Phosphatidylserine and Phosphatidylethanolamine Asymmetry Have a Negligible Effect on the Global Structure, Dynamics, and Interactions of the KRAS Lipid Anchor.

The intrinsically disordered C-terminus of the prominent oncogenic protein KRAS-4B (KRAS) selectively interacts and clusters with phosphatidylserine (PS) lipids in the plasma membrane (PM). This 11-residue segment, called tK, contains a polybasic domain (PBD) of six contiguous lysine residues and a farnesylated cysteine. Previous molecular dynamics (MD) simulation studies of tK in phosphatidylcholine (PC)/PS bilayers have suggested that backbone conformational dynamics modulate tK-PS interactions. These simulations have been conducted in symmetric membranes whereas the PM is compositionally asymmetric, with the inner leaflet, where KRAS is localized, being enriched with PS and phosphatidylethanolamine (PE) lipids. To examine if bilayer asymmetry affects tK conformational dynamics and interaction with lipids, we conducted two 10 μs long MD simulations of tK bound to a PC/PS and a PC/PS/PE bilayer in which the PS and PE lipids are distributed in one leaflet. We found that, first, these compositional asymmetries caused differences in acyl chain dynamics between leaflets, but the equilibrium structural and dynamic properties of the two asymmetric bilayers are similar; second, in both systems tK is highly dynamic and samples at least two distinct conformational states; third, PS-tK hydrogen-bonding interactions vary with peptide backbone conformations, and lysine side chains in the PBD predominantly interact with the serine oxygens of PS. These results are in good agreement with previous observations of tK in symmetric membranes. The effects of POPS asymmetry or the presence of POPE on tK are limited to modulating the relative contribution of individual side chains to interactions with lipids and redistributing conformational substates. Additional observations include the larger flexibility of tK in the current simulations, which we attribute to the longer duration of the simulations and the use of the CHARMM36m force field, which more accurately models intrinsically disordered peptides such as tK.