Physiologically Based Pharmacokinetic (PBPK) Modeling and Simulation Approaches: A Systematic Review of Published Models, Applications, and Model Verification

Physiologically Based Pharmacokinetic Modeling of Cannabidiol, Delta-9-Tetrahydrocannabinol, and Their Metabolites in Healthy Adults After Administration by Multiple Routes

The two most extensively studied cannabinoids, cannabidiol (CBD) and delta-9-tetrahydrocannabinol (THC), are used for myriad conditions. THC is predominantly eliminated via the cytochromes P450 (CYPs), whereas CBD is eliminated through both CYPs and UDP-glucuronosyltransferases (UGTs). The fractional contributions of these enzymes to cannabinoid metabolism have shown conflicting results among studies. Physiologically based pharmacokinetic (PBPK) models for CBD and THC and for drug-drug interaction studies involving CBD or THC as object drugs were developed and verified to improve estimates of these contributions. First, physicochemical and pharmacokinetic parameters for CBD, THC, and their metabolites (7-OH-CBD, 11-OH-THC, and 11-COOH-THC) were obtained from the literature or optimized. Second, PBPK base models were developed for CBD and THC after intravenous administration. Third, beginning with the intravenous models, absorption models were developed for CBD after oral and oromucosal spray administration and for THC after oral, inhalation, and oromucosal spray administration. The full models well-captured the area under the concentration-time curve (AUC) and peak concentration (Cmax) of CBD and THC from the verification dataset. Predicted AUC and Cmax for CBD and 7-OH-CBD were within two-fold of the observed data. For THC, 11-OH-THC, and 11-COOH-THC, 100%, 100%, and 83% of the predicted AUC values were within two-fold, respectively, of the observed values; 100%, 92%, and 94% of the predicted Cmax values, respectively, were within two-fold of the observed values. The verified models could be used to help address critical public health needs, including assessing potential drug interaction risks involving CBD and THC.

Applications of pharmacometrics in drug development

The last two decades have witnessed profound changes in how advanced computational tools can help leverage tons of data to improve our knowledge, and ultimately reduce cost and increase productivity in drug development. Pharmacometrics has demonstrated its impact through model-informed drug development (MIDD) approaches. It is now an indispensable component throughout the whole continuum of drug discovery, development, regulatory review, and approval. Today, applications of pharmacometrics are common in designing better trials and accelerating evidence-based decisions. Newly emerging technologies, especially those from data and computer sciences, are being integrated with existing computational tools used in the pharmaceutical industry at a remarkably fast pace. The new challenges faced by the pharmacometrics community are not what or how to contribute, but which optimal MIDD strategy should be adopted to maximize its value in the decision-making process. While we are embracing new innovative approaches and tools, this article discusses how a variety of existing modeling tools, with differentiated advantages and focus, can work in concert to inform drug development.

Overcoming Challenges in Small-Molecule Drug Bioavailability: A Review of Key Factors and Approaches

The bioavailability of small-molecule drugs remains a critical challenge in pharmaceutical development, significantly impacting therapeutic efficacy and commercial viability. This review synthesizes recent advances in understanding and overcoming bioavailability limitations, focusing on key physicochemical and biological factors influencing drug absorption and distribution. We examine cutting-edge strategies for enhancing bioavailability, including innovative formulation approaches, rational structural modifications, and the application of artificial intelligence in drug design. The integration of nanotechnology, 3D printing, and stimuli-responsive delivery systems are highlighted as promising avenues for improving drug delivery. We discuss the importance of a holistic, multidisciplinary approach to bioavailability optimization, emphasizing early-stage consideration of ADME properties and the need for patient-centric design. This review also explores emerging technologies such as CRISPR-Cas9-mediated personalization and microbiome modulation for tailored bioavailability enhancement. Finally, we outline future research directions, including advanced predictive modeling, overcoming biological barriers, and addressing the challenges of emerging therapeutic modalities. By elucidating the complex interplay of factors affecting bioavailability, this review aims to guide future efforts in developing more effective and accessible small-molecule therapeutics.

Prediction of human pharmacokinetic parameters incorporating SMILES information

This study aimed to develop a model incorporating natural language processing analysis for the simplified molecular-input line-entry system (SMILES) to predict clearance (CL) and volume of distribution at steady state (Vd,ss) in humans. The construction of CL and Vd,ss prediction models involved data from 435 to 439 compounds, respectively. In machine learning, features such as animal pharmacokinetic data, in vitro experimental data, molecular descriptors, and SMILES were utilized, with XGBoost employed as the algorithm. The ChemBERTa model was used to analyze substance SMILES, and the last hidden layer embedding of ChemBERTa was examined as a feature. The model was evaluated using geometric mean fold error (GMFE), r2, root mean squared error (RMSE), and accuracy within 2- and 3-fold error. The model demonstrated optimal performance for CL prediction when incorporating animal pharmacokinetic data, in vitro experimental data, and SMILES as features, yielding a GMFE of 1.768, an r2 of 0.528, an RMSE of 0.788, with accuracies within 2-fold and 3-fold error reaching 75.8% and 81.8%, respectively. The model's performance in Vd,ss prediction was optimized by leveraging animal pharmacokinetic data and in vitro experimental data as features, yielding a GMFE of 1.401, an r2 of 0.902, an RMSE of 0.413, with accuracies within 2-fold and 3-fold error reaching 93.8% and 100%, respectively. This study has developed a highly predictive model for CL and Vd,ss. Specifically, incorporating SMILES information into the model has predictive power for CL.

Causal, predictive or observational? Different understandings of key event relationships for adverse outcome pathways and their implications on practice

The Adverse Outcome Pathways (AOPs) framework is pivotal in toxicology, but the, terminology describing Key Event Relationships (KERs) varies within AOP guidelines.This study examined the usage of causal, observational and predictive terms in AOP, documentation and their adaptation in AOP development. A literature search and text, analysis of key AOP guidance documents revealed nuanced usage of these terms, with KERs often described as both causal and predictive. The adaptation of, terminology varies across AOP development stages. Evaluation of KER causality often, relies targeted blocking experiments and weight-of-evidence assessments in the, putative and qualitative stages. Our findings highlight a potential mismatch between,terminology in guidelines and methodologies in practice, particularly in inferring,causality from predictive models. We argue for careful consideration of terms like, causal and essential to facilitate interdisciplinary communication. Furthermore, integrating known causality into quantitative AOP models remains a challenge.

Drug dosing optimization in critically ill children under continuous renal replacement therapy: from basic concepts to the bedside model informed precision dosing

INTRODUCTION

Optimizing drug dosage in critically ill children undergoing Continuous Renal Replacement Therapy (CRRT) is mandatory and challenging, given the many factors impacting pharmacokinetics and pharmacodynamics coupled with the vulnerability of this population.

AREAS COVERED

A good understanding of the mechanisms that determine drug elimination via the CRRT technique is useful to avoid prescription pitfalls, however limited by the high between and within subject variability. The developments of population pharmacokinetic and physiologically based pharmacokinetic models derived from in-vivo and in-vitro studies, are challenging, but remain the most appropriate tool to suggest adjusted dosage regimens for every patient, throughout treatment. We searched PubMed using the search string: 'pediatrics OR children' AN 'continuous renal replacement therapy' AND 'pharmacokinetics' AND 'model informed precision dosing' AND, 'physiologically based pharmacokinetics,' AND 'therapeutic drug monitoring' until January 2024, regardless of language or publication status.

EXPERT OPINION

Familiarizing the pediatric intensivists with the therapeutic drug monitoring and providing clinicians the individualized prescribing software such as Model Informed Precision Dosing would be a significant step forward. The clinical benefit for patients remains to be demonstrated.

Advancements in Virtual Bioequivalence: A Systematic Review of Computational Methods and Regulatory Perspectives in the Pharmaceutical Industry

BACKGROUND/OBJECTIVES

The rise of virtual bioequivalence studies has transformed the pharmaceutical landscape, enabling more efficient drug development processes. This systematic review aims to explore advancements in physiologically based pharmacokinetic (PBPK) modeling, its regulatory implications, and its role in achieving virtual bioequivalence, particularly for complex drug formulations.

METHODS

We conducted a systematic review of clinical trials using computational methods, particularly PBPK modeling, to carry out bioequivalence assessments. Eligibility criteria are emphasized during in silico modeling and pharmacokinetic simulations. Comprehensive literature searches were performed across databases such as PubMed, Scopus, and the Cochrane Library. A search strategy using key terms and Boolean operators ensured that extensive coverage was achieved. We adhered to the PRISMA guidelines in regard to the study selection, data extraction, and quality assessment, focusing on key characteristics, methodologies, outcomes, and regulatory perspectives from the FDA and EMA.

RESULTS

Our findings indicate that PBPK modeling significantly enhances the prediction of pharmacokinetic profiles, optimizing dosing regimens, while minimizing the need for extensive clinical trials. Regulatory agencies have recognized this utility, with the FDA and EMA developing frameworks to integrate in silico methods into drug evaluations. However, challenges such as study heterogeneity and publication bias may limit the generalizability of the results.

CONCLUSIONS

This review highlights the critical need for standardized protocols and robust regulatory guidelines to facilitate the integration of virtual bioequivalence methodologies into pharmaceutical practices. By embracing these advancements, the pharmaceutical industry can improve drug development efficiency and patient outcomes, paving the way for innovative therapeutic solutions. Continued research and adaptive regulatory frameworks will be essential in navigating this evolving field.

Physiologically Based Pharmacokinetic Modeling to Assess Ritonavir-Digoxin Interactions and Recommendations for Co-Administration Regimens

BACKGROUND

Digoxin is a commonly used cardiac glycoside drug in clinical practice, primarily transported by P-glycoprotein (P-gp) and susceptible to the influence of P-gp inhibitors/inducers. Concurrent administration of ritonavir and digoxin may significantly increase the plasma concentration of digoxin. Due to the narrow therapeutic window of digoxin, combined use may lead to severe toxic effects.

PURPOSE

Utilize a Physiology-Based Pharmacokinetic (PBPK) model to simulate and predict the impact of the interaction between ritonavir and digoxin on the pharmacokinetics (PK) of digoxin, and provide recommendations for the combined medication regimen.

METHODS

Using PK-Sim®, develop individual PBPK models for ritonavir and digoxin. Simulate the exposure in a drug-drug interaction (DDI) scenario by implementing ritonavir's inhibition of P-glycoprotein (P-gp) on digoxin. Evaluate the performance of the models by comparing the predicted and observed plasma concentration-time curves and predicted versus observed PK parameter values. Finally, adjust the dosing regimen for the combined therapy based on the changes in exposure.

RESULTS

According to the model simulations, the steady-state exposure of digoxin increased by 86.5% and 90.2% for oral administration, and 80.2% and 90.2% for intravenous administration, respectively, when 0.25 mg or 0.5 mg of digoxin was administered concurrently with ritonavir. By reducing the dose of digoxin by 45% or doubling the oral administration interval, similar steady-state concentrations can be achieved compared to when the drugs are not co-administered.

CONCLUSIONS

In clinical practice, the influence of drug interactions on the plasma concentration changes of digoxin within the body should be considered to ensure the safety and effectiveness of treatment.

Physiologically based pharmacokinetic modelling to predict potential drug-drug interactions of dersimelagon (MT-7117)

AIMS

Dersimelagon is a novel, investigational, orally administered, selective agonist of the melanocortin-1 receptor that has demonstrated efficacy at increasing symptom-free light exposure and an acceptable safety profile in patients with protoporphyria. A phase 1 drug-drug interaction (DDI) study demonstrated that dersimelagon 300 mg has the potential for clinically relevant DDIs with drugs that are substrates for breast cancer resistance protein, such as atorvastatin and rosuvastatin. This study uses physiologically based pharmacokinetic (PBPK) modelling to further investigate the DDI effects at lower doses of dersimelagon with substrate drugs.

METHODS

The data from in silico, in vitro and in vivo studies were used to construct a PBPK model for dersimelagon to assess the DDI potential between dersimelagon and substrate drugs for cytochrome P450 3A, P-glycoprotein, organic anion transporting polypeptide 1B1/1B3, organic anion transporter 3 and breast cancer resistance protein, including atorvastatin and rosuvastatin.

RESULTS

The systemic exposure of atorvastatin based on the maximum plasma concentration and area under the plasma concentration-time curve was predicted to increase 1.21-fold and 1.25-fold, respectively, if coadministered with dersimelagon 100 mg, and 1.42-fold and 1.45-fold with dersimelagon 200 mg. The systemic exposure of rosuvastatin followed trends similar to atorvastatin (1.67-fold and 1.34-fold increase in maximum plasma concentration and area under the plasma concentration-time curve, respectively, with dersimelagon 100 mg, and 2.40-fold and 1.69-fold with dersimelagon 200 mg).

CONCLUSION

Overall, PBPK modelling results indicate that the simulated changes in plasma exposure of atorvastatin and rosuvastatin following coadministration with dersimelagon 100 or 200 mg are not clinically significant, but caution and appropriate clinical monitoring should be recommended.

A physiologically based pharmacokinetic model of voriconazole in human CNS-Integrating time-dependent inhibition of CYP3A4, genetic polymorphisms of CYP2C19 and possible transporter mechanisms

OBJECTIVES

Voriconazole is a classical antifungal drug that is often used to treat CNS fungal infections due to its permeability through the BBB. However, its clinical use remains challenging because of its narrow therapeutic window and wide inter-individual variability. In this study, we proposed an optimised and validated PBPK model by integrating in vitro, in vivo and clinical data to simulate the distribution and PK process of voriconazole in the CNS, providing guidance for clinical individualised treatment.

METHODS

The model structure was optimised and tissue-to-plasma partition coefficients were obtained through animal experiments. Using the allometric relationships, the distribution of voriconazole in the human CNS was predicted. The model integrated factors affecting inter-individual variation and drug interactions of voriconazole-polymorphisms in the CYP2C19 gene and auto-inhibition and then was validated using real clinical data.

RESULTS

The overall AFE value showing model predicted differences was 1.1420 in the healthy population; and in the first prediction of plasma and CSF in actual clinical patients, 89.5% of the values were within the 2-fold error interval, indicating good predictive performance of the model. The bioavailability of voriconazole varied at different doses (39%-86%), and the optimised model conformed to this pattern (46%-83%).

CONCLUSIONS

Combined with the relevant pharmacodynamic indexes, the PBPK model provides a feasible way for precise medication in patients with CNS infection and improve the treatment effect and prognosis.

Physiologically based pharmacokinetic model to predict drug-drug interactions with the antibody-drug conjugate enfortumab vedotin

Enfortumab vedotin is an antibody-drug conjugate (ADC) comprised of a Nectin-4-directed antibody and monomethyl auristatin E (MMAE), which is primarily eliminated through P-glycoprotein (P-gp)-mediated excretion and cytochrome P450 3A4 (CYP3A4)-mediated metabolism. A physiologically based pharmacokinetic (PBPK) model was developed to predict effects of combined P-gp with CYP3A4 inhibitor/inducer (ketoconazole/rifampin) on MMAE exposure when coadministered with enfortumab vedotin and study enfortumab vedotin with CYP3A4 (midazolam) and P-gp (digoxin) substrate exposure. A PBPK model was built for enfortumab vedotin and unconjugated MMAE using the PBPK simulator ADC module. A similar model was developed with brentuximab vedotin, an ADC with the same valine-citrulline-MMAE linker as enfortumab vedotin, for MMAE drug-drug interaction (DDI) verification using clinical data. The DDI simulation predicted a less-than-2-fold increase in MMAE exposure with enfortumab vedotin plus ketoconazole (MMAE geometric mean ratio [GMR] for maximum concentration [Cmax], 1.15; GMR for area under the time-concentration curve from time 0 to last quantifiable concentration [AUClast], 1.38). Decreased MMAE exposure above 50% but below 80% was observed with enfortumab vedotin plus rifampin (MMAE GMR Cmax, 0.72; GMR AUClast, 0.47). No effect of enfortumab vedotin on midazolam or digoxin systemic exposure was predicted. Results suggest that combination enfortumab vedotin, P-gp, and a CYP3A4 inhibitor may result in increased MMAE exposure and patients should be monitored for potential adverse effects. Combination P-gp and a CYP3A4 inducer may result in decreased MMAE exposure. No exposure change is expected for CYP3A4 or P-gp substrates when combined with enfortumab vedotin.ClinicalTrials.gov identifier Not applicable.

Seminar: Functional Exposomics and Mechanisms of Toxicity-Insights from Model Systems and NAMs

BACKGROUND

Significant progress has been made over the past decade in measuring the chemical components of the exposome, providing transformative population-scale frameworks in probing the etiologic link between environmental factors and disease phenotypes. While the analytical technologies continue to evolve with reams of data being generated, there is an opportunity to complement exposome-wide association studies (ExWAS) with functional analyses to advance etiologic search at organismal, cellular, and molecular levels.

OBJECTIVES

Exposomics is a transdisciplinary field aimed at enabling discovery-based analysis of the nongenetic factors that contribute to disease, including numerous environmental chemical stressors. While advances in exposure assessment are enhancing population-based discovery of exposome-wide effects and chemical exposure agents, functional screening and elucidation of biological effects of exposures represent the next logical step toward precision environmental health and medicine. In this work, we focus on the use, strategies, and prospects of alternative approaches and model systems to enhance the current human exposomics framework in biomarker search and causal understanding, spanning from bench-based nonmammalian organisms and cell culture to computational new approach methods (NAMs).

DISCUSSION

We visit the definition of the functional exposome and exposomics and discuss a need to leverage alternative models as opposed to mammalian animals for delineating exposome-wide health effects. Under the "three Rs" principle of reduction, replacement, and refinement, model systems such as roundworms, fruit flies, zebrafish, and induced pluripotent stem cells (iPSCs) are advantageous over mammals (e.g., rodents or higher vertebrates). These models are cost-effective, and cell-specific genetic manipulations in these models are easier and faster, compared to mammalian models. Meanwhile, in silico NAMs enhance hazard identification and risk assessment in humans by bridging the translational gaps between toxicology data and etiologic inference, as represented by in vitro to in vivo extrapolation (IVIVE) and integrated approaches to testing and assessment (IATA) under the adverse outcome pathway (AOP) framework. Together, these alternatives offer a strong toolbox to support functional exposomics to study toxicity and causal mediators underpinning exposure-disease links. https://doi.org/10.1289/EHP13120.

Contribution of Continued Dermal Exposure of PFAS-Containing Sunscreens to Internal Exposure: Extrapolation from In Vitro and In Vivo Tests to Physiologically Based Toxicokinetic Models

Per- and polyfluoroalkyl substances (PFASs) are widely present in sunscreen products as either active ingredients or impurities. They may penetrate the human skin barrier and then pose potential health risks. Herein, we aimed to develop a physiologically based toxicokinetic (PBTK) model capable of predicting the body loading of PFASs after repeated, long-term dermal application of commercial sunscreens. Ten laboratory-prepared sunscreens, generally falling into two categories of water-in-oil (W/O) and oil-in-water (O/W) sunscreens, were subject to in vitro percutaneous penetration test to assess the impacts of four sunscreen ingredients on PFAS penetration. According to the results, two sunscreen formulas representing W/O and O/W types that mostly enhanced PFAS dermal absorption were then selected for a subsequent 30 day in vivo exposure experiment in mice. PBTK models were successfully established based on the time-dependent PFAS concentrations in mouse tissues (R2 = 0.885-0.947) and validated through another 30 day repeated exposure experiment in mice using two commercially available sunscreens containing PFASs (R2 = 0.809-0.835). The PBTK model results suggest that applying sunscreen of the same amount on a larger skin area is more conducive to PFAS permeation, thus enhancing the exposure risk. This emphasizes the need for caution in practical sunscreen application scenarios, particularly during the summer months.

Model-informed drug development in pediatric, pregnancy and geriatric drug development: States of the art and future

The challenges of drug development in pediatric, pregnant and geriatric populations are a worldwide concern shared by regulatory authorities, pharmaceutical companies, and healthcare professionals. Model-informed drug development (MIDD) can integrate and quantify real-world data of physiology, pharmacology, and disease processes by using modeling and simulation techniques to facilitate decision-making in drug development. In this article, we reviewed current MIDD policy updates, reflected on the integrity of physiological data used for MIDD and the effects of physiological changes on the drug PK, as well as summarized current MIDD strategies and applications, so as to present the state of the art of MIDD in pediatric, pregnant and geriatric populations. Some considerations are put forth for the future improvements of MIDD including refining regulatory considerations, improving the integrity of physiological data, applying the emerging technologies, and exploring the application of MIDD in new therapies like gene therapies for special populations.

Advances in Modeling Approaches for Oral Drug Delivery: Artificial Intelligence, Physiologically-Based Pharmacokinetics, and First-Principles Models

Oral drug absorption is the primary route for drug administration. However, this process hinges on multiple factors, including the drug's physicochemical properties, formulation characteristics, and gastrointestinal physiology. Given its intricacy and the exorbitant costs associated with experimentation, the trial-and-error method proves prohibitively expensive. Theoretical models have emerged as a cost-effective alternative by assimilating data from diverse experiments and theoretical considerations. These models fall into three categories: (i) data-driven models, encompassing classical pharmacokinetics, quantitative-structure models (QSAR), and machine/deep learning; (ii) mechanism-based models, which include quasi-equilibrium, steady-state, and physiologically-based pharmacokinetics models; and (iii) first principles models, including molecular dynamics and continuum models. This review provides an overview of recent modeling endeavors across these categories while evaluating their respective advantages and limitations. Additionally, a primer on partial differential equations and their numerical solutions is included in the appendix, recognizing their utility in modeling physiological systems despite their mathematical complexity limiting widespread application in this field.

High-Resolution Mass Spectrometry for Human Exposomics: Expanding Chemical Space Coverage

In the modern "omics" era, measurement of the human exposome is a critical missing link between genetic drivers and disease outcomes. High-resolution mass spectrometry (HRMS), routinely used in proteomics and metabolomics, has emerged as a leading technology to broadly profile chemical exposure agents and related biomolecules for accurate mass measurement, high sensitivity, rapid data acquisition, and increased resolution of chemical space. Non-targeted approaches are increasingly accessible, supporting a shift from conventional hypothesis-driven, quantitation-centric targeted analyses toward data-driven, hypothesis-generating chemical exposome-wide profiling. However, HRMS-based exposomics encounters unique challenges. New analytical and computational infrastructures are needed to expand the analysis coverage through streamlined, scalable, and harmonized workflows and data pipelines that permit longitudinal chemical exposome tracking, retrospective validation, and multi-omics integration for meaningful health-oriented inferences. In this article, we survey the literature on state-of-the-art HRMS-based technologies, review current analytical workflows and informatic pipelines, and provide an up-to-date reference on exposomic approaches for chemists, toxicologists, epidemiologists, care providers, and stakeholders in health sciences and medicine. We propose efforts to benchmark fit-for-purpose platforms for expanding coverage of chemical space, including gas/liquid chromatography-HRMS (GC-HRMS and LC-HRMS), and discuss opportunities, challenges, and strategies to advance the burgeoning field of the exposome.

Clinical Ocular Exposure Extrapolation for a Complex Ophthalmic Suspension Using Physiologically Based Pharmacokinetic Modeling and Simulation

The development of generic ophthalmic drug products with complex formulations is challenging due to the complexity of the ocular system and a lack of sensitive testing to evaluate the interplay of its physiology with ophthalmic drugs. New methods are needed to facilitate the development of ophthalmic generic drug products. Ocular physiologically based pharmacokinetic (O-PBPK) models can provide insight into drug partitioning in eye tissues that are usually not accessible and/or are challenging to sample in humans. This study aims to demonstrate the utility of an ocular PBPK model to predict human exposure following the administration of ophthalmic suspension. Besifloxacin (Bes) suspension is presented as a case study. The O-PBPK model for Bes ophthalmic suspension (Besivance® 0.6%) accounts for nasolacrimal drainage, suspended particle dissolution in the tears, ocular absorption, and distribution in the rabbit eye. A topical controlled release formulation was used to integrate the effect of Durasite® on Bes ocular retention. The model was subsequently used to predict Bes exposure after its topical administration in humans. Drug-specific parameters were used as validated for rabbits. The physiological parameters were adjusted to match human ocular physiology. Simulated human ocular pharmacokinetic profiles were compared with the observed ocular tissue concentration data to assess the OCAT models' ability to predict human ocular exposure. The O-PBPK model simulations adequately described the observed concentrations in the eye tissues following the topical administration of Bes suspension in rabbits. After adjustment of physiological parameters to represent the human eye, the extrapolation of clinical ocular exposure following a single ocular administration of Bes suspension was successful.

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.

Advancing Drug Safety in Drug Development: Bridging Computational Predictions for Enhanced Toxicity Prediction

The attrition rate of drugs in clinical trials is generally quite high, with estimates suggesting that approximately 90% of drugs fail to make it through the process. The identification of unexpected toxicity issues during preclinical stages is a significant factor contributing to this high rate of failure. These issues can have a major impact on the success of a drug and must be carefully considered throughout the development process. These late-stage rejections or withdrawals of drug candidates significantly increase the costs associated with drug development, particularly when toxicity is detected during clinical trials or after market release. Understanding drug-biological target interactions is essential for evaluating compound toxicity and safety, as well as predicting therapeutic effects and potential off-target effects that could lead to toxicity. This will enable scientists to predict and assess the safety profiles of drug candidates more accurately. Evaluation of toxicity and safety is a critical aspect of drug development, and biomolecules, particularly proteins, play vital roles in complex biological networks and often serve as targets for various chemicals. Therefore, a better understanding of these interactions is crucial for the advancement of drug development. The development of computational methods for evaluating protein-ligand interactions and predicting toxicity is emerging as a promising approach that adheres to the 3Rs principles (replace, reduce, and refine) and has garnered significant attention in recent years. In this review, we present a thorough examination of the latest breakthroughs in drug toxicity prediction, highlighting the significance of drug-target binding affinity in anticipating and mitigating possible adverse effects. In doing so, we aim to contribute to the development of more effective and secure drugs.

Recent Progress on Physiologically Based Pharmacokinetic (PBPK) Model: A Review Based on Bibliometrics

Physiologically based pharmacokinetic/toxicokinetic (PBPK/PBTK) models are designed to elucidate the mechanism of chemical compound action in organisms based on the physiological, biochemical, anatomical, and thermodynamic properties of organisms. After nearly a century of research and practice, good results have been achieved in the fields of medicine, environmental science, and ecology. However, there is currently a lack of a more systematic review of progress in the main research directions of PBPK models, especially a more comprehensive understanding of the application in aquatic environmental research. In this review, a total of 3974 articles related to PBPK models from 1996 to 24 March 2024 were collected. Then, the main research areas of the PBPK model were categorized based on the keyword co-occurrence maps and cluster maps obtained by CiteSpace. The results showed that research related to medicine is the main application area of PBPK. Four major research directions included in the medical field were "drug assessment", "cross-species prediction", "drug-drug interactions", and "pediatrics and pregnancy drug development", in which "drug assessment" accounted for 55% of the total publication volume. In addition, bibliometric analyses indicated a rapid growth trend in the application in the field of environmental research, especially in predicting the residual levels in organisms and revealing the relationship between internal and external exposure. Despite facing the limitation of insufficient species-specific parameters, the PBPK model is still an effective tool for improving the understanding of chemical-biological effectiveness and will provide a theoretical basis for accurately assessing potential risks to ecosystems and human health. The combination with the quantitative structure-activity relationship model, Bayesian method, and machine learning technology are potential solutions to the previous research gaps.

Investigating clinical pharmacokinetics of brivaracetam by using a pharmacokinetic modeling approach

The development of technology and the processing speed of computing machines have facilitated the evaluation of advanced pharmacokinetic (PK) models, making modeling processes simple and faster. The present model aims to analyze the PK of brivaracetam (BRV) in healthy and diseased populations. A comprehensive literature review was conducted to incorporate the BRV plasma concentration data and its input parameters into PK-Sim software, leading to the creation of intravenous (IV) and oral models for both populations. The developed physiologically based pharmacokinetic (PBPK) model of BRV was then assessed using the visual predictive checks, mean observed/predicted ratios (Robs/pre), and average fold error for PK parameters including the maximum systemic concentration (Cmax), the area under the curve at time 0 to t (AUC0-∞), and drug clearance (CL). The PBPK model of BRV demonstrated that mean Robs/pre ratios of the PK parameters remained within the acceptable limits when assessed against a twofold error margin. Furthermore, model predictions were carried out to assess how AUC0-∞ is affected following the administration of BRV in individuals with varying degrees of liver cirrhosis, ranging from different child-pugh (CP) scores like A, B, and C. Moreover, dose adjustments were recommended by considering the variations in Cmax and CL in various kidney disease stages (mild to severe).

Investigating pharmacokinetic profiles of Centella asiatica using machine learning and PBPK modelling

Physiologically based pharmacokinetic (PBPK) modeling serves as a valuable tool for determining the distribution and disposition of substances in the body of an organism. It involves a mathematical representation of the interrelationships among crucial physiological, biochemical, and physicochemical parameters. A lack of the values of pharmacokinetic parameters can be challenging in constructing a PBPK model. Herein, we propose an artificial intelligence framework to evaluate a key pharmacokinetic parameter, the intestinal effective permeability (Peff). The publicly available Peff dataset was utilized to develop regression machine learning models. The XGBoost model demonstrates the best test accuracy of R-squared (R2, coefficient of determination) of 0.68. The model is then applied to compute the Peff of asiaticoside and madecassoside, the parent compounds found in Centella asiatica. Subsequently, PBPK modeling was conducted to evaluate the biodistribution of the herbal substances following oral administration in a rat model. The simulation results were evaluated and validated, which agreed with the existing in vivo studies in rats. This in silico pipeline presents a potential approach for investigating the pharmacokinetic parameters and profiles of drugs or herbal substances, which can be used independently or integrated into other modeling systems.

Another string to your bow: machine learning prediction of the pharmacokinetic properties of small molecules

INTRODUCTION

Prediction of pharmacokinetic (PK) properties is crucial for drug discovery and development. Machine-learning (ML) models, which use statistical pattern recognition to learn correlations between input features (such as chemical structures) and target variables (such as PK parameters), are being increasingly used for this purpose. To embed ML models for PK prediction into workflows and to guide future development, a solid understanding of their applicability, advantages, limitations, and synergies with other approaches is necessary.

AREAS COVERED

This narrative review discusses the design and application of ML models to predict PK parameters of small molecules, especially in light of established approaches including in vitro-in vivo extrapolation (IVIVE) and physiologically based pharmacokinetic (PBPK) models. The authors illustrate scenarios in which the three approaches are used and emphasize how they enhance and complement each other. In particular, they highlight achievements, the state of the art and potentials of applying machine learning for PK prediction through a comphrehensive literature review.

EXPERT OPINION

ML models, when carefully crafted, regularly updated, and appropriately used, empower users to prioritize molecules with favorable PK properties. Informed practitioners can leverage these models to improve the efficiency of drug discovery and development process.

DMPK perspective on quantitative model analysis for chimeric antigen receptor cell therapy: Advances and challenges

Chimeric antigen receptor (CAR) cells are genetically engineered immune cells that specifically target tumor-associated antigens and have revolutionized cancer treatment, particularly in hematological malignancies, with ongoing investigations into their potential applications in solid tumors. This review provides a comprehensive overview of the current status and challenges in drug metabolism and pharmacokinetics (DMPK) for CAR cell therapy, specifically emphasizing on quantitative modeling and simulation (M&S). Furthermore, the recent advances in quantitative model analysis have been reviewed, ranging from clinical data characterization to mechanism-based modeling that connects in vitro and in vivo nonclinical and clinical study data. Additionally, the future perspectives and areas for improvement in CAR cell therapy translation have been reviewed. This includes using formulation quality considerations, characterization of appropriate animal models, refinement of in vitro models for bottom-up approaches, and enhancement of quantitative bioanalytical methodology. Addressing these challenges within a DMPK framework is pivotal in facilitating the translation of CAR cell therapy, ultimately enhancing the patients' lives through efficient CAR cell therapies.

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.

Warfarin-A natural anticoagulant: A review of research trends for precision medication

BACKGROUND

Warfarin is a widely prescribed anticoagulant in the clinic. It has a more considerable individual variability, and many factors affect its variability. Mathematical models can quantify the quantitative impact of these factors on individual variability.

PURPOSE

The aim is to comprehensively analyze the advanced warfarin dosing algorithm based on pharmacometrics and machine learning models of personalized warfarin dosage.

METHODS

A bibliometric analysis of the literature retrieved from PubMed and Scopus was performed using VOSviewer. The relevant literature that reported the precise dosage of warfarin calculation was retrieved from the database. The multiple linear regression (MLR) algorithm was excluded because a recent systematic review that mainly reviewed this algorithm has been reported. The following terms of quantitative systems pharmacology, mechanistic model, physiologically based pharmacokinetic model, artificial intelligence, machine learning, pharmacokinetic, pharmacodynamic, pharmacokinetics, pharmacodynamics, and warfarin were added as MeSH Terms or appearing in Title/Abstract into query box of PubMed, then humans and English as filter were added to retrieve the literature.

RESULTS

Bibliometric analysis revealed important co-occuring MeShH and index keywords. Further, the United States, China, and the United Kingdom were among the top countries contributing in this domain. Some studies have established personalized warfarin dosage models using pharmacometrics and machine learning-based algorithms. There were 54 related studies, including 14 pharmacometric models, 31 artificial intelligence models, and 9 model evaluations. Each model has its advantages and disadvantages. The pharmacometric model contains biological or pharmacological mechanisms in structure. The process of pharmacometric model development is very time- and labor-intensive. Machine learning is a purely data-driven approach; its parameters are more mathematical and have less biological interpretation. However, it is faster, more efficient, and less time-consuming. Most published models of machine learning algorithms were established based on cross-sectional data sourced from the database.

CONCLUSION

Future research on personalized warfarin medication should focus on combining the advantages of machine learning and pharmacometrics algorithms to establish a more robust warfarin dosage algorithm. Randomized controlled trials should be performed to evaluate the established algorithm of warfarin dosage. Moreover, a more user-friendly and accessible warfarin precision medicine platform should be developed.

Construction of a physiologically based pharmacokinetic model of paclobutrazol and exposure estimation in the human body

Paclobutrazol (PBZ) is a plant growth regulator that can delay plant growth and improve plant resistance and yield. Although it has been widely used in the growth of medicinal plants, human beings may take it by taking traditional Chinese medicine. There are no published studies on PBZ exposure in humans or standardized limits for PBZ in medicinal plants. We measured the solubility, oil-water partition coefficient (logP), and pharmacokinetics of PBZ in rats and established a physiologically based pharmacokinetic (PBPK) model of PBZ in rats. This was followed by extrapolation to healthy Chinese adult males as a theoretical foundation for future risk assessment of PBZ. The results showed that PBZ had low solubility and high fat solubility. Pharmacokinetic experiments showed that PBZ was absorbed rapidly but eliminated slowly in rats. On this basis, the rat PBPK model was successfully constructed and extrapolated to healthy Chinese adult males to predict the plasma concentration-time curve and exposure of PBZ in humans. The construction of the PBPK model of PBZ in this study facilitates the determination of the standard formulation limits and risk assessment of PBZ residues in medicinal plants.

Development and verification of mechanistic vaginal absorption and metabolism model to predict systemic exposure after vaginal ring and gel application

AIMS

The current work describes the development of mechanistic vaginal absorption and metabolism model within Simcyp Simulator to predict systemic concentrations following vaginal application of ring and gel formulations.

METHODS

Vaginal and cervix physiology parameters were incorporated in the model development. The study highlights the model assumptions including simulation results comparing systemic concentrations of 5 different compounds, namely, dapivirine, tenofovir, lidocaine, ethinylestradiol and etonogestrel, administered as vaginal ring or gel. Due to lack of data, the vaginal absorption parameters were calculated based on assumptions or optimized. The model uses release rate/in vitro release profiles with formulation characteristics to predict drug mass transfer across vaginal tissue into the systemic circulation.

RESULTS

For lidocaine and tenofovir vaginal gel, the predicted to observed AUC0-t and Cmax ratios were well within 2-fold error limits. The average fold error (AFE) and absolute AFE indicating bias and precision of predictions range from 0.62 to 1.61. For dapivirine, the pharmacokinetic parameters are under and overpredicted in some studies due to lack of formulation composition details and relevance of release rate used in ring model. The predicted to observed AUC0-t and Cmax ratios were well within 2-fold error limits for etonogestrel and ethinylestradiol vaginal ring (AFEs and absolute AFEs from 0.84 to 1.83).

CONCLUSION

The current study provides first of its kind physiologically based pharmacokinetic framework integrating physiology, population and formulation data to carry out in silico mechanistic vaginal absorption studies, with the potential for virtual bioequivalence assessment in the future.

In Vitro In Vivo Extrapolation and Bioequivalence Prediction for Immediate-Release Capsules of Cefadroxil Based on a Physiologically-Based Pharmacokinetic ACAT Model

Physiologically based pharmacokinetic (PBPK) modeling is a mechanistic concept, which helps to judge the effects of biopharmceutical properties of drug product such as in vitro dissolution on its pharmacokinetic and in vivo performance. With the application of virtual bioequivalence (VBE) study, the drug product development using model-based approach can help in evaluating the possibility of extending BCS-based biowaiver. Therefore, the current study was intended to develop PBPK model as well as in vitro in vivo extrapolation (IVIVE) for BCS class III drug i.e. cefadroxil. A PBPK model was created in GastroPlus™ 9.8.3 utilizing clinical data of immediate-release cefadroxil formulations. By the examination of simulated and observed plasma drug concentration profiles, the predictability of the proposed model was assessed for the prediction errors. Furthermore, mechanistic deconvolution was used to create IVIVE, and the plasma drug concentration profiles and pharmacokinetic parameters were predicted for different virtual formulations with variable cefadroxil in vitro release. Virtual bioequivalence study was also executed to assess the bioequivalence of the generic verses the reference drug product (Duricef®). The developed PBPK model satisfactorily predicted Cmax and AUC0-t after cefadroxil single and multiple oral dose administrations, with all individual prediction errors within the limits except in a few cases. Second order polynomial correlation function obtained accurately predict in vivo drug release and plasma concentration profile of cefadroxil test and reference (Duricef®) formulation. The VBE study also proved test formulation bioequivalent to reference formulation and the statistical analysis on pharmacokinetic parameters reported 90% confidence interval for Cmax and AUC0-t in the FDA acceptable limits. The analysis found that a validated and verified PBPK model with a mechanistic background is as a suitable approach to accelerate generic drug development.

Physiologically based pharmacokinetic modeling to predict the pharmacokinetics of codeine in different CYP2D6 phenotypes

Objectives

Codeine, a prodrug used as an opioid agonist, is metabolized to the active product morphine by CYP2D6. This study aimed to establish physiologically based pharmacokinetic (PBPK) models of codeine and morphine and explore the influence of CYP2D6 genetic polymorphisms on the pharmacokinetics of codeine and morphine.

Methods

An initial PBPK modeling of codeine in healthy adults was established using PK-Sim® software and subsequently extrapolated to CYP2D6 phenotype-related PBPK modeling based on the turnover frequency (Kcat) of CYP2D6 for different phenotype populations (UM, EM, IM, and PM). The mean fold error (MFE) and geometric mean fold error (GMFE) methods were used to compare the differences between the predicted and observed values of the pharmacokinetic parameters to evaluate the accuracy of PBPK modeling. The validated models were then used to support dose safety for different CYP2D6 phenotypes.

Results

The developed and validated CYP2D6 phenotype-related PBPK model successfully predicted codeine and morphine dispositions in different CYP2D6 phenotypes. Compared with EMs, the predicted AUC0-∞ value of morphine was 98.6% lower in PMs, 60.84% lower in IMs, and 73.43% higher in UMs. Morphine plasma exposure in IMs administered 80 mg of codeine was roughly comparable to that in EMs administered 30 mg of codeine. CYP2D6 UMs may start dose titration to achieve an optimal individual regimen and avoid a single dose of over 20 mg. Codeine should not be used in PMs for pain relief, considering its insufficient efficacy.

Conclusion

PBPK modeling can be applied to explore the dosing safety of codeine and can be helpful in predicting the effect of CYP2D6 genetic polymorphisms on drug-drug interactions (DDIs) with codeine in the future.

Physiologically Based Pharmacokinetic Modeling for Maribavir to Inform Dosing in Drug-Drug Interaction Scenarios with CYP3A4 Inducers and Inhibitors

Maribavir, an orally available antiviral agent, has been approved in multiple countries for the treatment of patients with refractory post-transplant cytomegalovirus (CMV) infection and/or disease. Maribavir is primarily metabolized by CYP3A4; coadministration with CYP3A4 inducers and inhibitors may significantly alter maribavir exposure, thereby affecting its efficacy and safety. The effect of CYP3A4 inducers and inhibitors on maribavir exposure was evaluated based on a drug-drug interaction (DDI) study and physiologically-based pharmacokinetic (PBPK) modeling. The effect of rifampin (a strong inducer of CYP3A4 and moderate inducer of CYP1A2), administered at a 600 mg dose once daily, on maribavir pharmacokinetics was assessed in a clinical phase 1 DDI study in healthy participants. A full PBPK model for maribavir was developed and verified using in vitro and clinical pharmacokinetic data from phase 1 studies. The verified PBPK model was then used to simulate maribavir DDI interactions with various CYP3A4 inducers and inhibitors. The DDI study results showed that coadministration with rifampin decreased the maribavir maximum plasma concentration (Cmax), area under the plasma concentration-time curve (AUC), and trough concentration (Ctrough) by 39%, 60%, and 82%, respectively. Based on the results from the clinical DDI study, the coadministration of maribavir with rifampin is not recommended. The PBPK model did not predict a clinically significant effect of CYP3A4 inhibitors on maribavir exposure; however, it predicted that strong or moderate CYP3A4 inducers, including carbamazepine, efavirenz, phenobarbital, and phenytoin, may reduce maribavir exposure to a clinically significant extent, and may prompt the consideration of a maribavir dosing increase, in accordance with local approved labels and/or regulations.

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.

Physiologically-based pharmacokinetic modelling guided dose evaluations of nirmatrelvir/ritonavir in renal impairment for the management of COVID-19

We aimed to address factors contributing to the pharmacokinetic changes of nirmatrelvir/ritonavir in renal impaired (RI) patients and recommend dosing adjustment via a physiologically-based pharmacokinetic (PBPK) modelling approach. A PBPK model of nirmatrelvir/ritonavir was developed via Simcyp® Simulator. Sensitivity analysis of the influence of hepatic CYP3A4 intrinsic clearance and abundance, as well as hepatic non-CYP3A4 metabolism (other human liver microsomes [HLM] CLint) was performed to evaluate the effects of RI on oral clearance of nirmatrelvir. Other HLM CLint, the most sensitive parameter, was adjusted, and the simulated plasma concentration profiles of nirmatrelvir in severe RI subjects were within the therapeutic index of 292-10 000 ng/mL for dosing regimens of loading doses of 300/100 mg followed by 150/100 mg or 75/100 mg twice daily of nirmatrelvir/ritonavir. Considering that nirmatrelvir is available as a 150 mg tablet, we recommend 300/100 mg followed by 150/100 mg twice daily as the dosing regimen to be investigated in severe RI.

The Effect of Inflammatory Bowel Disease and Irritable Bowel Syndrome on Pravastatin Oral Bioavailability: In vivo and in silico evaluation using bottom-up wbPBPK modeling

The common disorders irritable bowel syndrome (IBS) and inflammatory bowel disease (IBD) can modify the drugs' pharmacokinetics via their induced pathophysiological changes. This work aimed to investigate the impact of these two diseases on pravastatin oral bioavailability. Rat models for IBS and IBD were used to experimentally test the effects of IBS and IBD on pravastatin pharmacokinetics. Then, the observations made in rats were extrapolated to humans using a mechanistic whole-body physiologically-based pharmacokinetic (wbPBPK) model. The rat in vivo studies done herein showed that IBS and IBD decreased serum albumin (> 11% for both), decreased PRV binding in plasma, and increased pravastatin absolute oral bioavailability (0.17 and 0.53 compared to 0.01) which increased plasma, muscle, and liver exposure. However, the wbPBPK model predicted muscle concentration was much lower than the pravastatin toxicity thresholds for myotoxicity and rhabdomyolysis. Overall, IBS and IBD can significantly increase pravastatin oral bioavailability which can be due to a combination of increased pravastatin intestinal permeability and decreased pravastatin gastric degradation resulting in higher exposure. This is the first study in the literature investigating the effects of IBS and IBD on pravastatin pharmacokinetics. The high interpatient variability in pravastatin concentrations as induced by IBD and IBS can be reduced by oral administration of pravastatin using enteric-coated tablets. Such disease (IBS and IBD)-drug interaction can have more drastic consequences for narrow therapeutic index drugs prone to gastric degradation, especially for drugs with low intestinal permeability.

Nano- and microplastic PBK modeling in the context of human exposure and risk assessment

Insufficient data on nano- and microplastics (NMP) hinder robust evaluation of their potential health risks. Methodological disparities and the absence of established toxicity thresholds impede the comparability and practical application of research findings. The diverse attributes of NMP, such as variations in sizes, shapes, and compositions, complicate human health risk assessment. Although probability density functions (PDFs) show promise in capturing this diversity, their integration into risk assessment frameworks is limited. Physiologically based kinetic (PBK) models offer a potential solution to bridge the gap between external exposure and internal dosimetry for risk evaluation. However, the heterogeneity of NMP poses challenges for accurate biodistribution modeling. A literature review, encompassing both experimental and modeling studies, was conducted to examine biodistribution studies of monodisperse micro- and nanoparticles. The literature search in PubMed and Scopus databases yielded 39 studies that met the inclusion criteria. Evaluation criteria were adapted from previous Quality Assurance and Quality Control (QA-QC) studies, best practice guidelines from WHO (2010), OECD guidance (2021), and additional criteria specific to NMP risk assessment. Subsequently, a conceptual framework for a comprehensive NMP-PBK model was developed, addressing the multidimensionality of NMP particles. Parameters for an NMP-PBK model are presented. QA-QC evaluations revealed that most experimental studies scored relatively well (>0) in particle characterizations and environmental settings but fell short in criteria application for biodistribution modeling. The evaluation of modeling studies revealed that information regarding the model type and allometric scaling requires improvement. Three potential applications of PDFs in PBK modeling of NMP are identified: capturing the multidimensionality of the NMP continuum, quantifying the probabilistic definition of external exposure, and calculating the bio-accessibility fraction of NMP in the human body. A framework for an NMP-PBK model is proposed, integrating PDFs to enhance the assessment of NMP's impact on human health.

Research Landscape of Physiologically Based Pharmacokinetic Model Utilization in Different Fields: A Bibliometric Analysis (1999-2023)

PURPOSE

The physiologically based pharmacokinetic (PBPK) modeling has received increasing attention owing to its excellent predictive abilities. However, there has been no bibliometric analysis about PBPK modeling. This research aimed to summarize the research development and hot points in PBPK model utilization overall through bibliometric analysis.

METHODS

We searched for publications related to the PBPK modeling from 1999 to 2023 in the Web of Science Core Collection (WoSCC) database. The Microsoft Office Excel, CiteSpace and VOSviewers were used to perform the analyses.

RESULTS

A total of 4,649 records from 1999 to 2023 were identified, and the largest number of publications focused in the period 2018-2023. The United States was the leading country, and the Environmental Protection Agency (EPA) was the leading institution. The journal Drug Metabolism and Disposition published and co-cited the most articles. Drug-drug interactions, special populations, and new drug development are the main topics in this research field.

CONCLUSION

We first visualize the research landscape and hotspots of the PBPK modeling through bibliometric methods. Our study provides a better understanding for researchers, especially beginners about the dynamization of PBPK modeling and presents the relevant trend in the future.

Physiologically Based Pharmacokinetic Modeling to Unravel the Drug-gene Interactions of Venlafaxine: Based on Activity Score-dependent Metabolism by CYP2D6 and CYP2C19 Polymorphisms

BACKGROUND

Venlafaxine (VEN) is a commonly utilized medication for alleviating depression and anxiety disorders. The presence of genetic polymorphisms gives rise to considerable variations in plasma concentrations across different phenotypes. This divergence in phenotypic responses leads to notable differences in both the efficacy and tolerance of the drug.

PURPOSE

A physiologically based pharmacokinetic (PBPK) model for VEN and its metabolite O-desmethylvenlafaxine (ODV) to predict the impact of CYP2D6 and CYP2C19 gene polymorphisms on VEN pharmacokinetics (PK).

METHODS

The parent-metabolite PBPK models for VEN and ODV were developed using PK-Sim® and MoBi®. Leveraging prior research, derived and implemented CYP2D6 and CYP2C19 activity score (AS)-dependent metabolism to simulate exposure in the drug-gene interactions (DGIs) scenarios. The model's performance was evaluated by comparing predicted and observed values of plasma concentration-time (PCT) curves and PK parameters values.

RESULTS

In the base models, 91.1%, 94.8%, and 94.6% of the predicted plasma concentrations for VEN, ODV, and VEN + ODV, respectively, fell within a twofold error range of the corresponding observed concentrations. For DGI scenarios, these values were 81.4% and 85% for VEN and ODV, respectively. Comparing CYP2D6 AS = 2 (normal metabolizers, NM) populations to AS = 0 (poor metabolizers, PM), 0.25, 0.5, 0.75, 1.0 (intermediate metabolizers, IM), 1.25, 1.5 (NM), and 3.0 (ultrarapid metabolizers, UM) populations in CYP2C19 AS = 2.0 group, the predicted DGI AUC0-96 h ratios for VEN were 3.65, 3.09, 2.60, 2.18, 1.84, 1.56, 1.34, 0.61, and for ODV, they were 0.17, 0.35, 0.51, 0.64, 0.75, 0.83, 0.90, 1.11, and the results were similar in other CYP2C19 groups. It should be noted that PK differences in CYP2C19 phenotypes were not similar across different CYP2D6 groups.

CONCLUSIONS

In clinical practice, the impact of genotyping on the in vivo disposition process of VEN should be considered to ensure the safety and efficacy of treatment.

Challenging the status quo: a framework for mechanistic and human-relevant cardiovascular safety screening

Traditional approaches to preclinical drug safety assessment have generally protected human patients from unintended adverse effects. However, these assessments typically occur too late to make changes in the formulation or in phase 1 and beyond, are highly dependent on animal studies and have the potential to lead to the termination of useful drugs due to liabilities in animals that are not applicable in patients. Collectively, these elements come at great detriment to both patients and the drug development sector. This phenomenon is particularly problematic in the area of cardiovascular safety assessment where preclinical attrition is high. We believe that a more efficient and translational approach can be defined. A multi-tiered assessment that leverages our understanding of human cardiovascular biology, applies human cell-based in vitro characterizations of cardiovascular responses to insult, and incorporates computational models of pharmacokinetic relationships would enable earlier and more translational identification of human-relevant liabilities. While this will take time to develop, the ultimate goal would be to implement such assays both in the lead selection phase as well as through regulatory phases.

Towards More Robust Evaluation of the Predictive Performance of Physiologically Based Pharmacokinetic Models: Using Confidence Intervals to Support Use of Model-Informed Dosing in Clinical Care

BACKGROUND AND OBJECTIVE

With the rise in the use of physiologically based pharmacokinetic (PBPK) modeling over the past decade, the use of PBPK modeling to underpin drug dosing for off-label use in clinical care has become an attractive option. In order to use PBPK models for high-impact decisions, thorough qualification and validation of the model is essential to gain enough confidence in model performance. Currently, there is no agreed method for model acceptance, while clinicians demand a clear measure of model performance before considering implementing PBPK model-informed dosing. We aim to bridge this gap and propose the use of a confidence interval for the predicted-to-observed geometric mean ratio with predefined boundaries. This approach is similar to currently accepted bioequivalence testing procedures and can aid in improved model credibility and acceptance.

METHODS

Two different methods to construct a confidence interval are outlined, depending on whether individual observations or aggregate data are available from the clinical comparator data sets. The two testing procedures are demonstrated for an example evaluation of a midazolam PBPK model. In addition, a simulation study is performed to demonstrate the difference between the twofold criterion and our proposed method.

RESULTS

Using midazolam adult pharmacokinetic data, we demonstrated that creating a confidence interval yields more robust evaluation of the model than a point estimate, such as the commonly used twofold acceptance criterion. Additionally, we showed that the use of individual predictions can reduce the number of required test subjects. Furthermore, an easy-to-implement software tool was developed and is provided to make our proposed method more accessible.

CONCLUSIONS

With this method, we aim to provide a tool to further increase confidence in PBPK model performance and facilitate its use for directly informing drug dosing in clinical care.

Physiologically Based Pharmacokinetic Modeling of Lacosamide in Patients With Hepatic and Renal Impairment and Pediatric Populations to Support Pediatric Dosing Optimization

PURPOSE

Lacosamide (LCM) is a new-generation anti-seizure medication that is efficacious in patients with focal seizures with or without secondary generalization. Until now, the efficacy, safety, and tolerability of LCM are still lacking in Chinese epilepsy patients, particularly for pediatric populations and patients with renal or hepatic impairment.

METHODS

This study was conducted to develop a physiologically based pharmacokinetic (PBPK) model to characterize the pharmacokinetics of LCM in Chinese populations and predict the pharmacokinetics of LCM in Chinese pediatric populations and patients with renal or hepatic impairment. Using data from clinical investigations, the developed PBPK model was validated by comparing predicted and observed blood concentration data.

FINDINGS

Doses should be reduced to approximately 82%, 75%, 63%, and 76% of the Chinese healthy adult dose in patients with mild, moderate, and severe renal impairment and end-stage renal disease; and approximately 89%, 72%, and 36% of the Chinese healthy adult dose in patients with Child Pugh-A, B, and C hepatic impairment. For pediatric populations, intravenous doses should be adjusted to 1.75 mg/kg for newborns, 2.5 mg/kg for toddlers, 2.2 mg/kg mg for preschool and school age, and 2 mg/kg mg for adolescents to achieve an equivalent plasma exposure of 2 mg/kg LCM in adults. The oral doses should be adjusted to 20 mg for toddlers, 32 mg for preschool, 45 mg for school age, and 95 mg for adolescents to achieve an approximately equivalent plasma exposure of 100 mg LCM in adults.

IMPLICATIONS

The PBPK model of LCM can be utilized to optimize dosage regimens for special populations.

Multi-channel Small Animal Drug Metabolism Real-Time Monitoring Fluorescence System

PURPOSE

The data acquisition of drug metabolism analysis requires a lot of time and animal resources. However, there are often many deviations in the results of pharmacokinetic analysis. Conventional methods cannot measure the blood drug concentration data in multiple tissues at the same time, and the data is obtained by in vitro measurement, which produces time errors, in vitro data errors, and individual differences between animals. In the analysis of pharmacokinetic parameters, it will seriously affect the pass rate of clinical trials of R&D drugs and the accuracy of the dosing schedule. To the best of our knowledge, we have not found the study of in vivo blood drug concentration using multi-channel equipment. Therefore, the purpose of this paper is to build a set of multi-organ monitoring and analysis instruments for synchronously monitoring the metabolism of drugs in various tissues of small animals, so as to obtain real in vivo data of blood drug concentration in real time.

PROCEDURES

Using the fluorescence properties and laser-induced fluorescence principle of drugs, we designed six channels to monitor the changes of fluorescence-labeled drugs in their main metabolic organs, a multi-channel calibration method was proposed to improve the accuracy of the time-division multiplexing, the real-time collection of drug concentration in vivo is realized, and the drug metabolism curve in vivo can be observed.

RESULTS

The instrument satisfies the collection of small doses of drugs such as microgram; the detection sensitivity can reach 10 ng/ml, and can monitor and collect the drug metabolism of multiple small animal tissues at the same time, which greatly reduces the use of animals, reduces the differences between individuals, and reduces consumption cost and improve the detection efficiency of parameters, and obtain data information that is closer to the real biology.

CONCLUSION

The real-time continuous monitoring and data collection of the drug metabolism in the plasma of living small animals and the important organs such as kidney, liver, and spleen were realized. The research and development of new drugs and clinical research have higher practical value.

Physiologically based pharmacokinetic modeling of apixaban to predict exposure in populations with hepatic and renal impairment and elderly populations

BACKGROUND

Apixaban is a factor Xa inhibitor with a limited therapeutic index that belongs to the family of oral direct anticoagulants. The pharmacokinetic (PK) behavior of apixaban may be altered in elderly populations and populations with renal or hepatic impairment, necessitating dosage adjustments.

METHODS

This study was conducted to examine how the physiologically based pharmacokinetic (PBPK) model describes the PKs of apixaban in adult and elderly populations and to determine the PKs of apixaban in elderly populations with renal and hepatic impairment. After PBPK models were constructed using the reported physicochemical properties of apixaban and clinical data, they were validated using data from clinical studies involving various dose ranges. Comparing predicted and observed blood concentration data and PK parameters was utilized to evaluate the model's fit performance.

RESULTS

Doses should be reduced to approximately 70% of the healthy adult population for the healthy elderly population to achieve the same PK exposure; approximately 88%, 71%, and 89% of that for the elderly populations with mild, moderate, and severe renal impairment, respectively; and approximately 96%, 81%, and 58% of that for the Child Pugh-A, Child Pugh-B, and Child Pugh-C hepatic impairment elderly populations, respectively to achieve the same PK exposure.

CONCLUSION

The findings indicate that the renal and hepatic function might be considered for apixaban therapy in Chinese elderly patients and the PBPK model can be used to optimize dosage regimens for specific populations.

Model-Informed Drug Development: In Silico Assessment of Drug Bioperformance following Oral and Percutaneous Administration

The pharmaceutical industry has faced significant changes in recent years, primarily influenced by regulatory standards, market competition, and the need to accelerate drug development. Model-informed drug development (MIDD) leverages quantitative computational models to facilitate decision-making processes. This approach sheds light on the complex interplay between the influence of a drug's performance and the resulting clinical outcomes. This comprehensive review aims to explain the mechanisms that control the dissolution and/or release of drugs and their subsequent permeation through biological membranes. Furthermore, the importance of simulating these processes through a variety of in silico models is emphasized. Advanced compartmental absorption models provide an analytical framework to understand the kinetics of transit, dissolution, and absorption associated with orally administered drugs. In contrast, for topical and transdermal drug delivery systems, the prediction of drug permeation is predominantly based on quantitative structure-permeation relationships and molecular dynamics simulations. This review describes a variety of modeling strategies, ranging from mechanistic to empirical equations, and highlights the growing importance of state-of-the-art tools such as artificial intelligence, as well as advanced imaging and spectroscopic techniques.

Physiologically based radiopharmacokinetic (PBRPK) modeling to simulate and analyze radiopharmaceutical therapies: studies of non-linearities, multi-bolus injections, and albumin binding

BACKGROUND

We aimed to develop a publicly shared computational physiologically based pharmacokinetic (PBPK) model to reliably simulate and analyze radiopharmaceutical therapies (RPTs), including probing of hot-cold ligand competitions as well as alternative injection scenarios and drug designs, towards optimal therapies.

RESULTS

To handle the complexity of PBPK models (over 150 differential equations), a scalable modeling notation called the "reaction graph" is introduced, enabling easy inclusion of various interactions. We refer to this as physiologically based radiopharmacokinetic (PBRPK) modeling, fine-tuned specifically for radiopharmaceuticals. As three important applications, we used our PBRPK model to (1) study the effect of competition between hot and cold species on delivered doses to tumors and organs at risk. In addition, (2) we evaluated an alternative paradigm of utilizing multi-bolus injections in RPTs instead of prevalent single injections. Finally, (3) we used PBRPK modeling to study the impact of varying albumin-binding affinities by ligands, and the implications for RPTs. We found that competition between labeled and unlabeled ligands can lead to non-linear relations between injected activity and the delivered dose to a particular organ, in the sense that doubling the injected activity does not necessarily result in a doubled dose delivered to a particular organ (a false intuition from external beam radiotherapy). In addition, we observed that fractionating injections can lead to a higher payload of dose delivery to organs, though not a differential dose delivery to the tumor. By contrast, we found out that increased albumin-binding affinities of the injected ligands can lead to such a differential effect in delivering more doses to tumors, and this can be attributed to several factors that PBRPK modeling allows us to probe.

CONCLUSIONS

Advanced computational PBRPK modeling enables simulation and analysis of a variety of intervention and drug design scenarios, towards more optimal delivery of RPTs.

Computational Nanotoxicology Models for Environmental Risk Assessment of Engineered Nanomaterials

Although engineered nanomaterials (ENMs) have tremendous potential to generate technological benefits in numerous sectors, uncertainty on the risks of ENMs for human health and the environment may impede the advancement of novel materials. Traditionally, the risks of ENMs can be evaluated by experimental methods such as environmental field monitoring and animal-based toxicity testing. However, it is time-consuming, expensive, and impractical to evaluate the risk of the increasingly large number of ENMs with the experimental methods. On the contrary, with the advancement of artificial intelligence and machine learning, in silico methods have recently received more attention in the risk assessment of ENMs. This review discusses the key progress of computational nanotoxicology models for assessing the risks of ENMs, including material flow analysis models, multimedia environmental models, physiologically based toxicokinetics models, quantitative nanostructure-activity relationships, and meta-analysis. Several challenges are identified and a perspective is provided regarding how the challenges can be addressed.

Minimally sufficient experimental design using identifiability analysis

Mathematical models are increasingly being developed and calibrated in tandem with data collection, empowering scientists to intervene in real time based on quantitative model predictions. Well-designed experiments can help augment the predictive power of a mathematical model but the question of when to collect data to maximize its utility for a model is non-trivial. Here we define data as model-informative if it results in a unique parametrization, assessed through the lens of practical identifiability. The framework we propose identifies an optimal experimental design (how much data to collect and when to collect it) that ensures parameter identifiability (permitting confidence in model predictions), while minimizing experimental time and costs. We demonstrate the power of the method by applying it to a modified version of a classic site-of-action pharmacokinetic/pharmacodynamic model that describes distribution of a drug into the tumor microenvironment (TME), where its efficacy is dependent on the level of target occupancy in the TME. In this context, we identify a minimal set of time points when data needs to be collected that robustly ensures practical identifiability of model parameters. The proposed methodology can be applied broadly to any mathematical model, allowing for the identification of a minimally sufficient experimental design that collects the most informative data.

Quetiapine-Related Deaths: In Search of a Surrogate Endpoint

Quetiapine is a second-generation antipsychotic drug available for two and half decades. Due to increased misuse, prescription outside the approved indications, and availability on the black market, it is being encountered in medicolegal autopsies more frequently. For instance, it has been linked to increased mortality rates, most likely due to its adverse effects on the cardiovascular system. Its pharmacokinetic features and significant postmortem redistribution challenge traditional sampling in forensic toxicology. Therefore, a systematic literature review was performed, inclusive of PubMed, the Web of Science-core collection, and the Scopus databases; articles were screened for the terms "quetiapine", "death", and "autopsy" to reevaluate each matrix used as a surrogate endpoint in the forensic toxicology of quetiapine-related deaths. Ultimately, this review considers the results of five studies that were well presented (more than two matrices, data available for all analyses, for instance). The highest quetiapine concentrations were usually measured in the liver tissue. As interpreted by their authors, the results of the considered studies showed a strong correlation between some matrices, but, unfortunately, the studies presented models with poor goodness of fit. The distribution of quetiapine in distinct body compartments/tissues showed no statistically significant relationship with the length of the postmortem interval. Furthermore, this study did not confirm the anecdotal correlation of peripheral blood concentrations with skeletal muscle concentrations. Otherwise, there was no consistency regarding selecting an endpoint for analysis.

Solubility-enabling formulations for oral delivery of lipophilic drugs: considering the solubility-permeability interplay for accelerated formulation development

INTRODUCTION

Tackling low water solubility of drug candidates is a major challenge in today's pharmaceutics/biopharmaceutics, especially by means of modern solubility-enabling formulations. However, drug absorption from these formulations oftentimes remains unchanged or even decreases, despite substantial solubility enhancement.

AREAS COVERED

In this article, we overview the simultaneous effects of the formulation on the solubility and the apparent permeability of the drug, and analyze the contribution of this solubility-permeability interplay to the success/failure of the formulation to increase the overall absorption and bioavailability. Three different patterns of interplay were identified: (1) solubility-permeability tradeoff in which every solubility gain comes with a price of concomitant permeability loss; (2) an advantageous interplay pattern in which the permeability remains unchanged alongside the solubility gain; and (3) an optimal interplay pattern in which the formulation increases both the solubility and the permeability. Passive vs. active intestinal permeability considerations in the context of the solubility-permeability interplay are also thoroughly discussed.

EXPERT OPINION

The solubility-permeability interplay pattern of a given formulation has a critical effect on its overall success/failure, and hence, taking into account both parameters in solubility-enabling formulation development is prudent and highly recommended.

Chronic kidney disease and physiologically based pharmacokinetic modeling: a critical review of existing models

INTRODUCTION

Physiologically based pharmacokinetic (PBPK) modeling is a paradigm shift in this era for determining the exposure of drugs in pediatrics, geriatrics, and patients with chronic diseases where clinical trials are difficult to conduct.

AREAS COVERED

This review has collated data regarding published PBPK models on chronic kidney disease (CKD), including the drug and system-specific input model parameters and model evaluation criteria. Four databases were used from 13th June 2023 to 10th July 2023 for identifying the relevant studies that met the inclusion/exclusion criteria. Alterations in plasma protein (albumin/alpha-1 acid glycoprotein), gastric emptying time, hematocrit, small intestinal transit time, the abundance of cytochrome (CYP) 450 enzymes, glomerular filtration rate, and physicochemical parameters for different drugs were explicitly elaborated from earlier reported studies. Moreover, model evaluation depicted that models in CKD for most of the included drugs were within the allowed two-fold error range.

EXPERT OPINION

This review will provide insights for researchers on applying PBPK models in managing patients with different levels of CKD to prevent undesirable side effects and increase the effectiveness of drug therapy.

In Silico and In Vivo Pharmacokinetic Evaluation of 84-B10, a Novel Drug Candidate against Acute Kidney Injury and Chronic Kidney Disease

Acute kidney injury (AKI) and chronic kidney disease (CKD) have become public health problems due to high morbidity and mortality. Currently, drugs recommended for patients with AKI or CKD are extremely limited, and candidates based on a new mechanism need to be explored. 84-B10 is a novel 3-phenylglutaric acid derivative that can activate the mitochondrial protease, Lon protease 1 (LONP1), and may protect against cisplatin-induced AKI and unilateral ureteral obstruction- or 5/6 nephrectomy [5/6Nx]-induced CKD model. Preclinical studies have shown that 84-B10 has a good therapeutic effect, low toxicity, and is a good prospect for further development. In the present study, the UHPLC-MS/MS method was first validated then applied to the pharmacokinetic study and tissue distribution of 84-B10 in rats. Physicochemical properties of 84-B10 were then acquired in silico. Based on these physicochemical and integral physiological parameters, a physiological based pharmacokinetic (PBPK) model was developed using the PK-Sim platform. The fitting accuracy was estimated with the obtained experimental data. Subsequently, the validated model was employed to predict the pharmacokinetic profiles in healthy and chronic kidney injury patients to evaluate potential clinical outcomes. Cmax in CKD patients was about 3250 ng/mL after a single dose of 84-B10 (0.41 mg/kg), and Cmax,ss was 1360 ng/mL after multiple doses. This study may serve in clinical dosage setting in the future.