Modelling membrane proteins using structural restraints.

Here we present a procedure for modelling membrane proteins which employs molecular dynamics simulations incorporating target restraints derived from low-resolution structures alongside distance restraints derived from mutagenesis data. The application of the modelling procedure to the closed conformation of the pore domain of the nicotinic acetylcholine receptor is described. This domain is formed by a parallel bundle of five M2 helices. Each M2 helix is kinked due to cumulative distortions of backbone (phi, psi) values. The central region of M2 may adopt a more distorted conformation. This would enable a ring of conserved leucine residues (one from each M2 helix) to pack together, occluding the central pore and thus preventing ion permeation. Molecular dynamics simulations on isolated helices that kink formation is not an inherent property of M2.