| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Biophysical Journal 71: 1350-1355 (1996)
© 1996 the Biophysical Society
Biophysics Laboratory, Food and Drug Administration, Rockville, Maryland 20852-1448, USA. rwpastor@deimos.cber.nih.gov
ABSTRACT
As sketched in Fig. 1, a current molecular dynamics computer simulation of a lipid bilayer fails to capture significant features of the macroscopic system, including long wavelength undulations. Such fluctuations are intrinsically connected to the value of the macroscopic (or thermodynamic) surface tension (cf. Eqs. 1 and 9; for a related treatment, see Brochard et al., 1975, 1976). Consequently, the surface tension that might be evaluated in an MD simulation should not be expected to equal the surface tension obtained from macroscopic measurements. Put another way, the largest of the three simulations presented here contained over 16,000 atoms and required substantial computer time to complete, but modeled a system of only 36 lipids per side. From this perspective it is not surprising that the system is not at the thermodynamic limit. An important practical consequence of this effect is that simulations with fluctuating area should be carried out with a nonzero applied surface tension (gamma 0 of Fig. 2) even when the macroscopic tension is zero, or close to zero. Computer simulations at fixed surface area, which can explicitly determine pressure anisotropy at the molecular level, should ultimately lend insight into the value of gamma 0, including its dependence on lipid composition and other membrane components. As we have noted and will describe further in separate publications (Feller et al., 1996; Feller et al., manuscript in preparation), surface tensions obtained from simulations can be distorted by inadequate initial conditions and convergence, and are sensitive to potential energy functions, force truncation methods, and system size; it is not difficult, in fact, to tune terms in the potential energy function so as to yield surface tensions close to zero. This is why parameters should be tested extensively on simpler systems, for example, monolayers. The estimates of gamma 0 that we have presented here should be regarded as qualitative, and primarily underscore the assertion that the surface tension of a microscopically flat, simulation-sized patch is significantly greater than zero. As the simulation cell length increases, the surface tension that would be evaluated (or should be applied) decreases; in the limit of micrometer-sized simulation cells, gamma would approach zero or its appropriate thermodynamic value. The theories presented here also imply that the estimation of bilayer surface tension from monolayer data should take the degree of flatness into account. These conclusions are independent of the precise values of parameters such as bending constants. In conclusion, from the simulator's perspective, the question "What is the surface tension of a bilayer?" is better phrased as "What is the value of the applied surface tension necessary to simulate a particular experimental system with a given number of lipids?". As we have shown, the answer to the second question varies, but it should not be assumed a priori to equal zero.
This article has been cited by other articles:
 |
|
 |
|
  H. Jang, J. Zheng, and R. Nussinov Models of {beta}-Amyloid Ion Channels in the Membrane Suggest That Channel Formation in the Bilayer Is a Dynamic Process Biophys. J., September 15, 2007; 93(6): 1938 - 1949. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
J. Sonne, M. O. Jensen, F. Y. Hansen, L. Hemmingsen, and G. H. Peters Reparameterization of All-Atom Dipalmitoylphosphatidylcholine Lipid Parameters Enables Simulation of Fluid Bilayers at Zero Tension Biophys. J., June 15, 2007; 92(12): 4157 - 4167. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
J. B. Klauda, N. Kucerka, B. R. Brooks, R. W. Pastor, and J. F. Nagle Simulation-Based Methods for Interpreting X-Ray Data from Lipid Bilayers Biophys. J., April 15, 2006; 90(8): 2796 - 2807. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
A. Skibinsky, R. M. Venable, and R. W. Pastor A Molecular Dynamics Study of the Response of Lipid Bilayers and Monolayers to Trehalose Biophys. J., December 1, 2005; 89(6): 4111 - 4121. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
M. Venturoli, B. Smit, and M. M. Sperotto Simulation Studies of Protein-Induced Bilayer Deformations, and Lipid-Induced Protein Tilting, on a Mesoscopic Model for Lipid Bilayers with Embedded Proteins Biophys. J., March 1, 2005; 88(3): 1778 - 1798. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
R. W. Benz, F. Castro-Roman, D. J. Tobias, and S. H. White Experimental Validation of Molecular Dynamics Simulations of Lipid Bilayers: A New Approach Biophys. J., February 1, 2005; 88(2): 805 - 817. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
J. Wohlert and O. Edholm The Range and Shielding of Dipole-Dipole Interactions in Phospholipid Bilayers Biophys. J., October 1, 2004; 87(4): 2433 - 2445. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
D. Bemporad, C. Luttmann, and J. W. Essex Computer Simulation of Small Molecule Permeation across a Lipid Bilayer: Dependence on Bilayer Properties and Solute Volume, Size, and Cross-Sectional Area Biophys. J., July 1, 2004; 87(1): 1 - 13. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
F. Y. Jiang, Y. Bouret, and J. T. Kindt Molecular Dynamics Simulations of the Lipid Bilayer Edge Biophys. J., July 1, 2004; 87(1): 182 - 192. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
P. M. Kasson and V. S. Pande Molecular Dynamics Simulation of Lipid Reorientation at Bilayer Edges Biophys. J., June 1, 2004; 86(6): 3744 - 3749. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
T. Huber, A. V. Botelho, K. Beyer, and M. F. Brown Membrane Model for the G-Protein-Coupled Receptor Rhodopsin: Hydrophobic Interface and Dynamical Structure Biophys. J., April 1, 2004; 86(4): 2078 - 2100. [Abstract] [Full Text] [PDF]
|
|
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
