Search

The mechanisms responsible for the pathogenesis and progression of Amyotrophic Lateral Sclerosis (ALS) remain poorly understood, making the diagnosis of ALS challenging. We aimed to find the novel gene biomarkers via computationally analyzing microarray expression studies, in three different cell lineages, namely myotube cells, astrocyte cells and oligodendrocyte cells. Microarray gene expression profiles were obtained and analyzed for three cell types: myotube cell lineage (GSE122261), astrocyte, and oligodendrocyte cell lineage (GSE87385). A comprehensive computational pipeline, tailored explicitly for microarray gene expression profiling studies, was devised to analyze the sample groups, wherein the myotube sample group comprised of six control (GSM3462697, GSM3462698, GSM3462699, GSM3462700, GSM3462701, GSM3462702) & six diseased (GSM3462691, GSM3462692, GSM3462693, GSM3462694, GSM3462695, GSM3462696) samples were considered. Similarly, for the astrocyte sample group two samples each for the control (GSM2330040, GSM2330042) and the diseased (GSM2330039, GSM2330041), and for the oligodendrocyte sample group, 2 control (GSM2330043, GSM2330045) samples and two diseased (GSM2330044, GSM2330046) samples were considered for the current study. The in-depth interaction of these DEGs was studied using MCODE and subjected to preliminary functional analysis using ClueGO/CluePedia plug-in. Qiagen's IPA software was employed for enrichment analysis, which generated the key canonical pathways and a list of potential biomarker molecules specific to each sample group. The preliminary analysis yielded 512 DEGs across all 3-sample groups, wherein 139 DEGs belonged to the myotube sample group, 216 DEGs for the astrocyte sample group, and 157 DEGs for the oligodendrocytes sample group. The data suggests growth hormone signaling and its activity, ErbB signaling pathway, and JAK/STAT signaling pathway are some of the pathways that are significantly dysregulated and play a crucial role in the development and progression of ALS. KISS1R and CSHL1 are potential genes that could act as diagnostic biomarkers in myotube cell types. Also, KRAS, TGFB2, JUN, and SMAD6 genes may be used as prognostic biomarkers to differentiate between early and late-stage ALS-diseased patients.

The number of confirmed COVID-19 cases is rapidly increasing with no direct treatment for the disease. Few repurposed drugs, such as Remdesivir, Chloroquine, Hydroxychloroquine, Lopinavir, and Ritonavir, are being tested against SARS-CoV-2. Remdesivir is the drug of choice for Ebola virus disease and has been authorized for emergency use. This drug acts against SARS-CoV-2 by inhibiting the RNA-dependent-RNA-polymerase (RdRp) of SARS-CoV-2. RdRp of viruses is prone to mutations that confer drug resistance. A recent study by Pachetti et al. in 2020 identified the P323L mutation in the RdRp protein of SARS-CoV-2. In this study, we aimed to determine the potency of lead compounds similar to Remdesivir, which can be used as an alternative when variants of SARS-CoV-2 develop resistance due to RdRp mutations. The initial screening yielded 704 compounds that were 90% similar to the control drug, Remdesivir. On further evaluation through drugability and antiviral inhibition percentage analyses, we shortlisted 32 and seven compounds, respectively. These seven compounds were further analyzed for their molecular interactions, which revealed that all seven compounds interacted with RdRp with higher affinity than Remdesivir under native conditions. However, three compounds failed to interact with the mutant protein with higher affinity than Remdesivir. Dynamic cross-correlation matrix (DCCM) and vector field collective motions analyses were performed to identify the precise movements of docked complexes' residues. Furthermore, the compound SCHEMBL20144212 showed a high affinity for native and mutant proteins and might provide an alternative against SARS-CoV-2 variants that might confer resistance to Remdesivir. Further validations by in vitro and in vivo studies are needed to confirm the efficacy of our lead compounds for their inhibition against SARS-CoV-2.

Sjögren-Larsson syndrome (SLS) is an autoimmune disorder inherited in an autosomal recessive pattern. To date, 80 missense mutations have been identified in association with the Aldehyde Dehydrogenase 3 Family Member A2 (ALDH3A2) gene causing SLS. Disruption of the function of ALDH3A2 leads to excessive accumulation of fat in the cells, which interferes with the normal function of protective membranes or materials that are necessary for the body to function normally. We retrieved 54 missense mutations in the ALDH3A2 from the OMIM, UniProt, dbSNP, and HGMD databases that are known to cause SLS. These mutations were examined with various in silico stability tools, which predicted that the mutations p.S308N and p.R423H that are located at the protein-protein interaction domains are the most destabilizing. Furthermore, to determine the atomistic-level differences within the protein-protein interactions owing to mutations, we performed macromolecular simulation (MMS) using GROMACS to validate the motion patterns and dynamic behavior of the biological system. We found that both mutations (p.S380N and p.R423H) had significant effects on the protein-protein interaction and disrupted the dimeric interactions. The computational pipeline provided in this study helps to elucidate the potential structural and functional differences between the ALDH3A2 native and mutant homodimeric proteins, and will pave the way for drug discovery against specific targets in the SLS patients.

Blau syndrome (BS), which affects the eyes, skin, and joints, is an autosomal dominant genetic inflammatory disorder. BS is caused by mutations in the NOD2 gene. However, there are no direct treatments, and treatment with conventional anti-inflammatory drugs such as adrenal glucocorticoids, anti-metabolites, and biological agents such as anti-TNF and infliximab have all been attempted with varying degrees of success. In this study, we tried to identify all the reported mutations in the NOD2 protein that cause BS. Collectively, 114 missense mutations were extracted from the UniProt, ClinVar, and HGMD databases. The mutations were further subjected to pathogenic, stability, and conservation analyses. According to these computational analyses, six missense mutations (R334Q, R334W, E383G, E383K, R426H, and T605P) were found to be highly deleterious, destabilizing, and positioned in the conserved position. ADP to ATP conversion plays a crucial role in switching the closed-form of NOD2 protein to the open-form, thus activating the protein. Accordingly, the mutations in the ADP binding sites have received more attention in comparison to the mutations in the non-ADP binding positions. Interestingly, the W490L mutation is positioned in the ADP binding site and exhibits highly deleterious and destabilizing properties. Additionally, W490L was also found to be conserved, with a ConSurf score of 7. Therefore, we further performed homology modeling to determine the 3D structure of native NOD2 and the W490L mutant. Molecular docking analysis was carried out to understand the change in the interaction of ADP with the NOD2 protein. We observed that ADP had a stronger interaction with the native NOD2 protein compared to the W490L mutant. Finally, ADP complexed with native NOD2 and W490L mutant were subjected to molecular dynamics simulations, and the trajectories were analyzed. In the simulations, we observed decreased deviation and fluctuations in native NOD2, whereas decreased compactness and inter- and intramolecular hydrogen bonds were observed in the W490L mutant. This study is expected to serve as a platform for developing targeted drug therapy for BS.

Type I galactosemia is a very rare autosomal recessive genetic metabolic disorder that occurs because of the mutations present in the galactose-1-phosphate uridyl transferase (GALT) gene, resulting in a deficiency of the GALT enzyme. The action of the GALT enzyme is to convert galactose-1-phosphate and uridine diphosphate glucose into glucose-1-phosphate (G1P) and uridine diphosphate-galactose, a crucial second step of the Leloir pathway. A missense mutation in the GALT enzyme leads to variable galactosemia's clinical presentations, ranging from mild to severe. Our study aimed to employ a comprehensive computational pipeline to analyze the most prevalent missense mutations (p.S135L, p.K285 N, p.Q188R, and p.N314D) responsible for galactosemia; these genes could serve as potential targets for chaperone therapy. We analyzed the four mutations through different computational analyses, including amino acid conservation, in silico pathogenicity and stability predictions, and macromolecular simulations (MMS) at 50 ns The stability and pathogenicity predictors showed that the p.Q188R and p.S135L mutants are the most pathogenic and destabilizing. In agreement with these results, MMS analysis demonstrated that the p.Q188R and p.S135L mutants possess higher deviation patterns, reduced compactness, and intramolecular H-bonds of the protein. This could be due to the physicochemical modifications that occurred in the mutants p.S135L and p.Q188R compared to the native. Evolutionary conservation analysis revealed that the most prevalent mutations positions were conserved among different species except N314. The proposed research study is intended to provide a basis for the therapeutic development of drugs and future treatment of classical galactosemia and possibly other genetic diseases using chaperone therapy.

Alcohol use disorder (AUD) is a multifactorial psychiatric behavior disorder. Disulfiram is the first approved drug by the Food and Drug Administration for alcohol-dependent patients, which targets the ALDH2 enzyme. Several genes are known to be involved in alcohol metabolism; mutations in any of these genes are known to be associated with AUD. The E504K mutation in the ALDH2 of the precursor protein or the E487K of the mature protein (E504K/E487K; ALDH2*2 allele) is carried by approximately 8% of the world population. In this study, we aimed to test the known inactive allele ALDH2*2, to validate the use of our extensive computational pipeline (in silico tools, molecular modeling, and molecular docking) for testing the interaction between the ALDH2*2 allele, NAD+, and Disulfiram. In silico predictions showed that the E504K variant of ALDH2 to be pathogenic and destabilizing with the maximum number of prediction in silico tools. Consequently, we studied the effect of this mutation mainly on the interaction between NAD+ -E504K and Disulfiram-E504K complexes using molecular docking technique, and molecular dynamics (MD) analysis. From the molecular docking analysis with NAD+ , we observed that the interaction affinity of the NAD+ decreases with the impact of E504K variant. On the other hand, the drug Disulfiram showed similar interaction in both the native and mutant ALDH2 proteins. Further, the comprehensive MD analysis predicted that the E504K destabilizes the protein and influences the NAD+ and Disulfiram interactions. Our findings reveal that the interaction of NAD+ to the protein is disturbed by the E504K/E487K variant whereas the drug Disulfiram has a similar effect as both native ALDH2 and ALDH2 bearing E504K/E487K variant. This study provides a platform to understand the effect of E504K/E487K on the molecular interaction with NAD+ and Disulfiram.

The DNA repair system is crucial to repair the error resulting in DNA replication. MSH2-MSH6 protein complex plays a significant role in maintaining the mismatch repair mechanism. Mutations in the interface between the two proteins compromise their function in the repair process. The present study aims to understand the impact of missense mutations in the interacting sites of the MSH2-MSH6 protein complex. MSH6 is unstable due to the disordered N-terminal domain. This is stabilized by the MSH2 hetero-dimerization. We used pathogenicity and stability predictors to identify the missense mutations that could be more pathogenic with the destabilizing property. The mutations W764C of MSH2, and L1201F and G1316E of MSH6 were predicted to be highly deleterious and destabilizing by all the in silico predictors. The dynamic motion of the native and mutant (W764C) MSH2-MSH6 protein complexes was further investigated using Molecular Dynamics Simulations of the GROMACS package. The Root Mean Square Deviation (RMSD), Radius of Gyration (Rg), and change in a number of intramolecular hydrogen bonds (H-bonds) were analyzed using the embedded packages of GROMACS. From the simulation studies, we observed higher deviation, lower protein compactness, and a decrease in the number of intramolecular hydrogen bonds in the mutant W764C MSH2-MSH6 protein complex. The observed results from the computational methods suggest the involvement of higher structural impact on the MSH2-MSH6 protein complex upon W764C mutation could affect the DNA repair mechanism.

DNA repair mechanism is a process through which the cell repairs its damaged DNA. Although there are several mechanisms involved in the DNA repair mechanisms, the direct reversal method is the simplest and does not require a reference template, in which the guanine bases are often methylated, and the methyl guanine methyl transferase protein (MGMT) reverses them. The mutations occurring in the MGMT protein might result in dysfunction of such DNA repair mechanism. In this study, we attempted to evaluate the impact of six missense mutations (Y114E, Y114A, R128G, R128A, R128K, and C145A) at three active-site positions (Y114, C145, and R128) as this might hinder the DNA binding to the protein. These six mutations were subjected to pathogenicity, stability, and conservation analysis using online servers such as PredictSNP, iStable, and ConSurf, respectively. From the predictions, all the six mutations were almost predicted to be significant. Considering true positives, true negatives, false positives, and false negatives, three mutations (Y114E, R128G, and C145A) showed "loss of DNA repair activity," and were analyzed further using molecular dynamics simulations (MDS) using GROMACS for 50ns. MDS run showed that the C145A mutant demonstrated higher structural deviation, decreased compactness, and the binding patterns. The Y114E mutant showed almost a null effect from the structural analysis. Finally, the R128G mutant showed structural variations in between the C145A and Y114E mutations of MGMT protein. We believe that the observed findings in this computational approach might further pave a way of providing better treatment measures by understanding the DNA repair mechanisms.