Team:Slovenia/ModelingPK

From 2012.igem.org

(Difference between revisions)
 
(43 intermediate revisions not shown)
Line 2: Line 2:
<head>
<head>
<meta http-equiv="X-UA-Compatible" content="IE=edge" />
<meta http-equiv="X-UA-Compatible" content="IE=edge" />
 +
 +
<!-- back to top -->
 +
<div style="position:fixed; bottom:45px; right:30px; width:100px; height:66px; background-color:transparent;">
 +
<a href="#topofthepage">
 +
<table style="background-color:transparent;" onclick="window.location = '#topofthepage'" class="invisible" style="height:100%;">
 +
<tr class="invisible" style="background-color:transparent;">
 +
<td class="invisible" style="background-color:transparent;" valign="center">
 +
<img width="100px" src ="https://static.igem.org/mediawiki/2012/1/14/Svn12_hp_btt.png"></img>
 +
</td></tr></table>
 +
</a>
 +
</div>
 +
 +
<style type="text/css">
<style type="text/css">
Line 201: Line 214:
#cssmenu ul li > ul li{display:block; list-style:inside none; padding:0; margin:0; position:relative;}  
#cssmenu ul li > ul li{display:block; list-style:inside none; padding:0; margin:0; position:relative;}  
#cssmenu ul li > ul li a{ outline:none; display:block; position:relative; margin:0; padding:8px 20px; font:10pt Arial, Helvetica, sans-serif; color:#fff; text-decoration:none; text-shadow:1px 1px 0 rgba(0,0,0, 0.5); }  
#cssmenu ul li > ul li a{ outline:none; display:block; position:relative; margin:0; padding:8px 20px; font:10pt Arial, Helvetica, sans-serif; color:#fff; text-decoration:none; text-shadow:1px 1px 0 rgba(0,0,0, 0.5); }  
 +
#cssmenu ul li > ul li a table tr td span{ outline:none; display:block; position:relative; margin:0; padding:0px 0px; font:10pt Arial, Helvetica, sans-serif; color:#fff; text-decoration:none; text-shadow:1px 1px 0 rgba(0,0,0, 0.5); }
#cssmenu, #cssmenu > ul > li > ul > li a:hover
#cssmenu, #cssmenu > ul > li > ul > li a:hover
{ background:#043A6B;  
{ background:#043A6B;  
Line 216: Line 230:
#cssmenu > ul > li > a:hover{background:#0C5DA5;}  
#cssmenu > ul > li > a:hover{background:#0C5DA5;}  
/* end CSS navigation menu (blue) */
/* end CSS navigation menu (blue) */
 +
 +
/*new table start*/
 +
table.newtable {background-color:transparent;}
 +
td.newtable, th.newtable {background-color:transparent;}
 +
thead.newtable{ }
 +
tbody .newtable{}
 +
/*new table*/
Line 292: Line 313:
<body>
<body>
 +
<a name="topofthepage" style="background-color:transparent;"></a>
<div id="banner">
<div id="banner">
-
<a style="position:absolute; top:0px; left:490px;" href="https://2012.igem.org/Main_Page"><b>iGEM 2012</b></a>
+
<a style="position:absolute; top:0px; left:490px;" href="https://2012.igem.org/Main_Page"><b></b></a>
</div>
</div>
 +
<!--<div id="container">-->
<!--<div id="container">-->
<div id="mainmenu">
<div id="mainmenu">
 +
 +
<!-- start main menu -->
<!-- start main menu -->
Line 317: Line 342:
<li><a href='https://2012.igem.org/Team:Slovenia/TheSwitchDesignedTALregulators'><span>Designed TAL regulators</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/TheSwitchDesignedTALregulators'><span>Designed TAL regulators</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/TheSwitchMutualRepressorSwitch'><span>Mutual repressor switch</span></a></li>  
<li><a href='https://2012.igem.org/Team:Slovenia/TheSwitchMutualRepressorSwitch'><span>Mutual repressor switch</span></a></li>  
-
<li><a href='https://2012.igem.org/Team:Slovenia/TheSwitchPositiveFeedbackLoopSwitch'><span>Positive feedback loop switch</span></a></li>  
+
<li><a href='https://2012.igem.org/Team:Slovenia/TheSwitchPositiveFeedbackLoopSwitch'><table onclick="window.location = 'https://2012.igem.org/Team:Slovenia/TheSwitchPositiveFeedbackLoopSwitch';" class="newtable"><tr class="newtable"><td class="newtable"><span>Positive feedback loop switch</span></td><td class="newtable"><img style="margin-right:-15px;" width="25px" src="https://static.igem.org/mediawiki/2012/e/ee/Svn12_hp_new.png"></img></td></tr></table></a></li>
 +
    <li><a href='https://2012.igem.org/Team:Slovenia/TheSwitchControls'><table onclick="window.location = 'https://2012.igem.org/Team:Slovenia/TheSwitchControls';" class="newtable"><tr class="newtable"><td class="newtable"><span>Controls</span></td><td class="newtable"><img style="margin-right:-81px;" width="25px" src="https://static.igem.org/mediawiki/2012/e/ee/Svn12_hp_new.png"></img></td></tr></table></a></li>  
  </ul>
  </ul>
</li>
</li>
Line 325: Line 351:
<li><a href='https://2012.igem.org/Team:Slovenia/SafetyMechanismsEscapeTag'><span>Escape tag</span></a></li>  
<li><a href='https://2012.igem.org/Team:Slovenia/SafetyMechanismsEscapeTag'><span>Escape tag</span></a></li>  
<li><a href='https://2012.igem.org/Team:Slovenia/SafetyMechanismsTermination'><span>Termination</span></a></li>  
<li><a href='https://2012.igem.org/Team:Slovenia/SafetyMechanismsTermination'><span>Termination</span></a></li>  
-
<li><a href='https://2012.igem.org/Team:Slovenia/SafetyMechanismsMicrocapsuleDegradation'><span>Microcapsule degradation</span></a></li>  
+
    <li><a href="https://2012.igem.org/Team:Slovenia/SafetyMechanismsMicrocapsuleDegradation"><table  onclick="window.location = 'https://2012.igem.org/Team:Slovenia/SafetyMechanismsMicrocapsuleDegradation';" class="newtable"><tr class="newtable"><td class="newtable"><span>Microcapsule degradation</span></td><td class="newtable"><img style="margin-right:-15px;" width="25px" src="https://static.igem.org/mediawiki/2012/e/ee/Svn12_hp_new.png"></img></td></tr></table></a></li>  
  </ul>
  </ul>
</li>
</li>
Line 333: Line 359:
<li><a href='https://2012.igem.org/Team:Slovenia/ImplementationHepatitisC'><span>Hepatitis C</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/ImplementationHepatitisC'><span>Hepatitis C</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/ImplementationIschaemicHeartDisease'><span>Ischaemic heart disease</span></a></li>  
<li><a href='https://2012.igem.org/Team:Slovenia/ImplementationIschaemicHeartDisease'><span>Ischaemic heart disease</span></a></li>  
 +
    <li><a href='https://2012.igem.org/Team:Slovenia/ImplementationImpact'><table onclick="window.location = 'https://2012.igem.org/Team:Slovenia/ImplementationImpact';" class="newtable"><tr class="newtable"><td class="newtable"><span>Impact</span></td><td class="newtable"><img style="margin-right:-86px;" width="25px" src="https://static.igem.org/mediawiki/2012/e/ee/Svn12_hp_new.png"></img></td></tr></table></a></li>
 
 
  </ul>
  </ul>
Line 340: Line 367:
  <ul>
  <ul>
<li><a href='https://2012.igem.org/Team:Slovenia/Modeling'><span>Overview</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/Modeling'><span>Overview</span></a></li>
-
<li><a href='https://2012.igem.org/Team:Slovenia/ModelingPK'><span>Pharmacokinetics</span></a></li>
+
    <li><a href='https://2012.igem.org/Team:Slovenia/ModelingPK'><table onclick="window.location = 'https://2012.igem.org/Team:Slovenia/ModelingPK';" class="newtable"><tr class="newtable"><td class="newtable"><span>Pharmacokinetics</span></td><td class="newtable"><img style="margin-right:-15px;" width="25px" src="https://static.igem.org/mediawiki/2012/e/ee/Svn12_hp_new.png"></img></td></tr></table></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/ModelingMethods'><span>Modeling methods</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/ModelingMethods'><span>Modeling methods</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/ModelingMutualRepressorSwitch'><span>Mutual repressor switch</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/ModelingMutualRepressorSwitch'><span>Mutual repressor switch</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/ModelingPositiveFeedbackLoopSwitch'><span>Positive feedback loop switch</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/ModelingPositiveFeedbackLoopSwitch'><span>Positive feedback loop switch</span></a></li>
-
<li><a href='https://2012.igem.org/Team:Slovenia/ModelingQuantitativeModel'><span>Quantitative and stability model</span></a></li>  
+
<li><a href='https://2012.igem.org/Team:Slovenia/ModelingQuantitativeModel'><table onclick="window.location = 'https://2012.igem.org/Team:Slovenia/ModelingQuantitativeModel';" class="newtable"><tr class="newtable"><td class="newtable"><span>Experimental model</span></td><td class="newtable"><img style="margin-right:-15px;" width="25px" src="https://static.igem.org/mediawiki/2012/e/ee/Svn12_hp_new.png"></img></td></tr></table></a></li>  
-
<li><a href='https://2012.igem.org/Team:Slovenia/ModelingInteractiveSimulations'><span>Interactive simulations</span></a></li>
+
    <li><a href='https://2012.igem.org/Team:Slovenia/ModelingInteractiveSimulations'><table onclick="window.location = 'https://2012.igem.org/Team:Slovenia/ModelingInteractiveSimulations';" class="newtable"><tr class="newtable"><td class="newtable"><span>Interactive simulations</span></td><td class="newtable"><img style="margin-right:-15px;" width="25px" src="https://static.igem.org/mediawiki/2012/e/ee/Svn12_hp_new.png"></img></td></tr></table></a></li>
  </ul>
  </ul>
</li>
</li>
Line 356: Line 383:
  <ul>
  <ul>
<li><a href='https://2012.igem.org/Team:Slovenia/Notebook'><span>Experimental methods</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/Notebook'><span>Experimental methods</span></a></li>
-
<li><a href='https://2012.igem.org/Team:Slovenia/NotebookLablog'><span>Lablog</span></a></li>
+
    <li><a href='https://2012.igem.org/Team:Slovenia/NotebookLablog'><table onclick="window.location = 'https://2012.igem.org/Team:Slovenia/NotebookLablog';" class="newtable"><tr class="newtable"><td class="newtable"><span>Lablog</span></td><td class="newtable"><img style="margin-right:-90px;" width="25px" src="https://static.igem.org/mediawiki/2012/e/ee/Svn12_hp_new.png"></img></td></tr></table></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/NotebookLabSafety'><span>Lab safety</span></a></li>  
<li><a href='https://2012.igem.org/Team:Slovenia/NotebookLabSafety'><span>Lab safety</span></a></li>  
  </ul>
  </ul>
Line 379: Line 406:
<li><a href='https://2012.igem.org/Team:Slovenia/Team'><span>Team members</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/Team'><span>Team members</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/TeamAttributions'><span>Attributions</span></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/TeamAttributions'><span>Attributions</span></a></li>
 +
<li><a href='https://2012.igem.org/Team:Slovenia/TeamCollaborations'><table  onclick="window.location = 'https://2012.igem.org/Team:Slovenia/TeamCollaborations';" class="newtable"><tr class="newtable"><td class="newtable"><span>Collaborations</span></td><td class="newtable"><img style="margin-right:-20px;" width="25px" src="https://static.igem.org/mediawiki/2012/e/ee/Svn12_hp_new.png"></img></td></tr></table></a></li>
<li><a href='https://2012.igem.org/Team:Slovenia/TeamGallery'><span>Gallery</span></a></li>  
<li><a href='https://2012.igem.org/Team:Slovenia/TeamGallery'><span>Gallery</span></a></li>  
<li><a href='https://2012.igem.org/Team:Slovenia/TeamSponsors'><span>Sponsors</span></a></li>  
<li><a href='https://2012.igem.org/Team:Slovenia/TeamSponsors'><span>Sponsors</span></a></li>  
Line 386: Line 414:
</div>
</div>
<!-- end main menu -->
<!-- end main menu -->
 +
Line 392: Line 421:
<div id="main">
<div id="main">
<br/>
<br/>
-
<h1>Pharmacokinetic model</h1>
 
 +
<!-- dummy link na bannerju -->
 +
<a href="https://2012.igem.org/Main_Page">
 +
<div id="dummy" style="background-color:transparent; position:absolute; left:870px; top:25px; width:115px; height:80px; z-index:100; opacity:0.0;">
 +
</div>
 +
</a>
 +
 +
<h1>Pharmacokinetic model</h1>
<table class="summary">
<table class="summary">
Line 399: Line 434:
<td class="summary" style="font-size:110%;">
<td class="summary" style="font-size:110%;">
<p><b>What is pharmacokinetic model?</b><br>Pharmacokinetic model is a quantitative description of drug absorption, distribution, metabolism and elimination from the body.</p><br>
<p><b>What is pharmacokinetic model?</b><br>Pharmacokinetic model is a quantitative description of drug absorption, distribution, metabolism and elimination from the body.</p><br>
-
<p><b>Why do we need it?</b><br>- to compare conventional and our therapy<br>- to calculate required drug production in microencapsulated cells</p><br>
+
<p style="margin-bottom:0px;"><b>Why do we need it?</b>
-
<p><b>How was our pharmacokinetic model constructed?</b><br>- physiologically based<br>- organs and tissues are presented as compartments<br>- processes in the body are described with ODEs<br>- customized for chosen drug and it's administration</p><br>
+
<ul style="margin-left:15px;">
-
<p><b>For which drugs (diseases) did we test our model?</b><br>- interferon (hepatitis C)<br>- anakinra (ischemic heart disease)</p><br>
+
<li>To compare conventional and our therapy.</li>
-
<p><b>3D animation</b><br>3D animation serves as visual represention of how drug concentrations change throughout the course of the therapy in various compartments. For the implementation, we used C++ programming language and OpenGL for 3D rendering.</p><br>
+
<li>To calculate required drug production at target site.</li>
-
<p><b>Results</b><br>- our therapy does not produce concentration fluctuations; absence of concentration peaks and lows<br>- with local drug production, drug concentration rises to desired level within a day (for both tested drugs)<br>- for each drug we got an estimation of the degree of drug localization</p><br>
+
</ul>
-
<p><b>What did results tell us?</b><br>Results show us the drug is more localized; the degree of localization is dependant on characteristics of tissue (e.g. perfusion), where microencapsulated cells would be applied. Additionally, absence of high concentration peaks suggests that our therapy would have less systemic side effects than conventional methods. Since the drug would be constantly produced, drug concentration would not drop enough for the drug to be therapeutically ineffective. According to these findings, steady drug levels over time therefore imply an optimized treatment. With known drug production rate of microencapsulated cells, it is possible to adjust dosage to reach the desired drug concentration in target tissue. Therapeutic production was estimated according to activity assays with mammalian cells. After obtaining the results of pharmacokinetic simulation, we were able to calculate number of microencapsulated cells needed for the therapy.</p><br>
+
</p><br/>
 +
 
 +
 
 +
<p style="margin-bottom:0px;"><b>How was our pharmacokinetic model constructed?</b>
 +
<ul style="margin-left:15px;">
 +
<li>The model is physiologically based.</li>
 +
<li>Organs and tissues are presented as compartments.</li>
 +
<li>Processes in the body are described with ordinary differential equations (ODEs).</li>
 +
<li>The model is customized for chosen drug and its type of administration.</li>
 +
</ul></p><br/>
 +
 
 +
<p style="margin-bottom:0px;"><b>For which drugs (diseases) did we test our model?</b>
 +
<ul style="margin-left:15px;">
 +
<li>interferon (hepatitis C)</li>
 +
<li>anakinra (ischemic heart disease)</li>
 +
</ul>
 +
</p><br/>
 +
 
 +
 
 +
<p><b>3D animation</b><br>3D animation serves as a real-time visual representation of how drug concentrations change throughout the course of the therapy in various compartments. The animation was implemented in C++ programming language and OpenGL was used for 3D rendering.</p><br>
 +
<img style="width:70%;"src="https://static.igem.org/mediawiki/2012/2/2e/Svn12_hcv_cmp_-_Copy.png"></img><br/>
 +
 
 +
 
 +
<p style="margin-bottom:0px;"><b>Results</b>
 +
<ul style="margin-left:15px;">
 +
<li>Our therapy does not produce concentration fluctuations - absence of concentration peaks and lows.</li>
 +
<li>With local drug production, drug concentration rises to a desired level within a day (for both tested drugs).</li>
 +
<li>For each drug, we obtained an estimation of the degree of drug localization.</li>
 +
</ul>
 +
</p><br/>
 +
 
 +
 
 +
<p><b>What did the results tell us?</b>
 +
<br>Results show us the drug is more localized; the degree of localization depends on characteristics of tissue (e.g. perfusion) where microencapsulated cells would be applied. Additionally, absence of high concentration peaks suggests that our therapy would have less systemic side effects than conventional methods. Since the drug would be constantly produced, drug concentration would not drop enough for the drug to be therapeutically ineffective. According to these findings, steady drug levels over time therefore imply an optimized treatment. With known drug production rate of microencapsulated cells, it is possible to adjust the dosage to reach the desired drug concentration in the target tissue. Therapeutic production was estimated according to activity assays with mammalian cells. After obtaining the results of pharmacokinetic simulation, we were able to calculate the number of microencapsulated cells needed for the therapy.</p>
</td>
</td>
</tr>
</tr>
</table>
</table>
-
<img style="width:70%;"src="https://static.igem.org/mediawiki/2012/2/2e/Svn12_hcv_cmp_-_Copy.png"></img>
 
-
<p>We developed a 3D visualization of pharmacokinetic model data representing our findings in that aspect.</p>
+
 
 +
<p></p>
<h2>Introduction</h2>
<h2>Introduction</h2>
-
<p>The major downsides of standard interferon treatment are substantial side effects that to a large degree are the consequence of very high drug concentrations, occurring shortly after drug administration. If high concentration peaks could be avoided and lower levels maintained steadily over time, this would result in reduced side effects without compromising therapeutic effectiveness. This idea is currently being tested in clinical trials using the interferon infusion pump (<i>COPE-HCV: Phase 2, randomized, open-label, active-control, dose-ranging study of interferon alfa-2b given via continuous sub-Q infusion; trial by Medronic Inc.</i>)</p>
+
<p>The major downsides of standard interferon treatment are substantial side effects that to a large degree are the consequence of very high drug concentrations, occurring shortly after drug administration. If high concentration peaks could be avoided and lower levels maintained steadily over time, this would result in reduced side effects without compromising therapeutic effectiveness. This idea is currently being tested in clinical trials using the interferon infusion pump (<i>COPE-HCV: Phase 2, randomized, open-label, active-control, dose-ranging study of interferon alfa-2b given via continuous sub-Q infusion; trial by Medronic Inc.</i>).</p>
<p>One of our goals is to show that therapy with drug-producing cells is more beneficial than standard treatments based on drug injections, which are in use today.  
<p>One of our goals is to show that therapy with drug-producing cells is more beneficial than standard treatments based on drug injections, which are in use today.  
Line 421: Line 489:
<p>We tested this hypothesis and compared standard therapies with our proposed treatment.
<p>We tested this hypothesis and compared standard therapies with our proposed treatment.
<br>  
<br>  
-
Because of complex physiological mechanisms and an extensive set of biological parameters that can not be accurately measured for either ethical or technical reasons, we developed a model that covers the most crucial aspects and processes. At the same time, computer simulation provides simpler and faster option than <i>in vivo</i> research.<p/>
+
Because of complex physiological mechanisms and an extensive set of biological parameters that cannot be accurately measured for either ethical or technical reasons, we developed a model that covers the most crucial aspects and processes. At the same time, computer simulation provides simpler and faster option than <i>in vivo</i> research.<p/>
 +
 
 +
<h2>3D visualization</h2>
 +
<p>We implemented a real-time 3D visualization of the pharmacokinetic simulation for hepatitis C and ischemia. This way, we were able to visually represent how drug concentrations change throughout the course of the therapy in various compartments, such as liver. In our 3D visualizations, drug concentration levels are represented as red color intensities - the compartments with high concentrations are closer to saturated red color, while the compartments with low concentrations are closer to white color. Red color intensity of a compartment was calculated as a ratio of the compartment's drug concentration at a given time to the maximal reached concentration of any compartment during the course of the therapy, or therapies compared.</p>
 +
<p>For implementation, we used C++ programming language and OpenGL for 3D rendering. The 3D anatomical models, obtained from <a href="http://dbarchive.biosciencedbc.jp/en/bodyparts3d/desc.html">BodyParts3D</a> database, were optimized for real-time rendering using <a href="http://meshlab.sourceforge.net/">MeshLab</a> and rendered as vertex arrays in OpenGL. Compartment concentration values as a function of time were obtained as vectors from the pharmacokinetic simulation and integrated into the application.</p><br/>
<h2>A physiologically based pharmacokinetic model</h2>
<h2>A physiologically based pharmacokinetic model</h2>
<p>A pharmacokinetic model is a quantitative description of drug absorption, distribution, metabolism and elimination from the body. A model is defined by a system of ordinary differential equations to represent essential drug kinetics.
<p>A pharmacokinetic model is a quantitative description of drug absorption, distribution, metabolism and elimination from the body. A model is defined by a system of ordinary differential equations to represent essential drug kinetics.
<br>
<br>
-
We took into account the neccessary biological parameters and based the simulated processes on actual physiological mechanisms to construct a physiologically plausible model.</p>
+
We took into account the necessary biological parameters and based the simulated processes on actual physiological mechanisms to construct a physiologically plausible model.</p>
<h2>Selecting the optimal physiological model</h2>
<h2>Selecting the optimal physiological model</h2>
-
<p>A physiologically based model is composed of multiple compartments which represent organs of the body. Parts were chosen in accordance with drug and tissue specifics, so that the relevant organs are represented as separate compartments, while other tissues were merged on the basis of common characteristics.</p>
+
<p>A physiologically based model is composed of multiple compartments which represent organs of the body. Parts were chosen in accordance with drug and tissue specifics, so that the relevant organs are represented as separate compartments, while other tissues were merged on the basis of common characteristics (Levitt et al., 2006).</p>
<h3>Determing compartments</h3>
<h3>Determing compartments</h3>
-
<p>We used the following characteristics as a criteria for splitting and merging organs:
+
<p>We used the following characteristics as criteria for splitting and merging organs:
<ul style="margin-left:30px;">
<ul style="margin-left:30px;">
<li>Liver is the target organ of therapy.</li>
<li>Liver is the target organ of therapy.</li>
Line 441: Line 513:
<li>Skin, musle and adipose tissue have similar, slow blood perfusion.</li>
<li>Skin, musle and adipose tissue have similar, slow blood perfusion.</li>
<li>Gut, spleen and heart are all rapidly perfused tissues.</li>
<li>Gut, spleen and heart are all rapidly perfused tissues.</li>
-
<li>Interferon is mainly eliminated via renal catabolism, while hepatic metabolism accounts only for a minor pathway of elimination.</li></ul></p>
+
<li>Interferon is mainly eliminated via renal catabolism, while hepatic metabolism accounts only for a minor pathway of elimination. </li></ul></p>
-
<p>We decided to define separate compartments for the liver, kidney and lungs. All other rapidly perfused tissues are grouped together as one part. Since venous blood enters the lungs and arterial blood flows into all other organs, we separately simulate venous and arterial blood. Because interferon does not cross the blood-brain barrier it is not necessary for the brain to be modelled separately. Skin, muscle, fat and other slowly perfused tissues are merged together into one compartment.</p>
+
<p>We decided to define separate compartments for the liver, kidney and lungs. All other rapidly perfused tissues are grouped together as one part. Since venous blood enters the lungs and arterial blood flows into all other organs, we separately simulate venous and arterial blood. Because interferon does not cross the blood-brain barrier it is not necessary for the brain to be modeled separately. Skin, muscle, fat and other slowly perfused tissues are merged together into one compartment.</p>
-
<p>We constructed three final models - two for standard interferon treatments and a third for a prospective therapy with drug-producing microencapsulated cells. The fundamental design is the same in all models, they are only modified for specific entry points of the drug and the corresponding absorption or production processes. On the diagram, blocks representing different administrations are shown in distinct colors: blue for the intravenous bolus, green for the subcutaneous injection and red for the interferon production by microencapsulated cells that are implanted into liver.</p>
+
<p>We constructed three final models - two for standard interferon treatments and a third for a prospective therapy with drug-producing microencapsulated cells. The fundamental design is the same in all models; they are only modified for specific entry points of the drug and the corresponding absorption or production processes. On the diagram, blocks representing different administrations are shown in distinct colors: blue for the intravenous bolus, green for the subcutaneous injection and red for the interferon production by microencapsulated cells that are implanted into liver (Figure 1).</p>
<p>
<p>
</br>
</br>
Line 451: Line 523:
<center><img src="https://static.igem.org/mediawiki/2012/6/6b/SVN12_PK_1.png" style="height:800px; width:auto;"/></center>
<center><img src="https://static.igem.org/mediawiki/2012/6/6b/SVN12_PK_1.png" style="height:800px; width:auto;"/></center>
</br></p>
</br></p>
-
<center><b>Figure 1.</b> Multi-compartmental physiologically based pharmacokinetic model, customized for different Interferon therapies. Main difference is the route of administration - for each administration, there is different compartment from which drug enters the body. Blue blocks are used in modeling for intravenous administration, green blocks are specific for subcutaneous administration and red blocks represent continuous production of the drug, in this case in the liver compartment. Therapies of course don't differ only in drug entry points, but also in absorption or production processes, percent of drug that enters the blood stream and so on.</center><br><br>
+
<center><b>Figure 1.</b> Multi-compartment physiologically based pharmacokinetic model, customized for different Interferon therapies. Main difference is the route of administration - for each administration, there is different compartment from which drug enters the body. Blue blocks are used in modeling for intravenous administration, green blocks are specific for subcutaneous administration and red blocks represent continuous production of the drug, in this case in the liver compartment. Therapies of course don't differ only in drug entry points, but also in absorption or production processes, percent of drug that enters the blood stream and so on.</center><br><br>
</br>
</br>
Line 458: Line 530:
<p>Types of parameters used:</p>
<p>Types of parameters used:</p>
-
<h3>1.) Species specific</h3>
+
<h3>1.) Species specific parameters</h3>
<p>Q<sub>i</sub> – blood flows to tissues
<p>Q<sub>i</sub> – blood flows to tissues
<br>
<br>
Line 479: Line 551:
                 </tr>
                 </tr>
                 <thead class="normal">
                 <thead class="normal">
-
<tr class="normal"><td class="normal">Comparment</td>
+
<tr class="normal"><td class="normal">Compartment</td>
                                   <td class="normal">Parameter name</td>
                                   <td class="normal">Parameter name</td>
                                   <td class="normal">Percent cardiac output</td>
                                   <td class="normal">Percent cardiac output</td>
Line 540: Line 612:
                 </tr>
                 </tr>
                 <thead class="normal">
                 <thead class="normal">
-
<tr class="normal"><td class="normal">Comparment</td>
+
<tr class="normal"><td class="normal">Compartment</td>
                                   <td class="normal">Parameter name</td>
                                   <td class="normal">Parameter name</td>
                                   <td class="normal">Percent body volume</td>
                                   <td class="normal">Percent body volume</td>
Line 587: Line 659:
<p>* average body density = 1 kg /L</p>
<p>* average body density = 1 kg /L</p>
-
<h3>2.) Individual specific</h3>
+
<h3>2.) Individual specific parameters</h3>
<p>BW - body weight
<p>BW - body weight
<br>
<br>
Line 601: Line 673:
Values used in our model present an average adult male, weighing 70 kg, with mean cardiac output and normal renal function. </p>
Values used in our model present an average adult male, weighing 70 kg, with mean cardiac output and normal renal function. </p>
-
<h3>3.) Chemical specific </h3>
+
<h3>3.) Chemical specific parameters</h3>
<p>t<sub>1/2</sub>  - drug half life
<p>t<sub>1/2</sub>  - drug half life
<br>
<br>
Line 608: Line 680:
P<sub>blood:tissue</sub> - partition coefficient</p>
P<sub>blood:tissue</sub> - partition coefficient</p>
-
<p>Interferon is expected to be found only in the plasma and not in red blood cells, therefore we can conclude that the amount of interferon found in plasma is equal to the amount in blood. Since blood is comprised of four parts plasma and three parts red blood cells, we can calculate the blood to plasma ratio to be:</p>
+
<p>Interferon is expected to be found only in the plasma and not in red blood cells; therefore we can conclude that the amount of interferon found in plasma is equal to the amount in blood. Since blood is comprised of four parts plasma and three parts red blood cells, we can calculate the blood to plasma ratio to be:</p>
<p><center><img src="https://static.igem.org/mediawiki/2012/0/0a/Svn12_ENACBE_Pplasma.gif"/></center></p>
<p><center><img src="https://static.igem.org/mediawiki/2012/0/0a/Svn12_ENACBE_Pplasma.gif"/></center></p>
Line 618: Line 690:
<b>Uncertainty of partition coefficients</b>
<b>Uncertainty of partition coefficients</b>
-
<p>A critical element in human PBPK modelling is the uncertainty of values of the partition coefficients. Partition coefficients are an important aspect of pharmacokinetic modeling, because they denote how the drug distributes throughout body tissues. The value of each coefficient has a complex dependence on solubility, permeability, pH, binding affinity of the drug to the receptor receptor and many other factors. The difference in a few amino acids between subtypes of interferon alpha can impact values of coefficients quite noticably. Even so, there are deviations of evaluated drug distribution of the same subtype of interferon alpha, depending on the detection method (radioactivity, ELISA). </p>
+
<p>A critical element in human physiologically based pharmacokinetic modeling is the uncertainty of values of the partition coefficients. Partition coefficients are an important aspect of pharmacokinetic modeling, because they denote how the drug distributes throughout body tissues. The value of each coefficient has a complex dependence on solubility, permeability, pH, binding affinity of the drug to the receptor receptor and many other factors. The difference in a few amino acids between subtypes of interferon alpha can impact values of coefficients quite noticeably. Even so, there are deviations of evaluated drug distribution of the same subtype of interferon alpha, depending on the detection method (radioactivity, ELISA). </p>
-
<p>For legal and ethical reasons these values cannot be directly measured in humans. We had to rely on various studies of interferon tissue distribution in rodents to calculate partition coefficients. It is generally assumed that animal and human partition coefficients are similiar for the same kind of tissue.</p>
+
<p>For legal and ethical reasons these values cannot be directly measured in humans. We had to rely on various studies of interferon tissue distribution in rodents to calculate partition coefficients (Johns et al, 1990). It is generally assumed that animal and human partition coefficients are similar for the same kind of tissue.</p>
-
<b>Chemical and ROA* specific</b>
+
<b>Chemical and ROA* specific parameters</b>
<p>* ROA - route of administration</p>
<p>* ROA - route of administration</p>
<p>D – dose
<p>D – dose
Line 635: Line 707:
<h2>General mass-balance equations</h2>
<h2>General mass-balance equations</h2>
-
<p>The equation below describes the change in concentration over time in non-eliminating tissues. The equation has the same form for both rapidly and slowly perfused tissue. Each compartment is then described with it's distinctive values for blood flow, concentration and partition coefficient.</p>
+
<p>The equation below describes the change in concentration over time in non-eliminating tissues. The equation has the same form for both rapidly and slowly perfused tissue. Each compartment is then described with its distinctive values for blood flow, concentration and partition coefficient.</p>
<p><center><img src="https://static.igem.org/mediawiki/2012/a/a3/Svn12_ENACBE_tissue.gif"/></center><br></p>
<p><center><img src="https://static.igem.org/mediawiki/2012/a/a3/Svn12_ENACBE_tissue.gif"/></center><br></p>
Line 656: Line 728:
<h2>Modeling pharmacokinetic processes</h2>
<h2>Modeling pharmacokinetic processes</h2>
-
<h3>Absorption</h3>
+
<h3>1.) Absorption</h3>
-
<p>Absorption processes depend on the route of administration, </p>
+
<p>Absorption processes depend on the route of administration: </p>
-
<b>Intravenous administration</b>
+
<b> a.) Intravenous administration</b><br>
-
<p>The drug is injected directly into the blood stream, so there is no special absorption process. The dose is bolus and enters the system completely.</p>
+
<p>The drug is injected directly into the blood stream, so there is no special absorption process. The dose is bolus and enters the system completely.</p><br>
-
<b>Subcutaneous administration</b>
+
<b> b.) Subcutaneous administration</b><br>
<p>In subcutaneous administration, more than 80% of the initial dose enters systemic circulation. This percentage is defined as bioavailability and can reach up to 95% for interferon alpha given subcutaneously.
<p>In subcutaneous administration, more than 80% of the initial dose enters systemic circulation. This percentage is defined as bioavailability and can reach up to 95% for interferon alpha given subcutaneously.
<br>
<br>
Line 680: Line 752:
-
<p>Constant k<sub>a</sub> was taken from literature and the value of k<sub>0</sub> was calculated using paramteres from the same source. Fz denotes the portion of the drug that is absorbed by zero order process.</p>
+
<p>Constant k<sub>a</sub> was taken from literature and the value of k<sub>0</sub> was calculated using parameters from the same source. Fz denotes the portion of the drug that is absorbed by zero order process.</p>
<p><center><img src="https://static.igem.org/mediawiki/2012/a/ac/Svn12_ENACBE_k0.gif"/></center><br></p>
<p><center><img src="https://static.igem.org/mediawiki/2012/a/ac/Svn12_ENACBE_k0.gif"/></center><br></p>
-
<h3>Distribution</h3>
+
<h3>2.) Distribution</h3>
<b>Estimating partition coefficients</b><br>
<b>Estimating partition coefficients</b><br>
-
<p>Distribution process is usually modelled indirectly with the use of partition coefficients, which are drug and tissue specific. We derived partition coefficients according to the available literature. However, these values do differ from study to study (e.g. that kidney concentration is 1-9x times that of plasma), but they overall agree that interferon concentration is highest in the kidney, very high in the liver and lungs and to some extent in other rapidly perfused tissues, but lower in the adipose and other slowly perfused tissues.</p>
+
<p>Distribution process is usually modeled indirectly with the use of partition coefficients, which are drug and tissue specific. We derived partition coefficients according to the tissue distribution data in available literature (Thitinian et al., 2009; Johns et al., 1990). However, these values do differ from study to study (e.g. that kidney concentration is 1-9x times that of plasma (Bohoslawec et al., 1986), but they overall agree that after  subcutaneous or intravenous administration, interferon concentration is highest in the kidney, very high in the liver and lungs and to some extent in other rapidly perfused tissues, but lower in the adipose and other slowly perfused tissues.</p>
<p><center><img src="https://static.igem.org/mediawiki/2012/9/95/Svn12_enacbeKONC.gif"/></center><br></p>
<p><center><img src="https://static.igem.org/mediawiki/2012/9/95/Svn12_enacbeKONC.gif"/></center><br></p>
<table class="normal" style="font-size:90%; width:60%; text-align:center;">
<table class="normal" style="font-size:90%; width:60%; text-align:center;">
<tbody  class="normal">
<tbody  class="normal">
-
<thead class="normal"><td class="normal">Comparment</td>
+
<thead class="normal"><td class="normal">Compartment</td>
                                   <td class="normal">Parameter name</td>
                                   <td class="normal">Parameter name</td>
                                   <td class="normal">Value</td>  
                                   <td class="normal">Value</td>  
Line 724: Line 796:
</table>  
</table>  
-
<h3>Elimination</h3>
+
<h3>3.) Elimination</h3>
<p>In the literature the terminal half-life for interferon alpha is from 3-8 h, with a mean around 5 h. Interferon alpha is detectable in the plasma for 4-8 h after the rapid intravenous injection and for 16-30h after subcutaneous administration.
<p>In the literature the terminal half-life for interferon alpha is from 3-8 h, with a mean around 5 h. Interferon alpha is detectable in the plasma for 4-8 h after the rapid intravenous injection and for 16-30h after subcutaneous administration.
<br>
<br>
Line 734: Line 806:
<p></p>
<p></p>
-
<h3>Modeling subcutaneous therapy</h3>
+
<h3>Interferon alpha subcutaneous therapy modeling</h3>
<p>We simulated common type of interferon administration:
<p>We simulated common type of interferon administration:
 +
<ul style="margin-left:30px;">
<li> 3 times per week, subcutaneous administration </li>
<li> 3 times per week, subcutaneous administration </li>
<li> 3 MU = 11100 ng, dose </li>
<li> 3 MU = 11100 ng, dose </li>
<li> F = 93%, bioavailability</li>
<li> F = 93%, bioavailability</li>
<li> Ko, Ka - zero and first order absorption</li>
<li> Ko, Ka - zero and first order absorption</li>
-
<li> t<sub>1/2</sub> = 5h, half life</li>
+
<li> t<sub>1/2</sub> = 5h, half life</li></ul>
-
<p><br>We constructed the model according to data from the literature. Results are shown in graph below.</p><br>
+
<p><br>We constructed the model according to data from the literature (Roferon A, Hoffman-La Roche; Chatelut, 1999; Thitinian et al., 2009; Marcellin et al., 2003). Results are shown in graph below (Figure 2).</p><br>
<p><center><img src="https://static.igem.org/mediawiki/2012/4/4c/Sv12_Hepatitis_SC_v1.png"/></center></p>
<p><center><img src="https://static.igem.org/mediawiki/2012/4/4c/Sv12_Hepatitis_SC_v1.png"/></center></p>
<center><b>Figure 2.</b> Standard subcutaneous therapy: high peaks right after administration, followed by rapid concentration decrease. </center><br><br>
<center><b>Figure 2.</b> Standard subcutaneous therapy: high peaks right after administration, followed by rapid concentration decrease. </center><br><br>
Line 747: Line 820:
<h3>Calculating desired drug production</h3>
<h3>Calculating desired drug production</h3>
-
<p>Our model assumes no cell divisions and therefore a constant therapeutic cell count in the body. In this case, cell production of a drug is mathematically equivallent to a continuous infusion of a drug into the tissue, where cells are located.
+
<p>Our model assumes no cell divisions and therefore a constant therapeutic cell count in the body. In this case, cell production of a drug is mathematically equivalent to a continuous infusion of a drug into the tissue, where cells are located.
-
To calculate the desired concentration level in the target tissue, we calculated AUC * of liver concentrations for subcutaneous administration over a period of one week. This allowed us to estimate the average concentration in the liver which corresponds to the therapeutic level we wish to reach and maintain with cellular drug production. The combination of kidney elimination and constant production leads to stable drug levels in tissues. After about four half-lives, drug concentrations rise to final levels and remain steady as long as the production rate stays the same.  </p><br>
+
To calculate the desired concentration level in the target tissue, we calculated AUC (area under the curve) of liver concentrations for subcutaneous administration over a period of one week. This allowed us to estimate the average concentration in the liver which corresponds to the therapeutic level we wish to reach and maintain with cellular drug production. The combination of kidney elimination and constant production leads to stable drug levels in tissues. After about four half-lives, drug concentrations rise to final levels and remain steady as long as the production rate stays the same (Figure 3).  </p><br>
<p><center><img src="https://static.igem.org/mediawiki/2012/d/dc/Svn12_Hepatitis_HEK_v1.png"/></center><br></p>
<p><center><img src="https://static.igem.org/mediawiki/2012/d/dc/Svn12_Hepatitis_HEK_v1.png"/></center><br></p>
<center><b>Figure 3.</b> Drug distribution troughout tissues with Switch-IT therapy.</center><br><br>
<center><b>Figure 3.</b> Drug distribution troughout tissues with Switch-IT therapy.</center><br><br>
Line 756: Line 829:
<p>If we compare result of different treatments it clearly shows that local therapy with drug producing cells does not result in concentration fluctuations, which are present in common therapies.
<p>If we compare result of different treatments it clearly shows that local therapy with drug producing cells does not result in concentration fluctuations, which are present in common therapies.
<br>
<br>
-
Concentrations never reach levels as high as in subcutaneous or intravenous administration. On the other hand, concentrations do not drop to a point wich would prevent the drug from becoming therapeutically inefficient.</p>
+
Concentrations never reach levels as high as in subcutaneous or intravenous administration. On the other hand, concentrations do not drop to a point, which would prevent the drug from becoming therapeutically inefficient.</p>
-
<p>When the drug is produced at the target site, concentration levels between target tissues and other organs are inclined in favor of the target tissue - although kidneys show a substantial uptake of the drug in subcutaneous administration, the relative difference between kidney and liver concentrations are smaller with localized therapy.</p><br>
+
<p>When the drug is produced at the target site, concentration levels between target tissues and other organs are inclined in favor of the target tissue - although kidneys show a substantial uptake of the drug in subcutaneous administration, the relative difference between kidney and liver concentrations are smaller with localized therapy (Figure 4).</p><br>
<p><center><img src="https://static.igem.org/mediawiki/2012/c/c5/Sv12_Hepatitis_SCinHEK_v1.png"/></center><br></p>
<p><center><img src="https://static.igem.org/mediawiki/2012/c/c5/Sv12_Hepatitis_SCinHEK_v1.png"/></center><br></p>
Line 766: Line 839:
<p>When we calculate the desired steady concentration in the target tissue, it was possible to estimate how many cells  would need to be implanted. From the required production rate, which is thought to be linearly proportional to the cell count, we can calculate the total number of cells needed. Considering the average number of cells per microcapsule and mean capsule size, we can calculate the amount of microcapsules we would need to inject in order to attain the wanted therapeutic effect. </p>
<p>When we calculate the desired steady concentration in the target tissue, it was possible to estimate how many cells  would need to be implanted. From the required production rate, which is thought to be linearly proportional to the cell count, we can calculate the total number of cells needed. Considering the average number of cells per microcapsule and mean capsule size, we can calculate the amount of microcapsules we would need to inject in order to attain the wanted therapeutic effect. </p>
-
<h2>Modeling anakinra</h2>
+
<h2>Anakinra therapy modeling</h2>
-
<p>For modeling the local application of anakinra-producing microencapsulated cells in the ischaemic muscle tissue we adopted slightly different approach. Liver is perfused by sinusoid capillaries which allow for passage of macromolecules to the tissue, therefore even protein therapeutics can be rapidly equilibrated with blood and distributed around the organism. Muscle tissue, however, is perfused with continuous capillaries, which prevents fast passage of macromolecules from the tissue to the blood. Protein molecules pass the continuous capillaries only through slow trans-pinocytosis. Consequently this increases the local concentration of therapeutics permitting the achievement of high local concentration with relatively little loss to the system and elimination by kidneys.</p>
+
<p>For modeling the local application of anakinra-producing microencapsulated cells in the ischaemic muscle tissue we adopted slightly different approach. Liver is perfused by sinusoid capillaries which allow for passage of macromolecules to the tissue, therefore even protein therapeutics can be rapidly equilibrated with blood and distributed around the organism. Muscle tissue, however, is perfused with continuous capillaries, which prevents fast passage of macromolecules from the tissue to the blood (Figure 5). Protein molecules pass the continuous capillaries only through slow trans-pinocytosis. Consequently this increases the local concentration of therapeutics permitting the achievement of high local concentration with relatively little loss to the system and elimination by kidneys.</p>
<br>
<br>
<p><center><img src="https://static.igem.org/mediawiki/2012/5/54/Zile.jpg"/></center></p>
<p><center><img src="https://static.igem.org/mediawiki/2012/5/54/Zile.jpg"/></center></p>
<center><b>Figure 5. </b>Different capillary types in muscle tissue (continuous capillary, left) and in liver (sinusoid capillary, right). Image from: (Mescher, 2009)</center><br><br>
<center><b>Figure 5. </b>Different capillary types in muscle tissue (continuous capillary, left) and in liver (sinusoid capillary, right). Image from: (Mescher, 2009)</center><br><br>
-
<p>Therefore we took this feature of the tissue into account when calculating partition coefficients between tissue and blood compartment. Myoglobin has similar molecular weight to anakinra, therefore we assumed its partition coefficient to be the same(Hall, 2010). We have calculated the therapeutic concentration of anakinra to be 25 ng/mL, which we derived from the anakinra IC_50 value (Dahlén et al., 2008).</p>
+
<p>Therefore we took this feature of the tissue into account when calculating partition coefficients between tissue and blood compartment. Myoglobin has similar molecular weight as anakinra; therefore we assumed its partition coefficient to be the same(Hall, 2010). We calculated the therapeutic concentration of anakinra to be 25 ng/mL, which we derived from the anakinra IC_50 value (Dahlén et al., 2008).</p>
-
<h3>Compartmental model and parameters</h3>
+
<h3>Compartmental model and simulation parameters</h3>
-
<p>According to a study of anakinra distribution and elimination after subcutaneous injection, we calculated relevant compartments for physiological model. At first, we modeled standard regime for anakinra therapy: 100 mg per day subcutaneously with next parameters:<br>
+
<p>According to a study of anakinra distribution and elimination after subcutaneous injection, we calculated relevant compartments for physiological model. At first, we modeled standard regime for anakinra therapy: 100 mg per day subcutaneously with the following parameters:<br>
 +
<ul style="margin-left:30px;">
<li>bioavailability: F = 95%</li>
<li>bioavailability: F = 95%</li>
<li>half life: T<sub>1/2</sub> = 5h</li>
<li>half life: T<sub>1/2</sub> = 5h</li>
<li>time to maximum concentration: T<sub>max</sub> = 4h</li>
<li>time to maximum concentration: T<sub>max</sub> = 4h</li>
-
<li>maximum concentration: C<sub>max</sub> = 780 ng/ml</li>
+
<li>maximum concentration: C<sub>max</sub> = 780 ng/mL</li>
-
<li>average concentration: C<sub>avg</sub> = 350 ng/ml</li><br><br>
+
<li>average concentration: C<sub>avg</sub> = 350 ng/mL</li></ul><br><br>
<p><center><img src="https://static.igem.org/mediawiki/2012/e/e1/Svn12_modelISHEMIC.png" style="height:800px; width:auto;"/></center></p>
<p><center><img src="https://static.igem.org/mediawiki/2012/e/e1/Svn12_modelISHEMIC.png" style="height:800px; width:auto;"/></center></p>
-
<center><b>Figure 6.</b> Multi-compartmental model for anakinra. For local therapy model, we added one additional compartment (shown in purple color) inside the heart compartment. This represents the affected area, where microencapsulated cells would be administered to. In conventional therapy model, the whole heart, both infarcted and undamaged tissue, are represented as one, since standard treatments can not target selected area specifically. </center><br><br>
+
<center><b>Figure 6.</b> Multi-compartment model for anakinra. For local therapy model, we added one additional compartment (shown in purple color) inside the heart compartment. This represents the affected area, where microencapsulated cells would be administered to. In conventional therapy model, the whole heart, both infarcted and undamaged tissue, are represented as one, since standard treatments can not target selected area specifically. </center><br><br>
-
<p>We calculated following partition coefficients according to a study on drug distribution in rodents.</p>
+
<p>We calculated following partition coefficients according to a study on drug distribution in rodents (Kim et al., 1995).</p>
Line 814: Line 888:
<tr class="normal">
<tr class="normal">
   <td class="normal">plasma</td>
   <td class="normal">plasma</td>
-
   <td class="normal">26%</td>
+
   <td class="normal">26 %</td>
</tr>
</tr>
</tbody>
</tbody>
Line 820: Line 894:
<img style="width:75%;" src="https://static.igem.org/mediawiki/2012/7/77/Doses_ISH_values1.png"></img>
<img style="width:75%;" src="https://static.igem.org/mediawiki/2012/7/77/Doses_ISH_values1.png"></img>
-
<center><b>Figure 7.</b> Dose distribution plot, fitted to literature data, which was used to determine the partition coefficients for the listed tissues.</center><br><br>
+
<center><b>Figure 7.</b> Dose distribution plot, fitted to the literature data, which was used to determine the partition coefficients for the listed tissues.</center><br><br>
-
<p>Partition coefficients of tissues are proportionate to drug concentrations (for non-eliminating tissue). From drug amount in each tissue and tissue size/weight data, we could find candidates for coefficient values. Parameters are further tuned to fit measured data as much as pssible.
+
<p>Partition coefficients of tissues are proportionate to drug concentrations (for non-eliminating tissue). From drug amount in each tissue and tissue volumes, we could find calculate concentration ratios and then estimate coefficient values. Parameters are further tuned to fit measured data as much as possible.
<table class="normal" style="font-size:90%; width:60%; text-align:left;">
<table class="normal" style="font-size:90%; width:60%; text-align:left;">
<tbody  class="normal">
<tbody  class="normal">
-
<thead class="normal"><td class="normal">Comparment</td>
+
<thead class="normal"><td class="normal">Compartment</td>
                                   <td class="normal">Parameter name</td>
                                   <td class="normal">Parameter name</td>
                                   <td class="normal">Value</td>  
                                   <td class="normal">Value</td>  
Line 875: Line 949:
<p><center><img style="width:75%;" src="https://static.igem.org/mediawiki/2012/e/ed/Anakinra_22sc.png"/></center></p><br>
<p><center><img style="width:75%;" src="https://static.igem.org/mediawiki/2012/e/ed/Anakinra_22sc.png"/></center></p><br>
-
 
+
<center><b>Figure 8. </b>Concentration distribution throughout tissues after subcutaneous administration. In this kind of administration, affected myocard cannot be targeted specifically; drug concentration is equal in both infarcted and undamaged heart tissue. With this kind of administration, the drug is also distributed throughout the whole body. Drug concentration is even higher in plasma or kidneys than concentration in heart.</center><br><br><br><br>
<!--p><center><img style="width:75%;" src="https://static.igem.org/mediawiki/2012/e/ed/Anakinra_22sc.png"/></center></p-->
<!--p><center><img style="width:75%;" src="https://static.igem.org/mediawiki/2012/e/ed/Anakinra_22sc.png"/></center></p-->
-
<p><center><img style="width:75%;" src="https://static.igem.org/mediawiki/2012/e/ea/Anakinra_30mcgNaDan.png"/></center></p><br><br>
+
<p><center><img style="width:75%;" src="https://static.igem.org/mediawiki/2012/e/ea/Anakinra_30mcgNaDan.png"/></center></p><br>
 +
<center><b>Figure 9. </b>Concentration distribution throughout tissues in local therapy. Here, the drug would be constantly produced in microencapsulated cells, which would be implanted into the damaged part of heart. Because heart muscle tissue is differently perfused than liver, very small amount of drug enters systemic circulation and localization of drug is even more evident.</center><br><br>
<p>We achieved significantly better results with our proposed therapy comparing with conventional treatment.</p><br>
<p>We achieved significantly better results with our proposed therapy comparing with conventional treatment.</p><br>
Line 894: Line 969:
<h2 style="color:grey;">References</h2>
<h2 style="color:grey;">References</h2>
<p style="color:grey;">
<p style="color:grey;">
-
Cawthorne, C., Prenant, C., Smigova, a, Julyan, P., Maroy, R., Herholz, K., Rothwell, N., et al. 2011. Biodistribution, pharmacokinetics and metabolism of interleukin-1 receptor antagonist (IL-1RA) using [18F]-IL1RA and PET imaging in rats. British journal of pharmacology 162,659–72.<br/><br/>
+
Bohoslawec, O., Trown, P. W., and Wills, R. J. (1986) Pharmacokinetics and tissue distribution of recombinant human alpha A, D, A/D(Bgl), and I interferons and mouse alpha-interferon in mice. <i>Journal of Interferon Research</i> <b>6</b>, 207-13.<br/><br/>
 +
 
 +
Cawthorne, C., Prenant, C., Smigova, a, Julyan, P., Maroy, R., Herholz, K., Rothwell, N., et al. (2011) Biodistribution, pharmacokinetics and metabolism of interleukin-1 receptor antagonist (IL-1RA) using [18F]-IL1RA and PET imaging in rats. <i>British journal of pharmacology</i> <b>162</b>, 659–72.<br/><br/>
 +
 
 +
Chatelut, E., Rostaing, L., Grégoire, N., Payen, J. L., Pujol, a, Izopet, J., Houin, G., et al. (1999) A pharmacokinetic model for alpha interferon administered subcutaneously. <i>British journal of clinical pharmacology</i> <b>47</b>, 365–71.<br/><br/>
 +
 
 +
Dahlén, E., Barchan, K., Herrlander, D., Höjman, P., Karlsson, M., Ljung, L., Andersson, M., et al. (2008) Development of Interleukin-1 Receptor Antagonist Mutants with Enhanced Antagonistic Activity In Vitro and Improved Therapeutic Efficacy in Collagen-Induced Arthritis. <i>J. Immunotoxicol.</i> <b>5</b>, 189-99.<br/><br/>
 +
 
 +
Hall, J. E. (2010) Guyton and Hall Textbook of Medical Physiology: with STUDENT CONSULT Online Access, 12e (Guyton Physiology) (p. 1120). Saunders.<br/><br/>
 +
 
 +
Johns, T. G., Kerry, J. A., Veitch, B. A. J., Johns, T. G., Kerry, J. A., Veitch, B. A. J., Mackay, I. R., et al. (1990) Pharmacokinetics, Tissue Distribution , and Cell Localization of [ 35 S ] Methionine-labeled Recombinant Human and Murine α Interferons in Mice. <i>Cancer Res.</i> <b>50</b>, 4718–23.<br/><br/>
 +
 
 +
Kim, D.C., Reitz, B., Carmichael, D. F., Bloedow, D. C. (1995) Kidney as a major clearance organ for recombinant human interleukin-1 receptor antagonist. <i>Journal of Pharmaceutical Sciences</i> <b>84</b>, 575–580.<br/><br/>
 +
 
 +
Lavoie, T. B., Kalie, E., Crisafulli-Cabatu, S., Abramovich, R., DiGioia, G., Moolchan, K., Pestka, S., et al. (2011) Binding and activity of all human alpha interferon subtypes. <i>Cytokine</i> <b>56</b>, 282–9.<br/><br/>
-
Chatelut, E., Rostaing, L., Grégoire, N., Payen, J. L., Pujol, a, Izopet, J., Houin, G., et al. 1999. A pharmacokinetic model for alpha interferon administered subcutaneously. British journal of clinical pharmacology 47,365–71.<br/><br/>
+
Levitt, D. G., Schoemaker, R. C. (2006) Human physiologically based pharmacokinetic model for ACE inhibitors: ramipril and ramiprilat. <i>BMC clinical pharmacology</i> <b>6</b>, 1.<br/><br/>
-
Johns, T. G., Kerry, J. A., Veitch, B. A. J., Johns, T. G., Kerry, J. A., Veitch, B. A. J., Mackay, I. R., et al. 1990. Pharmacokinetics, Tissue Distribution , and Cell Localization of [ 35 S ] Methionine-labeled Recombinant Human and Murine α Interferons in Mice 4718–4723.<br/><br/>
+
Marcellin, P., Boyer, N. (2003) Transition of care between paediatric and adult gastroenterology. Chronic viral hepatitis. <i>Best Pract Res Clin Gastroenterol.</i> <b>17</b>, 259-75. <br/><br/>
-
Kim, D.-C., Reitz, B., Carmichael, D. F., Bloedow, D. C. 1995. Kidney as a major clearance organ for recombinant human interleukin-1 receptor antagonist. Journal of Pharmaceutical Sciences 84,575–580.<br/><br/>
+
Mescher, A. (2009) Junqueira’s Basic Histology: Text and Atlas, 12th Edition (p. 480). McGraw-Hill Medical.<br/><br/>
-
Lavoie, T. B., Kalie, E., Crisafulli-Cabatu, S., Abramovich, R., DiGioia, G., Moolchan, K., Pestka, S., et al. 2011. Binding and activity of all human alpha interferon subtypes. Cytokine 56,282–9.<br/><br/>
+
Mrhar A. (1983) Računalniška simulacija nelinearnih farmakokinetičnih modelov pri razvoju optimalnih zdravil (Computer simulation of nonlinear pharmacokinetic models for development of optimal drugs), <i>Ph.D. Thesis</i><br/><br/>
-
Levitt, D. G., Schoemaker, R. C. 2006. Human physiologically based pharmacokinetic model for ACE inhibitors: ramipril and ramiprilat. BMC clinical pharmacology 6,1.<br/><br/>
+
Roferon-A (Interferon alfa-2a, recombinant), Hoffmann-La Roche Inc.<br/><br/>
-
Marcellin, P., Boyer, N. 2003. Chronic viral hepatitis 17,259–275.<br/><br/>
+
Roy, A., Georgopoulos, P. G. (1997) Mechanistic Modeling of Transport and Metabolism in Physiological Systems Towards a Framework for an Exposure and Dose Modeling and Analysis System.<br/><br/>
-
ROFERON -A. n.d.<br/><br/>
+
Thitinan, S., McConville, J. T. (2009) Interferon alpha delivery systems for the treatment of hepatitis C. <i>International journal of pharmaceutics</i> <b>369</b>, 121–35.<br/><br/>
-
Roy, A., Georgopoulos, P. G. 1997. Mechanistic Modeling of Transport and Metabolism in Physiological Systems Towards a Framework for an Exposure and Dose Modeling and Analysis System.<br/><br/>
+
Yang, B.B., Baughman, S., Sullivan, J. T. (2003) Pharmacokinetics of anakinra in subjects with different levels of renal function. <i>Clinical pharmacology and therapeutics</i> <b>74</b>, 85–94.
-
Thitinan, S., McConville, J. T. 2009. Interferon alpha delivery systems for the treatment of hepatitis C. International journal of pharmaceutics 369,121–35.
 
-
Yang, B.-B., Baughman, S., Sullivan, J. T. 2003. Pharmacokinetics of anakinra in subjects with different levels of renal function. Clinical pharmacology and therapeutics 74,85–94.
 

Latest revision as of 03:46, 27 October 2012


Pharmacokinetic model

What is pharmacokinetic model?
Pharmacokinetic model is a quantitative description of drug absorption, distribution, metabolism and elimination from the body.


Why do we need it?

  • To compare conventional and our therapy.
  • To calculate required drug production at target site.


How was our pharmacokinetic model constructed?

  • The model is physiologically based.
  • Organs and tissues are presented as compartments.
  • Processes in the body are described with ordinary differential equations (ODEs).
  • The model is customized for chosen drug and its type of administration.


For which drugs (diseases) did we test our model?

  • interferon (hepatitis C)
  • anakinra (ischemic heart disease)


3D animation
3D animation serves as a real-time visual representation of how drug concentrations change throughout the course of the therapy in various compartments. The animation was implemented in C++ programming language and OpenGL was used for 3D rendering.



Results

  • Our therapy does not produce concentration fluctuations - absence of concentration peaks and lows.
  • With local drug production, drug concentration rises to a desired level within a day (for both tested drugs).
  • For each drug, we obtained an estimation of the degree of drug localization.


What did the results tell us?
Results show us the drug is more localized; the degree of localization depends on characteristics of tissue (e.g. perfusion) where microencapsulated cells would be applied. Additionally, absence of high concentration peaks suggests that our therapy would have less systemic side effects than conventional methods. Since the drug would be constantly produced, drug concentration would not drop enough for the drug to be therapeutically ineffective. According to these findings, steady drug levels over time therefore imply an optimized treatment. With known drug production rate of microencapsulated cells, it is possible to adjust the dosage to reach the desired drug concentration in the target tissue. Therapeutic production was estimated according to activity assays with mammalian cells. After obtaining the results of pharmacokinetic simulation, we were able to calculate the number of microencapsulated cells needed for the therapy.

Introduction

The major downsides of standard interferon treatment are substantial side effects that to a large degree are the consequence of very high drug concentrations, occurring shortly after drug administration. If high concentration peaks could be avoided and lower levels maintained steadily over time, this would result in reduced side effects without compromising therapeutic effectiveness. This idea is currently being tested in clinical trials using the interferon infusion pump (COPE-HCV: Phase 2, randomized, open-label, active-control, dose-ranging study of interferon alfa-2b given via continuous sub-Q infusion; trial by Medronic Inc.).

One of our goals is to show that therapy with drug-producing cells is more beneficial than standard treatments based on drug injections, which are in use today. We predict that if the drug is constantly produced inside the body, it could reach a steady concentration at almost any desired level. We believe localization of therapeutic cells would also decrease the proportionate drug concentrations in non-target tissues, thus further reducing the side effects.

We tested this hypothesis and compared standard therapies with our proposed treatment.
Because of complex physiological mechanisms and an extensive set of biological parameters that cannot be accurately measured for either ethical or technical reasons, we developed a model that covers the most crucial aspects and processes. At the same time, computer simulation provides simpler and faster option than in vivo research.

3D visualization

We implemented a real-time 3D visualization of the pharmacokinetic simulation for hepatitis C and ischemia. This way, we were able to visually represent how drug concentrations change throughout the course of the therapy in various compartments, such as liver. In our 3D visualizations, drug concentration levels are represented as red color intensities - the compartments with high concentrations are closer to saturated red color, while the compartments with low concentrations are closer to white color. Red color intensity of a compartment was calculated as a ratio of the compartment's drug concentration at a given time to the maximal reached concentration of any compartment during the course of the therapy, or therapies compared.

For implementation, we used C++ programming language and OpenGL for 3D rendering. The 3D anatomical models, obtained from BodyParts3D database, were optimized for real-time rendering using MeshLab and rendered as vertex arrays in OpenGL. Compartment concentration values as a function of time were obtained as vectors from the pharmacokinetic simulation and integrated into the application.


A physiologically based pharmacokinetic model

A pharmacokinetic model is a quantitative description of drug absorption, distribution, metabolism and elimination from the body. A model is defined by a system of ordinary differential equations to represent essential drug kinetics.
We took into account the necessary biological parameters and based the simulated processes on actual physiological mechanisms to construct a physiologically plausible model.

Selecting the optimal physiological model

A physiologically based model is composed of multiple compartments which represent organs of the body. Parts were chosen in accordance with drug and tissue specifics, so that the relevant organs are represented as separate compartments, while other tissues were merged on the basis of common characteristics (Levitt et al., 2006).

Determing compartments

We used the following characteristics as criteria for splitting and merging organs:

  • Liver is the target organ of therapy.
  • Interferon alpha is widely distributed into body tissues – highest concentrations occur in kidney, liver and lung.
  • Interferon is a water-soluble molecule; it is only poorly distributed in adipose tissue.
  • Interferon alpha does not cross the blood-brain barrier.
  • Skin and muscle tissues do not seem to have much higher concentrations of the drug in comparison to adipose tissue.
  • Skin, musle and adipose tissue have similar, slow blood perfusion.
  • Gut, spleen and heart are all rapidly perfused tissues.
  • Interferon is mainly eliminated via renal catabolism, while hepatic metabolism accounts only for a minor pathway of elimination.

We decided to define separate compartments for the liver, kidney and lungs. All other rapidly perfused tissues are grouped together as one part. Since venous blood enters the lungs and arterial blood flows into all other organs, we separately simulate venous and arterial blood. Because interferon does not cross the blood-brain barrier it is not necessary for the brain to be modeled separately. Skin, muscle, fat and other slowly perfused tissues are merged together into one compartment.

We constructed three final models - two for standard interferon treatments and a third for a prospective therapy with drug-producing microencapsulated cells. The fundamental design is the same in all models; they are only modified for specific entry points of the drug and the corresponding absorption or production processes. On the diagram, blocks representing different administrations are shown in distinct colors: blue for the intravenous bolus, green for the subcutaneous injection and red for the interferon production by microencapsulated cells that are implanted into liver (Figure 1).




Figure 1. Multi-compartment physiologically based pharmacokinetic model, customized for different Interferon therapies. Main difference is the route of administration - for each administration, there is different compartment from which drug enters the body. Blue blocks are used in modeling for intravenous administration, green blocks are specific for subcutaneous administration and red blocks represent continuous production of the drug, in this case in the liver compartment. Therapies of course don't differ only in drug entry points, but also in absorption or production processes, percent of drug that enters the blood stream and so on.



Parameters

Types of parameters used:

1.) Species specific parameters

Qi – blood flows to tissues
Vi – organ volumes

Blood flows
Cardiac output Qc Qc = 5.58 L/min = 335L/h
Compartment Parameter name Percent cardiac output
lung Qlung 100% Qc
venous blood Qven 100% Qc
arterial blood Qart 100% Qc
liver Qliver 25% Qc
kidney Qkidney 19% Qc
rapidly perfused tissue Qrpt 18% Qc
slowly perfused tissue Qspt 38% Qc

Sum of blood flows through liver, kidney, rapidly and slowly perfused tissue must always equal total cardiac output.

Tissue Volumes
Body Weight BW BW = 70 kg => *BV = 70 L
Compartment Parameter name Percent body volume
lung Vlung 0.8% BV
venous blood Vven 5.57% BV
arterial blood Vart 2.43% BV
liver Vliver 2.60% BV
kidney Vkidney 0.44% BV
rapidly perfused tissue Vrpt 5.16% BV
slowly perfused tissue Vspt 83% BV

* average body density = 1 kg /L

2.) Individual specific parameters

BW - body weight
Qc - cardiac output
varying percents of tissue volumes (e.g. percent body fat)

The variability of parameters

Parameter values can range significantly between individuals, depending on factors such as age, sex, renal function, activity level and diet.
For instance, cardiac output can vary significantly even in one individual, depending on the current activity (sleeping, sitting, running etc.). There can be substantial differences in, for example, the percent of body fat (accounted for in slowly perfused tissue) comparing individuals with an otherwise similar profile (same age, sex, etc.).
Values used in our model present an average adult male, weighing 70 kg, with mean cardiac output and normal renal function.

3.) Chemical specific parameters

t1/2 - drug half life
kel - elimination rate (kidney)
Pblood:tissue - partition coefficient

Interferon is expected to be found only in the plasma and not in red blood cells; therefore we can conclude that the amount of interferon found in plasma is equal to the amount in blood. Since blood is comprised of four parts plasma and three parts red blood cells, we can calculate the blood to plasma ratio to be:

The equation has the same form for each kind of tissue:

Uncertainty of partition coefficients

A critical element in human physiologically based pharmacokinetic modeling is the uncertainty of values of the partition coefficients. Partition coefficients are an important aspect of pharmacokinetic modeling, because they denote how the drug distributes throughout body tissues. The value of each coefficient has a complex dependence on solubility, permeability, pH, binding affinity of the drug to the receptor receptor and many other factors. The difference in a few amino acids between subtypes of interferon alpha can impact values of coefficients quite noticeably. Even so, there are deviations of evaluated drug distribution of the same subtype of interferon alpha, depending on the detection method (radioactivity, ELISA).

For legal and ethical reasons these values cannot be directly measured in humans. We had to rely on various studies of interferon tissue distribution in rodents to calculate partition coefficients (Johns et al, 1990). It is generally assumed that animal and human partition coefficients are similar for the same kind of tissue.

Chemical and ROA* specific parameters

* ROA - route of administration

D – dose
k0 – absorption rate constant (zero-order process)
ka – absorption rate constant (first-order process)
F – bioavailability (percent of absorbed dose)
kprod – production rate constant

General mass-balance equations

The equation below describes the change in concentration over time in non-eliminating tissues. The equation has the same form for both rapidly and slowly perfused tissue. Each compartment is then described with its distinctive values for blood flow, concentration and partition coefficient.


We can use the same form for the liver compartment, since the metabolism of interferon is negligible.


Although the lungs also represent a non-eliminating tissue in this model, the equation is slightly different, since venous and not arterial blood flows into the tissue.


The kidneys represent a site of elimination, so we had to include this process as well.


In the equation describing the change of concentration in venous blood, there is a sum of blood flows which flows in from multiple compartments. These include liver, kidney, rapidly and slowly perfused tissue.


Blood flow from the lungs is accounted for in the equation for the arterial blood compartment.


Modeling pharmacokinetic processes

1.) Absorption

Absorption processes depend on the route of administration:

a.) Intravenous administration

The drug is injected directly into the blood stream, so there is no special absorption process. The dose is bolus and enters the system completely.


b.) Subcutaneous administration

In subcutaneous administration, more than 80% of the initial dose enters systemic circulation. This percentage is defined as bioavailability and can reach up to 95% for interferon alpha given subcutaneously.
Time to peak Tmax
The absorption half-life is approximately 2.3 h

Relying on the available studies, we specified the absorption as a two phase process:

1.) At the beginning, when the drug concentration at the injection site is highest, only a limited amount of drug can be absorbed into the blood - the rate of absorption is maximal. This is described by a zero order process from t=0 to t=tk.
Zero order process equation:


2.) At t=tk the concentration drops to a level where rate of absorption is proportional to local concentration. Absorption follows first order kinetics from t= tk until the bioavailable drug is completely absorbed.
First order process equation:


Constant ka was taken from literature and the value of k0 was calculated using parameters from the same source. Fz denotes the portion of the drug that is absorbed by zero order process.


2.) Distribution

Estimating partition coefficients

Distribution process is usually modeled indirectly with the use of partition coefficients, which are drug and tissue specific. We derived partition coefficients according to the tissue distribution data in available literature (Thitinian et al., 2009; Johns et al., 1990). However, these values do differ from study to study (e.g. that kidney concentration is 1-9x times that of plasma (Bohoslawec et al., 1986), but they overall agree that after subcutaneous or intravenous administration, interferon concentration is highest in the kidney, very high in the liver and lungs and to some extent in other rapidly perfused tissues, but lower in the adipose and other slowly perfused tissues.


Compartment Parameter name Value
lung Pblood:lung 2.5
liver Pblood:liver 2.8
kidney Pblood:kidney 8.5
rapidly perfused tissue Pblood:rpt 2.2
slowly perfused tissue Pblood:spt 0.7

3.) Elimination

In the literature the terminal half-life for interferon alpha is from 3-8 h, with a mean around 5 h. Interferon alpha is detectable in the plasma for 4-8 h after the rapid intravenous injection and for 16-30h after subcutaneous administration.
The elimination constant for the kidneys was determined by curve fitting to data from the obtained studies.

Simulation data for different therapies

Interferon alpha subcutaneous therapy modeling

We simulated common type of interferon administration:

  • 3 times per week, subcutaneous administration
  • 3 MU = 11100 ng, dose
  • F = 93%, bioavailability
  • Ko, Ka - zero and first order absorption
  • t1/2 = 5h, half life


We constructed the model according to data from the literature (Roferon A, Hoffman-La Roche; Chatelut, 1999; Thitinian et al., 2009; Marcellin et al., 2003). Results are shown in graph below (Figure 2).


Figure 2. Standard subcutaneous therapy: high peaks right after administration, followed by rapid concentration decrease.


Calculating desired drug production

Our model assumes no cell divisions and therefore a constant therapeutic cell count in the body. In this case, cell production of a drug is mathematically equivalent to a continuous infusion of a drug into the tissue, where cells are located. To calculate the desired concentration level in the target tissue, we calculated AUC (area under the curve) of liver concentrations for subcutaneous administration over a period of one week. This allowed us to estimate the average concentration in the liver which corresponds to the therapeutic level we wish to reach and maintain with cellular drug production. The combination of kidney elimination and constant production leads to stable drug levels in tissues. After about four half-lives, drug concentrations rise to final levels and remain steady as long as the production rate stays the same (Figure 3).



Figure 3. Drug distribution troughout tissues with Switch-IT therapy.


Results and conclusions

If we compare result of different treatments it clearly shows that local therapy with drug producing cells does not result in concentration fluctuations, which are present in common therapies.
Concentrations never reach levels as high as in subcutaneous or intravenous administration. On the other hand, concentrations do not drop to a point, which would prevent the drug from becoming therapeutically inefficient.

When the drug is produced at the target site, concentration levels between target tissues and other organs are inclined in favor of the target tissue - although kidneys show a substantial uptake of the drug in subcutaneous administration, the relative difference between kidney and liver concentrations are smaller with localized therapy (Figure 4).



Figure 4. Comparison of subcutaneous and Switch-IT therapy. We can see the new therapy results in stable concentration levels, without peaks and lows.


From model to implementation

When we calculate the desired steady concentration in the target tissue, it was possible to estimate how many cells would need to be implanted. From the required production rate, which is thought to be linearly proportional to the cell count, we can calculate the total number of cells needed. Considering the average number of cells per microcapsule and mean capsule size, we can calculate the amount of microcapsules we would need to inject in order to attain the wanted therapeutic effect.

Anakinra therapy modeling

For modeling the local application of anakinra-producing microencapsulated cells in the ischaemic muscle tissue we adopted slightly different approach. Liver is perfused by sinusoid capillaries which allow for passage of macromolecules to the tissue, therefore even protein therapeutics can be rapidly equilibrated with blood and distributed around the organism. Muscle tissue, however, is perfused with continuous capillaries, which prevents fast passage of macromolecules from the tissue to the blood (Figure 5). Protein molecules pass the continuous capillaries only through slow trans-pinocytosis. Consequently this increases the local concentration of therapeutics permitting the achievement of high local concentration with relatively little loss to the system and elimination by kidneys.


Figure 5. Different capillary types in muscle tissue (continuous capillary, left) and in liver (sinusoid capillary, right). Image from: (Mescher, 2009)


Therefore we took this feature of the tissue into account when calculating partition coefficients between tissue and blood compartment. Myoglobin has similar molecular weight as anakinra; therefore we assumed its partition coefficient to be the same(Hall, 2010). We calculated the therapeutic concentration of anakinra to be 25 ng/mL, which we derived from the anakinra IC_50 value (Dahlén et al., 2008).

Compartmental model and simulation parameters

According to a study of anakinra distribution and elimination after subcutaneous injection, we calculated relevant compartments for physiological model. At first, we modeled standard regime for anakinra therapy: 100 mg per day subcutaneously with the following parameters:

  • bioavailability: F = 95%
  • half life: T1/2 = 5h
  • time to maximum concentration: Tmax = 4h
  • maximum concentration: Cmax = 780 ng/mL
  • average concentration: Cavg = 350 ng/mL


Figure 6. Multi-compartment model for anakinra. For local therapy model, we added one additional compartment (shown in purple color) inside the heart compartment. This represents the affected area, where microencapsulated cells would be administered to. In conventional therapy model, the whole heart, both infarcted and undamaged tissue, are represented as one, since standard treatments can not target selected area specifically.


We calculated following partition coefficients according to a study on drug distribution in rodents (Kim et al., 1995).

tissuepercent of dose
kidney 11 %
liver 9 %
gut 6 %
lung 2 %
muscle 43 %
plasma 26 %
Figure 7. Dose distribution plot, fitted to the literature data, which was used to determine the partition coefficients for the listed tissues.


Partition coefficients of tissues are proportionate to drug concentrations (for non-eliminating tissue). From drug amount in each tissue and tissue volumes, we could find calculate concentration ratios and then estimate coefficient values. Parameters are further tuned to fit measured data as much as possible.

Compartment Parameter name Value
lung Pblood:lung 0.4
liver Pblood:liver 0.8
kidney Pblood:kidney 9
rapidly perfused tissue Pblood:rpt 2.2
slowly perfused tissue Pblood:spt 0.23
muscle Pblood:muscle 0.25
heart Pblood:heart 0.2
target tissue Pblood:target 33



Figure 8. Concentration distribution throughout tissues after subcutaneous administration. In this kind of administration, affected myocard cannot be targeted specifically; drug concentration is equal in both infarcted and undamaged heart tissue. With this kind of administration, the drug is also distributed throughout the whole body. Drug concentration is even higher in plasma or kidneys than concentration in heart.





Figure 9. Concentration distribution throughout tissues in local therapy. Here, the drug would be constantly produced in microencapsulated cells, which would be implanted into the damaged part of heart. Because heart muscle tissue is differently perfused than liver, very small amount of drug enters systemic circulation and localization of drug is even more evident.


We achieved significantly better results with our proposed therapy comparing with conventional treatment.


References

Bohoslawec, O., Trown, P. W., and Wills, R. J. (1986) Pharmacokinetics and tissue distribution of recombinant human alpha A, D, A/D(Bgl), and I interferons and mouse alpha-interferon in mice. Journal of Interferon Research 6, 207-13.

Cawthorne, C., Prenant, C., Smigova, a, Julyan, P., Maroy, R., Herholz, K., Rothwell, N., et al. (2011) Biodistribution, pharmacokinetics and metabolism of interleukin-1 receptor antagonist (IL-1RA) using [18F]-IL1RA and PET imaging in rats. British journal of pharmacology 162, 659–72.

Chatelut, E., Rostaing, L., Grégoire, N., Payen, J. L., Pujol, a, Izopet, J., Houin, G., et al. (1999) A pharmacokinetic model for alpha interferon administered subcutaneously. British journal of clinical pharmacology 47, 365–71.

Dahlén, E., Barchan, K., Herrlander, D., Höjman, P., Karlsson, M., Ljung, L., Andersson, M., et al. (2008) Development of Interleukin-1 Receptor Antagonist Mutants with Enhanced Antagonistic Activity In Vitro and Improved Therapeutic Efficacy in Collagen-Induced Arthritis. J. Immunotoxicol. 5, 189-99.

Hall, J. E. (2010) Guyton and Hall Textbook of Medical Physiology: with STUDENT CONSULT Online Access, 12e (Guyton Physiology) (p. 1120). Saunders.

Johns, T. G., Kerry, J. A., Veitch, B. A. J., Johns, T. G., Kerry, J. A., Veitch, B. A. J., Mackay, I. R., et al. (1990) Pharmacokinetics, Tissue Distribution , and Cell Localization of [ 35 S ] Methionine-labeled Recombinant Human and Murine α Interferons in Mice. Cancer Res. 50, 4718–23.

Kim, D.C., Reitz, B., Carmichael, D. F., Bloedow, D. C. (1995) Kidney as a major clearance organ for recombinant human interleukin-1 receptor antagonist. Journal of Pharmaceutical Sciences 84, 575–580.

Lavoie, T. B., Kalie, E., Crisafulli-Cabatu, S., Abramovich, R., DiGioia, G., Moolchan, K., Pestka, S., et al. (2011) Binding and activity of all human alpha interferon subtypes. Cytokine 56, 282–9.

Levitt, D. G., Schoemaker, R. C. (2006) Human physiologically based pharmacokinetic model for ACE inhibitors: ramipril and ramiprilat. BMC clinical pharmacology 6, 1.

Marcellin, P., Boyer, N. (2003) Transition of care between paediatric and adult gastroenterology. Chronic viral hepatitis. Best Pract Res Clin Gastroenterol. 17, 259-75.

Mescher, A. (2009) Junqueira’s Basic Histology: Text and Atlas, 12th Edition (p. 480). McGraw-Hill Medical.

Mrhar A. (1983) Računalniška simulacija nelinearnih farmakokinetičnih modelov pri razvoju optimalnih zdravil (Computer simulation of nonlinear pharmacokinetic models for development of optimal drugs), Ph.D. Thesis

Roferon-A (Interferon alfa-2a, recombinant), Hoffmann-La Roche Inc.

Roy, A., Georgopoulos, P. G. (1997) Mechanistic Modeling of Transport and Metabolism in Physiological Systems Towards a Framework for an Exposure and Dose Modeling and Analysis System.

Thitinan, S., McConville, J. T. (2009) Interferon alpha delivery systems for the treatment of hepatitis C. International journal of pharmaceutics 369, 121–35.

Yang, B.B., Baughman, S., Sullivan, J. T. (2003) Pharmacokinetics of anakinra in subjects with different levels of renal function. Clinical pharmacology and therapeutics 74, 85–94.


Next: Modeling methods >>