Parameterization of Physiologically Based Biopharmaceutics Models: Workshop Summary Report.

This Article shares the proceedings from the August 29th, 2023 (day 1) workshop "Physiologically Based Biopharmaceutics Modeling (PBBM) Best Practices for Drug Product Quality: Regulatory and Industry Perspectives". The focus of the day was on model parametrization; regulatory authorities from Canada, the USA, Sweden, Belgium, and Norway presented their views on PBBM case studies submitted by industry members of the IQ consortium. The presentations shared key questions raised by regulators during the mock exercise, regarding the PBBM input parameters and their justification. These presentations also shed light on the regulatory assessment processes, content, and format requirements for future PBBM regulatory submissions. In addition, the day 1 breakout presentations and discussions gave the opportunity to share best practices around key questions faced by scientists when parametrizing PBBMs. Key questions included measurement and integration of drug substance solubility for crystalline vs amorphous drugs; impact of excipients on apparent drug solubility/supersaturation; modeling of acid-base reactions at the surface of the dissolving drug; choice of dissolution methods according to the formulation and drug properties with a view to predict the in vivo performance; mechanistic modeling of in vitro product dissolution data to predict in vivo dissolution for various patient populations/species; best practices for characterization of drug precipitation from simple or complex formulations and integration of the data in PBBM; incorporation of drug permeability into PBBM for various routes of uptake and prediction of permeability along the GI tract.

Predicting Human Dermal Drug Concentrations Using PBPK Modeling and Simulation: Clobetasol Propionate Case Study.

Quantitative in silico tools may be leveraged to mechanistically predict the dermato-pharmacokinetics of compounds delivered from topical and transdermal formulations by integrating systems of rate equations that describe permeation through the formulation and layers of skin and pilo-sebaceous unit, and exchange with systemic circulation via local blood flow. Delivery of clobetasol-17 propionate (CP) from DermovateTM cream was simulated using the Transdermal Compartmental Absorption & Transit (TCATTM) Model in GastroPlus®. The cream was treated as an oil-in-water emulsion, with model input parameters estimated from publicly available information and quantitative structure-permeation relationships. From the ranges of values available for model input parameters, a set of parameters was selected by comparing model outputs to CP dermis concentration-time profiles measured by dermal open-flow microperfusion (Bodenlenz et al. Pharm Res. 33(9):2229-38, 2016). Predictions of unbound dermis CP concentrations were reasonably accurate with respect to time and skin depth. Parameter sensitivity analyses revealed considerable dependence of dermis CP concentration profiles on drug solubility in the emulsion, relatively less dependence on dispersed phase volume fraction and CP effective diffusivity in the continuous phase of the emulsion, and negligible dependence on dispersed phase droplet size. Effects of evaporative water loss from the cream and corticosteroid-induced vasoconstriction were also assessed. This work illustrates the applicability of computational modeling to predict sensitivity of dermato-pharmacokinetics to changes in thermodynamic and transport properties of a compound in a topical formulation, particularly in relation to rate-limiting steps in skin permeation. Where these properties can be related to formulation composition and processing, such a computational approach may support the design of topically applied formulations.

Bacterial metabolism rescues the inhibition of intestinal drug absorption by food and drug additives.

Food and drug products contain diverse and abundant small-molecule additives (excipients) with unclear impacts on human physiology, drug safety, and response. Here, we evaluate their potential impact on intestinal drug absorption. By screening 136 unique compounds for inhibition of the key intestinal transporter OATP2B1 we identified and validated 24 potent OATP2B1 inhibitors, characterized by higher molecular weight and hydrophobicity compared to poor or noninhibitors. OATP2B1 inhibitors were also enriched for dyes, including 8 azo (R-N=N-R') dyes. Pharmacokinetic studies in mice confirmed that FD&C Red No. 40, a common azo dye excipient and a potent inhibitor of OATP2B1, decreased the plasma level of the OATP2B1 substrate fexofenadine, suggesting that FD&C Red No. 40 has the potential to block drug absorption through OATP2B1 inhibition in vivo. However, the gut microbiomes of multiple unrelated healthy individuals as well as diverse human gut bacterial isolates were capable of inactivating the identified azo dye excipients, producing metabolites that no longer inhibit OATP2B1 transport. These results support a beneficial role for the microbiome in limiting the unintended effects of food and drug additives in the intestine and provide a framework for the data-driven selection of excipients. Furthermore, the ubiquity and genetic diversity of gut bacterial azoreductases coupled to experiments in conventionally raised and gnotobiotic mice suggest that variations in gut microbial community structure may be less important to consider relative to the high concentrations of azo dyes in food products, which have the potential to saturate gut bacterial enzymatic activity.

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

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

A Bayesian population physiologically based pharmacokinetic absorption modeling approach to support generic drug development: application to bupropion hydrochloride oral dosage forms.

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.

Interactions of Oral Molecular Excipients with Breast Cancer Resistance Protein, BCRP.

Mechanistic-understanding-based selection of excipients may improve formulation development strategies for generic drug products and potentially accelerate their approval. Our study aimed at investigating the effects of molecular excipients present in orally administered FDA-approved drug products on the intestinal efflux transporter, BCRP (ABCG2), which plays a critical role in drug absorption with potential implications on drug safety and efficacy. We determined the interactions of 136 oral molecular excipients with BCRP in isolated membrane vesicles and identified 26 excipients as BCRP inhibitors with IC50 values less than 5 µM using 3H-cholecystokinin octapeptide (3H-CCK8). These BCRP inhibitors belonged to three functional categories of excipients: dyes, surfactants, and flavoring agents. Compared with noninhibitors, BCRP inhibitors had significantly higher molecular weights and SLogP values. The inhibitory effects of excipients identified in membrane vesicles were also evaluated in BCRP-overexpressing HEK293 cells at similar concentrations. Only 1 of the 26 inhibitors of BCRP identified in vesicles inhibited BCRP-mediated 3H-oxypurinol uptake by more than 50%, consistent with the notion that BCRP inhibition depends on transmembrane or intracellular availability of the inhibitors. Collectively, the results of this study provide new information on excipient selection during the development of drug products with active pharmaceutical ingredients that are BCRP substrates.

Population pharmacokinetic analysis of AR-67, a lactone stable camptothecin analogue, in cancer patients with solid tumors.

Background AR-67 is a novel camptothecin analogue at early stages of drug development. The phase 1 clinical trial in cancer patients with solid tumors was completed and a population pharmacokinetic model (POP PK) was developed to facilitate further development of this investigational agent. Methods Pharmacokinetic data collected in the phase 1 clinical trial were utilized for the development of a population POP PK by implementing the non-linear mixed effects approach. Patient characteristics at study entry were evaluated as covariates in the model. Subjects (N = 26) were treated at nine dosage levels (1.2-12.4 mg/m2/day) on a daily × 5 schedule. Hematological toxicity data were modeled against exposure metrics. Results A two-compartment POP PK model best described the disposition of AR-67 by fitting a total of 328 PK observations from 25 subjects. Following covariate model selection, age remained as a significant covariate on central volume. The final model provided a good fit for the concentration versus time data and PK parameters were estimated with good precision. Clearance, inter-compartmental clearance, central volume and peripheral volume were estimated to be 32.2 L/h, 28.6 L/h, 6.83 L and 25.0 L, respectively. Finally, exposure-pharmacodynamic analysis using Emax models showed that plasma drug concentration versus time profiles are better predictors of AR-67-related hematologic toxicity were better predictors of leukopenia and thrombocytopenia, as compared to total dose. Conclusions A POP PK model was developed to characterize AR-67 pharmacokinetics and identified age as a significant covariate. Exposure PK metrics Cmax and AUC were shown to predict hematological toxicity. Further efforts to identify clinically relevant determinants of AR-67 disposition and effects in a larger patient population are warranted.

Leveraging Modeling and Simulation to Enhance the Efficiency of Bioequivalence Approaches for Generic Drugs: Highlights from the 2023 Generic Drug Science and Research Initiatives Public Workshop.

The 2023 Generic Drug Science and Research Initiative Public Workshop organized by the U.S. Food and Drug Administration (FDA) discussed the research needs to improve and enhance bioequivalence (BE) approaches for generic drug development. FDA takes such research needs and panel discussions into account to develop its Generic Drug User Fee Amendments III (GDUFA III) Science and Research Initiatives specific to generics. During the five workshop sessions, presentations and panel discussions focused on identifying and addressing scientific gaps and research needs related to nitrosamine impurity issues, BE assessment for oral products, innovative BE approaches for long-acting injectable products, alternative BE approaches for orally inhaled products, and advanced BE methods for topical products. Specifically, this report highlights the discussions on how to improve BE assessment for developing generic drug products based on research priorities for leveraging quantitative methods and modeling, as well as artificial intelligence/machine learning (AI/ML).

Experience Learned and Perspectives on Using Model-Integrated Evidence in the Regulatory Context for Generic Drug Products-a Meeting Report.

This report summarizes relevant insights and discussions from a 2022 FDA public workshop titled Best Practices for Utilizing Modeling Approaches to Support Generic Product Development which illustrated how model-integrated evidence has been used and can be leveraged further to inform generic drug product development and regulatory decisions during the assessment of generic drug applications submitted to the FDA. The workshop attendees discussed that model-integrated evidence (MIE) approaches for generics are being applied in the space of long-acting injectable (LAI) products to develop shorter and more cost-effective alternative study designs for LAI products. Modeling and simulation approaches are utilized to support virtual BE assessments at the site of action for locally acting drug products and to assess the impact of food on BE assessments for oral dosage forms. The factors contributing to the success of the model-informed drug development program under PDUFA VI were discussed. The generic drug industry shared that decisions on formulation candidate/formulation variant selection, on pilot in vivo bioavailability studies, and on alternative study designs for BE assessment are informed by modeling and simulation approaches. There was agreement that interactions between the regulatory agencies and the industry are desirable because they improve the industry's understanding of scientific and other regulatory considerations on implementing modeling and simulation approaches in drug development and regulatory submissions.

The effect of breast cancer resistance protein, multidrug resistant protein 1, and organic anion-transporting polypeptide 1B3 on the antitumor efficacy of the lipophilic camptothecin 7-t-butyldimethylsilyl-10-hydroxycamptothecin (AR-67) in vitro.

AR-67 (7-t-butyldimethylsilyl-10-hydroxycamptothecin) is a lipophilic camptothecin analog, currently under early stage clinical trials. Transporters are known to have an impact on the disposition of camptothecins and on the response to chemotherapeutics in general due to their expression in tumor tissues. Therefore, we investigated the interplay between the breast cancer resistance protein (BCRP), multidrug resistant protein 1 (MDR1), and organic anion-transporting polypeptide (OATP) 1B1/1B3 transporters and AR-67 and their impact on the toxicity profile of AR-67. Using cell lines expressing the aforementioned transporters, we showed that the lipophilic AR-67 lactone form is a substrate for efflux transporters BCRP and MDR1. Additionally, OATP1B1 and OATP1B3 facilitated the uptake of AR-67 carboxylate in SLCO1B1- and SLCO1B3-transfected cell systems compared with the mock-transfected ones. Notably, both BCRP and MDR1 conferred resistance to AR-67 lactone. Prompted by recent studies showing increased OATP1B3 expression in certain cancer types, we investigated the effect of OATP1B3 expression on cell viability after exposure to AR-67 carboxylate. OATP1B3-expressing cells had increased carboxylate uptake as compared with mock-transfected cells but were not sensitized because the intracellular amount of lactone was 50-fold higher than that of carboxylate and comparable between OATP1B3-expressing and OATP1B3-nonexpressing cells. In conclusion, BCRP- and MDR1-mediated efflux of AR-67 lactone confers resistance to AR-67, but OATP1B3-mediated uptake of the AR-67 carboxylate does not sensitize OATP1B3-expressing tumor cells.

Applications of Modeling and Simulation Approaches in Support of Drug Product Development of Oral Dosage Forms and Locally Acting Drug Products: a Symposium Summary.

The number of modeling and simulation applications, including physiologically based pharmacokinetic (PBPK) models, physiologically based biopharmaceutics modeling (PBBM), and empirical models, has been constantly increasing along with the regulatory acceptance of these methodologies. While aiming at minimizing unnecessary human testing, these methodologies are used today to support the development and approval of novel drug products and generics. Modeling approaches are leveraged today for assessing drug-drug interaction, informing dose adjustments in renally or hepatically impaired patients, perform dose selection in pediatrics and pregnant women and diseased populations, and conduct biopharmaceutics-related assessments such as establish clinically relevant specifications for drug products and achieve quality assurance throughout the product life cycle. In the generics space, PBPK analyses are utilized toward virtual bioequivalence assessments within the scope of alternative bioequivalence approaches, product-specific guidance development, and food effect assessments among others. Case studies highlighting the evolving and expanding role of modeling and simulation approaches within the biopharmaceutics space were presented at the symposium titled "Model Informed Drug Development (MIDD): Role in Dose Selection, Vulnerable Populations, and Biowaivers - Chemical Entities" and Prologue "PBPK/PBBM to inform the Bioequivalence Safe Space, Food Effects, and pH-mediated DDIs" at the American Association of Pharmaceutical Scientists (AAPS) PharmSci 360 Annual Meeting in Boston, MA, on October 16-19, 2022, and are summarized here.

Mechanistic modeling of drug products applied to the skin: A workshop summary report.

The development of a generic drug product involves demonstrating that there is no significant difference in the rate and extent to which the active ingredient becomes available at the site of action, relative to the reference listed drug product. This remains challenging for many locally acting topical dermatological products because measuring the concentration of the active ingredient at the site of action in the skin may not be straightforward, and, in most instances, there are no established relationships between skin and plasma pharmacokinetic profiles. In recent years, the Office of Generic Drugs of the US Food and Drug Administration (FDA) established scientific research programs with the goal of enhancing patient access to high quality, affordable topical dermatological generics. A key strategy of these research programs was to leverage modeling and simulation methodologies that accelerate the development of these generics by facilitating alternative bioequivalence approaches for dermatological drug products. This report summarizes relevant insights and discussions from a 2021 FDA public workshop titled "Regulatory Utility of Mechanistic Modeling to Support Alternative Bioequivalence Approaches," which illustrated how mechanistic modeling and simulation approaches can be utilized (and have been used) to inform generic drug product development and regulatory decisions during the assessment of generic drug applications submitted to the FDA.

Regulatory utility of mechanistic modeling to support alternative bioequivalence approaches: A workshop overview.

On September 30 and October 1, 2021, the US Food and Drug Administration (FDA) and the Center for Research on Complex Generics cosponsored a live virtual workshop titled "Regulatory Utility of Mechanistic Modeling to Support Alternative Bioequivalence Approaches." The overall aims of the workshop included (i) engaging the generic drug industry and other involved stakeholders regarding how mechanistic modeling and simulation can support their product development and regulatory submissions; (ii) sharing the current state of mechanistic modeling for bioequivalence (BE) assessment through case studies; (iii) establishing a consensus on best practices for using mechanistic modeling approaches, such as physiologically based pharmacokinetic modeling and computational fluid dynamics modeling, for BE assessment; and (iv) introducing the concept of a Model Master File to improve model sharing between model developers, industry, and the FDA. More than 1500 people registered for the workshop. Based on a postworkshop survey, the majority of participants reported that their fundamental scientific understanding of mechanistic models was enhanced, there was greater consensus on model validation and verification, and regulatory expectations for mechanistic modeling submitted in abbreviated new drug applications were clarified by the workshop.

Drug release patterns and cytotoxicity of PEG-poly(aspartate) block copolymer micelles in cancer cells.

PURPOSE: To test physicochemical and biological properties of PEG-poly(aspartate) [PEG-p(Asp)] block copolymer micelles entrapping doxorubicin hydrochloride (DOX) through ionic interaction. METHODS: PEG-p(Asp) was synthesized from 5 kDa PEG and 20 Asp units. Carboxyl groups of p(Asp) were present as benzyl ester [PEG-p(Asp/Bz)], sodium salt [PEG-p(Asp/Na)] or free acid [PEG-p(Asp/H)]. Block copolymers and DOX were mixed at various ratios to prepare polymer micelles, which were subsequently characterized to determine particle size, drug loading and release patterns, and cytotoxicity against prostate (PC3 and DU145) and lung (A549) cancer cell lines. RESULTS: PEG-p(Asp/Bz), Na- and H-micelles entrapped 1.1, 56.8 and 40.6 wt.% of DOX, respectively. Na- and H-micelles (<100 nm) showed time-dependent DOX release at pH 7.4, which was accelerated at pH 5.0. Na-micelles were most stable at pH 7.4, retaining 31.8% of initial DOX for 48 h. Cytotoxicity of Na-micelles was 23.2% (A549), 28.5% (PC3) and 45.9% (DU145) more effective than free DOX. CONCLUSION: Ionic interaction appeared to entrap DOX efficiently in polymer micelles from PEG-p(Asp) block copolymers. Polymer micelles possessing counter ions (Na) of DOX in the core were the most stable, releasing drugs for prolonged time in a pH-dependent manner, and suppressing cancer cells effectively.

Physiologically-Based Pharmacokinetic Modeling to Support Determination of Bioequivalence for Dermatological Drug Products: Scientific and Regulatory Considerations.

Physiologically-based pharmacokinetic (PBPK) modeling and simulation provides mechanism-based predictions of the pharmacokinetics of an active ingredient following its administration in humans. Dermal PBPK models describe the skin permeation and disposition of the active ingredient following the application of a dermatological product on the skin of virtual healthy and diseased human subjects. These models take into account information on product quality attributes, physicochemical properties of the active ingredient and skin (patho)physiology, and their interplay with each other. Regulatory and product development decision makers can leverage these quantitative tools to identify factors impacting local and systemic exposure. In the realm of generic drug products, the number of US Food and Drug Administratioin (FDA) interactions that use dermal PBPK modeling to support alternative bioequivalence (BE) approaches is increasing. In this report, we share scientific considerations on the development, verification and validation (V&V), and application of PBPK models within the context of a virtual BE assessment for dermatological drug products. We discuss the challenges associated with model V&V for these drug products stemming from the fact that target-site active ingredient concentrations are typically not measurable. Additionally, there are no established relationships between local and systemic PK profiles, when the latter are quantifiable. To that end, we detail a multilevel model V&V approach involving validation for the model of the drug product of interest coupled with the overall assessment of the modeling platform in use while leveraging in vitro and in vivo data related to local and systemic bioavailability.

A Critical Overview of the Biological Effects of Excipients (Part I): Impact on Gastrointestinal Absorption.

Traditionally, excipients have been considered in drug development from the perspective of their influence on drug solubility, manufacturability, and ability to control in vitro and in vivo drug release. These effects have been largely evaluated through studies involving in vitro dissolution methods. However, there is a growing awareness that what had previously been considered biologically inert excipients can exert numerous in vivo effects. This includes the potential to change gastrointestinal (GI) transit time, enterocyte passive transcellular or paracellular permeability, active transport activity, or presystemic drug metabolism. In this critical overview of the biological effects of excipients (Part I), we provide a summary of select published studies that explore these various in vivo factors. We also include a table that points readers to published reviews that list a range of excipients known to have biological activity. A subsequent discussion on in vitro, in vivo, and in silico methods that can be used to explore these excipient effects is provided in a separate (Part 2) continuation of this critical overview.

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

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

Percutaneous absorption of steroids from finite doses: Predicting urinary excretion from in vitro skin permeation testing.

Pharmacokinetic (PK) models are widely used to describe drug permeation across the epidermal membrane barrier, the stratum corneum (SC). Here, we extend our previously reported diffusion and compartment-in-series models to describe plasma concentrations, urinary excretion-time profiles and exposure estimates after topically applied finite doses of solvent deposited solids. In vivo models were derived by convolution of a skin absorption input function for finite dosing with that for in vivo disposition PK. In vitro skin permeation test (IVPT) and in vivo urinary excretion data for cortisone, desoxycorticosterone, and testosterone were extracted from literature for model validation and establishment of in vitro - in vivo relationships (IVIVR). Both SC diffusion and SC 3-compartment-in-series PK models adequately described experimental in vitro and in vivo permeation data, with similar model parameter estimates for SC diffusion time and bioavailability. A satisfactory IVIVR was generated for cortisone, whereas testosterone and desoxycorticosterone showed higher bioavailability in vitro compared to in vivo. In recognising that future prospective studies need to both have an adequate sampling schedule and be harmonized for robust IVIVRs, we developed expressions for predicting extent of absorption and time for peak absorption for both in vitro and in vivo studies. Other study parameters, such as application site, applied dose, and application techniques, can also affect drug permeability through skin during dosage form metamorphosis after finite dose application, and a lack of correlation may result if these are poorly matched.

Physiologically-based pharmacokinetic modeling to support bioequivalence and approval of generic products: A case for diclofenac sodium topical gel, 1.

Establishing bioequivalence (BE) for dermatological drug products by conducting comparative clinical end point studies can be costly and the studies may not be sufficiently sensitive to detect certain formulation differences. Quantitative methods and modeling, such as physiologically-based pharmacokinetic (PBPK) modeling, can support alternative BE approaches with reduced or no human testing. To enable PBPK modeling for regulatory decision making, models should be sufficiently verified and validated (V&V) for the intended purpose. This report illustrates the US Food and Drug Administration (FDA) approval of a generic diclofenac sodium topical gel that was based on a totality of evidence, including qualitative and quantitative sameness and physical and structural similarity to the reference product, an in vivo BE study with PK end points, and, more importantly, for the purposes of this report, a virtual BE assessment leveraging dermal PBPK modeling and simulation instead of a comparative clinical end point study in patients. The modeling approach characterized the relationship between systemic (plasma) and local (skin and synovial fluid) diclofenac exposure and demonstrated BE between the generic and reference products at the presumed site of action. Based on the fit-for-purpose modeling principle, the V&V process involved assessing observed data of diclofenac concentrations in skin tissues and plasma, and the overall performance of the modeling platform for relevant products. Using this case as an example, this report provides current scientific considerations on good practices for model V&V and the establishment of BE for dermatological drug products when leveraging PBPK modeling and simulation for regulatory decision making.