Use of Model-Informed Drug Development to Streamline Development of Long-Acting Products: Can These Successes Be Translated to Long-Acting Hormonal Contraceptives?

Long-acting contraceptives are the most effective reversible contraceptive methods. Increasing patients' access to these contraceptives may translate into fewer unintended pregnancies and lead to substantial individual and public health benefits. However, development of long-acting products can be complex and challenging. This review provides (a) an overview of representative development programs for long-acting antipsychotics as cases for conceptual translation to long-acting contraceptives, (b) several case examples on how modeling and simulation have been used to streamline the development of long-acting products, and (c) examples of challenges andopportunities in developing long-acting contraceptives and information on how exposure-response relationships of commonly used progestins may enable regulators and developers to rely on prior findings of effectiveness and safety from an approved contraceptive to streamline the development of long-acting contraceptives. The US Food and Drug Administration is seeking assistance from stakeholders to provide data from studies in which pharmacokinetic and pharmacodynamic or clinical outcomes of hormonal contraceptives were evaluated and not previously submitted.

Asymmetric Hydroboration Approach to the Scalable Synthesis of ((1R,3S)-1-Amino-3-((R)-6-hexyl-5,6,7,8-tetrahydronaphthalen-2-yl)cyclopentyl)methanol (BMS-986104) as a Potent S1P1 Receptor Modulator.

We describe a highly efficient route for the synthesis of 4a (BMS-986104). A key step in the synthesis is the asymmetric hydroboration of trisubstituted alkene 6. Particularly given the known difficulties involved in this type of transformation (6 → 7), the current methodology provides an efficient approach to prepare this class of compounds. In addition, we disclose the efficacy of 4a in a mouse EAE model, which is comparable to 4c (FTY720). Mechanistically, 4a exhibited excellent remyelinating effects on lysophosphatidylcholine (LPC) induced demyelination in a three-dimensional brain cell culture assay.

Identification of Tricyclic Agonists of Sphingosine-1-phosphate Receptor 1 (S1P1) Employing Ligand-Based Drug Design.

Fingolimod (1) is the first approved oral therapy for the treatment of relapsing remitting multiple sclerosis. While the phosphorylated metabolite of fingolimod was found to be a nonselective S1P receptor agonist, agonism specifically of S1P1 is responsible for the peripheral blood lymphopenia believed to be key to its efficacy. Identification of modulators that maintain activity on S1P1 while sparing activity on other S1P receptors could offer equivalent efficacy with reduced liabilities. We disclose in this paper a ligand-based drug design approach that led to the discovery of a series of potent tricyclic agonists of S1P1 with selectivity over S1P3 and were efficacious in a pharmacodynamic model of suppression of circulating lymphocytes. Compound 10 had the desired pharmacokinetic (PK) and pharmacodynamic (PD) profile and demonstrated maximal efficacy when administered orally in a rat adjuvant arthritis model.

Discovery and Structure-Activity Relationship (SAR) of a Series of Ethanolamine-Based Direct-Acting Agonists of Sphingosine-1-phosphate (S1P1).

Sphingosine-1-phosphate (S1P) is a bioactive sphingolipid metabolite that regulates a multitude of physiological processes such as lymphocyte trafficking, cardiac function, vascular development, and inflammation. Because of the ability of S1P1 receptor agonists to suppress lymphocyte egress, they have great potential as therapeutic agents in a variety of autoimmune diseases. In this article, the discovery of selective, direct acting S1P1 agonists utilizing an ethanolamine scaffold containing a terminal carboxylic acid is described. Potent S1P1 agonists such as compounds 18a and 19a which have greater than 1000-fold selectivity over S1P3 are described. These compounds efficiently reduce blood lymphocyte counts in rats through 24 h after single doses of 1 and 0.3 mpk, respectively. Pharmacodynamic properties of both compounds are discussed. Compound 19a was further studied in two preclinical models of disease, exhibiting good efficacy in both the rat adjuvant arthritis model (AA) and the mouse experimental autoimmune encephalomyelitis model (EAE).

Identification and synthesis of potent and selective pyridyl-isoxazole based agonists of sphingosine-1-phosphate 1 (S1P1).

The synthesis and structure-activity relationship (SAR) of a series of pyridyl-isoxazole based agonists of S1P1 are discussed. Compound 5b provided potent in vitro activity with selectivity, had an acceptable pharmacokinetic profile, and demonstrated efficacy in a dose dependent manner when administered orally in a rodent model of arthritis.

Potent and Selective Agonists of Sphingosine 1-Phosphate 1 (S1P1): Discovery and SAR of a Novel Isoxazole Based Series.

Sphingosine 1-phosphate (S1P) is the endogenous ligand for the sphingosine 1-phosphate receptors (S1P1-5) and evokes a variety of cellular responses through their stimulation. The interaction of S1P with the S1P receptors plays a fundamental physiological role in a number of processes including vascular development and stabilization, lymphocyte migration, and proliferation. Agonism of S1P1, in particular, has been shown to play a significant role in lymphocyte trafficking from the thymus and secondary lymphoid organs, resulting in immunosuppression. This article will detail the discovery and SAR of a potent and selective series of isoxazole based full agonists of S1P1. Isoxazole 6d demonstrated impressive efficacy when administered orally in a rat model of arthritis and in a mouse experimental autoimmune encephalomyelitis (EAE) model of multiple sclerosis.

Identification and Preclinical Pharmacology of BMS-986104: A Differentiated S1P1 Receptor Modulator in Clinical Trials.

Clinical validation of S1P receptor modulation therapy was achieved with the approval of fingolimod (Gilenya, 1) as the first oral therapy for relapsing remitting multiple sclerosis. However, 1 causes a dose-dependent reduction in the heart rate (bradycardia), which occurs within hours after first dose. We disclose the identification of clinical compound BMS-986104 (3d), a novel S1P1 receptor modulator, which demonstrates ligand-biased signaling and differentiates from 1 in terms of cardiovascular and pulmonary safety based on preclinical pharmacology while showing equivalent efficacy in a T-cell transfer colitis model.

Development of novel, 384-well high-throughput assay panels for human drug transporters: drug interaction and safety assessment in support of discovery research.

Transporter proteins are known to play a critical role in affecting the overall absorption, distribution, metabolism, and excretion characteristics of drug candidates. In addition to efflux transporters (P-gp, BCRP, MRP2, etc.) that limit absorption, there has been a renewed interest in influx transporters at the renal (OATs, OCTs) and hepatic (OATPs, BSEP, NTCP, etc.) organ level that can cause significant clinical drug-drug interactions (DDIs). Several of these transporters are also critical for hepatobiliary disposition of bilirubin and bile acid/salts, and their inhibition is directly implicated in hepatic toxicities. Regulatory agencies took action to address transporter-mediated DDI with the goal of ensuring drug safety in the clinic and on the market. To meet regulatory requirements, advanced bioassay technology and automation solutions were implemented for high-throughput transporter screening to provide structure-activity relationship within lead optimization. To enhance capacity, several functional assay formats were miniaturized to 384-well throughput including novel fluorescence-based uptake and efflux inhibition assays using high-content image analysis as well as cell-based radioactive uptake and vesicle-based efflux inhibition assays. This high-throughput capability enabled a paradigm shift from studying transporter-related issues in the development space to identifying and dialing out these concerns early on in discovery for enhanced mechanism-based efficacy while circumventing DDIs and transporter toxicities.

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.

Characterization of efflux transporters involved in distribution and disposition of apixaban.

The studies reported here were conducted to investigate the transport characteristics of apixaban (1-(4-methoxyphenyl)-7-oxo-6-(4-(2-oxopiperidin-1-yl)phenyl)-4,5,6,7-tetrahydro-1H-pyrazolo[3,4-c]pyridine-3-carboxamide) and to understand the impact of transporters on apixaban distribution and disposition. In human permeability glycoprotein (P-gp)- and breast cancer resistance protein (BCRP)-cDNA-transfected cell monolayers as well as Caco-2 cell monolayers, the apparent efflux ratio of basolateral-to-apical (PcB-A) versus apical-to-basolateral permeability (PcA-B) of apixaban was >10. The P-gp- and BCRP-facilitated transport of apixaban was concentration- and time-dependent and did not show saturation over a wide range of concentrations (1-100 µM). The efflux transport of apixaban was also demonstrated by the lower mucosal-to-serosal permeability than that of the serosal-to-mucosal direction in isolated rat jejunum segments. Apixaban did not inhibit digoxin transport in Caco-2 cells. Ketoconazole decreased the P-gp-mediated apixaban efflux in Caco-2 and the P-gp-cDNA-transfected cell monolayers, but did not affect the apixaban efflux to a meaningful extent in the BCRP-cDNA-transfected cell monolayers. Coincubation of a P-gp inhibitor (ketoconazole or cyclosporin A) and a BCRP inhibitor (Ko134) provided more complete inhibition of apixaban efflux in Caco-2 cells than separate inhibition by individual inhibitors. Naproxen inhibited apixaban efflux in Caco-2 cells but showed only a minimal effect on apixaban transport in the BCRP-transfected cells. Naproxen was the first nonsteroidal antiinflammatory drug that was demonstrated as a weak P-gp inhibitor. These results demonstrate that apixaban is a substrate for efflux transporters P-gp and BCRP, which can help explain its low brain penetration, and low fetal exposures and high milk excretion in rats.

Cynomolgus monkey as a potential model to assess drug interactions involving hepatic organic anion transporting polypeptides: in vitro, in vivo, and in vitro-to-in vivo extrapolation.

Organic anion-transporting polypeptides (OATP) 1B1, 1B3, and 2B1 can serve as the loci of drug-drug interactions (DDIs). In the present work, the cynomolgus monkey was evaluated as a potential model for studying OATP-mediated DDIs. Three cynomolgus monkey OATPs (cOATPs), with a high degree of amino acid sequence identity (91.9, 93.5, and 96.6% for OATP1B1, 1B3, and 2B1, respectively) to their human counterparts, were cloned, expressed, and characterized. The cOATPs were stably transfected in human embryonic kidney cells and were functionally similar to the corresponding human OATPs (hOATPs), as evident from the similar uptake rate of typical substrates (estradiol-17ß-d-glucuronide, cholecystokinin octapeptide, and estrone-3-sulfate). Moreover, six known hOATP inhibitors exhibited similar IC(50) values against cOATPs. To further evaluate the appropriateness of the cynomolgus monkey as a model, a known hOATP substrate [rosuvastatin (RSV)]-inhibitor [rifampicin (RIF)] pair was examined in vitro; the monkey-derived parameters (RSV K(m) and RIF IC(50)) were similar (within 3.5-fold) to those obtained with hOATPs and human primary hepatocytes. In vivo, the area under the plasma concentration-time curve of RSV (3 mg/kg, oral) given 1 hour after a single RIF dose (15 mg/kg, oral) was increased 2.9-fold in cynomolgus monkeys, consistent with the value (3.0-fold) reported in humans. A number of in vitro-in vivo extrapolation approaches, considering the fraction of the pathways affected and free versus total inhibitor concentrations, were also explored. It is concluded that the cynomolgus monkey has the potential to serve as a useful model for the assessment of OATP-mediated DDIs in a nonclinical setting.

Generation and characterization of a novel multidrug resistance protein 2 humanized mouse line.

The multidrug resistance protein (MRP) 2 is predominantly expressed in liver, intestine, and kidney, where it plays an important role in the excretion of a range of drugs and their metabolites or endogenous compounds into bile, feces, and urine. Mrp knockout [Mrp2(-/-)] mice have been used recently to study the role of MRP2 in drug disposition. Here, we describe the first generation and initial characterization of a mouse line humanized for MRP2 (huMRP2), which is nulled for the mouse Mrp2 gene and expresses the human transporter in the organs and cell types where MRP2 is normally expressed. Analysis of the mRNA expression for selected cytochrome P450 and transporter genes revealed no major changes in huMRP2 mice compared with wild-type controls. We show that human MRP2 is able to compensate functionally for the loss of the mouse transporter as demonstrated by comparable bilirubin levels in the humanized mice and wild-type controls, in contrast to the hyperbilirubinemia phenotype that is observed in MRP2(-/-) mice. The huMRP2 mouse provides a model to study the role of the human transporter in drug disposition and in assessing the in vivo consequences of inhibiting this transporter by compounds interacting with human MRP2.

Validation of in vitro cell-based human blood-brain barrier model using clinical positron emission tomography radioligands to predict in vivo human brain penetration.

We have evaluated a novel in vitro cell-based human blood-brain barrier (BBB) model that could predict in vivo human brain penetration for compounds with different BBB permeabilities using the clinical positron emission tomography (PET) data. Comparison studies were also performed to demonstrate that the in vitro cell-based human BBB model resulted in better predictivity over the traditional permeability model in discovery organizations, Caco-2 cells. We evaluated the in vivo BBB permeability of [(18)F] and [(11)C]-compounds in humans by PET imaging. The in vivo plasma-brain exchange parameters used for comparison were determined in humans by PET using a kinetic analysis of the radiotracer binding. For each radiotracer, the parameters were determined by fitting the brain kinetics of the radiotracer using a two-tissue compartment model of the ligand-receptor interaction. Bidirectional transport studies with the same compounds as in in vivo studies were carried out using the in vitro cell-based human BBB model as well as Caco-2 cells. The in vitro cell-based human BBB model has important features of the BBB in vivo and is suitable for discriminating between CNS and non-CNS marketed drugs. A very good correlation (r(2) = 0.90; P < 0.001) was demonstrated between in vitro BBB permeability and in vivo permeability coefficient. In contrast, a poor correlation (r(2) = 0.17) was obtained between Caco-2 data and in vivo human brain penetration. This study highlights the potential of this in vitro cell-based human BBB model in drug discovery and shows that it can be an extremely effective screening tool for CNS programs.

Small molecule antagonist of leukocyte function associated antigen-1 (LFA-1): structure-activity relationships leading to the identification of 6-((5S,9R)-9-(4-cyanophenyl)-3-(3,5-dichlorophenyl)-1-methyl-2,4-dioxo-1,3,7-triazaspiro[4.4]nonan-7-yl)nicotinic acid (BMS-688521).

Leukocyte function-associated antigen-1 (LFA-1), also known as CD11a/CD18 or alpha(L)beta(2), belongs to the beta(2) integrin subfamily and is constitutively expressed on all leukocytes. The major ligands of LFA-1 include three intercellular adhesion molecules 1, 2, and 3 (ICAM 1, 2, and 3). The interactions between LFA-1 and the ICAMs are critical for cell adhesion, and preclinical animal studies and clinical data from the humanized anti-LFA-1 antibody efalizumab have provided proof-of-concept for LFA-1 as an immunological target. This article will detail the structure-activity relationships (SAR) leading to a novel second generation series of highly potent spirocyclic hydantoin antagonists of LFA-1. With significantly enhanced in vitro and ex vivo potency relative to our first clinical compound (1), as well as demonstrated in vivo activity and an acceptable pharmacokinetic and safety profile, 6-((5S,9R)-9-(4-cyanophenyl)-3-(3,5-dichlorophenyl)-1-methyl-2,4-dioxo-1,3,7-triazaspiro-[4.4]nonan-7-yl)nicotinic acid (2e) was selected to advance into clinical trials.

SAR of PXR transactivation in benzimidazole-based IGF-1R kinase inhibitors.

The SAR of PXR transactivation by 3-(benzimidazol-2-yl)-pyridine-2-one based ATP competitive inhibitors of Insulin-like Growth Factor 1 Receptor kinase (IGF-1R) is discussed. Compounds without PXR transactivation, with in vivo antitumor activity, reduced protein binding and improved oral exposure are presented.

Discovery of a 2,4-disubstituted pyrrolo[1,2-f][1,2,4]triazine inhibitor (BMS-754807) of insulin-like growth factor receptor (IGF-1R) kinase in clinical development.

This report describes the biological activity, characterization, and SAR leading to 9d (BMS-754807) a small molecule IGF-1R kinase inhibitor in clinical development.

P-gp inhibition potential in cell-based models: which "calculation" method is the most accurate?

The objective was to directly compare the four different "calculation" methods of assessing P-gp inhibition potential using experimental data obtained from approximately 60 structurally diverse internal research and marketed compounds. Bidirectional studies for digoxin (probe for P-gp substrate) were performed with and without test compounds (at 10 microM). Four different calculation methods were applied to the same dataset (raw bidirectional permeability values) to obtain the "percent inhibition of P-gp" for these compounds using the different methods. Significantly different inhibition potential was obtained with the "exact" same experimental dataset depending on the calculation method used. Subsequently, entirely different conclusions regarding the "inhibition potential" of test compound was reached due to the different calculation methods. Based on the direct comparison of these methods, method no. 3 (i.e., inhibition of B to A permeability of digoxin) is recommended as the calculation method ideal during screening stages due to its high throughput amenability. The methodology is capable of rapidly screening compounds with adequate reliability for early stage drug discovery. Method no. 3 provides an abridged version of a bidirectional study that is fully capable of identifying all non-inhibitors (0-20%), moderate inhibitors (20-60%), and potent inhibitors (>60%) and demonstrates high correlation with method no. 1 (inhibition based on both A to B and B to A permeability of digoxin). Nevertheless, method no. 1 might be appropriate for more detailed mechanistic studies required in late stage discovery and development.

Discovery and evaluation of 4-(2-(4-chloro-1H-pyrazol-1-yl)ethylamino)-3-(6-(1-(3-fluoropropyl)piperidin-4-yl)-4-methyl-1H-benzo[d]imidazol-2-yl)pyridin-2(1H)-one (BMS-695735), an orally efficacious inhibitor of insulin-like growth factor-1 receptor kinase with broad spectrum in vivo antitumor activity.

We previously reported that 1 (BMS-536924), a benzimidazole inhibitor of the insulin-like growth factor-1 receptor, had demonstrated in vivo antitumor activity. This lead compound was found to have potent CYP3A4 inhibition, CYP3A4 induction mediated by PXR transactivation, poor aqueous solubility, and high plasma protein binding. Herein we disclose the evolution of this chemotype to address these issues. This effort led to 10 (BMS-695735), which exhibits improved ADME properties, a low risk for drug-drug interactions, and in vivo efficacy in multiple xenograft models.