Development of PBPK/PD Model to Predict Impact of Genetic Polymorphisms on the Pharmacokinetics and Pharmacodynamicsrepresented by receptor/transporter occupancy of CNS drugs

Saeed Alqahtani, Amal Kaddoumi

Department of Basic Pharmaceutical Sciences, School of Pharmacy, University of Louisiana at Monroe, Monroe, LA, 71201, USA

Corresponding author:

Dr. Amal Kaddoumi, School of Pharmacy, University of Louisiana at Monroe, 1800 Bienville Dr., Monroe, LA 71201. Phone: (318) 342-1460; Fax: (318) 342-1737;

E-mail:

General workflow of PBPK model development and verification criterion

The general workflow of PBPK modeling and simulation of test compounds in normal patients and patients with genetic variations consisted of the following steps (see Fig. 1).First,the PBPK model was initially developed using the physicochemical, biopharmaceutical, and pharmacokinetic parameters obtained from literature or estimated by ADMET Predictor(Simulation Plus Inc.); simulated mean plasma concentration-time profiles were compared with those obtained from clinical studies including single and multiple oral dosing.Population-dependent physiological parameters in human PBPK models were obtained using the Population Estimates for Age-Related Physiology™ module in GastroPlus, where the generic demographic and physiological parameters (e.g., age, tissue volume, and perfusion rate) for a normal American population were based on the National Health and Nutrition Examination Survey as described previously[1].Each simulation trail included 50 virtual subjects with demographics in terms of age, ethnicity, and percentage of women similar to those reported in the corresponding clinical trails[2-12].Human pharmacokinetics profiles of test compounds were taken from literature by scanning with GetData Graph Digitizer version 2.26 and were then used in the modeling. Second, the drug-specific parameters underwent refinement if predicted PK profileand parameters deviate significantly (<0.8-fold or >1.2-fold) from the observed data. Third, the PBPK model was used to predict the plasma concentration-time profiles and PK parameters in individuals with genetic variations by incorporating these alterations and compared to observed genetic polymorphisms data. Fourth, the PBPK model was linked to a PD model to evaluate the relationship between brain drug concentration and drug effect (Figure 1).

Verification of the established PBPK model was primarily based on AUC and Cmax. Mean AUC, Cmax, and Cminwere predicted and compared with published data. The predicted mean population PK parameters of the drug should fall within 80–120% of the observed values (i.e. the predicted/observed ratio should fall within 0.8-1.2).

Table S1. Tissue-to-Plasma Partition Coefficients (Kp) of quetiapine and fluvoxamine used in PBPK models.

Tissue / Quetiapine / Fluvoxamine
Lung / 1.15 / 60
Adipose / 5.16 / 2.03
Muscle / 5.04 / 22
Liver / 12.55 / 38
Spleen / 8.63 / 20
Heart / 6.75 / 10.48
Brain / 7.56 / 20
Kidney / 12.62 / 9.0
Skin / 4.92 / 6.19
Red Bone Marrow / 5.85 / 3.26
Yellow Bone Marrow / 5.16 / 2.03
Rest of Body / 5.64 / 3.26
Reproductive organ / 12.63 / 14.78

References

  1. Xia B, Heimbach T, Gollen R, Nanavati C, He H. A simplified PBPK modeling approach for prediction of pharmacokinetics of four primarily renally excreted and CYP3A metabolized compounds during pregnancy. AAPS J. 2013;15(4):1012-24. doi:10.1208/s12248-013-9505-3.
  2. DeVane CL, Nemeroff CB. Clinical pharmacokinetics of quetiapine: an atypical antipsychotic. Clin Pharmacokinet. 2001;40(7):509-22. doi:10.2165/00003088-200140070-00003.
  3. Grimm SW, Richtand NM, Winter HR, Stams KR, Reele SB. Effects of cytochrome P450 3A modulators ketoconazole and carbamazepine on quetiapine pharmacokinetics. Br J Clin Pharmacol. 2006;61(1):58-69. doi:10.1111/j.1365-2125.2005.02507.x.
  4. Kim KA, Joo HJ, Lee HM, Park JY. Influence of ABCB1 and CYP3A5 genetic polymorphisms on the pharmacokinetics of quetiapine in healthy volunteers. Pharmacogenet Genomics. 2014;24(1):35-42. doi:10.1097/FPC.0000000000000020.
  5. Nikisch G, Baumann P, Oneda B, Kiessling B, Weisser H, Mathe AA, et al. Cytochrome P450 and ABCB1 genetics: association with quetiapine and norquetiapine plasma and cerebrospinal fluid concentrations and with clinical response in patients suffering from schizophrenia. A pilot study. J Psychopharmacol. 2011;25(7):896-907. doi:10.1177/0269881110389208.
  6. Nord M, Nyberg S, Brogren J, Jucaite A, Halldin C, Farde L. Comparison of D(2) dopamine receptor occupancy after oral administration of quetiapine fumarate immediate-release and extended-release formulations in healthy subjects. Int J Neuropsychopharmacol. 2011;14(10):1357-66. doi:10.1017/S1461145711000514.
  7. van der Weide K, van der Weide J. The influence of the CYP3A4*22 polymorphism on serum concentration of quetiapine in psychiatric patients. J Clin Psychopharmacol. 2014;34(2):256-60. doi:10.1097/JCP.0000000000000070.
  8. Carrillo JA, Dahl ML, Svensson JO, Alm C, Rodriguez I, Bertilsson L. Disposition of fluvoxamine in humans is determined by the polymorphic CYP2D6 and also by the CYP1A2 activity. Clin Pharmacol Ther. 1996;60(2):183-90. doi:10.1016/S0009-9236(96)90134-4.
  9. Spigset O, Granberg K, Hagg S, Soderstrom E, Dahlqvist R. Non-linear fluvoxamine disposition. Br J Clin Pharmacol. 1998;45(3):257-63.
  10. Strauss WL, Layton ME, Dager SR. Brain elimination half-life of fluvoxamine measured by 19F magnetic resonance spectroscopy. Am J Psychiatry. 1998;155(3):380-4.
  11. Takano A, Suhara T, Ichimiya T, Yasuno F, Suzuki K. Time course of in vivo 5-HTT transporter occupancy by fluvoxamine. J Clin Psychopharmacol. 2006;26(2):188-91. doi:10.1097/01.jcp.0000203201.34323.d3.
  12. Watanabe J, Suzuki Y, Fukui N, Sugai T, Ono S, Inoue Y, et al. Dose-dependent effect of the CYP2D6 genotype on the steady-state fluvoxamine concentration. Ther Drug Monit. 2008;30(6):705-8. doi:10.1097/FTD.0b013e31818d73b3.

1