Nonequilibrium Molecular Dynamics Method to Generate Poiseuille-Like Flow between Lipid Bilayers

There are various flows inside and outside cells in vivo. Nonequilibrium molecular dynamics (NEMD) simulation is a useful tool for understanding the effects of these flows on the dynamics of biomolecules. We propose an NEMD method to generate a Poiseuille-like flow between lipid bilayers. We extended the conventional equilibrium MD method to produce a flow by adding constant external force terms to the water molecules. Using the Lagrange multiplier method, the center of mass of the lipid bilayer is constrained so that the flow does not sweep away the lipid bilayer, but the individual lipid molecules fluctuate. The temperature of the system is controlled properly in the solution and membrane by using the Nosé-Hoover thermostat. We found that the flow velocity increases linearly as the applied external force term increases. It is possible to estimate the appropriate value of acceleration to generate a flow with an arbitrary velocity using this proportional relation once a single short MD simulation is performed. We also found that the flow between two lipid bilayers is slower than the analytical solution of the Navier-Stokes equations between rigid parallel plates due to the interactions between water molecules and the membrane. This method can be applied not only to a flow on lipid membranes but also to a flow on soft surfaces generally.