 Research article
 Open Access
 Published:
Periodicity and dosage optimization of an RNAi model in eukaryotes cells
BMC Bioinformatics volume 20, Article number: 340 (2019)
Abstract
Background
As a highly efficient and specific gene regulation technology, RNAi has broad application fields and good prospects. The effect of RNAi enhances as the dosage of siRNA increases, while an exorbitant siRNA dosage will inhibit the RNAi effect. So it is crucial to formulate a doseeffect model to describe the degradation effects of the target mRNA at different siRNA dosages.
Results
In this work, a simple RNA interference model with hill kinetic function (Giulia Cuccato et al. (2011)) is extended. Firstly, by introducing both the degradation time delay τ_{1} of mRNA caused by siRNA and the transportation time delay τ_{2} of mRNA from the nucleus to the cytoplasm during protein translation, one acquires a novel delay differential equations (DDEs) model with physiology lags. Secondly, qualitative analyses are executed to identify regions of stability of the positive equilibrium and to determine the corresponding parameter scales. Next, the approximate period of the limit cycle at Hopf bifurcation points is computed. Furthermore we analyze the parameter sensitivity of the limit cycle. Finally, we propose an optimal strategy to select siRNA dosage which arouses significant silencing efficiency.
Conclusions
Our researches indicate that when the dosage of siRNA is large, oscillating periods are identical for disparate number of siRNA target sites even if it greatly impacts the critical siRNA dosage which is the switch of oscillating behavior. Furthermore, parametric sensitivity analyses of limit cycle disclose that both of degradation lag and maximum degradation rate of mRNA due to RNAi are principal elements on determining periodic oscillation. Our explorations will provide evidence for gene regulation and RNAi.
Background
The mechanism for sequencespecific posttranscriptional gene silencing that is induced by doublestranded RNA (dsRNA), leading to the regression of the target messenger RNA (mRNA) [1]. This common phenomenon in many eukaryotes, including insects, is named RNA interference (RNAi) by Fire et al. RNAi in animals [2] and in plants [3], is an evolutionarily conservative defense against transgenic or exotic virus infringement mechanism [4]. The process of RNAi can be divided into four stages:

Step 1. Double stranded RNA (dsRNA) expressed in or introduced into the cell is cleaved into fragments of 2123 base pairs (called small interfering RNA, abbreviated as siRNA) by the Dicer enzyme.

Step 2. siRNAs are firstly adhered to RNA Induced Silencing Complex (RISC), and whereafter split into the sense strands which are deserted [5], and the antisense strands which are still roped to RISC.

Step 3. An available siRNARISC complex, includes the siRNA loaded to the Ago protein, is packaged by antisense strand. Then it identifies and unites target mRNAs via the principle of complementary base pairing.

Step 4. The antisense strand commands a endonuclease bound to RISC (an Argonaute protein called ‘slicer’) to operate the degradation of the target mRNA. And next the complex is liberated to dispose further mRNA targets.
In recent years, some studies have shown that synthetic siRNAs can effectively trigger RNAi in eukaryotes [6] and the siRNAs seem to avoid offtarget effects prompted by longer doublestranded RNAs in mammalian cells [7]. This discovery makes the application of RNAi technology more convenient. The high efficiency and specificity of RNAi make it become a powerful tool for researching gene function. RNAi also provides a novel idea to schedule synthetic biological circuits for synthetic biology [8]. In the treatment of certain genetic diseases, for example viral infections [9], cancer [10] and inherited genetic disorders [11], RNAi has the potential to become a new type of therapeutic tool. In the field of pest management, RNAi also shows its talents [12]. And RNAi technology first approved by the US Environmental Protection Agency as a pesticide in 2017.
Because excessive siRNAs not only affect its efficiency [13], but attract offtarget effect [7]. For RNAi application, it is necessary to find a quantitative mathematical model that can describe the relationship between the dosage of siRNA and the RNAi effect. Giulia Cuccato et al. (2011), according to vitro experimental data and squared error measure, capture the most efficient mathematical model of RNA interference in [13].
However, for the model, we consider that there are two important time delays that cannot be ignored during the entire RNAi process. First, degradation of mRNA due to RNAi. Here, we use τ_{1} to describe this time delay. Next, carriage of mRNA from nucleus to cytoplasm. Thus, we introduce τ_{2} to represent this time delay. In our work, we start from the model proposed in [13] and then modify it. First, we conduct a qualitative analysis of the model with delay. Our result show that the stability of the only positive equilibrium has changed: it is stable while the original model without time delays, as the time delay increases, it will turn into damped oscillation and lose its stability via a Hopf Bifurcation. Therefore, time delay plays an important role in dynamics of RNAi model and should not be ignored in the modeling of genetic regulation. Next, we introduce the solution to the periodic value of the periodic solution of the system with the limit cycle. And we analyze the parameter sensitivity of the amplitude and period of a periodic solution for a system with a limit cycle. Finally, we give optimal control for quantitative RNAi model by optimization theory.
Results
Qualitative analysis
When the delays are finite, the characteristic equations are functions of delays. As values of the delays change, the stability of the trivial solution may also changes. Such phenomena is often refereed to as stability switches. Next the qualitative analysis of model (20) will be conducted.
Stability and Hopf Bifurcation
In this section, we discuss the local asymptotic stability of the unique positive equilibrium \(Q^{*}(\tilde M, \tilde P)\) and the existence of Hopf bifurcation. Setting β=rS^{n}/(θ^{n}+S^{n}), \(\tilde M \) and \(\tilde P\) are denoted by
For τ_{1,2}>0, characteristic equation of model (20) is given by
Obviously, λ_{1}=−d_{p} is a negative root of the Eq. (1). Next let the first item of the left side of the Eq. (1) be
Lemma 1
For ω∈[π/(2τ_{1}),π/τ_{1}], let \(\phantom {\dot {i}\!}\beta _{0}=e^{d_{m}\tau _{1}1}/\tau _{1},\) β_{1}=−d_{m}/ cos(ωτ_{1})>0. Then the following results hold.
(a) If β<β_{0}, f(λ) has two real negative roots.
(b) If β=β_{0}, f(λ) has one real negative root.
(c) If β_{0}<β<β_{1}, f(λ) has two complex conjugate roots with Re(λ)<0.
(d) If β=β_{1}, f(λ) has two complex conjugate roots with Re(λ)=0.
(e) If β>β_{1}, f(λ) has two complex conjugate roots with Re(λ)>0.
Proof
Function (2) implies f(−∞)=+∞, f(+∞)=+∞ and f(0)=d_{m}+β>0. Then, letting \(\phantom {\dot {i}\!}f^{\prime }(\lambda)=1\beta \tau _{1}e^{\lambda \tau _{1}}=0\) yields
So, f(λ) maybe has negative root only if βτ_{1}<1. In addition, because f(λ^{∗}) is minimum of f(λ) for every \(\lambda \in \mathbb {R}\), thus the function (2) has one real negative root λ^{∗} if f(λ^{∗})=0, namely, β=β_{0}. Hence (b) is proved. If β<β_{0}, we obtain f(λ^{∗})<0 and the function (2) has two negative real roots, then (a) is proved too.
For β>β_{0}, the Eq. (2) may has two complex roots. For that, we assume that there exists a solution of the characteristic equation of the form λ=iω(ω>0). Putting it into f(λ), it follows
Comparing real and imaginary parts we get,
Squaring and adding the first and the second of (3), we get \((d_{m}^{2}+\omega ^{2})/\beta ^{2}=1\), that is, \(\omega ^{2}=\beta ^{2}d_{m}^{2}\). Hence positive solution \(\omega _{0}=\sqrt {\beta ^{2}d_{m}^{2}}\) exists if β>d_{m}. And, corresponding to λ=iω_{0} and the first equation of (3), there exists \(\tau _{1}^{*}>0\) such that,
and (c) and (d) are proved. When β>β_{1}, Re(λ)>0, so (e) is proved. □
Theorem 1
For model (20), the following results hold.
(a) If β≤β_{0}, then the equilibrium Q^{∗} is asymptotically stable.
(b) If β_{0}<β<β_{1}, then the equilibrium Q^{∗} is oscillatory stable.
(c) If β>β_{1}, then the equilibrium Q^{∗} is unstable. Furthermore if β=β_{1}, Hopf bifurcation occurs.
Proof
(a) and (b) are apparently valid by Lemma 1. Now differentiating (1) with respect to τ_{1} gives
Then at λ=iω_{0}, one gets
Thus
□
So, Re(λ)>0 providing \(\tau _{1}>\tau _{1}^{*}\). By the Hopf bifurcation theorem, the condition with \(\tau _{1}>\tau _{1}^{*}\) and Re(λ)>0 guarantees that the Hopf bifurcation at β=β_{1} is supercritical. The result (c) is proved.
The bifurcation diagram of Eq. (20) as a function of the delay τ_{1} and of the parameter β is shown in Fig. 1a. Using some parameter values suggested in [14], we take k_{m}=10, d_{m}=0.05, k_{p}=1, d_{p}=0.01, other parameter values are set by θ=10, n=4, S=30. It is always possible to choose values of r, θ, n and S such that β>β_{1}(τ_{1}), where β_{1}(τ_{1}) is the parameter that determines a supercritical Hopf Bifurcation. In this case, the delay model (20) has asymptotically stable oscillatory solutions (limit cycle solutions) in Fig. 1a, and the time evolution of protein is shown in Fig. 1c. The phase diagrams for system (20) in damped oscillation region and at the limit cycle are shown in Fig. 1b and d, respectively.
Period of the bifurcating oscillatory solution
We knew in the above subsection how the delay differential equation model was capable of generating limit cycle periodic solutions. One indication of their existence is if the steady state is unstable by growing oscillations, although this is certainly not conclusive. From the analysis of the previous section, we found that τ_{2} has no effect on the stability of the system, and the first equation of (20) is independent. So, resorting to the periodicity of the first equation, we intent to analyze the period of the whole system as the delay occurs. Hence pull out the first equation of (20) separately,
Linearising (4) at the first component of the steady state, \(\tilde {M}=k_{m}/(d_{m}+\beta)\), that is, writing \(M(t)\tilde {M}=m(t)\), yields
By looking for solutions m(t) in the form m(t)=ce^{λt}, we get
where c is a constant and the eigenvalues λ are solutions of (6), a transcendental equation in which τ_{1}>0. It is not easy to find the analytical solutions of (6). However, all we really want to know from a stability point of view is whether there are any solutions with Re(λ)>0 which from the form of m(t) implies instability since in this case m(t) grows exponentially with time.
Putting λ=μ+iω, in (6), and now take the real and imaginary parts of the transcendental equation in (6), namely,
We knew that the steady state Q^{∗} is stable if \(0<\tau _{1}<\tau _{1}^{*}\) and the delay Eq. (20) has an stable periodic solution for \(\tau _{1}=\tau _{1}^{*}\). In the latter case we expect the solution to exhibit stable limit cycle behaviour. The critical value \(\tau _{1}=\tau _{1}^{*}\) is the bifurcation value. The effect of delay in models is usually to increase the potential for instability. Here as τ_{1} is increased beyond the bifurcation value \(\tau _{1}^{*}\), the steady state becomes unstable.
Near the bifurcation value we can get a estimate of the period of the bifurcating oscillatory solution as follows. Consider the dimensionless form and let
The solution λ=μ+iω, of (7), with the Re(λ)=0 when \(\tau _{1}=\tau _{1}^{*}\) is μ=0, \(\omega =\omega _{0}=\sqrt {\beta ^{2}d_{m}^{2}}\). For ε small we expect μ and ω to differ from μ=0 and ω=ω_{0} also by small quantities so let
where δ and σ are to be determined. Substituting these into the second of (7) and expanding for small δ, σ and ε gives
while the first of (7) gives
Thus on solving these simultaneously
Considering that the imaginary part of λ is the root cause of m(t) periodicity, near the bifurcation, the first of (5) with (6) gives
This shows that the delay Eq. (20) has an stable periodic solution due to the occurrence of the Hopf bifurcation with period
T describes the relations of the Hopf bifurcation period with the maximal regression rate r, siRNA dosage S, the half saturation coefficient θ and the number of siRNA target sites n. Here, Fig. 2 reveals that the levels of mRNA and the protein are persistence when the dosage of siRNA are small, otherwise the periodic oscillating happens. Meanwhile, it indicates that when the dosage of siRNA is large, oscillating periods are identical for disparate number of siRNA target sites even though it greatly impacts the critical siRNA dosage S_{n} which is the switch of oscillating behavior.
Our delayed differential equations are applied to model gene regulatory network due to RNAi. The periodic solutions of delayed differential equations are subjected to parameters. So, it is necessary that a parametric sensitivity analysis for amplitude and period of periodic solutions.
Parametric Sensitivity
In this section, we use a sensitivity analysis method proposed in [15], and focus on sensitivities of amplitude and of the period when our delay model (20) possesses a periodic solution. Then define the sensitivity equation for parameter d_{m}:
In similar way, we get the sensitivity equation with respect to τ_{1}:
Analogously, the other sensitivity equations on the rest parameters can be captured, I won’t list them here. Solving the there equations, and according to the circumscription of sensitivities of the limit cycle in [15], we obtain the relative sensitivities of the amplitude and of the period shown in Fig. 3.
We observe that τ_{1}, RNAi process delay, has a effective impact on both amplitude and period, while τ_{2}, the mRNA translational delay, has inappreciable influence. Because, the occurrence of the limit cycle is only related to the value of τ_{1}, and τ_{2} does not affect the stability of the equilibrium point of model (20). Moreover, parameter r, the maximal degradation rate of the mRNA due to RNAi, has a important affection on period too. This is because the value of r determines the satisfaction of β=β_{1}. When β=β_{1}, the system (20) will have a limit cycle, where β is the degradation rate of mRNA due to RNAi. In other words, in eukaryotic cells, if the rate of degradation of mRNA due to RNAi is greater than the rate of degradation of the mRNA itself, τ_{1} and r will be important parameters in the quantitative delay system (20).
Optimizing the dosage of siRNA in RNAi
During the RNAi, excessive siRNAs not only affect the efficiency, but also attracts offtarget effect. So the rational dosage of siRNA is crucial for both enhancing RNAi efficiency and reducing cost.
Optimal control for model without delay
Define a cost function as
where P_{s} is the cost of per unit siRNA, T is the terminal time. The first part of (10) represents the cost of siRNA consumed in [0,T], and the second part shows the accumulation of protein (denoted by PA) in [0,T]. The aim of this work is to minimize the cost of siRNA and the accumulation of protein.
Problem(Q1). For model (19), choose S∈[0,200] (according to the experiments in [13]) to minimize the cost function (10).
Since the constraint of the cost function (10) is only the state equation, it must be observed during [0,T]. Then the Lagrange multiplier vector can be used to introduce the equality constraint into the integrand part of the definite integral, thus transforming the constrained optimization problem into an unconstrained optimization problem. Then there is
Introduce a Hamiltonian function H as follows:
Then corresponding costate equations is determined by
and the corresponding gradients of the cost function (11) with respect to S is
Optimal control for model with delay
Problem (Q2) For model (20), opt for S∈[0,200] (according to the experiments in [13]) to minimize the cost function (10).
Rewriting the cost function (10) in the same way, yields
Define a Hamiltonian function H by
The corresponding costate equations is dominated by
with jump conditions
and boundary conditions
The corresponding gradients of the cost function (14) with respect to S is governed by
According to [14], we take k_{m}=10, d_{m}=0.05, k_{p}=1, d_{p}=0.01, and set other parameter values r=0.02, θ=10, n=4, P_{s}=1, T=60, M(0)=160, P(0)=10000. Then, we solve two optimal problems with these parameter values by using Matlab programs.
Simulation 1. Comparison of the optimal value and not optimal value about model without delay
The solution obtained by the optimizer is S=43.01. Substituting S into β, one gets β=0.0199≈r=0.02. This shows that the degradation rate of mRNA due to RNAi has reached the maximum. When S>43.01, correspondingly, β is almost unanimously close to r. It implies that the amounts of protein accumulated are identical and the degradation of mRNA is subject to saturation effects when siRNA dosage is larger. Meanwhile, by employing the optimal result S=43.01 together with a nooptimal value S=10 of siRNA dosage and remaining parameters are as above, we make a comparison about the time evolution of mRNA and protein dosages of model without delay under different siRNA dosages controls (see Fig. 4).
Simulation 2. Comparison of the optimal value and not optimal value about model with delay
For this case, take τ_{1}=2.5, τ_{2}=1 and the other parameters are taken as those in simulation 1. The solution obtained by the optimizer is S=30.37. After substituting, one gets β=0.0198≈r=0.02. Although the two results differ by 12.64, the corresponding β is almost the same. This shows that the degradation rate of mRNA due to RNAi is almost the highest under the optimal conditions. Similarly, when S>30.37, the RNAimediated degradation of mRNA is subject to saturation effects, we also make a comparison like simulation 1 at S=30.37 and S=10 (see Fig. 5). In addition, Table 1 gives the value of the optimal siRNA dosage, protein accumulation (PA) and the cost function value J. Obviously, with the participation of time delays, the accumulation of protein is much lower than when there is no time delays, although both of S are taken at the best value.
Discussion
What we interest in is a mathematical model that reflects the relationship between the RNAi effect and the siRNA dose, which is called the doseeffect model. The study had three primary goals. The first was to depict and forecast the evolution rules of mRNA and protein by the dynamic analysis. The second was to study the effect of parameters on periodic oscillation. The third was to explore the optimal dosage for the significant silencing efficiency. Our work provides a theoretical basis for more precise and economical RNAi experiments and applications. Even so, there are some questions worth exploring further. One is that the degradation and amplification process of siRNA should be considered in RNAi model. The second is that the stochastic effects and variable siRNA dosage should be involved in our model. These factors will result in more complicated dynamic behaviors and reveal more mechanisms of RNAi.
Conclusions
In this paper, we reference a simple Hill kinetic model proposed by [13] and consider the potential effect of two time delays. One is degradation of mRNA due to RNAi, other one is carriage of mRNA from nucleus to cytoplasm. For the improved timedelay system, the role of time delays and the dynamic behavior of system are discussed. Qualitative analyses indicate that the introduction of time delays changes the dynamic behaviors of the system. In detail, as delays increase, the unique positive equilibrium firstly is oscillatory stable and then loses its stability via a Hopf Bifurcation. Furthermore, we give the corresponding parameter scales for these results. Meanwhile, the period of the oscillation solution shows that when the dosage of siRNA is large, oscillating periods are identical for disparate number of siRNA target sites in spite of it greatly impacts the critical siRNA dosage which is the switch of oscillating behavior. And then, parametric sensitivities of the limit cycle is determined. The results indicate that both of degradation lag and maximum degradation rate of mRNA due to RNAi are principal elements on determining periodic oscillation. After that, we propose and solve a simple optimization problem for ODEs model (19) and DDEs model (20) based on the optimization theory. The rational dosage of siRNA is given for both enhancing RNAi efficiency and reducing cost by a Matlab program. The results imply that the optimal dosage of siRNA with delay effects is less than one without time delay.
Methods
In this section, we apply and expand the model recommended in [13]. This model well describes the mRNA and protein level in RNA interference process for different dosages of siRNA in mammalian cells in vitro, and great predicts the saturation effect observed experimentally of the RNAi process [13]. The RNAi process caused by siRNA (S) is encapsulated into a whole, and the degradation of the target mRNA (M) due to RNAi is expressed in the form of a functional reaction. In addition, the protein corresponding to the target mRNA is denoted as P. The time evolution of the dosages of mRNA and protein can be described by the ordinary differential equations (ODEs) as follows:
where M is transcribed at a rate k_{m} from the promoter; d_{m} and d_{p} are the degradation rates of M and P, respectively. P is translated at a rate k_{p} form M. The extra degradation rate of M as a result of RNAi is the third segment of the first equation of (19), which is a Hillkinetic model. Positive integer n is a Hill coefficient, representing the number of siRNA bounded on the target mRNA (or the number of siRNA target sites). r and θ tie to the potency of RNAi induced by siRNA [16]: r denotes the maximal regression rate of M because of RNAi, θ is the dosage of S required to reach half of the maximal degeneration rate r.
Time delay plays an important role in many biological dynamical systems. There are two important biological delays that must be considered when modeling RNAi. One is the RNAi process caused by siRNA, using τ_{1} to describe it. The other one is the transportation process of mRNA from nucleus to cytoplasm, introducing τ_{2} to represent it. Then, the time evolution of the dosages of mRNA and protein can be described by the following delay differential equations (DDEs):
with the initial condition: M(t)=M(0) and P(t)=P(0) for −max{τ_{1},τ_{2}}≤t≤0. It is assumed that all the parameters of model (20) are positive.
In real RNAi experiments and applications, the biological time delays are ubiquitous, such as inhibiting the expression of chitinase of migratory locust, gene knockout in animal and inhibiting cancer proliferation. Therefore, our improved time delay model is more convincing in describing the relationship between siRNA measurement and RNAi efficiency in eukaryotic cells.
Availability of data and materials
Data sharing is not applicable to this article as no datasets were generated or analysed during the current study.
Abbreviations
 DDEs:

Delay differential equations
 dsRNA:

Doublestranded RNA
 mRNA:

Messenger RNA
 ODEs:

Ordinary differential equations
 RISC:

RNA induced silencing complex
 RNAi:

RNA interference
 siRNA:

Small interfering RNA
References
 1
Fire A, Xu S, Montgomery MK, Kostas SA, Driver SE, Mello CC. Potent and specific genetic interference by doublestranded rna in caenorhabditis elegans. Nature. 1998; 391(6669):806–811.
 2
Hannon GJ. Rna interference. Nature. 2002; 418(6894):244–251.
 3
Baulcombe DC. Rna silencing in plants. Nature. 2004; 431(7006):356–363.
 4
Chen YMZ. Rna interference. Chinese Journal of Biological Engineering. 2003; 23(3):39–43. in Chinese.
 5
Filipowicz W. Rnai: the nuts and bolts of the risc machine. Cell. 2005; 122(1):17–20.
 6
Brummelkamp TR, Bernards R, Agami R. A system for stable expression of short interfering rnas in mammalian cells. Science. 2002; 296(5567):550–553.
 7
Caplen NJ, Parrish S, Imani F, Fire A, Morgan RA. Specific inhibition of gene expression by small doublestranded rnas in invertebrate and vertebrate systems. Proceedings of the National Academy of Sciences of the United States of America. 2001; 98(17):9742–9747.
 8
Deans TL, Cantor CR, Collins JJ. A tunable genetic switch based on rnai and repressor proteins for regulating gene expression in mammalian cells. Cell. 2007; 130(2):363–372.
 9
Barik S, Bitko V. Prospects of rna interference therapy in respiratory viral diseases: update 2006. Expert Opinion on Biological Therapy. 2006; 6(11):1151–1160.
 10
Takeshita F, Ochiya T. Therapeutic potential of rna interference against cancer. Cancer Science. 2006; 97(8):689–696.
 11
Aagaard L, Rossi JJ. Rnai therapeutics: Principles, prospects and challenges. Advanced Drug Delivery Reviews. 2007; 59(2):75–86.
 12
Price DRG, Gatehouse JA. Rnaimediated crop protection against insects. Trends in Biotechnology. 2008; 26(7):393–400.
 13
Cuccato G, Polynikis A, Siciliano V, Graziano M, Bernardo MD, Bernardo DD. Modeling rna interference in mammalian cells. BMC Systems Biology. 2011; 5(1):19–19.
 14
Zhou P, Cai S, Liu Z, Wang R. Mechanisms generating bistability and oscillations in micrornamediated motifs. Physical Review E. 2012; 85(4):041916.
 15
Ingalls B, Mincheva M, Roussel MR. Parametric sensitivity analysis of oscillatory delay systems with an application to gene regulation. Bulletin of Mathematical Biology. 2017; 79(7):1539–1563.
 16
Khanin R, Vinciotti V. Computational modeling of posttranscriptional gene regulation by micrornas. Journal of Computational Biology. 2008; 15(3):305–316.
Acknowledgements
The authors thank the referees for their careful reading of the original manuscript and many valuable comments and suggestions, which greatly improved the presentation of this paper.
Funding
Funding bodies did not play any role in the design of the study and in writing this manuscript.
Author information
Affiliations
Contributions
YP presented the ideas and designed the frame of this paper; TM and MZ finished the proofs, computes and writing of the first draft; CL polished, revised the last draft. All authors read and approved the final manuscript.
Corresponding author
Correspondence to Yongzhen Pei.
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver(http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 RNA interference
 Delay
 Oscillation period
 Sensitivity analyse
 Optimal control