Computational Investigation of Transmural Differences in Left Ventricular Contractility

Biventricular intraventricular mechanical and electrical dyssynchrony in pulmonary arterial hypertension

Background

Pulmonary arterial hypertension (PAH) leads to myocardial remodeling, manifesting as mechanical dyssynchrony (M-dys) and electrical dyssynchrony (E-dys), in both right (RV) and left ventricles (LV). However, the impacts of layer-specific intraventricular M-dys on biventricular functions and its association with E-dys in PAH remain unclear.

Methods

Seventy-nine newly diagnosed patients with PAH undergoing cardiac magnetic resonance scanning were consecutively recruited between January 2011 and December 2017. The biventricular volumetric and layer-specific intraventricular M-dys were analyzed. The QRS duration z-scores were calculated after adjusting for age and sex.

Results

77.22 % of patients were female (mean age 30.30 ± 9.79 years; median follow-up 5.53 years). Further, 29 (36.71 %) patients succumbed to all-cause mortality by the end of the study. At the baseline, LV layer-specific intraventricular M-dys had apparent transmural gradients compared with RV in the radial and circumferential directions. However, deceased patients lost the transmural gradients. The LV longitudinal strain rate time to late diastolic peak in the myocardial region (LVmyoLSRTTLDPintra) predicted long-term survival. The Kaplan-Meier curve revealed that patients with PAH with LVmyoLSRTTLDPintra <20.01 milliseconds had a worse prognosis. Larger right ventricle (RV) intraventricular M-dys resulted in worse RV ejection fraction. However, larger LV intraventricular M-dys in the late diastolic phase indicated remarkable exercise capacity and higher LV stroke volume index. E-dys and intraventricular M-dys had no direct correlations.

Conclusions

The layer-specific intraventricular M-dys had varying impacts on biventricular functions in PAH. PAH patients with LVmyoLSRTTLDPintra <20.01 milliseconds had a worse prognosis. LV intraventricular M-dys in the late diastolic phase needs more attention to precisely evaluate LV function.

Quantifying the microstructural and biomechanical changes in the porcine ventricles during growth and remodelling

Cardiac tissue growth and remodelling (G & R) occur in response to the changing physiological demands of the heart after birth. The early shift to pulmonary circulation produces an immediate increase in ventricular workload, causing microstructural and biomechanical changes that serve to maintain overall physiological homoeostasis. Such cardiac G & R continues throughout life. Quantifying the tissue's mechanical and microstructural changes because of G & R is of increasing interest, dovetailing with the emerging fields of personalised and precision solutions. This study aimed to determine equibiaxial, and non-equibiaxial extension, stress-relaxation, and the underlying microstructure of the passive porcine ventricles tissue at four time points spanning from neonatal to adulthood. The three-dimensional microstructure was investigated via two-photon excited fluorescence and second-harmonic generation microscopy on optically cleared tissues, describing the 3D orientation, rotation and dispersion of the cardiomyocytes and collagen fibrils. The results revealed that during biomechanical testing, myocardial ventricular tissue possessed non-linear, anisotropic, and viscoelastic behaviour. An increase in stiffness and viscoelasticity was noted for the left and right ventricular free walls from neonatal to adulthood. Microstructural analyses revealed concomitant increases in cardiomyocyte rotation and dispersion. This study provides baseline data, describing the biomechanical and microstructural changes in the left and right ventricular myocardial tissue during G & R, which should prove valuable to researchers in developing age-specific, constitutive models for more accurate computational simulations. STATEMENT OF SIGNIFICANCE: There is a dearth of experimental data describing the growth and remodelling of left and right ventricular tissue. The published literature is fragmented, with data reported via different experimental techniques using tissues harvested from a variety of animals, with different gender and ages. This prevents developing a continuum of data spanning birth to death, so limiting the potential that can be leveraged to aid computational modelling and simulations. In this study, equibiaxial, non-equibiaxial, and stress-relaxation data are presented, describing directional-dependent material responses. The biomechanical data is consolidated with equivalent microstructural data, an important element for the development of future material models. Combined, these data describe microstructural and biomechanical changes in the ventricles, spanning G &R from neonatal to adulthood.

Multiscale characterization of left ventricle active behavior in the mouse

The myocardium possesses an intricately designed microarchitecture to produce an optimal cardiac contraction. The contractile behavior of the heart is generated at the sarcomere level and travels across several length scales to manifest as the systolic function at the organ level. While passive myocardial behavior has been studied extensively, the translation of active tension produced at the fiber level to the organ-level function is not well understood. Alterations in cardiac systolic function are often key sequelae in structural heart diseases, such as myocardial infarction and systolic heart failure; thus, characterization of the contractile behavior of the heart across multiple length scales is essential to improve our understanding of mechanisms collectively leading to depressed systolic function. In this study, we present a methodology to characterize the active behavior of left ventricle free wall (LVFW) myocardial tissues in mice. Combined with active tests in papillary muscle fibers and conventional in vivo contractility measurement at the organ level in an animal-specific manner, we establish a multiscale active characterization of the heart from fiber to organ. In addition, we quantified myocardial architecture from histology to shed light on the directionality of the contractility at the tissue level. The LVFW tissue activation-relaxation behavior under isometric conditions was qualitatively similar to that of the papillary muscle fiber bundle. However, the maximum stress developed in the LVFW tissue was an order of magnitude lower than that developed by a fiber bundle, and the time taken for active forces to plateau was 2-3 orders of magnitude longer. Although the LVFW tissue exhibited a slightly stiffer passive response in the circumferential direction, the tissues produced significantly larger active stresses in the longitudinal direction during active testing. Also, contrary to passive viscoelastic stress relaxation, active stresses relaxed faster in the direction with larger peak stresses. The multiscale experimental pipeline presented in this work is expected to provide crucial insight into the contractile adaptation mechanisms of the heart with impaired systolic function. STATEMENT OF SIGNIFICANCE: Heart failure cause significant alterations to the contractile-relaxation behavior of the yocardium. Multiscale characterization of the contractile behavior of the myocardium is essential to advance our understanding of how contractility translates from fiber to organ and to identify the multiscale mechanisms leading to impaired cardiac function. While passive myocardial behavior has been studied extensively, the investigation of tissue-level contractile behavior remains critically scarce in the literature. To the best of our knowledge, our study here is the first to investigate the contractile behavior of the left ventricle at multiple length scales in small animals. Our results indicate that the active myocardial wall is a function of transmural depth and relaxes faster in the direction with larger peak stresses.

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.

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.

In-silico human electro-mechanical ventricular modelling and simulation for drug-induced pro-arrhythmia and inotropic risk assessment

Human-based computational modelling and simulation are powerful tools to accelerate the mechanistic understanding of cardiac patho-physiology, and to develop and evaluate therapeutic interventions. The aim of this study is to calibrate and evaluate human ventricular electro-mechanical models for investigations on the effect of the electro-mechanical coupling and pharmacological action on human ventricular electrophysiology, calcium dynamics, and active contraction. The most recent models of human ventricular electrophysiology, excitation-contraction coupling, and active contraction were integrated, and the coupled models were calibrated using human experimental data. Simulations were then conducted using the coupled models to quantify the effects of electro-mechanical coupling and drug exposure on electrophysiology and force generation in virtual human ventricular cardiomyocytes and tissue. The resulting calibrated human electro-mechanical models yielded active tension, action potential, and calcium transient metrics that are in agreement with experiments for endocardial, epicardial, and mid-myocardial human samples. Simulation results correctly predicted the inotropic response of different multichannel action reference compounds and demonstrated that the electro-mechanical coupling improves the robustness of repolarisation under drug exposure compared to electrophysiology-only models. They also generated additional evidence to explain the partial mismatch between in-silico and in-vitro experiments on drug-induced electrophysiology changes. The human calibrated and evaluated modelling and simulation framework constructed in this study opens new avenues for future investigations into the complex interplay between the electrical and mechanical cardiac substrates, its modulation by pharmacological action, and its translation to tissue and organ models of cardiac patho-physiology.

Relationship of Transmural Variations in Myofiber Contractility to Left Ventricular Ejection Fraction: Implications for Modeling Heart Failure Phenotype With Preserved Ejection Fraction

The pathophysiological mechanisms underlying preserved left ventricular (LV) ejection fraction (EF) in patients with heart failure and preserved ejection fraction (HFpEF) remain incompletely understood. We hypothesized that transmural variations in myofiber contractility with existence of subendocardial dysfunction and compensatory increased subepicardial contractility may underlie preservation of LVEF in patients with HFpEF. We quantified alterations in myocardial function in a mathematical model of the human LV that is based on the finite element method. The fiber-reinforced material formulation of the myocardium included passive and active properties. The passive material properties were determined such that the diastolic pressure-volume behavior of the LV was similar to that shown in published clinical studies of pressure-volume curves. To examine changes in active properties, we considered six scenarios: (1) normal properties throughout the LV wall; (2) decreased myocardial contractility in the subendocardium; (3) increased myocardial contractility in the subepicardium; (4) myocardial contractility decreased equally in all layers, (5) myocardial contractility decreased in the midmyocardium and subepicardium, (6) myocardial contractility decreased in the subepicardium. Our results indicate that decreased subendocardial contractility reduced LVEF from 53.2 to 40.5%. Increased contractility in the subepicardium recovered LVEF from 40.5 to 53.2%. Decreased contractility transmurally reduced LVEF and could not be recovered if subepicardial and midmyocardial contractility remained depressed. The computational results simulating the effects of transmural alterations in the ventricular tissue replicate the phenotypic patterns of LV dysfunction observed in clinical practice. In particular, data for LVEF, strain and displacement are consistent with previous clinical observations in patients with HFpEF, and substantiate the hypothesis that increased subepicardial contractility may compensate for subendocardial dysfunction and play a vital role in maintaining LVEF.

Differential Effects of Isoproterenol on Regional Myocardial Mechanics in Rat using 3D cine DENSE Cardiovascular Magnetic Resonance

The present study assessed the acute effects of isoproterenol on left ventricular (LV) mechanics in healthy rats with the hypothesis that ß-adrenergic stimulation influences the mechanics of different myocardial regions of the LV wall in different ways. To accomplish this, magnetic resonance images were obtained in the LV of healthy rats with or without isoproterenol infusion. The LV contours were divided into basal, mid-ventricular, and apical regions. Additionally, the mid-ventricular myocardium was divided into three transmural layers with each layer partitioned into four segments (i.e., septal, inferior, lateral, and anterior). Peak systolic strains and torsion were quantified for each region. Isoproterenol significantly increased peak systolic radial strain and circumferential-longitudinal shear strain, as well as ventricular torsion, throughout the basal, mid-ventricle, and apical regions. In the mid-ventricle, isoproterenol significantly increased peak systolic radial strain, and induced significant increases in peak systolic circumferential strain and longitudinal strain in the septum. Isoproterenol consistently increased peak systolic circumferential-longitudinal shear strain in all mid-ventricular segments. Ventricular torsion was significantly increased in nearly all segments except the inferior sub-endocardium. The effects of isoproterenol on LV systolic mechanics (i.e., 3D strains and torsion) in healthy rats depend on the region. This region-dependency is also strain component-specific. These results provide insight into the regional response of LV mechanics to ß-adrenergic stimulation in rats, and could act as a baseline for future studies on subclinical abnormalities associated with the inotropic response in heart disease.

Evaluation of a Novel Finite Element Model of Active Contraction in the Heart

Finite element (FE) modeling is becoming a widely used approach for the investigation of global heart function. In the present study, a novel model of cellular-level systolic contraction, which includes both length- and velocity-dependence, was implemented into a 3D non-linear FE code. To validate this new FE implementation, an optimization procedure was used to determine the contractile parameters, associated with sarcomeric function, by comparing FE-predicted pressure and strain to experimental measures collected with magnetic resonance imaging and catheterization in the ventricles of five healthy rats. The pressure-volume relationship generated by the FE models matched well with the experimental data. Additionally, the regional distribution of end-systolic strains and circumferential-longitudinal shear angle exhibited good agreement with experimental results overall, with the main deviation occurring in the septal region. Moreover, the FE model predicted a heterogeneous distribution of sarcomere re-lengthening after ventricular ejection, which is consistent with previous in vivo studies. In conclusion, the new FE active contraction model was able to predict the global performance and regional mechanical behaviors of the LV during the entire cardiac cycle. By including more accurate cellular-level mechanisms, this model could provide a better representation of the LV and enhance cardiac research related to both systolic and diastolic dysfunction.