Predicting the most deleterious missense nsSNPs of the protein isoforms of the human HLA-G gene and in silico evaluation of their structural and functional consequences

Computational exploration of SLC14A1 genetic variants through structure modeling, protein-ligand docking, and molecular dynamics simulation

The urea transporter UT-B1, encoded by the SLC14A1 gene, has been hypothesized to be a significant protein whose deficiency and dysfunction contribute to the pathogenesis of bladder cancer and many other diseases. Several studies reported the association of genetic alterations in the SLC14A1 (UT-B1) gene with bladder carcinogenesis, suggesting a need for thorough characterization of the UT-B1 protein's coding and non-coding variants. This study used various computational techniques to investigate the commonly occurring germ-line missense and non-coding SNPs (ncSNPs) of the SLC14A1 gene (UT-B1) for their structural, functional, and molecular implications for disease susceptibility and dysfunctionality. SLC14A1 missense variants, primarily identified from the ENSEMBL genome browser, were screened through twelve functionality prediction tools leading to two variants D280Y (predicted detrimental by maximum tools) and D280N (high global MAF) for rs1058396. Subsequently, the ConSurf and NetSurf tools revealed the D280 residue to be in a variable site and exposed on the protein surface. According to I-Mutant2.0 and MUpro, both variants are predicted to cause a significant effect on protein stability. Analysis of molecular docking anticipated these two variants to decrease the binding affinity of UT-B1 protein for the examined ligands to a significant extent. Molecular dynamics also disclosed the possible destabilization of the UT-B1 protein due to single nucleotide polymorphism compared to wild-type protein which may result in impaired protein function. Furthermore, several non-coding SNPs were estimated to affect transcription factor binding and regulation of SLC14A1 gene expression. Additionally, two ncSNPs were found to affect miRNA-based post-transcriptional regulation by creating new seed regions for miRNA binding. This comprehensive in-silico study of SLC14A1 gene variants may serve as a springboard for future large-scale investigations examining SLC14A1 polymorphisms.

Investigating the Influence of ANTXR2 Gene Mutations on Protective Antigen Binding for Heightened Anthrax Resistance

Bacillus anthracis is the bacterium responsible for causing the zoonotic disease called anthrax. The disease presents itself in different forms like gastrointestinal, inhalation, and cutaneous. Bacterial spores are tremendously adaptable, can persist for extended periods and occasionally endanger human health. The Anthrax Toxin Receptor-2 (ANTXR2) gene acts as membrane receptor and facilitates the entry of the anthrax toxin into host cells. Additionally, mutations in the ANTXR2 gene have been linked to various autoimmune diseases, including Hyaline Fibromatosis Syndrome (HFS), Ankylosing Spondylitis (AS), Juvenile Hyaline Fibromatosis (JHF), and Infantile Systemic Hyalinosis (ISH). This study delves into the genetic landscape of ANTXR2, aiming to comprehend its associations with diverse disorders, elucidate the impacts of its mutations, and pinpoint minimal non-pathogenic mutations capable of reducing the binding affinity of the ANTXR2 gene with the protective antigen. Recognizing the pivotal role of single-nucleotide polymorphisms (SNPs) in shaping genetic diversity, we conducted computational analyses to discern highly deleterious and tolerated non-synonymous SNPs (nsSNPs) in the ANTXR2 gene. The Mutpred2 server determined that the Arg465Trp alteration in the ANTXR2 gene leads to altered DNA binding (p = 0.22) with a probability of a deleterious mutation of 0.808; notably, among the identified deleterious SNPs, rs368288611 (Arg465Trp) stands out due to its significant impact on altering the DNA-binding ability of ANTXR2. We propose these SNPs as potential candidates for hypertension linked to the ANTXR2 gene, which is implicated in blood pressure regulation. Noteworthy among the tolerated substitutions is rs200536829 (Ala33Ser), recognized as less pathogenic; this highlights its potential as a valuable biomarker, potentially reducing side effects on the host while also reducing binding with the protective antigen protein. Investigating these SNPs holds the potential to correlate with several autoimmune disorders and mitigate the impact of anthrax disease in humans.

Investigation of a novel TBC1D24 variation causing autosomal dominant non-syndromic hearing loss

Hearing loss is considered one of the most common sensory neurological defects, with approximately 60% of cases attributed to genetic factors. Human pathogenic variants in the TBC1D24 gene are associated with various clinical phenotypes, including dominant nonsyndromic hearing loss DFNA65, characterized by progressive hearing loss after the development of language. This study provides an in-depth analysis of the causative gene and mutations in a family with hereditary deafness. We recruited a three-generation family with autosomal dominant nonsyndromic hearing loss (ADNSHL) and conducted detailed medical histories and relevant examinations. Next-generation sequencing (NGS) was used to identify genetic variants in the proband, which were then validated using Sanger sequencing. Multiple computational software tools were employed to predict the impact of the variant on the function and structure of the TBC1D24 protein. A series of bioinformatics tools were applied to determine the conservation characteristics of the sequence, establish a three-dimensional structural model, and investigate changes in molecular dynamics. A detailed genotype and phenotype analysis were carried out. The family exhibited autosomal dominant, progressive, postlingual, and nonsyndromic sensorineural hearing loss. A novel heterozygous variant, c.1459C>T (p.His487Tyr), in the TBC1D24 gene was identified and confirmed to be associated with the hearing loss phenotype in this family. Conservation analysis revealed high conservation of the amino acid affected by this variant across different species. The mutant protein showed alterations in thermodynamic stability, elasticity, and conformational dynamics. Molecular dynamics simulations indicated changes in RMSD, RMSF, Rg, and SASA of the mutant structure. We computed the onset age of non-syndromic hearing loss associated with mutations in the TBC1D24 gene and identified variations in the hearing progression time and annual threshold deterioration across different frequencies. The identification of a new variant associated with rare autosomal dominant nonsyndromic hereditary hearing loss in this family broadens the range of mutations in the TBC1D24 gene. This variant has the potential to influence the interaction between the TLDc domain and TBC domain, thereby affecting the protein's biological function.

Unraveling Extremely Damaging IRAK4 Variants and Their Potential Implications for IRAK4 Inhibitor Efficacy

Interleukin-1-receptor-associated kinase 4 (IRAK4) possesses a crucial function in the toll-like receptor (TLR) signaling pathway, and the dysfunction of this molecule could lead to various infectious and immune-related diseases in addition to cancers. IRAK4 genetic variants have been linked to various types of diseases. Therefore, we conducted a comprehensive analysis to recognize the missense variants with the most damaging impacts on IRAK4 with the employment of diverse bioinformatics tools to study single-nucleotide polymorphisms' effects on function, stability, secondary structures, and 3D structure. The residues' location on the protein domain and their conservation status were investigated as well. Moreover, docking tools along with structural biology were engaged in analyzing the SNPs' effects on one of the developed IRAK4 inhibitors. By analyzing IRAK4 gene SNPs, the analysis distinguished ten variants as the most detrimental missense variants. All variants were situated in highly conserved positions on an important protein domain. L318S and L318F mutations were linked to changes in IRAK4 secondary structures. Eight SNPs were revealed to have a decreasing effect on the stability of IRAK4 via both I-Mutant 2.0 and Mu-Pro tools, while Mu-Pro tool identified a decreasing effect for the G198E SNP. In addition, detrimental effects on the 3D structure of IRAK4 were also discovered for the selected variants. Molecular modeling studies highlighted the detrimental impact of these identified SNP mutant residues on the druggability of the IRAK4 ATP-binding site towards the known target inhibitor, HG-12-6, as compared to the native protein. The loss of important ligand residue-wise contacts, altered protein global flexibility, increased steric clashes, and even electronic penalties at the ligand-binding site interfaces were all suggested to be associated with SNP models for hampering the HG-12-6 affinity towards IRAK4 target protein. This given model lays the foundation for the better prediction of various disorders relevant to IRAK4 malfunction and sheds light on the impact of deleterious IRAK4 variants on IRAK4 inhibitor efficacy.

Discovering Deleterious Single Nucleotide Polymorphisms of Human AKT1 Oncogene: An In Silico Study

BACKGROUND

AKT1 is a serine/threonine kinase necessary for the mediation of apoptosis, angiogenesis, metabolism, and cell proliferation in both normal and cancerous cells. The mutations in the AKT1 gene have been associated with different types of cancer. Further, the AKT1 gene mutations are also reported to be associated with other diseases such as Proteus syndrome and Cowden syndromes. Hence, this study aims to identify the deleterious AKT1 missense SNPs and predict their effect on the function and structure of the AKT1 protein using various computational tools.

METHODS

Extensive in silico approaches were applied to identify deleterious SNPs of the human AKT1 gene and assessment of their impact on the function and structure of the AKT1 protein. The association of these highly deleterious missense SNPs with different forms of cancers was also analyzed. The in silico approach can help in reducing the cost and time required to identify SNPs associated with diseases.

RESULTS

In this study, 12 highly deleterious SNPs were identified which could affect the structure and function of the AKT1 protein. Out of the 12, four SNPs-namely, G157R, G159V, G336D, and H265Y-were predicted to be located at highly conserved residues. G157R could affect the ligand binding to the AKT1 protein. Another highly deleterious SNP, R273Q, was predicted to be associated with liver cancer.

CONCLUSIONS

This study can be useful for pharmacogenomics, molecular diagnosis of diseases, and developing inhibitors of the AKT1 oncogene.

Identification of differentially expressed genes and SNPs linked to harvest body weight of genetically improved rohu carp, Labeo rohita

In most of the aquaculture selection programs, harvest body weight has been a preferred performance trait for improvement. Molecular interplay of genes linked to higher body weight is not elucidated in major carp species. The genetically improved rohu carp with 18% average genetic gain per generation with respect to harvest body weight is a promising candidate for studying genes' underlying performance traits. In the present study, muscle transcriptome sequencing of two groups of individuals, with significant difference in breeding value, belonging to the tenth generation of rohu carp was performed using the Illumina HiSeq 2000 platform. A total of 178 million paired-end raw reads were generated to give rise to 173 million reads after quality control and trimming. The genome-guided transcriptome assembly and differential gene expression produced 11,86,119 transcripts and 451 upregulated and 181 downregulated differentially expressed genes (DEGs) between high-breeding value and low-breeding value (HB & LB) groups, respectively. Similarly, 39,158 high-quality coding SNPs were identified with the Ts/Tv ratio of 1.23. Out of a total of 17 qPCR-validated transcripts, eight were associated with cellular growth and proliferation and harbored 13 SNPs. The gene expression pattern was observed to be positively correlated with RNA-seq data for genes such as myogenic factor 6, titin isoform X11, IGF-1 like, acetyl-CoA, and thyroid receptor hormone beta. A total of 26 miRNA target interactions were also identified to be associated with significant DETs (p-value < 0.05). Genes such as Myo6, IGF-1-like, and acetyl-CoA linked to higher harvest body weight may serve as candidate genes in marker-assisted breeding and SNP array construction for genome-wide association studies and genomic selection.

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.

Identification of mutations in porcine STAT5A that contributes to the transcription of CISH

Identification of causative genes or genetic variants associated with phenotype traits benefits the genetic improvement of animals. CISH plays a role in immunity and growth, however, the upstream transcriptional factors of porcine CISH and the genetic variations in these factors remain unclear. In this study, we firstly identified the minimal core promoter of porcine CISH and confirmed the existence of STATx binding sites. Overexpression and RT-qPCR demonstrated STAT5A increased CISH transcriptional activity (P < 0.01) and mRNA expression (P < 0.01), while GATA1 inhibited CISH transcriptional activity (P < 0.01) and the following mRNA expression (P < 0.05 or P < 0.01). Then, the putative functional genetic variations of porcine STAT5A were screened and a PCR-SSCP was established for genotype g.508A>C and g.566C>T. Population genetic analysis showed the A allele frequency of g.508A>C and C allele frequency of g.566C>T was 0.61 and 0.94 in Min pigs, respectively, while these two alleles were fixed in the Landrace population. Statistical analysis showed that Min piglets with CC genotype at g.566C>T or Hap1: AC had higher 28-day body weight, 35-day body weight, and ADG than TC or Hap3: CT animals (P < 0.05, P < 0.05). Further luciferase activity assay demonstrated that the activity of g.508A>C in the C allele was lower than the A allele (P < 0.05). Collectively, the present study demonstrated that STAT5A positively regulated porcine CISH transcription, and SNP g.566C>T in the STAT5A was associated with the Min piglet growth trait.

HBD-2 variants and SARS-CoV-2: New insights into inter-individual susceptibility

Background

A deep understanding of the causes of liability to SARS-CoV-2 is essential to develop new diagnostic tests and therapeutics against this serious virus in order to overcome this pandemic completely. In the light of the discovered role of antimicrobial peptides [such as human b-defensin-2 (hBD-2) and cathelicidin LL-37] in the defense against SARS-CoV-2, it became important to identify the damaging missense mutations in the genes of these molecules and study their role in the pathogenesis of COVID-19.

Methods

We conducted a comprehensive analysis with multiple in silico approaches to identify the damaging missense SNPs for hBD-2 and LL-37; moreover, we applied docking methods and molecular dynamics analysis to study the impact of the filtered mutations.

Results

The comprehensive analysis reveals the presence of three damaging SNPs in hBD-2; these SNPs were predicted to decrease the stability of hBD-2 with a damaging impact on hBD-2 structure as well. G51D and C53G mutations were located in highly conserved positions and were associated with differences in the secondary structures of hBD-2. Docking-coupled molecular dynamics simulation analysis revealed compromised binding affinity for hBD-2 SNPs towards the SARS-CoV-2 spike domain. Different protein-protein binding profiles for hBD-2 SNPs, in relation to their native form, were guided through residue-wise levels and differential adopted conformation/orientation.

Conclusions

The presented model paves the way for identifying patients prone to COVID-19 in a way that would guide the personalization of both the diagnostic and management protocols for this serious disease.

Analyzing PKC Gamma (+ 19,506 A/G) polymorphism as a promising genetic marker for HCV-induced hepatocellular carcinoma

BACKGROUND

HCC is a major health concern worldwide. PKC gamma, a member of the conventional PKC subclass, is involved in many cancer types, but the protein has received little attention in the context of single nucleotide polymorphisms and HCC. Therefore, the study aims to investigate the association of PKC gamma missense SNP with HCV-induced hepatocellular carcinoma.

METHODS

The PKC gamma nsSNPs were retrieved from the ENSEMBL genome browser and the deleterious nsSNPs were filtered out through involvingPredictSNP2, CADD, DANN, FATHMM, FunSeq2 and GWAVA. Among the filtered nsSNPs, nsSNP rs1331262028 was identified to be the most pathogenic one. Through involving I-TASSER, ProjectHOPE, I-Mutant, MUpro, mCSM, SDM, DynaMut and MutPred, the influence of SNP rs1331262028 on protein structure, function and stability was estimated. A molecular Dynamic simulation was run to determine the conformational changes in mutant protein structure compared to wild. The blood samples were collected for genotyping analysis and for assessing ALT levels in the blood.

RESULTS

The study identified for the first time an SNP (rs1331262028) of PRKCG to strongly decrease protein stability and induce HCC. The RMSD, RMSF, and Rg values of mutant and wild types found were significantly different. Based on OR and RR values of 5.194 and 2.287, respectively, genotype analysis revealed a higher correlation between the SNP homozygous wild Typeform, AA, and the disease while patients with genotype AG have higher viral load.

CONCLUSION

Outcomes of the current study delineated PKC gamma SNP rs1331262028 as a genetic marker for HCV-induced HCC that could facilitate disease management after further validation.

Computational Analysis of Deleterious SNPs in NRAS to Assess Their Potential Correlation With Carcinogenesis

The NRAS gene is a well-known oncogene that acts as a major player in carcinogenesis. Mutations in the NRAS gene have been linked to multiple types of human tumors. Therefore, the identification of the most deleterious single nucleotide polymorphisms (SNPs) in the NRAS gene is necessary to understand the key factors of tumor pathogenesis and therapy. We aimed to retrieve NRAS missense SNPs and analyze them comprehensively using sequence and structure approaches to determine the most deleterious SNPs that could increase the risk of carcinogenesis. We also adopted structural biology methods and docking tools to investigate the behavior of the filtered SNPs. After retrieving missense SNPs and analyzing them using six in silico tools, 17 mutations were found to be the most deleterious mutations in NRAS. All SNPs except S145L were found to decrease NRAS stability, and all SNPs were found on highly conserved residues and important functional domains, except R164C. In addition, all mutations except G60E and S145L showed a higher binding affinity to GTP, implicating an increase in malignancy tendency. As a consequence, all other 14 mutations were expected to increase the risk of carcinogenesis, with 5 mutations (G13R, G13C, G13V, P34R, and V152F) expected to have the highest risk. Thermodynamic stability was ensured for these SNP models through molecular dynamics simulation based on trajectory analysis. Free binding affinity toward the natural substrate, GTP, was higher for these models as compared to the native NRAS protein. The Gly13 SNP proteins depict a differential conformational state that could favor nucleotide exchange and catalytic potentiality. A further application of experimental methods with all these 14 mutations could reveal new insights into the pathogenesis and management of different types of tumors.

NELL-1 in Genome-Wide Association Studies across Human Diseases

Neural epidermal growth factor-like (EGFL)-like protein (NELL)-1 is a potent and key osteogenic factor in the development and regeneration of skeletal tissues. Intriguingly, accumulative data from genome-wide association studies (GWASs) have started unveiling potential broader roles of NELL-1 beyond its functions in bone and cartilage. With exploration of the genetic variants of the entire genome in large-scale disease cohorts, GWASs have been used for establishing the connection between specific single-nucleotide polymorphisms of NELL1, in addition to osteoporosis, metabolic diseases, inflammatory conditions, neuropsychiatric diseases, neurodegenerative disorders, and malignant tumors. This review summarizes the findings from GWASs on the manifestation, significance level, implications on function, and correlation of specific NELL1 single-nucleotide polymorphisms in various disorders in humans. By offering a unique and comprehensive correlation between genetic variants and plausible functions of NELL1 in GWASs, this review illustrates the wide range of potential effects of a single gene on the pathogenesis of multiple disorders in humans.

A Novel Mutation of the KLK6 Gene in a Family With Knee Osteoarthritis

To investigate the correlation between gene mutation and knee osteoarthritis (KOA), a whole-exome sequencing (WES) was applied to analyze blood samples of four KOA patients and two normal subjects in a family. Gene mutations were identified by gene-trapping and high-throughput sequencing analysis across the differences between the patients and normal subjects. The interactive gene network analysis on the retrieval of interacting genes (STRING) database and the KOA-related genes expression data sets was performed. A possibly detrimental and nonsynonymous mutation at the kallikrein-related peptidase 6 (KLK6) gene (rs201586262, c. C80A, P27H) was identified and attracted our attention. KLK6 belongs to the kallikrein family of serine proteases and its serum level is known as a prevalent biomarker in inflammatory and malignant diseases. KLK6 expresses in the extracellular compartment for matrix degradation, highlighting that KLK6 plays a role in the pathogenesis of KOA. By using the gene databases, the KOA-related genes were mined after de-duplication and IL6 was selected as the most relevant gene through interactive analysis of protein-protein interaction (PPI) network. The data suggested that KLK6 gene mutation and the related expression alteration of IL6 gene might determine the occurrence of hereditary KOA. The is the first study discovering the gene mutation of KLK6 as a factor of pathogenesis of KOA, especially the hereditary KOA.

Appropriate fixative for MEM-G/9 staining of cultured human HLA-G-positive JEG-3 trophoblast tumor cells

Human leukocyte antigen (HLA-G) participates in immunosuppression and is useful for prenatal diagnostics. Isolation of fetal cells positive for HLA-G by HLA-G antibody conjugated nanoparticles from the cervix of pregnant women is the basis for non-invasive prenatal testing. Endocervical specimens are fixed in transport medium before isolation using antibody conjugated nanoparticles. Staining of HLA-G using MEM-G/9 antibody, however, is restricted to unfixed cells. We investigated the effect of several fixatives on the interaction of HLA-G with MEM-G/9 in the HLA-G-positive cell line, JEG-3. We investigated absolute methanol, 1:1 acetate buffer:methanol, Pap solution and paraformaldehyde. The effects of these fixatives were evaluated using immunofluorescence. We found no MEM-G/9 surface staining of methanol fixed cells. Approximately 40% of JEG-3 cells fixed with paraformaldehyde failed to stain. Nearly all cells were stained with MEM-G/9 following fixation with acetate buffer:methanol or Pap solution. Our findings indicate the importance of using an appropriate fixative for preserving HLA-G cell surface antigen for studies using the MEM-G/9 antibody.