Robustness analysis of the detailed kinetic model of an ErbB signaling network by using dynamic sensitivity.
Ontology highlight
ABSTRACT:
SUBMITTER: Masunaga H
PROVIDER: S-EPMC5443533 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Intracellular signaling/synthetic pathways are being increasingly extensively characterized. However, while these pathways can be displayed in static diagrams, in reality they exist with a degree of dynamic complexity that is responsible for heterogeneous cellular behavior. Multiple parallel pathways exist and interact concurrently, limiting the ability to integrate the various identified mechanisms into a cohesive whole. Computational methods have been suggested as a means of concatenating this knowledge to aid in the understanding of overall system dynamics. Since the eventual goal of biomedical research is the identification and development of therapeutic modalities, computational representation must have sufficient detail to facilitate this 'engineering' process. Adding to the challenge, this type of representation must occur in a perpetual state of incomplete knowledge. We present a modeling approach to address this challenge that is both detailed and qualitative. This approach is termed 'dynamic knowledge representation,' and is intended to be an integrated component of the iterative cycle of scientific discovery.BioNetGen (BNG), a software platform for modeling intracellular signaling pathways, was used to model the toll-like receptor 4 (TLR-4) signal transduction cascade. The informational basis of the model was a series of reference papers on modulation of (TLR-4) signaling, and some specific primary research papers to aid in the characterization of specific mechanistic steps in the pathway. This model was detailed with respect to the components of the pathway represented, but qualitative with respect to the specific reaction coefficients utilized to execute the reactions. Responsiveness to simulated lipopolysaccharide (LPS) administration was measured by tumor necrosis factor (TNF) production. Simulation runs included evaluation of initial dose-dependent response to LPS administration at 10, 100, 1000 and 10,000, and a subsequent examination of preconditioning behavior with increasing LPS at 10, 100, 1000 and 10,000 and a secondary dose of LPS at 10,000 administered at approximately 27h of simulated time. Simulations of 'knockout' versions of the model allowed further examination of the interactions within the signaling cascade.The model demonstrated a dose-dependent TNF response curve to increasing stimulus by LPS. Preconditioning simulations demonstrated a similar dose-dependency of preconditioning doses leading to attenuation of response to subsequent LPS challenge - a 'tolerance' dynamic. These responses match dynamics reported in the literature. Furthermore, the simulated 'knockout' results suggested the existence and need for dual negative feedback control mechanisms, represented by the zinc ring-finger protein A20 and inhibitor kappa B proteins (IkappaB), in order for both effective attenuation of the initial stimulus signal and subsequent preconditioned 'tolerant' behavior.We present an example of detailed, qualitative dynamic knowledge representation using the TLR-4 signaling pathway, its control mechanisms and overall behavior with respect to preconditioning. The intent of this approach is to demonstrate a method of translating the extensive mechanistic knowledge being generated at the basic science level into an executable framework that can provide a means of 'conceptual model verification.' This allows for both the 'checking' of the dynamic consequences of a mechanistic hypothesis and the creation of a modular component of an overall model directed at the engineering goal of biomedical research. It is hoped that this paper will increase the use of knowledge representation and communication in this fashion, and facilitate the concatenation and integration of community-wide knowledge.
| S-EPMC2652645 | biostudies-literature
Project description:Ligand-induced homo- and hetero-dimer formation of ErbB receptors results in different biological outcomes irrespective of recruitment and activation of similar effector proteins. Earlier experimental research indicated that cells expressing both EGFR (epidermal growth factor receptor) and the ErbB4 receptor (E1/4 cells) induced E1/4 cell-specific B-Raf activation and higher extracellular signal-regulated kinase (ERK) activation, followed by cellular transformation, than cells solely expressing EGFR (E1 cells) in Chinese hamster ovary (CHO) cells. Since our experimental data revealed the presence of positive feedback by ERK on upstream pathways, it was estimated that the cross-talk/feedback pathway structure of the Raf-MEK-ERK cascade might affect ERK activation dynamics in our cell system. To uncover the regulatory mechanism concerning the ERK dynamics, we used topological models and performed parameter estimation for all candidate structures that possessed ERK-mediated positive feedback regulation of Raf. The structure that reliably reproduced a series of experimental data regarding signal amplitude and duration of the signaling molecules was selected as a solution. We found that the pathway structure is characterized by ERK-mediated positive feedback regulation of B-Raf and B-Raf-mediated negative regulation of Raf-1. Steady-state analysis of the estimated structure indicated that the amplitude of Ras activity might critically affect ERK activity through ERK-B-Raf positive feedback coordination with sustained B-Raf activation in E1/4 cells. However, Rap1 that positively regulates B-Raf activity might be less effective concerning ERK and B-Raf activity. Furthermore, we investigated how such Ras activity in E1/4 cells can be regulated by EGFR/ErbB4 heterodimer-mediated signaling. From a sensitivity analysis of the detailed upstream model for Ras activation, we concluded that Ras activation dynamics is dominated by heterodimer-mediated signaling coordination with a large initial speed of dimerization when the concentration of the ErbB4 receptor is considerably high. Such characteristics of the signaling cause the preferential binding of the Grb2-SOS complex to heterodimer-mediated signaling molecules.
| S-EPMC2262155 | biostudies-literature
Project description:BackgroundIntegrative analysis between dynamical modeling of metabolic networks and data obtained from high throughput technology represents a worthy effort toward a holistic understanding of the link among phenotype and dynamical response. Even though the theoretical foundation for modeling metabolic network has been extensively treated elsewhere, the lack of kinetic information has limited the analysis in most of the cases. To overcome this constraint, we present and illustrate a new statistical approach that has two purposes: integrate high throughput data and survey the general dynamical mechanisms emerging for a slightly perturbed metabolic network.Methodology/principal findingsThis paper presents a statistic framework capable to study how and how fast the metabolites participating in a perturbed metabolic network reach a steady-state. Instead of requiring accurate kinetic information, this approach uses high throughput metabolome technology to define a feasible kinetic library, which constitutes the base for identifying, statistical and dynamical properties during the relaxation. For the sake of illustration we have applied this approach to the human Red blood cell metabolism (hRBC) and its capacity to predict temporal phenomena was evaluated. Remarkable, the main dynamical properties obtained from a detailed kinetic model in hRBC were recovered by our statistical approach. Furthermore, robust properties in time scale and metabolite organization were identify and one concluded that they are a consequence of the combined performance of redundancies and variability in metabolite participation.Conclusions/significanceIn this work we present an approach that integrates high throughput metabolome data to define the dynamic behavior of a slightly perturbed metabolic network where kinetic information is lacking. Having information of metabolite concentrations at steady-state, this method has significant relevance due its potential scope to analyze others genome scale metabolic reconstructions. Thus, I expect this approach will significantly contribute to explore the relationship between dynamic and physiology in other metabolic reconstructions, particularly those whose kinetic information is practically nulls. For instances, I envisage that this approach can be useful in genomic medicine or pharmacogenomics, where the estimation of time scales and the identification of metabolite organization may be crucial to characterize and identify (dis)functional stages.
| S-EPMC2654918 | biostudies-literature
Project description:Integrin signaling network is responsible for regulating a wide variety of fundamental biological processes ranging from cell survival to cell death. While individual components of the network have been studied through experimental and computational methods, the network robustness and the flow of information through the network have not been characterized in a quantitative framework. Using a probability based model implemented through GRID computing, we approach the reduced signaling network and show that the network is highly robust and the final stable steady state is independent of the initial configurations. However, the path from the initial and the final state is intrinsically dependent on the state of the input nodes. Our results demonstrate a rugged funnel-like landscape for the signaling network where the final state is unique, but the paths are dependent on initial conditions.
| S-EPMC3982507 | biostudies-other
Project description:The ErbB signaling pathways, which regulate diverse physiological responses such as cell survival, proliferation and motility, have been subjected to extensive molecular analysis. Nonetheless, it remains poorly understood how different ligands induce different responses and how this is affected by oncogenic mutations. To quantify signal flow through ErbB-activated pathways we have constructed, trained and analyzed a mass action model of immediate-early signaling involving ErbB1-4 receptors (EGFR, HER2/Neu2, ErbB3 and ErbB4), and the MAPK and PI3K/Akt cascades. We find that parameter sensitivity is strongly dependent on the feature (e.g. ERK or Akt activation) or condition (e.g. EGF or heregulin stimulation) under examination and that this context dependence is informative with respect to mechanisms of signal propagation. Modeling predicts log-linear amplification so that significant ERK and Akt activation is observed at ligand concentrations far below the K(d) for receptor binding. However, MAPK and Akt modules isolated from the ErbB model continue to exhibit switch-like responses. Thus, key system-wide features of ErbB signaling arise from nonlinear interaction among signaling elements, the properties of which appear quite different in context and in isolation.
| S-EPMC2644173 | biostudies-literature
Project description:Molecular dynamics simulations contain detailed kinetic information related to the functional states of proteins and macromolecules, but this information is obscured by the high dimensionality of configurational space. Markov state models and transition network models are widely applied to extract kinetic descriptors from equilibrium molecular dynamics simulations. In this study, we developed the Directed Kinetic Transition Network (DKTN)-a graph representation of a master equation which is appropriate for describing nonequilibrium kinetics. DKTN models the transition rate matrix among different states under detailed balance. Adopting the mixing time from the Markov chain, we use the half mixing time as the criterion to identify critical state transition regarding the protein conformational change. The similarity between the master equation and the Kolmogorov equation suggests that the DKTN model can be reformulated into the continuous-time Markov chain model, which is a general case of the Markov chain without a specific lag time. We selected a photo-sensitive protein, vivid, as a model system to illustrate the usage of the DKTN model. Overall, the DKTN model provides a graph representation of the master equation based on chemical kinetics to model the protein conformational change without the underlying assumption of the Markovian property.
| S-EPMC6800283 | biostudies-literature
Project description:Biological systems are known to be both robust and evolvable to internal and external perturbations, but what causes these apparently contradictory properties? We used Boolean network modeling and attractor landscape analysis to investigate the evolvability and robustness of the human signaling network. Our results show that the human signaling network can be divided into an evolvable core where perturbations change the attractor landscape in state space, and a robust neighbor where perturbations have no effect on the attractor landscape. Using chemical inhibition and overexpression of nodes, we validated that perturbations affect the evolvable core more strongly than the robust neighbor. We also found that the evolvable core has a distinct network structure, which is enriched in feedback loops, and features a higher degree of scale-freeness and longer path lengths connecting the nodes. In addition, the genes with high evolvability scores are associated with evolvability-related properties such as rapid evolvability, low species broadness, and immunity whereas the genes with high robustness scores are associated with robustness-related properties such as slow evolvability, high species broadness, and oncogenes. Intriguingly, US Food and Drug Administration-approved drug targets have high evolvability scores whereas experimental drug targets have high robustness scores.
| S-EPMC4117429 | biostudies-literature
Project description:Cell state transition is at the core of biological processes in metazoan, which includes cell differentiation, epithelial-to-mesenchymal transition (EMT) and cell reprogramming. In these cases, it is important to understand the molecular mechanism of cellular stability and how the transitions happen between different cell states, which is controlled by a gene regulatory network (GRN) hard-wired in the genome. Here we use Boolean modeling of GRN to study the cell state transition of EMT and systematically compare four available methods to calculate the cellular stability of three cell states in EMT in both normal and genetically mutated cases. The results produced from four methods generally agree but do not totally agree with each other. We show that distribution of one-degree neighborhood of cell states, which are the nearest states by Hamming distance, causes the difference among the methods. From that, we propose a new method based on one-degree neighborhood, which is the simplest one and agrees with other methods to estimate the cellular stability in all scenarios of our EMT model. This new method will help the researchers in the field of cell differentiation and cell reprogramming to calculate cellular stability using Boolean model, and then rationally design their experimental protocols to manipulate the cell state transition.
| S-EPMC6089891 | biostudies-other
Project description:Deregulation of ErbB signaling plays a key role in the progression of multiple human cancers. To help understand ErbB signaling quantitatively, in this work we combine traditional experiments with computational modeling, building a model that describes how stimulation of all four ErbB receptors with epidermal growth factor (EGF) and heregulin (HRG) leads to activation of two critical downstream proteins, extracellular-signal-regulated kinase (ERK) and Akt. Model analysis and experimental validation show that (i) ErbB2 overexpression, which occurs in approximately 25% of all breast cancers, transforms transient EGF-induced signaling into sustained signaling, (ii) HRG-induced ERK activity is much more robust to the ERK cascade inhibitor U0126 than EGF-induced ERK activity, and (iii) phosphoinositol-3 kinase is a major regulator of post-peak but not pre-peak EGF-induced ERK activity. Sensitivity analysis leads to the hypothesis that ERK activation is robust to parameter perturbation at high ligand doses, while Akt activation is not.
| S-EPMC2132449 | biostudies-literature
Project description:Robustness and sensitivity of responses generated by cell signaling networks has been associated with survival and evolvability of organisms. However, existing methods analyzing robustness and sensitivity of signaling networks ignore the experimentally observed cell-to-cell variations of protein abundances and cell functions or contain ad hoc assumptions. We propose and apply a data-driven maximum entropy based method to quantify robustness and sensitivity of Escherichia coli (E. coli) chemotaxis signaling network. Our analysis correctly rank orders different models of E. coli chemotaxis based on their robustness and suggests that parameters regulating cell signaling are evolutionary selected to vary in individual cells according to their abilities to perturb cell functions. Furthermore, predictions from our approach regarding distribution of protein abundances and properties of chemotactic responses in individual cells based on cell population averaged data are in excellent agreement with their experimental counterparts. Our approach is general and can be used to evaluate robustness as well as generate predictions of single cell properties based on population averaged experimental data in a wide range of cell signaling systems.
| S-EPMC3955891 | biostudies-literature