The nicotinic acetylcholine receptor: from molecular model to single-channel conductance

Electrochemical impedance spectroscopy of tethered bilayer membranes

The electrochemical impedance spectra (EIS) of tethered bilayer membranes (tBLMs) were analyzed, and the analytical solution for the spectral response of membranes containing natural or artificially introduced defects was derived. The analysis carried out in this work shows that the EIS features of an individual membrane defect cannot be modeled by conventional electrical elements. The primary reason for this is the complex nature of impedance of the submembrane ionic reservoir separating the phospholipid layer and the solid support. We demonstrate that its EIS response, in the case of radially symmetric defects, is described by the Hankel functions of a complex variable. Therefore, neither the impedance of the submembrane reservoir nor the total impedance of tBLMs can be modeled using the conventional elements of the equivalent electrical circuits of interfaces. There are, however, some limiting cases in which the complexity of the EIS response of the submembrane space reduces. In the high frequency limit, the EIS response of a submembrane space that surrounds the defect transforms into a response of a constant phase element (CPE) with the exponent (α) value of 0.5. The onset of this transformation is, beside other parameters, dependent on the defect size. Large-sized defects push the frequency limit lower, therefore, the EIS spectra exhibiting CPE behavior with α ≈ 0.5, can serve as a diagnostic criterion for the presence of such defects. In the low frequency limit, the response is dependent on the density of the defects, and it transforms into the capacitive impedance if the area occupied by a defect is finite. The higher the defect density, the higher the frequency edge at which the onset of the capacitive behavior is observed. Consequently, the presented analysis provides practical tools to evaluate the defect density in tBLMs, which could be utilized in tBLM-based biosensor applications. Alternatively, if the parameters of the defects, e.g., ion channels, such as the diameter and the conductance are known, the EIS data analysis provides a possibility to estimate other physical parameters of the system, such as thickness of the submembrane reservoir and its conductance. Finally, current analysis demonstrates a possibility to discriminate between the situations, in which the membrane defects are evenly distributed or clustered on the surface of tBLMs. Such sensitivity of EIS could be used for elucidation of the mechanisms of interaction between the proteins and the membranes.

Bookshelf: a simple curation system for the storage of biomolecular simulation data

Molecular dynamics simulations can now routinely generate data sets of several hundreds of gigabytes in size. The ability to generate this data has become easier over recent years and the rate of data production is likely to increase rapidly in the near future. One major problem associated with this vast amount of data is how to store it in a way that it can be easily retrieved at a later date. The obvious answer to this problem is a database. However, a key issue in the development and maintenance of such a database is its sustainability, which in turn depends on the ease of the deposition and retrieval process. Encouraging users to care about meta-data is difficult and thus the success of any storage system will ultimately depend on how well used by end-users the system is. In this respect we suggest that even a minimal amount of metadata if stored in a sensible fashion is useful, if only at the level of individual research groups. We discuss here, a simple database system which we call ‘Bookshelf’, that uses python in conjunction with a mysql database to provide an extremely simple system for curating and keeping track of molecular simulation data. It provides a user-friendly, scriptable solution to the common problem amongst biomolecular simulation laboratories; the storage, logging and subsequent retrieval of large numbers of simulations.

Download URL:http://sbcb.bioch.ox.ac.uk/bookshelf/

Mammalian nicotinic acetylcholine receptors: from structure to function

The classical studies of nicotine by Langley at the turn of the 20th century introduced the concept of a "receptive substance," from which the idea of a "receptor" came to light. Subsequent studies aided by the Torpedo electric organ, a rich source of muscle-type nicotinic receptors (nAChRs), and the discovery of alpha-bungarotoxin, a snake toxin that binds pseudo-irreversibly to the muscle nAChR, resulted in the muscle nAChR being the best characterized ligand-gated ion channel hitherto. With the advancement of functional and genetic studies in the late 1980s, the existence of nAChRs in the mammalian brain was confirmed and the realization that the numerous nAChR subtypes contribute to the psychoactive properties of nicotine and other drugs of abuse and to the neuropathology of various diseases, including Alzheimer's, Parkinson's, and schizophrenia, has since emerged. This review provides a comprehensive overview of these findings and the more recent revelations of the impact that the rich diversity in function and expression of this receptor family has on neuronal and nonneuronal cells throughout the body. Despite these numerous developments, our understanding of the contributions of specific neuronal nAChR subtypes to the many facets of physiology throughout the body remains in its infancy.

Dynamics of heteropentameric nicotinic acetylcholine receptor: implications of the gating mechanism

The dynamics characteristics of the currently available structure of Torpedo nicotinic acetylcholine receptor (nAChR), including the extracellular, transmembrane, and intracellular domains (ICDs), were analyzed using the Gaussian Network Model (GNM) and Anisotropic Network Model (ANM). We found that a symmetric quaternary twist motion, reported previously in the literature in a homopentameric receptor (Cheng et al. J Mol Biol 2006;355:310-324; Taly et al. Biophys J 2005;88:3954-3965), occurred also in the heteropentameric Torpedo nAChR. We believe, however, that the symmetric twist alone is not sufficient to explain a large body of experimental data indicating asymmetry and subunit nonequivalence during gating. Here we report our results supporting the hypothesis that a combination of symmetric and asymmetric motions opens the gate. We show that the asymmetric motion involves tilting of the TM2 helices. Furthermore, our study reveals three additional aspects of channel dynamics: (1) loop A serves as an allosteric mediator between the ligand binding loops and those at the domain interface, particularly the linker between TM2 and TM3; (2) the ICD can modulate the pore dynamics and thus should not be neglected in gating studies; and (3) the F loops, which are peculiarly longer and poorly-conserved in non-alpha-subunits, have important dynamical implications.

Blocking of the Nicotinic Acetylcholine Receptor Ion Channel by Chlorpromazine, a Noncompetitive Inhibitor: A Molecular Dynamics Simulation Study

A large series of pharmacological agents, distinct from the typical competitive antagonists, block in a noncompetitive manner the permeability response of the nicotinic acetylcholine receptor (nAChR) to the neurotransmitter acetylcholine. Taking the neuroleptic chlorpromazine (CPZ) as an example of such agents, the blocking mechanism of noncompetitive inhibitors to the ion channel pore of the nAChR has been explored at the atomic level using both conventional and steered molecular dynamics (MD) simulations. Repeated steered MD simulations have permitted calculation of the free energy (approximately 36 kJ/mol) of CPZ binding and identification of the optimal site in the region of the serine and leucine rings, at approximately 4 A from the pore entrance. Coulomb and the Lennard-Jones interactions between CPZ and the ion channel as well as the conformational fluctuations of CPZ were examined to assess the contribution of each to the binding of CPZ to the nAChR. The MD simulations disclose a dynamic interaction of CPZ binding to the nAChR ionic channel. The cationic ammonium head of CPZ forms strong hydrogen bonds with Glu262 (alpha), Asp268 (beta), Glu272 (beta), Ser276 (beta), Glu280 (delta), Gln271 (gamma), Glu275 (gamma), and Asn279 (gamma) nAChR residues. Finally, the conventional MD simulation of CPZ at its identified binding site demonstrates that the binding of CPZ not only blocks ion transport through the channel but also markedly inhibits the conformational transitions of the channel, necessary for nAChR to carry out its biological function.

Computed pore potentials of the nicotinic acetylcholine receptor

Electrostatic surface potentials in the vestibule of the nicotinic acetylcholine receptor (nAChR) were computed from structural models using the University of Houston Brownian Dynamics program to determine their effect on ion conduction and ionic selectivity. To further determine whether computed potentials accurately reflect the electrostatic environment of the channel, the potentials were used to predict the rate constants for diffusion-enhanced fluorescence energy transfer; the calculated energy transfer rates are directly comparable with those determined experimentally (see companion article by Meltzer et al. in this issue). To include any effects on the local potentials by the bound acceptor fluorophore crystal violet, its binding site was first localized within the pore by fluorescence energy transfer measurements from dansyl-C6-choline bound to the agonist sites and also by simulations of binding using Autodock. To compare the computed potentials with those determined experimentally, we used the predicted energy transfer rates from Tb3+ chelates of varying charge to calculate an expected potential using the Boltzmann relationship. This expected potential (from -20 to -40 mV) overestimates the values determined experimentally (from -10 to -25 mV) by two- to fourfold at similar conditions of ionic strength. Although the results indicate a basic discrepancy between experimental and computed surface potentials, both methods demonstrate that the vestibular potential has a relatively small effect on conduction and selectivity.

Electrolytic transport through a synthetic nanometer-diameter pore

We have produced single, synthetic nanometer-diameter pores by using a tightly focused, high-energy electron beam to sputter atoms in 10-nm-thick silicon nitride membranes. Subsequently, we measured the ionic conductance as a function of time, bath concentration, and pore diameter to infer the conductivity and ionic mobility through the pores. The pore conductivity is found to be much larger than the bulk conductivity for dilute bath concentrations, where the Debye length is larger than the pore radius, whereas it is comparable with or less than the bulk for high bath concentrations. We interpret these observations by using multiscale simulations of the ion transport through the pores. Molecular dynamics is used to estimate the ion mobility, and ion transport in the pore is described by the coupled Poisson-Nernst-Planck and the Stokes equations that are solved self-consistently for the ion concentration and velocity and electrical potential. We find that the measurements are consistent with the presence of fixed negative charge in the pore wall and a reduction of the ion mobility because of the fixed charge and the ion proximity to the pore wall.

Conformational dynamics of the nicotinic acetylcholine receptor channel: a 35-ns molecular dynamics simulation study

The nicotinic acetylcholine receptor (AChR) is the paradigm of ligand-gated ion channels, integral membrane proteins that mediate fast intercellular communication in response to neurotransmitters. A 35-ns molecular dynamics simulation has been performed to explore the conformational dynamics of the entire membrane-spanning region, including the ion channel pore of the AChR. In the simulation, the 20 transmembrane (TM) segments that comprise the whole TM domain of the receptor were inserted into a large dipalmitoylphosphatidylcholine (DPPC) bilayer. The dynamic behavior of individual TM segments and their corresponding AChR subunit helix bundles was examined in order to assess the contribution of each to the conformational transitions of the whole channel. Asymmetrical and asynchronous motions of the M1-M3 TM segments of each subunit were revealed. In addition, the outermost ring of five M4 TM helices was found to convey the effects exerted by the lipid molecules to the central channel domain. Remarkably, a closed-to-open conformational shift was found to occur in one of the channel ring positions in the time scale of the present simulations, the possible physiological significance of which is discussed.

Molecular dynamics simulation of the M2 helices within the nicotinic acetylcholine receptor transmembrane domain: structure and collective motions

Multiple nanosecond duration molecular dynamics simulations were performed on the transmembrane region of the Torpedo nicotinic acetylcholine receptor embedded within a bilayer mimetic octane slab. The M2 helices and M2-M3 loop regions were free to move, whereas the outer (M1, M3, M4) helix bundle was backbone restrained. The M2 helices largely retain their hydrogen-bonding pattern throughout the simulation, with some distortions in the helical end and loop regions. All of the M2 helices exhibit bending motions, with the hinge point in the vicinity of the central hydrophobic gate region (corresponding to residues alphaL251 and alphaV255). The bending motions of the M2 helices lead to a degree of dynamic narrowing of the pore in the region of the proposed hydrophobic gate. Calculations of Born energy profiles for various structures along the simulation trajectory suggest that the conformations of the M2 bundle sampled correspond to a closed conformation of the channel. Principal components analyses of each of the M2 helices, and of the five-helix M2 bundle, reveal concerted motions that may be relevant to channel function. Normal mode analyses using the anisotropic network model reveal collective motions similar to those identified by principal components analyses.

Asymmetric structural motions of the homomeric alpha7 nicotinic receptor ligand binding domain revealed by molecular dynamics simulation

A homology model of the ligand binding domain of the alpha7 nicotinic receptor is constructed based on the acetylcholine-binding protein crystal structure. This structure is refined in a 10 ns molecular dynamics simulation. The modeled structure proves fairly resilient, with no significant changes at the secondary or tertiary structural levels. The hypothesis that the acetylcholine-binding protein template is in the activated or desensitized state, and the absence of a bound agonist in the simulation suggests that the structure may also be relaxing from this state to the activatable state. Candidate motions that take place involve not only the side chains of residues lining the binding sites, but also the subunit positions that determine the overall shape of the receptor. In particular, two nonadjacent subunits move outward, whereas their partners counterclockwise to them move inward, leading to a marginally wider interface between themselves and an overall asymmetric structure. This in turn affects the binding sites, producing two that are more open and characterized by distinct side-chain conformations of W54 and L118, although motions of the side chains of all residues in every binding site still contribute to a reduction in binding site size, especially the outward motion of W148, which hinders acetylcholine binding. The Cys loop at the membrane interface also displays some flexibility. Although the short simulation timescale is unlikely to sample adequately all the conformational states, the pattern of observed motions suggests how ligand binding may correlate with larger-scale subunit motions that would connect with the transmembrane region that controls the passage of ions. Furthermore, the shape of the asymmetry with binding sites of differing affinity for acetylcholine, characteristic of other nicotinic receptors, may be a natural property of the relaxed, activatable state of alpha7.

Emerging structure of the nicotinic acetylcholine receptors

The conversion of acetylcholine binding into ion conduction across the membrane is becoming more clearly understood in terms of the structure of the receptor and its transitions. A high-resolution structure of a protein that is homologous to the extracellular domain of the receptor has revealed the binding sites and subunit interfaces in great detail. Although the structures of the membrane and cytoplasmic domains are less well determined, the channel lining and the determinants of selectivity have been mapped. The location and structure of the gates, and the coupling between binding sites and gates, remain to be established.

Side-chain ionization states in a potassium channel

KcsA is a bacterial K+ channel that is gated by pH. Continuum dielectric calculations on the crystal structure of the channel protein embedded in a low dielectric slab suggest that side chains E71 and D80 of each subunit, which lie adjacent to the selectivity filter region of the channel, form a proton-sharing pair in which E71 is neutral (protonated) and D80 is negatively charged at pH 7. When K+ ions are introduced into the system at their crystallographic positions the pattern of proton sharing is altered. The largest perturbation is for a K+ ion at site S3, i.e., interacting with the carbonyls of T75 and V76. The presence of multiple K+ ions in the filter increases the probability of E71 being ionized and of D80 remaining neutral (i.e., protonated). The ionization states of the protein side chains influence the potential energy profile experienced by a K+ ion as it is translated along the pore axis. In particular, the ionization state of the E71-D80 proton-sharing pair modulates the shape of the potential profile in the vicinity of the selectivity filter. Such reciprocal effects of ion occupancy on side-chain ionization states, and of side-chain ionization states on ion potential energy profiles will complicate molecular dynamics simulations and related studies designed to calculate ion permeation energetics.

Simulations of ion channels--watching ions and water move

Ion channels mediate electrical excitability in neurons and muscle. Three-dimensional structures for model peptide channels and for a potassium (K+) channel have been combined with computer simulations to permit rigorous exploration of structure-function relations of channels. Water molecules and ions within transbilayer pores tend to diffuse more slowly than in bulk solutions. In the narrow selectivity filter of the bacterial K+ channel (i.e. the region of the channel that discriminates between different species of ions) a column of water molecules and K+ ions moves in a concerted fashion. By combining atomistic simulations (in which all atoms of the channel molecule, water and ions are treated explicitly) with continuum methods (in which the description of the channel system is considerably simplified) it is possible to simulate some of the physiological properties of channels.