Thermodynamic description of polymorphism in Q- and N-rich peptide aggregates revealed by atomistic simulation

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.

Identifying Biological and Biophysical Features of Different Maturation States of α-Synuclein Amyloid Fibrils

Protein aggregates, hereunder amyloid fibrils, can undergo a maturation process, whereby early formed aggregates undergo a structural and physicochemical transition leading to more mature species. In the case of amyloid-related diseases, such maturation confers distinctive biological properties of the aggregates, which may account for a range of diverse pathological subtypes. Here, we present a protocol for the preparation of α-synuclein amyloid fibrils differing in the level of their maturation. We utilize widely accessible biophysical techniques to characterize the structure and morphology and a simple thermal treatment procedure to test their thermodynamic stability. Their biological properties are probed by means of binding to native plasma membrane sheets originating from mammalian cell lines.

A focused review on three-dimensional bioprinting technology for artificial organ fabrication

Three-dimensional (3D) bioprinting technology has attracted a great deal of interest because it can be easily adapted to many industries and research sectors, such as biomedical, manufacturing, education, and engineering. Specifically, 3D bioprinting has provided significant advances in the medical industry, since such technology has led to significant breakthroughs in the synthesis of biomaterials, cells, and accompanying elements to produce composite living tissues. 3D bioprinting technology could lead to the immense capability of replacing damaged or injured tissues or organs with newly dispensed cell biomaterials and functional tissues. Several types of bioprinting technology and different bio-inks can be used to replicate cells and generate supporting units as complex 3D living tissues. Bioprinting techniques have undergone great advancements in the field of regenerative medicine to provide 3D printed models for numerous artificial organs and transplantable tissues. This review paper aims to provide an overview of 3D-bioprinting technologies by elucidating the current advancements, recent progress, opportunities, and applications in this field. It highlights the most recent advancements in 3D-bioprinting technology, particularly in the area of artificial organ development and cancer research. Additionally, the paper speculates on the future progress in 3D-bioprinting as a versatile foundation for several biomedical applications.

Quantification of amyloid fibril polymorphism by nano-morphometry reveals the individuality of filament assembly

Amyloid fibrils are highly polymorphic structures formed by many different proteins. They provide biological function but also abnormally accumulate in numerous human diseases. The physicochemical principles of amyloid polymorphism are not understood due to lack of structural insights at the single-fibril level. To identify and classify different fibril polymorphs and to quantify the level of heterogeneity is essential to decipher the precise links between amyloid structures and their functional and disease associated properties such as toxicity, strains, propagation and spreading. Employing gentle, force-distance curve-based AFM, we produce detailed images, from which the 3D reconstruction of individual filaments in heterogeneous amyloid samples is achieved. Distinctive fibril polymorphs are then classified by hierarchical clustering, and sample heterogeneity is objectively quantified. These data demonstrate the polymorphic nature of fibril populations, provide important information regarding the energy landscape of amyloid self-assembly, and offer quantitative insights into the structural basis of polymorphism in amyloid populations.

Amyloid Evolution: Antiparallel Replaced by Parallel

Several atomic structures have now been found for micrometer-scale amyloid fibrils or elongated microcrystals using a range of methods, including NMR, electron microscopy, and X-ray crystallography, with parallel β-sheet appearing as the most common secondary structure. The etiology of amyloid disease, however, indicates nanometer-scale assemblies of only tens of peptides as significant agents of cytotoxicity and contagion. By combining solution X-ray with molecular dynamics, we show that antiparallel structure dominates at the first stages of aggregation for a specific set of peptides, being replaced by parallel at large length scales only. This divergence in structure between small and large amyloid aggregates should inform future design of molecular therapeutics against nucleation or intercellular transmission of amyloid. Calculations and an overview from the literature argue that antiparallel order should be the first appearance of structure in many or most amyloid aggregation processes, regardless of the endpoint. Exceptions to this finding should exist, depending inevitably on the sequence and on solution conditions.

Early Stage Alpha-Synuclein Amyloid Fibrils are Reservoirs of Membrane-Binding Species

The presence of αSN fibrils indisputably associates with the development of synucleinopathies. However, while certain fibril morphologies have been linked to downstream pathological phenotypes, others appear less harmful, leading to the concept of fibril strains, originally described in relation to prion disease. Indeed, the presence of fibrils does not associate directly with neurotoxicity. Rather, it has been suggested that the toxic compounds are soluble amyloidogenic oligomers, potentially co-existing with fibrils. Here, combining synchrotron radiation circular dichroism, transmission electron microscopy and binding assays on native plasma membrane sheets, we reveal distinct biological and biophysical differences between initial and matured fibrils, transformed within the timespan of few days. Immature fibrils are reservoirs of membrane-binding species, which in response to even gentle experimental changes release into solution in a reversible manner. In contrast, mature fibrils, albeit macroscopically indistinguishable from their less mature counterparts, are structurally robust, shielding the solution from the membrane active soluble species. We thus show that particular biological activity resides transiently with the fibrillating sample, distinct for one, but not the other, spontaneously formed fibril polymorph. These results shed new light on the principles of fibril polymorphism with consequent impact on future design of assays and therapeutic development.

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.

Peptide strand length controls the energetics of self-assembly and morphology of β-sheet fibrils

Self-assembling peptides can be used as versatile, natural, and multifunctional building blocks to produce a variety of well-defined nanostructures, materials and devices for applications in medicine and nanotechnology. Here, we concentrate on the 1D self-assembly of de novo designed Px-2 peptide β-strands into anti-parallel β-sheet tapes and higher order aggregates. We study six members of the Px-2 family, ranging from 3 amino acids (aa) to 13 aa in length, using a range of complementary experimental techniques, computer simulation and theoretical statistical mechanics. The critical concentration for self-assembly (c*) is found to increase systematically with decreasing peptide length. The shortest peptide found to self-assemble into soluble β-tapes in water is a 5 amino acid residue peptide. These investigations help decipher the role of the peptide length in controlling self-assembly, aggregate morphology, and material properties. By extracting free energies from these data using a statistical mechanical analysis and combining the results with computer simulations at the atomistic level, we can extract the entropy of association for individual β-strands.

Self-consistent implementation of meta-GGA functionals for the ONETEP linear-scaling electronic structure package

Accurate and computationally efficient exchange-correlation functionals are critical to the successful application of linear-scaling density functional theory (DFT). Local and semi-local functionals of the density are naturally compatible with linear-scaling approaches, having a general form which assumes the locality of electronic interactions and which can be efficiently evaluated by numerical quadrature. Presently, the most sophisticated and flexible semi-local functionals are members of the meta-generalized-gradient approximation (meta-GGA) family, and depend upon the kinetic energy density, τ, in addition to the charge density and its gradient. In order to extend the theoretical and computational advantages of τ-dependent meta-GGA functionals to large-scale DFT calculations on thousands of atoms, we have implemented support for τ-dependent meta-GGA functionals in the ONETEP program. In this paper we lay out the theoretical innovations necessary to implement τ-dependent meta-GGA functionals within ONETEP's linear-scaling formalism. We present expressions for the gradient of the τ-dependent exchange-correlation energy, necessary for direct energy minimization. We also derive the forms of the τ-dependent exchange-correlation potential and kinetic energy density in terms of the strictly localized, self-consistently optimized orbitals used by ONETEP. To validate the numerical accuracy of our self-consistent meta-GGA implementation, we performed calculations using the B97M-V and PKZB meta-GGAs on a variety of small molecules. Using only a minimal basis set of self-consistently optimized local orbitals, we obtain energies in excellent agreement with large basis set calculations performed using other codes. Finally, to establish the linear-scaling computational cost and applicability of our approach to large-scale calculations, we present the outcome of self-consistent meta-GGA calculations on amyloid fibrils of increasing size, up to tens of thousands of atoms.

Applications of large-scale density functional theory in biology

Density functional theory (DFT) has become a routine tool for the computation of electronic structure in the physics, materials and chemistry fields. Yet the application of traditional DFT to problems in the biological sciences is hindered, to a large extent, by the unfavourable scaling of the computational effort with system size. Here, we review some of the major software and functionality advances that enable insightful electronic structure calculations to be performed on systems comprising many thousands of atoms. We describe some of the early applications of large-scale DFT to the computation of the electronic properties and structure of biomolecules, as well as to paradigmatic problems in enzymology, metalloproteins, photosynthesis and computer-aided drug design. With this review, we hope to demonstrate that first principles modelling of biological structure-function relationships are approaching a reality.

Assembly of Triblock Amphiphilic Peptides into One-Dimensional Aggregates and Network Formation

Peptide assembly plays a key role in both neurological diseases and development of novel biomaterials with well-defined nanostructures. Synthetic model peptides provide a unique platform to explore the role of intermolecular interactions in the assembly process. A triblock peptide architecture designed by the Hartgerink group is a versatile system which relies on Coulomb interactions, hydrogen bonding, and hydrophobicity to guide these peptides' assembly at three different length scales: β-sheets, double-wall ribbon-like aggregates, and finally a highly porous network structure which can support gels with ≤1% by weight peptide concentration. In this study, by using molecular dynamics simulations of a structure based implicit solvent coarse grained model, we analyzed this hierarchical assembly process. Parametrization of our CG model is based on multiple-state points from atomistic simulations, which enables this model to represent the conformational adaptability of the triblock peptide molecule based on the surrounding medium. Our results indicate that emergence of the double-wall β-sheet packing mechanism, proposed in light of the experimental evidence, strongly depends on the subtle balance of the intermolecular forces. We demonstrate that, even though backbone hydrogen bonding dominates the early nucleation stages, depending on the strength of the hydrophobic and Coulomb forces, alternative structures such as zero-dimensional aggregates with two β-sheets oriented orthogonally (which we refer to as a cross-packed structure) and β-sheets with misoriented hydrophobic side chains are also feasible. We discuss the implications of these competing structures for the three different length scales of assembly by systematically investigating the influence of density, counterion valency, and hydrophobicity.

Polyglutamine Fibrils: New Insights into Antiparallel β-Sheet Conformational Preference and Side Chain Structure

Understanding the structure of polyglutamine (polyQ) amyloid-like fibril aggregates is crucial to gaining insights into the etiology of at least ten neurodegenerative disorders, including Huntington's disease. Here, we determine the structure of D2Q10K2 (Q10) fibrils using ultraviolet resonance Raman (UVRR) spectroscopy and molecular dynamics (MD). Using UVRR, we determine the fibril peptide backbone Ψ and glutamine (Gln) side chain χ3 dihedral angles. We find that most of the fibril peptide bonds adopt antiparallel β-sheet conformations; however, a small population of peptide bonds exist in parallel β-sheet structures. Using MD, we simulate three different potential fibril structural models that consist of either β-strands or β-hairpins. Comparing the experimentally measured Ψ and χ3 angle distributions to those obtained from the MD simulated models, we conclude that the basic structural motif of Q10 fibrils is an extended β-strand structure. Importantly, we determine from our MD simulations that Q10 fibril antiparallel β-sheets are thermodynamically more stable than parallel β-sheets. This accounts for why polyQ fibrils preferentially adopt antiparallel β-sheet conformations instead of in-register parallel β-sheets like most amyloidogenic peptides. In addition, we directly determine, for the first time, the structures of Gln side chains. Our structural data give new insights into the role that the Gln side chains play in the stabilization of polyQ fibrils. Finally, our work demonstrates the synergistic power and utility of combining UVRR measurements and MD modeling to determine the structure of amyloid-like fibrils.

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.

Orientations of residues along the β-arch of self-assembled amylin fibril-like structures lead to polymorphism

Amylin is an endocrine hormone peptide that consists of 37 residues and is the main component of extracellular amyloid deposits found in the pancreas of most type 2 diabetes patients. Amylin peptides are self-assembled to form oligomers and fibrils. So far, four different molecular structures of the self-assembled amylin fibrils have been observed experimentally: two ssNMR models and two crystal models. This study reveals, for the first time, that there are four self-assembled amylin forms that differ in the orientations of the side chains along the β-arch and are all derived from the two ssNMR models. The two ssNMR models are composed of these four different self-assembled forms of amylin, and the two crystal models are composed of two different self-assembled forms of amylin. This study illustrates at the atomic level the differences among the four experimental models and proposes eight new models of self-assembled amylin that are also composed of the four different self-assembled forms of amylin. Our results show polymorphism of the self-assembled fibril-like amylin, with a slight preference of some of the newly constructed models over the experimental models. Finally, we propose that two different self-assembled fibril-like forms of amylin can interact to form a new fibril-like amylin. We investigated this argument and found that some fibril-like amylin prefers to interact to form stable fibril-like structures, whereas others disfavor it. Our work provides new insights that may suggest strategies for future pharmacological studies that aim to find ways to ameliorate the interactions between polymorphic oligomers and fibrils of amylin.

DNAJB6 is a peptide-binding chaperone which can suppress amyloid fibrillation of polyglutamine peptides at substoichiometric molar ratios

Expanded polyglutamine (polyQ) stretches lead to protein aggregation and severe neurodegenerative diseases. A highly efficient suppressor of polyQ aggregation was identified, the DNAJB6, when molecular chaperones from the HSPH, HSPA, and DNAJ families were screened for huntingtin exon 1 aggregation in cells (Hageman et al. in Mol Cell 37(3):355-369, 2010). Furthermore, also aggregation of polyQ peptides expressed in cells was recently found to be efficiently suppressed by co-expression of DNAJB6 (Gillis et al. in J Biol Chem 288:17225-17237, 2013). These suppression effects can be due to an indirect effect of DNAJB6 on other cellular components or to a direct interaction between DNAJB6 and polyQ peptides that may depend on other cellular components. Here, we have purified the DNAJB6 protein to investigate the suppression mechanism. The purified DNAJB6 protein formed large heterogeneous oligomers, in contrast to the more canonical family member DNAJB1 which is dimeric. Purified DNAJB6 protein, at substoichiometric molar ratios, efficiently suppressed fibrillation of polyQ peptides with 45°Q in a thioflavin T fibrillation. No suppression was obtained with DNAJB1, but with the closest homologue to DNAJB6, DNAJB8. The suppression effect was independent of HSPA1 and ATP. These data, based on purified proteins and controlled fibrillation in vitro, strongly suggest that the fibrillation suppression is due to a direct protein-protein interaction between the polyQ peptides and DNAJB6 and that the DNAJB6 has unique fibrillation suppression properties lacking in DNAJB1. Together, the data obtained in cells and in vitro support the view that DNAJB6 is a peptide-binding chaperone that can interact with polyQ peptides that are incompletely degraded by and released from the proteasome.

Role of sequence and structural polymorphism on the mechanical properties of amyloid fibrils

Amyloid fibrils playing a critical role in disease expression, have recently been found to exhibit the excellent mechanical properties such as elastic modulus in the order of 10 GPa, which is comparable to that of other mechanical proteins such as microtubule, actin filament, and spider silk. These remarkable mechanical properties of amyloid fibrils are correlated with their functional role in disease expression. This suggests the importance in understanding how these excellent mechanical properties are originated through self-assembly process that may depend on the amino acid sequence. However, the sequence-structure-property relationship of amyloid fibrils has not been fully understood yet. In this work, we characterize the mechanical properties of human islet amyloid polypeptide (hIAPP) fibrils with respect to their molecular structures as well as their amino acid sequence by using all-atom explicit water molecular dynamics (MD) simulation. The simulation result suggests that the remarkable bending rigidity of amyloid fibrils can be achieved through a specific self-aggregation pattern such as antiparallel stacking of β strands (peptide chain). Moreover, we have shown that a single point mutation of hIAPP chain constituting a hIAPP fibril significantly affects the thermodynamic stability of hIAPP fibril formed by parallel stacking of peptide chain, and that a single point mutation results in a significant change in the bending rigidity of hIAPP fibrils formed by antiparallel stacking of β strands. This clearly elucidates the role of amino acid sequence on not only the equilibrium conformations of amyloid fibrils but also their mechanical properties. Our study sheds light on sequence-structure-property relationships of amyloid fibrils, which suggests that the mechanical properties of amyloid fibrils are encoded in their sequence-dependent molecular architecture.

A kinetic approach to the sequence-aggregation relationship in disease-related protein assembly

It is generally accepted that oligomers of aggregating proteins play an important role in the onset of neurodegenerative diseases. While in silico aggregation studies of full length amyloidogenic proteins are computationally expensive, the assembly of short protein fragments derived from these proteins with similar aggregating properties has been extensively studied. In the present work, molecular dynamics simulations are performed to follow peptide aggregation on the microsecond time scale. By defining aggregation states, we identify transition networks, disconnectivity graphs, and first passage time distributions to describe the kinetics of the assembly process. This approach unravels differences in the aggregation into hexamers of two peptides with different primary structures. The first is GNNQQNY, a hydrophilic fragment from the prion protein Sup35, and the second is KLVFFAE, a fragment from amyloid-β protein, with a hydrophobic core delimited by two charged amino acids. The assembly of GNNQQNY suggests a mechanism of monomer addition, with a bias toward parallel peptide pairs and a gradual increase in the amount of β-strand content. For KLVFFAE, a mechanism involving dimers rather than monomers is revealed, involving a generally higher β-strand content and a transition toward a larger number of antiparallel peptide pairs during the rearrangement of the hexamer. The differences observed for the aggregation of the two peptides suggests the existence of a sequence-aggregation relationship.

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.

Influence of temperature on formation of perfect tau fragment fibrils using PRIME20/DMD simulations

We investigate the fibrillization process for amyloid tau fragment peptides (VQIVYK) by applying the discontinuous molecular dynamics method to a system of 48 VQIVYK peptides modeled using a new protein model/force field, PRIME20. The aim of the article is to ascertain which factors are most important in determining whether or not a peptide system forms perfect coherent fibrillar structures. Two different directional criteria are used to determine when a hydrogen bond occurs: the original H-bond constraints and a parallel preference H-bond constraint that imparts a slight bias towards the formation of parallel versus antiparallel strands in a β-sheet. Under the original H-bond constraints, the resulting fibrillar structures contain a mixture of parallel and antiparallel pairs of strands within each β-sheet over the whole fibrillization temperature range. Under the parallel preference H-bond constraints, the β-sheets within the fibrillar structures are more likely to be parallel and indeed become perfectly parallel, consistent with X-ray crystallography, at a high temperature slightly below the fibrillization temperature. The high temperature environment encourages the formation of perfect fibril structures by providing enough time and space for peptides to rearrange during the aggregation process. There are two different kinetic mechanisms, template assembly with monomer addition at high temperature and merging/rearrangement without monomer addition at low temperature, which lead to significant differences in the final fibrillar structure. This suggests that the diverse fibril morphologies generally observed in vitro depend more on environmental conditions than has heretofore been appreciated.

Self-assembly of amphipathic β-sheet peptides: insights and applications

Amphipathic peptides composed of alternating polar and nonpolar residues have a strong tendency to self-assemble into one-dimensional, amyloid-like fibril structures. Fibrils derived from peptides of general (XZXZ)(n) sequence in which X is hydrophobic and Z is hydrophilic adopt a putative β-sheet bilayer. The bilayer configuration allows burial of the hydrophobic X side chain groups in the core of the fibril and leaves the polar Z side chains exposed to solvent. This architectural arrangement provides fibrils that maintain high solubility in water and has facilitated the recent exploitation of self-assembled amphipathic peptide fibrils as functional biomaterials. This article is a critical review of the development and application of self-assembling amphipathic peptides with a focus on the fundamental insight these types of peptides provide into peptide self-assembly phenomena.

Protein aggregation: mechanisms and functional consequences

Understanding the mechanisms underlying protein misfolding and aggregation has become a central issue in biology and medicine. Compelling evidence show that the formation of amyloid aggregates has a negative impact in cell function and is behind the most prevalent human degenerative disorders, including Alzheimer's Parkinson's and Huntington's diseases or type 2 diabetes. Surprisingly, the same type of macromolecular assembly is used for specialized functions by different organisms, from bacteria to human. Here we address the conformational properties of these aggregates, their formation pathways, their role in human diseases, their functional properties and how bioinformatics tools might be of help to study these protein assemblies.

Spontaneous formation of twisted Aβ(16-22) fibrils in large-scale molecular-dynamics simulations

Protein aggregation is associated with fatal neurodegenerative diseases, including Alzheimer's and Parkinson's. Mapping out kinetics along the aggregation pathway could provide valuable insights into the mechanisms that drive oligomerization and fibrillization, but that is beyond the current scope of computational research. Here we trace out the full kinetics of the spontaneous formation of fibrils by 48 Aβ(16-22) peptides, following the trajectories in molecular detail from an initial random configuration to a final configuration of twisted protofilaments with cross-β-structure. We accomplish this by performing large-scale molecular-dynamics simulations based on an implicit-solvent, intermediate-resolution protein model, PRIME20. Structural details such as the intersheet distance, perfectly antiparallel β-strands, and interdigitating side chains analogous to a steric zipper interface are explained by and in agreement with experiment. Two characteristic fibrillization mechanisms - nucleation/templated growth and oligomeric merging/structural rearrangement - emerge depending on the temperature.

Alternative packing modes leading to amyloid polymorphism in five fragments studied with molecular dynamics

Amyloid aggregates have been implicated in the pathogenesis of diseases such as type 2 diabetes, Alzheimer's, Parkinson's, and prion disease. Recently determined microcrystal structures of several short peptide segments derived from fibril-forming proteins revealed coexistence of alternative aggregation modes (amyloid polymorphism) formed by the same segment. This polymorphism may help in understanding the influence of the side chain packing on the amyloid stability. Here we use molecular dynamics (MD) simulation to analyze the stability of five pairs of polar and nonpolar polymorphic oligomers. MD simulation shows polymorphs with steric zipper interface containing large polar and/or aromatic side chains (GNNQQNY, and NNQNTF) are more stable than steric zipper interfaces made of small or hydrophobic residues (SSTNGVG, VQIVYK, and MVGGVV). Several geometric analyses revealed that larger sheet to sheet interface of the dry steric zipper through polar Q/N rich side chains holds the sheets together. Mutant simulations (Q/N→G) show substitutions with glycine disrupt the steric zipper, leading to unstable oligomers. Stability of Q/N rich oligomers was found to result from the large average number of hydrogen bonds. The molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) method reports the nonpolar component of free energy to be favorable, while electrostatic solvation is unfavorable for β-sheet association. Knowledge of structural properties of these fibrils might be useful for developing therapeutic agents against amyloidoses as well as for developing biomaterials. © 2011 Wiley Periodicals, Inc. Biopolymers (Pept Sci) 98: 131-144, 2012.

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.

Systematic examination of polymorphism in amyloid fibrils by molecular-dynamics simulation

Amyloid fibrils often exhibit polymorphism. Polymorphs are formed when proteins or peptides with identical sequences self-assemble into fibrils containing substantially different arrangements of the β-strands. We used atomistic molecular-dynamics simulation to examine the thermodynamic stability of a amyloid fibrils in different polymorphic forms by performing a systematic investigation of sequence and symmetry space for a series of peptides with a range of physicochemical properties. We show that the stability of fibrils depends on both sequence and the symmetry because these factors determine the availability of favorable interactions between the peptide strands within a sheet and in intersheet packing. By performing a detailed analysis of these interactions as a function of symmetry, we obtained a series of simple design rules that can be used to determine which polymorphs of a given sequence are most likely to form thermodynamically stable fibrils. These rules can potentially be employed to design peptide sequences that aggregate into a preferred polymorphic form for nanotechnological purposes.

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.

Molecular dynamic simulation of wild type and mutants of the polymorphic amyloid NNQNTF segments of elk prion: structural stability and thermodynamic of association

A hexapeptide with amino acid sequence NNQNTF from the elk prion protein forms amyloid fibrils. Here we use molecular dynamic simulations of the oligomers and their single point glycine mutants to study their stability. In an effort to probe the structural stability and association thermodynamic in a realistic environment, all wildtype of NNQNTF polymorphic forms with different size and their corresponding double layer 5 strands single point glycine mutants were subjected to a total of 500 ns of explicit-solvent molecular dynamics (MD) simulation. Our results show that the structural stability of the NNQNTF oligomers increases with increasing the number of β-strands for double layers. Our results also demonstrated that hydrophobic interaction is the principle driving force to stabilize the adjacent β-strands while the steric zipper is responsible for holding the neighboring β-sheet layers together. We used MM-PBSA approach free energy calculations to determine the role of 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 observed trend in the entropy difference between wildtype and glycine mutant. 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 may provide valuable insight for amyloid aggregation inhibitor design efforts.

Can molecular dynamics simulations assist in design of specific inhibitors and imaging agents of amyloid aggregation? Structure, stability and free energy predictions for amyloid oligomers of VQIVYK, MVGGVV and LYQLEN

The aggregation modes of hexapeptide fragments of Tau, Insulin and Aβ peptide (VQIVYK, MVGGVV and LYQLEN) were found from their microcrystalline structures that had been recently resolved by X-ray analysis. The atomic structures reveal a dry self-complementary interface between the neighboring β-sheet layers, termed "steric zipper". In this study we perform several all-atom molecular dynamics simulations with explicit water to analyze stability of the crystalline fragments of 2-10 hexapeptides each and their analogs with single glycine replacement mutations to investigate the structural stability, aggregation behavior and thermodynamic of the amyloid oligomers. Upon comparing single and double layer models, our results reveal that additional strands contribute significantly to the structural stability of the peptide oligomers for double layer model, while in the case of single layer model the stability decreases (or remains the same in the case of LYQLEN). This is in agreement with the previous studies performed on different types of amyloid models. We also replaced the side-chains participating in the steric zipper interfaces with glycine. None of the mutants were structurally stable compared to the respective wild type model, except for mutants V2G and V6G in MVGGVV2 case. The exception can be explained by structural features of this particular polymorph. The double layer decamer and dodecamer aggregates of the wild type hexapeptides appear to be stable at 300K, which is confirmed by the conservation of high anti-parallel β-sheet content throughout the whole simulation time. Deletions of the side chains resulted in decline of secondary structure content compared to corresponding wild type indicating that the role of the replaced amino acid in stabilizing the structure. Detailed analysis of the binding energy reveals that stability of these peptide aggregates is determined mainly by the van der Waals and hydrophobic forces that can serve as quantitative measure of shape complementarities between the side chains. This observation implies that interactions among side chains forming the dehydrated steric zipper, rather than among those exposed to water, are the major structural determinant. The electrostatic repulsion destabilizes the studied double layer aggregates in two cases, while stabilizes the other two. Negative total binding free energy indicates that both wild type and mutants complex formation is favorable. However, the mutants complexation is less favorable than the wild type's. The present study provides the atomic level understanding of the aggregation behavior and the driving force for the amyloid aggregates, and could be useful for rational design of amyloid inhibitors and amyloid-specific biomarkers for diagnostic purposes.

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.

GNNQQNY--investigation of early steps during amyloid formation

Protein aggregation has been implicated in the pathology of several neurodegenerative diseases, and a better understanding of how it proceeds is essential for the development of therapeutic strategies. Recently, the amyloidogenic heptapeptide GNNQQNY has emerged as a molecule of choice for fundamental studies of protein aggregation. A number of experimental and computational studies have examined the structure of the GNNQQNY aggregate. Less work, however, has been aimed at understanding its aggregation pathway. In this study, we present a detailed computational analysis of such a pathway. To that end, transition path sampling Monte Carlo simulations are used to examine the dimerization process. A statistical analysis of the reaction pathways shows that the dimerization reaction proceeds via a zipping mechanism, initiated with the formation of distinct contacts at the third residue (N). Asparagine residues are found to play a key role in the early stages of aggregation. And, contrary to previous belief, it is also shown that the tyrosine terminal group is not required to stabilize the dimer. In fact, an asparagine residue leads to faster aggregation of the peptide.

Characterizing the assembly of the Sup35 yeast prion fragment, GNNQQNY: structural changes accompany a fiber-to-crystal switch

Amyloid-like fibrils can be formed by many different proteins and peptides. The structural characteristics of these fibers are very similar to those of amyloid fibrils that are deposited in a number of protein misfolding diseases, including Alzheimer's disease and the transmissible spongiform encephalopathies. The elucidation of two crystal structures from an amyloid-like fibril-forming fragment of the yeast prion, Sup35, with sequence GNNQQNY, has contributed to knowledge regarding side-chain packing of amyloid-forming peptides. Both structures share a cross-beta steric zipper arrangement but vary in the packing of the peptide, particularly in terms of the tyrosine residue. We investigated the fibrillar and crystalline structure and assembly of the GNNQQNY peptide using x-ray fiber diffraction, electron microscopy, intrinsic and quenched tyrosine fluorescence, and linear dichroism. Electron micrographs reveal that at concentrations between 0.5 and 10 mg/mL, fibers form initially, followed by crystals. Fluorescence studies suggest that the environment of the tyrosine residue changes as crystals form. This is corroborated by linear dichroism experiments that indicate a change in the orientation of the tyrosine residue over time, which suggests that a structural rearrangement occurs as the crystals form. Experimental x-ray diffraction patterns from fibers and crystals also suggest that these species are structurally distinct. A comparison of experimental and calculated diffraction patterns contributes to an understanding of the different arrangements accessed by the peptide.

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.