Applications of Quantitative Systems Pharmacology in Model-Informed Drug Discovery: Perspective on Impact and Opportunities.

Quantitative systems pharmacology (QSP) approaches have been increasingly applied in the pharmaceutical since the landmark white paper published in 2011 by a National Institutes of Health working group brought attention to the discipline. In this perspective, we discuss QSP in the context of other modeling approaches and highlight the impact of QSP across various stages of drug development and therapeutic areas. We discuss challenges to the field as well as future opportunities.

Mathematical and Experimental Validation of Flux Dialysis Method: An Improved Approach to Measure Unbound Fraction for Compounds with High Protein Binding and Other Challenging Properties.

A flux dialysis method to measure unbound fraction (fu) of compounds with high protein binding and other challenging properties was tested and validated. This method is based on the principle that the initial flux rate of a compound through a size-excluding dialysis membrane is proportional to the product of the compound initial concentration, fu, and unbound dialysis membrane permeability (Pmem). Therefore, fu can be determined from the initial concentration and flux rate, assuming membrane Pmem is known. Compound initial flux rates for 14 compounds were determined by dialyzing human plasma containing compound (donor side) versus compound-free plasma (receiver side) and measuring the rate of compound appearance into the receiver side. Eleven compounds had known fu values obtained from conventional methods (ranging from 0.000013 to 0.22); three compounds (bedaquiline, lapatinib, and pibrentasvir) had previously qualified fu values (e.g., <0.001).Pmem estimated from flux rates and known fu values did not meaningfully differ among the compounds and were consistent with previously published values, indicating that Pmem is a constant for the dialysis membrane. This Pmem constant and the individual compound flux rates were used to calculate fu values. The flux dialysis fu values for the 11 compounds were in good agreement with their reported fu values (all within 2.5-fold; R2 = 0.980), confirming the validity of the method. Furthermore, the flux dialysis method allowed discrete fu to be estimated for the three compounds with previously qualified fu Theoretical and experimental advantages of the flux dialysis method over other dialysis-based protein binding methods are discussed.

Preclinical QSP Modeling in the Pharmaceutical Industry: An IQ Consortium Survey Examining the Current Landscape.

A cross-industry survey was conducted to assess the landscape of preclinical quantitative systems pharmacology (QSP) modeling within pharmaceutical companies. This article presents the survey results, which provide insights on the current state of preclinical QSP modeling in addition to future opportunities. Our results call attention to the need for an aligned definition and consistent terminology around QSP, yet highlight the broad applicability and benefits preclinical QSP modeling is currently delivering.

Translational pharmacokinetic-pharmacodynamic analysis in the pharmaceutical industry: an IQ Consortium PK-PD Discussion Group perspective.

With inadequate efficacy being the primary cause for the attrition of drug candidates in clinical development, the need to better predict clinical efficacy earlier in the drug development process has increased in importance in the pharmaceutical industry. Here, we review current applications of translational pharmacokinetic-pharmacodynamic (PK-PD) modeling of preclinical data in the pharmaceutical industry, including best practices. Preclinical translational PK-PD modeling has been used in many therapeutic areas and has been impactful to drug development. The role of preclinical translational PK-PD modeling in drug discovery and development will continue to evolve and broaden, given that its broad implementation in the pharmaceutical industry is relatively recent and many opportunities still exist for its further application.

Mechanisms and Predictions of Drug-Drug Interactions of the Hepatitis C Virus Three Direct-Acting Antiviral Regimen: Paritaprevir/Ritonavir, Ombitasvir, and Dasabuvir.

To assess drug-drug interaction (DDI) potential for the three direct-acting antiviral (3D) regimen of ombitasvir, dasabuvir, and paritaprevir, in vitro studies profiled drug-metabolizing enzyme and transporter interactions. Using mechanistic static and dynamic models, DDI potential was predicted for CYP3A, CYP2C8, UDP-glucuronosyltransferase (UGT) 1A1, organic anion-transporting polypeptide (OATP) 1B1/1B3, breast cancer resistance protein (BCRP), and P-glycoprotein (P-gp). Perpetrator static model DDI predictions for metabolizing enzymes were within 2-fold of the clinical observations, but additional physiologically based pharmacokinetic modeling was necessary to achieve the same for drug transporters. When perpetrator interactions were assessed, ritonavir was responsible for the strong increase in exposure of sensitive CYP3A substrates, whereas paritaprevir (an OATP1B1/1B3 inhibitor) greatly increased the exposure of sensitive OATP1B1/1B3 substrates. The 3D regimen drugs are UGT1A1 inhibitors and are predicted to moderately increase plasma exposure of sensitive UGT1A1 substrates. Paritaprevir, ritonavir, and dasabuvir are BCRP inhibitors. Victim DDI predictions were qualitatively in line with the clinical observations. Plasma exposures of the 3D regimen were reduced by strong CYP3A inducers (paritaprevir and ritonavir; major CYP3A substrates) but were not affected by strong CYP3A4 inhibitors, since ritonavir (a CYP3A inhibitor) is already present in the regimen. Strong CYP2C8 inhibitors increased plasma exposure of dasabuvir (a major CYP2C8 substrate), OATP1B1/1B3 inhibitors increased plasma exposure of paritaprevir (an OATP1B1/1B3 substrate), and P-gp or BCRP inhibitors (all compounds are substrates of P-gp and/or BCRP) increased plasma exposure of the 3D regimen. Overall, the comprehensive mechanistic assessment of compound disposition along with mechanistic and PBPK approaches to predict victim and perpetrator DDI liability may enable better clinical management of nonstudied drug combinations with the 3D regimen.

5-sulfonyl-benzimidazoles as selective CB2 agonists-part 2.

In a previous communication, the SAR of a series of potent and selective 5-sulfonyl-benzimidazole CB2-receptor agonists was described. The lack of in vivo activity of compounds from this series was attributed to their poor solubility and metabolic stability. In this Letter, we report on the further optimization of this series, leading to the relatively polar and peripherically acting CB2 agonists 41 and 49. Although both compounds were not active in acute pain models, the less selective compound 41 displayed good, sustained activity in a chronic model of neuropathic pain without the tolerance observed with morphine. In addition, both 41 and 49 delayed the onset of clinical symptoms in an experimental model for Multiple sclerosis.

Cardio-vascular safety beyond hERG: in silico modelling of a guinea pig right atrium assay.

As chemists can easily produce large numbers of new potential drug candidates, there is growing demand for high capacity models that can help in driving the chemistry towards efficacious and safe candidates before progressing towards more complex models. Traditionally, the cardiovascular (CV) safety domain plays an important role in this process, as many preclinical CV biomarkers seem to have high prognostic value for the clinical outcome. Throughout the industry, traditional ion channel binding data are generated to drive the early selection process. Although this assay can generate data at high capacity, it has the disadvantage of producing high numbers of false negatives. Therefore, our company applies the isolated guinea pig right atrium (GPRA) assay early-on in discovery. This functional multi-channel/multi-receptor model seems much more predictive in identifying potential CV liabilities. Unfortunately however, its capacity is limited, and there is no room for full automation. We assessed the correlation between ion channel binding and the GPRA's Rate of Contraction (RC), Contractile Force (CF), and effective refractory frequency (ERF) measures assay using over six thousand different data points. Furthermore, the existing experimental knowledge base was used to develop a set of in silico classification models attempting to mimic the GPRA inhibitory activity. The Naïve Bayesian classifier was used to built several models, using the ion channel binding data or in silico computed properties and structural fingerprints as descriptors. The models were validated on an independent and diverse test set of 200 reference compounds. Performances were assessed on the bases of their overall accuracy, sensitivity and specificity in detecting both active and inactive molecules. Our data show that all in silico models are highly predictive of actual GPRA data, at a level equivalent or superior to the ion channel binding assays. Furthermore, the models were interpreted in terms of the descriptors used to highlight the undesirable areas in the explored chemical space, specifically regions of low polarity, high lipophilicity and high molecular weight. In conclusion, we developed a predictive in silico model of a complex physiological assay based on a large and high quality set of experimental data. This model allows high throughput in silico safety screening based on chemical structure within a given chemical space.

5-Sulfonyl-benzimidazoles as selective CB2 agonists.

A novel series of benzimidazole CB2-receptor agonists was synthesized and the structure-activity relationship explored. The results showed agonistic activities with an EC(50) up to 0.5 nM and excellent selectivity (>4000-fold) over the CB1 receptor. The size of the substituent on the 2-position determined the level of agonism, ranging from inverse agonism to partial agonism to full agonism, which was more pronounced for the rat CB2 receptor. A wide variation of sulfonyl substituents at the benzimidazole 5-position was tolerated, which was used to optimize the drug-like properties. This resulted into lead compound 14j that can be used to investigate the potential of a selective, peripherically acting CB2 agonist. The in vitro profile of key compounds is displayed using pie bar charts (VlaaiVis).

Prediction of human pharmacokinetics using physiologically based modeling: a retrospective analysis of 26 clinically tested drugs.

The aim of this study was to evaluate different physiologically based modeling strategies for the prediction of human pharmacokinetics. Plasma profiles after intravenous and oral dosing were simulated for 26 clinically tested drugs. Two mechanism-based predictions of human tissue-to-plasma partitioning (P(tp)) from physicochemical input (method Vd1) were evaluated for their ability to describe human volume of distribution at steady state (V(ss)). This method was compared with a strategy that combined predicted and experimentally determined in vivo rat P(tp) data (method Vd2). Best V(ss) predictions were obtained using method Vd2, providing that rat P(tp) input was corrected for interspecies differences in plasma protein binding (84% within 2-fold). V(ss) predictions from physicochemical input alone were poor (32% within 2-fold). Total body clearance (CL) was predicted as the sum of scaled rat renal clearance and hepatic clearance projected from in vitro metabolism data. Best CL predictions were obtained by disregarding both blood and microsomal or hepatocyte binding (method CL2, 74% within 2-fold), whereas strong bias was seen using both blood and microsomal or hepatocyte binding (method CL1, 53% within 2-fold). The physiologically based pharmacokinetics (PBPK) model, which combined methods Vd2 and CL2 yielded the most accurate predictions of in vivo terminal half-life (69% within 2-fold). The Gastroplus advanced compartmental absorption and transit model was used to construct an absorption-disposition model and provided accurate predictions of area under the plasma concentration-time profile, oral apparent volume of distribution, and maximum plasma concentration after oral dosing, with 74%, 70%, and 65% within 2-fold, respectively. This evaluation demonstrates that PBPK models can lead to reasonable predictions of human pharmacokinetics.

The prediction of drug metabolism, tissue distribution, and bioavailability of 50 structurally diverse compounds in rat using mechanism-based absorption, distribution, and metabolism prediction tools.

The aim of this study was to assess a physiologically based modeling approach for predicting drug metabolism, tissue distribution, and bioavailability in rat for a structurally diverse set of neutral and moderate-to-strong basic compounds (n = 50). Hepatic blood clearance (CL(h)) was projected using microsomal data and shown to be well predicted, irrespective of the type of hepatic extraction model (80% within 2-fold). Best predictions of CL(h) were obtained disregarding both plasma and microsomal protein binding, whereas strong bias was seen using either blood binding only or both plasma and microsomal protein binding. Two mechanistic tissue composition-based equations were evaluated for predicting volume of distribution (V(dss)) and tissue-to-plasma partitioning (P(tp)). A first approach, which accounted for ionic interactions with acidic phospholipids, resulted in accurate predictions of V(dss) (80% within 2-fold). In contrast, a second approach, which disregarded ionic interactions, was a poor predictor of V(dss) (60% within 2-fold). The first approach also yielded accurate predictions of P(tp) in muscle, heart, and kidney (80% within 3-fold), whereas in lung, liver, and brain, predictions ranged from 47% to 62% within 3-fold. Using the second approach, P(tp) prediction accuracy in muscle, heart, and kidney was on average 70% within 3-fold, and ranged from 24% to 54% in all other tissues. Combining all methods for predicting V(dss) and CL(h) resulted in accurate predictions of the in vivo half-life (70% within 2-fold). Oral bioavailability was well predicted using CL(h) data and Gastroplus Software (80% within 2-fold). These results illustrate that physiologically based prediction tools can provide accurate predictions of rat pharmacokinetics.