Location: Dixit, Perelson, 2004 @ e66b6263c530 / dixit_perelson_2004.cellml

Author:
pmr2.import <nobody@models.cellml.org>
Date:
2006-09-03 23:34:04+12:00
Desc:
committing version01 of dixit_perelson_2004
Permanent Source URI:
https://models.cellml.org/workspace/dixit_perelson_2004/rawfile/e66b6263c530c891e1c92cb8c989bad13296ddf2/dixit_perelson_2004.cellml

<?xml version='1.0' encoding='utf-8'?>
<!--  FILE :  dixit_model_I_2004.xml

CREATED :  12th December 2003

LAST MODIFIED : 12th December 2003

AUTHOR :  Catherine Lloyd
          Bioengineering Institute
          The University of Auckland
          
MODEL STATUS :  This model conforms to the CellML 1.0 Specification released on
10th August 2001, and the 16/1/02 CellML Metadata 1.0 Specification.

DESCRIPTION :  This file contains a CellML description of Dixit and Perelson's 1st 2004 mathematical model of viral dynamics in the presence of reverse transcriptase inhibitors.
CHANGES:  
  
--><model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" cmeta:id="dixit_perelson_2004_version01" name="dixit_perelson_2004_version01">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Modelling the Complex Patterns of Viral Load Decay under Antiretroviral Therapy</title>
  <author>
    <firstname>Catherine</firstname>
          <surname>Lloyd</surname>
    <affiliation>
      <shortaffil>Bioengineering Institute, University of Auckland</shortaffil>
    </affiliation>
  </author>
</articleinfo>
  <section id="sec_status">
    <title>Model Status</title>
    <para>
            This is the original unchecked version of the model imported from the previous
            CellML model repository, 24-Jan-2006.
          </para>
  </section>
  <sect1 id="sec_structure">
<title>Model Structure</title>

<para>
When treated with antiretroviral therapy, the plasma viral loads in HIV patients decline in three distinct phases:
</para>
<itemizedlist>
  <listitem>
            <para>When the treatment is initially administered there is a 6 to 72 hour delay before the viral loads begin to decline;</para>
          </listitem>
  <listitem>
            <para>following this <emphasis>shoulder</emphasis> period viral load decline rapidly for about one week;</para>
          </listitem>
  <listitem>
            <para>finally, viral loads continue to decrease, but at a slower rate.</para>
          </listitem>
</itemizedlist>
<para>
Mathematical models which describe this process have been developed.  For more details please see:
</para>
<itemizedlist>
  <listitem>
            <para>
              <ulink url="${HTML_EXMPL_HERZ_MODEL}">Herz <emphasis>et al.</emphasis>, Modelling Viral Dynamics <emphasis>In Vivo</emphasis>, 1996</ulink>;</para>
          </listitem>
  <listitem>
            <para>
              <ulink url="${HTML_EXMPL_NOWAK_MODEL}">Nowak and Bangham, Modelling the Population Dynamics of Immune Responses to Persistent Viruses, 1996</ulink>;</para>
          </listitem>
  <listitem>
            <para>
              <ulink url="${HTML_EXMPL_PERELSON_MODEL}">Perelson <emphasis>et al.</emphasis>, Modelling HIV-1 Dynamics <emphasis>In Vivo</emphasis>, 1996</ulink>;</para>
          </listitem>
  <listitem>
            <para>
              <ulink url="${HTML_EXMPL_GROSSMAN_MODEL}">Grossman <emphasis>et al.</emphasis>, Modelling the Kinetics of Plasma Virus Following the Initiation of Therapy, 1998</ulink>;</para>
          </listitem>
  <listitem>
            <para>
              <ulink url="${HTML_EXMPL_MITTLER_MODEL}">Mittler <emphasis>et al.</emphasis>, Influence of Delayed Viral Production on Viral Dynamics in HIV-1 Infected Patients, 1998</ulink>;</para>
          </listitem>
  <listitem>
            <para>
              <ulink url="${HTML_EXMPL_WODARZ_MODEL}">Wodarz and Nowak, Modelling the Interaction Between HIV and the Immune System, 1999</ulink>; and </para>
          </listitem>
  <listitem>
            <para>
              <ulink url="${HTML_EXMPL_BONHOEFFER_MODEL}">Bonhoeffer <emphasis>et al.</emphasis> Modelling the Population Dynamics of Virus Infected Cells, 2000</ulink>.</para>
          </listitem>
</itemizedlist>

<para>
Despite these models the shoulder phase, where viral loads remain virtually unchanged by the drug therapy, is poorly understood.  Two factors are known to contribute to this delay: Firstly, a pharmacological delay which is the period between the drug administration and the drug action.  This delay is due to the time it takes for the drug to be absorbed from the gut into the blood stream, and then carried to and taken up by the target cells.  The second factor is an intracellular delay which is due to the finite time required for an infected cell to replicate the virus.   
</para>

<para>
In the Dixit and Perelson 2004 publication described here, the authors present a mathematical model that combines pharmacokinetics and viral dynamics and includes intracellular delay.  The model of viral dynamics is based on that published by <ulink url="${HTML_EXMPL_PERELSON_MODEL}">Perelson <emphasis>et al.</emphasis>
          </ulink> (see <xref linkend="fig_cell_diagram"/> below).  To accurately represent the pharmacokinetics they assume a two compartment model consisting of blood and cells.  The models of viral dynamics and pharmacokinetics are connected via a parameter which represents drug specific efficacy.  Because reverse transcriptase inhibitors (RTIs) and protease inhibitors (PIs) have very different effects on the duration of the intracellular delay, they have been modelled separately (see <xref linkend="fig_reaction_diagram1"/> and <xref linkend="fig_reaction_diagram2"/>).  </para>

<para>
The models have been described here in CellML (the raw CellML description of the Dixit and Perelson 2004 models can be downloaded in various formats as described in <xref linkend="sec_download_this_model"/>).     
</para>

<para>
The complete original paper reference is cited below:
</para>

<para>
<ulink url="http://www.sciencedirect.com/science?_ob=ArticleURL&amp;_udi=B6WMD-4B1WXVT-6&amp;_user=140507&amp;_coverDate=01%2F07%2F2004&amp;_alid=134704688&amp;_rdoc=1&amp;_fmt=summary&amp;_orig=search&amp;_qd=1&amp;_cdi=6932&amp;_sort=d&amp;_docanchor=&amp;view=c&amp;_acct=C000011498&amp;_version=1&amp;_urlVersion=0&amp;_userid=140507&amp;md5=574df08b1a72cfc52ed9a2bb186e39a5">Complex patterns of viral load decay under antiretroviral therapy: influence of pharmacokinetics and intracelllar delay</ulink>, Narendra M. Dixit and Alan S. Perelson, 2004, <ulink url="http://www.sciencedirect.com/science?_ob=JournalURL&amp;_cdi=6932&amp;_auth=y&amp;_acct=C000011498&amp;_version=1&amp;_urlVersion=0&amp;_userid=140507&amp;md5=c8e64d20ff203cda2255f4b6b7affe4f">
            <emphasis>Journal of Theoretical Biology</emphasis>
          </ulink>, 226, 95-109.  (<ulink url="http://www.sciencedirect.com/science?_ob=ArticleURL&amp;_udi=B6WMD-4B1WXVT-6&amp;_coverDate=01%2F07%2F2004&amp;_alid=134704688&amp;_rdoc=1&amp;_fmt=&amp;_orig=search&amp;_qd=1&amp;_cdi=6932&amp;_sort=d&amp;view=c&amp;_acct=C000011498&amp;_version=1&amp;_urlVersion=0&amp;_userid=140507&amp;md5=d77dd62b0a4028ce41dd953455d4bf3c">Full text (HTML)</ulink> and PDF versions of the article are available on the <emphasis>Journal of Theoretical Biology</emphasis> website.)  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&amp;db=PubMed&amp;list_uids=14637059&amp;dopt=Abstract">PubMed ID: 14637059</ulink>
</para>

<informalfigure float="0" id="fig_cell_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>cell diagram</title>
    </objectinfo>
    <imagedata fileref="cell_diagram.gif"/>
  </imageobject>
</mediaobject>
<caption>Schematic summary of the dynamics of HIV-1 infection <emphasis>in vivo</emphasis>.  This model is based on that published by Perelson <emphasis>et al.</emphasis> in 1996.</caption>
</informalfigure>

<informalfigure float="0" id="fig_reaction_diagram1">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>reaction diagram1</title>
    </objectinfo>
    <imagedata fileref="reaction_diagram.gif"/>
  </imageobject>
</mediaobject>
<caption>Pharmacokinetics - the two compartment model for protease inhibitors.</caption>
</informalfigure>

<informalfigure float="0" id="fig_reaction_diagram2">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>reaction diagram2</title>
    </objectinfo>
    <imagedata fileref="reaction_diagram2.gif"/>
  </imageobject>
</mediaobject>
<caption>Pharmacokinetics - the two compartment model for reverse transcriptase inhibitors.  The main difference between the models is that RTIs must be phosphorylated within the cell in order to be in their active form.</caption>
</informalfigure>

</sect1>
</article>
</documentation>
  
  
  
  <units name="ml">
    <unit units="litre" prefix="milli"/>
  </units>
  
  <units name="mg">
    <unit units="gram" prefix="milli"/>
  </units>
  
  <units name="per_ml">
    <unit units="litre" prefix="milli" exponent="-1"/>
  </units>
  
  <units name="mg_per_ml">
    <unit units="gram" prefix="milli"/>
    <unit units="litre" prefix="milli" exponent="-1"/>
  </units>
  
  <units name="day">
    <unit units="second" multiplier="86400.0"/>
  </units>
  
  <units name="flux">
    <unit units="ml"/> 
    <unit units="day" exponent="-1"/>
  </units>
  
  <units name="first_order_rate_constant">
    <unit units="day" exponent="-1"/>
  </units>
  
  <units name="second_order_rate_constant">
    <unit units="ml" exponent="-1"/>
    <unit units="day" exponent="-1"/>
  </units>

  
  <component name="environment">
    <variable units="day" public_interface="out" name="time"/>
  </component>
  
  <component cmeta:id="T" name="T">
    <rdf:RDF>
      <rdf:Description rdf:about="T">
        <dc:title>T</dc:title>
        <dcterms:alternative>uninfected target CD4 cells</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="per_ml" public_interface="out" name="T" initial_value="1E6"/>
    
    <variable units="second_order_rate_constant" name="lamda" initial_value="1E4"/>
    <variable units="first_order_rate_constant" name="d" initial_value="0.01"/>
    
    <variable units="flux" public_interface="in" name="k"/>
    <variable units="per_ml" public_interface="in" name="V"/>
    <variable units="dimensionless" public_interface="in" name="epsilon_RTI"/>
    <variable units="day" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>T</ci> 
        </apply>
        <apply>
          <minus/>
          <ci> lamda </ci>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> d </ci>
              <ci> T </ci>
            </apply>
            <apply>
              <times/>
              <apply>
                <minus/>
                <cn cellml:units="dimensionless"> 1.0 </cn>
                <ci> epsilon_RTI </ci>
              </apply>
              <ci> k </ci> 
              <ci> T </ci> 
              <ci> V </ci> 
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="T_" name="T_">
    <rdf:RDF>
      <rdf:Description rdf:about="T_">
        <dc:title>T_</dc:title>
        <dcterms:alternative>productively infected CD4 cells</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="per_ml" public_interface="out" name="T_"/>
    
    <variable units="first_order_rate_constant" name="tau" initial_value="1.5"/>
    <variable units="first_order_rate_constant" name="m" initial_value="0.01"/>
    
    <variable units="flux" public_interface="in" name="k"/>
    <variable units="per_ml" public_interface="in" name="V"/>
    <variable units="per_ml" public_interface="in" name="T"/>
    <variable units="dimensionless" public_interface="in" name="epsilon_RTI"/>
    <variable units="first_order_rate_constant" public_interface="in" name="delta"/>
    <variable units="day" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>T_</ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <apply>
              <minus/>
              <cn cellml:units="dimensionless"> 1.0 </cn>
              <apply>
                <times/>
                <ci> epsilon_RTI </ci>
                <apply>
                  <minus/>
                  <ci> time </ci>
                  <ci> tau </ci>
                </apply>
              </apply>
            </apply>
            <ci> k </ci> 
            <ci> T </ci>
            <apply>
              <minus/>
              <ci> time </ci>
              <ci> tau </ci>
            </apply> 
            <ci> V </ci>
            <apply>
              <minus/>
              <ci> time </ci>
              <ci> tau </ci>
            </apply>
            <apply>
              <exp/>
              <apply>
                <times/>
                <apply>
                  <minus/>
                  <ci> m </ci>  
                </apply>
                <ci> tau </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> delta </ci> 
            <ci> T_ </ci>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="V" name="V">
    <rdf:RDF>
      <rdf:Description rdf:about="V">
        <dc:title>V</dc:title>
        <dcterms:alternative>viral load</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="per_ml" public_interface="out" name="V"/>
    
    <variable units="dimensionless" name="N" initial_value="2500.0"/>
    <variable units="first_order_rate_constant" name="c" initial_value="23.0"/>
    
    <variable units="per_ml" public_interface="in" name="T_"/>
    <variable units="first_order_rate_constant" public_interface="in" name="delta"/>
    <variable units="day" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>V</ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> N </ci>
            <ci> delta </ci>
            <ci> T_ </ci>
          </apply>
          <apply>
            <times/>
            <ci> c </ci>
            <ci> V </ci> 
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component cmeta:id="epsilon_RTI" name="epsilon_RTI">
    <rdf:RDF>
      <rdf:Description rdf:about="epsilon_RTI">
        <dc:title>instantaneous drug efficacy</dc:title>
        <dcterms:alternative>non-infectious viral load</dcterms:alternative>  
      </rdf:Description>
    </rdf:RDF>  
    <variable units="dimensionless" public_interface="out" name="epsilon_RTI"/>
    
    <variable units="mg_per_ml" name="IC50" initial_value="0.54"/>
    
    <variable units="mg_per_ml" public_interface="in" name="Cc"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>epsilon_RTI</ci> 
        <apply>
          <divide/>
          <ci> Cc </ci>
          <apply>
            <plus/>
            <ci> IC50 </ci>
            <ci> Cc </ci>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="Mg">
    <variable units="mg_per_ml" public_interface="out" name="Mg"/>
    
    <variable units="first_order_rate_constant" public_interface="in" name="ka"/>
    <variable units="day" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>Mg</ci> 
        </apply>
        <apply>
          <times/>
          <apply>
            <minus/>
            <ci> ka </ci>
          </apply>
          <ci> Mg </ci>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="Mb">
    <variable units="mg_per_ml" public_interface="out" name="Mb"/>
    
    <variable units="first_order_rate_constant" public_interface="in" name="ka"/>
    <variable units="first_order_rate_constant" public_interface="in" name="ke"/>
    <variable units="day" public_interface="in" name="time"/>
    <variable units="mg_per_ml" public_interface="in" name="Mg"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>Mb</ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> ka </ci>
            <ci> Mg </ci>
          </apply>
          <apply>
            <times/>
            <ci> ke </ci>
            <ci> Mb </ci>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="Cb">
    <variable units="mg_per_ml" public_interface="out" name="Cb"/>
    
    <variable units="ml" name="Vd" initial_value="87500.0"/>
    <variable units="dimensionless" name="F" initial_value="0.39"/>
    <variable units="mg" name="D" initial_value="300.0"/>
    
    <variable units="first_order_rate_constant" public_interface="in" name="ka"/>
    <variable units="first_order_rate_constant" public_interface="in" name="ke"/>
    <variable units="mg_per_ml" public_interface="in" name="Mb"/>
    <variable units="day" public_interface="in" name="time"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>Cb</ci> 
        <apply>
          <times/>
          <apply>
            <divide/>
            <apply>
              <times/>
              <ci> F </ci>
              <ci> D </ci>
            </apply>
            <ci> Vd </ci>
          </apply>
          <apply>
            <divide/>
            <ci> ka </ci>
            <apply>
              <minus/>
              <ci> ke </ci>
              <ci> ka </ci>
            </apply>
          </apply>
          <apply>
            <minus/>
            <apply>
              <exp/>
              <apply>
                <times/>
                <apply>
                  <minus/>
                  <ci> ka </ci>
                </apply>
                <ci> time </ci>
              </apply>
            </apply>
            <apply>
              <exp/>
              <apply>
                <times/>
                <apply>
                  <minus/>
                  <ci> ke </ci>
                </apply>
                <ci> time </ci>
              </apply>
            </apply>
          </apply>
        </apply>
      </apply>
      
      <apply>
        <eq/> 
        <ci> Vd </ci> 
        <apply>
          <divide/> 
          <ci>Mb</ci>  
          <ci>Cb</ci>
        </apply>
      </apply>   
    </math>
  </component>
  
  <component name="Cx">
    <variable units="mg_per_ml" public_interface="out" name="Cx"/>
    
    <variable units="dimensionless" name="H" initial_value="1800.0"/>
    <variable units="dimensionless" name="fb" initial_value="0.07"/>
    
    <variable units="mg_per_ml" public_interface="in" name="Cb"/>
    <variable units="mg_per_ml" public_interface="in" name="Cc"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <ci>Cx</ci> 
        <piecewise>
          <piece>
            <apply>
              <minus/>
              <apply>
                <times/>
                <apply>
                  <minus/>
                  <cn cellml:units="dimensionless"> 1.0 </cn>
                  <ci> fb </ci>
                </apply>
                <ci> H </ci>
                <ci> Cb </ci>
              </apply>
              <ci> Cc </ci>
            </apply>
            <apply>
              <gt/>
              <apply>
                <minus/>
                <apply>
                  <times/>
                  <apply>
                    <minus/>
                    <cn cellml:units="dimensionless"> 1.0 </cn>
                    <ci> fb </ci>
                  </apply>
                  <ci> H </ci>
                  <ci> Cb </ci>
                </apply>
                <ci> Cc </ci>
              </apply>
              <cn cellml:units="mg_per_ml"> 0.0 </cn>
            </apply>
          </piece>
          <otherwise>
            <cn cellml:units="mg_per_ml"> 0.0 </cn>
          </otherwise>
        </piecewise>
      </apply>
    </math>
  </component>
  
  <component name="Cc">
    <variable units="mg_per_ml" public_interface="out" name="Cc"/>
    
    <variable units="first_order_rate_constant" public_interface="in" name="k1f"/>
    <variable units="first_order_rate_constant" public_interface="in" name="k1b"/>
    <variable units="first_order_rate_constant" public_interface="in" name="kacell"/>
    <variable units="first_order_rate_constant" public_interface="in" name="kecell"/>
    <variable units="day" public_interface="in" name="time"/>
    <variable units="mg_per_ml" public_interface="in" name="Cx"/>
    <variable units="mg_per_ml" public_interface="in" name="Ccp"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>Cc</ci> 
        </apply>
        <apply>
          <plus/>
          <apply>
            <minus/>
            <apply>
              <times/>
              <ci> kacell </ci>
              <ci> Cx </ci>
            </apply>
            <apply>
              <plus/>
              <apply>
                <times/>
                <ci> kecell </ci>
                <ci> Cc </ci>
              </apply>
              <apply>
                <times/>
                <ci> k1f </ci>
                <ci> Cc </ci>
              </apply>
            </apply>
          </apply>
          <apply>
            <times/>
            <ci> k1b </ci>
            <ci> Ccp </ci>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="Ccp">
    <variable units="mg_per_ml" public_interface="out" name="Ccp"/>
    
    <variable units="first_order_rate_constant" public_interface="in" name="k1f"/>
    <variable units="first_order_rate_constant" public_interface="in" name="k1b"/>
    <variable units="first_order_rate_constant" public_interface="in" name="kecell"/>
    <variable units="first_order_rate_constant" public_interface="in" name="k2f"/>
    <variable units="first_order_rate_constant" public_interface="in" name="k2b"/>
    <variable units="day" public_interface="in" name="time"/>
    <variable units="mg_per_ml" public_interface="in" name="Cc"/>
    <variable units="mg_per_ml" public_interface="in" name="Ccpp"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>Ccp</ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> k1f </ci>
              <ci> Cc </ci>
            </apply>
            <apply>
              <times/>
              <ci> k2b </ci>
              <ci> Ccpp </ci>
            </apply>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> kecell </ci>
              <ci> Ccp </ci>
            </apply>
            <apply>
              <times/>
              <ci> k1b </ci>
              <ci> Ccp </ci>
            </apply>
            <apply>
              <times/>
              <ci> k2f </ci>
              <ci> Ccp </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
  
  <component name="Ccpp">
    <variable units="mg_per_ml" public_interface="out" name="Ccpp"/>
    
    <variable units="first_order_rate_constant" public_interface="in" name="kecell"/>
    <variable units="first_order_rate_constant" public_interface="in" name="k2f"/>
    <variable units="first_order_rate_constant" public_interface="in" name="k2b"/>
    <variable units="day" public_interface="in" name="time"/>
    <variable units="mg_per_ml" public_interface="in" name="Ccp"/>
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci>Ccpp</ci> 
        </apply>
        <apply>
          <minus/>
          <apply>
            <times/>
            <ci> k2f </ci>
            <ci> Ccp </ci>
          </apply>
          <apply>
            <plus/>
            <apply>
              <times/>
              <ci> kecell </ci>
              <ci> Ccpp </ci>
            </apply>
            <apply>
              <times/>
              <ci> k2b </ci>
              <ci> Ccpp </ci>
            </apply>
          </apply>
        </apply>
      </apply>    
    </math>
  </component>
 
  <component name="kinetic_parameters"> 
    <variable units="flux" public_interface="out" name="k" initial_value="2.4E-8"/>
    <variable units="first_order_rate_constant" public_interface="out" name="delta"/>
    <variable units="first_order_rate_constant" public_interface="out" name="ka" initial_value="14.64"/>
    <variable units="first_order_rate_constant" public_interface="out" name="ke" initial_value="6.86"/>
    <variable units="first_order_rate_constant" public_interface="out" name="kacell" initial_value="24000.0"/>
    <variable units="first_order_rate_constant" public_interface="out" name="kecell" initial_value="1.1"/>
    <variable units="first_order_rate_constant" public_interface="out" name="k1f" initial_value="9.6"/>
    <variable units="first_order_rate_constant" public_interface="out" name="k1b" initial_value="30.3"/>
    <variable units="first_order_rate_constant" public_interface="out" name="k2f" initial_value="270.7"/>
    <variable units="first_order_rate_constant" public_interface="out" name="k2b" initial_value="95.5"/>
  </component>
  
  <connection>
    <map_components component_2="environment" component_1="T"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="T_"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="V"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="Mg"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="Mb"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="Cb"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="Cc"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="Ccp"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="environment" component_1="Ccpp"/>
    <map_variables variable_2="time" variable_1="time"/>
  </connection>
  
  <connection>
    <map_components component_2="T_" component_1="T"/>
    <map_variables variable_2="T" variable_1="T"/>
  </connection>
  
  <connection>
    <map_components component_2="V" component_1="T"/>
    <map_variables variable_2="V" variable_1="V"/>
  </connection>
  
  <connection>
    <map_components component_2="V" component_1="T_"/>
    <map_variables variable_2="V" variable_1="V"/>
    <map_variables variable_2="T_" variable_1="T_"/>
  </connection>
  
  <connection>
    <map_components component_2="kinetic_parameters" component_1="T"/>
    <map_variables variable_2="k" variable_1="k"/>
  </connection>
  
  <connection>
    <map_components component_2="kinetic_parameters" component_1="T_"/>
    <map_variables variable_2="delta" variable_1="delta"/>
    <map_variables variable_2="k" variable_1="k"/>
  </connection>
  
  <connection>
    <map_components component_2="epsilon_RTI" component_1="T"/>
    <map_variables variable_2="epsilon_RTI" variable_1="epsilon_RTI"/>
  </connection>
  
  <connection>
    <map_components component_2="epsilon_RTI" component_1="T_"/>
    <map_variables variable_2="epsilon_RTI" variable_1="epsilon_RTI"/>
  </connection>
  
  <connection>
    <map_components component_2="Cc" component_1="epsilon_RTI"/>
    <map_variables variable_2="Cc" variable_1="Cc"/>
  </connection>
  
  <connection>
    <map_components component_2="kinetic_parameters" component_1="V"/>
    <map_variables variable_2="delta" variable_1="delta"/>
  </connection>
  
  <connection>
    <map_components component_2="kinetic_parameters" component_1="Mg"/>
    <map_variables variable_2="ka" variable_1="ka"/>
  </connection>
  
  <connection>
    <map_components component_2="kinetic_parameters" component_1="Mb"/>
    <map_variables variable_2="ka" variable_1="ka"/>
    <map_variables variable_2="ke" variable_1="ke"/>
  </connection>
  
  <connection>
    <map_components component_2="Mg" component_1="Mb"/>
    <map_variables variable_2="Mg" variable_1="Mg"/>
  </connection>
  
  <connection>
    <map_components component_2="Mb" component_1="Cb"/>
    <map_variables variable_2="Mb" variable_1="Mb"/>
  </connection>
  
  <connection>
    <map_components component_2="kinetic_parameters" component_1="Cb"/>
    <map_variables variable_2="ka" variable_1="ka"/>
    <map_variables variable_2="ke" variable_1="ke"/>
  </connection>
  
  <connection>
    <map_components component_2="Cx" component_1="Cc"/>
    <map_variables variable_2="Cx" variable_1="Cx"/>
    <map_variables variable_2="Cc" variable_1="Cc"/>
  </connection>
  
  <connection>
    <map_components component_2="Cb" component_1="Cx"/>
    <map_variables variable_2="Cb" variable_1="Cb"/>
  </connection>
  
  <connection>
    <map_components component_2="Ccp" component_1="Cc"/>
    <map_variables variable_2="Ccp" variable_1="Ccp"/>
    <map_variables variable_2="Cc" variable_1="Cc"/>
  </connection>
  
  <connection>
    <map_components component_2="Ccpp" component_1="Ccp"/>
    <map_variables variable_2="Ccp" variable_1="Ccp"/>
    <map_variables variable_2="Ccpp" variable_1="Ccpp"/>
  </connection>
  
  <connection>
    <map_components component_2="kinetic_parameters" component_1="Cc"/>
    <map_variables variable_2="kacell" variable_1="kacell"/>
    <map_variables variable_2="kecell" variable_1="kecell"/>
    <map_variables variable_2="k1f" variable_1="k1f"/>
    <map_variables variable_2="k1b" variable_1="k1b"/>
  </connection>
  
  <connection>
    <map_components component_2="kinetic_parameters" component_1="Ccp"/>
    <map_variables variable_2="kecell" variable_1="kecell"/>
    <map_variables variable_2="k1f" variable_1="k1f"/>
    <map_variables variable_2="k1b" variable_1="k1b"/>
    <map_variables variable_2="k2f" variable_1="k2f"/>
    <map_variables variable_2="k2b" variable_1="k2b"/>
  </connection>
  
  <connection>
    <map_components component_2="kinetic_parameters" component_1="Ccpp"/>
    <map_variables variable_2="kecell" variable_1="kecell"/>
    <map_variables variable_2="k2f" variable_1="k2f"/>
    <map_variables variable_2="k2b" variable_1="k2b"/>
  </connection>
  





<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
  <rdf:Seq rdf:about="rdf:#citationAuthorsSeq">
    <rdf:li rdf:resource="rdf:#author1Vcard"/>
    <rdf:li rdf:resource="rdf:#author2Vcard"/>
  </rdf:Seq>
  <rdf:Description rdf:about="rdf:#author2VcardN">
    <vCard:Given>Alan</vCard:Given>
    <vCard:Family>Perelson</vCard:Family>
    <vCard:Other>S</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#11c31409-a448-48fc-b64a-f546ee0b7345">
    <dcterms:W3CDTF>2004-01-07</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#f9bc7305-5eef-430d-acdf-2448bc821e28">
    <dc:subject rdf:resource="rdf:#4c233164-4c77-4298-9096-8d07bebaa418"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b0ce441e-f8af-4f2d-b839-41ca0bca262d">
    <vCard:Orgname>The University of Auckland</vCard:Orgname>
    <vCard:Orgunit>The Bioengineering Institute</vCard:Orgunit>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#9f3cd882-1bce-47fc-bb76-e60f562b033a">
    <vCard:ORG rdf:resource="rdf:#b0ce441e-f8af-4f2d-b839-41ca0bca262d"/>
    <vCard:EMAIL rdf:resource="rdf:#8b960d83-8355-463d-b79d-20313b2b7172"/>
    <vCard:N rdf:resource="rdf:#678e8bb9-f403-4ddf-b4ac-07100e721b14"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#950f1196-c576-4697-96fe-82ed03d765ce">
    <dc:creator rdf:resource="rdf:#b0e4ad90-748a-440d-bfc4-4d5496267620"/>
    <rdf:value>This is the CellML description of Dixit and Perelson's 1st 2004 mathematical model of viral dynamics in the presence of reverse transcriptase inhibitors.</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#cd5072f4-573e-410e-ad9e-defc35c67619">
    <dc:creator rdf:resource="rdf:#citationAuthorsSeq"/>
    <dc:title>Complex patterns of viral load decay under antiretroviral therapy: influence of pharmacokinetics and intracellular delay</dc:title>
    <bqs:volume>226</bqs:volume>
    <bqs:first_page>95</bqs:first_page>
    <bqs:Journal rdf:resource="rdf:#ca3e8aa9-ce20-45b0-b377-5f3702bd9f43"/>
    <dcterms:issued rdf:resource="rdf:#11c31409-a448-48fc-b64a-f546ee0b7345"/>
    <bqs:last_page>109</bqs:last_page>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#8b960d83-8355-463d-b79d-20313b2b7172">
    <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
    <rdf:value>c.lloyd@auckland.ac.nz</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b0e4ad90-748a-440d-bfc4-4d5496267620">
    <vCard:FN>Catherine Lloyd</vCard:FN>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#author2Vcard">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#author2VcardN"/>
  </rdf:Description>
  <rdf:Description rdf:about="#dixit_perelson_2004_version01">
    <dc:title>
        Dixit and Perelson's 1st 2004 mathematical model of viral dynamics in 
        the presence of reverse transcriptase inhibitors.
      </dc:title>
    <cmeta:comment rdf:resource="rdf:#950f1196-c576-4697-96fe-82ed03d765ce"/>
    <bqs:reference rdf:resource="rdf:#b50e3c9e-a81b-4bb8-ad58-e1073512b923"/>
    <bqs:reference rdf:resource="rdf:#f9bc7305-5eef-430d-acdf-2448bc821e28"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#author1Vcard">
    <rdf:type rdf:resource="http://www.cellml.org/bqs/1.0#Person"/>
    <vCard:N rdf:resource="rdf:#author1VcardN"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#ca3e8aa9-ce20-45b0-b377-5f3702bd9f43">
    <dc:title>Journal of Theoretical Biology</dc:title>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#88e11d82-46f4-4a2f-96c8-3bb095f3f60e">
    <dcterms:W3CDTF>2003-12-12</dcterms:W3CDTF>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#author1VcardN">
    <vCard:Given>Narendra</vCard:Given>
    <vCard:Family>Dixit</vCard:Family>
    <vCard:Other>M</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#678e8bb9-f403-4ddf-b4ac-07100e721b14">
    <vCard:Given>Catherine</vCard:Given>
    <vCard:Family>Lloyd</vCard:Family>
    <vCard:Other>May</vCard:Other>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#4c233164-4c77-4298-9096-8d07bebaa418">
    <bqs:subject_type>keyword</bqs:subject_type>
    <rdf:value>immunology</rdf:value>
  </rdf:Description>
  <rdf:Description rdf:about="">
    <dc:publisher>
        The University of Auckland, Bioengineering Institute
      </dc:publisher>
    <dcterms:created rdf:resource="rdf:#88e11d82-46f4-4a2f-96c8-3bb095f3f60e"/>
    <dc:creator rdf:resource="rdf:#9f3cd882-1bce-47fc-bb76-e60f562b033a"/>
  </rdf:Description>
  <rdf:Description rdf:about="rdf:#b50e3c9e-a81b-4bb8-ad58-e1073512b923">
    <bqs:Pubmed_id>14637059</bqs:Pubmed_id>
    <bqs:JournalArticle rdf:resource="rdf:#cd5072f4-573e-410e-ad9e-defc35c67619"/>
  </rdf:Description>
</rdf:RDF>
</model>