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

Measuring and Modelling the Effect of Inorganic Phosphate on Cross-bridge Mechanics in Human Cardiac Muscle

Many common chronic diseases operate at the intersection of metabolic and cardiovascular dysfunction. In order to model the effects of these diseases and investigate underlying causes we are developing a cardiomyocyte model which incorporates both the mechanics and metabolic factors that underlie work done by the heart. In this paper we present the first experimental results from our study measuring mechanical properties in human cardiac trabeculae, including the effect of inorganic phosphate (Pi) on the complex modulus at 37 °C. Extending our previous mathematical model, we have developed a computationally efficient model of cardiac cross-bridge mechanics which is sensitive to changes in cellular Pi. This extended model was parameterised with human cardiac complex modulus data. It captured the changes to cardiac mechanics following an increase in Pi concentration that we measured experimentally, including a reduced elastic modulus and a right-shift in frequency. The human cardiac trabecula we studied had a low sensitivity to Pi compared to what has been previously reported in mammalian cardiac tissue, which suggests that the muscle may have cellular compensatory mechanisms to cope with elevated Pi levels. This study demonstrates the feasibility of our experimental-modelling pipeline for future investigation of mechanical and metabolic effects in the diseased human heart.Clinical Relevance- This study presents the first measurement of the effect of Pi on the stiffness frequency response of human cardiac tissue and extends an experimental-modelling framework appropriate for investigating effects of disease on the human heart.

Calcium Overload and Mitochondrial Metabolism

Mitochondria calcium is a double-edged sword. While low levels of calcium are essential to maintain optimal rates of ATP production, extreme levels of calcium overcoming the mitochondrial calcium retention capacity leads to loss of mitochondrial function. In moderate amounts, however, ATP synthesis rates are inhibited in a calcium-titratable manner. While the consequences of extreme calcium overload are well-known, the effects on mitochondrial function in the moderately loaded range remain enigmatic. These observations are associated with changes in the mitochondria ultrastructure and cristae network. The present mini review/perspective follows up on previous studies using well-established cryo-electron microscopy and poses an explanation for the observable depressed ATP synthesis rates in mitochondria during calcium-overloaded states. The results presented herein suggest that the inhibition of oxidative phosphorylation is not caused by a direct decoupling of energy metabolism via the opening of a calcium-sensitive, proteinaceous pore but rather a separate but related calcium-dependent phenomenon. Such inhibition during calcium-overloaded states points towards mitochondrial ultrastructural modifications, enzyme activity changes, or an interplay between both events.

Effects of altered cellular ultrastructure on energy metabolism in diabetic cardiomyopathy: an in silico study

Diabetic cardiomyopathy is a leading cause of heart failure in diabetes. At the cellular level, diabetic cardiomyopathy leads to altered mitochondrial energy metabolism and cardiomyocyte ultrastructure. We combined electron microscopy (EM) and computational modelling to understand the impact of diabetes-induced ultrastructural changes on cardiac bioenergetics. We collected transverse micrographs of multiple control and type I diabetic rat cardiomyocytes using EM. Micrographs were converted to finite-element meshes, and bioenergetics was simulated over them using a biophysical model. The simulations also incorporated depressed mitochondrial capacity for oxidative phosphorylation (OXPHOS) and creatine kinase (CK) reactions to simulate diabetes-induced mitochondrial dysfunction. Analysis of micrographs revealed a 14% decline in mitochondrial area fraction in diabetic cardiomyocytes, and an irregular arrangement of mitochondria and myofibrils. Simulations predicted that this irregular arrangement, coupled with the depressed activity of mitochondrial CK enzymes, leads to large spatial variation in adenosine diphosphate (ADP)/adenosine triphosphate (ATP) ratio profile of diabetic cardiomyocytes. However, when spatially averaged, myofibrillar ADP/ATP ratios of a cardiomyocyte do not change with diabetes. Instead, average concentration of inorganic phosphate rises by 40% owing to lower mitochondrial area fraction and dysfunction in OXPHOS. These simulations indicate that a disorganized cellular ultrastructure negatively impacts metabolite transport in diabetic cardiomyopathy.

This article is part of the theme issue ‘The cardiomyocyte: new revelations on the interplay between architecture and function in growth, health, and disease’.

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.

Sensing local energetics to acutely regulate mitophagy in skeletal muscle

The energetic requirements of skeletal muscle to sustain movement, as during exercise, is met largely by mitochondria, which form an intricate, interconnected reticulum. Maintenance of a healthy mitochondrial reticulum is essential for skeletal muscle function, suggesting quality control pathways are spatially governed. Mitophagy, the process by which damaged and/or dysfunctional regions of the mitochondrial reticulum are removed and degraded, has emerged as an integral part of the molecular response to exercise. Upregulation of mitophagy in response to acute exercise is directly connected to energetic sensing mechanisms through AMPK. In this review, we discuss the connection of mitophagy to muscle energetics and how AMPK may spatially control mitophagy through multiple potential means.

The Cell Physiome: What Do We Need in a Computational Physiology Framework for Predicting Single-Cell Biology?

Modern biology and biomedicine are undergoing a big data explosion, needing advanced computational algorithms to extract mechanistic insights on the physiological state of living cells. We present the motivation for the Cell Physiome project: a framework and approach for creating, sharing, and using biophysics-based computational models of single-cell physiology. Using examples in calcium signaling, bioenergetics, and endosomal trafficking, we highlight the need for spatially detailed, biophysics-based computational models to uncover new mechanisms underlying cell biology. We review progress and challenges to date toward creating cell physiome models. We then introduce bond graphs as an efficient way to create cell physiome models that integrate chemical, mechanical, electromagnetic, and thermal processes while maintaining mass and energy balance. Bond graphs enhance modularization and reusability of computational models of cells at scale. We conclude with a look forward at steps that will help fully realize this exciting new field of mechanistic biomedical data science. Expected final online publication date for the Annual Review of Biomedical Data Science, Volume 5 is August 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

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.

Paradoxes of Hymenoptera flight muscles, extreme machines

In the Carboniferous, insects evolved flight. Intense selection drove for high performance and approximately 100 million years later, Hymenoptera (bees, wasps and ants) emerged. Some species had proportionately small wings, with apparently impossible aerodynamic challenges including a need for high frequency flight muscles (FMs), powered exclusively off aerobic pathways and resulting in extreme aerobic capacities. Modern insect FMs are the most refined and form large dense blocks that occupy 90% of the thorax. These can beat wings at 200 to 230 Hz, more than double that achieved by standard neuromuscular systems. To do so, rapid repolarisation was circumvented through evolution of asynchronous stimulation, stretch activation, elastic recoil and a paradoxically slow Ca2+ reuptake. While the latter conserves ATP, considerable ATP is demanded at the myofibrils. FMs have diminished sarcoplasmic volumes, and ATP is produced solely by mitochondria, which pack myocytes to maximal limits and have very dense cristae. Gaseous oxygen is supplied directly to mitochondria. While FMs appear to be optimised for function, several unusual paradoxes remain. FMs lack any significant equivalent to the creatine kinase shuttle, and myofibrils are twice as wide as those of within cardiomyocytes. The mitochondrial electron transport systems also release large amounts of reactive oxygen species (ROS) and respiratory complexes do not appear to be present at any exceptional level. Given that the loss of the creatine kinase shuttle and elevated ROS impairs heart function, we question how do FM shuttle adenylates at high rates and tolerate oxidative stress conditions that occur in diseased hearts?

Mitophagy-regulated mitochondrial health strongly protects the heart against cardiac dysfunction after acute myocardial infarction

Autophagy including mitophagy serves as an important regulatory mechanism in the heart to maintain the cellular homeostasis and to protect against heart damages caused by myocardial infarction (MI). The current study aims to dissect roles of general autophagy and specific mitophagy in regulating cardiac function after MI. By using Beclin1^+/- , Fundc1 knockout (KO) and Fundc1 transgenic (TG) mouse models, combined with starvation and MI models, we found that Fundc1 KO caused more severe mitochondrial and cardiac dysfunction damages than Beclin1^+/- after MI. Interestingly, Beclin1^+/- caused notable decrease of total autophagy without detectable change to mitophagy, and Fundc1 KO markedly suppressed mitophagy but did not change the total autophagy activity. In contrast, starvation increased total autophagy without changing mitophagy while Fundc1 TG elevated total autophagy and mitophagy in mouse hearts. As a result, Fundc1 TG provided much stronger protective effects than starvation after MI. Moreover, Beclin1^+/- /Fundc1 TG showed increased total autophagy and mitophagy to a level comparable to Fundc1 TG per se, and completely reversed Beclin1^+/- -caused aggravation of mitochondrial and cardiac injury after MI. Our results reveal that mitophagy but not general autophagy contributes predominantly to the cardiac protective effect through regulating mitochondrial function.

Isolation and reconstruction of cardiac mitochondria from SBEM images using a deep learning-based method

Mitochondrial morphological defects are a common feature of diseased cardiac myocytes. However, quantitative assessment of mitochondrial morphology is limited by the time-consuming manual segmentation of electron micrograph (EM) images. To advance understanding of the relation between morphological defects and dysfunction, an efficient morphological reconstruction method is desired to enable isolation and reconstruction of mitochondria from EM images. We propose a new method for isolating and reconstructing single mitochondria from serial block-face scanning EM (SBEM) images. CDeep3M, a cloud-based deep learning network for EM images, was used to segment mitochondrial interior volumes and boundaries. Post-processing was performed using both the predicted interior volume and exterior boundary to isolate and reconstruct individual mitochondria. Series of SBEM images from two separate cardiac myocytes were processed. The highest F1-score was 95% using 50 training datasets, greater than that for previously reported automated methods and comparable to manual segmentations. Accuracy of separation of individual mitochondria was 80% on a pixel basis. A total of 2315 mitochondria in the two series of SBEM images were evaluated with a mean volume of 0.78 µm³. The volume distribution was very broad and skewed; the most frequent mitochondria were 0.04-0.06 µm³, but mitochondria larger than 2.0 µm³ accounted for more than 10% of the total number. The average short-axis length was 0.47 µm. Primarily longitudinal mitochondria (0-30 degrees) were dominant (54%). This new automated segmentation and separation method can help quantitate mitochondrial morphology and improve understanding of myocyte structure-function relationships.

Mitochondria-localized AMPK responds to local energetics and contributes to exercise and energetic stress-induced mitophagy

Mitochondria form a complex, interconnected reticulum that is maintained through coordination among biogenesis, dynamic fission, and fusion and mitophagy, which are initiated in response to various cues to maintain energetic homeostasis. These cellular events, which make up mitochondrial quality control, act with remarkable spatial precision, but what governs such spatial specificity is poorly understood. Herein, we demonstrate that specific isoforms of the cellular bioenergetic sensor, 5' AMP-activated protein kinase (AMPKα1/α2/β2/γ1), are localized on the outer mitochondrial membrane, referred to as mitoAMPK, in various tissues in mice and humans. Activation of mitoAMPK varies across the reticulum in response to energetic stress, and inhibition of mitoAMPK activity attenuates exercise-induced mitophagy in skeletal muscle in vivo. Discovery of a mitochondrial pool of AMPK and its local importance for mitochondrial quality control underscores the complexity of sensing cellular energetics in vivo that has implications for targeting mitochondrial energetics for disease treatment.

Cell Volume (3D) Correlative Microscopy Facilitated by Intracellular Fluorescent Nanodiamonds as Multi-Modal Probes

Three-dimensional correlative light and electron microscopy (3D CLEM) is attaining popularity as a potential technique to explore the functional aspects of a cell together with high-resolution ultrastructural details across the cell volume. To perform such a 3D CLEM experiment, there is an imperative requirement for multi-modal probes that are both fluorescent and electron-dense. These multi-modal probes will serve as landmarks in matching up the large full cell volume datasets acquired by different imaging modalities. Fluorescent nanodiamonds (FNDs) are a unique nanosized, fluorescent, and electron-dense material from the nanocarbon family. We hereby propose a novel and straightforward method for executing 3D CLEM using FNDs as multi-modal landmarks. We demonstrate that FND is biocompatible and is easily identified both in living cell fluorescence imaging and in serial block-face scanning electron microscopy (SB-EM). We illustrate the method by registering multi-modal datasets.

Impaired Myocardial Energetics Causes Mechanical Dysfunction in Decompensated Failing Hearts

Cardiac mechanical function is supported by ATP hydrolysis, which provides the chemical-free energy to drive the molecular processes underlying cardiac pumping. Physiological rates of myocardial ATP consumption require the heart to resynthesize its entire ATP pool several times per minute. In the failing heart, cardiomyocyte metabolic dysfunction leads to a reduction in the capacity for ATP synthesis and associated free energy to drive cellular processes. Yet it remains unclear if and how metabolic/energetic dysfunction that occurs during heart failure affects mechanical function of the heart. We hypothesize that changes in phosphate metabolite concentrations (ATP, ADP, inorganic phosphate) that are associated with decompensation and failure have direct roles in impeding contractile function of the myocardium in heart failure, contributing to the whole-body phenotype. To test this hypothesis, a transverse aortic constriction (TAC) rat model of pressure overload, hypertrophy, and decompensation was used to assess relationships between metrics of whole-organ pump function and myocardial energetic state. A multiscale computational model of cardiac mechanoenergetic coupling was used to identify and quantify the contribution of metabolic dysfunction to observed mechanical dysfunction. Results show an overall reduction in capacity for oxidative ATP synthesis fueled by either fatty acid or carbohydrate substrates as well as a reduction in total levels of adenine nucleotides and creatine in myocardium from TAC animals compared to sham-operated controls. Changes in phosphate metabolite levels in the TAC rats are correlated with impaired mechanical function, consistent with the overall hypothesis. Furthermore, computational analysis of myocardial metabolism and contractile dynamics predicts that increased levels of inorganic phosphate in TAC compared to control animals kinetically impair the myosin ATPase crossbridge cycle in decompensated hypertrophy/heart failure.

Precision remodeling: how exercise improves mitochondrial quality in myofibers

The profound energetic demand of prolonged exercise imposed upon skeletal muscle and the heart is met by oxidation of substrate within mitochondria. As such, several coordinated events are initiated in order to maintain mitochondria, collectively known as mitochondrial quality control. In this review, we discuss how mitochondrial quality control functions to maintain the integrity of the reticulum and energy production in response to prolonged exercise, as well as the relevant signaling events that dictate these responses. Based upon the prevailing data in the field, we propose a model where exercise-mediated quality control may be chiefly regulated through local mechanisms, thus allowing for the remarkable precision in mitochondrial quality control events.

Quantifying drug-induced structural toxicity in hepatocytes and cardiomyocytes derived from hiPSCs using a deep learning method

Cardiac and hepatic toxicity result from induced disruption of the functioning of cardiomyocytes and hepatocytes, respectively, which is tightly related to the organization of their subcellular structures. Cellular structure can be analyzed from microscopy imaging data. However, subtle or complex structural changes that are not easily perceived may be missed by conventional image-analysis techniques. Here we report the evaluation of PhenoTox, an image-based deep-learning method of quantifying drug-induced structural changes using human hepatocytes and cardiomyocytes derived from human induced pluripotent stem cells. We assessed the ability of the deep learning method to detect variations in the organization of cellular structures from images of fixed or live cells. We also evaluated the power and sensitivity of the method for detecting toxic effects of drugs by conducting a set of experiments using known toxicants and other methods of screening for cytotoxic effects. Moreover, we used PhenoTox to characterize the effects of tamoxifen and doxorubicin-which cause liver toxicity-on hepatocytes. PhenoTox revealed differences related to loss of cytochrome P450 3A4 activity, for which it showed greater sensitivity than a caspase 3/7 assay. Finally, PhenoTox detected structural toxicity in cardiomyocytes, which was correlated with contractility defects induced by doxorubicin, erlotinib, and sorafenib. Taken together, the results demonstrated that PhenoTox can capture the subtle morphological changes that are early signs of toxicity in both hepatocytes and cardiomyocytes.

Impact of global structure on diffusive exploration of organelle networks

We investigate diffusive search on planar networks, motivated by tubular organelle networks in cell biology that contain molecules searching for reaction partners and binding sites. Exact calculation of the diffusive mean first-passage time on a spatial network is used to characterize the typical search time as a function of network connectivity. We find that global structural properties — the total edge length and number of loops — are sufficient to largely determine network exploration times for a variety of both synthetic planar networks and organelle morphologies extracted from living cells. For synthetic networks on a lattice, we predict the search time dependence on these global structural parameters by connecting with percolation theory, providing a bridge from irregular real-world networks to a simpler physical model. The dependence of search time on global network structural properties suggests that network architecture can be designed for efficient search without controlling the precise arrangement of connections. Specifically, increasing the number of loops substantially decreases search times, pointing to a potential physical mechanism for regulating reaction rates within organelle network structures.