Acid reducing agents (ARAs) reduce the dissolution rate of weakly basic drugs in the stomach potentially leading to lower bioavailability. Formulating the API as a rapidly dissolving salt is one strategy employed to reduce the impact of ARAs on dissolution of such drugs. In the present work, a model drug was selected with an immediate release formulation of the free base dosed in both the absence and presence of the ARA famotidine. In the latter case, bioavailability is restricted and several salt formulations were investigated. To simulate these drug products a mechanistic physiologically based pharmacokinetic (PBPK) model was built using the Simcyp Simulator, which illustrates the advantage of formulating an API as a salt compared to the free base form. The simulations use a mechanistic salt model utilising knowledge of the solubility product which was applied to predict the salt advantage. The developed PBPK model exemplifies that it can be critical to account for the surface pH and solubility when modelling the dissolution of low pKa bases and their salts in the gastric environment. In particular, the mechanistic salt model can be used to aid in screening and salt form selection where the aim is to mitigate effects of ARAs.

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. Study Highlights WHAT IS THE CURRENT KNOWLEDGE ON THE TOPIC?No single software platform allows practitioners to compare parametric and nonparametric estimates of calibrated physiologically-based pharmacokinetic (PBPK) model parameters. WHAT QUESTION DID THIS STUDY ADDRESS?How do a maximum likelihood parametric (quasi-random parametric expectation maximization) and nonparametric (nonparametric adaptive grid estimation) algorithms, and two Bayesian numerical methods (Hamiltonian Markov Chain Monte Carlo and MH) compare in results and timing for calibration of a PBPK model?

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, the similarity, or otherwise, of profiles from two different formulation in the gastrointestinal tract (GIT) may not be associated with same conclusions based on systemic circulation. Here, we have explored the concordance of BE conclusions (or otherwise) for a set of hypothetical formulations based on drug concentration-profiles in various segments of gut vs those in systemic circulation using virtual trials powered by physiologically-based (PBPK) models. Moreover, by changing physiological and biological attributes of the GIT, the impact of Crohn’s Disease on the BE conclusions was explored. The virtual BE (VBE) was applied to an Entocort EC, a pH/time-dependent release formulation of budesonide, as model drug that is targeted at ileum and colon. Substantial “discordance” between the 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 with systemic-derived VBE. Nonetheless, a higher f2 value than the default similarity standard (50) was required to achieve local BE. In the case of Crohn’s Disease, a product-specific f2 value might be needed for drugs such as Entocort EC to ensure local BE, knowing the insensitivity of systemic data to show lack of BE. The results, being limited to the single model drug product, Entocort EC, may not provide a a general and conclusive picture. However, for the first time, it demonstrated the potential discordances between the local gut vs systemic BE.