Molecular Insights into the Role of Pathogenic nsSNPs in GRIN2B Gene Provoking Neurodevelopmental Disorders

Structural and dynamics insights into the GBA variants associated with Parkinson's disease

The GBA1 gene encodes for the lysosomal enzyme glucocerebrosidase (GCase), which maintains glycosphingolipid homeostasis and regulates the autophagy process. Genomic variants of GBA1 are associated with Goucher disease; however, several heterozygous variants of GBA (E326K, T369M, N370S, L444P) are frequent high-risk factors for Parkinson's disease (PD). The underlying mechanism of these variants has been revealed through functional and patient-centered research, but the structural and dynamical aspects of these variants have not yet been thoroughly investigated. In the current study, we used a thorough computational method to pinpoint the structural changes that GBA underwent because of genomic variants and drug binding mechanisms. According to our findings, PD-linked nsSNP variants of GBA showed structural variation and abnormal dynamics when compared to wild-typ. The docking analysis demonstrated that the mutants E326K, N370S, and L444P have higher binding affinities for Ambroxol. Root means square deviation (RMSD), Root mean square fluctuation analysis (RMSF), and MM-GBSA analysis confirmed that the Ambroxol are more stable in the binding site of N370S and L444P, and that their binding affinities are stronger as compared to the wild-type and T369M variants of GBA. The evaluation of hydrogen bonds and the calculation of the free binding energy provided additional evidence in favor of this conclusion. When docked with Ambroxol, GBA demonstrated an increase in binding affinity and catalytic activity. Understanding the therapeutic efficacy and potential against the aforementioned changes in the GBA will be beneficial in order to use more efficient methods for developing novel drugs.

The deleterious variants of N-acetylgalactosamine-6-sulfatase (GalN6S) enzyme trigger Morquio a syndrome by disrupting protein foldings

Lysosomal enzymes degrade cellular macromolecules, while their inactivation causes human hereditary metabolic disorders. Mucopolysaccharidosis IVA (MPS IVA; Moquio A syndrome) is one of the lysosomal storage disorders caused by a defective Galactosamine-6-sulfatase (GalN6S) enzyme. In several populations, disease incidence is elevated due to missense mutations brought on by non-synonymous allelic variation in the GalN6S enzyme. Here, we studied the effect of non-synonymous single nucleotide polymorphism (nsSNPs) on the structural dynamics of the GalN6S enzyme and its binding with N-acetylgalactosamine (GalNAc) using all-atom molecular dynamics simulation and an essential dynamics approach. Consequently, in this study, we have identified three functionally disruptive mutations in domain-I and domain-II, that is, S80L, R90W, and S162F, which presumably contribute to post-translational modifications. The study delineated that both domains work cooperatively, and alteration in domain II (S80L, R90W) leads to conformational changes in the catalytic site in domain-I, while mutation S162F mainly provokes higher residual flexibility of domain II. These results show that these mutations impair the hydrophobic core, implying that Morquio A syndrome is caused by misfolding of the GalN6S enzyme. The results also show the instability of the GalN6S-GalNAc complex upon substitution. Overall, the structural dynamics resulting from point mutations give the molecular rationale for Moquio A syndrome and, more importantly, the Mucopolysaccharidoses (MPS) family of diseases, re-establishing MPS IVA as a protein-folding disease.Communicated by Ramaswamy H. Sarma.

Single nucleotide polymorphisms of ANKK1, DDR4, and GRIN2B genes predict behavior in a prospective cohort of Mexican children and adolescents

Numerous studies have established associations between single nucleotide polymorphisms (SNPs) and various behavioral and neurodevelopmental conditions. This study explores the links between SNPs in candidate genes involved in central nervous system (CNS) physiology and their implications for the behavioral and emotional aspects in children and teenagers. A total of 590 participants, aged 7-15 years, from the Early Life Exposures In Mexico To Environmental Toxicants (ELEMENT) cohort study in Mexico City, underwent genotyping for at least one of 15 CNS gene-related SNPs at different timepoints. We employed multiple linear regression models to assess the potential impact of genetic variations on behavioral and cognitive traits, as measured by the Behavioral Assessment System for Children (BASC) and Conners parent rating scales. Significant associations were observed, including the rs1800497 TC genotype (ANKK1) with the Cognitive Problems/Inattention variable (p value = 0.003), the rs1800955 CT genotype (DDR4) with the Emotional Lability Global index variable (p value = 0.01), and the rs10492138 GA and rs7970177 TC genotypes (GRIN2B) with the Depression variable (p values 0.007 and 0.012, respectively). These finds suggest potential genetic profiles associated with "risk" and "protective" behaviors for these SNPs. Our results provide valuable insights into the role of genetic variations in neurobehavior and highlight the need for further research in the early identification and intervention in individuals at risk for these conditions.

Identification of potential natural products derived from fungus growing termite, inhibiting Pseudomonas aeruginosa quorum sensing protein LasR using molecular docking and molecular dynamics simulation approach

Pseudomonas aeruginosa, the most common opportunistic pathogen, is becoming antibiotic-resistant worldwide. The fate of P. aeruginosa, a multidrug-resistant strain, can be determined by multidrug efflux pumps, enzyme synthesis, outer membrane protein depletion, and target alterations. Microbial niches have long used quorum sensing (QS) to synchronize virulence gene expression. Computational methods can aid in the development of novel P. aeruginosa drug-resistant treatments. The tripartite symbiosis in termites that grow fungus may help special microbes find new antimicrobial drugs. To find anti-quorum sensing natural products that could be used as alternative therapies, a library of 376 fungal-growing termite-associated natural products (NPs) was screened for their physicochemical properties, pharmacokinetics, and drug-likeness. Using GOLD, the top 74 NPs were docked to the QS transcriptional regulator LasR protein. The five lead NPs with the highest gold score and drug-like properties were chosen for a 200-ns molecular dynamics simulation to test the competitive activity of different compounds against negative catechin. Fridamycin and Daidzein had stable conformations, with mean RMSDs of 2.48 and 3.67 Å, respectively, which were similar to Catechin's 3.22 Å. Fridamycin and Daidzein had absolute binding energies of -71.186 and -52.013 kcal/mol, respectively, which were higher than the control's -42.75 kcal/mol. All the compounds within the active site of the LasR protein were kept intact by Trp54, Arg55, Asp67, and Ser123. These findings indicate that termite gut and fungus-associated NPs, specifically Fridamycin and Daidzein, are potent QS antagonists that can be used to treat P. aeruginosa's multidrug resistance.Communicated by Ramaswamy H. Sarma.

Identification of novel peptide inhibitors for the KRas-G12C variant to prevent oncogenic signaling

Kirsten rat sarcoma viral oncogene homolog (KRas) activating mutations are common in solid tumors, accounting for 90%, 45%, and 35% of pancreatic, colorectal, and lung cancers (LC), respectively. Each year, nearly 150k new cases (both men and women) of KRas-mutated malignancies are reported in the United States. NSCLC (non-small cell lung cancer) accounts for 80% of all LC cases. KRas mutations are found in 15% to 25% of NSCLC patients. The main cause of NSCLC is the KRas-G12C mutation. The drugs Sotorasib and Adagrasib were recently developed to treat advanced NSCLC caused by the KRas-G12C mutation. Most patients do not respond to KRas-G12C inhibitors due to cellular, molecular, and genetic resistance. Because of their safety, efficacy, and selectivity, peptide inhibitors have the potential to treat newly developing KRas mutations. Based on the KRas mutations, peptide inhibitors that are highly selective and specific to individual lung cancers can be rationally designed. The current study uses an alanine and residue scanning approach to design peptide inhibitors for KRas-G12C based on the known peptide. Our findings show that substitution of F3K, G11T, L8C, T14C, K13D, G11S, and G11P considerably enhances the binding affinity of the novel peptides, whereas F3K, G11T, L8C, and T14C peptides have higher stability and favorable binding to the altered peptides. Overall, our study paves the road for the development of potential therapeutic peptidomimetics that target the KRas-G12C complex and may inhibit the KRas and SOS complex from interacting.Communicated by Ramaswamy H. Sarma.

Autosomal Dominant Intellectual Development Disorder-6 (MRD6) Without Seizures Linked to a De Novo Mutation in the grin2b Gene Revealed by Exome Sequencing: A Case Report of a Moroccan Child

Autosomal dominant intellectual development disorder-6 (MRD6) arises from a grin2b gene mutation, inducing neurodevelopmental issues. The effects of MRD6 encompass cognitive disabilities, seizures, muscle tone decline, and autism-like traits. Its severity ranges from mild impairment to severe epilepsy. The disorder's rarity is emphasized by roughly 100 reported GRIN2B-related cases, spotlighting the gene's significance in brain development. We present the case of a three-year-old Moroccan boy who was referred to a neuropediatric department for a molecular diagnosis. Initial genetic testing yielded inconclusive results, and subsequent tests for Angelman syndrome and metabolic diseases showed no abnormalities. Given the complexity of the disorder, exome sequencing was employed to identify the underlying genetic cause. Exome sequencing identified a nonsense (STOP) mutation c.3912C>G (p.Tyr1304Ter) in the grin2b gene in the heterozygous state known to be present in MRD6 (Online Mendelian Inheritance in Man (OMIM) 613970). The family segregation study shows that this is a de novo variant, which is confirmed by Sanger sequencing. This variant has not been previously reported in the GnomAD database. Based on current scientific knowledge, the variant is considered pathogenic (PVS1, PS2, PM2, PP3, PP5) according to the criteria of the American College of Medical Genetics and Genomics (ACMG). The mutation in the grin2b gene (p.Tyr1304Ter) was predicted to be deleterious through bioinformatics analysis tools. This study highlights the crucial role of the grin2b gene in normal brain development and communication within the nervous system. It also sheds light on the impact of a novel genetic mutation, identified through exome sequencing, on causing an intellectual developmental disorder in a child patient from Morocco.

A comprehensive protein interaction map and druggability investigation prioritized dengue virus NS1 protein as promising therapeutic candidate

Dengue Virus (DENV) is a serious threat to human life worldwide and is one of the most dangerous vector-borne diseases, causing thousands of deaths annually. We constructed a comprehensive PPI map of DENV with its host Homo sapiens and performed various bioinformatics analyses. We found 1195 interactions between 858 human and 10 DENV proteins. Pathway enrichment analysis was performed on the two sets of gene products, and the top 5 human proteins with the maximum number of interactions with dengue viral proteins revealed noticeable results. The non-structural protein NS1 in DENV had the maximum number of interactions with the host protein, followed by NS5 and NS3. Among the human proteins, HBA1 and UBE2I were associated with 7 viral proteins, and 3 human proteins (CSNK2A1, RRP12, and HSP90AB1) were found to interact with 6 viral proteins. Pharmacophore-based virtual screening of millions of compounds in the public databases was performed to identify potential DENV-NS1 inhibitors. The lead compounds were selected based on RMSD values, docking scores, and strong binding affinities. The top ten hit compounds were subjected to ADME profiling which identified compounds C2 (MolPort-044-180-163) and C6 (MolPort-001-742-737) as lead inhibitors against DENV-NS1. Molecular dynamics trajectory analysis and intermolecular interactions between NS1 and the ligands displayed the molecular stability of the complexes in the cellular environment. The in-silico approaches used in this study could pave the way for the development of potential specie-specific drugs and help in eliminating deadly viral infections. Therefore, experimental and clinical assays are required to validate the results of this study.

In Silico Characterization and Analysis of Clinically Significant Variants of Lipase-H (LIPH Gene) Protein Associated with Hypotrichosis

Hypotrichosis is an uncommon type of alopecia (hair loss) characterized by coarse scalp hair caused by the reduced or fully terminated activity of the Lipase-H (LIPH) enzyme. LIPH gene mutations contribute to the development of irregular or non-functional proteins. Because several cellular processes, including cell maturation and proliferation, are inhibited when this enzyme is inactive, the hair follicles become structurally unreliable, undeveloped, and immature. This results in brittle hair, as well as altered hair shaft development and structure. Because of these nsSNPs, the protein's structure and/or function may be altered. Given the difficulty in discovering functional SNPs in genes associated with disease, it is possible to assess potential functional SNPs before conducting broader population investigations. As a result, in our in silico analysis, we separated potentially hazardous nsSNPs of the LIPH gene from benign representatives using a variety of sequencing and architecture-based bioinformatics approaches. Using seven prediction algorithms, 9 out of a total of 215 nsSNPs were shown to be the most likely to cause harm. In order to distinguish between potentially harmful and benign nsSNPs of the LIPH gene, in our in silico investigation, we employed a range of sequence- and architecture-based bioinformatics techniques. Three nsSNPs (W108R, C246S, and H248N) were chosen as potentially harmful. The present findings will likely be helpful in future large population-based studies, as well as in drug discovery, particularly in the creation of personalized medicine, since this study provides an initial thorough investigation of the functional nsSNPs of LIPH.

In Silico Analysis of nsSNPs of Human KRAS Gene and Protein Modeling Using Bioinformatic Tools

The KRAS gene belongs to the RAS family and codes for 188 amino acid residues of KRAS protein, with a molecular mass of 21.6 kD. Non-synonymous single-nucleotide polymorphisms (nsSNPs) have been identified within the coding region in which some are associated with different diseases. However, structural changes are not well defined yet. In this study, we first categorized SNPs in the KRAS coding area and then used computational methods to determine their impact on the protein structure and stability. In addition, the three-dimensional model of KRAS was taken from the Protein Data Bank for structural modeling. Furthermore, genomic data were extracted from a variety of sources, including the 1000 Genome Project, dbSNPs, and ENSEMBLE, and assessed through in silico methods. Based on various tools used in this study, 10 out of 48 missense SNPs with rsIDs were found deleterious. The substitution of alanine for proline at position 146 pushed several residues toward the center of the protein. Arginine instead of leucine has a minor effect on protein structure and stability. In addition, the substitution of proline for leucine at the 34th position disrupted the structure and led to a bigger size than the wild-type protein, hence interrupting the protein interaction. Using the well-intended computational approach and applying several bioinformatic tools, we characterized and identified most damaging nsSNPs and further explored the structural dynamics and stability of KRAS protein.

Comprehensive in-silico analysis of deleterious SNPs in APOC2 and APOA5 and their differential expression in cancer and cardiovascular diseases conditions

Genetic variations in APOC2 and APOA5 genes involve activating lipoprotein lipase (LPL), responsible for the hydrolysis of triglycerides (TG) in blood and whose impaired functions affect the TG metabolism and are associated with metabolic diseases. In this study, we investigate the biological significance of genetic variations at the DNA sequence and structural level using various computational tools. Subsequently, 8 (APOC2) and 17 (APOA5) non-synonymous SNPs (nsSNPs) were identified as high-confidence deleterious SNPs based on the effects of the mutations on protein conservation, stability, and solvent accessibility. Furthermore, based on our docking results, the interaction of native and mutant forms of the corresponding proteins with LPL depicts differences in root mean square deviation (RMSD), and binding affinities suggest that these mutations may affect their function. Furthermore, in vivo, and in vitro studies have shown that differential expression of these genes in disease conditions due to the influence of nsSNPs abundance may be associated with promoting the development of cancer and cardiovascular diseases. Preliminary screening using computational methods can be a helpful start in understanding the effects of mutations in APOC2 and APOA5 on lipid metabolism; however, further wet-lab experiments would further strengthen the conclusions drawn from the computational study.

Fragment-Based Approaches Identified Tecovirimat-Competitive Novel Drug Candidate for Targeting the F13 Protein of the Monkeypox Virus

Monkeypox is a serious public health issue in tropical and subtropical areas. Antivirals that target monkeypox proteins might lead to more effective and efficient therapy. The F13 protein is essential for the growth and maturation of the monkeypox virus. F13 inhibition might be a viable therapeutic target for monkeypox. The in silico fragment-based drug discovery method for developing antivirals may provide novel therapeutic options. In this study, we generated 800 compounds based on tecovirimat, an FDA-approved drug that is efficacious at nanomolar quantities against monkeypox. These compounds were evaluated to identify the most promising fragments based on binding affinity and pharmacological characteristics. The top hits from the chemical screening were docked into the active site of the F13 protein. Molecular dynamics simulations were performed on the top two probable new candidates from molecular docking. The ligand-enzyme interaction analysis revealed that the C2 ligand had lower binding free energy than the standard ligand tecovirimat. Water bridges, among other interactions, were shown to stabilize the C2 molecule. Conformational transitions and secondary structure changes in F13 protein upon C2 binding show more native three-dimensional folding of the protein. Prediction of pharmacological properties revealed that compound C2 may be promising as a drug candidate for monkeypox fever. However, additional in vitro and in vivo testing is required for validation.

Identification of fungus-growing termite-associated halogenated-PKS maduralactomycin a as a potential inhibitor of MurF protein of multidrug-resistant Acinetobacter baumannii

Multidrug-resistant Acinetobacter baumannii infections have become a major public health concern globally. Inhibition of its essential MurF protein has been proposed as a potential target for broad-spectrum drugs. This study aimed to evaluate the potential of a novel ecological niche of 374 fungus-growing termite associated Natural Products (NPs). The molecular docking and computational pharmacokinetics screened four compounds, i.e., Termstrin B, Fridamycin A, Maduralactomycin A, and Natalenamide C, as potential compounds that have higher binding affinities and favourable protein-ligand interactions. The compound Maduralactomycin A induced more stability based on its lowest average RMSD value (2.31 Å) and low standard deviation (0.35) supported by the consistent flexibility and β-factor during the protein's time-dependent motion. While hydrogen bond analysis indicated that Termstrin B has formed the strongest intra-protein interaction, solvent accessibility was in good agreement with Maduralactomycin A compactness. Maduralactomycin A has the strongest binding energy among all the compounds (-348.48 kcal/mol) followed by Termstrin B (-321.19 kcal/mol). Since these findings suggest Maduralactomycin A and Termstrin B as promising candidates for inhibition of MurF protein, the favourable binding energies of Maduralactomycin A make it a more important compound to warrant further investigation. However, experimental validation using animal models and clinical trials is recommended before reaching any final conclusions.

Structural Evaluation and Conformational Dynamics of ZNF141T474I Mutation Provoking Postaxial Polydactyly Type A

Postaxial Polydactyly (PAP) is a congenital disorder of limb abnormalities characterized by posterior extra digits. Mutations in the N-terminal region of the Zinc finger protein 141 (ZNF141) gene were recently linked with PAP type A. Zinc finger proteins exhibit similarity at their N-terminal regions due to C2-H2 type Zinc finger domains, but their functional preferences vary significantly by the binding patterns of DNA. Methods: This study delineates the pathogenic association, miss-fold aggregation, and conformational paradigm of a missense variant (c.1420C > T; p.T474I) in ZNF141 gene segregating PAP through a molecular dynamics simulations approach. Results: In ZNF141 protein, helices play a crucial role by attaching three specific target DNA base pairs. In ZNF141T474I protein, H1, H3, and H6 helices attain more flexibility by acquiring loop conformation. The outward disposition of the proximal portion of H9-helix in mutant protein occurs due to the loss of prior beta-hairpins at the C terminal region of the C2-H2 domain. The loss of hydrogen bonds and exposure of hydrophobic residues to solvent and helices turning to loops cause dysfunction of ZNF141 protein. These significant changes in the stability and conformation of the mutant protein were validated using essential dynamics and cross-correlation maps, which revealed that upon point mutation, the overall motion of the proteins and the correlation between them were completely different, resulting in Postaxial polydactyly type A. Conclusions: This study provides molecular insights into the structural association of ZNF141 protein with PAP type A. Identification of active site residues and legends offers new therapeutic targets for ZNF141 protein. Further, it reiterates the functional importance of the last residue of a protein.