The pore domain of the nicotinic acetylcholine receptor: molecular modeling, pore dimensions, and electrostatics

Homologous mutations on different subunits cause unequal but additive effects on n-alcohol block in the nicotinic receptor pore

Hydrophobic antagonists of the nicotinic acetylcholine receptor inhibit channel activity by binding within the transmembrane pore formed by the second of four transmembrane domains (M2) on each of the receptor's subunits. Hydrophobic mutagenesis near the middle (10' locus) of the alpha-subunit M2 domain results in channels that are much more sensitive to block by long-chain alcohols and general anesthetics, indicating that the inhibitory site on wild-type receptors is nearby. To determine whether other receptor subunits also contribute to the blocker site, the hydrophobic mutagenesis strategy was extended to all four subunits at 10' loci. alpha S10'l causes the largest increase in apparent hexanol binding (4.3-fold compared to wild type), approximately twice the size of the change caused by beta T10'l (2.2-fold). gamma A10'l and delta A10'l mutations cause much smaller changes in apparent hexanol binding affinity (about 1.2-fold each), even when corrected for their smaller degree of side-chain hydrophobicity changes. When 10'l mutant subunits are coexpressed, the change from wild type in apparent hexanol binding energy (delta delta Gmixture) is roughly equal to the sum of hexanol binding energy changes for the constituent mutant subunits (sigma delta delta Gsubunits). The simplest model consistent with these results is one in which hydrophobic blockers make simultaneous contact with all five M2 10' residues, but the extent of contact is much greater for the alpha and beta than for gamma and delta side chains.

Alamethicin channels - modelling via restrained molecular dynamics simulations

Alamethicin channels have been modelled as approximately parallel bundles of transbilayer helices containing between N = 4 and 8 helices per bundle. Initial models were generated by in vacuo restrained molecular dynamics (MD) simulations, and were refined by 60 ps MD simulations with water molecules present within and at the mouths of the central pore. The helix bundles were stabilized by networks of H-bonds between intra-pore water molecules and Gln-7 side-chains. Channel conductances were predicted on the basis of pore radius profiles, and suggested that the N = 4 bundle formed an occluded pore, whereas pores with N > or = 5 helices per bundle were open. Continuum electrostatics calculations suggested that the N = 6 pore is cation-selective, whereas pores with N > or = 7 helices per bundle were predicted to be somewhat less ion-selective.

A novel method for structure-based prediction of ion channel conductance properties

A rapid and easy-to-use method of predicting the conductance of an ion channel from its three-dimensional structure is presented. The method combines the pore dimensions of the channel as measured in the HOLE program with an Ohmic model of conductance. An empirically based correction factor is then applied. The method yielded good results for six experimental channel structures (none of which were included in the training set) with predictions accurate to within an average factor of 1.62 to the true values. The predictive r2 was equal to 0.90, which is indicative of a good predictive ability. The procedure is used to validate model structures of alamethicin and phospholamban. Two genuine predictions for the conductance of channels with known structure but without reported conductances are given. A modification of the procedure that calculates the expected results for the effect of the addition of nonelectrolyte polymers on conductance is set out. Results for a cholera toxin B-subunit crystal structure agree well with the measured values. The difficulty in interpreting such studies is discussed, with the conclusion that measurements on channels of known structure are required.

Simulation studies of alamethicin-bilayer interactions

Alamethicin is an alpha-helical peptide that forms voltage-activated ion channels. Experimental data suggest that channel formation occurs via voltage-dependent insertion of alamethicin helices into lipid bilayers, followed by self-assembly of inserted helices to form a parallel helix bundle. Changes in the kink angle of the alamethicin helix about its central proline residue have also been suggested to play a role in channel gating. Alamethicin helices generated by simulated annealing and restrained molecular dynamics adopt a kink angle similar to that in the x-ray crystal structure, even if such simulations start with an idealized unkinked helix. This suggests that the kinked helix represents a stable conformation of the molecule. Molecular dynamics simulations in the presence of a simple bilayer model and a transbilayer voltage difference are used to explore possible mechanisms of helix insertion. The bilayer is represented by a hydrophobicity potential. An alamethicin helix inserts spontaneously in the absence of a transbilayer voltage. Application of a cis positive voltage decreases the time to insertion. The helix kink angle fluctuates during the simulations. Insertion of the helix is associated with a decrease in the mean kink angle, thus helping the alamethicin molecule to span the bilayer. The simulation results are discussed in terms of models of alamethicin channel gating.

HOLE: a program for the analysis of the pore dimensions of ion channel structural models

A method (HOLE) that allows the analysis of the dimensions of the pore running through a structural model of an ion channel is presented. The algorithm uses a Monte Carlo simulated annealing procedure to find the best route for a sphere with variable radius to squeeze through the channel. Results can be displayed in a graphical fashion or visualized with most common molecular graphical packages. Advances include a method to analyze the anisotropy within a pore. The method can also be used to predict the conductance of channels using a simple empirically corrected ohmic model. As an example the program is applied to the cholera toxin B-subunit pentamer. The compatibility of the crystal structure and conductance data is established.