key: cord-0725864-x52u8w28
authors: Dutta, Mycal; Tareq, Abu Montakim; Rakib, Ahmed; Mahmud, Shafi; Sami, Saad Ahmed; Mallick, Jewel; Islam, Mohammad Nazmul; Majumder, Mohuya; Uddin, Md. Zia; Alsubaie, Abdullah; Almalki, Abdulraheem S. A.; Khandaker, Mayeen Uddin; Bradley, D.A.; Rana, Md. Sohel; Emran, Talha Bin
title: Phytochemicals from Leucas zeylanica Targeting Main Protease of SARS-CoV-2: Chemical Profiles, Molecular Docking, and Molecular Dynamics Simulations
date: 2021-08-17
journal: Biology (Basel)
DOI: 10.3390/biology10080789
sha: 2e3a6a009e59bbfe77699a6e46aaca2e3053b856
doc_id: 725864
cord_uid: x52u8w28

SIMPLE SUMMARY: Molecular docking in conjunction with molecular dynamics simulation was accomplished as they extend an ample opportunity to screen plausible inhibitors of the main protease from Leucas zeylanica. The preferential phytochemicals were identified from L. zeylanica through gas chromatography–mass spectrometry (GC-MS). The pre-eminent three identified phytochemicals exhibited toxicity by no means during the scrutinization of ADME/T prominences. Moreover, pharmacologically distinguishing characteristics and the biological activity of the lead phytochemicals were satisfying as an antiviral drug contender. Additionally, the molecular dynamics simulation exhibited thermal stability and a stable binding affinity of the protein–compound complex that referred to the appreciable efficacy of lead optimization. Therefore, the preferable phytochemicals are worth substantial evaluation in the biological laboratory to recommend plausible antiviral drug candidates. ABSTRACT: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a contemporary coronavirus, has impacted global economic activity and has a high transmission rate. As a result of the virus’s severe medical effects, developing effective vaccinations is vital. Plant-derived metabolites have been discovered as potential SARS-CoV-2 inhibitors. The SARS-CoV-2 main protease (M(pro)) is a target for therapeutic research because of its highly conserved protein sequence. Gas chromatography–mass spectrometry (GC-MS) and molecular docking were used to screen 34 compounds identified from Leucas zeylanica for potential inhibitory activity against the SARS-CoV-2 M(pro). In addition, prime molecular mechanics–generalized Born surface area (MM-GBSA) was used to screen the compound dataset using a molecular dynamics simulation. From molecular docking analysis, 26 compounds were capable of interaction with the SARS-CoV-2 M(pro), while three compounds, namely 11-oxa-dispiro[4.0.4.1]undecan-1-ol (−5.755 kcal/mol), azetidin-2-one 3,3-dimethyl-4-(1-aminoethyl) (−5.39 kcal/mol), and lorazepam, 2TMS derivative (−5.246 kcal/mol), exhibited the highest docking scores. These three ligands were assessed by MM-GBSA, which revealed that they bind with the necessary M(pro) amino acids in the catalytic groove to cause protein inhibition, including Ser144, Cys145, and His41. The molecular dynamics simulation confirmed the complex rigidity and stability of the docked ligand–M(pro) complexes based on the analysis of mean radical variations, root-mean-square fluctuations, solvent-accessible surface area, radius of gyration, and hydrogen bond formation. The study of the postmolecular dynamics confirmation also confirmed that lorazepam, 11-oxa-dispiro[4.0.4.1]undecan-1-ol, and azetidin-2-one-3, 3-dimethyl-4-(1-aminoethyl) interact with similar M(pro) binding pockets. The results of our computerized drug design approach may assist in the fight against SARS-CoV-2.

In recent decades, the world has witnessed an unprecedented number of life-threatening human disease outbreaks caused by an array of pathogenic organisms, including several notable viral diseases, such as influenza, chikungunya, Nipah, Zika, and Ebola [1, 2] . However, the ongoing spread of coronavirus disease 2019 (COVID-19) has been exponential. It has already surpassed most previous viral infections in terms of infectivity and has become the center of global attention. Wuhan, a populous Chinese city located in the Hubei province was the first location where this acute respiratory infection was identified in late December 2019 [3] . COVID-19 has taken a significant toll on people worldwide and on 11 March 2020 was declared by the World Health Organization (WHO) a pandemic. This highly contagious infection has had a detrimental impact on the global healthcare management system and, as of 1 February 2021, >100 million confirmed cases have been reported, including more than 2 million estimated deaths worldwide [4] .

SARS-CoV-2 is a pleomorphic, enveloped, nonsegmented, single-stranded RNA betacoronavirus belonging to the Coronaviridae family and features a large genome (27-32 kb) that encodes both structural and nonstructural proteins. SARS-CoV-2 is associated with a higher transmission rate than other well-known human beta-coronaviruses, such as SARS-CoV and Middle Eastern respiratory syndrome coronavirus (MERS-CoV) [5] . The coronavirus main protease (M pro ) is a nonstructural protein that plays a crucial role in protein translation, viral replication, and maturation [6, 7] . In a recent study, Liu et al. confirmed the existence of the M pro (also known as 3CL Pro or chymotrypsin-like protease) enzyme in SARS-CoV-2 [8] . The genome of SARS-CoV-2 encodes pp1a and pp1ab, two large polyproteins, similar to other Coronaviridae genomes [9] . The resulting polyproteins, pp1a and pp1ab, must be cleaved to generate mature nonstructural proteins (nsps) [10] . The large pp1a (replicase 1ab) is generated inside the cell via genomic RNA transcription. Therefore, the inhibition of M pro activity is anticipated to result in the prevention of viral replication, as similar cleavage specificity has not been identified in any human proteases [11] , indicating that the polypeptide cannot be properly cleaved in the absence of M pro . Additionally, M pro has very low cytotoxicity and low similarity with human proteases [12] . Proteins produced by this pathogenic organism have been demonstrated to intervene the host immune response, and M pro enzyme-specific T cells have been encountered in SARS-CoV-2 patients [13, 14] . Therefore, M pro is considered to represent a promising drug target for antiviral drug development.

In general, cytokine production, cell death, inflammation, and other pathophysiological processes are commonly associated with disruptions in redox balance, resulting in oxidative stress during viral infections, which can negatively affect the respiratory tract. Previously, viral replication was strongly correlated with the excessive production of reactive oxygen species (ROS) and a reduction in the components of antioxidant mecha-nisms [15] [16] [17] [18] . Inflammatory reactions are triggered by the COVID-19 infection, resulting in the subsequent release of proinflammatory cytokines, which can cause acute lung damage [19] . Oxidative stress also plays a critical role in the perpetuation of the cytokine storm cycle and is important for blood clotting mechanisms [20] . The observed increase in COVID-19 infection severity in patients diagnosed with chronic diseases has been linked with the poor performance of the antioxidant system, suggesting that antioxidants may represent a prospective therapeutic option for COVID-19 infection [21] .

A close connection exists between the innate immune response and the thrombotic response, and recent COVID-19 clinical data have revealed a correlation between this infection and thrombotic complications, which might result in increased incidence of microvascular thrombosis, venous thromboembolic illness, and stroke. Markers of COVID-19 include thrombotic complications, which are often associated with multiorgan failure and increased fatality [22] .

Recently, phytochemicals have been investigated against different target proteins of SARS-CoV-2 to find appropriate lead compounds for COVID-19 infection. Baicalin, baicalein, 25-hydroxycholesterol, chrysosplenetin, shikonin, panduratin A, and quercetin are some of the examples of plant compounds that exhibited a potential effect against SARS-CoV-2 during in vitro studies. In silico studies were conducted on a broad spectrum for a plethora of medicinal plants and a huge number of compounds were screened. This not only helped indicate the potentiality of natural plant compounds but also reduced the number of tedious and costly wet-laboratory experiments [23] . Therefore, in our current study, we endeavored to assess the roles of phytoconstituents identified from Leucas zeylanica, a medicinal plant belonging to the Lamiaceae family, in the management of COVID-19 infection by employing computational biology approaches. Phytocompounds possess a wide range of pharmacological activities, and traditional healers have employed plants belonging to the Leucas genus to treat various disease states, indicating an immense potential for the discovery of new lead compounds [24] . L. zeylanica is a weed commonly referred to as "Ceylon slitwort" and locally known as "Kusha" [25] . The plant is widely distributed throughout China, India, Bangladesh, Sri Lanka, Thailand, Indonesia, Philippines, Vietnam, Cambodia, Nepal, Myanmar, Malaysia, and New Guinea [26] . A phytochemical screening of an L. zeylanica methanol extract confirmed the presence of alkaloids, flavonoids, tannins, steroids, and glycosides, which contribute to the traditional medicinal properties of the plant [27] . This plant is traditionally used as a vermifuge ingredient in addition to the treatment of burning sensations during urination, scabies, convulsion, fever, jaundice, scorpion and snake bites, colds, rheumatism, roundworm, psoriasis, anorexia, flatulence, colic, and malaria [24, 25, 28, 29] . The antimalarial drugs chloroquine and hydroxychloroquine have been suggested as potential anti-COVID-19 therapies in recent studies; therefore, the traditional use of L. zeylanica against malaria may be significant [30] . The ethnopharmacological activities of this plant, including anti-inflammatory, antidiarrheal, antimicrobial, antioxidant, thrombolytic, hepatoprotective, analgesic, larvicidal, and insecticidal activities, have been reported in previous studies [25, 26, 31] . Importantly, the plant exhibits significant antioxidant and thrombolytic properties, which may be useful against these components of the SARS-CoV-2 pathogenesis. This study computationally investigated the roles played by compounds identified in L. zeylanica to combat SARS-CoV-2.

The aerial parts of L. zeylanica were collected from the forest area of Chittagong Hill Tract in November 2015, which was acknowledged by a prominent botanist from Bangladesh Forest Research Institute, Chittagong, Bangladesh. The Bangladesh National Herbarium has stored this plant for future reference with a voucher specimen for identification (accession no. BFRI-107). 

The aerial parts of the L. zeylanica plant were cleaned and cut into tiny pieces (0.4-0.5 mm), dried in air and ground (Moulinex Blender AK-241, Moulinex, France) into a powder (40-80 mesh, 355 g). The powder was immersed in 2 L of methanol for approximately 5 min, and the decoction was allowed to stand for 30 min before being filtered through Whatman No. 1 filter paper. At a temperature of less than 50 • C, filtrate from cheesecloth and Whatman filter paper No. 1 was condensed using a rotating evaporator (RE 200, Bibby Sterling Ltd., Staffordshire, UK). The extracts were placed into glass Petri dishes (90 × 15 mm, Pyrex, Germany) and allowed to air dry to evaporate the solvent completely, resulting in a final yield of 4.4%-5.6% w/w.

The aerial parts of L. zeylanica (methanol extract) were inspected with a mass spectrometer (TQ 8040, Shimadzu Corporation, Kyoto, Japan), using an electron-impact ionization process, combined with a gas chromatograph (GC-17A, Shimadzu Corporation), using a silica capillary (Rxi-5 ms; 0.25 m film, 30 m long and internal diameter 0.32 mm). The oven temperature was set at 70 • C (0 min); 10 • C, 150 • C (5 min); 12 • C, 200 • C (15 min); 12 • C, 220 • C (5 min), with a hold time of 10 min. The inlet temperature was 260 • C. A rate of 0.6 mL/min was used for the flux at a constant pressure of 90 kPa using helium gas. The temperature interface between the GC and MS was maintained at a constant 280 • C. The MS was performed in scan mode, with a range from 40 to 350 amu. The sample was injected at 1 µL, and the entire GC-MS process lasted for 50 min [32] . The results were compared against the National Institute of Standards and Technology (NIST) GC-MS library version 08-S for the identification of compounds in the peak areas.

From the GC-MS analysis, 34 compounds (a total of 35 with the standard) were downloaded in .sdf two-dimensional (2D) format from the PubChem database (https: /pubchem.ncbi.nlm.nih.gov/, accessed on 12 March 2021). For the preparation of the ligand, the LigPrep tool (Maestro v 11.1) was used. The pH 7.0 ± 2.0 was used for the generation of ionization states of the compounds with Epik 2.2 (Force field: OPLS3) in Schrödinger ver.11.1. Up to 32 possible stereoisomers per ligand were retained.

The three-dimensional (3D) structure of the SARS-CoV-2 M pro (PDB ID: 6LU7) [33] was retrieved from the RCSB Protein Data Bank (https://www.rcsb.org/structure/, accessed on 12 March 2021) in PDB format [34] . The Protein Preparation Wizard (Schrödinger ver.11.1) was used to prepare the 6LU7 receptor using the following processes: optimization, removal of water molecules, and minimization (Force field: OPLS3).

The grid generation (Schrödinger ver.11.1) for the selected receptor was performed using the default parameters (Force field: OPLS3). Receptor grids were calculated for the prepared proteins for the observation of poses by various ligands bound within the active predicted site during the docking procedure. The van der Waals radius scaling factor and the partial atomic charge were 1.00 and 0.25, respectively. A cubic box of specific dimensions centered on the centroid of the active site residues was obtained for the receptor. The bounding box was set to 14 × 14 × 14 Å for docking experiments. Ligand docking was followed by the flexible standard precision (Schrödinger ver.11.1), and the docking score and the interactions of the ligand docking were recorded [35] . The results were represented as negative scores in kcal/mol, the final scoring was done on energy-minimized poses and shown as a Glide score. Discovery Studio (DS) version 4.5 was utilized to generate the 2D and 3D representations of the compounds [36] . The figures were generated using Microsoft PowerPoint 2019.

The pharmacokinetic characteristics of all identified phytocompounds were evaluated and screened for drug candidacy using Lipinski's rule of five (RO5) [37] . According to Lipinski's RO5, a compound may exhibit optimal drug-like behavior if the selected parameters are within the specified limit and do not violate more than one of the following five criteria: molecular weight <500 g/mol; ≤5 hydrogen bond donors; ≤10 hydrogen bond acceptors; lipophilicity <5; and molar refractivity between 40 and 130. The web-based tool SwissADME, which is considered to be a convenient drug discovery tool, was used to analyze the druglikeness criteria of the detected biological compounds [38] . Compounds that pass Lipinski's RO5 can be considered suitable candidates for the development of new drug entities.

The dynamic motion of the drug-protein complex was evaluated through a molecular dynamics simulation study. The simulation study was conducted using the YASARA software (version 20.1.1) [39] with the aid of the AMBER14 force field [40] . The N3 inhibitor complex was used as a control system to be compared against the other three complexes. The periodic boundary condition was maintained, and a cubic simulation cell was created that was 20 Å larger than the biological systems in all cases. The NPT ensemble method was used and a Berendsen thermostat was applied to control the temperature of four systems. For the calculation of long-range electrostatic interactions, the particle mesh Ewald method [41] was applied, and the short-range van der Waals and Coulomb interactions were analyzed using a cutoff radius of 8 Å. The TIP3P or transferable intramolecular potential 3 points were used to add Na and Cl ions [42] . The total physiological conditions of the system were set to a temperature of 298 K, pH 7.4, and 0.9% NaCl. For the initial energy minimization process, the steepest descent gradient approach was used with a simulated annealing method. The molecular dynamics simulation trajectories were saved after every 100 ps using a time step of 1.25 fs [43] . The simulation study was conducted for 100 ns to analyze the root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), radius of gyration (Rg), solvent-accessible surface area (SASA), secondary structure, and the number of hydrogen bonds [44] [45] [46] . The molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method was applied to calculate the binding free energy. The YASARA macro file was edited for this calculation [47, 48] . The 1000 trajectory files were considered for MM-PBSA calculation.

A total of 34 compounds were identified from the aerial parts of L. zeylanica using GC-MS, which are listed in Figure 1 and Table 1 , along with their chemical compositions. The total ionic chromatogram (TIC) is shown in Figure 1 . Thirty-four compounds were selected for molecular docking analyses because the specific biological activities of interest for these compounds have not yet been established.

The interactions between various identified compounds from L. zeylanica and the SARS-CoV-2 receptor (PDB ID: 6LU7) are presented in Table 2 . Of the 34 compounds analyzed, 26 compounds interacted with the SARS-CoV-2 receptor, and the three compounds with the highest docking scores were identified as 11-oxa-dispiro Figure S1 . Importantly, the docking experiment delineated that standard inhibitor N3 showed the highest docking score (−7.013 kcal/mol) for the SARS-CoV-2 M pro compared to the other studied compounds. 

A total of 34 compounds were identified from the aerial parts of L. zeylanica using GC-MS, which are listed in Figure 1 and Table 1 , along with their chemical compositions. The total ionic chromatogram (TIC) is shown in Figure 1 . Thirty-four compounds were selected for molecular docking analyses because the specific biological activities of interest for these compounds have not yet been established. 

Hexadecanoic acid, 2-hydroxy- 

Lipinski's rules of five was employed to evaluate the various pharmacokinetic parameters of the identified phytocompounds. Eleven phytocompounds, namely hexadecamethylcyclooctasiloxane, octadecamethylcyclononasiloxane, cyclodecasiloxane, eicosamethyl, N,N -methylenebis(oleamide), (Z,Z); squalene, α-tocopheryl acetate, bis(heptamethylcyclotetrasiloxy)hexamethyltrisiloxane, stigmasterol, γ-sitosterol, 9,19cyclolanost-24-en-3-ol, acetate, (3.beta.), and stigmast-4-en-3-one failed to fulfill Lipinski's RO5, as these compounds contravened more than one rule. However, 23 other compounds conformed with Lipinski's RO5 and may demonstrate optimal drug-like behavior. The result of the absorption, distribution, metabolism, and excretion/toxicity (ADME/T) analysis is shown in Table 3 .

In this simulation study, lorazepam, 11-oxa-dispiro[4.0.4.1]undecan-1-ol, azetidin-2one-3, 3-dimethyl-4-(1-aminoethyl), and inhibitor N3 were denoted as D1, D2, D3, and control, respectively. The RMSD from Figure 6A illustrated that the control drug-protein complex had a higher RMSD trend compared with those of the other three complexes. Initially, none of the four complexes displayed large fluctuations, and they generally remained in a steady state. However, after 40 ns, all complexes had slightly higher and lower RMSD trends, which indicated complex flexibility. Eventually, the complexes returned to a steady trend again for the remainder of the simulation time, exhibiting rigid profiles. The degree of mobility in a biological system can be indicated by the Rg profile. Figure 6B shows that the control drug-protein complex had a lower Rg profile than the experimental drugs, indicating the compacted nature of the protein complex, whereas higher Rg values, which correlate with the repeated folding and unfolding protein behavior, were observed for the protein complexes containing both D1 and D2. By contrast, the complex containing D3 demonstrated initial stability until 40 ns, after which the Rg value increased, which may represent the loose packaging of the system. However, after 60 ns the D3 complex regained a steady Rg pattern, similar to those observed for the other complexes.

The surface area of the biological systems and their corresponding binding patterns with ligand molecules can be assessed through SASA analysis. The D3 complex showed an increasing SASA value until 40 ns, after which the SASA stabilized ( Figure 6C ). This increasing trend in SASA represents protein expansion and comparatively loose binding. By contrast, D1, D2, and the control complex presented with stable SASA profiles and did not deviate. Therefore, these complexes experienced no changes in the surface area and formed more rigid profiles compared with the D3 complex. The quantitative measurement of hydrogen bonds in the drug-protein complex represents the constant nature and molecular recognition of the complexes. The D1 and D3 complexes formed more hydrogen bonds than the control and D2 complexes. More hydrogen bonds indicate an increasingly stable nature for the complex ( Figure 6D ).

The flexibility across the amino acid residues of the protein-drug system can be evaluated through RMSF descriptors. All of the complexes and their respective amino acid residues had lower RMSF values, indicating reduced flexibility. However, some residues, such as Ser1 (helix-strand), Gly2 (helix-strand), Leu50 (beta-turn), Arg60 (helixstrand), Asn72 (beta-turn), Lys97 (beta-turn), Tyr154 (beta-turn), Phe223 (beta-turn), His246 (helix-strand), Ser301 (beta-turn), Gly302 (beta-turn), Val303 (beta-turn), Thr304 (beta-turn), Phe305 (beta-turn), and Gln306 (beta-turn), presented with more flexibility associated with higher RMSF profiles in the molecular dynamics simulation ( Figure 6E ). 

With an increasing number of cases worldwide, the COVID-19 situation con worsen on a daily basis, especially in developed countries such as the USA. Th giousness of SARS-CoV-2 is above and beyond that of SARS-CoV or MERS-C In molecular modeling and computational drug design schemes, the PBSA system is a widely used solvation model for the estimation of the binding free energy of the drugprotein complex. Better binding and more compact interactions are indicated by higher binding energy values in the MM-PBSA calculations [44] [45] [46] . The reference control protein structure of the protein had more binding energy, which indicated a positive interaction pattern. The complex molecules D2 and D3 had similar binding energy patterns to that observed for the control complex, whereas D1 had slightly less binding energy than the control complex, as demonstrated in Figure 6F .

With an increasing number of cases worldwide, the COVID-19 situation continues to worsen on a daily basis, especially in developed countries such as the USA. The contagiousness of SARS-CoV-2 is above and beyond that of SARS-CoV or MERS-CoV, two other members of the beta-coronavirus family. However, no specific treatment has been developed to treat this infectious disease thus far. As a result, the mortality rate continues to increase rapidly. Due to the absence of specific therapeutic drugs, current treatment approaches primarily involve symptom relief and supportive care. Several clinical trials have been conducted for several drug candidates, including remdesivir, hydroxychloroquine, and lopinavir/ritonavir [49] . However, these drugs have not yet amassed sufficient evidence to support their clinical applications. Scientists and researchers worldwide are working together to identify treatment strategies for this deadly coronavirus. Recently, a research study suggested a role for immunopathological considerations in the treatment of SARS-CoV-2 [50] . Wang et al. showed that human monoclonal antibodies could have a neutralizing effect against SARS-CoV-2, primarily by targeting the spike glycoprotein of the virus [51] . However, traditional drug development processes are tedious, and the development of suitable drug candidates can take as long as 15 years, which is not a feasible approach for identifying appropriate cures for COVID-19. Computer-aided drug design may represent a potential method for identifying lead compounds to treat SARS-CoV-2, which can result in both rapid and accurate results. A plethora of studies that have applied computational biology techniques have successfully predicted novel lead compounds for combating SARS-CoV-2 in addition to designing in silico epitope-based vaccines against SARS-CoV-2 [2, 5, 30, 52, 53] . However, these types of studies continue to require further wet-lab verification to be developed into effective therapeutic strategies for COVID-19.

Medicinal plants have been used as an ideal source of therapeutic agents for thousands of years [54, 55] . These therapeutic agents can be developed into medicines and can often be extracted from the crude extracts of several plant parts, including the leaves, stems, bark, rhizomes, fruits, and whole-plant materials [56, 57] . The plant kingdoms feature several biomolecules that possess numerous and varied biological activities [58, 59] . A recent review from Yang et al. reported a role for traditional medicinal practices in the management of SARS-CoV-2 infection [60] . However, due to the lack of comprehensive studies, the effects of biologically active molecules of plant origins against COVID-19 remain relatively unexplored. Through computational studies, researchers can predict the activities of several phytocompounds that might be effective for curbing the activity of the SARS-CoV-2 M pro [30, 61] . The SARS-CoV-2 M pro plays important roles in viral replication and enzymatic cleavage, including the processing of the pp1a and pp1ab polyproteins. We have designed in silico studies to examine the inhibition of SARS-CoV-2 M pro activity using phytocompounds identified from the methanolic extract of the aerial parts of L. zeylanica. The early detection of SARS-CoV-2 infections by the human immune system occurs through pattern recognition receptors (PRRs), which trigger a nuclear factor kappa B (NF-κB) response [62] . NF-κB activation regulates the inflammatory response through the overproduction of proinflammatory cytokines and chemokines. NF-κB has been considered to serve as a potential biomarker for oxidative stress, and the activation of NF-κB can be regulated by antioxidants [63] . Moreover, NF-κB activation leads to tissue factor expression, which ultimately induces a hypercoagulable state. Several previous studies have reported the antioxidant and thrombolytic attributes of ethanolic extracts of L. zeylanica [64, 65] . Therefore, in the present study, we rationalize that several identified phytochemicals derived from the aerial parts of L. zeylanica could potentially play roles against SARS-CoV-2 M pro activation. Molecular docking is the most common computer-aided technique, which has been widely used during the last century [66] . Different algorithms have been developed to increase the accuracy of molecular docking analyses [67] . Molecular docking studies simulate the characteristics of small molecules during binding with the active site of a particular receptor protein. In the current study, we docked several identified phytochemicals identified in the aerial parts of L. zeylanica with the SARS-CoV-2 M pro (PDB ID: 6LU7). In addition, to validate the molecular docking approach, we performed molecular docking analysis using a known molecule inhibitor, N3 (positive control), which was bound to the SARS-CoV-2 M pro crystal structure (PDB ID: 6LU7). Although lorazepam, 2TMS derivative; azetidin-2-one 3, 3-dimethyl-4-(1-aminoethyl); and 11-oxa-dispiro[4.0.4.1]undecan-1-ol possessed the highest docking scores among the targeted compounds, their binding affinities were all weaker than that of the positive control N3. However, unlike inhibitor N3, the lorazepam, 2TMS derivative interacted with the Cys145 residue through both pi-alkyl and pi-sulfur interactions. Moreover, 3-butynoic acid; undecanoic acid, 11-bromo; methyl ester, 6-octadecenoic acid; methyl ester, (Z); hexadecanoic acid-2-hydroxy-1-(hydroxymethyl)ethyl ester; octadecanoic acid-2,3-dihydroxypropyl ester; 13-docosenamide, (Z); α-tocopheryl acetate; campesterol; stigmasterol; γ-sitosterol; 4-campestene-3-one-9; 19-cyclolanost-24-en-3-ol, acetate, (3-beta.); and stigmast-4-en-3-one also interacted with the Cys145 residue, and octadecanoic acid-2, 3-dihydroxypropyl ester; 13-docosenamide, (Z); and stigmasterol all interacted with Cys145 through the formation of a hydrogen bond. Additionally, azetidin-2-one 3, 3-dimethyl-4-(1-aminoethyl) and eight other compounds formed hydrophobic interactions with the His41 residue. The Cys145 and His41 residues have recently been identified as components of the SARS-CoV-2 M pro catalytic dyads. 3-Butynoic acid and octadecanoic acid, 2,3dihydroxypropyl ester each interacted with both the Gly143 and Ser144 residues through the formation of hydrogen bonds. Lorazepam, 2TMS derivative, undecanoic acid, 11bromo-methyl ester, 6-octadecenoic acid, and methyl ester, (Z) also interacted with Gly143, whereas stigmasterol, 9,19-cyclolanost-24-en-3-ol, acetate, (3-beta.), and 4,22-cholestadien-3-one interacted with Ser144 through hydrogen bonding. In addition, 3-butynoic acid, undecanoic acid, 11-bromo-methyl ester, 3,7,11,15-tetramethyl-2-hexadecen-1-ol/phytol, 6-octadecenoic acid, methyl ester, (Z); octadecanoic acid, 2,3-dihydroxypropyl ester, 13docosenamide, (Z); and inhibitor N3 all interacted with the His163 residue, but only octadecanoic acid, 2,3-dihydroxypropyl ester, and 13-docosenamide, (Z) formed hydrogen bonds with the His163 residue, whereas the remaining compounds, including inhibitor N3, yielded only hydrophobic interactions. Azetidin-2-one 3,3-dimethyl-4-(1-aminoethyl), 11-oxa-dispiro[4.0.4.1]undecan-1-ol, 3-azabicyclo[3.2.2]nonane, phytol acetate; phytol, hexadecanoic acid, methyl ester, 9,12-octadecadienoic acid (Z,Z), methyl ester, α-tocopheryl acetate; squalene, campesterol, γ-sitosterol; and 4-campestene-3-one interacted with the Met165 residue through hydrophobic interactions. The analysis of the docked compounds' intermolecular interactions depicted the potentiality of the identified compounds' ability to form interactions with the substrate binding cleft and essential residues in the active pocket of the SARS-CoV-2 M pro . The residues Gly143, Ser144, His163, His164, Met165, Glu166, Leu167, Asp187, Arg188, Gln189, Thr190, Ala191, and Gln192 have been indicated as playing pivotal roles by forming strong hydrogen bonds and substantial hydrophobic bond interactions with the SARS-CoV-2 M pro [68, 69] . Therefore, the predicted interactions of our identified plant compounds delineated as potential proteolytic function inhibitors of the SARS-CoV-2 M pro . Previous work has demonstrated both the antiviral and antioxidative properties of γ-sitosterol [70] . In addition, stigmasterol and its synthetic derivatives have been demonstrated to possess antiviral attributes [71, 72] . Furthermore, a study by Okoro et al. showed that stigmasterol and bauerenol significantly inhibited the HIV-I integrase [73] .

Although drug research has indicated a high level of interest in the investigation of natural products, realistically, the synthesis and purification of vast arrays of new compounds represent a research bottleneck [74] . To compensate for the time-consuming and costly nature of new product development, high-throughput screening techniques have been developed as effective methods for the identification of new hit compounds [75] . However, these hit compounds often result in pharmacokinetic failures that result in elimination following Phase II clinical trials [76] . Therefore, current research works have begun to involve key investigations of pharmacokinetics and pharmacodynamics parameters, which are often referred to as drug-likeness properties [77] . Several rules have been established to facilitate the acceptance of promising molecules, and Lipinski's RO5 is the most commonly applied set of drug-likeness attributes [78] . In the current study, we analyzed the druglikeness properties of the identified compounds based on Lipinski's RO5 by utilizing the SWISS-ADME server. Our analysis indicated that most of the identified phytocompounds derived from the aerial parts of L. zeylanica followed Lipinski's RO5. Furthermore, it was reported that, the antiviral drugs approved by the FDA in the last five years involve a molecular weight of 513.97, hydrogen bond donors of 2.95, and hydrogen bond acceptors of 9.13 [23] . The targeted natural compounds are also within the reported criteria, indicating that these compounds may be evaluated as potential drug molecules.

We used a methodology for computer-aided drug design to detect potent SARS-CoV-2 M pro inhibitors following the L. zeylanica phytochemical analyses. Our research showed that several of the identified phytocompounds from L. zeylanica followed Lipinski's RO5 and were therefore assessed as promising therapeutic molecules. The catalytic residues M pro , Ser 144, Cys 145, and His41 were shown to be linked in postmolecular dynamic structures to the investigated phytochemicals. The molecular dynamics simulations conducted for the docked complexes revealed more insight into the rigidity and binding stability of these protein-drug complexes. Following additional investigations in biological laboratories, these compounds may possibly be developed into effective SARS-CoV-2 therapeutic candidates. 

A comprehensive review on global contributions and recognition of pharmacy professionals amidst COVID-19 pandemic: Moving from present to future. Futur

Epitope-Based Immunoinformatics Approach on Nucleocapsid Protein of Severe Acute Respiratory Syndrome-Coronavirus-2

Outbreak of pneumonia of unknown etiology in Wuhan, China: The mystery and the miracle

COVID-19 in early 2021: Current status and looking forward

Immunoinformatics-guided design of an epitope-based vaccine against severe acute respiratory syndrome coronavirus 2 spike glycoprotein

Remdesivir and chloroquine effectively inhibit the recently emerged novel coronavirus (2019-nCoV) in vitro

Antiviral Drug Discovery Against SARS-CoV

Crystal structure of SARS-CoV-2 main protease provides a basis for design of improved α-ketoamide inhibitors

Mechanisms and enzymes involved in SARS coronavirus genome expression

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

Diverse Immunological Factors Influencing Pathogenesis in Patients with COVID-19: A Review on Viral Dissemination, Immunotherapeutic Options to Counter Cytokine Storm and Inflammatory Responses

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

The origin, transmission and clinical therapies on coronavirus disease 2019 (COVID-19) outbreak-An update on the status

Characterization of anti-viral immunity in recovered individuals infected by SARS-CoV-2

Influenza virus replication in lung epithelial cells depends on redox-sensitive pathways activated by NOX4 -derived ROS

Inhibition of Nox2 Oxidase Activity Ameliorates Influenza A Virus-Induced Lung Inflammation

Oxidative stress favours herpes virus infection in vertebrates: A meta-analysis

Molecular immune pathogenesis and diagnosis of COVID-19

SARS-CoV-2 infection pathogenesis is related to oxidative stress as a response to aggression

Oxidative Stress Associated with SARS-Cov-2 (COVID-19) Increases the Severity of the Lung Disease-A Systematic Review

The Emerging Threat of (Micro)Thrombosis in COVID-19 and Its Therapeutic Implications

Recent efforts for drug identification from phytochemicals against SARS-CoV-2: Exploration of the chemical space to identify druggable leads

A review of plants of genus Leucas

Leucas zeylanica (L.) R. Br. protects ethanol and hydrogen peroxideinduced oxidative stress on hepatic tissue of rats

Mosquito Larvicidal Activities and Chemical Compositions of the Essential Oils of Leucas zeylanica Growing Wild in Vietnam

Plant-based phytochemical screening by targeting main protease of SARS-CoV-2 to design effective potent inhibitors

Antioxidative, antimicrobial and cytotoxic effects of the phenolics of Leea indica leaf extract

A new triterpenoid glucoside from Leucas zeylanica

Biochemical and Computational Approach of Selected Phytocompounds from Tinospora crispa in the Management of COVID-19

Potential of Leucas zeylanica extract to eliminate E. coli and S. aureus in Corbicula fluminea

Virtual screening and molecular dynamics simulation study of plant derived compounds to identify potential inhibitor of main protease from SARS-CoV-2

The Protein Data Bank

Protein and ligand preparation: Parameters, protocols, and influence on virtual screening enrichments

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

Phytochemical screening, antioxidant, thrombolytic, α-amylase inhibition and cytotoxic activities of ethanol extract of Steudnera colocasiifolia K. Koch leaves

Making optimal use of empirical energy functions: Force-field parameterization in crystal space

The Amber Lipid Force Field

Fast empirical pKa prediction by Ewald summation

YASARA-yet another scientific artificial reality application. YASARA Org

New ways to boost molecular dynamics simulations

Exploring the potent inhibitors and binding modes of phospholipase A2 through in silico investigation

Comparative molecular investigation of the potential inhibitors against SARS-CoV-2 main protease: A molecular docking study

Assessment of structurally and functionally high-risk nsSNPs impacts on human bone morphogenetic protein receptor type IA (BMPR1A) by computational approach

In vivo neuroprotective, antinociceptive, anti-inflammatory potential in Swiss albino mice and in vitro antioxidant and clot lysis activities of fractionated Holigarna longifolia Roxb. bark extract

Identification of COX-2 inhibitors via structure-based virtual screening and molecular dynamics simulation

Potential covalent drugs targeting the main protease of the SARS-CoV-2 coronavirus

COVID-19: Immunopathology and its implications for therapy

Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: A descriptive study

A Sequence Homology and Bioinformatic Approach Can Predict Candidate Targets for Immune Responses to SARS-CoV-2

Appraisals of the Bangladeshi Medicinal Plant Calotropis gigantea Used by Folk Medicine Practitioners in the Management of COVID-19: A

GC-MS Phytochemical Profiling, Pharmacological Properties, and In Silico Studies of Chukrasia velutina Leaves: A Novel Source for Bioactive Agents

Anxiolytic, Antinociceptive, Anti-inflammatory and Antipyretic Effects of a Chloroform Extract from the Leaves of Urena sinuata in Rodents

In vivo and in vitro pharmacological activities of Tacca integrifolia rhizome and investigation of possible lead compounds against breast cancer through in silico approaches

Immunoinformatics-guided design of multi-epitope vaccine from structural proteins of severe acute respiratory syndrome-coronavirus-2

Pharmacological studies on the antinociceptive, anxiolytic and antidepressant activity of Tinospora crispa

Bin Emran, T. Antipyretic and hepatoprotective potential of Tinospora crispa and investigation of possible lead compounds through in silico approaches

Traditional Chinese Medicine in the Treatment of Patients Infected with 2019-New Coronavirus (SARS-CoV-2): A Review and Perspective

A molecular modelling approach for identifying antiviral selenium-containing heterocyclic compounds that inhibit the main protease of SARS-CoV-2: An in silico investigation

Coronavirus COV-19/SARS-CoV-2 affects women less than men: Clinical response to viral infection

Transcription factor NF-κB as a potential biomarker for oxidative stress

Phytochemical and Pharmacological Activities of Leucas zeylanica (L.) R. Br. (Family: Lamiaceae)

Drug Repurposing Approach against Novel Coronavirus Disease (COVID-19) through Virtual Screening Targeting SARS-CoV-2 Main Protease

A geometric approach to macromolecule-ligand interactions

Molecular Docking: A Powerful Approach for Structure-Based Drug Discovery

SARS-CoV 3CL protease cleaves its C-terminal autoprocessing site by novel subsite cooperativity

Mechanism of the Maturation Process of SARS-CoV 3CL Protease

Potential of plant bioactive compounds as SARS-CoV-2 main protease (Mpro) and spike (S) glycoprotein inhibitors: A molecular docking study

Antiviral Action of Synthetic Stigmasterol Derivatives on Herpes Simplex Virus Replication in Nervous CellsIn Vitro

Antiviral activity of 7-keto-stigmasterol obtained from green Antarctic algae Prasiola crispa against equine herpesvirus 1. Environ. Boil

Potential of leaf extract of Pangium edule Reinw as HIV-1 protease inhibitor: A computational biology approach

Efficacy of Phytochemicals Derived from Avicennia officinalis for the Management of COVID-19: A Combined In Silico and Biochemical Study

A brief history of novel drug discovery technologies

A practical view of 'druggability'

Assessing drug-likeness-What are we missing?

Identification of potential phytochemicals from Citrus limon against main protease of SARS-CoV-2: Molecular docking, molecular dynamic simulations and quantum computations