Investigating Hydrophilic Pores in Model Lipid Bilayers Using Molecular Simulations: Correlating Bilayer Properties with Pore-Formation Thermodynamics

Cell-penetrating and antimicrobial peptides show a remarkable ability to translocate across physiological membranes. Along with factors such as electric-potential-induced perturbations of membrane structure and surface tension effects, experiments invoke porelike membrane configurations during the solute transfer process into vesicles and cells. The initiation and formation of pores are associated with a nontrivial free-energy cost, thus necessitating a consideration of the factors associated with pore formation and the attendant free energies. Because of experimental and modeling challenges related to the long time scales of the translocation process, we use umbrella sampling molecular dynamics simulations with a lipid-density-based order parameter to investigate membrane-pore-formation free energy employing Martini coarse-grained models. We investigate structure and thermodynamic features of the pore in 18 lipids spanning a range of headgroups, charge states, acyl chain lengths, and saturation. We probe the dependence of pore-formation barriers on the area per lipid, lipid bilayer thickness, and membrane bending rigidities in three different lipid classes. The pore-formation free energy in pure bilayers and peptide translocating scenarios are significantly coupled with bilayer thickness. Thicker bilayers require more reversible work to create pores. The pore-formation free energy is higher in peptide lipid systems than in peptide-free lipid systems due to penalties to maintain the solvation of charged hydrophilic solutes within the membrane environment.