1. Introduction
Signal Transducer and Activator of Transcription (STAT) proteins are a family of transcription factors latently present in the cytoplasm and participate in a variety of cellular events following cytokines and growth factors signaling [1,2]. STAT proteins are involved in intracellular signaling downstream of the type I and type II cytokine receptors. Upon activation, translocation to the nucleus, binding to their specific promoter regions of target genes and regulation of their transcription subsequently takes place [3,4]. Seven proteins have been identified (STAT1, -2, -3, -4, -5a, -5b, and -6) and share a common structure consisting of an SH2 domain that mediates STAT interactions through homo- or heterodimers, a coiled-coil domain, which is important for dimer nuclear localization, a DNA-binding domain, which leads to target gene transcription, and a transactivation domain [5,6].
The Signal Transducer and Activator of Transcription 1 (STAT1) gene is composed of 25 exons and 7 domains, located on chromosome 2q32.2 [7,8,9]. STAT1 is an essential mediator of the JAK-STAT signaling pathway in response to interferons [8,10,11,12]. It plays a crucial role in the biological immune response against intracellular mycobacterial infection as well as viral infections [8,13,14]. Upon type I IFN-gamma (IFN-γ) binding to cell surface receptors, there is a signaling pathway through protein kinases then activation of Jak kinases (TYK2 and JAK1) with tyrosine phosphorylation of STAT1, dimerization of phosphorylated STAT1, and association with ISGF3G/IRF-9 forming ISGF3 transcription factor [15]. ISGF3 enters the nucleus and binds to the IFN-stimulated response element (ISRE) to activate the transcription of IFN-stimulated genes (ISG), which bring the cell into an antiviral state [16]. Moreover, in response to type II IFN, STAT1 is tyrosine- and serine-phosphorylated; it then forms a homodimer termed IFN-gamma-activated factor (GAF) [17] that migrates into the nucleus and binds to the IFN-gamma-activated sequence (GAS) to drive the expression of the target genes, inducing a cellular antiviral state [18].
Genetic variants within STAT1 gene lead to loss-of-function (LOF) and gain-of-function (GOF) phenotypes, with a wide range of clinical presentations, including autoimmunity and life-threatening mycobacterial, severe viral, and bacterial infections [19,20,21]. STAT1 amorphic alleles cause severe viral and bacterial infections, while hypomorphic alleles cause mild disseminated mycobacterial disease [22]. Moreover, hypermorphic mutations are responsible for a variety of clinical presentations such as chronic mucocutaneous candidiasis (CMC), arterial aneurysms, autoimmunity, and squamous cell cancers [23]. STAT1 gain-of-function (GOF) mutation, mostly located at coiled-coil (CCD) and DNA-binding domains (DBD) causing hyper-phosphorylation of STAT1 protein, thus enhanced STAT1-dependent responses to interferons (IFNs) and IL-27, with sequential impairment of Th17 cell development [24,25,26]. GOF mutation is associated with chronic mucocutaneous candidiasis [10,27,28], while patients with LOF mutations display an increased susceptibility to intracellular bacteria, including a Mendelian susceptibility to mycobacterial disease (MSMD) [10,22].
Single-nucleotide polymorphisms (SNPs) constitute a common form of genetic variation in humans [29]. The nonsynonymous SNPs (nsSNPs) cause alteration in the amino acid residues because of variation in the sequence of DNA at a single position of a nucleotide (A, T, C, or G), which contributes to the functional diversity of the related proteins [30,31,32].
Recently, bioinformatics tools have played a significant role in the prediction of damaging SNPs and their relationship with diseases [33]. The influence of STAT1 nsSNPs on protein structure and function has not been thoroughly investigated, despite their potential importance; this indicates a substantial scientific gap. Nonetheless, limited published articles have systematically examined STAT1 SNPs by bioinformatics approaches.
The objective of this study is to define the structural and functional characterization of the most pathogenic variations of the STAT1 gene. We performed a comprehensive STAT1-SNPs analysis using bioinformatics prediction tools combined with artificial intelligence models to identify the pathogenic and deleterious SNPs, providing novel insights into their involvement in immune dysregulation and establishing a foundation for subsequent functional and clinical research.
2. Materials and Method
An overview of the complete methodological approach is shown in Figure 1.
2.1. Data Retrieval
We gathered the data for the human STAT1 gene from the National Center for Biological Information (NCBI) website (
2.2. Phenotype Prediction of Deleterious ns SNPs
We predicted the deleterious nsSNPs by using eight different tools. Sorting Intolerant from Tolerant (SIFT) (
Polyphen-2 (
Provean (
SNAP2 (
PHD-SNP (
SNP and GO (
P-Mut is a web-based tool for the annotation of pathological variants on proteins. It allows fast and accurate prediction of the pathological properties of single-point amino acid mutations based on the use of a neural network. It is available at (
Protein Analysis through Evolutionary Relationships (PANTHER) (
2.3. Predicting Functional and Structural Effects of the nsSNP
MutPred v1.2 (
2.4. Protein Stability Analysis of Predicted STAT1 nsSNPs
I-Mutant 3.0 is available at (
MUpro, a group of machine learning methods, predicts the effects of single amino acid substitutions on protein stability [45]. It uses both support vector machines and neural networks; the output is either increased or decreased stability [45]. MUpro also interprets the result based on Gibbs free energy (ΔΔG), with a confidence score between −1 and 11. It is available at
DDMUT (
2.5. Prediction of Missense Variant Pathogenicity
Alpha Missense is an adaptation of alphafold fine-tuned on human and primate variant population frequency databases to predict missense variant pathogenicity. It works by combining structural context and evolutionary conservation. This model achieves state-of-the-art results across a wide range of genetic and experimental benchmarks, all without explicitly training on such data [47].
2.6. Three-Dimensional Structure Prediction and Visualization
We predicted the 3D structure using an artificial intelligence system, AlphaFold (
UCSF ChimeraX 1.9 is a robust application that enables interactive viewing and analysis of various molecular structures and related data, including density maps, sequence alignments, and supramolecular assemblies [49]. It allows the mapping and visualization of amino acid substitutions. Chimera X is available at
2.7. Phenotypic Effects Prediction
Project Hope (version 1.0) is an online web server used to analyze the structural and conformational variations that have resulted from single amino acid substitutions [50]. We uploaded STAT1 protein sequence, wild-type amino acids, and mutants. The results provided describe the change in the physiochemical properties of the amino acid in the given SNPs. It is available at (
DDMUT can also detect changes in the biological interactions between wild-type amino acids and neighborhood residues in comparison with mutant residues [46].
2.8. Conservational Analysis and Surface Accessibility Prediction of STAT1
The ConSurf bioinformatics tool (
2.9. Identification of nsSNPs in STAT1 Protein Domains
We submitted the FASTA sequence of the STAT1 protein to the InterPro server (
2.10. Prediction of Protein–Protein Interactions
A precomputed database, STRING (
3. Results
3.1. Distribution of STAT1 Gene SNP Datasets
The total number of SNPs was 10,989. There were 888 frame shift mutations; 480 SNPs located in the coding region, of which 247 were nsSNPs and 233 were synonymous SNPs (sSNPs), while 9.621 SNPs were in noncoding regions, of which 375 occurred in the 3′UTR, 131 in the 5′UTR region, and the rest (9115) were in the intronic region, as shown in Figure 2. We chose nonsynonymous coding SNPs for our investigation.
3.2. Identification of Deleterious Missense Mutation
All 247 nsSNPs were retrieved and subjected to pathogenicity prediction web servers. Sixty-four nsSNPs were found to be deleterious by SIFT and were further subjected to crosschecking by using three different tools (Poly-Phen-2, PROVEAN, and SNAP2).
The shortlisted 33 nsSNPs passed the first four tools, presented in Table 1, then were submitted to another set of four tools: P Mut, PhD-SNP, SNPs and GO, and PANTHER. In total, 29 SNPs out of the 33 predicted by the first set of tools are disease-causing by P mut, 21 out of 33 are disease-causing by Panther, 20 are disease-causing by PhD-SNP, and 14 out of 33 by SNP and GO. A final nine nsSNPs passed all eight tools shown in Table 2. We further analyzed the final set of SNPs for the functional and structural modifications.
3.3. MutPred Prediction for Functional and Structural Modifications
We submitted the shortlisted nine nsSNPs to the MutPred server, along with the resultant probability scores and their p values in Table 3. The structural and functional alterations predicted include loss of disorder, catalytic residue, glycosylation, gain of phosphorylation, solvent accessibility, ubiquitination, and molecular recognition features (MoRF) binding. According to these predictions, several nsSNPs might be the reason behind any possible structural and functional modifications of STAT1 protein.
3.4. Prediction of Change in STAT1 Stability Due to Mutation
We used I mutant, MUpro, and DDMUT servers to predict the effect of the nsSNPs on protein stability. The result revealed that six variants destabilized the STAT 1 protein, namely (I648T) rs759271255, (V642D) rs752542806, (R602W) rs 1209841496, (L600P) rs137852678, (I578N) rs767475430, and (W504C) rs916580554. The results are presented in Table 4.
3.5. Pathogenicity Prediction Results
We analyzed STAT1 nsSNP by Alpha-Missense, and we found that all the pathogenic nsSNPs that were predicted by the previous tools were also classified as pathogenic in Alpha-Missense, presented in Table 5. The heat map represented the mutations in STAT1, as shown in Figure 3.
3.6. The Conservational Status and Surface Accessibility Analysis of STAT1 Protein
Highly conserved residues are most likely to be involved in proteins’ structural integrity and functions. We evaluated the conservational profile for the STAT1 protein. The ConSurf algorithm represented the structural and functional conservation levels of all the amino acid residues of the STAT1 protein. Four SNPs (I648T, L600P, W504C, and I578N) are predicted to be located in a conserved region. L600P and I578N are predicted to be structural residues (highly conserved and buried). V642D is predicted as buried, and R602W is predicted as a functional residue (highly conserved and exposed), presented in Table 6.
3.7. Three-Dimensional Structure Prediction by AlphaFold and SNP Visualization by ChimeraX
An individual residue confidence score (pLDDT) between 0 and 100 is generated by the AlphaFold algorithm. Alphafold produces a per residue confidence score (pLDDT) 1–100. Regions with low pLDDT may be unstructured in isolation. The majority of the 3D structural region corresponds to alpha-helical domains and has extremely high confidence (pLDDT > 90). The remaining components of the model are depicted as unresolved loops with low (70 > pLDDT > 50) and extremely low (pLDDT > 50) scores, as in Figure 4.
We used ChimeraX to visualize the 3D structures of the wild-type amino acids in blue and the mutant residues in red, as shown in Figure 5.
3.8. The Physical Outcome of Predicted SNPs
We examined the impact of the generated damaging SNPs on the three-dimensional structure of STAT1 using the HOPE server. The server predicted that all the mutated amino acids were different in size; one had a different charge, and six had different hydrophobicity. The results are in Table 7.
Loss of the interactions between the wild-type amino acid and other amino acids in the protein and/or development of new interactions or bonds between the mutant residue of the protein and the other amino acids in the protein were predicted by DDMUT, as presented in Figure 6.
3.9. Domain Identification of the STAT1 Protein by the InterPro Server
The InterPro tool predicted the domain regions of the STAT1 protein. The STAT1, SH2 domain (a phosphotyrosine binding pocket) at position (557–707), STAT transcription factor, DNA binding domain at (323–458), and STAT1_TAZ2-binding domain (715–739) are conserved sites. Src homology 2 (SH2) domain profile (573–670), SH2 domain (578–638), STAT1 transcription factor, all alpha domain (144–305), and STAT transcription factor protein interaction (2–12) are as in Table 8.
3.10. STAT1–Protein Interaction
Analysis of protein–protein interaction using the STRING network revealed that STAT1 interacts with 10 proteins, which include other proteins of the same STAT family (STAT2 and STAT3), proteins of the JAK family (JAK1 and JAK2), IFR1, IFR9, IFNGR1, CREBBP, KBNA1, and PIAS1, as presented in Figure 7.
4. Discussion
We evaluated the functional and pathogenic sequences of missense SNPs of the human STAT1 gene, utilizing 12 diverse in silico prediction tools (SIFT, PolyPhen2, PROVEAN, PANTHER, P MUT, PhD-SNP, SNPs&GO, SNAP2, and MutPred2). In silico prediction analysis identified six variants (I648T, V642D, R602W, L600P, I578N, and W504C) considered pathogenic and deleterious. These mutations have a major impact on the protein’s physicochemical characteristics, such as its size and charge hydrophobicity, which ultimately affect the protein’s stability and function and may have an impact on disease. Furthermore, we assessed the effect of missense SNPs on the stability of the STAT1 structure utilizing three stability prediction algorithms: I-Mutant3, MUpro, and DDMUT. All the variants revealed a reduction in stability by the three stability prediction tools (I-Mutant3, MUpro, and DDMUT). In general, we assumed that all missense SNPs in the STAT1 gene were highly unstable in their protein structures, so they were selected for further structural bioinformatics analysis utilizing various tools to explore the consequences of tentatively destructive missense SNPs on STAT1 protein function. To evaluate the conservation profile, we used the ConSurf algorithm to represent the structural and functional conservation levels of all the amino acid residues of STAT1 protein. The ConSurf analysis revealed that the variant in position 602 is a functional residue in a highly conserved and exposed position. Structural residues in highly conserved and buried positions were identified in positions 600 and 578. The identified variants were found in a highly conserved region; this finding suggests that they might be involved in modifications of molecular mechanisms such as bond gain or loss.
STAT1 GOF mutations with CMC were first described in 2001 and 2011, respectively; later, studies confirmed that STAT1-GOF mutations cause immunodeficiency and immune dysregulation, with a wide clinical spectrum [54].
Among the six SNPs identified linked to STAT1 gene mutations in this study, some of these SNPs have been associated with diseases in previous studies, while others were projected to be so in this study using various computational tools. Population genetics and clinical studies are crucial to verifying the results of such research, even though utilizing computational techniques to analyze the impact of the SNPs may aid in identifying disease-related SNPs.
One mutation, namely L600P, has already been previously reported as a mutation in the STAT1 gene in an infant who died of a viral-like illness associated with complete STAT1 deficiency and carried a homozygous nucleotide substitution (T→C) in exon 20, resulting in the substitution of a proline for a leucine at amino acid position 600 [55]. This mutation was found to be pathogenic using all the bioinformatics tools. I648T, V642D, R602W, I578N, and W504C were not reported previously.
Three mutations, namely L706S (rsRCV000009610), Q463H (VAR_065817), and E320Q (VAR_065816), have been reported as mutations in the STAT1 gene. The two previously reported types of autosomal-dominant (AD) Mendelian susceptibility to mycobacterial disease (AD-MSMD) causing STAT1 mutations are located in the tail segment domain (p.L706S) or in the DNA-binding domain (p.E320Q and p.Q463H) [56]. These mutations were not available in the dbSNP database. Two other SNPs (K637E) and (K673R) affecting the SH2 domain, which has been previously reported in two cases with AD-STAT1 deficiency in two unrelated patients from Japan and Saudi Arabia, were also not available in the dbSNPs database at the time of the analysis [56].
Two mutations linked to chronic mucocutaneous candidiasis are (T437I) and (Q271P). Q271P occurred within a specific pocket of the STAT1 coiled-coil domain, near residues essential for dephosphorylation, and was identified in a German patient who presented at 1 year of age with autosomal dominant chronic mucocutaneous candidiasis, showed signs of thyroid autoimmunity, and died at age 41 from squamous cell carcinoma [57,58]. These mutations were not available in the dbSNPs database.
The A267V variant in STAT1 has been reported in >10 individuals with chronic mucocutaneous candidiasis (CMC) and segregated with disease in 16 individuals from nine families [59]. This mutation was not present in the dbSNP database.
Interestingly, nsSNPs in the STAT1 gene will ultimately affect and may disturb the normal function of other interacting genes. As our study was in detail, it provides all the information and analysis needed for the identification of the most damaging nsSNPs. Like ours, there are certain limitations in every study. Utilizing in silico technologies is now a crucial method for identifying disease-related SNPs. In this study, the STAT1 gene underwent a thorough analysis utilizing 18 genetics analysis tools (10 computational tools and 8 AI-based methods) to determine the impact of nsSNPs on the protein’s structure and function.
Our study is based on computer tools and web servers, which are based on mathematical and statistical algorithms. Therefore, to confirm these results, experimental investigation is necessary.
5. Conclusions
Our study provides an insight about nsSNPs of the STAT1 gene, its protein 3D structure, and its interactions with other genes, which might be helpful in future studies of STAT1 in order to better understand its role in immunity and all related diseases.
Conceptualization, E.K., L.A.K. and M.A.; Data curation, E.K., A.A. and M.A.; Funding acquisition, E.K.; Investigation, E.K., M.A., L.A.K. and A.A.; Methodology, E.K. and M.A.; Software, E.K. and L.A.K.; Supervision, A.A.; Validation, E.K. and M.A.; Writing—original draft, E.K., L.A.K., M.A. and A.A.; Writing—review and editing, E.K., L.A.K., M.A. and A.A. All authors have read and agreed to the published version of the manuscript.
Not applicable.
Not applicable.
The original contributions presented in this study are included in the article/
The authors declare no conflicts of interest.
Footnotes
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.
Figure 3. Heat map generated by Alpha-Missense shows the variations in STAT1 gene.
Figure 5. Effect of the six most deleterious nsSNPs on the STAT1 protein structure. ChimeraX software was used to visualize the 3D structure of the wild-type (blue), mutant residues (red) and gold ion (yellow).
Figure 5. Effect of the six most deleterious nsSNPs on the STAT1 protein structure. ChimeraX software was used to visualize the 3D structure of the wild-type (blue), mutant residues (red) and gold ion (yellow).
Figure 5. Effect of the six most deleterious nsSNPs on the STAT1 protein structure. ChimeraX software was used to visualize the 3D structure of the wild-type (blue), mutant residues (red) and gold ion (yellow).
Figure 5. Effect of the six most deleterious nsSNPs on the STAT1 protein structure. ChimeraX software was used to visualize the 3D structure of the wild-type (blue), mutant residues (red) and gold ion (yellow).
Figure 6. Difference in ionic interactions between the wild-type (A) and mutant residues (B).
Figure 6. Difference in ionic interactions between the wild-type (A) and mutant residues (B).
Figure 6. Difference in ionic interactions between the wild-type (A) and mutant residues (B).
Figure 6. Difference in ionic interactions between the wild-type (A) and mutant residues (B).
List of nsSNPs that were predicted to have deleterious effect by SIFT, PolyPhen-2, Provean, and SNAP.
SNP ID | Amino Acid Change | SIFT | Poly-Phen-2 | Provean | SNAP2 | |||||
---|---|---|---|---|---|---|---|---|---|---|
Prediction | TI | Effect | Score | Effect | Score | Prediction | Score | |||
1 | rs1173266737 | P728A | Deleterious | 0 | PB | 0.974 | Deleterious | −3.14 | E | 15 |
2 | rs1374373369 | D674V | Deleterious | 0.01 | PB | 0.989 | Deleterious | −7.327 | E | 47 |
3 | rs771679419 | Y668F | Deleterious | 0 | PS | 0.609 | Deleterious | −3.583 | E | 67 |
4 | rs759271255 | I648T | Deleterious | 0.01 | PB | 1 | Deleterious | −4.197 | E | 56 |
5 | rs752542806 | V642D | Deleterious | 0 | PB | 0.994 | Deleterious | −5.149 | E | 60 |
6 | rs1387961263 | L639F | Deleterious | 0 | PB | 1 | Deleterious | −3.489 | E | 7 |
7 | rs1209841496 | R602W | Deleterious | 0 | PB | 1 | Deleterious | −7.4 | E | 87 |
8 | rs137852678 | L600P | Deleterious | 0 | PB | 1 | Deleterious | −6.472 | E | 91 |
9 | rs1398307167 | P596Q | Deleterious | 0.01 | PS | 0.951 | Deleterious | −3.974 | E | 57 |
10 | rs1398307167 | P596L | Deleterious | 0.01 | PB | 0.966 | Deleterious | −5.853 | E | 58 |
11 | rs767475430 | I578N | Deleterious | 0 | PB | 0.1 | Deleterious | −6.392 | E | 80 |
12 | rs113988352 | I561T | Deleterious | 0 | PB | 0.981 | Deleterious | −4.143 | E | 38 |
13 | rs1803838 | P538L | Deleterious | 0.03 | PS | 0.588 | Deleterious | −3.097 | E | 10 |
14 | rs916580554 | W504C | Deleterious | 0 | PB | 0.999 | Deleterious | −11,937 | E | 48 |
15 | rs1185249247 | S503N | Deleterious | 0 | PS | 0.949 | Deleterious | −2.621 | E | 50 |
16 | rs866554932 | P481R | Deleterious | 0.04 | PB | 0.1 | Deleterious | −6.451 | E | 6 |
17 | rs935654762 | V455A | Deleterious | 0 | PB | 0.997 | Deleterious | −3.255 | E | 35 |
18 | rs527393923 | T450M | Deleterious | 0 | PB | 0.1 | Deleterious | −4.499 | E | 56 |
19 | rs760409880 | L448F | Deleterious | 0 | PS | 0.816 | Deleterious | −3.135 | E | 33 |
20 | rs776192196 | P326L | Deleterious | 0 | PB | 0.996 | Deleterious | −6.947 | E | 29 |
21 | rs763976174 | R304H | Deleterious | 0.02 | PS | 0.850 | Deleterious | −2.809 | E | 52 |
22 | rs751403509 | R304C | Deleterious | 0 | PB | 0.1 | Deleterious | −4.767 | E | 39 |
23 | rs779371351 | I248N | Deleterious | 0 | PB | 0.1 | Deleterious | −5.877 | E | 38 |
24 | rs779371351 | I248T | Deleterious | 0 | PB | 0.1 | Deleterious | −4.218 | E | 42 |
25 | rs1017740241 | C247Y | Deleterious | 0 | PB | 0.1 | Deleterious | −9.192 | E | 49 |
26 | rs763588438 | V149G | Deleterious | 0.01 | PB | 0.987 | Deleterious | −4.942 | E | 48 |
27 | rs1482374494 | A119T | Deleterious | 0 | PB | 0.1 | Deleterious | −3.113 | E | 26 |
28 | rs756147217 | P98S | Deleterious | 0 | PB | 0.1 | Deleterious | −5.885 | E | 48 |
29 | rs865962653 | S51L | Deleterious | 0.01 | PS | 0.883 | Deleterious | −4.4 | E | 34 |
30 | rs781389511 | A46T | Deleterious | 0 | PB | 0.1 | Deleterious | −2.543 | E | 1 |
31 | rs34255470 | I30T | Deleterious | 0.02 | PB | 0.1 | Deleterious | −3.391 | E | 17 |
32 | rs11549696 | P27T | Deleterious | 0 | PB | 0.1 | Deleterious | −6.503 | E | 64 |
33 | rs1233778383 | W4C | Deleterious | 0 | PB | 1 | Deleterious | −10.563 | E | 23 |
PB: probably damaging, PS: possibly damaging, E: effect.
List of pathological nsSNPs predicted by PhD-SNP, SNPs and GO, P Mut, and PANTHER.
nsSNP | Amino Acid | P MUT | PhD-SNP | SNP&GO | PANTHER | |||||
---|---|---|---|---|---|---|---|---|---|---|
Prediction | Score | Prediction | RI | Prediction | RI | Effect | Preservation | |||
1 | rs1374373369 | D674V | Disease | 90% | Disease | 9 | Disease | 8 | probably damaging | 1036 |
2 | rs759271255 | I648T | Disease | 92% | Disease | 6 | Disease | 6 | probably damaging | 842 |
3 | rs752542806 | V642D | Disease | 87% | Disease | 5 | Disease | 8 | probably damaging | 455 |
4 | rs1209841496 | R602W | Disease | 93% | Disease | 8 | Disease | 3 | probably damaging | 1237 |
5 | rs137852678 | L600P | Disease | 93% | Disease | 9 | Disease | 8 | probably damaging | 1036 |
6 | rs767475430 | I578N | Disease | 93% | Disease | 8 | Disease | 6 | probably damaging | 1237 |
7 | rs916580554 | W504C | Disease | 91% | Disease | 6 | Disease | 8 | probably damaging | 750 |
8 | rs527393923 | T450M | Disease | 86% | Disease | 1 | Disease | 2 | probably damaging | 1036 |
9 | rs865962653 | S51L | Disease | 88% | Disease | 5 | Disease | 1 | probably damaging | 750 |
MutPred probability values of deleterious and pathogenic nsSNPs identified in STAT1.
SNP ID | Amino Acid | MutPred 2 Score | Affected PROSITE and ELM Motifs | Molecular Mechanisms | Probability | p-Value | |
---|---|---|---|---|---|---|---|
1 | rs1374373369 | D674V | 0.813 | Gain of Strand | 0.26 | 0.04 | |
2 | rs759271255 | I648T | 0.893 | Altered Stability | 0.16 | 0.02 | |
3 | rs752542806 | V642D | 0.867 | ELME000063 ELME000085 | Altered Ordered interface | 0.35 | 4.2 × 10−3 |
Gain of Relative solvent | 0.30 | 7.3 × 10−3 | |||||
Altered Transmembrane protein | 0.18 | 8.6 × 10−3 | |||||
Altered DNA binding | 0.15 | 0.04 | |||||
4 | rs1209841496 | R602W | 0.896 | - | Gain of Strand | 0.27 | 0.02 |
Altered Stability | 0.09 | 0.05 | |||||
5 | rs137852678 | L600P | 0.965 | ELME000052 | Gain of Intrinsic disorder | 0.31 | 0.04 |
Altered Stability | 0.28 | 6.6 × 10−3 | |||||
6 | rs767475430 | I578N | 0.936 | PS00008 | |||
7 | rs916580554 | W504C | 0.807 | ELME000197 | |||
8 | rs527393923 | T450M | 0.373 | ||||
9 | rs865962653 | S51L | 0.665 | ELME000063 | Altered transmembrane protein | 0.23 | 2.4 × 10−3 |
p-values ≤ 0.05.
Deleterious and pathogenic nsSNPs predicted to have a significant decrease on protein stability by I-MUTANT 3.0 algorithm, MUpro, and DDMUT.
SNP ID | Amino Acid | I Mutant 3 | MUpro | DDMUT | |||||
---|---|---|---|---|---|---|---|---|---|
Stability | RI | DDG | Stability | DDG (kcal/mol) | Stability | DDG (kcal/mol) | |||
1 | rs759271255 | I648T | Decrease | 9 | −2.43 | Decrease | −2.4802937 | Destabilizing | −2.93 |
2 | rs752542806 | V642D | Decrease | 8 | −1.85 | Decrease | −1.8071037 | Destabilizing | −1.11 |
3 | rs1209841496 | R602W | Decrease | 3 | −0.20 | Decrease | −1.0486884 | Destabilizing | −0.19 |
4 | rs137852678 | L600P | Decrease | 3 | −1.54 | Decrease | −1.6074419 | Destabilizing | −3.06 |
5 | rs767475430 | I578N | Decrease | 5 | −1.92 | Decrease | −0.98144877 | Destabilizing | −0.84 |
6 | rs916580554 | W504C | Decrease | 8 | −1.41 | Decrease | −0.86533645 | Destabilizing | −0.73 |
Alpha-Missense prediction of the pathogenic nsSNPs in STAT1.
SNP ID | Substitution | Alpha-Missense Pathogenicity | Alpha-Missense Prediction | |
---|---|---|---|---|
1 | rs759271255 | I648T | 0.9875 | Likely Pathogenic |
2 | rs752542806 | V642D | 0.9916 | Likely Pathogenic |
3 | rs1209841496 | R602W | 0.9982 | Likely Pathogenic |
4 | rs137852678 | L600P | 0.9998 | Likely Pathogenic |
5 | rs767475430 | I578N | 0.9986 | Likely Pathogenic |
6 | rs916580554 | W504C | 0.9815 | Likely Pathogenic |
Conservation profile of most damaging nsSNPs of STAT1.
SNP ID | Amino Acid | Conservation | Prediction | |
---|---|---|---|---|
1 | rs759271255 | I648T | 8 | Conserved and buried |
2 | rs752542806 | V642D | 6 | Buried |
3 | rs120984149 | R602W | 9 | (functional residues), highly conserved and exposed |
4 | rs137852678 | L600P | 9 | (structural residues), highly conserved and buried |
5 | rs767475430 | I578N | 9 | (structural residues), highly conserved and buried |
Changes in physical properties between wild-type and mutant residues predicted by project hope.
SNPs | Difference in Size | Difference in Charge | Difference in Hydrophobicity | Disrupt Hydrogen Bond | Affect Contact with Ligand Molecules | |
---|---|---|---|---|---|---|
1 | I648T | Yes | No | Yes | No | Yes |
2 | V642D | Yes | No | Yes | No | Yes |
3 | R602W | Yes | No | Yes | No | Yes |
4 | L600P | Yes | No | Yes | No | No |
5 | I578N | Yes | No | Yes | Yes | Yes |
6 | W504C | Yes | No | Yes | No | Yes |
Domain regions of the selected most damaging nsSNPs in STAT1.
STAT1 Domains (Position) | SNPs |
---|---|
STAT1, SH2 domain (557–707) | Y668F, I648T, V642D, R602W, and L600P |
STAT1 transcription factor, DNA binding domain (323–458) | R304C |
SH2 domain (578–638) | I578N |
Src homology 2 (SH2) domain profile (573–670) | I578N |
Supplementary Materials
The following supporting information can be downloaded at:
References
1. Awasthi, N.; Liongue, C.; Ward, A.C. STAT proteins: A kaleidoscope of canonical and non-canonical functions in immunity and cancer. J. Hematol. Oncol.; 2021; 14, 198. [DOI: https://dx.doi.org/10.1186/s13045-021-01214-y] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34809691]
2. Calò, V.; Migliavacca, M.; Bazan, V.; Macaluso, M.; Buscemi, M.; Gebbia, N.; Russo, A. STAT proteins: From normal control of cellular events to tumorigenesis. J. Cell. Physiol.; 2003; 197, pp. 157-168. [DOI: https://dx.doi.org/10.1002/jcp.10364] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/14502555]
3. Zhong, M.; Henriksen, M.A.; Takeuchi, K.; Schaefer, O.; Liu, B.; ten Hoeve, J.; Ren, Z.; Mao, X.; Chen, X.; Shuai, K. et al. Implications of an antiparallel dimeric structure of nonphosphorylated STAT1 for the activation-inactivation cycle. Proc. Natl. Acad. Sci. USA; 2005; 102, pp. 3966-3971. [DOI: https://dx.doi.org/10.1073/pnas.0501063102]
4. Mao, X.; Ren, Z.; Parker, G.N.; Sondermann, H.; Pastorello, M.A.; Wang, W.; McMurray, J.S.; Demeler, B.; Darnell, J.E.; Chen, X. Structural bases of unphosphorylated STAT1 association and receptor binding. Mol. Cell; 2005; 17, pp. 761-771. [DOI: https://dx.doi.org/10.1016/j.molcel.2005.02.021] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15780933]
5. Metwally, H.; Kishimoto, T. Distinct Phosphorylation of STAT1 Confers Distinct DNA Binding and Gene-regulatory Properties. J. Cell. Signal.; 2020; 1, pp. 50-55. [DOI: https://dx.doi.org/10.33696/Signaling.1.012]
6. Lorenzini, T.; Dotta, L.; Giacomelli, M.; Vairo, D.; Badolato, R. STAT mutations as program switchers: Turning primary immunodeficiencies into autoimmune diseases. J. Leukoc. Biol.; 2017; 101, pp. 29-38. [DOI: https://dx.doi.org/10.1189/jlb.5RI0516-237RR]
7. Asano, T.; Utsumi, T.; Kagawa, R.; Karakawa, S.; Okada, S. Inborn errors of immunity with loss- and gain-of-function germline mutations in STAT1. Clin. Exp. Immunol.; 2023; 212, pp. 96-106. [DOI: https://dx.doi.org/10.1093/cei/uxac106]
8. Mizoguchi, Y.; Okada, S. Inborn errors of STAT1 immunity. Curr. Opin. Immunol.; 2021; 72, pp. 59-64. [DOI: https://dx.doi.org/10.1016/j.coi.2021.02.009]
9. Verhoeven, Y.; Tilborghs, S.; Jacobs, J.; De Waele, J.; Quatannens, D.; Deben, C.; Prenen, H.; Pauwels, P.; Trinh, X.B.; Wouters, A. et al. The potential and controversy of targeting STAT family members in cancer. Semin. Cancer Biol.; 2020; 60, pp. 41-56. [DOI: https://dx.doi.org/10.1016/j.semcancer.2019.10.002]
10. Liongue, C.; Sobah, M.L.; Ward, A.C. Signal transducer and activator of transcription proteins at the nexus of immunodeficiency, autoimmunity and cancer. Biomedicines; 2023; 12, 45. [DOI: https://dx.doi.org/10.3390/biomedicines12010045]
11. Reich, N.C. STATs get their move on. Jak-stat; 2013; 2, e27080. [DOI: https://dx.doi.org/10.4161/jkst.27080] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24470978]
12. de Prati, A.C.; Ciampa, A.R.; Cavalieri, E.; Zaffini, R.; Darra, E.; Menegazzi, M.; Suzuki, H.; Mariotto, S. STAT1 as a new molecular target of anti-inflammatory treatment. Curr. Med. Chem.; 2005; 12, pp. 1819-1828. [DOI: https://dx.doi.org/10.2174/0929867054546645] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16101503]
13. Tolomeo, M.; Cavalli, A.; Cascio, A. STAT1 and its crucial role in the control of viral infections. Int. J. Mol. Sci.; 2022; 23, 4095. [DOI: https://dx.doi.org/10.3390/ijms23084095] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35456913]
14. Asano, T.; Noma, K.; Mizoguchi, Y.; Karakawa, S.; Okada, S. Human STAT1 gain of function with chronic mucocutaneous candidiasis: A comprehensive review for strengthening the connection between bedside observations and laboratory research. Immunol. Rev.; 2024; 322, pp. 81-97. [DOI: https://dx.doi.org/10.1111/imr.13300]
15. Shuai, K.; Schindler, C.; Prezioso, V.R.; Darnell, J.E. Activation of transcription by IFN-γ: Tyrosine phosphorylation of a 91-kD DNA binding protein. Science; 1992; 258, pp. 1808-1812. [DOI: https://dx.doi.org/10.1126/science.1281555]
16. Heim, M.H. The Jak-STAT pathway: Cytokine signalling from the receptor to the nucleus. J. Recept. Signal Transduct.; 1999; 19, pp. 75-120. [DOI: https://dx.doi.org/10.3109/10799899909036638]
17. Eilers, A.; Georgellis, D.; Klose, B.; Schindler, C.; Ziemiecki, A.; Harpur, A.G.; Wilks, A.F.; Decker, T. Differentiation-regulated serine phosphorylation of STAT1 promotes GAF activation in macrophages. Mol. Cell. Biol.; 1995; 15, pp. 3579-3586. [DOI: https://dx.doi.org/10.1128/MCB.15.7.3579]
18. Decker, T.; Lew, D.J.; Mirkovitch, J.; Darnell, J.E. Cytoplasmic activation of GAF, an IFN-gamma-regulated DNA-binding factor. EMBO J.; 1991; 10, pp. 927-932. [DOI: https://dx.doi.org/10.1002/j.1460-2075.1991.tb08026.x]
19. Meesilpavikkai, K.; Hirankarn, N.; Dalm, V.A.S.H.; van Hagen, P.M.; Dik, W.A.; IJspeert, H. Unraveling the Immunogenetics of STAT Proteins: Clinical Perspectives on Gain-of-Function and Loss-of-Function Variants. Asian Pac. J. Allergy Immunol.; 2024; 42, pp. 105-122. [DOI: https://dx.doi.org/10.12932/AP-270124-1776]
20. Chen, X.; Chen, J.; Chen, R.; Mou, H.; Sun, G.; Yang, L.; Jia, Y.; Zhao, Q.; Wen, W.; Zhou, L. et al. Genetic and Functional Identifying of Novel STAT1 Loss-of-Function Mutations in Patients with Diverse Clinical Phenotypes. J. Clin. Immunol.; 2022; 42, pp. 1778-1794. [DOI: https://dx.doi.org/10.1007/s10875-022-01339-w]
21. Boisson-Dupuis, S.; Kong, X.-F.; Okada, S.; Cypowyj, S.; Puel, A.; Abel, L.; Casanova, J.-L. Inborn errors of human STAT1: Allelic heterogeneity governs the diversity of immunological and infectious phenotypes. Curr. Opin. Immunol.; 2012; 24, pp. 364-378. [DOI: https://dx.doi.org/10.1016/j.coi.2012.04.011] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22651901]
22. Tsumura, M.; Okada, S.; Sakai, H.; Yasunaga, S.; Ohtsubo, M.; Murata, T.; Obata, H.; Yasumi, T.; Kong, X.-F.; Abhyankar, A. et al. Dominant-negative STAT1 SH2 domain mutations in unrelated patients with Mendelian susceptibility to mycobacterial disease. Hum. Mutat.; 2012; 33, pp. 1377-1387. [DOI: https://dx.doi.org/10.1002/humu.22113] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22573496]
23. Uzel, G.; Sampaio, E.P.; Lawrence, M.G.; Hsu, A.P.; Hackett, M.; Dorsey, M.J.; Noel, R.J.; Verbsky, J.W.; Freeman, A.F.; Janssen, E. et al. Dominant gain-of-function STAT1 mutations in FOXP3 wild-type immune dysregulation-polyendocrinopathy-enteropathy-X-linked-like syndrome. J. Allergy Clin. Immunol.; 2013; 131, pp. 1611-1623. [DOI: https://dx.doi.org/10.1016/j.jaci.2012.11.054] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23534974]
24. Hartono, S.P.; Vargas-Hernández, A.; Ponsford, M.J.; Chinn, I.K.; Jolles, S.; Wilson, K.; Forbes, L.R. Novel STAT1 Gain-of-Function Mutation Presenting as Combined Immunodeficiency. J. Clin. Immunol.; 2018; 38, pp. 753-756. [DOI: https://dx.doi.org/10.1007/s10875-018-0554-3]
25. Henrickson, S.E.; Dolan, J.G.; Forbes, L.R.; Vargas-Hernández, A.; Nishimura, S.; Okada, S.; Kersun, L.S.; Brodeur, G.M.; Heimall, J.R. Gain-of-Function STAT1 Mutation With Familial Lymphadenopathy and Hodgkin Lymphoma. Front. Pediatr.; 2019; 7, 160. [DOI: https://dx.doi.org/10.3389/fped.2019.00160]
26. Okada, S.; Asano, T.; Moriya, K.; Boisson-Dupuis, S.; Kobayashi, M.; Casanova, J.-L.; Puel, A. Human STAT1 Gain-of-Function Heterozygous Mutations: Chronic Mucocutaneous Candidiasis and Type I Interferonopathy. J. Clin. Immunol.; 2020; 40, pp. 1065-1081. [DOI: https://dx.doi.org/10.1007/s10875-020-00847-x]
27. Wang, X.; Zhao, W.; Chen, F.; Zhou, P.; Yan, Z. Chinese Pedigree of Chronic Mucocutaneous Candidiasis Due to STAT1 Gain-of-Function Mutation: A Case Study and Literature Review. Mycopathologia; 2023; 188, pp. 87-97. [DOI: https://dx.doi.org/10.1007/s11046-022-00685-y]
28. Egri, N.; Esteve-Solé, A.; Deyà-Martínez, À.; Ortiz de Landazuri, I.; Vlagea, A.; García, A.P.; Cardozo, C.; Garcia-Vidal, C.; Bartolomé, C.S.; Español-Rego, M. et al. Primary immunodeficiency and chronic mucocutaneous candidiasis: Pathophysiological, diagnostic, and therapeutic approaches. Allergol. Immunopathol.; 2021; 49, pp. 118-127. [DOI: https://dx.doi.org/10.15586/aei.v49i1.20]
29. Collins, F.S.; Brooks, L.D.; Chakravarti, A. A DNA polymorphism discovery resource for research on human genetic variation. Genome Res.; 1998; 8, pp. 1229-1231. [DOI: https://dx.doi.org/10.1101/gr.8.12.1229]
30. Arshad, S.; Ishaque, I.; Mumtaz, S.; Rashid, M.U.; Malkani, N. In-Silico Analyses of Nonsynonymous Variants in the BRCA1 Gene. Biochem. Genet.; 2021; 59, pp. 1506-1526. [DOI: https://dx.doi.org/10.1007/s10528-021-10074-7]
31. Yazar, M.; Özbek, P. In Silico Tools and Approaches for the Prediction of Functional and Structural Effects of Single-Nucleotide Polymorphisms on Proteins: An Expert Review. OMICS J. Integr. Biol.; 2021; 25, pp. 23-37. [DOI: https://dx.doi.org/10.1089/omi.2020.0141] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33058752]
32. Allemailem, K.S.; Almatroudi, A.; Alrumaihi, F.; Makki Almansour, N.; Aldakheel, F.M.; Rather, R.A.; Afroze, D.; Rah, B. Single nucleotide polymorphisms (SNPs) in prostate cancer: Its implications in diagnostics and therapeutics. Am. J. Transl. Res.; 2021; 13, pp. 3868-3889. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34017579]
33. Clifford, R.J.; Edmonson, M.N.; Nguyen, C.; Scherpbier, T.; Hu, Y.; Buetow, K.H. Bioinformatics tools for single nucleotide polymorphism discovery and analysis. Ann. N. Y. Acad. Sci.; 2004; 1020, pp. 101-109. [DOI: https://dx.doi.org/10.1196/annals.1310.011]
34. Artimo, P.; Jonnalagedda, M.; Arnold, K.; Baratin, D.; Csardi, G.; de Castro, E.; Duvaud, S.; Flegel, V.; Fortier, A.; Gasteiger, E. et al. ExPASy: SIB bioinformatics resource portal. Nucleic Acids Res.; 2012; 40, pp. W597-W603. [DOI: https://dx.doi.org/10.1093/nar/gks400]
35. 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, pp. 1073-1081. [DOI: https://dx.doi.org/10.1038/nprot.2009.86]
36. Adzhubei, I.; Jordan, D.M.; Sunyaev, S.R. Predicting functional effect of human missense mutations using PolyPhen-2. Curr. Protoc. Hum. Genet.; 2013; 76, pp. 7.20.1-7.20.41. [DOI: https://dx.doi.org/10.1002/0471142905.hg0720s76]
37. Choi, Y.; Chan, A.P. PROVEAN web server: A tool to predict the functional effect of amino acid substitutions and indels. Bioinformatics; 2015; 31, pp. 2745-2747. [DOI: https://dx.doi.org/10.1093/bioinformatics/btv195]
38. Bromberg, Y.; Yachdav, G.; Rost, B. SNAP predicts effect of mutations on protein function. Bioinformatics; 2008; 24, pp. 2397-2398. [DOI: https://dx.doi.org/10.1093/bioinformatics/btn435]
39. Capriotti, E.; Calabrese, R.; Casadio, R. Predicting the insurgence of human genetic diseases associated to single point protein mutations with support vector machines and evolutionary information. Bioinformatics; 2006; 22, pp. 2729-2734. [DOI: https://dx.doi.org/10.1093/bioinformatics/btl423]
40. 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, (Suppl. S3), S6. [DOI: https://dx.doi.org/10.1186/1471-2164-14-S3-S6]
41. 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, pp. W222-W228. [DOI: https://dx.doi.org/10.1093/nar/gkx313] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28453649]
42. Mi, H.; Muruganujan, A.; Thomas, P.D. PANTHER in 2013: Modeling the evolution of gene function, and other gene attributes, in the context of phylogenetic trees. Nucleic Acids Res.; 2013; 41, pp. D377-D386. [DOI: https://dx.doi.org/10.1093/nar/gks1118] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23193289]
43. Pejaver, V.; Urresti, J.; Lugo-Martinez, J.; Pagel, K.A.; Lin, G.N.; Nam, H.-J.; Mort, M.; Cooper, D.N.; Sebat, J.; Iakoucheva, L.M. et al. Inferring the molecular and phenotypic impact of amino acid variants with MutPred2. Nat. Commun.; 2020; 11, 5918. [DOI: https://dx.doi.org/10.1038/s41467-020-19669-x]
44. Capriotti, E.; Fariselli, P.; Calabrese, R.; Casadio, R. Predicting protein stability changes from sequences using support vector machines. Bioinformatics; 2005; 21, (Suppl. S2), pp. ii54-ii58. [DOI: https://dx.doi.org/10.1093/bioinformatics/bti1109]
45. Cheng, J.; Randall, A.; Baldi, P. Prediction of protein stability changes for single-site mutations using support vector machines. Proteins Struct. Funct. Bioinform.; 2006; 62, pp. 1125-1132. [DOI: https://dx.doi.org/10.1002/prot.20810]
46. Zhou, Y.; Pan, Q.; Pires, D.E.V.; Rodrigues, C.H.M.; Ascher, D.B. DDMut: Predicting effects of mutations on protein stability using deep learning. Nucleic Acids Res.; 2023; 51, pp. W122-W128. [DOI: https://dx.doi.org/10.1093/nar/gkad472]
47. Cheng, J.; Novati, G.; Pan, J.; Bycroft, C.; Žemgulytė, A.; Applebaum, T.; Pritzel, A.; Wong, L.H.; Zielinski, M.; Sargeant, T. et al. Accurate proteome-wide missense variant effect prediction with AlphaMissense. Science; 2023; 381, eadg7492. [DOI: https://dx.doi.org/10.1126/science.adg7492]
48. Varadi, M.; Anyango, S.; Deshpande, M.; Nair, S.; Natassia, C.; Yordanova, G.; Yuan, D.; Stroe, O.; Wood, G.; Laydon, A. et al. AlphaFold Protein Structure Database: Massively expanding the structural coverage of protein-sequence space with high-accuracy models. Nucleic Acids Res.; 2022; 50, pp. D439-D444. [DOI: https://dx.doi.org/10.1093/nar/gkab1061]
49. Meng, E.C.; Goddard, T.D.; Pettersen, E.F.; Couch, G.S.; Pearson, Z.J.; Morris, J.H.; Ferrin, T.E. UCSF ChimeraX: Tools for structure building and analysis. Protein Sci.; 2023; 32, e4792. [DOI: https://dx.doi.org/10.1002/pro.4792]
50. Venselaar, H.; Te Beek, T.A.H.; Kuipers, R.K.P.; Hekkelman, M.L.; Vriend, G. Protein structure analysis of mutations causing inheritable diseases. An e-Science approach with life scientist friendly interfaces. BMC Bioinform.; 2010; 11, 548. [DOI: https://dx.doi.org/10.1186/1471-2105-11-548]
51. Ashkenazy, H.; Erez, E.; Martz, E.; Pupko, T.; Ben-Tal, N. ConSurf 2010: Calculating evolutionary conservation in sequence and structure of proteins and nucleic acids. Nucleic Acids Res.; 2010; 38, pp. W529-W533. [DOI: https://dx.doi.org/10.1093/nar/gkq399] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20478830]
52. Mulder, N.J.; Apweiler, R.; Attwood, T.K.; Bairoch, A.; Bateman, A.; Binns, D.; Biswas, M.; Bradley, P.; Bork, P.; Bucher, P. et al. InterPro Consortium InterPro: An integrated documentation resource for protein families, domains and functional sites. Brief Bioinform.; 2002; 3, pp. 225-235. [DOI: https://dx.doi.org/10.1093/bib/3.3.225] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/12230031]
53. Szklarczyk, D.; Franceschini, A.; Kuhn, M.; Simonovic, M.; Roth, A.; Minguez, P.; Doerks, T.; Stark, M.; Muller, J.; Bork, P. et al. The STRING database in 2011: Functional interaction networks of proteins, globally integrated and scored. Nucleic Acids Res.; 2011; 39, pp. D561-D568. [DOI: https://dx.doi.org/10.1093/nar/gkq973] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21045058]
54. Zhang, W.; Chen, X.; Gao, G.; Xing, S.; Zhou, L.; Tang, X.; Zhao, X.; An, Y. Clinical Relevance of Gain- and Loss-of-Function Germline Mutations in STAT1: A Systematic Review. Front. Immunol.; 2021; 12, 654406. [DOI: https://dx.doi.org/10.3389/fimmu.2021.654406]
55. Dupuis, S.; Jouanguy, E.; Al-Hajjar, S.; Fieschi, C.; Al-Mohsen, I.Z.; Al-Jumaah, S.; Yang, K.; Chapgier, A.; Eidenschenk, C.; Eid, P. et al. Impaired response to interferon-α/β and lethal viral disease in human STAT1 deficiency. Nat. Genet.; 2003; 33, pp. 388-391. [DOI: https://dx.doi.org/10.1038/ng1097]
56. Chapgier, A.; Boisson-Dupuis, S.; Jouanguy, E.; Vogt, G.; Feinberg, J.; Prochnicka-Chalufour, A.; Casrouge, A.; Yang, K.; Soudais, C.; Fieschi, C. et al. Novel STAT1 alleles in otherwise healthy patients with mycobacterial disease. PLoS Genet.; 2006; 2, e131. [DOI: https://dx.doi.org/10.1371/journal.pgen.0020131]
57. Liu, L.; Okada, S.; Kong, X.-F.; Kreins, A.Y.; Cypowyj, S.; Abhyankar, A.; Toubiana, J.; Itan, Y.; Audry, M.; Nitschke, P. et al. Gain-of-function human STAT1 mutations impair IL-17 immunity and underlie chronic mucocutaneous candidiasis. J. Exp. Med.; 2011; 208, pp. 1635-1648. [DOI: https://dx.doi.org/10.1084/jem.20110958]
58. Wang, X.; Zhang, R.; Wu, W.; Wang, A.; Wan, Z.; van de Veerdonk, F.L.; Li, R. New and recurrent STAT1 mutations in seven Chinese patients with chronic mucocutaneous candidiasis. Int. J. Dermatol.; 2017; 56, pp. e30-e33. [DOI: https://dx.doi.org/10.1111/ijd.13427]
59. Breuer, O.; Daum, H.; Cohen-Cymberknoh, M.; Unger, S.; Shoseyov, D.; Stepensky, P.; Keller, B.; Warnatz, K.; Kerem, E. Autosomal dominant gain of function STAT1 mutation and severe bronchiectasis. Respir. Med.; 2017; 126, pp. 39-45. [DOI: https://dx.doi.org/10.1016/j.rmed.2017.03.018]
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2025 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
Background: The Signal Transducer and Activator of Transcription 1 (STAT1) gene is an essential component of the JAK-STAT signaling pathway. This pathway plays a pivotal role in the regulation of different cellular processes, including immune responses, cell growth, and apoptosis. Mutations in the STAT1 gene contribute to a variety of immune system dysfunctions. Objectives: We aim to identify disease-susceptible single-nucleotide polymorphisms (SNPs) in STAT1 gene and predict structural changes associated with the mutations that disrupt normal protein–protein interactions using different computational algorithms. Methods: Several in silico tools, such as SIFT, Polyphen v2, PROVEAN, SNAP2, PhD-SNP, SNPs&GO, Pmut, and PANTHER, were used to determine the deleterious nsSNPs of the STAT1. Further, we evaluated the potentially deleterious SNPs for their effect on protein stability using I-Mutant, MUpro, and DDMUT. Additionally, we predicted the functional and structural effects of the nsSNPs using MutPred. We used Alpha-Missense to predict missense variant pathogenicity. Moreover, we predicted the 3D structure of STAT1 using an artificial intelligence system, alphafold, and the visualization of the 3D structures of the wild-type amino acids and the mutant residues was performed using ChimeraX 1.9 software. Furthermore, we analyzed the structural and conformational variations that have resulted from SNPs using Project Hope, while changes in the biological interactions between wild type, mutant amino acids, and neighborhood residues was studied using DDMUT. Conservational analysis and surface accessibility prediction of STAT1 was performed using ConSurf. We predicted the protein–protein interaction using STRING database. Results: In the current study, we identified six deleterious nsSNPs (R602W, I648T, V642D, L600P, I578N, and W504C) and their effect on protein structure, function, and stability. Conclusions: These findings highlight the potential of approaches to pinpoint pathogenic SNPs, providing a time- and cost-effective alternative to experimental approaches. To the best of our knowledge, this is the first comprehensive study in which we analyze STAT1 gene variants using both bioinformatics and artificial-intelligence-based model tools.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details
1 Department of Basic Medical Sciences, College of Medicine, Prince Sattam bin Abdulaziz University, Al Kharj 16278, Saudi Arabia
2 Department of Physiology, Faculty of Medicine, King Abdul-Aziz University, Rabigh 25724, Saudi Arabia
3 Plastic Surgery, Department of Surgery, College of Medicine, Prince Sattam bin Abdulaziz University, Al Kharj 16278, Saudi Arabia