Department of Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, Minnesota 55455-0132, USA.
J Phys Chem B. 2010 Mar 4;114(8):2790-7. doi: 10.1021/jp909640g.
The free energies of adsorption of the monomer or dimer of the cationic beta-hairpin antimicrobial peptide protegrin-1 (PG1) in a specific binding orientation on a lipid bilayer are determined using molecular dynamics (MD) simulations and Poisson-Boltzmann calculations. The bilayer is composed of anionic palmitoyl-oleoyl-phosphatidylglycerol (POPG) and palmitoyl-oleoyl-phosphatidylethanolamine (POPE) with ratio 1:3 (POPG/POPE). PG1 is believed to kill bacteria by binding on their membranes. There, it forms pores that lyse the bacteria. Herein we focus on the thermodynamics of binding. In particular, we explore the role of counterion release from the lipid bilayer upon adsorption of either the monomeric or the dimeric form of PG1. Twenty-two 4-ns-long MD trajectories of equilibrated systems are generated to determine the free energy profiles for the monomer and dimer as a function of the distance between the peptide(s) and the membrane surface. The MD simulations are conducted at 11 different separations from the membrane for each of the two systems, one with PG1, the second with a PG1 dimer of only a specific orientation of the monomer and dimer without taking into account the change of entropy for the peptide. To calculate the potential of mean force for each peptide/membrane system, a variant of constrained MD and thermodynamic integration is used. We observed that PG1 dimer binds more favorably to the POPG/POPE membrane. A simple method for relating the free energy profile to the PG1-membrane binding constant is employed to predict a free energy of adsorption of -2.4 +/- 0.8 kcal/mol. A corresponding PG1-dimer-membrane binding constant is calculated as -3.5 +/- 1.1 kcal/mol. Free energy profiles from MD simulation were extensively analyzed and compared with results of Poisson-Boltzmann theory. We find the peptide-membrane attraction to be dominated by the entropy increase due to the release of counterions in a POPG/POPE lipid bilayer.
采用分子动力学(MD)模拟和泊松-玻尔兹曼计算,确定阳离子β发夹抗菌肽 PG1 的单体或二聚体在特定结合方向上在脂质双层上的吸附自由能。双层由阴离子棕榈酰-油酰-磷脂酰甘油(POPG)和棕榈酰-油酰-磷脂酰乙醇胺(POPE)组成,比例为 1:3(POPG/POPE)。PG1 被认为通过结合在其膜上而杀死细菌。在那里,它形成孔,使细菌裂解。本文关注的是结合的热力学。特别是,我们探讨了在吸附单体或二聚体形式的 PG1 时,来自脂质双层的抗衡离子释放对吸附的作用。生成了 22 个 4-ns 长的平衡系统 MD 轨迹,以确定单体和二聚体的自由能曲线作为肽与膜表面之间距离的函数。对于每个系统,在 11 个不同的膜分离处进行 MD 模拟,一个系统为 PG1,第二个系统为 PG1 二聚体,仅考虑单体和二聚体的特定取向,而不考虑肽的熵变。为了计算每个肽/膜系统的平均力势能,使用了约束 MD 和热力学积分的变体。我们观察到 PG1 二聚体更有利于与 POPG/POPE 膜结合。采用一种将自由能曲线与 PG1-膜结合常数相关联的简单方法来预测吸附自由能为-2.4 +/- 0.8 kcal/mol。相应的 PG1-二聚体-膜结合常数计算为-3.5 +/- 1.1 kcal/mol。对 MD 模拟的自由能曲线进行了广泛的分析,并与泊松-玻尔兹曼理论的结果进行了比较。我们发现,在 POPG/POPE 脂质双层中,由于抗衡离子的释放,肽-膜的吸引力主要由熵的增加所主导。