Translocation thermodynamics of linear and cyclic nonaarginine into model DPPC bilayer via coarse-grained molecular dynamics simulation: implications of pore formation and nonadditivity

Structural mechanisms and underlying thermodynamic determinants of efficient internalization of charged cationic Peptides (cell-penetrating Peptides, CPPs) such as TAT, polyarginine, and their variants, into cells, cellular constructs, and model membrane/lipid bilayers (large and giant unilamellar or multilamelar vesicles) continue to garner significant attention. Two widely held views on the translocation mechanism center on endocytotic and nonendocytotic (diffusive) processes. Espousing the view of a purely diffusive internalization process (supported by recent experimental evidence, [Säälik, P.; et al. J. Controlled Release 2011, 153, 117-125]), we consider the underlying free energetics of the translocation of a nonaarginine peptide (Arg9) into a model DPPC bilayer. In the case of the Arg9 cationic peptide, recent experiments indicate a higher internalization efficiency of the cyclic structure (cyclic Arg9) relative to the linear conformer. Furthermore, recent all-atom resolution molecular dynamics simulations of cyclic Arg9 [Huang, K.; et al. Biophys. J., 2013, 104, 412-420] suggested a critical stabilizing role of water- and lipid-constituted pores that form within the bilayer as the charged Arg9 translocates deep into the bilayer center. Herein, we use umbrella sampling molecular dynamics simulations with coarse-grained Martini lipids, polarizable coarse-grained water, and peptide to explore the dependence of translocation free energetics on peptide structure and conformation via calculation of potentials of mean force along preselected reaction paths allowing and preventing membrane deformations that lead to pore formation. Within the context of the coarse-grained force fields we employ, we observe significant barriers for Arg9 translocation from bulk aqueous solution to bilayer center. Moreover, we do not find free-energy minima in the headgroup-water interfacial region, as observed in simulations using all-atom force fields. The pore-forming paths systematically predict lower free-energy barriers (ca. 90 kJ/mol lower) than the non pore-forming paths, again consistent with all-atom force field simulations. The current force field suggests no preference for the more compact or covalently cyclic structures upon entering the bilayer. Decomposition of the PMF into the system's components indicates that the dominant stabilizing contribution along the pore-forming path originates from the membrane as both layers of it deformed due to the formation of pore. Furthermore, our analysis revealed that although there is significant entropic stabilization arising from the enhanced configurational entropy exposing more states as the peptide moves through the bilayer, the enthalpic loss (as predicted by the interactions of this coarse-grained model) far outweighs any former stabilization, thus leading to significant barrier to translocation. Finally, we observe reduction in the translocation free-energy barrier for a second Arg9 entering the bilayer in the presence of an initial peptide restrained at the center, again, in qualitative agreement with all-atom force fields.