Cryo-EM structure of the human cardiac myosin filament.

Pumping of the heart is powered by filaments of the motor protein myosin that pull on actin filaments to generate cardiac contraction. In addition to myosin, the filaments contain cardiac myosin-binding protein C (cMyBP-C), which modulates contractility in response to physiological stimuli, and titin, which functions as a scaffold for filament assembly1. Myosin, cMyBP-C and titin are all subject to mutation, which can lead to heart failure. Despite the central importance of cardiac myosin filaments to life, their molecular structure has remained a mystery for 60 years2. Here we solve the structure of the main (cMyBP-C-containing) region of the human cardiac filament using cryo-electron microscopy. The reconstruction reveals the architecture of titin and cMyBP-C and shows how myosin's motor domains (heads) form three different types of motif (providing functional flexibility), which interact with each other and with titin and cMyBP-C to dictate filament architecture and function. The packing of myosin tails in the filament backbone is also resolved. The structure suggests how cMyBP-C helps to generate the cardiac super-relaxed state3; how titin and cMyBP-C may contribute to length-dependent activation4; and how mutations in myosin and cMyBP-C might disturb interactions, causing disease5,6. The reconstruction resolves past uncertainties and integrates previous data on cardiac muscle structure and function. It provides a new paradigm for interpreting structural, physiological and clinical observations, and for the design of potential therapeutic drugs.

Integrated multi-omic characterization of congenital heart disease.

The heart, the first organ to develop in the embryo, undergoes complex morphogenesis that when defective results in congenital heart disease (CHD). With current therapies, more than 90% of patients with CHD survive into adulthood, but many suffer premature death from heart failure and non-cardiac causes1. Here, to gain insight into this disease progression, we performed single-nucleus RNA sequencing on 157,273 nuclei from control hearts and hearts from patients with CHD, including those with hypoplastic left heart syndrome (HLHS) and tetralogy of Fallot, two common forms of cyanotic CHD lesions, as well as dilated and hypertrophic cardiomyopathies. We observed CHD-specific cell states in cardiomyocytes, which showed evidence of insulin resistance and increased expression of genes associated with FOXO signalling and CRIM1. Cardiac fibroblasts in HLHS were enriched in a low-Hippo and high-YAP cell state characteristic of activated cardiac fibroblasts. Imaging mass cytometry uncovered a spatially resolved perivascular microenvironment consistent with an immunodeficient state in CHD. Peripheral immune cell profiling suggested deficient monocytic immunity in CHD, in agreement with the predilection in CHD to infection and cancer2. Our comprehensive phenotyping of CHD provides a roadmap towards future personalized treatments for CHD.

Cardiac length-dependent activation driven by force-dependent thick-filament dynamics.

The length-dependent activation (LDA) of maximum force and calcium sensitivity are established features of cardiac muscle contraction but the dominant underlying mechanisms remain to be fully clarified. Alongside the well-documented regulation of contraction via the thin filaments, experiments have identified an additional force-dependent thick-filament activation, whereby myosin heads parked in a so-called off state become available to generate force. This process produces a feedback effect that may potentially drive LDA. Using biomechanical modeling of a human left-ventricular myocyte, this study investigates the extent to which the off-state dynamics could, by itself, plausibly account for LDA, depending on the specific mathematical formulation of the feedback. We hypothesized four different models of the off-state regulatory feedback based on (A) total force, (B) active force, (C) sarcomere strain, and (D) passive force. We tested if these models could reproduce the isometric steady-state and dynamic LDA features predicted by an earlier published model of a human left-ventricle myocyte featuring purely phenomenological length dependences. The results suggest that only total-force feedback (A) is capable of reproducing the expected behaviors, but that passive tension could provide a length-dependent signal on which to initiate the feedback. Furthermore, by attributing LDA to off-state dynamics, our proposed model also qualitatively reproduces experimentally observed effects of the off-state-stabilizing drug mavacamten. Taken together, these results support off-state dynamics as a plausible primary mechanism underlying LDA.

Right Ventricular Sarcomere Contractile Depression and the Role of Thick Filament Activation in Human Heart Failure With Pulmonary Hypertension.

BACKGROUND: Right ventricular (RV) contractile dysfunction commonly occurs and worsens outcomes in patients with heart failure with reduced ejection fraction and pulmonary hypertension (HFrEF-PH). However, such dysfunction often goes undetected by standard clinical RV indices, raising concerns that they may not reflect aspects of underlying myocyte dysfunction. We thus sought to characterize RV myocyte contractile depression in HFrEF-PH, identify those components reflected by clinical RV indices, and uncover underlying biophysical mechanisms. METHODS: Resting, calcium-, and load-dependent mechanics were prospectively studied in permeabilized RV cardiomyocytes isolated from explanted hearts from 23 patients with HFrEF-PH undergoing cardiac transplantation and 9 organ donor controls. RESULTS: Unsupervised machine learning using myocyte mechanical data with the highest variance yielded 2 HFrEF-PH subgroups that in turn mapped to patients with decompensated or compensated clinical RV function. This correspondence was driven by reduced calcium-activated isometric tension in decompensated clinical RV function, whereas surprisingly, many other major myocyte contractile measures including peak power and myocyte active stiffness were similarly depressed in both groups. Similar results were obtained when subgroups were first defined by clinical indices, and then myocyte mechanical properties in each group compared. To test the role of thick filament defects, myofibrillar structure was assessed by x-ray diffraction of muscle fibers. This revealed more myosin heads associated with the thick filament backbone in decompensated clinical RV function, but not compensated clinical RV function, as compared with controls. This corresponded to reduced myosin ATP turnover in decompensated clinical RV function myocytes, indicating less myosin in a crossbridge-ready disordered-relaxed (DRX) state. Altering DRX proportion (%DRX) affected peak calcium-activated tension in the patient groups differently, depending on their basal %DRX, highlighting potential roles for precision-guided therapeutics. Last, increasing myocyte preload (sarcomere length) increased %DRX 1.5-fold in controls but only 1.2-fold in both HFrEF-PH groups, revealing a novel mechanism for reduced myocyte active stiffness and by extension Frank-Starling reserve in human heart failure. CONCLUSIONS: Although there are many RV myocyte contractile deficits in HFrEF-PH, commonly used clinical indices only detect reduced isometric calcium-stimulated force, which is related to deficits in basal and recruitable %DRX myosin. Our results support use of therapies to increase %DRX and enhance length-dependent recruitment of DRX myosin heads in such patients.

Hypertrophic cardiomyopathy ß-cardiac myosin mutation (P710R) leads to hypercontractility by disrupting super relaxed state.

Hypertrophic cardiomyopathy (HCM) is the most common inherited form of heart disease, associated with over 1,000 mutations, many in ß-cardiac myosin (MYH7). Molecular studies of myosin with different HCM mutations have revealed a diversity of effects on ATPase and load-sensitive rate of detachment from actin. It has been difficult to predict how such diverse molecular effects combine to influence forces at the cellular level and further influence cellular phenotypes. This study focused on the P710R mutation that dramatically decreased in vitro motility velocity and actin-activated ATPase, in contrast to other MYH7 mutations. Optical trap measurements of single myosin molecules revealed that this mutation reduced the step size of the myosin motor and the load sensitivity of the actin detachment rate. Conversely, this mutation destabilized the super relaxed state in longer, two-headed myosin constructs, freeing more heads to generate force. Micropatterned human induced pluripotent derived stem cell (hiPSC)-cardiomyocytes CRISPR-edited with the P710R mutation produced significantly increased force (measured by traction force microscopy) compared with isogenic control cells. The P710R mutation also caused cardiomyocyte hypertrophy and cytoskeletal remodeling as measured by immunostaining and electron microscopy. Cellular hypertrophy was prevented in the P710R cells by inhibition of ERK or Akt. Finally, we used a computational model that integrated the measured molecular changes to predict the measured traction forces. These results confirm a key role for regulation of the super relaxed state in driving hypercontractility in HCM with the P710R mutation and demonstrate the value of a multiscale approach in revealing key mechanisms of disease.

Predominant myosin superrelaxed state in canine myocardium with naturally occurring dilated cardiomyopathy.

Dilated cardiomyopathy (DCM) is a naturally occurring heart failure condition in humans and dogs, notably characterized by a reduced contractility and ejection fraction. As the identification of its underlying cellular and molecular mechanisms remain incomplete, the aim of the present study was to assess whether the molecular motor myosin and its known relaxed conformational states are altered in DCM. For that, we dissected and skinned thin cardiac strips from left ventricle obtained from six DCM Doberman Pinschers and six nonfailing (NF) controls. We then used a combination of Mant-ATP chase experiments and X-ray diffraction to assess both energetic and structural changes of myosin. Using the Mant-ATP chase protocol, we observed that in DCM dogs, the amount of myosin molecules in the ATP-conserving conformational state, also known as superrelaxed (SRX), is significantly increased when compared with NF dogs. This alteration can be rescued by applying EMD-57033, a small molecule activating myosin. Conversely, with X-ray diffraction, we found that in DCM dogs, there is a higher proportion of myosin heads in the vicinity of actin when compared with NF dogs (1,0 to 1,1 intensity ratio). Hence, we observed an uncoupling between energetic (Mant-ATP chase) and structural (X-ray diffraction) data. Taken together, these results may indicate that in the heart of Doberman Pinschers with DCM, myosin molecules are potentially stuck in a nonsequestered but ATP-conserving SRX state, that can be counterbalanced by EMD-57033 demonstrating the potential for a myosin-centered pharmacological treatment of DCM.NEW & NOTEWORTHY The key finding of the present study is that, in left ventricles of dogs with a naturally occurring dilated cardiomyopathy, relaxed myosin molecules favor a nonsequestered superrelaxed state potentially impairing sarcomeric contractility. This alteration is rescuable by applying a small molecule activating myosin known as EMD-57033.

History-dependent muscle resistance to stretch remains high after small, posturally relevant pre-movements.

The contributions of intrinsic muscle fiber resistance during mechanical perturbations to standing and other postural behaviors are unclear. Muscle short-range stiffness is known to vary depending on the current level and history of the muscle's activation, as well as the muscle's recent movement history; this property has been referred to as history dependence or muscle thixotropy. However, we currently lack sufficient data about the degree to which muscle stiffness is modulated across posturally relevant characteristics of muscle stretch and activation. We characterized the history dependence of muscle's resistance to stretch in single, permeabilized, activated, muscle fibers in posturally relevant stretch conditions and activation levels. We used a classic paired muscle stretch paradigm, varying the amplitude of a 'conditioning' triangular stretch-shorten cycle followed by a 'test' ramp-and-hold imposed after a variable inter-stretch interval. We tested low (<15%), intermediate (15-50%) and high (>50%) muscle fiber activation levels, evaluating short-range stiffness and total impulse in the test stretch. Muscle fiber resistance to stretch remained high at conditioning amplitudes of <1% optimal fiber length, L0, and inter-stretch intervals of >1 s, characteristic of healthy standing postural sway. An ∼70% attenuation of muscle resistance to stretch was reached at conditioning amplitudes of >3% L0 and inter-stretch intervals of <0.1 s, characteristic of larger, faster postural sway in balance-impaired individuals. The thixotropic changes cannot be predicted solely on muscle force at the time of stretch. Consistent with the disruption of muscle cross-bridges, muscle resistance to stretch during behavior can be substantially attenuated if the prior motion is large enough and/or frequent enough.

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

FiberSim: A flexible open-source model of myofilament-level contraction.

FiberSim is a flexible open-source model of myofilament-level contraction. The code uses a spatially explicit technique, meaning that it tracks the position and status of each contractile molecule within the lattice framework. This allows the model to simulate some of the mechanical effects modulated by myosin-binding protein C, as well as the dose dependence of myotropes and the effects of varying isoform expression levels. This paper provides a short introduction to FiberSim and presents simulations of tension-pCa curves with and without regulation of thick and thin filament activation by myosin-binding protein C. A myotrope dose-dependent response as well as slack/re-stretch maneuvers to assess rates of tension recovery are also presented. The software was designed to be flexible (the user can define their own model and/or protocol) and computationally efficient (simulations can be performed on a regular laptop). We hope that other investigators will use FiberSim to explore myofilament level mechanisms and to accelerate research focusing on the contractile properties of sarcomeres.

Reduced cardiac muscle power with low ATP simulating heart failure.

For patients with heart failure, myocardial ATP level can be reduced to one-half of that observed in healthy controls. This marked reduction (from ≈8 mM in healthy controls to as low as 3-4 mM in heart failure) has been suggested to contribute to impaired myocardial contraction and to the decreased pump function characteristic of heart failure. However, in vitro measures of maximum myofilament force generation, maximum shortening velocity, and the actomyosin ATPase activity show effective KM values for MgATP ranging from ≈10 µM to 150 µM, well below the intracellular ATP level in heart failure. Thus, it is not clear that the fall of myocardial ATP observed in heart failure is sufficient to impair the function of the contractile proteins. Therefore, we tested the effect of low MgATP levels on myocardial contraction using demembranated cardiac muscle preparations that were exposed to MgATP levels typical of the range found in non-failing and failing hearts. Consistent with previous studies, we found that a 50% reduction in MgATP level (from 8 mM to 4 mM) did not reduce maximum force generation or maximum velocity of shortening. However, we found that a 50% reduction in MgATP level caused a 20%-25% reduction in maximal power generation (measured during muscle shortening against a load) and a 20% slowing of cross-bridge cycling kinetics. These results suggest that the decreased cellular ATP level occurring in heart failure contributes to the impaired pump function of the failing heart. Since the ATP-myosin ATPase dissociation constant is estimated to be submillimolar, these findings also suggest that MgATP concentration affects cross-bridge dynamics through a mechanism that is more complex than through the direct dependence of MgATP concentration on myosin ATPase activity. Finally, these studies suggest that therapies targeted to increase adenine nucleotide pool levels in cardiomyocytes might be beneficial for treating heart failure.

Chaperone-mediated autophagy protects cardiomyocytes against hypoxic-cell death.

Chaperone-mediated autophagy (CMA) is a chaperone-dependent process of selective cytosolic protein turnover that targets specific proteins to lysosomes for degradation. Enhancing protein degradation mechanisms has been shown to be beneficial in multiple models of cardiac disease, including myocardial infarction (MI) and ischemia-reperfusion (I/R) injury. However, the causal role of CMA in cardiomyocyte injury and death is largely unknown. Hypoxia is an important contributor to both MI and I/R damage, which are major, precedent causes of heart failure. Upregulating CMA was hypothesized to protect against hypoxia-induced cardiomyocyte death. Lysosome-associated membrane protein 2a (Lamp2a) overexpression and knockdown were used to causally study CMA's role in hypoxically stressed cardiomyocytes. LAMP2a protein levels were used as both a primary indicator and driver of CMA function. Hypoxic stress was stimulated by CoCl2 treatment, which increased LAMP2a protein levels (+1.4-fold) and induced cardiomyocyte apoptosis (+3.2-4.0-fold). Lamp2a siRNA knockdown (-3.2-fold) of control cardiomyocytes increased apoptosis (+1.8-fold) suggesting that loss of CMA is detrimental for cardiomyocyte survival. However, there was neither an additive nor a synergistic effect on cell death when Lamp2a-silenced cells were treated with CoCl2. Conversely, Lamp2a overexpression (+3.0-fold) successfully reduced hypoxia-induced apoptosis by ∼50%. LAMP2a was also significantly increased (+1.7-fold) in ischemic heart failure patient samples, similar to hypoxically stressed cardiomyocytes. The failing ischemic hearts may have had insufficient CMA activation. To our knowledge, this study for the first time establishes a protective role for CMA (via Lamp2a overexpression) against hypoxia-induced cardiomyocyte loss and reveals the intriguing possibility that CMA activation may offer a cardioprotective treatment for ischemic heart disease.

Renal Angiotensinogen Is Predominantly Liver Derived in Nonhuman Primates.

[Figure: see text].

MgADP Promotes Myosin Head Movement toward Actin at Low [Ca2+] to Increase Force Production and Ca2+-Sensitivity of Contraction in Permeabilized Porcine Myocardial Strips.

Myosin cross-bridges dissociate from actin following Mg2+-adenosine triphosphate (MgATP) binding. Myosin hydrolyses MgATP into inorganic phosphate (Pi) and Mg2+-adenosine diphosphate (ADP), and release of these hydrolysis products drives chemo-mechanical energy transitions within the cross-bridge cycle to power muscle contraction. Some forms of heart disease are associated with metabolic or enzymatic dysregulation of the MgATP-MgADP nucleotide pool, resulting in elevated cytosolic [MgADP] and impaired muscle relaxation. We investigated the mechanical and structural effects of increasing [MgADP] in permeabilized myocardial strips from porcine left ventricle samples. Sarcomere length was set to 2.0 µm at 28 °C, and all solutions contained 3% dextran T-500 to compress myofilament lattice spacing to near-physiological values. Under relaxing low [Ca2+] conditions (pCa 8.0, where pCa = -log10[Ca2+]), tension increased as [MgADP] increased from 0-5 mM. Complementary small-angle X-ray diffraction measurements show that the equatorial intensity ratio, I1,1/I1,0, also increased as [MgADP] increased from 0 to 5 mM, indicating myosin head movement away from the thick-filament backbone towards the thin-filament. Ca2+-activated force-pCa measurements show that Ca2+-sensitivity of contraction increased with 5 mM MgADP, compared to 0 mM MgADP. These data show that MgADP augments tension at low [Ca2+] and Ca2+-sensitivity of contraction, suggesting that MgADP destabilizes the quasi-helically ordered myosin OFF state, thereby shifting the cross-bridge population towards the disordered myosin ON state. Together, these results indicate that MgADP enhances the probability of cross-bridge binding to actin due to enhancement of both thick and thin filament-based activation mechanisms.

Cardiac myosin regulatory light chain kinase modulates cardiac contractility by phosphorylating both myosin regulatory light chain and troponin I.

Heart muscle contractility and performance are controlled by posttranslational modifications of sarcomeric proteins. Although myosin regulatory light chain (RLC) phosphorylation has been studied extensively in vitro and in vivo, the precise role of cardiac myosin light chain kinase (cMLCK), the primary kinase acting upon RLC, in the regulation of cardiomyocyte contractility remains poorly understood. In this study, using recombinantly expressed and purified proteins, various analytical methods, in vitro and in situ kinase assays, and mechanical measurements in isolated ventricular trabeculae, we demonstrate that human cMLCK is not a dedicated kinase for RLC but can phosphorylate other sarcomeric proteins with well-characterized regulatory functions. We show that cMLCK specifically monophosphorylates Ser23 of human cardiac troponin I (cTnI) in isolation and in the trimeric troponin complex in vitro and in situ in the native environment of the muscle myofilament lattice. Moreover, we observed that human cMLCK phosphorylates rodent cTnI to a much smaller extent in vitro and in situ, suggesting species-specific adaptation of cMLCK. Although cMLCK treatment of ventricular trabeculae exchanged with rat or human troponin increased their cross-bridge kinetics, the increase in sensitivity of myofilaments to calcium was significantly blunted by human TnI, suggesting that human cTnI phosphorylation by cMLCK modifies the functional consequences of RLC phosphorylation. We propose that cMLCK-mediated phosphorylation of TnI is functionally significant and represents a critical signaling pathway that coordinates the regulatory states of thick and thin filaments in both physiological and potentially pathophysiological conditions of the heart.

Fast-relaxing cardiomyocytes exert a dominant role in the relaxation behavior of heterogeneous myocardium.

Substantial variation in relaxation rate exists among cardiomyocytes within small volumes of myocardium; however, it is unknown how this variability affects the overall relaxation mechanics of heart muscle. In this study, we sought to modulate levels of cellular heterogeneity in a computational model, then validate those predictions using an engineered heart tissue platform. We formulated an in silico tissue model composed of half-sarcomeres with varied relaxation rates, incorporating single-cell cardiomyocyte experimental data. These model tissues randomly sampled relaxation parameters from two offset distributions of fast- and slow-relaxing populations of half-sarcomeres. Isometric muscle twitch simulations predicted a complex relationship between relaxation time and the proportion of fast-versus slow-relaxing cells in heterogeneous tissues. Specifically, a 50/50 mixture of fast and slow cells did not lead to relaxation time that was the mean of the relaxation times associated with the two pure cases. Rather, the mean relaxation time was achieved at a ratio of 70:30 slow:fast relaxing cells, suggesting a disproportionate impact of fast-relaxing cells on overall tissue relaxation. To examine whether this behavior persists in vitro, we constructed engineered heart tissues from two lines of fast- and slow-relaxing human iPSC-derived cardiomyocytes. Cell tracking via fluorescent nanocrystals confirmed the presence of both cell populations in the 50/50 mixed tissues at the time of mechanical characterization. Isometric muscle twitch relaxation times of these mixed-population engineered heart tissues showed agreement with the predictions from the model, namely that the measured relaxation rate of 50/50 mixed tissues more closely resembled that of tissues made with 100% fast-relaxing cells. Our observations suggest that cardiomyocyte diversity can play an important role in determining tissue-level relaxation.

Diabetic microcirculatory disturbances and pathologic erythropoiesis are provoked by deposition of amyloid-forming amylin in red blood cells and capillaries.

In the setting of type-2 diabetes, there are declines of structural stability and functionality of blood capillaries and red blood cells (RBCs), increasing the risk for microcirculatory disturbances. Correcting hyperglycemia is not entirely effective at reestablishing normal cellular metabolism and function. Therefore, identification of pathological changes occurring before the development of overt hyperglycemia may lead to novel therapeutic targets for reducing the risk of microvascular dysfunction. Here we determine whether RBC-capillary interactions are altered by prediabetic hypersecretion of amylin, an amyloid forming hormone co-synthesized with insulin, and is reversed by endothelial cell-secreted epoxyeicosatrienoic acids. In patients, we found amylin deposition in RBCs in association with type-2 diabetes, heart failure, cancer and stroke. Amylin-coated RBCs have altered shape and reduced functional (non-glycated) hemoglobin. Amylin-coated RBCs administered intravenously in control rats upregulated erythropoietin and renal arginase expression and activity. We also found that diabetic rats expressing amyloid-forming human amylin in the pancreas (the HIP rat model) have increased tissue levels of hypoxia-inducible transcription factors, compared to diabetic rats that express non-amyloid forming rat amylin (the UCD rat model). Upregulation of erythropoietin correlated with lower hematocrit in the HIP model indicating pathologic erythropoiesis. In the HIP model, pharmacological upregulation of endogenous epoxyeicosatrienoic acids protected the renal microvasculature against amylin deposition and also reduced renal accumulation of HIFs. Thus, prediabetes induces dysregulation of amylin homeostasis and promotes amylin deposition in RBCs and the microvasculature altering RBC-capillary interaction leading to activation of hypoxia signaling pathways and pathologic erythropoiesis. Hence, dysregulation of amylin homeostasis could be a therapeutic target for ameliorating diabetic vascular complications.

A short history of the development of mathematical models of cardiac mechanics.

Cardiac mechanics plays a crucial role in atrial and ventricular function, in the regulation of growth and remodelling, in the progression of disease, and the response to treatment. The spatial scale of the critical mechanisms ranges from nm (molecules) to cm (hearts) with the fastest events occurring in milliseconds (molecular events) and the slowest requiring months (growth and remodelling). Due to its complexity and importance, cardiac mechanics has been studied extensively both experimentally and through mathematical models and simulation. Models of cardiac mechanics evolved from seminal studies in skeletal muscle, and developed into cardiac specific, species specific, human specific and finally patient specific calculations. These models provide a formal framework to link multiple experimental assays recorded over nearly 100 years into a single unified representation of cardiac function. This review first provides a summary of the proteins, physiology and anatomy involved in the generation of cardiac pump function. We then describe the evolution of models of cardiac mechanics starting with the early theoretical frameworks describing the link between sarcomeres and muscle contraction, transitioning through myosin-level models to calcium-driven systems, and ending with whole heart patient-specific models.

Differential Effects of Isoproterenol on Regional Myocardial Mechanics in Rat Using Three-Dimensional 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, midventricular, and apical regions. Additionally, the midventricular 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 (CL) shear strain, as well as ventricular torsion, throughout the basal, midventricle, and apical regions. In the midventricle, 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 CL shear strain in all midventricular segments. Ventricular torsion was significantly increased in nearly all segments except the inferior subendocardium. The effects of isoproterenol on LV systolic mechanics (i.e., three-dimensional (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.

Force-Dependent Recruitment from the Myosin Off State Contributes to Length-Dependent Activation.

Cardiac muscle develops more force when it is activated at longer lengths. The concentration of Ca2+ required to develop half-maximal force also decreases. These effects are known as length-dependent activation and are thought to play critical roles in the Frank-Starling relationship and cardiovascular homeostasis. The molecular mechanisms underpinning length-dependent activation remain unclear, but recent experiments suggest that they may include recruitment of myosin heads from the off (sometimes called super-relaxed) state. This manuscript presents a mathematical model of muscle contraction that was developed to investigate this hypothesis. Myosin heads in the model transitioned between an off state (that could not interact with actin), an on state (that could bind to actin), and a single attached state. Simulations were fitted to experimental data using multidimensional parameter optimization. Statistical analysis showed that a model in which the rate of the off-to-on transition increased linearly with force reproduced the length-dependent behavior of chemically permeabilized myocardium better than a model with a constant off-to-on transition rate (F-test, p < 0.001). This result suggests that the thick-filament transitions are modulated by force. Additional calculations showed that the model incorporating a mechanosensitive thick filament could also reproduce twitch responses measured in a trabecula stretched to different lengths. A final set of simulations was then used to test the model. These calculations predicted how reducing passive stiffness would impact the length dependence of the calcium sensitivity of contractile force. The prediction (a 60% reduction in ΔpCa50) mimicked the 58% reduction in ΔpCa50 in myocardium from rats that expressed a giant isoform of titin and had low resting tension. Together, these computational results suggest that force-dependent recruitment of myosin heads from the thick-filament off state contributes to length-dependent activation and the Frank-Starling relationship.