This Article Full Text Full Text (PDF) All Versions of this Article: biophysj.104.045716v1 87/5/3242    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by McWhirter, J. L. Articles by Voth, G. A. Search for Related Content PubMed PubMed Citation Articles by McWhirter, J. L. Articles by Voth, G. A.

Coupling Field Theory with Mesoscopic Dynamical Simulations of Multicomponent Lipid Bilayers

Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah

Correspondence: Address reprint requests to Professor Gregory A. Voth, Dept. of Chemistry, University of Utah, 315 S. 1400 E. Rm 2020, Salt Lake City, UT 84112-0850. Tel.: 801-581-7272; E-mail: voth{at}chem.utah.edu.

A method for simulating a two-component lipid bilayer membrane in the mesoscopic regime is presented. The membrane is modeled as an elastic network of bonded points; the spring constants of these bonds are parameterized by the microscopic bulk modulus estimated from earlier atomistic nonequilibrium molecular dynamics simulations for several bilayer mixtures of DMPC and cholesterol. The modulus depends on the composition of a point in the elastic membrane model. The dynamics of the composition field is governed by the Cahn-Hilliard equation where a free energy functional models the coupling between the composition and curvature fields. The strength of the bonds in the elastic network are then modulated noting local changes in the composition and using a fit to the nonequilibrium molecular dynamics simulation data. Estimates for the magnitude and sign of the coupling parameter in the free energy model are made treating the bending modulus as a function of composition. A procedure for assigning the remaining parameters in the free energy model is also outlined. It is found that the square of the mean curvature averaged over the entire simulation box is enhanced if the strength of the bonds in the elastic network are modulated in response to local changes in the composition field. We suggest that this simulation method could also be used to determine if phase coexistence affects the stress response of the membrane to uniform dilations in area. This response, measured in the mesoscopic regime, is already known to be conditioned or renormalized by thermal undulations.

This article has been cited by other articles:

				G. S. Ayton, P. D. Blood, and G. A. Voth Membrane Remodeling from N-BAR Domain Interactions: Insights from Multi-Scale Simulation Biophys. J., May 15, 2007; 92(10): 3595 - 3602. [Abstract] [Full Text] [PDF]

				Q. Shi and G. A. Voth Multi-Scale Modeling of Phase Separation in Mixed Lipid Bilayers Biophys. J., October 1, 2005; 89(4): 2385 - 2394. [Abstract] [Full Text] [PDF]

				G. S. Ayton, J. L. McWhirter, P. McMurtry, and G. A. Voth Coupling Field Theory with Continuum Mechanics: A Simulation of Domain Formation in Giant Unilamellar Vesicles Biophys. J., June 1, 2005; 88(6): 3855 - 3869. [Abstract] [Full Text] [PDF]