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

Structural impact, ligand-protein interactions, and molecular phenotypic effects of TGF-β1 gene variants: In silico analysis with implications for idiopathic pulmonary fibrosis

BACKGROUND

Idiopathic Pulmonary Fibrosis (IPF) is a chronic interstitial lung disease resulting in progressively deteriorating lung function. Transforming growth factor-β1 (TGF-β1) belongs to the TGF superfamily and exerts a profibrotic role in promoting lung fibrosis by facilitating fibroblast infiltration and activity, extracellular matrix deposition, and inhibition of collagen breakdown, thus promoting tissue remodelling and IPF.

MATERIALS AND METHODS

We evaluated the link between pathogenic TGF-β1 SNPs and IPF pathogenesis and the structure-activity functional consequences of those SNPs on the TGF-β1 protein. Several computational algorithms were merged to address the functional consequences of TGF-β1 gene mutations to protein stability, putative post-translational modification sites, ligand-protein interactions, and molecular phenotypic effects. These included FATHMM, POLYPHEN2, PROVEAN, and SIFT tools (identifying deleterious nsSNPs in the TGF-β1 gene), along with Pmut, PhD-SNP, SNAP, MutPred and the related TMHMM, MARCOIL, and DisProt algorithms (predicting structural disorders). INPS-MD was also used to evaluate the mutation-induced TGF-β1 protein's stability and MODPRED for recognition of post-translational TGF-β1 modification.

RESULTS

In total, 14 major pathogenic variants markedly impact the destabilization of the TGF-β1 protein, with most of these high-risk mutations associated with decreased stability of the TGF-β1 protein as per the I-Mutant, MUpro, and INPS-MD tools. R205W, R185W, R180Q, D86Y, and I300T variants were proposed to participate in the post-translational modifications, thus affecting affect protein-ligand interactions. Furthermore, at-risk genetic variants appear to target conserved regions in the alpha helices, random coils, and extracellular loops, resulting in a varied composition of amino acids, charge, hydrophobicity, and spatial architecture.

CONCLUSIONS

This study manuscript comprehensively analyzes gene variants within the TGF-β1 gene, offering novel insights into their structural and functional implications in interacting with target sites. This study is significant for the development of targeted therapeutic strategies and personalized treatment approaches for patients with inflammatory lung diseases such as IPF.

Computational insights into NIMA-related kinase 6: unraveling mutational effects on structure and function

The NEK6 (NIMA-related kinase 6) serine/threonine kinase is a pivotal player in a multitude of cellular processes, including the regulation of the cell cycle and the response to DNA damage. Its significance extends to disease pathogenesis, as changes in NEK6 activity have been linked to the development of cancer. Non-synonymous single nucleotide polymorphisms (nsSNPs) in NEK6 have been linked to cancer as they alter the protein's native structure and function. The association between NEK6 activity and cancer development has prompted researchers to explore the effects of genetic variations within the NEK6 gene. Therefore, we utilized advanced computational tools to analyze 155 high-confidence nsSNPs in the NEK6 gene. From this analysis, 21 nsSNPs were identified as potentially harmful, raising concerns about their impact on NEK6 activity and cancer risk. These 21 mutations were then examined for structural alterations, and eight of nsSNPs (I51M, V76A, I134N, Y152D, R171Q, V186G, L237R, and C285S) were found to destabilize the protein. Among the destabilizing mutations screened, a specific mutation, R171Q, stood out due to its conserved nature. To understand its impact on the protein and conformation, all-atom molecular dynamics simulations (MDS) for 100 ns were performed for both Wildtype NEK6 (WT-NEK6) and R171Q. The simulations revealed that the R171Q variant was unstable and led to significant conformational changes in NEK6. This study provides valuable insights into NEK6 dysfunction caused by single amino acid alterations, offering a novel understanding of the molecular mechanisms underlying NEK6-related cancer progression.

Exploring the Structural and Functional Effects of Nonsynonymous SNPs in the Human Serotonin Transporter Gene Through In Silico Approaches

The sodium-dependent serotonin transporter SLC6A4 (solute carrier family 6 member 4) gene encodes an intrinsic membrane protein that transmits the serotonin neurotransmitter from synaptic clefts into presynaptic neurons. The product of the SLC6A4 gene is related to the regulation of mood and social behavior, sleep, appetite, memory, digestion, and sexual desire. This protein is a target for antidepressant and psychostimulant drugs, thus prolonged neurotransmitter signaling remains blocked. In this study, the functional consequences of nsSNPs in the human SLC6A4 gene were explored through computational tools: PhD-SNP, SIFT, Align GVGD, PROVEAN, PMut, nsSNP Analyzer, SNPs&GO, SNAP2, PolyPhen2, and PANTHER to identify the most deleterious and damaging nsSNPs. Then the mutant protein stabilities were assessed using I-Mutant, MUpro, and MutPred2; amino acid conservation using ConSurf, and posttranslational modification analysis using MusiteDEEP and PROSPER. Furthermore, the 3-dimensional (3D) model of the mutated proteins was predicted and validated using SPARKS-X, Verify3D, and PROCHECK. The protein–ligand binding sites were analyzed using the COACH meta-server. Results from this study predicted that T192M, G342E, R607C, W282S, R104C, P131L, P156L, and N351S were the most structurally and functionally significant nsSNPs in the human SLC6A4 gene. Arg607 and Pro156 were the predicted sites for posttranslational modifications, and Thr192 and Try282 were the ligand-binding sites in the human SLC6A4 gene. The analyzed data also suggested that R104C, P131L, P156L, T192M, G342E, and W282S mutants might affect the binding of sodium ions with this protein. Taken together, this study provided important information on structurally and functionally important nsSNPs of the human SLC6A4 gene for further experimental validation. In the future, these damaging nsSNPs of the SLC6A4 gene have the potential to be evaluated as prognostic biomarkers for SLC6A4-related disorder diagnosis and research.

Scrutinizing Deleterious Nonsynonymous SNPs and Their Effect on Human POLD1 Gene

POLD1 (DNA polymerase delta 1, catalytic subunit) is a protein-coding gene that encodes the large catalytic subunit (POLD1/p125) of the DNA polymerase delta (Polδ) complex. The consequence of missense or nonsynonymous SNPs (nsSNPs), which occur in the coding region of a specific gene, is the replacement of single amino acid. It may also change the structure, stability, and/or functions of the protein. Mutation in the POLD1 gene is associated with autosomal dominant predisposition to colonic adenomatous polyps, colon cancer, endometrial cancer (EDMC), breast cancer, and brain tumors. These de novo mutations in the POLD1 gene also result in autosomal dominant MDPL syndrome (mandibular hypoplasia, deafness, progeroid features, and lipodystrophy). In this study, genetic variations of POLD1 which may affect the structure and/or function were analyzed using different types of bioinformatics tools. A total of 17038 nsSNPs for POLD1 were collected from the NCBI database, among which 1317 were missense variants. Out of all missense nsSNPs, 28 were found to be deleterious functionally and structurally. Among these deleterious nsSNPs, 23 showed a conservation scale of >5, 2 were predicted to be associated with binding site formation, and one acted as a posttranslational modification site. All of them were involved in coil, extracellular structures, or helix formation, and some cause the change in size, charge, and hydrophobicity.

Structural Consequence of Non-Synonymous Single-Nucleotide Variants in the N-Terminal Domain of LIS1

Disruptive neuronal migration during early brain development causes severe brain malformation. Characterized by mislocalization of cortical neurons, this condition is a result of the loss of function of migration regulating genes. One known neuronal migration disorder is lissencephaly (LIS), which is caused by deletions or mutations of the LIS1 (PAFAH1B1) gene that has been implicated in regulating the microtubule motor protein cytoplasmic dynein. Although this class of diseases has recently received considerable attention, the roles of non-synonymous polymorphisms (nsSNPs) in LIS1 on lissencephaly progression remain elusive. Therefore, the present study employed combined bioinformatics and molecular modeling approach to identify potential damaging nsSNPs in the LIS1 gene and provide atomic insight into their roles in LIS1 loss of function. Using this approach, we identified three high-risk nsSNPs, including rs121434486 (F31S), rs587784254 (W55R), and rs757993270 (W55L) in the LIS1 gene, which are located on the N-terminal domain of LIS1. Molecular dynamics simulation highlighted that all variants decreased helical conformation, increased the intermonomeric distance, and thus disrupted intermonomeric contacts in the LIS1 dimer. Furthermore, the presence of variants also caused a loss of positive electrostatic potential and reduced dimer binding potential. Since self-dimerization is an essential aspect of LIS1 to recruit interacting partners, thus these variants are associated with the loss of LIS1 functions. As a corollary, these findings may further provide critical insights on the roles of LIS1 variants in brain malformation.

Targeting the chemokine network in atherosclerosis

Chemokines and their receptors represent a potential target for immunotherapy in chronic inflammation. They comprise a large family of cytokines with chemotactic activity, and their cognate receptors are expressed on all cells of the body. This network dictates leukocyte recruitment and activation, angiogenesis, cell proliferation and maturation. Dysregulation of chemokine and chemokine receptor expression as well as function participates in many pathologies including cancer, autoimmune diseases and chronic inflammation. In atherosclerosis, a lipid-driven chronic inflammation of middle-sized and large arteries, chemokines and their receptors participates in almost all stages of the disease from initiation of fatty streaks to mature atherosclerotic plaque formation. Atherosclerosis and its complications are the main driver of mortality and morbidity in cardiovascular diseases (CVD). Hence, exploring new fields of therapeutic targeting of atherosclerosis is of key importance. This review gives an overview of the recent advances on the role of key chemokines and chemokine receptors in atherosclerosis, addresses chemokine-based biomarkers at biochemical, imaging and genetic level in human studies, and highlights the clinial trials targeting atherosclerosis.

Histological Analysis, Bioinformatics Profile, and Expression of Methylenetetrahydrofolate Reductase (MTHFR) in Bovine Testes

Simple Summary

To date, several genes have been sequenced but their corresponding protein characteristics remain unknown. This study highlights the histological structure of bovine (yellow-cattle and yak) testes as a build-up to exploring the bioinformatics profile and expression characteristics of methylenetetrahydrofolate reductase (MTHFR) in bovine testes. Our findings suggest that yellow-cattle testis have similar anatomical characteristics with that of yak, except for the weight or size, for which that of yellow-cattle is significantly higher or greater than yak. We also found that the secondary and 3D protein structures of MTHFR were similar to that of humans, with differences in the number of nucleotides, amino acids, and some physico-chemical characteristics. Moreover, MTHFR mRNA expression was higher in adult yellow-cattle and yak compared to their juvenile ones, however, its protein expression was higher but not statistically significant in adult yellow-cattle and yak compared to the juvenile ones. This provides a basis for further investigations into the regulatory function of MTHFR in bovine testes.

Abstract

Methylenetetrahydrofolate reductase (MTHFR), an enzyme expressed in mammalian testes, exerts a direct effect on spermatogenesis; however, its protein characteristics in bovine testes remain unknown. Here, we analysed bovine testicular structure, MTHFR bioinformatics profile, mRNA, and protein expression characteristics in yellow-cattle (y-c) and yak testis using histological procedures, bioinformatics analysis, qRT-PCR, and western blot. Testes from 13 bovines, ≤2 years juvenile (y-c, n = 3; yak, n = 3) and ≥3 years adult (y-c, n = 3; yak, n = 4) were collected and analysed. Anatomical characteristics of testis in y-c and yak were similar except the weight or size for which that of y-c was significantly higher or greater than yak. In y-c, an open reading frame (ORF) for 2600 nucleotides sequence, encoding 655 amino acids showed high homology with zebu cattle (99.51%) and wild yak (98.68%). Secondary and 3D protein structures were similar to that of humans with differences in the number of nucleotides, amino acids, and some physico-chemical characteristics. MTHFR mRNA expression in y-c and yak were significantly higher in adult testes compared with juvenile ones. However, its protein expression was higher, but not statistically significant, in adult y-c and yak compared to the juvenile ones. The highlights and inferences of these and other findings are discussed.

Computational analysis of functional single nucleotide polymorphisms associated with SLC26A4 gene

Single Nucleotide Polymorphisms (SNPs) are the most common candidate mutations in human beings that play a vital role in the genetic basis of certain diseases. Previous studies revealed that Solute Carrier Family 26 Member 4 (SLC26A4) being an essential gene of the multi-faceted transporter family SLC26 facilitates reflexive movement of Iodide into follicular lumen through apical membrane of thyrocyte. SLC26A4 gene encodes Pendred protein, a membrane glycoprotein, highly hydrophobic in nature, present at the apical membrane of thyrocyte functioning as transporter of iodide for thyroid cells. A minor genetic variation in SLC26A4 can cause Pendred syndrome, a syndrome associated with thyroid glands and deafness. In this study, we performed in-silico analysis of 674 missense SNPs of SLC26A4 using different computational platforms. The bunch of tools including SNPNEXUS, SNAP-2, PhD-SNP, SNPs&GO, I-Mutant, ConSurf, and ModPred were used to predict 23 highly confident damaging and disease causing nsSNPs (G209V, G197R, L458P, S427P, Q101P, W472R, N392Y, V359E, R409C, Q235R, R409P, G139V, G497S, H723R, D87G, Y127H, F667C, G334A, G95R, S427C, R291W, Q383H and E384G) that could potentially alter the SLC26A4 gene. Moreover, protein structure prediction, protein-ligand docking and Molecular Dynamics simulation were performed to confirm the impact of two evident alterations (Y127H and G334A) on the protein structure and function.

Computational analysis of high-risk SNPs in human CHK2 gene responsible for hereditary breast cancer: A functional and structural impact

Nowadays CHK2 mutation is studied frequently in hereditary breast and ovarian cancer patients in addition to BRCA1/BRCA2. CHK2 is a tumor suppressor gene that encodes a serine/threonine kinase, also involved in pathways such as DNA repair, cell cycle regulation and apoptosis in response to DNA damage. CHK2 is a well-studied moderate penetrance gene that correlates with third high risk susceptibility gene with an increased risk for breast cancer. Hence before planning large population study, it is better to scrutinize putative functional SNPs of CHK2 using different computational tools. In this study, we have used various computational approaches to identify nsSNPs which are deleterious to the structure and/or function of CHK2 protein that might be causing this disease. Computational analysis was performed by different in silico tools including SIFT, Align GVGD, SNAP-2, PROVEAN, Poly-Phen-2, PANTHER, PhD-SNP, MUpro, iPTREE-STAB, Consurf, InterPro, NCBI Conserved Domain Search tool, ModPred, SPARKS-X, RAMPAGE, Verify-3D, FT Site, COACH and PyMol. Out of 78 nsSNP of human CHK2 gene, seven nsSNPs were predicted functionally most significant SNPs. Among these seven nsSNP, p.Arg160Gly, p.Gly210Arg and p.Ser415Phe are highly conserved residues with conservation score of 9 and three nsSNP were predicted to be involved in post translational modification. The p.Arg160Gly and p.Gly210Arg may interfere in phosphopeptide binding site on FHA conserved domain. The p.Ser415Phe may interfere in formation of activation loop of protein-kinase domain and might interfere in interactions of CHK2 with ligand. The study concludes that mutation of serine to phenylalanine at position 415 is a major mutation in native CHK2 protein which might contribute to its malfunction, ultimately causing disease. This is the first comprehensive study, where CHK2 gene variants are analyzed using in silico tools hence it will be of great help while considering large scale studies and also in developing precision medicines related to these polymorphisms in the era of personalized medicine.

Computational and artificial neural network based study of functional SNPs of human LEPR protein associated with reproductive function

Genetic polymorphisms are mostly associated with inherited diseases, detecting and analyzing the biological significance of functional single-nucleotide polymorphisms (SNPs) using wet laboratory experiments is an arduous task hence the computational analysis of putative SNPs is essential before conducting a study on a large population. SNP in the leptin receptor (LEPR) could result in the retention of intracellular signalling due to the structural and functional instability of the receptor causing abnormal reproductive function in human. In this first comprehensive computational analysis of LEPR gene mutation, we have identified and analyzed the functional consequence and structural significance of the SNPs in LEPR using recently developed several computational algorithms. Thirteen deleterious mutations such as W13C, S93G, I232R, Q307H, Y354C, E497A, Q571H, R612H, K656N, T690A, T699M V741M, and L760R were identified in the LEPR gene coding region. Backpropagation algorithm has been developed to forestall the deleterious nature of SNP and to validate the outcome of the tested computational tools. From ConSurf prediction three SNPs (Q571H, R612H, and T699M) were highly conserved on LEPR protein and the most deleterious variant R612H had one hydrogen bond abolished and severely reduced protein stability. Molecular docking suggested that the mutant (R612H) LEPR had lowest binding energy than native LEPR with the ligand molecule. Thus the energetically destructive changeover of ARG to HIS in R612H could possibly affect the LEPR protein structural stability and functional constancy due to interruption in the amino acid interactions and could result in reproductive disorders in human and increases the complication in obstetric and pregnancy outcome.

Identification of deleterious SNPs and their effects on BCL11A, the master regulator of fetal hemoglobin expression

The B-cell lymphoma/leukemia 11A protein (encoded by BCL11A gene) is a key regulator of fetal-to-adult hemoglobin switching as seen in post-natal life. Although genetic polymorphisms like SNPs in BCL11A gene are expected to affect fetal hemoglobin (HbF) expression levels, yet their implications are poorly studied. This study utilizes a computational approach to identify the deleterious SNPs which may affect the structure and function of BCL11A protein. The study also generated a 3D structure of native and mutants. The analysis identified two SNPs in BCL11A as highly deleterious: N391K and C414S which are expected to affect structure and stability of the protein. According to conservation analysis, both residues N391 and C414 were identified as highly conserved. Additionally, post-translational modification sites were predicted at both sites. Ligand binding sites were also predicted in N391 and C414. Therefore, N391K and C414S in BCL11A can considered as important candidates to mediate HbF variation.

In Silico Analysis of Coding/Noncoding SNPs of Human RETN Gene and Characterization of Their Impact on Resistin Stability and Structure

Resistin (RETN) is a gene coding for proinflammatory adipokine called resistin secreted by macrophages in humans. Single nucleotide polymorphisms (SNPs) in RETN are linked to obesity and insulin resistance in various populations. Using dbSNP, 78 nonsynonymous SNPs (nsSNPs) were retrieved and tested on a PredictSNP 1.0 megaserver. Among these, 15 nsSNPs were predicted as highly deleterious and thus subjected to further analyses, such as conservation, posttranscriptional modifications, and stability. The 3D structure of human resistin was generated by homology modeling using Swiss model. Root-mean-square deviation (RMSD), hydrogen bonds (h-bonds), and interactions were estimated. Furthermore, UTRscan served to identify UTR functional SNPs. Among the 15 most deleterious nsSNPs, 13 were predicted to be highly conserved including variants in posttranslational modification sites. Stability analysis predicted 9 nsSNPs (I32S, C51Y, G58E, G58R, C78S, G79C, W98C, C103G, and C104Y) which can decrease protein stability with at least three out of the four algorithms used in this study. These nsSNPs were chosen for structural analysis. Both variants C51Y and C104Y showed the highest RMS deviations (1.137 Å and 1.308 Å, respectively) which were confirmed by the important decrease in total h-bonds. The analysis of hydrophobic and hydrophilic interactions showed important differences between the native protein and the 9 mutants, particularly I32S, G79C, and C104Y. Six SNPs in the 3'UTR (rs920569876, rs74176247, rs1447199134, rs943234785, rs76346269, and rs78048640) were predicted to be implicated in polyadenylation signal. This study revealed 9 highly deleterious SNPs located in the human RETN gene coding region and 6 SNPs within the 3'UTR that may alter the protein structure. Interestingly, these SNPs are worth to be analyzed in functional studies to further elucidate their effect on metabolic phenotype occurrence.