Impact of I30T and I30M substitution in MPZ gene associated with Dejerine–Sottas syndrome type B (DSSB): A molecular modeling and dynamics

Identification and structural characterization of a pathogenic ARSA missense variant in two consanguineous families from Jammu and Kashmir (India) with late infantile metachromatic leukodystrophy

BACKGROUND

Metachromatic leukodystrophy (MLD) is a rare lysosomal storage disorder caused by a deficiency of Arylsulfatase A (ARSA) enzyme activity. Its clinical manifestations include progressive motor and cognitive decline. ARSA gene mutations are frequent in MLD.

METHODS AND RESULTS

In the present study, whole exome sequencing (WES) was employed to decipher the genetic cause of motor and cognitive decline in proband's of two consanguineous families from J&K (India). Clinical investigations using radiological and biochemical analysis revealed MLD-like features. WES confirmed a pathogenic variant in the ARSA gene. Molecular simulation dynamics was applied for structural characterization of the variant.

CONCLUSION

We report the identification of a pathogenic missense variant (c.1174 C > T; p.Arg390Trp) in the ARSA gene in two cases of late infantile MLD from consanguineous families in Jammu and Kashmir, India. Our study utilized genetic analysis and molecular dynamics simulations to identify and investigate the structural consequences of this mutation. The molecular dynamics simulations revealed significant alterations in the structural dynamics, residue interactions, and stability of the ARSA protein harbouring the p.Arg390Trp mutation. These findings provide valuable insights into the molecular mechanisms underlying the pathogenicity of this variant in MLD.

Transcriptome Profiling of miRNA-mRNA Interactions and Associated Mechanisms in Chemotherapy-Induced Neuropathic Pain

Chemotherapy-induced neuropathic pain (CINP) is a dose-limiting adverse event affecting 40% of chemotherapy patients. MiRNA-mRNA interaction plays an important role in various processes. However, detailed profiling of miRNA-mRNA interactions in CINP remains unclear. Here, a rat-based CINP model was established using paclitaxel, followed by nociceptive behavioral tests related to mechanical allodynia, thermal hyperalgesia, and cold allodynia. The landscape of miRNA-mRNA interaction in the spinal dorsal horn was investigated through mRNA transcriptomics and small RNA sequencing. Under CINP condition, 86 differentially expressed mRNAs and 56 miRNAs were identified. Gene Set Enrichment Analysis (GSEA), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses indicated the activity of Odorant binding, postsynaptic specialization and synaptic density, extracellular matrix, mitochondrial matrix, retrograde endocannabinoid signaling, and GTPase activity. Protein-protein interaction (PPI), networks of circRNA-miRNA-mRNA, lncRNA-miRNA-mRNA, and TF-genes were demonstrated. We next explored the immune infiltration microenvironment and found a higher infiltration abundance of Th17 and a lower abundance of MDSC in CINP. RT-qPCR and dual-luciferase assays were used to verify the sequencing results, and single-cell analysis based on the SekSeeq database was conducted. Combined with bioinformatics analyses and experimental validations, Mpz, a protein-coding gene specifically expressed in Schwann cells, was found critical in maintaining CINP under miRNA regulation. Therefore, these data highlight the expression patterns of miRNA-mRNA, and the underlying mechanism in the spinal dorsal horn under CINP condition, and Mpz may serve as a promising therapeutic target for patients with CINP.

Clinical and mutational spectrum of paediatric Charcot-Marie-Tooth disease in a large cohort of Chinese patients

Background: Charcot-Marie-Tooth disease (CMT) is the most common inherited neurological disorder suffered in childhood. To date, the disease features have not been extensively characterized in the Chinese paediatric population. In this study, we aimed to analyse the clinical profiles and genetic distributions of a paediatric CMT cohort in China. Methods: A total of 181 paediatric CMT patients were enrolled. After preexcluding PMP22 duplication/deletion by multiplex ligation-dependent probe amplification (MLPA), Sanger sequencing, targeted next-generation sequencing (NGS) or whole-exome sequencing (WES) was performed to obtain a genetic diagnosis. Detailed information was collected to explore the spectrum of subtypes and genotype-phenotype correlations. Results: Pathogenic mutations were identified in 68% of patients in this study; with PMP22 duplication, MFN2 and GJB1 were the most frequent disease-causing genes. Of note, respect to the higher prevalence worldwide, CMT1A (18.2%) was relatively lower in our cohort. Besides, the mean age at onset (8.3 ± 5.7 years) was significantly older in our series. In genotype-phenotype analyse, PMP22 point mutations were considered the most severe genotypes and were mostly de novo. In addition, the de novo mutations were identified in up to 12.7% of all patients, which was higher than that in other studies. Conclusion: We identified a relatively lower detection rate of PMP22 duplication and a higher frequency of de novo variants among paediatric patients in China. We also identified the genetic and phenotypic heterogeneity of this cohort, which may provide clues for clinicians in directing genetic testing strategies for Chinese patients with early-onset CMT.

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.

In silico and in vitro effects of the I30T mutation on myelin protein zero instability in the cell membrane

Charcot-Marie-Tooth (CMT) diseases are a heterogeneous group of genetic peripheral neuropathies caused by mutations in a variety of genes, which are involved in the development and maintenance of peripheral nerves. Myelin protein zero (MPZ) is expressed by Schwann cells, and MPZ mutations can lead to primarily demyelinating polyneuropathies including CMT type 1B. Different mutations demonstrate various forms of disease pathomechanisms, which may be beneficial in understanding the disease cellular pathology. Our molecular dynamics simulation study on the possible impacts of I30T mutation on the MPZ protein structure suggested a higher hydrophobicity and thus lower stability in the membranous structures. A study was also conducted to predict native/mutant MPZ interactions. To validate the results of the simulation study, the native and mutant forms of the MPZ protein were separately expressed in a cellular model, and the protein trafficking was chased down in a time course pattern. In vitro studies provided more evidence on the instability of the MPZ protein due to the mutation. In this study, qualitative and quantitative approaches were adopted to confirm the instability of mutant MPZ in cellular membranes.

The impact of missense mutation in PIGA associated to paroxysmal nocturnal hemoglobinuria and multiple congenital anomalies-hypotonia-seizures syndrome 2: A computational study

Paroxysmal nocturnal hemoglobinuria (PNH) is an acquired clonal blood disorder that manifests with hemolytic anemia, thrombosis, and peripheral blood cytopenias. The disease is caused by the deficiency of two glycosylphosphatidylinositols (GPI)-anchored proteins (CD55 and CD59) in the hemopoietic stem cells. The deficiency of GPI-anchored proteins has been associated with the somatic mutations in phosphatidylinositol glycan class A (PIGA). However, the mutations that do not cause PNH is associated with the multiple congenital anomalies-hypotonia-seizures syndrome 2 (MCAHS2). To best of our knowledge, no computational study has been performed to explore at an atomistic level the impact of PIGA missense mutations on the structure and dynamics of the protein. Therefore, we focused our study to provide molecular insights into the changes in protein structural dynamics upon mutation. In the initial step, screening for the most pathogenic mutations from the pool of publicly available mutations was performed. Further, to get a better understanding, pathogenic mutations were mapped to the modeled structure and the resulting protein was subjected to 100 ns molecular dynamics simulation. The residues close to C- and N-terminal regions of the protein were found to exhibit greater flexibility upon mutation. Our study suggests that four mutations are highly effective in altering the structural conformation and stability of the PIGA protein. Among them, mutant G48D was found to alter protein's structural dynamics to the greatest extent, both on a local and a global scale.

Molecular insights of the G2019S substitution in LRRK2 kinase domain associated with Parkinson's disease: A molecular dynamics simulation approach

The G2019S substitution in the Leucine-rich repeat kinase 2 (LRRK2) is significantly associated with Parkinson's disease (PD). This substitution was identified in both familial and sporadic forms of PD with a higher frequency. Few computational studies have reported the impact of G2019S substitution on inhibitors of the kinase domain of LRRK2. However, no computational study deeply investigated the possible impact of the G2019S substitution on the kinase domain in its Apo conformation. Therefore, in this study, we used 200 ns molecular dynamic simulation using the GROMACS 5.1.4 package software to investigate the impact of the G2019S substitution on the structure of the kinase domain of LRRK2. Our results indicate that the G2019S substitution affects the dynamics and stability of LRRK2 by decreasing the flexibility and increasing the compactness of the kinase domain and showing its tendency to be in an active conformation for long time interval because of the high energy barrier between active and inactive conformation. This study predicts the molecular pathogenicity mechanism of the G2019S on patients with PD and provides a potential platform for developing therapeutics for patients with PD that harbor this amino acid substitution.

Understanding the structure-function relationship of HPRT1 missense mutations in association with Lesch-Nyhan disease and HPRT1-related gout by in silico mutational analysis

The nucleotide salvage pathway is used to recycle degraded nucleotides (purines and pyrimidines); one of the enzymes that helps to recycle purines is hypoxanthine guanine phosphoribosyl transferase 1 (HGPRT1). Therefore, defects in this enzyme lead to the accumulation of DNA and nucleotide lesions and hence replication errors and genetic disorders. Missense mutations in hypoxanthine phosphoribosyl transferase 1 (HPRT1) are associated with deficiencies such as Lesch-Nyhan disease and chronic gout, which have manifestations such as arthritis, neurodegeneration, and cognitive disorders. In the present study, we collected 88 non-synonymous single nucleotide polymorphisms (nsSNPs) from the UniProt, dbSNP, ExAC, and ClinVar databases. We used a series of sequence-based and structure-based in silico tools to prioritize and characterize the most pathogenic and stabilizing or destabilizing nsSNPs. Moreover, to obtain the structural impact of the pathogenic mutations, we mapped the mutations to the crystal structure of the HPRT protein. We further subjected these mutant proteins to a 50 ns molecular dynamics simulation (MDS). The MDS trajectory showed that all mutant proteins altered the structural conformation and dynamic behavior of the HPRT protein and corroborated its association with LND and gout. This study provides essential information regarding the use of HPRT protein mutants as potential targets for therapeutic development.

Computational insights of K1444N substitution in GAP-related domain of NF1 gene associated with neurofibromatosis type 1 disease: a molecular modeling and dynamics approach

The NF1 gene encodes for neurofibromin protein, which is ubiquitously expressed, but most highly in the central nervous system. Non-synonymous SNPs (nsSNPs) in the NF1 gene were found to be associated with Neurofibromatosis Type 1 disease, which is characterized by the growth of tumors along nerves in the skin, brain, and other parts of the body. In this study, we used several in silico predictions tools to analyze 16 nsSNPs in the RAS-GAP domain of neurofibromin, the K1444N (K1423N) mutation was predicted as the most pathogenic. The comparative molecular dynamic simulation (MDS; 50 ns) between the wild type and the K1444N (K1423N) mutant suggested a significant change in the electrostatic potential. In addition, the RMSD, RMSF, Rg, hydrogen bonds, and PCA analysis confirmed the loss of flexibility and increase in compactness of the mutant protein. Further, SASA analysis revealed exchange between hydrophobic and hydrophilic residues from the core of the RAS-GAP domain to the surface of the mutant domain, consistent with the secondary structure analysis that showed significant alteration in the mutant protein conformation. Our data concludes that the K1444N (K1423N) mutant lead to increasing the rigidity and compactness of the protein. This study provides evidence of the benefits of the computational tools in predicting the pathogenicity of genetic mutations and suggests the application of MDS and different in silico prediction tools for variant assessment and classification in genetic clinics.

A profound computational study to prioritize the disease-causing mutations in PRPS1 gene

Charcot-Marie-Tooth disease (CMT) is one of the most commonly inherited congenital neurological disorders, affecting approximately 1 in 2500 in the US. About 80 genes were found to be in association with CMT. The phosphoribosyl pyrophosphate synthetase 1 (PRPS1) is an essential enzyme in the primary stage of de novo and salvage nucleotide synthesis. The mutations in the PRPS1 gene leads to X-linked Charcot-Marie-Tooth neuropathy type 5 (CMTX5), PRS super activity, Arts syndrome, X-linked deafness-1, breast cancer, and colorectal cancer. In the present study, we obtained 20 missense mutations from UniProt and dbSNP databases and applied series of comprehensive in silico prediction methods to assess the degree of pathogenicity and stability. In silico tools predicted four missense mutations (D52H, M115 T, L152P, and D203H) to be potential disease causing mutations. We further subjected the four mutations along with native protein to 50 ns molecular dynamics simulation (MDS) using Gromacs package. The resulting trajectory files were analyzed to understand the stability differences caused by the mutations. We used the Root Mean Square Deviation (RMSD), Radius of Gyration (Rg), solvent accessibility surface area (SASA), Covariance matrix, Principal Component Analysis (PCA), Free Energy Landscape (FEL), and secondary structure analysis to assess the structural changes in the protein upon mutation. Our study suggests that the four mutations might affect the PRPS1 protein function and stability of the structure. The proposed study may serve as a platform for drug repositioning and personalized medicine for diseases that are caused by the PRPS1 deficiency.

Substitution impact of highly conserved arginine residue at position 75 in GJB1 gene in association with X-linked Charcot-Marie-tooth disease: A computational study

X-linked Charcot-Marie-Tooth type 1 X (CMTX1) disease is a subtype of Charcot-Marie-Tooth (CMT), which is mainly caused by mutations in the GJB1 gene. It is also known as connexin 32 (Cx32) that leads to Schwann cell abnormalities and peripheral neuropathy. CMTX1 is considered as the second most common form of CMT disease. The aim of this study is to computationally predict the potential impact of different single amino acid substitutions at position 75 of Cx32, from arginine (R) to proline (P), glutamine (Q) and tryptophan (W). This position is known to be highly conserved among the family of connexin. To understand the structural and functional changes due to these single amino acid substitutions, we employed a homology-modeling technique to build the three-dimensional structure models for the native and mutant proteins. The protein structures were further embedded into a POPC lipid bilayer, inserted into a water box, and subjected to molecular dynamics simulation for 50 ns. Our results show that the mutants R75P, R75Q and R75W display variable structural conformation and dynamic behavior compared to the native protein. Our data proves useful in predicting the potential pathogenicity of the mutant proteins and is expected to serve as a platform for drug discovery for patients with CMT.