Department of Chemistry , The University of Western Ontario , London , Ontario , Canada N6A 3K7.
Departamento de Química, FFCLRP , Universidade de São Paulo , Avenida Bandeirantes 3900 , 14040-901 Ribeirão Preto , SP , Brazil.
J Chem Inf Model. 2020 Feb 24;60(2):522-536. doi: 10.1021/acs.jcim.9b00790. Epub 2019 Nov 26.
The ELBA coarse-grained force field was originally developed for lipids, and its water model is described as a single-site Lennard-Jones particle with electrostatics modeled by an embedded point-dipole, while other molecules in this force field have a three (or four)-to-one mapping scheme. Here, ELBA was applied to investigate the self-assembly processes of (DPC) micelle, 1,2-/ bicelles, and DPPC/cyclohexane/water reverse micelles through coarse-grained molecular dynamics (MD) simulations. New parameters were obtained using a simplex algorithm-based calibration procedure to determine the Lennard-Jones parameters for cyclohexane, dodecane, and cyclohexane-dodecane cross-interactions. Density, self-diffusion coefficient, surface tension, and mixture excess volume were found to be in fair agreement with experimental data. These new parameters were used in the simulations, and the obtained structures were analyzed for shape, size, volume, and surface area. Except for the shape of DPC micelles, all other properties match well with available experimental data and all-atom simulations. Remarkably, in agreement with experiments the rodlike shape of the DPPC reverse micelle is well described by ELBA, while all-atom data in the literature predicts a disclike shape. To further check the consistency of the force field in reproducing the correct shapes of reverse micelles, additional simulations were performed doubling the system size. Two distinct reverse micelles were obtained both presenting the rodlike shape and correct aggregation number.
ELBA 粗粒化力场最初是为脂质开发的,其水模型被描述为具有静电作用的单个 Lennard-Jones 粒子,静电作用由嵌入的点偶极子建模,而该力场中的其他分子具有三(或四)对一的映射方案。在这里,ELBA 被应用于通过粗粒化分子动力学 (MD) 模拟研究 (DPC) 胶束、1,2-/双胶束和 DPPC/环己烷/水反胶束的自组装过程。通过基于单纯形算法的校准程序获得了新的参数,以确定环己烷、十二烷和环己烷-十二烷交叉相互作用的 Lennard-Jones 参数。密度、自扩散系数、表面张力和混合超额体积与实验数据吻合良好。这些新参数用于模拟,并且对获得的结构进行了形状、大小、体积和表面积的分析。除了 DPC 胶束的形状外,所有其他性质都与可用的实验数据和全原子模拟很好地匹配。值得注意的是,与实验一致,ELBA 很好地描述了 DPPC 反胶束的棒状形状,而文献中的全原子数据预测出了盘状形状。为了进一步检查力场在复制反胶束正确形状方面的一致性,对系统尺寸进行了加倍的额外模拟。获得了两种不同的反胶束,均呈现棒状形状和正确的聚集数。