Shaker pore structure as predicted by annealed atomic simulation using symmetry and novel geometric restraints.

Recent studies making use of channel-blocking peptides as molecular calipers have revealed the architecture of the pore-forming region of Shaker-type potassium channels. Here we show that the low-resolution, experimentally derived geometric information can be incorporated as restraints within the context of an annealed molecular dynamics simulation to predict an atomic structure for the channel pore which, by virtue of restraints, conforms to the experimental evidence. The simulation is reminiscent of the computational method employed by nuclear magnetic resonance (NMR) spectroscopists to resolve solution structures of biological macromolecules, but in lieu of restraints conventionally derived from NMR spectra, novel restraints are developed that include side-chain orientation of amino acid residues and assumed symmetry of protein subunits. The method presented here offers the possibility of expanding cooperation between simulation and experiment in developing structural models, especially for systems such as ion channels whose three-dimensional structures may not be amenable to determination by direct methods at the present time.