An integrated electromechanical-growth heart model for simulating cardiac therapies

Simulated Effects of Acute Left Ventricular Myocardial Infarction on Mitral Regurgitation in an Ovine Model

Ischemic mitral regurgitation (IMR) occurs from incomplete coaptation of the mitral valve (MV) after myocardial infarction (MI), typically worsened by continued remodeling of the left ventricular (LV). The importance of LV remodeling is clear as IMR is induced by the post-MI dual mechanisms of mitral annular dilation and leaflet tethering from papillary muscle (PM) distension via the MV chordae tendineae (MVCT). However, the detailed etiology of IMR remains poorly understood, in large part due to the complex interactions of the MV and the post-MI LV remodeling processes. Given the patient-specific anatomical complexities of the IMR disease processes, simulation-based approaches represent an ideal approach to improve our understanding of this deadly disease. However, development of patient-specific models of left ventricle-mitral valve (LV-MV) interactions in IMR are complicated by the substantial variability and complexity of the MR etiology itself, making it difficult to extract underlying mechanisms from clinical data alone. To address these shortcomings, we developed a detailed ovine LV-MV finite element (FE) model based on extant comprehensive ovine experimental data. First, an extant ovine LV FE model (Sci. Rep. 2021 Jun 29;11(1):13466) was extended to incorporate the MV using a high fidelity ovine in vivo derived MV leaflet geometry. As it is not currently possible to image the MVCT in vivo, a functionally equivalent MVCT network was developed to create the final LV-MV model. Interestingly, in pilot studies, the MV leaflet strains did not agree well with known in vivo MV leaflet strain fields. We then incorporated previously reported MV leaflet prestrains (J. Biomech. Eng. 2023 Nov 1;145(11):111002) in the simulations. The resulting LV-MV model produced excellent agreement with the known in vivo ovine MV leaflet strains and deformed shapes in the normal state. We then simulated the effects of regional acute infarctions of varying sizes and anatomical locations by shutting down the local myocardial contractility. The remaining healthy (noninfarcted) myocardium mechanical behaviors were maintained, but allowed to adjust their active contractile patterns to maintain the prescribed pressure-volume loop behaviors in the acute post-MI state. For all cases studied, the LV-MV simulation demonstrated excellent agreement with known LV and MV in vivo strains and MV regurgitation orifice areas. Infarct location was shown to play a critical role in resultant MV leaflet strain fields. Specifically, extensional deformations of the posterior leaflets occurred in the posterobasal and laterobasal infarcts, while compressive deformations of the anterior leaflet were observed in the anterobasal infarct. Moreover, the simulated posterobasal infarct induced the largest MV regurgitation orifice area, consistent with experimental observations. The present study is the first detailed LV-MV simulation that reveals the important role of MV leaflet prestrain and functionally equivalent MVCT for accurate predictions of LV-MV interactions. Importantly, the current study further underscored simulation-based methods in understanding MV function as an integral part of the LV.

Multiscale Finite Element Modeling of Left Ventricular Growth in Simulations of Valve Disease

Multiscale models of the cardiovascular system are emerging as effective tools for investigating the mechanisms that drive ventricular growth and remodeling. These models can predict how molecular-level mechanisms impact organ-level structure and function and could provide new insights that help improve patient care. MyoFE is a multiscale computer framework that bridges molecular and organ-level mechanisms in a finite element model of the left ventricle that is coupled with the systemic circulation. In this study, we extend MyoFE to include a growth algorithm, based on volumetric growth theory, to simulate concentric growth (wall thickening/thinning) and eccentric growth (chamber dilation/constriction) in response to valvular diseases. Specifically in our model, concentric growth is controlled by time-averaged total stress along the fiber direction over a cardiac cycle while eccentric growth responds to time-averaged intracellular myofiber passive stress over a cardiac cycle. The new framework correctly predicted different forms of growth in response to two types of valvular diseases, namely aortic stenosis and mitral regurgitation. Furthermore, the model predicted that LV size and function are nearly restored (reversal of growth) when the disease-mimicking perturbation was removed in the simulations for each valvular disorder. In conclusion, the simulations suggest that time-averaged total stress along the fiber direction and time-averaged intracellular myofiber passive stress can be used to drive concentric and eccentric growth in simulations of valve disease.

Review of cardiac-coronary interaction and insights from mathematical modeling

Cardiac-coronary interaction is fundamental to the function of the heart. As one of the highest metabolic organs in the body, the cardiac oxygen demand is met by blood perfusion through the coronary vasculature. The coronary vasculature is largely embedded within the myocardial tissue which is continually contracting and hence squeezing the blood vessels. The myocardium-coronary vessel interaction is two-ways and complex. Here, we review the different types of cardiac-coronary interactions with a focus on insights gained from mathematical models. Specifically, we will consider the following: (1) myocardial-vessel mechanical interaction; (2) metabolic-flow interaction and regulation; (3) perfusion-contraction matching, and (4) chronic interactions between the myocardium and coronary vasculature. We also provide a discussion of the relevant experimental and clinical studies of different types of cardiac-coronary interactions. Finally, we highlight knowledge gaps, key challenges, and limitations of existing mathematical models along with future research directions to understand the unique myocardium-coronary coupling in the heart. This article is categorized under: Cardiovascular Diseases > Computational Models Cardiovascular Diseases > Biomedical Engineering Cardiovascular Diseases > Molecular and Cellular Physiology.

A multiscale finite element model of left ventricular mechanics incorporating baroreflex regulation

Cardiovascular function is regulated by a short-term hemodynamic baroreflex loop, which tries to maintain arterial pressure at a normal level. In this study, we present a new multiscale model of the cardiovascular system named MyoFE. This framework integrates a mechanistic model of contraction at the myosin level into a finite-element-based model of the left ventricle pumping blood through the systemic circulation. The model is coupled with a closed-loop feedback control of arterial pressure inspired by a baroreflex algorithm previously published by our team. The reflex loop mimics the afferent neuron pathway via a normalized signal derived from arterial pressure. The efferent pathway is represented by a kinetic model that simulates the net result of neural processing in the medulla and cell-level responses to autonomic drive. The baroreflex control algorithm modulates parameters such as heart rate and vascular tone of vessels in the lumped-parameter model of systemic circulation. In addition, it spatially modulates intracellular Ca2+ dynamics and molecular-level function of both the thick and the thin myofilaments in the left ventricle. Our study demonstrates that the baroreflex algorithm can maintain arterial pressure in the presence of perturbations such as acute cases of altered aortic resistance, mitral regurgitation, and myocardial infarction. The capabilities of this new multiscale model will be utilized in future research related to computational investigations of growth and remodeling.

An updated Lagrangian constrained mixture model of pathological cardiac growth and remodelling

Progressive left ventricular (LV) growth and remodelling (G&R) is often induced by volume and pressure overload, characterized by structural and functional adaptation through myocyte hypertrophy and extracellular matrix remodelling, which are dynamically regulated by biomechanical factors, inflammation, neurohormonal pathways, etc. When prolonged, it can eventually lead to irreversible heart failure. In this study, we have developed a new framework for modelling pathological cardiac G&R based on constrained mixture theory using an updated reference configuration, which is triggered by altered biomechanical factors to restore biomechanical homeostasis. Eccentric and concentric growth, and their combination have been explored in a patient-specific human LV model under volume and pressure overload. Eccentric growth is triggered by overstretching of myofibres due to volume overload, i.e. mitral regurgitation, whilst concentric growth is driven by excessive contractile stress due to pressure overload, i.e. aortic stenosis. Different biological constituent's adaptations under pathological conditions are integrated together, which are the ground matrix, myofibres and collagen network. We have shown that this constrained mixture-motivated G&R model can capture different phenotypes of maladaptive LV G&R, such as chamber dilation and wall thinning under volume overload, wall thickening under pressure overload, and more complex patterns under both pressure and volume overload. We have further demonstrated how collagen G&R would affect LV structural and functional adaption by providing mechanistic insight on anti-fibrotic interventions. This updated Lagrangian constrained mixture based myocardial G&R model has the potential to understand the turnover processes of myocytes and collagen due to altered local mechanical stimuli in heart diseases, and in providing mechanistic links between biomechanical factors and biological adaption at both the organ and cellular levels. Once calibrated with patient data, it can be used for assessing heart failure risk and designing optimal treatment therapies. STATEMENT OF SIGNIFICANCE: Computational modelling of cardiac G&R has shown high promise to provide insight into heart disease management when mechanistic understandings are quantified between biomechanical factors and underlying cellular adaptation processes. The kinematic growth theory has been dominantly used to phenomenologically describe the biological G&R process but neglecting underlying cellular mechanisms. We have developed a constrained mixture based G&R model with updated reference by taking into account different mechanobiological processes in the ground matrix, myocytes and collagen fibres. This G&R model can serve as a basis for developing more advanced myocardial G&R models further informed by patient data to assess heart failure risk, predict disease progression, select the optimal treatment by hypothesis testing, and eventually towards a truly precision cardiology using in-silico models.

Computational modelling of multi-temporal ventricular-vascular interactions during the progression of pulmonary arterial hypertension

A computational framework is developed to consider the concurrent growth and remodelling (G&R) processes occurring in the large pulmonary artery (PA) and right ventricle (RV), as well as ventricular-vascular interactions during the progression of pulmonary arterial hypertension (PAH). This computational framework couples the RV and the proximal PA in a closed-loop circulatory system that operates in a short timescale of a cardiac cycle, and evolves over a long timescale due to G&R processes in the PA and RV. The framework predicts changes in haemodynamics (e.g. 68.2% increase in mean PA pressure), RV geometry (e.g. 38% increase in RV end-diastolic volume) and PA tissue microstructure (e.g. 90% increase in collagen mass) that are consistent with clinical and experimental measurements of PAH. The framework also predicts that a reduction in RV contractility is associated with long-term RV chamber dilation, a common biomarker observed in the late-stage PAH. Sensitivity analyses on the G&R rate constants show that large PA stiffening (both short and long term) is affected by RV remodelling more than the reverse. This framework can serve as a foundation for the future development of a more predictive and comprehensive cardiovascular G&R model with realistic heart and vascular geometries.

Computational cardiac physiology for new modelers: Origins, foundations, and future

Mathematical models of the cardiovascular system have come a long way since they were first introduced in the early 19th century. Driven by a rapid development of experimental techniques, numerical methods, and computer hardware, detailed models that describe physical scales from the molecular level up to organs and organ systems have been derived and used for physiological research. Mathematical and computational models can be seen as condensed and quantitative formulations of extensive physiological knowledge and are used for formulating and testing hypotheses, interpreting and directing experimental research, and have contributed substantially to our understanding of cardiovascular physiology. However, in spite of the strengths of mathematics to precisely describe complex relationships and the obvious need for the mathematical and computational models to be informed by experimental data, there still exist considerable barriers between experimental and computational physiological research. In this review, we present a historical overview of the development of mathematical and computational models in cardiovascular physiology, including the current state of the art. We further argue why a tighter integration is needed between experimental and computational scientists in physiology, and point out important obstacles and challenges that must be overcome in order to fully realize the synergy of experimental and computational physiological research.

Computational models of ventricular mechanics and adaptation in response to right-ventricular pressure overload

Pulmonary arterial hypertension (PAH) is associated with substantial remodeling of the right ventricle (RV), which may at first be compensatory but at a later stage becomes detrimental to RV function and patient survival. Unlike the left ventricle (LV), the RV remains understudied, and with its thin-walled crescent shape, it is often modeled simply as an appendage of the LV. Furthermore, PAH diagnosis is challenging because it often leaves the LV and systemic circulation largely unaffected. Several treatment strategies such as atrial septostomy, right ventricular assist devices (RVADs) or RV resynchronization therapy have been shown to improve RV function and the quality of life in patients with PAH. However, evidence of their long-term efficacy is limited and lung transplantation is still the most effective and curative treatment option. As such, the clinical need for improved diagnosis and treatment of PAH drives a strong need for increased understanding of drivers and mechanisms of RV growth and remodeling (G&R), and more generally for targeted research into RV mechanics pathology. Computational models stand out as a valuable supplement to experimental research, offering detailed analysis of the drivers and consequences of G&R, as well as a virtual test bench for exploring and refining hypotheses of growth mechanisms. In this review we summarize the current efforts towards understanding RV G&R processes using computational approaches such as reduced-order models, three dimensional (3D) finite element (FE) models, and G&R models. In addition to an overview of the relevant literature of RV computational models, we discuss how the models have contributed to increased scientific understanding and to potential clinical treatment of PAH patients.

Numerical evaluation of cardiac mechanical markers as estimators of the electrical activation time

Recent advances in the development of noninvasive cardiac imaging technologies have made it possible to measure longitudinal and circumferential strains at a high spatial resolution also at intramural level. Local mechanical activation times derived from these strains can be used as noninvasive estimates of electrical activation, in order to determine, eg, the origin of premature ectopic beats during focal arrhythmias or the pathway of reentrant circuits. The aim of this work is to assess the reliability of mechanical activation time markers derived from longitudinal and circumferential strains, denoted by AT_ell and AT_ecc , respectively, by means of three-dimensional cardiac electromechanical simulations. These markers are compared against the electrical activation time (AT_v ), computed from the action potential waveform, and the reference mechanical activation markers derived from the active tension and fiber strain waveforms, denoted by AT_ta and AT_eff , respectively. Our numerical simulations are based on a strongly coupled electromechanical model, including bidomain representation of the cardiac tissue, mechanoelectric (ie, stretch-activated channels) and geometric feedbacks, transversely isotropic strain energy function for the description of passive mechanics and detailed membrane and excitation-contraction coupling models. The results have shown that, during endocardial and epicardial ectopic stimulations, all the mechanical markers considered are highly correlated with AT_v , exhibiting correlation coefficients larger than 0.8. However, during multiple endocardial stimulations, mimicking the ventricular sinus rhythm, the mechanical markers are less correlated with the electrical activation time, because of the more complex resulting excitation sequence. Moreover, the inspection of the endocardial and epicardial isochrones has shown that the AT_ell and AT_ecc mechanical activation sequences reproduce only some qualitative features of the electrical activation sequence, such as the areas of early and late activation, but in some cases, they might yield wrong excitation sources and significantly different isochrones patterns.

Multiscale simulations of left ventricular growth and remodeling

Cardiomyocytes can adapt their size, shape, and orientation in response to altered biomechanical or biochemical stimuli. The process by which the heart undergoes structural changes-affecting both geometry and material properties-in response to altered ventricular loading, altered hormonal levels, or mutant sarcomeric proteins is broadly known as cardiac growth and remodeling (G&R). Although it is likely that cardiac G&R initially occurs as an adaptive response of the heart to the underlying stimuli, prolonged pathological changes can lead to increased risk of atrial fibrillation, heart failure, and sudden death. During the past few decades, computational models have been extensively used to investigate the mechanisms of cardiac G&R, as a complement to experimental measurements. These models have provided an opportunity to quantitatively study the relationships between the underlying stimuli (primarily mechanical) and the adverse outcomes of cardiac G&R, i.e., alterations in ventricular size and function. State-of-the-art computational models have shown promise in predicting the progression of cardiac G&R. However, there are still limitations that need to be addressed in future works to advance the field. In this review, we first outline the current state of computational models of cardiac growth and myofiber remodeling. Then, we discuss the potential limitations of current models of cardiac G&R that need to be addressed before they can be utilized in clinical care. Finally, we briefly discuss the next feasible steps and future directions that could advance the field of cardiac G&R.

Sex Differences in Drug-Induced Arrhythmogenesis

The electrical activity in the heart varies significantly between men and women and results in a sex-specific response to drugs. Recent evidence suggests that women are more than twice as likely as men to develop drug-induced arrhythmia with potentially fatal consequences. Yet, the sex-specific differences in drug-induced arrhythmogenesis remain poorly understood. Here we integrate multiscale modeling and machine learning to gain mechanistic insight into the sex-specific origin of drug-induced cardiac arrhythmia at differing drug concentrations. To quantify critical drug concentrations in male and female hearts, we identify the most important ion channels that trigger male and female arrhythmogenesis, and create and train a sex-specific multi-fidelity arrhythmogenic risk classifier. Our study reveals that sex differences in ion channel activity, tissue conductivity, and heart dimensions trigger longer QT-intervals in women than in men. We quantify the critical drug concentration for dofetilide, a high risk drug, to be seven times lower for women than for men. Our results emphasize the importance of including sex as an independent biological variable in risk assessment during drug development. Acknowledging and understanding sex differences in drug safety evaluation is critical when developing novel therapeutic treatments on a personalized basis. The general trends of this study have significant implications on the development of safe and efficacious new drugs and the prescription of existing drugs in combination with other drugs.

Computational models of cardiac hypertrophy

Cardiac hypertrophy, defined as an increase in mass of the heart, is a complex process driven by simultaneous changes in hemodynamics, mechanical stimuli, and hormonal inputs. It occurs not only during pre- and post-natal development but also in adults in response to exercise, pregnancy, and a range of cardiovascular diseases. One of the most exciting recent developments in the field of cardiac biomechanics is the advent of computational models that are able to accurately predict patterns of heart growth in many of these settings, particularly in cases where changes in mechanical loading of the heart play an import role. These emerging models may soon be capable of making patient-specific growth predictions that can be used to guide clinical interventions. Here, we review the history and current state of cardiac growth models and highlight three main limitations of current approaches with regard to future clinical application: their inability to predict the regression of heart growth after removal of a mechanical overload, inability to account for evolving hemodynamics, and inability to incorporate known growth effects of drugs and hormones on heart growth. Next, we outline growth mechanics approaches used in other fields of biomechanics and highlight some potential lessons for cardiac growth modeling. Finally, we propose a multiscale modeling approach for future studies that blends tissue-level growth models with cell-level signaling models to incorporate the effects of hormones in the context of pregnancy-induced heart growth.

Role of coronary flow regulation and cardiac-coronary coupling in mechanical dyssynchrony associated with right ventricular pacing

Mechanical dyssynchrony (MD) affects left ventricular (LV) mechanics and coronary perfusion. To understand the multifactorial effects of MD, we developed a computational model that bidirectionally couples the systemic circulation with the LV and coronary perfusion with flow regulation. In the model, coronary flow in the left anterior descending (LAD) and left circumflex (LCX) arteries affects the corresponding regional contractility based on a prescribed linear LV contractility-coronary flow relationship. The model is calibrated with experimental measurements of LV pressure and volume, as well as LAD and LCX flow rate waveforms acquired under regulated and fully dilated conditions from a swine under right atrial (RA) pacing. The calibrated model is applied to simulate MD. The model can simultaneously reproduce the reduction in mean LV pressure (39.3%), regulated flow (LAD: 7.9%; LCX: 1.9%), LAD passive flow (21.6%), and increase in LCX passive flow (15.9%). These changes are associated with right ventricular pacing compared with RA pacing measured in the same swine only when LV contractility is affected by flow alterations with a slope of 1.4 mmHg/mL² in a contractility-flow relationship. In sensitivity analyses, the model predicts that coronary flow reserve (CFR) decreases and increases in the LAD and LCX with increasing delay in LV free wall contraction. These findings suggest that asynchronous activation associated with MD impacts 1) the loading conditions that further affect the coronary flow, which may explain some of the changes in CFR, and 2) the coronary flow that reduces global contractility, which contributes to the reduction in LV pressure.NEW & NOTEWORTHY A computational model that couples the systemic circulation of the left ventricular (LV) and coronary perfusion with flow regulation is developed to study the effects of mechanical dyssynchrony. The delayed contraction in the LV free wall with respect to the septum has a significant effect on LV function and coronary flow reserve.

Multiscale Modeling of Cardiovascular Function Predicts That the End-Systolic Pressure Volume Relationship Can Be Targeted via Multiple Therapeutic Strategies

Most patients who develop heart failure are unable to elevate their cardiac output on demand due to impaired contractility and/or reduced ventricular filling. Despite decades of research, few effective therapies for heart failure have been developed. In part, this may reflect the difficulty of predicting how perturbations to molecular-level mechanisms that are induced by drugs will scale up to modulate system-level properties such as blood pressure. Computer modeling might help with this process and thereby accelerate the development of better therapies for heart failure. This manuscript presents a new multiscale model that uses a single contractile element to drive an idealized ventricle that pumps blood around a closed circulation. The contractile element was formed by linking an existing model of dynamically coupled myofilaments with a well-established model of myocyte electrophysiology. The resulting framework spans from molecular-level events (including opening of ion channels and transitions between different myosin states) to properties such as ejection fraction that can be measured in patients. Initial calculations showed that the model reproduces many aspects of normal cardiovascular physiology including, for example, pressure-volume loops. Subsequent sensitivity tests then quantified how each model parameter influenced a range of system level properties. The first key finding was that the End Systolic Pressure Volume Relationship, a classic index of cardiac contractility, was ∼50% more sensitive to parameter changes than any other system-level property. The second important result was that parameters that primarily affect ventricular filling, such as passive stiffness and Ca²⁺ reuptake via sarco/endoplasmic reticulum Ca²⁺-ATPase (SERCA), also have a major impact on systolic properties including stroke work, myosin ATPase, and maximum ventricular pressure. These results reinforce the impact of diastolic function on ventricular performance and identify the End Systolic Pressure Volume Relationship as a particularly sensitive system-level property that can be targeted using multiple therapeutic strategies.

Biomechanics of infarcted left ventricle: a review of modelling

Mathematical modelling in biomechanics of infarcted left ventricle (LV) serves as an indispensable tool for remodelling mechanism exploration, LV biomechanical property estimation and therapy assessment after myocardial infarction (MI). However, a review of mathematical modelling after MI has not been seen in the literature so far. In the paper, a systematic review of mathematical models in biomechanics of infarcted LV was established. The models include comprehensive cardiovascular system model, essential LV pressure-volume and stress-stretch models, constitutive laws for passive myocardium and scars, tension models for active myocardium, collagen fibre orientation optimization models, fibroblast and collagen fibre growth/degradation models and integrated growth-electro-mechanical model after MI. The primary idea, unique characteristics and key equations of each model were identified and extracted. Discussions on the models were provided and followed research issues on them were addressed. Considerable improvements in the cardiovascular system model, LV aneurysm model, coupled agent-based models and integrated electro-mechanical-growth LV model are encouraged. Substantial attention should be paid to new constitutive laws with respect to stress-stretch curve and strain energy function for infarcted passive myocardium, collagen fibre orientation optimization in scar, cardiac rupture and tissue damage and viscoelastic effect post-MI in the future.

Computational modeling of cardiac growth and remodeling in pressure overloaded hearts-Linking microstructure to organ phenotype

Cardiac growth and remodeling (G&R) refers to structural changes in myocardial tissue in response to chronic alterations in loading conditions. One such condition is pressure overload where elevated wall stresses stimulate the growth in cardiomyocyte thickness, associated with a phenotype of concentric hypertrophy at the organ scale, and promote fibrosis. The initial hypertrophic response can be considered adaptive and beneficial by favoring myocyte survival, but over time if pressure overload conditions persist, maladaptive mechanisms favoring cell death and fibrosis start to dominate, ultimately mediating the transition towards an overt heart failure phenotype. The underlying mechanisms linking biological factors at the myocyte level to biomechanical factors at the systemic and organ level remain poorly understood. Computational models of G&R show high promise as a unique framework for providing a quantitative link between myocardial stresses and strains at the organ scale to biological regulatory processes at the cellular level which govern the hypertrophic response. However, microstructurally motivated, rigorously validated computational models of G&R are still in their infancy. This article provides an overview of the current state-of-the-art of computational models to study cardiac G&R. The microstructure and mechanosensing/mechanotransduction within cells of the myocardium is discussed and quantitative data from previous experimental and clinical studies is summarized. We conclude with a discussion of major challenges and possible directions of future research that can advance the current state of cardiac G&R computational modeling. STATEMENT OF SIGNIFICANCE: The mechanistic links between organ-scale biomechanics and biological factors at the cellular size scale remain poorly understood as these are largely elusive to investigations using experimental methodology alone. Computational G&R models show high promise to establish quantitative links which allow more mechanistic insight into adaptation mechanisms and may be used as a tool for stratifying the state and predict the progression of disease in the clinic. This review provides a comprehensive overview of research in this domain including a summary of experimental data. Thus, this study may serve as a basis for the further development of more advanced G&R models which are suitable for making clinical predictions on disease progression or for testing hypotheses on pathogenic mechanisms using in-silico models.

Multisite pacing and myocardial scars: a computational study

Cardiac resynchronization therapy (CRT) is a frequently effective treatment modality for dyssynchronous heart failure, however, 30% of patients do not respond, usually due to suboptimal activation of the left ventricle (LV). Multisite pacing (MSP) may increase the response rate, but its effect in the presence of myocardial scars is not fully understood. We use a computational model to study the outcome of MSP in an LV with scars in two different locations and of two different sizes. The LV was stimulated from anterior, posterior and lateral locations individually and in pairs, while a septal stimulation site represented right ventricular (RV) pacing. Intraventricular pressures were measured, and outcomes evaluated in terms of maximum LV pressure gradient (dP/dt_max)- change compared to isolated RV pacing. The best result obtained using various LV pacing locations included a combination of sites remote from scars and the septum. The highest dP/dt_max increase was achieved, regardless of scar size, using MSP with one pacing site located on the LV free wall opposite to the scar and one site opposite to the septum. These in silico modelling results suggest that making placement of pacing electrodes dependent on location of scarring, may alter acute haemodynamics and that such modelling may contribute to future CRT optimization.

Multiscale Models of Cardiac Muscle Biophysics and Tissue Remodeling in Hypertrophic Cardiomyopathies

Myocardial hypertrophy is the result of sustained perturbations to the mechanical and/or neurohormonal homeostasis of cardiac cells and is driven by integrated, multiscale biophysical and biochemical processes that are currently not well defined. In this brief review, we highlight recent computational and experimental models of cardiac hypertrophy that span mechanisms from the molecular level to the tissue level. Specifically, we focus on: (i) molecular-level models of the structural dynamics of sarcomere proteins in hypertrophic hearts, (ii) cellular-level models of excitation-contraction coupling and mechanosensitive signaling in disease-state myocytes, and (iii) organ-level models of myocardial growth kinematics and predictors thereof. Finally, we discuss how spanning these scales and combining multiple experimental/computational models will provide new information about the processes governing hypertrophy and potential methods to prevent or reverse them.

Model of Anisotropic Reverse Cardiac Growth in Mechanical Dyssynchrony

Based on recent single-cell experiments showing that longitudinal myocyte stretch produces both parallel and serial addition of sarcomeres, we developed an anisotropic growth constitutive model with elastic myofiber stretch as the growth stimuli to simulate long-term changes in biventricular geometry associated with alterations in cardiac electromechanics. The constitutive model is developed based on the volumetric growth framework. In the model, local growth evolutions of the myocyte’s longitudinal and transverse directions are driven by the deviations of maximum elastic myofiber stretch over a cardiac cycle from its corresponding local homeostatic set point, but with different sensitivities. Local homeostatic set point is determined from a simulation with normal activation pattern. The growth constitutive model is coupled to an electromechanics model and calibrated based on both global and local ventricular geometrical changes associated with chronic left ventricular free wall pacing found in previous animal experiments. We show that the coupled electromechanics-growth model can quantitatively reproduce the following: (1) Thinning and thickening of the ventricular wall respectively at early and late activated regions and (2) Global left ventricular dilation as measured in experiments. These findings reinforce the role of elastic myofiber stretch as a growth stimulant at both cellular level and tissue-level.

Interactions between structural remodeling and volumetric growth in right ventricle in response to pulmonary arterial hypertension

Pulmonary arterial hypertension (PAH) exerts substantial pressure overload on the right ventricle (RV). The associated RV free wall (RVFW) adaptation could consist of myocardial hypertrophy, augmented intrinsic contractility, collagen fibrosis, and structural remodeling in an attempt to cope with pressure overload. If RVFW adaptation cannot maintain the RV stroke volume, RV dilation will prevail as an exit mechanism which usually decompensates the RV function leading to RV failure. Our knowledge of the factors determining the transition from the upper limit of RVFW adaptation to RV decompensation and the role of fiber remodeling events in this transition remains very limited. Computational heart models that connect the growth and remodeling (G\&R) events at the fiber and tissue levels with alterations in the organ-level function are essential to predict the temporal order and the compensatory level of the underlying mechanisms. In this work, building upon our recent rodent heart models (RHM) of PAH, we integrated mathematical models that describe time-evolution volumetric growth of the RV and structural remodeling of the RVFW. Results suggest that augmentation of the intrinsic contractility of myofibers accompanied by an increase in passive stiffness of RVFW is among the first remodeling events through which the RV strives to maintain the cardiac output. Interestingly, we found that the observed reorientation of the myofibers towards the longitudinal (apex-to-base) direction was a maladaptive mechanism that impaired the contractile pattern of RVFW and advanced along with RV dilation at later stages of PAH development.

Coupled agent-based and hyperelastic modelling of the left ventricle post-myocardial infarction

Understanding the healing and remodelling processes induced by myocardial infarction (MI) of the heart is important, and the mechanical properties of the myocardium post-MI can be indicative for effective treatments aimed at avoiding eventual heart failure. MI remodelling is a multiscale feedback process between the mechanical loading and cellular adaptation. In this paper, we use an agent-based model to describe collagen remodelling by fibroblasts regulated by chemical and mechanical cues after acute MI, and upscale into a finite element 3D left ventricular model. We model the dispersed collagen fibre structure using the angular integration method and have incorporated a collagen fibre tension-compression switch in the left ventricle (LV) model. This enables us to study the scar healing (collagen deposition, degradation, and reorientation) of a rat heart post-MI. Our results, in terms of collagen accumulation and alignment, compare well with published experimental data. In addition, we show that different shapes of the MI region can affect the collagen remodelling, and in particular, the mechanical cue plays an important role in the healing process.

A fiduciary marker-based framework to assess heterogeneity and anisotropy of right ventricular epicardial strains in the beating ovine heart

Quantifying ventricular deformation in health and disease is critical to our understanding of normal heart function, heart disease mechanisms, and the effect of medical treatments. Imaging modalities have been developed that can measure ventricular deformation non-invasively. However, because of the small thickness, complex shape, and anatomic position of the right ventricle, using these technologies to determine its deformation remains challenging. Here we develop a first fiduciary marker-based method to assess heterogeneity and anisotropy of right ventricular epicardial strain across the entire free wall. To this end, we combine a high-density array of sonomicrometry crystals implanted across the entire right ventricular epicardial surface with a subdivision surface algorithm and a large deformation kinematics framework. We demonstrate our approach on four beating ovine hearts and present a preliminary regional analysis of circumferential, longitudinal, and areal strain. Moreover, we illustrate maps of the same strains across the entire right ventricular epicardial surface to highlight their spatial heterogeneity and anisotropy. We observe in these animals that RV epicardial strains vary throughout the cardiac cycle, are heterogeneous across the RV free wall, and are anisotropic with larger compressive strains, i.e., contraction, in the longitudinal direction than in the circumferential direction. Average peak compressive strains vary by region between -3.34% and -8.29% in circumferential direction, and -4.02% and -10.57% in longitudinal direction. In summary, we introduce an experimental framework that will allow us to study disease- and device-induced deformations, and long-term consequences of these deformations, including heterogeneous and anisotropic effects.

Organ-level validation of a cross-bridge cycling descriptor in a left ventricular finite element model: effects of ventricular loading on myocardial strains

Although detailed cell‐based descriptors of cross‐bridge cycling have been applied in finite element (FE) heart models to describe ventricular mechanics, these multiscale models have never been tested rigorously to determine if these descriptors, when scaled up to the organ‐level, are able to reproduce well‐established organ‐level physiological behaviors. To address this void, we here validate a left ventricular (LV) FE model that is driven by a cell‐based cross‐bridge cycling descriptor against key organ‐level heart physiology. The LV FE model was coupled to a closed‐loop lumped parameter circulatory model to simulate different ventricular loading conditions (preload and afterload) and contractilities. We show that our model is able to reproduce a linear end‐systolic pressure volume relationship, a curvilinear end‐diastolic pressure volume relationship and a linear relationship between myocardial oxygen consumption and pressure–volume area. We also show that the validated model can predict realistic LV strain‐time profiles in the longitudinal, circumferential, and radial directions. The predicted strain‐time profiles display key features that are consistent with those measured in humans, such as having similar peak strains, time‐to‐peak‐strain, and a rapid change in strain during atrial contraction at late‐diastole. Our model shows that the myocardial strains are sensitive to not only LV contractility, but also to the LV loading conditions, especially to a change in afterload. This result suggests that caution must be exercised when associating changes in myocardial strain with changes in LV contractility. The methodically validated multiscale model will be used in future studies to understand human heart diseases.

Predicting drug-induced arrhythmias by multiscale modeling

Drugs often have undesired side effects. In the heart, they can induce lethal arrhythmias such as torsades de pointes. The risk evaluation of a new compound is costly and can take a long time, which often hinders the development of new drugs. Here, we establish a high-resolution, multiscale computational model to quickly assess the cardiac toxicity of new and existing drugs. The input of the model is the drug-specific current block from single cell electrophysiology; the output is the spatio-temporal activation profile and the associated electrocardiogram. We demonstrate the potential of our model for a low-risk drug, ranolazine, and a high-risk drug, quinidine: For ranolazine, our model predicts a prolonged QT interval of 19.4% compared with baseline and a regular sinus rhythm at 60.15 beats per minute. For quinidine, our model predicts a prolonged QT interval of 78.4% and a spontaneous development of torsades de pointes both in the activation profile and in the electrocardiogram. Our model reveals the mechanisms by which electrophysiological abnormalities propagate across the spatio-temporal scales, from specific channel blockage, via altered single cell action potentials and prolonged QT intervals, to the spontaneous emergence of ventricular tachycardia in the form of torsades de pointes. Our model could have important implications for researchers, regulatory agencies, and pharmaceutical companies on rationalizing safe drug development and reducing the time-to-market of new drugs.

Computational Modeling for Cardiac Resynchronization Therapy

Cardiac resynchronization therapy (CRT) is an effective treatment for heart failure (HF) patients with an electrical substrate pathology causing ventricular dyssynchrony. However 40-50% of patients do not respond to treatment. Cardiac modeling of the electrophysiology, electromechanics, and hemodynamics of the heart has been used to study mechanisms behind HF pathology and CRT response. Recently, multi-scale dyssynchronous HF models have been used to study optimal device settings and optimal lead locations, investigate the underlying cardiac pathophysiology, as well as investigate emerging technologies proposed to treat cardiac dyssynchrony. However the breadth of patient and experimental data required to create and parameterize these models and the computational resources required currently limits the use of these models to small patient numbers. In the future, once these technical challenges are overcome, biophysically based models of the heart have the potential to become a clinical tool to aid in the diagnosis and treatment of HF.

Preconditioned augmented Lagrangian formulation for nearly incompressible cardiac mechanics

Computational modeling of the heart is a subject of substantial medical and scientific interest, which may contribute to increase the understanding of several phenomena associated with cardiac physiological and pathological states. Modeling the mechanics of the heart have led to considerable insights, but it still represents a complex and a demanding computational problem, especially in a strongly coupled electromechanical setting. Passive cardiac tissue is commonly modeled as hyperelastic and is characterized by quasi-incompressible, orthotropic, and nonlinear material behavior. These factors are known to be very challenging for the numerical solution of the model. The near-incompressibility is known to cause numerical issues such as the well-known locking phenomenon and ill-conditioning of the stiffness matrix. In this work, the augmented Lagrangian method is used to handle the nearly incompressible condition. This approach can potentially improve computational performance by reducing the condition number of the stiffness matrix and thereby improving the convergence of iterative solvers. We also improve the performance of iterative solvers by the use of an algebraic multigrid preconditioner. Numerical results of the augmented Lagrangian method combined with a preconditioned iterative solver for a cardiac mechanics benchmark suite are presented to show its improved performance.

High Spatial Resolution Multi-Organ Finite Element Modeling of Ventricular-Arterial Coupling

While it has long been recognized that bi-directional interaction between the heart and the vasculature plays a critical role in the proper functioning of the cardiovascular system, a comprehensive study of this interaction has largely been hampered by a lack of modeling framework capable of simultaneously accommodating high-resolution models of the heart and vasculature. Here, we address this issue and present a computational modeling framework that couples finite element (FE) models of the left ventricle (LV) and aorta to elucidate ventricular-arterial coupling in the systemic circulation. We show in a baseline simulation that the framework predictions of (1) LV pressure-volume loop, (2) aorta pressure-diameter relationship, (3) pressure-waveforms of the aorta, LV, and left atrium (LA) over the cardiac cycle are consistent with the physiological measurements found in healthy human. To develop insights of ventricular-arterial interactions, the framework was then used to simulate how alterations in the geometrical or, material parameter(s) of the aorta affect the LV and vice versa. We show that changing the geometry and microstructure of the aorta model in the framework led to changes in the functional behaviors of both LV and aorta that are consistent with experimental observations. On the other hand, changing contractility and passive stiffness of the LV model in the framework also produced changes in both the LV and aorta functional behaviors that are consistent with physiology principles.

Interpreting Activation Mapping of Atrial Fibrillation: A Hybrid Computational/Physiological Study

Atrial fibrillation is the most common rhythm disorder of the heart associated with a rapid and irregular beating of the upper chambers. Activation mapping remains the gold standard to diagnose and interpret atrial fibrillation. However, fibrillatory activation maps are highly sensitive to far-field effects, and often disagree with other optical mapping modalities. Here we show that computational modeling can identify spurious non-local components of atrial fibrillation electrograms and improve activation mapping. We motivate our approach with a cohort of patients with potential drivers of persistent atrial fibrillation. In a computational study using a monodomain Maleckar model, we demonstrate that in organized rhythms, electrograms successfully track local activation, whereas in atrial fibrillation, electrograms are sensitive to spiral wave distance and number, spiral tip trajectories, and effects of fibrosis. In a clinical study, we analyzed n = 15 patients with persistent atrial fibrillation that was terminated by limited ablation. In five cases, traditional activation maps revealed a spiral wave at sites of termination; in ten cases, electrogram timings were ambiguous and activation maps showed incomplete reentry. By adjusting electrogram timing through computational modeling, we found rotational activation, which was undetectable with conventional methods. Our results demonstrate that computational modeling can identify non-local deflections to improve activation mapping and explain how and where ablation can terminate persistent atrial fibrillation. Our hybrid computational/physiological approach has the potential to optimize map-guided ablation and improve ablation therapy in atrial fibrillation.

Predicting the cardiac toxicity of drugs using a novel multiscale exposure-response simulator

A common but serious side effect of many drugs is torsades de pointes, a rhythm disorder that can have fatal consequences. Torsadogenic risk has traditionally been associated with blockage of a specific potassium channel and an increased recovery period in the electrocardiogram. However, the mechanisms that trigger torsades de pointes remain incompletely understood. Here we establish a computational model to explore how drug-induced effects propagate from the single channel, via the single cell, to the whole heart level. Our mechanistic exposure-response simulator translates block-concentration characteristics for arbitrary drugs into three-dimensional excitation profiles and electrocardiogram recordings to rapidly assess torsadogenic risk. For the drug of dofetilide, we show that this risk is highly dose-dependent: at a concentration of 1x, QT prolongation is 55% but the heart maintains its regular sinus rhythm; at 5.7x, QT prolongation is 102% and the heart spontaneously transitions into torsades de points; at 30x, QT prolongation is 132% and the heart adapts a quasi-depolarized state with numerous rapidly flickering local excitations. Our simulations suggest that neither potassium channel blockage nor QT interval prolongation alone trigger torsades de pointes. The underlying mechanism predicted by our model is early afterdepolarization, which translates into pronounced U waves in the electrocardiogram, a signature that is correctly predicted by our model. Beyond the risk assessment of existing drugs, our exposure-response simulator can become a powerful tool to optimize the co-administration of drugs and, ultimately, guide the design of new drugs toward reducing life threatening drug-induced rhythm disorders in the heart.

Structural and Functional Differences Between Porcine Aorta and Vena Cava

Elastin and collagen fibers are the major load-bearing extracellular matrix (ECM) constituents of the vascular wall. Arteries function differently than veins in the circulatory system; however as a result from several treatment options, veins are subjected to sudden elevated arterial pressure. It is thus important to recognize the fundamental structure and function differences between a vein and an artery. Our research compared the relationship between biaxial mechanical function and ECM structure of porcine thoracic aorta and inferior vena cava. Our study suggests that aorta contains slightly more elastin than collagen due to the cyclical extensibility, but vena cava contains almost four times more collagen than elastin to maintain integrity. Furthermore, multiphoton imaging of vena cava showed longitudinally oriented elastin and circumferentially oriented collagen that is recruited at supraphysiologic stress, but low levels of strain. However in aorta, elastin is distributed uniformly, and the primarily circumferentially oriented collagen is recruited at higher levels of strain than vena cava. These structural observations support the functional finding that vena cava is highly anisotropic with the longitude being more compliant and the circumference stiffening substantially at low levels of strain. Overall, our research demonstrates that fiber distributions and recruitment should be considered in addition to relative collagen and elastin contents. Also, the importance of accounting for the structural and functional differences between arteries and veins should be taken into account when considering disease treatment options.

Patient-Specific Computational Analysis of Ventricular Mechanics in Pulmonary Arterial Hypertension

Patient-specific biventricular computational models associated with a normal subject and a pulmonary arterial hypertension (PAH) patient were developed to investigate the disease effects on ventricular mechanics. These models were developed using geometry reconstructed from magnetic resonance (MR) images, and constitutive descriptors of passive and active mechanics in cardiac tissues. Model parameter values associated with ventricular mechanical properties and myofiber architecture were obtained by fitting the models with measured pressure-volume loops and circumferential strain calculated from MR images using a hyperelastic warping method. Results show that the peak right ventricle (RV) pressure was substantially higher in the PAH patient (65 mmHg versus 20 mmHg), who also has a significantly reduced ejection fraction (EF) in both ventricles (left ventricle (LV): 39% versus 66% and RV: 18% versus 64%). Peak systolic circumferential strain was comparatively lower in both the left ventricle (LV) and RV free wall (RVFW) of the PAH patient (LV: -6.8% versus -13.2% and RVFW: -2.1% versus -9.4%). Passive stiffness, contractility, and myofiber stress in the PAH patient were all found to be substantially increased in both ventricles, whereas septum wall in the PAH patient possessed a smaller curvature than that in the LV free wall. Simulations using the PAH model revealed an approximately linear relationship between the septum curvature and the transseptal pressure gradient at both early-diastole and end-systole. These findings suggest that PAH can induce LV remodeling, and septum curvature measurements may be useful in quantifying transseptal pressure gradient in PAH patients.

Improved identifiability of myocardial material parameters by an energy-based cost function

Myocardial stiffness is a valuable clinical biomarker for the monitoring and stratification of heart failure (HF). Cardiac finite element models provide a biomechanical framework for the assessment of stiffness through the determination of the myocardial constitutive model parameters. The reported parameter intercorrelations in popular constitutive relations, however, obstruct the unique estimation of material parameters and limit the reliable translation of this stiffness metric to clinical practice. Focusing on the role of the cost function (CF) in parameter identifiability, we investigate the performance of a set of geometric indices (based on displacements, strains, cavity volume, wall thickness and apicobasal dimension of the ventricle) and a novel CF derived from energy conservation. Our results, with a commonly used transversely isotropic material model (proposed by Guccione et al.), demonstrate that a single geometry-based CF is unable to uniquely constrain the parameter space. The energy-based CF, conversely, isolates one of the parameters and in conjunction with one of the geometric metrics provides a unique estimation of the parameter set. This gives rise to a new methodology for estimating myocardial material parameters based on the combination of deformation and energetics analysis. The accuracy of the pipeline is demonstrated in silico, and its robustness in vivo, in a total of 8 clinical data sets (7 HF and one control). The mean identified parameters of the Guccione material law were [Formula: see text] and [Formula: see text] ([Formula: see text], [Formula: see text], [Formula: see text]) for the HF cases and [Formula: see text] and [Formula: see text] ([Formula: see text], [Formula: see text], [Formula: see text]) for the healthy case.

Improved identifiability of myocardial material parameters by an energy-based cost function

Myocardial stiffness is a valuable clinical biomarker for the monitoring and stratification of heart failure (HF). Cardiac finite element models provide a biomechanical framework for the assessment of stiffness through the determination of the myocardial constitutive model parameters. The reported parameter intercorrelations in popular constitutive relations, however, obstruct the unique estimation of material parameters and limit the reliable translation of this stiffness metric to clinical practice. Focusing on the role of the cost function (CF) in parameter identifiability, we investigate the performance of a set of geometric indices (based on displacements, strains, cavity volume, wall thickness and apicobasal dimension of the ventricle) and a novel CF derived from energy conservation. Our results, with a commonly used transversely isotropic material model (proposed by Guccione et al.), demonstrate that a single geometry-based CF is unable to uniquely constrain the parameter space. The energy-based CF, conversely, isolates one of the parameters and in conjunction with one of the geometric metrics provides a unique estimation of the parameter set. This gives rise to a new methodology for estimating myocardial material parameters based on the combination of deformation and energetics analysis. The accuracy of the pipeline is demonstrated in silico, and its robustness in vivo, in a total of 8 clinical data sets (7 HF and one control). The mean identified parameters of the Guccione material law were [Formula: see text] and [Formula: see text] ([Formula: see text], [Formula: see text], [Formula: see text]) for the HF cases and [Formula: see text] and [Formula: see text] ([Formula: see text], [Formula: see text], [Formula: see text]) for the healthy case.

Simulating the effects of growth and fiber dispersion on the electromechanical response of a cardiac ventricular wedge affected from concentric hypertrophy

In this paper, we analyze the epicardial electromechanical response of an in silico cardiac ventricular wedge under both healthy and concentric hypertrophic conditions. This is achieved by taking into account the growth of the wedge thickness and the fiber dispersion that may follow. The electromechanical response is described in terms of some macroscopic measures, i.e. the action potential duration, the conduction velocity, the contractility and the contraction force. Our results suggest that growth reduces the action potential duration and conduction velocity, whilst it increases the contractility and contraction force, yielding an overall negative effect. In presence of fiber dispersion, the action potential duration and conduction velocity are not affected further, whilst the effect on the contractility and contraction force is enhanced.

Multi-scale Modeling of the Cardiovascular System: Disease Development, Progression, and Clinical Intervention

Cardiovascular diseases (CVDs) are the leading cause of death in the western world. With the current development of clinical diagnostics to more accurately measure the extent and specifics of CVDs, a laudable goal is a better understanding of the structure-function relation in the cardiovascular system. Much of this fundamental understanding comes from the development and study of models that integrate biology, medicine, imaging, and biomechanics. Information from these models provides guidance for developing diagnostics, and implementation of these diagnostics to the clinical setting, in turn, provides data for refining the models. In this review, we introduce multi-scale and multi-physical models for understanding disease development, progression, and designing clinical interventions. We begin with multi-scale models of cardiac electrophysiology and mechanics for diagnosis, clinical decision support, personalized and precision medicine in cardiology with examples in arrhythmia and heart failure. We then introduce computational models of vasculature mechanics and associated mechanical forces for understanding vascular disease progression, designing clinical interventions, and elucidating mechanisms that underlie diverse vascular conditions. We conclude with a discussion of barriers that must be overcome to provide enhanced insights, predictions, and decisions in pre-clinical and clinical applications.

Mathematical modeling of cardiac growth and remodeling

This review provides an overview of the current state of mathematical models of cardiac growth and remodeling (G&R). We concisely describe the experimental observations associated with cardiac G&R and discuss existing mathematical models that describe this process. To facilitate the discussion, we have organized the G&R models in terms of (1) the physical focus (biochemical vs mechanical) and (2) the process that they describe (myocyte hypertrophy vs extracellular matrix remodeling). The review concludes with a discussion of some possible directions that can advance the existing state of cardiac G&R mathematical modeling. WIREs Syst Biol Med 2016, 8:211-226. doi: 10.1002/wsbm.1330 For further resources related to this article, please visit the WIREs website.