Search

The separation of lipid mixtures into thermodynamically stable phase-separated domains is dependent on lipid composition, temperature, and system size. Using molecular dynamics simulations, the line tension between thermodynamically stable lipid domains formed from ternary mixtures of di-C16:0 PC:di-C18:2 PC:cholesterol at 40:40:20 mol. % ratio was investigated via two theoretical approaches. The line tension was found to be 3.1 ± 0.2 pN by capillary wave theory and 4.7 ± 3.7 pN by pressure tensor anisotropy approaches for coarse-grained models based on the Martini force field. Using an all-atom model of the lipid membrane based on the CHARMM36 force field, the line tension was found to be 3.6 ± 0.9 pN using capillary wave theory and 1.8 ± 2.2 pN using pressure anisotropy approaches. The discrepancy between estimates of the line tension based on capillary wave theory and pressure tensor anisotropy methods is discussed. Inclusion of protein in Martini membrane lipid mixtures was found to reduce the line tension by 25%–35% as calculated by the capillary wave theory approach. To further understand and predict the behavior of proteins in phase-separated membranes, we have formulated an analytical Flory-Huggins model and parameterized it against the simulation results. Taken together these results suggest a general role for proteins in reducing the thermodynamic cost associated with domain formation in lipid mixtures and quantifies the thermodynamic driving force promoting the association of proteins to domain interfaces.

Reverse micelles (RMs) are recognized as a paradigm of molecular self-assembly and used in a variety of applications, such as chemical synthesis and molecular structure refinement. Nevertheless, many fundamental properties including their equilibrium size distribution, internal structure, and mechanism of self-assembly remain poorly understood. To provide an enhanced microscopic understanding of the assembly process and resulting structural distribution, we perform multiple nonequilibrium molecular dynamics simulations of dioctyl sulfosuccinate sodium salt (AOT) RM assembly, quantifying RM size, water core structure, and dynamics. Rapid assembly of smaller RM from a random mixture is observed to establish a constant AOT water loading within a nanosecond consistent with a diffusion-adsorption mechanism validated through the Monte-Carlo simulation of a model system. The structure of RM water cores and RM molecular volume during RM assembly is characterized during the AOT assembly process. A moment-closure equation is developed from a novel master equation model to elucidate the elementary events underlying the AOT self-assembly process. The resulting kinetic model is used to explore the role of monomer addition and dissociation, RM association and dissociation, and RM collision-induced exchange, all dependent on average RM size, which provides fundamental insight regarding the mechanisms and time scales for AOT RM self-assembly. The nascent dynamics that rapidly establish water loading, intermediate time scales of RM fusion, and longer time scale dynamics of inter-RM exchange essential in establishing the equilibrium condition are quantified through these kinetic models. Overall, this work provides insight into AOT RM self-assembly and provides a general theoretical framework for the analysis of the molecular self-assembly dynamics and mechanism.

Model cellular membranes are known to form micro- and macroscale lipid domains dependent on molecular composition. The formation of macroscopic lipid domains by lipid mixtures has been the subject of many simulation investigations. We present a critical study of system size impact on lipid domain phase separation into liquid-ordered and liquid-disordered macroscale domains in ternary lipid mixtures. In the popular di-C16:0 PC:di-C18:2 PC:cholesterol at 35:35:30 ratio mixture, we find systems with a minimum of 1480 lipids to be necessary for the formation of macroscopic phase separated domains and systems of 10 000 lipids to achieve structurally converged conformations similar to the thermodynamic limit. To understand these results and predict the behavior of any mixture forming two phases, we develop and investigate an analytical Flory-Huggins model which is recursively validated using simulation and experimental data. We find that micro- and macroscale domains can coexist in ternary mixtures. Additionally, we analyze the distributions of specific lipid-lipid interactions in each phase, characterizing domain structures proposed based on past experimental studies. These findings offer guidance in selecting appropriate system sizes for the study of phase separations and provide new insights into the nature of domain structure for a popular ternary lipid mixture.