 Research
 Open access
 Published:
Stability analysis of genetic regulatory network with additive noises
BMC Genomics volumeÂ 9, ArticleÂ number:Â S21 (2008)
Abstract
Background
Genetic regulatory networks (GRN) can be described by differential equations with SUM logic which has been found in many natural systems. Identification of the network components and transcriptional rates are critical to the output behavior of the system. Though transcriptional rates cannot be measured in vivo, biologists have shown that they are alterable through artificial factors in vitro.
Results
This study presents the theoretical research work on a novel nonlinear control and stability analysis of genetic regulatory networks. The proposed control scheme can drive the genetic regulatory network to desired levels by adjusting transcriptional rates. Asymptotic stability proof is conducted with Lyapunov argument for both noisefree and additive noises cases. Computer simulation results show the effectiveness of the control design and robustness of the regulation scheme with additive noises.
Conclusions
With the knowledge of interaction between transcriptional factors and gene products, the research results can be applied in the design of modelbased experiments to regulate gene expression profiles.
Background
Genetic networks regulate sophisticated biological functions by interacting genes and proteins and support homeostasis in metabolism and coordinate events during the developmental program. Research on stability analysis and regulation/control of these genetic networks are particularly important. Pioneer experimental studies in construction of genetic networks to manipulate protein levels or even to construct gene circuits with repressor functions have been carried out [1â€“9]. These experiments have demonstrated interesting properties of GRNs of E. coli in the presence of specified repressors. With different repressors, these GRNs include single or multiple interactions between genes and proteins. In a single gene regulatory network [1], the negative feedback that is integrated in the system decreases celltocell fluctuations in protein concentration measurements. Distribution of the regulated protein concentration is proportional to the degradation rate of the gene network. In a twogene regulation network [3], bistability is shown by coupling two proteins with negative regulation in the synthesis of each other. Stability analysis of this bigene regulation network is also presented based on parameter bifurcations. The significance of this experiment is that the transition between two stable states of the GRN is much sharper with respect to the intracellular stimuli, i.e., performance can be adjusted by special input. In a trigene regulatory model [2], three genes are regulated with cyclic repressibility. The system exhibited selfsustained oscillations over the entire growth phase of the host E. coli cells for certain biochemical parameters. These biological experimental results have shown that genetic networks can be regulated by a scheme of local promotor control, i.e. the number, type and placement of regulatory protein binding sites. However, quantitative analysis of the regulation function has not been studied.
The objective of our current study is to develop a mathematical model of the trigene regulation network and extend the theoretical stability analysis to the case with measurement noises. A novel control scheme is proposed to change the state of a genetic regulatory network by adjusting transcriptional rates. The proposed control scheme can provide biologists useful design techniques in modelbased experiments to predict protein levels in genetic regulatory networks by adjusting specific regulatory factors. We will use the biological scheme of adjusting regulatory factors reported in research articles [9â€“14]. These regulatory factors will be adjusted based on errors between the measured gene products, mRNA and protein, concentration levels and their desired values to regulate the gene expressions profiles. The proposed control is based on the concentrations of m RNAs and proteins which can be easily measured by current molecular biology techniques. Therefore, it's useful when some of the transcription rates are unmeasurable.
We will introduce the general model of GRNs to explain the work, present the control design and the stability analysis , and then give a simulation example and conclusion.
Mathematical model of genetic regulatory networks
A gene regulatory network for a eukaryote is shown in figure 1. Potential inputs of the system include a variety of developmental and environmental stimuli. System outputs are the synthesized proteins. Inputs activate a complex chain of intracellular reactions that activate a regulatory molecule, transcriptional factor, to translocate from the cytoplasm into the nucleus. In the nucleus, active transcription factors recognize a specific segment of DNA, termed a promotor. The promotor informs RNA polymerase, which binds closely to the promotor, where to start transcribing genetic information on DNA to mRNA. The mRNA molecules then leave the cell nucleus and enter the cytoplasm where proteins are synthesized in the presence of transfer RNAs. When the translated protein is capable of interacting with its own or other gene's transcription factor (denoted by the dashed arrow in figure 1, a regulatory or feedback loop is formed. Such transcriptional regulation is the typical method utilized by cells to control gene expression [9, 15]. Feedback can occur in either a positive (activator) or a negative (repressor) direction. In a gene regulatory network, natural regulatory molecules are either activated transcription factors or proteins that activate transcription factors. Gene expression is very sensitive to changes in the composition of regulatory factors, which limits attempts to control gene expression. If artificial regulatory factors [10â€“14] can be used to change the activity of the transcription factors and the affinity of DNA binding, we can affect the transcriptional rate in GRNs and control the output behavior of the GRNs. Thus, in order to derive novel agents that regulate gene expressions, we must first understand GRN output behavior.
Two types of mathematical models have been developed to understand the working scheme of genetic networks: (1) Boolean networks and (2) sets of differential equations [9, 16â€“19]. Boolean network model describes the expression of each gene with two states: ON or OFF, and the state of a gene is determined by a Boolean function of the states of other related genes. The differential equation model describes concentrations of mRNAs and proteins as continuous functions, which can provide more accurate and detailed dynamic information. Since genetic networks are high dimensional and nonlinear in nature, it is logical to consider such genetic networks from a nonlinear dynamic viewpoint, i.e., nonlinear differential equations. From a control point of view, main purpose of our mathematical model is to predict and manipulate the dynamic behaviors by analyzing available measurements. Researches of control design for Boolean networks have been carried out [20, 21], while for continuous differential model, there are comparatively few literature references available [22].
In this section, the studied genetic network model is described by nonlinear differential equations
where m_{ i }(t) and p_{ j }(t) are the concentrations of m RNAs and proteins in a genetic regulatory network, i = 1, 2, â€¦ , n, and j = 1, 2, â€¦ , m. Parameters k_{ mi } and k_{ pj } are degradation rates of the m RNAs and proteins, k_{ dji } are assumed to be constants describing the link between m RNAs and proteins. If the jth protein regulates the ith gene as illustrated by the dashed arrow in figure 1, there is a regulatory link in the network. Such regulation effects of proteins to gene expressions are described by wrapped nonlinear terms b_{ i }(p_{1}(t), p_{2}(t), â€¦ , p_{ m }(t)), which are nonlinear functions of p_{1}(t), p_{2}(t), â€¦ , p_{ m }(t), and can be described by SUM logic (additive effect) of each protein to the specified gene. Detailed description of the relationship between Î±_{ i }s, b_{ i }s and transcriptional factors have been previously published [9, 17, 18]. Adjustment of the regulatory factors will change parameters Î±_{ ij } in the mathematical depiction written in equation (2).
where case 1 represents that transcription factor (protein) j is an activator of gene i, and in case 2, transcription factor j is a repressor of gene i. Parameters H in equation (2) is the Hill coefficient, Î² is a positive constant, and Î±_{ ij } is the dimensionless transcription rate of transcription factor j to gene i. Equation (1) can be rewritten as the following form by applying (2)
where m(t)=[{m}_{1}(t),\xe2\u20ac\xa6{m}_{i}(t)\xe2\u20ac\xa6{m}_{n}(t)]\xe2\u02c6\u02c6{\mathrm{\xe2\u201e\u0153}}^{n},p(t)=[{p}_{1}(t),\xe2\u20ac\xa6{p}_{j}(t)\xe2\u20ac\xa6{p}_{m}(t)]\xe2\u02c6\u02c6{\mathrm{\xe2\u201e\u0153}}^{m},A=diag\{\xe2\u02c6\u2019{k}_{{m}_{i}}\}\xe2\u02c6\u02c6{\mathrm{\xe2\u201e\u0153}}^{n\xc3\u2014n},C=diag\{\xe2\u02c6\u2019{k}_{{p}_{j}}\}\xe2\u02c6\u02c6{\mathrm{\xe2\u201e\u0153}}^{m\xc3\u2014m},D=diag\left\{{k}_{dji}\right\}\xe2\u02c6\u02c6{\mathrm{\xe2\u201e\u0153}}^{m\xc3\u2014n}. Vector {l}_{0}\xe2\u02c6\u02c6{\mathrm{\xe2\u201e\u0153}}^{n} is determined by
within Î©_{ ij }, which is a set of repressors for gene i. Matrix B\xe2\u02c6\u02c6{\mathrm{\xe2\u201e\u0153}}^{n\xc3\u2014m} describes the interconnection of transcription factors and genes in the network with element Î±_{ ij },
and g\left(p\left(t\right)\right)={\left[g\left({p}_{1}\left(t\right)\right)\xe2\u20ac\xa6g\left({p}_{j}\left(t\right)\right)\xe2\u20ac\xa6g\left({p}_{n}\left(t\right)\right)\right]}^{T}\xe2\u02c6\u02c6{\mathrm{\xe2\u201e\u0153}}^{m},g\left({p}_{j}\left(t\right)\right)=\frac{{\left({p}_{j}\left(t\right)/\mathrm{\xce\xb2}\right)}^{H}}{1+{\left({p}_{j}\left(t\right)/\mathrm{\xce\xb2}\right)}^{H}}. Some special properties of the nonlinear function g(p(t)) should be pointed out:

g(p_{ j }) > 0 always holds with p_{ j } > 0 and the equal sign holds only when p_{ j } =0;

g(p_{ j }) is a monotonic increasing function, i.e. \frac{{\xe2\u02c6\u201a}_{g}({p}_{j})}{\xe2\u02c6\u201a{p}_{j}}\xe2\u2030\yen 0;

g(p_{ j }) satisfies a sector condition, which can be obtained from the meanvalue theorem, g\left({e}_{j}+{p}_{j}\right)\xe2\u02c6\u2019g\left({p}_{j}\right)=\frac{\xe2\u02c6\u201ag}{\xe2\u02c6\u201a{p}_{j}}{p}_{j}=\mathrm{\xce\xbe}\xe2\u2039\dots {e}_{j},\mathrm{\xce\xbe}\xe2\u02c6\u02c6\left[{p}_{j},{p}_{j}+{e}_{j}\right], and thus, 0\xe2\u2030\xa4{e}_{j}[g\left({e}_{j}+{p}_{j}\right)\xe2\u02c6\u2019g\left({p}_{j}\right)]\xe2\u2030\xa4{\mathrm{\xce\xb3}}_{j}{e}_{j}^{2}, where {\mathrm{\xce\xb3}}_{j}=\mathrm{max}\left\{\frac{H{\left({p}_{j}/\mathrm{\xce\xb2}\right)}^{H\xe2\u02c6\u20191}}{\mathrm{\xce\xb2}{\left(1+{\left({p}_{j}/\mathrm{\xce\xb2}\right)}^{H}\right)}^{2}}{p}_{j}\xe2\u02c6\u02c6\left[{p}_{j},{p}_{j}+e\right]\right\}. The vector case of the above equation can be expressed as
0\xe2\u2030\xa4{e}^{T}\left[g\left(e+p\right)\xe2\u02c6\u2019g\left(p\right)\right]\xe2\u2030\xa4\mathrm{\xce\xb3}{e}^{T}e,(5)
where Î³ = max {Î³_{1}, Î³_{ j } â€¦, Î³_{ m }}.
Combine the state m(t), p(t) together, the state space model for control design and stability analysis is written as follows:
where {x}^{T}=\left[{m}^{T}{p}^{T}\right],{l}^{T}=\left[{l}_{0}^{T}0\right],x,l\xe2\u02c6\u02c6{\mathrm{\xe2\u201e\u0153}}^{n+m},{C}_{c}=[{0}_{m\xc3\u2014n}{I}_{m\xc3\u2014m}]\xe2\u02c6\u02c6{\mathrm{\xe2\u201e\u0153}}^{m\xc3\u2014(n+m)}
Methods: adaptive control of genetic networks in a noisefree case
Nonlinear adaptive control has been applied to many systems to improve performance. The control objective here is to drive the current state of genetic networks to desired values m^{*} and p^{*}. In order to make the control design biologically meaningful, the following assumptions are introduced.
Assumption 1: Transcription rates of mRNAs in the studied genetic network are adjustable.
This assumption guarantees the possibility of adjusting the transcription rates and drive the current states to the desired values.
Assumption 2: State x^{*T} = [m^{*T} p^{*T}] is the stable steady state generated by the same genetic network in equation (3) with desired transcription rates.
This assumption guarantees that m^{*} and p^{*} are achievable and has biological meaning for a real genetic network. The stability also guarantees that once the state is driven to x^{*}, it will stay there.
The control design proposed here includes two parts: i) stability analysis of x^{*} generated by the genetic network with the desired degradation, transcription rate, the given sector parameter Î³; ii) control design to drive the current states to the desired x^{*}. Assume x^{*} is generated by
where {B}_{c}^{*} is the desired transcription rates. Stability of the genetic network means: (1) State x^{*} is the equilibrium of equation (7), i.e. {A}_{c}{x}^{*}\xe2\u02c6\u2019{B}_{c}^{*}g\left({C}_{c}{x}^{*}\right)+l=0; (2) Starting at any initial states x_{0} close to x^{*}, the trajectory will converge to x^{*} as time goes to infinity. If we define e^{*} = x âˆ’ x^{*}, the error e â†’ 0 as t â†’ âˆž.
Lemma 1: The system (7) is globally exponentially stable if the transfer function matrix G(s) = C_{ c }(sI âˆ’ A_{ C })^{âˆ’1}{B}_{c}^{*}is strictly positive real, i.e. with the controllable and observable pair\{{A}_{c},{B}_{c}^{*},{C}_{c}\}, there exist positive matrices P and Q, P = P^{T}, such that
Proof: With constant x^{*}, the error dynamics can be expressed as
Choosing the Lyapunov function {V}_{1}(e)=\frac{1}{2}{e}^{T}Pe, its time derivative {\stackrel{\xcb\u2122}{V}}_{1}(e) is obtained by the following equations.
which is negative definite since (C_{ c }e)^{T} {g[C_{ c }e + C_{ c }x^{*})]â€”g(C_{ c }x^{*})} > 0 by applying equation (5). Thus, the error dynamics will converge to zero exponentially by satisfying the strict positive real condition given as equation (8) on the system parameters.
Remark 1: Lemma 1 gives a sufficient condition on parameter settings for a stable genetic network. The stability is determined by system parameters and the sector condition of the nonlinear feedback function.
Lemma 1 provides an easy way to check the stability of a genetic network, since it is easy to get the transfer function for linear time invariant system and check whether the transfer function is strictly positive real or not.
Remark 2: While applying the small gain theorem, with the consideration of the sector parameter Î³, we can get the sufficient and necessary condition: {\mathrm{sup}}_{\mathrm{\xcf\u2030}}\leftG\left(j\mathrm{\xcf\u2030}\right)\right<1/\mathrm{\xce\xb3},\mathrm{\xcf\u2030}\xe2\u02c6\u02c6\mathrm{\xe2\u201e\u0153}. Since the necessary condition is not involved in the following control design, we ignore the proof here. Related information of it can be found in [23â€“25].
Based on the above stability analysis, the following theorem gives the control design that drives the current state to the desired state by adjusting the transcription rates of the system.
Theorem 1: Assume transcription rates of gene i in system (7) can be adjusted by the control law
where i = 1, 2, â€¦ , n, j = 1, 2, â€¦ , m, e_{ k } is the k th element in the error vector between the current state and the desired state, i.e., e = xâ€”x^{*}, g_{ j } = g(p_{ j }(t)), P_{ ij } is the element in the positive definite matrix P defined in equation (8). The system (6) will converge to the desired state x^{*} asymptotically as time goes to infinity.
Proof: From equations (6, 7), the error dynamics is
Choose Lyapunov candidate V\left(e,{\stackrel{\xcb\u0153}{\mathrm{\xce\pm}}}_{ij}\right)=\frac{1}{2}{e}^{T}Pe+{\displaystyle {\xe2\u02c6\u2018}_{i,j=1}^{i=n,j=m}\frac{1}{2{k}_{ij}}{\stackrel{\xcb\u0153}{\mathrm{\xce\pm}}}_{ij}^{2}}, where {\stackrel{\xcb\u0153}{\mathrm{\xce\pm}}}_{ij}\text{s} are defined as {\stackrel{\xcb\u0153}{\mathrm{\xce\pm}}}_{ij}={B}_{c}\xe2\u02c6\u2019{B}_{c}^{*}={\mathrm{\xce\pm}}_{ij}\xe2\u02c6\u2019{\mathrm{\xce\pm}}_{ij}^{*}. Recall that the bottom block with m Ã— n dimension in both matrices are zeroes, {\stackrel{\xcb\u0153}{\mathrm{\xce\pm}}}_{ij} is only related to the transcription rates in the system. Thus, the time derivative of the Lyapunov function can be derived as follows.
with the adaption control chosen as {\stackrel{\xcb\u2122}{\stackrel{\xcb\u0153}{\mathrm{\xce\pm}}}}_{ij}=\xe2\u02c6\u2019{\stackrel{\xcb\u2122}{\mathrm{\xce\pm}}}_{ij}={k}_{ij}{\displaystyle {\xe2\u02c6\u2018}_{i,j,k=1}^{i=n,j=m,k=n+m}{e}_{k}{P}_{ki}{g}_{j}}.
From the above Lyapunov argument, with a positive definite V and negative definite \stackrel{\xcb\u2122}{V} all the errors decrease to zero. This concludes the proof that the tracking error of the genetic network from current state to the desired state is globally asymptotically stable with the adaptive control given in (11), i.e., the states are driven to desired levels. As the tracking errors converge to zero, the adjustable transcription rates Î±_{ ij } will converge to constants.
Boundedness with additive noises
When signals are sensed, signal distortion, transmission delay and noise are unavoidable. In this section an additive measurement noise will be considered, and the distorted measurements of mRNA and protein concentration levels will be used in the adaption law.
Considering the aforementioned system in equation (6) with additive noises
where f:{\mathrm{\xe2\u201e\u0153}}^{+}\xc3\u2014{\mathrm{\xe2\u201e\u0153}}^{n}\xe2\u2020\u2019{\mathrm{\xe2\u201e\u0153}}^{n},d:{\mathrm{\xe2\u201e\u0153}}^{+}\xc3\u2014{\mathrm{\xe2\u201e\u0153}}^{n}\xe2\u2020\u2019{\mathrm{\xe2\u201e\u0153}}^{n} are piecewise continuous in t, the systems can be viewed as the nominal system
with perturbation term d(t, x). The variable u is a function of x and is omitted here for simplicity.
Assume the nominal system has an equilibrium point at the origin, if d(t, x) = 0 as x = 0, then the origin is still an equilibrium point with the disturbance d(t, x). However, in most cases, the disturbance does not satisfy this assumption, i.e. d(t, 0) â‰ 0. In this case, the origin is not an equilibrium point and no conclusion can be drawn about stability of the origin. The following Theorem shows that the best result we can expect is the uniform boundedness of the disturbed system when the origin is exponentially stable.
Theorem 2: Consider the system and adaptation laws given by equations (6), and (11), the system is uniformly bounded with disturbance \xe2\u20ac\u2013d\left(t,x\right)\xe2\u20ac\u2013\xe2\u2030\xa4\mathrm{\xce\xb4} where Î´ is a positive constant. The upper bound b_{ u } and lower bound b_{ l } of the system is determined by
where 0 < Î¼ < 1, respectively. Î±_{ i }s are positive constants defined as follows.
proof From the proof of Theorem 1, the error dynamics of the whole system is exponentially stable at the origin. Thus, there exist a Lyapunov function V(\mathrm{\xce\xa6}),\mathrm{\xce\xa6}={[e,{\stackrel{\xcb\u0153}{\mathrm{\xce\pm}}}_{ij}]}^{T} satisfies the following condition:
With the additive noise disturbance (t, Î¦) considered, the time derivative of V_{2} is rewritten as
for any \xe2\u20ac\u2013\mathrm{\xce\xa6}\xe2\u20ac\u2013\xe2\u2030\yen \frac{\mathrm{\xce\xb4}{\mathrm{\xce\pm}}_{4}}{\mathrm{\xce\xbc}{\mathrm{\xce\pm}}_{3}}. This gives the lower bound as bl=\frac{\mathrm{\xce\xb4}{\mathrm{\xce\pm}}_{4}}{\mathrm{\xce\xbc}{\mathrm{\xce\pm}}_{3}}. From BIBO Theorem from the reference [26], the upper bound can be determined as bu=\frac{{\mathrm{\xce\pm}}_{4}\mathrm{\xce\xb4}}{{\mathrm{\xce\pm}}_{3}\mathrm{\xce\xbc}}\sqrt{\frac{{\mathrm{\xce\pm}}_{2}}{{\mathrm{\xce\pm}}_{1}}}.
From theorem 2, it can be shown that with bounded noises, both the state estimation errors and adaptation control are bounded. With disturbance bound Î´ â†’ 0, differences between the adaptation control and the desired constant transcriptional rates shrink to zero also.
Simulation results and discussion
To show the effectiveness of the proposed control, an example is simulated with Matlab. We consider the dynamics of the trigene network described in [2].
where i sequence is lacl, tetR, cl and j sequence is cl, lacl, tetR, k_{ mi } = 1, k_{ pi } = 1, n = 2 and Î± is the control parameter. With parameters chosen as Î±_{ i } = 2.5, the genetic network has a globally exponential equilibrium at {m}_{i}^{*}={p}_{i}^{*}=1.1150, which is the steady state. Evolution trajectories of m RNA and protein concentration levels are shown in figure 2.
The control objective is to alter the current state x = [3, 2, 1, 1, 3, 2]^{T} to the desired states x^{*} by adjusting transcription rates Î±_{i}. Assume there is no regulation at the initial state, the control law is chosen to be {\stackrel{\xcb\u2122}{\mathrm{\xce\pm}}}_{i}=\xe2\u02c6\u201910{e}_{i}{g}_{j}. Evolution trajectories of m_{ i } and p_{ i } with control design are shown in figure 3. The control variables are shown in figure 4. It can be seen that the steady states converge to x^{*} and the transcription rates converge to the desired value 2.5 in the simulation. The simulation also shows that the convergence rate of the state and the parameter with the proposed control is faster than the original genetic network.
In order to show the robust of the control scheme, an additive white noise is injected to the measurements of the mRNA and protein concentration levels. Magnitude of the noises are bounded by 0.5. All parameter settings are the same as the ideal case illustrated above. Simulation results of the regulated mRNA and protein concentration levels are shown in figure 5 and the control of the transcription rates are shown in figure 6.
Conclusions and future research
In the current study, a general regulatory network is presented and a new control scheme is proposed. This scheme controls the current states of the genetic network to desired values. We have obtained global convergence of the tracking error by online adjustment of the transcription rates. We have applied Lyapunov argument to the convergence analysis of the states and boundedness of the transcription rates. A trigene regulatory is simulated with the proposed algorithm. Effectiveness of the proposed controller design is verified by Matlab simulation for noisefree measurement and bounded noises.
In this research, interaction between the proteins and mRNAs are described by Hill functions with SUM logic. More detailed studies need to be carried out to determine the structure of a real GRN and extend the control scheme to more general cases. In addition, the estimation algorithms are developed with continuous measurers of the state. In real biological systems, such measures will be collected at specific time points. Estimation error caused by sampling effects should also be considered. Together, our results demonstrated that GRN could be regulated by artificially changing the transcription rates, to approach the desired gene expression levels.
References
Becskei A, Seraphin B, Serrano L: Positive feedback in eukaryotic gene networks: cell differentiation by graded to binary response conversion. EMBO J. 2001, 20: 25282535. 10.1093/emboj/20.10.2528.
Elowitz MB, Leibler S: A synthetic oscillatory network of transcriptional regulators. Nature. 2000, 403 (6767): 335338. 10.1038/35002125.
Gardner TS, Cantor CR, Collins JJ: Construction of a genetic toggle switch in Escherichia coli. Nature. 2000, 403 (6767): 339342. 10.1038/35002131.
Hasty J, Issacs F: Designer gene networks: Towards fundamental cellular control. Chaos: An Interdisciplinary Journal of Nonlinear Science. 2001, 11: 207220. 10.1063/1.1345702.
Hasty J, McMillen D, Isaacs F, Collins JJ: Computational studies of gene regulatory networks: in numero molecular biology. Nature Reviews Genetics. 2001, 2 (4): 26810.1038/35066056.
Hasty J, McMillen D, Collins JJ: Engineered gene circuits. Nature. 2002, 420 (6912): 22410.1038/nature01257.
Hasty J, Pradines J, Dolnik M, Collins JJ: Noisebased switches and amplifiers for gene expression. PNAS. 2000, 97 (5): 20752080. 10.1073/pnas.040411297.
McAdams H, Arkin A: Stochastic mechanisms in gene expression. PNAS. 1997, 94 (3): 814819. 10.1073/pnas.94.3.814.
Wolf DM, Eeckman FH: On the relationship between genomic regulatory element organization and gene regulatory dynamics. Journal of Theoretical Biology. 1998, 195 (2): 167186. 10.1006/jtbi.1998.0790.
Ansari AZ: Fingers reach for the genome. Nat. Biotechnol. 2003, 21: 242243. 10.1038/nbt0303242.
Beerli RR, Barbas III CF: Engineering polydactyl zincfinger transcription factors. Nat. Biotechnol. 2002, 20: 135141. 10.1038/nbt0202135.
Belshaw PJ, Ho SN, Crabtree GR, Schreiber SL: Controlling protein association and subcellular localization with a synthetic ligand that induces heterodimerization of proteins. Proc Natl Acad Sci U S A. 1996, 93: 46044607. 10.1073/pnas.93.10.4604.
Park KS, Lee DK, Lee H, Lee Y, Jang YS, Kim YH, Yang HY, Lee SI, Seol W, Kim JS: Phenotypic alteration of eukaryotic cells using randomized libraries of artificial transcription factors. Nat. Biotechnol. 2003, 21: 12081214. 10.1038/nbt868.
Schmitz K, Schepers U: Polymides as artificial transcription factors: novel tools for molecular medcine?. Angewandte Chemistry International Edition. 2004, 43: 24722475. 10.1002/anie.200301745.
Dickson AJBWaRW R: Genetic regulation: the lac control region. Science. 1975, 187: 2735. 10.1126/science.1088926.
de Jong H: Modeling and simulation of genetic regulatory systems: A literature review. J Comput Biol. 2002, 9: 67103. 10.1089/10665270252833208.
Keller AD: Model genetic circuits encoding autoregulatory transcription factors. Journal of Theoretical Biology. 1995, 172 (2): 169185. 10.1006/jtbi.1995.0014.
Likhoshvai MY VA, Fadeev S: Relationship between a gene network graph and qualitative modes of its functioning. Molecular Biology. 2001, 35 (6): 926932. 10.1023/A:1013206906557.
Smolen P, Baxter DA, Byrne JH: Mathematical modeling of gene networks. Neuron. 2000, 26 (3): 567580. 10.1016/S08966273(00)811940.
External control in Markovian genetic regulatory networks. Machine Learning. 2003, 52: 169191. 10.1023/A:1023909812213. Aniruddha Datta MLBERD Ashish Choudhary
External control in Markovian genetic regulatory networks: the imperfect information case. Bioinformatics. 2004, 20 (6): 924930. 10.1093/bioinformatics/bth008. Aniruddha Datta MLBERD Ashish Choudhary
Chen L, Aihara K: Stability of genetic regulatory networks with time delay. Circuits and Systems I: Fundamental Theory and Applications, IEEE Transactions on [see also Circuits and Systems I: Regular Papers, IEEE Transactions on]. 2002, 49 (5): 602608. 10.1109/TCSI.2002.1001949. [10577122].
Khalil H: Nonlinear Systems. 2002, Prentice Hall, Upper Saddle River, NJ
Kristic KM, Kanellahopoulos I, Kokotovic P: Nonlinear and Adaptive Control Design. 1995, John Wiley & Sons
Vidyasagar M: Nonlinear Systems Analysis. 2002, Prentice Hall, Englewood Cliffs, NJ
Qu Z: Adaptive and robust controls of uncertain systems with nonlinear parameterization. Automatic Control, IEEE Transactions on. 2003, 48 (10): 18171824. 10.1109/TAC.2003.817931. [00189286].
Acknowledgements
The authors acknowledge support from the NSF (0649172YJ) and NIH (R01 HL75360 MLL). The authors also appreciate Dr. Yufei Huang for his helpful discussion.
This article has been published as part of BMC Genomics Volume 9 Supplement 1, 2008: The 2007 International Conference on Bioinformatics & Computational Biology (BIOCOMP'07). The full contents of the supplement are available online at http://www.biomedcentral.com/14712164/9?issue=S1.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
This study is carried out by the authors together. Dr. Lindsey contributes to the biology background and the discussion part. Dr. Jin contributes to the control and analysis part.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Jin, Y., Lindsey, M. Stability analysis of genetic regulatory network with additive noises. BMC Genomics 9 (Suppl 1), S21 (2008). https://doi.org/10.1186/147121649S1S21
Published:
DOI: https://doi.org/10.1186/147121649S1S21