key: cord-0709064-xpnlq1t2 authors: Palanisamy, Kandhan; Rubavathy, S. M. Esther; Prakash, Muthuramalingam; Thilagavathi, Ramasamy; Hosseini-Zare, Maryam S.; Selvam, Chelliah title: Antiviral activities of natural compounds and ionic liquids to inhibit the Mpro of SARS-CoV-2: a computational approach date: 2022-01-28 journal: RSC advances DOI: 10.1039/d1ra08604a sha: 88571192d71b9f7993b15757b7171a92c797ead0 doc_id: 709064 cord_uid: xpnlq1t2 The recalcitrant spread of the COVID-19 pandemic produced by the novel coronavirus SARS-CoV-2 is one of the most destructive occurrences in history. Despite the availability of several effective vaccinations and their widespread use, this line of immunization often faces questions about its long-term efficacy. Since coronaviruses rapidly change, and multiple SARS-CoV-2 variants have emerged around the world. Therefore, finding a new target-based medication became a priority to prevent and control COVID-19 infections. The main protease (Mpro) is a salient enzyme in coronaviruses that plays a vital role in viral replication, making it a fascinating therapeutic target for SARS-CoV-2. We screened 0.2 million natural products against the Mpro of SARS-CoV-2 using the Universal Natural Product Database (UNPD). As well, we studied the role of ionic liquids (ILs) on the structural stabilization of Mpro. Cholinium-based ILs are biocompatible and used for a variety of biomedical applications. Molecular docking was employed for the initial screening of natural products and ILs against Mpro. To predict the drug-likeness features of lead compounds, we calculated the ADMET properties. We performed MD simulations for the selected complexes based on the docking outcomes. Using MM/PBSA approaches, we conclude that compounds NP-Hit2 (−25.6 kcal mol(−1)) and NP-Hit3 (−25.3 kcal mol(−1)) show stronger binding affinity with Mpro. The hotspot residues of Thr25, Leu27, His41, Met49, Cys145, Met165, and Gln189 strongly interacted with the natural compounds. Furthermore, naproxenate, ketoprofenate, and geranate, cholinium-based ILs strongly interact with Mpro and these ILs have antimicrobial properties. Our findings will aid in the development of effective Mpro inhibitors. The emergence of the highly transmissible disease SARS-CoV-2 has given rise to the pandemic COVID-19, causing widespread dread and a severe impact on the world economy. 1 An unexpected outbreak of COVID-19, triggered by the highly contagious SARS-CoV-2 virus, has affected millions of people worldwide. A large number of deaths are reported every day at the hands of this devastating disease. It is the seventh coronavirus to afflict humans, and because of its transmissible nature, it has a greater fatality rate. [2] [3] [4] The World Health Organization (WHO) proclaimed the COVID-19 outbreak a Public Health Emergency of International Concern (PHEIC) in March 2020. 5 Elderly persons with diabetes, pulmonary issues, and lung injury who require hospitalization or succumb to COVID-19 can suffer severe sickness, dry cough, and respiratory-related diseases. 6 According to a recent study, viral infection leads to long-term detrimental consequences such as kidney impairment and myocardial ischemia. Currently, seven COVID-19 vaccines have received approval for the Emergency Use Authorization (EUA) by WHO. The global immunization campaign is underway, but the infection rate continues to rise due to the emergence of COVID-19 mutations. 7 The most promising treatment option for patients with COVID-19 is convalescent plasma therapy. Nonetheless, numerous COVID-19 variations have emerged due to mutation, posing a new challenge to humankind. 8 Alternative immunization strategies are urgently required. COVID-19 strains are more infectious and have a higher fatality rate. As a result, new effective antiviral medications are essential to control the rising number of COVID-19 patients. Previous studies have shown that people with low immunity and pre-existing health problems like cancer, diabetes, hypertension, and respiratory disorders, among other things, are more vulnerable to severe infection. 6 The coronavirus latches onto the human spike protein and interacts with the human angiotensin-converting enzyme 2 (hACE2). SARS-CoV-2 is a beta coronavirus belonging to the Coronaviridae family, single-strand positive RNA virus ((+) ssRNA), capable of infecting animals to humans (i.e., zoonotic). 9 It contains two cysteine proteases and four structural proteins. The structural proteins of spike glycoprotein (S protein) are involved in virus generation and replication, as well as membrane (M), envelope (E), and nucleocapsid (N) proteins. 10 The acid-base forms of quercetin and its metabolites potential inhibitors against the targets of furin and S protein of SARS-CoV-2. 11 Furthermore, the replication processes are mediated by the two proteases 3C-like protease (3CLpro) or main protease (Mpro) and Papain-like protease (PLpro). 12 The viral polyproteins are cleaved into 11 nonstructural proteins (NPs) by Mpro, while three nonstructural proteins are cleaved by PLpro. Other proteolytic enzymes, which are more similar to human proteases, are more complicated to target. SARS CoV-2 Mpro is a promising drug target for developing inhibitors for COVID-19 infections. 13 15 Mutations occurred in 282 out of 306 residues in the Mpro. Remaining 24 residues had no reported mutations, these nonmutational residues are called as cold spot residues. 16 According to recent research, vaccination options for COVID-19 mutations are limited. To counteract and reduce COVID-19 infection rates, we need to focus on antiviral medications from a wide range of sources that can effectively inhibit Mpro. Antiviral drugs must meet the following criteria: high binding affinity for the target protein, bioavailability, low manufacturing costs, and minimal risk of side effects. 17 Natural products and herbs have been extensively examined for their antiviral properties and shown to have signicant activity against coronavirus. 18 Michael acceptor inhibitors, a-ketoamides, aldehydes, ketones, and their derivatives, among others, have a strong interaction with Mpro's active sites. 19, 20 Recently, Marković and his co-workers studied, adamantane-isothiourea derivatives have strong inhibition potency towards Mpro and mutation of SARS-CoV-2 Spike protein D614G. 21 Coumarin derivatives and their Pd(II) complexes have stronger binding affinity with Mpro of SARS-CoV-2 compared to the FDA approved drugs. 22 These studies explored that natural compounds and their derivatives potent inhibitors for Mpro of SARS-CoV-2. With human cells, these natural compounds have potent antifungal, antibacterial, anti-inammatory, and antiviral properties. 23 We have screened the natural products from the universal natural product database with the Mpro of SARS-CoV-2. Ionic liquids (ILs) have recently gained prominence in biological applications as solvents, stabilizers, and destabilizers. 24, 25 Because of their antibacterial properties, ILs are used in biomedical applications. 26 The physicochemical attributes of designer solvent ILs, such as viscosity, solubility, toxicity, and hydrophobicity, can all be calibrated. Juhi and her co-workers investigated the antiviral activity of ionic liquids, herbs, and natural compounds against Mpro. As a result, ILs have a higher inhibitory action against SARS-CoV-2. 27 Bio-degradation, toxicity, inadequate physicochemical qualities, and expense have been issued with classic ILs like imidazolium, 28 pyrrolidinium, and pyridinium. 29 Cholinium-based ILs have several specic characteristics, including nontoxicity and biodegradability. 30 Recent developments have utilized cholinium-based ILs for many drug delivery applications due to their biocompatibility and benign nature. 31 Insulin structure is strongly stabilized by cholinium-based ILs which are used in biomedical applications. 32, 33 Cholinium-based ILs are chosen against Mpro based on their antibacterial and anti-inammatory properties. The structure, stability and interactions of protein-ligand complexes were determined using integrated computational methods of molecular docking, atomistic molecular dynamics (MD) simulations, MM/PBSA binding free energy calculations, and ADMET properties. Virtual screening methods are widely used for drug designs and development. [34] [35] [36] [37] Large library of natural products (UNPD) and cholinium-based ILs screened using molecular docking. As a result, we chose the most stable complexes of MD simulations to interpret the stability and interaction mechanism between natural products and ILs. Our research will contribute further insights into the inhibitory potential of natural compounds and ILs. Based on these ndings, we can look for potential inhibitors of the SARS-CoV-2 Mpro. Natural products. The Universal Natural Products Database (UNPD), which was established for virtual screening, consists of 229 358 molecules that were conformationally analyzed. The Xray crystal structure of Mpro was retrieved from the protein data bank (PDB ID: 6xr3) which is shown in Fig. 1 . The receptor grid was created with the OEDOCKING's pdb2receptor tool. The FRED (Fast Rigid Exhaustive Docking) module of OpenEye was used to construct the virtual screen. 38 FRED performs a rigid docking of ligand into the protein active site. To generate an ensemble of poses, FRED performs a systematic and comprehensive search. Omega2 was used to generate numerous conformers for the specied compounds. To exhibit the various conformations of selected library compounds, OMEGA (Open-Eye) was employed. It is a rapid conformer generation tool developed by OpenEye that optimizes geometry using a modi-ed Merck Molecular Force Field (MMFF). 39 During the conformation generation, the OMEGA default parameters were employed. The binding poses that resulted from docking were examined. We chose four compounds (UNPD68644, UNPD19759, UNPD75460 and UNPD23673) from the top hits and performed molecular dynamics simulations on them. The natural products of UNPD68644, UNPD19759, UNPD75460 and UNPD75473 are denoted as NP-Hit1, NP-Hit2, NP-Hit3 and NP-Hit4, respectively. The resulting binding modes of selected compounds were analyzed and discussed. The virtual screening of natural compounds are displayed in Fig. 2 . Ionic liquids (ILs). We have chosen 30 cholinium-based ionic liquids, depending on the different biological activities. The three-dimensional coordinates of these ionic liquids are taken from the PubChem database. 40 The structures were optimized employing the B3LYP and M05-2X methods with 6-31+G(d,p) electronic level of theory. These methods are more suitable for geometry optimization and different kinds of noncovalent interactions. [41] [42] [43] In addition, frequency calculations were used to conrm the global minimum in the potential energy surface (PES). All the ILs were optimized using the Gaussian16 program. 44 The optimized geometries and their energies were tabulated in ESI Table S1 . † We calculated the atomic partial charges before docking using the restrained electrostatic potential (RESP). Molecular docking was performed using Autodock 4.2. 45 The structures of most stable compounds of natural products and cholinium-based ILs are displayed in Fig. 3 . We screened the topmost stable compounds for MD simulation studies based on their binding score. ADMET characteristic analysis used to predict the pharmacokinetic properties of all selected classes of compounds were discussed in this work, such as adsorption, distribution, metabolism, excretion, and toxicity (ADMET). SwissADME was used to predict the pharmacokinetic characteristics of a selected class of ionic liquids and natural compounds. 46 Veber's rule, Ghose's rule, and Lipinski's rule of ve (ROF) were calculated for our selected compounds. The druglikeness can be determined using ROF, which states that the molecular weight, hydrophobicity, number of hydrogen bond (nHB) acceptors, and nHB donors should not exceed 500 kDa, 5, 10, and 5, respectively, to be considered a possible drug-like molecule. Veber's rule said that the topological polar surface area (TPSA) and the number of rotatable bonds should not exceed 140Å 2 and 10, respectively, to be considered a prospective drug-like molecule. Another foremost aspect of drug design is determining the toxicity of a specic class of chemicals. The computational approach to toxicity evaluation aids in calculating the toxic dose level in an animal model. It decreases the number of animal model trials and reduces the probability of failure in the experimental technique. We ran MD simulations for Mpro complexes of natural products and cholinium-based ILs using Gromacs 2016.3. 47 The complexes of Mpro and natural products were placed in a cubic box with a volume of 9.8 nm 3 . The system is solvated with the SPC/E water model and four Na + ions were added to maintain neutrality. 48 For Mpro and ligands, force led parameters of Amber ff99SB-ILDN and general amber force eld (GAFF) were used, respectively. 49, 50 The atomic partial charges for ligands were calculated the B3LYP/6-31+G** level of electronic calculations using the RESP method. 51 For all three dimensions, periodic boundary conditions (PBC) were used. The steepest descent algorithm was used to achieve energy minimization of the system. At 300 K and 1 bar pressure, the system was equilibrated with isothermal-isochoric (NVT) and isothermalisobaric (NPT) ensembles using velocity-rescaling and Berendsen coupling, respectively. For long-range electrostatic interactions, the Particle Mesh Ewald technique has been utilized with a cut-off distance of 1.2 nm. 52 The LINCS algorithm was utilized to limit all hydrogen and heavy atom bonds. 53 The integration time step of 2 fs was applied using the leap-Frog algorithm. 54 The molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) approaches were used for binding free energy calculations. 55 We solvated Mpro with a 3 M concentration of choliniumbased ILs to better understand its stability in the presence of ILs medium. At 3 M concentration, the effect of stabilization and destabilization of protein was determined. 56 The GAFF parameters were used for cholinium-based ILs. Further, the above-mentioned MD protocol was applied to the system. The basic analysis was done with the gromacs suite of programs. VMD and pyMOL tools were used to visualize and analyze the trajectories. 57, 58 The production run was performed for 200 ns, and the time series with energies were plotted to establish the equilibration. The importance of natural compounds as antiviral agents has already been explored. We performed molecular docking investigations involving naturally occurring compounds (UNPD), and ILs in addition to the control drugs for treating SARS-CoV-2 in this work. On protein targets, the binding methods and binding sites of ligand were determined. The most stable conformations were sorted according to their highest binding energy, and the results were then used to interpret the results individually. Docking studies Natural compounds. We performed molecular docking for a large library (i.e., 229 358 molecules) of natural compounds against Mpro. Table 1 shows the results of the top hits and their interactions. NP-Hit1, NP-Hit2, NP-Hit3 and NP-Hit4 are the natural compounds with a binding score range of À11.87 to À10.91 kcal mol À1 . Complex formation is associated with negative binding energy, implying that the complex formation is spontaneous. These four compounds exhibit a higher affinity for binding than lopinavir (À9.9 kcal mol À1 ), a potent Mpro inhibitor that has been explored in vitro and preclinical studies. Molecular docking, ADMET and MD simulations on Mprolopinavir complex is shown in ESI. † In molecular docking, FRED and Autodock 4.2 were used for natural compounds and cholinium-based ILs, respectively. Natural compounds interact with the active sites (His41 and Cys145) of the Mpro, as shown by the docked poses. The docked conformations of Mpro with natural compounds are shown in Fig. 4 . NP-Hit1 compound forms ve H-bonding interactions with Cys44, Glu143, His164, Asp187, and Gln189. NP-Hit2 also has six H-bonding interactions with Tyr54, Phe140, Ser144, His163, Glu166, and Asp187. The amino acids involved in the H-bonding interactions of NP-Hit3 are Phe140, Leu141, Asn142, His163, and Arg188. The residues in contact with the NP-Hit4 through H-bond are Thr25, Cys44, Phe140, Ser144, His163, Glu166, His172, and Arg188. Aside from this, the compounds NP-Hit1, NP-Hit2, NP-Hit3, NP-Hit4 forms p-p interactions with His41. NP-Hit4 shows anion-p stacking interactions with the amino acid residue Gln189. Ionic liquids (ILs). We have selected the biocompatible cholinium-based ILs for the screening against Mpro of SARS-CoV-2. The docking energies for all the 30 compounds were shown in ESI (Table S3 †) . The conformers for four ILs with various anion moieties, choline naproxenate, choline ketoprofenate, choline linoleate, We have calculated the ADMET properties for the stable compounds. The ADMET analysis evaluates drug-likeness properties of natural compounds as well ILs which is displayed in Table 2 . The drug-likeness parameter is inuenced by the way the drug is administered. Aqueous solubility and intestinal permeability are associated with these pharmacokinetic parameters. Lipinski's rule of ve (RO5) aids in predicting the medicinal and combinatorial chemistry of a given class of compounds. The predicted value of various properties such as size, lipophilicity, polarity, etc. The value of molecular weight describes the size of the molecule. The log P value represents the lipophilicity of the molecules, which is connected to the drug molecule's solubility in an aqueous environment. The RO5 implies that when there are more than 5 H-bond donors, 10 H-bond acceptors, the molecular weight is greater than 500 Da, and the computed log P is greater than 5, poor absorption is more likely. TPSA, number of rotatable bonds, and molar refractivity were calculated and applied Veber and Ghose rules to assess drug likeness. PAINS (Pan-assay interference) analysis was carried out on the lead compounds generated by the virtual screening and has shown distinct structural features, which can result in false-positive results in some cases when used in virtual screenings and can interact non-specically with random targets rather than specic targets. As a result, it's critical to rule out compounds having PAINS characteristics to avoid false positives. Based on the ADMET results, natural compounds have high GI (Gastro-Intestinal tract) absorption. Furthermore, ionic liquids follow the Lipinski rule of ve and have high absorption in the GI system. From Table 2 , we conclude that our natural compounds and cholinium-based ILs have good drug-like properties. Hepatotoxicity and carcinogenicity are found inactive in both natural compounds and ionic liquids. Hepatotoxicity refers to liver damage caused by chemicals. The liver is a vital organ for converting and removing toxins, and it is vulnerable to their toxicity. Structural stability of Mpro with natural compounds and ILs. To understand the structural stability and binding interactions of the top-ranked docked complexes, MD simulations were performed. The structural stability and compactness of Mpro complexes were calculated using root mean square deviation (RMSD), and radius of gyration (R g ) analysis. The RMSD of C-alpha atoms was calculated for Mpro with the presence of a water medium for the reference structure. In Fig. 6 RMSD plots reveal that all the complexes are stabilized with the crystal structure of Mpro during the MD simulations. The average RMSD values of Mpro with natural compounds and ILs were displayed in ESI (Table S4 and Table S5 †). The natural compounds NP-Hit3 and NP-Hit4 have observed the lowest RMSD values of 1.77 and 1.78Å, respectively. While the other compound NP-Hit1 has a higher RMSD value of 2.04Å. We considered IL medium for the stabilization of Mpro structure, cholinium-based ILs used for various drug delivery applications. The RMSD plots of Mpro structure in the presence of IL medium are shown in Fig. 6C . The average RMSD values of Mpro in the IL medium are displayed in ESI (Table S5 †) . Cholinium-based ILs were tightly bound with the protein and conserved the structural stability of Mpro. The system density is represented by the radius of gyration (Rg), which has an impact on the folding rate and protein compactness. Rg was used to determine the compactness of the complexes, and it was found that Rg of all systems was consistent with the RMSD system. This shows that the protein remained stable and compacted throughout 200 ns. The radius of gyration (Rg) of protein and natural compound complexes was found to be between 2.20 and 2.21 nm. Also, the radius of gyration for the protein and ionic liquids complexes is found to be between 2.16 and 2.19 nm. Fig. 6B and D reveal the compactness and unfolding of Mpro structure binding with the natural products and ILs. There is no large difference in the average Rg values. Thereby, we conclude that both natural compounds and ILs maintain structural stability. Also, the exibility of C-alpha atoms of Mpro is plotted in ESI (Fig. S3 †) . Thereby, the exibility of the Mpro structure is minimized in the presence of natural compounds and ILs. All these results conclude that both natural products and ILs equally stabilize the Mpro structure. Further, the average number of H-bonding interactions calculated for Mpro-natural compounds and Mpro-ILs are is shown in ESI Table S6 . † NP-Hit4 compound form $2 H-bonding interactions with Mpro, whereas compound NP-Hit3 has very less probability of forming H-bonding interactions. Compound NP-Hit1 and NP-Hit2 are stabilized by one H-bonding interaction with Mpro. On the other hand, 3 M concentration of ILs forms more number of Hbonded interactions. Choline ketoprofenate and naproxenate form more H-bonding interactions, while linoleate forms less number of H-bonds. It is interesting to note that the stability of Mpro-ligand complex highly depend on the number of Hbonded interactions with natural compounds and ILs. We computed the spatial distribution functions (SDF) of the individual ion species surrounding the protein to investigate local ion coordination. We calculated the SDF of choline cation with naproxenate, ketoprofenate, linoleate, and geranate anions to better understand the coordination of cholinium-based ILs around Mpro, as shown in Fig. 7 . It is found from calculations that the naproxenate anions are strongly accumulated around the protein surface than the choline cations. The close inspection reveals that the anions form H-bonding interactions with the protein structure and choline cation interacted with the negatively charged amino acids. The density of all the other anions is higher around the Mpro except linoleate. It has fewer interactions with the structure of the Mpro. shows that the binding energy of the compounds NP-Hit2 (À25.6 kcal mol À1 ) and NP-Hit3 (À25.3 kcal mol À1 ) revealed more binding affinity than the other compounds, whereas the binding energy of NP-Hit4 is À24.1 kcal mol À1 . The binding free energies of natural products and cholinium-based ILs with Mpro were shown in ESI Table S7 and S8. † The binding energies with time series are plotted in ESI (Fig. S4 †) . A comparison of MM-PBSA calculations to molecular docking revealed greater binding energies, indicating the compounds will likely be potential hits. The binding free energy of NP-Hit1 is À19.2 kcal mol À1 , which is very less compared to other compounds. Also, the GI absorption is very less. It concludes that the other three compounds NP-Hit2, NP-Hit3, and NP-Hit4 have a higher binding affinity, and hence, we can consider these are potential inhibitors against Mpro of SARS-CoV-2. Using MM-PBSA, we performed a per-residue decomposition analysis to understand the impact of appropriate amino acid interactions on binding free energy. The hot spot residues that contribute to protein-ligand binding interactions were determined using a free energy decomposition approach. For each system, we extracted 2000 snapshots from the MD trajectories and estimated the average ligand-residue interaction energies (DE lig-res ). The hot spot residues have (DE lig-res ) greater than À1.0 kcal mol À1 which is shown in Table 3 . These residues interact with the ligands and form stronger H-bond and hydrophobic interactions. The most signicant hotspot residue is Thr25, Leu27, His41, Met49, Cys145, Met165, and Gln189. Natural compounds NP-Hit2 and NP-Hit3 have strongly interacted with the hot spot residues. These results indicate that NP-Hit2 and NP-Hit3 are potent inhibitors for the Mpro of SARS-CoV-2. SARS-CoV-2 is a highly contagious disease that has spread over the world, increasing the fatality rate. Seeking new antiviral drugs for SARS-CoV-2 is a top priority. Computer-aided drug design techniques can assist in identifying a viable COVID-19 medication and management at a lower cost and in lesser time. Existing drugs have adverse effects, and the mechanisms of their interactions are yet unknown. We have selected topranked natural compounds and a series of biocompatible ILs for the drug-likeness study. The top-ranked natural compounds and ionic liquids were screened using docking studies based on their binding energies. For both natural compounds and ionic liquids, ADMET properties are used to study the drug-likeness attributes. For the four natural compounds whose binding affinity is À11 kcal mol À1 and also for the ionic liquids the binding affinity ranges from À10 to À9 kcal mol À1 are chosen for further screening in which MD simulation was used. As a result, protein-ligand complexes remain stable throughout the simulation time. The binding free energy for the top-ranked natural compounds was calculated using the MM-PBSA method. The most signicant hotspot residues Thr-25, Leu-27, His-41, Met-49, Cys-145, Met-165, and Gln-189 have made signicant contributions to the protein-ligand binding, which may be useful in drug development for Mpro target interactions. In conclusion, the natural compounds of NP-Hit2 and NP-Hit3 have higher binding affinity and obey the ADMET properties. Additionally, NP-Hit2 and NP-Hit3 interact with the cold spot residues, which is resist the mutations of Mpro. The active site residues His-41, Phe-140, Cys-145, Glu-166, and His-172 have low mutation frequencies, cold spot residue interactions with natural products will be studied in the future. Also, the ionic liquids choline ketoprofenate and choline geranate strongly enhance the stability of Mpro. In conclusion, we screened topranked natural compounds among 0.2 million compounds and top-ranked ionic liquids have made signicant contributions to the protein-ligand binding, which may be useful in drug development for COVID-19. There are no conicts to declare. Evidence-Based Complementary Altern Biological Activity of Ionic Liquids and Their Application in Pharmaceutics and Medicine KP and SMER thank the SRM Institute of Science and Technology (SRM-IST) Research Fellowship for his research work. MP thanks the Department of Science and Technology-Science and Engineering Research Board (DST-SERB), New Delhi, India for the nancial support. The authors also thank SRM-IST for providing the supercomputing facility and nancial support.