Multistationary and Oscillatory Modes of Free Radicals Generation by the Mitochondrial Respiratory Chain Revealed by a Bifurcation Analysis

The mitochondrial electron transport chain transforms energy satisfying cellular demand and generates reactive oxygen species (ROS) that act as metabolic signals or destructive factors. Therefore, knowledge of the possible modes and bifurcations of electron transport that affect ROS signaling provides insight into the interrelationship of mitochondrial respiration with cellular metabolism. Here, a bifurcation analysis of a sequence of the electron transport chain models of increasing complexity was used to analyze the contribution of individual components to the modes of respiratory chain behavior. Our algorithm constructed models as large systems of ordinary differential equations describing the time evolution of the distribution of redox states of the respiratory complexes. The most complete model of the respiratory chain and linked metabolic reactions predicted that condensed mitochondria produce more ROS at low succinate concentration and less ROS at high succinate levels than swelled mitochondria. This prediction was validated by measuring ROS production under various swelling conditions. A numerical bifurcation analysis revealed qualitatively different types of multistationary behavior and sustained oscillations in the parameter space near a region that was previously found to describe the behavior of isolated mitochondria. The oscillations in transmembrane potential and ROS generation, observed in living cells were reproduced in the model that includes interaction of respiratory complexes with the reactions of TCA cycle. Whereas multistationarity is an internal characteristic of the respiratory chain, the functional link of respiration with central metabolism creates oscillations, which can be understood as a means of auto-regulation of cell metabolism. © 2012 Selivanov et al

Unveiling a key role of oxaloacetate-glutamate interaction in regulation of respiration and ROS generation in nonsynaptic brain mitochondria using a kinetic model.

Glutamate plays diverse roles in neuronal cells, affecting cell energetics and reactive oxygen species (ROS) generation. These roles are especially vital for neuronal cells, which deal with high amounts of glutamate as a neurotransmitter. Our analysis explored neuronal glutamate implication in cellular energy metabolism and ROS generation, using a kinetic model that simulates electron transport details in respiratory complexes, linked ROS generation and metabolic reactions. The analysis focused on the fact that glutamate attenuates complex II inhibition by oxaloacetate, stimulating the latter's transformation into aspartate. Such a mechanism of complex II activation by glutamate could cause almost complete reduction of ubiquinone and deficiency of oxidized form (Q), which closes the main stream of electron transport and opens a way to massive ROS generating transfer in complex III from semiquinone radicals to molecular oxygen. In this way, under low workload, glutamate triggers the respiratory chain (RC) into a different steady state characterized by high ROS generation rate. The observed stepwise dependence of ROS generation on glutamate concentration experimentally validated this prediction. However, glutamate's attenuation of oxaloacetate's inhibition accelerates electron transport under high workload. Glutamate-oxaloacetate interaction in complex II regulation underlies the observed effects of uncouplers and inhibitors and acceleration of Ca2+ uptake. Thus, this theoretical analysis uncovered the previously unknown roles of oxaloacetate as a regulator of ROS generation and glutamate as a modifier of this regulation. The model predicted that this mechanism of complex II activation by glutamate might be operative in situ and responsible for excitotoxicity. Spatial-time gradients of synthesized hydrogen peroxide concentration, calculated in the reaction-diffusion model with convection under a non-uniform local approximation of nervous tissue, have shown that overproduction of H2O2 in a cell causes excess of its level in neighbor cells

Integration of enzyme kinetic models and isotopomer distribution analysis for studies of in situ cell operation

A current trend in neuroscience research is the use of stable isotope tracers in order to address metabolic processes in vivo. The tracers produce a huge number of metabolite forms that differ according to the number and position of labeled isotopes in the carbon skeleton (isotopomers) and such a large variety makes the analysis of isotopomer data highly complex. On the other hand, this multiplicity of forms does provide sufficient information to address cell operation in vivo. By the end of last millennium, a number of tools have been developed for estimation of metabolic flux profile from any possible isotopomer distribution data. However, although well elaborated, these tools were limited to steady state analysis, and the obtained set of fluxes remained disconnected from their biochemical context. In this review we focus on a new numerical analytical approach that integrates kinetic and metabolic flux analysis. The related computational algorithm estimates the dynamic flux based on the time-dependent distribution of all possible isotopomers of metabolic pathway intermediates that are generated from a labeled substrate. The new algorithm connects specific tracer data with enzyme kinetic characteristics, thereby extending the amount of data available for analysis: it uses enzyme kinetic data to estimate the flux profile, and vice versa, for the kinetic analysis it uses in vivo tracer data to reveal the biochemical basis of the estimated metabolic fluxes

Edelfosine-induced metabolic changes in cancer cells that precede the overproduction of reactive oxygen species and apoptosis

Background: Metabolic flux profiling based on the analysis of distribution of stable isotope tracer in metabolites is an important method widely used in cancer research to understand the regulation of cell metabolism and elaborate new therapeutic strategies. Recently, we developed software Isodyn, which extends the methodology of kinetic modeling to the analysis of isotopic isomer distribution for the evaluation of cellular metabolic flux profile under relevant conditions. This tool can be applied to reveal the metabolic effect of proapoptotic drug edelfosine in leukemia Jurkat cell line, uncovering the mechanisms of induction of apoptosis in cancer cells. Results: The study of 13C distribution of Jukat cells exposed to low edelfosine concentration, which induces apoptosis in ¿5% of cells, revealed metabolic changes previous to the development of apoptotic program. Specifically, it was found that low dose of edelfosine stimulates the TCA cycle. These metabolic perturbations were coupled with an increase of nucleic acid synthesis de novo, which indicates acceleration of biosynthetic and reparative processes. The further increase of the TCA cycle fluxes, when higher doses of drug applied, eventually enhance reactive oxygen species (ROS) production and trigger apoptotic program. Conclusion: The application of Isodyn to the analysis of mechanism of edelfosine-induced apoptosis revealed primary drug-induced metabolic changes, which are important for the subsequent initiation of apoptotic program. Initiation of such metabolic changes could be exploited in anticancer therapy

Compartmentation of glycogen metabolism revealed from 13C isotopologue distributions

Background: Stable isotope tracers are used to assess metabolic flux profiles in living cells. The existing methods of measurement average out the isotopic isomer distribution in metabolites throughout the cell, whereas the knowledge of compartmental organization of analyzed pathways is crucial for the evaluation of true fluxes. That is why we accepted a challenge to create a software tool that allows deciphering the compartmentation of metabolites based on the analysis of average isotopic isomer distribution. Results: The software Isodyn, which simulates the dynamics of isotopic isomer distribution in central metabolic pathways, was supplemented by algorithms facilitating the transition between various analyzed metabolic schemes, and by the tools for model discrimination. It simulated 13C isotope distributions in glucose, lactate, glutamate and glycogen, measured by mass spectrometry after incubation of hepatocytes in the presence of only labeled glucose or glucose and lactate together (with label either in glucose or lactate). The simulations assumed either a single intracellular hexose phosphate pool, or also channeling of hexose phosphates resulting in a different isotopic composition of glycogen. Model discrimination test was applied to check the consistency of both models with experimental data. Metabolic flux profiles, evaluated with the accepted model that assumes channeling, revealed the range of changes in metabolic fluxes in liver cells. Conclusions: The analysis of compartmentation of metabolic networks based on the measured 13C distribution was included in Isodyn as a routine procedure. The advantage of this implementation is that, being a part of evaluation of metabolic fluxes, it does not require additional experiments to study metabolic compartmentation. The analysis of experimental data revealed that the distribution of measured 13C-labeled glucose metabolites is inconsistent with the idea of perfect mixing of hexose phosphates in cytosol. In contrast, the observed distribution indicates the presence of a separate pool of hexose phosphates that is channeled towards glycogen synthesis

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

Systems Medicine: from molecular features and models to the clinic in COPD

Background and hypothesis Chronic Obstructive Pulmonary Disease (COPD) patients are characterized by heterogeneous clinical manifestations and patterns of disease progression. Two major factors that can be used to identify COPD subtypes are muscle dysfunction/wasting and co-morbidity patterns. We hypothesized that COPD heterogeneity is in part the result of complex interactions between several genes and pathways. We explored the possibility of using a Systems Medicine approach to identify such pathways, as well as to generate predictive computational models that may be used in clinic practice. Objective and method Our overarching goal is to generate clinically applicable predictive models that characterize COPD heterogeneity through a Systems Medicine approach. To this end we have developed a general framework, consisting of three steps/objectives: (1) feature identification, (2) model generation and statistical validation, and (3) application and validation of the predictive models in the clinical scenario. We used muscle dysfunction and co-morbidity as test cases for this framework. Results: In the study of muscle wasting we identified relevant features (genes) by a network analysis and generated predictive models that integrate mechanistic and probabilistic models. This allowed us to characterize muscle wasting as a general de-regulation of pathway interactions. In the co-morbidity analysis we identified relevant features (genes/pathways) by the integration of gene-disease and disease-disease associations. We further present a detailed characterization of co-morbidities in COPD patients that was implemented into a predictive model. In both use cases we were able to achieve predictive modeling but we also identified several key challenges, the most pressing being the validation and implementation into actual clinical practice. Conclusions: The results confirm the potential of the Systems Medicine approach to study complex diseases and generate clinically relevant predictive models. Our study also highlights important obstacles and bottlenecks for such approaches (e.g. data availability and normalization of frameworks among others) and suggests specific proposals to overcome them

Systems Medicine: from molecular features and models to the clinic in COPD

BACKGROUND AND HYPOTHESIS: Chronic Obstructive Pulmonary Disease (COPD) patients are characterized by heterogeneous clinical manifestations and patterns of disease progression. Two major factors that can be used to identify COPD subtypes are muscle dysfunction/wasting and co-morbidity patterns. We hypothesized that COPD heterogeneity is in part the result of complex interactions between several genes and pathways. We explored the possibility of using a Systems Medicine approach to identify such pathways, as well as to generate predictive computational models that may be used in clinic practice. OBJECTIVE AND METHOD: Our overarching goal is to generate clinically applicable predictive models that characterize COPD heterogeneity through a Systems Medicine approach. To this end we have developed a general framework, consisting of three steps/objectives: (1) feature identification, (2) model generation and statistical validation, and (3) application and validation of the predictive models in the clinical scenario. We used muscle dysfunction and co-morbidity as test cases for this framework. RESULTS: In the study of muscle wasting we identified relevant features (genes) by a network analysis and generated predictive models that integrate mechanistic and probabilistic models. This allowed us to characterize muscle wasting as a general de-regulation of pathway interactions. In the co-morbidity analysis we identified relevant features (genes/pathways) by the integration of gene-disease and disease-disease associations. We further present a detailed characterization of co-morbidities in COPD patients that was implemented into a predictive model. In both use cases we were able to achieve predictive modeling but we also identified several key challenges, the most pressing being the validation and implementation into actual clinical practice. CONCLUSIONS: The results confirm the potential of the Systems Medicine approach to study complex diseases and generate clinically relevant predictive models. Our study also highlights important obstacles and bottlenecks for such approaches (e.g. data availability and normalization of frameworks among others) and suggests specific proposals to overcome them