Phase separation behavior of mixed lipid systems at neutral and low pH: coarse-grained simulations with DMD/LIME

Development of a coarse-grained lipid model, LIME 2.0, for DSPE using multistate iterative Boltzmann inversion and discontinuous molecular dynamics simulations

We suggest an improved version of the intermediate resolution implicit solvent model for lipids, LIME, that was previously developed for use with discontinuous molecular dynamics (DMD) simulations. LIME gets its geometrical and the energy parameters between bonded and nonbonded pairs of coarse-grained (CG) sites from atomistic simulations. The improved model, LIME 2.0, uses multiple square wells rather than the single square well used in original LIME to obtain intermolecular interactions that more faithfully mimic those from atomistic simulations. The multi-state iterative Boltzmann inversion (MS-IBI) scheme is used to determine the interaction parameters. This means that a single set of interaction parameters between coarse-grained sites can be used to represent the lipid bilayers at different temperatures. The physical properties of CG DSPE lipid bilayer are calculated using CG simulations and compared to atomistic simulations results to verify the improved model. The phase transition temperature of the lipid bilayer is measured accurately and the lipid translocation phenomenon, " flip-flop" is observed through CG simulation. These results suggest that CG parameterization using multiple square-well and the MS-IBI scheme is well suited to the study of lipid bilayers cross a range of temperatures with DMD simulations.

Critical size dependence of domain formation observed in coarse-grained simulations of bilayers composed of ternary lipid mixtures

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.

Modeling nanoparticle wrapping or translocation in bilayer membranes

The spontaneous wrapping of nanoparticles by membranes is of increasing interest as nanoparticles become more prevalent in consumer products and hence more likely to enter the human body. We introduce a simulations-based tool that can be used to visualize the molecular level interaction between nanoparticles and bilayer membranes. By combining LIME, an intermediate resolution, implicit solvent model for phospholipids, with discontinuous molecular dynamics (DMD), we are able to simulate the wrapping or embedding of nanoparticles by 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) bilayer membranes. Simulations of hydrophilic nanoparticles with diameters from 10 Å to 250 Å show that hydrophilic nanoparticles with diameters greater than 20 Å become wrapped while the nanoparticle with a diameter of 10 Å does not. Instead this smaller particle became embedded in the bilayer surface where it can interact with the hydrophilic head groups of the lipid molecules. We also investigate the interaction between a DPPC bilayer and hydrophobic nanoparticles with diameters 10 Å to 40 Å. These nanoparticles do not undergo the wrapping process; instead they directly penetrate the membrane and embed themselves within the inner hydrophobic core of the bilayers.