key: cord-0768218-gf9inovy authors: Mahmud, Shafi; Uddin, Mohammad Abu Raihan; Paul, Gobindo Kumar; Shimu, Mst Sharmin Sultana; Islam, Saiful; Rahman, Ekhtiar; Islam, Ariful; Islam, Md Samiul; Promi, Maria Meha; Emran, Talha Bin; Saleh, Md Abu title: Virtual screening and molecular dynamics simulation study of plant-derived compounds to identify potential inhibitors of main protease from SARS-CoV-2 date: 2021-02-01 journal: Brief Bioinform DOI: 10.1093/bib/bbaa428 sha: 3ddfc507b4a26bc1bf542dd3c77da49a221d2a63 doc_id: 768218 cord_uid: gf9inovy The new coronavirus (SARS-CoV-2) halts the world economy and caused unbearable medical emergency due to high transmission rate and also no effective vaccine and drugs has been developed which brought the world pandemic situations. The main protease (M(pro)) of SARS-CoV-2 may act as an effective target for drug development due to the conservation level. Herein, we have employed a rigorous literature review pipeline to enlist 3063 compounds from more than 200 plants from the Asian region. Therefore, the virtual screening procedure helps us to shortlist the total compounds into 19 based on their better binding energy. Moreover, the Prime MM-GBSA procedure screened the compound dataset further where curcumin, gartanin and robinetin had a score of (−59.439, −52.421 and − 47.544) kcal/mol, respectively. The top three ligands based on binding energy and MM-GBSA scores have most of the binding in the catalytic groove Cys145, His41, Met165, required for the target protein inhibition. The molecular dynamics simulation study confirms the docked complex rigidity and stability by exploring root mean square deviations, root mean square fluctuations, solvent accessible surface area, radius of gyration and hydrogen bond analysis from simulation trajectories. The post-molecular dynamics analysis also confirms the interactions of the curcumin, gartanin and robinetin in the similar binding pockets. Our computational drug designing approach may contribute to the development of drugs against SARS-CoV-2. top three ligands based on binding energy and MM-GBSA scores have most of the binding in the catalytic groove Cys145, The world is currently facing a miserable situation due to the outbreak of the COVID-19 (coronavirus diseases 2019) viral pandemic which occurred by a novel coronavirus namely SARS-CoV-2 (severe acute respiratory syndrome coronavirus-2) [1] . Since the appearance of the first confirmed case in late December of 2019 at Wuhan, capital of Central China's Hubei Province [2] , this viral pneumonia has infected 35 659 007 people with 1 044 269 confirmed dead till 7 October 2020 [3] . Due to the pandemic, most of the countries were forced to adopt a lockdown mode along with social distancing and quarantine strategy triggering dilapidated economic fallout and human suffering [4] . This SARS-CoV-2 is a genus of beta (β)-coronavirus, which enveloped with a single-stranded RNA virus [5] . Six types of coronaviruses are pathogenic in the human body, most of them including, HCoV-229E, HCoV-HKU1, HCoVNL63 and HCoV-OC43 are less pathogenic and responsible for causing the common cold, but in China and Saudi Arabia, the severe infectious SARS-CoV and MERS-CoV have shown higher pathogenicity in the year of 2002 and 2012. The case fatality rate was relatively higher the SARS (∼10%) and MERS (∼35%) than COVID-19 (0.2-7.7%) [6] where its basic reproduction number (R0) (∼5.7) of COVID-19 [7] is higher compared to the SARS (∼2-3) and MERS (<1) [8] . The SARS-CoV-2 genome has 79% similarity with SARS-CoV-1 and 50% with MERS-CoV [9] , which declared a successor of SARS-CoV-1 by the U.S. National Institutes of Health [10] . The genomic RNA of Coronaviruses ranges from 26-32 kb containing a minimum of six open reading frames [8] and two overlapping polyproteins (pp1a and pp1ab) are encoded by the first ORF (ORF1a/b) of its two-thirds of the genome length [11] . The main protease (M pro ) also is known as 3CL pro (3chymotrypsin-like cysteine protease CCP), which along with the aid of papain-like protease (PL pro ) proteolytically cleaves these viral polyproteins into 16 non-structural proteins (nsp1-16) [12] . The M pro cleaves polypeptide sequence after glutamate residues which makes it an ideal candidate due to its substrate specificity [8] . The polyprotein cleavage sites for M pro in MERS-CoV, SARS-CoV, SARS-CoV-2 also exhibit similar kinds of substrate specificity and recognition [13] . This is a vital step in the interim of viral replication; the enzymes like RdRp-(RNA dependent RNA polymerase) or nsp13 are essential for the replication process, but it cannot fully function without prior photolytic release [8] . These above-mentioned non-structural proteins enjoin the production of sub-genomic RNAs that eventually translated into an envelope, membrane, spike and nucleo-capsid proteins as well as other accessory proteins [11] . Active SARS-CoV-2 M pro is a cysteine protease homodimer that combines two protomers stationing nearly perpendicular to one another [11] . Each one of the protomers contains three domains, and features a substrate-binding noncanonical catalytic dyad (His41 and Cys145) that is located in a cleft within domains I and II (residues 10−99 and 100−182, respectively) [14] . Domains I and II consist of chymotrypsin (which like double βbarrel fold), whereas domain III (198−303 amino acid residues) is mainly comprised of five antiparallel α-helices [12] . The M pro is conserved among coronaviruses [11] . It shows ∼99% identity with BatCoV RaTG13 M pro , ∼96% with the previous SARS-CoV M pro and ∼ 50% only with MERS-CoV M pro through amino acid sequence alignment [8] , in contrast, SARS-CoV-2 M pro plays a pivotal role in arbitrating viral replication and transcription, and inhibition of its activity is expected to block the viral replication, and maturation by providing a key enzyme [11] . SARS-CoV-2 M pro has a particular cleavage specificity unlike any other human proteases [5] . Therefore, this enzyme is an attractive therapeutic target for CoVs than others [8] . Moreover, synthesis and design of new antiviral drugs can be aided through plant-based phytochemicals with more efficacy and specificity. For example, isoquinoline types alkaloid, emetine is widely used as an amoebicidal drug, quinone as well as cancer drug paclitaxel [15, 16] . From the early 1980s to till the day, most drug candidates developed from plant-based natural compounds [17] . Moreover, the use of traditional Chinese medicine Chonqing, in Chinese COVID-19 patients, suggests that plant-based drug molecules can act as an effective solution as complementary or addition with therapeutic drugs [18] . Furthermore, effective and greater antiviral activity was observed for numerous plant-derived compounds, for example, SARS [19] , SARS-CoV-2 [20] and Chikungunya [21] . Computer-aided drug design would be a great weapon in this regard [22] . It is an efficient tool in the searching tool for promising drug candidates in a very short time and costeffective way [23] . Here an in-depth literature review was done to make a library of plant-derived natural compounds from commonly used medicinal plants in the Asian region. We further used a virtual screening technique to shortlist the potent phytochemicals with inhibiting the activity of SARS-CoV-2 M pro . Top hit molecules derived from these were then docked against the target enzyme. A molecular dynamics study was preferred due to the lack of reliance on providing the precise binding mode, binding energy, entropy energy, constant motion effect and solvation effect through molecular docking [24] . We propose that these compounds could be further potential therapeutic candidates for in vitro and in vivo antiviral studies followed by a clinical trial of SARS-CoV-2. The three-dimensional structure of the M pro enzyme of SARS-CoV-2 (PDB ID: 6LU7) [25] was taken from the protein data bank [26] . Protein preparation wizard of the Schrödinger suite version 2020-3 was subjected for preparation [27] . The water molecules were removed beyond 5.0 Å from het groups. The protein structure was preprocessed by adding the hydrogen with assigned bond orders and creating disulfide bonds, including the conversion of selenomethionines to methionines whereby it is applicable. H-bond networks were optimized, and protonation states were yields at pH 7.0. Finally, by applying the OPLS3e force field, the energy minimization process was done [28] . Afterward, keeping the active site residues His41 and Cys145 specific as centroid the grid was generated using the 'Receptor Grid Generation' of the Schrödinger suite. By literature search, about 3063 compounds (Supplementary file 1) were enlisted from the plant-based phytochemicals. The Pub-Chem database [29] was utilized for compound extraction and the Ligprep module of the Schrödinger suite [30] was employed by applying default parameters for ligand molecule preparation. Furthermore, Epik version v5.3 was used to obtain multiple states of the ligand molecules at pH 7.0 ± 2, and for the likelihood of reliability in the biological condition, the high energy ionization/tautomer states were removed [31, 32] . The Qikprop version 6.5 program was run [33] before employing structure-based virtual screening, to screen the plant compounds by applying Lipinski's rules of five [34] . The Glide program and its virtual screening workflow process were applied, including three docking protocols; high throughput virtual screening or HTVS, Standard Precision or SP module, and Extra Precision or XP module [35, 36] . To obtain the best compounds, we follow the previous literature as earlier published [37] . The HTVS was used to dock each ligand to the receptor, which generates one pose. Though the SP docking protocol provides a good scoring function retaining the good scoring states [38] , about 50% of the total plant-derived compounds were transferred from HTVS to SP which helps to find out the false-positive results. Moreover, about 30% of SPs total ligand molecules were processed to XP, where the XP provides the best scoring states [35, 36] . Three poses for each ligand were generated by the XP. Furthermore, only 10% of total compounds were specified postprocessed through MM-GBSA (molecular mechanicsgeneralized born surface area) for the accuracy of pose ranking and final selectivity [22, 39] . The Prime MM-GBSA module from the Schrodinger software package was utilized to calculate the binding affinity. The higher degree of rigidity of the ligand attached protein is indicated by the higher negative MM-GBSA value. The Prime MM-GBSA process consists of three different approaches: OPLS molecular mechanics energies, an SGB solvation model and a non-polar solvent. The binding free energies were calculated from the following equations: where G = EMM + VSGB + GNP Therefore, to perceive their rigidity along with motion and structural stability in simulation conditions, the best three ligands are selected for further processing. The pharmacological and carcinogenic properties of the compounds were assessed with the aid of admetSAR [40] , Swissadme [41] and PkCSM [42] webserver. The canonical smiles of the screened complex were used as an entry system of the complex. The dynamics simulation of the screened ligand molecules was conducted to analyze the conformational behavior and protein stability of the complex. The YASARA software package version 20.1.1 [43] was used to conduct simulation for three complexes where the AMBER14 force field was used [44] . The NPT ensemble method was used in this simulation system and also, the Berendsen thermostat method was applied to control the temperature of the systems. To calculate long-range electrostatic interactions, the particle-mesh Ewald method [45] was employed, whereas a cut-off radius of 8 Å was considered [46] for short-range van der Waals and Coulomb interactions. The cubic cell, which was employed for simulation was 20 Å bigger than the drug-protein complexes in all cases and a periodic boundary condition was maintained. The system was neutralized by the addition of 0.9% NaCl at 298 k temperature. The initial energy minimization process of the systems was incorporated by the simulated annealing method by applying the steepest gradient algorithms [47] . The simulation was carried out with a time step of 1.25 fs and the simulation trajectory was saved for every 100 ps. Finally, the molecular dynamics simulation was carried out for 50 ns, and trajectories were used to analyze root mean square deviations (RMSD), root mean square fluctuations (RMSF), solvent accessible surface area (SASA), radius of gyration (Rg), and hydrogen bond number of the system [22, 48, 49] . The result of virtual screening revealed the identification of about three compounds as a potential inhibitor as they exhibit non-covalent interactions with the amino acid derived from the catalytic residue of the M pro of SARS-CoV-2. The large phytochemical dataset was initially screened in Qikprop module where 2117 compounds were passed to the next passed, following Lipinski filtering where 1432 compounds were screened. The high throughput virtual screening module then filter and 716 phytochemicals were again filtered in SP or standard precision program. The XP or extra precision program was allowed for more accurate screening (Figure 1 ). The glide ligand efficiency score, SP GScore, Glide evdw, Glide ecoul, Glide energy, Glide emodel score were taken into consideration for screening the ligand molecules from the dataset. The lower energy score from these descriptors utilized to rank the top 19 molecules from the dataset. The XP Gscore was found higher in Quercetin as Moreover, for better approximation, the poses of the ligand molecules generated from the XP module was additionally processed with MM-GBSA. Hence, among 19 compounds, only three compounds have been selected by observing the binding energy ( G bind) with the catalytic region of the M pro . Although, quercetin, myrecetin, taxifolin, trans-caftaric acid, luteolin 7 methyl ether exhibit better energy score but they had higher energy score MM-GBSA. So these molecules were not taken Table 1 , binding interaction in Figure 2 and the interacting residues of the M pro with ligand molecules with the active site residues are shown in Table 2 . The pharmacokinetics and toxicity properties of the ligand need to be assessed to ensure the safety and efficacy level of the hit molecules. Moreover, CNS permeability, blood-brain barrier absorption, p-glycoprotein inhibition, hepatotoxicity, carcinogenicity, CYP inhibition of the lead molecules were checked. The CNS permeability determines the ability to permeable through the blood-brain barrier where CNS > −2 considered to penetrate the central nervous system. Among the three compounds, no carcinogenic and toxicity profiles were observed in carcinogenicity and AMES toxicity assessment. The selected compounds exhibit a positive response in Lipinski rule of five where the number of hydrogen bond donors, acceptor and surface area of the ligand molecules were explored. Also, the molecular weight (MW) of the compounds was (368.38, 396.439 and 302.238) g/mol, respectively, for curcumin, gartanin and robinetin which was good for heat molecules consideration as higher MW in compounds violates the rule of Lipinski. The hydrogen bond donors were found as 2, 4 and 5 whereas hydrogen bond acceptors were found as 6, 6 and 7, respectively, for curcumin, gartanin and robinetin. Although, some carcinogenic properties and the violation of Lipinski ' s rule of five was observed for some cancer drugs but some deviation in the ADMET properties might be acceptable if the drug molecules exhibit desired pharmacological properties Table 3 . The molecular dynamics simulation of the docked complex was employed to validate the docking study, and also the dynamic motion of the docked complex was analyzed to understand their degree of stability. The RMSF of the c-alpha atoms is illustrated in Figure 3 , where curcumin and M pro complex had initial stability and fluctuated a little bit after 10 ns. This complex maintained a stable RMSD profile till the rest simulation periods. Also, the M pro protein from SARS-CoV-2 and robinetin complex had a similar RMSD trend till 35 ns and thereafter increased a little bit but these complexes did not over-fluctuate which demonstrates rigid conformation. Therefore, gartanin and M pro complex had a sharp increase in 30-35 ns, followed by a lower trend like the starting phase. Interestingly, all three complexes did not disclose RMSD descriptors over 2.5 Å which validates the rigid conformation of the drug complexes. From Figure 3 , it was also observed that robinetin complex had a similar Rg profile from 0 to 25 ns and thereafter, decreased Rg descriptor which indicates the tight packaging system of the complex as the Rg of the drug complex indicates the compact nature of the protein. Moreover, the degree of protein folding and unfolding greatly depends on the value of the Rgs. Moreover, gartanin and the M pro complex had lower Rg profile from 0 to 25 ns and a higher rise of Rg was observed from 25 to 45 ns, which indicates loose packaging of the system. The curcumin complexes had a moderate Rg profile compared to the other two complexes, which illustrates the less mobile nature. Moreover, the solvent-accessible surface area of the complexes was analyzed from the simulation trajectories to assess the complex volume change through the simulation trajectory. Figure 3 demonstrates that the gartanin complex had a lower SASA profile at the starting phase and thereafter increased the SASA descriptors, which indicate expansion in the protein surface area. Therefore, this complex truncated its volume by a significant degree for the rest of the simulation trajectory. Interestingly, robinetin and curcumin complexes had a similar SASA trend which is higher than gartanin although little fluctuations were observed from 35-40 ns simulation time. On the other hand, the hydrogen bond has a crucial role in contributing stability in drug-protein complex and molecular recognition. All three drug-protein complexes remained stable in the simulation trajectory which indicates the rigidness of the complexes (Figure 3 ). The flexibility among the amino acid residues can be illustrated through the root mean square fluctuation or RMSF profile. From Figure 3 , it can be observed that most of the amino acid residues from three complex had lower flexibility as they did not have higher RMSF values except, Ser1(helix-strand), Gly2(helix-strand), Phe3(helixstrand), Lys12(helix-strand), Thr24(beta-turn), Glu47(betaturn), Glu55(helix-strand), Ile59(helix-strand), Asn72(betaturn), Asn142(beta-turn), Phe140(beta-turn), Tyr154(beta-turn), Ala193(gamma-turn), Gln189(gamma-turn), Arg222(beta-turn), Gln244(helix-strand), Arg279(helix-strand), Ser301(beta-turn), Gly302(beta-turn), Val303(beta-turn), Thr304(beta-turn), Phe305 (beta-turn) and Gln306(beta-turn) residue. The last snapshot from the molecular dynamics simulation trajectory was analyzed again to find out a change in interaction dynamics of the docked complex. The results were tabulated in Table 4 , where several interactions ( Figure 4 ) were observed the same as docked complex which includes less conformational variation along with structural compactness. The curcumin and M pro of the SARS-CoV-2 structure were stabilized by four hydrogen bonds at Arg188, Gln189, His163 and Thr190. Among these four hydrogen bonds, Thr190 was present in the pre-MD structure and this hydrogen bonding at the active groove of the M pro enzyme may be responsible for favorable binding energy in docking. Moreover, one pi-sulfur bond at Cys145 (active site), one alkyl bond at Ala191 (active site), two pi-alkyl bonds at His163 (active site), Leu167 (active site) were also observed in post MD structure of the curcumin complex. Like curcumin, gartanin complex also had four hydrogen bonds but they are positioned among different residues at Gln189 (active site), Gln192 (active site), His164 (active site), Arg188 residues. Also, two alkyl bonds at Met165 (active site), Cys145 (active site), and five pi alkyl bonds at His41 (active site), His163(active site), His164(active site), Phe181, Pro168 (active site) stabilized the gartanin complex. Maximum non-bonded interaction at the active sites was also observed for robinetin complex where four hydrogen bonds at Cys44, Asn142, Thr26, Asp187, and two hydrophobic bonds at Met49 and His41 was present. The novel coronavirus namely SARS-CoV-2 has created a pandemic situation due to the high rate of mortality, and no effective drugs or vaccine to treat against SARS-CoV-2 shaped the new global disaster. Although several clinical trials have been undergoing, the drug development methods are time-consuming and costly, there is a need for fast and effective development of active antiviral agents. Conversely, computer-aided drug design may assist the researcher to find the new therapeutic agent against SARS-CoV-2 due to its rapid and accurate screening capability from a vast small molecule library [50] [51] [52] . Due to having a large impact on the function of SARS-CoV-2, the M pro has become the best target for different therapeutic tactics. It contains three different domain regions; Domain I which covers 1-99 residues, Domain II consists of 100-182 amino acid sequences and Domain III covers 198 to the last residues 303 [11] . It has been observed that the activation through dimerization mechanism through Cys145 and His 41. It is significant that through the drug development technique those catalytic residues would be the best stage for the development of strong inhibitors of M pro [11] . In conjunction with the catalytic sites, there are two more subsites known as S1 and S2, and three more shallow subsites; S3, S4 and S5, where the S1 consist of His163, Glu166, Cys145, Gly143, His172, Phe140 residues, and S2 contains Cys145, His41, and Thr25; whereas the other three subsites consist of Met49, His41, Met165, Glu166 and Gln189 residues [12] . The identification of potent inhibitors from virtual screening workflow allows screening from the diverse large compound library. Additionally, depreciation of the false-positive result in screening modules combinatorial approaches were considered (e.g., HTVS, SP, XP, along with MM-GBSA scoring). These lead compounds identification criteria followed with the addition of binding in the active groove of the targeted protein to explore the possible binding sites in the active points and as well as the interaction dynamics [53] [54] [55] . Based on our analysis, we can suggest that the top three screened compounds with better energy interacted with a catalytic residue which is a prerequisite for the inhibition. The first candidate curcumin surrounds curcumin moiety consists of feruloyl chromophores which are joined by a methyl group [56] . The curcumin has chemo-preventive and chemotherapeutic activity [57] , anti-inflammation [58] , anti-parasite [59] , and carcinogenicity suppression [60] . This compound binds with the M pro of SARS-CoV-2 with numerous non-covalent interactions, where it creates one hydrophobic bond at the active cavity (Cys145). Also, it shows seven hydrogen bonds with the active site residues; Gly143, Gln189, Thr190, Pro168, Leu141 and Glu166. Additionally, it also forms two hydrophobic bonds with the targeted M pro at the active region; Met165 and Pro168 [61, 62] . Besides, candidate gartanin is widely known as plant metabolites and antineoplastic agents [63] . The gartanin has activity toward Alzheimer's diseases [64] , mTOR pathway as well as anti-cancer agents [65] . This phytochemical compound also showed one hydrogen bond and catalytic amino acid residue from the M pro has hydrophobic bond interaction at His41, one hydrophobic bond, and one pi-Sulfur with the active site residue of Cys145. More significantly, it created six hydrogen bonds in the active part of this protein also; Gly143, Gln189, Asn142. Additionally, two more hydrophobic interactions stabilized the complex by making contacts at Met49 and Met165 [61, 66] . The robinetin has anti-mutagenesis effect [67] , antitumorgenicity [68] , atheroprotective effect [69] in enzymatic and protein assay. The compound robinetin showed one strong hydrogen bond with 3.05 Å at His41 residue along with multiple hydrophobic bonds at the active residue of Cys145. Moreover, multiple hydrogen bonds (Thr26, Asp187 and Met165) at the active groove of the M pro gives the stability of the complex. It also formed two hydrophobic bonds with Met49 correspondingly [61, 70] . Moreover, to understand and endorse the molecular docking study, dynamics simulation was conducted. The RMSD profile from curcumin, gartanin and robinetin did not exceed 2.5 Å as their average were 1.72, 1.873 1.69 Å, respectively, which specifies their overall integrity. The robinetin and curcumin complex were comparatively more inflexible compared to gartanin as this complex had inflexibility in the middle phase of the simulation. On the other hand, RMSF analysis established that Domain II had strict nature whereas amino acid residues from Domain I and Domain III had more flexibility in the helical and turn region, although almost every residue was fixed. Also, the surface area of these complexes was not altered by surpassing the SASA values too much as their average of curcumin, gartanin and robinetin were 14175.2 Å2, 14063.98 Å2 and 14161.77 Å2, respectively. The hydrogen bonding assessment also aligned with the result of other descriptors from molecular dynamics simulation as they did not over-fluctuate across the simulation trajectory. Furthermore, superimposition between the docked complex and post-MD docked complex after 100 ns molecular dynamics simulation revealed that both structures had a similar binding position in the active points as the top three complex curcumin, gartanin and robinetin had RMSD value of 1.76, 2.03 and 1.43 Å, respectively ( Figure 5) . We also took a snapshot (Figures 6-8 ) from 0, 25, 50, 75 and 100 ns from molecular dynamics simulation trajectory for the top three phytochemicals but no drastic change was observed for their binding pose. The combinatorial docking and molecular dynamics approach suggest three lead compounds may interfere with the function of the M pro of SARS-CoV-2; however, these plant-derived phytochemicals need to be tested more in the lab to check the efficacy along with inhibitory potential in in vitro condition. In this word, we have employed a computational drug design workflow to recognize potent inhibitors of the M pro from SARS-CoV-2. We have combined the phytochemical dataset consist of over 3000 compounds from Asian plants to investigate their inhibitor potentiality. The virtual screening procedure along with MM-GBSA approaches aids in shortening the list from over 3000 to three potential lead molecules, curcumin, gartanin and robinetin. Furthermore, binding pose and interactions from the docking study were further evaluated through a molecular dynamics simulation study where multiple descriptors from simulation trajectories confirm their stability. We also found that the catalytic residue of the M pro , Cys145 and His41 binds with the drug molecules and their existence was also confirmed in the post-MD structures. However, toxicity and pharmacological estimation of the drug molecules confirms better absorption and metabolism profile along with no toxicity probabilities. Since this study is solely based on multiple computational tools and simulation studies, it requires further evaluation in the lab, and also this study may be helpful for future researchers to work with precise target molecules from an extensive library in search of effective drug development against COVID-19. • The main protease from SARS-CoV-2 can be targeted as an attractive key protein as it has a central role in activating viral replicase through posttranslational modification. The target protein was employed to search for potent inhibitors through virtual screening and molecular dynamics. • The compound dataset from Asian plants was retrieved by data mining and literature review. These compound lists were used to screen with the combination of molecular docking and MM-GBSA approaches in the virtual screening process. • The top three drug candidates were further assessed in ADMET filtering where no toxicity and positive pharmacological properties were found. • The molecular dynamics simulation of the docked complex gives insights into their inflexibility and stability of binding. These potential compounds may be used as a potent drug candidate after further evaluation in the biological lab. Supplementary data are available online at Briefings in Bioinformatics. Identification of SARS-CoV-2 cell entry inhibitors by drug repurposing using in silico structure-based virtual screening approach The continuing 2019-nCoV epidemic threat of novel coronaviruses to global healththe latest 2019 novel coronavirus outbreak in Wuhan, China Cell entry mechanisms of SARS-CoV-2 An investigation into the identification of potential inhibitors of SARS-CoV-2 main protease using molecular docking study RESEARCH high contagiousness and rapid spread of severe acute respiratory syndrome coronavirus 2 The SARS-CoV-2 main protease as drug target SARS-CoV-2, SARS-CoV, and MERS-CoV: a comparative overview Infectious (National Institute of Allergy and, Diseases), et al. Aerosol and surface stability of SARS-CoV-2 as compared with SARS-CoV-1 Structure-based design of antiviral drug candidates targeting the SARS-CoV-2 main protease Identification of chymotrypsin-like protease inhibitors of SARS-CoV-2 via integrated computational approach Substrate specificity profiling of SARS-CoV-2 Mpro protease provides basis for anti-COVID-19 drug design Targeting the dimerization of the main protease of coronaviruses: a potential broad-Spectrum therapeutic strategy Natural product-derived phytochemicals as potential agents against coronaviruses: a review Herbal plants and plant preparations as remedial approach for viral diseases Natural products as sources of new drugs from 1981 to 2014 Clinical features and treatment of COVID-19 patients in Northeast Chongqing In vitro susceptibility of 10 clinical isolates of SARS coronavirus to selected antiviral compounds Scutellaria baicalensis extract and baicalein inhibit replication of SARS-CoV-2 and its 3C-like protease in vitro Deciphering the potential of baicalin as an antiviral agent for chikungunya virus infection Exploring the potent inhibitors and binding modes of phospholipase A2 through in silico investigation Role of computeraided drug design in modern drug discovery Molecular docking in modern drug discovery: principles and recent applications Structure of Mpro from SARS-CoV-2 and discovery of its inhibitors The protein data Bank Protein and ligand preparation: parameters, protocols, and influence on virtual screening enrichments OPLS3e: extending force field coverage for drug-like small molecules PubChem substance and compound databases Epik: a software program for pKaprediction and protonation state generation for drug-like molecules Quantum chemical calculation and binding modes of H1R; a combined study of molecular docking and DFT for suggesting therapeutically potent H1R antagonist Lead-and drug-like compounds: the rule-of-five revolution Flexible ligand docking with glide Extra precision glide: docking and scoring incorporating a model of hydrophobic enclosure for protein−ligand complexes Molecular docking and dynamics study of natural compound for potential inhibition of main protease of SARS-CoV-2 Online publication Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy Assessing the performance of the MM/PBSA and MM/GBSA methods. 6. Capability to predict protein-protein binding free energies and re-rank binding poses generated by protein-protein docking AdmetSAR: a comprehensive source and free tool for assessment of chemical ADMET properties SwissADME: a free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules pkCSM: predicting smallmolecule pharmacokinetic and toxicity properties using graph-based signatures YASARA-yet another scientific artificial reality application The amber biomolecular simulation programs An efficient linear-scaling Ewald method for long-range electrostatic interactions in combined QM/MM calculations Fast empirical p K a prediction by Ewald summation New ways to boost molecular dynamics simulations Assessment of structurally and functionally high-risk nsSNPs impacts on human bone morphogenetic protein receptor type IA (BMPR1A) by computational approach Computational discovery of plant-based inhibitors against human carbonic anhydrase IX and molecular dynamics simulation Online publication Comparative molecular investigation of the potential inhibitors against SARS-CoV-2 main protease: a molecular docking study A dynamic simulation study of FDA drug from zinc database against COVID-19 main protease receptor Repurposing of the approved small molecule drugs in order to inhibit SARS-CoV-2 S protein and human ACE2 interaction through virtual screening approaches Docking and virtual screening strategies for GPCR drug discovery Identification of potential Mpro inhibitors for the treatment of COVID-19 by using systematic virtual screening approach Virtual screening approaches towards the discovery of toll-like receptor modulators The essential medicinal chemistry of curcumin In vivo evaluation of curcumin-loaded nanoparticles in a A549 xenograft mice model Properties of Lewis lung carcinoma cells surviving curcumin toxicity Validity of silver, chitosan, and curcumin nanoparticles as anti-giardia agents Determination of minimum effective dose and optimal dosing schedule for liposomal curcumin in a xenograft human pancreatic cancer model In silico prediction of potential inhibitors for the main protease of SARS-CoV-2 using molecular docking and dynamics simulation based drug-repurposing Identification of high-affinity inhibitors of SARS-CoV-2 main protease: towards the development of effective COVID-19 therapy Gartanin, an isoprenylated xanthone from the mangosteen fruit (Garcinia mangostana), is an androgen receptor degradation enhancer Natural Xanthones from Garcinia mangostana with multifunctional activities for the therapy of Alzheimer's disease The effect of gartanin, a naturally occurring xanthone in mangosteen juice, on the mTOR pathway, autophagy, apoptosis, and the growth of human urinary bladder cancer cell lines Identification of potential inhibitors of SARS-CoV-2 main protease from Aloe vera compounds: a molecular docking study Anti-mutagenesis and anti-promotion by apigenin, robinetin and indole-3-carbinol Effect of ellagic acid and hydroxylated flavonoids on the tumorigenicity of benzo[a]pyrene and (±)-7β, 8α-dihydroxy-9α Atheroprotective effects of antioxidants through inhibition of mitogenactivated protein kinases Crystal structure of SARS-CoV-2 main protease provides a basis for design of improved αketoamide inhibitors The authors receive no funding from an external source. The authors declare no conf lict of interest.