Journal of Environmental Science and Health, Part A Toxic/Hazardous Substances and Environmental Engineering

ISSN: (Print) (Online) Journal homepage: https://www.tandfonline.com/loi/lesa20

The inhibitory effect of some natural bioactive compounds against SARS-CoV-2 main protease: insights from molecular docking analysis and molecular dynamic simulation Doaa A. Abdelrheem , Shimaa A. Ahmed , H. R. Abd El-Mageed , Hussein S. Mohamed , Aziz A. Rahman , Khaled N. M. Elsayed & Sayed A. Ahmed To cite this article: Doaa A. Abdelrheem , Shimaa A. Ahmed , H. R. Abd El-Mageed , Hussein S. Mohamed , Aziz A. Rahman , Khaled N. M. Elsayed & Sayed A. Ahmed (2020): The inhibitory effect of some natural bioactive compounds against SARS-CoV-2 main protease: insights from molecular docking analysis and molecular dynamic simulation , Journal of Environmental Science and Health, Part A, DOI: 10.1080/10934529.2020.1826192 To link to this article: https://doi.org/10.1080/10934529.2020.1826192

Published online: 01 Oct 2020.

Submit your article to this journal

Article views: 84

View related articles

View Crossmark data

Full Terms & Conditions of access and use can be found at https://www.tandfonline.com/action/journalInformation?journalCode=lesa20

JOURNAL OF ENVIRONMENTAL SCIENCE AND HEALTH, PART A https://doi.org/10.1080/10934529.2020.1826192

The inhibitory effect of some natural bioactive compounds against SARS-CoV-2 main protease: insights from molecular docking analysis and molecular dynamic simulation Doaa A. Abdelrheema, Shimaa A. Ahmeda, H. R. Abd El-Mageedb, Hussein S. Mohamedc Khaled N. M. Elsayede, and Sayed A. Ahmeda

, Aziz A. Rahmand,

a Department of Chemistry, Faculty of Science, Beni-Suef University, Beni-Suef, Egypt; bMicro-Analysis and Environmental Research and Community Services Center, Faculty of Science, Beni-Suef University, Beni-Suef, Egypt; cResearch Institute of Medicinal and Aromatic Plants (RIMAP), Beni-Suef University, Beni-Suef, Egypt; dDepartment of Pharmacy, University of Rajshahi, Rajshahi, Bangladesh; eDepartment of Botany, Faculty of Science, Beni-Suef University, Beni-Suef, Egypt

ABSTRACT

ARTICLE HISTORY

This work aimed at evaluating the inhibitory effect of ten natural bioactive compounds (1–10) as potential inhibitors of SARS-CoV-2-3CL main protease (PDB ID: 6LU7) and SARS-CoV main proteases (PDB IDs: 2GTB and 3TNT) by molecular docking analysis. The inhibitory effect of all studied compounds was studied with compared to some proposed antiviral drugs which currently used in COVID-19 treatment such as chloroquine, hydroxychloroquine, azithromycin, remdesivir, baloxvir, lopinavir, and favipiravir. Homology modeling and sequence alignment was computed to evaluate the similarity between the SARS-CoV-2-3CL main protease and other SARS-CoV receptors. ADMET properties of all studied compounds were computed and reported. Also, molecular dynamic (MD) simulation was performed on the compound which has the highest binding affinity inside 6LU7 obtained from molecular docking analysis to study it is stability inside receptor in explicit water solvent. Based on molecular docking analysis, we found that caulerpin has the highest binding affinity inside all studied receptors compared to other bioactive compounds and studied drugs. Our homology modeling and sequence alignment showed that SARS-CoV main protease (PDB ID: 3TNT) shares high similarity with 3CLpro (96.00%). Also, ADMET properties confirmed that caulerpin obeys Lipinski’s rule and passes ADMET property, which make it a promising compound to act as a new safe natural drug against SARS-CoV-2-3CL main protease. Finally, MD simulation confirmed that the complex formed between caulerpin and 3CLpro is stable in water explicit and had no major effect on the flexibility of the protein throughout the simulations and provided a suitable basis for our study. Also, binding free energy between caulerpin and 6LU7 confirmed the efficacy of the caulerpin molecule against SARS-CoV-2 main protease. So, this study suggested that caulerpin could be used as a potential candidate in COVID-19 treatment.

Received 4 May 2020 Accepted 21 August 2020

Introduction Coronaviruses (CoV) are a great family of viruses varying from the common cold virus which is the lowest serious disease in this family to the more serious diseases such as MERS, CWID, and SARS. Also, their structure has a characteristic RNA genome. However, CoV are more familiar with animals; seven of them can induce the human respiratory system.[1] From December 2019, a new coronavirus called COVID19 that was first appeared in Wuhan, Hubei Province of China, and now, is spreading rapidly across China and other parts of the world. COVID-19 has become a serious threat to the world public health, causing 244,792 deaths from 3,484,640 cases as of 3 May 2020 from the entire world. The current COVID-19 outbreak is caused by SARS-CoV-2 that has been known as the seventh member of the family of coronaviruses.[1] SARS-CoV-2 was established to have higher sequence homology toward SARS-CoV than that of CONTACT Sayed A. Ahmed ß 2020 Taylor & Francis Group, LLC

[email protected]

KEYWORDS

COVID-19 virus protease; natural products; caulerpin; molecular docking; MD simulation

MERS-CoV according to the whole genome sequence alignment analysis in different studies.[2] There are four structural proteins characterize coronaviruse genome, spike glycoprotein (S), matrix glycoprotein (M), nucleocapsid protein (N), and small envelope protein (E) were reported.[3] Additionally, 3CLpro is also known as Nsp5 of SARS-Cov-2 required for the maturation of coronaviruses. 3CLpro is first automatically cleaved from poly-proteins to produce mature enzymes, and then, further cleaves downstream Nsps at 11 sites to release Nsp4–Nsp1631. 3CLpro directly mediates the maturation of Nsps. So, 3CLpro is critical for the viral life cycle, making it an attractive target of anti-coronavirus drug evolution.[4,5] There are different studies that had been performed on the main proteases of SARS-Coronavirus (2GTB) to develop antiviral treatments of COVID-19 virus because it shares 96% similarity with 3CLpro of COVID-19.[6]

Department of Chemistry, Faculty of Science, Beni-Suef University 62511, Beni-Suef, Egypt.

2

D. A. ABDELRHEEM ET AL.

Table 1. Source of studied compounds belong to natural sources and their biological activities. Compounds Oleic acid (1)

Saringosterols (2,3)

b-sitosterol (4)

Sources Terrestrial plant Red macroalgae Green macroalgae Brown macroalgae Microalgae Endophytic fungi form plant Vegetable oils Terrestrial plant Red macroalgae Green macroalgae Brown macroalgae

Species name Zaleya decandra[8] Ceramium virgatum[18] Ulva intestinalis[18] Fucus sp.[19] Chlorella vulgaris[20] Torreya Grandis[21] Hazelnut[22] Strychnos spinosa[25] Acanthophora spicifera[26] Cladophora fascicularis[27] Sargassum muticum[28]

Terrestrial plant Red macroalgae Green macroalgae Brown macroalgae Microalgae Terrestrial plant Red macroalgae Brown macroalgae Cyanobacteria Brown macroalgae

Synadenium glaucescens[31] Eucheuma cottonii[32] Ulva fasciata[33] Sargassum glaucescens[34] Nannochloropsis[35] Soybean[39] Exophyllum wentii[40] Sargassum horneri[41] Phormidium sp.[42] Sargassum naozhouense[44]

Loliolide (7)

Terrestrial plant Brown macroalgae

Canscora decussata[46] Sargassum naozhouense[44]

Hexadecanoic acid (8)

Terrestrial plant

Canthium parviflorum[49]

Caulerpin (10)

Green macroalgae Red macroalgae Brown macroalgae

Caulerpa racemosa[11] Chondria armata[12] Sargassum platycarpum[13]

Glycoglycerolipids (5,9) Kjellmanianone(6)

The utilization of natural compounds alternatives or complementary therapies has received growing attention due to their low toxicity, fewer side effects and in some cases as the only effective treatment. Chemical drugs are shown different side effects and ineffectiveness in some cases for longterm use. The World Health Organization (WHO) evaluated that about 80% of people use natural compounds in the treatment field.[7] Actually, natural compounds have various and effective biological activities such as antimicrobial, anticancer, anti-inflammatory, and anti-diabetic.[8] Medicinal plants are an important source of molecules with various pharmacological properties and effective biological activities like antimicrobial, anticancer, anti-inflammatory, antiviral, and antidiabetic. So, bioactive compounds isolated from medicinal plants have attracted significant attention because there are safe and effective against different pathogens since ancient times.[9,10] Caulerpin, the low toxic bisindole alkaloid is a more common compound of the genus Caulerpa of green macroalgae and it was isolated from Caulerpa racemosa,[11] red alga Chondria armata,[12] and brown alga Sargassum platycarpum.[13] Caulerpin has various pharmacological properties and effective biological activities such as antitumor, anti-diabetic, anticancer, anti-larvicidal, anticorrosion, antitubercular, antimicrobial, spasmolytic, antinociceptive, plant growth regulatorory activity, and anti-inflammatory.[13,14] Caulerpin exhibited antiviral activities against Chikungunya

Biological activity Antimicrobial,[22] anticancer,[23] and antiviral activities[24]

Antitrypanosomal[25]

Anti-obesity,[28] a novel selective LXRb agonist,[29] and anticancer activities[30] Antimicrobial,[36] antiviral activity against hepatitis B virus,[37] antioxidant and anticancer activities[38] Accumulation inhibition,[39] antitumor, antiviral, anti-inflammatory and anticancr activities[30,43] Antibacterial effect against E. coli and Bacillus subtilis[45] Antioxidant and a cell protective effect on a monkey kidney fibroblast cell line,[47] and anticancer, antibacterial and antifungal activities[48] Antimicrobial, anti-inflammatory, antioxidant, hypocholesterolemia, pesticide, hemolytic and 5-Alpha reductase inhibitor,[50] and antiviral activity[51] Antitumor, anti-diabetic, anticancer, antilarvicidal, anticorrosion, anti-herpes, antitubercular, antimicrobial, cytotoxic, antiviral, spasmolytic, antinociceptive, antiAlzheimer’s disease, plant growth regulatorory activity and anti-inflammatory activities[13–17]

virus[15] and herpes simplex virus type 1.[16] Also, caulerpin showed good inhibition activity against Alzheimer’s disease.[17] Due to the low-toxicity and biological activity of all studied compounds especially caulerpin in different biological applications as summarized in Table 1. We expected that they could be used to develop formal inhibitors against COVID-19. In this study, our homology modeling and sequence alignment of 2019-nCoV main protease show that SARSCoV main protease (PDB ID: 3TNT) shares 96% similarity with 3CLpro of COVID-19. Also, homology modeling and sequence alignment of 2019-nCoV main protease show that SARS-CoV main protease (PDB ID: 2GTB) shares 96.76% similarity with 3CLpro of COVID-19 as described in this study.[6] So, we study the inhibitory effect of some bioactive compounds obtained from natural sources against SARSCoV-2-3CLpro and SARS-CoV main proteases (PDB IDs: 2GTB and 3TNT). Due to their high similarity of 2GTB and 3TNT with 3CLpro, inhibition of these receptors could be effective in 3CLpro treatment. Based on molecular docking analysis, we found that caulerpin has the highest binding affinity against 6LU7, 3TNT, and 2GTB receptors compared to other bioactive compounds and the proposed drugs. The pharmacogenic and toxicity properties of all studied compounds also are computed and confirmed that most of all studied bioactive compounds especially, caulerpin obey

JOURNAL OF ENVIRONMENTAL SCIENCE AND HEALTH, PART A

3

Figure 1. All studied bioactive compounds under investigation in this work.

Lipinski’s rule and pass ADME property. Finally, molecular dynamic (MD) simulation confirmed that the complex formed between caulerpin and 6LU7 is stable in water explicit and had no major effect on the flexibility of the protein throughout the simulations and provided a suitable basis for our study. Based on this work, we believe that caulerpin could be used to develop an antiviral molecule against COVID-19. Thus, this compound could be a potential candidate for in vitro and in vivo antiviral studies followed by clinical COVID-19 treatment. It is worth mentioning that there is not any work in literature discusses the inhibition effect of all studied compounds against COVID-19. Also, most of all related work in this field studied the inhibition effect against only COVID-19 receptors, but in this work we study the inhibition activity of all studied compounds against COVID-19 main protease and two SARS-CoV main proteases receptors. Finally, we also build homology modeling and sequence alignment on SARS-CoV-2-3CL main protease to show its similarity with other SARS-CoV receptors. While, most of all related work in this field did not study it.

Material and methods Structures and biological activity of all studied bioactive compounds In our study, we choose ten public bioactive compounds that already are abundant and isolated from different species of plants or marine algae to study the inhibitory effect of them against main proteases of CoV-2 (3CLpro) with (PDB ID: 6LU7) and main proteases of SARS-CoV (Mpro) (PDB IDs: 2GTB and 3TNT). According to literature survey, these compounds showed a broad spectrum of biological activities like antimicrobial, antioxidant, anticancer, antiviral, and anti-inflammatory that making them an attractive target to evaluate their potential to become potential candidate

inhibitors against SARS-CoV-2-3CLpro. All studied bioactive compounds in this work are shown in Figure 1. Also, their sources, species names and biological activities of all studied compounds are illustrated in Table 1. Preparation of the SARS-Cov-2-3CLpro structure and inhibitors compounds The crystal structures of SARS-CoV-2-3CLpro (PDB code: 6LU7) and main proteases of SARS-Coronavirus (Mpro) with (PDB IDs: 2GTB and 3TNT) were downloaded from the Protein Data Bank (www.pdb.org), and any heteroatoms and water molecules were removed before molecular docking studies. The 3-dimensional (3D) structures of all studied compounds and some proposed antiviral drugs (bioactive compounds, chloroquine, hydroxychloroquine, azithromycin, remdesivir, baloxvir, lopinavir, and favipiravir) were obtained from PubChem (https://pubchem.ncbi.nlm.nih.gov/), in .sdf format. PubChem is a chemical substance and biological activities repository consisting of three databases, including substance, compound, and bioassay databases. Determination of the amino acids in the active site of a protein was determined using the Biovia Discovery Studio software[52] to analyze the Grid box and docking evaluation results. Homology modeling and sequence alignment of 2019-nCoV main protease Homology modeling is performed for the 2019-nCoV main protease. The crystal structure of SARS-CoV-2-3CLpro was used as a template (PDB ID: 6LU7). Therefore, sequence similarity searches were performed by using BLASTp analysis which showed suitable templates for the homology modeling. Target sequences of the SARS-CoV-2-3CLpro were created according to the build homology models

4

D. A. ABDELRHEEM ET AL.

Figure 2. The alignment sequence between SARS-CoV-2-3CLpro (PDB ID: 6LU7) and SARS sister (PDB ID: 3TNT). The a helices are shown as red blocks. The b sheets are shown as blue blocks.

Figure 3. The final Sequence alignment model between SARS-CoV-2-3CLpro (PDB ID: 6LU7) and SARS sister (PDB ID: 3TNT).

implicit in the Discovery Studio software as described in this study.[53] Molecular docking All Ligands in this study (bioactive compounds, chloroquine, hydroxychloroquine, azithromycin, remdesivir, baloxvir, lopinavir, and favipiravir) were optimized before docking by Avogadro version 1.2 with Force Field type MMFF94 and saved in .pdb format. Discovery Studio used for protein optimization, by removing water and other atoms to prepare protein for the docking process. Molecular docking between ligands and 6LU7 was carried out using Auto Dock Tools (ADT) graphical user interface supported by MGL tools.

Polar hydrogen was added and atomic charges were computed by Kollman and Gasteiger method. We defined a grid size with 60 Å
60 Å
60 Å for two receptors and the Lamarckian Genetic Algorithm (LGA) was assigned to carry out the molecular docking process, as described in this study.[54] The output of the Auto Dock was further analyzed with discovery studio program and Pymol version 1.7.4.5 software package. Analysis of drug likeness and ADMET properties of all studied bioactive compounds The drug likeness prediction of all studied bioactive compounds was carried out by Lipinski filter (http://www.scfbio-

JOURNAL OF ENVIRONMENTAL SCIENCE AND HEALTH, PART A

5

Figure 4. The Ramachandran plot of 6LU7.

Table 2. Protein target structures of 6LU7 and 3TNT and active site amino acids. PDB ID

Macromolecule

Active site

6LU7

GLN 107, ARG 107, GLN 127, THR 169, VAL125, LYS 5, ARG 105, GLN 127 THR 111, GLY 143, SER 144, CYS 145, GLN 192, GLU 166 LYS 137, ARG 298, GLN 110, LYS 137, GLU 288

3TNT

VAL125, TYR126, GLY127, PHE140, GLN 107, ARG 107, GLN 127, THR 169, VAL125, LYS 5, ARG 105, GLN 127 THR 111, GLY 143, SER 144, CYS 145, GLN 192, GLU 166 LYS 137, ARG 298, GLN 110, LYS 137, GLU 288ARG188, GLN189, GLN192, ALA198, LYS236, TYR237, GLN273

iitd.res.in/software/drugdesign/lipinski.jsp), according to which an orally active drug should comply with a minimum of four of the five laid down criteria for drug likeness namely: molecular mass, cLogP, hydrogen donor and acceptor, and molar refractive index.[55] Furthermore, the pharmacokinetic properties like absorption, distribution, metabolism, excretion, and the toxicity of all studied compounds were predicted utilizing admetSAR database (http:// lmmd.ecust.edu.cn/admetsar1/predict).[56]

MD simulation The structure of the highest binding complex obtained from molecular docking study between compound 10 (caulerpin) and 6LU7 was prepared for MD simulation using slandered dynamic cascade implicit in discovery studio. The MD simulation of SARS-CoV-2-3CLpro- inhibitor complex was carried out at 30 ns using CHARMm force field for all atoms in the complex. The simulation started by solvating the complex in triclinic box using TIP3P water model. The

6

D. A. ABDELRHEEM ET AL.

Table 3. Molecular docking analysis of studied compounds (1–10) and some antiviral drugs against 6LU7. Ligand

Binding energy (DG) k.cal/mol

Ligand efficiency

Inhibition constant (mM)

Intermolecular energy

VDW-H bond desolvation energy

–7.24 –7.47 –7.55 –8.02 –8.03 –9.22 –9.02 –8.10 –9.25 –9.28 –8.95 –9.23 –8.55 –9.24 –9.18 –9.26 –9.23

–0.26 –0.19 –0.32 –0.25 –0.29 –0.18 –0.24 –0.40 –0.28 –0.20 –0.20 –0.23 –0.22 –0.24 –0.18 –0.29 –0.27

940 820 765 523 420 160 220 427 315 189 525 418 330 380 128 320 515

–11.03 –11.21 –11.43 –12.54 –13.42 –15.65 –14.34 –12.13 –14.65 –16.32 –16.43 –15.32 –15.54 –16.32 –14.65 –13.62 –17.43

–10.00 –9.43 –10.87 –11.65 –12.12 –14.43 –13.54 –11.31 –12.88 –15.43 –15.22 –14.43 –14.13 –15.43 –13.22 –12.86 –16.25

1 2 3 4 5 6 7 8 9 10 Chloroquine Hydroxychloroquine Azithromycin Remdesivir Baloxvir Lopinavir Favipiravir

counter ions were added to the neutralized the system. In which caulerpin–6LU7 complex with solvated in 6506 water molecules and neutralized by 20 sodium and 17 chlorides as counter ions. Periodic boundary conditions were used. Throughout the simulation, the studied system is maintained at a temperature of 300 K with constant pressure. Energy minimization was done for 50,000 steps. The trajectories were collected for every nanosecond to get insights into the interactions at the atomistic level. All MD protocol was carried out according to this study.[3] The complexes result from MD simulation was analyzed for root mean square deviation (RMSD), root mean square fluctuations (RMSF) and hydrogen bonds analysis. MM-PBSA binding free energy calculations The binding free energy of protein and ligand complexes can be calculated by combining the molecular Mechanic/ Poisson–Boltzmann Surface Area (MM-PBSA) with MD. The MD scripts were extracted to perform MM-PBSA calculations. The binding free energy provides an overview of the biomolecular interactions between protein and ligand. The binding energy constitutes of potential energy, polar and non-polar solvation energies. The MM-PBSA binding free energies were calculated by utilizing nanoscale molecular dynamics (NAMD) and visual molecular dynamics (VMD) software programs according to this study.[57]

Results and discussions Homology modeling and sequence alignment of 2019nCoV main protease It is noticed that SARS-CoV and SARS-CoV-2-3CLpro share remarkable 96.00% sequence alignment among all other human coronaviruses as shown in Figure 2, based on our homology modeling and sequence alignment of 2019-nCoV main protease. The crystal structure of SARS-CoV-2-3CLpro (PDB ID: 6LU7) is highly similar to its SARS-CoV sister (PDB ID: 3TNT) at high resolution with 1.59 Å. multiple additional sequencing studies have been performed for

SARS-CoV-2, including a landmark preprint, which suggested renaming 2019-nCoV to SARS-CoV-2 on based on results similar to ours.[58] The highly conserved region of the structural elements was found, the least PDF total energy of 6LU7 and 3TNT are 2728.98 and 2881.43, respectively, that is, a reliable statistical potential to assess the quality of homology models in protein structure prediction. Also, DOPE score of 6LU7 and 3TNT are 70754.34 and 70644.22, respectively. SARS-CoV-2 and SARS-Cov-3CLpro have nine a-helices and 13 b-strands which make up three distinct domains, i.e., domain I, domain II, and domain III. All CoV proteases family consist of three domain, in which Domains I (residues 8–101) and II (residues 102–184) have one antiparallel b-barrel, that resemble the trypsin-like serine proteases structure. Domain III (residues 201–306) consists of 5 a-helices (a5-a9), that are connected by a long loop (residues185–200) with domain II as reported in this study.[6] There are only four mutated residues (Phe 13, Asn 65, Ala 94, and Val 35) between SARS-CoV-2 and SARS-CoV-23CLpro as shown in Figure 3. The Ramachandran plot built also, by discovery studio shows that 100% of the residues in the allowed regions, 97.0% in the most favored region as shown in Figure 4. Additionally, 88.2% of the residues have averaged 3D–1D score  0.3 based on the Verify 3D software, while the overall quality factor of ERRAT is 96.0%. Our homology modeling and sequence alignment of 2019nCoV main protease is the best agreement with this study,[58,59] which confirmed that the structure of SARSCoV (PDB Code 3TNT) and SARS-CoV-2-3CLpro (PDB Code 6LU7) share remarkable 96.75% sequence alignment. Due to the high sequence alignment between 3TNT and 2GTB with 6LU7, 3TNT, and 2GTB could be a potential drug target for 2019-nCoV and the inhibition of 3TNT and 2GTB protease would help to restrict the viral maturation thereby decreasing the SARS-CoV-2 infection in humans. Molecular docking The amino acids found in the active site pockets of 6LU7, 3TNT, and 2GTB and SARS-CoV-2, were elucidated from

JOURNAL OF ENVIRONMENTAL SCIENCE AND HEALTH, PART A

7

Figure 5. The best docking pose of the studied compounds generated using Pymol software. The SARS-CoV-2-3CLpro is shown in cartoon representation while the studied compounds are shown as red color surrounding by all sticks and byres of 6LU7 within 5 Å of compound.

PDB files by discovery studio and summarized in Table 2. Table 3 displays the molecular docking analysis of all studied bioactive compounds and some proposed antiviral drugs against 6LU7, including binding energy, ligand efficiency, inhibition constant, intermolecular energy, and van der Waals (VDW)-H Bond desolvation energy. Figure 5 is shown the best docking poses of the studied bioactive compounds (1–10) inside 6LU7. Also, Figure 6 displays docking analysis visualization of 6LU7 binding with all studied compounds in which the yellow dots show H-bonds formed between bioactive compounds and 6LU7 residues. The Number of H-bonds, H-bonding Residues and H-bonding distance produced from docking for all studied bioactive

compounds against 6LU7 shown in Table 4. Chloroquine, hydroxychloroquine, azithromycin, remdesivir, baloxvir, lopinavir, and favipiravir also docked inside 6LU7 using the same protocol applied on all studied compounds. An in silico analysis study showed that only compound 10 can inhibit 6LU7 which have the higher binding energies and inhibition constants and higher Number of H-bonds with 6LU7 amino acid residues with compared to chloroquine, hydroxychloroquine, azithromycin, remdesivir, baloxvir, lopinavir, and favipiravir inhibitors as shown in Tables 3 and 4. Compound 10 exhibits the highest binding energy (–9.28 k.cal/mol) with compared to other compounds and formed three hydrogen bond interactions with LYS 137, GLU 288, and LYS 5 of 6LU7 amino

8

D. A. ABDELRHEEM ET AL.

Figure 6. Docking analysis visualization of 6LU7 binding with all studied compounds generated using Pymol software. The yellow dots show H-bonds.

acid residues with distance 2.25, 2.27, and 2.38 Å, respectively. Also, all studied bioactive compounds are docked inside 3TNT and 2GTB and the obtained binding energy are summarized in Table 5. Also, compound 10 shows the highest binding energy compared to all studied bioactive compounds and antiviral drugs. So, compound 10 may act as a potential inhibitor of SARS-CoV-2-3CLpro and SARS-CoV-2. Analysis of drug likeness and ADMET properties of all studied bioactive compounds Lipinski’s rule of five is commonly utilized in development and drug design to expect oral bioavailability of drug

molecules. Lipinski’s rule was established based on five rules to compute the ability of the compound to act as an orally active drug was calculated and shown in Table 6. So, orally active drugs must have no more than one violation of the following standards: (i) octanol/water partition coefficient (log P) which measured the lipophilicity of a molecule must be not greater than five. (ii) A molecular weight (MW) less than 500 Da. (iii) not more than five hydrogen bond donors (nON). (iv) not more than 10 hydrogen bond acceptors (nOHN). The topological polar surface area (TPSA) is measured the bioavailability of the drug molecule. TPSA is closely related to the hydrogen bonding potential of a compound. TPSA of studied compounds was noticed in the

JOURNAL OF ENVIRONMENTAL SCIENCE AND HEALTH, PART A

9

Table 4. The Number of H-bonds, H-bonding Residues, and H-bonding distance produced from docking for all studied bioactive compounds against 6LU7. Compounds

N. of H-bonds

1

2

2

2

3

2

4 5

1 2

6

3

7

3

8 9

1 3

10

3

H-bonding residues

H-bonding distance (Å)

GLN 107 ARG 107 GLN 127 THR 169 VAL125 LYS 5 ARG 105 GLN 127 THR 111 GLY 143 SER 144 CYS 145 GLN 192 GLU 166 CYS 145 LYS 137 THR 111 ARG 298 GLN 110 LYS 137 GLU 288 LYS 5

2.24 2.09 2.32 2.33 2.24 2.35 2.41 2.17 2.36 2.35 2.39 2.29 2.27 2.30 2.33 2.35 2.30 2.41 2.33 2.25 2.27 2.38

Table 5. Molecular docking analysis of studied compounds (1–10) and some antiviral drugs against 3TNT and 2GTB. Ligand

3TNT Binding energy (DG) k.cal/mol

2GTB Binding energy (DG) k.cal/mol

–6.14 –7.22 –7.25 –7.88 –7.82 –8.13 –8.32 –8.44 –9.55 –9.62 –9.12 –9.34 –9.24 –9.53 –9.45 –9.50 –9.58

–7.24 –7.37 –7.44 –7.65 –7.64 –8.09 –8.14 –8.37 –9.44 –9.77 –9.32 –9.38 –9.42 –9.67 –9.51 –9.56 –9.73

1 2 3 4 5 6 7 8 9 10 Chloroquine Hydroxychloroquine Azithromycin Remdesivir Baloxvir Lopinavir favipiravir

Table 6. Prediction of molecular properties descriptors of all studied bioactive compounds. Bioactive compounds 1 2 3 4 5 6 7 8 9 10

Mass

Hydrogen bond donor

Hydrogen bond acceptors

cLogP

Molar refractivity

TPSA

282.00 428.00 428.00 368.00 750.00 186.00 196.00 256.00 712.00 398.00

1 2 2 0 4 1 1 1 4 2

2 2 2 1 10 5 3 2 10 4

6.10 7.11 7.16 0 9.20 –0.60 –1.40 5.55 3.39 3.94

87.08 135.68 138.78 0 207.42 41.68 51.60 77.94 187.92 114.13

32.41 150.55 21.26 35.49 38.42 46.00 150.88 36.19 25.87 36.76

Table 7. Prediction of ADMET descriptors of all studied bioactive compounds. Compounds 1 2 3 4 5 6 7 8 9 10

BBB

HIA

Caco-2 permeability

CYP inhibitory promiscuity

AMES toxicity

Carcinogenicity

Rat acute toxicity LD50 mol/kg

BBBþ BBBþ BBBþ BBBþ BBBBBBþ BBBþ BBBþ BBBþ BBBþ

HIAþ HIAþ HIAþ HIAþ HIAHIAþ HIAþ HIAþ HIAþ HIAþ

Caco21 Caco21 Caco21 Caco21 Caco2 Caco21 Caco21 Caco21 Caco21 Caco21

Low Low Low Low Low Low Low Low Low Low

Non-toxic Non-toxic Non-toxic Non-toxic Non-toxic Non-toxic Non-toxic Non-toxic Non-toxic Non-toxic

Non-carcinogenic Non-carcinogenic Non-carcinogenic Non-carcinogenic Non-carcinogenic Non-carcinogenic Non-carcinogenic Non-carcinogenic Non-carcinogenic Non-carcinogenic

1.95 2.29 2.38 2.60 2.17 2.18 2.21 2.75 2.70 2.23

10

D. A. ABDELRHEEM ET AL.

Figure 7. MD snapdhot of solvated compound10-6LU7 system obtained at 30 nm.

range of 25.87–153.50 Å and is well below the limit of 160 Å. It can be predicted that all studied bioactive compounds obeyed Lipinski’s rule of five and are likely to be orally active except compounds 5 and 9 as shown in Table 6. The database supports ADMET profiles which involve some features to study the ability of the studied compounds to act as drug leads such as Blood–brain barrier (BBB) penetration, human intestinal absorption (HIA), Caco-2 cell permeability, CYP inhibitory promiscuity, AMES toxicity, carcinogenicity, and rat acute toxicity LD50 are calculated and displayed in Table 7. As shown in Table 6, all studied compounds may cross blood brain barrier (BBB) and absorb in the human intestine (HIA) along are permeable for Cacoe2 cells, whereas, compound 5 showed a negative result for BBB, HIA, and Cacoe2 cell permeability. Cytochrome P450 (CYP) is a group of isozymes containing the

metabolism of drugs, steroids, fatty acids, bile acids, and carcinogens. The results indicate that these studied compounds are non-substrate and non-inhibitor of CYP enzymes.[13,60] In terms of AMES toxicity, all studied compounds were observed to be non-toxic. Carcinogenicity model indicated non-Carcinogenic nature of all studied compounds. Rat Acute Toxicity LD50 of all studied compounds was found between 1.95 and 2.75 mol/kg. The finding strongly provides the ability of most of all studied compounds to act as a drug, except compound 5 as shown in Table 7.

MD simulation To confirm the docking results and get more insight into the stability of the ligand–protein complex, MD simulations

JOURNAL OF ENVIRONMENTAL SCIENCE AND HEALTH, PART A

11

Figure 8. Plot of root mean square deviation (RMSD) of C–Ca–N backbone vs. simulation time for solvated 6LU7 in the absence and presence of compound 10 during 30 ns molecular dynamic simulations.

Figure 9. The root mean square fluctuation (RMSF) values of 6LU7 in the absence and presence of compound 10 at 30 nm were plotted against residue numbers.

were carried out for the compound which has the highest binding affinity (compound 10) inside 6LU7 complex in the solvated states at 30 ns as shown in Figure 7. The results of MD simulations have been examined based on RMSD, RMSF, and the number of hydrogen bonds as a function of time.

receptors (6LU7) and the complexes with a low average RMSF. This trend in the RMSF plot for the complex indicates that binding caulerpin to the receptor was stable and had no major effect on the flexibility of the protein throughout the simulations. To explore more insights on the local protein flexibility, the time average of RMSF values of the 300 amino acids of 6LU7 in the absence and presence of the caulerpin over simulation period were calculated. The RMSF values for the three complexes suggested that the following residues (ALA206, VAL204, and LEU205) showed less fluctuation in the presence and presence of the caulerpin. The average RMSF values of these residues were 0.33 ± 0.06, 0.40 ± 0.07, and 0.43 ± 0.08 Å for 6LU7, and 0.37 ± 0.07, 0.48 ± 0.08, and 0.52 ± 0.09 Å for caulerpin–6LU7 complex.

Root mean square deviation To examine the change in the protein dynamics and the conformational stability of the protein–ligand complexes, the protein complexed with compound 10 was subjected to 30 ns MD simulations. Standard dynamics cascade module implicit in DS software was employed to measure the RMSD and RMSF. The RMSD measures the direct changes in the protein from the initial coordinates. The RMSD values of the protein backbone in complex with the potential inhibitors were computed concerning the initial structure as a frame reference (0 to 10 ns). The RMSD values steadily increased from 0 to 5 ns, and reached equilibration after that throughout the simulation period. The RMSD values for the studied complex showed oscillations between 3 to 5 ns indicating that the studied compound was adapting another conformation within the binding pocket as shown in Figure 8. The average RMSD values for the last 1 ns for the 6LU7 and caulerpin þ 6LU7 complex, were 1.13 ± 0.07 and 1. 22 ± 0.05, respectively. Lower RMSD value of the caulerpin in the complex with 6LU7 indicates its stability and provided a suitable basis for our study. Root mean square fluctuation RMSF was measured based on the backbone atom of each amino acid residues and the plot of RMSF was used to depict the fluctuations at the residue level. RMSF plot as shown in Figure 9 of solvated 6LU7 and 6LU7 in complex with caulerpin during 30 ns MDs simulations exhibited a similar trend of residue fluctuation profile for both free

Hydrogen bond analysis and MM-PBSA binding free energy Hydrogen bonds are one of the essential elements answerable for the molecular interactions in biological systems. Hydrogen bonds provide the basis for molecular recognition and selectivity by imparting directionality and explicitness to molecular interactions. The protein–ligand interactions were guided by the changes in the secondary structures, which in turn were regulated by the hydrogen bonds. MD simulation provided different conformations in which a protein could be found in actual biological conditions. Each conformation of a protein is supposed to have its interaction pattern with the ligand. We calculated the number of hydrogen bonds formed during the complete run of MD simulations for selected complexes, as presented in Figure 10. As shown in Figure 10, in the caulerpin–6LU7 complex, the most number of conformations formed up to 10 hydrogen bonds during the simulation. A very few conformations showed less than 5 and greater than 10 hydrogen bonds. The bioactive caulerpin molecule was able to maintain strong interaction with the binding pocket of 6LU7 throughout the simulation period. We measured binding free energy at last 10 ns from MD simulation. The final binding energy is a cumulative

12

D. A. ABDELRHEEM ET AL.

Figure 10. Hydrogen bond profiles of 6LU7-Compound 10 complex at 30 nm.

sum of van der Wall, electrostatic, polar solvation, and SASA energy. Except for the polar solvation energy, all other forms of energy contributed favorably to the interaction between caulerpin molecule and 6LU7. The caulerpin molecule showed the binding free energy equal to 256.44 ± 27.55 kJ/mol, SASA equal to 28.58 ± 4.32 kJ/mol, DEpolar solvation equal to 283.58 ± 44.53 kJ/mol, DEElectrostatic equal to 127.76 ± 40.50 kJ/mol, and DEVan der Waal equal to 384.44 ± 34.40 kJ/mol inside 6LU7. This confirmed the efficacy of the caulerpin molecule against SARS-CoV-2 main protease.

to maintain strong interaction with the binding pocket of 6LU7 throughout the simulation period. Finally, the caulerpin molecule showed the binding free energy inside 6LU7 equal to 256.44 ± 27.55 kJ/mol which confirmed the efficacy of this molecule against SARS-CoV-2 main protease. Hence, caulerpin is highly effective against SARS-CoV-2 main protease and can be explored further for drug repurposing against the successful inhibition of COVID-19.

ORCID Hussein S. Mohamed

https://orcid.org/0000-0001-6524-9649

Conclusion Recently, the rapidly spreading outbreak of COVID-19 has challenged the healthcare sector of the world. Numerous antiviral or other conventional drugs are being examined against COVID-19, but not yet positive results have come. To contribute to this fight against COVID-19, virtual screening-based molecular docking and MD simulation were investigated to identify novel and potential to inhibit SARSCoV-2-3CLpro main protease. Based on our molecular docking analysis we found that among all studied compounds, caulerpin has the highest binding affinity against all studied receptors 6LU7, 3TNT, and 2GTB with compared to some proposed antiviral drug currently used in COVID-19 treatment. The ADMET and molecular descriptors properties strongly provide that most of all studied compounds obey Lipinski’s rule of five and exhibit significant biological activities especially, caulerpin compound. MD simulation confirmed the stability of caulerpin-6LU7 receptor in the presence of explicit water solvent. RMSD of the receptor–ligands complex has maintained the stability at around 1.8 Å and RMSD of caulerpin complexed with the protein are in the favorable range of within 1.2 Å and has remained stable during the simulations. The backbone atoms of the complex and free receptor show similar RMSF, indicating the stability of the caulerpin inside 6LU7 receptor. Also, the most number of conformations formed up to 10 hydrogen bonds during the simulation in caulerpin–6LU7 complex. A very few conformations showed less than 5 and greater than 10 hydrogen bonds. The bioactive caulerpin molecule was able

References [1]

[2]

[3]

[4]

[5]

[6]

[7]

[8]

Li, Y.; Zhang, J.; Wang, N.; Li, H.; Shi, Y.; Guo, G.; Zou, Q. Therapeutic Drugs Targeting 2019-nCoV Main Protease by High-Throughput Screening. bioRxiv 2020. Zhu, N.; Zhang, D.; Wang, W.; Li, X.; Yang, B.; Song, J.; Zhao, X.; Huang, B.; Shi, W.; Lu, R.; et al. A Novel Coronavirus from Patients with Pneumonia in China, 2019. N. Engl. J. Med. 2020, 382, 727–733. Xu, X.; Chen, P.; Wang, J.; Feng, J.; Zhou, H.; Li, X.; Zhong, W.; Hao, P. Evolution of the Novel Coronavirus from the Ongoing Wuhan Outbreak and Modeling of Its Spike Protein for Risk of Human Transmission. Sci. China Life Sci. 2020, 63, 457–460. DOI: 10.1007/s11427-020-1637-5. Xia, B.; Kang, X. Activation and Maturation of SARS-CoV Main Protease. Protein Cell 2011, 2, 282–290. DOI: 10.1007/ s13238-011-1034-1. Lu, I.-L.; Mahindroo, N.; Liang, P.-H.; Peng, Y.-H.; Kuo, C.-J.; Tsai, K.-C.; Hsieh, H.-P.; Chao, Y.-S.; Wu, S.-Y. StructureBased Drug Design and Structural Biology Study of Novel Nonpeptide Inhibitors of Severe Acute Respiratory Syndrome Coronavirus Main Protease. J. Med. Chem. 2006, 49, 5154–5161. DOI: 10.1021/jm060207o. Zhang, H.; Saravanan, K. M.; Yang, Y.; Hossain, M. T.; Li, J.; Ren, X.; Wei, Y. Deep Learning Based Drug Screening for Novel Coronavirus 2019-nCov. Preprints. 2020. Baell, J. B. Feeling Nature’s PAINS: Natural Products, Natural Product Drugs, and Pan Assay Interference Compounds (PAINS). J. Nat. Prod. 2016, 79, 616–628. DOI: 10.1021/acs. jnatprod.5b00947. Malarvizhi, D.; Anusooriya, P.; Meenakshi, P.; Sundaram, S.; Oirere, E.; Gopalakrishnan, V. K. Isolation, Structural Characterization of Oleic Acid from Zaleya decandra Root

JOURNAL OF ENVIRONMENTAL SCIENCE AND HEALTH, PART A

[9] [10]

[11]

[12] [13]

[14]

[15]

[16]

[17]

[18]

[19]

[20]

[21]

[22]

[23]

[24]

Extract. Anal. Chem. Lett. 2016, 6, 669–677. DOI: 10.1080/ 22297928.2016.1238319. Boxall, A. B. A. The Environmental Side Effects of Medication. EMBO Rep. 2004, 5, 1110–1116. Angiola, L.; Sacchetti, G.; Efferth, T. Antimicrobial and Antioxidant Activities of Natural Compounds. Evid-Based. Complem. Altern. Med. 2018, 2018, 1–3. Nagappan, T.; Vairappan, C. S. Nutritional and Bioactive Properties of Three Edible Species of Green Algae, Genus Caulerpa (Caulerpaceae). J. Appl. Phycol. 2014, 26, 1019–1027. DOI: 10.1007/s10811-013-0147-8. Mangala, S. W.; Govenkar, B. Constituents of Chondria Armata. Phytochemistry 2000, 54, 979–981. Abdelrheem, D. A.; Abd El-Mageed, H. R.; Mohamed, H. S.; Rahman, A. A.; Elsayed, K. N.; Ahmed, S. A. Bis-Indole Alkaloid Caulerpin from a New Source Sargassum platycarpum: Isolation, Characterization, In Vitro Anticancer Activity, Binding with Nucleobases by DFT Calculations and MD Simulation. J. Biomol. Struct. Dyn. 2020, 1–11. DOI: 10.1080/ 07391102.2020.1784285. Lunagariya, J.; Bhadja, P.; Zhong, S.; Vekariya, R.; Xu, S. Marine Natural Product Bis-Indole Alkaloid Caulerpin: Chemistry and Biology. MRMC. 2019, 19, 751–761. DOI: 10. 2174/1389557517666170927154231. Esteves, P. O.; de Oliveira, M. C.; de Souza Barros, C.; CirneSantos, C. C.; Laneuvlille, V. T.; Palmer Paix~ao, I. C. Antiviral Effect of Caulerpin against Chikungunya. Nat. Prod. Commun. 2019, 14, 1934578X19878295. Macedo, N. R. P. V.; Ribeiro, M. S.; Villac¸a, R. C.; Ferreira, W.; Pinto, A. M.; Teixeira, V. L.; Cirne-Santos, C.; Paix~ao, I. C. N. P.; Giongo, V. Caulerpin as a Potential Antiviral Drug against Herpes Simplex Virus Type 1. Rev. Bras. Farmacogn. 2012, 22, 861–867. DOI: 10.1590/S0102-695X2012005000072. Lorenzo, V. P.; Barbosa Filho, J. M.; Scotti, L.; Scotti, M. T. Combined Structure-and Ligand-Based Virtual Screening to Evaluate Caulerpin Analogs with Potential Inhibitory Activity against Monoamine Oxidase B. Rev. Bras. Farmacogn. Braz. J. Pharmacogn. 2015, 25, 690–697. DOI: 10.1016/j.bjp.2015.08.005. Horincar, V. B.; Parfene, G.; Tyagi, A. K.; Gottardi, D.; Dinica, R.; Guerzoni, M. E.; Bahrim, G. Extraction and Characterization of Volatile Compounds and fatty acids from Red and Green Macroalgae from the Romanian Black Sea in Order to Obtain Valuable Bioadditives and Biopreservatives. J. Appl. Phycol. 2014, 26, 551–559. DOI: 10.1007/s10811-013-0053-0. Otero, P.; Quintana, S. E.; Reglero, G.; Fornari, T.; GarcıaRisco, M. R. Pressurized Liquid Extraction (PLE) as an Innovative Green Technology for the Effective Enrichment of Galician Algae Extracts with High Quality Fatty Acids and Antimicrobial and Antioxidant Properties. Mar. Drugs 2018, 16, 156. DOI: 10.3390/md16050156. Mahmoud, E. A.; Farahat, L. A.; Aziz, Z. K. A.; Fatthallah, N. A. R. A. S.; El Din, R. A. S. Evaluation of the Potential for Some Isolated Microalgae to Produce Biodiesel. Egypt J. Pet. 2015, 24, 97–101. DOI: 10.1016/j.ejpe.2015.02.010. Yang, Y.; Jin, Z.; Jin, Q.; Dong, M. Isolation and Fatty Acid Analysis of Lipid-Producing Endophytic Fungi from Wild Chinese Torreya grandis. Microbiology 2015, 84, 710–716. DOI: 10.1134/S0026261715050173. € ¸en, O.; € Kiran, I.; € Dag, I. _ Ali, O.; _ Antimicrobial Evaluation of Ozs Oleic Acid Biotransformation Metabolites. International Congress on Applied Biological Sciences, 2015. Mericli, F.; Becer, E.; Kabadayı, H.; Hanoglu, A.; Yigit Hanoglu, D.; Ozkum Yavuz, D.; Ozek, T.; Vatansever, S. Fatty Acid Composition and Anticancer Activity in Colon Carcinoma Cell Lines of Prunus dulcis Seed Oil. Pharm. Biol. 2017, 55, 1239–1248. DOI: 10.1080/13880209.2017.1296003. Sands, J. A. Inactivation and Inhibition of Replication of the Enveloped Bacteriophage u6 by Fatty Acids. Antimicrob. Agents. Chemother. 1977, 12, 523–528. DOI: 10.1128/AAC.12.4. 523.

[25]

[26]

[27]

[28]

[29]

[30]

[31]

[32]

[33]

[34]

[35]

[36]

[37]

[38]

[39]

[40]

13

Hoet, S.; Pieters, L.; Muccioli, G. G.; Habib-Jiwan, J. L.; Opperdoes, F. R.; Quetin-Leclercq, J. Antitrypanosomal Activity of Triterpenoids and Sterols from the Leaves of Strychnos spinosa and Related Compounds. J. Nat. Prod. 2007, 70, 1360–1363. DOI: 10.1021/np070038q. Shi, D.; Guo, S.; Fan, X. A New Ketosteroid from Red Alga Acanthophora spicifera. Chin. J. Ocean. Limnol. 2011, 29, 674–678. DOI: 10.1007/s00343-011-0183-7. Huang, X.; Li, M.; Xu, B.; Zhu, X.; Deng, Z.; Lin, W. Proteasome and NF-KappaB inhibiting phaeophytins from the green alga Cladophora fascicularis. Molecules 2007, 12, 582–592. DOI: 10.3390/12030582. Lee, J. A.; Cho, Y.; Hong, S. S.; Ahn, E. Anti-Obesity Activity of Saringosterol Isolated from Sargassum muticum (Yendo) Fensholt Extract in 3T3-L1 Cells . Phytother. Res. 2017, 31, 1694–1701. DOI: 10.1002/ptr.5892. Chen, Z.; Liu, J.; Fu, Z.; Ye, C.; Zhang, R.; Song, Y.; Zhang, Y.; Li, H.; Ying, H.; Liu, H. 24(S)-Saringosterol from Edible Marine Seaweed Sargassum fusiforme Is a Novel Selective LXRb Agonist . J. Agric. Food Chem. 2014, 62, 6130–6137. DOI: 10. 1021/jf500083r. Ahmed, S. A.; Rahman, A. A.; Elsayed, K. N.; Abd El-Mageed, H. R.; Mohamed, H. S.; Ahmed, S. A. Cytotoxic Activity, Molecular Docking, Pharmacokinetic Properties and Quantum Mechanics Calculations of the Brown Macroalga Cystoseira trinodis Compounds. J. Biomol. Struct. Dyn. 2020, 1–31. Nyigo, V. A.; Peter, X.; Mabiki, F.; Malebo, H. M.; Mdegela, R. H.; Fouche, G. Isolation and Identification of Euphol and b-Sitosterol from the Dichloromethane Extracts of Synadenium glaucescens. J. Phytopharm. 2016, 5, 100–104. Fasya, A. G.; Baderos, A.; Madjid, A. D. R.; Amalia, S.; Megawati, D. S. Isolation, Identification and Bioactivity of Steroids Compounds from Red Algae Eucheuma cottonii Petroleum Ether Fraction. AIP Conf. Proc. 2019, 2120, 30025. El-Sherif, M. F.; S, S. H.; M, G. Studies on Fatty Acids and Sterols Content of Some Chlorophycean Members. Egypt J. Phycol. 2004, 5, 35–50. Payghami, N.; Jamili, S.; Rustaiyan, A.; Saeidnia, S.; Nikan, M.; Gohari, A. R. Alpha-Amylase Inhibitory Activity and Sterol Composition of the Marine Algae, Sargassum glaucescens. Pharmacogn. Res. 2015, 7, 314. Fithriani, D.; Ambarwaty, D. Identification of Bioactive Compounds from Nannochloropsis sp. IOP Conf. Ser.: Earth Environ. Sci. 2020, 404, 12064. Ododo, M. M.; Choudhury, M. K.; Dekebo, A. H. Structure Elucidation of b-Sitosterol with Antibacterial Activity from the Root Bark of Malva parviflora. Springerplus 2016, 5, 1210DOI: 10.1186/s40064-016-2894-x. Parvez, M. K.; Alam, P.; Arbab, A. H.; Al-Dosari, M. S.; Alhowiriny, T. A.; Alqasoumi, S. I. Analysis of Antioxidative and Antiviral Biomarkers b-Amyrin, b-Sitosterol, Lupeol, Ursolic Acid in Guiera senegalensis Leaves Extract by Validated HPTLC Methods. Saudi Pharm. J. 2018, 26, 685–693. DOI: 10. 1016/j.jsps.2018.02.022. Shahat, A. A.; Hidayathulla, S.; Khan, A. A.; Alanazi, A. M.; Al Meanazel, O. T.; Alqahtani, A. S.; Alsaid, M. S.; Hussein, A. A. Phytochemical Profiling, Antioxidant and Anticancer Activities of Gastrocotyle hispida Growing in Saudi Arabia. Acta Trop. 2019, 191, 243–247. DOI: 10.1016/j.actatropica.2019.01.013. Shiota, A.; Hada, T.; Baba, T.; Sato, M.; Yamanaka-Okumura, H.; Yamamoto, H.; Taketani, Y.; Takeda, E. Protective Effects of Glycoglycerolipids Extracted from Spinach on 5-Fluorouracil Induced Intestinal Mucosal Injury. J. Med. Invest. 2010, 57, 314–320. DOI: 10.2152/jmi.57.314. Illijas, M. I.; Indy, I. R.; Yasui, H.; Itabashi, Y. Lipid Class and Fatty Acid Composition of a Little-Known and Rarely Collected Alga Exophyllum wentii Weber-Van Bosse from Bali Island, Indonesia. J. Oleo Sci. 2009, 58, 103–110. DOI: 10.5650/jos.58. 103.

14

[41]

[42]

[43]

[44]

[45]

[46]

[47]

[48]

[49]

[50]

D. A. ABDELRHEEM ET AL.

Plouguerne, E.; da Gama, B. A. P.; Pereira, R. C.; BarretoBergter, E. Glycolipids from Seaweeds and Their Potential Biotechnological Applications. Front. Cell. Infect. Microbiol. 2014, 4, 174. DOI: 10.3389/fcimb.2014.00174. Marcolongo, G.; De Appolonia, F.; Venzo, A.; Berrie, C. P.; Carofiglio, T.; Ceschi Berrini, C. Diacylglycerolipids Isolated from a Thermophile Cyanobacterium from the Euganean Hot springs. Nat. Prod. Res. 2006, 20, 766–774. DOI: 10.1080/ 14786410500176393. Zhang, J.; Li, C.; Yu, G.; Guan, H. Total Synthesis and Structure-Activity Relationship of Glycoglycerolipids from Marine Organisms. Mar. Drugs 2014, 12, 3634–3659. DOI: 10. 3390/md12063634. Peng, Y.; Huang, R. M.; Lin, X. P.; Liu, Y. H. Norisoprenoids from the Brown Alga Sargassum naozhouense Tseng et Lu. Molecules 2018, 23, 348. DOI: 10.3390/molecules23020348. Baidya, M.; Griffin, K. A.; Yamamoto, H. Catalytic Enantioselective O-Nitrosocarbonyl Aldol Reaction of b-Dicarbonyl Compounds. J. Am. Chem. Soc. 2012, 134, 18566–18569. DOI: 10.1021/ja309311z. Singh, A. K.; Chaudhuri, P. K.; Ghosal, S. Chemical Constituents of Genianaceae XX: Natural Occurrence of (-)-Loliolide in Canscora decussata. J. Pharm. Sci. 1976, 65, 1549–1551. DOI: 10.1002/jps.2600651038. Yang, X.; Kang, M. C.; Lee, K. W.; Kang, S. M.; Lee, W. W.; Jeon, Y. J. Antioxidant Activity and Cell Protective Effect of Loliolide Isolated from Sargassum ringgoldianum Subsp. Coreanum. Algae 2011, 26, 201–208. DOI: 10.4490/algae.2011. 26.2.201. Grabarczyk, M.; Wi nska, K.; Ma˛czka, W.; Potaniec, B.; Anioł, M. Loliolide–the Most Ubiquitous Lactone. Biologica 2015, 11, 1–8. DOI: 10.1515/fobio-2015-0001. Krishnan, K. R.; James, F.; Mohan, A. Isolation and Characterization of n-Hexadecanoic Acid from Canthium parviflorum Leaves. J. Chem. Pharm. Res. 2016, 8, 614–617. Abubakar, M.; Majinda, R. GC-MS Analysis and Preliminary Antimicrobial Activity of Albizia adianthifolia (Schumach) and Pterocarpus angolensis (DC). Medicines 2016, 3, 3. DOI: 10. 3390/medicines3010003.

[51]

[52]

[53]

[54]

[55]

[56]

[57]

[58]

[59]

[60]

Libran-Perez, M.; Pereiro, P.; Figueras, A.; Novoa, B. Antiviral Activity of Palmitic Acid via Autophagic Flux Inhibition in Zebrafish (Danio rerio). Fish Shellfish Immunol. 2019, 95, 595–605. DOI: 10.1016/j.fsi.2019.10.055. Shen, M. Y.; Sali, A. Statistical Potential for Assessment and Prediction of Protein Structures. Protein Sci. 2006, 15, 2507–2524. DOI: 10.1110/ps.062416606. Chen, X. R.; Wang, X. T.; Hao, M. Q.; Zhou, Y. H.; Cui, W. Q.; Xing, X. X.; Xu, C. G.; Bai, J. W.; Li, Y. H. Homology Modeling and Virtual Screening to Discover Potent Inhibitors Targeting the Imidazole Glycerophosphate Dehydratase Protein in Staphylococcus xylosus. Front. Chem. 2017, 5, 98. DOI: 10. 3389/fchem.2017.00098. Sang, P.; Tian, S. H.; Meng, Z. H.; Yang, L. Q. Insight Derived from Molecular Docking and Molecular Dynamics Simulations into the Binding Interactions Between HIV-1 Protease Inhibitors and SARS-CoV-2 3CLpro. Preprint. 2020. Lipinski, C. A. Lead- and Drug-Like Compounds: The Rule-ofFive Revolution. Drug Discov. Today Technol. 2004, 1, 337–341. DOI: 10.1016/j.ddtec.2004.11.007. Cheng, F.; Li, W.; Zhou, Y.; Shen, J.; Wu, Z.; Liu, G.; Lee, P. W.; Tang, Y. admetSAR: A Comprehensive Source and Free Tool for Assessment of Chemical ADMET Properties. J. Chem. Inf. Model. 2012, 52, 3099–3105. DOI: 10.1021/ci300367a PMID: 23092397. Bhardwaj, V. K.; Singh, R.; Sharma, J.; Rajendran, V.; Purohit, R.; Kumar, S. Identification of Bioactive Molecules from Tea Plant as SARS-CoV-2 Main Protease Inhibitors. J. Biomol. Struct. Dyn. 2020, 1–13. DOI: 10.1080/07391102.2020.1766572. Hilgenfeld, R. From SARS to MERS: Crystallographic Studies on Coronaviral Proteases Enable Antiviral Drug Design. FEBS J. 2014, 281, 4085–4096. DOI: 10.1111/febs.12936. Liu, X.; Zhang, B.; Jin, Z.; Yang, H.; Rao, Z. The Crystal Structure of COVID-19 Main Protease in Complex with an Inhibitor N3. PDB. 2020. Wang, J. Fast Identification of Possible Drug Treatment of Coronavirus Disease -19 (COVID-19) through Computational Drug Repurposing Study. J. Chem. Inf. Model. 2020.