Crystal Structure of SARS-CoV-2 Main Protease in Complex with the Non-Covalent Inhibitor ML188

Abstract

Viral proteases are critical enzymes for the maturation of many human pathogenic viruses and thus are key targets for direct acting antivirals (DAAs). The current viral pandemic caused by SARS-CoV-2 is in dire need of DAAs. The Main protease (M^pro) is the focus of extensive structure-based drug design efforts which are mostly covalent inhibitors targeting the catalytic cysteine. ML188 is a non-covalent inhibitor designed to target SARS-CoV-1 M^pro, and provides an initial scaffold for the creation of effective pan-coronavirus inhibitors. In the current study, we found that ML188 inhibits SARS-CoV-2 M^pro at 2.5 µM, which is more potent than against SAR-CoV-1 M^pro. We determined the crystal structure of ML188 in complex with SARS-CoV-2 M^pro to 2.39 Å resolution. Sharing 96% sequence identity, structural comparison of the two complexes only shows subtle differences. Non-covalent protease inhibitors complement the design of covalent inhibitors against SARS-CoV-2 main protease and are critical initial steps in the design of DAAs to treat CoVID 19.

1. Introduction

Coronavirus disease 2019 (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) [1]. Zoonotic transmissions of coronaviruses to humans have caused at least three major outbreaks during the past two decades; SARS in 2002 [2,3], MERS in 2012 [4], and SARS-2 in 2019–2020 [5,6]. The previous outbreaks were contained with rigorous public health interventions, and basic research and vaccine/drug developments had largely been suspended. As is now evident from the current global pandemic [7], coronaviruses pose a major threat to human health, and treatment options need to be developed not only for the current pandemic but also for future coronavirus outbreaks. In addition to the major efforts for developing vaccines around the world, direct acting antivirals (DAAs) are critical to treat vulnerable patients to decrease morbidity and mortality.

Coronaviruses are non-segmented positive-sense single-stranded RNA viruses [8]. The viral genome encodes multiple open reading frames (ORF), including ORF1a and ORF1b, which are translated by the host machinery to generate two polyproteins (pp): pp1a and pp1ab [9]. These polyproteins are cleaved by two viral proteases, Main protease (M^pro) and the Papain-Like protease (PL^pro), liberating 16 individual proteins (nsp1–16). These proteolysis events are essential for the SARS-CoV-2 viral life cycle [10]. Thus, these two proteases represent excellent targets for the development of DAAs. Consistent with this notion, viral proteases are well-established targets for DAAs, with FDA-approved inhibitors as part of the mainstay of treatment for both HIV-1 and Hepatitis C virus infections. Given that the entire genomes of SARS-CoV-1 and SARS-CoV-2 are 79% identical [11], as well as the fact that SARS-CoV-2 M^pro (SARS2-M^pro) and SARS-CoV-2 PL^pro are 96% and 83% identical to their SARS-CoV-1 orthologs, respectively (Figure S1) [12], indicates that DAAs targeting either of these two proteases will likely show efficacy for a range of coronaviruses.

The current study focuses on SARS2-M^pro, a chymotrypsin-like cysteine protease that resides in the viral non-structural protein 5 (nsp5) [13]. M^pro cleaves 11 sites in pp1a and pp1ab, including autoproteolysis. The 4% difference in sequence with SARS1 corresponds to 12 changes in the 306-residue enzyme. Whereas most changes occur distal to the active site, we previously established that distal mutations can significantly affect inhibitor potency [14,15,16,17,18]. Potent, selective inhibition of SARS2-M^pro could be a successful therapeutic for Covid-19 and future coronaviruses. Structure-based drug design is already playing a pivotal role in targeting this viral enzyme. Fragment screens have placed hundreds of small fragments around SARS2-M^pro, in the active site and allosteric sites [19]. Continuing this quest, several publications include crystal structures of covalent inhibitors that target the catalytic cysteine [20,21,22], including one that has entered into clinical trials [23].

While covalent inhibitors have demonstrated high potency and efficacy, there is a lack of potent non-covalent inhibitors of coronavirus proteases. One potential scaffold is ML188 (Figure 1a), a non-covalent inhibitor designed to inhibit SARS1-M^pro [24]. ML188 was reported to bind SARS1-M^pro with an IC₅₀ of 1.5 ± 0.3 µM and an EC₅₀ of 12.9 ± 0.7 µM in cellular assays [24]. Pure ML188 has not been tested against SARS2-M^pro, but racemic ML188 inhibits with an IC₅₀ of 3.14 µM [25], which is comparable to racemic ML188 against SARS1-M^pro (4.8 ± 0.8 µM) [24]. Showing this scaffold’s versatility, ML188 was also found to inhibit M^pro from porcine epidemic diarrhea virus (PEDV) with only 17-fold less potency (IC₅₀ of 25.4 ± 1.4 μM), despite SARS1-M^pro and PEDV-M^pro sharing only 45.4% sequence identity [26]. This indicates that ML188 may provide a scaffold for a robust non-covalent pan-coronavirus inhibitor.

In this study, we characterize the complex of SARS-CoV-2 M^pro with the noncovalent inhibitor ML188. We determine that ML188 has enhanced binding potency to SARS-CoV-2 compared with SARS-CoV-1 and solved the cocrystal structure of ML188 in complex with SARS2-M^pro to 2.39 Å resolution. Overall, the complexes of M^pro with ML188 are very similar in both proteases, but subtle differences likely contribute to the higher potency of ML188 against SARS2-M^pro. 2. Materials and Methods 2.1. Protein Expression and Purification of SARS2-M^pro

The SARS2-M^pro plasmid, kindly provided by Rolf Hilgenfeld [21], was transformed into Escherichia coli strain HI-Control™ BL21(DE3) (Lucigen, Middleton, WI, USA). The transformed cells were pre-cultured at 37 °C in LB medium with ampicillin (100 μg/mL) overnight, and the cell culture was inoculated into TB medium containing 50 mM sodium phosphate (pH 7.0) and ampicillin (100 μg/mL). When OD₆₀₀ value reached ~2.0, 0.5 mM IPTG was added to induce SARS2-M^pro expression and the cell culture was further incubated overnight at 20 °C. Cells were harvested by centrifugation at 5000 rpm for 20 min, resuspended in lysis buffer (50 mM Tris–HCl (pH 8.0), 400 mM NaCl, 1 mM TCEP) and lysed by a cell disruptor. The lysate was clarified by ultracentrifugation at 18,000 rpm for 50 min. The supernatant was loaded onto a HisTrap FF column (Cytiva, Marlborough, MA, USA) equilibrated with lysis buffer, washed with lysis buffer and followed by elution using elution buffer (50 mM Tris–HCl pH 8.0, 400 mM NaCl, 500 mM imidazole, 1 mM TCEP) with a linear gradient of imidazole ranging from 0 mM to 500 mM. The fractions of M^pro-His tag were mixed with GST-PreScission protease-His-tag at a molar ratio of 5:1 to remove the C-terminal His tag. The PreScission-treated M^pro was applied to nickel column to remove the GST-PreScission protease-His-tag and protein with uncleaved His-tag. The His-tag cleaved M^pro in the flow-through was further purified by size-exclusion chromatography (HiLoad™ 16/60 Superdex 75 (Cytiva, Marlborough, MA, USA)) and stored in 20 mM HEPES pH 7.5, 150 mM NaCl, 1 mM TCEP.

2.2. M^Pro Inhibition Assay The M^Pro peptide substrate Dabcyl-KTSAVLQSGFRKM-E(Edans)-NH₂, was purchased from GenScript (Piscataway, NJ, USA). His-tagged SARS1-M^Pro was purchased from Sino Biological Inc. (Wayne, PA, USA). All assays were done in a 96-well half area plate (Corning, Corning, NY, USA). Peptide cleavage was measured using 50 nM enzyme. Assays were done in 50 mM HEPES pH 7.5, 150 mM NaCl, 1 mM EDTA, 1 mM DTT. SARS1-M^pro and SARS2-M^pro were incubated with either buffer or 0–50 µM of ML188 for 20 min. ML188 was purchased from MedChemExpress (Monmouth Junction, NJ, USA, CAT# HY-136259) with 98.35% purity. The reaction was initiated by adding 25 µM peptide substrate, followed by 30 min incubation at 25 °C. Fluorescence was measured at 485 nm with excitation at 340 nm with EnVision 2105 plate reader (Perkin Elmer, Waltham, MA, USA). Experiment was performed in duplicate and the error from global fit with variable hill slope to obtain IC₅₀ value is reported. 2.3. Protein Crystallization All crystallization screens tested provided conditions that produced M^pro cocrystals. A condition producing large crystals was discovered using the PACT Premier crystal screen (Molecular Dimensions, Maumee, OH, USA), Well E9, containing 20% (w/v) PEG 3350 and 0.2 M Potassium Sodium Tartrate Tetrahydrate. The SARS2-M^pro-ML188 cocrystal was grown at room temperature by hanging drop vapor diffusion method in a 24-well VDX hanging-drop tray (Hampton Research, Journey Aliso Viejo, CA, USA) with a protease concentration of 6.0 mg/mL with 6-fold molar excess of ML188 (10% DMSO) and mixed with the precipitant solution at a 1:1 ratio (1 µL:1 µL) and micro-seeded (1:100–1:10,000 dilution) with a cat whisker. Crystals appeared overnight and grew to diffraction quality after 3 days. As data was collected at 100 K, cryogenic conditions consisted of the precipitant solution supplemented with 25% glycerol. 2.4. Data Collection and Structure Determination

Diffraction quality crystals were flash frozen under a cryostream when mounted on our in-house Rigaku_Saturn944 X-ray system (Rigaku, The Woodlands, TX, USA). Co-crystal diffraction intensities were indexed, integrated, and scaled using HKL3000 [27]. The structure was solved using molecular replacement with PHASER [28] using an M^pro monomer (PDB: 6M03 by Zhang et al. DOI: 10.2210/pdb6M03/pdb). Model building and refinement was performed using Coot [29] and Phenix [30]. During refinement, optimized stereochemical weights were utilized. The ML188 ligand was designed in Maestro and the output SDF file was used in the Phenix program eLBOW [31] to generate the CIF file containing atomic positions and constraints necessary for ligand refinement. Iterative rounds of crystallographic refinement were carried out until convergence was achieved. To limit bias throughout the refinement process, five percent of the data was reserved for the free R-value calculation [32]. MolProbity [33] was applied to evaluate the final structure before deposition in the PDB [34,35]. Structure analysis, superposition and figure generation was done using PyMOL [36]. X-ray data collection and crystallographic refinement statistics are presented in the Supporting Information (Table S1).

2.5. Intermolecular vdW Contact Analysis of Crystal Structures

To calculate the intermolecular vdW interaction energies, the crystal structures were prepared using the Schrödinger Protein Preparation Wizard [37]. Hydrogen atoms were added, protonation states were determined, and the structures were minimized. Subsequently, force field parameters were assigned using the OPLS3 force field [38]. Interaction energies between the inhibitor and protease were estimated using a simplified Lennard-Jones potential, as previously described in detail [39]. Briefly, the vdW energy was calculated for pairwise interactions depending on the types of atoms interacting and the distance between them.

3. Results

To compare the potency of ML188, we purchased SARS1-M^pro and expressed and purified SARS2-M^pro for enzyme inhibition assays. Using a FRET-based enzymatic assay, ML188 inhibits SARS1-M^pro with an IC₅₀ of 4.5 ± 0.5 µM and inhibits SARS2-M^pro with an IC₅₀ of 2.5 ± 0.3 µM. Therefore, ML188 is approximately twice as potent in SARS2 compared to SARS1 (Figure 1b).

The SARS-M^pro is a functional homodimer (Figure 2A) that binds and cleaves 11 different substrates. The protease has well defined subsites at P1 and P2, which recognize glutamine and large hydrophobic residues, respectively. To elucidate structural differences in ML188 binding, the cocrystal structure of ML188 in complex with SARS2-M^pro was determined to 2.39 Å resolution (Table S1; Figure 2). The SARS2-M^pro-ML188 complex was solved in the same space group (C2) as the SARS1-M^pro-ML188 complex (PDB: 3V3M). Both structures contain one M^pro subunit in the asymmetric unit, but the complexes had different cell dimensions with vastly different crystal packing (Figure S2). The electron density permitted ML188 to be modeled unambiguously in the same orientation as found in the SARS1-M^pro-ML188 complex (Figure 2D). Overall, ML188 binds similarly in both complexes and makes all the same interactions (Figure 2E), with subtle differences.

Aligning the SARS1-M^pro-ML188 and SARS2-M^pro-ML188 complexes in PyMOL finds the protease alpha carbons have a root-mean-square deviation (RMSD) of 0.53 Å. The largest relevant C-α backbone differences occurred at residues 45–49 in the loop above the S2 subsite (Figure S3). This may relate to the substitution in this loop and the overall 12-residue difference between SARS1 and SARS2 M^pro (Figure 3A and Figure S1). The mutation closest to the active site is A46S (Figure 3B,C), which was found to affect the dynamics of that pocket [12]. Changes in protease dynamics can affect the overall binding potency of inhibitors.

The structures were used to calculate protease-inhibitor van der Waals (vdW) contacts that showed minor differences in interactions, SARS1-M^pro-ML188 with −59.3 kcal/mol versus SARS2-M^pro-ML188 with −62.0 kcal/mol. This difference of 2.3 kcal/mol is directly due to the conformers modeled into the electron density of three side chains, specifically M49, N142, and Q189 (Figure S4). The vdW calculations also show that ML188 interacts more with the catalytic histidine H41. All other interactions are within the margin of error for structures of this resolution, but the overall trends suggest ML188 packs deeper into the S2 subsite of SARS2-M^pro, but less deep into the S1 subsite compared to SARS1-M^pro (Figure S4).

4. Conclusions In this study, we determine that the non-covalent inhibitor ML188 has enhanced binding potency to SARS2-M^pro compared to SARS1-M^pro. We characterize the complex of SARS2-M^pro with ML188 and compare it to the SARS1-M^pro-ML188 complex. Overall, the complexes are very similar in both proteases, but subtle differences likely contribute to the higher potency of ML188 against SARS2-M^pro. This indicates ML188 may provide a scaffold for a robust non-covalent pan-coronavirus inhibitor.

Figure 1.(a) Chemical structure of ML188. (b) Dose-response curves and IC50 values of compound ML188 against SARS2-MPro and SARS1-MPro.

Figure 2.(A) SARS2-Mpro dimer and single subunit shown in surface representation (rotated on the plane of the page). (B) SARS2-Mpro active site with subsites labeled and catalytic residues colored yellow. (C) ML188 in the SARS2-Mpro active site with catalytic residues labeled. (D) Electron density around ML188 and catalytic residues. The 2Fo-Fc direct maps are depicted as grey mesh contoured at 1.0 σ while the Fo-Fc difference maps have positive density depicted as green mesh contoured at 3.0 σ and negative density as red mesh contoured at -3.0 σ. (E) Alignment of ML188 in SARS2-Mpro (teal protein and green ligand; PDB: 7L0D) and SARS1-Mpro (magenta protein and orange ligand; PDB: 3V3M).

Figure 3.(A) The twelve amino acid differences between SARS2 and SARS1-Mpro shown as a dimer, single subunit, and active site. SARS2 (in cyan) and SARS1 (in magenta) are shown with transparent surfaces and amino acid differences shown as sticks. The dimer subunits (in orange) and the catalytic residues (in yellow) are included. (B) Stereo pair of SARS2-Mpro-ML188 complex active site. (C) Stereo pair of SARS1-Mpro-ML188 complex active site.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4915/13/2/174/s1, Figure S1: Sequence alignment of SARS2-Mpro and SARS1-Mpro. Table S1: X-ray data collection and crystallographic refinement statistics. Figure S2: Comparison of crystallographic symmetry mates. Figure S3: C-alpha distance differences between SARS2-Mpro ML188 complex and SARS1-Mpro ML188 complex. Figure S4: Protease-Inhibitor vdW differences.

Author Contributions

Conceptualization, G.J.L. and A.A.; Investigation, G.J.L., J.M.L., R.T., E.A.N. and A.C.R.; Data Curation, G.J.L. and A.C.R.; Writing-Original Draft Preparation, G.J.L.; Writing-Review and Editing, N.K.Y., C.A.S. and P.R.T.; Visualization, G.J.L.; Supervision, N.K.Y., C.A.S. and P.R.T.; Funding Acquisition, P.R.T. and C.A.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Institute of General Medical Sciences, grant numbers R01 GM135919 and R35 GM118112; National Institute of Allergy and Infectious Diseases grant number R21 AI149716.

Data Availability Statement

The structural data presented in this study are openly available in the Protein Data Bank, identification number 7L0D. Inhibition data presented in this study are available on request from the corresponding author.

Acknowledgments

This manuscript was written in honor of Stephan Oroszlan, with his expertise in retroviruses, realizing the significance viral protease (or proteinases) early on as therapeutic targets [40]. In this quest he characterized their complex substrate specificity for HTLV-1 [41], HIV-1 [42,43,44,45], and HIV-2 [46,47] and many other retroviral proteases. As we now face a new global viral pandemic, characterizing and targeting the viral proteases would be a clear pathway forward in the legacy of Oroszlan.

Conflicts of Interest

The authors declare no conflict of interest.

References

1. Guan, W.J.; Ni, Z.Y.; Hu, Y.; Liang, W.H.; Ou, C.Q.; He, J.X.; Liu, L.; Shan, H.; Lei, C.L.; Hui, D.S.C.; et al. China Medical Treatment Expert Group for C. Clinical Characteristics of Coronavirus Disease 2019 in China. N. Engl. J. Med. 2020.

2. Drosten, C.; Gunther, S.; Preiser, W.; van der Werf, S.; Brodt, H.R.; Becker, S.; Rabenau, H.; Panning, M.; Kolesnikova, L.; Fouchier, R.A.; et al. Identification of a novel coronavirus in patients with severe acute respiratory syndrome. N. Engl. J. Med. 2003, 348, 1967-1976.

3. Ksiazek, T.G.; Erdman, D.; Goldsmith, C.S.; Zaki, S.R.; Peret, T.; Emery, S.; Tong, S.; Urbani, C.; Comer, J.A.; Lim, W.; et al. A novel coronavirus associated with severe acute respiratory syndrome. N. Engl. J. Med. 2003, 348, 1953-1966.

4. Zaki, A.M.; van Boheemen, S.; Bestebroer, T.M.; Osterhaus, A.D.; Fouchier, R.A. Isolation of a novel coronavirus from a man with pneumonia in Saudi Arabia. N. Engl. J. Med. 2012, 367, 1814-1820.

5. Wu, F.; Zhao, S.; Yu, B.; Chen, Y.M.; Wang, W.; Song, Z.G.; Hu, Y.; Tao, Z.W.; Tian, J.H.; Pei, Y.Y.; et al. A new coronavirus associated with human respiratory disease in China. Nature 2020, 579, 265-269.

6. Zhu, N.; Zhang, D.; Wang, W.; Li, X.; Yang, B.; Song, J.; Zhao, X.; Huang, B.; Shi, W.; Lu, R.; et al. A Novel Coronavirus from Patients with Pneumonia in China, 2019. N. Engl. J. Med. 2020, 382, 727-733.

7. Gabriel, S.E. The epidemiology of rheumatoid arthritis. Rheum. Dis. Clin. N. Am. 2001, 27, 269-281.

8. Fehr, A.R.; Perlman, S. Coronaviruses: An overview of their replication and pathogenesis. Methods Mol. Biol. 2015, 1282, 1-23.

9. Su, S.; Wong, G.; Shi, W.; Liu, J.; Lai, A.C.K.; Zhou, J.; Liu, W.; Bi, Y.; Gao, G.F. Epidemiology, Genetic Recombination, and Pathogenesis of Coronaviruses. Trends. Microbiol. 2016, 24, 490-502.

10. Zumla, A.; Chan, J.F.; Azhar, E.I.; Hui, D.S.; Yuen, K.Y. Coronaviruses-Drug discovery and therapeutic options. Nat. Rev. Drug Discov. 2016, 15, 327-347.

11. Lu, R.; Zhao, X.; Li, J.; Niu, P.; Yang, B.; Wu, H.; Wang, W.; Song, H.; Huang, B.; Zhu, N.; et al. Genomic characterisation and epidemiology of 2019 novel coronavirus: Implications for virus origins and receptor binding. Lancet 2020, 395, 565-574.

12. Bzowka, M.; Mitusinska, K.; Raczynska, A.; Samol, A.; Tuszynski, J.A.; Gora, A. Structural and Evolutionary Analysis Indicate That the SARS-CoV-2 Mpro Is a Challenging Target for Small-Molecule Inhibitor Design. Int. J. Mol. Sci. 2020, 21, 3099.

13. Wang, H.; He, S.; Deng, W.; Zhang, Y.; Li, G.; Sun, J.; Zhao, W.; Guo, Y.; Yin, Z.; Li, D.; et al. Comprehensive Insights into the Catalytic Mechanism of Middle East Respiratory Syndrome 3C-Like Protease and Severe Acute Respiratory Syndrome 3C-Like Protease. ACS Catal. 2020, 10, 5871-5890.

14. Ragland, D.A.; Nalivaika, E.A.; Nalam, M.N.; Prachanronarong, K.L.; Cao, H.; Bandaranayake, R.M.; Cai, Y.; Kurt-Yilmaz, N.; Schiffer, C.A. Drug resistance conferred by mutations outside the active site through alterations in the dynamic and structural ensemble of HIV-1 protease. J. Am. Chem. Soc. 2014, 136, 11956-11963.

15. Ragland, D.A.; Whitfield, T.W.; Lee, S.K.; Swanstrom, R.; Zeldovich, K.B.; Kurt-Yilmaz, N.; Schiffer, C.A. Elucidating the Interdependence of Drug Resistance from Combinations of Mutations. J. Chem. Theory Comput. 2017, 13, 5671-5682.

16. Henes, M.; Lockbaum, G.J.; Kosovrasti, K.; Leidner, F.; Nachum, G.S.; Nalivaika, E.A.; Lee, S.K.; Spielvogel, E.; Zhou, S.; Swanstrom, R.; et al. Picomolar to Micromolar: Elucidating the Role of Distal Mutations in HIV-1 Protease in Conferring Drug Resistance. ACS Chem. Biol. 2019, 14, 2441-2452.

17. Leidner, F.; Kurt Yilmaz, N.; Schiffer, C.A. Target-Specific Prediction of Ligand Affinity with Structure-Based Interaction Fingerprints. J. Chem. Inf. Model. 2019, 59, 3679-3691.

18. Whitfield, T.W.; Ragland, D.A.; Zeldovich, K.B.; Schiffer, C.A. Characterizing Protein-Ligand Binding Using Atomistic Simulation and Machine Learning: Application to Drug Resistance in HIV-1 Protease. J. Chem. Theory Comput. 2020, 16, 1284-1289.

19. Douangamath, A.; Fearon, D.; Gehrtz, P.; Krojer, T.; Lukacik, P.; Owen, C.D.; Resnick, E.; Strain-Damerell, C.; Aimon, A.; Abranyi-Balogh, P.; et al. Crystallographic and electrophilic fragment screening of the SARS-CoV-2 main protease. Nat. Commun. 2020, 11, 5047.

20. Hoffman, R.L.; Kania, R.S.; Brothers, M.A.; Davies, J.F.; Ferre, R.A.; Gajiwala, K.S.; He, M.; Hogan, R.J.; Kozminski, K.; Li, L.Y.; et al. Discovery of Ketone-Based Covalent Inhibitors of Coronavirus 3CL Proteases for the Potential Therapeutic Treatment of COVID-19. J. Med. Chem. 2020, 63, 12725-12747.

21. Zhang, L.; Lin, D.; Sun, X.; Curth, U.; Drosten, C.; Sauerhering, L.; Becker, S.; Rox, K.; Hilgenfeld, R. Crystal structure of SARS-CoV-2 main protease provides a basis for design of improved alpha-ketoamide inhibitors. Science 2020, 368, 409-412.

22. Dai, W.; Zhang, B.; Jiang, X.M.; Su, H.; Li, J.; Zhao, Y.; Xie, X.; Jin, Z.; Peng, J.; Liu, F.; et al. Structure-based design of antiviral drug candidates targeting the SARS-CoV-2 main protease. Science 2020, 368, 1331-1335.

23. Boras, B.; Jones, R.M.; Anson, B.J.; Arenson, D.; Aschenbrenner, L.; Bakowski, M.A.; Beutler, N.; Binder, J.; Chen, E.; Eng, H.; et al. Discovery of a Novel Inhibitor of Coronavirus 3CL Protease as a Clinical Candidate for the Potential Treatment of COVID-19. bioRxiv 2020.

24. Jacobs, J.; Grum-Tokars, V.; Zhou, Y.; Turlington, M.; Saldanha, S.A.; Chase, P.; Eggler, A.; Dawson, E.S.; Baez-Santos, Y.M.; Tomar, S.; et al. Discovery, synthesis, and structure-based optimization of a series of N-(tert-butyl)-2-(N-arylamido)-2-(pyridin-3-yl) acetamides (ML188) as potent noncovalent small molecule inhibitors of the severe acute respiratory syndrome coronavirus (SARS-CoV) 3CL protease. J. Med. Chem. 2013, 56, 534-546.

25. Zaidman, D.; Gehrtz, P.; Filep, M.; Fearon, D.; Prilusky, J.; Duberstein, S.; Cohen, G.; Owen, D.; Resnick, E.; Strain-Damerell, C.; et al. An automatic pipeline for the design of irreversible derivatives identifies a potent SARS-CoV-2 Mpro inhibitor. bioRxiv 2020.

26. St John, S.E.; Anson, B.J.; Mesecar, A.D. X-Ray Structure and Inhibition of 3C-like Protease from Porcine Epidemic Diarrhea Virus. Sci. Rep. 2016, 6, 25961.

27. Otwinowski, Z.; Minor, W. Processing of X-ray diffraction data collected in oscillation mode. Methods Enzymol. 1997, 276, 307-326.

28. McCoy, A.J.; Grosse-Kunstleve, R.W.; Adams, P.D.; Winn, M.D.; Storoni, L.C.; Read, R.J. Phaser crystallographic software. J. Appl. Crystallogr. 2007, 40 Pt 4, 658-674.

29. Emsley, P.; Cowtan, K. Coot: Model-building tools for molecular graphics. Acta Crystallogr. D Biol. Crystallogr. 2004, 60 Pt 12, 2126-2132.

30. Adams, P.D.; Afonine, P.V.; Bunkoczi, G.; Chen, V.B.; Davis, I.W.; Echols, N.; Headd, J.J.; Hung, L.W.; Kapral, G.J.; Grosse-Kunstleve, R.W.; et al. PHENIX: A comprehensive Python-based system for macromolecular structure solution. Acta Crystallogr. D Biol. Crystallogr. 2010, 66 Pt 2, 213-221.

31. Moriarty, N.W.; Grosse-Kunstleve, R.W.; Adams, P.D. electronic Ligand Builder and Optimization Workbench (eLBOW): A tool for ligand coordinate and restraint generation. Acta Crystallogr. D Biol. Crystallogr. 2009, 65 Pt 10, 1074-1080.

32. Brunger, A.T. Free R value: A novel statistical quantity for assessing the accuracy of crystal structures. Nature 1992, 355, 472-475.

33. Davis, I.W.; Leaver-Fay, A.; Chen, V.B.; Block, J.N.; Kapral, G.J.; Wang, X.; Murray, L.W.; Arendall, W.B., 3rd; Snoeyink, J.; Richardson, J.S.; et al. MolProbity: All-atom contacts and structure validation for proteins and nucleic acids. Nucleic Acids Res. 2007, 35, W375-W383.

34. Berman, H.M.; Westbrook, J.; Feng, Z.; Gilliland, G.; Bhat, T.N.; Weissig, H.; Shindyalov, I.N.; Bourne, P.E. The Protein Data Bank. Nucleic Acids Res. 2000, 28, 235-242.

35. Berman, H.; Henrick, K.; Nakamura, H. Announcing the worldwide Protein Data Bank. Nat. Struct. Biol. 2003, 10, 980.

36. DeLano, W.L. Pymol: An open-source molecular graphics tool. CCP4 Newsl. Protein. Crystallogr. 2002, 40, 82-92.

37. Schrödinger. Schrödinger Release 2019-2; Schrödinger, LLC: New York, NY, USA, 2019.

38. Harder, E.; Damm, W.; Maple, J.; Wu, C.; Reboul, M.; Xiang, J.Y.; Wang, L.; Lupyan, D.; Dahlgren, M.K.; Knight, J.L.; et al. OPLS3: A Force Field Providing Broad Coverage of Drug-like Small Molecules and Proteins. J. Chem. Theory Comput. 2016, 12, 281-296.

39. Nalam, M.N.; Ali, A.; Altman, M.D.; Reddy, G.S.; Chellappan, S.; Kairys, V.; Ozen, A.; Cao, H.; Gilson, M.K.; Tidor, B.; et al. Evaluating the Substrate-Envelope Hypothesis: Structural Analysis of Novel HIV-1 Protease Inhibitors Designed to be Robust against Drug Resistance. J. Virol. 2010, 84, 368-378.

40. Tozser, J.; Oroszlan, S. Proteolytic events of HIV-1 replication as targets for therapeutic intervention. Curr. Pharm. Des. 2003, 9, 1803-1815.

41. Tozser, J.; Zahuczky, G.; Bagossi, P.; Louis, J.M.; Copeland, T.D.; Oroszlan, S.; Harrison, R.W.; Weber, I.T. Comparison of the substrate specificity of the human T-cell leukemia virus and human immunodeficiency virus proteinases. Eur. J. Biochem. 2000, 267, 6287-6295.

42. Tozser, J.; Bagossi, P.; Weber, I.T.; Louis, J.M.; Copeland, T.D.; Oroszlan, S. Studies on the symmetry and sequence context dependence of the HIV-1 proteinase specificity. J. Biol. Chem. 1997, 272, 16807-16814.

43. Tozser, J.; Weber, I.T.; Gustchina, A.; Blaha, I.; Copeland, T.D.; Louis, J.M.; Oroszlan, S. Kinetic and modeling studies of S3-S3' subsites of HIV proteinases. Biochemistry 1992, 31, 4793-4800.

44. Roberts, M.M.; Copeland, T.D.; Oroszlan, S. In situ processing of a retroviral nucleocapsid protein by the viral proteinase. Protein Eng. 1991, 4, 695-700.

45. Tozser, J.; Gustchina, A.; Weber, I.T.; Blaha, I.; Wondrak, E.M.; Oroszlan, S. Studies on the role of the S4 substrate binding site of HIV proteinases. FEBS Lett. 1991, 279, 356-360.

46. Bagossi, P.; Cheng, Y.S.; Oroszlan, S.; Tozser, J. Comparison of the specificity of homo- and heterodimeric linked HIV-1 and HIV-2 proteinase dimers. Protein Eng. 1998, 11, 439-445.

47. Tozser, J.; Blaha, I.; Copeland, T.D.; Wondrak, E.M.; Oroszlan, S. Comparison of the HIV-1 and HIV-2 proteinases using oligopeptide substrates representing cleavage sites in Gag and Gag-Pol polyproteins. FEBS Lett. 1991, 281, 77-80.