Stochastic Transcriptional Regulatory Systems with Time Delays: A Mean-Field Approximation

NONLINEAR DYNAMICS OF CELL CYCLES WITH STOCHASTIC MATHEMATICAL MODELS

Cell cycles are fundamental components of all living organisms and their systematic studies extend our knowledge about the interconnection between regulatory, metabolic, and signaling networks, and therefore open new opportunities for our ultimate efficient control of cellular processes for disease treatments, as well as for a wide variety of biomedical and biotechnological applications. In the study of cell cycles, nonlinear phenomena play a paramount role, in particular in those cases where the cellular dynamics is in the focus of attention. Quantification of this dynamics is a challenging task due to a wide range of parameters that require estimations and the presence of many stochastic effects. Based on the originally deterministic model, in this paper we develop a hierarchy of models that allow us to describe the nonlinear dynamics accounting for special events of cell cycles. First, we develop a model that takes into account fluctuations of relative concentrations of proteins during special events of cell cycles. Such fluctuations are induced by varying rates of relative concentrations of proteins and/or by relative concentrations of proteins themselves. As such fluctuations may be responsible for qualitative changes in the cell, we develop a new model that accounts for the effect of cellular dynamics on the cell cycle. Finally, we analyze numerically nonlinear effects in the cell cycle by constructing phase portraits based on the newly developed model and carry out a parametric sensitivity analysis in order to identify parameters for an efficient cell cycle control. The results of computational experiments demonstrate that the metabolic events in gene regulatory networks can qualitatively influence the dynamics of the cell cycle.

Modeling and analysis of the dynamic behavior of the XlnR regulon in Aspergillus niger

Background

In this paper the dynamics of the transcription-translation system for XlnR regulon in Aspergillus niger is modeled. The model is based on Hill regulation functions and uses ordinary differential equations. The network response to a trigger of D-xylose is considered and stability analysis is performed. The activating, repressive feedback, and the combined effect of the two feedbacks on the network behavior are analyzed.

Results

Simulation and systems analysis showed significant influence of activating and repressing feedback on metabolite expression profiles. The dynamics of the D-xylose input function has an important effect on the profiles of the individual metabolite concentrations. Variation of the time delay in the feedback loop has no significant effect on the pattern of the response. The stability and existence of oscillatory behavior depends on which proteins are involved in the feedback loop.

Conclusions

The dynamics in the regulation properties of the network are dictated mainly by the transcription and translation degradation rate parameters, and by the D-xylose consumption profile. This holds true with and without feedback in the network. Feedback was found to significantly influence the expression dynamics of genes and proteins. Feedback increases the metabolite abundance, changes the steady state values, alters the time trajectories and affects the response oscillatory behavior and stability conditions. The modeling approach provides insight into network behavioral dynamics particularly for small-sized networks. The analysis of the network dynamics has provided useful information for experimental design for future in vitro experimental work.

Planning interventions in biological networks

Modeling the dynamics of biological processes has recently become an important research topic in computational biology and systems engineering. One of the most important reasons to model a biological process is to enable high-throughput in-silico experiments that attempt to predict or intervene in the process. These experiments can help accelerate the design of therapies through their rapid and inexpensive replication and alteration. While some techniques exist for reasoning with biological processes, few take advantage of the flexible and scalable algorithms popular in AI research. In reasoning about interventions in biological processes, where scalability is crucial for feasible application, we apply AI planning-based search techniques and demonstrate their advantage over existing enumerative methods. We also present a novel formulation of intervention planning that relies on models that characterize and attempt to change the phenotype of a system. We study three biological systems: the yeast cell cycle, a model of the human aging process, and the Wnt5a network governing the metastasis of melanoma in humans. The contribution of our investigation is in demonstrating that: (i) prior approaches, based on dynamic programming, cannot scale as well as heuristic search, and (ii) the newly found scalability enables us to plan previously unknown sequences of interventions that reveal novel and biologically significant responses in the systems which are consistent with biological knowledge in the literature.

Quantitative model for inferring dynamic regulation of the tumour suppressor gene p53

BACKGROUND

The availability of various "omics" datasets creates a prospect of performing the study of genome-wide genetic regulatory networks. However, one of the major challenges of using mathematical models to infer genetic regulation from microarray datasets is the lack of information for protein concentrations and activities. Most of the previous researches were based on an assumption that the mRNA levels of a gene are consistent with its protein activities, though it is not always the case. Therefore, a more sophisticated modelling framework together with the corresponding inference methods is needed to accurately estimate genetic regulation from "omics" datasets.

RESULTS

This work developed a novel approach, which is based on a nonlinear mathematical model, to infer genetic regulation from microarray gene expression data. By using the p53 network as a test system, we used the nonlinear model to estimate the activities of transcription factor (TF) p53 from the expression levels of its target genes, and to identify the activation/inhibition status of p53 to its target genes. The predicted top 317 putative p53 target genes were supported by DNA sequence analysis. A comparison between our prediction and the other published predictions of p53 targets suggests that most of putative p53 targets may share a common depleted or enriched sequence signal on their upstream non-coding region.

CONCLUSIONS

The proposed quantitative model can not only be used to infer the regulatory relationship between TF and its down-stream genes, but also be applied to estimate the protein activities of TF from the expression levels of its target genes.

Stochastic models and numerical algorithms for a class of regulatory gene networks

Regulatory gene networks contain generic modules, like those involving feedback loops, which are essential for the regulation of many biological functions (Guido et al. in Nature 439:856-860, 2006). We consider a class of self-regulated genes which are the building blocks of many regulatory gene networks, and study the steady-state distribution of the associated Gillespie algorithm by providing efficient numerical algorithms. We also study a regulatory gene network of interest in gene therapy, using mean-field models with time delays. Convergence of the related time-nonhomogeneous Markov chain is established for a class of linear catalytic networks with feedback loops.

Steady-state expression of self-regulated genes

MOTIVATION

Regulatory gene networks contain generic modules such as feedback loops that are essential for the regulation of many biological functions. The study of the stochastic mechanisms of gene regulation is instrumental for the understanding of how cells maintain their expression at levels commensurate with their biological role, as well as to engineer gene expression switches of appropriate behavior. The lack of precise knowledge on the steady-state distribution of gene expression requires the use of Gillespie algorithms and Monte-Carlo approximations.

METHODOLOGY

In this study, we provide new exact formulas and efficient numerical algorithms for computing/modeling the steady-state of a class of self-regulated genes, and we use it to model/compute the stochastic expression of a gene of interest in an engineered network introduced in mammalian cells. The behavior of the genetic network is then analyzed experimentally in living cells.

RESULTS

Stochastic models often reveal counter-intuitive experimental behaviors, and we find that this genetic architecture displays a unimodal behavior in mammalian cells, which was unexpected given its known bimodal response in unicellular organisms. We provide a molecular rationale for this behavior, and we implement it in the mathematical picture to explain the experimental results obtained from this network.