Thermodynamics-based Flux Analysis Research Articles

Growth in fluctuating light buffers plants against photorespiratory perturbations

Photorespiration (PR) is the pathway that detoxifies the product of the oxygenation reaction of Rubisco. It has been hypothesized that in dynamic light environments, PR provides a photoprotective function. To test this hypothesis, we characterized plants with varying PR enzyme activities under fluctuating and non-fluctuating light conditions. Contrasting our expectations, growth of mutants with decreased PR enzyme levels was least affected in fluctuating light compared with wild type. Results for growth, photosynthesis and metabolites combined with thermodynamics-based flux analysis revealed two main causal factors for this unanticipated finding: reduced rates of photosynthesis in fluctuating light and complex re-routing of metabolic fluxes. Only in non-fluctuating light, mutants lacking the glutamate:glyoxylate aminotransferase 1 re-routed glycolate processing to the chloroplast, resulting in photooxidative damage through H2O2 production. Our results reveal that dynamic light environments buffer plant growth and metabolism against photorespiratory perturbations.

Metabolic control analysis of L-tryptophan producing Escherichia coli applying targeted perturbation with shikimate

L-tryptophan production from glycerol with Escherichia coli was analysed by perturbation studies and metabolic control analysis. The insertion of a non-natural shikimate transporter into the genome of an Escherichia coli L-tryptophan production strain enabled targeted perturbation within the product pathway with shikimate during parallelised short-term perturbation experiments with cells withdrawn from a 15 L fed-batch production process. Expression of the shikimate/H+-symporter gene (shiA) from Corynebacterium glutamicum did not alter process performance within the estimation error. Metabolic analyses and subsequent extensive data evaluation were performed based on the data of the parallel analysis reactors and the production process. Extracellular rates and intracellular metabolite concentrations displayed evident deflections in cell metabolism and particularly in chorismate biosynthesis due to the perturbations with shikimate. Intracellular flux distributions were estimated using a thermodynamics-based flux analysis method, which integrates thermodynamic constraints and intracellular metabolite concentrations to restrain the solution space. Feasible flux distributions, Gibbs reaction energies and concentration ranges were computed simultaneously for the genome-wide metabolic model, with minimum bias in relation to the direction of metabolic reactions. Metabolic control analysis was applied to estimate elasticities and flux control coefficients, predicting controlling sites for L-tryptophan biosynthesis. The addition of shikimate led to enhanced deviations in chorismate biosynthesis, revealing a so far not observed control of 3-dehydroquinate synthase on L-tryptophan formation. The relative expression of the identified target genes was analysed with RT-qPCR. Transcriptome analysis revealed disparities in gene expression and the localisation of target genes to further improve the microbial L-tryptophan producer by metabolic engineering.

MultiTFA: a Python package for multi-variate thermodynamics-based flux analysis.

MotivationWe achieve a significant improvement in thermodynamic-based flux analysis (TFA) by introducing multivariate treatment of thermodynamic variables and leveraging component contribution, the state-of-the-art implementation of the group contribution methodology. Overall, the method greatly reduces the uncertainty of thermodynamic variables.ResultsWe present multiTFA, a Python implementation of our framework. We evaluated our application using the core Escherichia coli model and achieved a median reduction of 6.8 kJ/mol in reaction Gibbs free energy ranges, while three out of 12 reactions in glycolysis changed from reversible to irreversible.Availability and implementationOur framework along with documentation is available on https://github.com/biosustain/multitfa.Supplementary informationSupplementary data are available at Bioinformatics online.

Physicochemical and metabolic constraints for thermodynamics-based stoichiometric modelling under mesophilic growth conditions.

Metabolic engineering in the post-genomic era is characterised by the development of new methods for metabolomics and fluxomics, supported by the integration of genetic engineering tools and mathematical modelling. Particularly, constraint-based stoichiometric models have been widely studied: (i) flux balance analysis (FBA) (in silico), and (ii) metabolic flux analysis (MFA) (in vivo). Recent studies have enabled the incorporation of thermodynamics and metabolomics data to improve the predictive capabilities of these approaches. However, an in-depth comparison and evaluation of these methods is lacking. This study presents a thorough analysis of two different in silico methods tested against experimental data (metabolomics and 13C-MFA) for the mesophile Escherichia coli. In particular, a modified version of the recently published matTFA toolbox was created, providing a broader range of physicochemical parameters. Validating against experimental data allowed the determination of the best physicochemical parameters to perform the TFA (Thermodynamics-based Flux Analysis). An analysis of flux pattern changes in the central carbon metabolism between 13C-MFA and TFA highlighted the limited capabilities of both approaches for elucidating the anaplerotic fluxes. In addition, a method based on centrality measures was suggested to identify important metabolites that (if quantified) would allow to further constrain the TFA. Finally, this study emphasised the need for standardisation in the fluxomics community: novel approaches are frequently released but a thorough comparison with currently accepted methods is not always performed.

Large-scale kinetic metabolic models of Pseudomonas putida KT2440 for consistent design of metabolic engineering strategies

BackgroundPseudomonas putida is a promising candidate for the industrial production of biofuels and biochemicals because of its high tolerance to toxic compounds and its ability to grow on a wide variety of substrates. Engineering this organism for improved performances and predicting metabolic responses upon genetic perturbations requires reliable descriptions of its metabolism in the form of stoichiometric and kinetic models.ResultsIn this work, we developed kinetic models of P. putida to predict the metabolic phenotypes and design metabolic engineering interventions for the production of biochemicals. The developed kinetic models contain 775 reactions and 245 metabolites. Furthermore, we introduce here a novel set of constraints within thermodynamics-based flux analysis that allow for considering concentrations of metabolites that exist in several compartments as separate entities. We started by a gap-filling and thermodynamic curation of iJN1411, the genome-scale model of P. putida KT2440. We then systematically reduced the curated iJN1411 model, and we created three core stoichiometric models of different complexity that describe the central carbon metabolism of P. putida. Using the medium complexity core model as a scaffold, we generated populations of large-scale kinetic models for two studies. In the first study, the developed kinetic models successfully captured the experimentally observed metabolic responses to several single-gene knockouts of a wild-type strain of P. putida KT2440 growing on glucose. In the second study, we used the developed models to propose metabolic engineering interventions for improved robustness of this organism to the stress condition of increased ATP demand.ConclusionsThe study demonstrates the potential and predictive capabilities of the kinetic models that allow for rational design and optimization of recombinant P. putida strains for improved production of biofuels and biochemicals. The curated genome-scale model of P. putida together with the developed large-scale stoichiometric and kinetic models represents a significant resource for researchers in industry and academia.

PyTFA and matTFA: a Python package and a Matlab toolbox for Thermodynamics-based Flux Analysis.

SummarypyTFA and matTFA are the first published implementations of the original TFA paper. Specifically, they include explicit formulation of Gibbs energies and metabolite concentrations, which enables straightforward integration of metabolite concentration measurements.MotivationHigh-throughput analytic technologies provide a wealth of omics data that can be used to perform thorough analyses for a multitude of studies in the areas of Systems Biology and Biotechnology. Nevertheless, most studies are still limited to constraint-based Flux Balance Analyses (FBA), neglecting an important physicochemical constraint: thermodynamics. Thermodynamics-based Flux Analysis (TFA) in metabolic models enables the integration of quantitative metabolomics data to study their effects on the net-flux directionality of reactions in the network. In addition, it allows us to estimate how far each reaction operates from thermodynamic equilibrium, which provides critical information for guiding metabolic engineering decisions.ResultsWe present a Python package (pyTFA) and a Matlab toolbox (matTFA) that implement TFA. We show an example of application on both a reduced and a genome-scale model of E. coli., and demonstrate TFA and data integration through TFA reduce the feasible flux space with respect to FBA.Availability and implementationDocumented implementation of TFA framework both in Python (pyTFA) and Matlab (matTFA) are available on www.github.com/EPFL-LCSB/.Supplementary information Supplementary data are available at Bioinformatics online.

Bioenergetics-based modeling of Plasmodium falciparum metabolism reveals its essential genes, nutritional requirements, and thermodynamic bottlenecks

Novel antimalarial therapies are urgently needed for the fight against drug-resistant parasites. The metabolism of malaria parasites in infected cells is an attractive source of drug targets but is rather complex. Computational methods can handle this complexity and allow integrative analyses of cell metabolism. In this study, we present a genome-scale metabolic model (iPfa) of the deadliest malaria parasite, Plasmodium falciparum, and its thermodynamics-based flux analysis (TFA). Using previous absolute concentration data of the intraerythrocytic parasite, we applied TFA to iPfa and predicted up to 63 essential genes and 26 essential pairs of genes. Of the 63 genes, 35 have been experimentally validated and reported in the literature, and 28 have not been experimentally tested and include previously hypothesized or novel predictions of essential metabolic capabilities. Without metabolomics data, four of the genes would have been incorrectly predicted to be non-essential. TFA also indicated that substrate channeling should exist in two metabolic pathways to ensure the thermodynamic feasibility of the flux. Finally, analysis of the metabolic capabilities of P. falciparum led to the identification of both the minimal nutritional requirements and the genes that can become indispensable upon substrate inaccessibility. This model provides novel insight into the metabolic needs and capabilities of the malaria parasite and highlights metabolites and pathways that should be measured and characterized to identify potential thermodynamic bottlenecks and substrate channeling. The hypotheses presented seek to guide experimental studies to facilitate a better understanding of the parasite metabolism and the identification of targets for more efficient intervention.

Thermodynamics-based Metabolite Sensitivity Analysis in metabolic networks

The increasing availability of large metabolomics datasets enhances the need for computational methodologies that can organize the data in a way that can lead to the inference of meaningful relationships. Knowledge of the metabolic state of a cell and how it responds to various stimuli and extracellular conditions can offer significant insight in the regulatory functions and how to manipulate them. Constraint based methods, such as Flux Balance Analysis (FBA) and Thermodynamics-based flux analysis (TFA), are commonly used to estimate the flow of metabolites through genome-wide metabolic networks, making it possible to identify the ranges of flux values that are consistent with the studied physiological and thermodynamic conditions. However, unless key intracellular fluxes and metabolite concentrations are known, constraint-based models lead to underdetermined problem formulations. This lack of information propagates as uncertainty in the estimation of fluxes and basic reaction properties such as the determination of reaction directionalities. Therefore, knowledge of which metabolites, if measured, would contribute the most to reducing this uncertainty can significantly improve our ability to define the internal state of the cell. In the present work we combine constraint based modeling, Design of Experiments (DoE) and Global Sensitivity Analysis (GSA) into the Thermodynamics-based Metabolite Sensitivity Analysis (TMSA) method. TMSA ranks metabolites comprising a metabolic network based on their ability to constrain the gamut of possible solutions to a limited, thermodynamically consistent set of internal states. TMSA is modular and can be applied to a single reaction, a metabolic pathway or an entire metabolic network. This is, to our knowledge, the first attempt to use metabolic modeling in order to provide a significance ranking of metabolites to guide experimental measurements.