Emanuele Perola

A detailed comparison of current docking and scoring methods on systems of pharmaceutical relevance.

A thorough evaluation of some of the most advanced docking and scoring methods currently available is described, and guidelines for the choice of an appropriate protocol for docking and virtual screening are defined. The generation of a large and highly curated test set of pharmaceutically relevant protein–ligand complexes with known binding affinities is described, and three highly regarded docking programs (Glide, GOLD, and ICM) are evaluated on the same set with respect to their ability to reproduce crystallographic binding orientations. Glide correctly identified the crystallographic pose within 2.0 Å in 61% of the cases, versus 48% for GOLD and 45% for ICM. In general Glide appears to perform most consistently with respect to diversity of binding sites and ligand flexibility, while the performance of ICM and GOLD is more binding site–dependent and it is significantly poorer when binding is predominantly driven by hydrophobic interactions. The results also show that energy minimization and reranking of the top N poses can be an effective means to overcome some of the limitations of a given docking function. The same docking programs are evaluated in conjunction with three different scoring functions for their ability to discriminate actives from inactives in virtual screening. The evaluation, performed on three different systems (HIV‐1 protease, IMPDH, and p38 MAP kinase), confirms that the relative performance of different docking and scoring methods is to some extent binding site–dependent. GlideScore appears to be an effective scoring function for database screening, with consistent performance across several types of binding sites, while ChemScore appears to be most useful in sterically demanding sites since it is more forgiving of repulsive interactions. Energy minimization of docked poses can significantly improve the enrichments in systems with sterically demanding binding sites. Overall Glide appears to be a safe general choice for docking, while the choice of the best scoring tool remains to a larger extent system‐dependent and should be evaluated on a case‐by‐case basis. Proteins 2004.

Novel Dual-Targeting Benzimidazole Urea Inhibitors of DNA Gyrase and Topoisomerase IV Possessing Potent Antibacterial Activity: Intelligent Design and Evolution through the Judicious Use of Structure-Guided Design and Stucture−Activity Relationships

The discovery of new antibacterial agents with novel mechanisms of action is necessary to overcome the problem of bacterial resistance that affects all currently used classes of antibiotics. Bacterial DNA gyrase and topoisomerase IV are well-characterized clinically validated targets of the fluoroquinolone antibiotics which exert their antibacterial activity through inhibition of the catalytic subunits. Inhibition of these targets through interaction with their ATP sites has been less clinically successful. The discovery and characterization of a new class of low molecular weight, synthetic inhibitors of gyrase and topoisomerase IV that bind to the ATP sites are presented. The benzimidazole ureas are dual targeting inhibitors of both enzymes and possess potent antibacterial activity against a wide spectrum of relevant pathogens responsible for hospital- and community-acquired infections. The discovery and optimization of this novel class of antibacterials by the use of structure-guided design, modeling, and structure-activity relationships are described. Data are presented for enzyme inhibition, antibacterial activity, and in vivo efficacy by oral and intravenous administration in two rodent infection models.

An Analysis of the Binding Efficiencies of Drugs and Their Leads in Successful Drug Discovery Programs

In order to investigate the evolution of binding efficiency in successful drug discovery programs, a data set of 60 lead/drug pairs with known binding affinities has been compiled and analyzed. Low-end thresholds for the binding efficiencies of viable leads and drugs have been derived. On average, the drugs in the set are significantly larger and more potent but have similar lipophilicity relative to their originating leads, suggesting that the ability to maintain low levels of lipophilicity while increasing molecular weight is one of the keys to a successful drug discovery program. A number of examples demonstrate that large increases in binding efficiency from leads to more elaborate drugs sharing the same scaffold can be achieved. The importance of dissecting a lead structure to identify the most efficient fragments and the option of sacrificing binding efficiency to optimize other properties are discussed, and relevant examples are highlighted.

Novel inhibitors of cholesterol degradation in Mycobacterium tuberculosis reveal how the bacterium's metabolism is constrained by the intracellular environment.

Mycobacterium tuberculosis (Mtb) relies on a specialized set of metabolic pathways to support growth in macrophages. By conducting an extensive, unbiased chemical screen to identify small molecules that inhibit Mtb metabolism within macrophages, we identified a significant number of novel compounds that limit Mtb growth in macrophages and in medium containing cholesterol as the principle carbon source. Based on this observation, we developed a chemical-rescue strategy to identify compounds that target metabolic enzymes involved in cholesterol metabolism. This approach identified two compounds that inhibit the HsaAB enzyme complex, which is required for complete degradation of the cholesterol A/B rings. The strategy also identified an inhibitor of PrpC, the 2-methylcitrate synthase, which is required for assimilation of cholesterol-derived propionyl-CoA into the TCA cycle. These chemical probes represent new classes of inhibitors with novel modes of action, and target metabolic pathways required to support growth of Mtb in its host cell. The screen also revealed a structurally-diverse set of compounds that target additional stage(s) of cholesterol utilization. Mutants resistant to this class of compounds are defective in the bacterial adenylate cyclase Rv1625/Cya. These data implicate cyclic-AMP (cAMP) in regulating cholesterol utilization in Mtb, and are consistent with published reports indicating that propionate metabolism is regulated by cAMP levels. Intriguingly, reversal of the cholesterol-dependent growth inhibition caused by this subset of compounds could be achieved by supplementing the media with acetate, but not with glucose, indicating that Mtb is subject to a unique form of metabolic constraint induced by the presence of cholesterol.

Discovery of a Novel, First-in-Class, Orally Bioavailable Azaindole Inhibitor (VX-787) of Influenza PB2.

In our effort to develop agents for the treatment of influenza, a phenotypic screening approach utilizing a cell protection assay identified a series of azaindole based inhibitors of the cap-snatching function of the PB2 subunit of the influenza A viral polymerase complex. Using a bDNA viral replication assay (Wagaman, P. C., Leong, M. A., and Simmen, K. A. Development of a novel influenza A antiviral assay. J. Virol. Methods 2002, 105, 105-114) in cells as a direct measure of antiviral activity, we discovered a set of cyclohexyl carboxylic acid analogues, highlighted by VX-787 (2). Compound 2 shows strong potency versus multiple influenza A strains, including pandemic 2009 H1N1 and avian H5N1 flu strains, and shows an efficacy profile in a mouse influenza model even when treatment was administered 48 h after infection. Compound 2 represents a first-in-class, orally bioavailable, novel compound that offers potential for the treatment of both pandemic and seasonal influenza and has a distinct advantage over the current standard of care treatments including potency, efficacy, and extended treatment window.

Preclinical Activity of VX-787, a First-in-Class, Orally Bioavailable Inhibitor of the Influenza Virus Polymerase PB2 Subunit

ABSTRACT VX-787 is a novel inhibitor of influenza virus replication that blocks the PB2 cap-snatching activity of the influenza viral polymerase complex. Viral genetics and X-ray crystallography studies provide support for the idea that VX-787 occupies the 7-methyl GTP (m7GTP) cap-binding site of PB2. VX-787 binds the cap-binding domain of the PB2 subunit with a KD (dissociation constant) of 24 nM as determined by isothermal titration calorimetry (ITC). The cell-based EC50 (the concentration of compound that ensures 50% cell viability of an uninfected control) for VX-787 is 1.6 nM in a cytopathic effect (CPE) assay, with a similar EC50 in a viral RNA replication assay. VX-787 is active against a diverse panel of influenza A virus strains, including H1N1pdm09 and H5N1 strains, as well as strains with reduced susceptibility to neuraminidase inhibitors (NAIs). VX-787 was highly efficacious in both prophylaxis and treatment models of mouse influenza and was superior to the neuraminidase inhibitor, oseltamivir, including in delayed-start-to-treat experiments, with 100% survival at up to 96 h postinfection and partial survival in groups where the initiation of therapy was delayed up to 120 h postinfection. At different doses, VX-787 showed a 1-log to >5-log reduction in viral load (relative to vehicle controls) in mouse lungs. Overall, these favorable findings validate the PB2 subunit of the viral polymerase as a drug target for influenza therapy and support the continued development of VX-787 as a novel antiviral agent for the treatment of influenza infection.

Minimizing false positives in kinase virtual screens

In spite of recent improvements in docking and scoring methods, high false‐positive rates remain a common issue in structure‐based virtual screening. In this study, the distinctive features of false positives in kinase virtual screens were investigated. A series of retrospective virtual screens on kinase targets was performed on specifically designed test sets, each combining true ligands and experimentally confirmed inactive compounds. A systematic analysis of the docking poses generated for the top‐ranking compounds highlighted key aspects differentiating true hits from false positives. The most recurring feature in the poses of false positives was the absence of certain key interactions known to be required for kinase binding. A systematic analysis of 444 crystal structures of ligand‐bound kinases showed that at least two hydrogen bonds between the ligand and the backbone protein atoms in the kinase hinge region are present in 90% of the complexes, with very little variability across targets. Closer inspection showed that when the two hydrogen bonds are present, one of three preferred hinge‐binding motifs is involved in 96.5% of the cases. Less than 10% of the false positives satisfied these two criteria in the minimized docking poses generated by our standard protocol. Ligand conformational artifacts were also shown to contribute to the occurrence of false positives in a number of cases. Application of this knowledge in the form of docking constraints and post‐processing filters provided consistent improvements in virtual screening performance on all systems. The false‐positive rates were significantly reduced and the enrichment factors increased by an average of twofold. On the basis of these results, a generalized two‐step protocol for virtual screening on kinase targets is suggested. Proteins 2006.

Development of a Rule-Based Method for the Assessment of Protein Druggability

Target selection is a critical step in the majority of modern drug discovery programs. The viability of a drug target depends on two components: biological relevance and chemical tractability. The concept of druggability was introduced to describe the second component, and it is defined as the ability of a target to bind a drug-like molecule with a therapeutically useful level of affinity. To investigate the rules that govern druggability, we developed an algorithm to isolate and characterize the binding pockets of protein targets. Using this algorithm, we performed a comparative analysis between the relevant pockets of 60 targets of approved drugs and a diverse set of 440 ligand-binding pockets. As a result, we defined a preferred property space for druggable pockets based on five key properties (volume, depth, enclosure, percentage of charged residues and hydrophobicity), and we represented it with a set of simple rules. These rules may be applicable in the future to evaluate the chemical tractability of prospective targets.

Second-generation antibacterial benzimidazole ureas: discovery of a preclinical candidate with reduced metabolic liability.

Compound 3 is a potent aminobenzimidazole urea with broad-spectrum Gram-positive antibacterial activity resulting from dual inhibition of bacterial gyrase (GyrB) and topoisomerase IV (ParE), and it demonstrates efficacy in rodent models of bacterial infection. Preclinical in vitro and in vivo studies showed that compound 3 covalently labels liver proteins, presumably via formation of a reactive metabolite, and hence presented a potential safety liability. The urea moiety in compound 3 was identified as being potentially responsible for reactive metabolite formation, but its replacement resulted in loss of antibacterial activity and/or oral exposure due to poor physicochemical parameters. To identify second-generation aminobenzimidazole ureas devoid of reactive metabolite formation potential, we implemented a metabolic shift strategy, which focused on shifting metabolism away from the urea moiety by introducing metabolic soft spots elsewhere in the molecule. Aminobenzimidazole urea 34, identified through this strategy, exhibits similar antibacterial activity as that of 3 and did not label liver proteins in vivo, indicating reduced/no potential for reactive metabolite formation.

A Novel Inhibitor of Gyrase B Is a Potent Drug Candidate for Treatment of Tuberculosis and Nontuberculosis Mycobacterial Infections

ABSTRACT New drugs to treat drug-resistant tuberculosis are urgently needed. Extensively drug-resistant and probably the totally drug-resistant tuberculosis strains are resistant to fluoroquinolones like moxifloxacin, which target gyrase A, and most people infected with these strains die within a year. In this study, we found that a novel aminobenzimidazole, VXc-486, which targets gyrase B, potently inhibits multiple drug-sensitive isolates and drug-resistant isolates of Mycobacterium tuberculosis in vitro (MICs of 0.03 to 0.30 μg/ml and 0.08 to 5.48 μg/ml, respectively) and reduces mycobacterial burdens in lungs of infected mice in vivo. VXc-486 is active against drug-resistant isolates, has bactericidal activity, and kills intracellular and dormant M. tuberculosis bacteria in a low-oxygen environment. Furthermore, we found that VXc-486 inhibits the growth of multiple strains of Mycobacterium abscessus, Mycobacterium avium complex, and Mycobacterium kansasii (MICs of 0.1 to 2.0 μg/ml), as well as that of several strains of Nocardia spp. (MICs of 0.1 to 1.0 μg/ml). We made a direct comparison of the parent compound VXc-486 and a phosphate prodrug of VXc-486 and showed that the prodrug of VXc-486 had more potent killing of M. tuberculosis than did VXc-486 in vivo. In combination with other antimycobacterial drugs, the prodrug of VXc-486 sterilized M. tuberculosis infection when combined with rifapentine-pyrazinamide and bedaquiline-pyrazinamide in a relapse infection study in mice. Furthermore, the prodrug of VXc-486 appeared to perform at least as well as the gyrase A inhibitor moxifloxacin. These findings warrant further development of the prodrug of VXc-486 for the treatment of tuberculosis and nontuberculosis mycobacterial infections.