The novel coronavirus SARS-CoV-2, the causative agent of COVID-19 respiratory disease, has infected over 2.3 million people, killed over 160,000, and caused worldwide social and economic disruption 1,2 . There are currently no antiviral drugs with proven clinical efficacy, nor are there vaccines for its prevention, and these efforts are hampered by limited knowledge of the molecular details of SARS-CoV-2 infection. To address this, we cloned, tagged and expressed 26 of the 29 SARS-CoV-2 proteins in human cells and identified the human proteins physically associated with each using affinity-purification mass spectrometry (AP-MS), identifying 332 high-confidence SARS-CoV-2-human protein-protein interactions (PPIs). Among these, we identify 66 druggable human proteins or host factors targeted by 69 compounds (29 FDA-approved drugs, 12 drugs in clinical trials, and 28 preclinical compounds). Screening a subset of these in multiple viral assays identified two sets of pharmacological agents that displayed antiviral activity: inhibitors of mRNA translation and predicted regulators of the Sigma1 and Sigma2 receptors. Further studies of these host factor targeting agents, including their combination with drugs that directly target viral enzymes, could lead to a therapeutic regimen to treat COVID-19.

A critical barrier to entry into structure-based virtual screening is the lack of a suitable, easy to access database of purchasable compounds. We have therefore prepared a library of 727 842 molecules, each with 3D structure, using catalogs of compounds from vendors (the size of this library continues to grow). The molecules have been assigned biologically relevant protonation states and are annotated with properties such as molecular weight, calculated LogP, and number of rotatable bonds. Each molecule in the library contains vendor and purchasing information and is ready for docking using a number of popular docking programs. Within certain limits, the molecules are prepared in multiple protonation states and multiple tautomeric forms. In one format, multiple conformations are available for the molecules. This database is available for free download (http://zinc.docking.org) in several common file formats including SMILES, mol2, 3D SDF, and DOCK flexibase format. A Web-based query tool incorporating a molecular drawing interface enables the database to be searched and browsed and subsets to be created. Users can process their own molecules by uploading them to a server. Our hope is that this database will bring virtual screening libraries to a wide community of structural biologists and medicinal chemists.

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 () 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 .

Ligand enrichment among top-ranking hits is a key metric of molecular docking. To avoid bias, decoys should resemble ligands physically, so that enrichment is not simply a separation of gross features, yet be chemically distinct from them, so that they are unlikely to be binders. We have assembled a directory of useful decoys (DUD), with 2950 ligands for 40 different targets. Every ligand has 36 decoy molecules that are physically similar but topologically distinct, leading to a database of 98,266 compounds. For most targets, enrichment was at least half a log better with uncorrected databases such as the MDDR than with DUD, evidence of bias in the former. These calculations also allowed forty-by-forty cross docking, where the enrichments of each ligand set could be compared for all 40 targets, enabling a specificity metric for the docking screens. DUD is freely available online as a benchmarking set for docking at http://blaster.docking.org/dud/.

The identification of protein function based on biological information is an area of intense research. Here we consider a complementary technique that quantitatively groups and relates proteins based on the chemical similarity of their ligands. We began with 65,000 ligands annotated into sets for hundreds of drug targets. The similarity score between each set was calculated using ligand topology. A statistical model was developed to rank the significance of the resulting similarity scores, which are expressed as a minimum spanning tree to map the sets together. Although these maps are connected solely by chemical similarity, biologically sensible clusters nevertheless emerged. Links among unexpected targets also emerged, among them that methadone, emetine and loperamide (Imodium) may antagonize muscarinic M3, alpha2 adrenergic and neurokinin NK2 receptors, respectively. These predictions were subsequently confirmed experimentally. Relating receptors by ligand chemistry organizes biology to reveal unexpected relationships that may be assayed using the ligands themselves.

Highlights d Phosphoproteomics analysis of SARS-CoV-2-infected cells uncovers signaling rewiring d Infection promotes host p38 MAPK cascade activity and shutdown of mitotic kinases d Infection stimulates CK2-containing filopodial protrusions with budding virus d Kinase activity analysis identifies potent antiviral drugs and compounds

Whereas drugs are intended to be selective, at least some bind to several physiologic targets, explaining both side effects and efficacy. As many drug-target combinations exist, it would be useful to explore possible interactions computationally. Here, we compared 3,665 FDA-approved and investigational drugs against hundreds of targets, defining each target by its ligands. Chemical similarities between drugs and ligand sets predicted thousands of unanticipated associations. Thirty were tested experimentally, including the antagonism of the β1 receptor by the transporter inhibitor Prozac, the inhibition of the 5-HT transporter by the ion channel drug Vadilex, and antagonism of the histamine H4 receptor by the enzyme inhibitor Rescriptor. Overall, 23 new drug-target associations were confirmed, five of which were potent (< 100 nM). The physiological relevance of one such, the drug DMT on serotonergic receptors, was confirmed in a knock-out mouse. The chemical similarity approach is systematic and comprehensive, and may suggest side-effects and new indications for many drugs.

High-throughput and virtual screening are widely used to discover novel leads for drug design. On examination, many screening hits appear non-drug-like: they act noncompetitively, show little relationship between structure and activity, and have poor selectivity. Attempts to develop these peculiar molecules into viable leads are often futile, and much time can be wasted on the characterization of these "phony" hits. Despite their common occurrence, the mechanism of action of these promiscuous molecules remains unknown. To investigate this problem, 45 diverse screening hits were studied. Fifteen of these were previously reported as inhibitors of various receptors, including beta-lactamase, malarial protease, dihydrofolate reductase, HIV Tar RNA, thymidylate synthase, kinesin, insulin receptor, tyrosine kinases, farnesyltransferase, gyrase, prions, triosephosphate isomerase, nitric oxide synthase, phosphoinositide 3-kinase, and integrase; 30 were from an in-house screening library of a major pharmaceutical company. In addition to their original targets, 35 of these 45 compounds were shown to inhibit several unrelated model enzymes. These 35 screening hits included compounds, such as fullerenes, dyes, and quercetin, that have repeatedly shown activity against diverse targets. When tested against the model enzymes, the compounds showed time-dependent but reversible inhibition that was dramatically attenuated by albumin, guanidinium, or urea. Surprisingly, increasing the concentration of the model enzymes 10-fold largely eliminated inhibition, despite a 1000-fold excess of inhibitor; a well-behaved competitive inhibitor did not show this behavior. One model to explain these observations was that the active form of the promiscuous inhibitors was an aggregate of many individual molecules. To test this hypothesis, light scattering and electron microscopy experiments were performed. The nonspecific inhibitors were observed to form particles of 30-400 nm diameter by both techniques. In control experiments, a well-behaved competitive inhibitor and an inactive dye-like molecule were not observed to form aggregates. Consistent with the hypothesis that the aggregates are the inhibitory species, the particle size and IC(50) values of the promiscuous inhibitors varied monotonically with ionic strength; a competitive inhibitor was unaffected by changes in ionic strength. Unexpectedly, aggregate formation appears to explain the activity of many nonspecific inhibitors and may account for the activity of many promiscuous screening hits. Molecules acting via this mechanism may be widespread in drug discovery screening databases. Recognition of these compounds may improve screening results in many areas of pharmaceutical interest.