Amyloids are protein aggregates involved in various protein condensation diseases. Our study aims to investigate the influence of glycine on the fibrillization mechanism of -lactoglobulin (BLG), a model protein known to form amyloid fibrils from hydrolysed peptides in low pH aqueous solutions. We conducted atomistic molecular dynamics simulations of aqueous solutions of native and unfolded BLG in glycine buffer at pH 2.0. During the simulations we put our focus on analysing protein-protein/buffer interactions, structural electrostatic potential mapping, and the residence times of glycine and glycinium near specific amino acid residues. Glycinium cations were found to preferentially interact with specific protein residues potentially masking the outer disulfide bonds, affecting thiol deprotonation and influencing disulfide scrambling equilibrium. These interactions can potentially hinder hydrolysis and change the fibrillization pathway. Further investigations, such as constant pH MD simulations, simulations on disulfide bounded oligomers are warranted to validate these findings and deepen our understanding of protein aggregation mechanisms.