Ultra-large library docking for discovering new chemotypes.

Despite intense interest in expanding chemical space, libraries containing hundreds-of-millions to billions of diverse molecules have remained inaccessible. Here we investigate structure-based docking of 170 million make-on-demand compounds from 130 well-characterized reactions. The resulting library is diverse, representing over 10.7 million scaffolds that are otherwise unavailable. For each compound in the library, docking against AmpC ß-lactamase (AmpC) and the D4 dopamine receptor were simulated. From the top-ranking molecules, 44 and 549 compounds were synthesized and tested for interactions with AmpC and the D4 dopamine receptor, respectively. We found a phenolate inhibitor of AmpC, which revealed a group of inhibitors without known precedent. This molecule was optimized to 77 nM, which places it among the most potent non-covalent AmpC inhibitors known. Crystal structures of this and other AmpC inhibitors confirmed the docking predictions. Against the D4 dopamine receptor, hit rates fell almost monotonically with docking score, and a hit-rate versus score curve predicted that the library contained 453,000 ligands for the D4 dopamine receptor. Of 81 new chemotypes discovered, 30 showed submicromolar activity, including a 180-pM subtype-selective agonist of the D4 dopamine receptor.

Energy penalties enhance flexible receptor docking in a model cavity.

Protein flexibility remains a major challenge in library docking because of difficulties in sampling conformational ensembles with accurate probabilities. Here, we use the model cavity site of T4 lysozyme L99A to test flexible receptor docking with energy penalties from molecular dynamics (MD) simulations. Crystallography with larger and smaller ligands indicates that this cavity can adopt three major conformations: open, intermediate, and closed. Since smaller ligands typically bind better to the cavity site, we anticipate an energy penalty for the cavity opening. To estimate its magnitude, we calculate conformational preferences from MD simulations. We find that including a penalty term is essential for retrospective ligand enrichment; otherwise, high-energy states dominate the docking. We then prospectively docked a library of over 900,000 compounds for new molecules binding to each conformational state. Absent a penalty term, the open conformation dominated the docking results; inclusion of this term led to a balanced sampling of ligands against each state. High ranked molecules were experimentally tested by Tm upshift and X-ray crystallography. From 33 selected molecules, we identified 18 ligands and determined 13 crystal structures. Most interesting were those bound to the open cavity, where the buried site opens to bulk solvent. Here, highly unusual ligands for this cavity had been predicted, including large ligands with polar tails; these were confirmed both by binding and by crystallography. In docking, incorporating protein flexibility with thermodynamic weightings may thus access new ligand chemotypes. The MD approach to accessing and, crucially, weighting such alternative states may find general applicability.

Property-Unmatched Decoys in Docking Benchmarks.

Enrichment of ligands versus property-matched decoys is widely used to test and optimize docking library screens. However, the unconstrained optimization of enrichment alone can mislead, leading to false confidence in prospective performance. This can arise by over-optimizing for enrichment against property-matched decoys, without considering the full spectrum of molecules to be found in a true large library screen. Adding decoys representing charge extrema helps mitigate over-optimizing for electrostatic interactions. Adding decoys that represent the overall characteristics of the library to be docked allows one to sample molecules not represented by ligands and property-matched decoys but that one will encounter in a prospective screen. An optimized version of the DUD-E set (DUDE-Z), as well as Extrema and sets representing broad features of the library (Goldilocks), is developed here. We also explore the variability that one can encounter in enrichment calculations and how that can temper one's confidence in small enrichment differences. The new tools and new decoy sets are freely available at http://tldr.docking.org and http://dudez.docking.org.

Discovery of Lysine-Targeted eIF4E Inhibitors through Covalent Docking.

Eukaryotic translation initiation factor 4E (eIF4E) binds the m7GTP cap structure at the 5'-end of mRNAs, stimulating the translation of proteins implicated in cancer cell growth and metastasis. eIF4E is a notoriously challenging target, and most of the reported inhibitors are negatively charged guanine analogues with negligible cell permeability. To overcome these challenges, we envisioned a covalent targeting strategy. As there are no cysteines near the eIF4E cap binding site, we developed a covalent docking approach focused on lysine. Taking advantage of a "make-on-demand" virtual library, we used covalent docking to identify arylsulfonyl fluorides that target a noncatalytic lysine (Lys162) in eIF4E. Guided by cocrystal structures, we elaborated arylsulfonyl fluoride 2 to 12, which to our knowledge is the first covalent eIF4E inhibitor with cellular activity. In addition to providing a new tool for acutely inactivating eIF4E in cells, our computational approach may offer a general strategy for developing selective lysine-targeted covalent ligands.

Testing inhomogeneous solvation theory in structure-based ligand discovery.

Binding-site water is often displaced upon ligand recognition, but is commonly neglected in structure-based ligand discovery. Inhomogeneous solvation theory (IST) has become popular for treating this effect, but it has not been tested in controlled experiments at atomic resolution. To do so, we turned to a grid-based version of this method, GIST, readily implemented in molecular docking. Whereas the term only improves docking modestly in retrospective ligand enrichment, it could be added without disrupting performance. We thus turned to prospective docking of large libraries to investigate GIST's impact on ligand discovery, geometry, and water structure in a model cavity site well-suited to exploring these terms. Although top-ranked docked molecules with and without the GIST term often overlapped, many ligands were meaningfully prioritized or deprioritized; some of these were selected for testing. Experimentally, 13/14 molecules prioritized by GIST did bind, whereas none of the molecules that it deprioritized were observed to bind. Nine crystal complexes were determined. In six, the ligand geometry corresponded to that predicted by GIST, for one of these the pose without the GIST term was wrong, and three crystallographic poses differed from both predictions. Notably, in one structure, an ordered water molecule with a high GIST displacement penalty was observed to stay in place. Inclusion of this water-displacement term can substantially improve the hit rates and ligand geometries from docking screens, although the magnitude of its effects can be small and its impact in drug binding sites merits further controlled studies.

Novel K-Ras G12C Switch-II Covalent Binders Destabilize Ras and Accelerate Nucleotide Exchange.

The success of targeted covalent inhibitors in the global pharmaceutical industry has led to a resurgence of covalent drug discovery. However, covalent inhibitor design for flexible binding sites remains a difficult task due to a lack of methodological development. Here, we compared covalent docking to empirical electrophile screening against the highly dynamic target K-RasG12C. While the overall hit rate of both methods was comparable, we were able to rapidly progress a docking hit to a potent irreversible covalent binder that modifies the inactive, GDP-bound state of K-RasG12C. Hydrogen-deuterium exchange mass spectrometry was used to probe the protein dynamics of compound binding to the switch-II pocket and subsequent destabilization of the nucleotide-binding region. SOS-mediated nucleotide exchange assays showed that, contrary to prior switch-II pocket inhibitors, these new compounds appear to accelerate nucleotide exchange. This study highlights the efficiency of covalent docking as a tool for the discovery of chemically novel hits against challenging targets.

Homologous ligands accommodated by discrete conformations of a buried cavity.

Conformational change in protein-ligand complexes is widely modeled, but the protein accommodation expected on binding a congeneric series of ligands has received less attention. Given their use in medicinal chemistry, there are surprisingly few substantial series of congeneric ligand complexes in the Protein Data Bank (PDB). Here we determine the structures of eight alkyl benzenes, in single-methylene increases from benzene to n-hexylbenzene, bound to an enclosed cavity in T4 lysozyme. The volume of the apo cavity suffices to accommodate benzene but, even with toluene, larger cavity conformations become observable in the electron density, and over the series two other major conformations are observed. These involve discrete changes in main-chain conformation, expanding the site; few continuous changes in the site are observed. In most structures, two discrete protein conformations are observed simultaneously, and energetic considerations suggest that these conformations are low in energy relative to the ground state. An analysis of 121 lysozyme cavity structures in the PDB finds that these three conformations dominate the previously determined structures, largely modeled in a single conformation. An investigation of the few congeneric series in the PDB suggests that discrete changes are common adaptations to a series of growing ligands. The discrete, but relatively few, conformational states observed here, and their energetic accessibility, may have implications for anticipating protein conformational change in ligand design.

Customizable de novo design strategies for DOCK: Application to HIVgp41 and other therapeutic targets.

De novo design can be used to explore vast areas of chemical space in computational lead discovery. As a complement to virtual screening, from-scratch construction of molecules is not limited to compounds in pre-existing vendor catalogs. Here, we present an iterative fragment growth method, integrated into the program DOCK, in which new molecules are built using rules for allowable connections based on known molecules. The method leverages DOCK's advanced scoring and pruning approaches and users can define very specific criteria in terms of properties or features to customize growth toward a particular region of chemical space. The code was validated using three increasingly difficult classes of calculations: (1) Rebuilding known X-ray ligands taken from 663 complexes using only their component parts (focused libraries), (2) construction of new ligands in 57 drug target sites using a library derived from ∼13M drug-like compounds (generic libraries), and (3) application to a challenging protein-protein interface on the viral drug target HIVgp41. The computational testing confirms that the de novo DOCK routines are robust and working as envisioned, and the compelling results highlight the potential utility for designing new molecules against a wide variety of important protein targets. © 2017 Wiley Periodicals, Inc.

DOCK 6: Impact of new features and current docking performance.

This manuscript presents the latest algorithmic and methodological developments to the structure-based design program DOCK 6.7 focused on an updated internal energy function, new anchor selection control, enhanced minimization options, a footprint similarity scoring function, a symmetry-corrected root-mean-square deviation algorithm, a database filter, and docking forensic tools. An important strategy during development involved use of three orthogonal metrics for assessment and validation: pose reproduction over a large database of 1043 protein-ligand complexes (SB2012 test set), cross-docking to 24 drug-target protein families, and database enrichment using large active and decoy datasets (Directory of Useful Decoys [DUD]-E test set) for five important proteins including HIV protease and IGF-1R. Relative to earlier versions, a key outcome of the work is a significant increase in pose reproduction success in going from DOCK 4.0.2 (51.4%) → 5.4 (65.2%) → 6.7 (73.3%) as a result of significant decreases in failure arising from both sampling 24.1% → 13.6% → 9.1% and scoring 24.4% → 21.1% → 17.5%. Companion cross-docking and enrichment studies with the new version highlight other strengths and remaining areas for improvement, especially for systems containing metal ions. The source code for DOCK 6.7 is available for download and free for academic users at http://dock.compbio.ucsf.edu/.

Computer-aided identification, synthesis, and biological evaluation of novel inhibitors for botulinum neurotoxin serotype A.

Botulinum neurotoxins (BoNTs) are among the most potent biological toxin known to humans, and are classified as Category A bioterrorism agents by the Centers for Disease Control and prevention (CDC). There are seven known BoNT serotypes (A-G) which have been thus far identified in literature. BoNTs have been shown to block neurotransmitter release by cleaving proteins of the soluble NSF attachment protein receptor (SNARE) complex. Disruption of the SNARE complex precludes motor neuron failure which ultimately results in flaccid paralysis in humans and animals. Currently, there are no effective therapeutic treatments against the neurotoxin light chain (LC) after translocation into the cytosols of motor neurons. In this work, high-throughput in silico screening was employed to screen a library of commercially available compounds from ZINC database against BoNT/A-LC. Among the hit compounds from the in silico screening, two lead compounds were identified and found to have potent inhibitory activity against BoNT/A-LC in vitro, as well as in Neuro-2a cells. A few analogs of the lead compounds were synthesized and their potency examined. One of these analogs showed an enhanced activity than the lead compounds.

Considerations Around Structure-Based Drug Discovery for KRAS Using DOCK.

Molecular docking is a popular computational tool in drug discovery. Leveraging structural information, docking software predicts binding poses of small molecules to cavities on the surfaces of proteins. Virtual screening for ligand discovery is a useful application of docking software. In this chapter, using the enigmatic KRAS protein as an example system, we endeavor to teach the reader about best practices for performing molecular docking with UCSF DOCK. We discuss methods for virtual screening and docking molecules on KRAS. We present the following six points to optimize our docking setup for prosecuting a virtual screen: protein structure choice, pocket selection, optimization of the scoring function, modification of sampling spheres and sampling procedures, choosing an appropriate portion of chemical space to dock, and the choice of which top scoring molecules to pick for purchase.

Structural dynamics of RAF1-HSP90-CDC37 and HSP90 complexes reveal asymmetric client interactions and key structural elements.

RAF kinases are integral to the RAS-MAPK signaling pathway, and proper RAF1 folding relies on its interaction with the chaperone HSP90 and the cochaperone CDC37. Understanding the intricate molecular interactions governing RAF1 folding is crucial for comprehending this process. Here, we present a cryo-EM structure of the closed-state RAF1-HSP90-CDC37 complex, where the C-lobe of the RAF1 kinase domain binds to one side of the HSP90 dimer, and an unfolded N-lobe segment of the RAF1 kinase domain threads through the center of the HSP90 dimer. CDC37 binds to the kinase C-lobe, mimicking the N-lobe with its HxNI motif. We also describe structures of HSP90 dimers without RAF1 and CDC37, displaying only N-terminal and middle domains, which we term the semi-open state. Employing 1 µs atomistic simulations, energetic decomposition, and comparative structural analysis, we elucidate the dynamics and interactions within these complexes. Our quantitative analysis reveals that CDC37 bridges the HSP90-RAF1 interaction, RAF1 binds HSP90 asymmetrically, and that HSP90 structural elements engage RAF1's unfolded region. Additionally, N- and C-terminal interactions stabilize HSP90 dimers, and molecular interactions in HSP90 dimers rearrange between the closed and semi-open states. Our findings provide valuable insight into the contributions of HSP90 and CDC37 in mediating client folding.

Grid-based molecular footprint comparison method for docking and de novo design: application to HIVgp41.

Scoring functions are a critically important component of computer-aided screening methods for the identification of lead compounds during early stages of drug discovery. Here, we present a new multigrid implementation of the footprint similarity (FPS) scoring function that was recently developed in our laboratory which has proven useful for identification of compounds which bind to a protein on a per-residue basis in a way that resembles a known reference. The grid-based FPS method is much faster than its Cartesian-space counterpart, which makes it computationally tractable for on-the-fly docking, virtual screening, or de novo design. In this work, we establish that: (i) relatively few grids can be used to accurately approximate Cartesian space footprint similarity, (ii) the method yields improved success over the standard DOCK energy function for pose identification across a large test set of experimental co-crystal structures, for crossdocking, and for database enrichment, and (iii) grid-based FPS scoring can be used to tailor construction of new molecules to have specific properties, as demonstrated in a series of test cases targeting the viral protein HIVgp41. The method is available in the program DOCK6.

Evaluation of DOCK 6 as a pose generation and database enrichment tool.

In conjunction with the recent American Chemical Society symposium titled "Docking and Scoring: A Review of Docking Programs" the performance of the DOCK6 program was evaluated through (1) pose reproduction and (2) database enrichment calculations on a common set of organizer-specified systems and datasets (ASTEX, DUD, WOMBAT). Representative baseline grid score results averaged over five docking runs yield a relatively high pose identification success rate of 72.5 % (symmetry corrected rmsd) and sampling rate of 91.9 % for the multi site ASTEX set (N = 147) using organizer-supplied structures. Numerous additional docking experiments showed that ligand starting conditions, symmetry, multiple binding sites, clustering, and receptor preparation protocols all affect success. Encouragingly, in some cases, use of more sophisticated scoring and sampling methods yielded results which were comparable (Amber score ligand movable protocol) or exceeded (LMOD score) analogous baseline grid-score results. The analysis highlights the potential benefit and challenges associated with including receptor flexibility and indicates that different scoring functions have system dependent strengths and weaknesses. Enrichment studies with the DUD database prepared using the SB2010 preparation protocol and native ligand pairings yielded individual area under the curve (AUC) values derived from receiver operating characteristic curve analysis ranging from 0.29 (bad enrichment) to 0.96 (good enrichment) with an average value of 0.60 (27/38 have AUC ≥ 0.5). Strong early enrichment was also observed in the critically important 1.0-2.0 % region. Somewhat surprisingly, an alternative receptor preparation protocol yielded comparable results. As expected, semi-random pairings yielded poorer enrichments, in particular, for unrelated receptors. Overall, the breadth and number of experiments performed provide a useful snapshot of current capabilities of DOCK6 as well as starting points to guide future development efforts to further improve sampling and scoring.

Implementation and evaluation of a docking-rescoring method using molecular footprint comparisons.

A docking-rescoring method, based on per-residue van der Waals (VDW), electrostatic (ES), or hydrogen bond (HB) energies has been developed to aid discovery of ligands that have interaction signatures with a target (footprints) similar to that of a reference. Biologically useful references could include known drugs, inhibitors, substrates, transition states, or side-chains that mediate protein-protein interactions. Termed footprint similarity (FPS) score, the method, as implemented in the program DOCK, was validated and characterized using: (1) pose identification, (2) crossdocking, (3) enrichment, and (4) virtual screening. Improvements in pose identification (6­12%) were obtained using footprint-based (FPS(VDW+ES)) vs. standard DOCK (DCE(VDW+ES)) scoring as evaluated on three large datasets (680­775 systems) from the SB2010 database. Enhanced pose identification was also observed using FPS (45.4% or 70.9%) compared with DCE (17.8%) methods to rank challenging crossdocking ensembles from carbonic anhydrase. Enrichment tests, for three representative systems, revealed FPSVDW+ES scoring yields significant early fold enrichment in the top 10% of ranked databases. For EGFR, top FPS poses are nicely accommodated in the molecular envelope defined by the reference in comparison with DCE, which yields distinct molecular weight bias toward larger molecules. Results from a representative virtual screen of ca. 1 million compounds additionally illustrate how ligands with footprints similar to a known inhibitor can readily be identified from within large commercially available databases. By providing an alternative way to rank ligand poses in a simple yet directed manner we anticipate that FPS scoring will be a useful tool for docking and structure-based design.

RAS Nanoclusters: Dynamic Signaling Platforms Amenable to Therapeutic Intervention.

RAS proteins are mutated in approximately 20% of all cancers and are generally associated with poor clinical outcomes. RAS proteins are localized to the plasma membrane and function as molecular switches, turned on by partners that receive extracellular mitogenic signals. In the on-state, they activate intracellular signal transduction cascades. Membrane-bound RAS molecules segregate into multimers, known as nanoclusters. These nanoclusters, held together through weak protein-protein and protein-lipid associations, are highly dynamic and respond to cellular input signals and fluctuations in the local lipid environment. Disruption of RAS nanoclusters results in downregulation of RAS-mediated mitogenic signaling. In this review, we discuss the propensity of RAS proteins to display clustering behavior and the interfaces that are associated with these assemblies. Strategies to therapeutically disrupt nanocluster formation or the stabilization of signaling incompetent RAS complexes are discussed.

A practical guide to large-scale docking.

Structure-based docking screens of large compound libraries have become common in early drug and probe discovery. As computer efficiency has improved and compound libraries have grown, the ability to screen hundreds of millions, and even billions, of compounds has become feasible for modest-sized computer clusters. This allows the rapid and cost-effective exploration and categorization of vast chemical space into a subset enriched with potential hits for a given target. To accomplish this goal at speed, approximations are used that result in undersampling of possible configurations and inaccurate predictions of absolute binding energies. Accordingly, it is important to establish controls, as are common in other fields, to enhance the likelihood of success in spite of these challenges. Here we outline best practices and control docking calculations that help evaluate docking parameters for a given target prior to undertaking a large-scale prospective screen, with exemplification in one particular target, the melatonin receptor, where following this procedure led to direct docking hits with activities in the subnanomolar range. Additional controls are suggested to ensure specific activity for experimentally validated hit compounds. These guidelines should be useful regardless of the docking software used. Docking software described in the outlined protocol (DOCK3.7) is made freely available for academic research to explore new hits for a range of targets.

Origins of resistance to the HIVgp41 viral entry inhibitor T20.

Peptide T20, which targets the HIV protein gp41, represents the first approved member of the class of HIV drugs known as membrane fusion inhibitors. However, mechanisms which lead to resistance through clinical use of T20 are not well-understood because the structure of the bound complex remains undetermined. In this report, an atomic-level model of a T20-gp41 complex embedded in an explicit DOPC membrane was constructed, and molecular dynamics simulations, followed by binding energy analysis (MM-GBSA method), were performed to delineate structural and energetic features that contribute to drug resistance. Per-residue binding footprints for T20 with wild-type gp41 reveal strong intermolecular van der Waals, Coulombic, and H-bond interactions in striking agreement with clinically observed resistance patterns. In addition, seven deleterious gp41 point mutations (L33Q, L33S, G36V, I37K, V38E, Q40H, and Q40K) were simulated, and all correctly exhibited decreases in the level of binding, including the fact that L33Q and Q40K are most detrimental. Six of the seven simulations yield good quantitative agreement (r(2) = 0.72; N = 6) with available experimental fold resistance data. Results from energy decomposition, heat map analysis, and differential (mutant minus wild-type) footprinting indicate the following. (1) Mutations disrupt intermolecular H-bonding and reduce the level of favorable contact with gp41 at M19. (2) Charged mutations (I37K, Q40K, and V38E) lead to significant Coulombic changes that weaken favorable van der Waals interactions. (3) Q40K is more detrimental than I37K because of interaction differences with a polar/charged patch on T20 in the initial (wild-type) state. (4) Resistance for L33S versus L33Q likely involves side chain packing differences in the final (mutated) state. A valuable finding of the work involves identification of favorable interactions among the C-terminal end of T20 (WNWF motif), residues on gp41 (including the fusion peptide), and headgroups in the adjacent membrane. The results suggest a complete T20 binding site would contribute to a stable complex, which could help to explain why prior studies, which employed truncated gp41 constructs, reported that C-terminal T20 residues may not interact with gp41. A hypothesis resulting from this study is that peptides could be designed to increase the level of favorable contact with both the membrane and gp41 which would lead to enhanced activity.

Vorinostat increases carboplatin and paclitaxel activity in non-small-cell lung cancer cells.

We observed a 53% response rate in non-small cell lung cancer (NSCLC) patients treated with vorinostat plus paclitaxel/carboplatin in a Phase I trial. Studies were undertaken to investigate the mechanism (s) underlying this activity. Growth inhibition was assessed in NSCLC cells by MTT assay after 72 hr of continuous drug exposure. Vorinostat (1 microM) inhibited growth by: 17% +/- 7% in A549, 28% +/- 6% in 128-88T, 39% +/- 8% in Calu1 and 41% +/- 7% in 201T cells. Vorinostat addition to carboplatin or paclitaxel led to significantly greater growth inhibition than chemotherapy alone in all 4 cell lines. Vorinostat (1 microM) synergistically increased the growth inhibitory effects of carboplatin/paclitaxel in 128-88T cells. When colony formation was measured after drug withdrawal, vorinostat significantly increased the effects of carboplatin but not paclitaxel. The % colony formation was control 100%; 1 microM vorinostat, 83% +/- 10%; 5 microM carboplatin, 41% +/- 11%; carboplatin/vorinostat, 8% +/- 4%; 2 nM paclitaxel, 53% +/- 11%; paclitaxel/vorinostat, 46% +/- 21%. In A549 and 128-88T, vorinostat potentiated carboplatin induction of gamma-H2AX (a DNA damage marker) and increased alpha-tubulin acetylation (a marker for stabilized mictrotubules). In A549, combination of vorinostat with paclitaxel resulted in a synergistic increase in alpha-tubulin acetylation, which reversed upon drug washout. We conclude that vorinostat interacts favorably with carboplatin and paclitaxel in NSCLC cells, which may explain the provocative response observed in our clinical trial. This likely involves a vorinostat-mediated irreversible increase in DNA damage in the case of carboplatin and a reversible increase in microtubule stability in the case of paclitaxel.