Biopharmaceutic In Vitro In Vivo Extrapolation (IVIV_E) Informed Physiologically-Based Pharmacokinetic Model of Ritonavir Norvir Tablet Absorption in Humans Under Fasted and Fed State Conditions

The Use of Global Sensitivity Analysis to Assess the Oral Absorption of Weakly Basic Compounds: A Case Example of Dipyridamole

OBJECTIVE

To utilize the global system analysis (GSA) in oral absorption modeling to gain a deeper understanding of system behavior, improve model accuracy, and make informed decisions during drug development.

METHODS

GSA was utilized to give insight into which drug substance (DS), drug product (DP), and/or physiological parameter would have an impact on peak plasma concentration (Cmax) and area under the curve (AUC) of dipyridamole as a model weakly basic compound. GSA guided the design of in vitro experiments and oral absorption risk assessment using FormulatedProducts v2202.1.0. The solubility and precipitation profiles of dipyridamole in different bile salt concentrations were measured. The results were then used to build a mechanistic oral absorption model.

RESULTS

GSA warranted further investigation into the precipitation kinetics and its link to the levels of bile salt concentrations. Mechanistic modeling studies demonstrated that a precipitation-integrated modeling approach appropriately predicted the mean plasma profiles, Cmax, and AUC from the clinical studies.

CONCLUSIONS

This work shows the value of GSA utilization in early development to guide in vitro experimentation and build more confidence in identifying the critical parameters for the mathematical models.

Physiological Dynamics in the Upper Gastrointestinal Tract and the Development of Gastrointestinal Absorption Models for the Immediate-Release Oral Dosage Forms in Healthy Adult Human

This review is a revisit of various oral drug absorption models developed in the past decades, focusing on how to incorporate the physiological dynamics in the upper gastrointestinal (GI) tract. For immediate-release oral drugs, GI absorption is a critical input of drug exposure and subsequent human body response, yet difficult to model largely due to the complex GI environment. One of the biggest hurdles lies at capturing the high within-subject variability (WSV) of bioavailability measures, which can be mechanistically explained by the GI physiological dynamics. A thorough summary of how GI dynamics is handled in the absorption models would promote the development of mechanism-based oral drug absorption models, aid in the design of clinical studies regarding dosing regimens and bioequivalence studies based on WSV, and advance the decision-making on formulation selection.

Dose decision of HSK7653 oral immediate release tablets in specific populations clinical trials based on mechanistic physiologically-based pharmacokinetic model

HSK7653, an oral dipeptidyl peptidase-4 inhibitor administered every 2 weeks, is a candidate for the treatment of type 2 diabetes mellitus. The major elimination pathway of HSK7653 in vivo is renal excretion, and hepatic metabolism and fecal excretion of unchanged compound contribute less to the systemic clearance of HSK7653. Considering the disposition characteristics and the potential indication population of HSK7653, evaluating the HSK7653 exposure in patients with renal impairment and geriatric populations is a prerequisite for bringing more benefits to the patients. Here, a PBPK model was developed based on in vitro experimental results, such as dissolution, permeability, and metabolism, and the in vivo renal clearance, to evaluate the effects of physiological factors and food on HSK7653 exposure in specific populations, including adult and elder individuals with renal impairment and geriatric populations. Simulation results showed that the AUC of HSK7653 increased by 46%, 82%, and 129% in adult patients with mild, moderate, and severe renal impairment, and by 56%, 78%, and 101% in patients aged 65-75, 75-85 and 85-95 years, respectively. The AUC increased in the range of 62%-83%, 98%-133%, and 153%-195% in elderly patients (65-95 years) with mild, moderate, and severe renal impairment, respectively. Moreover, two different absorption model development methods (dissolution profile method and the diffusion layer model method) predicted that food had no effect on the exposure of the same simulated population. Since the predicted AUC of HSK7653 at the 10 mg dose in various specific populations was still within the relatively flat results of the exposure-response analysis, the 10 mg dose of HSK7653 was first used to explore the exposure in the renal impairment population (CTR20221952).

Leveraging the use of in vitro and computational methods to support the development of enabling oral drug products: An InPharma commentary

Due to the strong tendency towards poorly soluble drugs in modern development pipelines, enabling drug formulations such as amorphous solid dispersions, cyclodextrins, co-crystals and lipid-based formulations are frequently applied to solubilize or generate supersaturation in gastrointestinal fluids, thus enhancing oral drug absorption. Although many innovative in vitro and in silico tools have been introduced in recent years to aid development of enabling formulations, significant knowledge gaps still exist with respect to how best to implement them. As a result, the development strategy for enabling formulations varies considerably within the industry and many elements of empiricism remain. The InPharma network aims to advance a mechanistic, animal-free approach to the assessment of drug developability. This commentary focuses current status and next steps that will be taken in InPharma to identify and fully utilize 'best practice' in vitro and in silico tools for use in physiologically based biopharmaceutic models.

Prediction of gastric pH-mediated drug exposure using physiologically-based pharmacokinetic modeling: A case study of itraconazole

Abnormal gastric acidity, including achlorhydria, can act as a significant source of variability in orally administered drugs especially with pH-sensitive solubility profiles, such as weak bases, potentially resulting in an undesirable therapeutic response. This study aimed to evaluate the utility of physiologically-based pharmacokinetic (PBPK) modeling in the prediction of gastric pH-mediated drug exposure by using itraconazole, a weak base, as a case. An itraconazole PBPK model was developed on the mechanistic basis of its absorption kinetics in a middle-out manner from a stepwise in vitro-in vivo extrapolation to in vivo refinement. Afterward, an independent prospective clinical study evaluating gastric pH and itraconazole pharmacokinetics (PKs) under normal gastric acidity and esomeprazole-induced gastric hypoacidity was conducted for model validation. Validation was performed by comparing the predicted data with the clinical observations, and the valid model was subsequently applied to predict PK changes under achlorhydria. The developed itraconazole PBPK model showed reasonable reproducibility for gastric pH-mediated exposure observed in the clinical investigation. Based on the model-based simulations, itraconazole exposure was expected to be decreased up to 65% under achlorhydria, and furthermore, gastric pH-mediated exposure could be mechanistically interpreted according to sequential variation in total solubility, dissolution, and absorption. This study suggested the utility of PBPK modeling in the prediction of gastric pH-mediated exposure, especially for drugs whose absorption is susceptible to gastric pH. Our findings will serve as a leading model for further mechanistic assessment of exposure depending on gastric pH for various drugs, ultimately contributing to personalized pharmacotherapy.

Assessment of food effects during clinical development

Food-drug interactions frequently hamper oral drug development due to various physicochemical, physiological and formulation-dependent mechanisms. This has stimulated the development of a range of promising biopharmaceutical assessment tools which, however, lack standardized settings and protocols. Hence, this manuscript aims to provide an overview of the general approach and the methodology used in food effect assessment and prediction. For in vitro dissolution-based predictions, the expected food effect mechanism should be carefully considered when selecting the level of complexity of the model, together with its drawbacks and advantages. Typically, in vitro dissolution profiles are then incorporated into physiologically based pharmacokinetic models, which can estimate the impact of food-drug interactions on bioavailability within 2-fold prediction error, at least. Positive food effects related to drug solubilization in the GI tract are easier to predict than negative food effects. Preclinical animal models also provide a good level of food effect prediction, with beagle dogs remaining the gold standard. When solubility-related food-drug interactions have large clinical impact, advanced formulation approaches can be used to improve fasted state pharmacokinetics, hence decreasing the fasted/fed difference in oral bioavailability. Finally, the knowledge from all studies should be combined to secure regulatory approval of the labelling instructions.

Integrating In Vitro Biopharmaceutics into Physiologically Based Biopharmaceutic Model (PBBM) to Predict Food Effect of BCS IV Zwitterionic Drug (GSK3640254)

A strategy followed to integrate in vitro solubility and permeability data into a PBBM model to predict the food effect of a BCS IV zwitterionic drug (GSK3640254) observed in clinical studies is described. The PBBM model was developed, qualified and verified using clinical data of an immediate release (IR)-tablet (10-320 mg) obtained in healthy volunteers under fasted and fed conditions. The solubility of GSK3640254 was a function of its ionization state, the media composition and pH, whereas its permeability determined using MDCK cell lines was enhanced by the presence of mixed micelles. In vitro data alongside PBBM modelling suggested that the positive food effect observed in the clinical studies was attributed to micelle-mediated enhanced solubility and permeability. The biorelevant media containing oleic acid and cholesterol in fasted and fed levels enabled the model to appropriately capture the magnitude of the food effect. Thus, by using Simcyp^® v20 software, the PBBM model accurately predicted the results of the food effect and predicted data were within a two-fold error with 70% being within 1.25-fold. The developed model strategy can be effectively adopted to increase the confidence of using PBBM models to predict the food effect of BCS class IV drugs.

Mechanistic Modeling of In Vitro Skin Permeation and Extrapolation to In Vivo for Topically Applied Metronidazole Drug Products Using a Physiologically Based Pharmacokinetic Model

Physiologically based pharmacokinetic (PBPK) modeling has increasingly been employed in dermal drug development and regulatory assessment, providing a framework to integrate relevant information including drug and drug product attributes, skin physiology parameters, and population variability. The current study aimed to develop a stepwise modeling workflow with knowledge gained from modeling in vitro skin permeation testing (IVPT) to describe in vivo exposure of metronidazole locally in the stratum corneum following topical application of complex semisolid drug products. The initial PBPK model of metronidazole in vitro skin permeation was developed using infinite and finite dose aqueous metronidazole solution. Parameters such as stratum corneum lipid-water partition coefficient (K_sclip/water) and stratum corneum lipid diffusion coefficient (D_sclip) of metronidazole were optimized using IVPT data from simple aqueous solutions (infinite) and MetroGel (10 mg/cm² dose application), respectively. The optimized model, when parameterized with physical and structural characteristics of the drug products, was able to accurately predict the mean cumulative amount permeated (cm²/h) and flux (μg/cm²/h) profiles of metronidazole following application of different doses of MetroGel and MetroCream. Thus, the model was able to capture the impact of differences in drug product microstructure and metamorphosis of the dosage form on in vitro metronidazole permeation. The PBPK model informed by IVPT study data was able to predict the metronidazole amount in the stratum corneum as reported in clinical studies. In summary, the proposed model provides an enhanced understanding of the potential impact of drug product attributes in influencing in vitro skin permeation of metronidazole. Key kinetic parameters derived from modeling the metronidazole IVPT data improved the predictions of the developed PBPK model of in vivo local metronidazole concentrations in the stratum corneum. Overall, this work improves our confidence in the proposed workflow that accounts for drug product attributes and utilizes IVPT data toward improving predictions from advanced modeling and simulation tools.

Multi-phase multi-layer mechanistic dermal absorption (MPML MechDermA) model to predict local and systemic exposure of drug products applied on skin

Physiologically-based pharmacokinetic models combine knowledge about physiology, drug product properties, such as physicochemical parameters, absorption, distribution, metabolism, excretion characteristics, formulation attributes, and trial design or dosing regimen to mechanistically simulate drug pharmacokinetics (PK). The current work describes the development of a multiphase, multilayer mechanistic dermal absorption (MPML MechDermA) model within the Simcyp Simulator capable of simulating uptake and permeation of drugs through human skin following application of drug products to the skin. The model was designed to account for formulation characteristics as well as body site- and sex- population variability to predict local and systemic bioavailability. The present report outlines the structure and assumptions of the MPML MechDermA model and includes results from simulations comparing absorption at multiple body sites for two compounds, caffeine and benzoic acid, formulated as solutions. Finally, a model of the Feldene (piroxicam) topical gel, 0.5% was developed and assessed for its ability to predict both plasma and local skin concentrations when compared to in vivo PK data.

Guide to development of compound files for PBPK modeling in the Simcyp population-based simulator

The Simcyp Simulator is a software platform for population physiologically‐based pharmacokinetic (PBPK) modeling and simulation. It links in vitro data to in vivo absorption, distribution, metabolism, excretion and pharmacokinetic/pharmacodynamic outcomes to explore clinical scenarios and support drug development decisions, including regulatory submissions and drug labels. This tutorial describes the different input parameters required, as well as the considerations needed when developing a PBPK model within the Simulator, for a small molecule intended for oral administration. A case study showing the development and application of a PBPK model for ondansetron is herein used to aid the understanding of different PBPK model development concepts.

A Mechanistic Absorption and Disposition Model of Ritonavir to Predict Exposure and Drug-Drug Interaction Potential of CYP3A4/5 and CYP2D6 Substrates

BACKGROUND AND OBJECTIVES

Due to health authority warnings and the recommended limited use of ketoconazole as a model inhibitor of cytochrome P450 (CYP) 3A4 in clinical drug-drug interaction (DDI) studies, there is a need to search for alternatives. Ritonavir is a strong inhibitor for CYP3A4/5-mediated DDIs and has been proposed as a suitable alternative to ketoconazole. It can also be used as a weak inhibitor for CYP2D6-mediated DDIs. Most of the currently available physiologically based pharmacokinetic (PBPK) inhibitor models developed for predicting DDIs use first-order absorption models, which do not mechanistically capture the effect of formulations on the systemic exposure of the inhibitor. Thus, the main purpose of the current study was to verify the predictive performance of a mechanistic absorption and disposition model of ritonavir when it was applied to the inhibition of CYP2D6 and CYP3A4/5 by ritonavir.

METHODS

A PBPK model that incorporates formulation characteristics and enzyme kinetic parameters for post-absorptive pharmacokinetic processes of ritonavir was constructed. Key absorption-related parameters in the model were determined using mechanistic modelling of in vitro biopharmaceutics experiments. The model was verified for systemic exposure and DDI risk assessment using clinical observations from 13 and 18 studies, respectively.

RESULTS

Maximal inhibition of hepatic (3.53% of the activity remaining) and gut (5.16% of the activity remaining) CYP3A4 activity was observed when ritonavir was orally administered in doses of 100 mg or higher. The PBPK model accurately described the concentrations of ritonavir in the different simulated studies. The prediction accuracy for maximum concentration (C_max) and area under the plasma concentration versus time curve (AUC) were assessed. The bias (average fold error, AFE) for the prediction of C_max and AUC was 0.92 and 1.06, respectively, and the precision (absolute average fold error, AAFE) was 1.29 and 1.23, respectively. The PBPK model predictions for all C_max and AUC ratios when ritonavir was used as an inhibitor of CYP metabolism fell within twofold of the clinical observations. The prediction accuracy for C_max and AUC ratios had a bias (AFE) of 0.85 and 0.99, respectively, and a precision (AAFE) of 1.21 and 1.33, respectively.

CONCLUSIONS

The current model, which incorporates formulation characteristics and mechanistic disposition parameters, can be used to assess the DDI potential of CYP3A4/5 and CYP2D6 substrates administered with a twice-daily dose of 100 mg of ritonavir for 14 days.

Characterizing the Physicochemical Properties of Two Weakly Basic Drugs and the Precipitates Obtained from Biorelevant Media

Generally, some weakly basic insoluble drugs will undergo precipitate and redissolution after emptying from the stomach to the small intestinal, resulting in the limited ability to predict the absorption characteristics of compounds in advance. Absorption is determined by the solubility and permeability of compounds, which are related to physicochemical properties, while knowledge about the absorption of redissolved precipitate is poorly documented. Considering that biorelevant media have been widely used to simulate gastrointestinal fluids, sufficient precipitates can be obtained in biorelevant media in vitro. Herein, the purpose of this manuscript is to evaluate the physicochemical properties of precipitates obtained from biorelevant media and active pharmaceutical ingredients (API), and then to explore the potential absorption difference between API and precipitates. Precipitates can be formed by the interaction between compounds and intestinal fluid contents, leading to changes in the crystal structure, melting point, and melting process. However, the newly formed crystals have some advantageous properties compared with the API, such as the improved dissolved rate and the increased intrinsic dissolution rate. Additionally, the permeability of some precipitates obtained from biorelevant media was different from API. Meanwhile, the permeability of rivaroxaban and Drug-A was decreased by 1.92-fold and 3.53-fold, respectively, when the experiments were performed in a biorelevant medium instead of a traditional medium. Therefore, the absorption of precipitate may differ from that of API, and the permeability assay in traditional medium may be overestimated. Based on the research results, it is crucial to understand the physicochemical properties of precipitates and API, which can be used as the departure point to improve the prediction performance of absorption.

Biorelevant dissolution testing and physiologically based absorption modeling to predict in vivo performance of supersaturating drug delivery systems

Supersaturating drug delivery systems (SDDS) enhance the oral absorption of poorly water-soluble drugs by achieving a supersaturated state in the gastrointestinal tract. The maintenance of a supersaturated state is decided by the complex interplay among inherent properties of drug, excipients and physiological conditions of gastrointestinal tract. The biopharmaceutical advantage through SDDS can be mechanistically investigated by coupling biopredictive dissolution testing with physiologically based absorption modeling (PBAM). However, the development of biopredictive dissolution methods possess challenges due to concurrent dissolution, supersaturation, precipitation, and possible redissolution of precipitates during gastrointestinal transit of SDDS. In this comprehensive review, our effort is to critically assess the current state-of-knowledge and provide future directions for PBAM of SDDS. The review outlines various methods used to retrieve physiologically relevant values for input parameters like solubility, dissolution, precipitation, lipid-digestion and permeability of SDDS. SDDS-specific parameterization includes solubility values corresponding to apparent physical form, dissolution in physiologically relevant volumes with biorelevant media, and transfer experiments to incorporate precipitation kinetics. Interestingly, the lack of experimental permeability values and modification of absorption flux through SDDS possess the additional challenge for its PBAM. Supersaturation triggered permeability modifications are reported to fit the observed plasma concentration-time profile. Hence, the experimental insights on good fitting with modified permeability can be potential area of future research for the development of in vitro methods to reliably predict oral absorption of SDDS.

Physiologically Based Pharmacokinetic/Pharmacodynamic Modeling to Support Waivers of In Vivo Clinical Studies: Current Status, Challenges, and Opportunities

Physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) modeling has been extensively applied to quantitatively translate in vitro data, predict the in vivo performance, and ultimately support waivers of in vivo clinical studies. In the area of biopharmaceutics and within the context of model-informed drug discovery and development (MID3), there is a rapidly growing interest in applying verified and validated mechanistic PBPK models to waive in vivo clinical studies. However, the regulatory acceptance of PBPK analyses for biopharmaceutics and oral drug absorption applications, which is also referred to variously as "PBPK absorption modeling" [Zhang et al. CPT: Pharmacometrics Syst. Pharmacol. 2017, 6, 492], "physiologically based absorption modeling", or "physiologically based biopharmaceutics modeling" (PBBM), remains rather low [Kesisoglou et al. J. Pharm. Sci. 2016, 105, 2723] [Heimbach et al. AAPS J. 2019, 21, 29]. Despite considerable progress in the understanding of gastrointestinal (GI) physiology, in vitro biopharmaceutic and in silico tools, PBPK models for oral absorption often suffer from an incomplete understanding of the physiology, overparameterization, and insufficient model validation and/or platform verification, all of which can represent limitations to their translatability and predictive performance. The complex interactions of drug substances and (bioenabling) formulations with the highly dynamic and heterogeneous environment of the GI tract in different age, ethnic, and genetic groups as well as disease states have not been yet fully elucidated, and they deserve further research. Along with advancements in the understanding of GI physiology and refinement of current or development of fully mechanistic in silico tools, we strongly believe that harmonization, interdisciplinary interaction, and enhancement of the translational link between in vitro, in silico, and in vivo will determine the future of PBBM. This Perspective provides an overview of the current status of PBBM, reflects on challenges and knowledge gaps, and discusses future opportunities around PBPK/PD models for oral absorption of small and large molecules to waive in vivo clinical studies.

Bottom-Up Physiologically Based Oral Absorption Modeling of Free Weak Base Drugs

In this study, we systematically evaluated "bottom-up" physiologically based oral absorption modeling, focusing on free weak base drugs. The gastrointestinal unified theoretical framework (the GUT framework) was employed as a simple and transparent model. The oral absorption of poorly soluble free weak base drugs is affected by gastric pH. Alternation of bulk and solid surface pH by dissolving drug substances was considered in the model. Simple physicochemical properties such as pKa, the intrinsic solubility, and the bile micelle partition coefficient were used as input parameters. The fraction of a dose absorbed (Fa) in vivo was obtained by reanalyzing the pharmacokinetic data in the literature (15 drugs, a total of 85 Fa data). The AUC ratio with/without a gastric acid-reducing agent (AUCr) was collected from the literature (22 data). When gastric dissolution was neglected, Fa was underestimated (absolute average fold error (AAFE) = 1.85, average fold error (AFE) = 0.64). By considering gastric dissolution, predictability was improved (AAFE = 1.40, AFE = 1.04). AUCr was also appropriately predicted (AAFE = 1.54, AFE = 1.04). The Fa values of several drugs were slightly overestimated (less than 1.7-fold), probably due to neglecting particle growth in the small intestine. This modeling strategy will be of great importance for drug discovery and development.

Hybridized nanoamorphous micellar dispersion using a QbD-DM3 linked rational product design strategy for ritonavir: A BCS IV drug

A QbD-DM³ linked rational product design strategy was adopted to create a hybridized ritonavir (RTV, BCS Class IV) nanoamorphous micellar dispersion (RTV-NAD). A DM³ research strategy was employed in conjunction with the quality-by-design spaces, and quality target product profile to link the critical material attributes and critical process parameters to the quality target product profile's critical product attributes QbD elements. A Box-Behnken design and multivariate analysis using multiple linear regression and partial least squares provided data analysis. The hybridized strategy leveraged three different mechanisms to increase RTV's solubility and four mechanisms to increase its dissolution rate. Statistically significant models were generated for critical product attributes: particle size (p = 0.0000, R² adjusted = 0.9513), polydispersity index (p = 0.0002, R² adjusted = 0.6398), zeta potential (p = 0.0000, R² adjusted = 0.9744), and drug loading on a dry basis (p = 0.0000, R² adjusted = 0.9951). The impact of drug concentration, Soluplus® concentration, and solvent:antisolvent ratio, their interactions and square effects on the critical product attributes were assessed by multivariate analysis. The QbD optimal formulation was determined for RTV-NAD. Multiple linear regression and partial least squares computational predictability was evaluated using three verification batches. The prediction error for critical product attributes was <5%. RTV-NAD and ritonavir microsuspension were characterized by x-ray diffraction and in-vitro dissolution studies. X-ray diffraction confirmed the amorphous nature of the RTV-NAD. RTV-NAD exhibited a 'spring-hover' dissolution profile at pH 4.5. At pH 6.8, a classic 'spring-parachute' dissolution behavior was observed.