Using all-atom simulations in explicit solvent to study aggregation of amphipathic peptides into amyloid-like fibrils

Thermal unfolding of alpha-lactalbumin at acidic pH: Insights from molecular dynamics simulations

The transformation of globular proteins into fibrils passes through several stages, including the formation of partially expanded conformational states different from the native or fully unfolded forms. Here we used molecular dynamics simulations to characterize the thermal unfolding of alpha-lactalbumin on the microsecond timescale in the range of temperatures of 300-440 K. Comparative analysis of structural changes, mobility of different parts of protein, and pathways through the free energy landscape during the unfolding of alpha-lactalbumin at different temperatures reveals the existence of several intermediate states separated by small energy barriers. The lifetime of these intermediates depends on temperature and varies from nanoseconds to microseconds.

Pore Formation by Amyloid-like Peptides: Effects of the Nonpolar-Polar Sequence Pattern

One of the mechanisms accounting for the toxicity of amyloid peptides in diseases like Alzheimer's and Parkinson's is the formation of pores on the plasma membrane of neurons. Here, we perform unbiased all-atom simulations of the full membrane damaging pathway, which includes adsorption, aggregation, and perforation of the lipid bilayer accounting for pore-like structures. Simulations are performed using four peptides made with the same amino acids. Differences in the nonpolar-polar sequence pattern of these peptides prompt them to adsorb into the membrane with the extended conformations oriented either parallel [peptide labeled F1, Ac-(FKFE)2-NH2], perpendicular (F4, Ac-FFFFKKEE-NH2), or with an intermediate orientation (F2, Ac-FFKKFFEE-NH2, and F3, Ac-FFFKFEKE-NH2) in regard to the membrane surface. At the water-lipid interface, only F1 fully self-assembles into β-sheets, and F2 peptides partially fold into an α-helical structure. The β-sheets of F1 emerge as electrostatic interactions attract neighboring peptides to intermediate distances where nonpolar side chains can interact within the dry core of the bilayer. This complex interplay between electrostatic and nonpolar interactions is not observed for the other peptides. Although β-sheets of F1 peptides are mostly parallel to the membrane, some of their edges penetrate deep inside the bilayer, dragging water molecules with them. This precedes pore formation, which starts with the flow of two water layers through the membrane that expand into a stable cylindrical pore delimited by polar faces of β-sheets spanning both leaflets of the bilayer.

Peptide Self-Assembly into Amyloid Fibrils: Unbiased All-Atom Simulations

Protein self-assembly plays an important role in biological systems, accounting for the formation of mesoscopic structures that can be highly symmetric as in the capsid of viruses or disordered as in molecular condensates or exhibit a one-dimensional fibrillar morphology as in amyloid fibrils. Deposits of the latter in tissues of individuals with degenerative diseases like Alzheimer's and Parkinson's has motivated extensive efforts to understand the sequence of molecular events accounting for their formation. These studies aim to identify on-pathway intermediates that may be the targets for therapeutic intervention. This detailed knowledge of fibril formation remains obscure, in part due to challenges with experimental analyses of these processes. However, important progress is being achieved for short amyloid peptides due to advances in our ability to perform completely unbiased all-atom simulations of the self-assembly process. This perspective discusses recent developments, their implications, and the hurdles that still need to be overcome to further advance the field.

Resolving the Nanoscale Structure of β-Sheet Peptide Self-Assemblies Using Single-Molecule Orientation-Localization Microscopy

Synthetic peptides that self-assemble into cross-β fibrils are versatile building blocks for engineered biomaterials due to their modularity and biocompatibility, but their structural and morphological similarities to amyloid species have been a long-standing concern for their translation. Further, their polymorphs are difficult to characterize by using spectroscopic and imaging techniques that rely on ensemble averaging to achieve high resolution. Here, we utilize Nile red (NR), an amyloidophilic fluorogenic probe, and single-molecule orientation-localization microscopy (SMOLM) to characterize fibrils formed by the designed amphipathic enantiomers KFE8L and KFE8D and the pathological amyloid-beta peptide Aβ42. Importantly, NR SMOLM reveals the helical (bilayer) ribbon structure of both KFE8 and Aβ42 and quantifies the precise tilt of the fibrils' inner and outer backbones in relevant buffer conditions without the need for covalent labeling or sequence mutations. SMOLM also distinguishes polymorphic branched and curved morphologies of KFE8, whose backbones exhibit much more heterogeneity than those of typical straight fibrils. Thus, SMOLM is a powerful tool to interrogate the structural differences and polymorphism between engineered and pathological cross-β-rich fibrils.

Nucleation and Growth of Amyloid Fibrils

The formation of amyloid fibrils is a complex phenomenon that remains poorly understood at the atomic scale. Herein, we perform extended unbiased all-atom simulations in explicit solvent of a short amphipathic peptide to shed light on the three mechanisms accounting for fibril formation, namely, nucleation via primary and secondary mechanisms, and fibril growth. We find that primary nucleation takes place via the formation of an intermediate state made of two laminated β-sheets oriented perpendicular to each other. The amyloid fibril spine subsequently emerges from the rotation of these β-sheets to account for peptides that are parallel to each other and perpendicular to the main axis of the fibril. Growth of this spine, in turn, takes place via a dock-and-lock mechanism. We find that peptides dock onto the fibril tip either from bulk solution or after diffusing on the fibril surface. The latter docking pathway contributes significantly to populate the fibril tip with peptides. We also find that side chain interactions drive the motion of peptides in the lock phase during growth, enabling them to adopt the structure imposed by the fibril tip with atomic fidelity. Conversely, the docked peptide becomes trapped in a local free energy minimum when docked-conformations are sampled randomly. Our simulations also highlight the role played by nonpolar fibril surface patches in catalyzing and orienting the formation of small cross-β structures. More broadly, our simulations provide important new insights into the pathways and interactions accounting for primary and secondary nucleation as well as the growth of amyloid fibrils.

Effects of ion type and concentration on the structure and aggregation of the amyloid peptide A β 16 - 22 $$ {\boldsymbol{\beta}}_{16-22} $$

Among the various factors controlling the amyloid aggregation process, the influences of ions on the aggregation rate and the resulting structures are important aspects to consider, which can be studied by molecular simulations. There is a wide variety of protein force fields and ion models, raising the question of which model to use in such studies. To address this question, we perform molecular dynamics simulations of Aβ16-22 , a fragment of the Alzheimer's amyloid β peptide, using different protein force fields, AMBER99SB-disp (A99-d) and CHARMM36m (C36m), and different ion parameters. The influences of NaCl and CaCl2 at various concentrations are studied and compared with the systems without the addition of ions. Our results indicate a sensitivity of the peptide-ion interactions to the different ion models. In particular, we observe a strong binding of Ca2+ to residue E22 with C36m and also with the Åqvist ion model used together with A99-d, which slightly affects the monomeric Aβ16-22 structures and the aggregation rate, but significantly affects the oligomer structures formed in the aggregation simulations. For example, at high Ca2+ concentrations, there was a switch from an antiparallel to a parallel β-sheet. Such ionic influences are of biological relevance because local ion concentrations can change in vivo and could help explain the polymorphism of amyloid fibrils.

Resolving the nanoscale structure of β-sheet assemblies using single-molecule orientation-localization microscopy

Synthetic peptides that self-assemble into cross-β fibrils have remarkable utility as engineered biomaterials due to their modularity and biocompatibility, but their structural and morphological similarity to amyloid species has been a long-standing concern for their translation. Further, their polymorphs are difficult to characterize using spectroscopic and imaging techniques that rely on ensemble averaging to achieve high resolution. Here, we utilize single-molecule orientation-localization microscopy (SMOLM) to characterize fibrils formed by the designed amphipathic enantiomers, KFE8L and KFE8D, and the pathological amyloid-beta peptide Aβ42. SMOLM reveals that the orientations of Nile red, as it transiently binds to both KFE8 and Aβ42, are consistent with a helical (bilayer) ribbon structure and convey the precise tilt of the fibrils' inner and outer backbones. SMOLM also finds polymorphic branched and curved morphologies of KFE8 whose backbones exhibit much more heterogeneity than those of more typical straight fibrils. Thus, SMOLM is a powerful tool to interrogate the structural differences and polymorphism between engineered and pathological cross β-rich fibrils.

Aggregation of Amyloidogenic Peptide Uperin-Molecular Dynamics Simulations

Uperin 3.5 is a remarkable natural peptide obtained from the skin of toadlets comprised of 17 amino acids which exhibits both antimicrobial and amyloidogenic properties. Molecular dynamics simulations were performed to study the β-aggregation process of uperin 3.5 as well as two of its mutants, in which the positively charged residues Arg7 and Lys8 have been replaced by alanine. All three peptides rapidly underwent spontaneous aggregation and conformational transition from random coils to beta-rich structures. The simulations reveal that the initial and essential step of the aggregation process involves peptide dimerization and the formation of small beta-sheets. A decrease in positive charge and an increase in the number of hydrophobic residues in the mutant peptides lead to an increase in the rate of their aggregation.

Peptide-Membrane Binding: Effects of the Amino Acid Sequence

An understanding of how the amino acid sequence affects the interaction of peptides with lipid membranes remains mostly unknown. This type of knowledge is required to rationalize membrane-induced toxicity of amyloid peptides and to design peptides that can interact with lipid bilayers. Here, we perform a systematic study of how variations in the sequence of the amphipathic Ac-(FKFE)₂-NH₂ peptide affect its interaction with zwitterionic lipid bilayers using extensive all-atom molecular dynamics simulations in explicit solvent. Our results show that peptides with a net positive charge bind more frequently to the lipid bilayer than neutral or negatively charged sequences. Moreover, neutral amphipathic peptides made with the same numbers of phenylalanine (F), lysine (K), and glutamic (E) amino acids at different positions in the sequence differ significantly in their frequency of binding to the membrane. We find that peptides bind with a higher frequency to the membrane if their positive lysine side chains are more exposed to the solvent, which occurs if they are located at the extremity (as opposed to the middle) of the sequence. Non-polar residues play an important role in accounting for the adsorption of peptides onto the membrane. In particular, peptides made with less hydrophobic non-polar residues (e.g., valine and alanine) are significantly less adsorbed to the membrane compared to peptides made with phenylalanine. We also find that sequences where phenylalanine residues are located at the extremities of the peptide have a higher tendency to be adsorbed.

Atomic Insights into Amyloid-Induced Membrane Damage

Amphipathic peptides can cause biological membranes to leak either by dissolving their lipid content via a detergent-like mechanism or by forming pores on the membrane surface. These modes of membrane damage have been related to the toxicity of amyloid peptides and to the activity of antimicrobial peptides. Here, we perform the first all-atom simulations in which membrane-bound amphipathic peptides self-assemble into β-sheets that subsequently either form stable pores inside the bilayer or drag lipids out of the membrane surface. An analysis of these simulations shows that the acyl tail of lipids interact strongly with non-polar side chains of peptides deposited on the membrane. These strong interactions enable lipids to be dragged out of the bilayer by oligomeric structures accounting for detergent-like damage. They also disturb the orientation of lipid tails in the vicinity of peptides. These distortions are minimized around pore structures. We also show that membrane-bound β-sheets become twisted with one of their extremities partially penetrating the lipid bilayer. This allows peptides on opposite leaflets to interact and form a long transmembrane β-sheet, which initiates poration. In simulations, where peptides are deposited on a single leaflet, the twist in β-sheets allows them to penetrate the membrane and form pores. In addition, our simulations show that fibril-like structures produce little damage to lipid membranes, as non-polar side chains in these structures are unavailable to interact with the acyl tail of lipids.

β-Sheet to Random Coil Transition in Self-Assembling Peptide Scaffolds Promotes Proteolytic Degradation

One of the most desirable properties that biomaterials designed for tissue engineering or drug delivery applications should fulfill is biodegradation and resorption without toxicity. Therefore, there is an increasing interest in the development of biomaterials able to be enzymatically degraded once implanted at the injury site or once delivered to the target organ. In this paper, we demonstrate the protease sensitivity of self-assembling amphiphilic peptides, in particular, RAD16-I (AcN-RADARADARADARADA-CONH2), which contains four potential cleavage sites for trypsin. We detected that when subjected to thermal denaturation, the peptide secondary structure suffers a transition from β-sheet to random coil. We also used Matrix-Assisted Laser Desorption/Ionization-Time-Of-Flight (MALDI-TOF) to detect the proteolytic breakdown products of samples subjected to incubation with trypsin as well as atomic force microscopy (AFM) to visualize the effect of the degradation on the nanofiber scaffold. Interestingly, thermally treated samples had a higher extent of degradation than non-denatured samples, suggesting that the transition from β-sheet to random coil leaves the cleavage sites accessible and susceptible to protease degradation. These results indicate that the self-assembling peptide can be reduced to short peptide sequences and, subsequently, degraded to single amino acids, constituting a group of naturally biodegradable materials optimal for their application in tissue engineering and regenerative medicine.

Peptide Self-Assembled Nanostructures: From Models to Therapeutic Peptides

Self-assembly is the most suitable approach to obtaining peptide-based materials on the nano- and mesoscopic scales. Applications span from peptide drugs for personalized therapy to light harvesting and electron conductive media for solar energy production and bioelectronics, respectively. In this study, we will discuss the self-assembly of selected model and bioactive peptides, in particular reviewing our recent work on the formation of peptide architectures of nano- and mesoscopic size in solution and on solid substrates. The hierarchical and cooperative characters of peptide self-assembly will be highlighted, focusing on the structural and dynamical properties of the peptide building blocks and on the nature of the intermolecular interactions driving the aggregation phenomena in a given environment. These results will pave the way for the understanding of the still-debated mechanism of action of an antimicrobial peptide (trichogin GA IV) and the pharmacokinetic properties of a peptide drug (semaglutide) currently in use for the therapy of type-II diabetes.