In silico and in vivo characterization of cabralealactone, solasodin and salvadorin in a rat model: potential anti-inflammatory agents
Abstract
Background
The present study investigates the hepato- and DNA-protective effects of standardized extracts of Cleome brachycarpa (cabralealactone), Solanum incanum (solasodin), and Salvadora oleioides (salvadorin) in rats.
Materials and methods
Hepatotoxicity was induced with intraperitoneal injection of carbon tetrachloride (CCl4) (1 mL/kg b.wt.) once a week for 12 weeks. The hepato- and DNA protective effects of the extracts in different combinations were compared with that of a standard drug Clavazin (200 mg/kg b.wt.). Tissue alanine aminotransferase, alpha-fetoprotein, tumor necrosis factor alpha (TNF-α), isoprostanes-2α, malondialdehyde, and 8-hydroxydeoxyguanosine, the significant hallmarks of oxidative stress, were studied.
Results
Histopathological findings of the liver sections from the rat group which received CCl4+cabralealactone, solasodin, and salvadorin demonstrated improved centrilobular hepatocyte regeneration with moderate areas of congestion and infiltration comparable with Clavazin. For in silico study, the identified compounds were subjected to molecular docking with cyclooxygenase-2 and TNF-α followed by a molecular dynamics study, which indicated their potential as anti-inflammatory agents.
Conclusion
Cabralealactone, solasodin, and salvadorin confer some hepatoprotective and DNA-damage protective effects against CCl4-induced toxicity. They successfully restored the normal architecture of hepatocytes and have the potential to be used as inhibitor to main culprits, that is, cyclooxygenase-2 and TNF-α. They can combat oxidative stress and liver injuries both as mono and combinational therapies. However, combination therapy has more ameliorating effects.
Introduction
Currently, more than 10% of humans are affected by continuous exposure to hepatotoxic mediators like alcohol, viruses, toxic substances, and biotransformed metabolites, which may contribute to the development of liver cirrhosis and hepatocellular carcinoma.1 Liver is an organ, which is at higher risk for toxic metabolic reactions that lead to chronic liver inflammation, mainly associated with the innate defense system of the biological organization for elimination of toxic metabolites. Elevation of liver enzymes (alanine aminotransferase [ALT], aspartate aminotransferase, gamma-glutamyl transferase, and alkaline phosphatase) is a useful indicator of an abnormal hepatic functioning in addition to the clinical manifestations seen.2 Since hepatotoxicity is associated with inflammation and oxidative stress, herbal extracts with potent anti-inflammatory and anti-oxidative properties are known to have therapeutic values in the treatment of fibrosis and toxic injuries. Therefore, the herbal extracts may be potential sources of new chemical metabolites with good therapeutic potential in xenobiotic-induced liver injuries.1,2
Cleome brachycarpa belongs to the Cleomaceae family and is locally known as Ponwar. It grows widely in sandy regions of Africa (especially in North Africa) and Asia (particularly in Egypt, Saudi Arabia, Afghanistan, Pakistan, India, and Iran). The plant is also known as spider flower because of its flower shape.3C. brachycarpa is reported to have some amazing biochemical activities in the biological system such as antimicrobial, hepatoprotective, anticancer, antiviral, analgesic, antiemetic, immunomodulatory, and antioxidant properties.4 It is also used in the Iranian traditional medicine in the treatment of jaundice, diarrhea, bronchitis, scabies, malaria, heart failure, immune system deficiency, seizure and rheumatoid arthritis, fever, chills, and stomach disorders, although in high doses, it has been reported to induce stomach swelling.5,6 Moreover, in a recent study, the ethanolic extract of its stems, roots, and leaves shows some antioxidant effects.7 The total phenolic and flavonoid contents of these parts are also reported to be active against 1,1-diphenyl-2-picrylhydrazyl (DPPH) radical scavenging and ferric reducing ability of plasma activities.8 Another research established in vitro antibacterial effects of aerial parts,9 memory-boosting effects in Alzheimer,10 and anxiolytic effects in the central nervous system.7 Aerial parts of plant may show cytotoxic effects because of its constituents i.e., flavonoid and sesqiterpene.11 Spectroscopic analysis reveals that 90% of its ethanolic extract contains active phytochemical constituents, including cabralealactone (5.1%), trinortriterpenoid dilactone (0.6%), deacetoxybrachycarpone (1.1%), and ursolic acid (0.2%).7,12 In 2010, the anticancer activity of ethanolic extract of C. brachycarpa was evaluated by using albino mice against Ehrlich Ascites Carcinoma cell line in 2 doses (200 and 400 mg/kg) for 9 consecutive days. The anticancer effect was evaluated by observing the tumor volume, viable tumor cell count, and hematological parameters. All the above-mentioned parameters showed a significant (p<0.01) decrease and the life span of the tumor-bearing host showed an increase, indicating that the plant extract has a potential anticancer activity.11
Solanum incanum is locally known as thorn apple or bitter apple, an important plant utilized in Pakistan. It grows abundantly in Africa and is also found in the Middle East and India. All the parts of the plant are valued since they are mainly used to combat liver pain, ulcer, and benign tumors.13 Isolation of solasodin metabolite from S. incanum berries indicates the presence of a chemical structure related to that of a steroidal glycol-alkaloid, which can confer a protective effect against carbon tetrachloride (CCL4)-induced hepatic insult. Cytotoxic investigation of solasodin indicates the suppressive and apoptotic characteristic against hepatitis (specifically HepG2 and Hep3B) cell lines.14 The anti-cancer effect of the herb is mainly attributed due to its glycosidal alkaloids.15–17 It has been reported that solasodin causes necrosis of cancer cells via upregulation of tumour necrosis factor (TNF) receptor and activation of intrinsic mitochondrial apoptotic signaling pathways.18,19 Studies have reported the hepatoprotective effect of fresh tuber roots of S. incanum on CCl4-induced liver damage.20
Salvadora oleoides or locally known as “Mitha Jall or Peelu” is an evergreen plant found in India and Pakistan, as well as in Southern Iran. It is traditionally used in the treatment of enlarged spleen, rheumatism, low fever, snake bites, and cough, alongside its widespread use as a purgative.21 Two compounds, octadecanoyl heptanoate and δ-lactone of 3 isobutyl-5 hydroxy-19-oxononadeconic acid have been isolated from its stem bark, and salvadorin, a new dimeric dihydroisocoumarin has been isolated from chloroform fraction of S. oleoides.
Gas chromatography and mass spectrometry analysis of essential oils based on analysis of hydro-distilled extracts indicated the presence of 35 (94.0%) and 25 (91.1%) compounds in the leaf and stem, respectively. The leaf was found to contain high concentrations of 2-methoxy-4-vinylphenol (25.4%), cis-3-hexenyl benzoate (16.8%), phytol (13.9%), n-hexadecanoic acid (6.9%), and trans-β-damascenone (2.1%). On the other hand, the stem oil contains high concentrations of 2-methoxy-4-vinylphenol (21.6%), phytol (12.9%), n-hexadecanoic acid (3.6%), octacosane (7.9%), nonacosane (7.3%), 1-octadecene (5.8%), heptacosane (5.9%), hexacosane (4.5%), pentacosane (3.4%), squalene (3.9%), and trans-β damascenone (2.3%).22
A study was conducted to evaluate the hepatoprotective activity induced in albino rats following administration of S. oleioides aqueous extract. The extract significantly elevated the high-density lipoprotein (HDL) cholesterol level and decreased serum glucose, total cholesterol (TC), triglycerides, low-density lipoprotein (LDL) cholesterol, TC/HDL, and LDL/HDL ratios while the body weight, liver function indices, and selected tumor markers were approximately comparable with the control animals. These findings suggest that S. oleioides extract is very potent in ameliorating liver injury.23 New technological advances have been employed to evaluate the pharmacological properties of natural products using state-of-the-art computer-based approach, the use of which has become inevitable in drug discovery pipelines due to molecular docking and molecular dynamics simulations.24–26 The combination of such molecular modeling techniques alongside investigation of biological activity is an important approach in elucidating novel plant-based drugs.27 Additionally, although many natural products are widely utilized in folk medicine, their toxic potentials have not been reported. Considering these factors, the aim of the current study is to evaluate the protective effects of C. brachycarpa (cabralealactone), S. incanum (salvadorin), and S. oleioides (solasodin) against DNA damage and confirm their hepatoprotective effects.
Materials and methods
All chemicals and reagents used were of analytical grades and were purchased from Sigma Chemical Co. (St Louis, MO, USA). Standardized 90% ethanolic extracts of cabralealactone, 80% ethanolic extracts of both salvadorin and solasodin were purchased from Sigma Chemical Co. Ethical approval was obtained from the Institutional Review Board of The University of Lahore (Approval No: USM/Animal Ethics approval/2009/[45] [140]) and all animal experiments were performed in accordance with the Institutional Guidelines for the Care and Use of Animals for Scientific Purposes. A p-value of <0.05 was considered as statistically significant.
The animals were divided into several groups and received different treatments (Table 1). Subcutaneous injection of CCl4 (1:1 v/v) was administered with saline (1 mL/kg b.wt. per week) to albino rats (Rattus norvegicus) every week for 12 consecutive weeks to induce hepatocellular injury as confirmed by the liver function tests.23 At the end of the experiment, the animals were sacrificed by an intraperitoneal injection of ketamine (Arevipharma GmbH, Radebeul, Germany) (1 mL/g).
Following sacrifice, the liver lobes were separated and collected in glass tubes containing 5% formalin for histopathology examination. For homogenate preparation, freshly dissected liver tissue was homogenized in phosphate buffer (50 mM, pH 7.2), centrifuged at 3,000 rpm for 10 min and the supernatant separated to yield 25% (w/v) homogenate. It was stored at −80°C until biochemical analysis. The estimation of alanine aminotransferase (ALT) was conducted by using a commercially available Bio Merux and Randox kit.28 The levels of α-fetoprotein (AFP), 8-hydroxy-2-deoxyguanosine (8-OHdG), tumor necrosis factor alpha (TNF-α), and isoprostanes-2 alpha (isoP-2α) were determined by using commercially available ELISA kits (Abcam, Cambridge, MA, USA).
Lipid peroxidation in the blood samples was estimated calorimetrically by measuring thiobarbituric acid reactive substances (TBARS).29 Briefly, sample (0.2 mL), sodium dodecyl sulfate (8.1%, 0.2 mL), acetic acid (20%, 1.5 mL), and thiobarbituric acid (0.8%, 1.5 mL) were mixed, followed by centrifugation (3,000 rpm for 10 min). The upper organic layer was aspirated, and the optical density was read at 532 nm using a spectrophotometer (Echelle; LTB Lasertechnik Berlin Gmbh, Berlin, Germany) against the blank. The levels of lipid peroxides were expressed as millimoles of TBARS/g of liver tissue based on a standard curve. DPPH scavenging activity of plant extracts was measured to investigate their antioxidant status.30 Ethanolic DPPH was added to the prepared stirred utions with the plant extracts and the mixture was kept in the dark for 30 min. Absorbance was measured at 517 nm against a blank whereas, ascorbic acid was used as a positive control.
Statistical analysis
The data was subjected to analysis of variance (ANOVA) using COSTAT computer package (Cohort software, Berkeley, CA, USA). The mean values were compared with the least significance difference (at <0.05) test. One-way ANOVA were performed followed by Duncan S’ Multiple Range Test. The values are mean ± SD and p-values (<0.05) were considered as significant.
In silico studies
Dataset collection and molecular docking
The 3 active compounds (cabralealactone, solasodine, and salvadorin) were retrieved from ChemSpider (IDs: 23281688, 391288, and 28282881 respectively) and then converted to 3D Mol2 format using an Open Babel.31,32 Since the crystal structures of Rattus norvegicus cyclooxygenase-2 (COX2) and TNF-α were not present in the protein data bank (PDB) repository, a homology modeling has to be conducted by employing a SwissModel.33 The crystal structures of murine COX2 (PDB ID: 4FM5) and TNF-α (PDB ID: 2TNF) were selected to model the rat proteins.34,35 To generate a dimeric structure of TNF-α, individual modeled chains were later recombined as a single PDB coordinate file.
The predicted models were further evaluated based on Procheck, ERRAT, Z-score, and VERIFY 3D.36–38 The structural alignment of COX-2 and TNF-α homology models on their templates are presented in Figure 1. Subsequently, the AutoDock Vina standard molecular docking procedure was executed as previously documented and analysis of the binding site was performed before docking.39–41 COX2 showed high homology 96% with the template, considering a 100% query cover and an E-value of 0.0, while TNF-α was predicted to have a 95% homology, 87% query cover, and an E-value of 6e−109. Superimposition of the co-crystallized structure with inhibitor led to the identification of binding pocket residues. For each protein, docking grid of 30Å × 30Å × 30Å was made covering the binding site residues. Details of the predocking procedure, including protein and ligand preparation, optimization, and minimization were elaborated elsewhere.42
Molecular dynamics simulations
The molecular dynamics simulation is performed using an AMBER12 package to gain insight into the conformational dynamics of the protein and its binding site during a particular time, and evaluate the validity of the obtained binding mode of ligands from docking calculations.27,43 In addition, molecular dynamics simulations were used to rid unfavorable forces on the protein to reach equilibrium. The ligand geometry obtained from the docking calculations was used as a starting conformation for molecular dynamics study. Partial charges for ligands were calculated and the force field parameters (force constants for angles and bonds) were used based on the proposed parameters of ANTECHAMBER and General Amber Force Field. Amber ff99SB was also used to provide protein parameters for the modeling.41 Hydrogen atoms were added by using the t-leap module in Amber16.
The system was solvated to a 10 Å octagonal box of TIP3P water molecules and neutralized by adding an appropriate number of counter ions. The resultant simulation setup was subjected to the MD simulations in the ligand-bound states using the Particle-Mesh Ewald (PME) molecular docking (MD) module in the AMBER 16 MD package. Before the MD simulation of 20 ns, the system was relaxed by performing an energy minimization to remove the largest strains in the system. Then, after heating the system from 0 to 300 K, isotropic position scaling was used to keep the system at a constant pressure (1 atm) while the Berendsen thermostat method was used to maintain a constant temperature (300 K) during the MD simulation. The PME method was used to investigate the contribution of long-range electrostatic interactions. A 10 Å cut-off was used for both PME and van der Waals interactions. To examine the convergence of the system to equilibrium, the temperature, density, and the total energy was checked to ensure post-MD analysis. MD production run of 20 ns was executed for each complex while the Combined Cluster Analysis (CPPTRAJ) module of Amber 16 was used for post-MD analysis.44
Results
Significant restoration of the investigated biochemical variables in rat groups that received the herb treatments were seen (Table 2). The biochemical parameters of rats that received monotherapy of salvadorin and solasodin extracts (group H) were significantly ameliorated (41.55±4.34 IU/L, 17.71±1.20 ng/L, 18.01±2.38 pg/mL, 15.20±2.58 ng/mL, 15.41±3.03 pg/mL, 3.65±0.53 nmol/mL, and 1.55±0.386 ng/mL) compared with the control group. In addition, the levels of ALT, AFP, 8-OHdG, TNF-α, IsoP-2a, malondialdehyde (MDA), and COX-2 were significantly improved compared with the control. There remained a significant change in the variables of groups treated with plant extracts. Also, combination therapies of these plant extracts caused significant and remarkable restoration of all parameters compared with the standard drug Clavazin. It is plausible that the antioxidant effect of the plant extracts activate the proteolytic enzymes involved in DNA repair that supports its traditional use as anti-cancer agent.
Histopathological findings
The histopathological findings of the liver sections in the control group (Figure 2A) indicated the presence of visible central veins, a terminal hepatic venule, and an interlobular portal vein. The hepatocytes were in contact with the sinusoids, which appear as distensible vascular channels in connection with highly fenestrated endothelial cells surrounded by a population of phagocytic Kuppfer cells. However, CCl4-treated liver sections indicated the presence of prominent necrosis, vacuolization, central vein congestion, and fat vesicles (Figure 2B). The sinusoids, which are in dilated form, are characteristic of widened hepatic capillaries suggestive of hepatic tumors occurring due to sinusoidal infiltration by malignant or benign cells. The prominent architectural alterations in the CCl4-treated toxic liver are karyomegaly (enlargement of the nucleus) and populated activated Kuppfer cells. Examination of the liver sections in animals treated with standardized extracts of cabralealactone (100 mg/kg) revealed that the liver almost regained up to 30% of its normal architecture with reduced intensity of liver fibrosis (Figure 2C). Also, some liver recovery was observed in animals exposed to solasodin (Figure 2D) and salvadorin (Figure 2E) monotherapies as characterized by the lower extent of vacuolization seen.
Nevertheless, compared with the monotherapies, a significantly better protective effect was observed in the liver sections of rats that received combined treatment of standardized extracts of cabralealactone and salvadorin (Figure 2F), cabralealactone and solasodin (Figure 2G), and solasodin and salvadorin (Figure 2H), which were characterized by the presence of normal hepatic cords, absence of fatty infiltration, reduced intensity of sinusoidal dilation, and central vein congestion. However, the best recovery was seen in rat groups exposed to the 3 polyherbal treatment of cabralealactone, salvadorin, and solasodin (Figure 2I), indicating the potentials of these plants as anti-inflammatory agents.
DPPH bscavenging activity
Significant DPPH scavenging activity was observed with salvadorın (94.26%) followed by cabralealactone (89.26%) and solasodın (87.06%) compared with ascorbic acid (91.26%) as standards as in Figure 3.
In silico studies
Molecular docking and MD simulations
The 3 active compounds were used for molecular docking with COX-2 and TNF-α. The chemical 2D structure representations of the active compounds are shown in Figure 4. The binding energy along with the pharmacokinetics properties of each compound is tabulated in Table 3 while the molecular interactions are diagrammatically represented in Figures 5 and 6. Relative to the COX-2 enzyme, salvadorin exhibited a higher binding energy of −9.8 kcal/mol compared with solasodine and cabralealactone, which showed −7.9 and −6.8 kcal/mol, respectively. In relation to TNF-α, salvadorin exhibited a potentially higher binding energy of −11.4 kcal/mol, while the binding energies corresponding to solasodine and cabralealactone were −8.1 and −7.4 kcal/mol, respectively. Moreover, all ligands were docked in the lowest energy conformation inside the binding pocket as shown in Figures 5 and 6.
MD simulation of 20 nanosecond (ns) was performed on each lowest docking energy conformation that was chosen as a starting complex for COX-2 and TNF-α. The root mean square deviation (RMSD) of all backbone atoms of each target protein on the bound ligand was computed as a function of time in 1 ns as displayed in Figure 7. Relative to COX-2-salvadorin complex, the RMSD initially increased until 6 ns and then gradually became stable between the mean value of 1.75Å and 2.51Å until the end of the simulation (red), while the RMSD continued to fluctuate between 1.1Å and 2.2Å in the final 10 ns and continued to increase in both cases where COX-2 complexed with solasodine (yellow) and cabralealactone (green). In relation to the TNF-α-salvadorin complex, the backbone atoms RMSD remained stable throughout the 20 ns simulations between 1.15Å and 1.76Å and oscillated around the mean value of 0.51Å. On the other hand, TNF-α-solasodine showed stability in RMSD with small fluctuations while the backbone atoms of TNF-α complexed with cabralealactone showed higher fluctuations and continued to rise until the end; the ligand remained inside the binding pocket of both proteins throughout 20 ns simulation as was evident from the RMSD of salvadorin (black), solasodine (brown), and cabralealactone (orange).
Discussion
To our knowledge, this is the first study to confirm that cabralealactone, solasodin, and salvadorin have DNA-damage protective effects against CCl4-induced genotoxicity and may have the potential to be investigated as a COX-2 inhibitor and anti-inflammatory agent. The emergence of oxidative stress as a central player in chronic metabolic diseases has opened new avenues for research in the pathology of a wide range of diseases and their treatments to increase the survival rate.45 The mechanism of action of CCl4-induced liver damage is a complex mechanism that indicates its absorption in parenchymal hepatocytes and its metabolism results in the formation of free radicals in the presence of cytochrome P450 -dependent monooxygenases.44 The resultant free radicals form a covalent bond with biological macromolecules preferably with fatty acids and phospholipids of the plasma membrane.46 Once free radicals are produced, lipid peroxidation is initiated, resulting in cellular membrane peroxidation, chloromethylation reaction, polyunsaturated fatty acids (PUFA) deterioration, and ultimately disturbance of the structural integrity. The destruction of phospholipid-associated PUFA leads to physiological alterations in the mitochondria, endoplasmic reticulum, and is mostly prevalent in cellular plasma membrane.47 In addition to this, the metabolism of calcium-associated homeostasis is disturbed, leading to the loss or leakage of cellular calcium. The resulting end products of PUFA form a covalent association with proteins and mimic the catalytic activity of enzymes.48,49
CCl4 disrupts the cellular integrity by activating inflammatory cytokines such as TNF-α, interleukin (IL)-6, ALT, 8-OHdG (the hallmark of DNA-based oxidative damage, MDA, and isoP-2α). A study reported that patients with raised levels of AFP, histologically increased steatosis, and inflammatory infiltration point to DNA-based oxidative stress.48 Another research demonstrated that high levels of oxidative stress contribute to dysfunctioning of adipocytokines, leading to dysregulated production of reactive oxygen species from the adipose tissues that caused increased oxidative damage in the blood. This ultimately affected the normal functioning of the liver, skeletal muscles, and aorta of the patients with chronic liver disease.50,51 The response of cells to the DNA damage is demonstrated via the production of tumor suppressor gene p53 that allows genetic fidelity by inducing cell cycle arrest, ultimately ameliorating the cell damage due to apoptosis as shown in Figure 8.52,53
TNF-α is a cytokine with a molecular mass of 55 kDa, which mediates its cytotoxic activities through the expression of 2 receptors (p55 and p75). The induction of reactive oxygen species production by TNF-α occurs via the endothelial mitochondria and nicotinamide adenine dinucleotide phosphate in the cell membrane. Thus, the cytotoxic activity of TNF-α induced by CCl4 disrupts the hepatocyte mitochondrion function, which stimulates the production of reactive oxygen species in the mitochondria, particularly at the ubiquinone site while TNF-α disrupts the mitochondrial electron transport chain at the complex III region, ultimately resulting in the accumulation of free oxygen radicals inside the mitochondrion. In addition, the increased production of TNF-α also stimulates the nuclear factor-kappa beta activation and upregulated IL-6 expression.
The antioxidant-rich standardized extracts of cabralealactone, salvadorin, and solasodin are useful in repairing the DNA damage and ameliorating the TNF-α-induced cytotoxicity in hepatocytes and malignant transformation of hepatocytes as evidenced by the histological findings. IsoP-2α is recognized as potent indicator of increased oxidative stress in hepatic cirrhosis conditions. The formation of isoprostanes (the isomer of prostaglandin F-2α) is non-enzymatic when free oxygen radicals attack arachidonic acid in the lipid component of the cell membrane. Therefore, isoP-2α is extensively employed as a gold standard index of oxidative stress and DNA damage in liver tissue (Table 4). In our study, high significant correlation exists between some of the aforesaid parameters, which reveals that increased release of ALT from the hepatocytes is responsible for elevated AFP levels in the tissue due to increased production of reactive oxygen species associated with cellular damage, ultimately leading to the DNA damage. Elevated levels of necrosis and regeneration due to CCl4-induced liver injury. Upregulated cellular turnover in the inflammation cycle and oxidative DNA damage results in the appearance of genetic and epigenetic changes, which may be due to the increased expression of cellular oncogenes and growth and angiogenic factors; activation of telomerase activity; and downregulation of tumor suppressor genes.54 Nevertheless, the result should be interpretated with caution since AFP has been reported to yield some false-positive findings (20%) in chronic hepatitis C patients and tends to be higher in the presence of liver cirrhosis (20%–50%).1,55 The current findings corroborated with similar trends seen in the literatures where ALT, MDA, 8-OHdG, AFP, TNF-α, IsoP-2α, and COX-2 levels were elevated. In addition, significant correlation between MDA and 8-OHdG concentrations indicates covalent binding of DNA products with the reactive oxygen species resulting in oxidative DNA damage.
The pyrimidine double bond between C4-C5 is very susceptible to reactive oxygen species, which respond by generating a spectrum of oxidative pyrimidine damage products such as 5-hydroxydeoxyguanosine, thymine glycol, 5-hydroxydeoxycytidine, 8-hydroxydeoxyadenosine, and other purine oxidative products. The production and accumulation of 8-OHdG have been related to carcinogenesis, which causes a decrease in the gene expression and loss of functional integrity of the hepatocytes.56 Variations in the antioxidant capacity of the standardized extracts of cabralealactone, salvadorin, and solasodin conferred significant amelioration on liver regeneration, survival, and decrease in endogenous DNA damage.
The 3 active compounds were finally subjected to in silico studies to determine the binding potential and stability of protein’s all-backbone atoms57,58 via molecular docking, targeting TNF-α and COX2 enzyme followed by 20 ns MD simulation to analyze the atomistic behavior of protein’s backbone with ligand bound inside the binding pocket. Between the 3 compounds, salvadorin yielded the most stable RMSD with minimum fluctuation seen in the final 10–12 ns, which inferred better stability during the formation of the complex (Figure 7). Nevertheless, both in silico and in vivo studies suggested that all 3 active compounds have good and desirable affinity parameters, yielding stable molecular interactions.59–61 Thus, these compounds can be configured as promising agents for the development of COX-2 and TNF-α antagonists.
Limitations
The limitation of the current study is the relatively low number of samples and variables due to financial constraints. Efficacy of the results can be improved by increasing the number of samples and taken into consideration its molecular level. Hence, further studies in these areas would be helpful in understanding and raising new sets of ideas in the field.62–64
Conclusion
It is concluded that cabralealactone, solasodin, and salvadorin showed hepatoprotective and DNA-damage protective effects against CCl4-induced toxicity. The herbs successfully restored the normal architecture of hepatocytes and had the potential to be used as inhibitor to COX-2 and TNF-α. They can combat oxidative stress and liver injuries when used both as mono and combination therapies, with the later showing better findings perhaps due to synergistic effects.
Three-dimensional structural representations of COX-2 and TNF-α, and Ramachandran plot assessment.
Notes: Structural alignment of homology models (in red) of TNF-α (A) and COX-2 (B) with templates (green) are presented in ribbon representation. Ramachandran plots are shown below for TNF-α and COX-2 in the left and right panels, respectively.
Abbreviations: COX-2, cyclooxygenase-2; TNF-α, tumor necrosis factor alpha.
Histopathological findings in rat model. (A) Control liver (×40). (B) CCl4-treated liver (×40). (C) CCl4+cabralealactone-treated liver (×40). (D) CCl4+solasodin treated liver (H&E ×40). (E) CCl4+salvadorin treated liver (×40). (F) CCl4+cabralealactone+solasodin-treated liver (H&E ×40). (G) CCl4+solasodin+salvadorin treated liver (×40). (H) CCl4+cabralealactone+salvadorin treated liver (×40). (I) CCl4+cabralealactone+salvadorin+solasodin treated liver (×40).
Abbreviation: CCl4, carbon tetrachloride.
DPPH scavanging activity (%).
Abbreviation: DPPH, 1,1-diphenyl-2-picrylhydrazyl.
Structure representation of cross-sectional area of COX-2 with docked compounds.
Notes: Cross-sectional area of COX-2 with docked compounds, cabralealactone (pink), solasodine (light brown), and salvadorin (green) in stick representation inside the binding pocket. Molecular interactions in different colors are displayed alongside using a DS visualizer (Accelrys Inc., San Diego, CA, USA).
Abbreviation: COX-2, cyclooxygenase-2.
Molecular surface representation of dimeric TNF-α.
Notes: Molecular surface representation of dimeric TNF-α (in 2 different colors for each chain) with docked compounds, cabralealactone (pink), solasodine (light brown), and salvadorin (green) in stick representation inside the binding pocket. Molecular interactions in different colors are displayed alongside using a DS visualizer (Accelrys Inc., San Diego, CA, USA).
Abbreviation: TNF-α, tumor necrosis factor alpha.
Simulation results of protein-ligand complexes.
Notes: RMSDs during MD simulation of COX-2 and TNF-α of protein-ligand complexes. The RMSDs are plotted as a function of time in 1 ns of the protein backbone of COX-2 and TNF-α with docked ligands during MD simulation at 20 ns. Histopathological findings confirm carbon tetrachloride (CCl4) induces liver injury by binding to the extracellular receptor which results from the toxic metabolites of CCl4 causing centrilobular hepatic necrosis and steatosis that disturb the essential cellular processes. CCl4 is biotransformed by cytochrome P450 system in the hepatic microsomes, producing trichloromethyl free radical (CCl−°3) that reacts with cellular molecules and perturb the cellular processes. In the presence of oxygen, CCl°3 is converted to trichloromethyl peroxy radical (CCl3OO−°), resulting in oxidative stress, lipid peroxidation, mitochondrial dysfunction, and loss of membrane integrity. The described model of liver fibrosis indicates that macrophages progress the differentiation and proper functioning of fibroblasts, and act as potent therapeutic agent. Standardized extract of Cleome brachycarpa (cabralealactone) suppressed the activation of p53, decreased the expression level of anti-apoptotic protein Bax, and increased the expression level of apoptotic Bcl-2, which resulted in reduced cytochrome C release from the mitochondria into the cytoplasm and inactivated caspases 3 and 9. Ultimately, cabralealactone suppressed fibrosis and necrosis of the hepatocytes. Acute injection of CCl4 induces hepatocyte necrosis or apoptosis. Damaged hepatocyte release lipid antigens (represented by HSC) resulting in activation of Nk cells. The activated HSC produce a variety of pro-inflammatory cytokines, that is, TNF-α, NF-κβ, IL-6, and COX-2, contributing to hepatocyte inflammation. Treatment of hepatocyte via oral and intraperitoneal injection of standardized extracts of Solanum incanum (solasodin) and Salvadora oleoides (salvadorin) can attenuate the production of pro-inflammatory cytokines via inhibition of HSC activation and inhibition of the release of cytochrome c oxidase from the mitochondria. Molecular docking studies were carried out to assess the binding affinity of the compounds to proteins, thus mediating tumor growth.
Abbreviations: CCl4, carbon tetrachloride; COX-2, cyclooxygenase-2; IL, interleukin; MD, molecular docking; HSC, hepatic stellate cell; NF-κB, nuclear factor-kappa beta; RMSD, root mean square deviation; TNF-α, tumor necrosis factor alpha.
Molecular mechanism involved in the liver damage induced with the carbon tetrachloride (CCl4).
Abbreviations: CCl4, carbon tetrachloride; COX-2, cyclooxygenase-2; JNK, Jun N-terminal kinase; MAPK, mitogen-activated protein kinase; NF-κβ, nuclear factor-kappa beta; ROS, reactive oxygen species; TNF-α, tumor necrosis factor alpha.
| Groups (n=5) | Treatments |
|---|---|
| A | Control |
| B | CCl4 alone |
| C | CCl4+cabralealactone |
| D | CCl4+solasodın |
| E | CCl4+salvadorın |
| F | CCl4+cabralealactone+solasodın |
| G | CCl4+cabralealactone+salvadorın |
| H | CCl4+solasodın+salvadorın |
| I | CCl4+cabralealactone+solasodın+salvadorın |
| J | CCl4+Clavazin |
Notes: Dose of CCl4 (1 mL/kg b.wt. per week). Dose of cabralealactone, solasodın, salvadorın (100 mg/kg b.wt.). Dose of Clavazin (200 mg/kg b.wt.).
Abbreviation: CCl4, carbon tetrachloride.
Biochemical response of cabralealactone, solasodin, salvadorin following CCl4-induced liver insult in rats
| Groups | Mean ± SD (n=10) | ||||||
|---|---|---|---|---|---|---|---|
| ALT (IU/L) | AFP (ng/L) | 8-OHdG (pg/mL) | TNF-α (ng/mL) | IsoP-2α (pg/mL) | MDA (nmol/mL) | COX-2 (ng/mL) | |
| A | 25.11±4.32f,g | 7.32±3.60f | 2.88±0.33e | 20.52±3.95d | 29.67±2.3e,f | 1.99±0.60g | 0.51±0.0016a |
| B | 179.82±20.03a | 92.26±10.60a | 39.06±10.51a | 163.55±28.36a | 171.96±24.44a | 12.74±1.37a | 2.09±0.035b |
| C | 74.24±11.97b | 40.87±10.19b | 29.15±1.90b | 56.33±6.78b | 103.56±6.25b | 6.35±1.52b,c | 1.60±0.135d |
| D | 63.50±5.58b,c | 37.65±10.82b | 28.89±5.40b | 34.30±7.09c | 87.24±17.81c | 6.77±1.26b | 1.75±0.196d |
| E | 53.07±7.48c,d | 43.15±8.80b | 23.10±3.27c | 22.79±3.11c,d | 56.24±6.40d | 5.41±1.35b–d | 1.55±0.283d |
| F | 42.07±3.25d,e | 27.31±5.58c | 18.02±3.30c | 17.51±1.17d | 32.47±7.03e | 4.05±1.30d–f | 1.67±0.335d |
| G | 40.13±2.75e | 20.494±3.97c,d | 20.72±2.24c | 15.41±2.51d | 19.60±3.72e–g | 5.09±0.85c–e | 1.76±0.462d |
| H | 41.55±4.34d,e | 17.71±1.20d,e | 18.01±2.38c | 15.20±2.58d | 15.41±3.03f,g | 3.65±0.53e,f | 1.55±0.386d |
| I | 20.45±5.05g | 9.36±1.38e,f | 8.69±2.08d | 12.40±1.69d | 9.84±0.97g | 2.63±0.59f,g | 0.62±0.014a |
| J | 35.10±8.12e,f | 10.83±2.56e,f | 11.92±1.95d | 12.69±4.00d | 16.19±2.28f,g | 3.25±1.22f,g | 0.82±0.058c |
| LSD (0.05) | 11.70 | 8.87 | 5.48 | 12.52 | 13.89 | 1.43 | 1.25 |
| p-value | 0.0413 | 0.011 | 0.006 | 0.018 | 0.007 | 0.014 | 0.000 |
Notes: CCl4 @ 1 mL/kg b.wt. per week. Cabralealactone, solasodin, salvadorin @ 100 mg/kg b.wt. Clavazin @ 200 mg/kg b.wt. Superscript letters indicate that values not sharing a common letter differ significantly.
Abbreviations: 8-OHdG, 8-hydroxy-2-deoxyguanosine; ALT, alanine aminotransferase; AFP, α-fetoprotein; CCl4, carbon tetrachloride; COX-2, cyclooxygenase-2; isoP-2α, isoprostanes-2 alpha; LSD, least significance difference; MDA, malondialdehyde; TNF-α, tumor necrosis factor alpha; A, Control; B, CCl4; C, cabralealactone; D, solasodin; E, salvadorin; F, cabralealactone+solasodin; G, cabralealactone+salvadorin; H, solasodin+salvadorin; I, cabralealactone+solasodin+salvadorin; J, Clavazin.
| Properties | Cabralealactone | Solasodine | Salvadorin |
|---|---|---|---|
| Mass | 552.70 | 429.68 | 414.62 |
| Log-p-value | 8.69 | 5.57 | 6.34 |
| H-bond acceptors | 3.00 | 3.00 | 3.00 |
| H-bonds donors | 0.00 | 1.00 | 0.00 |
| Rotatable bonds | 5.00 | 0.00 | 1.00 |
| PSA | 35.53 | 32.70 | 43.37 |
| RO5 violations | 2.00 | 1.00 | 1.00 |
| Atoms | 78.00 | 78.00 | 72.00 |
| Rings | 7.00 | 6.00 | 5.00 |
| Binding affinity of compounds (kcal/mol) with COX-2 | −6.80 | −7.90 | −9.80 |
| Binding affinity of compounds (kcal/mol) with TNF-α | −7.40 | −8.10 | −11.40 |
Abbreviations: COX-2, cyclooxygenase-2; PSA, polar surface area; RO5, Lipinski’s rule of 5; TNF-α, tumor necrosis factor alpha.
Pearson correlation coeeficients of different variables in liver tissues of rats under CCL4 stress receiving cabralealactone, solasodin, salvadorin, and when used in combinations
| Variables | ALT | AFP | 8-OHdG | TNF-α | IsoP-2α | MDA | COX-2 |
|---|---|---|---|---|---|---|---|
| ALT (IU/L) | 1.000 | 0.913** | 0.751** | 0.953** | 0.905** | 0.912** | 0.594** |
| AFP (ng/L) | 1.000 | 0.817** | 0.900** | 0.889** | 0.880** | 0.746** | |
| 8-OHdG (pg/mL) | 1.000 | 0.714** | 0.824** | 0.815** | 0.699** | ||
| TNF-α (ng/mL) | 1.000 | 0.912** | 0.884** | 0.749** | |||
| IsoP-2α (pg/mL) | 1.000 | 0.886** | 0.648** | ||||
| MDA (nmol/mL) | 1.000 | 0.669** | |||||
| COX-2 (ng/mL) | 1.000 |
**Significant at the 0.01 level (two-tailed).
Abbreviations: 8-OHdG, 8-hydroxy-2-deoxyguanosine; ALT, alanine aminotransferase; AFP, α-fetoprotein; COX-2, cyclooxygenase-2; isoP-2α, isoprostanes-2 alpha; MDA, malondialdehyde; TNF-α, tumor necrosis factor alpha.
Footnotes
Disclosure
The authors report no conflicts of interest in this work.
Acknowledgments
The authors are grateful for the valuable contribution of MH Qazi, Director Centre for Research in Molecular Medicine (CRiMM), University of Lahore, Pakistan, for the financial support and for critically reviewing the manuscript.
References
- 1. Hepatocellular carcinoma. A worldwide problem and the major risk factorsDig Dis Sci199136962972[PubMed][Google Scholar]
- 2. Biochemistry of oxidative stressAngew Chem Int Engl19862510581071[Google Scholar]
- 3. Phytochemical and biological ınvestigation of Cleome brachycarpa Vhel. growing in EgyptInt J Pharm Sci Res20145940084014[Google Scholar]
- 4. Problem plants in ngoro ngoro conservation area
Final Report to the NCAA 2002 - 5. Chemical composition and biological activity of the essential oil from Jamaican Cleome serrataNat Prod Commun20127912311232[PubMed][Google Scholar]
- 6. Biochemical aspects of cellular oxidative stressArch Biol Med Exp (Santiago)1988218592[PubMed][Google Scholar]
- 7. Hepatoprotective activity of Cleome viscosa Linn. extract against thioacetamide-induced hepatotoxicity in ratsNat Prod Res2009231412891297[PubMed][Google Scholar]
- 8. In vivo investigation of hepatoprotective activity of Cleome viscosa L. in albino ratsDer Pharmacia Lettre201683308313[Google Scholar]
- 9. In vitro antibacterial activity and total phenolic and flavonoid content of Cleome brachycarpa Vahl ex DC methanol extract and its subfractionsIJGHC20143312381245[Google Scholar]
- 10. Evaluation of antioxidant activity of Cleome brachycarpa Vahl ex DC, an under-exploited desert plant of United Arab EmiratesPharmacology20123167173[Google Scholar]
- 11. Phytochemical and biological investigation of Cleome brachycarpa Vahl. growing in EgyptInt J Pharm Sci Res20159024[Google Scholar]
- 12. Chemistry-based studies on oxidative DNA damage: formation, repair, and mutagenesisFree Radic Biol Med200233450456[PubMed][Google Scholar]
- 13. Evaluation of anticancer activity of Cleome gynandra on Ehrlich’s Ascites Carcinoma treated miceJ Ethnopharmacol20101291131134[PubMed][Google Scholar]
- 14. Ripe fruit of Solanum nigrum L. inhibits cell growth and induces apoptosis in MCF-7 cellsFood Chem Toxicol2003411014211428[PubMed][Google Scholar]
- 15. Isolation of steroidal glycoalkaloids from Solanum incanum by two counters current chromatographic methodsPhytochemistry1991302685687[PubMed][Google Scholar]
- 16. Solamargine induces apoptosis and sensitizes breast cancer cells to cisplatinFood Chem Toxicol2007451121552164[PubMed][Google Scholar]
- 17. Oxidative damage to the lipids and proteins of the lungs, testis and kidney of rats during carbon tetrachloride intoxicationClin Chem Acta1999289177179[Google Scholar]
- 18. Spectrophotometric method for estimation of alkaloids precipitable with Dragendorffs reagent in plant materialJ AOAC Int200386611241127[PubMed][Google Scholar]
- 19. The effect of kefir on the activities of GSH-Px, GST, CAT, GSH and LPO levels in carbon tetrachloride-induced mice tissuesJ Vet Med B Infect Dis Vet Public Health200350412416[PubMed][Google Scholar]
- 20. Hepatoprotective activity of Solanum incanum (thorn apple) against carbon tetrachloride intoxication in mouse liverJ Ethnopharmacol19954726974[PubMed][Google Scholar]
- 21. Reduction of carbon tetrachloride-induced hepatotoxic effect by oral administration of betaine in male wistar rats: a morphometric histolgic studyVet Pathol2000373231238[PubMed][Google Scholar]
- 22. Pharmacognostical and phytochemical investigation of Salvadora oleoides decne. stemIJPPS20135Suppl 2128130[Google Scholar]
- 23. Antidiabetic and antihyperlipidemic effects of ethanolic extract of Salvadora persica L. on Alloxan-induced diabetic ratsPelagia Res Lib201343178182[Google Scholar]
- 24. Molecular docking and in silico ADMET studies of silibinin and glycyrrhetic acid anti-inflammatory activityTrop J Pharm Res20171616774[Google Scholar]
- 25. Medicinal plant phytochemicals and their inhibitory activities against pancreatic lipase: molecular docking combined with molecular dynamics simulation approachNat Prod Res201717[Google Scholar]
- 26. Structure-based virtual screening and molecular docking for the identification of potential multi-targeted inhibitors against breast cancerBreast Cancer (Dove Med Press)20179447459[PubMed][Google Scholar]
- 27. Integrated computational approach for virtual hit ıdentification against ebola viral proteins VP35 and VP40Int J Mol Sci201617111748[Google Scholar]
- 28. Bioinformatics opportunities for identification and study of medicinal plantsBrief Bioinform2012142238250[PubMed][Google Scholar]
- 29. Protective effect of n-acetyl cysteine in carbon tetrachloride-induced hepatotoxicity in ratsIran J Pharmacol Ther20054118123[Google Scholar]
- 30. A colorimetric method for the determination of serum glutamic oxalacetic and glutamic pyruvic transaminassAm J Clin Pathol19572815683[PubMed][Google Scholar]
- 31. Assay for lipid peroxidation in animal tissues by thiobarbituric acid reactionAnal Biochem197995351358[PubMed][Google Scholar]
- 32. Antioxidative properties of xanthan on the anti-oxidation of soybean oil in cyclodextrin emulsionJ Agric Food Chem1992401992945948[Google Scholar]
- 33. ChemSpider: an online chemical information resourceJ Chem Educ2010871111231124[Google Scholar]
- 34. Open babel: an open chemical toolboxJ Cheminform20113133[PubMed][Google Scholar]
- 35. SWISS-MODEL: an automated protein homology-modelling serverNucleic Acids Res2003311333813385[PubMed][Google Scholar]
- 36. Substrate-selective inhibition of cyclooxygenase-2: development and evaluation of achiral profen probesACS Med Chem Lett201239759763[PubMed][Google Scholar]
- 37. The structure of mouse tumour-necrosis factor at 1.4 Å resolution: towards modulation of its selectivity and trimerizationActa Crystallogr D Biol Crystallogr1999554772778[PubMed][Google Scholar]
- 38. PROCHECK: a program to check the stereochemical quality of protein structuresJ Appl Crystallogr1993262283291[Google Scholar]
- 39. ProSA-web: interactive web service for the recognition of errors in three-dimensional structures of proteinsNucleic Acids Res2007352407410[Google Scholar]
- 40. VERIFY3D: assessment of protein models with three-dimensional profilesMethods Enzymol1997277396404[PubMed][Google Scholar]
- 41. Glycyrrhetinic acid and E. resveratroloside act as potential plant derived compounds against dopamine receptor D3 for Parkinson’s disease: a pharmacoinformatics studyDrug Des Devel Ther20159187[Google Scholar]
- 42. Pharmacoinformatics approach for investigation of alternative potential hepatitis C virus nonstructural protein 5B inhibitorsDrug Des Devel Ther201591825[Google Scholar]
- 43. AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreadingJ Comput Chem2010312455461[PubMed][Google Scholar]
- 44. AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of moleculesComput Phys Commun19959113141[Google Scholar]
- 45. Towards peptide vaccines against Zika virus: ımmunoinformatics combined with molecular dynamics simulations to predict antigenic epitopes of Zika viral proteinsSci Rep2016637313733[Google Scholar]
- 46. Structural survey of zinc-containing proteins and development of the zinc AMBER force field (ZAFF)J Chem Theory Comput20106929352947[PubMed][Google Scholar]
- 47. PTRAJ and CPPTRAJ: software for processing and analysis of molecular dynamics trajectory dataJ Chem Theory Comput20139730843095[PubMed][Google Scholar]
- 48. Carbon tetrachloride-induced nephrotoxicity and protective effect of betaine in Sprage-Dawley ratsUrology2003622353356[PubMed][Google Scholar]
- 49. Octacosanol attenuates disrupted hepatic reactive oxygen species metabolism associated with acute liver injury progression in rats intoxicated with carbon tetrachlorideJ Clin Biochem Nutr2008422118125[PubMed][Google Scholar]
- 50. The prevention of carbon tetrachloride-induced hepatotoxicity in mice by α-hederin: inhibition of cytochrome P450 2E1 expressionBiochem Mol Biol Int199845163170[PubMed][Google Scholar]
- 51. Improved preservation of erythrocyte antioxidant capacity with ascorbic acid reperfusion after ischemia: a comparative study in rat hind limb modelInt Journal Altern Med200521[Google Scholar]
- 52. Erythrocyte antioxidant enzymes in patients with cataractAnn Clin Lab Sci200636201204[PubMed][Google Scholar]
- 53. Hepatocellular carcinoma: epidemiology and molecular carcinogenesisGastroenterology200713225572576[PubMed][Google Scholar]
- 54. Utility of lentil lectin affinity of alpha-fetoprotein in the diagnosis of hepatocellular carcinomaJ Hepatol199625166171[PubMed][Google Scholar]
- 55. Hepatoprotective and alterative properties of kombucha tea against carbon tetrachloride-induced toxicityJ Microbiol Biotechnol2009194397402[PubMed][Google Scholar]
- 56. DNA damage-induced phosphorylation of p53 alleviates inhibition by MDM2J Cell199791325334[Google Scholar]
- 57. High α-fetoprotein level in HCV-related nodular liver cell dysplasiaAm J Gastroenterol20119422962297[Google Scholar]
- 58. Clinical significance of elevated AFP in patients with chronic hepatitis C, but not hepatocellular carcinomaAm J Gastroenterol200499860865[PubMed][Google Scholar]
- 59. A novel mechanism of cyclooxygenase-2 inhibition involving interactions with Ser-530 and Tyr-385J Biol Chem2003278464576345769[PubMed][Google Scholar]
- 60. The structure of ibuprofen bound to cyclooxygenase-2J Struct Biol201518916266[PubMed][Google Scholar]
- 61. COX-2 inhibition potentiates antiangiogenic cancer therapy and prevents metastasis in preclinical modelsSci Transl Med20146242242[Google Scholar]
- 62. Antibacterial properties of some essential oils against fine enterobacteriaceae speciesAL-Bututh20071121119[Google Scholar]
- 63. Anxiolytic and CNS depressant Effects of ethanolic extract of Cleome brachycarpa revealed after Neuropharma-cological ScreeningWJPS2014217605610[Google Scholar]
- 64. Chemical constituents of essential oil of Salvadora oleoidesJ Pharmacy Res2012523662367[Google Scholar]