Phi-value analysis of a linear, sequential reaction mechanism: theory and application to ion channel 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.

STRUCTURE, GATING, AND REGULATION OF THE CFTR ANION CHANNEL

The cystic fibrosis transmembrane conductance regulator (CFTR) belongs to the ATP binding cassette (ABC) transporter superfamily but functions as an anion channel crucial for salt and water transport across epithelial cells. CFTR dysfunction, because of mutations, causes cystic fibrosis (CF). The anion-selective pore of the CFTR protein is formed by its two transmembrane domains (TMDs) and regulated by its cytosolic domains: two nucleotide binding domains (NBDs) and a regulatory (R) domain. Channel activation requires phosphorylation of the R domain by cAMP-dependent protein kinase (PKA), and pore opening and closing (gating) of phosphorylated channels is driven by ATP binding and hydrolysis at the NBDs. This review summarizes available information on structure and mechanism of the CFTR protein, with a particular focus on atomic-level insight gained from recent cryo-electron microscopic structures and on the molecular mechanisms of channel gating and its regulation. The pharmacological mechanisms of small molecules targeting CFTR's ion channel function, aimed at treating patients suffering from CF and other diseases, are briefly discussed.

The conformational wave in capsaicin activation of transient receptor potential vanilloid 1 ion channel

The capsaicin receptor TRPV1 has been intensively studied by cryo-electron microscopy and functional tests. However, though the apo and capsaicin-bound structural models are available, the dynamic process of capsaicin activation remains intangible, largely due to the lack of a capsaicin-induced open structural model and the low occupancy of the transition states. Here we report that reducing temperature toward the freezing point substantially increased channel closure events even in the presence of saturating capsaicin. We further used a combination of fluorescent unnatural amino acid (fUAA) incorporation, computational modeling, and rate-equilibrium linear free-energy relationships analysis (Φ-analysis) to derive the fully open capsaicin-bound state model, and reveal how the channel transits from the apo to the open state. We observed that capsaicin initiates a conformational wave that propagates through the S4–S5 linker towards the S6 bundle and finally reaching the selectivity filter. Our study provides a temporal mechanism for capsaicin activation of TRPV1.

The capsaicin receptor TRPV1 has been structurally characterized, but the capsaicin activation dynamics remain elusive. Here authors use fluorescent unnatural amino acid incorporation, computational modeling and Φ-analysis to derive the capsaicin-bound open state model and reveal the capsaicin induced conformational changes.

Asymmetry of movements in CFTR's two ATP sites during pore opening serves their distinct functions

CFTR, the chloride channel mutated in cystic fibrosis (CF) patients, is opened by ATP binding to two cytosolic nucleotide binding domains (NBDs), but pore-domain mutations may also impair gating. ATP-bound NBDs dimerize occluding two nucleotides at interfacial binding sites; one site hydrolyzes ATP, the other is inactive. The pore opens upon tightening, and closes upon disengagement, of the catalytic site following ATP hydrolysis. Extent, timing, and role of non-catalytic-site movements are unknown. Here we exploit equilibrium gating of a hydrolysis-deficient mutant and apply Φ value analysis to compare timing of opening-associated movements at multiple locations, from the cytoplasmic ATP sites to the extracellular surface. Marked asynchrony of motion in the two ATP sites reveals their distinct roles in channel gating. The results clarify the molecular mechanisms of functional cross-talk between canonical and degenerate ATP sites in asymmetric ABC proteins, and of the gating defects caused by two common CF mutations.

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

Gupta et al. use single-channel electrophysiology to investigate the gating mechanism of acetylcholine receptor ion channels. They propose that channel opening starts at the M2–M3 linker and ligand-binding sites and proceeds through four brief intermediate conformations before ending with the collapse of a gate bubble.

Nicotinic acetylcholine receptors are allosteric proteins that generate membrane currents by isomerizing (“gating”) between resting and active conformations under the influence of neurotransmitters. Here, to explore the mechanisms that link the transmitter-binding sites (TBSs) with the distant gate, we use mutant cycle analyses to measure coupling between residue pairs, phi value analyses to sequence domain rearrangements, and current simulations to reproduce a microsecond shut component (“flip”) apparent in single-channel recordings. Significant interactions between amino acids separated by >15 Å are rare; an exception is between the αM2–M3 linkers and the TBSs that are ∼30 Å apart. Linker residues also make significant, local interactions within and between subunits. Phi value analyses indicate that without agonists, the linker is the first region in the protein to reach the gating transition state. Together, the phi pattern and flip component suggest that a complete, resting↔active allosteric transition involves passage through four brief intermediate states, with brief shut events arising from sojourns in all or a subset. We derive energy landscapes for gating with and without agonists, and propose a structure-based model in which resting→active starts with spontaneous rearrangements of the M2–M3 linkers and TBSs. These conformational changes stabilize a twisted extracellular domain to promote transmembrane helix tilting, gate dilation, and the formation of a “bubble” that collapses to initiate ion conduction. The energy landscapes suggest that twisting is the most energetically unfavorable step in the resting→active conformational change and that the rate-limiting step in the reverse process is bubble formation.

Molecular simulation study of the unbinding of α-conotoxin [ϒ4E]GID at the α7 and α4β2 neuronal nicotinic acetylcholine receptors

The α7 and α4β2 neuronal nicotinic receptors belonging to the family of ligand-gated ion channels are most prevalent in the brain, and are implicated in various neurodegenerative disorders. α-conotoxin GID (and its analogue [ϒ4E]GID) specifically inhibits these subtypes, with more affinity towards the human α7 (hα7) subtype, and is valuable in understanding the physiological roles of these receptors. In this study, we use umbrella-sampling molecular dynamics simulations to understand the mechanism of interaction between [ϒ4E]GID and the agonist binding pockets of the α4β2 and the hα7 receptors, and to estimate their relative binding affinities (ΔG_bind). The obtained ΔG_bind values indicate stronger interaction with the hα7 receptor, in agreement with previous experimental studies. Simulations also revealed different unbinding pathways between the two receptor subtypes, enabling identification of a number of interactions at locations far from the orthosteric binding site which may explain the difference in [ϒ4E]GID potency. The pathways identified will help in the design of novel conotoxins with increased potency at α4β2, for which there is currently no known highly potent conotoxin inhibitor. Computational mutational free energy analyses also revealed a number of possible single-site mutations to GID which might enhance its selective binding to α4β2 over α7.

Timing of CFTR pore opening and structure of its transition state

In CFTR, the chloride ion channel mutated in cystic fibrosis (CF) patients, pore opening is coupled to ATP-binding-induced dimerization of two cytosolic nucleotide binding domains (NBDs) and closure to dimer disruption following ATP hydrolysis. CFTR opening rate, unusually slow because of its high-energy transition state, is further slowed by CF mutation ΔF508. Here, we exploit equilibrium gating of hydrolysis-deficient CFTR mutant D1370N and apply rate-equilibrium free-energy relationship analysis to estimate relative timing of opening movements in distinct protein regions. We find clear directionality of motion along the longitudinal protein axis and identify an opening transition-state structure with the NBD dimer formed but the pore still closed. Thus, strain at the NBD/pore-domain interface, the ΔF508 mutation locus, underlies the energetic barrier for opening. Our findings suggest a therapeutic opportunity to stabilize this transition-state structure pharmacologically in ΔF508-CFTR to correct its opening defect, an essential step toward restoring CFTR function.

Insights into the gating mechanism of the ryanodine-modified human cardiac Ca2+-release channel (ryanodine receptor 2)

Ryanodine receptors (RyRs) are intracellular membrane channels playing key roles in many Ca(2+) signaling pathways and, as such, are emerging novel therapeutic and insecticidal targets. RyRs are so named because they bind the plant alkaloid ryanodine with high affinity and although it is established that ryanodine produces profound changes in all aspects of function, our understanding of the mechanisms underlying altered gating is minimal. We address this issue using detailed single-channel gating analysis, mathematical modeling, and energetic evaluation of state transitions establishing that, with ryanodine bound, the RyR pore adopts an extremely stable open conformation. We demonstrate that stability of this state is influenced by interaction of divalent cations with both activating and inhibitory cytosolic sites and, in the absence of activating Ca(2+), trans-membrane voltage. Comparison of the conformational stability of ryanodine- and Imperatoxin A-modified channels identifies significant differences in the mechanisms of action of these qualitatively similar ligands.

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.

Hydrophobic interactions between the voltage sensor and pore mediate inactivation in Kv11.1 channels

Kv11.1 channels are critical for the maintenance of a normal heart rhythm. The flow of potassium ions through these channels is controlled by two voltage-regulated gates, termed "activation" and "inactivation," located at opposite ends of the pore. Crucially in Kv11.1 channels, inactivation gating occurs much more rapidly, and over a distinct range of voltages, compared with activation gating. Although it is clear that the fourth transmembrane segments (S4), within each subunit of the tetrameric channel, are important for controlling the opening and closing of the activation gate, their role during inactivation gating is much less clear. Here, we use rate equilibrium free energy relationship (REFER) analysis to probe the contribution of the S4 "voltage-sensor" helix during inactivation of Kv11.1 channels. Contrary to the important role that charged residues play during activation gating, it is the hydrophobic residues (Leu529, Leu530, Leu532, and Val535) that are the key molecular determinants of inactivation gating. Within the context of an interconnected multi-domain model of Kv11.1 inactivation gating, our REFER analysis indicates that the S4 helix and the S4-S5 linker undergo a conformational rearrangement shortly after that of the S5 helix and S5P linker, but before the S6 helix. Combining REFER analysis with double mutant cycle analysis, we provide evidence for a hydrophobic interaction between residues on the S4 and S5 helices. Based on a Kv11.1 channel homology model, we propose that this hydrophobic interaction forms the basis of an intersubunit coupling between the voltage sensor and pore domain that is an important mediator of inactivation gating.

Agonists binding nicotinic receptors elicit specific channel-opening patterns at αγ and αδ sites

'Embryonic' muscle-type nicotinic acetylcholine receptor channels (nAChRs) bind ligands at interfaces of α- and γ- or δ-subunits. αγ and αδ sites differ in affinity, but their contributions to opening the channel have remained elusive. We compared high-resolution patch clamp currents evoked by epibatidine (Ebd), carbamylcholine (CCh) and acetylcholine (ACh). Ebd binds with 75-fold higher affinity at αγ than at αδ sites, whereas CCh and ACh prefer αδ sites. Similar short (τO1), intermediate (τO2) and long (τO3) types of opening were observed with all three agonists. τO2 openings were maximally prevalent at low Ebd concentrations, binding at αγ sites. By contrast, τO1 openings appear to be generated at αδ sites. In addition, two types of burst appeared: short bursts of an average of 0.75 ms (τB1) that should arise from the αγ site, and long bursts of 12-25 ms (τB2) in duration arising from double liganded receptors. Limited by the temporal resolution, the closings within bursts were invariant at 3 μs. Corrected for missed closings, in the case of ACh the openings within long bursts lasted 170 μs and those in short bursts about 30 μs. Blocking αδ sites with α-conotoxin M1 (CTx) eliminated both τO1 and τB2 and left only τO2 and the short τB1 bursts, as expected. Furthermore we found desensitization when the receptors bound ACh only at the αγ site. When CTx was applied to 'embryonic' mouse endplates, monoquantal current rise times were increased, and amplitude and decay time constants were reduced, as expected. Thus the αγ and αδ sites of nAChRs elicit specific channel-opening patterns.

Conformational changes in acetylcholine binding protein investigated by temperature accelerated molecular dynamics

Despite the large number of studies available on nicotinic acetylcholine receptors, a complete account of the mechanistic aspects of their gating transition in response to ligand binding still remains elusive. As a first step toward dissecting the transition mechanism by accelerated sampling techniques, we study the ligand-induced conformational changes of the acetylcholine binding protein (AChBP), a widely accepted model for the full receptor extracellular domain. Using unbiased Molecular Dynamics (MD) and Temperature Accelerated Molecular Dynamics (TAMD) simulations we investigate the AChBP transition between the apo and the agonist-bound state. In long standard MD simulations, both conformations of the native protein are stable, while the agonist-bound structure evolves toward the apo one if the orientation of few key sidechains in the orthosteric cavity is modified. Conversely, TAMD simulations initiated from the native conformations are able to produce the spontaneous transition. With respect to the modified conformations, TAMD accelerates the transition by at least a factor 10. The analysis of some specific residue-residue interactions points out that the transition mechanism is based on the disruption/formation of few key hydrogen bonds. Finally, while early events of ligand dissociation are observed already in standard MD, TAMD accelerates the ligand detachment and, at the highest TAMD effective temperature, it is able to produce a complete dissociation path in one AChBP subunit.

Energetic contributions to channel gating of residues in the muscle nicotinic receptor β1 subunit

In the pentameric ligand-gated ion channel family, transmitter binds in the extracellular domain and conformational changes result in channel opening in the transmembrane domain. In the muscle nicotinic receptor and other heteromeric members of the family one subunit does not contribute to the canonical agonist binding site for transmitter. A fundamental question is whether conformational changes occur in this subunit. We used records of single channel activity and rate-equilibrium free energy relationships to examine the β1 (non-ACh-binding) subunit of the muscle nicotinic receptor. Mutations to residues in the extracellular domain have minimal effects on the gating equilibrium constant. Positions in the channel lining (M2 transmembrane) domain contribute strongly and relatively late during gating. Positions thought to be important in other subunits in coupling the transmitter-binding to the channel domains have minimal effects on gating. We conclude that the conformational changes involved in channel gating propagate from the binding-site to the channel in the ACh-binding subunits and subsequently spread to the non-binding subunit.

Conserved N-terminal negative charges support optimally efficient N-type inactivation of Kv1 channels

N-type inactivation is produced by the binding of a potassium channel's N-terminus within the open pore, blocking conductance. Previous studies have found that introduction of negative charges into N-terminal inactivation domains disrupts inactivation; however, the Aplysia AKv1 N-type inactivation domain contains two negatively charged residues, E2 and E9. Rather than being unusual, sequence analysis shows that this N-terminal motif is highly conserved among Kv1 sequences across many phyla. Conservation analysis shows some tolerance at position 9 for other charged residues, like D9 and K9, whereas position 2 is highly conserved as E2. To examine the functional importance of these residues, site directed mutagenesis was performed and effects on inactivation were recorded by two electrode voltage clamp in Xenopus oocytes. We find that inclusion of charged residues at positions 2 and 9 prevents interactions with non-polar sites along the inactivation pathway increasing the efficiency of pore block. In addition, E2 appears to have additional specific electrostatic interactions that stabilize the inactivated state likely explaining its high level of conservation. One possible explanation for E2's unique importance, consistent with our data, is that E2 interacts electrostatically with a positive charge on the N-terminal amino group to stabilize the inactivation domain at the block site deep within the pore. Simple electrostatic modeling suggests that due to the non-polar environment in the pore in the blocked state, even a 1 Å larger separation between these charges, produced by the E2D substitution, would be sufficient to explain the 65× reduced affinity of the E2D N-terminus for the pore. Finally, our studies support a multi-step, multi-site N-type inactivation model where the N-terminus interacts deep within the pore in an extended like structure placing the most N-terminal residues 35% of the way across the electric field in the pore blocked state.

An outline of desensitization in pentameric ligand-gated ion channel receptors

Pentameric ligand-gated ion channel (pLGIC) receptors exhibit desensitization, the progressive reduction in ionic flux in the prolonged presence of agonist. Despite its pathophysiological importance and the fact that it was first described over half a century ago, surprisingly little is known about the structural basis of desensitization in this receptor family. Here, we explain how desensitization is defined using functional criteria. We then review recent progress into reconciling the structural and functional basis of this phenomenon. The extracellular-transmembrane domain interface is a key locus. Activation is well known to involve conformational changes at this interface, and several lines of evidence suggest that desensitization involves a distinct conformational change here that is incompatible with activation. However, major questions remain unresolved, including the structural basis of the desensitization-induced agonist affinity increase and the mechanism of pore closure during desensitization.

Pore helices play a dynamic role as integrators of domain motion during Kv11.1 channel inactivation gating

Proteins that form ion-selective pores in the membrane of cells are integral to many rapid signaling processes, including regulating the rhythm of the heartbeat. In potassium channels, the selectivity filter is critical for both endowing an exquisite selectivity for potassium ions, as well as for controlling the flow of ions through the pore. Subtle rearrangements in the complex hydrogen-bond network that link the selectivity filter to the surrounding pore helices differentiate conducting (open) from nonconducting (inactivated) conformations of the channel. Recent studies suggest that beyond the selectivity filter, inactivation involves widespread rearrangements of the channel protein. Here, we use rate equilibrium free energy relationship analysis to probe the structural changes that occur during selectivity filter gating in Kv11.1 channels, at near atomic resolution. We show that the pore helix plays a crucial dynamic role as a bidirectional interface during selectivity filter gating. We also define the molecular bases of the energetic coupling between the pore helix and outer helix of the pore domain that occurs early in the transition from open to inactivated states, as well as the coupling between the pore helix and inner helix late in the transition. Our data demonstrate that the pore helices are more than just static structural elements supporting the integrity of the selectivity filter; instead they play a crucial dynamic role during selectivity filter gating.

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.

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.

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.

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.

Discrete-state representation of ion permeation coupled to fast gating in a model of CLC-chloride channels: analytic estimation of the state-to-state rate constants

Analytical estimation of state-to-state rate constants is carried out for a recently developed discrete state model of chloride ion motion in a CLC chloride channel (Coalson and Cheng, J. Phys. Chem. B 2010, 114, 1424). In the original presentation of this model, the same rate constants were evaluated via three-dimensional Brownian dynamics simulations. The underlying dynamical theory is an appropriate single- or multiparticle three-dimensional Smoluchowski equation. Taking advantage of approximate geometric symmetries (based on the details of the model channel geometry), well-known formulas for state-to-state transition rates are appealed to herein and adapted as necessary to the problem at hand. Rates of ionic influx from a bulk electrolyte reservoir to the nearest binding site within the channel pore are particularly challenging to compute analytically because they reflect multi-ion interactions (as opposed to single-ion dynamics). A simple empirical correction factor is added to the single-ion rate constant formula in this case to account for the saturation of influx rate constants with increasing bulk Cl(-) concentration. Overall, the agreement between all analytically estimated rate constants is within a factor of 2 of those computed via three-dimensional Brownian dynamics simulations, and often better than this. Current-concentration curves obtained using rate constants derived from these two different computational approaches agree to within 25%.

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.

Mapping the sequence of conformational changes underlying selectivity filter gating in the K(v)11.1 potassium channel

The potassium channel selectivity filter both discriminates between K(+) and sodium ions and contributes to gating of ion flow. Static structures of conducting (open) and nonconducting (inactivated) conformations of this filter are known; however, the sequence of protein rearrangements that connect these two states is not. We show that closure of the selectivity filter gate in the human K(v)11.1 K(+) channel (also known as hERG, for ether-a-go-go-related gene), a key regulator of the rhythm of the heartbeat, is initiated by K(+) exit, followed in sequence by conformational rearrangements of the pore domain outer helix, extracellular turret region, voltage sensor domain, intracellular domains and pore domain inner helix. In contrast to the simple wave-like sequence of events proposed for opening of ligand-gated ion channels, a complex spatial and temporal sequence of widespread domain motions connect the open and inactivated states of the K(v)11.1 K(+) channel.

Relationships between structural dynamics and functional kinetics in oligomeric membrane receptors

Recent efforts to broaden understanding of the molecular mechanisms of membrane receptors in signal transduction make use of rate-equilibrium free-energy relationships (REFERs), previously applied to chemical reactions, enzyme kinetics, and protein folding. For oligomeric membrane receptors, we distinguish between a), the Leffler parameter α_L, to characterize the global transition state for the interconversion between conformations; and b), the Fersht parameter, ϕ_F, to assign the degree of progression of individual residue positions at the transition state. For both α_L and ϕ_F, insights are achieved by using harmonic energy profiles to reflect the dynamic nature of proteins, as illustrated with single-channel results reported for normal and mutant nicotinic receptors. We also describe new applications of α_L based on published results. For large-conductance calcium-activated potassium channels, data are satisfactorily fit with an α_L value of 0.65, in accord with REFERs. In contrast, results reported for the flip conformational state of glycine and nicotinic receptors are in disaccord with REFERs, since they yield α_L values outside the usual limits of 0–1. Concerning published ϕ_F values underlying the conformational wave hypothesis for nicotinic receptors, we note that interpretations may be complicated by variations in the width of harmonic energy profiles.

Interdependence of receptor activation and ligand binding in HCN2 pacemaker channels

HCN pacemaker channels are tetramers mediating rhythmicity in neuronal and cardiac cells. The activity of these channels is controlled by both membrane voltage and the ligand cAMP, binding to each of the four channel subunits. The molecular mechanism underlying channel activation and the relationship between the two activation stimuli are still unknown. Using patch-clamp fluorometry and a fluorescent cAMP analog, we show that full ligand-induced activation appears already with only two ligands bound to the tetrameric channel. Kinetic analysis of channel activation and ligand binding suggests direct interaction between the voltage sensor and the cyclic nucleotide-binding domain, bypassing the pore. By exploiting the duality of activation in HCN2 channels by voltage and ligand binding, we quantify the increase of the binding affinity and overall free energy for binding upon channel activation, proving thus the principle of reciprocity between ligand binding and conformational change in a receptor protein.

Energetics of gating at the apo-acetylcholine receptor transmitter binding site

Acetylcholine receptor channels switch between conformations that have a low versus high affinity for the transmitter and conductance for ions (R↔R*; gating). The forward isomerization, which begins at the transmitter binding sites and propagates ∼50 Å to the narrow region of the pore, occurs by approximately the same sequence of molecular events with or without agonists present at the binding sites. To pinpoint the forces that govern the R versus R* agonist affinity ratio, we measured single-channel activation parameters for apo-receptors having combinations of mutations of 10 transmitter binding site residues in the α (Y93, G147, W149, G153, Y190, C192, and Y198), ε (W55 and P121), or δ (W57) subunit. Gating energy changes were largest for the tryptophan residues. The αW149 energy changes were coupled with those of the other aromatic amino acids. Mutating the aromatic residues to Phe reduces the R/R* equilibrium dissociation constant ratio, with αY190 and αW149 being the most sensitive positions. Most of the mutations eliminated long-lived spontaneous openings. The results provide a foundation for understanding how ligands trigger protein conformational change.

Discrete-state representation of ion permeation coupled to fast gating in a model of ClC chloride channels: comparison to multi-ion continuous space Brownian dynamics simulations

A discrete-state model of chloride ion motion in a ClC chloride channel is constructed, following a previously developed multi-ion continuous space model of the same system (Cheng, M. H.; Mamonov, A. B.; Dukes, J. W.; Coalson, R. D. J. Phys. Chem. B 2007, 111, 5956) that included a simplistic representation of the fast gate in this channel. The reducibility of the many-body continuous space to the eight discrete-state model considered in the present work is examined in detail by performing three-dimensional Brownian dynamics simulations of each allowed state-to-state transition in order to extract the appropriate rate constant for this process, and then inserting the pairwise rate constants thereby obtained into an appropriate set of kinetic master equations. Experimental properties of interest, including the rate of Cl(-) ion permeation through the open channel and the average rate of closing of the fast gate as a function of bulk Cl(-) ion concentrations in the intracellular and extracellular electrolyte reservoirs are computed. Good agreement is found between the results obtained via the eight discrete-state model versus the multi-ion continuous space model, thereby encouraging continued development of the discrete-state model to include more complex behaviors observed experimentally in these channels.

Binding, activation and modulation of Cys-loop receptors

It is over forty years since the major neurotransmitters and their protein receptors were identified, and over twenty years since determination of the first amino-acid sequences of the Cys-loop receptors that recognize acetylcholine, serotonin, GABA and glycine. The last decade has seen the first structures of these proteins (and related bacterial and molluscan homologues) determined to atomic resolution. Hopefully over the next decade, more detailed molecular structures of entire Cys-loop receptors in drug-bound and drug-free conformations will become available. These, together with functional studies, will provide a clear picture of how these receptors participate in neurotransmission and how structural variations between receptor subtypes impart their unique characteristics. This insight should facilitate the design of novel and improved therapeutics to treat neurological disorders. This review considers our current understanding about the processes of agonist binding, receptor activation and channel opening, as well as allosteric modulation of the Cys-loop receptor family.

The gating isomerization of neuromuscular acetylcholine receptors

Acetylcholine receptor-channels are allosteric proteins that isomerize ('gate') between conformations that have a low vs. high affinity for the transmitter and conductance for ions. In order to comprehend the mechanism by which the affinity and conductance changes are linked it is of value to know the magnitude, timing and distribution of energy flowing through the system. Knowing both the di- and unliganded gating equilibrium constants (E(2) and E(0)) is a foundation for understanding the AChR gating mechanism and for engineering both the ligand and the protein to operate in predictable ways. In adult mouse neuromuscular receptors activated by acetylcholine, E(2) = 28 and E(0) approximately 6.5 x 10(7). At each (equivalent) transmitter binding site acetylcholine provides approximately 5.2 kcal mol(1) to motivate the isomerization. The partial agonist choline provides approximately 3.3 kcal mol(1). The relative time of a residue's gating energy change is revealed by the slope of its rate-equilibrium constant relationship. A map of this parameter suggests that energy propagates as a conformational cascade between the transmitter binding sites and the gate region. Although gating energy changes are widespread throughout the protein, some residues are particularly sensitive to perturbations. Several specific proposals for the structural events that comprise the gating conformational cascade are discussed.

Application of rate-equilibrium free energy relationship analysis to nonequilibrium ion channel gating mechanisms

Rate-equilibrium free energy relationship (REFER) analysis provides information on transition-state structures and has been applied to reveal the temporal sequence in which the different regions of an ion channel protein move during a closed-open conformational transition. To date, the theory used to interpret REFER relationships has been developed only for equilibrium mechanisms. Gating of most ion channels is an equilibrium process, but recently several ion channels have been identified to have retained nonequilibrium traits in their gating cycles, inherited from transporter-like ancestors. So far it has not been examined to what extent REFER analysis is applicable to such systems. By deriving the REFER relationships for a simple nonequilibrium mechanism, this paper addresses whether an equilibrium mechanism can be distinguished from a nonequilibrium one by the characteristics of their REFER plots, and whether information on the transition-state structures can be obtained from REFER plots for gating mechanisms that are known to be nonequilibrium cycles. The results show that REFER plots do not carry information on the equilibrium nature of the underlying gating mechanism. Both equilibrium and nonequilibrium mechanisms can result in linear or nonlinear REFER plots, and complementarity of REFER slopes for opening and closing transitions is a trivial feature true for any mechanism. Additionally, REFER analysis provides limited information about the transition-state structures for gating schemes that are known to be nonequilibrium cycles.

Gating of nicotinic ACh receptors: latest insights into ligand binding and function

Nicotinic acetylcholine receptors (nAChRs) are in the superfamily of cys-loop receptors, and are widely expressed in the nervous system where they participate in a variety of physiological functions, including regulating excitability and neurotransmitter release, as well as neuromuscular contraction. Members of the cys-loop family of receptors, which also includes the molluscan ACh-binding protein (AChBP), a soluble protein that is analogous to the extracellular ligand-binding domain of the cys-loop receptors, are pentameric assemblies of five subunits, with each subunit arranged around a central pore. The binding of ACh to the extracellular interface between two subunits induces channel opening. With the recent 4 A resolution of the Torpedo nAChR, and the crystal structure of the AChBP, much has been learned about the structure of the ligand-binding domain and the channel pore, as well as major structural rearrangements that may confer channel opening, including a major rearrangement of the C-loop within the ligand binding pocket, and perhaps other regions including the F-loop (the beta8-beta9 linker), the beta1-beta2 linker and the cys-loop. Here I will review the latest findings from my lab aimed at a further understanding of the function of the neuronal nAChR channels (and in particular the role of desensitization), and our search for novel AChBP species that may lead to a further understanding of the function of the cys-loop receptor family.

Structural basis of activation of cys-loop receptors: the extracellular-transmembrane interface as a coupling region

Cys-loop receptors mediate rapid transmission throughout the nervous system by converting a chemical signal into an electric one. They are pentameric proteins with an extracellular domain that carries the transmitter binding sites and a transmembrane region that forms the ion pore. Their essential function is to couple the binding of the agonist at the extracellular domain to the opening of the ion pore. How the structural changes elicited by agonist binding are propagated through a distance of 50 A to the gate is therefore central for the understanding of the receptor function. A step forward toward the identification of the structures involved in gating has been given by the recently elucidated high-resolution structures of Cys-loop receptors and related proteins. The extracellular-transmembrane interface has attracted attention because it is a structural transition zone where beta-sheets from the extracellular domain merge with alpha-helices from the transmembrane domain. Within this zone, several regions form a network that relays structural changes from the binding site toward the pore, and therefore, this interface controls the beginning and duration of a synaptic response. In this review, the most recent findings on residues and pairwise interactions underlying channel gating are discussed, the main focus being on the extracellular-transmembrane interface.

Effects of jamming on nonequilibrium transport times in nanochannels

Many biological channels perform highly selective transport without direct input of metabolic energy and without transitions from a "closed" to an "open" state during transport. Mechanisms of selectivity of such channels serve as an inspiration for creation of artificial nanomolecular sorting devices and biosensors. To elucidate the transport mechanisms, it is important to understand the transport on the single molecule level in the experimentally relevant regime when multiple particles are crowded in the channel. In this Letter we analyze the effects of interparticle crowding on the nonequilibrium transport times through a finite-length channel by means of analytical theory and computer simulations.

Energy and structure of the M2 helix in acetylcholine receptor-channel gating

We studied single-channel currents from neuromuscular acetylcholine receptor-channels with mutations in the pore-lining, M2 helix of the epsilon-subunit. Three parameters were quantified: 1), the diliganded gating equilibrium constant (E(2)), which reflects the energy difference between C(losed) and O(pen) conformations; 2), the correlation between the opening rate constant and E(2) on a log-log scale (Phi), which illuminates the energy character of the residue (C- versus O-like) within the C<-->O isomerization process; and 3), the open-channel current amplitude (i(0)), which reports whether a mutation alters the energetics of ion permeation. The largest E(2) changes were observed in the cytoplasmic half of epsilonM2 (5', 9', 12', 13', and 16'), with smaller changes apparent for residues > or =17'. Phi was approximately 0.54 for most epsilonM2 residues, but was approximately 0.32 at the positions that had largest E(2) changes. An arginine substitution reduced i(0) significantly at six positions, with the magnitude of the reduction increasing, 16'-->2'. The measurements suggest that the 9', 12', and 13' residues experience large and late free-energy changes in the channel-opening process. We speculate that in the gating isomerization the pore-facing residues >6' and <16' experience multiple energy perturbations associated with changes in protein structure and, perhaps, hydration.

Probability distributions of Markovian sodium channel states during propagating axonal impulses with or without recovery supernormality

This study addressed a macroscopic neurophysiological phenomenon - supernormality during the recovery phase of propagating axonal impulses - in explicit chemical terms. Excitation was reconstructed numerically using the kinetic scheme of multiple-state probabilistic transitions within a population of voltage-dependent sodium channels (NaCh) derived by Vandenberg and Bezanilla ("PC" scheme). Each NaCh transition was characterized as a reversible Markov process with voltage-dependent rate constants associated with each respective directional transition. While recovery reconstructed with the Hodgkin-Huxley formalism included a supernormal period, the PC scheme did not. The present analysis showed that the occurrence and degree of supernormality with the PC scheme was determined by the relative speed of the transitions within the closed loop of the kinetic scheme; supernormality was promoted by speeding these kinetics. The analysis also showed that concurrent with supernormality, the faster loop kinetics caused (1) an elevation in the C(1) --> C(2) transitions, and (2) a reduction in the I(4) --> I(5) transitions. Thus, macroscopic functionality in information processing could be expressed in terms of probabilistic interstate transitions among a population of NaCh molecules.

Recent advances in understanding the structure of nicotinic acetylcholine receptors

Nicotinic acetylcholine receptors (nAChRs), members of the Cys-loop ligand-gated ion channels (LGICs) superfamily, are involved in signal transduction upon binding of the neurotransmitter acetylcholine or exogenous ligands, such as nicotine. nAChRs are pentameric assemblies of homologous subunits surrounding a central pore that gates cation flux, and are expressed at the neuromuscular junction and in the nervous system and several nonneuronal cell types. The 17 known nAChR subunits assemble into a variety of pharmacologically distinct receptor subtypes. nAChRs are implicated in a range of physiological functions and pathophysiological conditions related to muscle contraction, learning and memory, reward, motor control, arousal, and analgesia, and therefore present an important target for drug research. Such studies would be greatly facilitated by knowledge of the high-resolution structure of the nAChR. Although this information is far from complete, important progress has been made mainly based on electron microscopy studies of Torpedo nAChR and the high-resolution X-ray crystal structures of the homologous molluscan acetylcholine-binding proteins, the extracellular domain of the mouse nAChR alpha1 subunit, and two prokaryotic pentameric LGICs. Here, we review some of the latest advances in our understanding of nAChR structure and gating.

Multiple intermediate states precede pore block during N-type inactivation of a voltage-gated potassium channel

N-type inactivation of voltage-gated potassium channels is an autoinhibitory process that occurs when the N terminus binds within the channel pore and blocks conduction. N-type inactivation and recovery occur with single-exponential kinetics, consistent with a single-step reaction where binding and block occur simultaneously. However, recent structure-function studies have suggested the presence of a preinactivated state whose formation and loss regulate inactivation and recovery kinetics. Our studies on N-type inactivation of the Shaker-type AKv1 channel support a multiple-step inactivation process involving a series of conformational changes in distinct regions of the N terminus that we have named the polar, flex, and latch regions. The highly charged polar region forms interactions with the surface of the channel leading up to the side window openings between the T1 domain and the channel transmembrane domains, before the rate-limiting step occurs. This binding culminates with a specific electrostatic interaction between R18 and EDE161-163 located at the entrance to the side windows. The latch region appears to work together with the flex region to block the pore after polar region binding occurs. Analysis of tail currents for a latch region mutant shows that both blocked and unblocked states exist after the rate-limiting transition is passed. Our results suggest that at least two intermediate states exist for N-type inactivation: a polar region-bound state that is formed before the rate-limiting step, and a pre-block state that is formed by the flex and latch regions during the rate-limiting step.

Relationship between nucleotide binding and ion channel gating in cystic fibrosis transmembrane conductance regulator

We have employed rate-equilibrium free energy relationship (REFER) analysis to characterize the dynamic events involved in the allosteric regulation of cystic fibrosis transmembrane conductance regulator (CFTR) function. A wide range of different hydrolysable and poorly hydrolysable nucleoside triphosphates were used to elucidate the role of ATP hydrolysis in CFTR function. The linearity of the REFER plots and Phi values near unity for all ligands tested implies that CFTR channel gating is a reversible thermally driven process with all structural reorganization in the binding site(s) completed prior to channel opening. This is consistent with the requirement for nucleotide binding for channel opening. However, the channel structural transition from the open to the closed state occurs independently of any events in the binding sites. Similar results were obtained on substitution of amino acids at coupling joints between both nucleotide binding domains (NBD) and cytoplasmic loops (CL) in opposite halves of the protein, indicating that any structural reorganization there also had occurred in the channel closed state. The fact that fractional Phi values were not observed in either of these distant sites suggests that there may not be a deterministic 'lever-arm' mechanism acting between nucleotide binding sites and the channel gate. These findings favour a stochastic coupling between binding and gating in which all structural transitions are thermally driven processes. We speculate that increase of channel open state probability is due to reduction of the number of the closed state configurations available after physical interaction between ligand bound NBDs and the channel.

Effects of multiple occupancy and interparticle interactions on selective transport through narrow channels: theory versus experiment

Many biological and artificial transport channels function without direct input of metabolic energy during a transport event and without structural rearrangements involving transitions from a closed to an open state. Nevertheless, such channels are able to maintain efficient and selective transport. It has been proposed that attractive interactions between the transported molecules and the channel can increase the transport efficiency and that the selectivity of such channels can be based on the strength of the interaction of the specifically transported molecules with the channel. Herein, we study the transport through narrow channels in a framework of a general kinetic theory, which naturally incorporates multiparticle occupancy of the channel and non-single-file transport. We study how the transport efficiency and the probability of translocation through the channel are affected by interparticle interactions in the confined space inside the channel, and establish conditions for selective transport. We compare the predictions of the model with the available experimental data and find good semiquantitative agreement. Finally, we discuss applications of the theory to the design of artificial nanomolecular sieves.

Aromatic Residues {epsilon}Trp-55 and {delta}Trp-57 and the Activation of Acetylcholine Receptor Channels

The two transmitter binding sites of the neuromuscular acetylcholine (ACh) receptor channel contain several aromatic residues, including a tryptophan located on the complementary, negative face of each binding pocket. These two residues, Trp-55 in the epsilon subunit and Trp-57 in the delta subunit, were mutated (AEFHILRVY), and for most constructs the rate constants for acetylcholine binding and channel gating were estimated by using single channel kinetic analyses. The rate constants for unliganded channel opening and closing were also estimated for some mutants. From these measurements we calculated all of the equilibrium constants of the "allosteric" cycle as follows: diliganded gating, unliganded gating, dissociation from the C(losed) conformation, and dissociation from the O(pen) conformation. The results indicate the following. (i) These aromatic side chains play a relatively minor role in ACh receptor channel activation. (ii) The main consequence of mutations is to reduce the affinity of the O conformation of the binding site for ACh, with the effect being greater at the epsilon subunit. (iii) In epsilon (but not delta) the aromatic nature of the side chain is important in determining affinity, to a slightly greater degree in the O conformation. Phi value analyses (of both tryptophan residues) show Phi approximately 1 for both the ACh binding and diliganded gating reactions. (iv) This suggests that the structural boundaries of the dynamic elements of the gating conformational change may not be subunit-delimited, and (v) the mutated tryptophan residues experience energy changes that occur relatively early in both the ligand-binding and channel-gating reactions.

Unliganded gating of acetylcholine receptor channels

We estimated the unliganded opening and closing rate constants of neuromuscular acetylcholine receptor-channels (AChRs) having mutations that increased the gating equilibrium constant. For some mutant combinations, spontaneous openings occurred in clusters. For 25 different constructs, the unliganded gating equilibrium constant (E(0)) was correlated with the product of the predicted fold-increase in the diliganded gating equilibrium constant caused by each mutation alone. We estimate that (i) E(0) for mouse, wild-type alpha(2)beta delta epsilon AChRs is approximately 1.15 x 10(-7); (ii) unliganded AChRs open for approximately 80 micros, once every approximately 15 min; (iii) the affinity for ACh of the O(pen) conformation is approximately 10 nM, or approximately 15,600 times greater than for the C(losed) conformation; (iv) the ACh-monoliganded gating equilibrium constant is approximately 1.7 x 10(-3); (v) the C-->O isomerization reduces substantially ACh dissociation, but only slightly increases association; and (vi) ACh provides only approximately 0.9 k(B)T more binding energy per site than carbamylcholine but approximately 3.1 k(B)T more than choline, mainly because of a low O conformation affinity. Most mutations of binding site residue alphaW149 increase E(0). We estimate that the mutation alphaW149F reduces the ACh affinity of C only by 13-fold, but of O by 190-fold. Rate-equilibrium free-energy relationships for different regions of the protein show similar slopes (Phi values) for un- vs. diliganded gating, which suggests that the conformational pathway of the gating structural change is fundamentally the same with and without agonists. Agonist binding is a perturbation that (like most mutations) changes the energy, but not the mechanism, of the gating conformational change.

Observations on rate theory for rugged energy landscapes

The potential energy profile for many complex reactions of proteins, such as folding or allosteric conformational change, involves many different scales of molecular motion along the reaction coordinate. Although it is natural to model the dynamics of motion along such rugged energy landscapes as diffusional (the Smoluchowski equation; SE), problems arise because the frictional forces generated by the molecular surround are typically not strong enough to justify the use of the SE. Here, we discuss the fundamental theory behind the SE and note that it may be justified through a master equation when reduced to its continuum limit. However, the SE cannot be used for rough energy landscapes, where the continuum limit is ill defined. Instead, we suggest that one should use a mean first passage time expression derived from a master equation, and show how this approach can be used to glean information about the underlying dynamics of barrier crossing. We note that the potential profile in the SE is that of the microbarriers between conformational substates, and that there is a temperature-dependent, effective friction associated with the long residence time in the microwells that populate the rough landscape. The number of recrossings of the overall barrier is temperature-dependent, governed by the microbarriers and not by the effective friction. We derive an explicit expression for the mean number of recrossings and its temperature dependence. Finally, we note that the mean first passage time can be used as a departure point for measuring the roughness of the landscape.

On the nature of partial agonism in the nicotinic receptor superfamily

Partial agonists are ligands that bind to receptors but produce only a small maximum response even at concentrations where all receptors are occupied. In the case of ligand-activated ion channels, it has been supposed since 1957 that partial agonists evoke a small response because they are inefficient at eliciting the change of conformation between shut and open states of the channel. We have investigated partial agonists for two members of the nicotinic superfamily-the muscle nicotinic acetylcholine receptor and the glycine receptor-and find that the open-shut reaction is similar for both full and partial agonists, but the response to partial agonists is limited by an earlier conformation change ('flipping') that takes place while the channel is still shut. This has implications for the interpretation of structural studies, and in the future, for the design of partial agonists for therapeutic use.

Gating at the mouth of the acetylcholine receptor channel: energetic consequences of mutations in the alphaM2-cap

Gating of nicotinic acetylcholine receptors from a C(losed) to an O(pen) conformation is the initial event in the postsynaptic signaling cascade at the vertebrate nerve-muscle junction. Studies of receptor structure and function show that many residues in this large, five-subunit membrane protein contribute to the energy difference between C and O. Of special interest are amino acids located at the two transmitter binding sites and in the narrow region of the channel, where C↔O gating motions generate a low↔high change in the affinity for agonists and in the ionic conductance, respectively. We have measured the energy changes and relative timing of gating movements for residues that lie between these two locations, in the C-terminus of the pore-lining M2 helix of the α subunit (‘αM2-cap’). This region contains a binding site for non-competitive inhibitors and a charged ring that influences the conductance of the open pore. αM2-cap mutations have large effects on gating but much smaller effects on agonist binding, channel conductance, channel block and desensitization. Three αM2-cap residues (αI260, αP265 and αS268) appear to move at the outset of channel-opening, about at the same time as those at the transmitter binding site. The results suggest that the αM2-cap changes its secondary structure to link gating motions in the extracellular domain with those in the channel that regulate ionic conductance.

Structural answers and persistent questions about how nicotinic receptors work

The electron diffraction structure of nicotinic acetylcholine receptor (nAChR) from Torpedo marmorata and the X-ray crystallographic structure of acetylcholine binding protein (AChBP) are providing new answers to persistent questions about how nAChRs function as biophysical machines and as participants in cellular and systems physiology. New high-resolution information about nAChR structures might come from advances in crystallography and NMR, from extracellular domain nAChRs as high fidelity models, and from prokaryotic nicotinoid proteins. At the level of biophysics, structures of different nAChRs with different pharmacological profiles and kinetics will help describe how agonists and antagonists bind to orthosteric binding sites, how allosteric modulators affect function by binding outside these sites, how nAChRs control ion flow, and how large cytoplasmic domains affect function. At the level of cellular and systems physiology, structures of nAChRs will help characterize interactions with other cellular components, including lipids and trafficking and signaling proteins, and contribute to understanding the roles of nAChRs in addiction, neurodegeneration, and mental illness. Understanding nAChRs at an atomic level will be important for designing interventions for these pathologies.

Acetylcholine receptor gating at extracellular transmembrane domain interface: the "pre-M1" linker

Charged residues in the beta10-M1 linker region ("pre-M1") are important in the expression and function of neuromuscular acetylcholine receptors (AChRs). The perturbation of a salt bridge between pre-M1 residue R209 and loop 2 residue E45 has been proposed as being a principle event in the AChR gating conformational "wave." We examined the effects of mutations to all five residues in pre-M1 (positions M207-P211) plus E45 in loop 2 in the mouse alpha(1)-subunit. M207, Q208, and P211 mutants caused small (approximately threefold) changes in the gating equilibrium constant (K(eq)), but the changes for R209, L210, and E45 were larger. Of 19 different side chain substitutions at R209 on the wild-type background, only Q, K, and H generated functional channels, with the largest change in K(eq) (67-fold) from R209Q. Various R209 mutants were functional on different E45 backgrounds: H, Q, and K (E45A), H, A, N, and Q (E45R), and K, A, and N (E45L). Phi values for R209 (on the E45A background), L210, and E45 were 0.74, 0.35, and 0.80, respectively. Phi values for R209 on the wt and three other backgrounds could not be estimated because of scatter. The average coupling energy between 209/45 side chains (six different pairs) was only -0.33 kcal/mol (for both alpha subunits, combined). Pre-M1 residues are important for expression of functional channels and participate in gating, but the relatively modest changes in closed- vs. open-state energy caused mutations, the weak coupling energy between these residues and the functional activity of several unmatched-charge pairs are not consistent with the perturbation of a salt bridge between R209 and E45 playing the principle role in gating.

Acetylcholine receptor gating at extracellular transmembrane domain interface: the cys-loop and M2-M3 linker

Acetylcholine receptor channel gating is a propagated conformational cascade that links changes in structure and function at the transmitter binding sites in the extracellular domain (ECD) with those at a “gate” in the transmembrane domain (TMD). We used Φ-value analysis to probe the relative timing of the gating motions of α-subunit residues located near the ECD–TMD interface. Mutation of four of the seven amino acids in the M2–M3 linker (which connects the pore-lining M2 helix with the M3 helix), including three of the four residues in the core of the linker, changed the diliganded gating equilibrium constant (K_eq) by up to 10,000-fold (P272 > I274 > A270 > G275). The average Φ-value for the whole linker was ∼0.64. One interpretation of this result is that the gating motions of the M2–M3 linker are approximately synchronous with those of much of M2 (∼0.64), but occur after those of the transmitter binding site region (∼0.93) and loops 2 and 7 (∼0.77). We also examined mutants of six cys-loop residues (V132, T133, H134, F135, P136, and F137). Mutation of V132, H134, and F135 changed K_eq by 2800-, 10-, and 18-fold, respectively, and with an average Φ-value of 0.74, similar to those of other cys-loop residues. Even though V132 and I274 are close, the energetic coupling between I and V mutants of these positions was small (≤0.51 kcal mol⁻¹). The M2–M3 linker appears to be the key moving part that couples gating motions at the base of the ECD with those in TMD. These interactions are distributed along an ∼16-Å border and involve about a dozen residues.

Gating of nicotinic ACh receptors; new insights into structural transitions triggered by agonist binding that induce channel opening

Nicotinic acetylcholine receptors (nAChRs) are in the superfamily of Cys-loop ligand-gated ion channels, and are pentameric assemblies of five subunits, with each subunit arranged around the central ion-conducting pore. The binding of ACh to the extracellular interface between two subunits induces channel opening. With the recent 4 A resolution of the Torpedo nAChR, and the crystal structure of the related molluscan ACh binding protein, much has been learned about the structure of the ligand binding domain and the channel pore, as well as major structural rearrangements that may confer channel opening. For example, the putative pathway coupling agonist binding to channel gating may include a major rearrangement of the C-loop within the ligand binding pocket, and the disruption of a salt bridge between an arginine residue at the end of the beta10 strand and a glutamate residue in the beta1-beta2 linker. Here we will review and discuss the latest structural findings aiming to further refine the transduction pathway linking binding to gating for the nAChR channels, and discuss similarities and differences among the different members of this Cys-loop superfamily of receptors.