Application of Physiologically Based Pharmacokinetic (PBPK) Modeling to Support Dose Selection: Report of an FDA Public Workshop on PBPK

Physiologically Based Pharmacokinetic Model Qualification and Reporting Procedures for Regulatory Submissions: A Consortium Perspective

This work provides a perspective on the qualification and verification of physiologically based pharmacokinetic (PBPK) platforms/models intended for regulatory submission based on the collective experience of the Simcyp Consortium members. Examples of regulatory submission of PBPK analyses across various intended applications are presented and discussed. European Medicines Agency (EMA) and US Food and Drug Administration (FDA) recent draft guidelines regarding PBPK analyses and reporting are encouraging, and to advance the use and acceptability of PBPK analyses, more clarity and flexibility are warranted.

In vitro studies and in silico predictions of fluconazole and CYP2C9 genetic polymorphism impact on siponimod metabolism and pharmacokinetics

Purpose

The purpose of the study is to investigate the enzyme(s) responsible for siponimod metabolism and to predict the inhibitory effects of fluconazole as well as the impact of cytochrome P450 (CYP) 2C9 genetic polymorphism on siponimod pharmacokinetics (PK) and metabolism.

Methods

In vitro metabolism studies were conducted using human liver microsomes (HLM), and enzyme phenotyping was assessed using a correlation analysis method. SimCYP, a physiologically based PK model, was developed and used to predict the effects of fluconazole and CYP2C9 genetic polymorphism on siponimod metabolism. Primary PK parameters were generated using the SimCYP and WinNonlin software.

Results

Correlation analysis suggested that CYP2C9 is the main enzyme responsible for siponimod metabolism in humans. Compared with the CYP2C9*1/*1 genotype, HLM incubations from CYP2C9*3/*3 and CYP2C9*2/*2 donors showed ~ 10- and 3-fold decrease in siponimod metabolism, respectively. Simulations of enzyme contribution predicted that in the CYP2C9*1/*1 genotype, CYP2C9 is predominantly responsible for siponimod metabolism (~ 81%), whereas in the CYP2C9*3/*3 genotype, its contribution is reduced to 11%. The predicted exposure increase of siponimod with fluconazole 200 mg was 2.0–2.4-fold for CYP2C9*1/*1 genotype. In context of single dosing, the predicted mean area under the curve (AUC) is 2.7-, 3.0- and 4.5-fold higher in the CYP2C9*2/*2, CYP2C9*2/*3 and CYP2C9*3/*3 genotypes, respectively, compared with the CYP2C9*1/*1 genotype.

Conclusion

.Enzyme phenotyping with correlation analysis confirmed the predominant role of CYP2C9 in the biotransformation of siponimod and demonstrated the functional consequence of CYP2C9 genetic polymorphism on siponimod metabolism. Simulation of fluconazole inhibition closely predicted a 2-fold AUC change (ratio within ~ 20% deviation) to the observed value. In silico simulation predicted a significant reduction in siponimod clearance in the CYP2C9*2/*2 and CYP2C9*3/*3 genotypes based on the in vitro metabolism data; the predicted exposure was close (within 30%) to the observed results for the CYP2C9*2/*3 and CYP2C9*3/*3 genotypes.

Electronic supplementary material

The online version of this article (10.1007/s00228-017-2404-2) contains supplementary material, which is available to authorized users.

Ribociclib Bioavailability Is Not Affected by Gastric pH Changes or Food Intake: In Silico and Clinical Evaluations

Ribociclib (KISQALI), a cyclin‐dependent kinase 4/6 inhibitor approved for the first‐line treatment of HR+/HER2– advanced breast cancer with an aromatase inhibitor, is administered with no restrictions on concomitant gastric pH‐elevating agents or food intake. The influence of proton pump inhibitors (PPIs) on ribociclib bioavailability was assessed using 1) biorelevant media solubility, 2) physiologically based pharmacokinetic (PBPK) modeling, 3) noncompartmental analysis (NCA) of clinical trial data, and 4) population PK (PopPK) analysis. This multipronged approach indicated no effect of gastric pH changes on ribociclib PK and served as a platform for supporting ribociclib labeling language, stating no impact of gastric pH‐altering agents on the absorption of ribociclib, without a dedicated drug–drug interaction trial. The bioequivalence of ribociclib exposure with or without a high‐fat meal was demonstrated in a clinical trial. Lack of restrictions on ribociclib dosing may facilitate better patient compliance and therefore clinical benefit.

In vitro to in vivo extrapolation for high throughput prioritization and decision making

In vitro chemical safety testing methods offer the potential for efficient and economical tools to provide relevant assessments of human health risk. To realize this potential, methods are needed to relate in vitro effects to in vivo responses, i.e., in vitro to in vivo extrapolation (IVIVE). Currently available IVIVE approaches need to be refined before they can be utilized for regulatory decision-making. To explore the capabilities and limitations of IVIVE within this context, the U.S. Environmental Protection Agency Office of Research and Development and the National Toxicology Program Interagency Center for the Evaluation of Alternative Toxicological Methods co-organized a workshop and webinar series. Here, we integrate content from the webinars and workshop to discuss activities and resources that would promote inclusion of IVIVE in regulatory decision-making. We discuss properties of models that successfully generate predictions of in vivo doses from effective in vitro concentration, including the experimental systems that provide input parameters for these models, areas of success, and areas for improvement to reduce model uncertainty. Finally, we provide case studies on the uses of IVIVE in safety assessments, which highlight the respective differences, information requirements, and outcomes across various approaches when applied for decision-making.

Predictive Performance of Physiologically Based Pharmacokinetic Models for the Effect of Food on Oral Drug Absorption: Current Status

A comprehensive search in literature and published US Food and Drug Administration reviews was conducted to assess whether physiologically based pharmacokinetic (PBPK) modeling could be prospectively used to predict clinical food effect on oral drug absorption. Among the 48 resulted food effect predictions, ∼50% were predicted within 1.25‐fold of observed, and 75% within 2‐fold. Dissolution rate and precipitation time were commonly optimized parameters when PBPK modeling was not able to capture the food effect. The current work presents a knowledgebase for documenting PBPK experience to predict food effect.

Prediction of drug-drug interactions using physiologically-based pharmacokinetic models of CYP450 modulators included in Simcyp software

In recent years, physiologically based PharmacoKinetic (PBPK) modeling has received growing interest as a useful tool for the assessment of drug pharmacokinetics. It has been demonstrated to be informative and helpful to quantify the modification in drug exposure due to specific physio-pathological conditions, age, genetic polymorphisms, ethnicity and particularly drug-drug interactions (DDIs). In this paper, the prediction success of DDIs involving various cytochrome P450 isoenzyme (CYP) modulators namely ketoconazole (a competitive inhibitor of CYP3A), itraconazole (a competitive inhibitor of CYP3A), clarithromycin (a mechanism-based inhibitor of CYP3A), quinidine (a competitive inhibitor of CYP2D6), paroxetine (a mechanism-based inhibitor of CYP2D6), ciprofloxacin (a competitive inhibitor of CYP1A2), fluconazole (a competitive inhibitor of CYP2C9/2C19) and rifampicin (an inducer of CYP3A) were assessed using Simcyp^® software. The aim of this report was to establish confidence in each CYP-specific modulator file so they can be used in the future for the prediction of DDIs involving new victim compounds. Our evaluation of these PBPK models suggested that they can be successfully used to evaluate DDIs in untested scenarios. The only noticeable exception concerned a quinidine inhibitor model that requires further improvement. Additionally, other important aspects such as model validation criteria were discussed.

Predictive Performance of Physiologically Based Pharmacokinetic (PBPK) Modeling of Drugs Extensively Metabolized by Major Cytochrome P450s in Children

The accuracy of physiologically based pharmacokinetic (PBPK) model prediction in children, especially those younger than 2 years old, has not been systematically evaluated. The aim of this study was to characterize the pediatric predictive performance of the PBPK approach for 10 drugs extensively metabolized by CYP1A2 (theophylline), CYP2C8 (desloratidine, montelukast), CYP2C9 (diclofenac), CYP2C19 (esomeprazole, lansoprazole), CYP2D6 (tramadol), and CYP3A4 (itraconazole, ondansetron, sufentanil). Model performance in children was evaluated by comparing simulated plasma concentration-time profiles with observed clinical results for each drug and age group. PBPK models reasonably predicted the pharmacokinetics of desloratadine, diclofenac, itraconazole, lansoprazole, montelukast, ondansetron, sufentanil, theophylline, and tramadol across all age groups. Collectively, 58 out of 67 predictions were within 2-fold and 43 out of 67 predictions within 1.5-fold of observed values. Developed PBPK models can reasonably predict exposure in children age 1 month and older for an array of predominantly CYP metabolized drugs.

Combining "Bottom-up" and "Top-down" Approaches to Assess the Impact of Food and Gastric pH on Pictilisib (GDC-0941) Pharmacokinetics

Pictilisib, a weakly basic compound, is an orally administered, potent, and selective pan-inhibitor of phosphatidylinositol 3-kinases for oncology indications. To investigate the significance of high-fat food and gastric pH on pictilisib pharmacokinetics (PK) and enable label recommendations, a dedicated clinical study was conducted in healthy volunteers, whereby both top-down (population PK, PopPK) and bottom-up (physiologically based PK, PBPK) approaches were applied to enhance confidence of recommendation and facilitate the clinical development through scenario simulations. The PopPK model identified food (for absorption rate constant (Ka )) and proton pump inhibitors (PPI, for relative bioavailability (Frel ) and Ka ) as significant covariates. Food and PPI also impacted the variability of Frel . The PBPK model accounted for the supersaturation tendency of pictilisib, and gastric emptying physiology successfully predicted the food and PPI effect on pictilisib absorption. Our research highlights the importance of applying both quantitative approaches to address critical drug development questions.

Physiologically Based Pharmacokinetic Model of the CYP2D6 Probe Atomoxetine: Extrapolation to Special Populations and Drug-Drug Interactions

Physiologically based pharmacokinetic (PBPK) modeling of drug disposition and drug-drug interactions (DDIs) has become a key component of drug development. PBPK modeling has also been considered as an approach to predict drug disposition in special populations. However, whether models developed and validated in healthy populations can be extrapolated to special populations is not well established. The goal of this study was to determine whether a drug-specific PBPK model validated using healthy populations could be used to predict drug disposition in specific populations and in organ impairment patients. A full PBPK model of atomoxetine was developed using a training set of pharmacokinetic (PK) data from CYP2D6 genotyped individuals. The model was validated using drug-specific acceptance criteria and a test set of 14 healthy subject PK studies. Population PBPK models were then challenged by simulating the effects of ethnicity, DDIs, pediatrics, and renal and hepatic impairment on atomoxetine PK. Atomoxetine disposition was successfully predicted in 100% of healthy subject studies, 88% of studies in Asians, 79% of DDI studies, and 100% of pediatric studies. However, the atomoxetine area under the plasma concentration versus time curve (AUC) was overpredicted by 3- to 4-fold in end stage renal disease and hepatic impairment. The results show that validated PBPK models can be extrapolated to different ethnicities, DDIs, and pediatrics but not to renal and hepatic impairment patients, likely due to incomplete understanding of the physiologic changes in these conditions. These results show that systematic modeling efforts can be used to further refine population models to improve the predictive value in this area.

Prediction of a Therapeutic Dose for Buagafuran, a Potent Anxiolytic Agent by Physiologically Based Pharmacokinetic/Pharmacodynamic Modeling Starting from Pharmacokinetics in Rats and Human

Physiologically based pharmacokinetic (PBPK)/pharmacodynamic (PD) models can contribute to animal-to-human extrapolation and therapeutic dose predictions. Buagafuran is a novel anxiolytic agent and phase I clinical trials of buagafuran have been completed. In this paper, a potentially effective dose for buagafuran of 30 mg t.i.d. in human was estimated based on the human brain concentration predicted by a PBPK/PD modeling. The software GastroPlus^TM was used to build the PBPK/PD model for buagafuran in rat which related the brain tissue concentrations of buagafuran and the times of animals entering the open arms in the pharmacological model of elevated plus-maze. Buagafuran concentrations in human plasma were fitted and brain tissue concentrations were predicted by using a human PBPK model in which the predicted plasma profiles were in good agreement with observations. The results provided supportive data for the rational use of buagafuran in clinic.

Humans Vary, So Cardiac Models Should Account for That Too!

The utilization of mathematical modeling and simulation in drug development encompasses multiple mathematical techniques and the location of a drug candidate in the development pipeline. Historically speaking they have been used to analyze experimental data (i.e., Hill equation) and clarify the involved physical and chemical processes (i.e., Fick laws and drug molecule diffusion). In recent years the advanced utilization of mathematical modeling has been an important part of the regulatory review process. Physiologically based pharmacokinetic (PBPK) models identify the need to conduct specific clinical studies, suggest specific study designs and propose appropriate labeling language. Their application allows the evaluation of the influence of intrinsic (e.g., age, gender, genetics, disease) and extrinsic [e.g., dosing schedule, drug-drug interactions (DDIs)] factors, alone or in combinations, on drug exposure and therefore provides accurate population assessment. A similar pathway has been taken for the assessment of drug safety with cardiac safety being one the most advanced examples. Mechanistic mathematical model-informed safety evaluation, with a focus on drug potential for causing arrhythmias, is now discussed as an element of the Comprehensive in vitro Proarrhythmia Assay. One of the pillars of this paradigm is the use of an in silico model of the adult human ventricular cardiomyocyte to integrate in vitro measured data. Existing examples (in vitro—in vivo extrapolation with the use of PBPK models) suggest that deterministic, epidemiological and clinical data based variability models can be merged with the mechanistic models describing human physiology. There are other methods available, based on the stochastic approach and on population of models generated by randomly assigning specific parameter values (ionic current conductance and kinetic) and further pruning. Both approaches are briefly characterized in this manuscript, in parallel with the drug-specific variability.

Investigating the state of physiologically based kinetic modelling practices and challenges associated with gaining regulatory acceptance of model applications

Physiologically based kinetic (PBK) models are used widely throughout a number of working sectors, including academia and industry, to provide insight into the dosimetry related to observed adverse health effects in humans and other species. Use of these models has increased over the last several decades, especially in conjunction with emerging alternative methods to animal testing, such as in vitro studies and data-driven in silico quantitative-structure-activity-relationship (QSAR) predictions. Experimental information derived from these new approach methods can be used as input for model parameters and allows for increased confidence in models for chemicals that did not have in vivo data for model calibration. Despite significant advancements in good modelling practice (GMP) for model development and evaluation, there remains some reluctance among regulatory agencies to use such models during the risk assessment process. Here, the results of a survey disseminated to the modelling community are presented in order to inform the frequency of use and applications of PBK models in science and regulatory submission. Additionally, the survey was designed to identify a network of investigators involved in PBK modelling and knowledgeable of GMP so that they might be contacted in the future for peer review of PBK models, especially in regards to vetting the models to such a degree as to gain a greater acceptance for regulatory purposes.

•

Physiologically Based kinetic (PBK) models are used widely in academia, industry, and government.

•

Good modelling practice (GMP) for model development and evaluation continues to expand.

•

Further guidance for establishing GMP is called for.

•

There remains some reluctance among regulatory agencies to use PBK models.

•

The next generation of PBK models could be developed using only data from in vitro and in silico methods.

The Application of Physiologically Based Pharmacokinetic Modeling to Predict the Role of Drug Transporters: Scientific and Regulatory Perspectives

Transporters play an important role in drug absorption, disposition, and drug action. The evaluation of drug transporters requires a comprehensive understanding of transporter biology and pharmacology. Physiologically based pharmacokinetic (PBPK) models may offer an integrative platform to quantitatively evaluate the role of drug transporters and its interplay with other drug disposition processes such as passive drug diffusion and elimination by metabolizing enzymes. To date, PBPK modeling and simulations integrating drug transporters lag behind that for drug-metabolizing enzymes. In addition, predictive performance of PBPK has not been well established for predicting the role of drug transporters in the pharmacokinetics of a drug. To enhance overall predictive performance of transporter-based PBPK models, it is necessary to have a detailed understanding of transporter biology for proper representation in the models and to have a quantitative understanding of the contribution of transporters in the absorption and metabolism of a drug. This article summarizes PBPK-based submissions evaluating the role of drug transporters to the Office of Clinical Pharmacology of the US Food and Drug Administration.

Evaluation of transporters in drug development: Current status and contemporary issues

Transporters govern the access of molecules to cells or their exit from cells, thereby controlling the overall distribution of drugs to their intracellular site of action. Clinically relevant drug-drug interactions mediated by transporters are of increasing interest in drug development. Drug transporters, acting alone or in concert with drug metabolizing enzymes, can play an important role in modulating drug absorption, distribution, metabolism and excretion, thus affecting the pharmacokinetics and/or pharmacodynamics of a drug. The drug interaction guidance documents from regulatory agencies include various decision criteria that may be used to predict the need for in vivo assessment of transporter-mediated drug-drug interactions. Regulatory science research continues to assess the prediction performances of various criteria as well as to examine the strength and limitations of each prediction criterion to foster discussions related to harmonized decision criteria that may be used to facilitate global drug development. This review discusses the role of transporters in drug development with a focus on methodologies in assessing transporter-mediated drug-drug interactions, challenges in both in vitro and in vivo assessments of transporters, and emerging transporter research areas including biomarkers, assessment of tissue concentrations, and effect of diseases on transporters.

Physiologically Based Pharmacokinetic Modeling for Predicting the Effect of Intrinsic and Extrinsic Factors on Darunavir or Lopinavir Exposure Coadministered With Ritonavir

Management of comorbidities and medications is complex in HIV-1-infected patients. The overall objective of this project was to develop separate physiologically based pharmacokinetic (PBPK) substrate models for the protease inhibitors darunavir and lopinavir. These protease inhibitors are used in the treatment of HIV infection. Both darunavir and lopinavir are coadministered with another medication that inhibits cytochrome (CYP) 3A. The current project focused on PBPK modeling for darunavir and lopinavir coadministered with ritonavir. Darunavir and lopinavir PBPK models that accounted for ritonavir CYP3A inhibition effects (linked PBPK models) were developed. The linked PBPK models were then used to predict the effect on darunavir or lopinavir exposure from CYP modulators. In the next step, the predicted effect of hepatic impairment was evaluated. Additional exploratory analyses predicted CYP3A inhibition effects on darunavir or lopinavir exposure in simulated hepatically impaired subjects. The linked PBPK models reasonably predicted darunavir or lopinavir exposure based on simulations with CYP inhibitors or inducers. Exploratory simulations using the linked darunavir or lopinavir PBPK models indicated CYP3A inhibition may further increase darunavir or lopinavir exposure in patients with hepatic impairment.

Quantitative Modeling and Simulation in PMDA: A Japanese Regulatory Perspective

In Japan in October 2016, the Pharmaceuticals and Medical Devices Agency (PMDA) began to receive electronic data in new drug applications (NDAs). These electronic data are useful to conduct regulatory assessment of sponsors’ submissions and contribute to the PMDA's research. In this article, we summarize the number of submissions of quantitative modeling and simulation (M&S) documents in NDAs in Japan, and we describe our current thinking and activities about quantitative M&S in PMDA.

Development of a paediatric physiologically based pharmacokinetic model to assess the impact of drug-drug interactions in tuberculosis co-infected malaria subjects: A case study with artemether-lumefantrine and the CYP3A4-inducer rifampicin

The fixed dosed combination of artemether and lumefantrine (AL) is widely used for the treatment of malaria in adults and children in sub-Sahara Africa, with lumefantrine day 7 concentrations being widely used as a marker for clinical efficacy. Both are substrates for CYP3A4 and susceptible to drug-drug interactions (DDIs); indeed, knowledge of the impact of these factors is currently sparse in paediatric population groups. Confounding malaria treatment is the co-infection of patients with tuberculosis. The concomitant treatment of AL with tuberculosis chemotherapy, which includes the CYP3A4 inducer rifampicin, increases the risk of parasite recrudescence and malaria treatment failure. This study developed a population-based PBPK model for AL in adults capable of predicting the pharmacokinetics of AL under non-DDI and DDI conditions, as well as predicting AL pharmacokinetics in paediatrics of 2-12years of age. The validated model was utilised to assess the concomitant treatment of rifampicin and lumefantrine under standard body-weight based treatment regimens for 2-5year olds, and demonstrated that no subjects attained the target day 7 concentration (C_d7) of 280ng/mL, highlighting the importance of this DDI and the potential risk of malaria-TB based DDIs. An adapted 7-day treatment regimen was simulated and resulted in 63% and 74.5% of subjects attaining the target C_d7 for 1-tablet and 2-tablet regimens respectively.

Development and Qualification of Physiologically Based Pharmacokinetic Models for Drugs With Atypical Distribution Behavior: A Desipramine Case Study

Desipramine is a secondary tricyclic amine, which is primarily metabolized by cytochrome 2D6. It shows a high volume of distribution (Vss) (10–50 L/kg) due to its high lipophilicity, unspecific phospholipid binding, and lysosomal trapping. The objective of this study was to develop and qualify a physiologically based pharmacokinetic (PBPK) model for desipramine, which accounts for the high Vss of the drug following intravenous and oral administration of doses up to 100 mg. The model also accounts for the extended time to reach maximum concentration after oral dosing due to enterocyte trapping. Once developed and qualified in adults, we characterized the dynamic changes in metabolism and pharmacokinetics of desipramine after birth by scaling the system‐specific parameters of the model from adults to pediatrics. The developed modeling strategy provides a prototypical workflow that can also be applied to other drugs with similar properties and a high volume of distribution.

Physiologically Based Pharmacokinetic Modeling of Drug Transporters to Facilitate Individualized Dose Prediction

Physiologically based pharmacokinetic modeling is a commonly used strategy in the drug development and regulatory submissions. This commentary provides a critical overview of the current status of physiologically based pharmacokinetic methodologies to predict transporter-mediated pharmacokinetics, in addition to the impact of disease and genetics with respect to local and systemic concentration.

Impact of physiologically based pharmacokinetic models on regulatory reviews and product labels: Frequent utilization in the field of oncology

Physiologically based pharmacokinetic (PBPK) modeling can be used to predict drug pharmacokinetics in virtual populations using models that integrate understanding of physiological systems. PBPK models have been widely utilized for predicting pharmacokinetics in clinically untested scenarios during drug applications and regulatory reviews in recent years. Here, we provide a comprehensive review of the application of PBPK in new drug application (NDA) review documents from the US Food and Drug Administration (FDA) in the past 4 years.

Delineating the Role of Various Factors in Renal Disposition of Digoxin through Application of Physiologically Based Kidney Model to Renal Impairment Populations

Development of submodels of organs within physiologically-based pharmacokinetic (PBPK) principles and beyond simple perfusion limitations may be challenging because of underdeveloped in vitro-in vivo extrapolation approaches or lack of suitable clinical data for model refinement. However, advantage of such models in predicting clinical observations in divergent patient groups is now commonly acknowledged. Mechanistic understanding of altered renal secretion in renal impairment is one area that may benefit from such models, despite knowledge gaps in renal pathophysiology. In the current study, a PBPK kidney model was developed for digoxin, accounting for the roles of organic anion transporting peptide 4C1 (OATP4C1) and P-glycoprotein (P-gp) in its tubular secretion, with the aim to investigate the impact of age and renal impairment (moderate to severe) on renal drug disposition. Initial PBPK simulations based on changes in glomerular filtration rate (GFR) underestimated the observed reduction in digoxin renal excretion clearance (CLR) in subjects with moderately impaired renal function relative to healthy. Reduction in either proximal tubule cell number or the OATP4C1 abundance in the mechanistic kidney model successfully predicted 59% decrease in digoxin CLR, in particular when these changes were proportional to reduction in GFR. In contrast, predicted proximal tubule concentration of digoxin was only sensitive to changes in the transporter expression/ million proximal tubule cells. Based on the mechanistic modeling, reduced proximal tubule cellularity and OATP4C1 abundance, and inhibition of OATP4C1-mediated transport, are proposed as possible causes of reduced digoxin renal secretion in renally impaired patients.

Report from the EMA workshop on qualification and reporting of physiologically based pharmacokinetic (PBPK) modeling and simulation

On Nov 21, 2016, the European Medicines Agency (EMA) hosted a workshop to discuss its draft guideline on qualification and reporting of physiologically based pharmacokinetic (PBPK) analysis.¹ Published on July 21, 2016, the draft PBPK guideline is currently under the period of public comments.

Evaluation of Drug-Drug Interaction Potential Between Sacubitril/Valsartan (LCZ696) and Statins Using a Physiologically Based Pharmacokinetic Model

Sacubitril/valsartan (LCZ696) has been approved for the treatment of heart failure. Sacubitril is an in vitro inhibitor of organic anion-transporting polypeptides (OATPs). In clinical studies, LCZ696 increased atorvastatin C_max by 1.7-fold and area under the plasma concentration-time curve by 1.3-fold, but had little or no effect on simvastatin or simvastatin acid exposure. A physiologically based pharmacokinetics modeling approach was applied to explore the underlying mechanisms behind the statin-specific LCZ696 drug interaction observations. The model incorporated OATP-mediated clearance (CL_int,T) for simvastatin and simvastatin acid to successfully describe the pharmacokinetic profiles of either analyte in the absence or presence of LCZ696. Moreover, the model successfully described the clinically observed drug effect with atorvastatin. The simulations clarified the critical parameters responsible for the observation of a low, yet clinically relevant, drug-drug interaction DDI between sacubitril and atorvastatin and the lack of effect with simvastatin acid. Atorvastatin is administered in its active form and rapidly achieves C_max that coincide with the low C_max of sacubitril. In contrast, simvastatin requires a hydrolysis step to the acid form and therefore is not present at the site of interactions at sacubitril concentrations that are inhibitory. Similar models were used to evaluate the drug-drug interaction risk for additional OATP-transported statins which predicted to maximally result in a 1.5-fold exposure increase.

A four-compartment PBPK heart model accounting for cardiac metabolism - model development and application

In the field of cardiac drug efficacy and safety assessment, information on drug concentration in heart tissue is desirable. Because measuring drug concentrations in human cardiac tissue is challenging in healthy volunteers, mathematical models are used to cope with such limitations. With a goal of predicting drug concentration in cardiac tissue, we have developed a whole-body PBPK model consisting of seventeen perfusion-limited compartments. The proposed PBPK heart model consisted of four compartments: the epicardium, midmyocardium, endocardium, and pericardial fluid, and accounted for cardiac metabolism using CYP450. The model was written in R. The plasma:tissues partition coefficients (Kp) were calculated in Simcyp Simulator. The model was fitted to the concentrations of amitriptyline in plasma and the heart. The estimated parameters were as follows: 0.80 for the absorption rate [h-1], 52.6 for Kprest, 0.01 for the blood flow through the pericardial fluid [L/h], and 0.78 for the P-parameter describing the diffusion between the pericardial fluid and epicardium [L/h]. The total cardiac clearance of amitriptyline was calculated as 0.316 L/h. Although the model needs further improvement, the results support its feasibility, and it is a first attempt to provide an active drug concentration in various locations within heart tissue using a PBPK approach.

Physiologically based pharmacokinetic modeling in regulatory decision-making at the European Medicines Agency

Physiologically based pharmacokinetic (PBPK) modeling is a valuable tool in drug development and regulatory assessment, as it offers the opportunity to simulate the pharmacokinetics of a compound, with a mechanistic understanding, in a variety of populations and situations. This work reviews the use and impact of such modeling in selected regulatory procedures submitted to the European Medicines Agency (EMA) before the end of 2015, together with its subsequent reflection in public documents relating to the assessment of these procedures. It is apparent that the reference to PBPK modeling in regulatory public documents underrepresents its use. A positive trend over time of the number of PBPK models submitted is shown, and in a number of cases the results of these may impact the decision-making process or lead to recommendations in the product labeling. These results confirm the need for regulatory guidance in this field, which is currently under development by the EMA.

Applied Concepts in PBPK Modeling: How to Build a PBPK/PD Model

The aim of this tutorial is to introduce the fundamental concepts of physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) modeling with a special focus on their practical implementation in a typical PBPK model building workflow. To illustrate basic steps in PBPK model building, a PBPK model for ciprofloxacin will be constructed and coupled to a pharmacodynamic model to simulate the antibacterial activity of ciprofloxacin treatment.

Physiologically-based pharmacokinetic models: approaches for enabling personalized medicine

Personalized medicine strives to deliver the 'right drug at the right dose' by considering inter-person variability, one of the causes for therapeutic failure in specialized populations of patients. Physiologically-based pharmacokinetic (PBPK) modeling is a key tool in the advancement of personalized medicine to evaluate complex clinical scenarios, making use of physiological information as well as physicochemical data to simulate various physiological states to predict the distribution of pharmacokinetic responses. The increased dependency on PBPK models to address regulatory questions is aligned with the ability of PBPK models to minimize ethical and technical difficulties associated with pharmacokinetic and toxicology experiments for special patient populations. Subpopulation modeling can be achieved through an iterative and integrative approach using an adopt, adapt, develop, assess, amend, and deliver methodology. PBPK modeling has two valuable applications in personalized medicine: (1) determining the importance of certain subpopulations within a distribution of pharmacokinetic responses for a given drug formulation and (2) establishing the formulation design space needed to attain a targeted drug plasma concentration profile. This review article focuses on model development for physiological differences associated with sex (male vs. female), age (pediatric vs. young adults vs. elderly), disease state (healthy vs. unhealthy), and temporal variation (influence of biological rhythms), connecting them to drug product formulation development within the quality by design framework. Although PBPK modeling has come a long way, there is still a lengthy road before it can be fully accepted by pharmacologists, clinicians, and the broader industry.

Virtual Clinical Trial Toward Polytherapy Safety Assessment: Combination of Physiologically Based Pharmacokinetic/Pharmacodynamic-Based Modeling and Simulation Approach With Drug-Drug Interactions Involving Terfenadine as an Example

A Quantitative Systems Pharmacology approach was utilized to predict the cardiac consequences of drug-drug interaction (DDI) at the population level. The Simcyp in vitro-in vivo correlation and physiologically based pharmacokinetic platform was used to predict the pharmacokinetic profile of terfenadine following co-administration of the drug. Electrophysiological effects were simulated using the Cardiac Safety Simulator. The modulation of ion channel activity was dependent on the inhibitory potential of drugs on the main cardiac ion channels and a simulated free heart tissue concentration. ten Tusscher's human ventricular cardiomyocyte model was used to simulate the pseudo-ECG traces and further predict the pharmacodynamic consequences of DDI. Consistent with clinical observations, predicted plasma concentration profiles of terfenadine show considerable intra-subject variability with recorded C_max values below 5 ng/mL for most virtual subjects. The pharmacokinetic and pharmacodynamic effects of inhibitors were predicted with reasonable accuracy. In all cases, a combination of the physiologically based pharmacokinetic and physiology-based pharmacodynamic models was able to differentiate between the terfenadine alone and terfenadine + inhibitor scenario. The range of QT prolongation was comparable in the clinical and virtual studies. The results indicate that mechanistic in vitro-in vivo correlation can be applied to predict the clinical effects of DDI even without comprehensive knowledge on all mechanisms contributing to the interaction.

PBPK modeling and simulation in drug research and development

Physiologically based pharmacokinetic (PBPK) modeling and simulation can be used to predict the pharmacokinetic behavior of drugs in humans using preclinical data. It can also explore the effects of various physiologic parameters such as age, ethnicity, or disease status on human pharmacokinetics, as well as guide dose and dose regiment selection and aid drug-drug interaction risk assessment. PBPK modeling has developed rapidly in the last decade within both the field of academia and the pharmaceutical industry, and has become an integral tool in drug discovery and development. In this mini-review, the concept and methodology of PBPK modeling are briefly introduced. Several case studies were discussed on how PBPK modeling and simulation can be utilized through various stages of drug discovery and development. These case studies are from our own work and the literature for better understanding of the absorption, distribution, metabolism and excretion (ADME) of a drug candidate, and the applications to increase efficiency, reduce the need for animal studies, and perhaps to replace clinical trials. The regulatory acceptance and industrial practices around PBPK modeling and simulation is also discussed.

Translational PK/PD modeling to increase probability of success in drug discovery and early development

In this review we present ways in which translational PK/PD modeling can address opportunities to enhance probability of success in drug discovery and early development. This is achieved by impacting efficacy and safety-driven attrition rates, through increased focus on the quantitative understanding and modeling of translational PK/PD. Application of the proposed principles early in the discovery and development phases is anticipated to bolster confidence of successfully evaluating proof of mechanism in humans and ultimately improve Phase II success. The present review is centered on the application of predictive modeling and simulation approaches during drug discovery and early development, and more specifically of mechanism-based PK/PD modeling. Case studies are presented, focused on the relevance of M&S contributions to real-world questions and the impact on decision making.

Perspectives in regulatory science: translational and clinical pharmacology

This paper focuses on the role of clinical and translational pharmacology in the drug development and the regulatory process. Contemporary regulatory issues faced by FDA's Office of Clinical Pharmacology (OCP) in fulfilling its mission to advance the science of drug response and translate patient diversity into optimal drug therapy are discussed. Specifically current focus of the following key aspects of the drug development and regulatory science processes are discussed: the OCP vision and mission, two key OCP initiatives (i.e. guidance modernization, labeling and health communications), and translational and clinical pharmacology related regulatory science issues in (i.e. uncertainty, breakthrough therapies, individualization).

Predictive Performance of Physiologically Based Pharmacokinetic and Population Pharmacokinetic Modeling of Renally Cleared Drugs in Children

Predictive performance of physiologically based pharmacokinetic (PBPK) and population pharmacokinetic (PopPK) models of drugs predominantly eliminated through kidney in the pediatric population was evaluated. After optimization using adult clinical data, the verified PBPK models can predict 33 of 34 drug clearance within twofold of the observed values in children 1 month and older. More specifically, 10 of 11 of predicted clearance values were within 1.5‐fold of those observed in children between 1 month and 2 years old. The PopPK approach also predicted 19 of 21 drug clearance within twofold of the observed values in children. In summary, our analysis demonstrated both PBPK and PopPK adult models, after verification with additional adult pharmacokinetic (PK) studies and incorporation of known ontogeny of renal filtration, could be applied for dosing regimen recommendation in children 1 month and older for renally eliminated drugs in a first‐in‐pediatric study.

Key to Opening Kidney for In Vitro-In Vivo Extrapolation Entrance in Health and Disease: Part II: Mechanistic Models and In Vitro-In Vivo Extrapolation

It is envisaged that application of mechanistic models will improve prediction of changes in renal disposition due to drug-drug interactions, genetic polymorphism in enzymes and transporters and/or renal impairment. However, developing and validating mechanistic kidney models is challenging due to the number of processes that may occur (filtration, secretion, reabsorption and metabolism) in this complex organ. Prediction of human renal drug disposition from preclinical species may be hampered by species differences in the expression and activity of drug metabolising enzymes and transporters. A proposed solution is bottom-up prediction of pharmacokinetic parameters based on in vitro-in vivo extrapolation (IVIVE), mediated by recent advances in in vitro experimental techniques and application of relevant scaling factors. This review is a follow-up to the Part I of the report from the 2015 AAPS Annual Meeting and Exhibition (Orlando, FL; 25th-29th October 2015) which focuses on IVIVE and mechanistic prediction of renal drug disposition. It describes the various mechanistic kidney models that may be used to investigate renal drug disposition. Particular attention is given to efforts that have attempted to incorporate elements of IVIVE. In addition, the use of mechanistic models in prediction of renal drug-drug interactions and potential for application in determining suitable adjustment of dose in kidney disease are discussed. The need for suitable clinical pharmacokinetics data for the purposes of delineating mechanistic aspects of kidney models in various scenarios is highlighted.

A Physiologically Based Pharmacokinetic Model to Describe Artemether Pharmacokinetics in Adult and Pediatric Patients

Artemether is co-administered with lumefantrine as part of a fixed-dose combination therapy for malaria in both adult and pediatric patients. However, artemether exposure is higher in younger infants (1-3 months) with a lower body weight (<5 kg) as compared to older infants (3-6 months) with a higher body weight (≥5 to <10 kg), children, and adults. In contrast, lumefantrine exposure is similar in all age groups. This article describes the clinically observed artemether exposure data in pediatric populations across various age groups (1 month to 12 years) and body weights (<5 or ≥5 kg) using physiologically based pharmacokinetic (PBPK) mechanistic models. A PBPK model was developed using artemether physicochemical, biopharmaceutic, and metabolic properties together with known enzyme ontogeny and pediatric physiology. The model was verified using clinical data from adult patients after multiple doses of oral artemether, and was then applied to simulate the exposure in children and infants. The simulated PBPK concentration-time profiles captured observed clinical data. Consistent with the clinical data, the PBPK model simulations indicated a higher artemether exposure for younger infants with lower body weight. A PBPK model developed for artemether reliably described the clinical data from adult and pediatric patients.

Physiologically Based Absorption Modeling to Explore the Impact of Food and Gastric pH Changes on the Pharmacokinetics of Alectinib

Alectinib, a lipophilic, basic, anaplastic lymphoma kinase (ALK) inhibitor with very low aqueous solubility, has received Food and Drug Administration-accelerated approval for the treatment of patients with ALK+ non-small-cell lung cancer. This paper describes the application of physiologically based absorption modeling during clinical development to predict and understand the impact of food and gastric pH changes on alectinib absorption. The GastroPlus^™ software was used to develop an absorption model integrating in vitro and in silico data on drug substance properties. Oral pharmacokinetics was simulated by linking the absorption model to a disposition model fit to pharmacokinetic data obtained after an intravenous infusion. Simulations were compared to clinical data from a food effect study and a drug-drug interaction study with esomeprazole, a gastric acid-reducing agent. Prospective predictions of a positive food effect and negligible impact of gastric pH elevation were confirmed with clinical data, although the exact magnitude of the food effect could not be predicted with confidence. After optimization of the absorption model with clinical food effect data, a refined model was further applied to derive recommendations on the timing of dose administration with respect to a meal. The application of biopharmaceutical absorption modeling is an area with great potential to further streamline late stage drug development and with impact on regulatory questions.

Role of Quantitative Clinical Pharmacology in Pediatric Approval and Labeling

Dose selection is one of the key decisions made during drug development in pediatrics. There are regulatory initiatives that promote the use of model-based drug development in pediatrics. Pharmacometrics or quantitative clinical pharmacology enables development of models that can describe factors affecting pharmacokinetics and/or pharmacodynamics in pediatric patients. This manuscript describes some examples in which pharmacometric analysis was used to support approval and labeling in pediatrics. In particular, the role of pharmacokinetic (PK) comparison of pediatric PK to adults and utilization of dose/exposure-response analysis for dose selection are highlighted. Dose selection for esomeprazole in pediatrics was based on PK matching to adults, whereas for adalimumab, exposure-response, PK, efficacy, and safety data together were useful to recommend doses for pediatric Crohn's disease. For vigabatrin, demonstration of similar dose-response between pediatrics and adults allowed for selection of a pediatric dose. Based on model-based pharmacokinetic simulations and safety data from darunavir pediatric clinical studies with a twice-daily regimen, different once-daily dosing regimens for treatment-naïve human immunodeficiency virus 1-infected pediatric subjects 3 to <12 years of age were evaluated. The role of physiologically based pharmacokinetic modeling (PBPK) in predicting pediatric PK is rapidly evolving. However, regulatory review experiences and an understanding of the state of science indicate that there is a lack of established predictive performance of PBPK in pediatric PK prediction. Moving forward, pharmacometrics will continue to play a key role in pediatric drug development contributing toward decisions pertaining to dose selection, trial designs, and assessing disease similarity to adults to support extrapolation of efficacy.

Physiologically Based Pharmacokinetic Modeling in Pediatric Oncology Drug Development

Childhood cancer represents more than 100 rare and ultra-rare diseases, with an estimated 12,400 new cases diagnosed each year in the United States. As such, this much smaller patient population has led to pediatric oncology drug development lagging behind that for adult cancers. Developing drugs for pediatric malignancies also brings with it a number of unique trial design considerations, including flexible enrollment approaches, age-appropriate formulation, acceptable sampling schedules, and balancing the need for age-stratified dosing regimens, given the smaller patient populations. The regulatory landscape for pediatric pharmacotherapy has evolved with U.S. Food and Drug Administration (FDA) legislation such as the 2012 FDA Safety and Innovation Act. In parallel, regulatory authorities have recommended the application of physiologically based pharmacokinetic (PBPK) modeling, for example, in the recently issued FDA Strategic Plan for Accelerating the Development of Therapies for Pediatric Rare Diseases. PBPK modeling provides a quantitative and systems-based framework that allows the effects of intrinsic and extrinsic factors on drug exposure to be modeled in a mechanistic fashion. The application of PBPK modeling in drug development for pediatric cancers is relatively nascent, with several retrospective analyses of cytotoxic therapies, and latterly for targeted agents such as obatoclax and imatinib. More recently, we have employed PBPK modeling in a prospective manner to inform the first pediatric trials of pinometostat and tazemetostat in genetically defined populations (mixed lineage leukemia-rearranged and integrase interactor-1-deficient sarcomas, respectively). In this review, we evaluate the application of PBPK modeling in pediatric cancer drug development and discuss the important challenges that lie ahead in this field.

In Silico Modeling of Gastrointestinal Drug Absorption: Predictive Performance of Three Physiologically Based Absorption Models

Gastrointestinal (GI) drug absorption is a complex process determined by formulation, physicochemical and biopharmaceutical factors, and GI physiology. Physiologically based in silico absorption models have emerged as a widely used and promising supplement to traditional in vitro assays and preclinical in vivo studies. However, there remains a lack of comparative studies between different models. The aim of this study was to explore the strengths and limitations of the in silico absorption models Simcyp 13.1, GastroPlus 8.0, and GI-Sim 4.1, with respect to their performance in predicting human intestinal drug absorption. This was achieved by adopting an a priori modeling approach and using well-defined input data for 12 drugs associated with incomplete GI absorption and related challenges in predicting the extent of absorption. This approach better mimics the real situation during formulation development where predictive in silico models would be beneficial. Plasma concentration-time profiles for 44 oral drug administrations were calculated by convolution of model-predicted absorption-time profiles and reported pharmacokinetic parameters. Model performance was evaluated by comparing the predicted plasma concentration-time profiles, Cmax, tmax, and exposure (AUC) with observations from clinical studies. The overall prediction accuracies for AUC, given as the absolute average fold error (AAFE) values, were 2.2, 1.6, and 1.3 for Simcyp, GastroPlus, and GI-Sim, respectively. The corresponding AAFE values for Cmax were 2.2, 1.6, and 1.3, respectively, and those for tmax were 1.7, 1.5, and 1.4, respectively. Simcyp was associated with underprediction of AUC and Cmax; the accuracy decreased with decreasing predicted fabs. A tendency for underprediction was also observed for GastroPlus, but there was no correlation with predicted fabs. There were no obvious trends for over- or underprediction for GI-Sim. The models performed similarly in capturing dependencies on dose and particle size. In conclusion, it was shown that all three software packages are useful to guide formulation development. However, as a consequence of the high fraction of inaccurate predictions (prediction error >2-fold) and the clear trend toward decreased accuracy with decreased predicted fabs observed with Simcyp, the results indicate that GI-Sim and GastroPlus perform better than Simcyp in predicting the intestinal absorption of the incompletely absorbed drugs when a higher degree of accuracy is needed. In addition, this study suggests that modeling and simulation research groups should perform systematic model evaluations using their own input data to maximize confidence in model performance and output.

Atomoxetine: A Review of Its Pharmacokinetics and Pharmacogenomics Relative to Drug Disposition

Atomoxetine is a selective norepinephrine (NE) reuptake inhibitor approved for the treatment of attention-deficit/hyperactivity disorder (ADHD) in children (≥6 years of age), adolescents, and adults. Its metabolism and disposition are fairly complex, and primarily governed by cytochrome P450 (CYP) 2D6 (CYP2D6), whose protein expression varies substantially from person to person, and by race and ethnicity because of genetic polymorphism. These differences can be substantial, resulting in 8-10-fold differences in atomoxetine exposure between CYP2D6 poor metabolizers and extensive metabolizers. In this review, we have attempted to revisit and analyze all published clinical pharmacokinetic data on atomoxetine inclusive of public access documents from the new drug application submitted to the United States Food and Drug Administration (FDA). The present review focuses on atomoxetine metabolism, disposition, and genetic polymorphisms of CYP2D6 as they specifically relate to atomoxetine, and provides an in-depth discussion of the fundamental pharmacokinetics of the drug including its absorption, distribution, metabolism, and excretion in pediatric and adult populations. Further, a summary of relationships between genetic variants of CYP2D6 and to some degree, CYP2C19, are provided with respect to atomoxetine plasma concentrations, central nervous system (CNS) pharmacokinetics, and associated clinical implications for pharmacotherapy. Lastly, dosage adjustments based on pharmacokinetic principles are discussed.

Recent Advances in Development and Application of Physiologically-Based Pharmacokinetic (PBPK) Models: a Transition from Academic Curiosity to Regulatory Acceptance

There is a renewed surge of interest in applications of physiologically-based pharmacokinetic (PBPK) models by the pharmaceutical industry and regulatory agencies. Developing PBPK models within a systems pharmacology context allows separation of the parameters pertaining to the animal or human body (the system) from that of the drug and the study design which is essential to develop generic drug-independent models used to extrapolate PK/PD properties in various healthy and patient populations. This has expanded the classical paradigm to a ‘predict-learn-confirm-apply’ concept. Recently, a number of drug labels are informed by simulation results generated using PBPK models. These cases show that either the simulations are used in lieu of conducting clinical studies or have informed the drug label that otherwise would have been silent in some specific situations. It will not be surprising to see applications of these models in implementing precision dosing at the point of care in the near future.