Department of Biophysics, Faculty of Biological Sciences, Tarbiat Modares University, Tehran, Iran.
Department of Mathematics, Tarbiat Modares University, Tehran, Iran.
PLoS One. 2019 Apr 5;14(4):e0214673. doi: 10.1371/journal.pone.0214673. eCollection 2019.
Lamellar and hexagonal lipid structures are of particular importance in the biological processes such as membrane fusion and budding. Atomistic simulations of formation of these phases and transitions between them are computationally prohibitive, hence development of coarse-grained models is an important part of the methodological development in this area. Here we apply systematic bottom-up coarse-graining to model different phase structures formed by 1,2-dioleoylphosphatidylethanolamine (DOPE) lipid molecules. We started from atomistic simulations of DOPE lipids in water carried out at two different water/lipid molar ratio corresponding to the lamellar Lα and inverted hexagonal HII structures at low and high lipid concentrations respectively. The atomistic trajectories were mapped to coarse-grained trajectories, in which each lipid was represented by 14 coarse-grained sites. Then the inverse Monte Carlo method was used to compute the effective coarse-grained potentials which for the coarse-grain model reproduce the same structural properties as the atomistic simulations. The potentials derived from the low concentration atomistic simulation were only able to form a bilayer structure, while both Lα and HII lipid phases were formed in simulations with potentials obtained at high concentration. The typical atomistic configurations of lipids at high concentration combine fragments of both lamellar and non-lamellar structures, that is reflected in the extracted coarse-grained potentials which become transferable and can form a wide range of structures including the inverted hexagonal, bilayer, tubule, vesicle and micellar structures.
层状和六方脂质结构在膜融合和出芽等生物过程中具有重要意义。这些相的形成和它们之间的转变的原子模拟在计算上是不可行的,因此粗粒化模型的发展是该领域方法学发展的重要组成部分。在这里,我们应用系统的自下而上的粗粒化方法来模拟由 1,2-二油酰基磷脂酰乙醇胺 (DOPE) 脂质分子形成的不同相结构。我们从在两种不同的水/脂质摩尔比下进行的 DOPE 脂质在水中的原子模拟开始,这两种摩尔比分别对应于低脂质浓度下的层状 Lα 和反相六方 HII 结构,以及高脂质浓度下的 Lα 和反相六方 HII 结构。原子轨迹被映射到粗粒轨迹上,其中每个脂质由 14 个粗粒位点表示。然后,逆蒙特卡罗方法被用来计算有效的粗粒化势,对于粗粒化模型,这些势可以再现与原子模拟相同的结构特性。从低浓度原子模拟得到的势只能形成双层结构,而在高浓度势下的模拟中则形成了 Lα 和 HII 脂质相。高浓度下脂质的典型原子构型结合了层状和非层状结构的片段,这反映在提取的粗粒化势中,这些势变得可转移,并可以形成包括反相六方、双层、管状、囊泡和胶束等多种结构。