Beta-adrenergic stimulation maintains cardiac function in Serca2 knockout mice

Myocardial SERCA2 Protects Against Cardiac Damage and Dysfunction Caused by Inhaled Bromine

Myocardial sarcoendoplasmic reticulum calcium ATPase 2 (SERCA2) activity is critical for heart function. We have demonstrated that inhaled halogen (chlorine or bromine) gases inactivate SERCA2, impair calcium homeostasis, increase proteolysis, and damage the myocardium ultimately leading to cardiac dysfunction. To further elucidate the mechanistic role of SERCA2 in halogen-induced myocardial damage, we used bromine-exposed cardiac-specific SERCA2 knockout (KO) mice [tamoxifen-administered SERCA2 (flox/flox) Tg (αMHC-MerCreMer) mice] and compared them to the oil-administered controls. We performed echocardiography and hemodynamic analysis to investigate cardiac function 24 hours after bromine (600 ppm for 30 minutes) exposure and measured cardiac injury markers in plasma and proteolytic activity in cardiac tissue and performed electron microscopy of the left ventricle (LV). Cardiac-specific SERCA2 knockout mice demonstrated enhanced toxicity to bromine. Bromine exposure increased ultrastructural damage, perturbed LV shape geometry, and demonstrated acutely increased phosphorylation of phospholamban in the KO mice. Bromine-exposed KO mice revealed significantly enhanced mean arterial pressure and sphericity index and decreased LV end diastolic diameter and LV end systolic pressure when compared with the bromine-exposed control FF mice. Strain analysis showed loss of synchronicity, evidenced by an irregular endocardial shape in systole and irregular vector orientation of contractile motion across different segments of the LV in KO mice, both at baseline and after bromine exposure. These studies underscore the critical role of myocardial SERCA2 in preserving cardiac ultrastructure and function during toxic halogen gas exposures. SIGNIFICANCE STATEMENT: Due to their increased industrial production and transportation, halogens such as chlorine and bromine pose an enhanced risk of exposure to the public. Our studies have demonstrated that inhalation of these halogens leads to the inactivation of cardiopulmonary SERCA2 and results in calcium overload. Using cardiac-specific SERCA2 KO mice, these studies further validated the role of SERCA2 in bromine-induced myocardial injury. These studies highlight the increased susceptibility of individuals with pathological loss of cardiac SERCA2 to the effects of bromine.

Myofilament-based physiological regulatory compensation preserves diastolic function in failing hearts with severe Ca2+ handling deficits

Severe dysfunction in cardiac muscle intracellular Ca2+ handling is a common pathway underlying heart failure. Here we used an inducible genetic model of severe Ca2+ cycling dysfunction by the targeted temporal gene ablation of the cardiac Ca2+ ATPase, SERCA2, in otherwise normal adult mice. In this model, in vivo heart performance was minimally affected initially, even though Serca2a protein was markedly reduced. The mechanism underlying the sustained in vivo heart performance in the weeks prior to complete heart pump failure and death is not clear and is important to understand. Studies were primarily focused on understanding how in vivo diastolic function could be relatively normal under conditions of marked Serca2a deficiency. Interestingly, data show increased cardiac troponin I (cTnI) serine 23/24 phosphorylation content in hearts upon Serca2a ablation in vivo. We report that hearts isolated from the Serca2-deficient mice retained near normal heart pump functional responses to β-adrenergic stimulation. Unexpectedly, using genetic complementation models, in concert with inducible Serca2 ablation, data show that Serca2a-deficient hearts that also lacked the central β-adrenergic signaling-dependent Serca2a negative regulator, phospholamban (PLN), had severe diastolic dysfunction that could still be corrected by β-adrenergic stimulation. Notably, integrating a serines 23/24-to-alanine PKA-refractory sarcomere incorporated cTnI molecular switch complex in mice deficient in Serca2 showed blunting of β-adrenergic stimulation-mediated enhanced diastolic heart performance. Taken together, these data provide evidence of a compensatory regulatory role of the myofilaments as a critical physiological bridging mechanism to aid in preserving heart diastolic performance in failing hearts with severe Ca2+ handling deficits.

A compartmentalized mathematical model of the β1- and β2-adrenergic signaling systems in ventricular myocytes from mouse in heart failure

Mouse models of heart failure are extensively used to research human cardiovascular diseases. In particular, one of the most common is the mouse model of heart failure resulting from transverse aortic constriction (TAC). Despite this, there are no comprehensive compartmentalized mathematical models that describe the complex behavior of the action potential, [Ca²⁺]_i transients, and their regulation by β₁- and β₂-adrenergic signaling systems in failing mouse myocytes. In this paper, we develop a novel compartmentalized mathematical model of failing mouse ventricular myocytes after TAC procedure. The model describes well the cell geometry, action potentials, [Ca²⁺]_i transients, and β₁- and β₂-adrenergic signaling in the failing cells. Simulation results obtained with the failing cell model are compared with those from the normal ventricular myocytes. Exploration of the model reveals the sarcoplasmic reticulum Ca²⁺ load mechanisms in failing ventricular myocytes. We also show a larger susceptibility of the failing myocytes to early and delayed afterdepolarizations and to a proarrhythmic behavior of Ca²⁺ dynamics upon stimulation with isoproterenol. The mechanisms of the proarrhythmic behavior suppression are investigated and sensitivity analysis is performed. The developed model can explain the existing experimental data on failing mouse ventricular myocytes and make experimentally testable predictions of a failing myocyte's behavior.

Sarcoplasmic Reticulum Calcium Release Is Required for Arrhythmogenesis in the Mouse

Dysfunctional sarcoplasmic reticulum Ca²⁺ handling is commonly observed in heart failure, and thought to contribute to arrhythmogenesis through several mechanisms. Some time ago we developed a cardiomyocyte-specific inducible SERCA2 knockout mouse, which is remarkable in the degree to which major adaptations to sarcolemmal Ca²⁺ entry and efflux overcome the deficit in SR reuptake to permit relatively normal contractile function. Conventionally, those adaptations would also be expected to dramatically increase arrhythmia susceptibility. However, that susceptibility has never been tested, and it is possible that the very rapid repolarization of the murine action potential (AP) allows for large changes in sarcolemmal Ca²⁺ transport without substantially disrupting electrophysiologic stability. We investigated this hypothesis through telemetric ECG recording in the SERCA2-KO mouse, and patch-clamp electrophysiology, Ca²⁺ imaging, and mathematical modeling of isolated SERCA2-KO myocytes. While the SERCA2-KO animals exhibit major (and unique) electrophysiologic adaptations at both the organ and cell levels, they remain resistant to arrhythmia. A marked increase in peak L-type calcium (I_CaL) current and slowed I_CaL decay elicited pronounced prolongation of initial repolarization, but faster late repolarization normalizes overall AP duration. Early afterdepolarizations were seldom observed in KO animals, and those that were observed exhibited a mechanism intermediate between murine and large mammal dynamical properties. As expected, spontaneous SR Ca²⁺ sparks and waves were virtually absent. Together these findings suggest that intact SR Ca²⁺ handling is an absolute requirement for triggered arrhythmia in the mouse, and that in its absence, dramatic changes to the major inward currents can be resisted by the substantial K⁺ current reserve, even at end-stage disease.

Cardiac muscle regulatory units are predicted to interact stronger than neighboring cross-bridges

Strong interactions between cross-bridges (XB) and regulatory units (RU) lead to a steep response of cardiac muscle to an increase in intracellular calcium. We developed a model to quantitatively assess the influence of different types of interactions within the sarcomere on the properties of cardiac muscle. In the model, the ensembles consisting of cross-bridge groups connected by elastic tropomyosin are introduced, and their dynamics is described by a set of partial differential equations. Through large scans in the free energy landscape, we demonstrate the different influence of RU-RU, XB-XB, and XB-RU interactions on the cooperativity coefficient of calcium binding, developed maximal force, and calcium sensitivity. The model solution was fitted to reproduce experimental data on force development during isometric contraction, shortening in physiological contraction, and ATP consumption by acto-myosin. On the basis of the fits, we quantified the free energy change introduced through RU-RU and XB-XB interactions and showed that RU-RU interaction leads to ~ 5 times larger change in the free energy profile of the reaction than XB-XB interaction. Due to the deterministic description of muscle contraction and its thermodynamic consistency, we envision that the developed model can be used to study heart muscle biophysics on tissue and organ levels.

Mathematical model for β1-adrenergic regulation of the mouse ventricular myocyte contraction

The β1-adrenergic regulation of cardiac myocyte contraction plays an important role in regulating heart function. Activation of this system leads to an increased heart rate and stronger myocyte contraction. However, chronic stimulation of the β1-adrenergic signaling system can lead to cardiac hypertrophy and heart failure. To understand the mechanisms of action of β1-adrenoceptors, a mathematical model of cardiac myocyte contraction that includes the β1-adrenergic system was developed and studied. The model was able to simulate major experimental protocols for measurements of steady-state force-calcium relationships, cross-bridge release rate and force development rate, force-velocity relationship, and force redevelopment rate. It also reproduced quite well frequency and isoproterenol dependencies for intracellular Ca2+ concentration ([Ca2+]i) transients, total contraction force, and sarcomere shortening. The mathematical model suggested the mechanisms of increased contraction force and myocyte shortening on stimulation of β1-adrenergic receptors is due to phosphorylation of troponin I and myosin-binding protein C and increased [Ca2+]i transient resulting from activation of the β1-adrenergic signaling system. The model was used to simulate work-loop contractions and estimate the power during the cardiac cycle as well as the effects of 4-aminopyridine and tedisamil on the myocyte contraction. The developed mathematical model can be used further for simulations of contraction of ventricular myocytes from genetically modified mice and myocytes from mice with chronic cardiac diseases.

NEW & NOTEWORTHY A new mathematical model of mouse ventricular myocyte contraction that includes the β1-adrenergic system was developed. The model simulated major experimental protocols for myocyte contraction and predicted the effects of 4-aminopyridine and tedisamil on the myocyte contraction. The model also allowed for simulations of work-loop contractions and estimation of the power during the cardiac cycle.

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.

Isolated perfused working hearts provide valuable additional information during phenotypic assessment of the diabetic mouse heart

Although murine models for studying the development of cardiac dysfunction in diabetes mellitus are well established, their reported cardiac phenotypes vary. These reported divergences may, in addition to the severity of different models, also be linked to the methods used for cardiac functional assessment. In the present study, we examined the functional changes using conventional transthoracic echocardiography (in vivo) and isolated heart perfusion techniques (ex vivo), in hearts from two mouse models; one with an overt type 2 diabetes (the db/db mouse) and one with a prediabetic state, where obesity was induced by a high-fat diet (HFD). Analysis of left ventricular function in the isolated working hearts from HFD-fed mice, suggested that these hearts develop diastolic dysfunction with preserved systolic function. Accordingly, in vivo examination demonstrated maintained systolic function, but we did not find parameters of diastolic function to be altered. In db/db mice, ex vivo working hearts showed both diastolic and systolic dysfunction. Although in vivo functional assessment revealed signs of diastolic dysfunction, the hearts did not display reduced systolic function. The contrasting results between ex vivo and in vivo function could be due to systemic changes that may sustain in vivo function, or a lack of sensitivity using conventional transthoracic echocardiography. Thus, this study demonstrates that the isolated perfused working heart preparation provides unique additional information related to the development of cardiomyopathy, which might otherwise go unnoticed when only using conventional echocardiographic assessment.

Distinct physiological effects of β1- and β2-adrenoceptors in mouse ventricular myocytes: insights from a compartmentalized mathematical model

The β1- and β2-adrenergic signaling systems play different roles in the functioning of cardiac cells. Experimental data show that the activation of the β1-adrenergic signaling system produces significant inotropic, lusitropic, and chronotropic effects in the heart, whereas the effects of the β2-adrenergic signaling system is less apparent. In this paper, a comprehensive compartmentalized experimentally based mathematical model of the combined β1- and β2-adrenergic signaling systems in mouse ventricular myocytes is developed to simulate the experimental findings and make testable predictions of the behavior of the cardiac cells under different physiological conditions. Simulations describe the dynamics of major signaling molecules in different subcellular compartments; kinetics and magnitudes of phosphorylation of ion channels, transporters, and Ca2+ handling proteins; modifications of action potential shape and duration; and [Ca2+]i and [Na+]i dynamics upon stimulation of β1- and β2-adrenergic receptors (β1- and β2-ARs). The model reveals physiological conditions when β2-ARs do not produce significant physiological effects and when their effects can be measured experimentally. Simulations demonstrated that stimulation of β2-ARs with isoproterenol caused a marked increase in the magnitude of the L-type Ca2+ current, [Ca2+]i transient, and phosphorylation of phospholamban only upon additional application of pertussis toxin or inhibition of phosphodiesterases of type 3 and 4. The model also made testable predictions of the changes in magnitudes of [Ca2+]i and [Na+]i fluxes, the rate of decay of [Na+]i concentration upon both combined and separate stimulation of β1- and β2-ARs, and the contribution of phosphorylation of PKA targets to the changes in the action potential and [Ca2+]i transient.

Simulation of the effects of moderate stimulation/inhibition of the β1-adrenergic signaling system and its components in mouse ventricular myocytes

The β1-adrenergic signaling system is one of the most important protein signaling systems in cardiac cells. It regulates cardiac action potential duration, intracellular Ca2+concentration ([Ca2+]i) transients, and contraction force. In this paper, a comprehensive experimentally based mathematical model of the β1-adrenergic signaling system for mouse ventricular myocytes is explored to simulate the effects of moderate stimulations of β1-adrenergic receptors (β1-ARs) on the action potential, Ca2+and Na+dynamics, as well as the effects of inhibition of protein kinase A (PKA) and phosphodiesterase of type 4 (PDE4). Simulation results show that the action potential prolongations reach saturating values at relatively small concentrations of isoproterenol (∼0.01 μM), while the [Ca2+]itransient amplitude saturates at significantly larger concentrations (∼0.1–1.0 μM). The differences in the response of Ca2+and Na+fluxes to moderate stimulation of β1-ARs are also observed. Sensitivity analysis of the mathematical model is performed and the model limitations are discussed. The investigated model reproduces most of the experimentally observed effects of moderate stimulation of β1-ARs, PKA, and PDE4 inhibition on the L-type Ca2+current, [Ca2+]itransients, and the sarcoplasmic reticulum Ca2+load and makes testable predictions for the action potential duration and [Ca2+]itransients as functions of isoproterenol concentration.

A Spatially Detailed Model of Isometric Contraction Based on Competitive Binding of Troponin I Explains Cooperative Interactions between Tropomyosin and Crossbridges

Biophysical models of cardiac tension development provide a succinct representation of our understanding of force generation in the heart. The link between protein kinetics and interactions that gives rise to high cooperativity is not yet fully explained from experiments or previous biophysical models. We propose a biophysical ODE-based representation of cross-bridge (XB), tropomyosin and troponin within a contractile regulatory unit (RU) to investigate the mechanisms behind cooperative activation, as well as the role of cooperativity in dynamic tension generation across different species. The model includes cooperative interactions between regulatory units (RU-RU), between crossbridges (XB-XB), as well more complex interactions between crossbridges and regulatory units (XB-RU interactions). For the steady-state force-calcium relationship, our framework predicts that: (1) XB-RU effects are key in shifting the half-activation value of the force-calcium relationship towards lower [Ca²⁺], but have only small effects on cooperativity. (2) XB-XB effects approximately double the duty ratio of myosin, but do not significantly affect cooperativity. (3) RU-RU effects derived from the long-range action of tropomyosin are a major factor in cooperative activation, with each additional unblocked RU increasing the rate of additional RU’s unblocking. (4) Myosin affinity for short (1–4 RU) unblocked stretches of actin of is very low, and the resulting suppression of force at low [Ca²⁺] is a major contributor in the biphasic force-calcium relationship. We also reproduce isometric tension development across mouse, rat and human at physiological temperature and pacing rate, and conclude that species differences require only changes in myosin affinity and troponin I/troponin C affinity. Furthermore, we show that the calcium dependence of the rate of tension redevelopment k_tr is explained by transient blocking of RU’s by a temporary decrease in XB-RU effects.

Force generation in cardiac muscle cells is driven by changes in calcium concentration. Relatively small changes in the calcium concentration over the course of a heart beat lead to the large changes in force required to fully contract and relax the heart. This is known as ‘cooperative activation’, and involves a complex interaction of several proteins involved in contraction. Current computer models which reproduce force generation often do not represent these processes explicitly, and stochastic approaches that do tend to require large amounts of computational power to solve, which limit the range of investigations in which they can be used. We have created an new computational model that captures the underlying physiological processes in more detail, and is more efficient than stochastic approaches, while still being able to run a large range of simulations. The model is able to explain the biological processes leading to the cooperative activation of muscle. In addition, the model reproduces how this cooperative activation translates to normal muscle function to generate force from changes in calcium across three different species.

Quantifying inter-species differences in contractile function through biophysical modelling

Animal models and measurements are frequently used to guide and evaluate clinical interventions. In this context, knowledge of inter-species differences in physiology is crucial for understanding the limitations and relevance of animal experimental assays for informing clinical applications. Extensive effort has been put into studying the structure and function of cardiac contractile proteins and how differences in these translate into the functional properties of muscles. However, integrating this knowledge into a quantitative description, formalising and highlighting inter-species differences both in the kinetics and in the regulation of physiological mechanisms, remains challenging. In this study we propose and apply a novel approach for the quantification of inter-species differences between mouse, rat and human. Assuming conservation of the fundamental physiological mechanisms underpinning contraction, biophysically based computational models are fitted to simulate experimentally recorded phenotypes from multiple species. The phenotypic differences between species are then succinctly quantified as differences in the biophysical model parameter values. This provides the potential of quantitatively establishing the human relevance of both animal-based experimental and computational models for application in a clinical context. Our results indicate that the parameters related to the sensitivity and cooperativity of calcium binding to troponin C and the activation and relaxation rates of tropomyosin/crossbridge binding kinetics differ most significantly between mouse, rat and human, while for example the reference tension, as expected, shows only minor differences between the species. Hence, while confirming expected inter-species differences in calcium sensitivity due to large differences in the observed calcium transients, our results also indicate more unexpected differences in the cooperativity mechanism. Specifically, the decrease in the unbinding rate of calcium to troponin C with increasing active tension was much lower for mouse than for rat and human. Our results also predicted crossbridge binding to be slowest in human and fastest in mouse.

Improving the stability of cardiac mechanical simulations

In the field of cardiac modeling, the mechanical action of the heart is often simulated using finite element methods. These simulations are becoming increasingly challenging as the computational domain is customized to a patient's anatomy, within which large heterogeneous tension gradients are generated via biophysical cell models which drive simulations of the cardiac pump cycle. The convergence of nonlinear solvers in simulations of large deformation mechanics depends on many factors. When extreme stress or irregular deformations are modeled, commonly used numerical methods can often fail to find a solution, which can prevent investigation of interesting parameter variations or use of models in a clinical context with high standards for robustness. This paper outlines a novel numerical method that is straightforward to implement and which significantly improves the stability of these simulations. The method involves adding a compressibility penalty to the standard incompressible formulation of large deformation mechanics. We compare the method's performance when used with both a direct discretization of the equations for incompressible solid mechanics, as well as the formulation based on an isochoric/deviatoric split of the deformation gradient. The addition of this penalty decreases the tendency for solutions to deviate from the incompressibility constraint, and significantly improves the ability of the Newton solver to find a solution. Additionally, our method maintains the expected order of convergence under mesh refinement, has nearly identical solutions for the pressure-volume relations, and stabilizes the solver to allow challenging simulations of both diastolic and systolic function on personalized patient geometries.

Computational modeling of Takotsubo cardiomyopathy: effect of spatially varying β-adrenergic stimulation in the rat left ventricle

In Takotsubo cardiomyopathy, the left ventricle shows apical ballooning combined with basal hypercontractility. Both clinical observations in humans and recent experimental work on isolated rat ventricular myocytes suggest the dominant mechanisms of this syndrome are related to acute catecholamine overload. However, relating observed differences in single cells to the capacity of such alterations to result in the extreme changes in ventricular shape seen in Takotsubo syndrome is difficult. By using a computational model of the rat left ventricle, we investigate which mechanisms can give rise to the typical shape of the ventricle observed in this syndrome. Three potential dominant mechanisms related to effects of β-adrenergic stimulation were considered: apical-basal variation of calcium transients due to differences in L-type and sarco(endo)plasmic reticulum Ca²⁺-ATPase activation, apical-basal variation of calcium sensitivity due to differences in troponin I phosphorylation, and apical-basal variation in maximal active tension due to, e.g., the negative inotropic effects of p38 MAPK. Furthermore, we investigated the interaction of these spatial variations in the presence of a failing Frank-Starling mechanism. We conclude that a large portion of the apex needs to be affected by severe changes in calcium regulation or contractile function to result in apical ballooning, and smooth linear variation from apex to base is unlikely to result in the typical ventricular shape observed in this syndrome. A failing Frank-Starling mechanism significantly increases apical ballooning at end systole and may be an important additional factor underpinning Takotsubo syndrome.

A computational pipeline for quantification of mouse myocardial stiffness parameters

The mouse is an important model for theoretical-experimental cardiac research, and biophysically based whole organ models of the mouse heart are now within reach. However, the passive material properties of mouse myocardium have not been much studied. We present an experimental setup and associated computational pipeline to quantify these stiffness properties. A mouse heart was excised and the left ventricle experimentally inflated from 0 to 1.44kPa in eleven steps, and the resulting deformation was estimated by echocardiography and speckle tracking. An in silico counterpart to this experiment was built using finite element methods and data on ventricular tissue microstructure from diffusion tensor MRI. This model assumed a hyperelastic, transversely isotropic material law to describe the force-deformation relationship, and was simulated for many parameter scenarios, covering the relevant range of parameter space. To identify well-fitting parameter scenarios, we compared experimental and simulated outcomes across the whole range of pressures, based partly on gross phenotypes (volume, elastic energy, and short- and long-axis diameter), and partly on node positions in the geometrical mesh. This identified a narrow region of experimentally compatible values of the material parameters. Estimation turned out to be more precise when based on changes in gross phenotypes, compared to the prevailing practice of using displacements of the material points. We conclude that the presented experimental setup and computational pipeline is a viable method that deserves wider application.

Insight into model mechanisms through automatic parameter fitting: a new methodological framework for model development

Background

Striking a balance between the degree of model complexity and parameter identifiability, while still producing biologically feasible simulations using modelling is a major challenge in computational biology. While these two elements of model development are closely coupled, parameter fitting from measured data and analysis of model mechanisms have traditionally been performed separately and sequentially. This process produces potential mismatches between model and data complexities that can compromise the ability of computational frameworks to reveal mechanistic insights or predict new behaviour. In this study we address this issue by presenting a generic framework for combined model parameterisation, comparison of model alternatives and analysis of model mechanisms.

Results

The presented methodology is based on a combination of multivariate metamodelling (statistical approximation of the input–output relationships of deterministic models) and a systematic zooming into biologically feasible regions of the parameter space by iterative generation of new experimental designs and look-up of simulations in the proximity of the measured data. The parameter fitting pipeline includes an implicit sensitivity analysis and analysis of parameter identifiability, making it suitable for testing hypotheses for model reduction. Using this approach, under-constrained model parameters, as well as the coupling between parameters within the model are identified. The methodology is demonstrated by refitting the parameters of a published model of cardiac cellular mechanics using a combination of measured data and synthetic data from an alternative model of the same system. Using this approach, reduced models with simplified expressions for the tropomyosin/crossbridge kinetics were found by identification of model components that can be omitted without affecting the fit to the parameterising data. Our analysis revealed that model parameters could be constrained to a standard deviation of on average 15% of the mean values over the succeeding parameter sets.

Conclusions

Our results indicate that the presented approach is effective for comparing model alternatives and reducing models to the minimum complexity replicating measured data. We therefore believe that this approach has significant potential for reparameterising existing frameworks, for identification of redundant model components of large biophysical models and to increase their predictive capacity.

Impaired left ventricular mechanical and energetic function in mice after cardiomyocyte-specific excision of Serca2

Sarco(endo)plasmic reticulum Ca2+ -ATPase (SERCA)2 transports Ca2+ from the cytosol into the sarcoplasmic reticulum of cardiomyocytes and is essential for maintaining myocardial Ca2+ handling and thus the mechanical function of the heart. SERCA2 is a major ATP consumer in excitation-contraction coupling but is regarded to contribute to energetically efficient Ca2+ handling in the cardiomyocyte. Previous studies using cardiomyocyte-specific SERCA2 knockout (KO) mice have demonstrated that decreased SERCA2 activity reduces the Ca2+ transient amplitude and induces compensatory Ca2+ transport mechanisms that may lead to more inefficient Ca2+ transport. In this study, we examined the relationship between left ventricular (LV) function and myocardial O2 consumption (MVo2) in ex vivo hearts from SERCA2 KO mice to directly measure how SERCA2 elimination influences mechanical and energetic features of the heart. Ex vivo hearts from SERCA2 KO hearts developed mechanical dysfunction at 4 wk and demonstrated virtually no working capacity at 7 wk. In accordance with the reported reduction in Ca2+ transient amplitude in cardiomyocytes from SERCA2 KO mice, work-independent MVo2 was decreased due to a reduced energy cost of excitation-contraction coupling. As these hearts also showed a marked impairment in the efficiency of chemomechanical energy transduction (contractile efficiency, i.e, work-dependent MVo2), hearts from SERCA2 KO mice were found to be mechanically inefficient. This ex vivo evaluation of mechanical and energetic function in hearts from SERCA2 KO mice brings together findings from previous experimental and mathematical modeling-based studies and demonstrates that reduced SERCA2 activity not only leads to mechanical dysfunction but also to energetic dysfunction.

Integrating multi-scale data to create a virtual physiological mouse heart

While the virtual physiological human (VPH) project has made great advances in human modelling, many of the tools and insights developed as part of this initiative are also applicable for facilitating mechanistic understanding of the physiology of a range of other species. This process, in turn, has the potential to provide human relevant insights via a different scientific path. Specifically, the increasing use of mice in experimental research, not yet fully complemented by a similar increase in computational modelling, is currently missing an important opportunity for using and interpreting this growing body of experimental data to improve our understanding of cardiac function. This overview describes our work to address this issue by creating a virtual physiological mouse model of the heart. We describe the similarities between human- and mouse-focused modelling, including the reuse of VPH tools, and the development of methods for investigating parameter sensitivity that are applicable across species. We show how previous results using this approach have already provided important biological insights, and how these can also be used to advance VPH heart models. Finally, we show an example application of this approach to test competing multi-scale hypotheses by investigating variations in length-dependent properties of cardiac muscle.