Myosin in autoinhibited off state(s), stabilized by mavacamten, can be recruited in response to inotropic interventions.

Mavacamten is a FDA-approved small-molecule therapeutic designed to regulate cardiac function at the sarcomere level by selectively but reversibly inhibiting the enzymatic activity of myosin. It shifts myosin toward ordered off states close to the thick filament backbone. It remains elusive whether these myosin heads in the off state(s) can be recruited in response to physiological stimuli when required to boost cardiac output. We show that cardiac myosins stabilized in these off state(s) by mavacamten are recruitable by 1) Ca2+, 2) increased chronotropy [heart rate (HR)], 3) stretch, and 4) ß-adrenergic (ß-AR) stimulation, all known physiological inotropic interventions. At the molecular level, we show that Ca2+ increases myosin ATPase activity by shifting mavacamten-stabilized myosin heads from the inactive super-relaxed state to the active disordered relaxed state. At the myofilament level, both Ca2+ and passive lengthening can shift mavacamten-ordered off myosin heads from positions close to the thick filament backbone to disordered on states closer to the thin filaments. In isolated rat cardiomyocytes, increased stimulation rates enhanced shortening fraction in mavacamten-treated cells. This observation was confirmed in vivo in telemetered rats, where left-ventricular dP/dtmax, an index of inotropy, increased with HR in mavacamten-treated animals. Finally, we show that ß-AR stimulation in vivo increases left-ventricular function and stroke volume in the setting of mavacamten. Our data demonstrate that the mavacamten-promoted off states of myosin in the thick filament are at least partially activable, thus preserving cardiac reserve mechanisms.

Using Multiscale Simulations as a Tool to Interpret Equatorial X-ray Fiber Diffraction Patterns from Skeletal Muscle.

Synchrotron small-angle X-ray diffraction is the method of choice for nm-scale structural studies of striated muscle under physiological conditions and on millisecond time scales. The lack of generally applicable computational tools for modeling X-ray diffraction patterns from intact muscles has been a significant barrier to exploiting the full potential of this technique. Here, we report a novel "forward problem" approach using the spatially explicit computational simulation platform MUSICO to predict equatorial small-angle X-ray diffraction patterns and the force output simultaneously from resting and isometrically contracting rat skeletal muscle that can be compared to experimental data. The simulation generates families of thick-thin filament repeating units, each with their individually predicted occupancies of different populations of active and inactive myosin heads that can be used to generate 2D-projected electron density models based on known Protein Data Bank structures. We show how, by adjusting only a few selected parameters, we can achieve a good correspondence between experimental and predicted X-ray intensities. The developments presented here demonstrate the feasibility of combining X-ray diffraction and spatially explicit modeling to form a powerful hypothesis-generating tool that can be used to motivate experiments that can reveal emergent properties of muscle.

Effect of Myosin Isoforms on Cardiac Muscle Twitch of Mice, Rats and Humans.

To understand how pathology-induced changes in contractile protein isoforms modulate cardiac muscle function, it is necessary to quantify the temporal-mechanical properties of contractions that occur under various conditions. Pathological responses are much easier to study in animal model systems than in humans, but extrapolation between species presents numerous challenges. Employing computational approaches can help elucidate relationships that are difficult to test experimentally by translating the observations from rats and mice, as model organisms, to the human heart. Here, we use the spatially explicit MUSICO platform to model twitch contractions from rodent and human trabeculae collected in a single laboratory. This approach allowed us to identify the variations in kinetic characteristics of α- and ß-myosin isoforms across species and to quantify their effect on cardiac muscle contractile responses. The simulations showed how the twitch transient varied with the ratio of the two myosin isoforms. Particularly, the rate of tension rise was proportional to the fraction of α-myosin present, while the ß-isoform dominated the rate of relaxation unless α-myosin was >50%. Moreover, both the myosin isoform and the Ca2+ transient contributed to the twitch tension transient, allowing two levels of regulation of twitch contraction.

The effect of variable troponin C mutation thin filament incorporation on cardiac muscle twitch contractions.

One of the complexities of understanding the pathology of familial forms of cardiac diseases is the level of mutation incorporation in sarcomeres. Computational models of the sarcomere that are spatially explicit offer an approach to study aspects of mutational incorporation into myofilaments that are more challenging to get at experimentally. We studied two well characterized mutations of cardiac TnC, L48Q and I61Q, that decrease or increase the release rate of Ca2+ from cTnC, k-Ca, resulting in HCM and DCM respectively [1]. Expression of these mutations in transgenic mice was used to provide experimental data for incorporation of 30 and 50% (respectively) into sarcomeres. Here we demonstrate that fixed length twitch contractions of trabeculae from mice containing mutant differ from WT; L48Q trabeculae have slower relaxation while I61Q trabeculae have markedly reduced peak tension. Using our multiscale modelling approach [2] we were able to describe the tension transients of WT mouse myocardium. Tension transients for the mutant cTnCs were simulated with changes in k-Ca, measured experimentally for each cTnC mutant in whole troponin complex, a change in the affinity of cTnC for cTnI, and a reduction in the number of detached crossbridges available for binding. A major advantage of the multiscale explicit 3-D model is that it predicts the effects of variable mutation incorporation, and the effects of variations in mutation distribution within thin filaments in sarcomeres. Such effects are currently impossible to explore experimentally. We explored random and clustered distributions of mutant cTnCs in thin filaments, as well as distributions of individual thin filaments with only WT or mutant cTnCs present. The effects of variable amounts of incorporation and non-random distribution of mutant cTnCs are more marked for I61Q than L48Q cTnC. We conclude that this approach can be effective for study on mutations in multiple proteins of the sarcomere. SUMMARY: A challenge in experimental studies of diseases is accounting for the effect of variable mutation incorporation into myofilaments. Here we use a spatially explicit computational approach, informed by experimental data from transgenic mice expressing one of two mutations in cardiac Troponin C that increase or decrease calcium sensitivity. We demonstrate that the model can accurately describe twitch contractions for the data and go on to explore the effect of variable mutant incorporation and localization on simulated cardiac muscle twitches.

Estimation of Forces on Actin Filaments in Living Muscle from X-ray Diffraction Patterns and Mechanical Data.

Many biological processes are triggered or driven by mechanical forces in the cytoskeletal network, but these transducing forces have rarely been assessed. Striated muscle, with its well-organized structure provides an opportunity to assess intracellular forces using small-angle X-ray fiber diffraction. We present a new methodology using Monte Carlo simulations of muscle contraction in an explicit 3D sarcomere lattice to predict the fiber deformations and length changes along thin filaments during contraction. Comparison of predicted diffraction patterns to experimental meridional X-ray reflection profiles allows assessment of the stepwise changes in intermonomer spacings and forces in the myofilaments within living muscle cells. These changes along the filament length reflect the effect of forces from randomly attached crossbridges. This approach enables correlation of the molecular events, such as the current number of attached crossbridges and the distributions of crossbridge forces to macroscopic measurements of force and length changes during muscle contraction. In addition, assessments of fluctuations in local forces in the myofilaments may reveal how variations in the filament forces acting on signaling proteins in the sarcomere M-bands and Z-discs modulate gene expression, protein synthesis and degradation, and as well to mechanisms of adaptation of muscle in response to changes in mechanical loading.

Computational Modeling on Drugs Effects for Left Ventricle in Cardiomyopathy Disease.

Cardiomyopathy is associated with structural and functional abnormalities of the ventricular myocardium and can be classified in two major groups: hypertrophic (HCM) and dilated (DCM) cardiomyopathy. Computational modeling and drug design approaches can speed up the drug discovery and significantly reduce expenses aiming to improve the treatment of cardiomyopathy. In the SILICOFCM project, a multiscale platform is developed using coupled macro- and microsimulation through finite element (FE) modeling of fluid-structure interactions (FSI) and molecular drug interactions with the cardiac cells. FSI was used for modeling the left ventricle (LV) with a nonlinear material model of the heart wall. Simulations of the drugs' influence on the electro-mechanics LV coupling were separated in two scenarios, defined by the principal action of specific drugs. We examined the effects of Disopyramide and Dygoxin which modulate Ca2+ transients (first scenario), and Mavacamten and 2-deoxy adenosine triphosphate (dATP) which affect changes of kinetic parameters (second scenario). Changes of pressures, displacements, and velocity distributions, as well as pressure-volume (P-V) loops in the LV models of HCM and DCM patients were presented. Additionally, the results obtained from the SILICOFCM Risk Stratification Tool and PAK software for high-risk HCM patients closely followed the clinical observations. This approach can give much more information on risk prediction of cardiac disease to specific patients and better insight into estimated effects of drug therapy, leading to improved patient monitoring and treatment.

Myosin in autoinhibited off state(s), stabilized by mavacamten, can be recruited via inotropic effectors.

Mavacamten is a novel, FDA-approved, small molecule therapeutic designed to regulate cardiac function by selectively but reversibly inhibiting the enzymatic activity of myosin. It shifts myosin towards ordered off states close to the thick filament backbone. It remains unresolved whether mavacamten permanently sequesters these myosin heads in the off state(s) or whether these heads can be recruited in response to physiological stimuli when required to boost cardiac output. We show that cardiac myosins stabilized in these off state(s) by mavacamten are recruitable by Ca2+, increased heart rate, stretch, and ß-adrenergic (ß-AR) stimulation, all known physiological inotropic effectors. At the molecular level, we show that, in presence of mavacamten, Ca2+ increases myosin ATPase activity by shifting myosin heads from the reserve super-relaxed (SRX) state to the active disordered relaxed (DRX) state. At the myofilament level, both Ca2+ and passive lengthening can shift ordered off myosin heads from positions close to the thick filament backbone to disordered on states closer to the thin filaments in the presence of mavacamten. In isolated rat cardiomyocytes, increased stimulation rates enhanced shortening fraction in mavacamten-treated cells. This observation was confirmed in vivo in telemetered rats, where left-ventricular dP/dtmax, an index of inotropy, increased with heart rate in mavacamten treated animals. Finally, we show that ß-AR stimulation in vivo increases left-ventricular function and stroke volume in the setting of mavacamten. Our data demonstrate that the mavacamten-promoted off states of myosin in the thick filament are activable, at least partially, thus leading to preservation of cardiac reserve mechanisms.

SILICOFCM platform, multiscale modeling of left ventricle from echocardiographic images and drug influence for cardiomyopathy disease.

BACKGROUND AND OBJECTIVE: In silico clinical trials are the future of medicine and virtual testing and simulation are the future of medical engineering. The use of a computational platform can reduce costs and time required for developing new models of medical devices and drugs. The computational platform, which is one of the main results of the SILICOFCM project, was developed using state-of-the-art finite element modeling for macro simulation of fluid-structure interaction with micro modeling at the molecular level for drug interaction with the cardiac cells. SILICOFCM platform is using for risk prediction and optimal drug therapy of familial cardiomyopathy in a specific patient. METHODS: In order to obtain 3D image reconstruction, the U-net architecture was used to determine geometric parameters for the left ventricle which were extracted from the echocardiographic apical and M-mode views. A micro-mechanics cellular model which includes three kinetic processes of sarcomeric proteins interactions was developed. It allows simulation of the drugs which are divided into three major groups defined by the principal action of each drug. Fluid-solid coupling for the left ventricle was presented. A nonlinear material model of the heart wall that was developed by using constitutive curves which include the stress-strain relationship was used. RESULTS: The results obtained with the parametric model of the left ventricle where pressure-volume (PV) diagrams depend on the change of Ca2+ were presented. It directly affects the ejection fraction. The presented approach with the variation of the left ventricle (LV) geometry and simulations which include the influence of different parameters on the PV diagrams are directly interlinked with drug effects on the heart function. It includes different drugs such as Entresto and Digoxin that directly affect the cardiac PV diagrams and ejection fraction. CONCLUSIONS: Computational platforms such as the SILICOFCM platform are novel tools for risk prediction of cardiac disease in a specific patient that will certainly open a new avenue for in silico clinical trials in the future.

Multiscale modeling of twitch contractions in cardiac trabeculae.

Understanding the dynamics of a cardiac muscle twitch contraction is complex because it requires a detailed understanding of the kinetic processes of the Ca2+ transient, thin-filament activation, and the myosin-actin cross-bridge chemomechanical cycle. Each of these steps has been well defined individually, but understanding how all three of the processes operate in combination is a far more complex problem. Computational modeling has the potential to provide detailed insight into each of these processes, how the dynamics of each process affect the complexity of contractile behavior, and how perturbations such as mutations in sarcomere proteins affect the complex interactions of all of these processes. The mechanisms involved in relaxation of tension during a cardiac twitch have been particularly difficult to discern due to nonhomogeneous sarcomere lengthening during relaxation. Here we use the multiscale MUSICO platform to model rat trabecular twitches. Validation of computational models is dependent on being able to simulate different experimental datasets, but there has been a paucity of data that can provide all of the required parameters in a single experiment, such as simultaneous measurements of force, intracellular Ca2+ transients, and sarcomere length dynamics. In this study, we used data from different studies collected under similar experimental conditions to provide information for all the required parameters. Our simulations established that twitches either in an isometric sarcomere or in fixed-length, multiple-sarcomere trabeculae replicate the experimental observations if models incorporate a length-tension relationship for the nonlinear series elasticity of muscle preparations and a scheme for thick-filament regulation. The thick-filament regulation assumes an off state in which myosin heads are parked onto the thick-filament backbone and are unable to interact with actin, a state analogous to the super-relaxed state. Including these two mechanisms provided simulations that accurately predict twitch contractions over a range of different conditions.

X-ray diffraction from nonuniformly stretched helical molecules.

The fibrous proteins in living cells are exposed to mechanical forces interacting with other subcellular structures. X-ray fiber diffraction is often used to assess deformation and movement of these proteins, but the analysis has been limited to the theory for fibrous molecular systems that exhibit helical symmetry. However, this approach cannot adequately interpret X-ray data from fibrous protein assemblies where the local strain varies along the fiber length owing to interactions of its molecular constituents with their binding partners. To resolve this problem a theoretical formulism has been developed for predicting the diffraction from individual helical molecular structures nonuniformly strained along their lengths. This represents a critical first step towards modeling complex dynamical systems consisting of multiple helical structures using spatially explicit, multi-scale Monte Carlo simulations where predictions are compared with experimental data in a 'forward' process to iteratively generate ever more realistic models. Here the effects of nonuniform strains and the helix length on the resulting magnitude and phase of diffraction patterns are quantitatively assessed. Examples of the predicted diffraction patterns of nonuniformly deformed double-stranded DNA and actin filaments in contracting muscle are presented to demonstrate the feasibly of this theoretical approach.