Modeling stochastic noise in gene regulatory systems

Unveiling the Dynamics behind Glioblastoma Multiforme Single-Cell Data Heterogeneity

Glioblastoma Multiforme is a brain tumor distinguished by its aggressiveness. We suggested that this aggressiveness leads single-cell RNA-sequence data (scRNA-seq) to span a representative portion of the cancer attractors domain. This conjecture allowed us to interpret the scRNA-seq heterogeneity as reflecting a representative trajectory within the attractor's domain. We considered factors such as genomic instability to characterize the cancer dynamics through stochastic fixed points. The fixed points were derived from centroids obtained through various clustering methods to verify our method sensitivity. This methodological foundation is based upon sample and time average equivalence, assigning an interpretative value to the data cluster centroids and supporting parameters estimation. We used stochastic simulations to reproduce the dynamics, and our results showed an alignment between experimental and simulated dataset centroids. We also computed the Waddington landscape, which provided a visual framework for validating the centroids and standard deviations as characterizations of cancer attractors. Additionally, we examined the stability and transitions between attractors and revealed a potential interplay between subtypes. These transitions might be related to cancer recurrence and progression, connecting the molecular mechanisms of cancer heterogeneity with statistical properties of gene expression dynamics. Our work advances the modeling of gene expression dynamics and paves the way for personalized therapeutic interventions.

Feasibility and reliability of sequential logic with gene regulatory networks

Gene regulatory networks exhibiting Boolean behaviour, e.g. AND, OR or XOR, have been routinely designed for years. However, achieving more sophisticated functions, such as control or computation, usually requires sequential circuits or so-called state machines. For such a circuit, outputs depend both on inputs and the current state of the system. Although it is still possible to design such circuits by analogy with digital electronics, some particularities of biology make the task trickier. The impact of two of them, namely the stochasticity of biological processes and the inhomogeneity in the response of regulation mechanisms, are assessed in this paper. Numerical simulations performed in two use cases point out high risks of malfunctions even for designed GRNs functional from a theoretical point of view. Several solutions to improve reliability of such systems are also discussed.

External Noise and External Signal Induced Transition of Gene Switch and Coherence Resonance in the Genetic Regulatory System

The transition of gene switch induced by external noises (multiplicative external noise and additive external noise) and external signals is investigated in the genetic regulatory system. Results show that the state-to-state transition of gene switch as well as resonant behaviors, such as the explicit coherence resonance (ECR), implicit coherence resonance (ICR) and control parameter coherence biresonance (CPCBR), can appear when noises are injected into the genetic regulatory system. The ECR is increased with the increase of the control parameter value when starting from the supercritical Hopf bifurcation parameter point, and there exists a critical control parameter value for the occurrence of ECR. However, the ICR is decreased as the control parameter value is increased when starting from the subcritical Hopf bifurcation point. In particular, the coherence of ECR is higher and more sensitive to noise than that of ICR. When an external signal is introduced into the system, the enhancement or suppression of the CPCBR and the number of peaks strongly depend on the frequency and amplitude of the external signal. Furthermore, the gene regulation system can selectively enhance or decrease the noise-induced oscillation signals at preferred frequency and amplitude of an external signal.

Splitting nodes and linking channels: A method for assembling biocircuits from stochastic elementary units

Akin to electric circuits, we construct biocircuits that are manipulated by cutting and assembling channels through which stochastic information flows. This diagrammatic manipulation allows us to create a method which constructs networks by joining building blocks selected so that (a) they cover only basic processes; (b) it is scalable to large networks; (c) the mean and variance-covariance from the Pauli master equation form a closed system; and (d) given the initial probability distribution, no special boundary conditions are necessary to solve the master equation. The method aims to help with both designing new synthetic signaling pathways and quantifying naturally existing regulatory networks.

Stochastic modeling suggests that noise reduces differentiation efficiency by inducing a heterogeneous drug response in glioma differentiation therapy

Background

Glioma differentiation therapy is a novel strategy that has been used to induce glioma cells to differentiate into glia-like cells. Although some advances in experimental methods for exploring the molecular mechanisms involved in differentiation therapy have been made, a model-based comprehensive analysis is still needed to understand these differentiation mechanisms and improve the effects of anti-cancer therapeutics. This type of analysis becomes necessary in stochastic cases for two main reasons: stochastic noise inherently exists in signal transduction and phenotypic regulation during targeted therapy and chemotherapy, and the relationship between this noise and drug efficacy in differentiation therapy is largely unknown.

Results

In this study, we developed both an additive noise model and a Chemical-Langenvin-Equation model for the signaling pathways involved in glioma differentiation therapy to investigate the functional role of noise in the drug response. Our model analysis revealed an ultrasensitive mechanism of cyclin D1 degradation that controls the glioma differentiation induced by the cAMP inducer cholera toxin (CT). The role of cyclin D1 degradation in human glioblastoma cell differentiation was then experimentally verified. Our stochastic simulation demonstrated that noise not only renders some glioma cells insensitive to cyclin D1 degradation during drug treatment but also induce heterogeneous differentiation responses among individual glioma cells by modulating the ultrasensitive response of cyclin D1. As such, the noise can reduce the differentiation efficiency in drug-treated glioma cells, which was verified by the decreased evolution of differentiation potential, which quantified the impact of noise on the dynamics of the drug-treated glioma cell population.

Conclusion

Our results demonstrated that targeting the noise-induced dynamics of cyclin D1 during glioma differentiation therapy can increase anti-glioma effects, implying that noise is a considerable factor in assessing and optimizing anti-cancer drug interventions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-016-0316-x) contains supplementary material, which is available to authorized users.

Fuzzy Stochastic Petri Nets for Modeling Biological Systems with Uncertain Kinetic Parameters

Stochastic Petri nets (SPNs) have been widely used to model randomness which is an inherent feature of biological systems. However, for many biological systems, some kinetic parameters may be uncertain due to incomplete, vague or missing kinetic data (often called fuzzy uncertainty), or naturally vary, e.g., between different individuals, experimental conditions, etc. (often called variability), which has prevented a wider application of SPNs that require accurate parameters. Considering the strength of fuzzy sets to deal with uncertain information, we apply a specific type of stochastic Petri nets, fuzzy stochastic Petri nets (FSPNs), to model and analyze biological systems with uncertain kinetic parameters. FSPNs combine SPNs and fuzzy sets, thereby taking into account both randomness and fuzziness of biological systems. For a biological system, SPNs model the randomness, while fuzzy sets model kinetic parameters with fuzzy uncertainty or variability by associating each parameter with a fuzzy number instead of a crisp real value. We introduce a simulation-based analysis method for FSPNs to explore the uncertainties of outputs resulting from the uncertainties associated with input parameters, which works equally well for bounded and unbounded models. We illustrate our approach using a yeast polarization model having an infinite state space, which shows the appropriateness of FSPNs in combination with simulation-based analysis for modeling and analyzing biological systems with uncertain information.