Real-world application of PBPK in drug discovery

The utility of PBPK models in support of drug development has been well documented. During the discovery stage, PBPK has increasingly been applied for early risk assessment, prediction of human dose, toxicokinetic dose projection and early formulation assessment. Previous review articles have proposed model building and application strategies for PBPK-based first in human predictions with comprehensive descriptions of the individual components of PBPK models. This includes the generation of decision trees, based on comprehensive literature reviews, to guide the application of PBPK in the discovery setting. The goal of this mini review is to provide additional guidance on the real-world application of PBPK, in support of the discovery stage of drug development. In this mini review, our goal is to provide guidance on the typical steps involved in the development and application of a PBPK model during drug discovery to assist in decision making. We have illustrated our recommended approach through description of case examples, where PBPK has been successfully applied to aid in human PK projection, candidate selection and prediction of drug interaction liability for parent and metabolite. Through these case studies, we have highlighted fundamental issues, including pre-verification in preclinical species, the application of empirical scalars in the prediction of in vivo clearance from in vitro systems, in silico prediction of permeability and the exploration of aqueous and biorelevant solubility data to predict dissolution. In addition, current knowledge gaps have been highlighted and future directions proposed. Significance Statement Through description of three case studies, we have highlighted the fundamental principles of PBPK application during drug discovery. These include pre-verification of the model in preclinical species, application of empirical scalars where necessary in the prediction of clearance, in silico prediction of permeability, and the exploration of aqueous and biorelevant solubility data to predict dissolution. In addition, current knowledge gaps have been highlighted and future directions proposed.

Physiologically-Based Pharmacokinetic Modeling of PAXLOVID™ with First-Order Absorption Kinetics

PURPOSE

PAXLOVID™ is nirmatrelvir tablets co-packaged with ritonavir tablets. Ritonavir is used as a pharmacokinetics (PK) enhancer to reduce metabolism and increase exposure of nirmatrelvir. This is the first disclosure of Paxlovid physiologically-based pharmacokinetic (PBPK) model.

METHODS

Nirmatrelvir PBPK model with first-order absorption kinetics was developed using in vitro, preclinical, and clinical data of nirmatrelvir in the presence and absence of ritonavir. Clearance and volume of distribution were derived from nirmatrelvir PK obtained using a spray-dried dispersion (SDD) formulation where it is considered to be dosed as an oral solution, and absorption is near complete. The fraction of nirmatrelvir metabolized by CYP3A was estimated based on in vitro and clinical ritonavir drug-drug interaction (DDI) data. First-order absorption parameters were established for both SDD and tablet formulation using clinical data. Nirmatrelvir PBPK model was verified with both single and multiple dose human PK data, as well as DDI studies. Simcyp® first-order ritonavir compound file was also verified with additional clinical data.

RESULTS

The nirmatrelvir PBPK model described the observed PK profiles of nirmatrelvir well with predicted AUC and Cmax values within ± 20% of the observed. The ritonavir model performed well resulting in predicted values within twofold of observed.

CONCLUSIONS

Paxlovid PBPK model developed in this study can be applied to predict PK changes in special populations, as well as model the effect of victim and perpetrator DDI. PBPK modeling continues to play a critical role in accelerating drug discovery and development of potential treatments for devastating diseases such as COVID-19. NCT05263895, NCT05129475, NCT05032950 and NCT05064800.

Development and validation of two bioanalysis methods for the determination of etimicin in human serum and urine by liquid chromatography-tandem mass spectrometry: Applications to a human pharmacokinetic and breakpoint study

Etimicin is a fourth-generation aminoglycoside antibiotic. It has potent activity and low toxicity when employed for the treatment of Gram-negative and Gram-positive bacterial infections. The pharmacokinetics of etimicin in humans have not been elucidated completely. Two liquid chromatography-tandem mass spectrometry (LC-MS/MS) bioanalytical methods, without the use of any ion-pairing reagents, were developed and validated for the quantification of etimicin in human samples of serum and urine. Using a deuterated reagent as the internal standard, analytes in serum and urine samples were extracted by protein precipitation and dilution before LC-MS/MS analysis, respectively. For the two methods, chromatographic separations were undertaken under isocratic elution of water-ammonia solution-acetic acid (96:3.6:0.2, v/v/v) and methanol at 50%:50% and a flow rate of 0.35 ml/min within 5 min. A Waters XTerra MS C18 column (2.1 × 150 mm, 3.5 μm) and a column temperature of 40°C were chosen. A Sciex Qtrap 5500 mass spectrometer equipped with an electrospray ion source was used in both methods under multiple-reaction monitoring in positive-ion mode. The two methods showed good linearity, accuracy, and precision with high recovery and a minimal matrix effect in the range of 50.0-20000 ng/ml for serum samples and 50.0-10000 ng/ml for urine samples, respectively. Carry-over effects were not observed. Etimicin remained stable in human samples of serum or urine under the storage, preparation, and analytical conditions of the two methods. These two simple and reliable methods were applied successfully to a dose-escalation, phase I clinical trial of etimicin in Chinese healthy volunteers after intravenous administration of single and multiple doses. Based on these two methods we ascertained, for the first time, the comprehensive pharmacokinetics of etimicin in humans, which will be used for the exploration of the breakpoint research further.

CYP2C8-Mediated Formation of a Human Disproportionate Metabolite of the Selective NaV1.7 Inhibitor DS-1971a, a Mixed Cytochrome P450 and Aldehyde Oxidase Substrate

Predicting human disproportionate metabolites is difficult, especially when drugs undergo species-specific metabolism mediated by cytochrome P450s (P450s) and/or non-P450 enzymes. This study assessed human metabolites of DS-1971a, a potent Na_v1.7-selective blocker, by performing human mass balance studies and characterizing DS-1971a metabolites, in accordance with the Metabolites in Safety Testing (MIST) guidance. In addition, we investigated the mechanism by which the major human disproportionate metabolite (M1) was formed. After oral administration of radiolabeled DS-1971a, the major metabolites in human plasma were P450-mediated monoxidized metabolites M1 and M2 with area under the curve ratios of 27% and 10% of total drug-related exposure, respectively; the minor metabolites were dioxidized metabolites produced by aldehyde oxidase and P450s. By comparing exposure levels of M1 and M2 between humans and safety assessment animals, M1 but not M2 was found to be a human disproportionate metabolite, requiring further characterization under the MIST guidance. Incubation studies with human liver microsomes indicated that CYP2C8 was responsible for the formation of M1. Docking simulation indicated that, in the formation of M1 and M2, there would be hydrogen bonding and/or electrostatic interactions between the pyrimidine and sulfonamide moieties of DS-1971a and amino acid residues Ser100, Ile102, Ile106, Thr107, and Asn217 in CYP2C8, and that the cyclohexane ring of DS-1971a would be located near the heme iron of CYP2C8. These results clearly indicate that M1 is the predominant metabolite in humans and a human disproportionate metabolite due to species-specific differences in metabolism. Significance Statement This report is the first to show a human disproportionate metabolite generated by CYP2C8-mediated primary metabolism. We clearly demonstrate that DS-1971a, a mixed aldehyde oxidase and cytochrome P450 substrate, was predominantly metabolized by CYP2C8 to form M1, a human disproportionate metabolite. Species differences in the formation of M1 highlight the regio- and stereoselective metabolism by CYP2C8, and the proposed interaction between DS-1971a and CYP2C8 provides new knowledge of CYP2C8-mediated metabolism of cyclohexane-containing substrates.

Machine Learning Models for Human In Vivo Pharmacokinetic Parameters with In-House Validation

Prior to clinical development, a comprehensive pharmacokinetic characterization of a novel drug is required to understand its exposure at the site of action and elimination. Accordingly, in vitro assays and animal pharmacokinetic studies are regularly employed to predict drug exposure in humans, which is often costly and time-consuming. For this reason, the prediction of human pharmacokinetics at the point of design would be of high value for drug discovery. Therefore, we have established a comprehensive data curation protocol that enables machine learning evaluation of 12 human in vivo pharmacokinetic parameters using only chemical structure information and available doses for 1001 unique compounds. These machine learning models were thoroughly investigated and validated using both an independent hold-out test set and AstraZeneca clinical data. In addition, the availability of preclinical predictions for a subset of internal clinical candidates allowed us to compare our in silico approach with state-of-the-art pharmacokinetic predictions. Based on this evaluation, three fit-for-purpose models for AUC PO (R_test² = 0.63; RMSE_test = 0.76), C_max PO (R_test² = 0.68; RMSE_test = 0.62), and Vd_ss IV (R_test² = 0.47; RMSE_test = 0.50) were identified. Based on the findings, our machine learning models have considerable potential for practical applications in drug discovery, such as influencing decision-making in drug discovery projects and progression of drug candidates toward the clinic.

Human Variability in Carboxylesterases and carboxylesterase-related Uncertainty Factors for Chemical Risk Assessment

Carboxylesterases (CES) are an important class of enzymes involved in the hydrolysis of a range of chemicals and show large inter-individual variability in vitro. An extensive literature search was performed to identify in vivo probe substrates for CES1 and CES2 together with their protein content and enzymatic activity. Human pharmacokinetic (PK) data on Cmax, clearance, and AUC were extracted from 89 publications and Bayesian meta-analysis was performed using a hierarchical model to derive CES-related variability distributions and related uncertainty factors (UF). The CES-related variability indicated that 97.5% of healthy adults are covered by the kinetic default UF (3.16), except for clopidogrel and dabigatran etexilate. Clopidogrel is metabolised for a small amount by the polymorphic CYP2C19, which can have an impact on the overall pharmacokinetics, while the variability seen for dabigatran etexilate might be due to differences in the absorption, since this can be influenced by food intake. The overall CES-related variability was moderate to high in vivo (<CV 50%), which might be due to possible polymorphism in the enzyme but also to the small sample size available per chemical. The presented CES-related variability can be used in combination with in vitro data to derive pathway-specific distributions.

Comparative pharmacokinetic study on three formulations of Astragali Radix by an LC-MS/MS method for determination of formononetin in human plasma

Astragali Radix (AR) is a widely used traditional Chinese medicine for healing the cardiovascular, liver and immune systems. Recently, superfine pulverizing technology has been applied to developing novel formulations to improve bioavailability of the active constituents in herbs, such as ultrafine granular powder of AR. In this study, a universal and sensitive quantitative method based on LC-MS/MS was employed for determining formononetin, the main flavonoid in AR, in human plasma for comparative pharmacokinetics of three oral formulations of AR. Formononetin and IS (quercetin) were extracted by ethyl acetate from human plasma and were separated on a C₁₈ column with a mobile phase consisting of acetonitrile and 0.1% formic acid. Positive-ion electrospray-ionization mode was applied in mass spectrometric detection. The quantitative method was validated with regards to selectivity, linearity, accuracy and precision, matrix effect, extraction recovery and stability, and was applied to comparing the pharmacokinetics of ultrafine granular powder (UGP), ultrafine powder (UP) and traditional decoction pieces (TDP) of AR after oral administration. The peak concentration and areas under the concentration-time curve of formononetin in UGP and UP were significantly higher than those of TDP. UGP and UP could significantly improve the bioavailability of AR in human compared with TDP after oral administration.

Human pharmacokinetics of ginkgo terpene lactones and impact of carboxylation in blood on their platelet-activating factor antagonistic activity

Terpene lactones are a class of bioactive constituents of standardized preparations of Ginkgo biloba leaf extract, extensively used as add-on therapies in patients with ischemic cardiovascular and cerebrovascular diseases. This investigation evaluated human pharmacokinetics of ginkgo terpene lactones and impact of their carboxylation in blood. Human subjects received oral YinXing-TongZhi tablet or intravenous ShuXueNing, two standardized ginkgo preparations. Their plasma protein-binding and platelet-activating factor antagonistic activity were assessed in vitro. Their carboxylation was assessed in phosphate-buffered saline (pH 7.4) and in human plasma. After dosing YinXing-TongZhi tablet, ginkgolides A and B and bilobalide exhibited significantly higher systemic exposure levels than ginkgolides C and J; after dosing ShuXueNing, ginkgolides A, B, C, and J exhibited high exposure levels. The compounds' unbound fractions in plasma were 45-92%. Apparent oral bioavailability of ginkgolides A and B was mostly >100%, while that of ginkgolides C and J was 6-15%. Bilobalide's bioavailability was probably high but lower than that of ginkgolides A/B. Terminal half-lives of ginkgolides A, B, and C (4-7 h) after dosing ShuXueNing were shorter than their respective values (6-13 h) after dosing YinXing-TongZhi tablet. Half-life of bilobalide after dosing the tablet was around 5 h. Terpene lactones were roughly evenly distributed in various body fluids and tissues; glomerular-filtration-based renal excretion was the predominant elimination route for the ginkgolides and a major route for bilobalide. Terpene lactones circulated as trilactones and monocarboxylates. Carboxylation reduced platelet-activating factor antagonistic activity of ginkgolides A, B, and C. Ginkgolide J, bilobalide, and ginkgo flavonoids exhibited no such bioactivity. Collectively, differences in terpene lactones' exposure between the two preparations and influence of their carboxylation in blood should be considered in investigating the relative contributions of terpene lactones to ginkgo preparations' therapeutic effects. The results here will inform rational clinical use of ginkgo preparations.

A mechanistic PK/PD model for two anti-IL13 antibodies explains the difference in total IL-13 accumulation observed in clinical studies

IMA-638 and IMA-026 are humanized IgG1 monoclonal antibodies (mAbs) that target non-overlapping epitopes of IL-13. Separate first-in-human single ascending dose studies were conducted for each mAb. These studies had similar study designs, but mild to moderate asthmatics were recruited for the IMA-638 study and healthy subjects were recruited for the IMA-026 study. IMA-638 and IMA-026 showed similar pharmacokinetic (PK) profiles, but very different total IL-13 (free and drug bound IL-13) profiles; free IL13 was not measured. IMA-026 treatment induced a dose-dependent accumulation of total IL-13, while IMA-638 treatment led to a much smaller accumulation without any clear dose-response. To understand the differences between the two total IL-13 profiles and to predict the free IL-13 profiles for each mAb treatment, a mechanistic PK/pharmacodynamic model was developed. PK-related parameters were first fit to the mean PK profiles of each mAb separately; thereafter, the target-related parameters were fit to both total IL-13 profiles simultaneously. The IL-13 degradation rate was assumed to be the same for asthma patients and healthy subjects. The model suggests that an approximately 100× faster elimination of IL-13-IMA-638 complex than IL-13-IMA-026 complex could be responsible for the differences observed in total IL-13 profiles for the two mAbs. Furthermore, the model predicts that IMA-638 administration results in greater and more prolonged free IL-13 inhibition than equivalent dosing of IMA-026 despite similar binding KD and PK profile. In conclusion, joint modeling of two similar molecules provided mechanistic insight that the elimination rate of mAb-target complex can regulate the degree of free target inhibition.

Setipiprant, a selective oral antagonist of human CRTH2: relative bioavailability of a capsule and a tablet formulation in healthy female and male subjects

BACKGROUND

CRTH2 is a prostaglandin D2 receptor that plays an important role in allergic inflammation. Setipiprant is a potent CRTH2 antagonist under development for the treatment of allergic diseases.

OBJECTIVE

The aim of this study was to evaluate the tolerability and pharmacokinetics of a single oral dose of a setipiprant capsule (reference) and a tablet formulation.

METHODS

This was an open-label, 2-period, 2-way crossover, randomized study in which 20 healthy women and men (1:1 ratio) received either 2 250-mg capsules or a 500-mg tablet of setipiprant. Subjects were between 18 and 45 years old, with a body mass index of 18.0 to 28.0 kg/m(2). Differences in pharmacokinetics of setipiprant formulations were explored overall and by sex.

RESULTS

All subjects completed the study. Both formulations were well tolerated, with headache the most frequently reported adverse event (25% of subjects), followed by flatulence (15%) and somnolence and fatigue (10%). The adverse event profile in men and women and between formulations was similar. The ratios of geometric means for Cmax (0.94; 95% CI, 0.79-1.12) and AUC0-∞ (1.01; 95% CI, 0.92-1.12) were mostly within the limits of 0.80 to 1.25. When corrected for weight, the differences observed between sexes, within each treatment, for Cmax (capsules: 1.01; 95% CI, 0.71-1.44; tablet: 0.89; 95% CI, 0.62-1.26) and AUC0-∞ (capsules: 1.12; 95% CI, 0.86-1.47; tablet: 0.96; 95% CI, 0.73-1.25) were minor.

CONCLUSION

Both the setipiprant formulations were well tolerated. Setipiprant pharmacokinetics were similar between formulations, overall, and between sexes. The new tablet formulation may constitute a valid alternative to the capsule formulation for later clinical development phases. ClinicalTrials.gov identifier: NCT01877629.

Pharmacokinetics and metabolism of brotizolam in humans

Pharmacokinetic studies were performed in healthy young volunteers and in elderly patients after oral administration of single doses (0.5 mg), increasing doses (0.5-1.5 mg), and multiple doses (1.0 mg) of brotizolam. Brotizolam was absorbed quickly from the gastro-intestinal tract. Elimination half-lives were in the range of 3.6-7.9 h. In healthy young volunteers as well as in elderly patients, there was neither a tendency for brotizolam to accumulate nor was there any indication of enzyme induction. Brotizolam was metabolized almost completely into hydroxylated compounds which were conjugated prior to renal excretion. After oral administration of [14C]-brotizolam, two-thirds of excretion of radioactivity was renal and was completed within 4 days.