Analysis of the kinetics and bistability of ubiquinol:cytochrome c oxidoreductase

MtDNA Copy Number in Oral Epithelial Cells Serves as a Potential Biomarker of Mitochondrial Damage by Neonicotinoid Exposure: A Cross-Sectional Study

As the mitochondrial DNA copy number (mtDNAcn) has been reported to be a biomarker for mtDNA damage in honeybees when exposed to sublethal neonicotinoids, the feasibility of using human mitochondria as a predictor upon neonicotinoid exposure remains elusive. This study investigated the association between the urinary neonicotinoid and the relative mtDNAcn (RmtDNAcn) of oral epithelial cells collected in a cross-sectional study with repeated measurements over 6 weeks. The molecular mechanism underlying neonicotinoid-caused mitochondrial damage was also examined by in vitro assay. Herein, the average integrated urinary neonicotinoid (IMIRPF) concentration ranged from 8.01 to 13.70 μg/L (specific gravity-adjusted) during the sampling period. Concomitantly, with an increase in the urinary IMIRPF, the RmtDNAcn significantly increased from 1.20 (low group) to 1.93 (high group), indicating potential dose-dependent mitochondrial damage. Furthermore, the linear regression analysis confirmed the significant correlation between the IMIRPF and RmtDNAcn. Results from in vitro assays demonstrated that neonicotinoid exposure led to the inhibition of the genes encoding mitochondrial oxidative phosphorylation (OXPHOS) complexes I and III (e.g., ND2, ND6, CytB, and CYC1), accompanied by increased reactive oxygen species production in SH-SY5Y cells. Conjointly, neonicotinoid exposure led to mitochondrial dysfunction and a resulting increase in the RmtDNAcn, which may serve as a plausible biomarker in humans.

Effects of ROS pathway inhibitors and NADH and FADH2 linked substrates on mitochondrial bioenergetics and ROS emission in the heart and kidney cortex and outer medulla

Mitochondria are major sources of reactive oxygen species (ROS), which play important roles in both physiological and pathological processes. However, the specific contributions of different ROS production and scavenging components in the mitochondria of metabolically active tissues such as heart and kidney cortex and outer medulla (OM) are not well understood. Therefore, the goal of this study was to determine contributions of different ROS production and scavenging components and provide detailed comparisons of mitochondrial respiration, bioenergetics, ROS emission between the heart and kidney cortex and OM using tissues obtained from the same Sprague-Dawley rat under identical conditions and perturbations. Specifically, data were obtained using both NADH-linked substrate pyruvate + malate and FADH2-linked substrate succinate followed by additions of inhibitors of different components of the electron transport chain (ETC) and oxidative phosphorylation (OxPhos) and other ROS production and scavenging systems. Currently, there is limited data available for the mitochondria of kidney cortex and OM, the two major energy-consuming tissues in the body only next to the heart, and scarce quantitative information on the interplay between mitochondrial ROS production and scavenging systems in the three tissues. The findings from this study demonstrate significant differences in mitochondrial respiratory and bioenergetic functions and ROS emission among the three tissues. The results quantify the rates of ROS production from different complexes of the ETC, identify the complexes responsible for variations in mitochondrial membrane depolarization and regulations of ROS production, and quantify the contributions of ROS scavenging enzymes towards overall mitochondrial ROS emission. These findings advance our fundamental knowledge of tissue-specific and substrate-dependent mitochondrial respiratory and bioenergetic functions and ROS emission. This is important given the critical role that excess ROS production, oxidative stress, and mitochondrial dysfunction in the heart and kidney cortex and OM play in the pathogenesis of cardiovascular and renal diseases, including salt-sensitive hypertension.

AW, Dash RK. Computational Modeling of Substrate-Dependent Mitochondrial Respiration and Bioenergetics in the Heart and Kidney Cortex and Outer Medulla

Integrated computational modeling provides a mechanistic and quantitative framework to characterize alterations in mitochondrial respiration and bioenergetics in response to different metabolic substrates in-silico. These alterations play critical roles in the pathogenesis of diseases affecting metabolically active organs such as heart and kidney. Therefore, the present study aimed to develop and validate thermodynamically constrained integrated computational models of mitochondrial respiration and bioenergetics in the heart and kidney cortex and outer medulla (OM). The models incorporated the kinetics of major biochemical reactions and transport processes as well as regulatory mechanisms in the mitochondria of these tissues. Intrinsic model parameters such as Michaelis-Menten constants were fixed at previously estimated values, while extrinsic model parameters such as maximal reaction and transport velocities were estimated separately for each tissue. This was achieved by fitting the model solutions to our recently published respirometry data measured in isolated rat heart and kidney cortex and OM mitochondria utilizing various NADH- and FADH2-linked metabolic substrates. The models were validated by predicting additional respirometry and bioenergetics data, which were not used for estimating the extrinsic model parameters. The models were able to predict tissue-specific and substrate-dependent mitochondrial emergent metabolic system properties such as redox states, enzyme and transporter fluxes, metabolite concentrations, membrane potential, and respiratory control index under diverse physiological and pathological conditions. The models were also able to quantitatively characterize differential regulations of NADH- and FADH2-linked metabolic pathways, which contribute differently toward regulations of oxidative phosphorylation and ATP synthesis in the heart and kidney cortex and OM mitochondria.

The evolution of the human mitochondrial bc1 complex- adaptation for reduced rate of superoxide production?

The mitochondrial bc1 complex is a major source of mitochondrial superoxide. While bc1-generated superoxide plays a beneficial signaling role, excess production of superoxide lead to aging and degenerative diseases. The catalytic core of bc1 comprises three peptides -cytochrome b, Fe-S protein, and cytochrome c1. All three core peptides exhibit accelerated evolution in anthropoid primates. It has been suggested that the evolution of cytochrome b in anthropoids was driven by a pressure to reduce the production of superoxide. In humans, the bc1 core peptides exhibit anthropoid-specific substitutions that are clustered near functionally critical sites that may affect the production of superoxide. Here we compare the high-resolution structures of bovine, mouse, sheep and human bc1 to identify structural changes that are associated with human-specific substitutions. Several cytochrome b substitutions in humans alter its interactions with other subunits. Most significantly, there is a cluster of seven substitutions, in cytochrome b, the Fe-S protein, and cytochrome c1 that affect the interactions between these proteins at the tether arm of the Fe-S protein and may alter the rate of ubiquinone oxidation and the rate of superoxide production. Another cluster of substitutions near heme bH and the ubiquinone reduction site, Qi, may affect the rate of ubiquinone reduction and thus alter the rate of superoxide production. These results are compatible with the hypothesis that cytochrome b in humans (and other anthropoid primates) evolve to reduce the rate of production of superoxide thus enabling the exceptional longevity and exceptional cognitive ability of humans.

Kinetic Mathematical Modeling of Oxidative Phosphorylation in Cardiomyocyte Mitochondria

Oxidative phosphorylation (OXPHOS) is an oxygen-dependent process that consumes catabolized nutrients to produce adenosine triphosphate (ATP) to drive energy-dependent biological processes such as excitation-contraction coupling in cardiomyocytes. In addition to in vivo and in vitro experiments, in silico models are valuable for investigating the underlying mechanisms of OXPHOS and predicting its consequences in both physiological and pathological conditions. Here, we compare several prominent kinetic models of OXPHOS in cardiomyocytes. We examine how their mathematical expressions were derived, how their parameters were obtained, the conditions of their experimental counterparts, and the predictions they generated. We aim to explore the general landscape of energy production mechanisms in cardiomyocytes for future in silico models.

Mechanisms and mathematical modelling of ROS production by the mitochondrial electron transport chain

Reactive oxygen species (ROS) are recognised both as damaging molecules and intracellular signalling entities. In addition to its role in ATP generation, the mitochondrial electron transport chain (ETC) constitutes a relevant source of mitochondrial ROS, in particular during pathological conditions. Mitochondrial ROS homeostasis depends on species- and site-dependent ROS production, their bioreactivity, diffusion, and scavenging. However, our quantitative understanding of mitochondrial ROS homeostasis has thus far been hampered by technical limitations, including lack of truly site- and/or ROS-specific reporter molecules. In this context, the use of computational models is of great value to complement and interpret empirical data, as well as to predict variables that are difficult to assess experimentally. During the last decades, various mechanistic models of ETC-mediated ROS production have been developed. Although these often-complex models have generated novel insights, their parameterisation, analysis, and integration with other computational models is not straightforward. In contrast, phenomenological (sometimes termed "minimal") models use a relatively small set of equations to describe empirical relationship(s) between ROS-related and other parameters, and generally aim to explore system behaviour and generate hypotheses for experimental validation. In this review, we first discuss ETC-linked ROS homeostasis and introduce various detailed mechanistic models. Next, we present how bioenergetic parameters (e.g. NADH/NAD⁺ ratio, mitochondrial membrane potential) relate to site-specific ROS production within the ETC and how these relationships can be used to design minimal models of ROS homeostasis. Finally, we illustrate how minimal models have been applied to explore pathophysiological aspects of ROS.

What is the Role of Lipid Membrane-embedded Quinones in Mitochondria and Chloroplasts? Chemiosmotic Q-cycle versus Murburn Reaction Perspective

Quinones are found in the lipid membranes of prokaryotes like E. coli and cyanobacteria, and are also abundant in eukaryotic mitochondria and chloroplasts. They are intricately involved in the reaction mechanism of redox phosphorylations. In the Mitchellian chemiosmotic school of thought, membrane-lodged quinones are perceived as highly mobile conveyors of two-electron equivalents from the first leg of Electron Transport Chain (ETC) to the 'second pit-stop' of Cytochrome bc₁ or b₆f complex (CBC), where they undergo a regenerative 'Q-cycle'. In Manoj's murburn mechanism, the membrane-lodged quinones are perceived as relatively slow-moving one- or two- electron donors/acceptors, enabling charge separation and the CBC resets a one-electron paradigm via 'turbo logic'. Herein, we compare various purviews of the two mechanistic schools with respect to: constraints in mobility, protons' availability, binding of quinones with proteins, structural features of the protein complexes, energetics of reaction, overall reaction logic, etc. From various perspectives, the murburn mechanism appeals as a viable alternative explanation well-rooted in thermodynamics/kinetics and one which lends adequate structure-function correlations for the roles of quinones, lipid membrane and associated proteins.

Hysteresis and bistability in the succinate-CoQ reductase activity and reactive oxygen species production in the mitochondrial respiratory complex II

The mitochondrial respiratory Complex II (CII) is one of key enzymes of cell energy metabolism, linking the tricarboxylic acid (TCA) cycle and the electron transport chain (ETC). CII reversibly oxidizes succinate to fumarate in the TCA cycle and transfers the electrons, produced by this reaction to the membrane quinone pool, providing ubiquinol QH₂ to ETC. CII is also known as a generator of reactive oxygen species (ROS).

It was shown experimentally that succinate can serve as not only a substrate in the forward succinate-quinone oxidoreductase (SQR) direction, but also an enzyme activator. Molecular and kinetic mechanisms of this property of CII are still unclear.

In order to account for activation of CII by succinate in the forward SQR direction, we developed and analyzed a computational mechanistic model of electron transfer and ROS formation in CII. It was found that re-binding of succinate to the unoccupied dicarboxylate binding site when FAD is reduced with subsequent oxidation of FADH₂ creates a positive feedback loop in the succinate oxidation. The model predicts that this positive feedback can result in hysteresis and bistable switches in SQR activity and ROS production in CII. This requires that the rate constant of re-binding of succinate has to be higher than the rate constant of the initial succinate binding to the active center when FAD is oxidized.

Hysteresis and bistability in the SQR activity and ROS production in CII can play an important physiological role. In the presence of hysteresis with two stable branches with high and low SQR activity, high SQR activity is maintained even with a very strong drop in the succinate concentration, which may be necessary in the process of cell functioning in stressful situations. For the same reason, a high stationary rate of ROS production in CII can be maintained at low succinate concentrations.

•

Computational model of electron flows in the respiratory Complex II was developed.

•

The Complex II model predicts a positive feedback loop in the succinate oxidation.

•

Complex II can operate as a bistable switch between two alternative stable states.

•

ROS production in Complex II may have a hysteretic behaviour.

Modelling mitochondrial ROS production by the respiratory chain

ROS (superoxide and oxygen peroxide in this paper) play a dual role as signalling molecules and strong oxidizing agents leading to oxidative stress. Their production mainly occurs in mitochondria although they may have other locations (such as NADPH oxidase in particular cell types). Mitochondrial ROS production depends in an interweaving way upon many factors such as the membrane potential, the cell type and the respiratory substrates. Moreover, it is experimentally difficult to quantitatively assess the contribution of each potential site in the respiratory chain. To overcome these difficulties, mathematical models have been developed with different degrees of complexity in order to analyse different physiological questions ranging from a simple reproduction/simulation of experimental results to a detailed model of the possible mechanisms leading to ROS production. Here, we analyse experimental results concerning ROS production including results still under discussion. We then critically review the three models of ROS production in the whole respiratory chain available in the literature and propose some direction for future modelling work.

Atoms to Phenotypes: Molecular Design Principles of Cellular Energy Metabolism

We report a 100-million atom-scale model of an entire cell organelle, a photosynthetic chromatophore vesicle from a purple bacterium, that reveals the cascade of energy conversion steps culminating in the generation of ATP from sunlight. Molecular dynamics simulations of this vesicle elucidate how the integral membrane complexes influence local curvature to tune photoexcitation of pigments. Brownian dynamics of small molecules within the chromatophore probe the mechanisms of directional charge transport under various pH and salinity conditions. Reproducing phenotypic properties from atomistic details, a kinetic model evinces that low-light adaptations of the bacterium emerge as a spontaneous outcome of optimizing the balance between the chromatophore's structural integrity and robust energy conversion. Parallels are drawn with the more universal mitochondrial bioenergetic machinery, from whence molecular-scale insights into the mechanism of cellular aging are inferred. Together, our integrative method and spectroscopic experiments pave the way to first-principles modeling of whole living cells.

Analysis of a Functional Dimer Model of Ubiquinol Cytochrome c Oxidoreductase

Ubiquinol cytochrome c oxidoreductase (bc₁ complex) serves as an important electron junction in many respiratory systems. It funnels electrons coming from NADH and ubiquinol to cytochrome c, but it is also capable of producing significant amounts of the free radical superoxide. In situ and in other experimental systems, the enzyme exists as a dimer. But until recently, it was believed to operate as a functional monomer. Here we show that a functional dimer model is capable of explaining both kinetic and superoxide production rate data. The model consists of six electronic states characterized by the number of electrons deposited on the complex. It is fully reversible and strictly adheres to the thermodynamics governing the reactions. A total of nine independent data sets were used to parameterize the model. To explain the data with a consistent set of parameters, it was necessary to incorporate intramonomer Coulombic effects between hemes b_L and b_H and intermonomer Coulombic effects between b_L hemes. The fitted repulsion energies fall within the theoretical range of electrostatic calculations. In addition, model analysis demonstrates that the Q pool is mostly oxidized under normal physiological operation but can switch to a more reduced state when reverse electron transport conditions are in place.

Feedback Regulation and Time Hierarchy of Oxidative Phosphorylation in Cardiac Mitochondria

To determine how oxidative ATP synthesis is regulated in the heart, the responses of cardiac mitochondria oxidizing pyruvate to alterations in [ATP], [ADP], and inorganic phosphate ([Pi]) were characterized over a range of steady-state levels of extramitochondrial [ATP], [ADP], and [Pi]. Evolution of the steady states of the measured variables with the flux of respiration shows that: (1) a higher phosphorylation potential is achieved by mitochondria at higher [Pi] for a given flux of respiration; (2) the time hierarchy of oxidative phosphorylation is given by phosphorylation subsystem, electron transport chain, and substrate dehydrogenation subsystems listed in increasing order of their response times; (3) the matrix ATP hydrolysis mass action ratio [ADP] × [Pi]/[ATP] provides feedback to the substrate dehydrogenation flux over the entire range of respiratory flux examined in this study; and finally, (4) contrary to previous models of regulation of oxidative phosphorylation, [Pi] does not modulate the activity of complex III.

Catalytic Coupling of Oxidative Phosphorylation, ATP Demand, and Reactive Oxygen Species Generation

Competing models of mitochondrial energy metabolism in the heart are highly disputed. In addition, the mechanisms of reactive oxygen species (ROS) production and scavenging are not well understood. To deepen our understanding of these processes, a computer model was developed to integrate the biophysical processes of oxidative phosphorylation and ROS generation. The model was calibrated with experimental data obtained from isolated rat heart mitochondria subjected to physiological conditions and workloads. Model simulations show that changes in the quinone pool redox state are responsible for the apparent inorganic phosphate activation of complex III. Model simulations predict that complex III is responsible for more ROS production during physiological working conditions relative to complex I. However, this relationship is reversed under pathological conditions. Finally, model analysis reveals how a highly reduced quinone pool caused by elevated levels of succinate is likely responsible for the burst of ROS seen during reperfusion after ischemia.

Global Kinetic Analysis of Mammalian E3 Reveals pH-dependent NAD+/NADH Regulation, Physiological Kinetic Reversibility, and Catalytic Optimum

Mammalian E3 is an essential mitochondrial enzyme responsible for catalyzing the terminal reaction in the oxidative catabolism of several metabolites. E3 is a key regulator of metabolic fuel selection as a component of the pyruvate dehydrogenase complex (PDHc). E3 regulates PDHc activity by altering the affinity of pyruvate dehydrogenase kinase, an inhibitor of the enzyme complex, through changes in reduction and acetylation state of lipoamide moieties set by the NAD(+)/NADH ratio. Thus, an accurate kinetic model of E3 is needed to predict overall mammalian PDHc activity. Here, we have combined numerous literature data sets and new equilibrium spectroscopic experiments with a multitude of independently collected forward and reverse steady-state kinetic assays using pig heart E3. The latter kinetic assays demonstrate a pH-dependent transition of NAD(+) activation to inhibition, shown here, to our knowledge, for the first time in a single consistent data set. Experimental data were analyzed to yield a thermodynamically constrained four-redox-state model of E3 that simulates pH-dependent activation/inhibition and active site redox states for various conditions. The developed model was used to determine substrate/product conditions that give maximal E3 rates and show that, due to non-Michaelis-Menten behavior, the maximal flux is different compared with the classically defined kcat.

Determination of the catalytic mechanism for mitochondrial malate dehydrogenase

The kinetics of malate dehydrogenase (MDH) catalyzed oxidation/reduction of L-malate/oxaloacetate is pH-dependent due to the proton generated/taken up during the reaction. Previous kinetic studies on the mitochondrial MDH did not yield a consensus kinetic model that explains both substrate and pH dependency of the initial velocity. In this study, we propose, to our knowledge, a new kinetic mechanism to explain kinetic data acquired over a range of pH and substrate concentrations. Progress curves in the forward and reverse reaction directions were obtained under a variety of reactant concentrations to identify associated kinetic parameters. Experiments were conducted at physiologically relevant ionic strength of 0.17 M, pH ranging between 6.5 and 9.0, and at 25 °C. The developed model was built on the prior observation of proton uptake upon binding of NADH to MDH, and that the MDH-catalyzed oxidation of NADH may follow an ordered bi-bi mechanism with NADH/NAD binding to the enzyme first, followed by the binding of oxaloacetate/L-malate. This basic mechanism was expanded to account for additional ionic states to explain the pH dependency of the kinetic behavior, resulting in what we believe to be the first kinetic model explaining both substrate and pH dependency of the reaction velocity.

Computational modeling analysis of mitochondrial superoxide production under varying substrate conditions and upon inhibition of different segments of the electron transport chain

A computational mechanistic model of superoxide (O2•-) formation in the mitochondrial electron transport chain (ETC) was developed to facilitate the quantitative analysis of factors controlling mitochondrial O2•- production and assist in the interpretation of experimental studies. The model takes into account all individual electron transfer reactions in Complexes I and III. The model accounts for multiple, often seemingly contradictory observations on the effects of ΔΨ and ΔpH, and for the effects of multiple substrate and inhibitor conditions, including differential effects of Complex III inhibitors antimycin A, myxothiazol and stigmatellin. Simulation results confirm that, in addition to O2•- formation in Complex III and at the flavin site of Complex I, the quinone binding site of Complex I is an additional superoxide generating site that accounts for experimental observations on O2•- production during reverse electron transfer. However, our simulation results predict that, when cytochrome c oxidase is inhibited during oxidation of succinate, ROS production at this site is eliminated and almost all superoxide in Complex I is generated by reduced FMN, even when the redox pressure for reverse electron transfer from succinate is strong. In addition, the model indicates that conflicting literature data on the kinetics of electron transfer in Complex III involving the iron-sulfur protein-cytochrome bL complex can be resolved in favor of a dissociation of the protein only after electron transfer to cytochrome bH. The model predictions can be helpful in understanding factors driving mitochondrial superoxide formation in intact cells and tissues.

A pH-dependent kinetic model of dihydrolipoamide dehydrogenase from multiple organisms

Dihydrolipoamide dehydrogenase is a flavoenzyme that reversibly catalyzes the oxidation of reduced lipoyl substrates with the reduction of NAD(+) to NADH. In vivo, the dihydrolipoamide dehydrogenase component (E3) is associated with the pyruvate, α-ketoglutarate, and glycine dehydrogenase complexes. The pyruvate dehydrogenase (PDH) complex connects the glycolytic flux to the tricarboxylic acid cycle and is central to the regulation of primary metabolism. Regulation of PDH via regulation of the E3 component by the NAD(+)/NADH ratio represents one of the important physiological control mechanisms of PDH activity. Furthermore, previous experiments with the isolated E3 component have demonstrated the importance of pH in dictating NAD(+)/NADH ratio effects on enzymatic activity. Here, we show that a three-state mechanism that represents the major redox states of the enzyme and includes a detailed representation of the active-site chemistry constrained by both equilibrium and thermodynamic loop constraints can be used to model regulatory NAD(+)/NADH ratio and pH effects demonstrated in progress-curve and initial-velocity data sets from rat, human, Escherichia coli, and spinach enzymes. Global fitting of the model provides stable predictions to the steady-state distributions of enzyme redox states as a function of lipoamide/dihydrolipoamide, NAD(+)/NADH, and pH. These distributions were calculated using physiological NAD(+)/NADH ratios representative of the diverse organismal sources of E3 analyzed in this study. This mechanistically detailed, thermodynamically constrained, pH-dependent model of E3 provides a stable platform on which to accurately model multicomponent enzyme complexes that implement E3 from a variety of organisms.

Determining the origins of superoxide and hydrogen peroxide in the mammalian NADH:ubiquinone oxidoreductase

NADH:ubiquinone oxidoreductase (complex I) is a proton pump in the electron transport chain that can produce a significant amounts of superoxide and hydrogen peroxide. While the flavin mononucleotide (FMN) is the putative site for hydrogen peroxide generation, sites responsible for superoxide are less certain. Here, data on complex I kinetics and ROS generation are analyzed using a computational model to determine the sites responsible for superoxide. The analysis includes all the major redox centers: the FMN, iron-sulfur cluster N2, and semiquinone. Analysis reveals that the fully reduced FMN and semiquinone are the primary sources of superoxide, and the iron-sulfur cluster N2 produces none. The FMN radical only produces ROS when the quinone reductase site is blocked. Model simulations reveal that ROS generation is maximized during reverse electron transport with both the FMN and the semiquinone producing similar amounts of superoxide. In addition, the model successfully predicts the increase in ROS generation when the membrane potential is high and matrix pH is alkaline. Of the total ROS produced by complex I, the majority originates from the FMN.

Mathematical models of mitochondrial aging and dynamics

Research on the role of mitochondria in aging and disease is rapidly growing. Furthermore, in recent years, it also became clear that mitochondria are dynamic structures undergoing constant and rapid cycles of fusion and fission. The involvement of mitochondria in multiple complex processes makes them a prime target for mathematical and computational modeling. This review consists of two parts. In the first (Section 2), we provide a detailed introduction to the underlying concepts of mathematical modeling to help the reader who is not so familiar with these techniques to judge the requirements and results that can be obtained through modeling. In the second part (Section 3), we review existing mathematical and computational models that investigate mitochondrial dynamics and the role of mitochondria for the aging process.