The energetic cost of burying charged groups in the hydrophobic core of lipid bilayers has been controversial, with simulations giving higher estimates than certain experiments. Implicit membrane approaches are usually deemed too simplistic for this problem. Here we challenge this view. The free energy of transfer of amino acid side chains from water to the membrane center predicted by IMM1 is reasonably close to all-atom free energy calculations. The shape of the free energy profile, however, for the charged side chains needs to be modified to reflect the all-atom simulation findings (IMM1-LF). Membrane thinning is treated by combining simulations at different membrane widths with an estimate of membrane deformation free energy from elasticity theory. This approach is first tested on the voltage sensor and the isolated S4 helix of potassium channels. The voltage sensor is stably inserted in a transmembrane orientation for both the original and the modified model. The transmembrane orientation of the isolated S4 helix is unstable in the original model, but a stable local minimum in IMM1-LF, slightly higher in energy than the interfacial orientation. Peptide translocation is addressed by mapping the effective energy of the peptide as a function of vertical position and tilt angle, which allows identification of minimum energy pathways and transition states. The barriers computed for the S4 helix and other experimentally studied peptides are low enough for an observable rate. Thus, computational results and experimental studies on the membrane burial of peptide charged groups appear to be consistent.

Pore formation by membrane-active peptides, naturally encountered in innate immunity and infection, could have important medical and technological applications. Recently, the well-studied lytic peptide melittin has formed the basis for the development of combinatorial libraries from which potent pore-forming peptides have been derived, optimized to work under different conditions. We investigate three such peptides, macrolittin70, which is most active at neutral pH; pHD15, which is active only at low pH; and MelP5_D6, which was rationally designed to be active at low pH but formed only small pores. There are three, six, and six acidic residues in macrolittin70, pHD15, and MelP5_D6, respectively. We perform multi-microsecond simulations in 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) of hexamers of these peptides starting from transmembrane orientations at neutral pH (all residues at standard protonation), low pH (acidic residues and His protonated), and highly acidic environments in which C-termini are also protonated. Previous simulations of the parent peptides melittin and MelP5 in 1,2-dimyristoylsn-glycero-3-phosphocholine (DMPC) are repeated in POPC. We find that the most potent pore-forming peptides exhibit strong interpeptide interactions, including salt bridges, H-bonds, and polar interactions. Protonation of the C-terminus promotes helicity and pore size. The proximity of the peptides allows fewer lipid headgroups to line the pores than in previous simulations, making the pores intermediate between barrel stave and toroidal. Based on these structures and geometrical arguments, we attempt to rationalize the factors that under different conditions can increase or decrease pore stability and propose mutations that could be tested experimentally.

No abstract