Temperature dependence of acetylcholine receptor channels activated by different agonists

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.

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.

Incomplete removal of extracellular glutamate controls synaptic transmission and integration at a cerebellar synapse

Synapses of glutamatergic mossy fibers (MFs) onto cerebellar unipolar brush cells (UBCs) generate slow excitatory (ON) or inhibitory (OFF) postsynaptic responses dependent on the complement of glutamate receptors expressed on the UBC's large dendritic brush. Using mouse brain slice recording and computational modeling of synaptic transmission, we found that substantial glutamate is maintained in the UBC synaptic cleft, sufficient to modify spontaneous firing in OFF UBCs and tonically desensitize AMPARs of ON UBCs. The source of this ambient glutamate was spontaneous, spike-independent exocytosis from the MF terminal, and its level was dependent on activity of glutamate transporters EAAT1-2. Increasing levels of ambient glutamate shifted the polarity of evoked synaptic responses in ON UBCs and altered the phase of responses to in vivo-like synaptic activity. Unlike classical fast synapses, receptors at the UBC synapse are virtually always exposed to a significant level of glutamate, which varies in a graded manner during transmission.

Effects of hyperthermia on the effective concentration of rocuronium and sugammadex-mediated reversal in isolated phrenic nerve hemidiaphragm preparations of rats

Background

Hyperthermia is relatively rare during general anesthesia; however, a few studies have been conducted on hyperthermia and the neuromuscular blockade (NMB) induced by rocuronium, and the reversal of NMB by sugammadex. We investigated the effect of hyperthermia status on the NMB induced by rocuronium, and its reversal by sugammadex, in isolated phrenic nerve hemidiaphragm (PNHD) preparations of the rat.

Methods

Thirty-three male Sprague-Dawley rat PNHD preparations were randomly assigned to three groups at different temperatures (36 °C, 38 °C, and 40 °C; each group, n = 11, in Krebs solution). The train-of-four (TOF) and twitch height responses were checked mechanomyographically. The PNHD were treated with progressively increasing doses of rocuronium and three effective concentrations (ECs), EC50, EC90, and EC95, of rocuronium were analyzed in each group via nonlinear regression analysis. Then, sugammadex was administered in doses equimolar to rocuronium. Thereafter, the T1 height (%), TOFR (%) and the duration index were measured.

Results

The EC of rocuronium (EC50, EC90, and EC95) decreased significantly in accordance with increasing temperature. The groups at 36 °C and 40 °C showed clear differences in all areas (all P < 0.001). Moreover, the T1 height (%) and the duration index upon sugammadex administration showed faster recovery results in the36 °C than the 38 °C and 40 °C groups.

Conclusion

A rise of temperature from 38 °C to 40 °C in rat PNHD preparations proportionally enhanced the NMB induced by rocuronium. In addition, equimolar doses of sugammadex to the administered rocuronium showed a slower recovery time as the temperature rises.

Hyperthermia exacerbates the acute effects of psychoactive substances on neuronal activity measured using microelectrode arrays (MEAs) in rat primary cortical cultures in vitro

Hyperthermia is a well-known, potentially life-threatening, side effect of stimulant psychoactive substances that worsens the neurological outcome of hospitalized patients. However, current in vitro methods to assess the hazard of psychoactive substances do not account for hyperthermia. Therefore, this study determined the potency of five psychoactive substances (cocaine, MDMA (3,4-methylenedioxymethamphetamine), methamphetamine, 3-MMC (3-methylmethcathinone) and TFMPP (3-trifluoromethylphenylpiperazine)) to affect neuronal activity at physiological and hyperthermic conditions. Neuronal activity of rat cortical cultures grown on microelectrode arrays (MEAs) was recorded at 37 °C before exposure. Following 30 min and 4.5 h drug exposure (1-1000 μM) at 37 °C or 41 °C, neuronal activity was measured at either 37 °C or 41 °C. Without drug exposure, hyperthermia induced a modest decrease in neuronal activity. Following acute (30 min) exposure at 37 °C, all drugs concentration-dependently inhibited neuronal activity. Increasing the temperature to 41 °C significantly exacerbated the reduction of neuronal activity ~ 2-fold for all drugs compared to 37 °C. Prolonged (4.5 h) exposure at 41 °C decreased neuronal activity comparable to 37 °C. Neuronal activity (partly) recovered following drug exposure at both temperatures, although recovery from exposure at 41 °C was less pronounced for most drugs. None of the exposure conditions affected viability. Since acute exposure at hyperthermic conditions exacerbates the decrease in neuronal activity induced by psychoactive substances, effects of hyperthermia should be included in future hazard assessment of illicit drugs and new psychoactive substances (NPS).

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.

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/.

Cyclic activation of endplate acetylcholine receptors

Agonists turn on receptors because they have a higher affinity for active versus resting conformations of the protein. Activation can occur by either of two pathways that connect to form a cycle: Agonists bind to resting receptors that then become active, or resting receptors activate and then bind agonists. We used mutations to construct endplate acetylcholine receptors (AChRs) having only one functional neurotransmitter-binding site and single-channel electrophysiology to measure independently binding constants for four different agonists, to both resting and active conformations of each site. For all agonists and sites, the total free energy change in each pathway was the same, confirming the activation cycle without external energy. Other results show that (i) there is no cooperativity between sites; (ii) agonist association is slower than diffusion in resting receptors but nearly diffusional in active receptors; (iii) whereas resting affinity is determined mainly by agonist association, active affinity is determined mainly by agonist dissociation; and (iv) at each site and for all agonists, receptor activation approximately doubles the agonist-binding free energy. We discuss a two-step mechanism for binding that involves diffusion and a local conformational change ("catch") that is modulated by receptor activation. The results suggest that binding to a resting site and the switch to high affinity are both integral parts of a single allosteric transition. We hypothesize that catch ensures proper signal recognition in complex chemical environments and that binding site compaction is a determinant of both resting and active affinity.

The Kinetics and the Permeation Properties of Piezo Channels

Piezo channels are eukaryotic, cation-selective mechanosensitive channels (MSCs), which show rapid activation and voltage-dependent inactivation. The kinetics of these channels are largely consistent across multiple cell types and different stimulation paradigms with some minor variability. No accessory subunits that associate with Piezo channels have been reported. They are homotrimers and each ∼300kD monomer has an N-terminal propeller blade-like mechanosensing module, which can confer mechanosensing capabilities on ASIC-1 (the trimeric non-MSC, acid-sensing ion channel-1) and a C-terminal pore module, which influences conductance, selectivity, and channel inactivation. Repeated stimulation can cause domain fracture and diffusion of these channels leading to synchronous loss of inactivation. The reconstituted channels spontaneously open only in asymmetric bilayers but lack inactivation. Mutations that cause hereditary xerocytosis alter PIEZO1 kinetics. The kinetics of the wild-type PIEZO1 and alterations thereof in mutants (M2225R, R2456K, and DhPIEZO1) are summarized in the form of a quantitative model and hosted online. The pore is permeable to alkali ions although Li⁺ permeates poorly. Divalent cations, notably Ca²⁺, traverse the channel and inhibit the flux of monovalents. The large monovalent organic cations such as tetramethyl ammonium and tetraethyl ammonium can traverse the channel, but slowly, suggesting a pore diameter of ∼8Å, and the estimated in-plane area change upon opening is around 6-20nm². Ruthenium red can enter the channel only from the extracellular side and seems to bind in a pocket close to residue 2496.

Structural correlates of affinity in fetal versus adult endplate nicotinic receptors

Adult-type nicotinic acetylcholine receptors (AChRs) mediate signalling at mature neuromuscular junctions and fetal-type AChRs are necessary for proper synapse development. Each AChR has two neurotransmitter binding sites located at the interface of a principal and a complementary subunit. Although all agonist binding sites have the same core of five aromatic amino acids, the fetal site has ∼30-fold higher affinity for the neurotransmitter ACh. Here we use molecular dynamics simulations of adult versus fetal homology models to identify complementary-subunit residues near the core that influence affinity, and use single-channel electrophysiology to corroborate the results. Four residues in combination determine adult versus fetal affinity. Simulations suggest that at lower-affinity sites, one of these unsettles the core directly and the others (in loop E) increase backbone flexibility to unlock a key, complementary tryptophan from the core. Swapping only four amino acids is necessary and sufficient to exchange function between adult and fetal AChRs.

Adult and fetal nicotinic acetylcholine receptors (AChRs) have different functional requirements and affinity for ACh. Here, the authors use molecular dynamics and electrophysiology to investigate this affinity, and identify four amino acids that when swapped exchange function between adult and fetal AChRs.

Exploring the positive allosteric modulation of human α7 nicotinic receptors from a single-channel perspective

Enhancement of α7 nicotinic receptor (nAChR) function by positive allosteric modulators (PAMs) is a promising therapeutic strategy to improve cognitive deficits. PAMs have been classified only on the basis of their macroscopic effects as type I, which only enhance agonist-induced currents, and type II, which also decrease desensitization and reactivate desensitized nAChRs. To decipher the molecular basis underlying these distinct activities, we explored the effects on single-α7 channel currents of representative members of each type and of less characterized compounds. Our results reveal that all PAMs enhance open-channel lifetime and produce episodes of successive openings, thus indicating that both types affect α7 kinetics. Different PAM types show different sensitivity to temperature, suggesting different mechanisms of potentiation. By using a mutant α7 receptor that is insensitive to the prototype type II PAM (PNU-120596), we show that some though not all type I PAMs share the structural determinants of potentiation. Overall, our study provides novel information on α7 potentiation, which is key to the ongoing development of therapeutic compounds.

Agonist activation of a nicotinic acetylcholine receptor

How does an agonist activate a receptor? In this article I consider the activation process in muscle nicotinic acetylcholine receptors (AChRs), a prototype for understanding the energetics of binding and gating in other ligand-gated ion channels. Just as movements that generate gating currents activate voltage-gated ion channels, movements at binding sites that generate an increase in affinity for the agonist activate ligand-gated ion channels. The main topics are: i) the schemes and intermediate states of AChR activation, ii) the energy changes of each of the steps, iii) the sources of the energies, iv) the three kinds of AChR agonist binding site and v) the correlations between binding and gating energies. The binding process is summarized as sketches of different conformations of an agonist site. The results suggest that agonists lower the free energy of the active conformation of the protein in stages by establishing favorable, local interactions at each binding site, independently. This article is part of the Special Issue entitled 'The Nicotinic Acetylcholine Receptor: From Molecular Biology to Cognition'.

Activation of endplate nicotinic acetylcholine receptors by agonists

The interaction of a small molecule made in one cell with a large receptor made in another is the signature event of cell signaling. Understanding the structure and energy changes associated with agonist activation is important for engineering drugs, receptors and synapses. The nicotinic acetylcholine receptor (AChR) is a ∼300kD ion channel that binds the neurotransmitter acetylcholine (ACh) and other cholinergic agonists to elicit electrical responses in the central and peripheral nervous systems. This mini-review is in two sections. First, general concepts of skeletal muscle AChR operation are discussed in terms of energy landscapes for conformational change. Second, adult vs. fetal AChRs are compared with regard to interaction energies between ACh and agonist-site side chains, measured by single-channel electrophysiology and molecular dynamics simulations. The five aromatic residues that form the core of each agonist binding site can be divided into two working groups, a triad (led by αY190) that behaves similarly at all sites and a coupled pair (led by γW55) that has a large influence on affinity only in fetal AChRs. Each endplate AChR has 5 homologous subunits, two of α(1) and one each of β, δ, and either γ (fetal) or ϵ (adult). These nicotinic AChRs have only 2 functional agonist binding sites located in the extracellular domain, at αδ and either αγ or αϵ subunit interfaces. The receptor undergoes a reversible, global isomerization between structures called C and O. The C shape does not conduct ions and has a relatively low affinity for ACh, whereas O conducts cations and has a higher affinity. When both agonist sites are empty (filled only with water) the probability of taking on the O conformation (PO) is low, <10(-6). When ACh molecules occupy the agonist sites the C→O opening rate constant and C↔O gating equilibrium constant increase dramatically. Following a pulse of ACh at the nerve-muscle synapse, the endplate current rises rapidly to reach a peak that corresponds to PO ∼0.96.

Quasithermodynamic contributions to the fluctuations of a protein nanopore

Proteins undergo thermally activated conformational fluctuations among two or more substates, but a quantitative inquiry on their kinetics is persistently challenged by numerous factors, including the complexity and dynamics of various interactions, along with the inability to detect functional substates within a resolvable time scale. Here, we analyzed in detail the current fluctuations of a monomeric β-barrel protein nanopore of known high-resolution X-ray crystal structure. We demonstrated that targeted perturbations of the protein nanopore system, in the form of loop-deletion mutagenesis, accompanying alterations of electrostatic interactions between long extracellular loops, produced modest changes of the differential activation free energies calculated at 25 °C, ΔΔG(⧧), in the range near the thermal energy but substantial and correlated modifications of the differential activation enthalpies, ΔΔH(⧧), and entropies, ΔΔS(⧧). This finding indicates that the local conformational reorganizations of the packing and flexibility of the fluctuating loops lining the central constriction of this protein nanopore were supplemented by changes in the single-channel kinetics. These changes were reflected in the enthalpy-entropy reconversions of the interactions between the loop partners with a compensating temperature, TC, of ∼300 K, and an activation free energy constant of ∼41 kJ/mol. We also determined that temperature has a much greater effect on the energetics of the equilibrium gating fluctuations of a protein nanopore than other environmental parameters, such as the ionic strength of the aqueous phase as well as the applied transmembrane potential, likely due to ample changes in the solvation activation enthalpies. There is no fundamental limitation for applying this approach to other complex, multistate membrane protein systems. Therefore, this methodology has major implications in the area of membrane protein design and dynamics, primarily by revealing a better quantitative assessment on the equilibrium transitions among multiple well-defined and functionally distinct substates of protein channels and pores.

Functional differences between neurotransmitter binding sites of muscle acetylcholine receptors

A muscle acetylcholine receptor (AChR) has two neurotransmitter binding sites located in the extracellular domain, at αδ and either αε (adult) or αγ (fetal) subunit interfaces. We used single-channel electrophysiology to measure the effects of mutations of five conserved aromatic residues at each site with regard to their contribution to the difference in free energy of agonist binding to active versus resting receptors (ΔGB1). The two binding sites behave independently in both adult and fetal AChRs. For four different agonists, including ACh and choline, ΔGB1 is ∼-2 kcal/mol more favorable at αγ compared with at αε and αδ. Only three of the aromatics contribute significantly to ΔGB1 at the adult sites (αY190, αY198, and αW149), but all five do so at αγ (as well as αY93 and γW55). γW55 makes a particularly large contribution only at αγ that is coupled energetically to those contributions of some of the α-subunit aromatics. The hydroxyl and benzene groups of loop C residues αY190 and αY198 behave similarly with regard to ΔGB1 at all three kinds of site. ACh binding energies estimated from molecular dynamics simulations are consistent with experimental values from electrophysiology and suggest that the αγ site is more compact, better organized, and less dynamic than αε and αδ. We speculate that the different sensitivities of the fetal αγ site versus the adult αε and αδ sites to choline and ACh are important for the proper maturation and function of the neuromuscular synapse.

Functional anatomy of an allosteric protein

Synaptic receptors are allosteric proteins that switch on and off to regulate cell signalling. Here, we use single-channel electrophysiology to measure and map energy changes in the gating conformational change of a nicotinic acetylcholine receptor. Two separated regions in the α-subunits--the transmitter-binding sites and αM2-αM3 linkers in the membrane domain--have the highest ϕ-values (change conformation the earliest), followed by the extracellular domain, most of the membrane domain and the gate. Large gating-energy changes occur at the transmitter-binding sites, α-subunit interfaces, the αM1 helix and the gate. We hypothesize that rearrangements of the linkers trigger the global allosteric transition, and that the hydrophobic gate unlocks in three steps. The mostly local character of side-chain energy changes and the similarly high ϕ-values of separated domains, both with and without ligands, suggest that gating is not strictly a mechanical process initiated by the affinity change for the agonist.

Catch-and-hold activation of muscle acetylcholine receptors having transmitter binding site mutations

Agonists turn on receptors because their target sites have a higher affinity in the active versus resting conformation of the protein. We used single-channel electrophysiology to measure the lower-affinity (LA) and higher-affinity (HA) equilibrium dissociation constants for acetylcholine in adult-type muscle mouse nicotinic receptors (AChRs) having mutations of agonist binding site amino acids. For a series of agonists and for all mutations of αY93, αG147, αW149, αY190, αY198, εW55, and δW57, the change in LA binding energy was approximately half that in HA binding energy. The results were analyzed as a linear free energy relationship between LA and HA agonist binding, the slope of which (κ) gives the fraction of the overall binding chemical potential where the LA complex is established. The linear correlation between LA and HA binding energies suggests that the overall binding process is by an integrated mechanism (catch-and-hold). For the agonist and the above mutations, κ ∼ 0.5, but side-chain substitutions of two residues had a slope that was significantly higher (0.90; αG153) or lower (0.25; εP121). The results suggest that backbone rearrangements in loop B, loop C, and the non-α surface participate in both LA binding and the LA ↔ HA affinity switch. It appears that all of the intermediate steps in AChR activation comprise a single, energetically coupled process.

SOLVING ION CHANNEL KINETICS WITH THE QuB SOFTWARE

The ability to record the currents from single ion channels led to the need to extract the underlying kinetic model from such data. This inverse hidden Markov problem is difficult but led to the creation of a software suite called QuB utilizing likelihood optimization. This review presents the software. The software is open source and, in addition to solving kinetic models, has many generic database operations including report generation with publishable graphics, function fitting and scripting for new and repeated processing and AD/DA I/O. The core algorithms allow for constraints such as fixed rates or maintaining detailed balance in the model. All rate constants can be driven by a stimulus and the system can analyze nonstationary data. QuB also can analyze the kinetics of multichannel data where individual events cannot be discriminated, but the fitting algorithms utilize the signal variance as well as the mean to fit models. QuB can be applied to any data appropriately modeled with Markov kinetics and has been utilized to solve ion channels but also the movement of motor proteins, the sleep cycles in mice, and physics processes.

[Formula: see text]Special Issue Comment: This is a review about the software QuB that can extract a model from the trajectory. It is connected with the review about treatments when solving single molecules,60 and the reviews about enzymes.61,62

Asymmetric transmitter binding sites of fetal muscle acetylcholine receptors shape their synaptic response

Neuromuscular acetylcholine receptors (AChRs) have two transmitter binding sites: at α-δ and either α-γ (fetal) or α-ε (adult) subunit interfaces. The γ-subunit of fetal AChRs is indispensable for the proper development of neuromuscular synapses. We estimated parameters for acetylcholine (ACh) binding and gating from single channel currents of fetal mouse AChRs expressed in tissue-cultured cells. The unliganded gating equilibrium constant is smaller and less voltage-dependent than in adult AChRs. However, the α-γ binding site has a higher affinity for ACh and provides more binding energy for gating compared with α-ε; therefore, the diliganded gating equilibrium constant at -100 mV is comparable for both receptor subtypes. The -2.2 kcal/mol extra binding energy from α-γ compared with α-δ and α-ε is accompanied by a higher resting affinity for ACh, mainly because of slower transmitter dissociation. End plate current simulations suggest that the higher affinity and increased energy from α-γ are essential for generating synaptic responses at low pulse [ACh].

Action of nicotine and analogs on acetylcholine receptors having mutations of transmitter-binding site residue αG153

A primary target for nicotine is the acetylcholine receptor channel (AChR). Some of the ability of nicotine to activate differentially AChR subtypes has been traced to a transmitter-binding site amino acid that is glycine in lower affinity and lysine in higher affinity AChRs. We studied the effects of mutations of this residue (αG153) in neuromuscular AChRs activated by nicotine and eight other agonists including nornicotine and anabasine. All of the mutations increased the unliganded gating equilibrium constant. The affinity of the resting receptor (K_d) and the net binding energy from the agonist for gating (ΔG_B) were estimated by cross-concentration fitting of single-channel currents. In all but one of the agonist/mutant combinations there was a moderate decrease in K_d and essentially no change in ΔG_B. The exceptional case was nicotine plus lysine, which showed a large, >8,000-fold decrease in K_d but no change in ΔG_B. The extraordinary specificity of this combination leads us to speculate that AChRs with a lysine at position αG153 may be exposed to a nicotine-like compound in vivo.

Loop C and the mechanism of acetylcholine receptor-channel gating

Agonist molecules at the two neuromuscular acetylcholine (ACh) receptor (AChR) transmitter-binding sites increase the probability of channel opening. In one hypothesis for AChR activation (“priming”), the capping of loop C at each binding site transfers energy independently to the distant gate over a discrete structural pathway. We used single-channel analyses to examine the experimental support for this proposal with regard to brief unliganded openings, the effects of loop-C modifications, the effects of mutations to residues either on or off the putative pathway, and state models for describing currents at low [ACh]. The results show that (a) diliganded and brief unliganded openings are generated by the same essential, global transition; (b) the radical manipulation of loop C does not prevent channel opening but impairs agonist binding; (c) both on- and off-pathway mutations alter gating by changing the relative stability of the open-channel conformation by local interactions rather than by perturbing a specific site–gate communication link; and (d) it is possible to estimate directly the rate constants for agonist dissociation from and association to both the low and high affinity forms of the AChR-binding site by using a cyclic kinetic model. We conclude that the mechanism of energy transfer between the binding sites and the gate remains an open question.

The energy and work of a ligand-gated ion channel

Ligand-gated ion channels are allosteric membrane proteins that isomerize between C(losed) and O(pen) conformations. A difference in affinity for ligands in the two states influences the C↔O "gating" equilibrium constant. The energies associated with adult-type mouse neuromuscular nicotinic acetylcholine receptor (AChR) channel gating have been measured by using single-channel electrophysiology. Without ligands, the free energy, enthalpy and entropy of gating are ΔG0=+8.4, ΔH0=+10.9 and TΔS0=+2.5kcal/mol (-100mV, 23°C). Many mutations throughout the protein change ΔG0, including natural ones that cause disease. Agonists and most mutations change approximately independently the ground-state energy difference; thus, it is possible to forecast and engineer AChR responses simply by combining perturbations. The free energy of the low↔high affinity change for the neurotransmitter at each of two functionally equivalent binding sites is ΔGB(ACh)=-5.1kcal/mol. ΔGB(ACh) is set mainly by interactions of ACh with just three binding site aromatic groups. For a series of structurally related agonists, there is a correlation between the energies of low- and high-affinity binding, which implies that gating commences with the formation of the low-affinity complex. Brief, intermediate states in binding and gating have been detected. Several proposals for the nature of the gating transition-state energy landscape and the isomerization mechanism are discussed.

Structure and pharmacology of pentameric receptor channels: from bacteria to brain

Orthologs of the pentameric receptor channels that mediate fast synaptic transmission in the central and peripheral nervous systems have been found in several bacterial species and in a single archaea genus. Recent X-ray structures of bacterial and invertebrate pentameric receptors point to a striking conservation of the structural features within the whole family, even between distant prokaryotic and eukaryotic members. These structural data reveal general principles of molecular organization that allow allosteric membrane proteins to mediate chemoelectric transduction. Notably, several conformations have been solved, including open and closed channels with distinct global tertiary and quaternary structure. The data reveal features of the ion channel architecture and of diverse categories of binding sites, such as those that bind orthosteric ligands, including neurotransmitters, and those that bind allosteric modulators, such as general anesthetics, ivermectin, or lipids. In this review, we summarize the most recent data, discuss insights into the mechanism of action in these systems, and elaborate on newly opened avenues for drug design.

An integrated catch-and-hold mechanism activates nicotinic acetylcholine receptors

In neuromuscular acetylcholine (ACh) receptor channels (AChRs), agonist molecules bind with a low affinity (LA) to two sites that can switch to high affinity (HA) and increase the probability of channel opening. We measured (by using single-channel kinetic analysis) the rate and equilibrium constants for LA binding and channel gating for several different agonists of adult-type mouse AChRs. Almost all of the variation in the equilibrium constants for LA binding was from differences in the association rate constants. These were consistently below the limit set by diffusion and were substantially different even though the agonists had similar sizes and the same charge. This suggests that binding to resting receptors is not by diffusion alone and, hence, that each binding site can undergo two conformational changes ("catch" and "hold") that connect three different structures (apo-, LA-bound, and HA-bound). Analyses of ACh-binding protein structures suggest that this binding site, too, may adopt three discrete structures having different degrees of loop C displacement ("capping"). For the agonists we tested, the logarithms of the equilibrium constants for LA binding and LA↔HA gating were correlated. Although agonist binding and channel gating have long been considered to be separate processes in the activation of ligand-gated ion channels, this correlation implies that the catch-and-hold conformational changes are energetically linked and together comprise an integrated process having a common structural basis. We propose that loop C capping mainly reflects agonist binding, with its two stages corresponding to the formation of the LA and HA complexes. The catch-and-hold reaction coordinate is discussed in terms of preopening states and thermodynamic cycles of activation.

Entropy-enthalpy Compensation of Biomolecular Systems in Aqueous Phase: a Dry Perspective

We survey thermodynamic measurements on processes involving biological macromolecules in aqueous solution, which illustrate well the ubiquitous phenomenon of entropy-enthalpy compensation. The processes include protein folding/unfolding and ligand binding/unbinding, with compensation temperatures varying by about 50 K around an average near 293 K. We show that incorporating both near-exact entropy-enthalpy compensation (due to solvent relaxation) and multi-excitation entropy (from vibrational quanta) leads to a compensation temperature in water of about 230 K. We illustrate a general procedure for subtracting solvent and environment-related terms to determine the bare Gibbs free energy changes of chemical processes.

The intrinsic energy of the gating isomerization of a neuromuscular acetylcholine receptor channel

Nicotinic acetylcholine receptor (AChR) channels at neuromuscular synapses rarely open in the absence of agonists, but many different mutations increase the unliganded gating equilibrium constant (E₀) to generate AChRs that are active constitutively. We measured E₀ for two different sets of mutant combinations and by extrapolation estimated E₀ for wild-type AChRs. The estimates were 7.6 and 7.8 × 10⁻⁷ in adult-type mouse AChRs (−100 mV at 23°C). The values are in excellent agreement with one obtained previously by using a completely different method (6.5 × 10⁻⁷, from monoliganded gating). E₀ decreases with depolarization to the same extent as does the diliganded gating equilibrium constant, e-fold with ∼60 mV. We estimate that at −100 mV the intrinsic energy of the unliganded gating isomerization is +8.4 kcal/mol (35 kJ/mol), and that in the absence of a membrane potential, the intrinsic chemical energy of this global conformational change is +9.4 kcal/mol (39 kJ/mol). Na⁺ and K⁺ in the extracellular solution have no measureable effect on E₀, which suggests that unliganded gating occurs with only water occupying the transmitter binding sites. The results are discussed with regard to the energy changes in receptor activation and the competitive antagonism of ions in agonist binding.

Drying transition in the hydrophobic gate of the GLIC channel blocks ion conduction

The theoretical prediction of water drying transitions near nonpolar surfaces has stimulated an intensive search for biological processes exploiting this extreme form of hydrophobicity. Here we quantitatively demonstrate that drying of a hydrophobic constriction is the major determinant of ion conductance in the GLIC pentameric ion channel. Molecular-dynamics simulations show that in the closed state, the channel conductance is ∼12 orders-of-magnitude lower than in the open state. This large drop in conductance is remarkable because even in the functionally closed conformation the pore constriction remains wide enough for the passage of sodium ions, aided by a continuous bridge of ∼12 water molecules. However, we find that the free energy cost of hydrating the hydrophobic gate is large, accounting almost entirely for the energetic barrier blocking ion passage. The free energies of transferring a sodium ion into a prehydrated gate in functionally closed and open states differ by only 1.2 kcal/mol, compared to an 11 kcal/mol difference in the costs of hydrating the hydrophobic gate. Conversely, ion desolvation effects play only minor roles in GLIC ion channel gating. Our simulations help rationalize experiments probing the gating kinetics of the nicotinic acetylcholine receptor in response to mutations of pore-lining residues. The molecular character and phase behavior of water should thus be included in quantitative descriptions of ion channel gating.

An outer membrane protein undergoes enthalpy- and entropy-driven transitions

β-Barrel membrane proteins often fluctuate among various open substates, yet the nature of these transitions is not fully understood. Using temperature-dependent, single-molecule electrophysiology analysis, along with rational protein design, we show that OccK1, a member of the outer membrane carboxylate channel from Pseudomonas aeruginosa, features a discrete gating dynamics comprising both enthalpy-driven and entropy-driven current transitions. OccK1 was chosen for the analysis of these transitions, because it is a monomeric transmembrane β-barrel of a known high-resolution crystal structure and displays three distinguishable, time-resolvable open substates. Native and loop-deletion OccK1 proteins showed substantial changes in the activation enthalpies and entropies of the channel transitions, but modest alterations in the equilibrium free energies, confirming that the system never departs from equilibrium. Moreover, some current fluctuations of OccK1 indicated a counterintuitive, negative activation enthalpy, which was compensated by a significant decrease in the activation entropy. Temperature scanning of the single-channel properties of OccK1 exhibited a thermally induced switch of the energetically most favorable open substate at the lowest examined temperature of 4 °C. Therefore, such a semiquantitative assessment of the current fluctuation dynamics not only demonstrates the complexity of channel gating but also reveals distinct functional traits of a β-barrel outer membrane protein under different temperature circumstances.

Sources of energy for gating by neurotransmitters in acetylcholine receptor channels

Nicotinic acetylcholine receptors (AChRs) mediate signaling in the central and peripheral nervous systems. The AChR gating conformational change is powered by a low- to high-affinity change for neurotransmitters at two transmitter binding sites. We estimated (from single-channel currents) the components of energy for gating arising from binding site aromatic residues in the α-subunit. All mutations reduced the energy (TyrC1>>TrpB≈TyrC2>TyrA), with TyrC1 providing ~40% of the total. Considered one at a time, the fractional energy contributions from the aromatic rings were TrpB ~35%, TyrC1 ~28%, TyrC2 ~28%, and TyrA ~10%. Together, TrpB, TyrC1, and TyrC2 comprise an "aromatic triad" that provides much of the total energy from the transmitter for gating. Analysis of mutant pairs suggests that the energy contributions from some residues are nearly independent. Mutations of TyrC1 cause particularly large energy reductions because they remove two favorable and approximately equal interactions between the aromatic ring and the quaternary amine of the agonist and between the hydroxyl and αLysβ7.

Changes in temperature have opposing effects on current amplitude in α7 and α4β2 nicotinic acetylcholine receptors

We have examined the effect of temperature on the electrophysiological properties of three neuronal nicotinic acetylcholine receptor (NACHR) subtypes: the rapidly desensitizing homomeric α7 nAChR, the more slowly desensitizing heteromeric α4β2 nAChR and on α7 nAChRs containing a transmembrane mutation (L247T) that results in dramatically reduced desensitization. In all cases, the functional properties of receptors expressed in Xenopus oocytes at room temperature (RT; 21°C) were compared to those recorded at either physiological temperature (37°C) or at lower temperature (4°C). Alterations in temperature had dramatically differing effects on the amplitude of whole-cell responses detected with these three nAChR subtypes. Compared to responses at RT, the amplitude of agonist-evoked responses with α4β2 nAChRs was increased at high temperature (125±9%, n = 6, P<0.01) and reduced at low temperature (47±5%, n = 6, P<0.01), whereas the amplitude of α7 responses was reduced at high temperature (27±7%, n = 11, P<0.001) and increased at low temperatures (224±16%, n = 10, P<0.001). In contrast to the effects of temperature on α4β2 and wild type α7 nAChRs, the amplitude of α7 nAChRs containing the L247T mutation was unaffected by changes in temperature. In addition, changes in temperature had little or no effect on current amplitude when α7 nAChRs were activated by the largely non-desensitizing allosteric agonist 4BP-TQS. Despite these differing effects of temperature on the amplitude of agonist-evoked responses in different nAChRs, changes in temperature had a consistent effect on the rate of receptor desensitization on all subtypes examined. In all cases, higher temperature resulted in increased rates of desensitization. Thus, it appears that the differing effects of temperature on the amplitudes of whole-cell responses cannot be explained by temperature-induced changes in receptor desensitization rates.

Mapping heat exchange in an allosteric protein

Nicotinic acetylcholine receptors (AChRs) are synaptic ion channels that spontaneously isomerize (i.e., gate) between resting and active conformations. We used single-molecule electrophysiology to measure the temperature dependencies of mouse neuromuscular AChR gating rate and equilibrium constants. From these we estimated free energy, enthalpy, and entropy changes caused by mutations of amino acids located between the transmitter binding sites and the middle of the membrane domain. The range of equilibrium enthalpy change (13.4 kcal/mol) was larger than for free energy change (5.5 kcal/mol at 25°C). For two residues, the slope of the rate-equilibrium free energy relationship (Φ) was approximately constant with temperature. Mutant cycle analysis showed that both free energies and enthalpies are additive for energetically independent mutations. We hypothesize that changes in energy associated with changes in structure mainly occur close to the site of the mutation, and, hence, that it is possible to make a residue-by-residue map of heat exchange in the AChR gating isomerization. The structural correlates of enthalpy changes are discussed for 12 different mutations in the protein.

Design and control of acetylcholine receptor conformational change

Allosteric proteins use energy derived from ligand binding to promote a global change in conformation. The "gating" equilibrium constant of acetylcholine receptor-channels (AChRs) is influenced by ligands, mutations, and membrane voltage. We engineered AChRs to have specific values of this constant by combining these perturbations, and then calculated the corresponding values for a reference condition. AChRs were designed to have specific rate and equilibrium constants simply by adding multiple, energetically independent mutations with known effects on gating. Mutations and depolarization (to remove channel block) changed the diliganded gating equilibrium constant only by changing the unliganded gating equilibrium constant (E(0)) and did not alter the energy from ligand binding. All of the tested perturbations were approximately energetically independent. We conclude that naturally occurring mutations mainly adjust E(0) and cause human disease because they generate AChRs that have physiologically inappropriate values of this constant. The results suggest that the energy associated with a structural change of a side chain in the gating isomerization is dissipated locally and is mainly independent of rigid body or normal mode motions of the protein. Gating rate and equilibrium constants are estimated for seven different AChR agonists using a stepwise engineering approach.