Nicotinic Receptor Fourth Transmembrane Domain: Hydrogen Bonding by Conserved Threonine Contributes to Channel Gating Kinetics

Cecilia Bouzata, Francisco Barrantesa, and Steven Sineb
a From the Instituto de Investigaciones Bioquímicas, Universidad Nacional del Sur-Consejo Nacional de Investigaciones Científicas y Técnicas, 8000 Bahía Blanca, Argentina
b Receptor Biology Laboratory, Department of Physiology and Biophysics, Mayo Foundation, Rochester, Minnesota 55905

Correspondence to: Cecilia Bouzat, Instituto de Investigaciones Bioquímicas, UNS-CONICET, Camino La Carrindanga Km 7, 8000 Bahía Blanca, Argentina. Fax:54-291-4861201 E-mail:inbouzat{at}criba.edu.ar.


  Abstract
Top
Abstract
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
Acknowledgements
References

The fourth transmembrane domain (M4) of the nicotinic acetylcholine receptor (AChR) contributes to the kinetics of activation, yet its close association with the lipid bilayer makes it the outermost of the transmembrane domains. To investigate mechanistic and structural contributions of M4 to AChR activation, we systematically mutated {alpha}T422, a conserved residue that has been labeled by hydrophobic probes, and evaluated changes in rate constants underlying ACh binding and channel gating steps. Aromatic and nonpolar mutations of {alpha}T422 selectively affect the channel gating step, slowing the rate of opening two- to sevenfold, and speeding the rate of closing four- to ninefold. Additionally, kinetic modeling shows a second doubly liganded open state for aromatic and nonpolar mutations. In contrast, serine and asparagine mutations of {alpha}T422 largely preserve the kinetics of the wild-type AChR. Thus, rapid and efficient gating of the AChR channel depends on a hydrogen bond involving the side chain at position 422 of the M4 transmembrane domain.

Key Words: patch clamp, kinetic analysis, nicotinic acetylcholine receptor channel gating, fourth transmembrane domain, hydrogen bond


  INTRODUCTION
Top
Abstract
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
Acknowledgements
References

The nicotinic acetylcholine receptor (AChR)1 is a pentamer of homologous subunits with composition {alpha}2ß{gamma}{delta} in fetal muscle and {alpha}2ß{epsilon}{delta} in the adult. Each subunit contains an amino-terminal extracellular domain of ~210 amino acids, four transmembrane domains (M1–M4), and a short extracellular carboxy-terminal tail. The M2 domain of each subunit contributes to the cation-selective channel, and agonist binding triggers motion of M2 to initiate ion flow (Unwin 1995 ). The locations and functional roles of the M1, M3, and M4 transmembrane domains are not as well understood. The M1 domain appears closely associated with both the ion conducting pathway and the bilayer, as revealed by substituted cysteine accessibility and labeling by hydrophobic reagents (Akabas and Karlin 1995 ; Blanton and Cohen 1994 ). The M3 domain has been shown to have some contact with lipid (Blanton and Cohen 1994 ) and to contribute to channel opening and closing rates (Wang et al. 1999 ). The M4 domain is the least conserved among the transmembrane domains, is the most hydrophobic, and has been extensively labeled by hydrophobic probes (Blanton and Cohen 1992 , Blanton and Cohen 1994 ), suggesting that it is the outermost of the transmembrane domains. Labeling by the photoactivatable hydrophobic probe [3H]TID {[3H]3-(trifluoromethyl)3-m-(iodophenyl)diazirine} occurs at C412, M415, C418, T422, and V425 of the {alpha} subunit of Torpedo californica AChR, making these the best candidates for residues in contact with the lipid bilayer (Blanton and Cohen 1994 ). Mutagenesis studies suggest that the M4 domain is essential for proper AChR activation (Lee et al. 1994 ; Ortiz-Miranda et al. 1997 ). By alanine scanning of the lipid-exposed residues in the M4 domain of the {alpha} subunit, we recently showed that T422 affects channel open duration (Bouzat et al. 1998 ).

Here we examine the structural and mechanistic contributions of the M4 domain to AChR activation by systematically mutating {alpha}T422 of the mouse muscle AChR and evaluating changes in activation kinetics. Our results show that {alpha}T422 contributes through a hydrogen bond to both opening and closing steps.


  MATERIALS AND METHODS
Top
Abstract
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
Acknowledgements
References

Construction of Mutant Subunits
Mouse cDNAs were used subcloned into the cytomegalovirus-based expression vector pRBG4 (Sine 1993 ). Mutant {alpha} subunits were constructed by bridging the naturally occurring sites BstX-1 and BspM-1 with synthetic double-stranded oligonucleotides (Bio-Synthesis Inc.), essentially as described previously (Bouzat et al. 1998 ). Single-stranded oligonucleotides were purified by PAGE and annealed before ligation. Restriction mapping and dideoxy sequencing on polyacrylamide gels confirmed all constructs.

Expression of AChR
Human embryonic kidney cells (HEK293) were transfected with {alpha} (wild-type or mutant), ß, {delta}, and {epsilon} cDNAs using calcium phosphate precipitation at a subunit ratio of 2:1:1:1 for {alpha}:ß:{delta}:{epsilon}, respectively, essentially as described previously (Bouzat et al. 1994 , Bouzat et al. 1998 ). For transfections, cells at 40–50% confluence were incubated for 8–12 h at 37°C with the calcium phosphate precipitate containing the cDNAs in DMEM plus 10% fetal bovine serum. Cells were used for single-channel measurements 1 or 2 d after transfection.

Patch-Clamp Recordings
Recordings were obtained in the cell-attached configuration (Hamill et al. 1981 ) at a membrane potential of -70 mV and at 20°C. The bath and pipette solutions contained 142 mM KCl, 5.4 mM NaCl, 1.8 mM CaCl2, 1.7 mM MgCl2, and 10 mM HEPES, pH 7.4. Patch pipettes were pulled from 7052 capillary tubes (Garner Glass), to produce final resistances of 5–7 M{Omega}, and were coated with Sylgard (Dow Corning Corp.). Pipette concentrations of ACh ranged from 10 to 1000 µM. Single channel currents were recorded using an Axopatch 200 B patch-clamp amplifier (Axon Instruments, Inc.) and digitized at 94 kHz with a PCM adapter (VR-10B; Instrutech). Data were transferred to a computer using the program Acquire (Bruxton Corp.) and detected by the half-amplitude threshold criterion using the program TAC 3.0 (Bruxton Corp.) at a final bandwidth of 10 kHz. Open- and closed-time histograms were plotted using a logarithmic abscissa and a square root ordinate (Sigworth and Sine 1987 ) and fitted to the sum of exponentials by maximum likelihood using the program TACFit (Bruxton Corp.). For each patch corresponding to a given ACh concentration, the number of opening events ranged from 4,000 to 10,000.

Open probability within clusters (Popen) was experimentally determined at each ACh concentration by calculating the mean fraction of time the channel is open within a cluster. The experimental Popen determinations were compared with theoretical dose–response curves calculated from either Scheme 1 or Scheme 2, using the fitted rate constants in Table 2.


 
View this table:
[in this window]
[in a new window]
 
Table 1. Open Probability, Mean Open Duration, and Mean Closed Duration of Clusters of AChRs Containing {alpha} wild-type, {alpha}T422A, and {alpha}T422S Subunits


 
View this table:
[in this window]
[in a new window]
 
Table 2. Kinetic Parameters for Mouse Wild-Type and {alpha}T422 Mutant AChRs

Kinetic Analysis
At ACh concentrations >10 µM, clusters of openings corresponding to a single channel were identified as a series of closely spaced events preceded and followed by closed intervals greater than a specified duration; this duration was taken as the point of intersection of the predominant closed-time component and the succeeding component in the closed-time histogram. To minimize errors in assigning cluster boundaries at the lower ACh concentrations, we analyzed only recordings from patches with low channel activity. Only clusters containing >10 openings were considered for further analysis. In addition, any clusters showing double openings were rejected. For each recording, kinetically homogeneous clusters were selected based on their mean channel open duration and open probability distributions (Wang et al. 1997 ). For each cluster within a recording, we calculated the Popen and mean open duration and plotted their distributions. Histograms were constructed and fitted by a Gaussian function, which yielded a mean and a standard deviation for each distribution. Typically, the distributions contained a dominant, approximately Gaussian component and minor contributions of clusters with different properties. Clusters showing mean open duration and open probability values within 2 SD of the mean of the major component were selected and retained for the kinetic analysis. Typically, >80% of the clusters were selected. For each recording, mean values of open probability, and mean open and closed durations changed little after the selection procedure (see examples of selection results in Table 1). The similarity of the evaluated parameters before and after selection indicates no qualitative changes in the dwell-time distributions due to the selection procedure. The resulting open and closed intervals, from single patches at several ACh concentrations, were analyzed according to a kinetic scheme using the program MIL (QuB suite, www.qub.buffalo.edu). In brief, MIL simultaneously fits recordings obtained at different ACh concentrations, and estimates rate constants in a kinetic scheme using a maximum likelihood method that corrects for missed events (Qin et al. 1996 ). We typically included recordings obtained at 10, 30, 50, 100, 200, and 300 µM ACh in the analysis of each mutant AChR. The dead time was typically 30 µs, but the fitted rate constants were stable over a range of dead times ranging from 18 to 40 µs (data not shown). Probability density functions of open and closed durations were calculated from the fitted rate constants and instrumentation dead time and superimposed on the experimental dwell-time histogram, as described by Qin et al. 1996 . Calculated rate constants were accepted only if the resulting probability density functions reasonably superimposed upon the experimental open- and closed-duration histograms at all ACh concentrations. To check the final set of fitted rate constants, open and closed intervals were simulated according to Scheme 1 or Scheme 2, using the program SIMU (QuB suite). Simulated channel events were analyzed and the resulting dwell times were fitted with the program MIL. The original rate constants were recovered with 10–15% accuracy. Comparison of the fits to different schemes was evaluated by the log likelihood ratio (LLR). For example, a LLR of 50 indicates that the probability of the model with the higher LL is e50 times more than that of the scheme with the lower LL (Salamone et al. 1999 ).


  RESULTS
Top
Abstract
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
Acknowledgements
References

We previously examined the M4 transmembrane domain by alanine-scanning mutagenesis, and found that {alpha}T422 was unique in that its mutation markedly affected the rate of channel closing (Bouzat et al. 1998 ). To further examine the structural contribution of {alpha}T422 to AChR activation, we replaced it by polar (serine and asparagine), aromatic (tryptophan and tyrosine), and nonpolar amino acids (alanine, valine, and cysteine). Then we transfected HEK293 cells with either wild-type or mutant {alpha} plus wild-type ß, {epsilon}, and {delta} subunit cDNAs and recorded single channel currents.

AChR channels were activated by a range of desensitizing concentrations of ACh (10–1,000 µM) to produce clear clusters of events corresponding to a single channel (Sakmann et al. 1980 ). Wild-type as well as mutant AChRs open in clusters of well-defined activation episodes at ACh concentrations >10 µM. Each activation episode begins with the transition of a single receptor from the desensitized to the activatable state and terminates by returning to the desensitized state (Fig 1). For both wild-type and mutant AChRs, closed intervals within clusters become progressively briefer with increasing ACh concentration (Fig 2). However, compared with wild-type, most of the mutant AChRs show both briefer openings and prolonged intracluster closings (Fig 1).



View larger version (31K):
[in this window]
[in a new window]
 
Figure 1. Clusters of single-channel currents from wild-type and mutant AChRs. Channels activated by 100 µM ACh were recorded from HEK cells expressing mouse wild-type ({alpha}2ß{epsilon}{delta}) and mutant {alpha}T422A, {alpha}T422V, {alpha}T422W, {alpha}T422Y, {alpha}T422C, {alpha}T422N, and {alpha}T422S AChRs. Currents are displayed at a bandwidth of 9 kHz with channel openings as upward deflections. Membrane potential, -70 mV.



View larger version (36K):
[in this window]
[in a new window]
 
Figure 2. Kinetics of activation of wild-type and mutant AChRs. (Left) Channel traces corresponding to AChRs containing the indicated {alpha} subunit and activated by 30 and 300 µM ACh. Currents are displayed at a bandwidth of 9 kHz with channel openings as upward deflections. (Center and right) Closed- and open-time histograms with the fit for Scheme 1 (wild-type and {alpha}T422S) or Scheme 2 ({alpha}T422A) superimposed. Histograms were constructed with the selected clusters. The ordinate corresponds to the square root of the fraction of events per bin. Closed-time histograms are truncated to eliminate slow components due to desensitization.

For both wild-type and mutant AChRs, the predominant peak of the closed-time histograms moves to shorter duration with increasing ACh concentration (Fig 2). However, for mutant AChRs containing alanine, closed intervals are clearly longer than those of wild-type AChRs at all ACh concentrations (Fig 2). Similar results were obtained for receptors containing nonpolar (valine and cysteine) or aromatic (tryptophan and tyrosine) substitutions. The predominant open-time component of either wild-type or mutant AChRs does not change with changes in ACh concentration (Fig 2). However for {alpha}T422A (Fig 2), as well as for {alpha}T422V, {alpha}T422W, {alpha}T422Y, and {alpha}T422C AChRs (not shown), open-time distributions shift to briefer durations compared with wild-type AChRs. In contrast, when T422 is replaced by serine (Fig 2) or asparagine (not shown), open- and closed-duration histograms are similar to those of wild type.

The observed changes in the kinetics of AChR activation could be due to changes in rate constants underlying either ACh binding or channel gating steps. To identify the kinetic steps affected by each M4 mutation, we fitted kinetic schemes to the open and closed dwell times. For wild-type AChRs, we used the classical activation scheme (Scheme 1), where two agonists (A) bind to the receptor (R) in the resting state with association rates k+1 and k+2 and dissociate with rates k-1 and k-2. Receptors occupied by one agonist open with rate ß1 and close with rate {alpha}1, and AChRs occupied by two agonist molecules open with rate ß2 and close with rate {alpha}2. At high agonist concentrations (>100 µM ACh) channel blockade is evident, requiring addition of the blocked state A2B.

To estimate the set of rate constants, Scheme 1 was fitted to the data using the program MIL (Qin et al. 1996 ). MIL computes the likelihood of the experimental series of open and closed times given a set of trial constants, and then changes the rate constants to maximize the likelihood. We simultaneously analyzed recordings obtained at multiple ACh concentrations (10–300 µM), with the aim of representing all the states in Scheme 1 in the fitting. The procedure for kinetic analysis requires selection of clusters to produce kinetically homogeneous data (MIL, QUB suite; Salamone et al. 1999 ). For each recording, kinetic homogeneity was determined by computing the mean channel open duration and open probability of each cluster, and then plotting their distributions. Only clusters showing Popen and mean open duration values within the mean ±2 SD were included in the subsequent analysis. Comparison of distributions before and after selection (see examples in Table 1) indicates that the selected clusters are representative of the entire population of events for each recording. In addition, open- and closed-time histograms constructed with selected data are very similar to those constructed with the original events (not shown). Data of open and closed intervals corresponding to the selected clusters were subsequently used to fit Scheme 1.

For wild-type AChR, ß2 was constrained to its known value (Sine et al. 1995 ; Wang et al. 1997 ; Salamone et al. 1999 ) because brief closings due to gating and channel blocking became indistinguishable at high ACh concentrations (Wang et al. 1997 ; Salamone et al. 1999 ). When ß2 was allowed to vary freely, MIL failed to converge to a well-defined set of rate constants and approached a value of ~100,000 s-1 (data not shown). Similar observations were recently reported by Salamone et al. 1999 . Although channel blocking is too brief to be completely resolved with our instrumentation dead time of 30 µs, we retain channel block in Scheme 1 because significant numbers of blockages are resolved at our highest ACh concentrations of 300 µM. Moreover, the unblocking rate constant obtained under the present conditions is similar to values previously reported (Maconochie and Steinbach 1995 ; Wang et al. 1997 ). Finally, association and dissociation rate constants were assumed to be equal at both binding sites, as previously described for adult mouse AChRs (Akk and Auerbach 1996 ; Wang et al. 1997 ; Salamone et al. 1999 ).

Rate constant estimates obtained for wild-type AChR, shown in Table 2, agree with those previously reported for mouse AChR (Wang et al. 1997 ; Salamone et al. 1999 ). The fitting analysis confirms that singly occupied AChRs open slowly and with low efficiency, but doubly occupied receptors open rapidly and efficiently.

We next applied kinetic analysis to AChRs containing mutations of {alpha}T422 using clusters selected as described for wild-type AChRs (Table 1). Scheme 1 could not adequately describe data obtained from AChRs containing the hydrophobic/aromatic mutations {alpha}T422A, {alpha}T422V, {alpha}T422Y, {alpha}T422W, and {alpha}T422C. Scheme 1 predicts a minor open-time component, associated with AR*, that decreases in relative area with increasing ACh concentration, but these mutant channels show two open components whose relative areas do not change with changes in ACh concentration. To account for two concentration-independent open-time components, we fitted Scheme 2 to dwell times from mutant AChRs. The classical activation scheme (Scheme 1) is a subset of Scheme 2, which has the opening step separated into two sequential steps. An alternative branched scheme, containing two open states connected to the A2R closed state, also described the data, but was found to be less likely for all mutant AChRs, as judged by maximum likelihood analysis. For example, for the {alpha}T422A AChR, Scheme 2 was e110 times more likely than the branched scheme.

In fitting Scheme 2 to data from AChRs containing hydrophobic or aromatic mutations, the two binding sites were assumed to have equal microscopic rate constants, as found for wild-type AChR. In contrast to fitting Scheme 1 to data from the wild type, it was possible to allow ß2 to vary freely. For the hydrophobic and aromatic mutations, the unblocking rate constant (k-b) was constrained to the value used for wild-type AChRs because M4 does not contribute to the channel pore (Blanton and Cohen 1994 ), and because we found the channel conductance to remain constant for all M4 mutations. The resulting fitting analysis establishes that the T422V, T422A, T422Y, T422W, and T422C mutations markedly affect the kinetics of channel gating (Table 2). The steps primarily affected are channel opening and closing, with changes of two- to sevenfold in the opening rate (ß2 in Scheme 2) and four- to ninefold in the closing rate ({alpha}2) (Table 2). Therefore, the doubly occupied mutant AChR opens with greater latency and reduced efficiency compared with wild-type AChRs. However, after opening, the mutant channels rapidly return to the closed state or, with low probability, reach a more stable open state, represented by A2R** in Scheme 2.

In contrast to the hydrophobic and aromatic mutations, AChRs containing the serine or asparagine mutations could be well fitted by Scheme 1 (Fig 2). Using the same constraints as in wild-type AChRs (ß2 = 49,000 s-1 and equivalent binding sites), the best-fit rate constants indicate that the primary effect of the serine and asparagine mutations is a twofold increase in the channel closing rate, {alpha}2 (Table 2). Because, as described above for wild-type AChRs, it was necessary to constrain ß2 in fitting data from {alpha}T422S and {alpha}T422N AChRs, and because the hydrophobic and aromatic mutations affected ß2, we systematically varied ß2 from 19,000 to 49,000 s-1. Based on likelihood, the best description of the {alpha}T422S and {alpha}T422N data was obtained when ß2 was 49,000 s-1, as found for wild-type AChR. For example, for the {alpha}T422S AChR, the fit using ß2 = 49,000 was e24, e81, and e235 times more likely than with ß2 equal to 39,000, 29,000, and 19,000 s-1, respectively.

Although dwell times from the hydrophobic and aromatic mutants were not well described by Scheme 1, dwell times from wild-type, {alpha}T422S, and {alpha}T422N AChRs could be well described by either Scheme 1 or Scheme 2 (Table 2). To show that changes in rate constants due to the mutations are not scheme dependent, we also fitted Scheme 2 to data from wild-type as well as {alpha}T422S and {alpha}T422N AChRs (Table 2). Comparison of log likelihoods, obtained after fitting both schemes to wild-type and {alpha}T422N data, reveals that Scheme 1 provides the better description; Scheme 1 is e75 and e170 times more likely than Scheme 2 for wild-type and T422N AChRs, respectively. On the other hand, for the {alpha}T422S mutation, Scheme 2 is e72 more likely than Scheme 1. Fitting Scheme 2 to data from the serine and asparagine mutations reveals only slight changes in the gating steps, similar to results obtained with Scheme 1 (Table 2). Comparison of wild-type and M4 mutant rate constants, calculated on the basis of Scheme 2, also shows that {alpha}T422 contributes primarily to channel opening and closing, and that these steps are largely affected by the presence of hydrophobic or aromatic residues (Table 2). Thus, residues preserving the hydrogen bonding ability of the original threonine support nearly normal gating kinetics.

The channel gating equilibrium constant, {theta}2, calculated as ß2/{alpha}2, significantly decreases in most of the mutants (Table 2). For each side chain at position 422, we calculated the net free energy change for the gating equilibrium, as well as the activation free energy changes for opening and closing transitions (Table 3). Free energy changes of ~2 kcal/mol in the gating equilibrium are observed when nonpolar or aromatic amino acids are substituted into position 422, whereas only slight changes occur when serine or asparagine are substituted. When threonine is replaced by either aromatic or nonpolar amino acids, the activation free energy for opening increases, whereas, for all mutants, activation free energy for closing decreases (Table 3).


 
View this table:
[in this window]
[in a new window]
 
Table 3. Differences in Free Energy of Channel Gating, and Open and Closing Reactions

Rates of agonist association and dissociation (Table 2) and dissociation equilibrium constants (Table 4) are not significantly affected by these mutations. Thus, T422 selectively affects channel opening and closing transitions.


 
View this table:
[in this window]
[in a new window]
 
Table 4. Equilibrium Constants for Mouse Wild-Type and {alpha}T422 Mutant AChRs

To determine the overall consequences of each mutation for AChR activation, and to evaluate consistency of the fitted rate constants, we computed the open probability from the fitted rate constants and compared it with the experimentally determined open probability (Fig 3). For wild-type AChRs, open probability increases with increasing ACh concentration, reaching a maximum of ~0.9 at 100 µM and showing an EC50 of 40 µM (Fig 3). For most mutant AChRs, both the EC50 and maximum open probability change. The open probability curves are shifted to higher ACh concentrations, owing to decreased efficiency of channel gating. For the mutants T422A, T422V, T422W, T422Y, and T422C, open probability is <0.6 at 1,000 µM, and the EC50 values range from 200 to 300 µM (Fig 3). Thus, mutations of {alpha}T422 affect both efficacy and EC50 for AChR activation. As expected, the profile for the T422N mutant is similar to that of wild-type AChR, and only slight changes are observed when threonine is replaced by serine. The EC50 and maximum open probability values for the T422S mutant are 80 µM and 0.8, respectively (Fig 3). The dose–response curves computed from the fitted rate constants closely follow the experimentally determined open probability determinations, demonstrating overall consistency of the data.



View larger version (25K):
[in this window]
[in a new window]
 
Figure 3. Agonist concentration dependence of the channel open probability. For AChRs containing the indicated {alpha} subunits, the mean fraction of time the channel is open during a cluster (Popen) was experimentally determined at the indicated concentrations of ACh. Each point corresponds to the mean ± SD of three to four patches. The curves were calculated from either Scheme 1 ({alpha}WT, {alpha}T422S, and {alpha}T422N) or Scheme 2 ({alpha}T422V, {alpha}T422A, {alpha}T422W, {alpha}T422Y, and {alpha}T422C), using the fitted rate constants in Table 2, and superimposed on the Popen measurements. Popen curves corresponding to wild type, T422S, and T422N AChRs calculated from Scheme 2 are almost identical to those shown in the figure. •, wild-type; bolded open circle, {alpha}T422N; {circ}, {alpha}T422S; {diamondsuit}, {alpha}T422C; {blacksquare}, {alpha}T422V; {square}, {alpha}T422A; {blacktriangleup}, {alpha}T422W; and {triangleup}, {alpha}T422Y.


  DISCUSSION
Top
Abstract
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
Acknowledgements
References

T422 of the M4 domain of the {alpha} subunit is presumed to be located at the lipid–protein interface, owing to its susceptibility to labeling by a hydrophobic probe (Blanton and Cohen 1994 ). Although M4 is the least conserved of the four transmembrane domains, the high degree of conservation of T422 among subunits and species indicates that it might play a structural or functional role. We previously reported that among residues in M4, {alpha}T422 is the most sensitive to mutation, with the alanine substitution decreasing the mean open time fivefold (Bouzat et al. 1998 ). Here we elucidated both mechanistic and structural bases for the contribution of {alpha}T422 to channel gating. Kinetic analysis of mutant {alpha}T422 AChRs disclose relatively small contributions to steps governing ACh binding to the resting, closed state of the AChR, but reveals significant contributions to steps underlying channel opening and closing.

Kinetic analysis requires that open and closed dwell times correspond to a single AChR channel. Clusters of opening events, corresponding to activation of the AChR between desensitized periods, were clearly distinguished at ACh concentrations >10 and 30 µM for wild-type and mutant AChRs, respectively. Because only clusters containing a single open conductance level were considered for kinetic analysis, it is possible to unequivocally assign the dwell times to sojourns in open and closed states of a single AChR. In addition, the clustering behavior observed for all mutant AChRs clearly demonstrates that, like wild-type, they can undergo cycles of desensitization and resensitization in the continued presence of ACh.

Based on gatings kinetics, one can distinguish the polar serine and asparagine from the nonpolar or aromatic amino acid mutations at position 422 as follows. (a) Kinetic changes due to replacement of T422 by serine or asparagine: recordings obtained from these mutant AChRs are similar to those from wild-type AChRs (Fig 1). Slight changes in the open- and closed-time histograms are observed (Fig 2). Data can be well fitted by the classical description of activation for wild-type AChR. The only detectable change is a twofold increase in the channel closing rate. (b) Kinetic changes due to replacement of {alpha}T422 by nonpolar or aromatic residues: recordings show very brief openings and prolonged closings. The classical activation scheme does not fit the experimental data. However, expanding the classical scheme to include two sequentially connected open states (Scheme 2) satisfactorily describes the data. Kinetic analysis based on Scheme 2 reveals that the rate of opening to the first open state is slower than the analogous step in wild-type AChRs, but the closing rate from this state ({alpha}2) is much faster. Because {alpha}2 is faster than the competing rate ß3, the open channel (A2R*) has a much greater probability of closing than of reaching the second, more stable open state (A2R**). Thus, nonpolar or aromatic mutations dramatically impair the initial opening step in the gating pathway. Interestingly, {alpha}T422C shows a different behavior within this group of mutants; it is the only mutant in which the opening rate is much more affected than the closing rate. Comparison of rate constants of wild-type and M4 mutant AChRs, calculated on the basis of both Scheme 1 and Scheme 2, confirms that {alpha}T422 contributes to channel opening and closing regardless of the kinetic scheme used for analysis. Scheme 2 was recently used to explain the activation of an AChR mutated in the M3 domain associated with a congenital myasthenic syndrome (Wang et al. 1999 ). This mutation, {alpha}V285I, also causes abnormally slow opening and rapid closing rates. Thus, Scheme 2 may represent a general mechanism that is not distinguished in wild-type AChR because the intermediate open state may be too short-lived to be detected, and the channel rapidly reaches the final open state. Thus, certain mutations that disrupt gating may have the ability to unmask intermediate open states.

The results also reveal the structural basis of the contribution of T422 to channel gating. Threonine and serine both contain polar hydroxyl groups, but differ by the presence of a methyl group attached to the beta carbon. Although the lack of a methyl group in serine slightly affects gating, {alpha}T422S AChRs are kinetically similar to wild-type AChRs. However, if only the hydroxyl group of threonine is replaced by a methyl group, as in valine, gating is substantially disrupted. Moreover, {alpha}T422V AChRs behave similarly to those containing structurally unrelated amino acids, such as tryptophan. The asparagine side chain provides similarly electronegative atoms to the original threonine side chain, and the polar amide group is similar to the hydroxyl group in its propensity for donating hydrogen bonds; T422N AChRs are almost kinetically identical to wild-type AChRs. Although cysteine can form hydrogen bonds, the kinetics of T422C AChRs differ from those of wild type. This finding can be explained by the fact that the thiol group is a poor hydrogen bond donor and forms weaker hydrogen bonds compared with hydroxyl or amide groups. Because hydroxyl and amide groups can equally well accept hydrogen bonds, and because, in contrast, cysteine is a poor hydrogen bond acceptor (Pal and Chakrabarti 1998 ), an alternative interpretation would be that the side chain at position 422 acts as a hydrogen bond acceptor. In addition, the lower pKa of cysteine compared with threonine, or the greater hydrophobicity of cysteine compared with threonine, serine, or asparagine may contribute to the impaired channel gating. The observed gating kinetics of T422Y AChRs indicate that tyrosine does not provide the correct hydrogen bonding capability, perhaps owing to delocalization of electron density from the hydroxyl oxygen to the aromatic ring, or to the bulk of the aromatic ring. Interestingly, the tyrosine hydroxyl group is a poorer hydrogen-bond acceptor than an aliphatic or water hydroxyl group (Fersht et al. 1985 ). In summary, our results show that correct gating requires a hydrogen bond involving the side chain at position 422 in the M4 domain of the {alpha} subunit.

T422 is located in the last third of the M4 domain, closest to the outer leaflet of the bilayer (Fig 4). Whether the M4 domain is an {alpha} helix (Blanton and Cohen 1994 ) or a ß-sheet (Unwin 1993 ) remains controversial. In addition, current models of the AChR place M4 with a tilt of ~30° with respect to the membrane (Ortells et al. 1997 ). Thus, atomic structure resolution of the AChR will be needed to unequivocally identify the hydrogen-bonding partner for {alpha}T422. Although T422 is labeled by the hydrophobic probe TID, and therefore is probably exposed to the lipid, its hydroxyl group could interact with a hydrogen bond acceptor or donor within a transmembrane domain of the same or another subunit. The more energetic hydrogen bonding expected in a hydrophobic environment may account for the substantial functional importance of the hydrogen bond with which T422 is involved. Thus, the present work demonstrates that subtle interactions originating far away from the two main functional domains of the AChR, the ion pore and the binding site, significantly affect channel gating.



View larger version (18K):
[in this window]
[in a new window]
 
Figure 4. Sequence alignments of the M4 domains. The conserved T422 is shown as a shaded column.

Calculation of free energy changes for channel opening and closing disclosed only slight decreases in free energy for closing between serine and asparagine mutant AChRs, but revealed significant changes in free energy for both opening and closing steps for all other mutant AChRs. The free energy differences are on the order of those calculated for hydrogen bonds between uncharged residues in an aqueous solution (0.5–1.5 kcal/mol; Fersht et al. 1985 ), or for other weak interactions such as dipole–dipole interactions. Our observation that a hydrogen bond is required for proper gating agrees with the calculated free energy changes. However, the correlation between the free energy of a hydrogen bond and free energy of activation for a state transition is not necessarily strict, given that disruption of a chemical bond can have global consequences.

The significant role of a hydrogen bond at position 422 in channel gating demonstrated in the current work explains the high degree of conservation of T422 among subunits and species (Fig 4). Abnormal activation of AChR has been shown to underlie congenital myasthenic syndromes (CMS) (Engel et al. 1998 ). Although mutations in the M4 domain have not been identified in CMS, our findings indicate that if they naturally occurred they could lead to a CMS.

Given the present understanding of the topology of the M4 domain, our results demonstrate that a lipid-exposed residue involved in hydrogen bonding is necessary for proper gating of the AChR.


Scheme 1.


Scheme 2.


  Footnotes

1 Abbreviations used in this paper: AChR, nicotinic acetylcholine receptor; M4, fourth transmembrane domain.


  Acknowledgements
Top
Abstract
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
Acknowledgements
References

We especially thank Dr. H.-L. Wang for advice and helpful discussion, Ana M. Roccamo, Dora Ortiz, and Horacio De Genaro for their expert technical assistance, and Nina Bren for providing two mutant subunits.

This work was supported by grants from Universidad Nacional del Sur, CONICET, and Agencia Nacional de Promoción Científica (C. Bouzat and F.J. Barrantes), National Institutes of Health grant NS-31744 (S. Sine), and Fogarty International Center grant 1R03 TW01185-01 (C. Bouzat and S. Sine).

Submitted: 6 December 1999
Revised: 27 March 2000
Accepted: 28 March 2000


  References
Top
Abstract
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
Acknowledgements
References

    Akabas, M., Karlin, A. 1995. Identification of acetylcholine receptor channel-lining residues in the M1 segment of the {alpha}-subunit. Biochemistry. 34:12496-12500[Medline].

    Akk, G., Auerbach, A. 1996. Inorganic, monovalent cations compete with agonists for the transmitter binding site of nicotinic acetylcholine receptors. Biophys. J. 70:2652-2658[Abstract].

    Blanton, M.P., Cohen, J.B. 1992. Mapping the lipid-exposed regions in the Torpedo californica nicotinic acetylcholine receptor. Biochemistry. 31:3738-3750[Medline].

    Blanton, M.P., Cohen, J.B. 1994. Identifying the lipid-protein interface of the Torpedo nicotinic acetylcholine receptor: secondary structure implications. Biochemistry. 33:2859-2872[Medline].

    Bouzat, C., Bren, N., Sine, S.M. 1994. Structural basis of the different gating kinetics of fetal and adult nicotinic acetylcholine receptors. Neuron. 13:1395-1402[Medline].

    Bouzat, C., Roccamo, A.M., Garbus, I., Barrantes, F.J. 1998. Mutations at lipid-exposed residues of the acetylcholine receptor affect its gating kinetics. Mol. Pharmacol 54:146-153[Abstract].

    Engel, A.G., Ohno, K., Milone, M., Sine, S.M. 1998. Congenital myasthenic syndromes. New insights from molecular genetics and patch-clamp studies. Ann. NY Acad. Sci. USA. 841:140-156[Free Full Text].

    Fersht, A.R., Shi, J.-P., Knill-Jone, J., Lowe, D.M., Wilkinson, A.J., Blow, D.M., Brick, P., Carter, P., Waye, M.M.Y., Winter, G. 1985. Hydrogen bonding and biological specificity analyzed by protein engineering. Nature. 314:235-238[Medline].

    Hamill, O.P., Marty, A., Neher, E., Sakmann, B., Sigworth, F.J. 1981. Improved patch-clamp techniques for high-resolution current recording from cells and cell-free membrane patches. Pflügers Arch. 391:85-100[Medline].

    Lee, Y.H., Li, L., Lasalde, J., Rojas, L., McNamee, M., Ortiz-Miranda, S.I., Pappone, P. 1994. Mutations in the M4 domain of Torpedo californica acetylcholine receptor dramatically alter ion channel fuinction. Biophys. J. 66:646-653[Abstract].

    Maconochie, D.J., Steinbach, J.H. 1995. Block by acetylcholine of mouse muscle nicotinic receptors, stably expressed in fibroblasts. J. Gen. Physiol. 106:113-147[Abstract].

    Ortells, M.O., Barrantes, G.E., Wood, C., Lunt, G.G., Barrantes, F.J. 1997. Molecular modelling of the nicotinic acetylcholine receptor transmembrane region in the open state. Prot. Eng. 10:511-517[Abstract].

    Ortiz-Miranda, S.I., Lasalde, J.A., Pappone, P.A., McNamee, M.G. 1997. Mutations in the M4 domain of the Torpedo californica nicotinic acetylcholine receptor alter channel opening and closing. J. Membr. Biol. 158:17-30[Medline].

    Pal, D., Chakrabarti, P. 1998. Different types of interactions involving cysteine sulfhydril group in proteins. J. Biomol. Struct. Dyn. 15:1059-1071[Medline].

    Qin, F.A., Auerbach, A., Sachs, F. 1996. Estimating single-channel kinetic parameters from idealized patch clamp data containing missed events. Biophys. J 70:264-280[Abstract].

    Sakmann, B.J., Patlak, J., Neher, E. 1980. Single acetylcholine-activated channels show burst-kinetics in the presence of desensitizing concentrations of agonist. Nature. 286:71-73[Medline].

    Salamone, F.N., Zhou, M., Auerbach, A. 1999. A re-examination of adult mouse nicotinic acetylcholine receptor channel activation kinetics. J. Physiol 516:315-330[Abstract/Free Full Text].

    Sigworth, F., Sine, S.M. 1987. Data transformation for improved display and fitting of single-channel dwell time histograms. Biophys. J. 52:1047-1052[Abstract].

    Sine, S.M. 1993. Molecular dissection of subunit interfaces in the acetylcholine receptor: identification of residues that determine curare selectivity. Proc. Natl. Acad. Sci. USA. 90:9436-9440[Abstract].

    Sine, S.M., Ohno, K., Bouzat, C., Auerbach, A., Milone, M., Pruitt, J.N., Engel, A.G. 1995. Mutation of the acetylcholine receptor alpha subunit causes a slow-channel myathenic syndrome by enhancing agonist binding affinity. Neuron. 15:229-239[Medline].

    Unwin, N. 1993. The nicotinic acetylcholine receptor at 9 Å resolution. J. Membr. Biol. 229:1101-1124.

    Unwin, N. 1995. Acetylcholine receptor channel imaged in the open state. Nature. 373:37-43[Medline].

    Wang, H.-L., Auerbach, A., Bren, N., Ohno, K., Engel, A.G., Sine, S.M. 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, H.-L., Milone, M., Ohno, K., Shen, X.-M., Tsujino, A., Batocchi, A.P., Tonali, P., Brengman, J., Engel, A.G., Sine, S.M. 1999. Acetylcholine receptor M3 domain: stereochemical and volume contributions to channel gating. Nat. Neurosci. 2:226-233[Medline].