Modeling delay in genetic networks: from delay birth-death processes to delay stochastic differential equations.
Ontology highlight
ABSTRACT:
SUBMITTER: Gupta C
PROVIDER: S-EPMC4039735 | biostudies-other | 2014 May
REPOSITORIES: biostudies-other
Similar Datasets
Project description:Biological regulatory systems, such as cell signaling networks, nervous systems and ecological webs, consist of complex dynamical interactions among many components. Network motif models focus on small sub-networks to provide quantitative insight into overall behavior. However, such models often overlook time delays either inherent to biological processes or associated with multi-step interactions. Here we systematically examine explicit-delay versions of the most common network motifs via delay differential equation (DDE) models, both analytically and numerically. We find many broadly applicable results, including parameter reduction versus canonical ordinary differential equation (ODE) models, analytical relations for converting between ODE and DDE models, criteria for when delays may be ignored, a complete phase space for autoregulation, universal behaviors of feedforward loops, a unified Hill-function logic framework, and conditions for oscillations and chaos. We conclude that explicit-delay modeling simplifies the phenomenology of many biological networks and may aid in discovering new functional motifs.
| S-EPMC7979834 | biostudies-literature
Project description:Contaminated sites are recognized as the "hotspot" for the development and spread of antibiotic resistance in environmental bacteria. It is very challenging to understand mechanism of development of antibiotic resistance in polluted environment in the presence of different anthropogenic pollutants. Uncertainties in the environmental processes adds complexity to the development of resistance. This study attempts to develop mathematical model by using stochastic partial differential equations for the transport of fluoroquinolone and its resistant bacteria in riverine environment. Poisson's process is assumed for the diffusion approximation in the stochastic partial differential equations (SPDE). Sensitive analysis is performed to evaluate the parameters and variables for their influence over the model outcome. Based on their sensitivity, the model parameters and variables are chosen and classified into environmental, demographic, and anthropogenic categories to investigate the sources of stochasticity. Stochastic partial differential equations are formulated for the state variables in the model. This SPDE model is then applied to the 100 km stretch of river Musi (South India) and simulations are carried out to assess the impact of stochasticity in model variables on the resistant bacteria population in sediments. By employing the stochasticity in model variables and parameters we came to know that environmental and anthropogenic variations are not able to affect the resistance dynamics at all. Demographic variations are able to affect the distribution of resistant bacteria population uniformly with standard deviation between 0.087 and 0.084, however, is not significant to have any biological relevance to it. The outcome of the present study is helpful in simplifying the model for practical applications. This study is an ongoing effort to improve the model for the transport of antibiotics and transport of antibiotic resistant bacteria in polluted river. There is a wide gap between the knowledge of stochastic resistant bacterial growth dynamics and the knowledge of transport of antibiotic resistance in polluted aquatic environment, this study is one step towards filling up that gap.
| S-EPMC7494867 | biostudies-literature
Project description:Time series analysis with nonlinear delay differential equations (DDEs) reveals nonlinear as well as spectral properties of the underlying dynamical system. Here, global DDE models were used to analyze 5 min data segments of electrocardiographic (ECG) recordings in order to capture distinguishing features for different heart conditions such as normal heart beat, congestive heart failure, and atrial fibrillation. The number of terms and delays in the model as well as the order of nonlinearity of the model have to be selected that are the most discriminative. The DDE model form that best separates the three classes of data was chosen by exhaustive search up to third order polynomials. Such an approach can provide deep insight into the nature of the data since linear terms of a DDE correspond to the main time-scales in the signal and the nonlinear terms in the DDE are related to nonlinear couplings between the harmonic signal parts. The DDEs were able to detect atrial fibrillation with an accuracy of 72%, congestive heart failure with an accuracy of 88%, and normal heart beat with an accuracy of 97% from 5 min of ECG, a much shorter time interval than required to achieve comparable performance with other methods.
| S-EPMC3710263 | biostudies-other
Project description:Network models extend evolutionary game theory to settings with spatial or social structure and have provided key insights on the mechanisms underlying the evolution of cooperation. However, network models have also proven sensitive to seemingly small details of the model architecture. Here we investigate two popular biologically motivated models of evolution in finite populations: Death-Birth (DB) and Birth-Death (BD) processes. In both cases reproduction is proportional to fitness and death is random; the only difference is the order of the two events at each time step. Although superficially similar, under DB cooperation may be favoured in structured populations, while under BD it never is. This is especially troubling as natural populations do not follow a strict one birth then one death regimen (or vice versa); such constraints are introduced to make models more tractable. Whether structure can promote the evolution of cooperation should not hinge on a simplifying assumption. Here, we propose a mixed rule where in each time step DB is used with probability ? and BD is used with probability 1-?. We derive the conditions for selection favouring cooperation under the mixed rule for all social dilemmas. We find that the only qualitatively different outcome occurs when using just BD (? = 0). This case admits a natural interpretation in terms of kin competition counterbalancing the effect of kin selection. Finally we show that, for any mixed BD-DB update and under weak selection, cooperation is never inhibited by population structure for any social dilemma, including the Snowdrift Game.
| S-EPMC3557300 | biostudies-literature
Project description:The use of RNA-seq as the preferred method for the discovery and validation of small RNA biomarkers is hindered by high variability and biased sequence counts. In this paper we develop a statistical model for sequence counts that accounts for ligase bias and stochastic variation in library amplification steps and sequencing depth variation. Our analytical contributions are the description of the Linear Quadratic (LQ) relation between the mean and variance of the sequence counts in an RNA-seq experiment and the derivation of the Poisson truncated mixture as the underlying probability distribution for RNA-seq data. Using a large number of sequencing datasets, we demonstrate here how one can use this modeling framework to calculate empirical correction factors for ligase bias, while accounting for random variation in sequence counts. Bias correction may remove the majority of bias in the absence of differential expression and more than 40% of the bias in the presence of variable expression of miRNAs. Empirical bias correction factors appear to be nearly constant over at least one and up to four orders of magnitude of total RNA input and independent of sample composition.
2017-03-15 | GSE93399 | GEO
Project description:Many models in Systems Biology are described as a system of Ordinary Differential Equations, which allows for transient, steady-state or bifurcation analysis when kinetic information is available. Complementary structure-related qualitative analysis techniques have become increasingly popular in recent years, like qualitative model checking or pathway analysis (elementary modes, invariants, flux balance analysis, graph-based analyses, chemical organization theory, etc.). They do not rely on kinetic information but require a well-defined structure as stochastic analysis techniques equally do. In this article, we look into the structure inference problem for a model described by a system of Ordinary Differential Equations and provide conditions for the uniqueness of its solution. We describe a method to extract a structured reaction network model, represented as a bipartite multigraph, for example, a continuous Petri net (CPN), from a system of Ordinary Differential Equations (ODEs). A CPN uniquely defines an ODE, and each ODE can be transformed into a CPN. However, it is not obvious under which conditions the transformation of an ODE into a CPN is unique, that is, when a given ODE defines exactly one CPN. We provide biochemically relevant sufficient conditions under which the derived structure is unique and counterexamples showing the necessity of each condition. Our method is implemented and available; we illustrate it on some signal transduction models from the BioModels database. A prototype implementation of the method is made available to modellers at http://contraintes.inria.fr/~soliman/ode2pn.html, and the data mentioned in the "Results" section at http://contraintes.inria.fr/~soliman/ode2pn_data/. Our results yield a new recommendation for the import/export feature of tools supporting the SBML exchange format.
| S-EPMC3008708 | biostudies-literature
Project description:BACKGROUND: Various computational models have been of interest due to their use in the modelling of gene regulatory networks (GRNs). As a logical model, probabilistic Boolean networks (PBNs) consider molecular and genetic noise, so the study of PBNs provides significant insights into the understanding of the dynamics of GRNs. This will ultimately lead to advances in developing therapeutic methods that intervene in the process of disease development and progression. The applications of PBNs, however, are hindered by the complexities involved in the computation of the state transition matrix and the steady-state distribution of a PBN. For a PBN with n genes and N Boolean networks, the complexity to compute the state transition matrix is O(nN22n) or O(nN2n) for a sparse matrix. RESULTS: This paper presents a novel implementation of PBNs based on the notions of stochastic logic and stochastic computation. This stochastic implementation of a PBN is referred to as a stochastic Boolean network (SBN). An SBN provides an accurate and efficient simulation of a PBN without and with random gene perturbation. The state transition matrix is computed in an SBN with a complexity of O(nL2n), where L is a factor related to the stochastic sequence length. Since the minimum sequence length required for obtaining an evaluation accuracy approximately increases in a polynomial order with the number of genes, n, and the number of Boolean networks, N, usually increases exponentially with n, L is typically smaller than N, especially in a network with a large number of genes. Hence, the computational efficiency of an SBN is primarily limited by the number of genes, but not directly by the total possible number of Boolean networks. Furthermore, a time-frame expanded SBN enables an efficient analysis of the steady-state distribution of a PBN. These findings are supported by the simulation results of a simplified p53 network, several randomly generated networks and a network inferred from a T cell immune response dataset. An SBN can also implement the function of an asynchronous PBN and is potentially useful in a hybrid approach in combination with a continuous or single-molecule level stochastic model. CONCLUSIONS: Stochastic Boolean networks (SBNs) are proposed as an efficient approach to modelling gene regulatory networks (GRNs). The SBN approach is able to recover biologically-proven regulatory behaviours, such as the oscillatory dynamics of the p53-Mdm2 network and the dynamic attractors in a T cell immune response network. The proposed approach can further predict the network dynamics when the genes are under perturbation, thus providing biologically meaningful insights for a better understanding of the dynamics of GRNs. The algorithms and methods described in this paper have been implemented in Matlab packages, which are attached as Additional files.
| S-EPMC3532238 | biostudies-literature
Project description:BackgroundThis work explores the quantitative characteristics of the local transcriptional regulatory network based on the availability of time dependent gene expression data sets. The dynamics of the gene expression level are fitted via a stochastic differential equation model, yielding a set of specific regulators and their contribution.ResultsWe show that a beta sigmoid function that keeps track of temporal parameters is a novel prototype of a regulatory function, with the effect of improving the performance of the profile prediction. The stochastic differential equation model follows well the dynamic of the gene expression levels.ConclusionWhen adapted to biological hypotheses and combined with a promoter analysis, the method proposed here leads to improved models of the transcriptional regulatory networks.
| S-EPMC1892092 | biostudies-literature
Project description:We use the polynomial least squares method (PLSM), which allows us to compute analytical approximate polynomial solutions for a very general class of strongly nonlinear delay differential equations. The method is tested by computing approximate solutions for several applications including the pantograph equations and a nonlinear time-delay model from biology. The accuracy of the method is illustrated by a comparison with approximate solutions previously computed using other methods.
| S-EPMC4134833 | biostudies-other
Project description:Diffusion processes are governed by external triggers and internal dynamics in complex systems. Timely and cost-effective control of infectious disease spread critically relies on uncovering underlying diffusion mechanisms, which is challenging due to invisible infection pathways and time-evolving intensity of infection cases. Here, we propose a new diffusion framework for stochastic processes, which models disease spread across metapopulations by incorporating human mobility as topological pathways in a heterogeneous social system. We apply Bayesian inference with the stochastic Expectation-Maximization algorithm to quantify underlying diffusion dynamics in terms of exogeneity and endogeneity and estimate cross-regional infection flow based on Granger causality. The effectiveness of our proposed model is shown by using comprehensive simulation procedures (robustness tests with noisy data considering missing or delayed human case reporting in real situations) and by applying the model to real data from 15-y dengue outbreaks in Australia.
| S-EPMC6329989 | biostudies-literature