 |
The Journal of Neuroscience, March 11, 2009, 29(10):3189-3199; doi:10.1523/JNEUROSCI.6185-08.2009
 Previous Article | Next Article 
Cellular/Molecular
Binding to Gating Transduction in Nicotinic Receptors: Cys-Loop Energetically Couples to Pre-M1 and M2–M3 Regions
Won Yong Lee,1,2
Chris R. Free,1,2 and
Steven M. Sine1,2,3
1Receptor Biology Laboratory and Departments of 2Physiology and Biomedical Engineering and 3Neurology, Mayo Clinic College of Medicine, Rochester, Minnesota 55905
 |
Abstract
|
|---|
The nicotinic acetylcholine receptor (AChR) transduces binding of nerve-released ACh into opening of an intrinsic ion channel, yet the intraprotein interactions behind transduction remain to be fully elucidated. Attention has focused on the region of the AChR in which the β1–β2 and Cys-loops from the extracellular domain project into a cavity framed by residues preceding the first transmembrane domain (pre-M1) and the linker spanning transmembrane domains M2 and M3. Previous studies identified a principal transduction pathway in which the pre-M1 domain is coupled to the M2–M3 linker through the β1–β2 loop. Here we identify a parallel pathway in which the pre-M1 domain is coupled to the M2–M3 linker through the Cys-loop. Mutagenesis, single-channel kinetic analyses and thermodynamic mutant cycle analyses reveal energetic coupling among  Leu 210 from the pre-M1 domain, Phe 135 and Phe 137 from the Cys-loop, and Leu 273 from the M2–M3 linker. Residues at equivalent positions of non--subunits show negligible coupling, indicating these interresidue couplings are specific to residues in the -subunit. Thus, the extracellular β1–β2 and Cys-loops bridge the pre-M1 domain and M2–M3 linker to transduce agonist binding into channel gating.
|
Introduction
|
|---|
Rapid synaptic transmission comprises quantal release of neurotransmitter, followed by activation of postsynaptic receptors. The core function of postsynaptic receptors is to bind the neurotransmitter and then transduce binding into opening of an ion channel several tens of angstroms from the binding site. Atomic-scale insight into the transduction process is now possible because of high-resolution x-ray structures of prokaryotic receptor ancestors (Hilf and Dutzler, 2008 ; Bocquet et al., 2009), the isolated binding domains (Armstrong et al., 1998; Dellisanti et al., 2007), homologs of the binding domains (Brejc et al., 2001; Hansen et al., 2005), and a 4-Å-resolution cryo-electron microscopic structure of the Torpedo AChR (Unwin, 2005). Given residue locations within functionally crucial regions of synaptic receptors, structural, electrophysiological, and computational approaches provide powerful means to identify intraprotein structures that transduce neurotransmitter binding into channel gating.
The AChR from the motor endplate has served as a prototype to understand receptor-coupled ion channels from the Cys-loop receptor superfamily. The AChR contains five homologous subunits arranged as barrel-staves around a central ion-conducting pore. Each subunit contains a large extracellular domain composed primarily of β-sheets, a pore domain composed of four transmembrane -helices, and a cytoplasmic domain composed of partially defined -helical structure (Unwin, 2005). Transduction begins with binding of ACh to each of two sites located at interfaces between extracellular regions of - and either  - or  -subunits. After ACh enters each binding site, a hairpin structure called the C-loop changes from an open to a closed conformation that traps ACh within the site. One β-strand from each C-loop extends to the junction of the extracellular and pore domains where it forms the pre-M1 region that covalently links the two domains. Two loops from the extracellular domain, the β1–β2- and Cys-loops, also project into the interdomain junction (Fig. 1). The β1–β2 loop is a key element of a principal pathway that transduces ACh binding into channel gating, and through a series of interresidue contacts, functionally couples the pre-M1 domain to the linker spanning the M2 and M3 -helices from the pore domain (Lee and Sine, 2005). The Cys-loop is situated adjacent to the β1–β2 loop, where it also bridges the pre-M1 domain and M2–M3 linker (Unwin, 2005). Found in all subunits in the superfamily of Cys-loop receptors, the Cys-loop is essential for AChR assembly (Fu and Sine, 1996; Green and Wanamaker, 1997) and rapid and efficient gating of the AChR channel (Jha et al., 2007).
Here we show that contributions of the Cys-loop to channel gating depend on residues in the pre-M1 domain and the M2–M3 linker. Site-directed mutagenesis, single-channel kinetic analyses, and thermodynamic mutant cycle analyses reveal strong energetic coupling among Phe 135 and Phe 137 from the Cys-loop, Leu 210 from the pre-M1 domain, and Leu 273 from the M2–M3 linker. The overall findings show that the β1–β2 and Cys-loops form functionally parallel transduction pathways linking the pre-M1 domain to the M2–M3 linker.
|
Materials and Methods
|
|---|
Construction of wild-type and mutant AChRs
Human -, β-, -, and -subunit cDNAs were subcloned in the CMV-based mammalian expression vector pRBG4 (Lee et al., 1991) as described (Ohno et al., 1996). Site-directed mutations were made using the QuikChange Site-Directed Mutagenesis Kit (Stratagene). The presence of each mutation and the absence of unwanted mutations were confirmed by sequencing the entire cDNA insert.
Mammalian cell expression
All experiments were performed using the BOSC 23 cell line (CRL-11270, ATCC), a variant of the HEK 293 cell line. Cells were maintained in the DMEM containing FBS (10% v/v) at 37°C until they reached  50–70% confluence. Wild-type or mutant AChR cDNAs plus the cDNA encoding green fluorescent protein were transfected by calcium-phosphate precipitation using final concentrations of each cDNA of 0.68 µg/ml. Patch-clamp recordings were performed 1–2 d after transfection.
Patch-clamp single-channel recordings
To record single-channel currents, cells transfected with wild-type or mutant AChR cDNAs were rinsed with and maintained in the following bath solution (in mM): KCl 142, NaCl 5.4, CaCl2 1.8, MgCl2 1.7, and HEPES 10 (pH was adjusted to 7.4 with NaOH). The same solution was used to fill patch pipettes. ACh (Sigma-Aldrich) was kept as a 10 mM stock in bath solution and stored at –80°C until use. Glass micropipettes (type 7052; Garner Glass) were coated with Sylgard 184 (Dow Corning) and heat polished to yield resistances of 5–8 M . After identifying transfected cells under fluorescence optics, single-channel currents were recorded in the cell-attached configuration using the Axopatch 200A (Molecular Devices) at a holding membrane potential of –70 mV at 22°C. Data were collected from two to four different patches for each ACh concentration, choosing only recordings in which channel activity was low enough to allow identification of activation episodes resulting from a single AChR channel. The current signal was low-pass filtered at 50 kHz and recorded to hard disk at 200 kHz using the program Acquire (Bruxton). Recordings of single-channel currents spanning a range of ACh concentrations from each mutant receptor were typically performed over a period of 12 h.
Single-channel kinetic analysis
Detailed methods of single-channel kinetic analyses were described previously (Lee and Sine, 2004). Briefly, the digitized current signal was filtered using a 10 kHz digital Gaussian filter (Colquhoun and Sigworth, 1983). Channel events were detected by the half-amplitude threshold criterion using the program TAC (Bruxton) with an imposed dead time of 10 µs. Precise determination of each dwell time at threshold was achieved by cubic spline interpolation of the digital signal, and the measured dwell time at threshold was corrected for the effects of the Gaussian filter (Colquhoun and Sigworth, 1983). Open and closed time histograms were fitted by the sum of exponentials by maximum likelihood using the program TACFit (Bruxton). Openings corresponding to a single receptor channel were identified by assigning a critical closed time defined as the point of intersection of the closed time component that depended on ACh concentration with the succeeding concentration-independent component (Ohno et al., 1996), presumably due to fast desensitization. Because this method removes most but not all closed dwell times due to fast desensitization, analyses were also performed in which the critical closed time was defined as the point of intersection of the closed time component due to fast desensitization with the succeeding component, and the scheme was modified to include a desensitized state connected to the doubly occupied open state. To obtain kinetically homogeneous data, defined single-channel episodes containing at least five openings were analyzed for channel open probability, mean open time, and mean closed time; episodes within two SDs of the mean were accepted for further analysis (Wang et al., 1997). This processing removes a small population of single brief openings flanked by long closed periods (Sine and Steinbach, 1987). For receptors examined in this study, the overall selection retained >80% of the detected single-channel events. The kinetic scheme was fitted simultaneously to all the data obtained over a range of ACh concentrations using MIL software (QuB suite, State University of New York, Buffalo), which employs a maximum likelihood optimization criterion, corrects for missed events, and computes error estimates of the fitted parameters (Qin et al., 1996). An instrument dead time of 22 µs was uniformly applied to all recordings. For a given wild-type or mutant receptor, our analysis comprised data from two to four patches for each ACh concentration, with the ACh concentrations spaced at half log-unit intervals over a two to three log-unit range of concentration. For each receptor mutant, the range of ACh concentrations spanned from the minimum to the maximum open channel probability.
The following kinetic scheme was fitted to the single-channel open and closed dwell times:  Two ACh (A) molecules associate with the receptor (R) with rate constants k+1 and k+2, and dissociate with rate constants k–1 and k–2. Closed but occupied receptors (AR and A2R) produce singly (AR*) and doubly (A2R*) occupied open states that close with rate constants 1 and 2, respectively. Rate constants for open channel blocking (k+b) and unblocking (k–b) by ACh are included.
Double mutant cycle analysis
To detect interresidue energetic coupling and estimate its magnitude, we applied mutant cycle analysis (Carter et al., 1984) to the diliganded gating equilibrium constant  2, computed from the ratio of rate constants β2/2 (Lee and Sine, 2005). The following mutant cycle examines energetic coupling between two residues, 1 and 2, for which substitutions alter a functional parameter such as the channel gating equilibrium constant examined in the current study.  For each mutation (mt), the change in 2 relative to that of the wild-type (wt) AChR is expressed as  G = –RTln(2mt/2wt), where R is the gas constant and T the absolute temperature. The free energy change along either pathway from wt1–wt2 to mt1–mt2 is equal, expressed as G1 + G2' = G1' + G2. Rearranging gives G1 – G1' = G2 – G2', which gives the first-order coupling free energy for the mutant cycle, Gint. If the two residues are independent, G1 = G1' and G2 = G2', whereas if they are interdependent, G1  G1' and G2 G2'.
|
Results
|
|---|
Key residues from the Cys-loop, pre-M1 domain, and M2–M3 linker
At the border of the AChR extracellular and pore domains, an assembly of alternating aromatic and aliphatic residues lodges between the pre-M1 domain and the M2–M3 linker (Fig. 1A) (PDB code 2bg9) (Unwin, 2005). This assembly comprises Leu 210 from the pre-M1 domain, Phe 135 and Phe 137 from the Cys-loop, and Leu 273 from the M2–M3 linker (Fig. 1B). Although residue side chains are depicted with atomic precision, their precise placement is tentative because the resolution of the structure is 4 Å. Among subunits of the Cys-loop receptor superfamily, residues at positions equivalent to these are not conserved, but are hydrophobic or contain hydrophobic elements (Fig. 2). Furthermore, the four residues are conserved among -subunits that form heteromeric AChRs (1–6). Thus, location of the residues between the extracellular and pore domains and their sequence conservation, hydrophobicity, and apparent spatial continuity suggest that they form a functionally coupled network that contributes to the rapid and efficient channel gating characteristic of the AChR.
Contributions of the four residues to channel gating
Previous structure–function studies showed that Phe 135, Phe 137, Leu 210, and Leu 273 contribute to channel gating (Chakrapani et al., 2004; Jha et al., 2007; Purohit and Auerbach, 2007), but how the contributions depend on surrounding residues is not known. Furthermore, the contributions of Phe 137 and Leu 210 to channel gating have not been examined with ACh as the agonist. To provide a foundation to investigate functional interdependence of these residues, we constructed mutations of each residue, transfected BOSC 23 cells with cDNAs encoding each mutant -subunit, together with wild-type β-, -, and -subunit cDNAs, and recorded single-channel currents elicited by ACh (concentration range 3 µM to 1 mM). To identify mutations optimal for studying interresidue energetic coupling, we engineered different side chains at each position, recorded single-channel currents activated by low and high concentrations of ACh, and chose mutations that produced large and readily measurable changes in single-channel closed and open dwell times (supplemental Fig. S1, available at www.jneurosci.org as supplemental material). Recordings of single-channel currents show qualitatively that mutation of Phe 135 to Leu decreases, whereas mutations of Phe 137 to Leu, Leu 210 to Gln, and Leu 273 to Phe increase the probability of channel opening (Fig. 3). Thus, we chose these mutations for further study.
View larger version (31K):
[in this window]
[in a new window]
|
Figure 3. Single-channel currents and dwell time histograms from AChRs containing the indicated mutations. Currents elicited by 30 µM ACh are shown at a bandwidth of 10 kHz, with channel openings as upward deflections; a submaximal concentration of ACh is chosen for illustration to best highlight the effects of the mutations on open channel probability. Corresponding histograms of dwell times within identified clusters of single-channel events are shown on logarithmic time axes with overlaid probability density functions generated from fitting a kinetic scheme simultaneously to data obtained over the entire range of ACh concentrations (see Materials and Methods; fitted rate constants are given in Table 1).
|
|
Changes in the probability of channel opening can arise from changes in rate constants underlying agonist binding or channel gating. To quantify the effect of each mutation on channel gating rate constants, we analyzed sequences of open and closed dwell times using a kinetic scheme in which two agonists bind sequentially followed by gating of the channel (see Materials and Methods). For each recording, long closed times corresponding to dwells in desensitized states were systematically removed, as described previously (Lee and Sine, 2004), and the kinetic scheme was fitted simultaneously to dwell time sequences elicited by the full range of ACh concentrations using the maximum likelihood optimization criterion (Qin et al., 1996). For both wild-type and mutant AChRs, probability density functions computed from the fitted rate constants superimpose well on histograms of closed and open dwell times (Fig. 3). Further, the fitted rate constants for channel opening and closing show that the mutation F135L decreases the channel gating equilibrium constant by 30-fold, while the mutations F137L, L210Q, and L273F increase the gating equilibrium constant 1.5-fold to fourfold (Table 1). Rate constants underlying agonist binding are also affected, with the greatest changes occurring for mutations that enhance channel gating.
Energetic coupling between Phe 135, Leu 210, and Leu 273
The cryo-electron microscopic structure of the Torpedo AChR shows that Leu 210 and Leu 273 are predicted to bracket Phe 135 (Fig. 1), suggesting these residues contribute to channel gating in an interdependent manner. To determine whether the residues' contributions are interdependent, we mutated each residue individually and in pairs, determined the channel gating equilibrium constant for each mutant AChR, and cast the results as two-dimensional mutant cycles (Fig. 4). The mutation F135L decreases the channel gating equilibrium constant by a free energy change of 1.9 kcal/mol relative to that of the wild-type AChR. However, if the mutation L210Q is present in the same AChR, F135L has a negligible effect on channel gating. Conversely, the mutation L210Q increases the channel gating equilibrium constant by a free energy change of –0.7 kcal/mol, but when F135L is present in the same AChR, L210Q enhances channel gating by a free energy change of –2.6 kcal/mol. The first-order coupling free energy, given by the difference between free energy changes for any two parallel lines of the F135L/L210Q mutant cycle, is large, –1.9 kcal/mol; the negative sign indicates the double mutant enhances channel gating relative to that expected from additive contributions of the two single mutations. Analogously, the mutant cycle comprising the mutations F135L and L273F yields a first-order coupling free energy of –1.1 kcal/mol. Thus, the contribution of Phe 135 to channel gating depends on both Leu 210 and Leu 273.
To determine whether equivalent residues in the - and -subunits contribute to channel gating, and further, to determine whether any detected contributions are interdependent, we mutated each residue individually and in pairs, measured channel gating equilibrium constants for each mutant AChR, and generated two-dimensional mutant cycles (Fig. 4). Mutations of residues in the - and -subunits equivalent to Phe 135 decrease channel gating free energy up to 1.5 kcal/mol, but mutations of residues equivalent to Leu 210 or Leu 273 alter gating free energy by 0.4 kcal/mol or less (Table 2). Moreover, the set of two-dimensional mutant cycles reveals negligible first-order coupling free energies (Fig. 4, Table 3), indicating residues at equivalent positions of the - and -subunits are functionally independent. Mutations of equivalent residues in the β-subunit have only minor effects on channel gating free energy (Table 2), so mutant cycle analyses were not performed for these residues. Thus, energetic coupling among Phe 135, Leu 210, and Leu 273 is unique to the -subunit.
Context dependence of pairwise interresidue coupling
Interresidue coupling has been widely observed to be context dependent (Zhang et al., 1996; Albeck et al., 2000), particularly for residues within 7 Å of each other (Schreiber and Fersht, 1995). Thus, to determine whether energetic coupling between any two residues of the triad Phe 135, Leu 210, and Leu 273 depends on the third residue, we generated a cubic mutant cycle (Fig. 5A). The left vertical plane, F135L/L210Q, exhibits a first-order coupling free energy of –1.9 kcal/mol, but the opposing plane, in which the mutation L273F is present in all four mutant AChRs, exhibits a decreased coupling free energy of –0.67 kcal/mol. Analogously, the top plane, F135L/L273F, exhibits a coupling free energy of 1.1 kcal/mol, but the opposing plane, in which the mutation L210Q is present, exhibits a decreased coupling free energy of 0.16 kcal/mol. Finally, the front plane, L210Q/L273F, exhibits a coupling free energy of 0.47 kcal/mol, but the opposing plane, in which the mutation F135L is present, exhibits an increased coupling energy of 1.7 kcal/mol. For the overall cubic cycle, the second-order coupling free energy, given by the difference in first-order coupling free energies between any pair of parallel planes, is 1.2 kcal/mol. Thus, Phe 135, Leu 210, and Leu 273 contribute to channel gating as an interdependent triad.
View larger version (38K):
[in this window]
[in a new window]
|
Figure 5. Energetic coupling between Phe 135, Phe 137, Leu 210, and Leu 273. For a cubic mutant cycle, each plane depicts the first-order coupling free energy (Gint) for the indicated residue pair, as in Figure 4. The free energy difference between any pair of parallel planes gives the second-order coupling free energy (Gint) for the residue triad. For each mutant AChR, single-channel currents elicited by 300 µM ACh are shown at a bandwidth of 10 kHz. A, Energetic coupling of Phe 135, Leu 210, and Leu 273. For the front plane, Gint = 0.47 kcal/mol, and for the back plane, Gint = 1.7 kcal/mol. For this cubic cycle, Gint = 1.2 kcal/mol. B, Energetic coupling of Phe 135, Phe 137, and Leu 210. For the front plane, Gint = 0.11 kcal/mol, and for the back plane, Gint = 1.2 kcal/mol. For this cubic cycle, Gint = 1.1 kcal/mol. Insets show residues examined in each cubic cycle.
|
|
Phe 137 is the fourth member of a coupled tetrad
To determine whether Phe 137 is energetically coupled to residues of the triad just identified, we generated a cubic mutant cycle composed of mutations of Phe 137, Phe 135, and Leu 210 (Fig. 5B). The front plane, F137L/L210Q, exhibits a first-order coupling free energy of 0.11 kcal/mol, but the opposing plane in which the mutation F135L is present in all four mutant AChRs, exhibits an increased coupling free energy of 1.2 kcal/mol. A possible structural interpretation is that coupling between F135 and L210 is so strong that mutation of Phe 137 does not affect the contribution of Leu 210 to channel gating. However if Leu is substituted for Phe 135, Phe 137 remains the nearest aromatic neighbor for Leu 210, and coupling shifts to this residue pair. The left vertical plane, F135L/F137L, exhibits a first-order coupling free energy of –1.1 kcal/mol, but the opposing plane in which the mutation L210Q is present, exhibits a negligible coupling energy of 0.01 kcal/mol. Finally, the top plane, F135L/L210Q, exhibits a coupling energy of –1.9 kcal/mol, but the opposing plane, in which the mutation F137L is present, exhibits a decreased coupling energy of –0.77 kcal/mol. For the overall cubic cycle, the second-order coupling free energy of 1.1 kcal/mol indicates that Phe 137, Leu 210, Phe 135 contribute to channel gating as a coupled triad. Because both Leu 210 and Phe 135 couple strongly to Leu 273, we infer that Phe 137 constitutes the fourth member of a coupled tetrad.
Context dependence of residues in the - and -subunits
Although only weak energetic coupling is observed among pairs of residues in the - and -subunits equivalent to Phe 135, Leu 210, and Leu 273, our analyses of cubic mutant cycles raise the possibility that a neighboring residue may mask interresidue coupling. We therefore generated cubic mutant cycles to examine this possibility, focusing on residues in the - and -subunits at positions equivalent to Phe 135, Leu 210, and Leu 273. For the cubic mutant cycles generated from the equivalent residue triads, second-order coupling free energies are low, –0.28 for the triad in the -subunit and –0.13 kcal/mol for the triad in the -subunit (supplemental Fig. S2, available at www.jneurosci.org as supplemental material). Thus, energetic coupling between residues in the Cys-loop, pre-M1 domain, and M2–M3 linker is specific to the -subunit.
Tyr 277 does not couple to residues in the tetrad
Tyr 277 is predicted to project from the distal end of the M2–M3 linker toward the coupled tetrad (Fig. 1B), and is conserved as an aromatic residue at equivalent positions in all subunits of the Cys-loop receptor superfamily (Fig. 2). Furthermore, mutation of Tyr 277 to Leu decreases free energy of channel gating by 1.1 kcal/mol (Table 1). Therefore, to examine energetic coupling among Tyr 277, Phe 135, Phe 137, and Leu 210, we generated two cubic mutant cycles. In the first cubic cycle (Fig. 6A), the front plane, Y277L/L210Q, exhibits a first-order coupling free energy of –0.9 kcal/mol, while the opposing plane, containing the mutation F135L, exhibits a coupling energy of –0.6 kcal/mol. In the second cubic cycle (Fig. 6B), the front plane, Y277L/F137L, exhibits a coupling free energy of –0.35 kcal/mol, while the opposing plane, containing the mutation F135L, maintains a low coupling energy of 0.16 kcal/mol. For the overall cubic cycles in Figure 6, A and B, the second-order coupling free energies are low, 0.33 and 0.51 kcal/mol, respectively. Thus, although Tyr 277 contributes to channel gating, its contribution depends only weakly on Leu 210, Phe 135, and Phe 137.
View larger version (31K):
[in this window]
[in a new window]
|
Figure 6. Energetic coupling between Phe 135, Phe 137, Leu 210, and Tyr 277. A, Energetic coupling of Phe 135, Leu 210, and Tyr 277. For the front plane, Gint = –0.89 kcal/mol, and for the back plane, Gint = –0.56 kcal/mol. For this cubic cycle, Gint = –0.33 kcal/mol. B, Energetic coupling of Phe 135, Phe 137, and Tyr 277. For the front plane, Gint = –0.35 kcal/mol, and for the back plane, Gint = 0.16 kcal/mol. For this cubic cycle, Gint = –0.51 kcal/mol. Insets show residues examined in each cubic cycle.
|
|
Functional separation of Cys-loop and principal coupling pathways
The preceding findings show that the Cys-loop is analogous to the β1–β2 loop in functionally linking the pre-M1 domain to the M2–M3 linker. To determine whether the Cys-loop and principal coupling pathways overlap functionally, we asked whether Ile 274, which is predicted to contact residues from both pathways, couples functionally to residues within either pathway. In particular, Ile 274 is one of four apparently interlocked aliphatic residues that physically join the Cys-loop to the M2–M3 linker (Fig. 1C), and among Cys-loop receptors, is conserved as Ile or Met (Fig. 2). Furthermore, the mutation I274A decreases free energy of channel gating by 2.5 kcal/mol, primarily by decreasing the rate constant for channel opening (Table 1); of the mutations examined in this study, I274A affects channel gating to the greatest degree, in broad agreement with previous studies (Jha et al., 2007).
To determine whether the functional contribution of Ile 274 depends on residues within the Cys-loop coupling pathway, we generated a cubic mutant cycle comprising mutations of Ile 274 and the two residues in the coupled tetrad predicted to be most proximal, Phe 135 and Leu 273 (Fig. 7A). The top and left planes, I274A/F135L and I274A/L273F, exhibit first-order coupling free energies of 0.1 and 0.02 kcal/mol, respectively, while the opposing planes, containing a mutation of the third residue, maintain low coupling energies of 0.06 and –0.02 kcal/mol. The front plane, F135L/L273F, exhibits a coupling free energy of –1.06 kcal/mol, while the opposing plane, containing the mutation I274A, exhibits a coupling energy of –1.1 kcal/mol. For the overall cubic cycle, the second-order coupling free energy is –0.04 kcal/mol, showing that although Ile 274 contributes strongly to channel gating, its contribution does not depend on Phe 135 and Leu 273 within the Cys-loop coupling pathway.
View larger version (34K):
[in this window]
[in a new window]
|
Figure 7. Energetic coupling between Ile 274 and residues in the Cys-loop or principal coupling pathways. A, Energetic coupling of Phe 135, Leu 210, and Ile 274. For the front plane, Gint = –1.06 kcal/mol, and for the back plane, Gint = –1.1 kcal/mol. For this cubic cycle, Gint = –0.04 kcal/mol. B, Energetic coupling of Pro 272, Val 132, and Ile 274. For the front plane, Gint = –0.01 kcal/mol, and for the back plane, Gint = 1.2 kcal/mol. For this cubic cycle, Gint = –1.2 kcal/mol. Insets show residues examined in each cubic cycle.
|
|
To determine whether the functional contribution of Ile 274 depends on residues within the principal coupling pathway, we generated a cubic mutant cycle including mutations of Val 132 and Pro 272, which contribute to the principal coupling pathway (Lee and Sine, 2005; Lee et al., 2008) (Fig. 7B). The front plane, I274A/V132A, exhibits a negligible first-order coupling free energy, but the opposing plane, containing the mutation P272A, exhibits a coupling energy of 1.2 kcal/mol. Conversely, the left plane, I274A/P272A, exhibits a large first-order coupling energy of –1.8 kcal/mol, but the opposing plane, containing the mutation V132A, exhibits a reduced coupling energy of –0.67 kcal/mol. For the overall cubic cycle, the second-order coupling free energy of –1.2 kcal/mol shows that the strong contribution of Ile 274 to channel gating depends on Val 132 and Pro 272 of the principal coupling pathway.
|
Discussion
|
|---|
The present findings reveal pairwise energetic coupling among residues of the Cys-loop and those in the neighboring pre-M1 and M2–M3 regions, and show that these interresidue couplings are essential for rapid and efficient gating of the AChR channel. The emerging network of coupled residues is present in all -subunits that form heteromeric AChRs (1–6), suggesting that it is generally important in transducing agonist binding into channel gating. Energetic coupling is observed only among residues in the -subunits, analogous to interresidue coupling in the previously identified principal coupling pathway (Lee and Sine, 2005). The contribution of the Cys-loop to transduction is analogous to that of the β1–β2 loop of the principal coupling pathway in that both loops link the pre-M1 to the M2–M3 region. Ileu 274 in the M2–M3 linker does not couple energetically to residues in the Cys-loop pathway, but it couples strongly to residues in the principal pathway, suggesting that the two pathways are functionally separate. The strongest interresidue coupling free energies are observed between residues predicted to be nearest neighbors based on the cryo-electron microscopic Torpedo AChR structure, and as widely observed in proteins, the coupling energies are context dependent. Moreover, in light of recent high-resolution structures of bacterial Cys-loop receptors in candidate closed and open conformations (Bocquet et al., 2009; Hilf and Dutzler, 2009), our findings suggest the principal and Cys-loop pathways contribute jointly in transducing agonist binding to gating of the AChR channel.
Our experimental strategy centers on the use of double mutant cycle analyses to assess interresidue interactions (Horovitz and Fersht, 1990; Horovitz et al., 1990). Mutant cycle analysis is founded on the principle that mutation of a residue causes a free energy change, G1, which depends on other residues within the protein. If mutation of a second residue alters the effect of mutation of the first residue, yielding a different free energy change, G2, the two residues are deemed to interact with a first-order coupling free energy, G = G2 – G1. Interresidue coupling free energy is a thermodynamic parameter, and does not disclose whether coupling arises through a direct or a propagated interaction. However, if the two residues are known to make contact, the most likely interpretation is that in the context of the surrounding structure, the coupling arises through direct contact. Because resolution of the Torpedo AChR is not sufficient to unequivocally establish direct residue contact, the interresidue coupling free energies could have arisen either through direct contact or through interresidue interactions within an 7 Å radius. However, it is worth noting that the interresidue couplings observed here support residue placement in the cryo-electron microscopic structure of the Torpedo AChR (Unwin, 2005).
The cryo-electron microscopic structure of the Torpedo AChR nevertheless establishes that the Cys-loop inserts between the pre-M1 domain and the M2–M3 loop. Residues stemming from these structures (Leu 210, Phe 137, Phe 135, and Leu 273) are predicted to project into the space underneath the Cys-loop and face the pore, forming a row of alternating aliphatic and aromatic side chains approximately parallel to the plane of the membrane. The high-resolution structures of the bacterial homologs ELIC and GLIC confirm the overall placement of the Cys-loop in the lower resolution Torpedo AChR structure, but the detailed shape, orientation relative to the membrane, and pairwise residue contacts of the Cys-loop differ between Torpedo AChR and bacterial counterparts. Nevertheless, in both Torpedo and bacterial structures, the C-terminal half of the Cys-loop associates with the pre-M1 domain, the N-terminal half associates with the M2–M3 linker, and the prevalent interresidue contacts are hydrophobic. In both the inactive ELIC and active GLIC structures, the Cys-loop maintains close association with surrounding residues (Bocquet et al., 2009; Hilf and Dutzler, 2009), suggesting that the interresidue couplings detected here originate from interresidue contacts that are maintained between closed and open states, rather than from contacts that form or break upon transition between the two states.
Throughout this work, the kinetics of AChR activation was analyzed using a kinetic model in which two agonists bind sequentially, the channel gates, and ACh blocks the open channel (Materials and Methods). A recent study, however, provided evidence for a transient closed state, called flipped, that followed agonist binding but preceded channel opening (Lape et al., 2008). The flipped state was suggested to be a major determinant of agonist efficacy, forming with much greater efficiency for full agonists compared with partial agonists. It is possible that mutations examined in this study, rather than altering elementary steps underlying agonist binding or channel gating, altered elementary steps underlying flipping. However, introducing a flipped state into the kinetic model increases the number of free parameters and thus substantially increases the difficulty of obtaining well defined rate constants. In recognition of this potential limitation, we emphasize that the rate constants underlying agonist binding and channel gating presented here should be considered apparent rate constants.
The chief coupling free energies identified herein correspond to the residue pairs Phe 135/Leu 210, Phe135/Leu273, and Phe135/Phe137. The coupling free energies likely arise from optimal tertiary packing of side chains that buries large hydrophobic surfaces. In model peptides that can form -helices, pairs of aromatic or pairs of aromatic and Leu residues stabilize the helix when placed at positions i and i + 4 of the peptide chain, in contrast to lack of stabilization by pairs of aliphatic residues; the stabilization was suggested to arise from decreased loss of conformational entropy for aromatic–aliphatic compared with aliphatic–aliphatic pairs (Padmanabhan and Baldwin, 1994; Creamer and Rose, 1995). The cubic cycle F135L/L210Q/L273F exhibits a second-order coupling free energy of 1.2 kcal/mol and provides evidence that the interresidue triad is energetically coupled. This coupling energy is about half of that determined for residue triads from the principal coupling pathway, V46A/P272G/S269L (–2.0 kcal/mol) (Lee and Sine, 2005) and V132A/V46A/P272A (2.3 kcal/mol) (Lee et al., 2008). Tertiary packing of hydrophobic side chains is likely a major source of interresidue energetic coupling in both the Cys-loop and principal transduction pathways.
Conformational rearrangements underlying transduction of binding to gating have been detected in two distinct regions of the representative Cys-loop receptor: the agonist binding site and the structural transition zone between extracellular and pore domains. The collective evidence suggests a structural mechanism of the transduction process, described as follows.
In the absence of agonist, the peripheral C-loop tilts away from the binding site, maintaining an opened up conformation that allows small molecules to enter the site. Supporting evidence comes from x-ray structures of the homologous AChBP (Celie et al., 2004; Hansen et al., 2005), the bacterial channel ELIC (Hilf and Dutzler, 2008), and the cryo-electron microscopic structure of the Torpedo AChR (Unwin, 2005), all in the absence of activators. After the activator binds, the C-loop changes from an uncapped to a capped conformation, restricting exit or entry of small molecules, again as demonstrated by x-ray structures of AChBP and the bacterial channel GLIC in the presence of activators. Notably, the only comparison of liganded and unliganded structures for the same protein was achieved for AChBP from Aplysia (Hansen et al., 2005). In solution, activator-induced changes of the C-loop have been observed for Lymnaea AChBP through measurements of intrinsic tryptophan fluorescence (Gao et al., 2005), changes in chemical shift by NMR (Gao et al., 2006), and molecular dynamics simulations (Gao et al., 2005; Cheng et al., 2006b). Thus, the accumulated evidence indicates that the initial trigger for channel gating is a capping motion of the C-loop (Mukhtasimova et al., 2005).
Attention has also focused on the junction between the extracellular and pore domains, because this region features an abrupt transition where β-sheets meet -helices. Within this transition zone of each AChR subunit, five segments from different parts of the linear sequence merge: three loops from the extracellular domain, the covalent link between the extracellular and pore domains, and the linker joining the M2 and M3 -helices. A wealth of structure–function studies show that residues from all five intrazone segments are crucial for coupling agonist binding to channel gating (for review, see Sine and Engel, 2006). In particular, a highly conserved salt bridge within the hydrophobic core of the -subunits joins, both physically and functionally, the pre-M1 domain to the tip of the β1–β2 loop (Lee and Sine, 2005; Unwin, 2005). The β1–β2 loop, in turn, forms a hydrophobic pin-in-socket connection to the M2–M3 linker at the top of the pore (Unwin, 2005). This triad of interlinked segments, called the principal coupling pathway, was originally identified in light of the 4 Å resolution cryo-electron microscopic structure of the Torpedo AChR. However the positioning of key elements of the principal pathway, such as the buried salt bridge formed by Arg 209 and Glu 45, have been verified by high-resolution structures of the isolated -subunit extracellular domain (Dellisanti et al., 2007) and the bacterial channel GLIC (Hilf and Dutzler, 2008) and by structure–function studies of the homomeric  1 GABA receptor (Wang et al., 2007).
Comparison of x-ray structures of ELIC in the inactive state and GLIC in the active state suggest that agonist binding causes translation of the β1–β2 loop against the M2–M3 linker, displacing the tops of the M2 and M3 -helices away from an axis through the center of the channel, thus initiating ion conduction (Bocquet et al., 2009; Hilf and Dutzler, 2009). Although the β-strands stemming from the β1–β2 loop do not connect directly to the agonist binding site, interstrand hydrogen bonds hold these strands within the overall β-barrel structure, which also encompasses elements that form the principal face of the agonist binding site. Thus, when the C-loop caps the entrance to the binding site, a quaternary twist of the overall β-barrel (Taly et al., 2006; Cheng et al., 2006a) could cause translation of the β1–β2 loop against the M2–M3 linker and consequently open the channel.
If the link between the pre-M1, β1–β2, and M2–M3 regions is a principal trigger of channel gating, how does the Cys-loop contribute to gating? Structures of the Torpedo AChR and bacterial channels show that the β1–β2 loop overlies the M2–M3 linker and is proximal to an axis through the center of the channel, whereas the larger Cys-loop inserts between the M2–M3 linker and the pre-M1 domain, and is farther from the central axis. The four β-strands that emanate from the β1–β2 and Cys-loops are hydrogen bonded within the overall β-barrel structure, so a quaternary twist of the barrel would likely move both loops in concert. The two loops impinge on opposite sides and at opposite ends of the M2–M3 linker, which extends radially from the central axis. Studies of linear free energy relationships of channel gating reveal similar slopes of the plot of channel opening rate against channel opening equilibrium constant for the β1–β2 and Cys-loops (Jha et al., 2007), suggesting the two loops change conformation at the same time. A concerted movement of the two loops would displace the linker away from the central axis, moving the tops of the M2 and M3 -helices in the same direction, dilating the top of the pore (Bocquet et al., 2009; Hilf and Dutzler, 2009). Thus, through specific interresidue couplings, the β1–β2 and Cys-loops loops appear to act jointly on M2–M3 linker to open the pore for ion conduction. Upon relaxation of the twisted β-barrel, return of the β1–β2 and Cys-loops to their original positions would constrict the top of the pore, closing the channel.
|
Footnotes
|
|---|
Received Dec. 30, 2008;
revised Feb. 5, 2009;
accepted Feb. 6, 2009.
W. Y. Lee's present address: Department of Biology, University of Utah, Salt Lake City, UT 84112.
This work was supported by National Institutes of Health Grant NS031744 (S.M.S.).
Correspondence should be addressed to Steven M. Sine at the above address. Email: sine{at}mayo.edu
Copyright © 2009 Society for Neuroscience 0270-6474/09/293189-11$15.00/0
|
References
|
|---|
Albeck S, Unger R, Schreiber G (2000) Evaluation of direct and cooperative contributions towards the strength of buried hydrogen bonds and salt bridges. J Mol Biol 298:503–520.[CrossRef][Web of Science][Medline] Armstrong N, Sun Y, Chen GQ, Gouaux E (1998) Structure of a glutamate-receptor ligand-binding core in complex with kainate. Nature 395:913–917.[CrossRef][Medline] Bocquet N, Nury H, Baaden M, Le Poupon C, Changeux JP, Delarue M, Corringer PJ (2009) X-ray structure of a pentameric ligand-gated ion channel in an apparently open conformation. Nature 457:111–114.[CrossRef][Web of Science][Medline] Brejc K, van Dijk WJ, Klaassen RV, Schuurmans M, van Der Oost J, Smit AB, Sixma TK (2001) Crystal structure of an ACh-binding protein reveals the ligand-binding domain of nicotinic receptors. Nature 411:269–276.[CrossRef][Medline] Carter PJ, Winter G, Wilkinson AJ, Fersht AR (1984) The use of double mutants to detect structural changes in the active site of the tyrosyl-tRNA synthetase (Bacillus stearothermophilus). Cell 38:835–840.[CrossRef][Web of Science][Medline] Celie PH, van Rossum-Fikkert SE, van Dijk WJ, Brejc K, Smit AB, Sixma TK (2004) Nicotine and carbamylcholine binding to nicotinic acetylcholine receptors as studied in AChBP crystal structures. Neuron 41:907–914.[CrossRef][Web of Science][Medline] Chakrapani S, Bailey TD, Auerbach A (2004) Gating dynamics of the acetylcholine receptor extracellular domain. J Gen Physiol 123:341–356.[Abstract/Free Full Text] Cheng X, Lu B, Grant B, Law RJ, McCammon JA (2006a) Channel opening motion of alpha7 nicotinic acetylcholine receptor as suggested by normal mode analysis. J Mol Biol 355:310–324.[CrossRef][Web of Science][Medline] Cheng X, Wang H, Grant B, Sine SM, McCammon JA (2006b) Targeted molecular dynamics study of C-loop closure and channel gating in nicotinic receptors. PLoS Comput Biol 2:e134.[CrossRef][Medline] Colquhoun D, Sigworth F (1983) Fitting and statistical analysis of single channel records. In: Single-channel recording (Sakmann B, Neher E, eds), pp 191–264. New York: Plenum. Creamer TP, Rose GD (1995) Interactions between hydrophobic side chains within a-helices. Protein Sci 4:1305–1314.[Web of Science][Medline] Dellisanti CD, Yao Y, Stroud JC, Wang ZZ, Chen L (2007) Crystal structure of the extracellular domain of nAChR 1 bound to -bungarotoxin at 1.94 Å resolution. Nat Neurosci 10:953–962.[CrossRef][Web of Science][Medline] Fu DX, Sine SM (1996) Asymmetric contribution of the conserved disulfide loop to subunit oligomerization and assembly of the nicotinic acetylcholine receptor. J Biol Chem 271:31479–31484.[Abstract/Free Full Text] Gao F, Bren N, Burghardt TP, Hansen S, Henchman RH, Taylor P, McCammon JA, Sine SM (2005) Acetylcholine-mediated conformational changes in acetylcholine-binding protein revealed by simulation and intrinsic tryptophan fluorescence. J Biol Chem 280:8443–8451.[Abstract/Free Full Text] Gao F, Mer G, Tonelli M, Hansen SB, Burghardt TP, Taylor P, Sine SM (2006) Solution NMR of acetylcholine binding protein reveals agonist-mediated conformational change of the C-loop. Mol Pharmacol 70:1230–1235.[Abstract/Free Full Text] Green WN, Wanamaker CP (1997) The role of the cystine loop in acetylcholine receptor assembly. J Biol Chem 272:20945–20953.[Abstract/Free Full Text] Hansen SB, Sulzenbacher G, Huxford T, Marchot P, Taylor P, Bourne Y (2005) Structures of Aplysia AChBP complexes with nicotinic agonists and antagonists reveal distinctive binding interfaces and conformations. EMBO J 24:3635–3646.[CrossRef][Web of Science][Medline] Hilf RJ, Dutzler R (2008) X-ray structure of a prokaryotic pentameric ligand-gated ion channel. Nature 452:375–379.[CrossRef][Medline] Hilf RJ, Dutzler R (2009) Structure of a potentially open state of a proton-activated pentameric ligand-gated ion channel. Nature 457:115–118.[CrossRef][Web of Science][Medline] Horovitz A, Fersht AR (1990) Strategy for analysing the co-operativity of intramolecular interactions in peptides and proteins. J Mol Biol 214:613–617.[CrossRef][Web of Science][Medline] Horovitz A, Serrano L, Avron B, Bycrogt M, Fersht AR (1990) Strength and cooperativity of contributions of surface salt bridges to protein stability. J Mol Biol 216:1031–1044.[Web of Science][Medline] Jha A, Cadugan DJ, Purohit P, Auerbach A (2007) Acetylcholine receptor gating at extracellular transmembrane domain interface: the cys-loop and m2–m3 linker. J Gen Physiol 130:547–558.[Abstract/Free Full Text] Lape R, Colquhoun D, Sivilotti LG (2008) On the nature of partial agonism in the nicotinic receptor superfamily. Nature 454:722–727.[Web of Science][Medline] Lee BS, Gunn RB, Kopito RR (1991) Functional differences among nonerythroid anion exchangers expressed in a transfected human cell line. J Biol Chem 266:11448–11454.[Abstract/Free Full Text] Lee WY, Sine SM (2004) Invariant aspartic acid in muscle nicotinic receptor contributes selectively to the kinetics of agonist binding. J Gen Physiol 124:555–567.[Abstract/Free Full Text] Lee WY, Sine SM (2005) Principal pathway coupling agonist binding to channel gating in nicotinic receptors. Nature 438:243–247.[CrossRef][Medline] Lee WY, Free C, Sine SM (2008) Nicotinic receptor inter-loop proline anchors β1-β2 and Cys-loops in coupling agonist binding to channel gating. J Gen Physiol 132:265–278.[Abstract/Free Full Text] Mukhtasimova N, Free C, Sine SM (2005) Initial coupling of binding to gating mediated by conserved residues in muscle nicotinic receptor. J Gen Physiol 126:23–39.[Abstract/Free Full Text] Ohno K, Wang HL, Milone M, Bren N, Brengman JM, Nakano S, Quiram P, Pruitt JN, Sine SM, Engel AG (1996) Congenital myasthenic syndrome caused by decreased agonist binding affinity due to a mutation in the acetylcholine receptor epsilon subunit. Neuron 17:157–170.[CrossRef][Web of Science][Medline] Padmanabhan S, Baldwin RL (1994) Tests for helix-stabilizing interactions between various nonpolar side chains in alanine-based peptides. Protein Sci 3:1992–1997.[Web of Science][Medline] Purohit P, Auerbach A (2007) Acetylcholine receptor gating at extracellular transmembrane domain interface: the "pre-M1" linker. J Gen Physiol 130:559–568.[Abstract/Free Full Text] Qin F, Auerbach A, Sachs F (1996) Estimating single-channel kinetic parameters from idealized patch-clamp data containing missed events. Biophys J 70:264–280.[Web of Science][Medline] Schreiber G, Fersht AR (1995) Energetics of protein-protein interactions: analysis of the barnase-barstar interface by single mutations and double mutant cycles. J Mol Biol 248:478–486.[Web of Science][Medline] Sine SM, Engel AG (2006) Recent advances in Cys-loop structure and function. Nature 440:448–455.[CrossRef][Medline] Sine SM, Steinbach JH (1987) Activation of acetylcholine receptors on clonal mammalian BC3H-1 cells by high concentrations of agonist. J Physiol 385:325–359.[Abstract/Free Full Text] Taly A, Corringer PJ, Grutter T, Prado de Carvalho L, Karplus M, Changeux JP (2006) Implications of the quaternary twist allosteric model for the physiology and pathology of nicotinic receptors. Proc Natl Acad Sci U S A 103:16965–16970.[Abstract/Free Full Text] Unwin N (2005) Refined structure of the nicotinic acetylcholine receptor at 4 Å resolution. J Mol Biol 346:967–989.[CrossRef][Web of Science][Medline] Wang HL, Auerbach A, Bren N, Ohno K, Engel AG, Sine SM (1997) Mutation in the M1 domain of the acetylcholine receptor alpha subunit decreases the rate of agonist dissociation. J Gen Physiol 109:757–766.[Abstract/Free Full Text] Wang J, Lester HA, Dougherty DA (2007) Establishing an ion pair interaction in homomeric rho1 gamma-aminobutyric acid type A receptor that contributes to the gating pathway. J Biol Chem 282:26210–26216.[Abstract/Free Full Text] Zhang H, Skinner MM, Sandberg WS, Wang AH, Terwilliger TC (1996) Context dependence of mutational effects in a protein: the crystal structures of the V35I, I47V and V35I/I47V gene V protein core mutants. J Mol Biol 259:148–159.[CrossRef][Web of Science][Medline]
Related articles in J. Neurosci.:
- Novel Hypotheses from a Neuropsychological Case Study: Is the Visual Ventral Cortex Critical for Both Category-Generic and Category-Specific Form Perception?
- Johanna C. Goll
J. Neurosci. 2009 29: 11421-11423.
[Full Text]
|
|
Top
Abstract
Introduction
-subunits
