Druglikeness Analysis of Curcumin
The SMILES format of curcumin, (COC1=C(C=CC(=C1)C=CC(=O)CC(=O)C=CC2=CC(=C(C=C2)O)OC)O), was imported into swissame website to obtain relevant parameters.According to the Lipinski's rule of five, a drug-like compound should have a molecular weight of less than 500 g/mol, a polar surface area (PSA) of less than or equal to 140 Â, a computed octanol/water partition coefficient (XLogP3-AA) of less than 5, less than 10 rotatable bonds (RB), no more than 10 hydrogen bond acceptor (HBA), and no more than 5 hydrogen bond donors (HBD) (Chen et al. 2019). It can be seen from the obtained parameters that the properties of curcumin comply with the RO5, indicating that it has good drug-like properties (Table 1).
Composition-Targets Network
The related targets of curcumin were searched, 104 targets were obtained after removing the duplication, and 1911 CRC-related targets were obtained after removing the duplication. Next, 30 common targets were screened out, that were considered potential targets of curcumin in the treatment of CRC (Figure 3A). The composition-targets network was constructed by Cytoscape v3.7.2 (Figure 3B).
Construction of Protein-Protein Interaction Network (PPI)
We uploaded 30 common targets to the STRING database to determine their functional relationships and interactions. Then the protein interactions with the default confidence level of 0.4 were imported into Cytoscape v3.7.2 to generate a protein-protein interaction (PPI) network, which consisted of 26 nodes and 90 edges, as shown in Figure 4.
To identify the pivot nodes and essential proteins in the PPI network, the topology parameters of the node degree were calculated by the network analyzer, and the three centralities (betweenness, closeness and subgraph) were determined through the CytoNC as shown in Table 2.
GO and KEGG Pathway Enrichment Analysis
The GO and KEGG enrichment analysis were performed via David platform. The GO enrichment analysis is composed of biological process (BP), cellular component (CC) and molecular function (MF). A total of 140 items, BP: 93, CC:14, MF: 33, were obtained by GO functional analysis. The results showed that the effects of curcumin were related to protein kinase activity, ATP binding, negative regulation of apoptotic process and protein serine/threonine kinase activity,et al., as shown in Figure 5A.
In the enrichment analysis of KEGG pathway, 61 enrichment results were obtained. A total of 20 typical pathways were selected to make the visualizedbubble diagram after excluding irrelevant pathways (Figure 5B). The results showed that these pathways were mainly related to pathways in cancer, PI3K and Akt signaling pathway, FOXO signaling pathway, et al. Six targets (AKT1, RAF1, BRAF, EGFR, IKBKB, and STAT3) in the first 20 pathways participated in a high frequency (≥9 times), indicating that they played important roles in CRC. Ten representative signaling pathways are selected to construct a "pathways-targets" network, as shown in Figure 5C.
Molecular Docking
Curcumin is docked with three important targets AKT1, STAT3 and EGFR. These targets are selected not only because they are the key nodes of PPI network, but also they play important roles in KEGG enrichment pathways. The binding energies of AKT1, STAT3 and EGFR with curcumin were -9.9 kcal/mol, -8.7 kcal/mol, -8.5 kcal/mol, respectively. The binding energies of matrine to AKT1,STAT3 and EGFR were -7.8 kcal/mol, -8.7 kcal/mol, -7.6 kcal/mol, respectively(Table 3). It can be seen that curcumin has a strong binding force with key targets. The binding of curcumin with AKT1 is mainly through the hydrogen bonding with amino acid residues ASN53 and GLN79, hydrophobic interaction with TRP80, and π bonding with LEU210, LEU264, LYS268, VAL270, ILE84. The binding of curcumin with EGFR is mainly through the hydrophobic interaction with amino acid residues VAL762, PHE856, ALA743,MET790, CYS775, and π-bond interaction with LEU844. The binding of curcumin with STAT3 is mainly through the hydrogen bonding of amino acid residues ASP1021, ASN1008, ARG1007, GLU957, GLY962, hydrophobic interaction with LEU881, VAL889, ALA906 and π bonding with LEU1010 (Figure 6).
Experimental Verification
Curcumin Promoted Apoptosis of CRC Cells
The effect of curcumin on the apoptosis of HCT116 cells was evaluated throughflow cytometric analysis. After treatment with 0, 12.5, 25 and 50μmol/L of curcumin for 24 h, CRC cells were stained with Annexin V-FITC and PI to determine the degree of apoptosis. The results showed that the percentage of apoptotic cells increased significantly in a dose-dependent manner after curcumin treatment of HCT116 cells, indicating that curcumin induced apoptosis of HCT116 cells (Figure 7).
Validation of Targets
We further verified these targets were involved in curcumin-induced apoptosis in CRC cells by Western blotting. As shown in Figure 8, AKT1 protein levels decreased significantly in a dose-dependent manner. The difference was statistically significant (p<0.05). The results showed that AKT1 was an important target of curcumin in the treatment of CRC.