Driving β-Strands into 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.

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.

Cold Denaturation of Proteins: Where Bioinformatics Meets Thermodynamics to Offer a Mechanistic Understanding: Pea Protein As a Case Study

Protein structure can be altered with heat, but models which predict denaturation show that globular proteins also spontaneously unfold at low temperatures through cold denaturation. By an analysis of the primary structure of pea protein using bioinformatic modeling, a mechanism of pea protein cold denaturation is proposed. Pea protein is then fractionated into partially purified legumin and vicilin components, suspended in ethanol, and subjected to low temperatures (-10 to -20 °C). The structural characterizations of the purified fractions are conducted through FTIR, ζ potential, dynamic light scattering, and oil binding, and these are compared to the results of commercial protein isolates. The observed structural changes suggest that pea protein undergoes changes in structure as the result of low-temperature treatments, which could lead to innovative industrial processing techniques for functionalization by low-temperature processing.

Intrinsic Conformational Preferences and Interactions in α-Synuclein Fibrils: Insights from Molecular Dynamics Simulations

Amyloid formation by the intrinsically disordered α-synuclein protein is the hallmark of Parkinson's disease. We present atomistic Molecular Dynamics simulations of the core of α-synuclein using enhanced sampling techniques to describe the conformational and binding free energy landscapes of fragments implicated in fibril stabilization. The theoretical framework is derived to combine the free energy profiles of the fragments into the reaction free energy of a protein binding to a fibril. Our study shows that individual fragments in solution have a propensity toward attaining non-β conformations, indicating that in a fibril β-strands are stabilized by interactions with other strands. We show that most dimers of hydrogen-bonded fragments are unstable in solution, while hydrogen bonding stabilizes the collective binding of five fragments to the end of a fibril. Hydrophobic effects make further contributions to the stability of fibrils. This study is the first of its kind where structural and binding preferences of the five major fragments of the hydrophobic core of α-synuclein have been investigated. This approach improves sampling of intrinsically disordered proteins, provides information on the binding mechanism between the core sequences of α-synuclein, and enables the parametrization of coarse grained models.

Dynamics of Seeded Aβ40-Fibril Growth from Atomistic Molecular Dynamics Simulations: Kinetic Trapping and Reduced Water Mobility in the Locking Step

Filamentous β-amyloid aggregates are crucial for the pathology of Alzheimer's disease. Despite the tremendous biomedical importance, the molecular pathway of growth propagation is not completely understood and remains challenging to investigate by simulations due to the long time scales involved. Here, we apply extensive all-atom molecular dynamics simulations in explicit water to obtain free energy profiles and kinetic information from position-dependent diffusion profiles for three different Aβ9-40-growth processes: fibril elongation by single monomers at the structurally unequal filament tips and association of larger filament fragments. Our approach provides insight into the molecular steps of the kinetic pathway and allows close agreement with experimental binding free energies and macroscopic growth rates. Water plays a decisive role, and solvent entropy is identified as the main driving force for assembly. Fibril growth is disfavored energetically due to cancellation of direct peptide-peptide interactions and solvation effects. The kinetics of growth is consistent with the characteristic dock/lock mechanism, and docking is at least 2 orders of magnitude faster. During initial docking, interactions are mediated by transient non-native hydrogen bonds, which efficiently catch the incoming monomer or fragment already at separations of about 3 nm. In subsequent locking, the dynamics is much slower due to formation of kinetically trapped conformations caused by long-lived non-native hydrogen bonds. Fibril growth additionally requires collective motion of water molecules to create a dry binding interface. Fibril growth is further retarded due to reduced mobility of the involved hydration water, evident from a 2-fold reduction of the diffusion coefficient.

Thermodynamics of Aβ16-21 dissociation from a fibril: Enthalpy, entropy, and volumetric properties

Here, we provide insights into the thermodynamic properties of A β16-21 dissociation from an amyloid fibril using all-atom molecular dynamics simulations in explicit water. An umbrella sampling protocol is used to compute potentials of mean force (PMF) as a function of the distance ξ between centers-of-mass of the A β16-21 peptide and the preformed fibril at nine temperatures. Changes in the enthalpy and the entropic energy are determined from the temperature dependence of these PMF(s) and the average volume of the simulation box is computed as a function of ξ. We find that the PMF at 310 K is dominated by enthalpy while the entropic energy does not change significantly during dissociation. The volume of the system decreases during dissociation. Moreover, the magnitude of this volume change also decreases with increasing temperature. By defining dock and lock states using the solvent accessible surface area (SASA), we find that the behavior of the electrostatic energy is different in these two states. It increases (unfavorable) and decreases (favorable) during dissociation in lock and dock states, respectively, while the energy due to Lennard-Jones interactions increases continuously in these states. Our simulations also highlight the importance of hydrophobic interactions in accounting for the stability of A β16-21.

Fourier transform infrared spectroscopy provides an evidence of papain denaturation and aggregation during cold storage

Papain is a cysteine protease with wide substrate specificity and many applications. Despite its widespread applications, cold stability of papain has never been studied. Here, we used differential spectroscopy to monitor thermal denaturation process. Papain was the most stabile from 45 °C to 60 °C with ΔG°321 of 13.9±0.3 kJ/mol and Tm value of 84±1 °C. After cold storage, papain lost parts of its native secondary structures elements which gave an increase of 40% of intermolecular β-sheet content (band maximum detected at frequency of 1621 cm(-1) in Fourier transform infrared (FT-IR) spectrum) indicating the presence of secondary structures necessary for aggregation. The presence of protein aggregates after cold storage was also proven by analytical size exclusion chromatography. After six freeze-thaw cycles around 75% of starting enzyme activity of papain was lost due to cold denaturation and aggregation of unfolded protein. Autoproteolysis of papain did not cause significant loss of the protein activity. Upon the cold storage, papain underwent structural rearrangements and aggregation that correspond to other cold denatured proteins, rather than autoproteolysis which could have the commercial importance for the growing polypeptide based industry.

Role of side-chain interactions on the formation of α-helices in model peptides

The role played by side-chain interactions on the formation of α-helices is studied using extensive all-atom molecular dynamics simulations of polyalanine-like peptides in explicit TIP4P water. The peptide is described by the OPLS-AA force field except for the Lennard-Jones interaction between Cβ-Cβ atoms, which is modified systematically. We identify values of the Lennard-Jones parameter that promote α-helix formation. To rationalize these results, potentials of mean force (PMF) between methane-like molecules that mimic side chains in our polyalanine-like peptides are computed. These PMF exhibit a complex distance dependence where global and local minima are separated by an energy barrier. We show that α-helix propensity correlates with values of these PMF at distances corresponding to Cβ-Cβ of i-i+3 and other nearest neighbors in the α-helix. In particular, the set of Lennard-Jones parameters that promote α-helices is characterized by PMF that exhibit a global minimum at distances corresponding to i-i+3 neighbors in α-helices. Implications of these results are discussed.