Computational evaluation of bioactive compounds from Vitis vinifera as a novel β-catenin inhibitor for cancer treatment

β-catenin is an important unit of the Wnt/β-catenin signaling pathway, a conserved process involving several physiological activities, encompassing differentiation and cell proliferation, etc. The dysfunction or mutation in β-catenin causing the initiation and advancement of various neoplasm types, including colorectal cancer, breast cancer, etc., has been reported. Therefore, β-catenin is a therapeutic target. Hence, designing new inhibitors targeted against β-catenin will prevent cancerous cells’ involvement and eliminate the diseases. Studies showed that Vitis vinifera, a well-known grape species, contains different phytochemical substances, including aromatic acids, flavonoids, phenolic compounds, proanthocyanins, etc. V.vinifera exerts different anticancer properties such as apoptosis, cell proliferation, cell cycle arrest, and inhibition in cancerous cells. Structural bioinformatics methods, including molecular docking, molecular mechanics generalized Born surface area (MM/GBSA), absorption, distribution, metabolism, excretion studies (ADMET), and pharmacophore modeling approach, were used to determine the potential β-catenin inhibitors from V.vinifera bioactive compounds. Cis-astringin, rutin, caftaric acid, trans-caftaric acid, procyanidin B3, cis-Miyabenol C, and ampelopsin H are shown to be suitable inhibitors against β-catenin due to their binding affinity and interaction with the amino acids residues at the binding sites of β-catenin compared to Food and Drug Administration (FDA) approved drugs leucovorin Calcium and Xeloda prescribed to cure colorectal cancer. This study suggests that V. vinifera could be a good plant source for compounds that might treat cancer by inhibiting the Wnt/β-catenin signaling pathway.


Background
The Wnt/β-catenin pathway, otherwise named the canonical Wnt pathway, is a preserved mechanism that participates in different physiological activities, including apoptosis, cell proliferation, tissue homeostasis, and differentiation (Choi et al. 2020;Salik et al. 2020;Soleas et al. 2020). The dysfunction of β-catenin, a significant part of the Wnt signaling pathway, has been connected with tumors' initiation and progression. High level of β-catenin occurs when there is overexpression of the canonical WNT pathway. Hence, mutations in the Wnt/β-catenin signaling pathway have been involved in different cancer types, particularly colorectal cancer, osteoporosis, and other human diseases (Clevers & Nusse, 2012;Moon et al. 2002Moon et al. , 2004Nusse, 2005). The high mortality rate due to cancer can be prevented if this pathway is inhibited.
The Wnt signaling pathway can activate both the β-catenin-independent non-canonical pathway and the β-catenin-dependent canonical pathway. When Wnt signaling is absent, the degradation complex, comprising Axin scaffolding protein, tumor suppressor adenomatous polyposis coli (APC), and the (Glycogen synthase kinase 3) GSK 3β, degrades β-catenin, and this prevents the buildup of β-catenin. When Wnt signaling is present, the WNT receptor combines with low-density lipoproteinreceptor related protein 5/6 (LRP5/6) and Frizzled (FZD). This causes the interruption of the degrading complex, which induces the buildup of β-catenin and then moves to the nucleus and combines with the T cell factor/lymphoid enhancer-binding factor (TCF/LEF). Thus, promoting different β-catenin gene activation as shown in Fig 1. Mutations and dysfunction in Wnt/β-catenin signaling component genes such as tumor suppressor gene, β-catenin genes, and adenomatous polyposis coli (APC) are implicated in the constitutive activation of Wnt signaling in more than 90% of colorectal cancer cases (Novellasdemunt et al. 20152015). Hence, there is a need to introduce therapy targeted against the Wnt/β-catenin signaling pathway, which prevents cancerous cell involvement and gets rid of the diseases. Therefore, careful study of the genes involved or unregulated in colorectal cancer will give a clue for the generation of new sustainable therapeutic interventions. Phytochemicals from natural sources possess promising chemoprotective properties, which might serve as reliable and alternative cancer treatments. This study used bioactive compounds from Vitis vinifera to design new inhibitors targeted against β-catenin.
Vitis vinifera is a well-known grape species of the family Vitaceae. The different parts of Vitis vinifera contain phytochemical substances, including aromatic acids, flavonoids, phenolic compounds, proanthocyanins, stilbenoids, etc. (Filocamo et al. 2015;Goufo et al. 2020;Radulescu et al. 2020). Research showed that Vitis vinifera's bioactive constituents have anticancer, antioxidant, antidiabetic, antiviral, and anti-inflammatory potentials   (Chaudhary et al. 2021). Vitis vinifera exerts different anticancer properties such as inhibition, apoptosis, cell proliferation, and cell cycle arrest in cancerous cells (Nandakumar et al. 2008). This study employed a comprehensive computational approach to determine the therapeutic properties of different bioactive compounds from Vitis vinifera. The bioactive compounds were used to impede the Wnt/βcatenin pathway. The study includes molecular docking, ADMET screening, pharmacokinetics, and 3D pharmacophore modeling.

Ligand library generation and preparation
Secondary metabolites from V. vinifera that have been described were extracted in standard format (SDF) from the online database, Pubchem (Kim et al. 2016). Using the ligPrep tool (Release Schrödinger, 2017), the mined structures were converted into a three-dimensional structure by ionizing at pH (7.2 ± 0.2) and removing salt with Epik (Shelley et al. 2007;Schrodinger, 2021). The OPLS3 force field (Harder et al. 2016) was used for ionization and tautomeric state creation.

Preparation of target
X-ray crystalline structure of β-catenin coupled with an inhibitor (PDB ID: 1JDH) (Asthana et al. 2014) was recovered from Protein Data Bank. The protein was processed with the protein preparation wizard tool of Maestro, Schrodinger Suite. The protein was refined by optimizing the H-bond assignment and minimized using the OPLS3 force field.

Receptor grid generation
The receptor grid depicts the sector where the ligand and protein interact. The Receptor Grid Generation tool was used to create the processed protein grid on the binding site (Glide Grid). Selecting the co-crystallized ligand at the active site of 1JDH revealed the binding location. The ligand on the protein's crystal structure revealed information about the active site. A cubic grid box including all amino acid residues at the active site was automatically generated with the coordinate X = − 0.360, Y = 14.401, and Z = 47.310, respectively.

Docking
Docking was done on maestro 11.1 (Schrödinger Release 2017) with the Glide tool. Using extra precision (XP) and standard precision (SP) docking algorithm, the crystal structure of β-catenin was utilized to virtually screen the prepared compounds to predict compounds with the lowest docking score. The docking research was carried out with the protein expressed as a rigid body and the ligand's rotatable bonds set to be free.
To create comparison research, known cancer standards were docked using the same approach.

ADMET/Tox screening
To establish the pharmacokinetic profile, drug-likeness, and toxicity of the hit compounds, the swissADME (http:// www. swiss adme. ch) and Pro-Tox II online servers (https:// tox-new. chari te. de/ protox II) online servers were used.

MM/GBSA
The molecular mechanics/generalized Born surface area (MM/GBSA) continuum solvent model was employed to discover the docked protein-ligand complex binding free energy.
Rotamer search techniques from prime were used with the VSGB solvent approach and OPLS3 force field to complete this project.

Pharmacophore modeling
The lead compounds' receptor-ligand complexes were examined. A hypothesis (E-pharmacophore) was generated using the Schrodinger suite's phase interface to highlight the significant properties that actively contribute to the lead ligands' characteristic binding to the active sites of the target proteins.

Results
This study features a computational approach to screen phytochemicals of V. vinifera against β-catenin, a therapeutic target protein that takes part in carcinogenesis through the wnt pathway by employing molecular docking, pharmacokinetics screening, MM/GBSA, and pharmacophore modeling approach.
The molecular docking and binding energy study (Figs. 2 and 3) predicted seven compounds as lead compounds that show a better inhibitory potential when compared with the standard compounds, leucovorin Calcium and Xeloda.
A pharmacophore hypothesis which analyzes the characteristic binding of the top-scoring ligands to the target was analyzed and represented in Fig. 4 The lead compounds were subjected to ADME/Tox screening to study their physicochemical properties, pharmacokinetics profile, drug-likeness and toxicity. The results are recorded in Tables 1, 2 and 3.

Discussion
The molecular docking result identified seven compounds from V. vinifera that showed better inhibitory potential against the target enzyme by having better docking scores when compared with the standard compounds such as leucovorin Calcium and Xeloda currently used to manage cancer. A more negative docking score indicates more inhibitory potential. The docking scores, as represented in Fig. 2, show that cis-astringin, rutin, caftaric acid, trans-caftaric acid, procyanidin B3, cis-Miyabenol C, and ampelopsin H show better inhibitory potential by having docking scores of − 7.925, − 6.460, − 6.379, − 6.379, − 6.170, − 5.734 and − 5.548 kcal/ mol, respectively, compared to leucovorin Calcium and Xeloda which are standard compounds used to treat cancer with the docking score of − 5.490 and − 3.285 kcal/ mol, respectively. The differences in the lead compounds and standards docking scores are pictorially represented in Fig. 2.
Structural-based drug design is predominantly based on the interaction of protein-ligand, as significant inhibition is primarily determined by the interaction of the ligand with the combination of amino acid residues at the active site of the target enzyme. The interaction which contributes significantly to the inhibition of β-catenin in this study is shown in Fig. 3.
The 2D interaction, as shown in Fig. 3 The MM-GBSA technique is a significantly more accurate way of estimating protein-ligand complexes' binding free energies (dG) (Bandyopadhyay et al. 2021). It is one of the promising approaches used to increase virtual screening results. A negative dG value suggests that the complexes generated were stable in the target's binding pocket (Bathula et al. 2021). All lead compounds show a negative dG value, as shown in Fig. 2. The binding free energy of the docked complexes was − 54.285, − 42.282, − 31.078, − 31.078, − 43.693, − 42.348, and − 53.368 kcal/mol for cis-astringin, rutin, caftaric acid, trans-caftaric acid, procyanidin B3, cis-Miyabenol C, and ampelopsin H, respectively. This implies that all the hit compounds are more stable in the target's binding pocket than the standard compounds leucovorin Calcium and Xeloda, which have − 28.175, and − 24.100 kcal/mol as their free binding energy.
A pharmacophore hypothesis founded on the complex formed was developed to analyze the characteristic binding of the top-scoring ligands to the target. PHASE graphical user interface in Schrodinger's suite puts forward essential information on the molecular orientation  Interaction between the lead compounds and the amino acid residues at the active site of the target compared with standards. A-2D interaction between cis-astringin and the amino acid residues at the active site of the target. B-2D interaction between Rutin and the amino acid residues at the active site of the target. C-2D interaction between caftaric acid and the amino acid residues at the active site of the target. D-2D interaction between trans-caftaric acid and the amino acid residues at the active site of the target. E-2D interaction between procyanidin B3 and the amino acid residues at the active site of the target. F-2D interaction between cis-Miyabenol C and the amino acid residues at the active site of the target. G-2D interaction between ampelopsin H and the amino acid residues at the active site of the target. H-2D interaction between leucovorin Calcium and the amino acid residues at the active site of the target. I-2D interaction between Xeloda and the amino acid residues at the active site of the target  of vital functional groups predominantly involved in the characteristic binding of high-affinity ligands to the protein target (Dixon et al. 2006).
The generated E-pharmacophore hypothesis (Fig. 4) contains H-bond donor, H-bond acceptor interactions, and aromatic rings. This hypothesis can be used to form the basic structural foundation that will be general to all potential inhibitors of the protein target and can also be employed to create a basic skeleton of compounds with a specific angular distance that will bond firmly to the ligand-binding site of the targets.
The absorption, distribution, metabolism, and excretion (ADME) analysis has proven to be an essential component of the drug discovery process. It is a highly economical method for in silico prediction of drug action within a biological system (Kenakin et al. 2016). In this study, we use the SwissADME and Protox-II servers to predict the drug-like and toxicity profile of the selected lead and standard compounds.
The physicochemical properties are established on the Lipinski rule, which considers the five (5) main criteria to determine if a drug is orally active (Lipinski et al. 1997). The results are shown in Table 1. Cis-astringin showed very good physicochemical properties, low molecular weight, and a favorable bioavailability score and solubility. This makes it the top drug-like compound among the test compounds based on physicochemical profiling. While Calcium Xeloda showed a good physicochemical profile with no violation of the Lipinski's rule, leucovorin has a very poor physicochemical profile, and the molecular weight exceeds the Lipinski's standard (500 g/mol), the bioavailability score is low, and compared to Cisastringin, it is a poor compound.
The pharmacokinetic profile of the compounds is shown in Table 2. Ampelopsin H, rutin, and the standard compounds were P-glycoprotein substrates. P-glycoprotein is one of the ATP binding cassette (ABC) proteins involved in discharging molecules from the cell and preventing compounds from bioaccumulating and eliciting their response. (Lin et al. 2003) The low gastrointestinal absorption rate of the test compounds and leucovorin indicates their poor lipophilic nature. However, Calcium Xeloda showed a high rate of gastrointestinal absorption. The test and standard compounds do not permeate the blood-brain barrier. Cytochrome P450 (CYP) is a family of enzymes that catalyze the phase 1 metabolism of xenobiotics at large. Any compound that inhibits selected isoforms would induce a drug-drug interaction. (Esteves Table 2 SwissADME pharmacokinetic profile of the lead and standard compounds . Except for procyanidin B3, which was predicted to inhibit the isoform CYP34A, all of the other test compounds and the standard compounds were expected to be non-inhibitors of these enzymes and these cannot induce a drug-drug interaction. The Pro-tox II online server estimated the toxicity profile of compounds, and it is shown in Table 3. The least toxic compound is rutin. It has the highest predicted LD50 value of 5000 mg/kg and ranks 5 in the toxicity class. (It ranges from 1 to 6, i.e., most toxic to least toxic.) It is predicted to have no carcinogenic or hepatotoxic effect. Leucovorin has an expected LD50 value of 135 mg/kg and ranks 3 in toxicity, making it the most toxic compound. Although Calcium Xeloda is not predicted to be as toxic as leucovorin, it is a hepatotoxic compound.

Conclusions
This work used in silico analysis to predict seven compounds from V. vinifera phytoconstituents: cis-astringin, rutin, caftaric acid, trans-caftaric acid, procyanidin B3, cis-Miyabenol C, and ampelopsin H as lead compounds that can bind to β-catenin more robustly than the conventional compounds. When these lead compounds are subjected to ADME/Tox screening, it is discovered that they are somewhat drug-like, with cis-astringin being the most drug-like compound compared to the co-crystallized compound.
This shows that these compounds could be better druglike molecules, and V. vinifera could be a good plant source for a drug-like compound that could cure tumors by suppressing the Wnt//β-catenin pathway's target enzyme, β-catenin.
Lead compound optimization and additional in vitro and in vitro analyses are recommended to validate this experiment.