Neurodegenerative diseases, such as Alzheimer's, Parkinson's, and Huntington's, affect millions worldwide and share a common feature: the aggregation of intrinsically disordered proteins into toxic oligomers that interact with cell membranes. In Alzheimer's disease (AD), amyloid-beta (Aβ) peptides accumulate and bind to plasma membranes, potentially disrupting cellular function. The complex interplay between amyloidogenic peptides and lipid membranes, particularly the role of anionic lipids, is crucial in disease pathogenesis but challenging to characterize experimentally. The literature presents conflicting results on the influence of anionic lipids on peptide aggregation kinetics, highlighting a knowledge gap. To address this, we used coarse-grained molecular dynamics (CG-MD) simulations to study interactions between a model amyloidogenic peptide, amyloid-β's K16LVFFAE22 fragment (Aβ16-22), and mixed lipid bilayers. We used phosphatidylserine (PS) and phosphatidylcholine (PC) as representative anionic and zwitterionic lipids, respectively, examining the mixed bilayer compositions of 0% PS-100% PC, 10% PS-90% PC, and 30% PS-70% PC. Our simulations revealed that membranes enriched in anionic lipids enhance peptide adsorption and interaction kinetics. The aggregation dynamics was modulated by two competing factors: increased local peptide concentration near negatively charged membranes, which promoted aggregation, and peptide-lipid interactions, which slowed it down. Higher percentages of anionic lipids led to smaller and more ordered aggregates and enhanced lipid demixing, leading to the formation of PS clusters. These findings contribute to understanding membrane-mediated peptide aggregation in neurodegenerative disorders, potentially guiding new therapeutic strategies targeting the early stages of protein aggregation in various neurodegenerative diseases.