Wild type and mutants of the HET‐s(218–289) prion show different flexibility at fibrillar ends: A simulation study

The C‐terminal segment (residues 218–289) of the HET‐s protein of the filamentous fungus Podosporina anserina is a prion‐forming domain. The structural model of the HET‐s(218–289) amyloid fibril based on solid‐state nuclear magnetic resonance (NMR) restraints shows a β solenoid topology which is comprised of a β‐sheet core and interconnecting loops. For the single‐point mutants Phe286Ala and Trp287Ala, slower aggregation rates in vitro and loss of prionic infectivity have been reported recently. Here we have used molecular dynamics to compare the flexibility of the mutants and wild type. The simulations, initiated from a trimeric aggregate extracted from the NMR structural model, show structural stability on a 100‐ns time scale for wild type and mutants. Analysis of the fluctuations along the simulations reveals that the mutants are less flexible than the wild type in the C‐terminal segment at only one of the two external monomers. Analysis of interaction energy and buried accessible surface indicates that residue Phe286 in particular is stabilized in the Trp287Ala mutant. The simulation results provide an atomistic explanation of the suggestion (based on indirect experimental evidence) that flexibility at the protofibril end(s) is required for fibril elongation. Moreover, they provide further evidence that the growth of the HET‐s amyloid fibril is directional. Proteins 2014; 82:399–404. © 2013 Wiley Periodicals, Inc.     

The N-terminal propeptide of lung surfactant protein C is necessary for biosynthesis and prevents unfolding of a metastable alpha-helix

The lung surfactant-associated protein C (SP-C) consists mainly of a polyvaline alpha-helix, which is stable in a lipid membrane. However, in agreement with the predicted beta-strand conformation of a polyvaline segment, helical SP-C unfolds and transforms into beta-sheet aggregates and amyloid fibrils within a few days in aqueous organic solvents. SP-C fibril formation and aggregation have been associated with lung disease. Here, we show that in a recently isolated biosynthetic precursor of SP-C (SP-Ci), a 12 residue N-terminal propeptide locks the metastable polyvaline part in a helical conformation. The SP-Ci helix does not aggregate or unfold during several weeks of incubation, as judged by hydrogen/deuterium exchange and mass spectrometry. Hydrogen/deuterium exchange experiments further indicate that the propeptide reduces exchange in parts corresponding to mature SP-C. Finally, in an acidic environment, SP-Ci unfolds and aggregates into amyloid fibrils like SP-C. These data suggest a direct role of the N-terminal propeptide in SP-C biosynthesis.     

Molecular dynamics simulation of the SH3 domain aggregation suggests a generic amyloidogenesis mechanism

We use molecular dynamics simulation to study the aggregation of Src SH3 domain proteins. For the case of two proteins, we observe two possible aggregation conformations: the closed form dimer and the open aggregation state. The closed dimer is formed by "domain swapping"-the two proteins exchange their RT-loops. All the hydrophobic residues are buried inside the dimer so proteins cannot further aggregate into elongated amyloid fibrils. We find that the open structure-stabilized by backbone hydrogen bond interactions-packs the RT-loops together by swapping the two strands of the RT-loop. The packed RT-loops form a beta-sheet structure and expose the backbone to promote further aggregation. We also simulate more than two proteins, and find that the aggregate adopts a fibrillar double beta-sheet structure, which is formed by packing the RT-loops from different proteins. Our simulations are consistent with a possible generic amyloidogenesis scenario.     

Molecular dynamics simulation of the nanofibrils formed by amyloid-based peptide amphiphiles

Abstract

Atomistic molecular dynamics simulations have been performed on the peptide amphiphiles (PAs) with four amyloid beta peptide fragments as head groups. The stable structures were monitored by the root mean square deviation with respect to the energy minimised initial structures. Random coil and β-sheet structures with hydrogen bonds along and perpendicular to the long axis of the nanofibre were obtained due to the different nature of the head groups. Influences of pH and capping ends on the nanofibre structures were investigated through variation of the protonation states of the ionic amino acids in the peptides. The peptides with opposite charges on both sides were found to have the fewest β-sheet structures, and the charges on the outer terminal tended to destruct the β-sheets while those at the inner side did not. The isolated charge in the centre of peptides was found to be able to promote the formation of regular β-sheets, while multiple charged residues could not support ordered β-sheet structures. When charge neutralisation occurred between adjacent residues, regular β-sheet laminates might also occur for systems with charges at the outer terminal. With the increase of β-sheet structures formed, the original twisted structures found for random coil structures of the PAs could be diminished by the hydrogen bonds.     

Pathways and intermediates of amyloid fibril formation

The lack of understanding of amyloid fibril formation at the molecular level is a major obstacle in devising strategies to interfere with the pathologies linked to peptide or protein aggregation. In particular, little is known on the role of intermediates and fibril elongation pathways as well as their dependence on the intrinsic tendency of a polypeptide chain to self-assembly by β-sheet formation (β-aggregation propensity). Here, coarse-grained simulations of an amphipathic polypeptide show that a decrease in the β-aggregation propensity results in a larger heterogeneity of elongation pathways, despite the essentially identical structure of the final fibril. Protofibrillar intermediates that are thinner, shorter and less structured than the final fibril accumulate along some of these pathways. Moreover, the templated formation of an additional protofilament on the lateral surface of a protofibril is sometimes observed as a collective transition. Conversely, for a polypeptide model with a high β-aggregation propensity, elongation proceeds without protofibrillar intermediates. Therefore, changes in intrinsic β-aggregation propensity modulate the relative accessibility of parallel routes of aggregation.     

Atomistic mechanism of polyphenol amyloid aggregation inhibitors: molecular dynamics study of Curcumin,Exifone, and Myricetin interaction with the segment of tau peptide oligomer

Amyloid fibrils are highly ordered protein aggregates associated with many diseases affecting millions of people worldwide. Polyphenols such as Curcumin, Exifone, and Myricetin exhibit modest inhibition toward fibril formation of tau peptide which is associated with Alzheimer’s disease. However, the molecular mechanisms of this inhibition remain elusive. We investigated the binding of three polyphenol molecules to the protofibrils of an amyloidogenic fragment VQIVYK of tau peptide by molecular dynamics simulations in explicit solvent. We find that polyphenols induce conformational changes in the oligomer aggregate. These changes disrupt the amyloid H bonding, perturbing the aggregate. While the structural evolution of the control oligomer with no ligand is limited to the twisting of the β-sheets without their disassembly, the presence of polyphenol molecule pushes the β-sheets apart, and leads to a loosely packed structure where two of four β-sheets dissociate in each of the three cases considered here. The H-bonding capacity of polyphenols is responsible for the observed behavior. The calculated binding free energies and its individual components enabled better understanding of the binding. Results indicated that the contribution from Van der Waals interactions is more significant than electrostatic contribution to the binding. The findings from this study are expected to assist in the development of aggregation inhibitors. Significant binding between polyphenols and aggregate oligomer identified in our simulations confirms the previous experimental observations in which polyphenols refold the tau peptide without forming covalent bonds.     

Folding and dimerization of the ionic peptide EAK 16-IV

Folding and dimerization of an ionic polyalanine-based peptide chain (EAK16-IV) are simulated with nonspecific interactions. It is found that there is a competition between two kinds of structural motifs under different strengths of electrostatic interactions. The dominance of hairpin-like structures would be realized with a strong electrostatic interaction both thermodynamically and kinetically, showing the importance of the electrostatic interaction on the formation of hairpin-like structures. Simulations on the dimerization with strong electrostatic interaction are also carried out. It is found that the concentration contributes essentially to the shape of the dimers. These studies demonstrate that the strong interactions and kinetic factors might be important for the ordered amyloid aggregates.     

Heme binding inhibits the fibrillization of amyloidogenic apomyoglobin and determines lack of aggregate cytotoxicity

Myoglobin is an alpha-helical globular protein containing two highly conserved tryptophanyl residues at positions 7 and 14 in the N-terminal region. The double W/F replacement renders apomyoglobin highly susceptible to aggregation and amyloid-like fibril formation under physiological conditions. In this work we analyze the early stage of W7FW14F apomyoglobin aggregation following the time dependence of the process by far-UV CD, Fourier-transform infrared (FTIR) spectroscopy, and heme-binding properties. The results show that the aggregation of W7FW14F apomyoglobin starts from a native-like globin state able to bind the prosthetic group with spectroscopic properties similar to those observed for wild-type apoprotein. Nevertheless, it rapidly aggregates, forming amyloid fibrils. However, when the prosthetic group is added before the beginning of aggregation, amyloid fibrillization is inhibited, although the aggregation process is not prevented. Moreover, the apomyoglobin aggregates formed in these conditions are not cytotoxic differently from what is observed for all amyloidogenic proteins. These results open new insights into the relationship between the structure adopted by the protein into the aggregates and their ability to trigger the impairment of cell viability.     

Chain-length dependence of alpha-helix to beta-sheet transition in polylysine: model of protein aggregation studied by temperature-tuned FTIR spectroscopy

The chain-length dependence of the alpha-helix to beta-sheet transition in poly(L-lysine) is studied by temperature-tuned FTIR spectroscopy. This study shows that heterogeneous samples of poly(L-lysine), comprising polypeptide chains with various lengths, undergo the alpha-beta transition at an intermediate temperature compared to homogeneous ingredients. This holds true as long as each individual fraction of the polypeptide is capable of adopting an antiparallel beta-sheet structure. The tendency is that the longer chain is, the lower the alpha-beta transition temperature is, which has been linked to the presence of distorted or solvated helices with turns or beta sheets in elongating chains of poly(L-lysine). As such helical structures are apparently conducive to the alpha-beta transition, this draws a comparison to the hypothesis of metastable protein conformational states being a common stage in amyloid-formation pathways. The antiparallel architecture of the beta sheet is likely to reflect the pretransition interhelical interactions in poly(L-lysine). Namely, the chains are arranged in an antiparallel manner because of energetically favored antiparallel pre-assembly of dipolar alpha helices.     

Molecular dynamics simulations of Alzheimer's beta-amyloid protofilaments

Filamentous amyloid aggregates are central to the pathology of Alzheimer's disease. We use all-atom molecular dynamics (MD) simulations with explicit solvent and multiple force fields to probe the structural stability and the conformational dynamics of several models of Alzheimer's beta-amyloid fibril structures, for both wild-type and mutated amino acid sequences. The structural models are based on recent solid state NMR data. In these models, the peptides form in-register parallel beta-sheets along the fibril axis, with dimers of two U-shaped peptides located in layers normal to the fibril axis. Four different topologies are explored for stacking the beta-strand regions against each other to form a hydrophobic core. Our MD results suggest that all four NMR-based models are structurally stable, and we find good agreement with dihedral angles estimated from solid-state NMR experiments. Asp23 and Lys28 form buried salt-bridges, resulting in an alternating arrangement of the negatively and positively charged residues along the fibril axis that is reminiscent of a one-dimensional ionic crystal. Interior water molecules are solvating the buried salt-bridges. Based on data from NMR measurements and MD simulations of short amyloid fibrils, we constructed structural models of long fibrils. Calculated X-ray fiber diffraction patterns show the characteristics of packed beta-sheets seen in experiments, and suggest new experiments that could discriminate between various fibril topologies.     

Interpreting the aggregation kinetics of amyloid peptides

Amyloid fibrils are insoluble mainly beta-sheet aggregates of proteins or peptides. The multi-step process of amyloid aggregation is one of the major research topics in structural biology and biophysics because of its relevance in protein misfolding diseases like Alzheimer's, Parkinson's, Creutzfeld-Jacob's, and type II diabetes. Yet, the detailed mechanism of oligomer formation and the influence of protein stability on the aggregation kinetics are still matters of debate. Here a coarse-grained model of an amphipathic polypeptide, characterized by a free energy profile with distinct amyloid-competent (i.e. beta-prone) and amyloid-protected states, is used to investigate the kinetics of aggregation and the pathways of fibril formation. The simulation results suggest that by simply increasing the relative stability of the beta-prone state of the polypeptide, disordered aggregation changes into fibrillogenesis with the presence of oligomeric on-pathway intermediates, and finally without intermediates in the case of a very stable beta-prone state. The minimal-size aggregate able to form a fibril is generated by collisions of oligomers or monomers for polypeptides with unstable or stable beta-prone state, respectively. The simulation results provide a basis for understanding the wide range of amyloid-aggregation mechanisms observed in peptides and proteins. Moreover, they allow us to interpret at a molecular level the much faster kinetics of assembly of a recently discovered functional amyloid with respect to the very slow pathological aggregation.     

Amyloid-like fibril formation of co-chaperonin GroES: nucleation and extension prefer different degrees of molecular compactness

The molecular chaperone GroES, together with GroEL from Escherichia coli, is the best characterized protein of the molecular chaperone family. Here, we report on the in vitro formation of GroES amyloid-like fibrils and the mechanism of formation. When incubated for several weeks at neutral pH in the presence of the denaturant guanidine hydrochloride, GroES formed a typical amyloid fibril; unbranched, twisted, and extended filaments stainable by thioflavin T and Congo red. GroES fibril formation was accelerated by the addition of preformed fibril seeds, in accordance with a nucleation-extension mechanism. Interestingly, whereas the spontaneous formation of GroES fibrils was favored in the structural transition region of GroES dissociation/unfolding, the extension of fibrils from preformed fibril seeds was favored in the region corresponding to an expanded molecular state. We concluded that the two stages of GroES fibril formation prefer different molecular states of the same protein. The significance of this preference is discussed.     

Interpeptide interactions induce helix to strand structural transition in Aβ peptides

Replica exchange molecular dynamics and all‐atom implicit solvent model are used to compute the structural propensities in Aβ monomers, dimers, and Aβ peptides bound to the edge of amyloid fibril. These systems represent, on an approximate level, different stages in Aβ aggregation. Aβ monomers are shown to form helical structure in the N‐terminal (residues 13 to 21). Interpeptide interactions in Aβ dimers and, especially, in the peptides bound to the fibril induce a dramatic shift in the secondary structure, from helical states toward β‐strand conformations. The sequence region 10–23 in Aβ peptide is found to form most of interpeptide interactions upon aggregation. Simulation results are tested by comparing the chemical shifts in Aβ monomers computed from simulations and obtained experimentally. Possible implications of our simulations for designing aggregation‐resistant variants of Aβ are discussed. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Amyloidogenic sequences in native protein structures

Numerous short peptides have been shown to form β‐sheet amyloid aggregates in vitro. Proteins that contain such sequences are likely to be problematic for a cell, due to their potential to aggregate into toxic structures. We investigated the structures of 30 proteins containing 45 sequences known to form amyloid, to see how the proteins cope with the presence of these potentially toxic sequences, studying secondary structure, hydrogen‐bonding, solvent accessible surface area and hydrophobicity. We identified two mechanisms by which proteins avoid aggregation: Firstly, amyloidogenic sequences are often found within helices, despite their inherent preference to form β structure. Helices may offer a selective advantage, since in order to form amyloid the sequence will presumably have to first unfold and then refold into a β structure. Secondly, amyloidogenic sequences that are found in β structure are usually buried within the protein. Surface exposed amyloidogenic sequences are not tolerated in strands, presumably because they lead to protein aggregation via assembly of the amyloidogenic regions. The use of α‐helices, where amyloidogenic sequences are forced into helix, despite their intrinsic preference for β structure, is thus a widespread mechanism to avoid protein aggregation.     

Stability of Aβ‐fibril fragments in the presence of fatty acids

We consider the effect of lauric acid on the stability of various fibril‐like assemblies of Aβ peptides. For this purpose, we have performed molecular dynamics simulations of these assemblies either in complex with lauric acid or without presence of the ligand. While we do not observe a stabilizing effect on Aβ₄₀‐fibrils, we find that addition of lauric acid strengthens the stability of fibrils built from the triple‐stranded S‐shaped Aβ₄₂‐peptides considered to be more toxic. Or results may help to understand how the specifics of the brain‐environment modulate amyloid formation and propagation.     

Effects of interfaces on aggregates of peptides derived from pancreatic islet amyloid polypeptide

Recent experimental studies indicate that oligomeric complexes of misfolded proteins and peptides are the primary agents of cytotoxicity in amyloid-related diseases. Given the prevalence of mixed-polarity interfaces in physiological environments, an understanding of the mechanisms of interactions between amorphous (pre-fibrillar) aggregates and surfaces is important for completing our knowledge of the behaviour of peptide aggregation phenomena. We have employed fully solvated molecular dynamics simulations to study the morphology, interactions and peptide conformations of disordered aggregates of the amyloidogenic NFGAIL (derived from human islet amyloid polypeptide) and non-amyloidogenic AGAIL peptides upon adsorption to vapour–water, decane–water, bilayer and solid–water interfaces. All of the interfaces studied promote elongation and surface-spreading of both peptide aggregates, with the liquid–liquid interface being particularly efficient at altering the gross morphology of disordered aggregates. NFGAIL aggregates are more effective at disrupting lipid bilayers compared to AGAIL. Additionally, the interfaces studied cause greater changes in peptide conformations within the NFGAIL aggregates compared to AGAIL. We propose that simulations may elucidate the capability of interfaces to effect changes in the behaviour of disordered peptide aggregates, which may also serve to provide measures of the intrinsic fibrillogenicity of a given peptide sequence.     

Computational studies of the structure, dynamics and native content of amyloid-like fibrils of ribonuclease A

It is a common belief that some residues of a protein are more important than others. In some cases, point mutations of some residues make butterfly effect on the protein structure and function, but in other cases they do not. In addition, the residues important for the protein function tend to be not only conserved but also coevolved with other interacting residues in a protein. Motivated by these observations, the authors propose that there is a network composed of the residues, the residue-residue coevolution network (RRCN), where nodes are residues and links are set when the coevolutionary interaction strengths between residues are sufficiently large. The authors build the RRCN for the 44 diverse protein families. The interaction strengths are calculated by using McBASC algorithm. After constructing the RRCN, the authors identify residues that have high degree of connectivity (hub nodes), and residues that play a central role in network flow of information (C(I) nodes). The authors show that these residues are likely to be functionally important residues. Moreover, the C(I) nodes appear to be more relevant to the function than the hub nodes. Unlike other similar methods, the method described in this study is solely based on sequences. Therefore, the method can be applied to the function annotation of a wider range of proteins.     

Effects of chain length on the aggregation of model polyglutamine peptides: molecular dynamics simulations

Aggregation in the brain of polyglutamine-containing proteins is either a cause or an associated symptom of nine hereditary neurodegenerative disorders including Huntington's disease. The molecular level mechanisms by which these proteins aggregate are still unclear. In an effort to shed light on this important phenomenon, we are investigating the aggregation of model polyglutamine peptides using molecular-level computer simulation with a simplified model of polyglutamine that we have developed. This model accounts for the most important types of intra- and inter-molecular interactions-hydrogen bonding and hydrophobic interactions-while allowing the folding process to be simulated in a reasonable time frame. The model is used to examine the folding of isolated polyglutamine peptides 16, 32, and 48 residues long and the folding and aggregation of systems of 24 model polyglutamine peptides 16, 24, 32, 36, 40, and 48 residues long. Although the isolated polyglutamine peptides did form some alpha and beta backbone-backbone hydrogen bonds they did not have as many of these bonds as they would have if they had folded into a complete alpha helix or beta sheet. In one of the simulations on the isolated polyglutamine peptide 48 residues long, we observed a structure that resembles a beta helix. In the multi-chain simulations we observed amorphous aggregates at low temperatures, ordered aggregates with significant beta sheet character at intermediate temperatures, and random coils at high temperatures. We have found that the temperature at which the model peptides undergo the transition from amorphous aggregates to ordered aggregates and the temperature at which the model peptides undergo the transition from ordered aggregates to random coils increase with increasing chain length. Our finding that the stability of the ordered aggregates increases as the peptide chain length increases may help to explain the experimentally observed relation between polyglutamine tract length and aggregation in vitro and disease progression in vivo. We have also observed in our simulations that the optimal temperature for the formation of beta sheets increases with chain length up to 36 glutamine residues but not beyond. Equivalently, at fixed temperature we find a transition from a region dominated by random coils at chain lengths less than 36 to a region dominated by relatively ordered beta sheet structures at chain lengths greater than 36. Our finding of this critical chain length of 36 glutamine residues is interesting because a critical chain length of 37 glutamine residues has been observed experimentally.