Coronaviruses can evolve and spread rapidly to cause severe disease morbidity and mortality, as exemplified by SARS-CoV-2 variants of the COVID-19 pandemic. Although currently available vaccines remain mostly effective against SARS-CoV-2 variants, additional treatment strategies are needed. Inhibitors that target essential viral enzymes, such as proteases and polymerases, represent key classes of antivirals. However, clinical use of antiviral therapies inevitably leads to emergence of drug resistance. In this study we implemented a strategy to pre-emptively address drug resistance to protease inhibitors targeting the main protease (Mpro) of SARS-CoV-2, an essential enzyme that promotes viral maturation. We solved nine high-resolution cocrystal structures of SARS-CoV-2 Mpro bound to substrate peptides and six structures with cleavage products. These structures enabled us to define the substrate envelope of Mpro, map the critical recognition elements, and identify evolutionarily vulnerable sites that may be susceptible to resistance mutations that would compromise binding of the newly developed Mpro inhibitors. Our results suggest strategies for developing robust inhibitors against SARS-CoV-2 that will retain longer-lasting efficacy against this evolving viral pathogen.

Herein we provide a living summary of the data generated during the COVID Moonshot project focused on the development of SARS-CoV-2 main protease (Mpro) inhibitors. Our approach uniquely combines crowdsourced medicinal chemistry insights with high throughput crystallography, exascale computational chemistry infrastructure for simulations, and machine learning in triaging designs and predicting synthetic routes. This manuscript describes our methodologies leading to both covalent and non-covalent inhibitors displaying protease IC50 values under 150 nM and viral inhibition under 5 uM in multiple different viral replication assays. Furthermore, we provide over 200 crystal structures of fragment-like and lead-like molecules in complex with the main protease. Over 1000 synthesized and ordered compounds are also reported with the corresponding activity in Mpro enzymatic assays using two different experimental setups. The data referenced in this document will be continually updated to reflect the current experimental progress of the COVID Moonshot project, and serves as a citable reference for ensuing publications. All of the generated data is open to other researchers who may find it of use.

Toxic dipeptide-repeat (DPR) proteins are produced from expanded G4C2 repeats in the C9ORF72 gene, the most common genetic cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Two DPR proteins, poly-PR and poly-GR, repress cellular translation but the molecular mechanism remains unknown. Here we show that poly-PR and poly-GR of ≥20 repeats inhibit the ribosome’s peptidyl-transferase activity at nanomolar concentrations, comparable to specific translation inhibitors. High-resolution cryogenic electron microscopy (cryo-EM) reveals that poly-PR and poly-GR block the polypeptide tunnel of the ribosome, extending into the peptidyl-transferase center (PTC). Consistent with these findings, the macrolide erythromycin, which binds in the tunnel, competes with poly-PR and restores peptidyl-transferase activity. Our results demonstrate that strong and specific binding of poly-PR and poly-GR in the ribosomal tunnel blocks translation, revealing the structural basis of their toxicity in C9ORF72-ALS/FTD.

Drugs that target the main protease (Mpro) of SARS-CoV-2 are effective therapeutics that have entered clinical use. Wide-scale use of these drugs will apply selection pressure for the evolution of resistance mutations. To understand resistance potential in Mpro, we performed comprehensive surveys of amino acid changes that can cause resistance in a yeast screen to nirmatrelvir (contained in the drug Paxlovid), and ensitrelvir (Xocova) that is currently in phase III trials. The most impactful resistance mutation (E166V) recently reported in multiple viral passaging studies with nirmatrelvir showed the strongest drug resistance score for nirmatrelvir, while P168R had the strongest resistance score for ensitrelvir. Using a systematic approach to assess potential drug resistance, we identified 142 resistance mutations for nirmatrelvir and 177 for ensitrelvir. Among these mutations, 99 caused apparent resistance to both inhibitors, suggesting a strong likelihood for the evolution of cross-resistance. Many mutations that exhibited inhibitor-specific resistance were consistent with distinct ways that each inhibitor protrudes beyond the substrate envelope. In addition, mutations with strong drug resistance scores tended to have reduced function. Our results indicate that strong pressure from nirmatrelvir or ensitrelvir will select for multiple distinct resistant lineages that will include both primary resistance mutations that weaken interactions with drug while decreasing enzyme function and secondary mutations that increase enzyme activity. The comprehensive identification of resistance mutations enables the design of inhibitors with reduced potential of developing resistance and aids in the surveillance of drug resistance in circulating viral populations.

Coronaviruses, as exemplified by SARS-CoV-2, can evolve and spread rapidly to cause severe disease morbidity and mortality. Direct acting antivirals (DAAs) are highly effective in decreasing disease burden especially when they target essential viral enzymes, such as proteases and polymerases, as demonstrated in HIV-1 and HCV and most recently SARS-CoV-2. Optimization of these DAAs through iterative structure-based drug design has been shown to be critical. Particularly, the evolutionarily conserved molecular mechanisms underlying viral replication can be leveraged to develop robust antivirals against rapidly evolving viral targets. The main protease (Mpro) of SARS-CoV-2, which is evolutionarily constrained to recognize and cleave 11 specific sites to promote viral maturation, exemplifies one such target. In this study we define the substrate envelope of Mpro by determining the molecular basis of substrate recognition, through nine high-resolution cocrystal structures of SARS-CoV-2 Mpro with the viral cleavage sites. These structures enable identification of evolutionarily vulnerable sites beyond the substrate envelope that may be susceptible to drug resistance and compromise binding of the newly developed Mpro inhibitors.

With the spread of SARS-CoV-2 throughout the globe causing the COVID-19 pandemic, the threat of zoonotic transmissions of coronaviruses (CoV) has become even more evident. As human infections have been caused by alpha- and beta-CoVs, structural characterization and inhibitor design mostly focused on these two genera. However, viruses from the delta and gamma genera also infect mammals and pose a potential zoonotic transmission threat. Here, we determined the inhibitor-bound crystal structures of the main protease (Mpro) from the delta-CoV porcine HKU15 and gamma-CoV SW1 from the beluga whale. A comparison with the apo structure of SW1 Mpro, which is also presented here, enabled the identification of structural arrangements upon inhibitor binding at the active site. The cocrystal structures reveal binding modes and interactions of two covalent inhibitors, PF-00835231 (active form of lufotrelvir) bound to HKU15, and GC376 bound to SW1 Mpro. These structures may be leveraged to target diverse coronaviruses and toward the structure-based design of pan-CoV inhibitors.

With the spread of SARS-CoV-2 throughout the globe to cause the COVID-19 pandemic, the threat of zoonotic transmissions of coronaviruses (CoV) has become even more evident. As human infections have been caused by alpha- and beta-CoVs, structural characterization and inhibitor design mostly focused on these two genera. However, viruses from the delta and gamma genera also infect mammals and pose potential zoonotic transmission threat. Here, we determined the inhibitor-bound crystal structures of the main protease (Mpro) from the delta-CoV porcine HKU15 and gamma-CoV SW1 from beluga whale. Comparison with the apo structure of SW1 Mpro, which we also present here, enabled identifying structural arrangements upon inhibitor binding at the active site. The binding modes and interactions of two covalent inhibitors, PF-00835231 (lufotrelvir) bound to HKU15 and GC376 bound to SW1 Mpro, reveal features that may be leveraged to target diverse coronaviruses and toward structure-based design of pan-CoV inhibitors.

No abstract