Search

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., (© 2022. This is a U.S. government work and not under copyright protection in the U.S.; foreign copyright protection may apply.)

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., Competing Interests: The authors have declared that no competing interests exist.

Stable isotope-resolved metabolomics (SIRM), based on the analysis of biological samples from living cells incubated with artificial isotope enriched substrates, enables mapping the rates of biochemical reactions (metabolic fluxes). We developed software supporting a workflow of analysis of SIRM data obtained with mass spectrometry (MS). The evaluation of fluxes starting from raw MS recordings requires at least three steps of computer support: first, extraction of mass spectra of metabolites of interest, then correction of the spectra for natural isotope abundance, and finally, evaluation of fluxes by simulation of the corrected spectra using a corresponding mathematical model. A kinetic model based on ordinary differential equations (ODEs) for isotopomers of metabolites of the corresponding biochemical network supports the final part of the analysis, which provides a dynamic flux map.

Altered metabolism is a hallmark of cancer, but little is still known about its regulation. In this study, we measure transcriptomic, proteomic, phospho-proteomic and fluxomics data in a breast cancer cell-line (MCF7) across three different growth conditions. Integrating these multiomics data within a genome scale human metabolic model in combination with machine learning, we systematically chart the different layers of metabolic regulation in breast cancer cells, predicting which enzymes and pathways are regulated at which level. We distinguish between two types of reactions, directly and indirectly regulated. Directly-regulated reactions include those whose flux is regulated by transcriptomic alterations (~890) or via proteomic or phospho-proteomics alterations (~140) in the enzymes catalyzing them. We term the reactions that currently lack evidence for direct regulation as (putative) indirectly regulated (~930). Many metabolic pathways are predicted to be regulated at different levels, and those may change at different media conditions. Remarkably, we find that the flux of predicted indirectly regulated reactions is strongly coupled to the flux of the predicted directly regulated ones, uncovering a tiered hierarchical organization of breast cancer cell metabolism. Furthermore, the predicted indirectly regulated reactions are predominantly reversible. Taken together, this architecture may facilitate rapid and efficient metabolic reprogramming in response to the varying environmental conditions incurred by the tumor cells. The approach presented lays a conceptual and computational basis for mapping metabolic regulation in additional cancers.

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)., Competing Interests: The authors have declared that no competing interests exist.

Background: Tracing stable isotopes, such as 13 C 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 13 C labeled, and this isotopic effect also can be corrected automatically, if the user provides a measured spectrum of the substance in which the 13 C labeling is known a priori., Conclusion: Using the developed program improves the reliability of stable isotope tracer data analysis.

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., (Copyright © 2016 Elsevier B.V. All rights reserved.)

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.

Background: It has been suggested that the adipokine resistin links obesity and insulin resistance, although how resistin acts on muscle metabolism is controversial. We aimed to quantitatively analyse the effects of resistin on the glucose metabolic flux profile and on insulin response in L6E9 myotubes at the metabolic level using a tracer-based metabolomic approach and our in-house developed software, Isodyn., Results: Resistin significantly increased glucose uptake and glycolysis, altering pyruvate utilisation by the cell. In the presence of resistin, insulin only slightly increased glucose uptake and glycolysis, and did not alter the flux profile around pyruvate induced by resistin. Resistin prevented the increase in gene expression in pyruvate dehydrogenase-E1 and the sharp decrease in gene expression in cytosolic phosphoenolpyruvate carboxykinase-1 induced by insulin., Conclusions: These data suggest that resistin impairs the metabolic activation of insulin. This impairment cannot be explained by the activity of a single enzyme, but instead due to reorganisation of the whole metabolic flux distribution.

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.

Catalytic activity has been demonstrated for holotransketolase in the absence of free bivalent cations in the medium. The two active centers of the enzyme are equivalent in both the catalytic activity and the affinity for the substrates. In the presence of free Ca²⁺ (added to the medium from an external source), this equivalence is lost: negative cooperativity is induced on binding of either xylulose 5-phosphate (donor substrate) or ribose 5-phosphate (acceptor substrate), whereupon the catalytic conversion of the bound substrates causes the interaction between the centers to become positively cooperative. Moreover, the enzyme total activity increase is observed.

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.

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"., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

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.

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.

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.

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.

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.

Catalytic activity of two active sites of transketolase and their affinity towards the substrates (xylulose-5-phosphate and ribose-5-phosphate) has been studied in the presence of Ca2+ and Mg2+. In the presence of Ca2+, the active sites exhibit negative cooperativity in binding both xylulose-5-phosphate (donor substrate) and ribose-5-phosphate (acceptor substrate) and positive cooperativity in the catalytic transformation of the substrates. In the presence of Mg2+, nonequivalence of the active sites is not observed.

Reactive oxygen species (ROS) generation in mitochondria as a side product of electron and proton transport through the inner membrane is important for normal cell operation as well as development of pathology. Matrix and cytosol alkalization stabilizes semiquinone radical, a potential superoxide producer, and we hypothesized that proton deficiency under the excess of electron donors enhances reactive oxygen species generation. We tested this hypothesis by measuring pH dependence of reactive oxygen species released by mitochondria. The experiments were performed in the media with pH varying from 6 to 8 in the presence of complex II substrate succinate or under more physiological conditions with complex I substrates glutamate and malate. Matrix pH was manipulated by inorganic phosphate, nigericine, and low concentrations of uncoupler or valinomycin. We found that high pH strongly increased the rate of free radical generation in all of the conditions studied, even when DeltapH=0 in the presence of nigericin. In the absence of inorganic phosphate, when the matrix was the most alkaline, pH shift in the medium above 7 induced permeability transition accompanied by the decrease of ROS production. ROS production increase induced by the alkalization of medium was observed with intact respiring mitochondria as well as in the presence of complex I inhibitor rotenone, which enhanced reactive oxygen species release. The phenomena revealed in this report are important for understanding mechanisms governing mitochondrial production of reactive oxygen species, in particular that related with uncoupling proteins.

The biochemical effects of training programmes have been studied with a kinetic model of central metabolism, using enzyme activities and metabolite concentrations measured at rest and after 30 s maximum-intensity exercise, collected before and after long and short periods of training, which differed only by the duration of the rest intervals. After short periods of training the glycolytic flux at rest was three times higher than it had been before training, whereas during exercise the flux and energy consumption remained the same as before training. Long periods of training had less effect on the glycolytic flux at rest, but increased it in response to exercise, increasing the contribution of oxidative phosphorylation.

Although a high ATP diffusion rate implies homogeneous distribution of the principal energetic currency in the cytosol, local diffusion barriers represented by macromolecular structures can render ATP concentrations to be inhomogeneous. A method is presented here that provides apparent diffusion coefficient values in local intracellular regions and allows the estimation of spatial metabolite distribution. The apparent local diffusion coefficient for ATP in cardiac myofibrils was determined from the analysis of diffusion-dependent rightward shift of the substrate dependence for actomyosin ATPase activity using the reaction-diffusion model, which accounted for the properties of phosphotransfer reactions. This functional analysis, which took into account the local diffusional ATP delivery to the active sites, provided an apparent value that was three orders of magnitude lower than that defined by direct methods for the cytosol. The low value of the diffusion coefficient was shown to define unusual properties of the intracellular space in working heart, where small reductions in ATP levels in the surrounding cytosol result in a large drop in [ATP] inside myofibrils. This drop is critical for vital cellular functions, and the analysis presented here defines its physical basis. The diffusion barriers thus defined explain the coexistence of pathological energy deficit with almost normal average ATP levels.

Motivation: Metabolic flux analysis of biochemical reaction networks using isotope tracers requires software tools that can analyze the dynamics of isotopic isomer (isotopomer) accumulation in metabolites and reveal the underlying kinetic mechanisms of metabolism regulation. Since existing tools are restricted by the isotopic steady state and remain disconnected from the underlying kinetic mechanisms, we have recently developed a novel approach for the analysis of tracer-based metabolomic data that meets these requirements. The present contribution describes the last step of this development: implementation of (i) the algorithms for the determination of the kinetic parameters and respective metabolic fluxes consistent with the experimental data and (ii) statistical analysis of both fluxes and parameters, thereby lending it a practical application., Results: The C++ applications package for dynamic isotopomer distribution data analysis was supplemented by (i) five distinct methods for resolving a large system of differential equations; (ii) the 'simulated annealing' algorithm adopted to estimate the set of parameters and metabolic fluxes, which corresponds to the global minimum of the difference between the computed and measured isotopomer distributions; and (iii) the algorithms for statistical analysis of the estimated parameters and fluxes, which use the covariance matrix evaluation, as well as Monte Carlo simulations. An example of using this tool for the analysis of (13)C distribution in the metabolites of glucose degradation pathways has demonstrated the evaluation of optimal set of parameters and fluxes consistent with the experimental pattern, their range and statistical significance, and also the advantages of using dynamic rather than the usual steady-state method of analysis., Availability: Software is available free from http://www.bq.ub.es/bioqint/selivanov.htm

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.

The main object of this work is to present simplified and three-dimensional views of the pentose phosphate and Calvin cycles, emphasizing their functional and chemical similarities., (Copyright © 2006 International Union of Biochemistry and Molecular Biology, Inc.)

Motivation: Addition of labeled substrates and the measurement of the subsequent distribution of the labels in isotopomers in reaction networks provide a unique method for assessing metabolic fluxes in whole cells. However, owing to insufficiency of information, attempts to quantify the fluxes often yield multiple possible sets of solutions that are consistent with a given experimental pattern of isotopomers. In the study of the pentose phosphate pathways, the need to consider isotope exchange reactions of transketolase (TK) and transaldolase (TA) (which in past analyses have often been ignored) magnifies this problem; but accounting for the interrelation between the fluxes known from biochemical studies and kinetic modeling solves it. The mathematical relationships between kinetic and equilibrium constants restrict the domain of estimated fluxes to the ones compatible not only with a given set of experimental data, but also with other biochemical information., Method: We present software that integrates kinetic modeling with isotopomer distribution analysis. It solves the ordinary differential equations for total concentrations (accounting for the kinetic mechanisms) as well as for all isotopomers in glycolysis and the pentose phosphate pathway (PPP). In the PPP the fluxes created in the TK and TA reactions are expressed through unitary rate constants. The algorithms that account for all the kinetic and equilbrium constant constraints are integrated with the previously developed algorithms, which have been further optimized. The most time-consuming calculations were programmed directly in assembly language; this gave an order of magnitude decrease in the computation time, thus allowing analysis of more complex systems. The software was developed as C-code linked to a program written in Mathematica (Wolfram Research, Champaign, IL), and also as a C++ program independent from Mathematica., Results: Implementing constraints imposed by kinetic and equilibrium constants in the isotopomer distribution analysis in the data from the cancer cells eliminated estimates of fluxes that were inconsistent with the kinetic mechanisms of TK and TA. Fluxes measured experimentally in cells can be used to estimate better the kinetics of TK and TA as they operate in situ. Thus, our approach of integrating various methods for in situ flux analysis opens up the possibility of designing new types of experiments to probe metabolic interrelationships, including the incorporation of additional biochemical information., Availability: Software is available freely at: http://www.bq.ub.es/bioqint/selivanov.htm, Contact: martacascante@ub.edu

Motivation: Analysis of the conversion of (13)C glucose within the metabolic network allows the evaluation of the biochemical fluxes in interconnecting metabolic pathways. Such analyses require solving hundreds of equations with respect to individual isotopomer concentrations, and this assumes applying special software even for constructing the equations. The algorithm, proposed by others could be improved., Method: A C-code linked to the program written in Mathematica (Wolfram Research Inc.), constructs and solves differential equations for all isotopomer concentrations, using the general enzyme characteristics (K(m), equilibrium constant, etc.). This code uses innovative algorithm of determination for the isotopomers-products, thus essentially decreasing the computation time. Feasible metabolic fluxes are provided by the parameters of enzyme kinetics found from the data fitting., Results: The software effectively evaluates metabolic fluxes based on the measured isotopomer distribution, as was illustrated by the analysis of glycolysis and pentose phosphate cycle. The mechanism of transketolase and transaldolase catalysis was shown to induce a specific kind of isotopomer re-distribution, which, despite the significance of its effect, usually is not taken into account., Availability: The software could be freely downloaded from the site: http://bq.ub.es/bioqint/label_distribution/.

Data from site-directed mutagenesis and X-ray crystallography show that His103 of holotransketolase (holoTK) does not come into contact with thiamin diphosphate (ThDP) but stabilizes the transketolase (TK) reaction intermediate, alpha,beta-dihydroxyethyl-thiamin diphosphate, by forming a hydrogen bond with the oxygen of its beta-hydroxyethyl group [Eur. J. Biochem. 233 (1995) 750; Proc. Natl. Acad. Sci. USA 99 (2002) 591]. We studied the influence of His103 mutation on ThDP-binding and enzymatic activity. It was found that mutation does not affect the affinity of the coenzyme to apotransketolase (apoTK) in the presence of Ca(2+) (a cation found in the native holoenzyme) but changes all the kinetic parameters of the ThDP-apoTK interaction in the presence of Mg(2+) (a cation commonly used in ThDP-dependent enzymes studies). It was concluded that the structures of TK active centers formed in the presence of Mg(2+) and Ca(2+) are not identical. Mutation of His103 led to a significant acceleration of the one-substrate reaction but a slow down of the two-substrate reaction so that the rates of both types of catalysis became equal. Our results provide evidence for the intermediate-stabilizing function of His103.

Medical records of 320 persons, who were injured in terrorist attacks in 1999-2003, were analyzed. The main errors in rendering the anesthetic-and-resuscitation care were pointed out. An algorithm of actions is elaborated for the anesthesiologist working at the reception, surgical and intensive care departments of a hospital. Owing to the mentioned algorithm, the lethality rate of wounded persons dropped from 18 to 10%.

Stress tolerance of the heart requires high-fidelity metabolic sensing by ATP-sensitive potassium (K(ATP)) channels that adjust membrane potential-dependent functions to match cellular energetic demand. Scanning of genomic DNA from individuals with heart failure and rhythm disturbances due to idiopathic dilated cardiomyopathy identified two mutations in ABCC9, which encodes the regulatory SUR2A subunit of the cardiac K(ATP) channel. These missense and frameshift mutations mapped to evolutionarily conserved domains adjacent to the catalytic ATPase pocket within SUR2A. Mutant SUR2A proteins showed aberrant redistribution of conformations in the intrinsic ATP hydrolytic cycle, translating into abnormal K(ATP) channel phenotypes with compromised metabolic signal decoding. Defective catalysis-mediated pore regulation is thus a mechanism for channel dysfunction and susceptibility to dilated cardiomyopathy.

Transmission of energetic signals to membrane sensors, such as the ATP-sensitive K+ (KATP) channel, is vital for cellular adaptation to stress. Yet, cell compartmentation implies diffusional hindrances that hamper direct reception of cytosolic energetic signals. With high intracellular ATP levels, KATP channels may sense not bulk cytosolic, but rather local submembrane nucleotide concentrations set by membrane ATPases and phosphotransfer enzymes. Here, we analyzed the role of adenylate kinase and creatine kinase phosphotransfer reactions in energetic signal transmission over the strong diffusional barrier in the submembrane compartment, and translation of such signals into a nucleotide response detectable by KATP channels. Facilitated diffusion provided by creatine kinase and adenylate kinase phosphotransfer dissipated nucleotide gradients imposed by membrane ATPases, and shunted diffusional restrictions. Energetic signals, simulated as deviation of bulk ATP from its basal level, were amplified into an augmented nucleotide response in the submembrane space due to failure under stress of creatine kinase to facilitate nucleotide diffusion. Tuning of creatine kinase-dependent amplification of the nucleotide response was provided by adenylate kinase capable of adjusting the ATP/ADP ratio in the submembrane compartment securing adequate KATP channel response in accord with cellular metabolic demand. Thus, complementation between creatine kinase and adenylate kinase systems, here predicted by modeling and further supported experimentally, provides a mechanistic basis for metabolic sensor function governed by alterations in intracellular phosphotransfer fluxes.

Transduction of metabolic signals is essential in preserving cellular homeostasis. Yet, principles governing integration and synchronization of membrane metabolic sensors with cell metabolism remain elusive. Here, analysis of cellular nucleotide fluxes and nucleotide-dependent gating of the ATP-sensitive K+ (K(ATP)) channel, a prototypic metabolic sensor, revealed a diffusional barrier within the submembrane space, preventing direct reception of cytosolic signals. Creatine kinase phosphotransfer, captured by 18O-assisted 31P NMR, coordinated tightly with ATP turnover, reflecting the cellular energetic status. The dynamics of high energy phosphoryl transfer through the creatine kinase relay permitted a high fidelity transmission of energetic signals into the submembrane compartment synchronizing K(ATP) channel activity with cell metabolism. Knock-out of the creatine kinase M-CK gene disrupted signal delivery to K(ATP) channels and generated a cellular phenotype with increased electrical vulnerability. Thus, in the compartmentalized cell environment, phosphotransfer systems shunt diffusional barriers and secure regimented signal transduction integrating metabolic sensors with the cellular energetic network.

Two substrates of the transketolase reaction are known to bind with the enzyme according to a ping-pong mechanism [1]. It is shown in this work that high concentrations of ribose-5-phosphate (acceptor substrate) compete with xylulose-5-phosphate (donor substrate), suppressing the transketolase activity (Ki = 3.8 mM). However, interacting with the donor-substrate binding site on the protein molecule, the acceptor substrate, unlike the donor substrate, does not cause any change in the active site of the enzyme. The data are interesting in terms of studying the regulatory mechanism of the transketolase activity and the structure of the enzyme-substrate complex.

The two-step mechanism of coenzyme (thiamine diphosphate, ThDP) binding with two initially identical active sites of apotransketolase has been examined with a kinetic model. Cooperativity between sites in the primary ThDP binding and in the following conformational transition has been analyzed. The only reliable difference between sites is shown to be the tenfold difference in the backward rate constants of the conformational transition; this means that the cooperative interaction between sites takes place only after termination of both steps of ThDP binding in both sites.

Recent evidence underlines a key role of mitochondrial Ca2+ fluxes in cell Ca2+ signalling. We present here a kinetic model simulating the Ca2+ fluxes generated by mitochondria during mitochondrial Ca(2+)-induced Ca2+ release (mCICR) resulting from the operation of the permeability transition pore (PTP). Our model connects the Ca2+ fluxes through the ruthenium redsensitive Ca2+ uniporter, the respiration-dependent and passive H+ fluxes, the rate of oxygen consumption, the movements of weak acids across the mitochondrial membrane, the electrical transmembrane potential (delta psi), and operation of the PTP. We find that two factors are crucial to account for the various mCICR profiles that can be observed experimentally: (i) the dependence of PTP opening and closure on matrix pH (pHi), and (ii) the relative inhibition of the respiratory rate consecutive to PTP opening. The resulting model can simulate irreversible Ca2+ efflux from mitochondria, as well as the genesis of damped or sustained Ca2+ oscillations, and of single Ca2+ spikes. The model also simulates the main features of mCICR, i.e. the threshold-dependence of mCICR triggering, and the all-or-nothing nature of mCICR operation. Our model should appear useful to further mathematically address the consequences of mCICR on the spatiotemporal organisation of Ca2+ signals, as a 'plug-in' module for the existing models of cell Ca2+ signalling.

The two-step mechanism of coenzyme (TDP) binding to apotransketolase has been examined by kinetic modeling, and the rate and equilibrium constants for each binding step for two active sites have been determined. The dissociation constants for the primary fast binding step and the forward rate constants for the secondary slow binding step have been shown to be similar for two active sites. The backward rate constants for the secondary binding step are different for two active sites, providing the kinetic mechanism of their non-equivalence in TDP binding.

A kinetic model of bisubstrate reaction catalyzed by baker's yeast transketolase is proposed. The model considers individual stages of substrates reversible primary binding. The model corresponds to the observed kinetics of product accumulation within a wide range of initial substrate concentrations. Kinetic parameters for the best simulation of the experimental data are defined. The equilibrium constants of the primary binding of both the initial and produced ketose and also the initial aldose were unequivocally determined by varying the initial substrate concentrations. The dissociation constants of the primary enzyme-substrate complex for the initial ketose (xylulose 5-phosphate) and the reaction product (sedoheptulose 7-phosphate) were found to differ by more than by two orders of magnitude. The result is discussed in the context of the hypothesis of flip-flop functioning of the transketolase active sites.

A model of myosin complex with two active sites is considered. Each of the sites is supposed to pass through four states during a cycle of ATP hydrolysis. The model describes the distribution of the complex between states and the ATP hydrolysis rate as dependent on ATP concentration and reaction constants. Some of the constants were determined experimentally [6], but these values reproduce experimental curves [5] only with the assumption of cooperativity. If the first site is bind to actin, the rate of binding for the second one is shown to increase tenfold. If one of the sites is bind to actin by a "pulling" bridge, and the second site is in "rigor" state, then the ATP binding for the second site is about ten times faster. The transition into the rigor state proceeds much faster if both sites form pulling bridges. The rate of binding to actin is the same for one and for two sites.

A kinetic model for the pyruvate dehydrogenase complex is analyzed. The model takes into account intermediate channeling through the lipoyl network attached to the complex core, as well as inter-related regulatory effects of protein X acetylation and enzyme phosphorylation. The model predicts undamped oscillations of enzyme activity.

A novel regulation type, which may be observed as an unusual kinetic 'cooperativity', is predicted for the alpha-keto acid dehydrogenase complexes. The inter-relationship of this regulation with the well-known regulatory effect of enzyme phosphorylation is discussed.