cammey e. cole meredith college march 11, 2005
DESCRIPTION
Some Deterministic Models in Mathematical Biology: Physiologically Based Pharmacokinetic Models for Toxic Chemicals. Cammey E. Cole Meredith College March 11, 2005. Outline. Introduction to compartment models Research examples Linear model Analytics Graphics Nonlinear model - PowerPoint PPT PresentationTRANSCRIPT
![Page 1: Cammey E. Cole Meredith College March 11, 2005](https://reader035.vdocument.in/reader035/viewer/2022062802/56814565550346895db2365f/html5/thumbnails/1.jpg)
1
Some Deterministic Models in Mathematical Biology: Physiologically Based
Pharmacokinetic Models for Toxic Chemicals
Cammey E. Cole
Meredith College
March 11, 2005
![Page 2: Cammey E. Cole Meredith College March 11, 2005](https://reader035.vdocument.in/reader035/viewer/2022062802/56814565550346895db2365f/html5/thumbnails/2.jpg)
2
Outline• Introduction to compartment models• Research examples• Linear model
– Analytics– Graphics
• Nonlinear model• Exploration
![Page 3: Cammey E. Cole Meredith College March 11, 2005](https://reader035.vdocument.in/reader035/viewer/2022062802/56814565550346895db2365f/html5/thumbnails/3.jpg)
3
Physiologically Based Pharmacokinetic (PBPK)
ModelsA physiologically based
pharmacokinetic (PBPK) model for the uptake and elimination of a chemical in rodents is developed to relate the amount of IV and orally administered chemical to the tissue doses of the chemical and its metabolite.
![Page 4: Cammey E. Cole Meredith College March 11, 2005](https://reader035.vdocument.in/reader035/viewer/2022062802/56814565550346895db2365f/html5/thumbnails/4.jpg)
4
Characteristics of PBPK Models
• Compartments are to represent the amount or concentration of the chemical in a particular tissue.
• Model incorporates known tissue volumes and blood flow rates; this allows us to use the same model across multiple species.
• Similar tissues are grouped together.• Compartments are assumed to be well-mixed.
![Page 5: Cammey E. Cole Meredith College March 11, 2005](https://reader035.vdocument.in/reader035/viewer/2022062802/56814565550346895db2365f/html5/thumbnails/5.jpg)
5
Example of Compartment in PBPK Model
• QK is the blood flow into the kidney.
• CVK is the concentration of drug in the venous blood leaving the kidney.
QK CVKKidney
![Page 6: Cammey E. Cole Meredith College March 11, 2005](https://reader035.vdocument.in/reader035/viewer/2022062802/56814565550346895db2365f/html5/thumbnails/6.jpg)
6
Example of Compartment in PBPK Model
• CK is the concentration of drug in the kidney at time t.
• CVK is the concentration of drug in the venous blood leaving the kidney at time t.
• QK is the blood flow into the kidney.
• VK is the volume of the kidney.
K
KKKK
V
CVCQ
dt
dC
![Page 7: Cammey E. Cole Meredith College March 11, 2005](https://reader035.vdocument.in/reader035/viewer/2022062802/56814565550346895db2365f/html5/thumbnails/7.jpg)
7
BenzeneAveloar Space
Lung Blood
Fat
Slowly Perfused
Rapidly Perfused
Liver
StomachMuconic
Acid
ExhaledInhaled
VenousBlood
CV
ArterialBlood
CA
SlowlyPerfused
Blood
Liver
Blood
Liver
Blood
Liver
BOBO
BO
BO
PH
PH
HQ
HQ
BZ
BZ
BZ
BZ
BZ
ConjugatesPH Catechol
THB
ConjugatesHQ
PMA
Benzene
Benzene Oxide Phenol Hydroquinone
RapidlyPerfusedFat
BO PH
SlowlyPerfused
PH
Fat
PH
RapidlyPerfused Fat
SlowlyPerfused
RapidlyPerfused
HQ HQ HQ
![Page 8: Cammey E. Cole Meredith College March 11, 2005](https://reader035.vdocument.in/reader035/viewer/2022062802/56814565550346895db2365f/html5/thumbnails/8.jpg)
8
Benzene Plot
![Page 9: Cammey E. Cole Meredith College March 11, 2005](https://reader035.vdocument.in/reader035/viewer/2022062802/56814565550346895db2365f/html5/thumbnails/9.jpg)
9
Benzene Plot
![Page 10: Cammey E. Cole Meredith College March 11, 2005](https://reader035.vdocument.in/reader035/viewer/2022062802/56814565550346895db2365f/html5/thumbnails/10.jpg)
10
4-Methylimidazole (4-MI)
Metabolism
Adipose
Other Tissues
Kidney
Urine
Liver
B
l
o
o
d
Initial condition:IV exposure
Other Tissues
Leaving System
Blood
Stomach
Initial condition: gavage exposure
Feces
4-MI
Metabolite
![Page 11: Cammey E. Cole Meredith College March 11, 2005](https://reader035.vdocument.in/reader035/viewer/2022062802/56814565550346895db2365f/html5/thumbnails/11.jpg)
11
4-MI Female Rat Data (NTP TK)
![Page 12: Cammey E. Cole Meredith College March 11, 2005](https://reader035.vdocument.in/reader035/viewer/2022062802/56814565550346895db2365f/html5/thumbnails/12.jpg)
12
Linear Model ExampleA drug or chemical enters the body via the stomach.
Where does it go?Assume we can think about the body as three
compartments:– Stomach (where drug enters)– Liver (where drug is metabolized)– All other tissues
Assume that once the drug leaves the stomach, it can not return to the stomach.
![Page 13: Cammey E. Cole Meredith College March 11, 2005](https://reader035.vdocument.in/reader035/viewer/2022062802/56814565550346895db2365f/html5/thumbnails/13.jpg)
13
Schematic of Linear ModelStomach
Other Tissues
Liver
• x1, x2 , and x3 represent amounts of the drug.• a, b, and c represent flow rates.
x2
x1
x3
a
b c
![Page 14: Cammey E. Cole Meredith College March 11, 2005](https://reader035.vdocument.in/reader035/viewer/2022062802/56814565550346895db2365f/html5/thumbnails/14.jpg)
14
Linear Model EquationsLet’s look at the change of amounts in each
compartment, assuming the mass balance principle is applied.
323
3212
11
cxbxdt
dx
cxbxaxdt
dx
axdt
dx
![Page 15: Cammey E. Cole Meredith College March 11, 2005](https://reader035.vdocument.in/reader035/viewer/2022062802/56814565550346895db2365f/html5/thumbnails/15.jpg)
15
Linear Model (continued)Let’s now write the system in matrix form.
3
2
1
3
2
1
0
00
x
x
x
cb
cba
a
x
x
x
dt
d
![Page 16: Cammey E. Cole Meredith College March 11, 2005](https://reader035.vdocument.in/reader035/viewer/2022062802/56814565550346895db2365f/html5/thumbnails/16.jpg)
16
Linear Model (continued)
• Find the eigenvectors and eigenvalues.
• Write general solution of the differential equation.
• Use initial conditions of the system to determine particular solution.
![Page 17: Cammey E. Cole Meredith College March 11, 2005](https://reader035.vdocument.in/reader035/viewer/2022062802/56814565550346895db2365f/html5/thumbnails/17.jpg)
17
Linear Model (continued)For our example, let a=3, b=4, and c=11. Then, our general
solution would be given by:
tt ekekk
x
x
x15
33
21
3
2
1
1
1
0
1
2
3
4
11
0
![Page 18: Cammey E. Cole Meredith College March 11, 2005](https://reader035.vdocument.in/reader035/viewer/2022062802/56814565550346895db2365f/html5/thumbnails/18.jpg)
18
Initial ConditionsUsing the initial conditions of
we are representing the fact that the drug began in the stomach and there were no background levels of the drug in the system.
,0)0(
0)0(
9)0(
3
2
1
x
x
x
![Page 19: Cammey E. Cole Meredith College March 11, 2005](https://reader035.vdocument.in/reader035/viewer/2022062802/56814565550346895db2365f/html5/thumbnails/19.jpg)
19
Graphical Results
![Page 20: Cammey E. Cole Meredith College March 11, 2005](https://reader035.vdocument.in/reader035/viewer/2022062802/56814565550346895db2365f/html5/thumbnails/20.jpg)
20
Schematic of Nonlinear Model
x1, x2 , x3, and x4 represent amounts of the drug.
Stomach
Other Tissues
Liverx2
x1
x3
a
b c Metabolitex4
2
2
xK
Vx
![Page 21: Cammey E. Cole Meredith College March 11, 2005](https://reader035.vdocument.in/reader035/viewer/2022062802/56814565550346895db2365f/html5/thumbnails/21.jpg)
21
Nonlinear Model Equations
2
24
323
2
2321
2
11
xK
xV
dt
dx
cxbxdt
dx
xK
xVcxbxax
dt
dx
axdt
dx
![Page 22: Cammey E. Cole Meredith College March 11, 2005](https://reader035.vdocument.in/reader035/viewer/2022062802/56814565550346895db2365f/html5/thumbnails/22.jpg)
22
Nonlinear Model a=0.2, b=0.4, c=0.1, V=0.3, K=4
![Page 23: Cammey E. Cole Meredith College March 11, 2005](https://reader035.vdocument.in/reader035/viewer/2022062802/56814565550346895db2365f/html5/thumbnails/23.jpg)
23
Exploration
• What would happen if the parameters were changed?
• What would happen if the initial conditions were changed?
We will now use Phaser to explore these questions. Phaser files are on the website:
www.meredith.edu/math/faculty/cole/maaworkshop/pbpkmodels.htm