Fragment-Based Analysis of Ligand Dockings Improves Classification of Actives.

We describe ADChemCast, a method for using results from virtual screening to create a richer representation of a target binding site, which may be used to improve ranking of compounds and characterize the determinants of ligand-receptor specificity. ADChemCast clusters docked conformations of ligands based on shared pairwise receptor-ligand interactions within chemically similar structural fragments, building a set of attributes characteristic of binders and nonbinders. Machine learning is then used to build rules from the most informational attributes for use in reranking of compounds. In this report, we use ADChemCast to improve the ranking of compounds in 11 diverse proteins from the Database of Useful Decoys-Enhanced (DUD-E) and demonstrate the utility of the method for characterizing relevant binding attributes in HIV reverse transcriptase.

Proteome-wide covalent ligand discovery in native biological systems.

Small molecules are powerful tools for investigating protein function and can serve as leads for new therapeutics. Most human proteins, however, lack small-molecule ligands, and entire protein classes are considered 'undruggable'. Fragment-based ligand discovery can identify small-molecule probes for proteins that have proven difficult to target using high-throughput screening of complex compound libraries. Although reversibly binding ligands are commonly pursued, covalent fragments provide an alternative route to small-molecule probes, including those that can access regions of proteins that are difficult to target through binding affinity alone. Here we report a quantitative analysis of cysteine-reactive small-molecule fragments screened against thousands of proteins in human proteomes and cells. Covalent ligands were identified for >700 cysteines found in both druggable proteins and proteins deficient in chemical probes, including transcription factors, adaptor/scaffolding proteins, and uncharacterized proteins. Among the atypical ligand-protein interactions discovered were compounds that react preferentially with pro- (inactive) caspases. We used these ligands to distinguish extrinsic apoptosis pathways in human cell lines versus primary human T cells, showing that the former is largely mediated by caspase-8 while the latter depends on both caspase-8 and -10. Fragment-based covalent ligand discovery provides a greatly expanded portrait of the ligandable proteome and furnishes compounds that can illuminate protein functions in native biological systems.

Computational protein-ligand docking and virtual drug screening with the AutoDock suite.

Computational docking can be used to predict bound conformations and free energies of binding for small-molecule ligands to macromolecular targets. Docking is widely used for the study of biomolecular interactions and mechanisms, and it is applied to structure-based drug design. The methods are fast enough to allow virtual screening of ligand libraries containing tens of thousands of compounds. This protocol covers the docking and virtual screening methods provided by the AutoDock suite of programs, including a basic docking of a drug molecule with an anticancer target, a virtual screen of this target with a small ligand library, docking with selective receptor flexibility, active site prediction and docking with explicit hydration. The entire protocol will require ∼5 h.

A virtual screen discovers novel, fragment-sized inhibitors of Mycobacterium tuberculosis InhA.

Isoniazid (INH) is usually administered to treat latent Mycobacterium tuberculosis (Mtb) infections and is used in combination therapy to treat active tuberculosis (TB). Unfortunately, resistance to this drug is hampering its clinical effectiveness. INH is a prodrug that must be activated by Mtb catalase-peroxidase (KatG) before it can inhibit InhA (Mtb enoyl-acyl-carrier-protein reductase). Isoniazid-resistant cases of TB found in clinical settings usually involve mutations in or deletion of katG, which abrogate INH activation. Compounds that inhibit InhA without requiring prior activation by KatG would not be affected by this resistance mechanism and hence would display continued potency against these drug-resistant isolates of Mtb. Virtual screening experiments versus InhA in the GO Fight Against Malaria (GO FAM) project were designed to discover new scaffolds that display base-stacking interactions with the NAD cofactor. GO FAM experiments included targets from other pathogens, including Mtb, when they had structural similarity to a malaria target. Eight of the 16 soluble compounds identified by docking against InhA plus visual inspection were modest inhibitors and did not require prior activation by KatG. The best two inhibitors discovered are both fragment-sized compounds and displayed Ki values of 54 and 59 µM, respectively. Importantly, the novel inhibitors discovered have low structural similarity to known InhA inhibitors and thus help expand the number of chemotypes on which future medicinal chemistry efforts can be focused. These new fragment hits could eventually help advance the fight against INH-resistant Mtb strains, which pose a significant global health threat.

Distinguishing binders from false positives by free energy calculations: fragment screening against the flap site of HIV protease.

Molecular docking is a powerful tool used in drug discovery and structural biology for predicting the structures of ligand-receptor complexes. However, the accuracy of docking calculations can be limited by factors such as the neglect of protein reorganization in the scoring function; as a result, ligand screening can produce a high rate of false positive hits. Although absolute binding free energy methods still have difficulty in accurately rank-ordering binders, we believe that they can be fruitfully employed to distinguish binders from nonbinders and reduce the false positive rate. Here we study a set of ligands that dock favorably to a newly discovered, potentially allosteric site on the flap of HIV-1 protease. Fragment binding to this site stabilizes a closed form of protease, which could be exploited for the design of allosteric inhibitors. Twenty-three top-ranked protein-ligand complexes from AutoDock were subject to the free energy screening using two methods, the recently developed binding energy analysis method (BEDAM) and the standard double decoupling method (DDM). Free energy calculations correctly identified most of the false positives (≥83%) and recovered all the confirmed binders. The results show a gap averaging ≥3.7 kcal/mol, separating the binders and the false positives. We present a formula that decomposes the binding free energy into contributions from the receptor conformational macrostates, which provides insights into the roles of different binding modes. Our binding free energy component analysis further suggests that improving the treatment for the desolvation penalty associated with the unfulfilled polar groups could reduce the rate of false positive hits in docking. The current study demonstrates that the combination of docking with free energy methods can be very useful for more accurate ligand screening against valuable drug targets.

Crystallographic fragment-based drug discovery: use of a brominated fragment library targeting HIV protease.

A library of 68 brominated fragments was screened against a new crystal form of inhibited HIV-1 protease in order to probe surface sites in soaking experiments. Often, fragments are weak binders with partial occupancy, resulting in weak, difficult-to-fit electron density. The use of a brominated fragment library addresses this challenge, as bromine can be located unequivocally via anomalous scattering. Data collection was carried out in an automated fashion using AutoDrug at SSRL. Novel hits were identified in the known surface sites: 3-bromo-2,6-dimethoxybenzoic acid (Br6) in the flap site and 1-bromo-2-naphthoic acid (Br27) in the exosite, expanding the chemistry of known fragments for development of higher affinity potential allosteric inhibitors. At the same time, mapping the binding sites of a number of weaker binding Br-fragments provides further insight into the nature of these surface pockets.

Virtual screening against acetylcholine binding protein.

The nicotinic acetylcholine receptors (nAChRs) are a member of the ligand-gated ion channel family and play a key role in the transfer of information across neurological networks. The X-ray crystal structure of agonist-bound α(7) acetylcholine binding protein (AChBP) has been recognized as the most appropriate template to model the ligand-binding domain of nAChR for studying the molecular mechanism of the receptor-ligand interactions. Virtual screening of the National Cancer Institute diversity set, a library of 1990 compounds with nonredundant pharmacophore profiles, using AutoDock against AChBPs revealed 51 potential candidates. In vitro radioligand competition assays using [(3)H] epibatidine against the AChBPs from the freshwater snails, Lymnaea stagnalis, and from the marine species, Aplysia californica and the mutant (AcY55W), revealed seven compounds from the list of candidates that had micromolar to nanomolar affinities for the AChBPs. Further investigation on α(7)nAChR expressing in Xenopus oocytes and on the recombinant receptors with fluorescence resonance energy transfer (FRET)-based calcium sensor expressing in HEK cells showed that seven compounds were antagonists of α(7)nAChR, only one compound (NSC34352) demonstrated partial agonistic effect at low dose (10 µM), and two compounds (NSC36369 and NSC34352) were selective antagonists on α(7)nAchR with moderate potency. These hits serve as novel templates/scaffolds for development of more potent and specific in the AChR systems.

Virtual screening against alpha-cobratoxin.

alpha-Cobratoxin (Cbtx), the neurotoxin isolated from the venom of the Thai cobra Naja kaouthia , causes paralysis by preventing acetylcholine (ACh) binding to nicotinic acetylcholine receptors (nAChRs). In the current study, the region of the Cbtx molecule that is directly involved in binding to nAChRs is used as the target for anticobratoxin drug design. The crystal structure (1YI5) of Cbtx in complex with the acetylcholine binding protein (AChBP), a soluble homolog of the extracellular binding domain of nAChRs, was selected to prepare an alpha-cobratoxin active binding site for docking. The amino acid residues (Ser182-Tyr192) of the AChBP structure, the binding site of Cbtx, were used as the positive control to validate the prepared Cbtx active binding site (root mean square deviation < 1.2 A). Virtual screening of the National Cancer Institute diversity set, a library of 1990 compounds with nonredundant pharmacophore profiles, using AutoDock against the Cbtx active site, revealed 39 potential inhibitor candidates. The adapted in vitro radioligand competition assays using [(3)H]epibatidine and [(125)I]bungarotoxin against the AChBPs from the marine species, Aplysia californica (Ac), and from the freshwater snails, Lymnaea stagnalis (Ls) and Bolinus truncates (Bt), revealed 4 compounds from the list of inhibitor candidates that had micromolar to nanomolar interferences for the toxin binding to AChBPs. Three hits (NSC42258, NSC121865, and NSC134754) can prolong the survival time of the mice if administered 30 min before injection with Cbtx, but only NSC121865 and NSC134754 can prolong the survival time if injected immediately after injection with Cbtx. These inhibitors serve as novel templates/scaffolds for the development of more potent and specific anticobratoxin.

Pursuing aldose reductase inhibitors through in situ cross-docking and similarity-based virtual screening.

Aldose reductase (ALR2) is a critical enzyme in the development of the major complications of diabetes mellitus. Herein, new molecular entities active against ALR2 were discovered through an integrated receptor- and ligand-based virtual screening campaign. Twelve candidates were found to inhibit this enzyme in the micromolar range including two ligands having an IC(50) below 3 muM. Six new compounds, structurally unrelated to the known ARIs, have been identified, opening up opportunity for lead optimization.

Structure-based virtual screening and biological evaluation of Mycobacterium tuberculosis adenosine 5'-phosphosulfate reductase inhibitors.

Tuberculosis is among the world's deadliest infectious diseases. APS reductase catalyzes the first committed step in bacterial sulfate reduction and is a validated drug target against latent tuberculosis infection. We performed a virtual screening to identify APSR inhibitors. These inhibitors represent the first non-phosphate-based molecules to inhibit APSR. Common chemical features lay the foundation for the development of agents that could shorten the duration of chemotherapy by targeting the latent stage of TB infection.