Annular structures as intermediates in fibril formation of Alzheimer Abeta17-42

Single-Molecule Maps of Membrane Insertion by Amyloid-β Oligomers Predict Their Toxicity

The interaction of small Amyloid-β (Aβ) oligomers with the lipid membrane is an important component of the pathomechanism of Alzheimer's disease (AD). However, oligomers are heterogeneous in size. How each type of oligomer incorporates into the membrane, and how that relates to their toxicity, is unknown. Here, we employ a single molecule technique called Q-SLIP (Quencher-induced Step Length Increase in Photobleaching) to measure the membrane insertion of each monomeric unit of individual oligomers of Aβ42, Aβ40, and Aβ40-F19-Cyclohexyl alanine (Aβ40-F19Cha), and correlate it with their toxicity. We observe that the N-terminus of Aβ42 inserts close to the center of the bilayer, the less toxic Aβ40 inserts to a shallower depth, and the least toxic Aβ40-F19Cha has no specific distribution. This oligomer-specific map provides a mechanistic representation of membrane-mediated Aβ toxicity and should be a valuable tool for AD research.

Early aggregation mechanism of Aβ16-22 revealed by Markov state models

Aβ_16-22 is believed to have critical role in early aggregation of full length amyloids that are associated with the Alzheimer's disease and can aggregate to form amyloid fibrils. However, the early aggregation mechanism is still unsolved. Here, multiple long-term molecular dynamics simulations combining with Markov state model were used to probe the early oligomerization mechanism of Aβ_16-22 peptides. The identified dimeric form adopted either globular random-coil or extended β-strand like conformations. The observed dimers of these variants shared many overall conformational characteristics but differed in several aspects at detailed level. In all cases, the most common type of secondary structure was intermolecular antiparallel β-sheets. The inter-state transitions were very frequent ranges from few to hundred nanoseconds. More strikingly, those states which contain fraction of β secondary structure and significant amount of extended coiled structures, therefore exposed to the solvent, were majorly participated in aggregation. The assembly of low-energy dimers, in which the peptides form antiparallel β sheets, occurred by multiple pathways with the formation of an obligatory intermediates. We proposed that these states might facilitate the Aβ_16-22 aggregation through a significant component of the conformational selection mechanism, because they might increase the aggregates population by promoting the inter-chain hydrophobic and the hydrogen bond contacts. The formation of early stage antiparallel β sheet structures is critical for oligomerization, and at the same time provided a flat geometry to seed the ordered β-strand packing of the fibrils. Our findings hint at reorganization of this part of the molecule as a potentially critical step in Aβ aggregation and will insight into early oligomerization for large β amyloids.

Defining the Neuropathological Aggresome across in Silico, in Vitro, and ex Vivo Experiments

The loss of proteostasis over the life course is associated with a wide range of debilitating degenerative diseases and is a central hallmark of human aging. When left unchecked, proteins that are intrinsically disordered can pathologically aggregate into highly ordered fibrils, plaques, and tangles (termed amyloids), which are associated with countless disorders such as Alzheimer's disease, Parkinson's disease, type II diabetes, cancer, and even certain viral infections. However, despite significant advances in protein folding and solution biophysics techniques, determining the molecular cause of these conditions in humans has remained elusive. This has been due, in part, to recent discoveries showing that soluble protein oligomers, not insoluble fibrils or plaques, drive the majority of pathological processes. This has subsequently led researchers to focus instead on heterogeneous and often promiscuous protein oligomers. Unfortunately, significant gaps remain in how to prepare, model, experimentally corroborate, and extract amyloid oligomers relevant to human disease in a systematic manner. This Review will report on each of these techniques and their successes and shortcomings in an attempt to standardize comparisons between protein oligomers across disciplines, especially in the context of neurodegeneration. By standardizing multiple techniques and identifying their common overlap, a clearer picture of the soluble neuropathological aggresome can be constructed and used as a baseline for studying human disease and aging.

Modulating amyloid fibrillation in a minimalist model peptide by intermolecular disulfide chemical reduction

Peptide structural transformation and aggregation is associated with a large number of outsider aetiology diseases, and it is intrinsically linked to amyloid peptide aggregation. Diphenylalanine self-assembled structures are used as robust minimalist beta amyloids not only to elucidate protein aggregation but also to generate hydrogels. Herein, we employed a neutral model peptide Ac-Phe-Phe-Cys-NH2 (Ac-FFC-NH2) to elucidate the role of intermolecular disulfide bonds in protein fibrillation. The Ac-FFC-NH2 peptide initially self-assembles into nanospheres that evolve to amyloid type fibrils under mild oxidative conditions. Incubation of the peptide in the presence of the chemical reduction agent TCEP inhibits the formation of the fibrils, detecting only spherical nanostructures with no secondary structure. Importantly, we triggered the transformation of the preformed linear straight amyloid fibrils to non-fibrillar structures by TCEP treatment. Under this condition, the amyloid bundles are transformed into rings, which evolve to a new spherical microstructure. We showed that the chemical reduction of intermolecular S-S in internal amyloid sequences might favour the off-path intermediates of amyloid fibril growth, even when the fibrils are formed. Our findings demonstrated that in internal amyloid sequences, the formation of intermolecular S-S promotes the formation of amyloid type fibrils; meanwhile, its reduction stabilises non-fibrillar structures. Altogether, this work provides fundamental understanding at the molecular and supramolecular level, thus facilitating the rational design of therapeutic tools for protein aggregation diseases.

Solvent Composition Effects on the Structural Properties of the Aβ42 Monomer from the 3D-RISM-KH Molecular Theory of Solvation

Structural characterization of amyloid (A)β peptides implicated in Alzheimer's disease is a challenging problem due to their intrinsically disordered nature and their high propensity for aggregation. Only limited information is currently available from experiments on conformational properties and aggregation pathways of the peptides in cellular environments. In silico modeling complements experimental information, providing atomistic insight into structure and dynamics of different Aβ species. All-atom explicit solvent molecular dynamics (MD) simulations with a properly selected force field can deliver reliable structural and dynamic information. In the case of intrinsically disordered Aβ peptides, enhanced sampling simulations beyond the nanosecond time scale are required to obtain statistically meaningful results even for simple solvent conditions. To overcome the challenges of conformational sampling in crowded cellular environments, alternative approaches have to be used, including postprocessing of MD data. In this study, we employ the statistical-mechanical, three-dimensional reference interaction site model with the Kovalenko-Hirata closure integral equation molecular theory of solvation to describe solvent composition effects on the conformational equilibrium in a structural ensemble of the Aβ42 (covering residues 1-42) monomer based on a statistical reweighting technique. The methodology enables a computationally efficient prediction on how different factors in the cellular environment, such as solvent composition, nonpolar solvation, and macromolecular crowding, affect the structural properties of the monomer. Similarities have been identified between changes in the structural ensemble caused by nonpolar solvation and crowded environments modeled by ionic solution with large negative ions. In particular, both solvent conditions reduce the random coil content and enhance the helical structure content of the monomer. In contrast to the previous studies, which reported increased α-helical content of peptides in crowded environments, this work attributes these structural features to the difference in solvent exposure of hydrophilic residues of the monomer for different secondary structure elements, rather than to (entropic) excluded volume effects.

Structure of amyloid β25-35 in lipid environment and cholesterol-dependent membrane pore formation

The amyloid β (Aβ) peptide and its shorter variants, including a highly cytotoxic Aβ_25–35 peptide, exert their neurotoxic effect during Alzheimer’s disease by various mechanisms, including cellular membrane permeabilization. The intrinsic polymorphism of Aβ has prevented the identification of the molecular basis of Aβ pore formation by direct structural methods, and computational studies have led to highly divergent pore models. Here, we have employed a set of biophysical techniques to directly monitor Ca²⁺-transporting Aβ_25–35 pores in lipid membranes, to quantitatively characterize pore formation, and to identify the key structural features of the pore. Moreover, the effect of membrane cholesterol on pore formation and the structure of Aβ_25–35 has been elucidated. The data suggest that the membrane-embedded peptide forms 6- or 8-stranded β-barrel like structures. The 8-stranded barrels may conduct Ca²⁺ ions through an inner cavity, whereas the tightly packed 6-stranded barrels need to assemble into supramolecular structures to form a central pore. Cholesterol affects Aβ_25–35 pore formation by a dual mechanism, i.e., by direct interaction with the peptide and by affecting membrane structure. Collectively, our data illuminate the molecular basis of Aβ membrane pore formation, which should advance both basic and clinical research on Alzheimer’s disease and membrane-associated pathologies in general.

Dihydrochalcone molecules destabilize Alzheimer's amyloid-β protofibrils through binding to the protofibril cavity

Dihydrochalcone molecules destabilize Aβ_17–42protofibrils by disrupting the N-terminal β1 region and the turn region through binding to the protofibril cavity.

Investigation of the interaction of amyloid β peptide (11–42) oligomers with a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) membrane using molecular dynamics simulation

The mechanism of pore formation in model neural cell membranes by β amyloid (Aβ) peptides was investigated using molecular dynamics simulation which indicated that Aβ oligomers of size equal or greater than 3 has a higher tendency for pore formation than monomers and that cholesterol tends to retard Aβ binding and insertion into the membrane.

Role of the N-terminus for the stability of an amyloid-β fibril with three-fold symmetry

A key player in Alzheimer’s disease is the peptide amyloid-beta (Aβ), whose aggregation into small soluble oligomers, protofilaments, and fibrils finally leads to plaque deposits in human brains. The aggregation behavior of Aβ is strongly modulated by the nature and composition of the peptide’s environment and by its primary sequence properties. The N-terminal residues of Aβ play an important role, because they are known to change the peptide’s aggregation propensity. Since these residues are for the first time completely resolved at the molecular level in a three-fold symmetric fibril structure derived from a patient, we chose that system as template for a systematic investigation of the influence of the N-terminus upon structural stability. Using atomistic molecular dynamics simulations, we examined several fibrillar systems comprising three, six, twelve and an infinite number of layers, both with and without the first eight residues. First, we found that three layers are not sufficient to stabilize the respective Aβ topology. Second, we observed a clear stabilizing effect of the N-terminal residues upon the overall fibril fold: truncated Aβ systems were less stable than their full-length counterparts. The N-terminal residues Arg5, Asp7, and Ser8 were found to form important interfilament contacts stabilizing the overall fibril structure of three-fold symmetry. Finally, similar structural rearrangements of the truncated Aβ species in different simulations prompted us to suggest a potential mechanism involved in the formation of amyloid fibrils with three-fold symmetry.

Probing oligomerization of amyloid beta peptide in silico

Aggregation of amyloid β (Aβ) peptide is implicated in fatal Alzheimer's disease, for which no cure is available. Understanding the mechanisms responsible for this aggregation is required in order for therapies to be developed. In an effort to better understand the molecular mechanisms involved in spontaneous aggregation of Aβ peptide, extensive molecular dynamics simulations are reported, and the results are analyzed through a combination of structural biology tools and a novel essential collective dynamics method. Several model systems composed of ten or twelve Aβ_17-42 chains in water are investigated, and the influence of metal ions is probed. The results suggest that Aβ monomers tend to aggregate into stable globular-like oligomers with 13-23% of β-sheet content. Two stages of oligomer formation have been identified: quick collapse within the first 40 ns of the simulation, characterized by a decrease in inter-chain separation and build-up of β-sheets, and the subsequent slow relaxation of the oligomer structure. The resulting oligomers comprise a stable, coherently moving sub-aggregate of 6-9 strongly inter-correlated chains. Cu²⁺ and Fe²⁺ ions have been found to develop coordination bonds with carboxylate groups of E22, D23 and A42, which remain stable during 200 ns simulations. The presence of Fe²⁺, and particularly Cu²⁺ ions, in negatively charged cavities has been found to cause significant changes in the structure and dynamics of the oligomers. The results indicate, in particular, that formation of non-fibrillar oligomers might be involved in early template-free aggregation of Aβ_17-42 monomers, with charged species such as Cu²⁺ or Fe²⁺ ions playing an important role.

Influence of Hyperglycemic Conditions on Self-Association of the Alzheimer's Amyloid β (Aβ1-42) Peptide

Clinical studies have identified a correlation between type-2 diabetes mellitus and cognitive decrements en route to the onset of Alzheimer's disease (AD). Recent studies have established that post-translational modifications of the amyloid β (Aβ) peptide occur under hyperglycemic conditions; particularly, the process of glycation exacerbates its neurotoxicity and accelerates AD progression. In view of the assertion that macromolecular crowding has an altering effect on protein self-assembly, it is crucial to characterize the effects of hyperglycemic conditions via crowding on Aβ self-assembly. Toward this purpose, fully atomistic molecular dynamics simulations were performed to study the effects of glucose crowding on Aβ dimerization, which is the smallest known neurotoxic species. The dimers formed in the glucose-crowded environment were found to have weaker associations as compared to that of those formed in water. Binding free energy calculations show that the reduced binding strength of the dimers can be mainly attributed to the overall weakening of the dispersion interactions correlated with substantial loss of interpeptide contacts in the hydrophobic patches of the Aβ units. Analysis to discern the differential solvation pattern in the glucose-crowded and pure water systems revealed that glucose molecules cluster around the protein, at a distance of 5-7 Å, which traps the water molecules in close association with the protein surface. This preferential exclusion of glucose molecules and resulting hydration of the Aβ peptides has a screening effect on the hydrophobic interactions, which in turn diminishes the binding strength of the resulting dimers. Our results imply that physical effects attributed to crowded hyperglycemic environments are incapable of solely promoting Aβ self-assembly, indicating that further mechanistic studies are required to provide insights into the self-assembly of post-translationally modified Aβ peptides, known to possess aggravated toxicity, under these conditions.

A Protocol for the Design of Protein and Peptide Nanostructure Self-Assemblies Exploiting Synthetic Amino Acids

In recent years there has been increasing interest in nanostructure design based on the self-assembly properties of proteins and polymers. Nanodesign requires the ability to predictably manipulate the properties of the self-assembly of autonomous building blocks, which can fold or aggregate into preferred conformational states. The design includes functional synthetic materials and biological macromolecules. Autonomous biological building blocks with available 3D structures provide an extremely rich and useful resource. Structural databases contain large libraries of protein molecules and their building blocks with a range of sizes, shapes, surfaces, and chemical properties. The introduction of engineered synthetic residues or short peptides into these building blocks can greatly expand the available chemical space and enhance the desired properties. Herein, we summarize a protocol for designing nanostructures consisting of self-assembling building blocks, based on our recent works. We focus on the principles of nanostructure design with naturally occurring proteins and synthetic amino acids, as well as hybrid materials made of amyloids and synthetic polymers.

Polymorphism of fibrillar structures depending on the size of assembled Aβ17-42 peptides

The size of assembled Aβ_17-42 peptides can determine polymorphism during oligomerization and fibrillization, but the mechanism of this effect is unknown. Starting from separate random monomers, various fibrillar oligomers with distinct structural characteristics were identified using discontinuous molecular dynamics simulations based on a coarse-grained protein model. From the structures observed in the simulations, two characteristic oligomer sizes emerged, trimer and paranuclei, which generated distinct structural patterns during fibrillization. A majority of the simulations for trimers and tetramers formed non-fibrillar oligomers, which primarily progress to off-pathway oligomers. Pentamers and hexamers were significantly converted into U-shape fibrillar structures, meaning that these oligomers, called paranuclei, might be potent on-pathway intermediates in fibril formation. Fibrillar oligomers larger than hexamers generated substantial polymorphism in which hybrid structures were readily formed and homogeneous fibrillar structures appeared infrequently.

Computational Methods for Structural and Functional Studies of Alzheimer's Amyloid Ion Channels

Aggregation can be studied by a range of methods, experimental and computational. Aggregates form in solution, across solid surfaces, and on and in the membrane, where they may assemble into unregulated leaking ion channels. Experimental probes of ion channel conformations and dynamics are challenging. Atomistic molecular dynamics (MD) simulations are capable of providing insight into structural details of amyloid ion channels in the membrane at a resolution not achievable experimentally. Since data suggest that late stage Alzheimer's disease involves formation of toxic ion channels, MD simulations have been used aiming to gain insight into the channel shapes, morphologies, pore dimensions, conformational heterogeneity, and activity. These can be exploited for drug discovery. Here we describe computational methods to model amyloid ion channels containing the β-sheet motif at atomic scale and to calculate toxic pore activity in the membrane.

Electric Field Effects on Short Fibrils of Aβ Amyloid Peptides

Amyloid fibrils are highly ordered protein aggregates, which are associated with many neurodegenerative diseases. The assembling dynamics of monomeric beta-amyloid peptides, Aβ, into small aggregates (and then into long fibrils) is still debated and has become a hot topic. In this study, we conducted molecular dynamics simulations in explicit water of small Aβ protofibrils (from monomer to pentamer) under the perturbation of an externally applied electric field with the aim of investigating the fundamental molecular interactions involved in the aggregation mechanism. Dynamics of small adducts of Aβ(16-42) in the presence of an electric field, which was shown before to accelerate the conformational change of a single molecule, indicate that the structural resilience increases with the number of molecules in the aggregate. In particular, for 50 ns, the pentamer shows an enhanced stability in secondary structure, number of hydrogen bonds, and number of salt bridges, even in the presence of the field perturbation. The resilience to the field perturbation is linked to the variation of the induced dipole moment of the aggregates that tends to level off very rapidly with the growing number of molecules, thereby reducing the energy available per molecule to produce structural changes. The results also show that in the presence of the field the stability of the hydrophobic second β-sheet (β2, residues 31-42) is higher than that of the first one (β1, residues 18-26). In particular, we identify Gly33, Gly37, and Met35 as the most important residues that stabilize the intermolecular packing and may act as nucleation sites for fibrillization. Furthermore, dynamics of the full-length Aβ(1-42) pentameric aggregate, which include the highly charged random coil residues 1-15, confirmed the key role of the second hydrophobic core in the protofibril structure.

Structural Conversion of Aβ17-42 Peptides from Disordered Oligomers to U-Shape Protofilaments via Multiple Kinetic Pathways

Discovering the mechanisms by which proteins aggregate into fibrils is an essential first step in understanding the molecular level processes underlying neurodegenerative diseases such as Alzheimer's and Parkinson's. The goal of this work is to provide insights into the structural changes that characterize the kinetic pathways by which amyloid-β peptides convert from monomers to oligomers to fibrils. By applying discontinuous molecular dynamics simulations to PRIME20, a force field designed to capture the chemical and physical aspects of protein aggregation, we have been able to trace out the entire aggregation process for a system containing 8 Aβ17-42 peptides. We uncovered two fibrillization mechanisms that govern the structural conversion of Aβ17-42 peptides from disordered oligomers into protofilaments. The first mechanism is monomeric conversion templated by a U-shape oligomeric nucleus into U-shape protofilament. The second mechanism involves a long-lived and on-pathway metastable oligomer with S-shape chains, having a C-terminal turn, en route to the final U-shape protofilament. Oligomers with this C-terminal turn have been regarded in recent experiments as a major contributing element to cell toxicity in Alzheimer's disease. The internal structures of the U-shape protofilaments from our PRIME20/DMD simulation agree well with those from solid state NMR experiments. The approach presented here offers a simple molecular-level framework to describe protein aggregation in general and to visualize the kinetic evolution of a putative toxic element in Alzheimer's disease in particular.

Molecular understanding of a potential functional link between antimicrobial and amyloid peptides

Antimicrobial and amyloid peptides do not share common sequences, typical secondary structures, or normal biological activity but both the classes of peptides exhibit membrane-disruption ability to induce cell toxicity. Different membrane-disruption mechanisms have been proposed for antimicrobial and amyloid peptides, individually, some of which are not exclusive to either peptide type, implying that certain common principles may govern the folding and functions of different cytolytic peptides and associated membrane disruption mechanisms. Particularly, some antimicrobial and amyloid peptides have been identified to have dual complementary amyloid and antimicrobial properties, suggesting a potential functional link between amyloid and antimicrobial peptides. Given that some similar structural and membrane-disruption characteristics exist between the two classes of peptides, this review summarizes major findings, recent advances, and future challenges related to antimicrobial and amyloid peptides and strives to illustrate the similarities, differences, and relationships in the sequences, structures, and membrane interaction modes between amyloid and antimicrobial peptides, with a special focus on direct interactions of the peptides with the membranes. We hope that this review will stimulate further research at the interface of antimicrobial and amyloid peptides - which has been studied less intensively than either type of peptides - to decipher a possible link between both amyloid pathology and antimicrobial activity, which can guide drug design and peptide engineering to influence peptide-membrane interactions important in human health and diseases.

Structural diversity of Alzheimer's disease amyloid-β dimers and their role in oligomerization and fibril formation

Alzheimer's disease (AD) is associated with the formation of toxic amyloid-β (Aβ)42 oligomers, and recent evidence supports a role for Aβ dimers as building blocks for oligomers. Molecular dynamics simulation studies have identified clans for the dominant conformations of Aβ42 forming dimers; however, it is unclear if a larger spectrum of dimers is involved and which set(s) of dimers might evolve to oligomers verse fibrils. Therefore, for this study we generated multiple structural conformations of Aβ42, using explicit all-atom molecular dynamics, and then clustering the different structures based on key conformational similarities. Those matching a selection threshold were then used to model a process of oligomerization. Remarkably, we showed a greater diversity in Aβ dimers than previously described. Depending on the clan family, different types of Aβ dimers were obtained. While some had the tendency to evolve into oligomeric rings, others formed fibrils of diverse characteristics. Then we selected the dimers that would evolve to membranephilic annular oligomers. Nearly one third of the 28 evaluated annular oligomers had the dimer interfaces between the neighboring Aβ42 monomers with possible salt bridges between the residue K28 from one side and either residue E22 or D23 on the other. Based on these results, key amino acids were identified for point mutations that either enhanced or suppressed the formation and toxicity of oligomer rings. Our studies suggest a greater diversity of Aβ dimers. Understanding the structure of Aβ dimers might be important for the rationale design of small molecules that block formation of toxic oligomers.

Biomimetic design of affinity peptide ligand for capsomere of virus-like particle

Virus-like particle (VLP) of murine polyomavirus (MPV) is a T = 7d icosahedral capsid that self-assembles from 72 capsomeres (Caps), each of which is a pentamer of major coat protein VP1. VLP has great potential in vaccinology, gene therapy, drug delivery, and materials science. However, its application is hindered by high cost downstream processes, leading to an urgent demand of a highly efficient affinity ligand for the separation and purification of Cap by affinity chromatography. Herein a biomimetic design strategy of an affinity peptide ligand of Cap has been developed on the basis of the binding structure of the C-terminus of minor coat protein (VP2-C) on the inner surface of Cap. The molecular interactions between VP2-C and Cap were first examined using all-atom molecular dynamics (MD) simulations coupled with the molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) method, where V283, P285, D286, W287, L289, and Y296 of VP2-C were identified as the hot spots. An affinity peptide library (DWXLXLXY, X denotes arbitrary amino acids except cysteine) was then constructed for virtual screening sequently by docking with AUTODOCK VINA, binding structure comparison, and final docking with ROSETTA FlexPepDock. Ten peptide candidates were selected and further confirmed by MD simulations and MM/PBSA, where DWDLRLLY was found to have the highest affinity to Cap. In DWDLRLLY, six residues are favorable for the binding, including W2, L4, L6 and Y8 inheriting from VP2-C, and R5 and L7 selected in the virtual screening. This confirms the high efficiency and accuracy of the biomimetic design strategy. DWDLRLLY was then experimentally validated by a one-step purification of Cap from crude cell lysate using affinity chromatography with the octapeptide immobilized on Sepharose gel. The purified Caps were observed to self-assemble into VLP with consistent structure of authentic MPV.

In silico investigation and targeting of amyloid β oligomers of different size

Aggregation of amyloid β (Aβ) peptides into fibrils has been implicated in the pathogenesis of Alzheimer's disease (AD). As a result, in recent years, substantial efforts have been expended in the study of the mechanism of aggregation of the Aβ peptide as well as of its inhibition by potential drug molecules. In this context, we have built a model of the Aβ(17-42) deca-oligomer using the solid-state NMR (ssNMR) structure of the Aβ(17-42) penta-oligomer as a reference. Both the penta- and deca-oligomer systems have been studied by all-atom molecular dynamics (MD) simulations and used as target systems for the investigation of the mechanism of action of a trehalose-derived Aβ aggregation inhibitor. In the deca-oligomer all the main structural features of the putative fibrillar state are retained. Moreover, the simulations reveal a remarkable gain in stability as the oligomer grows. MD studies of the inhibitor in complex with the penta- and deca-oligomers indicate a significant destabilization of the structure beyond the hampering of the addition of successive Aβ peptides at the ends of the fibril due to the presence of the inhibitor molecule. Our work provides an easy and effective approach which could be useful for the in silico development of potential drug molecules acting at different stages of the progression of Aβ-related diseases.

Molecular interactions of Alzheimer amyloid-β oligomers with neutral and negatively charged lipid bilayers

Interaction of p3 (Aβ(17-42)) peptides with cell membranes is crucial for the understanding of amyloid toxicity associated with Alzheimer's disease (AD). Such p3-membrane interactions are considered to induce the disruption of membrane permeability and integrity, but the exact mechanisms of how p3 aggregates, particularly small p3 oligomers, induce receptor-independent membrane disruption are not yet completely understood. Here, we investigate the adsorption, orientation, and surface interaction of the p3 pentamer with lipid bilayers composed of both pure zwitterionic POPC (palmitoyl-oleoyl-phosphatidylcholine) and mixed anionic POPC-POPG (palmitoyl-oleoyl-phosphatidylglycerol) (3 : 1) lipids using explicit-solvent molecular dynamics (MD) simulations. MD simulation results show that the p3 pentamer has much stronger interactions with mixed POPC-POPG lipids than pure POPC lipids, consistent with experimental observation that Aβ adsorption and fibrillation are enhanced on anionic lipid bilayers. Although electrostatic interactions are main attractive forces to drive the p3 pentamer to adsorb on the bilayer surface, the adsorption of the p3 pentamer on the lipid bilayer with C-terminal β-strands facing toward the bilayer surface is a net outcome of different competitions between p3 peptides-lipid bilayer and ions-p3-bilayer interactions. More importantly, Ca(2+) ions are found to form ionic bridges to associate negatively charged residues of p3 with anionic headgroups of the lipid bilayer, resulting in Aβ-Ca(2+)-PO4(-) complexes. Intensive Ca(2+) bound to the lipid bilayer and Ca(2+) ionic bridges may lead to Ca(2+) hemostasis responsible for neuronal dysfunction and death. This work provides insights into the mutual structure, dynamics, and interactions of both Aβ peptides and lipid bilayers at the atomic level, which expand our understanding of the complex behavior of amyloid-induced membrane disruption.

Aggregation and fibril morphology of the Arctic mutation of Alzheimer's Aβ peptide by CD, TEM, STEM and in situ AFM

Morphology of aggregation intermediates, polymorphism of amyloid fibrils and aggregation kinetics of the "Arctic" mutant of the Alzheimer's amyloid β-peptide, Aβ((1-40))(E22G), in a physiologically relevant Tris buffer (pH 7.4) were thoroughly explored in comparison with the human wild type Alzheimer's amyloid peptide, wt-Aβ((1-40)), using both in situ atomic force and electron microscopy, circular dichroism and thioflavin T fluorescence assays. For arc-Aβ((1-40)) at the end of the 'lag'-period of fibrillization an abrupt appearance of ≈ 3 nm size 'spherical aggregates' with a homogeneous morphology, was identified. Then, the aggregation proceeds with a rapid growth of amyloid fibrils with a variety of morphologies, while the spherical aggregates eventually disappeared during in situ measurements. Arc-Aβ((1-40)) was also shown to form fibrils at much lower concentrations than wt-Aβ((1-40)): ≤ 2.5 μM and 12.5 μM, respectively. Moreover, at the same concentration, 50 μM, the aggregation process proceeds more rapidly for arc-Aβ((1-40)): the first amyloid fibrils were observed after c.a. 72 h from the onset of incubation as compared to approximately 7 days for wt-Aβ((1-40)). Amyloid fibrils of arc-Aβ((1-40)) exhibit a large variety of polymorphs, at least five, both coiled and non-coiled distinct fibril structures were recognized by AFM, while at least four types of arc-Aβ((1-40)) fibrils were identified by TEM and STEM and their mass-per-length statistics were collected suggesting supramolecular structures with two, four and six β-sheet laminae. Our results suggest a pathway of fibrillogenesis for full-length Alzheimer's peptides with small and structurally ordered transient spherical aggregates as on-pathway immediate precursors of amyloid fibrils.

Distinct annular oligomers captured along the assembly and disassembly pathways of transthyretin amyloid protofibrils

Background

Defects in protein folding may lead to severe degenerative diseases characterized by the appearance of amyloid fibril deposits. Cytotoxicity in amyloidoses has been linked to poration of the cell membrane that may involve interactions with amyloid intermediates of annular shape. Although annular oligomers have been detected in many amyloidogenic systems, their universality, function and molecular mechanisms of appearance are debated.

Methodology/Principal Findings

We investigated with high-resolution in situ atomic force microscopy the assembly and disassembly of transthyretin (TTR) amyloid protofibrils formed of the native protein by pH shift. Annular oligomers were the first morphologically distinct intermediates observed in the TTR aggregation pathway. Morphological analysis suggests that they can assemble into a double-stack of octameric rings with a 16±2 nm diameter, and displaying the tendency to form linear structures. According to light scattering data coupled to AFM imaging, annular oligomers appeared to undergo a collapse type of structural transition into spheroid oligomers containing 8–16 monomers. Disassembly of TTR amyloid protofibrils also resulted in the rapid appearance of annular oligomers but with a morphology quite distinct from that observed in the assembly pathway.

Conclusions/Significance

Our observations indicate that annular oligomers are key dynamic intermediates not only in the assembly but also in the disassembly of TTR protofibrils. The balance between annular and more compact forms of aggregation could be relevant for cytotoxicity in amyloidogenic disorders.

Thermodynamic analysis of the molecular interactions between amyloid β-protein fragments and (-)-epigallocatechin-3-gallate

(-)-Epigallocatechin-3-gallate (EGCG) has been proven effective in preventing the aggregation of amyloid β-protein 42 (Aβ42), and the thermodynamic interactions between Aβ42 and EGCG have been studied in our previous work ( J. Phys. Chem. B 2010, 114, 11576). Herein, to further probe the interactions between different regions of Aβ42 and EGCG, three Aβ42 fragments (i.e., Aβ1-16, Aβ1-30, and Aβ31-42) were synthesized, and the thermodynamic interactions between each of the fragments and EGCG at different EGCG and salt concentrations were investigated by isothermal titration calorimetry. The results indicate that, although hydrogen bonding and hydrophobic interaction are both involved in the interactions between Aβ42 and EGCG, hydrogen bonding mainly happens in Aβ1-16 while hydrophobic interaction mainly happens in Aβ17-42. It is found that when Aβ42 and its fragments are saturated by EGCG, their thermodynamic parameters have linear relationships. The saturated binding stoichiometry (N(s)) for Aβ42 is the sum of the N(s) values for Aβ1-30 and Aβ31-42, while ΔH(s), ΔS(s), and ΔG(s) for Aβ42 are half the sum of the values for Aβ1-30 and Aβ31-42. The result suggests that there are no specific interactions and binding sites in the Aβ42 and EGCG binding. The orders of ΔH(s) and TΔS(s) values for the Aβ fragments are determined as Aβ17-42 > Aβ31-42 > Aβ1-30 > Aβ1-16. Moreover, there is significant enthalpy-entropy compensation in the binding of EGCG to Aβ42 and its fragments, resulting in insignificant change of ΔG with the change of the solution environment. The research has shed new light on the molecular mechanisms of the interactions between EGCG and Aβ42.

Structures of Aβ17-42 trimers in isolation and with five small-molecule drugs using a hierarchical computational procedure

The amyloid-β protein (Aβ) oligomers are believed to be the main culprits in the cytoxicity of Alzheimer's disease (AD) and p3 peptides (Aβ17-42 fragments) are present in AD amyloid plaques. Many small-molecule or peptide-based inhibitors are known to slow down Aβ aggregation and reduce the toxicity in vitro, but their exact modes of action remain to be determined since there has been no atomic level of Aβ(p3)-drug oligomers. In this study, we have determined the structure of Aβ17-42 trimers both in aqueous solution and in the presence of five small-molecule inhibitors using a multiscale computational study. These inhibitors include 2002-H20, curcumin, EGCG, Nqtrp, and resveratrol. First, we used replica exchange molecular dynamics simulations coupled to the coarse-grained (CG) OPEP force field. These CG simulations reveal that the conformational ensemble of Aβ17-42 trimer can be described by 14 clusters with each peptide essentially adopting turn/random coil configurations, although the most populated cluster is characterized by one peptide with a β-hairpin at Phe19-Leu31. Second, these 14 dominant clusters and the less-frequent fibril-like state with parallel register of the peptides were subjected to atomistic Autodock simulations. Our analysis reveals that the drugs have multiple binding modes with different binding affinities for trimeric Aβ17-42 although they interact preferentially with the CHC region (residues 17-21). The compounds 2002-H20 and Nqtrp are found to be the worst and best binders, respectively, suggesting that the drugs may interfere at different stages of Aβ oligomerization. Finally, explicit solvent molecular dynamics of two predicted Nqtrp-Aβ17-42 conformations describe at atomic level some possible modes of action for Nqtrp.

Introduction and technical survey: protein aggregation and fibrillogenesis

In this chapter we provided the overall background to the subject of protein aggregation and fibrillogenesis in amyloidogenesis, with introduction and brief discussion of the various topics that are included with the coming chapters. The division of the book into basic science and clinical science sections enables correlation of the topics to be made. The many proteins and peptides that have currently been found to undergo fibrillogenesis are tabulated. A broad technical survey is made, to indicate the vast array of techniques currently available to study aspects of protein oligomerization, aggregation and fibrillogenesis. These are split into three groups and tabulated, as the microscopical techniques, the analytical and biophysical methods, and the biochemical and cellular techniques. A few techniques are discussed, but in most cases only a link to relevant recent literature is provided.

Polymorphic structures of Alzheimer's β-amyloid globulomers

Background

Misfolding and self-assembly of Amyloid-β (Aβ) peptides into amyloid fibrils is pathologically linked to the development of Alzheimer's disease. Polymorphic Aβ structures derived from monomers to intermediate oligomers, protofilaments, and mature fibrils have been often observed in solution. Some aggregates are on-pathway species to amyloid fibrils, while the others are off-pathway species that do not evolve into amyloid fibrils. Both on-pathway and off-pathway species could be biologically relevant species. But, the lack of atomic-level structural information for these Aβ species leads to the difficulty in the understanding of their biological roles in amyloid toxicity and amyloid formation.

Methods and Findings

Here, we model a series of molecular structures of Aβ globulomers assembled by monomer and dimer building blocks using our peptide-packing program and explicit-solvent molecular dynamics (MD) simulations. Structural and energetic analysis shows that although Aβ globulomers could adopt different energetically favorable but structurally heterogeneous conformations in a rugged energy landscape, they are still preferentially organized by dynamic dimeric subunits with a hydrophobic core formed by the C-terminal residues independence of initial peptide packing and organization. Such structural organizations offer high structural stability by maximizing peptide-peptide association and optimizing peptide-water solvation. Moreover, curved surface, compact size, and less populated β-structure in Aβ globulomers make them difficult to convert into other high-order Aβ aggregates and fibrils with dominant β-structure, suggesting that they are likely to be off-pathway species to amyloid fibrils. These Aβ globulomers are compatible with experimental data in overall size, subunit organization, and molecular weight from AFM images and H/D amide exchange NMR.

Conclusions

Our computationally modeled Aβ globulomers provide useful insights into structure, dynamics, and polymorphic nature of Aβ globulomers which are completely different from Aβ fibrils, suggesting that these globulomers are likely off-pathway species and explaining the independence of the aggregation kinetics between Aβ globulomers and fibrils.

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.

Symmetry-based self-assembled nanotubes constructed using native protein structures: the key role of flexible linkers

We construct nanotubes using native protein structures and their native associations from structural databases. The construction is based on a shape-guided symmetric self-assembly concept. Our strategy involves fusing judiciously-selected oligomerization domains via peptide linkers. Linkers are inherently flexible, hence their choice is critical: they should position the domains in three-dimensional space in the desired orientation while retaining their own natural conformational tendencies; however, at the same time, retain the construct stability. Here we outline a design scheme which accounts for linker flexibility considerations, and present two examples. The first is HIV-1 capsid protein, which in vitro self-assembles into nanotubes and conical capsids, and its linker exists as a short flexible loop. The second involves novel nanotubes construction based on antimicrobial homodimer Magainin 2, employing linkers of distinct lengths and flexibility levels. Our strategy utilizes the abundance of unique shapes and sizes of proteins and their building blocks which can assemble into a vast number of combinations, and consequently, nanotubes of distinct morphologies and diameters. Computational design and assessment methodologies can help reduce the number of candidates for experimental validation. This is an invited paper for a special issue on protein dynamics, here focusing on flexibility in nanotube design based on protein building blocks.

Binding of nonsteroidal anti-inflammatory drugs to Abeta fibril

Nonsteroidal anti-inflammatory drugs are considered as potential therapeutic agents against Alzheimer's disease. Using replica exchange molecular dynamics and atomistic implicit solvent model, we studied the mechanisms of binding of naproxen and ibuprofen to the Abeta fibril derived from solid-state NMR measurements. The binding temperature of naproxen is found to be almost 40 K higher than of ibuprofen implicating higher binding affinity of naproxen. The key factor, which enhances naproxen binding, is strong interactions between ligands bound to the surface of the fibril. The naphthalene ring in naproxen appears to provide a dominant contribution to ligand-ligand interactions. In contrast, ligand-fibril interactions cannot explain differences in the binding affinities of naproxen and ibuprofen. The concave fibril edge with the groove is identified as the primary binding location for both ligands. We show that confinement of the ligands to the groove facilitates ligand-ligand interactions that lowers the energy of the ligands bound to the concave edge compared with those bound to the convex edge. Our simulations appear to provide microscopic rationale for the differing binding affinities of naproxen and ibuprofen observed experimentally.

An oligomeric equilibrium intermediate as the precursory nucleus of globular and fibrillar supramacromolecular assemblies in a PDZ domain

The equilibrium unfolding at neutral pH of the third PDZ domain of PSD95, as followed by DSC, is characterized by the presence of an equilibrium intermediate with clear signs of oligomerization. DLS and SEC measurements indicate that at 60-70 degrees C small oligomers populate, showing a typical beta-sheet far-UV CD spectrum. These intermediate species lead to the formation of rodlike particulates of approximately 12 nm, which remain in solution after 2 weeks incubation and grow until they adopt annular/spherical shapes of approximately 50 nm and protofibrils, which are subsequently fully transformed into fibrils. The fibrils can also disaggregate after the addition of 1:1 buffer dilution followed by cooling to room temperature, thus returning to the initial monomeric state. Growth kinetics, as shown by ThT and ANS fluorescence, show that the organization of the different supramacromolecular structures comes from a common nucleation unit, the small oligomers, which organize themselves before reaching the incubation temperature of 60 degrees C. Our experiments point toward the existence of a well-defined reversible, stepwise, and downhill organization of the processes involved in the association-dissociation of the intermediate. We estimate the enthalpy change accompanying the association-dissociation equilibria to be 130 kJ x mol(-1). Furthermore, the coalescence under essentially reversible conditions of different kinds of supramacromolecular assemblies renders this protein system highly interesting for biophysical studies aimed at our further understanding of amyloid pathological conditions.

Structural determination of Abeta25-35 micelles by molecular dynamics simulations

Amyloid-beta (Abeta) peptides and other amyloidogenic proteins can form a wide range of soluble oligomers of varied morphologies at the early aggregation stage, and some of these oligomers are biologically relevant to the pathogenesis of Alzheimer's disease. Spherical micelle-like oligomers have been often observed for many different types of amyloids. Here, we report a hybrid computational approach to systematically construct, search, optimize, and rank soluble micelle-like Abeta25-35 structures with different side-chain packings at the atomic level. Simulations reveal for the first time, to our knowledge, that two Abeta micelles with antiparallel peptide organization and distinct surface hydrophobicity display high structural stability. Stable micelles experience a slow secondary structural transition from turn to alpha-helix. Energetic analysis coupled with computational mutagenesis reveals that van der Waals and solvation energies play a more pronounced role in stabilizing the micelles, whereas the electrostatic energies present a stable but minor energetic contribution to peptide assemblies. Modeled Abeta micelles with shapes and dimensions similar to those of experimentally derived spherical structures also provide detailed information about the roles of structural dynamics and transition in the formation of amyloid fibrils. The strong binding affinity of our micelles to antibodies implies that micelles may be a biologically relevant species.

Structural convergence among diverse, toxic beta-sheet ion channels

Recent studies show that an array of beta-sheet peptides, including N-terminally truncated Abeta peptides (Abeta(11-42/17-42)), K3 (a beta(2)-microglobulin fragment), and protegrin-1 (PG-1) peptides form ion channel-like structures and elicit single channel ion conductance when reconstituted in lipid bilayers and induce cell damage through cell calcium overload. Striking similarities are observed in the dimensions of these toxic channels irrespective of their amino acid sequences. However, the intriguing question of preferred channel sizes is still unresolved. Here, exploiting ssNMR-based, U-shaped, beta-strand-turn-beta-strand coordinates, we modeled truncated Abeta peptide (p3) channels with different sizes (12- to 36-mer). Molecular dynamics (MD) simulations show that optimal channel sizes of the ion channels presenting toxic ionic flux range between 16- and 24-mer. This observation is in good agreement with channel dimensions imaged by AFM for Abeta(9-42), K3 fragment, and PG-1 channels and highlights the bilayer-supported preferred toxic beta-channel sizes and organization, regardless of the peptide sequence.

Structures of beta-amyloid peptide 1-40, 1-42, and 1-55-the 672-726 fragment of APP-in a membrane environment with implications for interactions with gamma-secretase

Aggregation of Amyloid beta (Abeta) peptide has been linked to the neurodegenerative Alzheimer's Disease and implicated in other amyloid diseases including cerebral amyloid angiopathy. Abeta peptide is generated by cleavage of the amyloid precursor protein (APP) by transmembrane proteases. It is crucial to determine the structures of beta-amyloid peptides in a membrane to provide a molecular basis for the cleavage mechanism. We report the structures of amyloid beta peptide (Abeta(1-40) and Abeta(1-42)) as well as the 672-726 fragment of APP (referred to as Abeta(1-55)) in a membrane environment determined by replica-exchange molecular dynamics simulation. Abeta(1-40) is found to have two helical domains A (13-22) and B(30-35) and a type I beta-turn at 23-27. The peptide is localized at the interface between membrane and solvent. Substantial fluctuations in domain A are observed. The dominant simulated tertiary structure of Abeta(1-40) is observed to be similar to the simulated Abeta(1-42) structure. However, there are differences observed in the overall conformational ensemble, as characterized by the two-dimensional free energy surfaces. The fragment of APP (Abeta(1-55)) is observed to have a long transmembrane helix. The position of the transmembrane region and ensemble of membrane structures are elucidated. The conformational transition between the transmembrane Abeta(1-55) structure, prior to cleavage, and the Abeta(1-40) structure, following cleavage, is proposed.

Molecular modeling of two distinct triangular oligomers in amyloid beta-protein

Amyloid-beta (Abeta) peptides exhibit many distinct structural morphology at the early aggregate stage, some of which are biological relevant to the pathogenesis of Alzheimer's disease (AD). Atomic-resolution structures of the early Abeta aggregates and their conformational changes in amyloid aggregation remain elusive. Here, we perform all-atom molecular modeling and dynamics simulations to obtain two stable triangular-like Abeta structures with the lowest packing energy, one corresponding to the Tycko's model (Paravastu, A.; Leapman, R.; Yau, W.; Tycko, R. Proc. Nat. Acad. Soc. U.S.A. 2008, 105, 18349-18354) (referred to C-WT model) and the other corresponding to computational model (N-WT model). Both models have the same 3-fold symmetry but distinct beta-sheet organizations in which three Abeta hexamers pack together via either C-terminal beta-strand residues or N-terminal beta-strand residues forming distinct hydrophobic cross section. Structural and energetic comparisons of two 3-fold Abeta oligomers, coupled with structural changes upon the mutations occurring at the interacting interfaces, reveal that although hydrophobic interactions are still dominant forces, electrostatic interactions are more favorable in the N-WT model due to the formation of more and stable intersheet salt bridges, while solvation energy is more favorable in the C-WT model due to more exposed hydrophilic residues to solvent. Both models display many common features similar to other amyloid oligomers and therefore are likely to be biologically relevant.

Misfolded amyloid ion channels present mobile beta-sheet subunits in contrast to conventional ion channels

In Alzheimer's disease, calcium permeability through cellular membranes appears to underlie neuronal cell death. It is increasingly accepted that calcium permeability involves toxic ion channels. We modeled Alzheimer's disease ion channels of different sizes (12-mer to 36-mer) in the lipid bilayer using molecular dynamics simulations. Our Abeta channels consist of the solid-state NMR-based U-shaped beta-strand-turn-beta-strand motif. In the simulations we obtain ion-permeable channels whose subunit morphologies and shapes are consistent with electron microscopy/atomic force microscopy. In agreement with imaged channels, the simulations indicate that beta-sheet channels break into loosely associated mobile beta-sheet subunits. The preferred channel sizes (16- to 24-mer) are compatible with electron microscopy/atomic force microscopy-derived dimensions. Mobile subunits were also observed for beta-sheet channels formed by cytolytic PG-1 beta-hairpins. The emerging picture from our large-scale simulations is that toxic ion channels formed by beta-sheets spontaneously break into loosely interacting dynamic units that associate and dissociate leading to toxic ionic flux. This sharply contrasts intact conventional gated ion channels that consist of tightly interacting alpha-helices that robustly prevent ion leakage, rather than hydrogen-bonded beta-strands. The simulations suggest why conventional gated channels evolved to consist of interacting alpha-helices rather than hydrogen-bonded beta-strands that tend to break in fluidic bilayers. Nature designs folded channels but not misfolded toxic channels.

Dynamics and stability of amyloid-like steric zipper assemblies with hydrophobic dry interfaces

Recent seminal investigations have suggested that the basic structural motif of amyloid fibers may be constituted by a tight association of two facing beta-sheets (steric zipper). Although this model has been derived from crystal structures of small peptide models, several theoretical investigations, essentially focused on steric zipper interface containing large polar and/or aromatic side chains, have confirmed the stability of this motif in a crystal-free context. To analyze the general validity of these findings, we carried out molecular dynamics (MD) simulations on aggregates stabilized by steric zipper interfaces made also of small or hydrophobic residues. In particular, we here characterized assemblies formed by the peptides SSTSAA and VQIVYK, whose structures have been recently solved at high resolution. In contrast to previous results obtained for polar/aromatic aggregates of the same size and with similar interface area, steric zipper assemblies composed of a pair of 10-stranded beta-sheets show high fluctuations and significant distortions in the simulation timescales (40-60 ns). Taking into account the crystal packing, the effect of the addition of an extra sheet to the assemblies was also evaluated. The MD results indicate that this addition does not provide extra-stabilization to the pair of sheet models. Although present data do not preclude the possibility that the steric zipper association identified in the crystal structure is the basic motif of SSTSAA and VQIVYK fibers, our findings highlight the importance of the nature of residues directly involved in the motif. Indeed, polar and aromatic residues that may form intrasheet and intersheet interactions likely provide a strong contribution to the steric zipper motif stability. Along this line, assemblies endowed with hydrophobic residues presumably require larger interfaces. In line with this suggestion, MD analysis of the HET-s(218-289) prion models composed of a similar number of strands shows that the assembly is endowed with a remarkable stability.