Peptide Sequence and Amyloid Formation

De novo designed aliphatic and aromatic peptides assemble into amyloid-like cytotoxic supramolecular nanofibrils

Peptides are very interesting biomolecules that upon self-association form a variety of thermodynamically stable supramolecular structures of nanometric dimension e.g. nanotubes, nanorods, nanovesicles, nanofibrils, nanowires and many others. Herein, we report six peptide molecules having a general chemical structure, H-Gaba-X-X-OH (Gaba: γ-aminobutyric acid, X: amino acid). Out of these six peptides, three are aromatic and the others are aliphatic. Atomic force microscopic (AFM) studies reveal that except peptide 6 (H-Gaba-Trp-Trp-OH), all the reported peptides adopt nanofibrillar morphology upon aggregation in aqueous medium. These supramolecular assemblies can recognize amyloid-specific molecular probe congo red (CR) and thioflavine t (ThT) and exhibit all the characteristic properties of amyloids. The MTT cell viability assay reveals that the toxicity of both aliphatic and aromatic peptides increases with increasing concentration of the peptides to both cancer (HeLa) and non-cancer (HEK 293) cells. Of note, the aromatic peptides show a slightly higher cytotoxic effect compared to the aliphatic peptides. Overall, the studies highlight the self-assembling nature of the de novo designed aliphatic and aromatic peptides and pave the way towards elucidating the intricacies of pathogenic amyloid assemblies.

Isotope-edited vibrational circular dichroism study reveals a flexible N-terminal structure of islet amyloid peptide (NFGAIL) in amyloid fibril form: A site-specific local structural analysis

The short peptide fragment NFGAIL (IAPf) is a well-known amyloidogenic peptide (22-27), derived from human islet amyloid polypeptide(hIAPP), whose fibrillar structure is often used to better understand the wild-type hIAPP amyloid fibrils, associated with type II diabetes. Despite an extensive study, the fibrillar structure of IAPf at the amino acid residue level is still unclear. Herein, the vibrational circular dichroism(VCD) spectroscopic technique coupled with isotope labelling strategy has been used to study the site-specific local structure of IAPf amyloid fibrils. Two ¹³C labeled IAPfs were designed and used along with unlabelled IAPf to achieve this. The ¹³C labelled (on -C=O) glycine(IAPf-G) and phenylalanine (IAPf-F) residues were introduced into the IAPf sequence separately by replacing natural glycine (residue 24) and phenylalanine (residue 23), respectively. VCD spectral analysis on IAPf-G suggests that IAPf fibrils adopt parallel β-sheet conformation with glycine residues are part of β-sheet and in-register. Unlike IAPf-G, VCD analysis on IAPf-F reveals that phenylalanine residues exist in the turn/hairpin conformation rather than β-sheet region. Both VCD results thus suggest that IAPf amyloid fibril consists of a mixture of β-sheet as a major conformation involving GAIL and turn/hairpin as a minor conformation involving NF rather than an idealized β-sheet involving all the amino acids. While previous studies speculated that the full NFGAIL sequence could participate in the β-sheet formation, the present site-specific structural analysis of IAPf amyloid fibrils at residue level using isotope-edited VCD has gained significant attention. Such residue level information has important implications for understanding the role of NFGAIL sequence in the amyloid fibrillation of hIAPP.

Tyrosine carbon dots inhibit fibrillation and toxicity of the human islet amyloid polypeptide

Misfolding and aggregation of the human islet amyloid polypeptide (hIAPP) are believed to play key roles in the pathophysiology of type-II diabetes. Here, we demonstrate that carbon dots (C-dots) prepared from the amino acid tyrosine inhibit fibrillation of hIAPP, reduce hIAPP-induced cell toxicity and block membrane disruption by the peptide. The pronounced inhibitory effect is traced to the display of ubiquitous aromatic residues upon the C-dots' surface, mimicking the anti-fibril and anti-toxic activity of natural polyphenolic compounds. Notably, spectroscopy and thermodynamics analysis demonstrated different hIAPP interactions and fibril inhibition effects induced by tyrosine-C-dots displaying phenolic residues and C-dots prepared from phenylalanine which exhibited phenyl units on their surface, underscoring the significance of hydrogen bonding mediated by the phenolic hydroxide moieties for the fibril modulation activity. The presented experiments attest to the potential of tyrosine-C-dots as a therapeutic vehicle for protein misfolding diseases, interfering in both π-π interactions as well as hydrogen bonding involving aromatic residues of amyloidogenic peptides.

Role of aromatic amino acids in amyloid self-assembly

Amyloids are proteins of a cross-β structure found as deposits in several diseases and also in normal tissues (nails, spider net, silk). Aromatic amino acids are frequently found in amyloid deposits. Although they are not indispensable, aromatic amino acids, phenylalanine, tyrosine and tryptophan, enhance significantly the kinetics of formation and thermodynamic stability, while tape or ribbon-like morphology is represented in systems with experimentally detected π-π interactions between aromatic rings. Analysis of geometries and energies of the amyloid PDB structures indicate the prevalence of aromatic-nonaromatic interactions and confirm that aromatic-aromatic interactions are not crucial for the amyloid formation.

Determination of the Aggregate Binding Site of Amyloid Protofibrils Using Electron Capture Dissociation Tandem Mass Spectrometry

Amyloid fibril formation is a hallmark in a range of human diseases. Analysis of the molecular details of amyloid aggregation, however, is limited by the difficulties in solubilizing, separating, and identifying the aggregated biomolecules. Additional labeling or protein modification is required in many current analytical techniques in order to provide molecular details of amyloid protein aggregation, but these modifications may result in protein structure disruption. Herein, ultrahigh resolution mass spectrometry (MS) with electron capture dissociation tandem MS (ECD MS/MS) has been applied to monitor the formation of early oligomers of human islet amyloid polypeptide (hIAPP), which aggregate rapidly in the pancreas of type II diabetes (T2D) patients. ECD MS/MS results show the aggregation region of the early oligomers is at the Ser-28/Ser-29 residue of a hIAPP unit and at the Asn-35 residue of another hIAPP unit near the C-terminus in the gas phase. These data contribute to the understanding of the binding site between hIAPP units which may help for specific target region therapeutic development in the future. Furthermore, MS has also been applied to quantify the amount of soluble amyloid protein remaining in the incubated solutions, which can be used to estimate the aggregation rate of amyloid protein during incubation (28 days). These data are further correlated with the results obtained using fluorescence spectroscopy and transmission electron microscopy (TEM) to generate a general overview of amyloid protein aggregation. The methods demonstrated in this article not only explore the aggregation site of hIAPP down to an amino acid residue level, but are also applicable to many amyloid protein aggregation studies.

Large-scale all-atom molecular dynamics alanine-scanning of IAPP octapeptides provides insights into the molecular determinants of amyloidogenicity

In order to investigate the early phase of the amyloid formation by the short amyloidogenic octapeptide sequence (‘NFGAILSS’) derived from IAPP, we carried out a 100ns all-atom molecular dynamics (MD) simulations of systems that contain 27 peptides and over 30,000 water molecules. The large-scale calculations were performed for the wild type sequence and seven alanine-scanned sequences using AMBER 8.0 on RIKEN’s special purpose MD-GRAPE3 supercomputer, using the all-atom point charge force field ff99, which do not favor β-structures. Large peptide clusters (size 18–26 mers) were observed for all simulations, and our calculations indicated that isoleucine at position 5 played important role in the formation of β-rich clusters. In the oligomeric state, the wild type and the S7A sequences had the highest β-structure content (~14%), as calculated by DSSP, in line with experimental observations, whereas I5A and G3A had the highest helical content (~20%). Importantly, the β-structure preferences of wild type IAPP originate from its association into clusters and are not intrinsic to its sequence. Altogether, the results of this first large-scale, multi-peptide all-atom molecular dynamics simulation appear to provide insights into the mechanism of amyloidogenic and non-amyloidogenic oligomers that mainly corroborate previous experimental observations.

Designing phenylalanine-based hybrid biological materials: controlling morphology via molecular composition

Harnessing the self-assembly of peptide sequences has demonstrated great promise in the domain of creating high precision shape-tunable biomaterials. The unique properties of peptides allow for a building block approach to material design. In this study, self-assembly of mixed systems encompassing two peptide sequences with identical hydrophobic regions and distinct polar segments is investigated. The two peptide sequences are diphenylalanine and phenylalanine-asparagine-phenylalanine. The study examines the impact of molecular composition (namely, the total peptide concentration and the relative tripeptide concentration) on the morphology of the self-assembled hybrid biological material. We report a rich polymorphism in the assemblies of these peptides and explain the relationship between the peptide sequence, concentration and the morphology of the supramolecular assembly.

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

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

Distinct C9orf72-Associated Dipeptide Repeat Structures Correlate with Neuronal Toxicity

Hexanucleotide repeat expansions in C9orf72 are the most common inherited cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). The expansions elicit toxicity in part through repeat-associated non-AUG (RAN) translation of the intronic (GGGGCC)n sequence into dipeptide repeat-containing proteins (DPRs). Little is known, however, about the structural characteristics and aggregation propensities of the dipeptide units comprising DPRs. To address this question, we synthesized dipeptide units corresponding to the three sense-strand RAN translation products, analyzed their structures by circular dichroism, electron microscopy and dye binding assays, and assessed their relative toxicity when applied to primary cortical neurons. Short, glycine-arginine (GR)3 dipeptides formed spherical aggregates and selectively reduced neuronal survival compared to glycine-alanine (GA)3 and glycine-proline (GP)3 dipeptides. Doubling peptide length had little effect on the structure of GR or GP peptides, but (GA)6 peptides formed β-sheet rich aggregates that bound thioflavin T and Congo red yet lacked the typical fibrillar morphology of amyloids. Aging of (GA)6 dipeptides increased their β-sheet content and enhanced their toxicity when applied to neurons. We also observed that the relative toxicity of each tested dipeptide was proportional to peptide internalization. Our results demonstrate that different C9orf72-related dipeptides exhibit distinct structural properties that correlate with their relative toxicity.

Different nanostructures caused by competition of intra- and inter-β-sheet interactions in hierarchical self-assembly of short peptides

To understand how molecular interactions lead to the self-assembly of twisted, helical and flat nanoribbons, we have compared the hierarchical self-assembly processes of three selected octapeptides with the same amino acid composition but different sequences by both experiments and molecular dynamics (MD) simulations. KE-F8 (NH2-KEFFFFKE-CONH2) and EK-F8 (NH2-KEFFFFEK-CONH2) have the same distribution of hydrophobic residues and only differ by swapping the positive and negative charged residues at their C-terminals, while KFE-8 (NH2-KFEFKFEF-CONH2) differs from KE-F8 and EK-F8 by having all hydrophobic and charged residues evenly distributed. MD simulations indicated that the competition between electrostatic and hydrophobic interactions at the molecular level results in different initial packing modes: KE-F8 monomers form completely matched anti-parallel β-sheets, EK-F8 monomers align with one residue shifting, and KFE-8 monomers pack β-sheets with two heterogeneous surfaces, consistent with previously suggested models. Driven by inter-strand and inter-sheet interactions, further growth of these molecular templates leads to larger oligomers with different twisting and stacking degrees, which are structurally consistent with the experimentally observed self-assembled morphologies. Further MD simulations showed that the competition between intra-β-sheet and inter-β-sheet interactions is responsible for the different twisting and stacking degrees of β-sheets and the subsequent formation of different nanostructures (twisted ribbons for KE-F8, helical ribbons/tubes for EK-F8 and flat ribbons for KFE-8). This study thus provided an important mechanistic insight into the fine tuning of molecular packing and interactions via peptide sequence variation leading to controllable self-assembly of twisted, helical and flat nanostructures.

Two-dimensional sum-frequency generation (2D SFG) spectroscopy: summary of principles and its application to amyloid fiber monolayers

By adding a mid-infrared pulse shaper to a sum-frequency generation (SFG) spectrometer, we have built a 2D SFG spectrometer capable of measuring spectra analogous to 2D IR spectra but with monolayer sensitivity and SFG selection rules. In this paper, we describe the experimental apparatus and provide an introduction to 2D SFG spectroscopy to help the reader interpret 2D SFG spectra. The main aim of this manuscript is to report 2D SFG spectra of the amyloid forming peptide FGAIL. FGAIL is a critical segment of the human islet amyloid polypeptide (hIAPP or amylin) that aggregates in people with type 2 diabetes. FGAIL is catalyzed into amyloid fibers by many types of surfaces. Here, we study the structure of FGAIL upon deposition onto a gold surface covered with a self-assembled monolayer of methyl-4-mercaptobenzoate (MMB) that produces an ester coating. FGAIL deposited on bare gold does not form ordered layers. The measured 2D SFG spectrum is consistent with amyloid fiber formation, exhibiting both the parallel (a+) and perpendicular (a-) symmetry modes associated with amyloid β-sheets. Cross peaks are observed between the ester stretches of the coating and the FGAIL peptides. Simulations are presented for two possible structures of FGAIL amyloid β-sheets that illustrate the sensitivity of the 2D SFG spectra to structure and orientation. These results provide some of the first molecular insights into surface catalyzed amyloid fiber structure.

Energy interactions in amyloid-like fibrils from NNQQNY

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

Full length amylin oligomer aggregation: insights from molecular dynamics simulations and implications for design of aggregation inhibitors

Amyloid oligomers are considered to play essential roles in the pathogenesis of amyloid-related degenerative diseases including type 2 diabetes. Using an explicit solvent all atomic MD simulation, we explored the stability, conformational dynamics and association force of different single-layer models of the full-length wild-type and glycine mutants of amylin (pentamer) obtained from a recent high resolution fibril model. The RMSF profile shows enhanced flexibility in the disorder (Lys1-Cys7) and turn region (Ser19-Gly23), along with smallest fluctuation at the residues (Asn14-Phe15-Leu16-Val17-His18) of β1 region and (Ala25-Ile26-Leu27-Ser28-Ser29) of the β2 region. We obtained a significant difference in backbone RMSD between the wild-type and the mutants, indicating that mutations affected the stability of the peptide. The RMSD and RMSF profiles indicate the edge and loop residues are the primary contributors to the overall conformational changes. The degree of structural similarity between the oligomers in the simulation and the fibril conformation is proposed as the possible explanation for experimentally observed shortening of the nucleation lag phase of amylin with oligomer seeding. On the basis of structure-stability findings, the β1 and β2 portions are optimal target for further anti-amyloid drug design. The MM-PBSA binding energy calculation reveals the binding of amylin: amylin strands in single layer is dominated by contributions from van der Waals interactions. The non-polar solvation term is also found to be favorable. While the electrostatic interactions and polar solvation energy was found to be favorable for the interaction for the larger aggregate and unfavorable for the smaller aggregates. A per-residue decomposition of the binding free energy has been performed to identify the residues contributing most to the self-association free energy. Residues found in the β-sheet regions were found to be key residue making the largest favorable contributions to the single-layer association. The result from our simulation could be used in rational design of new amylinomimetic agent, amylin aggregation inhibitors and amylin-specific biomarkers.

Observation of molecular inhibition and binding structures of amyloid peptides

Unveiling interactions between labeling molecules and amyloid fibrils is essential to develop new detection methods for studying amyloid structures under various conditions. This review endeavours to reflect the progress in studying interactions between molecular inhibitors and amyloid peptides using a series of experimental approaches, such as X-ray diffraction, nuclear magnetic resonance, scanning probe microscopy, and electron microscopy. The revealed binding mechanisms of anti-amyloid drugs and target proteins could benefit the rational design of drugs for prevention or treatment of amyloidal diseases.

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

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

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.

Complete phenotypic recovery of an Alzheimer's disease model by a quinone-tryptophan hybrid aggregation inhibitor

The rational design of amyloid oligomer inhibitors is yet an unmet drug development need. Previous studies have identified the role of tryptophan in amyloid recognition, association and inhibition. Furthermore, tryptophan was ranked as the residue with highest amyloidogenic propensity. Other studies have demonstrated that quinones, specifically anthraquinones, can serve as aggregation inhibitors probably due to the dipole interaction of the quinonic ring with aromatic recognition sites within the amyloidogenic proteins. Here, using in vitro, in vivo and in silico tools we describe the synthesis and functional characterization of a rationally designed inhibitor of the Alzheimer's disease-associated β-amyloid. This compound, 1,4-naphthoquinon-2-yl-L-tryptophan (NQTrp), combines the recognition capacities of both quinone and tryptophan moieties and completely inhibited Aβ oligomerization and fibrillization, as well as the cytotoxic effect of Aβ oligomers towards cultured neuronal cell line. Furthermore, when fed to transgenic Alzheimer's disease Drosophila model it prolonged their life span and completely abolished their defective locomotion. Analysis of the brains of these flies showed a significant reduction in oligomeric species of Aβ while immuno-staining of the 3^rd instar larval brains showed a significant reduction in Aβ accumulation. Computational studies, as well as NMR and CD spectroscopy provide mechanistic insight into the activity of the compound which is most likely mediated by clamping of the aromatic recognition interface in the central segment of Aβ. Our results demonstrate that interfering with the aromatic core of amyloidogenic peptides is a promising approach for inhibiting various pathogenic species associated with amyloidogenic diseases. The compound NQTrp can serve as a lead for developing a new class of disease modifying drugs for Alzheimer's disease.

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

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

Self-assembly of phenylalanine oligopeptides: insights from experiments and simulations

Studies of peptide-based nanostructures provide general insights into biomolecular self-assembly and can lead material engineering toward technological applications. The diphenylalanine peptide (FF) self-assembles into discrete, hollow, well ordered nanotubes, and its derivatives form nanoassemblies of various morphologies. Here we demonstrate for the first time, to our knowledge, the formation of planar nanostructures with beta-sheet content by the triphenylalanine peptide (FFF). We characterize these structures using various microscopy and spectroscopy techniques. We also obtain insights into the interactions and structural properties of the FF and FFF nanostructures by 0.4-micros, implicit-solvent, replica-exchange, molecular-dynamics simulations of aqueous FF and FFF solutions. In the simulations the peptides form aggregates, which often contain open or ring-like peptide networks, as well as elementary and network-containing structures with beta-sheet characteristics. The networks are stabilized by polar and nonpolar interactions, and by the surrounding aggregate. In particular, the charged termini of neighbor peptides are involved in hydrogen-bonding interactions and their aromatic side chains form "T-shaped" contacts, as in three-dimensional FF crystals. These interactions may assist the FF and FFF self-assembly at the early stage, and may also stabilize the mature nanostructures. The FFF peptides have higher network propensities and increased aggregate stabilities with respect to FF, which can be interpreted energetically.

Prediction of structural stability of short beta-hairpin peptides by molecular dynamics and knowledge-based potentials

BACKGROUND

The structural stability of peptides in solution strongly affects their binding affinities and specificities. Thus, in peptide biotechnology, an increase in the structural stability is often desirable. The present work combines two orthogonal computational techniques, Molecular Dynamics and a knowledge-based potential, for the prediction of structural stability of short peptides (< 20 residues) in solution.

RESULTS

We tested the new approach on four families of short beta-hairpin peptides: TrpZip, MBH, bhpW and EPO, whose structural stabilities have been experimentally measured in previous studies. For all four families, both computational techniques show considerable correlation (r > 0.65) with the experimentally measured stabilities. The consensus of the two techniques shows higher correlation (r > 0.82).

CONCLUSION

Our results suggest a prediction scheme that can be used to estimate the relative structural stability within a peptide family. We discuss the applicability of this predictive approach for in-silico screening of combinatorial peptide libraries.

Investigating the mechanism of peptide aggregation: insights from mixed monte carlo-molecular dynamics simulations

The early stages of peptide aggregation are currently not accessible by experimental techniques at atomic resolution. In this article, we address this problem through the application of a mixed simulation scheme in which a preliminary coarse-grained Monte Carlo analysis of the free-energy landscape is used to identify representative conformations of the aggregates and subsequent all-atom molecular dynamics simulations are used to analyze in detail possible pathways for the stabilization of oligomers. This protocol was applied to systems consisting of multiple copies of the model peptide GNNQQNY, whose detailed structures in the aggregated state have been recently solved in another study. The analysis of the various trajectories provides dynamical and structural insight into the details of aggregation. In particular, the simulations suggest a hierarchical mechanism characterized by the initial formation of stable parallel beta-sheet dimers and identify the formation of the polar zipper motif as a fundamental feature for the stabilization of initial oligomers. Simulation results are consistent with experimentally derived observations and provide an atomically detailed view of the putative initial stages of fibril formation.

Computational studies of the structure, dynamics and native content of amyloid‐like fibrils of ribonuclease A

The characterization at atomic resolution of amyloid-like protein aggregates is one of the fundamental problems of modern biology. In particular, the question whether native-like domains are retained or completely refolded in the amyloid state and the identification of possible mechanisms for macromolecular ordered aggregation represent major unresolved puzzles. To address these issues, in this article we examine the stability, dynamics, and conservation of native-like properties of several models of a previously designed amyloid-like fibril of RNase A (Sambashivan et al., Nature 2005; 437:266-269). Through the use of molecular dynamics (MD) simulations, we have provided molecular-level insights into the role of different parts of the sequence on the stability of fibrils, the collective properties of supramolecular complexes, and the presence of native-like conformations and dynamics in supramolecular aggregates. We have been able to show that within the fibrils the three-dimensional globular domain-swapped units preserve the conformational, dynamical, and hydration properties typical of the monomeric state, providing a rationalization for the experimentally observed catalytic activity of fibrils. The nativeness of the globular domains is not affected by the amyloidogenic stretches, which determine the molecular recognition process underlying aggregation through the formation of a stable steric zipper motif. Moreover, through the study of the hydration features of a single sheet model, we have been able to show that polyglutamine stretches of the domain-swapped ribonuclease tend to minimize the interaction with water in favor of sidechain-sidechain interactions, shedding light on the factors leading to the supramolecular assembly of beta-sheet layers into dry steric zippers.

Folding versus aggregation: polypeptide conformations on competing pathways

Protein aggregation has now become recognised as an important and generic aspect of protein energy landscapes. Since the discovery that numerous human diseases are caused by protein aggregation, the biophysical characterisation of misfolded states and their aggregation mechanisms has received increased attention. Utilising experimental techniques and computational approaches established for the analysis of protein folding reactions has ensured rapid advances in the study of pathways leading to amyloid fibrils and amyloid-related aggregates. Here we describe recent experimental and theoretical advances in the elucidation of the conformational properties of dynamic, heterogeneous and/or insoluble protein ensembles populated on complex, multidimensional protein energy landscapes. We discuss current understanding of aggregation mechanisms in this context and describe how the synergy between biochemical, biophysical and cell-biological experiments are beginning to provide detailed insights into the partitioning of non-native species between protein folding and aggregation pathways.

Structural and Pathway Complexity of β-Strand Reorganization within Aggregates of Human Transthyretin(105−115) Peptide

Interstrand conformational rearrangements of human transthyretin peptide (TTR(105-115)) within dimeric aggregates were simulated by means of molecular dynamics (MD) with implicit solvation model for a total length of 48 micros. The conformations sampled in the MD simulations were clustered to identify free energy minima without any projections of free energy surface. A connected graph was constructed with nodes (=clusters) and edges corresponding to free energy minima and transitions between nodes, respectively. This connected graph which reflects the complexity of the free energy surface was used to extract the transition disconnectivity graph, which reflects the overall free energy barriers between pairs of free energy minima but does not contain information on transition paths. The routes of transitions between important free energy minima were obtained by further processing the original graph and the MD data. We have found that both parallel and antiparallel aggregates are populated. The parallel aggregates with different alignment patterns are separated by nonnegligible free energy barriers. Multiroutes exist in the interstrand conformational reorganization. Most visited routes do not dominant the kinetics, while less visited routes contribute a little each but they are numerous and their total contributions are actually dominant. There are various kinds of reptation motions, including those through a beta-bulge, side-chain aided reptation, and flipping or rotation of a hairpin formed by one strand.

Changing the charge distribution of beta-helical-based nanostructures can provide the conditions for charge transfer

In this work we present a computational approach to the design of nanostructures made of structural motifs taken from left-handed beta-helical proteins. Previously, we suggested a structural model based on the self-assembly of motifs taken from Escherichia coli galactoside acetyltransferase (Protein Data Bank 1krr, chain A, residues 131-165, denoted krr1), which produced a very stable nanotube in molecular dynamics simulations. Here we modify this model by changing the charge distribution in the inner core of the system and testing the effect of this change on the structural arrangement of the construct. Our results demonstrate that it is possible to generate the proper conditions for charge transfer inside nanotubes based on assemblies of krr1 segment. The electronic transfer would be achieved by introducing different histidine ionization states in selected positions of the internal core of the construct, in addition to specific mutations with charged amino acids that altogether will allow the formation of coherent networks of aromatic ring stacking, salt-bridges, and hydrogen bonds.

Peptide self-assembly at the nanoscale: a challenging target for computational and experimental biotechnology

Self-assembly at the nanoscale is becoming increasingly important for the fabrication of novel supramolecular structures, with applications in the fields of nanobiotechnology and nanomedicine. Peptides represent the most favorable building blocks for the design and synthesis of nanostructures because they offer a great diversity of chemical and physical properties, they can be synthesized in large amounts, and can be modified and decorated with functional elements, which can be used in diverse applications. In this article, we review some of the most recent experimental advances in the use of nanoscale self-assembled peptide structures and the theoretical efforts aimed at understanding the microscopic determinants of their formation, stability and conformational properties. The combination of experimental observations and theoretical advances will be fundamental to fully realizing the biotechnological potential of peptide self-organization.

Use of constrained synthetic amino acids in beta-helix proteins for conformational control

A highly constrained amino acid has been introduced in the turn region of a beta-helix to increase the conformational stability of the native fold for nanotechnological purposes. The influence of this specific amino acid replacement in the final organization of beta-helix motifs has been evaluated by combining ab initio first-principles calculations on model systems and molecular dynamics simulations of entire peptide segments. The former methodology, which has been applied to a sequence containing three amino acids, has been used to develop adjusted templates. Calculations indicated that 1-amino-2,2-diphenylcyclopropanecarboxylic acid, a constrained cyclopropane analogue of phenylalanine, exhibits a strong tendency to form and promote folded conformations. On the other hand, molecular dynamics simulations are employed to probe the ability of such a synthetic amino acid to enhance the conformational stability of the beta-helix motif, which is the first requirement for further protein nanoengineering. A highly regular segment from a naturally occurring beta-helix protein was selected as a potential nanoconstruct module. Simulations of wild type and mutated segments revealed that the ability of the phenylalanine analogue to nucleate turn conformations enhances the conformational stability of the beta-helix motif in isolated peptide segments.

The molecular dynamics of assembly of the ubiquitous aortic medial amyloidal medin fragment

In recent years there is an increased understanding of the molecular conformation of amyloid fibrils. However, much less is known about the early events that lead to the formation of these medically important assemblies. The clarification of these very important mechanistic details on the process may indicate directions towards the inhibition of the early stages of the assembly, where harmful species are most likely to form. Here, we study the dynamics of assembly of short amyloidogenic peptide fragments from the medin polypeptide. This polypeptide is of unique interest since amyloid deposits composed of medin are found almost in all the population above the age of 50. Twelve independent 50 ns long molecular dynamics simulations in explicit water have been run on peptide NH2-NFGSVQFV-COOH, the minimal recognition hexapeptide element, NH2-NFGSVQ-COOH, and several single-point mutants. In all cases a three-stranded polymeric beta-sheet was used as the basic unit from which fibrils can be formed. Our results clearly indicate the need of well-defined sequence and stereochemical constraints to allow the formation of stable well-ordered aggregates. One of the key findings is the need for the presence of a phenylalanine residue, but not other hydrophobic amino acids, in specific positions within the peptide. Taken together, the results are consistent with recent high-resolution structures of amyloid assemblies and provide unique insights into the dynamics of these structures.

Probing amyloid fibril formation of the NFGAIL peptide by computer simulations

Amyloid fibril formation, as observed in Alzheimer's disease and type II diabetes, is currently described by a nucleation-condensation mechanism, but the details of the process preceding the formation of the nucleus are still lacking. In this study, using an activation-relaxation technique coupled to a generic energy model, we explore the aggregation pathways of 12 chains of the hexapeptide NFGAIL. The simulations show, starting from a preformed parallel dimer and ten disordered chains, that the peptides form essentially amorphous oligomers or more rarely ordered beta-sheet structures where the peptides adopt a parallel orientation within the sheets. Comparison between the simulations indicates that a dimer is not a sufficient seed for avoiding amorphous aggregates and that there is a critical threshold in the number of connections between the chains above which exploration of amorphous aggregates is preferred.

Aggregation of small peptides studied by molecular dynamics simulations

Peptides and proteins tend to aggregate under appropriate conditions. The amyloid fibrils that are ubiquitously found among these structures are associated with major human diseases like Alzheimer's disease, type II diabetes, and various prion diseases. Lately, it has been observed that even very short peptides like tetra and pentapeptides can form ordered amyloid structures. Here, we present aggregation studies of three such small polypeptide systems, namely, the two amyloidogenic peptides DFNKF and FF, and a control (nonamyloidogenic) one, the AGAIL. The respective aggregation process is studied by all-atom Molecular Dynamics simulations, which allow to shed light on the fine details of the association and aggregation process. Our analysis suggests that naturally aggregating systems exhibit significantly diverse overall cluster shape properties and specific intermolecular interactions. Additional analysis was also performed on the previously studied NFGAIL system.

Stability and structure of oligomers of the Alzheimer peptide Abeta16-22: from the dimer to the 32-mer

Several neurodegenerative diseases such as Alzheimer's, Parkinson's, and Huntington's diseases are associated with amyloid fibrils formed by different polypeptides. We probe the structure and stability of oligomers of different sizes of the fragment Abeta(16-22) of the Alzheimer beta-amyloid peptide using atomic-detail molecular dynamics simulations with explicit solvent. We find that only large oligomers form a stable beta-sheet aggregate, the minimum nucleus size being of the order of 8-16 peptides. This effect is attributed to better hydrophobic contacts and a better shielding of backbone-backbone hydrogen bonds from the solvent in bigger assemblies. Moreover, the observed stability of beta-sheet aggregates with a different number of layers can be explained on the basis of their solvent-accessible surface area. Depending on the stacking interface between the sheets, we observe straight or twisted structures, which could be linked to the experimentally observed polymorphism of amyloid fibrils. To compare our 32-mer structure to experimental data, we calculate its x-ray diffraction pattern. Good agreement is found between experimentally and theoretically determined reflections, suggesting that our model indeed closely resembles the structures found in vitro.