A mechanism for acetylcholine receptor gating based on structure, coupling, phi, and flip

The effect of unliganded gating on agonist response in nicotinic receptors

Understanding the agonist concentration-response curve (CRC) is the cornerstone in pharmacology. While CRC parameters, agonist potency (EC50) and efficacy (maximum response, Imax) are well-studied, the role of unliganded gating (minimum response, Imin) on CRC is often overlooked. This study explores the effect of unliganded gating on agonist response in muscle-type acetylcholine (ACh) receptors, focusing on the underexplored role of Imin in modulating EC50 and Imax. Three Gain-of-Function (GOF) mutations that increase, and two Loss-of-Function (LOF) mutations that decrease the unliganded gating equilibrium constant (L0) were studied using automated patch-clamp electrophysiology. GOF mutations enhanced agonist potency, whereas LOF mutations reduced it. The calculated CRC aligned well with empirical results, indicating that agonist CRC can be estimated from knowledge of L0. Reduction in agonist efficacy due to LOF mutations was calculated and subsequently validated using single-channel patch-clamp electrophysiology, a factor often obscured in normalized CRC. The study also evaluated the combined impact of mutations (L0) on CRC, confirming the predictive model. Further, no significant energetic coupling between distant residues (>15 Å) was found, indicating that the mutations' effects are localized and do not alter overall agonist affinity. These findings substantiate the role of unliganded gating in modulating agonist responses and establishes a predictive model for estimating CRC parameters from known changes in L0. The study highlights the importance of intrinsic activity in receptor theory.

The Nicotinic Acetylcholine Receptor and Its Pentameric Homologs: Toward an Allosteric Mechanism of Signal Transduction at the Atomic Level

The nicotinic acetylcholine receptor has served, since its biochemical identification in the 1970s, as a model of an allosteric ligand-gated ion channel mediating signal transition at the synapse. In recent years, the application of X-ray crystallography and high-resolution cryo-electron microscopy, together with molecular dynamic simulations of nicotinic receptors and homologs, have opened a new era in the understanding of channel gating by the neurotransmitter. They reveal, at atomic resolution, the diversity and flexibility of the multiple ligand-binding sites, including recently discovered allosteric modulatory sites distinct from the neurotransmitter orthosteric site, and the conformational dynamics of the activation process as a molecular switch linking these multiple sites. The model emerging from these studies paves the way for a new pharmacology based, first, upon the occurrence of an original mode of indirect allosteric modulation, distinct from a steric competition for a single and rigid binding site, and second, the design of drugs that specifically interact with privileged conformations of the receptor such as agonists, antagonists, and desensitizers. Research on nicotinic receptors is still at the forefront of understanding the mode of action of drugs on the nervous system.

Dynamics of receptor activation by agonists

How do agonists turn on receptors? The model system we have used to address this question is the adult-type skeletal muscle nicotinic acetylcholine receptor. This ligand-gated ion channel has two orthosteric sites (for neurotransmitters) in the extracellular domain linked to an allosteric site (a gate) in the transmembrane domain. The goal of this perspective is to summarize how measurements of agonist binding energy reveal the dynamics of the neurotransmitter sites and the fundamental link between binding and gating.

A release of local subunit conformational heterogeneity underlies gating in a muscle nicotinic acetylcholine receptor

Synaptic receptors respond to neurotransmitters by opening an ion channel across the post-synaptic membrane to elicit a cellular response. Here we use recent Torpedo acetylcholine receptor structures and functional measurements to delineate a key feature underlying allosteric communication between the agonist-binding extracellular and channel-gating transmembrane domains. Extensive mutagenesis at this inter-domain interface re-affirms a critical energetically coupled role for the principal α subunit β1-β2 and M2-M3 loops, with agonist binding re-positioning a key β1-β2 glutamate/valine to facilitate the outward motions of a conserved M2-M3 proline to open the channel gate. Notably, the analogous structures in non-α subunits adopt a locally active-like conformation in the apo state even though each L9' hydrophobic gate residue in each pore-lining M2 α-helix is closed. Agonist binding releases local conformational heterogeneity transitioning all five subunits into a conformationally symmetric open state. A release of conformational heterogeneity provides a framework for understanding allosteric communication in pentameric ligand-gated ion channels.

Mechanism of hydrophobic gating in the acetylcholine receptor channel pore

Neuromuscular acetylcholine receptors (AChRs) are hetero-pentameric, ligand-gated ion channels. The binding of the neurotransmitter acetylcholine (ACh) to two target sites promotes a global conformational change of the receptor that opens the channel and allows ion conduction through the channel pore. Here, by measuring free-energy changes from single-channel current recordings and using molecular dynamics simulations, we elucidate how a constricted hydrophobic region acts as a "gate" to regulate the channel opening in the pore of AChRs. Mutations of gate residues, including those implicated in congenital myasthenia syndrome, lower the permeation barrier of the channel substantially and increase the unliganded gating equilibrium constant (constitutive channel openings). Correlations between hydrophobicity and the observed free-energy changes, supported by calculations of water densities in the wild-type versus mutant channel pores, provide evidence for hydrophobic wetting-dewetting transition at the gate. The analysis of a coupled interaction network provides insight into the molecular mechanism of closed- versus open-state conformational changes at the gate. Studies of the transition state by "phi"(φ)-value analysis indicate that agonist binding serves to stabilize both the transition and the open state. Intersubunit interaction energy measurements and molecular dynamics simulations suggest that channel opening involves tilting of the pore-lining M2 helices, asymmetric outward rotation of amino acid side chains, and wetting transition of the gate region that lowers the barrier to ion permeation and stabilizes the channel open conformation. Our work provides new insight into the hydrophobic gate opening and shows why the gate mutations result in constitutive AChR channel activity.

Mechanistic basis of ligand efficacy in the calcium-activated chloride channel TMEM16A

Agonist binding in ligand-gated ion channels is coupled to structural rearrangements around the binding site, followed by the opening of the channel pore. In this process, agonist efficacy describes the equilibrium between open and closed conformations in a fully ligand-bound state. Calcium-activated chloride channels in the TMEM16 family are important sensors of intracellular calcium signals and are targets for pharmacological modulators, yet a mechanistic understanding of agonist efficacy has remained elusive. Using a combination of cryo-electron microscopy, electrophysiology, and autocorrelation analysis, we now show that agonist efficacy in the ligand-gated channel TMEM16A is dictated by the conformation of the pore-lining helix α6 around the Ca2+ -binding site. The closure of the binding site, which involves the formation of a π-helix below a hinge region in α6, appears to be coupled to the opening of the inner pore gate, thereby governing the channel's open probability and conductance. Our results provide a mechanism for agonist binding and efficacy and a structural basis for the design of potentiators and partial agonists in the TMEM16 family.

Agonist efficiency links binding and gating in a nicotinic receptor

Receptors signal by switching between resting (C) and active (O) shapes ('gating') under the influence of agonists. The receptor's maximum response depends on the difference in agonist binding energy, O minus C. In nicotinic receptors, efficiency (η) represents the fraction of agonist binding energy applied to a local rearrangement (an induced fit) that initiates gating. In this receptor, free energy changes in gating and binding can be interchanged by the conversion factor η. Efficiencies estimated from concentration-response curves (23 agonists, 53 mutations) sort into five discrete classes (%): 0.56 (17), 0.51(32), 0.45(13), 0.41(26), and 0.31(12), implying that there are 5 C versus O binding site structural pairs. Within each class efficacy and affinity are corelated linearly, but multiple classes hide this relationship. η unites agonist binding with receptor gating and calibrates one link in a chain of coupled domain rearrangements that comprises the allosteric transition of the protein.

Mutation of valine 53 at the interface between extracellular and transmembrane domains of the β2 principal subunit affects the GABAA receptor gating

GABA_A receptors (gamma-aminobutyric acid type A receptors) are pentameric ligand-gated ion channels mediating inhibition in adult mammalian brains. Their static structure has been intensely studied in the past years but the underlying molecular activatory mechanisms remain obscure. The interface between extracellular and transmembrane domains has been recognized as a key player in the receptor gating. However, the role of the valine 53 in the β1-β2 loop of the principal subunit (β₂) remains controversial showing differences compared to homologous residues in some cys-loop counterparts such as nAChR. To address the role of the β₂V53 residue in the α₁β₂γ_2L receptor gating, we performed high resolution macroscopic and single-channel recordings. To explore underlying molecular mechanisms a variety of substituting amino acids were investigated: Glutamate and Lysine (different electric charge), Alanine (aliphatic, larger than Valine) and Histidine (same residue as in homologous α₁H55). We report that mutation of the β₂V53 residue results in alterations of nearly all gating transitions including opening/closing, preactivation and desensitization. A dramatic gating impairment was observed for glutamate substitution (β₂V53E) but β₂V53K mutation had a weak effect. The impact of histidine substitution was also small while β₂V53A markedly affected the receptor but to a smaller extent than β₂V53E. Considering available structures in desensitized and bicuculline blocked shut states we propose that strongly detrimental effect of β₂V53E mutation on receptor activation results from electrostatic interaction between the glutamate and β₂K274 on the loop M2-M3 which stabilizes the receptor in the shut state. We conclude that β₂V53 is strongly involved in mechanisms underlying the receptor gating.

Illumination of a progressive allosteric mechanism mediating the glycine receptor activation

Pentameric ligand-gated ion channel mediate signal transduction at chemical synapses by transiting between resting and open states upon neurotransmitter binding. Here, we investigate the gating mechanism of the glycine receptor fluorescently labeled at the extracellular-transmembrane interface by voltage-clamp fluorometry (VCF). Fluorescence reports a glycine-elicited conformational change that precedes pore opening. Low concentrations of glycine, partial agonists or specific mixtures of glycine and strychnine trigger the full fluorescence signal while weakly activating the channel. Molecular dynamic simulations of a partial agonist bound-closed Cryo-EM structure show a highly dynamic nature: a marked structural flexibility at both the extracellular-transmembrane interface and the orthosteric site, generating docking properties that recapitulate VCF data. This work illuminates a progressive propagating transition towards channel opening, highlighting structural plasticity within the mechanism of action of allosteric effectors.

Recent Insight into Lipid Binding and Lipid Modulation of Pentameric Ligand-Gated Ion Channels

Pentameric ligand-gated ion channels (pLGICs) play a leading role in synaptic communication, are implicated in a variety of neurological processes, and are important targets for the treatment of neurological and neuromuscular disorders. Endogenous lipids and lipophilic compounds are potent modulators of pLGIC function and may help shape synaptic communication. Increasing structural and biophysical data reveal sites for lipid binding to pLGICs. Here, we update our evolving understanding of pLGIC–lipid interactions highlighting newly identified modes of lipid binding along with the mechanistic understanding derived from the new structural data.

Antagonistic Mechanism of α-Conotoxin BuIA toward the Human α3β2 Nicotinic Acetylcholine Receptor

Nicotinic acetylcholine receptors (nAChRs) are pentameric ligand-gated ion channels that are abundantly expressed in the central and peripheral nervous systems, playing an important role in mediating neurotransmitter release and inter-synaptic signaling. Dysfunctional nAChRs are associated with neurological disorders, and studying the structure and function of nAChRs is essential for development of drugs or strategies for treatment of related diseases. α-Conotoxins are selective antagonists of the nAChR and are an important class of drug leads. So far, the antagonistic mechanism of α-conotoxins toward the nAChRs is still unclear. In this study, we built an α3β2 nAChR homology model and investigated its conformational transition mechanism upon binding with a highly potent inhibitor, α-conotoxin BuIA, through μs molecular dynamic simulations and site-directed mutagenesis studies. The results suggested that the α3β2 nAChR underwent global conformational transitions and was stabilized into a closed state with three hydrophobic gates present in the transmembrane domain by BuIA. Finally, the probable antagonistic mechanism of BuIA was proposed. Overall, the closed-state model of the α3β2 nAChR bound with BuIA is not only essential for understanding the antagonistic mechanism of α-conotoxins but also particularly valuable for development of therapeutic inhibitors in future.

Antagonistic postsynaptic and presynaptic actions of cyclohexanol on neuromuscular synaptic transmission and function

Intentional ingestion of agricultural organophosphorus insecticides is a significant public health issue in rural Asia, causing thousands of deaths annually. Some survivors develop a severe, acute or delayed myasthenic syndrome. In animal models, similar myasthenia has been associated with increasing plasma concentration of one insecticide solvent metabolite, cyclohexanol. We investigated possible mechanisms using voltage and current recordings from mouse neuromuscular junctions (NMJs) and transfected human cell lines. Cyclohexanol (10-25 mM) reduced endplate potential (EPP) amplitudes by 10-40% and enhanced depression during repetitive (2-20 Hz) stimulation by up to 60%. EPP decay was prolonged more than twofold. Miniature EPPs were attenuated by more than 50%. Cyclohexanol inhibited whole-cell currents recorded from CN21 cells expressing human postjunctional acetylcholine receptors (hnAChR) with an IC₅₀ of 3.74 mM. Cyclohexanol (10-20 mM) also caused prolonged episodes of reduced-current, multi-channel bursting in outside-out patch recordings from hnAChRs expressed in transfected HEK293T cells, reducing charge transfer by more than 50%. Molecular modelling indicated cyclohexanol binding (-6 kcal/mol) to a previously identified alcohol binding site on nicotinic AChR α-subunits. Cyclohexanol also increased quantal content of evoked transmitter release by ∼50%. In perineurial recordings, cyclohexanol selectively inhibited presynaptic K⁺ currents. Modelling indicated cyclohexanol binding (-3.8 kcal/mol) to voltage-sensitive K⁺ channels at the same site as tetraethylammonium (TEA). TEA (10 mM) blocked K⁺ channels more effectively than cyclohexanol but EPPs were more prolonged in 20 mM cyclohexanol. The results explain the pattern of neuromuscular dysfunction following ingestion of organophosphorus insecticides containing cyclohexanol precursors and suggest that cyclohexanol may facilitate investigation of mechanisms regulating synaptic strength at NMJs. KEY POINTS: Intentional ingestion of agricultural organophosphorus insecticides is a significant public health issue in rural Asia, causing thousands of deaths annually. Survivors may develop a severe myasthenic syndrome or paralysis, associated with increased plasma levels of cyclohexanol, an insecticide solvent metabolite. Analysis of synaptic transmission at neuromuscular junctions in isolated mouse skeletal muscle, using isometric tension recording and microelectrode recording of endplate voltages and currents, showed that cyclohexanol reduced postsynaptic sensitivity to acetylcholine neurotransmitter (reduced quantal size) while simultaneously enhancing evoked transmitter release (increased quantal content). Patch recording from transfected cell lines, together with molecular modelling, indicated that cyclohexanol causes selective, allosteric antagonism of postsynaptic nicotinic acetylcholine receptors and block of presynaptic K⁺ -channel function. The data provide insight into the cellular and molecular mechanisms of neuromuscular weakness following intentional ingestion of agricultural organophosphorus insecticides. Our findings also extend understanding of the effects of alcohols on synaptic transmission and homeostatic synaptic function.

Towards the Idea of Molecular Brains

How can single cells without nervous systems perform complex behaviours such as habituation, associative learning and decision making, which are considered the hallmark of animals with a brain? Are there molecular systems that underlie cognitive properties equivalent to those of the brain? This review follows the development of the idea of molecular brains from Darwin’s “root brain hypothesis”, through bacterial chemotaxis, to the recent discovery of neuron-like r-protein networks in the ribosome. By combining a structural biology view with a Bayesian brain approach, this review explores the evolutionary labyrinth of information processing systems across scales. Ribosomal protein networks open a window into what were probably the earliest signalling systems to emerge before the radiation of the three kingdoms. While ribosomal networks are characterised by long-lasting interactions between their protein nodes, cell signalling networks are essentially based on transient interactions. As a corollary, while signals propagated in persistent networks may be ephemeral, networks whose interactions are transient constrain signals diffusing into the cytoplasm to be durable in time, such as post-translational modifications of proteins or second messenger synthesis. The duration and nature of the signals, in turn, implies different mechanisms for the integration of multiple signals and decision making. Evolution then reinvented networks with persistent interactions with the development of nervous systems in metazoans. Ribosomal protein networks and simple nervous systems display architectural and functional analogies whose comparison could suggest scale invariance in information processing. At the molecular level, the significant complexification of eukaryotic ribosomal protein networks is associated with a burst in the acquisition of new conserved aromatic amino acids. Knowing that aromatic residues play a critical role in allosteric receptors and channels, this observation suggests a general role of π systems and their interactions with charged amino acids in multiple signal integration and information processing. We think that these findings may provide the molecular basis for designing future computers with organic processors.

Nicotinic receptors: From protein allostery to computational neuropharmacology

We propose an extension and further development of the Monod-Wyman-Changeux model for allosteric transitions of regulatory proteins to brain communications and specifically to neurotransmitters receptors, with the nicotinic acetylcholine receptor (nAChR) as a model of ligand-gated ion channels. The present development offers an expression of the change of the gating isomerization constant caused by pharmacological ligand binding in terms of its value in the absence of ligands and several "modulation factors", which vary with orthosteric ligand binding (agonists/antagonists), allosteric ligand binding (positive allosteric modulators/negative allosteric modulators) and receptor desensitization. The new - explicit - formulation of such "modulation factors", provides expressions for the pharmacological attributes of potency, efficacy, and selectivity for the modulatory ligands (including endogenous neurotransmitters) in terms of their binding affinity for the active, resting, and desensitized states of the receptor. The current formulation provides ways to design neuroactive compounds with a controlled pharmacological profile, opening the field of computational neuro-pharmacology.

Mutational analysis to explore long-range allosteric couplings involved in a pentameric channel receptor pre-activation and activation

Pentameric ligand-gated ion channels (pLGICs) mediate chemical signaling through a succession of allosteric transitions that are yet not completely understood as intermediate states remain poorly characterized by structural approaches. In a previous study on the prototypic bacterial proton-gated channel GLIC, we generated several fluorescent sensors of the protein conformation that report a fast transition to a pre-active state, which precedes the slower process of activation with pore opening. Here, we explored the phenotype of a series of allosteric mutations, using simultaneous steady-state fluorescence and electrophysiological measurements over a broad pH range. Our data, fitted to a three-state Monod-Wyman-Changeux model, show that mutations at the subunit interface in the extracellular domain (ECD) principally alter pre-activation, while mutations in the lower ECD and in the transmembrane domain principally alter activation. We also show that propofol alters both transitions. Data are discussed in the framework of transition pathways generated by normal mode analysis (iModFit). It further supports that pre-activation involves major quaternary compaction of the ECD, and suggests that activation involves principally a reorganization of a ‘central gating region’ involving a contraction of the ECD β-sandwich and the tilt of the channel lining M2 helix.

Pathways for nicotinic receptor desensitization

Nicotinic acetylcholine receptors (AChRs) are ligand-gated ion channels that generate transient currents by binding agonists and switching rapidly between closed- and open-channel conformations. Upon sustained exposure to ACh, the cell response diminishes slowly because of desensitization, a process that shuts the channel even with agonists still bound. In liganded receptors, the main desensitization pathway is from the open-channel conformation, but after agonists dissociate the main recovery pathway is to the closed-channel conformation. In this Viewpoint, I discuss two mechanisms that can explain the selection of different pathways, a question that has puzzled the community for 60 yr. The first is based on a discrete-state model (the “prism”), in which closed, open, and desensitized conformational states interconnect directly. This model predicts that 5% of unliganded AChRs are desensitized. Different pathways are taken with versus without agonists because ligands have different energy properties (φ values) at the transition states of the desensitization and recovery reactions. The second is a potential energy surface model (the “monkey saddle”), in which the states connect indirectly at a shared transition state region. Different pathways are taken because agonists shift the position of the gating transition state relative to the point where gating and desensitization conformational trajectories intersect. Understanding desensitization pathways appears to be a problem of kinetics rather than of thermodynamics. Other aspects of the two mechanisms are considered, as are experiments that may someday distinguish them.

Agonist efficiency from concentration-response curves: Structural implications and applications

Agonists are evaluated by a concentration-response curve (CRC), with a midpoint (EC₅₀) that indicates potency, a high-concentration asymptote that indicates efficacy, and a low-concentration asymptote that indicates constitutive activity. A third agonist attribute, efficiency (η), is the fraction of binding energy that is applied to the conformational change that activates the receptor. We show that η can be calculated from EC₅₀ and the asymptotes of a CRC derived from either single-channel or whole-cell responses. For 20 agonists of skeletal muscle nicotinic receptors, the distribution of η-values is bimodal with population means at 51% (including acetylcholine, nornicotine, and dimethylphenylpiperazinium) and 40% (including epibatidine, varenicline, and cytisine). The value of η is related inversely to the size of the agonist's headgroup, with high- versus low-efficiency ligands having an average volume of 70 vs. 102 Å³. Most binding site mutations have only a small effect on acetylcholine efficiency, except for αY190A (35%), αW149A (60%), and those at αG153 (42%). If η is known, the EC₅₀ and high-concentration asymptote can be calculated from each other. Hence, an entire CRC can be estimated from the response to a single agonist concentration, and efficacy can be estimated from EC₅₀ of a CRC that has been normalized to 1. Given η, the level of constitutive activity can be estimated from a single CRC.

A critical residue in the α1M2-M3 linker regulating mammalian GABAA receptor pore gating by diazepam

Benzodiazepines (BZDs) are a class of widely prescribed psychotropic drugs that modulate activity of GABA_A receptors (GABA_ARs), neurotransmitter-gated ion channels critical for synaptic transmission. However, the physical basis of this modulation is poorly understood. We explore the role of an important gating domain, the α₁M2–M3 linker, in linkage between the BZD site and pore gate. To probe energetics of this coupling without complication from bound agonist, we use a gain of function mutant (α₁L9'Tβ₂γ_2L) directly activated by BZDs. We identify a specific residue whose mutation (α₁V279A) more than doubles the energetic contribution of the BZD positive modulator diazepam (DZ) to pore opening and also enhances DZ potentiation of GABA-evoked currents in a wild-type background. In contrast, other linker mutations have little effect on DZ efficiency, but generally impair unliganded pore opening. Our observations reveal an important residue regulating BZD-pore linkage, thereby shedding new light on the molecular mechanism of these drugs.

Loss of Choline Agonism in the Inner Ear Hair Cell Nicotinic Acetylcholine Receptor Linked to the α10 Subunit

The α9α10 nicotinic acetylcholine receptor (nAChR) plays a fundamental role in inner ear physiology. It mediates synaptic transmission between efferent olivocochlear fibers that descend from the brainstem and hair cells of the auditory sensory epithelium. The α9 and α10 subunits have undergone a distinct evolutionary history within the family of nAChRs. Predominantly in mammalian vertebrates, the α9α10 receptor has accumulated changes at the protein level that may ultimately relate to the evolutionary history of the mammalian hearing organ. In the present work, we investigated the responses of α9α10 nAChRs to choline, the metabolite of acetylcholine degradation at the synaptic cleft. Whereas choline is a full agonist of chicken α9α10 receptors it is a partial agonist of the rat receptor. Making use of the expression of α9α10 heterologous receptors, encompassing wild-type, heteromeric, homomeric, mutant, chimeric, and hybrid receptors, and in silico molecular docking, we establish that the mammalian (rat) α10 nAChR subunit underscores the reduced efficacy of choline. Moreover, we show that whereas the complementary face of the α10 subunit does not play an important role in the activation of the receptor by ACh, it is strictly required for choline responses. Thus, we propose that the evolutionary changes acquired in the mammalian α9α10 nAChR resulted in the loss of choline acting as a full agonist at the efferent synapse, without affecting the triggering of ACh responses. This may have accompanied the fine-tuning of hair cell post-synaptic responses to the high-frequency activity of efferent medial olivocochlear fibers that modulate the cochlear amplifier.

Structure of the Native Muscle-type Nicotinic Receptor and Inhibition by Snake Venom Toxins

The nicotinic acetylcholine receptor, a pentameric ligand-gated ion channel, converts the free energy of binding of the neurotransmitter acetylcholine into opening of its central pore. Here we present the first high-resolution structure of the receptor type found in muscle-endplate membrane and in the muscle-derived electric tissues of fish. The native receptor was purified from Torpedo electric tissue and functionally reconstituted in lipids optimal for cryo-electron microscopy. The receptor was stabilized in a closed state by the binding of α-bungarotoxin. The structure reveals the binding of a toxin molecule at each of two subunit interfaces in a manner that would block the binding of acetylcholine. It also reveals a closed gate in the ion-conducting pore, formed by hydrophobic amino acid side chains, located ∼60 Å from the toxin binding sites. The structure provides a framework for understanding gating in ligand-gated channels and how mutations in the acetylcholine receptor cause congenital myasthenic syndromes.

Simulating ion channel activation mechanisms using swarms of trajectories

Atomic-level studies of protein activity represent a significant challenge as a result of the complexity of conformational changes occurring on wide-ranging timescales, often greatly exceeding that of even the longest simulations. A prime example is the elucidation of protein allosteric mechanisms, where localized perturbations transmit throughout a large macromolecule to generate a response signal. For example, the conversion of chemical to electrical signals during synaptic neurotransmission in the brain is achieved by specialized membrane proteins called pentameric ligand-gated ion channels. Here, the binding of a neurotransmitter results in a global conformational change to open an ion-conducting pore across the nerve cell membrane. X-ray crystallography has produced static structures of the open and closed states of the proton-gated GLIC pentameric ligand-gated ion channel protein, allowing for atomistic simulations that can uncover changes related to activation. We discuss a range of enhanced sampling approaches that could be used to explore activation mechanisms. In particular, we describe recent application of an atomistic string method, based on Roux's "swarms of trajectories" approach, to elucidate the sequence and interdependence of conformational changes during activation. We illustrate how this can be combined with transition analysis and Brownian dynamics to extract thermodynamic and kinetic information, leading to understanding of what controls ion channel function. © 2019 Wiley Periodicals, Inc.

Molecular dynamics simulations of dihydro-β-erythroidine bound to the human α4β2 nicotinic acetylcholine receptor

BACKGROUND AND PURPOSE

The heteromeric α4β2 nicotinic acetylcholine receptor (nAChR) is abundant in the human brain and is associated with a range of CNS disorders. This nAChR subtype has been recently crystallised in a conformation that was proposed to represent a desensitised state. Here, we investigated the conformational transition mechanism of this nAChR from a desensitised to a closed/resting state.

EXPERIMENTAL APPROACH

The competitive antagonist dihydro-β-erythroidine (DHβE) was modelled by replacement of the agonist nicotine in the α4β2 nAChR experimental structure. DHβE is used both in vitro and in vivo for its ability to block α4β2 nAChRs. This system was studied by three molecular dynamics simulations with a combined simulation time of 2.6 μs. Electrophysiological studies of mutated receptors were performed to validate the simulation results.

KEY RESULTS

The relative positions of the extracellular and transmembrane domains in the models are distinct from those of the desensitised state structure and are compatible with experimental structures of Cys-loop receptors captured in a closed/resting state.

CONCLUSIONS AND IMPLICATIONS

Our model suggests that the side chains of α4 L257 (9') and α4 L264 (16') are the main constrictions in the transmembrane pore. The involvement of position 9' in channel gating is well established, but position 16' was only previously identified as a gate for the bacterial channels, ELIC and GLIC. L257 but not L264 was found to influence the slow component of desensitisation. The structure of the antagonist-bound state proposed here should be valuable for the development of therapeutic or insecticide compounds.

A structural perspective on GABAA receptor pharmacology

GABA_A receptors are pentameric ligand-gated chloride channels of crucial importance for the vertebrate nervous system physiology. They typically modulate the fast inhibitory neurotransmission, and represent the target receptors for major classes of drugs used in the clinic, such as benzodiazepines and general anesthetics. Recent technological progress in structural biology, in particular single-particle cryo-electron microscopy, has led to fundamental advances in understanding the detailed organization and signalling mechanisms of major GABA_A receptor subtypes. This effort culminated with the high-resolution structural analysis of an intact, full-length human heteropentameric receptor, α1β3γ2, in a lipid bilayer and in complex with small molecule ligands including the commonly used benzodiazepines diazepam (Valium) and alprazolam (Xanax). These structures reveal multiple aspects of receptor activation and provide a path for rational design of subunit-specific GABA_A receptor modulators.

A single molecular distance predicts agonist binding energy in nicotinic receptors

Agonists turn on receptors because they bind more strongly to active (R*) versus resting (R) conformations of their target sites. Here, to explore how agonists activate neuromuscular acetylcholine receptors, we built homology models of R and R* neurotransmitter binding sites, docked ligands to those sites, ran molecular dynamics simulations to relax ("equilibrate") the structures, measured binding site structural parameters, and correlated them with experimental agonist binding energies. Each binding pocket is a pyramid formed by five aromatic amino acids and covered partially by loop C. We found that in R* versus R, loop C is displaced outward, the pocket is smaller and skewed, the agonist orientation is reversed, and a key nitrogen atom in the agonist is closer to the pocket center (distance dx) and a tryptophan pair but farther from αY190. Of these differences, the change in dx shows the largest correlation with experimental binding energy and provides a good estimate of agonist affinity, efficacy, and efficiency. Indeed, concentration-response curves can be calculated from just dx values. The contraction and twist of the binding pocket upon activation resemble gating rearrangements of the extracellular domain of related receptors at a smaller scale.

Efficiency measures the conversion of agonist binding energy into receptor conformational change

Receptors alternate between resting↔active conformations that bind agonists with low↔high affinity. Here, we define a new agonist attribute, energy efficiency (η), as the fraction of ligand-binding energy converted into the mechanical work of the activation conformational change. η depends only on the resting/active agonist-binding energy ratio. In a plot of activation energy versus binding energy (an "efficiency" plot), the slope gives η and the y intercept gives the receptor's intrinsic activation energy (without agonists; ΔG₀). We used single-channel electrophysiology to estimate η for eight different agonists and ΔG₀ in human endplate acetylcholine receptors (AChRs). From published equilibrium constants, we also estimated η for agonists of K_Ca1.1 (BK channels) and muscarinic, γ-aminobutyric acid, glutamate, glycine, and aryl-hydrocarbon receptors, and ΔG₀ for all of these except K_Ca1.1. Regarding AChRs, η is 48-56% for agonists related structurally to acetylcholine but is only ∼39% for agonists related to epibatidine; ΔG₀ is 8.4 kcal/mol in adult and 9.6 kcal/mol in fetal receptors. Efficiency plots for all of the above receptors are approximately linear, with η values between 12% and 57% and ΔG₀ values between 2 and 12 kcal/mol. Efficiency appears to be a general attribute of agonist action at receptor binding sites that is useful for understanding binding mechanisms, categorizing agonists, and estimating concentration-response relationships.

Electrostatics, proton sensor, and networks governing the gating transition in GLIC, a proton-gated pentameric ion channel

The pentameric ligand-gated ion channel (pLGIC) from Gloeobacter violaceus (GLIC) has provided insightful structure-function views on the permeation process and the allosteric regulation of the pLGICs family. However, GLIC is activated by pH instead of a neurotransmitter and a clear picture for the gating transition driven by protons is still lacking. We used an electrostatics-based (finite difference Poisson-Boltzmann/Debye-Hückel) method to predict the acidities of all aspartic and glutamic residues in GLIC, both in its active and closed-channel states. Those residues with a predicted pK_a close to the experimental pH₅₀ were individually replaced by alanine and the resulting variant receptors were titrated by ATR/FTIR spectroscopy. E35, located in front of loop F far away from the orthosteric site, appears as the key proton sensor with a measured individual pK_a at 5.8. In the GLIC open conformation, E35 is connected through a water-mediated hydrogen-bond network first to the highly conserved electrostatic triad R192-D122-D32 and then to Y197-Y119-K248, both located at the extracellular domain-transmembrane domain interface. The second triad controls a cluster of hydrophobic side chains from the M2-M3 loop that is remodeled during the gating transition. We solved 12 crystal structures of GLIC mutants, 6 of them being trapped in an agonist-bound but nonconductive conformation. Combined with previous data, this reveals two branches of a continuous network originating from E35 that reach, independently, the middle transmembrane region of two adjacent subunits. We conclude that GLIC's gating proceeds by making use of loop F, already known as an allosteric site in other pLGICs, instead of the classic orthosteric site.

The neural γ2α1β2α1β2 gamma amino butyric acid ion channel receptor: structural analysis of the effects of the ivermectin molecule and disulfide bridges

While ~30% of the human genome encodes membrane proteins, only a handful of structures of membrane proteins have been resolved to high resolution. Here, we studied the structure of a member of the Cys-loop ligand gated ion channel protein superfamily of receptors, human type A γ₂α₁β₂α₁β₂ gamma amino butyric acid receptor complex in a lipid bilayer environment. Studying the correlation between the structure and function of the gamma amino butyric acid receptor may enhance our understanding of the molecular basis of ion channel dysfunctions linked with epilepsy, ataxia, migraine, schizophrenia and other neurodegenerative diseases. The structure of human γ₂α₁β₂α₁β₂ has been modeled based on the X-ray structure of the Caenorhabditis elegans glutamate-gated chloride channel via homology modeling. The template provided the first inhibitory channel structure for the Cys-loop superfamily of ligand-gated ion channels. The only available template structure before this glutamate-gated chloride channel was a cation selective channel which had very low sequence identity with gamma aminobutyric acid receptor. Here, our aim was to study the effect of structural corrections originating from modeling on a more reliable template structure. The homology model was analyzed for structural properties via a 100 ns molecular dynamics (MD) study. Due to the structural shifts and the removal of an open channel potentiator molecule, ivermectin, from the template structure, helical packing changes were observed in the transmembrane segment. Namely removal of ivermectin molecule caused a closure around the Leu 9 position along the ion channel. In terms of the structural shifts, there are three potential disulfide bridges between the M1 and M3 helices of the γ₂ and 2 α₁ subunits in the model. The effect of these disulfide bridges was investigated via monitoring the differences in root mean square fluctuations (RMSF) of individual amino acids and principal component analysis of the MD trajectory of the two homology models-one with the disulfide bridge and one with protonated Cys residues. In all subunit types, RMSF of the transmembrane domain helices are reduced in the presence of disulfide bridges. Additionally, loop A, loop F and loop C fluctuations were affected in the extracellular domain. In cross-correlation analysis of the trajectory, the two model structures displayed different coupling in between the M2-M3 linker region, protruding from the membrane, and the β1-β2/D loop and cys-loop regions in the extracellular domain. Correlations of the C loop, which collapses directly over the bound ligand molecule, were also affected by differences in the packing of transmembrane helices. Finally, more localized correlations were observed in the transmembrane helices when disulfide bridges were present in the model. The differences observed in this study suggest that dynamic coupling at the interface of extracellular and ion channel domains differs from the coupling introduced by disulfide bridges in the transmembrane region. We hope that this hypothesis will be tested experimentally in the near future.

Nicotinic acetylcholine receptors at the single-channel level

Over the past four decades, the patch clamp technique and nicotinic ACh (nACh) receptors have established an enduring partnership. Like all good partnerships, each partner has proven significant in its own right, while their union has spurred innumerable advances in life science research. A member and prototype of the superfamily of pentameric ligand-gated ion channels, the nACh receptor is a chemo-electric transducer, binding ACh released from nerves and rapidly opening its channel to cation flow to elicit cellular excitation. A subject of a Nobel Prize in Physiology or Medicine, the patch clamp technique provides unprecedented resolution of currents through single ion channels in their native cellular environments. Here, focusing on muscle and α7 nACh receptors, we describe the extraordinary contribution of the patch clamp technique towards understanding how they activate in response to neurotransmitter, how subtle structural and mechanistic differences among nACh receptor subtypes translate into significant physiological differences, and how nACh receptors are being exploited as therapeutic drug targets.

LINKED ARTICLES

This article is part of a themed section on Nicotinic Acetylcholine Receptors. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v175.11/issuetoc/.

Spontaneous activity, singly bound states and the impact of alpha1Phe64 mutation on GABAAR gating in the novel kinetic model based on the single-channel recordings

GABA_A receptor is the primary mediator of inhibition in the adult mammalian brain. Our recent studies revealed that a classic gating scheme for GABA_AR needed to be updated with an intermediate step (flipping) and that the α₁Phe64 mutation at the GABA binding site affects this transition. However, description of flipping at the single-channel level remains incomplete. In particular, its role in singly-bound and spontaneous activity remains unknown. We have performed thus single-channel recordings over wide range of agonist concentration for wild-type α₁β₂γ_2L receptors and α₁Phe64 mutants. For WT receptors we observed relatively frequent brief spontaneous openings which were also present at low [GABA]. However, closed times distributions for spontaneous activity and at low [GABA] were clearly different indicating that a proportion of short-lived openings were due to liganded, most likely singly bound receptors. Increasing [GABA] resulted in prolongation of bursts and increased occurrence of bursts with long openings and short closures. Mutations of α₁Phe64 residue dramatically affected the open and closed time distributions at high and saturating [GABA], especially in the case of cysteine mutants. However, this mutation weakly affected spontaneous or singly bound activity. Model fitting of our single-channel data led us to propose a novel and, to our knowledge, most complete GABA_AR kinetic model in which flipping occurs in singly and doubly bound states. However, spontaneous activity did not reveal involvement of flipping. Moreover, we report that α₁Phe64 mutation affects not only the flipping but also the opening/closing transitions indicating its generalized impact on the receptor gating.

GABA Type A Receptor Activation in the Allosteric Coagonist Model Framework: Relationship between EC50 and Basal Activity

The concerted transition model for multimeric proteins is a simple formulation for analyzing the behavior of transmitter-gated ion channels. We used the model to examine the relationship between the EC₅₀ for activation of the GABA type A (GABA_A) receptor by the transmitter GABA and basal activity employing concatemeric ternary GABA_A receptors expressed in Xenopus oocytes. Basal activity, reflecting the receptor function in the absence of the transmitter, can be changed either by mutation to increase constitutive activity or by the addition of a second agonist (acting at a different site) to increase background activity. The model predicts that either mechanism for producing a change in basal activity will result in identical effects on the EC₅₀ We examined receptor activation by GABA while changing the level of basal activity with the allosterically acting anesthetics propofol, pentobarbital, or alfaxalone. We found that the relationship between EC₅₀ and basal activity was well described by the concerted transition model. Changes in the basal activity by gain-of-function mutations also resulted in predictable changes in the EC₅₀ Finally, we altered the number of GABA-binding sites by a mutation and again found that the relationship could be well described by the model. Overall, the results support the idea that interactions between the transmitter GABA and the allosteric agonists propofol, pentobarbital, or alfaxalone can be understood as reflecting additive and independent free energy changes, without assuming any specific interactions.

Un-gating and allosteric modulation of a pentameric ligand-gated ion channel captured by molecular dynamics

Pentameric ligand-gated ion channels (pLGICs) mediate intercellular communication at synapses through the opening of an ion pore in response to the binding of a neurotransmitter. Despite the increasing availability of high-resolution structures of pLGICs, a detailed understanding of the functional isomerization from closed to open (gating) and back is currently missing. Here, we provide the first atomistic description of the transition from open to closed (un-gating) in the glutamate-gated chloride channel (GluCl) from Caenorhabditis Elegans. Starting with the active-state structure solved in complex with the neurotransmitter L-glutamate and the positive allosteric modulator (PAM) ivermectin, we analyze the spontaneous relaxation of the channel upon removal of ivermectin by explicit solvent/membrane Molecular Dynamics (MD) simulations. The μs-long trajectories support the conclusion that ion-channel deactivation is mediated by two distinct quaternary transitions, i.e. a global receptor twisting followed by the radial expansion (or blooming) of the extracellular domain. At variance with previous models, we show that pore closing is exclusively regulated by the global twisting, which controls the position of the β1-β2 loop relative to the M2-M3 loop at the EC/TM domain interface. Additional simulations with L-glutamate restrained to the crystallographic binding mode and ivermectin removed indicate that the same twisting isomerization is regulated by agonist binding at the orthosteric site. These results provide a structural model for gating in pLGICs and suggest a plausible mechanism for the pharmacological action of PAMs in this neurotransmitter receptor family. The simulated un-gating converges to the X-ray structure of GluCl resting state both globally and locally, demonstrating the predictive character of state-of-art MD simulations.

A chimeric prokaryotic-eukaryotic pentameric ligand gated ion channel reveals interactions between the extracellular and transmembrane domains shape neurosteroid modulation

Pentameric ligand-gated ion channels (pLGICs) are the targets of several clinical and endogenous allosteric modulators including anesthetics and neurosteroids. Molecular mechanisms underlying allosteric drug modulation are poorly understood. Here, we constructed a chimeric pLGIC by fusing the extracellular domain (ECD) of the proton-activated, cation-selective bacterial channel GLIC to the transmembrane domain (TMD) of the human ρ1 chloride-selective GABAAR, and tested the hypothesis that drug actions are regulated locally in the domain that houses its binding site. The chimeric channels were proton-gated and chloride-selective demonstrating the GLIC ECD was functionally coupled to the GABAρ TMD. Channels were blocked by picrotoxin and inhibited by pentobarbital, etomidate and propofol. The point mutation, ρ TMD W328M, conferred positive modulation and direct gating by pentobarbital. The data suggest that the structural machinery mediating general anesthetic modulation resides in the TMD. Proton-activation and neurosteroid modulation of the GLIC-ρ chimeric channels, however, did not simply mimic their respective actions on GLIC and GABAρ revealing that across domain interactions between the ECD and TMD play important roles in determining their actions. Proton-induced current responses were biphasic suggesting that the chimeric channels contain an additional proton sensor. Neurosteroid modulation of the GLIC-ρ chimeric channels by the stereoisomers, 5α-THDOC and 5β-THDOC, were swapped compared to their actions on GABAρ indicating that positive versus negative neurosteroid modulation is not encoded solely in the TMD nor by neurosteroid isomer structure but is dependent on specific interdomain connections between the ECD and TMD. Our data reveal a new mechanism for shaping neurosteroid modulation.

Back to the future: Rational maps for exploring acetylcholine receptor space and time

Global functions of nicotinic acetylcholine receptors, such as subunit cooperativity and compatibility, likely emerge from a network of amino acid residues distributed across the entire pentameric complex. Identification of such networks has stymied traditional approaches to acetylcholine receptor structure and function, likely due to the cryptic interdependency of their underlying amino acid residues. An emerging evolutionary biochemistry approach, which traces the evolutionary history of acetylcholine receptor subunits, allows for rational mapping of acetylcholine receptor sequence space, and offers new hope for uncovering the amino acid origins of these enigmatic properties.

String method solution of the gating pathways for a pentameric ligand-gated ion channel

Pentameric ligand-gated ion channels control synaptic neurotransmission by converting chemical signals into electrical signals. Agonist binding leads to rapid signal transduction via an allosteric mechanism, where global protein conformational changes open a pore across the nerve cell membrane. We use all-atom molecular dynamics with a swarm-based string method to solve for the minimum free-energy gating pathways of the proton-activated bacterial GLIC channel. We describe stable wetted/open and dewetted/closed states, and uncover conformational changes in the agonist-binding extracellular domain, ion-conducting transmembrane domain, and gating interface that control communication between these domains. Transition analysis is used to compute free-energy surfaces that suggest allosteric pathways; stabilization with pH; and intermediates, including states that facilitate channel closing in the presence of an agonist. We describe a switching mechanism that senses proton binding by marked reorganization of subunit interface, altering the packing of β-sheets to induce changes that lead to asynchronous pore-lining M2 helix movements. These results provide molecular details of GLIC gating and insight into the allosteric mechanisms for the superfamily of pentameric ligand-gated channels.

Orthosteric and allosteric potentiation of heteromeric neuronal nicotinic acetylcholine receptors

Heteromeric nicotinic ACh receptors (nAChRs) were thought to have two orthodox agonist-binding sites at two α/β subunit interfaces. Highly selective ligands are hard to develop by targeting orthodox agonist sites because of high sequence similarity of this binding pocket among different subunits. Recently, unorthodox ACh-binding sites have been discovered at some α/α and β/α subunit interfaces, such as α4/α4, α5/α4 and β3/α4. Targeting unorthodox sites may yield subtype-selective ligands, such as those for (α4β2)₂ α5, (α4β2)₂ β3 and (α6β2)₂ β3 nAChRs. The unorthodox sites have unique pharmacology. Agonist binding at one unorthodox site is not sufficient to activate nAChRs, but it increases activation from the orthodox sites. NS9283, a selective agonist for the unorthodox α4/α4 site, was initially thought to be a positive allosteric modulator (PAM). NS9283 activates nAChRs with three engineered α4/α4 sites. PAMs, on the other hand, act at allosteric sites where ACh cannot bind. Known PAM sites include the ACh-homologous non-canonical site (e.g. morantel at β/α), the C-terminus (e.g. Br-PBTC and 17β-estradiol), a transmembrane domain (e.g. LY2087101) or extracellular and transmembrane domain interfaces (e.g. NS206). Some of these PAMs, such as Br-PBTC and 17β-estradiol, require only one subunit to potentiate activation of nAChRs. In this review, we will discuss differences between activation from orthosteric and allosteric sites, their selective ligands and clinical implications. These studies have advanced understanding of the structure, assembly and pharmacology of heteromeric neuronal nAChRs.

LINKED ARTICLES

This article is part of a themed section on Nicotinic Acetylcholine Receptors. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v175.11/issuetoc.