AlphaFold2 structures guide prospective ligand discovery.

AlphaFold2 (AF2) models have had wide impact, but they have had mixed success in retrospective ligand recognition. We prospectively docked large libraries against unrefined AF2 models of the σ2 and 5-HT2A receptors, testing hundreds of new molecules and comparing results to docking against the experimental structures. Hit rates were high and similar for the experimental and the AF2 structures, as were affinities. The success of docking against the AF2 models was achieved despite differences in orthosteric residue conformations versus the experimental structures. Determination of the cryo-electron microscopy structure for one of the more potent 5HT2A 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 ligand discovery.

Discovery of Protease-Activated Receptor 4 (PAR4)-Tethered Ligand Antagonists Using Ultralarge Virtual Screening.

Here, we demonstrate a structure-based small molecule virtual screening and lead optimization pipeline using a homology model of a difficult-to-drug G-protein-coupled receptor (GPCR) target. Protease-activated receptor 4 (PAR4) is activated by thrombin cleavage, revealing a tethered ligand that activates the receptor, making PAR4 a challenging target. A virtual screen of a make-on-demand chemical library yielded a one-hit compound. From the single-hit compound, we developed a novel series of PAR4 antagonists. Subsequent lead optimization via simultaneous virtual library searches and structure-based rational design efforts led to potent antagonists of thrombin-induced activation. Interestingly, this series of antagonists was active against PAR4 activation by the native protease thrombin cleavage but not the synthetic PAR4 agonist peptide AYPGKF.

Profiling the proximal proteome of the activated µ-opioid receptor.

The µ-opioid receptor (µOR) represents an important target of therapeutic and abused drugs. So far, most understanding of µOR activity has focused on a subset of known signal transducers and regulatory molecules. Yet µOR signaling is coordinated by additional proteins in the interaction network of the activated receptor, which have largely remained invisible given the lack of technologies to interrogate these networks systematically. Here we describe a proteomics and computational approach to map the proximal proteome of the activated µOR and to extract subcellular location, trafficking and functional partners of G-protein-coupled receptor (GPCR) activity. We demonstrate that distinct opioid agonists exert differences in the µOR proximal proteome mediated by endocytosis and endosomal sorting. Moreover, we identify two new µOR network components, EYA4 and KCTD12, which are recruited on the basis of receptor-triggered G-protein activation and might form a previously unrecognized buffering system for G-protein activity broadly modulating cellular GPCR signaling.

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.

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.

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.

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.

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.

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.

Preparing for the next pandemic.

New lead drugs to treat COVID-19 are beginning to emerge.

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.

Colloidal aggregation confounds cell-based Covid-19 antiviral screens.

Colloidal aggregation is one of the largest contributors to false-positives in early drug discovery and chemical biology. Much work has focused on its impact on pure-protein screens; here we consider aggregations role in cell-based infectivity assays in Covid-19 drug repurposing. We began by investigating the potential aggregation of 41 drug candidates reported as SARs-CoV-2 entry inhibitors. Of these, 17 formed colloidal-particles by dynamic light scattering and exhibited detergent-dependent enzyme inhibition. To evaluate antiviral efficacy of the drugs in cells we used spike pseudotyped lentivirus and pre-saturation of the colloids with BSA. The antiviral potency of the aggregators was diminished by at least 10-fold and often entirely eliminated in the presence of BSA, suggesting antiviral activity can be attributed to the non-specific nature of the colloids. In confocal microscopy, the aggregates induced fluorescent puncta of labeled spike protein, consistent with sequestration of the protein on the colloidal particles. Addition of either non-ionic detergent or of BSA disrupted these puncta. These observations suggest that colloidal aggregation is common among cell-based anti-viral drug repurposing, and perhaps cell-based assays more broadly, and offers rapid counter-screens to detect and eliminate these artifacts, allowing the community invest resources in compounds with true potential as a Covid-19 therapeutic.

Cryo-EM structure of Alzheimer's disease tau filaments with PET ligand MK-6240.

Positron Emission Tomography (PET) ligands have advanced Alzheimer's disease (AD) diagnosis and treatment. Using autoradiography and cryo-EM, we identified AD brain tissue with elevated tau burden, purified filaments, and determined the structure of second-generation high avidity PET ligand MK-6240 at 2.31 Å resolution, which bound at a 1:1 ratio within the cleft of tau paired-helical filament (PHF), engaging with glutamine 351, lysine K353, and isoleucine 360. This information elucidates the basis of MK-6240 PET in quantifying PHF deposits in AD and may facilitate the structure-based design of superior ligands against tau amyloids.

Proteomic and genetic analyses of influenza A viruses identify pan-viral host targets.

Influenza A Virus (IAV) is a recurring respiratory virus with limited availability of antiviral therapies. Understanding host proteins essential for IAV infection can identify targets for alternative host-directed therapies (HDTs). Using affinity purification-mass spectrometry and global phosphoproteomic and protein abundance analyses using three IAV strains (pH1N1, H3N2, H5N1) in three human cell types (A549, NHBE, THP-1), we map 332 IAV-human protein-protein interactions and identify 13 IAV-modulated kinases. Whole exome sequencing of patients who experienced severe influenza reveals several genes, including scaffold protein AHNAK, with predicted loss-of-function variants that are also identified in our proteomic analyses. Of our identified host factors, 54 significantly alter IAV infection upon siRNA knockdown, and two factors, AHNAK and coatomer subunit COPB1, are also essential for productive infection by SARS-CoV-2. Finally, 16 compounds targeting our identified host factors suppress IAV replication, with two targeting CDK2 and FLT3 showing pan-antiviral activity across influenza and coronavirus families. This study provides a comprehensive network model of IAV infection in human cells, identifying functional host targets for pan-viral HDT.

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.

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

Synthetic Ionizable Colloidal Drug Aggregates Enable Endosomal Disruption.

Colloidal drug aggregates enable the design of drug-rich nanoparticles; however, the efficacy of stabilized colloidal drug aggregates is limited by entrapment in the endo-lysosomal pathway. Although ionizable drugs are used to elicit lysosomal escape, this approach is hindered by toxicity associated with phospholipidosis. It is hypothesized that tuning the pKa of the drug would enable endosomal disruption while avoiding phospholipidosis and minimizing toxicity. To test this idea, 12 analogs of the nonionizable colloidal drug fulvestrant are synthesized with ionizable groups to enable pH-dependent endosomal disruption while maintaining bioactivity. Lipid-stabilized fulvestrant analog colloids are endocytosed by cancer cells, and the pKa of these ionizable colloids influenced the mechanism of endosomal and lysosomal disruption. Four fulvestrant analogs-those with pKa values between 5.1 and 5.7-disrupted endo-lysosomes without measurable phospholipidosis. Thus, by manipulating the pKa of colloid-forming drugs, a tunable and generalizable strategy for endosomal disruption is established.

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.

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.