<?xml version="1.0"?>
<model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" xml:base="file:///home/jlaw060/cellml_repository/pmr2/overgaard_2007/overgaard_pharmacodynamic_2007.cellml" cmeta:id="overgaard_pharmacodynamic_2007" name="overgaard_pharmacodynamic_2007">
<documentation xmlns="http://cellml.org/tmp-documentation">
<article>
<articleinfo>
<author>
<firstname>Catherine</firstname>
<surname>Lloyd</surname>
<affiliation>
<shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
</affiliation>
</author>
</articleinfo>
<section id="sec_status">
<title>Model Status</title>
<para>
This CellML model represents the "Pharmacodynamic Model" in the published paper. The model runs in PCEnv to recreate figure3 in the published paper. The units have been chaecked and they balance. Note that this model will not run in COR due to the presence of a "remainder" function. Please also note that we have chosen "high dose levels" (deltahigh_dose=1), and IL-21 is administered on days 1, 3 and 5, at a dose of 3mg/kg. The f2 function in equation 3 of the paper has been ignored as this is only relevant for drugs directly affecting heat conduction, which IL-21 is not expected to. Thank you to the original model author Rune Viig Overgaard for providing us with the NONMEM code and for clarifying certain aspects of the model.
</para>
</section>
<sect1 id="sec_structure">
<title>Model Structure</title>
<para>
ABSTRACT: PURPOSE: To describe the pharmacodynamic effects of recombinant human interleukin-21 (IL-21) on core body temperature in cynomolgus monkeys using basic mechanisms of heat regulation. A major effort was devoted to compare the use of ordinary differential equations (ODEs) with stochastic differential equations (SDEs) in pharmacokinetic pharmacodynamic (PKPD) modelling. METHODS: A temperature model was formulated including circadian rhythm, metabolism, heat loss, and a thermoregulatory set-point. This model was formulated as a mixed-effects model based on SDEs using NONMEM. RESULTS: The effects of IL-21 were on the set-point and the circadian rhythm of metabolism. The model was able to describe a complex set of IL-21 induced phenomena, including 1) disappearance of the circadian rhythm, 2) no effect after first dose, and 3) high variability after second dose. SDEs provided a more realistic description with improved simulation properties, and further changed the model into one that could not be falsified by the autocorrelation function. CONCLUSIONS: The IL-21 induced effects on thermoregulation in cynomolgus monkeys are explained by a biologically plausible model. The quality of the model was improved by the use of SDEs.
</para>
<informalfigure float="0" id="fig_reaction_diagram">
<mediaobject>
<imageobject>
<objectinfo>
<title>model diagram</title>
</objectinfo>
<imagedata fileref="overgaard_2007.png"/>
</imageobject>
</mediaobject>
<caption>Schematic diagram of the model for IL-21 induced regulation of core body temperature.</caption>
</informalfigure>
<para>
The original paper reference is cited below:
</para>
<para>
PKPD model of interleukin-21 effects on thermoregulation in monkeys - application and evaluation of stochastic differential equations, Rune Viig Overgaard, Nick Holford, Klaus A. Rytved and Henrik Madsen, 2007, Pharmaceutical Research, 24, (2), 298-309. <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&cmd=Retrieve&dopt=AbstractPlus&list_uids=17009101&query_hl=1&itool=pubmed_docsum">PubMed ID: 17009101</ulink>
</para>
</sect1>
</article>
</documentation>
<units name="hour">
<unit multiplier="3600" units="second"/>
</units>
<units name="day">
<unit multiplier="86400" units="second"/>
</units>
<units name="per_hour">
<unit exponent="-1" units="hour"/>
</units>
<units name="per_day">
<unit exponent="-1" units="day"/>
</units>
<units name="hours_per_day">
<unit units="hour"/>
<unit exponent="-1" units="day"/>
</units>
<units name="W">
<unit units="watt"/>
</units>
<units name="kg_per_day_mg">
<unit units="kilogram"/>
<unit exponent="-1" units="day"/>
<unit exponent="-1" prefix="milli" units="gram"/>
</units>
<units name="mg_per_kg">
<unit prefix="milli" units="gram"/>
<unit exponent="-1" units="kilogram"/>
</units>
<units name="W_per_kg">
<unit units="W"/>
<unit exponent="-1" units="kilogram"/>
</units>
<units name="W_per_kg_C">
<unit units="W"/>
<unit exponent="-1" units="kilogram"/>
<unit exponent="-1" units="celsius"/>
</units>
<units name="W_per_kg_C2">
<unit units="W"/>
<unit exponent="-1" units="kilogram"/>
<unit exponent="-2" units="celsius"/>
</units>
<units name="kJ_per_kg_C">
<unit prefix="kilo" units="joule"/>
<unit exponent="-1" units="kilogram"/>
<unit exponent="-1" units="celsius"/>
</units>
<units name="seconds_per_hour">
<unit units="second"/>
<unit exponent="-1" units="hour"/>
</units>
<component name="environment">
<variable cmeta:id="environment_time" name="time" public_interface="out" units="hour"/>
</component>
<component name="model_parameters">
<variable initial_value="21.0" name="T_a" public_interface="out" units="celsius"/>
<variable initial_value="38.0" name="T_b" public_interface="out" units="celsius"/>
<variable initial_value="1.57" name="delta_T" public_interface="out" units="celsius"/>
<variable initial_value="0.0258" name="kinc" public_interface="out" units="W_per_kg_C2"/>
<variable initial_value="24.0" name="tdose1" public_interface="out" units="hour"/>
<variable initial_value="72.0" name="tdose2" public_interface="out" units="hour"/>
<variable initial_value="120.0" name="tdose3" public_interface="out" units="hour"/>
</component>
<component name="M_c">
<variable name="M_c" public_interface="out" units="W_per_kg"/>
<variable initial_value="17.5" name="t_day" units="hour"/>
<variable initial_value="6.73" name="t_night" units="hour"/>
<variable name="tprime" units="second"/>
<variable initial_value="86400" name="day_length" units="second"/>
<variable name="M_day" public_interface="in" units="W_per_kg"/>
<variable name="M_night" public_interface="in" units="W_per_kg"/>
<variable name="time" public_interface="in" units="hour"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="1">
<apply>
<eq/>
<ci>tprime</ci>
<apply>
<rem/>
<apply>
<times/>
<ci>time</ci>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_hour">3600</cn>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="second">1</cn>
</apply>
<ci>day_length</ci>
</apply>
</apply>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="1a">
<apply>
<eq/>
<ci>M_c</ci>
<piecewise>
<piece>
<ci>M_night</ci>
<apply>
<and/>
<apply>
<geq/>
<apply>
<divide/>
<ci>tprime</ci>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="seconds_per_hour">3600</cn>
</apply>
<ci>t_night</ci>
</apply>
<apply>
<lt/>
<apply>
<divide/>
<ci>tprime</ci>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="seconds_per_hour">3600</cn>
</apply>
<ci>t_day</ci>
</apply>
</apply>
</piece>
<otherwise>
<ci>M_day</ci>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="M">
<variable cmeta:id="M_M" initial_value="3.5" name="M" public_interface="out" units="W_per_kg"/>
<variable initial_value="1.1375" name="km" units="per_hour"/>
<variable name="M_c" public_interface="in" units="W_per_kg"/>
<variable name="time" public_interface="in" units="hour"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="2">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>M</ci>
</apply>
<apply>
<times/>
<apply>
<minus/>
<ci>km</ci>
</apply>
<apply>
<minus/>
<ci>M</ci>
<ci>M_c</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="T">
<variable cmeta:id="T_T" initial_value="38.785" name="T" public_interface="out" units="celsius"/>
<variable initial_value="3.47" name="c" units="kJ_per_kg_C"/>
<variable name="k" public_interface="in" units="W_per_kg_C"/>
<variable name="M" public_interface="in" units="W_per_kg"/>
<variable name="T_a" public_interface="in" units="celsius"/>
<variable name="time" public_interface="in" units="hour"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="3">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>T</ci>
</apply>
<apply>
<times/>
<apply>
<power/>
<ci>c</ci>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">-1</cn>
</apply>
<apply>
<minus/>
<ci>M</ci>
<apply>
<times/>
<ci>k</ci>
<apply>
<minus/>
<ci>T</ci>
<ci>T_a</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="k">
<variable name="k" public_interface="out" units="W_per_kg_C"/>
<variable initial_value="0.0" name="BR" units="dimensionless"/>
<variable initial_value="0.144" name="pEtot" units="dimensionless"/>
<variable initial_value="5.35" name="kR" units="per_day"/>
<variable initial_value="3.0" name="AMT_dose" units="mg_per_kg"/>
<variable initial_value="1.0" name="pEf1" units="per_day"/>
<variable initial_value="0.2" name="pEs1" units="kg_per_day_mg"/>
<variable initial_value="3.57" name="pEf2" units="per_day"/>
<variable initial_value="2.43" name="pEs2" units="kg_per_day_mg"/>
<variable initial_value="8.0" name="pEf3" units="per_day"/>
<variable initial_value="50.0" name="pEs3" units="kg_per_day_mg"/>
<variable cmeta:id="k_E_slow" name="E_slow" units="per_day"/>
<variable cmeta:id="k_E_fast" name="E_fast" units="per_day"/>
<variable name="f2_drug" units="W_per_kg_C"/>
<variable name="kinc" public_interface="in" units="W_per_kg_C2"/>
<variable name="kb" public_interface="in" units="W_per_kg_C"/>
<variable name="T" public_interface="in" units="celsius"/>
<variable name="T_b" public_interface="in" units="celsius"/>
<variable name="f_prime" public_interface="in" units="dimensionless"/>
<variable name="gNsTs1" public_interface="in" units="dimensionless"/>
<variable name="gNsTs2" public_interface="in" units="dimensionless"/>
<variable name="gNsTs3" public_interface="in" units="dimensionless"/>
<variable name="gNfTf1" public_interface="in" units="dimensionless"/>
<variable name="gNfTf2" public_interface="in" units="dimensionless"/>
<variable name="gNfTf3" public_interface="in" units="dimensionless"/>
<variable name="time" public_interface="in" units="hour"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="4">
<apply>
<eq/>
<ci>k</ci>
<apply>
<plus/>
<ci>kb</ci>
<apply>
<times/>
<ci>kinc</ci>
<apply>
<minus/>
<ci>T</ci>
<apply>
<times/>
<ci>T_b</ci>
<apply>
<plus/>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
<apply>
<times/>
<ci>pEtot</ci>
<ci>BR</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>f2_drug</ci>
</apply>
</apply>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="4a">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>BR</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<apply>
<times/>
<ci>f_prime</ci>
<apply>
<plus/>
<ci>E_slow</ci>
<ci>E_fast</ci>
</apply>
</apply>
<apply>
<minus/>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
<ci>BR</ci>
</apply>
</apply>
<apply>
<times/>
<ci>kR</ci>
<ci>BR</ci>
</apply>
</apply>
</apply>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="4b">
<apply>
<eq/>
<ci>E_slow</ci>
<apply>
<times/>
<ci>AMT_dose</ci>
<ci>pEs2</ci>
<apply>
<plus/>
<ci>gNsTs1</ci>
<ci>gNsTs2</ci>
<ci>gNsTs3</ci>
</apply>
</apply>
</apply>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="4c">
<apply>
<eq/>
<ci>E_fast</ci>
<apply>
<times/>
<ci>pEf2</ci>
<apply>
<plus/>
<ci>gNfTf1</ci>
<ci>gNfTf2</ci>
<ci>gNfTf3</ci>
</apply>
</apply>
</apply>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="4d">
<apply>
<eq/>
<ci>f2_drug</ci>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="W_per_kg_C">0</cn>
</apply>
</math>
</component>
<component name="T_day">
<variable name="T_day" public_interface="out" units="celsius"/>
<variable name="T_b" public_interface="in" units="celsius"/>
<variable name="delta_T" public_interface="in" units="celsius"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="5">
<apply>
<eq/>
<ci>T_day</ci>
<apply>
<plus/>
<ci>T_b</ci>
<apply>
<divide/>
<ci>delta_T</ci>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</math>
</component>
<component name="T_night">
<variable name="T_night" public_interface="out" units="celsius"/>
<variable name="T_b" public_interface="in" units="celsius"/>
<variable name="delta_T" public_interface="in" units="celsius"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="6">
<apply>
<eq/>
<ci>T_night</ci>
<apply>
<minus/>
<ci>T_b</ci>
<apply>
<divide/>
<ci>delta_T</ci>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</math>
</component>
<component name="kb">
<variable name="kb" public_interface="out" units="W_per_kg_C"/>
<variable initial_value="3.0" name="M_b" units="W_per_kg"/>
<variable name="T_b" public_interface="in" units="celsius"/>
<variable name="T_a" public_interface="in" units="celsius"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="7">
<apply>
<eq/>
<ci>kb</ci>
<apply>
<divide/>
<ci>M_b</ci>
<apply>
<minus/>
<ci>T_b</ci>
<ci>T_a</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="M_day">
<variable name="M_day" public_interface="out" units="W_per_kg"/>
<variable name="kinc" public_interface="in" units="W_per_kg_C2"/>
<variable name="kb" public_interface="in" units="W_per_kg_C"/>
<variable name="T_b" public_interface="in" units="celsius"/>
<variable name="T_a" public_interface="in" units="celsius"/>
<variable name="T_day" public_interface="in" units="celsius"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="8">
<apply>
<eq/>
<ci>M_day</ci>
<apply>
<times/>
<apply>
<plus/>
<ci>kb</ci>
<apply>
<times/>
<ci>kinc</ci>
<apply>
<minus/>
<ci>T_day</ci>
<ci>T_b</ci>
</apply>
</apply>
</apply>
<apply>
<minus/>
<ci>T_day</ci>
<ci>T_a</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="M_night">
<variable name="M_night" public_interface="out" units="W_per_kg"/>
<variable cmeta:id="M_night_f_prime" name="f_prime" public_interface="out" units="dimensionless"/>
<variable initial_value="45.12" name="t_prime" units="hour"/>
<variable initial_value="0.2229166" name="alpha" units="per_hour"/>
<variable initial_value="1.0" name="delta_high_dose" units="dimensionless"/>
<variable name="M_night_baseline" units="W_per_kg"/>
<variable name="tdose1" public_interface="in" units="hour"/>
<variable name="kinc" public_interface="in" units="W_per_kg_C2"/>
<variable name="kb" public_interface="in" units="W_per_kg_C"/>
<variable name="T_b" public_interface="in" units="celsius"/>
<variable name="T_a" public_interface="in" units="celsius"/>
<variable name="T_night" public_interface="in" units="celsius"/>
<variable name="M_day" public_interface="in" units="W_per_kg"/>
<variable name="time" public_interface="in" units="hour"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="9">
<apply>
<eq/>
<ci>f_prime</ci>
<apply>
<times/>
<ci>delta_high_dose</ci>
<apply>
<power/>
<apply>
<plus/>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci>alpha</ci>
</apply>
<apply>
<minus/>
<ci>time</ci>
<apply>
<plus/>
<ci>tdose1</ci>
<ci>t_prime</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">-1</cn>
</apply>
</apply>
</apply>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="9a">
<apply>
<eq/>
<ci>M_night</ci>
<apply>
<plus/>
<apply>
<times/>
<apply>
<minus/>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
<ci>f_prime</ci>
</apply>
<ci>M_night_baseline</ci>
</apply>
<apply>
<times/>
<ci>f_prime</ci>
<ci>M_day</ci>
</apply>
</apply>
</apply>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="9b">
<apply>
<eq/>
<ci>M_night_baseline</ci>
<apply>
<times/>
<apply>
<plus/>
<ci>kb</ci>
<apply>
<times/>
<ci>kinc</ci>
<apply>
<minus/>
<ci>T_night</ci>
<ci>T_b</ci>
</apply>
</apply>
</apply>
<apply>
<minus/>
<ci>T_night</ci>
<ci>T_a</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="gNT">
<variable name="gNsTs1" public_interface="out" units="dimensionless"/>
<variable name="gNsTs2" public_interface="out" units="dimensionless"/>
<variable name="gNsTs3" public_interface="out" units="dimensionless"/>
<variable name="gNfTf1" public_interface="out" units="dimensionless"/>
<variable name="gNfTf2" public_interface="out" units="dimensionless"/>
<variable name="gNfTf3" public_interface="out" units="dimensionless"/>
<variable initial_value="4.0" name="Ns" units="dimensionless"/>
<variable initial_value="4.0" name="Nf" units="dimensionless"/>
<variable initial_value="2.45" name="Ts" units="day"/>
<variable initial_value="0.368" name="Tf" units="day"/>
<variable name="X1" units="day"/>
<variable name="X2" units="day"/>
<variable name="X3" units="day"/>
<variable name="Kf" units="per_day"/>
<variable name="Ks" units="per_day"/>
<variable name="tdose1" public_interface="in" units="hour"/>
<variable name="tdose2" public_interface="in" units="hour"/>
<variable name="tdose3" public_interface="in" units="hour"/>
<variable name="time" public_interface="in" units="hour"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="10">
<apply>
<eq/>
<ci>gNsTs1</ci>
<piecewise>
<piece>
<apply>
<times/>
<apply>
<divide/>
<apply>
<power/>
<ci>Ks</ci>
<ci>Ns</ci>
</apply>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_day">6</cn>
</apply>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci>Ks</ci>
</apply>
<ci>X1</ci>
</apply>
</apply>
<apply>
<power/>
<ci>X1</ci>
<apply>
<minus/>
<ci>Ns</ci>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<gt/>
<ci>X1</ci>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="day">0</cn>
</apply>
</piece>
<otherwise>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0</cn>
</otherwise>
</piecewise>
</apply>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="11">
<apply>
<eq/>
<ci>gNsTs2</ci>
<piecewise>
<piece>
<apply>
<times/>
<apply>
<divide/>
<apply>
<power/>
<ci>Ks</ci>
<ci>Ns</ci>
</apply>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_day">6</cn>
</apply>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci>Ks</ci>
</apply>
<ci>X2</ci>
</apply>
</apply>
<apply>
<power/>
<ci>X2</ci>
<apply>
<minus/>
<ci>Ns</ci>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<gt/>
<ci>X2</ci>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="day">0</cn>
</apply>
</piece>
<otherwise>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0</cn>
</otherwise>
</piecewise>
</apply>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="12">
<apply>
<eq/>
<ci>gNsTs3</ci>
<piecewise>
<piece>
<apply>
<times/>
<apply>
<divide/>
<apply>
<power/>
<ci>Ks</ci>
<ci>Ns</ci>
</apply>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_day">6</cn>
</apply>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci>Ks</ci>
</apply>
<ci>X3</ci>
</apply>
</apply>
<apply>
<power/>
<ci>X3</ci>
<apply>
<minus/>
<ci>Ns</ci>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<gt/>
<ci>X3</ci>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="day">0</cn>
</apply>
</piece>
<otherwise>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0</cn>
</otherwise>
</piecewise>
</apply>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="13">
<apply>
<eq/>
<ci>gNfTf1</ci>
<piecewise>
<piece>
<apply>
<times/>
<apply>
<divide/>
<apply>
<power/>
<ci>Kf</ci>
<ci>Nf</ci>
</apply>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_day">6</cn>
</apply>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci>Kf</ci>
</apply>
<ci>X1</ci>
</apply>
</apply>
<apply>
<power/>
<ci>X1</ci>
<apply>
<minus/>
<ci>Nf</ci>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<gt/>
<ci>X1</ci>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="day">0</cn>
</apply>
</piece>
<otherwise>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0</cn>
</otherwise>
</piecewise>
</apply>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="14">
<apply>
<eq/>
<ci>gNfTf2</ci>
<piecewise>
<piece>
<apply>
<times/>
<apply>
<divide/>
<apply>
<power/>
<ci>Kf</ci>
<ci>Nf</ci>
</apply>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_day">6</cn>
</apply>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci>Kf</ci>
</apply>
<ci>X2</ci>
</apply>
</apply>
<apply>
<power/>
<ci>X2</ci>
<apply>
<minus/>
<ci>Nf</ci>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<gt/>
<ci>X2</ci>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="day">0</cn>
</apply>
</piece>
<otherwise>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0</cn>
</otherwise>
</piecewise>
</apply>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="15">
<apply>
<eq/>
<ci>gNfTf3</ci>
<piecewise>
<piece>
<apply>
<times/>
<apply>
<divide/>
<apply>
<power/>
<ci>Kf</ci>
<ci>Nf</ci>
</apply>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="per_day">6</cn>
</apply>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci>Kf</ci>
</apply>
<ci>X3</ci>
</apply>
</apply>
<apply>
<power/>
<ci>X3</ci>
<apply>
<minus/>
<ci>Nf</ci>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<gt/>
<ci>X3</ci>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="day">0</cn>
</apply>
</piece>
<otherwise>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">0</cn>
</otherwise>
</piecewise>
</apply>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="16">
<apply>
<eq/>
<ci>X1</ci>
<apply>
<divide/>
<apply>
<minus/>
<ci>time</ci>
<ci>tdose1</ci>
</apply>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="hours_per_day">24</cn>
</apply>
</apply>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="17">
<apply>
<eq/>
<ci>X2</ci>
<apply>
<divide/>
<apply>
<minus/>
<ci>time</ci>
<ci>tdose2</ci>
</apply>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="hours_per_day">24</cn>
</apply>
</apply>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="18">
<apply>
<eq/>
<ci>X3</ci>
<apply>
<divide/>
<apply>
<minus/>
<ci>time</ci>
<ci>tdose3</ci>
</apply>
<cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="hours_per_day">24</cn>
</apply>
</apply>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="19">
<apply>
<eq/>
<ci>Kf</ci>
<apply>
<divide/>
<ci>Nf</ci>
<ci>Tf</ci>
</apply>
</apply>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="20">
<apply>
<eq/>
<ci>Ks</ci>
<apply>
<divide/>
<ci>Ns</ci>
<ci>Ts</ci>
</apply>
</apply>
</math>
</component>
<connection>
<map_components component_1="M_c" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="M" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="T" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="M_night" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="k" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="gNT" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="M_day" component_2="M_c"/>
<map_variables variable_1="M_day" variable_2="M_day"/>
</connection>
<connection>
<map_components component_1="M_night" component_2="M_c"/>
<map_variables variable_1="M_night" variable_2="M_night"/>
</connection>
<connection>
<map_components component_1="M_c" component_2="M"/>
<map_variables variable_1="M_c" variable_2="M_c"/>
</connection>
<connection>
<map_components component_1="M" component_2="T"/>
<map_variables variable_1="M" variable_2="M"/>
</connection>
<connection>
<map_components component_1="model_parameters" component_2="T"/>
<map_variables variable_1="T_a" variable_2="T_a"/>
</connection>
<connection>
<map_components component_1="k" component_2="T"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="k" variable_2="k"/>
</connection>
<connection>
<map_components component_1="kb" component_2="k"/>
<map_variables variable_1="kb" variable_2="kb"/>
</connection>
<connection>
<map_components component_1="M_night" component_2="k"/>
<map_variables variable_1="f_prime" variable_2="f_prime"/>
</connection>
<connection>
<map_components component_1="gNT" component_2="k"/>
<map_variables variable_1="gNsTs1" variable_2="gNsTs1"/>
<map_variables variable_1="gNsTs2" variable_2="gNsTs2"/>
<map_variables variable_1="gNsTs3" variable_2="gNsTs3"/>
<map_variables variable_1="gNfTf1" variable_2="gNfTf1"/>
<map_variables variable_1="gNfTf2" variable_2="gNfTf2"/>
<map_variables variable_1="gNfTf3" variable_2="gNfTf3"/>
</connection>
<connection>
<map_components component_1="model_parameters" component_2="k"/>
<map_variables variable_1="T_b" variable_2="T_b"/>
<map_variables variable_1="kinc" variable_2="kinc"/>
</connection>
<connection>
<map_components component_1="model_parameters" component_2="T_day"/>
<map_variables variable_1="T_b" variable_2="T_b"/>
<map_variables variable_1="delta_T" variable_2="delta_T"/>
</connection>
<connection>
<map_components component_1="model_parameters" component_2="T_night"/>
<map_variables variable_1="T_b" variable_2="T_b"/>
<map_variables variable_1="delta_T" variable_2="delta_T"/>
</connection>
<connection>
<map_components component_1="model_parameters" component_2="kb"/>
<map_variables variable_1="T_b" variable_2="T_b"/>
<map_variables variable_1="T_a" variable_2="T_a"/>
</connection>
<connection>
<map_components component_1="kb" component_2="M_day"/>
<map_variables variable_1="kb" variable_2="kb"/>
</connection>
<connection>
<map_components component_1="T_day" component_2="M_day"/>
<map_variables variable_1="T_day" variable_2="T_day"/>
</connection>
<connection>
<map_components component_1="model_parameters" component_2="M_day"/>
<map_variables variable_1="T_a" variable_2="T_a"/>
<map_variables variable_1="T_b" variable_2="T_b"/>
<map_variables variable_1="kinc" variable_2="kinc"/>
</connection>
<connection>
<map_components component_1="kb" component_2="M_night"/>
<map_variables variable_1="kb" variable_2="kb"/>
</connection>
<connection>
<map_components component_1="M_day" component_2="M_night"/>
<map_variables variable_1="M_day" variable_2="M_day"/>
</connection>
<connection>
<map_components component_1="T_night" component_2="M_night"/>
<map_variables variable_1="T_night" variable_2="T_night"/>
</connection>
<connection>
<map_components component_1="model_parameters" component_2="M_night"/>
<map_variables variable_1="T_a" variable_2="T_a"/>
<map_variables variable_1="T_b" variable_2="T_b"/>
<map_variables variable_1="kinc" variable_2="kinc"/>
<map_variables variable_1="tdose1" variable_2="tdose1"/>
</connection>
<connection>
<map_components component_1="model_parameters" component_2="gNT"/>
<map_variables variable_1="tdose1" variable_2="tdose1"/>
<map_variables variable_1="tdose2" variable_2="tdose2"/>
<map_variables variable_1="tdose3" variable_2="tdose3"/>
</connection>
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"><rdf:Description rdf:nodeID="n1"><Given xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Klaus</Given><Family xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Rytved</Family><Other xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">A</Other></rdf:Description><rdf:Description rdf:about=""><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">PKPD model of interleukin-21 effects on thermoregulation in monkeys (Pharmacodynamic)</title><creator xmlns="http://purl.org/dc/elements/1.1/"><rdf:Description rdf:nodeID="n2"/></creator><created xmlns="http://purl.org/dc/terms/"><rdf:Description rdf:nodeID="n3"/></created></rdf:Description><rdf:Description rdf:nodeID="n4"><Given xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Catherine</Given><Family xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Lloyd</Family><Other xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">May</Other></rdf:Description><rdf:Description rdf:nodeID="n5"><rdf:first><rdf:Description rdf:nodeID="n6"/></rdf:first><rdf:rest><rdf:Description rdf:about="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/></rdf:rest></rdf:Description><rdf:Description rdf:about="#overgaard_pharmacodynamic_2007"><reference xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:nodeID="n7"/></reference><reference xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:nodeID="n8"/></reference><cmeta:comment><rdf:Description rdf:nodeID="n9"/></cmeta:comment><cmeta:species xml:lang="en">Monkey</cmeta:species><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">
The Overgaard et al. 2007 PKPD model of the effects of IL_21 on thermoregulation in monkeys
</title><simulation xmlns="http://www.cellml.org/metadata/simulation/1.0#"><rdf:Description rdf:nodeID="n10"/></simulation></rdf:Description><rdf:Description rdf:nodeID="n11"><volume xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">24</volume><creator xmlns="http://purl.org/dc/elements/1.1/"><rdf:Description rdf:nodeID="n12"/></creator><Journal xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:nodeID="n13"/></Journal><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">
PKPD model of interleukin-21 effects on thermoregulation in monkeys--application and evaluation of stochastic differential equations
</title><first_page xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">298</first_page><issued xmlns="http://purl.org/dc/terms/"><rdf:Description rdf:nodeID="n14"/></issued><last_page xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">309</last_page></rdf:Description><rdf:Description rdf:nodeID="n15"><N xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:nodeID="n1"/></N></rdf:Description><rdf:Description rdf:nodeID="n16"><Orgname xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">The University of Auckland</Orgname><Orgunit xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Auckland Bioengineering Institute</Orgunit></rdf:Description><rdf:Description rdf:nodeID="n17"><Person xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:nodeID="n15"/></Person></rdf:Description><rdf:Description rdf:nodeID="n18"><rdf:type><rdf:Description rdf:about="http://www.example.org/foo/bar/bazfoo"/></rdf:type></rdf:Description><rdf:Description rdf:nodeID="n19"><Person xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:nodeID="n20"/></Person></rdf:Description><rdf:Description rdf:nodeID="n21"><rdf:type><rdf:Description rdf:about="http://imc.org/vCard/3.0#internet"/></rdf:type><rdf:value xml:lang="en">c.lloyd@auckland@auckland.ac.nz</rdf:value></rdf:Description><rdf:Description rdf:nodeID="n3"><W3CDTF xmlns="http://purl.org/dc/terms/" xml:lang="en">2009-10-20</W3CDTF></rdf:Description><rdf:Description rdf:nodeID="n8"><Pubmed_id xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">17009101</Pubmed_id><JournalArticle xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:nodeID="n11"/></JournalArticle></rdf:Description><rdf:Description rdf:nodeID="n22"><N xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:nodeID="n23"/></N></rdf:Description><rdf:Description rdf:nodeID="n24"><Person xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:nodeID="n22"/></Person></rdf:Description><rdf:Description rdf:nodeID="n23"><Given xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Henrik</Given><Family xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Madsen</Family></rdf:Description><rdf:Description rdf:nodeID="n9"><creator xmlns="http://purl.org/dc/elements/1.1/"><rdf:Description rdf:nodeID="n25"/></creator><rdf:value xml:lang="en">
This is the CellML description of Overgaard et al.'s PKPD model of the effects of IL_21 on thermoregulation in monkeys
</rdf:value></rdf:Description><rdf:Description rdf:nodeID="n12"><rdf:_4><rdf:Description rdf:nodeID="n24"/></rdf:_4><rdf:type><rdf:Description rdf:about="http://www.w3.org/1999/02/22-rdf-syntax-ns#Seq"/></rdf:type><rdf:_1><rdf:Description rdf:nodeID="n19"/></rdf:_1><rdf:_3><rdf:Description rdf:nodeID="n17"/></rdf:_3><rdf:_2><rdf:Description rdf:nodeID="n26"/></rdf:_2></rdf:Description><rdf:Description rdf:nodeID="n13"><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">Pharmaceutical Research</title></rdf:Description><rdf:Description rdf:nodeID="n27"><Given xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Nick</Given><Family xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Holford</Family></rdf:Description><rdf:Description rdf:nodeID="n14"><W3CDTF xmlns="http://purl.org/dc/terms/" xml:lang="en">2007-02</W3CDTF></rdf:Description><rdf:Description rdf:nodeID="n2"><EMAIL xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:nodeID="n21"/></EMAIL><N xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:nodeID="n4"/></N><ORG xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:nodeID="n16"/></ORG></rdf:Description><rdf:Description rdf:nodeID="n28"><N xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:nodeID="n27"/></N></rdf:Description><rdf:Description rdf:nodeID="n25"><FN xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Catherine Lloyd</FN></rdf:Description><rdf:Description rdf:nodeID="n6"><endingValue xmlns="http://www.cellml.org/metadata/simulation/1.0#" xml:lang="en">200</endingValue><pointDensity xmlns="http://www.cellml.org/metadata/simulation/1.0#nonstandard-" xml:lang="en">10000</pointDensity><maximumStepSize xmlns="http://www.cellml.org/metadata/simulation/1.0#" xml:lang="en">0.1</maximumStepSize></rdf:Description><rdf:Description rdf:nodeID="n7"><subject xmlns="http://purl.org/dc/elements/1.1/"><rdf:Description rdf:nodeID="n29"/></subject></rdf:Description><rdf:Description rdf:nodeID="n30"><Given xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Rune</Given><Family xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Overgaard</Family><Other xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Viig</Other></rdf:Description><rdf:Description rdf:nodeID="n20"><N xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:nodeID="n30"/></N></rdf:Description><rdf:Description rdf:nodeID="n29"><subject_type xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">keyword</subject_type><rdf:value><rdf:Description rdf:nodeID="n31"/></rdf:value></rdf:Description><rdf:Description rdf:nodeID="n10"><boundIntervals xmlns="http://www.cellml.org/metadata/simulation/1.0#"><rdf:Description rdf:nodeID="n5"/></boundIntervals></rdf:Description><rdf:Description rdf:nodeID="n26"><Person xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:nodeID="n28"/></Person></rdf:Description><rdf:Description rdf:nodeID="n31"><rdf:_4 xml:lang="en">immunomodulation</rdf:_4><rdf:type><rdf:Description rdf:about="http://www.w3.org/1999/02/22-rdf-syntax-ns#Bag"/></rdf:type><rdf:_1 xml:lang="en">PKPD</rdf:_1><rdf:_3 xml:lang="en">thermoregulation</rdf:_3><rdf:_2 xml:lang="en">metabolism</rdf:_2></rdf:Description></rdf:RDF></model>