Oligomerization of amyloid Abeta16-22 peptides using hydrogen bonds and hydrophobicity forces

Uncovering supramolecular chirality codes for the design of tunable biomaterials

In neurodegenerative diseases, polymorphism and supramolecular assembly of β-sheet amyloids are implicated in many different etiologies and may adopt either a left- or right-handed supramolecular chirality. Yet, the underlying principles of how sequence regulates supramolecular chirality remains unknown. Here, we characterize the sequence specificity of the central core of amyloid-β 42 and design derivatives which enable chirality inversion at biologically relevant temperatures. We further find that C-terminal modifications can tune the energy barrier of a left-to-right chiral inversion. Leveraging this design principle, we demonstrate how temperature-triggered chiral inversion of peptides hosting therapeutic payloads modulates the dosed release of an anticancer drug. These results suggest a generalizable approach for fine-tuning supramolecular chirality that can be applied in developing treatments to regulate amyloid morphology in neurodegeneration as well as in other disease states.

Modulating Nanodroplet Formation En Route to Fibrillization of Amyloid Peptides with Designed Flanking Sequences

Soluble oligomers populating early amyloid aggregation can be regarded as nanodroplets of liquid-liquid phase separation (LLPS). Amyloid peptides typically contain hydrophobic aggregation-prone regions connected by hydrophilic linkers and flanking sequences, and such a sequence hydropathy pattern drives the formation of supramolecular structures in the nanodroplets and modulates subsequent fibrillization. Here, we studied LLPS and fibrillization of coarse-grained amyloid peptides with increasing flanking sequences. Nanodroplets assumed lamellar, cylindrical micellar, and spherical micellar structures with increasing peptide hydrophilic/hydrophobic ratios, and such morphologies governed subsequent fibrillization processes. Adding glycine-serine repeats as flanking sequences to Aβ16-22, the amyloidogenic core of amyloid-β, our computational predictions of morphological transitions were corroborated experimentally. The uncovered inter-relationships between the peptide sequence pattern, oligomer/nanodroplet morphology, and fibrillization pathway, kinetics, and structure may contribute to our understanding of pathogenic amyloidosis in aging, facilitate future efforts ameliorating amyloidosis through peptide engineering, and aid in the design of novel amyloid-based functional nanobiomaterials and nanocomposites.

Kinetic Mechanism of Amyloid-β-(16-22) Peptide Fibrillation

The kinetic mechanism of amyloid fibril formation by a peptide fragment containing seven residues of the amyloid-β protein Aβ-(16-22) was investigated. We found that the N- and C-terminal unprotected Aβ-(16-22), containing no aggregation nuclei, showed rapid fibrillation within seconds to minutes in a neutral aqueous buffer solution. The fibrillation kinetics were well described by the nucleation-elongation model, suggesting that primary nucleation was the rate-limiting step. On the basis of both experimental and theoretical analyses, the aggregated nucleus was estimated to be composed of 6-7 peptide molecules, wherein the two β-sheets were associated with their hydrophobic surfaces. Thin fibers with widths of 10-20 nm were formed, which increased their length and thickness, attaining a width of >20 nm over several tens of minutes, probably owing to the lateral association of the fibers. Electrostatic and hydrophobic interactions play important roles in aggregation. These results provide a basis for understanding the fibrillation of short peptides.

Development of a Hybrid-Resolution Force Field for Peptide Self-Assembly Simulations: Optimizing Peptide-Peptide and Peptide-Solvent Interactions

Atomic descriptions of peptide self-assembly are crucial to an understanding of disease-related peptide aggregation and the design of peptide-assembled materials. Obtaining these descriptions through computer simulation is challenging because current force fields, which were not designed for this process and are often unable to describe correctly peptide self-assembly behavior and the sequence dependence. Here, we developed a framework using dipeptide aggregation as a model system to improve force fields for simulations of self-assembly. Aggregation-related structural properties were designed and used to guide the optimization of peptide-peptide and peptide-solvent interactions. With this framework, we developed a self-assembly force field, termed PACE-ASM, by reoptimizing a hybrid-resolution force field that was originally developed for folding simulation. With its applicability in folding simulations, the new PACE was used to simulate the self-assembly of two disease-related short peptides, Aβ_16-21 and PHF6, into β-sheet-rich cross-β amyloids. These simulations reproduced the crystal structures of Aβ_16-21 and PHF6 amyloids at near-atomic resolution and captured the difference in packing orientations between the two sequences, a task which is challenging even with all-atom force fields. Apart from cross-β amyloids, the self-assembly of emerging helix-rich cross-α amyloids by another peptide PSMα3 can also be correctly described with the new PACE, manifesting the versatility of the force field. We demonstrated that the ability of the PACE-ASM to model peptide self-assembly is based largely on its improved description of peptide-peptide and peptide-solvent interactions. This was achieved with our optimization framework that can readily identify and address the deficiency in describing these interactions.

Simulations of Cross-Amyloid Aggregation of Amyloid-β and Islet Amyloid Polypeptide Fragments

Amyloid-beta (Aβ) and islet amyloid polypeptide (IAPP) are small peptides, classified as amyloids, that have the potential to self-assemble and form cytotoxic species, such as small soluble oligomers and large insoluble fibrils. The formation of Aβ aggregates facilitates the progression of Alzheimer's disease (AD), while IAPP aggregates induce pancreatic β-cell apoptosis, leading to exacerbation of Type 2 diabetes (T2D). Cross-amyloid interactions between Aβ and IAPP have been described both in vivo and in vitro, implying the role of Aβ or IAPP as modulators of cytotoxic self-aggregation of each species, and suggesting that Aβ-IAPP interactions are a potential molecular link between AD and T2D. Using molecular dynamics simulations, "hot spot" regions of the two peptides were studied to understand the formation of hexamers in a heterogenous and homogenous peptide-containing environment. Systems of only Aβ_(16-22) peptides formed antiparallel, β-barrel-like structures, while systems of only IAPP_(20-29) peptides formed stacked, parallel beta sheets and had relatively unstable aggregation structures after 2 μs of simulation time. Systems containing both Aβ and IAPP (1:1 ratio) hexamers showed antiparallel, β-barrel-like structures, with an interdigitated arrangement of Aβ_(16-22) and IAPP_(20-29). These β-barrel structures have features of cytotoxic amyloid species identified in previous literature. Ultimately, this work seeks to provide atomistic insight into both the mechanism behind cross-amyloid interactions and structural morphologies of these toxic amyloid species.

Computer Simulations Aimed at Exploring Protein Aggregation and Dissociation

Protein aggregation can lead to well-defined structures that are functional, but is also the cause of the death of neuron cells in many neurodegenerative diseases. The complexity of the molecular events involved in the aggregation kinetics of amyloid proteins and the transient and heterogeneous characters of all oligomers prevent high-resolution structural experiments. As a result, computer simulations have been used to determine the atomic structures of amyloid proteins at different association stages as well as to understand fibril dissociation. In this chapter, we first review the current computer simulation methods used for aggregation with some atomistic and coarse-grained results aimed at better characterizing the early formed oligomers and amyloid fibril formation. Then we present the applications of non-equilibrium molecular dynamics simulations to comprehend the dissociation of protein assemblies.

Design of carboxylated single-walled carbon nanotubes as highly efficient inhibitors against Aβ40 fibrillation based on the HyBER mechanism

Misfolding and the subsequent self-assembly of amyloid-β protein (Aβ) is very important in the occurrence of Alzheimer's disease (AD). Thus, inhibition of Aβ aggregation is currently an effective method to alleviate and treat AD. Herein, a carboxylated single-walled carbon nanotube (SWCNT-COOH) was rationally designed based on the hydrophobic binding-electrostatic repulsion (HyBER) mechanism. The inhibitory effect of SWCNT-COOH on Aβ fibrillogenesis was first studied. Based on the results of thioflavin T fluorescence and atomic force microscopy imaging assays, it was shown that SWCNT-COOH can not only effectively inhibit Aβ aggregation, but also depolymerize the mature fibrils of Aβ. In addition, its inhibitory action will be affected by the content of carboxyl groups. Moreover, the influence of SWCNT-COOH on cytotoxicity induced by Aβ was investigated by the MTT method. It was found that SWCNT-COOH can produce an anti-Aβ neuroprotective effect in vitro. Molecular dynamics simulations showed that SWCNT-COOH significantly destroyed the overall and internal structural stability of an Aβ40 trimer. Moreover, SWCNT-COOH interacted strongly with the N-terminal region, turn region and C-terminal region of the Aβ40 trimer via hydrogen bonds, salt bridges and π-π interactions, which triggered a large structural disturbance of the Aβ40 trimer, reduced the β-sheet content of the Aβ40 trimer and led to more disorder in these regions. All the above data not only reveal the suppressive effect of SWCNT-COOH on Aβ aggregation, but also reveal its inhibitory mechanism, which provides a useful clue to exploit anti-Aβ drugs in the future.

Modulating Aβ Fibrillogenesis with 'Trojan' peptides

Abnormal aggregation of beta-amyloid (Aβ) peptide into amyloid plaques in the brain has been identified as one of the key factors in instigating AD pathogenesis. Inhibition of Aβ aggregation can be an important therapeutic strategy in disease management. In this work, we demonstrate the application of structure-based design of short peptides ('trojan peptides'), intended to intervene in the aggregation of the core recognition domain of amyloid-beta peptide, a known malefactor in Alzheimer's disease. The modulatory effect of trojan peptides has been assessed using ThT fluorescence assay, FETEM imaging, IR, and toxicity assays on model neuronal cell lines. Experimental results suggest that designed trojan peptides could impede the aggregation of the core amyloid fibril forming segment of Aβ peptide, arrest the formation of toxic fibrillar assemblies, and reduce cytotoxicity of the neuronal cell lines.

Molecular insights into the inhibitory mechanisms of gallate moiety on the Aβ1-40 amyloid aggregation: A molecular dynamics simulation study

Alzheimer's disease is the most common form of neurodegenerative disease and the formation of Aβ amyloid aggregates has been widely demonstrated to be the principal cause of Alzheimer's disease. Our previous study and other studies suggested that the gallate moiety played an obligatory role in the inhibition process of naturally occurring polyphenols on Aβ amyloid fibrils formation. However, the detailed mechanisms were still unknown. Thus, in the present study, the gallic acid (GA) was specially selected and the molecular recognition mechanisms between GA molecules and Aβ_1-40 monomer were examined and analyzed by molecular dynamics simulation. The in silico experiments revealed that GA significantly prevented the conformational changes of Aβ_1-40 monomer with no β-sheet structure during the whole 100 ns. By analyzing the binding sites of GA molecules to Aβ_1-40 monomer, we found that both hydrophilic and hydrophobic amino acid residues were participated in the binding of GA molecules to Aβ_1-40 monomer. Moreover, results from the binding free energy analysis further demonstrated that the strength of polar interactions was significantly stronger than that of nonpolar interactions. We believed that our results could help to elucidate the underlying mechanisms of gallate moiety on the anti-amyloidogenic effects of polyphenols at the atomic level.

Conformational dynamics of amyloid-β (16–22) peptide in aqueous ionic liquids

Molecular dynamics simulations of amyloid-β (16–22) peptide dimer in water as well as at two different experimentally studied concentrations of hydrated ionic liquids (ILs), ethylammonium mesylate (EAM), ethylammonium nitrate (EAN), and triethylammonium mesylate (TEAM), were carried out employing an umbrella sampling method. We used the average Ψ angle of the peptide backbone as the reaction coordinate to observe the conformational changes of a peptide dimer. Secondary structural element values were calculated for the peptide dimer along the reaction coordinate to see the transition of the peptide dimer between β-sheet and α-helix conformations. We observe the β-sheet conformation as the global minimum on the free energy surfaces in both EAM and EAN ILs at both the concentrations and at a low concentration of TEAM. However, we observe α-helix conformation as the global minimum at a high concentration of TEAM. Our results are in good correlation with the experimental findings. We calculated the average number of intramolecular and intermolecular hydrogen bonds of α-helix and β-sheet conformations in all solutions, and they are in correlation with the secondary structure element values. To understand the peptide–IL interactions, atom–atom radial distribution functions of cation, anion, and water around amide oxygen and hydrogen atoms were calculated. The solvent-accessible surface area of the peptide dimer was calculated to understand the exposure of the peptide towards the solvent during conformational changes. Finally, van der Waals (vdW) and Coulomb interaction energies were calculated between peptide–cation, peptide–anion, and peptide–water to understand the stability of conformations in different concentrations. We find that the TEA cation has more vdW interaction energy compared to Coulomb interaction energy with peptide in 70% (w/w) TEAM, which mimics a membrane-like environment to induce α-helix conformation rather than β-sheet conformation.

Molecular dynamics simulations of amyloid-β (16–22) peptide dimer at two different experimentally studied concentrations of hydrated ethylammonium mesylate, ethylammonium nitrate, and triethylammonium mesylate were carried out employing an umbrella sampling method.

ATP Controls the Aggregation of Aβ16-22 Peptides

The oligomerization of Aβ_16-22 peptide, which is the hydrophobic core region of full-length Aβ_1-42, causes Alzheimer's disease (AD). This progressive neurodegenerative disease affects over 44 million people worldwide. However, very few synthesized drug molecules are available to inhibit the aggregation of Aβ. Recently, experimental studies have shown that the biological ATP molecule prevents Aβ fibrillation at the millimolar scale; however, the significance of ATP molecules on Aβ fibrillation and the mechanism behind it remain elusive. We have carried out a total of 7.5 μs extensive all-atom molecular dynamics and 8.82 μs of umbrella sampling in explicit water using AMBER14SB, AMBER99SB-ILDN, and AMBER-FB15 force fields for Aβ_16-22 peptide, to investigate the role of ATP on the disruption of Aβ_16-22 prefibrils. From various analyses, such as secondary structure analysis, residue-wise contact map, SASA, and interaction energies, we have observed that, in the presence of ATP, the aggregation of Aβ_16-22 peptide is very unfavorable. Moreover, the biological molecule ATP interacts with the Aβ_16-22 peptide via hydrogen bonding, π-π stacking, and NH-π interactions which, ultimately, prevent the aggregation of Aβ_16-22 peptide. Hence, we assume that the deficiency of ATP may cause Alzheimer's disease (AD).

Effects of All-Atom Molecular Mechanics Force Fields on Amyloid Peptide Assembly: The Case of Aβ16-22 Dimer

We investigated the effects of 17 widely used atomistic molecular mechanics force fields (MMFFs) on the structures and kinetics of amyloid peptide assembly. To this end, we performed large-scale all-atom molecular dynamics simulations in explicit water on the dimer of the seven-residue fragment of the Alzheimer's amyloid-β peptide, Aβ16-22, for a total time of 0.34 ms. We compared the effects of these MMFFs by analyzing various global reaction coordinates, secondary structure contents, the fibril population, the in-register and out-of-register architectures, and the fibril formation time at 310 K. While the AMBER94, AMBER99, and AMBER12SB force fields do not predict any β-sheets, the seven force fields, AMBER96, GROMOS45a3, GROMOS53a5, GROMOS53a6, GROMOS43a1, GROMOS43a2, and GROMOS54a7, form β-sheets rapidly. In contrast, the following five force fields, AMBER99-ILDN, AMBER14SB, CHARMM22*, CHARMM36, and CHARMM36m, are the best candidates for studying amyloid peptide assembly, as they provide good balances in terms of structures and kinetics. We also investigated the assembly mechanisms of dimeric Aβ16-22 and found that the fibril formation rate is predominantly controlled by the total β-strand content.

Thermodynamic phase diagram of amyloid-β (16-22) peptide

The aggregation of monomeric amyloid β protein (Aβ) peptide into oligomers and amyloid fibrils in the mammalian brain is associated with Alzheimer's disease. Insight into the thermodynamic stability of the Aβ peptide in different polymeric states is fundamental to defining and predicting the aggregation process. Experimental determination of Aβ thermodynamic behavior is challenging due to the transient nature of Aβ oligomers and the low peptide solubility. Furthermore, quantitative calculation of a thermodynamic phase diagram for a specific peptide requires extremely long computational times. Here, using a coarse-grained protein model, molecular dynamics (MD) simulations are performed to determine an equilibrium concentration and temperature phase diagram for the amyloidogenic peptide fragment Aβ_16-22 Our results reveal that the only thermodynamically stable phases are the solution phase and the macroscopic fibrillar phase, and that there also exists a hierarchy of metastable phases. The boundary line between the solution phase and fibril phase is found by calculating the temperature-dependent solubility of a macroscopic Aβ_16-22 fibril consisting of an infinite number of β-sheet layers. This in silico determination of an equilibrium (solubility) phase diagram for a real amyloid-forming peptide, Aβ_16-22, over the temperature range of 277-330 K agrees well with fibrillation experiments and transmission electron microscopy (TEM) measurements of the fibril morphologies formed. This in silico approach of predicting peptide solubility is also potentially useful for optimizing biopharmaceutical production and manufacturing nanofiber scaffolds for tissue engineering.

Mechanistic insight into E22Q-mutation-induced antiparallel-to-parallel β-sheet transition of Aβ16−22fibrils: an all-atom simulation study

E22Q mutation of Aβ₁₆₋₂₂fibrils facilitates parallel β-sheet formation by enhancing Q22–Q22 hydrogen-bonding interaction and A21–A21, F20–F20, F19–F19 and V18–V18 hydrophobic interaction.

On the Applicability of Force Fields To Study the Aggregation of Amyloidogenic Peptides Using Molecular Dynamics Simulations

Molecular dynamics simulations play an essential role in understanding biomolecular processes such as protein aggregation at temporal and spatial resolutions which are not attainable by experimental methods. For a correct modeling of protein aggregation, force fields must accurately represent molecular interactions. Here, we study the effect of five different force fields on the oligomer formation of Alzheimer's Aβ_16-22 peptide and two of its mutants: Aβ_16-22(F19V,F20V), which does not form fibrils, and Aβ_16-22(F19L) which forms fibrils faster than the wild type. We observe that while oligomer formation kinetics depends strongly on the force field, structural properties, such as the most relevant protein-protein contacts, are similar between them. The oligomer formation kinetics obtained with different force fields differ more from each other than the kinetics between aggregating and nonaggregating peptides simulated with a single force field. We discuss the difficulties in comparing atomistic simulations of amyloid oligomer formation with experimental observables.

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

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

Self-Assembly Pathways of β-Sheet-Rich Amyloid-β(1-40) Dimers: Markov State Model Analysis on Millisecond Hybrid-Resolution Simulations

Early oligomerization during amyloid-β (Aβ) aggregation is essential for Aβ neurotoxicity. Understanding how unstructured Aβs assemble into oligomers, especially those rich in β-sheets, is essential but remains challenging as the assembly process is too transient for experimental characterization and too slow for molecular dynamics simulations. So far, atomic simulations are limited only to studies of either oligomer structures or assembly pathways for short Aβ segments. To overcome the computational challenge, we combine in this study a hybrid-resolution model and adaptive sampling techniques to perform over 2.7 ms of simulations of formation of full-length Aβ40 dimers that are the earliest toxic oligomeric species. The Markov state model is further employed to characterize the transition pathways and associated kinetics. Our results show that for two major forms of β-sheet-rich structures reported experimentally, the corresponding assembly mechanisms are markedly different. Hairpin-containing structures are formed by direct binding of soluble Aβ in β-hairpin-like conformations. Formation of parallel, in-register structures resembling fibrils occurs ∼100-fold more slowly and involves a rapid encounter of Aβ in arbitrary conformations followed by a slow structural conversion. The structural conversion proceeds via diverse pathways but always requires transient unfolding of encounter complexes. We find that the transition kinetics could be affected differently by intra-/intermolecular interactions involving individual residues in a conformation-dependent manner. In particular, the interactions involving Aβ's N-terminal part promote the assembly into hairpin-containing structures but delay the formation of fibril-like structures, thus explaining puzzling observations reported previously regarding the roles of this region in the early assembly process.

Distinct oligomerization and fibrillization dynamics of amyloid core sequences of amyloid-beta and islet amyloid polypeptide

A direct observation of amyloid aggregation from isolated peptides to cross-β fibrils is crucial for understanding the nucleation-dependence process, but the corresponding macroscopic timescales impose a major computational challenge. Using rapid all-atom discrete molecular dynamics simulations, we capture the oligomerization and fibrillization dynamics of the amyloid core sequences of amyloid-β (Aβ) in Alzheimer's disease and islet amyloid polypeptide (IAPP) in type-2 diabetes, namely Aβ16-22 and IAPP22-28. Both peptides and their mixture spontaneously assemble into cross-β aggregates in silico, but follow distinct pathways. Aβ16-22 is highly aggregation-prone with a funneled free energy basin toward multi-layer β-sheet aggregates. IAPP22-28, on the other hand, features the accumulation of unstructured oligomers before the nucleation of β-sheets and growth into double-layer β-sheet aggregates. In the presence of Aβ16-22, the aggregation of IAPP22-28 is promoted by forming co-aggregated multi-layer β-sheets. Our study offers a detailed molecular insight to the long-postulated oligomerization-nucleation process in the amyloid aggregations.

The interaction with gold suppresses fiber-like conformations of the amyloid β (16-22) peptide

Inorganic surfaces and nanoparticles can accelerate or inhibit the fibrillation process of proteins and peptides, including the biomedically relevant amyloid β peptide. However, the microscopic mechanisms that determine such an effect are still poorly understood. By means of large-scale, state-of-the-art enhanced sampling molecular dynamics simulations, here we identify an interaction mechanism between the segments 16-22 of the amyloid β peptide, known to be fibrillogenic by itself, and the Au(111) surface in water that leads to the suppression of fiber-like conformations from the peptide conformational ensemble. Moreover, thanks to advanced simulation analysis techniques, we characterize the conformational selection vs. induced fit nature of the gold effect. Our results disclose an inhibition mechanism that is rooted in the details of the microscopic peptide-surface interaction rather than in general phenomena such as peptide sequestration from the solution.

Prediction of a stable associated liquid of short amyloidogenic peptides

Amyloid fibril formation is believed to be a nucleation-controlled process. Depending on the nature of peptide sequence, fibril nucleation can occur in one step, straight from a dilute solution, or in multiple steps via oligomers or disordered aggregates. What determines this process is poorly understood. Since the fibril formation kinetics is driven by thermodynamic forces, knowledge of the phase behavior is crucial. Here, we investigated the phase behavior of three short peptide sequences of varying side-chain hydrophobicity. Replica exchange molecular dynamics simulations of a mid-resolution model indicate that the weakly hydrophobic peptide forms fibrils directly from solution, whereas the most hydrophobic peptide forms a dense liquid phase before crystallizing into ordered fibrils at low temperatures. For the medium hydrophobic peptide we found evidence of a novel additional transition to a liquid phase consisting of clusters of aligned peptides, implying a three-step nucleation process. We tested the robustness of this prediction by applying Wertheim's theory and statistical associating fluid theory to a hard-sphere model dressed with isotropic and anisotropic attractions. We found that the ratio of interaction strengths strongly affects the phase behavior, and under certain conditions indeed gives rise to a stable polymerized liquid phase. The peptide clusters in the associated liquid tend to be slow and long-lived, which may give the oligomer droplet more time to act as a toxic oligomer, before turning into a fibril.

Amyloid Hydrogen Bonding Polymorphism Evaluated by (15)N{(17)O}REAPDOR Solid-State NMR and Ultra-High Resolution Fourier Transform Ion Cyclotron Resonance Mass Spectrometry

A combined approach, using Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) and solid-state NMR (Nuclear Magnetic Resonance), shows a high degree of polymorphism exhibited by Aβ species in forming hydrogen-bonded networks. Two Alzheimer's Aβ peptides, Ac-Aβ(16-22)-NH2 and Aβ(11-25), selectively labeled with (17)O and (15)N at specific amino acid residues were investigated. The total amount of peptides labeled with (17)O as measured by FTICR-MS enabled the interpretation of dephasing observed in (15)N{(17)O}REAPDOR solid-state NMR experiments. Specifically, about one-third of the Aβ peptides were found to be involved in the formation of a specific >C═(17)O···H-(15)N hydrogen bond with their neighbor peptide molecules, and we hypothesize that the rest of the molecules undergo ± n off-registry shifts in their hydrogen bonding networks.

Preformed template fluctuations promote fibril formation: insights from lattice and all-atom models

Fibril formation resulting from protein misfolding and aggregation is a hallmark of several neurodegenerative diseases such as Alzheimer's and Parkinson's diseases. Despite the fact that the fibril formation process is very slow and thus poses a significant challenge for theoretical and experimental studies, a number of alternative pictures of molecular mechanisms of amyloid fibril formation have been recently proposed. What seems to be common for the majority of the proposed models is that fibril elongation involves the formation of pre-nucleus seeds prior to the creation of a critical nucleus. Once the size of the pre-nucleus seed reaches the critical nucleus size, its thermal fluctuations are expected to be small and the resulting nucleus provides a template for sequential (one-by-one) accommodation of added monomers. The effect of template fluctuations on fibril formation rates has not been explored either experimentally or theoretically so far. In this paper, we make the first attempt at solving this problem by two sets of simulations. To mimic small template fluctuations, in one set, monomers of the preformed template are kept fixed, while in the other set they are allowed to fluctuate. The kinetics of addition of a new peptide onto the template is explored using all-atom simulations with explicit water and the GROMOS96 43a1 force field and simple lattice models. Our result demonstrates that preformed template fluctuations can modulate protein aggregation rates and pathways. The association of a nascent monomer with the template obeys the kinetics partitioning mechanism where the intermediate state occurs in a fraction of routes to the protofibril. It was shown that template immobility greatly increases the time of incorporating a new peptide into the preformed template compared to the fluctuating template case. This observation has also been confirmed by simulation using lattice models and may be invoked to understand the role of template fluctuations in slowing down fibril elongation in vivo.

Conformational features of the Aβ42 peptide monomer and its interaction with the surrounding solvent

Heterogeneous conformational flexibility of the Aβ monomers has been found to be correlated with the corresponding non-uniform entropy gains.

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.

Recent advances in the understanding of two-step nucleation of protein crystals

The two-step mechanism of nucleation of crystals in solutions posits that the formation of crystal nuclei occurs within structures of extended lifetimes, in which the nucleating solute is at high concentration. The validity of this mechanism has been demonstrated for proteins, small-molecule organic and inorganic materials, colloids, and polymers. Due to large molecule sizes, proteins are an ideal system to study the details of this nucleation pathway, in particular the formation mechanisms of the nucleation precursors and the associated physico-chemical rules. The precursors of protein crystal nuclei are protein-rich clusters of sizes ∼100 nm that contain 10 000–100 000 molecules and occupy less than 10⁻³of the total solution volume. Here we demonstrate, using oblique illumination microscopy, the liquid nature of the clusters of the protein lysozyme and reveal their inhomogeneous structure. We test a hypothesis put forth by theory that clusters primarily consist of transient protein oligomers. For this, we explore how varying the strength of the Coulomb interaction affects the cluster characteristics. We find that the cluster’s size is insensitive to variations of pH and ionic strength. In contrast, the addition of urea, a chaotropic agent that leads to protein unfolding, strongly decreases the cluster size. Shear stress, a known protein denaturant, induced by bubbling of the solutions with an inert gas, elicits a similar response. These observations support partial protein unfolding, followed by dimerization, as the mechanism of cluster formation. The amide hydrogen–deuterium exchange, monitored by nuclear magnetic resonance, highlights that lysozyme conformational flexibility is a condition for the formation of the protein-rich clusters and facilitates the nucleation of protein crystals.

Effects of hydroxylated carbon nanotubes on the aggregation of Aβ16-22 peptides: a combined simulation and experimental study

The pathogenesis of Alzheimer's disease (AD) is associated with the aggregation of amyloid-β (Aβ) peptides into toxic aggregates with ?-sheet character. In a previous computational study, we showed that pristine single-walled carbon nanotubes (SWCNTs) can inhibit the formation of β-sheet-rich oligomers in the central hydrophobic core fragment of Aβ (Aβ16-22). However, the poor solubility of SWCNTs in water hinders their use in biomedical applications and nanomedicine. Here, we investigate the influence of hydroxylated SWCNT, a water-soluble SWCNT derivative, on the aggregation of Aβ16-22 peptides using all-atom explicit-water replica exchange molecular dynamics simulations. Our results show that hydroxylated SWCNTs can significantly inhibit β-sheet formation and shift the conformations of Aβ16-22 oligomers from ordered β-sheet-rich structures toward disordered coil aggregates. Detailed analyses of the SWCNT-Aβ interaction reveal that the inhibition of β-sheet formation by hydroxylated SWCNTs mainly results from strong electrostatic interactions between the hydroxyl groups of SWCNTs and the positively charged residue K16 of Aβ16-22 and hydrophobic and aromatic stacking interactions between SWCNTs and F19 and F20. In addition, our atomic force microscopy and thioflavin T fluorescence experiments confirm the inhibitory effect of both pristine and hydroxylated SWCNTs on Aβ16-22 fibrillization, in support of our previous and present replica exchange molecular dynamics simulation results. These results demonstrate that hydroxylated SWCNTs efficiently inhibit the aggregation of Aβ16-22; in addition, they offer molecular insight into the inhibition mechanism, thus providing new clues for the design of therapeutic drugs against amyloidosis.

Investigation of the aggregation process of amyloid-β-(16-22) peptides and the dissolution of intermediate aggregates

The aggregation processes of amyloid-β-(16-22) peptides (Aβ16-22) are investigated by atomic force microscopy (AFM). It is found that Aβ16-22 peptides quickly aggregate from monomers to oligomers and flakelike structures and finally to fibrils. In particular, unusual morphology change is observed in an early stage of aggregation; that is, the originally formed flakelike structures would disappear in the following aggregation processes. To determine the evolution of the flakelike structures, in situ AFM imaging is carried out in liquid to reveal the real-time morphology change of Aβ16-22. The results provide clear evidence that the flakelike structures are in an unstable intermediate state, which would be dissolved into oligomers or short protofibrils for reorganization. Further fluorescence and attenuated total reflectance Fourier transform infrared (ATR-FTIR) experiments on thioflavin T(ThT) suggest that those flakelike structures contain β-sheet components.

Aβ(16-22) peptides can assemble into ordered β-barrels and bilayer β-sheets, while substitution of phenylalanine 19 by tryptophan increases the population of disordered aggregates

A recent experimental study reported that termini-uncapped Aβ(16-22) (with sequence KLVFFAE) peptides self-assembled into nanofibrils at pH 2.0. The oligomerization of this uncapped peptide at atomic level in acidic pH condition remains to be determined, as computational studies mainly focus on the self-assembly of capped Aβ(16-22) peptides at neutral pH condition. In this study, using replica exchange molecular dynamics (REMD) simulations with explicit solvent, we investigated the octameric structures of the uncapped Aβ(16-22) and its F19W variant at acidic pH condition. Our simulations reveal that the Aβ(16-22) octamers adopt various conformations, including closed β-barrels, bilayer β-sheets, and disordered aggregates. The closed β-barrel conformation is particularly interesting, as the cylindrical β-barrel has been reported recently as a cytotoxic species. Interpeptide contact probability analyses between all pairs of residues reveal that the hydrophobic and aromatic stacking interactions between F19 residues play an essential role in the formation of β-barrels and bilayer β-sheets. The importance of F19 and the steric effect on the structures of Aβ(16-22) octamers are further examined by REMD simulation of F19W mutant. This REMD run shows that substitution of F19 by W with a more bulky aromatic side chain significantly reduces the β-sheet content and in turn enhances the population of disordered aggregates, indicating that the steric effect significantly affect the self-assembly of low molecular weight Aβ(16-22) oligomers.

Oligomerization of Peptides LVEALYL and RGFFYT and Their Binding Affinity to Insulin

Recently it has been proposed a model for fibrils of human insulin in which the fibril growth proceeds via stacking LVEALYL (fragment 11-17 from chain B of insulin) into pairs of tightly interdigitated [Formula: see text]-sheets. The experiments have also shown that LVEALYL has high propensity to self-assembly and binding to insulin. This necessitates study of oligomerization of LVEALYL and its binding affinity to full-length insulin. Using the all-atom simulations with Gromos96 43a1 force field and explicit water it is shown that LVEALYL can aggregate. Theoretical estimation of the binding free energy of LVEALYL to insulin by the molecular mechanic Poisson-Boltzmann surface area method reveals its strong binding affinity to chain B, implying that, in agreement with the experiments, LVEALYL can affect insulin aggregation via binding mechanism. We predict that, similar to LVEALYL, peptide RGFFYT (fragment B22-27) can self-assemble and bind to insulin modulating its fibril growth process. The binding affinity of RGFFYT is shown to be comparable with that of LVEALYL.

Nucleation process of a fibril precursor in the C-terminal segment of amyloid-β

By extended atomistic simulations in explicit solvent and bias-exchange metadynamics, we study the aggregation process of 18 chains of the C-terminal segment of amyloid-β, an intrinsically disordered protein involved in Alzheimer's disease and prone to form fibrils. Starting from a disordered aggregate, we are able to observe the formation of an ordered nucleus rich in beta sheets. The rate limiting step in the nucleation pathway involves crossing a barrier of approximately 40 kcal/mol and is associated with the formation of a very specific interdigitation of the side chains belonging to different sheets. This structural pattern is different from the one observed experimentally in a microcrystal of the same system, indicating that the structure of a "nascent" fibril may differ from the one of an "extended" fibril.

The effects of ferulic acid on β-amyloid fibrillar structures investigated through experimental and computational techniques

BACKGROUND

Current research has indicated that small natural compounds could interfere with β-amyloid fibril growth and have the ability to disassemble preformed folded structures. Ferulic acid (FA), which possesses both hydrophilic and hydrophobic moieties and binds to peptides/proteins, is a potential candidate against amyloidogenesis. The molecular mechanisms connected to this action have not been elucidated in detail yet.

METHODS

Here the effects of FA on preformed fibrils are investigated by means of a concerted experimental-computational approach. Spectroscopic techniques, such as FTIR, fluorescence, size exclusion chromatography and confocal microscopy in combination with molecular dynamics simulations are used to identify those features which play a key role in the destabilization of the aggregates.

RESULTS

Experimental findings highlight that FA has disruptive effects on the fibrils. The computational analysis suggests that dissociation of peptides from the amyloid superstructures could take place along the fibril axis and be primarily determined by the cooperative rupture of the backbone hydrogen bonds and of the Asp-Lys salt bridges.

CONCLUSION

FA clusters could induce a sort of stabilization and tightening of the fibril structure in the short term and its disruption in the long term, inhibiting further fibril re-assembly through FA screening effects.

GENERAL SIGNIFICANCE

The combination of experimental and computational techniques could be successfully used to identify the disrupting action of FA on preformed Aβ fibrils in water solution.

Nanomaterials design and tests for neural tissue engineering

Nanostructured scaffolds recently showed great promise in tissue engineering: nanomaterials can be tailored at the molecular level and scaffold morphology may more closely resemble features of extracellular matrix components in terms of porosity, framing and biofunctionalities. As a consequence, both biomechanical properties of scaffold microenvironments and biomaterial-protein interactions can be tuned, allowing for improved transplanted cell engraftment and better controlled diffusion of drugs. Easier said than done, a nanotech-based regenerative approach encompasses different fields of know-how, ranging from in silico simulations, nanomaterial synthesis and characterization at the nano-, micro- and mesoscales to random library screening methods (e.g. phage display), in vitro cellular-based experiments and validation in animal models of the target injury. All of these steps of the "assembly line" of nanostructured scaffolds are tightly interconnected both in their standard analysis techniques and in their most recent breakthroughs: indeed their efforts have to jointly provide the deepest possible analyses of the diverse facets of the challenging field of neural tissue engineering. The purpose of this review is therefore to provide a critical overview of the recent advances in and drawbacks and potential of each mentioned field, contributing to the realization of effective nanotech-based therapies for the regeneration of peripheral nerve transections, spinal cord injuries and brain traumatic injuries. Far from being the ultimate overview of such a number of topics, the reader will acknowledge the intrinsic complexity of the goal of nanotech tissue engineering for a conscious approach to the development of a regenerative therapy and, by deciphering the thread connecting all steps of the research, will gain the necessary view of its tremendous potential if each piece of stone is correctly placed to work synergically in this impressive mosaic.

Exploring the inter-molecular interactions in amyloid-β protofibril with molecular dynamics simulations and molecular mechanics Poisson-Boltzmann surface area free energy calculations

Aggregation of amyloid-β (Aβ) peptides correlates with the pathology of Alzheimer's disease. However, the inter-molecular interactions between Aβ protofibril remain elusive. Herein, molecular mechanics Poisson-Boltzmann surface area analysis based on all-atom molecular dynamics simulations was performed to study the inter-molecular interactions in Aβ(17-42) protofibril. It is found that the nonpolar interactions are the important forces to stabilize the Aβ(17-42) protofibril, while electrostatic interactions play a minor role. Through free energy decomposition, 18 residues of the Aβ(17-42) are identified to provide interaction energy lower than -2.5 kcal/mol. The nonpolar interactions are mainly provided by the main chain of the peptide and the side chains of nine hydrophobic residues (Leu17, Phe19, Phe20, Leu32, Leu34, Met35, Val36, Val40, and Ile41). However, the electrostatic interactions are mainly supplied by the main chains of six hydrophobic residues (Phe19, Phe20, Val24, Met35, Val36, and Val40) and the side chains of the charged residues (Glu22, Asp23, and Lys28). In the electrostatic interactions, the overwhelming majority of hydrogen bonds involve the main chains of Aβ as well as the guanidinium group of the charged side chain of Lys28. The work has thus elucidated the molecular mechanism of the inter-molecular interactions between Aβ monomers in Aβ(17-42) protofibril, and the findings are considered critical for exploring effective agents for the inhibition of Aβ aggregation.

Dynamics of metastable β-hairpin structures in the folding nucleus of amyloid β-protein

The amyloid β-protein (Aβ), which is present predominately as a 40- or 42-residue peptide, is postulated to play a seminal role in the pathogenesis of Alzheimer's disease (AD). Folding of the Aβ(21-30) decapeptide region is a critical step in the aggregation of Aβ. We report results of constant temperature all-atom molecular dynamics simulations in explicit water of the dynamics of monomeric Aβ(21-30) and its Dutch [Glu22Gln], Arctic [Glu22Gly], and Iowa [Asp23Asn] isoforms that are associated with familial forms of cerebral amyloid angiopathy and AD. The simulations revealed a variety of loop conformers that exhibited a hydrogen bond network involving the Asp23 and Ser26 amino acids. A population of conformers, not part of the loop population, was found to form metastable β-hairpin structures with the highest probability in the Iowa mutant. At least three β-hairpin structures were found that differed in their hydrogen bonding register, average number of backbone hydrogen bonds, and lifetimes. Analysis revealed that the Dutch mutant had the longest β-hairpin lifetime (≥500 ns), closely followed by the Iowa mutant (≈500 ns). Aβ(21-30) and the Arctic mutant had significantly lower lifetimes (≈200 ns). Hydrophobic packing of side chains was responsible for enhanced β-hairpin lifetimes in the Dutch and Iowa mutants, whereas lifetimes in Aβ(21-30) and its Arctic mutant were influenced by the backbone hydrogen bonding. The data suggest that prolonged β-hairpin lifetimes may impact peptide pathogenicity in vivo.