Ancestral acetylcholine receptor β-subunit forms homopentamers that prime before opening spontaneously

  1. Christian JG Tessier
  2. Raymond M Sturgeon
  3. Johnathon R Emlaw
  4. Gregory D McCluskey
  5. F Javier Pérez-Areales
  6. Corrie JB daCosta  Is a corresponding author
  1. Department of Chemistry and Biomolecular Sciences, University of Ottawa, Canada
  2. Centre for Chemical and Synthetic Biology, University of Ottawa, Canada

Abstract

Human adult muscle-type acetylcholine receptors are heteropentameric ion channels formed from two α-subunits, and one each of the β-, δ-, and ε-subunits. To form functional channels, the subunits must assemble with one another in a precise stoichiometry and arrangement. Despite being different, the four subunits share a common ancestor that is presumed to have formed homopentamers. The extent to which the properties of the modern-day receptor result from its subunit complexity is unknown. Here, we discover that a reconstructed ancestral muscle-type β-subunit can form homopentameric ion channels. These homopentamers open spontaneously and display single-channel hallmarks of muscle-type acetylcholine receptor activity. Our findings attest to the homopentameric origin of the muscle-type acetylcholine receptor, and demonstrate that signature features of its function are both independent of agonist and do not necessitate the complex heteropentameric architecture of the modern-day protein.

Editor's evaluation

This compelling study outlines how phylogenetic reconstruction can yield important insights into evolutionary aspects of assembly and function in acetylcholine receptors (AChRs). The authors elegantly demonstrate that ancestral AChR β subunits can form homomeric channels that share important functional hallmarks with their modern cousins, including unliganded gating. The work provides an intriguing framework to evaluate the evolution of the broader family of pentameric ligand-gated ion channels.

https://doi.org/10.7554/eLife.76504.sa0

Introduction

Ligand-gated ion channels convert chemical signals into electrical impulses by coupling the binding of small molecules to the opening of an ion-conducting transmembrane pore (Hille, 2001). Of the many types of ligand-gated channels, the superfamily of pentameric ligand-gated ion channels (pLGICs) is the largest and most structurally and functionally diverse (Jaiteh et al., 2016). Formed from five identical or homologous subunits arranged around a central ion-conducting pore, pLGICs are found in almost all forms of life (Corringer et al., 2012). In prokaryotes, pLGICs appear to be exclusively homopentameric, while in eukaryotes both homo- and heteropentameric channels are common. Humans express more than 40 different pLGIC subunits, which are subdivided based on whether they form cation-selective channels activated by acetylcholine or serotonin (5-hydroxytryptamine [5-HT3]), or anion-selective channels activated by γ-aminobutyric acid or glycine. This repertoire of pLGIC subunits, combined with their ability to form homo- and heteropentamers, leads to a wealth of structural and functional diversity, presumably to meet a variety of synaptic needs (Sine and Engel, 2006).

The archetypal pLGIC is the heteropentameric muscle-type nicotinic acetylcholine receptor (AChR), with the human adult form composed of two α-subunits, and one each of the β-, δ-, and ε-subunits (Figure 1; Changeux, 2020a; Changeux, 2020b). To gain insight into the structure, function, and evolution of the AChR, we have employed an ancestral reconstruction approach (Emlaw et al., 2021; Emlaw et al., 2020; Prinston et al., 2017; Tessier et al., 2017). Previously we resurrected a putative ancestral AChR β-subunit. Referred to as ‘βAnc’, this subunit differed from its human counterpart by 132 amino acids (i.e. approximately 30% of the total amino acid sequence), yet was able to substitute for the human β-subunit, and also supplant the human δ-subunit, forming functional hybrid ancestral/human AChRs (Emlaw et al., 2021). These hybrid AChRs were ternary mixtures, containing two βAnc subunits, two human α-subunits, and one human ε-subunit. A concatameric construct confirmed that the two βAnc subunits resided next to each other, occupying positions in the heteropentameric complex usually reserved for the human β- and δ-subunits (Figure 1; Emlaw et al., 2021). Regardless of whether they bind agonist or not, all pLGIC subunits have both principal (+) and complementary (−) subunit interfaces. To sit next to each other, the (+) and (−) interfaces of any two neighbouring subunits must be structurally compatible. Thus, for βAnc to replace both the human β- and δ-subunits, the (+) and (−) interfaces of βAnc must be compatible with each other, raising the possibility that multiple βAnc subunits could coassemble to form βAnc homomers (Figure 1).

Subunit composition of heterologously expressed acetylcholine receptors (AChR).

Subunit stoichiometry and arrangement of the human adult muscle-type AChR (left), where the agonist-binding sites at the α–δ and α–ε subunit interfaces are indicated with asterisks (*). A reconstructed ancestral β-subunit (βAnc; purple) forms hybrid AChRs (middle) where βAnc substitutes for the human β-subunit (β; orange) and supplants the human δ-subunit (δ; green). The principal (+) and complementary (−) interfaces of βAnc must be compatible for two βAnc subunits to sit side-by-side (red highlight), which predicts that homomers formed from multiple βAnc subunits should be possible (right, boxed).

Here, using single-channel measurements, we demonstrate that βAnc readily forms homopentameric channels that open spontaneously in the absence of acetylcholine. This spontaneous activity displays hallmarks of the muscle-type AChR, including steady-state single-channel burst behaviour. These findings demonstrate how fundamental characteristics of AChR activation are independent of agonist, and not a result of the complex heteropentameric architecture of the muscle-type AChR. Finally, an alternate ancestral β-subunit, reconstructed using a phylogenetic tree that matched the accepted species relationships (Emlaw et al., 2021), and which only shared ~85% identity with βAnc, revealed that these unexpected characteristics of βAnc are robust to phylogenetic uncertainty, and thus deeply embedded within AChR β-subunit structure and evolutionary history.

Results

Our first inkling that βAnc may be able to form homomeric channels came from heterogeneity in single-channel recordings acquired after alterations to our original heterologous expression/transfection protocol. Typical cotransfection of cDNAs encoding human α-, δ-, and ε-subunits with a cDNA encoding βAnc leads to robust cell surface expression of βAnc-containing hybrid AChRs (Prinston et al., 2017). In an attempt to lower overall AChR expression with the purpose of facilitating single-channel analysis, we reduced the total amount of subunit cDNA in our transfections (~sixfold), while maintaining the same 2:1:1:1 subunit cDNA ratio (by weight; α:βAnc:δ:ε). Reducing the amount of cDNA lowered overall AChR expression as expected, but also led to heterogeneity in our patches (Figure 2A), which was not present in original single-channel recordings of βAnc-containing AChRs (Prinston et al., 2017). Instead of a single population of channels with a uniform amplitude of ~10 pA, and a burst behaviour indicative of βAnc-containing AChRs (Figure 2A; left inset), we also observed a second class of channels with a different kinetic signature, and an increased amplitude (Figure 2A; right inset). A similar trend was not observed when cDNA encoding the wild-type human β-subunit was cotransfected instead of βAnc (Figure 2—figure supplement 1), indicating that βAnc was the source of the heterogeneity. Consistent with this, lowering the proportion of βAnc cDNA in the transfection mixture reduced the fraction of high amplitude channels (Figure 2B and D), while transfecting exclusively with βAnc cDNA resulted in patches where all channel openings had a uniformly high amplitude (Figure 2C and D). This demonstrated that when transfected alone, βAnc forms functional ion channels. The present work stems from this unexpected observation and describes characterisation of these previously unobserved channels.

Figure 2 with 1 supplement see all
Single-channel recordings of βAnc-containing channels.

(A) Representative continuous recording from a cell-attached patch where cells were transfected with cDNAs encoding human muscle-type α-, δ-, and ε-subunits, and an additional cDNA encoding βAnc at a cDNA ratio of 2:1:1:1 (α:βAnc:δ:ε). (B) Same as in (A), but where cells were transfected with an altered 2:0.25:1:1 cDNA ratio, making the βAnc subunit limiting, or (C) where cells were transfected with only the cDNA encoding βAnc. In all cases openings are upward deflections, in the presence of 30 µM acetylcholine, and with an applied voltage of –120 mV. Continuous recordings are digitally filtered to 5 kHz, and the scale bar (2 s, 10 pA) in (A) applies to (B) and (C). Insets are digitally filtered to 10 kHz, with boxes representing scale bars (300 ms, 25 pA). (D) Event-based amplitude histograms for single-channel bursts from each of the patches shown in (A), (B), and (C). In each case, the height of the bins was normalised to the total number of bursts in each patch (A: 40; B: 50; C: 41), with the scale bar representing the indicated fraction (0.20) of the total bursts.

The traces in Figure 2 were recorded in the presence of agonist (30 μM acetylcholine). In the human adult muscle-type AChR, the agonist-binding sites are located at the α–δ and α–ε interfaces, and the β-subunit is the only subunit that does not participate directly in agonist binding (Figure 1; Rahman et al., 2022; Rahman et al., 2020; Zarkadas et al., 2022). We were therefore surprised to see single-channel activity in patches from cells transfected exclusively with βAnc, as channels formed from a muscle-type β-subunit alone would not be expected to have intact agonist-binding sites. To determine if the activity of βAnc-alone channels was dependent upon acetylcholine, we recorded single-channel activity in the absence of acetylcholine (Figure 3). When no acetylcholine was present, patches from cells that were transfected exclusively with βAnc cDNA still displayed single-channel activity, indicating that βAnc-alone channels open spontaneously under these conditions (Figure 3A). Furthermore, spontaneous activity occurred as bursts of closely spaced openings, separated by brief closings (Figure 3B; Colquhoun and Hawkes, 1982). The briefest of these intervening closings were reminiscent of classic ‘nachschlag shuttings’ (Figure 3B, ‘i’ in inset), observed in early patch clamp recordings from frog end-plate nicotinic receptors, and originally thought to relate to agonist efficacy (Colquhoun and Sakmann, 1981). Thus, despite being homomeric and lacking agonist-binding sites, βAnc-alone channels display single-channel hallmarks of the muscle-type AChR.

Figure 3 with 1 supplement see all
Spontaneous single-channel openings of βAnc homomers.

(A) Representative continuous recording of a cell-attached patch from cells transfected with a single cDNA encoding βAnc. Recording was made in the absence of acetylcholine and at an applied voltage of –120 mV. Data was digitally filtered to 5 kHz (scale bar = 2 s, 10 pA). (B) Single burst of openings from a homomeric βAnc channel, shown digitally filtered to 10 kHz (scale bar = 25 ms, 10 pA). Inset depicts (i) brief and (ii) long closings within bursts, where the former (i) are reminiscent of ‘nachschlag shuttings’ (scale bar = 1 ms, 5 pA). (C) Open and closed dwell duration histograms for the representative patch depicted in (B). Individual exponential components determined manually (dashed lines) and kinetic fits from MIL (solid lines) are overlaid. Global kinetic fitting was performed on three individual recordings, from two separate transfections. (D) The single-channel data fit a three-state scheme (Scheme 1), where C, C′, and O′ correspond to closed, closed-primed, and open-primed states. Rate constants with units s–1 are shown above and below corresponding arrow, with error estimates provided in Table 1.

Figure 3—source data 1

Source data for Figure 3.

Detected single-channel event durations of spontaneously opening βAnc homomers. Compressed file includes three TAC 4.3.3 event files (*.evt format) of the single-channel detections for the three recordings used in the presented kinetic analysis, as well as the associated R scripts (*.txt format) for defining and sorting bursts.

https://cdn.elifesciences.org/articles/76504/elife-76504-fig3-data1-v2.zip
Figure 3—source data 2

Unrasterized version of Figure 3.

https://cdn.elifesciences.org/articles/76504/elife-76504-fig3-data2-v2.eps

To gain insight into the spontaneous activity of βAnc-alone channels, we performed kinetic analysis of our single-channel data. First, we determined a critical closed duration (τcrit) to define bursts arising from a single ion channel. Then we determined the minimum number of components in our apparent open and closed dwell duration histograms by fitting each with a sum of exponentials. Open duration histograms were fit well by a single exponential component, while closed duration histograms required at least two components (Figure 3C). This suggested that a minimal scheme with a single open state and two closed states is necessary and sufficient to describe the spontaneous activity of βAnc-alone channels. Based on this, we then fit the sequence of single-channel dwells using the three possible kinetic schemes, two linear and one cyclic, comprising a single open state and two closed states (Figure 3D; Figure 3—figure supplement 1). As a control, we also fit a simplified two-state scheme, where a single open state was connected to a single closed state, which, based on the relatively poor fit of the closed durations, confirmed that inclusion of a second closed state was justified (Figure 3—figure supplement 1). Overlaying the resulting fits on top of duration histograms revealed that each of the possible three-state schemes fits the observed dwells equally well, thus discriminating between the possible kinetic schemes is not trivial (Figure 3—figure supplement 1). We settled upon the simple linear scheme, where βAnc-alone channels transition from a closed state, C, to an intermediate closed state, C′, before opening to O′. The form of this scheme, with an intermediate closed state that precedes channel opening, is guided by models of AChR activation that include a single ‘flipping’ or multiple ‘priming’ steps (Lape et al., 2008; Mukhtasimova et al., 2016; Mukhtasimova et al., 2009). Given that for βAnc-alone channels there is a single intermediate closed state that precedes channel opening, we refer to this state in our scheme as ‘flipped’ or ‘primed’.

In the presence of acetylcholine, the single-channel current traces appeared different (Figure 4A). As the concentration of acetylcholine increased from 10 to 100 μM, there was a progressive decrease in the duration of openings, as well as a reciprocal increase in the number of short-lived closings within each burst (Figure 4A). This can be observed as a leftward shift in the open duration histograms as the concentration of acetylcholine is increased. This single-channel behaviour is a hallmark of open-channel block, a ubiquitous property of AChR agonists, including acetylcholine (Lape et al., 2009; Mukhtasimova et al., 2016; Ogden and Colquhoun, 1985; Sine and Steinbach, 1984). Consistent with this blockage profile, the same trend, albeit with longer-lived blocking events, was observed with the well-characterized AChR open-channel blocker, 2-[(2,6-dimethylphenyl)amino]-N,N,N-trimethyl-2-oxoethaniminium chloride (QX-222) (Figure 4B; Charnet et al., 1990; Leonard et al., 1988; Pascual and Karlin, 1998).

Figure 4 with 2 supplements see all
Open-channel block of βAnc homomers by acetylcholine and 2-[(2,6-dimethylphenyl)amino]-N,N,N-trimethyl-2-oxoethaniminium chloride (QX-222).

Representative single-channel activity of βAnc homomers in the presence of increasing concentrations of (A) acetylcholine and (B) QX-222. Openings are upward deflections. Recordings were obtained with an applied voltage of –120 mV. Data were filtered to 10 kHz (scale bars = 25 ms, 10 pA; applies to (A) and (B)). The sequence of dwells from each dataset, encompassing the full concentration range of the blocker, was globally fit to the same three-state scheme used for βAnc, where an additional fourth state corresponding to the open/blocked channel was added (Scheme 2). Global kinetic fits were performed on three individual recordings for each concentration of blocker, from at least two separate transfections, corresponding to 15 total patches for each global fit. Note that the recordings in the absence of blocker are the same for each dataset. Rate constants are overlaid on the scheme below each dataset, with error estimates presented in Table 1.

Figure 4—source data 1

Source data for Figure 4.

Detected single-channel event durations of spontaneously opening βAnc homomers in the absence (agonist-free) and presence of increasing concentrations of acetylcholine (ACh; 10 µM, 30 µM, 60 µM, and 100 µM) or 2-[(2,6-dimethylphenyl)amino]-N,N,N-trimethyl-2-oxoethaniminium chloride (QX-222; QX; 10 µM, 30 µM, 60 µM, and 100µM). Compressed file includes 27 TAC 4.3.3 event files (*.evt format) of the single-channel detections for the three recordings for each condition (fileA, fileB, and fileC in each case) in the presented kinetic analysis, as well as the associated R scripts (*.txt format) for defining and sorting bursts.

https://cdn.elifesciences.org/articles/76504/elife-76504-fig4-data1-v2.zip
Figure 4—source data 2

Unrasterized version of Figure 4.

https://cdn.elifesciences.org/articles/76504/elife-76504-fig4-data2-v2.eps

The kinetics of open-channel block are determined by interactions between the blocking molecule and residues that line the channel pore in the open state, and therefore provide indirect structural insight into the open state. To compare the open state structures of βAnc-alone homomers with wild-type AChRs, we determined the kinetics of acetylcholine and QX-222 block for both types of channels (Figure 4; Figure 4—figure supplements 1 and 2). To fit our βAnc-alone single-channel data recorded in the presence of a blocker, we introduced an additional open, but blocked (i.e. non-conducting) state connected to our open state, where the forward rate of block was dependent upon the concentration of the blocking molecule (Figure 4A and B; Neher and Steinbach, 1978). We then globally fit each of our βAnc-alone datasets encompassing between 0 and 100 μM acetylcholine or QX-222. Initially, we restricted the rates of the core (C-C′-O′) scheme to those inferred in the absence of blocker; however, allowing all parameters to be estimated, led to negligible changes in the inferred rates of block. For βAnc-alone and wild-type channels, the rates of acetylcholine and QX-222 block were comparable (see Table 1 and Table 2). Of note, while the forward rate of QX-222 block (k+B) was almost the same for both types of channels, the reverse rate of QX-222 unblocking (k–B) was nearly twofold faster for βAnc-alone channels, indicating that the open pore of βAnc-alone channels has a slightly reduced affinity for QX-222. Regardless of this nuance, the similarity in the profiles of acetylcholine and QX-222 block suggests that the structure of the open pore in the two types of channels is similar.

Table 1
Single-channel kinetics of spontaneously opening βAnc homomers.
Homomerk+1k-1Kβ1α1Θ1k+Bk-BKB(μM)
No agonist
(3 patches)
3900
(110)
7600
(450)
0.5110,000
(300)
310
(3)
32.26N/AN/AN/A
Acetylcholine
(15 patches)
3500
(85)
5400
(265)
0.658000
(200)
330
(3)
24.24170*
(1.8)
62,000
(450)
364.71
QX-222
(15 patches)
4000
(110)
8000
(440)
0.5010,500
(300)
320
(3)
32.8195*
(0.5)
1400
(7)
14.74
Acetylcholine (Constrained)
(15 patches)
390076000.5110,00031032.26170*
(1.5)
60,800
(410)
357.65
QX-222 (Constrained)
(15 patches)
390076000.5110,00031032.2695*
(0.5)
1400
(6.5)
14.74
  1. Note: Rate constants were estimated from fitting Scheme 1 or Scheme 2 presented in Figures 3 and 4, respectively. Data were globally fit (number of patches indicated in each case) over a range of acetylcholine/2-[(2,6-dimethylphenyl)amino]-N,N,N-trimethyl-2-oxoethaniminium chloride (QX-222) concentrations, with rate constants and associated errors (parentheses) estimated by MIL (see Materials and methods). Priming (K′), gating (θ1), and blocking (KB) equilibrium constants represent k+1′/k-1′, β1/α1, and k-B/k+B, respectively. Association rate constants (*); (k+B) are presented in units of μM–1·s–1, while remaining rate constants are presented in units of s–1. Constrained rates (presented in italicized) were held constant at the values determined from the agonist-free dataset, while the blocking rate constants were estimated.

Table 2
Kinetics of acetylcholine (ACh) activation and 2-[(2,6-dimethylphenyl)amino]-N,N,N-trimethyl-2-oxoethaniminium chloride (QX-222) block of human adult muscle-type acetylcholine receptors (AChRs).
WTk+1k-1K1(μM)k+2k-2K2(μM)β1α1Θ1β2α2Θ2k+B (ACh)k-B(ACh)KB(ACh)(μM)k+B(QX-222)k-B(QX-222)KB(QX-222)(μM)
ACh
(24 pat.)
650*
(35)
14,400
(1000)
22.12325*
(20)
26,500
(400)
81.5433
(2.5)
8750
(650)
3.77 E-0314,000
(450)
1000
(10)
14215*
(3)
110,000
(800)
511.63N/AN/AN/A
QX-222 (10μM ACh)
(15 pat.)
650014,4002.21325026,5008.153387503.77 E-0314,0001000142150110,00051.16100*
(0.5)
845
(5.5)
8.45
QX-222 (30μM ACh)
(15 pat.)
19,50014,4000.74975026,5002.723387503.77 E-0314,0001000146450110,00017.05110*
(0.5)
845
(4.5)
7.68
  1. Note: Rate constants were estimated from fitting Scheme 3 or Scheme 4 presented in Figure 4—figure supplements 1 and 2, respectively. Where ‘A’ represents agonist, and ‘R’, ‘R*’, and ‘R*B’ represent the closed, open, and open-blocked states of the human adult muscle-type acetylcholine receptor (AChR). Data were globally fit (number of patches indicated in each case) over a range of 2-[(2,6-dimethylphenyl)amino]-N,N,N-trimethyl-2-oxoethaniminium chloride (QX-222) concentrations with fixed concentration of acetylcholine (ACh: 10 μM or 30 μM), with rate constants and associated errors (parentheses) estimated with MIL (see Materials and methods). Apparent binding (Kn), apparent gating (θn), and blocking (KB) equilibrium constants represent k-n/k+n, βn/αn, and k-B/k+B, respectively. Association rate constants (k+1, k+2, and k+B) are presented in units of μM–1·s–1, while remaining rate constants are presented in units of s–1. Constrained rates (italicized) were held constant to the rates estimated from fitting of wild-type (WT) in the absence of QX-222, allowing the QX-222 blocking rates to be estimated independently.

Given that βAnc-alone channels are expressed in the absence of other AChR subunits, a reasonable hypothesis is that they are homopentamers. To determine the subunit stoichiometry of βAnc-alone channels, we employed a single-channel electrical fingerprinting strategy, where mutations altering unitary conductance were used to count the number of individual βAnc subunits in βAnc-alone channels. A similar strategy has been employed with tetrameric potassium channels (Niu and Magleby, 2002), and other pLGICs (Andersen et al., 2011), including both the homopentameric α7 AChR (Andersen et al., 2013; daCosta et al., 2015; daCosta and Sine, 2013) and the heteropentameric muscle-type AChR (Emlaw et al., 2021). The approach relies on identifying high-conductance (HC) and low-conductance (LC) variants of the βAnc subunit, and then co-expressing them to reveal a number of amplitude classes. Openings in each amplitude class originate from channels incorporating the same ratio of HC to LC subunits, and based on the total number of amplitude classes, the number of βAnc subunits within βAnc-alone channels can be inferred.

When βAnc is expressed alone, the resulting channels exhibit a single, uniform amplitude, distributed around a mean of ~16 pA (Figure 5A and D), making the wild-type βAnc subunit an ideal HC subunit for electrical fingerprinting. To identify a LC variant of βAnc, we took advantage of a structural feature inherent to eukaryotic pLGICs: as conducting ions exit the channel’s transmembrane pore, they are obliged to pass through one of five portals in the cytoplasmic domain (Rahman et al., 2020; Unwin, 2005). Framed by charged or polar residues from each subunit, these portals influence single-channel conductance. Mimicking the homologous 5-HT3A receptor, which has an unusually low single-channel conductance (Hales et al., 2006; Kelley et al., 2003; Peters et al., 2004), we substituted three arginine residues (E420R, D424R, and E428R) into this region of βAnc. When βAnc harbouring three arginines in this region was transfected by itself, the resulting channels exhibited a reduced single-channel amplitude centred around~1–2pA (Figure 5B and D). With its markedly reduced amplitude, βAnc harbouring three arginine residues is a suitable LC subunit for electrical fingerprinting.

Figure 5 with 1 supplement see all
Electrical fingerprinting to determine subunit stoichiometry of βAnc homomers.

Representative single-channel activity from cells transfected with (A) cDNA encoding the wild-type high-conductance (HC) βAnc subunit, or (B) a mutant low-conductance (LC) βAnc variant harbouring substitutions that reduce single-channel amplitude. (C) Cotransfection of cDNAs encoding HC and LC βAnc variants led to patches (two shown) with heterogeneous amplitudes. (D) The amplitudes segregate into six well-defined amplitude classes (total of 495 bursts combined from the two patches in (C)), where the highest and lowest amplitude classes match that of the all-HC (1569 bursts) and all-LC classes (883 bursts), respectively. Plot of the mean amplitude of each class as a function of the presumed number of incorporated HC subunits (error bars = standard deviations of the mean but are smaller than the points themselves). Recordings were obtained with an applied voltage of –120 mV, and traces were digitally filtered to 1 kHz to facilitate amplitude detection (scale bar = 50 ms, 5 pA; applies to (A), (B), and (C)).

When cDNAs encoding HC and LC variants of βAnc were transfected together, a variety of single-channel amplitudes were observed in each patch (Figure 5C). The relative proportion of channels with high versus low amplitude could, to some degree, be tuned by the ratio of HC to LC βAnc cDNA used for transfection (Figure 5—figure supplement 1). Constructing event-based amplitude histograms, and pooling amplitudes from more than one recording, revealed that the amplitudes segregated into as many as six amplitude classes, with the highest and lowest amplitude classes matching that of the HC and LC forms of βAnc-alone channels. The difference in amplitude between successive classes was somewhat regular, demonstrating five approximately equal contributions to single-channel conductance (Figure 5D), consistent with the hypothesis that βAnc-alone channels are homopentamers.

As noted previously, reconstruction of βAnc was based upon a molecular phylogeny that diverged from the accepted species phylogeny (Emlaw et al., 2021; Emlaw et al., 2020; Prinston et al., 2017; Tessier et al., 2017). Reconstruction of ancestral protein sequences is based upon a best-fit model of amino acid evolution, a multiple sequence alignment, as well as a phylogenetic tree relating the sequences within the alignment (Thornton, 2004). We therefore wondered if the ability of βAnc to form spontaneously opening homomers was an artefact of the discordant tree used to reconstruct it. To test this, we took advantage of an alternate ancestral β-subunit, called ‘βAncS’, whose reconstruction was based upon a molecular phylogeny that matched the accepted species phylogeny (Emlaw et al., 2021). Despite 67 substitutions and 6 indels relative to βAnc, when expressed alone, βAncS still formed homomeric channels that opened in bursts in the absence of acetylcholine (Figure 6). This demonstrated that the ability of βAnc to form homomers that spontaneously open is not an artefact of the phylogeny used to reconstruct it. Instead, this surprising ability of βAnc is robust to the phylogenetic uncertainties inherent in ancestral sequence reconstruction, as well as substantial variation in the amino acid sequence of the reconstructed ancestral β-subunits.

Spontaneous single-channel openings of homomers formed from an alternate ancestral β-subunit (βAncS).

(A) Representative continuous recording of a cell-attached patch from cells transfected with a single cDNA encoding βAncS. Recording was made in the absence of acetylcholine, at an applied voltage of –120 mV, where spontaneous openings are upward deflections. Data was digitally filtered to 5 kHz (scale bar = 2 s, 10 pA). (B) Bursts from homomeric βAncS channels, each exhibiting one of three different types (i, ii, iii) of openings (scale bar = 25 ms, 10 pA). The boxed burst at the bottom is an example of a single burst that contains more than one type of opening (ii/i). (C) Open and closed dwell duration histograms for the representative patch depicted in (B). Individual exponential components determined manually (dashed lines) and kinetic fits from MIL (solid lines) are overlaid. Global kinetic fitting was performed on three individual recordings, from two separate transfections. The exponential components (i, ii, iii) in the open duration histogram correspond to the different types of openings observed within the bursts in panel (B). (D) The three-state scheme in Figure 3 (Scheme 1) was expanded to include additional priming steps (‘singly’, ‘doubly’, and ‘triply’ primed), each with their own connected open state (Scheme 5). Rate constants are shown, with error estimates provided in Table 3.

Figure 6—source data 1

Source data for Figure 6.

Detected single-channel event durations of spontaneously opening βAncS homomers. Compressed file includes three TAC 4.3.3 event files (*.evt format) of the single-channel detections for the three recordings used in the presented kinetic analysis, as well as the associated R scripts (*.txt format) for defining and sorting bursts.

https://cdn.elifesciences.org/articles/76504/elife-76504-fig6-data1-v2.zip
Figure 6—source data 2

Unrasterized version of Figure 6.

https://cdn.elifesciences.org/articles/76504/elife-76504-fig6-data2-v2.eps

While βAncS forms homomers that open spontaneously, inspection of the βAncS single-channel activity revealed additional complexity not seen with βAnc. For βAncS, single-channel bursts were heterogeneous, displaying at least three distinct kinetic behaviours (Figure 6B). In some cases, the kinetic behaviour changed within a burst, demonstrating that the different kinetics were possible within the same channel (Figure 6B; boxed). This heterogeneity was also reflected in apparent open and closed duration histograms, with each displaying a minimum of three exponential components (Figure 6C). The increased number of exponential components indicated additional open and closed states relative to βAnc, and thus that the three-state scheme used to fit βAnc was insufficient to describe the spontaneous activity of βAncS. To account for the additional states, we expanded our original scheme to include two additional ‘priming’ steps, where openings could occur from one of three primed states (Figure 6D). This scheme, with multiple priming steps, builds directly upon the one used to fit muscle-type AChRs that had been engineered to open spontaneously (Mukhtasimova et al., 2009). Given the heterogeneity of the βAncS single-channel activity, and the complexity of this scheme, we caution against over interpretation of the inferred rates. Nevertheless, we note that in accord with the muscle-type AChR, the equilibrium gating constants appear to increase (Table 3; compare Θ1, Θ2, and Θ3), and thus the open states become more and more favoured, for each successive priming step. In any case, the fits suggest that a scheme of this form, with multiple stages of priming, is adequate to describe the complex spontaneous single-channel activity of βAncS homomers.

Table 3
Single-channel kinetics of spontaneously opening βAncS homomers.
k+1k-1β1α1Θ1k+2′′k-2′′β2α2Θ2k+3′′′k-3′′′β3α3Θ3
βAncS
(3 pat.)
2000
(80)
2500
(140)
11,500
(450)
17,000
(800)
0.67613,650
(400)
14,800
(2000)
175,000
(17,000)
6800
(1400)
25.742000
(300)
5000
(450)
40,000
(1550)
500
(15)
80
  1. Note: Rate constants were estimated from fitting Scheme 5 in Figure 6 (and below). Data were globally fit (three individual patches from two separate transfections) with rate constants and associated errors (parentheses) estimated within MIL (see Materials and methods). Gating equilibrium (θn) constants represent βn/αn. Rate constants are presented as s–1.

Discussion

The goal of this work was to identify the source of heterogeneity in single-channel recordings from cells transfected with our previously described, reconstructed ancestral AChR β-subunit (Prinston et al., 2017). The heterogeneity in single-channel activity was only evident upon modification of our original transfection protocol. Specifically, when the total amount of transfected AChR subunit cDNA was reduced ~sixfold, a second class of single-channel openings with increased amplitude and a different kinetic signature was observed. Evidently, βAnc can both incorporate into heteropentameric AChRs (Emlaw et al., 2021), and self-associate to form homopentameric βAnc channels. Reducing the overall expression of AChR subunits was sufficient to tip the balance between the incorporation of βAnc into hetero- versus homopentamers, and thus unmask this previously unobserved ability of βAnc to form homomeric channels.

At first glance, the ability of βAnc to form homopentamers is unexpected. Modern-day muscle-type β-subunits do not appear to form homopentamers, and instead appear fully entrenched within heteropentameric muscle-type AChRs. Since reconstruction of βAnc was informed almost exclusively by modern muscle-type β-subunits, there was no reason to expect that βAnc would behave differently (Prinston et al., 2017). However, as mentioned previously, βAnc can replace both the human muscle-type β- and δ-subunits in hybrid ancestral/human AChRs (Emlaw et al., 2021). For this to be possible the principal (+) and complementary (−) subunit interfaces of βAnc must be compatible with each other (Figure 1), leading to the logical hypothesis, confirmed here, that βAnc can form homopentamers. From an evolutionary perspective, this ability to revert the muscle-type β-subunit, which is entrenched in a heteropentamer, back to a subunit capable of forming homopentamers, attests to the presumed homopentameric origins of pentameric ligand-gated ion channel subunits (Ortells and Lunt, 1995).

We have also shown that βAnc homopentamers open spontaneously. This is surprising. What is even more surprising is that the spontaneous single-channel activity of βAnc homopentamers resembles that of the agonist-activated muscle-type AChR. Some of the first single-channel recordings of frog AChRs activated by agonists revealed that agonist-activated openings occur in quick succession, as bursts of activity originating from the same channel (Sakmann et al., 1980). The main effect of increasing the concentration of agonist was to increase the open probability within bursts, with saturating concentrations of full agonists, such as acetylcholine, leading to bursts where the open probability exceeded 0.90 (i.e. the channel was open for more than 90% of the duration of the burst). Spontaneous openings of βAnc homopentamers also occur in bursts, and despite the absence of agonist, activation appears efficient, with the probability of being open within a burst also exceeding 0.90. The kinetic structure of bursts of spontaneous βAnc openings also resembles that of the agonist-activated AChR, with bursts containing several types of closings, the briefest of which are reminiscent of classic ‘nachschlag shuttings’. This is significant because it was originally proposed that the duration of nachschlag shuttings was related to agonist efficacy (Colquhoun and Sakmann, 1981). However, subsequent work showed that nachschlag shuttings were independent of agonist, which necessitated the introduction of an additional closed state appended to the agonist-activated open state in early schemes of AChR activation (Sine and Steinbach, 1986). This latter finding was some of the impetus for refined schemes, where the additional closed states preceded channel opening and were referred to as ‘flipped’ or ‘primed’ (Lape et al., 2008; Mukhtasimova et al., 2009). Our kinetic analysis has shown that the spontaneous activity of βAnc fits an analogous scheme, containing an intermediate closed state that also precedes channel opening, but where the agonist binding steps have been omitted due to the absence of agonist. Evidently, these functional hallmarks of AChR activation do not arise from the complex heteropentameric architecture of the muscle-type receptor, nor do they depend upon the presence of agonist. Instead, they are fundamental properties preserved and encoded in the reconstructed ancestral amino acid sequence of βAnc (and βAncS).

Relative to the human β-subunit, which does not form spontaneously opening homopentamers, βAnc contains 125 individual substitutions and 7 deletions (Prinston et al., 2017). These 132 amino acid differences are scattered throughout the entire β-subunit, and presumably not all of them are required for βAnc to form spontaneously opening homopentamers. Thus, although the differences in amino acid sequence between βAnc and the human β-subunit are sufficient to convert the human β-subunit into a subunit capable of (1) forming homopentamers and (2) spontaneously opening, the large number of residues involved, as well as their delocalised nature, obscures the subset of residues strictly necessary for either (or both) functions. Nevertheless, several of the residues substituted between the human β-subunit and βAnc align with known determinants of homomeric assembly and spontaneous activity in various pLGICs (Appendix 1—figure 1).

Compared to the human wild-type and βAnc-containing AChRs, βAnc homopentamers exhibit increased single-channel amplitude and thus conductance (Figure 2D; Figure 5; Figure 2—figure supplement 1). This increased conductance can be explained by amino acid differences in the second transmembrane segment (i.e. M2) of each subunit, which is one of the more conserved regions across all pLGICs. M2 lines the channel pore at its narrowest constriction and is a well-known determinant of single-channel conductance (Imoto et al., 1986). With exception of the 2′ and 6′ positions, M2 is conserved between βAnc and the human β-subunit (Emlaw et al., 2020). However, to reconcile the differences in conductance, it is important to consider the M2 segments of not just the β-subunits, but all the subunits that make up the various channels. Owing to the symmetric (homopentamers) or pseudosymmetric (heteropentamers) architecture of pLGICs, analogous M2 residues from each subunit occupy similar positions in the pore, forming concentric rings of amino acids. Several of these rings influence conductance, with the –4′, –1′, 2′, and 20′ positions being particularly important, and referred to as the ‘cytoplasmic’, ‘intermediate’, ‘central’, and ‘extracellular’ rings, respectively (Imoto et al., 1991; Imoto et al., 1988). In βAnc homopentamers, with five identical βAnc subunits, the amino acids contributed from individual subunits are identical. This is significant because βAnc contains acidic residues at the –4′ (Asp), –1′ (Glu), and 20′ (Asp) positions, meaning that all 15 residues forming these three decisive rings are acidic, effectively maximising their negative potential (Appendix 1—figure 2). In heteropentamers, such as human wild-type and βAnc-containing AChRs, the residues contributed by individual subunits are in some cases different, with polar/neutral (Gln) and basic (Lys) residues contributed by the ε- and δ-subunits, respectively (Appendix 1—figure 2). Substitution of neutral or basic residues into any of these three rings has been shown to reduce conductance (Imoto et al., 1988). Thus, unlike for βAnc homopentamers that have higher conductance, the negative potential of these rings is not maximised in human wild-type and βAnc-containing heteropentamers.

Based upon the rates of open-channel block by acetylcholine and the well-known open-channel blocker, QX-222, the structure of the open pore in βAnc homopentamers and the wild-type AChR appear to be similar. Nevertheless, kinetic fitting revealed that compared to wild-type, βAnc homopentamers exhibit an almost twofold reduced affinity for QX-222. While we have employed the simplest model for open-channel block (Neher and Steinbach, 1978), which neglects potential increased complexity at higher QX-222 concentrations (Neher, 1983), the simplest interpretation is that, like conductance, this reduced apparent affinity can be explained by amino acid differences in M2. Specifically, the 6′ and 10′ sites are known determinants of QX-222 block (Charnet et al., 1990; Leonard et al., 1988; Pascual and Karlin, 1998), with a more polar environment at 6′, and a less polar environment at 10′, favouring interaction with QX-222 (Charnet et al., 1990). In the wild-type AChR, the polarity of the 6′ ring is limited by a hydrophobic phenylalanine residue contributed by the human β-subunit (Appendix 1—figure 2). In βAnc homopentamers, the 6′ ring does not contain this phenylalanine residue and is instead composed of five polar serine residues, predicting an increased affinity for QX-222. Conversely, the 10′ ring in βAnc homopentamers contains five threonine residues making it more polar, and thus a less favourable environment for QX-222, than in the wild-type AChR, which contains two nonpolar alanine residues (from the δ- and ε-subunits) (Appendix 1—figure 2). Evidently, the energetic contributions of these two rings to QX-222 blockage balance out such that βAnc homopentamers exhibit only a slightly reduced apparent affinity for QX-222 compared to the wild-type AChR.

Phylogenetic analyses of pLGICs suggests that present-day heteropentameric subunits, such as those entrenched in the muscle-type AChR, evolved from ancestral subunits capable of self-associating to form homopentamers (Dent, 2010; Jaiteh et al., 2016; Le Novère and Changeux, 1995; Ortells, 2016; Ortells and Lunt, 1995; Pedersen et al., 2019; Tsunoyama and Gojobori, 1998). Within this framework, strictly entrenched subunits have lost compatibility between their (+) and (–) interfaces, and thus the ability to self-associate. This is significant because the pentameric architecture of pLGICs places a constraint upon when the ability to self-associate became dispensable, and thus when strict subunit entrenchment could have evolved. For the simplest heteropentamers, made up of two types of subunits, at least one of the subunit types will have to occupy a minimum of three positions within the pentamer. In these cases, two copies of the same subunit must sit next to each other, and thus retain the ability to self-associate. Indeed, in neuronal α4β2 heteropentameric AChRs, both the α4- and β2-subunits have retained the ability to self-associate and form α4–α4 or β2–β2 interfaces, meaning that in principle, a variety of subunit stoichiometries and arrangements are possible. While additional factors such as subunit packing, relative expression, trafficking and assembly, and the presence of cellular chaperones provide additional constraints that limit the stoichiometries and arrangements observed (Mazzaferro et al., 2021; Moroni et al., 2006; Nelson et al., 2003; Son et al., 2009; Walsh et al., 2018), the requirement for self-association must remain for at least one of the subunit types. Thus, while the ability to self-associate is not sufficient to permit all possible heteropentameric stoichiometries and arrangements, it is necessary.

Unlike neuronal heteropentameric channels, such as α4β2 above, whose subunit stoichiometries and arrangements exhibit a degree of plasticity, the subunits in the muscle-type AChR are strictly entrenched. In this case, the definition of strict entrenchment neglects the γ- for ε-subunit swap that occurs during development (Mishina et al., 1986), because it occurs at the same location within the heteropentamer, and thus does not alter the stoichiometry and arrangement of the α-, β-, or δ-subunits. The above structural considerations imply that such strict entrenchment should only be expected in heteropentamers composed of three or more subunit types, and where each of the subunits has also lost the ability to self-associate. Muscle-type AChRs fulfill both criteria, in that they are made up of four types of subunits that are each incapable of self-associating and forming homopentamers. The finding that reconstructed ancestral β-subunits retain the ability to self-associate and form homomers, suggests that, unlike their extant human and Torpedo counterparts, the ancestral β-subunits were not strictly entrenched, and that their entrenchment occurred in parallel lineages after the last common ancestor shared by humans and Torpediniformes (i.e. Gnathostomata).

A degree of uncertainty is inherent in the reconstruction of any ancestral protein, and to solidify evolutionary conclusions it is important to assess whether the functions of putative ancestors are robust to these uncertainties. We have shown that an alternate ancestral β-subunit (βAncS), with 67 substitutions and 6 indels relative to βAnc, was still able to form homomers that opened spontaneously. Furthermore, spontaneous openings of βAncS homomers also occurred in bursts that contained brief closings (i.e. ‘nachschlag shuttings’) and had an open probability that exceeded 0.90. These findings demonstrate that the ability of βAnc to form spontaneously opening homomers is robust to substantial variations in the inferred ancestral β-subunit amino acid sequence. Thus, the capacity of these reconstructed ancestral β-subunits to form spontaneously opening homomers is deeply embedded in their structure and evolutionary history.

The kinetic behaviour of βAncS homomers was more complex than originally observed with βAnc. This increased complexity necessitated the expansion of the original scheme used to fit βAnc to include two additional priming steps. Both the original βAnc and expanded βAncS schemes parallel the scheme used to describe the muscle-type AChR, which included two priming steps (i.e. ‘singly’ and ‘doubly’ primed) that each correlated with conformational changes around the two agonist-binding sites (Mukhtasimova et al., 2009). In the case of βAncS, which lacks agonist-binding sites, the simplest interpretation is that the different levels of priming correlate with conformational changes occurring within individual subunits. Although the states in Scheme 5 are labelled as ‘singly’, ‘doubly’, and ‘triply’ primed, the three priming steps could represent the three most terminal priming steps, where (assuming that βAncS also forms homopentamers) three, four, or five βAncS subunits are primed. Within this framework, openings from βAncS channels with zero, one, or two primed subunits are presumably unstable, and thus not observed in our single-channel recordings. In an alternate scenario, ‘singly’ and ‘doubly’ primed could refer to βAncS channels with one or two primed subunits, respectively. While ‘triply’ primed could refer to channels with three or more primed subunits, but where openings from βAncS channels with three, four, or five primed subunits are indistinguishable. Applied to βAnc, these interpretations suggest that either (1) openings from only the terminal priming step, where all five βAnc subunits are primed, are stable enough to be observed in βAnc homopentamers, or (2) openings can occur with fewer than five primed subunits, but these openings are kinetically indistinguishable. Regardless of the interpretation, priming is an important step in the activation of both βAncS and βAnc homopentamers.

Modern mechanisms of AChR activation include intermediate closed states, which place the roots of agonism at a stage in the activation process that precedes channel opening (Lape et al., 2008; Mukhtasimova et al., 2016; Mukhtasimova et al., 2009). A consequence of these mechanisms is that the ultimate opening and closing rates of the channel are independent of agonist. Here, we have shown that additional single-channel hallmarks of AChR function are independent of agonist, as they occur in spontaneously opening homopentameric channels that are devoid of agonist-binding sites, and which are formed from reconstructed ancestral AChR β-subunits. Often overlooked, the β-subunit is the least conserved of the four AChR subunits, and is the only subunit that does not contribute residues to the two AChR agonist-binding sites. Despite these considerations, hallmarks of AChR function remain deeply embedded in β-subunit sequence, structure, and evolutionary history. Given that these functional hallmarks are independent of agonist, it is tempting to speculate that they predate agonism, and thus that agonism evolved subsequently as an additional layer of regulation in this family of pentameric ion channels.

Materials and methods

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Cell line (Homo sapiens)BOSC 23ATCCCRL11270
(discontinued)
Modified Homo sapiens embryonic kidney cells
Recombinant DNA reagentpRBG4 – AChR α1Provided by Steven M. Sine (Mayo Clinic)Homo sapiens CHRNA1 (Accession: NM_000079.4)
Recombinant DNA reagentpRBG4 – AChR β1Provided by Steven M. Sine (Mayo Clinic)Homo sapiens CHRNB1
(Accession: NM_000747.3)
Recombinant DNA reagentpRBG4 – AChR δProvided by Steven M. Sine (Mayo Clinic)Homo sapiens CHRND
(Accession: NM_000751.3)
Recombinant DNA reagentpRBG4 – AChR εProvided by Steven M. Sine (Mayo Clinic)Homo sapiens CHRNE
(Accession: NM_000080.4)
Recombinant DNA reagentpRBG4 – AChR βAncCustom gene synthesisConstruct originating from: PMID:28689969
Recombinant DNA reagentpRBG4 – AChR βAncSCustom gene synthesisConstruct originating from: PMID:33579823
Recombinant DNA reagentpGreenLanternProvided by Steven M. Sine (Mayo Clinic)
Commercial assay or kitQ5 DNA polymeraseNew England Biolabs, incM0491PCR
Sequence-based reagentSDM_AChRBAncLC_FThis paperMutagenesis primerAGAACGCTGAAGAGAGACTGGCAGTACGTGGCCAT
Sequence-based reagentSDM_AChRBAncLC_RThis paperMutagenesis primerATAGTCCTCTCTTTTCTGCAGCTGCTCAGCGAT
Chemical compound, drugAcetylcholine ChlorideSigmaA9101-10VLPurity: 99%
Chemical compound, drugQX-222Tocris1043/10Purity:>98%
Software, algorithmTAC 4.3.3Bruxton (https://www.bruxton.com/legacy.html)Single-channel recording, detection, and analysis
Software, algorithmRhttps://www.r-project.org/Open-source statistical computing software
Software, algorithmscburstshttps://cran.r-project.org/web/packages/scbursts/index.htmlR Package – single-channel burst analysis
Software, algorithmextreme-valueshttps://cran.r-project.org/web/packages/extremevalues/index.htmlR Package – outlier detection
Software, algorithmMASShttps://cran.r-project.org/web/packages/MASS/index.htmlR Package – function and statistical analysis
Software, algorithmxlsxhttps://cran.r-project.org/web/packages/xlsx/index.htmlR Package – read and write excel files

Materials

The 2-[(2,6-dimethylphenyl)amino]-N,N,N-trimethyl-2-oxoethaniminium chloride (QX-222) was purchased from Tocris Bioscience. All other chemicals, including acetylcholine chloride, were purchased from Sigma-Aldrich.

Molecular Biology

Request a detailed protocol

cDNAs of human muscle-type AChR subunits (⍺1, β1, δ, and ε) in the pRBG4 plasmid were provided by Steven M. Sine (Mayo Clinic), while cDNAs encoding βAnc and βAncS were reconstructed and cloned into pRBG4 as described previously (Emlaw et al., 2021; Prinston et al., 2017). Mutations to produce the LC variant of βAnc (E420R, D424R, E428R) were introduced by inverse PCR (Silva et al., 2017). Sanger sequencing confirmed the entire reading frame for all constructs.

Mammalian cell expression

Request a detailed protocol

Combinations of human and ancestral subunit cDNAs were transfected into BOSC 23 cells (Pear et al., 1993), originally from ATCC (CRL11270), but provided by Steven M. Sine (Mayo Clinic) (RRID:CVCL_4401). Cells were maintained in Dulbecco’s modified Eagle’s medium (DMEM; Corning) containing 10% (vol/vol) fetal bovine serum (Gibco) at 37°C, until they reached 50–70% confluency. Cells were then transfected using calcium phosphate precipitation, and transfections terminated after 3–4 h by exchanging the medium. All experiments were performed one day post transfection (between 16 and 24 h after exchanging the medium). A separate plasmid encoding green fluorescent protein was included in all transfections to facilitate identification of transfected cells.

Cell line authentication and mycoplasma testing

Request a detailed protocol

Approximately 5 million confluent cells were harvested and their total DNA isolated (E.Z.N.A. Tissue DNA Kit), and then submitted to The Centre for Applied Genomics Genetic Analysis Facility (The Hospital for Sick Children, Toronto, Canada) for STR profiling using Promega’s GenePrint 24 System. A similarity search on the 8,159 human cell lines with STR profiles in Cellosaurus release 42.0 was conducted on the resulting STR profile, which revealed that the cell line shares closest identity (88%, CLASTR 1.4.4 STR Similarity Search Tool score) with cell line Anjou 65 (CVCL_3645). Anjou 65 is a child of CVCL_1926 (HEK293T/17) and is itself a parent line of CVCL_X852 (Bartlett 96). Bartlett 96 is the parent line of BOSC 23 (Pear et al., 1993). PCR tests confirmed that the cells were free from detectable mycoplasma contamination (Uphoff and Drexler, 2011, 2002).

Single-channel patch clamp recordings

Request a detailed protocol

Single-channel patch clamp recordings were performed as previously described (Mukhtasimova et al., 2016). Recordings from BOSC 23 cells transiently transfected with cDNAs encoding wild-type, ancestral, or LC subunits, were obtained in a cell-attached patch configuration. All recordings were obtained with a membrane potential of –120 mV, with room temperature maintained between 20 and 22°C. The external bath solution contained 142 mM KCl, 5.4 mM NaCl, 0.2 mM CaCl2, and 10mM 4-(2-hydroxyethyl)–1-piperazineethanesulfonic acid (HEPES), adjusted to pH 7.40 with KOH. The pipette solution contained 80 mM KF, 20 mM KCl, 40 mM K-aspartate, 2 mM MgCl2, 1 mM ethylene glycol-bis(β-aminoethyl ether)-N,N,N′,N′-tetraacetic acid, and 10 mM HEPES, adjusted to a pH of 7.40 with KOH. Acetylcholine and QX-222 were added to pipette solutions at their desired final concentrations and stored at –80°C. Patch pipettes were fabricated from type 7052 or 8250 non-filamented glass (King Precision Glass) with inner and outer diameters of 1.15 and 1.65 mm, respectively, and coated with SYLGARD 184 (Dow Corning). Prior to recording, electrodes were heat polished to yield a resistance of 5–8 MΩ. Single-channel currents were recorded on an Axopatch 200B patch clamp amplifier (Molecular Devices), with a gain of 100 mV/pA and an internal Bessel filter at 100 kHz. Data were sampled at 1.0 μs intervals using a BNC-2090 A/D converter with a National Instruments PCI 6111e acquisition card and recorded by the program Acquire (Bruxton).

Dwell time and kinetic analysis

Request a detailed protocol

Single-channel detections were performed using the program TAC 4.3.3 (Bruxton). Data were analysed with an applied 10 kHz digital Gaussian filter. Opening and closing transitions were detected using the 50% threshold crossing criterion, and open and closed dwell duration histograms were generated within the program TACfit 4.3.3 (Bruxton). Histograms were visually fit with a minimum sum of exponential components. From the closed duration histograms, the intersection of the slowest activation and fastest inactivation/desensitisation components was taken as the critical closed duration (τcrit) (Colquhoun and Hawkes, 1982; Sine et al., 1990). Closings longer than τcrit (corresponding to inter-burst closings) were removed from analysis. Events were imported into R using scbursts (Drummond et al., 2019), where individual durations were corrected for instrument risetime (Colquhoun and Sigworth, 1995), and bursts were defined by a τcrit (Sine et al., 1990). Bursts with fewer than three events were omitted from further analysis. The probability of being open within a burst (i.e. burst Popen) was calculated for each burst, and bursts with a Popen that did not fit within the normal distribution were removed using extremevalues (van der Loo, 2010). The distribution of burst Popen was then fit with a Gaussian distribution, and bursts within two standard deviations from the mean were used for further kinetic analysis (Drummond et al., 2019). A user-defined kinetic scheme (see Figure 3—figure supplement 1 for βAnc, and the modified del Castillo and Katz scheme (Colquhoun and Sakmann, 1985) for the human adult AChR in Figure 4—figure supplement 1) was fit to the sequence of single-channel dwells in the global dataset using maximum likelihood implemented within MIL (QUB suite, State University of New York, Buffalo, NY). With a user-defined dead time of 18.83 µs, MIL corrected for missed events, estimated model parameters by maximum likelihood, and gave standard errors of the estimated parameters (see Tables 1–3; Qin et al., 1996).

Electrical fingerprinting

Request a detailed protocol

The HC and LC variants of βAnc were transfected at 1:4 and 4:1 (HC:LC) cDNA ratios. Transfections and single-channel recordings were performed as described above. For detections, data were filtered to 1 kHz, and bursts defined by a uniform τcrit of 2 ms imposed upon all recordings. Using the program TAC 4.3.3 (Bruxton), amplitudes of single-channel bursts were measured as the difference between open- and closed-channel currents. Amplitudes of individual bursts were pooled from separate recordings to generate event-based amplitude histograms (EBAHs; Figure 5), which were fit with Gaussian distributions within TACfit (Bruxton).

Schemes

Request a detailed protocol
Kinetic scheme describing spontaneous single-channel activity of βAnc homomers.
Kinetic scheme describing single-channel activity of βAnc homomers in the presence of acetylcholine or 2-[(2,6-dimethylphenyl)amino]-N,N,N-trimethyl-2-oxoethaniminium chloride (QX-222).
Kinetic scheme describing single-channel activity of the human adult muscle-type acetylcholine (ACh) receptor.
Kinetic scheme describing single-channel activity of the human adult muscle-type acetylcholine (ACh) receptor in the presence of 2-[(2,6-dimethylphenyl)amino]-N,N,N-trimethyl-2-oxoethaniminium chloride (QX-222).
Kinetic scheme describing spontaneous single-channel activity of βAncS homomers.

Appendix 1

Multiple sequence alignments

Provided is the multiple sequence alignment (Supplementary file 1) used to identify residues that differ between βAnc, βAncS, and the human β-subunit, and that also align with residues shown to be important for homomeric pLGIC assembly (Alves et al., 2011; Hannan and Smart, 2018; Sexton et al., 2021; Taylor et al., 1999) and spontaneous activity (Beckstead, 2002; Engel et al., 1996; Grosman and Auerbach, 2000; Miller et al., 2008; Mukhtasimova et al., 2009; Nayak et al., 2012; Ohno et al., 1995; Pan et al., 1997; Purohit and Auerbach, 2009; Torres and Weiss, 2002; Appendix 1—figure 1). Also shown is a local alignment of the second transmembrane region (M2) of the adult human muscle-type subunits (α-, β-, δ-, and ε-subunits), as well as βAnc and βAncS (Appendix 1—figure 2) highlighting differences in well-characterized determinants of single-channel conductance (Imoto et al., 1991; Imoto et al., 1988) and open-channel block (Charnet et al., 1990; Leonard et al., 1988; Pascual and Karlin, 1998).

Appendix 1—figure 1
Multiple sequence alignment of βAnc and βAncS with the human β-subunit (βHs).

Positions aligning with those previously reported as being important for spontaneous openings (blue) and homomeric assembly (orange) in other pentameric ligand-gated ion channels (pLGICs) are highlighted. Sites implicated in both spontaneous opening and homomeric assembly are highlighted in yellow. In each case, if the residue implicated is conserved across all three subunits the highlight is faded. Corresponding references for each position are indicated (a–n), where a: Torres and Weiss, 2002; b: Beckstead, 2002; c: Miller et al., 2008; d: Nayak et al., 2012; e:Purohit and Auerbach, 2009; f: Taylor et al., 1999; g: Sexton et al., 2021; h: Alves et al., 2011; i: Hannan and Smart, 2018; j: Pan et al., 1997; k: Mukhtasimova et al., 2009; l: Ohno et al., 1995; m: Grosman and Auerbach, 2000; n: Engel et al., 1996.

Appendix 1—figure 2
Primary and tertiary structure of the M2 segments from the adult human muscle-type AChR subunits, as well the two reconstructed ancestral β-subunits (βAnc and βAncS).

(A) Multiple sequence alignment of residues comprising the pore-lining M2 helix, highlighting positions discussed in the main text and known to influence single-channel conductance (black boxes; –4′, –1′, 2′, and 20′), or open-channel block (grey boxes; 6′ and 10′). (B) Tertiary structure of the AChR pore. M2 segments corresponding to the residues shown in panel “A” are mapped onto the Cα backbone of the Torpedo AChR (PDB: 6UWZ). The two subunits in the front (α and ε) have been hidden for illustrative purposes, and positions of highlighted residues in “A” are mapped onto the β-subunit with their Cα shown as black and grey spheres. (C) Analogous residues from each subunit form concentric rings, where the chemical properties of amino acids in each ring influence single-channel conductance or open-channel block. Residues contributed to each ring by individual subunits in the heteropentameric human adult muscle-type AChR (left) and βAnc-containing hybrid AChR (middle), as well as in βAnc homompentamers (right) are shown. Acidic, basic, polar, and nonpolar residues are colored red, blue, green, and yellow, respectively.

Data availability

All data generated and analysed during this study are included in the manuscript and as source data.

References

  1. Book
    1. Colquhoun D
    2. Sigworth FJ
    (1995) Fitting and Statistical Analysis of Single-Channel Records
    In: Colquhoun D, editors. Single-Channel Recording. Springer. pp. 483–587.
    https://doi.org/10.1007/978-1-4419-1229-9
  2. Book
    1. Hille B
    (2001)
    Ion Channels of Excitable Membranes
    Sinauer Associates, Inc.
  3. Book
    1. Silva D
    2. Santos G
    3. Barroca M
    4. Collins T
    (2017) Inverse PCR for point mutation introduction
    In: Domingues L, editors. PCR: Methods and Protocols. New York: Springer. pp. 87–100.
    https://doi.org/10.1007/978-1-4939-7060-5_5
    1. Taylor PM
    2. Thomas P
    3. Gorrie GH
    4. Connolly CN
    5. Smart TG
    6. Moss SJ
    (1999)
    Identification of amino acid residues within GABA(A) receptor beta subunits that mediate both homomeric and heteromeric receptor expression
    The Journal of Neuroscience 19:6360–6371.
  4. Software
    1. van der Loo M
    (2010)
    extremevalues, an R package for outlier detection in univariate data
    Extremevalues.

Article and author information

Author details

  1. Christian JG Tessier

    1. Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Canada
    2. Centre for Chemical and Synthetic Biology, University of Ottawa, Ottawa, Canada
    Contribution
    CJGT acquired and analyzed all electrophysiological data, except the data for the electrical fingerprinting, which was recorded and analyzed by RMS, CJGT, and CJBdC. CJGT and CJBdC interpreted the data and wrote the manuscript., Conceptualization, Formal analysis, Investigation, Writing – original draft, Writing – review and editing, Methodology
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-6006-7755
  2. Raymond M Sturgeon

    1. Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Canada
    2. Centre for Chemical and Synthetic Biology, University of Ottawa, Ottawa, Canada
    Contribution
    Conceptualization, Formal analysis, Investigation, Methodology, RMS recorded and analyzed electrophysiological data for the electrical fingerprinting experiments., Writing – original draft, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-0043-0241
  3. Johnathon R Emlaw

    1. Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Canada
    2. Centre for Chemical and Synthetic Biology, University of Ottawa, Ottawa, Canada
    Contribution
    Conceptualization, JRE acquired some of the initial recordings containing βAnc homopentamers., Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-7255-182X
  4. Gregory D McCluskey

    1. Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Canada
    2. Centre for Chemical and Synthetic Biology, University of Ottawa, Ottawa, Canada
    Contribution
    Conceptualization, Formal analysis, GDM reconstructed βAncS and acquired βAncS electrophysiological data., Investigation, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-6594-4665
  5. F Javier Pérez-Areales

    1. Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Canada
    2. Centre for Chemical and Synthetic Biology, University of Ottawa, Ottawa, Canada
    Present address
    Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
    Contribution
    FJPA acquired some of the initial recordings containing βAnc homopentamers., Investigation, Methodology, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-9525-9346
  6. Corrie JB daCosta

    1. Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Canada
    2. Centre for Chemical and Synthetic Biology, University of Ottawa, Ottawa, Canada
    Contribution
    CJBdC interpreted the data and wrote the manuscript, and supervised the project., Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Supervision, Writing – original draft, Writing – review and editing
    For correspondence
    cdacosta@uottawa.ca
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-9546-5331

Funding

Natural Sciences and Engineering Research Council of Canada (Discovery Grant RGPIN-2016-04801)

  • Corrie JB daCosta

Canadian Institutes of Health Research (Project Grant 377068)

  • Corrie JB daCosta

Ontario Council on Graduate Studies, Council of Ontario Universities (Graduate Scholarship)

  • Christian JG Tessier

New Frontiers in Research Fund (NFRFE-2018-453 00064)

  • Corrie JB daCosta

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

CJGT was funded in part by an Ontario Graduate Scholarship, while JRE was funded in part by a Natural Sciences and Engineering Research Council of Canada (NSERC) CREATE Scholarship. CJBdC acknowledges grants from NSERC (RGPIN-2016–04801), the Canada Foundation for Innovation (34475), the New Frontiers in Research Fund (NFRFE-2018–453 00064), and the Canadian Institutes of Health Research (377068).

Copyright

© 2022, Tessier et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 966
    views
  • 268
    downloads
  • 7
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Christian JG Tessier
  2. Raymond M Sturgeon
  3. Johnathon R Emlaw
  4. Gregory D McCluskey
  5. F Javier Pérez-Areales
  6. Corrie JB daCosta
(2022)
Ancestral acetylcholine receptor β-subunit forms homopentamers that prime before opening spontaneously
eLife 11:e76504.
https://doi.org/10.7554/eLife.76504

Share this article

https://doi.org/10.7554/eLife.76504

Further reading

    1. Structural Biology and Molecular Biophysics
    Maxim Oleynikov, Samie R Jaffrey
    Research Article

    The functional effects of an RNA can arise from complex three-dimensional folds known as tertiary structures. However, predicting the tertiary structure of an RNA and whether an RNA adopts distinct tertiary conformations remains challenging. To address this, we developed BASH MaP, a single-molecule dimethyl sulfate (DMS) footprinting method and DAGGER, a computational pipeline, to identify alternative tertiary structures adopted by different molecules of RNA. BASH MaP utilizes potassium borohydride to reveal the chemical accessibility of the N7 position of guanosine, a key mediator of tertiary structures. We used BASH MaP to identify diverse conformational states and dynamics of RNA G-quadruplexes, an important RNA tertiary motif, in vitro and in cells. BASH MaP and DAGGER analysis of the fluorogenic aptamer Spinach reveals that it adopts alternative tertiary conformations which determine its fluorescence states. BASH MaP thus provides an approach for structural analysis of RNA by revealing previously undetectable tertiary structures.

    1. Biochemistry and Chemical Biology
    2. Structural Biology and Molecular Biophysics
    Lingzhi Gao, Dian Chen, Yu Liu
    Research Article

    Riboswitches represent a class of non-coding RNA that possess the unique ability to specifically bind ligands and, in response, regulate gene expression. A recent report unveiled a type of riboswitch, known as the guanidine-IV riboswitch, which responds to guanidine levels to regulate downstream genetic transcription. However, the precise molecular mechanism through which the riboswitch senses its target ligand and undergoes conformational changes remain elusive. This gap in understanding has impeded the potential applications of this riboswitch. To bridge this knowledge gap, our study investigated the conformational dynamics of the guanidine-IV riboswitch RNA upon ligand binding. We employed single-molecule fluorescence resonance energy transfer (smFRET) to dissect the behaviors of the aptamer, terminator, and full-length riboswitch. Our findings indicated that the aptamer portion exhibited higher sensitivity to guanidine compared to the terminator and full-length constructs. Additionally, we utilized Position-specific Labelling of RNA (PLOR) combined with smFRET to observe, at the single-nucleotide and single-molecule level, the structural transitions experienced by the guanidine-IV riboswitch during transcription. Notably, we discovered that the influence of guanidine on the riboswitch RNA’s conformations was significantly reduced after the transcription of 88 nucleotides. Furthermore, we proposed a folding model for the guanidine-IV riboswitch in the absence and presence of guanidine, thereby providing insights into its ligand-response mechanism.