Plasticity of acetylcholine receptor gating motions via rate-energy relationships

Controlling Allosteric Networks in Proteins

Allosteric transition, defined as conformational changes induced by ligand binding, is one of the fundamental properties of proteins. Allostery has been observed and characterized in many proteins, and has been recently utilized to control protein function via regulation of protein activity. Here, we review the physical and evolutionary origin of protein allostery, as well as its importance to protein regulation, drug discovery, and biological processes in living systems. We describe recently developed approaches to identify allosteric pathways, connected sets of pairwise interactions that are responsible for propagation of conformational change from the ligand-binding site to a distal functional site. We then present experimental and computational protein engineering approaches for control of protein function by modulation of allosteric sites. As an example of application of these approaches, we describe a synergistic computational and experimental approach to rescue the cystic-fibrosis-associated protein cystic fibrosis transmembrane conductance regulator, which upon deletion of a single residue misfolds and causes disease. This example demonstrates the power of allosteric manipulation in proteins to both elucidate mechanisms of molecular function and to develop therapeutic strategies that rescue those functions. Allosteric control of proteins provides a tool to shine a light on the complex cascades of cellular processes and facilitate unprecedented interrogation of biological systems.

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.

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.

Highly conserved tyrosine 37 stabilizes desensitized states and restricts calcium permeability of ATP-gated P2X3 receptor

Tyrosine 37 in the first transmembrane (TM1) domain is highly conserved in ATP-gated P2X receptors suggesting its fundamental role. We tested whether Y37 contributes to the desensitization of P2X3 receptors, which is currently not well understood. By combining electrophysiological, imaging and modeling approaches, we studied desensitization of various Y37 P2X3 mutants and potential partners of Y37. Unlike the membrane current of the WT receptor, which desensitized in seconds, Y37A mutant current did not fully desensitize even after minutes-long applications of β,γ-meATP, α,β-meATP, ATP or 2MeS-ATP. The fractional calcium current was enhanced in the Y37A mutant. Y37F did not rescue the native P2X3 phenotype indicating a role for the hydroxyl group of Y37 for the WT receptor. Homology modeling indicated I318 or I319 in TM2 as potential partners for Y37 in the receptor closed state. We tested this hypothesis by creating a permanent interaction between the two residues via disulfide bond. Whereas single Y37C, I318C and I319C mutants were functional, the double mutants Y37C-I318C and Y37C-I319C were non-functional. Using a cyclic model of receptor operation, we suggest that the conserved tyrosine 37 links TM1 to TM2 of adjacent subunit to stabilize desensitized states and restricts calcium permeability through the ion channel.

Role of the ectodomain serine 275 in shaping the binding pocket of the ATP-gated P2X3 receptor

ATP-activated P2X3 receptors expressed in nociceptive sensory neurons play an important role in pain signaling. Basic properties of this receptor subtype, including very strong desensitization, depend on the rate of dissociation of the agonist from the binding site. Even though the rough structure of the ATP binding site has been proposed on the basis of the X-ray structure of the zebrafish P2X4 receptor and mutagenesis studies, the fine subunit-specific structural properties predisposing the receptor to tight capture of the agonist inside the binding pocket have not been elucidated. In this work, by exploring in silico the functional role for the left flipper located in the ectodomain region, we identified within this loop a candidate residue S275, which could contribute to the closure of the agonist-binding pocket. Testing of the S275 mutants using the patch-clamp technique revealed a crucial role for S275 in agonist binding and receptor desensitization. The S275A mutant showed a reduced rate of onset of desensitization and accelerated resensitization and was weakly inhibited by nanomolar agonist. Extracellular calcium application produced inhibition instead of facilitation of membrane currents. Moreover, some full agonists became only partial agonists when applied to the S275A receptor. These effects were stronger with the more hydrophobic mutants S275C and S275V. Taken together, our data suggest that S275 contributes to the closure of the agonist-binding pocket and that effective capture of the agonist provided by the left flipper in calcium-dependent manner determines the high rate of desensitization, slow recovery, and sensitivity to nanomolar agonist of the P2X3 receptor.

Temperature dependence of acetylcholine receptor channels activated by different agonists

The temperature dependence of agonist binding and channel gating were measured for wild-type adult neuromuscular acetylcholine receptors activated by acetylcholine, carbamylcholine, or choline. With acetylcholine, temperature changed the gating rate constants (Q(10) ≈ 3.2) but had almost no effect on the equilibrium constant. The enthalpy change associated with gating was agonist-dependent, but for all three ligands it was approximately equal to the corresponding free-energy change. The equilibrium dissociation constant of the resting conformation (K(d)), the slope of the rate-equilibrium free-energy relationship (Φ), and the acetylcholine association and dissociation rate constants were approximately temperature-independent. In the mutant αG153S, the choline association and dissociation rate constants were temperature-dependent (Q(10) ≈ 7.4) but K(d) was not. By combining two independent mutations, we were able to compensate for the catalytic effect of temperature on the decay time constant of a synaptic current. At mouse body temperature, the channel-opening and -closing rate constants are ∼400 and 16 ms(-1). We hypothesize that the agonist dependence of the gating enthalpy change is associated with differences in ligand binding, specifically to the open-channel conformation of the protein.

Conformational plasticity of the influenza A M2 transmembrane helix in lipid bilayers under varying pH, drug binding, and membrane thickness

Membrane proteins change their conformations to respond to environmental cues, thus conformational plasticity is important for function. The influenza A M2 protein forms an acid-activated proton channel important for the virus lifecycle. Here we have used solid-state NMR spectroscopy to examine the conformational plasticity of membrane-bound transmembrane domain of M2 (M2TM). (13)C and (15)N chemical shifts indicate coupled conformational changes of several pore-facing residues due to changes in bilayer thickness, drug binding, and pH. The structural changes are attributed to the formation of a well-defined helical kink at G34 in the drug-bound state and in thick lipid bilayers, nonideal backbone conformation of the secondary-gate residue V27 in the presence of drug, and nonideal conformation of the proton-sensing residue H37 at high pH. The chemical shifts constrained the (ϕ, ψ) torsion angles for three "basis" states, the equilibrium among which explains the multiple resonances per site in the NMR spectra under different combinations of bilayer thickness, drug binding, and pH conditions. Thus, conformational plasticity is important for the proton conduction and inhibition of M2TM. The study illustrates the utility of NMR chemical shifts for probing the structural plasticity and folding of membrane proteins.

Molecular determinants of coupling between the domain III voltage sensor and pore of a sodium channel

In a voltage-dependent sodium channel, the activation of voltage sensors upon depolarization leads to the opening of the pore gates. To elucidate the principles underlying this conformational coupling, we investigated a putative gating interface in domain III of the sodium channel using voltage-clamp fluorimetry and tryptophan-scanning mutagenesis. Most mutations have similar energetic effects on voltage-sensor activation and pore opening. However, several mutations stabilized the activated voltage sensor while concurrently destabilizing the open pore. When mapped onto a homology model of the sodium channel, most localized to hinge regions of the gating interface. Our analysis shows that these residues are involved in energetic coupling of the voltage sensor to the pore when both are in resting and when both are in activated conformations, supporting the notion that electromechanical coupling in a voltage-dependent ion channel involves the movement of rigid segments connected by elastic hinges.

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.

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.

Toolbox for analyzing finite two-state trajectories

In many experiments, the aim is to deduce an underlying multisubstate on-off kinetic scheme (KS) from the statistical properties of a two-state trajectory. However, a two-state trajectory that is generated from an on-off KS contains only partial information about the KS, and so, in many cases, more than one KS can be associated with the data. We recently showed that the optimal way to solve this problem is to use canonical forms of reduced dimensions (RDs). RD forms are on-off networks with connections only between substates of different states, where the connections can have nonexponential waiting time probability density functions (WT-PDFs). In theory, only a single RD form can be associated with the data. To utilize RD forms in the analysis of the data, a RD form should be associated with the data. Here, we give a toolbox for building a RD form from a finite time, noiseless, two-state trajectory. The methods in the toolbox are based on known statistical methods in data analysis, combined with statistical methods and numerical algorithms designed specifically for the current problem. Our toolbox is self-contained-it builds a mechanism based only on the information it extracts from the data, and its implementation is fast (analyzing a 10;{6}cycle trajectory from a 30-parameter mechanism takes a couple of hours on a PC with a 2.66GHz processor). The toolbox is automated and is freely available for academic research upon electronic request.

A residue at the cytoplasmic entrance of BK-type channels regulating single-channel opening by its hydrophobicity

Single large-conductance calcium-activated K(+) (BK) channels encoded by the mSlo gene usually have synchronous gating, but a Drosophila dSlo (A2/C2/E2/G5/10) splice variant (dSlo1A) exhibits very flickery openings. To probe this difference in gating, we constructed a mutant I323T. This channel exhibits four subconductance levels similar to those of dSlo1A. Rectification of the single-channel current-voltage relation of I323T decreased as [Ca(2+) ](in) increased from 10 to 300 microM. Mutagenesis suggests that the hydrophobicity of the residue at the position is important for the wild-type gating; i.e., increasing hydrophobicity prolongs open duration. Molecular dynamics simulation suggests that four hydrophobic pore-lining residues at position 323 of mSlo act cooperatively in a "shutter-like" mechanism gating the permeation of K(+) ions. Rate-equilibrium free energy relations analysis shows that the four I323 residues in an mSlo channel have a conformation 65% similar to the closed conformation during gating. Based on these observations, we suggest that the appearance of rectification and substates of BK-type channels arise from a reduction of the cooperativity among these four residues and a lower probability of being open.

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.

Phenylalanine-508 mediates a cytoplasmic-membrane domain contact in the CFTR 3D structure crucial to assembly and channel function

Deletion of phenylalanine-508 (Phe-508) from the N-terminal nucleotide-binding domain (NBD1) of the cystic fibrosis transmembrane conductance regulator (CFTR), a member of the ATP-binding cassette (ABC) transporter family, disrupts both its folding and function and causes most cystic fibrosis. Most mutant nascent chains do not pass quality control in the ER, and those that do remain thermally unstable, only partially functional, and are rapidly endocytosed and degraded. Although the lack of the Phe-508 peptide backbone diminishes the NBD1 folding yield, the absence of the aromatic side chain is primarily responsible for defective CFTR assembly and channel gating. However, the site of interdomain contact by the side chain is unknown as is the high-resolution 3D structure of the complete protein. Here we present a 3D structure of CFTR, constructed by molecular modeling and supported biochemically, in which Phe-508 mediates a tertiary interaction between the surface of NBD1 and a cytoplasmic loop (CL4) in the C-terminal membrane-spanning domain (MSD2). This crucial cytoplasmic membrane interface, which is dynamically involved in regulation of channel gating, explains the known sensitivity of CFTR assembly to many disease-associated mutations in CL4 as well as NBD1 and provides a sharply focused target for small molecules to treat CF. In addition to identifying a key intramolecular site to be repaired therapeutically, our findings advance understanding of CFTR structure and function and provide a platform for focused biochemical studies of other features of this unique ABC ion channel.

Properties of the generalized master equation: Green's functions and probability density functions in the path representation

The Green's function for the master equation and the generalized master equation in path representation is an infinite sum over the length of path probability density functions (PDFs). In this paper, the properties of path PDFs are studied both qualitatively and quantitatively. The results are used in building efficient approximations for Green's function in 1D, and are relevant in modeling and in data analysis.

Solid-state NMR characterization of conformational plasticity within the transmembrane domain of the influenza A M2 proton channel

Membrane protein function within the membrane interstices is achieved by mechanisms that are not typically available to water-soluble proteins. The whole balance of molecular interactions that stabilize three-dimensional structure in the membrane environment is different from that in an aqueous environment. As a result interhelical interactions are often dominated by non-specific van der Waals interactions permitting dynamics and conformational heterogeneity in these interfaces. Here, solid-state NMR data of the transmembrane domain of the M2 protein from influenza A virus are used to exemplify such conformational plasticity in a tetrameric helical bundle. Such data lead to very high resolution structural restraints that can identify both subtle and substantial structural differences associated with various states of the protein. Spectra from samples using two different preparation protocols, samples prepared in the presence and absence of amantadine, and spectra as a function of pH are used to illustrate conformational plasticity.

Conformational dynamics of the alphaM3 transmembrane helix during acetylcholine receptor channel gating

Muscle acetylcholine receptors are synaptic ion channels that "gate" between closed- and open-channel conformations. We used Phi-value analysis to probe the transition state of the diliganded gating reaction with regard to residues in the M3, membrane-spanning helix of the muscle acetylcholine receptor alpha-subunit. Phi (a fraction between 1 and 0) parameterizes the extent to which a mutation changes the opening versus the closing rate constant and, for a linear reaction mechanism, the higher the Phi-value, the "earlier" the gating motion. In the upper half of alphaM3 the gating motions of all five tested residues were temporally correlated (Phi approximately 0.30) and serve to link structural changes occurring at the middle of the M2, pore-lining helix with those occurring at the interface of the extracellular and transmembrane domains. alphaM3 belongs to a complex and diverse set of synchronously moving parts that change structure relatively late in the channel-opening process. The propagation of the gating Brownian conformational cascade has a complex spatial distribution in the transmembrane domain.