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

Comparison of Left Ventricular Function Derived from Subject-Specific Inverse Finite Element Modeling Based on 3D ECHO and Magnetic Resonance Images

Three-dimensional echocardiography (3D ECHO) and magnetic resonance (MR) imaging are frequently used in patients and animals to evaluate heart functions. Inverse finite element (FE) modeling is increasingly applied to MR images to quantify left ventricular (LV) function and estimate myocardial contractility and other cardiac biomarkers. It remains unclear, however, as to whether myocardial contractility derived from the inverse FE model based on 3D ECHO images is comparable to that derived from MR images. To address this issue, we developed a subject-specific inverse FE model based on 3D ECHO and MR images acquired from seven healthy swine models to investigate if there are differences in myocardial contractility and LV geometrical features derived using these two imaging modalities. We showed that end-systolic and end-diastolic volumes derived from 3D ECHO images are comparable to those derived from MR images (R2=0.805 and 0.969, respectively). As a result, ejection fraction from 3D ECHO and MR images are linearly correlated (R2=0.977) with the limit of agreement (LOA) ranging from -17.95% to 45.89%. Using an inverse FE modeling to fit pressure and volume waveforms in subject-specific LV geometry reconstructed from 3D ECHO and MR images, we found that myocardial contractility derived from these two imaging modalities are linearly correlated with an R2 value of 0.989, a gradient of 0.895, and LOA ranging from -6.11% to 36.66%. This finding supports using 3D ECHO images in image-based inverse FE modeling to estimate myocardial contractility.

Preconditioning with selective autoretroperfusion: In vivo and in silico evidence of washout hypothesis

Introduction

Prompt reperfusion of coronary artery after acute myocardial infarction (AMI) is crucial for minimizing heart injury. The myocardium, however, may experience additional injury due to the flow restoration itself (reperfusion injury, RI). The purpose of this study was to demonstrate that short preconditioning (10 min) with selective autoretroperfusion (SARP) ameliorates RI, based on a washout hypothesis.

Methods

AMI was induced in 23 pigs (3 groups) by occluding the left anterior descending (LAD) artery. In SARP-b (SARP balloon inflated) and SARP-nb (SARP balloon deflated) groups, arterial blood was retroperfused for 10 min via the great cardiac vein before releasing the arterial occlusion. A mathematical model of coronary circulation was used to simulate the SARP process and evaluate the potential washout effect.

Results

SARP restored left ventricular function during LAD occlusion. Ejection fraction in the SARP-b group returned to baseline levels, compared to SARP-nb and control groups. Infarct area was significantly larger in the control group than in the SARP-b and SARP-nb groups. End-systolic wall thickness was preserved in the SARP-b compared to the SARP-nb and control groups. Analyte values (pH, lactate, glucose, and others), measured every 2 min during retroperfusion, suggest a "washout" effect as one important mechanism of action of SARP in reducing infarct size. With SARP, the values progressively approached baseline levels. The mathematical model also confirmed a possible washout effect of tracers.

Discussion

RI can be ameliorated by delaying restoration of arterial flow for a brief period of time while pretreating the infarction with SARP to restore homeostasis via a washout mechanism.

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.

Investigation into the two-way interaction of coronary flow and heart function in coronary artery disease predicted by a computational model of autoregulation of coronary flow

This study presents a closed-loop computational model to investigate the interplay between heart function, coronary flow, and systemic circulation during exercise, with a specific focus on the impact of coronary artery stenosis. The model incorporates a lumped representation of the heart, main arteries, and coronary arteries, establishing a closed circulatory system. The simulation investigates the autoregulation of coronary flow in response to myocardial oxygen demands during physical exercise by incorporating sympathetic and parasympathetic functions. This study establishes a closed supply-demand loop and investigates the effect of coronary flow deficiency on heart function and systemic circulation in coronary artery diseases during exercise. In coronary artery diseases with low stenosis, heart function and systemic flow resemble those of a healthy person. However, as stenosis intensifies with physical exercise, an additional regulatory mechanism (reg2) is activated. This mechanism adjusts coronary flow by reducing myocardial contractility (E) and increasing heart rate (HR) while maintaining cardiac output (CO). The study results indicate that, at the highest exercise intensity for a healthy individual (HR = 150), the value of E increases from 6 to 8.65mmHg/ml. Meanwhile, for a patient with 85 % coronary artery stenosis in the same exercise intensity, the HR increases to 200, and the value of E decreases to 3.45mmHg/ml. The results also demonstrate that the initiation of the (reg2) mechanism at rest occurs at 83 % stenosis, while at the highest exercise intensity, this mechanism commences at 67 % stenosis.

Simulation of coronary capillary transit time based on full vascular model of the heart

Capillary transit time (CTT) is a fundamental determinant of gas exchange between blood and tissues in the heart and other organs. Despite advances in experimental techniques, it remains difficult to measure coronary CTT in vivo. Here, we developed a novel computational framework that couples coronary microcirculation with cardiac mechanics in a closed-loop system that enables prediction of hemodynamics in the entire coronary network, including arteries, veins, and capillaries. We also developed a novel "particle-tracking" approach for computing CTT where "virtual tracers" are individually tracked as they traverse the capillary network. Model predictions compare well with blood pressure and flow rate distributions in the arterial network reported in previous studies. Model predictions of transit times in the capillaries (1.21 ± 1.5 s) and entire coronary network (11.8 ± 1.8 s) also agree with measurements. We show that, with increasing coronary artery stenosis (as quantified by fractional flow reserve, FFR), intravascular pressure and flow rate downstream are reduced but remain non-stationary even at 100 % stenosis because some flow (∼3 %) is redistributed from the non-occluded to the occluded territories. Importantly, the model predicts that occlusion of a large artery results in higher CTT. For moderate stenosis (FFR > 0.6), the increase in CTT (from 1.21 s without stenosis to 2.23 s at FFR=0.6) is caused by a decrease in capillary flow rate. In severe stenosis (FFR = 0.1), the increase in CTT to 14.2 s is due to both a decrease in flow rate and an increase in path length taken by "virtual tracers" in the capillary network.

An in silico study of the effects of left ventricular assist device on right ventricular function and inter-ventricular interaction

BACKGROUND

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.

METHODS

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 byV 60 derived from its end-systolic PV relationship, and (2) RV diastolic function indexed byV 20 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.

RESULTS

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.

CONCLUSION

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.

Myocardial perfusion and flow reserve in the asynchronous heart: mechanistic insight from a computational model

The tight coupling between myocardial oxygen demand and supply has been recognized for decades, but it remains controversial whether this coupling persists under asynchronous activation, such as during left bundle branch block (LBBB). Furthermore, it is unclear whether the amount of local cardiac wall growth, following longer-lasting asynchronous activation, can explain differences in myocardial perfusion distribution between subjects. For a better understanding of these matters, we built upon our existing modeling framework for cardiac mechanics-to-perfusion coupling by incorporating coronary autoregulation. Regional coronary flow was regulated with a vasodilator signal based on regional demand, as estimated from regional fiber stress-strain area. Volume of left ventricular wall segments was adapted with chronic asynchronous activation toward a homogeneous distribution of myocardial oxygen demand per tissue weight. Modeling results show that 1) both myocardial oxygen demand and supply are decreased in early activated regions and increased in late-activated regions; 2) but that regional hyperemic flow remains unaffected; while 3) regional myocardial flow reserve (the ratio of hyperemic to resting myocardial flow) decreases with increases in absolute regional myocardial oxygen demand as well as with decreases in wall thickness. These findings suggest that septal hypoperfusion in LBBB represents an autoregulatory response to reduced myocardial oxygen demand. Furthermore, oxygen demand-driven remodeling of wall mass can explain asymmetric hypertrophy and the related homogenization of myocardial perfusion and flow reserve. Finally, the inconsistent observations of myocardial perfusion distribution can primarily be explained by the degree of dyssynchrony, the degree of asymmetric hypertrophy, and the imaging modality used.NEW & NOTEWORTHY This versatile modeling framework couples myocardial oxygen demand to oxygen supply and myocardial growth, enabling simulation of resting and hyperemic myocardial flow during acute and chronic asynchronous ventricular activation. Model-based findings suggest that reported inconsistencies in myocardial perfusion and flow reserve responses with asynchronous ventricular activation between patients can primarily be explained by the degree of dyssynchrony and wall mass remodeling, which together determine the heterogeneity in regional oxygen demand and, hence, supply with autoregulation.

Mechanism of exercise intolerance in heart diseases predicted by a computer model of myocardial demand-supply feedback system

BACKGROUND AND OBJECTIVE

The myocardial demand-supply feedback system plays an important role in augmenting blood supply in response to exercise-induced increased myocardial demand. During this feedback process, the myocardium and coronary blood flow interact bidirectionally at many different levels.

METHODS

To investigate these interactions, a novel computational framework that considers the closed myocardial demand-supply feedback system was developed. In the framework coupling the systemic circulation of the left ventricle and coronary perfusion with regulation, myocardial work affects coronary perfusion via flow regulation mechanisms (e.g., metabolic regulation) and myocardial-vessel interactions, whereas coronary perfusion affects myocardial contractility in a closed feedback system. The framework was calibrated based on the measurements from healthy subjects under graded exercise conditions, and then was applied to simulate the effects of graded exercise on myocardial demand-supply under different physiological and pathological conditions.

RESULTS

We found that the framework can recapitulate key features found during exercise in clinical and animal studies. We showed that myocardial blood flow is increased but maximum hyperemia is reduced during exercise, which led to a reduction in coronary flow reserve. For coronary stenosis and myocardial inefficiency, the model predicts that an increase in heart rate is necessary to maintain the baseline cardiac output. Correspondingly, the resting coronary flow reserve is exhausted and the range of heart rate before exhaustion of coronary flow reserve is reduced. In the presence of metabolic regulation dysfunction, the model predicts that the metabolic vasodilator signal is higher at rest, saturates faster during exercise, and as a result, causes quicker exhaustion of coronary flow reserve.

CONCLUSIONS

Model predictions showed that the coronary flow reserve deteriorates faster during graded exercise, which in turn, suggests a decrease in exercise tolerance for patients with stenosis, myocardial inefficiency and metabolic flow regulation dysfunction. The findings in this study may have clinical implications in diagnosing cardiovascular diseases.

A Closed-Loop Modeling Framework for Cardiac-to-Coronary Coupling

The mechanisms by which cardiac mechanics effect coronary perfusion (cardiac-to-coronary coupling) remain incompletely understood. Several coronary models have been proposed to deepen our understanding of coronary hemodynamics, but possibilities for in-depth studies on cardiac-to-coronary coupling are limited as mechanical properties like myocardial stress and strain are most often neglected. To overcome this limitation, a mathematical model of coronary mechanics and hemodynamics was implemented in the previously published multi-scale CircAdapt model of the closed-loop cardiovascular system. The coronary model consisted of a relatively simple one-dimensional network of the major conduit arteries and veins as well as a lumped parameter model with three transmural layers for the microcirculation. Intramyocardial pressure was assumed to arise from transmission of ventricular cavity pressure into the myocardial wall as well as myocardial stiffness, based on global pump mechanics and local myofiber mechanics. Model-predicted waveforms of global epicardial flow velocity, as well as of intramyocardial flow and diameter were qualitatively and quantitatively compared with reported data. Versatility of the model was demonstrated in a case study of aortic valve stenosis. The reference simulation correctly described the phasic pattern of coronary flow velocity, arterial flow impediment, and intramyocardial differences in coronary flow and diameter. Predicted retrograde flow during early systole in aortic valve stenosis was in agreement with measurements obtained in patients. In conclusion, we presented a powerful multi-scale modeling framework that enables realistic simulation of coronary mechanics and hemodynamics. This modeling framework can be used as a research platform for in-depth studies of cardiac-to-coronary coupling, enabling study of the effect of abnormal myocardial tissue properties on coronary hemodynamics.

Optimization of cardiac resynchronization therapy based on a cardiac electromechanics-perfusion computational model

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.

Transmural Distribution of Coronary Perfusion and Myocardial Work Density Due to Alterations in Ventricular Loading, Geometry and Contractility

Myocardial supply changes to accommodate the variation of myocardial demand across the heart wall to maintain normal cardiac function. A computational framework that couples the systemic circulation of a left ventricular (LV) finite element model and coronary perfusion in a closed loop is developed to investigate the transmural distribution of the myocardial demand (work density) and supply (perfusion) ratio. Calibrated and validated against measurements of LV mechanics and coronary perfusion, the model is applied to investigate changes in the transmural distribution of passive coronary perfusion, myocardial work density, and their ratio in response to changes in LV contractility, preload, afterload, wall thickness, and cavity volume. The model predicts the following: (1) Total passive coronary flow varies from a minimum value at the endocardium to a maximum value at the epicardium transmurally that is consistent with the transmural distribution of IMP; (2) Total passive coronary flow at different transmural locations is increased with an increase in either contractility, afterload, or preload of the LV, whereas is reduced with an increase in wall thickness or cavity volume; (3) Myocardial work density at different transmural locations is increased transmurally with an increase in either contractility, afterload, preload or cavity volume of the LV, but is reduced with an increase in wall thickness; (4) Myocardial work density-perfusion mismatch ratio at different transmural locations is increased with an increase in contractility, preload, wall thickness or cavity volume of the LV, and the ratio is higher at the endocardium than the epicardium. These results suggest that an increase in either contractility, preload, wall thickness, or cavity volume of the LV can increase the vulnerability of the subendocardial region to ischemia.