Next Article in Journal
The Entomopathogenic Nematodes H. bacteriophora and S. carpocapsae Inhibit the Activation of proPO System of the Nipa Palm Hispid Octodonta nipae (Coleoptera: Chrysomelidae)
Previous Article in Journal
Identification and Functional Analysis of SabHLHs in Santalum album L.
Previous Article in Special Issue
Detection of Specific Immune Cell Subpopulation Changes Associated with Systemic Immune Inflammation–Index Level in Germ Cell Tumors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

In-Silico Analysis of Deleterious SNPs of FGF4 Gene and Their Impacts on Protein Structure, Function and Bladder Cancer Prognosis

1
Faculty of Applied Sciences, UCSI University, Kuala Lumpur 56000, Malaysia
2
School of Pharmacy, Monash University Malaysia, Jalan Lagoon Selatan, Bandar Sunway 47500, Malaysia
3
Faculty of Health and Life Sciences, INTI International University, Persiaran Perdana BBN, Putra Nilai, Nilai 71800, Malaysia
4
Health Sciences Division, Abu Dhabi Women’s College, Higher Colleges of Technology, Abu Dhabi P.O. Box 41012, United Arab Emirates
5
Centre of Research for Advanced Aquaculture (CORAA), UCSI University, Kuala Lumpur 56000, Malaysia
6
He & Ni Academy, Office Tower B, Northpoint Mid Valley City, Kuala Lumpur 59200, Malaysia
*
Authors to whom correspondence should be addressed.
Life 2022, 12(7), 1018; https://doi.org/10.3390/life12071018
Submission received: 16 June 2022 / Revised: 3 July 2022 / Accepted: 5 July 2022 / Published: 9 July 2022

Abstract

:
Dysregulation of fibroblast growth factors is linked to the pathogenesis of bladder cancer. The role of FGF1 and FGF3 is evident in bladder cancer; however, the role of FGF4 is vague. Despite being reported that FGF4 interacts with FGF1 and FGF3 in MAPK pathways, its pathogenesis and mechanism of action are yet to be elucidated. Therefore, this study aimed to elucidate pathogenic nsSNPs and their role in the prognosis of bladder cancer by employing in-silico analysis. The nsSNPs of FGF4 were retrieved from the NCBI database. Different in silico tools, PROVEAN, SIFT, PolyPhen-2, SNPs&GO, and PhD-SNP, were used for predicting the pathogenicity of the nsSNPs. Twenty-seven nsSNPs were identified as “damaging”, and further stability analysis using I-Mutant 2.0 and MUPro indicated 22 nsSNPs to cause decreased stability (DDG scores < −0.5). Conservation analysis predicted that Q97K, G106V, N164S, and N167S were highly conserved and exposed. Biophysical characterisation indicated these nsSNPs were not tolerated, and protein-protein interaction analysis showed their involvement in the GFR-MAPK signalling pathway. Furthermore, Kaplan Meier bioinformatics analyses indicated that the FGF4 gene deregulation affected the overall survival rate of patients with bladder cancer, leading to prognostic significance. Thus, based on these analyses, our study suggests that the reported nsSNPs of FGF4 may serve as potential targets for diagnoses and therapeutic interventions focusing on bladder cancer.

Graphical Abstract

1. Introduction

Bladder cancer is the top ten most common cancer in the world and is reported to cause 573,000 new cases and 213,000 deaths worldwide in 2020 [1,2]. It affects the urothelial cells that line the urinary bladder. The non-muscle invasive bladder cancer (NMIBC) is prominent in 80% of patients, and muscle-invasive bladder cancer (MIBC) was reported in others [3]. Despite aggressive therapy, up to 50% of NMIBC patients return, and up to 30% develop MIBC [2]. This necessitates a lifetime surveillance cystoscopy to determine the original grade and stage of the disease due to its high recurrence rate.
The pathogenesis of bladder cancer has been studied in different signalling pathways. Recent studies report that upregulation of FGF19-FGFR4 signalling is crucial for carcinogenesis and cancer development. FGFR4 is an attractive target for developing a novel therapeutic, focusing on the FGF19-FGFR4 pathway [4]. MAPK activation due to FGF mutations was reported in 85% of NMIBC cases, emphasising the relevance of this pathway in bladder cancer aetiology [5]. FGFR is a highly conserved receptor tyrosine kinases family that mediates cellular proliferation, differentiation, and death [6]. It has been reported that aberrant FGF signalling promotes tumour development, promoting cancer cell proliferation and survival, and tumour angiogenesis [7]. The mammalian FGF family comprises highly conserved transmembrane tyrosine kinase receptors; FGFR1, FGFR2, FGFR3, and FGFR4. The most prevalent mutation in BC has been reported for FGFR1 and FGF3 with 50–60% NMIBC and 10–15% MIBC, respectively [8,9]. The FGFR mutation is associated with a better prognosis for BC patients, including enhanced survival and a lower chance of recurrence and progression [10]. A recent clinical study reported a 40% response rate on BC treated with a FGFR inhibitor, underscoring the receptor’s potential value as a therapeutic target.
Unlike the other FGF family, the Fibroblast Growth Factor 4 (FGF4) gene does not increase the incidence of cancer, however, it has been reported to be associated with poor prognosis in multiple cancer types [11,12]. The poor prognosis of cancer has been associated with a burden of fatigue, shortness of breath, weight loss, muscle wasting, and pain [13]. This causes cancer progression, side effects of treatment, high levels of systemic inflammation, and malnutrition, leading to poor health-related quality of life (HRQoL).
The prognostic value of FGF4 in bladder cancer remains to be elucidated, and the exact mechanism is yet to be deliberated. Here we first investigated the association between FGF4 and poor prognosis in bladder cancer, and we have also determined the “high-risk” nsSNPs towards the pathogenesis of bladder cancer by employing an in-silico approach. The various bioinformatics steps involved in this study are shown in Figure 1.

2. Materials and Methods

2.1. Retrieving nsSNPs

The SNPs of the FGF4 gene were retrieved from National Center for Biological Information (NCBI) dbSNP database (CRCh37.p13) and the FASTA sequence was obtained from the UniProtKB database. The data were then subjected to multiple bioinformatic tools for SNPs analysis (Figure 1).
A total of 2642 SNPs of FGF4 (Gene ID: 2249) were retrieved from the NCBI dbSNP (https://www.ncbi.nlm.nih.gov/snp/ (accessed on 1 July 2021)), of which 186 were nsSNPs. Information on the allele change, global minor allele frequency (MAF), and residue changes were also retrieved from the database for compilation. The amino acid sequence of the FGF4n FASTA format was obtained from the UniProtKB (UniProt ID: P08620) (https://www.uniprot.org/uniprot/ (accessed on 1 July 2021)).

2.2. Identifying Deleterious nsSNPs

Deleterious nsSNPs were identified using five different in silico tools: PROVEAN (Protein Variation Effect Analyzer) [14] (http://provean.jcvi.org/index.php (accessed on 12 July 2021)) which is embedded with SIFT (Sorting Intolerant from Tolerance) [15,16] (https://sift.bii.a-star.edu.sg/ (accessed on 12 July 2021)), Polyphen-2 (Polymorphism Phenotyping v2) [17] (http://genetics.bwh.harvard.edu/pph2/ (accessed on 15 July 2021)), SNPs&GO (Single Nucleotide Polymorphism and Gene Ontology) [18] (https://snps.biofold.org/snps-and-go//snps-and-go.html (accessed on 19 July 2021)) and PhD-SNPs (Predictor of human Del Single Nucleotide Polymorphism) [18] (https://snps.biofold.org/phd-snp/PhD-SNP.html (accessed on 19 July 2021)). Only nsSNPs predicted to be damaging by all 5 tools were proceeded to further downstream analysis.

2.3. Validating the High-Risk nsSNPs

PMut [19] (http://mmb.irbbarcelona.org/PMut/ (accessed on 26 July 2021)) was used to screen the pathogenicity of the high-risk nsSNPs. This server utilises a predictor engine that was trained using the SwissVar 2016 entries database which contains the mutations that were considered pathogenic and also neutral which were found in over 12,000 proteins. PMut allows its user to train a custom predictor for a more precise prediction using their own datasets but for this study, the default predictor provided by the PMut was used. PMut’s predictor will give the result in the form of prediction scores that ranged from 0 to 1. Mutations with scores between 0 to 0.5 will be predicted as neutral and those scores greater than 0.5 will be predicted as pathogenic.

2.4. Determining Protein Stability

Protein stability was determined using I-Mutant 2.0 [20] (https://folding.biofold.org/i-mutant/i-mutant2.0.html (accessed on 3 August 2021)). The website predicts the effect of the amino acid substitution on the stability of the protein as well as the free energy value (DDG value) and the reliability index (RI) value of the amino acids. The RI value reveals the reliability of the prediction, where 0 indicates the least reliable while a 10 indicates the most reliable result. The DDG value measures the energy changes between a folded and unfolded structure. When the DDG value is higher than zero, the mutation is said to be able to increase the protein stability whereas a negative DDG value will decrease the stability of the protein.

2.5. Analyzing Protein Evolutionary Conservation

ConSurf [21] (https://consurf.tau.ac.il/ (accessed on 16 August 2021)) was used for determining the evolutionary conservation of each amino acid residue. It obtains the results by first running a BLAST analysis on the query protein sequence against the UNIREF-90 database and the resulting sequences were aligned using MAFFT which is a multiple sequence alignment (MSA) program that was used for creating MSA of the amino acids sequences. The MSA sequences generated were then used for constructing a phylogenetic tree and calculating the conservation rate based on an empirical Bayesian methodology. The output of this tool was given in the form of a conservation score and the prediction of the amino acid residue on whether it is exposed, buried, functional or structural. The conservation rates were then normalized and grouped into 9 different grades ranging from 1 to 9, with 1 being the most rapidly evolving position, 5 being the position of intermediate rates, and 9 being the most evolutionary conserved positions, meaning that the residue has a slower evolution rate compared to others.

2.6. Biophysical Characteristic Analysis with Align-GVGD

Align-GVGD (A-GVGD) [22] (http://agvgd.hci.utah.edu/ (accessed on 23 August 2021)) tool was used to predict the biological effect of missense substitutions. It predicts the transactivity by measuring the Grantham Variation (GV) which is the degree of biochemical variation among amino acids and Grantham Deviation (GD), which reflects the ‘biochemical distance’ of the mutant amino acid from the observed amino acid at a particular position [23]. The server then uses the GD value to predict whether the nsSNPS is neutral, deleterious, or unclassified. A GD = 0, is predicted as neutral, GV > 61.3 and 0 < GD ≤ 61.3, then the residue is considered as a mutant that has its composition, polarity, and volume fall close to the observed range of variation, and it will also be predicted as neutral by the server. If GV = 0 and GD > 0, this indicates the nsSNP as deleterious.

2.7. Analyzing Protein Interacting Network with Cytoscape

To analyze the protein interacting network, Cytoscape 3.8.2 [24] (https://cytoscape.org/ (accessed on 31 August 2021)) was employed to visualize the biological interaction between molecular complexes, modules, or pathways that possess different biological functions of the FGF4 domain. The STRING (Search Tool for Retrieval of Interacting Genes/Proteins) interaction database [25] (https://string-db.org/ (accessed on 31 August 2021)) was used as the network source for the analysis.

2.8. Prediction of Structural Alteration in FGF4 Domains Using SWISS-Model

SWISS-Model [26] (https://swissmodel.expasy.org/ (accessed on 14 September 2021)) was employed to generate 3D models to predict the structural alterations in the FGF4 domain. The nsSNPs that were identified as functional (exposed and conserved) in ConSurf and present in the FGF domain were selected for analysis.

2.9. Prognosis Analysis

The prognosis was evaluated using the Kaplan-Kaplan-Meier Plotter (http://www.kmplot.com// (accessed on 4 May 2022)) [27]. The survival curves and log-p values were obtained for bladder cancer.

3. Results

3.1. Predicting Deleterious nsSNPs of FGF4

Out of 186 nsSNPs, 68 nsSNPs were predicted as “deleterious” by PROVEAN, 79 “damaging” by SIFT, 88 nsSNPs were predicted as “probably damaging and possibly damaging” by PolyPhen-2, 35 nsSNPs were predicted as “disease” by PhD-SNP and 28 nsSNPs were predicted as “disease” by SNPs&GO. Analysis using PMut showed that all 27 nsSNPs scored > 0.5, suggesting pathogenicity. The results of all the six tools were integrated and 27 nsSNPs predicted as “high-risk” by all tools were proceeded with further analysis (Table 1).

3.2. Predicting Effects of High-Risk nsSNPs on Protein Stability

Protein stability of the nsSNPs was determined by I-Mutant 2.0 and MUPro by comparing free energy. The 27 nsSNPs predicted as “damaging” and “high-risk” were submitted to I-Mutant 2.0 and MUPro for the prediction of change in protein stability due to mutation. Twenty-three nsSNPs were predicted to have decreased protein stability by I-Mutant while MUPro determined 26 nsSNPs to have decreased stability. A total of 22 nsSNPs were predicted by both tools to decrease the protein stability of FGF4 (Table 2).

3.3. Evolutionary Conservation Analysis

The evolutionary conservation of the proteins was determined using the ConSurf tool. Out of the 22high risk nsSNPs associated with decreased protein stability, four nsSNPs; rs1259280329, rs1363460000, rs1413186512, and rs930844659 were predicted to be exposed and functional. These nsSNPs possess a highly conserved sequence, indicating a slower evolutionary rate compared to other residues, therefore involved in important functional roles [28].

3.4. Biophysical Characteristic Analysis

Based on Align-GVGD results, the twentytwo FGF4 nsSNPs fall within class C65 (n = 14), class C55 (n = 3), class C45 (n = 3), class C15 (n = 1), and class C0 (n = 1). Classes of C45, C55 and C65 indicate mutations with functional impact on protein, whereas class C35 (intermediate class) along with classes of C0, C15, and C25 denote mutations with no apparent effect on protein function. The rs1259280329, rs1363460000, rs1413186512, and rs930844659 were predicted to cause functional effect on the protein.

3.5. Construction of Protein-Protein Interaction Network

The interaction network between proteins was constructed by STRING and visualised by Cytoscape (Figure 2). The interaction consisted of 11 nodes and 49 edges. It is predicted that the FGF4 is associated with FGFR2, FGFR1, FGFR3, FGFR4, MAPK1, KRAS, MAPK3, HRAS, NRAS, and KAL1. The degree, average shortest path length, closeness centrality, betweenness centrality, and neighbourhood connectivity of these 10 proteins in interaction with FGF4 protein were also analysed using Cytoscape.

3.6. Prediction of Structural Alteration in FGF4 Domains Using SWISS-Model

Mutant 3D models of rs1259280329, rs1363460000, rs1413186512, and rs930844659 localized in the FGF domain were predicted and generated using SWISS-Model (Figure 3). As hydrophobicity has a major contribution to protein function and structure, the hydrophobicity of wild-type and mutant residues were analyzed in SWISS-Model to investigate their physicochemical properties. The polarity of Q97K rs1259280329 (hydrophilic) remained unchanged throughout the mutation. Changes in polarity were observed for G106V rs1363460000 (neutral to hydrophobic), N164S rs1413186512 (hydrophilic to neutral), and N167S rs930844659 (hydrophilic to neutral).

3.7. Prognosis of FGF4 in Malignancies

A Kaplan-Meier plotter was used to determine the prognostic value of the FGF4 gene by combining gene expression and cancer patient survival. Hazard ratio (HR) with 95% confidence intervals (CI) and logrank p-value were calculated. FGF4 gene showed a hazard ratio (HR) = 8.21 (95% CI, 1.13–9.58) and logrank p-value = 0.012 for bladder carcinoma indicating that the result was statistically significant (the relation between the high expression of the HLA-G gene and more survival rate) (Figure 4).

4. Discussion

Dysregulation of fibroblast growth factors, especially FGF1 and FGF3, have been associated with the risk of bladder cancer [29]. However, to the best of our knowledge, the role of FGF4 in bladder cancer is limited, despite its prominent role in MAPK pathway activation linked with FGF1 and FGF3. This has created an intrigue to investigate the role of FGF4 in bladder cancer and further elucidate its role in prognosis. In this study, we have developed a pipeline (Figure 1) to determine the pathogenic nsSNPs associated with bladder cancer and its prognosis.
In this study, the nsSNPs of FGF4 were subjected to different bioinformatics tools to determine their structural and functional effect on the protein [30]. Damaging nsSNPs were predicted using five different tools, resulting in 27 nsSNPs as “highly damaging”. In order to further narrow down the number of possible pathogenic nsSNPs, P-Mut, I-Mutant, Biophysics analysis, and ConSurf tools were used to predict protein stability, the evolutionary conservation of amino acids, the physical and chemical properties, and changes in protein structure after mutations. This resulted in four “high-risk” nsSNPs; rs1259280329, rs1363460000, rs1413186512, and rs930844659, which were predicted (i) pathogenic by all five predicting tools; (ii) reduced protein stability; and (iii) evolutionary conservation showed that these nsSNPs as highly conserved. This indicates that these nsSNPs with altered protein stability may cause misfolding, degradation, or aberrant conglomeration of proteins [31]. Moreover, we also found that these highly deleterious nsSNPs with high conservation scores could increase the risk of tumorigenesis by inactivating FGF4.
Protein-protein interaction analysis using Cytoscape 3.8.2 indicated FGF4 associated with FGFR2, FGFR1, FGFR3, FGFR4, MAPK1, KRAS, MAPK3, HRAS, NRAS, and KAL1. Based on the shortest average path length and the highest closeness and betweenness centrality, FGFR1, FGFR2, and FGFR3 have the most substantial network interaction with FGF4. FGF4 binds to the tyrosine kinase receptors encoded by FGFR1 (FGFR1c), FGFR2 (FGFR2c), FGFR3 (FGFR3c), and FGFR4 to begin the initiation of signalling. Loss of FGFR protein activity regulation is reported to be oncogenic leading to the overexpression of FGFR protein [32]. Another study has shown that the G388R in the protein kinase domain on FGFR4 has been proven to increase the potential of promoting cancer cells [33].
The Q97K rs1259280329, G106V rs1363460000, N164S rs1413186512, and N167S rs930844659 predicted deleterious are located within the FGF domain, which acts as a binding site for heparan sulphate (HS) proteoglycans to trigger the activation for FGF receptors [34]. The G106V (neutral to hydrophobic), N164S (hydrophilic to neutral), N167S (hydrophilic to neutral), and G190E (neutral to hydrophobic) were predicted to change in polarity. Hydrophobic residues have a significant role in the folding of a protein chain, thus a protein needs to fold accurately to carry out its functions efficiently. Therefore, the change in polarity of G106V could potentially distort the conformation of the resulting protein. This might lead to improper binding of FGF4 to the receptors, causing incomplete signalling transduction and also the production of mutagenic proteins. According to Bellosta (2001) [35], mutated FGF4 tends to interact with FGFR1 and produces ligands with reduced receptor binding potential. Studies have also shown that the gain in the hydrophobicity in G190E nsSNP could disrupt the MAPK pathways that are responsible for cell differentiation, as well as affect the regulations of the states of actin filaments that allow cell migration and division [31,36].
We also evaluated the FGF4 against different types of cancer using the Kaplan-Meier bioinformatics analyses. The results indicated that the FGF4 gene deregulation might affect the overall survival rate of patients with bladder cancer and thus affecting prognosis significance. This finding agrees with Zaharieva and colleagues (2003) who reported gene amplification of FGF4 among bladder cancer patients.

5. Conclusions

In conclusion, a bioinformatics pipeline was developed to efficiently predict the deleterious effect of nsSNPs of FGF4 on bladder cancer. We have also analysed the prognosis of FGF4. We believe that a similar approach can be utilised in annotating and predicting the functional and structural effect of nsSNPs on other diseases.

Author Contributions

Conceptualization, M.S. and W.-S.Y.; methodology, E.C.L., S.W.L. and K.J.T.; formal analysis, E.C.L. and M.S.; writing—original draft preparation, E.C.L.; writing—review and editing, W.H.C. and K.-S.L.; supervision, J.-Y.L. and W.-S.Y.; funding acquisition, J.-Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Program Inisiatif Dana bagi Penggerak Program SDG-WWF (Grant code: MA01MOF1)-UCSI Grant code: IND-FAS-2022/02.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data available in this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lenis, A.T.; Lec, P.M.; Chamie, K.; MSHS, M. Bladder Cancer: A Review. JAMA 2020, 324, 1980–1991. [Google Scholar] [CrossRef] [PubMed]
  2. Xiao, J.-F.; Caliri, A.W.; Duex, J.E.; Theodorescu, D. Targetable Pathways in Advanced Bladder Cancer: FGFR Signaling. Cancers 2021, 13, 4891. [Google Scholar] [CrossRef]
  3. Ng, K.; Stenzl, A.; Sharma, A.; Vasdev, N. Urinary biomarkers in bladder cancer: A review of the current landscape and future directions. Urol. Oncol. Semin. Orig. Investig. 2021, 39, 41–51. [Google Scholar] [CrossRef] [PubMed]
  4. Lang, L.; Teng, Y. Fibroblast Growth Factor Receptor 4 Targeting in Cancer: New Insights into Mechanisms and Therapeutic Strategies. Cells 2019, 8, 31. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Hernández, S.; López-Knowles, E.; Lloreta, J.; Kogevinas, M.; Amorós, A.; Tardón, A.; Carrato, A.; Serra, C.; Malats, N.; Real, F.X. Prospective study of FGFR3 mutations as a prognostic factor in nonmuscle invasive urothelial bladder carcinomas. J. Clin. Oncol. 2006, 24, 3664–3671. [Google Scholar] [CrossRef] [PubMed]
  6. Liu, X.; Zhang, W.; Geng, D.; He, J.; Zhao, Y.; Yu, L. Clinical significance of fibroblast growth factor receptor-3 mutations in bladder cancer: A systematic review and meta-analysis. Genet. Mol. Res. 2014, 13, 1109–1120. [Google Scholar] [CrossRef] [PubMed]
  7. Turner, N.; Grose, R. Fibroblast growth factor signalling: From development to cancer. Nat. Rev. Cancer 2010, 10, 116–129. [Google Scholar] [CrossRef] [PubMed]
  8. Simon, R.; Richter, J.; Wagner, U.; Fijan, A.; Bruderer, J.; Schmid, U.; Ackermann, D.; Maurer, R.; Alund, G.; Knönagel, H. High-throughput tissue microarray analysis of 3p25 (RAF1) and 8p12 (FGFR1) copy number alterations in urinary bladder cancer. Cancer Res. 2001, 61, 4514–4519. [Google Scholar]
  9. Touat, M.; Ileana, E.; Postel-Vinay, S.; André, F.; Soria, J.-C. Targeting FGFR signaling in cancer. Clin. Cancer Res. 2015, 21, 2684–2694. [Google Scholar] [CrossRef] [Green Version]
  10. Loriot, Y.; Necchi, A.; Park, S.H.; Garcia-Donas, J.; Huddart, R.A.; Burgess, E.F.; Fleming, M.T.; Rezazadeh, A.; Mellado, B.; Varlamov, S. Erdafitinib (ERDA.; JNJ-42756493), a Pan-Fibroblast Growth Factor Receptor (FGFR) Inhibitor, in Patients (pts) with Metastatic or Unresectable Urothelial Carcinoma (mUC) and FGFR Alterations (FGFRa): Phase 2 Continuous Versus Intermittent Dosing. J. Clin. Oncol. 2018, 36, 411. [Google Scholar] [CrossRef]
  11. Bange, J.; Prechtl, D.; Cheburkin, Y.; Specht, K.; Harbeck, N.; Schmitt, M.; Knyazeva, T.; Müller, S.; Gärtner, S.; Sures, I. Cancer progression and tumor cell motility are associated with the FGFR4 Arg388 allele. Cancer Res. 2002, 62, 840–847. [Google Scholar] [PubMed]
  12. Spinola, M.; Leoni, V.P.; Tanuma, J.-i.; Pettinicchio, A.; Frattini, M.; Signoroni, S.; Agresti, R.; Giovanazzi, R.; Pilotti, S.; Bertario, L. FGFR4 Gly388Arg polymorphism and prognosis of breast and colorectal cancer. Oncol. Rep. 2005, 14, 415–419. [Google Scholar] [CrossRef]
  13. Peddle-McIntyre, C. Poor prognosis cancers–The role of exercise. J. Sci. Med. Sport 2018, 21, S74. [Google Scholar] [CrossRef]
  14. Choi, Y.; Chan, A.P. PROVEAN web server: A tool to predict the functional effect of amino acid substitutions and indels. Bioinformatics 2015, 31, 2745–2747. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Hu, J.; Ng, P.C. Predicting the effects of frameshifting indels. Genome Biol. 2012, 13, R9. [Google Scholar] [CrossRef] [Green Version]
  16. Kumar, P.; Henikoff, S.; Ng, P.C. Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat. Protoc. 2009, 4, 1073–1081. [Google Scholar] [CrossRef]
  17. Adzhubei, I.A.; Schmidt, S.; Peshkin, L.; Ramensky, V.E.; Gerasimova, A.; Bork, P.; Kondrashov, A.S.; Sunyaev, S.R. A method and server for predicting damaging missense mutations. Nat. Methods 2010, 7, 248–249. [Google Scholar] [CrossRef] [Green Version]
  18. Capriotti, E.; Calabrese, R.; Fariselli, P.; Martelli, P.L.; Altman, R.B.; Casadio, R. WS-SNPs&GO: A web server for predicting the deleterious effect of human protein variants using functional annotation. BMC Genom. 2013, 14, S6. [Google Scholar] [CrossRef] [Green Version]
  19. López-Ferrando, V.; Gazzo, A.; de la Cruz, X.; Orozco, M.; Gelpí, J.L. PMut: A web-based tool for the annotation of pathological variants on proteins, 2017 update. Nucleic Acids Res. 2017, 45, W222–W228. [Google Scholar] [CrossRef]
  20. Capriotti, E.; Fariselli, P.; Casadio, R. I-Mutant2.0: Predicting stability changes upon mutation from the protein sequence or structure. Nucleic Acids Res. 2005, 33, W306–W310. [Google Scholar] [CrossRef] [Green Version]
  21. Ashkenazy, H.; Abadi, S.; Martz, E.; Chay, O.; Mayrose, I.; Pupko, T.; Ben-Tal, N. ConSurf 2016: An improved methodology to estimate and visualize evolutionary conservation in macromolecules. Nucleic Acids Res. 2016, 44, W344–W350. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Tavtigian, S.V.; Deffenbaugh, A.M.; Yin, L.; Judkins, T.; Scholl, T.; Samollow, P.B.; de Silva, D.; Zharkikh, A.; Thomas, A. Comprehensive statistical study of 452 BRCA1 missense substitutions with classification of eight recurrent substitutions as neutral. J. Med. Genet. 2006, 43, 295–305. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Mathe, E.; Olivier, M.; Kato, S.; Ishioka, C.; Hainaut, P.; Tavtigian, S.V. Computational approaches for predicting the biological effect of p53 missense mutations: A comparison of three sequence analysis based methods. Nucleic Acids Res. 2006, 34, 1317–1325. [Google Scholar] [CrossRef]
  24. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef] [PubMed]
  25. Szklarczyk, D.; Gable, A.L.; Nastou, K.C.; Lyon, D.; Kirsch, R.; Pyysalo, S.; Doncheva, N.T.; Legeay, M.; Fang, T.; Bork, P. The STRING database in 2021: Customizable protein–protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 2021, 49, D605–D612. [Google Scholar] [CrossRef] [PubMed]
  26. Waterhouse, A.; Bertoni, M.; Bienert, S.; Studer, G.; Tauriello, G.; Gumienny, R.; Heer, F.T.; de Beer, T.A.P.; Rempfer, C.; Bordoli, L. SWISS-MODEL: Homology modelling of protein structures and complexes. Nucleic Acids Res. 2018, 46, W296–W303. [Google Scholar] [CrossRef] [Green Version]
  27. Lánczky, A.; Győrffy, B. Web-based survival analysis tool tailored for medical research (KMplot): Development and implementation. J. Med. Internet Res. 2021, 23, e27633. [Google Scholar] [CrossRef]
  28. Ben Chorin, A.; Masrati, G.; Kessel, A.; Narunsky, A.; Sprinzak, J.; Lahav, S.; Ashkenazy, H.; Ben-Tal, N. ConSurf-DB: An accessible repository for the evolutionary conservation patterns of the majority of PDB proteins. Protein Sci. 2020, 29, 258–267. [Google Scholar] [CrossRef]
  29. Subbiah, V.; Iannotti, N.; Gutierrez, M.; Smith, D.; Féliz, L.; Lihou, C.; Tian, C.; Silverman, I.; Ji, T.; Saleh, M. FIGHT-101, a first-in-human study of potent and selective FGFR 1–3 inhibitor pemigatinib in pan-cancer patients with FGF/FGFR alterations and advanced malignancies. Ann. Oncol. 2022, 33, 522–533. [Google Scholar] [CrossRef]
  30. Lim, S.W.; Tan, K.J.; Azuraidi, O.M.; Sathiya, M.; Lim, E.C.; Lai, K.S.; Yap, W.-S.; Afizan, N.A.R.N.M. Functional and structural analysis of non-synonymous single nucleotide polymorphisms (nsSNPs) in the MYB oncoproteins associated with human cancer. Sci. Rep. 2021, 11, 24206. [Google Scholar] [CrossRef]
  31. Cooper, G.M.; Hausman, R. A molecular approach. In The Cell, 2nd ed.; Sinauer Associates: Sunderland, MA, USA, 2000. [Google Scholar]
  32. Haugsten, E.M.; Wiedlocha, A.; Olsnes, S.; Wesche, J. Roles of fibroblast growth factor receptors in carcinogenesis. Mol. Cancer Res. 2010, 8, 1439–1452. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Kim, J.H.; Jeong, S.Y.; Jang, H.J.; Park, S.T.; Kim, H.S. FGFR4 Gly388Arg Polymorphism Reveals a Poor Prognosis, Especially in Asian Cancer Patients: A Meta-Analysis. Front. Oncol. 2021, 11, 762528. [Google Scholar] [CrossRef] [PubMed]
  34. Anjay, A. National Center for Biotechnology Information (ncBi). 2012. Available online: https://www.ncbi.nlm.nih.gov/ (accessed on 16 September 2021).
  35. Bellosta, P.; Iwahori, A.; Plotnikov, A.N.; Eliseenkova, A.V.; Basilico, C.; Mohammadi, M. Identification of receptor and heparin binding sites in fibroblast growth factor 4 by structure-based mutagenesis. Mol. Cell. Biol. 2001, 21, 5946–5957. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Zhang, W.; Liu, H.T. MAPK signal pathways in the regulation of cell proliferation in mammalian cells. Cell Res. 2002, 12, 9–18. [Google Scholar] [CrossRef]
Figure 1. Schematic workflow for predicting high-risk deleterious nsSNPs and cancer prognosis.
Figure 1. Schematic workflow for predicting high-risk deleterious nsSNPs and cancer prognosis.
Life 12 01018 g001
Figure 2. Protein–protein interaction network of FGF4 with 10 partners.
Figure 2. Protein–protein interaction network of FGF4 with 10 partners.
Life 12 01018 g002
Figure 3. Comparison of wild-type FGF4 protein structure with its mutant forms. (A) Superimposed structures of wild-type FGF4 protein with mutant protein Q97K, Glutamine into a Lysine at position 97, (B) Superimposed structures of wild-type FGF4 protein with mutant protein G106V, Glycine into a Valine at position 106. (C) Superimposed structures of wild-type FGF4 protein with mutant protein N164S, Asparagine into a Serine at position 164, and (D) Superimposed structures of wild-type FGF4 protein with mutant protein N167S, Asparagine into a Serine at position 167.
Figure 3. Comparison of wild-type FGF4 protein structure with its mutant forms. (A) Superimposed structures of wild-type FGF4 protein with mutant protein Q97K, Glutamine into a Lysine at position 97, (B) Superimposed structures of wild-type FGF4 protein with mutant protein G106V, Glycine into a Valine at position 106. (C) Superimposed structures of wild-type FGF4 protein with mutant protein N164S, Asparagine into a Serine at position 164, and (D) Superimposed structures of wild-type FGF4 protein with mutant protein N167S, Asparagine into a Serine at position 167.
Life 12 01018 g003
Figure 4. Kaplan-Meier plot showing correlation of deregulation of FGF4 gene and overall survival rate of (A) bladder carcinoma and (B) lung adenocarcinoma patients.
Figure 4. Kaplan-Meier plot showing correlation of deregulation of FGF4 gene and overall survival rate of (A) bladder carcinoma and (B) lung adenocarcinoma patients.
Life 12 01018 g004
Table 1. High risk nsSNPs predicted by the six different in silico tools.
Table 1. High risk nsSNPs predicted by the six different in silico tools.
nsSNPs IDAA
Change
PROVEANSIFTPolyphen-2 PhD-SNPSNPs&GOPMut
Pred (Cut Off = −2.5)ScPred (Cut Off = 0.05)ScPredScPred (Cut Off = 0.5)RIProbPred (Cut Off = 0.5)RIProbPred (Cut Off = 0.5)Sc
rs1383383982D75VDel−3.9Dmg0.002Pro.dmg0.993Dis10.610Dis20.542Dis0.81
rs922987433D75YDel−3.92Dmg0.002Pro.dmg1Dis10.680Dis40.538Dis0.81
rs760825703R85WDel−5Dmg0.016Pro.dmg0.988Dis50.889Dis80.769Dis0.77
rs1266598072G91DDel−5.56Dmg0.001Pro.dmg1Dis10.878Dis80.554Dis0.79
rs1194178508G93DDel−5.74Dmg0Pro.dmg1Dis60.917Dis80.825Dis0.83
rs1250040489G93RDel−6.53Dmg0Pro.dmg1Dis60.898Dis80.803Dis0.79
rs775542907F94SDel−6.44Dmg0.001Pro.dmg0.980Dis70.888Dis80.853Dis0.83
rs1259280329Q97KDel−3.5Dmg0.001Pro.dmg0.998Dis60.854Dis70.791Dis0.76
rs1469284144I104NDel−6.05Dmg0Pro.dmg1Dis70.801Dis60.871Dis0.89
rs1363460000G106VDel−8.14Dmg0Pro.dmg1Dis50.739Dis50.755Dis0.86
rs1432374845L118RDel−4.61Dmg0.001Pro.dmg0.996Dis10.714Dis40.543Dis0.59
rs1245810774G124VDel−8.58Dmg0Pro.dmg1Dis50.685Dis60.639Dis0.91
rs539419605G124SDel−5.74Dmg0Pro.dmg1Dis30.793Dis40.756Dis0.77
rs374997743I128FDel−2.91Dmg0.004Pro.dmg0.96Dis10.862Dis70.534Dis0.79
rs979866825G130SDel−5.68Dmg0Pro.dmg1Dis20.688Dis40.578Dis0.74
rs966807008S133IDel−5.32Dmg0Pro.dmg1Dis40.781Dis60.712Dis0.88
rs781699363A138TDel−3.51Dmg0.002Pro.dmg1Dis40.696Dis50.52Dis0.87
rs757487910M139LDel−2.86Dmg0Pro.dmg0.992Dis40.764Dis50.697Dis0.83
rs1283278927L145PDel−6.29Dmg0Pro.dmg1Dis50.836Dis70.741Dis0.78
rs764426431Y146CDel−7.41Dmg0.001Pro.dmg1Dis30.760Dis50.651Dis0.85
rs779058257E154GDel−5.75Dmg0.001Pro.dmg1Dis60.773Dis50.785Dis0.82
rs756008893C155SDel−9.47Dmg0.006Pro.dmg1Dis60.864Dis70.782Dis0.9
rs1413186512N164SDel−4.73Dmg0Pro.dmg0.999Dis10.732Dis50.552Dis0.86
rs986306143Y166HDel−4.73Dmg0Pro.dmg1Dis20.755Dis50.583Dis0.77
rs930844659N167SDel−4.73Dmg0Pro.dmg1Dis20.562Dis10.581Dis0.83
rs1182350769S171YDel−5.31Dmg0Pro.dmg1Dis10.794Dis60.53Dis0.78
rs866953016G190EDel−7.27Dmg0Pro.dmg1Dis40.680Dis40.705Dis0.9
Table 2. Overall results returned by I-Mutant 2.0, MUPro, ConSurf, and Align-GVGD.
Table 2. Overall results returned by I-Mutant 2.0, MUPro, ConSurf, and Align-GVGD.
nsSNPs IDAA ChangeI-Mutant 2.0MUProConSurfAlign-GVGD
StabRIDDGStabDDGPredScGVGDPred
rs1383383982D75VDecrease1−1.23Increase0.15056579Ex40152.01Class C65
rs922987433D75YDecrease1−0.99Decrease−0.24113434 Ex40159.94Class C65
rs760825703R85WDecrease6−0.34Decrease−0.69402177 Ex70101.29Class C65
rs1266598072G91DDecrease6−0.72Decrease−0.36065593Bu7093.77Class C65
rs1194178508G93DDecrease8−1.16Decrease−0.17477258Ex7093.77Class C65
rs1250040489G93RDecrease7−0.96Decrease−0.28775219 Ex70125.13Class C65
rs775542907F94SDecrease8−2.77Decrease−1.5641414 Ex40154.81Class C65
rs1259280329Q97KDecrease4−0.58Decrease−0.62838088 Ex & Fn8053.23Class C45
rs1469284144I104NDecrease5−0.54Decrease−1.1107512Bu80148.91Class C65
rs1363460000G106VDecrease3−1.32Decrease−0.73393095 Ex & Fn90108.79Class C65
rs1432374845L118RDecrease8−2.07Decrease−1.6708666 Bu70101.88Class C65
rs1245810774G124VIncrease1−0.16Decrease−0.49592065 Ex70108.79Class C65
rs539419605G124SDecrease7−0.95Decrease−1.2395925Ex7055.27Class C55
rs374997743I128FDecrease7−1.88Decrease−1.3022244 Bu8021.28Class C15
rs979866825G130SDecrease7−1.02Decrease−0.90820561 Bu7055.27Class C55
rs966807008S133IIncrease1−0.37Decrease−0.05605553 Bu70141.8Class C65
rs781699363A138TDecrease7−1.15Decrease−1.0148323 Bu8058.02Class C55
rs757487910M139LDecrease6−0.35Decrease−0.67563291 Bu & St9014.3Class C0
rs1283278927L145PDecrease8−1.92Decrease−1.9056486Bu8097.78Class C65
rs764426431Y146CDecrease10.15Decrease−0.46046358Ex70193.72Class C65
rs779058257E154GDecrease4−0.61Decrease−0.9707573Ex7097.85Class C65
rs756008893C155SDecrease6−1.37Decrease−1.4335656 Bu & St90111.67Class C65
rs1413186512N164SDecrease6−0.69Decrease−0.86604139Ex & Fn9046.24Class C45
rs986306143Y166HDecrease8−1.68Decrease−1.3498789 Bu & St9083.33Class C65
rs930844659N167SDecrease5−1.06Decrease−0.72801042 Ex & Fn9046.24Class C45
rs1182350769S171YIncrease40.96Decrease−1.0600288 Ex & Fn90143.11Class C65
rs866953016G190EIncrease30.64Decrease−0.34251717Ex & Fn9097.85Class C65
AA: Amino acid; Pred: Prediction; Sc: Score; RI: Reliability Index; Stab: Stability prediction; DDG: Free energy change; RMSD: Root mean square deviation; Dis: Disease; Ex: Exposed; Fn: Functional; Bu: Buried; St: Structural; GV: Grantham Variation; GD: Grantham Deviation.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lim, E.C.; Lim, S.W.; Tan, K.J.; Sathiya, M.; Cheng, W.H.; Lai, K.-S.; Loh, J.-Y.; Yap, W.-S. In-Silico Analysis of Deleterious SNPs of FGF4 Gene and Their Impacts on Protein Structure, Function and Bladder Cancer Prognosis. Life 2022, 12, 1018. https://doi.org/10.3390/life12071018

AMA Style

Lim EC, Lim SW, Tan KJ, Sathiya M, Cheng WH, Lai K-S, Loh J-Y, Yap W-S. In-Silico Analysis of Deleterious SNPs of FGF4 Gene and Their Impacts on Protein Structure, Function and Bladder Cancer Prognosis. Life. 2022; 12(7):1018. https://doi.org/10.3390/life12071018

Chicago/Turabian Style

Lim, Ee Chen, Shu Wen Lim, Kenneth JunKai Tan, Maran Sathiya, Wan Hee Cheng, Kok-Song Lai, Jiun-Yan Loh, and Wai-Sum Yap. 2022. "In-Silico Analysis of Deleterious SNPs of FGF4 Gene and Their Impacts on Protein Structure, Function and Bladder Cancer Prognosis" Life 12, no. 7: 1018. https://doi.org/10.3390/life12071018

APA Style

Lim, E. C., Lim, S. W., Tan, K. J., Sathiya, M., Cheng, W. H., Lai, K. -S., Loh, J. -Y., & Yap, W. -S. (2022). In-Silico Analysis of Deleterious SNPs of FGF4 Gene and Their Impacts on Protein Structure, Function and Bladder Cancer Prognosis. Life, 12(7), 1018. https://doi.org/10.3390/life12071018

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