Mean field stochastic boundary molecular dynamics simulation of a phospholipid in a membrane.

Computer simulations of phospholipid membranes have been carried out by using a combined approach of molecular and stochastic dynamics and a mean field based on the Marcelja model. First, the single-chain mean field simulations of Pastor et al. [(1988) J. Chem. Phys. 89, 1112-1127] were extended to a complete dipalmitoylphosphatidylcholine molecule; a 102-ns Langevin dynamics simulation is presented and compared with experiment. Subsequently, a hexagonally packed seven-lipid array was simulated with Langevin dynamics and a mean field at the boundary and with molecular dynamics (and no mean field) in the center. This hybrid method, mean field stochastic boundary molecular dynamics, reduces bias introduced by the mean field and eliminates the need for periodic boundary conditions. As a result, simulations extending to tens of nanoseconds may be carried out by using a relatively small number of molecules to model the membrane environment. Preliminary results of a 20-ns simulation are reported here. A wide range of motions, including overall reorientation with a nanosecond decay time, is observed in both simulations, and good agreement with NMR, IR, and neutron diffraction data is found.

[1]  M. Karplus,et al.  Active site dynamics in protein molecules: A stochastic boundary molecular‐dynamics approach , 1985, Biopolymers.

[2]  Scott H. Northrup,et al.  Molecular mechanics of kink formation in lipid monolayers , 1984 .

[3]  Herman J. C. Berendsen,et al.  MOLECULAR-DYNAMICS SIMULATION OF A BILAYER-MEMBRANE , 1982 .

[4]  U. Singh,et al.  A NEW FORCE FIELD FOR MOLECULAR MECHANICAL SIMULATION OF NUCLEIC ACIDS AND PROTEINS , 1984 .

[5]  E. Helfand,et al.  Brownian dynamics study of transitions in a polymer chain of bistable oscillators , 1978 .

[6]  Michael L. Klein,et al.  Molecular dynamics study of a sodium octanoate micelle in aqueous solution , 1988 .

[7]  R. Mendelsohn,et al.  Quantitative determination of conformational disorder in the acyl chains of phospholipid bilayers by infrared spectroscopy. , 1989, Biochemistry.

[8]  J J Wendoloski,et al.  Molecular dynamics simulation of a phospholipid micelle. , 1989, Science.

[9]  J. Seelig,et al.  The dynamic structure of fatty acyl chains in a phospholipid bilayer measured by deuterium magnetic resonance. , 1974, Biochemistry.

[10]  J. Seelig,et al.  Deuterium order parameters in relation to thermodynamic properties of a phospholiped bilayer. A statistical mechanical interpretation. , 1975, Biochemistry.

[11]  P. Meier,et al.  Proton spin relaxation dispersion studies of phospholipid membranes , 1988 .

[12]  Robert Zwanzig,et al.  STATISTICAL ERROR DUE TO FINITE TIME AVERAGING IN COMPUTER EXPERIMENTS. , 1969 .

[13]  S. Harvey,et al.  Molecular mechanics model of supercoiled DNA. , 1989, Journal of molecular biology.

[14]  A. Saupe,et al.  Eine einfache molekular-statistische Theorie der nematischen kristallinflüssigen Phase. Teil II , 1960 .

[15]  H. Berendsen,et al.  Molecular dynamics simulation of a smectic liquid crystal with atomic detail , 1988 .

[16]  D. Engelman,et al.  Lipid bilayer thickness varies linearly with acyl chain length in fluid phosphatidylcholine vesicles. , 1983, Journal of molecular biology.

[17]  Scott H. Northrup,et al.  Molecular dynamics simulation of disorder transitions in lipid monolayers , 1985 .

[18]  Gerald D. Williams,et al.  New view of lipid bilayer dynamics from 2H and 13C NMR relaxation time measurements. , 1983, Proceedings of the National Academy of Sciences of the United States of America.

[19]  B. Brooks,et al.  An analysis of the accuracy of Langevin and molecular dynamics algorithms , 1988 .

[20]  D. Gruen.  A model for the chains in amphiphilic aggregates. 1. Comparison with a molecular dynamics simulation of a bilayer , 1985 .

[21]  S. Marčelja,et al.  Chain ordering in liquid crystals. II. Structure of bilayer membranes. , 1974, Biochimica et biophysica acta.

[22]  M. Karplus,et al.  Stochastic boundary conditions for molecular dynamics simulations of ST2 water , 1984 .

[23]  E Helfand.  Dynamics of Conformational Transitions in Polymers , 1984, Science.

[24]  Egbert Egberts.  Molecular dynamics simulation of multibilayer membranes. , 1988 .

[25]  A. Szabó,et al.  Model-free approach to the interpretation of nuclear magnetic resonance relaxation in macromolecules. 1. Theory and range of validity , 1982 .

[26]  L. Verlet.  Computer "Experiments" on Classical Fluids. I. Thermodynamical Properties of Lennard-Jones Molecules , 1967 .

[27]  G Büldt,et al.  Neutron diffraction studies on phosphatidylcholine model membranes. II. Chain conformation and segmental disorder. , 1979, Journal of molecular biology.

[28]  S. Marčelja.  Molecular Model for Phase Transition in Biological Membranes , 1973, Nature.

[29]  J. Mccammon,et al.  Molecular dynamics with stochastic boundary conditions , 1982 .

[30]  Martin Karplus,et al.  A simulation based model of NMR T1 relaxation in lipid bilayer vesicles , 1988 .

[31]  P. Hitchcock,et al.  A refinement analysis of the crystallography of the phospholipid, 1,2-dilauroyl-DL-phosphatidylethanolamine, and some remarks on lipid—lipid and lipid-protein interactions , 1977, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[32]  J. Seelig,et al.  Lipid conformation in model membranes and biological membranes , 1980, Quarterly Reviews of Biophysics.

[33]  Jean-Paul Ryckaert,et al.  Molecular dynamics of liquid n-butane near its boiling point , 1975 .

[34]  Martin Karplus,et al.  Brownian dynamics simulation of a lipid chain in a membrane bilayer , 1988 .

[35]  Herman J. C. Berendsen,et al.  MOLECULAR-DYNAMICS OF A BILAYER-MEMBRANE , 1983 .

[36]  M. Karplus,et al.  Inertial effects in butane stochastic dynamics , 1989 .

[37]  S. Rice,et al.  A molecular dynamics study of the structure of a model Langmuir monolayer of amphiphile molecules , 1988 .