Quantitative Characterization of Protein-Lipid Interactions by Free Energy Simulation between Binary Bilayers
Soohyung Park1, Min Sun Yeom2, Olaf S Andersen3
1Departments of Biological Sciences and Bioengineering , Lehigh University , Bethlehem , Pennsylvania 18015 , United States.
Journal of Chemical Theory and Computation
|September 28, 2019
Summary
This study introduces a binary bilayer system (BBS) to accurately measure protein-lipid interactions using molecular dynamics (MD) simulations. The novel method reveals gramicidin A
Area of Science:
- Biophysics
- Computational Biology
- Membrane Protein Dynamics
Background:
- Understanding protein-lipid interactions is crucial for membrane biology.
- Traditional models often simplify the complex environment of lipid bilayers.
- Accurate characterization of protein behavior within membranes requires advanced simulation techniques.
Purpose of the Study:
- To quantitatively characterize protein-lipid interactions using a novel binary bilayer system (BBS).
- To determine the free energy difference for transferring a model protein (gramicidin A channel) between lipid bilayers of different thicknesses.
- To reconcile discrepancies between simulation results and continuum elastic models for membrane protein behavior.
Main Methods:
- Development and application of a binary bilayer system (BBS) comprising DLPC and DMPC lipid bilayers.
- Utilizing umbrella sampling molecular dynamics (MD) simulations for free energy calculations.
- Embedding a gramicidin A (gA) channel within the BBS to study its interactions.
Main Results:
- Calculated a free energy difference of -2.2 ± 0.7 kcal/mol for gA transfer from DLPC to DMPC bilayers.
- Determined an effective hydrophobic length of ~26 Å for the gA channel, resolving discrepancies with previous models.
- Identified an interfacial barrier in the free energy profile due to line tension between bilayers of different thicknesses.
Conclusions:
- The BBS provides an efficient and accurate method for all-atom resolution studies of protein-lipid interactions.
- The effective hydrophobic length of gA is longer than previously assumed, aligning simulation results with continuum models.
- Interfacial barriers play a significant role in peptide translocation between different membrane regions, impacting cellular processes.


