Determining the bistability parameter ranges of artificially induced lac operon using the root locus method

Root locus-based stability analysis for biological systems

Background: The first objective for realizing and handling biological systems is to choose a suitable model prototype and then perform structure and parameter identification. Afterwards, a theoretical analysis is needed to understand the characteristics, abilities, and limitations of the underlying systems. Generalized Michaelis-Menten kinetics (MM) and S-systems are two well-known biochemical system theory-based models. Research on steady-state estimation of generalized MM systems is difficult because of their complex structure. Further, theoretical analysis of S-systems is still difficult because of the power-law structure, and even the estimation of steady states can be easily achieved via algebraic equations. Aim: We focus on how to flexibly use control technologies to perform deeper biological system analysis. Methods: For generalized MM systems, the root locus method (proposed by Walter R. Evans) is used to predict the direction and rate (flux) limitations of the reaction and to estimate the steady states and stability margins (relative stability). Mode analysis is additionally introduced to discuss the transient behavior and the setting time. For S-systems, the concept of root locus, mode analysis, and the converse theorem are used to predict the dynamic behavior, to estimate the setting time and to analyze the relative stability of systems. Theoretical results were examined via simulation in a Simulink/MATLAB environment. Results: Four kinds of small functional modules (a system with reversible MM kinetics, a system with a singular or nearly singular system matrix and systems with cascade or branch pathways) are used to describe the proposed strategies clearly. For the reversible MM kinetics system, we successfully predict the direction and the rate (flux) limitations of reactions and obtain the values of steady state and net flux. We observe that theoretically derived results are consistent with simulation results. Good prediction is observed ([Formula: see text]% accuracy). For the system with a (nearly) singular matrix, we demonstrate that the system is neither globally exponentially stable nor globally asymptotically stable but globally semistable. The system possesses an infinite gain margin (GM denoting how much the gain can increase before the system becomes unstable) regardless of how large or how small the values of independent variables are, but the setting time decreases and then increases or always decreases as the values of independent variables increase. For S-systems, we first demonstrate that the stability of S-systems can be determined by linearized systems via root loci, mode analysis, and block diagram-based simulation. The relevant S-systems possess infinite GM for the values of independent variables varying from zero to infinity, and the setting time increases as the values of independent variables increase. Furthermore, the branch pathway maintains oscillation until a steady state is reached, but the oscillation phenomenon does not exist in the cascade pathway because in this system, all of the root loci are located on real lines. The theoretical predictions of dynamic behavior for these two systems are consistent with the simulation results. This study provides a guideline describing how to choose suitable independent variables such that systems possess satisfactory performance for stability margins, setting time and dynamic behavior. Conclusion: The proposed root locus-based analysis can be applied to any kind of differential equation-based biological system. This research initiates a method to examine system dynamic behavior and to discuss operating principles.

Influence of far-red light coherence on the functional state of plants

The influence of the coherence of far-red (730 nm) light on the functional activity of plants was studied. Blackberry explants cultivated in vitro on an artificial nutrient medium served as a biological model. The explants were irradiated with light beams with different spatial and temporal coherence. The average cell size D was taken as the discrimination threshold for the coherence length L_{coh} and the correlation radius r_{cor}. The results of irradiation were judged by the length and number of shoots formed on each explant. The greatest photoinduced effect was observed when the conditions L_{coh}, r_{cor}>D were fulfilled, i.e., when the cell fit completely in the coherence volume of the light wave field. Significant differences in growth parameters were also observed in the variants of the experiment with a constant frequency spectrum of radiation (fixed L_{coh}), but different r_{cor}. It is concluded that the correlation properties of radiation affect photoregulatory processes.

Bifurcation analysis of bistable and oscillatory dynamics in biological networks using the root-locus method

Most of the biological systems including gene regulatory networks can be described well by ordinary differential equation models with rational non-linearities. These models are derived either based on the reaction kinetics or by curve fitting to experimental data. This study demonstrates the applicability of the root-locus-based bifurcation analysis method for studying the complex dynamics of such models. The effectiveness of the bifurcation analysis in determining the exact parameter regions in each of which the system shows a certain dynamical behaviour, such as bistability, oscillation, and asymptotically equilibrium dynamics is shown by considering two mostly studied gene regulatory networks, namely Gardner's genetic toggle switch and p53 gene network possessing two-phase (mono-stable/oscillation) dynamics.

Two-dimensional polynomial type canonical relaxation oscillator model for p53 dynamics

p53 network, which is responsible for DNA damage response of cells, exhibits three distinct qualitative behaviours; low state, oscillation and high state, which are associated with normal cell cycle progression, cell cycle arrest and apoptosis, respectively. The experimental studies demonstrate that these dynamics of p53 are due to the ATM and Wip1 interaction. This paper proposes a simple two-dimensional canonical relaxation oscillator model based on the identified topological structure of ATM and Wip1 interaction underlying these qualitative behaviours of p53 network. The model includes only polynomial terms that have the interpretability of known ATM and Wip1 interaction. The introduced model is useful for understanding relaxation oscillations in gene regulatory networks. Through mathematical analysis, we investigate the roles of ATM and Wip1 in forming of these three essential behaviours, and show that ATM and Wip1 constitute the core mechanism of p53 dynamics. In agreement with biological findings, we show that Wip1 degradation term is a highly sensitive parameter, possibly related to mutations. By perturbing the corresponding parameters, our model characterizes some mutations such as ATM deficiency and Wip1 overexpression. Finally, we provide intervention strategies considering our observation that Wip1 seems to be an important target to conduct therapies for these mutations.

Aggregation for Computing Multi-Modal Stationary Distributions in 1-D Gene Regulatory Networks

This paper proposes aggregation-based, three-stage algorithms to overcome the numerical problems encountered in computing stationary distributions and mean first passage times for multi-modal birth-death processes of large state space sizes. The considered birth-death processes which are defined by Chemical Master Equations are used in modeling stochastic behavior of gene regulatory networks. Computing stationary probabilities for a multi-modal distribution from Chemical Master Equations is subject to have numerical problems due to the probability values running out of the representation range of the standard programming languages with the increasing size of the state space. The aggregation is shown to provide a solution to this problem by analyzing first reduced size subsystems in isolation and then considering the transitions between these subsystems. The proposed algorithms are applied to study the bimodal behavior of the lac operon of E. coli described with a one-dimensional birth-death model. Thus, the determination of the entire parameter range of bimodality for the stochastic model of lac operon is achieved.