Garcinia cambogia Phenolics as Potent Anti-COVID-19 Agents: Phytochemical Profiling, Biological Activities, and Molecular Docking

COVID-19 is a disease caused by the coronavirus SARS-CoV-2 and became a pandemic in a critically short time. Phenolic secondary metabolites attracted much attention from the pharmaceutical industries for their easily accessible natural sources and proven antiviral activity. In our mission, a metabolomics study of the Garcinia cambogia Roxb. fruit rind was performed using LC-HRESIMS to investigate its chemical profile, especially the polar aspects, followed by a detailed phytochemical analysis, which led to the isolation of eight known compounds. Using spectrometric techniques, the isolated compounds were identified as quercetin, amentoflavone, vitexin, rutin, naringin, catechin, p-coumaric, and gallic acids. The antiviral activities of the isolated compounds were investigated using two assays; the 3CL-Mpro enzyme showed that naringin had a potent effect with IC50 16.62 μg/mL, followed by catechin and gallic acid (IC50 26.2, 30.35 μg/mL, respectively), while the direct antiviral inhibition effect of naringin confirmed the potency with an EC50 of 0.0169 μM. To show the molecular interaction, in situ molecular docking was carried out using a COVID-19 protease enzyme. Both biological effects and docking studies showed the hydrophobic interactions with Gln 189 or Glu 166, per the predicated binding pose of the isolated naringin.


Introduction
COVID-19 is a novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). It is a contagious strain of coronavirus that appeared in Wuhan, China in December 2019. The virus is believed to be of animal origin, emerging from the animal market through human−animal contact and causing respiratory and gastrointestinal tract infections with various symptoms starting from fever and general fatigue to shortness of breath, respiratory issues, renal failures, and death [1].
The SARS-CoV-2 virus infection spreads by person-to-person contact with another contaminated person. Serious complications could happen to immune-compromised individuals and the older population, who can develop a critical condition [2]. According to WORLDMETER records (12 September 2022), more than 614 million cases have been reported since the start of this pandemic [3]. WHO also recorded over 5.7 million new Plants 2022, 11, 2521 2 of 10 cases in the first week of July 2022, and over 6.3 million deaths have been reported globally [4]. Billions of COVID-19 vaccines were aligned with drug protocols to fight this global pandemic.
Coronaviruses cause multiple respiratory and intestinal infections in both humans and animals. The 3-chymotrypsin-like cysteine protease (3CL Pro ), known as main protease (M pro ), plays an essential role in the development of the polyproteins that result from the viral RNA translation [5]. Consequently, 3CL Pro inhibitors are a promising choice that can block viral replication and be used to treat COVID-19 patients. The 3CL protease assay used has the advantage in that it is a homogeneous assay with no time-consuming washing steps.
Nowadays, the importance of phytomedicine against COVID-19 has gained much attention by using either plants' crude extracts or their pure isolates. Many publications revealed the role of some plants against different respiratory viruses, including COVID-19, through many mechanisms, such as blocking virus release or even inhibiting enzymes responsible for replication and viral entry inhibition. Examples of these plants and their active metabolites include Pelargonium sidoides, Isatis indigotica, Artemisia annua, and glycyrrhizin from Glycyrrhiza glabra roots [6].
The Guttiferae (Clusiaceae) is a tropical family comprising about 37 genera and 1610 species of trees and shrubs [7]. The genus Garcinia contains a large number of bioactive secondary metabolites. It has received considerable attention for its chemical content in the last few years as its extracts are rich in derivatives of polyphenols, polyisoprenylated benzophenones, xanthones, and bioflavonoids [8]. Phenolic compounds have important functional properties, such as antioxidant, anticancer, and anti-inflammatory capacities, among others [9]. In a recent publication, the garcinoic acid in Garcinia kola suppressed the induced hyper-inflammatory SARS-CoV-2 spike glycoprotein S1 by suppressing the activation of NF-κBin human PBMCs [10].
Garcinia cambogia Roxb. (synonymous to Garcinia cowa Roxb.) is native to southeast Asia and commonly known as Malabar tamarind. It was reported in a recent publication [8] that its fruits are rich in hydroxycitric acid, which has an antiobesity effect, as well as oxy-guttiferone-K, M, and I with their benzophenones' analogues, which are topoisomerase II inhibitors. Benzophenones garcinol, rheediaxanthone-A, and isogarcinol were also identified from its peel and volkensiflavone and fukugetin bioflavonoids from its pericarp.
Previous investigations of the antioxidant, antileukemic, antifungal, antiviral, antibacterial, hypolipidemic, and antiadipogenic activities of the Garcinia cambogia Roxb. plant extracts have been conducted [8]. Biological studies on seed extracts showed that they have vasodilatory, antiulcerogenic, and antihistaminic effects. The gastroprotective effects were assumed to diminish acidity and improve the mucosal defenses [11]. A photograph of the fruit is provided in Figure S1.
The objective of this work was to assess the potential of Garcinia cambogia Roxb. as a source of phenolic compounds, including its chemical characterization and evaluation of antiviral activities against COVID-19.

Metabolomic Analysis
Chemical profiling of the secondary metabolites of the aerial parts of crude ethanolic extract of the Garcinia cambogia Roxb. fruit rind by the use of LC-HRESIMS led to the identification and characterization of 35 various metabolites (1-35) belonging to various chemical classes, among which phenolics, such as flavonoids, catechins, and phenolic acids, dominated. The identification of the resulted metabolites was processed by applying algorithms and macros, which coupled MZmine with both in-house and online databases (METLIN and DNP databases for plant natural products) (Table S1, Figure S2 in the supplementary file).

Biological Activities
The inhibition of SARS-CoV-2 3CL Pro ; IC50 activity by the eight isolated compounds was tested at a final concentration of 50 μM. It was noticed, however, that the incubation of the 3CL Pro with the compound 5 produced a significant reduction of the protease activity compared to the plumbagin positive standard followed by both compounds 6 and 8, while other compounds showed moderate activities (Table 1, Figure S3). Detailed results and standard curves can be found in the supplementary file (Table S2, Figure S4).   (6), p-coumaric acid (7), gallic acid (8).

Biological Activities
The inhibition of SARS-CoV-2 3CL Pro ; IC 50 activity by the eight isolated compounds was tested at a final concentration of 50 µM. It was noticed, however, that the incubation of the 3CL Pro with the compound 5 produced a significant reduction of the protease activity compared to the plumbagin positive standard followed by both compounds 6 and 8, while other compounds showed moderate activities (Table 1, Figure S3). Detailed results and standard curves can be found in the supplementary file (Table S2, Figure S4).

Human Coronavirus (COVID-19) Antiviral Assay
An antiviral evaluation of the isolated compounds against the coronavirus was performed with two concentrations (1 mM and 10 µM); it was noted that compound 5 possessed a potent inhibition activity of COVID-19 in both test concentrations compared to the remdesivir positive standard. This potency forced us to repeat this test on it to evaluate its EC 50 with a serial concentration (10, 1, 0.1, 0.01, and 0.001 µM), which confirmed its anti-COVID efficacy ( Table 2, Figures 2 and S5). Detailed results and standard curves can be found in the supplementary file (Tables S3 and S4, Figure S6).

. Human Coronavirus (COVID-19) Antiviral Assay
An antiviral evaluation of the isolated compounds against the coronavirus was performed with two concentrations (1 mM and 10 μM); it was noted that compound 5 possessed a potent inhibition activity of COVID-19 in both test concentrations compared to the remdesivir positive standard. This potency forced us to repeat this test on it to evaluate its EC50 with a serial concentration (10, 1, 0.1, 0.01, and 0.001 μM), which confirmed its anti-COVID efficacy ( Table 2, Figures 2 and S5). Detailed results and standard curves can be found in the supplementary file (Tables S3 and S4, Figure S6).

Molecular Docking Studies
Docking scores for the natural isolated compounds were provided with the ligand to

Molecular Docking Studies
Docking scores for the natural isolated compounds were provided with the ligand to draw a predictive picture of where these target compounds are binding with the COVID-19 protease (Table 3). A literature survey of the SARS-CoV-2 main proteases revealed that a relative hydrophobic group is essential for the activity to bind to Gln 189 residue, and the presence of Gln 191 and Glu 166 can form the hydrophobic subsite in other SARS-CoV-2 main proteases, which may create extra van der Waal interactions in addition to one hydrogen bond interaction with Thr 190 [14]. In this study, most of the target compounds showed hydrophobic interactions with Gln 189 or Glu 166. As per the predicted binding pose of compound 5, one hydrogen bond donor between the hydroxyl phenolic group and the Thr 190 amino acid residue was recognized in addition to a hydrophobic interaction with the Gln 189 residue. Moreover, compound 8 formed the traditional hydrophobic interaction with Gln 189 in addition to two hydrogen bond interactions with Gln 192 and Glu 166. The extra binding interaction of compound 8 may explain the partial increase of vial inhibition. Compound 7 showed two hydrogen bond interactions via binding hydroxyl and carbonyl groups with Glu 166 and Gln 192 residues, respectively. Finally, compound 6 bound to Thr 26 and His 41 via hydrogen bond and hydrophobic interactions (Figures 3 and 4).

Discussion
Several reports of phenolic natural compounds with anti-SARS-CoV-2 activity indicated metabolites that have been or are currently being studied. Flavonoids are considered natural biologically active phenolic compounds, their antioxidant power due to the presence of hydroxyl functional groups, which enable them to be used as anticancer and antiviral agents by preventing the infection and enhancing immunity. They affect coronaviruses at all stages during their life cycle, starting with cell penetration and the entry inside it, replication of the coronavirus's nucleic acid, and finally, release of the virion outside the cell; they also could have a role with cellular targets of the host [15]. These phenolics could be an essential source in the battle against coronaviruses with advantageous safe administration without systemic toxicity. This work highlights the significance of phenolics, including flavonoids and their potential on various key SARS-CoV-2 targets with molecular docking analysis.
A previous study explored the in vitro activities of over 69 flavonoids and their efficacy against SARS-CoV-2 using many SARS-CoV-2 targets, such as the SARS-CoV-2 3C-like protease (3CL pro ), nucleocapsid (N) protein, the spike (S) protein-ACE2 interaction, papainlike protease (PL pro ), and helicase [15]. The results detailed SARS-CoV-2 3C-like protease (3CL pro ) as the most promising target, explaining the large number of flavonoids' inhibitory activities against this protease. A recent study was performed on various flavonoids for screening their effect against coronavirus infection including quercetin, which showed the capacity to block the entry of SARS-CoV-2 in host cells [16]. Naringin aglycone (Naringenin) previously proved its activity against the hepatitis C virus (HCV) infection [17] and the alphavirus [18], while in our studies, we confirmed a direct role of naringin in the inhibition of the replication of the virus in lung cells pre-infection and even post-infection more than quercetin and other isolated phenolics. In SARS-CoV-2, molecular docking assumed that naringin has the power to inhibit SARS-CoV-2's 3CL-M pro and, thus, prohibit viral replication compared to the standard remdesivir, which is reported as a promising antiviral drug, similar to the ritonavir drug [19]. Many reports of compounds of natural sources with anti-SARS-CoV-2 potential are presently being investigated. Naringin flavonoid is the diglucoside of naringenin aglycone, which had been studied previously and proved to have potential effects against SARS-CoV-2.
Needless to say, all the prospective isolated compounds require in vivo anti-COVID-19 assessment and pharmacokinetic evaluation before going for clinical trial to find useful plant products for COVID-19 therapy, as elaborated in this article. Additionally, a species of Garcinia genus, Garcinia kola, was recently recorded for the treatment of respiratory disorders.

Conclusions
As shown in our study on Garcinia cambogia Roxb. fruit rind, it is rich in secondary metabolites, mainly phenolic compounds. Eight of them were isolated and subjected to in vitro anti-COVID-19 and in silico studies. Naringin was the most promising compound, followed by catechin and gallic acid; these results were further confirmed by the in silico experiment.

Plant Material
Garcinia cambogia Roxb. fruit rinds were purchased from a folk medicine market in Haram, Egypt in February 2022, confirmed by Asmaa Khamis, lecturer of Plant Taxonomy and the Flora-Department of Botany-Faculty of Science-Fayoum University. Specimen no. FuPD-2 was deposited at the department of Pharmacognosy, Faculty of Pharmacy, Fayoum University.

Preparation of Extracts
The dried Garcinia cambogia fruit rind (500 g) was well cleaned and milled by Moulinex caba (type 843, code 243, 220 vac 50 Hz 750 W) at 3000 RCF for 6 min at room temperature to obtain a homogenate powder, followed by extraction by cold maceration with 70% ethanol. The solvent was evaporated under reduced pressure at 40 • C to yield 42 g of dried crude total 70% EtOH extract, which was further used for metabolomic, phytochemical, and biological studies. The total EtOH extract of Garcinia cambogia Roxb. (40 g) was suspended in 250 mL water and successively fractionated using a separating funnel and solvents of increasing polarity: n-hexane (3.5 g), dichloromethane (1.1 g), ethyl acetate (1.2 g), n-butanol (2.2 g), and water (4.25 g). The remaining crude extract (2 g) was used for both metabolomic and biological studies.

Metabolomic Analysis
The total 70% EtOH extract of Garcinia cambogia Roxb. fruit rind under investigation was subjected to metabolomic analysis using LC-HRESIMS with the conditions previously described in the supplementary file, page S2 [20].

Phytochemical Study
A detailed chemical investigation from the plant materials to various chromatographic techniques has been applied to obtain the phenolic compounds to correlate phenolic compounds and biological activity. Details are mentioned in the supplementary file, page S3.

Solvent Systems
The following solvent systems were used for developing the TLC chromatography chromatograms, based on several trials: S 1 : CH 2 Cl 2 /MeOH; 9:1; S 2 : CH 2 Cl 2 /MeOH; 7.5:2.5; S 3 : CH 2  According to the protocol of the manufacturer, the assay performed used SARS-CoV-2 proteases assay kits (M Pro , Catalog: 79955-1, BPS Bioscience, Inc., Allentown, PA, USA), using GC376 (MW = 507.5 Dalton) as the positive control. The assay principle is when the enzyme cleaves the substrate and leads to fluorescence production, measured using a microplate reader at wavelengths 460 for emission and 360 nm for excitation (Tecan Spark, Männedorf, Switzerland). In a 96-well plate, different concentrations of selected test compounds (10 µL) were mixed with 30 µL of diluted protease (15 µg/mL) were added, followed by the incubation of these mixtures at room temperature for half an hour. Afterward, a mixture of 10 µL of the substrate and the reaction buffer was added to the wells to attain a final volume of 50 µL with a final concentration of 40 µM. This final mixture was subjected to incubation at 20 • C for four hours using a fluorimeter (TECAN spark microplate-reading) in order to measure the produced fluorescence.

Coronavirus COVID-19 (Genesig ® Real-Time PCR Assay)
Genesig ® Real-Time PCR assay is an in vitro diagnostic test developed to detect SARSal-Tim viral RNA. A sample was taken from nasopharyngeal and oropharyngeal swabs and sputum, processed under the World Health Organization (WHO) recommendations and according to the manufacturer's instructions (Coronavirus COVID-19 Genesig ® Real-Time PCR Assay Issue 5.0, 2020, Primerdesign Ltd., Hampshire, United Kingdom). The extraction of samples with a volume of 700 µL each was performed using a CE IVD extraction system for the isolation of RNA with the aid of GXT DNA ⁄RNA extraction kits, followed by mixing the resuspended components with the reagents of the PCR reaction platform in an appropriate 96-well plate to perform the experiment. PCR assay results were interpreted and evaluated using the Applied Biosystem ® 7500 Real-Time PCR system.

Molecular Docking Studies
Crystallographic structures of the COVID-19 protease enzyme in combination with its ligand are available in the Protein DataBank (PDB file 5R7Z) [21]. The revision of protein errors and protein structure creation were performed by the structure arrangement process on default rules on the molecular operating environment (MOE). Finally, protein partial charges were calculated by the Gasteiger methodology. Molecular docking analyses were performed with MOE2022 software [22]. The ligand coordinates were designed using the ChemDraw Ultra 11.0 program. Then, their protonation, the correction of the atom, and bond types were defined, followed by the hydrogen atoms addition, and finally minimization performance (MMFF94x, gradient: 0.01) [23]. The COVID-19 protease enzyme docking experiment was also accomplished by redocking the ligand in the PDB file 5R7Z, which was then deleted.
The default Triangle Matcher placement method for docking was processed by choosing a GBVI/WSA dG scoring function to determine the binding free energy of the ligand from a given pose to rank the final poses. The ligand complex with the enzyme possessing the lowest score was the best selection. The ligand redocking with its target discovered an RMSD 0.840 A • , which assured that the ligand bound to a similar pocket and confirmed the reliability of docking parameters.

Statistical Analysis
Biochemical data results have been recorded as means ± SE and then submitted to one-way ANOVA, comparing the means between groups using both the least-significant difference (LSD) and Duncan's multiple range tests. The outcome results could be considered as statistically significant when p-value < 0.05 (Grapg Pad Prism 5).