[论文解读] Response of the hydrophilic part of lipid membranes to changing conditions - a critical comparison of simulations to experiments

O. H. Samuli Ollila|arXiv (Cornell University)|Sep 9, 2013

Lipid Membrane Structure and Behavior参考文献 4被引用 2

一句话总结

本研究通过将14种原子尺度脂质膜模型的模拟结果与完全水化双分子层中甘油和胆碱头基C-H键向量序参数的实验数据进行对比，评估了这些模型的性能。尽管当前模型尚不足以精确解析原子结构，但CHARMM36、GAFFlipid和MacRog模型在二面角分布得到改进后展现出潜力，且所有模型均正确预测了脱水条件下P-N向量的倾斜行为。

ABSTRACT

Phospholipids are essential building blocks of biological membranes. Despite of vast amount of accurate experimental data the atomistic resolution structures sampled by the glycerol backbone and choline headgroup in phoshatidylcholine bilayers are not known. Atomistic resolution molecular dynamics simulation model would automatically resolve the structures giving an interpretation of experimental results, if the model would reproduce the experimental data. In this work we compare the C-H bond vector order parameters for glycerol backbone and choline headgroup between 14 different atomistic resolution models and experiments in fully hydrated lipid bilayer. The current models are not accurately enough to resolve the structure. However, closer inspection of three best performing models (CHARMM36, GAFFlipid and MacRog) suggest that improvements in the sampled dihedral angle distributions would potentilly lead to the model which would resolve the structure. Despite of the inaccuracy in the fully hydrated structures, the response to the dehydration, i.e. P-N vector tilting more parallel to membrane normal, is qualitatively correct in all models. The CHARMM36 and MacRog models describe the interactions between lipids and cholesterol better than Berger/Holtje model. This work has been, and continues to be, progressed and discussed through the blog: this http URL. Everyone is invited to join the discussion and make contributions through the blog. The manuscript will be eventually submitted to an appropriate scientific journal. Everyone who has contributed to the work through the blog will be offered coauthorship. For more details see: this http URL.

研究动机与目标

• 评估14种原子尺度脂质模型在完全水化磷脂酰胆碱双分子层中重现实验C-H键向量序参数的准确性，重点关注甘油和胆碱头基部分。
• 识别哪些模型最能捕捉脂质头基在原子尺度下的结构动力学行为。
• 确定当前模型是否能可靠预测脂质膜在脱水条件下的响应行为。
• 通过揭示二面角分布中的缺陷，为未来模型开发提供指导。
• 通过开放博客形式的协作与共同作者机会，推动社区驱动的模型改进。

提出的方法

• 将14种原子尺度脂质模型的模拟C-H键向量序参数与完全水化磷脂酰胆碱双分子层的实验数据进行对比。
• 分析各模型中甘油骨架和胆碱头基区域的二面角分布。
• 通过测量P-N向量相对于膜法线的倾斜程度，评估模型对脱水的响应。
• 将CHARMM36、MacRog和Berger/Holtje模型中的脂质-胆固醇相互作用与实验趋势进行基准测试。
• 利用公开博客记录、讨论并迭代优化模型性能，整合社区反馈。
• 将社区贡献者纳入最终论文，并提供共同作者资格。

实验结果

研究问题

• RQ1 在完全水化双分子层中，哪些原子尺度脂质模型最准确地再现了甘油和胆碱头基的实验C-H键向量序参数？
• RQ2 当前模型在预测脂质膜脱水响应方面表现如何，特别是P-N向量重新取向的行为？
• RQ3 现有模型在二面角分布上的哪些具体缺陷限制了其对原子尺度头基结构的解析能力？
• RQ4 不同模型在模拟脂质-胆固醇相互作用方面表现如何，特别是在结构有序性方面的差异？
• RQ5 通过开放讨论实现的社区驱动模型优化，能否提升脂质膜模拟的预测准确性？

主要发现

• 在测试的14种原子尺度脂质模型中，没有一种能准确再现完全水化磷脂酰胆碱双分子层中甘油骨架和胆碱头基的原子结构。
• CHARMM36、GAFFlipid和MacRog模型表现最佳，若能改进二面角采样，有望实现原子尺度结构的精确解析。
• 所有模型在定性上均正确再现了脱水条件下的实验响应，包括P-N向量向膜法线方向倾斜的增加。
• CHARMM36和MacRog模型在描述脂质-胆固醇相互作用方面优于Berger/Holtje模型，表明其对膜横向组织的表征更优。
• 本研究明确了未来改进路径：通过优化表现最佳模型的二面角参数，实现原子尺度精度。
• 已建立一种开放、协作的博客工作流，用于迭代改进模型，贡献者可获得共同作者资格。