The acetylcholine-binding protein (AChBP) is homologous to the ligand-binding domain of the nicotinic acetylcholine receptor (nAChR) and other members of the Cys-loop family of neurotransmitter receptors. The high-resolution X-ray structures of AChBP mean it has been used as a model from which to understand agonist and antagonist binding to nAChRs. We present here a molecular dynamics (MD) study of AChBP with nicotine and carbamylcholine bound. Our results suggest that the ligand imposes rigidity on the binding pocket residues. The simulations also suggest that the protein undergoes breathing motions with respect to the five-fold axis, a motion that has been postulated to be related to gating in the nAChR. We analyzed the behaviour of the water molecules in and around the binding site and found that they occupied five distinct sites within the binding pocket. Water occupied these sites in the absence of ligand, but the presence of ligand increased the probability that a water molecule would be found in these sites. Finally, we demonstrate how the positions of these waters might be used in the design of new ligands by comparing the positions of these sites with other recent structures.

The structure of a homopentameric alpha7 nicotinic acetylcholine receptor is modelled by combining structural information from two sources: the X-ray structure of a water soluble acetylcholine binding protein from Lymnea stagnalis, and the electron microscopy derived structure of the transmembrane domain of the Torpedo nicotinic receptor. The alpha7 nicotinic receptor model is generated by simultaneously optimising: (i) chain connectivity, (ii) avoidance of stereochemically unfavourable contacts, and (iii) contact between the beta1-beta2 and M2-M3 loops that have been suggested to be involved in transmission of conformational change between the extracellular and transmembrane domains. A Gaussian network model was used to predict patterns of residue mobility in the alpha7 model. The results of these calculations suggested a flexibility gradient along the transmembrane domain, with the extracellular end of the domain more flexible that the intracellular end. Poisson-Boltzmann (PB) energy calculations and atomistic (molecular dynamics) simulations were used to estimate the free energy profile of a Na+ ion as a function of position along the axis of the pore-lining M2 helix bundle of the transmembrane domain. Both types of calculation suggested a significant energy barrier to exist in the centre of the (closed) pore, consistent with a "hydrophobic gating" model. Estimations of the PB energy profile as a function of ionic strength suggest a role of the extracellular domain in determining the cation selectivity of the alpha7 nicotinic receptor. These studies illustrate how molecular models of members of the nicotinic receptor superfamily of channels may be used to study structure-function relationships.

Nicotinic acetylcholine receptors (nAChRs) are ligand-gated ion channels mediating fast cholinergic synaptic transmission in the brain and at neuromuscular junctions. We used the structure of the acetylcholine binding protein from Lymnaea stagnalis to model the chicken ␣7 agonist-binding domain. The initial models and a preliminary docking study suggested that position Leu118 may play an important role in determining agonist actions on ␣7. A prediction from these in silico studies, that L118E and L118D would retain binding to acetylcholine but L118K and L118R would not, was confirmed in electrophysiological studies on functional recombinant mutant receptors expressed in Xenopus laevis oocytes. The functional studies also demonstrated that residues at position 118 have a dramatic effect on the actions of imidacloprid (a partial agonist of wild-type ␣7 receptors) and its des-nitro derivative. Molecular dynamics simulations confirmed that Leu118 can strongly influence agonist binding and that the model was robust in terms of its prediction for acetylcholine binding. Together, the results indicate a role for Leu118 in influencing agonist actions on ␣7 nAChRs.