Dataset Information

Are Quasi-Steady-State Approximated Models Suitable for Quantifying Intrinsic Noise Accurately?

ABSTRACT: Large gene regulatory networks (GRN) are often modeled with quasi-steady-state approximation (QSSA) to reduce the huge computational time required for intrinsic noise quantification using Gillespie stochastic simulation algorithm (SSA). However, the question still remains whether the stochastic QSSA model measures the intrinsic noise as accurately as the SSA performed for a detailed mechanistic model or not? To address this issue, we have constructed mechanistic and QSSA models for few frequently observed GRNs exhibiting switching behavior and performed stochastic simulations with them. Our results strongly suggest that the performance of a stochastic QSSA model in comparison to SSA performed for a mechanistic model critically relies on the absolute values of the mRNA and protein half-lives involved in the corresponding GRN. The extent of accuracy level achieved by the stochastic QSSA model calculations will depend on the level of bursting frequency generated due to the absolute value of the half-life of either mRNA or protein or for both the species. For the GRNs considered, the stochastic QSSA quantifies the intrinsic noise at the protein level with greater accuracy and for larger combinations of half-life values of mRNA and protein, whereas in case of mRNA the satisfactory accuracy level can only be reached for limited combinations of absolute values of half-lives. Further, we have clearly demonstrated that the abundance levels of mRNA and protein hardly matter for such comparison between QSSA and mechanistic models. Based on our findings, we conclude that QSSA model can be a good choice for evaluating intrinsic noise for other GRNs as well, provided we make a rational choice based on experimental half-life values available in literature.

SUBMITTER: Sengupta D

PROVIDER: S-EPMC4556639 | biostudies-literature | 2015

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Are Quasi-Steady-State Approximated Models Suitable for Quantifying Intrinsic Noise Accurately?

Sengupta Dola D Kar Sandip S

PloS one 20150901 9

Large gene regulatory networks (GRN) are often modeled with quasi-steady-state approximation (QSSA) to reduce the huge computational time required for intrinsic noise quantification using Gillespie stochastic simulation algorithm (SSA). However, the question still remains whether the stochastic QSSA model measures the intrinsic noise as accurately as the SSA performed for a detailed mechanistic model or not? To address this issue, we have constructed mechanistic and QSSA models for few frequentl ...[more]

PMID: 26327626

Similar Datasets

Project description:During tissue development, patterns of gene expression determine the spatial arrangement of cell types. In many cases, gradients of secreted signalling molecules-morphogens-guide this process by controlling downstream transcriptional networks. A mechanism commonly used in these networks to convert the continuous information provided by the gradient into discrete transitions between adjacent cell types is the genetic toggle switch, composed of cross-repressing transcriptional determinants. Previous analyses have emphasised the steady state output of these mechanisms. Here, we explore the dynamics of the toggle switch and use exact numerical simulations of the kinetic reactions, the corresponding Chemical Langevin Equation, and Minimum Action Path theory to establish a framework for studying the effect of gene expression noise on patterning time and boundary position. This provides insight into the time scale, gene expression trajectories and directionality of stochastic switching events between cell states. Taking gene expression noise into account predicts that the final boundary position of a morphogen-induced toggle switch, although robust to changes in the details of the noise, is distinct from that of the deterministic system. Moreover, the dramatic increase in patterning time close to the boundary predicted from the deterministic case is substantially reduced. The resulting stochastic switching introduces differences in patterning time along the morphogen gradient that result in a patterning wave propagating away from the morphogen source with a velocity determined by the intrinsic noise. The wave sharpens and slows as it advances and may never reach steady state in a biologically relevant time. This could explain experimentally observed dynamics of pattern formation. Together the analysis reveals the importance of dynamical transients for understanding morphogen-driven transcriptional networks and indicates that gene expression noise can qualitatively alter developmental patterning.

Project description:Oxygen levels in cancerous tissue can have a significant effect on treatment response: hypoxic tissue is both more radioresistant and more chemoresistant than well-oxygenated tissue. While recent advances in medical imaging have facilitated real-time observation of macroscopic oxygenation, the underlying physics limits the resolution to the millimetre domain, whereas oxygen tension varies over a micrometre scale. If the distribution of oxygen in the tumour micro-environment can be accurately estimated, then the effect of potential dose escalation to these hypoxic regions could be better modelled, allowing more realistic simulation of biologically adaptive treatments. Reaction-diffusion models are commonly used for modelling oxygen dynamics, with a variety of functional forms assumed for the dependence of oxygen consumption rate (OCR) on cellular status and local oxygen availability. In this work, we examine reaction-diffusion models of oxygen consumption in spherically and cylindrically symmetric geometries. We consider two different descriptions of oxygen consumption: one in which the rate of consumption is constant and one in which it varies with oxygen tension in a hyperbolic manner. In each case, we derive analytic approximations to the steady-state oxygen distribution, which are shown to closely match the numerical solutions of the equations and accurately predict the extent to which oxygen can diffuse. The derived expressions relate the limit to which oxygen can diffuse into a tissue to the OCR of that tissue. We also demonstrate that differences between these functional forms are likely to be negligible within the range of literature estimates of the hyperbolic oxygen constant, suggesting that the constant consumption rate approximation suffices for modelling oxygen dynamics for most values of OCR. These approximations also allow the rapid identification of situations where hyperbolic consumption forms can result in significant differences from constant consumption rate models, and so can reduce the computational workload associated with numerical solutions, by estimating both the oxygen diffusion distances and resultant oxygen profile. Such analysis may be useful for parameter fitting in large imaging datasets and histological sections, and allows easy quantification of projected differences between functional forms of OCR.

Dataset Information

Are Quasi-Steady-State Approximated Models Suitable for Quantifying Intrinsic Noise Accurately?

Publications

Are Quasi-Steady-State Approximated Models Suitable for Quantifying Intrinsic Noise Accurately?

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets