<?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="shrestha_2010" name="shrestha_2010">
<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 runs in both COR and OpenCel to replicate the published results (figures 7 and 10). The units have been checked and they are consistent. This particular version of the model represents subject 1.
</para>
</section>
<sect1 id="sec_structure">
<title>Model Structure</title>
<para>
ABSTRACT: A complex bio-mechanism, commonly referred to as calcium homeostasis, regulates plasma ionized calcium (Ca(2+)) concentration in the human body within a narrow range which is crucial for maintaining normal physiology and metabolism. Taking a step towards creating a complete mathematical model of calcium homeostasis, we focus on the short-term dynamics of calcium homeostasis and consider the response of the parathyroid glands to acute changes in plasma Ca(2+) concentration. We review available models, discuss their limitations, then present a two-pool, linear, time-varying model to describe the dynamics of this calcium homeostasis subsystem, the Ca-PTH axis. We propose that plasma PTH concentration and plasma Ca(2+) concentration bear an asymmetric reverse sigmoid relation. The parameters of our model are successfully estimated based on clinical data corresponding to three healthy subjects that have undergone induced hypocalcemic clamp tests. In the first validation of this kind, with parameters estimated separately for each subject we test the model's ability to predict the same subject's induced hypercalcemic clamp test responses. Our results demonstrate that a two-pool, linear, time-varying model with an asymmetric reverse sigmoid relation characterizes the short-term dynamics of the Ca-PTH axis.
</para>
<para>
The original paper reference is cited below:
</para>
<para>
A mathematical model of parathyroid hormone response to acute changes in plasma ionized calcium concentration in humans, Shrestha RP, Hollot CV, Chipkin SR, Schmitt CP, and Chait Y, 2010, <emphasis>Mathematical Biosciences</emphasis>, 226, 46-57. <ulink url="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&list_uids=20406649&dopt=Abstract">PubMed ID: 20406649</ulink>
</para>
<informalfigure float="0" id="fig_pathway_diagram">
<mediaobject>
<imageobject>
<objectinfo>
<title>model diagram</title>
</objectinfo>
<imagedata fileref="shrestha_2010.png"/>
</imageobject>
</mediaobject>
<caption>A schematic diagram of the model.</caption>
</informalfigure>
</sect1>
</article>
</documentation>
<units name="minute">
<unit units="second" multiplier="60.0"/>
</units>
<units name="picomole">
<unit units="mole" prefix="pico"/>
</units>
<units name="millimole">
<unit prefix="milli" units="mole"/>
</units>
<units name="millimolar">
<unit prefix="milli" units="mole"/>
<unit exponent="-1" units="litre"/>
</units>
<units name="flux">
<unit units="picomole"/>
<unit units="minute" exponent="-1"/>
</units>
<units name="per_minute">
<unit units="minute" exponent="-1"/>
</units>
<units name="litre_per_millimole">
<unit units="litre"/>
<unit units="millimole" exponent="-1"/>
</units>
<component name="environment">
<variable name="time" units="minute" public_interface="out"/>
</component>
<component cmeta:id="x1" name="x1">
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
<rdf:Description rdf:about="#x1">
<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">x1</dc:title>
<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">PTH in PTG pool</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable initial_value="0.00" name="x1" public_interface="out" units="picomole"/>
<variable name="k" public_interface="in" units="flux"/>
<variable name="lambda_Ca" public_interface="in" units="per_minute"/>
<variable name="lambda1" public_interface="in" units="per_minute"/>
<variable name="time" public_interface="in" units="minute"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="1">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> x1 </ci>
</apply>
<apply>
<minus/>
<ci> k </ci>
<apply>
<plus/>
<apply>
<times/>
<ci> lambda_Ca </ci>
<ci> x1 </ci>
</apply>
<apply>
<times/>
<ci> lambda1 </ci>
<ci> x1 </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component cmeta:id="lambda_Ca" name="lambda_Ca">
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
<rdf:Description rdf:about="#lambda_Ca">
<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">lambda_Ca</dc:title>
<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">secretion rate function of PTH in PTG pool</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable name="lambda_Ca" public_interface="out" units="per_minute"/>
<variable name="A" public_interface="in" units="per_minute"/>
<variable name="B" public_interface="in" units="per_minute"/>
<variable name="Ca" public_interface="in" units="millimolar"/>
<variable name="S" public_interface="in" units="millimolar"/>
<variable name="m_Ca" public_interface="in" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="2">
<apply>
<eq/>
<ci> lambda_Ca </ci>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<minus/>
<ci> A </ci>
<ci> B </ci>
</apply>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<power/>
<apply>
<divide/>
<ci> Ca </ci>
<ci> S </ci>
</apply>
<ci> m_Ca </ci>
</apply>
</apply>
</apply>
<ci> B </ci>
</apply>
</apply>
</math>
</component>
<component name="m_Ca">
<variable name="m_Ca" public_interface="out" units="dimensionless"/>
<variable name="m1" public_interface="in" units="dimensionless"/>
<variable name="m2" public_interface="in" units="dimensionless"/>
<variable name="beta" public_interface="in" units="litre_per_millimole"/>
<variable name="Ca" public_interface="in" units="millimolar"/>
<variable name="R" public_interface="in" units="millimolar"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="3">
<apply>
<eq/>
<ci> m_Ca </ci>
<apply>
<plus/>
<apply>
<divide/>
<ci> m1 </ci>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci> beta </ci>
</apply>
<apply>
<minus/>
<ci> R </ci>
<ci> Ca </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci> m2 </ci>
</apply>
</apply>
</math>
</component>
<component cmeta:id="x2" name="x2">
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
<rdf:Description rdf:about="#x2">
<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">x2</dc:title>
<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">PTH in plasma</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable initial_value="0.00" name="x2" public_interface="out" units="picomole"/>
<variable name="PTH" units="picomole"/>
<variable name="x1" public_interface="in" units="picomole"/>
<variable name="lambda_Ca" public_interface="in" units="per_minute"/>
<variable name="lambda2" public_interface="in" units="per_minute"/>
<variable name="time" public_interface="in" units="minute"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="4">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci> time </ci>
</bvar>
<ci> x2 </ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> lambda_Ca </ci>
<ci> x1 </ci>
</apply>
<apply>
<times/>
<ci> lambda2 </ci>
<ci> x2 </ci>
</apply>
</apply>
</apply>
</math>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="4a">
<apply>
<eq/>
<ci> PTH </ci>
<apply>
<divide/>
<ci> x2 </ci>
<cn cellml:units="dimensionless">2.75</cn>
</apply>
</apply>
</math>
</component>
<component cmeta:id="Ca" name="Ca">
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
<rdf:Description rdf:about="#Ca">
<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">Ca</dc:title>
<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">Calcium in plasma</dcterms:alternative>
</rdf:Description>
</rdf:RDF>
<variable name="Ca" public_interface="out" units="millimolar"/>
<variable name="Ca_0" public_interface="in" units="millimolar"/>
<variable name="Ca_1" public_interface="in" units="millimolar"/>
<variable name="alpha" public_interface="in" units="per_minute"/>
<variable name="t0" public_interface="in" units="minute"/>
<variable name="time" public_interface="in" units="minute"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="5">
<apply>
<eq/>
<ci> Ca </ci>
<piecewise>
<piece>
<ci> Ca_0 </ci>
<apply><lt/>
<ci> time </ci>
<ci> t0 </ci>
</apply>
</piece>
<otherwise>
<apply><minus/>
<ci> Ca_0 </ci>
<apply>
<times/>
<ci> Ca_1 </ci>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci> alpha </ci>
</apply>
<apply>
<minus/>
<ci> time </ci>
<ci> t0 </ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
</math>
</component>
<component name="k">
<variable name="k" public_interface="out" units="flux"/>
<variable name="lambda1" public_interface="in" units="per_minute"/>
<variable name="lambda2" public_interface="in" units="per_minute"/>
<variable name="x1_n" public_interface="in" units="picomole"/>
<variable name="x2_n" public_interface="in" units="picomole"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="6">
<apply>
<eq/>
<ci> k </ci>
<apply>
<plus/>
<apply>
<times/>
<ci> lambda2 </ci>
<ci> x2_n </ci>
</apply>
<apply>
<times/>
<ci> lambda1 </ci>
<ci> x1_n </ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="S">
<variable name="S" public_interface="out" units="millimolar"/>
<variable name="Ca_0" public_interface="in" units="millimolar"/>
<variable name="A" public_interface="in" units="per_minute"/>
<variable name="B" public_interface="in" units="per_minute"/>
<variable name="x1_n" public_interface="in" units="picomole"/>
<variable name="x2_n" public_interface="in" units="picomole"/>
<variable name="lambda2" public_interface="in" units="per_minute"/>
<variable name="m_Ca" public_interface="in" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="7">
<apply>
<eq/>
<ci> S </ci>
<apply>
<times/>
<ci> Ca_0 </ci>
<apply>
<power/>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci> x1_n </ci>
<ci> B </ci>
</apply>
<apply>
<times/>
<ci> lambda2 </ci>
<ci> x2_n </ci>
</apply>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci> x1_n </ci>
<ci> A </ci>
</apply>
<apply>
<times/>
<ci> lambda2 </ci>
<ci> x2_n </ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci> m_Ca </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="A">
<variable name="A" public_interface="out" units="per_minute"/>
<variable name="lambda1" public_interface="in" units="per_minute"/>
<variable name="lambda2" public_interface="in" units="per_minute"/>
<variable name="x2_max" public_interface="in" units="picomole"/>
<variable name="k" public_interface="in" units="flux"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="8">
<apply>
<eq/>
<ci> A </ci>
<apply>
<divide/>
<apply>
<times/>
<ci> lambda1 </ci>
<ci> lambda2 </ci>
<ci> x2_max </ci>
</apply>
<apply>
<minus/>
<ci> k </ci>
<apply>
<times/>
<ci> lambda2 </ci>
<ci> x2_max </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="B">
<variable name="B" public_interface="out" units="per_minute"/>
<variable name="lambda1" public_interface="in" units="per_minute"/>
<variable name="lambda2" public_interface="in" units="per_minute"/>
<variable name="x2_min" public_interface="in" units="picomole"/>
<variable name="k" public_interface="in" units="flux"/>
<math xmlns="http://www.w3.org/1998/Math/MathML" id="9">
<apply>
<eq/>
<ci> B </ci>
<apply>
<divide/>
<apply>
<times/>
<ci> lambda1 </ci>
<ci> lambda2 </ci>
<ci> x2_min </ci>
</apply>
<apply>
<minus/>
<ci> k </ci>
<apply>
<times/>
<ci> lambda2 </ci>
<ci> x2_min </ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="model_parameters">
<variable name="lambda1" public_interface="out" initial_value="0.0125" units="per_minute"/>
<variable name="lambda2" public_interface="out" initial_value="0.5595" units="per_minute"/>
<variable name="m1" public_interface="out" initial_value="112.5200" units="dimensionless"/>
<variable name="m2" public_interface="out" initial_value="15.00" units="dimensionless"/>
<variable name="R" public_interface="out" initial_value="1.2162" units="millimolar"/>
<variable name="beta" public_interface="out" initial_value="1e6" units="litre_per_millimole"/>
<variable name="x1_n" public_interface="out" initial_value="490.7800" units="picomole"/>
<variable name="x2_n" public_interface="out" initial_value="6.6290" units="picomole"/>
<variable name="x2_min" public_interface="out" initial_value="0.6697" units="picomole"/>
<variable name="x2_max" public_interface="out" initial_value="14.0430" units="picomole"/>
<variable name="t0" public_interface="out" initial_value="575.0" units="minute"/>
<variable name="alpha" public_interface="out" initial_value="0.0442" units="per_minute"/>
<variable name="Ca_0" public_interface="out" initial_value="1.255" units="millimolar"/>
<variable name="Ca_1" public_interface="out" initial_value="0.1817" units="millimolar"/>
</component>
<connection>
<map_components component_1="x1" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="x2" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Ca" component_2="environment"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="x1" component_2="k"/>
<map_variables variable_1="k" variable_2="k"/>
</connection>
<connection>
<map_components component_1="x1" component_2="lambda_Ca"/>
<map_variables variable_1="lambda_Ca" variable_2="lambda_Ca"/>
</connection>
<connection>
<map_components component_1="x1" component_2="model_parameters"/>
<map_variables variable_1="lambda1" variable_2="lambda1"/>
</connection>
<connection>
<map_components component_1="lambda_Ca" component_2="A"/>
<map_variables variable_1="A" variable_2="A"/>
</connection>
<connection>
<map_components component_1="lambda_Ca" component_2="B"/>
<map_variables variable_1="B" variable_2="B"/>
</connection>
<connection>
<map_components component_1="lambda_Ca" component_2="Ca"/>
<map_variables variable_1="Ca" variable_2="Ca"/>
</connection>
<connection>
<map_components component_1="lambda_Ca" component_2="S"/>
<map_variables variable_1="S" variable_2="S"/>
</connection>
<connection>
<map_components component_1="lambda_Ca" component_2="m_Ca"/>
<map_variables variable_1="m_Ca" variable_2="m_Ca"/>
</connection>
<connection>
<map_components component_1="m_Ca" component_2="Ca"/>
<map_variables variable_1="Ca" variable_2="Ca"/>
</connection>
<connection>
<map_components component_1="m_Ca" component_2="model_parameters"/>
<map_variables variable_1="m1" variable_2="m1"/>
<map_variables variable_1="m2" variable_2="m2"/>
<map_variables variable_1="beta" variable_2="beta"/>
<map_variables variable_1="R" variable_2="R"/>
</connection>
<connection>
<map_components component_1="x2" component_2="lambda_Ca"/>
<map_variables variable_1="lambda_Ca" variable_2="lambda_Ca"/>
</connection>
<connection>
<map_components component_1="x2" component_2="x1"/>
<map_variables variable_1="x1" variable_2="x1"/>
</connection>
<connection>
<map_components component_1="x2" component_2="model_parameters"/>
<map_variables variable_1="lambda2" variable_2="lambda2"/>
</connection>
<connection>
<map_components component_1="Ca" component_2="model_parameters"/>
<map_variables variable_1="Ca_0" variable_2="Ca_0"/>
<map_variables variable_1="Ca_1" variable_2="Ca_1"/>
<map_variables variable_1="alpha" variable_2="alpha"/>
<map_variables variable_1="t0" variable_2="t0"/>
</connection>
<connection>
<map_components component_1="k" component_2="model_parameters"/>
<map_variables variable_1="lambda1" variable_2="lambda1"/>
<map_variables variable_1="lambda2" variable_2="lambda2"/>
<map_variables variable_1="x1_n" variable_2="x1_n"/>
<map_variables variable_1="x2_n" variable_2="x2_n"/>
</connection>
<connection>
<map_components component_1="S" component_2="A"/>
<map_variables variable_1="A" variable_2="A"/>
</connection>
<connection>
<map_components component_1="S" component_2="B"/>
<map_variables variable_1="B" variable_2="B"/>
</connection>
<connection>
<map_components component_1="S" component_2="m_Ca"/>
<map_variables variable_1="m_Ca" variable_2="m_Ca"/>
</connection>
<connection>
<map_components component_1="S" component_2="model_parameters"/>
<map_variables variable_1="Ca_0" variable_2="Ca_0"/>
<map_variables variable_1="x1_n" variable_2="x1_n"/>
<map_variables variable_1="x2_n" variable_2="x2_n"/>
<map_variables variable_1="lambda2" variable_2="lambda2"/>
</connection>
<connection>
<map_components component_1="A" component_2="k"/>
<map_variables variable_1="k" variable_2="k"/>
</connection>
<connection>
<map_components component_1="A" component_2="model_parameters"/>
<map_variables variable_1="x2_max" variable_2="x2_max"/>
<map_variables variable_1="lambda1" variable_2="lambda1"/>
<map_variables variable_1="lambda2" variable_2="lambda2"/>
</connection>
<connection>
<map_components component_1="B" component_2="k"/>
<map_variables variable_1="k" variable_2="k"/>
</connection>
<connection>
<map_components component_1="B" component_2="model_parameters"/>
<map_variables variable_1="x2_min" variable_2="x2_min"/>
<map_variables variable_1="lambda1" variable_2="lambda1"/>
<map_variables variable_1="lambda2" variable_2="lambda2"/>
</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#">
<rdf:Description rdf:about="">
<dc:creator rdf:parseType="Resource">
<vCard:N rdf:parseType="Resource">
<vCard:Family>Lloyd</vCard:Family>
<vCard:Given>Catherine</vCard:Given>
<vCard:Other>May</vCard:Other>
</vCard:N>
<vCard:EMAIL rdf:parseType="Resource">
<rdf:value>c.lloyd@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>
</rdf:Description>
<rdf:Description rdf:about="#shrestha_2010">
<bqs:reference rdf:parseType="Resource">
<dc:subject rdf:parseType="Resource">
<bqs:subject_type>keyword</bqs:subject_type>
<rdf:value>
<rdf:Bag>
<rdf:li>endocrine</rdf:li>
<rdf:li>parathyroid hormone</rdf:li>
<rdf:li>PTH axis</rdf:li>
</rdf:Bag>
</rdf:value>
</dc:subject>
</bqs:reference>
<bqs:reference rdf:parseType="Resource">
<bqs:Pubmed_id>20406649</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>Shrestha</vCard:Family>
<vCard:Given>R</vCard:Given>
<vCard:Other>P</vCard:Other>
</vCard:N>
</bqs:Person>
</rdf:li>
<rdf:li rdf:parseType="Resource">
<bqs:Person rdf:parseType="Resource">
<vCard:N rdf:parseType="Resource">
<vCard:Family>Hollot</vCard:Family>
<vCard:Given>C</vCard:Given>
<vCard:Other>V</vCard:Other>
</vCard:N>
</bqs:Person>
</rdf:li>
<rdf:li rdf:parseType="Resource">
<bqs:Person rdf:parseType="Resource">
<vCard:N rdf:parseType="Resource">
<vCard:Family>Chipkin</vCard:Family>
<vCard:Given>S</vCard:Given>
<vCard:Other>R</vCard:Other>
</vCard:N>
</bqs:Person>
</rdf:li>
<rdf:li rdf:parseType="Resource">
<bqs:Person rdf:parseType="Resource">
<vCard:N rdf:parseType="Resource">
<vCard:Family>Schmitt</vCard:Family>
<vCard:Given>C</vCard:Given>
<vCard:Other>P</vCard:Other>
</vCard:N>
</bqs:Person>
</rdf:li>
<rdf:li rdf:parseType="Resource">
<bqs:Person rdf:parseType="Resource">
<vCard:N rdf:parseType="Resource">
<vCard:Family>Chait</vCard:Family>
<vCard:Given>Y</vCard:Given>
</vCard:N>
</bqs:Person>
</rdf:li>
</rdf:Seq>
</dc:creator>
<dc:title>
A mathematical model of parathyroid hormone response to acute changes in plasma ionized calcium concentration in humans
</dc:title>
<dcterms:issued rdf:parseType="Resource">
<dcterms:W3CDTF>2010-07</dcterms:W3CDTF>
</dcterms:issued>
<bqs:Journal rdf:parseType="Resource">
<dc:title>Mathematical Biosciences</dc:title>
</bqs:Journal>
<bqs:volume>226</bqs:volume>
<bqs:first_page>46</bqs:first_page>
<bqs:last_page>57</bqs:last_page>
</bqs:JournalArticle>
</bqs:reference>
</rdf:Description>
</rdf:RDF>
</model>
