Difference between revisions of "Team:NJU-China/RNAi"

Line 269: Line 269:
 
</TD>
 
</TD>
 
     <TD width="73%" bgColor=#FFFFFF style="vertical-align:top;padding-left:80px;padding-right:80px;padding-top:50px;padding-bottom:50px;word-wrap:break-word;">
 
     <TD width="73%" bgColor=#FFFFFF style="vertical-align:top;padding-left:80px;padding-right:80px;padding-top:50px;padding-bottom:50px;word-wrap:break-word;">
−
<h1> 2. RNAi module </h1>
+
<h1> 2. RNAi module </h1>
  
 
<h2> 2.1 Introduction </h2>
 
<h2> 2.1 Introduction </h2>
  
−
&nbsp;&nbsp;&nbsp;
 
  
 
RNA interference (RNAi) is a major tool for transiently suppressing the expression of  
 
RNA interference (RNAi) is a major tool for transiently suppressing the expression of  
Line 290: Line 289:
 
<br><br>
 
<br><br>
  
−
&nbsp;&nbsp;&nbsp;
 
  
 
A challenge for the successful clinical application of RNAi-based drugs is determining  
 
A challenge for the successful clinical application of RNAi-based drugs is determining  
Line 310: Line 308:
 
<h2> 2.2 Model methods </h2>
 
<h2> 2.2 Model methods </h2>
  
−
&nbsp;&nbsp;&nbsp;
 
  
 
This model is inspired by the paper written by Bartlett and Davis [1]. The system uses  
 
This model is inspired by the paper written by Bartlett and Davis [1]. The system uses  
Line 338: Line 335:
 
these processes are included in the model equations.
 
these processes are included in the model equations.
  
−
&nbsp;&nbsp;&nbsp;
+
<br><br>
−
+
 
 
Applying the usual mass action to the reaction network, we can easily obtain the  
 
Applying the usual mass action to the reaction network, we can easily obtain the  
  
Line 345: Line 342:
  
 
<br><br>
 
<br><br>
−
 
−
&nbsp;&nbsp;&nbsp;
 
  
 
The RISC complex, derived from endosomal escape, may disassociate into free RISC and  
 
The RISC complex, derived from endosomal escape, may disassociate into free RISC and  
Line 352: Line 347:
 
siRNA, form an activated mRNA-RISC complex or be degraded, as represented by  
 
siRNA, form an activated mRNA-RISC complex or be degraded, as represented by  
  
−
<I>kdisRISC</I>, <I>kformRISCm and <I>kdegRISC</I>, respectively. The free RISC and  
+
<I>kdisRISC</I>, <I>kformRISCm</I> and <I>kdegRISC</I>, respectively. The free RISC and  
  
 
siRNA may again form a RISC complex, which is represented by <I>kformRISC</I>. The  
 
siRNA may again form a RISC complex, which is represented by <I>kformRISC</I>. The  
Line 381: Line 376:
 
Equation_RNAi_2.jpg"> <br><br>
 
Equation_RNAi_2.jpg"> <br><br>
  
−
&nbsp;&nbsp;&nbsp;
 
  
 
Activated RISC complex bound to mRNA induces the cleavage of target mRNA  
 
Activated RISC complex bound to mRNA induces the cleavage of target mRNA  
Line 387: Line 381:
 
(<I>kcleavage</I>). Additionally, activated RISC complex may undergo degradation  
 
(<I>kcleavage</I>). Additionally, activated RISC complex may undergo degradation  
  
−
(<I>kdegRISC</I>) or disassociation (<I>kdisRISC</I>m).
+
(<I>kdegRISC</I>) or disassociation (<I>kdisRISCm</I>).
  
 
<br><br>
 
<br><br>
Line 395: Line 389:
 
Equation_RNAi_3.jpg"> <br><br>
 
Equation_RNAi_3.jpg"> <br><br>
  
−
&nbsp;&nbsp;&nbsp;
 
  
 
The balance of formation (<I>kformmRNA</I>) and degradation (<I>kdegmRNA</I>) of mRNA  
 
The balance of formation (<I>kformmRNA</I>) and degradation (<I>kdegmRNA</I>) of mRNA  
Line 409: Line 402:
 
Equation_RNAi_4.jpg"> <br><br>
 
Equation_RNAi_4.jpg"> <br><br>
  
−
&nbsp;&nbsp;&nbsp;
 
  
 
The variables and parameters of this model can be accessed here. All the parameters we  
 
The variables and parameters of this model can be accessed here. All the parameters we  
Line 419: Line 411:
 
<h2> 2.3 Results </h2>
 
<h2> 2.3 Results </h2>
  
−
&nbsp;&nbsp;&nbsp;
 
  
 
The simulation results demonstrate the effect of siRNA on MOR knockdown in vivo.  
 
The simulation results demonstrate the effect of siRNA on MOR knockdown in vivo.  
Line 435: Line 426:
 
<br><br>
 
<br><br>
  
−
&nbsp;&nbsp;&nbsp;
 
  
 
As depicted in the figure, the expression level of MOR protein exhibits a rapid  
 
As depicted in the figure, the expression level of MOR protein exhibits a rapid  
Line 457: Line 447:
 
<br><br>
 
<br><br>
  
−
&nbsp;&nbsp;&nbsp;
 
  
 
To explore the dose effect on MOR knockdown, we set different initial conditions and  
 
To explore the dose effect on MOR knockdown, we set different initial conditions and  
Line 471: Line 460:
 
<br><br>
 
<br><br>
  
−
&nbsp;&nbsp;&nbsp;
 
  
 
To optimize the dosing schedule, the lasting times for efficiency and recovery needs to  
 
To optimize the dosing schedule, the lasting times for efficiency and recovery needs to  
Line 523: Line 511:
  
 
effects of heroin. Neuroscience, 158, 474-483. <br>
 
effects of heroin. Neuroscience, 158, 474-483. <br>
 +
 +
 +
  
  

Revision as of 18:10, 17 September 2015

model

2. RNAi module

2.1 Introduction

RNA interference (RNAi) is a major tool for transiently suppressing the expression of genes. Many mathematical models have been constructed to elucidate the mechanism of RNA interference and provide accurate predictions. Nevertheless, most of the current models focus merely on RNAi and fail to consider the delivery process. We modeled the delivery process and the input variant in this module should be the output result of the delivery module.

A challenge for the successful clinical application of RNAi-based drugs is determining the dosing schedule required for efficacy. Patients may ask the following questions: How soon will RNAi-based drugs take to exert efficacy after injection? How long will the efficacy last? What is the dose I need to take, and will it be too costly? How soon will the level of Mu opioid receptor (MOR) protein recover? Is RNAi therapy safe enough? Mathematical modeling using simple kinetic equations for each step in the RNAi process can shed light on many of these questions.

2.2 Model methods

This model is inspired by the paper written by Bartlett and Davis [1]. The system uses the presence of the RISC complex, which is formed in exosomes and escaped from endosome, as a stable source to provide silencing power. Then, the RISC units are targeted to mRNA having the same sequence as the siRNA that triggers this process, binding with mRNA to form an activated RISC-mRNA complex. Once bound to complementary mRNA, activated RISC may induce the degradation of mRNA and further silence protein expression.



Figure 5. Schematic diagram of RNA interference pathway. Degradation of the RISC complex, siRNA, mRNA and protein is not shown here for clear illustration. However, these processes are included in the model equations.

Applying the usual mass action to the reaction network, we can easily obtain the following model equations:

The RISC complex, derived from endosomal escape, may disassociate into free RISC and siRNA, form an activated mRNA-RISC complex or be degraded, as represented by kdisRISC, kformRISCm and kdegRISC, respectively. The free RISC and siRNA may again form a RISC complex, which is represented by kformRISC. The amount of free RISC proteins available for the formation of activated complex is rtot (free RISC protein) – R – C – kdisRISC*R (disassociated RISC protein derived from endosomes). Thus, the total numbers of siRNA- RISC complexes can be modeled using the equations below.



    The number of free siRNA in the cytosol is governed by the equation below.



Activated RISC complex bound to mRNA induces the cleavage of target mRNA (kcleavage). Additionally, activated RISC complex may undergo degradation (kdegRISC) or disassociation (kdisRISCm).



The balance of formation (kformmRNA) and degradation (kdegmRNA) of mRNA and protein (kformprot kdegprot) is interrupted by RISC-induced cleavage of mRNA.



The variables and parameters of this model can be accessed here. All the parameters we used in this module are reported in the literature [1].

2.3 Results

The simulation results demonstrate the effect of siRNA on MOR knockdown in vivo.



Figure 6. Effect of siRNA on MOR mRNA (B) and protein (A) knockdown in vivo. The quantity of total exosomes injected is 300 μg which contains 3 nmol siRNA.

As depicted in the figure, the expression level of MOR protein exhibits a rapid exponential decay and reaches lowest level 12 hours after exosome injection, following a similar pattern observed with the relative level of MOR mRNA.



Figure 7. Effect of dose on MOR mRNA (A) and protein (B) knockdown in vivo. The initial quantity of total exosome injected was set at 50 μg, 100 μg, 200 μg, 400μg and 600 μg, containing 0.5 nmol, 1 nmol, 2 nmol, 4 nmol and 6 nmol siRNA, respectively.

To explore the dose effect on MOR knockdown, we set different initial conditions and ran simulations. The result shows that the concentrations of exosomes and siRNA have significant impact on knockdown efficiency and recovery time. A high dosing schedule leads to more complete knockdown of MOR protein and takes longer for MOR protein levels to recover.

To optimize the dosing schedule, the lasting times for efficiency and recovery needs to be considered. Elongating the lasting time of drug efficiency while shortening the recovery time seems paradoxical based on the simulation data. The literature has reported that 3 nmol of siRNA is adequate for repressing reward effects after 7 days of injection with the relative level of MOR protein reaching approxiamately 80%. If we assume that the threshold of relative level of MOR protein below which opioid reward effects are repressed, is 80% [2], then injecting 400 μg exosome (4 nmol siRNA) might be the best choice. The efficacy of the drug could last for about one week, and another week would be required for MOR protein levels to absolutely recover. Increasing the frequency of dosing may also help to lengthen the drug efficacy time.

2.4 Model Variables

*********************这里插第一幅表格*********************************************

2.5 Model Parameters

*********************这里插第二幅表格*********************************************

References:
1.Bartlett, D.W. and Davis, M.E. (2006) Insights into the kinetics of siRNA-mediated gene silencing from live-cell and live-animal bioluminescent imaging. Nucleic Acids Res, 34, 322-333.
2.Zhang, Y., Landthaler, M., Schlussman, S.D., Yuferov, V., Ho, A., Tuschl, T. and Kreek, M.J. (2009) Mu opioid receptor knockdown in the substantia nigra/ventral tegmental area by synthetic small interfering RNA blocks the rewarding and locomotor effects of heroin. Neuroscience, 158, 474-483.