Team:Freiburg Bioware/Modeling/Virus Infection

From 2010.igem.org

(Difference between revisions)
 
(38 intermediate revisions not shown)
Line 1: Line 1:
{{:Team:Freiburg_Bioware/Head}}
{{:Team:Freiburg_Bioware/Head}}
 +
{{:Team:Freiburg_Bioware/jquery}}
{{:Team:Freiburg_Bioware/menu_home}}
{{:Team:Freiburg_Bioware/menu_home}}
<html>
<html>
-
 
<h1>Model for Virus Infection</h1>
<h1>Model for Virus Infection</h1>
-
As in the previous model for the virus production we established a ODE model based on the law of mass action. The following paragraph explains the reaction scheme and our model assumptions. In the subsequent paragraphs the system of differential equations is specified and the implementation in MathWorks® MATLAB is discussed.<br>
+
<p style="text-align: justify;">
-
The last section deals with our modeling results.
+
As described in the previous model for the virus production we established a
 +
model of <span style="font-weight: bold;">ordinary differential
 +
equation</span> (<span style="font-style: italic;">ODE</span>) based on
 +
the law of mass action. The following paragraph depicts
 +
the reaction scheme and our model assumptions. In the subsequent
 +
paragraphs the system of differential equations is specified and the
 +
implementation in MathWorks® MATLAB is discussed. Furthermore
 +
reasonable model extensions are presented.<br>
 +
The last section deals with our modeling results and gives an outlook
 +
to further modeling applications.
<br>
<br>
<br>
<br>
-
 
+
</p>
<h2>Reaction Scheme</h2>
<h2>Reaction Scheme</h2>
-
To simplify the mathematical description of the reaction scheme we divide the cell into four compartments: the <b>extracellular matrix</b> (all quantities with the index <i>ext</i>), the space in <b>endosomes</b> (<i>end</i>), the <b>cytoplasm</b> (<i>cyt</i>) and the <b>nucleus</b> (<i>nuc</i>).<br>  
+
<p style="text-align: justify;">
-
A target cell is transduced by <b>viral particles</b> (<i>V</i>) in the extracellular matrix. Depending on their <b>degree of modification</b> (<i>m</i>) and thus their specificity they can bind to <b>receptors</b> (<i>R</i>) on the cell surface. Once a receptor has formed a <b>complex with a virus particle</b> (<i>VR</i>) receptor dimerization (<i>R<sub>2</sub></i>) occurs and the whole <b>complex</b> (<i>VR<sub>2</sub></i>) is invaginated and endosomes are formed. The virus is released from the endosome to the cytoplasm and is transported to the nucleus where uncoating of the capsid is initiated and the <b>single stranded DNA</b> (<i>ssDNA</i>) is released.<br>
+
To simplify the mathematical description of the reaction scheme we
-
Finally viral mRNA is processed and transported into the cytoplasm and the <b>enzyme for therapeutic approach</b> (<i>E</i>) is produced.  
+
divide the cell into four compartments: the <b>extracellular matrix</b>
-
<br><br>
+
(all quantities with the index <i>ext</i>), the space in <b>endosomes</b>
-
 
+
(<i>end</i>), the <b>cytoplasm</b> (<i>cyt</i>) and the <b>nucleus</b>
-
 
+
(<i>nuc</i>).<br>
-
<table>
+
A target cell is transduced by <b>viral particles</b> (<i>V</i>) in
-
<td>
+
the extracellular matrix depending on their <b>degree of modification</b>
-
<center>
+
(<i>m</i>) and thus their specificity they can bind to <b>receptors</b>
-
<img width="380" src="https://static.igem.org/mediawiki/2010/e/e1/Freiburg10_VirusInfectionScheme01.png" alt="Reaction scheme for the virus production" />
+
(<i>R</i>) on the cell surface. Once a receptor has formed a <b>complex
 +
with a virus particle</b> (<i>VR</i>) receptor dimerization (<i>R<sub>2</sub></i>)
 +
occurs and the whole <b>complex</b> (<i>VR<sub>2</sub></i>) is
 +
internalized and
 +
endosomes are formed. The virus is released from the
 +
endosome to the cytoplasm and is transported to the nucleus where
 +
uncoating of the capsid is initiated and <b>single stranded DNA</b>
 +
(<i>ssDNA</i>) is released.<br>
 +
Finally viral mRNA is processed and transported into the cytoplasm and
 +
the <b>enzyme for therapeutic approach</b> (<i>E</i>) is produced. <br>
 +
<br>
 +
</p>
 +
<table style="width: 905px;">
 +
<tbody>
 +
<tr>
 +
<td style="width: 60%;">
 +
<center style="width: 380px;"><img
 +
src="https://static.igem.org/mediawiki/2010/e/e1/Freiburg10_VirusInfectionScheme01.png"
 +
alt="Reaction scheme for the virus production" height="800"><br>
 +
<br>
 +
<div style="text-align: left;"><span style="font-weight: bold;">Figure
 +
1:</span> Schematical overview of virus infection<br>
 +
</div>
</center>
</center>
</td>
</td>
<td>
<td>
-
<center>
+
<center style="width: 40%;"><img
-
<img width="359" height="335" src="https://static.igem.org/mediawiki/2010/9/90/Freiburg10_VirusInfection01.png" alt="Reaction scheme for the virus production" />
+
src="https://static.igem.org/mediawiki/2010/9/90/Freiburg10_VirusInfection01.png"
-
</center>
+
alt="Reaction scheme for the virus production" width="399"></center>
</td>
</td>
 +
</tr>
 +
</tbody>
</table>
</table>
-
<br><br>
+
<br>
-
 
+
<br>
<h2>Reduced Reaction Scheme</h2>
<h2>Reduced Reaction Scheme</h2>
-
 
+
<p style="text-align: justify;">
-
To reduce the number of equations and parameters we considered two different approaches for the reaction dynamics.
+
This modeling approach neglects the fastest linear transport processes
-
 
+
and
-
<h3>Model 1 (LMA)</h3>
+
describes the receptor binding and dimerization in terms of the <b>law
-
 
+
of mass action</b> (<i>LMA</i>).
-
The first approach only neglects linear transport processes and describes the receptor binding and dimerization in terms of the <b>law of mass action</b> (<i>LMA</i>).
+
<br>
-
 
+
<br>
 +
</p>
<center>
<center>
-
<img width="359" height="227" src="https://static.igem.org/mediawiki/2010/8/83/Freiburg10_VirusInfection02.png" alt="Reaction scheme for the virus production" />
+
<img
-
</center>
+
src="https://static.igem.org/mediawiki/2010/8/83/Freiburg10_VirusInfection02.png"
-
<br><br>
+
alt="Reaction scheme for the virus production" height="227" width="359"></center>
-
 
+
<br>
-
<h3>Model 2 (MM)</h3>
+
<br>
-
The second approach simplifies the scheme to 3 reactions assuming <b>Michaelis-Menten kinetics</b> (<i>MM</i>) in the quasi steady state approximation (Briggs and Haldane 1925) without explicit receptor dimerization.
+
-
 
+
-
<center>
+
-
<img width="429" height="134" src="https://static.igem.org/mediawiki/2010/3/3f/Freiburg10_VirusInfection03.png" alt="Reaction scheme for the virus production" />
+
-
</center>
+
-
<br><br>
+
-
 
+
<h2>Differential Equations</h2>
<h2>Differential Equations</h2>
 +
<p style="text-align: justify;">The ODE model consists of 7 equations
 +
containing 9 rate constants. To
 +
the equation for the extracellular virus concentration a degradation
 +
term with rate constant <span style="font-style: italic;">(k</span><sub
 +
style="font-style: italic;">7,1</sub><span style="font-style: italic;">)</span>
 +
is added corresponding to the immune response of the target
 +
system. The temporal behavior of <i>ssDNA</i> is supplemented by a
 +
linear
 +
degradation term with rate constant <span style="font-style: italic;">(k</span><sub
 +
style="font-style: italic;">6,1</sub><span style="font-style: italic;">)</span>.<br>
 +
</p>
 +
<center><img
 +
src="https://static.igem.org/mediawiki/2010/8/85/Freiburg10_VirusInfection04.png"
 +
alt="Reaction scheme for the virus production" width="633"></center>
 +
<br>
 +
<br>
 +
<h3>Model Extensions</h3>
 +
<p style="text-align: justify;">Since internalization depends on the <b>degree
 +
of modification</b> (<i>m</i>) the model is extended by a non-linear
 +
dependency in the rate constants. Additionally a decreasing production
 +
efficiency reduces the amount
 +
of available <i>ssDNA</i> in the nucleus. The estimated functional
 +
dependencies are shown in the figure below.<br>
 +
</p>
 +
<table style="text-align: left; width: 905px; height: 481px;" border="0"
 +
cellpadding="2" cellspacing="2">
 +
<tbody>
 +
<tr>
 +
<td style="vertical-align: top; width: 538px;"><img
 +
src="https://static.igem.org/mediawiki/2010/5/54/Freiburg10_ModificationEfficiency.png"
 +
alt="" width="500"><br>
 +
<span style="font-weight: bold;">Figure 2:</span> Assumed
 +
functions for internalization (<span style="font-style: italic;">k(m)</span>,
 +
blue) and fraction of infectious virus particles (<span
 +
style="font-style: italic;">r(m)</span>, red) dependant on
 +
modification degree (<span style="font-style: italic;">m</span>).<br>
 +
</td>
 +
<td
 +
style="vertical-align: top; width: 507px; text-align: justify;">Regarding the <i>k(m)</i>
 +
we assumed that a minimal modification degree of 10%
 +
is necessary for virus receptor complex formation due to the fact that
 +
at least two binding sites are required for receptor dimerization. For
 +
higher <i>m</i>-values the targeting efficiency saturates because
 +
enough binding sites exists that every dimerized receptor can build a
 +
complex with a virus.<br>
 +
<br>
 +
The values of <i>r(m)</i> are normalized to the ratio of infectious
 +
particles to empty virus capsids and therefore start at one.
 +
Increasing <i>m</i> results in a decreasing amount of infectious
 +
particles until a minimal value is reached for the completely modified
 +
virus (m=100%).<br>
 +
<br>
 +
Future experiments could provide a better understanding of those
 +
functional dependencies and therefore improve the model predictions
 +
concerning the targeting efficiency discussed later in this chapter.<br>
 +
</td>
 +
</tr>
 +
</tbody>
 +
</table>
 +
<br>
 +
<br>
 +
<h2>Methods and Simulation</h2>
 +
<p style="text-align: justify;">
 +
The ODE model was implemented in MathWorks® MATLAB R2010b. Integration
 +
of the differential equations was achieved using the stiff integration
 +
routine <i>ode15s</i> with automatic integration step size management.<br>
 +
For initial conditions we took virus concentrations used in the
 +
experiments. Rate constants were estimated according to published
 +
values. The entire internalization process of viral particles being
 +
transported from the extracellular matrix through the cytoplasm into the
 +
nuleus happens within seconds (Seisenberger et al. 2001; Bartlett et
 +
al. 2000) whereas protein synthetization, capsid formation and virus
 +
packaging is taking much more time.<br>
 +
Furthermore experiments have been conducted in order to
 +
determine typical time scales corresponding to biological dynamics.<br>
 +
<br>
 +
The used parameters are given in the table below.
 +
<br>
 +
</p>
<center>
<center>
-
<img width="655" height="314" src="https://static.igem.org/mediawiki/2010/8/85/Freiburg10_VirusInfection04.png" alt="Reaction scheme for the virus production" />
+
<table style="text-align: left; width: 432px; height: 239px;" border="0"
 +
cellpadding="2" cellspacing="2">
 +
<tbody>
 +
<tr>
 +
<td style="vertical-align: top;">
 +
<div style="text-align: center;"><img
 +
src="https://static.igem.org/mediawiki/2010/c/c7/Freiburg10_RateConstants02.png"
 +
alt="" width="405"><br>
 +
</div>
 +
<br>
 +
<span style="font-weight: bold;">Table 1:</span> rate constants
 +
for the virus infection model<br>
 +
</td>
 +
</tr>
 +
</tbody>
 +
</table>
</center>
</center>
-
<br><br>
 
-
<center>
 
-
<img width="800" height="285" src="https://static.igem.org/mediawiki/2010/8/84/Freiburg10_VirusInfection05.png" alt="Reaction scheme for the virus production" />
 
-
</center>
 
-
 
-
 
-
<br><br>
 
-
<h2>Methods and Simulation</h2>
 
-
The ODE model was implemented in MathWorks® MATLAB R2010b. Integration of the differential equations was achieved using the stiff integrator <i>ode15s</i> with automatic integration step size management.<br>
 
-
To calibrate the dynamics of the mathematical model to those of biological system we used time lapse data of fluorescence experiments as well as published values for the rate constants.
 
-
<br><br>
 
-
The parameters used are given in the table below. Also you can download the MATLAB source code.
 
<br>
<br>
-
 
+
<br>
-
<h3>Model 1 (LMA)</h3>
+
<a
-
<img width="405" height="337" src="https://static.igem.org/mediawiki/2010/c/c7/Freiburg10_RateConstants02.png" alt="" />
+
href="https://static.igem.org/mediawiki/2010/f/f8/Freiburg10_VirusInfectionCode.m">get
-
<br><br>
+
the m-File (MATLAB source code)!</a><br>
-
 
+
<br>
-
<h3>Model 2 (MM)</h3>
+
<br>
-
<img width="405" height="337" src="https://static.igem.org/mediawiki/2010/0/0e/Freiburg10_RateConstants03.png" alt="" />
+
-
<br><br>
+
-
 
+
-
<a href="https://2010.igem.org/Wiki/images/8/87/Freiburg10_VirusInfectionCode.m">get the .m-File (MATLAB source code)</a><br>  
+
-
 
+
-
 
+
-
<br><br>
+
<h2>Results and Discussion</h2>
<h2>Results and Discussion</h2>
-
 
+
<p style="text-align: justify;">
-
<h3>Model 1 (LMA)</h3>
+
Figure 3 shows the time course of the virus infection model. Starting
 +
with an initial virus concentration of 10<sup>7</sup> viral particles
 +
per mililiter
 +
the curve decays fast as the infectious particles are bound to
 +
receptors and imported to endosomes in the cytoplasm. Since
 +
recepor-virus complexes are built the amout of free recepors on the
 +
cell surface decreases while the complex concentrations are rising. <br>
 +
As soon as the first particles reach the nucleus and release their
 +
viral DNA the intranuclear concentration of <span
 +
style="font-style: italic;">ssDNA</span> rises and enzyme
 +
transcription is initialized. Therefore the enzyme concentration in
 +
cytoplasm is increased until a steady state is reached due to the fact
 +
that <span style="font-style: italic;">ssDNA</span> is degraded inside
 +
the nucleus.<br>
 +
</p>
 +
<br>
 +
<br>
<center>
<center>
-
<img width=800 src="https://static.igem.org/mediawiki/2010/8/8c/Freiburg10_VirusInfectionPlot01.png" alt="" />
+
<table style="text-align: left; width: 862px; height: 533px;" border="0"
 +
cellpadding="2" cellspacing="2">
 +
<tbody>
 +
<tr>
 +
<td style="vertical-align: top;"><img
 +
src="https://static.igem.org/mediawiki/2010/8/8c/Freiburg10_VirusInfectionPlot01.png"
 +
alt="" width="800"><br>
 +
<span style="font-weight: bold;">Figure 3:</span> Time course of
 +
concentrations. Extracellular virus particles (blue), free receptors on
 +
cell surface (green), virus-receptor complex (red), dimerized
 +
receptor-virus complex (turquoise), intranuclear virus (purple),
 +
single-stranded DNA (yellow) and enzyme in cytoplasm (black).<br>
 +
</td>
 +
</tr>
 +
</tbody>
 +
</table>
</center>
</center>
-
<br><br>
+
<br>
-
 
+
<p style="text-align: justify;">
 +
As mentioned in the previous paragraph concerning the model extensions
 +
the internalization rate as well as the amount of empty and for this
 +
reason non-functional virus particles strongly depends on the degree of
 +
modification. Hence one of the major goals of simulation was to
 +
calculate the optimal ratio between wild type and modified capsid
 +
proteins on the virus surface.<br>
 +
The simulation result can be seen in figure 4 B where, depending on the
 +
assumed functional behavior, the best targeting efficiency was reached
 +
with a 60% modified capsid. <br>
 +
</p>
 +
<br>
 +
<div style="text-align: center;">
<center>
<center>
-
<img width=800 src="https://static.igem.org/mediawiki/2010/9/9e/Freiburg10_VirusInfectionModPlot01.png" alt="" />
+
<table
 +
style="text-align: left; width: 903px; margin-left: auto; margin-right: auto; height: 880px;"
 +
border="0" cellpadding="2" cellspacing="2">
 +
<tbody>
 +
<tr>
 +
<td style="vertical-align: top;"><img
 +
src="https://static.igem.org/mediawiki/2010/9/9e/Freiburg10_VirusInfectionModPlot01.png"
 +
alt="" width="800"><br>
 +
<span style="font-weight: bold;">Figure 4: A:</span> Enzyme
 +
concentration for different&nbsp; degrees of modification (<span
 +
style="font-style: italic;">m</span>) - the steady state depends on <span
 +
style="font-style: italic;">m</span>. <span style="font-weight: bold;">B:</span>
 +
Enyme concentration as a function of <span style="font-style: italic;">m</span>.
 +
The optimum is reached for <span style="font-style: italic;">m=60%.</span><br>
 +
</td>
 +
</tr>
 +
</tbody>
 +
</table>
</center>
</center>
-
<br><br>
+
</div>
-
 
+
<br>
-
<h3>Model 2 (MM)</h3>
+
<br>
-
<center>
+
-
<img width=800 src="https://static.igem.org/mediawiki/2010/2/24/Freiburg10_VirusInfectionPlot02.png" alt="" />
+
-
</center>
+
-
<br><br>
+
-
 
+
-
 
+
-
<center>
+
-
<img width=800 src="https://static.igem.org/mediawiki/2010/a/a7/Freiburg10_VirusInfectionModPlot02.png" alt="" />
+
-
</center>
+
-
 
+
-
 
+
</html>
</html>
-
 
{{:Team:Freiburg_Bioware/Footer}}
{{:Team:Freiburg_Bioware/Footer}}

Latest revision as of 02:40, 28 October 2010

Model for Virus Infection

As described in the previous model for the virus production we established a model of ordinary differential equation (ODE) based on the law of mass action. The following paragraph depicts the reaction scheme and our model assumptions. In the subsequent paragraphs the system of differential equations is specified and the implementation in MathWorks® MATLAB is discussed. Furthermore reasonable model extensions are presented.
The last section deals with our modeling results and gives an outlook to further modeling applications.

Reaction Scheme

To simplify the mathematical description of the reaction scheme we divide the cell into four compartments: the extracellular matrix (all quantities with the index ext), the space in endosomes (end), the cytoplasm (cyt) and the nucleus (nuc).
A target cell is transduced by viral particles (V) in the extracellular matrix depending on their degree of modification (m) and thus their specificity they can bind to receptors (R) on the cell surface. Once a receptor has formed a complex with a virus particle (VR) receptor dimerization (R2) occurs and the whole complex (VR2) is internalized and endosomes are formed. The virus is released from the endosome to the cytoplasm and is transported to the nucleus where uncoating of the capsid is initiated and single stranded DNA (ssDNA) is released.
Finally viral mRNA is processed and transported into the cytoplasm and the enzyme for therapeutic approach (E) is produced.

Reaction scheme for the virus production

Figure 1: Schematical overview of virus infection
Reaction scheme for the virus production


Reduced Reaction Scheme

This modeling approach neglects the fastest linear transport processes and describes the receptor binding and dimerization in terms of the law of mass action (LMA).

Reaction scheme for the virus production


Differential Equations

The ODE model consists of 7 equations containing 9 rate constants. To the equation for the extracellular virus concentration a degradation term with rate constant (k7,1) is added corresponding to the immune response of the target system. The temporal behavior of ssDNA is supplemented by a linear degradation term with rate constant (k6,1).

Reaction scheme for the virus production


Model Extensions

Since internalization depends on the degree of modification (m) the model is extended by a non-linear dependency in the rate constants. Additionally a decreasing production efficiency reduces the amount of available ssDNA in the nucleus. The estimated functional dependencies are shown in the figure below.


Figure 2: Assumed functions for internalization (k(m), blue) and fraction of infectious virus particles (r(m), red) dependant on modification degree (m).
Regarding the k(m) we assumed that a minimal modification degree of 10% is necessary for virus receptor complex formation due to the fact that at least two binding sites are required for receptor dimerization. For higher m-values the targeting efficiency saturates because enough binding sites exists that every dimerized receptor can build a complex with a virus.

The values of r(m) are normalized to the ratio of infectious particles to empty virus capsids and therefore start at one. Increasing m results in a decreasing amount of infectious particles until a minimal value is reached for the completely modified virus (m=100%).

Future experiments could provide a better understanding of those functional dependencies and therefore improve the model predictions concerning the targeting efficiency discussed later in this chapter.


Methods and Simulation

The ODE model was implemented in MathWorks® MATLAB R2010b. Integration of the differential equations was achieved using the stiff integration routine ode15s with automatic integration step size management.
For initial conditions we took virus concentrations used in the experiments. Rate constants were estimated according to published values. The entire internalization process of viral particles being transported from the extracellular matrix through the cytoplasm into the nuleus happens within seconds (Seisenberger et al. 2001; Bartlett et al. 2000) whereas protein synthetization, capsid formation and virus packaging is taking much more time.
Furthermore experiments have been conducted in order to determine typical time scales corresponding to biological dynamics.

The used parameters are given in the table below.



Table 1: rate constants for the virus infection model


get the m-File (MATLAB source code)!


Results and Discussion

Figure 3 shows the time course of the virus infection model. Starting with an initial virus concentration of 107 viral particles per mililiter the curve decays fast as the infectious particles are bound to receptors and imported to endosomes in the cytoplasm. Since recepor-virus complexes are built the amout of free recepors on the cell surface decreases while the complex concentrations are rising.
As soon as the first particles reach the nucleus and release their viral DNA the intranuclear concentration of ssDNA rises and enzyme transcription is initialized. Therefore the enzyme concentration in cytoplasm is increased until a steady state is reached due to the fact that ssDNA is degraded inside the nucleus.




Figure 3: Time course of concentrations. Extracellular virus particles (blue), free receptors on cell surface (green), virus-receptor complex (red), dimerized receptor-virus complex (turquoise), intranuclear virus (purple), single-stranded DNA (yellow) and enzyme in cytoplasm (black).

As mentioned in the previous paragraph concerning the model extensions the internalization rate as well as the amount of empty and for this reason non-functional virus particles strongly depends on the degree of modification. Hence one of the major goals of simulation was to calculate the optimal ratio between wild type and modified capsid proteins on the virus surface.
The simulation result can be seen in figure 4 B where, depending on the assumed functional behavior, the best targeting efficiency was reached with a 60% modified capsid.



Figure 4: A: Enzyme concentration for different  degrees of modification (m) - the steady state depends on m. B: Enyme concentration as a function of m. The optimum is reached for m=60%.