Factors That Affect the Degree of Twist in β-Sheet Structures: A Molecular Dynamics Simulation Study of a Cross-β Filament of the GNNQQNY Peptide

Molecular Insights into the Dynamics of Amyloid Fibril Growth: Elongation and Lateral Assembly of GNNQQNY Protofibrils

The self-assembly of peptides and proteins into β-sheet rich amyloid fibrils is linked to both functional and pathological states. In this study, the growth of fibrillar structures of the short peptide GNNQQNY, a fragment from the yeast prion Sup35 protein, was examined. Molecular dynamics simulations were used to study alternative mechanisms of fibril growth, including elongation through binding of monomers as well as fibril self-assembly into larger, more mature structures. It was found that after binding, monomers diffused along preformed fibrils toward the ends, supporting the mechanism of fibril growth via elongation. Lateral assembly of protofibrils was found to occur readily, suggesting that this could be the key to transitioning from isolated fibrils to mature multilayer structures. Overall, the work provides mechanistic insights into the competitive pathways that govern amyloid fibril growth.

Programming co-assembled peptide nanofiber morphology via anionic amino acid type: Insights from molecular dynamics simulations

Co-assembling peptides can be crafted into supramolecular biomaterials for use in biotechnological applications, such as cell culture scaffolds, drug delivery, biosensors, and tissue engineering. Peptide co-assembly refers to the spontaneous organization of two different peptides into a supramolecular architecture. Here we use molecular dynamics simulations to quantify the effect of anionic amino acid type on co-assembly dynamics and nanofiber structure in binary CATCH(+/-) peptide systems. CATCH peptide sequences follow a general pattern: CQCFCFCFCQC, where all C's are either a positively charged or a negatively charged amino acid. Specifically, we investigate the effect of substituting aspartic acid residues for the glutamic acid residues in the established CATCH(6E-) molecule, while keeping CATCH(6K+) unchanged. Our results show that structures consisting of CATCH(6K+) and CATCH(6D-) form flatter β-sheets, have stronger interactions between charged residues on opposing β-sheet faces, and have slower co-assembly kinetics than structures consisting of CATCH(6K+) and CATCH(6E-). Knowledge of the effect of sidechain type on assembly dynamics and fibrillar structure can help guide the development of advanced biomaterials and grant insight into sequence-to-structure relationships.

Simulations of Amyloid-Forming Peptides in the Crystal State

There still is little treatment available for amyloid diseases, despite their significant impact on individuals and the social and economic implications for society. One reason for this is that the physical nature of amyloid formation is not understood sufficiently well. Therefore, fundamental research at the molecular level remains necessary to support the development of therapeutics. A few structures of short peptides from amyloid-forming proteins have been determined. These can in principle be used as scaffolds for designing aggregation inhibitors. Attempts to this end have often used the tools of computational chemistry, in particular molecular simulation. However, few simulation studies of these peptides in the crystal state have been presented so far. Hence, to validate the capability of common force fields (AMBER19SB, CHARMM36m, and OPLS-AA/M) to yield insight into the dynamics and structural stability of amyloid peptide aggregates, we have performed molecular dynamics simulations of twelve different peptide crystals at two different temperatures. From the simulations, we evaluate the hydrogen bonding patterns, the isotropic B-factors, the change in energy, the Ramachandran plots, and the unit cell parameters and compare the results with the crystal structures. Most crystals are stable in the simulations but for all force fields there is at least one that deviates from the experimental crystal, suggesting more work is needed on these models.

Functionalized Mesoporous Silicas Direct Structural Polymorphism of Amyloid-β Fibrils

The aggregation of amyloid-β (Aβ) is associated with the onset of Alzheimer's disease (AD) and involves a complex kinetic pathway as monomers self-assemble into fibrils. A central feature of amyloid fibrils is the existence of multiple structural polymorphs, which complicates the development of disease-relevant structure-function relationships. Developing these relationships requires new methods to control fibril structure. In this work, we evaluated the effect that mesoporous silicas (SBA-15) functionalized with hydrophobic (SBA-PFDTS) and hydrophilic groups (SBA-PEG) have on the aggregation kinetics and resulting structure of Aβ_1-40 fibrils. The hydrophilic SBA-PEG had little effect on amyloid kinetics, while as-synthesized and hydrophobic SBA-PFDTS accelerated aggregation kinetics. Subsequently, we quantified the relative population of fibril structures formed in the presence of each material using electron microscopy. Fibrils formed from Aβ_1-40 exposed to SBA-PEG were structurally similar to control fibrils. In contrast, Aβ_1-40 incubated with SBA-15 or SBA-PFDTS formed fibrils with shorter crossover distances that were more structurally representative of fibrils found in AD patient derived samples. Overall, our results suggest that mesoporous silicas and other exogenous materials are promising scaffolds for the de novo production of specific fibril polymorphs of Aβ_1-40 and other amyloidogenic proteins.

Growth, Polymorphism, and Spatially Controlled Surface Immobilization of Biotinylated Variants of IAPP21-27 Fibrils

The islet amyloid polypeptide (IAPP) is a regulatory peptide that can aggregate into fibrillar structures associated with type 2 diabetes. In this study, the IAPP_21-27 segment was modified with a biotin linker at the N-terminus (Btn-GNNFGAIL) to immobilize peptide fibrils on streptavidin-coated surfaces. Key residues for fibril formation of the N-terminal biotinylated IAPP_21-27 segment were identified by using an alanine scanning approach combined with molecular dynamics simulations, thioflavin T fluorescence measurements, and scanning electron microscopy. Significant contributions of phenylalanine (F23), leucine (L27), and isoleucine (I26) for the fibrillation of the short peptide segment were identified. The fibril morphologies of the peptide variants differed depending on their primary sequence, ranging from flexible and semiflexible to stiff and crystal-like structures. These insights could advance the design of new functional hybrid bionanomaterials and fibril-engineered surface coatings using short peptide segments. To validate this concept, the biotinylated fibrils were immobilized on streptavidin-coated surfaces under spatial control.

Twisted Ribbon Aggregates in a Model Peptide System

The model peptides A₈K and A₁₀K self-assemble in water into ca. 100 nm long ribbon-like aggregates. These structures can be described as β-sheets laminated into a ribbon structure with a constant elliptical cross-section of 4 by 8 nm, where the longer axis corresponds to a finite number, N ≈ 15, of laminated sheets, and 4 nm corresponds to a stretched peptide length. The ribbon cross-section is strikingly constant and independent of the peptide concentration. High-contrast transmission electron microscopy shows that the ribbons are twisted with a pitch λ ≈ 15 nm. The self-assembly is analyzed within a simple model taking into account the interfacial free energy of the hydrophobic β-sheets and a free energy penalty arising from an increased stretching of hydrogen bonds within the laminated β-sheets, arising from the twist of the ribbons. The model predicts an optimal value N, in agreement with the experimental observations.

Amino acid conformations control the morphological and chiral features of the self-assembled peptide nanostructures: Young investigators perspective

HYPOTHESIS

A variety of nanostructures with different chiral features can be self-assembled from short peptides with highly similar sequences. We hypothesize that these supramolecular nanostructures are ruled by the constituent amino acid residues which adopt their conformations under the influence of intra-/inter-molecular interactions during peptide self-assembly.

APPROACH

Through reviewing recent advances in the self-assembly of short peptides and focusing on the relationship between amino acid conformations, peptide secondary structures and intra-/inter-molecular interactions within the supramolecular architectures, we aim to rationalize the complex interactive processes involved in the self-assembly of short, designed peptides.

RESULTS

Given the highly complexing interactive processes, the adoption of amino acid conformations and their control over peptide self-assembly consist of 4 main steps: (1) Each amino acid residue adopts its unique conformation in a specific sequence; (2) The sequence exhibits its own main chain geometry and determines the propensity of the intermolecular alignment within the building block; (3) The structural propensity of the building block and the packing mode between them determine the self-assembled structural features such as twisting, growth and chirality; (4) In addition to intra-/inter-molecular interactions, inter-sheet and inter-building block interactions could also affect the residue conformations and nanostructures, causing structural readjustment.

Carbamylation promotes amyloidogenesis and induces structural changes in Tau-core hexapeptide fibrils

BACKGROUND

Carbamylation is a non-enzymatic post-translational modification (PTM), which involves the covalent modification of N-terminus of protein or ε-amino group of Lys. The role of carbamylation in several age-related disorders is well documented, however, the relationship between carbamylation and neurodegenerative disorders including Alzheimer's disease remains uncharted.

METHODS

In the present study, using aggregation-prone tau-core hexapeptide fragments ³⁰⁶VQIVYK³¹¹ (PHF6) and ²⁷⁵VQIINK²⁸⁰ (PHF6*) as models, we have elucidated the effect of carbamylation on aggregation kinetics and the changes occurring in the 3-dimensional architecture of fibrils using biophysical assays and molecular dynamics simulations.

RESULTS

We found that carbamylation aids in amyloid formation and can convert the unstructured off-pathway aggregates into robust amyloids, which were toxic to cells. Electron microscopy images and molecular dynamics simulations of PHF6 fibrils showed that carbamylated peptides can form excess hydrogen bonds and modulate the pitch length and twist of peptides fibrils. We have also compared N-terminal carbamylation to acetylation and further extended our finding to full length tau that exhibits aggregation upon carbamylation even in the absence of any external inducer.

CONCLUSION

Our in vitro and in silico results together suggest that carbamylation can modulate the aggregation pathway of the amyloidegenic sequences and cause structural changes in fibril assemblies.

GENERAL SIGNIFICANCE

Carbamylation acts as a switch, which triggers the aggregation in short amyloidogenic peptide fragments and modulate the structural changes in resulting amyloid fibrils.

Positive effect of strong acidity on the twist of Aβ42 fibrils and the counteraction of Aβ42 N-terminus

pH is a crucial factor in terms of affecting the aggregation and morphology of β-Amyloid and hence a focus of study. In this study, structural and mechanical properties of a series of models (5, 6, …, 30 layer) of one-fold Aβ₄₂ fibrils at pH 1.5, 3.0 and 7.5, have been computed by using all-atom molecular dynamics simulations. 12, 14, and 15 layers are established to be the smallest realistic models for Aβ₄₂ fibrils at pH 1.5, 3.0 and 7.5, with twist angles of 0.40°, 0.34°, 0.31° respectively, disclosing the favorable effect of strong acidity on fibril twist. However, these angles are all lower than that (0.48°) determined for the truncated Aβ_17-42 fibril at pH 7.5, indicating that the disordered N-terminal depresses greatly the fibril twist and the lower pH disfavors the depression. Three commonly used indices to measure the fibril properties, namely number of H-bonds, interstrand distance and β-sheet content have imperceptible changes with the pH alternation, therefore changes in fibril twist can be taken as a probe to monitor fibril properties. By contrast, N-terminus is determined not only to inhibit the U-shaped fibril twist by hampering the stagger between β1 and β2 strands, but also to play a vital carrier role in feeling solution (i.e., pH, salt) changes. These results can help design the nextgeneration of amyloid materials for state-of-the-art bio-nano-med applications by changing the solution pH or modifying chain length.

Cardiolipin dynamics and binding to conserved residues in the mitochondrial ADP/ATP carrier

Cardiolipin in eukaryotes is found in the mitochondrial inner membrane, where it interacts with membrane proteins and, although not essential, is necessary for the optimal activity of a number of proteins. One of them is the mitochondrial ADP/ATP carrier, which imports ADP into the mitochondrion and exports ATP. In the crystal structures, cardiolipin is bound to three equivalent sites of the ADP/ATP carrier, but its role is unresolved. Conservation of residues at these cardiolipin binding sites across other members of the mitochondrial carrier superfamily indicates cardiolipin binding is likely to be important for the function of all mitochondrial carriers. Multiscale simulations were performed in a cardiolipin-containing membrane to investigate the dynamics of cardiolipin around the yeast and bovine ADP/ATP carriers in a lipid bilayer and the properties of the cardiolipin-binding sites. In coarse-grain simulations, cardiolipin molecules bound to the carriers for longer periods of time than phosphatidylcholine and phosphatidylethanolamine lipids—with timescales in the tens of microseconds. Three long-lived cardiolipin binding sites overlapped with those in the crystal structures of the carriers. Other shorter-lived cardiolipin interaction sites were identified in both membrane leaflets. However, the timescales of the interactions were of the same order as phosphatidylcholine and phosphatidylethanolamine, suggesting that these sites are not specific for cardiolipin binding. The calculation of lipid binding times and the overlap of the cardiolipin binding sites between the structures and simulations demonstrate the potential of multiscale simulations to investigate the dynamics and behavior of lipids interacting with membrane proteins.

•

Coarse-grained models of AAC in mixed lipid bilayers were simulated.

•

Three long-lived cardiolipin sites correspond to those in the crystal structures.

•

No other long-lived binding sites were observed for cardiolipin or other phospholipids.

•

Trimethylation of Lys-51 of AAC had no effect on cardiolipin interactions.

Energetics Underlying Twist Polymorphisms in Amyloid Fibrils

Amyloid fibrils are highly ordered protein aggregates associated with more than 40 human diseases. The exact conditions under which the fibrils are grown determine many types of reported fibril polymorphism, including different twist patterns. Twist-based polymorphs display unique mechanical properties in vitro, and the relevance of twist polymorphism in amyloid diseases has been suggested. We present transmission electron microscopy images of Aβ42-derived (amyloid β) fibrils, which are associated with Alzheimer's disease, demonstrating the presence of twist variability even within a single long fibril. To better understand the molecular underpinnings of twist polymorphism, we present a structural and thermodynamics analysis of molecular dynamics simulations of the twisting of β-sheet protofilaments of a well-characterized cross-β model: the GNNQQNY peptide from the yeast prion Sup35. The results show that a protofilament model of GNNQQNY is able to adopt twist angles from -11° on the left-hand side to +8° on the right-hand side in response to various external conditions, keeping an unchanged peptide structure. The potential of mean force (PMF) of this cross-β structure upon twisting revealed that only ∼2kBT per peptide are needed to stabilize a straight conformation with respect to the left-handed free-energy minimum. The PMF also shows that the canonical structural core of β-sheets, i.e., the hydrogen-bonded backbone β-strands, favors the straight conformation. However, the concerted effects of the side chains contribute to twisting, which provides a rationale to correlate polypeptide sequence, environmental growth conditions and number of protofilaments in a fibril with twist polymorphisms.

N-Terminus Binding Preference for Either Tanshinone or Analogue in Both Inhibition of Amyloid Aggregation and Disaggregation of Preformed Amyloid Fibrils-Toward Introducing a Kind of Novel Anti-Alzheimer Compounds

Amyloid-β (Aβ₄₀/Aβ₄₂) peptide with a length of 40 or 42 residues is naturally secreted as cleavage product of the amyloid precursor protein, and formation of Aβ aggregates in a patient's brain is a hallmark of Alzheimer's disease (AD). Therefore, disaggregation and disruption provide potential therapeutic approaches to reduce, inhibit, and even reverse Aβ aggregation. The disaggregation/inhibition effect of the inhibitors applies generally to both Aβ₄₀ and Aβ₄₂ aggregations. Here we capture the atomic-level details of the interaction between Aβ₄₀/Aβ₄₂ and either natural tanshinone compound TS1 or its derivative TS0, and observe novel results by using molecular dynamics simulations. We observe that the natural TS1 indeed inhibits the monomolecular Aβ₄₂ (mAβ₄₂) aggregation and disaggregates Aβ₄₂ amyloid fibrils, being in good agreement with the experimental results. TS1 is favorable to stabilize mAβ₄₀ and even Aβ₄₀ fibril, playing an opposite role to that in the Aβ₄₂ counterpart, however. TS0 can inhibit the misfolding of either mAβ₄₀ or mAβ₄₂ and disaggregate Aβ₄₂ fibril but stabilize the Aβ₄₀ fibril. Using a combination of secondary structural analysis, MM-PBSA binding energy calculations, and radial distribution functions computations, we find that both TS0 and TS1, especially the former, prefer to bind at the charged residues within disordered N-terminus with a scarce positive binding energy and disappear the characteristic C-terminal bend region of Aβ₄₂ fibril, as well as twist the Aβ₄₂ fibril seriously. It turns out to destabilize the Aβ₄₂ fibril and enable the conversion of U-shaped Aβ₄₂ fibril from the onefold to the twofold morphologies. The N-terminal binding preference helps us to identify N-terminal region as the specific epitope for specific inhibitors/drugs (such as TS0 and analogues), heralding unusual inhibition/disaggregation or stabilization mechanisms, and offering an alternative direction in engineering new inhibitors to treat AD.

Cardiolipin binds selectively but transiently to conserved lysine residues in the rotor of metazoan ATP synthases

The anionic lipid cardiolipin is an essential component of active ATP synthases. In metazoans, their rotors contain a ring of eight c-subunits consisting of inner and outer circles of N- and C-terminal α-helices, respectively. The beginning of the C-terminal α-helix contains a strictly conserved and fully trimethylated lysine residue in the lipid head-group region of the membrane. Larger rings of known structure, from c9-c15 in eubacteria and chloroplasts, conserve either a lysine or an arginine residue in the equivalent position. In computer simulations of hydrated membranes containing trimethylated or unmethylated bovine c8-rings and bacterial c10- or c11-rings, the head-groups of cardiolipin molecules became associated selectively with these modified and unmodified lysine residues and with adjacent polar amino acids and with a second conserved lysine on the opposite side of the membrane, whereas phosphatidyl lipids were attracted little to these sites. However, the residence times of cardiolipin molecules with the ring were brief and sufficient for the rotor to turn only a fraction of a degree in the active enzyme. With the demethylated c8-ring and with c10- and c11-rings, the density of bound cardiolipin molecules at this site increased, but residence times were not changed greatly. These highly specific but brief interactions with the rotating c-ring are consistent with functional roles for cardiolipin in stabilizing and lubricating the rotor, and, by interacting with the enzyme at the inlet and exit of the transmembrane proton channel, in participation in proton translocation through the membrane domain of the enzyme.

FRET evidence for untwisting of amyloid fibrils on the surface of model membranes

Apolipoprotein A-I (apoA-I) is an amyloid-forming protein whose amyloidogenic properties are attributed mainly to its N-terminal fragment. Cell membranes are thought to be the primary target for the toxic amyloid aggregates. In the present study Förster resonance energy transfer (FRET) between the membrane fluorescent probe Laurdan as a donor and amyloid-specific dye Thioflavin T (ThT) as an acceptor was employed to explore the interactions of amyloid fibrils from apoA-I variants 1-83/G26R and 1-83/G26R/W@8 with the model membranes composed of phosphatidylcholine and its mixture with cholesterol. The changes in FRET efficiency upon fibril-lipid binding were found to correlate with the extent of protein fibrillization. AFM imaging revealed the presence of two polymorphic states of fibrillar 1-83/G26R/W@8 with the helical and twisted ribbon morphologies. The simulation-based analysis of the experimental FRET profiles provided the arguments in favor of untwisting of fibrillar assemblies upon their interaction with the model membranes. Evidence for the face-on orientation and superficial bilayer location of the membrane-bound fragments of 1-83/G26R/W@8 fibrils was obtained.

The architecture of amyloid-like peptide fibrils revealed by X-ray scattering, diffraction and electron microscopy

Structural analysis of protein fibrillation is inherently challenging. Given the crucial role of fibrils in amyloid diseases, method advancement is urgently needed. A hybrid modelling approach is presented enabling detailed analysis of a highly ordered and hierarchically organized fibril of the GNNQQNY peptide fragment of a yeast prion protein. Data from small-angle X-ray solution scattering, fibre diffraction and electron microscopy are combined with existing high-resolution X-ray crystallographic structures to investigate the fibrillation process and the hierarchical fibril structure of the peptide fragment. The elongation of these fibrils proceeds without the accumulation of any detectable amount of intermediate oligomeric species, as is otherwise reported for, for example, glucagon, insulin and α-synuclein. Ribbons constituted of linearly arranged protofilaments are formed. An additional hierarchical layer is generated via the pairing of ribbons during fibril maturation. Based on the complementary data, a quasi-atomic resolution model of the protofilament peptide arrangement is suggested. The peptide structure appears in a β-sheet arrangement reminiscent of the β-zipper structures evident from high-resolution crystal structures, with specific differences in the relative peptide orientation. The complexity of protein fibrillation and structure emphasizes the need to use multiple complementary methods.

Universal behavior in the mesoscale properties of amyloid fibrils

Amyloid fibrils are ubiquitous proteinaceous aggregates occurring in vivo and in vitro, with an invariant structural fingerprint at the molecular length scale. However, interpretation of their mesoscopic architectures is complicated by diverse observable polymorphic states. We here present a constitutive model for amyloid fibrils based on the minimization of the total energy per fibril. The model is benchmarked on real amyloid fibrils studied by atomic force microscopy. We use multistranded β-lactoglobulin amyloid fibrils as a model system exhibiting a rich polymorphism. The constitutive model quantitatively recapitulates the main mesoscopic topological features of amyloid fibrils, that is, the evolution of fibril periodicity as a function of the ionic strength of the solution and of the fibril width. A universal mesoscopic structural signature of the fibrils emerges from this picture, predicting a general, parameter-free law for the periodicity of the fibrils, that depends solely on the number of protofilaments per fibril. These predictions are validated experimentally and conclusively highlight the role of competing electrostatic and elastic contributions as the main players in the establishment of amyloid fibrils structure.

Energy interactions in amyloid-like fibrils from NNQQNY

We use large-scale MP2 calculations to analyze the interactions appearing in amyloid fibers, which are difficult to determine experimentally. To this end, dimers and trimers of the hexapeptide NNQQNY from the yeast prion-like protein Sup35 were considered as model systems. We studied the energy interactions present in the three levels of organization in which the formation of amyloid fibrils is structured. The structural changes in the hydrogen bonds were studied too. It was found that the most energetic process is the formation of the β-sheet, which is equally due to both hydrogen bonds and van der Waals interactions. The aromatic rings help stabilize these aggregates through stacking of the aromatic rings of tyrosine, the stability produced by the aromatics residues increasing with their aromaticity. The formation of the basic unit of the assembled proto-fiber, the steric zipper, is less energetic and is associated to both dispersion forces and hydrogen bonds. The interactions between pair of β-sheets across the peptide-to-peptide contact through the tyrosine rings are cooperative and due to dispersion effects. Moreover, the strength of this interaction can rationalize the variation of mobility of the aromatic ring in the tyrosine units found in solid NMR experiments.

Interplay of sequence, topology and termini charge in determining the stability of the aggregates of GNNQQNY mutants: a molecular dynamics study

This study explores the stabilities of single sheet parallel systems of three sequence variants of ¹GNNQQNY⁷, N2D, N2S and N6D, with variations in aggregate size (5–8) and termini charge (charged or neutral). The aggregates were simulated at 300 and 330 K. These mutations decrease amyloid formation in the yeast prion protein Sup35. The present study finds that these mutations cause instability even in the peptide context. The protonation status of termini is found to be a key determinant of stabilities; other determinants are sequence, position of mutation and aggregate size. All systems with charged termini are unstable, whereas both stable and unstable systems are found when the termini are neutral. When termini are charged, the largest stable aggregate for the N2S and N6D systems has 3 to 4 peptides whereas N2D mutation supports oligomers of larger size (5-and 6-mers) as well. Mutation at 2^nd position (N2S and N2D) results in fewer H-bonds at the mutated as well as neighboring (Gly1/Gln4) positions. However, no such effect is found if mutation is at 6^th position (N6D). The effect of Asn→Asp mutation depends on the position and termini charge: it is more destabilizing at the 2^nd position than at the 6^th in case of neutral termini, however, the opposite is true in case of charged termini. Appearance of twist in stable systems and in smaller aggregates formed in unstable systems suggests that twist is integral to amyloid arrangement. Disorder, dissociation or rearrangement of peptides, disintegration or collapse of aggregates and formation of amorphous aggregates observed in these simulations are likely to occur during the early stages of aggregation also. The smaller aggregates formed due to such events have a variety of arrangements of peptides. This suggests polymorphic nature of oligomers and presence of a heterogeneous mixture of oligomers during early stages of aggregation.

Factors that drive peptide assembly and fibril formation: experimental and theoretical analysis of Sup35 NNQQNY mutants

Residue mutations have substantial effects on aggregation kinetics and propensities of amyloid peptides and their aggregate morphologies. Such effects are attributed to conformational transitions accessed by various types of oligomers such as steric zipper or single β-sheet. We have studied the aggregation propensities of six NNQQNY mutants: NVVVVY, NNVVNV, NNVVNY, VIQVVY, NVVQIY, and NVQVVY in water using a combination of ion-mobility mass spectrometry, transmission electron microscopy, atomic force microscopy, and all-atom molecular dynamics simulations. Our data show a strong correlation between the tendency to form early β-sheet oligomers and the subsequent aggregation propensity. Our molecular dynamics simulations indicate that the stability of a steric zipper structure can enhance the propensity for fibril formation. Such stability can be attained by either hydrophobic interactions in the mutant peptide or polar side-chain interdigitations in the wild-type peptide. The overall results display only modest agreement with the aggregation propensity prediction methods such as PASTA, Zyggregator, and RosettaProfile, suggesting the need for better parametrization and model peptides for these algorithms.

Effect of Protonation State on the Stability of Amyloid Oligomers Assembled from TTR(105-115)

Amyloid fibrils are self-assembled aggregates of polypeptides that are implicated in the development of several human diseases. A peptide derived from amino acids 105-115 of the human plasma protein transthyretin forms homogeneous and well-defined fibrils and, as a model system, has been the focus of a number of studies investigating the formation and structure of this class of aggregates. Self-assembly of TTR(105-115) occurs at low pH, and this work explores the effect of protonation on the growth and stability of small cross-β aggregates. Using molecular dynamics simulations of structures up to the decamer in both protonated and deprotonated states, we find that, whereas hexamers are more stable for protonated peptides, higher order oligomers are more stable when the peptides are deprotonated. Our findings imply a change in the acid pK of the protonated C-terminal group during the formation of fibrils, which leads to stabilization of higher-order oligomers through electrostatic interactions.

The self-assembly, aggregation and phase transitions of food protein systems in one, two and three dimensions

The aggregation of proteins is of fundamental relevance in a number of daily phenomena, as important and diverse as blood coagulation, medical diseases, or cooking an egg in the kitchen. Colloidal food systems, in particular, are examples that have great significance for protein aggregation, not only for their importance and implications, which touches on everyday life, but also because they allow the limits of the colloidal science analogy to be tested in a much broader window of conditions, such as pH, ionic strength, concentration and temperature. Thus, studying the aggregation and self-assembly of proteins in foods challenges our understanding of these complex systems from both the molecular and statistical physics perspectives. Last but not least, food offers a unique playground to study the aggregation of proteins in three, two and one dimensions, that is to say, in the bulk, at air/water and oil/water interfaces and in protein fibrillation phenomena. In this review we will tackle this very ambitious task in order to discuss the current understanding of protein aggregation in the framework of foods, which is possibly one of the broadest contexts, yet is of tremendous daily relevance.

Driving forces and structural determinants of steric zipper peptide oligomer formation elucidated by atomistic simulations

Understanding the structural and energetic requirements of non-fibrillar oligomer formation harbors the potential to decipher an important yet still elusive part of amyloidogenic peptide and protein aggregation. Low-molecular-weight oligomers are described to be transient and polymorphic intermediates in the nucleated self-assembly process to highly ordered amyloid fibers and were additionally found to exhibit a profound cytotoxicity. However, detailed structural information on the oligomeric species involved in the nucleation cannot be readily inferred from experiments. Here, we study the spontaneous assembly of steric zipper peptides from the tau protein, insulin and α-synuclein with atomistic molecular dynamics simulations on the microsecond timescale. Detailed analysis of the forces driving the oligomerization reveals a common two-step process akin to a general condensation-ordering mechanism and thus provides a rational understanding of the molecular basis of peptide self-assembly. Our results suggest that the initial formation of partially ordered peptide oligomers is governed by the solvation free energy, whereas the dynamical ordering and emergence of β-sheets are mainly driven by optimized inter-peptide interactions in the collapsed state. A novel mapping technique based on collective coordinates is employed to highlight similarities and differences in the conformational ensemble of small oligomer structures. Elucidating the dynamical and polymorphic β-sheet oligomer conformations at atomistic detail furthermore suggests complementary sheet packing characteristics similar to steric zipper structures, but with a larger heterogeneity in the strand alignment pattern and sheet-to-sheet arrangements compared to the cross-β motif found in the fibrillar or crystalline states.

Structural complexity of a composite amyloid fibril

The molecular structure of amyloid fibrils and the mechanism of their formation are of substantial medical and biological importance, but present an ongoing experimental and computational challenge. An early high-resolution view of amyloid-like structure was obtained on amyloid-like crystals of a small fragment of the yeast prion protein Sup35p: the peptide GNNQQNY. As GNNQQNY also forms amyloid-like fibrils under similar conditions, it has been theorized that the crystal's structural features are shared by the fibrils. Here we apply magic-angle-spinning (MAS) NMR to examine the structure and dynamics of these fibrils. Previously multiple NMR signals were observed for such samples, seemingly consistent with the presence of polymorphic fibrils. Here we demonstrate that peptides with these three distinct conformations instead assemble together into composite protofilaments. Electron microscopy (EM) of the ribbon-like fibrils indicates that these protofilaments combine in differing ways to form striations of variable widths, presenting another level of structural complexity. Structural and dynamic NMR data reveal the presence of highly restricted side-chain conformations involved in interfaces between differently structured peptides, likely comprising interdigitated steric zippers. We outline molecular interfaces that are consistent with the observed EM and NMR data. The rigid and uniform structure of the GNNQQNY crystals is found to contrast distinctly with the more complex structural and dynamic nature of these "composite" amyloid fibrils. These results provide insight into the fibril-crystal distinction and also indicate a necessary caution with respect to the extrapolation of crystal structures to the study of fibril structure and formation.

Self-organizing bioinspired oligothiophene-oligopeptide hybrids

In this minireview, we survey recent advances in the synthesis, characterization, and modeling of new oligothiophene-oligopeptide hybrids capable of forming nanostructured fibrillar aggregates in solution and on solid substrates. Compounds of this class are promising for applications because their self-assembly and stimuli-responsive properties, provided by the peptide moieties combined with the semiconducting properties of the thiophene blocks, can result in novel opportunities for the design of advanced smart materials. These bio-inspired molecular hybrids are experimentally shown to form stable fibrils as visualized by AFM and TEM. While the experimental evidence alone is not sufficient to reveal the exact molecular organization of the fibrils, theoretical approaches based on quantum chemistry calculations and large-scale atomistic molecular dynamics simulations are attempted in an effort to reveal the structure of the fibrils at the nanoscale. Based on the combined theoretical and experimental analysis, the most likely models of fibril formation and aggregation are suggested.

Nanomechanics of functional and pathological amyloid materials

Amyloid or amyloid-like fibrils represent a general class of nanomaterials that can be formed from many different peptides and proteins. Although these structures have an important role in neurodegenerative disorders, amyloid materials have also been exploited for functional purposes by organisms ranging from bacteria to mammals. Here we review the functional and pathological roles of amyloid materials and discuss how they can be linked back to their nanoscale origins in the structure and nanomechanics of these materials. We focus on insights both from experiments and simulations, and discuss how comparisons between functional protein filaments and structures that are assembled abnormally can shed light on the fundamental material selection criteria that lead to evolutionary bias in multiscale material design in nature.

Unique example of amyloid aggregates stabilized by main chain H-bond instead of the steric zipper: molecular dynamics study of the amyloidogenic segment of amylin wild-type and mutants

Most proteins do not aggregate while in their native functional states. However, they may be disturbed from their native conformation by certain change in the environment, and form unwanted oligomeric or polymeric aggregates. Recent experimental data demonstrate that soluble oligomers of amyloidogenic proteins are responsible for amyloidosis and its cytotoxicity. Human islet amyloid polypeptide (IAPP or amylin) is a 37-residue hormone found as fibrillar deposits in pancreatic extracts of nearly all type II diabetics. In this study we performed in silico mutation analysis to examine the stability of the double layer five strand aggregates formed by heptapeptide NNFGAIL segment from amyline peptide. This segment is one of the shortest fragments that can form amyloid fibrils similar to those formed by the full length peptide. The mutants obtained by single glycine replacement were also studied to investigate the specificity of the dry self-complementary interface between the neighboring β-sheet layers. The molecular dynamics simulations of the aggregates run for 20 ns at 330 K, the degree of the aggregate disassembly was investigated using several geometry analysis tools: the root mean square deviations of the C(α) atoms, root mean square fluctuations per residue, twist angles, interstrand distances, fraction of the secondary structure elements, and number of H-bonds. The analysis shows that most mutations make the aggregates unstable, and their stabilities were dependent to a large extent on the position of replaced residues. Our mutational simulations are in agreement with the pervious experimental observations. We also used free binding energy calculations to determine the role of different components: nonpolar effects, electrostatics and entropy in binding. Nonpolar effects remained consistently more favorable in wild type and mutants reinforcing the importance of hydrophobic effects in protein-protein binding. While entropy systematically opposed binding in all cases, there was no clear trend in the entropy difference between wildtype and glycine mutants. Free energy decomposition shows residues situated at the interface were found to make favorable contributions to the peptide-peptide association. The study of the wild type and mutants in an explicit solvent could provide valuable insight into the future computer guided design efforts for the amyloid aggregation inhibitor.

A multiscale approach to characterize the early aggregation steps of the amyloid-forming peptide GNNQQNY from the yeast prion sup-35

The self-organization of peptides into amyloidogenic oligomers is one of the key events for a wide range of molecular and degenerative diseases. Atomic-resolution characterization of the mechanisms responsible for the aggregation process and the resulting structures is thus a necessary step to improve our understanding of the determinants of these pathologies. To address this issue, we combine the accelerated sampling properties of replica exchange molecular dynamics simulations based on the OPEP coarse-grained potential with the atomic resolution description of interactions provided by all-atom MD simulations, and investigate the oligomerization process of the GNNQQNY for three system sizes: 3-mers, 12-mers and 20-mers. Results for our integrated simulations show a rich variety of structural arrangements for aggregates of all sizes. Elongated fibril-like structures can form transiently in the 20-mer case, but they are not stable and easily interconvert in more globular and disordered forms. Our extensive characterization of the intermediate structures and their physico-chemical determinants points to a high degree of polymorphism for the GNNQQNY sequence that can be reflected at the macroscopic scale. Detailed mechanisms and structures that underlie amyloid aggregation are also provided.

The formation of amyloid fibrils is associated with many neurodegenerative diseases such as Alzheimer's, Creutzfeld-Jakob, Parkinson's, the Prion disease and diabetes mellitus. In all cases, proteins misfold to form highly ordered insoluble aggregates called amyloid fibrils that deposit intra- and extracellularly and are resistant to proteases. All proteins are believed to have the instrinsic capability of forming amyloid fibrils that share common specific structural properties that have been observed by X-ray crystallography and by NMR. However, little is known about the aggregation dynamics of amyloid assemblies, and their toxicity mechanism is therefore poorly understood. It is believed that small amyloid oligomers, formed on the aggregation pathway of full amyloid fibrils, are the toxic species. A detailed atomic characterization of the oligomerization process is thus necessary to further our understanding of the amyloid oligomer's toxicity. Our approach here is to study the aggregation dynamics of a 7-residue amyloid peptide GNNQQNY through a combination of numerical techniques. Our results suggest that this amyloid sequence can form fibril-like structures and is polymorphic, which agrees with recent experimental observations. The ability to fully characterize and describe the aggregation pathway of amyloid sequences numerically is key to the development of future drugs to target amyloid oligomers.

Structural characterization of GNNQQNY amyloid fibrils by magic angle spinning NMR

Several human diseases are associated with the formation of amyloid aggregates, but experimental characterization of these amyloid fibrils and their oligomeric precursors has remained challenging. Experimental and computational analysis of simpler model systems has therefore been necessary, for instance, on the peptide fragment GNNQQNY7−13 of yeast prion protein Sup35p. Expanding on a previous publication, we report here a detailed structural characterization of GNNQQNY fibrils using magic angle spinning (MAS) NMR. On the basis of additional chemical shift assignments we confirm the coexistence of three distinct peptide conformations within the fibrillar samples, as reflected in substantial chemical shift differences. Backbone torsion angle measurements indicate that the basic structure of these coexisting conformers is an extended β-sheet. We structurally characterize a previously identified localized distortion of the β-strand backbone specific to one of the conformers. Intermolecular contacts are consistent with each of the conformers being present in its own parallel and in-register sheet. Overall the MAS NMR data indicate a substantial difference between the structure of the fibrillar and crystalline forms of these peptides, with a clearly increased complexity in the GNNQQNY fibril structure. These experimental data can provide guidance for future work, both experimental and theoretical, and provide insights into the distinction between fibril growth and crystal formation.

Alzheimer's abeta(1-40) amyloid fibrils feature size-dependent mechanical properties

Amyloid fibrils are highly ordered protein aggregates that are associated with several pathological processes, including prion propagation and Alzheimer's disease. A key issue in amyloid science is the need to understand the mechanical properties of amyloid fibrils and fibers to quantify biomechanical interactions with surrounding tissues, and to identify mechanobiological mechanisms associated with changes of material properties as amyloid fibrils grow from nanoscale to microscale structures. Here we report a series of computational studies in which atomistic simulation, elastic network modeling, and finite element simulation are utilized to elucidate the mechanical properties of Alzheimer's Abeta(1-40) amyloid fibrils as a function of the length of the protein filament for both twofold and threefold symmetric amyloid fibrils. We calculate the elastic constants associated with torsional, bending, and tensile deformation as a function of the size of the amyloid fibril, covering fibril lengths ranging from nanometers to micrometers. The resulting Young's moduli are found to be consistent with available experimental measurements obtained from long amyloid fibrils, and predicted to be in the range of 20-31 GPa. Our results show that Abeta(1-40) amyloid fibrils feature a remarkable structural stability and mechanical rigidity for fibrils longer than approximately 100 nm. However, local instabilities that emerge at the ends of short fibrils (on the order of tens of nanometers) reduce their stability and contribute to their disassociation under extreme mechanical or chemical conditions, suggesting that longer amyloid fibrils are more stable. Moreover, we find that amyloids with lengths shorter than the periodicity of their helical pitch, typically between 90 and 130 nm, feature significant size effects of their bending stiffness due the anisotropy in the fibril's cross section. At even smaller lengths (50 nm), shear effects dominate lateral deformation of amyloid fibrils, suggesting that simple Euler-Bernoulli beam models fail to describe the mechanics of amyloid fibrils appropriately. Our studies reveal the importance of size effects in elucidating the mechanical properties of amyloid fibrils. This issue is of great importance for comparing experimental and simulation results, and gaining a general understanding of the biological mechanisms underlying the growth of ectopic amyloid materials.

HRMAS 1H NMR conformational study of the resin-bound amyloid-forming peptide GNNQQNY from the yeast prion Sup35

The conversion of soluble proteins to insoluble amyloid fibrils is associated with numerous human diseases. The peptide GNNQQNY is a short segment of the yeast prion protein Sup35 that previously has been found to form amyloid fibrils in a similar manner to the protein itself. The approach taken in this work was to attach this peptide sequence to an insoluble polymer matrix through solid phase peptide synthesis and give it the internal freedom to fold into its local conformation in an organic solvent. Observation of its monomeric structure, free from the effects of aggregation, entropic solvent effects, and neighboring molecules, was possible by two-dimensional high-resolution magic angle spinning (1)H NMR spectroscopy. Analysis of the through-bond correlations and through-space interactions observed in the spectra, combined with global energy minimization via computational studies, led to the observation that the peptide chain adopts a compact beta-like turn at the central hydrophilic residues. The technique of peptide attachment to a polymer resin and observation by NMR may allow for future study of single peptide fragments prone to aggregation.

In silico study of full-length amyloid beta 1-42 tri- and penta-oligomers in solution

Amyloid oligomers are considered to play causal roles in the pathogenesis of amyloid-related degenerative diseases including Alzheimer's disease. Using MD simulation techniques, we explored the contributions of the different structural elements of trimeric and pentameric full-length Abeta1-42 aggregates in solution to their stability and conformational dynamics. We found that our models are stable at a temperature of 310 K, and converge toward an interdigitated side-chain packing for intermolecular contacts within the two beta-sheet regions of the aggregates: beta1 (residues 18-26) and beta2 (residues 31-42). MD simulations reveal that the beta-strand twist is a characteristic element of Abeta-aggregates, permitting a compact, interdigitated packing of side chains from neighboring beta-sheets. The beta2 portion formed a tightly organized beta-helix, whereas the beta1 portion did not show such a firm structural organization, although it maintained its beta-sheet conformation. Our simulations indicate that the hydrophobic core comprising the beta2 portion of the aggregate is a crucial stabilizing element in the Abeta aggregation process. On the basis of these structure-stability findings, the beta2 portion emerges as an optimal target for further antiamyloid drug design.

Atomic-scale simulations confirm that soluble beta-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 beta-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 beta-sheet-rich assemblies. We simulate PSAMs with different-sized beta-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 beta-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 beta-sheets to some extent, implying that a twisted and bent beta-sheet is a general motif of beta-rich assemblies. Thus, our studies indicate that soluble beta-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.