Intrinsic lipolysis rate for systematic design of lipid-based formulations.

Lipid-based formulations (LBFs) are used by the pharmaceutical industry in oral delivery systems for both poorly water-soluble drugs and biologics. Digestibility is key for the performance of LBFs and in vitro lipolysis is commonly used to compare the digestibility of LBFs. Results from in vitro lipolysis experiments depend highly on the experimental conditions and formulation characteristics, such as droplet size (which defines the surface area available for digestion) and interfacial structure. This study introduced the intrinsic lipolysis rate (ILR) as a surface area-independent approach to compare lipid digestibility. Pure acylglycerol nanoemulsions, stabilized with polysorbate 80 at low concentration, were formulated and digested according to a standardized pH-stat lipolysis protocol. A methodology originally developed to calculate the intrinsic dissolution rate of poorly water-soluble drugs was adapted for the rapid calculation of ILR from lipolysis data. The impact of surfactant concentration on the apparent lipolysis rate and lipid structure on ILR was systematically investigated. The surfactant polysorbate 80 inhibited lipolysis of tricaprylin nanoemulsions in a concentration-dependent manner. Coarse-grained molecular dynamics simulations supported these experimental observations. In the absence of bile and phospholipids, tricaprylin was shielded from lipase at 0.25% polysorbate 80. In contrast, the inclusion of bile salt and phospholipid increased the surfactant-free area and improved the colloidal presentation of the lipids to the enzyme, especially at 0.125% polysorbate 80. At a constant and low surfactant content, acylglycerol digestibility increased with decreasing acyl chain length, decreased esterification, and increasing unsaturation. The calculated ILR of pure acylglycerols was successfully used to accurately predict the IRL of binary lipid mixtures. The ILR measurements hold great promise as an efficient method supporting pharmaceutical formulation scientists in the design of LBFs with specific digestion profiles.

A Tribute to Professor Per Artursson - Scientist, Explorer, Mentor, Innovator, and Giant in Pharmaceutical Research.

This issue of the Journal of Pharmaceutical Sciences is dedicated to Professor Per Artursson and the groundbreaking contributions he has made and continues to make in the Pharmaceutical Sciences. Per is one of the most cited researchers in his field, with more than 30,000 citations and an h-index of 95 as of September 2020. Importantly, these citations are distributed over the numerous fields he has explored, clearly showing the high impact the research has had on the discipline. We provide a short portrait of Per, with emphasis on his personality, driving forces and the inspirational sources that shaped his career as a world-leading scientist in the field. He is a curious scientist who deftly moves between disciplines and has continued to innovate, expand boundaries, and profoundly impact the pharmaceutical sciences throughout his career. He has developed new tools and provided insights that have significantly contributed to today's molecular and mechanistic approaches to research in the fields of intestinal absorption, cellular disposition, and exposure-efficacy relationships of pharmaceutical drugs. We want to celebrate these important contributions in this special issue of the Journal of Pharmaceutical Sciences in Per's honor.

Input Estimation for Extended-Release Formulations Exemplified with Exenatide.

Estimating the in vivo absorption profile of a drug is essential when developing extended-release medications. Such estimates can be obtained by measuring plasma concentrations over time and inferring the absorption from a model of the drug's pharmacokinetics. Of particular interest is to predict the bioavailability-the fraction of the drug that is absorbed and enters the systemic circulation. This paper presents a framework for addressing this class of estimation problems and gives advice on the choice of method. In parametric methods, a model is constructed for the absorption process, which can be difficult when the absorption has a complicated profile. Here, we place emphasis on non-parametric methods that avoid making strong assumptions about the absorption. A modern estimation method that can address very general input-estimation problems has previously been presented. In this method, the absorption profile is modeled as a stochastic process, which is estimated using Markov chain Monte Carlo techniques. The applicability of this method for extended-release formulation development is evaluated by analyzing a dataset of Bydureon, an injectable extended-release suspension formulation of exenatide, a GLP-1 receptor agonist for treating diabetes. This drug is known to have non-linear pharmacokinetics. Its plasma concentration profile exhibits multiple peaks, something that can make parametric modeling challenging, but poses no major difficulties for non-parametric methods. The method is also validated on synthetic data, exploring the effects of sampling and noise on the accuracy of the estimates.

Tenapanor administration and the activity of the H+ -coupled transporter PepT1 in healthy volunteers.

AIM: Tenapanor (RDX5791/AZD1722), an inhibitor of gastrointestinal Na+ /H+ exchanger NHE3, is being evaluated for the treatment of patients with constipation-predominant irritable bowel syndrome and the treatment of hyperphosphataemia in patients with chronic kidney disease on dialysis. By reducing intestinal H+ secretion, inhibition of NHE3 by tenapanor could indirectly affect H+ -coupled transporter activity, leading to drug-drug interactions. We investigated the effect of tenapanor on the activity of the H+ -coupled peptide transporter PepT1 via assessment of the pharmacokinetics of cefadroxil - a compound transported by PepT1 - in healthy volunteers. METHODS: In this open-label, two-period crossover, phase 1 study (NCT02140281), 28 volunteers received in random order: a single dose of cefadroxil 500 mg for 1 day; and tenapanor 15 mg twice daily over 4 days followed by single doses of both cefadroxil 500 mg and tenapanor 15 mg on day 5. There was a 4-day washout between treatment periods. RESULTS: Cefadroxil exposure was similar when administered alone or in combination with tenapanor {geometric least-squares mean ratios [(cefadroxil + tenapanor)/cefadroxil] (90% confidence interval): area under the concentration-time curve 93.3 (90.6-96.0)%; maximum concentration in plasma 95.9 (89.8-103)%}. Tenapanor treatment caused a softening of stool consistency and an increase in stool frequency, consistent with its expected pharmacodynamic effect. No safety concerns were identified and tenapanor was not detected in plasma. CONCLUSIONS: These results suggest that tenapanor 15 mg twice daily does not have a clinically relevant impact on the activity of the H+ -coupled transporter PepT1 in humans. This may guide future research on drug-drug interactions involving NHE3 inhibitors.

Substrate and method dependent inhibition of three ABC-transporters (MDR1, BCRP, and MRP2).

Drug transport and drug-drug interactions (DDI) with human ABC transporters are generally investigated in mammalian cell lines or inverted membrane vesicles from insect cells (Sf9) overexpressing the transporter of interest. In this study, we instead used membrane vesicles from human embryonic kidney cells (HEK293) overexpressing wild type MDR1/Pgp (ABCB1), BCRP (ABCG2), and MRP2 (ABCC2) with the aim to study the concentration dependent inhibition of shared and prototypic probe substrates. We first investigated 15 substrates and identified estrone-17-beta-glucorinide (E17G) as shared substrate. Nine specific and general inhibitors were then studied using E17G and prototypic probe substrates. The results were compared with those previously obtained in Sf9 vesicles and cell lines of canine (MDCKII) and human (Saos-2) origin. For the majority of inhibitors, Ki values differed <10-fold between E17G and probe substrates. Significant differences in Ki values were observed for about one third of the inhibitors. The transport inhibition potencies in HEK293 vesicles were in good agreement with those obtained in Sf9 vesicles. Large differences were found in the inhibition potencies observed in the vesicular systems compared to the cellular systems. Nevertheless, the rank order correlations between the different experimental systems were generally good. Our study provides further information on substrate dependent inhibition of ABC-transporters, and suggests that simple ranking of compounds can be used as a tier one approach to bridge results obtained in different experimental systems.

Combined in vitro-in vivo approach to assess the hepatobiliary disposition of a novel oral thrombin inhibitor.

Two clinical trials and a large set of in vitro transporter experiments were performed to investigate if the hepatobiliary disposition of the direct thrombin inhibitor prodrug AZD0837 is the mechanism for the drug-drug interaction with ketoconazole observed in a previous clinical study. In Study 1, [(3)H]AZD0837 was administered to healthy male volunteers (n = 8) to quantify and identify the metabolites excreted in bile. Bile was sampled directly from the jejunum by duodenal aspiration via an oro-enteric tube. In Study 2, the effect of ketoconazole on the plasma and bile pharmacokinetics of AZD0837, the intermediate metabolite (AR-H069927), and the active form (AR-H067637) was investigated (n = 17). Co-administration with ketoconazole elevated the plasma exposure to AZD0837 and the active form approximately 2-fold compared to placebo, which may be explained by inhibited CYP3A4 metabolism and reduced biliary clearance, respectively. High concentrations of the active form was measured in bile with a bile-to-plasma AUC ratio of approximately 75, indicating involvement of transporter-mediated excretion of the compound. AZD0837 and its metabolites were further investigated as substrates of hepatic uptake and efflux transporters in vitro. Studies in MDCK-MDR1 cell monolayers and P-glycoprotein (P-gp) expressing membrane vesicles identified AZD0837, the intermediate, and the active form as substrates of P-gp. The active form was also identified as a substrate of the multidrug and toxin extrusion 1 (MATE1) transporter and the organic cation transporter 1 (OCT1), in HEK cells transfected with the respective transporter. Ketoconazole was shown to inhibit all of these three transporters; in particular, inhibition of P-gp and MATE1 occurred in a clinically relevant concentration range. In conclusion, the hepatobiliary transport pathways of AZD0837 and its metabolites were identified in vitro and in vivo. Inhibition of the canalicular transporters P-gp and MATE1 may lead to enhanced plasma exposure to the active form, which could, at least in part, explain the clinical interaction with ketoconazole.

A P-gp vesicular transport inhibition assay - optimization and validation for drug-drug interaction testing.

Accurate determination of potential drug-drug interaction mediated by efflux transporters (tDDI) is crucial to assess the risk of pharmacokinetic interaction and toxicity of drugs. Passive permeability and uptake transporter mediated transport are important covariates of cell-based inhibition assays that need to be taken into consideration when performing kinetic analysis of data. Vesicular uptake inhibition has been considered by regulatory agencies as a viable alternative for testing tDDI potential of low passive permeability drugs in particular. Membranes prepared from a P-gp overexpressing human cell line has superior transport properties over membranes prepared from Sf9 cells and cholesterol enriched Sf9 membranes. P-gp expressed in this membrane effluxes N-methyl-quinidine (NMQ) with high affinity (K(m) is 3.65 µM) and a high rate (V(max) is 656 pmol/mg protein/min). Digoxin, vinblastine and paclitaxel, established P-gp substrates inhibited transport of NMQ with estimated K(i) values of 250, 0.1 and 0.6 µM, respectively. A panel of 11 drugs that have been listed by regulatory agencies as reference inhibitors were used to validate the assay to predict clinical inhibition potential. All the drugs that have been implicated in P-gp mediated DDI were found to be inhibitors in a relevant concentration range.

Variability in P-glycoprotein inhibitory potency (IC50) using various in vitro experimental systems: implications for universal digoxin drug-drug interaction risk assessment decision criteria.

A P-glycoprotein (P-gp) IC50 working group was established with 23 participating pharmaceutical and contract research laboratories and one academic institution to assess interlaboratory variability in P-gp IC50 determinations. Each laboratory followed its in-house protocol to determine in vitro IC50 values for 16 inhibitors using four different test systems: human colon adenocarcinoma cells (Caco-2; eleven laboratories), Madin-Darby canine kidney cells transfected with MDR1 cDNA (MDCKII-MDR1; six laboratories), and Lilly Laboratories Cells--Porcine Kidney Nr. 1 cells transfected with MDR1 cDNA (LLC-PK1-MDR1; four laboratories), and membrane vesicles containing human P-glycoprotein (P-gp; five laboratories). For cell models, various equations to calculate remaining transport activity (e.g., efflux ratio, unidirectional flux, net-secretory-flux) were also evaluated. The difference in IC50 values for each of the inhibitors across all test systems and equations ranged from a minimum of 20- and 24-fold between lowest and highest IC50 values for sertraline and isradipine, to a maximum of 407- and 796-fold for telmisartan and verapamil, respectively. For telmisartan and verapamil, variability was greatly influenced by data from one laboratory in each case. Excluding these two data sets brings the range in IC50 values for telmisartan and verapamil down to 69- and 159-fold. The efflux ratio-based equation generally resulted in severalfold lower IC50 values compared with unidirectional or net-secretory-flux equations. Statistical analysis indicated that variability in IC50 values was mainly due to interlaboratory variability, rather than an implicit systematic difference between test systems. Potential reasons for variability are discussed and the simplest, most robust experimental design for P-gp IC50 determination proposed. The impact of these findings on drug-drug interaction risk assessment is discussed in the companion article (Ellens et al., 2013) and recommendations are provided.

Application of receiver operating characteristic analysis to refine the prediction of potential digoxin drug interactions.

In the 2012 Food and Drug Administration (FDA) draft guidance on drug-drug interactions (DDIs), a new molecular entity that inhibits P-glycoprotein (P-gp) may need a clinical DDI study with a P-gp substrate such as digoxin when the maximum concentration of inhibitor at steady state divided by IC50 ([I1]/IC50) is ≥0.1 or concentration of inhibitor based on highest approved dose dissolved in 250 ml divide by IC50 ([I2]/IC50) is ≥10. In this article, refined criteria are presented, determined by receiver operating characteristic analysis, using IC50 values generated by 23 laboratories. P-gp probe substrates were digoxin for polarized cell-lines and N-methyl quinidine or vinblastine for P-gp overexpressed vesicles. Inhibition of probe substrate transport was evaluated using 15 known P-gp inhibitors. Importantly, the criteria derived in this article take into account variability in IC50 values. Moreover, they are statistically derived based on the highest degree of accuracy in predicting true positive and true negative digoxin DDI results. The refined criteria of [I1]/IC50 ≥ 0.03 and [I2]/IC50 ≥ 45 and FDA criteria were applied to a test set of 101 in vitro-in vivo digoxin DDI pairs collated from the literature. The number of false negatives (none predicted but DDI observed) were similar, 10 and 12%, whereas the number of false positives (DDI predicted but not observed) substantially decreased from 51 to 40%, relative to the FDA criteria. On the basis of estimated overall variability in IC50 values, a theoretical 95% confidence interval calculation was developed for single laboratory IC50 values, translating into a range of [I1]/IC50 and [I2]/IC50 values. The extent by which this range falls above the criteria is a measure of risk associated with the decision, attributable to variability in IC50 values.

In vitro and in silico strategies to identify OATP1B1 inhibitors and predict clinical drug-drug interactions.

PURPOSE: To establish in vitro and in silico models that predict clinical drug-drug interactions (DDIs) with the OATP1B1 (SLCO1B1) transporter. METHODS: The inhibitory effect of 146 drugs and drug-like compounds on OATP1B1-mediated transport was studied in HEK293 cells. A computational model was developed to predict OATP1B1 inhibition. Concentration-dependent effects were investigated for six compounds; clinical DDIs were predicted by calculating change in exposure (i.e. R-values) in eight different ways. RESULTS: Sixty-five compounds were identified as OATP1B1 inhibitors at 20 µM. The computational model predicted the test set with 80% accuracy for inhibitors and 91% for non-inhibitors. In vitro-in vivo comparisons underscored the importance of using drugs with known clinical effects as references. Thus, reference drugs, cyclosporin A, gemfibrozil, and fenofibrate, provided an inhibition interval to which three antiviral drugs, atazanavir, lopinavir, and amprenavir, could be compared and their clinical DDIs with OATP1B1 classified. CONCLUSIONS: Twenty-two new OATP1B1 inhibitors were identified, a predictive OATP1B1 inhibition in silico model was developed, and successful predictions of clinical DDIs were obtained with OATP1B1.

Development of a highly sensitive method using liquid chromatography-multiple reaction monitoring to quantify membrane P-glycoprotein in biological matrices and relationship to transport function.

The quantification of P-glycoprotein [P-gp, ABCB1, multidrug resistance 1 (MDR1)] protein in biological matrices is considered a key factor missing for useful translation of in vitro functional data to the in vivo situation and for comparison of transporter data among different in vitro models. In the present study a liquid chromatography (LC)-mass spectrometry method was developed to quantify P-gp membrane protein levels in different biological matrices. The amount of P-gp transporter protein was measured in Caco-2 cell monolayers and in inside-out human embryonic kidney (HEK)-MDR1 vesicles. From both in vitro systems, two preparations with different functionality were used. Transporter function was determined as digoxin efflux in Caco-2 cell monolayers and N-methylquinidine (NMQ) uptake in membrane vesicles, and, in addition, mRNA expression in the Caco-2 monolayers was measured. The results showed an excellent relationship between NMQ uptake functionality in inside-out HEK-MDR1 vesicles and protein contents. Similar concordance between the digoxin efflux and P-gp content in different Caco-2 cell cultures was observed, whereas mRNA levels are indicative of increased P-gp content and activity in older Caco-2 cultures, however, not yielding the same quantitative relationship. The results from both Caco-2 and HEK-MDR1 membrane vesicles confirm that the protein content is directly related to the level of activity in the respective system. The method presented here to quantify P-gp protein by LC-multiple reaction monitoring will facilitate the development of future methodologies to bridge between expression systems and cell/tissue models and to scale from in vitro models to whole organs.

Effects of ketoconazole on the in vivo biotransformation and hepatobiliary transport of the thrombin inhibitor AZD0837 in pigs.

Ketoconazole has been shown in clinical trials to increase the plasma exposure of the oral anticoagulant prodrug AZD0837 [(2S)-N-{4- [(Z)-amino(methoxyimino)methyl]benzyl}-1-{(2R)-2-[3-chloro-5-(difluoromethoxy)phenyl]-2-hydroxyethanoyl}-azetidine-2-carboxamide] and its active metabolite, AR-H067637 [(2S)-N-{4-[amino(imino)methyl]benzyl}-1-{(2R)-2-[3-chloro-5-(difluoromethoxy)phenyl]-2-hydroxyethanoyl}-azetidine-2-carboxamide]. To investigate the biotransformation of AZD0837 and the effect of ketoconazole on this process, we used an experimental model in pigs that allows repeated sampling from three blood vessels, the bile duct, and a perfused intestinal segment. The pigs received AZD0837 (500 mg) given enterally either alone (n = 5) or together with single-dose ketoconazole (600 mg) (n = 6). The prodrug (n = 2) and its active metabolite (n = 2) were also administered intravenously to provide reference doses. The plasma data revealed considerable interindividual variation in the exposure of the prodrug, intermediate metabolite, and active metabolite. However, AR-H067637 was detected at very high concentrations in the bile with low variability (Ae(bile) = 53 ± 6% of the enteral dose), showing that the compound had indeed been formed in all of the animals and efficiently transported into the bile canaliculi. Concomitant dosing with ketoconazole increased the area under the plasma concentration-time curve for AZD0837 (by 99%) and for AR-H067637 (by 51%). The effect on the prodrug most likely arose from inhibited CYP3A-mediated metabolism. Reduced metabolism also seemed to explain the increased plasma exposure of the active compound because ketoconazole prolonged the terminal half-life with no apparent effect on the extensive biliary excretion and biliary clearance. These in vivo results were supported by in vitro depletion experiments for AR-H067637 in pig liver microsomes with and without the addition of ketoconazole.