Application of PBPK modelling in drug discovery and development at Pfizer

Considerations from an International Regulatory and Pharmaceutical Industry (IQ MPS Affiliate) Workshop on the Standardization of Complex In Vitro Models in Drug Development

In May 2022, there is an International Regulatory and Pharmaceutical Industry (Innovation and Quality [IQ] Microphysiological Systems [MPS] Affiliate) Workshop on the standardization of complex in vitro models (CIVMs) in drug development. This manuscript summarizes the discussions and conclusions of this joint workshop organized and executed by the IQ MPS Affiliate and the United States Food and Drug Administration (FDA). A key objective of the workshop is to facilitate discussions around opportunities and/or needs for standardization of MPS and chart potential pathways to increase model utilization in the context of regulatory decision making. Participation in the workshop included 200 attendees from the FDA, IQ MPS Affiliate, and 26 global regulatory organizations and affiliated parties representing Europe, Japan, and Canada. It is agreed that understanding global perspectives regarding the readiness of CIVM/MPS models for regulatory decision making and potential pathways to gaining acceptance is useful to align on globally. The obstacles are currently too great to develop standards for every context of use (COU). Instead, it is suggested that a more tractable approach may be to think of broadly applicable standards that can be applied regardless of COU and/or organ system. Considerations and next steps for this effort are described.

Prediction of pharmacokinetics of an anaplastic lymphoma kinase inhibitor in rat and monkey: application of physiologically based pharmacokinetic model as an alternative tool to minimise animal studies

The pharmacokinetic (PK) and toxicokinetic profile of a drug from its preclinical evaluation helps the researcher determine whether the drug should be tested in humans based on its safety and toxicity.Preclinical studies require time and resources and are prone to error. Moreover, according to the United States Food and Drug Administration Modernisation Act 2, animal testing is no longer mandatory for new drug development, and an animal-free alternative, such as cell-based assay and computer models, can be used.Different physiologically based PK models were developed for an anaplastic lymphoma kinase inhibitor in rats and monkeys after intravenous and oral administration using its physicochemical properties and in vitro characterisation data.The developed model was validated against the in vivo data available in the literature, and the validation results were found within the acceptable limit. A parameter sensitivity analysis was performed to identify the properties of the compound influencing the PK profile.This work demonstrates the application of the physiologically based PK model to predict the PKs of a drug, which will eventually assist in reducing the number of animal studies and save time and cost of drug discovery and development.

Progressive tools and critical strategies for development of best fit PBPK model aiming better in vitro-in vivo correlation

Nowadays, conducting discriminative dissolution experiments employing physiologically based pharmacokinetic modeling (PBPK) or physiologically based biopharmaceutical modeling (PBBM) is gaining significant importance in quantitatively predicting oral absorption of drugs. Mechanistic understanding of each process involved in drug absorption and its impact on the performance greatly facilitates designing a formulation with high confidence. Unfortunately, the biggest challenge scientists are facing in current days is the lack of standardized protocol for integrating dissolution experiment data during PBPK modeling. However, in vitro-in vivo drug release interrelation can be improved with the consideration and development of appropriate biorelevant dissolution media that closely mimic physiological conditions. Multiple reported dissolution models have described nature and functionality of different regions of the gastrointestinal tract (GI) to more accurately design discriminative dissolution media. Dissolution experiment data can be integrated either mechanistically or without a mechanism depending primarily on the formulation type, biopharmaceutics classification system (BCS) class and particle size of the drug substance. All such parameters are required to be considered for selecting the appropriate functions during PBPK modeling to produce a best fit model. The primary focus of this review is to critically discuss various progressive dissolution models and tools, existing challenges and approaches for establishing best fit PBPK model aiming better in vitro-in vivo correlation (IVIVC). Strategies for proper selection of dissolution models as an input function in PBPK/PBBM modeling have also been critically discussed. Logical and scientific pathway for selection of different type of functions and integration events in the commercially available in silico software has been described through case studies.

In Silico Studies to Support Vaccine Development

The progress that has been made in computer science positioned in silico studies as an important and well-recognized methodology in the drug discovery and development process. It has numerous advantages in terms of costs and also plays a huge impact on the way the research is conducted since it can limit the use of animal models leading to more sustainable research. Currently, human trials are already being partly replaced by in silico trials. EMA and FDA are both endorsing these studies and have been providing webinars and guidance to support them. For instance, PBPK modeling studies are being used to gather data on drug interactions with other drugs and are also being used to support clinical and regulatory requirements for the pediatric population, pregnant women, and personalized medicine. This trend evokes the need to understand the role of in silico studies in vaccines, considering the importance that these products achieved during the pandemic and their promising hope in oncology. Vaccines are safer than other current oncology treatments. There is a huge variety of strategies for developing a cancer vaccine, and some of the points that should be considered when designing the vaccine technology are the following: delivery platforms (peptides, lipid-based carriers, polymers, dendritic cells, viral vectors, etc.), adjuvants (to boost and promote inflammation at the delivery site, facilitating immune cell recruitment and activation), choice of the targeted antigen, the timing of vaccination, the manipulation of the tumor environment, and the combination with other treatments that might cause additive or even synergistic anti-tumor effects. These and many other points should be put together to outline the best vaccine design. The aim of this article is to perform a review and comprehensive analysis of the role of in silico studies to support the development of and design of vaccines in the field of oncology and infectious diseases. The authors intend to perform a literature review of all the studies that have been conducted so far in preparing in silico models and methods to support the development of vaccines. From this point, it was possible to conclude that there are few in silico studies on vaccines. Despite this, an overview of how the existing work could support the design of vaccines is described.

Current State and Challenges of Physiologically Based Biopharmaceutics Modeling (PBBM) in Oral Drug Product Development

Physiologically based biopharmaceutics modeling (PBBM) emphasizes the integration of physicochemical properties of drug substance and formulation characteristics with system physiological parameters to predict the absorption and pharmacokinetics (PK) of a drug product. PBBM has been successfully utilized in drug development from discovery to postapproval stages and covers a variety of applications. The use of PBBM facilitates drug development and can reduce the number of preclinical and clinical studies. In this review, we summarized the major applications of PBBM, which are classified into six categories: formulation selection and development, biopredictive dissolution method development, biopharmaceutics risk assessment, clinically relevant specification settings, food effect evaluation and pH-dependent drug-drug-interaction risk assessment. The current state of PBBM applications is illustrated with examples from published studies for each category of application. Despite the variety of PBBM applications, there are still many hurdles limiting the use of PBBM in drug development, that are associated with the complexity of gastrointestinal and human physiology, the knowledge gap between the in vitro and the in vivo behavior of drug products, the limitations of model interfaces, and the lack of agreed model validation criteria, among other issues. The challenges and essential considerations related to the use of PBBM are discussed in a question-based format along with the scientific thinking on future research directions. We hope this review can foster open discussions between the pharmaceutical industry and regulatory agencies and encourage collaborative research to fill the gaps, with the ultimate goal to maximize the applications of PBBM in oral drug product development.

In Vitro/In Vivo Translation of Synergistic Combination of MDM2 and MEK Inhibitors in Melanoma Using PBPK/PD Modelling: Part III

The development of in vitro/in vivo translational methods and a clinical trial framework for synergistically acting drug combinations are needed to identify optimal therapeutic conditions with the most effective therapeutic strategies. We performed physiologically based pharmacokinetic-pharmacodynamic (PBPK/PD) modelling and virtual clinical trial simulations for siremadlin, trametinib, and their combination in a virtual representation of melanoma patients. In this study, we built PBPK/PD models based on data from in vitro absorption, distribution, metabolism, and excretion (ADME), and in vivo animals' pharmacokinetic-pharmacodynamic (PK/PD) and clinical data determined from the literature or estimated by the Simcyp simulator (version V21). The developed PBPK/PD models account for interactions between siremadlin and trametinib at the PK and PD levels. Interaction at the PK level was predicted at the absorption level based on findings from animal studies, whereas PD interaction was based on the in vitro cytotoxicity results. This approach, combined with virtual clinical trials, allowed for the estimation of PK/PD profiles, as well as melanoma patient characteristics in which this therapy may be noninferior to the dabrafenib and trametinib drug combination. PBPK/PD modelling, combined with virtual clinical trial simulation, can be a powerful tool that allows for proper estimation of the clinical effect of the above-mentioned anticancer drug combination based on the results of in vitro studies. This approach based on in vitro/in vivo extrapolation may help in the design of potential clinical trials using siremadlin and trametinib and provide a rationale for their use in patients with melanoma.

In-Depth Analysis of Physiologically Based Pharmacokinetic (PBPK) Modeling Utilization in Different Application Fields Using Text Mining Tools

In the past decade, only a small number of papers have elaborated on the application of physiologically based pharmacokinetic (PBPK) modeling across different areas. In this review, an in-depth analysis of the distribution of PBPK modeling in relation to its application in various research topics and model validation was conducted by text mining tools. Orange 3.32.0, an open-source data mining program was used for text mining. PubMed was used for data retrieval, and the collected articles were analyzed by several widgets. A total of 2699 articles related to PBPK modeling met the predefined criteria. The number of publications per year has been rising steadily. Regarding the application areas, the results revealed that 26% of the publications described the use of PBPK modeling in early drug development, risk assessment and toxicity assessment, followed by absorption/formulation modeling (25%), prediction of drug-disease interactions (20%), drug-drug interactions (DDIs) (17%) and pediatric drug development (12%). Furthermore, the analysis showed that only 12% of the publications mentioned model validation, of which 51% referred to literature-based validation and 26% to experimentally validated models. The obtained results present a valuable review of the state-of-the-art regarding PBPK modeling applications in drug discovery and development and related fields.

Towards best use and regulatory acceptance of generic physiologically based kinetic (PBK) models for in vitro-to-in vivo extrapolation (IVIVE) in chemical risk assessment

With an increasing need to incorporate new approach methodologies (NAMs) in chemical risk assessment and the concomitant need to phase out animal testing, the interpretation of in vitro assay readouts for quantitative hazard characterisation becomes more important. Physiologically based kinetic (PBK) models, which simulate the fate of chemicals in tissues of the body, play an essential role in extrapolating in vitro effect concentrations to in vivo bioequivalent exposures. As PBK-based testing approaches evolve, it will become essential to standardise PBK modelling approaches towards a consensus approach that can be used in quantitative in vitro-to-in vivo extrapolation (QIVIVE) studies for regulatory chemical risk assessment based on in vitro assays. Based on results of an ECETOC expert workshop, steps are recommended that can improve regulatory adoption: (1) define context and implementation, taking into consideration model complexity for building fit-for-purpose PBK models, (2) harmonise physiological input parameters and their distribution and define criteria for quality chemical-specific parameters, especially in the absence of in vivo data, (3) apply Good Modelling Practices (GMP) to achieve transparency and design a stepwise approach for PBK model development for risk assessors, (4) evaluate model predictions using alternatives to in vivo PK data including read-across approaches, (5) use case studies to facilitate discussions between modellers and regulators of chemical risk assessment. Proof-of-concepts of generic PBK modelling approaches are published in the scientific literature at an increasing rate. Working on the previously proposed steps is, therefore, needed to gain confidence in PBK modelling approaches for regulatory use.

Evaluation of the Success of High-Throughput Physiologically Based Pharmacokinetic (HT-PBPK) Modeling Predictions to Inform Early Drug Discovery

Minimizing in vitro and in vivo testing in early drug discovery with the use of physiologically based pharmacokinetic (PBPK) modeling and machine learning (ML) approaches has the potential to reduce discovery cycle times and animal experimentation. However, the prediction success of such an approach has not been shown for a larger and diverse set of compounds representative of a lead optimization pipeline. In this study, the prediction success of the oral (PO) and intravenous (IV) pharmacokinetics (PK) parameters in rats was assessed using a “bottom-up” approach, combining in vitro and ML inputs with a PBPK model. More than 240 compounds for which all of the necessary inputs and PK data were available were used for this assessment. Different clearance scaling approaches were assessed, using hepatocyte intrinsic clearance and protein binding as inputs. In addition, a novel high-throughput PBPK (HT-PBPK) approach was evaluated to assess the scalability of PBPK predictions for a larger number of compounds in drug discovery. The results showed that bottom-up PBPK modeling was able to predict the rat IV and PO PK parameters for the majority of compounds within a 2- to 3-fold error range, using both direct scaling and dilution methods for clearance predictions. The use of only ML-predicted inputs from the structure did not perform well when using in vitro inputs, likely due to clearance miss predictions. The HT-PBPK approach produced comparable results to the full PBPK modeling approach but reduced the simulation time from hours to seconds. In conclusion, a bottom-up PBPK and HT-PBPK approach can successfully predict the PK parameters and guide early discovery by informing compound prioritization, provided that good in vitro assays are in place for key parameters such as clearance.

Artificial Intelligence in Drug Safety and Metabolism

The use of artificial intelligence methods in drug safety began in the early 2000s with applications such as predicting bacterial mutagenicity and hERG inhibition. The field has been endlessly expanding ever since and the models have become more complex. These approaches are now integrated into molecule risk assessment processes along with in vitro and in vivo methods. Today, artificial intelligence can be used in every phase of drug discovery and development, from profiling chemical libraries in early discovery, to predicting off-target effects in the mid-discovery phase, to assessing potential mutagenic impurities in development and degradants as part of life cycle management. This chapter provides an overview of artificial intelligence in drug safety and describes its application throughout the entire discovery and development process.

A Mechanistic Site-Of-Action Model: A Tool for Informing Right Target, Right Compound, And Right Dose for Therapeutic Antagonistic Antibody Programs

Quantitative modeling is increasingly utilized in the drug discovery and development process, from the initial stages of target selection, through clinical studies. The modeling can provide guidance on three major questions–is this the right target, what are the right compound properties, and what is the right dose for moving the best possible candidate forward. In this manuscript, we present a site-of-action modeling framework which we apply to monoclonal antibodies against soluble targets. We give a comprehensive overview of how we construct the model and how we parametrize it and include several examples of how to apply this framework for answering the questions postulated above. The utilities and limitations of this approach are discussed.

Physiologically based pharmacokinetic (PBPK) modeling for prediction of celecoxib pharmacokinetics according to CYP2C9 genetic polymorphism

Celecoxib is a non-steroidal anti-inflammatory drug (NSAID) and a representative selective cyclooxygenase (COX)-2 inhibitor, which is commonly prescribed for osteoarthritis, rheumatoid arthritis, ankylosing spondylitis, acute pain, and primary dysmenorrhea. It is mainly metabolized by CYP2C9 and partly by CYP3A4 after oral administration. Many studies reported that CYP2C9 genetic polymorphism has significant effects on the pharmacokinetics of celecoxib and the occurrence of adverse drug reactions. The aim of this study was to develop a physiologically based pharmacokinetic (PBPK) model of celecoxib according to CYP2C9 genetic polymorphism for personalized pharmacotherapy. Initially, a clinical pharmacokinetic study was conducted where a single dose (200 mg) of celecoxib was administered to 39 healthy Korean subjects with CYP2C9*1/*1 or CYP2C9*1/*3 genotypes to obtain data for PBPK development. Based on the conducted pharmacokinetic study and a previous pharmacokinetic study involving subjects with CYP2C9*1/*13 and CYP2C9*3/*3 genotype, PBPK model for celecoxib was developed. A PBPK model for CYP2C9*1/*1 genotype group was developed and then scaled to other genotype groups (CYP2C9*1/*3, CYP2C9*1/*13 and CYP2C9*3/*3). After model development, model validation was performed with comparison of five pharmacokinetic studies. As a result, the developed PBPK model of celecoxib successfully described the pharmacokinetics of each CYP2C9 genotype group and its predicted values were within the acceptance criterion. Additionally, all the predicted values were within two-fold error range in comparison to the previous pharmacokinetic studies. This study demonstrates the possibility of determining the appropriate dosage of celecoxib for each individual through the PBPK modeling with CYP2C9 genomic information. This approach could contribute to the reduction of adverse drug reactions of celecoxib and enable precision medicine.

Intestinal Excretion, Intestinal Recirculation, and Renal Tubule Reabsorption Are Underappreciated Mechanisms That Drive the Distribution and Pharmacokinetic Behavior of Small Molecule Drugs

Drug reabsorption following biliary excretion is well-known as enterohepatic recirculation (EHR). Renal tubular reabsorption (RTR) following renal excretion is also common but not easily assessed. Intestinal excretion (IE) and enteroenteric recirculation (EER) have not been recognized as common disposition mechanisms for metabolically stable and permeable drugs. IE and intestinal reabsorption (IR:EHR/EER), as well as RTR, are governed by dug concentration gradients, passive diffusion, active transport, and metabolism, and together they markedly impact disposition and pharmacokinetics (PK) of small molecule drugs. Disruption of IE, IR, or RTR through applications of active charcoal (AC), transporter knockout (KO), and transporter inhibitors can lead to changes in PK parameters. The impacts of intestinal and renal reabsorption on PK are under-appreciated. Although IE and EER/RTR can be an intrinsic drug property, there is no apparent strategy to optimize compounds based on this property. This review seeks to improve understanding and applications of IE, IR, and RTR mechanisms.

Physiologically-Based Pharmacokinetic (PBPK) Modeling of Buprenorphine in Adults, Children and Preterm Neonates

Buprenorphine plays a crucial role in the therapeutic management of pain in adults, adolescents and pediatric subpopulations. However, only few pharmacokinetic studies of buprenorphine in children, particularly neonates, are available as conducting clinical trials in this population is especially challenging. Physiologically-based pharmacokinetic (PBPK) modeling allows the prediction of drug exposure in pediatrics based on age-related physiological differences. The aim of this study was to predict the pharmacokinetics of buprenorphine in pediatrics with PBPK modeling. Moreover, the drug-drug interaction (DDI) potential of buprenorphine with CYP3A4 and P-glycoprotein perpetrator drugs should be elucidated. A PBPK model of buprenorphine and norbuprenorphine in adults has been developed and scaled to children and preterm neonates, accounting for age-related changes. One-hundred-percent of the predicted AUC_last values in adults (geometric mean fold error (GMFE): 1.22), 90% of individual AUC_last predictions in children (GMFE: 1.54) and 75% in preterm neonates (GMFE: 1.57) met the 2-fold acceptance criterion. Moreover, the adult model was used to simulate DDI scenarios with clarithromycin, itraconazole and rifampicin. We demonstrate the applicability of scaling adult PBPK models to pediatrics for the prediction of individual plasma profiles. The novel PBPK models could be helpful to further investigate buprenorphine pharmacokinetics in various populations, particularly pediatric subgroups.

Prediction of active human dose: learnings from 20 years of Merck KGaA experience, illustrated by case studies

High-quality dose predictions based on a good understanding of target engagement is one of the main translational goals in drug development. Here, we systematically evaluate active human dose predictions for 15 Merck KGaA/EMD Serono assets spanning several modalities and therapeutic areas. Using case studies, we illustrate the value of adhering to the translational best practices of having an exposure-response relationship in an appropriate animal model; having validated, translatable pharmacodynamic (PD) biomarkers measurable in Phase I populations in the right tissue; having a deeper understanding of biology; and capturing uncertainties in predictions. Given the gap in publications on the subject, we believe that the learnings from this unique diverse data set, which are generic to the industry, will trigger actions to improve future predictions.

In Vitro Dissolution and in Silico Modeling Shortcuts in Bioequivalence Testing

Purpose: To review in vitro testing and simulation platforms that are in current use to predict in vivo performances of generic products as well as other situations to provide evidence for biowaiver and support drug formulations development. Methods: Pubmed and Google Scholar databases were used to review published literature over the past 10 years. The terms used were “simulation AND bioequivalence” and “modeling AND bioequivalence” in the title field of databases, followed by screening, and then reviewing. Results: A total of 22 research papers were reviewed. Computer simulation using software such as GastroPlus™, PK-Sim^® and SimCyp^® find applications in drug modeling. Considering the wide use of optimization for in silico predictions to fit observed data, a careful review of publications is required to validate the reliability of these platforms. For immediate release (IR) drug products belonging to the Biopharmaceutics Classification System (BCS) classes I and III, difference factor (ƒ₁) and similarity factor (ƒ₂) are calculated from the in vitro dissolution data of drug formulations to support biowaiver; however, this method can be more discriminatory and may not be useful for all dissolution profiles. Conclusions: Computer simulation platforms need to improve their mechanistic physiologically based pharmacokinetic (PBPK) modeling, and if prospectively validated within a small percentage of error from the observed clinical data, they can be valuable tools in bioequivalence (BE) testing and formulation development.

An in vitro toolbox to accelerate anti-malarial drug discovery and development

Background

Modelling and simulation are being increasingly utilized to support the discovery and development of new anti-malarial drugs. These approaches require reliable in vitro data for physicochemical properties, permeability, binding, intrinsic clearance and cytochrome P450 inhibition. This work was conducted to generate an in vitro data toolbox using standardized methods for a set of 45 anti-malarial drugs and to assess changes in physicochemical properties in relation to changing target product and candidate profiles.

Methods

Ionization constants were determined by potentiometric titration and partition coefficients were measured using a shake-flask method. Solubility was assessed in biorelevant media and permeability coefficients and efflux ratios were determined using Caco-2 cell monolayers. Binding to plasma and media proteins was measured using either ultracentrifugation or rapid equilibrium dialysis. Metabolic stability and cytochrome P450 inhibition were assessed using human liver microsomes. Sample analysis was conducted by LC–MS/MS.

Results

Both solubility and fraction unbound decreased, and permeability and unbound intrinsic clearance increased, with increasing Log D_7.4. In general, development compounds were somewhat more lipophilic than legacy drugs. For many compounds, permeability and protein binding were challenging to assess and both required the use of experimental conditions that minimized the impact of non-specific binding. Intrinsic clearance in human liver microsomes was varied across the data set and several compounds exhibited no measurable substrate loss under the conditions used. Inhibition of cytochrome P450 enzymes was minimal for most compounds.

Conclusions

This is the first data set to describe in vitro properties for 45 legacy and development anti-malarial drugs. The studies identified several practical methodological issues common to many of the more lipophilic compounds and highlighted areas which require more work to customize experimental conditions for compounds being designed to meet the new target product profiles. The dataset will be a valuable tool for malaria researchers aiming to develop PBPK models for the prediction of human PK properties and/or drug–drug interactions. Furthermore, generation of this comprehensive data set within a single laboratory allows direct comparison of properties across a large dataset and evaluation of changing property trends that have occurred over time with changing target product and candidate profiles.

Sampling Site Has a Critical Impact on Physiologically Based Pharmacokinetic Modeling

It has been shown that arterial (central) and venous (peripheral) plasma drug concentrations can be very different. While pharmacokinetic studies typically measure drug concentrations from the peripheral vein such as the arm vein, physiologically based pharmacokinetic (PBPK) models generally output simulated concentrations from the central venous compartment that physiologically represents the right atrium, a merge of the superior and inferior vena cava. In this study, a physiologically based peripheral forearm sampling site model was developed and verified using nicotine, ketamine, lidocaine, and fentanyl as model drugs. This verified model allows output of simulated peripheral venous concentrations that can be meaningfully compared with observed pharmacokinetic data from the arm vein. The generalized effect of PBPK model sampling site on simulation output was investigated. Drugs and metabolites with large volumes of distribution showed considerable concentration discrepancy between the simulated central venous compartment and the peripheral arm vein after intravenous or oral administration, resulting in significant differences in values for C _max and time taken to reach C _max (t _max ) In addition, the simulated central venous metabolite profile showed an unexpected profile that was not observed in the peripheral arm vein. Using fentanyl as a model compound, we show that using the wrong sampling site in PBPK models can lead to biased model evaluation and subsequent erroneous model parameter optimization. Such an error in model parameters along with the discrepant sampling site could dramatically mislead the pharmacokinetic prediction in unstudied clinical scenarios, affecting the assessment of drug safety and efficacy. Overall, this study shows that PBPK model publications should specify the model sampling sites and match them with those employed in clinical studies. SIGNIFICANCE STATEMENT: Our study shows that sampling from the central venous compartment (right atrium) during physiologically based pharmacokinetic model development gives rise to biased model evaluation and erroneous model parameterization when observed data are collected from the peripheral arm vein. This can lead to a clinically significant error in predictions of plasma concentration-time profiles in unstudied scenarios. To address this error, we developed and verified a novel peripheral sampling site model to simulate arm vein drug concentrations that can be applied to different drug dosing scenarios.

Quantitative Analysis of Tozadenant Using Liquid Chromatography-Mass Spectrometric Method in Rat Plasma and Its Human Pharmacokinetics Prediction Using Physiologically Based Pharmacokinetic Modeling

Tozadenant is one of the selective adenosine A2a receptor antagonists with a potential to be a new Parkinson's disease (PD) therapeutic drug. In this study, a liquid chromatography-mass spectrometry based bioanalytical method was qualified and applied for the quantitative analysis of tozadenant in rat plasma. A good calibration curve was observed in the range from 1.01 to 2200 ng/mL for tozadenant using a quadratic regression. In vitro and preclinical in vivo pharmacokinetic (PK) properties of tozadenant were studied through the developed bioanalytical methods, and human PK profiles were predicted using physiologically based pharmacokinetic (PBPK) modeling based on these values. The PBPK model was initially optimized using in vitro and in vivo PK data obtained by intravenous administration at a dose of 1 mg/kg in rats. Other in vivo PK data in rats were used to validate the PBPK model. The human PK of tozadenant after oral administration at a dose of 240 mg was simulated by using an optimized and validated PBPK model. The predicted human PK parameters and profiles were similar to the observed clinical data. As a result, optimized PBPK model could reasonably predict the PK in human.

Variance based global sensitivity analysis of physiologically based pharmacokinetic absorption models for BCS I-IV drugs

Regulatory agencies have a strong interest in sensitivity analysis for the evaluation of physiologically-based pharmacokinetic (PBPK) models used in pharmaceutical research and drug development and regulatory submissions. One of the applications of PBPK is the prediction of fraction absorbed and bioavailability for drugs following oral administration. In this context, we performed a variance based global sensitivity analysis (GSA) on in-house PBPK models for drug absorption, with the aim of identifying key parameters that influence the predictions of the fraction absorbed and the bioavailability for neutral, acidic and basic compounds. This analysis was done for four different classes of drugs, defined according to the Biopharmaceutics Classification System, differentiating compounds by permeability and solubility. For class I compounds (highly permeable, highly soluble), the parameters that mainly influence the fraction absorbed are related to the formulation properties, for class II compounds (highly permeable, lowly soluble) to the dissolution process, for class III (lowly permeable, highly soluble) to both absorption process and formulation properties and for class IV (lowly permeable, lowly soluble) to both absorption and dissolution processes. Considering the bioavailability, the results are similar to those for the fraction absorbed, with the addition that parameters related to gut wall and liver clearance influence as well the predictions. This work aimed to give a demonstration of the GSA methodology and highlight its importance in improving our understanding of PBPK absorption models and in guiding the choice of parameters that can safely be assumed, estimated or require data generation to allow informed model prediction.

Application of Physiologically Based Pharmacokinetic Modeling to Predict Pharmacokinetics in Healthy Japanese Subjects

Pharmacokinetics (PKs) in Japanese healthy subjects were simulated for nine compounds using physiologically based PK (PBPK) models parameterized with physicochemical properties, preclinical absorption, distribution, metabolism, and excretion (ADME) data, and clinical PK data from non‐Japanese subjects. For each dosing regimen, 100 virtual trials were simulated and predicted/observed ratios for peak plasma concentration (C_max) and area under the curve (AUC) were calculated. As qualification criteria, it was prespecified that >80% of simulated trials should demonstrate ratios to observed data ranging from 0.5–2.0. Across all compounds and dose regimens studied, 93% of simulated C_max values in Japanese subjects fulfilled the criteria. Similarly, for AUC, 77% of single‐dosing regimens and 100% of multiple‐dosing regimens fulfilled the criteria. In summary, mechanistically incorporating the appropriate ADME properties into PBPK models, followed by qualification using non‐Japanese clinical data, can predict PKs in the Japanese population and lead to efficient trial design and conduct of Japanese phase I studies.

Whole-body physiology-based pharmacokinetics of caspofungin for general patients, intensive care unit patients and hepatic insufficiency patients

Caspofungin is an echinocandin antifungal agent licensed as a first-line therapy for invasive candidiasis in patients with moderate to severe illness or recent exposure to azoles. In this study we developed a whole-body physiology-based pharmacokinetics (WB-PBPK) model to predict the pharmacokinetics (PK) of caspofungin, and combined with Monte Carlo simulation (MCS) to optimize clinical dosage regimens of caspofungin in different kinds of patients. A WB-PBPK model of caspofungin was built and validated with raw data from 4 previous trials of general patients, intensive care unit (ICU) patients with Child-Pugh B, ICU patients on continuous renal replacement therapy, mild and moderate hepatic insuffciency (HI) patients. MCS was used to optimize clinical dosage regimens of caspofungin in these patients. A cumulative fraction of response (CFR) value of ≥90% was considered to be the minimum for achieving optimal empirical therapy. The simulated results of the WB-PBPK model were in good agreement with observed values of all trials. For general and ICU patients with caspofungin 70/50 mg, AUC and Cmax were decreased with the increase of body weight (BW) and showed great variation. MCS showed all general patients achieved CFR≥90% regardless of BW. But not all ICU patients with higher BW (≥70 kg) could achieve CFR≥90%. Compared with standard dosage regimens in general patients, caspofungin 70/35 mg in ICU patients with Child-Pugh B achieved significantly decreased AUC and Cmax, but obtained similar AUC and Cmax in moderate HI patients with Child-Pugh B. The WB-PBPK model of caspofungin is able to predict PK of all populations correctly. The combined WB-PBPK model with MCS can successfully optimize clinical dosage regimens of caspofungin in all patient populations.

Physiologically based pharmacokinetic modelling of atomoxetine with regard to CYP2D6 genotypes

Atomoxetine is a norepinephrine reuptake inhibitor indicated in the treatment of attention-deficit/hyperactivity disorder. It is primarily metabolized by CYP2D6 to its equipotent metabolite, 4-hydroxyatomoxetine, which promptly undergoes further glucuronidation to an inactive 4-HAT-O-glucuronide. Clinical trials have shown that decreased CYP2D6 activity leads to substantially elevated atomoxetine exposure and increase in adverse reactions. The aim of this study was to to develop a pharmacologically based pharmacokinetic (PBPK) model of atomoxetine in different CYP2D6 genotypes. A single 20 mg dose of atomoxetine was given to 19 healthy Korean individuals with CYP2D6*wt/*wt (*wt = *1 or *2) or CYP2D6*10/*10 genotype. Based on the results of this pharmacokinetic study, a PBPK model for CYP2D6*wt/*wt individuals was developed. This model was scaled to those with CYP2D6*10/*10 genotype, as well as CYP2D6 poor metabolisers. We validated this model by comparing the predicted pharmacokinetic parameters with diverse results from the literature. The presented PBPK model describes the pharmacokinetics after single and repeated oral atomoxetine doses with regard to CYP2D6 genotype and phenotype. This model could be utilized for identification of appropriate dosages of atomoxetine in patients with reduced CYP2D6 activity to minimize the adverse events, and to enable personalised medicine.

Absolute Oral Bioavailability of Creatine Monohydrate in Rats: Debunking a Myth

Creatine is an ergogenic compound used by athletes to enhance performance. Supplementation with creatine monohydrate (CM) has been suggested for musculoskeletal and neurological disorders. Until now, little is known about its pharmacokinetic profile. Our objective was to determine the oral bioavailability of CM and the influence of dose on oral absorption. Rats were dosed orally with low dose (10 mg/kg) or high dose (70 mg/kg) 13C-labeled CM. Blood samples were removed at various time points. Muscle and brain tissue were collected at the conclusion of the study. Plasma and tissue levels of 13C-labeled creatine were determined using liquid chromatography-tandem mass spectrometry (LC-MS/MS). Physiologically based pharmacokinetic (PBPK) models of CM were built using GastroPlus™. These models were used to predict the plasma concentration-time profiles of creatine hydrochloride (CHCL), which has improved aqueous solubility compared to CM. Absolute oral bioavailability for low dose CM was 53% while high dose CM was only 16%. The simulated Cmax of 70 mg/kg CHCL was around 35 μg/mL compared to 14 μg/mL for CM with a predicted oral bioavailability of 66% with CHCL compared to 17% with CM. Our results suggest that the oral bioavailability of CM is less than complete and subject to dose and that further examination of improved dosage formulations of creatine is warranted.

Prediction of Drug-Drug Interactions with Bupropion and Its Metabolites as CYP2D6 Inhibitors Using a Physiologically-Based Pharmacokinetic Model

The potential of inhibitory metabolites of perpetrator drugs to contribute to drug-drug interactions (DDIs) is uncommon and underestimated. However, the occurrence of unexpected DDI suggests the potential contribution of metabolites to the observed DDI. The aim of this study was to develop a physiologically-based pharmacokinetic (PBPK) model for bupropion and its three primary metabolites—hydroxybupropion, threohydrobupropion and erythrohydrobupropion—based on a mixed “bottom-up” and “top-down” approach and to contribute to the understanding of the involvement and impact of inhibitory metabolites for DDIs observed in the clinic. PK profiles from clinical researches of different dosages were used to verify the bupropion model. Reasonable PK profiles of bupropion and its metabolites were captured in the PBPK model. Confidence in the DDI prediction involving bupropion and co-administered CYP2D6 substrates could be maximized. The predicted maximum concentration (C_max) area under the concentration-time curve (AUC) values and C_max and AUC ratios were consistent with clinically observed data. The addition of the inhibitory metabolites into the PBPK model resulted in a more accurate prediction of DDIs (AUC and C_max ratio) than that which only considered parent drug (bupropion) P450 inhibition. The simulation suggests that bupropion and its metabolites contribute to the DDI between bupropion and CYP2D6 substrates. The inhibitory potency from strong to weak is hydroxybupropion, threohydrobupropion, erythrohydrobupropion, and bupropion, respectively. The present bupropion PBPK model can be useful for predicting inhibition from bupropion in other clinical studies. This study highlights the need for caution and dosage adjustment when combining bupropion with medications metabolized by CYP2D6. It also demonstrates the feasibility of applying the PBPK approach to predict the DDI potential of drugs undergoing complex metabolism, especially in the DDI involving inhibitory metabolites.

Whole body physiologically based modelling of β-blockers in the rat: events in tissues and plasma following an i.v. bolus dose

BACKGROUND AND PURPOSE

Whole body physiologically based pharmacokinetic (PBPK) models have been increasingly applied in drug development to describe kinetic events of therapeutic agents in animals and humans. The advantage of such modelling is the ability to incorporate vast amounts of physiological information, such as organ blood flow and volume, to ensure that the model is as close to reality as possible.

EXPERIMENTAL APPROACH

Previous PBPK model development of enantiomers of a series of seven racemic β-blockers, namely, acebutolol, betaxolol, bisoprolol, metoprolol, oxprenolol, pindolol and propranolol, together with S-timolol in rat was based on tissue and blood concentration data at steady state. Compounds were administered in several cassettes with the composition mix and blood and tissue sampling times determined using a D-optimal design.

KEY RESULTS

Closed-loop PBPK models were developed initially based on the application of open loop forcing function models to individual tissues and compounds. For the majority of compounds and tissues, distribution kinetics was adequately characterized by perfusion rate-limited models. For some compounds in the testes and gut, a permeability rate-limited distribution model was required to best fit the data. Parameter estimates of the tissue-to-blood partition coefficient through fitting of individual enantiomers and of racemic pair were generally in agreement and also concur with those from previous steady-state experiments.

CONCLUSIONS AND IMPLICATIONS

PBPK modelling is a very powerful tool to aid drug discovery and development of therapeutic agents in animals and humans. However, careful consideration of the assumptions made during the modelling exercise is essential.

The Constraints, Construction, and Verification of a Strain-Specific Physiologically Based Pharmacokinetic Rat Model

The use of in vitro-in vivo extrapolation (IVIVE) techniques, mechanistically incorporated within physiologically based pharmacokinetic (PBPK) models, can harness in vitro drug data and enhance understanding of in vivo pharmacokinetics. This study's objective was to develop a user-friendly rat (250 g, male Sprague-Dawley) IVIVE-linked PBPK model. A 13-compartment PBPK model including mechanistic absorption models was developed, with required system data (anatomical, physiological, and relevant IVIVE scaling factors) collated from literature and analyzed. Overall, 178 system parameter values for the model are provided. This study also highlights gaps in available system data required for strain-specific rat PBPK model development. The model's functionality and performance were assessed using previous literature-sourced in vitro properties for diazepam, metoprolol, and midazolam. The results of simulations were compared against observed pharmacokinetic rat data. Predicted and observed concentration profiles in 10 tissues for diazepam after a single intravenous (i.v.) dose making use of either observed i.v. clearance (CL_iv) or in vitro hepatocyte intrinsic clearance (CL_int) for simulations generally led to good predictions in various tissue compartments. Overall, all i.v. plasma concentration profiles were successfully predicted. However, there were challenges in predicting oral plasma concentration profiles for metoprolol and midazolam, and the potential reasons and according solutions are discussed.

On the Nature of Physiologically-Based Pharmacokinetic Models -A Priori or A Posteriori? Mechanistic or Empirical?

Physiologically-based pharmacokinetic (PBPK) models explicitly incorporate tissue-specific blood flows, partition coefficients, and metabolic processes. Since PBPK models are derived using physiologic parameters and interactions of the compound with tissue components, these models are considered to be "bottom up" as opposed to "top down". Modeling approaches can be characterized as either a posteriori (observational) or a priori (based solely on theory). Furthermore, approaches can be mechanistic (structure and components based on mechanisms) or empirical (based on observations alone). Both "bottom up" and "top down" approaches can incorporate either empirical or mechanistic components. In this perspective, we discuss some of the methods and assumptions of current PBPK modeling approaches. Specifically, we discuss drug partitioning into phospholipids and neutral lipids, use of blood-plasma ratios to estimate basic drug tissue partitioning, and clearance of neutral and acidic drugs. Based on these discussions, we believe that current PBPK models are mechanistic but a posteriori and semi-empirical.

Prediction of pharmacokinetics and drug-drug interaction potential using physiologically based pharmacokinetic (PBPK) modeling approach: A case study of caffeine and ciprofloxacin

Over the last decade, physiologically based pharmacokinetics (PBPK) application has been extended significantly not only to predicting preclinical/human PK but also to evaluating the drug-drug interaction (DDI) liability at the drug discovery or development stage. Herein, we describe a case study to illustrate the use of PBPK approach in predicting human PK as well as DDI using in silico, in vivo and in vitro derived parameters. This case was composed of five steps such as: simulation, verification, understanding of parameter sensitivity, optimization of the parameter and final evaluation. Caffeine and ciprofloxacin were used as tool compounds to demonstrate the “fit for purpose” application of PBPK modeling and simulation for this study. Compared to caffeine, the PBPK modeling for ciprofloxacin was challenging due to several factors including solubility, permeability, clearance and tissue distribution etc. Therefore, intensive parameter sensitivity analysis (PSA) was conducted to optimize the PBPK model for ciprofloxacin. Overall, the increase in C_max of caffeine by ciprofloxacin was not significant. However, the increase in AUC was observed and was proportional to the administered dose of ciprofloxacin. The predicted DDI and PK results were comparable to observed clinical data published in the literatures. This approach would be helpful in identifying potential key factors that could lead to significant impact on PBPK modeling and simulation for challenging compounds.

Lipid-associated oral delivery: Mechanisms and analysis of oral absorption enhancement

The majority of newly discovered oral drugs are poorly water soluble, and co-administration with lipids has proven effective in significantly enhancing bioavailability of some compounds with low aqueous solubility. Yet, lipid-based delivery technologies have not been widely employed in commercial oral products. Lipids can impact drug transport and fate in the gastrointestinal (GI) tract through multiple mechanisms including enhancement of solubility and dissolution kinetics, enhancement of permeation through the intestinal mucosa, and triggering drug precipitation upon lipid emulsion depletion (e.g., by digestion). The effect of lipids on drug absorption is currently not quantitatively predictable, in part due to the multiple complex dynamic processes that can be impacted by lipids. Quantitative mechanistic analysis of the processes significant to lipid system function and overall impact on drug absorption can aid in the understanding of drug-lipid interactions in the GI tract and exploitation of such interactions to achieve optimal lipid-based drug delivery. In this review, we discuss the impact of co-delivered lipids and lipid digestion on drug dissolution, partitioning, and absorption in the context of the experimental tools and associated kinetic expressions used to study and model these processes. The potential benefit of a systems-based consideration of the concurrent multiple dynamic processes occurring upon co-dosing lipids and drugs to predict the impact of lipids on drug absorption and enable rational design of lipid-based delivery systems is presented.

An S-warfarin and AZD1981 interaction: in vitro and clinical pilot data suggest the N-deacetylated amino acid metabolite as the primary perpetrator

AIM

AZD1981 is an orally bioavailable chemoattractant receptor-homologous molecule expressed on Th2 cells (CRTh2) receptor antagonist progressed to phase II trials for the treatment of allergic asthma. Previously performed in vitro human hepatocyte incubations identified N-deacetylated AZD1981 as a primary metabolite. We report on metabolite exposure from a clinical excretion balance, on in vitro studies performed to determine the likelihood of a metabolite-dependent drug-drug interaction (DDI) and on a clinical warfarin DDI study. The aim was to demonstrate that N-deacetylated AZD1981 is responsible for the observed interaction.

METHODS

The excretion and biotransformation of [¹⁴ C]-AZD1981 were studied in healthy male volunteers, and subsequently in vitro cytochrome P450 (CYP) inhibition and hepatocyte uptake investigations were carried out with metabolites and the parent drug. A clinical DDI study using coadministered twice-daily 100 mg and 400 mg AZD1981 with 25 mg warfarin was performed.

RESULTS

The excretion balance study showed N-deacetylated AZD1981 to be the most abundant metabolite present in plasma. In vitro data revealed the metabolite to be a weak CYP2C9 time-dependent inhibitor, subject to more active hepatic uptake than the parent molecule. Clinically, the S-warfarin area under the plasma concentration-time curve increased, on average, 1.4-fold [95% confidence interval (CI) 1.22, 1.50] and 2.4-fold (95% CI 2.11, 2.64) after 100 mg (n = 13) and 400 mg (n = 11) AZD1981 administration, respectively. In vitro CYP inhibition and hepatocyte uptake data were used to explain the interaction.

CONCLUSIONS

N-deacetylated AZD1981 can be added to the small list of drug metabolites reported as sole contributors to clinical drug-drug interactions, with weak time-dependent inhibition exacerbated by efficient hepatic uptake being the cause.

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.

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.

Ontogeny of plasma proteins, albumin and binding of diazepam, cyclosporine, and deltamethrin

BACKGROUND

To characterize the ontogeny of plasma albumin and total proteins, due to the lack of a comprehensive pediatric database. Secondly, to establish the magnitude and duration of maturational changes in binding of highly-bound drugs/chemicals.

METHODS

Anonymized plasma samples from 296 donors were pooled in 6 age brackets from birth to adolescence. Total protein and albumin levels were measured in each age group, as was the age-dependency of plasma binding of diazepam (DZP), cyclosporine (CYC), and deltamethrin (DLM), a pyrethroid insecticide.

RESULTS

Plasma levels of albumin and total proteins steadily increased for the first 1-3 y of life. Unbound DZP and CYC fractions were elevated three- to fourfold in neonates, but decreased to adult levels after 1 and 3 y, respectively. Unbound DLM levels exceeded those in adults for just 1 mo.

CONCLUSION

Neonates and infants under 1-3 y may be at risk from increased amounts of free drug, when given standard doses of some highly-bound drugs. Pyrethroid insecticides might be anticipated to pose increased risk for 1 mo.

Preclinical pharmacokinetics of TPN729MA, a novel PDE5 inhibitor, and prediction of its human pharmacokinetics using a PBPK model

AIM

TPN729MA is a novel selective PDE5 inhibitor currently under clinical development in China for the treatment of erectile dysfunction. In this study we characterized its preclinical pharmacokinetics (PK) and predict its human PK using a physiologically based pharmacokinetic (PBPK) model.

METHODS

The preclinical PK of TPN729MA was studied in rats and dogs. Human clearance (CL) values for TPN729MA were predicted from various allometric methods and from intrinsic CL determined in human liver microsomes. Human PK and plasma concentration versus time profiles of TPN729MA were predicted by using a PBPK model in GastroPlus. Considering the uncertainties in the prediction, a preliminary human study was conducted in 3 healthy male volunteers with an oral dose of 25 mg.

RESULTS

After a single intravenous administration of TPN729MA at a dose of 1 mg/kg in rats and 3 mg/kg in dogs, the plasma CL was 69.7 mL·min(-1)·kg(-1) in rats and 26.3 mL·min(-1)·kg(-1) in dogs, and the steady-state volumes of distribution (V(ss)) were 7.35 L/kg in rats and 6.48 L/kg in dogs. The oral bioavailability of TPN729MA was 10% in rats and above 34% in dogs. Profiles of predicted plasma concentration versus time were similar to those observed in humans at 25 mg, and the predicted T(max), C(max) and AUC values were within 2-fold of the observed values.

CONCLUSION

TPN729MA demonstrates good preclinical PK. This compound is a valuable candidate for further clinical development. This study shows the benefits of using a PBPK model to predict PK in humans.

Are Physiologically Based Pharmacokinetic Models Reporting the Right C(max)? Central Venous Versus Peripheral Sampling Site

Physiologically based pharmacokinetic (PBPK) models can over-predict maximum plasma concentrations (C_max) following intravenous administration. A proposed explanation is that invariably PBPK models report the concentration in the central venous compartment, rather than the site where the samples are drawn. The purpose of this study was to identify and validate potential corrective models based on anatomy and physiology governing the blood supply at the site of sampling and incorporate them into a PBPK platform. Four models were developed and scrutinised for their corrective potential. All assumed the peripheral sampling site concentration could be described by contributions from surrounding tissues and utilised tissue-specific concentration-time profiles reported from the full-PBPK model within the Simcyp Simulator. Predicted concentrations for the peripheral site were compared to the observed C_max. The models results were compared to clinical data for 15 studies over seven compounds (alprazolam, imipramine, metoprolol, midazolam, omeprazole, rosiglitazone and theophylline). The final model utilised tissue concentrations from adipose, skin, muscle and a contribution from artery. Predicted C_max values considering the central venous compartment can over-predict the observed values up to 10-fold whereas the new sampling site predictions were within 2-fold of observed values. The model was particularly relevant for studies where traditional PBPK models over-predict early time point concentrations. A successful corrective model for C_max prediction has been developed, subject to further validation. These models can be enrolled as built-up modules into PBPK platforms and potentially account for factors that may affect the initial mixing of the blood at the site of sampling.

Effect of multiple intravenous doses of lanicemine (AZD6765) on the pharmacokinetics of midazolam in healthy subjects

The objectives of the present study were to evaluate safety and tolerability as well as the effects of multiple doses of lanicemine on the pharmacokinetics of a CYP3A substrate, midazolam. A total of 46 healthy volunteers were enrolled in the open-label, fixed-sequence, nonrandomized study. All volunteers received an oral dose of 5 mg of midazolam alone or after 6 days of 150 mg daily intravenous infusion of lanicemine. Lanicemine reached a plasma Cmax of 1.51 μg/mL after 150 mg daily dosing to steady state. The geometric mean CL, Vss, and t1/2 of lanicemine were 8.1 L/h, 122.0 L, and 10.4 hours, respectively. The geometric least-squares mean ratios and 90% confidence intervals for midazolam AUC0- ∞ , and Cmax were within the 80% to 125% limits when lanicemine plus midazolam treatment was compared with midazolam alone, demonstrating that daily dosing with 150 mg of lanicemine for 6 days had no effect on CYP3A activity. Comprehensive physiologically based pharmacokinetic modeling using in vitro and in silico findings also indicated lanicemine would have little impact on the pharmacokinetics of CYP3A substrate, such as midazolam. In addition, lanicemine and midazolam administered alone or in combination were generally safe and well tolerated.

Rationalizing underprediction of drug clearance from enzyme and transporter kinetic data: from in vitro tools to mechanistic modeling

Over the years, there has been an increase in the number and quality of available in vitro tools for the assessment of clearance. Complexity of data analysis and modelling of corresponding in vitro data has increased in an analogous manner, in particular for the simultaneous characterization of transporter and metabolism kinetics, together with intracellular binding and passive diffusion. In the current chapter, the impact of different factors on the in vitro-in vivo extrapolation of clearance will be addressed in a stepwise manner, from the selection of the most adequate in vitro system and experimental design/condition to the corresponding modelling of data generated. The application of static or physiologically based pharmacokinetic models in the prediction of clearance will be discussed, highlighting limitations and current challenges of some of the approaches. Particular focus will be on the ability of in vitro and in silico predictive tools to overcome the trend of clearance underprediction. Improvements made as a result of inclusion of extrahepatic metabolism and consideration of transporter-metabolism interplay across different organs will be discussed.

Ubiquity: a framework for physiological/mechanism-based pharmacokinetic/pharmacodynamic model development and deployment

Practitioners of pharmacokinetic/pharmacodynamic modeling routinely employ various software packages that enable them to fit differential equation based mechanistic or empirical models to biological/pharmacological data. The availability and choice of different analytical tools, while enabling, can also pose a significant challenge in terms of both, implementation and transferability. A package has been developed that addresses these issues by creating a simple text-based format, which provides methods to reduce coding complexity and enables the modeler to describe the components of the model based on the underlying physiochemical processes. A Perl script builds the system for multiple formats (ADAPT, MATLAB, Berkeley Madonna, etc.), enabling analysis across several software packages and reducing the chance for transcription error. Workflows can then be built around this package, which can increase efficiency and model availability. As a proof of concept, tools are included that allow models constructed in this format to be run with MATLAB both at the scripting level and through a generic graphical application that can be compiled and run as a stand-alone application.

Impact of physiologically based pharmacokinetic modeling and simulation in drug development

Physiologically based pharmacokinetic modeling and simulation can be used to predict the pharmacokinetics of drugs in human populations and to explore the effects of varying physiologic parameters that result from aging, ethnicity, or disease. In addition, the effects of concomitant medications on drug exposure can be investigated; prediction of the magnitude of drug interactions can impact regulatory communications or internal decision-making regarding the requirement for a clinical drug interaction study. Modeling and simulation can also help to inform the design and timings of clinical drug interaction studies, resulting in more efficient use of limited resources and improved planning in addition to promoting mechanistic understanding of observed drug interactions. These approaches have been used in GlaxoSmithKline from drug discovery to registration and have been applied to 41 drugs from a number of therapeutic areas. This report highlights the variety of questions that can be addressed by prospective or retrospective application of modeling and simulation and the impact this can have on clinical drug development (from candidate selection through clinical development to regulatory submissions).

Drug metabolites as cytochrome p450 inhibitors: a retrospective analysis and proposed algorithm for evaluation of the pharmacokinetic interaction potential of metabolites in drug discovery and development

Understanding drug-drug interactions (DDIs) is a key component of clinical practice ensuring patient safety and efficacy of medicines. The role of drug metabolites in DDIs is a developing area of science, and has been recently highlighted in a draft regulatory guidance. The guidance states that metabolites representing ≥25% of the parent drug's area under the plasma concentration/time curve and/or >10% of exposure of total drug-related material should trigger in vitro characterization of metabolites for cytochrome P450 inhibition and propensity for DDIs. The relationship between in vitro cytochrome P450 inhibitory potency, systemic exposure, and DDI potential of drug metabolites was examined using the Pfizer development database to identify compounds with pre-existing in vivo biotransformation data, where circulating metabolites were identified in humans. The database yielded 33 structurally diverse compounds with collectively 115 distinct circulating metabolites. Of these, 52% (60/115) achieved exposures >25% of parent drug levels as judged from mass balance/metabolite identification studies. It was noted that 14 metabolite standards for 12 parent drugs had been synthesized, monitored in clinical studies, and examined for cytochrome P450 inhibition. For the 14 metabolite/parent drug pairs, no clinically relevant DDIs were expected to occur against the major human cytochrome P450 isoforms. A review of the literature for parent/metabolite DDI information was also conducted to examine trends using a larger data set. Leveraging the analysis of both internal and literature-based data sets, an algorithm was devised for use in drug discovery/early development to assess cytochrome P450 inhibitory potential of drug metabolites and the propensity to cause a clinically relevant DDI.