Aggregation of the HIV-1 Gag protein onto the plasma membrane (PM) enables viral budding and infection propagation. Gag assembly at the membrane interface is mediated by its matrix domain (MA), the Myristoylated (Myr) N-terminus. MA targets the PM through electrostatic interactions, mainly at its highly-basic-region (HBR). The mechanism of Myr insertion and its role in proteinmembrane dynamics remains unclear. Using all-atom molecular dynamics, we examined an MA unit in the vicinity of lipid bilayers that model different characteristics of the PM. Interaction with PIP2 and PS lipids is highly favored around the HBR, and is enough to keep the protein bound. Additionally, we simulated three MA units near our bilayers and quantified the collective effects of free monomers vs. formed trimers on Myr insertion events. Micro-second-long trajectories allowed us to observe Myr insertion, propose a mechanism, quantify specific interactions with lipids, and examine the response of the local membrane environment.