Early studies have shown that moderate levels of calcium overload can cause lower oxidative phosphorylation rates. However, the mechanistic interpretations of these findings were inadequate. And while the effect of excessive calcium overload on mitochondrial function is well appreciated, there has been little to no reports on the consequences of low to moderate calcium overload. To resolve this inadequacy, mitochondrial function from guinea pig hearts was quantified using several well-established methods including high-resolution respirometry and spectrofluorimetry and analyzed using mathematical modeling. We measured key mitochondrial variables such as respiration, mitochondrial membrane potential, buffer calcium, and substrate effects for a range of mitochondrial calcium loads from near zero to levels approaching mitochondrial permeability transition. In addition, we developed a computer model closely mimicking the experimental conditions and used this model to design experiments capable of eliminating many hypotheses generated from the data analysis. We subsequently performed those experiments and determined why mitochondrial ADP-stimulated respiration is significantly lowered during calcium overload. We found that when calcium phosphate levels, not matrix free calcium, reached sufficient levels, complex I activity is inhibited, and the rate of ATP synthesis is reduced. Our findings suggest that calcium phosphate granules form physical barriers that isolate complex I from NADH, disrupt complex I activity, or destabilize cristae and inhibit NADH-dependent respiration.

Mitochondria have a remarkable ability to uptake and store massive amounts of calcium. However, the consequences of massive calcium accumulation remain enigmatic. In the present study, we analyzed a series of time-course experiments to identify the sequence of events that occur in a population of guinea pig cardiac mitochondria exposed to excessive calcium overload that cause mitochondrial permeability transition (MPT). By analyzing coincident structural and functional data, we determined that excessive calcium overload is associated with large calcium phosphate granules and inner membrane fragmentation, which explains the extent of mitochondrial dysfunction. This data also reveals a novel mechanism for cyclosporin A, an inhibitor of MPT, in which it preserves cristae despite the presence of massive calcium phosphate granules in the matrix. Overall, these findings establish a mechanism of calcium-induced mitochondrial dysfunction and the impact of calcium regulation on mitochondrial structure and function.

Mitochondrial permeability transition (PT) is a phenomenon of stress-induced increase in nonspecific permeability of the mitochondrial inner membrane that leads to disruption of oxidative phosphorylation and cell death. Quantitative measurement of the membrane permeability increase during PT is critically important for understanding the PT’s impact on mitochondrial function. The elementary unit of PT is a PT pore (PTP), a single channel presumably formed by either ATP synthase or adenine nucleotide translocator (ANT). It is not known how many channels are open in a single mitochondrion during PT, which makes it difficult to quantitatively estimate the overall degree of membrane permeability. Here, we used wide-field microscopy to record mitochondrial swelling and quantitatively measure rates of single-mitochondrion volume increase during PT-induced high-amplitude swelling. PT was quantified by calculating the rates of water flux responsible for measured volume changes. The total water flux through the mitochondrial membrane of a single mitochondrion during PT was in the range of (2.5 ± 0.4) × 10−17 kg/s for swelling in 2 mM Ca2+ and (1.1 ± 0.2) × 10−17 kg/s for swelling in 200 µM Ca2+. Under these experimental conditions, a single PTP channel with ionic conductance of 1.5 nS could allow passage of water at the rate of 0.65 × 10−17 kg/s. Thus, we estimate the integral ionic conductance of the whole mitochondrion during PT to be 5.9 ± 0.9 nS for 2 mM concentration of Ca2+ and 2.6 ± 0.4 nS for 200 µM of Ca2+. The number of PTPs per mitochondrion ranged from one to nine. Due to the uncertainties in PTP structure and model parameters, PTP count results may be slightly underestimated. However, taking into account that each mitochondrion has ∼15,000 copies of ATP synthases and ANTs, our data imply that PTP activation is a rare event that occurs only in a small subpopulation of these proteins.

Mitochondrial cristae are dynamic invaginations of the inner membrane and play a key role in its metabolic capacity to produce ATP. Structural alterations caused by either genetic abnormalities or detrimental environmental factors impede mitochondrial metabolic fluxes and lead to a decrease in their ability to meet metabolic energy requirements. While some of the key proteins associated with mitochondrial cristae are known, very little is known about how the inner membrane dynamics are involved in energy metabolism. In this study, we present a computational strategy to understand how cristae are formed using a phase-based separation approach of both the inner membrane space and matrix space, which are explicitly modeled using the Cahn–Hilliard equation. We show that cristae are formed as a consequence of minimizing an energy function associated with phase interactions which are subject to geometric boundary constraints. We then extended the model to explore how the presence of calcium phosphate granules, entities that form in calcium overload conditions, exert a devastating inner membrane remodeling response that reduces the capacity for mitochondria to produce ATP. This modeling approach can be extended to include arbitrary geometrical constraints, the spatial heterogeneity of enzymes, and electrostatic effects to mechanize the impact of ultrastructural changes on energy metabolism.

Mitochondrial reactive oxygen species (ROS) play important roles in cellular signaling; however, certain pathological conditions such as ischemia/reperfusion (I/R) injury disrupt ROS homeostasis and contribute to cell death. A major impediment to developing therapeutic measures against oxidative stress induced cellular damage is the lack of a quantitative framework to identify the specific sources and regulatory mechanisms of mitochondrial ROS production. We developed a thermodynamically consistent, mass-and-charge balanced, kinetic model of mitochondrial ROS homeostasis focused on redox sites of electron transport chain complexes I, II, and III. The model was calibrated and corroborated using comprehensive data sets relevant to ROS homeostasis. The model predicts that complex I ROS production dominates other sources under conditions favoring a high membrane potential with elevated NADH and QH2 levels. In general, complex I contributes to significant levels of ROS production under pathological conditions, while complexes II and III are responsible for basal levels of ROS production, especially when QH2 levels are elevated. The model also reveals that hydrogen peroxide production by complex I underlies the non-linear relationship between ROS emission and O2 at low O2 concentrations. Lastly, the model highlights the need to quantify scavenging system activity under different conditions to establish a complete picture of mitochondrial ROS homeostasis. In summary, we describe the individual contributions of the ETS complex redox sites to total ROS emission in mitochondria respiring under various combinations of NADH- and Q-linked respiratory fuels under varying workloads.