Screening and insilico analysis of deleterious nsSNPs (missense) in human CSF3 for their effects on protein structure,stability and function

Human granulocyte colony stimulating factor (hG-CSF), known as CSF3, plays an important role in the growth, differentiation, proliferation, survival, and activation of the granulocyte cell lineage such as neutrophils and their precursors. Functional reduction in native CSF3 protein results in reduced proliferation and activation of neutrophils and leads to neutropenia. Single nucleotide polymorphisms (SNPs) in the CSF3 gene may have deleterious effects on the CSF3 protein function. This study was undertaken to find the functional SNPs in the human CSF3 gene. Results suggest that 18.9% of all the SNPs in the dbSNP database for CSF3 gene were present in the coding region. Out of 59 non-synonymous SNPs (nsSNPs), 26 nsSNPs were predicted to be non-tolerable by SIFT whereas 18 and 7 nsSNPs were predicted as probably damaging and possibly damaging, respectively by PolyPhen. Among this 31 nsSNPs, 16 nsSNPs were identified to be potentially deleterious by PANTHER server, and 4 nsSNPs were found to be neutral by PROVEAN. SNPAnalyzer predicted 7 nsSNPs to be neutral phenotype and the remaining 24 nsSNPs to be associated with diseases. Among the predicted nsSNPs, rs144408123, rs144408123, rs145136406, rs145311241, rs373191696, rs762945096, rs763688260, rs767572172, rs775326370, rs777777864, rs777983866, rs781596455, rs139072004, rs757612684, rs772911210, rs139072004, rs746634544, rs749993200, rs763426127, rs772466210 were identified as deleterious and potentially damaging. I-Mutant analysis revealed that th 20 nsSNPs were important for protein stability of CSF3. Therefore, th 20 nsSNPs may be used for the wider population-based genetic studies and also should be taken into account while engineering the recombinant CSF3 protein for clinical use.     

Computational analysis for the determination of deleterious nsSNPs in human MTHFD1 gene

Single nucleotide polymorphisms (SNPs) are the most common genetic polymorphisms and play a major role in many inherited diseases. Methylenetetrahydrofolate dehydrogenase 1 (MTHFD1) is one of the enzymes involved in folate metabolism. In the present study, the functional and structural consequences of nsSNPs of human MTHFD1 gene was analyzed using various computational tools like SIFT, PolyPhen2, PANTHER, PROVEAN, SNAP2, nsSNPAnalyzer, PhD-SNP, SNPs&GO, I-Mutant, MuPro, ConSurf, InterPro, NCBI Conserved Domain Search tool, ModPred, SPARKS-X, RAMPAGE, FT Site and PyMol. Out of 327 nsSNPs form human MTHFD1 gene, total 45 SNPs were predicted as functionally most significant SNPs, among which 17 were highly conserved and functional, 17 were highly conserved and structural residues. Among 45 most significant SNPs, 15 were predicted to be involved in post translational modifications. The p.Gly165Arg may interfere in homodimer interface formation. The p.Asn439Lys and p.Asp445Asn may interfere in binding interactions of MTHFD1 protein with cesium cation and potassium. The two SNPs (p.Asp562Gly and p.Gly637Cys) might interfere in interactions of MTHFD1 with ligand.     

Computational algorithmic and molecular dynamics study of functional and structural impacts of non-synonymous single nucleotide polymorphisms in human DHFR gene

Human dihydrofolate reductase (DHFR) is a conserved enzyme that is central to folate metabolism and is widely targeted in pathogenic diseases as well as cancers. Although studies have reported the fact that genetic mutations in DHFR leads to a rare autosomal recessive inborn error of folate metabolism and drug resistance, there is a lack of an extensive study on how the deleterious non-synonymous SNPs (nsSNPs) disrupt its phenotypic effects. In this study, we aim at discovering the structural and functional consequences of nsSNPs in DHFR by employing a combined computational approach consisting of ten recently developed in silico tools for identification of damaging nsSNPs and molecular dynamics (MD) simulation for getting deeper insights into the magnitudes of damaging effects. Our study revealed the presence of 12 most deleterious nsSNPs affecting the native phenotypic effects, with three (R71T, G118D, Y122D) identified in the co-factor and ligand binding active sites. MD simulations also suggested that these three SNPs particularly Y122D, alter the overall structural flexibility and dynamics of the native DHFR protein which can provide more understandings into the crucial roles of these mutants in influencing the loss of DHFR function.     

Computational analyses prioritize and reveal the deleterious nsSNPs in human angiotensinogen gene

Angiotensinogen (AGT) is a key component of renin-angiotensin-aldosterone system (RAAS), which plays central role in blood pressure homeostasis. Association of AGT polymorphisms have been investigated in different ethnic populations in variety of cardiovascular and non-cardiovascular conditions. In this study, 354 non-synonymous SNPs (nsSNPs) of AGT were evaluated to predict damaging and structurally important variants. Majority of the deleterious nsSNPs occurred in the evolutionary conserved regions. Several of these nsSNPs were found to affect post-translational modifications like methylation, glycosylation, phosphorylation, ubiquitination etc. Structural evaluations predicted 19 variants as destabilizing and some of them were also predicted to destabilize the renin-AGT interaction. Therefore, the present computational investigation predicted pathogenic and functionally important variants of human AGT gene. The study has also shown that AGT deregulation is associated with survival outcome in patients with gastric and breast cancer, using microarray gene expression profile. Furthermore, the computationally screened nsSNPs can be analyzed in population based genotyping studies and may help futuristic drug development in the area of AGT pharmacogenomics.     

Identification and structural characterization of deleterious non-synonymous single nucleotide polymorphisms in the human SKP2 gene

In SCF (Skp, Cullin, F-box) ubiquitin-protein ligase complexes, S-phase kinase 2 (SKP2) is one of the major players of F-box family, that is responsible for the degradation of several important cell regulators and tumor suppressor proteins. Despite of having significant evidence for the role of SKP2 on tumorgenesis, there is a lack of available data regarding the effect of non-synonymous polymorphisms. In this communication, the structural and functional consequences of non-synonymous single nucleotide polymorphisms (nsSNPs) of SKP2 have been reported by employing various computational approaches and molecular dynamics simulation. Initially, several computational tools like SIFT, PolyPhen-2, PredictSNP, I-Mutant 2.0 and ConSurf have been implicated in this study to explore the damaging SNPs. In total of 172 nsSNPs, 5 nsSNPs were identified as deleterious and 3 of them were predicted to be decreased the stability of protein. Guided from ConSurf analysis, P101L (rs761253702) and Y346C (rs755010517) were categorized as the highly conserved and functional disrupting mutations. Therefore, these mutations were subjected to three dimensional model building and molecular dynamics simulation study for the detailed structural consequences upon the mutations. The study revealed that P101L and Y346C mutations increased the flexibility and changed the structural dynamics. As both these mutations are located in the most functional regions of SKP2 protein, these computational insights might be helpful to consider these nsSNPs for wet-lab confirmatory analysis as well as in rationalizing future population based studies and structure based drug design against SKP2.     

Seven nonsynonymous SNPs in the gene encoding human deoxyribonuclease II may serve as a functional SNP potentially implicated in autoimmune dysfunction

Many nonsynonymous SNPs in the human DNase II gene (DNASE2), potentially relevant to autoimmunity in conditions such as rheumatoid arthritis, have been identified, but only limited population data are available and no studies have evaluated whether such SNPs are functional. Genotyping of all the 15 nonsynonymous human DNase II SNPs was performed in three ethnic groups including 16 different populations using the PCR‐restriction fragment length polymorphism technique. A series of constructs corresponding to each SNP was examined. Fifteen nonsynonymous SNPs in the gene, except for p.Val206Ile in a Korean population, exhibited a mono‐allelic distribution in all of the populations. On the basis of alterations in the activity levels resulting from the corresponding amino acid substitutions, four activity‐abolishing and five activity‐reducing SNPs were confirmed to be functional. The amino acid residues in activity‐abolishing SNPs were conserved in animal DNase II. All the nonsynonymous SNPs that affected the catalytic activity of human DNase II showed extremely low genetic heterogeneity. However, a minor allele of seven SNPs producing a loss‐of‐function or extremely low activity‐harboring variant could serve as a genetic risk factor for autoimmune dysfunction. These functional SNPs in DNASE2 may have clinical implications in relation to the prevalence of autoimmune diseases.     

Analysis of the variations in IL-28RA gene and their association with allergic rhinitis

IL-28RA is one of the important candidate genes for complex trait of genetic diseases, but there is no published information of the genetic variation in this gene. We scanned the seven exons and their boundary introns sequence of IL-28RA including the promoter regions to analyze genetic variation sites, and identified eighteen single nucleotide polymorphisms (SNPs) and two variation sites. We chose seven SNPs (g.-1193 A>C, g.-30 C>T, g.17654 C>T, g.27798 A>G, g.31265 C>T, g.31911 C>T and g.32349 G>A) of them for large sample size genotyping, and assessed the association of genotype and allele frequencies of these SNPs between allergic rhinitis patients and non-allergic rhinitis controls. We also compared the genotype frequencies between Korean controls and Han Chinese control or Korean Chinese control. We investigated the frequencies of haplotype constructed by these SNPs between allergic rhinitis patients and non-allergic rhinitis controls. Our results suggested that the g.32349 G>A polymorphism of IL-28RA might be associated with susceptibility to allergic rhinitis (P=0.032), but seems to have no relationship with serum total IgE levels. The haplotype frequencies by these SNPs also show significant association between controls and allergic rhinitis patients.     

Assessment of structurally and functionally high-risk nsSNPs impacts on human bone morphogenetic protein receptor type IA (BMPR1A) by computational approach

BMPR1A (BMP type 1 receptor) is a transmembrane cell-surface receptor also known as ALK3 (activin-like kinases-3) encodes for a type I serine/threonine kinase receptor and a member of the transforming growth-factor β–receptor (TGF-β) super family. The BMPR1A has a significant interaction with BMP-2 for protein activity and also has a low affinity with growth and differentiation factor 5 (GDF5); positively regulates chondrocyte differentiation. The genetic variations can alter the structure and function of the BMPR1A gene that causes several diseases such as juvenile polyposis syndrome or hereditary cancer-predisposing syndrome. The current study was carried out to identify potential deleterious non-synonymous single nucleotide polymorphisms (nsSNPs) in BMPR1A by implementing different computational algorithms such as SIFT, PolyPhen2, SNAP2, PROVEAN, PhD-SNP, SNPs&GO, nsSNPAnalyzer, and P-Mut. From 205 nsSNPs in BMPR1A, 7 nsSNPs (C76Y, C124R, C124Y, C376Y, R443C, R480W, and W487R) were predicted as deleterious in 8 prediction algorithms. The Consurf analysis showed that selected 7 nsSNPs were present in the highly conserved regions. Molecular dynamics simulation analysis also performed to explore conformational changes in the variant structure with respect to its native structure. According to the MDS result, all variants flexibility and rigidity were unbalanced, which may alter the structural and functional behavior of the native protein. Although, three nsSNPs i.e., C124R, C376Y, and R443C have already been reported in patients associated with JPS, but their structural and functional molecular studies remain uncharacterized. Therefore, the findings of this study can provide a better understanding of uncharacterized nsSNPS and to find their association with disease susceptibility and also facilitate to the researchers for designing or developing the target dependent drugs.     

In-Silico analysis of missense SNPs in Human HPPD gene associated with Tyrosinemia type III and Hawkinsinuria

HPPD gene codes a dioxygenase enzyme involved in catalysis of different molecules such as tyrosine and phenylalanine by oxidizing them to produce energy. A single change in protein can trigger serious genetic disorders like Tyrosinemia type III and Hawkinsinuria. This study aims to identify the functional missense SNPs of the HPPD gene by using multiple computational tools. All deleterious missense SNPs retrieved from Ensembl and OMIM database were evaluated through six different software. Ultimately, out of 148 missense SNPs, only 27 were confirmed as diseasecausing SNPs by developing a consensus approach. These damaging SNPs were further examined to evaluate their impact on protein stability and energy including their evolutionary conservation. Native and mutated proteins structures were also designed and superimposed by I-TASSER and PyMol respectively. This work results in narrowing down missense SNPs which are still not confirmed experimentally and demands the confirmation by GWAS data. Thus, these missense SNPs could directly or indirectly destabilize the amino acid interactions causing functional deviations of protein.     

The Novel Role of MIF in the Secretion of IL-25, IL-31, and IL-33 from PBMC of Patients with Rheumatoid Arthritis

Rheumatoid arthritis (RA) is an autoimmune inflammatory joint disease with complex pathogenesis associated with cytokine dysregulation. Macrophage migration inhibitory factor (MIF) plays a role in systemic inflammation and joint destruction in RA and could be associated with the secretion of other immune-modulatory cytokines such as IL-25, IL-31, and IL-33. For the above, our main aim was to evaluate the IL-25, IL-31, and IL-33 secretion from recombinant human MIF (rhMIF)-stimulated peripheral blood mononuclear cells (PBMC) of RA patients. The rhMIF and lipopolysaccharide (LPS) plus rhMIF stimuli promote the secretion of IL-25, IL-31, and IL-33 (p < 0.05) from PBMC of RA patients. The study groups, the different stimuli, and the interaction between both showed a statistically significant effect on the secretion of IL-25 (p < 0.05) and IL-31 (p < 0.01). The study of the effect of the RA patient treatments and their interaction with the effect of stimuli did not show an interaction between them. In conclusion, our study generates new evidence for the role of MIF in the secretion of IL-25, IL-31, and IL-33 and its immunomodulatory effect on RA.     

Genetic and expression analysis of SNPs in the human deoxyribonuclease II: SNPs in the promoter region reduce its in vivo activity through decreased promoter activity

Five SNPs in the human DNase II gene have been reported to be associated with rheumatoid arthritis (RA). Genotype and haplotype analysis of 14 SNPs, nine SNPs of which reported in the NCBI dbSNP database in addition to these five SNPs, was performed in healthy subjects. The enzymatic activities of the amino acid substituted DNase II corresponding to each SNP and serum DNase II in healthy Japanese, and promoter activities derived from each haplotype of the RA‐related SNPs were measured. Significant correlations between genotype in each RA‐related SNP and enzymatic activity levels were found; alleles associated with RA exhibited a reduction in serum DNase II activity. Furthermore, the promoter activities of each reporter construct corresponding to predominant haplotypes in three SNPs in the promoter region of the gene exhibited significant correlation with levels of serum DNase II activity. These findings indicate these three SNPs could alter the promoter activity of DNASE2, leading to a decline in DNase II activity in the serum through gene expression. Since the three SNPs in the promoter region of the DNase II gene could affect in vivo DNase II activity through reduction of the promoter activity, it is feasible to identify these SNPs susceptible to RA.     

Rational Design,Synthesis and Evaluation of Oxazolo[4,5-c]-quinolinone Analogs as Novel Interleukin-33 Inhibitors

Interleukin-33 (IL-33) is an epithelial-derived cytokine that plays an important role in immune-mediated diseases such as asthma, atopic dermatitis, and rheumatoid arthritis. Although IL-33 is considered a potential target for the treatment of allergy-related diseases, no small molecule that inhibits IL-33 has been reported. Based on the structure-activity relationship and in vitro 2D NMR studies employing ¹⁵N-labeled IL-33, we identified that the oxazolo[4,5-c]-quinolinone analog 7 c binds to the interface region of IL-33 and IL-33 receptor (ST2), an orphan receptor of the IL-1 receptor family. Compound 7 c effectively inhibited the production of IL-6 in human mast cells in a dose-dependent manner. Compound 7 c is the first low molecular weight IL-33 inhibitor and may be used as a prototype molecule for structural optimization and investigation of the IL-33/ST2 signaling pathway.     

Confirmation that SNPs in the high mobility group‐A2 gene (HMGA2) are associated with adult height in the Japanese population; wide‐ranging population survey of height‐related SNPs in HMGA2

Adult height is a highly heritable trait in that multiple genes are involved. Recent genome‐wide association studies have identified a novel single‐nucleotide polymorphism (SNP) rs1042725 in the high mobility group‐A2 gene (HMGA2) and shown it to be associated with human height in Caucasian populations. We performed a replication study to examine the associations between SNPs in HMGA2 and adult height in the Japanese population based on autopsy cases. Although we could not confirm a significant association between rs1042725 in HMGA2 and adult height, another SNP, rs7968902, in the gene achieved significance for its association in the same populations, and the effect was the same as that documented previously. These findings permit us to conclude that the SNPs in HMGA2 are common variants influencing human height across different populations. Moreover, a worldwide population study of these SNPs using 14 different populations including Asians, Africans and Caucasians demonstrated that both haplotypes and genotypes for three height‐related SNPs (rs1042725, rs7968682 and rs7968902) in HMGA2 were distributed in an ethnicity‐dependent manner. This information will be useful for clarifying the genetic basis of human height.     

Computational prediction of the effects of non-synonymous single nucleotide polymorphisms on the GPI-anchor transamidase subunit GPI8p of Plasmodium falciparum

Drug resistance is increasingly evolving in malaria parasites; hence, it is important to discover and establish alternative drug targets. In this context, GPI-anchor transamidase (GPI-T) is a potential drug target primarily of its crucial role in the development and survival of the parasite in the GPI anchor biosynthesis pathway. The present investigation was undertaken to explore the plausible effects of nsSNP on the structure and functions of GPI-T subunit GPI8p of Plasmodium falciparum. The GPI8p (PF3D7_1128700) was analyzed using various sequence-based and structure-based computational tools such as SIFT, PROVEAN, PredictSNP, SNAP2, I-Mutant, MuPro, ConSurf, NetSurfP, MUSTER, COACH server and STRING server. Of the 34 nsSNPs submitted for functional analysis, 18 nsSNPs (R124 L, N143 K, Y145 F, V157I, T195S, K379E, I392 K, I437 T, Y438H, N439D, Y441H, N442D, N448D, N451D, D457A, D457Y, I458 L and N460 K) were predicted to have deleterious effects on the protein GPI8p. Additionally, I-Mutant 2.0 and MuPro both showed a decrease in stability after mutation as a result of these nsSNPs, suggesting the destabilization of protein. ConSurf findings suggest that most of the regions were highly conserved. In addition, COACH server was used to predict the ligand binding sites. It was found that no mutation was present at the predicted ligand binding site. The results of the STRING database showed that the protein GPI8p interacts with those proteins which either involve the biosynthetic process of attaching GPI anchor to protein or GPI anchor. The present study suggested that the GPI8p could be a novel target for anti-malarial drugs, which provides significant details for further experimentation.     

Global genetic analysis of all single nucleotide polymorphisms in exons of the human deoxyribonuclease I-like 3 gene and their effect on its catalytic activity

Deoxyribonucleases (DNases) have been suggested to be implicated in the pathophysiology of autoimmune diseases. In the DNASE1L3 gene encoding human DNase I‐like 3 (DNase 1L3), a member of the DNase I family, only two non‐synonymous (R178 H and R206C) single nucleotide polymorphisms (SNPs) have been examined [Ueki et al., Clin. Chim. Acta 2009, 407, 20–24]. Three other non‐synonymous (G82R, K96N, and I243M) and four synonymous (S17S, T84T, R92R, and A181A) SNPs, in addition to R206C and R178H, have been identified in DNASE1L3. We investigated the distribution of all these SNPs in exons of the gene in eight Asian, three African, and three Caucasian populations worldwide using newly devised genotyping methods. SNP T84T showed polymorphism in all the populations, and R92R was polymorphic in the three African and three Caucasian populations; R206C was distributed only in Caucasian populations. In contrast, no minor allele was found in five SNPs (S17S, G82R, K96N, A181A, and I243M) in DNASE1L3. Generally, the DNase 1L3 gene shows relatively low genetic diversity with regard to exonic SNPs. When the effect of amino acid/nucleotide substitutions resulting from the SNPs on DNase 1L3 activity was examined, none of the synonymous SNPs had any effect on the DNase 1L3 activity, whereas among non‐synonymous SNPs, SNP G82R diminished the activity of the enzyme, being similar to R206C. These findings permit us to assume that, although only R206 exhibits polymorphisms in a Caucasian‐specific manner, at least SNPs G82R and R206C in DNASE1L3 might be potential risk factors for autoimmune disease.     

Five non‐synonymous SNPs in the gene encoding human deoxyribonuclease I‐like 2 implicated in terminal differentiation of keratinocytes reduce or abolish its activity

Several non‐synonymous SNPs in the human deoxyribonuclease I‐like 2 (DNase 1L2) gene responsible for DNA degradation during terminal differentiation of epidermal keratinocytes have been identified. However, only limited population data are available, and furthermore the effect of these SNPs on the DNase 1L2 activity remains unknown. Genotyping of all of the 17 SNPs was performed using the PCR‐RFLP method in three ethnic groups including 14 different populations. A series of amino acid‐substituted DNase 1L2 corresponding to each SNP was expressed, and its activity was measured. All of the six non‐synonymous SNPs exhibited a mono‐allelic distribution, whereas the distribution of some SNPs other than exonic ones was ethnicity‐dependent. Each of the minor alleles in SNPs, p.Ala20Asp, p.Val104Leu, p.Asp197Ala, p.Glu274Lys and p.Asp287Asn, among the non‐synonymous SNPs produced low or no activity‐harbouring DNase 1L2. DNase 1L2 is well conserved, retaining full levels of enzymatic activity, with regard to these exonic SNPs in human populations. It seems plausible to assume that these SNPs affecting the activity may be one of the factors responsible for a genetic pre‐disposition for failure of differentiation‐associated cell death in various keratinocyte lineages, thereby leading to the development of parakeratosis. Our results may have clinical implications in relation to the pathogenesis of parakeratosis.     

Identification of the polymorphisms in IFITM3 gene and their association in a Korean population with ulcerative colitis

Interferons play critical roles in tumor pathogenesis by controlling apoptosis and through cellular anti-proliferative and differentiation activities. Interferon inducible transmembrane protein (IFITM) family genes have been implicated in several cellular processes such as the homotypic cell adhesion functions of IFN and cellular anti-proliferative activities. Expression levels of IFITM genes have been found to be up-regulated in gastric cancer cells and colorectal tumors. IFITM3 (also known as 1-8U) is a member of the IFITM family, and has been described as a key player in specification of germ cell fate. IFITM3 was first isolated from a genetic screen aimed at identifying genes involved in acquisition of germ cell competence. It has been proposed that epiblast cells have the highest expression of IFITM3 initiated germ cell specification and that homotypic association can discriminate germ cells from their somatic neighbors. In an attempt to better understand the genetic influences of IFITM3 on ulcerative colitis, we have identified possible variation sites and single nucleotide polymorphisms (SNPs) through two exons and their boundary IFITM3 intron sequences including the ~2.1 kb promoter regions. To determine whether or not these IFITM3 SNPs are associated with susceptibility to ulcerative colitis, frequencies of the genotype and allele of IFITM3 polymorphisms were analyzed on genomic DNAs isolated from patients with ulcerative colitis and from healthy controls. We also investigated the haplotype frequencies constructed by these SNPs in both groups. In this study, we also showed that expression level of IFITM3 mRNA was significantly higher in tissues of the ileum and cecum of the digestive system. We identified a total of seven SNPs and multiple variation regions in the IFITM3 gene. The genotype frequency of the g.-204T>G polymorphism in patients with ulcerative colitis was significantly different from that of the control group. Our results strongly suggest that polymorphisms of the IFITM3 gene may be associated with susceptibility to ulcerative colitis.