The enhanced cell penetration ability of arginine-rich peptides, such as nonaarginine (R 9 R₉), compared to their lysine-rich counterparts, remains incompletely understood. Atomistic simulations reveal that R 9 R₉ binds significantly stronger (≈ 20 kJ / mol 20 kJ/mol) and penetrates deeper into the anionic lipid headgroup region than its lysine equivalent. This enhanced interaction translates into a stronger induction of negative membrane curvature by R 9 R₉. We introduce an integrative modeling workflow to extract and incorporate material properties from molecular simulations into a continuum membrane model that includes peptide binding and curvature induction. Our model predicts that stable membrane invaginations, as observed in studies of cell penetration, require excess membrane and are stable only for R 9 R₉. By analyzing lipid and protein sorting coupled to the membrane structure, we explain the interplay of Gaussian and mean curvature in providing a mechanistic basis for the initial membrane deformation events potentially involved in "Arginine Magic" cell entry pathways.
Baxová et al. (2026) studied this question.