Molecular dynamics simulation of a phospholipid membrane
Tóm tắt
We present the results of molecular dynamics (MD) simulations of a phospholipid membrane in water, including full atomic detail. The goal of the simulations was twofold: first we wanted to set up a simulation system which is able to reproduce experimental results and can serve as a model membrane in future simulations. This goal being reached it is then further possible to gain insight in to those properties that are experimentally more difficult to access. The system studied is dipalmitoylphosphatidylcholine/water, consisting of 5408 atoms. Using original force field parameters the membrane turned out to approach a gel-like state. With slight changes of the parameters, the system adopted a liquid-crystalline state. Separate 80 ps runs were performed on both the gel and liquid-crystalline systems. Comparison of MD results with reliable experimental data (bilayer repeat distance, surface area per lipid, tail order parameters, atom distributions) showed that our simulations, especially the one in the liquid-crystalline phase, can serve as a realistic model for a phospholipid membrane. Further analysis of the trajectories revealed valuable information on various properties. In the liquid-crystalline phase, the interface turns out to be quite diffuse, with water molecules penetrating into the bilayer to the position of the carbonyl groups. The 10–90% width of the interface turns out to be 1.3 nm and the width of the hydrocarbon interior 3.0 nm. The headgroup dipoles are oriented at a small angle with respect to the bilayer plane. The resulting charge distribution is almost completely cancelled by the water molecules. The electron density distribution shows a large dip in the middle of the membrane. In this part the tails are more flexible. The mean life time between dihedral transitions is 20 ps. The average number of gauche angles per tail is 3.5. The occurrence of kinks is not a significant feature.
Tài liệu tham khảo
Akutsu H (1981) Direct determination by raman scattering of the conformation of the choline group in phospholipid bilayers. Biochemistry 20:7359–7366
Akutsu H, Kyogoku Y (1977) Conformational difference in the polar groups of phosphatidylcholine and phosphatidylethanolamine in aqueous phase. Chem Phys Lipids 18:285–303
Bareman JP, Klein ML (1990) Collective tilt behavior in dense, substrate-supported monolayers of long-chain molecules: a molecular dynamics study. J Phys Chem 94:5202–5205
Berendsen HJC, Postma JPM, Gunsteren WF van, Hermans J (1981) Interaction models for water in relation to protein hydration. In: Pullman B (ed) Intermolecular forces. Reidel, Dordrecht, pp 331–342
Berendsen HJC, Postma JPM, Gunsteren WF van, Dinola A, Haak JR (1984) Molecular dynamics with coupling to an external bath. J Chem Phys 81:3684–3690
Braganza LF, Worcester DL (1986) Hydrostatic pressure induces hydrocarbon chain interdigitation in single-component phospholipid bilayers. Biochemistry 25:2591–2596
Brown MF, Seelig J, Häberlen U (1979) Structural dynamics in phospholipid bilayers from deuterium spin-lattice relaxation time measurements. J Chem Phys 70:5045–5053
Büldt G, Gally HU, Seelig A, Seelig J, Zaccai G (1978) Neutron diffraction studies on selectively deuterated phospholipid bilayers. Nature 271:182–184
Büldt G, Gally HU, Seelig J (1979) Neutron diffraction studies on phosphatidylcholine model membranes. J Mol Biol 134:673–691
Cevc G (1987) How membrane chain melting properties are regulated by the polar surface of the lipid bilayer. Biochemistry 26:6305–6310
Chandrasekhar I(1992) In: Gaber BP, Easwaran KRK (eds) Biomembrane structure and function. Adenine Press, New York, p 353
Chapman D, Williams RW, Ladbrooke BD (1967) Physical studies of phospholipids. VI. Thermotropic and lyotropic mesomorphism of some 1,2-diacylphosphatidylcholines (lecithins). Chem Phys Lipids 1:445–475
Egberts E (1988) Molecular dynamics simulation of multibilayer membranes. Ph.D.Thesis, Rijksuniversiteit Groningen
Egberts E, Berendsen HJC (1988) Molecular dynamics simulation of a smectic liquid crystal with atomic detail. J Chem Phys 89:3718–3732
Fahey PF, Webb WW (1978) Lateral diffusion in phospholipid bilayer membranes and multilamellar liquid crystals. Biochemistry 17:3046–3053
Frischleder H (1980) Quantum-chemical and empirical calculations on phospholipids. VII The conformational behavior of the phosphatidylcholine headgroup in layer systems. Chem Phys Lipids 27:83–92
Frischleder H, Peinel G (1982) Quantum-chemical and statistical calculations on phospholipids. Chem Phys Lipids 30:121–158
Gally HU, Niederberger W, Seelig J (1975) Conformation and motion of the choline head group in bilayers of dipalmitoyl-e-snphosphatidylcholine. Biochemistry 14:3647–3652
Gruen DWR (1981) A mean-field model of the alkane-saturated lipid bilayer above the phase transition. I. Development of the model. Biophys J 33:149–166
Gunsteren WF van, Berendsen HJC (1977) Algorithms for macromolecular dynamics and constraint dynamics. Mol Phys 34:1311–1327
Gunsteren WF van, Berendsen HJC (1987) Groningen molecular simulation (GROMOS) library manual, Biomos, Nijenborgh 4, 9747 AG Groningen, The Netherlands
Hauser H, Pascher I, Pearson RH, Sundell S (1981) Preferred conformation and molecular packing of phosphatidylethanolamine and phosphatidylcholine. Biochim Biophys Acta 650:21–51
Hockney RW Eastwood JW (1981) Computer simulation using particles. McGraw-Hill, New York
Horvath LI, Cirak J, Vigh L (1980) Relation of raman order parameters to spin labeling parameters. Chem Phys Lipids 27:237–250
Inoko Y, Mitsui T (1978) Structural parameters of dipalmitoyl phosphatidylcholine lamellar phases and bilayer phase transitions. J Phys Soc Japan 44:1918–1924
Janiak MJ, Small DM, Shipley GG (1976) Nature of the thermal pretransiton of synthetic phospholipids: dimyristoyl- and dipalmitoyllecithin. Biochemistry 15:4575–4580
Janiak MJ, Small DM, Shipley GG (1979) Temperature and compositional dependence of the structure of hydrated dimyristoyl lecithin. J Biol Chem 254:6068–6078
Jönsson B, Edholm O, Teleman O (1986) Molecular dynamics simulations of a sodium octanoate micelle in aqueous solution. J Chem Phys 85:2259–2271
Levine YK, Wilkins MHF (1971) Structure of oriented lipid bilayers. Nature New Biol 230:69–72
Lewis BA, Engelman DM (1983) Lipid bilayer thickness varies linearly with acyl chain length in fluid phosphatidylcholine vesicles. J Mol Biol 166:211–217
Lewis RNAH, Mannock DA, McElhaney RN (1989) Effect of fatty acyl chain length and structure on the lamellar gel to liquid-crystalline and lamellar to reversed hexagonal phase transitions of aqueous phosphatidylethanolamine dispersions. Biochemistry 28:541–548
Lindblom G, Wennerström H (1977) Biological and model membranes studies by nuclear magnetic resonance of spin one half nuclei. Q Rev Biophys 10:67–96
Lis LJ, McAlister M, Fuller N, Rand RP, Parsegian VA (1982) Interactions between neutral phospholipid bilayer membranes. Biophys J 37:657–665
London RE, Walker TE, Wilson DM, Matwiyoff NA (1979) Application of doubly decoupled carbon-13{proton, nitrogen-14-NMR spectroscopy to studies of phospholipids}. Chem Phys Lipids 25:7–14
Loof H de, Harvey SC, Segrest JP, Pastor RW (1991) Mean field stochastic boundary molecular dynamics simulation of a phospholipid in a membrane. Biochemistry 30:2099–2113
Marrink SJ, Berkowitz M, Berendsen HJC (1993) Molecular dynamics simulation of a membrane/water interface: the ordering of water and its relation to the hydration force. Langmuir (in press)
Marsh D (1974) Statistical mechanics of the fluidity of phospholipid bilayers and membranes. J Membr Biol 18:145–162
McConnell HM, Devaux P, Scandella C (1971) In: Fox CF (ed) Membrane research. Academic Press, New York
Meier P, Ohmes E, Kothe G (1986) Multipulse dynamic nuclear magnetic resonance of phospholipid membranes. J Chem Phys 85:3598–3614
Mendelsohn R, Davies MA, Brauner JW, Schuster HF, Dluhy RA (1989) Quantitative determination of conformational disorder in the acyl chains of phospholipid bilayers by infrared spectroscopy. Biochemistry 28:8934–8939
Meraldi JP, Schlitter J (1981) A statistical mechanical treatment of fatty acyl chain order in phospholipid bilayers and correlation with experimental data. B. Dipalmitoyl-3-sn-phosphatidylcholine. Biochim Biophys Acta 645:193–210
Moller MA, Tildesley DJ, Kim KS, Quirke N (1991) Molecular dynamics simulation of a Langmuir-Blodgett film. J Chem Phys 94:8390–8401
Nagle JF, Wiener MC (1988) Structure of fully hydrated bilayer dispersions. Biochim Biophys Acta 942:1–10
Phillips MC, Williams RM, Chapman D (1968) Hydrocarbon chain motions in lipid liquid crystals. Chem Phys Lipids 3:234–244
Pink DA, Green TJ, Chapman D (1980) Raman scattering in bilayers of saturated phosphatidylcholines. Experiment and theory. Biochemistry 19:349–356
Ploeg P van der, Berendsen HJC (1982) Molecular dynamics simulation of a bilayer membrane. J Chem Phys 76:3271–3276
Ploeg P van der, Berendsen HJC (1983) Molecular dynamics of a bilayer membrane. Mol Phys 49:233–248
Reiss-Husson J (1967) Structure des phases liquide-cristallines de différents phospholipides, monoglycérides, sphingolipides, anhydres on en presence d'eau. J Mol Biol 25:363–382
Rubinstein JLR, Smith BA, McConnell HM (1979) Lateral diffusion in binary mixtures of cholesterol and phosphatidylcholines. Proc Natl Acad Sci, USA 76:15–18
Ruocco MJ, Shipley GG (1982) Characterization of the sub-transition of hydrated dipalmitoylphosphatidylcholine bilayers. Biochim Biophys Acta 691:309–320
Ryckaert JP, Bellemans A (1975) Molecular dynamics of liquid n-butane near its boiling point. Chem Phys Lett 30:123–125
Ryckaert JP, Bellemans A (1978) Faraday Discuss Chem. Soc. 66:95
Ryckaert JP, Ciccotti G, Berendsen HJC (1977) Numerical integration of the cartesian equations of motion of a system with constraints: Molecular dynamics of n-alkanes. J Comp Phys 23:327–341
Schindler H, Seelig J (1975) Deuterium order parameters in relation to thermodynamic properties of a phospholipid bilayer. A statistical mechanical interpretation. Biochemistry 14:2283–2278
Seelig J (1977) Deuterium magnetic resonance: theory and application to lipid membranes. Q Rev Biophys 10:353–418
Seelig J (1978) Nuclear magnetic resonance and the head group structure of phospholipids in membranes. Biochim Biophys Acta 515:105–140
Seelig J, Gaily HU (1976) Investigation of phosphatidylethanolamine bilayers by deuterium and phosphorus-31 nuclear magnetic resonance. Biochemistry 15:5199–5204
Seelig A, Seelig J (1974) The dynamic structure of fatty acyl chains in a phospholipid bilayer measured by deuterium magnetic resonance. Biochemistry 13:4839–4845
Seelig A, Seelig J (1975) Bilayers of dipalmitoyl-3-sn-phosphatidylcholine conformational differences between the fatty acyl chains. Biochim Biophys Acta 406:1–5
Seelig J, Seelig A (1980) Lipid conformation in model membranes and biological membranes. Q Rev Biophys 13:19–61
Seelig J, Gaily HU, Wohlgemuth R (1977) Orientation and flexibility of the choline head group in phosphatidylcholine bilayers. Biochim Biophys Acta 467:109–119
Sheats JR, McConnell HM (1978) A photochemical technique for measuring lateral diffusion of spin-labeled phospholipids in membranes. Proc Natl Acad Sci, USA 75:4661–4663
Smith BA, McConnell HM (1978) Determination of molecular motion in membranes using periodic pattern photobleaching. Proc Natl Acad Sci, USA 15:2759–2763
Sundaralingam M (1972) Molecular structures and conformations of the phospholipids and sphingomyelins. Ann NY Acad Sci, USA 195:324–355
Tamm LK, Seelig J (1983) Lipid solvation of cytochrome c oxidase. Deuterium, nitrogen-14, and phosphorus-31 nuclear magnetic resonance studies on the phosphocholine head group and on cis-unsaturated fatty acyl chains. Biochemistry 22:1474–1483
Tardieu A, Luzzati V, Reman FC (1973) Structure and polymorphism of the hydrocarbon chains of lipids: a study of lecithinwater phases. J Mol Biol 75:711–733
Tenchov BG, Lis LJ, Quinn PJ (1987) Mechanism and kinetics of the subtransition in hydrated L-lipalmitoylphosphatidylcholine. Biochim Biophys Acta 897:143–151
Träuble H, Haynes DH (1974) Volume change in lipid bilayer lamellae at the crystalline-liquid crystalline phase transition. Chem Phys Lipids 7:324
Ulmius J, Wennerström H, Lindblom G, Arvidson G (1977) Deuteron nuclear magnetic resonance studies of phase equilibria in a lecithin-water system. Biochemistry 16:5742–5745
Vaz WLC, Stümpel J, Hallmann D, Gambacorta A, Rosa M de (1987) Bounding fluid viscosity and translational diffusion in a fluid lipid bilayer. Eur Biophys J 15:111–115
Wiener MC, White SH (1992a) Structure of a fluid dioleoylphosphatidylcholine bilayer determined by joint refinement of x-ray and neutron diffraction data. II. Distribution and packing of terminal methyl groups. Biophys J 61:428–433
Wiener MC, White SH (1992b) Structure of a fluid dioleoylphosphatidylcholine bilayer determined by joint refinement of x-ray and neutron diffraction data. I. Complete structure. Biophys J 61:434
Wu E, Jacobson K, Papahadjopoulos D (1977) Lateral diffusion in phospholipid multibilayers measured by fluorescence recovery after photobleaching. Biochemistry 16:3936–3941
Zaccai G, Büldt G, Seelig A, Seelig J (1979) Neutron diffraction studies on phosphatidylcholine model membranes. II Chain conformation and segmental disorder. J Mol Biol 134:693–706
Zaccai G, Blasie JK, Schoenborn BP (1980) Neutron diffraction studies on the location of water in lecithin bilayer model membranes. Proc Natl Acad Sci, USA 72:376–380