Ancestral acetylcholine receptor β-subunit forms homopentamers that prime before opening spontaneously
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.sa0Introduction
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).
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.
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.
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).
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.
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.
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.
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.
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
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 protocolcDNAs 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 protocolCombinations 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 protocolApproximately 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 protocolSingle-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 protocolSingle-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 protocolThe 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 protocolAppendix 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).
Data availability
All data generated and analysed during this study are included in the manuscript and as source data.
References
-
Functional relationships between agonist binding sites and coupling regions of homomeric Cys-loop receptorsThe Journal of Neuroscience 31:3662–3669.https://doi.org/10.1523/JNEUROSCI.5940-10.2011
-
On the stochastic properties of bursts of single ion channel openings and of clusters of burstsPhilosophical Transactions of the Royal Society of London. Series B, Biological Sciences 300:1–59.https://doi.org/10.1098/rstb.1982.0156
-
BookFitting and Statistical Analysis of Single-Channel RecordsIn: Colquhoun D, editors. Single-Channel Recording. Springer. pp. 483–587.https://doi.org/10.1007/978-1-4419-1229-9
-
Stoichiometry for α-bungarotoxin block of α7 acetylcholine receptorsNature Communications 6:8057.https://doi.org/10.1038/ncomms9057
-
The evolution of pentameric ligand-gated ion channelsAdvances in Experimental Medicine and Biology 683:11–23.https://doi.org/10.1007/978-1-4419-6445-8_2
-
Contingency between Historical Substitutions in the Acetylcholine Receptor PoreACS Chemical Neuroscience 11:2861–2868.https://doi.org/10.1021/acschemneuro.0c00410
-
Kinetic, Mechanistic, and Structural Aspects of Unliganded Gating of Acetylcholine Receptor ChannelsJournal of General Physiology 115:621–635.https://doi.org/10.1085/jgp.115.5.621
-
Cell surface expression of homomeric GABAA receptors depends on single residues in subunit transmembrane domainsThe Journal of Biological Chemistry 293:13427–13439.https://doi.org/10.1074/jbc.RA118.002792
-
Agonist and blocking actions of choline and tetramethylammonium on human muscle acetylcholine receptorsThe Journal of Physiology 587:5045–5072.https://doi.org/10.1113/jphysiol.2009.176305
-
Molecular evolution of the nicotinic acetylcholine receptor: an example of multigene family in excitable cellsJournal of Molecular Evolution 40:155–172.https://doi.org/10.1007/BF00167110
-
NACHO and 14-3-3 promote expression of distinct subunit stoichiometries of the α4β2 acetylcholine receptorCellular and Molecular Life Sciences 78:1565–1575.https://doi.org/10.1007/s00018-020-03592-x
-
Mapping a molecular link between allosteric inhibition and activation of the glycine receptorNature Structural & Molecular Biology 15:1084–1093.https://doi.org/10.1038/nsmb.1492
-
Improved resolution of single channel dwell times reveals mechanisms of binding, priming, and gating in muscle AChRThe Journal of General Physiology 148:43–63.https://doi.org/10.1085/jgp.201611584
-
The intrinsic energy of the gating isomerization of a neuromuscular acetylcholine receptor channelThe Journal of General Physiology 139:349–358.https://doi.org/10.1085/jgp.201110752
-
Local anaesthetics transiently block currents through single acetylcholine-receptor channelsThe Journal of Physiology 277:153–176.https://doi.org/10.1113/jphysiol.1978.sp012267
-
Alternate stoichiometries of alpha4beta2 nicotinic acetylcholine receptorsMolecular Pharmacology 63:332–341.https://doi.org/10.1124/mol.63.2.332
-
Ion channel block by acetylcholine, carbachol and suberyldicholine at the frog neuromuscular junctionProceedings of the Royal Society of London, Series B, Biological Sciences 225:329–355.https://doi.org/10.1098/rspb.1985.0065
-
Evolutionary history of the ligand-gated ion-channel superfamily of receptorsTrends in Neurosciences 18:121–127.https://doi.org/10.1016/0166-2236(95)93887-4
-
Delimiting the binding site for quaternary ammonium lidocaine derivatives in the acetylcholine receptor channelThe Journal of General Physiology 112:611–621.https://doi.org/10.1085/jgp.112.5.611
-
Evolution of vertebrate nicotinic acetylcholine receptorsBMC Evolutionary Biology 19:38.https://doi.org/10.1186/s12862-018-1341-8
-
The 5-hydroxytryptamine type 3 (5-HT3) receptor reveals a novel determinant of single-channel conductanceBiochemical Society Transactions 32:547–552.https://doi.org/10.1042/BST0320547
-
Structural mechanism of muscle nicotinic receptor desensitization and block by curareNature Structural & Molecular Biology 29:386–394.https://doi.org/10.1038/s41594-022-00737-3
-
BookInverse PCR for point mutation introductionIn: Domingues L, editors. PCR: Methods and Protocols. New York: Springer. pp. 87–100.https://doi.org/10.1007/978-1-4939-7060-5_5
-
Agonists block currents through acetylcholine receptor channelsBiophysical Journal 46:277–283.https://doi.org/10.1016/S0006-3495(84)84022-9
-
Identification of amino acid residues within GABA(A) receptor beta subunits that mediate both homomeric and heteromeric receptor expressionThe Journal of Neuroscience 19:6360–6371.
-
Back to the future: Rational maps for exploring acetylcholine receptor space and timeBiochimica et Biophysica Acta. Proteins and Proteomics 1865:1522–1528.https://doi.org/10.1016/j.bbapap.2017.08.006
-
Resurrecting ancient genes: experimental analysis of extinct moleculesNature Reviews. Genetics 5:366–375.https://doi.org/10.1038/nrg1324
-
Evolution of nicotinic acetylcholine receptor subunitsMolecular Biology and Evolution 15:518–527.https://doi.org/10.1093/oxfordjournals.molbev.a025951
-
Refined structure of the nicotinic acetylcholine receptor at 4A resolutionJournal of Molecular Biology 346:967–989.https://doi.org/10.1016/j.jmb.2004.12.031
-
Softwareextremevalues, an R package for outlier detection in univariate dataExtremevalues.
Article and author information
Author details
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
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)
Further reading
-
- Structural Biology and Molecular Biophysics
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.
-
- Biochemistry and Chemical Biology
- Structural Biology and Molecular Biophysics
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.