Exploring the effect of disease causing mutations in metal binding sites of human ARSA in metachromatic leukodystrophy.

The arylsulfatase A (ARSA) gene is observed to be deficient in patients with metachromatic leukodystrophy (MLD), a type of lysosomal storage disease. MLD is a severe neurodegenerative disorder characterized by an autosomal recessive inheritance pattern. This study aimed to map the most deleterious mutations at the metal binding sites of ARSA and the amino acids in proximity to the mutated positions. We utilized an array of computational tools, including PredictSNP, MAPP, PhD-SNP, PolyPhen-1, PolyPhen-2, SIFT, SNAP, and ConSurf, to identify the most detrimental mutations potentially implicated in MLD collected from UniProt, ClinVar, and HGMD. Two mutations, D29N and D30H, as being extremely deleterious based on assessments of pathogenicity, conservation, biophysical characteristics, and stability analysis. The D29 and D30 are located at the metal-interacting regions of ARSA and found to undergo post-translational modification, specifically phosphorylation. Henceforth, the in-depth effect of metal binding upon mutation was examined using molecular dynamics simulations (MDS) before and after phosphorylation. The MDS results exhibited high deviation for the D29N and D30H mutations in comparison to the native, and the same was confirmed by significant residue fluctuation and reduced compactness. These structural alterations suggest that such mutations may influence protein functionality, offering potential avenues for personalized therapeutic and providing a basis for potential mutation-specific treatments for severe MLD patients.

A computational approach to analyzing the functional and structural impacts of Tripeptidyl-Peptidase 1 missense mutations in neuronal ceroid lipofuscinosis.

Neuronal ceroid-lipofuscinosis (NCLs) are a group of severe neurodegenerative conditions, most likely present in infantile, late infantile, juvenile, and adult-onset forms. Their phenotypic characteristics comprise eyesight damage, reduced motor activity and cognitive function, and sometimes tend to die in the initial stage. In recent studies, NCLs have been categorized into at least 14 genetic collections (CLN1-14). CLN2 gene encodes Tripeptidyl peptidase 1 (TPP1), which affects late infantile-onset form. In this study, we retrieved a mutational dataset screening for TPP1 protein from various databases (ClinVar, UniProt, HGMD). Fifty-six missense mutants were enumerated with computational methods to perceive the significant mutants (G475R and G501C) and correlated with clinical and literature data. A structure-based screening method was initiated to understand protein-ligand interaction and dynamic simulation. The docking procedure was performed for the native (3EDY) and mutant (G473R and G501C) structures with Gemfibrozil (gem), which lowers the lipid level, decreases the triglycerides amount in the blood circulation, and controls hyperlipidemia. The Native had an interaction score of -5.57 kcal/mol, and the mutants had respective average binding scores of -6.24 (G473R) and - 5.17 (G501C) kcal/mol. Finally, molecular dynamics simulation showed that G473R and G501C mutants had better flexible and stable orientation in all trajectory analyses. Therefore, this work gives an extended understanding of both functional and structural levels of influence for the mutant form that leads to NCL disorder.

Elucidating the functional impact of G137V and G144R variants in Maroteaux Lamy's Syndrome by Molecular Dynamics Simulation.

Mucopolysaccharidoses VI (Maroteaux Lamy syndrome) is a metabolic disorder due to the loss of enzyme activity of N-acetyl galactosamine-4-sulphatase arising from mutations in the ARSB gene. The mutated ARSB is the origin for the accumulation of GAGs within the lysosome leading to severe growth deformities, causing lysosomal storage disease. The main focus of this study is to identify the deleterious variants by applying bioinformatics tools to predict the conservation, pathogenicity, stability, and effect of the ARSB variants. We examined 170 missense variants, of which G137V and G144R were the resultant variants predicted detrimental to the progression of the disease. The native along with G137V and G144R structures were fixed as the receptors and subjected to Molecular docking with the small molecule Odiparcil to analyze the binding efficiency and the varied interactions of the receptors towards the drug. The interaction resulted in similar docking scores of - 7.3 kcal/mol indicating effective binding and consistent interactions of the drug with residues CYS117, GLN118, THR182, and GLN517 for native, along with G137V and G144R structures. Molecular Dynamics were conducted to validate the stability and flexibility of the native and variant structures on ligand binding. The overall study indicates that the drug has similar therapeutic towards the native and variant based on the higher binding affinity and also the complexes show stability with an average of 0.2 nm RMS value. This can aid in the future development therapeutics for the Maroteaux Lamy syndrome.

Controlling cell proliferation by targeting cyclin-dependent kinase 6 using drug repurposing approach.

Cyclin-dependent kinase 6 (CDK6) is an essential kinase in cell cycle progression, which is a viable target for inhibitors in various malignancies, including breast cancer. This study aimed to virtually screen efficient compounds as new leads in treating breast cancer using a drug repurposing approach. Apoptosis regulatory compounds were taken from the seleckchem database. Molecular docking experiments were carried out in the presence of abemaciclib, a routinely used FDA drug. Compared to conventional drugs, the two compounds demonstrated a higher binding affinity for CDK6. Compounds (N-benzyl-6-[(4-hydroxyphenyl)methyl]-8-(naphthalen-1-ylmethyl)-4,7-dioxo-3,6,9,9a-tetrahydro-2H-pyrazino[1,2-a]pyrimidine-1-carboxamide) and (1'-[4-[1-(4-fluorophenyl)indol-3-yl]butyl]spiro[1H-2-benzofuran-3,4'-piperidine]) were discovered to have an inhibitory effect against CDK6 at -8.49 and -6.78kcal/mol, respectively, compared to -8.09kcal/mol of the control molecule, the interacting residues of these two new compounds were found to fall within the binding site of the CDK6 molecule. Both compounds exhibited equal ADME features compared with abemaciclib and would be well distributed and metabolized by the body with an appropriate druglikeness range. Lastly, molecular dynamics was initiated for 200ns for the selected potent inhibitors and abemaciclib as complexed with CDK6. The RMSD, RMSF, Rg, H-Bond interactions, SASA, PCA, FEL, and MM/PBSA analysis were performed for the complexes to assess the stability, fluctuations, radius of gyration, hydrogen bond interaction, solvent accessibility, essential dynamics, free energy landscape, and MM/PBSA. The selected two compounds are small molecules in the appropriate druglikeness range. The results observed in molecular docking and molecular dynamics simulations were most promising for two compounds, suggesting their potent inhibitory effect against CDK6. We propose that these candidate compounds can undergo in vitro validation and in vivo testing for their further use against cancer.

Deciphering the effect of mutations in MMAA protein causing methylmalonic acidemia-A computational approach.

Methylmalonic acidemia (MMA) is a rare genetic disorder affecting multiple body systems. We aimed to investigate the pathogenic mutations in MMAA that are associated with isolated methylmalonic acidemia to identify the structural behavior of MMAA upon mutation. The algorithms such as PredictSNP, iStable, ConSurf, and Align GVGD were employed to analyze the consequence of the mutations. Molecular docking was carried out for the native MMAA, L89P, G274D, and R359G to interpret its interactions with the GDP substrate. The docked complexes were simulated for 200ns aiding GROMACS in apprehending the behavior of MMAA upon mutation and GDP binding. After simulation, cα disruptions were observed using the RMSF plot, which indicated that several regions of mutant MMAAs have highly fluctuated. The gyration and H-bond plots were used to understand the compactness and intermolecular interaction with the GDP molecule. The MDS analysis showed that the mutations L89P, G274D, and R359G are highly unstable even after GDP binding, with the least compactness, fewer H-bonds, and larger conformational cα motions. Our study provided structural and dynamic insights into MMAA protein, which further helps to characterize these mutants and provide potential treatment strategies for MMA patients.

An integrative analysis to distinguish between emphysema (EML) and alpha-1 antitrypsin deficiency-related emphysema (ADL)-A systems biology approach.

Lung Emphysema is an abnormal enlargement of the air sacs followed by the destruction of alveolar walls without any prominent fibrosis. This study primarily identifies the differentially expressed genes (DEGs), interactions between them, and their significant involvement in the activated signaling cascades. The dataset with ID GSE1122 (five normal lung tissue samples, five of usual emphysema, and five of alpha-1 antitrypsin deficiency-related emphysema) from the gene expression omnibus (GEO) was analyzed using the GEO2R tool. The physical association between the DEGs were mapped using the STRING tool and was visualized in the Cytoscape software. The enriched functional processes were identified with the ClueGO plugin's help from Cytoscape. Further integrative functional annotation was performed by implying the GeneGo Metacore™ to distinguish the enriched pathway maps, process networks, and GO processes. The results from this analysis revealed the critical signaling cascades that have been either activated or inhibited due to identified DEGs. We found the activated pathways such as immune response IL-1 signaling pathway, positive regulation of smooth muscle migration, BMP signaling pathway, positive regulation of leukocyte migration, NIK/NF-kappB signaling, and cytochrome-c oxidase activity. Finally, we mapped four crucial genes (CCL5, ALK, TAC1, CD74, and HLA-DOA) by comparing the functional annotations that could be significantly influential in emphysema molecular pathogenesis. Our study provides insights into the pathogenesis of emphysema and helps in developing potential drug targets against emphysema.

A computational approach to analyse the amino acid variants of GLB1 protein causing GM1 Gangliosidosis.

Lysosomal storage diseases comprise different forms of autosomal recessive disorders from which GM1 gangliosidosis has categorized by the accumulation of complex glycolipids associated with a range of progressive neurologic phenotypes. GM1 gangliosidosis is an inherited disorder that progressively destroys nerve cells (neurons) in the brain and spinal cord. GM1 has three main types of onsets, namely infantile (type I), juvenile (type II), and adult (type III) forms. This study provides a series of computational methods that examine the mutations that occurred in GLB1 protein. Initially, the mutational analysis started with 689 amino acid variants for a sequence-based screening and it was done with quite a few In-silico tools to narrow down the most significant variants by utilizing the standard tools; namely, Evolutionary analysis (77 variants), Pathogenicity prediction (44 variants), Stability predictions (30 variants), Biophysical functions (19 variants) and according to the binding site of protein structure with PDB ID 3THC, seven variants (Y83D, Y83H, Y270S, Y270D, W273R, W273D, and Y333H) were narrowed down. Structure based analysis was performed to understand the interacting profile of the native protein and variants with Miglustat; which is the currently used FDA drug as an alternative to enzyme replacement therapy. Molecular Docking study was done to analyze the protein interaction with Miglustat (ligand), as a result native (3THC) structure had a binding affinity of -8.18 kcal/mol and two variant structures had an average binding affinities of -2.61 kcal/mol (Y83D) and - 7.63 kcal/mol (Y270D). Finally, Molecular Dynamics Simulation was performed to know the mutational activity of the protein structures on Miglustat for 50,000 ps. The Y83D variant showed higher deviation than native protein and Y270D in all trajectory analysis. The analysis was done to the protein structures to check the structural variations happened through simulations. This study aids to understand the most deleterious mutants, the activity of the drug to the protein structure and also gives an insight on the stability of the drug with the native and selected variants.