Cardiac resynchronization therapy (CRT) is an established treatment for left bundle branch block (LBBB) resulting in mechanical dyssynchrony. Approximately 1/3 of patients with CRT, however, are ...non-responders. To understand factors affecting CRT response, an electromechanics-perfusion computational model based on animal-specific left ventricular (LV) geometry and coronary vascular networks located in the septum and LV free wall is developed. The model considers contractility-flow and preload-activation time relationships, and is calibrated to simultaneously match the experimental measurements in terms of the LV pressure, volume waveforms and total coronary flow in the left anterior descending and left circumflex territories from 2 swine models under right atrium and right ventricular pacing. The model is then applied to investigate the responses of CRT indexed by peak LV pressure and (dP/dt)max at multiple pacing sites with different degrees of perfusion in the LV free wall. Without the presence of ischemia, the model predicts that basal-lateral endocardial region is the optimal pacing site that can best improve (dP/dt)max by 20%, and is associated with the shortest activation time. In the presence of ischemia, a non-ischemic region becomes the optimal pacing site when coronary flow in the ischemic region fell below 30% of its original value. Pacing at the ischemic region produces little response at that perfusion level. The optimal pacing site is associated with one that optimizes the LV activation time. These findings suggest that CRT response is affected by both pacing site and coronary perfusion, which may have clinical implication in improving CRT responder rates.
•A novel electromechanics-perfusion computational model based on animal-specific geometry and measurements is developed.•The relationship between contractility and myocardial perfusion is considered.•The calibrated model is applied to optimize cardiac resynchronization therapy (CRT) response.•The effects of pacing site and degree of coronary perfusion on CRT response are investigated.
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
Overview of cardiac–coronary interaction mechanisms. M→C and C→M denote effects of myocardium on coronary vasculature and vice versa, respectively. M↔C denotes bi‐directional cardiac–coronary interactions.
Left ventricular assist device (LVAD) is associated with a high incidence of right ventricular (RV) failure, which is hypothesized to be caused by the occurring inter-ventricular interactions when ...the LV is unloaded. Factors contributing to these interactions are unknown.
We used computer modeling to investigate the impact of the HeartMate 3 LVAD on RV functions. The model was first calibrated against pressure-volume (PV) loops associated with a heart failure (HF) patient and validated against measurements of inter-ventricular interactions in animal experiments. The model was then applied to investigate the effects of LVAD on (1) RV chamber contractility indexed by
derived from its end-systolic PV relationship, and (2) RV diastolic function indexed by
derived from its end-diastolic PV relationship. We also investigated how septal wall thickness and regional contractility affect the impact of LVAD on RV function.
The impact of LVAD on RV chamber contractility is small at a pump speed lower than 4k rpm. At a higher pump speed between 4k and 9k rpm, however, RV chamber contractility is reduced (by ~3% at 6k rpm and ~10% at 9k rpm). The reduction of RV chamber contractility is greater with a thinner septal wall or with a lower myocardial contractility at the LV free wall, septum, or RV free wall.
RV chamber contractility is reduced at a pump speed higher than 4k rpm, and this reduction is greater with a thinner septal wall or lower regional myocardial contractility. Findings here may have clinical implications in identifying LVAD patients who may suffer from RV failure.
Abstract Even when entirely unloaded, biological structures are not stress-free, as shown by Y.C. Fung׳s seminal opening angle experiment on arteries and the left ventricle. As a result of this ...prestrain, subject-specific geometries extracted from medical imaging do not represent an unloaded reference configuration necessary for mechanical analysis, even if the structure is externally unloaded. Here we propose a new computational method to create physiological residual stress fields in subject-specific left ventricular geometries using the continuum theory of fictitious configurations combined with a fixed-point iteration. We also reproduced the opening angle experiment on four swine models, to characterize the range of normal opening angle values. The proposed method generates residual stress fields which can reliably reproduce the range of opening angles between 8.7±1.8 and 16.6±13.7 as measured experimentally. We demonstrate that including the effects of prestrain reduces the left ventricular stiffness by up to 40%, thus facilitating the ventricular filling, which has a significant impact on cardiac function. This method can improve the fidelity of subject-specific models to improve our understanding of cardiac diseases and to optimize treatment options.
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.
Critical aortic stenosis (AS) of the fetal heart causes a drastic change in the cardiac biomechanical environment. Consequently, a substantial proportion of such cases will lead to a ...single-ventricular birth outcome. However, the biomechanics of the disease is not well understood. To address this, we performed Finite Element (FE) modelling of the healthy fetal left ventricle (LV) based on patient-specific 4D ultrasound imaging, and simulated various disease features observed in clinical fetal AS to understand their biomechanical impact. These features included aortic stenosis, mitral regurgitation (MR) and LV hypertrophy, reduced contractility, and increased myocardial stiffness. AS was found to elevate LV pressures and myocardial stresses, and depending on severity, can drastically decrease stroke volume and myocardial strains. These effects are moderated by MR. AS alone did not lead to MR velocities above 3 m/s unless LV hypertrophy was included, suggesting that hypertrophy may be involved in clinical cases with high MR velocities. LV hypertrophy substantially elevated LV pressure, valve flow velocities and stroke volume, while reducing LV contractility resulted in diminished LV pressure, stroke volume and wall strains. Typical extent of hypertrophy during fetal AS in the clinic, however, led to excessive LV pressure and valve velocity in the FE model, suggesting that reduced contractility is typically associated with hypertrophy. Increased LV passive stiffness, which might represent fibroelastosis, was found to have minimal impact on LV pressures, stroke volume, and wall strain. This suggested that fibroelastosis could be a by-product of the disease progression and does not significantly impede cardiac function. Our study demonstrates that FE modelling is a valuable tool for elucidating the biomechanics of congenital heart disease and can calculate parameters which are difficult to measure, such as intraventricular pressure and myocardial stresses.
Finite element (FE) modeling is becoming increasingly prevalent in the world of cardiac mechanics; however, many existing FE models are phenomenological and thus do not capture cellular-level ...mechanics. This work implements a cellular-level contraction scheme into an existing nonlinear FE code to model ventricular contraction. Specifically, this contraction model incorporates three myosin states: OFF-, ON-, and an attached force-generating state. It has been speculated that force-dependent transitions from the OFF- to ON-state may contribute to length-dependent activation at the cellular level. The current work investigates the contribution of force-dependent recruitment out of the OFF-state to ventricular-level function, specifically the Frank–Starling relationship, as seen through the end-systolic pressure–volume relationship (ESPVR). Five FE models were constructed using geometries of rat left ventricles obtained via cardiac magnetic resonance imaging. FE simulations were conducted to optimize parameters for the cellular contraction model such that the differences between FE predicted ventricular pressures for the models and experimentally measured pressures were minimized. The models were further validated by comparing FE predicted end-systolic strain to experimentally measured strain. Simulations mimicking vena cava occlusion generated descending pressure volume loops from which ESPVRs were calculated. In simulations with the inclusion of the OFF-state, using a force-dependent transition to the ON-state, the ESPVR calculated was steeper than in simulations excluding the OFF-state. Furthermore, the ESPVR was also steeper when compared to models that included the OFF-state without a force-dependent transition. This suggests that the force-dependent recruitment of thick filament heads from the OFF-state at the cellular level contributes to the Frank–Starling relationship observed at the organ level.
Abstract Injection of biomaterials into diseased myocardium has been associated with decreased myofiber stress, restored left ventricular (LV) geometry and improved LV function. However, its exact ...mechanism(s) of action remained unclear. In this work, we present the first patient-specific computational model of biomaterial injection that accounts for the possibility of residual strain and stress introduced by this treatment. We show that the presence of residual stress can create more heterogeneous regional myofiber stress and strain fields. Our simulation results show that the treatment generates low stress and stretch areas between injection sites, and high stress and stretch areas between the injections and both the endocardium and epicardium. Globally, these local changes are translated into an increase in average myofiber stress and its standard deviation (from 6.9±4.6 to 11.2±48.8 kPa and 30±15 to 35.1±50.9 kPa at end-diastole and end-systole, respectively). We also show that the myofiber stress field is sensitive to the void-to-size ratio. For a constant void size, the myofiber stress field became less heterogeneous with decreasing injection volume. These results suggest that the residual stress and strain possibly generated by biomaterial injection treatment can have large effects on the regional myocardial stress and strain fields, which may be important in the remodeling process.
Individually personalized computational models of heart mechanics can be used to estimate important physiological and clinically‐relevant quantities that are difficult, if not impossible, to directly ...measure in the beating heart. Here, we present a novel and efficient framework for creating patient‐specific biventricular models using a gradient‐based data assimilation method for evaluating regional myocardial contractility and estimating myofiber stress. These simulations can be performed on a regular laptop in less than 2 h and produce excellent fit between measured and simulated volume and strain data through the entire cardiac cycle. By applying the framework using data obtained from 3 healthy human biventricles, we extracted clinically important quantities as well as explored the role of fiber angles on heart function. Our results show that steep fiber angles at the endocardium and epicardium are required to produce simulated motion compatible with measured strain and volume data. We also find that the contraction and subsequent systolic stresses in the right ventricle are significantly lower than that in the left ventricle. Variability of the estimated quantities with respect to both patient data and modeling choices are also found to be low. Because of its high efficiency, this framework may be applicable to modeling of patient specific cardiac mechanics for diagnostic purposes.
Here we present a novel and efficient framework for creating patient‐specific biventricular models using a gradient‐based data assimilation method for evaluating regional myocardial contractility and estimating myofiber stress. By applying the framework using data obtained from 3 healthy human biventricles, we extracted clinically important quantities and explored the role of fiber angles on heart function. Because of its high efficiency, this framework may be applicable to modeling patient‐specific cardiac mechanics for diagnostic purposes.