Dataset Information

Mathematical identification of critical reactions in the interlocked feedback model.

ABSTRACT: Dynamic simulations are necessary for understanding the mechanism of how biochemical networks generate robust properties to environmental stresses or genetic changes. Sensitivity analysis allows the linking of robustness to network structure. However, it yields only local properties regarding a particular choice of plausible parameter values, because it is hard to know the exact parameter values in vivo. Global and firm results are needed that do not depend on particular parameter values. We propose mathematical analysis for robustness (MAR) that consists of the novel evolutionary search that explores all possible solution vectors of kinetic parameters satisfying the target dynamics and robustness analysis. New criteria, parameter spectrum width and the variability of solution vectors for parameters, are introduced to determine whether the search is exhaustive. In robustness analysis, in addition to single parameter sensitivity analysis, robustness to multiple parameter perturbation is defined. Combining the sensitivity analysis and the robustness analysis to multiple parameter perturbation enables identifying critical reactions. Use of MAR clearly identified the critical reactions responsible for determining the circadian cycle in the Drosophila interlocked circadian clock model. In highly robust models, while the parameter vectors are greatly varied, the critical reactions with a high sensitivity are uniquely determined. Interestingly, not only the per-tim loop but also the dclk-cyc loop strongly affect the period of PER, although the dclk-cyc loop hardly changes its amplitude and it is not potentially influential. In conclusion, MAR is a powerful method to explore wide parameter space without human-biases and to link a robust property to network architectures without knowing the exact parameter values. MAR identifies the reactions critically responsible for determining the period and amplitude in the interlocked feedback model and suggests that the circadian clock intensively evolves or designs the kinetic parameters so that it creates a highly robust cycle.

SUBMITTER: Kurata H

PROVIDER: S-EPMC2040204 | biostudies-literature | 2007 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Mathematical identification of critical reactions in the interlocked feedback model.

Kurata Hiroyuki H Tanaka Takayuki T Ohnishi Fumitaka F

PloS one 20071031 10

Dynamic simulations are necessary for understanding the mechanism of how biochemical networks generate robust properties to environmental stresses or genetic changes. Sensitivity analysis allows the linking of robustness to network structure. However, it yields only local properties regarding a particular choice of plausible parameter values, because it is hard to know the exact parameter values in vivo. Global and firm results are needed that do not depend on particular parameter values. We pro ...[more]

PMID: 17971866

Similar Datasets

Project description:Maintaining normal potassium (K+) concentrations in the extra- and intracellular fluid is critical for cell function. K+ homeostasis is achieved by ensuring proper distribution between extra- and intracellular fluid compartments and by matching K+ excretion with intake. The Na+-K+-ATPase pump facilitates K+ uptake into the skeletal muscle, where most K+ is stored. Na+-K+-ATPase activity is stimulated by insulin and aldosterone. The kidneys regulate long term K+ homeostasis by controlling the amount of K+ excreted through urine. Renal handling of K+ is mediated by a number of regulatory mechanisms, including an aldosterone-mediated feedback control, in which high extracellular K+ concentration stimulates aldosterone secretion, which enhances urine K+ excretion, and a gastrointestinal feedforward control mechanism, in which dietary K+ intake increases K+ excretion. Recently, a muscle-kidney cross talk signal has been hypothesized, where the K+ concentration in skeletal muscle cells directly affects urine K+ excretion without changes in extracellular K+ concentration. To understand how these mechanisms coordinate under different K+ challenges, we have developed a compartmental model of whole-body K+ regulation. The model represents the intra- and extracellular fluid compartments in a human (male) as well as a detailed kidney compartment. We included (i) the gastrointestinal feedforward control mechanism, (ii) the effect of insulin and (iii) aldosterone on Na+-K+-ATPase K+ uptake, and (iv) aldosterone stimulation of renal K+ secretion. We used this model to investigate the impact of regulatory mechanisms on K+ homeostasis. Model predictions showed how the regulatory mechanisms synthesize to ensure that the extra- and intracelluller fluid K+ concentrations remain in normal range in times of K+ loading and fasting. Additionally, we predict that without the hypothesized muscle-kidney cross talk signal, the model was unable to predict a return to normal extracellular K+ concentration after a period of high K+ loading or depletion.

Project description:BackgroundTransforming growth factor β (TGF-β) signalling regulates the development of embryos and tissue homeostasis in adults. In conjunction with other oncogenic changes, long-term perturbation of TGF-β signalling is associated with cancer metastasis. Although TGF-β signalling can be complex, many of the signalling components are well defined, so it is possible to develop mathematical models of TGF-β signalling using reduction and scaling methods. The parameterization of our TGF-β signalling model is consistent with experimental data.ResultsWe developed our mathematical model for the TGF-β signalling pathway, i.e. the RF- model of TGF-β signalling, using the "rapid equilibrium assumption" to reduce the network of TGF-β signalling reactions based on the time scales of the individual reactions. By adding time-delayed positive feedback to the inherent time-delayed negative feedback for TGF-β signalling. We were able to simulate the sigmoidal, switch-like behaviour observed for the concentration dependence of long-term (> 3 hours) TGF-β stimulation. Computer simulations revealed the vital role of the coupling of the positive and negative feedback loops on the regulation of the TGF-β signalling system. The incorporation of time-delays for the negative feedback loop improved the accuracy, stability and robustness of the model. This model reproduces both the short-term and long-term switching responses for the intracellular signalling pathways at different TGF-β concentrations. We have tested the model against experimental data from MEF (mouse embryonic fibroblasts) WT, SV40-immortalized MEFs and Gp130 F/F MEFs. The predictions from the RF- model are consistent with the experimental data.ConclusionsSignalling feedback loops are required to model TGF-β signal transduction and its effects on normal and cancer cells. We focus on the effects of time-delayed feedback loops and their coupling to ligand stimulation in this system. The model was simplified and reduced to its key components using standard methods and the rapid equilibrium assumption. We detected differences in short-term and long-term signal switching. The results from the RF- model compare well with experimental data and predict the dynamics of TGF-β signalling in cancer cells with different mutations.

Project description:BackgroundFeedback loops in gene regulatory networks play pivotal roles in governing functional dynamics of cells. Systems approaches demonstrated characteristic dynamical features, including multistability and oscillation, of positive and negative feedback loops. Recent experiments and theories have implicated highly interconnected feedback loops (high-feedback loops) in additional nonintuitive functions, such as controlling cell differentiation rate and multistep cell lineage progression. However, it remains challenging to identify and visualize high-feedback loops in complex gene regulatory networks due to the myriad of ways in which the loops can be combined. Furthermore, it is unclear whether the high-feedback loop structures with these potential functions are widespread in biological systems. Finally, it remains challenging to understand diverse dynamical features, such as high-order multistability and oscillation, generated by individual networks containing high-feedback loops. To address these problems, we developed HiLoop, a toolkit that enables discovery, visualization, and analysis of several types of high-feedback loops in large biological networks.ResultsHiLoop not only extracts high-feedback structures and visualize them in intuitive ways, but also quantifies the enrichment of overrepresented structures. Through random parameterization of mathematical models derived from target networks, HiLoop presents characteristic features of the underlying systems, including complex multistability and oscillations, in a unifying framework. Using HiLoop, we were able to analyze realistic gene regulatory networks containing dozens to hundreds of genes, and to identify many small high-feedback systems. We found more than a 100 human transcription factors involved in high-feedback loops that were not studied previously. In addition, HiLoop enabled the discovery of an enrichment of high feedback in pathways related to epithelial-mesenchymal transition.ConclusionsHiLoop makes the study of complex networks accessible without significant computational demands. It can serve as a hypothesis generator through identification and modeling of high-feedback subnetworks, or as a quantification method for motif enrichment analysis. As an example of discovery, we found that multistep cell lineage progression may be driven by either specific instances of high-feedback loops with sparse appearances, or generally enriched topologies in gene regulatory networks. We expect HiLoop's usefulness to increase as experimental data of regulatory networks accumulate. Code is freely available for use or extension at https://github.com/BenNordick/HiLoop .

Dataset Information

Mathematical identification of critical reactions in the interlocked feedback model.

Publications

Mathematical identification of critical reactions in the interlocked feedback model.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets