Search

In the intact heart, the distribution of electromechanical delay (EMD), the time interval between local depolarization and myocyte shortening onset, depends on the loading conditions. The distribution of EMD throughout the heart remains, however, unknown because current experimental techniques are unable to evaluate three-dimensional cardiac electromechanical behavior. The goal of this study was to determine the three-dimensional EMD distributions in the intact ventricles for sinus rhythm (SR) and epicardial pacing (EP) by using a new, to our knowledge, electromechanical model of the rabbit ventricles that incorporates a biophysical representation of myofilament dynamics. Furthermore, we aimed to ascertain the mechanisms that underlie the specific three-dimensional EMD distributions. The results revealed that under both conditions, the three-dimensional EMD distribution is nonuniform. During SR, EMD is longer at the epicardium than at the endocardium, and is greater near the base than at the apex. After EP, the three-dimensional EMD distribution is markedly different; it also changes with the pacing rate. For both SR and EP, late-depolarized regions were characterized with significant myofiber prestretch caused by the contraction of the early-depolarized regions. This prestretch delays myofiber-shortening onset, and results in a longer EMD, giving rise to heterogeneous three-dimensional EMD distributions., (2010 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

We develop a point model of the cardiac myofilament (MF) to simulate a wide variety of experimental muscle characterizations including Force-Ca relations and twitches under isometric, isosarcometric, isotonic, and auxotonic conditions. Complex MF behaviors are difficult to model because spatial interactions cannot be directly implemented as ordinary differential equations. We therefore allow phenomenological approximations with careful consideration to the relationships with the underlying biophysical mechanisms. We describe new formulations that avoid mean-field approximations found in most existing MF models. To increase the scope and applicability of the model, we include length- and temperature-dependent effects that play important roles in MF responses. We have also included a representation of passive restoring forces to simulate isolated cell shortening protocols. Possessing both computational efficiency and the ability to simulate a wide variety of muscle responses, the MF representation is well suited for coupling to existing cardiac cell models of electrophysiology and Ca-handling mechanisms. To illustrate this suitability, the MF model is coupled to the Chicago rabbit cardiomyocyte model. The combined model generates realistic appearing action potentials, intracellular Ca transients, and cell shortening signals. The combined model also demonstrates that the feedback effects of force on Ca binding to troponin can modify the cytosolic Ca transient.

We have developed a model of cardiac thin filament activation using an Ising model approach from equilibrium statistical physics. This model explicitly represents nearest-neighbor interactions between 26 troponin/tropomyosin units along a one-dimensional array that represents the cardiac thin filament. With transition rates chosen to match experimental data, the results show that the resulting force-pCa (F-pCa) relations are similar to Hill functions with asymmetries, as seen in experimental data. Specifically, Hill plots showing (log(F/(1-F)) vs. log [Ca]) reveal a steeper slope below the half activation point (Ca(50)) compared with above. Parameter variation studies show interplay of parameters that affect the apparent cooperativity and asymmetry in the F-pCa relations. The model also predicts that Ca binding is uncooperative for low [Ca], becomes steeper near Ca(50), and becomes uncooperative again at higher [Ca]. The steepness near Ca(50) mirrors the steep F-pCa as a result of thermodynamic considerations. The model also predicts that the correlation between troponin/tropomyosin units along the one-dimensional array quickly decays at high and low [Ca], but near Ca(50), high correlation occurs across the whole array. This work provides a simple model that can account for the steepness and shape of F-pCa relations that other models fail to reproduce.

A model of the functional release unit (FRU) in rat cardiac muscle consisting of one dihydropyridine receptor (DHPR) and eight ryanodine receptor (RyR) channels, and the volume surrounding them, is formulated. It is assumed that no spatial [Ca2+] gradients exist in this volume, and that each FRU acts independently. The model is amenable to systematic parameter studies in which FRU dynamics are simulated at the channel level using Monte Carlo methods with Ca2+ concentrations simulated by numerical integration of a coupled system of differential equations. Using stochastic methods, Ca(2+)-induced Ca2+ release (CICR) shows both high gain and graded Ca2+ release that is robust when parameters are varied. For a single DHPR opening, the resulting RyR Ca2+ release flux is insensitive to the DHPR open duration, and is determined principally by local sarcoplasmic reticulum (SR) Ca2+ load, consistent with experimental data on Ca2+ sparks. In addition, single RyR openings are effective in triggering Ca2+ release from adjacent RyRs only when open duration is long and SR Ca2+ load is high. This indicates relatively low coupling between RyRs, and suggests a mechanism that limits the regenerative spread of RyR openings. The results also suggest that adaptation plays an important modulatory role in shaping Ca2+ release duration and magnitude, but is not solely responsible for terminating Ca2+ release. Results obtained with the stochastic model suggest that high gain and gradedness can occur by the recruitment of independent FRUs without requiring spatial [Ca2+] gradients within a functional unit or cross-coupling between adjacent functional units.

We construct a detailed mathematical model for Ca2+ regulation in the ventricular myocyte that includes novel descriptions of subcellular mechanisms based on recent experimental findings: 1) the Keizer-Levine model for the ryanodine receptor (RyR), which displays adaptation at elevated Ca2+; 2) a model for the L-type Ca2+ channel that inactivates by mode switching; and 3) a restricted subspace into which the RyRs and L-type Ca2+ channels empty and interact via Ca2+. We add membrane currents from the Luo-Rudy Phase II ventricular cell model to our description of Ca2+ handling to formulate a new model for ventricular action potentials and Ca2+ regulation. The model can simulate Ca2+ transients during an action potential similar to those seen experimentally. The subspace [Ca2+] rises more rapidly and reaches a higher level (10-30 microM) than the bulk myoplasmic Ca2+ (peak [Ca2+]i approximately 1 microM). Termination of sarcoplasmic reticulum (SR) Ca2+ release is predominately due to emptying of the SR, but is influenced by RyR adaptation. Because force generation is roughly proportional to peak myoplasmic Ca2+, we use [Ca2+]i in the model to explore the effects of pacing rate on force generation. The model reproduces transitions seen in force generation due to changes in pacing that cannot be simulated by previous models. Simulation of such complex phenomena requires an interplay of both RyR adaptation and the degree of SR Ca2+ loading. This model, therefore, shows improved behavior over existing models that lack detailed descriptions of subcellular Ca2+ regulatory mechanisms.