This Article Full Text Full Text (PDF) All Versions of this Article: biophysj.106.084483v1 91/10/3617    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Benz, R. W. Articles by Tobias, D. J. Search for Related Content PubMed PubMed Citation Articles by Benz, R. W. Articles by Tobias, D. J.

Diffraction-Based Density Restraints for Membrane and Membrane-Peptide Molecular Dynamics Simulations

Ryan W. Benz ^*, Hirsh Nanda  , Francisco Castro-Román , Stephen H. White  and Douglas J. Tobias ^*

^* Department of Chemistry, and Department of Physiology and Biophysics, University of California, Irvine, California; and NIST Center for Neutron Research, National Institute of Standards and Technology, Gaithersburg, Maryland

Correspondence: Address reprint requests to Douglas J. Tobias, E-mail: dtobias{at}uci.edu; or Stephen H. White, E-mail: stephen.white{at}uci.edu.

We have recently shown that current molecular dynamics (MD) atomic force fields are not yet able to produce lipid bilayer structures that agree with experimentally-determined structures within experimental errors. Because of the many advantages offered by experimentally validated simulations, we have developed a novel restraint method for membrane MD simulations that uses experimental diffraction data. The restraints, introduced into the MD force field, act upon specified groups of atoms to restrain their mean positions and widths to values determined experimentally. The method was first tested using a simple liquid argon system, and then applied to a neat dioleoylphosphatidylcholine (DOPC) bilayer at 66% relative humidity and to the same bilayer containing the peptide melittin. Application of experiment-based restraints to the transbilayer double-bond and water distributions of neat DOPC bilayers led to distributions that agreed with the experimental values. Based upon the experimental structure, the restraints improved the simulated structure in some regions while introducing larger differences in others, as might be expected from imperfect force fields. For the DOPC-melittin system, the experimental transbilayer distribution of melittin was used as a restraint. The addition of the peptide caused perturbations of the simulated bilayer structure, but which were larger than observed experimentally. The melittin distribution of the simulation could be fit accurately to a Gaussian with parameters close to the observed ones, indicating that the restraints can be used to produce an ensemble of membrane-bound peptide conformations that are consistent with experiments. Such ensembles pave the way for understanding peptide-bilayer interactions at the atomic level.

This article has been cited by other articles:

				C. L. Wee, K. Balali-Mood, D. Gavaghan, and M. S. P. Sansom The Interaction of Phospholipase A2 with a Phospholipid Bilayer: Coarse-Grained Molecular Dynamics Simulations Biophys. J., August 15, 2008; 95(4): 1649 - 1657. [Abstract] [Full Text] [PDF]

				S. H. White Membrane Protein Insertion: The Biology-Physics Nexus J. Gen. Physiol., April 30, 2007; 129(5): 363 - 369. [Full Text] [PDF]