AlphaFold2 structures guide prospective ligand discovery.

AlphaFold2 (AF2) models have had wide impact but mixed success in retrospective ligand recognition. We prospectively docked large libraries against unrefined AF2 models of the σ2 and serotonin 2A (5-HT2A) receptors, testing hundreds of new molecules and comparing results with those obtained from docking against the experimental structures. Hit rates were high and similar for the experimental and AF2 structures, as were affinities. Success in docking against the AF2 models was achieved despite differences between orthosteric residue conformations in the AF2 models and the experimental structures. Determination of the cryo-electron microscopy structure for one of the more potent 5-HT2A ligands from the AF2 docking revealed residue accommodations that resembled the AF2 prediction. AF2 models may sample conformations that differ from experimental structures but remain low energy and relevant for ligand discovery, extending the domain of structure-based drug design.

Novel Fragment Inhibitors of PYCR1 from Docking-Guided X-ray Crystallography.

The proline biosynthetic enzyme Δ1-pyrroline-5-carboxylate (P5C) reductase 1 (PYCR1) is one of the most consistently upregulated enzymes across multiple cancer types and central to the metabolic rewiring of cancer cells. Herein, we describe a fragment-based, structure-first approach to the discovery of PYCR1 inhibitors. Thirty-seven fragment-like carboxylic acids in the molecular weight range of 143-289 Da were selected from docking and then screened using X-ray crystallography as the primary assay. Strong electron density was observed for eight compounds, corresponding to a crystallographic hit rate of 22%. The fragments are novel compared to existing proline analog inhibitors in that they block both the P5C substrate pocket and the NAD(P)H binding site. Four hits showed inhibition of PYCR1 in kinetic assays, and one has lower apparent IC50 than the current best proline analog inhibitor. These results show proof-of-concept for our inhibitor discovery approach and provide a basis for fragment-to-lead optimization.

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.

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.

AI-Powered Virtual Screening of Large Compound Libraries Leads to the Discovery of Novel Inhibitors of Sirtuin-1.

The discovery of new scaffolds and chemotypes via high-throughput screening is tedious and resource intensive. Yet, there are millions of small molecules commercially available, rendering comprehensive in vitro tests intractable. We show how smart algorithms reduce large screening collections to target-specific sets of just a few hundred small molecules, allowing for a much faster and more cost-effective hit discovery process. We showcase the application of this virtual screening strategy by preselecting 434 compounds for Sirtuin-1 inhibition from a library of 2.6 million compounds, corresponding to 0.02% of the original library. Multistage in vitro validation ultimately confirmed nine chemically novel inhibitors. When compared to a competitive benchmark study for Sirtuin-1, our method shows a 12-fold higher hit rate. The results demonstrate how AI-driven preselection from large screening libraries allows for a massive reduction in the number of small molecules to be tested in vitro while still retaining a large number of hits.

Magnet for the Needle in Haystack: "Crystal Structure First" Fragment Hits Unlock Active Chemical Matter Using Targeted Exploration of Vast Chemical Spaces.

Fragment-based drug discovery (FBDD) has successfully led to approved therapeutics for challenging and "undruggable" targets. In the context of FBDD, we introduce a novel, multidisciplinary method to identify active molecules from purchasable chemical space. Starting from four small-molecule fragment complexes of protein kinase A (PKA), a template-based docking screen using Enamine's multibillion REAL Space was performed. A total of 93 molecules out of 106 selected compounds were successfully synthesized. Forty compounds were active in at least one validation assay with the most active follow-up having a 13,500-fold gain in affinity. Crystal structures for six of the most promising binders were rapidly obtained, verifying the binding mode. The overall success rate for this novel fragment-to-hit approach was 40%, accomplished in only 9 weeks. The results challenge the established fragment prescreening paradigm since the standard industrial filters for fragment hit identification in a thermal shift assay would have missed the initial fragments.

A multi-pronged approach targeting SARS-CoV-2 proteins using ultra-large virtual screening.

The unparalleled global effort to combat the continuing severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic over the last year has resulted in promising prophylactic measures. However, a need still exists for cheap, effective therapeutics, and targeting multiple points in the viral life cycle could help tackle the current, as well as future, coronaviruses. Here, we leverage our recently developed, ultra-large-scale in silico screening platform, VirtualFlow, to search for inhibitors that target SARS-CoV-2. In this unprecedented structure-based virtual campaign, we screened roughly 1 billion molecules against each of 40 different target sites on 17 different potential viral and host targets. In addition to targeting the active sites of viral enzymes, we also targeted critical auxiliary sites such as functionally important protein-protein interactions.

A Multi-Pronged Approach Targeting SARS-CoV-2 Proteins Using Ultra-Large Virtual Screening.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), previously known as 2019 novel coronavirus (2019-nCoV), has spread rapidly across the globe, creating an unparalleled global health burden and spurring a deepening economic crisis. As of July 7th, 2020, almost seven months into the outbreak, there are no approved vaccines and few treatments available. Developing drugs that target multiple points in the viral life cycle could serve as a strategy to tackle the current as well as future coronavirus pandemics. Here we leverage the power of our recently developed in silico screening platform, VirtualFlow, to identify inhibitors that target SARS-CoV-2. VirtualFlow is able to efficiently harness the power of computing clusters and cloud-based computing platforms to carry out ultra-large scale virtual screens. In this unprecedented structure-based multi-target virtual screening campaign, we have used VirtualFlow to screen an average of approximately 1 billion molecules against each of 40 different target sites on 17 different potential viral and host targets in the cloud. In addition to targeting the active sites of viral enzymes, we also target critical auxiliary sites such as functionally important protein-protein interaction interfaces. This multi-target approach not only increases the likelihood of finding a potent inhibitor, but could also help identify a collection of anti-coronavirus drugs that would retain efficacy in the face of viral mutation. Drugs belonging to different regimen classes could be combined to develop possible combination therapies, and top hits that bind at highly conserved sites would be potential candidates for further development as coronavirus drugs. Here, we present the top 200 in silico hits for each target site. While in-house experimental validation of some of these compounds is currently underway, we want to make this array of potential inhibitor candidates available to researchers worldwide in consideration of the pressing need for fast-tracked drug development.