Search

Several factors can affect drug release from polylactide coglycolide (PLGA)-based formulations, including polymer and drug properties, formulation components, manufacturing processes, and environmental in vitro or in vivo conditions. To achieve optimal release profiles for specific drug delivery applications, it is crucial to understand the mechanistic processes that determine drug release from PLGA-based formulations. In the current study, we developed a mechanistic model for the in vitro drug release of PLGA-based solid implants. The model accounts for all known critical quality attributes (CQAs) and considers the most important release rate processes, including water or dissolution medium influx into the porous structure of the implant, initial noncatalytic hydrolysis of PLGA, autocatalytic hydrolysis, dissolution of oligomers and monomers into the aqueous medium, the liberation of the trapped solid drug from the polymer matrix, dissolution of the solid drug into the wetted pore network, diffusion of the dissolved drug out of the implant, and distribution of the dissolved drug into the dissolution medium. The model has been validated using in vitro release data obtained from implants of four drugs (buserelin, afamelanotide, brimonidine, and nafarelin). The model presented in this manuscript provides valuable insights into the kinetics and mechanism of drug release from PLGA-based solid implants and has demonstrated the potential for optimizing formulation design. The in vitro release model, coupled with physiologically based pharmacokinetic (PBPK) modeling, can predict the in vivo performance of implants and can be used to support bioequivalence studies in a drug development program.

Background/Objectives: Plasma levels of 4β-hydroxycholesterol (4β-OHC), a CYP3A-specific metabolite of cholesterol, are elevated after administration of CYP3A inducers like rifampicin and carbamazepine. To simulate such plasma 4β-OHC increase, we developed a physiologically based pharmacokinetic (PBPK) model of cholesterol and 4β-OHC in the Simcyp PBPK Simulator (Version 23, Certara UK Ltd.) using a middle-out approach. Methods: Relevant physicochemical properties and metabolic pathway data for CYP3A and CYP27A1 was incorporated in the model. Results: The PBPK model recovered the observed baseline plasma 4β-OHC levels in Caucasian, Japanese, and Korean populations. The model also captured the higher baseline 4β-OHC levels in females compared to males, indicative of sex-specific differences in CYP3A abundance. More importantly, the model recapitulated the increased 4β-OHC plasma levels after multiple-dose rifampicin treatment in six independent studies, indicative of hepatic CYP3A induction. The verified model also captured the altered 4β-OHC levels in CYP3A4/5 polymorphic populations and with other CYP3A inducers. The model is limited by scant data on relative contributions of CYP3A and CYP27A1 pathways and does not account for regulatory mechanisms that control plasma cholesterol and 4β-OHC levels. Conclusion: This study provides a quantitative fit-for-purpose and framed-for-future modelling framework for an endogenous biomarker to evaluate the DDI risk with hepatic CYP3A induction.

The Simcyp Simulator is a software platform widely used in the pharmaceutical industry to conduct stochastic physiologically-based pharmacokinetic (PBPK) modeling. This approach has the advantage of combining routinely generated in vitro data on drugs and drug products with knowledge of biology and physiology parameters to predict a priori potential pharmacokinetic changes in absorption, distribution, metabolism, and excretion for populations of interest. Combining such information with pharmacodynamic knowledge of drugs enables planning for potential dosage adjustment when clinical studies are feasible. Although the conduct of dedicated clinical studies in some patient groups (e.g., with hepatic or renal diseases) is part of the regulatory path for drug development, clinical studies for all permutations of covariates potentially affecting pharmacokinetics are impossible to perform. The role of PBPK in filling the latter gap is becoming more appreciated. This tutorial describes the different input parameters required for the creation of a virtual population giving robust predictions of likely changes in pharmacokinetics. It also highlights the considerations needed to qualify the models for such contexts of use. Two case studies showing the step-by-step development and application of population files for obese or morbidly obese patients and individuals with Crohn's disease are provided as the backbone of our tutorial to give some hands-on and real-world examples., (© 2024 The Author(s). CPT: Pharmacometrics & Systems Pharmacology published by Wiley Periodicals LLC on behalf of American Society for Clinical Pharmacology and Therapeutics.)

Introduction: In virtual bioequivalence (VBE) assessments, pharmacokinetic models informed with in vitro data and verified with small clinical trials' data are used to simulate otherwise unfeasibly large trials. Simulated VBE trials are assessed in a frequentist framework as if they were real despite the unlimited number of virtual subjects they can use. This may adequately control consumer risk but imposes unnecessary risks on producers. We propose a fully Bayesian model-integrated VBE assessment framework that circumvents these limitations., Methods: We illustrate our approach with a case study on a hypothetical paliperidone palmitate (PP) generic long-acting injectable suspension formulation using a validated population pharmacokinetic model published for the reference formulation. BE testing, study power, type I and type II error analyses or their Bayesian equivalents, and safe-space analyses are demonstrated., Results: The fully Bayesian workflow is more precise than the frequentist workflow. Decisions about bioequivalence and safe space analyses in the two workflows can differ markedly because the Bayesian analyses are more accurate., Discussion: A Bayesian framework can adequately control consumer risk and minimize producer risk . It rewards data gathering and model integration to make the best use of prior information. The frequentist approach is less precise but faster to compute, and it can still be used as a first step to narrow down the parameter space to explore in safe-space analyses., Competing Interests: FB and CB are employees and potential shareholders of Certara Inc., (Copyright © 2024 Bois and Brochot.)

Background: Understanding the variability across the human population with respect to toxicodynamic responses after exposure to chemicals, such as environmental toxicants or drugs, is essential to define safety factors for risk assessment to protect the entire population. Activation of cellular stress response pathways are early adverse outcome pathway (AOP) key events of chemical-induced toxicity and would elucidate the estimation of population variability of toxicodynamic responses., Objectives: We aimed to map the variability in cellular stress response activation in a large panel of primary human hepatocyte (PHH) donors to aid in the quantification of toxicodynamic interindividual variability to derive safety uncertainty factors., Methods: High-throughput transcriptomics of over 8,000 samples in total was performed covering a panel of 50 individual PHH donors upon 8 to 24 h exposure to broad concentration ranges of four different toxicological relevant stimuli: tunicamycin for the unfolded protein response (UPR), diethyl maleate for the oxidative stress response (OSR), cisplatin for the DNA damage response (DDR), and tumor necrosis factor alpha ( TNF α ) for NF- κ B signaling. Using a population mixed-effect framework, the distribution of benchmark concentrations (BMCs) and maximum fold change were modeled to evaluate the influence of PHH donor panel size on the correct estimation of interindividual variability for the various stimuli., Results: Transcriptome mapping allowed the investigation of the interindividual variability in concentration-dependent stress response activation, where the average of BMCs had a maximum difference of 864-, 13-, 13-, and 259-fold between different PHHs for UPR, OSR, DDR, and NF- κ B signaling-related genes, respectively. Population modeling revealed that small PHH panel sizes systematically underestimated the variance and gave low probabilities in estimating the correct human population variance. Estimated toxicodynamic variability factors of stress response activation in PHHs based on this dataset ranged between 1.6 and 6.3., Discussion: Overall, by combining high-throughput transcriptomics and population modeling, improved understanding of interindividual variability in chemical-induced activation of toxicity relevant stress pathways across the human population using a large panel of plated cryopreserved PHHs was established, thereby contributing toward increasing the confidence of in vitro -based prediction of adverse responses, in particular hepatotoxicity. https://doi.org/10.1289/EHP11891.

The most common method for establishing bioequivalence (BE) is to demonstrate similarity of concentration-time profiles in the systemic circulation, as a surrogate to the site of action. However, similarity of profiles from two formulations in the systemic circulation does not imply similarity in the gastrointestinal tract (GIT) nor local BE. We have explored the concordance of BE conclusions for a set of hypothetical formulations based on budesonide concentration profiles in various segments of gut vs. those in systemic circulation using virtual trials powered by physiologically based pharmacokinetic (PBPK) models. The impact of Crohn's disease on the BE conclusions was explored by changing physiological and biological GIT attributes. Substantial 'discordance' between local and systemic outcomes of VBE was observed. Upper GIT segments were much more sensitive to formulation changes than systemic circulation, where the latter led to false conclusions for BE. The ileum and colon showed a lower frequency of discordance. In the case of Crohn's disease, a product-specific similarity factor might be needed for products such as Entocort ® EC to ensure local BE. Our results are specific to budesonide, but we demonstrate potential discordances between the local gut vs. systemic BE for the first time.

SimRFlow is a high-throughput physiologically based pharmacokinetic (PBPK) modelling tool which uses Certara's Simcyp® simulator. The workflow is comprised of three main modules: 1) a Data Collection module for automated curation of physicochemical (from ChEMBL and the Norman Suspect List databases) and experimental data (i.e.: clearance, plasma-protein binding, and blood-to-plasma ratio, from httk -R package databases), 2) a Simulation module which activates the Simcyp® simulator and runs Monte Carlo simulations on virtual subjects using the curated data, and 3) a Data Visualisation module for understanding the simulated compound-specific profiles and predictions. SimRFlow has three administration routes (oral, intravenous, dermal) and allows users to change some simulation parameters including the number of subjects, simulation duration, and dosing. Users are only expected to provide a file of the compounds they wish to simulate, and in return the workflow provides summary statistics, concentration-time profiles of various tissue types, and a database file (containing in-depth results) for each simulated compound. This is presented within a guided and easy-to-use R Shiny interface which provides many plotting options for the visualisation of concentration-time profiles, parameter distributions, trends between the different parameters, as well as comparison of predicted parameters across all batch-simulated compounds. The in-built R functions can be assembled in user-customised scripts which allows for the modification of the workflow for different purposes. SimRFlow proves to be a time-efficient tool for simulating a large number of compounds without any manual curation of physicochemical or experimental data necessary to run Simcyp® simulations., Competing Interests: All authors were employees of Certara UK Limited and may own shares in the company., (Copyright © 2022 Khalidi, Onasanwo, Islam, Jo, Fisher, Aidley, Gardner and Bois.)

Tizanidine, a centrally acting skeletal muscle relaxant, is predominantly metabolized by CYP1A2 and undergoes extensive hepatic first-pass metabolism after oral administration. As a highly extracted drug, the systemic exposure to tizanidine exhibits considerable interindividual variability and is altered substantially when coadministered with CYP1A2 inhibitors or inducers. The aim of the current study was to compare the performance of a permeability-limited multicompartment liver (PerMCL) model, which operates as an approximation of the dispersion model, and the well stirred model (WSM) for predicting tizanidine drug-drug interactions (DDIs). Physiologically based pharmacokinetic models were developed for tizanidine, incorporating the PerMCL model and the WSM, respectively, to simulate the interaction of tizanidine with a range of CYP1A2 inhibitors and inducers. Whereas the WSM showed a tendency to underpredict the fold change of tizanidine area under the plasma concentration-time curve (AUC ratio) in the presence of perpetrators, the use of PerMCL model increased precision (absolute average-fold error: 1.32-1.42 versus 1.58) and decreased bias (average-fold error: 0.97-1.25 versus 0.63) for the predictions of mean AUC ratios as compared with the WSM. The PerMCL model captured the observed range of individual AUC ratios of tizanidine as well as the correlation between individual AUC ratios and CYP1A2 activities without interactions, whereas the WSM was not able to capture these. The results demonstrate the advantage of using the PerMCL model over the WSM in predicting the magnitude and interindividual variability of DDIs for a highly extracted sensitive substrate tizanidine. SIGNIFICANCE STATEMENT: This study demonstrates the advantages of the PerMCL model, which operates as an approximation of the dispersion model, in mitigating the tendency of the WSM to underpredict the magnitude and variability of DDIs of a highly extracted CYP1A2 substrate tizanidine when it is administered with CYP1A2 inhibitors or inducers. The physiologically based pharmacokinetic modeling approach described herein is valuable to the understanding of drug interactions of highly extracted substrates and the source of its interindividual variability., (Copyright © 2022 by The Author(s).)

Adverse Outcome Pathways (AOPs) are increasingly used to support the integration of in vitro data in hazard assessment for chemicals. Quantitative AOPs (qAOPs) use mathematical models to describe the relationship between key events (KEs). In this paper, data obtained in three cell lines, LHUMES, HepG2 and RPTEC/TERT1, using similar experimental protocols, was used to calibrate a qAOP of mitochondrial toxicity for two chemicals, rotenone and deguelin. The objectives were to determine whether the same qAOP could be used for the three cell types, and to test chemical-independence by cross-validation with a dataset obtained on eight other chemicals in LHUMES cells. Repeating the calibration approach for both chemicals in three cell lines highlighted various practical difficulties. Even when the same readouts of KEs are measured, the mathematical functions used to describe the key event relationships may not be the same. Cross-validation in LHUMES cells was attempted by estimating chemical-specific potency at the molecular initiating events and using the rest of the calibrated qAOP to predict downstream KEs: toxicity of azoxystrobin, carboxine, mepronil and thifluzamide was underestimated. Selection of most relevant readouts and accurate characterization of the molecular initiating event for cross-validation are critical when designing in vitro experiments targeted at calibrating qAOPs., (Copyright © 2022 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Physiologically-based pharmacokinetic (PBPK) models usually include a large number of parameters whose values are obtained using in vitro to in vivo extrapolation. However, such extrapolations can be uncertain and may benefit from inclusion of evidence from clinical observations via parametric inference. When clinical interindividual variability is high, or the data sparse, it is essential to use a population pharmacokinetics inferential framework to estimate unknown or uncertain parameters. Several approaches are available for that purpose, but their relative advantages for PBPK modeling are unclear. We compare the results obtained using a minimal PBPK model of a canonical theophylline dataset with quasi-random parametric expectation maximization (QRPEM), nonparametric adaptive grid estimation (NPAG), Bayesian Metropolis-Hastings (MH), and Hamiltonian Markov Chain Monte Carlo sampling. QRPEM and NPAG gave consistent population and individual parameter estimates, mostly agreeing with Bayesian estimates. MH simulations ran faster than the others methods, which together had similar performance., (© 2022 The Authors. CPT: Pharmacometrics & Systems Pharmacology published by Wiley Periodicals LLC on behalf of American Society for Clinical Pharmacology and Therapeutics.)

Pharmacokinetics study the fate of xenobiotics in a living organism. Physiologically based pharmacokinetic (PBPK) models provide realistic descriptions of xenobiotics' absorption, distribution, metabolism, and excretion processes. They model the body as a set of homogeneous compartments representing organs, and their parameters refer to anatomical, physiological, biochemical, and physicochemical entities. They offer a quantitative mechanistic framework to understand and simulate the time-course of the concentration of a substance in various organs and body fluids. These models are well suited for performing extrapolations inherent to toxicology and pharmacology (e.g., between species or doses) and for integrating data obtained from various sources (e.g., in vitro or in vivo experiments, structure-activity models). In this chapter, we describe the practical development and basic use of a PBPK model from model building to model simulations, through implementation with an easily accessible free software., (© 2022. The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature.)

We propose a Bayesian population modeling and virtual bioequivalence assessment approach to establishing dissolution specifications for oral dosage forms. A generalizable semi-physiologically based pharmacokinetic absorption model with six gut segments and liver, connected to a two-compartment model of systemic disposition for bupropion hydrochloride oral dosage forms was developed. Prior information on model parameters for gut physiology, bupropion physicochemical properties, and drug product properties were obtained from the literature. The release of bupropion hydrochloride from immediate-, sustained- and extended-release oral dosage forms was described by a Weibull function. In vitro dissolution data were used to assign priors to the in vivo release properties of the three bupropion formulations. We applied global sensitivity analysis to identify the influential parameters for plasma bupropion concentrations and calibrated them. To quantify inter- and intra-individual variability, plasma concentration profiles in healthy volunteers that received the three dosage forms, each at two doses, were used. The calibrated model was in good agreement with both in vitro dissolution and in vivo exposure data. Markov Chain Monte Carlo samples from the joint posterior parameter distribution were used to simulate virtual crossover clinical trials for each formulation with distinct drug dissolution profiles. For each trial, an allowable range of dissolution parameters ("safe space") in which bioequivalence can be anticipated was established. These findings can be used to assure consistent product performance throughout the drug product life-cycle and to support manufacturing changes. Our framework provides a comprehensive approach to support decision-making in drug product development., (© 2021. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)