key: cord-1016854-lho0vil6
authors: El-Hawary, Seham S.; Rabeh, Mohamed A.; Raey, Mohamed A. El; El-Kadder, Essam M. Abd; Sobeh, Mansor; Abdelmohsen, Usama Ramadan; Albohy, Amgad; Andrianov, Alexander M.; Bosko, Ivan P.; Al-Sanea, Mohammad M.; El-Kolobby, Dalia G.
title: Metabolomic profiling of three Araucaria species, and their possible potential role against COVID-19
date: 2021-02-18
journal: Journal of biomolecular structure & dynamics
DOI: 10.1080/07391102.2021.1885494
sha: f47e5c08f2b7afe26f33df714d3cb0e0eecc7994
doc_id: 1016854
cord_uid: lho0vil6

The COVID-19 pandemic in Egypt is a part of the worldwide global crisis of coronavirus 2 (SARS-CoV-2). The contagious life-threatening condition causes acute respiratory syndrome. The present study aimed to assess the compounds identified by LC-MS of the methanolic leaves extracts from three conifers trees cultivated in Egypt (Araucaria bidwillii, Araucaria. cunninghamii and Araucaria heterophylla) via docking technique as potential inhibitor of COVID-19 virus on multiple targets; viral main protease (M(pro), 6LU7), non-structural protein-16 which is a methyl transferase (nsp16, 6W4H) and RNA dependent RNA polymerase (nsp12, 7BV2). Among the three targets, nsp16 was the best target recognized by the tested compounds as can be deduced from docking studies. Moreover, the methanolic extract of A. cunninghamii showed the highest radical-scavenging activity using (DPPH test) with 53.7 µg/mL comparable to ascorbic acid with IC(50) = 46 µg/mL The anti-inflammatory potential carried using enzyme linked immunoassay showed the highest activity for A. cunninghamii and A. bidwillii followed by A. heterophylla with IC(50) = 23.20 ± 1.17 µg/mL, 82.83 ± 3.21 µg/mL and 221.13 ± 6.7 µg/mL, respectively (Celecoxib was used as a standard drug with IC(50) = 141.92 ± 4.52 µg/mL). Moreover, a molecular docking study was carried for the LC-MS annotated metabolites to validate their anti-inflammatory inhibitory effect using Celecoxib as a reference compound and showed a high docking score (−7.7 kcal/mol) for Octadecyl (E) P-coumarate and (−7.3 kcal/mol) for secoisolariciresinol rhamnoside. Communicated by Ramaswamy H. Sarma

The COVID-19 pandemic in Egypt is a part of the worldwide global crisis of coronavirus 2 (SARS-CoV-2). It is a contagious condition causes acute respiratory syndrome (Chiappetta et al., 2020; Zabetakis et al., 2020) . The first case reported with that disease was a women in local market in Wuhan province, China in December 2019, the scientists believe it's acquired from animals and started spreading between humans (https://africacdc.org/covid-19/). The danger of that disease is the rapid spreading as it's highly contagious and its vigorous symptoms which might lead to death (Schett et al., 2020) . Up till now is available only a supportive treatment as a strategy in the treatment protocol, WHO welcomed the satisfactory results of using dexamethasone as a preliminary treatment for COVID-19 (World Health Organiztion, 2020) as it reduces the risk of mortality in infected patients for its role as anti-inflammatory reducing the immune-system cytokine storm reducing the exaggerated immune response caused by the COVID-19 (Mahmudpour et al., 2020) Inflammation is a biological body response as a result of exposure to stimuli which provokes the immune system to eliminate the harmful stimuli. Prostaglandins are responsible for the generation of any inflammatory response (Mart ınez-S anchez et al., 2020). Prostaglandins biosynthesis is markedly increased in inflamed tissue, arachidonic acid biosynthesis to prostaglandin H2 is catalyzed by cyclooxygenase (Lee et al., 2017) . As a result of the human defense mechanism during the inflammation, there is an excessive production of reactive oxygen species by phagocytic leucocytes (Jayawardena et al., 2020) , here appears the role of antioxidants to scavenge the reactive species, these auxiliary process aids tissue repair and healing (Ik et al., 2013) .

Plants used medicinally considered as a huge source for constituents with biological various activities. So, it was a substantial importance to seek a plant origin drug that might provide a safer therapy to control the covid-19 symptoms or as a curative treatment. Therefore, three Araucaria species were subjected for the following study to investigate their anti-inflammatory and antioxidant potential activities. Antiinflammatory activity for that genus was reported before using different techniques, exerting a significant results (Aslam et al., 2013; Elshamy et al., 2020) which was interesting to carry further investigations for the anti-inflammatory activity using different and accurate assay; as the enzyme immunoassay ELISA method which is performed here in the current research comparing three species.

It is well known now the role of inflammation in the pathogenesis of SARS-CoV-2 which suggests a promising effect of anti-inflammatory drugs as adjuvant therapy in COVID-19 patients to suppress the cytokine storm. Furthermore, dexamethasone was recently used to control SARS-CoV-2 progression (World Health Organiztion, 2020) . This information has encouraged us to virtually investigate our compounds to see if any of them could have a potential effect on the new COVID-19 virus possible targets. The active sites of three targets, viral main protease (M pro , 6LU7), nonstructural protein-16 which is a methyl transferase (nsp16, 6W4H) and RNA dependent RNA polymerase (nsp12, 7BV2) using a published recently high-resolution structure of COVID-19 main protease (Mpro ) (Jin et al., 2020; Zhang et al., 2020) allowed the chance for the creation of an inhibitor with an important role affecting the transcription and replication of virus (Aly, 2020) . The importance of Mpro in the life cycle of coronavirus in contrast it is not present as an identical human homologues, introduced Mpro as favoured antiviral targeted drug designed (Pillaiyar et al., 2016) .

A. bidwillii, A. cunninghamii and A. heterophylla leaves were collected and shed dried. The plants identity was confirmed at the Department of Flora and Phytotaxonomy, ARC, Cairo, Egypt.

The powdered leaves dried in the air of A. bidwillii, A. cunninghamii and A. heterophylla weighing 200 g for each species, was separately extracted with methanol by percolation. methanolic extracts were subjected for evaporation under reduced pressure at 45 C to a constant weight then were subjected to lyophilization and kept in the dark at 4 C until analysis.

Analytical grade reagents were purchased and used with no additional purification. Organic solvents: methanol, (Adwic, Nasr Pharma, Egypt), DPPH, standardized extract of Ascorbic acid., (Memphis Co., Egypt). (2, 2-diphenyl-1-picrylhydrazyl) free radical (Sigma Co., USA), COX-1 (human) Inhibitor Screening Assay Kit., by (Caymanchem, USA), Celecoxib (Pfizer Inc., Egypt).

Metabolic profiling was performed on methanolic extracts of A. bidwillii, A. heterophylla and A. cunninghamii according to (Sobeh et al., 2017) . HPLC-PDA-MS/MS was carried at Institute of Pharmacy and Molecular Biotechnology, Heidelberg University, Germany using a ThermoFinnigan LC system (Thermo Electron Corporation, Austin, TX, USA). A Zorbax Eclipse (XDB-C18), Rapid resolution of 4.6 Â 150 mm, 3.5 mm column was used (Agilent, Santa Clara, CA, USA).

Operating mass parameters were used in the negative mode. A combination of gradient water and acetonitrile (ACN), 0.1% formic acid have been used, through 60 min in 1 mL/min flow rate and 1:1 splitted before the ESI source, acetonitrile was raised from 5% ! 30%. The samples were injected separately using the autosampler. For MS analysis, LCQ-Duo ion trap with a ThermoQuest ESI source has been used. The Xcalibur software was used (Xcalibur TM 2.0.7, Thermo Scientific, Waltham, MA, USA) for system controlling.

The plant extract ability to quench DPPH free radicals was evaluated by a standard method (Takao et al., 1994) accredited with slight modifications (Kumarasamy et al., 2007) . Methanolic extracts of A. bidwillii, A. cunninghamii and A. heterophylla were dissolved in methanol to accomplish 15 mg/mL concentrations. Dilutions were carried to have the concentrations of 500, 250, 125, 62.5, 31.25 lg/mL. (100 lL each) of the solutions after dilution were added to 3 mL of methanolic solution of DPPH with the concentration of 0.002%. As, the standard ascorbic acid stock solution was dissolved in methanol for achieving a 1.5 mg/mL concentration, dilutions were carried to achieve concentrations 50, 25, 12.5, 6.25, 3.125 lg/mL incubation in the dark for 30 min, the absorbances were observed. The reaction undergoes between the antioxidant and free radical to produce 1,1-diphenyl-2-picryl hydrazine. The capacity to quench the free radical, DPPH was recorded at the absorbance of 517 nm.

The biosynthesis of AA to PGH2 is catalyzed by Cyclooxygenase. PGH2 produces PGF2a, via stannous chloride reduction to be assessed by enzyme immunoassay (ELISA) (Pradelles et al., 1985) . This is a competitive method among a PG tracer (PG-acetylcholinesterase conjugate) and PGs for the limited quantity of PG antiserum. The PG tracer amount that can bind to the PG antiserum is indirectly proportional to the amount of PGs in the wells, the PG tracer concentration is kept constant, but PG concentration varies. This complex of antibody-PG conjugate to an antibody (anti-IgG) which was before linked to the well. The buffer and the acetylcholinesterase substrate (Ellman's reagent) poured to the well. The reaction enzymatically produces a yellow colour which is measured via a spectrophotometer in a Microplate special Reader (BioRad, Japan) at 412 nm. The technique to have 100% COX activity was carried out in the presence and absence of DMSO as a solvent control. The inhibitory experiment was carried for extracts at different concentrations starting from 1 ending with 16 lg total phenolic compound/ mL) or of the celecoxib drug as a standard anti-inflammatory. After adding the enzyme, a pre-incubation period was allowed of 10 min before adding the inhibitor and 2 min of incubation in case of AA presence at 37 C. inactivation of COXs by keeping them in water and allowed to boil for 3 min as Enzyme control. 29 PG/mL was the detection limit and the coefficients of variations intra and interassay were 5 and 10%, respectively. The anti-inflammatory activity of the tested extracts was assessed by determining the percent of inhibition of PGE2 produced and the concentration of the test compounds leading to 50% inhibition of the release of PGE2 (IC 50 ) was detected from the curve of concentrationinhibition response by analytical regression.

Docking took place in the compounds active sites of the corresponding proteins including COX-1 (1EQG), main protease of SARS-CoV-2 (6LUV7), nsp16 (6W4H) and nsp12 (7BV2). Ligands and proteins were prepared as reported earlier (Ebada et al., 2020) . In short, structures were either drawn or used from the PubChem (www.pubchem.ncbi.nlm.nih.gov). Minimization for Structures via one thousand steps of steepest descent after that a minimization through one thousand steps of gradient conjugate. Available Proteins were used from the protein data bank (www.rcsb.org) then hydrogens were added and merged, water and non-bonded fragments were removed, and Gasteiger charges were calculated. AutoDock vina was used for carrying docking (Trott & Olson, 2010) with a grid box of 25 Â 25 Â 25 Å 3 centered on the internal ligand using exhaustiveness of 16. PyMOL was responsible for generating images (Schrodinger, 2010) .

Amber ff14SB force field was used to implement Amber18 to make the ligand-protein complexes in water. The module of antechamber was used to adjust the Gasteiger atomic partial charges (Case et al., 2020) . For preparation of the general Amber GAFF force field and the force field parameters (Wang et al., 2004) , Addition of hydrogen atoms was to target proteins via the tleap program of the AmberTools18 package (Case et al., 2020) . First, the complexes of ligand-protein were each laid in a box with cubical shape with conditions of periodic boundary, TIP3P water as a solvent was included in the box for the MD simulations of MD. (Jorgensen et al., 1983) , Cl À and Na þ ions supplied a salt of 0.10 M concentration. the system was set up, minimization of energy was carried via five hundred steps of steepest descent algorithm and five hundred steps of gradient conjugation method. The essential atoms for complex clustering were fixed by the additional harmonic potential with constant force of 2.0 kcal/mol subjecting the system to the phase of equilibration. Equilibration of the system was performed in three steps: (1) gradual heating of the system using a Langevin thermostat with a collision frequency of 2.0 ps À1 heated from 0 K to 310 K for 1 ns in NVT ensemble;

(2) equilibration of pressure was done at 1.0 bar in NPT for 1 ns ensemble via Berendsen barostat with a 2.0 ps (Case et al., 2020) ; (3) equilibration was achieved at 310 K over 0.5 ns in conditions of constant volume and the complex assembly constraints were eliminated. Once the equilibration was accomplished, the simulations of MD were performed at temperature T ¼ 310 K and p ¼ 1 bar and for 150 ns in NPT ensemble. SHAKE algorithm was used to constrained bonds participating in hydrogen atoms (Ryckaert, 1977) for accomplishing the integration time-step of 2 ps. Calculations of the electrostatic Long-range interactions were made via Particle Mesh Ewald algorithm (Essmann et al., 1995) . Truncations at 8 Å of van der Waals and Coulomb interactions and the calculations of energy values of binding were made with Amber18 (Case et al., 2020) via the MM/GBSA method (Genheden & Ryde, 2015) . 750 snapshots were calculated produced from the final 150 ns of the MD trajectories, by maintaining the snapshots every 0.2 ns. Computations of the energies of the polar solvation in continuum solvent via Poisson-Boltzmann continuum-solvation model with ionic strength of 0.10. The non-polar terms were evaluated via solvent accessible surface areas. Application of the Nmode module in Amber 18 was performed for the calculation of the entropy for the binding free energy (Case et al., 2020) . MD trajectories analysis was carried by the CPPTRAJ module of AmberTools 18 (Case et al., 2020) .

Secondary metabolites of crude methanolic extract of A. bidwillii, A. cunninghamii, and A. heterophylla dereplication revealed the annotation of various compounds of different classes. The annotated compounds detect the presence of diverse phenolic acids, biflavonoids and lignans Table 1 and Figure 1 .

The concentration inhibiting 50% of the release of PGE2 (IC 50 ) of the extracts under investigation, A.bidwillii, A.cunninghamii and A. heterophylla to examine the extracts anti-inflammatory activity was determined from the curve of concentration-inhibition response through regression analysis. The reported values are the means of three experiments and represented in Table 2 and Figure 2.

The DPPH and its % inhibition of methanolic extracts showed IC 50 values 93.32 lg/mL, 53.7 lg/mL and 120.226 lg/mL for A. bidwillii, A. cunninghamii and A. heterophylla respectively. The reference drug used was ascorbic acid which showed 46 lg/ mL Among these results methanolic extract of A. cunninghamii exhibited a potent antioxidant activity than the other extracts followed by A. bidwillii and A. heterophylla extracts. The total results of inhibition percentage as represented in Table 3 , respective to IC 50 values.

Absorbance of controlÞ Â 10

Standard Deviation (SD) ± means data set of n ¼ 3 was determined through the subsequent equation, whilst the values of IC 50 were calculated from plotting inhibition percentage against concentration, via a non-linear regression algorithm. Presenting data as mean values ± of (n ¼ 3) standard deviation.

The plant extracts of A. bidwillii, A. cunninghamii, and A. heterophylla showed anti-inflammatory effects as shown by their ability to inhibit the release of PGE2. The possible inhibitory effect of compounds in the extracts for COX-1 was studied using docking. The docking study was done using 1EQG PDB file which has co-crystalized ibuprofen. To establish the docking procedure, ibuprofen was redocked in its active site and the correct pose was predicted with high accuracy (RMSD of 0.381) and the docking energy was found to be À7.8 kcal/mol ( Figure 3A ). This indicates the validity of docking procedure to predict correct pose. Ibuprofen (and most of the anti-inflammatory non-steroidal drugs, NSAIDS) forms 2 hydrogen bonds with R120 and Y355. The hydrophobic part of ibuprofen extends in a hydrophobic pocket formed around I345, V349, L359, I523. Selective COX-2 inhibitors, on the other hand, do not inhibit COX-1 to the same extent.

Although that is the case, a selective COX-2 inhibitor, celecoxib, is known to bind to one of two subunits of COX-1 and the interference of the action of low dose aspirin (Rimon et al., 2010) . Binding site of celecoxib is close but different from that of known NSAIDs such as celecoxib ( Figure 3B ). Celecoxib sulphamoyl moiety forms hydrogen bonds with S516, Q192 and L352 while hydrophobic interactions are seen with V116, L531 and V349. Redocking of celecoxib in COX-1 resulted in an energy of À5.0 kcal/mol. Among the tested compounds lariciresinol glucoside, octadecyl (E)-p-coumarate and p-coumaroylquinic acid showed best docking scores which is comparable or even slightly higher than that of ibuprofen (Table 4 ). All the three compounds can form hydrogen bonds with the same residues similar to ibuprofen including R120 and Y355 as well as multiple hydrophobic interactions. In addition, lariciresinol glucoside and p-coumaroylquinic acid both formed hydrogen bonds with Y385 and p-coumaroylquinic acid has an extra hydrogen bond with S530 which is in close proximity to the binding pocket of celecoxib reported earlier. Proposed binding modes of these three compounds with COX-1 are illustrated in Figure 4 .

A recent study has shown that the anti-inflammatory corticosteroid dexamethasone has positive effect in improving survival rate of COVID-19 hospitalized patients as demonstrated by randomized clinical trials (Abdolahi et al., 2020) . This suggests promising effect of anti-inflammatory drugs as adjuvant therapy in COVID-19 patients for challenging the storm of cytokine which is a state of release of excessive of a diverse of inflammatory cytokines. The molecular mechanism of the cytokine storm is still unexplored. We were encouraged to theoretically study possible effect of our identified compounds against multiple targets of SARS-CoV-2. The compounds identified were docked in 3 targets active sites, viral main protease (M pro , 6LU7), non-structural protein-16 which is a methyl transferase (nsp16, 6W4H) and RNA dependent RNA polymerase (nsp12, 7BV2). Table 4 shows the docking results of tested compounds against the 3 viral targets. In general, bioflavonoids have shown better scores compared to other compounds in all the three targets. Among the three targets, nsp16 was the best target recognized by the tested compounds as can be deduced from docking energies of each compound with the 3 targets (Table 4 ). Three biflavonoids were among the top compounds bound to nsp16 with energies better than S-adenosylmethionine (SAM) which is the co-crystalized ligand. These biflavonoids include kayaflavone, amentoflavone and Bilobetin glucoside. A representative example of the binding of kayaflavone and amentoflavone is shown in Figure 5 . Rings A and B of biflavonoids overlapped with the purine ring of the SAM and form similar hydrogen bonds with D6912 and C6913. Other hydrogen bonds are also formed with D6897. The other side of the bioflavonoids extends into the KDKE methyl transferase motif and forms hydrogen bonds with E7001 in case of kayaflavone and with K6968 for amentoflavone. These results suggest possible effect of these compounds against SARS-CoV-2 and worth further investigation.

So, the molecular docking data revealed that all selected compounds effectively block the target proteins and that is confirmed by the binding energy values lower compared with those predicted for the molecules used as a positive control in the calculations. However, it is well-known that estimating the binding affinity profiles using molecular docking studies involves a number of approximations, resulting in the Lack of trustworthiness on the capability of scoring functions to produce accurate values of binding free energy (Meng et al., 2011; Salmaso & Moro, 2018) . In this context, the parameters of thermodynamic depending on the MD simulations are much more credible than those predicted from the analysis of the static docking models. contrary of the computational approaches, MD simulation moves each atom individually through the field of the rested atoms showing the flexibility of both the ligand and protein which is more effective than other algorithms (Meng et al., 2011; Salmaso & Moro, 2018) .

Nevertheless, analyzed compounds to the target proteins resulted from the post-modeling analysis of the static ligand/ protein structures their data, is confirmed by the MD simulations the high binding affinity. Analysis of the dynamic ligand/ protein models shows that, in both cases of interest, within the MD simulations, these complexes are relatively stable, as proved by the averages of binding free energies and the corresponding standard deviations (Table 5) . With the Given MM/ GBSA method errors of 1-3 kcal/mol (Genheden & Ryde, 2015) , we can propose that the dynamic ligand/protein structures having the averages of binding free energy lower than those calculated for the reference compounds using the same computational protocols (Table 5 ). The exception is compound B' targeting the NSP16-10 protein of SARS-Cov-2 and exhibiting the binding affinity comparable with that obtained for the control molecule SAM (Table 5) .

Regarding the observed stability of the dynamic ligand/ protein structures is also an evidence for the data depending on the time of the root-mean square deviations (RMSD) for the atomic positions of the static and dynamic models of the analyzed compounds bounding to the target proteins (Figures 6 and 7) . Analysis of Figures 6 and 7 reveals that these complexes do not support significant structural rearrangements on the MD trajectories, and this is confirmed by the calculated RMSD averages of the dynamic models of molecules A, B, C and A, B, C in the complexes with COX-1 and NSP16-10, respectively. For the ligand/COX-1 complexes, the RMSD mean values and standard deviations, which are 2.48 ± 0.58 Å (compound A), 2.74 ± 0.35 Å (compound B), 2.64 ± 0.37 (compound C), are close to those of 2.34 ± 0.36 Å calculated for the COX-1 inhibitor IBP ( Figure 6 ). In the case of the ligand/NSP16-10 complexes, these averages are 2.94 ± 0.74 Å (compound A'), 2.90 ± 0.33 Å (compound B'), 2.55 ± 0.26 Å (compound C'), and 3.31 ± 0.37 Å (SAM).

Comparison of the MD structures between every subsequent dynamic model with the previous one results in the averages of RMSD which are also evidence of the relative stability of the ligand/protein complexes at MD simulations. For the compounds bound to COX-1, these averages equal to 0.77 ± 0.05 Å (compound A), 0.78 ± 0.05 Å (compound B) and 0.78 ± 0.05 Å (compound C) are very close to those calculated for the control molecule IBP (0.76 ± 0.05 Å). Compounds À38.93 ± 4.02 À23.60 ± 5.42 À15.33 ± 6.44 B'

À27.84 ± 7.01 À19.21 ± 3.29 À8.63 ± 7.10 C'

À44.14 ± 3.26 À24.69 ± 5.97 À19.45 ± 5.84 SAM À27.30 ± 4.69 À18.60 ± 7.19 À8.70 ± 8.51

The compounds are designated as follows: A-3-Lariciresinol-glucoside, B-14-Octadecyl (E)-p-coumarate, C-16-p-coumaroylquinic acid, IBP-Ibuprofen. A'-5-Kayaflavon, B'-8-Amentoflavone, C'-11-Bilobetin-glucoside, SAM-S-Adenosylmethionine.

Values of mean of enthalpic and entropic components of free energy are <DH > and < TDS>, respectively; standard deviations corresponding to these values are (DH) STD and (TDS) STD À0.61 ± 0.27 À0.56 ± 0.27 À1.51 ± 0.55 a Data for the COX-1 residues with the binding enthalpy À0.4 kcal/mol are presented. b The averages of the residue contributions to the binding enthalpy and corresponding standard deviations are given. c The COX-1 residues dominating the ligand/protein interaction are highlighted by bold. A', B', C' show the mean values of RMSD in the complexes with NSP16-10 which are respectively 0.84 ± 0.06 Å, 0.84 ± 0.06 Å, 0.83 ± 0.06 Å and very close to the SAM/NSP16-10 structure (0.85 ± 0.06 Å).

Enthalpy components decomposition of the binding energy to the contributions of individual COX-1 amino acids exposes the residues dominating the ligand/COX-1 interactions in all cases of interest ( Table 6) . Analysis of Table 6 shows that these COX-1 residues are . The data obtained suggest that there are a number of the basic anchoring residues of COX-1 that provide strong attachment of the analyzed compounds at the active site of this protein. This supposition is supported by the values of the data of root-mean-square fluctuations (RMSF) for residues of COX-1 designating the flexibility of each amino acid throughout the MD simulations ( Figure 8 , Table 7) . Analysis of Table 7 shows that, according to the values of RMSF, the key anchoring residues and Ala-527 of COX-1 their position are restrained on the MD trajectories, in agreement with the data on their contributions into the binding enthalpy ( Table 6) .

The data on decomposition of binding energy into the contributions of individual NSP16-10 amino acids are given in Table 8 , and the data on the values of RMSF for the individual residues of this protein are presented in Figure 9 and Table 9 .

Intermolecular calculations of hydrogen bonds appeared in the MD trajectories of the ligand/COX-1 complexes indicate (Table 10 ) that the concerned compounds form hydrogen bonds with Arg-120 COX-1 (compound A), Leu-117 COX-1 (compound B), Ala-527 COX-1 and Leu-531 COX-1 (compound C). Notably, as compound C, the control molecule IBP is also involved in hydrogen bonding with Leu-531 COX-1 (Table 10) . Analysis of the dynamic ligand/NSP16-10 models reveals hydrogen bonds between compound A' and Asn-6899 NSP16-10 , compound B' and Lys-6933 NSP16-10 , as well as compound C' and Lys-6930 NSP16-10 , Asp-6931 NSP16-10 , Met-6929 NSP16-10 , and Cys-6913 NSP16-10 . At the same time, the reference compound SAM forms hydrogen bonds with the NSP16-10 residues Thr-4292, Asn-4293, Val-4295 (Table 10) . Among these H-bonds, it should be specially noted hydrogen bonds CÁÁÁAla-527 COX-1 and C'ÁÁÁLys-6930 NSP16-10 demonstrating the relatively high percentage occupancies on the MD trajectories (Table 10) .

Thus, the data of molecular modeling indicate that the analyzed compounds show strong attachment to the target proteins. This follows from the calculation of low binding free energy values for the ligand/protein complexes in terms of the scoring functions of molecular docking and molecular dynamics. In both cases, the predicted values of this À0.71 ± 0.91 -a Data for the NSP16-10 residues with the binding enthalpy À0.4 kcal/mol are presented. b The averages of the residue contributions to the binding enthalpy and corresponding standard deviations are given. c The NSP16-10 residues dominating the ligand/protein interaction are highlighted by bold. 3.12 --

The NSP16-10 residues presenting the dominant contributors to the ligand/ protein interaction are highlighted by bold. thermodynamic parameter testify to the higher binding affinity of the identified compounds to their targets compared with the reference molecules IBP and SAM.

In light of the results of the current investigation, we propose that the extracts of A. bidwillii, A. cunninghamii and A. heterophylla cultivated in Egypt, could be a possible source of potent natural anti-inflammatory and antioxidant metabolites. Those activities are of extensive importance as medicated agents for inhibiting and regress the progression of inflammation correlated with oxidative stress-relevant to viral infection and trying to avoid the side effects accompanied by the administration of dexamethasone such as elevated blood pressure, oedema and hyperglycemia. Those results were deduced from the virtual study carried for the structure relationship of LC-MS annotated compounds and the receptor interactions in three different binding sites of proteolytic enzyme in the COVID-19 as all the tested exhibited significant binding stability as observed for the amentoflavone and bilobetin glucoside with significant low binding energy. Our data suggested possible effect of these compounds against SARS-CoV-2 and worth further investigation. Several pathological conditions for instance as processes of inflammation are accompanied with reactive oxygen species generation. Thus, inhibitory activity on COX-1 enzyme of these plant extracts probably due to their capacity to reduce oxidative stress. Finally, we believe that the potential anti-COVID effects of the reported compounds required and worth further investigations.

No potential conflict of interest was reported by the authors. Donors and acceptors of the hydrogen bonds relating to the ligands are shown first, followed by the corresponding functional groups of the protein amino acids. The protein residues and percentage occupancies of hydrogen bonds are indicated in square brackets. The data on the hydrogen bonds exhibiting percentage occupancies >10%are given.

Letter to the editor: Efficacy of different methods of combination regimen administrations including dexamethasone, intravenous immunoglobulin, and interferon-beta to treat critically ill COVID-19 patients: A structured summary of a study protocol for a randomized controlled trial

Activity guided fractionation of Arum italicum miller tubers and the LC/MS-MS profiles

Molecular docking reveals the potential of aliskiren, dipyridamole, mopidamol, rosuvastatin, rolitetracycline and metamizole to inhibit COVID-19 virus main protease

Phytochemical and ethno-pharmacological review of the genus Araucaria-review

LC coupled to ion-trap MS for the rapid screening and detection of polyphenol antioxidants from Helichrysum stoechas

COVID-19 and the role of chronic inflammation in patients with obesity

In vivo antiulcer activity, phytochemical exploration, and molecular modelling of the polyphenolic-rich fraction of Crepis sancta extract

Identification of lignans by liquid chromatography-electrospray ionization ion-trap mass spectrometry

Investigation of some chemical constituents, cytotoxicity and antioxidant activities of Beta Vulgaris var. altissima cultivated in Egypt

Essential oil and its nanoemulsion of Araucaria heterophylla resin: Chemical characterization, anti-inflammatory, and antipyretic activities

A smooth particle mesh Ewald method

LC/MS/MS characterization of phenolic constituents in dried plums

The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities

Anti-inflammatory, antioxidant and antimicrobial activity of the stem bark extract and fractions of Ficus exasperata Vahl.(Moraceae)

Enhancing immunity in viral infections, with special emphasis on COVID-19: A review

Structure-based drug design, virtual screening and high-throughput screening rapidly identify antiviral leads targeting COVID-19

Comparison of simple potential functions for simulating liquid water

Screening seeds of some Scottish plants for free radical scavenging activity

Bioactivity-guided isolation of anti-inflammatory triterpenoids from the sclerotia of Poria cocos using LPS-stimulated Raw264.7 cells

COVID-19 cytokine storm: The anger of inflammation

Potential cytoprotective activity of ozone therapy in SARS-CoV-2/COVID-19

Phytochemical investigation of Cycas circinalis and Cycas revoluta leaflets: Moderately active antibacterial biflavonoids

An overview of severe acute respiratory syndrome-coronavirus (SARS-CoV) 3CL protease inhibitors: Peptidomimetics and small molecule chemotherapy

Enzyme immunoassays of eicosanoids using acetylcholine esterase as label: An alternative to radioimmunoassay

Flavonoids from Ephedra sinica stapf

Coxibs interfere with the action of aspirin by binding tightly to one monomer of cyclooxygenase-1

Ciccottiþ, G. & Berendsen, HJC numerical integration of the Cartesian equations of motion of a system with constraints: Molecular dynamics of N-alkanes

Qualitative analysis of phenolic compounds in apple pomace using liquid chromatography coupled to mass spectrometry in tandem mode

COVID-19 revisiting inflammatory pathways of arthritis

The PyMOL molecular graphics system

HPLC-PDA-MS/MS characterization of bioactive secondary metabolites from Turraea fischeri bark extract and its antioxidant and hepatoprotective activities in vivo

A simple screening method for antioxidants and isolation of several antioxidants produced by marine bacteria from fish and shellfish

AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading

Development and testing of a general amber force field

Analysis of the total biflavonoids extract from Selaginella doederleinii by HPLC-QTOF-MS and its in vitro and in vivo anticancer effects

COVID-19: The inflammation link and the role of nutrition in potential mitigation

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

Molecular docking: A powerful approach for structure-based drug discovery

Bridging molecular docking to molecular dynamics in exploring ligand-protein recognition process: An overview

The authors would like to extend their sincere appreciation to the central laboratory at Jouf University for support this study. M.M. Al-Sanea also extends his appreciation to the Korea Institute of Science and Technology (KIST) for funding this work through the Grant name "2021 KIST School Partnership Project".

The authors extend their appreciation to the Deputyship for Research and Innovation, Ministry of Education in Saudi Arabia for funding this work through the project number "375213500"