Creating a Structurally Realistic Finite Element Geometric Model of a Cardiomyocyte to Study the Role of Cellular Architecture in Cardiomyocyte Systems Biology

Computational Modeling and Imaging of the Intracellular Oxygen Gradient

Computational modeling can provide a mechanistic and quantitative framework for describing intracellular spatial heterogeneity of solutes such as oxygen partial pressure (pO₂). This study develops and evaluates a finite-element model of oxygen-consuming mitochondrial bioenergetics using the COMSOL Multiphysics program. The model derives steady-state oxygen (O₂) distributions from Fickian diffusion and Michaelis–Menten consumption kinetics in the mitochondria and cytoplasm. Intrinsic model parameters such as diffusivity and maximum consumption rate were estimated from previously published values for isolated and intact mitochondria. The model was compared with experimental data collected for the intracellular and mitochondrial pO₂ levels in human cervical cancer cells (HeLa) in different respiratory states and under different levels of imposed pO₂. Experimental pO₂ gradients were measured using lifetime imaging of a Förster resonance energy transfer (FRET)-based O₂ sensor, Myoglobin-mCherry, which offers in situ real-time and noninvasive measurements of subcellular pO₂ in living cells. On the basis of these results, the model qualitatively predicted (1) the integrated experimental data from mitochondria under diverse experimental conditions, and (2) the impact of changes in one or more mitochondrial processes on overall bioenergetics.

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.

Insights on the impact of mitochondrial organisation on bioenergetics in high-resolution computational models of cardiac cell architecture

Recent electron microscopy data have revealed that cardiac mitochondria are not arranged in crystalline columns but are organised with several mitochondria aggregated into columns of varying sizes spanning the cell cross-section. This raises the question—how does the mitochondrial arrangement affect the metabolite distributions within cardiomyocytes and what is its impact on force dynamics? Here, we address this question by employing finite element modeling of cardiac bioenergetics on computational meshes derived from electron microscope images. Our results indicate that heterogeneous mitochondrial distributions can lead to significant spatial variation across the cell in concentrations of inorganic phosphate, creatine (Cr) and creatine phosphate (PCr). However, our model predicts that sufficient activity of the creatine kinase (CK) system, coupled with rapid diffusion of Cr and PCr, maintains near uniform ATP and ADP ratios across the cell cross sections. This homogenous distribution of ATP and ADP should also evenly distribute force production and twitch duration with contraction. These results suggest that the PCr shuttle and associated enzymatic reactions act to maintain uniform force dynamics in the cell despite the heterogeneous mitochondrial organization. However, our model also predicts that under hypoxia activity of mitochondrial CK enzymes and diffusion of high-energy phosphate compounds may be insufficient to sustain uniform ATP/ADP distribution and hence force generation.

Mammalian cardiomyocytes contain a high volume of mitochondria, which maintains the continuous and bulk supply of ATP to sustain normal heart function. Previously, cardiac mitochondria were understood to be distributed in a regular, crystalline pattern, which facilitated a steady supply of ATP at different workloads. Using electron microscopy images of cell cross sections, we recently found that they are not regularly distributed inside cardiomyocytes. We created new spatially accurate computational models of cardiac cell bioenergetics and tested whether this heterogeneous distribution of mitochondria causes non-uniform energy supply and contractile force production in the cardiomyocyte. We found that ATP and ADP concentrations remain uniform throughout the cell because of the activity of creatine kinase (CK) enzymes that convert ATP produced in the mitochondria into creatine phosphate. Creatine phosphate rapidly diffuses to the myofibril region where it can be converted back to ATP for the contraction cycle in a timely manner. This mechanism is called the phosphocreatine shuttle (PCr shuttle). The PCr shuttle ensures that different areas of the cell produce the same amount of force regardless of the mitochondrial distribution. However, our model also shows that when the cellular oxygen supply is limited—as can be the case in conditions such as heart failure—the PCr shuttle cannot maintain uniform ATP and ADP concentrations across the cell. This causes a non-uniform acto-myosin force distribution and non-uniform twitch duration across the cell cross section. Our study suggests that mechanisms other than the PCr shuttle may be necessary to maintain uniform supply of ATP in a hypoxic environment.