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

Abstract 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.