TY - JOUR

T1 - Determining structural and mechanical properties from molecular dynamics simulations of lipid vesicles

AU - Braun, Anthony R.

AU - Sachs, Jonathan N.

PY - 2014/9/9

Y1 - 2014/9/9

N2 - We have developed an algorithm to determine membrane structure, area per lipid, and bending rigidity from molecular dynamics simulations of lipid vesicles. Current methods to extract structure from vesicle simulations define densities relative to the global center of mass of the vesicle. This approach ignores the long-wavelength fluctuations (undulations) that develop across the sphere and broaden the underlying structure. Our method establishes a local reference frame by defining a radially undulating reference surface (URS) and thereby removes the broadening effect of the undulations. Using an arc-length low-pass filter, we render the URS by defining the bilayer midplane on an equi-angular, ø-grid (colatitude, longitude). This surface is then expanded onto a truncated series of spherical harmonics. The spherical harmonic coefficients characterize the long-wavelength fluctuations that define both the local reference frame-used to determine the bilayers structure-and the area per lipid (AL) along the undulating surface. Additionally, the resulting power spectrum of spherical harmonic coefficients can be fit to a Helfrich continuum model for membrane bending in spherical geometry to extract bending rigidity (kc). kcvalues determined for both DMPC and DMPC + cholesterol (30 mol %) vesicles are consistent with values from corresponding flat-patch systems determined using an independent, previously published spectral method. These new tools to accurately extract structure, AL, and kcshould prove invaluable in evaluating the construction and equilibration of lipid vesicle simulations.

AB - We have developed an algorithm to determine membrane structure, area per lipid, and bending rigidity from molecular dynamics simulations of lipid vesicles. Current methods to extract structure from vesicle simulations define densities relative to the global center of mass of the vesicle. This approach ignores the long-wavelength fluctuations (undulations) that develop across the sphere and broaden the underlying structure. Our method establishes a local reference frame by defining a radially undulating reference surface (URS) and thereby removes the broadening effect of the undulations. Using an arc-length low-pass filter, we render the URS by defining the bilayer midplane on an equi-angular, ø-grid (colatitude, longitude). This surface is then expanded onto a truncated series of spherical harmonics. The spherical harmonic coefficients characterize the long-wavelength fluctuations that define both the local reference frame-used to determine the bilayers structure-and the area per lipid (AL) along the undulating surface. Additionally, the resulting power spectrum of spherical harmonic coefficients can be fit to a Helfrich continuum model for membrane bending in spherical geometry to extract bending rigidity (kc). kcvalues determined for both DMPC and DMPC + cholesterol (30 mol %) vesicles are consistent with values from corresponding flat-patch systems determined using an independent, previously published spectral method. These new tools to accurately extract structure, AL, and kcshould prove invaluable in evaluating the construction and equilibration of lipid vesicle simulations.

UR - http://www.scopus.com/inward/record.url?scp=84907221090&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84907221090&partnerID=8YFLogxK

U2 - 10.1021/ct500460u

DO - 10.1021/ct500460u

M3 - Article

AN - SCOPUS:84907221090

VL - 10

SP - 4160

EP - 4168

JO - Journal of Chemical Theory and Computation

JF - Journal of Chemical Theory and Computation

SN - 1549-9618

IS - 9

ER -