key: cord-253431-fjds5cdr
authors: Erukainure, Ochuko L.; Atolani, Olubunmi; Muhammad, Aliyu; Katsayal, Sanusi B.; Ebhuoma, Osadolor O.; Ibeji, Collins U.; Mesaik, M. Ahmed
title: Targeting the initiation and termination codons of SARS-CoV-2 spike protein as possible therapy against COVID-19: the role of novel harpagide 5-O-β-D-glucopyranoside from Clerodendrum volubile P Beauv. (Labiatae)
date: 2020-11-03
journal: Journal of biomolecular structure & dynamics
DOI: 10.1080/07391102.2020.1840439
sha: 
doc_id: 253431
cord_uid: fjds5cdr

The global spread of the coronavirus infections disease − 2019 (COVID-19) and the search for new drugs from natural products particularly from plants are receiving much attention recently. In this study, the therapeutic potential of a new iridoid glycoside isolated from the leaves of Clerodendrum volubile against COVID-19 was investigated. Harpagide 5-O-β-D-glucopyranoside (HG) was isolated, characterised and investigated for its druglikeness, optimized geometry, and pharmacokinetics properties. Its immunomodulatory was determined by chemiluminescence assay using polymorphonuclear neutrophils (PMNs) in addition to T-cell proliferation assay. In silico analysis was used in determining its molecular interaction with severe acute respiratory syndrome coronavirus-2 (SARS-COV-2). HG displayed potent druglikeness properties, with no inhibitory effect on cytochrome P(450) (1A2, 2C19, 2C9, 2D6 and 3A4) and a predicted LD(50) of 2000 mg/kg. Its (1)H-NMR chemical shifts showed a little deviation of 0.01 and 0.11 ppm for H-4 and H-9, respectively. HG significantly suppressed oxidative bursts in PMNs, while concomitantly inhibiting T-cell proliferation. It also displayed a very strong binding affinity with the translation initiation and termination sequence sites of spike (S) protein mRNA of SARS-COV-2, its gene product, and host ACE2 receptor. These results suggest the immunomodulatory properties and anti-SARS-COV-2 potentials of HG which can be explored in the treatment and management of COVID-19. Communicated by Ramaswamy H. Sarma

The recent pandemic of coronavirus disease 2019 has led to an unprecedented rate of mortality, with Italy, China, Spain and USA being the most hit (WHO., 2020) . It is caused by the novel b-coronavirus (enveloped non-segmented positive-sense RNA virus), SARS-CoV-2 which has a similar sequence identity to that of bat/human severe acute respiratory syndrome coronavirus-related coronavirus (SARSr-CoV) and bat SARS-like coronavirus (SLCoV) Schoeman & Fielding, 2019; Tian et al., 2020) . Since its origin in Wuhan, China as a cluster of pneumonia cases in December 2019 which has spread globally (da Silva et al., 2020) , there have been search for effective antiviral remedies against the diseases. Lack of effective treatment coupled with no valid vaccines increases the number of infections and mortalities on a daily basis (Mehta et al., 2020; Onder et al., 2020) . The common symptoms of COVID-19 include sore throat, cough, fatigue, slight dyspnoea, headache, conjunctivitis, diarrhea and fever (Pascarella et al., 2020) . The major cause of COVID-19 mortality is respiratory failure as present treatments are based on symptomatic management and respiratory sustenance (Pascarella et al., 2020) . Underlying ailments particularly diabetes mellitus and obesity have been implicated as high-risk factors of COVID-19 mortalities (Holman et al., 2020) . Boosting the immune system has also been reported as a preventive strategy against COVID-19, as cytokine release syndrome (CRS) characterized by exacerbated levels of macrophage inflammatory protein 1-a, tumour necrosis factor-a (TNF-a), interleukins (IL)-2 and À7, monocyte chemoattractant protein 1, and granulocytecolony stimulating factor have been implicated in patients with severe conditions (Mehta et al., 2020; Shi et al., 2020) . In vitro studies have reported the inhibitory effect of remdesivir and chloroquine on SARS-CoV-2 activity , and have shown some therapeutic effects in clinical trials . However, studies have reported deaths in COVID-19 patients treated with chloroquine and its derivatives owing to severe arrhythmias . It has also been shown to suppress the immune system as well as inhibit autophagy and xenophagy . Interestingly, the dynamics of the said virus largely depends on the molecular interactions between spike protein and human angiotensin converting enzyme 2 receptor (ACE2). Thus by implication, targeting these biomarkers have been reported as viable strategies for treating the disease (Gurwitz, 2020; Jiang et al., 2012) .

The antiviral and immunomodulatory activities of medicinal plants have been reported Vijayan et al., 2004) . This has been attributed to their ability to modulate inflammatory cytokines and phagocytic activities of monocytes, neutrophils and macrophages Harun et al., 2015) , with their phytochemical constituents playing an influential role. Clerodendrum volubile is among medicinal plants reported for their immunomodulatory activities (Erukainure, Hafizur, et al., 2017; Erukainure, Zaruwa, et al., 2017; Erukainure et al., 2016) .

Clerodendrum volubile is a leafy vegetable mostly utilized as food ingredients in Southern Nigeria (Erukainure, Sanni, et al., 2018; Erukainure et al., 2011) . Regarded as magic leaves, C. volubile is used in folk medicine for treating ulcers, arthritis, rheumatism, dropsy, and diabetes (Burkill, 1985) . Its antioxidant activities have been reported in diabetic rats, and cancerous cells which correlates with their antidiabetic and anticancer activities (Erukainure, Ashraf, et al., 2018; Erukainure, Hafizur, et al., 2017; . The flowers as well as protocatechuic acid and pectolinarigenin isolated from the leaves have been reported for their ability to modulate phagocytic oxidative burst in isolated neutrophils and macrophages which demonstrates their immunomodulatory activities (Erukainure, Hafizur, et al., 2017; Erukainure et al., 2016) . The dichloromethane (DCM) fraction of the leaves have been reported for its cytotoxic effect on human embryonic kidney (HEK293) cells (Erukainure, Narainpersad, et al., 2018) . The leaves are rich in phenolics contents which has been correlated with its reported biological activities (Erukainure, Narainpersad, et al., 2018; Oboh et al., 2017) .

To the best of our knowledge, this study reports for the isolation and characterization of a new iridoid glycoside, harpagide 5-O-b-D-glucopyranoside from the leaves of C. volubile. Based on previous reports on the immunomodulatory properties of the leaf extract, it is hypothesized that the isolated compound may contribute to the modulation of the immune system in favor of suppressing CRS. Thus, this study was aimed at investigating the antioxidative burst potentials of the compound on polymorpho-nuclear neutrophils (PMNs), as well as its antiproliferative effect on T-cells. The ability of the compound to suppress the activity of SARS-CoV-2 was investigated using combined in silico tools. This study also reports global epidemiology study of COVID-19 from January -August 2020.

The 1 H-and 13 C-NMR spectra were recorded at temperature 300 K on Bruker Avance 500 MHz spectrometer in MeOD. d (ppm) values were used for reporting chemical shifts. The 1 H-NMR spectrum was recorded for 128 scan using pulse sequence with an acquisition time 1.58 sec, relaxation delay of 1.5 sec, spectra width of 10,330 Hz and Intrinsic digital resolution of 0.3152 Hz with tetramethylsilane, TMS as internal standard. Hitachi UV-3200 spectrophotometer was used in obtaining UV spectra. A pre-coated silica gel 60 F 254 plates (E. Merck, 0.25 mm) was used for thin layer chromatography (TLC), which was detected under UV light (254 nm) after spraying with ceric sulphate reagent. JEOL JMS-HX-110 mass spectrometer was used for recording EI-MS.

Computational software used for this work includes ArcGIS 10.6 software, UCSF Chimera (v. 1.0.1), Discovery Studio 2017 R2 Client (v17.2.0.16349), Python prescription virtual screening software (PyRx, 8.0), ChemDraw Ultra v12.0.2, PatchDock server and other online computational tools available.

We obtained freely available daily data of confirmed cases and deaths for Coronavirus Disease (COVID-19) pandemic across the world from the COVID-19 -Statistics and Research website from January 2020 to March 2020 (Roser et al., 2020; ECDC, 2020) . The COVID-19 -Statistics and Research team retrieves daily global and national statistics of COVID-19 from two credible sources, the World Health Organization (WHO) and, the European Centre for Disease Prevention and Control (ECDC). We aggregated the daily confirmed data (cases and deaths) for each country into monthly data for January to March 2020, and subsequently performed the spatiotemporal distribution analysis of COVID-19 confirmed cases and deaths using the ArcGIS 10.6 software (ESRI, Redland, CA, USA).

Clerodendrum volubile leaves were obtained from local farmers at Ifon, Ondo State, Nigeria. The leaves were assigned a voucher number, UBH C284 after identification and authentication by Dr. Henry Akinbosun at the Herbarium, Department of Botany, University of Benin, Benin City, Nigeria. The leaves were air-dried to minimize loss of volatiles and light-sensitive compounds (Thamkaew et al., 2020) before blending to fine powder and stored in a zip-lock bag at room temperature until further analysis (2 days).

The blended sample was subjected to methanol (MeOH) extraction at room temperature. The solvent choice was based on the polarity of MeOH and reports as an optimal solvent for high content of phytochemical constituents (Truong et al., 2019) . The extract was concentrated in vacuo using a rotary evaporator. About 100 g of the extract was dissolved in MeOH/Distilled water (1:3), and thereafter subjected to gradient polarity based chemical fractionation using, n-hexane (Hex), n-dichloromethane (DCM), ethyl acetate (EtOAc), and butanol (BuOH). All the fractions were collected, concentrated in vacuo and stored in glass vials at À20˚C. The EtOAc fraction was parked into a silica gelloaded column chromatography for further fractionation. The column was subjected to solvent elution starting with 100% hexane, and thereafter mixtures of Hex and EtOAc (9:1) in increasing order of polarity up to 100% EtOAc. Each elute was collected in glass vials and its purity was confirmed by subjection to thin-layer chromatography on pre-coated silica gel 60 F 254 sheets. They were sprayed with ceric sulfate reagent, dried under hot-air (60 C), and viewed under UV light (254 nm). Elutes from Hex and EtOAc (8:2) showed a single spot. They were combined and concentrated in a fume compound, before subjection to NMR analysis for structural elucidation and identification. The compound was identified as harpagide 5-O-b-D-glucopyranoside.

The isolated compound (harpagide 5-O-b-D-glucopyranoside) was fully optimized using density functional theory (DFT) with M06-2X functional with 6-31 þ G(d,p) basis. The frequency calculation was computed confirming no negative frequencies. The scaling factor of 0.964 was applied to evaluate vibrational frequencies (Johnson III, 1999) . The 1 H magnetic shielding constants, by chemical shifts, obtained on a d-scale relative to the TMS, used as reference was calculated on complex 1 and 2 using the Gauge-Independent Atomic Orbital (GIAO) method developed by Wolinski et al. (1990) . Polarizable Continuum Model (PCM) through a single point (B3LYP/6-311 þ G(2d,p)) calculation was used for determining the 1 H NMR chemical shift (DMSO solvent: dielectric constant, e ¼ 46.826) (Cances et al., 1997; De Souza et al., 2017) .

To determine the druglikeness properties of the selected compounds, a supercomputing facility for Bioinformatics and Computational Biology (SCFBio platform; scfbio-iitd.res.in) was used. 3 D structures of the compounds were submitted through the platform and the prediction carried out based on Lipinski's rule. Lipinski's rule considered five basic properties including molecular weight, hydrogen bond acceptor and donor, lipophilicity index (LogP) and molar refractivity of a compound (Jayaram et al., 2012; Lipinski, 2004) .

The pharmacokinetics and oral lethal dose of the isolated compound was determined in silico via qualitative structure activity relationship (QSAR) vis-a-vis virtual molecular structure activity relationship studies (SARs) using the SwissADME 

Blood sample was collected from a consented apparently healthy volunteer, from which neutrophils were immediately isolated as previously described (Alyiu et al., 2014) . The isolated neutrophils were subjected to luminol-amplified chemiluminescence assay (Helfand et al., 1982) . Briefly, the isolated neutrophils (1 Â 10 6 /mL) were suspended in modified Hank's solution (MHS) and incubated with different concentrations (0.5, 5.00 and 50.00 mg/mL) of the isolated compound (harpagide 5-O-b-D-glucopyranoside) for 30 min. Zymosan (1 mg/ mL) was thereafter added, followed by 25 mL of luminol (10 À5 M). A luminometer (Lab system Luminoskan RS, Helsinki, Finland) was used in recording the total chemiluminescence (CL). Control (þC) consisted of MHS containing cells only.

Blood obtained from a consented apparently healthy volunteer was subjected to Ficoll-Hypaque gradient centrifugation and proliferation to extract T-cells as described previously (Nielsen et al., 2003) . The cells were cultured in a 96-well round bottom tissue culture plate at a concentration of 5 Â 10 5 /mL. They were thereafter stimulated with 1.25 mg/mL of phytohaemagglutinin and incubated with the isolated compound (0.5, 5.00 and 50.00 mg/mL) at 37 C in a 5% CO 2 incubator for 72 h. After incubation, the reaction was pulsed with titrated thymidine 0.5 mCi/well, and further incubated for 18 h. The cells were harvested with a cell harvester (SKATRON A.S. Flow Lab., Norway) into a glass fibre-filter (Cambridge Technology USA), and measured by counting with a liquid scintillation counter. After 120 s, the counts per minute (cpm) results were recorded . All studies were carried out under the ethical guidelines of the International Center for Chemical and Biological Sciences, University of Karachi, Karachi, Pakistan.

SARS-CoV-2 spike protein mRNA sequence SARS-CoV-2 spike protein mRNA sequence was obtained from the PubMed database from the coronavirus whole genome (NC_045512.2). About 42 nucleotide sequence from both initiation and termination codons (Table 1) were converted into 3 D single mRNA strand and prepared for docking.

The binding potentials of Harpagide 5-O-b-D-glucopyranoside, Azithromycin, Chloroquine, and Remdesivir was studied by means of molecular docking. The 3 D crystal structures of 2019-nCoV spike protein receptor-binding domain bound with ACE2 angiotensin converting enzyme 2 receptor (ACE2) were downloaded from the protein data bank (PDB) with PDB ID; 6M0J (Lan et al., 2020) with resolution 2.45 Å. The suitable grid box was determined using AutoDock tools 1.5.4 (Sanner, 1999) . The ligands (CID_2719, 121304016 and DB00207) were obtained from PUBMED, while the isolated compound (harpagide 5-O-b-D-glucopyranoside) was drawn using ChemDraw Ultra v12.0.2 and the file converted to pdb format. Initial optimization was done using Gaussian 09 (Frisch et al., 2009) , to obtain minimized conformations. The determined dimension was X ¼ 26 Y ¼ 26 Z ¼ 26 with 1.00 Å as the grid spacing. The optimum binding site for the ligand (Erukainure et al., 2020; Yang et al., 2013) was determined using Lamarckian genetic algorithm method. Gasteiger charges were computed using the AutoDock Tools graphical user interface supplied by MGL Tools (Morris et al., 2009 ).

To address biological variability, wet experiments were repeated at least thrice. Where appropriate, results were presented as mean ± SD and subjected to one-way analysis of variance (ANOVA). Tukey's HSD-multiple range post-hoc test was used in obtaining significant differences between means at p < 0.05. IBM Statistical Package for the Social Sciences (SPSS) for Windows, version 23.0 (IBM Corp., Armonk, NY, USA) was used for the analyses.

The spatiotemporal distribution of confirmed monthly COVID-19 cases and deaths from January to March 2020 is shown in Figure 1 . The maps suggest the impacts of the disease (cases and deaths) were initially felt in China in January 2020 before it spread to other countries where noticeable spikes in confirmed COVID-19 cases and deaths were reported in subsequent months. This suggests the epicentre of COVID-19 shifted from China to Europe (affecting mainly Italy, Spain, Germany, France, UK, Switzerland, Belgium, Netherlands, Turkey and Austria) at the end of March 2020. Other countries that reported high confirmed cases and deaths of COVID-19 in March 2020 were the USA, Iran, Canada, South Korea, Brazil and Chile. However, African countries reported fewer cases and deaths of COVID-19 compared to countries in other continents from January to March 2020. While low confirmed cases and deaths of COVID-19 are continuously reported across Africa, US and Europe (especially Spain, France and UK) are still characterized by high confirmed COVID-19 cases and death. Other countries currently experiencing a surge in COVID-19 cases and deaths are Russia, India, Indonesia, Philippines, Iran, Iraq and most South American Countries (especially Brazil, Columbia, Argentina, Peru, Mexico, Chile). This suggests, South America is currently the new epicentre of COVID-19 (Roser et al., 2020; ECDC, 2020).

As shown in Figure 2 , harpagide 5-O-b-D-glucopyranoside was isolated for the first time from the leaves of C. volubile. Its structure was elucidated using the combination of NMR, infrared and mass spectra data:

Harpagide 5 Table 2 . Major hydrogen atoms are designated in the calculated 1 H-NMR spectrum.

Drug likeness properties of the selected compounds were determined using criteria put forward by Lipinski. The result present in Table 3 and presents shown shows the parameter assessed, the scores and the reference value for each parameter, respectively. Harpagide 5-O-b-D-glucopyranoside showed strong compliance with respect to number of hydrogen bond donors, LogP and molar refractivity, except for hydrogen bond acceptors.

In silico prediction analysis predicted harpagide 5-O-b-D-glucopyranoside to be absorbed poorly in the gastrointestinal tract and not permeable across the blood brain barrier (BBB) as shown in Table 4 . Harpagide 5-O-b-D-glucopyranoside was also predicted as a P-glycoprotein inhibitor. The compound was predicted not to be an inhibitor of cytochromes 

As depicted in Figure 4A , incubation of harpagide 5-O-b-Dglucopyranoside with neutrophils significantly (p < 0.05) inhibited ROS production stimulated by zymosan in PMNs with an IC 50 value of 0.33 þ 0.01 mg/mL. The activity was dose dependent with increasing concentration. The compound also significantly (p < 0.05) suppressed the proliferation of T-cell proliferation, with a 26% inhibitory activity compared to the control as shown in Figure 4B .

The binding potentials of Harpagide 5-O-b-D-glucopyranoside and standard FDA approved drugs (Azithromycin, Chloroquine, and Remdesivir) were studied by means of molecular docking. Ligand-protein interactions between viral protein (SARS-CoV-2 spike protein), the host receptor target (ACE2) and Harpagide 5-O-b-D-glucopyranoside are presented in Figures 5-7 . Harpagide 5-O-b-D-glucopyranoside displayed a good binding in complex with the host receptor target, initiation and termination sequence of the viral spike protein messenger RNA compared to all studied standard drugs with binding affinities of À7.5, À6.4 and 6.3 kcal mol À1 respectively (Table 5) . Interestingly, it showed comparable binding affinity with Remdesivir, when in complex with the SARS-CoV-2 spike protein with a slight difference of 0.2 kcal mol À1 but better binding compared to Azithromycin and Chloroquine (Figures S2 and S3; and Table 5 ).

As COVID-19 pandemic feasts, scientists and policy makers are currently making tireless efforts to reduce transmission components of medicinal plants against COVID-19 is receiving a little or no attention at this time. This is without taking cognizance of the facts that even the so-called drugs (Chloroquine, Remdesivir, Azithromycin etc.) that are being repurposed against COVID-19 had traceable origin to medicinal plants. Boosting the host immune system has also been reported to be beneficial in ameliorating the effect of the virus and the severe progression of the disease (Rothan & Byrareddy, 2020; Shi et al., 2020) . In the present study, we investigated the epidemiology of COVID-19 and the potentials of harpagide 5-O-b-D-glucopyranoside, a new iridoid glycoside isolated from C. volubile leaves as a potential candidate for the inhibition of the proliferation of SARS-CoV-2. The notable initial incidence and spatial distribution of COVID-19 cases and deaths across China (Figures 1A and 1B) can be attributed to its novelty (Guan et al., 2020; Sohrabi et al., 2020) . This resulted in the disease being transmitted through hospitals and infected medical staff unknowingly. In addition, vulnerable cases, 60 years and above with underlying medical conditions such as diabetes, cardiovascular disease, hypertension, asthma and other respiratory diseases contributed to the majority of the cases and deaths in China (Zhang, 2020) . However, the downward trend of the disease in China (Figures 1C -1F) can be related to her timely amendment of governmental policies which included enforcing restricted personnel movements in the epicentre of the disease (Wuhan city) and other areas of China Kraemer et al., 2020; Zhu et al., 2020) , sharing vital prevention and control information (including wearing of masks, adequate washing and sanitizing of hands, symptoms and means of seek medical assistance), enhanced public health resources, providing improved and functional infectious disease surveillance, early detection systems and multi-sectoral rapid response systems to assist curtail the disease (Zhang, 2020) . China's overall initiative and strategies in response to COVID-19 can be attributed to the lessons and experiences of the SARS epidemic that occurred in 2003 (Zhang, 2020) . Contrarily, limited medical resources and the delayed response in implementing policies restricting movements of people may have resulted in the surge of cases and deaths across Europe and other countries with high confirmed cases and deaths ( Figures 1E and 1F) . In other words, countries reacted slowly in implementing measures and policies to reduce the chance of contact between infected and vulnerable people. While the low recorded COVID-19 cases and related deaths across Africa can be linked with few vulnerable populations because Africa's life expectancy at birth is 61 years (World Health Organisation, 2018). Thus, a global health emergency that involves continuous collaborative surveillance, information sharing, researching and implementing evidence-based public health and clinical practices and formulating effective control strategies cannot be overemphasised. Furthermore, in-depth studies aimed at assessing the vulnerability of young people with underlying medical conditions and amongst different blood groups is important. These measures would be instrumental in the implementation of targeted intervention policies and resources.

The 1 H NMR spectral of the isolated compound was compared to that of harpagide and allobentonicoside, compounds of close structural similarity as obtained in literatures (Jensen et al., 1989; Manguro et al., 2011; Venditti et al., 2017) . While harpagide has only one glucose unit attached at C-1 to the iridoid aglycone unit, the allobentonicoside has two glucose units attached at C-1 and C-5 as in harpagide 5- O-b-D-glucopyranoside. However, allobentonicoside has the hydroxyl converted to a carbonyl at C-6 with an unsaturation at C-7. The isolated compound has very similar 1 H and 13 C NMR data compared to the allobentonicoside ( Figure S1 ). However, the absence of carbonyl stretching vibration in the infrared spectrum confirms that the compound is not allobentonicoside. The molecular ion peak at m/z 549 obtained from FAB-MS for the isolated compound which corresponds to C 21 H 34 O 15 [M þ Na] further strengthens the claim of two attached glucose units as in allobentonicoside. The de-glucosidated aglycone was presented at m/z 387 calculated for C 15 H 26 O 10 in the FAB-MS spectrum. The isolated compound was thus presented as harpagide 5-O-b-D-glucopyranoside as depicted in Figure 2 . The DFT calculated result for the 1 H-NMR and the calculated vibrational frequencies obtained showed good agreement with experimental ( Figure 3) . The 1 H-NMR chemical shifts, with a small deviation of 0.01 and 0.11 ppm for H-4 and H-9 (Table 2) , respectively. Others displayed slightly increased deviations from the experimental values, though these deviations are in accordance with previously calculated DFT 1 H-NMR calculations (De Souza et al., 2017; Olanrewaju et al., 2020) .

Druglikeness properties of compounds screen-out promising ligands that are likely to be used as drugs based on some identified properties (Dong et al., 2018) . The properties comprise those considered crucial for both kinetic and dynamic interactions inside the biological system. Essentially, harpagide 5-O-b-D-glucopyranoside displayed a potential druggable character, with the exception of a very few properties (Table 3) . Molecular weight indicates the mass as well as the size of a ligand and plays a greater role in determining the solubility, absorption and distribution of a given drug candidate (Doak et al., 2014) . Basically, the molecular weight of the new ligands reported in this work appears to fit within the acceptable range of this parameter. Hydrogen bond acceptors and donors in drug structures play an important role in water solubility, membrane transport, distribution, and drug-receptor interactions (Cheng et al., 2007; Kumar et al., 2010) . Five hydrogen bond donors and ten hydrogen bond acceptors are required for a ligand to possess a satisfactory bioavailability and the optimum number of drug-target interactions via hydrogen bonds (Dong et al., 2018) . LogP indicates a relative concentration of the solute in octanol and water (octanol-water partition coefficient). It determines the degree of hydrophobicity of a ligand and plays a role in predicting drug absorption via the phospholipid bilayer of the cell membrane (Kumar et al., 2010) . Generally, a more polar compound typically encounters difficulties when passing through phospholipid bilayer of the cell membrane (Yusof & Segall, 2013) . Harpagide 5-O-b-D-glucopyranoside demonstrated approvingly a very good partition coefficient (Table 3 ). This suggests that it can be efficiently transported across cell membranes, as passive diffusion through a cell membrane is a key element of transcellular absorption (Doak et al., 2014; Yusof et al., 2014) .

The isolated compound was subjected to in silico pharmacokinetics and lethal toxicity analysis in order to predict its possible drug-drug interactions and toxicity when swallowed ( Table 4 ). The low GI absorption value insinuates harpagide 5-O-b-D-glucopyranoside may not be easily absorbed by the intestine. The small intestine is known for its important role in drug metabolism as the large surfaces areas of villi and microvilli determine the rate of drug absorption via the oral route (Pang, 2003) . The low absorption of the compound may be attributed to its inability to permeate the phospholipid bilayers of the intestinal membrane owing to the lipophobic properties of most glycosides (Chao et al., 2006) . This may also be responsible for the predicted impermeability of the compound across the BBB. The prediction of the compound as a substrate of P-gp, suggests it is translocated from the intracellular to the extracellular compartment via transmembrane efflux mediated by P-gp (Kim, 2002) . This may be responsible for the low GI absorption of the compound as P-gp has been implicated in reduced bioavailability of its substrate (Constantinides & Wasan, 2007) . The predicted inability of the compound to inhibit CYPs 1A2, 2C19, 2C9, 2D6 and 3A4 insinuates that there will be little or no drugdrug interactions when co-administered with other drugs which are metabolized by these enzymes (Obach et al., 2006) . The predicted LD 50 value of 2000 mg/kg and toxicity class of 4 suggest the compound is safe for administration .

Phagocytic respiratory oxidative burst has been implicated in innate immunity and inflammatory processes which are the major pathogenesis of COVID-19. Boosting the immune system particularly during the incubation and non-severe stages of the disease, has been recognized as a strategy to eliminate the virus and to suppress the progression to severe stages . Phagocytic respiratory oxidative burst involves the generation of ROS from increased nonmitochondrial oxygen consumption arising from dehydrogenation of glucose via the hexose monophosphate shunt (Allen, 1994; Erukainure, Hafizur, et al., 2017) . Continuous phagocytic-induced ROS production will lead to an imbalance in the innate antioxidant system, thereby resulting in oxidative stress and provocation of inflammatory cytokines (Ciz et al., 2012) . Exacerbated inflammatory cytokine production often regarded as cytokine release syndrome (CRS) has been implicated in the onset and progression of COVID-19 stages (Mehta et al., 2020; Shi et al., 2020) . Thus, arresting phagocytic respiratory oxidative bursts may be a therapeutic strategy in treating and managing the disease and its progression. In the present study, harpagide 5-O-b-D-glucopyranoside suppressed the production of ROS in PMNs as shown in Figure 4A . Thus, indicating a suppressive activity on phagocytic respiratory oxidative burst and may be a good drug candidate in the treatment and management of COVID-19 and its complications. This activity correlates with our previous reports on the ability of extracts and isolated compounds from C. voluble leaves and flower to modulate phagocytic oxidative burst in whole blood, PMNs and macrophages (Erukainure, Hafizur, et al., 2017; Erukainure et al., 2016) .

The suppressive activity on respiratory oxidative burst by harpagide 5-O-b-D-glucopyranoside is further corroborated by its ability to suppress T-cell proliferation ( Figure 4B ). Proliferation of T-cell has been implicated in the activation of CRS, and studies have suggested that targeting T-cell proliferation may suppress the severe progression of COVID-19 .

At molecular level, the viral envelope spike (S) protein of SARS-CoV-2 and angiotensin converting enzyme 2 receptor within the host are central to COVID-19 pathogenesis and response to therapeutic interventions among other biological factors . Basically, it has been reported that envelope spike (S) protein is important for coronavirus pathogenesis as it mediates receptor binding and membrane fusion implicated in host tropism and transmission capability . Mostly, the S protein is functionally divided into the S1 domain, responsible for receptor binding, and S2 domain, responsible for cell membrane fusion. Structure analysis suggests that the receptor-binding domain is composed of a core and an external subdomain. Angiotensin converting enzyme 2 (ACE2) has been reported to act as a cell receptor for the S protein (Kuster et al., 2020) . Therefore, targeting S protein and ACE2 at transcriptomic and proteomic levels could be a novel strategy for treatment and management of COVID-19. The docking simulations revealed that harpagide 5-O-b-D-glucopyranoside binds better with the evidence of hydrogen bonding interactions with amino acid residues around the binding pocket of the host receptor target as well as the nitrogenous base and ribose sugar of the mRNA sequence of initiation and termination sequence of the viral spike protein m RNA (Figures 5-7) . Remdesivir in complex with SARS-CoV-2 spike protein displayed stronger molecular interactions when compared with CQ and AM ( Figure S3 ). However, it is comparable with harpagide 5-O-b-D-glucopyranoside. This can be as a result of the four hydrogen bonding interactions with four amino acid residues in the bind pocket of the SARS-CoV-2 spike protein compared to three of hydrogen bonding interactions found in harpagide 5-O-b-D-glucopyranoside ( Figure S2 ). Remdesivir have been reported the recently to inhibit novel SARS-COV-2 in vitro, which invariably corroborated our findings (Table 5) . The good binding potential of Harpagide 5-O-b-D-glucopyranoside can prevent or impair the recruitment of initiating and termination factors around mRNA sequence, thus affecting the ultimate translation of S protein mRNA. Equally, the strong and robust interactions were visible between the compound and S protein as well as ACE2 receptor (Figures 6  and 7) . Thus, qualifying it as a potential inhibitor of S protein mRNA translation vis-a-vis barriers that affect S protein and ACE2 receptor interactions implicated in COVID-19 pathogenesis (Kuster et al., 2020) . However, further studies are needed to evaluate the actual efficacy of HG against SARS-COV-2 at in vitro and in vivo levels while targeting spike protein and its receptor at transcriptomic and proteomic levels.

Taken together, these results suggest the immunomodulatory properties of harpagide 5-O-b-D-glucopyranoside isolated from C. volubile leaves which may be explored in the treatment and management of COVID-19. This is further corroborated by its strong molecular interaction and binding with the translation initiation and termination sequence sites of S protein mRNA and S protein of SARS-COV-2 and its host receptor protein (ACE2 receptor). In addition, there is a need for implementation of targeted intervention policies and resources which should include a global health emergency involving continuous collaborative surveillance, information sharing, research and implementation of evidence-based public health and clinical practices.

Role of oxygen in phagocyte microbicidal action

Antioxidant, mitogenic and immunomodulatory potentials of acacia honey

ProTox-II: A webserver for the prediction of toxicity of chemicals

The useful plants of west tropical Africa

A new integral equation formalism for the polarizable continuum model: Theoretical background and applications to isotropic and anisotropic dielectrics

Bioavailability, metabolism, and pharmacokinetics of glycosides in Chinese herbs

Structure-based maximal affinity model predicts small-molecule druggability

Flavonoids inhibit the respiratory burst of neutrophils in mammals

Lipid formulation strategies for enhancing intestinal transport and absorption of P-glycoprotein (P-gp) substrate drugs: In vitro/in vivo case studies

Forecasting Brazilian and American COVID-19 cases based on artificial intelligence coupled with climatic exogenous variables

SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules

Structural analysis of flavonoids in solution through DFT 1H NMR chemical shift calculations: Epigallocatechin, Kaempferol and Quercetin

Oral druggable space beyond the rule of 5: Insights from drugs and clinical candidates

ADMETlab: A platform for systematic ADMET evaluation based on a comprehensively collected ADMET database

Fatty acids rich extract from clerodendrum volubile suppresses cell migration; abates oxidative stress; and regulates cell cycle progression in glioblastoma multiforme (U87 MG) cells

Anti-diabetic effect of the ethyl acetate fraction of Clerodendrum volubile: Protocatechuic acid suppresses phagocytic oxidative burst and modulates inflammatory cytokines

Suppressive effects of clerodendrum volubile P Beauv. [Labiatae] methanolic extract and its fractions on type 2 diabetes and its complications

Dacryodes edulis (G. Don) HJ Lam modulates glucose metabolism, cholinergic activities and Nrf2 expression, while suppressing oxidative stress and dyslipidemia in diabetic rats

Pectolinarigenin from the leaves of Clerodendrum volubile shows potent immunomodulatory activity by inhibiting T À cell proliferation and modulating respiratory oxidative burst in phagocytes

Flowers of Clerodendrum volubile exacerbate immunomodulation by suppressing phagocytic oxidative burst and modulation of COX-2 activity

Clerodendrum volubile inhibits key enzymes linked to type 2 diabetes but induces cytotoxicity in human embryonic kidney (HEK293) cells via exacerbated oxidative stress and proinflammation

Nutritional qualities and phytochemical constituents of Clerodendrum volubile, a tropical non-conventional vegetable

Clerodendrum volubile: Phenolics and applications to health

Suppression of phagocytic oxidative burst, cytotoxic effect, and computational prediction of oral toxicity of dietary fatty acids of Clerodendrum volubile stem

COVID-19 data

Epidemiological characteristics of 1 052 cases of COVID-19 in epidemic clusters

Clinical characteristics of coronavirus disease 2019 in China

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

Angiotensin receptor blockers as tentative SARS-CoV-2 therapeutics

Immunomodulatory effects of selected Malaysian plants on the CD18/11a expression and phagocytosis activities of leukocytes

Coronavirus disease 2019: What we know

Chemiluminescence response of human natural killer cells. I. The relationship between target cell binding, chemiluminescence, and cytolysis

Risk factors for COVID-19-related mortality in people with type 1 and type 2 diabetes in England: A population-based cohort study

Sanjeevini: A freely accessible web-server for target directed lead molecule discovery

Iridoids in Physostegia virginiana

Receptor-binding domains of spike proteins of emerging or re-emerging viruses as targets for development of antiviral vaccines

NIST 101. Computational chemistry comparison and benchmark database

Drugs as P-glycoprotein substrates, inhibitors, and inducers

The effect of human mobility and control measures on the COVID-19 epidemic in China

Analysis of physicochemical properties for drugs from nature

SARS-CoV2: Should inhibitors of the renin-angiotensin system be withdrawn in patients with COVID-19?

Structure of the SARS-CoV-2 spike receptor-binding domain bound to the ACE2 receptor

Lead-and drug-like compounds: The rule-of-five revolution

Iridoid glycosides from underground parts of Ajuga remota

COVID-19: Consider cytokine storm syndromes and immunosuppression

Use of antiviral drugs to reduce COVID-19 transmission. The Lancet. Global Health

AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility

Reliable prediction of T-cell epitopes using neural networks with novel sequence representations

The utility of in vitro cytochrome P450 inhibition data in the prediction of drug-drug interactions

Phenolic extracts from Clerodendrum volubile leaves inhibit cholinergic and monoaminergic enzymes relevant to the management of some neurodegenerative diseases

Synthesis, spectral, structure and computational studies of novel transition Metal (II) complexes of (Z)-((dimethylcarbamothioyl) thio)((1, 1, 1-trifluoro-4-(naphthalen-2-yl)-4-oxobut-2-en-2-yl) oxy)

Case-fatality rate and characteristics of patients dying in relation to COVID-19 in Italy

Modeling of intestinal drug absorption: Roles of transporters and metabolic enzymes (for the Gillette Review Series). Drug Metabolism and Disposition: The Biological Fate of Chemicals

COVID-19 diagnosis and management: A comprehensive review

The epidemiology and pathogenesis of coronavirus disease (COVID-19) outbreak

Coronavirus Pandemic (COVID-19). Our World In Data 2020

Python: A programming language for software integration and development

Coronavirus envelope protein: Current knowledge

COVID-19 infection: The perspectives on immune responses

World Health Organization declares global emergency: A review of the 2019 novel coronavirus (COVID-19)

Chloroquine in fighting COVID-19: Good, bad, or both?

A review of drying methods for improving the quality of dried herbs

Potent binding of 2019 novel coronavirus spike protein by a SARS coronavirus-specific human monoclonal antibody

Evaluation of the use of different solvents for phytochemical constituents, antioxidants, and in vitro anti-inflammatory activities of Severinia buxifolia

Reassessment of the polar fraction of Stachys alopecuros (L.) Benth. subsp. divulsa (Ten.) Grande (Lamiaceae) from the Monti Sibillini National Park: A potential source of bioactive compounds

Antiviral activity of medicinal plants of Nilgiris

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

Remdesivir in adults with severe COVID-19: A randomised, double-blind, placebo-controlled, multicenter trial

Coronavirus disease 2019 (COVID-19): Situation report

Efficient implementation of the gauge-independent atomic orbital method for NMR chemical shift calculations

Mortality and burden of disease

Binding of chrysoidine to catalase: Spectroscopy, isothermal titration calorimetry and molecular docking studies

Considering the impact drug-like properties have on the chance of success

Finding the rules for successful drug optimisation

Analysis of Epidemiological characteristics of new coronavirus pneumonia

COVID-19: A recommendation to examine the effect of hydroxychloroquine in preventing infection and progression

Network-based drug repurposing for novel coronavirus 2019-nCoV/ SARS-CoV-2

Epidemic trend of corona virus disease 2019 (COVID-19) in mainland China