Amino Acid Insertion Energetics in a POPC Bilayer from Unbiased Molecular Dynamics

Membrane association is governed by the thermodynamics of amino acid partitioning between water and the lipid bilayer. Here, we quantified amino acid side-chain insertion energetics in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer using unbiased molecular dynamics simulations. Equilibrium depth distributions of 28 analogs, including multiple protonation states, were converted into potentials of mean force (PMFs) by Boltzmann inversion. The resulting PMFs reproduced the main features of bilayer partitioning. Hydrophobic analogs favored the bilayer core, aromatic analogs were stabilized in interfacial regions, and polar or charged analogs remained unfavorable in the hydrophobic interior. A diglycine analog representing the peptide backbone behaved similarly to uncharged polar residues. Depth-dependent pKa profiles and orientational analyses further showed how protonation equilibria and aromatic-ring alignment influence insertion energetics. Agreement with experimental hydrophobicity scales supports the robustness of the approach. These results provide an efficient and internally consistent framework for characterizing bilayer insertion energetics and establish a reference for future studies in more complex lipid environments. O_FIG O_LINKSMALLFIG WIDTH=198 HEIGHT=200 SRC="FIGDIR/small/723583v1_ufig1.gif" ALT="Figure 1"> View larger version (79K): org.highwire.dtl.DTLVardef@127b12org.highwire.dtl.DTLVardef@14de924org.highwire.dtl.DTLVardef@53b27org.highwire.dtl.DTLVardef@16e8ee4_HPS_FORMAT_FIGEXP M_FIG C_FIG SIGNIFICANCEMembrane-associated proteins represent a large fraction of the proteome and include many major drug targets, yet quantitative understanding of their interactions with lipid bilayers remains limited. Here, we present an unbiased molecular dynamics framework for systematically determining amino acid side-chain insertion free energies in a model bilayer. By deriving potentials of mean force directly from equilibrium depth distributions, this approach enables internally consistent comparisons across residue classes and protonation states without biasing restraints. The resulting free-energy profiles reproduce established hydrophobicity trends and show how protonation equilibria and aromatic-ring orientation modulate bilayer partitioning. This scalable strategy provides a quantitative reference for residue-level membrane thermodynamics and establishes a foundation for extending insertion energetics to more diverse lipid compositions and more complex membrane-associated systems.