AlphaFold2 structures template ligand discovery

AlphaFold2 (AF2) and RosettaFold have greatly expanded the number of structures available for structure-based ligand discovery, even though retrospective studies have cast doubt on their direct usefulness for that goal. Here, we tested unrefined AF2 models prospectively, comparing experimental hit-rates and affinities from large library docking against AF2 models vs the same screens targeting experimental structures of the same receptors. In retrospective docking screens against the σ2 and the 5-HT2A receptors, the AF2 structures struggled to recapitulate ligands that we had previously found docking against the receptors' experimental structures, consistent with published results. Prospective large library docking against the AF2 models, however, yielded similar hit rates for both receptors versus docking against experimentally-derived structures; hundreds of molecules were prioritized and tested against each model and each structure of each receptor. The success of the AF2 models was achieved despite differences in orthosteric pocket residue conformations for both targets versus the experimental structures. Intriguingly, against the 5-HT2A receptor the most potent, subtype-selective agonists were discovered via docking against the AF2 model, not the experimental structure. To understand this from a molecular perspective, a cryoEM structure was determined for one of the more potent and selective ligands to emerge from docking against the AF2 model of the 5-HT2A receptor. Our findings suggest that AF2 models may sample conformations that are relevant for ligand discovery, much extending the domain of applicability of structure-based ligand discovery.

Structure-based discovery of CFTR potentiators and inhibitors

The cystic fibrosis transmembrane conductance regulator (CFTR) is a crucial ion channel whose loss of function leads to cystic fibrosis, while its hyperactivation leads to secretory diarrhea. Small molecules that improve CFTR folding (correctors) or function (potentiators) are clinically available. However, the only potentiator, ivacaftor, has suboptimal pharmacokinetics and inhibitors have yet to be clinically developed. Here we combine molecular docking, electrophysiology, cryo-EM, and medicinal chemistry to identify novel CFTR modulators. We docked ∼155 million molecules into the potentiator site on CFTR, synthesized 53 test ligands, and used structure-based optimization to identify candidate modulators. This approach uncovered novel mid-nanomolar potentiators as well as inhibitors that bind to the same allosteric site. These molecules represent potential leads for the development of more effective drugs for cystic fibrosis and secretory diarrhea, demonstrating the feasibility of large-scale docking for ion channel drug discovery.

Small vs. Large Library Docking for Positive Allosteric Modulators of the Calcium Sensing Receptor

Drugs acting as positive allosteric modulators (PAMs) to enhance the activation of the calcium sensing receptor (CaSR) and to suppress parathyroid hormone (PTH) secretion can treat hyperparathyroidism but suffer from side effects including hypocalcemia and arrhythmias. Seeking new CaSR modulators, we docked libraries of 2.7 million and 1.2 billion molecules against transforming pockets in the active-state receptor dimer structure. Consistent with simulations suggesting that docking improves with library size, billion-molecule docking found new PAMs with a hit rate that was 2.7-fold higher than the million-molecule library and with hits up to 37-fold more potent. Structure-based optimization of ligands from both campaigns led to nanomolar leads, one of which was advanced to animal testing. This PAM displays 100-fold the potency of the standard of care, cinacalcet, in ex vivo organ assays, and reduces serum PTH levels in mice by up to 80% without the hypocalcemia typical of CaSR drugs. Cryo-EM structures with the new PAMs show that they induce residue rearrangements in the binding pockets and promote CaSR dimer conformations that are closer to the G-protein coupled state compared to established drugs. These findings highlight the promise of large library docking for therapeutic leads, especially when combined with experimental structure determination and mechanism.

Large library docking for cannabinoid-1 receptor agonists with reduced side effects

Large library docking can reveal unexpected chemotypes that complement the structures of biological targets. Seeking new agonists for the cannabinoid-1 receptor (CB1R), we docked 74 million tangible molecules, prioritizing 46 high ranking ones for de novo synthesis and testing. Nine were active by radioligand competition, a 20% hit-rate. Structure-based optimization of one of the most potent of these (Ki = 0.7 uM) led to '4042, a 1.9 nM ligand and a full CB1R agonist. A cryo-EM structure of the purified enantiomer of '4042 ('1350) in complex with CB1R-Gi1 confirmed its docked pose. The new agonist was strongly analgesic, with generally a 5-10-fold therapeutic window over sedation and catalepsy and no observable conditioned place preference. These findings suggest that new cannabinoid chemotypes may disentangle characteristic cannabinoid side-effects from their analgesia, supporting the further development of cannabinoids as pain therapeutics.

DockOpt: A Tool for Automatic Optimization of Docking Models

Molecular docking is a widely used technique for leveraging protein structure for ligand discovery, but it remains difficult to utilize due to limitations that have not been adequately addressed. Despite some progress toward automation, docking still requires expert guidance, hindering its adoption by a broader range of investigators. To make docking more accessible, we developed a new utility called DockOpt, which automates the creation, evaluation, and optimization of docking models prior to their deployment in large-scale prospective screens. DockOpt outperforms our previous automated pipeline across all 43 targets in the DUDE-Z benchmark data set, and the generated models for 84% of targets demonstrate sufficient enrichment to warrant their use in prospective screens, with normalized LogAUC values of at least 15%. DockOpt is available as part of the Python package Pydock3 included in the UCSF DOCK 3.8 distribution, which is available for free to academic researchers at https://dock.compbio.ucsf.edu and free for everyone upon registration at https://tldr.docking.org.

Docking for Molecules That Bind in a Symmetric Stack with SymDOCK

Discovering ligands for amyloid fibrils, such as those formed by the tau protein, is an area of great current interest. In recent structures, ligands bind in stacks in the tau fibrils to reflect the rotational and translational symmetry of the fibril itself; in these structures, the ligands make few interactions with the protein but interact extensively with each other. To exploit this symmetry and stacking, we developed SymDOCK, a method to dock molecules that follow the protein's symmetry. For each prospective ligand pose, we apply the symmetry operation of the fibril to generate a self-interacting and fibril-interacting stack, checking that doing so will not cause a clash between the original molecule and its image. Absent a clash, we retain that pose and add the ligand-ligand van der Waals energy to the ligand's docking score (here using DOCK3.8). We can check these geometries and energies using an implementation of ANI, a neural-network-based quantum-mechanical evaluation of the ligand stacking energies. In retrospective calculations, symmetry docking can reproduce the poses of three tau PET tracers whose structures have been determined. More convincingly, in a prospective study, SymDOCK predicted the structure of the PET tracer MK-6240 bound in a symmetrical stack to AD PHF tau before that structure was determined; the docked pose was used to determine how MK-6240 fit the cryo-EM density. In proof-of-concept studies, SymDOCK enriched known ligands over property-matched decoys in retrospective screens without sacrificing docking speed and can address large library screens that seek new symmetrical stackers. Future applications of this approach will be considered.

Docking for EP4R antagonists active against inflammatory pain

The lipid prostaglandin E2 (PGE2) mediates inflammatory pain by activating G protein-coupled receptors, including the prostaglandin E2 receptor 4 (EP4R). Nonsteroidal anti-inflammatory drugs (NSAIDs) reduce nociception by inhibiting prostaglandin synthesis, however, the disruption of upstream prostanoid biosynthesis can lead to pleiotropic effects including gastrointestinal bleeding and cardiac complications. In contrast, by acting downstream, EP4R antagonists may act specifically as anti-inflammatory agents and, to date, no selective EP4R antagonists have been approved for human use. In this work, seeking to diversify EP4R antagonist scaffolds, we computationally dock over 400 million compounds against an EP4R crystal structure and experimentally validate 71 highly ranked, de novo synthesized molecules. Further, we show how structure-based optimization of initial docking hits identifies a potent and selective antagonist with 16 nanomolar potency. Finally, we demonstrate favorable pharmacokinetics for the discovered compound as well as anti-allodynic and anti-inflammatory activity in several preclinical pain models in mice.

Docking for molecules that bind in a symmetric stack with SymDOCK

Discovering ligands for amyloid fibrils, such as those formed by the tau protein, is an area of much current interest. In recent structures, ligands bind in stacks in the tau fibrils to reflect the rotational and translational symmetry of the fibril itself; in these structures the ligands make few interactions with the protein but interact extensively with each other. To exploit this symmetry and stacking, we developed SymDOCK, a method to dock molecules that follow the protein's symmetry. For each prospective ligand pose, we apply the symmetry operation of the fibril to generate a self-interacting and fibril-interacting stack, checking that doing so will not cause a clash between the original molecule and its image. Absent a clash, we retain that pose and add the ligand-ligand van der Waals energy to the ligand's docking score (here using DOCK3.8). We can check these geometries and energies using an implementation of ANI, a neural network-based quantum-mechanical evaluation of the ligand stacking energies. In retrospective calculations, symmetry docking can reproduce the poses of three tau PET tracers whose structures have been determined. More convincingly, in a prospective study SymDOCK predicted the structure of the PET tracer MK-6240 bound in a symmetrical stack to AD PHF tau before that structure was determined; the docked pose was used to determine how MK-6240 fit the cryo-EM density. In proof-of-concept studies, SymDOCK enriched known ligands over property-matched decoys in retrospective screens without sacrificing docking speed, and can address large library screens that seek new symmetrical stackers. Future applications of this approach will be considered.

Large library docking for novel SARS-CoV-2 main protease non-covalent and covalent inhibitors

Antiviral therapeutics to treat SARS-CoV-2 are needed to diminish the morbidity of the ongoing COVID-19 pandemic. A well-precedented drug target is the main viral protease (MPro ), which is targeted by an approved drug and by several investigational drugs. Emerging viral resistance has made new inhibitor chemotypes more pressing. Adopting a structure-based approach, we docked 1.2 billion non-covalent lead-like molecules and a new library of 6.5 million electrophiles against the enzyme structure. From these, 29 non-covalent and 11 covalent inhibitors were identified in 37 series, the most potent having an IC50 of 29 and 20 μM, respectively. Several series were optimized, resulting in low micromolar inhibitors. Subsequent crystallography confirmed the docking predicted binding modes and may template further optimization. While the new chemotypes may aid further optimization of MPro inhibitors for SARS-CoV-2, the modest success rate also reveals weaknesses in our approach for challenging targets like MPro versus other targets where it has been more successful, and versus other structure-based techniques against MPro itself.

Structure-Based Discovery of Inhibitors of the SARS-CoV-2 Nsp14 N7-Methyltransferase

An under-explored target for SARS-CoV-2 is the S-adenosyl methionine (SAM)-dependent methyltransferase Nsp14, which methylates the N7-guanosine of viral RNA at the 5'-end, allowing the virus to evade host immune response. We sought new Nsp14 inhibitors with three large library docking strategies. First, up to 1.1 billion lead-like molecules were docked against the enzyme's SAM site, leading to three inhibitors with IC50 values from 6 to 50 μM. Second, docking a library of 16 million fragments revealed 9 new inhibitors with IC50 values from 12 to 341 μM. Third, docking a library of 25 million electrophiles to covalently modify Cys387 revealed 7 inhibitors with IC50 values from 3.5 to 39 μM. Overall, 32 inhibitors encompassing 11 chemotypes had IC50 values < 50 μM and 5 inhibitors in 4 chemotypes had IC50 values < 10 μM. These molecules are among the first non-SAM-like inhibitors of Nsp14, providing starting points for future optimization.

Modeling the expansion of virtual screening libraries

Recently, 'tangible' virtual libraries have made billions of molecules readily available. Prioritizing these molecules for synthesis and testing demands computational approaches, such as docking. Their success may depend on library diversity, their similarity to bio-like molecules and how receptor fit and artifacts change with library size. We compared a library of 3 million 'in-stock' molecules with billion-plus tangible libraries. The bias toward bio-like molecules in the tangible library decreases 19,000-fold versus those 'in-stock'. Similarly, thousands of high-ranking molecules, including experimental actives, from five ultra-large-library docking campaigns are also dissimilar to bio-like molecules. Meanwhile, better-fitting molecules are found as the library grows, with the score improving log-linearly with library size. Finally, as library size increases, so too do rare molecules that rank artifactually well. Although the nature of these artifacts changes from target to target, the expectation of their occurrence does not, and simple strategies can minimize their impact.

Structure-based discovery of conformationally selective inhibitors of the serotonin transporter

The serotonin transporter (SERT) removes synaptic serotonin and is the target of anti-depressant drugs. SERT adopts three conformations: outward-open, occluded, and inward-open. All known inhibitors target the outward-open state except ibogaine, which has unusual anti-depressant and substance-withdrawal effects, and stabilizes the inward-open conformation. Unfortunately, ibogaine's promiscuity and cardiotoxicity limit the understanding of inward-open state ligands. We docked over 200 million small molecules against the inward-open state of the SERT. Thirty-six top-ranking compounds were synthesized, and thirteen inhibited; further structure-based optimization led to the selection of two potent (low nanomolar) inhibitors. These stabilized an outward-closed state of the SERT with little activity against common off-targets. A cryo-EM structure of one of these bound to the SERT confirmed the predicted geometry. In mouse behavioral assays, both compounds had anxiolytic- and anti-depressant-like activity, with potencies up to 200-fold better than fluoxetine (Prozac), and one substantially reversed morphine withdrawal effects.

Large-Scale Docking in the Cloud

Molecular docking is a pragmatic approach to exploit protein structures for new ligand discovery, but the growing size of available chemical space is increasingly challenging to screen on in-house computer clusters. We have therefore developed AWS-DOCK, a protocol for running UCSF DOCK in the AWS cloud. Our approach leverages the low cost and scalability of cloud resources combined with a low-molecule-cost docking engine to screen billions of molecules efficiently. We benchmarked our system by screening 50 million HAC 22 molecules against the DRD4 receptor with an average CPU time of around 1 s per molecule. We saw up to 3-fold variations in cost between AWS availability zones. Docking 4.5 billion lead-like molecules, a 7 week calculation on our 1000-core lab cluster, runs in about a week depending on accessible CPUs, in AWS for around $25,000, less than the cost of two new nodes. The cloud docking protocol is described in easy-to-follow steps and may be sufficiently general to be used for other docking programs. All the tools to enable AWS-DOCK are available free to everyone, while DOCK 3.8 is free for academic research.

ZINC-22─A Free Multi-Billion-Scale Database of Tangible Compounds for Ligand Discovery

Purchasable chemical space has grown rapidly into the tens of billions of molecules, providing unprecedented opportunities for ligand discovery but straining the tools that might exploit these molecules at scale. We have therefore developed ZINC-22, a database of commercially accessible small molecules derived from multi-billion-scale make-on-demand libraries. The new database and tools enable analog searching in this vast new space via a facile GUI, CartBlanche, drawing on similarity methods that scale sublinearly in the number of molecules. The new library also uses data organization methods, enabling rapid lookup of molecules and their physical properties, including conformations, partial atomic charges, c Log P values, and solvation energies, all crucial for molecule docking, which had become slow with older database organizations in previous versions of ZINC. As the libraries have continued to grow, we have been interested in finding whether molecular diversity has suffered, for instance, because certain scaffolds have come to dominate via easy analoging. This has not occurred thus far, and chemical diversity continues to grow with database size, with a log increase in Bemis-Murcko scaffolds for every two-log unit increase in database size. Most new scaffolds come from compounds with the highest heavy atom count. Finally, we consider the implications for databases like ZINC as the libraries grow toward and beyond the trillion-molecule range. ZINC is freely available to everyone and may be accessed at cartblanche22.docking.org, via Globus, and in the Amazon AWS and Oracle OCI clouds.

Iterative computational design and crystallographic screening identifies potent inhibitors targeting the Nsp3 macrodomain of SARS-CoV-2

The nonstructural protein 3 (NSP3) of the severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) contains a conserved macrodomain enzyme (Mac1) that is critical for pathogenesis and lethality. While small-molecule inhibitors of Mac1 have great therapeutic potential, at the outset of the COVID-19 pandemic, there were no well-validated inhibitors for this protein nor, indeed, the macrodomain enzyme family, making this target a pharmacological orphan. Here, we report the structure-based discovery and development of several different chemical scaffolds exhibiting low- to sub-micromolar affinity for Mac1 through iterations of computer-aided design, structural characterization by ultra-high-resolution protein crystallography, and binding evaluation. Potent scaffolds were designed with in silico fragment linkage and by ultra-large library docking of over 450 million molecules. Both techniques leverage the computational exploration of tangible chemical space and are applicable to other pharmacological orphans. Overall, 160 ligands in 119 different scaffolds were discovered, and 153 Mac1-ligand complex crystal structures were determined, typically to 1 Å resolution or better. Our analyses discovered selective and cell-permeable molecules, unexpected ligand-mediated conformational changes within the active site, and key inhibitor motifs that will template future drug development against Mac1.

Bespoke library docking for 5-HT2A receptor agonists with antidepressant activity

There is considerable interest in screening ultralarge chemical libraries for ligand discovery, both empirically and computationally1-4. Efforts have focused on readily synthesizable molecules, inevitably leaving many chemotypes unexplored. Here we investigate structure-based docking of a bespoke virtual library of tetrahydropyridines-a scaffold that is poorly sampled by a general billion-molecule virtual library but is well suited to many aminergic G-protein-coupled receptors. Using three inputs, each with diverse available derivatives, a one pot C-H alkenylation, electrocyclization and reduction provides the tetrahydropyridine core with up to six sites of derivatization5-7. Docking a virtual library of 75 million tetrahydropyridines against a model of the serotonin 5-HT2A receptor (5-HT2AR) led to the synthesis and testing of 17 initial molecules. Four of these molecules had low-micromolar activities against either the 5-HT2A or the 5-HT2B receptors. Structure-based optimization led to the 5-HT2AR agonists (R)-69 and (R)-70, with half-maximal effective concentration values of 41 nM and 110 nM, respectively, and unusual signalling kinetics that differ from psychedelic 5-HT2AR agonists. Cryo-electron microscopy structural analysis confirmed the predicted binding mode to 5-HT2AR. The favourable physical properties of these new agonists conferred high brain permeability, enabling mouse behavioural assays. Notably, neither had psychedelic activity, in contrast to classic 5-HT2AR agonists, whereas both had potent antidepressant activity in mouse models and had the same efficacy as antidepressants such as fluoxetine at as low as 1/40th of the dose. Prospects for using bespoke virtual libraries to sample pharmacologically relevant chemical space will be considered.

Structure-based discovery of nonopioid analgesics acting through the α2A-adrenergic receptor

Because nonopioid analgesics are much sought after, we computationally docked more than 301 million virtual molecules against a validated pain target, the α2A-adrenergic receptor (α2AAR), seeking new α2AAR agonists chemotypes that lack the sedation conferred by known α2AAR drugs, such as dexmedetomidine. We identified 17 ligands with potencies as low as 12 nanomolar, many with partial agonism and preferential Gi and Go signaling. Experimental structures of α2AAR complexed with two of these agonists confirmed the docking predictions and templated further optimization. Several compounds, including the initial docking hit '9087 [mean effective concentration (EC50) of 52 nanomolar] and two analogs, '7075 and PS75 (EC50 4.1 and 4.8 nanomolar), exerted on-target analgesic activity in multiple in vivo pain models without sedation. These newly discovered agonists are interesting as therapeutic leads that lack the liabilities of opioids and the sedation of dexmedetomidine.

Iterative computational design and crystallographic screening identifies potent inhibitors targeting the Nsp3 Macrodomain of SARS-CoV-2

The nonstructural protein 3 (NSP3) of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) contains a conserved macrodomain enzyme (Mac1) that is critical for pathogenesis and lethality. While small molecule inhibitors of Mac1 have great therapeutic potential, at the outset of the COVID-19 pandemic there were no well-validated inhibitors for this protein nor, indeed, the macrodomain enzyme family, making this target a pharmacological orphan. Here, we report the structure-based discovery and development of several different chemical scaffolds exhibiting low- to sub-micromolar affinity for Mac1 through iterations of computer-aided design, structural characterization by ultra-high resolution protein crystallography, and binding evaluation. Potent scaffolds were designed with in silico fragment linkage and by ultra-large library docking of over 450 million molecules. Both techniques leverage the computational exploration of tangible chemical space and are applicable to other pharmacological orphans. Overall, 160 ligands in 119 different scaffolds were discovered, and 152 Mac1-ligand complex crystal structures were determined, typically to 1 Å resolution or better. Our analyses discovered selective and cell-permeable molecules, unexpected ligand-mediated protein dynamics within the active site, and key inhibitor motifs that will template future drug development against Mac1.

Ligand Strain Energy in Large Library Docking

While small molecule internal strain is crucial to molecular docking, using it in evaluating ligand scores has remained elusive. Here, we investigate a technique that calculates strain using relative torsional populations in the Cambridge Structural Database, enabling fast precalculation of these energies. In retrospective studies of large docking screens of the dopamine D4 receptor and of AmpC β-lactamase, where close to 600 docking hits were tested experimentally, including such strain energies improved hit rates by preferentially reducing the ranks of strained high-scoring decoy molecules. In a 40-target subset of the DUD-E benchmark, we found two thresholds that usefully distinguished between ligands and decoys: one based on the total strain energy of the small molecules and another based on the maximum strain allowed for any given torsion within them. Using these criteria, about 75% of the benchmark targets had improved enrichment after strain filtering. Relying on precalculated population distributions, this approach is rapid, taking less than 0.04 s to evaluate a conformation on a standard core, making it pragmatic for precalculating strain in even ultralarge libraries. Since it is scoring function agnostic, it may be useful to multiple docking approaches; it is openly available at http://tldr.docking.org.

Fragment binding to the Nsp3 macrodomain of SARS-CoV-2 identified through crystallographic screening and computational docking

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) macrodomain within the nonstructural protein 3 counteracts host-mediated antiviral adenosine diphosphate-ribosylation signaling. This enzyme is a promising antiviral target because catalytic mutations render viruses nonpathogenic. Here, we report a massive crystallographic screening and computational docking effort, identifying new chemical matter primarily targeting the active site of the macrodomain. Crystallographic screening of 2533 diverse fragments resulted in 214 unique macrodomain-binders. An additional 60 molecules were selected from docking more than 20 million fragments, of which 20 were crystallographically confirmed. X-ray data collection to ultra-high resolution and at physiological temperature enabled assessment of the conformational heterogeneity around the active site. Several fragment hits were confirmed by solution binding using three biophysical techniques (differential scanning fluorimetry, homogeneous time-resolved fluorescence, and isothermal titration calorimetry). The 234 fragment structures explore a wide range of chemotypes and provide starting points for development of potent SARS-CoV-2 macrodomain inhibitors.

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.

ZINC20-A Free Ultralarge-Scale Chemical Database for Ligand Discovery

Identifying and purchasing new small molecules to test in biological assays are enabling for ligand discovery, but as purchasable chemical space continues to grow into the tens of billions based on inexpensive make-on-demand compounds, simply searching this space becomes a major challenge. We have therefore developed ZINC20, a new version of ZINC with two major new features: billions of new molecules and new methods to search them. As a fully enumerated database, ZINC can be searched precisely using explicit atomic-level graph-based methods, such as SmallWorld for similarity and Arthor for pattern and substructure search, as well as 3D methods such as docking. Analysis of the new make-on-demand compound sets by these and related tools reveals startling features. For instance, over 97% of the core Bemis-Murcko scaffolds in make-on-demand libraries are unavailable from "in-stock" collections. Correspondingly, the number of new Bemis-Murcko scaffolds is rising almost as a linear fraction of the elaborated molecules. Thus, an 88-fold increase in the number of molecules in the make-on-demand versus the in-stock sets is built upon a 16-fold increase in the number of Bemis-Murcko scaffolds. The make-on-demand library is also more structurally diverse than physical libraries, with a massive increase in disc- and sphere-like shaped molecules. The new system is freely available at zinc20.docking.org.

The activities of drug inactive ingredients on biological targets

Excipients, considered "inactive ingredients," are a major component of formulated drugs and play key roles in their pharmacokinetics. Despite their pervasiveness, whether they are active on any targets has not been systematically explored. We computed the likelihood that approved excipients would bind to molecular targets. Testing in vitro revealed 25 excipient activities, ranging from low-nanomolar to high-micromolar concentration. Another 109 activities were identified by testing against clinical safety targets. In cellular models, five excipients had fingerprints predictive of system-level toxicity. Exposures of seven excipients were investigated, and in certain populations, two of these may reach levels of in vitro target potency, including brain and gut exposure of thimerosal and its major metabolite, which had dopamine D3 receptor dissociation constant K _d values of 320 and 210 nM, respectively. Although most excipients deserve their status as inert, many approved excipients may directly modulate physiologically relevant targets.

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.

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.

Predicted Biological Activity of Purchasable Chemical Space

Whereas 400 million distinct compounds are now purchasable within the span of a few weeks, the biological activities of most are unknown. To facilitate access to new chemistry for biology, we have combined the Similarity Ensemble Approach (SEA) with the maximum Tanimoto similarity to the nearest bioactive to predict activity for every commercially available molecule in ZINC. This method, which we label SEA+TC, outperforms both SEA and a naïve-Bayesian classifier via predictive performance on a 5-fold cross-validation of ChEMBL’s bioactivity data set (version 21). Using this method, predictions for over 40% of compounds (>160 million) have either high significance (pSEA ≥ 40), high similarity (ECFP4MaxTc ≥ 0.4), or both, for one or more of 1382 targets well described by ligands in the literature. Using a further 1347 less-well-described targets, we predict activities for an additional 11 million compounds. To gauge whether these predictions are sensible, we investigate 75 predictions for 50 drugs lacking a binding affinity annotation in ChEMBL. The 535 million predictions for over 171 million compounds at 2629 targets are linked to purchasing information and evidence to support each prediction and are freely available via https://zinc15.docking.org and https://files.docking.org.

A Molecular Basis for Innovation in Drug Excipients

Excipients are ubiquitous in drug formulation, ensuring that active ingredient drugs are properly released on dosing, retain their properties over time, and are palatable, among other roles. Despite their crucial roles, surprisingly little is known about their systemic availability and activities on molecular targets. Here we review key excipient properties, introduce a public-accessible database that enumerates and categorizes them, and sketch a strategy for exploring their possible direct actions on molecular targets.

Identification of Novel Smoothened Ligands Using Structure-Based Docking

The seven transmembrane protein Smoothened is required for Hedgehog signaling during embryonic development and adult tissue homeostasis. Inappropriate activation of the Hedgehog signalling pathway leads to cancers such as basal cell carcinoma and medulloblastoma, and Smoothened inhibitors are now available clinically to treat these diseases. However, resistance to these inhibitors rapidly develops thereby limiting their efficacy. The determination of Smoothened crystal structures enables structure-based discovery of new ligands with new chemotypes that will be critical to combat resistance. In this study, we docked 3.2 million available, lead-like molecules against Smoothened, looking for those with high physical complementarity to its structure; this represents the first such campaign against the class Frizzled G-protein coupled receptor family. Twenty-one high-ranking compounds were selected for experimental testing, and four, representing three different chemotypes, were identified to antagonize Smoothened with IC₅₀ values better than 50 μM. A screen for analogs revealed another six molecules, with IC₅₀ values in the low micromolar range. Importantly, one of the most active of the new antagonists continued to be efficacious at the D473H mutant of Smoothened, which confers clinical resistance to the antagonist vismodegib in cancer treatment.

Docking Screens for Novel Ligands Conferring New Biology

It is now plausible to dock libraries of 10 million molecules against targets over several days or weeks. When the molecules screened are commercially available, they may be rapidly tested to find new leads. Although docking retains important liabilities (it cannot calculate affinities accurately nor even reliably rank order high-scoring molecules), it can often can distinguish likely from unlikely ligands, often with hit rates above 10%. Here we summarize the improvements in libraries, target quality, and methods that have supported these advances, and the open access resources that make docking accessible. Recent docking screens for new ligands are sketched, as are the binding, crystallographic, and in vivo assays that support them. Like any technique, controls are crucial, and key experimental ones are reviewed. With such controls, docking campaigns can find ligands with new chemotypes, often revealing the new biology that may be docking's greatest impact over the next few years.

ZINC 15--Ligand Discovery for Everyone

Many questions about the biological activity and availability of small molecules remain inaccessible to investigators who could most benefit from their answers. To narrow the gap between chemoinformatics and biology, we have developed a suite of ligand annotation, purchasability, target, and biology association tools, incorporated into ZINC and meant for investigators who are not computer specialists. The new version contains over 120 million purchasable “drug-like” compounds – effectively all organic molecules that are for sale – a quarter of which are available for immediate delivery. ZINC connects purchasable compounds to high-value ones such as metabolites, drugs, natural products, and annotated compounds from the literature. Compounds may be accessed by the genes for which they are annotated as well as the major and minor target classes to which those genes belong. It offers new analysis tools that are easy for nonspecialists yet with few limitations for experts. ZINC retains its original 3D roots – all molecules are available in biologically relevant, ready-to-dock formats. ZINC is freely available at http://zinc15.docking.org.

An Aggregation Advisor for Ligand Discovery

Colloidal aggregation of organic molecules is the dominant mechanism for artifactual inhibition of proteins, and controls against it are widely deployed. Notwithstanding an increasingly detailed understanding of this phenomenon, a method to reliably predict aggregation has remained elusive. Correspondingly, active molecules that act via aggregation continue to be found in early discovery campaigns and remain common in the literature. Over the past decade, over 12 thousand aggregating organic molecules have been identified, potentially enabling a precedent-based approach to match known aggregators with new molecules that may be expected to aggregate and lead to artifacts. We investigate an approach that uses lipophilicity, affinity, and similarity to known aggregators to advise on the likelihood that a candidate compound is an aggregator. In prospective experimental testing, five of seven new molecules with Tanimoto coefficients (Tc's) between 0.95 and 0.99 to known aggregators aggregated at relevant concentrations. Ten of 19 with Tc's between 0.94 and 0.90 and three of seven with Tc's between 0.89 and 0.85 also aggregated. Another three of the predicted compounds aggregated at higher concentrations. This method finds that 61 827 or 5.1% of the ligands acting in the 0.1 to 10 μM range in the medicinal chemistry literature are at least 85% similar to a known aggregator with these physical properties and may aggregate at relevant concentrations. Intriguingly, only 0.73% of all drug-like commercially available compounds resemble the known aggregators, suggesting that colloidal aggregators are enriched in the literature. As a percentage of the literature, aggregator-like compounds have increased 9-fold since 1995, partly reflecting the advent of high-throughput and virtual screens against molecular targets. Emerging from this study is an aggregator advisor database and tool ( http://advisor.bkslab.org ), free to the community, that may help distinguish between fruitful and artifactual screening hits acting by this mechanism.

Increasing chemical space coverage by combining empirical and computational fragment screens

Most libraries for fragment-based drug discovery are restricted to 1,000-10,000 compounds, but over 500,000 fragments are commercially available and potentially accessible by virtual screening. Whether this larger set would increase chemotype coverage, and whether a computational screen can pragmatically prioritize them, is debated. To investigate this question, a 1281-fragment library was screened by nuclear magnetic resonance (NMR) against AmpC β-lactamase, and hits were confirmed by surface plasmon resonance (SPR). Nine hits with novel chemotypes were confirmed biochemically with KI values from 0.2 to low mM. We also computationally docked 290,000 purchasable fragments with chemotypes unrepresented in the empirical library, finding 10 that had KI values from 0.03 to low mM. Though less novel than those discovered by NMR, the docking-derived fragments filled chemotype holes from the empirical library. Crystal structures of nine of the fragments in complex with AmpC β-lactamase revealed new binding sites and explained the relatively high affinity of the docking-derived fragments. The existence of chemotype holes is likely a general feature of fragment libraries, as calculation suggests that to represent the fragment substructures of even known biogenic molecules would demand a library of minimally over 32,000 fragments. Combining computational and empirical fragment screens enables the discovery of unexpected chemotypes, here by the NMR screen, while capturing chemotypes missing from the empirical library and tailored to the target, with little extra cost in resources.

Novel compounds lowering the cellular isoform of the human prion protein in cultured human cells

PURPOSE

Previous studies showed that lowering PrP(C) concomitantly reduced PrP(Sc) in the brains of mice inoculated with prions. We aimed to develop assays that measure PrP(C) on the surface of human T98G glioblastoma and IMR32 neuroblastoma cells. Using these assays, we sought to identify chemical hits, confirmed hits, and scaffolds that potently lowered PrP(C) levels in human brains cells, without lethality, and that could achieve drug concentrations in the brain after oral or intraperitoneal dosing in mice.

METHODS

We utilized HTS ELISA assays to identify small molecules that lower PrP(C) levels by ≥30% on the cell surface of human glioblastoma (T98G) and neuroblastoma (IMR32) cells.

RESULTS

From 44,578 diverse chemical compounds tested, 138 hits were identified by single point confirmation (SPC) representing 7 chemical scaffolds in T98G cells, and 114 SPC hits representing 6 scaffolds found in IMR32 cells. When the confirmed SPC hits were combined with structurally related analogs, >300 compounds (representing 6 distinct chemical scaffolds) were tested for dose-response (EC₅₀) in both cell lines, only studies in T98G cells identified compounds that reduced PrP(C) without killing the cells. EC₅₀ values from 32 hits ranged from 65 nM to 4.1 μM. Twenty-eight were evaluated in vivo in pharmacokinetic studies after a single 10 mg/kg oral or intraperitoneal dose in mice. Our results showed brain concentrations as high as 16.2 μM, but only after intraperitoneal dosing.

CONCLUSIONS

Our studies identified leads for future studies to determine which compounds might lower PrP(C) levels in rodent brain, and provide the basis of a therapeutic for fatal disorders caused by PrP prions.

Ligand pose and orientational sampling in molecular docking

Molecular docking remains an important tool for structure-based screening to find new ligands and chemical probes. As docking ambitions grow to include new scoring function terms, and to address ever more targets, the reliability and extendability of the orientation sampling, and the throughput of the method, become pressing. Here we explore sampling techniques that eliminate stochastic behavior in DOCK3.6, allowing us to optimize the method for regularly variable sampling of orientations. This also enabled a focused effort to optimize the code for efficiency, with a three-fold increase in the speed of the program. This, in turn, facilitated extensive testing of the method on the 102 targets, 22,805 ligands and 1,411,214 decoys of the Directory of Useful Decoys - Enhanced (DUD-E) benchmarking set, at multiple levels of sampling. Encouragingly, we observe that as sampling increases from 50 to 500 to 2000 to 5000 to 20000 molecular orientations in the binding site (and so from about 1×10¹⁰ to 4×10¹⁰ to 1×10¹¹ to 2×10¹¹ to 5×10¹¹ mean atoms scored per target, since multiple conformations are sampled per orientation), the enrichment of ligands over decoys monotonically increases for most DUD-E targets. Meanwhile, including internal electrostatics in the evaluation ligand conformational energies, and restricting aromatic hydroxyls to low energy rotamers, further improved enrichment values. Several of the strategies used here to improve the efficiency of the code are broadly applicable in the field.

Structure-based discovery of antagonists of nuclear receptor LRH-1

Liver receptor homolog 1 (nuclear receptor LRH-1, NR5A2) is an essential regulator of gene transcription, critical for maintenance of cell pluripotency in early development and imperative for the proper functions of the liver, pancreas, and intestines during the adult life. Although physiological hormones of LRH-1 have not yet been identified, crystallographic and biochemical studies demonstrated that LRH-1 could bind regulatory ligands and suggested phosphatidylinositols as potential hormone candidates for this receptor. No synthetic antagonists of LRH-1 are known to date. Here, we identify the first small molecule antagonists of LRH-1 activity. Our search for LRH-1 modulators was empowered by screening of 5.2 million commercially available compounds via molecular docking followed by verification of the top-ranked molecules using in vitro direct binding and transcriptional assays. Experimental evaluation of the predicted ligands identified two compounds that inhibit the transcriptional activity of LRH-1 and diminish the expression of the receptor's target genes. Among the affected transcriptional targets are co-repressor SHP (small heterodimer partner) as well as cyclin E1 (CCNE1) and G0S2 genes that are known to regulate cell growth and proliferation. Treatments of human pancreatic (AsPC-1), colon (HT29), and breast adenocarcinoma cells T47D and MDA-MB-468 with the LRH-1 antagonists resulted in the receptor-mediated inhibition of cancer cell proliferation. Our data suggest that specific antagonists of LRH-1 could be used as specific molecular probes for elucidating the roles of the receptor in different types of malignancies.

Directory of useful decoys, enhanced (DUD-E): better ligands and decoys for better benchmarking

A key metric to assess molecular docking remains ligand enrichment against challenging decoys. Whereas the directory of useful decoys (DUD) has been widely used, clear areas for optimization have emerged. Here we describe an improved benchmarking set that includes more diverse targets such as GPCRs and ion channels, totaling 102 proteins with 22886 clustered ligands drawn from ChEMBL, each with 50 property-matched decoys drawn from ZINC. To ensure chemotype diversity, we cluster each target’s ligands by their Bemis–Murcko atomic frameworks. We add net charge to the matched physicochemical properties and include only the most dissimilar decoys, by topology, from the ligands. An online automated tool (http://decoys.docking.org) generates these improved matched decoys for user-supplied ligands. We test this data set by docking all 102 targets, using the results to improve the balance between ligand desolvation and electrostatics in DOCK 3.6. The complete DUD-E benchmarking set is freely available at http://dude.docking.org.

ZINC: a free tool to discover chemistry for biology

ZINC is a free public resource for ligand discovery. The database contains over twenty million commercially available molecules in biologically relevant representations that may be downloaded in popular ready-to-dock formats and subsets. The Web site also enables searches by structure, biological activity, physical property, vendor, catalog number, name, and CAS number. Small custom subsets may be created, edited, shared, docked, downloaded, and conveyed to a vendor for purchase. The database is maintained and curated for a high purchasing success rate and is freely available at zinc.docking.org.

Virtual ligand screening against comparative protein structure models

Virtual ligand screening uses computation to discover new ligands of a protein by screening one or more of its structural models against a database of potential ligands. Comparative protein structure modeling extends the applicability of virtual screening beyond the atomic structures determined by X-ray crystallography or NMR spectroscopy. Here, we describe an integrated modeling and docking protocol, combining comparative modeling by MODELLER and virtual ligand screening by DOCK.

Chemical informatics and target identification in a zebrafish phenotypic screen

Target identification is a core challenge in chemical genetics. Here we use chemical similarity to predict computationally the targets of 586 compounds active in a zebrafish behavioral assay. Of 20 predictions tested, 11 had activities ranging from 1 to 10,000nM on the predicted targets. The role of two of these targets was tested in the original zebrafish phenotype. Prediction of targets from chemotype is rapid and may be generally applicable.

Statistical potential for modeling and ranking of protein-ligand interactions

Applications in structural biology and medicinal chemistry require protein-ligand scoring functions for two distinct tasks: (i) ranking different poses of a small molecule in a protein binding site and (ii) ranking different small molecules by their complementarity to a protein site. Using probability theory, we developed two atomic distance-dependent statistical scoring functions: PoseScore was optimized for recognizing native binding geometries of ligands from other poses and RankScore was optimized for distinguishing ligands from nonbinding molecules. Both scores are based on a set of 8,885 crystallographic structures of protein-ligand complexes but differ in the values of three key parameters. Factors influencing the accuracy of scoring were investigated, including the maximal atomic distance and non-native ligand geometries used for scoring, as well as the use of protein models instead of crystallographic structures for training and testing the scoring function. For the test set of 19 targets, RankScore improved the ligand enrichment (logAUC) and early enrichment (EF(1)) scores computed by DOCK 3.6 for 13 and 14 targets, respectively. In addition, RankScore performed better at rescoring than each of seven other scoring functions tested. Accepting both the crystal structure and decoy geometries with all-atom root-mean-square errors of up to 2 Å from the crystal structure as correct binding poses, PoseScore gave the best score to a correct binding pose among 100 decoys for 88% of all cases in a benchmark set containing 100 protein-ligand complexes. PoseScore accuracy is comparable to that of DrugScore(CSD) and ITScore/SE and superior to 12 other tested scoring functions. Therefore, RankScore can facilitate ligand discovery, by ranking complexes of the target with different small molecules; PoseScore can be used for protein-ligand complex structure prediction, by ranking different conformations of a given protein-ligand pair. The statistical potentials are available through the Integrative Modeling Platform (IMP) software package (http://salilab.org/imp) and the LigScore Web server (http://salilab.org/ligscore/).

Ligand discovery from a dopamine D3 receptor homology model and crystal structure

G protein-coupled receptors (GPCRs) are intensely studied as drug targets and for their role in signaling. With the determination of the first crystal structures, interest in structure-based ligand discovery increased. Unfortunately, for most GPCRs no experimental structures are available. The determination of the D(3) receptor structure and the challenge to the community to predict it enabled a fully prospective comparison of ligand discovery from a modeled structure versus that of the subsequently released crystal structure. Over 3.3 million molecules were docked against a homology model, and 26 of the highest ranking were tested for binding. Six had affinities ranging from 0.2 to 3.1 μM. Subsequently, the crystal structure was released and the docking screen repeated. Of the 25 compounds selected, five had affinities ranging from 0.3 to 3.0 μM. One of the new ligands from the homology model screen was optimized for affinity to 81 nM. The feasibility of docking screens against modeled GPCRs more generally is considered.

The presynaptic component of the serotonergic system is required for clozapine's efficacy

Clozapine, by virtue of its absence of extrapyramidal side effects and greater efficacy, revolutionized the treatment of schizophrenia, although the mechanisms underlying this exceptional activity remain controversial. Combining an unbiased cheminformatics and physical screening approach, we evaluated clozapine's activity at >2350 distinct molecular targets. Clozapine, and the closely related atypical antipsychotic drug olanzapine, interacted potently with a unique spectrum of molecular targets. This distinct pattern, which was not shared with the typical antipsychotic drug haloperidol, suggested that the serotonergic neuronal system was a key determinant of clozapine's actions. To test this hypothesis, we used pet1(-/-) mice, which are deficient in serotonergic presynaptic markers. We discovered that the antipsychotic-like properties of the atypical antipsychotic drugs clozapine and olanzapine were abolished in a pharmacological model that mimics NMDA-receptor hypofunction in pet1(-/-) mice, whereas haloperidol's efficacy was unaffected. These results show that clozapine's ability to normalize NMDA-receptor hypofunction, which is characteristic of schizophrenia, depends on an intact presynaptic serotonergic neuronal system.

The chemical basis of pharmacology

Molecular biology now dominates pharmacology so thoroughly that it is difficult to recall that only a generation ago the field was very different. To understand drug action today, we characterize the targets through which they act and new drug leads are discovered on the basis of target structure and function. Until the mid-1980s the information often flowed in reverse: investigators began with organic molecules and sought targets, relating receptors not by sequence or structure but by their ligands. Recently, investigators have returned to this chemical view of biology, bringing to it systematic and quantitative methods of relating targets by their ligands. This has allowed the discovery of new targets for established drugs, suggested the bases for their side effects, and predicted the molecular targets underlying phenotypic screens. The bases for these new methods, some of their successes and liabilities, and new opportunities for their use are described.

Complementarity between a docking and a high-throughput screen in discovering new cruzain inhibitors

Virtual and high-throughput screens (HTS) should have complementary strengths and weaknesses, but studies that prospectively and comprehensively compare them are rare. We undertook a parallel docking and HTS screen of 197861 compounds against cruzain, a thiol protease target for Chagas disease, looking for reversible, competitive inhibitors. On workup, 99% of the hits were eliminated as false positives, yielding 146 well-behaved, competitive ligands. These fell into five chemotypes: two were prioritized by scoring among the top 0.1% of the docking-ranked library, two were prioritized by behavior in the HTS and by clustering, and one chemotype was prioritized by both approaches. Determination of an inhibitor/cruzain crystal structure and comparison of the high-scoring docking hits to experiment illuminated the origins of docking false-negatives and false-positives. Prioritizing molecules that are both predicted by docking and are HTS-active yields well-behaved molecules, relatively unobscured by the false-positives to which both techniques are individually prone.