Molecular mechanisms underlying the impact of mutations in SOD1 on its conformational properties associated with amyotrophic lateral sclerosis as revealed with molecular modelling

Molecular dynamics analysis of superoxide dismutase 1 mutations suggests decoupling between mechanisms underlying ALS onset and progression

Mutations in the superoxide dismutase 1 (SOD1) gene are the second most common known cause of ALS. SOD1 variants express high phenotypic variability and over 200 have been reported in people with ALS. It was previously proposed that variants can be broadly classified in two groups, 'wild-type like' (WTL) and 'metal binding region' (MBR) variants, based on their structural location and biophysical properties. MBR variants, but not WTL variants, were associated with a reduction of SOD1 enzymatic activity. In this study we used molecular dynamics and large clinical datasets to characterise the differences in the structural and dynamic behaviour of WTL and MBR variants with respect to the wild-type SOD1, and how such differences influence the ALS clinical phenotype. Our study identified marked structural differences, some of which are observed in both variant groups, while others are group specific. Moreover, collecting clinical data of approximately 500 SOD1 ALS patients carrying variants, we showed that the survival time of patients carrying an MBR variant is generally longer (∼6 years median difference, p < 0.001) with respect to patients with a WTL variant. In conclusion, our study highlighted key differences in the dynamic behaviour between WTL and MBR SOD1 variants, and between variants and wild-type SOD1 at an atomic and molecular level, that could be further investigated to explain the associated phenotypic variability. Our results support the hypothesis of a decoupling between mechanisms of onset and progression of SOD1 ALS, and an involvement of loss-of-function of SOD1 with the disease progression.

Neuroprotective and nephroprotective effects of Ircinia sponge in polycyclic aromatic hydrocarbons (PAHs) induced toxicity in animal model: a pharmacological and computational approach

The present study investigated the neuroprotective and nephroprotective effects of the sponge Ircinia sp. ethyl acetate extract (ISPE) against persistent aromatic pollutants in vitro and in vivo. Different exponential experimental assays were applied to this study. An in vitro study to investigate the potential therapeutic effect of ISPE using antioxidants (for example, ABTS and DPPH) and anti-Alzheimer assays (inhibition of acetylcholinesterase); the in-vivo study was designed to evaluate the protective effect of ISPE as neuroprotective and nephroprotective against the destructive effect of PAH. Several assays included oxidative assays (LPO), antioxidant biomarkers (GSH, GST), and inflammatory and neurodegenerative biomarkers (PTK,SAA). Additionally, the results were confirmed using histopathological examination. The in silico screening study improved the in vitro and in vivo findings through interaction between the aryl hydrocarbon receptor (AHR) and the polyphenolic content of ISPE extract, which was determined using LCMSM. The results and discussion showed that ISPE exhibited a promising antioxidant and anti-acetylcholinesterase activity as evidenced by IC50 values of 49.74, 28.25, and 0.18 µg/mL in DPPH, ABTS, and acetylcholinesterase inhibition assays, respectively. In vivo, the study showed that animals receiving ISPE before poly aromatic hydrocarbons administration PAHs (Prot, ISPE) showed significant amelioration in kidney functions manifested by the reduction of serum urea, uric acid, and creatinine by 40.6%, 66.4%, and 134.8%, respectively, concerning PAH-injected mice (HAA). Prot, ISPE revealed a decline in malondialdehyde (MDA) and total proteins (TP) in kidney and brain tissues by 73.63% and 50.21%, respectively, for MDA and 59.82% and 80.41%, respectively, for TP with respect to HAA. Prot, ISPE showed significant elevation in reduced glutathione (GSH) and glutathione transferase (GST) in kidney and brain tissues and reduction in the inflammatory and pre-cancerous biomarkers, namely, serum protein tyrosine kinases (PTKs) and serum amyloid A (SAA). These findings were further supported by histopathological examination of kidney and brain tissues, which revealed normal structure approaching normal control. Metabolic profiling of ISPE using LC-MS-MS showed the presence of fourteen polyphenolic compounds belonging mainly to phenolic acids and flavonoids. In silico study revealed that all the tested compounds exerted certain binding with the aryl hydrocarbon receptor, where rutin showed the best fitting (ΔG =  - 7.6 kcal/mol-1) with considerable pharmacokinetic and pharmacodynamic properties revealed from in silico ADME (Absorption, Distribution, Metabolism, and Excretion) study. Hence, it can be concluded that the Ircinia sponge showed a promising protective effect versus kidney and brain toxicity triggered by PAHs.

Amyotrophic lateral sclerosis disease-related mutations disrupt the dimerization of superoxide dismutase 1 - A comparative molecular dynamics simulation study

More than 150 genes are involved in amyotrophic lateral sclerosis (ALS), with superoxide dismutase 1 (SOD1) being one of the most studied. Mutations in SOD1 gene, which encodes the enzyme SOD1 is the second most prevalent and studied cause of familial ALS. SOD1 is a ubiquitous, homodimeric metalloenzyme that forms a critical component of the cellular defense against reactive oxygen species. Several mutations in the SOD1 enzyme cause misfolding, dimerization instability, and increased aggregate formation in ALS. However, there is a lack of information on the dimerization of SOD1 monomers and the mechanistic underpinnings on how the pathogenic mutations disrupt the dimerization mechanism. Here, we presented microsecond-scale molecular dynamics (MD) simulations to unravel how interface-based mutations compromise SOD1 dimerization and provide mechanistic understanding into the corresponding process using WT and three interface-based mutant systems (A4V, T54R, and I113T). Structural stability analysis showed that the mutant systems displayed disparate variations in the catalytic sites which may directly alter the stability and activity of the SOD1 enzyme. Based on the dynamic network analysis and principal component analysis, it has been identified that the mutations weakened the correlated motions along the dimer interface and altered the protein conformational behavior, thus weakening the stability of dimer formation. Moreover, the simulation results identified crucial residues such as G51, D52, G114, I151, and Q153 in establishing the dimerization interaction network, which were weakened or absent in the presence of interfacial mutants. Surface potential analysis on mutant systems also displayed changes in the dimerization potential, thus showing the unfavorable dimer formation. Furthermore, network analysis identified the hotspot residues necessary for SOD1 signal transduction which were surprisingly found in the catalytic sites rather than the anticipated dimerization interface.

Yeast Proteins may Reversibly Aggregate like Amphiphilic Molecules

More than a hundred proteins in yeast reversibly aggregate and phase-separate in response to various stressors, such as nutrient depletion and heat shock. We know little about the protein sequence and structural features behind this ability, which has not been characterized on a proteome-wide level. To identify the distinctive features of aggregation-prone protein regions, we apply machine learning algorithms to genome-scale limited proteolysis-mass spectrometry (LiP-MS) data from yeast proteins. LiP-MS data reveals that 96 proteins show significant structural changes upon heat shock. We find that in these proteins the propensity to phase separate cannot be solely driven by disordered regions, because their aggregation-prone regions (APRs) are not significantly disordered. Instead, the phase separation of these proteins requires contributions from both disordered and structured regions. APRs are significantly enriched in aliphatic residues and depleted in positively charged amino acids. Aggregator proteins with longer APRs show a greater propensity to aggregate, a relationship that can be explained by equilibrium statistical thermodynamics. Altogether, our observations suggest that proteome-wide reversible protein aggregation is mediated by sequence-encoded properties. We propose that aggregating proteins resemble supra-molecular amphiphiles, where APRs are the hydrophobic parts, and non-APRs are the hydrophilic parts.

Computer analysis of the relation between hydrogen bond stability in SOD1 mutants and the survival time of amyotrophic lateral sclerosis patients

BACKGROUND AND OBJECTIVE

Mutations in the SOD1 protein can lead to the death of motor neurons, which, in turn, causes an incurable disease called amyotrophic lateral sclerosis (ALS). At the same time, the mechanism of the onset and development of this disease is not fully understood and is often contradictory.

METHODS

Accelerated Molecular Dynamics as implemented in the OpenMM library, principal component analysis, regression analysis, random forest method.

RESULTS

The stability of hydrogen bonds in 72 mutants of the SOD1 protein was calculated. Principal component analysis was carried out. Based on ten principal components acting as predictors, a multiple linear regression model was constructed. An analysis of the correlation of these ten principal components with the initial values of the stability of hydrogen bonds made it possible to characterize the contribution of known structurally and functionally important sites in the SOD1 to the scatter of ALS patients' survival time.

CONCLUSION

Such an analysis made it possible to put forward hypotheses about the relationship between the stabilizing and destabilizing effects of mutations in different structurally and functionally important regions of SOD1 with the patients's survival time.

Superoxide Dismutase 1 in Health and Disease: How a Frontline Antioxidant Becomes Neurotoxic

Cu/Zn superoxide dismutase (SOD1) is a frontline antioxidant enzyme catalysing superoxide breakdown and is important for most forms of eukaryotic life. The evolution of aerobic respiration by mitochondria increased cellular production of superoxide, resulting in an increased reliance upon SOD1. Consistent with the importance of SOD1 for cellular health, many human diseases of the central nervous system involve perturbations in SOD1 biology. But far from providing a simple demonstration of how disease arises from SOD1 loss-of-function, attempts to elucidate pathways by which atypical SOD1 biology leads to neurodegeneration have revealed unexpectedly complex molecular characteristics delineating healthy, functional SOD1 protein from that which likely contributes to central nervous system disease. This review summarises current understanding of SOD1 biology from SOD1 genetics through to protein function and stability.

Comprehensive in silico analysis and molecular dynamics of the superoxide dismutase 1 (SOD1) variants related to amyotrophic lateral sclerosis

Amyotrophic Lateral Sclerosis (ALS) is the most frequent motor neuron disorder, with a significant social and economic burden. ALS remains incurable, and the only drugs approved for its treatments confers a survival benefit of a few months for the patients. Missense mutations in superoxide dismutase 1 (SOD1), a major cytoplasmic antioxidant enzyme, has been associated with ALS development, accounting for 23% of its familial cases and 7% of all sporadic cases. This work aims to characterize in silico the structural and functional effects of SOD1 protein variants. Missense mutations in SOD1 were compiled from the literature and databases. Twelve algorithms were used to predict the functional and stability effects of these mutations. ConSurf was used to estimate the evolutionary conservation of SOD1 amino-acids. GROMACS was used to perform molecular dynamics (MD) simulations of SOD1 wild-type and variants A4V, D90A, H46R, and I113T, which account for approximately half of all ALS-SOD1 cases in the United States, Europe, Japan, and United Kingdom, respectively. 233 missense mutations in SOD1 protein were compiled from the databases and literature consulted. The predictive analyses pointed to an elevated rate of deleterious and destabilizing predictions for the analyzed variants, indicating their harmful effects. The ConSurf analysis suggested that mutations in SOD1 mainly affect conserved and possibly functionally essential amino acids. The MD analyses pointed to flexibility and essential dynamics alterations at the electrostatic and metal-binding loops of variants A4V, D90A, H46R, and I113T that could lead to aberrant interactions triggering toxic protein aggregation. These alterations may have harmful implications for SOD1 and explain their association with ALS. Understanding the effects of SOD1 mutations on protein structure and function facilitates the design of further experiments and provides relevant information on the molecular mechanism of pathology, which may contribute to improvements in existing treatments for ALS.

Learning the changes of barnase mutants thermostability from structural fluctuations obtained using anisotropic network modeling

In biotechnology applications, rational design of new proteins with improved physico-chemical properties includes a number of important tasks. One of the greatest practical and fundamental challenges is the design of highly thermostable protein enzymes that maintain catalytic activity at high temperatures. This problem may be solved by introducing mutations into the wild-type enzyme protein. In this work, to predict the impact of such mutations in barnase protein we applied the anisotropic network modeling approach, revealing atomic fluctuations in structural regions that are changed in mutants compared to the wild-type protein. A regression model was constructed based on these structural features that can allow one to predict the thermal stability of new barnase mutants. Moreover, the analysis of regression model provides a mechanistic explanation of how the structural features can contribute to the thermal stability of barnase mutants.

Critical Sites of DNA Backbone Integrity for Damaged Base Removal by Formamidopyrimidine-DNA Glycosylase

DNA glycosylases, the enzymes that initiate base excision DNA repair, recognize damaged bases through a series of precisely orchestrated movements. Most glycosylases sharply kink the DNA axis at the lesion site and extrude the target base from the DNA double helix into the enzyme's active site. Little attention has been paid so far to the role of the physical continuity of the DNA backbone in allowing the required conformational distortion. Here, we analyze base excision by formamidopyrimidine-DNA glycosylase (Fpg) from substrates keeping all phosphates but containing a nick within three nucleotides of the lesion in either DNA strand. Four phosphoester linkages at the damaged nucleotide and two nucleotides 3' to it were essential for Fpg activity, while the breakage of the others, even at the same critical phosphates, had no effect or even stimulated the reaction. Reduction of the likelihood of hydrogen bonding at the nicks by using dideoxynucleotides as their 3'-terminal groups was more detrimental for the activity. All phosphoester bonds in the complementary strand were dispensable for base excision, but nicks close to the orphaned nucleotide caused early termination of damaged strand cleavage. Elastic network analysis of Fpg-DNA structures showed that the vibrational motions of the critical phosphates are strongly correlated, in part due to the presence of the protein. Overall, our results suggest that mechanical forces propagating along the DNA backbone play a critical role in the correct conformational distortion of DNA by Fpg and possibly by other target base-everting DNA glycosylases.

Improved regression model to predict an impact of SOD1 mutations on ALS patients survival time based on analysis of hydrogen bond stability

Amyotrophic lateral sclerosis (ALS) is an incurable neurodegenerative disease characterised by the inevitable degeneration of central and peripheral motor neurons. Aggregation of mutant SOD1 is one of the molecular mechanisms underlying the onset of the disease. There are a number of regression models designed to predict the survival of patients based on an analysis of experimental data on thermostability, heterodimerisation energy, and changes in the hydrophobicity of SOD1 mutants. Previously, we proposed regression models linking the change in the stability of hydrogen bonds in mutant SOD1 calculated using molecular dynamics and elastic networks with patients survival time. In this study, these models were improved in terms of accuracy of survival time prediction by taking into account the variance of survival time values relative to the mean, the number of patients carrying each specific mutation, and the use of random forest regression as a regression method. The accuracy of the previous models was roughly 5.2 years while the accuracy of the new ones are up to 4 years. The model is also superior to those published by other authors. It was found that the hydrogen bonds important for prediction of survival time are formed by residues at positions located in the regions of the protein responsible for aggregation as well as in structural and functionally important sites.