Next Article in Journal
Cigarette Smoking Alters the Expression of Circulating microRNAs and Its Potential Diagnostic Value in Female Lung Cancer Patients
Next Article in Special Issue
The Role of Acidosis in the Pathogenesis of Severe Forms of COVID-19
Previous Article in Journal
A Comparison of Acute Pharmacological Effects of Methylone and MDMA Administration in Humans and Oral Fluid Concentrations as Biomarkers of Exposure
Previous Article in Special Issue
Alcohol, Cannabis and Crossfading: Concerns for COVID-19 Disease Severity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Phytochemicals from Leucas zeylanica Targeting Main Protease of SARS-CoV-2: Chemical Profiles, Molecular Docking, and Molecular Dynamics Simulations

1
Department of Pharmacy, BGC Trust University Bangladesh, Chittagong 4381, Bangladesh
2
Department of Pharmacy, Jahangirnagar University, Savar, Dhaka 1342, Bangladesh
3
Department of Pharmacy, International Islamic University Chittagong, Chittagong 4318, Bangladesh
4
Department of Pharmacy, Faculty of Biological Sciences, University of Chittagong, Chittagong 4331, Bangladesh
5
Microbiology Laboratory, Department of Genetic Engineering and Biotechnology, University of Rajshahi, Rajshahi 6205, Bangladesh
6
Drug Discovery, GUSTO A Research Group, Chittagong 4203, Bangladesh
7
Department of Physics, College of Khurma, Taif University, Taif 21944, Saudi Arabia
8
Department of Chemistry, Faculty of Science, Taif University, Taif 21974, Saudi Arabia
9
Centre for Applied Physics and Radiation Technologies, School of Engineering and Technology, Sunway University, Bandar Sunway 47500, Malaysia
10
Department of Physics, University of Surrey, Guilford GU2 7XH, UK
*
Authors to whom correspondence should be addressed.
Biology 2021, 10(8), 789; https://doi.org/10.3390/biology10080789
Submission received: 14 June 2021 / Revised: 5 August 2021 / Accepted: 15 August 2021 / Published: 17 August 2021
(This article belongs to the Special Issue Coronavirus Disease 2019 (COVID-19))

Abstract

:

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 (Mpro) 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 Mpro. 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 Mpro, 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 Mpro 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–Mpro 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 Mpro binding pockets. The results of our computerized drug design approach may assist in the fight against SARS-CoV-2.

1. Introduction

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 beta-coronavirus 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 (Mpro) 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 Mpro (also known as 3CLPro 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 Mpro 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 Mpro. Additionally, Mpro 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 Mpro enzyme-specific T cells have been encountered in SARS-CoV-2 patients [13,14]. Therefore, Mpro 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 mechanisms [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.

2. Materials and Methods

2.1. Collection and Identification of Plant Material

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).

2.2. Preparation of Plant Extract and Decoction Preparation

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.

2.3. Gas Chromatography–Mass Spectrometry (GC-MS) Analysis

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.

2.4. Molecular Docking Study

2.4.1. Ligand Preparation

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.

2.4.2. Protein Preparation

The three-dimensional (3D) structure of the SARS-CoV-2 Mpro (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).

2.4.3. Receptor Grid Generation and Glide Molecular Docking

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.

2.5. ADME/T Properties Analysis

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 drug-likeness 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.

2.6. Molecular Dynamics Simulation

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.

3. Results

3.1. GC-MS Analysis

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.

3.2. Molecular Docking Study

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[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), as represented in Figure 2, Figure 3, Figure 4 and Figure 5 and Supplementary 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 Mpro compared to the other studied compounds.

3.3. ADME/T Properties Analysis

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,19-cyclolanost-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.

3.4. Molecular Dynamics Simulation

In this simulation study, lorazepam, 11-oxa-dispiro[4.0.4.1]undecan-1-ol, azetidin-2-one-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 (helix-strand), 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).
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 drug–protein 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.

4. Discussion

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 Mpro [30,61]. The SARS-CoV-2 Mpro 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 Mpro 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 Mpro 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 Mpro (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 Mpro 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 Mpro catalytic dyads. 3-Butynoic acid and octadecanoic acid, 2,3-dihydroxypropyl ester each interacted with both the Gly143 and Ser144 residues through the formation of hydrogen bonds. Lorazepam, 2TMS derivative, undecanoic acid, 11-bromo- 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, 13-docosenamide, (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 Mpro. 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 Mpro [68,69]. Therefore, the predicted interactions of our identified plant compounds delineated as potential proteolytic function inhibitors of the SARS-CoV-2 Mpro. 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 drug-likeness 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.

5. Conclusions

We used a methodology for computer-aided drug design to detect potent SARS-CoV-2 Mpro 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 Mpro, 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.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/biology10080789/s1, Supplementary File S1: A 2D representation of 23 compounds against the SARS-CoV-2 receptor (PDB: 6LU7). (a) 3-Butynoic acid; (b) cyanoacetic acid; (c) methyl 11-bromoundecanoate; (d) 3-azabicyclo[3.2.2]nonane; (e) phytol acetate; (f) 3,7,11,15-tetramethyl-2-hexadecen-1-ol; (g) hexadecanoic acid, methyl ester; (h) 9,12-octadecadienoic acid (Z,Z)-, methyl ester; (i) 6-octadecenoic acid, methyl ester, (Z)-; (j) methyl stearate; (k) hexadecanoic acid, 2-hydroxy-1-(hydroxymethyl)ethyl ester; (l) octadecanoic acid, 2,3-dihydroxypropyl ester; (m) 13-docosenamide, (Z)-; (n) N,N’-methylenebis(oleamide), (Z, Z)-; (o) squalene; (p) α-tocopheryl acetate; (q) campesterol; (r) stigmasterol; (s) γ-sitosterol; (t) 4-campestene-3-one; (u) 9,19-cyclolanost-24-en-3-ol, acetate, (3.beta.); (v) 4,22-cholestadien-3-one; and (w) stigmast-4-en-3-one.

Author Contributions

Conceptualization: M.D., A.M.T., A.R., S.M. and T.B.E.; methodology: M.D., A.M.T., A.R., S.M., S.A.S., J.M., M.N.I., M.M., M.Z.U., M.S.R. and T.B.E.; software: M.D., A.M.T., A.R., S.M., A.S.A.A. and T.B.E.; validation: A.A., A.S.A.A., M.U.K., D.A.B., M.S.R. and T.B.E.; formal analysis: A.A., A.S.A.A., M.U.K., D.A.B., M.S.R. and T.B.E.; investigation: M.D., A.M.T., A.R., S.M., S.A.S., J.M., M.N.I., M.M., M.Z.U., M.S.R. and T.B.E.; resources: M.S.R. and T.B.E.; data curation: M.D., A.M.T., A.R., S.M. and T.B.E.; writing—original draft preparation: M.D., A.M.T., A.R., S.M. and T.B.E.; writing—review and editing: M.N.I., A.A., A.S.A.A., M.U.K., D.A.B., M.S.R. and T.B.E.; visualization: K M.D., A.M.T., A.R., S.M. and T.B.E.; supervision: M.S.R. and T.B.E.; project administration: M.S.R. and T.B.E.; funding acquisition: A.A., A.S.A.A., M.U.K., D.A.B., M.S.R. and T.B.E. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Available data are presented in the manuscript.

Acknowledgments

This research was supported by the Taif University Researchers Supporting Project number (TURSP-2020/163), Taif 21944, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sami, S.A.; Marma, K.K.S.; Chakraborty, A.; Singha, T.; Rakib, A.; Uddin, G.; Hossain, M.K.; Uddin, S.M.N. A comprehensive review on global contributions and recognition of pharmacy professionals amidst COVID-19 pandemic: Moving from present to future. Futur. J. Pharm. Sci. 2021, 7, 1–16. [Google Scholar] [CrossRef]
  2. Rakib, A.; Sami, S.; Islam, A.; Ahmed, S.; Faiz, F.; Khanam, B.; Marma, K.; Rahman, M.; Uddin, M.; Nainu, F.; et al. Epitope-Based Immunoinformatics Approach on Nucleocapsid Protein of Severe Acute Respiratory Syndrome-Coronavirus-2. Molecules 2020, 25, 5088. [Google Scholar] [CrossRef]
  3. Lu, H.; Stratton, C.W.; Tang, Y. Outbreak of pneumonia of unknown etiology in Wuhan, China: The mystery and the miracle. J. Med. Virol. 2020, 92, 401–402. [Google Scholar] [CrossRef] [Green Version]
  4. Wang, C.; Wang, Z.; Wang, G.; Lau, J.Y.-N.; Zhang, K.; Li, W. COVID-19 in early 2021: Current status and looking forward. Signal Transduct. Target. Ther. 2021, 6, 1–14. [Google Scholar] [CrossRef]
  5. Rakib, A.; Sami, S.A.; Mimi, N.J.; Chowdhury, M.; Eva, T.A.; Nainu, F.; Paul, A.; Shahriar, A.; Tareq, A.M.; Emon, N.U.; et al. Immunoinformatics-guided design of an epitope-based vaccine against severe acute respiratory syndrome coronavirus 2 spike glycoprotein. Comput. Biol. Med. 2020, 124, 103967. [Google Scholar] [CrossRef]
  6. Wang, M.; Cao, R.; Zhang, L.; Yang, X.; Liu, J.; Xu, M.; Shi, Z.; Hu, Z.; Zhong, W.; Xiao, G. Remdesivir and chloroquine effectively inhibit the recently emerged novel coronavirus (2019-nCoV) in vitro. Cell Res. 2020, 30, 269–271. [Google Scholar] [CrossRef]
  7. Wu, Y.-S.; Lin, W.-H.; Hsu, J.T.-A.; Hsieh, H.-P. Antiviral Drug Discovery Against SARS-CoV. Curr. Med. Chem. 2006, 13, 2003–2020. [Google Scholar] [CrossRef] [PubMed]
  8. Zhang, L.; Lin, D.; Sun, X.; Curth, U.; Drosten, C.; Sauerhering, L.; Becker, S.; Rox, K.; Hilgenfeld, R. Crystal structure of SARS-CoV-2 main protease provides a basis for design of improved α-ketoamide inhibitors. Science 2020, 368, 409–412. [Google Scholar] [CrossRef] [Green Version]
  9. Thiel, V.; Ivanov, K.A.; Putics, A.; Hertzig, T.; Schelle, B.; Bayer, S.; Weißbrich, B.; Snijder, E.; Rabenau, H.; Doerr, H.W.; et al. Mechanisms and enzymes involved in SARS coronavirus genome expression. J. Gen. Virol. 2003, 84, 2305–2315. [Google Scholar] [CrossRef]
  10. Hilgenfeld, R. From SARS to MERS: Crystallographic studies on coronaviral proteases enable antiviral drug design. FEBS J. 2014, 281, 4085–4096. [Google Scholar] [CrossRef] [Green Version]
  11. Rabaan, A.A.; Al-Ahmed, S.H.; Garout, M.A.; Al-Qaaneh, A.M.; Sule, A.A.; Tirupathi, R.; Mutair, A.A.; Alhumaid, S.; Al-Omari, A.; Hasan, A.; et al. Diverse Immunological Factors Influencing Pathogenesis in Patients with COVID-19: A Review on Viral Dissemination, Immunotherapeutic Options to Counter Cytokine Storm and Inflammatory Responses. Pathogens 2021, 10, 565. [Google Scholar] [CrossRef]
  12. Anand, K. Ziebuhr J Fau-Wadhwani P.; Wadhwani P Fau-Mesters JR; Mesters Jr Fau-Hilgenfeld R.; Hilgenfeld R. Coronavirus main proteinase (3CLpro) structure: Basis for design of anti-SARS drugs. Science 2003, 300, 1763–1767. [Google Scholar] [CrossRef] [Green Version]
  13. Guo, Y.-R.; Cao, Q.-D.; Hong, Z.-S.; Tan, Y.-Y.; Chen, S.-D.; Jin, H.-J.; Tan, K.-S.; Wang, D.-Y.; Yan, Y. The origin, transmission and clinical therapies on coronavirus disease 2019 (COVID-19) outbreak—An update on the status. Mil. Med. Res. 2020, 7, 1–10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Dong, C.; Ni, L.; Ye, F.; Chen, M.-L.; Feng, Y.; Deng, Y.-Q.; Zhao, H.; Wei, P.; Ge, J.; Li, X. Characterization of anti-viral immunity in recovered individuals infected by SARS-CoV-2. MedRxiv 2020. [Google Scholar] [CrossRef] [Green Version]
  15. Amatore, D.; Sgarbanti, R.; Aquilano, K.; Baldelli, S.; Limongi, D.; Civitelli, L.; Nencioni, L.; Garaci, E.; Ciriolo, M.R.; Palamara, A.T. Influenza virus replication in lung epithelial cells depends on redox-sensitive pathways activated by NOX4 -derived ROS. Cell. Microbiol. 2015, 17, 131–145. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Vlahos, R.; Stambas, J.; Bozinovski, S.; Broughton, B.R.S.; Drummond, G.; Selemidis, S. Inhibition of Nox2 Oxidase Activity Ameliorates Influenza A Virus-Induced Lung Inflammation. PLoS Pathog. 2011, 7, e1001271. [Google Scholar] [CrossRef] [Green Version]
  17. Sebastiano, M.; Chastel, O.; de Thoisy, B.; Eens, M.; Costantini, D. Oxidative stress favours herpes virus infection in vertebrates: A meta-analysis. Curr. Zoöl. 2016, 62, 325–332. [Google Scholar] [CrossRef] [Green Version]
  18. Khomich, O.A.; Kochetkov, S.N.; Bartosch, B.; Ivanov, A.V. Redox Biology of Respiratory Viral Infections. Viruses 2018, 10, 392. [Google Scholar] [CrossRef] [Green Version]
  19. Li, X.; Geng, M.; Peng, Y.; Meng, L.; Lu, S. Molecular immune pathogenesis and diagnosis of COVID-19. J. Pharm. Anal. 2020, 10, 102–108. [Google Scholar] [CrossRef]
  20. Cecchini, R.; Cecchini, A.L. SARS-CoV-2 infection pathogenesis is related to oxidative stress as a response to aggression. Med Hypotheses 2020, 143, 110102. [Google Scholar] [CrossRef]
  21. Samir, D. Oxidative Stress Associated with SARS-Cov-2 (COVID-19) Increases the Severity of the Lung Disease–A Systematic Review. J. Infect. Dis. Epidemiology 2020, 6, 121. [Google Scholar] [CrossRef]
  22. McFadyen, J.D.; Stevens, H.; Peter, K. The Emerging Threat of (Micro)Thrombosis in COVID-19 and Its Therapeutic Implications. Circ. Res. 2020, 127, 571–587. [Google Scholar] [CrossRef] [PubMed]
  23. Joshi, G.; Sindhu, J.; Thakur, S.; Rana, A.; Sharma, G.; Mayank; Poduri, R. Recent efforts for drug identification from phytochemicals against SARS-CoV-2: Exploration of the chemical space to identify druggable leads. Food Chem. Toxicol. 2021, 152, 112160. [Google Scholar] [CrossRef]
  24. Chouhan, H.S.; Singh, S.K. A review of plants of genus Leucas. J. Pharmacogn. Phytother. 2011, 3, 13–26. [Google Scholar] [CrossRef]
  25. Hossain, S.; Rahman, M.; Fatima, N.; Haque, M.; Islam, J. Leucas zeylanica (L.) R. Br. protects ethanol and hydrogen peroxide-induced oxidative stress on hepatic tissue of rats. Int. Curr. Pharm. J. 2013, 2, 148–151. [Google Scholar] [CrossRef] [Green Version]
  26. Hung, N.H.; Chuong, N.T.H.; Satyal, P.; Hieu, H.V.; Dai, D.N.; Huong, L.T.; Sinh, L.H.; Ngoc, N.T.B.; Hien, V.T.; Setzer, W.N. Mosquito Larvicidal Activities and Chemical Compositions of the Essential Oils of Leucas zeylanica Growing Wild in Vietnam. Nat. Prod. Commun. 2019, 14, 1–7. [Google Scholar] [CrossRef]
  27. Mahmud, S.; Biswas, S.; Paul, G.K.; Mita, M.A.; Promi, M.M.; Afrose, S.; Hasan, R.; Zaman, S.; Uddin, M.S.; Dhama, K.; et al. Plant-based phytochemical screening by targeting main protease of SARS-CoV-2 to design effective potent inhibitors. Biology 2021, 10, 589. [Google Scholar] [CrossRef] [PubMed]
  28. Rahman, M.A.; bin Imran, T.; Islam, S. Antioxidative, antimicrobial and cytotoxic effects of the phenolics of Leea indica leaf extract. Saudi J. Biol. Sci. 2013, 20, 213–225. [Google Scholar] [CrossRef] [Green Version]
  29. Chen, G.-Y.; Zhang, B.; Zhao, T.; Nidhal, N.; Jia-Li, W.; Zhou, X.-M.; Chun-Yan, D. A new triterpenoid glucoside from Leucas zeylanica. Nat. Prod. Res. 2020, 34, 1874–1878. [Google Scholar] [CrossRef]
  30. Rakib, A.; Paul, A.; Chy, N.U.; Sami, S.A.; Baral, S.K.; Majumder, M.; Tareq, A.M.; Amin, M.N.; Shahriar, A.; Uddin, Z.; et al. Biochemical and Computational Approach of Selected Phytocompounds from Tinospora crispa in the Management of COVID-19. Molecules 2020, 25, 3936. [Google Scholar] [CrossRef]
  31. Abdullah, F.; Nasir, S.N.A.M.; Han, D.K.; Appalasamy, S.; Nor, M.M.; Rak, A.E. Potential of Leucas zeylanica extract to eliminate E. coli and S. aureus in Corbicula fluminea (“Etak”) tissue. Malays. J. Fundam. Appl. Sci. 2019, 15, 597–599. [Google Scholar] [CrossRef]
  32. Uddin, Z.; Paul, A.; Rakib, A.; Sami, S.; Mahmud, S.; Rana, S.; Hossain, S.; Tareq, A.; Dutta, M.; Emran, T.; et al. Chemical Profiles and Pharmacological Properties with In Silico Studies on Elatostema papillosum Wedd. Molecules 2021, 26, 809. [Google Scholar] [CrossRef]
  33. Mahmud, S.; Uddin, M.A.R.; Paul, G.K.; Shimu, M.S.S.; Islam, S.; Rahman, E.; Islam, A.; Islam, M.S.; Promi, M.M.; Emran, T.B.; et al. Virtual screening and molecular dynamics simulation study of plant derived compounds to identify potential inhibitor of main protease from SARS-CoV-2. Brief. Bioinform. 2021, 22, 1402–1414. [Google Scholar] [CrossRef]
  34. Berman, H.M.; Battistuz, T.; Bhat, T.N.; Bluhm, W.F.; Bourne, P.E.; Burkhardt, K.; Feng, Z.; Gilliland, G.L.; Iype, L.; Jain, S.; et al. The Protein Data Bank. Acta Crystallogr. Sect. D Biol. Crystallogr. 2002, 58, 899–907. [Google Scholar] [CrossRef]
  35. Sastry, G.M.; Adzhigirey, M.; Day, T.; Annabhimoju, R.; Sherman, W. Protein and ligand preparation: Parameters, protocols, and influence on virtual screening enrichments. J. Comput. Mol. Des. 2013, 27, 221–234. [Google Scholar] [CrossRef]
  36. Studio, D. Dassault Systemes BIOVIA, Discovery Studio Modelling Environment, Release 4.5; Accelrys Softw. Inc.: San Diego, CA, USA, 2015. [Google Scholar]
  37. Lipinski, C.A.; Lombardo, F.; Dominy, B.W.; Feeney, P.J. Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv. Drug Deliv. Rev. 1997, 23, 3–25. [Google Scholar] [CrossRef]
  38. Kabir, M.S.H.; Hossain, M.M.; Kabir, M.I.; Rahman, M.M.; Hasanat, A.; Emran, T.B.; Rahman, M.A. Phytochemical screening, antioxidant, thrombolytic, α-amylase inhibition and cytotoxic activities of ethanol extract of Steudnera colocasiifolia K. Koch leaves. J. Young Pharmacists 2016, 8, 391–397. [Google Scholar] [CrossRef] [Green Version]
  39. Krieger, E.; Darden, T.; Nabuurs, S.B.; Finkelstein, A.; Vriend, G. Making optimal use of empirical energy functions: Force-field parameterization in crystal space. Proteins: Struct. Funct. Bioinform. 2004, 57, 678–683. [Google Scholar] [CrossRef] [PubMed]
  40. Dickson, C.J.; Madej, B.D.; Skjevik, Å.A.; Betz, R.M.; Teigen, K.; Gould, I.R.; Walker, R.C. Lipid14: The Amber Lipid Force Field. J. Chem. Theory Comput. 2014, 10, 865–879. [Google Scholar] [CrossRef] [PubMed]
  41. Krieger, E.; Nielsen, J.E.; Spronk, C.A.; Vriend, G. Fast empirical pKa prediction by Ewald summation. J. Mol. Graph. Model. 2006, 25, 481–486. [Google Scholar] [CrossRef]
  42. Krieger, E.; Vriend, G.; Spronk, C. YASARA–yet another scientific artificial reality application. YASARA Org. 2013, 993, 51–78. [Google Scholar]
  43. Krieger, E.; Vriend, G. New ways to boost molecular dynamics simulations. J. Comput. Chem. 2015, 36, 996–1007. [Google Scholar] [CrossRef]
  44. Mahmud, S.; Parves, R.; Riza, Y.; Sujon, K.M.; Ray, S.; Alam Tithi, F.; Zaoti, Z.F.; Alam, S.; Absar, N. Exploring the potent inhibitors and binding modes of phospholipase A2 through in silico investigation. J. Biomol. Struct. Dyn. 2020, 38, 4221–4231. [Google Scholar] [CrossRef]
  45. Khan, A.; Mahmud, S.; Alam, A.S.M.R.U.; Rahman, E.; Ahmed, F.; Rahmatullah, M. Comparative molecular investigation of the potential inhibitors against SARS-CoV-2 main protease: A molecular docking study. J. Biomol. Struct. Dyn. 2020, 1–7. [Google Scholar] [CrossRef]
  46. Islam, J.; Parves, R.; Mahmud, S.; Alam Tithi, F.; Reza, A. Assessment of structurally and functionally high-risk nsSNPs impacts on human bone morphogenetic protein receptor type IA (BMPR1A) by computational approach. Comput. Biol. Chem. 2019, 80, 31–45. [Google Scholar] [CrossRef]
  47. Uddin, M.Z.; Rana, M.S.; Hossain, S.; Dutta, E.; Ferdous, S.; Dutta, M.; Emran, T.B. 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. J. Complement. Integr. Med. 2019, 17, 1–10. [Google Scholar] [CrossRef]
  48. Razzaghi-Asl, N.; Mirzayi, S.; Mahnam, K.; Sepehri, S. Identification of COX-2 inhibitors via structure-based virtual screening and molecular dynamics simulation. J. Mol. Graph. Model. 2018, 83, 138–152. [Google Scholar] [CrossRef] [PubMed]
  49. Liu, S.; Zheng, Q.; Wang, Z. Potential covalent drugs targeting the main protease of the SARS-CoV-2 coronavirus. Bioinformatics 2020, 36, 3295–3298. [Google Scholar] [CrossRef] [PubMed]
  50. Cao, X. COVID-19: Immunopathology and its implications for therapy. Nat. Rev. Immunol. 2020, 20, 269–270. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Chen, N.; Zhou, M.; Dong, X.; Qu, J.; Gong, F.; Han, Y.; Qiu, Y.; Wang, J.; Liu, Y.; Wei, Y. Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: A descriptive study. Lancet 2020, 395, 507–513. [Google Scholar] [CrossRef] [Green Version]
  52. Grifoni, A.; Sidney, J.; Zhang, Y.; Scheuermann, R.H.; Peters, B.; Sette, A. A Sequence Homology and Bioinformatic Approach Can Predict Candidate Targets for Immune Responses to SARS-CoV-2. Cell Host Microbe 2020, 27, 671–680. [Google Scholar] [CrossRef]
  53. Dutta, M.; Nezam, M.; Chowdhury, S.; Rakib, A.; Paul, A.; Sami, S.A.; Uddin, M.Z.; Rana, M.S.; Hossain, S.; Effendi, Y.; et al. Appraisals of the Bangladeshi Medicinal Plant Calotropis gigantea Used by Folk Medicine Practitioners in the Management of COVID-19: A Biochemical and Computational Approach. Front. Mol. Biosci. 2021, 8, 625391. [Google Scholar] [CrossRef]
  54. Jahan, I.; Tona, M.R.; Sharmin, S.; Sayeed, M.A.; Tania, F.Z.; Paul, A.; Chy, N.U.; Rakib, A.; Bin Emran, T.; Simal-Gandara, J. GC-MS Phytochemical Profiling, Pharmacological Properties, and In Silico Studies of Chukrasia velutina Leaves: A Novel Source for Bioactive Agents. Molecules 2020, 25, 3536. [Google Scholar] [CrossRef]
  55. Bin Emran, T.; Mowla, T.-E.; Ahmed, S.; Zahan, S.; Rakib, A.; Hasan, M.S.; Amin, M.N.; Mow, T.R.; Uddin, M.M.N. Sedative, Anxiolytic, Antinociceptive, Anti-inflammatory and Antipyretic Effects of a Chloroform Extract from the Leaves of Urena sinuata in Rodents. J. Appl. Life Sci. Int. 2018, 16, 1–19. [Google Scholar] [CrossRef]
  56. Ahmed, S.; Rakib, A.; Islam, A.; Khanam, B.H.; Faiz, F.B.; Paul, A.; Chy, N.U.; Alam Bhuiya, N.M.M.; Uddin, M.M.N.; Ullah, S.M.A.; et al. In vivo and in vitro pharmacological activities of Tacca integrifolia rhizome and investigation of possible lead compounds against breast cancer through in silico approaches. Clin. Phytosci. 2019, 5, 1–13. [Google Scholar] [CrossRef]
  57. Obaidullah, A.J.; Alanazi, M.A.; Alsaif, N.A.; Albassam, H.; Almehizia, A.A.; Alqahtani, A.A.; Mahmud, S.; Sami, S.A.; Emran, T.B. Immunoinformatics-guided design of multi-epitope vaccine from structural proteins of severe acute respiratory syndrome-coronavirus-2. RSC Adv. 2021, 11, 18103–18121. [Google Scholar] [CrossRef]
  58. Rakib, A.; Ahmed, S.; Islam, A.; Uddin, M.M.N.; Paul, A.; Chy, N.U.; Bin Emran, T.; Seidel, V. Pharmacological studies on the antinociceptive, anxiolytic and antidepressant activity of Tinospora crispa. Phytother. Res. 2020, 34, 2978–2984. [Google Scholar] [CrossRef]
  59. Rakib, A.; Ahmed, S.; Islam, A.; Haye, A.; Uddin, S.M.N.; Uddin, M.M.N.; Hossain, M.K.; Paul, A.; Bin Emran, T. Antipyretic and hepatoprotective potential of Tinospora crispa and investigation of possible lead compounds through in silico approaches. Food Sci. Nutr. 2020, 8, 547–556. [Google Scholar] [CrossRef] [Green Version]
  60. Yang, Y.; Islam, S.; Wang, J.; Li, Y.; Chen, X. Traditional Chinese Medicine in the Treatment of Patients Infected with 2019-New Coronavirus (SARS-CoV-2): A Review and Perspective. Int. J. Biol. Sci. 2020, 16, 1708–1717. [Google Scholar] [CrossRef]
  61. Rakib, A.; Nain, Z.; Islam, M.A.; Sami, S.A.; Mahmud, S.; Islam, A.; Ahmed, S.; Siddiqui, A.B.F.; Babu, S.M.O.F.; Hossain, P.; et al. A molecular modelling approach for identifying antiviral selenium-containing heterocyclic compounds that inhibit the main protease of SARS-CoV-2: An in silico investigation. Brief. Bioinform. 2021, 22, 1476–1498. [Google Scholar] [CrossRef] [PubMed]
  62. Conti, P.; Younes, A. Coronavirus COV-19/SARS-CoV-2 affects women less than men: Clinical response to viral infection. J. Biol. Regul. Homeost. Agents 2020, 34, 339–343. [Google Scholar]
  63. Berg, R.V.D.; Haenen, G.; Bast, A. Transcription factor NF-κB as a potential biomarker for oxidative stress. Br. J. Nutr. 2001, 86, S121–S127. [Google Scholar] [CrossRef] [Green Version]
  64. Islam, M.A.; Biswas, N.N. Phytochemical and Pharmacological Activities of Leucas zeylanica (L.) R. Br. (Family: Lamiaceae). Hamdard Medicus 2010, 53, 19–29. [Google Scholar]
  65. Chowdhury, K.H.; Chowdhury, M.R.; Mahmud, S.; Tareq, A.M.; Hanif, N.B.; Banu, N.; Reza, A.S.M.A.; Emran, T.B.; Simal-Gandara, J. Drug Repurposing Approach against Novel Coronavirus Disease (COVID-19) through Virtual Screening Targeting SARS-CoV-2 Main Protease. Biology 2020, 10, 2. [Google Scholar] [CrossRef]
  66. Kuntz, I.D.; Blaney, J.M.; Oatley, S.J.; Langridge, R.; Ferrin, T.E. A geometric approach to macromolecule-ligand interactions. J. Mol. Biol. 1982, 161, 269–288. [Google Scholar] [CrossRef]
  67. Meng, X.-Y.; Zhang, H.-X.; Mezei, M.; Cui, M. Molecular Docking: A Powerful Approach for Structure-Based Drug Discovery. Curr. Comput. Drug Des. 2011, 7, 146–157. [Google Scholar] [CrossRef]
  68. Muramatsu, T.; Takemoto, C.; Kim, Y.-T.; Wang, H.; Nishii, W.; Terada, T.; Shirouzu, M.; Yokoyama, S. SARS-CoV 3CL protease cleaves its C-terminal autoprocessing site by novel subsite cooperativity. Proc. Natl. Acad. Sci. USA 2016, 113, 12997–13002. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Hsu, M.-F.; Kuo, C.-J.; Chang, K.-T.; Chang, H.-C.; Chou, C.-C.; Ko, T.-P.; Shr, H.-L.; Chang, G.-G.; Wang, A.H.-J.; Liang, P.-H. Mechanism of the Maturation Process of SARS-CoV 3CL Protease. J. Biol. Chem. 2005, 280, 31257–31266. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Tallei, T.E.; Tumilaar, S.G.; Niode, N.J.; Fatimawali; Kepel, B.J.; Idroes, R.; Effendi, Y.; Sakib, S.A.; Emran, T.B. Potential of plant bioactive compounds as SARS-CoV-2 main protease (Mpro) and spike (S) glycoprotein inhibitors: A molecular docking study. Scientifica 2020, 2020, 6307457. [Google Scholar] [CrossRef] [PubMed]
  71. Petrera, E.; Níttolo, A.G.; Alché, L.E. Antiviral Action of Synthetic Stigmasterol Derivatives on Herpes Simplex Virus Replication in Nervous CellsIn Vitro. BioMed Res. Int. 2014, 2014, 1–9. [Google Scholar] [CrossRef] [Green Version]
  72. Marinho, R.D.S.S.; Ramos, C.J.B.; Leite, J.P.G.; Teixeira, V.L.; Paixão, I.C.N.D.P.; Belo, C.A.D.; Pereira, A.B.; Pinto, A.M.V. Antiviral activity of 7-keto-stigmasterol obtained from green Antarctic algae Prasiola crispa against equine herpesvirus 1. Environ. Boil. Fishes 2017, 29, 555–562. [Google Scholar] [CrossRef]
  73. Tumilaar, S.G.; Siampa, J.P.; Fatimawali; Kepel, B.J.; Niode, N.J.; Idroes, R.; Rakib, A.; Emran, T.B.; Tallei, T.E. Potential of leaf extract of Pangium edule Reinw as HIV-1 protease inhibitor: A computational biology approach. J. Appl. Pharma. Sci. 2021, 11, 101–110. [Google Scholar] [CrossRef]
  74. Mahmud, S.; Paul, G.K.; Afroze, M.; Islam, S.; Gupt, S.B.R.; Razu, M.H.; Biswas, S.; Zaman, S.; Uddin, M.S.; Khan, M.; et al. Efficacy of Phytochemicals Derived from Avicennia officinalis for the Management of COVID-19: A Combined In Silico and Biochemical Study. Molecules 2021, 28, 2210. [Google Scholar] [CrossRef] [PubMed]
  75. Gershell, L.J.; Atkins, J.H. A brief history of novel drug discovery technologies. Nat. Rev. Drug Discov. 2003, 2, 321–327. [Google Scholar] [CrossRef] [PubMed]
  76. Keller, T.H.; Pichota, A.; Yin, Z. A practical view of ‘druggability’. Curr. Opin. Chem. Biol. 2006, 10, 357–361. [Google Scholar] [CrossRef] [PubMed]
  77. Vistoli, G.; Pedretti, A.; Testa, B. Assessing drug-likeness–What are we missing? Drug Discov. Today 2008, 13, 285–294. [Google Scholar] [CrossRef]
  78. Khan, J.; Sakib, S.A.; Mahmud, S.; Khan, Z.; Islam, M.N.; Sakib, M.A.; Emran, T.B.; Simal-Gandara, J. Identification of potential phytochemicals from Citrus limon against main protease of SARS-CoV-2: Molecular docking, molecular dynamic simulations and quantum computations. J. Biomol. Struct. Dyn. 2021, 2021, 1–12. [Google Scholar] [CrossRef]
Figure 1. Total ionic chromatogram (TIC) for the methanolic extract of the aerial parts of L. zeylanica via gas chromatography–mass spectrometry (GC-MS).
Figure 1. Total ionic chromatogram (TIC) for the methanolic extract of the aerial parts of L. zeylanica via gas chromatography–mass spectrometry (GC-MS).
Biology 10 00789 g001
Figure 2. Shown here is 11-oxa-dispiro[4.0.4.1]undecan-1-ol binding the SARS-CoV-2 receptor (PDB ID: 6LU7). (A) A 3D representation and (B) 2D representation. Hydrogen bonds are displayed as green balls and sticks, hydrophobic bonds (Pi–Pi/Pi–sigma/amide–Pi interactions) are displayed as violet balls and sticks, hydrophobic bonds (Pi–alkyl/alkyl interactions stacking) are displayed as pink balls and sticks, hydrophobic bonds (Pi–sulfur) are displayed as gold balls and sticks, and carbon–hydrogen bonds are displayed as white balls and sticks.
Figure 2. Shown here is 11-oxa-dispiro[4.0.4.1]undecan-1-ol binding the SARS-CoV-2 receptor (PDB ID: 6LU7). (A) A 3D representation and (B) 2D representation. Hydrogen bonds are displayed as green balls and sticks, hydrophobic bonds (Pi–Pi/Pi–sigma/amide–Pi interactions) are displayed as violet balls and sticks, hydrophobic bonds (Pi–alkyl/alkyl interactions stacking) are displayed as pink balls and sticks, hydrophobic bonds (Pi–sulfur) are displayed as gold balls and sticks, and carbon–hydrogen bonds are displayed as white balls and sticks.
Biology 10 00789 g002
Figure 3. Azetidin-2-one 3,3-dimethyl-4-(1-aminoethyl) binding the SARS-CoV-2 receptor (PDB ID: 6LU7). (A) A 3D representation and (B) 2D representation. Hydrogen bonds are displayed as green balls and sticks, hydrophobic bonds (Pi–Pi/Pi–sigma/amide–Pi interactions) are displayed as violet balls and sticks, hydrophobic bonds (Pi–alkyl/alkyl interactions stacking) are displayed as pink balls and sticks, hydrophobic bonds (Pi–sulfur) are displayed as gold balls and sticks, and carbon–hydrogen bonds are displayed as white balls and sticks.
Figure 3. Azetidin-2-one 3,3-dimethyl-4-(1-aminoethyl) binding the SARS-CoV-2 receptor (PDB ID: 6LU7). (A) A 3D representation and (B) 2D representation. Hydrogen bonds are displayed as green balls and sticks, hydrophobic bonds (Pi–Pi/Pi–sigma/amide–Pi interactions) are displayed as violet balls and sticks, hydrophobic bonds (Pi–alkyl/alkyl interactions stacking) are displayed as pink balls and sticks, hydrophobic bonds (Pi–sulfur) are displayed as gold balls and sticks, and carbon–hydrogen bonds are displayed as white balls and sticks.
Biology 10 00789 g003
Figure 4. Lorazepam, 2TMS derivative binding the SARS-CoV-2 receptor (PDB ID: 6LU7). (A) A 3D representation and (B) 2D representation. Hydrogen bonds are displayed as green balls and sticks, hydrophobic bonds (Pi–Pi/Pi–sigma/amide–Pi interactions) are displayed as violet balls and sticks, hydrophobic bonds (Pi–alkyl/alkyl interactions stacking) are displayed as pink balls and sticks, hydrophobic bonds (Pi–sulfur) are displayed as gold balls and sticks, and carbon–hydrogen bonds are displayed as white balls and sticks.
Figure 4. Lorazepam, 2TMS derivative binding the SARS-CoV-2 receptor (PDB ID: 6LU7). (A) A 3D representation and (B) 2D representation. Hydrogen bonds are displayed as green balls and sticks, hydrophobic bonds (Pi–Pi/Pi–sigma/amide–Pi interactions) are displayed as violet balls and sticks, hydrophobic bonds (Pi–alkyl/alkyl interactions stacking) are displayed as pink balls and sticks, hydrophobic bonds (Pi–sulfur) are displayed as gold balls and sticks, and carbon–hydrogen bonds are displayed as white balls and sticks.
Biology 10 00789 g004
Figure 5. Standard inhibitor N3 binding the SARS-CoV-2 receptor (PDB ID: 6LU7). (A) A 3D representation and (B) 2D representation. Hydrogen bonds are displayed as green balls and sticks, hydrophobic bonds (Pi–Pi/Pi–sigma/amide–Pi interactions) are displayed as violet balls and sticks, hydrophobic bonds (Pi–alkyl/alkyl interactions stacking) are displayed as pink balls and sticks, hydrophobic bonds (Pi–sulfur) are displayed as gold balls and sticks, and carbon–hydrogen bonds are displayed as white balls and sticks.
Figure 5. Standard inhibitor N3 binding the SARS-CoV-2 receptor (PDB ID: 6LU7). (A) A 3D representation and (B) 2D representation. Hydrogen bonds are displayed as green balls and sticks, hydrophobic bonds (Pi–Pi/Pi–sigma/amide–Pi interactions) are displayed as violet balls and sticks, hydrophobic bonds (Pi–alkyl/alkyl interactions stacking) are displayed as pink balls and sticks, hydrophobic bonds (Pi–sulfur) are displayed as gold balls and sticks, and carbon–hydrogen bonds are displayed as white balls and sticks.
Biology 10 00789 g005
Figure 6. The molecular dynamics simulation of the four docked complexes. (AF) Present, RMSD, Rg, SASA, hydrogen bonds, RMSF, and MM-PBSA analyses, respectively. Here, lorazepam, 11-oxa-dispiro[4.0.4.1]undecan-1-ol, azetidin-2-one-3, 3-dimethyl-4-(1-aminoethyl), and inhibitor N3 were denoted as D1, D2, D3, and control, respectively.
Figure 6. The molecular dynamics simulation of the four docked complexes. (AF) Present, RMSD, Rg, SASA, hydrogen bonds, RMSF, and MM-PBSA analyses, respectively. Here, lorazepam, 11-oxa-dispiro[4.0.4.1]undecan-1-ol, azetidin-2-one-3, 3-dimethyl-4-(1-aminoethyl), and inhibitor N3 were denoted as D1, D2, D3, and control, respectively.
Biology 10 00789 g006
Table 1. List of compounds identified from the L. zeylanica methanolic extract using gas chromatography–mass spectrometry analysis.
Table 1. List of compounds identified from the L. zeylanica methanolic extract using gas chromatography–mass spectrometry analysis.
Sl.NameRTm/zAreaConc.Peak Area (%)
1.3-Butynoic acid5.92740.0036,1350.0490.049299
2.Tetradecamethylcycloheptasiloxane8.22140.0023,0610.0310.031462
3.Trimethylsilyl 2,6-bis[(trimethylsilyl)oxy]benzoate10.13973.001,516,8952.0662.069501
4.Bis(heptamethylcyclotetrasiloxy)hexamethyltrisiloxane11.91473.001,205,6831.6421.644915
5.Lorazepam, 2TMS derivative11.91473.001,205,6831.6421.644915
6.Cyanoacetic acid10.66940.0049,2310.0670.067166
7.Methyl 11-bromoundecanoate11.07840.0045,4300.0620.06198
8.Azetidin-2-one 3,3-dimethyl-4-(1-aminoethyl)11.07840.0045,4300.0620.06198
9.3-Azabicyclo[3.2.2]nonane11.91640.0067,3270.0920.091854
10.Phytol acetate12.56768.002,737,1013.7283.73423
11.Hexadecanoic acid, methyl ester13.50274.0014,101,16119.20719.23823
12.11-Oxa-dispiro[4.0.4.1]undecan-1-ol13.50340.0091,6780.1250.125076
13.Hexadecamethylcyclooctasiloxane14.94673.00789,7721.0761.077487
14.9,12-Octadecadienoic acid (Z,Z)-, methyl ester15.22567.006,028,1388.2118.224195
15.6-Octadecenoic acid, methyl ester, (Z)-15.22555.002,986,8864.0684.075012
16.Phytol15.39371.0011,551,44015.73415.75964
17.Methyl stearate15.53074.002,412,1523.2863.290901
18.Octadecamethylcyclononasiloxane16.25773.00913,5211.2441.246318
19.Pseduosarsasapogenin-5,20-dien17.21983.00547,3130.7450.7467
20.Cyclodecasiloxane, eicosamethyl-19.54073.00929,8091.2661.268539
21.Hexadecanoic acid, 2-hydroxy-1-(hydroxymethyl)ethyl ester20.21198.00971,0871.3231.324855
22.Octadecanoic acid, 2,3-dihydroxypropyl ester23.85143.00356,2720.4850.486062
23.13-Docosenamide, (Z)-24.59659.0015,302,72520.84420.87752
24.N,N′-methylenebis(oleamide), (Z,Z)-24.595207.00247,5290.3370.337704
25.Squalene27.770207.00138,6190.1890.189118
26.α-Tocopheryl acetate28.073207.00203,3410.2770.277418
27.Campesterol29.656207.00201,5720.2750.275005
28.Stigmasterol30.546207.00156,0660.2130.212921
29.γ-Sitosterol31.581207.00165,3100.2250.225533
30.4-Campestene-3-one31.581207.00147,9600.2020.201862
31.9, 19-Cyclolanost-24-en-3-ol, acetate, (3.beta.)32.728207.00101,4950.1380.13847
32.4,22-Cholestadien-3-one32.728207.00101,4950.1380.13847
33.Stigmast-4-en-3-one33.655124.004,766,3236.4926.502699
34.1,2-Bis(trimethylsilyl)benzene35.243207.00163,5490.2230.22313
Table 2. Docking results (kcal/mol)) estimated for the top 34 compounds.
Table 2. Docking results (kcal/mol)) estimated for the top 34 compounds.
SL. No.NameDocking ScoreInteractions by H-BondHydrophobic Bonds (Pi–Alkyl/Alkyl Interaction)Hydrophobic Bonds (Pi–Pi/Pi–Sigma/ Amide–Pi Interaction)Pi–Sulfur Interaction
1.3-Butynoic acid−1.08Cys 145, Ser 144 (2), Gly 143His 172, His 163Phe 140
2.Tetradecamethylcycloheptasiloxane
3.Trimethylsilyl 2,6-bis[(trimethylsilyl)oxy]benzoate
4.Bis(heptamethylcyclotetrasiloxy)hexamethyltrisiloxane
5.Lorazepam, 2TMS derivative−5.246Gly 143Cys 145Cys 145, Met 49
6.Cyanoacetic acid−3.469Lys 61 (2), Arg 60, Asp 48
7.Methyl 11-bromoundecanoate−1.792Gly 143His 163, Cys 145
8.Azetidin-2-one 3,3-dimethyl-4-(1-aminoethyl)−5.39Gln 189, Tyr 54Met 165 (2)His 41
9.3-Azabicyclo[3.2.2]nonane−4.703His 41 (2), Met 49, Met 165
10.Phytol acetate−3.357Met 165, Glu 166Met 165, Cys 145 (2), Leu 27, His 41, Met 49, Met 165
11.Phytol−1.565Asn 142,His 163, His 172, Met 165 (2), His 41, Met 49
12.Hexadecanoic acid, methyl ester+0.991Asn 142Met 165,His 41
13.11-Oxa-dispiro[4.0.4.1]undecan-1-ol−5.755Glu 166Met 165
14.Hexadecamethylcyclooctasiloxane
15.9,12-Octadecadienoic acid (Z,Z)-, methyl ester−1.111Glu 166Met 165, Leu 167
16.6-Octadecenoic acid, methyl ester, (Z)-−0.399Ser 144, Gly 143Cys 145, His 163His 41
17.Methyl stearate−0.419Gln 189His 41, Leu 50, Met 49
18.Octadecamethylcyclononasiloxane
19.Pseduosarsasapogenin-5,20-dien
20.Cyclodecasiloxane, eicosamethyl-
21.Hexadecanoic acid, 2-hydroxy-1-(hydroxymethyl)ethyl ester−4.152Glu 166, Cys 145
22.Octadecanoic acid, 2,3-dihydroxypropyl ester−2.406His 163, Gln 143 (2), Ser 144, Cys 145
23.13-Docosenamide, (Z)-−4.46His 163, Cys 145His 41, Met 49
24.N,N’′methylenebis(oleamide), (Z, Z)-−4.057Phe 140His 41
25.Squalene−3.609Met 165 (2), Met 49, Cys 149, His 41
26.α-Tocopheryl acetate−4.871Leu 167, Met 165 (2), Met 49, His 41, Leu 27, Cys 145
27.Campesterol−3.776Thr 24Met 49 (2), His 41 (2), Met 165 (2), Cys 145
28.Stigmasterol−4.194Cys 145, Ser 144 (2)Ala 191, Pro 168
29.γ-Sitosterol−4.854Thr 26Cys 145 (3), His 41 (2), Met 49, Met 165 (2), Leu 167, Pro 168
30.4-Campestene-3-one−3.934Met 165, Met 49 (2), His 41 (2), Cys 145
31.9, 19-Cyclolanost-24-en-3-ol, acetate, (3.beta.)−3.105Asn 142, Ser 144Cys 145, Pro 168, Ala 191, Leu 50
32.4,22-Cholestadien-3-one−3.824Ser 144Ala 191,
33.Stigmast-4-en-3-one−3.696 His 41 (2), Cys 145, Leu 27
34.1,2-Bis(trimethylsilyl)benzene
35.Standard (inhibitor N3)−7.013Phe 140, Gly 143, His 164, Glu 166, Gln 189, Thr 190His 41, Met 49, Met 165, Leu 167, Pro 168, Ala 191
N.B.: Bold text indicates best docking scores.
Table 3. Pharmacological profile of the 23 ligand molecules passing the Lipinski’s RO5, as derived from the SwissADME webserver.
Table 3. Pharmacological profile of the 23 ligand molecules passing the Lipinski’s RO5, as derived from the SwissADME webserver.
CompoundsMolecular Weight a (g/mol)Hydrogen Bond Acceptors bHydrogen Bond Donors cMlogP dMolar Refractivity eNo. of Lipinski Violations f
<500≤10≤5<540–130≤1
3-Butynoic acid84.07210.3821.281
Tetradecamethylcycloheptasiloxane519.0870−1.54129.971
Trimethylsilyl 2,6-bis[(trimethylsilyl)oxy]benzoate370.66402.97103.150
Lorazepam, 2TMS derivative465.52404.66134.901
Hexadecanoic acid, 2-hydroxy-1-(hydroxymethyl)ethyl ester330.50423.1897.060
Cyanoacetic acid85.0631−0.9618.061
Methyl 11-bromoundecanoate279.21203.5668.950
Azetidin-2-one 3,3-dimethyl-4-(1-aminoethyl)-142.20220.7243.010
3-Azabicyclo[3.2.2]nonane125.21111.8343.060
Phytol, acetate338.57205.47108.681
Hexadecanoic acid, methyl ester270.45204.4485.120
11-Oxa-dispiro[4.0.4.1]undecan-1-ol168.23211.5246.160
9,12-Octadecadienoic acid (Z,Z)-, methyl ester294.47204.7093.780
6-Octadecenoic acid, methyl ester, (Z)-296.49204.8094.260
Phytol296.53115.2598.941
Methyl stearate298.50204.9194.730
Pseduosarsasapogenin-5,20-dien414.62324.42123.270
Octadecanoic acid, 2,3-dihydroxypropyl ester358.56423.63106.670
13-Docosenamide, (Z)-337.58115.06110.301
Campesterol400.68116.54128.421
4-Campestene-3-one398.66106.43127.461
4,22-Cholestadien-3-one382.62106.13122.181
1,2-Bis(trimethylsilyl)benzene222.47004.1372.400
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dutta, M.; Tareq, A.M.; Rakib, A.; Mahmud, S.; Sami, S.A.; Mallick, J.; Islam, M.N.; Majumder, M.; Uddin, M.Z.; Alsubaie, A.; et al. Phytochemicals from Leucas zeylanica Targeting Main Protease of SARS-CoV-2: Chemical Profiles, Molecular Docking, and Molecular Dynamics Simulations. Biology 2021, 10, 789. https://doi.org/10.3390/biology10080789

AMA Style

Dutta M, Tareq AM, Rakib A, Mahmud S, Sami SA, Mallick J, Islam MN, Majumder M, Uddin MZ, Alsubaie A, et al. Phytochemicals from Leucas zeylanica Targeting Main Protease of SARS-CoV-2: Chemical Profiles, Molecular Docking, and Molecular Dynamics Simulations. Biology. 2021; 10(8):789. https://doi.org/10.3390/biology10080789

Chicago/Turabian Style

Dutta, Mycal, Abu Montakim Tareq, Ahmed Rakib, Shafi Mahmud, Saad Ahmed Sami, Jewel Mallick, Mohammad Nazmul Islam, Mohuya Majumder, Md. Zia Uddin, Abdullah Alsubaie, and et al. 2021. "Phytochemicals from Leucas zeylanica Targeting Main Protease of SARS-CoV-2: Chemical Profiles, Molecular Docking, and Molecular Dynamics Simulations" Biology 10, no. 8: 789. https://doi.org/10.3390/biology10080789

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop