The α2 glycine receptor (GlyR) subunit, abundant in embryonic neurons, is replaced by α1 in the adult nervous system. The single-channel activity of homomeric α2 channels differs from that of α1-containing GlyRs, as even at the lowest glycine concentration (20 µM), openings occurred in long (>300-ms) groups with high open probability (Popen; 0.96; cell-attached recordings, HEK-expressed channels). Shut-time intervals within groups of openings were dominated by short shuttings of 5–10 µs. The lack of concentration dependence in the groups of openings suggests that they represent single activations, separated by very long shut times at low concentrations. Several putative mechanisms were fitted by maximizing the likelihood of the entire sequence of open and shut times, with exact missed-events allowance (program hjcfit). Records obtained at several glycine concentrations were fitted simultaneously. The adequacy of the different schemes was judged by the accuracy with which they predicted not only single-channel data but also the time course and concentration dependence of macroscopic responses elicited by rapid glycine applications to outside-out patches. The data were adequately described only with schemes incorporating a reaction intermediate in the activation, and the best was a flip mechanism with two binding sites and one open state. Fits with this mechanism showed that for α2 channels, the opening rate constant is very fast, ∼130,000 s−1, much as for α1β GlyRs (the receptor in mature synapses), but the estimated true mean open time is 20 times longer (around 3 ms). The efficacy for the flipping step and the binding affinity were lower for α2 than for α1β channels, but the overall efficacies were similar. As we previously showed for α1 homomeric receptors, in α2 glycine channels, maximum Popen is achieved when fewer than all five of the putative binding sites in the pentamer are occupied by glycine.
Out of the four known isoforms of the glycine receptor (GlyR) α subunit, α2 (Grenningloh et al., 1990; Akagi et al., 1991) may play a major role in the embryonic form of GlyRs (Legendre, 2001), as its transcripts, abundant in the embryonic and neonatal rat, are gradually replaced by α1 in the adult (Akagi et al., 1991; Malosio et al., 1991; Watanabe and Akagi, 1995). This is supported by a variety of functional observations. Takahashi et al. (1992) observed longer apparent open times in recombinant α2 versus α1 GlyRs, and in embryonic versus adult rat spinal GlyRs. The decay of glycinergic IPSCs also speeds up with development. The simplest explanation is that embryonic and neonatal GlyRs contain mostly α2 subunits, whereas adult receptors are α1β heteromers (Becker et al., 1988; Langosch et al., 1988; Hoch et al., 1989). Nevertheless, it is far from certain that the slow glycinergic synaptic currents observed early in development are mediated by α2 homomeric receptors. Homomeric GlyRs probably cannot cluster at the synapse, because only β subunits can interact with gephyrin (Meyer et al., 1995). Second, α2 is expressed in areas where no glycine-mediated synaptic transmission has been detected (Flint et al., 1998; Mangin et al., 2002) and by nonneuronal cells, such as glial or progenitor cells (Belachew et al., 1998; Nguyen et al., 2002). Mangin et al. (2003) showed that α2 homomeric receptors are unsuitable for fast synaptic transmission because they activate and deactivate slowly and have a low probability of opening in response to the rapid application of brief glycine pulses (Mangin et al., 2003). It has therefore been suggested that α2 homomeric receptors may be mainly extrasynaptic and that they can be activated by a paracrine-like release of glycine or taurine during neuronal development (Legendre, 2001; Clements, 2002; Mangin et al., 2003). In addition to that, synaptic currents recorded in neonatal hypoglossal motoneurons decay almost 10 times faster than currents elicited on recombinant α2 channels (Mangin et al., 2003), leading Singer et al. (1998) to suggest that they are mediated by α2β heteromeric GlyR.
Global fitting of kinetic models to single-channel data from recombinant α1-containing GlyRs led us to describe an activation mechanism (“flip”) that explicitly incorporates intermediate shut states between agonist binding and channel opening (Burzomato et al., 2004). Mechanisms that include pre-open shut states such as flip states, or for muscle nicotinic receptors, “primed” states (Mukhtasimova et al., 2009), are attractive because they offer a suggestive kinetic parallel to the events that are likely to happen in the extracellular domain before channel opening, such as domain closing and/or C-loop capping. In addition to that, they have been found to describe well channel activation for a variety of wild-type and mutant subunit combinations for glycine, muscle nicotinic, and GABAA receptors, notably with respect to the differences between full and partial agonists in the nicotinic superfamily and to the complex structure of unliganded openings in gain-of-function mutants (Burzomato et al., 2004; Plested et al., 2007; Lape et al., 2008; Mukhtasimova et al., 2009; Keramidas and Harrison, 2010). Our aim in the present work was to extend this approach to recombinant α2 GlyRs, to establish a quantitative activation mechanism for this channel, and to test whether a flip-type mechanism could provide an adequate description for both single-channel and macroscopic data from recombinant α2 homomers. For the first time in fitting GlyR channel activity, we found that a flip-type mechanism describes the data somewhat better than a less constrained mechanism incorporating distal shut states (such as the mechanism proposed by Jones and Westbrook, 1995).
MATERIALS AND METHODS
GlyR α2 subunit cDNA construct, HEK293 cell culture, and transfection
HEK293 cells (American Type Culture Collection and LGC Standards) were maintained in a humidified incubator at 37°C (95% air/5% CO2) in Dulbecco’s modified Eagle’s medium supplemented with: sodium pyruvate (0.11 g/liter), heat-inactivated fetal bovine serum (10% vol/vol), and penicillin G (100 U/ml)/streptomycin sulfate (100 µg/ml; all from Invitrogen). Cells were passaged every 2–3 d, up to 30 times. Cells were plated on glass coverslips (13-mm diameter, sterilized over flame), placed in 35-mm sterile Petri dishes containing 1.5–2 ml Dulbecco’s modified Eagle’s medium, 4–6 h before transfection, to allow cells to attach. Cells used for fast concentration jumps were plated onto coverslips coated with poly-l-lysine (Sigma-Aldrich) by incubation of the coverslips for ∼40 min, followed by washing with autoclaved water and oven drying at 60°C.
HEK293 cells were transiently transfected with the rat GlyR α2 cDNA (Isoform A; available from GenBank/EMBL/DDBJ under accession no. AJ310837; full-length sequence verified) in the pcDNA3.1 vector. Transfection was by calcium phosphate–DNA coprecipitation (Groot-Kormelink et al., 2002). In brief, cDNA was added to CaCl2 solution (340 mM in sterile water) at a volume ratio of 1:5, and this mixture was added drop-by-drop to an equal volume of 2× Hank’s buffered saline (280 mM NaCl, 2.8 mM Na2HPO4, and 50 mM HEPES, pH 7.2 with NaOH) to form the precipitate. The final mixture was then added dropwise over the plated cells. The cDNA used for transfection contained 3–20% α2/pcDNA3.1, 18% pEGFP-c1, and 62–79% pcDNA3.1 vector lacking the α2 coding sequence. Adding empty vector was found to be important for the health of the transfected cells and to ensure that the level of α2 GlyR expression was optimal for recording. The total amount of the final cDNA mixture per plate was kept constant at 3 µg. Recordings were performed 2–3 h after washing and no more than 3 d after transfection. We attempted to obtain single-channel recordings also from α2β heteromeric GlyR, but the expression level was consistently too low, despite considerable effort to improve it.
Cell-attached single-channel currents were recorded at 19–21°C with thick-walled borosilicate pipettes (with filament; Harvard Apparatus) coated near the tip with Sylgard (Dow Corning) and fire-polished before use to a final resistance of 5–12 MΩ. Electrodes were filled with extracellular medium (in mM: 20 Na-gluconate, 102.7 NaCl, 2 KCl, 2 CaCl2, 1.2 MgCl2, 10 HEPES, 20 TEA-Cl, 15 sucrose, and 14 glucose, pH 7.4 with NaOH; osmolarity 320 mOsm) plus the appropriate concentration of glycine. All solutions were prepared in HPLC-grade water (VWR International) to minimize contamination by glycine and filtered through a 0.2-µm Cyclopore track-etched membrane (GE Healthcare) to remove impurities that can block electrodes or affect the quality of the seal. The bath level was kept low to reduce noise.
Currents were recorded with an amplifier (Axopatch 200B; MDS Analytical Technologies) at a pipette potential of +100 mV, with no correction for junction potential, as this is calculated to be at most +1 mV (Clampex 9; MDS Analytical Technologies). Single-channel currents were prefiltered at 10 kHz by the amplifier’s four-pole Bessel filter and saved on a digital audio tape with a recorder (DRA-200; Bio-Logic). Recordings were replayed, filtered at 5 or 7 kHz, and sampled at 50 or 71.43 kHz, respectively, depending on the signal-to-noise ratio, with a Digidata 1322A and Clampex software (MDS Analytical Technologies).
Single-channel current analysis
Single-channel currents were idealized by time-course fitting (SCAN program). Segments showing drifts in the baseline, seal breakdowns, or dubious channel openings were excluded, and adjacent intervals were marked unusable as appropriate. After imposing a resolution of 20–25 µs (for both shut and open times), the number of resolved time intervals (sum of shut-time and open-period intervals) ranged from 498 to 2,368 in the 12 patches used in the maximum likelihood fits.
Idealized current amplitudes longer than twice the filter rise time (66 µs at 5 kHz) were plotted as an amplitude histogram (EKDIST program). The main conductance level accounted for >99% of all fitted open intervals, so we pooled together all openings for kinetic analysis. The average amplitude of openings to the main conductance level was 5.9 ± 0.2 pA (n = 25 patches) but varied considerably, ranging from 4.5 to 9.1 pA. In the 12 patches used in our analysis, the range of mean amplitudes was 4.9–8.0 pA (6.1 ± 0.3 pA). This variability in amplitude was expected for cell-attached records because it will be affected by the intracellular chloride concentration and the resting membrane potential, both of which vary from cell to cell. Variability in the intracellular chloride concentration and in the true transmembrane potential is a source of error in so far as it affects the rate constant values in receptor activation.
The distributions of apparent shut times and open periods were fitted with a mixture of exponential probability density functions. Each open period was defined as the duration of time for which the channel remained continuously open, regardless of amplitude; that is, the time between two adjacent shut times that were both longer than the imposed resolution.
The true number of the channels in each patch is unknown, but for kinetic analysis, it is essential to isolate stretches of openings that are likely to arise from one individual channel. This was done in the initial stage of the analysis by fitting shut-time distributions with EKDIST to determine a critical shut time (tcrit) that would allow us to divide the record into true one-channel stretches of openings. Fits to the distributions obtained at EKDIST are only descriptive and not mechanism dependent; hence, parameters estimated by EKDIST were not used further in the analysis. The criterion for a group of openings to be judged to have come from one individual channel is essentially that the probability of being open within the groups is high enough that we can be almost certain that a double opening would be seen if more than one channel were present (Colquhoun and Hawkes, 1990). Groups with fewer than 10 open periods were excluded. For fitting, we tested different values of tcrit (with or without the use of Colquhoun-Hawkes-Srodzinski vectors; see below). We found that the most consistent fits were obtained when we used the tcrit values reported in the Results.
For kinetic analysis, the postulated mechanisms were evaluated using maximum likelihood fitting and the program HJCFIT. For this purpose, the idealized data from all concentrations were grouped into three independent sets, with each set containing four patches at different agonist concentrations (20 or 30 µM, 100 µM, and 1 and 10 mM glycine). All data from one set were input into the program to be fitted simultaneously to a specified mechanism, together with the resolution of 20–25 µs (the duration of the fastest event that can be detected), a tcrit value and the glycine concentration for each patch. The program calculates a likelihood value for the dataset from the initial guesses for the rate constant values, taking into account the imposed resolution, which is used to correct for missed events (Hawkes et al., 1990, 1992). The rate constant values are then adjusted to maximize the likelihood until a maximum is reached.
Shut times longer than tcrit cannot be used directly in the fitting process because they are affected by the (unknowable) number of channels in the patch. Nevertheless, some information is recoverable from the long shut times. The true (one-channel) shut time between one activation and the next is not known, but it is known that it must be at least as long as tcrit. This information can be used to improve the fit by replacing the normal steady-state vectors as the start and end of each group of openings, with the Colquhoun-Hawkes-Srodzinski vectors described by Colquhoun et al. (1996; Eq. 5.11). This was done for low concentration records (in this case, all concentrations except 10 mM), and simulations show that it is useful (Colquhoun et al., 2003). It cannot be used at concentrations that are high enough that desensitization occurs, because it is model dependent and our models do not include desensitized states.
We tested whether our fits were robust to the choice of initial guesses for the rate constant values to exclude the possibility of convergence to local maxima in the likelihood surface. If the fit is good and the likelihood surface has only one well-defined maximum, changing the initial guesses should not affect the estimates for the rate constants. This was not always the case, and sometimes several equivalent solutions (with similar likelihood values) were possible, suggesting that the likelihood surface was rough with several maxima. In this case, the best solution was chosen on the basis of fit quality, physical plausibility of the rate constant estimates, and predictive power of the fitted rate constants. For example, we discarded fits that gave physically impossible association rate constant values (i.e., >109 M−1s−1). As a test of fit quality, we used several data displays: shut-time distributions, open-period distributions, and the predicted open probability (Popen)–concentration curve. We also used the results of the fits to simulate realizations of single-channel activity (SCSIM program) and to calculate macroscopic current relaxations to a glycine concentration jump (SCALCS program), and we compared them to our experimental observations. All analysis programs are available from http://www.ucl.ac.uk/Pharmacology/dcpr95.html.
Macroscopic currents: recording and analysis
Macroscopic currents evoked in outside-out patches by agonist concentration pulses were recorded with thick-walled borosilicate pipettes (with filament; Harvard Apparatus), fire-polished to a final resistance of 8–15 MΩ when filled with “low [Cl−] solution” (in mM): 121.1 K-gluconate, 1 CaCl2, 1 MgCl2, 10 HEPES, 11 EGTA, 6 TEA-Cl, and 2 MgATP, pH 7.2, adjusted with KOH; osmolarity 325 mOsm, adjusted with sucrose. Because of the very low Cl− concentration in this solution, the silver/silver chloride wire of the pipette holder was re-chlorided before each experiment to compensate for the increased flow of Cl− ions from the wire to the solution. The bath solution was the same as for single-channel recordings. All outside-out recordings were performed at a pipette holding potential of −100 mV. We did not correct for junction potential, and this gave a maximum error of +14.6 mV in low [Cl−] solution. During the 3–5-min intervals between glycine applications, patches were held at −60 mV to prolong the stability of the seal.
Glycine, dissolved in bath solution, was applied to outside-out patches with a theta tube (Hilgenberg GmbH) cut to a final diameter of ∼150 µm at the tip. The tube was driven by a piezo stepper (Burleigh Instruments, Inc.). The exchange time was measured by the application of 30% diluted bath solution before the experiment (to optimize the electrode position) and after the rupture of the patch. Macroscopic currents were recorded with an Axopatch 200B, digitized with Digidata 1322A, and saved directly on a computer via Clampex software (5 kHz filtering, 20–50 kHz sampling rate; all from MDS Analytical Technologies).
Different glycine concentrations were applied in random order for a duration chosen to allow the current response to reach a plateau or a peak. All currents were corrected for rundown using a standard glycine concentration (about EC30), which was applied every third response. Concentration–response curves were plotted using these corrected responses (averaging up to five responses), and each experiment was fitted to the Hill equation by least squares with equal weights, using the CVFIT program to get estimates of maximum current, Hill coefficient (nH), and EC50. Numbers in the text are means ± SD of the mean (SDM) from the individual fits of each experiment. For the purpose of display, responses were normalized to the fitted maximum response of the same experiment and pooled to give an average concentration–response curve.
To study the kinetics of macroscopic currents, 10–50 responses were recorded for each glycine concentration at intervals of at least 20 s to allow full recovery from desensitization (Gentet and Clements, 2002). Responses at each concentration for each patch were averaged, excluding failures or responses that contained patch breakdowns. Only experiments in which the rundown between the first and last three responses was <5% were included in the analysis. The rise time for average currents and tip potentials was measured as the time from 20 to 80% of the peak response using Clampfit 9.2 software; experiments in which the open tip response had a 20–80% exchange time slower than 300 µs were rejected. Both the time course of desensitization at prolonged glycine applications (200 ms) and the time course of deactivation at short applications (2 ms) were fitted with a sum of one or more exponentials.
For the purpose of clarity, we shall first define four terms that are used throughout to describe patterns of channel openings: group, activation, burst, and cluster.
The term “group” of openings is used as a generic description of any group of openings that is separated by short shut times, regardless of why they occur.
An “activation” is what happens between the start of the first opening of a channel until the end of the last opening before the receptor returns to the unoccupied resting state. The activation is the elementary event produced by agonist binding, consists of a burst of one or more openings that occur in quick succession, and is what underlies the deactivation rate after a macroscopic jump to zero concentration (see the appendix of Wyllie et al., 1998).
The term “burst” refers to the groups of openings that occur at very low agonist concentrations. At sufficiently low agonist concentrations, the observed bursts of openings are a close approximation to channel activations (Edmonds et al., 1995; see Discussion). The shut times between activations contain one or more sojourns in the unoccupied resting state, as well as bindings that happen not to produce an opening. They may also contain sojourns in desensitized states.
A “cluster” of openings, defined as in Sakmann et al. (1980), refers to the long groups of openings that are observed at high agonist concentrations. A cluster consists of many activations that occur in quick succession (because of the high agonist concentration), so one activation cannot be distinguished from the next one. Clusters are separated by shut times that have one or more sojourns in long-lived desensitized states. They are seen with most ligand-gated ion channels but are not easy to identify in α2 receptor records.
Properties of macroscopic currents from α2 GlyRs
To understand the α2 GlyR channel behavior in nonequilibrium conditions, we recorded macroscopic currents activated by fast application of glycine. We have recently shown (Pitt et al., 2008) that the gating of α1β GlyRs is modulated by internal Cl−. Channel deactivation after short glycine applications becomes faster in outside-out patches when intracellular chloride is lowered to concentrations similar to those likely to be present in the intact HEK cell (i.e., in cell-attached recordings). As this is likely to occur also for α2 homomeric receptors, we recorded the macroscopic jumps using a low (10-mM) intracellular chloride concentration to make sure that in the different experiments, channels were exposed to similar chloride concentrations so that the information from different patch configurations could be combined in our investigation of the activation mechanism. Traces in Fig. 1 A are the outside-out current responses of α2 homomeric GlyRs to the fast application of 2-ms pulses of 10 mM glycine. The results of fitting the time course of current responses to short pulses of 0.5–10 mM glycine are shown in Table I. Two exponential components were needed to fit the current decay, and the slower component was the most prominent by far.
We also recorded outside-out α2 GlyR macroscopic responses to long (200-ms) concentration jumps at saturating glycine concentrations (10 mM; −100 mV; low [Cl−] pipette solution). A typical current example is shown in Fig. 1 B. All patches (n = 6) revealed a complex desensitization phase that could be best fitted with mixtures of three exponential components: 1.9 ± 0.3 ms (0.6 ± 0.2% fractional area), 23 ± 3 ms (2.2 ± 0.6%), and 820 ± 180 ms (97.2 ± 0.7%). This complexity of desensitization was independent of current amplitude: there was no apparent relationship between the peak current amplitude of the different patches and the number of the components needed for fitting currents.
Finally, to characterize the concentration–response relation for α2 channels, we performed concentration jumps at a range of glycine concentrations. The traces in Fig. 1 C show examples of averaged currents recorded in the outside-out configuration using a low [Cl−] pipette solution (10 mM chloride; −100 mV). Concentration–response curves obtained from six patches gave an average EC50 of 770 ± 100 µM and an nH of 1.99 ± 0.25 (Fig. 1 D).
General features of α2 GlyR single-channel activity
Fig. 2 A shows representative glycine-activated single-channel currents recorded in the cell-attached configuration (pipette holding potential of +100 mV) from HEK293 cells transfected with rat α2 GlyR cDNA. It was immediately clear that the single-channel activity of α2 GlyRs is very different from that of α1 (Beato et al., 2004) or α1β GlyRs (Burzomato et al., 2004). Openings of α1 or α1β GlyRs occur in short bursts at low glycine concentrations (Fig. 2 B, bottom trace). These bursts are activations that, as glycine concentration is increased, become more closely spaced because the time spent in unliganded resting states gets shorter until, in high glycine, they become grouped into long, high Popen clusters, separated by long shut periods.
In apparent contrast to what is seen with α1 or α1β receptors, α2 activity appeared to be very similar at all glycine concentrations (Fig. 2, A and B). Even at 20 µM, the lowest glycine concentration at which single-channel currents could be detected, groups of α2 channel openings were long and separated by long shuttings (usually lasting for minutes). Because of the length of the shut times, each recording contained a relatively small number of these groups of openings and few long shut times. Another consequence of the small number of long shut times that could be recorded is that it was impossible to measure any changes in their length as glycine concentration was increased. Moreover, the “experimental” success rate was generally low, and only a minority of α2 patches displayed any channel activity (∼1 in 5 patches at glycine concentrations of 100 µM or greater and only ∼1 in 10 at 20–50 µM). We could find no single-channel activity in records at 10 µM glycine, despite recording from >70 good patches for 5–25 min.
After idealization of single-channel records by time-course fitting, the distributions of apparent open and shut times were plotted and fitted with mixtures of exponential probability density functions. The distributions of apparent open periods are shown in Fig. 2 A (middle column). These were fitted with up to three exponential components (Table II). The longest component was predominant at all concentrations, with an average area of 82–98%, and had a time constant of ∼30 ms at all concentrations other than 20 µM glycine, where it appeared to be about threefold shorter. Both of the two faster components represent a very small number of intervals and were poorly defined, and a single-exponential density fitted the distribution quite well in most cases. The absolute number of apparent openings or shuttings per patch was small, despite the length of the clusters and the length of the recordings. Apparent shut-time distributions (Fig. 2 A, right column) were fitted with a minimum of three components (Table III). At all glycine concentrations, most shuttings were contributed by the shortest component (area, 85–99%). Given its average time constant of 5–10 µs, the vast majority of the events in this component were shorter than the experimental resolution of 20–25 µs and thus were missed. The number of shut-time intervals that fell into the two longest components was very small, and consequently, these components were not very well defined. The time constants of these two components appeared to decrease about threefold when glycine concentration increased from 20 µM to 10 mM. The existence of more than one shut-time component at saturating concentrations (Fig. 2 A and Table III) is not compatible with mechanisms that have only one shut state when fully liganded. Additional fully bound shut states are therefore required, either as intermediates during channel opening or otherwise distal to the binding steps.
The records shown in Fig. 2 A show a superficial resemblance to the clusters of openings that are seen with other receptors at concentrations of agonist that are high enough to produce long desensitized intervals, but it seems that the appearance is deceptive. These long groups of openings appear to be unusually long individual activations of the channel, rather than clusters of many activations. This interpretation is confirmed by the analysis below.
The best way to display and characterize the concentration dependence of single-channel activity is to plot a one-channel concentration–Popen curve. This is done by identifying stretches of openings that are likely to be the expression of the activity of a single-channel molecule by setting for each idealized record a critical shut-time duration tcrit value (refer to Materials and methods). We attempted to construct a single-channel concentration–Popen curve for α2 receptors in this way, but it proved to be impossible. A typical display of the Popen values calculated for groups of openings, pooled and averaged, is shown in Fig. 2 C: the Popen was high, well above 90% at all agonist concentrations. Similar results were obtained at all the values of tcrit tested (1–100 ms). This is because (see Discussion) these groups represent channel activations, rather than clusters of several activations. Hence, the Popen data found in this way yielded no information about the concentration dependence of channel activity.
Maximum likelihood fit of kinetic mechanisms: initial considerations
Next, we performed maximum likelihood fitting of kinetic mechanisms to our single-channel data, using HJCFIT. Single-channel records were grouped into sets, and each set contained four patches that spanned the range of glycine concentrations (refer to Materials and methods). Because (as above) we were not able to assess the concentration dependence of the single-channel Popen (e.g., Burzomato et al., 2004; Lape et al., 2008), for α2 receptors we had to rely on the macroscopic concentration–response curve (Fig. 1 D) for the choice of concentrations. Ideally, to provide sufficient information on the different steps of the mechanism, each set should include records from low, intermediate, and saturating concentrations of glycine.
Initial attempts to fit a range of concentrations up to 1 mM gave inconsistent results. However, the macroscopic concentration–response curve (Fig. 1 D) suggests that 1 mM glycine is not a saturating concentration. Much better results were obtained when single-channel recordings at 10 mM glycine were included in the fits. This inclusion of 10 mM of data was essential to obtain good agreement of the macroscopic deactivation time constant predicted by the fit with our measurements in the concentration jump experiments and to predict the distinctive, unusually long single-channel activations, separated by long shut times.
For fitting, the idealized records are divided into groups of openings that are separated by a shut time that is longer than a specified tcrit. The aim of this is to obtain groups of openings that come from the same individual channel (refer to Materials and methods; Colquhoun et al., 1996, 2003). Because there was some ambiguity in the choice of tcrit (see about Popen curve above), we did several preliminary fits with different tcrit values. For most mechanisms tested, the fits were consistently more stable when the tcrit was set to 100 ms at 20–30 µM, 10 ms for bursts at higher concentrations (up to 1 mM), and 7–10 ms at 10 mM glycine.
To test the quality of a model, the predictions from a fit were compared with the experimental apparent shut-time and apparent open-period distributions. In this case, open and shut times were almost independent, so conditional distributions of open times were not a useful criterion for distinguishing mechanisms (Colquhoun and Hawkes, 1987; Colquhoun et al., 2003; Burzomato et al., 2004). Additional elements for evaluating the quality of a mechanism were the physical plausibility of the estimated rate constant values and their consistency across different sets. Finally, simulations of steady-state single-channel records and calculations of macroscopic currents (concentration jumps), on the basis of the model and rate constants that were fitted, were compared with our experimental data. In total, 44 different putative mechanisms were fitted (see Discussion). Most failed to describe the data well, and only a few of the most plausible candidates will be described here.
Scheme 1: a mechanism with two binding sites and two extra shut states
Kinetic Scheme 1 (Fig. 3 A) has been proposed by Mangin et al. (2003) as a minimal kinetic mechanism for the α2 homomeric GlyR on the basis of the analysis of macroscopic and single-channel currents activated by fast agonist application to outside-out patches. Scheme 1 contains two short-lived “desensitized” states, each connected to a bound shut state. This scheme resembles closely the one proposed for the activation of GABAA receptors by Jones and Westbrook (1995). Good fits of apparent open periods and shut times could be obtained, but only if the binding sites were allowed to interact (i.e., no constraints imposed on the binding rates), as for the analogous model for α1β receptors (Burzomato et al., 2004).
Examples of distributions of apparent open periods and shut times for a set fitted with Scheme 1 are plotted in Fig. 3 B. As shown from this example, the fit of Scheme 1 predicts well the shape of the apparent shut- and open-time distributions. The results from the fits to three datasets are summarized in Table IV. As seen from the table, the values for the gating rate constants (the shutting rate constant α and the opening rate constant β) were highly consistent across the three sets, with a very low coefficient of variation.
Rate constants other than those for opening and shutting were highly variable across the three sets (coefficient of variation [CVM] of 50% or more), suggesting that the scheme cannot estimate these rate constants reliably. Note also that the EC50 predicted from the fitted rate constants was almost sixfold lower than the EC50 measured in outside-out recordings. To test further the correctness of Scheme 1, we calculated the macroscopic current expected in response to 2-ms fast concentration jumps with 10 mM glycine. The current response was calculated for each set of fitted rate constants. The predicted rise time of macroscopic currents from the three sets fitted had a mean time constant of 10 ± 5 ms, which is ∼25-fold slower than the 20–80% rise time measured from currents evoked by 10 mM glycine (τ = 0.42 ± 0.08 ms; Table I). On the other hand, the mean deactivation time constant predicted by Scheme 1 was 120 ± 70 ms (n = 3 sets), which is very close to the experimental value of 130 ± 10 ms (the slow component; Table I). This corresponds to a mean burst length calculated from Scheme 1 of 120 ± 70 ms at 1 nM glycine (n = 3 sets).
We also checked whether simulated single-channel records (SCSIM program) resembled the pattern seen in experimental records. Single-channel simulations at 20 µM glycine using Scheme 1 failed to reproduce the most striking features we observed in our single-channel experimental records, namely the presence of groups of openings with duration of up to a few seconds long, separated by very long gaps lasting for minutes. For two sets, bursts (with lengths ranging from tens of milliseconds to ∼1 s) appeared continuously in the simulated records, and there was no sign of the very long shut times between them, as seen in experimental records (compare Fig. 3 C with Fig. 2 A). Only the simulation from the third set showed long gaps (Fig. 3 D), but even this realization did not show long, high Popen groups of openings. The mean burst length calculated from Scheme 1 at 20 µM glycine was only 300 ± 30 ms (n = 3 sets). These simulations show that Scheme 1 could not describe our experimental results adequately, despite a good match to the dwell-time distributions. Therefore, we proceeded to test alternative mechanisms for our data.
Scheme 2: a flip mechanism with two binding sites and one open state
Kinetic Scheme 2 (Fig. 4 A) is a variant of the flip mechanism, originally proposed by Burzomato et al. (2004) for α1β and α1 GlyRs. Scheme 2 is a modified version containing two binding sites (rather than the three needed to reach maximum Popen in α1 homomers). We have also chosen to omit from Scheme 2 the two partially liganded open states included in the α1/α1β flip mechanism. Including more than one open state in the mechanism barely improved the quality of the fits to the present data and resulted in poorly determined rate constant estimates for the entry and exit from partially liganded open states. This probably means that openings from partially liganded receptors are rare, so there are not enough transitions in our records to be able to describe these states accurately. The predominance of visits to a single open state is confirmed by the lack of any substantial correlation between open and shut times, and by the fact that the distribution of apparent open-period durations is not far from a single exponential.
An example of a dataset fit with Scheme 2 is shown in Fig. 4 B. There is a good agreement between the experimental distributions (histograms) and those that are predicted from the fits using the model (continuous line). The set shown in Fig. 4 B is the same as in Fig. 3 B for comparison. The average rate constants from three different sets are summarized in Table V. Note that these fits were done with the assumption that binding sites are independent; i.e., they do not interact, so all binding steps to a given conformation (R or F) are identical, regardless of how many sites are already occupied. This was done by imposing the constraint that k+1 = k+2 and k–1 = k–2.
As seen from the values in Table V, the fits with Scheme 2 are very consistent across the three sets, and most rate constants had low coefficients of variation. The least well-defined rate constants are those for entering and exiting the singly occupied flip state (δ1 and γ1), although they are clearly quite slow. This is not surprising given the evidence already cited that partially liganded states occur only rarely. Scheme 2 again predicts that the efficacy for the opening step is very high (350), and that the efficacy for the flipping steps increases as the number of bound glycine molecules increases (Table V). Also, the shut flipped states have a 40–60-fold higher affinity for glycine than the shut resting states. This is a well-characterized feature of the flip mechanism, which suggests that the binding of the agonist stabilizes the receptor in a higher affinity form reached by a conformational change that precedes channel opening (Burzomato et al., 2004). The average maximum Popen predicted by the fits to Scheme 2 was 0.99, and the predicted nH was 1.89, in good agreement with the value of 1.99 (nH; Fig. 1 D) obtained from our macroscopic dose–response curve experiments (see also Mangin et al., 2003) and with the maximum Popen of 0.99 from single-channel records. Note, however, that there was about a threefold difference between the experimental macroscopic EC50 (770 µM) and the one predicted by the model (240 µM; Table V).
Using the rate constants obtained from the fits described above, we calculated macroscopic current responses to a 2-ms fast concentration jump with 10 mM glycine. Scheme 2 predicted a time constant of 171 ± 6 ms (n = 3 sets) for the decay of these currents, which is similar to the 130 ± 10 ms we obtained from our macroscopic currents in outside-out patches. Scheme 2 predicted two significant components for the activation phase of the current macroscopic current (1.6 ± 0.2 and 0.19 ± 0.02 ms), but the contribution of the faster component was minimal. Consequently, the weighted calculated rise time (τ = 1.5 ± 0.2 ms) reflected mostly the slow component and was almost fourfold slower than the rise time measured from currents evoked by 10 mM glycine (20–80% rise time, 0.42 ± 0.08 ms; Table I).
Finally, steady-state single-channel records (at 20 µM glycine) simulated using the fitted rate constants showed groups of openings with a mean burst length of 490 ± 70 ms, separated by long shut times (>1 min; Fig. 4, C and D). These features of single channels are similar to our cell-attached recordings. It is not surprising that we saw no channel activity at 10 µM glycine, because calculations based on the fitted rate constants predict that activations should be separated by long silent periods of >4 min (see Discussion).
To summarize, based on the above observations and on preliminary fits to different models, Scheme 2 is the minimal model that can describe best both our single-channel and macroscopic data.
Scheme 3: a flip mechanism with three binding sites and one open state
The mechanism in Scheme 2 contains two binding sites. This is different from what we suggested earlier for the α1 homomeric GlyR, namely that maximal activation of the receptor requires the occupancy of three binding sites by glycine (Beato et al., 2004). It is of course possible that the number of functional binding sites is different between α2 and α1 homomers, but this is an important point to test. We therefore extended Scheme 2 and tested two related mechanisms (Scheme 3a and Scheme 3b), shown in Fig. 5 A. The two mechanisms differ only in that the monoliganded flipped state AF is omitted from Scheme 3b. As for Scheme 2, we assumed that binding sites are identical (within each form of the receptor) and do not interact. This was imposed by using the following constraints: k+1 = k+2 = k+3 and k–1 = k–2 = k–3 for both R and F states. Note that because of these constraints, the number of free parameters in Schemes 2 and 3b is the same, and both have nine free rate constants.
The apparent quality of the predictions of the dwell-time distributions was good and very similar for the two models; therefore, we show only the fits of Scheme 3a in Fig. 5 B, using the same set of data as in the previous figures for comparison. The predictions of the dwell-time distributions were almost identical for Schemes 2, 3a, and 3b (compare Fig. 4 B and 5 B), and the values of the rate constants obtained by fitting the three schemes (see Table V) were very similar, although the coefficient of variation was overall higher for the fits with Scheme 3b. The most variable rate constants are again those that describe the transitions between partially liganded flipped and resting states.
Scheme 2 is favored over Schemes 3a and 3b because of its better prediction of the Hill slope (1.89 ± 0.01 for Scheme 2, but 2.49 ± 0.06 for Scheme 3a and 2.49 ± 0.05 for Scheme 3b; see also Table V). Because we were not able to obtain a one-channel Popen curve, our only experimental estimate of a slope (nH = 1.99; Fig. 1 D) comes from the concentration–response curve obtained from the macroscopic current records in the outside-out configuration. We do not know how much this measurement is distorted by desensitization, so this is not a conclusive argument.
To check further the quality of Schemes 3a and 3b, we calculated macroscopic currents and simulated single-channel currents. Fig. 5 C shows a simulated single-channel trace (20 µM glycine), using the rate constants from a set fitted with Scheme 3a. In this record, openings occur in groups several hundred milliseconds long, separated by very long shut times, as in our experimental traces. Similar results were obtained from simulations with Scheme 3b (not depicted). Therefore, as with Scheme 2, Schemes 3a and 3b predicted well our single-channel data. The rise and deactivation time constants of macroscopic currents were calculated from the fitted rate constants. The average weighted rise time of a calculated current activated by 2 ms of 10 mM glycine was 1.5 ± 0.4 ms (n = 3 sets), and the deactivation time constant was 150 ± 30 ms. Both time constants were almost indistinguishable from those calculated for Scheme 2. The relationship between burst length and macroscopic deactivation is considered in the Discussion.
The main aim of this paper is to characterize the activation mechanism of homomeric α2 GlyRs. We used the method of maximum likelihood fitting of kinetic mechanisms to idealized single-channel data (Colquhoun et al., 2003) for the α2 GlyR. In brief: (a) single-channel data are idealized by time-course fitting, and patches of different agonist concentrations are grouped into sets; (b) a kinetic mechanism is postulated (usually taking into account exponential fits to experimental dwell-time distributions) and initial guesses are given for the rate constants; (c) the idealized records, the appropriate resolution, and the postulated mechanism are introduced into a program (HJCFIT) that can calculate the likelihood (i.e., the probability of the entire sequence of apparent open and shut times, given the current estimates of the rate constants) and maximize it by optimizing rate constants; and finally, (d) the validity of the model and of the final rate constants is checked by comparing their predictions with the experimental data (dwell-time histograms). For the α2 GlyR, we also examined how well the mechanism and rates fitted to single channels predicted the time course and concentration dependence of the macroscopic currents elicited by concentration jumps. The terms activation, burst, and cluster are defined at the beginning of the Results section.
Properties of α2 GlyR single-channel activity
The first kinetic analysis on GlyRs, performed on cultures of spinal motorneurons, showed a characteristic feature of these channels: the frequency of long and short openings is concentration dependent. In recordings at low glycine, increasing glycine concentration increases the proportion of long openings (without changing the time constants of the different opening components; Twyman and Macdonald, 1991). The same phenomenon was subsequently detected in recombinant homomeric and heteromeric GlyRs that contain the α1 subunit, where bursts at higher glycine concentrations are longer because they contain more openings and because the openings themselves appear to be longer (Fucile et al., 1999; Beato et al., 2002, 2004; Burzomato et al., 2004).
This is not the case for the activity of α2 GlyRs, where openings occur in long groups, regardless of concentration (Fig. 2 A). At first, the length of the groups of openings (Fig. 2 A) misled us into thinking that they were clusters (as defined above). Even at the lowest concentration of glycine at which channel activity could be obtained (20 µM), long groups of openings were observed (several hundreds of milliseconds), and the probability of being open within the group of openings was high (0.96; tcrit, 100 ms). However, the macroscopic concentration–response curve (Fig. 1 D) suggests that 20 µM glycine should produce only a very small response and therefore a very low overall single-channel Popen (<5%). This observation, together with the fact that we could not observe activity at concentrations lower than 20 µM, implies that the long groups of openings at 20 µM glycine must be (unusually long) individual activations of the channel (possibly lengthened a bit by rebinding). In other words, the true Popen must be low at 20 µM and should be measured including the very long shut periods that separate the groups of openings. This cannot be done because there is no way to tell whether two adjacent activations originate from the same individual channel.
The use of concentration jumps for fast agonist application provides the benefit of mimicking what occurs in vivo: a fast and brief exposure of the receptors to saturating concentrations of neurotransmitter in the cleft during synaptic transmission (Clements et al., 1992; Colquhoun et al., 1992). Moreover, when the application of agonist is short (1–2 ms), the decay of the current is known to reflect the duration of the burst at very low agonist concentration, i.e., the duration of channel activations (Wyllie et al., 1998). When jumps of 2 ms to three different glycine concentrations were applied, two deactivation components were identified. The time constant of the slow component, which is by far the predominant one, was ∼130 ms, which is ∼17 times slower than the deactivation of α1β channels (7.5 ms with 10 mM intracellular Cl−; Pitt et al., 2008). Note that Mangin et al. (2003) reported for α2 currents, using high chloride in the pipette, a deactivation time course similar to the one we observed (153 ms). We found (unpublished data) that changing intracellular chloride concentration to 130 mM had little effect on α2 currents, in contrast with our findings for α1β.
Choice of mechanism
Despite the fact that the single-channel records look very different from those seen with α1-containing homomers or heteromers, it was found eventually that the best description among the 44 mechanisms that we tested was provided by the flip model (Burzomato et al., 2004; Lape et al., 2008). An overview of the type of mechanisms tested is shown in Fig. 6.
The fit of the shut-time distributions was not perfect with this mechanism, and neither was the prediction of the rise time for concentration jumps. However, both were better than could be obtained with the main rival class of mechanisms, as represented here by the Mangin/Jones-Westbrook type of model (Scheme 1). In our original work on the flip model (Burzomato et al., 2004), the Jones-Westbrook–type mechanism fitted at least as well as flip, albeit with four more free parameters. The preference for the flip model then was based to a certain extent on physical grounds (it allowed a fit without having to postulate interaction between different binding sites), partly because the flip model had a plausible structural interpretation (C-loop capping, since confirmed in muscle nicotinic channels by Mukhtasimova et al., 2009). Mechanisms with preopening intermediate states work well for wild-type and mutant GlyR and muscle nicotinic channels and for GABAA channels (Burzomato et al., 2004; Plested et al., 2007; Lape et al., 2008; Mukhtasimova et al., 2009; Keramidas and Harrison, 2010). In the present case, there was more reason to discard a Mangin/Jones-Westbrook type of model in favor of flip than in our earlier work. For example, the fits with flip give more consistent rate constant values across experimental sets. Flip is also somewhat better at predicting EC50 values, the rise time of macroscopic currents, and the pattern of channel activity in single-channel records.
There are two respects in which the fit with the flip mechanism (Scheme 2), despite being the best we found, is less than perfect. The EC50 predicted by the fit to Scheme 2 was almost three times smaller than that observed for the macroscopic concentration–response curve. Also, the predicted rise time for responses to pulses of 10 mM glycine was over three times slower than the observed value (although that is a great deal closer than for Scheme 1, for which it was 25 times too slow). It could well be that the flip model is too simple to account for these details. The least plausible aspect of the flip mechanism is the postulate that flipping is a concerted conformation change. It might be better to allow flipping (or “priming”) to occur separately in each subunit (Mukhtasimova et al., 2009). Desirable though that might be, in its general form, the primed model has so many free parameters that they cannot all be estimated from the data. We have preferred an approximation that allows all of the rate constants to be estimated.
Features of the activation mechanism for α2 GlyRs: the number of open states and binding sites
The general flip model that describes the activation of homomeric α1 and heteromeric α1β GlyRs contains three open states (Burzomato et al., 2004). This is in line with early data on native and recombinant receptors that detected at least three distinguishable open states (Twyman and Macdonald, 1991; Beato et al., 2002, 2004; Lewis et al., 2003). In contrast to the α1-containing channels, α2 receptors could be fitted quite well with a flip model that contains only one open state. This suggests that partially liganded openings are either rare or too short to be detected. The latter explanation is supported by preliminary fits with a version of Scheme 2 that included monoliganded openings (not depicted). Partially liganded openings had a very high closing rate constant and were so short that, if they truly occurred, most of them would be missed with our resolution of 20 or 25 µs. For these reasons, we considered it more sensible to leave only one fully liganded open state. Another reason for this choice is that we could not detect open–shut correlations in our experiments with α2 receptors, and these are not expected to occur in mechanisms with only one open state (Colquhoun and Hawkes, 1987).
In the general flip scheme used to fit homomeric α1 channels, maximum Popen is reached when three (out of the possible five) binding sites are occupied by glycine. For the α2 channel, good, almost identical fits of single-channel data were obtained with flip mechanisms with two or three binding sites, and these resulted in very similar values for the rate constants. The biggest difference between the predictions of mechanisms with two or three binding sites was in the value of the Hill slope of the concentration–response curve. Scheme 2 (two binding sites) predicted a value of 1.9, in better agreement with our experimental data from macroscopic currents in the outside-out configuration (1.99; Fig. 1 D). The limitation is that we only have a macroscopic concentration–response curve, which is less satisfactory than a single-channel Popen curve because it is susceptible to errors from desensitization, as well as possible differences between outside-out patches and cell-attached patches. Even so, it is clear that occupancy of all the five potential agonist binding sites in a homomeric nicotinic-type receptor is not needed to achieve the maximum Popen, in agreement with our own conclusions for the α1 GlyR homomer (Beato et al., 2004) and with those of the elegant experiments of Rayes et al. (2009) on α7-5HT3 chimeric receptors.
Features of the activation mechanism of α2 GlyRs: efficacy and affinity
For the purpose of characterizing the kinetic behavior of homomeric α2 GlyRs, steady-state single-channel data were fitted to 44 different putative reaction schemes, only four of which are described here. These potential mechanisms fall into two distinct types. In Scheme 1, originally proposed by Mangin et al. (2003), the saturated receptor has only two states, open and shut, and the efficacy of the agonist depends only on the equilibrium constant for this step, E = β2/α2, in the classical way. The maximum Popen is
To get good fits, it was necessary to add two extra shut states, in the manner of Jones and Westbrook (1995). Although these states are labeled as “desensitized,” they are quite short-lived, and there is no reason to think that they have anything to do with macroscopic desensitization. The extra shut states are, therefore, entirely arbitrary. Furthermore, to obtain a good fit, it was necessary to assume that the two binding sites in Scheme 1 interact (the binding constant for the second binding is not the same as for the first).
For both of these reasons, we turned to flip mechanisms (Schemes 2, 3a, and 3b). In this class of mechanisms, the binding sites are supposed to be independent, and there is an extra shut state (the flipped conformation) between resting and open states. Therefore, the fully saturated receptor has three states (resting, flipped, and open) rather than two, and the maximum Popen depends on two different equilibrium constants. These are the equilibrium constant, E, for the open–shut reaction of the fully liganded receptor (the “opening efficacy”), and the equilibrium constant for the preceding conformation change from resting to flipped state, F (the “flipping efficacy”),
This can be written in the classical form if we define an “effective efficacy,” Eeff, as
With this definition, the maximum Popen for the flip type mechanism can be written in the form of Eq. 1 as:
The flip mechanism was first proposed for glycine α1β receptors by Burzomato et al. (2004). It follows from Eq. 3 that if the flipping efficacy, F, is small, the overall efficacy may be much lower than expected from the opening efficacy, E, and Lape et al. (2008) gave evidence that a reluctance to flip (small value of F) was the basis of partial agonism for both glycine and nicotinic receptors. A similar explanation for the partial agonism of choline at nicotinic receptors was proposed by Lape et al. (2009).
The fits of the flip mechanism (Scheme 2) suggested that glycine has a 46-fold higher affinity for the flipped conformation than for the resting conformation, little different from the 65-fold higher affinity for α1β GlyRs found by Burzomato et al. (2004). The results resemble α1β GlyRs also in that the increase in affinity is brought about almost entirely by an increase in the association rate constant for the flipped conformation, rather than a decrease in the dissociation rate. The association rate constant for the resting conformation, ∼105 M−1s−1, is unusually low. Even the association rate constant for the flipped conformation, ∼8 × 106 M−1s−1, is still slower than that for the binding of acetylcholine to muscle-type nicotinic receptors, which is >108 M−1s−1. Although the high affinity ratio (between resting and flipped states) results in a high degree of apparent cooperativity (the Hill slope is close to its maximum value of 2), the absolute value of the flipping equilibrium constant is only near 1, compared with 27 for α1β GlyRs (Lape et al., 2008). Nevertheless, the overall efficacy is high because of the high equilibrium constant (E2) for the shut-open step.
The values of the gating rate constants were the most consistent across experimental sets in our maximum likelihood fitting, regardless of which model was used for fitting. They ranged between 300 and 400 s−1 for the shutting rate α and 110,000 and 130,000 s−1 for the opening rate β. The values given in Table V for the fit of Scheme 2 have the following implications. The fast opening rate constant is very similar to that for the heteromeric α1β receptors (Burzomato et al., 2004), but the shutting rate constant is 20 times slower, so the mean length of the individual openings is ∼3 rather than 0.14 ms for heteromeric α1β receptors. The mean efficacy for the opening step is ∼350, much higher than the gating efficacy of either α1β heteromeric GlyR (α = 7,000 s−1; β = 129,000 s−1; E = 20) or α1 homomeric GlyR (α = 700 s−1; β = 28,000 s−1; E = 38; Burzomato et al., 2004). In α1-containing receptors, the lower (but still large) efficacy for the open–shut step stems from a faster shutting rate and implies a short mean open time for the α1 heteromer, but in the α1 homomer, it appears to result from a slower opening rate constant. For α2 receptors, the flipping equilibrium constant for saturated receptors (Table V) is around 1, so shut saturated receptors spend about half their time in the resting conformation (A2R) and half the time flipped (A2F). The overall effective efficacy is still very high, Eeff =171, from Eq. 3, which implies a maximum Popen of 99.4%, from Eq. 4.
The energy diagram in Fig. 7 plots the energy barriers for the transitions of fully liganded receptors from resting to flipped to open for the GlyRs we have so far analyzed (α1, α1β, and α2). It is interesting to note that the height of the energy barrier that the channel has to cross to go from resting to flipped seems to be linked in our series to whether the principal subunit in the receptor is α1 or α2. It is plausible that the subunit that provides the principal side of the binding site should play a major part in this transition to flipped, a transition that is likely to correspond to the conformational changes that follow agonist binding and spread through the extracellular domain to reach the channel gate (Grosman et al., 2000; Lee and Sine, 2005). Almost all of the sequence differences between α1 and α2 are in their extracellular domain, whereas there is only one amino acid difference in M2, the main domain that lines the channel and gates it shut. On the other hand, the M2 sequence is less conserved in the β subunit. Notably, the energy difference between the initial resting state and the final open state is very similar for the three receptors in the series, corresponding to the fact that the effective efficacy Eeff (Eq. 3) of glycine is very high for all three.
The mechanistic interpretation of α2 GlyR bursts
It was shown rigorously by Wyllie et al. (1998; appendix) that the time constants (and mean length) for bursts (activations) measured at sufficiently low agonist concentrations must be the same as the time constants for the decay of macroscopic concentration jumps to zero concentration. Obviously, bursts cannot be measured at zero concentration, so this theorem begs the question, What concentration is “sufficiently low”? For most sorts of receptor, it is possible to record bursts at a concentration that is low enough for there to be good agreement between the single-channel mean burst length and the macroscopic deactivation time constant. According to our interpretation, the α2 receptor is an exception to this rule.
The observed time constant for the predominant slow component of the macroscopic deactivation was 130 ms after a step from 10 mM glycine to zero (Table I). The fit of Scheme 2 predicts a predominant deactivation time constant of 171 ms, not far from that observed. On the other hand, this value does not agree very well with the observed mean burst length of 490 ms measured at the lowest concentration of glycine at which measurements could be made, 20 µM. Evidently, quite a lot of reopening from monoliganded channels must occur at 20 µM glycine, and this reopening prolongs the burst length to a value that is well above its low concentration limit, even at 20 µM glycine. Calculations from the fitted rate constants for Scheme 2 predicts that at a really low glycine concentration, 1 nM, the mean burst length should be 171 ms, in agreement with the predicted deactivation rate for Scheme 2. Clearly, the way the numbers work out in this case makes it impossible to record bursts at a sufficiently low concentration.
To make sense of the numbers, it is useful to define the relative probabilities that a channel moves next to each of the adjacent states in the scheme, regardless of how long it takes to do so. For example, in Scheme 2 (Fig. 8), a monoliganded resting channel AR (state 5) can move next to the resting state, R (state 6), it can flip to AF (state 3), or it can bind another glycine molecule to get to A2R (state 4). The probabilities are simply proportional to the rate constants for leaving state 5 (Colquhoun and Hawkes, 1995; p. 415), so, for example, the probability that a monoliganded channel dissociates is
where c is the concentration of glycine.
The concentration dependence of the length of an activation arises, in the case of α2 receptors, largely from the binding of glycine to the monoliganded flip state (AF). When 20 µM glycine is applied, opening will take place usually via A2R because the monoliganded flipping rate is small (π54 π42 is almost 20 times greater than π53 π32). After an opening, the channel must return to A2F (state 2), from where it is likely to reopen (π21 = 0.96), but if it does not reopen, it is three times less likely to go to AF (state 3) than to unflip to A2R.
If this occurs, even at the lowest concentration of glycine for which recordings could be made, 20 µM, the probability that from AF the channel can bind glycine again is high (π32 = 0.98). Because this binding returns the channel to A2F, and from there the channel is very likely to (re)open, this path accounts for the lengthening of the activation at 20 µM glycine. In contrast to that, the probability of a second binding occurring to the resting conformation, AR, is very low (π54 = 0.004 in 20 µM glycine), and the probability of flipping is even lower.
There are two reasons for the unusually long activations produced by the α2 receptor even at the lowest concentrations. The first is that individual openings are longer, on average, for the α2 receptor (3 ms; Table V) than for the α1β receptor (0.14 ms; Lape et al., 2008). The second reason is the large number of openings per activation seen with the α2 receptor, around 100, even at very low concentrations.
Shut times within an activation turn out to have a rather simple interpretation in the case of the α2 receptor. The distribution of shut times calculated from the fitted rate constants for Scheme 2 (Table V), with Eq. 3.90 from Colquhoun and Hawkes (1982), is dominated at all concentrations by a very short component with a mean of ∼8 µs, close to the mean duration of a sojourn in A2F (state 2). Clearly, most of the time during an activation is spent oscillating between A2F* and A2F. There is also a small component (a few percent in area) of shut times within activations that has a time constant that is close to the mean lifetime in A2R (state 4), ∼0.2 ms. This clearly represents channels that undergo transitions A2F*→A2F→A2R→A2F→A2F*. There will very rarely be oscillations between A2F and A2R, so this time constant of the shut-time distribution is close to the time spent in A2R (plus two brief sojourns in A2F). Still more rarely (at 20 µM glycine), there will be transitions A2F*→A2F→AF→A2F→A2F*. Again, there will be little oscillation, so these generate a small component of shut times with a mean life close to that in AF. At low concentrations, a channel that after opening reaches AR will very rarely rebind glycine and reopen, but this route will become increasingly common as glycine concentration is increased, and it will contribute to the increase in duration of activations with concentration.
Finally, the shut time between one activation and the next at a low concentration is predicted by the fit to be very long. At 20 µM glycine, the range is 25.6–92.9 s for the three sets used for fitting. The fact that it is so long explains why so few activations could be recorded, and it also explains why the Popen is so low (around 0.01) at 20 µM glycine. The shut time between activations is characterized by at least one sojourn in the resting state, R. However, even at 20 µM glycine, the mean lifetime of a single sojourn in R is quite short, a few hundred milliseconds. So, evidently, there is a lot of oscillation between shut states between one activation and the next. Consider oscillations between the resting state, R (state 6), and the monoliganded resting state, AR (state 5), and back. The distribution of the number of such oscillations can be found as follows. From R, the channel has to reach AR (state 5; i.e., π65 is 1). There will be no R–AR oscillations if the channel then moves straight to either state AF (state 3) or A2R (state 4) without returning to R. Therefore, the probability of there being no R–AR oscillations is P(0) = π53 + π54 = (1 − π56). The probability of one oscillation (R→AR→R→AR→state 3 or 4) must be P(1) = π56 (1 − π56). So, in general, the probability of r oscillations between R and AR before leaving for state 3 or state 4 is seen to be
These probabilities (from r = 1 to ∞) sum to 1, and the mean number of oscillations is
For example, in the fitted set that was used for the histograms, π56 = 0.996 at 20 µM glycine, so the mean number of R→AR→R oscillations is 229. The mean length of a sojourn in R was 251 ms, and the mean length of the monoliganded state AR is 2.17 ms. Thus, the mean length of a run of R→AR→R oscillations is 229 (251 + 2.17) = 58 s. This is almost as long as the mean shut time between activations (74 s in this case). Thus, the very long shut times that separate long activations are due to the many R→AR→R oscillations rather than to the long life of any individual shut state.
This work was funded by the Medical Research Council (Program Grant G0400869).
Christopher Miller served as editor.
P. Krashia and R. Lape contributed equally to this paper.
F. Lodesani’s present address is Advanced Center for Electronic Systems ARCES II School of Engineering, University of Bologna, 40125 Bologna, Italy.