Application of physiologically based pharmacokinetic modeling in the prediction of pharmacokinetics of bicyclol controlled-release formulation in human

Formulation Strategy of BCS-II Drugs by Coupling Mechanistic In-Vitro and Nonclinical In-Vivo Data with PBPK: Fundamentals of Absorption-Dissolution to Parameterization of Modelling and Simulation

BCS class II candidates pose challenges in drug development due to their low solubility and permeability. Researchers have explored various techniques; co-amorphous and solid dispersion are major approaches to enhance in-vitro drug solubility and dissolution. However, in-vivo oral bioavailability remains challenging. Physiologically based pharmacokinetic (PBPK) modeling with a detailed understanding of drug absorption, distribution, metabolism, and excretion (ADME) using a mechanistic approach is emerging. This review summarizes the fundamentals of the PBPK, dissolution-absorption models, parameterization of oral absorption for BCS class II drugs, and provides information about newly emerging artificial intelligence/machine learning (AI/ML) linked PBPK approaches with their advantages, disadvantages, challenges and areas of further exploration. Additionally, the fully integrated workflow for formulation design for investigational new drugs (INDs) and virtual bioequivalence for generic molecules falling under BCS-II are discussed.

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

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

SeDeM tool-driven full factorial design for osmotic drug delivery of tramadol HCl: Formulation development, physicochemical evaluation, and in-silico PBPK modeling for predictive pharmacokinetic evaluation using GastroPlus™

The study is based on using SeDeM expert system in developing controlled-release tramadol HCl osmotic tablets and its in-silico physiologically based pharmacokinetic (PBPK) modeling for in-vivo pharmacokinetic evaluation. A Quality by Design (QbD) based approach in developing SeDEM-driven full factorial osmotic drug delivery was applied. A 2⁴ Full-factorial design was used to make the trial formulations of tramadol HCl osmotic tablets using NaCl as osmogen, Methocel K4M as rate controlling polymer, and avicel pH 101 as diluent. The preformulation characteristics of formulations (F1-F16) were determined by applying SeDeM Expert Tool. The formulation was optimized followed by in-vivo predictive pharmacokinetic assessment using PBPK “ACAT” model of GastroPlus™. The FTIR results showed no interaction among the ingredients. The index of good compressibility (ICG) values of all trial formulation blends were ≥5, suggesting direct compression is the best-suited method. Formulation F3 and F4 were optimized based on drug release at 2, 10, and 16 h with a zero-order kinetic release (r² = 0.992 and 0.994). The SEM images confirmed micropores formation on the surface of the osmotic tablet after complete drug release. F3 and F4 were also stable (shelf life 29.41 and 23.46 months). The in vivo simulation of the pharmacokinetics of the PBPK in-silico model revealed excellent relative bioavailability of F3 and F4 with reference to tramadol HCl 50 mg IR formulations. The SeDeM expert tool was best utilized to evaluate the compression characteristics of selected formulation excipients and their blends for direct compression method in designing once-daily osmotically controlled-release tramadol HCl tablets. The in-silico GastroPlus™ PBPK modeling provided a thorough pharmacokinetic assessment of the optimized formulation as an alternative to tramadol HCl in vivo studies.

Opportunities and challenges of physiologically based pharmacokinetic modeling in drug delivery

Physiologically based pharmacokinetic (PBPK) modeling is an important in silico tool to bridge drug properties and in vivo PK behaviors during drug development. Over the recent decade, the PBPK method has been largely applied to drug delivery systems (DDS), including oral, inhaled, transdermal, ophthalmic, and complex injectable products. The related therapeutic agents have included small-molecule drugs, therapeutic proteins, nucleic acids, and even cells. Simulation results have provided important insights into PK behaviors of new dosage forms, which strongly support drug regulation. In this review, we comprehensively summarize recent progress in PBPK applications in drug delivery, which shows large opportunities for facilitating drug development. In addition, we discuss the challenges of applying this methodology from a practical viewpoint.

Physiologically based pharmacokinetic/pharmacodynamic modeling to evaluate the absorption of midazolam rectal gel

OBJECTIVE

We aimed to 1) develop physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) models of a novel midazolam rectal gel in healthy adults, 2) assess the contribution of different physiologically relevant factors in rectal absorption, and 3) to provide supports for future clinical studies of midazolam rectal gel.

METHODS

We developed the rectal PBPK model after built the intravenous and the oral PBPK model. Then, the physiological progress of rectal route was described in terms of the drug release, the rectal absorption and the particle first-pass elimination. Next, the validated PBPK model was combined with the sigmoid E_max PD model. This PBPK/PD model was used to identify the dose range and the critical parameters to ensure safety sedation.

RESULTS

Based on the simulations, the recommended maximum dose for adults' sedation was 15 mg. And the retention time of midazolam rectal gel should be longer than 3 h to reach over 80% pharmacokinetics and pharmacodynamics effects.

CONCLUSION

We successfully developed a PBPK/PD model for the midazolam rectal gel, which accurately described the PK/PD behavior in healthy adults and indicated the transit time of rectum was the most sensitive parameter for absorption. This PBPK/PD model would be expected to support the future clinical studies and pediatric application.

Integrated computer-aided formulation design: A case study of andrographolide/ cyclodextrin ternary formulation

Current formulation development strongly relies on trial-and-error experiments in the laboratory by pharmaceutical scientists, which is time-consuming, high cost and waste materials. This research aims to integrate various computational tools, including machine learning, molecular dynamic simulation and physiologically based absorption modeling (PBAM), to enhance andrographolide (AG) /cyclodextrins (CDs) formulation design. The lightGBM prediction model we built before was utilized to predict AG/CDs inclusion's binding free energy. AG/γ-CD inclusion complexes showed the strongest binding affinity, which was experimentally validated by the phase solubility study. The molecular dynamic simulation was used to investigate the inclusion mechanism between AG and γ-CD, which was experimentally characterized by DSC, FTIR and NMR techniques. PBAM was applied to simulate the in vivo behavior of the formulations, which were validated by cell and animal experiments. Cell experiments revealed that the presence of D-α-Tocopherol polyethylene glycol succinate (TPGS) significantly increased the intracellular uptake of AG in MDCK-MDR1 cells and the absorptive transport of AG in MDCK-MDR1 monolayers. The relative bioavailability of the AG-CD-TPGS ternary system in rats was increased to 2.6-fold and 1.59-fold compared with crude AG and commercial dropping pills, respectively. In conclusion, this is the first time to integrate various computational tools to develop a new AG-CD-TPGS ternary formulation with significant improvement of aqueous solubility, dissolution rate and bioavailability. The integrated computational tool is a novel and robust methodology to facilitate pharmaceutical formulation design.

Qualification and application of liquid chromatography-quadrupole time-of-flight mass spectrometric method for the determination of carisbamate in rat plasma and prediction of its human pharmacokinetics using physiologically based pharmacokinetic modeling

Carisbamate is an antiepileptic drug and it also has broad neuroprotective activity and anticonvulsant reaction. In this study, a liquid chromatography-quadrupole time-of-flight mass spectrometric (LC-qTOF-MS) method was developed and applied for the determination of carisbamate in rat plasma to support in vitro and in vivo studies. A quadratic regression (weighted 1/concentration²), with an equation y = ax² + bx + c, was used to fit calibration curves over the concentration range from 9.05 to 6,600 ng/mL for carisbamate in rat plasma. Preclinical in vitro and in vivo studies of carisbamate have been studied through the developed bioanalytical method. Based on these study results, human pharmacokinetic (PK) profile has been predicted using physiologically based pharmacokinetic (PBPK) modeling. The PBPK model was optimized and validated by using the in vitro and in vivo data. The human PK of carisbamate after oral dosing of 750 mg was simulated by using this validated PBPK model. The human PK parameters and profiles predicted from the validated PBPK model were similar to the clinical data. This PBPK model developed from the preclinical data for carisbamate would be useful for predicting the PK of carisbamate in various clinical settings.

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

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

Prediction of a Therapeutic Dose for Buagafuran, a Potent Anxiolytic Agent by Physiologically Based Pharmacokinetic/Pharmacodynamic Modeling Starting from Pharmacokinetics in Rats and Human

Physiologically based pharmacokinetic (PBPK)/pharmacodynamic (PD) models can contribute to animal-to-human extrapolation and therapeutic dose predictions. Buagafuran is a novel anxiolytic agent and phase I clinical trials of buagafuran have been completed. In this paper, a potentially effective dose for buagafuran of 30 mg t.i.d. in human was estimated based on the human brain concentration predicted by a PBPK/PD modeling. The software GastroPlus^TM was used to build the PBPK/PD model for buagafuran in rat which related the brain tissue concentrations of buagafuran and the times of animals entering the open arms in the pharmacological model of elevated plus-maze. Buagafuran concentrations in human plasma were fitted and brain tissue concentrations were predicted by using a human PBPK model in which the predicted plasma profiles were in good agreement with observations. The results provided supportive data for the rational use of buagafuran in clinic.

Scientific and Regulatory Considerations in Solid Oral Modified Release Drug Product Development

This review presents scientific and regulatory considerations for the development of solid oral modified release (MR) drug products. It includes a rationale for patient-focused development based on Quality-by-Design (QbD) principles. Product and process understanding of MR products includes identification and risk-based evaluation of critical material attributes (CMAs), critical process parameters (CPPs), and their impact on critical quality attributes (CQAs) that affect the clinical performance. The use of various biopharmaceutics tools that link the CQAs to a predictable and reproducible clinical performance for patient benefit is emphasized. Product and process understanding lead to a more comprehensive control strategy that can maintain product quality through the shelf life and the lifecycle of the drug product. The overall goal is to develop MR products that consistently meet the clinical objectives while mitigating the risks to patients by reducing the probability and increasing the detectability of CQA failures.

Evaluating a physiologically based pharmacokinetic model for predicting the pharmacokinetics of midazolam in Chinese after oral administration

AIM

To evaluate the SimCYP simulator ethnicity-specific population model for predicting the pharmacokinetics of midazolam, a typical CYP3A4/5 substrate, in Chinese after oral administration.

METHODS

The physiologically based pharmacokinetic (PBPK) model for midazolam was developed using a SimCYP population-based simulator incorporating Chinese population demographic, physiological and enzyme data. A clinical trial was conducted in 40 Chinese subjects (the half was females) receiving a single oral dose of 15 mg midazolam. The subjects were separated into 4 groups based on age (20-50, 51-65, 66-75, and above 76 years), and the pharmacokinetics profiles of each age- and gender-group were determined, and the results were used to verify the PBPK model.

RESULTS

Following oral administration, the simulated profiles of midazolam plasma concentrations over time in virtual Chinese were in good agreement with the observed profiles, as were AUC and Cmax. Moreover, for subjects of varying ages (20-80 years), the ratios of predicted to observed clearances were between 0.86 and 1.12.

CONCLUSION

The SimCYP PBPK model accurately predicted the pharmacokinetics of midazolam in Chinese from youth to old age. This study may provide novel insight into the prediction of CYP3A4/5-mediated pharmacokinetics in the Chinese population relative to Caucasians and other ethnic groups, which can support the rational design of bridging clinical trials.