Gramicidin A backbone and side chain dynamics evaluated by molecular dynamics simulations and nuclear magnetic resonance experiments. I: Molecular dynamics simulations

Helgi I. Ingolfsson, Yuhui Li, Vitaly V. Vostrikov, Hong Gu, James F. Hinton, Roger E. Koeppe, Benoît Roux, Olaf S. Andersen

Research output: Contribution to journalArticle

25 Scopus citations

Abstract

Gramicidin A (gA) channels provide an ideal system to test molecular dynamics (MD) simulations of membrane proteins. The peptide backbone lines a cation-selective pore, and due to the small channel size, the average structure and extent of fluctuations of all atoms in the peptide will influence ion permeation. This raises the question of how well molecular mechanical force fields used in MD simulations and potential of mean force (PMF) calculations can predict structure and dynamics as well as ion permeation. To address this question, we undertook a comparative study of nuclear magnetic resonance (NMR) observables predicted by fully atomistic MD simulations on a gA dimer embedded in a sodium dodecyl sulfate (SDS) micelle with measurements of the gA dimer backbone and tryptophan side chain dynamics using solution-state 15N NMR on gA dimers in SDS micelles (Vostrikov, V. V.; Gu, H.; Ingólfsson, H. I.; Hinton, J. F.; Andersen, O. S.; Roux, B.; Koeppe, R. E., II. J. Phys. Chem. B2011, DOI 10.1021/jp200906y, accompanying article). This comparison enables us to examine the robustness of the MD simulations done using different force fields as well as their ability to predict important features of the gA channel. We find that MD is able to predict NMR observables, including the generalized order parameters (S2), the 15N spin-lattice (T1) and spin-spin (T2) relaxation times, and the 1H-15N nuclear Overhauser effect (NOE), with remarkable accuracy. To examine further how differences in the force fields can affect the channel conductance, we calculated the PMF for K+ and Na+ permeation through a gA channel in a dimyristoylphosphatidylcholine (DMPC) bilayer. In this case, we find that MD is less successful in quantitatively predicting the single-channel conductance.

Original languageEnglish (US)
Pages (from-to)7417-7426
Number of pages10
JournalJournal of Physical Chemistry B
Volume115
Issue number22
DOIs
StatePublished - Jun 9 2011

Fingerprint Dive into the research topics of 'Gramicidin A backbone and side chain dynamics evaluated by molecular dynamics simulations and nuclear magnetic resonance experiments. I: Molecular dynamics simulations'. Together they form a unique fingerprint.

  • Cite this