Open AccessMolecular dynamics simulation of the gramicidin channel in a phospholipid bilayer.
Author(s) -
Thomas B. Woolf,
Benoı̂t Roux
Publication year - 1994
Publication title -
proceedings of the national academy of sciences of the united states of america
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 5.011
H-Index - 771
eISSN - 1091-6490
pISSN - 0027-8424
DOI - 10.1073/pnas.91.24.11631
Subject(s) - molecular dynamics , lipid bilayer , gramicidin , chemistry , van der waals force , phospholipid , chemical physics , bilayer , hydrogen bond , membrane , phosphatidylcholine , crystallography , biophysics , computational chemistry , molecule , organic chemistry , biochemistry , biology
A molecular dynamics simulation of the gramicidin A channel in an explicit dimyristoyl phosphatidylcholine bilayer was generated to study the details of lipid-protein interactions at the microscopic level. Solid-state NMR properties of the channel averaged over the 500-psec trajectory are in excellent agreement with available experimental data. In contrast with the assumptions of macroscopic models, the membrane/solution interface region is found to be at least 12 A thick. The tryptophan side chains, located within the interface, are found to form hydrogen bonds with the ester carbonyl groups of the lipids and with water, suggesting their important contribution to the stability of membrane proteins. Individual lipid-protein interactions are seen to vary from near 0 to -50 kcal/mol. The most strongly interacting conformations are short-lived and have a nearly equal contribution from both van der Waals and electrostatic energies. This approach for performing molecular dynamics simulations of membrane proteins in explicit phospholipid bilayers should help in studying the structure, dynamics, and energetics of lipid-protein interactions.