Simulations of membranes and other interfacial systems using P2(1) and Pc periodic boundary conditions

We demonstrate the ease and utility of simulating heterogeneous interfacial systems with P2(1) and Pc periodic boundary conditions which allow, for example, lipids in a membrane to switch leaflets. In preliminary tests, P2(1) was shown to yield equivalent results to P1 in simulations of bulk water, a water/vacuum interface, and pure DPPC bilayers with an equal number of lipids per leaflet; equivalence of Pc and P1 was also demonstrated for the former two systems. P2(1) was further tested in simulations involving the spreading of an octane film on water, and equilibration of a DPPC bilayer from an initial condition containing different numbers of lipids in the two leaflets. Lastly, a simulation in P2(1) of a DOPC/melittin membrane showed significant passage of lipids to the melittin-containing leaflet from the initial distribution, and lends insight into the condensation of lipids by melittin.