In silico analysis of functional nsSNPs in human TRPC6 gene associated with steroid resistant nephrotic syndrome

An inactivating human TRPC6 channel mutation without focal segmental glomerulosclerosis

Transient receptor potential cation channel-6 (TRPC6) gene mutations cause familial focal segmental glomerulosclerosis (FSGS), which is inherited as an autosomal dominant disease. In patients with TRPC6-related FSGS, all mutations map to the N- or C-terminal TRPC6 protein domains. Thus far, the majority of TRPC6 mutations are missense resulting in increased or decreased calcium influx; however, the fundamental molecular mechanisms causing cell injury and kidney pathology are unclear. We report a novel heterozygous TRPC6 mutation (V691Kfs*) in a large kindred with no signs of FSGS despite a largely truncated TRPC6 protein. We studied the molecular effects of V691Kfs* TRPC6 mutant using the tridimensional cryo-EM structure of the tetrameric TRPC6 protein. The results indicated that V691 is localized at the pore-forming transmembrane region affecting the ion conduction pathway, and predicted that V691Kfs* causes closure of the ion-conducting pathway leading to channel inactivation. We assessed the impact of V691Kfs* and two previously reported TRPC6 disease mutants (P112Q and G757D) on calcium influx in cells. Our data show that the V691Kfs* fully inactivated the TRCP6 channel-specific calcium influx consistent with a complete loss-of-function phenotype. Furthermore, the V691Kfs* truncation exerted a dominant negative effect on the full-length TRPC6 proteins. In conclusion, the V691Kfs* non-functional truncated TRPC6 is not sufficient to cause FSGS. Our data corroborate recently characterized TRPC6 loss-of-function and gain-of-function mutants suggesting that one defective TRPC6 gene copy is not sufficient to cause FSGS. We underscore the importance of increased rather than reduced calcium influx through TRPC6 for podocyte cell death.

Insilico prediction and functional analysis of nonsynonymous SNPs in human CTLA4 gene

The CTLA4 receptor is an immune checkpoint involved in the downregulation of T cells. Polymorphisms in this gene have been found to be associated with different diseases like rheumatoid arthritis, autosomal dominant immune dysregulation syndrome, juvenile idiopathic arthritis and autoimmune Addison's disease. Therefore, the identification of polymorphisms that have an effect on the structure and function of CTLA4 gene is important. Here we identified the most damaging missense or non-synonymous SNPs (nsSNPs) that might be crucial for the structure and function of CTLA4 using different bioinformatics tools. These in silico tools included SIFT, PROVEAN, PhD-SNP, PolyPhen-2 followed by MutPred2, I-Mutant 2.0 and ConSurf. The protein structures were predicted using Phyre2 and I-TASSER, while the gene-gene interactions were predicted by GeneMANIA and STRING. Our study identified three damaging missense SNPs rs1553657429, rs1559591863 and rs778534474 in coding region of CTLA4 gene. Among these SNPs the rs1553657429 showed a loss of potential phosphorylation site and was found to be highly conserved. The prediction of gene-gene interaction showed the interaction of CTlA4 with other genes and its importance in different pathways. This investigation of damaging nsSNPs can be considered in future while studying CTLA4 related diseases and can be of great importance in precision medicine.

Structural modeling of human AKAP3 protein and in silico analysis of single nucleotide polymorphisms associated with sperm motility

AKAP3 is a member of the A-kinase anchoring proteins and it is a constituent of the sperm fibrous sheath. AKAP3 is needed for the formation of sperm flagellum structure, sperm motility, and male fertility. This study aims to model the AKAP3 tertiary structure and identify the probable impact of four mutations characterized in infertile men on the AKAP3 structure. The T464S, I500T, E525K, and I661T substitutions were analyzed using in silico methods. The secondary structure and three-dimensional model of AKAP3 were determined using PSI-BLAST based secondary structure prediction and Robetta servers. The TM-score was used to quantitatively measure the structural similarities between native and mutated models. All of the desired substitutions were classified as benign. I-Mutant results showed all of the substitutions decreased AKAP3 stability; however, the I500T and I661T were more effective. Superposition and secondary structure comparisons between native and mutants showed no dramatic deviations. Our study provided an appropriate model for AKAP3. Destabilization of AKAP3 caused by these substitutions did not appear to induce structural disturbances. As AKAP3 is involved in male infertility, providing more structural insights and the impact of mutations that cause protein functional diversity could elucidate the etiology of male fertility problems at molecular level.

Mutational landscape of TRPC6, WT1, LMX1B, APOL1, PTPRO, PMM2, LAMB2 and WT1 genes associated with Steroid resistant nephrotic syndrome

BACKGROUND

Nephrotic syndrome appears as a group of symptoms like proteinuria, edema and hyperlipidemia. Identification of monogenic forms revealed the physiology and pathogenesis of the SRNS.

METHODS AND RESULTS

We performed Illumina panel sequencing of seven genes in 90 Indian patients to determine the role of these genetic mutations in nephrotic syndrome prognosis. Samtool was used for variants calling, and SnpEff and Snpsift did variants annotation. Clinical significance and variant classification were performed by the ClinVar database. In SSNS and SRNS patients, we found 0.78% pathogenic and 3.41% likely pathogenic mutations. Pathogenic mutations were found in LAMB2, LMX1B and WT1 genes, while likely pathogenic mutations were found in (6/13) LAMB2, (2/13) LMX1B, (2/13) TRPC6, (2/13) PTPRO and (1/13) PMM2 genes. Approximately 46% likely pathogenic mutations were contributed to the LAMB2 gene in SSNS and SRNS patients. We also detect 30 VUS (variants of uncertain significance), which were found (17/30) pathogenic and (13/30) likely pathogenic by different prediction tools.

CONCLUSIONS

Multigene panels were used for genetic screening of heterogeneous disorders like nephrotic syndrome in the Indian population. We found pathogenic, likely pathogenic and certain VUS, which were responsible for the pathogenesis of the disease. Therefore, mutational analysis of SSNS and SRNS is necessary to avoid adverse effects of corticosteroids, modify the intensity of immunosuppressing agents, and prevent the disease's progression.

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.

Evaluating the effect of BDNF Val66Met polymorphism on complex formation with HAP1 and Sortilin1 via structural modeling

Brain derived neurotrophic factor (BDNF) has a critical role in the neurogenesis, differentiation, survival of the neurons, regulation of the appetite, and energy homeostasis. Two key proteins, Huntingtin associated protein-1 (HAP1) and sortilin1, regulate the intracellular trafficking and stabilization of the precursor proBDNF through interaction with its prodomain region and mark it for secretion. Evidence suggests that the most frequent single nucleotide polymorphism (SNP) of BDNF gene (rs6265) has been associated with different psychiatric, neurodegenerative and eating disorders. In this study, structural bioinformatics and molecular dynamics (MD) simulations were applied, in order to get precise insights into the impact of Val66Met polymorphism on the proBDNF structure and its interaction with HAP1 and Sortilin1. Homology modeling, structure validation, refinement and also protein-protein docking were performed using appropriate servers. The stability, the fluctuations and the compactness of protein complexes were measured by MD simulation parameters including root mean square deviation (RMSD), root mean square fluctuation (RMSF) and Radius of gyration (Rg), respectively. The mutant proBDNF complexes with HAP1 and Sortilin1 revealed higher RMSD and RMSF values and also variable R_g over time compared with wild-type proBDNF. These computational results indicated that, wild-type proBDNF possessed more stable structure in binding with HAP1 and Sortilin1 compared with its mutant form. Therefore, Val66Met SNP could be deleterious due to making structural changes. It may cause a decrease in proBDNF secretion, which could possibly lead to different psychiatric, neurodegenerative and eating disorders. Further experimental lab studies are required for a more accurate conclusion.

Highly deleterious variations in COX1, CYTB, SCG5, FK2, PRL and PGF genes are the potential adaptation of the immigrated African ostrich population

Because of variable inconvenient living conditions in some places around the world, it is difficult to collect reliable physiological data for ostriches. Therefore, this study aims to provide a comprehensive in silico insight for the nature of polymorphism of important genetic loci that are related to physiological and reproductive traits. Sixty-nine mature ostriches ranging over half of Iraq were screened. Six exonic genetic loci, including cytochrome c oxidase I (COX1), cytochrome b (CYTB), secretogranin V (SCG5), feather keratin 2-like (FK2), prolactin (PRL) and placenta growth factor (PGF) were genotyped by PCR-single stranded conformation polymorphism (SSCP). Thirty-six novel SNPs, including seventeen nonsynonymous (ns) SNPs, were observed. Several computational software programs were utilized to assess the extent of the nsSNPs on their corresponding proteins structure, function and stability. The results showed several deleterious functional and stability changes in almost all the proteins studied. The total severity of each missense mutation was evaluated and compared with other nsSNPs accumulatively. It is evident from the extensive cumulative in silico computation that both p.E34D and p.E60K in PGF have the highest deleterious effect. The cumulative predictions from the present study are an impressive guide for the genotypes of African ostriches, which bypassed the expensive protocols for wet laboratory screening, to identify the effects of variants. To the best of our knowledge, this is the first investigation of its kind on the analyses and prediction outcome of missense mutations in African ostrich populations. The highly deleterious nsSNPs in the placenta growth factor are possible adaptive mutations which might be associated with adaptation in extreme and new environments. The flow and protocol of the computational predictions can be extended for various wild animals to identify the molecular nature of adaptations.

Characterization of NPHS2 gene polymorphisms associated to steroid resistance nephrotic syndrome in Indian children

Nephrotic syndrome (NS) is the common glomerular disease in children. These children are treated with steroids, depending upon their behavior. They are either steroid sensitive (SSNS) or steroid resistant (SRNS). NPHS2 gene mutants are linked to the risk of autosomal recessive SRNS and in some cases to SSNS. The present study has been performed to screen single nucleotide polymorphisms (SNPs) of the NPHS2 gene in a group of 90 Indian children suffering with NS (30 SSNS, 30 SRNS and 30 Controls) by PCR method followed by direct exon sequencing. Effect of SNPs on fold expression changes at transcript level of podocin was checked using quantitative real time PCR (qRT-PCR). SNPs identified through sequencing helps to carry out in-silico analysis. Overall 17 SNPs were identified in NPHS2 gene where 6 were found novel. Three missense SNPs p.R299Q, p.P20L and p.G35D were also identified in this population where SNP, p.G35D was found novel. In addition to sequencing analysis, results of in silico analysis shows that a mutant with these three missense SNPs has least ligand binding efficiency compared to native model. Moreover the significant observation of this study included two intronic SNPs c.451+23C>T and c.451+58A>T present in SRNS group of patients. These SNPs has shown high level of clinical significance within genomic and allelic frequency along with haplotypes and linkage disequilibrium count. The qRT-PCR analysis shows, down expression of podocin protein at transcript level in SRNS patients compared to SSNS patients. All these results support the fact that SNPs present in this population could affect the protein structural stability. Thus it is concluded that the polymorphisms predicted in this study might be disease causing in the NPHS2 gene and may have influence on the therapeutic response of NS patients.

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.

In Silico Analysis of nsSNPs of Carp TLR22 Gene Affecting its Binding Ability with Poly I:C

Immune response mediated by toll-like receptor 22 (TLR22), only found in teleost/amphibians, is triggered by double-stranded RNA binding to its LRR (leucine-rich repeats) ecto-domain. Accumulated evidences suggested that missense mutations in TLR genes affect its function. However, information on mutation linked pathogen recognition for TLR22 was lacking. The present study was commenced for predicting the effect of non-synonymous single-nucleotide polymorphisms (nsSNPs) on the pathogen recognizable LRR domain of TLR22 of farmed carp, Labeo rohita. The sequence-based algorithms (SIFT, PROVEAN and I-Mutant2.0) indicated that three SNPs (out of 27) such as p.L159F (rs76759876) and p.L529P (rs749355507) of LRR, and p.I836M (rs750758397) of intracellular motifs could potentially disrupt protein function. The 3D structure was generated using MODELLER 9.13 and further validated by SAVEs server. The simulated molecular docking of native TLR22 and mutants with poly I:C ligand indicated that mutations positioned at p.L159F and p.L529P of the LRR region affects the binding affinity significantly. This is the first kind of study of predicting nsSNPs of teleost TLR22 with disturbed ligand binding affinity with its extra-cellular LRR domain and thereby likely hindrance in subsequent signal transduction. This study serves as a guide for in vivo evaluation of impact of mutation on immune response mediated by teleost TLR22 gene.

Somatic SNPs of the BRCA2 gene at the fragments encoding RAD51 binding sites of canine mammary tumors

Mammary tumors are the most common tumor type both in women and in female dogs. In women, heritable breast cancers have been linked mutations in the breast cancer susceptibility gene BRCA2 and it contains eight BRC repeats in exon 11 that bind to RAD51. In this study, we investigated the sequence variations of BRC1-BRC8 and C-terminus of canine BRCA2 gene. From a total of 64 canine patients with mammary tumors, 31 mammary tumors with benign and malign carcinomas and the 3 normal mammary glands were used for the study. In this study, 19 SNPs of exon 11 of BRCA2 in canine mammary tumors were detected for the first time. The c.2383A>C (T1425P) SNP was found to be the most probable disease-associated nsSNP. Our findings suggest that T1425P variation in BRC3 to be the most probable disease-associated nsSNP and may affect RAD51 binding strength.