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

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.

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.

Mechanisms of coronary sinus reducer for treatment of myocardial ischemia: in silico study

The coronary sinus reducer (CSR) is an emerging medical device for treating patients with refractory angina, often associated with myocardial ischemia. Patients implanted with CSR have shown positive outcomes, but the underlying mechanisms are unclear. This study sought to understand the mechanisms of CSR by investigating its effects on coronary microcirculation hemodynamics that may help explain the therapy's efficacy. We applied a validated computer model of the coronary microcirculation to investigate how CSR affects hemodynamics under different degrees of coronary artery stenosis. With moderate coronary stenosis, an increase in capillary transit time (CTT) [up to 69% with near-complete coronary sinus (CS) occlusion] is the key change associated with CSR. Because capillaries in the microcirculation can still receive oxygenated blood from the upstream artery with moderate stenosis, the increase in CTT allows more time for the exchange of gases and nutrients, aiding tissue oxygenation. With severe coronary stenosis; however, the redistribution of blood draining from the nonischemic region to the ischemic region (up to 96% with near-complete CS occlusion) and the reduction in capillary flow heterogeneity are the key changes associated with CSR. Because blood draining from the nonischemic region is not completely devoid of O2, the redistribution of blood to the capillaries in the ischemic region by CSR is beneficial especially when little or no oxygenated blood reaches these capillaries. This simulation study provides insights into the mechanisms of CSR in improving clinical symptoms. The mechanisms differ with the severity of the upstream stenosis.NEW & NOTEWORTHY Emerging coronary venous retroperfusion treatments, particularly coronary sinus reducer (CSR) for refractory angina linked to myocardial ischemia, show promise; however, their mechanisms of action are not well understood. We find that CSR's effectiveness varies with the severity of coronary stenosis. In moderate stenosis, CSR improves tissue oxygenation by increasing capillary transit time, whereas in severe stenosis, it redistributes blood from nonischemic to ischemic regions and reduces capillary flow heterogeneity.

Systolic blood pressure ≤110 mm Hg is associated with severe coronary microvascular ischemia and higher risk for ventricular arrhythmias in hypertrophic cardiomyopathy

Background

Coronary microvascular dysfunction (CMD) and hypertension (HTN) occur frequently in hypertrophic cardiomyopathy (HCM), but whether blood pressure (BP) influences CMD and outcomes is unknown.

Objective

The purpose of this study was to test the hypothesis that HTN is associated with worse CMD and outcomes.

Methods

This retrospective study included 690 HCM patients. All patients underwent cardiac magnetic resonance imaging, echocardiography, and rhythm monitoring; 127 patients also underwent rest/vasodilator stress 13NH3 positron emission tomography myocardial perfusion imaging. Patients were divided into 3 groups based on their rest systolic blood pressure (SBP) (group 1 ≤110 mm Hg; group 2 111-140; group 3 >140 mm Hg) and were followed for development of ventricular tachycardia (VT)/ventricular fibrillation (VF), heart failure (HF), death, and composite outcome.

Results

Group 1 patients had the lowest age and left ventricular (LV) mass but the highest prevalence of nonobstructive hemodynamics and restrictive diastolic filling. LV scar was similar in the 3 groups. Group 1 had the lowest rest and stress myocardial blood flow (MBF) and highest SDS (summed difference score). Rest SBP was positively correlated with stress MBF and negatively correlated with SDS. Group 1 had the highest incidence of VT/VF, whereas the incidences of HF, death, and composite outcome were similar among the 3 groups. In multivariate analysis, rest SBP ≤110 mm Hg was independently associated with VT/VF (hazard ratio 2.6; 95% confidence interval 1.0-6.7; P = .04).

Conclusion

SBP ≤110 mm Hg is associated with greater severity of CMD and coronary microvascular ischemia and higher incidence of ventricular arrhythmias in HCM.

Computational analysis of ventricular mechanics in hypertrophic cardiomyopathy patients

Hypertrophic cardiomyopathy (HCM) is a genetic heart disease that is associated with many pathological features, such as a reduction in global longitudinal strain (GLS), myofiber disarray and hypertrophy. The effects of these features on left ventricle (LV) function are, however, not clear in two phenotypes of HCM, namely, obstructive and non-obstructive. To address this issue, we developed patient-specific computational models of the LV using clinical measurements from 2 female HCM patients and a control subject. Left ventricular mechanics was described using an active stress formulation and myofiber disarray was described using a structural tensor in the constitutive models. Unloaded LV configuration for each subject was first determined from their respective end-diastole LV geometries segmented from the cardiac magnetic resonance images, and an empirical single-beat estimation of the end-diastolic pressure volume relationship. The LV was then connected to a closed-loop circulatory model and calibrated using the clinically measured LV pressure and volume waveforms, peak GLS and blood pressure. Without consideration of myofiber disarray, peak myofiber tension was found to be lowest in the obstructive HCM subject (60 kPa), followed by the non-obstructive subject (242 kPa) and the control subject (375 kPa). With increasing myofiber disarray, we found that peak tension has to increase in the HCM models to match the clinical measurements. In the obstructive HCM patient, however, peak tension was still depressed (cf. normal subject) at the largest degree of myofiber disarray found in the clinic. The computational modeling workflow proposed here can be used in future studies with more HCM patient data.

The dependency of fetal left ventricular biomechanics function on myocardium helix angle configuration

The helix angle configuration of the myocardium is understood to contribute to the heart function, as finite element (FE) modeling of postnatal hearts showed that altered configurations affected cardiac function and biomechanics. However, similar investigations have not been done on the fetal heart. To address this, we performed image-based FE simulations of fetal left ventricles (LV) over a range of helix angle configurations, assuming a linear variation of helix angles from epicardium to endocardium. Results showed that helix angles have substantial influence on peak myofiber stress, cardiac stroke work, myocardial deformational burden, and spatial variability of myocardial strain. A good match between LV myocardial strains from FE simulations to those measured from 4D fetal echo images could only be obtained if the transmural variation of helix angle was generally between 110 and 130°, suggesting that this was the physiological range. Experimentally discovered helix angle configurations from the literature were found to produce high peak myofiber stress, high cardiac stroke work, and a low myocardial deformational burden, but did not coincide with configurations that would optimize these characteristics. This may suggest that the fetal development of myocyte orientations depends concurrently on several factors rather than a single factor. We further found that the shape, rather than the size of the LV, determined the manner at which helix angles influenced these characteristics, as this influence changed significantly when the LV shape was varied, but not when a heart was scaled from fetal to adult size while retaining the same shape. This may suggest that biomechanical optimality would be affected during diseases that altered the geometric shape of the LV.

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.

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.