280 likes | 487 Views
Human Body Drug Simulation. Nathan Liles Benjamin Munda. Presentation Outline. Objective Background Model Overview Organ and Body Theory Case Studies. Objective. Pharmacokinetics : Seeks to determine fate of substances administered externally to a living organism.
E N D
Human Body Drug Simulation Nathan Liles Benjamin Munda
Presentation Outline • Objective • Background • Model Overview • Organ and Body Theory • Case Studies
Objective Pharmacokinetics: Seeks to determine fate of substances administered externally to a living organism 1) Overall PBPK ADME model a) How to divide the body and describe movement? b) How to make mathematical model consistent with physiology/anatomy of the human body 2) Specific Component of a PBPK model a) How does a drug’s structure and charge affect its movement in capillaries? b) How does an enzyme within the liver interact with a drug? Are any intermediates created? Our project focuses on the first objective above: creating an overall PBPK model that is accurate, physiologically correct, and innovative Our focus, Oral Administration a) most common, practical b) model can be expanded
Background: Potential Applications On Pharmaceuticals every year: Companies R&D ≈ $70 billion b) Consumers spend≈ $200 billion • A Mathematical Body Model Could: • For Companies: • Accelerate R&D • For consumers: • Help doctors optimize doses • Lower costs of prescription drugs
Background: What Has Been Done Previous Work: Compartment models Well-mixed regions Few parameters <15 Few equations < 25 Our Work Began with current models of key organs 23 Tissues included, avoided well-mixed assumption Incorporated the tissues into a circulatory system with mass transfer on the capillary level ≈1,100 Simultaneous ODEs ≈200 constants ≈100,000 steps
Overview of Our Proposed Model • Path of drug • Enters stomach • Moves through small intestine • Enters the blood, first pass through the liver • Liver to the heart • Heart to the entire body • Interacts with the body • a) Reacts at intended site • b) Eliminated in liver, kidney • Returns to the Heart 1 2 3 4 6b 5 6b 7 6a
Absorption: Stomach • Stomach information: • 1) “Churning” creates a well mixed volume • 2) Exit flow rate of mass depends on the mass inside the stomach • 3) Little absorption of the drug into the bloodstream occurs in the stomach Governing Equations: Concentration: Dose (mg) enters stomach, which has some mass inside The drug exits with a semi-constant concentration and a flow rate that varies with: a) L or S? b) Mass c) Liberation Flow rate: 1) For liquids 2) For solids
Absorbtion: Small Intestine • Small Intestine Information: • Main site of drug absorption • ≈7 meters long with an average diameter of 2.5-3 cm • 3) Modeled as a PFR Governing Equation • Assumptions a) Radial variations are unimportant b) Diffusive flux term is negligible Final Form: • c) Superficial velocity is a variable
Absorbtion: Small Intestine A LaPlace transform was performed in the z-dimension. The equation became a linear homogenous ODE in the time dimension. This integral cannot be solved analytically, thus the inverse Laplace transform cannot be solved.
Numerical Methods Method of Lines 1) Discretize Space Variable Z 2) System of Equations to Solve
Numerical Methods 4th Order Runge-Kutta Method of lines means we have many ordinary differential equations=4RK Used to integrate a function: Where values for k (slope estimates) come from: Described by a 1st Order ODE: Given initial values for y estimates next y in time by:
Absorbtion: Small Intestine Amount absorbed through the gut wall for each time interval
Distribution • Estimations used in distribution: • 1) Time from SI to liver negligible • 2) Time from Inferior vena cava to heart negligible. • 3) 8 seconds to get from the right heart to the left heart • 4) 3 Seconds to get to the extremities. • 5) 3 seconds back to the heart 1 2 3 4 5
Distribution: Blood Concentration Blood concentration is the ‘heart’ of our model Heart chosen as site to track Governing Equations Heart Information: Well mixed tank Receives influx from entire body and outputs to the entire body 1) 2)
Distribution: Capillaries Governing Equations Capillary Information: Literature ≈ 40 billion capillaries Calculated (A, d, l) ≈ 20 billion capillaries 2) Number of capillaries in an organ estimated by percent of total body blood flow Example: Brain and kidney are small, but receive a large (≈30-35%) amount of blood, requires dense capillaries 1) 2)
Metabolism: The Liver Macrostructure Microstructure Governing Equations • Liver Information • Modeled as PFR Simultaneous with tissue • Blood Mixing produces convection • Movement slow enough for Dispersion to matter • Mass transfer between vascular and tissue compartments 1) 2)
Excretion: The Kidney Governing Equations Kidney Information: Blood enters kidney vascular system Some flow rate transferred to bladder by GFR Rest passes through capillaries where it can interact with tissue From tissue moves to bladder, where excreted 1) 2) 3)
Case Studies: Limitations Model Requires ≈200 Physiological Parameters: 1) Drug Differences 2) Human Differences 3) Literature Limitations Our Strategy: 1) Values have data/theory behind them 2) Human differences don’t matter 3) Drug parameters optimized to reproduce data
Case Study: Atenolol Important Information Acts to treat hypertension Acts in the brain 11.1% of the dose was absorbed into the brain The compartmental model does not predict the double peak
Case Study: Imatinib Mesylate Important Information Common anti-cancer drug Acts within tissue where tumors located 18.9% was absorbed into the esophagus and stomach The compartmental model and our model predict similar results for the blood concentration, however, the compartmental model would be unable to predict tissue concentrations.
Conclusion: • What did we accomplish? • A mathematical model that can accurately describe the way a drug moves through the body • Integrated all organs at the capillary level - A novel approach • Includes spatial variations in all body tissues • What should be done in the future? • Develop a method to determine the model parameters • Account for differences between people • Compare more extensively to simpler models