Simulation Studies on the Stabilities of Aggregates Formed by Fibril-Forming Segments of α-Synuclein

Abstract

We performed molecular dynamics simulations for various oligomers with different β-sheet conformations consisting of α-Synuclein 71–82 residues using an all atom force field and explicit water model. Tetramers of antiparallel β-sheet are shown to be stable, whereas parallel sheets are highly unstable due to the repulsive interactions between bulky and polar side chains as well as the weaker backbone hydrogen bonds. We also investigated the stabilities of double antiparallel β-sheets stacked with asymmetric and symmetric geometries. Our results show that this 12 amino acid residue peptide can form stable β-sheet conformers at 320K and higher temperatures. The backbone hydrogen bonds in β-sheet and the steric packing between hydrophobic side chains between β-sheets are shown to give conformational stabilities.     

Effects of Turn Stability and Side-Chain Hydrophobicity on the Folding of β-Structures

Key elements of β-structure folding include hydrophobic core collapse, turn formation, and assembly of backbone hydrogen bonds. In the present folding simulations of several β-hairpins and β-sheets (peptide 1, protein G B1 domain peptide, TRPZIP2, TRPZIP4, 20mer, and 20mer^DP6D), the folding free-energy landscape as a function of several reaction coordinates corresponding to the three key elements indicates apparent dependence on turn stability and side-chain hydrophobicity, which demonstrates different folding mechanisms of similar β-structures of varied sequences. Turn stability is found to be the key factor in determining the formation order of the three structural elements in the folding of β-structures. Moreover, turn stability and side-chain hydrophobicity both affect the stability of backbone hydrogen bonds. The three-stranded β-sheets fold through a three-state transition in which the formation of one hairpin always takes precedence over the other. The different stabilities of two anti-parallel hairpins in each three-stranded β-sheet are shown to correlate well with the different levels of their hydrophobic interactions.     

Tryptophan- and arginine-rich antimicrobial peptides: Structures and mechanisms of action

Antimicrobial peptides encompass a number of different classes, including those that are rich in a particular amino acid. An important subset are peptides rich in Arg and Trp residues, such as indolicidin and tritrpticin, that have broad and potent antimicrobial activity. The importance of these two amino acids for antimicrobial activity was highlighted through the screening of a complete combinatorial library of hexapeptides. These residues possess some crucial chemical properties that make them suitable components of antimicrobial peptides. Trp has a distinct preference for the interfacial region of lipid bilayers, while Arg residues endow the peptides with cationic charges and hydrogen bonding properties necessary for interaction with the abundant anionic components of bacterial membranes. In combination, these two residues are capable of participating in cation-π interactions, thereby facilitating enhanced peptide-membrane interactions. Trp sidechains are also implicated in peptide and protein folding in aqueous solution, where they contribute by maintaining native and nonnative hydrophobic contacts. This has been observed for the antimicrobial peptide from human lactoferrin, possibly restraining the peptide structure in a suitable conformation to interact with the bacterial membrane. These unique properties make the Arg- and Trp-rich antimicrobial peptides highly active even at very short peptide lengths. Moreover, they lead to structures for membrane-mimetic bound peptides that go far beyond regular α-helices and β-sheet structures. In this review, the structures of a number of different Trp- and Arg-rich antimicrobial peptides are examined and some of the major mechanistic studies are presented.     

Dimers of human β-defensins and their interactions with the POPG membrane

The stable dimeric structures of human β-defensin (HBD)-3 and -28 have been first computationally identified via a protein docking approach in conjunction with all-atom molecular dynamic simulation. We found that both HBD dimers contain an extended β-sheet platform stabilised mainly by the interaction of second β-sheets and further investigated interaction mechanisms of these dimers including HBD-2 against 1-palmitoyl-2-oleoyl-sn-phosphatidylglycerol membrane bilayer by using coarse-grained model combined with the ElNeDyn network. The extended β-sheet platform of the HBD dimer stayed over the bilayer due to the attachment of the amphipathic region located on one side of the β-sheet platform. The hydrophobic residues of HBDs on the surface interact with the hydrophobic tails of the lipids, whereas the positively charged residues interact with the lipid polar head groups. Finally, antimicrobial nature of HBD-2, HBD-3 and HBD-28 dimers is found to be kept because they are not detached in interacting with the membrane.     

Direct observations of shifts in the β-sheet register of a protein-peptide complex using explicit solvent simulations

Using explicit solvent molecular dynamics simulations, we were able to obtain direct observations of shifts in the hydrogen-bonding register of an intermolecular β-sheet protein-peptide complex. The β-sheet is formed between the FHA domain of cancer marker protein Ki67 (Ki67FHA) and a peptide fragment of the hNIFK signaling protein. Potential encounter complexes of the Ki67FHA receptor and hNIFK peptide are misregistered states of the β-sheet. Rearrangements of one of these misregistered states to the native state were captured in three independent simulations. All three rearrangements occurred by a common mechanism: an aromatic residue of the peptide (F263) anchors into a transient hydrophobic pocket of the receptor to facilitate the formation of native hydrogen bonds. To our knowledge, these simulations provide the first atomically detailed visualizations of a mechanism by which nature might correct for errors in the alignment of intermolecular β-sheets.     

The investigation of the secondary structure propensities and free-energy landscapes of peptide ligands by replica exchange molecular dynamics simulations

The conformational states of two peptide sequences that bind to staphylococcal enterotoxin B are sampled by replica exchange molecular dynamic (REMD) simulations in explicit water. REMD simulations were treated with 52 replicas in the range of 280–501 K for both peptides. The conformational ensembles of both peptides are dominated by random coil, bend and turn structures with a small amount of helical structures for each temperature. In addition, while an insignificant presence of β-bridge structures were observed for both peptides, the β-sheet structure was observed only for peptide 3. The results obtained from simulations at 300 K are consistent with the experimental results obtained from circular dichroism spectroscopy. From the analysis of REMD results, we also calculated hydrophobic and hydrophilic solvent accessible surface areas for both peptides, and it was observed that the hydrophobic segments of the peptides tend to form bend or turn structures. Moreover, the free-energy landscapes of both peptides were obtained by principal component analysis to understand how the secondary structural properties change according to their complex space. From the free-energy analysis, we have found several minima for both peptides at decreased temperature. For these obvious minima of both peptides, it was observed that the random coil, bend and turn structures are still dominant and the helix, β-bridge or β-sheet structures can appear or disappear with respect to minima. On the other hand, when we compare the results of REMD with conventional MD simulations for these peptides, the configurations of peptide 3 might be trapped in energy minima during the conventional MD simulations. Hence, it can be said that the REMD simulations have provided a sufficiently high sampling efficiency.     

Analysis and prediction of the packing of α-helices against a β-sheet in the tertiary structure of globular proteins

The packing of α-helices and β-sheets in six $\text{α}\text{β}$ proteins (e.g. flavodoxin) has been analysed. The results provide the basis for a computer algorithm to predict the tertiary structure of an $\text{α}\text{β}$ protein from its amino acid sequence and actual assignment of secondary structure.The packing of an individual α-helix against a β-sheet generally involves two adjacent ± 4 rows of non-polar residues on the α-helix at the positions i, i + 4, i + 8, i + 1, i + 5, i + 9. The pattern of interacting β-sheet residues results from the twisted nature of the sheet surface and the attendant rotation of the side-chains. At a more detailed level, four of the α-helical residues (i + 1, i + 4, i + 5 and i + 8) form a diamond that surrounds one particular β-sheet residue, generally isoleucine, leucine or valine. In general, the α-helix sits 10 Å above the sheet and lies parallel to the strand direction.The prediction follows a combinational approach. First, a list of possible β-sheet structures (10⁶ to 10¹⁴) is constructed by the generation of all β-sheet topologies and β-strand alignments. This list is reduced by constraints on topology and the location of non-polar residues to mediate the sheet/helix packing, and then rank-ordered on the extent of hydrogen bonding. This algorithm was uniformly applied to 16 $\text{α}\text{β}$ domains in 13 proteins. For every structure, one member of the reduced list was close to the crystal structure; the root-mean-square deviation between equivalenced C^α atoms averaged 5.6 Å for 100 residues. For the $\text{α}\text{β}$ proteins with pure parallel β-sheets, the total number of structures comparable to or better than the native in terms of hydrogen bonds was between 1 and 148. For proteins with mixed β-sheets, the worst case is glyceraldehyde-3-phosphate dehydrogenase, where as many as 3800 structures would have to be sampled. The evolutionary significance of these results as well as the potential use of a combinatorial approach to the protein folding problem are discussed.     

Aggregation of amphipathic peptides at an aqueous–organic interface using coarse-grained simulations

The effect of an aqueous/organic interface on the folding and aggregation of amphipathic peptides is examined by applying discontinuous molecular dynamics (DMD) simulations combined with an intermediate resolution protein model, PRIME20, to a peptide/interface system. The systems contain 48 (KLLK)₄ peptides in random coil or α-helical conformations interacting with both strong and weak interfaces. In the absence of an interface, most of the oligomers form helical bundles, a small fraction of which convert to β-sheets when the temperature is above the folding transition. Adding a weak interface decreases oligomer formation above the folding temperature and increases it below. Little monolayer formation is observed at the weak interface; instead reversible adsorption increases the local peptide concentration near the interface, promoting helical bundle formation in the aqueous phase below the folding temperature and β-sheet formation above the folding temperature. Introducing a strong interface leads to irreversible adsorption, promoting formation of helical monolayers below the folding temperature and mixed β-sheet/amorphous monolayers above the folding temperature. The (KLLK)₄ peptide is more likely to adsorb to the interface when it is in an α-helical conformation, as opposed to a random coil, because of its larger hydrophobic moment.     

Atomic-Scale Simulations Confirm that Soluble β-Sheet-Rich Peptide Self-Assemblies Provide Amyloid Mimics Presenting Similar Conformational Properties

The peptide self-assembly mimic (PSAM) from the outer surface protein A (OspA) can form highly stable but soluble β-rich self-assembly-like structures similar to those formed by native amyloid-forming peptides. However, unlike amyloids that predominantly form insoluble aggregates, PSAMs are highly water-soluble. Here, we characterize the conformations of these soluble β-sheet-rich assemblies. We simulate PSAMs with different-sized β-sheets in the presence and absence of end-capping proteins using all-atom explicit-solvent molecular dynamics, comparing the structural stability, conformational dynamics, and association force. Structural and free-energy comparisons among β-sheets with different numbers of layers and sequences indicate that in similarity to amyloids, the intersheet side chain-side chain interactions and hydrogen bonds combined with intrasheet salt bridges are the major driving forces in stabilizing the overall structural organization. A detailed structural analysis shows that in similarity to amyloid fibrils, all wild-type and mutated PSAM structures display twisted and bent β-sheets to some extent, implying that a twisted and bent β-sheet is a general motif of β-rich assemblies. Thus, our studies indicate that soluble β-sheet-rich peptide self-assemblies can provide good amyloid mimics, and as such confirm on the atomic scale that they are excellent systems for amyloid studies. These results provide further insight into the usefulness of such mimics for nanostructure design.     

Aromatic cross-strand ladders control the structure and stability of beta-rich peptide self-assembly mimics

Though β-rich self-assemblies comprise a major structural class of polypeptides, a detailed understanding of the determinants of their structure and stability is lacking. In particular, the roles of repetitive stretches of side chains running the long axis of these β-sheets, termed “cross-strand ladders,” remain poorly characterized due to the inherently insoluble and heterogeneous nature of self-assemblies. To overcome these experimental challenges, we have established a complementary experimental system termed “peptide self-assembly mimics” (PSAMs). The PSAMs capture a defined number of self-assembly-like peptide repeats within a soluble β-rich protein, making structural and energetic studies possible. In this work, we investigated the role of cross-strand ladders containing aromatic residues, which are prominent in self-assembling peptides. A combination of solution data and high-resolution crystal structures revealed that a single cross-strand ladder consisting solely of Tyr significantly stabilized, rigidified, and flattened the PSAM β-sheet. These characteristics would stabilize each β-sheet layer of a self-assembly and direct sheet conformations compatible with lamination. Our results therefore provide a rationale for the abundance of aromatic amino acids in fibril-forming peptides and establish important roles of cross-strand Tyr ladders in the structure and stability of β-rich peptide self-assemblies.     

Structure and alignment of the membrane-associated antimicrobial peptide arenicin by oriented solid-state NMR spectroscopy

The antimicrobial arenicin peptides are cationic amphipathic sequences that strongly interact with membranes. Through a cystine ring closure a cyclic β-sheet structure is formed in aqueous solution, which persists when interacting with model membranes. In order to investigate the conformation, interactions, dynamics, and topology of their bilayer-associated states, arenicin 1 and 2 were prepared by chemical solid-phase peptide synthesis or by bacterial overexpression, labeled selectively or uniformly with (15)N, reconstituted into oriented membranes, and investigated by proton-decoupled (31)P and (15)N solid-state NMR spectroscopy. Whereas the (31)P NMR spectra indicate that the peptide induces orientational disorder at the level of the phospholipid head groups, the (15)N chemical shift spectra agree well with a regular β-sheet conformation such as the one observed in micellar environments. In contrast, the data do not fit the twisted β-sheet structure found in aqueous buffer. Furthermore, the chemical shift distribution is indicative of considerable conformational and/or topological heterogeneity when at the same time the (15)N NMR spectra exclude alignments of the peptide where the β-sheet lies side ways on the membrane surface. The ensemble of experimental constraints, the amphipathic character of the peptide, and in particular the distribution of the six arginine residues are in agreement with a boatlike dimer structure, similar or related to the one observed in micellar solution, that floats on the membrane surface with the possibility to oligomerize into higher order structures and/or to insert in a transmembrane fashion.     

Structures and dynamics of β-barrel oligomer intermediates of amyloid-beta16-22 aggregation

Accumulating evidence suggests that soluble oligomers are more toxic than final fibrils of amyloid aggregations. Among the mixture of inter-converting intermediates with continuous distribution of sizes and secondary structures, oligomers in the β-barrel conformation – a common class of protein folds with a closed β-sheet – have been postulated as the toxic species with well-defined three-dimensional structures to perform pathological functions. A common mechanism for amyloid toxicity, therefore, implies that all amyloid peptides should be able to form β-barrel oligomers as the aggregation intermediates. Here, we applied all-atom discrete molecular dynamics (DMD) simulations to evaluate the formation of β-barrel oligomers and characterize their structures and dynamics in the aggregation of a seven-residue amyloid peptide, corresponding to the amyloid core of amyloid-β with a sequence of ¹⁶KLVFFAE²² (Aβ16-22). We carried out aggregation simulations with various numbers of peptides to study the size dependence of aggregation dynamics and assembly structures. Consistent with previous computational studies, we observed the formation of β-barrel oligomers in all-atom DMD simulations. Using a network-based approach to automatically identify β-barrel conformations, we systematically characterized β-barrels of various sizes. Our simulations revealed the conformational inter-conversion between β-barrels and double-layer β-sheets due to increased structural strains upon forming a closed β-barrel while maximizing backbone hydrogen bonds. The potential of mean force analysis further characterized the free energy barriers between these two states. The obtained structural and dynamic insights of β-barrel oligomers may help better understand the molecular mechanism of oligomer toxicities and design novel therapeutics targeting the toxic β-barrel oligomers. This article is part of a Special Issue entitled: Protein Aggregation and Misfolding at the Cell Membrane Interface edited by Ayyalusamy Ramamoorthy     

Investigating the effects of N-terminal acetylation on KFE8 self-assembly with 2D IR spectroscopy

Peptide self-assembly is an exciting and robust approach to create novel nanoscale materials for biomedical applications. However, the complex interplay between intra- and intermolecular interactions in peptide aggregation means that minor changes in peptide sequence can yield dramatic changes in supramolecular structure. Here, we use two-dimensional infrared spectroscopy to study a model amphiphilic peptide, KFE8, and its N-terminal acetylated counterpart, AcKFE8. Two-dimensional infrared spectra of isotope-labeled peptides reveal that AcKFE8 aggregates comprise two distinct β-sheet structures although KFE8 aggregates comprise only one of these structures. Using an excitonic Hamiltonian to simulate the vibrational spectra of model β-sheets, we determine that the spectra are consistent with antiparallel β-sheets with different strand alignments, specifically a two-residue shift in the register of the β-strands. These findings bring forth new insights into how N-terminal acetylation may subtly impact secondary structure, leading to larger effects on overall aggregate morphology. In addition, these results highlight the importance of understanding the residue-level structural differences that result from changes in peptide sequence to facilitate the rational design of peptide materials.     

Conformational entropy limits the transition from nucleation to elongation in amyloid aggregation

The formation of β-sheet-rich amyloid fibrils in Alzheimer’s disease and other neurodegenerative disorders is limited by a slow nucleation event. To understand the initial formation of β-sheets from disordered peptides, we used all-atom simulations to parameterize a lattice model that treats each amino acid as a binary variable with β- and non-β-sheet states. We show that translational and conformational entropy give the nascent β-sheet an anisotropic surface tension that can be used to describe the nucleus with 2D classical nucleation theory. Since translational entropy depends on concentration, the aspect ratio of the critical β-sheet changes with protein concentration. Our model explains the transition from the nucleation phase to elongation as the point where the β-sheet core becomes large enough to overcome the conformational entropy cost to straighten the terminal molecule. At this point the β-strands in the nucleus spontaneously elongate, which results in a larger binding surface to capture new molecules. These results suggest that nucleation is relatively insensitive to sequence differences in coaggregation experiments because the nucleus only involves a small portion of the peptide.     

Effects of arginine density on the membrane-bound structure of a cationic antimicrobial peptide from solid-state NMR

Solid-state NMR spectroscopy is used to determine the membrane-bound topological structure of a cationic β-hairpin antimicrobial peptide in which the number of Arg residues has been halved. The parent peptide, PG-1, was previously found to form transmembrane β-barrels in anionic membranes where the Arg residues complex with the lipid phosphate groups to cause toroidal pore defects in the membrane. In comparison, the charge-attenuated and less active mutant studied here forms β-sheets that lie on the surface of the zwitterionic membrane and only partially insert into the anionic membrane. The mutant also exhibits much looser contact with the lipid headgroups. These results indicate that transmembrane insertion and tight Arg-phosphate association are two important elements for strong antimicrobial activities of this class of peptides. Comparison with other β-hairpin antimicrobial peptides studied so far further suggests a relative potency scale for the various mechanisms of action for the β-sheet family of antimicrobial peptides. The transmembrane insertion-toroidal pore mechanism is the most potent in disrupting the lipid bilayer, followed by the large-amplitude in-plane motional mechanism. The carpet model, where peptides aggregate on the membrane surface to cause lateral expansion and eventual micellization of the membrane, is a weaker mechanism of action.     

Designed β-sheet-forming peptide 33mers with potent human bactericidal/permeability increasing protein-like bactericidal and endotoxin neutralizing activities

Novel peptide 33mers have been designed by incorporating β-conformation stabilizing residues from the β-sheet domains of α-chemokines and functionally important residues from the β-sheet domain of human neutrophil bactericidal protein (B/PI). B/PI is known for its ability to kill bacteria and to neutralize the action of bacterial endotoxin (lipopolysaccharide, LPS) which can induce septic shock leading to eventual death. Here, the goal was to make short linear peptides which demonstrate good β-sheet folding and maintain bioactivity as in native B/PI. A library of 24 peptide 33mers (βpep-1 to βpep-24) were synthesized with various amino acid substitutions. CD and NMR data acquired in aqueous solution indicate that βpep peptides form β-sheet structure to varying degrees and self-associate as dimers and tetramers like the α-chemokines. Bactericidal activity toward Gram-negative Pseudomonas aeruginosa was tested, and βpep-19 was found to be only about 5-fold less potent (62% kill at 1.2×10^?7 M) than native B/PI (80% kill at 2.9×10^?8 M). At LPS neutralization, βpep-2 and -23 were found to be most active (66–78% effective at 1.2×10^?6 M), being only about 50–100-fold less active than B/PI (50% at 1.5×10^?8 M). In terms of structure–activity relations, β-sheet structural stability correlates with the capacity to neutralize LPS, but not with bactericidal activity. Although a net positive charge is necessary for activity, it is not sufficient for optimal activity. Hydrophobic residues tend to influence activities indirectly by affecting structural stability. Furthermore, results show that sequentially and spatially related residues from the β-sheet domain of native B/PI can be designed into short linear peptides which show good β-sheet folding and retain much of the native activity. This research contributes to the development of solutions to the problem of multiple drug-resistant, opportunistic microorganisms like P. aeruginosa and of agents effective at neutralizing bacterial endotoxin.     

Homophilic adhesion mechanism of neurofascin, a member of the L1 family of neural cell adhesion molecules

The L1 family neural cell adhesion molecules play key roles in specifying the formation and remodeling of the neural network, but their homophilic interaction that mediates adhesion is not well understood. We report two crystal structures of a dimeric form of the headpiece of neurofascin, an L1 family member. The four N-terminal Ig-like domains of neurofascin form a horseshoe shape, akin to several other immunoglobulin superfamily cell adhesion molecules such as hemolin, axonin, and Dscam. The neurofascin dimer, captured in two crystal forms with independent packing patterns, reveals a pair of horseshoes in trans-synaptic adhesion mode. The adhesion interaction is mediated mostly by the second Ig-like domain, which features an intermolecular β-sheet formed by the joining of two individual GFC β-sheets and a large but loosely packed hydrophobic cluster. Mutagenesis combined with gel filtration assays suggested that the side chain hydrogen bonds at the intermolecular β-sheet are essential for the homophilic interaction and that the residues at the hydrophobic cluster play supplementary roles. Our structures reveal a conserved homophilic adhesion mode for the L1 family and also shed light on how the pathological mutations of L1 affect its structure and function.     

Secondary structure of proteins associated in thin films

The solid state secondary structure of myoglobin, RNase A, concanavalin A (Con A), poly(L -lysine), and two linear heterooligomeric peptides were examined by both far-uv CD spectroscopy¹ and by ir spectroscopy. The proteins associated from water solution on glass and mica surfaces into noncrystalline, amorphous films, as judged by transmission electron microscopy of carbon-platinum replicas of surface and cross-fractured layer. The association into the solid state induced insignificant changes in the amide CD spectra of all α-helical myoglobin, decreased the molar ellipticity of the α/β RNase A, and increased the molar ellipticity of all-β Con A with no change in the positions of the bands' maxima. High-temperature exposure of the films induced permanent changes in the conformation of all proteins, resulting in less α-helix and more β-sheet structure. The results suggest that the protein α-helices are less stable in films and that the secondary structure may rearrange into β-sheets at high temperature. Two heterooligomeric peptides and poly (L -lysine), all in solution at neutral pH with “random coil” conformation, formed films with variable degrees of their secondary structure in β-sheets or β-turns. The result corresponded to the protein-derived Chou-Fasman amino acid propensities, and depended on both temperature and solvent used. The ir and CD spectra correlations of the peptides in the solid state indicate that the CD spectrum of a “random” structure in films differs from random coil in solution. Formic acid treatment transformed the secondary structure of the protein and peptide films into a stable α-helix or β-sheet conformations. The results indicate that the proteins aggregate into a noncrystalline, glass-like state with preserved secondary structure. The solid state secondary structure may undergo further irreversible transformations induced by heat or solvent. © 1993 John Wiley & Sons, Inc.     

Secondary structure propensity and chirality of the amyloidophilic peptide p5 and its analogues impacts ligand binding - In vitro characterization

BackgroundPolybasic helical peptides, such as peptide p5, bind human amyloid extracts and synthetic amyloid fibrils. When radiolabeled, peptide p5 has been shown to specifically bind amyloid in vivo thereby allowing imaging of the disease. Structural requirements for heparin and amyloid binding have been studied using analogues of p5 that modify helicity and chirality.MethodsPeptide-ligand interactions were studied using CD spectroscopy and solution-phase binding assays with radiolabeled p5 analogues. The interaction of a subset of peptides was further studied by using molecular dynamics simulations.ResultsDisruption of the peptide helical structure reduced peptide binding to heparin and human amyloid extracts. The all-D enantiomer and the β-sheet-structured peptide bound all substrates as well as, or better than, p5. The interaction of helical and β-sheet structured peptides with Aβ fibrils was modeled and shown to involve both ionic and non-ionic interactions.ConclusionsThe α-helical secondary structure of peptide p5 is important for heparin and amyloid binding; however, helicity is not an absolute requirement as evidenced by the superior reactivity of a β-sheet peptide. The differential binding of the peptides with heparin and amyloid fibrils suggests that these molecular interactions are different. The all-D enantiomer of p5 and the β-sheet peptide are candidates for amyloid targeting reagents in vivo.General SignificanceEfficient binding of polybasic peptides with amyloid is dependent on the linearity of charge spacing in the context of an α-helical secondary structure. Peptides with an α-helix or β-sheet propensity and with similar alignment of basic residues is optimal.