Open science discovery of potent noncovalent SARS-CoV-2 main protease inhibitors.

We report the results of the COVID Moonshot, a fully open-science, crowdsourced, and structure-enabled drug discovery campaign targeting the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease. We discovered a noncovalent, nonpeptidic inhibitor scaffold with lead-like properties that is differentiated from current main protease inhibitors. Our approach leveraged crowdsourcing, machine learning, exascale molecular simulations, and high-throughput structural biology and chemistry. We generated a detailed map of the structural plasticity of the SARS-CoV-2 main protease, extensive structure-activity relationships for multiple chemotypes, and a wealth of biochemical activity data. All compound designs (>18,000 designs), crystallographic data (>490 ligand-bound x-ray structures), assay data (>10,000 measurements), and synthesized molecules (>2400 compounds) for this campaign were shared rapidly and openly, creating a rich, open, and intellectual property-free knowledge base for future anticoronavirus drug discovery.

Allosteric regulation and crystallographic fragment screening of SARS-CoV-2 NSP15 endoribonuclease.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of coronavirus disease 2019 (COVID-19). The NSP15 endoribonuclease enzyme, known as NendoU, is highly conserved and plays a critical role in the ability of the virus to evade the immune system. NendoU is a promising target for the development of new antiviral drugs. However, the complexity of the enzyme's structure and kinetics, along with the broad range of recognition sequences and lack of structural complexes, hampers the development of inhibitors. Here, we performed enzymatic characterization of NendoU in its monomeric and hexameric form, showing that hexamers are allosteric enzymes with a positive cooperative index, and with no influence of manganese on enzymatic activity. Through combining cryo-electron microscopy at different pHs, X-ray crystallography and biochemical and structural analysis, we showed that NendoU can shift between open and closed forms, which probably correspond to active and inactive states, respectively. We also explored the possibility of NendoU assembling into larger supramolecular structures and proposed a mechanism for allosteric regulation. In addition, we conducted a large fragment screening campaign against NendoU and identified several new allosteric sites that could be targeted for the development of new inhibitors. Overall, our findings provide insights into the complex structure and function of NendoU and offer new opportunities for the development of inhibitors.

Room-temperature crystallography reveals altered binding of small-molecule fragments to PTP1B.

Much of our current understanding of how small-molecule ligands interact with proteins stems from X-ray crystal structures determined at cryogenic (cryo) temperature. For proteins alone, room-temperature (RT) crystallography can reveal previously hidden, biologically relevant alternate conformations. However, less is understood about how RT crystallography may impact the conformational landscapes of protein-ligand complexes. Previously, we showed that small-molecule fragments cluster in putative allosteric sites using a cryo crystallographic screen of the therapeutic target PTP1B (Keedy et al., 2018). Here, we have performed two RT crystallographic screens of PTP1B using many of the same fragments, representing the largest RT crystallographic screens of a diverse library of ligands to date, and enabling a direct interrogation of the effect of data collection temperature on protein-ligand interactions. We show that at RT, fewer ligands bind, and often more weakly - but with a variety of temperature-dependent differences, including unique binding poses, changes in solvation, new binding sites, and distinct protein allosteric conformational responses. Overall, this work suggests that the vast body of existing cryo-temperature protein-ligand structures may provide an incomplete picture, and highlights the potential of RT crystallography to help complete this picture by revealing distinct conformational modes of protein-ligand systems. Our results may inspire future use of RT crystallography to interrogate the roles of protein-ligand conformational ensembles in biological function.

Discovery of novel druggable pockets on polyomavirus VP1 through crystallographic fragment-based screening to develop capsid assembly inhibitors.

Polyomaviruses are a family of ubiquitous double-stranded DNA viruses many of which are human pathogens. These include BK polyomavirus which causes severe urinary tract infection in immunocompromised patients and Merkel cell polyomavirus associated with aggressive cancers. The small genome of polyomaviruses lacks conventional drug targets, and no specific drugs are available at present. Here we focus on the main structural protein VP1 of BK polyomavirus which is responsible for icosahedral capsid formation. To provide a foundation towards rational drug design, we crystallized truncated VP1 pentamers and subjected them to a high-throughput screening for binding drug-like fragments through a direct X-ray analysis. To enable a highly performant screening, rigorous optimization of the crystallographic pipeline and processing with the latest generation PanDDA2 software were necessary. As a result, a total of 144 binding hits were established. Importantly, the hits are well clustered in six surface pockets. Three pockets are located on the outside of the pentamer and map on the regions where the 'invading' C-terminal arm of another pentamer is attached upon capsid assembly. Another set of three pockets is situated within the wide pore along the five-fold axis of the VP1 pentamer. These pockets are situated at the interaction interface with the minor capsid protein VP2 which is indispensable for normal functioning of the virus. Here we systematically analyse the three outside pockets which are highly conserved across various polyomaviruses, while point mutations in these pockets are detrimental for viral replication. We show that one of the pockets can accommodate antipsychotic drug trifluoperazine. For each pocket, we derive pharmacophore features which enable the design of small molecules preventing the interaction between VP1 pentamers and therefore inhibiting capsid assembly. Our data lay a foundation towards a rational development of first-in-class drugs targeting polyomavirus capsid.

Structure, mechanism and crystallographic fragment screening of the SARS-CoV-2 NSP13 helicase.

There is currently a lack of effective drugs to treat people infected with SARS-CoV-2, the cause of the global COVID-19 pandemic. The SARS-CoV-2 Non-structural protein 13 (NSP13) has been identified as a target for anti-virals due to its high sequence conservation and essential role in viral replication. Structural analysis reveals two "druggable" pockets on NSP13 that are among the most conserved sites in the entire SARS-CoV-2 proteome. Here we present crystal structures of SARS-CoV-2 NSP13 solved in the APO form and in the presence of both phosphate and a non-hydrolysable ATP analog. Comparisons of these structures reveal details of conformational changes that provide insights into the helicase mechanism and possible modes of inhibition. To identify starting points for drug development we have performed a crystallographic fragment screen against NSP13. The screen reveals 65 fragment hits across 52 datasets opening the way to structure guided development of novel antiviral agents.

Achieving Efficient Fragment Screening at XChem Facility at Diamond Light Source.

In fragment-based drug discovery, hundreds or often thousands of compounds smaller than ~300 Da are tested against the protein of interest to identify chemical entities that can be developed into potent drug candidates. Since the compounds are small, interactions are weak, and the screening method must therefore be highly sensitive; moreover, structural information tends to be crucial for elaborating these hits into lead-like compounds. Therefore, protein crystallography has always been a gold-standard technique, yet historically too challenging to find widespread use as a primary screen. Initial XChem experiments were demonstrated in 2014 and then trialed with academic and industrial collaborators to validate the process. Since then, a large research effort and significant beamtime have streamlined sample preparation, developed a fragment library with rapid follow-up possibilities, automated and improved the capability of I04-1 beamline for unattended data collection, and implemented new tools for data management, analysis and hit identification. XChem is now a facility for large-scale crystallographic fragment screening, supporting the entire crystals-to-deposition process, and accessible to academic and industrial users worldwide. The peer-reviewed academic user program has been actively developed since 2016, to accommodate projects from as broad a scientific scope as possible, including well-validated as well as exploratory projects. Academic access is allocated through biannual calls for peer-reviewed proposals, and proprietary work is arranged by Diamond's Industrial Liaison group. This workflow has already been routinely applied to over a hundred targets from diverse therapeutic areas, and effectively identifies weak binders (1%-30% hit rate), which both serve as high-quality starting points for compound design and provide extensive structural information on binding sites. The resilience of the process was demonstrated by continued screening of SARS-CoV-2 targets during the COVID-19 pandemic, including a 3-week turn-around for the main protease.

Bispecific repurposed medicines targeting the viral and immunological arms of COVID-19.

Effective agents to treat coronavirus infection are urgently required, not only to treat COVID-19, but to prepare for future outbreaks. Repurposed anti-virals such as remdesivir and human anti-inflammatories such as barcitinib have received emergency approval but their overall benefits remain unclear. Vaccines are the most promising prospect for COVID-19, but will need to be redeveloped for any future coronavirus outbreak. Protecting against future outbreaks requires the identification of targets that are conserved between coronavirus strains and amenable to drug discovery. Two such targets are the main protease (Mpro) and the papain-like protease (PLpro) which are essential for the coronavirus replication cycle. We describe the discovery of two non-antiviral therapeutic agents, the caspase-1 inhibitor SDZ 224015 and Tarloxotinib that target Mpro and PLpro, respectively. These were identified through extensive experimental screens of the drug repurposing ReFRAME library of 12,000 therapeutic agents. The caspase-1 inhibitor SDZ 224015, was found to be a potent irreversible inhibitor of Mpro (IC50 30 nM) while Tarloxotinib, a clinical stage epidermal growth factor receptor inhibitor, is a sub micromolar inhibitor of PLpro (IC50 300 nM, Ki 200 nM) and is the first reported PLpro inhibitor with drug-like properties. SDZ 224015 and Tarloxotinib have both undergone safety evaluation in humans and hence are candidates for COVID-19 clinical evaluation.

An automatic pipeline for the design of irreversible derivatives identifies a potent SARS-CoV-2 Mpro inhibitor.

Designing covalent inhibitors is increasingly important, although it remains challenging. Here, we present covalentizer, a computational pipeline for identifying irreversible inhibitors based on structures of targets with non-covalent binders. Through covalent docking of tailored focused libraries, we identify candidates that can bind covalently to a nearby cysteine while preserving the interactions of the original molecule. We found âˆ¼11,000 cysteines proximal to a ligand across 8,386 complexes in the PDB. Of these, the protocol identified 1,553 structures with covalent predictions. In a prospective evaluation, five out of nine predicted covalent kinase inhibitors showed half-maximal inhibitory concentration (IC50) values between 155 nM and 4.5 µM. Application against an existing SARS-CoV Mpro reversible inhibitor led to an acrylamide inhibitor series with low micromolar IC50 values against SARS-CoV-2 Mpro. The docking was validated by 12 co-crystal structures. Together these examples hint at the vast number of covalent inhibitors accessible through our protocol.

Exploring protein hotspots by optimized fragment pharmacophores.

Fragment-based drug design has introduced a bottom-up process for drug development, with improved sampling of chemical space and increased effectiveness in early drug discovery. Here, we combine the use of pharmacophores, the most general concept of representing drug-target interactions with the theory of protein hotspots, to develop a design protocol for fragment libraries. The SpotXplorer approach compiles small fragment libraries that maximize the coverage of experimentally confirmed binding pharmacophores at the most preferred hotspots. The efficiency of this approach is demonstrated with a pilot library of 96 fragment-sized compounds (SpotXplorer0) that is validated on popular target classes and emerging drug targets. Biochemical screening against a set of GPCRs and proteases retrieves compounds containing an average of 70% of known pharmacophores for these targets. More importantly, SpotXplorer0 screening identifies confirmed hits against recently established challenging targets such as the histone methyltransferase SETD2, the main protease (3CLPro) and the NSP3 macrodomain of SARS-CoV-2.

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.

Fragment Binding to the Nsp3 Macrodomain of SARS-CoV-2 Identified Through Crystallographic Screening and Computational Docking.

The SARS-CoV-2 macrodomain (Mac1) within the non-structural protein 3 (Nsp3) counteracts host-mediated antiviral ADP-ribosylation signalling. This enzyme is a promising antiviral target because catalytic mutations render viruses non-pathogenic. 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 diverse fragment libraries resulted in 214 unique macrodomain-binding fragments, out of 2,683 screened. An additional 60 molecules were selected from docking over 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 crystallographic and docking fragment hits were validated for solution binding using three biophysical techniques (DSF, HTRF, ITC). Overall, the 234 fragment structures presented explore a wide range of chemotypes and provide starting points for development of potent SARS-CoV-2 macrodomain inhibitors.

Crystallographic and electrophilic fragment screening of the SARS-CoV-2 main protease.

COVID-19, caused by SARS-CoV-2, lacks effective therapeutics. Additionally, no antiviral drugs or vaccines were developed against the closely related coronavirus, SARS-CoV-1 or MERS-CoV, despite previous zoonotic outbreaks. To identify starting points for such therapeutics, we performed a large-scale screen of electrophile and non-covalent fragments through a combined mass spectrometry and X-ray approach against the SARS-CoV-2 main protease, one of two cysteine viral proteases essential for viral replication. Our crystallographic screen identified 71 hits that span the entire active site, as well as 3 hits at the dimer interface. These structures reveal routes to rapidly develop more potent inhibitors through merging of covalent and non-covalent fragment hits; one series of low-reactivity, tractable covalent fragments were progressed to discover improved binders. These combined hits offer unprecedented structural and reactivity information for on-going structure-based drug design against SARS-CoV-2 main protease.

In crystallo-screening for discovery of human norovirus 3C-like protease inhibitors.

Outbreaks of human epidemic nonbacterial gastroenteritis are mainly caused by noroviruses. Viral replication requires a 3C-like cysteine protease (3CLpro) which processes the 200 kDa viral polyprotein into six functional proteins. The 3CLpro has attracted much interest due to its potential as a target for antiviral drugs. A system for growing high-quality crystals of native Southampton norovirus 3CLpro (SV3CP) has been established, allowing the ligand-free crystal structure to be determined to 1.3 Å in a tetrameric state. This also allowed crystal-based fragment screening to be performed with various compound libraries, ultimately to guide drug discovery for SV3CP. A total of 19 fragments were found to bind to the protease out of the 844 which were screened. Two of the hits were located at the active site of SV3CP and showed good inhibitory activity in kinetic assays. Another 5 were found at the enzyme's putative RNA-binding site and a further 11 were located in the symmetric central cavity of the tetramer.

Crystal structure of the catalytic D2 domain of the AAA+ ATPase p97 reveals a putative helical split-washer-type mechanism for substrate unfolding.

Several pathologies have been associated with the AAA+ ATPase p97, an enzyme essential to protein homeostasis. Heterozygous polymorphisms in p97 have been shown to cause neurological disease, while elevated proteotoxic stress in tumours has made p97 an attractive cancer chemotherapy target. The cellular processes reliant on p97 are well described. High-resolution structural models of its catalytic D2 domain, however, have proved elusive, as has the mechanism by which p97 converts the energy from ATP hydrolysis into mechanical force to unfold protein substrates. Here, we describe the high-resolution structure of the p97 D2 ATPase domain. This crystal system constitutes a valuable tool for p97 inhibitor development and identifies a potentially druggable pocket in the D2 domain. In addition, its P61 symmetry suggests a mechanism for substrate unfolding by p97. DATABASE: The atomic coordinates and structure factors have been deposited in the PDB database under the accession numbers 6G2V, 6G2W, 6G2X, 6G2Y, 6G2Z and 6G30.

Rapid Covalent-Probe Discovery by Electrophile-Fragment Screening.

Covalent probes can display unmatched potency, selectivity, and duration of action; however, their discovery is challenging. In principle, fragments that can irreversibly bind their target can overcome the low affinity that limits reversible fragment screening, but such electrophilic fragments were considered nonselective and were rarely screened. We hypothesized that mild electrophiles might overcome the selectivity challenge and constructed a library of 993 mildly electrophilic fragments. We characterized this library by a new high-throughput thiol-reactivity assay and screened them against 10 cysteine-containing proteins. Highly reactive and promiscuous fragments were rare and could be easily eliminated. In contrast, we found hits for most targets. Combining our approach with high-throughput crystallography allowed rapid progression to potent and selective probes for two enzymes, the deubiquitinase OTUB2 and the pyrophosphatase NUDT7. No inhibitors were previously known for either. This study highlights the potential of electrophile-fragment screening as a practical and efficient tool for covalent-ligand discovery.

Achieving a Good Crystal System for Crystallographic X-Ray Fragment Screening.

The XChem facility at Diamond Light Source offers fragment screening by X-ray crystallography as a general access user program. The main advantage of X-ray crystallography as a primary fragment screen is that it yields directly the location and pose of the fragment hits, whether within pockets of interest or merely on surface sites: this is the key information for structure-based design and for enabling synthesis of follow-up molecules. Extensive streamlining of the screening experiment at XChem has engendered a very active user program that is generating large amounts of data: in 2017, 36 academic and industry groups generated 35,000 datasets of uniquely soaked crystals. It has also generated a large number of learnings concerning the main remaining bottleneck, namely, obtaining a suitable crystal system that will support a successful fragment screen. Here we discuss the practicalities of generating screen-ready crystals that have useful electron density maps, and how to ensure they will be successfully reproduced and usable at a facility outside the home lab.

Gentle, fast and effective crystal soaking by acoustic dispensing.

The steady expansion in the capacity of modern beamlines for high-throughput data collection, enabled by increasing X-ray brightness, capacity of robotics and detector speeds, has pushed the bottleneck upstream towards sample preparation. Even in ligand-binding studies using crystal soaking, the experiment best able to exploit beamline capacity, a primary limitation is the need for gentle and nontrivial soaking regimens such as stepwise concentration increases, even for robust and well characterized crystals. Here, the use of acoustic droplet ejection for the soaking of protein crystals with small molecules is described, and it is shown that it is both gentle on crystals and allows very high throughput, with 1000 unique soaks easily performed in under 10 min. In addition to having very low compound consumption (tens of nanolitres per sample), the positional precision of acoustic droplet ejection enables the targeted placement of the compound/solvent away from crystals and towards drop edges, allowing gradual diffusion of solvent across the drop. This ensures both an improvement in the reproducibility of X-ray diffraction and increased solvent tolerance of the crystals, thus enabling higher effective compound-soaking concentrations. The technique is detailed here with examples from the protein target JMJD2D, a histone lysine demethylase with roles in cancer and the focus of active structure-based drug-design efforts.

The XChemExplorer graphical workflow tool for routine or large-scale protein-ligand structure determination.

XChemExplorer (XCE) is a data-management and workflow tool to support large-scale simultaneous analysis of protein-ligand complexes during structure-based ligand discovery (SBLD). The user interfaces of established crystallographic software packages such as CCP4 [Winn et al. (2011), Acta Cryst. D67, 235-242] or PHENIX [Adams et al. (2010), Acta Cryst. D66, 213-221] have entrenched the paradigm that a `project' is concerned with solving one structure. This does not hold for SBLD, where many almost identical structures need to be solved and analysed quickly in one batch of work. Functionality to track progress and annotate structures is essential. XCE provides an intuitive graphical user interface which guides the user from data processing, initial map calculation, ligand identification and refinement up until data dissemination. It provides multiple entry points depending on the need of each project, enables batch processing of multiple data sets and records metadata, progress and annotations in an SQLite database. XCE is freely available and works on any Linux and Mac OS X system, and the only dependency is to have the latest version of CCP4 installed. The design and usage of this tool are described here, and its usefulness is demonstrated in the context of fragment-screening campaigns at the Diamond Light Source. It is routinely used to analyse projects comprising 1000 data sets or more, and therefore scales well to even very large ligand-design projects.

Application of in situ diffraction in high-throughput structure determination platforms.

Macromolecular crystallography (MX) is the most powerful technique available to structural biologists to visualize in atomic detail the macromolecular machinery of the cell. Since the emergence of structural genomics initiatives, significant advances have been made in all key steps of the structure determination process. In particular, third-generation synchrotron sources and the application of highly automated approaches to data acquisition and analysis at these facilities have been the major factors in the rate of increase of macromolecular structures determined annually. A plethora of tools are now available to users of synchrotron beamlines to enable rapid and efficient evaluation of samples, collection of the best data, and in favorable cases structure solution in near real time. Here, we provide a short overview of the emerging use of collecting X-ray diffraction data directly from the crystallization experiment. These in situ experiments are now routinely available to users at a number of synchrotron MX beamlines. A practical guide to the use of the method on the MX suite of beamlines at Diamond Light Source is given.