Identification and structural characterization of deleterious non-synonymous single nucleotide polymorphisms in the human SKP2 gene

Field measures of strengthen plant-microbial remediation of PAHs-FQs compound pollution

In this study, five PAHs (benzo [b] fluoranthene (BbF), phenanthrene (Phe), fluoranthene (Flu), fluorene (Fl), benzo [A] pyrene (Bap)), and five FQs (ofloxacin (OFL), enrofloxacin (ENR), ciprofloxacin (CIP), norfloxacin (NOR), lomefloxacin (LOM)) were selected as ligands; peroxidase (1NML) was selected as receptor degrading protein. In the plant-microbial degradation, the factors with significant inhibitory effects are NOR, Bap, CIP, ENR, OFL, Flu, LOM, Phe, Fl, and BbF by the fractional factorial design experiment and molecular docking-assisted molecular dynamics methods. Using Taguchi experiment and molecular dynamics simulation methods, the main external field measures were designed and screened to effectively promote the degradation of PAHs-FQs under the combined pollution scenarios of Bap-CIP and BbF-NOR, respectively. The peroxidase mutation design plans with enhanced substrate affinity were then designed and screened using the DS software by predicting the virtual key amino acid of peroxidase. The novel biodegradable enzymes 2YCD-1, 2YCD-4, 2YCD-5, 2YCD-7, and 2YCD-9 had better structures and showed excellent degradability for PAHs and FQs. This study explored the degradation rules of the composite pollutants in the coexistence systems of multiple PAHs and FQs, providing the best external field measures for the control and treatment of the combined pollution effects of different PAHs and FQs. Overall, the current study has important practical significance for promoting the plant-microbial joint remediation of PAHs-FQs pollution and for reducing the combined pollution of PAHs and FQs in farmland systems.

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.

A novel causative functional mutation in GATA6 gene is responsible for familial dilated cardiomyopathy as supported by in silico functional analysis

Dilated cardiomyopathy (DCM), one of the most common types of cardiomyopathies has a heterogeneous nature and can be seen in Mendelian forms. Next Generation Sequencing is a powerful tool for identifying novel variants in monogenic disorders. We used whole-exome sequencing (WES) and Sanger sequencing techniques to identify the causative mutation of DCM in an Iranian pedigree. We found a novel variant in the GATA6 gene, leading to substituting Histidine by Tyrosine at position 329, observed in all affected family members in the pedigree, whereas it was not established in any of the unaffected ones. We hypothesized that the H329Y mutation may be causative for the familial pattern of DCM in this family. The predicted models of GATA6 and H329Y showed the high quality according to PROCHECK and ERRAT. Nonetheless, simulation results revealed that the protein stability decreased after mutation, while the flexibility may have been increased. Hence, the mutation led to the increased compactness of GATA6. Overall, these data indicated that the mutation could affect the protein structure, which may be related to the functional impairment of GATA6 upon H329Y mutation, likewise their involvement in pathologies. Further functional investigations would help elucidating the exact mechanism.

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.

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.

Whole proteome screening and identification of potential epitopes of SARS-CoV-2 for vaccine design-an immunoinformatic, molecular docking and molecular dynamics simulation accelerated robust strategy

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the most cryptic pandemic outbreak of the 21st century, has gripped more than 1.8 million people to death and infected almost eighty six million. As it is a new variant of SARS, there is no approved drug or vaccine available against this virus. This study aims to predict some promising cytotoxic T lymphocyte epitopes in the SARS-CoV-2 proteome utilizing immunoinformatic approaches. Firstly, we identified 21 epitopes from 7 different proteins of SARS-CoV-2 inducing immune response and checked for allergenicity and conservancy. Based on these factors, we selected the top three epitopes, namely KAYNVTQAF, ATSRTLSYY, and LTALRLCAY showing functional interactions with the maximum number of MHC alleles and no allergenicity. Secondly, the 3D model of selected epitopes and HLA-A*29:02 were built and Molecular Docking simulation was performed. Most interestingly, the best two epitopes predicted by docking are part of two different structural proteins of SARS-CoV-2, namely Membrane Glycoprotein (ATSRTLSYY) and Nucleocapsid Phosphoprotein (KAYNVTQAF), which are generally target of choice for vaccine designing. Upon Molecular Docking, interactions between selected epitopes and HLA-A*29:02 were further validated by 50 ns Molecular Dynamics (MD) simulation. Analysis of RMSD, Rg, SASA, number of hydrogen bonds, RMSF, MM-PBSA, PCA, and DCCM from MD suggested that ATSRTLSYY is the most stable and promising epitope than KAYNVTQAF epitope. Moreover, we also identified B-cell epitopes for each of the antigenic proteins of SARS CoV-2. Findings of our work will be a good resource for wet lab experiments and will lessen the timeline for vaccine construction.Communicated by Ramaswamy H. Sarma.

Unravelling the molecular basis of AM-6494 high potency at BACE1 in Alzheimer's disease: an integrated dynamic interaction investigation

β-amyloid precursor protein cleaving enzyme1 (BACE1) has prominently been an important drug design target implicated in Alzheimer's disease pathway. The failure rate of most of the already tested drugs at different clinical phases remains a major concern. Recently, AM-6494 was reported as a novel potent, highly selective, and orally effective inhibitor against BACE1. AM-6494 displayed no alteration of skin/fur colour in animal studies, an adverse effect common to previous BACE1 inhibitors. However, the atomistic molecular mechanism of BACE1 inhibition by AM-6494 remains unclear. To elucidate the binding mechanism of AM-6494 relative to umibecestat (CNP-520) as well as the structural changes when bound to BACE1, advanced computational techniques such as accelerated MD simulation and principal component analysis have been utilised. The results demonstrated higher binding affinity of AM-6494 at BACE1 with van der Waals as dominant energy contributor compared to umibecestat. Conformational monitoring of the β-hairpin flap covering the active site revealed an effective flap closure when bound with AM-6494 compared to CNP-520, which predominantly alternates between semi-open and closed conformations. The observed effective flap closure of AM-6494 explains its higher inhibitory power towards BACE1. Besides the catalytic Asp32/228 dyad, Tyr14, Leu30, Tyr71 and Gly230 represent critical residues in the potency of these inhibitors at BACE1 binding interface. The findings highlighted in this research provide a basis to explain AM-6494 high inhibitory potency and might assist in the design of new inhibitors with improved selectivity and potency for BACE1.

Exploring the structural basis and atomistic binding mechanistic of the selective antagonist blockade at D3 dopamine receptor over D2 dopamine receptor

More recently, there has been a paradigm shift toward selective drug targeting in the treatment of neurological disorders, including drug addiction, schizophrenia, and Parkinson's disease mediated by the different dopamine receptor subtypes. Antagonists with higher selectivity for D₃ dopamine receptor (D3DR) over D₂ dopamine receptor (D2DR) have been shown to attenuate drug-seeking behavior and associated side effects compared to non-subtype selective antagonists. However, high conservations among constituent residues of both proteins, particularly at the ligand-binding pockets, remain a challenge to therapeutic drug design. Recent studies have reported the discovery of two small-molecules R-VK4-40 and Y-QA31 which substantially inhibited D3DR with >180-fold selectivity over D2DR. Therefore, in this study, we seek to provide molecular and structural insights into these differential binding mechanistic using meta-analytic computational simulation methods. Findings revealed that R-VK4-40 and Y-QA31 adopted shallow binding modes and were more surface-exposed at D3DR while on the contrary, they exhibited deep hydrophobic pocket binding at D2DR. Also, two non-conserved residues; Tyr36^1.39 and Ser182^45.51 were identified in D3DR, based on their crucial roles and contributions to the selective binding of R-VK4-40 and Y-QA31. Importantly, both antagonists exhibited high affinities in complex with D3DR compared to D2DR, while van der Waals energies contributed majorly to their binding and stability. Structural analyses also revealed the distinct stabilizing effects of both compounds on D3DR secondary architecture relative to D2DR. Therefore, findings herein pinpointed the origin and mechanistic of selectivity of the compounds, which may assist in the rational design of potential small molecules of the D₂ -like dopamine family receptor subtype with improved potency and selectivity.

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.

Mechanistic insights into the deleterious roles of Nasu-Hakola disease associated TREM2 variants

Recently, the critical roles played by genetic variants of TREM2 (Triggering Receptor Expressed on Myeloid cells 2) in Alzheimer's disease have been aggressively highlighted. However, few studies have focused on the deleterious roles of Nasu-Hakola disease (NHD) associated TREM2 variants. In order to get insights into the contributions made by these variants to neurodegeneration, we investigated the influences of four NHD associated TREM2 mutations (Y38C, W50C, T66M, and V126G) on loss-of-function, and followed this with in silico prediction and conventional molecular dynamics simulation. NHD mutations were predicted to be highly deleterious by eight different in silico bioinformatics tools and found to induce conformational changes by molecular dynamics simulation. As compared with the wild-type, the four variants produced substantial differences in the collective motions of loop regions, which not only promoted structural remodeling in the CDR2 (complementarity-determining region 2) loop but also in the CDR1 loop, by changing inter- and intra-loop hydrogen bonding networks. In addition, structural studies in a free energy landscape analysis showed that Y38, T66, and V126 are crucial for maintaining the structural features of CDR1 and CDR2 loops, and that mutations in these positions produced steric clashes and loss of ligand binding. These results showed the presence of mutations in the TREM2 ectodomain induced flexibility and caused structural alterations. Dynamical scenarios, as provided by the present study, may be critical to our understanding of the roles of these TREM2 mutations in neurodegenerative diseases.

Structural and Dynamic Characterizations Highlight the Deleterious Role of SULT1A1 R213H Polymorphism in Substrate Binding

Sulfotransferase 1A1 (SULT1A1) is responsible for catalyzing various types of endogenous and exogenous compounds. Accumulating data indicates that the polymorphism rs9282861 (R213H) is responsible for inefficient enzymatic activity and associated with cancer progression. To characterize the detailed functional consequences of this mutation behind the loss-of-function of SULT1A1, the present study deployed molecular dynamics simulation to get insights into changes in the conformation and binding energy. The dynamics scenario of SULT1A1 in both wild and mutated types as well as with and without ligand showed that R213H induced local conformational changes, especially in the substrate-binding loop rather than impairing overall stability of the protein structure. The higher conformational changes were observed in the loop3 (residues, 235-263), turning loop conformation to A-helix and B-bridge, which ultimately disrupted the plasticity of the active site. This alteration reduced the binding site volume and hydrophobicity to decrease the binding affinity of the enzyme to substrates, which was highlighted by the MM-PBSA binding energy analysis. These findings highlight the key insights of structural consequences caused by R213H mutation, which would enrich the understanding regarding the role of SULT1A1 mutation in cancer development and also xenobiotics management to individuals in the different treatment stages.

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.