Next Article in Journal
Electroanalytical Techniques for the Detection of Selenium as a Biologically and Environmentally Significant Analyte—A Short Review
Next Article in Special Issue
Recent Advances in Discovery of Lead Structures from Microbial Natural Products: Genomics- and Metabolomics-Guided Acceleration
Previous Article in Journal
Assessing Green Methods for Pectin Extraction from Waste Orange Peels
Previous Article in Special Issue
Holospiniferoside: A New Antitumor Cerebroside from The Red Sea Cucumber Holothuria spinifera: In Vitro and In Silico Studies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Potential SARS-CoV-2 Main Protease and Spike Protein Inhibitors from the Genus Aloe: An In Silico Study for Drug Development

by
Mohamed E. Abouelela
1,
Hamdy K. Assaf
1,
Reda A. Abdelhamid
1,
Ehab S. Elkhyat
1,
Ahmed M. Sayed
2,
Tomasz Oszako
3,
Lassaad Belbahri
4,*,
Ahmed E. El Zowalaty
5,6,* and
Mohamed Salaheldin A. Abdelkader
7,*
1
Department of Pharmacognosy, Faculty of Pharmacy, Al-Azhar University, Assiut-Branch, Assiut 71524, Egypt
2
Department of Pharmacognosy, Faculty of Pharmacy, Nahda University, Beni-Suef 62513, Egypt
3
Department of Forest Protection, Forest Research Institute, 05-090 Sekocin Stary, Poland
4
Laboratory of Soil Biology, University of Neuchatel, 2000 Neuchatel, Switzerland
5
Sahlgrenska Center for Cancer Research, Department of Surgery, Institute of Clinical Sciences, University of Gothenburg, 405 30 Gothenburg, Sweden
6
Wallenberg Centre for Molecular and Translational Medicine, University of Gothenburg, 405 30 Gothenburg, Sweden
7
Department of Pharmacognosy, Faculty of Pharmacy, Sohag University, Nasr City, Sohag 82524, Egypt
*
Authors to whom correspondence should be addressed.
Molecules 2021, 26(6), 1767; https://doi.org/10.3390/molecules26061767
Submission received: 24 February 2021 / Revised: 16 March 2021 / Accepted: 18 March 2021 / Published: 21 March 2021

Abstract

:
Severe acute respiratory syndrome coronavirus (SARS-CoV-2) disease is a global rapidly spreading virus showing very high rates of complications and mortality. Till now, there is no effective specific treatment for the disease. Aloe is a rich source of isolated phytoconstituents that have an enormous range of biological activities. Since there are no available experimental techniques to examine these compounds for antiviral activity against SARS-CoV-2, we employed an in silico approach involving molecular docking, dynamics simulation, and binding free energy calculation using SARS-CoV-2 essential proteins as main protease and spike protein to identify lead compounds from Aloe that may help in novel drug discovery. Results retrieved from docking and molecular dynamics simulation suggested a number of promising inhibitors from Aloe. Root mean square deviation (RMSD) and root mean square fluctuation (RMSF) calculations indicated that compounds 132, 134, and 159 were the best scoring compounds against main protease, while compounds 115, 120, and 131 were the best scoring ones against spike glycoprotein. Compounds 120 and 131 were able to achieve significant stability and binding free energies during molecular dynamics simulation. In addition, the highest scoring compounds were investigated for their pharmacokinetic properties and drug-likeness. The Aloe compounds are promising active phytoconstituents for drug development for SARS-CoV-2.

1. Introduction

SARS-CoV-2, a novel coronavirus disease caused by Severe Acute Respiratory Syndrome Coronavirus 2, an RNA β-coronavirus, poses an increasing threat to human health. To date, SARS-CoV-2 has infected 105 million people worldwide (https://www.worldometers.info/, accessed on 20 February 2021). The disease is manifested by fever, cough, dyspnea and pneumonia with unknown etiology that worsens over time and can lead to death [1,2]. In addition, SARS-CoV-2 patients develop low levels of neutralizing antibodies leading to prolonged disease [3]. Entry of SARS-CoV-2 into host cells is a critical factor in its pathogenesis. The surface-anchored spike proteins of SARS-CoV-2 are key determinants of viral entry. They bind to surface receptors on host cells, then the virus enters endosomal pathway followed by fusion of viral and lysosomal membranes [4,5]. SARS-CoV-2 spike protein has N-terminus S1 domain, a receptor binding domain (RBD) that recognizes and binds angiotensin-converting enzyme 2 (ACE2) [6,7,8]. Proteolytic activation of SARS-CoV-2 spike proteins is mediated by the cell surface protease TMPRSS2, a process crucial for membrane fusion and viral entry [9]. Recent studies have shown that viral entry depends on the component of the head spike that recognizes the ACE2 receptor. In addition, structural and energetic analysis have shown that high-frequency contacts between ACE2 and SARS-CoV-2 spike protein lead to local conformational stability and large energetic cost was required for virus-cell collision at early stage facilitating cell entry [8,10,11].Coronavirus main protease (Mpro, also known as 3CLpro), is the best characterized drug target, with no known human protease having the same cleavage pattern, so its inhibition leads to specific blockade of viral replication [12]. Spike glycoprotein (S protein) and its RBD are important targets for therapeutic intervention that target host cell recognition and the membrane fusion process [13]. As a global health emergency, abundant collaborative efforts have rapidly emerged to investigate the effectiveness of different therapies as antiviral, monoclonal antibodies, immune-therapies, and vaccines [14]. Current antiviral therapies for other viruses as SARS-CoV-1, MERS-CoV and HIV as well as antimalarial drugs have been inspected for their activity against SARS-CoV-2. For instance, the antimalarial drug hydroxychloroquine blocks viral cell entry by inhibiting glycosylation of host receptors and proteolytic processing. In addition, the antiviral Favipiravir inhibits RNA polymerase and is involved in entry blocking. Both drugs showed potential in vitro activity against SARS-CoV-2 [15].
Although at least six SARS-CoV-2 vaccines have been developed and licensed for emergency use, the safety, efficacy, durability and availability to large populations have not been established, so it is too early to know if COVID-19 vaccines will provide long-term protection. In addition, there is still no effective drug therapy for SARS-CoV-2. The current therapeutic strategies depend on supportive therapy and symptomatic management. Natural products can serve as prophylactic agents, halt virus progression, inhibit inflammatory cytokines secretion, and reduce infection, complications and mortality of SARS-CoV-2 [16].
Natural products have been a valuable source of therapeutic agents, molecules with therapeutic potentials, and an important source of more efficient drugs that are based on the chemical structure of natural products. For example, flavonoids have shown significant antiviral activities [17]. Curcumin and luteolin also show therapeutic potential against HIV targeting viral protease and HIV-1 transactivator of transcription [18,19]. Kaempferol also exhibits anti HSV-1 and 2 activities [20,21] Aloe is an ancient common plant species used as a medicinal plant. The genus Aloe comprises about 581 species. Its pharmacological properties and phytochemical characteristics have been extensively studied and evaluated [22]. Previous studies have shown that natural products from Aloe possess anti-inflammatory, immunostimulant, anti-cancer, antioxidant, anti-ageing, wound healing, antifungal, antibacterial and antiviral activities [23]. Natural products from Aloe showed antiviral and inhibitory activities against HSV-1 and 2, human cytomegalovirus (HCMV), influenza A, polio and other hemagglutinating viruses [23,24,25]. Aloin, a major compound of Aloe species, significantly reduces influenza viruses replication including oseltamivir-resistant (H1N1) influenza virus [26]. This indicates that Aloe genus offer a rich source of potential anti-viral compounds.
Virtual screening and molecular modeling studies showed potential therapeutic activities of some natural products in inhibiting SARS-CoV-2 proteins including the main protease (Mpro), spike glycoprotein (S) and angiotensin converting enzyme-2 (ACE2) receptor which are promising potential therapeutic targets [22,27]. In the present study, we conducted computational screening and molecular dynamics study on a library of isolated molecules from Aloe genus, investigated the binding affinity of these compounds with SARS-CoV-2 main protease (Mpro), spike glycoprotein (S) through molecular docking analysis. We found six potential inhibitors from Aloe genus that effectively bind to SARS-CoV-2 main protease (Mpro), and three inhibitors that effectively bind spike glycoprotein receptor binding domain-ACE2 interface.

2. Results and Discussion

2.1. Phytochemical Constituents of Aloe

The tested library of active constituents from Aloe genus comprised phytochemicals that cover major classes of natural products (Figure 1). Phytochemical studies of the genus Aloe plants showed the presence of anthraquinones, chromones, coumarins, flavonoids, simple phenolic compounds, phenyl pyrans and phenyl pyrones, benzofurans, naphthalene derivatives, alkaloids and fatty acid derivatives (Table A1 and Figures S1–S18). Out of 237 compounds compiled in the library; anthraquinones were the most abundant constituents, with a percentage of 36.29%, followed by chromones (27.43%) and simple phenolic compounds (7.17%), while alkaloids, coumarin and fatty acid derivatives constituents were less abundant.

2.2. Structure-Based Virtual Screening and Molecular Docking of Aloe Phytochemicals on SARS-CoV-2 Spike Glycoprotein and Main Protease

High-throughput virtual screening of compounds from Aloe, was followed by molecular docking and MD simulation. Since ligand binding to a protein of interest is the first step in drug discovery, molecular docking is widely used to predict and identify ligands that fit into the binding pocket of a protein of interest [28]. Our screening was performed against two major drug discovery and therapeutic targets of SARS-CoV-2, spike glycoprotein and Mpro proteins [7,12]. SARS-CoV-2 main protease Mpro is critical for the life cycle of the virus. Approximately, two thirds of the SARS-CoV-2 genome is translated into polyproteins pp1a and pp1ab, that are cleaved with Mpro into nonstructural proteins that are involved in the production of viral membrane, spike and nucleocapsid proteins [29]. Mpro is a dimer that has cysteine and histidine in the active site which form a catalytic dyad, conserved among coronaviruses making it an ideal therapeutic target [12]. In molecular docking studies, the ligand-receptor interaction with protein active site residues is established by formation of some interactions including hydrogen bonds, Van der Waal force interaction, π-sigma bond, π–π interaction, electrostatic interaction, and many other hydrophobic interactions. Hydrogen bonds are essential for interaction, lowering the binding energy and stabilizing the ligand-receptor docked complex. Pharmacologically, it is well-known that blockade of a receptor active site by a ligand terminates its functional activity [30]. Our molecular docking approach was validated by docking of hydroxychloroquine, a potent inhibitor of SARS-CoV-2 Mpro. Hydroxychloroquine acts as a lysomotropic agent that inhibits viral entry and viral endocytosis. Viral entry and replication are highly dependent on the acidic pH of lysosomes and endosomes, and some host proteases which are also active in acidic pH (pH 5–5.5) [31]. Chloroquine and its analogues are diprotic weak bases that in their unprotonated forms, readily diffuse through cellular and organelle membranes such as lysosomes, endosomes and Golgi vesicles increasing pH from 6.3 to 6.7 [32,33,34]. In addition to disruption of endocytic pathway pH, chloroquine and hydroxychloroquine have been recently found to be potent inhibitors of SARS-CoV-2 Mpro but not viruses that belong to Rhabdoviridae [35]. In our study, the compounds previously isolated from Aloe plants were virtually screened against SARS-CoV-2 main protease Mpro (PDB ID: 6LU7) (Figure 2) and spike glycoprotein (PDB ID: 6M0J) (Figure 2) to find potential inhibitors for SARS-CoV-2. Using our docking approach, hydroxychloroquine interacted with SARS-CoV-2 protein Mpro and docked hydroxychloroquine bound to the active site with and RMSD of 1.2 Å. Molecular docking data were filtered to remove compounds with scores > −6.5 for both SARS-CoV-2 main protease Mpro (Figure 3 and Table A1) and spike glycoprotein (Figure 4 and Table A1). Molecular docking was performed by examining the interactions of these compounds with the active site residues of these proteins and analysis of results.
Compounds scoring lower than −5.00 kcal/mol are expected to be active. These compounds were then filtered by RMSD value [30], to evaluate experimental stability of the docked ligand conformers. RMSD values around 1.5 Å, are considered successful and stable while those beyond 2 Å indicate instability of ligand conformation and docking parameters [36]. For SARS-CoV-2 protein Mpro, the binding energy observed for these compounds ranged from−7.950 to −0.339 kcal/mol while for spike glycoprotein, binding energy ranged from −8.088 to −5.437 kcal/mol. The top three scoring compounds for SARS-CoV-2 protein Mpro were compound 132 (2′-oxo-2′-O-(3,4-dihydroxy-E-cinnamoyl)-(2′R) aloesinol-7-methyl ether), compound 134 (2′-oxo-2′-O-(4-hydroxy-3-methoxy-(E)-cinnamoyl)-(2′R)-aloesinol-7-methyl ether) and compound 159 (rutin), (Table 1 docking scores and Figure 5, top panel). These three compounds showed the strongest interaction with the active site of SARS-CoV-2 protein Mpro. Molecular 2D and 3D interactions complexes of compounds 132, 134 and 159 with SARS- SARS-CoV-2 protein Mpro are shown in Figure 6.
On the other hand, the top three scoring compounds for SARS-CoV-2 spike glycoprotein were compounds compound 115 (2″-O-(4-methoxycinnamoyl)-(S)-aloesinol), compound 120 (rabaichromone), and compound 131 (aloeribide), (Table 1 docking score, and Figure 5, lower panel).
These three compounds showed the strongest interaction with SARS-CoV-2 RBD. Molecular 2D and 3D interactions complexes of compounds 115, 120 and 131 with SARS- SARS-CoV-2 protein spike glycoprotein receptor binding domain are shown in Figure 7. In depth analysis showed that chromone derivatives 132 and 134 had high binding affinity as lead compounds for developing SARS-CoV-2 Mpro inhibitors. These compounds had a score of −7.683 kcal/mol (RSMD = 1.37) and −7.951 kcal/mol (RSMD = 1.72), respectively (Table 1).
The interacting residues of SARS-CoV-2 Mpro involved in interactions with compound 132 were ASN142A, ASN142A, HIS163A, GLN189A, GLU166A and GLN189A (Figure 6A1,B1 and Table 1), while compound 134 interacted with ASN142A, ASN142A and HIS163 (Figure 6A2,B2 and Table 1). In addition, the flavonoid compound 159 forms a hydrogen bond with the Mpro protein THR190A amino acid residue with −7.728 kcal/mol as scoring value (Figure 6A3,B3 and Table 1). Regarding the interaction of the compounds with SARS-CoV-2 spike glycoprotein, this was mainly supported by hydrogen bonds, π-H, ionic and hydrophobic interactions. The highest scoring compounds were the chromone derivatives 115 (2″-O-(4-methoxycinnamoyl)-(S)-aloesinol), 120 (rabaichromone) and 131 (aloeribide) with binding energies −8.057, −7.871 and −8.088 kcal/mol, respectively. Compound 115 interacted with SARS-CoV-2 RBD TRP566A, LYS562A, LYS562A and VAL209 (Figure 7A1,B1 and Table 1). Moreover, compound 131 interacted with GLN102A, ASN210A and ASP206A (Figure 7A2,B2 and Table 1) by π i- H bond while compound 120 formed hydrogen bonds with ALA396A, ASP206A and GLU208A amino acid residues (Figure 7A3,B3 and Table 1).

2.3. Molecular Dynamics Simulation

Conventional docking approaches do not account for the inherent protein binding site flexibility and the many protein conformational rearrangements [37]. Computational tools for drug discovery such as molecular dynamics take into account structural flexibility and entropic effects which produce accurate predictions of small molecule-protein binding thermodynamics and kinetics [38]. Hence dynamical docking considers flexibility of drug-protein binding and conformational changes, solvation of drug-protein complex and temperature [38,39]. Unbiased millisecond-long can predict spontaneous drug-protein entire binding [40]. In addition, recent developments in dynamical docking such as enhanced sampling for dynamical docking, path-based and alchemical transformations have greatly impacted drug discovery [38]. To validate molecular docking results, we subjected the top scoring compounds to unbiased molecular dynamics simulation experiments. The three top scoring Mpro inhibitor hits 132, 134, and 159 were able to achieve stable binding inside the active site with low deviations across the course of simulations (Average RMSD = 3.22, 3.32, and 3.86 Å, respectively) and convergent binding free energies (ΔG = −6.9, −6.8, and −6.5 kcal/mol, respectively), (Figure 8A).
With respect to SARS-CoV-2 spike glycoprotein, both compounds 120 and 131 were stable inside the binding site during MD simulation, with scoring average RMSDs of 2.81 Å and 3.96 Å, respectively, and ΔG of −7.4 and −6.8 kcal/mol, respectively (Figure 8B). On the other hand, compound 115 was significantly less stable (average RMSD = 6.2 Å) inside the SARS-CoV-2 spike glycoprotein binding site, and this instability was further translated into a low binding free energy (ΔG = −4.5 kcal/mol) compared to compounds 120 and 131 (Figure 8B). RMSF is an expression of the average residual mobility throughout simulation in a structure and a higher RMSF value indicates more flexibility during MD simulation. We calculated the RMSF value for the top scoring compounds from Aloe genus with SARS-CoV-2 Mpro and SARS-CoV-2 spike glycoprotein and plotted RMSF value versus residue number (Figure 8C,D). The results indicate that compounds 159 and 120 had high RMSF values compared to other compounds. The RMSD and RMSF values indicate that the top scoring compounds from Aloe genus were stable and had greater random motion during the simulation. The inhibitors identified in in our docking analysis that showed interaction with SARS CoV-2 spike protein and MPro are in agreement with previously reported results [41]. Arokiyaraj et al. found that several polyphenolic compounds from Geranii Herba, including geraniin, kaempferitrin, quercitin, gallic acid, and kaempferol interacted with amino acid residues in the SARS-CoV-2 RBD active site inhibiting the interaction of SARS-CoV-2 RBD with ACE2. Arokiyaraj et al. also reported that these polyphenolic compounds interacted strongly with amino acids in the active site of SARS-CoV-2 Mpro and its proximity leading to blockade of the nucleophilic attack toward His 41 and blockade of proteolytic activity. In agreement with this, we found that quercetin interacted with SARS CoV-2 RBD and Mpro with binding energies of −5 Kcal/mol and −5.5 Kcal/mol respectively, similar to the results reported by Arokiyaraj et al. for quercetin interaction with RBD and Mpro −5.71 Kcal/mol and −6.49 kcal/mol, respectively. In addition, we found that gallic acid interacted with SARS CoV-2 RBD and Mpro with binding energies of −4.19 Kcal/mol and −3.56 Kcal/mol, respectively, similar to the binding energies reported by Arokiyaraj et al. for the gallic acid interaction with SARS CoV-2 RBD and Mpro, −4.21 kcal/mol and −4.46 kcal/mol, respectively. These finding indicate that phenolic compounds from Aloe are potential inhibitors for SARS CoV-2 RBD and Mpro [41].

2.4. Drug like Properties, and Pharmacokinetic Prediction of the Ligands

Drug-like properties and pharmacokinetic properties are intrinsic characteristics of drugs that may need to be optimized independently from pharmacodynamics properties during drug development. It is a balance among molecular properties affecting pharmacodynamics and pharmacokinetics of small molecules. These molecular properties such as membrane permeability and bioavailability are always connected to some basic molecular descriptors such as lipophilicity log P, (Tendency of a compound to partition into an aqueous matrix versus lipid matrix), molecular weight (MW), topological polar surface area (TPSA), or hydrogen bond acceptors and donors count in a molecule. Lipophilicity impacts drug’s absorption, distribution, metabolism, elimination (ADME) and plasma protein binding properties. In addition, the number of hydrogen bond donors and hydrogen bond acceptors influence drug’s pKa (−log Ka). The solubility of small molecules impacts their bioavailability and the need for frequent dosing, hence we investigated the ADME properties, inhibition of cytochrome P450 (CYP), modulation of P-glycoprotein (Pgp), solubility, plasma protein binding and permeability of the top scoring compounds in our analysis. The best scoring compounds for both SARS-CoV-2 Mpro and spike glycoprotein were tested for obeying Lipinski’s rule of five parameters, which states that drugs having log P ranging from 0 to 5, have high possibility of oral absorption [42]. Data (Table 2) showed that the compounds have log P values that ranged from −1.06 to 2.8 that does not exceed 5.0 indicating reasonable probability of their good absorption. The number of hydrogen bond donors was variable and ranged from 4 to10 that is more than 5 and also hydrogen bond acceptors were 11–16 that is more than 10. All compounds have number of atoms that ranged from 40 to 43 which is within 20–70. In addition, the topological polar surface area (TPSA) of the compounds as parameter for the prediction of drug transport properties showed TPSA value greater than 140 Å2 tend to be poor at permeating cell membranes. Despite violation of some rules, approved anticancer and anti-infective drugs from natural products or their semisynthetic derivatives such as taxol and amphotericin B have also some violations but are biologically effective as drugs. Therefore, these results don’t interfere with the development of these compounds as potential SARS-CoV-2 therapeutic agents [43].
Data resulting from Prediction of ADME server (Table S1) revealed that compounds 115, 120, 131, 132, 134 may be better absorbed from the intestinal tract upon oral administration as they showed good human intestinal absorption (HIA) (77.110857, 57.614849, 82.803611, 66.184384 and 79.978975, respectively). Caco-2 cell permeability is a model for selecting drug candidates for oral administration [44]. All compounds showed medium Caco-2 predicted permeability and medium MDCK predicted cell permeability [45]. Moreover, all compounds showed moderate predicted plasma protein binding (PPB) (Table S1) except for compound 159, which showed weak predicted PPB, which indicates predicted decreased excretion and increased predicted half-life. It is important to consider drug’s interaction with plasma proteins, transporters and CYP450s for the successful selection of drug candidate. CYP2C19 and CYP2C9 inhibition leads to increased drug plasma concentration, leading to potential side effects [46,47]. All top scoring compounds were predicted to inhibit CYP2C19 and CYP2C9. CYP2D6 metabolizes many drugs and toxins [48]. None of the top scoring compounds showed predicted inhibitory activity to CYP2D6. CYP3A4 is also involved in metabolism of xenobiotics and is highly expressed in the liver and intestine [49]. The six top scoring compounds were predicted to inhibit CYP3A4. Drug resistance is a major concern in drug development. Multidrug resistance is regulated by a network of ATP-binding cassette (ABC) proteins that detoxify xenobiotics and act as drug transporters and efflux pumps. P glycoprotein (Pgp; ABCB1) is the most popular and well-studied efflux pump [50,51]. Pgp has intrinsic ATPase activity to drive active transport and generate a concentration gradient leading to transport of drugs to the extracellular space and inhibition of drug activity [51]. Pgp is highly expressed in blood-brain barrier cell, liver, intestine and kidney. Thus, it is important to predict drug’s binding to Pgp. Only compound 115 was predicted to inhibit Pgp and hence it may affect the activity or excretion of other Pgp substrates. Compound 159 had the highest water solubility (217.207 mg/L) while the other five compounds had low water solubility, hence this should be considered during drug development. In addition, skin permeability is an important factor to consider during drug development for the potential of dermal drug delivery and risk assessment of drugs that may contact skin [52]. Skin permeability also increases drug’s plasma concentration and activity. It has been reported that logP between −3 to +6 predict drug’s skin permeability [53]. SKlogD, SKlogP and SKlogS are related to drugs’ skin permeability and lipophilicity. All the six top scoring compounds had skin permeability values ranging from −4.6 to −3.6, indicating that they may not be absorbed through skin and thus should not pose a dermal exposure risk. Finally, all compounds did not have predicted ability to pass the blood brain barrier (BBB) and are not expected to be neurotoxic.

3. Materials and Methods

3.1. Phytochemical Review of Genus Aloe

Intensive review of the literatures in ScienceDirect, PubMed, SciFinder and has been conducted to identify compounds from Aloe genus.

3.2. Molecular Docking, Data Software and Visualization

3.2.1. Preparation of Protein and Active Site Prediction

In this study, two SARS-CoV-2 proteins which facilitate viral–host interaction and replication were selected from the RCSB protein databank (https://www.rscb.org/pdb, accessed on 20 February 2021). The proteins are SARS-CoV-2 main protease (PDB ID: 6LU7, resolution = 2.16 Å) [54] and spike glycoprotein (PDB ID: 6M0J, resolution = 2.45 Å) [42]. The 3D protein structures were prepared using Molecular Operating Environment software (MOE 2014.0901) Ligx option. The site finder function used to calculate and predict possible active potential site of selected proteins for ligand binding in the receptor. PyMol 2.3 software was used for visualization.

3.2.2. Preparation of Ligand

Reviewing the available literatures identified 237 phytochemical compounds that were isolated from genus Aloe (Table A1 and Figures S1–S18). All these molecular structures were imported to MOE and subjected to 3D protonation and energy minimization using MMFF94s force field and ligand database was constructed. Ligand coordinate files were extracted from PDB files and used as reference structures for root mean square deviation (RMSD) calculations.

3.2.3. Docking Analysis

Flexible ligand-rigid receptor docking was performed with MOE-DOCK for molecular docking. The parameters of scoring were Triangle Matcher, scoring was set at London dG with 30 output poses and rescoring at GBVI/WSA dG retaining 10 refined poses. The docking score, root mean square deviation (RSMD), 2D and 3D interactions were recorded [55]. The results of docked ligands are chosen based on the most negative docking score. The docking score represents the best-bound ligand conformations and relative binding affinities. The best-docked conformations for comparison of the binding of the drugs and targets of SARS-CoV2 were selected based on number of hydrogen bonds, binding energy (kcal/mol), upper and lower bound RMSD, number of interacting residues, and forces which stabilized the receptor-ligand complex. RMSD and RMSF of the ligand interaction with the target protein were calculated using the following formulas:
R M S D x =   1 N i = 1 n ( r i t ( t x ) )   r i ( t r e f   ) ) 2
R M S F i =   1 T t = 1 T < ( r i t ( t ) )   r i ( t r e f   ) ) 2 >
where N is the number of atoms, tref is the reference time, r′ is the position of the selected atoms in frame x after superimposing on the reference frame, frame x recorded at time tx, T is the trajectory time over which the RMSF was calculated, r is the position of an atom. Poses of docked compounds are automatically calculated by docking function in Molecular Operating Environment software.

3.3. Molecular Dynamics Simulation

MD simulation experiments were performed as previously described [43]. Briefly, the Molecular Dynamics (NAMD) 2.6 software [45], employing the CHARMM27 force field [56] was used for simulations. Hydrogen atoms were added to initial coordinates of proteins using the psfgen plugin included in the Visual Molecular Dynamics (VMD) 1.9 software [57]. Subsequently, the protein systems were solvated using TIP3P water particles and 0.15 M NaCl. The equilibration procedure comprised 1500 minimization steps followed by 30 ps of MD simulation at 10 K with fixed protein atoms. Then, the entire systems were minimized over 1500 steps at 0 K, followed by gradual heating from 10 to 310 K using temperature reassignment during the initial 60 ps of the 100 ps equilibration MD simulation. The final step involved NTP simulation (30 ps) using the Nose–Hoover Langevin piston pressure control at 310 K and 1.013 bars for density (volume) fitting [58]. Thereafter, the MD simulation experiments were continued for 25 ns for the entire systems (20 ns for the enzyme–ligand complexes). The trajectory was stored every 0.1 ns and further analyzed with the VMD 1.9 software. The MD simulation output over 25 ns provided several structural conformers that were sampled every 0.1 ns (250 poses) to evaluate conformational changes of the entire protein structure to analyze the RMSD. All parameters and topologies of the compounds selected for MD simulation were prepared using the online software Ligand Reader & Modeler [59] and the VMD Force Field Toolkit (ffTK) [57]. Binding free energy calculations (∆G) were performed using the free energy perturbation (FEP) method through the web-based software Absolute Ligand Binder [60] together with MD simulation software NAMD 2.6 [45]. Hydrogen bonds and hydrophobic interactions between protein and ligand were also analyzed using Protein-Ligand Interaction Profiler [61].

3.4. Drug Like Properties, and ADME Prediction of the Ligands

The drug likeliness of best pose scoring compounds is specified by the Lipinski’s rule and molecular properties prediction was calculated by the free access website https://www.molinspiration.com/cgi-bin/properties, accessed on 20 February 2021. ADME Prediction were determined by PreADMET estimation server website [62].

4. Conclusions

In recent years, advances in computational resources and software tools led to emergence of molecular dynamics (docking and scoring tool), as the first phase in drug screening and discovery. In addition, absence of new cell culture models for working safely with highly pathogenic viruses makes virtual screening, docking and dynamics of great importance. Aloe genus is a rich source of phytochemicals with a wide range of therapeutic activity. Several natural products from Aloe have shown strong antiviral activity, inhibiting replication and entry and of HSV-1 and 2, human cytomegalovirus (HCMV), influenza A and polio. Aloin significantly reduces replication of oseltamivir-resistant (H1N1) influenza virus. In our study, we applied computational screening of 237 natural product compounds from Aloe genus and identification of six compounds as stable potential inhibitors of SARS-CoV-2 main protease and spike glycoprotein. Our molecular docking analysis showed that theses six compounds are stable and safe. Compounds 132, 134 and 159 were the top three scoring potential inhibitors of for SARS-CoV-2 main protease. These compounds interacted strongly with amino acids in the active site of SARS-CoV-2 main protease. Rutin (154) is known to have antiviral activity against influenza virus [63]. Compounds 115, 120 and compound 131 were the top scoring potential inhibitors of SARS-CoV-2 spike glycoprotein. The results highlighted chromone derivatives as potential inhibitors for SARS-CoV-2 according to their best scores of binding affinities to the mentioned target proteins among the examined compounds. The results of this in-silico investigation (docking and molecular dynamics simulation) should have a great impact for drug repurposing studies. In the future, in vitro, in vivo and clinical studies shall be conducted to further validate the effectiveness of these compounds as potential treatments for COVID-19 and to identify compounds with best pharmacokinetic profiles. In addition, it will be of great importance to apply newly-developed algorithms and utilize the development of steered molecular dynamics for evaluating the binding of the top scoring compounds to SARS-CoV-2 target proteins [64].

Supplementary Materials

The following are available online. Figures S1–S18: Compounds isolated from genus Aloe, Table S1: Predicted pharmacokinetics of top scoring compounds.

Author Contributions

Conceptualization; M.E.A. and M.S.A.A.; Data curation, E.S.E., A.M.S., T.O., L.B. and A.E.E.Z.; Formal analysis, M.E.A., R.A.A., A.M.S., L.B. and A.E.E.Z.; Investigation, M.E.A., H.K.A., R.A.A., E.S.E., A.M.S., A.E.E.Z. and M.S.A.A.; Methodology, M.E.A., H.K.A., R.A.A. and A.M.S.; Project administration, E.S.E., A.E.E.Z. and M.S.A.A.; Resources, E.S.E., T.O., L.B., A.E.E.Z. and M.S.A.A.; Software, M.E.A., H.K.A., R.A.A. and A.M.S.; Supervision, E.S.E., A.E.E.Z. and M.S.A.A.; Validation, L.B., T.O.; Visualization, M.E.A., M.S.A.A. and A.E.E.Z.; Writing—original draft, M.E.A. and H.K.A.; Writing—review & editing, E.S.E., A.E.E.Z. and M.S.A.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare that they have no competing financial interests or personal relationships that could influence the work reported in this paper.

Sample Availability

Samples of the compounds are not available from the authors.

Appendix A

Table A1. Docking score of isolated compounds from Aloe plants against selected SARS CoV-2 proteins.
Table A1. Docking score of isolated compounds from Aloe plants against selected SARS CoV-2 proteins.
No.Name of the CompoundMain Protease (PDB ID: 6LU7)Spike Glycoprotein (PDB ID: 6M0J)Plant SourceRef.
ScoreRSMDScoreRSMD
I-Anthraquinones:
1Chrysophanol−5.45840.8734−5.44470.8709A. pulcherrima, A. dawei, A. megalacantha, A. vera[65,66,67,68]
28-O-Methylchrysophanol−5.06780.4878−4.84060.9352A. dawei[66]
3Aloe-emodin−5.33480.8782−6.05460.9449A. megalacantha, A. arborescens A. vera, A. ferox[67,68,69,70,71,72]
47-Hydroxy-aloe-emodin−5.36601.3842−4.89380.9332A. succotrina[73]
5Nataloe-emodin−5.50170.7836−5.25990.4924A. nyeriensis[74]
6Mono-O-methyl-nataloe-emodin−4.80591.2584−4.61130.9158A. speciosa[74]
7Emodin−4.40600.9758−4.85581.3623A. vera, A. spp[72,75,76,77]
8Saponarin II−4.75201.0190−5.11191.3501A. megalacantha, A. pulcherrima, A. dawei, A. saponaria.[65,66,67,78,79,80]
9Saponarin I−5.56360.8416−5.26251.3226A. megalacantha, A. pulcherrima, A. saponaria.[65,67,78,81]
10Saponarin III−5.59410.9280−5.99070.9256A. megalacantha[67]
11Helminthosporin−5.26980.6476−5.39631.3211A. megalacantha, A. dawei, A. Saponaria[66,67,79]
125-O-Methylziganein−5.60540.8884−3.32090.8923A. hijazensis[82]
13Isoxanthorin−4.77481.2234−3.90510.9124A. Saponaria[79,80]
14Deoxyerythrolaccin−4.79681.3812−5.25721.2640A. ferox, A. Saponaria[71]
15Laccaic acid D Methyl ester−6.16200.9956−5.37860.8172A. Saponaria, A. dawei[66,78]
16Madagascin−4.22121.1677−4.02690.9379A. vera[83]
173-Geranyloxyemodine−5.41611.4882−4.88611.1771A. vera[83]
18Aloetinic acid.−5.60061.6111−2.89681.3504A. vera[84]
19Nataloe-emodin-2-O-β-D-Glucopyranoside−5.39150.9107−4.41161.4211A. nyeriensis[85]
20Aloe-emodin-11-O-rhamnoside−6.55131.1442−6.32331.2130A. vera[72,86]
211,1′,8,8′-Tetrahydroxy-3,3′-dimethyl-4,7′-bianthracene-9,9′,10 (10′H)-trione−6.25190.8628−3.72871.2329A. Saponaria[81]
22Asphodelin−5.84740.9813−6.19661.2760A. megalacantha, A. Saponaria[67,81]
23(1,1′,8,8′,10-Pentahydroxy-3,3′-dimethyl-10,7′-bianthracene-9,9′,10′-trione)−5.25211.3390−6.44191.2682A. Saponaria[81]
2410-(chrysophanol-7′-yl)-10-hydroxychrysophanol-9-anthrone−6.09461.6603−5.74001.2765A. megalacantha[67]
25Chrysalodin−6.38191.4365−5.87871.3680A. megalacantha[67]
2610-O-Methylchrysalodin−5.11800.9839−5.16161.0078A. megalacantha[67]
27Elgonica-dimer A−6.69331.1067−2.58621.4349A. elgonica, A. vera[87,88,89]
28Elgonica-dimer B−6.44161.21425.43781.6013A. elgonica, A. vera[87,88,89]
291,4′,5′,8,9′-Pentahydroxy-2′,6-dimethyl[2,9′-bianthracene]-9,10′(9′H,10H)-dione−6.15951.1344−6.29831.4995A. Saponaria[81]
30Aloin A−6.17581.0225−5.80471.3500A. calidophila, A. schelpei, A. vera, A. perryi, A. ghibensis, A. gilbertii, A. trigonantha[68,69,90,91,92,93,94,95,96]
31Aloin B−7.27941.3710−6.32271.1594A. vera, A. perryi, A. ghibensis, A. gilbertii,A. trigonantha[68,69,90,91,92]
327-Hydroxy-8-O-methylaloin A−6.45551.4506−6.10341.1466A. vera[97]
337-Hydroxy-8-O-methylaloin B−3.35791.4722−5.11351.1263A. vera[97]
34Nataloin−3.12030.9706−7.10481.2243A. nyeriensis[85]
356′-O-Malonylnataloin−6.15411.3288−3.41421.0014A. ellenbeckii[98]
36Homonataloin A−3.46951.0697−5.44431.4948A. lateritia, A. distans, A. cremnophila, A. citrina, A. vera[92,99,100]
37Homonataloin B−4.93491.4914−4.98070.8487A. lateritia, A. excelsa, A. vera, A. perryi[92,99,101,102]
385-Hydroxyaloin A−4.87661.1457−6.61591.1690A. nobilis, A microstigma[99,103,104]
397-Hydroxyaloin A−3.01641.4360−5.55260.9072A. ghibensis, A. succotrina[73,91]
403′-Acetyl-5-hydroxyaloin A−6.67671.3118−5.18321.2296A. nobilis[104]
416′-Acetylglucosyl-5-hydroxyaloin A−5.35351.3456−5.49901.1070A. marlothii,Aloe rupestris[105]
427-Hydroxy (6′-acetylGlucosyl)-aloin−5.26291.2352−6.40081.4160A. succotrina[73,106]
437-hydroxy (6′-acetylglucosyl)-barbloine−4.25651.4076−7.22790.8115A. succotrina[73]
442′,6′-Diacetylglucosyl-5-hydroxyaloin A.−6.71571.5569−5.83131.3241A. nobilis[104]
454′,6′-Diacetylglucosyl-5-hydroxyaloin A−4.98251.7716−5.77461.0670A. nobilis[104]
468-O-Methoxy-7-hydroxyaloin A−5.26940.8441−1.26091.3127A. vera, A. trigonantha[92,93]
474′,6′-O-diacetate-7-Hydroxyaloin A−1.01691.1487−3.66641.1714A. succotrina[73]
484′,6′-O-diacetate-7-Hydroxyaloin B−5.92831.8131−3.81891.1443A. succotrina[73]
49Aloinoside A−5.53271.6795−7.37221.6129A. ferox, A. spp.[92,94,95,107]
50Aloinoside B−4.23862.0931−6.12250.9157A. vera, A. perryi[69,72,86,92,107,108]
51Aloinoside C−6.70351.3302−7.09751.2894A. spp.[69]
52Homonataloside B−7.45480.9307−7.09661.3986A. spp.[109]
53Microdontin A−5.38541.8687−3.05681.3618A. gilbertii, A. microdonta, A. calidophila, A. vera, A. perryi, A. schelpei[92,94,95,96,110]
54Microdontin B−6.44591.4075−6.20211.4041A. vera, A. perryi, A. microdonta[92,110]
55Microstigmin A−6.22121.2076−6.70731.1774A. microstigma, A. broomii[103]
56Desoxyaloin−5.13231.4371−5.45031.1224A. spp.[69]
578-O-Methoxy-7-hydroxyaloin B−4.90671.2985−6.66351.4981A. vera, A. trigonantha[92,93]
587-Hydroxyaloin B−5.71381.6502−6.22450.5540A. ghibensis, A. vera, A. succotrina[73,91]
596′-O-Acetyl-aloin B−6.28291.1810−6.27691.2832A. vera,A. trigonantha[93,111]
606′-O-Acetyl-aloin A−7.01701.3290−6.71741.4200A. vera, A. trigonantha[93,111]
616′-O-Acetyl-10-hydroxyaloin B−5.43051.3206−5.16751.1849A. claviflora[112]
6210-Hydroxyaloin A−5.52351.2517−5.39370.9974A. vera[68,72,86,97]
6310-Hydroxyaloin B−4.52071.1644−6.20371.4416A. vera,A. littoralis[68,72,97,113]
64Aloinoside D−4.02811.5439−5.70671.6872A. sp.[69]
65Deacetyllittoraloin−5.47660.8960−3.28651.2900A. littoralis[114]
66Littoraloin−6.70951.9848−6.64151.4619A. littoralis[114]
67Littoraloin−5.45681.3628−1.64501.3238A. littoralis[113]
68Deacetyllittoraloin−4.65361.4897−3.40451.2748A. littoralis[113]
69Littoraloside−5.13441.6895−5.76991.5576A. littoralis[114]
706,8-Dihydroxy-4-methylbenzanthrone−4.80930.9928−3.60171.1112A. vera[115]
71Anthrone; Enol-form−3.71491.7269−3.32771.7661A. vera[116]
723,4-Dihydro-3,5,7-trihydroxy-9-methyl-1(2H)-anthracenone−5.30770.5534−5.53671.5708A. vera[117]
73Aloesaponol II−5.55960.9169−4.08841.1997A. saponaria[78,80]
74Aloesaponol II-6-methyl ether−5.42780.8931−5.68290.8923A. dawei[66]
75Aloesaponol IV−5.31441.0041−5.87921.2443A. saponaria[79]
76Aloesaponol I−6.27310.6260−5.52681.1253A. megalacantha, A. dawei, A. saponaria[66,67,78,80]
77Aloesaponol III−5.23751.6252−2.82730.5258A. saponaria[79,80]
788-O-Methyl-aloesaponol III−4.75572.6159−5.40510.7606A. saponaria[79]
79Aloesaponol III-8-O-β-D-Glucopyranoside−7.19000.9527−4.93311.1900A. saponaria[79,80]
80Aloesaponol IV-8-O-β-D-Glucopyranoside−6.43290.9875−6.18781.0722A. saponaria[79]
81O-de-Methylaloesaponol IV-4-Epimer, 4-O-β-D-glucopyranoside−5.26211.6257−6.04351.6168A. vera[118,119]
82Aloesaponol IV-4-Epimer, 4-O-β-D-glucopyranoside−5.58781.4697−4.48751.2372A. vera[118,119]
83Aloesaponol II-6-glucoside−5.50971.6405−6.34301.7211A. saponaria[78,80]
84Aloesaponol I-6-O-β-D-Glucoside−6.40181.0686−6.47200.8722A. saponaria[78,80]
85Prechrysophanol−5.00971.3074−4.95030.9273A. graminicola[120]
86Aloechrysone−5.51831.0812−5.22031.1941A. berhana[121]
II-Chromones:
872,7-Dihydroxy-5-methylchromone−4.10441.5749−4.33802.1225A. arborescens[122]
88Altechromone A−4.39641.3753−4.01772.5253A. vera, A. ferox[76]
897-Hydroxy-5-(hydroxymethyl)-2-methylchromone−4.66541.0133−4.00341.0116A. vera, A. spp.[68,69,123]
905-(Hydroxymethyl)-7-methoxy-2-methylchromone−5.09421.7087−4.36040.6077A. vera[68]
91Aloesone−5.34340.8671−4.44591.6654A. spp.[124]
92Aloesol−5.48241.1217−4.51191.3491A. spp.[124]
93Saikochromone A−4.02461.1592−4.66750.8802A. vera[68]
942−Carboxyethenyl-5,7-dihydroxychromone−5.45600.8633−4.30331.4776A. cremnophila[125]
955-((4E)-2′-Oxopentenyl)-2-hydroxymethylchromone−5.24081.2458−4.65430.9525A. vera[68]
962-Acetonyl-7-hydroxy-8-(3-hydroxyacetonyl)-5-methylchromone−6.12951.1673−5.62940.8828A. ferox[126]
975-((S)-2′-Oxo-4′-hydroxypentyl)-2-hydroxymethylchromone−5.65330.8672−4.57421.0536A. spp., A. vera[68,69,123]
987-O-methyl-(R)-aloesinol−5.77491.3751−6.33491.0728A. capensis, A. rubroviolacea, A. spicata[69,123,127,128]
997-O-methyl-(S)-aloesinol−4.98621.0634−6.58261.3861A. vera[129,130,131]
100Aloesinol; (2′S)-form−6.54121.3314−5.92310.9830A. vera[129,130,131]
1018-C-Glucosyl-7-O-methylaloediol−6.56161.4702−6.19831.3370A. vera[130,131,132]
102C-2′-Decoumaroyl-aloeresin G−6.32032.2082−5.66740.6253A. vera, A. spp.[69,123,132]
103Deacetylaloesin−5.45801.4439−5.67690.7237A. vera var. chinensis[133]
104Isobiflorin−4.39171.0141−5.53611.3146A. vera[131]
1058-C-Glucosyl-noreugenin−6.01550.9498−5.33111.2871A. vera[132]
106Neoaloesin A−4.13141.4402−5.37271.2553A. vera[132,134]
107Aloesin (Aloeresin B)−5.19201.4721−6.00110.9641A. vera, A. monticola, A. trigonantha A. saponaria, A. arborescens[68,92,93,132,135,136,137]
1087-O-Methylaloesin−6.52111.0562−4.78730.7405A. vera, A. rupestris[92,105]
1092-Acetonyl-8-(2-furoylmethyl)-7-hydroxy-5-methylchromone−5.57851.2539−5.94861.4100A. ferox[126]
1108-C-Glucosyl-(R)-aloesol−6.98101.2199−6.58851.2379A. vera[68,92,132,136]
1118-C-Glucosyl-(S)-aloesol−5.99791.2224−3.01781.2511A. vera[132]
1128-C-Glucosyl-7-methoxy-(R)-aloesol−6.91221.5762−5.22510.8301A. vera[132]
1138-C-Glucosyl-7-methoxy-(S)-aloesol−6.28671.2856−4.99030.9918A. vera[68,132,136]
1142″-O-p-Coumaroyl-(S)-aloesinol−6.49491.1099−6.38801.7012A. nobilis[138]
1152″-O-(4-methoxycinnamoyl)-(S)-aloesinol−6.82220.9924−8.05741.4996A. nobilis[104]
1162″-O-Cinnamoyl-8-C-glucosyl-7-O-methyl-aloediol A−6.98420.9200−7.44911.1903A. vera[131]
1174″-Deoxyaloeresin D−1.96591.4660−7.39691.4416A. vera[139]
118Aloeresin D−7.38651.5512−6.68081.3432A. ferox, A. vera[69,123,132,140,141]
119Isoaloeresin D−5.82801.4838−7.55611.2248A. vera[68,92,129,130,131,132,136]
120Rabaichromone−7.39141.9603−7.87151.1703A. rabaiensis, A. vera[132,142]
121Isorabaichromone−7.54701.4108−7.26561.4428A. vera[129,130,131,132]
122Aloeresin J−6.97231.9908−6.43521.4979A. vera[132]
1238-[C-β-D-[2-O-(E)-cinnamoyl] glucopyranosyl]-2-[(R)-2-hydroxypropyl]7-methoxy-5-methylchromone−7.35160.9022−5.89351.2315A. vera[132]
124Aloeresin K−6.28571.5121−6.90201.8499A. vera[111,132]
125Aloeresin F−7.58141.9617−6.63201.4268A. peglerae[143,144]
1262″-O-E-cinnamoyl-2′-Ketone (2′R) aloesinol-7-methyl ether−6.52271.6510−5.79261.2582A. broomii[145]
1277-O-Methylaloeresin A−7.30611.1520−6.84820.9873A. vera, A. perryi, A. marlothii[92,105,132]
128Aloeresin A−7.02281.1297−6.17301.3563A. cremnophila, A. jacksonii, A. arborescens, A. vera[92,132,146]
129Isoaloeresin A−7.06981.8469−6.96011.6099A. ferox[147]
1306″-O-p-Coumaroylaloesin−6.10541.3654−5.55471.4791A. vera, A. castanea[143,144]
131Aloeribide−6.59881.1683−8.08861.1338A. vera, A. monticola[86]
1322″-O-(3,4-dihydroxy-E-cinnamoyl) 2′-Ketone, (2′R) Aloesinol-7-methyl ether−7.68301.3783−5.18570.8784A. broomii[145]
1332″-O-Feruloylaloesin−6.32461.4172−7.35711.3775A. arborescens[148]
1342″-O-(4-hydroxy-3-methoxy-E-cinnamoyl) 2′-Ketone, (2′R)-aloesinol-7-methyl ether−7.95101.7236−6.50961.2159A. africana[145]
1352″-O-tigloyl-2′-Ketone, (2′R)-aloesinol−7.00601.0927−1.48651.2366A. cremnophila, A. jacksonii[125]
1362″,6″-bis-O-(4-hydroxy-E-cinnamoyl) 2′-Ketone, (2′R)-aloesinol−7.64661.6059−6.85892.8077A. speciosa[145]
137Aloeresin E−6.84831.7166−3.39351.5214A. peglerae[143,144]
138Aloeresin C−3.19261.40364.65451.2671A. spp.[140,141]
139trans-p-coumaroyl-4′-O-Glucosyl-isoaloeresin D-I a−7.30701.8254−5.25101.8902A. vera[129,130,131,132]
140cis-p-coumaroyl-4′-O-Glucosylisoaloeresin D-II a−5.12141.9414−6.98041.5954A. vera[129,130,131,132]
1418-C-(2′-O-coumaroylglucosyl)-7-hydroxy-5-methyl-chromone-2-carboxylic acid−7.00931.2428−3.86211.2324A. vera, A. perryi[92]
1429-Dihydroxyl-2′-O-(Z)-cinnamoyl-7-methoxy-aloesin−7.66161.3432−6.55101.0865A. vera[132]
1435-((S)-2′-Oxo-4′-hydroxypentyl)-2-(β-D-glucopyranosyloxy-methyl) chromone−6.36041.6246−5.66220.9999A. vera[68,123]
1448-Glucosyl-(2′-O-cinnamoyl)-7-O-methyl-aloediol A−6.84591.3641−6.16061.1397A. vera[132]
145Aloeresin G−5.61361.0175−7.20211.2832A. vera[69,123,149]
146Aloeresin H−7.57381.6723−7.55251.2922A. ferox[150]
147Aloeresin I−7.43962.2729−7.43371.7717A. ferox[151]
148Iso-aloesin−6.58170.9104−5.32661.3373A. vera var. chinensis[152]
149Aloeveraside B−7.06581.3476−5.61131.4352A. vera[72,132]
150Aloeveraside A−7.13441.1636−4.83531.2241A. vera[72,132]
151Furoaloesone−5.27651.3977−5.30161.0583A. ferox[153]
III-Coumarins:
152Coumarin−4.35651.3443−3.51420.7049A. vera[86]
1537-Demethylsiderin−4.48571.2535−4.66411.3044A. vera, A. megalacantha[72,86]
IV-Flavonoids:
154Apigenin−4.88321.1056−4.94311.4340A. vera[154]
155Kaempferol−5.56220.6388−5.27311.3518A. vera[154]
156Quercetin−5.50450.8816−5.00911.3859A. vera[154]
157Myricetin−4.34581.1650−5.18721.7167A. vera[154]
158Quercitrin−5.98301.1253−6.67871.2189A. vera[154]
159Rutin−7.72851.4051−6.51751.7434A. vera[154]
1603′,5′,6,7-Tetra-methyl ether, 5-O-[α-L-rhamnopyranosyl-(1→6)-β-D-galactopyranoside], 3′,4′,5,5′,6,7-Hexahydroxyflavone−6.46801.4828−7.24751.4030A. vera[155]
161Catechin−5.52320.7426−5.07220.8989A. vera[154]
162Epicatechin−5.34311.3310−5.08291.3469A. vera[154]
1633′,4′,6-Tri-Me ether, 5-O-[α-L-rhamnopyranosyl-(1→6)-β-D-glucopyranoside], 3′,4′,5,6,7-Pentahydroxyisoflavone0.33900.9914−4.02571.8089A. vera[156]
V-Phenolic Compounds:
a-Simple phenolic compounds:
164Pyrocatechol−3.90370.5358−3.19881.0732A. ferox[157]
165Salicylaldehyde−2.67541.0693−3.21781.2797A. vera[72,86]
166p-Cresol−3.77801.1628−3.51101.1472A. vera[72,86]
167p-Hydroxyacetophenone−4.02241.8629−3.77730.7970A. ferox[157]
168p-Hydroxybenzaldehyde−3.10970.7939−3.54401.5859A. ferox[157]
169p-Anisaldehyde−3.91641.3669−3.44651.1778A. vera[72,86]
170Phloretic acid−4.38231.2905−4.11071.0512A. vera[72,86]
171Methyl-3-(4-hydroxyphenyl) propionate−4.74960.8674−3.91360.9518A. vera[72,86]
172Coumaric acid−4.50551.4107−3.65871.1690A. vera[154]
173Caffeic acid−4.19691.2736−3.88591.0162A. vera[154]
174Ferulic acid−4.38450.9040−4.46891.0173A. vera[154]
175Sinapic acid−5.10211.0097−4.27631.1207A. vera[154]
176Orcinol−3.93621.2430−3.76541.4027A. sp.[69]
177Gentisic acid−3.85951.6010−3.98281.3232A. vera[154]
178Protocatechuic acid−3.64471.5645−4.10491.3881A. vera[154]
179Vanillic acid−4.02561.0470−4.08270.9452A. vera[154]
180Gallic acid−4.19610.5757−3.56790.8059A. vera[154]
181Syringic acid−4.14671.2874−4.63081.2947A. vera[154]
1821-(2,4-Dihydroxy-6-methylphenyl) ethanone−3.81831.3698−3.93360.9642A. vera[72]
b-Phenyl-Pyran and Phenyl-pyrone derivatives:
183Aloenin aglycone.−4.84941.2878−3.89571.3574A. vera, A. spp.[68,69,123,158,159]
1846-(2,4-Dihydroxy-6-pentylphenyl)-4-hydroxy-2H-pyran-2-one−3.17491.0339−5.54281.0487A. arborescens[160]
1856-[[2,4-Dihydroxy-6-[2-(4-hydroxyphenyl) ethenyl]-4-hydroxy-2H-pyran-2-one−5.97421.4002−5.88521.4624A. arborescens[160]
186Aloenin A−5.37231.3544−4.56451.5586A. vera, A. arborescens[85,92,123,159]
1872″-O-trans-p-coumaroyl-aloenin−7.03281.6452−6.10261.2183A. vera, A. nyeriensis, A. spicata[159]
1884″,6″-Ethylidenealoenin−6.32212.1206−6.85801.0640A. arborescens, A. hijazensis[82]
189Aloenin C−4.05351.3070−4.77011.7252A. sp.[123]
19010-O-β-D-glucopyranosyl-aloenin−6.13781.0013−6.69761.3721A. sp., A. vera, A. spicata[127,136]
191Aloenin B−7.05731.9426−4.77661.3631A. hijazensis, A. spicata, A. vera[68,69,82,123,136]
1922,4-Dihydroxy-β-(4-hydroxyphenyl)-5-(4-methoxy-2-oxo-2H-pyran-6-yl)-6-methylbenzenepropanoic acid; Et ester−6.05581.2415−6.57471.4078A. vera[161]
1936−[[3,5-Dihydroxy-2-(1-oxohexyl) phenyl] methyl]-4-hydroxy-2H-pyran-2-one−5.51751.3328−5.19891.4473A. arborescens[160]
1946-[[3,5-Dihydroxy-2-[3-(4-hydroxyphenyl)-1-oxo-2-propenyl] phenyl] methyl]-4-hydroxy-2H-pyran-2-one−3.95371.4667−6.80471.2550A. arborescens[160]
1953,4-Dihydro-7-hydroxy-4-(4-hydroxyphenyl)-6-(4-methoxy-2-oxo-2H-pyran-6-yl)-5-methyl-2H-1-benzopyran-2-one−6.00221.3606−6.44751.1932A. vera[161]
196Feralolide−5.79801.8453−5.47811.4521A. arborescent, A. ferox, A. vera[70,82,86,162,163]
1975′-O-Methylferalolide−6.38081.3485−5.99901.3806A. vera[164]
198Feralolide-3′-O-β-D-Glucopyranoside−5.63641.3200−4.30412.3081A. arborescens, A. vera[70,165]
1993,3′-Bi(3,4-dihydro-6-methoxy-2H-1-benzopyran-4-ol). 3,3′-Bi(3,4-dihydro-4-hydroxy-6-methoxy-2H-1-benzopyran)−5.30931.7270−5.20611.4071A. vera[166]
c-Benzofurans:
2005-Hydroxy-3-methylnaphtho[2,3-c] furan-4(1H)-one−4.33831.1118−4.09580.7359A. ferox[71]
2015-Hydroxy-3-methylnaphtho[2,3-c] furan-4(9H)-one−4.13461.3867−4.09570.7362A. ferox[71]
202Isoeleutherol−4.75441.2527−4.52320.7749A. graminicola[167]
203Isoeleutherol glucoside−6.00441.3381−4.63251.4751A. saponaria[80]
2048-Hydroxy-1-methylnaphtho[2,3-c] furan-4,9-dione−5.04471.3119−4.32710.8506A. ferox[71]
d-Naphthalin derivatives:
205Droserone−4.55251.0154−4.14501.1420A. dawei[66]
206Droserone-5-methyl ether−4.88811.1213−4.28591.2524A. dawei[66]
207Hydroxydroserone−4.97121.0739−3.76561.9742A. dawei[66]
208Ancistroquinone C−4.88941.4452−4.21080.5425A. dawei[66]
2095,8-Dihydroxy-3-methoxy-2-methyl-1,4-naphthoquinone−5.05911.0203−3.59060.8949A. dawei[66]
210Malvone A−4.70030.6291−4.71581.0717A. dawei[66]
2116-Hydroxy-3,5-dimethoxy-2-methyl-1,4-naphthoquinone−4.86031.1151−4.21081.2529A. dawei[66]
2121,8-Dimethoxynepodinol−4.86251.0642−4.14780.9217A. megalacantha[67]
2133-Hydroxy-1-(1,7-dihydroxy-3,6-dimethoxynaphthalen-2-yl) propan-1-one−5.85040.7746−5.23761.0515A. vera[168]
214Plicataloside−3.41181.4389−5.67931.2349A. plicatilis[169]
215Kenyaloside−7.29101.6098−4.69541.3254a Kenyan A. spp.[170]
216Aloveroside A−7.58021.8561−6.86121.5133A. vera, A. spp.[69,123,171]
217Isoeleutherin−4.32171.3635−3.72340.6554A. graminicola[167]
2181-(4-Hydroxyphenyl)-6,9-dihydroxy-7-methyl-8-acetyl-1,2-dihydro-(3H)-naphtho[2,1-b] pyran-3-one−3.81171.5671−6.92880.9758A. ferox[172]
VI-Alkaloids:
2194,7-Dichloroquinoline−1.69760.2458−3.54151.2153A. hijazensis[82]
220N, N-Dimethyl-(+)-coniine−1.94200.6273−1.70661.3831A. sabaea[173]
221γ-Coniceine−4.23580.9461−3.68230.9519A. sp.[173,174]
VII-Fatty acid derivatives:
22210-Hydroxyoctadecanoic acid−6.14271.4991−4.46571.1518A. ferox[157]
22310-Oxooctadecanoic acid-6.33621.0628−5.80791.1759A. ferox[157]
224Methyl-26-O-feruloyl-oxyhexacosanate−7.60531.7939−7.23491.2702A. megalacantha[67]
VIII-Miscellaneous compounds:
225Nilic acid−4.00960.9979−2.39510.6370A. littoralis[113]
226N-(4-Chlorobutyl) butanamide−4.75960.6115−3.97580.9418A. sabaea[173]
2273-Furanmethanol−3.91770.9518−3.49390.7289A. arborescens[175]
2283, 6-Dioxo-3, 3a, 6, 6 a-tetrahydropyrrolo [3, 4-c] pyrrole-1, 4-dicarboxamide−4.10121.7165−3.80361.0232A. vera[176]
2291-(2,4-Dihydroxy-6-methylphenyl)-1-(4-hydroxyphenyl) ethane−4.77352.0838−2.09691.0198A. ferox[177]
230Chlorogenic acid−4.17221.3907−5.67611.2173A. vera, A. arborescens[154]
231Pluridone; (E)-form−4.97121.2855−3.81811.0980A. pluridens[178]
232Feroxidin−4.31491.0601−3.99581.0673A. ferox, A. arborescens[69,70,72,86,123,179,180]
233Feroxin A−5.53131.4467−5.28691.0687A. spp.[181]
234Feroxin B−6.64462.1091−6.29021.7466A. spp.[181]
235Veracylglucan A−5.78291.2199−4.76530.8308A. vera[182]
236Veracylglucan B−3.08361.0444−1.03681.0780A. vera[182]
237Veracylglucan C−6.34332.2092−7.25121.9263A. vera[182]

References

  1. EI Zowalaty, M.E.; Young, S.G.; Järhult, J.D. Environmental impact of the COVID-19 pandemic–a lesson for the future. Infect. Ecol. Epidemiol. 2020, 10, 1768023. [Google Scholar] [CrossRef]
  2. El Zowalaty, M.E.; Järhult, J.D. From SARS to COVID-19: A previously unknown SARS-related coronavirus (SARS-CoV-2) of pandemic potential infecting humans–Call for a One Health approach. One Health 2020, 9, 100124. [Google Scholar] [CrossRef]
  3. Wu, F.; Wang, A.; Liu, M.; Wang, Q.; Chen, J.; Xia, S.; Ling, Y.; Zhang, Y.; Xun, J.; Lu, L.; et al. Neutralizing antibody responses to SARS-CoV-2 in a COVID-19 recovered patient cohort and their implications. 2020. Available online: https://papers.ssrn.com/sol3/papers.cfm?abstract_id=3566211 (accessed on 15 January 2021).
  4. Li, F. Structure, Function, and Evolution of Coronavirus Spike Proteins. Annu. Rev. Virol. 2016, 3, 237–261. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Perlman, S.; Netland, J. Coronaviruses post-SARS: Update on replication and pathogenesis. Nat. Rev. Genet. 2009, 7, 439–450. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Li, F. Receptor Recognition Mechanisms of Coronaviruses: A Decade of Structural Studies. J. Virol. 2015, 89, 1954–1964. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Li, W.; Moore, M.J.; Vasilieva, N.; Sui, J.; Wong, S.K.; Berne, M.A.; Somasundaran, M.; Sullivan, J.L.; Luzuriaga, K.; Greenough, T.C.; et al. Angiotensin-converting enzyme 2 is a functional receptor for the SARS coronavirus. Nature 2003, 426, 450–454. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Boopathi, S.; Poma, A.B.; Kolandaivel, P. Novel 2019 coronavirus structure, mechanism of action, antiviral drug promises and rule out against its treatment. J. Biomol. Struct. Dyn. 2020, 1–10. [Google Scholar] [CrossRef] [Green Version]
  9. Heald-Sargent, T.; Gallagher, T. Ready, Set, Fuse! The Coronavirus Spike Protein and Acquisition of Fusion Competence. Viruses 2012, 4, 557–580. [Google Scholar] [CrossRef] [Green Version]
  10. Moreira, R.A.; Chwastyk, M.; Baker, J.L.; Guzman, H.A.V.; Poma, A.B. Quantitative determination of mechanical stability in the novel coronavirus spike protein. Nanoscale 2020, 12, 16409–16413. [Google Scholar] [CrossRef]
  11. Moreira, R.A.; Guzman, H.V.; Boopathi, S.; Baker, J.L.; Poma, A.B. Characterization of Structural and Energetic Differences between Conformations of the SARS-CoV-2 Spike Protein. Materials 2020, 13, 5362. [Google Scholar] [CrossRef]
  12. Anand, K. Coronavirus Main Proteinase (3CLpro) Structure: Basis for Design of Anti-SARS Drugs. Science 2003, 300, 1763–1767. [Google Scholar] [CrossRef] [Green Version]
  13. Gil, C.; Ginex, T.; Maestro, I.; Nozal, V.; Barrado-Gil, L.; Cuesta-Geijo, M.Á.; Urquiza, J.; Ramírez, D.; Alonso, C.; Campillo, N.E.; et al. COVID-19: Drug Targets and Potential Treatments. J. Med. Chem. 2020, 63, 12359–12386. [Google Scholar] [CrossRef] [PubMed]
  14. Heimfarth, L.; Serafini, M.R.; Martins-Filho, P.R.; Quintans, J.D.S.S.; Quintans-Júnior, L.J. Drug repurposing and cytokine management in response to COVID-19: A review. Int. Immunopharmacol. 2020, 88, 106947. [Google Scholar] [CrossRef] [PubMed]
  15. Chakraborty, R.; Parvez, S. COVID-19: An overview of the current pharmacological interventions, vaccines, and clinical trials. Biochem. Pharmacol. 2020, 180, 114184. [Google Scholar] [CrossRef]
  16. Thota, S.M.; Balan, V.; Sivaramakrishnan, V. Natural products as home-based prophylactic and symptom management agents in the setting of COVID-19. Phytother. Res. 2020, 34, 3148–3167. [Google Scholar] [CrossRef] [PubMed]
  17. Ninfali, P.; Antonelli, A.; Magnani, M.; Scarpa, E.S. Antiviral Properties of Flavonoids and Delivery Strategies. Nutrient 2020, 12, 2534. [Google Scholar] [CrossRef] [PubMed]
  18. Sui, Z.; Salto, R.; Li, J.; Craik, C.; de Montellano, P.R.O. Inhibition of the HIV-1 and HIV-2 proteases by curcumin and curcumin boron complexes. Bioorg. Med. Chem. 1993, 1, 415–422. [Google Scholar] [CrossRef]
  19. Mehla, R.; Bivalkar-Mehla, S.; Chauhan, A. A Flavonoid, Luteolin, Cripples HIV-1 by Abrogation of Tat Function. PLoS ONE 2011, 6, e27915. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Ürményi, F.G.G.; Saraiva, G.D.N.; Casanova, L.M.; Matos, A.D.S.; Camargo, L.M.D.M.; Romanos, M.T.V.; Costa, S.S. Anti-HSV-1 and HSV-2 Flavonoids and a New Kaempferol Triglycoside from the Medicinal Plant Kalanchoe daigremontiana. Chem. Biodivers. 2016, 13, 1707–1714. [Google Scholar] [CrossRef]
  21. Behbahani, M.; Sayedipour, S.; Pourazar, A.; Shanehsazzadeh, M. In vitro anti-HIV-1 activities of kaempferol and kaempferol-7-O-glucoside isolated from Securigera securidaca. Res. Pharm. Sci. 2015, 9, 463–469. [Google Scholar]
  22. Mpiana, P.T.; Tshibangu, D.S.; Kilembe, J.T.; Gbolo, B.Z.; Mwanangombo, D.T.; Inkoto, C.L.; Lengbiye, E.M.; Mbadiko, C.M.; Matondo, A.; Bongo, G.N.; et al. Identification of potential inhibitors of SARS-CoV-2 main protease from Aloe vera compounds: A molecular docking study. Chem. Phys. Lett. 2020, 754, 137751. [Google Scholar] [CrossRef]
  23. Rezazadeh, F.; Moshaverinia, M.; Motamedifar, M.; Alyaseri, M. Assessment of Anti HSV-1 Activity of Aloe Vera Gel Extract: An In Vitro Study. J. Dent. 2016, 17, 49–54. [Google Scholar]
  24. Abd-Alla, H.I.; Abu-Gabal, N.S.; Hassan, A.Z.; El-Safty, M.M.; Shalaby, N.M.M. Antiviral activity of Aloe hijazensis against some haemagglutinating viruses infection and its phytoconstituents. Arch. Pharmacal Res. 2012, 35, 1347–1354. [Google Scholar] [CrossRef]
  25. Saoo, K.; Miki, H.; Ohmori, M.; Winters, W.D. Antiviral Activity of Aloe Extracts against Cytomegalovirus. Phytother. Res. 1996, 10, 348–350. [Google Scholar] [CrossRef]
  26. Huang, C.-T.; Hung, C.-Y.; Hseih, Y.-C.; Chang, C.-S.; Velu, A.B.; He, Y.-C.; Huang, Y.-L.; Chen, T.-A.; Chen, T.-C.; Lin, C.-Y.; et al. Effect of aloin on viral neuraminidase and hemagglutinin-specific T cell immunity in acute influenza. Phytomedicine 2019, 64, 152904. [Google Scholar] [CrossRef] [PubMed]
  27. Subbaiyan, A.; Ravichandran, K.; Singh, S.V.; Sankar, M.; Thomas, P.; Dhama, K.; Malik, Y.S.; Singh, R.K.; Chaudhuri, P. In silico Molecular Docking Analysis Targeting SARS-CoV-2 Spike Protein and Selected Herbal Constituents. J. Pure Appl. Microbiol. 2020, 14, 989–998. [Google Scholar] [CrossRef]
  28. Jorgensen, W.L. The Many Roles of Computation in Drug Discovery. Science 2004, 303, 1813–1818. [Google Scholar] [CrossRef] [PubMed]
  29. Ramajayam, R.; Tan, K.-P.; Liang, P.-H. Recent development of 3C and 3CL protease inhibitors for anti-coronavirus and anti-picornavirus drug discovery. Biochem. Soc. Trans. 2011, 39, 1371–1375. [Google Scholar] [CrossRef]
  30. Salman, S.; Shah, F.H.; Idrees, J.; Idrees, F.; Velagala, S.; Ali, J.; Khan, A.A. Virtual screening of immunomodulatory medicinal compounds as promising anti-SARS-COV-2 inhibitors. Future Virol. 2020, 15, 267–275. [Google Scholar] [CrossRef]
  31. Simmons, G.; Gosalia, D.N.; Rennekamp, A.J.; Reeves, J.D.; Diamond, S.L.; Bates, P. Inhibitors of cathepsin L prevent severe acute respiratory syndrome coronavirus entry. Proc. Natl. Acad. Sci. USA 2005, 102, 11876–11881. [Google Scholar] [CrossRef] [Green Version]
  32. Al-Bari, M.A.A. Chloroquine analogues in drug discovery: New directions of uses, mechanisms of actions and toxic manifestations from malaria to multifarious diseases. J. Antimicrob. Chemother. 2015, 70, 1608–1621. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Yang, N.; Shen, H.-M. Targeting the Endocytic Pathway and Autophagy Process as a Novel Therapeutic Strategy in COVID-19. Int. J. Biol. Sci. 2020, 16, 1724–1731. [Google Scholar] [CrossRef]
  34. Rivinoja, A.; Hassinen, A.; Kokkonen, N.; Kauppila, A.; Kellokumpu, S. Elevated Golgi pH impairs terminalN-glycosylation by inducing mislocalization of Golgi glycosyltransferases. J. Cell. Physiol. 2009, 220, 144–154. [Google Scholar] [CrossRef]
  35. Li, Z.; Li, X.; Huang, Y.Y.; Wu, Y.; Liu, R.; Zhou, L.; Lin, Y.; Wu, D.; Zhang, L.; Liu, H.; et al. Identify potent SARS-CoV-2 main protease inhibitors via accelerated free energy perturbation-based virtual screening of existing drugs. Proc. Natl. Acad. Sci. USA 2020, 117, 27381–27387. [Google Scholar] [CrossRef] [PubMed]
  36. Velázquez-Libera, J.L.; Durán-Verdugo, F.; Valdés-Jiménez, A.; Núñez-Vivanco, G.; Caballero, J. LigRMSD: A web server for automatic structure matching and RMSD calculations among identical and similar compounds in protein-ligand docking. Bioinformatics 2020, 36, 2912–2914. [Google Scholar] [CrossRef]
  37. Totrov, M.; Abagyan, R. Flexible ligand docking to multiple receptor conformations: A practical alternative. Curr. Opin. Struct. Biol. 2008, 18, 178–184. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. De Vivo, M.; Masetti, M.; Bottegoni, G.; Cavalli, A. Role of Molecular Dynamics and Related Methods in Drug Discovery. J. Med. Chem. 2016, 59, 4035–4061. [Google Scholar] [CrossRef] [PubMed]
  39. De Vivo, M.; Cavalli, A. Recent advances in dynamic docking for drug discovery. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2017, 7, 1320. [Google Scholar] [CrossRef]
  40. Souza, P.C.T.; Thallmair, S.; Conflitti, P.; Ramírez-Palacios, C.; Alessandri, R.; Raniolo, S.; Limongelli, V.; Marrink, S.J. Protein–ligand binding with the coarse-grained Martini model. Nat. Commun. 2020, 11, 1–11. [Google Scholar] [CrossRef] [PubMed]
  41. Arokiyaraj, S.; Stalin, A.; Kannan, B.S.; Shin, H. Geranii Herba as a Potential Inhibitor of SARS-CoV-2 Main 3CLpro, Spike RBD, and Regulation of Unfolded Protein Response: An In Silico Approach. Antibiotics 2020, 9, 863. [Google Scholar] [CrossRef]
  42. Hussien, M.A.; Abdelaziz, A.E. Molecular docking suggests repurposing of brincidofovir as a potential drug targeting SARS-CoV-2 ACE2 receptor and main protease. Netw. Model. Anal. Health Inform. Bioinformatics 2020, 9, 1–18. [Google Scholar]
  43. Sayed, A.M.; Alhadrami, H.A.; El-Gendy, A.O.; Shamikh, Y.I.; Belbahri, L.; Hassan, H.M.; Abdelmohsen, U.R.; Rateb, M.E. Microbial natural products as potential inhibitors of SARS-CoV-2 main protease (Mpro). Microorganisms 2020, 8, 970. [Google Scholar] [CrossRef] [PubMed]
  44. Saurav, K.; Zhang, W.; Saha, S.; Zhang, H.; Li, S.; Zhang, Q.; Wu, Z.; Zhang, G.; Zhu, Y.; Verma, G. In silico molecular docking, preclinical evaluation of spiroindimicins A-D, lynamicin A and D isolated from deep marine sea derived Streptomyces sp. SCSIO 03032. Interdiscip. Sci. Comput. Life Sci. 2014, 6, 187–196. [Google Scholar] [CrossRef]
  45. Phillips, J.C.; Braun, R.; Wang, W.; Gumbart, J.; Tajkhorshid, E.; Villa, E.; Chipot, C.; Skeel, R.D.; Kalé, L.; Schulten, K. Scalable molecular dynamics with NAMD. J. Comput. Chem. 2005, 26, 1781–1802. [Google Scholar] [CrossRef] [Green Version]
  46. Foti, R.S.; Wahlstrom, J.L. CYP2C19 Inhibition: The Impact of Substrate Probe Selection on in Vitro Inhibition Profiles. Drug Metab. Dispos. 2007, 36, 523–528. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Van Booven, D.; Marsh, S.; McLeod, H.; Carrillo, M.W.; Sangkuhl, K.; Klein, T.E.; Altman, R.B. Cytochrome P450 2C9-CYP2C9. Pharm. Genom. 2010, 20, 277–281. [Google Scholar] [CrossRef] [PubMed]
  48. Bertilsson, L.; Dahl, M.-L.; Dalén, P.; Al-Shurbaji, A. Molecular genetics of CYP2D6: Clinical relevance with focus on psychotropic drugs. Br. J. Clin. Pharmacol. 2002, 53, 111–122. [Google Scholar] [CrossRef] [Green Version]
  49. Dai, D.; Tang, J.; Rose, R.; Hodgson, E.; Bienstock, R.J.; Mohrenweiser, H.W.; A Goldstein, J. Identification of variants of CYP3A4 and characterization of their abilities to metabolize testosterone and chlorpyrifos. J. Pharmacol. Exp. Ther. 2001, 299, 825–831. [Google Scholar]
  50. Wessler, J.D.; Grip, L.T.; Mendell, J.; Giugliano, R.P. The P-Glycoprotein Transport System and Cardiovascular Drugs. J. Am. Coll. Cardiol. 2013, 61, 2495–2502. [Google Scholar] [CrossRef] [Green Version]
  51. Sharom, F.J. Complex Interplay between the P-Glycoprotein Multidrug Efflux Pump and the Membrane: Its Role in Modulating Protein Function. Front. Oncol. 2014, 4, 41. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Singh, S.; Singh, J. Transdermal drug delivery by passive diffusion and iontophoresis: A review. Med. Res. Rev. 1993, 13, 569–621. [Google Scholar] [CrossRef] [PubMed]
  53. Mitragotri, S.; Anissimov, Y.G.; Bunge, A.L.; Frasch, H.F.; Guy, R.H.; Hadgraft, J.; Kasting, G.B.; Lane, M.E.; Roberts, M.S. Mathematical models of skin permeability: An overview. Int. J. Pharm. 2011, 418, 115–129. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Jin, Z.; Du, X.; Xu, Y.; Deng, Y.; Liu, M.; Zhao, Y.; Zhang, B.; Li, X.; Zhang, L.; Peng, C.; et al. Structure of M pro from SARS-CoV-2 and discovery of its inhibitors. Nature 2020, 582, 289–293. [Google Scholar] [CrossRef] [Green Version]
  55. Cozza, G.; Moro, S. Medicinal Chemistry and the Molecular Operating Environment (MOE): Application of QSAR and Molecular Docking to Drug Discovery. Curr. Top. Med. Chem. 2008, 8, 1555–1572. [Google Scholar] [CrossRef]
  56. MacKerell, A.D., Jr.; Bashford, D.; Bellott, M.L.D.R.; Dunbrack, R.L., Jr.; Evanseck, J.D.; Field, M.J.; Fischer, S.; Gao, J.; Guo, H.; Ha, S.; et al. All-atom empirical potential for molecular modeling and dynamics studies of proteins. J. Phys. Chem. B 1998, 102, 3586–3616. [Google Scholar] [CrossRef] [PubMed]
  57. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  58. Martyna, G.J.; Tobias, D.J.; Klein, M.L. Constant pressure molecular dynamics algorithms. J. Chem. Phys. 1994, 101, 4177–4189. [Google Scholar] [CrossRef]
  59. Jo, S.; Kim, T.; Iyer, V.G.; Im, W. CHARMM-GUI: A web-based graphical user interface for CHARMM. J. Comput. Chem. 2008, 29, 1859–1865. [Google Scholar] [CrossRef]
  60. Jo, S.; Jiang, W.; Lee, H.S.; Roux, B.; Im, W. CHARMM-GUI Ligand Binder for Absolute Binding Free Energy Calculations and Its Application. J. Chem. Inf. Model. 2013, 53, 267–277. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Salentin, S.; Schreiber, S.; Haupt, V.J.; Adasme, M.F.; Schroeder, M. PLIP: Fully automated protein–ligand interaction profiler. Nucleic Acids Res. 2015, 43, W443–W447. [Google Scholar] [CrossRef]
  62. Allam, A.E.; Assaf, H.K.; Hassan, H.A.; Shimizu, K.; Elshaier, Y.A.M.M. An in silico perception for newly isolated flavonoids from peach fruit as privileged avenue for a countermeasure outbreak of COVID-19. RSC Adv. 2020, 10, 29983–29998. [Google Scholar] [CrossRef]
  63. Ibrahim, A.K.; Youssef, A.I.; Arafa, A.S.; Ahmed, S.A. Anti-H5N1 virus flavonoids fromCapparis sinaicaVeill. Nat. Prod. Res. 2013, 27, 2149–2153. [Google Scholar] [CrossRef] [PubMed]
  64. Van Vuong, Q.; Nguyen, T.T.; Li, M.S. A New Method for Navigating Optimal Direction for Pulling Ligand from Binding Pocket: Application to Ranking Binding Affinity by Steered Molecular Dynamics. J. Chem. Inf. Model. 2015, 55, 2731–2738. [Google Scholar] [CrossRef] [PubMed]
  65. Abdissa, D.; Geleta, G.; Bacha, K.; Abdissa, N. Phytochemical investigation of Aloe pulcherrima roots and evaluation for its antibacterial and antiplasmodial activities. PLoS ONE 2017, 12, e0173882. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Abdissa, N.; Induli, M.; Fitzpatrick, P.; Alao, J.P.; Sunnerhagen, P.; Landberg, G.; Yenesew, A.; Erdélyi, M. Cytotoxic Quinones from the Roots of Aloe dawei. Molecules 2014, 19, 3264–3273. [Google Scholar] [CrossRef] [Green Version]
  67. Abdissa, N.; Gohlke, S.; Frese, M.; Sewald, N. Cytotoxic Compounds from Aloe megalacantha. Molecules 2017, 22, 1136. [Google Scholar] [CrossRef]
  68. Zhong, J.; Huang, Y.; Ding, W.; Wu, X.; Wan, J.; Luo, H. Chemical constituents of Aloe barbadensis Miller and their inhibitory effects on phosphodiesterase-4D. Fitoterapia 2013, 91, 159–165. [Google Scholar] [CrossRef]
  69. Sun, Y.N.; Li, L.Y.; Li, W.; Kang, J.S.; Hwang, I.; Kim, Y.H. Chemical Components from Aloe and their Inhibition of Indoleamine 2, 3-dioxygenase. Pharmacogn. Mag. 2017, 13, 58–63. [Google Scholar]
  70. Kurizaki, A.; Watanabe, T.; Devkota, H.P. Chemical Constituents from the Flowers of Aloe arborescens. Nat. Prod. Commun. 2019, 14, 1934578–19844135. [Google Scholar] [CrossRef] [Green Version]
  71. Koyama, J.; Ogura, T.; Tagahara, K. Naphtho[2,3-c]furan-4,9-dione and its derivatives from Aloe ferox. Phytochemistry 1994, 37, 1147–1148. [Google Scholar] [CrossRef]
  72. Rehman, N.U.; Al-Riyami, S.A.; Hussain, H.; Ali, A.; Khan, A.L.; Al-Harrasi, A. Secondary metabolites from the resins of Aloe vera and Commiphora mukul mitigate lipid peroxidation. Acta Pharm. 2019, 69, 433–441. [Google Scholar] [CrossRef] [Green Version]
  73. Sigler, A.; Rauwald, H.W. Aloe Plants Accumulate Anthrone-Type Anthranoids in Inflorescence and Leaves, and Tetrahydroanthracenes in Roots. Z. Nat. C 1994, 49, 286–292. [Google Scholar] [CrossRef]
  74. Van Oudtshoorn, M.v.R. Chemotaxonomic investigations in asphodeleae and aloineae (liliaceae). Phytochemistry 1964, 3, 383–390. [Google Scholar] [CrossRef]
  75. Kambizi, L.; Sultana, N.; Afolayan, A. Bioactive Compounds Isolated fromAloe ferox.: A Plant Traditionally Used for the Treatment of Sexually Transmitted Infections in the Eastern Cape, South Africa. Pharm. Biol. 2005, 42, 636–639. [Google Scholar] [CrossRef]
  76. Awe, W.; Kuemmell, H. On the occurrence of aloin in Aloe vera in addition to comparative studies with the fresh juice of Cape aloe (Aloe ferox) and the dried extract prepared from it. Arch. Pharm. 1962, 295, 819–822. [Google Scholar] [CrossRef]
  77. Hirata, T.; Suga, T. Biologically Active Constituents of Leaves and Roots of Aloe arborescens var. natalensis. Z. Nat. C 1977, 32, 731–734. [Google Scholar] [CrossRef]
  78. Yagi, A.; Makino, K.; Nishioka, I. Studies on the Constituents of Aloe sapnaria HAW. I. The Structures of Tetrahydroanthracene Derivatives and the Related Anthraquinones. Chem. Pharm. Bull. 1974, 22, 1159–1166. [Google Scholar] [CrossRef] [Green Version]
  79. Yagi, A.; Makino, K.; Nishioka, I. Studies on the constituents of Aloe saponaria Haw. II. The structures of tetrahydroanthracene derivatives, aloesaponol III and -IV. Chem. Pharm. Bull. 1977, 25, 1764–1770. [Google Scholar] [CrossRef] [Green Version]
  80. Yagi, A.; Makino, K.; Nishioka, I. Studies on the constituents of Aloe saponaria Haw. III. The structures of phenol glucosides. Chem. Pharm. Bull. 1977, 25, 1771–1776. [Google Scholar] [CrossRef] [Green Version]
  81. Yagi, A.; Makino, K.; Nishioka, I. Studies on the constituents of Aloe saponaria Haw. IV. The structures of bianthraquinoid pigments. Chem. Pharm. Bull. 1978, 26, 1111–1116. [Google Scholar] [CrossRef] [Green Version]
  82. Abd-Alla, H.I.; Shaaban, M.; Shaaban, K.A.; Abu-Gabal, N.S.; Shalaby, N.M.; Laatsch, H. New bioactive compounds fromAloe hijazensis. Nat. Prod. Res. 2009, 23, 1035–1049. [Google Scholar] [CrossRef]
  83. Epifano, F.; Fiorito, S.; Locatelli, M.; Taddeo, V.A.; Genovese, S. Screening for novel plant sources of prenyloxyanthraquinones: Senna alexandrina Mill. and Aloe vera (L.) Burm. F. Nat. Prod. Res. 2015, 29, 180–184. [Google Scholar] [CrossRef]
  84. Van Wyk, B.-E.; Wink, M. Medicinal Plants of the World; CABI: Wallingford, England, 2018. [Google Scholar]
  85. Conner, J.M.; Gray, A.I.; Reynolds, T.; Waterman, P.G. Anthraquinone, anthrone and phenylpyrone components of Aloe nyeriensis var. kedongensis leaf exudate. Phytochemistry 1987, 26, 2995–2997. [Google Scholar] [CrossRef]
  86. Rehman, N.U.; Khan, A.; Al-Harrasi, A.; Khiat, M.; Hussain, H.; Wadood, A.; Riaz, M. Natural urease inhibitors from Aloe vera resin and Lycium shawii and their structural-activity relationship and molecular docking study. Bioorg. Chem. 2019, 88, 102955. [Google Scholar] [CrossRef]
  87. Conner, J.M.; Gray, A.I.; Waterman, P.G.; Reynolds, T. Novel Anthrone-Anthraquinone Dimers from Aloe elgonica. J. Nat. Prod. 1990, 53, 1362–1364. [Google Scholar] [CrossRef]
  88. Shin, K.H.; Woo, W.S.; Lim, S.S.; Shim, C.S.; Chung, H.S.; Kennelly, E.J.; Kinghorn, A.D. Elgonica-Dimers A and B, Two Potent Alcohol Metabolism Inhibitory Constituents of Aloe arborescens. J. Nat. Prod. 1997, 60, 1180–1182. [Google Scholar] [CrossRef]
  89. Saleem, R.; Faizi, S.; Hussain, S.A.; Qazi, A.; Dar, A.; Ahmad, S.I.; Qazi, M.; Akhtar, S.; Hasnain, S.N.; Siddiqui, B.S.; et al. Hypotensive Effect of Chemical Constituents fromAloe barbadensis. Planta Med. 2001, 67, 757–760. [Google Scholar] [CrossRef] [PubMed]
  90. Auterhoff, H.; Graf, E.; Eurisch, G.; Alexa, M. Trennung des Aloins in Diastereomere und deren Charakterisierung. Arch. Pharm. 1980, 313, 113–120. [Google Scholar] [CrossRef]
  91. Tekassa, T.; Tewabe, Y.; Bisrat, D.; Hailu, A.; Asres, K. Antileishmanial activities of leaf latex and compound isolated from Aloe ghibensis Sebsebe & Friis. Ethiop. Pharm. J. 2020, 35, 51–58. [Google Scholar] [CrossRef]
  92. Aldayel, T.S.; Grace, M.H.; Lila, M.A.; Yahya, M.A.; Omar, U.M.; Alshammary, G. LC-MS characterization of bioactive metabolites from two Yemeni Aloe spp. with antioxidant and antidiabetic properties. Arab. J. Chem. 2020, 13, 5040–5049. [Google Scholar] [CrossRef]
  93. Megeressa, M.; Bisrat, D.; Mazumder, A.; Asres, K. Structural elucidation of some antimicrobial constituents from the leaf latex of Aloe trigonantha L.C. Leach. BMC Complement. Altern. Med. 2015, 15, 270. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  94. Teka, T.; Kassahun, H. Characterization and Evaluation of Antioxidant Activity of Aloe schelpei Reynolds. Drug Des. Dev. Ther. 2020, ume 14, 1003–1008. [Google Scholar] [CrossRef] [Green Version]
  95. Kassahun, A.; Bisrat, D. Free radical scavenging activities of three anthrones from Aloe gilbertii reynolds leaf latex. J. Nat. Prod. Plant Resour. 2017, 7, 40–44. [Google Scholar]
  96. Abeje, F.; Bisrat, D.; Hailu, A.; Asres, K. Phytochemistry and Antileishmanial Activity of the Leaf Latex ofAloe calidophilaReynolds. Phytother. Res. 2014, 28, 1801–1805. [Google Scholar] [CrossRef]
  97. Okamura, N.; Hine, N.; Harada, S.; Fujioka, T.; Mihashi, K.; Nishi, M.; Miyahara, K.; Yagi, A. Diastereomeric C-glucosylanthrones of Aloe vera leaves. Phytochemistry 1997, 45, 1519–1522. [Google Scholar] [CrossRef]
  98. Grace, O.; Kokubun, T.; Veitch, N.; Simmonds, M. Characterisation of a nataloin derivative from Aloe ellenbeckii, a maculate species from east Africa. South. Afr. J. Bot. 2008, 74, 761–763. [Google Scholar] [CrossRef]
  99. Rauwald, H.W.; Beil, A. 5-Hydroxyaloin A in the Genus Aloe Thin Layer Chromatographic Screening and High Performance Liquid Chromatographic Determination. Z. Nat. C 1993, 48, 1–4. [Google Scholar] [CrossRef]
  100. Asres, K.; Girma, B.; Bisrat, D. Antimalarial evaluation of the leaf latex of Aloe citrina and its major constituent. Anc. Sci. Life 2015, 34, 142–146. [Google Scholar] [CrossRef]
  101. Mebe, P.P. 2′-p-Methoxycoumaroylaloeresin, A C-glucoside from Aloe excelsa. Phytochemistry 1987, 26, 2646–2647. [Google Scholar] [CrossRef]
  102. Beaumont, J.; Reynolds, R.; Vaughan, J.G. Homonataloin in Aloe species. Planta Med. 1984, 50, 505–508. [Google Scholar] [CrossRef]
  103. Dagne, E.; Bisrat, D.; Van Wyk, B.-E.; Viljoen, A.; Hellwig, V.; Steglich, W. Anthrones from Aloe microstigma. Phytochemistry 1997, 44, 1271–1274. [Google Scholar] [CrossRef]
  104. Lv, L.; Tian, X.-Y.; Fang, W.-S. Three new antioxidantC-glucosylanthrones fromAloe nobilis. J. Asian Nat. Prod. Res. 2010, 12, 443–447. [Google Scholar] [CrossRef]
  105. Bisrat, D.; Dagne, E.; Van Wyk, B.-E.; Viljoen, A. Chromones and anthrones from Aloe marlothii and Aloe rupestris. Phytochemistry 2000, 55, 949–952. [Google Scholar] [CrossRef]
  106. Sigler, A.; Rauwald, H.W. First proof of anthrone aglycones and diastereomeric anthrone-C-glycosyls in flowers and bracts of Aloe species. Biochem. Syst. Ecol. 1994, 22, 287–290. [Google Scholar] [CrossRef]
  107. Hoerhammer, L.; Wagner, H.; Bittner, G. ALOINOSIDE B, A NEW GLYCOSIDE FROM ALOE. Z. Nat. B 1964, 19, 222–226. [Google Scholar]
  108. Gao, J.; Zhang, G.; Dai, R.; Bi, K. Isolation of Aloinoside B and Metabolism by Rat Intestinal Bacteria. Pharm. Biol. 2005, 42, 581–587. [Google Scholar] [CrossRef]
  109. Viljoen, A.M.; Van Wyk, B.-E.; Van Heerden, F.R. The chemotaxonomic value of the diglucoside anthrone homonataloside B in the genus Aloe. Biochem. Syst. Ecol. 2002, 30, 35–43. [Google Scholar] [CrossRef]
  110. Farah, M.; Andersson, R.; Samuelsson, G. Microdontin A and B: Two New Aloin Derivatives from Aloe microdonta. Planta Med. 1992, 58, 88–93. [Google Scholar] [CrossRef]
  111. Zhong, J.-S.; Huang, Y.-Y.; Zhang, T.-H.; Liu, Y.-P.; Ding, W.-J.; Wu, X.-F.; Xie, Z.-Y.; Luo, H.-B.; Wan, J.-Z. Natural phosphodiesterase-4 inhibitors from the leaf skin of Aloe barbadensis Miller. Fitoterapia 2015, 100, 68–74. [Google Scholar] [CrossRef]
  112. Dagne, E.; Bisrat, D.; Van Wyk, B.-E.; Viljoen, A. 10-Hydroxyaloin B 6’-O-Acetate, an Oxanthrone fromAloe claviflora. J. Nat. Prod. 1998, 61, 256–257. [Google Scholar] [CrossRef]
  113. Dagne, E.; Van Wyk, B.-E.; Stephenson, D.; Steglich, W. Three oxanthrones from Aloe littoralis. Phytochemistry 1996, 42, 1683–1687. [Google Scholar] [CrossRef]
  114. Dagne, E.; Bisrat, D.; Codina, C.; Bastida, J. A C,O-diglucosylated oxanthrone from Aloe littoralis. Phytochemistry 1998, 48, 903–905. [Google Scholar] [CrossRef]
  115. Wang, H.M.; Shi, W.; Xu, Y.K.; Wang, P.; Chen, W.; Liu, Y.; Lu, M.J.; Pa, J.Q. Isolation and spectral study of 4-methyl-6, 8-dihydroxy-7H-benz [de] anthracen-7-one. Magn. Reson. Chem. 2003, 41, 301–303. [Google Scholar] [CrossRef]
  116. Asamenew, G.; Bisrat, D.; Mazumder, A.; Asres, K. In Vitro Antimicrobial and Antioxidant Activities of Anthrone and Chromone from the Latex of Aloe harlana Reynolds. Phytother. Res. 2011, 25, 1756–1760. [Google Scholar] [CrossRef] [PubMed]
  117. Saleem, R.; Faizi, S.; Deeba, F.; Siddiqui, B.S.; Qazi, M.H. Anthrones from Aloe barbadensis. Phytochemistry 1997, 45, 1279–1282. [Google Scholar] [CrossRef]
  118. Yagi, A.; Shoyama, Y.; Nishioka, I. Formation of tetrahydroanthracene glucosides by callus tissue of Aloe saponaria. Phytochemistry 1983, 22, 1483–1484. [Google Scholar] [CrossRef]
  119. Yagi, A.; Hine, N.; Asai, M.; Nakazawa, M.; Tateyama, Y.; Okamura, N.; Fujioka, T.; Mihashi, K.; Shimomura, K. Tetrahydroanthracene glucosides in callus tissue from Aloe barbadensis leaves. Phytochemistry 1998, 47, 1267–1270. [Google Scholar] [CrossRef]
  120. Yenesew, A.; Ogur, J.; Duddeckt, H. (R)-Prechrysophanol from Aloe graminicola. Phytochemistry 1993, 34, 1442–1444. [Google Scholar] [CrossRef]
  121. Dagne, E.; Casser, I.; Steglich, W. Aloechrysone, a dihydroanthracenone from Aloe berhana. Phytochemistry 1992, 31, 1791–1793. [Google Scholar] [CrossRef]
  122. Abe, I.; Oguro, S.; Utsumi, Y.; Sano, Y.; Noguchi, H. Engineered Biosynthesis of Plant Polyketides: Chain Length Control in an Octaketide-Producing Plant Type III Polyketide Synthase. J. Am. Chem. Soc. 2005, 127, 12709–12716. [Google Scholar] [CrossRef]
  123. Sun, Y.N.; Li, W.; Yang, S.Y.; Kang, J.S.; Ma, J.Y.; Kim, Y.H. Isolation and identification of chromone and pyrone constituents from Aloe and their anti-inflammatory activities. J. Funct. Foods 2016, 21, 232–239. [Google Scholar] [CrossRef]
  124. Holdsworth, D. Chromones in Aloe species–Part II–Aloesone. Planta Med. 1972, 22, 54–58. [Google Scholar] [CrossRef]
  125. Conner, J.M.; Gray, A.I.; Reynolds, T.; Waterman, P.G. Anthrone and chromone components of Aloe cremnophila and A. jacksonii leaf exudates. Phytochemistry 1990, 29, 941–944. [Google Scholar] [CrossRef]
  126. Speranza, G.; Fontana, G.; Zanzola, S.; Di Meo, A. Studies on Aloe. 15.1Two New 5-Methylchromones from Cape Aloe. J. Nat. Prod. 1997, 60, 692–694. [Google Scholar] [CrossRef]
  127. Durı, L.; Morelli, C.; Crippa, S.; Speranza, G. 6-Phenylpyrones and 5-methylchromones from Kenya aloe. Fitoterapia 2004, 75, 520–522. [Google Scholar] [CrossRef] [PubMed]
  128. Blitzke, T.; Masaoud, M.; Schmidt, J. Constituents of Aloe rubroviolacea. Fitoterapia 2001, 72, 78–79. [Google Scholar] [CrossRef]
  129. Okamura, N.; Hine, N.; Harada, S.; Fujioka, T.; Mihashi, K.; Yagi, A. Three chromone components from Aloe vera leaves. Phytochemistry 1996, 43, 495–498. [Google Scholar] [CrossRef]
  130. Okamura, N.; Hine, N.; Tateyama, Y.; Nakazawa, M.; Fujioka, T.; Mirmhi, K.; Yagi, A. Three chromones of Aloe vera leaves. Phytochemistry 1997, 45, 1511–1513. [Google Scholar] [CrossRef]
  131. Okamura, N.; Hine, N.; Tateyama, Y.; Nakazawa, M.; Fujioka, T.; Mihashi, K.; Yagi, A. Five chromones from Aloe Vera leaves. Phytochemistry 1998, 49, 219–223. [Google Scholar] [CrossRef]
  132. Kahramanoğlu, I.; Chen, C.; Chen, J.; Wan, C. Chemical Constituents, Antimicrobial Activity, and Food Preservative Characteristics of Aloe vera Gel. Agronomy 2019, 9, 831. [Google Scholar] [CrossRef] [Green Version]
  133. Yuan, A.X.; Kang, S.H.; Qin, L.; Yuan, P.; Fan, Y.J. Chemical constituents of the leaves of Chinese aloe (Aloe vera var. chinensis). Zhongcaoyao 1994, 25, 339–341. [Google Scholar]
  134. Park, M.; Park, J.; Shin, Y.; Kim, W.; Lee, J.; Kim, K. Neoaloesin A: A NewC-Glucofuranosyl Chromone fromAloe barbadensis. Planta Med. 1996, 62, 363–365. [Google Scholar] [CrossRef] [PubMed]
  135. Haynes, L.; Holdsworth, D.; Russell, R. C-Glycosyl compounds. Part VI. Aloesin, a C-glucosylchromone from Aloe sp. J. Chem. Soc. C Org. 1970, 2581–2586. [Google Scholar] [CrossRef]
  136. Wu, X.-F.; Wan, J.-Z.; Luo, B.-J.; Yang, M.-R.; Ding, W.-J.; Zhong, J.-S. A novel naphthalene derivative from Aloe barbadensis. Yao xue xue bao = Acta Pharm. Sin. 2013, 48, 723–727. [Google Scholar]
  137. Hiruy, M.; Bisrat, D.; Mazumder, A.; Asres, K. Two chromones with antimicrobial activity from the leaf latex of Aloe monticola Reynolds. Nat. Prod. Res. 2021, 35, 1052–1056. [Google Scholar] [CrossRef]
  138. Lv, L.; Yang, Q.Y.; Zhao, Y.; Yao, C.S.; Sun, Y.; Yang, E.J.; Song, K.S.; Mook-Jung, I.; Fang, W.S. BACE1 (beta-secretase) inhibitory chromone glycosides from Aloe vera and Aloe nobilis. Planta Med. 2008, 74, 540–545. [Google Scholar] [CrossRef]
  139. Hutter, J.A.; Salman, M.; Stavinoha, W.B.; Satsangi, N.; Williams, R.F.; Streeper, R.T.; Weintraub, S.T. AntiinflammatoryC-Glucosyl Chromone fromAloe barbadensis. J. Nat. Prod. 1996, 59, 541–543. [Google Scholar] [CrossRef]
  140. Speranza, G.; Gramatica, P.; Dadá, G.; Manitto, P. Aloeresin C, a bitter C, O-diglucoside from Cape Aloe. Phytochemistry 1985, 24, 1571–1573. [Google Scholar] [CrossRef]
  141. Speranza, G.; Dada, G.; Lunazzi, L.; Gramatica, P.; Manitto, P. A C-glucosylated 5-methylchromone from Kenya aloe. Phytochemistry 1986, 25, 2219–2222. [Google Scholar] [CrossRef]
  142. Conner, J.M.; Gray, A.I.; Reynolds, T.; Waterman, P.G. Anthracene and chromone derivatives in the exudate of Aloe rabaiensis. Phytochemistry 1989, 28, 3551–3553. [Google Scholar] [CrossRef]
  143. Van Heerden, F.R.; Van Wyk, B.-E.; Viljoen, A.M. ChemInform Abstract: Aloeresins E (Ia) and F (Ib), Two Chromone Derivatives from Aloe Peglerae. ChemInform 2010, 28, 867–869. [Google Scholar] [CrossRef]
  144. Van Heerden, F.R.; Viljoen, A.M.; van Wyk, B.E. 6’-O-Coumaroylaloesin from Aloe castanea--a taxonomic marker for Aloe section Anguialoe. Phytochemistry 2000, 55, 117–120. [Google Scholar] [CrossRef]
  145. Holzapfel, C.W.; Wessels, P.L.; Van Wyk, B.-E.; Marais, W.; Portwig, M. Chromone and aloin derivatives from Aloe broomii, A. Africana and A. speciosa. Phytochemistry 1997, 45, 97–102. [Google Scholar] [CrossRef]
  146. Gramatica, P.; Monti, D.; Speranza, G.; Manitto, P. Aloe revisited the structure of aloeresin A. Tetrahedron Lett. 1982, 23, 2423–2424. [Google Scholar] [CrossRef]
  147. Speranza, G.; Martignoni, A.; Manitto, P. Iso-aloeresin A, a Minor Constituent of Cape Aloe. J. Nat. Prod. 1988, 51, 588–590. [Google Scholar] [CrossRef]
  148. Makino, K.; Yagi, A.; Nishioka, I. Studies on the constituents of Aloe arborescens Mill. var. natalensis Berger. II. The structures of two new aloesin esters. Chem. Pharm. Bull. 1974, 22, 1565–1570. [Google Scholar] [CrossRef] [Green Version]
  149. Xiao, Z.; Chen, D.; Si, J.; Tu, G.; Ma, L. Chemical constituents of Aloe vera. Yaoxue Xuebao 2000, 35, 120–123. [Google Scholar]
  150. Manitto, P.; Speranza, G.; De Tommasi, N.; Ortoleva, E.; Morelli, C.F. Aloeresin H, a new polyketide constituent of Cape aloe. Tetrahedron 2003, 59, 401–408. [Google Scholar] [CrossRef]
  151. Speranza, G.; Morelli, C.F.; Tubaro, A.; Altinier, G.; Durì, L.; Manitto, P. Aloeresin I, an Anti-Inflammatory 5-Methylchromone from Cape Aloe. Planta Med. 2005, 71, 79–81. [Google Scholar] [CrossRef]
  152. Yuan, A.X. The molecular structure of iso-aloesin isolated from the leaves of Aloe vera L. var. chinensis (Haw.) Berge. China J. Chin. Mater. Med. 1993, 18, 609–611. [Google Scholar]
  153. Bhaludra, C.S.S.; Bethapudi, R.R.; Murugulla, A.C.; Pullagummi, C.; Latha, T.; Venkatesh, K.; Bheemagani, A.J.; Pudutha, A.; Rani, A.R. Cultivation, phytochemical studies, biological activities and medicinal uses of Aloe ferox, grandfather of aloes an important amazing medicinal plant. Int. J. Pharmacol. 2013, 9, 405–415. [Google Scholar]
  154. López, A.; De Tangil, M.S.; Vega-Orellana, O.; Ramírez, A.S.; Rico, M. Phenolic Constituents, Antioxidant and Preliminary Antimycoplasmic Activities of Leaf Skin and Flowers of Aloe vera (L.) Burm. f. (syn. A. barbadensis Mill.) from the Canary Islands (Spain). Molecules 2013, 18, 4942–4954. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  155. Saxena, V.; Sharma, D. 5, 4’-dihydroxy 6, 7, 3’, 5’-tetramethoxy flavone 5-O-alpha-L-rhamno pyranosyl 16-O-beta-D-galactopyranoside from aloe barbadensis (leaves). J. Inst. Chem. 1998, 70, 179–182. [Google Scholar]
  156. Sexena, V.; Chourasia, S. 7-Hydroxy, 6, 3′,4′-Trimethoxy Isoflavone-5-O-alpha-L-Rhamnopy-ranosyl[1→6]-O-beta-D-Glucopyranoside of Aloe vera. J. Inst. Chem. 2000, 72, 195–197. [Google Scholar]
  157. Kametani, S.; Kojima-Yuasa, A.; Kikuzaki, H.; Kennedy, D.O.; Honzawa, M.; Matsui-Yuasa, I. Chemical Constituents of Cape Aloe and Their Synergistic Growth-Inhibiting Effect on Ehrlich Ascites Tumor Cells. Biosci. Biotechnol. Biochem. 2007, 71, 1220–1229. [Google Scholar] [CrossRef]
  158. Makino, K.; Yagi, A.; Nishioka, I. Studies on the Constituents of Aloe arborescens MILL. var. natalensis BERGER. I. The Structure of Aloearbonaside, a Glucoside of a New Type naturally Occurring Chromene. Chem. Pharm. Bull. 1973, 21, 149–156. [Google Scholar] [CrossRef] [Green Version]
  159. Speranza, G.; Dadá, G.; Lunazzi, L.; Gramatica, P.; Manitto, P. Aloenin B, a New Diglucosylated 6-Phenyl-2-pyrone from Kenya Aloe. J. Nat. Prod. 1986, 49, 800–805. [Google Scholar] [CrossRef]
  160. Shi, S.-P.; Wanibuchi, K.; Morita, H.; Endo, K.; Noguchi, H.; Abe, I. Enzymatic Formation of Unnatural Novel Chalcone, Stilbene, and Benzophenone Scaffolds by Plant Type III Polyketide Synthase. Org. Lett. 2009, 11, 551–554. [Google Scholar] [CrossRef] [PubMed]
  161. Zhang, X.-F.; Wang, H.-M.; Song, Y.-L.; Nie, L.-H.; Wang, L.-F.; Liu, B.; Shen, P.-P.; Liu, Y. Isolation, structure elucidation, antioxidative and immunomodulatory properties of two novel dihydrocoumarins from Aloe vera. Bioorg. Med. Chem. Lett. 2006, 16, 949–953. [Google Scholar] [CrossRef] [PubMed]
  162. Speranza, G.; Manitto, P.; Cassara’, P.; Monti, D. Feralolide, a dihydroisocoumarin from cape aloe. Phytochemistry 1993, 33, 175–178. [Google Scholar] [CrossRef]
  163. Rehman, N.U.; Hussain, H.; Khiat, M.; Khan, H.Y.; Abbas, G.; Green, I.R.; Al-Harrasi, A. Bioactive chemical constituents from the resin of Aloe vera. Z. Nat. B 2017, 72, 955–958. [Google Scholar] [CrossRef]
  164. Wang, H.M.; Shi, W.; Xu, Y.K.; Liu, Y.; Lü, M.J.; Pan, J.Q. Spectral study of a new dihydroisocoumarin. Magn. Reson. Chem. 2003, 41, 718–720. [Google Scholar] [CrossRef]
  165. Veitch, N.C.; Simmonds, M.S.; Blaney, W.M.; Reynolds, T. A dihydroisocoumarin glucoside from Aloe hildebrandtii. Phytochemistry 1994, 35, 1163–1166. [Google Scholar] [CrossRef]
  166. Saleem, R.; Faizi, S.; Deeba, F.; Siddiqui, B.S.; Qazi, M.H. A New Bisbenzopyran fromAloe barbadensisRoots. Planta Med. 1997, 63, 454–456. [Google Scholar] [CrossRef] [PubMed]
  167. Dagne, E.; Yenesew, A.; Asmellash, S.; Demissew, S.; Mavi, S. Anthraquinones, pre-anthraquinones and isoeleutherol in the roots of Aloe species. Phytochemistry 1994, 35, 401–406. [Google Scholar] [CrossRef]
  168. Kong, W.-S.; Li, J.; Liu, X.; Mi, Q.-L.; Chen, J.-H.; Li, X.-M.; Yang, G.-Y.; Hu, Q.-F.; Li, T.; Yang, Y.-K. A new naphthalene derivative from Aloe vera and its antibacterial activity. China J. Chin. Mater. Med. 2017, 42, 3761–3763. [Google Scholar]
  169. Wessels, P.L.; Holzapfel, C.W.; Van Wyk, B.-E.; Marais, W. Plicataloside, an O,O-diglycosylated naphthalene derivative from Aloe plicatilis. Phytochemistry 1996, 41, 1547–1551. [Google Scholar] [CrossRef]
  170. Speranza, G.; Monti, D.; Crippa, S.; Cairoli, P.; Morelli, C.F.; Manitto, P. Kenyaloside, a Novel O,O,O-Triglycosylated Naphthalene Derivative from the Exudate of Kenyan Aloe Species. Nat. Prod. Commun. 2006, 1, 1085–1088. [Google Scholar] [CrossRef]
  171. Yang, Q.-Y.; Yao, C.-S.; Fang, W.-S. A new triglucosylated naphthalene glycoside from Aloe vera L. Fitoterapia 2010, 81, 59–62. [Google Scholar] [CrossRef]
  172. Speranza, G.; Di Meo, A.; Manitto, P.; Monti, D.; Fontana, G. A New Benzochromanone Derivative from Cape Aloe†. J. Agric. Food Chem. 1996, 44, 274–277. [Google Scholar] [CrossRef]
  173. Blitzke, T.; Porzel, A.; Masaoud, M.; Schmidt, J. A chlorinated amide and piperidine alkaloids from Aloe sabaea. Phytochemistry 2000, 55, 979–982. [Google Scholar] [CrossRef]
  174. Hotti, H.; Häkkinen, S.T.; Seppänen-Laakso, T.; Rischer, H. Polyketide-Derived Alkaloids and Anthraquinones in Aloe Plants and Cell Cultures. J. Plant. Biotechnol. Res. 2019, 1, 1–15. [Google Scholar] [CrossRef] [PubMed]
  175. Dagne, E.; Bisrat, D.; Viljoen, A.; Van Wyk, B.-E. Chemistry of Aloe Species. Curr. Org. Chem. 2000, 4, 1055–1078. [Google Scholar] [CrossRef]
  176. Prasannaraja, C.; Kamalanathan, A.S.; Vijayalakshmi, M.A.; Venkataraman, K. A dipyrrole derivative from Aloe vera inhibits an anti-diabetic drug target Dipeptidyl Peptidase (DPP)-IV in vitro. Prep. Biochem. Biotechnol. 2020, 50, 511–520. [Google Scholar] [CrossRef]
  177. Speranza, G.; Corti, S.; Manitto, P. Isolation and Chemical Characterization of a New Constituent of Cape Aloe Having the 1,1-Diphenylethane Skeleton. J. Agric. Food Chem. 1994, 42, 2002–2006. [Google Scholar] [CrossRef]
  178. Confalone, P.N.; Huie, E.M.; Patel, N.G. ChemInform Abstract: The isolation, structure determination, and synthesis of pluridone, a novel insecticide from aloe pluridens. Chem. Inf. 1984, 15, 5563–5566. [Google Scholar] [CrossRef]
  179. Speranza, G.; Manitto, P.; Monti, D.; Lianza, F. Feroxidin, a novel 1-methyltetralin derivative isolated from cape aloe. Tetrahedron Lett. 1990, 31, 3077–3080. [Google Scholar] [CrossRef]
  180. Speranza, G.; Paolo, M.; Donata, P.; Diego, M. Absolute configuration of feroxidin: An experimental support to Snatzke’s helicity rules for tetralins. Chirality 1991, 3, 263–267. [Google Scholar] [CrossRef]
  181. Speranza, G.; Manitto, P.; Monti, D.; Pezzuto, D. Studies on Aloe, Part 10. Feroxins A and B, Two O-Glucosylated 1-Methyltetralins from Cape Aloe. J. Nat. Prod. 1992, 55, 723–729. [Google Scholar] [CrossRef]
  182. Esua, M.F.; Rauwald, J.-W. Novel bioactive maloyl glucans from Aloe vera gel: Isolation, structure elucidation and in vitro bioassays. Carbohydr. Res. 2006, 341, 355–364. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The percentage of different classes of phytochemicals reported from the genus Aloe. Anthraquinones 36.3%, chromones 27.4%, coumarin 0.8%, flavonoids 4.2%, simple phenolic compounds 8%, phenyl pyran and phenyl pyrone derivatives 7%, benzofurans 2%, naphthalene derivatives 5.9%, alkaloids 1.2%, fatty acid derivatives 1.2% and miscellaneous compounds 5.5%.
Figure 1. The percentage of different classes of phytochemicals reported from the genus Aloe. Anthraquinones 36.3%, chromones 27.4%, coumarin 0.8%, flavonoids 4.2%, simple phenolic compounds 8%, phenyl pyran and phenyl pyrone derivatives 7%, benzofurans 2%, naphthalene derivatives 5.9%, alkaloids 1.2%, fatty acid derivatives 1.2% and miscellaneous compounds 5.5%.
Molecules 26 01767 g001
Figure 2. Three-dimensional crystal structure of the therapeutic targets of SARS-CoV-2 Mpro main protease (PDB ID: 6LU7, cyan) and spike glycoprotein (PDB ID: 6M0J, magenta).
Figure 2. Three-dimensional crystal structure of the therapeutic targets of SARS-CoV-2 Mpro main protease (PDB ID: 6LU7, cyan) and spike glycoprotein (PDB ID: 6M0J, magenta).
Molecules 26 01767 g002
Figure 3. Docking scores and RMSD values of isolated compounds from Aloe genus against SARS-CoV-2 main protease Mpro.
Figure 3. Docking scores and RMSD values of isolated compounds from Aloe genus against SARS-CoV-2 main protease Mpro.
Molecules 26 01767 g003
Figure 4. Docking scores and RMSD values of isolated compounds from Aloe genus against SARS-CoV-2 spike glycoprotein.
Figure 4. Docking scores and RMSD values of isolated compounds from Aloe genus against SARS-CoV-2 spike glycoprotein.
Molecules 26 01767 g004
Figure 5. Chemical structures of top scoring active compounds from Aloe genus. Top panel, top scoring active compounds for SARS-CoV-2 protein Mpro. Lower panel, top scoring active compounds for SARS-CoV-2 protein spike glycoprotein.
Figure 5. Chemical structures of top scoring active compounds from Aloe genus. Top panel, top scoring active compounds for SARS-CoV-2 protein Mpro. Lower panel, top scoring active compounds for SARS-CoV-2 protein spike glycoprotein.
Molecules 26 01767 g005
Figure 6. 2D and 3D docking interactions complexes of compounds; 132 (A1,B1), 134 (A2,B2) and 159 (A3,B3) with SARS-CoV-2 main protease protein Mpro. Solid red line: Hydrogen bonds. Dashed pink line: Hydrophobic interactions.
Figure 6. 2D and 3D docking interactions complexes of compounds; 132 (A1,B1), 134 (A2,B2) and 159 (A3,B3) with SARS-CoV-2 main protease protein Mpro. Solid red line: Hydrogen bonds. Dashed pink line: Hydrophobic interactions.
Molecules 26 01767 g006
Figure 7. 2D and 3D docking interactions complexes of compounds 115 (A1,B1), 131 (A2,B2), and 120 (A3,B3) with SARS-CoV-2 spike glycoprotein. Solid red line: Hydrogen bonds. Dashed pink line: Hydrophobic interactions.
Figure 7. 2D and 3D docking interactions complexes of compounds 115 (A1,B1), 131 (A2,B2), and 120 (A3,B3) with SARS-CoV-2 spike glycoprotein. Solid red line: Hydrogen bonds. Dashed pink line: Hydrophobic interactions.
Molecules 26 01767 g007
Figure 8. Analysis of the molecular dynamics (MD) simulations for top scoring compounds from Aloe genus. (A) RMSD analysis for Mpro and the ligands 132, 134 and 159. (B) RMSD analysis for spike glycoprotein and the ligands 115, 120 and 131. (C) RMSF analysis for Mpro and the ligands 132, 134 and 159. (D) RMSF analysis for spike glycoprotein and the ligands 115, 120 and 131.
Figure 8. Analysis of the molecular dynamics (MD) simulations for top scoring compounds from Aloe genus. (A) RMSD analysis for Mpro and the ligands 132, 134 and 159. (B) RMSD analysis for spike glycoprotein and the ligands 115, 120 and 131. (C) RMSF analysis for Mpro and the ligands 132, 134 and 159. (D) RMSF analysis for spike glycoprotein and the ligands 115, 120 and 131.
Molecules 26 01767 g008
Table 1. Molecular docking results and interactions of the three top scoring compounds of Aloe species with SARS-CoV-2 proteins.
Table 1. Molecular docking results and interactions of the three top scoring compounds of Aloe species with SARS-CoV-2 proteins.
ProteinNo.Docking
Score
(kcal/mol)
RSMD 1
Refine
CLogPReceptorInteractionDistanceE (kcal/mol)
main Protease (PDB ID: 6LU7)132−7.681.370.25ASN142AH-donor3.01−1.9
ASN142AH-donor2.81−2.5
HIS163AH-acceptor3.35−0.7
GLN189AH-acceptor3.03−1.4
GLU166Api-H4.42−0.7
GLN189Api-H3.62−0.6
134−7.951.720.69ASN142AH-donor3.44−0.6
ASN142AH-donor2.78−2
HIS163AH-acceptor3.29−1.2
159−7.721.40−1.36THR190AH-donor2.94−0.8
Spike Glycoprotein (PDB ID: 6M0J)115−8.051.491.21TRP566AH-acceptor2.93−2.4
LYS562AH-acceptor3.11−12.2
LYS562AIonic3.11−3.8
VAL209pi-H4.17−0.6
VAL209pi-H4.27−0.6
131−8.081.131.02GLN102AH-donor3.01−1.2
ASN210AH-acceptor3.31−0.8
ASP206Api-H4.26−0.7
120−7.871.170.46ALA396AH-donor2.56−0.5
ASP206AH-donor2.78−3.0
GLU208AH-donor2.77−1.1
1 RMSD; Root mean square deviation.
Table 2. Drug-likeness and molecular properties of top scoring compounds predicted at molinspiration server. miLogP: Molinspiration LogP (Octanol-water partition coefficient); TPSA: Total polar surface area (drug transport properties); natoms-Number of atoms; Mol.Wt (g/mol): Molecular weight; nON: number of hydrogen bond acceptors; nOHHN: number of hydrogen bond donors; nviolations: Number of Lipinski’s rule of five parameters violations; nrotb: Number of Rotatable Bonds (molecular flexibility).
Table 2. Drug-likeness and molecular properties of top scoring compounds predicted at molinspiration server. miLogP: Molinspiration LogP (Octanol-water partition coefficient); TPSA: Total polar surface area (drug transport properties); natoms-Number of atoms; Mol.Wt (g/mol): Molecular weight; nON: number of hydrogen bond acceptors; nOHHN: number of hydrogen bond donors; nviolations: Number of Lipinski’s rule of five parameters violations; nrotb: Number of Rotatable Bonds (molecular flexibility).
115120131132134159
miLogP2.81.842.151.661.96−1.06
TPSA176.12196.35172.97193.19182.2269.43
natoms404140414243
MW556.56572.56554.55570.55584.57610.52
nON111211121216
nOHNH5645410
nviolations232223
nrotb9999106
volume486.85494.87480.99489.01506.53496.07
%ABS48.2341.2549.3242.3446.1416.04
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Abouelela, M.E.; Assaf, H.K.; Abdelhamid, R.A.; Elkhyat, E.S.; Sayed, A.M.; Oszako, T.; Belbahri, L.; El Zowalaty, A.E.; Abdelkader, M.S.A. Identification of Potential SARS-CoV-2 Main Protease and Spike Protein Inhibitors from the Genus Aloe: An In Silico Study for Drug Development. Molecules 2021, 26, 1767. https://doi.org/10.3390/molecules26061767

AMA Style

Abouelela ME, Assaf HK, Abdelhamid RA, Elkhyat ES, Sayed AM, Oszako T, Belbahri L, El Zowalaty AE, Abdelkader MSA. Identification of Potential SARS-CoV-2 Main Protease and Spike Protein Inhibitors from the Genus Aloe: An In Silico Study for Drug Development. Molecules. 2021; 26(6):1767. https://doi.org/10.3390/molecules26061767

Chicago/Turabian Style

Abouelela, Mohamed E., Hamdy K. Assaf, Reda A. Abdelhamid, Ehab S. Elkhyat, Ahmed M. Sayed, Tomasz Oszako, Lassaad Belbahri, Ahmed E. El Zowalaty, and Mohamed Salaheldin A. Abdelkader. 2021. "Identification of Potential SARS-CoV-2 Main Protease and Spike Protein Inhibitors from the Genus Aloe: An In Silico Study for Drug Development" Molecules 26, no. 6: 1767. https://doi.org/10.3390/molecules26061767

APA Style

Abouelela, M. E., Assaf, H. K., Abdelhamid, R. A., Elkhyat, E. S., Sayed, A. M., Oszako, T., Belbahri, L., El Zowalaty, A. E., & Abdelkader, M. S. A. (2021). Identification of Potential SARS-CoV-2 Main Protease and Spike Protein Inhibitors from the Genus Aloe: An In Silico Study for Drug Development. Molecules, 26(6), 1767. https://doi.org/10.3390/molecules26061767

Article Metrics

Back to TopTop