Functional and Structural Consequences of Damaging Single Nucleotide Polymorphisms in Human Prostate Cancer Predisposition Gene RNASEL

ASCARIS: Positional feature annotation and protein structure-based representation of single amino acid variations

Background

Genomic variations may cause deleterious effects on protein functionality and perturb biological processes. Elucidating the effects of variations is critical for developing novel treatment strategies for diseases of genetic origin. Computational approaches have been aiding the work in this field by modeling and analyzing the mutational landscape. However, new approaches are required, especially for accurate representation and data-centric analysis of sequence variations.

Method

In this study, we propose ASCARIS (Annotation and StruCture-bAsed RepresentatIon of Single amino acid variations), a method for the featurization (i.e., quantitative representation) of single amino acid variations (SAVs), which could be used for a variety of purposes, such as predicting their functional effects or building multi-omics-based integrative models. ASCARIS utilizes the direct and spatial correspondence between the location of the SAV on the sequence/structure and 30 different types of positional feature annotations (e.g., active/lipidation/glycosylation sites; calcium/metal/DNA binding, inter/transmembrane regions, etc.), along with structural features and physicochemical properties. The main novelty of this method lies in constructing reusable numerical representations of SAVs via functional annotations.

Results

We statistically analyzed the relationship between these features and the consequences of variations and found that each carries information in this regard. To investigate potential applications of ASCARIS, we trained variant effect prediction models that utilize our SAV representations as input. We carried out an ablation study and a comparison against the state-of-the-art methods and observed that ASCARIS has a competing and complementary performance against widely-used predictors. ASCARIS can be used alone or in combination with other approaches to represent SAVs from a functional perspective. ASCARIS is available as a programmatic tool at https://github.com/HUBioDataLab/ASCARIS and as a web-service at https://huggingface.co/spaces/HUBioDataLab/ASCARIS.

Non-synonymous SNPs variants of PRKCG and its association with oncogenes predispose to hepatocellular carcinoma

BACKGROUND

PRKCG encodes PKC γ, which is categorized under the classical protein kinase C family. No studies have specifically established the relationship between PRKCG nsSNPs with structural and functional variations in PKC γ in the context of hepatocellular carcinoma (HCC). The present study aims to uncover this link through in-silico and experimental studies.

METHODS

The 3D structure of PKC γ was predicted. Molecular Dynamic (MD) Simulations were run and estimates were made for interactions, stability, conservation and post-translational alterations between wild and mutant structures. The association of PRKCG levels with HCC survival rate was determined. Genotyping analyses were conducted to investigate the deleterious PRKCG nsSNP association with HCC. mRNA expression of PKC γ, HIF-1 alpha, AKT, SOCS3 and VEGF in the blood of controls and HCC patients was analyzed and a genetic cascade was constructed depicting these interactions.

RESULTS

The expression level of studied oncogenes was compared to tumour suppressor genes. Through Alphafold, the 3D structure of PKC γ was explored. Fifteen SNPs were narrowed down for in-silico analyses that were identified in exons 5, 10 and 18 and the regulatory and kinase domain of PKC γ. Root mean square deviation and fluctuation along with the radius of gyration unveiled potential changes between the wild and mutated variant structures. Mutant genotype AA (homozygous) corresponding to nsSNP, rs386134171 had more frequency in patients with OR (2.446), RR (1.564) and P-values (< 0.0029) that highlights its significant association with HCC compared to controls in which the wild genotype GG was found more prevalent.

CONCLUSION

nsSNP rs386134171 can be a genetic marker for HCC diagnosis and therapeutic studies. This study has laid down a road map for future studies to be conducted on HCC.

Impact of highly deleterious non-synonymous polymorphisms on GRIN2A protein's structure and function

GRIN2A is a gene that encodes NMDA receptors found in the central nervous system and plays a pivotal role in excitatory synaptic transmission, plasticity and excitotoxicity in the mammalian central nervous system. Changes in this gene have been associated with a spectrum of neurodevelopmental disorders such as epilepsy. Previous studies on GRIN2A suggest that non-synonymous single nucleotide polymorphisms (nsSNPs) can alter the protein's structure and function. To gain a better understanding of the impact of potentially deleterious variants of GRIN2A, a range of bioinformatics tools were employed in this study. Out of 1320 nsSNPs retrieved from the NCBI database, initially 16 were predicted as deleterious by 9 tools. Further assessment of their domain association, conservation profile, homology models, interatomic interaction, and Molecular Dynamic Simulation revealed that the variant I463S is likely to be the most deleterious for the structure and function of the protein. Despite the limitations of computational algorithms, our analyses have provided insights that can be a valuable resource for further in vitro and in vivo research on GRIN2A-associated diseases.

Computational screening and analysis of deleterious nsSNPs in human p14ARF (CDKN2A gene) protein using molecular dynamic simulation approach

Cyclin-dependent kinase inhibitor 2 A (CDKN2A) gene belongs to the cyclin-dependent kinase family that code for two transcripts (p16INK4A and p14ARF), both work as tumor suppressors proteins. The mutation that occurs in the p14ARF protein can lead to different types of cancers. Single nucleotide polymorphisms (SNPs) are an important type of genetic alteration that can lead to different types of diseases. In this study, we applied the computational strategy on human p14ARF protein to identify the potential deleterious nsSNPs and check their impact on the structure, function, and protein stability. We applied more than ten prediction tools to screen the retrieved 288 nsSNPs, consequently extracting four deleterious nsSNPs i.e., rs139725688 (R10G), rs139725688 (R21W), rs374360796 (F23L) and rs747717236 (L124R). Homology modeling, conservation and conformational analysis of mutant models were performed to examine the divergence of these variants from the native p14ARF structure. All-atom molecular dynamics simulation revealed a significant impact of these mutations on protein stability, compactness, globularity, solvent accessibility and secondary structure elements. Protein-protein interactions indicated that p14ARF operates as a hub linking clusters of different proteins and that changes in p14ARF may result in the disassociation of numerous signal cascades. Our current study is the first survey of computational analysis on p14ARF protein that determines the association of these nsSNPs with the altered function of p14ARF protein and leads to the development of various types of cancers. This research proposes the described functional SNPs as possible targets for proteomic investigations, diagnostic procedures, and treatments.Communicated by Ramaswamy H. Sarma.

Structural analysis and conformational dynamics of SOCS1 gene mutations involved in diffuse large B-cell lymphoma

OBJECTIVES

The SOCS1 gene is frequently mutated in primary Diffuse Large B-Cell Lymphoma (DLBCL) patients and is associated with a reduced survival rate. Using various computational techniques, the current study aims to identify Single Nucleotide Polymorphisms (SNPs) in the SOCS1 gene that are associated with the mortality rate of DLBCL patients. This study also evaluates the effects of SNPs on the structural instability of the SOCS1 protein in DLBCL patient.

METHODS

The cBioPortal webserver was used for mutations and determining how the SNP mutations affect the SOCS1 protein with various algorithms (PolyPhen-2.0, Provean, PhD-SNPg, SNPs&GO, SIFT, FATHMM, Predict SNP and SNAP). Five webservers (I-Mutant 2.0, MUpro, mCSM, DUET and SDM) were used for protein instability and the conserved status and were also predicted through different tools (ConSurf, Expasy, SOMPA). Lastly, MD simulations were run on the two chosen mutations (S116N and V128G) using GROMACS 5.0.1 to study how the mutations change the structure of SOCS1.

RESULTS

Among the 93 SOCS1 mutations detected in DLBCL patients, nine mutations were found to have a detrimental effect (damaging/deleterious/pathogenic/altered) on the SOCS1 protein. All the nine selected mutations are in the conserved region and four are on the extended strand site, four on the random coil site and one on the alpha helix position of the secondary protein structure. After anticipating the structural effects of these nine mutations, two were chosen (S116N and V128G) based on mutational frequency, location within the protein, structural effect (primary, secondary and tertiary) on stability and conservation status within the SOCS1 protein. The simulation of a 50 ns time interval revealed that the Rg value of S116N (2.17 nm) is higher than that of WT (1.98 nm), indicating a loss of structural compactness. In the case of the RMSD value, this mutated type (V128G) shows more deviation (1.54 nm) in comparison to the wild-type (2.14 nm) and another mutant type (S116N) (2.12 nm). The average RMSF values of wild-type and mutant types (V128G and S116N) were 0.88 nm, 0.49 nm, and 0.93 nm, respectively. The RMSF result shows that the mutant V128G structure is more stable than the wild-type and mutant S116N structures.

CONCLUSION

Based on all these computational predictions, this study finds that certain mutations, particularly S116N, have a destabilising and robust effect on the SOCS1 protein. These results can be used to learn more about the importance of SOCS1 mutations in DLBCL patients and to develop new ways to treat DLBCL.

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.

Pathogenic nsSNPs that increase the risks of cancers among the Orang Asli and Malays

Single-nucleotide polymorphisms (SNPs) are the most common genetic variations for various complex human diseases, including cancers. Genome-wide association studies (GWAS) have identified numerous SNPs that increase cancer risks, such as breast cancer, colorectal cancer, and leukemia. These SNPs were cataloged for scientific use. However, GWAS are often conducted on certain populations in which the Orang Asli and Malays were not included. Therefore, we have developed a bioinformatic pipeline to mine the whole-genome sequence databases of the Orang Asli and Malays to determine the presence of pathogenic SNPs that might increase the risks of cancers among them. Five different in silico tools, SIFT, PROVEAN, Poly-Phen-2, Condel, and PANTHER, were used to predict and assess the functional impacts of the SNPs. Out of the 80 cancer-related nsSNPs from the GWAS dataset, 52 nsSNPs were found among the Orang Asli and Malays. They were further analyzed using the bioinformatic pipeline to identify the pathogenic variants. Three nsSNPs; rs1126809 (TYR), rs10936600 (LRRC34), and rs757978 (FARP2), were found as the most damaging cancer pathogenic variants. These mutations alter the protein interface and change the allosteric sites of the respective proteins. As TYR, LRRC34, and FARP2 genes play important roles in numerous cellular processes such as cell proliferation, differentiation, growth, and cell survival; therefore, any impairment on the protein function could be involved in the development of cancer. rs1126809, rs10936600, and rs757978 are the important pathogenic variants that increase the risks of cancers among the Orang Asli and Malays. The roles and impacts of these variants in cancers will require further investigations using in vitro cancer models.

Emerging Diagnostic and Therapeutic Potentials of Human Hair Proteomics

The use of noninvasive human substrates to interrogate pathophysiological conditions has become essential in the post- Human Genome Project era. Due to its high turnover rate, and its long term capability to incorporate exogenous and endogenous substances from the circulation, hair testing is emerging as a key player in monitoring long term drug compliance, chronic alcohol abuse, forensic toxicology, and biomarker discovery, among other things. Novel high-throughput 'omics based approaches like proteomics have been underutilized globally in comprehending human hair morphology and its evolving use as a diagnostic testing substrate in the era of precision medicine. There is paucity of scientific evidence that evaluates the difference in drug incorporation into hair based on lipid content, and very few studies have addressed hair growth rates, hair forms, and the biological consequences of hair grooming or bleaching. It is apparent that protein-based identification using the human hair proteome would play a major role in understanding these parameters akin to DNA single nucleotide polymorphism profiling, up to single amino acid polymorphism resolution. Hence, this work seeks to identify and discuss the progress made thus far in the field of molecular hair testing using proteomic approaches, and identify ways in which proteomics would improve the field of hair research, considering that the human hair is mostly composed of proteins. Gaps in hair proteomics research are identified and the potential of hair proteomics in establishing a historic medical repository of normal and disease-specific proteome is also discussed.

Mutation-Structure-Function Relationship Based Integrated Strategy Reveals the Potential Impact of Deleterious Missense Mutations in Autophagy Related Proteins on Hepatocellular Carcinoma (HCC): A Comprehensive Informatics Approach

Autophagy, an evolutionary conserved multifaceted lysosome-mediated bulk degradation system, plays a vital role in liver pathologies including hepatocellular carcinoma (HCC). Post-translational modifications (PTMs) and genetic variations in autophagy components have emerged as significant determinants of autophagy related proteins. Identification of a comprehensive spectrum of genetic variations and PTMs of autophagy related proteins and their impact at molecular level will greatly expand our understanding of autophagy based regulation. In this study, we attempted to identify high risk missense mutations that are highly damaging to the structure as well as function of autophagy related proteins including LC3A, LC3B, BECN1 and SCD1. Number of putative structural and functional residues, including several sites that undergo PTMs were also identified. In total, 16 high-risk SNPs in LC3A, 18 in LC3B, 40 in BECN1 and 43 in SCD1 were prioritized. Out of these, 2 in LC3A (K49A, K51A), 1 in LC3B (S92C), 6 in BECN1 (S113R, R292C, R292H, Y338C, S346Y, Y352H) and 6 in SCD1 (Y41C, Y55D, R131W, R135Q, R135W, Y151C) coincide with potential PTM sites. Our integrated analysis found LC3B Y113C, BECN1 I403T, SCD1 R126S and SCD1 Y218C as highly deleterious HCC-associated mutations. This study is the first extensive in silico mutational analysis of the LC3A, LC3B, BECN1 and SCD1 proteins. We hope that the observed results will be a valuable resource for in-depth mechanistic insight into future investigations of pathological missense SNPs using an integrated computational platform.

Computational Analysis of High Risk Missense Variant in Human UTY Gene: A Candidate Gene of AZFa Sub-region

BACKGROUND

The human Ubiquitously transcribed tetratricopeptide repeat gene, Y-linked (UTY) gene encodes histone demethylase involved in protein-protein interactions. UTY protein evidence at protein level predicted intracellular and secreted protein. UTY is also involved in spermatogenesis process.

METHODS

The high-risk non-synonymous single nucleotide polymorphism in the coding region of the UTY gene was screened by SNP database and identified missense variants were subjected to computational analysis to understand the effect on protein function, stability and structure by SIFT, PolyPhen 2, PANTHER, PROVEAN, I-Mutant 2, iPTREE-STAB, ConSurf, ModPred, SPARKS-X, QMEAN, PROCHECK, project HOPE and STRING.

RESULTS

A total of 151 nsSNPs variants were retrieved in UTY gene out of which one missense variant (E18D) was predicted to be damaging or deleterious using SIFT, PolyPhen 2, PANTHER and PROVEAN. Additionally, E18D variant showed less stability, high conservation and having role in post translation modification using i-Mutant 2 and iPTREE-STAB, ConSurf and ModPred, respectively. The predicted 3D model of UTY using SPARKS-X with z-score of 15.16 was generated and validated via QMEAN (Z-score of 0.472) and PROCHECK which plots Ramachandran plot (85.3% residues in most favored regions, 12.3% in additionally allowed regions, 2.0% in generously allowed regions and 4.0% were in disallowed regions) and it indicates a good quality model. STRING showed that UTY interacts with ten different proteins.

CONCLUSION

This study revealed that SNP data available on database was deduced to find out the most damaging nsSNPs i.e. rs3212293 (E18D). Therefore, it provides useful information about functional SNPs for future prospects concerning infertility in men.