Mitochondrial energetics, pH regulation, and ion dynamics: a computational-experimental approach

Mitofusin 2-containing mitochondrial-reticular microdomains direct rapid cardiomyocyte bioenergetic responses via interorganelle Ca(2+) crosstalk

RATIONALE

Mitochondrial Ca(2+) uptake is essential for the bioenergetic feedback response through stimulation of Krebs cycle dehydrogenases. Close association of mitochondria to the sarcoplasmic reticulum (SR) may explain efficient mitochondrial Ca(2+) uptake despite low Ca(2+) affinity of the mitochondrial Ca(2+) uniporter. However, the existence of such mitochondrial Ca(2+) microdomains and their functional role are presently unresolved. Mitofusin (Mfn) 1 and 2 mediate mitochondrial outer membrane fusion, whereas Mfn2 but not Mfn1 tethers endoplasmic reticulum to mitochondria in noncardiac cells.

OBJECTIVE

To elucidate roles for Mfn1 and 2 in SR-mitochondrial tethering, Ca(2+) signaling, and bioenergetic regulation in cardiac myocytes.

METHODS AND RESULTS

Fruit fly heart tubes deficient of the Drosophila Mfn ortholog MARF had increased contraction-associated and caffeine-sensitive Ca(2+) release, suggesting a role for Mfn in SR Ca(2+) handling. Whereas cardiac-specific Mfn1 ablation had no effects on murine heart function or Ca(2+) cycling, Mfn2 deficiency decreased cardiomyocyte SR-mitochondrial contact length by 30% and reduced the content of SR-associated proteins in mitochondria-associated membranes. This was associated with decreased mitochondrial Ca(2+) uptake (despite unchanged mitochondrial membrane potential) but increased steady-state and caffeine-induced SR Ca(2+) release. Accordingly, Ca(2+)-induced stimulation of Krebs cycle dehydrogenases during β-adrenergic stimulation was hampered in Mfn2-KO but not Mfn1-KO myocytes, evidenced by oxidation of the redox states of NAD(P)H/NAD(P)(+) and FADH(2)/FAD.

CONCLUSIONS

Physical tethering of SR and mitochondria via Mfn2 is essential for normal interorganelle Ca(2+) signaling in the myocardium, consistent with a requirement for SR-mitochondrial Ca(2+) signaling through microdomains in the cardiomyocyte bioenergetic feedback response to physiological stress.

Toward an integrative computational model of the Guinea pig cardiac myocyte

The local control theory of excitation-contraction (EC) coupling asserts that regulation of calcium (Ca²⁺) release occurs at the nanodomain level, where openings of single L-type Ca²⁺ channels (LCCs) trigger openings of small clusters of ryanodine receptors (RyRs) co-localized within the dyad. A consequence of local control is that the whole-cell Ca²⁺ transient is a smooth continuous function of influx of Ca²⁺ through LCCs. While this so-called graded release property has been known for some time, its functional importance to the integrated behavior of the cardiac ventricular myocyte has not been fully appreciated. We previously formulated a biophysically based model, in which LCCs and RyRs interact via a coarse-grained representation of the dyadic space. The model captures key features of local control using a low-dimensional system of ordinary differential equations. Voltage-dependent gain and graded Ca²⁺ release are emergent properties of this model by virtue of the fact that model formulation is closely based on the sub-cellular basis of local control. In this current work, we have incorporated this graded release model into a prior model of guinea pig ventricular myocyte electrophysiology, metabolism, and isometric force production. The resulting integrative model predicts the experimentally observed causal relationship between action potential (AP) shape and timing of Ca²⁺ and force transients, a relationship that is not explained by models lacking the graded release property. Model results suggest that even relatively subtle changes in AP morphology that may result, for example, from remodeling of membrane transporter expression in disease or spatial variation in cell properties, may have major impact on the temporal waveform of Ca²⁺ transients, thus influencing tissue level electromechanical function.

Computational modeling of mitochondrial function

The advent of techniques with the ability to scan massive changes in cellular makeup (genomics, proteomics, etc.) has revealed the compelling need for analytical methods to interpret and make sense of those changes. Computational models built on sound physico-chemical mechanistic basis are unavoidable at the time of integrating, interpreting, and simulating high-throughput experimental data. Another powerful role of computational models is predicting new behavior provided they are adequately validated.Mitochondrial energy transduction has been traditionally studied with thermodynamic models. More recently, kinetic or thermo-kinetic models have been proposed, leading the path toward an understanding of the control and regulation of mitochondrial energy metabolism and its interaction with cytoplasmic and other compartments. In this work, we outline the methods, step-by-step, that should be followed to build a computational model of mitochondrial energetics in isolation or integrated to a network of cellular processes. Depending on the question addressed by the modeler, the methodology explained herein can be applied with different levels of detail, from the mitochondrial energy producing machinery in a network of cellular processes to the dynamics of a single enzyme during its catalytic cycle.