p13CMFA: Parsimonious 13C metabolic flux analysis.

Deciphering the mechanisms of regulation of metabolic networks subjected to perturbations, including disease states and drug-induced stress, relies on tracing metabolic fluxes. One of the most informative data to predict metabolic fluxes are 13C based metabolomics, which provide information about how carbons are redistributed along central carbon metabolism. Such data can be integrated using 13C Metabolic Flux Analysis (13C MFA) to provide quantitative metabolic maps of flux distributions. However, 13C MFA might be unable to reduce the solution space towards a unique solution either in large metabolic networks or when small sets of measurements are integrated. Here we present parsimonious 13C MFA (p13CMFA), an approach that runs a secondary optimization in the 13C MFA solution space to identify the solution that minimizes the total reaction flux. Furthermore, flux minimization can be weighted by gene expression measurements allowing seamless integration of gene expression data with 13C data. As proof of concept, we demonstrate how p13CMFA can be used to estimate intracellular flux distributions from 13C measurements and transcriptomics data. We have implemented p13CMFA in Iso2Flux, our in-house developed isotopic steady-state 13C MFA software. The source code is freely available on GitHub (https://github.com/cfoguet/iso2flux/releases/tag/0.7.2).

A key role for mitochondrial gatekeeper pyruvate dehydrogenase in oncogene-induced senescence.

In response to tenacious stress signals, such as the unscheduled activation of oncogenes, cells can mobilize tumour suppressor networks to avert the hazard of malignant transformation. A large body of evidence indicates that oncogene-induced senescence (OIS) acts as such a break, withdrawing cells from the proliferative pool almost irreversibly, thus crafting a vital pathophysiological mechanism that protects against cancer. Despite the widespread contribution of OIS to the cessation of tumorigenic expansion in animal models and humans, we have only just begun to define the underlying mechanism and identify key players. Although deregulation of metabolism is intimately linked to the proliferative capacity of cells, and senescent cells are thought to remain metabolically active, little has been investigated in detail about the role of cellular metabolism in OIS. Here we show, by metabolic profiling and functional perturbations, that the mitochondrial gatekeeper pyruvate dehydrogenase (PDH) is a crucial mediator of senescence induced by BRAF(V600E), an oncogene commonly mutated in melanoma and other cancers. BRAF(V600E)-induced senescence was accompanied by simultaneous suppression of the PDH-inhibitory enzyme pyruvate dehydrogenase kinase 1 (PDK1) and induction of the PDH-activating enzyme pyruvate dehydrogenase phosphatase 2 (PDP2). The resulting combined activation of PDH enhanced the use of pyruvate in the tricarboxylic acid cycle, causing increased respiration and redox stress. Abrogation of OIS, a rate-limiting step towards oncogenic transformation, coincided with reversion of these processes. Further supporting a crucial role of PDH in OIS, enforced normalization of either PDK1 or PDP2 expression levels inhibited PDH and abrogated OIS, thereby licensing BRAF(V600E)-driven melanoma development. Finally, depletion of PDK1 eradicated melanoma subpopulations resistant to targeted BRAF inhibition, and caused regression of established melanomas. These results reveal a mechanistic relationship between OIS and a key metabolic signalling axis, which may be exploited therapeutically.

MIDcor, an R-program for deciphering mass interferences in mass spectra of metabolites enriched in stable isotopes.

BACKGROUND: Tracing stable isotopes, such as 13C using various mass spectrometry (MS) methods provides a valuable information necessary for the study of biochemical processes in cells. However, extracting such information requires special care, such as a correction for naturally occurring isotopes, or overlapping mass spectra of various components of the cell culture medium. Developing a method for a correction of overlapping peaks is the primary objective of this study. RESULTS: Our computer program-MIDcor (free at https://github.com/seliv55/mid_correct) written in the R programming language, corrects the raw MS spectra both for the naturally occurring isotopes and for the overlapping of peaks corresponding to various substances. To this end, the mass spectra of unlabeled metabolites measured in two media are necessary: in a minimal medium containing only derivatized metabolites and chemicals for derivatization, and in a complete cell incubated medium. The MIDcor program calculates the difference (D) between the theoretical and experimentally measured spectra of metabolites containing only the naturally occurring isotopes. The result of comparison of D in the two media determines a way of deciphering the true spectra. (1) If D in the complete medium is greater than that in the minimal medium in at least one peak, then unchanged D is subtracted from the raw spectra of the labeled metabolite. (2) If D does not depend on the medium, then the spectrum probably overlaps with a derivatized fragment of the same metabolite, and D is modified proportionally to the metabolite labeling. The program automatically reaches a decision regarding the way of correction. For some metabolites/fragments in the case (2) D was found to decrease when the tested substance was 13C labeled, and this isotopic effect also can be corrected automatically, if the user provides a measured spectrum of the substance in which the 13C labeling is known a priori. CONCLUSION: Using the developed program improves the reliability of stable isotope tracer data analysis.

Restrictions in ATP diffusion within sarcomeres can provoke ATP-depleted zones impairing exercise capacity in chronic obstructive pulmonary disease.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) is characterized by the inability of patients to sustain a high level of ventilation resulting in perceived exertional discomfort and limited exercise capacity of leg muscles at average intracellular ATP levels sufficient to support contractility. METHODS: Myosin ATPase activity in biopsy samples from healthy and COPD individuals was implemented as a local nucleotide sensor to determine ATP diffusion coefficients within myofibrils. Ergometric parameters clinically measured during maximal exercise tests in both groups were used to define the rates of myosin ATPase reaction and aerobic ATP re-synthesis. The obtained parameters in combination with AK- and CK-catalyzed reactions were implemented to compute the kinetic and steady-state spatial ATP distributions within control and COPD sarcomeres. RESULTS: The developed reaction-diffusion model of two-dimensional sarcomeric space identified similar, yet extremely low nucleotide diffusion in normal and COPD myofibrils. The corresponding spatio-temporal ATP distributions, constructed during imposed exercise, predicted in COPD sarcomeres a depletion of ATP in the zones of overlap between actin and myosin filaments along the center axis at average cytosolic ATP levels similar to healthy muscles. CONCLUSIONS: ATP-depleted zones can induce rigor tension foci impairing muscle contraction and increase a risk for sarcomere damages. Thus, intra-sarcomeric diffusion restrictions at limited aerobic ATP re-synthesis can be an additional risk factor contributing to the muscle contractile deficiency experienced by COPD patients. GENERAL SIGNIFICANCE: This study demonstrates how restricted substrate mobility within a cellular organelle can provoke an energy imbalance state paradoxically occurring at abounding average metabolic resources.

HepatoDyn: A Dynamic Model of Hepatocyte Metabolism That Integrates 13C Isotopomer Data.

The liver performs many essential metabolic functions, which can be studied using computational models of hepatocytes. Here we present HepatoDyn, a highly detailed dynamic model of hepatocyte metabolism. HepatoDyn includes a large metabolic network, highly detailed kinetic laws, and is capable of dynamically simulating the redox and energy metabolism of hepatocytes. Furthermore, the model was coupled to the module for isotopic label propagation of the software package IsoDyn, allowing HepatoDyn to integrate data derived from 13C based experiments. As an example of dynamical simulations applied to hepatocytes, we studied the effects of high fructose concentrations on hepatocyte metabolism by integrating data from experiments in which rat hepatocytes were incubated with 20 mM glucose supplemented with either 3 mM or 20 mM fructose. These experiments showed that glycogen accumulation was significantly lower in hepatocytes incubated with medium supplemented with 20 mM fructose than in hepatocytes incubated with medium supplemented with 3 mM fructose. Through the integration of extracellular fluxes and 13C enrichment measurements, HepatoDyn predicted that this phenomenon can be attributed to a depletion of cytosolic ATP and phosphate induced by high fructose concentrations in the medium.

Relevance of the MEK/ERK signaling pathway in the metabolism of activated macrophages: a metabolomic approach.

The activation of immune cells in response to a pathogen involves a succession of signaling events leading to gene and protein expression, which requires metabolic changes to match the energy demands. The metabolic profile associated with the MAPK cascade (ERK1/2, p38, and JNK) in macrophages was studied, and the effect of its inhibition on the specific metabolic pattern of LPS stimulation was characterized. A [1,2-[(13)C](2)]glucose tracer-based metabolomic approach was used to examine the metabolic flux distribution in these cells after MEK/ERK inhibition. Bioinformatic tools were used to analyze changes in mass isotopomer distribution and changes in glucose and glutamine consumption and lactate production in basal and LPS-stimulated conditions in the presence and absence of the selective inhibitor of the MEK/ERK cascade, PD325901. Results showed that PD325901-mediated ERK1/2 inhibition significantly decreased glucose consumption and lactate production but did not affect glutamine consumption. These changes were accompanied by a decrease in the glycolytic flux, consistent with the observed decrease in fructose-2,6-bisphosphate concentration. The oxidative and nonoxidative pentose phosphate pathways and the ratio between them also decreased. However, tricarboxylic acid cycle flux did not change significantly. LPS activation led to the opposite responses, although all of these were suppressed by PD325901. However, LPS also induced a small decrease in pentose phosphate pathway fluxes and an increase in glutamine consumption that were not affected by PD325901. We concluded that inhibition of the MEK/ERK cascade interferes with central metabolism, and this cross-talk between signal transduction and metabolism also occurs in the presence of LPS.

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.

Compartmentation of membrane processes and nucleotide dynamics in diffusion-restricted cardiac cell microenvironment.

Orchestrated excitation-contraction coupling in heart muscle requires adequate spatial arrangement of systems responsible for ion movement and metabolite turnover. Co-localization of regulatory and transporting proteins into macromolecular complexes within an environment of microanatomical cell components raises intracellular diffusion barriers that hamper the mobility of metabolites and signaling molecules. Compared to substrate diffusion in the cytosol, diffusional restrictions underneath the sarcolemma are much larger and could impede ion and nucleotide movement by a factor of 10(3)-10(5). Diffusion barriers thus seclude metabolites within the submembrane space enabling rapid and vectorial effector targeting, yet hinder energy supply from the bulk cytosolic space implicating the necessity for a shunting transfer mechanism. Here, we address principles of membrane protein compartmentation, phosphotransfer enzyme-facilitated interdomain energy transfer, and nucleotide signal dynamics at the subsarcolemma-cytosol interface. This article is part of a Special Issue entitled "Local Signaling in Myocytes".

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.

MITODYN: An Open Source Software for Quantitative Modeling of Mitochondrial and Cellular Energy Metabolic Flux Dynamics in Health and Disease.

Mitochondrial respiratory chain (RC) transforms the reductive power of NADH or FADH2 oxidation into a proton gradient between the matrix and cytosolic sides of the inner mitochondrial membrane, that ATP synthase uses to generate ATP. This process constitutes a bridge between carbohydrates' central metabolism and ATP-consuming cellular functions. Moreover, the RC is responsible for a large part of reactive oxygen species (ROS) generation that play signaling and oxidizing roles in cells. Mathematical methods and computational analysis are required to understand and predict the possible behavior of this metabolic system. Here we propose a software tool that helps to analyze individual steps of respiratory electron transport in their dynamics, thus deepening understanding of the mechanism of energy transformation and ROS generation in the RC. This software's core is a kinetic model of the RC represented by a system of ordinary differential equations (ODEs). This model enables the analysis of complex dynamic behavior of the RC, including multistationarity and oscillations. The proposed RC modeling method can be applied to study respiration and ROS generation in various organisms and naturally extended to explore carbohydrates' metabolism and linked metabolic processes.