Bistability of Mitochondrial Respiration Underlies Paradoxical Reactive Oxygen Species Generation Induced by Anoxia

Increased production of reactive oxygen species (ROS) in mitochondria underlies major systemic diseases, and this clinical problem stimulates a great scientific interest in the mechanism of ROS generation. However, the mechanism of hypoxia-induced change in ROS production is not fully understood. To mathematically analyze this mechanism in details, taking into consideration all the possible redox states formed in the process of electron transport, even for respiratory complex III, a system of hundreds of differential equations must be constructed. Aimed to facilitate such tasks, we developed a new methodology of modeling, which resides in the automated construction of large sets of differential equations. The detailed modeling of electron transport in mitochondria allowed for the identification of two steady state modes of operation (bistability) of respiratory complex III at the same microenvironmental conditions. Various perturbations could induce the transition of respiratory chain from one steady state to another. While normally complex III is in a low ROS producing mode, temporal anoxia could switch it to a high ROS producing state, which persists after the return to normal oxygen supply. This prediction, which we qualitatively validated experimentally, explains the mechanism of anoxia-induced cell damage. Recognition of bistability of complex III operation may enable novel therapeutic strategies for oxidative stress and our method of modeling could be widely used in systems biology studies

Reactive Oxygen Species Production by Forward and Reverse Electron Fluxes in the Mitochondrial Respiratory Chain

Reactive oxygen species (ROS) produced in the mitochondrial respiratory chain (RC) are primary signals that modulate cellular adaptation to environment, and are also destructive factors that damage cells under the conditions of hypoxia/reoxygenation relevant for various systemic diseases or transplantation. The important role of ROS in cell survival requires detailed investigation of mechanism and determinants of ROS production. To perform such an investigation we extended our rule-based model of complex III in order to account for electron transport in the whole RC coupled to proton translocation, transmembrane electrochemical potential generation, TCA cycle reactions, and substrate transport to mitochondria. It fits respiratory electron fluxes measured in rat brain mitochondria fueled by succinate or pyruvate and malate, and the dynamics of NAD+ reduction by reverse electron transport from succinate through complex I. The fitting of measured characteristics gave an insight into the mechanism of underlying processes governing the formation of free radicals that can transfer an unpaired electron to oxygen-producing superoxide and thus can initiate the generation of ROS. Our analysis revealed an association of ROS production with levels of specific radicals of individual electron transporters and their combinations in species of complexes I and III. It was found that the phenomenon of bistability, revealed previously as a property of complex III, remains valid for the whole RC. The conditions for switching to a state with a high content of free radicals in complex III were predicted based on theoretical analysis and were confirmed experimentally. These findings provide a new insight into the mechanisms of ROS production in RC

13C labeling experiments at metabolic nonstationary conditions: An exploratory study

<p>Abstract</p> <p>Background</p> <p>Stimulus Response Experiments to unravel the regulatory properties of metabolic networks are becoming more and more popular. However, their ability to determine enzyme kinetic parameters has proven to be limited with the presently available data. In metabolic flux analysis, the use of ¹³C labeled substrates together with isotopomer modeling solved the problem of underdetermined networks and increased the accuracy of flux estimations significantly.</p> <p>Results</p> <p>In this contribution, the idea of increasing the information content of the dynamic experiment by adding ¹³C labeling is analyzed. For this purpose a small example network is studied by simulation and statistical methods. Different scenarios regarding available measurements are analyzed and compared to a non-labeled reference experiment. Sensitivity analysis revealed a specific influence of the kinetic parameters on the labeling measurements. Statistical methods based on parameter sensitivities and different measurement models are applied to assess the information gain of the labeled stimulus response experiment.</p> <p>Conclusion</p> <p>It was found that the use of a (specifically) labeled substrate will significantly increase the parameter estimation accuracy. An overall information gain of about a factor of six is observed for the example network. The information gain is achieved from the specific influence of the kinetic parameters towards the labeling measurements. This also leads to a significant decrease in correlation of the kinetic parameters compared to an experiment without ¹³C-labeled substrate.</p

On the effects of mechanical stress of biological membranes in modeling of swelling dynamics of biological systems

We highlight mechanical stretching and bending of membranes and the importance of membrane deformations in the analysis of swelling dynamics of biological systems, including cells and subcellular organelles. Membrane deformation upon swelling generates tensile stress and internal pressure, contributing to volume changes in biological systems. Therefore, in addition to physical (internal/external) and chemical factors, mechanical properties of the membranes should be considered in modeling analysis of cellular swelling. Here we describe an approach that considers mechanical properties of the membranes in the analysis of swelling dynamics of biological systems. This approach includes membrane bending and stretching deformations into the model, producing a more realistic description of swelling. We also discuss the effects of membrane stretching on swelling dynamics. We report that additional pressure generated by membrane bending is negligible, compared to pressures generated by membrane stretching, when both membrane surface area and volume are variable parameters. Note that bending deformations are reversible, while stretching deformation may be irreversible, leading to membrane disruption when they exceed a certain threshold level. Therefore, bending deformations need only be considered in reversible physiological swelling, whereas stretching deformations should also be considered in pathological irreversible swelling. Thus, the currently proposed approach may be used to develop a detailed biophysical model describing the transition from physiological to pathological swelling mode.National Aeronautics & Space Administration (NASA):80NSSC19M0049; PR Space Grant (NASA):NNX15AI11Hinfo:eu-repo/semantics/publishedVersio

A mathematical model of mitochondrial swelling

<p>Abstract</p> <p>Background</p> <p>The <it>permeabilization </it>of mitochondrial membranes is a decisive event in apoptosis or necrosis culminating in cell death. One fundamental mechanism by which such permeabilization events occur is the calcium-induced mitochondrial permeability transition. Upon Ca²⁺-uptake into mitochondria an increase in inner membrane permeability occurs by a yet unclear mechanism. This leads to a net water influx in the mitochondrial matrix, mitochondrial swelling, and finally the rupture of the outer membrane. Although already described more than thirty years ago, many unsolved questions surround this important biological phenomenon. Importantly, theoretical modeling of the mitochondrial permeability transition has only started recently and the existing mathematical models fail to characterize the swelling process throughout the whole time range.</p> <p>Results</p> <p>We propose here a new mathematical approach to the mitochondrial permeability transition introducing a specific delay equation and resulting in an optimized representation of mitochondrial swelling. Our new model is in accordance with the experimentally determined course of volume increase throughout the whole swelling process, including its initial lag phase as well as its termination. From this new model biological consequences can be deduced, such as the confirmation of a positive feedback of mitochondrial swelling which linearly depends on the Ca²⁺-concentration, or a negative exponential dependence of the average swelling time on the Ca²⁺-concentration. Finally, our model can show an initial shrinking phase of mitochondria, which is often observed experimentally before the actual swelling starts.</p> <p>Conclusions</p> <p>We present a model of the mitochondrial swelling kinetics. This model may be adapted and extended to diverse other inducing/inhibiting conditions or to mitochondria from other biological sources and thus may benefit a better understanding of the mitochondrial permeability transition.</p

The cystic fibrosis transmembrane recruiter the alter ego of CFTR as a multi-kinase anchor

This review focuses on a newly discovered interaction between protein kinases involved in cellular energetics, a process that may be disturbed in cystic fibrosis for unknown reasons. I propose a new model where kinase-mediated cellular transmission of energy provides mechanistic insight to a latent role of the cystic fibrosis transmembrane conductance regulator (CFTR). I suggest that CFTR acts as a multi-kinase recruiter to the apical epithelial membrane. My group finds that, in the cytosol, two protein kinases involved in cell energy homeostasis, nucleoside diphosphate kinase (NDPK) and AMP-activated kinase (AMPK), bind one another. Preliminary data suggest that both can also bind CFTR (function unclear). The disrupted role of this CFTR-kinase complex as ‘membrane transmitter to the cell’ is proposed as an alternative paradigm to the conventional ion transport mediated and CFTR/chloride-centric view of cystic fibrosis pathogenesis. Chloride remains important, but instead, chloride-induced control of the phosphohistidine content of one kinase component (NDPK, via a multi-kinase complex that also includes a third kinase, CK2; formerly casein kinase 2). I suggest that this complex provides the necessary near-equilibrium conditions needed for efficient transmission of phosphate energy to proteins controlling cellular energetics. Crucially, a new role for CFTR as a kinase controller is proposed with ionic concentration acting as a signal. The model posits a regulatory control relay for energy sensing involving a cascade of protein kinases bound to CFTR

Membrane potential and delta pH dependency of reverse electron transport-associated hydrogen peroxide production in brain and heart mitochondria

Succinate-driven reverse electron transport (RET) is one of the main sources of mitochondrial reactive oxygen species (mtROS) in ischemia-reperfusion injury. RET is dependent on mitochondrial membrane potential (Δψm) and transmembrane pH difference (ΔpH), components of the proton motive force (pmf); a decrease in Δψm and/or ΔpH inhibits RET. In this study we aimed to determine which component of the pmf displays the more dominant effect on RET-provoked ROS generation in isolated guinea pig brain and heart mitochondria respiring on succinate or α-glycerophosphate (α-GP). Δψm was detected via safranin fluorescence and a TPP+ electrode, the rate of H2O2 formation was measured by Amplex UltraRed, the intramitochondrial pH (pHin) was assessed via BCECF fluorescence. Ionophores were used to dissect the effects of the two components of pmf. The K+/H+ exchanger, nigericin lowered pHin and ΔpH, followed by a compensatory increase in Δψm that led to an augmented H2O2 production. Valinomycin, a K+ ionophore, at low [K+] increased ΔpH and pHin, decreased Δψm, which resulted in a decline in H2O2 formation. It was concluded that Δψm is dominant over ∆pH in modulating the succinate- and α-GP-evoked RET. The elevation of extramitochondrial pH was accompanied by an enhanced H2O2 release and a decreased ∆pH. This phenomenon reveals that from the pH component not ∆pH, but rather absolute value of pH has higher impact on the rate of mtROS formation. Minor decrease of Δψm might be applied as a therapeutic strategy to attenuate RET-driven ROS generation in ischemia-reperfusion injury