In silico analysis of nonsynonymous single-nucleotide polymorphisms (nsSNPs) of the SMPX gene

Computational Exploration of Single-Nucleotide Polymorphisms in the Human hRAS Gene: Implications and Insights

Background A group of genes called oncogenes includes the Harvey rat sarcoma virus (hRAS) gene. Along with hRAS, Kirsten rat sarcoma viral oncogene homolog (kRAS) and neuroblastoma RAS viral oncogene homolog (nRAS) genes belong to the Rat sarcoma (Ras) family of oncogenes. These three genes result in Rho guanosine triphosphate hydrolases (GTPases) as their protein product. Instructions for producing the protein hRAS, which is mainly involved in controlling cell division, are provided by the hRAS gene. The hRAS protein transfers signals from outside through a process called signal transduction. Because the hRAS protein is a GTPase, it changes the chemical guanosine-5'-triphosphate (GTP) into guanosine diphosphate (GDP). GTP and GDP molecules operate as switches to turn on and off the hRAS. This study aimed to anticipate the structure and stability of the protein resulting from missense single-nucleotide polymorphisms (SNPs) in the human hRAS genes. Methodology To investigate the possible negative effects associated with these SNPs, bioinformatic analysis is typically essential. The following tools were employed for forecasting harmful SNPs: Scale-Invariant Feature Transform (SIFT), Protein Analysis Through Evolutionary Relationships (PANTHER), non-synonymous SNP by Protein Variation Effect Analyzer (PROVEAN), and non-synonymous SNP by Single Nucleotide Polymorphism Annotation Platform (SNAP). Results The present study identified a total of 11 SNPs using the SIFT approach, which were shown to have functional significance. Only two of these 11 SNPs were determined to be tolerable, whereas nine were shown to be detrimental. Among the 11 SNPs analyzed, seven (Q61H, Q99H, K117R, A121D, A146V, R169W, R169Q) were classified as possibly damaging,and four (G13V, Q22K, Q61K, Q13V) were categorized as probably benign according to the predictions made by PANTHER tools. Therefore, the seven SNPs were identified as high-risk SNPs. Conclusions Given that SNPs have the potential to be candidates for cellular alterations brought on by mutations that are associated with cancer, this study provides vital information about how SNPs might be utilized as a diagnostic marker for cancer.

The pathogenic effect of SNPs on structure and function of human TLR4 using a computational approach

The human toll-like receptor (hTLR) 4 single nucleotide polymorphisms (SNPs) are interconnected with cancer, multiple genetic disorders and other immune-related diseases. The detrimental effect of SNPs in hTLR4 with respect to structure and function has not been explored in depth. The present study concatenates the biological consequences of the SNPs along with structural modifications predicted at the hTLR4 gene. A total of 7910 SNPs of hTLR4 were screened, and 21 damage-causing SNPs were identified. Out of 21, seven are present in the extracellular region, of which three were detected as deleterious and the fourth one as moderate. These three mutations are located in a highly conserved region and influence conformational change. The change leads to the widening of the Leucine-rich repeat (LRR) arc to a maximum of 16.9 Å and a minimum of 8.7 Å. Expansion/shortening of LRR arc, never discussed before, would cause loss of myeloid differentiation factor 2 (MD-2) interactions in the interior and diminish lipopolysaccharide (LPS) responses. Similarly, in all mutant structures, the binding region for HMGB1 and LPS is deflating or in an unsupportive conformation. Thus, SNPs affect the regular signaling cascade and might result in human sepsis, genetic disorders, cancer and other immunological related diseases.Communicated by Ramaswamy H. Sarma.

DeteX: A highly accurate software for detecting SNV and InDel in single and paired NGS data in cancer research

Background: Genetic testing is becoming more and more accepted in the auxiliary diagnosis and treatment of tumors. Due to the different performance of the existing bioinformatics software and the different analysis results, the needs of clinical diagnosis and treatment cannot be met. To this end, we combined Bayesian classification model (BC) and fisher exact test (FET), and develop an efficient software DeteX to detect SNV and InDel mutations. It can detect the somatic mutations in tumor-normal paired samples as well as mutations in a single sample. Methods: Combination of Bayesian classification model (BC) and fisher exact test (FET). Results: We detected SNVs and InDels in 11 TCGA glioma samples, 28 clinically targeted capture samples and 2 NCCL-EQA standard samples with DeteX, VarDict, Mutect, VarScan and GatkSNV. The results show that, among the three groups of samples, DeteX has higher sensitivity and precision whether it detects SNVs or InDels than other callers and the F1 value of DeteX is the highest. Especially in the detection of substitution and complex mutations, only DeteX can accurately detect these two kinds of mutations. In terms of single-sample mutation detection, DeteX is much more sensitive than the HaplotypeCaller program in Gatk. In addition, although DeteX has higher mutation detection capabilities, its running time is only .609 of VarDict, which is .704 and .343 longer than VarScan and MuTect, respectively. Conclusion: In this study, we developed DeteX to detect SNV and InDel mutations in single and paired samples. DeteX has high sensitivity and precision especially in the detection of substitution and complex mutations. In summary, DeteX from NGS data is a good SNV and InDel caller.

In silico design of epitope-based peptide vaccine against non-typhoidal Salmonella through immunoinformatic approaches

Non-typhoidal Salmonella (NTS) is one of the leading bacterial causes of many invasive human infections with a high antibiotic resistance profile. With this concern, the current study aimed to design an effective epitope-based peptide vaccine against NTS species as a successive and substitutive protective measure of invasive NTS disease. To design rationally, the current study considered a comprehensive in silico workflow combination of both immunoinformatics and molecular modeling approaches, including molecular docking and molecular dynamics (MD) simulation. We identified the two most promising T cell epitopes KVLYGIFAI and YGIFAITAL, and three B cell epitopes AAPVQVGEAAGS, TGGGDGSNT, and TGGGDGSNTGTTTT, in the outer membrane of NTS. Using these epitopes, a multiepitope vaccine was subsequently constructed along with appropriate adjuvant and linkers, which showed a good binding affinity and stability with toll-like receptor 2 (TLR2) in both molecular docking and MD simulation. Furthermore, in silico immune simulation described a strong immune response with a high number of antibodies, interferon-γ, and activated B and T cells. This study collectively suggests that predicted vaccine constructs could be considered potential vaccine candidates against common NTS species.Communicated by Ramaswamy H. Sarma.

Dynamic insights into the effects of nonsynonymous polymorphisms (nsSNPs) on loss of TREM2 function

Single nucleotide variations in Triggering Receptor Expressed on Myeloid Cells 2 (TREM2) are associated with many neurodegenerative diseases, including Nasu-Hakola disease (NHD), frontotemporal dementia (FTD), and late-onset Alzheimer's disease because they disrupt ligand binding to the extracellular domain of TREM2. However, the effects of nonsynonymous single nucleotide polymorphisms (nsSNPs) in TREM2 on disease progression remain unknown. In this study, we identified several high-risk nsSNPs in the TREM2 gene using various deleterious SNP predicting algorithms and analyzed their destabilizing effects on the ligand recognizing region of the TREM2 immunoglobulin (Ig) domain by molecular dynamics (MD) simulation. Cumulative prediction by all tools employed suggested the three most deleterious nsSNPs involved in loss of TREM2 function are rs549402254 (W50S), rs749358844 (R52C), and rs1409131974 (D104G). MD simulation showed that these three variants cause substantial structural alterations and conformational remodeling of the apical loops of the TREM2 Ig domain, which is responsible for ligand recognition. Detailed analysis revealed that these variants substantially increased distances between apical loops and induced conformation remodeling by changing inter-loop nonbonded contacts. Moreover, all nsSNPs changed the electrostatic potentials near the putative ligand-interacting region (PLIR), which suggested they might reduce specificity or loss of binding affinity for TREM2 ligands. Overall, this study identifies three potential high-risk nsSNPs in the TREM2 gene. We propose further studies on the molecular mechanisms responsible for loss of TREM2 function and the associations between TREM2 nsSNPs and neurodegenerative diseases.

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.

Identification of deleterious single nucleotide polymorphism (SNP)s in the human TBX5 gene & prediction of their structural & functional consequences: An in silico approach

T-box transcription factor 5 gene (TBX5) encodes the transcription factor TBX5, which plays a crucial role in the development of heart and upper limbs. Damaging single nucleotide variants in this gene alter the protein structure, disturb the functions of TBX5, and ultimately cause Holt-Oram Syndrome (HOS). By analyzing the available single nucleotide polymorphism information in the dbSNP database, this study was designed to identify the most deleterious TBX5 SNPs through insilico approaches and predict their structural and functional consequences.

Fifty-eight missense substitutions were found damaging by sequence homology-based tools: SIFT and PROVEAN, and structure homology-based tool PolyPhen-2. Various disease association meta-predictors further scrutinized these SNPs. Additionally, conservation profile of the amino acid residues, their surface accessibility, stability, and structural integrity of the native protein upon mutations were assessed. From these analyses, finally 5 SNPs were detected as the most damaging ones: [rs1565941579 (P85S), rs1269970792 (W121R), rs772248871 (V153D), rs769113870 (E208D), and rs1318021626 (I222N)]. Analyses of stop-lost, nonsense, UTR, and splice site SNPs were also conducted.

Through integrative bioinformatics analyses, this study has identified the SNPs that are deleterious to the TBX5 protein structure and have the potential to cause HOS. Further wet-lab experiments can validate these findings.

•

Deleterious SNPs in the human TBX5 gene responsible for Holt-Oram Syndrome have been identified.

•

58 missense and 2 nonsense SNPs were identified as deleterious.

•

86 3′ UTR SNPs were predicted to be located on miRNA target sites.

•

Possible effects of missense SNPs on the TBX5 protein structure have been studied.

Novel dynamic residue network analysis approaches to study allosteric modulation: SARS-CoV-2 Mpro and its evolutionary mutations as a case study

The rational search for allosteric modulators and the allosteric mechanisms of these modulators in the presence of mutations is a relatively unexplored field. Here, we established novel in silico approaches and applied them to SARS-CoV-2 main protease (Mpro) as a case study. First, we identified six potential allosteric modulators. Then, we focused on understanding the allosteric effects of these modulators on each of its protomers. We introduced a new combinatorial approach and dynamic residue network (DRN) analysis algorithms to examine patterns of change and conservation of critical nodes, according to five independent criteria of network centrality. We observed highly conserved network hubs for each averaged DRN metric on the basis of their existence in both protomers in the absence and presence of all ligands (persistent hubs). We also detected ligand specific signal changes. Using eigencentrality (EC) persistent hubs and ligand introduced hubs we identified a residue communication path connecting the allosteric binding site to the catalytic site. Finally, we examined the effects of the mutations on the behavior of the protein in the presence of selected potential allosteric modulators and investigated the ligand stability. One crucial outcome was to show that EC centrality hubs form an allosteric communication path between the allosteric ligand binding site to the active site going through the interface residues of domains I and II; and this path was either weakened or lost in the presence of some of the mutations. Overall, the results revealed crucial aspects that need to be considered in rational computational drug discovery.

Computational Insights into the Deleterious Impacts of Missense Variants on N-Acetyl-d-glucosamine Kinase Structure and Function

An enzyme of the mammalian amino-sugar metabolism pathway, N-acetylglucosamine kinase (NAGK), that synthesizes N-acetylglucosamine (GlcNAc)-6-phosphate, is reported to promote dynein functions during mitosis, axonal and dendritic growth, cell migration, and selective autophagy, which all are unrelated to its enzyme activity. As non-enzymatic structural functions can be altered by genetic variation, we made an effort in this study aimed at deciphering the pathological effect of nonsynonymous single-nucleotide polymorphisms (nsSNPs) in NAGK gene. An integrated computational approach, including molecular dynamics (MD) simulation and protein-protein docking simulation, was used to identify the damaging nsSNPs and their detailed structural and functional consequences. The analysis revealed the four most damaging variants (G11R, G32R, G120E, and A156D), which are highly conserved and functional, positioned in both small (G11R and G32R) and large (G120E and A156D) domains of NAGK. G11R is located in the ATP binding region, while variants present in the large domain (G120E and A156D) were found to induce substantial alterations in the structural organizations of both domains, including the ATP and substrate binding sites. Furthermore, all variants were found to reduce binding energy between NAGK and dynein subunit DYNLRB1, as revealed by protein-protein docking and MM-GBSA binding energy calculation supporting their deleteriousness on non-canonical function. We hope these findings will direct future studies to gain more insight into the role of these variants in the loss of NAGK function and their role in neurodevelopmental disorders.

Integrated System Pharmacology and In Silico Analysis Elucidating Neuropharmacological Actions of Withania somnifera in the Treatment of Alzheimer's Disease

BACKGROUND

Withania somnifera (WS), also referred to as Medhya Rasayana (nootropic or rejuvenating), has traditionally been prescribed for various neurological ailments, including dementia. Despite substantial evidence, pharmacological roles of WS, neither as nootropic nor as an antidementia agent, are well-understood at the cellular and molecular levels.

OBJECTIVES

We aimed at elucidating the pharmacological action mechanisms of WS root constituents against Alzheimer's Disease (AD) pathology.

METHODS

Various bioinformatics tools and resources, including DAVID, Cytoscape, NetworkAnalyst and KEGG pathway database were employed to analyze the interaction of WS root bioactive molecules with the protein targets of AD-associated cellular processes. We also used a molecular simulation approach to validate the interaction of compounds with selected protein targets.

RESULTS

Network analysis revealed that β-sitosterol, withaferin A, stigmasterol, withanolide A, and withanolide D are the major constituents of WS root that primarily target the cellular pathways such as PI3K/Akt signaling, neurotrophin signaling and toll-like receptor signaling and proteins such as Tropomyosin receptor Kinase B (TrkB), Glycogen Synthase Kinase-3β (GSK-3β), Toll-Like Receptor 2/4 (TLR2/4), and β-secretase (BACE-1). Also, the in silico analysis further validated the interaction patterns and binding affinity of the major WS compounds, particularly stigmasterol, withanolide A, withanolide D and β-sitosterol with TrkB, GSK-3β, TLR2/4, and BACE-1.

CONCLUSION

The present findings demonstrate that stigmasterol, withanolide A, withanolide D and β-sitosterol are the major metabolites that are responsible for the neuropharmacological action of WS root against AD-associated pathobiology, and TrkB, GSK-3β, TLR2/4, and BACE-1 could be the potential druggable targets.

Computational SNP Analysis and Molecular Simulation Revealed the Most Deleterious Missense Variants in the NBD1 Domain of Human ABCA1 Transporter

The ATP-binding cassette transporter A1 (ABCA1) is a membrane-bound exporter protein involved in regulating serum HDL level by exporting cholesterol and phospholipids to load up in lipid-poor ApoA-I and ApoE, which allows the formation of nascent HDL. Mutations in the ABCA1 gene, when presents in both alleles, disrupt the canonical function of ABCA1, which associates with many disorders related to lipid transport. Although many studies have reported the phenotypic effects of a large number of ABCA1 variants, the pathological effect of non-synonymous polymorphisms (nsSNPs) in ABCA1 remains elusive. Therefore, aiming at exploring the structural and functional consequences of nsSNPs in ABCA1, in this study, we employed an integrated computational approach consisting of nine well-known in silico tools to identify damaging SNPs and molecular dynamics (MD) simulation to get insights into the magnitudes of the damaging effects. In silico tools revealed four nsSNPs as being most deleterious, where the two SNPs (G1050V and S1067C) are identified as the highly conserved and functional disrupting mutations located in the NBD1 domain. MD simulation suggested that both SNPs, G1050V and S1067C, changed the overall structural flexibility and dynamics of NBD1, and induced substantial alteration in the structural organization of ATP binding site. Taken together, these findings direct future studies to get more insights into the role of these variants in the loss of the ABCA1 function.

Graphene-Based Nanoparticles as Potential Treatment Options for Parkinson's Disease: A Molecular Dynamics Study

INTRODUCTION

The study of abnormal aggregation of proteins in different tissues of the body has recently earned great attention from researchers in various fields of science. Concerning neurological diseases, for instance, the accumulation of amyloid fibrils can contribute to Parkinson's disease, a progressively severe neurodegenerative disorder. The most prominent features of this disease are the degeneration of neurons in the substantia nigra and accumulation of α-synuclein aggregates, especially in the brainstem, spinal cord, and cortical areas. Dopamine replacement therapies and other medications have reduced motor impairment and had positive consequences on patients' quality of life. However, if these medications are stopped, symptoms of the disease will recur even more severely. Therefore, the improvement of therapies targeting more basic mechanisms like prevention of amyloid formation seems to be critical. It has been shown that the interactions between monolayers like graphene and amyloids could prevent their fibrillation.

METHODS

For the first time, the impact of four types of last-generation graphene-based nanostructures on the prevention of α-synuclein amyloid fibrillation was investigated in this study by using molecular dynamics simulation tools.

RESULTS

Although all monolayers were shown to prevent amyloid fibrillation, nitrogen-doped graphene (N-Graphene) caused the most instability in the secondary structure of α-synuclein amyloids. Moreover, among the four monolayers, N-Graphene was shown to present the highest absolute value of interaction energy, the lowest contact level of amyloid particles, the highest number of hydrogen bonds between water and amyloid molecules, the highest instability caused in α-synuclein particles, and the most significant decrease in the compactness of α-synuclein protein.

DISCUSSION

Ultimately, it was concluded that N-Graphene could be the most effective monolayer to disrupt amyloid fibrillation, and consequently, prevent the progression of Parkinson's disease.

A comprehensive study on identifying the structural and functional SNPs of human neuronal membrane glycoprotein M6A (GPM6A)

Glycoprotein M6A, a stress related gene, plays an important role in synapse and filopodia formation. Filopodia formation is vital for development, immunity, angiogenesis, wound healing and metastasis. In this study, structural and functional analysis of high-risk SNPs associated with Glycoprotein M6-A were evaluated using six different bioinformatics tools. Results classified T210I, T134I, Y153H, I215T, F156L, T160I, I226T, R247W, R178C, W159R, N157S and P151L as deleterious mutants that are crucial for the structure and function of the protein causing malfunction of M6-a and ultimately leads to disease development. The three-dimensional structure of wild-type M6-a and mutant M6-a were also predicted. Furthermore, the effects of high risk substitutions were also analyzed with interaction with valproic acid. Based on structural models obtained, the binding pocket of ligand bound glycoprotein M6-A structure showed few core interacting residues which are different in the mutant models. Among all substitutions, F156L showed complete loss of binding pocket when interacting with valproic acid as compared to the wild type model. Up to the best of our knowledge this is the first comprehensive study where GPM6A mutations were analyzed. The mechanism of action of GPM6A is still not fully defined which limits the understanding of functional details encoding M6-A. Our results may help enlighten some molecular aspects underlying glycoprotein M6-A. Communicated by Ramaswamy H. Sarma.

Deciphering Molecular Mechanism of the Neuropharmacological Action of Fucosterol through Integrated System Pharmacology and In Silico Analysis

Fucosterol is an algae-derived unique phytosterol having several medicinal properties, including antioxidant, anti-inflammatory, anticholinesterase, neuroprotective, and so on. Accumulated evidence suggests a therapeutic promise of fucosterol in neurodegeneration; however, the in-depth pharmacological mechanism of its neuroprotection is poorly understood. Here, we employed system pharmacology and in silico analysis to elucidate the underlying mechanism of neuropharmacological action of fucosterol against neurodegenerative disorders (NDD). Network pharmacology revealed that fucosterol targets signaling molecules, receptors, enzymes, transporters, transcription factors, cytoskeletal, and various other proteins of cellular pathways, including tumor necrosis factor (TNF), mitogen-activated protein kinase (MAPK), phosphatidylinositol 3-kinase/protein kinase B (PI3K/Akt), neurotrophin, and toll-like receptor (TLR) signaling, which are intimately associated with neuronal survival, immune response, and inflammation. Moreover, the molecular simulation study further verified that fucosterol exhibited a significant binding affinity to some of the vital targets, including liver X-receptor-beta (LXR-β), glucocorticoid receptor (GR), tropomyosin receptor kinase B (TrkB), toll-like receptor 2/4 (TLR2/4), and β-secretase (BACE1), which are the crucial regulators of molecular and cellular processes associated with NDD. Together, the present system pharmacology and in silico findings demonstrate that fucosterol might play a significant role in modulating NDD-pathobiology, supporting its therapeutic application for the prevention and treatment of NDD.