Cooperative effects due to calcium binding by troponin and their consequences for contraction and relaxation of cardiac muscle under various conditions of mechanical loading

Optogenetic actuation in ChR2-transduced fibroblasts alter excitation-contraction coupling and mechano-electric feedback in coupled cardiomyocytes: a computational modeling study

With the help of the conventional electrical method and the growing optogenetic technology, cardiac fibroblasts (Fbs) have been verified to couple electrically with working myocytes and bring electrophysiological remodeling changes in them. The intrinsic properties of cardiac functional autoregulation represented by excitation-contraction coupling (ECC) and mechano-electric feedback (MEF) have also been extensively studied. However, the roles of optogenetic stimulation on the characteristics of ECC and MEF in cardiomyocytes (CMs) coupled with Fbs have been barely investigated. In this study, we proposed a combined model composed of three modules to explore these influences. Simulation results showed that (1) during ECC, an increased light duration (LD) strengthened the inflow of ChR2 current and prolonged action potential duration (APD), and extended durations of twitch and internal sarcomere deformation through the decreased dissociation of calcium with troponin C (CaTnC) complexes and the prolonged duration of Xb attachment-detachment; (2) during MEF, an increased LD was followed by a longer muscle twitch and deformation, and led to APD prolongation through the inward ChR2 current and its inward rectification kinetics, which far outweighed the effects of the delaying dissociation of CaTnC complexes and the prolonged reverse mode of Na⁺-Ca²⁺ exchange on AP shortening; (3) due to the ChR2 current's rectification feature, enhancing the light irradiance (LI) brought slight variations in peak or valley values of electrophysiological and mechanical parameters while did not change durations of AP and twitch and muscle deformation in both ECC and MEF. In conclusion, the inward ChR2 current and its inward rectification feature were found to affect significantly the durations of AP and twitch in both ECC and MEF. The roles of optogenetic actuation on both ECC and MEF should be considered in future cardiac computational optogenetics at the tissue and organ scale.

Insights From Computational Modeling Into the Contribution of Mechano-Calcium Feedback on the Cardiac End-Systolic Force-Length Relationship

In experimental studies on cardiac tissue, the end-systolic force-length relation (ESFLR) has been shown to depend on the mode of contraction: isometric or isotonic. The isometric ESFLR is derived from isometric contractions spanning a range of muscle lengths while the isotonic ESFLR is derived from shortening contractions across a range of afterloads. The ESFLR of isotonic contractions consistently lies below its isometric counterpart. Despite the passing of over a hundred years since the first insight by Otto Frank, the mechanism(s) underlying this protocol-dependent difference in the ESFLR remain incompletely explained. Here, we investigate the role of mechano-calcium feedback in accounting for the difference between these two ESFLRs. Previous studies have compared the dynamics of isotonic contractions to those of a single isometric contraction at a length that produces maximum force, without considering isometric contractions at shorter muscle lengths. We used a mathematical model of cardiac excitation-contraction to simulate isometric and force-length work-loop contractions (the latter being the 1D equivalent of the whole-heart pressure-volume loop), and compared Ca²⁺ transients produced under equivalent force conditions. We found that the duration of the simulated Ca²⁺ transient increases with decreasing sarcomere length for isometric contractions, and increases with decreasing afterload for work-loop contractions. At any given force, the Ca²⁺ transient for an isometric contraction is wider than that during a work-loop contraction. By driving simulated work-loops with wider Ca²⁺ transients generated from isometric contractions, we show that the duration of muscle shortening was prolonged, thereby shifting the work-loop ESFLR toward the isometric ESFLR. These observations are explained by an increase in the rate of binding of Ca²⁺ to troponin-C with increasing force. However, the leftward shift of the work-loop ESFLR does not superimpose on the isometric ESFLR, leading us to conclude that while mechano-calcium feedback does indeed contribute to the difference between the two ESFLRs, it does not completely account for it.

Mechano-calcium and mechano-electric feedbacks in the human cardiomyocyte analyzed in a mathematical model

Experiments on animal hearts (rat, rabbit, guinea pig, etc.) have demonstrated that mechano-calcium feedback (MCF) and mechano-electric feedback (MEF) are very important for myocardial self-regulation because they adjust the cardiomyocyte contractile function to various mechanical loads and to mechanical interactions between heterogeneous myocardial segments in the ventricle walls. In in vitro experiments on these animals, MCF and MEF manifested themselves in several basic classical phenomena (e.g., load dependence, length dependence of isometric twitches, etc.), and in the respective responses of calcium transients and action potentials. However, it is extremely difficult to study simultaneously the electrical, calcium, and mechanical activities of the human heart muscle in vitro. Mathematical modeling is a useful tool for exploring these phenomena. We have developed a novel model to describe electromechanical coupling and mechano-electric feedbacks in the human cardiomyocyte. It combines the ‘ten Tusscher–Panfilov’ electrophysiological model of the human cardiomyocyte with our module of myocardium mechanical activity taken from the ‘Ekaterinburg–Oxford’ model and adjusted to human data. Using it, we simulated isometric and afterloaded twitches and effects of MCF and MEF on excitation–contraction coupling. MCF and MEF were found to affect significantly the duration of the calcium transient and action potential in the human cardiomyocyte model in response to both smaller afterloads as compared to bigger ones and various mechanical interventions applied during isometric and afterloaded twitches.

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.

High tension in sarcomeres hinders myocardial relaxation: A computational study

Experiments have shown that the relaxation phase of cardiac sarcomeres during an isometric twitch is prolonged in muscles that reached a higher peak tension. However, the mechanism is not completely understood. We hypothesize that the binding of calcium to troponin is enhanced by the tension in the thin filament, thus contributing to the prolongation of contraction upon higher peak tension generation. To test this hypothesis, we developed a computational model of sarcomere mechanics that incorporates tension-dependence of calcium binding. The model was used to simulate isometric twitch experiments with time dependency in the form of a two-state cross-bridge cycle model and a transient intracellular calcium concentration. In the simulations, peak isometric twitch tension appeared to increase linearly by 51.1 KPa with sarcomere length from 1.9 μm to 2.2 μm. Experiments showed an increase of 47.3 KPa over the same range of sarcomere lengths. The duration of the twitch also increased with both sarcomere length and peak intracellular calcium concentration, likely to be induced by the inherently coupled increase of the peak tension in the thin filament. In the model simulations, the time to 50% relaxation (tR50) increased over the range of sarcomere lengths from 1.9 μm to 2.2 μm by 0.11s, comparable to the increased duration of 0.12s shown in experiments. Model simulated tR50 increased by 0.12s over the range of peak intracellular calcium concentrations from 0.87 μM to 1.45 μM. Our simulation results suggest that the prolongation of contraction at higher tension is a result of the tighter binding of Ca2+ to troponin in areas under higher tension, thus delaying the deactivation of the troponin.

A new myofilament contraction model with ATP consumption for ventricular cell model

A new contraction model of cardiac muscle was developed by combining previously described biochemical and biophysical models. The biochemical component of the new contraction model represents events in the presence of Ca2+-crossbridge attachment and power stroke following inorganic phosphate release, detachment evoked by the replacement of ADP by ATP, ATP hydrolysis, and recovery stroke. The biophysical component focuses on Ca2+ activation and force (F b) development assuming an equivalent crossbridge. The new model faithfully incorporates the major characteristics of the biochemical and biophysical models, such as F b activation by transient Ca2+ ([Ca2+]-F b), [Ca2+]-ATP hydrolysis relations, sarcomere length-F b, and F b recovery after jumps in length under the isometric mode and upon sarcomere shortening after a rapid release of mechanical load under the isotonic mode together with the load-velocity relationship. ATP consumption was obtained for all responses. When incorporated in a ventricular cell model, the contraction model was found to share approximately 60% of the total ATP usage in the cell model.

Effects of subchronic lead intoxication of rats on the myocardium contractility

Outbred male rats were repeatedly injected IP with sub-lethal doses of lead acetate 3 times a week during 5 weeks. They developed an explicit, even if moderate, lead intoxication characterized by typical hematological and some other features. The next day after the last injection the heart of each animal was excised, and the trabecules and papillary muscles from the right ventricle were used for modeling in vitro isometric (with varying starting length of the preparation) regimes of the contraction-relaxation cycle with different preloads. Several well-established parameters of this model were found changed compared with the preparations taken from the hearts of healthy control rats. Background in vivo calcium treatment attenuated both systemic and cardiotoxic effects of lead to an extent. We show for the first time that subchronic intoxication with lead caused myocardial preparations in a wide range of lengths to respond by a decrease in the time and speed parameters of the isometric contraction while maintaining its amplitude and by a decrease in the passive stiffness of trabecules. The responses of the various heart structures are outlined, and the isomyosin ratio is shown to have shifted towards the slow isoform. Mechanistic and toxicological inferences from the results are discussed.

Transmural cellular heterogeneity in myocardial electromechanics

Myocardial heterogeneity is an attribute of the normal heart. We have developed integrative models of cardiomyocytes from the subendocardial (ENDO) and subepicardial (EPI) ventricular regions that take into account experimental data on specific regional features of intracellular electromechanical coupling in the guinea pig heart. The models adequately simulate experimental data on the differences in the action potential and contraction between the ENDO and EPI cells. The modeling results predict that heterogeneity in the parameters of calcium handling and myofilament mechanics in isolated ENDO and EPI cardiomyocytes are essential to produce the differences in Ca2+ transients and contraction profiles via cooperative mechanisms of mechano-calcium-electric feedback and may further slightly modulate transmural differences in the electrical properties between the cells. Simulation results predict that ENDO cells have greater sensitivity to changes in the mechanical load than EPI cells. These data are important for understanding the behavior of cardiomyocytes in the intact heart.

The effects of load on transmural differences in contraction of isolated mouse ventricular cardiomyocytes

Mechanical properties of cardiomyocytes from different transmural regions are heterogeneous in the left ventricular wall. The cardiomyocyte mechanical environment affects this heterogeneity because of mechano-electric feedback mechanisms. In the present study, we investigated the effects of the mechanical load (preload and afterload) on transmural differences in contraction of subendocardial (ENDO) and subepicardial (EPI) single cells isolated from the murine left ventricle. Various preloads imposed via axial stretch and afterloads (unloaded and heavy loaded conditions) were applied to the cells using carbon fiber techniques for single myocytes. To simulate experimentally obtained results and to predict mechanisms underlying the cellular response to change in load, our mathematical models of the ENDO and EPI cells were used. Our major findings are the following. Our results show that ENDO and EPI cardiomyocytes have different mechanical responses to changes in preload to the cells. Under auxotonic contractions at low preload (unstretched cells), time to peak contraction (T_max) and the time constant of [Ca²⁺]_i transient decay were significantly longer in ENDO cells than in EPI cells. An increase in preload (stretched cells) prolonged T_max in both cell types; however, the prolongation was greater in EPI cells, resulting in a decrease in the transmural gradient in T_max at high preload. Comparing unloaded and heavy loaded (isometric) contractions of the cells we found that transmural gradient in the time course of contraction is independent of the loading conditions. Our mathematical cell models were able to reproduce the experimental results on the distinct cellular responses to changes in the mechanical load when we accounted for an ENDO/EPI difference in the parameters of cooperativity of calcium activation of myofilaments.

Effects of cellular electromechanical coupling on functional heterogeneity in a one-dimensional tissue model of the myocardium

Based on the experimental evidence, we developed a one-dimensional (1D) model of heterogeneous myocardial tissue consisting of in-series connected cardiomyocytes from distant transmural regions using mathematical models of subendocardial and subepicardial cells. The regional deformation patterns produced by our 1D model are consistent with the transmural regional strain patterns obtained experimentally in the normal heart in vivo. The modelling results suggest that the mechanical load may essentially affect the transmural gradients in the electrical and mechanical properties of interacting myocytes within a tissue, thereby regulating global myocardial output.

A simple model of cardiac muscle for multiscale simulation: Passive mechanics, crossbridge kinetics and calcium regulation

A simple model of cardiac muscle was designed for multiscale simulation of heart mechanics. Relaxed cardiac muscle was described as a transversally isotropic hyperelastic material. Active tension caused by actin-myosin crossbridges depends on the ensemble averaged strain of myosin heads bound to actin. Calcium activation was modeled by Ca²⁺ binding to troponin-C. To account for the dependence of troponin affinity for Ca²⁺ on myosin heads strongly bound to actin, the kinetics of troponin binding to Ca²⁺ in the overlap zone of the thin and thick filaments and outside it were separated. The changes in the length of these zones during muscle shortening or lengthening were accounted for explicitly. Simplified version of the model contains only 5 ordinary differential equations (ODE). Model parameters were estimated from a limited set of experiments with skeletal and cardiac muscle. Simulations have shown that model reproduces qualitatively a number of experimental observations: steady-state force-velocity and stiffness-velocity relations; mechanical responses to step changes in muscle length or load; steep Ca²⁺-tension relationship and its dependence on sarcomere length tension (the Frank-Starling mechanism); tension, shortening and Ca²⁺-transients in twitch isometric and isotonic contractions, tension development and redevelopment upon fast change in Ca²⁺ concentration or muscle release followed by re-stretch. We believe that the model can be effectively used for modeling contraction and relaxation of the heart.

Investigations of Molecular Mechanisms of Actin-Myosin Interactions in Cardiac Muscle

The functional characteristics of cardiac muscle depend on the composition of protein isoforms in the cardiomyocyte contractile machinery. In the ventricular myocardium of mammals, several isoforms of contractile and regulatory proteins are expressed - two isoforms of myosin (V1 and V3) and three isoforms of tropomyosin chains (α, β, and κ). Expression of protein isoforms depends on the animal species, its age and hormonal status, and this can change with pathologies of the myocardium. Mutations in these proteins can lead to cardiomyopathies. The functional significance of the protein isoform composition has been studied mainly on intact hearts or on isolated preparations of myocardium, which could not provide a clear comprehension of the role of each particular isoform. Present-day experimental techniques such as an optical trap and in vitro motility assay make it possible to investigate the phenomena of interactions of contractile and regulatory proteins on the molecular level, thus avoiding effects associated with properties of a whole muscle or muscle tissue. These methods enable free combining of the isoforms to test the molecular mechanisms of their participation in the actin-myosin interaction. Using the optical trap and the in vitro motility assay, we have studied functional characteristics of the cardiac myosin isoforms, molecular mechanisms of the calcium-dependent regulation of actin-myosin interaction, and the role of myosin and tropomyosin isoforms in the cooperativity mechanisms in myocardium. The knowledge of molecular mechanisms underlying myocardial contractility and its regulation is necessary for comprehension of cardiac muscle functioning, its disorders in pathologies, and for development of approaches for their correction.

Mechano-chemical Interactions in Cardiac Sarcomere Contraction: A Computational Modeling Study

We developed a model of cardiac sarcomere contraction to study the calcium-tension relationship in cardiac muscle. Calcium mediates cardiac contraction through its interactions with troponin (Tn) and subsequently tropomyosin molecules. Experimental studies have shown that a slight increase in intracellular calcium concentration leads to a rapid increase in sarcomeric tension. Though it is widely accepted that the rapid increase is not possible without the concept of cooperativity, the mechanism is debated. We use the hypothesis that there exists a base level of cooperativity intrinsic to the thin filament that is boosted by mechanical tension, i.e. a high level of mechanical tension in the thin filament impedes the unbinding of calcium from Tn. To test these hypotheses, we developed a computational model in which a set of three parameters and inputs of calcium concentration and sarcomere length result in output tension. Tension as simulated appeared in good agreement with experimentally measured tension. Our results support the hypothesis that high tension in the thin filament impedes Tn deactivation by increasing the energy required to detach calcium from the Tn. Given this hypothesis, the model predicted that the areas with highest tension, i.e. closest to the Z-disk end of the single overlap region, show the largest concentration of active Tn's.

Mechano-electric heterogeneity of the myocardium as a paradigm of its function

Myocardial heterogeneity is well appreciated and widely documented, from sub-cellular to organ levels. This paper reviews significant achievements of the group, led by Professor Vladimir S. Markhasin, Russia, who was one of the pioneers in studying and interpreting the relevance of cardiac functional heterogeneity.

Mechano-electric feedback in one-dimensional model of myocardium

We utilized our earlier developed 1D mathematical model of the heart muscle strand to study contribution of the bilateral interactions between excitation and contraction on the cellular and tissue levels to the local and global myocardium function. Numerical experiments on the model showed that an initially uniform strand, formed on the inherently identical cells, became functionally heterogeneous due to the asynchronous excitation via the electrical wave spread. Mechanical interactions between the cells and the mechano-electric feedback beat-to-beat affect the functional characteristics of coupled cardiomyocytes further, adjusting their electrical and mechanical heterogeneity to the activation timing. Model simulations showed that functional heterogeneity increases with an enlarged spatial extension of the myocardial strand (in terms of the longer slack length not a higher stretch of the strand), demonstrating a special role of the heart size in its function. Model analysis suggests that cooperative mechanisms of myofilament calcium activation contribute essentially to the generation of cellular functional heterogeneity in contracting cardiac tissue.

Modeling calcium regulation of contraction, energetics, signaling, and transcription in the cardiac myocyte

Calcium (Ca(2+)) plays many important regulatory roles in cardiac muscle cells. In the initial phase of the action potential, influx of Ca(2+) through sarcolemmal voltage-gated L-type Ca(2+) channels (LCCs) acts as a feed-forward signal that triggers a large release of Ca(2+) from the junctional sarcoplasmic reticulum (SR). This Ca(2+) drives heart muscle contraction and pumping of blood in a process known as excitation-contraction coupling (ECC). Triggered and released Ca(2+) also feed back to inactivate LCCs, attenuating the triggered Ca(2+) signal once release has been achieved. The process of ECC consumes large amounts of ATP. It is now clear that in a process known as excitation-energetics coupling, Ca(2+) signals exert beat-to-beat regulation of mitochondrial ATP production that closely couples energy production with demand. This occurs through transport of Ca(2+) into mitochondria, where it regulates enzymes of the tricarboxylic acid cycle. In excitation-signaling coupling, Ca(2+) activates a number of signaling pathways in a feed-forward manner. Through effects on their target proteins, these interconnected pathways regulate Ca(2+) signals in complex ways to control electrical excitability and contractility of heart muscle. In a process known as excitation-transcription coupling, Ca(2+) acting primarily through signal transduction pathways also regulates the process of gene transcription. Because of these diverse and complex roles, experimentally based mechanistic computational models are proving to be very useful for understanding Ca(2+) signaling in the cardiac myocyte.

Mechanical regulation of native and the recombinant calcium channel

L-type calcium channels are modulated by a host of mechanisms that include voltage, calcium ions (Ca(2+) dependent inactivation and facilitation), cytosolic proteins (CAM, CAMKII, PKA, PKC, etc.), and oxygen radicals. Here we describe yet another Ca(2+) channel regulatory mechanism that is induced by pressure-flow (PF) forces of ∼25dyn/cm(2) producing 35-60% inhibition of channel current. Only brief periods (300ms) of such PF pulses were required to suppress reversibly the current. Recombinant Ca(2+) channels (α1c77/β2a/α2δ and α1c77/β1/α2δ), expressed in HEK293 cells, were similarly suppressed by PF pulses. To examine whether Ca(2+) released by PF pulses triggered from different sub-cellular compartments (SR, ER, mitochondria) underlies the inhibitory effect of PF on the channel current, pharmacological agents and ionic substitutions were employed to probe this possibility. No significant difference in effectiveness of PF pulses to suppress ICa or IBa (used to inhibit CICR) was found between control cells and those exposed to U73122 and 2-APB (PLC and IP3R pathway modulators), thapsigargin and BAPTA (SERCA2a modulator), dinitrophenol, FCCP and Ru360 (mitochondrial inhibitors), l-NAME (NOS inhibitor signaling), cAMP and Pertussis toxin (Gi protein modulator). We concluded that the rapid and reversible modulation of the Ca(2+) channel by PF pulses is independent of intracellular release of Ca(2+) and Ca(2+) dependent inactivation of the channel and may represent direct mechanical regulatory effect on the channel protein in addition to previously reported Ca(2+)-release or entry dependent mechanism.

Cardiac electromechanical models: from cell to organ

The heart is a multiphysics and multiscale system that has driven the development of the most sophisticated mathematical models at the frontiers of computational physiology and medicine. This review focuses on electromechanical (EM) models of the heart from the molecular level of myofilaments to anatomical models of the organ. Because of the coupling in terms of function and emergent behaviors at each level of biological hierarchy, separation of behaviors at a given scale is difficult. Here, a separation is drawn at the cell level so that the first half addresses subcellular/single-cell models and the second half addresses organ models. At the subcellular level, myofilament models represent actin–myosin interaction and Ca-based activation. The discussion of specific models emphasizes the roles of cooperative mechanisms and sarcomere length dependence of contraction force, considered to be the cellular basis of the Frank–Starling law. A model of electrophysiology and Ca handling can be coupled to a myofilament model to produce an EM cell model, and representative examples are summarized to provide an overview of the progression of the field. The second half of the review covers organ-level models that require solution of the electrical component as a reaction–diffusion system and the mechanical component, in which active tension generated by the myocytes produces deformation of the organ as described by the equations of continuum mechanics. As outlined in the review, different organ-level models have chosen to use different ionic and myofilament models depending on the specific application; this choice has been largely dictated by compromises between model complexity and computational tractability. The review also addresses application areas of EM models such as cardiac resynchronization therapy and the role of mechano-electric coupling in arrhythmias and defibrillation.

Contribution of mechanical factors to arrhythmogenesis in calcium overloaded cardiomyocytes: model predictions and experiments

It is well-known that Ca²⁺ overload in cardiomyocytes may underlie arrhythmias. However, the possible contribution of mechanical factors to rhythm disturbances in Ca²⁺ overloaded myocytes has not been sufficiently investigated. We used a mathematical model of the electrical and mechanical activity of cardiomyocytes to reveal an essential role of the mechanisms of cardiac mechano-electric feedback in arrhythmogenesis in Ca²⁺ overloaded myocardium. In the model, the following mechanical factors increased Ca²⁺ overload in contracting cardiomyocytes and promoted rhythm disturbances: i) a decrease in the mechanical load for afterloaded contractions; and ii) a decrease in the initial length of sarcomeres for isometric twitches. In exact accordance with the model predictions, in experiments on papillary muscles from the right ventricle of guinea pigs with Ca²⁺ overloaded cardiomyocytes (using 0.5-1 μM of ouabain), we found that emergence of rhythm disturbances and extrasystoles depends on the mechanical conditions of muscle contraction.

Role of myocardial viscoelasticity in disturbances of electrical and mechanical activity in calcium overloaded cardiomyocytes: mathematical modeling

Cardiomyocyte Ca(2+) overload is closely linked to cardiac arrhythmias. We have earlier shown in a mathematical model that myocardium mechanical activity may contribute to rhythm disturbances induced by Ca(2+) overload in cardiomyocytes with reduced Na(+)-K(+) pump work (Sulman et al., 2008). The same model is used here to address possible contribution of the passive mechanical properties of cardiac muscle (i.e. myocardial viscous and elastic properties) to the arrhythmogenesis. In a series of contractions at regular pacing rate of 75 beats/min a model with higher viscosity demonstrated essentially earlier appearance of extrasystoles due to a faster cardiomyocyte Ca(2+) loading up to a level triggering spontaneous Ca(2+) releases from the sarcoplasmic reticulum. The model predicts that myocardial elasticity also may affect arrhythmogenesis in cardiomyocytes overloaded with Ca(2+). Contribution of the mechanical properties of the myocardial tissue to the arrhythmia has been analyzed for wide ranges of both viscosity and elasticity coefficients. The results suggest that myocardial viscoelastic properties may be a factor affecting Ca(2+) handling in cardiomyocytes and contributing to cardiac mechano-electric feedback in arrhythmogenesis.

Integrative modeling of the cardiac ventricular myocyte

Cardiac electrophysiology is a discipline with a rich 50-year history of experimental research coupled with integrative modeling which has enabled us to achieve a quantitative understanding of the relationships between molecular function and the integrated behavior of the cardiac myocyte in health and disease. In this paper, we review the development of integrative computational models of the cardiac myocyte. We begin with a historical overview of key cardiac cell models that helped shape the field. We then narrow our focus to models of the cardiac ventricular myocyte and describe these models in the context of their subcellular functional systems including dynamic models of voltage-gated ion channels, mitochondrial energy production, ATP-dependent and electrogenic membrane transporters, intracellular Ca dynamics, mechanical contraction, and regulatory signal transduction pathways. We describe key advances and limitations of the models as well as point to new directions for future modeling research. WIREs Syst Biol Med 2011 3 392-413 DOI: 10.1002/wsbm.122

Effects of cardiac myosin binding protein-C on the regulation of interaction of cardiac myosin with thin filament in an in vitro motility assay

Modulatory role of whole cardiac myosin binding protein-C (сMyBP-C) in regulation of cardiac muscle contractility was studied in the in vitro motility assay with rabbit cardiac myosin as a motor protein. The effects of cMyBP-C on the interaction of cardiac myosin with regulated thin filament were tested in both in vitro motility and ATPase assays. We demonstrate that the addition of cMyBP-C increases calcium regulated Mg-ATPase activity of cardiac myosin at submaximal calcium. The Hill coefficient for 'pCa-velocity' relation in the in vitro motility assay decreased and the calcium sensitivity increased when сMyBP-C was added. Results of our experiments testifies in favor of the hypothesis that сMyBP-C slows down cross-bridge kinetics when binding to actin.

Study of the interaction between rabbit cardiac contractile and regulatory proteins. An in vitro motility assay

A series of experiments was performed in an in vitro motility assay with reconstructed thin filaments to obtain pCa-force relationships for cardiac isomyosins V1 and V3. Two concentrations of each isomyosin (200 and 300 microg/ml) on the surface of a flow cell were tested. Isometric force was estimated as the amount of actin-binding protein, alpha-actinin, stopping thin filament movement. It was found that the amount of alpha-actinin stopping the movement at saturating calcium concentration for V3 was twice higher than for V1 at both concentrations of isoforms. Hill coefficients of cooperativity (h) were determined for pCa-force relationships. The value of h did not differ significantly for isoforms at 300 microg/ml of protein (h was 1.56 for V1 and 1.54 for V3). However, the Hill coefficient was higher for V3 isoform at 200 microg/ml (h = 2.00 and 1.76 for V3 and V1, respectively). Importantly, the Hill coefficient increased for both isoenzymes when their concentrations were decreased. The connection between Hill coefficient and cooperative interactions between cardiac contractile and regulatory proteins is analyzed in detail.

Mathematical modeling of mechanically modulated rhythm disturbances in homogeneous and heterogeneous myocardium with attenuated activity of na+ -k+ pump

A mathematical model of the cardiomyocyte electromechanical function is used to study contribution of mechanical factors to rhythm disturbances in the case of the cardiomyocyte calcium overload. Particular attention is paid to the overload caused by diminished activity of the sodium-potassium pump. It is shown in the framework of the model, where mechano-calcium feedback is accounted for that myocardium mechanics may significantly enhance arrhythmogenicity of the calcium overload. Specifically, a role of cross-bridge attachment/detachment processes, a role of mechanical conditions of myocardium contractions (length, load), and a role of myocardium viscosity in the case of simulated calcium overload have been revealed. Underlying mechanisms are analyzed. Several approaches are designed in the model and compared to each other for recovery of the valid myocardium electrical and mechanical performance in the case of the partially suppressed sodium-potassium pump.

Activation sequence as a key factor in spatio-temporal optimization of myocardial function

Using one-dimensional models of myocardial tissue, implemented as chains of virtual ventricular muscle segments that are kinematically connected in series, we studied the role of the excitation sequence in spatio-temporal organization of cardiac function. Each model element was represented by a well-verified mathematical model of cardiac electro-mechanical activity. We found that homogeneous chains, consisting of identical elements, respond to non-simultaneous stimulation by generation of complex spatio-temporal heterogeneities in element deformation. These are accompanied by the establishment of marked gradients in local electro-mechanical properties of the elements (heterogeneity in action potential duration, Ca2+ transient characteristics and sarcoplasmic reticulum Ca2+ loading). In heterogeneous chains, composed of elements simulating fast and slow contracting cardiomyocytes from different transmural layers, we found that only activation sequences where stimulation of the slower elements preceded that of faster ones gave rise to optimization of the system's electro-mechanical function, which was confirmed experimentally. Based on the results obtained, we hypothesize that the sequence of activation of cardiomyocytes in different ventricular layers is one of the key factors of spatio-temporal organization of myocardium. Moreover, activation sequence and regional differences in intrinsic electro-mechanical properties of cardiac muscle must be matched in order to optimize myocardial function.

Hybrid duplex: a novel method to study the contractile function of heterogeneous myocardium

In an earlier study, we experimentally mimicked the effects of mechanical interaction between different regions of the ventricular wall by allowing pairs of independently maintained cardiac muscle fibers to interact mechanically in series or in parallel. This simple physiological model of heterogeneous myocardium, which has been termed “duplex,” has provided new insight into basic effects of cardiac electromechanical heterogeneity. Here, we present a novel “hybrid duplex,” where one of the elements is an isolated cardiac muscle and the other a “virtual cardiac muscle.” The virtual muscle is represented by a computational model of cardiomyocyte electromechanical activity. We present in detail the computer-based digital control system that governs the mechanical interaction between virtual and biological muscle, the software used for data analysis, and working implementations of the model. Advantages of the hybrid duplex method are discussed, and experimental recordings are presented for illustration and as proof of the principle.

Influence of viscosity on myocardium mechanical activity: a mathematical model

We have previously proposed and validated a mathematical model of myocardium contraction-relaxation cycle based on current knowledge of regulatory role of Ca2+ and cross-bridge kinetics in cardiac cell. That model did not include viscous elements. Here we propose a modification of the model, in which two viscous elements are added, one in parallel to the contractile element, and one more in parallel to the series elastic element. The modified model allowed us to simulate and explain some subtle experimental data on relaxation velocity in isotonic twitches and on a mismatch between the time course of sarcomere shortening/lengthening and the time course of active force generation in isometric twitches. Model results were compared with experimental data obtained from 28 rat LV papillary muscles contracting and relaxing against various loads. Additional model analysis suggested contribution of viscosity to main inotropic and lusitropic characteristics of myocardium performance.

Mechano-electric interactions in heterogeneous myocardium: development of fundamental experimental and theoretical models

The heart is structurally and functionally a highly non-homogenous organ, yet its main function as a pump can only be achieved by the co-ordinated contraction of millions of ventricular cells. This apparent contradiction gives rise to the hypothesis that 'well-organised' inhomogeneity may be a pre-requisite for normal cardiac function. Here, we present a set of novel experimental and theoretical tools for the study of this concept. Heterogeneity, in its most condensed form, can be simulated using two individually controlled, mechanically interacting elements (duplex). We have developed and characterised three different types of duplexes: (i) biological duplex, consisting of two individually perfused biological samples (like thin papillary muscles or a trabeculae), (ii) virtual duplex, made-up of two interacting mathematical models of cardiac muscle, and (iii) hybrid duplex, containing a biological sample that interacts in real-time with a virtual muscle. In all three duplex types, in-series or in-parallel mechanical interaction of elements can be studied during externally isotonic, externally isometric, and auxotonic modes of contraction and relaxation. Duplex models, therefore, mimic (patho-)physiological mechano-electric interactions in heterogeneous myocardium at the multicellular level, and in an environment that allows one to control mechanical, electrical and pharmacological parameters. Results obtained using the duplex method show that: (i) contractile elements in heterogeneous myocardium are not 'independent' generators of tension/shortening, as their ino- and lusitropic characteristics change dynamically during mechanical interaction-potentially matching microscopic contractility to macroscopic demand, (ii) mechanical heterogeneity contributes differently to action potential duration (APD) changes, depending on whether mechanical coupling of elements is in-parallel or in-series, which may play a role in mechanical tuning of distant tissue regions, (iii) electro-mechanical activity of mechanically interacting contractile elements is affected by their activation sequence, which may optimise myocardial performance by smoothing intrinsic differences in APD. In conclusion, we present a novel set of tools for the experimental and theoretical investigation of cardiac mechano-electric interactions in healthy and/or diseased heterogeneous myocardium, which allows for the testing of previously inaccessible concepts.

Load dependence of ventricular performance explained by model of calcium-myofilament interactions

Although a simple concept of load-independent behavior of the intact heart evolved from early studies of isolated, intact blood-perfused hearts, more recent studies showed that, as in isolated muscle, the mode of contraction (isovolumic vs. ejection) impacts on end-systolic elastance. The purpose of the present study was to test whether a four-state model of myofilament interactions with length-dependent rate constants could explain the complex contractile behavior of the intact, ejecting heart. Studies were performed in isolated, blood-perfused canine hearts with intracellular calcium transients measured by macroinjected aequorin. Measured calcium transients were used as the driving function for the model, and length-dependent rate constants yielding the highest concordance between measured and model-predicted midwall stress at different isovolumic volumes were determined. These length-dependent rate constants successfully predicted contractile behavior on ejecting contractions. This, along with additional model analysis, suggests that length-dependent changes in calcium binding affinity may not be an important factor contributing to load-dependent contractile performance in the intact heart under physiological conditions.

Controlled auxotonic twitch in papillary muscle: a new computer-based control approach

Based on new advancements in digital technology, we developed a PC- and DSP-based measurement and control system for isolated papillary muscle experiments. High flexibility was obtained through a three level control. Length or force was controlled real-time with a sample frequency of 5000 Hz. Muscle length and up to three segment lengths were measured simultaneously and each of these lengths could be chosen as feedback variable. Individual algorithms were implemented for different twitch types. Batches of twitches were organized in experiment protocols. The system included a new twitch type, namely a controlled auxotonic twitch. In this twitch, the muscle acted against a virtual ideal spring, giving a proportional change in developed force and shortening. The value of the virtual spring constant could be set on-line or defined in the experiment protocol. An increasing virtual spring constant represented a smooth transition from isotonic to isometric conditions.