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

Mechanisms of myosin II force generation: insights from novel experimental techniques and approaches

Myosin II is a molecular motor that converts chemical energy derived from ATP hydrolysis into mechanical work. Myosin II isoforms are responsible for muscle contraction and a range of cell functions relying on the development of force and motion. When the motor attaches to actin, ATP is hydrolyzed and inorganic phosphate (Pi) and ADP are released from its active site. These reactions are coordinated with changes in the structure of myosin, promoting the so-called "power stroke" that causes the sliding of actin filaments. The general features of the myosin-actin interactions are well accepted, but there are critical issues that remain poorly understood, mostly due to technological limitations. In recent years, there has been a significant advance in structural, biochemical, and mechanical methods that have advanced the field considerably. New modeling approaches have also allowed researchers to understand actomyosin interactions at different levels of analysis. This paper reviews recent studies looking into the interaction between myosin II and actin filaments, which leads to power stroke and force generation. It reviews studies conducted with single myosin molecules, myosins working in filaments, muscle sarcomeres, myofibrils, and fibers. It also reviews the mathematical models that have been used to understand the mechanics of myosin II in approaches focusing on single molecules to ensembles. Finally, it includes brief sections on translational aspects, how changes in the myosin motor by mutations and/or posttranslational modifications may cause detrimental effects in diseases and aging, among other conditions, and how myosin II has become an emerging drug target.

Multiscale fiber remodeling in the infarcted left ventricle using a stress-based reorientation law

The organization of myofibers and extra cellular matrix within the myocardium plays a significant role in defining cardiac function. When pathological events occur, such as myocardial infarction (MI), this organization can become disrupted, leading to degraded pumping performance. The current study proposes a multiscale finite element (FE) framework to determine realistic fiber distributions in the left ventricle (LV). This is achieved by implementing a stress-based fiber reorientation law, which seeks to align the fibers with local traction vectors, such that contractile force and load bearing capabilities are maximized. By utilizing the total stress (passive and active), both myofibers and collagen fibers are reoriented. Simulations are conducted to predict the baseline fiber configuration in a normal LV as well as the adverse fiber reorientation that occurs due to different size MIs. The baseline model successfully captures the transmural variation of helical fiber angles within the LV wall, as well as the transverse fiber angle variation from base to apex. In the models of MI, the patterns of fiber reorientation in the infarct, border zone, and remote regions closely align with previous experimental findings, with a significant increase in fibers oriented in a left-handed helical configuration and increased dispersion in the infarct region. Furthermore, the severity of fiber reorientation and impairment of pumping performance both showed a correlation with the size of the infarct. The proposed multiscale modeling framework allows for the effective prediction of adverse remodeling and offers the potential for assessing the effectiveness of therapeutic interventions in the future. STATEMENT OF SIGNIFICANCE: The organization of muscle and collagen fibers within the heart plays a significant role in defining cardiac function. This organization can become disrupted after a heart attack, leading to degraded pumping performance. In the current study, we implemented a stress-based fiber reorientation law into a computer model of the heart, which seeks to realign the fibers such that contractile force and load bearing capabilities are maximized. The primary goal was to evaluate the effects of different sized heart attacks. We observed substantial fiber remodeling in the heart, which matched experimental observations. The proposed computational framework allows for the effective prediction of adverse remodeling and offers the potential for assessing the effectiveness of therapeutic interventions in the future.

Advancing clinical translation of cardiac biomechanics models: a comprehensive review, applications and future pathways

Cardiac mechanics models are developed to represent a high level of detail, including refined anatomies, accurate cell mechanics models, and platforms to link microscale physiology to whole-organ function. However, cardiac biomechanics models still have limited clinical translation. In this review, we provide a picture of cardiac mechanics models, focusing on their clinical translation. We review the main experimental and clinical data used in cardiac models, as well as the steps followed in the literature to generate anatomical meshes ready for simulations. We describe the main models in active and passive mechanics and the different lumped parameter models to represent the circulatory system. Lastly, we provide a summary of the state-of-the-art in terms of ventricular, atrial, and four-chamber cardiac biomechanics models. We discuss the steps that may facilitate clinical translation of the biomechanics models we describe. A well-established software to simulate cardiac biomechanics is lacking, with all available platforms involving different levels of documentation, learning curves, accessibility, and cost. Furthermore, there is no regulatory framework that clearly outlines the verification and validation requirements a model has to satisfy in order to be reliably used in applications. Finally, better integration with increasingly rich clinical and/or experimental datasets as well as machine learning techniques to reduce computational costs might increase model reliability at feasible resources. Cardiac biomechanics models provide excellent opportunities to be integrated into clinical workflows, but more refinement and careful validation against clinical data are needed to improve their credibility. In addition, in each context of use, model complexity must be balanced with the associated high computational cost of running these models.

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.

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.

Characterization of Exercise-Induced Myocardium Growth Using Finite Element Modeling and Bayesian Optimization

Cardiomyocyte growth can occur in both physiological (exercised-induced) and pathological (e.g., volume overload and pressure overload) conditions leading to left ventricular (LV) hypertrophy. Studies using animal models and histology have demonstrated the growth and remodeling process at the organ level and tissue-cellular level, respectively. However, the driving factors of growth and the mechanistic link between organ, tissue, and cellular growth remains poorly understood. Computational models have the potential to bridge this gap by using constitutive models that describe the growth and remodeling process of the myocardium coupled with finite element (FE) analysis to model the biomechanics of the heart at the organ level. Using subject-specific imaging data of the LV geometry at two different time points, an FE model can be created with the inverse method to characterize the growth parameters of each subject. In this study, we developed a framework that takes in vivo cardiac magnetic resonance (CMR) imaging data of exercised porcine model and uses FE and Bayesian optimization to characterize myocardium growth in the transverse and longitudinal directions. The efficacy of this framework was demonstrated by successfully predicting growth parameters of 18 synthetic LV targeted masks which were generated from three LV porcine geometries. The framework was further used to characterize growth parameters in 4 swine subjects that had been exercised. The study suggested that exercise-induced growth in swine is prone to longitudinal cardiomyocyte growth (58.0 ± 19.6% after 6 weeks and 79.3 ± 15.6% after 12 weeks) compared to transverse growth (4.0 ± 8.0% after 6 weeks and 7.8 ± 9.4% after 12 weeks). This framework can be used to characterize myocardial growth in different phenotypes of LV hypertrophy and can be incorporated with other growth constitutive models to study different hypothetical growth mechanisms.

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.

Force-dependent recruitment from myosin OFF-state increases end-systolic pressure-volume relationship in left ventricle

Finite element (FE) modeling is becoming increasingly prevalent in the world of cardiac mechanics; however, many existing FE models are phenomenological and thus do not capture cellular-level mechanics. This work implements a cellular-level contraction scheme into an existing nonlinear FE code to model ventricular contraction. Specifically, this contraction model incorporates three myosin states: OFF-, ON-, and an attached force-generating state. It has been speculated that force-dependent transitions from the OFF- to ON-state may contribute to length-dependent activation at the cellular level. The current work investigates the contribution of force-dependent recruitment out of the OFF-state to ventricular-level function, specifically the Frank-Starling relationship, as seen through the end-systolic pressure-volume relationship (ESPVR). Five FE models were constructed using geometries of rat left ventricles obtained via cardiac magnetic resonance imaging. FE simulations were conducted to optimize parameters for the cellular contraction model such that the differences between FE predicted ventricular pressures for the models and experimentally measured pressures were minimized. The models were further validated by comparing FE predicted end-systolic strain to experimentally measured strain. Simulations mimicking vena cava occlusion generated descending pressure volume loops from which ESPVRs were calculated. In simulations with the inclusion of the OFF-state, using a force-dependent transition to the ON-state, the ESPVR calculated was steeper than in simulations excluding the OFF-state. Furthermore, the ESPVR was also steeper when compared to models that included the OFF-state without a force-dependent transition. This suggests that the force-dependent recruitment of thick filament heads from the OFF-state at the cellular level contributes to the Frank-Starling relationship observed at the organ level.

Multiscale Models of Cardiac Muscle Biophysics and Tissue Remodeling in Hypertrophic Cardiomyopathies

Myocardial hypertrophy is the result of sustained perturbations to the mechanical and/or neurohormonal homeostasis of cardiac cells and is driven by integrated, multiscale biophysical and biochemical processes that are currently not well defined. In this brief review, we highlight recent computational and experimental models of cardiac hypertrophy that span mechanisms from the molecular level to the tissue level. Specifically, we focus on: (i) molecular-level models of the structural dynamics of sarcomere proteins in hypertrophic hearts, (ii) cellular-level models of excitation-contraction coupling and mechanosensitive signaling in disease-state myocytes, and (iii) organ-level models of myocardial growth kinematics and predictors thereof. Finally, we discuss how spanning these scales and combining multiple experimental/computational models will provide new information about the processes governing hypertrophy and potential methods to prevent or reverse them.

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.