Multiscale modeling in rodent ventricular myocytes

Multi-Scale Computational Modeling of Spatial Calcium Handling From Nanodomain to Whole-Heart: Overview and Perspectives

Regulation of intracellular calcium is a critical component of cardiac electrophysiology and excitation-contraction coupling. The calcium spark, the fundamental element of the intracellular calcium transient, is initiated in specialized nanodomains which co-locate the ryanodine receptors and L-type calcium channels. However, calcium homeostasis is ultimately regulated at the cellular scale, by the interaction of spatially separated but diffusively coupled nanodomains with other sub-cellular and surface-membrane calcium transport channels with strong non-linear interactions; and cardiac electrophysiology and arrhythmia mechanisms are ultimately tissue-scale phenomena, regulated by the interaction of a heterogeneous population of coupled myocytes. Recent advances in imaging modalities and image-analysis are enabling the super-resolution reconstruction of the structures responsible for regulating calcium homeostasis, including the internal structure of nanodomains themselves. Extrapolating functional and imaging data from the nanodomain to the whole-heart is non-trivial, yet essential for translational insight into disease mechanisms. Computational modeling has important roles to play in relating structural and functional data at the sub-cellular scale and translating data across the scales. This review covers recent methodological advances that enable image-based modeling of the single nanodomain and whole cardiomyocyte, as well as the development of multi-scale simulation approaches to integrate data from nanometer to whole-heart. Firstly, methods to overcome the computational challenges of simulating spatial calcium dynamics in the nanodomain are discussed, including image-based modeling at this scale. Then, recent whole-cell models, capable of capturing a range of different structures (such as the T-system and mitochondria) and cellular heterogeneity/variability are discussed at two different levels of discretization. Novel methods to integrate the models and data across the scales and simulate stochastic dynamics in tissue-scale models are then discussed, enabling elucidation of the mechanisms by which nanodomain remodeling underlies arrhythmia and contractile dysfunction. Perspectives on model differences and future directions are provided throughout.

Hybrid finite element and Brownian dynamics method for charged particles

Diffusion is often the rate-determining step in many biological processes. Currently, the two main computational methods for studying diffusion are stochastic methods, such as Brownian dynamics, and continuum methods, such as the finite element method. A previous study introduced a new hybrid diffusion method that couples the strengths of each of these two methods, but was limited by the lack of interactions among the particles; the force on each particle had to be from an external field. This study further develops the method to allow charged particles. The method is derived for a general multidimensional system and is presented using a basic test case for a one-dimensional linear system with one charged species and a radially symmetric system with three charged species.

Numerical analysis of the effect of T-tubule location on calcium transient in ventricular myocytes

Intracellular calcium (Ca2+) signaling in cardiac myocytes is vital for proper functioning of the heart. Understanding the intracellular Ca2+ dynamics would give an insight into the functions of normal and diseased hearts. In the current study, spatiotemporal Ca2+ dynamics is investigated in ventricular myocytes by considering Ca2+ release and re-uptake via sarcolemma and transverse tubules (T-tubules), Ca2+ diffusion and buffering in the cytosol, and the blockade of Ca2+ activities associated with the sarcoplasmic reticulum. This study is carried out using a three dimensional (3D) geometric model of a branch of T-tubule extracted from the electron microscopy (EM) images of a partial ventricular myocyte. Mathematical modeling is done by using a system of partial differential equations involving Ca2+, buffers, and membrane channels. Numerical simulation results suggest that a lack of T-tubule structure at the vicinity of the cell surface could increase the peak time of Ca2+ concentration in myocytes. The results also show that T-tubules and mobile buffers play an important role in the regulation of Ca2+ transient in ventricular myocytes.

Parallel acceleration for modeling of calcium dynamics in cardiac myocytes

Spatial-temporal calcium dynamics due to calcium release, buffering, and re-uptaking plays a central role in studying excitation-contraction (E-C) coupling in both healthy and defected cardiac myocytes. In our previous work, partial differential equations (PDEs) had been used to simulate calcium dynamics with realistic geometries extracted from electron microscopic imaging data. However, the computational costs of such simulations are very high on a single processor. To alleviate this problem, we have accelerated the numerical simulations of calcium dynamics by using graphics processing units (GPUs). Computational performance and simulation accuracy are compared with those based on a single CPU and another popular parallel computing technique, OpenMP.

Sensitivity of rabbit ventricular action potential and Ca²⁺ dynamics to small variations in membrane currents and ion diffusion coefficients

Little is known about how small variations in ionic currents and Ca²⁺ and Na⁺ diffusion coefficients impact action potential and Ca²⁺ dynamics in rabbit ventricular myocytes. We applied sensitivity analysis to quantify the sensitivity of Shannon et al. model (Biophys. J., 2004) to 5%-10% changes in currents conductance, channels distribution, and ion diffusion in rabbit ventricular cells. We found that action potential duration and Ca²⁺ peaks are highly sensitive to 10% increase in L-type Ca²⁺ current; moderately influenced by 10% increase in Na⁺-Ca²⁺ exchanger, Na⁺-K⁺ pump, rapid delayed and slow transient outward K⁺ currents, and Cl⁻ background current; insensitive to 10% increases in all other ionic currents and sarcoplasmic reticulum Ca²⁺ fluxes. Cell electrical activity is strongly affected by 5% shift of L-type Ca²⁺ channels and Na⁺-Ca²⁺ exchanger in between junctional and submembrane spaces while Ca²⁺-activated Cl⁻-channel redistribution has the modest effect. Small changes in submembrane and cytosolic diffusion coefficients for Ca²⁺, but not in Na⁺ transfer, may alter notably myocyte contraction. Our studies highlight the need for more precise measurements and further extending and testing of the Shannon et al. model. Our results demonstrate usefulness of sensitivity analysis to identify specific knowledge gaps and controversies related to ventricular cell electrophysiology and Ca²⁺ signaling.

New software developments for quality mesh generation and optimization from biomedical imaging data

In this paper we present a new software toolkit for generating and optimizing surface and volumetric meshes from three-dimensional (3D) biomedical imaging data, targeted at image-based finite element analysis of some biomedical activities in a single material domain. Our toolkit includes a series of geometric processing algorithms including surface re-meshing and quality-guaranteed tetrahedral mesh generation and optimization. All methods described have been encapsulated into a user-friendly graphical interface for easy manipulation and informative visualization of biomedical images and mesh models. Numerous examples are presented to demonstrate the effectiveness and efficiency of the described methods and toolkit.

Modeling effects of L-type ca(2+) current and na(+)-ca(2+) exchanger on ca(2+) trigger flux in rabbit myocytes with realistic T-tubule geometries

The transverse tubular system of rabbit ventricular myocytes consists of cell membrane invaginations (t-tubules) that are essential for efficient cardiac excitation-contraction coupling. In this study, we investigate how t-tubule micro-anatomy, L-type Ca²⁺ channel (LCC) clustering, and allosteric activation of Na⁺/Ca²⁺ exchanger by L-type Ca²⁺ current affects intracellular Ca²⁺ dynamics. Our model includes a realistic 3D geometry of a single t-tubule and its surrounding half-sarcomeres for rabbit ventricular myocytes. The effects of spatially distributed membrane ion-transporters (LCC, Na⁺/Ca²⁺ exchanger, sarcolemmal Ca²⁺ pump, and sarcolemmal Ca²⁺ leak), and stationary and mobile Ca²⁺ buffers (troponin C, ATP, calmodulin, and Fluo-3) are also considered. We used a coupled reaction-diffusion system to describe the spatio-temporal concentration profiles of free and buffered intracellular Ca²⁺. We obtained parameters from voltage-clamp protocols of L-type Ca²⁺ current and line-scan recordings of Ca²⁺ concentration profiles in rabbit cells, in which the sarcoplasmic reticulum is disabled. Our model results agree with experimental measurements of global Ca²⁺ transient in myocytes loaded with 50 μM Fluo-3. We found that local Ca²⁺ concentrations within the cytosol and sub-sarcolemma, as well as the local trigger fluxes of Ca²⁺ crossing the cell membrane, are sensitive to details of t-tubule micro-structure and membrane Ca²⁺ flux distribution. The model additionally predicts that local Ca²⁺ trigger fluxes are at least threefold to eightfold higher than the whole-cell Ca²⁺ trigger flux. We found also that the activation of allosteric Ca²⁺-binding sites on the Na⁺/Ca²⁺ exchanger could provide a mechanism for regulating global and local Ca²⁺ trigger fluxes in vivo. Our studies indicate that improved structural and functional models could improve our understanding of the contributions of L-type and Na⁺/Ca²⁺ exchanger fluxes to intracellular Ca²⁺ dynamics.

Hybrid finite element and Brownian dynamics method for diffusion-controlled reactions

Diffusion is often the rate determining step in many biological processes. Currently, the two main computational methods for studying diffusion are stochastic methods, such as Brownian dynamics, and continuum methods, such as the finite element method. This paper proposes a new hybrid diffusion method that couples the strengths of each of these two methods. The method is derived for a general multidimensional system, and is presented using a basic test case for 1D linear and radially symmetric diffusion systems.

Multi-Scale Continuum Modeling of Biological Processes: From Molecular Electro-Diffusion to Sub-Cellular Signaling Transduction

This article provides a brief review of multi-scale modeling at the molecular to cellular scale, with new results for heart muscle cells. A finite element-based simulation package (SMOL) was used to investigate the signaling transduction at molecular and sub-cellular scales (http://mccammon.ucsd.edu/smol/, http://FETK.org) by numerical solution of time-dependent Smoluchowski equations and a reaction-diffusion system. At the molecular scale, SMOL has yielded experimentally-validated estimates of the diffusion-limited association rates for the binding of acetylcholine to mouse acetylcholinesterase using crystallographic structural data. The predicted rate constants exhibit increasingly delayed steady-state times with increasing ionic strength and demonstrate the role of an enzyme's electrostatic potential in influencing ligand binding. At the sub-cellular scale, an extension of SMOL solves a non-linear, reaction-diffusion system describing Ca²⁺ ligand buffering and diffusion in experimentally-derived rodent ventricular myocyte geometries. Results reveal the important role for mobile and stationary Ca²⁺ buffers, including Ca²⁺ indicator dye. We found that the alterations in Ca²⁺-binding and dissociation rates of troponin C (TnC) and total TnC concentration modulate subcellular Ca²⁺ signals. Model predicts that reduced off-rate in whole troponin complex (TnC, TnI, TnT) versus reconstructed thin filaments (Tn, Tm, actin) alters cytosolic Ca²⁺ dynamics under control conditions or in disease-linked TnC mutations. The ultimate goal of these studies is to develop scalable methods and theories for integration of molecular-scale information into simulations of cellular-scale systems.

A localized meshless approach for modeling spatial-temporal calcium dynamics in ventricular myocytes

Spatial–temporal calcium dynamics due to calcium release, buffering and re-uptaking plays a central role in studying excitation–contraction (E–C) coupling in both normal and diseased cardiac myocytes. In this paper, we employ a meshless method, namely, the local radial basis function collocation method (LRBFCM), to model such calcium behaviors by solving a nonlinear system of reaction–diffusion partial differential equations. In particular, a simplified structural unit containing a single transverse tubule (T-tubule) and its surrounding half sarcomeres is investigated using the meshless method. Numerical results are compared with those generated by finite element methods, showing the capability and efficiency of the LRBFCM in modeling calcium dynamics in ventricular myocytes. The single T-tubule model is also extended to the whole-cell scale with T-tubules excluded to demonstrate the scalability of the proposed meshless method in handling very large domains. The experiments have shown that the LRBFCM is suitable to multiscale modeling of calcium dynamics in ventricular myocytes with high accuracy and efficiency.

Multiscale modeling of calcium dynamics in ventricular myocytes with realistic transverse tubules

Spatial-temporal Ca(2+) dynamics due to Ca(2+) release, buffering, and reuptaking plays a central role in studying excitation-contraction (E-C) coupling in both normal and diseased cardiac myocytes. In this paper, we employ two numerical methods, namely, the meshless method and the finite element method, to model such Ca(2+) behaviors by solving a nonlinear system of reaction-diffusion partial differential equations at two scales. In particular, a subcellular model containing several realistic transverse tubules (or t-tubules) is investigated and assumed to reside at different locations relative to the cell membrane. To this end, the Ca(2+) concentration calculated from the whole-cell modeling is adopted as part of the boundary constraint in the subcellular model. The preliminary simulations show that Ca(2+) concentration changes in ventricular myocytes are mainly influenced by calcium release from t-tubules.

Numerical analysis of Ca2+ signaling in rat ventricular myocytes with realistic transverse-axial tubular geometry and inhibited sarcoplasmic reticulum

The t-tubules of mammalian ventricular myocytes are invaginations of the cell membrane that occur at each Z-line. These invaginations branch within the cell to form a complex network that allows rapid propagation of the electrical signal, and hence synchronous rise of intracellular calcium (Ca²⁺). To investigate how the t-tubule microanatomy and the distribution of membrane Ca²⁺ flux affect cardiac excitation-contraction coupling we developed a 3-D continuum model of Ca²⁺ signaling, buffering and diffusion in rat ventricular myocytes. The transverse-axial t-tubule geometry was derived from light microscopy structural data. To solve the nonlinear reaction-diffusion system we extended SMOL software tool (http://mccammon.ucsd.edu/smol/). The analysis suggests that the quantitative understanding of the Ca²⁺ signaling requires more accurate knowledge of the t-tubule ultra-structure and Ca²⁺ flux distribution along the sarcolemma. The results reveal the important role for mobile and stationary Ca²⁺ buffers, including the Ca²⁺ indicator dye. In agreement with experiment, in the presence of fluorescence dye and inhibited sarcoplasmic reticulum, the lack of detectible differences in the depolarization-evoked Ca²⁺ transients was found when the Ca²⁺ flux was heterogeneously distributed along the sarcolemma. In the absence of fluorescence dye, strongly non-uniform Ca²⁺ signals are predicted. Even at modest elevation of Ca²⁺, reached during Ca²⁺ influx, large and steep Ca²⁺ gradients are found in the narrow sub-sarcolemmal space. The model predicts that the branched t-tubule structure and changes in the normal Ca²⁺ flux density along the cell membrane support initiation and propagation of Ca²⁺ waves in rat myocytes.

In cardiac muscle cells, calcium (Ca²⁺) is best known for its role in contraction activation. A remarkable amount of quantitative data on cardiac cell structure, ion-transporting protein distributions and intracellular Ca²⁺ dynamics has been accumulated. Various alterations in the protein distributions or cell ultra-structure are now recognized to be the primary mechanisms of cardiac dysfunction in a diverse range of common pathologies including cardiac arrhythmias and hypertrophy. Using a 3-D computational model, incorporating more realistic transverse-axial t-tubule geometry and considering geometric irregularities and inhomogeneities in the distribution of ion-transporting proteins, we analyze several important spatial and temporal features of Ca²⁺ signaling in rat ventricular myocytes. This study demonstrates that the computational models could serve as powerful tools for prediction and analyses of how the Ca²⁺ dynamics and cardiac excitation-contraction coupling are regulated under normal conditions or certain pathologies. The use of computational and mathematical approaches will help also to better understand aspects of cell functions that are not currently amenable to experimental investigation.

Cardiac cell modelling: observations from the heart of the cardiac physiome project

In this manuscript we review the state of cardiac cell modelling in the context of international initiatives such as the IUPS Physiome and Virtual Physiological Human Projects, which aim to integrate computational models across scales and physics. In particular we focus on the relationship between experimental data and model parameterisation across a range of model types and cellular physiological systems. Finally, in the context of parameter identification and model reuse within the Cardiac Physiome, we suggest some future priority areas for this field.