Structure of Gramicidin a in a Lipid Bilayer Environment Determined Using Molecular Dynamics Simulations and Solid-state NMR Data

Two different high-resolution structures recently have been proposed for the membrane-spanning gramicidin A channel: one based on solid-state NMR experiments in oriented phospholipid bilayers (Ketchem, R. R.; Roux, B.; Cross, T. A. Structure 1997, 5, 1655-1669; Protein Data Bank, PDB:1MAG); and one based on two-dimensional NMR in detergent micelles (Townsley, L. E.; Tucker, W. A.; Sham, S.; Hinton, J. F. Biochemistry 2001, 40, 11676-11686; PDB:1JNO). Despite overall agreement, the two structures differ in peptide backbone pitch and the orientation of several side chains; in particular that of the Trp at position 9. Given the importance of the peptide backbone and Trp side chains for ion permeation, we undertook an investigation of the two structures using molecular dynamics simulation with an explicit lipid bilayer membrane, similar to the system used for the solid-state NMR experiments. Based on 0.1 micros of simulation, both backbone structures converge to a structure with 6.25 residues per turn, in agreement with X-ray scattering, and broad agreement with SS backbone NMR observables. The side chain of Trp 9 is mobile, more so than Trp 11, 13, and 15, and undergoes spontaneous transitions between the orientations in 1JNO and 1MAG. Based on empirical fitting to the NMR results, and umbrella sampling calculations, we conclude that Trp 9 spends 80% of the time in the 1JNO orientation and 20% in the 1MAG orientation. These results underscore the utility of molecular dynamics simulations in the analysis and interpretation of structural information from solid-state NMR.