Basic concepts and principles of stoichiometric modeling of metabolic networks.
Ontology highlight
ABSTRACT:
SUBMITTER: Maarleveld TR
PROVIDER: S-EPMC4671265 | biostudies-literature | 2013 Sep
REPOSITORIES: biostudies-literature
Publications
Basic concepts and principles of stoichiometric modeling of metabolic networks.
Maarleveld Timo R TR Khandelwal Ruchir A RA Olivier Brett G BG Teusink Bas B Bruggeman Frank J FJ
Biotechnology journal 20130729 9
Metabolic networks supply the energy and building blocks for cell growth and maintenance. Cells continuously rewire their metabolic networks in response to changes in environmental conditions to sustain fitness. Studies of the systemic properties of metabolic networks give insight into metabolic plasticity and robustness, and the ability of organisms to cope with different environments. Constraint-based stoichiometric modeling of metabolic networks has become an indispensable tool for such studi ...[more]
Similar Datasets
Project description:Recent advances in metabolomics technologies have resulted in high-quality (time-resolved) metabolic profiles with an increasing coverage of metabolic pathways. These data profiles represent read-outs from often non-linear dynamics of metabolic networks. Yet, metabolic profiles have largely been explored with regression-based approaches that only capture linear relationships, rendering it difficult to determine the extent to which the data reflect the underlying reaction rates and their couplings. Here we propose an approach termed Stoichiometric Correlation Analysis (SCA) based on correlation between positive linear combinations of log-transformed metabolic profiles. The log-transformation is due to the evidence that metabolic networks can be modeled by mass action law and kinetics derived from it. Unlike the existing approaches which establish a relation between pairs of metabolites, SCA facilitates the discovery of higher-order dependence between more than two metabolites. By using a paradigmatic model of the tricarboxylic acid cycle we show that the higher-order dependence reflects the coupling of concentration of reactant complexes, capturing the subtle difference between the employed enzyme kinetics. Using time-resolved metabolic profiles from Arabidopsis thaliana and Escherichia coli, we show that SCA can be used to quantify the difference in coupling of reactant complexes, and hence, reaction rates, underlying the stringent response in these model organisms. By using SCA with data from natural variation of wild and domesticated wheat and tomato accession, we demonstrate that the domestication is accompanied by loss of such couplings, in these species. Therefore, application of SCA to metabolomics data from natural variation in wild and domesticated populations provides a mechanistic way to understanding domestication and its relation to metabolic networks.
| S-EPMC5741659 | biostudies-literature
Project description:Modeling is an important tool in drug development; population modeling is a complex process requiring robust underlying procedures for ensuring clean data, appropriate computing platforms, adequate resources, and effective communication. Although requiring an investment in resources, it can save time and money by providing a platform for integrating all information gathered on new therapeutic agents. This article provides a brief overview of aspects of modeling and simulation as applied to many areas in drug development.CPT: Pharmacometrics & Systems Pharmacology (2012) 1, e6; doi:10.1038/psp.2012.4; advance online publication 26 September 2012.
| S-EPMC3606044 | biostudies-other
Project description:The dynamic features of a genetic network's response to environmental fluctuations represent essential functional specifications and thus may constrain the possible choices of network architecture and kinetic parameters. To explore the connection between dynamics and network design, we have analyzed a general regulatory architecture that is commonly found in many metabolic pathways. Such architecture is characterized by a dual control mechanism, with end product feedback inhibition and transcriptional regulation mediated by an intermediate metabolite. As a case study, we measured with high temporal resolution the induction profiles of the enzymes in the leucine biosynthetic pathway in response to leucine depletion, using an automated system for monitoring protein expression levels in single cells. All the genes in the pathway are known to be coregulated by the same transcription factors, but we observed drastically different dynamic responses for enzymes upstream and immediately downstream of the key control point-the intermediate metabolite alpha-isopropylmalate (alphaIPM), which couples metabolic activity to transcriptional regulation. Analysis based on genetic perturbations suggests that the observed dynamics are due to differential regulation by the leucine branch-specific transcription factor Leu3, and that the downstream enzymes are strictly controlled and highly expressed only when alphaIPM is available. These observations allow us to build a simplified mathematical model that accounts for the observed dynamics and can correctly predict the pathway's response to new perturbations. Our model also suggests that transient dynamics and steady state can be separately tuned and that the high induction levels of the downstream enzymes are necessary for fast leucine recovery. It is likely that principles emerging from this work can reveal how gene regulation has evolved to optimize performance in other metabolic pathways with similar architecture.
| S-EPMC2429954 | biostudies-literature
Project description:Uncovering the general principles that govern the structure of metabolic networks is key to understanding the emergence and evolution of living systems. Artificial chemistries can help illuminate this problem by enabling the exploration of chemical reaction universes that are constrained by general mathematical rules. Here, we focus on artificial chemistries in which strings of characters represent simplified molecules, and string concatenation and splitting represent possible chemical reactions. We developed a novel Python package, ARtificial CHemistry NEtwork Toolbox (ARCHNET), to study string chemistries using tools from the field of stoichiometric constraint-based modeling. In addition to exploring the topological characteristics of different string chemistry networks, we developed a network-pruning algorithm that can generate minimal metabolic networks capable of producing a specified set of biomass precursors from a given assortment of environmental nutrients. We found that the composition of these minimal metabolic networks was influenced more strongly by the metabolites in the biomass reaction than the identities of the environmental nutrients. This finding has important implications for the reconstruction of organismal metabolic networks and could help us better understand the rise and evolution of biochemical organization. More generally, our work provides a bridge between artificial chemistries and stoichiometric modeling, which can help address a broad range of open questions, from the spontaneous emergence of an organized metabolism to the structure of microbial communities.
| S-EPMC8318951 | biostudies-literature
Project description:Population pharmacodynamic (PD) models describe the time course of drug effects, relating exposure to response, and providing a more robust understanding of drug action than single assessments. PD models can test alternative dose regimens through simulation, allowing for informed assessment of potential dose regimens and study designs. This is the third paper in a three-part series, providing an introduction into methods for developing and evaluating population PD models. Example files are available in the Supplementary Data.
| S-EPMC3917320 | biostudies-literature
Project description:Using laser capture microdissection to isolate over 100 specific cell populations, we report the profiling of prostate cancer progression from benign epithelium to metastatic disease. By analyzing these expression signatures in the context of over 15,000 “molecular concepts”, or sets of biologically connected genes, we generated an integrative model of prostate cancer progression. Keywords: disease state analysis
2006-12-17 | GSE6099 | GEO
Project description:Using laser capture microdissection to isolate over 100 specific cell populations, we report the profiling of prostate cancer progression from benign epithelium to metastatic disease. By analyzing these expression signatures in the context of over 15,000 â??molecular conceptsâ??, or sets of biologically connected genes, we generated an integrative model of prostate cancer progression. Keywords: disease state analysis Experiment Design: â?¢ Type of experiment: This study used microarray experiments to profile prostate cancer progression through the isolation of pure cell populations using laser capture microdissection (LCM) and OmniPlex Whole Transcriptome (WTA) Amplification. â?¢ Experimental factors: This study profiled various pure cell populations isolated by LCM from prostate tissue. â?¢ The number of hybridizations performed in the experiment: A total of 104 hybridizations were performed. â?¢ The type of reference used for the hybridizations: A single commercially available pool of benign prostate tissue (CPP, Clontech) was used as the reference in all hybridizations. For the reference sample, 10 ng of CPP total RNA was converted into a WTA cDNA library in a volume of 25 µl and five 5 µl aliquots were amplified by WTA PCR and purified. Five ng aliquots of the pooled products were subjected to a second WTA PCR amplification in the presence of amino-allyl dUTP and all products were pooled before labeling to ensure equal representation across all of the hybridizations. â?¢ Hybridization design: For all hybridizations, the experimental prostate sample isolated by LCM was labeled with Cy5 and the common reference of benign pooled prostate as described above was labeled with Cy3. â?¢ Quality control steps taken: In order to assess biological and technical reproducibility, we captured identical regions from two specimens (PCA_11 and MET_HR_3) and subjected one sample, PCA_30, to two hybridizations. For PCA_11, adjacent regions of the same section were captured and treated as separate samples. For MET_HR_3, two separate foci of the liver metastasis (from separate frozen blocks) were captured and treated as separate samples. For PCA_30, two aliquots of the primary WTA PCR products were separately subjected to a second round of WTA PCR amplification. Replicates from PCA_30 and MET_HR_3 were also hybridized to arrays from each of the print runs utilized in this study. Samples used, extract preparation and labeling: â?¢ The origin of the biological sample (for instance, name of the organism, the provider of the sample) and its characteristics: All samples used were of human origin. The table in the MIAME Checklist (Supplmentary Methods) summarizes the samples used. Hybridizations are named according to the class of sample profiled and the background color is the same as shown in Figure S3. Two print runs were used during the course of the study and the indicated run for each hybridization is given. For prostate cancer samples, the Gleason patterns present on the sample section used are indicated, as well as the Gleason patterns of the actual cells isolated by LCM. A unique case number for all profiled samples is given. Descriptions of samples are given where applicable, including the location of metastatic foci. â?¢ Manipulation of biological samples and protocols used: Laser Capture Microdissection (LCM) was performed from frozen tissue sections with the SL Microtest device using µCUT software (MMI). Approximately 10,000 cells were captured for each sample. Serial sections were used if cells could not be obtained from a single section. Total RNA was isolated from captured cells with the RNAqueous Micro kit (Ambion) and treated with DNAse I according to the manufacturer's instructions. RNA quantification was perfomed using Ribogreen (Molecular Probes). â?¢ Protocol for preparing the hybridization extract: For each sample, total RNA (~10ng) in a volume of 25 µl was converted into an OmniPlex WTA cDNA library and amplified by WTA PCR using reagents and protocols according to the beta-commercial TransPlex WTA kit (Rubicon Genomics, Ann Arbor, MI). For each sample, a single 10 µl aliquot of the WTA cDNA library was amplified by WTA PCR and products were purified. Four or five 5 ng aliquots of product were subjected to a second WTA PCR amplification in the presence of amino-allyl dUTP for post amplification labeling and products were pooled before proceeding with the hybridization. Yields after all WTA PCR amplifications were between 2 to 5 μg per reaction, and reaction progress was monitored in real time using SYBR green I (Molecular Probes) on the iCycler IQ (BioRad) or ABI 7300 Real Time PCR system (Applied Biosystems). Real-time PCR amplifications were terminated at or before plateau phase as measured by fluorescence incorporation to preserve maximum representation. â?¢ Labeling protocol(s): For all WTA amplified cDNA samples with amino-allyl dUTP incorporated (10-20 ï?g), nuclease free water was added to 500 ï?l, the sample was transferred to a Microcon YM-30 filter, and centrifuged at 13,000 rcf at RT for 12 minutes. The spin column was inverted in a new collection tube and centrifuged at 13,000 rcf for 2 min at RT. The amount of sample was measured and the volume was adjusted to 5 ï?l with nuclease free water. For labeling, 5 ï?l of 1 M sodium bicarbonate (pH 9.0) was added and allowed to incubate at RT for 15 min. Nine ï?l of anhydrous DMSO were added to Mono Reactive Cy3 and Cy5 dye packs (Amersham, Buckinghamshire, England), mixed thoroughly, and 1 ï?l of the proper dye was added to each sample. The labeling mixture was incubated in the dark at RT for 1 hour. The reaction was stopped by the addition of 9 ï?l of 4 M Hydroxylamine for 15 min. For each labeling mixture, RNase free water was added to 100 ï?l, 500 ï?l of PB buffer was added and Cy3 and Cy5 mixtures were added to separate Qiaquick (Qiagen, Valencia, CA) spin columns and centrifuged at 13,000 rcf. Columns were washed twice with 700 ï?l of buffer PE, and centrifuged dry at 13,000 rcf for 2 min. To elute, 60 ï?l of EB buffer was placed on the column, incubated for 5 min and centrifuged at 13,000 rcf for 1 min. The Cy3 labeled and Cy5 labeled cDNA for each hybridization were mixed and 1.5 ï?l were used for analysis on a ND-1000 spectrophotometer. â?¢ External controls (spikes): Not utilized Hybridization procedures and parameters: â?¢ The protocol and conditions used during hybridization, blocking and washing: For hybridization, 40 ï?g of human Cot-1 DNA (Invitrogen, Carlsbad, CA) was added to 120 ï?l of labeled cDNA, the probe was transferred to a Microcon YM-30 filter and centrifuged at 13,000 rcf for 6 min at RT. The spin column was inverted in a new collection tube and centrifuged at 13,000 rcf for 2 min at RT. The eluted probe volume was adjusted to 18.6 ï?l with nuclease free water and the following were added: 4 ï?l of yeast tRNA (10 ï?g/ï?l, Invitrogen, Carlsbad, CA), 4.9 ï?l of 20X SSC and 0.84 ï?l of 10% SDS. The probe was denatured at 100 Co for 3 min and centrifuged at 13,000 rcf for 45 sec. The probe was added directly to the microarray and a cover slip was added. Slides were hybridized overnight in DIE-TECH 1 or 5 slide hybridization chambers in a 65 Co water bath. After hybridization, cover slips were removed by incubation in 2x SSC / 0.05% SDS. Slides were then washed in 2x SSC / 0.05% SDS for 5 min and 0.2x SSC / 0.05% SDS for five minutes. Slides were then washed in 0.2x SSC for 10 sec and centrifuged dry at 500 rpm for 5 min. Measurement data and specifications: â?¢ Type of scanning hardware and software used: Microarrays were scanned using a GenePix 4000B scanner (Axon Instruments, Union City, CA). â?¢ Type of image analysis software used: Images were gridded and spots were quantified using GenePix Pro 4.0 software (Axon Instruments, Union City, CA). â?¢ A description of the measurements produced by the image-analysis software and a description of which measurements were used in the analysis: The measurements made by GenePix Pro 4.0 are described in the sample files. For all hybridizations, the log2 normalized Median of Ratios (as described below) was used. â?¢ Data selection and transformation procedures: Arrays were auto-gridded by GenePix 4.0. Features flagged by GenePix as not found during grid
2010-06-05 | E-GEOD-6099 | biostudies-arrayexpress
Project description:Gene regulation is intimately connected with metabolism, enabling the appropriate timing and tuning of biochemical pathways to substrate availability. In microorganisms, such as archaea and bacteria, transcription factors (TFs) often directly sense external cues such as nutrient substrates, metabolic intermediates, or redox status to regulate gene expression. Intense recent interest has characterized the functions of a large number of such regulatory TFs in archaea, which regulate a diverse array of unique archaeal metabolic capabilities. However, it remains unclear how the co-ordinated activity of the interconnected metabolic and transcription networks produces the dynamic flexibility so frequently observed in archaeal cells as they respond to energy limitation and intermittent substrate availability. In this review, we communicate the current state of the art regarding these archaeal networks and their dynamic properties. We compare the topology of these archaeal networks to those known for bacteria to highlight conserved and unique aspects. We present a new computational model for an exemplar archaeal network, aiming to lay the groundwork toward understanding general principles that unify the dynamic function of integrated metabolic-transcription networks across archaea and bacteria.
| S-EPMC7289023 | biostudies-literature
Project description:To develop and investigate detailed mathematical models of metabolic processes is one of the primary challenges in systems biology. However, despite considerable advance in the topological analysis of metabolic networks, kinetic modeling is still often severely hampered by inadequate knowledge of the enzyme-kinetic rate laws and their associated parameter values. Here we propose a method that aims to give a quantitative account of the dynamical capabilities of a metabolic system, without requiring any explicit information about the functional form of the rate equations. Our approach is based on constructing a local linear model at each point in parameter space, such that each element of the model is either directly experimentally accessible or amenable to a straightforward biochemical interpretation. This ensemble of local linear models, encompassing all possible explicit kinetic models, then allows for a statistical exploration of the comprehensive parameter space. The method is exemplified on two paradigmatic metabolic systems: the glycolytic pathway of yeast and a realistic-scale representation of the photosynthetic Calvin cycle.
| S-EPMC1524928 | biostudies-literature
Project description:Stochasticity is a key characteristic of intracellular processes such as gene regulation and chemical signalling. Therefore, characterizing stochastic effects in biochemical systems is essential to understand the complex dynamics of living things. Mathematical idealizations of biochemically reacting systems must be able to capture stochastic phenomena. While robust theory exists to describe such stochastic models, the computational challenges in exploring these models can be a significant burden in practice since realistic models are analytically intractable. Determining the expected behaviour and variability of a stochastic biochemical reaction network requires many probabilistic simulations of its evolution. Using a biochemical reaction network model to assist in the interpretation of time-course data from a biological experiment is an even greater challenge due to the intractability of the likelihood function for determining observation probabilities. These computational challenges have been subjects of active research for over four decades. In this review, we present an accessible discussion of the major historical developments and state-of-the-art computational techniques relevant to simulation and inference problems for stochastic biochemical reaction network models. Detailed algorithms for particularly important methods are described and complemented with Matlab® implementations. As a result, this review provides a practical and accessible introduction to computational methods for stochastic models within the life sciences community.
| S-EPMC6408336 | biostudies-literature