Dynamic properties of SOD1 mutants can predict survival time of patients carrying familial amyotrophic lateral sclerosis

Advances in the understanding of protein misfolding and aggregation through molecular dynamics simulation

Aberrant protein folding known as protein misfolding is counted as one of the striking factors of neurodegenerative diseases. The extensive range of pathologies caused by protein misfolding, aggregation and subsequent accumulation are mainly classified into either gain of function diseases or loss of function diseases. In order to seek for novel strategies for treatment and diagnosis of neurodegenerative diseases, insights into the mechanism of misfolding and aggregation is essential. A comprehensive knowledge on the factors influencing misfolding and aggregation is required as well. An extensive experimental study on protein aggregation is somewhat challenging due to the insoluble and noncrystalline nature of amyloid fibrils. Thus there has been a growing use of computational approaches including Monte Carlo simulation, docking simulation, molecular dynamics simulation in the study of protein misfolding and aggregation. The review presents a discussion on molecular dynamics simulation alone as to how it has emerged as a promising tool in the understanding of protein misfolding and aggregation in general, detailing upon three different aspects considering four misfold prone proteins in particular. It is noticeable that all four proteins considered in this review i.e prion, superoxide dismutase1, huntingtin and amyloid β are linked to chronic neurodegenerative diseases with debilitating effects. Initially the review elaborates on the factors influencing the misfolding and aggregation. Next, it addresses our current understanding of the amyloid structures and the associated aggregation mechanisms, finally, summarizing the contribution of this computational tool in the search for therapeutic strategies against the respective protein-deposition diseases.

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.

Studying ALS: Current Approaches, Effect on Potential Treatment Strategy

Amyotrophic lateral sclerosis (ALS) is one of the most common neurodegenerative diseases, characterized by inevitable progressive paralysis. To date, only two disease modifying therapeutic options are available for the patients with ALS, although they show very modest effect on disease course. The main reason of failure in the field of pharmacological correction of ALS is inability to untangle complex relationships taking place during ALS initiation and progression. Traditional methods of research, based on morphology or transgenic animal models studying provided lots of information about ALS throughout the years. However, translation of these results to humans was unsuccessful due to incomplete recapitulation of molecular pathology and overall inadequacy of the models used in the research.In this review we summarize current knowledge regarding ALS molecular pathology with depiction of novel methods applied recently for the studies. Furthermore we describe present and potential treatment strategies that are based on the recent findings in ALS disease mechanisms.

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.

Tryptophan 32 mediates SOD1 toxicity in a in vivo motor neuron model of ALS and is a promising target for small molecule therapeutics

SOD1 misfolding, toxic gain of function, and spread are proposed as a pathological basis of amyotrophic lateral sclerosis (ALS), but the nature of SOD1 toxicity has been difficult to elucidate. Uniquely in SOD1 proteins from humans and other primates, and rarely in other species, a tryptophan residue at position 32 (W32) is predicted to be solvent exposed and to participate in SOD1 misfolding. We hypothesized that W32 is influential in SOD1 acquiring toxicity, as it is known to be important in template-directed misfolding. We tested if W32 contributes to SOD1 cytotoxicity and if it is an appropriate drug target to ameliorate ALS-like neuromuscular deficits in a zebrafish model of motor neuron axon morphology and function (swimming). Embryos injected with human SOD1 variant with W32 substituted for a serine (SOD1^W32S) had reduced motor neuron axonopathy and motor deficits compared to those injected with wildtype or disease-associated SOD1. A library of FDA-approved small molecules was ranked with virtual screening based on predicted binding to W32, and subsequently filtered for analogues using a pharmacophore model based on molecular features of the uracil moiety of a small molecule previously predicted to interact with W32 (5'-fluorouridine or 5'-FUrd). Along with testing 5'-FUrd and uridine, a lead candidate from this list was selected based on its lower toxicity and improved blood brain barrier penetrance; telbivudine significantly rescued SOD1 toxicity in a dose-dependent manner. The mechanisms whereby the small molecules ameliorated motor neuron phenotypes were specifically mediated through human SOD1 and its residue W32, because these therapeutics had no measurable impact on the effects of UBQLN4^D90A, EtOH, or tryptophan-deficient human SOD1^W32S. By substituting W32 for a more evolutionarily conserved residue (serine), we confirmed the significant influence of W32 on human SOD1 toxicity to motor neuron morphology and function; further, we performed pharmaceutical targeting of the W32 residue for rescuing SOD1 toxicity. This unique residue offers future novel insights into SOD1 stability and toxic gain of function, and therefore poses an potential target for drug therapy.

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.

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

Background

So far, little is known about the molecular mechanisms of amyotrophic lateral sclerosis onset and progression caused by SOD1 mutations. One of the hypotheses is based on SOD1 misfolding resulting from mutations and subsequent deposition of its cytotoxic aggregates. This hypothesis is complicated by the fact that known SOD1 mutations of similar clinical effect could be distributed over the whole protein structure.

Results

In this work, a measure of hydrogen bond stability in conformational states was studied with elastic network analysis of 35 SOD1 mutants. Twenty-eight hydrogen bonds were detected in nine of 35 mutants with their stability being significantly different from that with the wild-type. These hydrogen bonds were formed by the amino acid residues known from the literature to be located in contact between SOD1 aggregates. Additionally, residues disposed between copper binding sites of both protein subunits were found from the models to form a stiff core, which can be involved in mechanical impulse transduction between these active centres.

Conclusions

The modelling highlights that both stability of the copper binding site and stability of the dimer can play an important role in ALS progression.

Electronic supplementary material

The online version of this article (10.1186/s12900-018-0080-9) contains supplementary material, which is available to authorized users.