The superfamily of pentameric ligand-gated ion channels (pLGICs) is unique among ionotropic receptors in that the same overall structure has evolved to generate multiple members with different combinations of agonist specificities and permeant-ion charge selectivities. However, aside from these differences, pLGICs have been typically regarded as having several invariant functional properties. These include pore blockade by extracellular quaternary-ammonium cations in the micromolar-to-millimolar concentration range (in the case of the cation-selective members), and a gain-of-function phenotype, which manifests as a slower deactivation time course, as a result of mutations that reduce the hydrophobicity of the transmembrane pore lining. Here, we tested this notion on three distantly related cation-selective members of the pLGIC superfamily: the mouse muscle nicotinic acetylcholine receptor (nAChR), and the bacterial GLIC and ELIC channels. Remarkably, we found that, whereas low millimolar concentrations of TMA+ and TEA+ block the nAChR and GLIC, neither of these two quaternary-ammonium cations blocks ELIC at such concentrations; instead, both carry measurable inward currents when present as the only cations on the extracellular side. Also, we found that, whereas lidocaine binding speeds up the current-decay time courses of the nAChR and GLIC in the presence of saturating concentrations of agonists, the binding of lidocaine to ELIC slows this time course down. Furthermore, whereas mutations that reduce the hydrophobicity of the side chains at position 9′ of the M2 α-helices greatly slowed the deactivation time course of the nAChR and GLIC, these mutations had little effect—or even sped up deactivation—when engineered in ELIC. Our data indicate that caution should be exercised when generalizing results obtained with ELIC to the rest of the pLGICs, but more intriguingly, they hint at the possibility that ELIC is a representative of a novel branch of the superfamily with markedly divergent pore properties despite a well-conserved three-dimensional architecture.
Structure–function studies of bacterial pentameric ligand-gated ion channels (pLGICs) have led to the notion that both structure and function are remarkably well-conserved in the superfamily despite little amino acid sequence conservation. Indeed, even in the absence of the eponymous cysteine loop or the long intracellular linker between the M3 and M4 transmembrane α-helices, the bacterial members studied thus far open and desensitize upon binding extracellular ligands, much like their homologues from the nervous system of animals (Bocquet et al., 2007; Gonzalez-Gutierrez and Grosman, 2010; Parikh et al., 2011; Zimmermann and Dutzler, 2011; Gonzalez-Gutierrez et al., 2012).
From studies performed largely on the animal members of the superfamily, it has been concluded that—aside from being activated by different ligands and having opposite charge selectivities—pLGICs form a group of ion channels that share several functional properties. For example: (a) the transmembrane pore of the cation-selective members is blocked by extracellular quaternary-ammonium cations in the micromolar-to-millimolar concentration range (Neher and Steinbach, 1978; Adler et al., 1979; Sine and Steinbach, 1984; Ogden and Colquhoun, 1985; Marshall et al., 1990; Cuevas and Adams, 1994; Zhang et al., 1995; Blanchet and Dulon, 2001; Akk and Steinbach, 2003; Purohit and Grosman, 2006), discriminates poorly among monovalent cations (Adams et al., 1980; Yang, 1990), and displays some permeability to Ca2+ (Adams et al., 1980; Decker and Dani, 1990; Yang 1990; Zhou and Neher, 1993; Elenes et al., 2009); (b) the transmembrane pore of the anion-selective members is blocked by picrotoxin (Chang and Weiss, 1998, 1999; Etter et al., 1999; Sedelnikova et al., 2006; Bali and Akabas, 2007; Wang et al., 2007); (c) the rearrangement of the loop between the extracellular domain β-strands 9 and 10 (the C-loop) is an integral part of the conformational changes that occur upon ligand binding or gating (Chen et al., 1995; Hansen et al., 2005; Mukhtasimova et al., 2009; Hibbs and Gouaux, 2011; Purohit and Auerbach, 2013; Yoluk et al., 2013); and (d) mutations that reduce the hydrophobicity of the transmembrane pore lining invariably lead to a gain-of-function phenotype (regardless of the channel’s charge selectivity) that results from a stabilized open-channel conformation and manifests, for example, as an increased sensitivity to agonists, an increased unliganded-gating activity, longer bursts and clusters of single-channel openings, a slower time course of deactivation, and a slower time course of desensitization (Revah et al., 1991; Filatov and White, 1995; Labarca et al., 1995; Kearney et al., 1996; Chang and Weiss, 1998, 1999; Thompson et al., 1999; Kosolapov et al., 2000; Bianchi and Macdonald, 2001; Cymes et al., 2002; Burzomato et al., 2003; Grosman, 2003; Shan et al., 2003; Cymes et al., 2005; Papke and Grosman, 2014). Importantly, the discovery of the more distantly related bacterial and archaeal pLGICs (Tasneem et al., 2005) has afforded us the possibility to challenge these concepts in the framework of an even more diverse group of homologues. Although it seems likely that all members of the superfamily—from bacterial and archaeal to human—form ion channels gated by extracellular ligands, the degree to which more detailed aspects of molecular function are conserved remains to be ascertained.
Not much is known about the relationship between structure and function in bacterial or archaeal pLGICs, but some intriguing differences have already begun to emerge. A case in point is the C-loop of the extracellular domain. Whereas mutations to this loop have been found to impair the activation of animal pLGICs profoundly (Chen et al., 1995; Shen et al., 2012), we have recently shown that the entire C-loop of the bacterial homologue ligand-gated ion channel from Gloeobacter violaceus (GLIC) can be deleted without compromising its function as a proton-gated channel (Gonzalez-Gutierrez et al., 2013). Another example is the atypical phenotype of leucine-to-alanine mutations at the pore-lining position 9′ of the M2 transmembrane α-helix of ligand-gated ion channel from Erwinia chrysanthemi (ELIC; Gonzalez-Gutierrez et al., 2012), a finding that we investigate further here with additional mutations. Yet another example is the larger number of exponential components observed for the distribution of open times of ELIC recorded at saturating concentrations of agonist (propylamine) compared with the single component observed under comparable experimental conditions in the cases of the (animal) muscle nicotinic acetylcholine receptor (nAChR) and the α1 glycine receptor (Marabelli et al., 2015).
In this paper, we set out to study some of the properties of the transmembrane pore. We compared the effects of pore blockers (tetramethylammonium [TMA+], tetraethylammonium [TEA+], and lidocaine) and mutations to the pore’s hydrophobic lining on the function of three distantly related cation-selective members of the superfamily: the mouse muscle nAChR, and the bacterial GLIC and ELIC channels. Strikingly, we found that, whereas millimolar concentrations of extracellular TMA+ and TEA+ block the nAChR and GLIC, neither of these two quaternary-ammonium cations blocks ELIC at a concentration as high as 50 mM, and instead, both carry measurable inward currents when present as the only cations on the extracellular side. Also, we found that, whereas lidocaine binding speeds up the time course of the current decay of the nAChR and GLIC on sustained application of saturating concentrations of agonist, the binding of lidocaine to ELIC slows this time course down, and the crystal structure of the ELIC–lidocaine complex shows the anesthetic bound to the cavity lined by the back of the M2 and the front of the M1 and M3 transmembrane α-helices. In marked contrast, in both the nAChR (Leonard et al., 1988; Charnet et al., 1990) and GLIC (Hilf et al., 2010), lidocaine (or its trimethyl-ammonium counterpart, QX-222) has been found to bind to the transmembrane pore, i.e., the cavity lined by the front side of the five M2 α-helices. Finally, whereas mutations that reduce the hydrophobicity of the pore-lining side chains greatly slowed down the time course of deactivation of the nAChR and GLIC—a phenomenon often regarded as a hallmark of the entire superfamily—we found that these mutations have little or even the opposite effect when engineered in ELIC.
MATERIALS AND METHODS
cDNA clones, mutagenesis, and heterologous expression
HEK-293 cells were transiently transfected with wild-type or mutant complementary DNAs (cDNAs) coding the mouse muscle adult-type nAChR (α1, β1, δ, and ε subunits), GLIC or ELIC using the calcium-phosphate precipitation method. nAChR cDNAs in the expression vector pRBG4 were provided by S. Sine (Mayo Clinic, Rochester, MN). GLIC and ELIC cDNAs in the pcDNA 3.1 (–) vector (Life Technologies) were prepared as described previously (Gonzalez-Gutierrez and Grosman, 2010; Gonzalez-Gutierrez et al., 2012). For the calculation of identities and similarities between pairs of amino acid sequences, the number of identical or similar residues was divided by the average number of total residues in the two aligned sequences. Point mutations were engineered using the QuikChange site-directed mutagenesis kit (Agilent Technologies) and were confirmed by dideoxy sequencing.
Electrophysiological recordings and analysis
Ensemble (macroscopic) currents were recorded from transfected HEK-293 cells using the outside-out configuration of the patch-clamp technique, whereas single-channel currents were recorded using either the outside-out or the cell-attached configuration, at ∼22°C. For outside-out recordings, the ligand was applied to the external aspect of excised patches as rapid jumps (solution-exchange time10–90% < 150 µs). These step-changes in the concentration of ligand were achieved by the rapid switching of two solutions flowing from either barrel of a piece of theta-type capillary glass mounted on a piezo-electric device (Burleigh-LSS-3100; discontinued). For all outside-out recordings, the patch-pipette solution consisted of (in mM) 110 KF, 40 KCl, 1 CaCl2, 11 EGTA, and 10 HEPES/KOH, pH 7.4; the solutions flowing through the two barrels of the perfusion tubing differed among experiments and are indicated in the text. For cell-attached recordings, the pipette solution consisted of (in mM) 142 KCl, 5.4 NaCl, 1.8 CaCl2, 1.7 MgCl2, 10 HEPES/KOH, pH 7.4, and the indicated agonist. The effective bandwidth for data analysis was DC–5 kHz for both macroscopic and single-channel currents. For display purposes, the macroscopic current traces were decimated, and the single-channel traces were filtered such that the cut-off frequency was 2–5 kHz. Both macroscopic and single-channel current recordings were analyzed using pClamp 9.0 (MDS Analytical Technologies).
To determine the effects of quaternary-ammonium cations and lidocaine on the values of macroscopic peak currents and decay time constants of the muscle nAChR and GLIC, we used paired experiments. In these, the response of each patch of membrane to a pulse of agonist was recorded sequentially in the absence, in the presence, and back again in the absence (that is, upon washout) of the compound under study with the goal of reducing the impact of the variability from patch to patch. Thus, each patch acted as its own control and generated a value for the effect of each added compound. These values (obtained from several patches) were averaged, and the corresponding mean and standard error are presented in the form of bar graphs throughout the paper as normalized peak currents and normalized time constants. In the case of ELIC, however, paired experiments proved difficult because the channel activity usually ran down as a function of time after patch excision. Therefore, to maximize the amount of data obtained from each patch, the number of maneuvers within each single experiment was kept to a minimum. Thus, responses in the presence and in the absence of the tested compounds were recorded from separate (i.e., unpaired) patches of membrane. The ratios between the means of these values were calculated, and their corresponding standard errors were propagated; it is these values of ELIC that are presented in the form of bar graphs throughout the paper.
Because of the slow desensitization time courses of wild-type GLIC and ELIC (with time constants, in some patches, as slow as 8 s), the agonist typically had to be applied for long intervals (∼1 min) for the currents to reach steady state, and hence, for accurate estimates of the exponential decay time constants to be obtained. These intervals proved to be too long, however, considering the average lifetime of a fast-perfused outside-out patch, on the one hand, and the multiple solution changes and consecutive agonist applications that we intended to perform on any given patch, on the other. Therefore, as a compromise, we applied 1-min pulses of agonist to only a few patches per tested condition and estimated the steady-state value of the current relative to the peak. With this information, we then proceeded to fix the steady-state current of fits to responses to shorter applications of agonist (5 s). In all cases, the steady-state values of the macroscopic currents on sustained application of agonist turned out to be nearly zero.
In an attempt to assess the statistical significance of the differences between the mean values of current decay time constants obtained under different conditions, we performed Student’s t tests (paired or unpaired, depending on the experimental design) and report the corresponding two-tailed p-values.
Protein purification, crystallization, and x-ray diffraction
ELIC was purified as described in previous work (Gonzalez-Gutierrez et al., 2012). Concentrated protein (∼10 mg/ml) was mixed with cysteamine (final concentration = 10 mM), Escherichia coli total polar lipids (Avanti Polar; final concentration = 0.5 mg/ml), and lidocaine (final concentration = 15 mM) or the bromo derivative N-1-(4-bromophenyl)-N-2,N-2 diethyl glycinamide (Chembridge; final concentration = 15 mM). The protein was crystallized by vapor diffusion in hanging drops at 4°C. The reservoir solution consisted of 200 mM (NH4)2SO4, 10–12% (wt/vol) PEG 4000 and 50 mM N-(2-acetamido) iminodiacetic acid/NaOH, pH 6.5–6.7. For cryoprotection, crystals were transferred for a few seconds into reservoir solution supplemented with 30% (vol/vol) ethylene-glycol and were subsequently flash frozen in liquid nitrogen. Diffraction data were collected at 100 K at the 21 ID–D/F/G beamline at the Advanced Photon Source (Argonne National Laboratory, Argonne, IL). The data were indexed, integrated, and scaled using HKL2000 (Otwinowski and Minor, 1997). The structures were solved by molecular replacement using PDB deposition 2VL0 as the search model. The Auto Build function in PHENIX (Adams et al., 2010) was used to generate a first model that was subsequently refined in PHENIX and Coot (Emsley et al., 2010). No attempts were made to refine regions of poor electron density as defined in the starting model. Cross-validation used 5% of the data in the calculation of the Rfree value. PyMOL was used to prepare molecular images (The PyMOL Molecular Graphics System, Version 1.3; Schrödinger). The structure of ELIC was solved in the presence of cysteamine and lidocaine to a resolution of 3.65 Å; the wavelength used was 0.97872 Å. Because only a weak density for lidocaine was observed, a second structure in the presence of cysteamine and the bromo derivative N-1-(4-bromophenyl)-N-2,N-2 diethyl glycinamide was also solved (to a resolution of 4.23 Å). To confirm the presence of the heavy atom, a third dataset was obtained using a wavelength of 0.91838 Å (solved to a resolution of 4.82 Å).
Effects of quaternary-ammonium cations and lidocaine on the nAChR and GLIC
To gain insight into the extent to which concepts derived from the study of the bacterial homologues of the pLGIC superfamily can be applied to their animal counterparts, we decided to compare some of their functional properties. We started by testing the effects of TEA+ and lidocaine on the (heteromeric) adult-type nAChR from mouse muscle, and the (homomeric) GLIC and ELIC channels from bacteria. Averaging the four types of subunit, the amino acid sequence of this nAChR is 16.7% identical (25.4% similar) to that of GLIC and 14.7% identical (23.1% similar) to that of ELIC. The amino acid sequence of GLIC is 20.4% identical (35.2% similar) to that of ELIC. Disregarding the highly variable M3–M4 linker, these values become 19.0% (identity) and 33.0% (similarity) for the muscle AChR and GLIC; 16.9% (identity) and 30.4% (similarity) for the muscle AChR and ELIC; and 20.5% (identity) and 35.8% (similarity) for GLIC and ELIC.
TEA+ applied to the extracellular side of outside-out patches at a concentration of 0.5–5.0 mM reduced the peak amplitude of the macroscopic currents flowing through the nAChR and GLIC with little (if any) effect on their time constants of desensitization (Figs. 1 and 2, and Table 1). The lower peak values of the nAChR-mediated macroscopic currents can be attributed to a combination of a decreased single-channel conductance and a lower open probability (Fig. 1 B), consistent with TEA+ acting as a pore blocker that enters and exits the pore fast enough to reduce the unitary current amplitude yet slow enough for the individual blocking events to be discernible as sojourns in the zero-current level. Although analogous recordings from GLIC were not clear enough to draw meaningful conclusions (owing to the low single-channel conductance of this bacterial channel even in the absence of TEA+), we surmise that the effect of this organic cation on GLIC is essentially the same. Regarding desensitization of these two channels, we would like to highlight the observation that the presence of blocking concentrations of TEA+ did not affect their time courses. Thus, the presence of a blocker inside the pore of a pLGIC’s crystal structure, for example, should not be taken to imply that the crystallized conformation represents the open state because pore block (at least by TEA+) does not prevent channel desensitization.
Extracellular lidocaine also reduced the peak currents mediated by the nAChR (Fig. 3, A and B), but unlike TEA+, this local anesthetic sped up the time course of current decay by a factor of ∼100 (Fig. 3 C and Table 1). The lower peak values of the nAChR-mediated macroscopic currents can be attributed to the dramatically reduced single-channel open probability elicited by 100-µM ACh in the presence of 0.5 mM lidocaine (Fig. 3 A, insets). Because lidocaine is a slower blocker of the nAChR than is TEA+ (note that, in the presence of lidocaine, the amplitude of the single-channel currents remained essentially unaltered), it may be asked whether the faster decay of the macroscopic currents reflects faster desensitization or merely the kinetics of lidocaine entry into the open-channel pore. To learn more about the effect of lidocaine, we coapplied it with ACh for short intervals (Fig. 3 D). If lidocaine acted as a fast-dissociating open-channel blocker that does not alter the kinetics of closed–open–desensitized state interconversions, for example, the current upon washout of both lidocaine and ACh would be expected to rebound and attain the same value as that observed at the end of ACh-only pulses. This is not what we observed, however. Upon lidocaine and ACh washout, the signal remained close to the zero-current baseline (Fig. 3 D, vertical dashed line), a finding that is entirely consistent with lidocaine strongly biasing the channel’s conformational landscape toward a nonconductive conformation (such as the desensitized state). However, it may be argued that this finding is also consistent with lidocaine dissociating slowly from the open-channel pore. As a result, we do not refer to this particular current decay time course as the time course of desensitization. Clearly, more experiments would be needed to elucidate the relative contributions of pore block and entry into desensitization to the current decay time course of the nAChR in the presence of both ACh and lidocaine. Obtaining such mechanistic insight was beyond the scope of this paper, however.
In the case of GLIC, the effect of extracellular 0.5 mM lidocaine on the macroscopic currents was qualitatively the same as that on the nAChR. Indeed, the peak currents were decreased to a large extent (Fig. 4, A and B), and the current-decay time course on continuous exposure to pH 4.5 extracellular solution was sped up, although in this case, only by a factor of ∼3 (Fig. 4 C and Table 1). In fact, the ranges of time constant values of GLIC with and without lidocaine overlapped partially (Table I), but the application of lidocaine to each individual patch (as illustrated in Fig. 4 A) sped up the time course of current decay in all 12 tested patches. The time constant fitted to each patch varied widely from patch to patch, but its value always became faster upon adding lidocaine (and thus, we cannot explain the slowing effect reported by Velisetty and Chakrapani (2012) for lidocaine on the current-decay time course of GLIC). Again, although single-channel recordings from GLIC in the presence of lidocaine were not clear enough to draw conclusions, it seems likely that this local anesthetic also lowers the channel’s open probability, much like it does on the nAChR.
Effects of quaternary-ammonium cations and lidocaine on ELIC
When tested on ELIC, extracellular TEA+ did not block the peak-current values elicited by saturating (10 mM) cysteamine, even when applied at a concentration as high as 50 mM (Fig. 5, A and B); neither did TEA+ block the channel when applied to the intracellular side of outside-out patches at a concentration of 5 mM (Fig. 5, A and B). Not surprisingly then, we found that ELIC’s kinetics of desensitization are essentially unaffected by the presence of TEA+ on either side of the membrane (Fig. 5 C and Table 1). Importantly, ELIC’s open probability between desensitized intervals in the presence of 10 mM cysteamine and the absence of TEA+ is similar (∼0.95; Gonzalez-Gutierrez et al., 2012) to that of the muscle nAChR in the presence of 100-µM ACh and the absence of TEA+ (∼0.85; Gonzalez-Gutierrez and Grosman, 2010), and therefore, the possibility that a lower open probability in the absence of blocker underlies the much lower extent of block observed for ELIC can be completely ruled out.
Application of extracellular 5 mM lidocaine to ELIC, on the other hand, reduced the peak currents elicited by pulses of 10 mM cysteamine by a factor of ∼16 (mean peak current in the absence of lidocaine = 466 ± 109 pA, n = 16 patches; mean peak current in the presence of lidocaine = 29 ± 8 pA, n = 14 patches), and single-channel recordings also revealed a markedly reduced open probability (Fig. 6, A and B). Note, however, that a concentration as high as 5 mM was needed for lidocaine to block ELIC to an extent that is similar to that achieved by 0.5 mM on the nAChR or GLIC (as judged from the reduction in the peak-current values). Regarding the decay of the macroscopic current in the presence of saturating agonist, application of extracellular lidocaine to ELIC had the opposite effect of that on the nAChR or GLIC: the time course was slowed down to the extent that no decline of the current was observed during 5-s pulses of 10 mM cysteamine in the continuous presence of 5 mM extracellular lidocaine (Fig. 6 A, inset, and Table 1). Puzzled by the strikingly different effects of lidocaine on the time courses of the current decay of the nAChR and GLIC, on the one hand, and that of ELIC, on the other, we set out to solve the crystal structure of ELIC in the presence of agonist (cysteamine) and either lidocaine or a brominated analogue. Although the structures were solved to a low resolution (≥3.65 Å; Table 2), the bromo derivative was clearly identified in the cavity lined by the back side of M2 and the front sides of M1 and M3 of several subunits (Fig. 7; the overall structure of ELIC, however, was not affected). Moreover, anomalous electron densities were found in these same regions (Fig. 7 D). Evidently, it would be important to determine whether and where lidocaine binds to the open-channel conformation of ELIC so as to be able to make a more direct comparison with the existing data on the binding of lidocaine (or its quaternary-ammonium derivatives) to the nAChR and GLIC (Leonard et al., 1988; Charnet et al., 1990; Hilf et al., 2010). It may well be that a given compound binds to different cavities of the channel in its different conformational states. Also, it is important to note that we have used lidocaine here only as a probe of structure–function relationships in a bacterial channel; our results should not be regarded as shedding light on the binding of anesthetics to their animal targets.
To probe the effect of structural perturbations to the lidocaine-binding site of ELIC, we engineered five single-point mutants (Y245A, Y248A, and T249A, in M2—corresponding to positions 14′, 17′, and 18′—and D263A and I266A, in M3; amino-acid numbering system as in Gonzalez-Gutierrez et al. ) and one double mutant (Y203A + Y204A, in M1), and recorded their responses to 5-s pulses of a saturating concentration of agonist (Fig. 8). Three of these (Y203A + Y204A, T18′A, and D263A) failed to give rise to functional expression, I266A was wild-type-like (Table 3), and Y14′A and Y17′A displayed faster desensitization time courses (by factors of ∼17 and ∼12, respectively; Table 3) and lower open probabilities. The combination of these two single-point mutations into a double mutant gave rise to, essentially, the same phenotype (Fig. 8 and Table 3). We also tested the effect of 5 mM lidocaine on two of these mutants. When applied to the wild type-like I266A ELIC, lidocaine reduced the peak-current amplitude and slowed down the time course of current decay much like it did on the wild-type channel (Fig. 8 I). When applied to the loss-of-function Y14′A mutant, lidocaine decreased the peak-current amplitude to the extent that its effect on the current-decay time course could not be determined (Fig. 8 J). In summary, no mutation to the lidocaine-binding sites mimicked or prevented all the effects of lidocaine binding on channel function. However, the marked effects of the Y14′A and Y17′A mutations on the kinetics of desensitization and on the open probability in the presence of a saturating concentration of agonist does indicate that perturbations to the M1–M2–M3 cavity—caused by mutations, and perhaps, the binding of small molecules—can tilt the conformational free-energy landscape of ELIC.
Because pore block by extracellular quaternary-ammonium cations is one of the most characteristic properties of the cation-selective pLGICs, we were intrigued by our finding that 50 mM TEA+ did not block the currents mediated by ELIC (Fig. 5, A and B). We reasoned that this lack of block could be the result of the open-channel pore of ELIC being too narrow to allow the access of a cation as large as TEA+ (∼6.9 Å diam; Conway et al., 1966; Edward, 1970; Ue, 1994). Alternatively, at the other end of the spectrum of possibilities, the pore of ELIC in the open state could be wide enough and have the necessary electrostatic properties for TEA+ to permeate nearly as fast as the mixture of K+, Na+, Ca2+, and Mg2+ of our extracellular solution in such a way that the competition between TEA+ and these inorganic cations would not result in a decreased conductance.
To test the possibility that the extracellular entrance to the transmembrane pore of ELIC is simply too narrow, we mutated the leucines at position 9′ and the phenylalanines at position 16′ to alanines and leucines, respectively, in an attempt to widen the pore cavity. The side chains at these two positions project into the lumen of the channel pore in the crystallized nonconductive conformation of ELIC (Hilf and Dutzler, 2008) and are likely to also do so in the open-channel conformation; in the muscle nAChR, they certainly do (Cymes et al., 2005). As was the case for wild-type ELIC, however, neither the peak value of the macroscopic currents nor the time course of the current response recorded from the L9′A + F16′L double mutant was affected by the continuous presence of extracellular TEA+ at a concentration as high as 50 mM (Fig. 9, A and B; and Table 1). The same was the case for TMA+ (a smaller cation, ∼5.7 Å diam; Conway et al., 1966; Edward, 1970; Ue, 1994) on both the wild-type and the L9′A + F16′L wider-pore mutant at a concentration of 50 mM (Fig. 9, C–E; and Table 1). To emphasize the atypical properties of ELIC, Fig. 10 illustrates the blocking effect of TMA+ on the nAChR at the macroscopic and single-channel levels when also applied to the extracellular side at a concentration of 50 mM. TMA+ is both a fast pore blocker and an agonist of the nAChR.
Regarding the open probability between desensitized intervals of the L9′A + F16′L mutant of ELIC in the presence of 10 mM cysteamine and the absence of quaternary-ammonium cations, we note that its value (∼0.45; Gonzalez-Gutierrez et al., 2012) is lower than the muscle-nAChR’s (∼0.85 in the presence of 100-µM ACh and the absence of quaternary-ammonium cations; Gonzalez-Gutierrez and Grosman, 2010) only by a factor of ∼1.9. In other words, although lower, this open probability is not low enough to account for the negligible block of this double mutant by 50 mM TEA+ or TMA+.
Quaternary-ammonium cations carry measurable currents through ELIC
To test the possibility that TMA+ can conduct measurable currents through ELIC, we replaced the extracellular mixture of K+, Na+, Ca2+, and Mg2+ present in the agonist-containing barrel of the theta-type tubing of our perfusion system with a solution containing TMA+ and agonist alone. More specifically, this solution consisted of (in mM) 142 TMACl, 10 cysteamine chloride, and 10 HEPES/TMAOH, pH 7.4; the solution flowing through the other barrel of the theta tube (that is, the agonist-free solution) remained unchanged (in mM: 142 KCl, 5.4 NaCl, 1.8 CaCl2, 1.7 MgCl2, and 10 HEPES/KOH, pH 7.4). To allow for a comparison across channel types, we also applied this type of pulse protocols to the nAChR and GLIC. Because TMA+ is a high-efficacy agonist of the nAChR, no additional ligands were added to the TMA+-containing solution when assaying this channel. In the case of GLIC, the pH of the 142-mM TMA+ solution was buffered with 10 mM acetic acid (instead of HEPES), and its value was adjusted to 4.5 with TMAOH.
The responses of the nAChR and GLIC are illustrated in Fig. 11, whereas that of ELIC is illustrated in Fig. 12. In the case of the nAChR, the ratio between the macroscopic current values before and after TMA+ removal was ∼0.08, thus suggesting that TMA+ carries a current that is ∼8.0% of that carried by the mixture of K+, Na+, Ca2+, and Mg2+ present in the agonist-free solution (Fig. 11 A). Because the fast deactivation of the nAChR (time constant = 1.1 ± 0.1 ms; n = 12) could have led to an underestimation of the peak-value reached upon TMA+ removal, we repeated this type of experiment with a mutant nAChR that deactivates more slowly (εT12′P; Ohno et al., 1995); in this case, the ratio was estimated to be ∼0.01 (Fig. 11 B). In other words, a value of ∼1% seems to be a more accurate estimate of the ratio of the current carried by TMA+ through the nAChR relative to the current carried by the mixture of K+, Na+, Ca2+, and Mg2+. In the case of GLIC, the response was very similar to that of the nAChR: we found the corresponding ratio to be ∼0.02 (Fig. 11 C), although, owing to the fast deactivation time course of GLIC (time constant = 0.91 ± 0.20 ms, in the absence of TMA+; n = 11), this number is likely to be an upper bound. Fig. 11 D summarizes the numerical values obtained from these experiments.
In striking contrast, the ratio of intra–pulse-to-tail current values corresponding to ELIC was much larger than those obtained for the nAChR or GLIC. Indeed, this ratio was ∼0.25 for wild-type ELIC and ∼0.35 for the wider-pore mutant (Fig. 12, A–C). We also measured the amplitude of inward single-channel currents elicited by 10 mM cysteamine at –80 mV in outside-out patches bathed by either the K+-Na+-Ca2+-Mg2+ solution or the 142-mM TMA+ solution on the extracellular side (Fig. 12 D). Under these conditions, we found the ratio of unitary-current amplitudes of wild-type ELIC to be 0.22 ± 0.01 (n = 13), a value that is very similar to that estimated from macroscopic current recordings (Fig. 12, A–C). Similar attempts to record single-channel currents carried by TMA+ alone through the nAChR or GLIC proved unsuccessful, which is hardly surprising given the small amplitude of the macroscopic currents shown in Fig. 11.
We also tested whether TEA+ permeates through the pore of ELIC. Using the same type of pulse protocol as that illustrated for TMA+ in Fig. 12 (A and B), we found the ratio between the current value reached at the end of the pulse of TEA+ plus agonist and the peak value reached upon switching the external solution to the mixture of K+, Na+, Ca2+, and Mg2+ to be ∼0.09 (Fig. 12 E). That is, with a solution containing (in mM) 142 TEACl, 10 cysteamine chloride, and 10 HEPES/TEAOH, pH 7.4 on the extracellular side, the inward currents were as large as ∼9% of their value when the extracellular solution consisted of (in mM) 142 KCl, 5.4 NaCl, 1.8 CaCl2, 1.7 MgCl2, and 10 HEPES/KOH, pH 7.4. Nearly the same value (∼0.10) was obtained with the L9′A + F16′L double mutant of ELIC (Fig. 12, F and G), indicating that the larger side chains of leucine at position 9′ and phenylalanine at position 16′ do not affect the rate of TEA+ permeation through the channel. Collectively, the results with TMA+ and TEA+ strongly suggest that, in the open-channel conformation, the pore of ELIC is wider at its narrowest constriction than are the pores of the muscle nAChR or GLIC.
Although TMA+ and TEA+ permeate faster through ELIC than they do through the nAChR or GLIC, as previously mentioned, the currents these cations carry through ELIC when present as the only extracellular cations are nevertheless smaller than the current carried by the mixture of K+, Na+, Ca2+, and Mg2+ (by a factor of ∼3–4 in the case of TMA+ and ∼11 in the case of TEA+). However, the addition of TMA+ or TEA+ (up to a concentration of 50 mM) to the K+-Na+-Ca2+-Mg2+ solution had hardly any effect on the size of the peak-current amplitudes carried by this mixture of inorganic cations alone (Figs. 5 and 9). Thus, this lack of blockade by TMA+ or TEA+ cannot be ascribed to the channel being too narrow to enter the pore (because these cations certainly conduct currents) or to TMA+ and TEA+ permeating as fast as the mixture of inorganic cations (because the currents carried by TMA+ and TEA+ are smaller). Thus, it seems as though the only explanation for the lack of current block illustrated in Figs. 5 and 9 is that TMA+ and TEA+ interact with the open-channel pore of ELIC with such low affinity that the other cations also present in the extracellular solution—that is, K+ (142 mM), Na+ (5.4 mM), Ca2+ (1.8 mM), and Mg2+ (1.7 mM)—outcompete them, at least when TMA+ and TEA+ are present at a concentration of 50 mM. To test this idea, we repeated the competition experiments now using a higher concentration of TEA+ and a simpler background solution (Fig. 13). More specifically, the TEA+-containing solution consisted of (in mM): 142 KCl, 142 TEACl, 10 cysteamine chloride, and 10 HEPES/KOH, pH 7.4; the solution flowing through the other barrel of the theta tube (that is, the TEA+-free solution) consisted of (in mM): 142 KCl and 10 HEPES/KOH, pH 7.4. Under these conditions, we found that the current carried by the equimolar TEA+-K+ solution was ∼53% of that recorded upon TEA+ removal, i.e., in the presence of 142 mM K+ alone, in the case of the wild type, and ∼54% in the case of the wider-pore mutant (Fig. 13 C). Evidently, a higher concentration of TEA+ is required to block K+ currents through ELIC than is required in the case of the nAChR or GLIC. Finally, we reasoned that—when present at the same concentration—a quaternary-ammonium cation larger than TEA+ would block K+ currents to a larger extent. As shown in Fig. 14, this expectation was fully borne out by experiments with 142 mM tetrapropylammonium (TPA+) chloride (TPACl; Fig. 14 A); in the presence of an equimolar mixture of KCl and TPACl (both at a concentration of 142 mM) applied to the extracellular side of ELIC, the current was ∼30% as large as that carried by the K+ solution alone. Fig. 14 also shows that 142 mM TPA+ carries a small but measurable current through ELIC when present as the only cation on the extracellular side (amounting to ∼3.5% of the current carried by the mixture of K+, Na+, Ca2+, and Mg2+ alone; Fig. 14 B), which is intriguing in light of the large diameter of this cation (∼7.8 Å; Conway et al., 1966; Edward, 1970; Ue, 1994).
It may be argued that the cation-conduction properties of ELIC could have been characterized much more directly and in more quantitative detail at the single-channel level. However, here, we specifically wanted to contrast these properties with those of the nAChR and GLIC, and given the low unitary conductance of the latter (even in the absence of blockers), this would have been nearly impossible to achieve using single-channel recordings. Thus, we felt that it was desirable to consider an experimental design that could be applied to all three channels under the same experimental conditions. The importance of including a second bacterial channel in the comparison seemed to fully justify the more indirect nature of macroscopic recordings.
Leucine-to-alanine mutations at position 9′ speed up the deactivation of ELIC
To gain insight into the interconversion among the closed, open, and desensitized states of ELIC, we compared the effects of reducing the hydrophobicity of the pore lining on the kinetics of deactivation. Channel deactivation is the decay of the macroscopic current that occurs on agonist removal, and thus, its kinetics depend on the rates of agonist dissociation and closed–open–desensitized state interconversions, with the contribution of each of these rates varying among different channels (Jones and Westbrook, 1995; Elenes et al., 2006). On the other hand, for channels that activate much faster than they desensitize, the kinetics of desensitization measured at saturating concentrations of agonist—that is, the kinetics of the macroscopic current decay that occurs during the application of agonist—depend on the rates at which the channel enters and leaves the desensitized state(s) while fully bound to the ligand. Thus, the time courses of deactivation and desensitization reflect different (although somewhat overlapping) sets of rates.
For these experiments, we chose the side chain at position 9′ (a leucine in all subunits of the muscle nAChR and ELIC, and an isoleucine in GLIC) because the effect of mutations at this position on the kinetics of deactivation—or, its single-molecule equivalent, the mean duration of bursts of single-channel openings (Wyllie et al., 1998)—has been extensively documented for the pLGICs from vertebrates (for example: Filatov and White, 1995; Labarca et al., 1995; Chang and Weiss, 1998; Kosolapov et al., 2000; Bianchi and Macdonald, 2001; Papke and Grosman, 2014). Fig. 15 illustrates the effect of shortening these side chains by mutation to alanine in the muscle nAChR and GLIC. The presence of a mutant alanine at this position of the two α subunits of the nAChR and in all five subunits of GLIC slowed down deactivation (by factors of ∼90 and ∼178, respectively; Table 4), a phenomenon shared by all other pLGICs for which the effect of this type of mutation on deactivation has been recorded. However, in ELIC, the leucine-to-alanine mutation sped up deactivation, with the time constant of the monoexponential fit being shortened by a factor of ∼3.5 (Fig. 16 and Table 4; see also Gonzalez-Gutierrez et al., 2012). To rule out the possibility that the phenylalanine at the pore-lining position 16′ of ELIC has anything to do with this unexpected behavior (indeed, five phenylalanines at this position seem to occur only in some bacterial and invertebrate homologues whose functional properties have not been extensively characterized yet), we also recorded the deactivation time course of the wider-pore L9′A + F16′L mutant that we used in our previously mentioned permeation studies; a leucine at position 16′ is present in four of the five subunits of the muscle nAChR, whereas an isoleucine at this position is present in all five subunits of GLIC. Remarkably, the deactivation time course of this double mutant of ELIC was faster too, also by a factor of ∼3.5 (Fig. 16, C and F; and Table 4; Gonzalez-Gutierrez et al., 2012). Mutating the leucine at position 9′ of ELIC to glutamine or glutamate—two polar amino acids—in the wild-type background gave rise to, essentially, the same result: these mutations did not slow down the deactivation time course of ELIC (Fig. 16, D–F; and Table 4).
As for desensitization, the effect of mutating the leucine at position 9′ of ELIC to alanine (both in the wild-type and the F16′L backgrounds) was unremarkable: the decay of the macroscopic current became much slower, to the extent of being negligible during the 5-s applications of 10 mM cysteamine (Fig. 9 A), largely as observed for the muscle nAChR (Filatov and White, 1995; Papke and Grosman, 2014), GLIC (Bocquet et al., 2007; Gonzalez-Gutierrez et al., 2012), and all other pLGICs examined to date (e.g., Bianchi and Macdonald, 2001).
Collectively, the effects of mutations at position 9′ of the nAChR and GLIC on the time courses of deactivation and desensitization suggest that reducing the hydrophobicity of the transmembrane pore-lining destabilizes the closed and desensitized states of these two channels relative to the open state. (Although changes in the agonist-dissociation rate constant would have contributed to affect the kinetics of deactivation, we deem it unlikely for a mutation in the middle of the M2 α-helix to affect the unbinding kinetics from the distant transmitter-binding sites.) This proposed relative stabilization of the open state is consistent with the longer clusters of single-channel openings and the higher closed-to-open gating equilibrium constant of nAChR 9′ mutants inferred from unitary currents recorded at saturating concentrations of agonist (Cymes et al., 2002). In marked contrast, the same experimental maneuvers applied to ELIC destabilize the open and the desensitized states, with the destabilization of the latter being larger. That the open state of ELIC becomes destabilized relative to the closed state upon mutation was also evident at the single-channel level, in recordings performed at saturating concentrations of cysteamine that showed an intracluster open probability lower than the wild-type’s (Gonzalez-Gutierrez et al., 2012).
Admittedly, we have not tested the effect of all other natural side chains at position 9′ of ELIC, but the results obtained with alanine, glutamine, and glutamate seem compelling enough to conclude that, in this channel, the relationship between pore hydrophobicity, pore-water content, and the relative stabilities of the closed, open, and desensitized states is quite unique. To our knowledge, ELIC is the only pLGIC studied to date that deviates from the rule that mutations that reduce the hydrophobicity of the transmembrane pore lining stabilize the open-channel conformation relative to both the closed and desensitized states.
A more polar lining increases the affinity of ELIC’s pore for cysteamine
The responses of the glutamine and glutamate mutants of ELIC to cysteamine pulses (Fig. 16, D and E) differed from those recorded from the other tested constructs (Fig. 16, A–C) in that the currents transiently rose upon cysteamine removal. The solutions alternately bathing the outside-out patches in these experiments were the K+-Na+-Ca2+-Mg2+ solution with or without 10 mM cysteamine, and thus, the observed current overshoots indicate that this agonist blocked the currents carried by these two mutants to some extent (with pKa1 and pKa2 values of ∼8.3 and ∼10.8 [Avdeef and Brown, 1984], a 10-mM cysteamine solution at pH 7.4 is expected to contain ∼8.9 mM of the doubly protonated, positively charged form). However, a careful examination of the electrophysiological responses recorded from the wild type, the L9′A single mutant, and the L9′A + F16′L double mutant failed to reveal such fast block by cysteamine, which is important in the context of this paper because these were the variants of ELIC that we used for the permeation studies described above. Therefore, our decision to disregard the presence of 10 mM cysteamine in the experiments that probed the permeation properties of ELIC is fully justified. Clearly, it seems as though the presence of a more polar lining in the L9′Q and L9′E mutants of ELIC increases the affinity of the pore for the, also polar, cysteamine.
Before this work, the characterization of two bacterial pLGICs had suggested that the functional properties of the noneukaryotic members of the superfamily fell well within the range of variability revealed by decades of work on their animal—both vertebrate and invertebrate—counterparts. Indeed, the differences identified between the bacterial and animal pLGICs were not much larger than those found between distantly related pLGICs from animals. For example, although the notion that GLIC is gated by an unusual agonist may have seemed surprising, it was already known that roundworms, mites, flies, and polychaete worms also have proton-gated pLGICs (Schnizler et al., 2005; Mounsey et al., 2007; Beg et al., 2008; Juneja et al., 2014), and that humans, dogs and cows have Zn2+-gated pLGICs (Davies et al., 2003; Houtani et al., 2005). Also, although the slow activation of ELIC upon fast exposure to a saturating concentration of a high-efficacy agonist (time constant = 38 ± 6 ms for 50 mM cysteamine in outside-out patches excised from HEK-293 cells; 7 patches) may have seemed uncharacteristic, it was well-known that the mouse serotonin type-3A receptor (5-HT3AR) also activates with a time course that is uncharacteristically slow for a neurotransmitter-gated channel (time constant ≅ 7.3 ms in outside-out patches excised from HEK-293 cells; Mott et al., 2001). Moreover, although the presence of a ring of five pore-occluding phenylalanines may have seemed a unique property of ELIC, examination of invertebrate genomes had revealed the existence of several pLGIC subunits from roundworms, sea squirts, and sea anemones that also contain a phenylalanine at position 16′ of the transmembrane M2 α-helix and that may also assemble as homopentamers in vivo (Beg and Jorgensen, 2003; Satoh et al., 2006; Putnam et al., 2007; Rufener et al., 2009). Furthermore, nothing seemed particularly remarkable about the single-channel conductance values of GLIC and ELIC. Indeed, GLIC’s value (∼8 pS; Bocquet et al., 2007) is similar to that of the 5-HT3AR–5-HT3BR heteromer (∼16 pS; Davies et al., 1999), whereas ELIC’s (∼84 pS; Zimmermann and Dutzler, 2011) is nearly the same as that of the adult-muscle nAChR (∼73 pS; Elenes et al., 2009) for inward currents recorded under comparable ion conditions.
In contrast, here, we presented experimental evidence for the notion that the superfamily of nicotinic receptor-like channels is functionally more diverse than previously thought. We found that ELIC displays a remarkably atypical response to several perturbations that have long been used to probe the pLGICs from animals. Specifically, ELIC is not blocked by quaternary-ammonium cations TMA+ or TEA+ applied extracellularly at a concentration as high as 50 mM; displays a slower macroscopic current-decay time course in the presence of lidocaine and a saturating concentration of agonist; binds lidocaine in the cavities delimited by the M1, M2, and M3 transmembrane α-helices; and deactivates with wild-type-like kinetics—or even faster—upon the introduction of mutations that reduce the hydrophobicity of the channel’s transmembrane pore lining. Instead, all other cation-selective pLGICs studied to date are blocked by extracellular TMA+ and TEA+ in the micromolar-to-millimolar concentration range (Adler et al., 1979; Marshall et al., 1990; Zhang et al., 1995; Blanchet and Dulon, 2001; Akk and Steinbach, 2003), their current-decay time course in the presence of a saturating concentration of agonist becomes faster upon lidocaine binding (Alberola-Die et al., 2011) to the M2-delimited pore cavity (Hilf et al., 2010), and deactivate much more slowly upon engineering mutations such as those tested here (Filatov and White, 1995; Labarca et al., 1995; Kosolapov et al., 2000; Papke and Grosman, 2014).
Thompson et al. (2012) reported that TMA+ and TEA+ at a concentration of ∼20 mM reduce the peak-amplitude of currents mediated by ELIC by a factor of ∼2, a result that may seem to contradict our findings presented above. However, it should be noted that, in these experiments, ELIC currents were elicited by 1.6 mM GABA, a concentration of agonist that elicits only a half-maximal response (EC50). In contrast, in our experiments, ELIC currents were elicited by a saturating concentration of cysteamine. Thus, we surmise that the reported inhibitory effect of TMA+ and TEA+ reflects the competition of these quaternary-ammonium cations (which, when applied alone, do not elicit any measurable current response; Thompson et al. ) with GABA for the channel’s agonist-binding sites rather than the blocking of the pore. Indeed, it has been reported that TMA+ competes with cysteamine (Zimmermann et al., 2012) for the agonist-binding sites of ELIC, but their relative affinities are such that, at a concentration of cysteamine as high as 10 mM (i.e., the value we used in our experiments), TMA+ is expected to be outcompeted even when present at a concentration of 50 mM. Therefore, there is no discrepancy between our results and those reported by Thompson et al. (2012).
No discrete feature of the amino-acid sequence of ELIC seemed “suspicious” to us so as to anticipate its unique pore properties, and thus, it seems impossible to predict ELIC-like behavior in newly discovered pLGIC subunits on the basis of sequence gazing alone. Furthermore, because only a structural model for a nonconductive conformation of ELIC is currently available (Hilf and Dutzler, 2008), identifying the three-dimensional structural basis for its unusual open-channel behavior does not seem straightforward, either. Perhaps, however, the lack of block by 50 mM TMA+ or TEA+ would be the simplest property to explain if the most intracellular turn of each of the five M2 α-helices moved little upon opening (relative to their position in the nonconductive crystallized conformation) and formed the narrowest constriction in the open-channel conformation. Indeed, the diameter of the pore at the level of the intracellular ends of the M2 α-helices in the crystallized nonconductive conformation of ELIC is ∼9 Å, a size that would be consistent with the relatively high permeation rate observed for TMA+, TEA+, and TPA+. For comparison, in a homology model of the open-channel conformation of the mouse muscle nAChR, this diameter is ∼7 Å (Harpole and Grosman, 2014), whereas in the open-channel model of GLIC (Protein Data Bank ID code 4HFI; Sauguet et al., 2013), it is 5.2 Å (forcing the side chains of the five pore-lining glutamates at position –2′ to move out of the way by adopting a rotamer with χ1 ≅ –60° and χ2 ≅ +60°; otherwise, the diameter would be even narrower). It is tempting to speculate that it is this putatively larger size of the transmembrane pore cavity in the open state of ELIC that also underlies the lower binding affinity for TMA+ and TEA+, and the atypical effect of mutations to pore-lining residues on the kinetics of deactivation.
The physiological role of ELIC in the bacterium Erwinia chrysanthemi (now reclassified as Dickeya dadantii) is not known. However, we reason that the low affinity of ELIC for quaternary-ammonium cations, for example, would be advantageous if the channel faced an external solution that is rich in hydrophobic cations, and yet, the bacterium required a fast influx of inorganic cations—such as Na+ and Ca2+—through the channel’s cation-selective pore. In addition, it is likely that permeant organic cations for which the pore of ELIC has a higher affinity (higher than for TMA+, TEA+ or TPA+) exist. If this were the case, and if such cations indeed existed in the solutions bathing the bacterial membrane, then ELIC could function as a channel for solutes other than small inorganic cations. Clearly, extensive future work will be required to elucidate how the peculiar properties of this channel give rise to its biological function. Moreover, there may well be other pLGICs that behave like ELIC among the largely unexplored bacterial, archaeal, and invertebrate members of the superfamily, and thus, ELIC’s unconventional features may turn out to be less unusual than we currently think.
As for its role as a model system for the entire superfamily, our data indicate that caution should be exercised when generalizing results obtained with ELIC to the rest of the pLGICs, especially as they pertain to the properties of the transmembrane pore domain. Nevertheless, we would like to emphasize the notion that a member with such different functional properties in the context of a well-conserved three-dimensional architecture offers an unparalleled opportunity to further our mechanistic understanding of the pLGICs from vertebrate animals.
We thank Keith Brister, Joseph Brunzelle, David Smith, Elena Kondrashkina, Spencer Anderson, and Zdzislaw Wawrzak (Life Sciences Collaborative Access Team, 21 ID–D/F/G at Argonne National Laboratory, Advanced Photon Source) for technical assistance during crystallographic data collection; S. Sine (Mayo Clinic, Rochester, MN) for wild-type muscle nAChR subunit cDNAs; Y. Bello, C. Johnson, N. Kowalczyk, M. Rigby, and G. Wang for assistance with mutagenesis, cell culturing and transfections; and T. Harpole for the critical reading of the manuscript.
This work was supported by National Institutes of Health Grant NS042169 (to C. Grosman) and the Richard and Margaret Romano Professorial Scholarship (to C. Grosman).
The authors declare no competing financial interests.
Sharona E. Gordon served as editor.
- Abbreviations used in this paper:
- ligand-gated ion channel from Erwinia chrysanthemi
- ligand-gated ion channel from Gloeobacter violaceus
- nicotinic acetylcholine receptor
- pentameric ligand-gated ion channel
- Submitted: 25 November 2014
- Accepted: 14 May 2015
This article is distributed under the terms of an Attribution–Noncommercial–Share Alike–No Mirror Sites license for the first six months after the publication date (see http://www.rupress.org/terms). After six months it is available under a Creative Commons License (Attribution–Noncommercial–Share Alike 3.0 Unported license, as described at http://creativecommons.org/licenses/by-nc-sa/3.0/).