Mechanistic signs of double-barreled structure in a fluoride ion channel
Abstract
The Fluc family of F− ion channels protects prokaryotes and lower eukaryotes from the toxicity of environmental F−. In bacteria, these channels are built as dual-topology dimers whereby the two subunits assemble in antiparallel transmembrane orientation. Recent crystal structures suggested that Fluc channels contain two separate ion-conduction pathways, each with two F− binding sites, but no functional correlates of this unusual architecture have been reported. Experiments here fill this gap by examining the consequences of mutating two conserved F−-coordinating phenylalanine residues. Substitution of each phenylalanine specifically extinguishes its associated F− binding site in crystal structures and concomitantly inhibits F− permeation. Functional analysis of concatemeric channels, which permit mutagenic manipulation of individual pores, show that each pore can be separately inactivated without blocking F− conduction through its symmetry-related twin. The results strongly support dual-pathway architecture of Fluc channels.
https://doi.org/10.7554/eLife.18767.001Introduction
The F− ion, ubiquitous in the aqueous biosphere since the dawn of life, has been generally considered an irrelevant nonparticipant in membrane biology. This view was recently upended by the discovery of F−-specific riboswitches in many bacterial genomes (Baker et al., 2012), a breakthrough that over the past few years has revealed F−-transporting membrane proteins underpinning a widespread microbial physiology of resistance to environmental F− toxicity (Baker et al., 2012; Stockbridge et al., 2012; Li S et al., 2013; Ji et al., 2014; Smith et al., 2015). These F− exporters fall into two phylogenetically unrelated classes: energy-consuming anion transporters of the CLC superfamily (Stockbridge et al., 2012; Brammer et al., 2014), and electrodiffusive, F−-specific ion channels of the Fluc family (Stockbridge et al., 2013). While CLC proteins have been studied for decades (Miller, 2015), Flucs represent a novel and idiosyncratic class of ion channels. Bacterial Fluc channels assemble as dimers of small subunits (120–130 residues, 4 transmembrane helices) arranged in antiparallel transmembrane topology. This dual-topology architecture, originally inferred from biochemical and electrophysiological behavior (Stockbridge et al., 2013; Stockbridge et al., 2014), was recently observed in crystal structures of two homodimeric Fluc homologues (Stockbridge et al., 2015), as depicted in Figure 1a. The hourglass-shaped channel is complexed with two 'monobody' crystallization chaperones, which plug each of the wide vestibules facing the two aqueous solutions. (Figure 1—figure supplement 1). These ~10 kDa engineered monobodies are also high-affinity channel blockers whose long-lived block-times are readily observed in single-channel recordings (Stockbridge et al., 2014).
The structures suggested an additional surprising feature of Fluc proteins: double-barreled architecture. In contrast to most ion channels, in which a single ion-permeation pore spans the membrane along a central axis, Flucs were proposed to contain two pores running along the sides of the complex, with each pore built by residues from both subunits (Stockbridge et al., 2015). However, this is a crystallographically 'dry' inference based on structural results alone, unaccompanied by any functional behavior that would add mechanistic vitality to support or refute this otherwise unprecedented picture. An additional crystallographic ambiguity is that the crevices housing the anion densities are too narrow to show themselves at 2.1–2.6 Å resolution as clearly delineated aqueous pores traversing the protein. For these reasons, we consider that the argument for dual-pathway construction is not yet consummated; it calls for tests based on functional properties of the pores manipulated individually. By solving crystal structures of ion-permeation mutants of a bacterial Fluc channel and analyzing the behavior of fused-dimer constructs, we now unambiguously establish the two-pore character of Fluc channels.
Results
All experiments here employ the previously described 'Ec2' Fluc homologue from an E. coli virulence plasmid (Stockbridge et al., 2013). The homodimeric channel's symmetric, dual-topology structure demands that the pores adopt antiparallel orientations and that the two chemically distinct F− ions in one pore, designated F1− and F2−, are mirrored symmetrically in the other. In each pore, the two ions are coordinated by dipolar H-bond donors situated along a transmembrane 'polar track' (Figure 1b) that constitutes a possible ion pathway through the channel. Moreover, the F− coordination shells are completed by conserved phenylalanine residues interacting edge-on with the bound anions (Figure 1b), Phe80 associated with the F1− pair and Phe83 with the F2− pair. These four aromatic side chains adopt a striking side-to-face 'Phe-box' arrangement that supports an unusual aromatic-halide coordination motif proposed to be essential for F−-specific binding and permeation (Stockbridge et al., 2015). Edge-on aromatic-anion coordination geometry makes chemical sense in light of the electropositive hydrogens of the quadrupolar aromatic ring (Jackson et al., 2007; Philip et al., 2011; Schwans et al., 2013). The two antiparallel polar tracks and their associated bound F− ions are well-separated in the crystal structure, with no plausible pathway for ions to transfer between tracks. These structural features led to the proposal, to be tested here, that the two polar tracks represent two separate pores.
Phe-box residues are essential for F− binding and transport
Mutation of the conserved Phe-box residues in a Fluc homologue from B. pertussis was previously shown to severely inhibit F− transport, an effect presumed to reflect the replacement of F−-friendly binding sites in the conduction pathways with kinetic barriers to ion movement (Stockbridge et al., 2015). This idea is now tested in Ec2, where we here examine occupancy of the F−-binding sites by solving crystal structures of Phe-box mutants F80I and F83I (Figure 1c, Table 1) and assessing their transport function. Substitution of each aromatic side chain by Ile specifically extinguishes the corresponding pair of F− densities, one in each pore, while leaving the other pair Phe-coordinated and unperturbed. The mutations cause only local structural disturbances such as alternate rotamers at nearby sidechains Ser84 and Ser110 (Cα rmsd 0.4 Å, 0.2 Å for F80I, F83I). Functionally, the mutations are dramatic, reducing F− flux in reconstituted liposomes to undetectable levels (Figure 1d), a >104-fold rate inhibition compared to the wildtype (WT) channel (Stockbridge et al., 2013). These facts imply that Phe80 and Phe83 lie along F− permeation pathways and that the functional defect in the mutants stems directly from the loss of the aromatic ring as a F−-binding partner in these pathways.
Inactivation of individual pores
The results above point to the importance of the Phe-box side chains in F− binding and permeation, and the physical organization of their associated F− ions makes single-pore architecture implausible. But the homodimeric nature of Ec2 limits our ability to examine individual pores functionally, since any such mutation appears twice in the channel. We therefore constructed a tandem dimer by connecting two copies of the WT Ec2 sequence with a linker containing a foreign transmembrane helix to force an antiparallel arrangement of the two domains (Figure 2a, Figure 2—figure supplement 1). This construct, denoted 'WT/WT', mimics eukaryotic Fluc genes (Smith et al., 2015; Stockbridge et al., 2013), and it reprises an engineered homologue previously used to infer dual-topology assembly of Fluc channels (Stockbridge et al., 2013). The concatemer expresses poorly (~40 μg/L culture), but well enough for purification and reconstitution in liposomes and planar lipid bilayers (Figure 2b–d). The concatemer is functionally indistinguishable from the homodimeric WT channel with respect to unitary current, high intrinsic open probability, brief, infrequent closings, F− selectivity, and submicromolar affinity block by monobody proteins (Figure 2—figure supplements 1–3). Moreover, the F80I / F80I and F83I / F83I double mutants with identical Phe substitutions in both domains - mimics of the mutant homodimers above - are likewise inactive (Figure 2c). Concatemeric Ec2 thus faithfully recapitulates the behavior of the corresponding homodimers.
This construct allows us to make channel-disruptive mutations along individual pores, and to assess whether channel function is retained by the unperturbed pore. The complete inhibition of F− flux by either Phe mutant in the homodimer suggests a simple experimental design to test whether Fluc contains two functional pores. We use a liposome-based F− efflux assay under 'Poisson-dilution' conditions (Stockbridge et al., 2013; Maduke et al., 1999; Wu et al., 2007), where most reconstituted liposomes carry only a single Fluc concatemer. Since efflux through a single channel is much faster than the time response of the F− detection system, the assay provides a straightforward all-or-none indicator for the presence of functional channels; liposomes containing active channels release F− immediately upon initiating efflux, regardless of whether one or two pores are active, while liposomes devoid of active channels, as with the functionally disabled Phe mutants (Figures 1d, 2c), show no F− efflux.
The flux behavior of the Phe-box mutants in the fused-domain channels meet all expectations of two-pore assembly, as Figure 3 illustrates. All four single Phe mutants, each of which leaves one pore unaltered, score active in these efflux experiments (Figure 3a). Since the two Phe residues from the same domain contribute to different pores, and thus each pore harbors Phe80 from one domain and Phe83 from the other, we can distinguish two classes of double mutants: trans, wherein the two mutations reside in different pores, or cis, with both mutations in the same pore. Two trans double mutants, F80I / F80I and F83I / F83I, are completely inactive, as shown above, while both cis double mutants (F80I / F83I, F83I / F80I) retain channel activity. (We were unable to test trans double mutants in which both Phe residues in the same domain were mutated, as these were biochemically intractable.) These results show that channels with both pathways mutated are inactive irrespective of the particular combination of mutations, while those mutated in only one pathway remain active, as expected for two-pore construction.
Single-channel behavior of single-pore mutants
Although F−-efflux experiments clearly distinguish active from inactive channel mutants, their limited time-resolution (~1 s) precludes quantitative examination of F− permeation rates. Single-channel electrical recording in planar lipid bilayers provides such information by directly comparing F− current through WT channels and functionally active mutants. The results (Figure 4, Figure 4—figure supplement 1) are unambiguous: all single Phe-box mutants show unitary conductance of 5.7 + 0.3 pS, roughly half that of WT/WT fused-domain or WT homodimeric channels (10.6 + 0.1 pS). Moreover, the single-pore channels all recapitulate the high open probability, and monobody block of WT. We additionally note that the cis double mutants also show clear single-channel conductance that is nonetheless lower (2.3 + 0.2 pS) than the single mutants, possibly due to subtle, non-local structural disturbance of the conducting pore arising from the double mutation in the inactive pore. This ~2.5-fold drop in conductance, however, pales in comparison to the >104-fold drop observed when adding a second mutation in the trans configuration (Figure 3).
Discussion
By examining the behavior of the Fluc Ec2 channel at a mechanistic level, these experiments now firmly establish the two-pore architecture suggested by crystal structures. We consider this to be a property of the entire Fluc family, a conclusion further supported by indications of genetic drift in non-symmetric Flucs and the asymmetrical effects of mutations in yeast Fluc channels (Smith et al., 2015; Stockbridge et al., 2015). The results also highlight the importance of the conserved Phe-box residues in crafting coordination spheres for which F− ions are willing to shed their extensive aqueous hydration shells and enter the water-depleted confines of the ion-conduction pathway.
Two unexpected features of the mutant channels invite further comment. First, WT Ec2 channels show infrequent, brief excursions to a substate of ~6 pS, approximately one-half the conductance of the open channel (Figure 2d). It is natural to envision these as full closings of a single pore in a two-pore complex (Stockbridge et al., 2013), a picture demanding that such partial closings be absent in the single-pore mutants. However, this idea is plainly contradicted by the kinetically similar substates of ~3 pS appearing in the single-pore mutants here (Figure 4a,c). Instead, the substates in both wildtype and mutant channels must reflect a rare conformation that acts simultaneously on both pores to reduce F− current in each, perhaps occurring in the wide, water-filled vestibules on the two ends of the channel complex. Second, the striking and aesthetically appealing edge-to-face arrangement of the Phe-box in the WT channel had originally suggested to us that these four sidechains mutually stabilize and orient each other so as to properly coordinate the F− ions (Stockbridge et al., 2015). However, the Phe-mutant crystal structures here, while validating F− coordination by the aromatic rings, refute the idea of a mutually stabilizing quartet, since the two Phe side chains remaining in each Ile substitution, though distant from each other, adopt identical orientations and crystallographic order as in the full Phe-box of WT (side chain rmsd 0.2 Å upon backbone alignment).
Finally, despite the well-defined locations of the F− ions and the span of the polar track residues, we acknowledge that the detailed trajectories of the two pores remain uncertain due to their very narrow bore. However, as indicated in the conjectural cartoon of Figure 5, it is likely that the pores originate and end in the two aqueous vestibules that give the channel a symmetrical hourglass shape. This suggestion follows from monobody interaction with Fluc channels. In crystal structures of both Fluc homologues, monobodies cap and intrude into the vestibules with long loops bearing their diversified sequences (Stockbridge et al., 2015); in channel recordings here, the S9 monobody inhibits both double- and single-pore channels fully to the zero-current level, as though both pores are occluded simultaneously by the blocker. This proposal thus envisions the vestibules as common regions of low selectivity from which F− ions gain access to the separate F−-specific pathways.
Materials and methods
Biochemical
Request a detailed protocolFluc and Ec2 monobody S9 were expressed as previously described (Stockbridge et al., 2014; Stockbridge et al., 2015). For functional studies, purified Fluc was run over a S200 size-exclusion column equilibrated with 25 mM 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid (HEPES), 100 mM NaCl, 5 mM n-decyl-β-D-maltoside (DM), pH 7.0. Fluc channels were reconstituted into E. coli polar extract liposomes, and detergent was removed by dialysis against 25 mM HEPES, 300 mM KF, pH 7.0. For Cl- efflux assays, these liposomes were then diluted 1:1 with 25 mM HEPES, 300 mM KCl, pH 7.0, sonicated to clarity, and then freeze-thawed 3 times. For crystallography, Fluc and monobody S9 were purified in 100 mM NaF, 10 mM HEPES, pH 7.0, with the Fluc solution also containing 5 mM DM. All Ec2 variants were constructed with standard PCR techniques.
Fluc concatemer design
Request a detailed protocolThe base Ec2 construct used here, denoted 'WT', contains a single mutation (R25K) that increases expression but does not affect functionality (Stockbridge et al., 2015). The WT/WT Fluc concatemer was synthetically constructed with a hexahistidine tag, and designed so that the two Fluc domains would differ in nucleotide sequence while maintaining identical protein sequence. A transmembrane inversion linker based on a non-dimerizing glycophorin A variant (Lemmon et al., 1992) was used to join the two flux domains, as described for a Fluc heterodimeric homologue (Stockbridge et al., 2013).
Anion efflux assays
Request a detailed protocolF− and Cl− efflux assays were performed to detect anions released from pre-loaded proteoliposomes, as described (Stockbridge et al., 2015). Liposomes were prepared for Poisson-dilution conditions at low protein density (0.2 µg protein/mg lipid), such that on average only a single channel protein per proteoliposome is incorporated (Stockbridge et al., 2013), and ~50% of the liposomes are protein-free. The liposomes, loaded with 300 mM F− or 150 mM F− + 150 mM Cl- were extruded 21 times through a 0.4 µm membrane, exchanged into 25 mM HEPES, 300 mM K-isethionate, pH 7.0 and diluted 20-fold into a flux buffer of the same composition also containing 1 mM of the appropriate halide. Efflux was initiated by adding 1 μM valinomycin and was followed electrochemically by monitoring anion appearance in the suspension using a F− or Cl− electrode. All traces are shown as relative anion efflux, calculated as the voltage output from the electrode normalized to the efflux signal upon dissolution of the liposomes with 30 mM octylglucoside. All efflux experiments were repeated 2–4 times, and representative traces are shown in figures.
Single-channel recording
Request a detailed protocolSingle-channel recording was performed via a Nanion Orbit-mini planar bilayer system, using 70% 1-palmitoyl-2-oleoyl-phosphatidylethanolamine / 30% 1-palmitoyl-2-oleoyl-phosphatidylglycerol (Avanti Polar Lipids), 5 mg/mL in n-nonane to form bilayers. Single channels were inserted by addition of Ec2-reconstituted liposomes to the 'cis' side of the bilayer, and current was recorded at +200 mV holding potential (cis side defined as zero voltage) in symmetrical solutions containing 300 mM NaF, 10 mM NaCl, 15 mM MOPS-NaOH pH 7.0, with 75–150 nM of monobody S9 added to the cis chamber unless otherwise noted. Recordings were low-pass filtered at 160 Hz, digitized at 1.25 kHz, and analyzed in Clampfit 10 after further digital filtering at 100 Hz. Channel conductance, averaged from 10–14 separate single-channel records under each condition, was measured from the difference between open vs blocked current in each recording. All recordings shown in figures are representative of behavior seen on many channels in multiple reconstitutions.
Crystallography
Request a detailed protocolEc2 in complex with monobody S9 were crystallized as previously described (Stockbridge et al., 2015). Fluc dimers and S9 monobodies were individually purified and mixed at a 1:1.2 molar ratio at 10 mg/mL total protein. Crystallization well solutions were 100 mM N-(2-Acetamido)iminodiacetic acid, pH 6.2 (F80I) or 6.5 (F83I), 50 mM LiNO3, and 31% (F80I) or 36% (F83I) polyethylene glycol 600. Crystals were formed in hanging drops, mixing 1 μL each of protein and well solution, and incubating at 22°C for 14–21 days. Crystals were frozen in liquid nitrogen.
Datasets were collected at Advanced Light Source beamline 8.2.2. Diffraction images were processed in iMosflm (Battye et al., 2011), and data were merged and scaled using Aimless (Evans, 2011). The complete F80I dataset consisted of three different datasets from different positions along a single long crystal, merged in Aimless. Initial structure solution of both Phe mutants was done via molecular replacement using Phaser (McCoy et al., 2007) with PDB #5A43 as search model. Refinement was done with Refmac5 (Winn et al., 2003) and Phenix (Adams et al., 2010), using TLS refinement and local NCS averaging. Re-refinement of the previously collated WT Ec2 structure (Stockbridge et al., 2015) additionally used SAD data directly for additional phase information. Real-space refinement was done in COOT (Emsley et al., 2010), and Molprobity (Chen et al., 2010) was used to carry out model validation. Unbiased maps shown in Figure 1 are mFo-DFc maps output by Refmac5 using a model that had never had any fluorides introduced at any point in the refinement. Datasets from second crystals for both F80I and F83I that diffracted to similar resolution reproduce all the major structural findings here, particularly with regard to the pattern of bound F− and the lack of perturbation of the overall Fluc structure.
References
-
PHENIX: a comprehensive python-based system for macromolecular structure solutionActa Crystallographica Section D Biological Crystallography 66:213–221.https://doi.org/10.1107/S0907444909052925
-
iMOSFLM: a new graphical interface for diffraction-image processing with MOSFLMActa Crystallographica Section D Biological Crystallography 67:271–281.https://doi.org/10.1107/S0907444910048675
-
F-/Cl- selectivity in CLCF-type F-/H+ antiportersJournal of General Physiology 144:129–136.https://doi.org/10.1085/jgp.201411225
-
MolProbity: all-atom structure validation for macromolecular crystallographyActa Crystallographica Section D Biological Crystallography 66:12–21.https://doi.org/10.1107/S0907444909042073
-
Features and development of CootActa Crystallographica Section D Biological Crystallography 66:486–501.https://doi.org/10.1107/S0907444910007493
-
An introduction to data reduction: space-group determination, scaling and intensity statisticsActa Crystallographica Section D Biological Crystallography 67:282–292.https://doi.org/10.1107/S090744491003982X
-
A preference for edgewise interactions between aromatic rings and carboxylate anions: the biological relevance of anion-quadrupole interactionsJournal of Physical Chemistry B 111:8242–8249.https://doi.org/10.1021/jp0661995
-
Bacterial fluoride resistance, Fluc channels, and the weak acid accumulation effectJournal of General Physiology 144:257–261.https://doi.org/10.1085/jgp.201411243
-
Glycophorin a dimerization is driven by specific interactions between transmembrane alpha-helicesJournal of Biological Chemistry 267:7683–7689.
-
High-level expression, functional reconstitution, and quaternary structure of a prokaryotic ClC-type chloride channelJournal of General Physiology 114:713–722.https://doi.org/10.1085/jgp.114.5.713
-
Phaser crystallographic softwareJournal of Applied Crystallography 40:658–674.https://doi.org/10.1107/S0021889807021206
-
Yeast Fex1p Is a constitutively expressed fluoride channel with functional asymmetry of its two homologous domainsJournal of Biological Chemistry 290:19874–19887.https://doi.org/10.1074/jbc.M115.651976
-
Proof of dual-topology architecture of Fluc F- channels with monobody blockersNature Communications 5:5120.https://doi.org/10.1038/ncomms6120
-
Macromolecular TLS refinement in REFMAC at moderate resolutionsMethods in Enzymology 374:300–321.https://doi.org/10.1016/S0076-6879(03)74014-2
-
Uncoupling and turnover in a Cl-/H+ exchange transporterJournal of General Physiology 129:317–329.https://doi.org/10.1085/jgp.200709756
Article and author information
Author details
Funding
Howard Hughes Medical Institute
- Ludmila Kolmakova-Partensky
National Institute of General Medical Sciences (RO1-GM107023)
- Ludmila Kolmakova-Partensky
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Dr Randy Stockbridge for designing the Fluc concatemer and for advice in the early stages of this project.
Copyright
© 2016, Last 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
-
- 1,957
- views
-
- 457
- downloads
-
- 30
- 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
Images taken by transmission electron microscopes are usually affected by lens aberrations and image defocus, among other factors. These distortions can be modeled in reciprocal space using the contrast transfer function (CTF). Accurate estimation and correction of the CTF is essential for restoring the high-resolution signal in cryogenic electron microscopy (cryoEM). Previously, we described the implementation of algorithms for this task in the cisTEM software package (Grant et al., 2018). Here we show that taking sample characteristics, such as thickness and tilt, into account can improve CTF estimation. This is particularly important when imaging cellular samples, where measurement of sample thickness and geometry derived from accurate modeling of the Thon ring pattern helps judging the quality of the sample. This improved CTF estimation has been implemented in CTFFIND5, a new version of the cisTEM program CTFFIND. We evaluated the accuracy of these estimates using images of tilted aquaporin crystals and eukaryotic cells thinned by focused ion beam milling. We estimate that with micrographs of sufficient quality CTFFIND5 can measure sample tilt with an accuracy of 3° and sample thickness with an accuracy of 5 nm.
-
- Structural Biology and Molecular Biophysics
Agonists enhance receptor activity by providing net-favorable binding energy to active over resting conformations, with efficiency (η) linking binding energy to gating. Previously, we showed that in nicotinic receptors, η-values are grouped into five structural pairs, correlating efficacy and affinity within each class, uniting binding with allosteric activation (Indurthi and Auerbach, 2023). Here, we use molecular dynamics (MD) simulations to investigate the low-to-high affinity transition (L→H) at the Torpedo α−δ nicotinic acetylcholine receptor neurotransmitter site. Using four agonists spanning three η-classes, the simulations reveal the structural basis of the L→H transition where: the agonist pivots around its cationic center (‘flip’), loop C undergoes staged downward displacement (‘flop’), and a compact, stable high-affinity pocket forms (‘fix’). The η derived from binding energies calculated in silico matched exact values measured experimentally in vitro. Intermediate states of the orthosteric site during receptor activation are apparent only in simulations, but could potentially be observed experimentally via time-resolved structural studies.