Location: Lockwood, Ewy, Hermann, Holford, 2006 @ c5bc10cf7952 / lockwood_ewy_hermann_holford_2006.cellml

Author:
Catherine Lloyd <c.lloyd@auckland.ac.nz>
Date:
2009-12-02 12:25:21+13:00
Desc:
Corrected units, corrected some equations, had to add a dummy ode to get the model to run in COR because it's purely algebraic. Also added documentation and metadata. Model now runs to almost match the paper but we need to add an equation describing drug input and clearance. This is not explicitly defined in an equation in the paper but it can hopefully be constructed from what is written in the text.
Permanent Source URI:
https://models.physiomeproject.org/workspace/lockwood_ewy_hermann_holford_2006/rawfile/c5bc10cf7952125929be52fc37a0a1eca256ba7c/lockwood_ewy_hermann_holford_2006.cellml

<?xml version="1.0" encoding="utf-8"?>


<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="lockwood_2006" name="lockwood_2006">

<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
  <articleinfo>
  <title>Application of clinical trial simulation to compare proof-of-concept study designs for drugs with a slow onset of effect; an example in Alzheimer's disease</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 CellML model represents version A of the published model (Table 1.A vPFK regulation by ATP-AMP concentrations) and runs in both PCEnv and COR to replicate the published results (Figure 2A).  The units have been checked and they are consistent.  We'd like to thank the original model author Mathieu Cloutier for his time spent curating the CellML model to get it matching the published model.
          </para>
  </section>
  <sect1 id="sec_structure">
  <title>Model Structure</title>

<para>
ABSTRACT: OBJECTIVE: Clinical trial simulation (CTS) was used to select a robust design to test the hypothesis that a new treatment was effective for Alzheimer's disease (AD). Typically, a parallel group, placebo controlled, 12-week trial in 200-400 AD patients would be used to establish drug effect relative to placebo (i.e., Ho: Drug Effect = 0). We evaluated if a crossover design would allow smaller and shorter duration trials. MATERIALS AND METHODS: A family of plausible drug and disease models describing the time course of the AD assessment scale (ADAS-Cog) was developed based on Phase I data and literature reports of other treatments for AD. The models included pharmacokinetic, pharmacodynamic, disease progression, and placebo components. Eight alternative trial designs were explored via simulation. One hundred replicates of each combination of drug and disease model and trial design were simulated. A 'positive trial' reflecting drug activity was declared considering both a dose trend test (p less than 0.05) and pair-wise comparisons to placebo (p less than 0.025). RESULTS: A 4 x 4 Latin Square design was predicted to have at least 80% power to detect activity across a range of drug and disease models. The trial design was subsequently implemented and the trial was completed. Based on the results of the actual trial, a conclusive decision about further development was taken. The crossover design provided enhanced power over a parallel group design due to the lower residual variability. CONCLUSION: CTS aided the decision to use a more efficient proof of concept trial design, leading to savings of up to US 4 M dollars in direct costs and a firm decision 8-12 months earlier than a 12-week parallel group trial.
</para>

<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
  <imageobject>
    <objectinfo>
      <title>model diagram</title>
    </objectinfo>
    <imagedata fileref="lockwood_2006.png"/>
  </imageobject>
</mediaobject>
<caption>Schematic diagram of the model.</caption>
</informalfigure>

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

<para>
Application of clinical trial simulation to compare proof-of-concept study designs for drugs with a slow onset of effect; an example in Alzheimer's disease, Peter Lockwood, Wayne Ewy, David Hermann and Nick Holford, 2006, <emphasis>Pharmaceutical Research</emphasis>, 23, (9), 2050-2059.  <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&amp;cmd=Retrieve&amp;dopt=AbstractPlus&amp;list_uids=16906456&amp;query_hl=1&amp;itool=pubmed_docsum">PubMed ID: 16906456</ulink>
</para>

</sect1>
</article>
</documentation>  
	
    <units name="units" base_units="yes" />
	
	<units name="units_ml_per_ng">
      <unit units="units"/>
	  <unit prefix="milli" units="litre"/>
      <unit exponent="-1" prefix="nano" units="gram"/>
    </units>
	
	<units name="ng_per_ml">
      <unit prefix="nano" units="gram"/>
      <unit exponent="-1" prefix="milli" units="litre"/>
    </units>
	
	<units name="day">
      <unit multiplier="86400" units="second"/>
    </units>
	
	<units name="per_day">
      <unit exponent="-1" units="day"/>
    </units>
	
	<units name="units_per_day">
      <unit units="units"/>
	  <unit exponent="-1" units="day"/>
    </units>
	
	<units name="year">
      <unit multiplier="31536000" units="second"/>
    </units>
	
	<units name="per_year">
      <unit exponent="-1" units="year"/>
    </units>
	
	<component name="environment">
      <variable name="time" public_interface="out" units="day"/>
    </component>
	
	<component name="population_pharmacodynamics_model">
      <variable name="S" units="units"/>
      <variable name="S0" units="units" initial_value="30"/>
      <variable name="alpha" units="units_per_day" initial_value="0.0164"/>
      <variable name="epsilon" units="units" initial_value="16"/>
	  
	  <variable name="time" units="day" public_interface="in"/>
      <variable name="ADAS_Cog_p" units="units" public_interface="in"/>
      <variable name="PD_CeA" units="units" public_interface="in"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>S</ci>
            <apply>
               <plus/>
               <ci>S0</ci>
               <apply>
                  <times/>
                  <ci>alpha</ci>
                  <ci>time</ci>
               </apply>
               <ci>ADAS_Cog_p</ci>
               <ci>PD_CeA</ci>
               <ci>epsilon</ci>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="placebo_response_model">
      <variable name="ADAS_Cog_p" units="units" public_interface="out"/>
      <variable name="beta_P" units="units" initial_value="-3"/>
      <variable name="Keq_p" units="per_day"/>
      <variable name="Kel_p" units="per_day"/>
      <variable name="t_half_el_p" units="day" initial_value="7"/>
      <variable name="t_half_eq_p" units="day" initial_value="6"/>
      
	  <variable name="time" units="day" public_interface="in"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>ADAS_Cog_p</ci>
            <apply>
               <times/>
               <apply>
                  <divide/>
                  <apply>
                     <times/>
                     <ci>beta_P</ci>
                     <ci>Keq_p</ci>
                  </apply>
                  <apply>
                     <minus/>
                     <ci>Keq_p</ci>
                     <ci>Kel_p</ci>
                  </apply>
               </apply>
               <apply>
                  <minus/>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <ci>Kel_p</ci>
                        </apply>
                        <ci>time</ci>
                     </apply>
                  </apply>
                  <apply>
                     <exp/>
                     <apply>
                        <times/>
                        <apply>
                           <minus/>
                           <ci>Keq_p</ci>
                        </apply>
                        <ci>time</ci>
                     </apply>
                  </apply>
               </apply>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>Keq_p</ci>
            <apply>
               <divide/>
               <apply>
                  <ln/>
                  <cn cellml:units="dimensionless">2</cn>
               </apply>
               <ci>t_half_eq_p</ci>
            </apply>
         </apply>
         
		 <apply>
            <eq/>
            <ci>Kel_p</ci>
            <apply>
               <divide/>
               <apply>
                  <ln/>
                  <cn cellml:units="dimensionless">2</cn>
               </apply>
               <ci>t_half_el_p</ci>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="pharmacokinetic_model">
      <variable name="CL" units="dimensionless"/>
      <variable name="smk" units="dimensionless" initial_value="1"/>
      <variable name="age" units="year" initial_value="40"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>CL</ci>
            <apply>
               <times/>
               <cn cellml:units="dimensionless">94.5</cn>
               <apply>
                  <exp/>
                  <apply>
                     <times/>
                     <apply>
                        <minus/>
						<cn cellml:units="per_year">0.0135</cn>
					 </apply>
                     <apply>
                        <minus/>
                        <ci>age</ci>
						<cn cellml:units="year">40</cn>
                     </apply>
                  </apply>
               </apply>
               <ci>smk</ci>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="drop_out_model">
      <variable name="Sv" units="dimensionless"/>
      
	  <variable name="time" units="day" public_interface="in"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>Sv</ci>
            <apply>
               <exp/>
               <apply>
                  <times/>
                  <apply>
                     <minus/>
                     <cn cellml:units="per_day">0.00145</cn>
                  </apply>
                  <ci>time</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="drug_response_model">
      <variable name="PD_CeA" units="units" public_interface="out"/>
      
	  <variable name="beta_a" units="units_ml_per_ng" initial_value="-0.047"/>
      <variable name="CeA" units="ng_per_ml" initial_value="25"/>
      
	  <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>PD_CeA</ci>
            <apply>
               <times/>
               <ci>beta_a</ci>
               <ci>CeA</ci>
            </apply>
         </apply>
      </math>
   </component>
   
   <component name="dummy_ode">
    <variable initial_value="0.2" name="dummy" units="dimensionless"/> 
    
    <variable name="time" public_interface="in" units="day"/> 
    
    <math xmlns="http://www.w3.org/1998/Math/MathML"> 
      <apply>
        <eq/> 
        <apply>
          <diff/> 
          <bvar>
            <ci>time</ci>
          </bvar> 
          <ci> dummy </ci> 
        </apply>
        <apply>
          <plus/> 
          <cn cellml:units="per_day">0.5</cn>
		  <cn cellml:units="per_day">0.6</cn>
        </apply>  
      </apply>
    </math> 
  </component>
   
   <connection>
      <map_components component_1="placebo_response_model" component_2="population_pharmacodynamics_model"/>
      <map_variables variable_1="ADAS_Cog_p" variable_2="ADAS_Cog_p"/>
   </connection>
   <connection>
      <map_components component_1="drug_response_model" component_2="population_pharmacodynamics_model"/>
      <map_variables variable_1="PD_CeA" variable_2="PD_CeA"/>
   </connection>
   <connection>
      <map_components component_1="drop_out_model" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="placebo_response_model" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="population_pharmacodynamics_model" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="dummy_ode" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   
   <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xmlns:bqs="http://www.cellml.org/bqs/1.0#" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:dcterms="http://purl.org/dc/terms/" xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">
    <foo xmlns="http://www.example.org/foo/bar/baz" />
    <!--
      The following RDF block contains metadata that applies to this document
      as a whole, as indicated by the empty about attribute on the 
      <rdf:Description> element.
    -->
    <rdf:Description rdf:about="">
      <!--
        The Model Builder Metadata.  The Dublin Core "creator" element is used  
        to indicate the person who translated the model into CellML.
      -->
      <dc:creator rdf:parseType="Resource">
        <vCard:N rdf:parseType="Resource">
          <vCard:Family>Lawson</vCard:Family>
          <vCard:Given>James</vCard:Given>
		</vCard:N>
        <vCard:EMAIL rdf:parseType="Resource">
          <rdf:value>j.lawson@auckland@auckland.ac.nz</rdf:value>
          <rdf:type rdf:resource="http://imc.org/vCard/3.0#internet"/>
        </vCard:EMAIL>
        <vCard:ORG rdf:parseType="Resource">
          <vCard:Orgname>The University of Auckland</vCard:Orgname>
          <vCard:Orgunit>Auckland Bioengineering Institute</vCard:Orgunit>
        </vCard:ORG>
      </dc:creator>
      
      <!--
        The Creation Date metadata. This is the date on which the model
        was translated into CellML.
      -->
      <dcterms:created rdf:parseType="Resource">
        <dcterms:W3CDTF>2009-12-02</dcterms:W3CDTF>
      </dcterms:created>
      
      <!--
        The Modification History metadata. This lists the changes that have been
        made to the document, who made the changes, and when they were made.
      -->
     
    </rdf:Description>
    
    <!--
      The following metadata refers to the model itself, as indicated by the
      reference to the ID "lockwood_2006", which is    
      declared on the <model> element.
    -->
    <rdf:Description rdf:about="#lockwood_2006">
   
      <!--  Keyword(s) -->
      <bqs:reference rdf:parseType="Resource">
        <dc:subject rdf:parseType="Resource">
          <bqs:subject_type>keyword</bqs:subject_type>
          <rdf:value>
            <rdf:Bag>
              <rdf:li>PKPD</rdf:li>
			  <rdf:li>pharmacokinetic pharmacodynamic model</rdf:li>
            </rdf:Bag>
          </rdf:value>
        </dc:subject>
      </bqs:reference>

      <bqs:reference rdf:parseType="Resource">
        <bqs:Pubmed_id>16906456</bqs:Pubmed_id>
        <bqs:JournalArticle rdf:parseType="Resource">
          <dc:creator>
            <rdf:Seq>
              <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N rdf:parseType="Resource">
                    <vCard:Family>Lockwood</vCard:Family>
                    <vCard:Given>P</vCard:Given>
				  </vCard:N>
                </bqs:Person>
              </rdf:li>
			  <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N rdf:parseType="Resource">
                    <vCard:Family>Ewy</vCard:Family>
                    <vCard:Given>W</vCard:Given>
				  </vCard:N>
                </bqs:Person>
              </rdf:li>
              <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N rdf:parseType="Resource">
                    <vCard:Family>Hermann</vCard:Family>
                    <vCard:Given>D</vCard:Given>
				  </vCard:N>
                </bqs:Person>
              </rdf:li>
			  <rdf:li rdf:parseType="Resource">
                <bqs:Person rdf:parseType="Resource">
                  <vCard:N rdf:parseType="Resource">
                    <vCard:Family>Holford</vCard:Family>
                    <vCard:Given>N</vCard:Given>
				  </vCard:N>
                </bqs:Person>
              </rdf:li>
            </rdf:Seq>
          </dc:creator>
          <dc:title>
            Application of clinical trial simulation to compare proof-of-concept study designs for drugs with a slow onset of effect; an example in Alzheimer's disease
          </dc:title>
          <dcterms:issued rdf:parseType="Resource">
            <dcterms:W3CDTF>2006-09</dcterms:W3CDTF>
          </dcterms:issued>
          <bqs:Journal rdf:parseType="Resource">
            <dc:title>Pharmaceutical Research</dc:title>
          </bqs:Journal>
		  <bqs:volume>23</bqs:volume>
          <bqs:first_page>2050</bqs:first_page>
          <bqs:last_page>2059</bqs:last_page>
        </bqs:JournalArticle>
      </bqs:reference>
    </rdf:Description>
  </rdf:RDF>
</model>