Molecular Dynamics Insights into the Aggregation Behavior of N-Terminal β-Lactoglobulin Peptides.

β-lactoglobulin (BLG) forms amyloid-like aggregates at high temperatures, low pH, and low ionic strengths. At a pH below 2, BLG undergoes hydrolysis into peptides, with N-terminal peptides 1-33 and 1-52 being prone to fibrillization, forming amyloid-like fibrils. Due to their good mechanical properties, BLG amyloids demonstrate great potential for diverse applications, including biosensors, nanocomposites, and catalysts. Consequently, further studies are essential to comprehensively understand the factors governing the formation of BLG amyloid-like morphologies. In this study, all-atom molecular dynamics simulations were employed to explore the aggregation of N-terminal 1-33 and 1-52 BLG peptides under conditions of pH 2 and at 10 mM NaCl concentration. The simulations revealed that the peptides spontaneously assembled into aggregates of varying sizes. The aggregation process was enabled by the low charge of peptides and the presence of hydrophobic residues within them. As the peptides associated into aggregates, there was a concurrent increase in β-sheet structures and the establishment of hydrogen bonds, enhancing the stability of the aggregates. Notably, on average, 1-33 peptides formed larger aggregates compared to their 1-52 counterparts, while the latter exhibited a slightly higher content of β-sheets and higher cluster orderliness. The applied approach facilitated insights into the early stages of amyloid-like aggregation and molecular-level insight into the formation of β-sheets, which serve as nucleation points for further fibril growth.