Editorial: mathematical modeling of cardiovascular systems: from physiology to the clinic

Editorial on the Research Topic
Mathematical Modeling of Cardiovascular Systems: From Physiology to the Clinic

Biomedical research has enjoyed an eruption of reductionism where organs are reduced to their respective tissues, cells, and molecules. Although reductionism (including proteomic and genomic) is a powerful approach to provide information on the components of the system, integration of the components and their interactions are necessary to reveal organ phenotype and function ( Kassab et al., 2016 ; Heikhmakhtiar and Lim, 2018 ). Biological integration is inherently complex and requires mathematical “ gluing” of numerous components of the system with numerous variables. Physics-based mathematical modeling is well-tailored for the task where the number of variables is so great that it eludes intuition. The past several decades have enjoyed enormous advancements in mathematical modeling as computational capabilities have advanced substantially and experimental databases to inform and calibrate the computational models have become available to minimize the number of ad hoc assumptions. Applications of mathematical models have ranged from understanding basic biological systems to generate new hypotheses to patient-specific modeling and simulation studies that can aid medical diagnosis and help design optimal treatments ( Wenk et al., 2009 ). Modeling and simulation studies have an important role here because they provide knowledge and insights not available from existing data modalities and make evidence-generation more efficient by reducing the number of real patients needed in clinical studies. Furthermore, the utility of mathematical models remains largely removed from the clinic. The objective of this issue of “ Frontiers in Physiology” is to place the spotlight on computational integration with emphasis on validation and clinical translation of mathematical models.

The cardiovascular system which includes the heart and blood vessels has enjoyed prodigious efforts on both reduction and integration for understanding physiology and pathophysiology. For example, it is now clinical practice to compute pressure drop in the coronary arteries using computerized tomography (to determine fractional flow reserve; Tan et al., 2017 ), and non-invasively measure regional heart wall motion and strain in patients by using echocardiography or cardiac magnetic resonance imaging, in order to solve the inverse problem and quantify regional myocardial diastolic compliance and systolic contractility. On the latter, there still is no reliable force transducer or strain gauge for directly measuring forces or stresses in the intact heart wall without interfering with heart function. Clinicians, engineers, and physiologists have been interested in measurement and/or computation of heart wall stress since the early 1900s. Myocardial stress is an important boundary condition for coronary blood flow, and it is directly related to myocardial oxygen consumption ( Yin, 1981 ). The working hypothesis is that early and late cardiovascular disease treatment outcomes will be improved by using surgical and/or percutaneous techniques that reduce forces or stresses to homeostatic levels ( Grossman, 1980 ). Since regional heart wall stress cannot be measured reliably, mathematical modeling based on the conservation laws of continuum mechanics is needed. Because heart wall geometry (including its fibrous architecture) is fully 3D and the mechanical properties of beating myocardium are non-linear, there are no exact solutions of the governing differential equations of motion. Thus, numerical methods are required to find numerical approximations. The most versatile numerical method for simulation of cardiovascular system is the finite element (FE) method ( Guccione et al., 2010 ), which is widely used in the automotive and aerospace industries. Most FE models concerned with heart disease include only the left ventricular (LV) heart chamber because it is under the greatest stress (highest pressures) and thus, most prone to failure.

Currently, numerous computational models can simulate the behavior of the LV, a few models can simulate both ventricles, and even fewer can simulate all full four-chambered heart ( Baillargeon et al., 2014 ). The response of FE model is governed by realistic electrical, structural, and fluid flow physics. Significant recent improvements have been made in determination of diastolic ( Sack et al.) and systolic ( Sack et al., 2016 ) myocardial material properties in FE models of normal human LV ( Genet et al., 2014 ). Relatively few computational studies have investigated the effects of mechanical circulatory support devices on cardiac function using realistic geometries ( Lim et al., 2012 ; McCormick et al., 2013 ). A recent study of the effects of an LV assist device on acute left heart failure ( Sack et al., 2016 ) presented stress results in both ventricles.

Demands for medical device evidence development addressing patient safety, therapeutic efficacy, and cost-benefit determination are increasing. These demands are well-documented drivers of medical device clinical trial size, complexity, timelines, and costs. At the same time, and seemingly paradoxically, patient and provider demands for faster treatment access, and more innovative treatment options are increasing. Given the legitimacy of these stakeholders’ data demands, combined with the realities of resource limitations (budget and time), historic data acquisition techniques will soon no longer meet patient and provider needs. While numerous solutions have been hypothesized, the adoption of clinical trial solutions which make clinical trials more efficient and provide the data necessary to show safety and efficacy has been scarce. Modeling and simulation add value to the medical device innovation ecosystem by providing knowledge and insights not available from existing data modalities and making evidence generation more efficient by reducing the need for more burdensome types of data and information. We can more fully utilize modeling and simulation to minimize burden on patients, investigators, manufacturers and regulators, with mechanisms that support efficient review of novel products and approaches without compromising safety or effectiveness. “ In the future, computer-based modeling may change the way we think about device validation in other ways, allowing for much smaller clinical trials, …, in that some ‘ clinical’ information may be derived from simulation ( Faris and Shuren, 2017 ).” Stated another way, digital evidence (evidence from well-controlled modeling and simulation studies) could serve as “ clinical” information to support the evaluation of medical devices.

It is our hope that this edition of “ Frontiers in Physiology” titled “ Mathematical Modeling of Cardiovascular Systems: From Physiology to the Clinic” will serve to highlight powerful analytical tools that may advance to the clinic to improve safety and efficacy of existing treatments and advance novel diagnostics and therapeutics. In this spirit, we are delighted to include 23 articles contributed by 101 authors to this special issue. The topics range from models of cardiac and vascular system in health and disease from molecular to the organ system including simulations of cardiovascular devices. To promote these and other translational and transformational efforts, it is important to encourage further computational modeling that range from molecular, to cellular, to tissues and organs in partnership with experimental validations. Validation is key to this effort as it will boost confidence in the computational simulations that will translate from the academic laboratories to the clinic.

Author Contributions

All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

Baillargeon, B., Rebelo, N., Fox, D. D., Taylor, R. L., and Kuhl, E. (2014). The Living Heart Project: a robust and integrative simulator for human heart function. Eur. J. Mech. A Solids 48, 38–47. doi: 10. 1016/j. euromechsol. 2014. 04. 001

PubMed Abstract | CrossRef Full Text | Google Scholar

Faris, O., and Shuren, J. (2017). An FDA viewpoint on unique considerations for medical-device clinical trials. N. Engl. J. Med . 376, 1350–1357. doi: 10. 1056/NEJMra1512592

PubMed Abstract | CrossRef Full Text | Google Scholar

Genet, M., Lee, L. C., Nguyen, R., Haraldsson, H., Acevedo-Bolton, G., Zhang, Z., et al. (2014). Distribution of normal human left ventricular myofiber stress at end diastole and end systole: a target for in silico design of heart failure treatments. J. Appl. Physiol. (1985). 117, 142–152. doi: 10. 1152/japplphysiol. 00255. 2014

PubMed Abstract | CrossRef Full Text | Google Scholar

Grossman, W. (1980). Cardiac hypertrophy: useful adaptation or pathologic process? Am. J. Med . 69, 576–584. doi: 10. 1016/0002-9343(80)90471-4

CrossRef Full Text | Google Scholar

Guccione, J. M., Kassab, G. S., and Ratcliffe, M. B. (2010). Computational Cardiovascular Mechanics: Modeling and Applications in Heart Failure . New York, NY: Springer, XVI, 436.

Google Scholar

Heikhmakhtiar, A. K., and Lim, K. M. (2018). Computational prediction of the combined effect of CRT and LVAD on cardiac electromechanical delay in LBBB and RBBB. Comput. Math. Methods Med . 2018: 4253928. doi: 10. 1155/2018/4253928

PubMed Abstract | CrossRef Full Text | Google Scholar

Kassab, G. S., An, G., Sander, E. A., Miga, M., Guccione, J., Ji, S., et al. (2016). Augmenting surgery via multi-scale modeling and translational systems biology in the era of precision medicine: a multidisciplinary perspective. Ann. Biomed. Eng. 44, 2611–2625. doi: 10. 1007/s10439-016-1596-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Lim, K. M., Constantino, J., Gurev, V., Zhu, R., Shim, E. B., and Trayanova, N. A. (2012). Comparison of the effects of continuous and pulsatile left ventricular-assist devices on ventricular unloading using a cardiac electromechanics model. J. Physiol. Sci . 62, 11–19. doi: 10. 1007/s12576-011-0180-9

PubMed Abstract | CrossRef Full Text | Google Scholar

McCormick, M., Nordsletten, D. A., Kay, D., and Smith, N. P. (2013). Simulating left ventricular fluid-solid mechanics through the cardiac cycle under LVAD support. J. Comput. Phys . 244, 80–96. doi: 10. 1016/j. jcp. 2012. 08. 008

CrossRef Full Text | Google Scholar

Sack, K. L., Baillargeon, B., Acevedo-Bolton, G., Genet, M., Rebelo, N., Kuhl, E., et al. (2016). Partial LVAD restores ventricular outputs and normalizes LV but not RV stress distributions in the acutely failing heart in silico . Int. J. Artif. Organs . 39, 421–430. doi: 10. 5301/ijao. 5000520

CrossRef Full Text | Google Scholar

Tan, X. W., Zheng, Q., Shi, L., Gao, F., Allen, J. C. Jr., Coenen, A., et al. (2017). Combined diagnostic performance of coronary computed tomography angiography and computed tomography derived fractional flow reserve for the evaluation of myocardial ischemia: A meta-analysis. Int. J. Cardiol . 236, 100–106. doi: 10. 1016/j. ijcard. 2017. 02. 053

PubMed Abstract | CrossRef Full Text | Google Scholar

Wenk, J. F., Wall, S. T., Peterson, R. C., Helgerson, S. L., Sabbah, H. N., Burger, M., et al. (2009). A method for automatically optimizing medical devices for treating heart failure: designing polymeric injection patterns. J. Biomech. Eng . 131: 121011. doi: 10. 1115/1. 4000165

PubMed Abstract | CrossRef Full Text | Google Scholar

Yin, F. C. (1981). Ventricular wall stress. Circ. Res . 49, 829–842. doi: 10. 1161/01. RES. 49. 4. 829

PubMed Abstract | CrossRef Full Text | Google Scholar