Critical assessment of popular biomolecular force fields for molecular dynamics simulations of folding and enzymatic activity of main protease of coronavirus SARS-CoV-2.

The main cysteine protease (Mpro) of coronavirus SARS-CoV-2 has become a promising target for computational development in anti-COVID-19 treatments. Here, we benchmarked the performance of six biomolecular molecular dynamics (MD) force fields (OPLS-AA, CHARMM27, CHARMM36, AMBER03, AMBER14SB and GROMOS G54A7) and three water models (TIP3P, TIP4P and SPC) for reproducing the native fold and the enzymatic activity of Mpro as monomeric and dimeric units. The MD sampling up to 1 µs suggested that the proper choice of the force fields and water models plays an essential role in reproducing the tertiary structure and the inter-residue distance between the catalytic dyad His41-Cys145. We found that while most benchmarked all-atom force fields reproduce well the native fold of Mpro, the CHARMM27/TIP3P and OPLS-AA/TIP4P setups revealed a good performance in reproducing the structure of the catalytic domain. In addition, these FF setups were also well-adopted for MD sampling of Mpro at the physiologic conditions by mimicking the presence of 100 mM NaCl and the elevated temperature of 310 K. Finally, both FFs were also performed well in reproducing the native fold of Mpro in a dimeric form. Therefore, comparing the preservation of the native fold of Mpro and the stability of its catalytic site architecture, our MD benchmarking suggests that the OPLS-AA/TIP4P and CHARMM27/TIP3P MD setups at the physiologic conditions may be well-suited for rapid in silico screening and developing broad-spectrum anti-coronaviral therapeutic agents.

Assessment of nucleotide/nucleoside analog intervention in primer-dependent viral RNA-dependent RNA polymerases.

Nucleotide/nucleoside analogs (NAs) are important compounds used in antiviral drug development. To understand the action mode of NA drugs, we present an enzymology protocol to initially evaluate the intervention mechanism of the NTP forms of NAs on a coronaviral RNA-dependent RNA polymerase (RdRP). We describe the preparation of SARS-CoV-2 RdRP proteins and RNA constructs, followed by a primer-dependent RdRP assay to assess NTP forms of NAs. Two representative NA drugs, sofosbuvir and remdesivir, are used for demonstration of this protocol. For complete details on the use and execution of this protocol, please refer to Wu et al. (2021).

A COVID moonshot: assessment of ligand binding to the SARS-CoV-2 main protease by saturation transfer difference NMR spectroscopy.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the etiological cause of the coronavirus disease 2019, for which no effective antiviral therapeutics are available. The SARS-CoV-2 main protease (Mpro) is essential for viral replication and constitutes a promising therapeutic target. Many efforts aimed at deriving effective Mpro inhibitors are currently underway, including an international open-science discovery project, codenamed COVID Moonshot. As part of COVID Moonshot, we used saturation transfer difference nuclear magnetic resonance (STD-NMR) spectroscopy to assess the binding of putative Mpro ligands to the viral protease, including molecules identified by crystallographic fragment screening and novel compounds designed as Mpro inhibitors. In this manner, we aimed to complement enzymatic activity assays of Mpro performed by other groups with information on ligand affinity. We have made the Mpro STD-NMR data publicly available. Here, we provide detailed information on the NMR protocols used and challenges faced, thereby placing these data into context. Our goal is to assist the interpretation of Mpro STD-NMR data, thereby accelerating ongoing drug design efforts.

SARS-CoV-2 Papain-Like Protease Potential Inhibitors-In Silico Quantitative Assessment.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) encodes the papain-like protease (PLpro). The protein not only plays an essential role in viral replication but also cleaves ubiquitin and ubiquitin-like interferon-stimulated gene 15 protein (ISG15) from host proteins, making it an important target for developing new antiviral drugs. In this study, we searched for novel, noncovalent potential PLpro inhibitors by employing a multistep in silico screening of a 15 million compound library. The selectivity of the best-scored compounds was evaluated by checking their binding affinity to the human ubiquitin carboxy-terminal hydrolase L1 (UCH-L1), which, as a deubiquitylating enzyme, exhibits structural and functional similarities to the PLpro. As a result, we identified 387 potential, selective PLpro inhibitors, from which we retrieved the 20 best compounds according to their IC50 values toward PLpro estimated by a multiple linear regression model. The selected candidates display potential activity against the protein with IC50 values in the nanomolar range from approximately 159 to 505 nM and mostly adopt a similar binding mode to the known, noncovalent SARS-CoV-2 PLpro inhibitors. We further propose the six most promising compounds for future in vitro evaluation. The results for the top potential PLpro inhibitors are deposited in the database prepared to facilitate research on anti-SARS-CoV-2 drugs.

Computational Selectivity Assessment of Protease Inhibitors against SARS-CoV-2.

The pandemic of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) poses a serious global health threat. Since no specific therapeutics are available, researchers around the world screened compounds to inhibit various molecular targets of SARS-CoV-2 including its main protease (Mpro) essential for viral replication. Due to the high urgency of these discovery efforts, off-target binding, which is one of the major reasons for drug-induced toxicity and safety-related drug attrition, was neglected. Here, we used molecular docking, toxicity profiling, and multiple molecular dynamics (MD) protocols to assess the selectivity of 33 reported non-covalent inhibitors of SARS-CoV-2 Mpro against eight proteases and 16 anti-targets. The panel of proteases included SARS-CoV Mpro, cathepsin G, caspase-3, ubiquitin carboxy-terminal hydrolase L1 (UCHL1), thrombin, factor Xa, chymase, and prostasin. Several of the assessed compounds presented considerable off-target binding towards the panel of proteases, as well as the selected anti-targets. Our results further suggest a high risk of off-target binding to chymase and cathepsin G. Thus, in future discovery projects, experimental selectivity assessment should be directed toward these proteases. A systematic selectivity assessment of SARS-CoV-2 Mpro inhibitors, as we report it, was not previously conducted.

Assessment of effective imidazole derivatives against SARS-CoV-2 main protease through computational approach.

Because of the fast increase in deaths due to Corona Viral Infection in majority region in the world, the detection of drugs potent of this infection is a major need. With this idea, docking study was executed on eighteen imidazole derivatives based on 7-chloro-4-aminoquinoline against novel Coronavirus (SARS-CoV-2). In this study, we carried out a docking study of these molecules in the active site of SARS-CoV-2 main protease. The result indicate that Molecules N° 3, 7 and 14 have more binding energy with SARS-CoV-2 main protease recently crystallized (pdb code 6LU7) in comparison with the other imidazole derivatives and the two drug; Chloroquine and hydroxychloroquine. Because of the best energy of interaction, these three molecules could have the most potential antiviral treatment of COVID-19 than the other studied compounds. The structures with best affinity in the binding site of the protease have more than 3 cycles and electronegative atoms in the structure. This may increase the binding affinity of these molecules because of formation of π-bonds, halogen interactions and/or Hydrogen bond interactions between compounds and the enzyme. So, compounds with more cycles and electronegative atoms could have a potent inhibition of SARS-CoV-2 main protease.