key: cord-315604-a6fvsd45
authors: Maurya, Santosh K.; Maurya, Akhilesh Kumar; Mishra, Nidhi; Siddique, Hifzur R.
title: Virtual screening, ADME/T, and binding free energy analysis of anti-viral, anti-protease, and anti-infectious compounds against NSP10/NSP16 methyltransferase and main protease of SARS CoV-2
date: 2020-06-01
journal: J Recept Signal Transduct Res
DOI: 10.1080/10799893.2020.1772298
sha: 
doc_id: 315604
cord_uid: a6fvsd45

Recently, a pathogen has been identified as a novel coronavirus (SARS-CoV-2) and found to trigger novel pneumonia (COVID-19) in human beings and some other mammals. The uncontrolled release of cytokines is seen from the primary stages of symptoms to last acute respiratory distress syndrome (ARDS). Thus, it is necessary to find out safe and effective drugs against this deadly coronavirus as soon as possible. Here, we downloaded the three-dimensional model of NSP10/NSP16 methyltransferase (PDB-ID: 6w6l) and main protease (PDB-ID: 6lu7) of COVID-19. Using these molecular models, we performed virtual screening with our anti-viral, inti-infectious, and anti-protease compounds, which are attractive therapeutics to prevent infection of the COVID-19. We found that top screened compound binds with protein molecules with good dock score with the help of hydrophobic interactions and hydrogen bonding. We observed that protease complexed with Cyclocytidine hydrochloride (anti-viral and anti-cancer), Trifluridine (anti-viral), Adonitol, and Meropenem (anti-bacterial), and Penciclovir (anti-viral) bound with a good docking score ranging from −6.8 to −5.1 (Kcal/mol). Further, NSP10/NSP16 methyltransferase complexed with Telbivudine, Oxytetracycline dihydrate (anti-viral), Methylgallate (anti-malarial), 2-deoxyglucose and Daphnetin (anti-cancer) from the docking score of −7.0 to −5.7 (Kcal/mol). In conclusion, the selected compounds may be used as a novel therapeutic agent to combat this deadly pandemic disease, SARS-CoV-2 infection, but needs further experimental research. HIGHLIGHTS: NSP10/NSP16 methyltransferase and main protease complex of SARS CoV-2 bind with selected drugs. NSP10/NSP16 methyltransferase and protease interacted with drugs by hydrophobic interactions. Compounds show good DG binging free energy with protein complexes. Ligands were found to follow the Lipinski rule of five.

A new coronavirus, SARS-CoV-2, was found to cause pulmonary disease in the city of Wuhan, China, in December 2019 [1] . The name SARS-CoV-2 is due to its genomic RNA, which is about 82% similar to the Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV) [2] . These viruses are found to link with clade b of the Betacoronavirus genus [3] . The disease caused by SARS-CoV-2 is known as coronavirus disease-19 . At the beginning of the COVID-19 outbreak/spread, most of the cases were found to connect to the wild animal, seafood, and animal market in Wuhan [4] . These human foods have been found efficient to enhance the number of transmissions among human and later infected humans promotes human to human transmission, which further leads to the exponential growth of these viruses [5] . The World Health Organization (WHO) declared this outbreak as pandemic on 11 March 2020. As of May 16, there are >4.59 million cases globally, with a $3.7% casefatality rate (https://www.worldometers.info/coronavirus).

Coronaviruses (CoVs) have been found to affect the respiratory tract and lead light to harsh respiratory tract infections in humans [6] . In the past two decades, the SARS-CoV (severe acute respiratory syndrome-coronavirus) and the MERS-CoV (middle east respiratory syndrome-coronavirus) have been found to emerge from animal reservoirs and lead to global epidemics with great anguish and fatality [7] .

The coronaviruses have been found to encode two proteases (translated nonstructural proteins), i.e. a papain-like protease (Plpro) and main proteases. The main protease is also known as 3-C-like protease (M-pro) [8] . It has been shown that the main Protease (Mpro, also called 3CLpro) is one of the best-characterized drug targets among coronaviruses [9] . Protease is necessary for the processing of polyproteins, which are translated from the SARS CoV-2 RNA [8] . The main Protease of SARS CoV-2 operates at about 11 sites of cleavage on the large polyprotein 1ab (replicase 1ab, size ¼ $790 kDa) and the recognition sequence is found to be Leu-Gln#Ser-Ala-Gly) (where # marks the site of cleavage). It has been seen that inhibiting the activity of the main protease enzyme could be an effective measure to check viral replication [10] .

Various forms of the viruses are found to perform replication in the cytoplasm of eukaryotes through evolving 2 0 -Omethyltransferases (2 0 -O-MTase) for modification of their mRNAs and carry a Cap-1 structure such as, m7GpppNm at the 5 0 end, which facilitate viral duplication and help to escape recognition against innate immune responses done by some immune cells in the hosts [11] . The stimulatory factor NSP10, along with NSP16, stabilizes the pocket of SAM-binding and continuing the RNA-binding groove of the substrate of NSP16, which is observed by crystallographic and biochemical analysis [12] . It has been observed that interference in the interaction between NSP16 and NSP10 of SARS-CoV by short peptides can inhibit the activity of 2 0 -O-MTase [13] . 2 0 -O-MTase was also found to catalyze the transfer of methyl group SAM or AdoMet to RNA substrate and produce S-adenosyl-L-homocysteine as a byproduct [14] . It was evidenced that the stimulation of NSP16 methyltransferase activity by NSP10 is a common and primary mechanism for the coronaviruses. It has been further demonstrated that NSP10 is synonymous in the stimulatory function among distinct coronaviruses, and the peptide, which is derived from the conserved domain of MHV NSP10, shows an inhibitory effect on 2 0 -O-methyltransferase activity during virus duplication. These results might use to find out specific anti-COVID-19 drugs to control infection [11] .

It was proved that the availability of FDA-approved anti-RNA-dependent RNA Polymerase (anti-RdRp) drugs could help in the treatment of patients and reduce the chances of the danger of the viral infection caused by SARS CoV2. These drugs can tightly bind to the RdRp and thus may be useful to treat the disease caused by the SARS CoV2. Furthermore, it has also been shown that Ribavirin, vitamin B12, and nicotinamide may be helpful against COVID 19 infections [12] . Moreover, it has also been shown that the anti-virals Velpatasvir and Ledipasvir are particularly attractive as therapeutics to combat the SARS CoV-2 with minimal side effects (fatigue and headache). The drugs Harvoni (Ledipasvir/Sofosbuvir) and Epclusa (Velpatasvir/Sofosbuvir) with their dual inhibitory actions on two viral enzymes could be very effective [15] .

The present study provided a comprehensive target of the resolved COVID19 structure of Mpro and NSP10/NSP16 Methyltransferase and found a suitable approach against this disease. We have performed in silico study on anti-viral, antiinfectious, and anti-protease compounds against methyltransferase-stimulatory factor complex and main protease.

The in-silico studies were performed on KBS Desktop having 12 GB RAM, Intel i5 generation with 4 cores in Indian Institute of Information Technology, Allahabad, UP, India. Grid-based Ligand Docking with Energetics (GLIDE) module of maestro 12.0 (Schrodinger LLC 2019, USA) was used for our study.

Anti-viral, anti-infectious, and anti-protease compounds were downloaded from the online database Selleckchem (https:// www.selleckchem.com/), in SDF format. LigPrep module of Maestro12.0 (Schrodinger) was used to prepare the ligand library for the Docking purpose. Ionization of ligands was retained in actual states with the realistic length of bonds as well as bond angles, ring conformation and tautomers were generated, using the OPLS-2005 force field.

The high-resolution structure of SARS-CoV2 main Protease (PDB ID: 6lu7) and NSP10/NSP16 Methyltransferase (PDB ID: 6w61) were downloaded online from the protein data bank (https://www.rcsb.org/). The protein preparation wizard of Maestro 12.0 was used for the preparation of the protein structure. The selected structure was processed for creating disulfide bonds, assigning proper bond orders, and the addition of missing hydrogens in the raw structure. Hydrogen bonds optimization was assigned using non-hydrogen atoms of protein, and the structure was energy minimized until the RMSD (root mean square deviation) reaches the value of 0.3 Å.

Protease and methyltransferase were subjected to the Sitemap module of Maestro 12.0 for the prediction of active binding sites. A total of 4 binding sites was generated, and the top one was selected for receptor grid generation. The grid was generated using the default setting of the receptor grid generation module of maestro 12.0, i.e. 0.25 was partial atomic charge and were 1.0 Å was Van der Waals radii of receptor atoms. The center of the sitemap position was selected to generate a grid box by using the receptor grid generation module of Maestro.

Virtual screening of prepared library was operated by using virtual screening workflow (VSW) of GLIDE (maestro 12.0, Schr€ odinger, LLC, NY, USA). VSW has many tools to calculate the properties of ligands, and filter ligands such as QuikProp and prefilters by Lipinski rule were selected. VWS of Maestro has three type docking precision: (1) . HTVS (High Throughput Virtual Screening) for fast screening of huge number ligand's docking, (2) standard-precision (SP) for screening and docking of ligands of unknown quality in large numbers, and (3) extra-precision (XP) is a powerful and discriminating procedure for docking ligands. Extra precision (XP) docking was operated on prepared ligands by keeping default parameters, i.e. scaling factor at 0.80 and partial charge at 0.15. The binding affinity of docked ligands for the receptor protease at the active site was generated from the docking binding energy.

Binding free energy of protease enzyme and docked ligands were calculated using the Prime MM-GBSA (Molecular Mechanics-Generalized Born Surface Area) modules of Maestro 12.0, which run on the OPLS_2005 molecular mechanic's energies, nonpolar solvation and VSGB solvent model [16] . Prime MM-GBSA module was operated using pose viewer file of ligand and receptor (generated after docking) to generate binding free energy. The MM-GBSA analysis employs molecular mechanics and solvent accessibility methods to depict free energies. This analysis is efficient and quite reliable for studying the energetics of biomolecular systems. Prediction of the binding free energies in biomolecular such as, protein, protein-drug complexes is the application of this approach. The following descriptors were used to generate changes in energy upon binding.

MM-GBSA DG bind ¼ Ligand binding free energy, MM-GBSA E complex ¼ complex free energy, MM-GBSA E protein ¼ free energy of the receptor without the ligand, MM-GBSA E Ligand ¼ unbound ligand-free energy.

Absorption, distribution, metabolism, excretion/toxicity (ADME/T) properties are used to eliminate inappropriate compounds before invested inappropriate compounds, time, and money. QuikProp module in VSW (Maestro 12.0) calculated the ADME information of ligands. QikProp module calculates various properties and provides values for descriptors, which predict drug-likeness of ligands and generated functional groups, rings, and a number of elements used for the prediction of descriptors of the ligand/s.

The crystal structure of the NSP10/NSP16 Methyltransferase (PDB ID: 6w6l) and main Protease (PDB ID: 6lu7) of SARS-CoV-2 was used for docking purposes. Molecular docking has been performed using the extra precision (XP) mode of gridbased ligand docking with energetics (GLIDE). We used known anti-viral, anti-infectious, and anti-protease inhibitors for the docking purpose (Figures 1 and 2) . We selected the top 5 compounds with good docking scores, which were docked with proteases and NSP10/NSP16 methyltransferase. Our result highlighted that; Cyclocytidine, Hydrochloride, Trifluridine, Adonitol Meropenem, and Penciclovir yielded a preeminent dock score with the main Protease of SARS-COV-2 with docking score of À6.8, À6.0, À5.5, À5.2, and À5.1 Kcal/mol respectively (Table 1) Protease and Telbivudine, Oxytetracycline dihydrate (antiviral), Methylgallate (anti-malarial), 2-deoxyglucose (anti-cancer), and Daphnetin (anti-inflammatory, and antitumor) showed good docking score with methyltransferase-stimulatory factor complex of NSP16 and NSP10 of SARS-CoV-2.

Protein-ligands interactions are generally highlighted by hydrophobic interactions and hydrogen bonds, which play a significant role in predicting the binding affinity of the ligand with proteins. All the selected compounds in the dataset were docked into the active site of the NSP10/NSP16 methyltransferase and the main Protease of SAR-CoV-2 (Figures 3  and 4 ) depicts the binding conformations of the selected compounds in the binding pocket of the NSP10/NSP16 methyltransferase and protease. The active site of the NSP10/NSP16 methyltransferase showed hydrophobic interaction with the help of Met248, Pro252, Phe251, Cys47, and Pro60 amino acids. The amino acids are present in the active site of Protease like Phe112, Phe8, and Phe294 are found to create hydrophobic interaction with top ligands. Besides Hydrophobic interactions, there are some hydrogen bonds present between amino acids and compounds. In the case of NSP10/NSP16 methyltransferase, mostly Gln88, Glu61, and thr48 formed hydrogen bonds with the atoms of ligands and in the active site of main protease complex Gln110, Thr111, Gln127, and Asn151 amino acids have been found to form a hydrogen bond with the atoms of ligands (Table 2) .

A good number of amino acids are found to appear in the hydrophobic interactions and hydrogen bonding. Moreover, the active site of the NSP10/NSP16 methyltransferase and main protease have mostly hydrophobic interactions with the ligands. Hydrophobic interaction successfully delineates specific functional groups that may be responsible for the hydrophobic generating effect of these compounds with strong binding affinity against target proteins and may play a highly influencing against SARS-CoV-2 infection

MM/GBSA calculation was performed to estimate the free binding energy of ligands with docked proteins. Docked top five compounds from ligands library with proteins showed the binding-free energy values. Compounds selected in virtual screening has been found to show the good result of binding free energy with their respective Proteins (Table 3) . All the compounds showed good binding free energy with their respective proteins of SARS-CoV-2.

ADME/T properties of the lead compounds ( Figure 4) were appraised by using the Qikprop application of Maestro 12.0). The most attractive aspect of compounds is their admirable QPlogo/w, QPlogS [17] molecular weight, H-bond donor and acceptor, QplogBB, and percentage of human oral value which satisfy the Lipinski rule of five [18] (Table 4 ). Moreover, polar surface area, high oral bioavailability, H-bond donors, and acceptors are being important criteria for the development of therapeutic agents. All these models unscrew the qualitative prediction and ranking of absorption, formulation effects on drug permeability, determining the mechanism(s) of permeability, and the potential for transporter-mediated drug-drug interactions. All the good scoring ligand has druglikeness properties by Lipinski's rule. It may also be predictable that the combination of these drugs/s may be a good idea to prevent COVID-19 infection.

Drug discovery has become an important area in the field of COVID 19 spreading worldwide. So, the identification of drug targets is a crucial phase to modulate the expression of any targets. Many scientists are involved in searching for a cure for COVID 19. Most of the researchers are targeting protease complex and NSP10/NSP16 methyltransferase, which play a crucial role in inhibits the infection of the virus of COVID 19 [10, 11] . It has been seen that various anti-viral, anti-malarial drugs target so many proteins of SARS CoV-2 with a good binding affinity [12, 19] . Our selected compounds are already approved by the regulatory agencies, present in the market, and easy to purchase. It has been observed in this study that protease complex with Cyclocytidine hydrochloride (anti-viral and anti-cancer), Trifluridine (anti-viral), Adonitol (anti-bacterial), Meropenem (antibiotic and anti-bacterial), and Penciclovir (anti-viral) bound with a good docking score ranging from À6.8 to À5.1 (Kcal/mol). Moreover, the NSP10/NSP16 methyltransferase complex was found to docked with Telbivudine, Oxytetracycline dihydrate (anti-viral), Methylgallate (anti-malarial), 2-deoxyglucose (anti-cancer), and Daphnetin (anti-inflammatory, and antitumor) from the docking score of À7.0 to À5.7 (Kcal/mol). Hydrophobic interactions and hydrogen bonding have been observed between the protein complexes and our selected compounds. These hydrophobic interactions and hydrogen bonding have been found to show a strong bonding between the atoms of compounds and protein molecules [20] . This ligand-protein interaction also shows good DG free binding energy and the strength of bonding between the protein and ligands. All the chosen compound also follows the criteria of Lipinski's rule of five and show properties to have the physicochemical and pharmacokinetic parameters of drugability [18] . Our ligands probably inhibit protease and NSP10/NSP16 methyltransferase molecules involved in the COVID-19. Moreover, it may be possible that applying inhibitors in combination will exert a synergetic effect against the SARS CoV-2.

In this study, it has been shown that selected screened compounds have good docking score, high DG binding free energy as well as strong hydrophobic interactions, and follows the Lipinski rule of five. Thus, these compounds might have the potential to be utilized against the NSP10/NSP16 MW: Molecular weight; logs: Predicted aqueous solubility; logo/w: Predicted octanol/water partition coefficient; AccptH: Estimated number of hydrogen bonds that would be accepted by the solute from water molecules in an aqueous solution; DonorH: Estimated number of hydrogen bonds that would be donated by the solute to water molecules in an aqueous solution; QPPCaco: Predicted apparent Caco-2 cell permeability in nm/sec (Caco2 cells are a model for the gut-blood barrier); QPlogBB: Predicted brain/blood partition coefficient; % Human Oral Abs: Predicted human oral absorption on 0 to 100% scale; QPlog HERG: Predicted IC50 value for blockage of HERG K þ channels; SASA: Solvent Accessible Surface Area.

A pneumonia outbreak associated with a new coronavirus of probable bat origin

The species Severe acute respiratory syndrome-related coronavirus: classifying 2019-nCoV and naming it SARS-CoV-2

Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study

COVID-19: a novel zoonotic disease caused by a coronavirus from China: what we know and what we don't

A new coronavirus associated with human respiratory disease in China

Coronaviruses -drug discovery and therapeutic options

Coronavirus infections-more than just the common cold

From SARS to MERS: crystallographic studies on coronaviral proteases enable antiviral drug design

Coronavirus main proteinase (3CLpro) structure: basis for design of anti-SARS drugs

a-Ketoamides as broad-spectrum inhibitors of coronavirus and enterovirus replication: Structurebased design, synthesis, and activity assessment

Coronavirus NSP10/NSP16 methyltransferase can be targeted by nsp10-derived peptide in vitro and in vivo to reduce replication and pathogenesis

Structural; physicochemical; biochemical; pharmacokinetics and toxicity properties of compound

Prediction of the SARS-CoV-2 (2019-nCoV) 3C-like protease (3CL pro) structure: virtual screening reveals velpatasvir, ledipasvir, and other drug repurposing candidates

Short peptides derived from the interaction domain of SARS coronavirus nonstructural protein NSP10 can suppress the 2'-O-methyltransferase activity of nsp10/nsp16 complex

RIG-I-mediated antiviral responses to single-stranded RNA bearing 5'-phosphates

Virtual screening and repurposing of FDA approved drugs against COVID-19 main protease

The VSGB 2.0 model: a next generation energy model for high resolution protein structure modeling

Prediction of drug solubility from structure

Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings

Hydroxychloroquine, a less toxic derivative of chloroquine, is effective in inhibiting SARS-CoV-2 infection in vitro

Novel potent inhibitors of Plasmodium vivax dihydrofolate reductase: an in silico antimalarial drug discovery. IJPER

The authors acknowledge the Indian Institute of Information technology and Aligarh Muslim University for providing infrastructure facilities.

The authors declare no conflict of interest.