Physiologically based pharmacokinetic model of brivaracetam to predict the exposure and dose exploration in hepatic impairment and elderly populations.

Brivaracetam (BRV) is a new third-generation antiseizure medication for the treatment of focal epileptic seizures. Its use has been increasing among epileptic populations in recent years, but pharmacokinetic (PK) behavior may change in hepatic impairment and the elderly populations. Due to ethical constraints, clinical trials are difficult to conduct and data are limited. This study used PK-Sim® to develop a physiologically based pharmacokinetic (PBPK) model for adults and extrapolate it to hepatic impairment and the elderly populations. The model was evaluated with clinical PK data, and dosage explorations were conducted. For the adult population with mild hepatic impairment, the dose is recommended to be adjusted to 70% of the recommended dose, and to 60% for moderate and severe hepatic impairment. For the elderly population with mild hepatic impairment under 80 years old, it is recommended that the dose be adjusted to 60% of the recommended dose and to 50% for moderate and severe conditions. The elderly population with hepatic impairment over 80 years old is adjusted to 50% of the recommended dose for all stages. Healthy elderly do not need to adjust. The BRV PBPK model was successfully developed, studying exposure in hepatic impairment and elderly populations and optimizing dosing regimens.

PBPK model-based gender-specific risk assessment of N-nitrosodimethylamine (NDMA) using human biomonitoring data.

Despite several screening levels for NDMA reported in water, soil, air, and drugs, the human risk assessment using biomonitoring concentrations has not been performed. In this study, gender-specific exposure guidance values were determined in humans, then biomonitoring measurements in healthy Korean subjects (32 men and 40 women) were compared to the exposure guidance values to evaluate the current exposure level to NDMA. For the human risk assessment of NDMA, the gender-specific physiologically based pharmacokinetic (PBPK) model was developed in humans using proper physiological parameters, partition coefficients, and biochemical parameters. Using the PBPK model, a Monte Carlo simulation was performed to describe the magnitudes of inter-individual variability and uncertainty on the single model predictions. The PBPK modeling and Monte Carlo simulation allowed the estimation of the relationship between external dose and blood concentration for the risk assessment. The procedure for the human risk assessment was summarized as follows: (1) estimating a steady-state blood concentration (Cavg) corresponding to the daily no observed adverse effect level (NOAEL) administration in rats; (2) applying uncertainty factors (UFs) for deriving the human Cavg; (3) determining the exposure guidance values as screening criteria; (4) interpreting the human biomonitoring measurements by forward and reverse dosimetry approaches. Using the biomonitoring concentrations, current daily exposures to NDMA were estimated to be 3.95 µg/day/kg for men and 10.60 µg/day/kg for women, respectively. The result of the study could be used as a basis for implementing further risk management and regulatory decision-making for NDMA.

Styrene lung cancer risk assessment: an alternative evaluation of human lung cancer risk assuming mouse lung tumors are potentially human relevant and operating by a threshold-based non-genotoxic mode of action.

Rodent inhalation studies indicate styrene is a mouse lung-specific carcinogen. Mode-of-action (MOA) analyses indicate that the lung tumors cannot be excluded as weakly quantitatively relevant to humans due to shared oxidative metabolites detected in rodents and humans. However, styrene also is not genotoxic following in vivo dosing. The objective of this review was to characterize occupational and general population cancer risks by conservatively assuming mouse lung tumors were relevant to humans but operating by a non-genotoxic MOA. Inhalation cancer values reference concentrations for respective occupational and general population exposures (RfCcar-occup and RfCcar-genpop) were derived from initial benchmark dose (BMD) modeling of mouse inhalation tumor dose-response data. An overall lowest BMDL10 of 4.7 ppm was modeled for lung tumors, which was further duration- and dose-adjusted by physiologically based pharmacokinetic (PBPK) modeling to derive RfCcar-occup/genpop values of 6.2 ppm and 0.8 ppm, respectively. With the exception of open-mold fiber reinforced composite workers not using personal protective equipment (PPE), the RfCcar-occup/genpop values are greater than typical occupational and general population human exposures, thus indicating styrene exposures represent a low potential for human lung cancer risk. Consistent with this conclusion, a review of styrene occupational epidemiology did not support a conclusion of an association between styrene exposure and lung cancer occurrence, and further supports a conclusion that the conservatively derived RfCcar-occup is lung cancer protective.

Beyond the basics: A deep dive into parameter estimation for advanced PBPK and QSP models.

Physiologically-based pharmacokinetic (PBPK) models and quantitative systems pharmacology (QSP) models have contributed to drug development strategies. The parameters of these models are commonly estimated by capturing observed values using the nonlinear least-squares method. Software packages for PBPK and QSP modeling provide a range of parameter estimation algorithms. To choose the most appropriate method, modelers need to understand the basic concept of each approach. This review provides a general introduction to the key points of parameter estimation with a focus on the PBPK and QSP models, and the respective parameter estimation algorithms. The latter part assesses the performance of five parameter estimation algorithms - the quasi-Newton method, Nelder-Mead method, genetic algorithm, particle swarm optimization, and Cluster Gauss-Newton method - using three examples of PBPK and QSP modeling. The assessment revealed that some parameter estimation results were significantly influenced by the initial values. Moreover, the choice of algorithms demonstrating good estimation results heavily depends on factors such as model structure and the parameters to be estimated. To obtain credible parameter estimation results, it is advisable to conduct multiple rounds of parameter estimation under different conditions, employing various estimation algorithms.

Preface to the special issue of Food and Chemical Toxicology on "New approach methodologies and machine learning in food safety and chemical risk assessment: Development of reproducible, open-source, and user-friendly tools for exposure, toxicokinetic, and toxicity assessments in the 21st century".

This Special Issue contains articles on applications of various new approach methodologies (NAMs) in the field of toxicology and risk assessment. These NAMs include in vitro high-throughput screening, quantitative structure-activity relationship (QSAR) modeling, physiologically based pharmacokinetic (PBPK) modeling, network toxicology analysis, molecular docking simulation, omics, machine learning, deep learning, and "template-and-anchor" multiscale computational modeling. These in vitro and in silico approaches complement each other and can be integrated together to support different applications of toxicology, including food safety assessment, dietary exposure assessment, chemical toxicity potency screening and ranking, chemical toxicity prediction, chemical toxicokinetic simulation, and to investigate the potential mechanisms of toxicities, as introduced further in selected articles in this Special Issue.

Applying Physiologically Based Pharmacokinetic Modeling to Interpret Carbamazepine's Nonlinear Pharmacokinetics and Its Induction Potential on Cytochrome P450 3A4 and Cytochrome P450 2C9 Enzymes.

Carbamazepine (CBZ) is commonly prescribed for epilepsy and frequently used in polypharmacy. However, concerns arise regarding its ability to induce the metabolism of other drugs, including itself, potentially leading to the undertreatment of co-administered drugs. Additionally, CBZ exhibits nonlinear pharmacokinetics (PK), but the root causes have not been fully studied. This study aims to investigate the mechanisms behind CBZ's nonlinear PK and its induction potential on CYP3A4 and CYP2C9 enzymes. To achieve this, we developed and validated a physiologically based pharmacokinetic (PBPK) parent-metabolite model of CBZ and its active metabolite Carbamazepine-10,11-epoxide in GastroPlus®. The model was utilized for Drug-Drug Interaction (DDI) prediction with CYP3A4 and CYP2C9 victim drugs and to further explore the underlying mechanisms behind CBZ's nonlinear PK. The model accurately recapitulated CBZ plasma PK. Good DDI performance was demonstrated by the prediction of CBZ DDIs with quinidine, dolutegravir, phenytoin, and tolbutamide; however, with midazolam, the predicted/observed DDI AUClast ratio was 0.49 (slightly outside of the two-fold range). CBZ's nonlinear PK can be attributed to its nonlinear metabolism caused by autoinduction, as well as nonlinear absorption due to poor solubility. In further applications, the model can help understand DDI potential when CBZ serves as a CYP3A4 and CYP2C9 inducer.

Translational two-pore PBPK model to characterize whole-body disposition of different-size endogenous and exogenous proteins.

Two-pore physiologically based pharmacokinetic (PBPK) modeling has demonstrated its potential in describing the pharmacokinetics (PK) of different-size proteins. However, all existing two-pore models lack either diverse proteins for validation or interspecies extrapolation. To fill the gap, here we have developed and optimized a translational two-pore PBPK model that can characterize plasma and tissue disposition of different-size proteins in mice, rats, monkeys, and humans. Datasets used for model development include more than 15 types of proteins: IgG (150 kDa), F(ab)2 (100 kDa), minibody (80 kDa), Fc-containing proteins (205, 200, 110, 105, 92, 84, 81, 65, or 60 kDa), albumin conjugate (85.7 kDa), albumin (67 kDa), Fab (50 kDa), diabody (50 kDa), scFv (27 kDa), dAb2 (23.5 kDa), proteins with an albumin-binding domain (26, 23.5, 22, 16, 14, or 13 kDa), nanobody (13 kDa), and other proteins (110, 65, or 60 kDa). The PBPK model incorporates: (i) molecular weight (MW)-dependent extravasation through large and small pores via diffusion and filtration, (ii) MW-dependent renal filtration, (iii) endosomal FcRn-mediated protection from catabolism for IgG and albumin-related modalities, and (iv) competition for FcRn binding from endogenous IgG and albumin. The finalized model can well characterize PK of most of these proteins, with area under the curve predicted within two-fold error. The model also provides insights into contribution of renal filtration and lysosomal degradation towards total elimination of proteins, and contribution of paracellular convection/diffusion and transcytosis towards extravasation. The PBPK model presented here represents a cross-modality, cross-species platform that can be used for development of novel biologics.

PBPK modeling to predict the pharmacokinetics of venlafaxine and its active metabolite in different CYP2D6 genotypes and drug-drug interactions with clarithromycin and paroxetine.

Venlafaxine, a serotonin-norepinephrine reuptake inhibitor (SNRI), is indicated for the treatment of major depressive disorder, social anxiety disorder, generalized anxiety disorder, and panic disorder. Venlafaxine is metabolized to the active metabolite desvenlafaxine mainly by CYP2D6. Genetic polymorphism of CYP2D6 and coadministration with other medications can significantly affect the pharmacokinetics and/or pharmacodynamics of venlafaxine and its active metabolite. This study aimed to establish the PBPK models of venlafaxine and its active metabolite related to CYP2D6 genetic polymorphism and to predict drug-drug interactions (DDIs) with clarithromycin and paroxetine in different CYP2D6 genotypes. Clinical pharmacogenomic data for venlafaxine and desvenlafaxine were collected to build the PBPK model. Physicochemical and absorption, distribution, metabolism, and excretion (ADME) characteristics of respective compounds were obtained from previously reported data, predicted by the PK-Sim® software, or optimized to capture the plasma concentration-time profiles. Model evaluation was performed by comparing the predicted pharmacokinetic parameters and plasma concentration-time profiles to the observed data. Predicted plasma concentration-time profiles of venlafaxine and its active metabolite were visually similar to the observed profiles and all predicted AUC and Cmax values for respective compounds were included in the twofold error range of observed values in non-genotyped populations and different CYP2D6 genotypes. When clarithromycin or clarithromycin plus paroxetine was concomitantly administered, predicted plasma concentration-time profiles of venlafaxine properly captured the observed profiles in two different CYP2D6 genotypes and all predicted DDI ratios for AUC and Cmax were included within the acceptance range. Consequently, the present model successfully captured the pharmacokinetic alterations of venlafaxine and its active metabolite according to CYP2D6 genetic polymorphism as well as the DDIs between venlafaxine and two CYP inhibitors. The present model can be used to predict the pharmacokinetics of venlafaxine and its active metabolite considering different races, ages, coadministered drugs, and CYP2D6 activity of individuals and it can contribute to individualized pharmacotherapy of venlafaxine.

Use of quantitative in vitro to in vivo extrapolation (QIVIVE) for the assessment of non-combustible next-generation product aerosols.

With the use of in vitro new approach methodologies (NAMs) for the assessment of non-combustible next-generation nicotine delivery products, new extrapolation methods will also be required to interpret and contextualize the physiological relevance of these results. Quantitative in vitro to in vivo extrapolation (QIVIVE) can translate in vitro concentrations into in-life exposures with physiologically-based pharmacokinetic (PBPK) modelling and provide estimates of the likelihood of harmful effects from expected exposures. A major challenge for evaluating inhalation toxicology is an accurate assessment of the delivered dose to the surface of the cells and the internalized dose. To estimate this, we ran the multiple-path particle dosimetry (MPPD) model to characterize particle deposition in the respiratory tract and developed a PBPK model for nicotine that was validated with human clinical trial data for cigarettes. Finally, we estimated a Human Equivalent Concentration (HEC) and predicted plasma concentrations based on the minimum effective concentration (MEC) derived after acute exposure of BEAS-2B cells to cigarette smoke (1R6F), or heated tobacco product (HTP) aerosol at the air liquid interface (ALI). The MPPD-PBPK model predicted the in vivo data from clinical studies within a factor of two, indicating good agreement as noted by WHO International Programme on Chemical Safety (2010) guidance. We then used QIVIVE to derive the exposure concentration (HEC) that matched the estimated in vitro deposition point of departure (POD) (MEC cigarette = 0.38 puffs or 11.6 µg nicotine, HTP = 22.9 puffs or 125.6 µg nicotine) and subsequently derived the equivalent human plasma concentrations. Results indicate that for the 1R6F cigarette, inhaling 1/6th of a stick would be required to induce the same effects observed in vitro, in vivo. Whereas, for HTP it would be necessary to consume 3 sticks simultaneously to induce in vivo the effects observed in vitro. This data further demonstrates the reduced physiological potency potential of HTP aerosol compared to cigarette smoke. The QIVIVE approach demonstrates great promise in assisting human health risk assessments, however, further optimization and standardization are required for the substantiation of a meaningful contribution to tobacco harm reduction by alternative nicotine delivery products.

Exploring the potential and challenges of developing physiologically-based toxicokinetic models to support human health risk assessment of microplastic and nanoplastic particles.

Microplastics (MPs) and nanoplastics (NPs) pollution has emerged as a significant and widespread environmental issue. Humans are inevitably exposed to MPs and NPs via ingestion, inhalation, and dermal contacts from various sources. However, mechanistic knowledge of their distribution, interaction, and potency in the body is still lacking. To address this knowledge gap, we have undertaken the task of elucidating the toxicokinetic (TK) behaviors of MPs and NPs, aiming to provide mechanistic information for constructing a conceptual physiologically based toxicokinetic (PBTK) model to support in silico modeling approaches. Our effort involved a thorough examination of the existing literature and data collation on the presence of MPs in the human body and in vitro/ex vivo/in vivo biodistribution across various cells and tissues. By comprehending the absorption, distribution, metabolism, and excretion mechanisms of MPs and NPs in relation to their physicochemical attributes, we established a foundational understanding of the link between external exposure and internal tissue dosimetry. We observed that particle size and surface chemistry have been thoroughly explored in previous experimental studies. However, certain attributes, such as polymer type, shape, and biofilm/biocorona, warrant attention and further examination. We discussed the fundamental disparities in TK properties of MPs/NPs from those of engineered nanoparticles. We proposed a preliminary PBTK framework with several possible modeling approaches and discussed existing challenges for further investigation. Overall, this article provides a comprehensive compilation of existing TK data of MPs/NPs, a critical overview of TK processes and mechanisms, and proposes potential PBTK modeling approaches, particularly regarding their applicability to the human system, and outlines future perspectives for developing PBTK models and their integration into human health risk assessment of MPs and NPs.

Physiologically based pharmacokinetic modeling for confirming the role of CYP3A1/2 and P-glycoprotein in detoxification mechanism between glycyrrhizic acid and aconitine in rats.

Fuzi, an effective common herb, is often combined with Gancao to treat disease in clinical practice with enhancing its efficacy and alleviating its toxicity. The major toxic and bioactive compounds in Fuzi and Gancao are aconitine (AC) and glycyrrhizic acid (GL), respectively. This study aims to elucidate detoxification mechanism between AC and GL from pharmacokinetic perspective using physiologically based pharmacokinetic (PBPK) model. In vitro experiments exhibited that AC was mainly metabolized by CYP3A1/2 in rat liver microsomes and transported by P-glycoprotein (P-gp) in Caco-2 cells. Kinetics assays showed that the Km and Vmax of AC towards CYP3A1/2 were 2.38 µM and 57.3 pmol/min/mg, respectively, whereas that of AC towards P-gp was 11.26 µM and 147.1 pmol/min/mg, respectively. GL markedly induced the mRNA expressions of CYP3A1/2 and MDR1a/b in rat primary hepatocytes. In vivo studies suggested that the intragastric and intravenous administration of GL significantly reduced systemic exposure of AC by 27% and 33%, respectively. Drug-drug interaction (DDI) model of PBPK predicted that co-administration of GL would decrease the exposure of AC by 39% and 45% in intragastric and intravenous dosing group, respectively. The consistency between predicted data and observed data confirmed that the upregulation of CYP3A1/2 and P-gp was the crucial detoxification mechanism between AC and GL. Thus, this study provides a demonstration for elucidating the compatibility mechanisms of herbal formula using PBPK modeling and gives support for the clinical co-medication of Fuzi and Gancao.

Docetaxel, cyclophosphamide, and epirubicin: application of PBPK modeling to gain new insights for drug-drug interactions.

The new adjuvant chemotherapy of docetaxel, epirubicin, and cyclophosphamide has been recommended for treating breast cancer. It is necessary to investigate the potential drug-drug Interactions (DDIs) since they have a narrow therapeutic window in which slight differences in exposure might result in significant differences in treatment efficacy and tolerability. To guide clinical rational drug use, this study aimed to evaluate the DDI potentials of docetaxel, cyclophosphamide, and epirubicin in cancer patients using physiologically based pharmacokinetic (PBPK) models. The GastroPlus™ was used to develop the PBPK models, which were refined and validated with observed data. The established PBPK models accurately described the pharmacokinetics (PKs) of three drugs in cancer patients, and the predicted-to-observed ratios of all the PK parameters met the acceptance criterion. The PBPK model predicted no significant changes in plasma concentrations of these drugs during co-administration, which was consistent with the observed clinical phenomenon. Besides, the verified PBPK models were then used to predict the effect of other Cytochrome P450 3A4 (CYP3A4) inhibitors/inducers on these drug exposures. In the DDI simulation, strong CYP3A4 modulators changed the exposure of three drugs by 0.71-1.61 fold. Therefore, patients receiving these drugs in combination with strong CYP3A4 inhibitors should be monitored regularly to prevent adverse reactions. Furthermore, co-administration of docetaxel, cyclophosphamide, or epirubicin with strong CYP3A4 inducers should be avoided. In conclusion, the PBPK models can be used to further investigate the DDI potential of each drug and to develop dosage recommendations for concurrent usage by additional perpetrators or victims.

Characterization of Pediatric Rectal Absorption, Drug Disposition, and Sedation Level for Midazolam Gel Using Physiologically Based Pharmacokinetic/Pharmacodynamic Modeling.

This study aims to explore and characterize the role of pediatric sedation via rectal route. A pediatric physiologically based pharmacokinetic-pharmacodynamic (PBPK/PD) model of midazolam gel was built and validated to support dose selection for pediatric clinical trials. Before developing the rectal PBPK model, an intravenous PBPK model was developed to determine drug disposition, specifically by describing the ontogeny model of the metabolic enzyme. Pediatric rectal absorption was developed based on the rectal PBPK model of adults. The improved Weibull function with permeability, surface area, and fluid volume parameters was used to extrapolate pediatric rectal absorption. A logistic regression model was used to characterize the relationship between the free concentrations of midazolam and the probability of sedation. All models successfully described the PK profiles with absolute average fold error (AAFE) < 2, especially our intravenous PBPK model that extended the predicted age to preterm. The simulation results of the PD model showed that when the free concentrations of midazolam ranged from 3.9 to 18.4 ng/mL, the probability of "Sedation" was greater than that of "Not-sedation" states. Combined with the rectal PBPK model, the recommended sedation doses were in the ranges of 0.44-2.08 mg/kg for children aged 2-3 years, 0.35-1.65 mg/kg for children aged 4-7 years, 0.24-1.27 mg/kg for children aged 8-12 years, and 0.20-1.10 mg/kg for adolescents aged 13-18 years. Overall, this model mechanistically quantified drug disposition and effect of midazolam gel in the pediatric population, accurately predicted the observed clinical data, and simulated the drug exposure for sedation that will inform dose selection for following pediatric clinical trials.

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.

Facing the Facts of Altered Plasma Protein Binding: Do Current Models Correctly Predict Changes in Fraction Unbound in Special Populations?

Accounting for variability in plasma protein binding of drugs is an essential input to physiologically-based pharmacokinetic (PBPK) models of special populations. Prediction of fraction unbound in plasma (fu) in such populations typically considers changes in plasma protein concentration while assuming that the binding affinity remains unchanged. A good correlation between predicted vs observed fu data reported for various drugs in a given special population is often used as a justification for such predictive methods. However, none of these analyses evaluated the prediction of the fold-change in fu in special populations relative to the reference population. This would be a more appropriate assessment of the predictivity, analogous to drug-drug interactions. In this study, predictive performance of the single protein binding model was assessed by predicting fu for alpha-1-acid glycoprotein and albumin bound drugs in hepatic impairment, renal impairment, paediatric, elderly, patients with inflammatory disease, and in different ethnic groups for a dataset of >200 drugs. For albumin models, the concordance correlation coefficients for predicted fu were >0.90 for 16 out of 17 populations with sub-groups, indicating strong agreement between predicted and observed values. In contrast, concordance correlation coefficients for predicted fold-change in fu for the same dataset were <0.38 for all populations and sub-groups. Trends were similar for alpha-1-acid glycoprotein models. Accordingly, the predictions of fu solely based on changes in protein concentrations in plasma cannot explain the observed values in some special populations. We recommend further consideration of the impact of changes in special populations to endogenous substances that competitively bind to plasma proteins, and changes in albumin structure due to posttranslational modifications. PBPK models of special populations for highly bound drugs should preferably use measured fu data to ensure reliable prediction of drug exposure or compare predicted unbound drug exposure between populations knowing that these will not be sensitive to changes in fu.

Prediction of ROS1 and TRKA/B/C occupancy in plasma and cerebrospinal fluid for entrectinib alone and in DDIs using physiologically based pharmacokinetic (PBPK) modeling approach.

PURPOSE: Entrectinib (ENT) is a potent c-ros oncogene 1(ROS1) and neurotrophic tyrosine receptor kinase (NTRKA/B/C) inhibitor. To determine the optimum dosage of ENT using ROS1 and NTRKA/B/C occupancy in plasma and cerebrospinal fluid (CSF) in drug-drug interactions (DDIs), physiologically-based pharmacokinetic (PBPK) models for healthy subjects and cancer population were developed for ENT and M5 (active metabolite). METHODS: The PBPK models were built using the modeling parameters of ENT and M5 that were mainly derived from the published paper on the ENT PBPK model, and then validated by the observed pharmacokinetics (PK) in plasma and CSF from healthy subjects and patients. RESULTS: The PBPK model showed that AUC, Cmax, and Ctrough ratios between predictions and observations are within the range of 0.5-2.0, except that the M5 AUC ratio is slightly above 2.0 (2.34). Based on the efficacy (> 75% occupancy for ROS1 and NTRKA/B/C) and safety (AUC < 160 µM·h and Cmax < 8.9 µM), the appropriate dosing regimens were identified. The appropriate dosage is 600 mg once daily (OD) when administered alone, reduced to 200 mg and 400 mg OD with itraconazole and fluconazole, respectively. ENT is not recommended for co-administration with rifampicin or efavirenz, but is permitted with fluvoxamine or dexamethasone. CONCLUSION: The PBPK models can serve as a powerful approach to predict ENT concentration as well as ROS1 and NTRKA/B/C occupancy in plasma and CSF.

Whole-Body Disposition and Physiologically Based Pharmacokinetic Modeling of Adeno-Associated Viruses and the Transgene Product.

To facilitate model-informed drug development (MIDD) of adeno-associated virus (AAV) therapy, here we have developed a physiologically based pharmacokinetic (PBPK) model for AAVs following preclinical investigation in mice. After 2E11 Vg/mouse dose of AAV8 and AAV9 encoding a monoclonal antibody (mAb) gene, whole-body disposition of both the vector and the transgene mAb was evaluated over 3 weeks. At steady-state, the following tissue-to-blood (T/B) concentration ratios were found for AAV8/9: ∼50 for liver; ∼10 for heart and muscle; ∼2 for brain, lung, kidney, adipose, and spleen; ≤1 for bone, skin, and pancreas. T/B values for mAb were compared with the antibody biodistribution coefficients, and five different clusters of organs were identified based on their transgene expression profile. All the biodistribution data were used to develop a novel AAV PBPK model that incorporates: (i) whole-body distribution of the vector; (ii) binding, internalization, and intracellular processing of the vector; (iii) transgene expression and secretion; and (iv) whole-body disposition of the secreted transgene product. The model was able to capture systemic and tissue PK of the vector and the transgene-produced mAb reasonably well. Pathway analysis of the PBPK model suggested that liver, muscle, and heart are the main contributors for the secreted transgene mAb. Unprecedented PK data and the novel PBPK model developed here provide the foundation for quantitative systems pharmacology (QSP) investigations of AAV-mediated gene therapies. The PBPK model can also serve as a quantitative tool for preclinical study design and preclinical-to-clinical translation of AAV-based gene therapies.

Physiologically based pharmacokinetic (PBPK) modeling of pitavastatin in relation to SLCO1B1 genetic polymorphism.

Pitavastatin, a potent 3-hydroxymethylglutaryl coenzyme A reductase inhibitor, is indicated for the treatment of hypercholesterolemia and mixed dyslipidemia. Hepatic uptake of pitavastatin is predominantly occupied by the organic anion transporting polypeptide 1B1 (OATP1B1) and solute carrier organic anion transporter family member 1B1 (SLCO1B1) gene, which is a polymorphic gene that encodes OATP1B1. SLCO1B1 genetic polymorphism significantly alters the pharmacokinetics of pitavastatin. This study aimed to establish the physiologically based pharmacokinetic (PBPK) model to predict pitavastatin pharmacokinetics according to SLCO1B1 genetic polymorphism. PK-Sim® version 10.0 was used to establish the whole-body PBPK model of pitavastatin. Our pharmacogenomic data and a total of 27 clinical pharmacokinetic data with different dose administration and demographic properties were used to develop and validate the model, respectively. Physicochemical properties and disposition characteristics of pitavastatin were acquired from previously reported data or optimized to capture the plasma concentration-time profiles in different SLCO1B1 diplotypes. Model evaluation was performed by comparing the predicted pharmacokinetic parameters and profiles to the observed data. Predicted plasma concentration-time profiles were visually similar to the observed profiles in the non-genotyped populations and different SLCO1B1 diplotypes. All fold error values for AUC and Cmax were included in the two fold range of observed values. Thus, the PBPK model of pitavastatin in different SLCO1B1 diplotypes was properly established. The present study can be useful to individualize the dose administration strategy of pitavastatin in individuals with various ages, races, and SLCO1B1 diplotypes.

PBPK modeling to predict the pharmacokinetics of pantoprazole in different CYP2C19 genotypes.

Pantoprazole is used to treat gastroesophageal reflux disease (GERD), maintain healing of erosive esophagitis (EE), and control symptoms related to Zollinger-Ellison syndrome (ZES). Pantoprazole is mainly metabolized by cytochrome P450 (CYP) 2C19, converting to 4'-demethyl pantoprazole. CYP2C19 is a genetically polymorphic enzyme, and the genetic polymorphism affects the pharmacokinetics and/or pharmacodynamics of pantoprazole. In this study, we aimed to establish the physiologically based pharmacokinetic (PBPK) model to predict the pharmacokinetics of pantoprazole in populations with various CYP2C19 metabolic activities. A comprehensive investigation of previous reports and drug databases was conducted to collect the clinical pharmacogenomic data, physicochemical data, and disposition properties of pantoprazole, and the collected data were used for model establishment. The model was evaluated by comparing the predicted plasma concentration-time profiles and/or pharmacokinetic parameters (AUC and Cmax) with the clinical observation results. The predicted plasma concentration-time profiles in different CYP2C19 phenotypes properly captured the observed profiles. All fold error values for AUC and Cmax were included in the two-fold range. Consequently, the minimal PBPK model for pantoprazole related to CYP2C19 genetic polymorphism was properly established and it can predict the pharmacokinetics of pantoprazole in different CYP2C19 phenotypes. The present model can broaden the insight into the individualized pharmacotherapy for pantoprazole.

Physiologically-based pharmacokinetic model for evaluating gender-specific exposures of N-nitrosodimethylamine (NDMA).

N-nitrosodimethylamine (NDMA) is classified as a human carcinogen and could be produced by both natural and industrial processes. Although its toxicity and histopathology have been well-studied in animal species, there is insufficient data on the blood and tissue exposures that can be correlated with the toxicity of NDMA. The purpose of this study was to evaluate gender-specific pharmacokinetics/toxicokinetics (PKs/TKs), tissue distribution, and excretion after the oral administration of three different doses of NDMA in rats using a physiologically-based pharmacokinetic (PBPK) model. The major target tissues for developing the PBPK model and evaluating dose metrics of NDMA included blood, gastrointestinal (GI) tract, liver, kidney, lung, heart, and brain. The predictive performance of the model was validated using sensitivity analysis, (average) fold error, and visual inspection of observations versus predictions. Then, a Monte Carlo simulation was performed to describe the magnitudes of inter-individual variability and uncertainty of the single model predictions. The developed PBPK model was applied for the exposure simulation of daily oral NDMA to estimate blood concentration ranges affecting health effects following acute-duration (≤ 14 days), intermediate-duration (15-364 days), and chronic-duration (≥ 365 days) intakes. The results of the study could be used as a scientific basis for interpreting the correlation between in vivo exposures and toxicological effects of NDMA.