Although computer simulation of biological molecules has seen widespread growth and is widely accepted as an important biochemical tool, it is hampered by lim ited computing resources. Biomolecular systems, by necessity, contain a large number of interaction sites. In many cases, these sites interact over quite large dis tances. Further, the time scales of biological interest are long, which requires that simulations of dynamical properties at the atomic level must be lengthy to ade quately probe these motions. We address these issues through discussions of atomic-level molecular dynam ics simulations of biological lipid bilayer membranes, which are key constructs in biochemistry. These simu lations reproduce many experimental observables and provide a degree of resolution currently unavailable experimentally. The lengths of these simulations, the longest of which was 2 nanoseconds, were sufficient to effectively sample many of the motions governing the behavior of biomembranes. Examples are given showing the importance of long-range interactions. The number of interaction sites required by these sim ulations is discussed, particularly the need for explicit representation of solvent molecules.
Get full access to this article
View all access options for this article.
References
1.
Aloisi, G., Foresti, M.L., and Guidelli, R.1990. A Monte Carlo simulation of water molecules near a charged wallJ. Chem. Phys.91:5592-5596.
2.
Alper, H.E., Bassolino, D., and Stouch, T.R.1993a. Computer simulation of a phospholipid monolayer-water system: The influence of long range forces on water structure and dynamicsJ. Chem. Phys.98(12):9798-9807.
3.
Alper, H.E., Bassolino, Donna, and Stouch, T.R.1993b. The limiting behavior of water hydrating a phospholipid monolayer : A computer simulation study. J. Chem. Phys.99(7):5547-5559.
4.
Alper, H.E., and Levy, R.M.1989. Computer simulations of the dielectric properties of water: Studies of the simple point charge and transferrable intermolecular potential modelsJ. Chem. Phys.91(2):1242.
5.
Alper, H.E., and Levy, R.M.1993. Dielectric and thermodynamics response of a generalized reaction field model for liquid state simulations. J. Chem. Phys.99(12):9847-9852.
6.
Applegate, Kenneth R., and Glomset, John A.1986. Computer-based modeling of the conformation and packing properties of docosahexaenoic acidJ. Lipid Res.27:658-680.
7.
Bareman, James P., Cardini, Gianni, and Klein, Michael L.1988. Characterization of structural and dynamical behavior in monolayer of long-chain molecules using molecular-dynamics calculationsPhys. Rev. Let.60(21):2152-2155.
8.
Bassolino-Klimas, Donna, Alper, H.E., and Stouch, T.R.1993. Solute diffusion in lipid bilayer membranes : an atomic level study by molecular dynamics simulationBio chemistry32:12624-12637.
9.
Belhadj, M., Alper, H.E., and Levy, R.M.1991. Molecular dynamics simulations of water with Ewald summation for the long range electrostatic interactionsChem. Phys. Lett.1:13.
10.
Berendsen, H.J.C., Postma, J.P.M., van Gunsteren, W.F., Di Nola, A., and Haah, J.R.1984. Molecular dynamics with coupling to an external bathJ. Chem. Phys.81(8):3684.
11.
Berendsen, H.J.C., Postma, J.P.M., van Gunsteren, W.F., and Hermans, J.1981. Interaction models for water in relation to protein hydration In Intermolecular Forces, edited by B. Pullman. Dordrecht: Reidel , pp. 331-342.
12.
Berkowitz, M.L., and Raghavan, K.1991. Computer simulation of a water/membrane interfaceLangmuir7:1442-1044.
13.
Binder, H., and Peinel, G.1985. Behaviour of water at membrane surfaces: A molecular dynamics study J. Molec. Struc. (Theochem), 123:155-163.
14.
Brasseur, R., DeLoof, H., Ruysschaert, J.M., and Rosseneu, M.1988. Conformational analysis of lipid-associating proteins in a lipid environmentBiochemica et Biophysica Acta943:95-102.
15.
Brasseur, R., Goormaghtigh, Erik, and Ruysschaert, Jean-Marie.1981. Theoretical conformational analysis of phospholipids bilayersBiochem. Biophys. Res. Commun.103(1):301-310.
16.
Brooks III, C.1987. The influence of long range force truncation on the thermodynamics of aqueous ionic solutionsJ. Chem. Phys.86(9):5156.
17.
Brooks III, C.L., Karplus, M., and Pettitt, B.M.1988. Prateins: A theoretical perspectiveNew York: John Wiley and Sons.
18.
Buldt, Georg, and Wohlgemuth, Roland.1981. The headgroup conformation of phospholipids in membranesJ. Membrane Biol.58:81-100.
19.
Cox, S.R., and Williams, D.E.1981. Representation of the molecular electrostatic potential by a net atomic charge modelJ. Comput. Chem.2(3):304-323.
20.
Cramer, Diego, Ben-Shaul, Avinoam, Chen, Zhong-Ying, and Gelbart, William M.1992. Monte Carlo and mean-field studies of successive phase transitions in rod monolayersJ. Chem. Phys.96(3):2236-2252.
21.
Damodaran, K.V., Merz, Kenneth M., and Gaber, Bruce Paul.1992. Structure and dynamics of the dilauroylphosphatidylethanolamine lipid bilayerBiochemistry31(33):7656-7664.
22.
Davis, Malcolm E., and McCammon, J. Andrew.1990. Electrostatics in biomolecular structure and dynamicsChemical Reviews90:509-521.
23.
De Loof, Hans, Harvey, Stephen C., Segrest, Jere P., and Pastor, Richard W.1991. Mean field stochastic boundary molecular dynamics simulation of a phospholipid in a membraneBiochemistiy30: 2099-2113.
24.
Dill, K.A., and Flory, P.J.1981. Molecular organization in micelles and vesiclesProc. Natl. Acad. Sci. USA78:676-680.
25.
Dill, Ken A., Koppel, Daniel E., Cantor, Robert S., Dill, James D., Bendedouch, Dalila, and Chen, Sow-Hsin.1984. Molecular conformations in surfactant micellesNature309:42-45.
26.
Din, Hon-Qiang, Karasawa, Naoki, and Goddard, William A., III.1992. Atomic level simulations on a million particles : The cell multipole method for Coulomb and London nonbond interactionsJ. Chem. Phys.97(6): 4309-4315.
27.
Drouffe, J.-M., Maggs, A.C., and Leibler, S.1991. Computer simulations of self-assembled membranesScience254: 1353-1356.
28.
Dufourc, Erick J., Mayer, Christian, Stohrer, Jurgen, Althoff, Gerhard, and Kothe, Gerd.1992. Dynamics of phosphate head groups in biomembranesBiophys. J.61:42-57.
29.
Edholm, O., Berendsen, H.J.C., and van der Ploeg, P.1983. Conformational entropy of a bilayer membrane derived from a molecular dynamics simulationMolecular Physics48(2):379-388.
30.
Edholm, O. and Johansson, J.1987. Lipid bilayer polypeptide interactions studied by molecular dynamics simulationEur. Biophys. J.14:203-209.
Egberts, E., and Berendsen, H.J.C.1988. Molecular dynamics simulation of a smectic liquid crystal with atomic detailJ. Chem. Phys.89:3718-3732.
33.
Engelmann, A.R., Llanos, C.M., Nyholm, P.G., Tapia, O., and Pascher, I.1987. Studies on membrane hydrationJ. Mol. Struc.151:81-102.
34.
Finean, J.B., Coleman, R., and Michell, R.H.1978. Membranes and their cellular functions2nd ed., Oxford: Blackwell.
35.
Franks, N.P.1976. Structural analysis of hydrated egg lecithin and cholesterol bilayersJ. Mol. Biol.100:345-358.
36.
French, Alfred E., and Brady, John W.1990. Computer modeling of carbohydrate moleculesWashington, D.C.: American Chemical Society .
37.
Gawrisch, K., Ruston, D., Zimmerberg, J., Parsegian, V.A., Rand, P.R., and Fuller, N.1992. Membrane dipole poten
38.
tials, hydration forces, and the ordering of water at membrane surfacesBiophys. J.61:1213.
39.
Gruen, D.W.R.1980. A statistical mechanical model of the lipid bilayer above its phase transitionBiochem. Biophys. Acta595: 161-183.
40.
Gunn, John R., and Dawson, K.A.1989. Microscopic model of amphiphilic assemblyJ. Chem. Phys.91(10):6393-6403.
41.
Hagler, A.T., Lifson, S., and Dauber, P.1979. Consistent force field studies of intermolecular forces in hydrogen bonded crystals. 2. A benchmark for the comparison of alternative force fieldsJ. Amer. Chem. Soc.101(18):5122-5230.
42.
Horwitz, A.F., and Klein, M.P.1973. The fluid state of lecithin bilayersJ. Supramolecular Structure (New York: Alan R. Liss, Inc.), 281-284.
43.
Hussin, A., and Scott, H.L.1987. Density and bonding profiles of inter-bilayer water as functions of bilayer separation: a Monte Carlo studyBio chimica et Biophysica Acta897:423-430.
44.
Israelachvili, Jacob N., and Wennerstrom, Hakan.1992. Entropic forces between amphiphilic surfaces in liquidsJ. Phys. Chem.96: 520-531.
45.
Jacobs, R.E., and White, S.H.1984. Behavior of hexane dissolved in DMPC bilayers: An NMR and calorimetric studyJ. Amer. Chem. Soc.106:915-920.
46.
Jonsson, B., Edholin, O., and Teleman, O.1986. Molecular dynamics simulations of a sodium octanoate micelle in aqueous solutionJ. Chem. Phys.85:2263-2271.
47.
Karaborni, S., and Toxvaerd, S.1992. Molecular dynamics simulations of Langmuir monolayers : A study of structure and thermodynamicsJ. Chem. Phys.96(7):5505-5515.
48.
Khalatur, Pavel G., and Pavlov, Alexander S.1987. Molecular motions in a liquid-crystalline lipid bilayer. Molecular dynamics simulationMakromol. Khim.188:3029-3040.
49.
King, Glen I., Jacobs, Russell E., and White, Stephen H.1985. Hexane dissolved in dioleoyllecithin bilayers has a partial molar volume of approximately zeroBiochemistry24:4637-4645.
50.
Klose, G. and Gawrisch, K.1981. Lipid-water interaction in model-membranesStudia Biophysica84(1):21-22.
51.
Knott, R.B., and Schoeborn, P.1986. Quantitation of water in membranes by neutron diffraction and X-ray techniquesMethods in Enzymology127 (Academic Press, Inc.).
52.
Kodama, Michiko.1986. Phase transition phenomena induced by the successive appearances of new type of aggregation states of water molecules in the Lα-DPPC-water systemThermochimica Acta109:81-89.
53.
Kodama, M., Kuwabara, M., and Seki, S. 1982.
54.
Successive phase-transition phenomena and phase diagram of the phosphatidylcholine-water system as revealed by differential scanning calorimetryBiochim. Biophys. Acta689:567-570.
55.
Kroll, D.M., and Gompper, G.1992. The conformation of fluid membranes: Monte Carlo simulationsScience255: 968-971.
56.
Larson, R.G.1992. Monte Carlo simulation of microstructural transitions in surfactant systemsJ. Chem. Phys.96: 7904-7918.
57.
Lewis, B.A., and Engelman, D.M.1983. Lipid bilayer thickness varies linearly with acyl chain length in fluid phosphatidylcholine vesiclesJ. Mol. Biol.166:211-217.
58.
Liang, Congxin, Ewig, Carl S., Stouch, Terry R., and Hagler, Arnold T.1993. ab initio studies of lipid model species. 1. Dimethylphosphate and methylpropylphosphate anionsJ. Amer. Chem. Soc.115(4):1537-1545.
59.
Liang, Congxin, Ewig, Carl S., Stouch, Terry R., and Hagler, Arnold T.1994. ab initio studies of lipid model species. 2. Inositol (In review).
60.
Loncharich, R.J., and Brooks, B.R.1989. The effects of truncating long-range forces on protein dynamicsProteins: Structure, Function, and Genetics.36:32.
61.
Marcelja, S.1974. Chain ordering in liquid crystals. II Structure of bilayer
62.
membranesBiochem. Biophys. Acta367:165.
63.
Marsh, D.1974. Statistical mechanics of the fluidity of phospholipid bilayers and membranesJ. Membrane Biol.18:145-162.
64.
McCammon, J.A., and Harvey, S.C.1987. Dynamics of proteins and nucleic acidsNew York: Cambridge University Press .
65.
Meraldi, Jean-Paul, and Schlitter, Jurgen.1981. A statistical mechanical treatment of fatty acyl chain order in phospholipid bilayers and correlation with experimental dataBiochimica et Biophysica Acta645:183-192.
66.
Milik, Mariusz, and Skolnick, Jeffrey.1993. Insertion of peptide chains into lipid membranes : An off-lattice Monte Carlo dynamics modelProteins: Structure, Function, and Genetics15: 10-25.
67.
Mit, B., Hubers, P.A.J., Esselink, K., Rupert, L.A.M., van Os, M.M., and Schlijper, A.G.1990. Computer simulations of a water/oil interface in the presence of micellesNa ture348:624-625.
68.
Mouritsen, Ole G.1990. Computer simulation of cooperative phenomena in lipid membranes In Molecular description of biological membranes by computer aided conformational analysis, Volume I, edited by Robert Brasseur. CRC Press, Boston.
69.
Muller-Plathe, F., Rogers, S.C., and van Gunsteren, W.F.1992. Computational evidence for anomalous diffusion of small molecules in amorphous polymersChem. Phys. Letts.199:237-243.
70.
Mulliken, R.S.1955. Electronic population analysis on LCAO-MO molecular wave functions. I. J. Chem. Phys.23:1833-1840.
71.
Murthy, N.S., and Worthington, C.R.1991. X-ray diffraction evidence for the presence of discrete water layers on the surface of membranesBiochimica et Biophysica Acta1062:172-176.
72.
Nagle, J.F.1976. Theory of lipid monolayer and bilayer phase transitions: Effect of headgroup interactionsJ. Membrane Biol.27:233-250.
73.
Nelander, J.C., and Blaurock, A.E.1978. Disorder in nerve myelin: Phasing the higher order reflections by means of the diffuse scatter. J. Mol. Biol.1118:497-532.
74.
Nicklas, K., Bocker, J., Schlenkrich, M., Brickmann, J., and Bopp, P.1991. Molecular dynamics studies of the interface between a model membrane and an aqueous solutionBiophys. J.60:261-272.
75.
Nilsson, Lennart, and Karplus, Martin.1986. Empirical energy functions for energy minimization and dynamics of nucleic acidsJ. Computational Chemistry5:591-616.
76.
Northrup, S.H.1984. Molecular mechanics of kink formation in lipid monolayersJ. Phys. Chem.88:3441-3446.
77.
Northrup, S.H., and Curvin, M.S.1985. Molecular dynamics simulation of disorder transitions in lipid monolayersJ. Phys. Chem.89:4707-4713.
78.
O'Leary, Timothy J.1987. Lateral diffusion of lipids in complex biological membranesProc. Natl. Acad. Sci. USA84:429-433.
79.
Parsegian, V.A., Rand, R.P., and Rau, D.C.1985. Hydration forces: What next?Chemica Scripta25:28-31.
80.
Pastor, R.W.1990. Determination of chain conformations in the membrane interior by Brownian dynamics simulations In Molecular description of biological membranes by computer aided conformational analysis , Volume I, edited by Robert Brasseur. Boston : CRC Press .
81.
Pastor, Richard W., Venable, Richard M. , and Karplus, Martin.1991. Model for the structure of the lipid bilayerProc. Natl. Acad. Sci. USA88: 892-896.
82.
Pastor, R.W., Venable, R.M., Karplus, M., and Szabo, A.1988. A simulation based model of NMR T1 relaxation in lipid bilayer vesiclesJ. Chem. Phys.89:1128.
83.
Pink, David A.1984. Theoretical studies of phospholipid bilayers and monolayers. Perturbing probes, monolayer phase transitions, and computer simulations of lipid-protein bilayersCan. J. Biochem. Cell Biol.62:760-777.
84.
Pink, D.A., and Hamboyan, H.1990. Effect of integral proteins upon bilayer permeability to small ionsEur. Biophys. J.18:245-253.
85.
van der Ploeg, P., and Berendsen, H.J.C.1982. Molecular dynamics simulation of a bilayer membraneJ. Chem. Phys.76: 3271-3276.
86.
van der Ploeg, P., and Berendsen, H.J.C.1983. Molecular dynamics of a bilayer membraneMol. Phys.49:233-248.
87.
Raghavan, K., Reddy, M., and Berkowitz, Max L.1992. A molecular dynamics study of the structure and dynamics of water between dilauroylphosphatidylethanolamine bilayersLangmuir8:233— 240.
88.
Rand, R.P., and Parsegian, V.A.1989. Hydration forces between phospholipid bilayersBiochemica et Biophysica Acta988:351-376.
89.
Ryckaert, Jean-Paul, McDonald, Ian R., and Klein, Michael L.1989. Disorder in the pseudohexagonal rotator phase of n-alkanes: Molecular-dynamics calculations for tircosaneMol. Phys.67(5):957-979.
90.
Saxton, Michael J.1992. Lateral diffusion and aggregationBiophys. J.61: 119-128.
91.
Scheringer, M., Hilfer, R., and Binder, K.1992. Orientational ordering in lipid monolayers: A two-dimensional model of rigid rods grafted onto a latticeJ. Chem. Phys.96(3):2269-2277.
92.
Schreiber, H., and Steinhauser, O.1992a. Taming cut-off induced artifacts in molecular dynamics studies of solvated polypeptides. The reaction field methodJ. Mol. Biol.228:909.
93.
Schreiber, H., and Steinhauser, O.1992b. Molecular dynamics simulations of solvated polypeptides : why the cut-off scheme does not workChem. Phys.168:75.
94.
Scott, H.L.1986. Monte Carlo calculations of order parameter profiles in models of lipid-protein interactions in bilayersBiochemistry25:6122-6126.
95.
Scott, H.L., and Kalaskar, S.1989. Lipid chains and cholesterol in model membranes: A Monte Carlo studyBiochemistry28:3687-3691.
96.
Scott, H.L., and Lee, C.Y.1980. The surface tension of lipid water interfaces : Monte Carlo simulationsJ. Chem. Phys.73(10):5351-5333.
97.
Seelig, J., and Seelig, A.1980. Lipid conformation in model membranes and biological membranesQ. Rev. Biophys. 13:19-61.
98.
Seelig, Joachim, Gally, Hans-Ulrich, and Wohlgemuth, Roland.1977. Orientation and flexibilitv of the choline head group in phosphatidylcholine bilayersBiochim. et Biophys. Acta467: 109-119.
99.
Simon, S.A., and McIntosh, T.J.1989. Magnitude of solvation pressure depends on dipole potentialProc. Natl. Acad. Sci. USA86:9263.
100.
Small, D.M.1967. Phase equilibria and structure of dry and hydrated egg lecithinJ. Lipid Researeh8: 551-557.
101.
Smith, P.E., and Pettitt, B.M.1991. Peptides in ionic solutions: A comparison of the Ewald and switching function techniquesJ. Chem. Phys.95(11):8430.
102.
Steinhauser, O.1981. Reaction field simulations of waterMol. Phys.45(2): 335.
103.
Stillinger, F.H., and Rahman, A.1974. Improved simulations of liquid water by molecular dynamicsJ. Chem. Phys.60(4): 1545.
104.
Stouch, T.R.1993. Lipid membrane structure and dynamics studied by all-atom molecular dynamics simulations of hydrated phospholipid bilayersMolecular simulations10(2-3):317-345.
105.
Stouch, T.R., and Williams, Donald E.1992. Conformational dependence of electrostatic potential derived charges of a lipid headgroup: GlycerylphosphorylcholineJ. Comput. Chem.13(5):622.
106.
Stouch, Terry R., and Williams, Donald E.1993. On the conformational dependence of electrostatic potential derived charges: Studies of the fitting procedureJ. Comput. Chem.14(7):858-866.
107.
Stouch, T.R., Ward, K.B., Altieri, A., and Hagler, A.T.1991. Simulations of lipid crystals: Characterization of potential energy functions and parameters for lecithin moleculesJ. Comput. Chem.12(8):1033-1046.
108.
Swindoll, R.D., and Haile, J.M.1984. A multiple time-step method for molecular dynamics simulations of fluids of chain moleculesJ. Comp. Phys.53:289.
109.
Tasaki, K., McDonald, S., and Brady, J.W.1993. Observations concerning the treatment of long-range interactions in molecular dynamics simulationsJ. Comp. Chem.14(3):278-284.
110.
Teleman, Olle, and Jonsson, Bo.1986. Vectorizing a general purpose M.D. simulation programJ. Comp. Chem.7:58.
111.
Tuckerman, Mark E., Martyna, Glenn J., and Berne, Bruce J.1990. Molecular dynamics algorithm for condensed systems with multiple time stepsJ. Chem. Phys.93(2): 1287-1291.
112.
Turner, G.L., and Oldfield, E.1979. Effect of local anaesthetics on hydrocarbon chain order in membranesNature277:669-670.
113.
Ulrich, A., Volke, F., and Watts, A.1990. The dependence of phospholipid head-group mobility on hydration as studied by deuterium-NMR spin-lattice relaxation time measurementsChem. Phys. Lipids55:61.
114.
Verlet, L.1967. Computer experiments on classical fluids. I. Thermodynamical properties of Lennard-Jones moleculesPhys. Rev.159(1): 98-103.
115.
Wallqvist, A.1991. Molecular dynamics study of a hydrophobic aggregate in an aqueous solution of methaneJ. Phys. Chem.95:8921-8927.
116.
Watanabe, K., Ferrario, M., and Klein, M.1988. Molecular dynamics study of a sodium octanoate micelle in aqueous solutionJ. Phys. Chem.92: 818-821.
117.
Watts, R.O.1974. Monte-Carlo studies of liquid waterMol. Phys.28(4): 1069.
118.
Weiner, M.C., King, G.I., and White, S.H.1991. Structure of a fluid dioleoylphosphatidylcholine bilayer determined by joint refinement of x-ray and neutron diffraction dataBiophys. J.60:568-576.
119.
Weiner, Scott J., Kollman, Peter A., Nguyen, Dzung T., and Case, David A.1986. An all atom force field for simulations of proteins and nucleic acidsJ. Comp. Chem.7(2):230-252.
120.
Wendoloski, I.J., Kimatian, S.J., Schutt, C.E., and Salemme, F.R.1989. Molecular dynamics simulation of a phospholipid micelleScience243: 636-638.
Williams, D.E., and Stouch, T.R.1993. Characterization of force fields for lipid molecules: Applications to crystal structuresJ. Comp. Chem.14(8):1066-1076.
123.
Williams, D.E., and Yan, Ji-Min.1988. Point-charge models for molecules derived from least-squares fitting of the electric potentialAdvances in Atomic and Molecular Physics23: 87-130.
124.
Woods, M.C., Haile, J.M., and O'Connell, J.P.1986. Internal structure of a model micelle via computer simulation. 2. Spherically confined aggregates with mobile head groupsJ. Phys. Chem.90: 1875-1885.
125.
Wong, C.F., and McCammon, J.A.1986. Isr. J. Chem.27:211.
126.
Yoon, D.Y., Smith, G.D., and Matsuda, T.1993. A comparison of a united atom and an explicit atom model in simulations of polymethylene.J. Chem. Phys.98(2):10037-10043.
