Quantifying Dynamic Regulation in Metabolic Pathways with Nonparametric Flux Inference.
Ontology highlight
ABSTRACT:
SUBMITTER: He F
PROVIDER: S-EPMC6531928 | biostudies-literature | 2019 May
REPOSITORIES: biostudies-literature
Publications
Quantifying Dynamic Regulation in Metabolic Pathways with Nonparametric Flux Inference.
He Fei F Stumpf Michael P H MPH
Biophysical journal 20190419 10
One of the central tasks in systems biology is to understand how cells regulate their metabolism. Hierarchical regulation analysis is a powerful tool to study this regulation at the metabolic, gene-expression, and signaling levels. It has been widely applied to study steady-state regulation, but analysis of the metabolic dynamics remains challenging because it is difficult to measure time-dependent metabolic flux. Here, we develop a nonparametric method that uses Gaussian processes to accurately ...[more]
Similar Datasets
Project description:Most metabolic pathways contain more reactions than metabolites and therefore have a wide stoichiometric matrix that corresponds to infinitely many possible flux distributions that are perfectly compatible with the dynamics of the metabolites in a given dataset. This under-determinedness poses a challenge for the quantitative characterization of flux distributions from time series data and thus for the design of adequate, predictive models. Here we propose a method that reduces the degrees of freedom in a stepwise manner and leads to a dynamic flux distribution that is, in a statistical sense, likely to be close to the true distribution.We applied the proposed method to the lignin biosynthesis pathway in switchgrass. The system consists of 16 metabolites and 23 enzymatic reactions. It has seven degrees of freedom and therefore admits a large space of dynamic flux distributions that all fit a set of metabolic time series data equally well. The proposed method reduces this space in a systematic and biologically reasonable manner and converges to a likely dynamic flux distribution in just a few iterations. The estimated solution and the true flux distribution, which is known in this case, show excellent agreement and thereby lend support to the method.The computational model was implemented in MATLAB (version R2014a, The MathWorks, Natick, MA). The source code is available at https://github.gatech.edu/VoitLab/Stepwise-Inference-of-Likely-Dynamic-Flux-Distributions and www.bst.bme.gatech.edu/research.php .mojdeh@gatech.edu or eberhard.voit@bme.gatech.edu.Supplementary data are available at Bioinformatics online.
| S-EPMC5860468 | biostudies-other
Project description:For an industrial fermentation process, it can be advantageous to decouple cell growth from product formation. This decoupling would allow for the rapid accumulation of biomass without inhibition from product formation, after which the fermentation can be switched to a mode where cells would grow minimally and primarily act as catalysts to convert substrate into desired product. The switch in fermentation mode should preferably be accomplished without the addition of expensive inducers. A common cell factory Saccharomyces cerevisiae is a Crabtree-positive yeast and is typically fermented at industrial scale under glucose-limited conditions to avoid the formation of ethanol. In this work, we aimed to identify and characterize promoters that depend on glucose concentration for use as dynamic control elements. Through analysis of mRNA data of S. cerevisiae grown in chemostats under glucose excess or limitation, we identified 34 candidate promoters that strongly responded to glucose presence or absence. These promoters were characterized in small-scale batch and fed-batch cultivations using a quickly maturing rapidly degrading green fluorescent protein yEGFP3-Cln2PEST as a reporter. Expressing 3-hydroxypropionic acid (3HP) pathway from a set of selected regulated promoters allowed for suppression of 3HP production during glucose-excess phase of a batch cultivation with subsequent activation in glucose-limiting conditions. Regulating the 3HP pathway by the ICL1 promoter resulted in 70% improvement of 3HP titer in comparison to PGK1 promoter.
| S-EPMC5972318 | biostudies-literature
Project description:Metabolic engineering of microorganisms to produce desirable products on an industrial scale can result in unbalanced cellular metabolic networks that reduce productivity and yield. Metabolic fluxes can be rebalanced using dynamic pathway regulation, but few broadly applicable tools are available to achieve this. We present a pathway-independent genetic control module that can be used to dynamically regulate the expression of target genes. We apply our module to identify the optimal point to redirect glycolytic flux into heterologous engineered pathways in Escherichia coli, resulting in titers of myo-inositol increased 5.5-fold and titers of glucaric acid increased from unmeasurable to >0.8 g/L, compared to the parent strains lacking dynamic flux control. Scaled-up production of these strains in benchtop bioreactors resulted in almost ten- and fivefold increases in specific titers of myo-inositol and glucaric acid, respectively. We also used our module to control flux into aromatic amino acid biosynthesis to increase titers of shikimate in E. coli from unmeasurable to >100 mg/L.
| S-EPMC5340623 | biostudies-literature
Project description:In this article, the existing concept of reversed percentile residual life, or percentile inactivity time, is recast to show that it can be used for routine analysis of time-to-event data under right censoring to summarize "life lost," which poses several advantages over the existing methods for survival analysis. An estimating equation approach is adopted to avoid estimation of the probability density function of the underlying time-to-event distribution to estimate the variance of the quantile estimator. Additionally a K-sample test statistic is proposed to test the ratio of the quantile lost lifespans. Simulation studies are performed to assess finite properties of the proposed K-sample statistic in terms of coverage probability and power. The proposed method is illustrated with a real data example from a breast cancer study.
| S-EPMC5724980 | biostudies-literature
Project description:BackgroundThe cellular signaling pathway (network) is one of the main topics of organismic investigations. The intracellular interactions between genes in a signaling pathway are considered as the foundation of functional genomics. Thus, what genes and how much they influence each other through transcriptional binding or physical interactions are essential problems. Under the synchronous measures of gene expression via a microarray chip, an amount of dynamic information is embedded and remains to be discovered. Using a systematically dynamic modeling approach, we explore the causal relationship among genes in cellular signaling pathways from the system biology approach.ResultsIn this study, a second-order dynamic model is developed to describe the regulatory mechanism of a target gene from the upstream causality point of view. From the expression profile and dynamic model of a target gene, we can estimate its upstream regulatory function. According to this upstream regulatory function, we would deduce the upstream regulatory genes with their regulatory abilities and activation delays, and then link up a regulatory pathway. Iteratively, these regulatory genes are considered as target genes to trace back their upstream regulatory genes. Then we could construct the regulatory pathway (or network) to the genome wide. In short, we can infer the genetic regulatory pathways from gene-expression profiles quantitatively, which can confirm some doubted paths or seek some unknown paths in a regulatory pathway (network). Finally, the proposed approach is validated by randomly reshuffling the time order of microarray data.ConclusionWe focus our algorithm on the inference of regulatory abilities of the identified causal genes, and how much delay before they regulate the downstream genes. With this information, a regulatory pathway would be built up using microarray data. In the present study, two signaling pathways, i.e. circadian regulatory pathway in Arabidopsis thaliana and metabolic shift pathway from fermentation to respiration in yeast Saccharomyces cerevisiae, are reconstructed using microarray data to evaluate the performance of our proposed method. In the circadian regulatory pathway, we identified mainly the interactions between the biological clock and the photoperiodic genes consistent with the known regulatory mechanisms. We also discovered the now less-known regulations between crytochrome and phytochrome. In the metabolic shift pathway, the casual relationship of enzymatic genes could be detected properly.
| S-EPMC555938 | biostudies-literature
Project description:Recently, approaches have been developed to sample the genetic content of heterogeneous environments (metagenomics). However, by what means these sequences link distinct environmental conditions with specific biological processes is not well understood. Thus, a major challenge is how the usage of particular pathways and subnetworks reflects the adaptation of microbial communities across environments and habitats-i.e., how network dynamics relates to environmental features. Previous research has treated environments as discrete, somewhat simplified classes (e.g., terrestrial vs. marine), and searched for obvious metabolic differences among them (i.e., treating the analysis as a typical classification problem). However, environmental differences result from combinations of many factors, which often vary only slightly. Therefore, we introduce an approach that employs correlation and regression to relate multiple, continuously varying factors defining an environment to the extent of particular microbial pathways present in a geographic site. Moreover, rather than looking only at individual correlations (one-to-one), we adapted canonical correlation analysis and related techniques to define an ensemble of weighted pathways that maximally covaries with a combination of environmental variables (many-to-many), which we term a metabolic footprint. Applied to available aquatic datasets, we identified footprints predictive of their environment that can potentially be used as biosensors. For example, we show a strong multivariate correlation between the energy-conversion strategies of a community and multiple environmental gradients (e.g., temperature). Moreover, we identified covariation in amino acid transport and cofactor synthesis, suggesting that limiting amounts of cofactor can (partially) explain increased import of amino acids in nutrient-limited conditions.
| S-EPMC2629784 | biostudies-literature
Project description:(+)-Catechin is an important antioxidant of green tea (Camelia sinensis (L.) O. Kuntze). Catechin is known for its positive role in anticancerous activity, extracellular matrix degradation, cell death regulation, diabetes, and other related disorders. As a result of enormous interest in and great demand for catechin, its biosynthesis using metabolic engineering has become the subject of concentrated research with the aim of enhancing (+)-catechin production. Metabolic flux is an essential concept in the practice of metabolic engineering as it helps in the identification of the regulatory element of a biosynthetic pathway. In the present study, an attempt was made to analyze the metabolic flux of the (+)-catechin biosynthesis pathway in order to decipher the regulatory element of this pathway. Firstly, a genetically encoded fluorescence resonance energy transfer (FRET)-based nanosensor (FLIP-Cat, fluorescence indicator protein for (+)-catechin) was developed for real-time monitoring of (+)-catechin flux. In vitro characterization of the purified protein of the nanosensor showed that the nanosensor was pH stable and (+)-catechin specific. Its calculated Kd was 139 µM. The nanosensor also performed real-time monitoring of (+)-catechin in bacterial cells. In the second step of this study, an entire (+)-catechin biosynthesis pathway was constructed and expressed in E. coli in two sets of plasmid constructs: pET26b-PT7-rbs-PAL-PT7-rbs-4CL-PT7-rbs-CHS-PT7-rbs-CHI and pET26b-T7-rbs-F3H-PT7-rbs- DFR-PT7-rbs-LCR. The E. coli harboring the FLIP-Cat was transformed with these plasmid constructs. The metabolic flux analysis of (+)-catechin was carried out using the FLIP-Cat. The FLIP-Cat successfully monitored the flux of catechin after adding tyrosine, 4-coumaric acid, 4-coumaroyl CoA, naringenin chalcone, naringenin, dihydroquercetin, and leucocyanidin, individually, with the bacterial cells expressing the nanosensor as well as the genes of the (+)-catechin biosynthesis pathway. Dihydroflavonol reductase (DFR) was identified as the main regulatory element of the (+)-catechin biosynthesis pathway. Information about this regulatory element of the (+)-catechin biosynthesis pathway can be used for manipulating the (+)-catechin biosynthesis pathway using a metabolic engineering approach to enhance production of (+)-catechin.
| S-EPMC7222200 | biostudies-literature
Project description:BackgroundAdvances in modern high-throughput techniques of molecular biology have enabled top-down approaches for the estimation of parameter values in metabolic systems, based on time series data. Special among them is the recent method of dynamic flux estimation (DFE), which uses such data not only for parameter estimation but also for the identification of functional forms of the processes governing a metabolic system. DFE furthermore provides diagnostic tools for the evaluation of model validity and of the quality of a model fit beyond residual errors. Unfortunately, DFE works only when the data are more or less complete and the system contains as many independent fluxes as metabolites. These drawbacks may be ameliorated with other types of estimation and information. However, such supplementations incur their own limitations. In particular, assumptions must be made regarding the functional forms of some processes and detailed kinetic information must be available, in addition to the time series data.ResultsThe authors propose here a systematic approach that supplements DFE and overcomes some of its shortcomings. Like DFE, the approach is model-free and requires only minimal assumptions. If sufficient time series data are available, the approach allows the determination of a subset of fluxes that enables the subsequent applicability of DFE to the rest of the flux system. The authors demonstrate the procedure with three artificial pathway systems exhibiting distinct characteristics and with actual data of the trehalose pathway in Saccharomyces cerevisiae.ConclusionsThe results demonstrate that the proposed method successfully complements DFE under various situations and without a priori assumptions regarding the model representation. The proposed method also permits an examination of whether at all, to what degree, or within what range the available time series data can be validly represented in a particular functional format of a flux within a pathway system. Based on these results, further experiments may be designed to generate data points that genuinely add new information to the structure identification and parameter estimation tasks at hand.
| S-EPMC3495652 | biostudies-literature
Project description:To quantify gene regulation, a function is required that relates transcription factor binding to DNA (input) to the rate of mRNA synthesis from a target gene (output). Such a 'gene regulation function' (GRF) generally cannot be measured because the experimental titration of inputs and simultaneous readout of outputs is difficult. Here we show that GRFs may instead be inferred from natural changes in cellular gene expression, as exemplified for the cell cycle in the yeast S. cerevisiae. We develop this inference approach based on a time series of mRNA synthesis rates from a synchronized population of cells observed over three cell cycles. We first estimate the functional form of how input transcription factors determine mRNA output and then derive GRFs for target genes in the CLB2 gene cluster that are expressed during G2/M phase. Systematic analysis of additional GRFs suggests a network architecture that rationalizes transcriptional cell cycle oscillations. We find that a transcription factor network alone can produce oscillations in mRNA expression, but that additional input from cyclin oscillations is required to arrive at the native behaviour of the cell cycle oscillator.
| S-EPMC5072840 | biostudies-literature
Project description:The ability to measure the properties of proteins at the single-molecule level offers an unparalleled glimpse into biological systems at the molecular scale. The interpretation of single-molecule time series has often been rooted in statistical mechanics and the theory of Markov processes. While existing analysis methods have been useful, they are not without significant limitations including problems of model selection and parameter nonidentifiability. To address these challenges, we introduce the use of nonparametric Bayesian inference for the analysis of single-molecule time series. These methods provide a flexible way to extract structure from data instead of assuming models beforehand. We demonstrate these methods with applications to several diverse settings in single-molecule biophysics. This approach provides a well-constrained and rigorously grounded method for determining the number of biophysical states underlying single-molecule data.
| S-EPMC4317543 | biostudies-literature