Experimental Error, Kurtosis, Activity Cliffs, and Methodology: What Limits the Predictivity of Quantitative Structure-Activity Relationship Models?

Given a particular descriptor/method combination, some quantitative structure-activity relationship (QSAR) datasets are very predictive by random-split cross-validation while others are not. Recent literature in modelability suggests that the limiting issue for predictivity is in the data, not the QSAR methodology, and the limits are due to activity cliffs. Here, we investigate, on in-house data, the relative usefulness of experimental error, distribution of the activities, and activity cliff metrics in determining how predictive a dataset is likely to be. We include unmodified in-house datasets, datasets that should be perfectly predictive based only on the chemical structure, datasets where the distribution of activities is manipulated, and datasets that include a known amount of added noise. We find that activity cliff metrics determine predictivity better than the other metrics we investigated, whatever the type of dataset, consistent with the modelability literature. However, such metrics cannot distinguish real activity cliffs due to large uncertainties in the activities. We also show that a number of modern QSAR methods, and some alternative descriptors, are equally bad at predicting the activities of compounds on activity cliffs, consistent with the assumptions behind "modelability." Finally, we relate time-split predictivity with random-split predictivity and show that different coverages of chemical space are at least as important as uncertainty in activity and/or activity cliffs in limiting predictivity.

Cutaneous Adverse Events Caused by Sulfonamide-Containing Drugs: Reality or Perception?

The term "sulfa allergy", originally coined for the sulfonamide antibiotics, has become associated with any drugs that contain a sulfonamide moiety. This raises a question: should medicinal chemists avoid making sulfonamide-containing compounds in drug discovery programs? The negative perception of sulfonamides is not supported by any systematic study or data analysis. To address this gap, an analysis of postmarketing reports of cutaneous AEs for drugs with and without a sulfonamide group was conducted. The analysis revealed no evidence of association between the presence of a sulfonamide moiety and a high reporting rate of cutaneous AEs and indicated that the risk of such AEs was associated with the presence of certain structural alerts and higher daily doses. These results strongly suggest that sulfonamide-containing compounds are not at a higher risk of presenting with a cutaneous adverse drug reaction compared with non-sulfonamides and that medicinal chemists should not avoid use of the sulfonamide group.

Moving beyond Binary Predictions of Human Drug-Induced Liver Injury (DILI) toward Contrasting Relative Risk Potential.

The hepatic risk matrix (HRM) was developed and used to differentiate lead clinical and back-up drug candidates against competitor/marketed drugs within the same pharmaceutical class for their potential to cause human drug-induced liver injury (DILI). The hybrid HRM scoring system blends physicochemical properties (Rule of Two Model: dose and lipophilicity or Partition Model: dose, ionization state, lipophilicity, and fractional carbon bond saturation) with common toxicity mechanisms (cytotoxicity, mitochondrial dysfunction, and bile salt export pump (BSEP) inhibition) that promote DILI. HRM scores are based on bracketed safety margins (<1, 1-10, 10-100, and >100× clinical Cmax,total). On the basis of well-established clinical safety experience of marketed/withdrawn drug candidates, the background analysis consists of 200 drugs from the Liver Toxicity Knowledge Base annotated as Most-DILI- (79), Less-DILI- (56), No-DILI- (47), and Ambiguous-DILI-concern (18) drugs. Scores were generated for over 21 internal and 7 external drug candidates discontinued for unacceptable incidence/magnitude of liver transaminase elevations during clinical trials or withdrawn for liver injury severity. Both hybrid scoring systems identified 70-80% Most-DILI-concern drugs, but more importantly, stratified successful/unsuccessful drug candidates for liver safety (incidence/severity of transaminase elevations and approved drug labels). Incorporating other mechanisms (reactive metabolite and cytotoxic metabolite generation and hepatic efflux transport inhibition, other than BSEP) to the HRM had minimal beneficial impact in DILI prediction/stratification. As is, the hybrid scoring system was positioned for portfolio assessments to contrast DILI risk potential of small molecule drug candidates in early clinical development. This stratified approach for DILI prediction aided decisions regarding drug candidate progression, follow-up mechanistic work, back-up selection, clinical dose selection, and due diligence assessments in favor of compounds with less implied clinical hepatotoxicity risk.

The Identification of Pivotal Transcriptional Factors Mediating Cell Responses to Drugs With Drug-Induced Liver Injury Liabilities.

Drug-induced liver injury (DILI) is a leading cause of drug attrition during drug development and a common reason for drug withdrawal from the market. The poor predictability of conventional animal-based approaches necessitates the development of alternative testing approaches. A body of evidence associates DILI with the induction of stress-response genes in liver cells. Here, we set out to identify signal transduction pathways predominantly involved in the regulation of gene transcription by DILI drugs. To this end, we employed ATTAGENE's cell-based multiplexed reporter assay, the FACTORIAL transcription factor (TF), that enables quantitative assessment of the activity of multiple stress-responsive TFs in a single well of cells. Homogeneous reporter system enables quantitative functional assessment of multiple transcription factors. Nat. Methods 5, 253-260). Using this assay, we assessed TF responses of the human hepatoma cell line HepG2 to a panel of 64 drug candidates, including 23 preclinical DILI and 11 clinical DILI compounds and 30 nonhepatotoxic compounds from a diverse physicochemical property space. We have identified 16 TF families that specifically responded to DILI drugs, including nuclear factor (erythroid-derived 2)-like 2 antioxidant response element, octamer, hypoxia inducible factor 1 alpha, farnesoid-X receptor, TCF/beta-catenin, aryl hydrocarbon receptor, activator protein-1, E2F, early growth response-1, metal-response transcription factor 1, sterol regulatory element-binding protein, paired box protein, peroxisome proliferator-activated receptor, liver X receptor, interferone regulating factor, and P53, and 2 promoters that responded to multiple TFs (cytomegalovirus and direct repeat 3/vitamin D receptor). Some of TFs identified here also have previously defined role in pathogenesis of liver diseases. These data demonstrate the utility of cost-effective, animal-free, TF profiling assay for detecting DILI potential of drug candidates at early stages of drug development.

Mechanisms of Skin Toxicity Associated with Metabotropic Glutamate Receptor 5 Negative Allosteric Modulators.

Cutaneous reactions represent one of the most common adverse drug effects observed in clinical trials leading to substantial compound attrition. Three negative allosteric modulators (NAMs) of metabotropic glutamate receptors (mGluRs), which represent an important target for neurological diseases, developed by Pfizer, were recently failed in preclinical development due to delayed type IV skin hypersensitivity observed in non-human primates (NHPs). Here we employed large-scale phenotypic profiling in standardized panels of human primary cell/co-culture systems to characterize the skin toxicity mechanism(s) of mGluR5 NAMs from two different series. Investigation of a database of chemicals tested in these systems and transcriptional profiling suggested that the mechanism of toxicity may involve modulation of nuclear receptor targets RAR/RXR, and/or VDR with AhR antagonism. The studies reported here demonstrate how phenotypic profiling of preclinical drug candidates using human primary cells can provide insights into the mechanisms of toxicity and inform early drug discovery and development campaigns.

Evaluating the Role of Multidrug Resistance Protein 3 (MDR3) Inhibition in Predicting Drug-Induced Liver Injury Using 125 Pharmaceuticals.

The role of bile salt export protein (BSEP) inhibition in drug-induced liver injury (DILI) has been investigated widely, while inhibition of the canalicular multidrug resistant protein 3 (MDR3) has received less attention. This transporter plays a pivotal role in secretion of phospholipids into bile and functions coordinately with BSEP to mediate the formation of bile acid-containing biliary micelles. Therefore, inhibition of MDR3 in human hepatocytes was examined across 125 drugs (70 of Most-DILI-concern and 55 of No-DILI-concern). Of these tested, 41% of Most-DILI-concern and 47% of No-DILI-concern drugs had MDR3 IC50 values of <50 µM. A better distinction across DILI classifications occurred when systemic exposure was considered where safety margins of 50-fold had low sensitivity (0.29), but high specificity (0.96). Analysis of physical chemical property space showed that basic compounds were twice as likely to be MDR3 inhibitors as acids, neutrals, and zwitterions and that inhibitors were more likely to have polar surface area (PSA) values of <100 Å2 and cPFLogD values between 1.5 and 5. These descriptors, with different cutoffs, also highlighted a group of compounds that shared dual potency as MDR3 and BSEP inhibitors. Nine drugs classified as Most-DILI-concern compounds (four withdrawn, four boxed warning, and one liver injury warning in their approved label) had intrinsic potency features of <20 µM in both assays, thereby reinforcing the notion that multiple inhibitory mechanisms governing bile formation (bile acid and phospholipid efflux) may confer additional risk factors that play into more severe forms of DILI as shown by others for BSEP inhibitors combined with multidrug resistance-associated protein (MRP2, MRP3, MRP4) inhibitory properties. Avoiding physical property descriptors that highlight dual BSEP and MDR3 inhibition or testing drug candidates for inhibition of multiple efflux transporters (e.g., BSEP, MDR3, and MRPs) may be an effective strategy for prioritizing drug candidates with less likelihood of causing clinical DILI.

Amine promiscuity and toxicology analysis.

Drug discovery programs often face challenges to obtain sufficient duration of action of the drug (i.e. seek longer half-lives). If the pharmacodynamic response is driven by free plasma concentration of the drug then extending the plasma drug concentration is a valid approach. Half-life is dependent on the volume of distribution, which in turn can be dependent upon the ionization state of the molecule. Basic compounds tend to have a higher volume of distribution leading to longer half-lives. However, it has been shown that bases may also have higher promiscuity. In this work, we describe an analysis of in vitro pharmacological profiling and toxicology data investigating the role of primary, secondary, and tertiary amines in imparting promiscuity and thus off-target toxicity. Primary amines are found to be less promiscuous in in vitro assays and have improved profiles in in vivo toxicology studies compared to secondary and tertiary amines.

Setting Clinical Exposure Levels of Concern for Drug-Induced Liver Injury (DILI) Using Mechanistic in vitro Assays.

Severe drug-induced liver injury (DILI) remains a major safety issue due to its frequency of occurrence, idiosyncratic nature, poor prognosis, and diverse underlying mechanisms. Numerous experimental approaches have been published to improve human DILI prediction with modest success. A retrospective analysis of 125 drugs (70 = most-DILI, 55 = no-DILI) from the Food and Drug Administration Liver Toxicity Knowledge Base was used to investigate DILI prediction based on consideration of human exposure alone or in combination with mechanistic assays of hepatotoxic liabilities (cytotoxicity, bile salt export pump inhibition, or mitochondrial inhibition/uncoupling). Using this dataset, human plasma Cmax,total ≥ 1.1 µM alone distinguished most-DILI from no-DILI compounds with high sensitivity/specificity (80/73%). Accounting for human exposure improved the sensitivity/specificity for each assay and helped to derive predictive safety margins. Compounds with plasma Cmax,total ≥ 1.1 µM and triple liabilities had significantly higher odds ratio for DILI than those with single/dual liabilities. Using this approach, a subset of recent pharmaceuticals with evidence of liver injury during clinical development was recognized as potential hepatotoxicants. In summary, plasma Cmax,total ≥ 1.1 µM along with multiple mechanistic liabilities is a major driver for predictions of human DILI potential. In applying this approach during drug development the challenge will be generating accurate estimates of plasma Cmax,total at efficacious doses in advance of generating true exposure data from clinical studies. In the meantime, drug candidates with multiple hepatotoxic liabilities should be deprioritized, since they have the highest likelihood of causing DILI in case their efficacious plasma Cmax,total in humans is higher than anticipated.

Chemotypes sensitivity and predictivity of in vivo outcomes for cytotoxic assays in THLE and HepG2 cell lines.

In the present study, we demonstrate the utility of in vitro ATP depletion assays in both THLE and HepG2 cells for predicting the toxicological outcome in Exploratory Toxicology Studies across 446 Pfizer proprietary compounds. Our results suggest a higher likelihood of selecting suitable compounds for in vivo safety studies by using cytotoxicity assays in multiple cell-lines over a single cell line. In addition, we demonstrate that different cell-lines have different sensitivities to compounds depending on their ionization state, that is, acid, base or neutral. HepG2 cells are more sensitive for basic compounds, whereas THLE cells have a relatively higher sensitivity for the acidic and neutral compounds. These in vitro cytotoxicity assays when combined with physicochemical properties (cLogP >3 and topological polar surface area (TPSA) <75Å(2)), are the most effective means to prioritize compounds having a lower probability of causing adverse events in vivo.

Synthesis, biological evaluation, hydration site thermodynamics, and chemical reactivity analysis of α-keto substituted peptidomimetics for the inhibition of Plasmodium falciparum.

A new series of peptidomimetic pseudo-prolyl-homophenylalanylketones were designed, synthesized and evaluated for inhibition of the Plasmodium falciparum cysteine proteases falcipain-2 (FP-2) and falcipain-3 (FP-3). In addition, the parasite killing activity of these compounds in human blood-cultured P. falciparum was examined. Of twenty-two (22) compounds synthesized, one peptidomimetic comprising a homophenylalanine-based α-hydroxyketone linked Cbz-protected hydroxyproline (39) showed the most potency (IC50 80 nM against FP-2 and 60 nM against FP-3). In silico analysis of these peptidomimetic analogs offered important protein-ligand structural insights including the role, by WaterMap, of water molecules in the active sites of these protease isoforms. The pseudo-dipeptide 39 and related compounds may serve as a promising direction forward in the design of competitive inhibitors of falcipains for the effective treatment of malaria.

Finding the rules for successful drug optimisation.

Drug discovery is a process of multiparameter optimisation, with the objective of finding compounds that achieve multiple, project-specific property criteria. These criteria are often based on the subjective opinion of the project team, but analysis of historical data can help to find the most appropriate profile. Computational 'rule induction' approaches enable an objective analysis of complex data to identify interpretable, multiparameter rules that distinguish compounds with the greatest likelihood of success for a project. Each property criterion highlights the most critical data that enable effective compound prioritisation decisions. We illustrate this with two applications: determining rules for simple, drug-like properties; and exploring experimental target inhibition data to find rules to reduce the risk of toxicity.

Analysis of Pfizer compounds in EPA's ToxCast chemicals-assay space.

The U.S. Environmental Protection Agency (EPA) launched the ToxCast program in 2007 with the goal of evaluating high-throughput in vitro assays to prioritize chemicals that need toxicity testing. Their goal was to develop predictive bioactivity signatures for toxic compounds using a set of in vitro assays and/or in silico properties. In 2009, Pfizer joined the ToxCast initiative by contributing 52 compounds with preclinical and clinical data for profiling across the multiple assay platforms available. Here, we describe the initial analysis of the Pfizer subset of compounds within the ToxCast chemical (n = 1814) and in vitro assay (n = 486) space. An analysis of the hit rate of Pfizer compounds in the ToxCast assay panel allowed us to focus our mining of assays potentially most relevant to the attrition of our compounds. We compared the bioactivity profile of Pfizer compounds to other compounds in the ToxCast chemical space to gain insights into common toxicity pathways. Additionally, we explored the similarity in the chemical and biological spaces between drug-like compounds and environmental chemicals in ToxCast and compared the in vivo profiles of a subset of failed pharmaceuticals having high similarity in both spaces. We found differences in the chemical and biological spaces of pharmaceuticals compared to environmental chemicals, which may question the applicability of bioactivity signatures developed exclusively based on the latter to drug-like compounds if used without prior validation with the ToxCast Phase-II chemicals. Finally, our analysis has allowed us to identify novel interactions for our compounds in particular with multiple nuclear receptors that were previously not known. This insight may help us to identify potential liabilities with future novel compounds.

Expression in yeast links field polymorphisms in PfATP6 to in vitro artemisinin resistance and identifies new inhibitor classes.

BACKGROUND: The mechanism of action of artemisinins against malaria is unclear, despite their widespread use in combination therapies and the emergence of resistance. RESULTS: Here, we report expression of PfATP6 (a SERCA pump) in yeast and demonstrate its inhibition by artemisinins. Mutations in PfATP6 identified in field isolates (such as S769N) and in laboratory clones (such as L263E) decrease susceptibility to artemisinins, whereas they increase susceptibility to unrelated inhibitors such as cyclopiazonic acid. As predicted from the yeast model, Plasmodium falciparum with the L263E mutation is also more susceptible to cyclopiazonic acid. An inability to knockout parasite SERCA pumps provides genetic evidence that they are essential in asexual stages of development. Thaperoxides are a new class of potent antimalarial designed to act by inhibiting PfATP6. Results in yeast confirm this inhibition. CONCLUSIONS: The identification of inhibitors effective against mutated PfATP6 suggests ways in which artemisinin resistance may be overcome.

Computer-aided drug design of falcipain inhibitors: virtual screening, structure-activity relationships, hydration site thermodynamics, and reactivity analysis.

Falcipains (FPs) are hemoglobinases of Plasmodium falciparum that are validated targets for the development of antimalarial chemotherapy. A combined ligand- and structure-based virtual screening of commercial databases was performed to identify structural analogs of virtual screening hits previously discovered in our laboratory. A total of 28 low micromolar inhibitors of FP-2 and FP-3 were identified and the structure-activity relationship (SAR) in each series was elaborated. The SAR of the compounds was unusually steep in some cases and could not be explained by a traditional analysis of the ligand-protein interactions (van der Waals, electrostatics, and hydrogen bonds). To gain further insights, a statistical thermodynamic analysis of explicit solvent in the ligand binding domains of FP-2 and FP-3 was carried out to understand the roles played by water molecules in binding of these inhibitors. Indeed, the energetics associated with the displacement of water molecules upon ligand binding explained some of the complex trends in the SAR. Furthermore, low potency of a subset of FP-2 inhibitors that could not be understood by the water energetics was explained in the context of poor chemical reactivity of the reactive centers of these compounds. The present study highlights the importance of considering energetic contributors to binding beyond traditional ligand-protein interactions.

Inhibitors of SARS-3CLpro: virtual screening, biological evaluation, and molecular dynamics simulation studies.

SARS-CoV from the coronaviridae family has been identified as the etiological agent of Severe Acute Respiratory Syndrome (SARS), a highly contagious upper respiratory disease that reached epidemic status in 2002. SARS-3CL(pro), a cysteine protease indispensible to the viral life cycle, has been identified as one of the key therapeutic targets against SARS. A combined ligand and structure-based virtual screening was carried out against the Asinex Platinum collection. Multiple low micromolar inhibitors of the enzyme were identified through this search, one of which also showed activity against SARS-CoV in a whole cell CPE assay. Furthermore, multinanosecond explicit solvent simulations were carried out using the docking poses of the identified hits to study the overall stability of the binding site interactions as well as identify important changes in the interaction profile that were not apparent from the docking study. Cumulative analysis of the evaluated compounds and the simulation studies led to the identification of certain protein-ligand interaction patterns which would be useful in further structure based design efforts.

Identification of novel malarial cysteine protease inhibitors using structure-based virtual screening of a focused cysteine protease inhibitor library.

Malaria, in particular that caused by Plasmodium falciparum , is prevalent across the tropics, and its medicinal control is limited by widespread drug resistance. Cysteine proteases of P. falciparum , falcipain-2 (FP-2) and falcipain-3 (FP-3), are major hemoglobinases, validated as potential antimalarial drug targets. Structure-based virtual screening of a focused cysteine protease inhibitor library built with soft rather than hard electrophiles was performed against an X-ray crystal structure of FP-2 using the Glide docking program. An enrichment study was performed to select a suitable scoring function and to retrieve potential candidates against FP-2 from a large chemical database. Biological evaluation of 50 selected compounds identified 21 diverse nonpeptidic inhibitors of FP-2 with a hit rate of 42%. Atomic Fukui indices were used to predict the most electrophilic center and its electrophilicity in the identified hits. Comparison of predicted electrophilicity of electrophiles in identified hits with those in known irreversible inhibitors suggested the soft-nature of electrophiles in the selected target compounds. The present study highlights the importance of focused libraries and enrichment studies in structure-based virtual screening. In addition, few compounds were screened against homologous human cysteine proteases for selectivity analysis. Further evaluation of structure-activity relationships around these nonpeptidic scaffolds could help in the development of selective leads for antimalarial chemotherapy.

Design, synthesis, and development of novel guaianolide-endoperoxides as potential antimalarial agents.

Design and synthesis of a guaianolide-endoperoxide (thaperoxide) 3 was pursued as a new antimalarial lead which was found to be noncytotoxic as compared to the natural product lead thapsigargin 2. Several analogues of 3 were successfully synthesized and found to be comparable to derivatives of artemisinin 1 in in vitro antimalarial assay. Among the synthesized compounds, 22 showed excellent in vitro potency against the cultured parasites (W2 IC(50) = 13 nM) without apparent cytotoxicity. Furthermore, SAR trends in thaperoxide analogues are presented and explained with the help of docking studies in the homology model of PfSERCA(PfATP6).

Computational approaches for the discovery of cysteine protease inhibitors against malaria and SARS.

Cysteine proteases are implicated in a variety of human physiological processes and also form an essential component of the life cycle of a number of pathogenic protozoa and viruses. The present review highlights the drug design approaches utilized to understand the mechanism of inhibition and discovery of inhibitors against protozoal cysteine protease, falcipain (a cysteine protease of P. falciparum which causes malaria), and viral cysteine protease, SARS-CoV M(pro) (a cysteine protease of severe acute respiratory syndrome corona virus). The article describes rational approaches for the design of inhibitors and focuses on a variety of structure as well as ligand-based modeling strategies adopted for the discovery of the inhibitors. Also, the key features of ligand recognition against these targets are accentuated. Although no apparent similarities exist between viral and protozoal cysteine proteases discussed here, the goal is to provide examples of rational drug design approaches adopted to design inhibitors against these proteases. The current review would be of interest to scientists engaged in the development of drug design strategies to target the cysteine proteases present in mammals and other lower order organisms.

Design, synthesis, and docking studies of novel benzimidazoles for the treatment of metabolic syndrome.

In addition to lowering blood pressure, telmisartan, an angiotensin (AT(1)) receptor blocker, has recently been shown to exert pleiotropic effects as a partial agonist of nuclear peroxisome proliferator-activated receptor gamma (PPAR gamma). On the basis of these findings and docking pose similarity between telmisartan and rosiglitazone in PPAR gamma active site, two classes of benzimidazole derivatives were designed and synthesized as dual PPAR gamma agonist/angiotensin II antagonists for the possible treatment of metabolic syndrome. Compound 4, a bisbenzimidazole derivative showed the best affinity for the AT(1) receptor with a K(i) = 13.4 nM, but it was devoid of PPAR gamma activity. On the other hand 9, a monobenzimidazole derivative, showed the highest activity in PPAR gamma transactivation assay (69% activation) with no affinity for the AT(1) receptor. Docking studies lead to the designing of a molecule with dual activity, 10, with moderate PPARgamma activity (29%) and affinity for the AT(1) receptor (K(i) = 2.5 microM).