Next Article in Journal
Comparison of the Differences between Two-Photon Excitation, Upconversion, and Conventional Photodynamic Therapy on Cancers in In Vitro and In Vivo Studies
Next Article in Special Issue
Unraveling the Impact of Six Pentacyclic Triterpenes Regulating Metabolic Pathways on Lung Carcinoma Cells
Previous Article in Journal
Small Molecule Drugs Targeting Viral Polymerases
Previous Article in Special Issue
The Therapeutic Potential of Essential Oils in Managing Inflammatory Skin Conditions: A Scoping Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Profibrotic Inflammatory Cytokines and Growth Factors Are Predicted as the Key Targets of Uncaria gambir (Hunter) Roxb. in Keloids: An Epistatic and Molecular Simulation Approach

by
Sri Suciati Ningsih
1,2,
Fadilah Fadilah
1,3,
Sri Widia A. Jusman
4,5,
Rahimi Syaidah
1,6,* and
Takashi Yashiro
7
1
Doctoral Program in Biomedical Sciences, Faculty of Medicine, Universitas Indonesia, Jakarta 10430, Indonesia
2
Faculty of Medicine, Universitas Muhammadiyah Prof. Dr. Hamka, Jakarta 12130, Indonesia
3
Department of Medical Chemistry, Faculty of Medicine, Universitas Indonesia, Jakarta 10430, Indonesia
4
Department of Biochemistry, Faculty of Medicine, Universitas Indonesia, Jakarta 10430, Indonesia
5
Center of Hypoxia and Oxidative Stress Studies, Faculty of Medicine, Universitas Indonesia, Jakarta 10430, Indonesia
6
Department of Histology, Faculty of Medicine, Universitas Indonesia, Jakarta 10430, Indonesia
7
Department of Histology, Jichi Medical University School of Medicine, Tochigi 329-0498, Japan
*
Author to whom correspondence should be addressed.
Pharmaceuticals 2024, 17(6), 662; https://doi.org/10.3390/ph17060662
Submission received: 17 April 2024 / Revised: 3 May 2024 / Accepted: 8 May 2024 / Published: 21 May 2024

Abstract

:
Keloid is characterized as the fibrotic tissue resulting from the increase of fibroblast activity. Uncaria gambir (Hunter) Roxb. possesses bioactive compounds that have potential as antifibrotic agents, while the mechanism of action in keloid has not yet been elucidated. The aim of this study was to investigate the interaction of gambir bioactive compounds with keloid target proteins using an epistatic and molecular simulation approach. The known bioactive compounds of gambir targets and keloid-related protein targets were screened using databases. The network was constructed and analyzed to obtain the core protein targets. The targets were enriched to describe the Gene Ontology (GO) and pathway related to the proteins. Eleven targets were defined as the main targets of gambir bioactive compounds related to keloid disease. Gambiriin C, Isogambirine, and Procyanidin B1 were identified as the most promising compounds with the highest binding energy to transforming growth factor beta 1 (TGFβ1), AKT serine/threonine kinase 1 (AKT1), and matrix metallopeptidase 1 (MMP1) as the target proteins. GO enrichment and pathway analysis found that gambir bioactive compounds may act on keloid-related target proteins to regulate cell proliferation, migration, transcription, and signal transduction activity via profibrotic cytokine and growth factor signaling pathways. This study provides a reference for potential targets, compounds, and pathways to explain the mechanism of gambir against keloid.

Graphical Abstract

1. Introduction

Keloids display exophytic overgrowth beyond the site of injury. The latest demographic study of keloid incidence reported that the incidence of keloids is higher in pigmented-skin individuals, particularly in Africa (8.5–16%). Even though this disease may not directly result in fatality, the manifestation may cause cosmetic, physical, and physiological effects that consequently affect the patient’s quality of life [1,2]. The pathomechanism of keloid has not been fully elucidated. Abnormality of the wound-healing process is addressed as the critical point of keloid formation. Robust inflammatory stimulation and over-sensitized fibroblasts driven by cytokines and growth factors during the wound-healing process contribute to excessive extracellular matrix (ECM) deposits in keloids. Over-sensitized fibroblasts are activated into myofibroblasts, which tend to have a higher level of proliferation, migration, and extracellular matrix synthesis capacity [3,4,5].
Several existing treatments for keloid, including surgery, steroid therapy, stem cell, radiation, and targeted therapy, are intended to inhibit the proliferation and activation of keloid fibroblasts. Unfortunately, none of these are considered a universal gold standard for keloid treatment due to the high recurrence rate, emergence of side effects, and requiring high cost and further study to reveal the effects in a clinical approach [2,6,7]. Thus, herb-based therapy is proposed as an option for an alternative way to treat keloid disease. Traditional herbal medicine possesses the characteristics of being multicomponent, multifunction, and multitarget; is of low toxicity; and has limited side effects, in general [8].
Gambir is an herbal medicine originating from Indonesia that has long been used by local residents as a traditional remedy. A majority of the phytochemical substances derived from gambir are categorized into the flavonoid, phenolic, and alkaloid groups [9]. Recent studies proved that the bioactive compounds of gambir have antibacterial, anticancer, anti-inflammation, as well as antifibrotic effects [10,11,12,13]. Desdiani et al. [14,15] found using in vivo study that gambir has antifibrotic effect against lung fibrosis through inhibition of nuclear factor kappa-B (NF-κB) and transforming growth factor receptor (TGF-βR). In line with these findings, Sriningsih [16] identified that gambir extract could inhibit the expression of tissue inhibitor matrix metalloproteinase (TIMP) and collagen-1 in rat liver fibrosis tissue. Keloid is defined as tissue fibrosis due to its effusive accumulation of ECM, particularly collagen [2]. Jusman et al. [17] have shown that some gambir bioactive compounds have potency to inhibit platelet-derived growth factor subunit A (PDGFA) as a target protein to reduce fibroblast proliferation through in silico study. Therefore, gambir is assumed to be a promising herb-based antifibrotic agent to deal with keloid disease.
Recent studies have highlighted the efficiency of bioinformatic analysis as a beneficial tool to explore the molecular mechanism of disease [18]; one form of this is epistatic study. Epistatic study is a beneficial tool to discover a new drug candidate for human disease by understanding the gene interaction and pathway that contribute to the disease risk. Furthermore, the researcher can identify potential targets that can disrupt the interaction and reduce the disease progression [19,20,21]. One branch of epistatic interaction study is network pharmacology, which can be used as an appropriate approach to analyze and model the interactions between drug–target networks [22]. To date, a comprehensive understanding of gambir’s mechanism of action against fibrotic tissue, particularly keloid, is still lacking. The conventional experimental method to explore all the proteins included in the mechanism is time-consuming, expensive, and so forth [23]. Therefore, this bioinformatic approach is suggested as an efficient pathway to analyze the interactions between keloid-related proteins and the bioactive compounds of gambir. This study aims to investigate the bioactive compounds of gambir and their interaction with keloid target proteins using a network pharmacology and molecular simulation approach. The result of this study is expected to provide a scientific explanation of the mechanism of gambir as an antikeloid candidate.

2. Results

2.1. Screening for Uncaria Gambir Bioactive Compounds and Its Target Proteins

A total 74 bioactive compounds of Uncaria gambir (Hunter) Roxb. were identified from the literature study [9,17,24], with 5048 protein-specific targets based on SWISS target prediction. A total of 277 targets from 29 bioactive compounds were predicted to be targeted to keloid-related target proteins. The drug-likeness nature, pharmacokinetic properties, and toxicity of the compounds are shown in Table 1. The drug-likeness properties identified from this study were the molecular weight (Mw), octanol–water partition coefficients (Log P), hydrogen-bond acceptor, and donor (HBA and HBD). The pharmacokinetic parameters included were the gastrointestinal absorption (GI Abs) and cytochrome (CYP) inhibition. The toxicity of the compounds was described by the toxicity class (1–6) according to the globally harmonized system of classification for the labeling of chemicals (GHS). Classes 1–2 mean fatal if swallowed, Class 3 is toxic if swallowed, Class 4 is harmful if swallowed, Class 5 may be harmful if swallowed, and Class 6 is non-toxic. All of the compounds were categorized in Classes 3 to 6.

2.2. Screening for Keloid-Related Target Proteins

A total of 96 out of 695 keloid-related target proteins were identified after combining the 2 databases. The 96 genes were then imported into Cytoscape 3.9.1 to plot the keloid-related target network diagram. A total of 96 nodes and 100 edges were obtained using the Cytoscape 3.9.1 platform, with a high confidence score (0.98). The top 11 of the 96 targets were selected based on the degree, betweenness centrality (BC), and closeness centrality (CC) score (Table 2). The keloid-related protein–protein interaction (PPI) network was visualized using GeneMANIA (Figure 1).

2.3. Interaction of Bioactive Compounds of Gambir on Keloid-Related Target Proteins

A total of 277 targets from 29 bioactive compounds were predicted to be targeted to keloid-related target proteins. All of the gambir bioactive compounds’ targets were merged with the keloid-related targets and keloid-related targets with a high degree value (HDV) to further discover the potential targets and mechanism (Figure 2A). Based on the diagram, there were 11 target proteins defined in the intersection area. All of these target proteins were defined as the main targets of gambir bioactive compounds in treating keloid disease. The drug–protein interaction (DPI) network of gambir’s bioactive compounds targeted to keloid-related target proteins is shown in Figure 2B. The pink dots represent the keloid protein targets, the orange dots represent gambir bioactive compounds, and the connecting lines represent the compound–target interactions. Gambiriin A1, Procyanidin B1, Procyanidin B3, Gambiriin A2, Epicatechin gallate, Gambiriin A3, Gambiriin C, Isogambirine, Procyanidin B2, and Procyanidin B4 are the top 10 bioactive compounds of gambir with high degree values resulting from the DPI network analysis (Table S3, Supplementary).

2.4. Gene Ontology (GO) and Pathway Enrichment Analysis

GO enrichment analysis of the 11 target proteins in the intersection area was performed using the DAVID database. The result discovered 111 biological processes (BPs), 12 cellular components (CCs), and 17 molecular functions (MFs). Table 3 shows the sorted BPs, CCs, and MFs by degree of significance, p < 0.05. The BP of the gambir compound–keloid–target is significantly enriched in the regulation of cell proliferation, transcription, protein phosphorylation, cell migration, apoptotic process, signal transduction, peptidyl-serine phosphorylation, and protein kinase B signaling. On the other hand, the MF is mainly enriched in identical protein and enzyme binding, growth factor and cytokine activity, protein (serine/threonine/tyrosine) kinase binding and activity, and protein homodimerization activity. The CC is primarily enriched in cytoplasm, nucleus, extracellular space, extracellular matrix, extracellular region, platelet alpha granule lumen, cell surface, collagen trimer, and secretory granule.
The pathway analysis obtained 211 hits from REACTOME and 146 hits from KEGG, respectively. The first five pathways were screened and sorted according to the adjusted p < 0.05, which most likely shared comparable signaling pathways. The top five signaling pathways discovered from REACTOME are signaling by interleukin (IL; IL-4 and IL-13), cytokine signaling in the immune system, transcriptional regulation of transcription factors or growth factors. Meanwhile, the AGE_RAGE signaling pathway in diabetic complication, pathway in cancer, HIF-1 signaling pathway, proteoglycan in cancer, FOXO signaling pathway are the top signaling pathways discovered from KEGG. The pathway prediction of the target proteins using KEGG is illustrated in Figure 3. This analysis showed results corresponding to the GO biological process and molecular function analyses. This indicates that the bioactive compounds of Uncaria gambir might treat keloid by acting on these pathways.

2.5. Molecular-Docking Analysis

Molecular docking was used to validate the interaction between the gambir compounds and targets. The top 10 gambir bioactive compounds with high degree values were docked with TGFβ1, AKT1, and MMP1. Some compounds had lower binding energy compared with TAC as the positive control, as shown in Table 4 and visualized in Figure 4. The lower the binding energy, the more stable the binding between the ligand and the protein [25]. It was assumed that the compounds may affect the keloid pathomechanism through these target proteins. Gambiriin C and TGFβ1, Isogambirine and AKT1, and Procyanidin B1 and MMP1 are the compound–protein complexes with the highest binding energy.

3. Discussion

Keloid is defined as a skin fibrosis-associated disorder that is characterized by excessive growth of fibroblasts, leading to aberrant accumulation of extracellular matrix (ECM) [2]. Herbal-based therapy is expected as an alternative treatment with low toxicity and side effects in treating keloid. Uncaria gambir (Hunter) Roxb. is an original herbal medicine from Indonesia that has potential as an antifibrosis agent. This study first identified the bioactive compounds of gambir related to keloid target proteins through an epistatic interaction and molecular simulation approach. Network pharmacology is one of the branches of epistatic study that is considered the next paradigm in discovering the pharmacology of certain herbal compounds and their relationship to target genes of the disease. The relationship between herbs, disease, and molecular targets is discovered using a network basis that represents the action of the compound in the disease pathomechanism. Therefore, this method is considered the ideal approach for studying an herbal-based medicine that has multiple compounds and multiple targets [26,27].
This study found that there are 10 bioactive compounds of gambir with high degree values (cut-off degree value ≥ 5). Molecular-docking validation was conducted on these bioactive compounds and target proteins (TGFβ1, AKT1, and MMP1). Gambiriin C and TGFβ1, Isogambirine and AKT1, and Procyanidin B1 and MMP1 were identified as the compound–protein complexes with the highest binding energy (Table 4 and Figure 4). Therefore, Gambiriin C, Isogambirine, and Procyanidin B1 were proposed to be the most promising gambir compounds to treat keloid disease. All these compounds are classified as flavonoids, particularly chalcane-flavan dimers [28]. Flavonoids have a number of medicinal benefits, including antioxidant and anti-inflammatory properties [29]. In addition, these bioactive compounds have approved drug-likeness, pharmacokinetics, and toxicity properties as a drug candidate. Gambir’s bioactive compounds show adverse properties to ideal categories, particularly molecular weight. In term of the pharmacokinetics properties, most of the compounds are highly absorbed in GI tracts and inhibit CYP3A4 enzyme. Most of the compounds are categorized in class 4 to 6 predicted toxicity classes, meaning that most of the compounds have middle to low toxic properties.
The keloid-related proteins were represented in the PPI network produced from the datamining of two databases. DisGeNet and GeneCards are sophisticated search engines that integrate and standardize data about disease-associated genes from multiple sources, including the scientific literature [30,31]. The other network was constructed by merging the previous PPI network and gambir–target network, resulting in the DPI network (Figure 2B). The networks were then analyzed to identify the most impactful proteins in the network based on three important values. The first value is the degree of a node; this referred to the number of edges between nodes in the network. The greater the degree, the more important the nodes in the network. The second is the BC; this represents the signal probability passed through a node. The higher the BC value, the more important the node. The third is the CC value, which reflects the network tightness. The tighter the network, the higher the efficiency [8]. Based on the analysis, 11 targets were identified as the main targets of gambir bioactive compounds related to keloid disease: EP300, AKT1, STAT3, TGFB1, VEGFA, EGFR, CCND1, IGF1R, TIMP1, MMP1, and PDGFA. Thus, it is speculated that the bioactive compounds of gambir may have pharmacological activity in keloid through these targets.
GO enrichment analysis of the 11 potential protein targets was performed to clarify the role of the gambir-targeted proteins in relation to the gene function. The biological process showed that regulation of cell proliferation, transcription, protein phosphorylation, cell migration, apoptotic process, signal transduction, peptidyl-serine phosphorylation, and protein kinase B signaling are the main aspects whereby gambir may play a crucial role in dealing with keloids. In line with this, the molecular function was mainly enriched in identical protein and enzyme binding, growth factor and cytokine activity, protein (serine/threonine/tyrosine) kinase activity, and protein homodimerization activity. Interestingly, the enrichment result shared the same main characteristic of keloid disease. The main characteristic of keloid is over-synthesis and accumulation of extracellular matrix (ECM) due to effusive fibroblast proliferation and differentiation. The fibroblast is the key player in tissue fibrosis, including keloid. The activation of keloid fibroblasts is stimulated by the paracrine effect of cytokines and growth factor in the keloid’s surrounding micro-environment [2,4,32]. Dysregulation of cytokines and growth factor at any stage of the wound-healing process may lead to prolonged inflammation as a favorable condition for keloid formation [33].
In parallel with the cellular component enrichment, the protein targets were primarily enriched in cytoplasm, nucleus, extracellular space, extracellular matrix, extracellular region, platelet alpha granule lumen, cell surface, collagen trimer, and secretory granule. The highly proliferated fibroblast and its activation into a myofibroblast leads to excessive ECM synthesis and deposition, particularly collagen [34]. Keloid scars have a 20-fold increase in collagen production compared with normal skin [3]. Overproduction of ECM is closely related to prolonged inflammation in keloid pathogenesis. Keloid formation is driven by multiple cytokines, growth factors, and transcription factors [2,35]. The study by Nangole et al. [36] revealed that the high expression of proinflammatory cytokines in keloid impacts the sensitivity and increases the production of collagen on keloid fibroblasts. This is in accordance with the result of the pathway enrichment analysis from REACTOME that showed signaling by cytokine, transcription factor, and growth factor were predicted as the main pathways of gambir on keloid. Some of the cytokines and growth factors counted as main players in tissue fibrosis are interleukin-4 (IL-4), IL-13, PDGF, and TGF-β1 [37].
The closest pathway prediction was generated from KEGG, which showed TGFB1 as one of the frontline proteins in the pathway (Figure 3). TGF-β is the key activator of fibroblasts in response to fibrosis. TGF-β could be released from epithelial or endothelial cells, resident fibroblasts, platelets, and macrophages [38]. TGF-β initiates a downstream signaling cascade through activation of the Smad2/3 protein and then forms a complex protein with Smad4. This complex protein migrates to the nucleus as the regulator of ECM organization, fibroblast proliferation, and differentiation [38,39]. From the same group of growth factors, PDGF regulates fibroblast activity through the phosphatidylinositol3-kinase–AKT (PI3K/Akt) pathway [40]. The activation of the catalytic subunit of PI3K undergoes enzymatic function to catalyze phosphorylation of phosphatidylinositol 4,5-bisphosphate to phosphatidylinositol 3,4,5-trisphosphate. These triphosphate molecules activate the Akt/PKB protein to induce cell proliferation and inhibit apoptosis [41]. Most of growth factor signaling pathways are categorized as protein kinase receptors with homodimerization activity, including PDGF, EGF, VEGF, and TGF-β receptor [42]. This is in accordance with the result of the molecular function enrichment analysis.
The TGF-β pathway is associated with epithelial mesenchymal transition (EMT) activation in the keloid pathomechanism. This is a process whereby epithelial cells lose adhesion and acquire expression of a mesenchymal component and manifest the increase of cell migration. EMT plays an important role in fibroblasts differentiation into myofibroblasts. There are four common sources of myofibroblasts: fibrocytes, resident fibroblasts, endothelial cells, and epithelial cells [37]. EMT may be induced by epidermal growth factor (EGF) and EGFR via activation of Smad2/3 [43]. Another study proved that the EGF/EFGR signaling pathway promotes EMT mediated by STAT3 [44,45]. Research by Lee et al. [46] on STAT3 found signaling was the most significantly enriched Gene Ontology term in keloid fibroblast. Along with the activation of EMT, fibroblast differentiation and activation impact the imbalance of ECM remodeling.
Abnormal ECM reconstruction during wound healing contributes to keloid formation. Unlike normal scars, a decreased level of MMP1 and an increased level of TIMP1 are found in keloids that are related to tissue fibrosis [47]. MMP1 is one of the type of collagenase that catalyzes the cleave of types I and III collagens. Collagen I and III are highly expressed in keloid tissue [48,49]. Meanwhile, TIMP1 is the endogenous inhibitor of MMP1. The imbalance between the expression of MMP and TIMP leads to excessive collagen I and III synthesis and deposition in fibrotic tissue like keloid [47]. The expression of MMP and TIMP is regulated by TGF-β. The inhibition of TGF-β and the subsequent downstream signaling contribute to a more balanced MMP1/TIMP1 ratio and may be a promising therapeutic approach for keloid treatments [50].
This recent study revealed the multicomponent, multitarget, and multipathway mechanism of Uncaria gambir bioactive compounds against keloid through network interaction and molecular-docking simulation. EP300, AKT1, STAT3, TGFB1, VEGFA, EGFR, CCND1, IGFR, TIMP, MMP, and PDGFA are defined as the core proteins targeted by gambir’s bioactive compounds through profibrotic cytokines and growth factor signaling pathways. On the other hand, other signaling pathways should not be ignored. Therefore, the other pathway enrichment analysis result is worthy of further study. In addition, further experimental research is required to analyze and validate the mechanism of gambir in keloid identified in this study.

4. Materials and Methods

4.1. Mining Databases of Gambir Bioactive Components and Target Protein

The bioactive compounds of Uncaria gambir were obtained from a literature study. The structure and canonical simplified molecular-input line-entry (SMILE) were accessed from PubChem (https://pubchem.ncbi.nlm.nih.gov) between 4–8 January 2023. The bioactive compounds that were predicted to be targeted to keloid-related target proteins were listed for further analysis. The data were then submitted to SwissADME (http://www.swissadme.ch/) on 14 July 2023 to predict the absorption, distribution, metabolism, and excretion (ADME) parameters, pharmacokinetics properties, and drug-likeness nature. The Mw, HBA and HBD and Log P were utilized as parameters to consider the drug-likeness assessment. The Log P value was presented as the consensus the partition coefficient between n-octanol and water (LogPo/w), which was generated from SwissADME [51]. The prediction of the toxicity profile of the compounds was analyzed using a toxicity database, Protox II (https://tox-new.charite.de/protox_II/) (accessed on 16 July 2023). Furthermore, SWISS target prediction (http://www.swisstargetprediction.ch/), searching for target proteins corresponding to each molecule with only human (Homo sapiens) genes, was set as data retrieval. The list of proteins was aligned using UniProt (http://www.uniprot.org/uploadlists/) to eliminate double-listed proteins (accessed on 20 June 2023).
The protein targets related to keloid were identified in the Gene Cards (https://www.genecards.org/) and DisGeNet (https://www.disgenet.org/) databases using “Keloid” as the keyword in the search engine column (accessed on 13 February 2023). The identified proteins were filtered to only Homo sapiens genes. We utilized a Venn script (https://bioinfogp.cnb.csic.es/tools/venny/) to identify keloid-related protein targets at the intersection of the databases [52]. This tool was also used for the gambir bioactive compound targets, keloid-related targets, and protein targets related to keloid with high degree values intersection lists.

4.2. Construction of Protein–Protein Interaction (PPI) and Drug–Protein Interaction (DPI) Network

The PPI and DPI networks were constructed, visualized, and analyzed using Cytoscape 3.9.1 software on 15 February 2023. The PPI network was constructed from the keloid-related protein list. The proteins were uploaded to the String database (https://string-db.org/) on 20 February 2023, with 0.98 as the high confidence score for protein interaction and Homo sapiens set as the only organism. STRING suggests a threshold of 0.15 to be considered as low confidence, 0.40 as medium confidence, and 0.90 as high confidence. The higher confidence score of the interaction reflects higher supporting evidence and reduced false positives. On the other hand, the higher confidence may reduce the sensitivity and potentially exclude valid interactions inadvertently [23,53].
Afterwards, the network was analyzed to discover the top rank of significance proteins based on the degree, betweenness centrality (BC), and closeness centrality (CC) scores. The network was then visualized using GeneMANIA (http://genemania.org). This website allows users to find genes likely to share functions based on the interaction [54]. The second network is drug–protein interaction (DPI). The network was created to visualize and analyze the interactions of the bioactive compounds of gambir with keloid-related target proteins. Furthermore, the networks were merged to determine a new network of compound–targets of gambir targeted to keloid disease. The screening of the key bioactive compounds was also subjected to this step. The bioactive compounds with high degree values were considered the key important compounds [55].

4.3. Gene Annotation (GO) and Pathway Enrichment Analyses

The core protein targets from the previous step were uploaded to perform enrichment analysis for Gene Ontology (GO), biological process (BP), cellular component (CC), and molecular function (MF), as well as the pathway. All of these were conducted to describe the gene set mechanism of action through the Database for Annotation, Visualization, and Integrated Discovery (DAVID, https://david.ncifcrf.gov/), REACTOME (https://reactome.org/), and Kyoto Encyclopedia of Genes and Genomes (KEGG) (https://www.genome.jp/kegg/). The list of genes was based on the adjusted p value of less than 0.05 as the cut-off value that indicates statistical significance [56,57]. All of these analysis were accessed between 20 June 2023 and 7 July 2023.

4.4. Molecular Docking

The interactions between the gambir bioactive compounds and target proteins were validated using molecular-docking analysis. The top 10 gambir bioactive compounds with HDVs were searched for and downloaded in SDF format for the 2D structure in PubChem (https://pubchem.ncbi.nlm.nih.gov/) on 25 July 2023. The 2D structures were then prepared and saved as pdb format using MarvinSketch 23.11 software. TGFB1, AKT1, and MMP1 were selected as the target proteins, considering that these proteins were in the up, middle, and downstream of the KEGG pathway prediction. Triamcinolone acetonide (TAC) was used as a positive control since this chemical has been widely used for keloid treatment due to its anti-inflammatory properties [58]. The pdb formats of the 3D structure were downloaded from the RCSB PDB (https://www.rcsb.org/) database on 28 July 2023. Moreover, 3KFD, 4GV1, and 1HFC were the PDB ID for TGFβ1, AKT1, and MMP1, respectively. All the protein and ligand structures were saved as pdbqt formats. Molecular docking between the gambir bioactive compounds and target proteins was conducted using AutoDockTools 1.5.6 software on 9 August 2023. The ligands were set as a rigid structure and positioned at the central of cavity, which was docked by the native ligand corresponding to each protein model. Afterwards, the grid box was adjusted for the x-, y-, and z-axis for each protein models and the grid box with the lowest binding energy was selected for further analysis. The Lamarckian Genetic Algorithm was employed for running the molecular-docking simulation. The absolute value of the docking score was identified as the binding affinity and stability between the compound and target. Discovery Studio 4.5. was used to visualize the result of the molecular-docking analyses on 15 August 2023 [59].

5. Conclusions

In conclusion, this study first predicted the interactions of the bioactive compounds of gambir with keloid target proteins through an epistatic interaction and molecular simulation approach. EP300, AKT1, STAT3, TGFB1, VEGFA, EGFR, CCND1, IGFR, TIMP, MMP, and PDGFA were defined as the core proteins targeted by gambir’s bioactive compounds with high degree values. Gambiriin C and TGFβ1, Isogambirine and AKT1, and Procyanidin B1 and MMP1 were identified as the compound–protein complexes with the highest binding energy. GO enrichment and pathway analysis found that gambir’s bioactive compounds may act on keloid-related target proteins to regulate cell proliferation, migration, transcription, and signal transduction activity via cytokine, transcription, and growth factor signaling pathways. This is proposed to be correlated to the dysregulation of cytokines and growth factors in prolonged inflammation as favorable conditions for keloid formation. This study provides a reference for potential targets and pathways to explain the molecular mechanism of gambir’s bioactive compounds against keloid disease.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ph17060662/s1, Table S1. Keloid-related target proteins, Table S2. Gambir’s bioactive compound target proteins. Table S3. Degree ranking of Gambir bioactive compound analysed by Cytoscape 3.9.1.

Author Contributions

Conceptualization, S.S.N. and F.F.; methodology, S.S.N. and F.F.; validation, F.F., S.W.A.J. and R.S.; formal analysis, S.S.N.; investigation, S.S.N.; funding acquisition, R.S. and T.Y.; writing—original draft preparation, S.S.N.; writing—review and editing, F.F., S.W.A.J., R.S. and T.Y.; visualization, S.S.N. and R.S.; supervision, F.F., S.W.A.J., R.S. and T.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by Universitas Indonesia (PUTI 2023 grant; grant number: NKB-400/UN2.RST/HKP.05.00/2023).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data are available in this manuscript.

Acknowledgments

The support from Universitas Indonesia, Centre or Education Services (Ministry of Education, Culture, Research, and Technology of the Republic of Indonesia), and Indonesia Endowment Funds for Education (LPDP) is gratefully acknowledged.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Elsaie, M.L. Update on management of keloid and hypertrophic scars: A systemic review. J. Cosmet. Dermatol. 2021, 20, 2729–2738. [Google Scholar] [CrossRef] [PubMed]
  2. Andrews, J.P.; Marttala, J.; Macarak, E.; Rosenbloom, J.; Uitto, J. Keloids: The paradigm of skin fibrosis—Pathomechanisms and treatment. Matrix Biol. 2016, 51, 37–46. [Google Scholar] [CrossRef] [PubMed]
  3. Berman, B.; Maderal, A.; Raphael, B. Keloids and hypertrophic scars: Pathophysiology, classification, and treatment. Dermatol. Surg. 2017, 43, S3–S18. [Google Scholar] [CrossRef] [PubMed]
  4. Ashcroft, K.J.; Syed, F.; Bayat, A. Site-specific keloid fibroblasts alter the behaviour of normal skin and normal scar fibroblasts through paracrine signalling. PLoS ONE 2013, 8, e75600. [Google Scholar] [CrossRef] [PubMed]
  5. Ud-Din, S.; Bayat, A. Keloid scarring or disease: Unresolved quasi-neoplastic tendencies in the human skin. Wound Repair Regen. 2020, 28, 422–426. [Google Scholar] [CrossRef]
  6. Ekstein, S.F.; Wyles, S.P.; Moran, S.L.; Meves, A. Keloids: A review of therapeutic management. Int. J. Dermatol. 2021, 60, 661–671. [Google Scholar] [CrossRef] [PubMed]
  7. Naik, P.P. Novel targets and therapies for keloid. Clin. Exp. Dermatol. 2022, 47, 507–515. [Google Scholar] [CrossRef] [PubMed]
  8. Kang, P.; Wu, Z.; Zhong, Y.; Wang, Z.; Zhou, C.; Huo, S.; Guo, H.; Li, S.; Xu, K.; Liu, L.; et al. A Network Pharmacology and Molecular Docking Strategy to Explore Potential Targets and Mechanisms Underlying the Effect of Curcumin on Osteonecrosis of the Femoral Head in Systemic Lupus Erythematosus. BioMed Res. Int. 2021, 2021, 5538643. [Google Scholar] [CrossRef] [PubMed]
  9. Munggari, I.P.; Kurnia, D.; Deawati, Y.; Julaeha, E. Current Research of Phytochemical, Medicinal and Non-Medicinal Uses of Uncaria gambir Roxb.: A Review. Molecules 2022, 27, 6551. [Google Scholar] [CrossRef] [PubMed]
  10. Dharsono, H.D.A.; Wibisono, L.; Hayati, A.T.; Apriyanti, E.; Satari, M.H.; Kurnia, D. Mode action prediction of catechin from Uncaria gambir Roxb. against UDP-N-acetylenolpyruvyl-glucosamine reductase (MurB enzyme) of Streptococcus mutans: In silico study. J. Adv. Pharm. Technol. Res. 2022, 13, 197–201. [Google Scholar] [CrossRef] [PubMed]
  11. Syarifah, S. Anticancer activity of Uncaria gambir roxb on T47D breast cancer cells. J. Phys. 2019, 1317, 012106. [Google Scholar] [CrossRef]
  12. Wardana, A.P.; Aminah, N.S.; Kristanti, A.N.; Zahrah, H.I.; Fahmi, M.Z.; Yuliati, A.; Suhaili; Choudhary, M.I. Anticancer Activity of Nano Uncaria gambir. 2022. Available online: https://www.researchsquare.com/article/rs-1497841/v1 (accessed on 19 September 2022).
  13. Yimam, M.; Lee, Y.C.; Kim, T.W.; Moore, B.; Jiao, P.; Hong, M.; Kim, H.J.; Nam, J.B.; Kim, M.R.; Oh, J.S.; et al. Analgesic and anti-Inflammatory effect of UP3005, a botanical composition Containing two standardized extracts of Uncaria gambir and Morus alba. Pharmacogn. Res. 2015, 7, S39–S46. [Google Scholar] [CrossRef] [PubMed]
  14. Desdiani, D.; Rengganis, I.; Djauzi, S.; Setiyono, A.; Sadikin, M.; Jusman, S.W.A.; Siregar, N.C.; Suradi, S.; Eyanoer, P.C.; Fadilah, F. In vitro assay and study interaction of Uncaria gambir (Hunter) Roxb. as anti-fibrotic activity against A549 cell line. Pharmacogn. J. 2020, 12, 1232–1240. [Google Scholar] [CrossRef]
  15. Desdiani, D.; Rengganis, I.; Djauzi, S.; Setiyono, A.; Sadikin, M.; Jusman, S.W.A.; Siregar, N.C.; Suradi, S.; Eyanoer, P.C. Fibropreventive and Antifibrotic Effects of Uncaria gambir on Rats with Pulmonary Fibrosis. Evid.-Based Complement. Altern. Med. 2022, 2022, 6721958. [Google Scholar] [CrossRef] [PubMed]
  16. Ningsih, S. Efek Hepatoprotektor Gambir (Uncaria gambir (Hunter) Roxb.) dalam Menghambat Pembentukan Pembentukan Kolagen dengan Menekan TiMP-1 In Vivo. Ph.D. Thesis, Universitas Indonesia, Kota Depok, Indonesia, 2015. [Google Scholar]
  17. Jusman, S.W.A.; Amalia, M.F.; Paramita, R.; Ningsih, S.; Fadilah, F. Structure-based screening of inhibitor platelet-derived growth factor from ethanol extract of Uncaria gambir (Hunter) Roxb. as an antifibrotic in keloid fibroblast cells. Appl. Pharm. Sci. 2023, 13, 36–43. [Google Scholar] [CrossRef]
  18. Wang, Q.; Zhong, Y.; Li, Z.; Zhu, D.; Lu, H.; Chen, P.; Li, C.; Peng, X.; Li, Q.; Zeng, K. Multitranscriptome analyses of keloid fibroblasts reveal the role of the HIF-1α/HOXC6/ERK axis in keloid development. Burn. Trauma 2022, 10, tkac013. [Google Scholar] [CrossRef] [PubMed]
  19. Nagel, R.L. Epistasis and the genetics of human diseases. Comptes Rendus Biol. 2005, 328, 606–615. [Google Scholar] [CrossRef] [PubMed]
  20. Quan, Y.; Liu, M.-Y.; Liu, Y.-M.; Zhu, L.-D.; Wu, Y.-S.; Luo, Z.-H.; Zhang, X.-Z.; Xu, S.-Z.; Yang, Q.-Y.; Zhang, H.-Y. Facilitating anti-cancer combinatorial drug discovery by targeting epistatic disease genes. Molecules 2018, 23, 736. [Google Scholar] [CrossRef] [PubMed]
  21. Wilkins, J.F.; Cannataro, V.L.; Shuch, B.; Townsend, J.P. Analysis of mutation, selection, and epistasis: An informed approach to cancer clinical trials. Oncotarget 2018, 9, 22243–22253. [Google Scholar] [CrossRef] [PubMed]
  22. Zhou, Y.; Lauschke, V.M. Pharmacogenomic network analysis of the gene-drug interaction landscape underlying drug disposition. Comput. Struct. Biotechnol. J. 2020, 18, 52–58. [Google Scholar] [CrossRef] [PubMed]
  23. Ding, Z.; Kihara, D. Computational identification of protein-protein interactions in model plant proteomes. Sci. Rep. 2019, 9, 8740. [Google Scholar] [CrossRef] [PubMed]
  24. Nur Sazwi, N.; Nalina, T.; Rahim, Z.H.A. Antioxidant and cytoprotective activities of Piper betle, Areca catechu, Uncaria gambir and betel quid with and without calcium hydroxide. BMC Complement. Altern. Med. 2013, 13, 351. [Google Scholar] [CrossRef] [PubMed]
  25. Ye, J.; Li, L.; Hu, Z. Exploring the Molecular Mechanism of Action of Yinchen Wuling Powder for the Treatment of Hyperlipidemia, Using Network Pharmacology, Molecular Docking, and Molecular Dynamics Simulation. BioMed Res. Int. 2021, 2021, 9965906. [Google Scholar] [CrossRef] [PubMed]
  26. Wang, Y.; Yang, S.H.; Zhong, K.; Jiang, T.; Zhang, M.; Kwan, H.Y.; Su, T. Network Pharmacology-Based Strategy for the Investigation of the Anti-Obesity Effects of an Ethanolic Extract of Zanthoxylum bungeanum Maxim. Front. Pharmacol. 2020, 11, 572387. [Google Scholar] [CrossRef] [PubMed]
  27. Lu, Z.; Huang, M.; Lin, H.; Wang, G.; Li, H. Network pharmacology and molecular docking approach to elucidate the mechanisms of Liuwei Dihuang pill in diabetic osteoporosis. J. Orthop. Surg. Res. 2022, 17, 314. [Google Scholar] [CrossRef] [PubMed]
  28. Taniguchi, S.; Kuroda, K.; Doi, K.-I.; Tanabe, M.; Shibata, T.; Yoshida, T.; Hatano, T. Revised structures of gambiriins A1, A2, B1, and B2, chalcane-flavan dimers from gambir (Uncaria gambir extract). Chem. Pharm. Bull. 2007, 55, 268–272. [Google Scholar] [CrossRef] [PubMed]
  29. Ullah, A.; Munir, S.; Badshah, S.L.; Khan, N.; Ghani, L.; Poulson, B.G.; Emwas, A.-H.; Jaremko, M. Important Flavonoids and Their Role as a Therapeutic Agent. Molecules 2020, 25, 5243. [Google Scholar] [CrossRef] [PubMed]
  30. Piñero, J.; Ramírez-Anguita, J.M.; Saüch-Pitarch, J.; Ronzano, F.; Centeno, E.; Sanz, F.; Furlong, L.I. The DisGeNET knowledge platform for disease genomics: 2019 update. Nucleic Acids Res. 2020, 48, D845–D855. [Google Scholar] [CrossRef] [PubMed]
  31. Safran, M.; Dalah, I.; Alexander, J.; Rosen, N.; Stein, T.I.; Shmoish, M.; Nativ, N.; Bahir, I.; Doniger, T.; Krug, H.; et al. GeneCards Version 3: The human gene integrator. Database 2010, 2010, baq020. [Google Scholar] [CrossRef] [PubMed]
  32. Macarak, E.J.; Wermuth, P.J.; Rosenbloom, J.; Uitto, J. Keloid disorder: Fibroblast differentiation and gene expression profile in fibrotic skin diseases. Exp. Dermatol. 2021, 30, 132–145. [Google Scholar] [CrossRef] [PubMed]
  33. Wang, Z.-C.; Zhao, W.-Y.; Cao, Y.; Liu, Y.-Q.; Sun, Q.; Shi, P.; Cai, J.-Q.; Shen, X.Z.; Tan, W.-Q. The Roles of Inflammation in Keloid and Hypertrophic Scars. Front Immunol. 2020, 11, 603187. [Google Scholar] [CrossRef] [PubMed]
  34. Mcginty, S.; Siddiqui, J.W. Keloid; StatPearls Publishing: St. Petersburg, FL, USA, 2022. [Google Scholar]
  35. Lee, C.C.; Tsai, C.H.; Chen, C.H.; Yeh, Y.C.; Chung, W.H.; Chen, C.B. An updated review of the immunological mechanisms of keloid scars. Front. Immunol. 2023, 14, 1117630. [Google Scholar] [CrossRef] [PubMed]
  36. Nangole, F.W.; Agak, G.W. Keloid pathophysiology: Fibroblast or inflammatory disorders?: Pathophysiology of Keloids. JPRAS Open 2019, 22, 44–54. [Google Scholar] [CrossRef] [PubMed]
  37. Borthwick, L.A.; Wynn, T.A.; Fisher, A.J. Cytokine mediated tissue fibrosis. Biochim. Biophys. Acta Mol. Basis Dis. 2013, 1832, 1049–1060. [Google Scholar] [CrossRef] [PubMed]
  38. Varmeh, S.; Egia, A.; McGrouther, D.; Tahan, S.R.; Bayat, A.; Pandolfi, P.P. Cellular Senescence as a Possible Mechanism for Halting Progression of Keloid Lesions. Genes Cancer 2011, 2, 1061–1066. [Google Scholar] [CrossRef] [PubMed]
  39. Deng, C.-C.; Hu, Y.-F.; Zhu, D.-H.; Cheng, Q.; Gu, J.-J.; Feng, Q.-L.; Zhang, L.-X.; Xu, Y.-P.; Wang, D.; Rong, Z.; et al. Single-cell RNA-seq reveals fibroblast heterogeneity and increased mesenchymal fibroblasts in human fibrotic skin diseases. Nat Commun. 2021, 12, 1–16. [Google Scholar] [CrossRef] [PubMed]
  40. Papadopoulos, N.; Lennartsson, J. The PDGF/PDGFR pathway as a drug target. Mol. Aspects Med. 2018, 62, 75–88. [Google Scholar] [CrossRef] [PubMed]
  41. Klinkhammer, B.M.; Floege, J.; Boor, P. PDGF in organ fibrosis. Mol. Aspects Med. 2018, 62, 44–62. [Google Scholar] [CrossRef] [PubMed]
  42. Dowling, C.M.; Kiely, P.A. Targeting protein Kinase C downstream of growth factor and adhesion signalling. Cancers 2015, 7, 1271–1291. [Google Scholar] [CrossRef] [PubMed]
  43. Kim, J.; Kong, J.; Chang, H.; Kim, H.; Kim, A. EGF induces epithelial-mesenchymal transition through phosphor-Smad2/3-Snail signaling pathway in breast cancer cells. Oncotarget 2016, 7, 85021–85032. [Google Scholar] [CrossRef] [PubMed]
  44. Lo, H.-W.; Hsu, S.-C.; Xia, W.; Cao, X.; Shih, J.-Y.; Wei, Y.; Abbruzzese, J.L.; Hortobagyi, G.N.; Hung, M.-C. Epidermal Growth Factor Receptor Cooperates with Signal Transducer and Activator of Transcription 3 to Induce Epithelial-Mesenchymal Transition in Cancer Cells via Up-regulation of TWIST Gene Expression. Cancer Res. 2008, 23, 9066–9076. [Google Scholar] [CrossRef] [PubMed]
  45. Lee, Y.I.; Shim, J.E.; Kim, J.; Lee, W.J.; Kim, J.W.; Nam, K.H.; Lee, J.H. WNT5A drives interleukin-6-dependent epithelial–mesenchymal transition via the JAK/STAT pathway in keloid pathogenesis. Burn. Trauma. 2022, 10, tkac023. [Google Scholar] [CrossRef] [PubMed]
  46. Lee, Y.S.; Liang, Y.C.; Wu, P.; Kulber, D.A.; Tanabe, K.; Chuong, C.M.; Ming, C.; Randall, W.; Tai-Lan, T. STAT3 Signaling Pathway is Implicated in Keloid Pathogenesis by Preliminary Transcriptome and Open Chromatin Analyses. Physiol. Behav. 2017, 174, 139–148. [Google Scholar] [CrossRef] [PubMed]
  47. Xue, M.; Jackson, C.J. Extracellular Matrix Reorganization During Wound Healing and Its Impact on Abnormal Scarring. Adv. Wound Care 2015, 4, 119–136. [Google Scholar] [CrossRef] [PubMed]
  48. Lee, D.E.; Trowbridge, R.M.; Ayoub, N.T.; Agrawal, D.K. High-mobility Group Box Protein-1, Matrix Metalloproteinases, and Vitamin D in Keloids and Hypertrophic Scars. PRS Glob. Open 2015, 3, e425. [Google Scholar] [CrossRef] [PubMed]
  49. Wulandari, E.; Jusman, S.W.A.; Moenadjat, Y.; Jusuf, A.A.; Sadikin, M. Expressions of Collagen I and III in Hypoxic Keloid Tissue. Kobe J. Med. Sci. 2016, 62, 58–69. [Google Scholar]
  50. Dohi, T.; Miyake, K.; Aoki, M.; Ogawa, R.; Akaishi, S.; Shimada, T.; Okada, T.; Hyakusoku, H. Tissue Inhibitor of Metalloproteinase-2 Suppresses Collagen Synthesis in Cultured Keloid Fibroblasts. PRS Glob. Open 2015, 3, e520. [Google Scholar] [CrossRef] [PubMed]
  51. Daina, A.; Michielin, O.; Zoete, V. SwissADME: A free web tool to evaluate pharmacokinetics, druglikeness and medicinal chemistry friendliness of small molecules. Sci. Rep. 2017, 7, 42717. [Google Scholar] [CrossRef] [PubMed]
  52. Venny, O.J. An Interactive Tool for Comparing Lists with Venn’s Diagrams. Available online: https://methods.sagepub.com/reference/dictionary-of-statistics-methodology/n2067.xml (accessed on 14 February 2023).
  53. Bozhilova, L.V.; Whitmore, A.V.; Wray, J.; Reinert, G.; Deane, C.M. Measuring rank robustness in scored protein interaction networks. BMC Bioinform. 2019, 20, 446. [Google Scholar] [CrossRef] [PubMed]
  54. Franz, M.; Rodriguez, H.; Lopes, C.; Zuberi, K.; Montojo, J.; Bader, G.D.; Morris, Q. GeneMANIA update 2018. Nucleic Acids Res. 2018, 46, 60–64. [Google Scholar] [CrossRef] [PubMed]
  55. Cui, Y.; Wang, H.; Wang, D.; Mi, J.; Chen, G.; Li, F.; Wang, Y.; Zhang, Y. Network Pharmacology Analysis on the Mechanism of Huangqi Sijunzi Decoction in Treating Cancer-Related Fatigue. J. Healthc. Eng. 2021, 2021, 9780677. [Google Scholar] [CrossRef] [PubMed]
  56. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Lempicki RA. Nat. Protoc. 2009, 4, 44–57. [Google Scholar] [CrossRef] [PubMed]
  57. Fabregat, A.; Sidiropoulos, K.; Viteri, G.; Marin-Garcia, P.; Ping, P.; Stein, L.; D’eustachio, P.; Hermjakob, H. Reactome diagram viewer: Data structures and strategies to boost performance. Bioinformatics 2018, 34, 1208–1214. [Google Scholar] [CrossRef] [PubMed]
  58. Coppola, M.M.; Salzillo, R.; Segreto, F.; Persichetti, P. Triamcinolone acetonide intralesional injection for the treatment of keloid scars: Patient selection and perspectives. Clin. Cosmet. Investig. Dermatol. 2018, 11, 387–396. [Google Scholar] [CrossRef] [PubMed]
  59. Dong, Y.; Zhao, Q.; Wang, Y. Network pharmacology-based investigation of potential targets of astragalus membranaceous-angelica sinensis compound acting on diabetic nephropathy. Sci. Rep. 2021, 11, 19496. [Google Scholar] [CrossRef] [PubMed]
Figure 1. PPI network for the candidate targets. Black node indicates target protein and the correlation is represented as connecting colors.
Figure 1. PPI network for the candidate targets. Black node indicates target protein and the correlation is represented as connecting colors.
Pharmaceuticals 17 00662 g001
Figure 2. (A) Intersection Venn diagram between the predicted target proteins of gambir bioactive compounds, keloid-related target proteins, and high degree of connection values of keloid-related target proteins. (B) Compound–target network of gambir bioactive compounds targeted to keloid disease.
Figure 2. (A) Intersection Venn diagram between the predicted target proteins of gambir bioactive compounds, keloid-related target proteins, and high degree of connection values of keloid-related target proteins. (B) Compound–target network of gambir bioactive compounds targeted to keloid disease.
Pharmaceuticals 17 00662 g002
Figure 3. Pathway prediction of the target proteins using KEGG. The red nodes are potential target proteins of Uncaria gambir, while the green nodes are relevant targets in the pathway.
Figure 3. Pathway prediction of the target proteins using KEGG. The red nodes are potential target proteins of Uncaria gambir, while the green nodes are relevant targets in the pathway.
Pharmaceuticals 17 00662 g003
Figure 4. Three best molecular docking results: TGFB1 and Gambiriin C (A), AKT1 and Isogambirine (B), and MMP1 and ProcyanidinB1 (C).
Figure 4. Three best molecular docking results: TGFB1 and Gambiriin C (A), AKT1 and Isogambirine (B), and MMP1 and ProcyanidinB1 (C).
Pharmaceuticals 17 00662 g004
Table 1. The drug-likeness nature, pharmacokinetics properties, and toxicity of Uncaria gambir bioactive compounds.
Table 1. The drug-likeness nature, pharmacokinetics properties, and toxicity of Uncaria gambir bioactive compounds.
NoMoleculeDrug-LikenessPharmacokineticsPredicted Toxicity Class
MwLog PHBAHBDGI AbsCYP Inhibitor
1Catechol110.110.9422HighNo3
2Gambiriin C562.521.84119LowCYP3A45
3Gambiriin A1580.541.691211LowCYP3A45
4Gambiriin A3580.541.451211LowCYP3A45
5Procyanidin B3578.521.531210LowCYP3A45
6Procyanidin B1578.521.531210LowCYP3A45
7Gambiriin A2580.541.691211LowCYP3A45
8Procyanidin B2578.521.531210LowCYP3A45
9Procyanidin B4578.521.531210LowCYP3A45
10Ursolic acid456.705.8832LowNo4
11Quercetin302.241.2375HighCYP2D6CYP3A4CYP3A23
12Quinic acid192.171.7565LowNo6
13Kaempferol286.241.5864HighCYP2D6CYP3A4CYP3A25
14Protocatechuic acid154.120.6543HighCYP3A44
15Epicatechin gallate442.371.25107LowNo4
16Dimeric proanthocyanidin576.501.61129LowCYP3A45
17Gambirflavan D1562.522.29119LowCYP2C9 CYP3A45
18Gambirflavan D2562.522.29119LowCYP2C9CYP3A45
19Hyperoside464.38−0.25128LowNo5
20Isoquercitrin464.38−0.25128LowNo5
21Pyrocatechol110.110.9722HighCYP3A43
22Phloroglucinol126.110.4533HighCYP3A43
234-methoxyphenol124.141.4121HighNo4
244-allyl-2-methoxyphenol164.202.2521HighCYP1A24
25Gambirine384.472.8552HighCYP2D63
26Isogambirine384.502.8252HighCYP2D63
27Tetrahydroalstonine352.432.6741HighCYP2D63
28Dihydrocorynantheine368.473.2241HighCYP2D6 CYP3A43
29Uncariagambiriine620.653.6196LowCYP2C94
Table 2. Top 11 keloid-related core protein targets in the PPI network.
Table 2. Top 11 keloid-related core protein targets in the PPI network.
GeneGene Full NameBCCCDegree
EP300E1A-binding protein p3000.190.4312
AKT1AKT serine/threonine kinase 10.120.4410
STAT3Signal transducer and activator of transcription 30.260.3910
TGFB1Transforming growth factor beta 10.110.448
VEGFAVascular endothelial growth factor A0.170.437
EGFREpidermal growth factor receptor0.180.446
CCND1Cyclin D10.050.414
IGF1RInsulin-like growth factor 1 receptor0.020.393
TIMP1TIMP metallopeptidase inhibitor 10.110.392
MMP1Matrix metallopeptidase 10.140.331
PDGFAPlatelet-derived growth factor alpha0.070.311
Table 3. GO analysis of core target proteins.
Table 3. GO analysis of core target proteins.
CategoryGene FunctionCountp-Value
BPPositive regulation of cell proliferation61.9 × 10−6
BPPositive regulation of transcription from RNA polymerase II promoter61.0 × 10−4
BPPositive regulation of protein phosphorylation51.7 × 10−6
BPPositive regulation of cell migration53.9 × 10−6
BPNegative regulation of apoptotic process55.8 × 10−5
BPPositive regulation of transcription, DNA-templated52.1 × 10−4
BPNegative regulation of transcription from RNA polymerase II promoter56.8 × 10−4
BPSignal transduction51.8 × 10−3
BPPositive regulation of peptidyl-serine phosphorylation48.5 × 10−6
BPPositive regulation of protein kinase B signaling42.3 × 10−5
MFIdentical protein binding65.4 × 10−4
MFEnzyme binding46.1 × 10−4
MFGrowth factor activity32.7 × 10−3
MFCytokine activity33.5 × 10−3
MFProtein kinase activity31.3 × 10−2
MFProtein serine/threonine/tyrosine kinase activity31.7 × 10−2
MFProtein kinase binding32.4 × 10−2
MFProtein homodimerization activity34.5 × 10−2
CCCytoplasm 82.5 × 10−3
CCNucleus 72.2 × 10−2
CCExtracellular space57.1 × 10−3
CCExtracellular matrix41.5 × 10−4
CCExtracellular region46.1 × 10−2
CCPlatelet alpha granule lumen33.7 × 10−4
CCCell surface33.1 × 10−2
CCCollagen trimer24.0 × 10−2
CCSecretory granule24.2 × 10−2
Table 4. Docking scores of Uncaria gambir bioactive compounds with target proteins.
Table 4. Docking scores of Uncaria gambir bioactive compounds with target proteins.
Compound/Affinity (kcal/mol)TGFβ1AKT1MMP1
Gambiriin A1−4.47−7.43−5.32
Procyanidin B1−6.76−3.66−9.54
Procyanidin B3−6.56−6.52−7.49
Procyanidin B2−6.03−5.98−8.38
Gambiriin A2−5.53−5.27−5.99
Epicatechin gallate−5.18−7.79−9.48
Gambiriin A3−5.45−7.68−6.24
Gambiriin C−6.90−7.99−8.46
Isogambirine−6.02−8.68−8.44
Procyanidin B4−6.62−8.47−8.84
TAC−6.19−8.88−6.69
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ningsih, S.S.; Fadilah, F.; Jusman, S.W.A.; Syaidah, R.; Yashiro, T. Profibrotic Inflammatory Cytokines and Growth Factors Are Predicted as the Key Targets of Uncaria gambir (Hunter) Roxb. in Keloids: An Epistatic and Molecular Simulation Approach. Pharmaceuticals 2024, 17, 662. https://doi.org/10.3390/ph17060662

AMA Style

Ningsih SS, Fadilah F, Jusman SWA, Syaidah R, Yashiro T. Profibrotic Inflammatory Cytokines and Growth Factors Are Predicted as the Key Targets of Uncaria gambir (Hunter) Roxb. in Keloids: An Epistatic and Molecular Simulation Approach. Pharmaceuticals. 2024; 17(6):662. https://doi.org/10.3390/ph17060662

Chicago/Turabian Style

Ningsih, Sri Suciati, Fadilah Fadilah, Sri Widia A. Jusman, Rahimi Syaidah, and Takashi Yashiro. 2024. "Profibrotic Inflammatory Cytokines and Growth Factors Are Predicted as the Key Targets of Uncaria gambir (Hunter) Roxb. in Keloids: An Epistatic and Molecular Simulation Approach" Pharmaceuticals 17, no. 6: 662. https://doi.org/10.3390/ph17060662

APA Style

Ningsih, S. S., Fadilah, F., Jusman, S. W. A., Syaidah, R., & Yashiro, T. (2024). Profibrotic Inflammatory Cytokines and Growth Factors Are Predicted as the Key Targets of Uncaria gambir (Hunter) Roxb. in Keloids: An Epistatic and Molecular Simulation Approach. Pharmaceuticals, 17(6), 662. https://doi.org/10.3390/ph17060662

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop