key: cord-0020344-cynugz51
authors: Sabiu, Saheed; Balogun, Fatai O.; Amoo, Stephen O.
title: Phenolics Profiling of Carpobrotus edulis (L.) N.E.Br. and Insights into Molecular Dynamics of Their Significance in Type 2 Diabetes Therapy and Its Retinopathy Complication
date: 2021-08-11
journal: Molecules
DOI: 10.3390/molecules26164867
sha: 5c49e8fe24cfce7e50d24f2e525ef68ce7066936
doc_id: 20344
cord_uid: cynugz51

Adverse effects associated with synthetic drugs in diabetes therapy has prompted the search for novel natural lead compounds with little or no side effects. Effects of phenolic compounds from Carpobrotus edulis on carbohydrate-metabolizing enzymes through in vitro and in silico methods were assessed. Based on the half-maximal inhibitory concentrations (IC(50)), the phenolic extract of the plant had significant (p < 0.05) in vitro inhibitory effect on the specific activity of alpha-amylase (0.51 mg/mL), alpha-glucosidase (0.062 mg/mL) and aldose reductase (0.75 mg/mL), compared with the reference standards (0.55, 0.72 and 7.05 mg/mL, respectively). Molecular interactions established between the 11 phenolic compounds identifiable from the HPLC chromatogram of the extract and active site residues of the enzymes revealed higher binding affinity and more structural compactness with procyanidin (−69.834 ± 6.574 kcal/mol) and 1,3-dicaffeoxyl quinic acid (−42.630 ± 4.076 kcal/mol) as potential inhibitors of alpha-amylase and alpha-glucosidase, respectively, while isorhamnetin-3-O-rutinoside (−45.398 ± 4.568 kcal/mol) and luteolin-7-O-beta-d-glucoside (−45.102 ± 4.024 kcal/mol) for aldose reductase relative to respective reference standards. Put together, the findings are suggestive of the compounds as potential constituents of C. edulis phenolic extract responsible for the significant hypoglycemic effect in vitro; hence, they could be exploited in the development of novel therapeutic agents for type-2 diabetes and its retinopathy complication.

Diabetes, one of the leading causes of death globally, is a chronic metabolic derangement leading to high levels of glucose in systemic circulation (hyperglycaemia) due to the inability of the body to manage available glucose levels, arising from ineffective or insensitive insulin secretion by islets of Langerhans beta cells of the pancreas [1] . The management of the disease is time-consuming, during which diverse secondary complications (nephropathy, neuropathy, cardiopathy, retinopathy, etc.) culminating in death may set in, if no viable treatment/management therapy is embraced. The prevalence of diabetes continues to rise with population growth, and alongside other non-communicable diseases accounted for 74% of the world's total deaths in 2019 [2] . Uncontrolled hyperglycaemia stimulates excessive free radical generation, culminating in oxidative stress associated with most diabetic complications such as neuropathy, nephropathy and retinopathy [3] .

The total phenolic content of the extract was 96.05 mg/g GAE. Based on the standards used for the HPLC analysis, 11 major phenolic compounds including sinapic acid, cacticin, hyperoside, 1,3-dicaffeoxyl quinnic acid, procyanidin, luteolin-7-O-beta-D-glucoside, rutin, epicatechin, isorhamnetin-3-O-rutinoside, chlorogenic acid and myricetin were identified from the chromatogram of the extract ( Figure 1 , Table 1 ). It is noteworthy that the prominent peaks as observed from the chromatogram with their highest relative abundances were for chlorogenic acid followed by rutin, luteolin-7-O-beta-D-glucoside and epicatechin, suggesting they are the major identifiable phenolic constituents of C. edulis (Figure 1 ). 

The antidiabetic effect of phenolic extract of C. edulis based on its half-maximal inhibitory concentrations (IC 50 ) of 0.06 and 0.75 mg/mL revealed significantly (p < 0.05) higher inhibitory activity against alpha-glucosidase and aldose reductase, respectively, than the standards [acarbose (0.72 mg/mL) and ranirestat (7.05 mg/mL)], except for alpha-amylase, where the extract (0.51 mg/mL) favourably competed with acarbose (0.55 mg/mL) ( Table 2) . Alpha-amylase and alpha-glucosidase found in the pancreas and intestine, respectively, are prominent enzymes involved in the uncontrolled or unending hydrolysis of carbohydrate to glucose. In a diabetic state, the control (with inhibitors) of these enzymes is necessary for the management of hyperglycaemia. Due to the available synthetic inhibitors such as acarbose and miglitol triggering gastro-intestinal dysfunctional events including bloating, flatulence and diarrhoea, the use of plant-based alternative therapeutic agents has been advocated [15] , provided they are able to strongly (low IC 50 score) and mildly (marginally high IC 50 value) inhibit the specific activity of alpha-glucosidase and alpha-amylase [16, 17] , as observed in this study with C. edulis phenolic extract. Although in a related study [14] , among the studied extracts (50% methanol, 70% acetone and aqueous) of the plant, only aqueous extract was tested against alpha-glucosidase, the result obtained in the current study fared comparably in terms of activity with the previously determined aqueous extract, which not only corroborates the antihyperglycemic effect of the plant but the use of polar solvent as medium of extraction or herbal preparation in indigenous medicine. Overall, the exhibited antihyperglycemic effect by the extract could be attributed to its phenolic compounds, as studies have reported phenolic compounds such as cacticin, epicatechin, luteolin glucosides, cyanidin (identified in this study) in plants such as Bergenia ciliata, Diospyros kaki and Eleusine coracana as inhibitors of alpha-amylase and alpha-glucosidase [7, 18] . Aldose reductase, an enzyme implicated in the complications of diabetes, particularly retinopathy, facilitates sorbitol production from glucose through polyol pathway during the glucose catabolism [19] . Drugs such as ranirestat, sorbinil and alrestatin act by improving the reduced glucose concentration, thus facilitating eventual absorption in tissues (such as neural, lens and glomeruli) [20] . However, like many other synthetic drugs, they exhibit side effects including nausea, fever and diarrhoea, hence, the need for alternative therapeutic molecules (inhibitors) from natural products in drug design. Based on its IC 50 value with respect to modulatory effect of inhibitors on aldose reductase, the elicited activity by C. edulis phenolic extract in this study could be indicative of its prospect in the management of long-term diabetic retinopathy complication [21] [22] [23] [24] [25] . Patel et al. [26] also advocated phenolic compounds such as luteolin-7-O-β-D-glucopyranoside and 4,5-di-O-caffeoylquinic acid which are similar compounds identified in this study, as inhibitors of aldose reductase with prospect for diabetes retinopathy. 

To gain insight into the probable interactions between the identified phenolic compounds (as revealed by the HPLC analysis) and the study enzymes in this study, computational evaluation was performed through molecular docking and MDS. Molecular docking, a measure of fitness and pose of a compound at the active site of an enzyme, normally gives high negative scores as a reflection of better pose of the compound [27] . In this study, phenolic compounds such as 1,3-dicaffeoxyl quanic acid, chlorogenic acid, epicatechin, luteolin-7-O-beta-D-glucoside, isorhamnetin-3-O-rutinoside, myricetin and rutin, had significant and better poses based on their scores than the reference standard, ranirestat, when docked with aldose reductase (Table 3) . Additionally, better poses were observed for epicatechin, luteolin-7-O-beta-D-glucoside, isorhamnetin-3-O-rutinoside, rutin, hyperoside and procyanidin with alpha-amylase compared to the resulting complexes with acarbose (Table 3) . While most of the identified compounds including 1,3-dicaffeoxyl quanic acid, chlorogenic acid, epicatechin, isorhamnetin-3-O-rutinoside, luteolin-7-O-beta-D-glucoside, myrcetin, rutin, cacticin, hyperoside and procyanidin showed good docking with alpha-glucosidase as depicted by the higher negative values than acarbose, other compounds (epicatechin, isorhamnetin-3-O-rutinoside, chlorogenic acid and rutin) had commendable binding at the active sites of the three enzymes (Table 3) , which is indicative of their prospective interaction with the enzymes [28] . However, since docking is only a preliminary reflection of the ligand's fitness within the binding pocket of a receptor, the binding orientations of the studied phenolics were subjected to further binding energy calculations and MDS. Typically looking at the thermodynamic calculations against alpha-amylase, procyanidin among other compounds had the highest (−69.834 kcal/mol) binding energy, which was better than the value for acarbose (−54.679 kcal/mol) and rutin (−46.826 kcal/mol) (Table 4) . Similarly, against alpha-glucosidase, 1,3-dicaffeoxyl quinic acid and hyperoside had higher binding energies than acarbose, while isorhamnetin-3-O-rutinoside by luteolin-7-O-beta-D-glucoside and rutin had higher binding energies than ranirestat against aldose reductase (Table 4 ). Higher negative values are indicative of stronger affinity of these compounds with the respective enzymes and hence possible better stability of the resulting complex [29] . While potential stronger affinities of phenolic compounds (over synthetic inhibitors) in complex with antidiabetic enzymes have been reported [30, 31] , Rasouli et al. [18] observed higher binding free energy of some of the phenolic compounds (rutin, cacticin and epicatechin) reported in this study against alpha-amylase as well as alpha-glucosidase. Table 3 . Molecular docking scores (kcal/mol) of the phenolic compounds with carbohydrate hydrolyzing enzymes and aldose reductase.

Aldose Reductase α-Amylase α-Glucosidase The enzyme-ligand complex is prone to conformational changes inducible by the binding ligand, and this could result in possible alteration of the biological activity of the enzyme [31] . Based on the affinity of the enzymes for the phenolics through 100 ns MDS in this study, a further probe into understanding the degree of stability and compactness of the resulting complexes was undertaken. To examine the stability while ensuring the equilibration of the enzyme-ligand, parameters including RMSD, RMSF and RoG evaluated over an extended period of 100 ns for unbound (apo-enzyme, i.e., alpha-amylase, alphaglucosidase and aldose reductase) and bound complexes (enzyme + inhibitor, i.e., standards and phenolic compounds) are presented in Figures 2-4 .

The nature/extent of stability and convergence of a ligand-receptor system is normally measured by RMSD [28] . In this study, the average RMSD values for procyanidin, rutin and acarbose (standard drug) were 1.62 Å, 1.68 Å and 1.82 Å, respectively, while it was 2.04 Å for α-amylase (apo-enzyme) ( Figure 2A ). With these findings, RMSD of the compounds and standard were lower than the apo-enzyme, meaning that both the compounds and standard possessed the ability to promote increased structural stability of the alpha-amylase. Additionally, the lower RMSD of the phenolic compounds compared to the standard is suggestive of their better advantage as prospective leads and possible novel inhibitors of the enzyme. Figure 2B shows higher average RMSD bound systems for the phenolic compounds, hyperoside (2.22 Å), 1,3-dicaffeoxyl quinic acid (1.73 Å) and luteolin-7-O-beta-D-glucoside (1.70 Å) relative to the apo-enzyme, alpha-glucosidase (1.48 Å). However, it must be noted that these compounds, particularly 1,3-dicaffeoxyl quinic acid and luteolin-7-O-beta-Dglucoside, had RSMD values that are not only within acceptable limits, but lower than that of the standard, acarbose (1.79 Å), pointing to their relative superiority over the reference standard. Interestingly, these results are in tandem with the thermodynamic profiles where 1,3-dicaffeoxyl quinic acid had the best binding free energy (Table 4) , synonymous with better compactness and structural stability of the complex. A similar trend was also observed where the RMSD value of the apo-enzyme was lower than the bound systems regarding aldose reductase. Typically, aldose reductase average RMSD (1.32 Å) was lower than ranirestat system (1.38 Å) which was not significant but only marginally lower than the value (1.43 Å) observed with luteolin-7-O-beta-D-glucoside ( Figure 2C ). This observation may not infer the inability of luteolin-7-O-beta-D-glucoside to promote the structural stability of the complex going by comparable RMSD value with ranirestat. Aside from the fact that the value is within acceptable limit, the effect elicited by luteolin-7-O-beta-D-glucoside with respect to its RMSD value corroborates its enhanced binding free energy in comparison with other compounds and the reference standard (Table 4 ). Radius of gyration (RoG) determines the total compactness of the enzyme-inhibitor binding [28, 32] . It is a measure of densification of the protease structure [33] , and the smaller the RoG value, the better the protease stability. In line with RMSD result, the lead compounds and standard drug revealed mean RoG values of 23.24 Å (procyanidin), 23.25 Å (rutin) and 23.37 Å (acarbose) lower than the apo-enzyme (23.54 Å), indicating that the binding of the compounds potentially stabilized alpha-amylase better than the control molecule ( Figure 3A) . However, the RoG results of compounds and standard drugs for alpha-glucosidase and aldose reductase do not follow the trend of RMSD, as there were reductions in RoG values of phenolic compounds such as 1,3-dicaffeoxyl quinic acid (27. The binding property of the inhibitor or ligand and the active site residues of each protein was further evaluated by RMSF. Increased or decreased fluctuations are sin qua non to high or low flexibility movement or interaction between ligands and the receptor amino acids residues [28] . In the finding for alpha-amylase system, rutin (2.79 Å) followed by acarbose (2.54 Å) exhibited the highest average RMSF values, while the lowest value was found with procyanidin (2.05 Å) among the studied interactions. While it was observed that compounds and the standard drug increased the enzyme (1.90 Å) fluctuation or amino acid residue flexibility, a kind of similar pattern of fluctuations was seen among the compounds, the standard drug and enzyme at 200, 325 and 350 residues ( Figure 4A ). Except for luteolin-7-O-beta-D-glucoside (1.88 Å), compounds including hyperoside (4.31 Å) and 1,3-dicaffeoxyl quinic acid (3.24 Å) were found to have higher average RMSF above the enzyme (3.06 Å). The observed fluctuations were seen around 350, 425 and 800 residues ( Figure 4B ). The highest RMSF in the aldose reductase system was 2.88 Å (standard drug), while the lowest for the studied interactions was 1. The interaction between the binding of molecules (ranirestat, acarbose) or compounds with the active site residues of the enzymes (alpha-amylase, alpha-glucosidase and aldose reductase) is represented by ligand-enzyme interaction plots (Figures 5-7) . The interactions between acarbose (standard), procyanidin and rutin on the active sites of alpha-amylase from the plots ( Figure 5A-C) were Van der Waals forces, hydrogen (to hydrogen) bonds, donor-donor interaction, C-H bond, π-π stacked interaction and π-alkyl bonds, though the number of these interactions differs between molecules and observed to be a consequence of their binding free energies. While acarbose Van der Waals forces (with Gln403, Phe405, Val400, Pro404, Thr332, Thr10, Phe334), H-H bonds (Arg397, Arg420, Gly333, Pro331, Trp279, Asn278) and C-H (Asp401, Lys277, Gly280, Ser288) formed, it had one interaction fewer compared with procyanidin (21 bonds) [Van Der Waal (Arg302, Ala306, Gly305, Leu236, Hie200, Leu161, Phe255, Glu232, Ala197), H-H bonds (Asp299, Asp196, Arg194, Trp58, Asp355, Hie298, Asn297), π-π (Trp57, Tyr61) and π-alkyl (Hie304, Ile234)], and this could explain why it had lower binding free energy in comparison with procyanidin (Table 4 ), which possessed more numbers of hydrogen bonds and presence of π-π stacked interaction and π-alkyl bonds. The binding free energy capacity of rutin (lower than acarbose and procyanidin) is corroborated by its number of molecular interactions [ (17) including Van Der Waals (Hie298, Hie200, Tyr61, Gly305, Leu164, Val97), H-H bonds (Gln62, Asp299, Asp196, Hie100, Hie304, Tyr150) and π-alkyl bonds (Leu161, Ala197, Trp58, Trp57)]. In terms of amino acid residues involved in the stability, it was observed that Trp57, Trp58, Try61, Leu162, Asp196, His201, Asp299 and Ala197 are the most important amino acid residues involved with compounds (procyanidin and rutin) at the active sites of alpha-amylase. Though these residues are absent in acarbose, our report agrees with the submission of Hashim et al. [34] , where Trp57, Trp58 and His201 have also been identified as important (catalytic) residues involved in alpha-amylase (1DHK) stability. 1,3-Dicaffeoxyl quinic acid [(Ala177, Asp511, Tyr186, Phe544, Tyr410, Ile339, Asp300, Trp272, Trp375, Lys449), (Asp175, Arg475, Asp412, Ile301) (Phe419), (Met413)] and hyperoside [(Arg613, Phe623, Phe625, Thr624, Pro626, Gly700, Gly664, Asn665, Ser727, Hie729), (Asp627, Glu244, Glu699, Arg642), (His698), (Val730) had the same number of interactions (17) with the active sites of alpha-glucosidase and are characterized by (contain the same number of) Van der Waal forces (10), H-bonds (4), π-π stacked interaction (1) and π-alkyl bonds (1); however the highest binding free energy found with 1,3-dicaffeoxyl quinnic acid may be attributed to unidentified carbon-H bonds (Ile176) and formed π-cation (Arg663) in hyperoside. In fact, the presence of π-cation in hyperoside may also be suggested to be the reason for lesser binding energy, as similarly witnessed in acarbose (Glu405) with far less binding energy and lacking π-π stacked, π-alkyl bonds and a lower number of Van der Waal forces (Gly157, Gly158, Ser177, Thr178, Cys176, Val407) ( Figure 6) . Similarly, the interactions [H-bonding (Leu303, Leu304, Leu305), vVn Der Waal forces (Lys224, Arg299, Val300, Ala302, Cys301, Cys306, Gly131, Tyr51), π-sulfur (Trp222), π-Alkyl (Phe125, Leu127) of ranirestat and standard molecule (14) with active sites of aldose reductase is lesser than those of isorhamnetin-3-O-rutinoside, rutin and luteolin-7-O-beta-D-glucoside exhibited in terms of number of interactions (20, 20 and 15 respectively) relative to the former (Figure 7) , and these interactions corroborated the findings from the binding free energies (Table 4) . It is interesting to note that although isorhamnetin-3-O-rutinoside and rutin revealed same number of interactions (20) , the presence of higher numbers of Van der Waal forces [(12) (Pro221, Leu304, Cys301, Ser305, Leu127, Tyr51, Tyr212, Ala48, Val50, Trp82, Phe124, Trp114)], hydrogen bonds [(5) (Lys24, Ala302, Val300, Trp23, Hie113)] and absence of π-cation bond for isorhamnetin-3-O-rutinoside as against 11 (Ser213, Val50, Trp82, Asn163, Phe125, Tyr51, Ala302, Val300, Lys80, Gly21, Asp46), 2 (Hie113, Trp114) and 1 (Lys24), respectively, for rutin could be responsible for the higher binding free energy witnessed in the former over the latter. Amino acids residues such as Leu303, Leu304, Cys301, Cys306, Val300, Phe125, Trp22, Ser305, Tyr51, Ala302 and Arg299 were found surrounding acarbose. While some of these amino acids residues, including Leu303, Cys301, Val300, Ala302 might be believed to be involved in the stability of the ligands-enzyme complexes because they are also shared by the promising compounds, others like Leu304 (rutin), Cys306 (all compounds), Tyr51, (luteolin-7-O-beta-D-glucoside), Phe125 (luteolin-7-O-beta-D-glucoside), Ser305 (luteolin-7-O-beta-D-glucoside) and Trp222 (rutin) are observed to be absent on the interaction plots of respective complexes. Additionally, and more importantly, the presence of π-π stacked interactions between amino acid residues [Trp23 and Tyr212 (rutin), Trp23 (luteolin-7-O-beta-D-glucoside) as well as Trp222 and Phe125 (isorhamnetin-3-O-rutinoside)] and aromatic rings on the compounds absent in ranirestat may be suggestive of its low free energy binding. The presence of π-π stacked interactions in receptor-ligand binding has been reported to be germane in the development of drugs [35] . Hence, it could be inferred that the presence of π-π stacked interaction in these compounds could be advantageous in the development of potential inhibitors against carbohydrate enzymes that trigger diabetes development and its related complication (retinopathy). 

Acarbose, ranirestat, porcine pancreatic α-amylase, rat intestine acetone powder, aldose reductase, extra pure starch, dinitrosalicylic acid (DNS), p-nitrophenyl glucopyranoside (pNPG) and 2-chloro-4-nitrophenyl α-D-maltotrioside were obtained from Sigma-Aldrich, St. Louis, MO, USA. All other chemicals and reagents used are of analytical grade. Carpobrotus edulis leaves collected from the Agricultural Research Council-Vegetables, Industrial and Medicinal Plants campus in Pretoria, South Africa with voucher specimen Mulaudzi RB# 200 deposited in Bews Herbarium, University of KwaZulu-Natal, as described by Mulaudzi et al. [14] were lyophilised and ground into fine powders using a rotor mill (Fritsh Pulverisette 14, Labotec, Midrand, South Africa).

The method of Mulaudzi et al. [14] with slight modification was adopted for the extraction of the powdered materials from the leaves of the plant (10 g) in 50% methanol (250 mL) under sonication in a cold water-containing water bath for 120 min. The extract was filtered through a Whatman No. 1 filter paper, centrifuged (3000 rpm, 15 min), and the resulting supernatant concentrated in a water bath at 40 • C. The dried extract obtained was kept airtight and refrigerated (10 • C) in a glass vial prior to use for total phenolic quantification, HPLC analysis and in vitro antidiabetic assays.

For total phenolic quantification, 50 µL (1 mg/mL) of the phenolic extract was added to 6.950 mL distilled water in a test tube, and gently shaken prior to the addition of Folin-Ciocalteau phenol reagent (0.5 mL) and sodium carbonate (1.5 mL; 20%). Subsequently, distilled water (1 mL) was added to the mixture, shaken vigorously and allowed to stand for 45 min prior to absorbance measurement at 760 nm. The standard (gallic acid) was prepared (0.8 mg/mL in 40% methanol) and its varying concentrations were treated in a similar manner as the phenolic extract [36] . The phenolic content of the extract was estimated from gallic acid standard curve and expressed as milligram per gram gallic acid equivalent (mg/g GAE).

The HPLC analysis was achieved based on the method of Peng et al. [37] with modifications. This was conducted using HPLC (Shimadzu Prominence-i LC-2030C 3D plus, Kyoto, Japan) coupled to a diode array UV detector (HPLC-DAD) and a high-resolution mass spectrometry (HPLC-HRMS) on an Ultimate 3000 RSL Cnano system (Thermo Scientific, Waltham, Massachusetts, United States of America). The mobile phase consisted of A (0.1% formic acid) and B (acetonitrile), and the flow rate was 0.25 mL/min with the temperature of the column (Sunfire C18, 5 µm, 4.6 mm × 150 mm, Waters Corporation, Milford, Massachusetts, United States of America) set at 35 • C and the sample volume maintained at 20 µL. The elution gradient varied from 1% A to 2% B linearly for 2 min and from 2-100% B in 50 min and thereafter, from 10% to 2% for 1 min and from 2% to 0% for 9 min. The chromatogram was based on photo diode array UV detector (DAD) with wavelengths spanning 190-800 nm based on the peak absorption of the analysed compounds. The identification of the compounds was achieved based on their individual retention times and MS fragment patterns compared with those of the standard phenolics (sinapic acid, cacticin, hyperoside, 1,3-dicaffeoxyl quinic acid, procyanidin, rutin, epicatechin, isorhamnetin-3-Orutinoside, chlorogenic acid, myricetin and luteolin-7-O-beta-D-glucoside) used in tandem with published data.

2.6. In vitro Assays 2.6.1. Alpha-Amylase Inhibitory Assay Using a previously reported protocol [38] , the activity of the extract against α-amylase was evaluated. Aliquots (0.1 mL) of either acarbose (reference standard) or the phenolic extract at varying concentrations (0.065-1.000 mg/mL) were added to α-amylase solution (0.1 mL of 0.5 mg/mL). Following a 10-min pre-incubation (25 • C) of the resulting solution, 1% starch solution in 0.02 M sodium phosphate buffer, was added and further incubated (25 • C, 10 min), before the reaction was halted by DNS (0.5 mL). The resulting mixtures were boiled (100 • C, 5 min) and subsequently cooled (25 • C) before final dilution (distilled water, 7.5 mL) and spectrophotometric absorbance reading (540 nm) (OPTIZEN POP, Apex Scientific, Yuseong-gu, Daejeon, Republic of Korea). The results presented as IC 50 (halfmaximal inhibitory concentration) value in each case was non-linearly extrapolated from maltose standard calibration curve.

For this assay, 50 µL of varying concentrations (0.065-1.000 mg/mL) of either acarbose or phenolic extract were added to 0.1 mL αglucosidase (1 M) before incubation (25 • C, 10 min). Thereafter, 0.05 mL of 5 mM p-NPG solution was added followed by a further 5 min incubation (25 • C) period before 0.05 mL Na 2 CO 3 (0.1 M). Following this treatment, a microplate reader (MULTISKAN GO 1519 Thermo Scientific, Vantaa, Finland) was used to take the absorbance readings (405 nm) and the IC 50 values were similarly non-linearly determined from a p-nitrophenol standard curve [38] .

In this assay, glyceraldehyde and NADPH were used as substrate and cofactor, respectively [39] . The final concentration of the dissolving solvent was kept equivalent among reaction mixtures in the presence of varying concentrations of either the phenolic extract or ranirestat (reference standard). The rates of reaction were monitored spectrophotometrically (340 nm) at 25 • C and compared with the control not containing the phenolic extract, and the IC 50 value was non-linearly determined from a calibration curve. Ranirestat was used as reference standard and the experiments were performed in triplicate.

The collection of crystal (X-ray) structure of the enzymes [PDB: 3RX3 (aldose reductase), 3W37 (α-glucosidase), and 1DHK (α-amylase)] were from the RSCB Protein Data Bank (https://www.rcsb.org/ accessed on 12 December 2020). The UCSF Chimera software V1.14 was used in the preparation of the enzymes in readiness for docking [40] , PubChem (https://pubchem.ncbi.nlm.nih.gov/ accessed on 15 December 2020) was used to retrieve the structures of the chromatogram-identified phenolic compounds (sinapic acid, cacticin, hyperoside, 1,3-dicaffeoxyl quinic acid, procyanidin, rutin, epicatechin, isorhamnetin-3-Orutinoside, chlorogenic acid, myricetin and luteolin-7-O-beta-D-glucoside) and standards (acarbose and ranirestat) and optimization of their three-dimensional structures executed using Avogadro software as previously reported [41] . The optimized compounds (ligands) and the enzymes were subsequently subjected to molecular docking.

The docking of the prepared phenolic compounds and standards into binding pockets of the enzymes (α-amylase, α-glucosidase, and aldose reductase) was by Autodock Vina Plugin on Chimera V1.14. Judging by the docking scores, complexes identified to have the best pose for each compound were ranked, selected and further analyzed through 100 ns molecular dynamics simulation (MDS).

The MDS was achieved as recently reported [28] , using the GPU (force fields) version obtainable in AMBER package, where the description of the system by FF18SB variant of the AMBER force field was carried out [42] . With the aid of Restrained Electrostatic Potential (RESP) and the General Amber Force Field (GAFF) methods of the ANTECHAMBER assisted with information on atomic partial charges for the compounds. Hydrogen atoms and Na+ and Cl-counter ions (to neutralize the system) were made possible with Leap module of AMBER 18. The residues were numbered 1-336, 913, and 496, respectively, for aldose reductase, α-glucosidase and α-amylase. The system in each case was then lowered implicitly within an orthorhombic box of TIP3P water molecules such that all atoms were within 8Å of any box edge. MDS total time carried-out were 100 ns. For each simulation, hydrogens atoms were constricted using the SHAKE algorithm. The step size of each simulation was 2 fs, and an SPFP precision model was used. The simulations align with the isobaric-isothermal ensemble (NPT), having randomized seeding, Berendsen barostat maintains 1 bar constant pressure, 2 ps pressure-coupling constant, 300 K temperature and Langevin thermostat with a collision frequency of 1.0 ps [43] .

Using PTRAJ, the systems were subsequently saved, and each trajectory analyzed every 1 ps, and the RoG, RMSF, and RMSD were analyzed with CPPTRAJ module (AMBER 18 suit).

Molecular Mechanics/GB Surface Area method (MM/GBSA) was adopted to assess the free binding energy while comparison of the systems binding affinity followed afterwards [44] . Binding free energy was averaged over 100,000 snapshots extracted from the 100 ns trajectory. The ∆G for each system (enzyme, complex and phenolics) was estimated as earlier reported [45] .

For the in vitro experiments, data analyses were carried out by Graph pad Prism version 3.0 using t-test (and nonparametric tests), supplemented with Mann-Whitney test. Results are expressed as mean ± standard error of the mean (SEM). Except otherwise stated, the raw data plots for the in silico evaluations were generated using the Origin data analysis software V18 (OriginLab, Northampton, MA, USA) (Seifert, 2014).

While the in vitro studies result gave an insight into possible antidiabetic potential of C. edulis, the HPLC analysis suggested and identified 11 phenolic compounds, which were further analysed as probable hypoglycaemic candidates through in silico studies. Based on the findings from the binding free energy, structural stability and compactness in this study, procyanidin was a better inhibitor of alpha-amylase, 1,3-dicaffeoxyl quinic acid against alpha-glucosidase while luteolin-7-O-beta-D-glucoside showed good inhibitory potentials of aldose reductase among other phenolic compounds. Thus, these molecules could be exploited in developing novel therapeutic candidates against postprandial hyperglycaemia and diabetic retinopathy. 

The data presented in this study are available in the article.

Pathophysiology of type 2 diabetes mellitus

Available online: www.who.int/news-room/fact-sheets/detail/the-top-10-causes-ofdeath

Membrane stabilization and probable mechanisms of hypoglycemic activity of fruit extract of Solanum incanum L. (Solanaceae)

The scope of phytotherapy in southern African antidiabetic healthcare

Inhibition of α-glucosidase, α-amylase, and aldose reductase by potato polyphenolic compounds

Aldose reductase inhibitory activity and antioxidant capacity of pomegranate extracts

Anti-diabetic potential of phenolic compounds: A review

Promising potential of dietary (poly) phenolic compounds in the prevention and treatment of diabetes mellitus

Effect of different commercial fertilizers, harvest date, and storage time on two organically grown blackberry cultivars: Physicochemical properties, antioxidant properties, and sugar profiles

Traditional uses, phytochemistry and pharmacological activity of Carpobrotus edulis: A global perspective

Polyphenols and their effects on diabetes management: A review

Carpobrotus edulis (L.) L.Bolus. SANBI PlantZAfrica

Antidiabetic, anti-inflammatory, anticholinesterase and cytotoxicity determination of two Carpobrotus species

In vitro antioxidant and antidiabetic activities of zinc oxide nanoparticles synthesized using different plant extracts

Aqueous roots extract of Dicoma anomala (Sond.) extenuates postprandial hyperglycaemia in vitro and its modulation against on the activities of carbohydrate-metabolizing enzymes in streptozotocin-induced diabetic Wistar Rats

Chemical group profiling, in vitro and in silico evaluation of Aristolochia ringens on α-amylase and α-glucosidase activity. Evid. Based Complement Altern

Differential α-amylase/α-glucosidase inhibitory activities of plant-derived phenolic compounds: A virtual screening perspective for the treatment of obesity and diabetes

Molecular docking evaluation of some natural phenolic compounds as aldose reductase inhibitors for diabetic complications

Aldose reductase inhibitors

Aldose reductase inhibitors from natural sources

Rat lens aldose reductase inhibitory constituents of Nelumbo nucifera stamens

Aldose reductase, still a compelling target for diabetic neuropathy

Extraction and identification of three major aldose reductase inhibitors from Artemisia montana

Antidiabetic properties of dietary phenolic compounds: Inhibition effects on α-amylase, aldose reductase, and α-glycosidase

Pharmacologically tested aldose reductase inhibitors isolated from plant sources-A concise report. Chin

Molecular dynamic mechanism(s) of inhibition of bioactive antiviral phytochemical compounds targeting cytochrome P450 3A4 and P-glycoprotein

Repurposing drugs and identification of inhibitors of integral proteins (spike protein and main protease) of SARS-CoV-2

Inhibitory evaluation of Curculigo latifolia on αglucosidase, DPP (IV) and in vitro studies in antidiabetic with molecular docking relevance to type 2 diabetes mellitus

Identification of phenolic compounds from Nettle as new candidate inhibitors of main enzymes responsible on type-II diabetes

Biological pretreatment of lignocellulosic biomass-An overview

Convergence of sampling in protein simulations

α-Glucosidase inhibitory effect of anthocyanins from Cinnamomum camphora fruit: Inhibition kinetics and mechanistic insights through in vitro and in silico studies

Antioxidant and α-amylase inhibitory property of Phyllanthus virgatus L.: An in vitro and molecular interaction study

Molecular recognition of protein-Ligand complexes: Applications to drug design

Analysis of total phenols and other oxidation substrates and antioxidants by means of Folin-Ciocalteu reagent

Rapid and sensitive detection of fipronil and its metabolites in edible oils by solid-phase extraction based on humic acid bonded silica combined with gas chromatography with electron capture detection

Potent α-glucosidase and α-amylase inhibitory activities of standardized 50% ethanolic extracts and sinensetin from Orthosiphon stamineus Benth as anti-diabetic mechanism

Kinetic and spectroscopic evidence for active site inhibition of human aldose reductase

Structure-based ligand design to overcome CYP inhibition in drug discovery projects

PubChem substance and compound databases

Molecular dynamics simulations: From structure function relationships to drug discovery

A quadratically convergent SHAKE in

MMGBSA as a tool to understand the binding affinities of filamin-peptide interactions

An insight on the nature of biochemical interactions between glycyrrhizin, myricetin and CYP3A4 isoform

The authors declare no conflict of interest.