A coarse-grained model for lipid bilayer formation, fusion, and its hydrodynamics*

In this paper a theoretical model for long-range, hydrophobic attraction between amphiphilic particles is developed to quantify the macroscopic assembly and mechanics of a lipid bilayer membrane in solvents. The non-local interactions between amphiphilic particles are obtained from the first domain variation of a hydrophobicity functional, giving rise to forces and torques (between particles) that dictate the motion of both particles and the surrounding solvent. The functional minimizer (that accounts for hydrophobicity at molecular-aqueous interfaces) is a solution to a boundary value problem of the screened Laplace equation. We reformulate the boundary value problem as a second-kind integral equation. Solving a mobility problem in Stokes flow is incorporated to obtain corresponding rigid body motion. The simulated fluid-particle systems exhibit a variety of multiscale behaviors over both time and length: Over short time scales, the numerical results show self-assembly for model lipid particles. For large system simulations, the particles form realistic configurations like micelles and bilayers. Over long time scales, the bilayer shapes emerging from the simulation appear to minimize a form of bending energy.