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

Applications of model simulation in pharmacological fields and the problems of theoretical reliability

The use of mathematical models has become increasingly prevalent in pharmacological fields, particularly in drug development processes. These models are instrumental in tasks such as designing clinical trials and assessing factors like efficacy, toxicity, and clinical practice. Various types of models have been developed and documented. Nevertheless, emphasizing the reliability of parameter values is crucial, as they play a pivotal role in shaping the behavior of the system. In some instances, parameter values reported previously are treated as fixed values, which can lead to convergence towards values that deviate substantially from those found in actual biological systems. This is especially true when parameter values are determined through fitting to limited observations. To mitigate this risk, the reuse of parameter values from previous reports should be approached with a critical evaluation of their validity. Currently, there is a proposal for a simultaneous search for plausible values for all parameters using comprehensive search algorithms in both pharmacokinetic and pharmacodynamic or systems pharmacological models. Implementing these methodologies can help address issues related to parameter determination. Furthermore, integrating these approaches with methods developed in the field of machine-learning field has the potential to enhance the reliability of parameter values and the resulting model outputs.

Modeling the complexity of drug-drug interactions: A physiologically-based pharmacokinetic study of Lenvatinib with Schisantherin A/Schisandrin A

BACKGROUND

Lenvatinib's efficacy as a frontline targeted therapy for radioactive iodine-refractory thyroid carcinoma and advanced hepatocellular carcinoma owes to its inhibition of multiple tyrosine kinases. However, as a CYP3A4 substrate, lenvatinib bears susceptibility to pharmacokinetic modulation by co-administered agents. Schisantherin A (STA) and schisandrin A (SIA) - bioactive lignans abundant in the traditional Chinese medicinal Wuzhi Capsule - act as CYP3A4 inhibitors, engendering the potential for drug-drug interactions (DDIs) with lenvatinib.

METHODS

To explore potential DDIs between lenvatinib and STA/SIA, we developed a physiologically-based pharmacokinetic (PBPK) model for lenvatinib and used it to construct a DDI model for lenvatinib and STA/SIA. The model was validated with clinical trial data and used to predict changes in lenvatinib exposure with combined treatment.

RESULTS

Following single-dose administration, the predicted area under the plasma concentration-time curve (AUC) and maximum plasma concentrations (Cmax) of lenvatinib increased 1.00- to 1.03-fold and 1.00- to 1.01-fold, respectively, in the presence of STA/SIA. Simulations of multiple-dose regimens revealed slightly greater interactions, with lenvatinib AUC0-t and Cmax increasing up to 1.09-fold and 1.02-fold, respectively.

CONCLUSION

Our study developed the first PBPK and DDI models for lenvatinib as a victim drug. STA and SIA slightly increased lenvatinib exposure in simulations, providing clinically valuable information on the safety of concurrent use. Given the minimal pharmacokinetic changes, STA/SIA are unlikely to interact with lenvatinib through pharmacokinetic alterations synergistically but rather may enhance efficacy through inherent anti-cancer efficacy of STA/ SIA.

Physiologically-based pharmacokinetic modeling of quinidine to establish a CYP3A4, P-gp, and CYP2D6 drug-drug-gene interaction network

The antiarrhythmic agent quinidine is a potent inhibitor of cytochrome P450 (CYP) 2D6 and P-glycoprotein (P-gp) and is therefore recommended for use in clinical drug-drug interaction (DDI) studies. However, as quinidine is also a substrate of CYP3A4 and P-gp, it is susceptible to DDIs involving these proteins. Physiologically-based pharmacokinetic (PBPK) modeling can help to mechanistically assess the absorption, distribution, metabolism, and excretion processes of a drug and has proven its usefulness in predicting even complex interaction scenarios. The objectives of the presented work were to develop a PBPK model of quinidine and to integrate the model into a comprehensive drug-drug(-gene) interaction (DD(G)I) network with a diverse set of CYP3A4 and P-gp perpetrators as well as CYP2D6 and P-gp victims. The quinidine parent-metabolite model including 3-hydroxyquinidine was developed using pharmacokinetic profiles from clinical studies after intravenous and oral administration covering a broad dosing range (0.1-600 mg). The model covers efflux transport via P-gp and metabolic transformation to either 3-hydroxyquinidine or unspecified metabolites via CYP3A4. The 3-hydroxyquinidine model includes further metabolism by CYP3A4 as well as an unspecific hepatic clearance. Model performance was assessed graphically and quantitatively with greater than 90% of predicted pharmacokinetic parameters within two-fold of corresponding observed values. The model was successfully used to simulate various DD(G)I scenarios with greater than 90% of predicted DD(G)I pharmacokinetic parameter ratios within two-fold prediction success limits. The presented network will be provided to the research community and can be extended to include further perpetrators, victims, and targets, to support investigations of DD(G)Is.

Experimental and computational models to investigate intestinal drug permeability and metabolism

Oral administration is the preferred route for drug administration that leads to better therapy compliance. The intestine plays a key role in the absorption and metabolism of oral drugs, therefore, new intestinal models are being continuously proposed, which contribute to the study of intestinal physiology, drug screening, drug side effects, and drug-drug interactions.Advances in pharmaceutical processes have produced more drug formulations, causing challenges for intestinal models. To adapt to the rapid evolution of pharmaceuticals, more intestinal models have been created. However, because of the complexity of the intestine, few models can take all aspects of the intestine into account, and some functions must be sacrificed to investigate other areas. Therefore, investigators need to choose appropriate models according to the experimental stage and other requirements to obtain the desired results.To help researchers achieve this goal, this review summarised the advantages and disadvantages of current commonly used intestinal models and discusses possible future directions, providing a better understanding of intestinal models.

Variations in the Frequencies of Polymorphisms in the CYP450s Genes in Eight Major Ethnicities of Iran: A Review of the Human Data

Genetic polymorphisms in cytochrome P450 genes can cause variation in metabolism. Thus, single nucleotide variants significantly impact drug pharmacokinetics, toxicity factors, and efficacy and safety of medicines. The distribution of CYP450 alleles varies drastically across ethnicities, with significant implications for personalized medicine and the healthcare system. We combined whole-genome and exome sequencing data to provide a review of CYP450 allele polymorphisms with clinical importance. Data were collected from 800 unrelated Iranians (100 subjects from 8 major ethnicities of Iran), more than 32,000 unrelated Europeans (other than Caucasian), and four Middle Eastern countries. We analyzed the frequencies and similarities of 17 CYP450 frequent alleles related to nine important CYP450 isoenzymes and homozygous and heterozygous genotypes based on these alleles in eight major Iranian ethnics by integrating these data with population-specific linkage information and compared these datasets with mentioned populations.

Candidate Proficiency Test Chemicals to Address Industrial Chemical Applicability Domains for in vitro Human Cytochrome P450 Enzyme Induction

Cytochrome P450 (CYP) enzymes play a key role in the metabolism of both xenobiotics and endogenous chemicals, and the activity of some CYP isoforms are susceptible to induction and/or inhibition by certain chemicals. As CYP induction/inhibition can bring about significant alterations in the level of in vivo exposure to CYP substrates and metabolites, CYP induction/inhibition data is needed for regulatory chemical toxicity hazard assessment. On the basis of available human in vivo pharmaceutical data, a draft Organisation for Economic Co-operation and Development Test Guideline (TG) for an in vitro CYP HepaRG test method that is capable of detecting the induction of four human CYPs (CYP1A1/1A2, 2B6, and 3A4), has been developed and validated for a set of pharmaceutical proficiency chemicals. However to support TG adoption, further validation data was requested to demonstrate the ability of the test method to also accurately detect CYP induction mediated by industrial and pesticidal chemicals, together with an indication on regulatory uses of the test method. As part of "GOLIATH", a European Union Horizon-2020 funded research project on metabolic disrupting chemical testing approaches, work is underway to generate supplemental validated data for an additional set of chemicals with sufficient diversity to allow for the approval of the guideline. Here we report on the process of proficiency chemical selection based on a targeted literature review, the selection criteria and considerations required for acceptance of proficiency chemical selection for OECD TG development (i.e. structural diversity, range of activity, relevant chemical sectors, global restrictions etc). The following 13 proposed proficiency chemicals were reviewed and selected as a suitable set for use in the additional validation experiments: tebuconazole, benfuracarb, atrazine, cypermethrin, chlorpyrifos, perfluorooctanoic acid, bisphenol A, N,N-diethyl-m-toluamide, benzo-[a]-pyrene, fludioxonil, malathion, triclosan, and caffeine. Illustrations of applications of the test method in relation to endocrine disruption and non-genotoxic carcinogenicity are provided.

Prediction of CYP‐mediated DDIs involving inhibition: Approaches to address the requirements for system qualification of the Simcyp Simulator

Physiologically‐based pharmacokinetic (PBPK) modeling is being increasingly used in drug development to avoid unnecessary clinical drug–drug interaction (DDI) studies and inform drug labels. Thus, regulatory agencies are recommending, or indeed requesting, more rigorous demonstration of the prediction accuracy of PBPK platforms in the area of their intended use. We describe a framework for qualification of the Simcyp Simulator with respect to competitive and mechanism‐based inhibition (MBI) of CYP1A2, CYP2D6, CYP2C8, CYP2C9, CYP2C19, and CYP3A4/5. Initially, a DDI matrix, consisting of a range of weak, moderate, and strong inhibitors and substrates with varying fraction metabolized by specific CYP enzymes that were susceptible to different degrees of inhibition, were identified. Simulations were run with 123 clinical DDI studies involving competitive inhibition and 78 clinical DDI studies involving MBI. For competitive inhibition, the overall prediction accuracy was good with an average fold error (AFE) of 0.91 and 0.92 for changes in the maximum plasma concentration (C_max) and area under the plasma concentration (AUC) time profile, respectively, as a consequence of the DDI. For MBI, an AFE of 1.03 was determined for both C_max and AUC. The prediction accuracy was generally comparable across all CYP enzymes, irrespective of the isozyme and mechanism of inhibition. These findings provide confidence in application of the Simcyp Simulator (V19 R1) for assessment of the DDI potential of drugs in development either as inhibitors or victim drugs of CYP‐mediated interactions. The approach described herein and the identified DDI matrix can be used to qualify subsequent versions of the platform.

Applications, Challenges, and Outlook for PBPK Modeling and Simulation: A Regulatory, Industrial and Academic Perspective

Several regulatory guidances on the use of physiologically based pharmacokinetic (PBPK) analyses and physiologically based biopharmaceutics model(s) (PBBM(s)) have been issued. Workshops are routinely held, demonstrating substantial interest in applying these modeling approaches to address scientific questions in drug development. PBPK models and PBBMs have remarkably contributed to model-informed drug development (MIDD) such as anticipating clinical PK outcomes affected by extrinsic and intrinsic factors in general and specific populations. In this review, we proposed practical considerations for a "base" PBPK model construction and development, summarized current status, challenges including model validation and gaps in system models, and future perspectives in PBPK evaluation to assess a) drug metabolizing enzyme(s)- or drug transporter(s)- mediated drug-drug interactions b) dosing regimen prediction, sampling timepoint selection and dose validation in pediatric patients from newborns to adolescents, c) drug exposure in patients with renal and/or and hepatic organ impairment, d) maternal-fetal drug disposition during pregnancy, and e) pH-mediated drug-drug interactions in patients treated with proton pump inhibitors/acid-reducing agents (PPIs/ARAs) intended for gastric protection. Since PBPK can simulate outcomes in clinical studies with enrollment challenges or ethical issues, the impact of PBPK models on waivers and how to strengthen study waiver is discussed.

Toxicity induced by ciprofloxacin and enrofloxacin: oxidative stress and metabolism

Ciprofloxacin (CIP) (human use) and enrofloxacin (ENR) (veterinary use) are synthetic anti-infectious medications that belong to the second generation of fluoroquinolones. They have a wide antimicrobial spectrum and strong bactericidal effects at very low concentrations via enzymatic inhibition of DNA gyrase and topoisomerase IV, which are required for DNA replication. They also have high bioavailability, rapid absorption with favorable pharmacokinetics and excellent tissue penetration, including cerebral spinal fluid. These features have made them the most applied antibiotics in both human and veterinary medicine. ENR is marketed exclusively for animal medicine and has been widely used as a therapeutic veterinary antibiotic, resulting in its residue in edible tissues and aquatic environments, as well as the development of resistance and toxicity. Estimation of the risks to humans due to antimicrobial resistance produced by CIP and ENR is important and of great interest. Moreover, in rare cases due to their overdose and/or prolonged administration, the development of CIP and ENR toxicity may occur. The toxicity of these fluoroquinolones antimicrobials is mainly related to reactive oxygen species (ROS) and oxidative stress (OS) generation, besides metabolism-related toxicity. Therefore, CIP is restricted in pregnant and lactating women, pediatrics and elderly similarly ENR do in the veterinary field. This review manuscript aims to identify the toxicity induced by ROS and OS as a common sequel of CIP and ENR. Furthermore, their metabolism and the role of metabolizing enzymes were reported.

Antifungal Drugs TDM: Trends and Update

PURPOSE

The increasing burden of invasive fungal infections results in growing challenges to antifungal (AF) therapeutic drug monitoring (TDM). This review aims to provide an overview of recent advances in AF TDM.

METHODS

We conducted a PubMed search for articles during 2016-2020 using "TDM" or "pharmacokinetics" or "drug-drug-interaction" with "antifungal," consolidated for each AF. Selection was limited to English language articles with human data on drug exposure.

RESULTS

More than 1000 articles matched the search terms. We selected 566 publications. The latest findings tend to confirm previous observations in real-life clinical settings. The pharmacokinetic variability related to special populations is not specific but must be considered. AF benefit-to-risk ratio, drug-drug interaction (DDI) profiles, and minimal inhibitory concentrations for pathogens must be known to manage at-risk situations and patients. Itraconazole has replaced ketoconazole in healthy volunteers DDI studies. Physiologically based pharmacokinetic modeling is widely used to assess metabolic azole DDI. AF prophylactic use was studied more for Aspergillus spp. and Mucorales in oncohematology and solid organ transplantation than for Candida (already studied). Emergence of central nervous system infection and severe infections in immunocompetent individuals both merit special attention. TDM is more challenging for azoles than amphotericin B and echinocandins. Fewer TDM requirements exist for fluconazole and isavuconazole (ISZ); however, ISZ is frequently used in clinical situations in which TDM is recommended. Voriconazole remains the most challenging of the AF, with toxicity limiting high-dose treatments. Moreover, alternative treatments (posaconazole tablets, ISZ) are now available.

CONCLUSIONS

TDM seems to be crucial for curative and/or long-term maintenance treatment in highly variable patients. TDM poses fewer cost issues than the drugs themselves or subsequent treatment issues. The integration of clinical pharmacology into multidisciplinary management is now increasingly seen as a part of patient care.

Predictive In Vitro-In Vivo Extrapolation for Time Dependent Inhibition of CYP1A2, CYP2C8, CYP2C9, CYP2C19 and CYP2D6 Using Pooled Human Hepatocytes, Human Liver Microsomes, and a Simple Mechanistic Static Model

Inactivation of Cytochrome P450 (CYP450) enzymes can lead to significant increases in exposure of co-medicants. The majority of reported in vitro to in vivo extrapolation (IVIVE) data have historically focused on CYP3A4 leaving the assessment of other CYP isoforms insubstantial. To this end, the utility of human hepatocytes (HHEP) and microsome (HLM) to predict clinically relevant DDIs was investigated with a focus on CYP1A2, CYP2C8, CYP2C9, CYP2C19 and CYP2D6. Evaluation of IVIVE for CYP2B6 was limited to only weak inhibition. A search of the University of Washington Drug-Drug Interaction Database was conducted to identify a clinically relevant weak, moderate and strong inhibitor for selective substrates of CYP1A2, CYP2C8, CYP2C9, CYP2C19 and CYP2D6, resulting in 18 inhibitors for in vitro characterization against 119 clinical interaction studies. Pooled human hepatocytes and HLM were pre-incubated with increasing concentrations of inhibitors for designated timepoints. Time dependent inhibition (TDI) was detected in HLM for four moderate/strong inhibitors suggesting that some optimization of incubation conditions (i.e. lower protein concentrations) is needed to capture weak inhibition. Clinical risk assessment was conducted by incorporating the in vitro derived kinetic parameters k_inact and K_I into static equations recommended by regulatory authorities. Significant overprediction was observed when applying the basic models recommended by regulatory agencies. Mechanistic static models (MSM), which consider the fraction of metabolism through the impacted enzyme, using the unbound hepatic inlet concentration lead to the best overall prediction accuracy with 92% and 85% of data from HHEPs and HLM, respectively, within 2-fold of the observed value. Significance Statement Collectively, the data demonstrate that coupling time-dependent inactivation parameters derived from pooled human hepatocytes and HLM with a mechanistic static model provides an easy and quantitatively accurate means to determine clinical DDI risk from in vitro data. Weak and moderate inhibitors did not show TDI under standard incubation conditions using HLM and optimization of incubation conditions is warranted. Recommendations are made with respect to input parameters for IVIVE of TDI with non-CYP3A enzymes using available data from HLM and HHEPs.

Prediction of Drug-Drug Interaction Potential of Tegoprazan Using Physiologically Based Pharmacokinetic Modeling and Simulation

This study aimed to develop a physiologically based pharmacokinetic (PBPK) model of tegoprazan and to predict the drug-drug interaction (DDI) potential between tegoprazan and cytochrome P450 (CYP) 3A4 perpetrators. The PBPK model of tegoprazan was developed using SimCYP Simulator^® and verified by comparing the model-predicted pharmacokinetics (PKs) of tegoprazan with the observed data from phase 1 clinical studies, including DDI studies. DDIs between tegoprazan and three CYP3A4 perpetrators were predicted by simulating the difference in tegoprazan exposure with and without perpetrators, after multiple dosing for a clinically used dose range. The final PBPK model adequately predicted the biphasic distribution profiles of tegoprazan and DDI between tegoprazan and clarithromycin. All ratios of the predicted-to-observed PK parameters were between 0.5 and 2.0. In DDI simulation, systemic exposure to tegoprazan was expected to increase about threefold when co-administered with the maximum recommended dose of clarithromycin or ketoconazole. Meanwhile, tegoprazan exposure was expected to decrease to ~30% when rifampicin was co-administered. Based on the simulation by the PBPK model, it is suggested that the DDI potential be considered when tegoprazan is used with CYP3A4 perpetrator, as the acid suppression effect of tegoprazan is known to be associated with systemic exposure.

Prediction of Fluoxetine and Norfluoxetine Pharmacokinetic Profiles Using Physiologically Based Pharmacokinetic Modeling

Fluoxetine is a selective serotonin reuptake inhibitor that is metabolized to norfluoxetine by cytochrome P450 (CYP) 2D6, CYP2C19, and CYP3A4. A physiologically based pharmacokinetic model for fluoxetine and norfluoxetine metabolism was developed to predict and investigate changes in concentration-time profiles according to fluoxetine dosage in the Korean population. The model was developed based on the Certara repository model and information gleaned from the literature. Digitally extracted clinical study data were used to develop and verify the model. Simulations for plasma concentrations of fluoxetine and norfluoxetine after a single dose of 60 or 80 mg fluoxetine were made based on 1000 virtual healthy Korean individuals using the SimCYP version 19 simulator. The mean ratios (simulated/observed) after a single administration of 80 mg fluoxetine for maximum plasma concentration, area under the plasma concentration-time curve, and apparent clearance were 1.12, 1.08, and 0.93 for fluoxetine; the ratios of maximum plasma concentration and area under the plasma concentration-time curve were 1.08 and 1.08, respectively, for norfluoxetine, indicating that the simulated concentration-time profiles of fluoxetine and norfluoxetine fitted the observed profiles well. The developed model was used to predict plasma fluoxetine and norfluoxetine concentration-time profiles after repeated administrations of fluoxetine in Korean volunteers. This physiologically based pharmacokinetic model could provide basic understanding of the pharmacokinetic profiles of fluoxetine and its metabolite under various situations.

A generic framework for the physiologically-based pharmacokinetic platform qualification of PK-Sim and its application to predicting cytochrome P450 3A4-mediated drug-drug interactions

The success of applications of physiologically‐based pharmacokinetic (PBPK) modeling in drug development and drug labeling has triggered regulatory agencies to demand rigorous demonstration of the predictive capability of the specific PBPK platform for a particular intended application purpose. The effort needed to comply with such qualification requirements exceeds the costs for any individual PBPK application. Because changes or updates of a PBPK platform would require (re‐)qualification, a reliable and efficient generic qualification framework is needed. We describe the development and implementation of an agile and sustainable technical framework for automatic PBPK platform (re‐)qualification of PK‐Sim^® embedded in the open source and open science GitHub landscape of Open Systems Pharmacology. The qualification approach enables the efficient assessment of all aspects relevant to the qualification of a particular purpose and provides transparency and traceability for all stakeholders. As a showcase example for the power and versatility of the qualification framework, we present the qualification of PK‐Sim^® for the intended purpose of predicting cytochrome P450 3A4 (CYP3A4)–mediated drug–drug interactions (DDIs). Several perpetrator PBPK models featuring various degrees of CYP3A4 modulation and different types of mechanisms (competitive inhibition, mechanism‐based inactivation, and induction) were coupled with a set of PBPK models of sensitive CYP3A4 victim drugs. Simulations were compared to a comprehensive data set of 135 observations from published clinical DDI studies. The platform's overall predictive performance showed reasonable accuracy and precision (geometric mean fold error of 1.4 for both area under the plasma concentration‐time curve ratios and peak plasma concentration ratios with/without perpetrator) and suggests that PK‐Sim^® can be applied to quantitatively assess CYP3A4‐mediated DDI in clinically untested scenarios.

The Prevalence of Selected Potential Drug-Drug Interactions of Analgesic Drugs and Possible Methods of Preventing Them: Lessons Learned From the Analysis of the Real-World National Database of 38 Million Citizens of Poland

Introduction: Drug-drug interactions may lead to poor health outcomes, as well as increased costs and utilization of healthcare services. Unfortunately, real-world data continuously prove high prevalence of potential drug-drug interactions (pDDIs) worldwide. Among identified drivers, ageing, multimorbidity and polypharmacy play a very important role. With these factors being widespread, the need for implementation of strategies minimizing the burden of pDDIs becomes an urgency. This, however, requires a better understanding of the prevalence of pDDIs and the underlying causative factors.

Aim of study: To assess the real-world prevalence of pDDIs and its characteristics in the general population of Poland, using analgesic drugs as a model, and to find out whether pDDIs are caused by prescribing coming from the very same prescribers (co-prescribing).

Methods: A retrospective analysis of the 2018 dispensation data of the National Health Fund (NHF) - the only Polish public healthcare payer organization with nationwide coverage. We searched for selected pDDIs of non-steroidal anti-inflammatory drugs (NSAIDs) with antihypertensives, other NSAIDs (double use), oral glucocorticoids, oral anticoagulants, selective serotonin reuptake inhibitors (SSRIs), serotonin–norepinephrine reuptake inhibitors (SNRIs), and antiplatelet drugs; as well as opioides with SSRIs, SNRIs, gabapentinoids, and benzodiazepines. A pDDI was deemed present if two drugs standing in a possible conflict were dispensed within the same calendar month.

Results: Out of 38.4 million citizens of Poland, 23.3 million were dispensed prescribed drugs reimbursed by NHF in 2018. In this cohort, we have identified 2,485,787 cases of analgesic drug pDDIs, corresponding with 6.47% of the Polish population. Out of these, the most prevalent pDDI was caused by “NSAIDs + antihypertensives” (1,583,575 cases, i.e., 4.12% of the Polish population), followed by “NSAIDs + NSAIDs” (538,640, 1.40%) and “NSAIDs + glucocorticoids” (213,504, 0.56%). The most persistent pDDIs among those studied were caused by “Opioids + Gabapentinoids” (2.19, 95%CI: 2.16–2.22 months). On average, 76.63% of all cases of pDDIs were caused by drugs prescribed by the very same prescribers.

Conclusion: Based on high-quality, nationwide data, we have found a high prevalence of analgesic drugs-related pDDIs in Poland. Over ¾ of the identified pDDIs were caused by co-prescribing, i.e., prescriptions issued by the same prescribers. The significance of the problem, illustrated with our findings on analgesic drugs-related pDDIs in Poland, deserves much more scientific and policymaker attention.

A quantitative model for metabolic intervention using gut microbes

As medicine shifts toward precision-based and personalized therapeutics, utilizing more complex biomolecules to treat increasingly difficult and rare conditions, microorganisms provide an avenue for realizing the production and processing necessary for novel drug pipelines. More so, probiotic microbes can be co-opted to deliver therapeutics by oral administration as living drugs, able to survive and safely transit the digestive tract. As living therapeutics are in their nascency, traditional pharmacokinetic-pharmacodynamic (PK-PD) models for evaluating drug candidates are not appropriate for this novel platform. Using a living therapeutic in late-stage clinical development for phenylketonuria (PKU) as a case study, we adapt traditional oral drug delivery models to properly evaluate and inform the engineering of living therapeutics. We develop the adapted for living therapeutics compartmental absorption and transit (ALT-CAT) model to provide metrics for drug efficacy across nine age groups of PKU patients and evaluate model parameters that are influenced by patient physiology, microbe selection and therapeutic production, and dosing formulations. In particular, the ALT-CAT model describes the mathematical framework to model the behavior of orally delivered engineered bacteria that act as living therapeutics by adapting similar methods that have been developed and widely-used for small molecular drug delivery and absorption.

Physiologically based pharmacokinetic modeling of altered tizanidine systemic exposure by CYP1A2 modulation: Impact of drug-drug interactions and cigarette consumption

Tizanidine is an alpha2-adrenergic agonist, used to treat spasticity associated with multiple sclerosis and spinal injury. Tizanidine is primarily metabolized by CYP1A2 and is considered a sensitive index substrate for this enzyme. The physiologically based pharmacokinetic (PBPK) modeling platform Simcyp® was used to evaluate the impact of CYP1A2 modulation on tizanidine exposure through drug-drug interactions (DDIs) and host-dependent habits (cigarette smoking). A PBPK model was developed to predict tizanidine disposition in healthy volunteers following oral administration. The model was verified based on agreement between model-simulated and clinically observed systemic exposure metrics (C_max, AUC). The model was then used to carry-out DDI simulations to predict alterations in tizanidine systemic exposure when co-administered with various CYP1A2 perpetrators including competitive inhibitors (fluvoxamine, ciprofloxacin), a mechanism-based inhibitor (rofecoxib), and an inducer (rifampin). Additional simulations were performed to evaluate the impact of cigarette smoking on systemic exposure. Under each scenario, the PBPK model was able to capture the observed fold changes in tizanidine C_max and AUC of tizanidine when coadministered with CYP1A2 inhibitors or inducers. These results add to the available research findings in the literature on PBPK predictions of drug-drug interactions and illustrate the potential application in drug development, specifically to support product labeling.

Physiologically-Based Pharmacokinetic Modeling for the Prediction of a Drug-Drug Interaction of Combined Effects on P-glycoprotein and Cytochrome P450 3A

Direct oral anticoagulants, such as apixaban and rivaroxaban, are important for the treatment and prophylaxis of venous thromboembolism and to reduce the risk of stroke and systemic embolism in patients with nonvalvular atrial fibrillation. Because apixaban and rivaroxaban are predominantly eliminated by cytochrome P450 (CYP) 3A and P‐glycoprotein (P‐gp), concomitant use of combined P‐gp and strong CYP3A4 inhibitors and inducers should be avoided. Physiologically‐based pharmacokinetic models for apixaban and rivaroxaban were developed to estimate the net effect of CYP3A induction, P‐gp inhibition, and P‐gp induction by rifampicin. The disposition of rivaroxaban is more complex compared with apixaban because both hepatic and renal P‐gp is considered to contribute to rivaroxaban elimination. Furthermore, organic anion transporter‐3, a renal uptake transporter, may also contribute the elimination of rivaroxaban from systemic circulation. The models were verified with observed clinical drug–drug interactions with CYP3A and P‐gp inhibitors. With the developed models, the predicted area under the concentration time curve and maximum concentration ratios were 0.43 and 0.48, respectively, for apixaban, and 0.50–0.52 and 0.72–0.73, respectively, for rivaroxaban when coadministered with 600 mg multiple doses of rifampicin and that were very close to observed data. The impact of each of the elimination pathways was assessed for rivaroxaban, and inhibition of CYP3A led to a larger impact over intestinal and hepatic P‐gp. Inhibition of renal organic anion transporter‐3 or P‐gp led to an overall modest interaction. The developed apixaban and rivaroxaban models can be further applied to the investigation of interactions with other P‐gp and/or CYP3A4 inhibitors and inducers.

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.

Recommendations for the Design of Clinical Drug-Drug Interaction Studies With Itraconazole Using a Mechanistic Physiologically-Based Pharmacokinetic Model

Regulatory agencies currently recommend itraconazole (ITZ) as a strong cytochrome P450 3A (CYP3A) inhibitor for clinical drug–drug interaction (DDI) studies. This work by an International Consortium for Innovation and Quality in Pharmaceutical Development working group (WG) is to develop and verify a mechanistic ITZ physiologically‐based pharmacokinetic model and provide recommendations for optimal DDI study design based on model simulations. To support model development and verification, in vitro and clinical PK data for ITZ and its metabolites were collected from WG member companies. The model predictions of ITZ DDIs with seven different CYP3A substrates were within the guest criteria for 92% of area under the concentration‐time curve ratios and 95% of maximum plasma concentration ratios, thus verifying the model for DDI predictions. The verified model was used to simulate various clinical DDI study scenarios considering formulation, duration of dosing, dose regimen, and food status to recommend the optimal design for maximal inhibitory effect by ITZ.

The Pharmacokinetic Prediction of Cyclosporin A after Coadministration with Wuzhi Capsule

We aim to describe the influence of principal ingredients of Wuzhi capsule, schisandrin A (SIA) and schisantherin A (STA), on the pharmacokinetics of cyclosporin A (CsA) and to quantify the herb-drug interactions (HDIs) between SIA, STA, and CsA. CsA is a first-line immunosuppressant for anti-rejection therapy after solid organ transplantation, while narrow therapeutic window associated with strong hepatotoxicity largely limited its use. Wuzhi capsule, a liver-protective drug, was approved for coadministration with CsA to reduce the hepatotoxicity. There are few studies exploring HDIs of CsA when coadministered with Wuzhi capsule. The essential adjusted physicochemical data and pharmacokinetic parameters of SIA, STA, and CsA were collected. Then physiologically based pharmacokinetic (PBPK) models of SIA, STA, and CsA were built and verified in healthy subjects using Simcyp respectively. The refined PBPK models were used to estimate potential HDIs between CsA and SIA, STA. The simulated plasma concentration-time curves of CsA, SIA, and STA were in good accordance with the observed profiles respectively. CsA pharmacokinetics were improved after coadministration. After a single dose and multiple doses, the area under the plasma concentration-time curve (AUC) of CsA was increased by 47% and 226% when coadministered with STA, respectively, and by 8% and 36% when coadministered with SIA, respectively. PBPK models sufficiently described the pharmacokinetics of CsA, SIA, and STA. Compared with SIA, STA inhibited CsA metabolism to a greater extent. Our result revealed the dose of CsA can be reduced to maintain similar profile when used concomitantly with Wuzhi capsule.

Physiologically-Based Pharmacokinetic Modeling for the Prediction of CYP2D6-Mediated Gene-Drug-Drug Interactions

The aim of this work was to predict the extent of Cytochrome P450 2D6 (CYP2D6)‐mediated drug–drug interactions (DDIs) in different CYP2D6 genotypes using physiologically‐based pharmacokinetic (PBPK) modeling. Following the development of a new duloxetine model and optimization of a paroxetine model, the effect of genetic polymorphisms on CYP2D6‐mediated intrinsic clearances of dextromethorphan, duloxetine, and paroxetine was estimated from rich pharmacokinetic profiles in activity score (AS)1 and AS2 subjects. We obtained good predictions for the dextromethorphan–duloxetine interaction (Ratio of predicted over observed area under the curve (AUC) ratio (R_pred/obs) 1.38–1.43). Similarly, the effect of genotype was well predicted, with an increase of area under the curve ratio of 28% in AS2 subjects when compared with AS1 (observed, 33%). Despite an approximately twofold underprediction of the dextromethorphan–paroxetine interaction, an R_pred/obs of 0.71 was obtained for the effect of genotype on the area under the curve ratio. Therefore, PBPK modeling can be successfully used to predict gene–drug–drug interactions (GDDIs). Based on these promising results, a workflow is suggested for the generic evaluation of GDDIs and DDIs that can be applied in other situations.

In silico models in drug development: where we are

The use and utility of computational models in drug development has significantly grown in the last decades, fostered by the availability of high throughput datasets and new data analysis strategies. These in silico approaches are demonstrating their ability to generate reliable predictions as well as new knowledge on the mode of action of drugs and the mechanisms underlying their side effects, altogether helping to reduce the costs of drug development. The aim of this review is to provide a panorama of developments in the field in the last two years.

Cryopreserved Human Intestinal Mucosal Epithelium: A Novel In Vitro Experimental System for the Evaluation of Enteric Drug Metabolism, Cytochrome P450 Induction, and Enterotoxicity

We report here a novel in vitro enteric experimental system, cryopreserved human intestinal mucosa (CHIM), for the evaluation of enteric drug metabolism, drug-drug interaction, drug toxicity, and pharmacology. CHIM was isolated from the small intestines of four human donors. The small intestines were first dissected into the duodenum, jejunum, and ileum, followed by collagenase digestion of the intestinal lumen. The isolated mucosa was gently homogenized to yield multiple cellular fragments, which were then cryopreserved in a programmable liquid cell freezer and stored in liquid nitrogen. After thawing and recovery, CHIM retained robust cytochrome P450 (P450) and non-P450 drug-metabolizing enzyme activities and demonstrated dose-dependent induction of transcription of CYP24A1 (approximately 300-fold) and CYP3A4 (approximately 3-fold) by vitamin D₃ as well as induction of CYP3A4 (approximately 3-fold) by rifampin after 24 hours of treatment. Dose-dependent decreases in cell viability quantified by cellular ATP content were observed for naproxen and acetaminophen, with higher enterotoxicity observed for naproxen, consistent with that observed in humans in vivo. These results suggest that CHIM may be a useful in vitro experimental model for the evaluation of enteric drug properties, including drug metabolism, drug-drug interactions, and drug toxicity.