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.

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 (α 2A AR), seeking new α 2A AR agonists chemotypes that lack the sedation conferred by known α 2A AR drugs, such as dexmedetomidine. We identified 17 ligands with potencies as low as 12 nanomolar, many with partial agonism and preferential G i and G o signaling. Experimental structures of α 2A AR 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 (EC 50 ) of 52 nanomolar] and two analogs, ‘7075 and PS75 (EC 50 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.

Antiviral therapeutics to treat SARS-CoV-2 are much desired for the on-going pandemic. A well-precedented viral enzyme is the main protease (MPro), which is now targeted by an approved drug and by several investigational drugs. With the inevitable liabilities of these new drugs, and facing viral resistance, there remains a call for new chemical scaffolds against MPro. We virtually docked 1.2 billion non-covalent and a new library of 6.5 million electrophilic molecules 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 μM and 20 μM, respectively. Several series were optimized, resulting in inhibitors active in the low micromolar range. Subsequent crystallography confirmed the docking predicted binding modes and may template further optimization. Together, these compounds reveal new chemotypes to aid in further discovery of MPro inhibitors for SARS-CoV-2 and other future coronaviruses.

Protein conformations are shaped by cellular environments, but how environmental changes alter the conformational landscapes of specific proteins in vivo remains largely uncharacterized, in part due to the challenge of probing protein structures in living cells. Here, we use deep mutational scanning to investigate how a toxic conformation of α-synuclein, a dynamic protein linked to Parkinson's disease, responds to perturbations of cellular proteostasis. In the context of a course for graduate students in the UCSF Integrative Program in Quantitative Biology, we screened a comprehensive library of α-synuclein missense mutants in yeast cells treated with a variety of small molecules that perturb cellular processes linked to α-synuclein biology and pathobiology. We found that the conformation of α-synuclein previously shown to drive yeast toxicity-an extended, membrane-bound helix-is largely unaffected by these chemical perturbations, underscoring the importance of this conformational state as a driver of cellular toxicity. On the other hand, the chemical perturbations have a significant effect on the ability of mutations to suppress α-synuclein toxicity. Moreover, we find that sequence determinants of αsynuclein toxicity are well described by a simple structural model of the membrane-bound helix.This model predicts that α-synuclein penetrates the membrane to constant depth across its length but that membrane affinity decreases toward the C terminus, which is consistent with orthogonal biophysical measurements. Finally, we discuss how parallelized chemical genetics experiments can provide a robust framework for inquiry-based graduate coursework.