Open AccessArticle

In Silico Discovery of Small-Molecule Inhibitors Targeting SARS-CoV-2 Main Protease

School of Pharmacy and Pharmaceutical Sciences & Institute of Materia Medica, Shandong First Medical University & Shandong Academy of Medical Sciences, 6699 Qingdao Road, Jinan 250117, China

NHC Key Laboratory of Biotechnology Drugs, Shandong Academy of Medical Sciences, 6699 Qingdao Road, Jinan 250117, China

Key Lab for Rare & Uncommon Diseases of Shandong Province, 6699 Qingdao Road, Jinan 250117, China

⁴

Department of Medicinal Chemistry, School of Pharmaceutical Sciences, Cheeloo College of Medicine, Shandong University, 44 West Culture Road, Jinan 250012, China

Authors to whom correspondence should be addressed.

Molecules 2023, 28(14), 5320; https://doi.org/10.3390/molecules28145320

Submission received: 30 March 2023 / Revised: 29 June 2023 / Accepted: 5 July 2023 / Published: 10 July 2023

(This article belongs to the Special Issue Computational Studies on the Development and Characterization of Pharmaceutical Materials)

Download

Browse Figures

Figure 1
Life cycle of SARS-CoV-2. "> Figure 2
Crystal structure of Mpro homodimer and the catalytic site. Ligand: green stick/sphere; key residues: pink stick; hydrogen bond: yellow dashed line; π-π stacking: purple dashed line; PDB code: 7KX5 [<a href="#B18-molecules-28-05320" class="html-bibr">18</a>]. "> Figure 3
Representative peptidic covalent Mpro inhibitors. "> Figure 4
Representative non-peptidic Mpro inhibitors. "> Figure 5
The flow chart of the discovery of inhibitors against SARS-CoV-2 Mpro via docking-based virtual screening. "> Figure 6
The RMSD and RMSF plots of AF-399/40713777 and AI-942/42301830 binding to Mpro generated via MD simulation. "> Figure 7
The protein–ligand binding conformations generated from the molecular dynamics simulation trajectory. Ligand: green stick; key residues: pink stick; hydrogen bond: yellow dashed line; π-π stacking: purple dashed line. The structure files (PDB files) can be found in the <a href="#app1-molecules-28-05320" class="html-app">Supplementary Materials</a>. "> Figure 8
The ADMET profile prediction of small molecules processed using ADMET lab 2.0. "> Figure 9
The toxicity prediction of representative ligands processed using the ProTox-II server. ">

Review Reports Versions Notes

Abstract

The COVID-19 pandemic has caused severe health threat globally, and novel SARS-Cov-2 inhibitors are urgently needed for antiviral treatment. The main protease (M^pro) of the virus is one of the most effective and conserved targets for anti-SARS-CoV-2 drug development. In this study, we utilized a molecular docking-based virtual screening approach against the conserved catalytic site to identify small-molecule inhibitors of SARS-CoV-2 M^pro. Further biological evaluation helped us identify two compounds, AF-399/40713777 and AI-942/42301830, with moderate inhibitory activity. Besides that, the in silico data, including molecular dynamics (MD) simulation, binding free energy calculations, and AMDET profiles, suggested that these two hits could serve as the starting point for the future development of COVID-19 intervention treatments.

Keywords:

SARS-CoV-2; main protease inhibitors; molecular docking-based virtual screening; MD simulation

1. Introduction

The COVID-19 pandemic, caused by SARS-CoV-2, has been a worldwide health emergency since 2019. SARS-CoV-2, which belongs to the Coronaviridae family, has a single-stranded positive RNA genome. The infection of SARS-CoV-2 could cause pneumonia, fatigue, diarrhea, and life-threatening cardiovascular complications or multiorgan failure. Two antiviral drugs (Remdesivir [1] and Nirmatrelvir [2]) have been approved by the FDA for emergency use. But further clinical data have shown that Remdesivir only has weak potency against SARS-CoV-2 infection [3], and the application of Nirmatrelvir can only prevent severe symptoms of the infection [4]. Novel anti-SARS-CoV-2 agents are still urgently needed.

The replication of SARS-CoV-2 can be divided into several steps, including fusion and entry, translation, polyprotein processing, RNA replication, viral assembly, and particle release (Figure 1). All of these steps are catalyzed by various enzymes of the virus (protease, RNA-dependent RNA polymerase, etc.) [5,6]. Theoretically, the inhibition of any enzyme could stop the replication of the virus. The fusion and entry of SARS-CoV-2 rely on the recognition of viral spike (S) protein and the host ACE2 (angiotensin-converting enzyme 2) [7]. The S protein is assembled as a homotrimer and is inserted in multiple copies into the phospholipid membrane of the virion, which gives it a crown-like appearance. The S protein consists of two non-covalently associated subunits: the S1 subunit binds ACE2 and the S2 subunit anchors the S protein to the viral membrane. The S2 subunit also includes a fusion peptide and other machinery necessary to mediate membrane fusion upon the infection of a new host cell. The receptor-binding domain (RBD) of the S1 subunit flips to an up conformation and binds to ACE2 [8]. Once the viral particle has fused and entered the cell, the viral positive-sense, single-stranded RNA genome (+ssRNA) is then exposed and translated into its polyprotein. The polyprotein is then cleaved into 16 non-structural proteins by two cysteine proteases, including the main protease (M^pro) and papain-like protease (PL^pro) [9,10]. M^pro cleaves the majority of non-structural proteins, which is crucial for viral replication, making it an ideal target for anti-viral agents. The viral genomic replication starts from the synthesis of full-length negative-sense RNA genome, which is the template for further production of +ssRNA. The newly produced +ssRNA is then used for the translation of polyprotein or for the packaging of new viral particles. The RNA replication relies on the replication and transcription complex (RTC), which contains RNA-dependent RNA polymerase (RdRp) and its cofactors, nsp7 and nsp8 [11,12]. Remdesivir, developed by Gilead, is approved by the FDA as novel nucleoside inhibitor of RdRp [1,3,13]. However, Remdesivir shows poor potency in inhibiting viral replication. Further investigation illustrated that RdRp possesses an in trans backtracking proofreading mechanism, which makes the inhibition potency of Remdesivir decreases badly [14]. The assembly and particle release of SARS-CoV-2 are also mediated by several enzymes from the virus and the host cell. The structural proteins (S protein, envelope (E) protein, membrane (M) protein, and nucleocapsid (N) protein) are translated, and the viral particle is assembled at the endoplasmic reticulum-to-Golgi compartment (ERGIC), where the viral RNA genome is encapsulated by the N protein and interacts with the M protein [15]. The M protein forms a scaffold for the virion and recruits the E and S proteins to the budding site. The E protein is a small membrane protein that regulates the membrane curvature and stability of the virion [16]. The S protein is a large glycoprotein that mediates the attachment and entry of the virus into the host cell by binding to the ACE2 receptor. The S protein undergoes proteolytic cleavage and conformational changes to expose its fusion peptide and facilitate membrane fusion. After assembly, the mature viral particles are transported from the ERGIC to the plasma membrane via vesicles. These vesicles fuse with the host cell membrane and release the virions via exocytosis [5,17].

M^pro is one of the most important enzymes that processes viral polyprotein into functional proteins. M^pro is a homodimer, where each monomer can be divided into two domains, the functional catalytic domain and the crucial dimerization domain. The catalytic site of M^pro consists of Cys145 and His41. The whole binding pocket can be identified as several sub-pockets (S1, S1’, S2, and S4). Occupation of the catalytic site could inhibit the normal function of M^pro and, thus, inhibit the replication of virus. The endogenous binding of ligands and the vital roles of M^pro make it an ideal target for drug design and development [17]. A shown in Figure 2, the inhibitor reported by Ma et al. [18] binds the catalytic site of M^pro and forms multiple hydrogen bonds with surrounding residues (Glu166, His163, Ser144, and Gly143). It also forms π-π stacking with His 41. The cocrystal structure of the M^pro–inhibitor complex can shed light on the development of novel M^pro inhibitors.

To date, several ligands targeting M^pro have been reported (Figure 3) [19,20,21]. Nirmatrelvir, developed by Pfizer, acts as an orally active ligand in combination with the CYP3A inhibitor Ritonavir [4,22]. GC-376 [23], PF-00835231 [24,25], and 11a [26] are currently under clinical trials. All of these peptidic inhibitors simulate the endogenous substrate of M^pro and possess warheads which could covalently bind with the vital catalytic residue Cys145.

Despite the high potency of peptidic covalent M^pro inhibitors, many of them show off-target side effects on host proteins. GC-376 shows inhibitory activity toward several cathepsins (IC₅₀ = 990, 74, and 0.56 nM to cathepsin L, cathepsin I, and cathepsin K, respectively). PF-00835231 also shows inhibitory activity toward several cathepsins (IC₅₀ = 146 nM and 1300 μM to cathepsin L and cathepsin B, respectively). 11a is potent inhibitor of cathepsin L (IC₅₀ = 210 nM) [18].

Compared to peptidic covalent M^pro inhibitors, non-peptidic inhibitors have attracted researchers’ attention due to their advantages, such as low off-target activity, higher oral bioavailability, and low toxicity. The non-peptidic inhibitors, 21, 23R, S-216722, and Wu04, display high antiviral activity and low toxicity (Figure 4) [27,28,29,30]. All of these findings suggest that non-peptidic M^pro inhibitors can be great potential candidate drugs for treating SARS-CoV-2 infection.

Computer-aided drug design (CADD) has been a general and practical technique highly used in the development of novel drugs; it has become a fundamental tool in the development of novel drugs for most pharmaceutical companies [31,32]. Several blockbuster drugs (Glivec, etc.) were developed using CADD [33]. In this context, virtual screening has been used to discover novel chemotypes of SARS-CoV-2 M^pro inhibitors. Approaches based on docking virtual screening using co-crystal structures, pharmacophore models, and quantitative structure–activity relationship (QSAR)-based screening are used to identify novel hits and guide the structure optimization of M^pro inhibitors. Moreover, molecular dynamics simulations have been exploited to investigate the binding mechanisms of these inhibitors [34,35].

Here, we report 17 compounds identified from molecular docking-based virtual screening targeting M^pro. The initial biological evaluations identified two novel scaffolds that possessed moderate potency toward M^pro. Molecular dynamics (MD) simulations of representative ligands targeting M^pro were performed to investigate the binding conformations of these ligands. The ADMET profiles of these ligands were also estimated to predict the potential draggability of representative molecules.

2. Results and Discussion

2.1. Molecular Docking-Based Virtual Screening

The cocrystal structure of M^pro with 23R (PDB ID: 7KX5) was selected as the structure complex since few cocrystal structures had been reported when this work started. The binding pocket was identified using the binding pose of 23R. Two commercial compound libraries (Enamine and Specs, total compound number > 3.5 M) were used for selection. A 3-step molecular docking-based virtual screening and a binding free energy (Molecular Mechanics/Generalized Born Surface Area, MM/GBSA) calculation were carried out. Finally, 17 compounds were selected based on the docking scores, MM/GBSA estimation, and manual selection (Figure 5, Table 1). The manual selection was based on the criteria reported by Fischer et al. [36], including fundamental interaction with key residues, binding conformation, and structural novelty.

2.2. Initial Biological Evaluation

These selected candidates were then purchased, and their enzyme inhibitory activities were evaluated using a fluorescence resonance energy transfer (FRET) assay (Table 2). Notably, AF-399/40713777 and AI-942/42301830 possess 47.11% and 37.67% of inhibitory activity toward M^pro at 100 μM, respectively. Moreover, AN-329/15538195, AN-655/14907067, AG-690/13705944, AK-968/37129380, and Z929753284 show >10% inhibition rate toward M^pro. These encouraging biological results indicate these diverse scaffolds can be used for further SAR studies. These novel scaffolds provide more possibility for developing good candidates for treating SARS-CoV-2 infection.

2.3. Molecular Dynamics (MD) Simulation

MD simulation is a powerful computational technique that can provide insights into the structure, dynamics, and interactions of biomolecules at the atomic level. The initial biological results suggest that AF-399/40713777 and AI-942/42301830 possess moderate inhibitory activity toward SARS-CoV-2 M^pro. The binding conformations of these two ligands with the pocket were then investigated. The MD simulation results show that both AF-399/40713777 and AI-942/42301830 bind to the catalytic site of M^pro. The RMSD plots of the ligands heavy atoms and protein backbones were generated (Figure 6). The RMSF values of the ligands were also estimated (Figure 6). The computer data indicate that both of AF-399/40713777 and AI-942/42301830 have stable binding conformation with M^pro, as shown through the 500 ns MD simulation.

The protein–ligand interaction diagram was further analyzed (Figure 7). AF-399/40713777 forms multiple interactions with surrounding residues. The carbonyl motifs interact with Gln189, Gln192, and Thr190 via hydrogen bond. The methoxy group forms hydrogen bond with Arg188. Besides these interactions, the phenyl ring forms π-π stacking with His41. These multiple interactions stabilize the binding conformation between the ligand and the pocket.

AI-942/42301830 mainly binds to M^pro through abundant hydrophobic interactions. His41 forms multiple π-π stacking with the two phenyl rings of the ligand. It also interacts with Met49, Leu167, and Pro168 via hydrophobic interactions. Moreover, the carbonyl motif of ester forms hydrogen bond with Glu166.

2.4. Binding Free Energy (MM/GBSA) Calculations

The binding free energy (MM/GBSA) [37] of AF-399/40713777 (−66.07 kcal/mol) and AI-942/42301830 (−78.32 kcal/mol) with SARS-CoV-2 M^pro were estimated. The calculated binding affinities did translate well into experimentally determined inhibitory potencies (Table 3). As illustrated in Figure 7, AF-399/40713777 forms multiple hydrogen bonds with surrounding residues. And the ΔG_Hbond value of AF-399/40713777 is relatively higher compared to that of AI-942/42301830. On the other hand, AI-942/42301830 forms abundant hydrophobic interactions with the binding pocket. The ΔG_Lipo of AI-942/42301830 contributes more compared to that of AF-399/40713777.

2.5. Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADMET) Prediction

ADMET is an important aspect of drug discovery and development as it determines the efficacy, safety, and optimal dosage of a drug. A poor ADMET profile is one of the main reasons for drug failure in clinical trials, resulting in wasted time, resources, and money. In order to eliminate compounds with undesirable pharmacokinetics and toxicity and reduce the risk of new drug discovery, we also evaluated the AMDET profiles of these two promising hits, AF-399/40713777 and AI-942/42301830, using ADMETlab 2.0 (Figure 8) and the ProTox-II server (Figure 9) [38,39]. ADMETlab 2.0 [38] is a web-based platform that provides comprehensive and accurate predictions of ADMET properties of molecules. ADMETlab 2.0 enables users to evaluate 17 physicochemical properties, 13 medicinal chemistry measures, 23 ADME endpoints, 27 toxicity endpoints, and 8 toxicophore rules of input molecules, using a multi-task graph attention framework. It also supports batch evaluation of molecular datasets and provides detailed explanation and optimal range of each property. ProTox-II [39] is a web server that provides predictions of oral toxicity for small molecules. It incorporates new features and functionalities, such as a more comprehensive database of toxic compounds, a more accurate prediction model, a more user-friendly interface, and a more informative output. ProTox-II aims to assist researchers and practitioners in the fields of drug discovery, chemical safety assessment, and environmental toxicology by providing reliable and interpretable toxicity predictions.

The prediction results showed that AF-399/40713777 has a desirable ADME profile, with proper logP and LogD, and acceptable LogS. AI-942/42301830 shows a relatively worse ADME profile, mainly because of the five phenyl rings in its structure. But both of them show desirable toxicity data. The predicted LD₅₀ of AF-399/40713777 and AI-942/42301830 is 1000 mg/kg and 2000 mg/kg, respectively.

3. Materials and Methods

3.1. Target Protein Structure and Ligand Preparation

The cocrystal complex structure of M^pro with 23R (PDB ID: 7KX5) obtained from a protein databank (rcsb.org) was used as the receptor structure. The protein preparation wizard, Schrödinger Suite 2022-2, was used to prepare the structure. Co-crystalized metals and ions and non-water solvents were deleted. The bond orders were reassigned using the Chemical Component Dictionary. The protein’s hydrogens were deleted and re-added. The terminal oxygen atoms were added, and water molecules beyond 8 Å of the binding ligand were deleted. The missing side chains of the whole structure were added and the hydrogen bonds were reassigned using PROPKA [40]. The whole structure was then energy minimized using the OPLS4 force field. The other parameters were set as default.

More than 3 million commercially available compounds from two vendors (Specs and Enamine) were used for screening. All ligands were prepared using the Ligprep module. The ionization state of all ligands was optimized with pH = 7.0 using Epik, and alternative stereoisomers were determined from the original 3D structures. The other parameters were set as default.

3.2. Molecular Docking-Based Virtual Screening

The binding pocket was identified using the receptor grid generation module. The cocrystal complex structure of M^pro with 23R (PDB ID: 7KX5) obtained from the protein databank (rcsb.org) was used to identify the binding pocket. All other parameters were set as default.

The virtual screening was carried out using the virtual screening workflow of Schrödinger Suite 2022-2. All prepared compound libraries were used for the screening. The binding site was previously generated using the receptor grid generation module of Schrödinger Suite 2022-2, which was identified using the position of 23R. The other parameters were set as default. With the binding grid box generated, Glide HTVS (high-throughput virtual screening) docking precision, standard precision (SP), and extra precision (XP) docking were used for screening in sequence. And the top 10% of best-scoring ligands for each procedure were selected for the next screening step. Eventually, the 10% best-scoring ligands (3208 ligands) after the extra precision (XP) docking were identified. The binding free energy (Molecular Mechanics/Generalized Born Surface Area, MM/GBSA) of the top 1000 ranking ligands with the receptor was estimated using the prime MMGBSA module [37]. The residues within 5 Å of the ligands were set as flexible, and all other parameters were set as default. The binding poses of these ligands were further checked manually. And the final candidates were selected based on the docking scores, binding free energy estimation, and manual selection (Table 1).

3.3. Initial Biological Evaluation of Selected Analogues

The enzyme inhibitory activity of selected candidates was evaluated using a previously reported well-developed fluorescence resonance energy transfer (FRET) assay [41].

A fluorescence resonance energy transfer (FRET) method was applied to measure M^pro inhibition of the tested compounds. Black 96-well plates, SARS-CoV-2 main protease solution (1 mg/mL, 500 mM of Tris, 150 mM of NaCl, 1 mM of EDTA, and 50% glycerol), and fluorescent substrate MCA-AVLQSGFR-Lys(Dnp)-Lys-NH₂ (20 mM in DMSO) were purchased from Beyotime Biotechnology (Shanghai, China). M^pro was diluted to 1.5 μM using an assay buffer (50 mM of Tris-HCl, 150 mM of NaCl, 20% glycerol, and pH = 7.3) and preserved at −20 °C. The fluorescent substrate was diluted to 500 μM in an assay buffer (50 mM of Tris-HCl, 150 mM of NaCl, 1 mM of EDTA, and pH = 7.3). Stock solutions (5 mM) of the test compounds were prepared in DMSO and diluted to the desired concentration using the assay buffer. Each well of a black 96-well plate was added 50 μL of assay buffer, 20 μL of diluted enzyme solution, and 20 μL of inhibitor solution. The control wells (no inhibitor) and blank wells (no enzyme) were measured in parallel experiments. Then, the 96-well plate was incubated at 37 °C in a shaking incubator for 10 min. After incubation, 10 µL of substrate solution was added per well to initiate the reaction. At the excitation wavelength of 320 nM and the emission wavelength of 405 nM, the fluorescence signal was measured every 10 s for 10 min using a SpectraMax iD5 multimode plate reader (Molecular Devices). Data from the first 60 s were linear fitted to calculate the slope values (V₀ for the control wells, and V_i for the test wells). %Inhibition (i%) in each well can, thus, be calculated using the following equation:

i% = 1 − V_i/V₀ × 100%

3.4. Molecular Dynamics Simulations for Interaction Analyses and Binding Free Energy Estimations

The binding poses of AF-399/40713777 and AI-942/42301830 with M^pro were investigated using molecular dynamics simulation. The initial docking poses of these two ligands were used to build the simulation system. The protein–ligand complex within the explicit solvent system with the OPLS4 force field was studied using the Desmond module of Schrödinger suite 2022-1. The atomic framework was solvated with a TIP3P water model with orthorhombic intermittent limit conditions for a 10 Å buffer region. The overlapping water molecules were eliminated, and Na⁺ or Cl⁻ was added as counter ions to neutralize the entire framework of atoms. An extra 0.15 M of NaCl was added into the system. The simulation was performed using an ensemble (NPT) of Nosé–Hoover thermostat and barostat to maintain a constant temperature of 300 K and a pressure of 1 bar in the system (Maestro, Schrödinger suits 2022-1). A hybrid energy minimization algorithm with 1000 steps of steepest descent, followed by conjugate gradient algorithms, was used. Then, 500 ns molecular dynamics simulations (Maestro, Schrödinger suits 2022-1) were performed and the post-dynamics simulation was analyzed using the simulation interaction diagram module. The binding free energy was calculated using the MM/GBSA module of Maestro based on the molecular dynamics simulations. The whole simulation time (500 ns) was sampled. The MM/GBSA calculations were estimated every 5 ns. The average of all estimations is presented as the results. The detailed data can be found in the Supplementary Materials.

3.5. ADMET Analysis

The ADMET profile of AF-399/40713777 and AI-942/42301830 were estimated using ADMET lab 2.0 (https://admetmesh.scbdd.com/ (accessed on 15 January 2023)). The molecular SIMLE strings were generated using Chemdraw and then submitted to the online system. An analysis report was then generated. The toxicity profile was evaluated using the ProTox-II server (https://tox-new.charite.de/ (accessed on 15 January 2023)).

4. Conclusions

In conclusion, we identified two novel non-peptidic inhibitors, AF-399/40713777 and AI-942/42301830, for SARS-CoV-2 M^pro using a molecular docking-based virtual screening approach. The initial enzyme inhibition assay indicated moderate inhibitory activities of both ligands toward M^pro. The molecular dynamics simulations investigated the possible binding poses of these two ligands. In addition, the calculated binding affinities (MM/GBSA) did translate well into experimentally determined inhibitory potencies. The predicted ADMET results showed that AF-399/40713777 possessed a desirable profile. The successfully identified hits suggest the feasibility of this hierarchical virtual screening approach in drug discovery. Moreover, these two hits provide novel scaffolds for exploring high-activity inhibitors targeting SARS-CoV-2 M^pro.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/molecules28145320/s1.

Author Contributions

M.G. conceptualized and carried out the work under the supervision of N.L. and Y.L.; D.K. and N.L. provided expert review throughout execution of the study and manuscript preparation. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the Academic Promotion Program of Shandong First Medical University (No. 2019LJ003).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Materials; further inquiries can be directed to the corresponding authors.

Acknowledgments

The authors acknowledge the technical support of Fabao Zhao and Lanlan Jing from Shandong University.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds used in the method part of the article are available from the authors.

References

Lee, T.C.; Murthy, S.; Del Corpo, O.; Senécal, J.; Butler-Laporte, G.; Sohani, Z.N.; Brophy, J.M.; McDonald, E.G. Remdesivir for the treatment of COVID-19: A systematic review and meta-analysis. Clin. Microbiol. Infect. 2022, 28, 1203–1210. [Google Scholar] [CrossRef] [PubMed]
Reina, J.; Iglesias, C. Nirmatrelvir plus ritonavir (Paxlovid) a potent SARS-CoV-2 3CLpro protease inhibitor combination. Rev. Esp. Quimioter. 2022, 35, 236–240. [Google Scholar] [CrossRef] [PubMed]
Jorgensen, S.C.J.; Kebriaei, R.; Dresser, L.D. Remdesivir: Review of Pharmacology, Pre-clinical Data, and Emerging Clinical Experience for COVID-19. Pharmacotherapy 2020, 40, 659–671. [Google Scholar] [CrossRef] [PubMed]
Reis, S.; Metzendorf, M.I.; Kuehn, R.; Popp, M.; Gagyor, I.; Kranke, P.; Meybohm, P.; Skoetz, N.; Weibel, S. Nirmatrelvir combined with ritonavir for preventing and treating COVID-19. Cochrane Database Syst. Rev. 2022, 9, Cd015395. [Google Scholar] [PubMed]
V’Kovski, P.; Kratzel, A.; Steiner, S.; Stalder, H.; Thiel, V. Coronavirus biology and replication: Implications for SARS-CoV-2. Nat. Rev. Microbiol. 2021, 19, 155–170. [Google Scholar] [CrossRef]
Yan, W.; Zheng, Y.; Zeng, X.; He, B.; Cheng, W. Structural biology of SARS-CoV-2: Open the door for novel therapies. Signal Transduct. Target. Ther. 2022, 7, 26. [Google Scholar] [CrossRef]
Scialo, F.; Daniele, A.; Amato, F.; Pastore, L.; Matera, M.G.; Cazzola, M.; Castaldo, G.; Bianco, A. ACE2: The Major Cell Entry Receptor for SARS-CoV-2. Lung 2020, 198, 867–877. [Google Scholar] [CrossRef]
Wrapp, D.; Wang, N.; Corbett, K.S.; Goldsmith, J.A.; Hsieh, C.L.; Abiona, O.; Graham, B.S.; McLellan, J.S. Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation. Science 2020, 367, 1260–1263. [Google Scholar] [CrossRef] [Green Version]
Jin, Z.; Du, X.; Xu, Y.; Deng, Y.; Liu, M.; Zhao, Y.; Zhang, B.; Li, X.; Zhang, L.; Peng, C.; et al. Structure of M(pro) from SARS-CoV-2 and discovery of its inhibitors. Nature 2020, 582, 289–293. [Google Scholar] [CrossRef] [Green Version]
Shin, D.; Mukherjee, R.; Grewe, D.; Bojkova, D.; Baek, K.; Bhattacharya, A.; Schulz, L.; Widera, M.; Mehdipour, A.R.; Tascher, G.; et al. Papain-like protease regulates SARS-CoV-2 viral spread and innate immunity. Nature 2020, 587, 657–662. [Google Scholar] [CrossRef]
Yin, W.; Mao, C.; Luan, X.; Shen, D.D.; Shen, Q.; Su, H.; Wang, X.; Zhou, F.; Zhao, W.; Gao, M.; et al. Structural basis for inhibition of the RNA-dependent RNA polymerase from SARS-CoV-2 by remdesivir. Science 2020, 368, 1499–1504. [Google Scholar] [CrossRef] [PubMed]
Wang, Y.; Anirudhan, V.; Du, R.; Cui, Q.; Rong, L. RNA-dependent RNA polymerase of SARS-CoV-2 as a therapeutic target. J. Med. Virol. 2021, 93, 300–310. [Google Scholar] [CrossRef] [PubMed]
Lin, H.X.J.; Cho, S.; Meyyur Aravamudan, V.; Sanda, H.Y.; Palraj, R.; Molton, J.S.; Venkatachalam, I. Remdesivir in Coronavirus Disease 2019 (COVID-19) treatment: A review of evidence. Infection 2021, 49, 401–410. [Google Scholar] [CrossRef] [PubMed]
Yan, L.; Yang, Y.; Li, M.; Zhang, Y.; Zheng, L.; Ge, J.; Huang, Y.C.; Liu, Z.; Wang, T.; Gao, S.; et al. Coupling of N7-methyltransferase and 3’-5’ exoribonuclease with SARS-CoV-2 polymerase reveals mechanisms for capping and proofreading. Cell 2021, 184, 3474–3485.e11. [Google Scholar] [CrossRef]
Yadav, R.; Chaudhary, J.K.; Jain, N.; Chaudhary, P.K.; Khanra, S.; Dhamija, P.; Sharma, A.; Kumar, A.; Handu, S. Role of Structural and Non-Structural Proteins and Therapeutic Targets of SARS-CoV-2 for COVID-19. Cells 2021, 10, 821. [Google Scholar] [CrossRef]
Cao, Y.; Yang, R.; Lee, I.; Zhang, W.; Sun, J.; Wang, W.; Meng, X. Characterization of the SARS-CoV-2 E Protein: Sequence, Structure, Viroporin, and Inhibitors. Protein Sci. 2021, 30, 1114–1130. [Google Scholar] [CrossRef]
Hu, Q.; Xiong, Y.; Zhu, G.H.; Zhang, Y.N.; Zhang, Y.W.; Huang, P.; Ge, G.B. The SARS-CoV-2 main protease (M(pro)): Structure, function, and emerging therapies for COVID-19. MedComm 2022, 3, e151. [Google Scholar] [CrossRef]
Ma, C.; Xia, Z.; Sacco, M.D.; Hu, Y.; Townsend, J.A.; Meng, X.; Choza, J.; Tan, H.; Jang, J.; Gongora, M.V.; et al. Discovery of Di- and Trihaloacetamides as Covalent SARS-CoV-2 Main Protease Inhibitors with High Target Specificity. J. Am. Chem. Soc. 2021, 143, 20697–20709. [Google Scholar] [CrossRef]
Chia, C.S.B.; Xu, W.; Shuyi Ng, P. A Patent Review on SARS Coronavirus Main Protease (3CL(pro)) Inhibitors. ChemMedChem 2022, 17, e202100576. [Google Scholar] [CrossRef]
Chen, J.; Ali, F.; Khan, I.; Zhu, Y.Z. Recent progress in the development of potential drugs against SARS-CoV-2. Curr. Res. Pharmacol. Drug Discov. 2021, 2, 100057. [Google Scholar] [CrossRef]
Macip, G.; Garcia-Segura, P.; Mestres-Truyol, J.; Saldivar-Espinoza, B.; Pujadas, G.; Garcia-Vallvé, S. A Review of the Current Landscape of SARS-CoV-2 Main Protease Inhibitors: Have We Hit the Bullseye Yet? Int. J. Mol. Sci. 2021, 23, 259. [Google Scholar] [CrossRef] [PubMed]
Li, J.; Lin, C.; Zhou, X.; Zhong, F.; Zeng, P.; Yang, Y.; Zhang, Y.; Yu, B.; Fan, X.; McCormick, P.J.; et al. Structural basis of main proteases of coronavirus bound to drug candidate PF-07321332. J. Virol. 2022, 96, e0201321. [Google Scholar] [CrossRef] [PubMed]
Vuong, W.; Fischer, C.; Khan, M.B.; van Belkum, M.J.; Lamer, T.; Willoughby, K.D.; Lu, J.; Arutyunova, E.; Joyce, M.A.; Saffran, H.A.; et al. Improved SARS-CoV-2 M(pro) inhibitors based on feline antiviral drug GC376: Structural enhancements, increased solubility, and micellar studies. Eur. J. Med. Chem. 2021, 222, 113584. [Google Scholar] [CrossRef] [PubMed]
Chen, W.; Shao, Y.; Peng, X.; Liang, B.; Xu, J.; Xing, D. Review of preclinical data of PF-07304814 and its active metabolite derivatives against SARS-CoV-2 infection. Front. Pharmacol. 2022, 13, 1035969. [Google Scholar] [CrossRef]
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. [Google Scholar] [CrossRef]
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. [Google Scholar] [CrossRef] [Green Version]
Zhang, C.H.; Stone, E.A.; Deshmukh, M.; Ippolito, J.A.; Ghahremanpour, M.M.; Tirado-Rives, J.; Spasov, K.A.; Zhang, S.; Takeo, Y.; Kudalkar, S.N.; et al. Potent Noncovalent Inhibitors of the Main Protease of SARS-CoV-2 from Molecular Sculpting of the Drug Perampanel Guided by Free Energy Perturbation Calculations. ACS Cent. Sci. 2021, 7, 467–475. [Google Scholar] [CrossRef]
Kitamura, N.; Sacco, M.D.; Ma, C.; Hu, Y.; Townsend, J.A.; Meng, X.; Zhang, F.; Zhang, X.; Ba, M.; Szeto, T.; et al. Expedited Approach toward the Rational Design of Noncovalent SARS-CoV-2 Main Protease Inhibitors. J. Med. Chem. 2022, 65, 2848–2865. [Google Scholar] [CrossRef]
Unoh, Y.; Uehara, S.; Nakahara, K.; Nobori, H.; Yamatsu, Y.; Yamamoto, S.; Maruyama, Y.; Taoda, Y.; Kasamatsu, K.; Suto, T.; et al. Discovery of S-217622, a Noncovalent Oral SARS-CoV-2 3CL Protease Inhibitor Clinical Candidate for Treating COVID-19. J. Med. Chem. 2022, 65, 6499–6512. [Google Scholar] [CrossRef]
Hou, N.; Shuai, L.; Zhang, L.; Xie, X.; Tang, K.; Zhu, Y.; Yu, Y.; Zhang, W.; Tan, Q.; Zhong, G.; et al. Development of Highly Potent Noncovalent Inhibitors of SARS-CoV-2 3CLpro. ACS Cent. Sci 2023, 9, 217–227. [Google Scholar] [CrossRef]
Gurung, A.B.; Ali, M.A.; Lee, J.; Farah, M.A.; Al-Anazi, K.M. An Updated Review of Computer-Aided Drug Design and Its Application to COVID-19. BioMed Res. Int. 2021, 2021, 8853056. [Google Scholar] [CrossRef] [PubMed]
Rahman, M.M.; Islam, M.R.; Akash, S.; Mim, S.A.; Rahaman, M.S.; Emran, T.B.; Akkol, E.K.; Sharma, R.; Alhumaydhi, F.A.; Sweilam, S.H.; et al. In silico investigation and potential therapeutic approaches of natural products for COVID-19: Computer-aided drug design perspective. Front. Cell. Infect. Microbiol. 2022, 12, 929430. [Google Scholar] [CrossRef] [PubMed]
Capdeville, R.; Buchdunger, E.; Zimmermann, J.; Matter, A. Glivec (STI571, imatinib), a rationally developed, targeted anticancer drug. Nat. Rev. Drug Discov. 2002, 1, 493–502. [Google Scholar] [CrossRef] [PubMed]
Shahwan, M.; Khan, M.S.; Husain, F.M.; Shamsi, A. Understanding binding between donepezil and human ferritin: Molecular docking and molecular dynamics simulation approach. J. Biomol. Struct. Dyn. 2022, 40, 3871–3879. [Google Scholar] [CrossRef] [PubMed]
Hirano, Y.; Okimoto, N.; Fujita, S.; Taiji, M. Molecular Dynamics Study of Conformational Changes of Tankyrase 2 Binding Subsites upon Ligand Binding. ACS Omega 2021, 6, 17609–17620. [Google Scholar] [CrossRef]
Fischer, A.; Smiesko, M.; Sellner, M.; Lill, M.A. Decision making in structure-based drug discovery: Visual inspection of docking results. J. Med. Chem. 2021, 64, 2489–2500. [Google Scholar] [CrossRef]
Ahinko, M.; Niinivehmas, S.; Jokinen, E.; Pentikäinen, O.T. Suitability of MMGBSA for the selection of correct ligand binding modes from docking results. Chem. Biol. Drug Des. 2019, 93, 522–538. [Google Scholar] [CrossRef] [Green Version]
Xiong, G.; Wu, Z.; Yi, J.; Fu, L.; Yang, Z.; Hsieh, C.; Yin, M.; Zeng, X.; Wu, C.; Lu, A.; et al. ADMETlab 2.0: An integrated online platform for accurate and comprehensive predictions of ADMET properties. Nucleic Acids Res. 2021, 49, W5–W14. [Google Scholar] [CrossRef]
Banerjee, P.; Eckert, A.O.; Schrey, A.K.; Preissner, R. ProTox-II: A webserver for the prediction of toxicity of chemicals. Nucleic Acids Res. 2018, 46, W257–W263. [Google Scholar] [CrossRef] [Green Version]
Rostkowski, M.; Olsson, M.H.; Søndergaard, C.R.; Jensen, J.H. Graphical analysis of pH-dependent properties of proteins predicted using PROPKA. BMC Struct. Biol. 2011, 11, 6. [Google Scholar] [CrossRef] [Green Version]
Gao, S.; Sylvester, K.; Song, L.; Claff, T.; Jing, L.; Woodson, M.; Weiße, R.H.; Cheng, Y.; Schäkel, L.; Petry, M.; et al. Discovery and Crystallographic Studies of Trisubstituted Piperazine Derivatives as Non-Covalent SARS-CoV-2 Main Protease Inhibitors with High Target Specificity and Low Toxicity. J. Med. Chem. 2022, 65, 13343–13364. [Google Scholar] [CrossRef] [PubMed]

Figure 1. Life cycle of SARS-CoV-2.

Figure 2. Crystal structure of M^pro homodimer and the catalytic site. Ligand: green stick/sphere; key residues: pink stick; hydrogen bond: yellow dashed line; π-π stacking: purple dashed line; PDB code: 7KX5 [18].

Figure 3. Representative peptidic covalent M^pro inhibitors.

Figure 4. Representative non-peptidic M^pro inhibitors.

Figure 5. The flow chart of the discovery of inhibitors against SARS-CoV-2 M^pro via docking-based virtual screening.

Figure 6. The RMSD and RMSF plots of AF-399/40713777 and AI-942/42301830 binding to M^pro generated via MD simulation.

Figure 7. The protein–ligand binding conformations generated from the molecular dynamics simulation trajectory. Ligand: green stick; key residues: pink stick; hydrogen bond: yellow dashed line; π-π stacking: purple dashed line. The structure files (PDB files) can be found in the Supplementary Materials.

Figure 8. The ADMET profile prediction of small molecules processed using ADMET lab 2.0.

Figure 9. The toxicity prediction of representative ligands processed using the ProTox-II server.

Table 1. Docking scores and MM/GBSA (kcal/mol) of selected ligands as anti-SARS-CoV-2 M^pro candidates.

ID	Docking Score	MM/GBSA
AN-329/15538195	−7.769	−102.03
AF-399/40713777	−7.929	−91.41
AN-655/14907067	−7.882	−86.30
AK-968/12101028	−7.911	−86.00
AG-690/13705944	−8.568	−85.86
AK-968/37129380	−7.829	−85.48
AG-205/36953218	−8.248	−85.06
AH-487/11927009	−7.764	−84.29
AI-942/42301830	−8.856	−87.32
AN-329/14726055	−8.798	−88.09
Z54217235	−8.670	−87.01
Z91218686	−9.095	−86.00
Z20007584	−10.038	−64.14
Z92376193	−9.139	−71.91
Z929753284	−9.018	−49.50
Z245966642	−8.850	−61.06
Z1603682175	−8.827	−51.59

Table 2. Enzyme inhibitory activities of selected candidates.

ID	% Inhibition @ 100 μM	% Inhibition @ 100 μM Average of 3 Samples
AN-329/15538195	18.9%	11.80% ± 5.65
AF-399/40713777	23.7%	47.11% ± 0.84
AN-655/14907067	13.3%	2.71% ± 7.02
AK-968/12101028	4.3%
AG-690/13705944	12.4%
AK-968/37129380	12.5%
AG-205/36953218	−13.8%
AH-487/11927009	−0.1%
AI-942/42301830	17.8%	37.67% ± 2.61
AN-329/14726055	−5.1%	−0.43% ± 3.97
Z54217235	−14.8%
Z91218686	−6.1%
Z20007584	4.5%
Z92376193	−5.0%
Z929753284	18.4%
Z245966642	3.8%
Z1603682175	7.2%

Table 3. MM/GBSA (kcal/mol) of AF-399/40713777 and AI-942/42301830 with SARS-CoV-2 M^{pro 1}.

Ligands	ΔG_Bind	ΔG_Coulomb	ΔG_Hbond	ΔG_Lipo	ΔG_vdW
AF-399/40713777	−66.07	−15.60	−0.66	−16.14	−51.01
AI-942/42301830	−78.32	−11.98	−0.51	−25.66	−62.19

¹ The whole simulations trajectory was sampled, and the estimation was calculated every 5 ns. The averages of all calculations are presented as the results. Complete data and standard deviations can be found in the Supplementary Materials.

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Gao, M.; Kang, D.; Liu, N.; Liu, Y. In Silico Discovery of Small-Molecule Inhibitors Targeting SARS-CoV-2 Main Protease. Molecules 2023, 28, 5320. https://doi.org/10.3390/molecules28145320

AMA Style

Gao M, Kang D, Liu N, Liu Y. In Silico Discovery of Small-Molecule Inhibitors Targeting SARS-CoV-2 Main Protease. Molecules. 2023; 28(14):5320. https://doi.org/10.3390/molecules28145320

Chicago/Turabian Style

Gao, Menghan, Dongwei Kang, Na Liu, and Yanna Liu. 2023. "In Silico Discovery of Small-Molecule Inhibitors Targeting SARS-CoV-2 Main Protease" Molecules 28, no. 14: 5320. https://doi.org/10.3390/molecules28145320

APA Style

Gao, M., Kang, D., Liu, N., & Liu, Y. (2023). In Silico Discovery of Small-Molecule Inhibitors Targeting SARS-CoV-2 Main Protease. Molecules, 28(14), 5320. https://doi.org/10.3390/molecules28145320

Article Menu