Using Simulation to Understand the Role of Titration on the Stability of a Peptide-Lipid Bilayer Complex

The pH-low insertion peptide (pHLIP) is an anionic membrane-active peptide with promising potential for applications in imaging of cancer tumors and targeted delivery of chemotherapeutics. The key advantage of pHLIP lies in its acid sensitivity: in acidic cellular environments, pHLIP can insert unidirectionally into the plasma membrane. Partitioning-folding coupling is triggered by titration of the acidic residues in pHLIP, transforming pHLIP from a hydrophilic to a hydrophobic peptide. Despite this knowledge, the reverse pathway that leads to exit of the peptide from the plasma membrane is poorly understood. Our hypothesis is that sequential deprotonation of pHLIP is a prerequisite for exit of the peptide from the plasma membrane. We carried out molecular dynamics (MD) simulations to characterize the effect that deprotonation of the acidic residues of pHLIP has on the stability of the peptide when inserted into a model lipid bilayer of 1-palmitoyl-2-oleoyl-sn-3-phosphocholine (POPC). Initiation of the exit mechanism is facilitated by a complex relationship between the peptide, bulk solvent, and the membrane environment. As the N-terminal acidic residues of pHLIP are deprotonated, localized loss of helicity drives unfolding of the peptide and more pronounced interactions with the bilayer at the lipid-water interface. Deprotonation of the C-terminal acidic residues (D25, D31, D33, and E34) leads to further loss of secondary structure distal from the C-terminus, as well as formation of a water channel that stabilizes the orientation of pHLIP parallel to the membrane normal. Together, these results help explain how stabilization of intermediates between the surface-bound and inserted states of pHLIP occur and provide insights into rational design of pHLIP variants with modified abilities of insertion.