Comparative structural and conformational studies on H43R and W32F mutants of copper-zinc superoxide dismutase by molecular dynamics simulation

Mechanistic insights into the conformational switch in profilin-1 subject to collective effects of mutation and histidine tautomerism

Mutations and histidine (His) tautomerism in profilin-1 (PFN1) are associated with amyotrophic lateral sclerosis (ALS). The conformational changes in PFN1 caused by the collective effects of G117V mutation and His tautomeric isomers εε, εδ, δε, and δδ were clarified using molecular dynamics (MD) simulations. The predominant structural variations were seen in α-helices, β-sheets, turns, and coils and the His tautomer's unique degree of disruption was seen in these conformations. The content of α-helices was 23.2 % in the εε and δδ isomers, but the observed α-helices in the isomers εδ and δε were 20.3 % and 21.7 % respectively. The percentage of β-sheet was found to be higher (34.1) in the εε isomer than in the εδ, δε, and δδ isomers, and the values were 30.4, 29.7, and 31.9, respectively. Intermolecular water dynamics analysis discloses that His 133 can form an intramolecular H-bond interaction (N^α-H---N^δ), confirming the experimental observations in the simulations of εε, δε, and δδ isomers of G117V PFN1 mutant. It was concluded that these solvent molecules are crucial for aggregation and must be considered in future research on the PFN1 associated with ALS. Overall, the study offers a thorough microscopic understanding of the pathogenic mechanisms behind conformational changes that cause aggregation illnesses like ALS.

The metal cofactor zinc and interacting membranes modulate SOD1 conformation-aggregation landscape in an in vitro ALS model

Aggregation of Cu-Zn superoxide dismutase (SOD1) is implicated in the motor neuron disease, amyotrophic lateral sclerosis (ALS). Although more than 140 disease mutations of SOD1 are available, their stability or aggregation behaviors in membrane environment are not correlated with disease pathophysiology. Here, we use multiple mutational variants of SOD1 to show that the absence of Zn, and not Cu, significantly impacts membrane attachment of SOD1 through two loop regions facilitating aggregation driven by lipid-induced conformational changes. These loop regions influence both the primary (through Cu intake) and the gain of function (through aggregation) of SOD1 presumably through a shared conformational landscape. Combining experimental and theoretical frameworks using representative ALS disease mutants, we develop a 'co-factor derived membrane association model' wherein mutational stress closer to the Zn (but not to the Cu) pocket is responsible for membrane association-mediated toxic aggregation and survival time scale after ALS diagnosis.

An Allosteric Pathway in Copper, Zinc Superoxide Dismutase Unravels the Molecular Mechanism of the G93A Amyotrophic Lateral Sclerosis-Linked Mutation

Several different mutations of the protein copper, zinc superoxide dismutase (SOD1) produce the neurodegenerative disorder amyotrophic lateral sclerosis (ALS). The molecular mechanism by which the diverse mutations converge to a similar pathology is currently unknown. The electrostatic loop (EL) of SOD1 is known to be affected in all of the studied ALS-linked mutations of SOD1. In this work, we employ a multiscale simulation approach to show that this perturbation corresponds to an increased probability of the EL detaching from its native position, exposing the metal site of the protein to water. From extensive atomistic and coarse-grained molecular dynamics (MD) simulations, we identify an allosteric pathway that explains the action of the distant G93A mutation on the EL. Finally, we employ quantum mechanics/molecular mechanics MD simulations to show that the opening of the EL decreases the Zn(II) affinity of the protein. As the loss of Zn(II) is at the center of several proposed pathogenic mechanisms in SOD1-linked ALS, the structural effect identified here not only is in agreement with the experimental data but also places the opening of the electrostatic loop as the possible main pathogenic effect for a significant number of ALS-linked SOD1 mutations.

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.

Molecular dynamics of a far positioned SOD1 mutant V14M reveals pathogenic misfolding behavior

Human superoxide dismutase (Cu/Zn SOD1) is a homodimeric enzyme. Mutations in Cu/Zn SOD1 causes a familial form of amyotrophic lateral sclerosis (fALS), and aggregation of mutant SOD1 has been proposed to play a role in neurodegeneration. Though a majority of the mutations are point substitutions, there are a few changes that result in amino acid deletions or truncations of the polypeptide. These pathogenic mutations are scattered throughout the three-dimensional structure of the dimeric enzyme, which creates a puzzling pattern to investigate the molecular determinants of fALS. The most common hypothesis proposed that the misfolding of SOD1 mutants are primarily triggered by decreased affinity for metal ions. However, this hypothesis is challenging, as a significant number of disease-causing mutations are located far away from the metal-binding site and dimer interface. So in the present study, we have investigated the influence of such a far positioned pathogenic mutation, V14M, in altering the stability and folding of the Cu/Zn SOD1. Though the location of Val14 is far positioned, it has a vital role in the stability of SOD1 by preserving its hydrophobic cluster at one end of the β barrel domain. We have performed MD simulations of the V14M mutant for 80 ns timescale. The results reveal the fact that irrespective of its location, V14M mutation triggers a conformational change that is more similar to that of the metal-deficient holo form and could resemble an intermediate state in the folding reaction which results in protein misfolding and aggregation.

Temperature-dependent conformational dynamics of cytochrome c: Implications in apoptosis

Heat, electric shock, and burn injuries induce apoptosis by releasing cytochrome c (cyt-c) from mitochondria and by subsequently activating the death protease, caspases-3. During apoptosis, cyt-c undergoes changes in the secondary structure that have been suggested to increase its peroxidase activity. Information about these structural changes will provide better understanding of the apoptotic mechanism. Hence, temperature-dependent conformational dynamics of cyt-c has been investigated through molecular dynamics (MD) simulations to explain the structural changes and to correlate them with its apoptotic behavior. We observe that, at lower temperatures (223, 248, and 300K), the secondary structure of cyt-c, remains stable, while at higher temperatures (323, 373, 423, and 473K), the secondary structural regions change significantly. Further, our MD results indicate that these structural changes are mainly localized on α-helices, turns, β-sheets, and important loops that were involved in the stabilization of the heme conformation. This conformational transition between specific regions of secondary structure of cyt-c directly affects the electron tunneling properties of the proteins as observed experimentally. We quantify and compare these changes and explain that the temperature plays a vital role in assuring the structural stability of cyt-c and thus its functions. Our findings from this MD study reproduce experimental results at high temperatures and provide evidence for the alteration of the heme through the disruption of the H-bonding interactions between specific regions of cyt-c, thereby enhancing its peroxidase activity which plays a crucial role in the apoptotic process.

Destabilization of the metal site as a hub for the pathogenic mechanism of five ALS-linked mutants of copper, zinc superoxide dismutase

Amyotrophic lateral sclerosis (ALS) is a lethal neurodegenerative disease, with no effective pharmacological treatment. Its pathogenesis is unknown, although a subset of the cases is linked to genetic mutations. A significant fraction of the mutations occur in one protein, copper, zinc superoxide dismutase (SOD1). The toxic function of mutant SOD1 has not been elucidated, but damage to the metal site of the protein is believed to play a major role. In this work, we study the electrostatic loop of SOD1, which we had previously proposed to work as a "solvent seal" isolating the metal site from water molecules. Out of the five contact points identified between the electrostatic loop and its dock in the rest of the protein, three points were found to be affected by ALS-linked mutations, with a total of five mutations identified. The effect of the five mutations was studied using methods of computational chemistry. We found that four of the mutations destabilize the proposed solvent seal, while the fifth mutation directly affects the metal-site stability. In the two contact points unaffected by ALS-linked mutations, the side chains of the residues were not found to play a stabilizing role. Our results show that the docking of the electrostatic loop to the rest of SOD1 plays a role in ALS pathogenesis, in support of that structure acting as a solvent barrier for the metal site. The results provide a unified pathogenic mechanism for five different ALS-linked mutations of SOD1.

Investigation of structural dynamics of Thrombocytopenia Cargeeg mutants of human apoptotic cytochrome c: A molecular dynamics simulation approach

Naturally occurring mutations to cytochrome c (cyt-c) have been identified recently in patients with mild autosomal dominant thrombocytopenia (low platelet levels), which yield cyt-c mutants with enhanced apoptotic activity. However, the molecular mechanism underlying this low platelet production and enhanced apoptosis remain unclear. Therefore, an attempt is made herein for the first time to investigate the effects of mutations of glycine 41 by serine (G41S) and tyrosine 48 by histidine (Y48H) on the conformational and dynamic changes of apoptotic (Fe³⁺) cyt-c using all atom molecular dynamics (MD) simulations in explicit water solvent. Our 30ns MD simulations demonstrate considerable structural differences in G41S and Y48H compared to wild type (WT) cyt-c, such as increasing distances between the critical electron transfer residues results in open conformation at the heme active site, large fluctuations in β-turns and α-helices. Additionally, although the β-sheets remain mostly unaffected in all the three cyt-c simulations, the α-helices undergo conformational switch to β-turns in both the mutant simulations. Importantly, this conformational switch of α-helix to β-turn around heme active site should attributes to the loss of intraprotein H-bonds in the mutant simulations especially between NE2 (His26) and O (Pro44) in agreement with the experimental report. Further, essential dynamics analysis reveals that overall motions of WT cyt-c is mainly involved only in the first eigenvector, but in G41S and Y48H the overall motions are mainly in three and two eigenvectors respectively. Overall, the detailed atomistic level information provide a unifying description for the molecular mechanism of structural destabilization, disregulation of platelet formation and enhanced peroxidase activity of the mutant cyt-c's in the pathology of intrinsic apoptosis.

Essential dynamics analysis captures the concerted motion of the integrin-binding site in jerdostatin, an RTS disintegrin

Disintegrins, small molecular weight proteins contained in the venom of vipers and rattlesnakes, are high-affinity and selectivity integrin antagonists. Disintegrins inhibitory epitope mainly consists in a tripeptide sequence localized in a mobile loop protruding from the protein core. RTS and/or KTS tripeptide characterizes the most recently discovered group of disintegrins that selectively block α1β1 integrin receptor. A NMR study dedicated to structure and dynamics properties of jerdostatin, an RTS disintegrin, demonstrated that the substitution of the native RTS with KTS motif impaired flexibility and inhibitory activity of the molecule. Here we add atomic details to the experimental profiles of jerdostatin and its R24K mutant by analyzing the dynamics behavior of the molecules through computational methods. For jerdostatin wild type, molecular dynamics simulations and essential dynamics analyses showed that Y31 residue acts as hinge element in the concerted motions involving the active loop and the C-terminal tail. R24 side chain ability to engage both cation-π and H-bond interactions with Y31 residue was found crucial for that breathing mechanism. Less significant loop-tail concerted motions were observed for the R24K mutant. The description at atomic resolution of jerdostatin dynamics is useful for decoding the influence of specific residues on disintegrin functional properties.

Critical Nucleus Structure and Aggregation Mechanism of the C-terminal Fragment of Copper-Zinc Superoxide Dismutase Protein

The aggregation of the copper-zinc superoxide dismutase (SOD1) protein is linked to familial amyotrophic lateral sclerosis, a progressive neurodegenerative disease. A recent experimental study has shown that the (147)GVIGIAQ(153) SOD1 C-terminal segment not only forms amyloid fibrils in isolation but also accelerates the aggregation of full-length SOD1, while substitution of isoleucine at site 149 by proline blocks its fibril formation. Amyloid formation is a nucleation-polymerization process. In this study, we investigated the oligomerization and the nucleus structure of this heptapeptide. By performing extensive replica-exchange molecular dynamics (REMD) simulations and conventional MD simulations, we found that the GVIGIAQ hexamers can adopt highly ordered bilayer β-sheets and β-barrels. In contrast, substitution of I149 by proline significantly reduces the β-sheet probability and results in the disappearance of bilayer β-sheet structures and the increase of disordered hexamers. We identified mixed parallel-antiparallel bilayer β-sheets in both REMD and conventional MD simulations and provided the conformational transition from the experimentally observed parallel bilayer sheets to the mixed parallel-antiparallel bilayer β-sheets. Our simulations suggest that the critical nucleus consists of six peptide chains and two additional peptide chains strongly stabilize this critical nucleus. The stabilized octamer is able to recruit additional random peptides into the β-sheet. Therefore, our simulations provide insights into the critical nucleus formation and the smallest stable nucleus of the (147)GVIGIAQ(153) peptide.

Interaction between dimer interface residues of native and mutated SOD1 protein: a theoretical study

Cu-Zn superoxide dismutase 1 (SOD1) is a highly conserved bimetallic protein enzyme, used for the scavenging the superoxide radicals (O2 (-)) produced due to aerobic metabolism in the mitochondrial respiratory chain. Over 100 mutations have been identified and found to be in the homodimeric structure of SOD1. The enzyme has to be maintained in its dimeric state for the structural stability and enzymatic activity. From our investigation, we found that the mutations apart from the dimer interface residues are found to affect the dimer stability of protein and hence enhancing the aggregation and misfolding tendency of mutated protein. The homodimeric state of SOD1 is found to be held together by the non-covalent interactions. The molecular dynamics simulation has been used to study the hydrogen bond interactions between the dimer interface residues of the monomers in native and mutated forms of SOD1 in apo- and holo-states. The results obtained by this analysis reveal the fact that the loss of hydrogen bond interactions between the monomers of the dimer is responsible for the reduced stability of the apo- and holo-mutant forms of SOD1. The conformers with dimer interface residues in native and mutated protein obtained by the molecular dynamics simulation is subjected to quantum mechanical study using M052X/6-31G(d) level of theory. The charge transfer between N-H···O interactions in the dimer interface residues were studied. The weak interaction between the monomers of the dimer accounts for the reduced dimerization and enhanced deformation energy in the mutated SOD1 protein.