Ion- and water-binding sites inside an occluded hourglass pore of a trimeric intracellular cation (TRIC) channel

Background Trimeric intracellular cation (TRIC) channels are crucial for Ca2+ handling in eukaryotes and are involved in K+ uptake in prokaryotes. Recent studies on the representative members of eukaryotic and prokaryotic TRIC channels demonstrated that they form homotrimeric units with the ion-conducting pores contained within each individual monomer. Results Here we report detailed insights into the ion- and water-binding sites inside the pore of a TRIC channel from Sulfolobus solfataricus (SsTRIC). Like the mammalian TRIC channels, SsTRIC is permeable to both K+ and Na+ with a slight preference for K+, and is nearly impermeable to Ca2+, Mg2+, or Cl–. In the 2.2-Å resolution K+-bound structure of SsTRIC, ion/water densities have been well resolved inside the pore. At the central region, a filter-like structure is shaped by the kinks on the second and fifth transmembrane helices and two nearby phenylalanine residues. Below the filter, the cytoplasmic vestibule is occluded by a plug-like motif attached to an array of pore-lining charged residues. Conclusions The asymmetric filter-like structure at the pore center of SsTRIC might serve as the basis for the channel to bind and select monovalent cations. A Velcro-like plug-pore interacting model has been proposed and suggests a unified framework accounting for the gating mechanisms of prokaryotic and eukaryotic TRIC channels. Electronic supplementary material The online version of this article (doi:10.1186/s12915-017-0372-8) contains supplementary material, which is available to authorized users.


Background
The regulated processes of Ca 2+ release from the intracellular stores and its uptake from the cytosol are vital for various biological processes including muscle contraction, neurotransmitter release, cell division, and apoptosis [1,2]. For instance, muscle contraction is initiated by membrane depolarization followed by opening of a ryanodine receptor (RyR) channel to release Ca 2+ from the lumen of the sarcoplasmic reticulum (SR) into the cytosol. The process is known as excitationcontraction (E-C) coupling [3]. Rapid efflux of Ca 2+ from the SR generates a transient negative potential inside the SR lumen and will hinder Ca 2+ release if the transmembrane potential remains unbalanced. Thus, efficient operation of E-C coupling requires not only the RyR to release Ca 2+ , but also counteracting ion channels to restore the balance of the SR membrane potential and maintain ion homeostasis within the SR lumen [4].
Two isoforms of SR/endoplasmic reticulum (ER) membrane proteins, called trimeric intracellular cation (TRIC) channels (TRIC-A and TRIC-B), presumably function as the counteracting ion channels facilitating the intracellular Ca 2+ handling processes [5]. Alternatively, they might serve to restore the balance of trans-SR K + after the RyRs close, instead of carrying countercurrent during Ca 2+ release [6]. TRIC channels are permeable to K + and Na + with moderate selectivity for K + over Na + and are impermeable to Ca 2+ , Mg 2+ , or anions [5,7]. TRIC-A and TRIC-B have distinct single-channel conductances as well as diverse regulatory mechanisms and physiological roles [8,9]. TRIC-A is regulated by transmembrane voltage [8] and may interact with the RyR functionally and physically [9,10]. TRIC-B channel is activated by micromolar Ca 2+ applied on the cytosolic side, but it is inhibited by Ca 2+ on the luminal side [8] and is also regulated by voltage [11]. TRIC-B may modulate the Ca 2+ -release channel activity of inositol trisphosphate (IP 3 ) receptor [9] and is essential for perinatal lung maturation [12]. Genetic mutations of the human TMEM38B gene (encoding the TRIC-B protein) are found in patients with a hereditary brittle bone disease called osteogenesis imperfecta [13][14][15][16]. TMEM38B knockout mice are deficient in producing collagen, and bone mineralization is impaired in the mutant animals [17,18]. Moreover, TRIC-A contributes to the maintenance of normal blood pressure and may serve as a potential pharmaceutical target for treating hypertension [19,20]. The association of TRIC channels with bone, pulmonary, and muscular diseases indicates that they have indispensable functions in the related physiological and developmental processes [9,17,18,21].
The structures of TRIC-B channels from Caenorhabditis elegans (CeTRIC-B1 and CeTRIC-B2) revealed a homotrimeric membrane protein-lipid complex with hourglass-shaped pores traversing through each individual monomer [22]. A phosphatidylinositol 4,5-biphosphate (PtdIns(4,5)P 2 , also known as PIP 2 ) lipid molecule binds specifically to each monomer at 1:1 stoichiometry, mediates trimerization of CeTRIC-B channels, and stabilizes the cytoplasmic gate of the channel. In addition to the members found in eukaryotic organisms, TRIC channel orthologs are widespread in prokaryotes including bacteria and archaea [23]. The prokaryotic TRIC members form the largest group of the TRIC family, outnumbering eukaryotic ones. A recent study suggested that they may mediate K + uptake, and the structures of two prokaryotic TRIC orthologs from Rhodobacter sphaeroides (RsTRIC) and Sulfolobus solfataricus (SsTRIC) were reported at 3.4and 2.6-Å resolution, respectively [24]. These prokaryotic orthologs form homotrimeric units resembling those of CeTRIC-B channels, and the ion-conducting pore is also contained within each individual monomer [24].
Although significant progress has been made on structural and functional studies of both eukaryotic and prokaryotic TRIC channels, the fundamental mechanistic questions regarding the molecular basis of ion selectivity remain largely open. The pore architectures of TRIC channels clearly do not resemble those of homotetrameric K + channels, such as the well-studied KcsA channel [25]. We do not know where the ion selectivity filter is located in TRIC channels or how TRIC channels selectively bind monovalent cations in their pore region. Therefore, it is indispensable to characterize the ionbinding sites along the permeation pathway of TRIC channels through high-resolution structural studies of TRIC channels. Moreover, the lack of the PIP 2 molecule in prokaryotic TRIC orthologs raises questions about how their gates are stabilized in the absence of PIP 2 in their structure. Here we describe the detailed ion-and water-binding sites inside the pore of each SsTRIC monomer, define the chemical basis of K + coordination inside a buried asymmetric filter-like structure, and provide new insights into a Velcro-like plug-pore interacting model accounting for the gating mechanism of TRIC channels.

Results
Ion-and water-binding sites within an hourglass-shaped pore As shown in Fig. 1, the SsTRIC protein reconstituted on giant unilamellar vesicles (GUVs) forms an ion channel permeable to K + ion and exhibits three major open states with conductances at 42.2 pS, 100.6 pS, and 163.7 pS, respectively. These conductances are similar to those of the CeTRIC-B1 channel (45.8 pS, 101.5 pS, and 153.2 pS for the first, second, and third open states, respectively) reported previously [22], suggesting that their open pores may have similar sizes. Furthermore, the SsTRIC channel is also permeable to Na + with a P K /P Na ratio at 1.21 (derived from the reversal potential at 4.84 mV under the 210 mM KCl/210 mM NaCl bi-ionic condition; see Fig. 2a, b), close to the P K /P Na value (1.5) of the mammalian TRIC-A channel reported previously [5]. Meanwhile, SsTRIC is essentially impermeable to Ca 2+ , Mg 2+ , or Cl -, as the reversal potential remains close to 0 mV when 75 mM CaCl 2 or MgCl 2 is added to the bath solution (the KCl concentration is the same on both sides; see Fig. 2c, e-g).
The structure of SsTRIC solved at 2.2-Å resolution shows that it forms a symmetrical homotrimer measuring~56 Å wide and~40 Å tall, and has flat surfaces on both the extracellular and cytoplasmic sides ( Fig. 3a and b). It superposes well with the previous structures of SsTRIC in complex with the Fab fragment of a monoclonal antibody (SsTRIC-Fab complex, root mean square deviation (RMSD) of C α atoms = 0.39 Å), and a bacterial TRIC homolog (RsTRIC) from Rhodobacter sphaeroides [24] (RMSD of C α atoms = 2.36 Å) (see Additional file 1: Figure S1). When compared to CeTRIC-B1 or the RyR (both are membrane proteins located on the ER/SR membrane), the membrane-embedded domain of the SsTRIC trimer appears to be~6 Å taller along the dimension parallel to the membrane normal (Additional file 2: Figure S2). Curiously, the tilt angles of transmembrane helices in SsTRIC with respect to the membrane normal are dramatically different from those in CeTRIC-B1/B2 (Table 1).
Except for M6, most transmembrane helices (M1-M5 and M7) in SsTRIC are tilted from the membrane normal at smaller angles compared to those in CeTRIC-B1/B2. This may arise from the adaptation of membrane proteins to different membrane environments, i.e., the plasma membrane of prokaryotic cells versus ER membranes of eukaryotic cells.
Both eukaryotic and prokaryotic TRIC channels harbor their ion-conducting pores within each subunit of the homotrimeric assemblies [22,24]. As shown in Fig. 3c, an hourglass-shaped pore runs through each SsTRIC monomer from cytoplasm to extracellular space. Along the pore, four bottleneck sites, including two (G1 and G2) on the cytoplasmic side, one at the central region (G3), and one close to the extracellular vestibule (G4), form highly constricted areas with extremely small widths at 0.7-2.0 Å (smaller than the diameter of a K + ion at~2.8 Å [26]). For the pore in RsTRIC, the G4 site around Phe103 appears to be slightly wider than the one in SsTRIC, although their overall profiles and pore-lining residues are very similar (Additional file 3: Figure S3a and b). For comparison, the pore within each CeTRIC-B monomer harbors two bottlenecks (instead of four), one on the cytoplasmic side and the other on the luminal side of the ER/SR [22] (Additional file 3: Figure S3c). The pore lumen surface of SsTRIC is mainly shaped by amino acid residues from M1, M2, M4, and M5 helices. Within the pore, a single file of eight well-resolved spherical densities was observed in the high-resolution electron density map of the K + -soaked crystal (Additional file 4: Figure S4a and b). The four ion/water molecules (P0-P3) on the cytoplasmic side are related to the other four at the P0′-P3′ sites on the extracellular side, and they approximately follow the internal pseudo-C 2 symmetry [22,24] between the M1-M2-M3 and M4-M5-M6 helices. Among these densities, two (P0 and P0′) are located at the pore center area between the G3 and G4 bottleneck sites (Fig. 3c), and they intercalate in the grooves of glycine-rich kinks around Ala44-Gly46 and Gly130-Gly134 regions on M2 and M5, respectively ( Fig. 3e and f). Previously, two spherical F o -F c densities at positions similar to P0 and P0′ were observed in the structure of the SsTRIC-Fab complex and assigned as water molecules [24].
To test whether any of these potential ion-/waterbinding sites inside the pore can be accessed by monovalent cations, we have soaked the SsTRIC crystal with Tl + , the electron-rich surrogate of the K + ion [27]. The anomalous difference Fourier peaks of Tl + ions are mostly detected on the P0′ and P3′ (outer) sites (Tl-1 and Tl-2 shown in Additional file 4: Figure S4c) inside . c The current-voltage relationships of the SsTRIC channel at O1, O2, and O3 states. The error bars indicate the standard errors of mean values (SEM). The number of measurements (n) is labeled above or below the data points the pore (and on numerous surface sites), but not on the P0-P3 (inner) sites. By referring to the positions of the Tl + peaks, the putative K + binding sites have been assigned in the structure of SsTRIC soaked with KCl. As shown in Fig. 3e and f, the K + ion bound on the P0′ site is mainly coordinated by the backbone carbonyl from Thr41, backbone amide from Gly45, and a water molecule nearby on the P0 site. Such a coordination mode resembles one half of the K + coordination observed previously in valinomycin (a K + -permeable ionophore, see Additional file 5: Figure S5a and b), but differs from the square antiprism-type coordination in KcsA (Additional file 5: Figure S5a and c). Due to steric hindrance of the bulky side chain of Phe104, three additional ligands below K + (as in valinomycin) are absent in SsTRIC. The water molecule on the P0 site is further ligated to the backbone carbonyl groups of Val129 and Gly130 and the backbone amide of Gly133. The coordination bond lengths between Representative recordings and current-voltage relationships of the wild-type SsTRIC measured under bi-ionic condition of 210 mM KCl (pipette) and 210 mM NaCl (bath). The I/V curve of wild-type SsTRIC measured under symmetrical 210 mM KCl (red) is included for comparison. The reversal potential (4.84 ± 0.86 mV for bi-ionic condition and −0.02 ± 0.63 mV for symmetrical condition) is derived by linear fitting of the I/V relationship data. The error bars indicate the standard errors of mean values (SEM, number of patches/n = 5 for bi-ionic condition, and n = 8 or 9 for symmetrical condition as shown in Fig. 1c). c, d, and e Representative recordings and current-voltage relationships of the wild-type SsTRIC (dark) and F104 A mutant (blue) measured under bi-ionic condition of 210 mM KCl (pipette)/ 210 mM KCl + 75 mM CaCl 2 (bath). The reversal potential is at 0.04 ± 0.09 mV (mean ± SEM, n = 4 except that n = 3 for the data point at −20 mV) for the wild type and 0.11 ± 1.77 mV (mean ± SEM, n = 5) for the F104A mutant. f and g Representative recordings and current-voltage relationships of the wild-type SsTRIC measured under bi-ionic condition of 210 mM KCl (pipette)/210 mM KCl + 75 mM MgCl 2 (bath). The reversal potential is at −0.04 ± 0.69 mV (mean ± SEM, n = 3). h and i Representative recordings and current-voltage relationships of F104A mutant measured under bi-ionic condition of 210 mM KCl (pipette)/210 mM NaCl (bath). For comparison, the I/V curve of wild-type tSsTRIC under the same condition is included i The reversal potential is at 1.92 ± 0.93 mV (mean ± SEM, n = 5 except that n = 3 for −20 mV and n = 4 for 40 mV) for F104A mutant or 4.84 ± 0.86 mV (mean ± SEM, n = 5) for the wild type. See Additional file 11 for the supporting data values K + and its ligands on P0′ site are consistent with those (2.7-2.9 Å) observed in the KcsA channels, although the coordination geometry of K + in SsTRIC is of trigonal shape instead of the square antiprism found in KcsA [25]. The introduction of a G45A/G47A or G132A/G133A mutation in RsTRIC (corresponding to A44/G46 and G133/ 134 in SsTRIC) abolishes its channel activity [24], indicating that the kinks on M2 and M5 helices are essential for the function of TRIC channels. In these two mutants, the kinks may have been distorted, leading to loss of the P0′ or P0 site (both crucial for K + binding). Similar kinks are also present on the M2 and M5 helices of CeTRIC-B channels, and a Rb + ion was found on a site nearby these regions [22]. Thus, such a double-kink structure with an inverted twofold pseudo-symmetry may represent a general characteristic feature of the TRIC channel family.
In SsTRIC, the aromatic rings of Phe16 from the M1 helix and Phe104 from the M4 helix enclose the P0 and P0′ sites from the other side opposing to the M5 kink and M2 kink, respectively ( Fig. 3e and f ). The side chains of Phe104 and Phe16 not only serve to constrict the pore center to form two bottlenecks (G3 and G4 in Fig. 3c) but may also contribute to binding K + through cation-π interactions [28]. Mutation of Phe104 to alanine decreases the P K /P Na ratio from 1.21 to 1.08 (as the reversal potential shifts from 4.84 to 1.92 mV under KCl/NaCl bi-ionic conditions, see Fig. 2h and i), and the F104A mutant is also nearly impermeable to Ca 2+ or Clas the wild type ( Fig. 2d and e). Therefore, Phe104 may have a crucial role in adjusting the permeability of SsTRIC to K + versus Na + , presumably by shaping the geometry of the P0′ site slightly more favorably for binding K + than Na + . Corresponding to the two phenylalanine residues in SsTRIC, CeTRICB1/B2 channels have their His34/35 and Lys136/Lys137 positioned nearby the kinks to enclose their central ion-binding sites [22]. Such a putative asymmetric cation-binding site at the pore center likely serves as the basis for TRIC channels to selectively bind and conduct monovalent cations (K + and Na + ) when the flanking bottleneck sites become wider at the open state. Besides the K + /Tl + ion on the P0′ site, the one on the P3′ site in the extracellular vestibule of the pore is coordinated by the side chain carboxyl of Asp50 and the backbone carbonyl of Gly46. This peripheral site is close to the extracellular surface, and it may serve to attract K + ions before they reach the pore center, or facilitate the release of K + to the extracellular space.

Gain-of-function mutants
To further verify whether the ion permeation pathway of the SsTRIC channel does traverse through each monomer as observed in the structure, we have introduced a series of alanine mutations to the key amino acid residues lining the pore lumen surface (Additional file 6: Figure S6a). Expression of the wild-type SsTRIC channel in Escherichia coli induces moderate inhibition of cell growth, while the empty vector or expression of a different membrane protein (the large-conductance mechanosensitive channel from E. coli, EcMscL) has little effect (See figure on previous page.) Fig. 3 Ion-and water-binding sites within the pores of SsTRIC channel. a Structure of a SsTRIC homotrimer at 2.2-Å resolution viewed along the membrane plane. The estimated position of membrane is indicated by the light gray boxes. The K + ions and water molecules within the pore region are shown as purple and red spheres, respectively. b The SsTRIC trimer viewed along the C3 axis from the cytoplasmic side. The green meshes indicate the 2F o -F c electron density (contoured at +1.0 × σ level) of a neutral lipid molecule. c The pore running through SsTRIC monomer is indicated by the tri-color dotted areas. Red areas, radius <1.15 Å; green areas, 1.15 Å < radius < 2.30 Å; blue areas, radius >2.30 Å. G1-G4 label the pore constriction sites (or gates) as shown by the red areas, while the green areas contain the binding sites for ion or water molecules labeled as P1′, P0′, P0, P1-P3 in d.
The key amino acid residues lining the pore lumen are shown as stick models. From middle to bottom, they are F104, F16, Y153, R137, D138, and M146. d Distribution of K + ions or water molecules within the pore region of an SsTRIC monomer. The sectional view of a monomer is shown as the electrostatic potential surface model overlaid with the cartoon model (color code: deep red, electronegative; deep blue, electropositive; silver, cartoon model of the protein backbone). The distances between two adjacent molecules are shown by the numbers (Å) labeled around the dashed lines. e A zoom-in view of the pore center with two ion-/water-binding sites. K + ions and water molecules are shown as silver and red bullets, respectively. Oxygen, nitrogen, and carbon atoms of the amino acid residues are colored red, blue, and yellow, respectively. The interactions between K + ion on P0′ site (or water molecule on P0 site) and the surrounding residues are indicated by the dashed lines with their distances (Å) labeled nearby. f The pore center viewed from a different angle to show the water molecule on P0 site and its surrounding residues The tilt angle (η) is defined as the angle between the axis of each individual transmembrane helix and the central C3 axis (parallel with the membrane normal) of the SsTRIC trimer or CeTRIC-B1/B2 trimer on cell growth (Additional file 6: Figure S6b and c). Such an inhibitory effect is likely due to the basal leaky activity of the SsTRIC channel stimulated by the resting membrane potential (Δψ) of E. coli cells. The activated SsTRIC channel may leak out intracellular K + ions, perturb the resting potential balance, and inhibit the growth of E. coli cells moderately (Additional file 6: Figure S6b). Mutation of the key residues along the pore lumen surface (such as F104A, D97A, R137A, D138A, M146A, and Y153A) leads to gain-of-function (GOF) phenotypes compared to the wild-type channel (Additional file 6: Figure S6b), resembling the phenotype of a wellestablished severe GOF mutant (G26H) of MscL [29,30]. Unlike the other mutants, the F16A protein is undetectable in the membrane fraction (Additional file 6: Figure S6d) or at the whole-cell level, indicating that the protein may have been degraded before forming a functional channel on the membrane due to its high toxicity. As a result, the cells hosting the F16A mutant (not expressing) grow nearly as normal as the empty vector. The control, the R187A mutant (at a site distant from the pore region), expresses normally and shows phenotypes similar to the wild type (Additional file 6: Figure S6c). Hence, the in vivo functional assay data are consistent with the observation of a potential ion permeation pathway contained within each SsTRIC monomer.

Neutral lipid instead of PIP 2 bound to SsTRIC
At the monomer-monomer interfaces of the SsTRIC trimer, lipid-like moieties have been observed ( Fig. 3a and b), but the identities of these interfacial cofactors are elusive [24]. They have elongated tubular features extending from the peripheral region toward the center of the SsTRIC trimer, and three fatty acyl group-like chains join at the center to form a triskelion-shaped structure enhancing the association of three SsTRIC monomers through hydrophobic interactions. Thin layer chromatography and mass spectrometry analyses of the lipid samples extracted from purified SsTRIC preparations suggest that they belong to highly hydrophobic neutral lipid species (likely a mixture of triacylglycerol and free fatty acids/derivatives from the E. coli membrane, Additional file 7: Figure S7a-c). A similar pattern of lipid bands is also observed in the sample extracted from the pure preparation of TRIC ortholog from E. coli (EcTRIC, Additional file 7: Figure S7a). These lipid molecules are not random species bound to SsTRIC. Instead, they are most likely specific lipids selected by the hydrophobic binding site with a well-defined triskelion shape. Electrospray ionization mass spectrometry analysis detected the presence of a molecule with an m/z value of 987.7, likely corresponding to an ionized triacylglycerol molecule, in the lipid samples extracted from purified SsTRIC protein preparation (Additional file 7: Figure S7b and c). Therefore, the triskelion-shaped density in the SsTRIC trimer is interpreted as a triacylglycerol molecule, and the model fits reasonably well with the electron density (Additional file 7: Figure S7d and e). The amino acid residues involved in binding the lipid are mainly hydrophobic residues from the M2 and M5 helices (Additional file 7: Figure S7e and f). In comparison, CeTRIC-B channels contain a PIP 2 lipid molecule per monomer, and the two fatty acyl chains of the PIP 2 molecule extend to the monomer-monomer interfaces (Additional file 7: Figure S7g and h), both contributing to trimerization simultaneously [22]. Like the lipid in SsTRIC, the fatty acyl chains of PIP 2 molecules in CeTRIC-B channels are also surrounded by hydrophobic residues from the M2 and M5 helices (Additional file 7: Figure S7h and i). Nevertheless, the bulky hydrophilic inositol 4,5-biphosphate head group (covalently linked to the 3-position of the glycerol backbone) of PIP 2 inserts through the gap between the M5 and M6 helices, and is wrapped at the center of each monomer to stabilize the cytoplasmic gate of CeTRIC-B channels (Additional file 7: Figure S7g and h). For the lipid molecule in SsTRIC, the three fatty acyl groups (instead of two) extend laterally through the fenestration between M6 and M2′ (from adjacent monomer) helices. Instead of bending toward the pore region within each monomer (like the PIP 2 head group), they have reached the external surface facing lipid bilayer (Additional file 7: Figure S7d and e). The two positively charged residues (Lys130 and Arg134 in CeTRIC-B2 shown in Additional file 7: Figure S7h) involved in binding the PIP 2 head group in CeTRIC-B are absent in SsTRIC or other prokaryotic orthologs. No bulky hydrophilic head groups are attached to the fatty acyl chains of the lipid found in the SsTRIC channel, suggesting that prokaryotic TRIC channels may adopt a different mechanism to stabilize the cytoplasmic gate. It is noteworthy that Sulfolobus solfataricus and other archaea contain etherlinked lipids (such as sn-2,3-diphytanylglycerol diether and glycerol-dialkyl-calditol-tetraether) [31] instead of the esterlinked lipids commonly found in bacterial or eukaryotic cells. Therefore, the SsTRIC from its native cellular environments may adopt the ether lipids in its structure.
A tethered plug-like motif at the cytoplasmic gate How does SsTRIC stabilize its cytoplasmic gate in the absence of a bulky lipid head group? The pore of the SsTRIC channel is more open on the extracellular side than the intracellular side (Fig. 4a). What is the cause of such dramatic differences between the outer and inner vestibules of the SsTRIC channel? Each SsTRIC monomer consists of three motifs, called motif A, motif B, and motif C, respectively (Additional file 8: Figure S8). Motif A contains M1, M2, M3, and two short amphipathic helices. Motif B comprises M4, M5, M6, and one short amphipathic helix in the M5-M6 loop region.
These two motifs correspond to the N-repeat and Crepeat described in the previous study [24]. Motif A is related to motif B by a pseudo-C2 axis running approximately parallel to the membrane plane (Additional file 8: Figure S8a and b). They superpose well with the root mean square deviation (RMSD) of α-carbon atoms at 1.42 Å (Additional file 8: Figure S8b), and their sequences share 25% identity (Additional file 8: Figure S8c). The side chains of three key residues from motif B, namely Tyr153, Met146, and Arg137, collectively contribute to occlusion of the inner vestibule of the pore (Additional file 8: Figure S8d). On the other side, their symmetry-related counterparts in motif A are either oriented away from the pore region (Tyr66 A versus Tyr153 B and Gln59 A versus Met146 B ) or replaced by Ala residue (Ala49 A versus Arg137 B ) (Additional file 8: Figure S8e). The insertion of Pro57 (in a cis configuration) within the M2-M3 loop region directs the side chains of Gln59 and Tyr66 away from the pore region to create an open cavity on the extracellular side (Additional file 8: Figure S8e). Thus, these differences between motifs A and B result in the asymmetric status of the SsTRIC channel with its outer vestibule open and inner vestibule closed. On the cytoplasmic side, the region between the M5 and M6 helices forms a loop-helix-loop motif (residues Asn142-Tyr153) (Fig. 4b). It plugs in the inner vestibule of the hourglass-shaped pore and contributes to the formation of the first two bottlenecks (G1 and G2 gates shown in Fig. 3c) within the cytoplasmic vestibule. This plug-like motif binds to Asp97, Arg137, and Asp138 residues lining the surface of the pore lumen, through hydrogen bonds and ionic interactions (Fig. 4c). Thereby, it is stabilized in a position to block the access of ions into the pore (Fig. 4a). Removal of the plug-like motif from the structure generates a wide-open cavity in the cytoplasmic vestibule region (Fig. 4d). In this model, the open Among the residues from the plug motif, Tyr153 is closest to the pore center near Phe16, and forms a strong hydrogen bond with Asp97 (Fig. 5a). Strikingly, the position of the Tyr153 bulky side chain coincides with that of the inositol 4,5-biphosphate head group of the PIP 2 molecule in CeTRIC-B1, when the structure of SsTRIC is superposed with that of CeTRIC-B1 (Fig. 5b). In CeTRIC-B1, the head group of the PIP 2 molecule stabilizes the cytoplasmic gate through binding to two basic residues on the pore lumen [22]. Likewise, Tyr153 in SsTRIC may have a similar role in controlling the cytoplasmic entrance through its bulky side chain, and may behave like a plug to the hourglass bottleneck. The electrophysiological analyses on the Y153A mutant of SsTRIC confirm its role as a crucial gating element of the channel. The mutant appears more active than the wild type (Fig. 5c). Remarkably, the right shift of the amplitude histogram peaks of Y153A in respect to those of the wild type indicates that the conductance of the Y153A mutant is larger than that of the wild type (Fig. 5d). Moreover, the open probability (NP o ) of Y153A is also higher than that of the wild type (Fig. 5e). The open state of the Y153A mutant also has a longer lifetime than that of wild type (Fig. 5f ), while the closedstate lifetime remains very similar (Fig. 5g). This indicates that the Y153A mutation gives the channel an increased stability at the open state. Taken together, the electrophysiological behaviors of the Y153A mutant and the structural observations both demonstrate that Tyr153 has a pivotal role in the gating process of the SsTRIC channel. It acts as a tethered plug being latched onto the surface of the pore lumen that occludes the pore from the cytoplasmic side.

Discussion
As shown in Additional file 5: Figure S5d and e, the site for TRIC channels to select and bind monovalent cations (K + /Na + ) during their permeation through the pore is most likely located at the central region around the two kinks on the M2 and M5 helices according to the following evidences. Firstly, two discrete water-/ionbinding sites (P0 and P0′) have been located in this region, and one of them (P0′) is accessible to a K + surrogate (Tl + ). Secondly, the coordination bond lengths between the ion on P0′ and the surrounding ligands are consistent with the K + -ligand bond length values observed in the KcsA or valinomycin structures (Additional file 5: Figure S5a-c). Thirdly, the kinks on M2 and M5 are general features conserved in the structures of both prokaryotic and eukaryotic TRIC channels. Such an asymmetric filter-like structure is evidently different from those of KcsA (or other homeo-tetrameric K + channels) and may represent a new type of monovalent cation-binding site for ion channels. In this filter, K + was bound in a monohydrated state and forms a trigonal geometry with its ligands, resembling half of the trigonal-antiprism K + coordinations in valinomycin (Additional file 5: Figure S5a and b). For comparison, the filter in KcsA contains four consecutive layers of discrete K + -binding sites with a square antiprism geometry, and the K + ion is bound in the center (Additional file 5: Figure S5f). Such a multilayer symmetrical filter in KcsA is more selective (for K + ) than the single-layer asymmetric filter in the TRIC channel. The P K /P Na permeability ratio for the KcsA channel is 166.7 or higher [32], whereas the P K /P Na ratio of the SsTRIC and mammalian TRIC channels is much lower at 1.21 and 1.5 [5], respectively. Lastly, the F104A mutant of SsTRIC has a lower P K /P Na ratio compared to the wild type, indicating that the mutant becomes less selective between K + and Na + . The result suggests that the geometry of the fourth pore constriction site (G4) shaped by the bulky side chain of Phe104 and M2/M5 kinks is likely crucial for tuning the relative selectivity of SsTRIC channel for K + versus Na + .
The plug regions in prokaryotic TRIC orthologs contain a consensus "PX 5-7 E(D/Q)XYA" motif, while the regions interacting with the plug comprise two consensus motifs, namely the "DA(T/S)XGL" motif on M4 and the "GGXXRD" motif on M5 (the bold residues are those involved in binding the plug). Such conservation suggests that the Velcro-like plug-pore interactions might serve as a general gating mechanism for the prokaryotic TRIC orthologs. These three consensus motifs have become diversified in the eukaryotic members. The comparison between SsTRIC and a CeTRIC-B channel reveals an unexpected diverse feature of the pore plugs among the prokaryotic and eukaryotic members of the family (Additional file 9: Figure S9a and b). While the prokaryotic orthologs utilize a built-in plug in the M5-M6 loop region to control their cytoplasmic gates, the eukaryotic members have evolved a different strategy by adopting a lipid molecule, namely phosphatidylinositol 4,5-bisphosphate (PIP 2 ), as the plug to their pores. The head group of PIP 2 latches onto two basic residues on the pore lumen surface and occludes the pore from the cytoplasmic side (Additional file 9: Figure S9b). The M5-M6 loop region in the eukaryotic TRIC proteins turns into an irregular structure capping over the PIP 2 head group, instead of forming a plug motif itself. In spite of these evident differences between prokaryotic and eukaryotic members of TRIC family, they share a common overall channel architecture and are probably unified by a general gating mechanism involving a Velcro-like plug-pore interaction. A model accounting for the potential gating mechanism of SsTRIC has been summarized in Additional file 10: Figure S10.
The previous bioinformatics study suggested potential roles of prokaryotic TRIC orthologs in the efflux of metabolites [23]. This prediction awaits verification by further experimental evidences to demonstrate the physiological relevance of TRIC channel function in the transport processes of metabolites. A recent work on RsTRIC and SsTRIC indicated that they are involved in K + uptake in E. coli cells [24]. Expression of SsTRIC on the E. coli membrane leads to moderate inhibition of cell growth, resembling the phenotype of a gain-of-function ion channel (Additional file 6: Figure S6b). It is likely constitutively active on the E. coli membrane when a resting potential at −220 to −140 mV is present at different growth phases [33]. Putatively, the prokaryotic TRIC members might serve as membrane-potential regulators to prevent imbalance of resting potential (by providing K + flux) and to maintain ion homeostasis in the cytoplasm of prokaryotic cells, instead of serving as a metabolite transporters themselves. Such a role may serve to regulate the active transport processes of metabolites across the membrane indirectly, as many transporters are driven by the electrochemical potential across the membrane [34].

Conclusions
In summary, the high-resolution view of ion-and water-binding sites in the SsTRIC channel unravels the presence of an asymmetric filter-like structure buried in the middle of an hourglass pore. The Velcro-like plug-pore interacting model derived in this study may offer a unified framework for understanding the gating mechanism of both prokaryotic and eukaryotic TRIC channels.

Cloning, protein expression, and purification
The gene encoding SsTRIC was synthesized (GenScript) with optimized codon usage for protein expression in Escherichia coli and inserted between the NdeI and XhoI sites in the pET21b vector. The C41 (DE3) E. coli strain was used for protein expression. For large-scale expression, the overnight starter culture was inoculated into Terrific Broth (TB) media at a 1:40 (v:v) ratio. The cells were grown to optical density OD 600 = 1.0 and then induced with 1 mM isopropyl β-D-1-thiogalactopyranoside (IPTG) at 37°C for 2 h before being harvested through centrifugation.
For the purification of SsTRIC protein, 10 g of defrosted cells were suspended in 100 ml of lysis buffer (50 mM Tris-HCl pH 8.0, 200 mM NaCl, 10 mM imidazole). After homogenization, 1.5 g of dodecyl-β-Dmaltoside (DDM) was added to the suspension to extract the membrane proteins, and the mixture was stirred on ice for 30 min. The preparation was sonicated for 2 min (1 s on and 5 s off ) in an ice-water bath. The cell lysate was centrifuged for 30 min at 16,000 RPM, 4°C in a JA-25.50 rotor (Beckman). The supernatant was collected and applied by gravity to a column with 2 ml Ni-NTA superflow resin (Qiagen). After all the sample had flowed through, the column was washed with 10 ml of equilibration buffer (25 mM Tris-HCl pH 7.5, 200 mM NaCl, 20 mM imidazole, 0.5% DM) and then 10 ml of washing buffer (25 mM Tris-HCl pH 7.5, 200 mM NaCl, 80 mM imidazole, 0.4% DM). After the elution buffer (25 mM Tris-HCl pH 7.5, 200 mM NaCl, 300 mM imidazole, 0.4% DM) was applied, the fractions with protein concentration above 0.3 mg ml -1 were pooled and concentrated to 10-15 mg ml -1 in a 50 kDa cutoff Amicon Ultra-4 centrifugal filter unit (Millipore). Next, trypsin was added to the concentrated SsTRIC protein sample at a 1:1000 (m:m) ratio and incubated at 37°C for 1 h for limited proteolysis treatment. The product was then applied to a Superdex 200 10/300 GL column (GE Healthcare) and eluted in a buffer containing 10 mM Tris HCl pH 7.5, 150 mM NaCl, and 1.1% OG. The major peak fraction of the eluted protein was pooled and concentrated to~10 mg ml -1 for the crystallization experiments.
Crystallization, crystal soaking, data collection, and processing The initial crystallization condition was identified through sparse matrix screening using the MemGold kit (Molecular Dimensions). The optimized recipe for crystallization of SsTRIC in OG involves setting up hanging drops with 1 μl 10-15 mg ml -1 protein sample, 1 μl well solution (22-26% PEG 3000, 0.1 M Tris-HCl pH 8.0, 0.2 M sodium acetate, and 0.2 M KCl), and 0.2 μl 30% ethylene glycol. The drops were equilibrated against 0.5 ml of well solution at 16°C. The best crystals were small hexagonal plates (0.1-0.3 mm in diagonal length) grown on the surface of the drops. To prepare crystal samples with only K + , the plate crystals initially grown in the mixture of NaAc and KCl were washed extensively in a stabilizing solution with KCl as the only salt (26% PEG 3000, 50 mM Tris-HCl pH 8.0, 1.1% OG, 0.1% DM, 3% ethylene glycol, 0.5 M KCl). These crystals were further soaked in the same solution at 16°C for~24 h before being harvested. The Tl + -containing crystals were prepared by washing in a stabilizing solution with 0.5 M KNO 3 to remove Clions. Otherwise, they would have formed a precipitate with Tl + ions because TlCl has low solubility in water. Next, the crystals were washed and soaked in a new stabilizing solution with 0.5 M TlNO 3 for 24 h.
For phasing, the crystals of A15C or A44C mutants were grown under the same conditions as for the wild type and derivatized with 1 mM CH 3 HgCl for 24-48 h. The mutation was introduced through the QuikChange Site-Directed Mutagenesis method (Stratagene). The high-resolution KCl dataset (Table 2) was collected at the BL17U of the Shanghai Synchrotron Radiation Facility (SSRF) equipped with an ADSC Q315 CCD detector. The anomalous diffraction data of the Tl + -containing crystals were collected at the home source with a Rigaku Micromax-007 rotating anode generator and a Rigaku IP IV++ detector. For data processing, the iMOSFLM [35] or HKL2000 [36] program was used.

Structure determination, refinement, and analysis
The structure was solved by the single-wavelength anomalous dispersion method using the Autosol program in the PHENIX suite [37]. The anomalous diffraction dataset of Hg-labeled A15C crystal collected at 1.00000 Å wavelength was used for phasing. Six Hg atoms were located in one asymmetric unit, and the initial set of selected area diffraction (SAD) phases has a figure of merit (FOM) of 0.258 which is further improved to 0.640 through density modification. The initial model automatically generated by the Autobuild program in PHENIX contains a nearly complete polypeptide structural model with R work = 0.250 and R free = 0.260. This model was further improved by iterative cycles of manual adjustment in Coot [38] and refinement in CNS 1.2 [39] using the high-resolution (2.2 Å) KCl-only data. The longest polypeptide chains were continuously traced from Met 1 to Pro 198. From the electron density map, it appears that the limited proteolysis by trypsin removed only seven residues at the carboxy-terminal region after Pro198 and the hexahistidine tag region. The transmembrane domain and the loop region remain intact after proteolysis. The K + ions were modeled and cross-validated by the anomalous difference Fourier peaks of the Tl + ion. The anomalous difference Fourier map of Tl + -containing crystals was computed by the FFT program in the CCP4 suite [40]. The data collection, phasing, and structure statistics are summarized in Table 2.
For the structural analysis, the PROMOTIF program [41] was used to analyze the secondary structures, Lsqman [42] was used to superpose different structures, HOLE [43] was used to probe the pores within the channel, PISA [44] was used to analyze protein interfaces and buried surfaces within the trimer, and the APBS [45] tool was used to calculate surface electrostatic potential. For sequence alignment, the output from the ClustalW [46] program was manually checked The R sym , R pim , I/σ, completeness, and redundancy are presented as the statistics for overall and outmost shell (in parentheses) a SSRF Shanghai Synchrotron Radiation Facility b R pim is the precision-indicating (multiplicity-weighted) R sym reported by the Scala program.
The HKL2000 program used to process TlNO 3 data did not report R pim statistics c Percentage of residues in most favored, additional allowed, and generously allowed regions in Ramachandran plot. No residues are observed in the disallowed region and readjusted. BOXSHADE 3.21 and ESPript programs [47] were used to generate sequence alignment figures. The cartoon structural figures were produced in PyMOL [48] or Chimera [49], and the electron density maps were displayed using Coot [38].

Lipid extraction and identification
To extract lipids, purified SsTRIC or E. coli TRIC (EcTRIC) protein samples (200 μl 10 mg ml -1 protein in β-DM) were mixed with 180 μl of solvent with chloroform, methanol, and concentrated HCl solution (1:2:0.02, v/v/v). Subsequently, 60 μl of chloroform and 60 μl of 2 M KCl (sigma) were added to each tube. The mixture was vortexed and then centrifuged for 5 min at 2400 g to separate the organic phase from the aqueous phase. The organic phase was then separated through thin layer chromatography (TLC). At first, 20 ul organic phases extracted from SsTRIC and EcTRIC protein samples were spotted on the TLC silica gel plate (Merck). As reference standards, triacylglycerol (TG, C18:1), diacylglycerol (DG, C16:0), monoacylglycerol (MG, C16:1), palmitic acid, and steric acid samples were applied on the same plate. The solvent used for TLC contains nhexane, diethyl ether, and acetic acid (70:30:1, v:v:v). After being separated on the TLC plate, the lipid fractions were visualized by spraying the plate with 0.5% iodinechloroform. To further identify the lipid fraction, mass spectroscopy was performed under the electrospray ionization (ESI) positive and negative scan modes. As a control, the organic phase extraction of the blank elution buffer used for protein purification was also examined through mass spectroscopy under the same modes.

Reconstitution of SsTRIC in small unilamellar vesicles (SUVs)
A lipid mixture contraining 90% 1,2-diphytanoyl-sn-glycero-3-phosphocholine (DPhPC) and 10% cholesterol (w/w) in chloroform was dried under vacuum in a CentriVap Concentrator (Labconco) for 4 h. The lipid sample was then suspended at 10 mg ml -1 in a low-salt buffer (1 mM 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid (HEPES) pH 7.2, 5 mM KCl). SUVs were formed by tip sonication (50 Hz, 1 s on, 1 s off for 1 min) and then presolubilized by 10 mM DM for 30 min at RT. Subsequently, the SsTRIC protein sample was added to the presolubilized SUVs to achieve a protein:lipid ratio of 1:10 (m:m). More DM was added to a final concentration of 17.5 mM, and the resulting mixture was gently agitated for 1 h at room temperature. Detergent was removed by dialysis in the low-salt buffer (1 mM HEPES pH 7.2, 5 mM KCl). The external buffer was changed every 12 h for 3 days. After dialysis, the resulting SUVs were aliquoted, flash-frozen in liquid nitrogen, and stored at −80°C.

Preparation of giant unilamellar vesicles (GUVs) for electrophysiology
The GUV samples were generated by the electroformation technique using the Nanion Vesicle Prep Pro device (Nanion). Before electroformation, trehalose was added to the preformed SUV solution to a final concentration of 10 mM in order to protect the SsTRIC protein during the partial dehydration process. About 10 μl of SUV solution was applied in small droplets (~0.2 μl/droplet) on the indium tin oxide (ITO)-treated glass slide. The droplets were left to dry in the open air for approximately 30 min at room temperature. Subsequently, 300 μl of 1 M sorbitol solution was carefully added onto the lipid film, and a cassette sandwiching the sample in the middle was assembled. During assembly, we ensured that the ITO layers of the slides were facing and that they touched the sample. For the electroformation process, the protocol was set as 0.1 to 1.0 V at 12 Hz frequency for 3 h. For the next step, the frequency was lowered to 4 Hz and the voltage was raised to 2 V for 30 min to detach the GUVs from the glass slides. The temperature was kept constant at 36°C throughout the electroformation process. Alternatively, the purified protein can be reconstituted directly on the GUVs (protein:lipid = 1:200-250, w:w; lipid: 95% azolectin + 5% cholesterol, w:w) prepared through a modified sucrose method [50]. In some cases when the GUVs did not attach well to the bottom of the sample chamber, the lipid composition was adjusted by adding DPhPC to the mixture at a cholesterol:DPhPC:azolectin (w:w:w) ratio of 1:5:17 or 1:10:15.

Electrophysiology
All recordings were performed with the inside-out configuration. The intracellular side of the GUV membrane was exposed to the bath solution, and the extracellular side was exposed to the pipette solution. The bath and pipette solutions contained 210 mM KCl and 10 mM HEPES (pH 7.2). Patch pipettes with resistances of 8-9 MΩ were used, and the patch resistance increased to~2 GΩ after the pipette sealed tightly with the GUV membrane. For recording the data under bi-ionic conditions, the pipette solution was kept constant with 210 mM KCl and 10 mM HEPES (pH 7.2), while the bath solutions were: 210 mM NaCl and 10 mM HEPES (pH 7.2) for P K /P Na analysis; 210 mM KCl, 75 mM CaCl 2 or MgCl 2 , and 10 mM HEPES (pH 7.2) for the test of Ca 2+ /Mg 2+ and Clpermeability. Single-channel recordings were made with an EPC-10 amplifier (HEKA, Lambrecht, Germany) under different voltage settings at 50 kHz with a 0.5-kHz filter and a 50-Hz notch filter. All experiments were done at room temperature (21-24°C).
The Clampfit Version 9.0 (Axon Instruments, Foster City, CA) was used for data analysis, Excel Version 2010 (Microsoft) and OriginPro 8 were used for statistical analysis, and Igor Pro 6.32A (WaveMetrics, USA) was used for graphics. The single-channel conductance was obtained through linear fitting of the currentvoltage plots. The statistical data are reported as mean value ± SEM. The Student's t test was used to assess statistical significance; n represented the number of experiments analyzed.

Cell-based functional assay
The plasmids carrying SsTRIC mutants were transformed into the C41 (DE3) E. coli strain for protein expression and in vivo functional assay. For the measurement of cell growth curves, single colonies of the E. coli transformants were used to inoculate 3 ml of liquid LB media plus 100 μg/ml ampicillin and grown at 37°C for 2-3 h. The cell densities of the cultures were normalized to OD 600 = 0.1 with LB media. Subsequently, 75 μl of diluted cell culture was added into the wells of a 96-well plate and mixed with 75 μl LB media with 1 mM IPTG added. The plate was then sealed with the CyclerSeal film (Platemax), and the growth of cells at 37°C within each well was monitored continuously in the Thermo Varioskan Flash Plate Reader. The OD 600 was taken every 20 min while the plate was shaken at 300 cycles of horizontal shakes per minute (SPM).
The protein expression was analyzed through western blots on the membrane fractions of the E. coli cells harvested after being induced by IPTG. To prepare the membrane fractions, 1 g of defrosted cells was suspended in 5 ml lysis buffer (50 mM Tris-HCl pH 8.0, 200 mM NaCl, 10 mM imidazole). Lysozyme was added to a final concentration of 1 mg ml -1 , and the suspension was stirred at 4°C for 1 h. The preparation was sonicated for 2 min in an ice-water bath with a program setting of 1 s on and 5 s off. The cell lysate was spun at 14,800 RPM, 4°C for 30 min. The supernatant was collected and centrifuged at 14,800 RPM, 4°C for another 30 min to remove large insoluble debris. Subsequently, the supernatant (1 ml) was ultracentrifuged at 38,000 RPM (~100,000 × g) in an S140AT rotor (Hitachi) at 4°C for 2 h. The pellet was resuspended in 1 ml of washing buffer (50 mM Tris-HCl pH 8.0, 200 mM NaCl, 10 mM imidazole) and ultracentrifuged again. The membrane pellets were collected and weighted, followed by resuspension in the washing buffer to a final concentration of 50 mg ml -1 . The samples were mixed with 5 × SDS loading buffer and then separated by electrophoresis on 12% SDS-PAGE gel. For western blot analysis, the protein bands on the gel were transferred to a nitrocellulose membrane. After blocking and washing, the blot was probed by HRP-conjugated anti-Histag mAb antibody (GenScript, cat. number: A00612, lot number: 16B001004, RRID: AB_915573) and then developed with the Western Lightning ULTRA substrate (PerkinElmer).

Additional files
Additional file 1: Figure S1. Superposition of SsTRIC structure with two previous structures of prokaryotic TRIC orthologs. a and b SsTRIC trimer and monomer superposed on the structure of SsTRIC-Fab complex (Protein Data Bank (PDB):5H35). The Fab antibody fragment is omitted for clarity. c and d SsTRIC trimer and monomer superposed on the structure of RsTRIC (PDB:5H36). Protein backbones are represented as ribbon cartoon models. Additional file 3: Figure S3. The pore structures of SsTRIC compared to those of RsTRIC and CeTRIC-B1. a SsTRIC, b RsTRIC, c CeTRIC-B1. PDB codes: RsTRIC, 5H36; CeTRIC-B1, 5EGI. Color codes for the pore profiles based on the data output by HOLE program: red, low radius surface (<1.5 Å); green areas, normal pore surface (1.5-2.3 Å); light blue areas, high radius surface (>2.3 Å). The protein backbones are shown as ribbon models. The key amino acid residues shaping the constriction areas along the pore are highlighted as green stick models, while the PIP 2 in CeTRIC-B1 is shown as yellow and red sticks. (TIF 891 kb) Additional file 4: Figure S4. Electron density maps of the pore region in SsTRIC channel. a 2F o -F c electron densities of the pore region contoured at 1.5 × σ level. The refined structures are shown as stick models superposed on the map. The water molecules and K + ions are displayed as red and silver bullets, respectively. b The same map zoomed in around the pore center to show the local view of the P0 and P0′ sites encircled by the M2 kink, M5 kink, Phe16, and Phe104. c The anomalous difference Fourier peaks of the Tl + ions bound to the SsTRIC channel. The map at 3.4-Å resolution is contoured at 3 × σ level and shown as magenta meshes superposed on the C α -trace model of a SsTRIC trimer. The images shown are stereo pairs, and the view is approximately perpendicular to the C3 axis. The Tl-1 and Tl-2 peaks are located at the sites corresponding to P0′ and P3′ sites shown in a. The relatively weaker densities of Tl-1 and Tl-2 peaks compared to those on surfaces indicate that the occupancy of Tl + ions on these internal sites is much lower than on the surface ones. Note that only one monomer in the trimer contains a peak above 3 × σ level at P0′ site (Tl-1), suggesting that this site at the current conformational state is not yet fully accessible to Tl + /K + . Stereo pairs are shown in a-c. (TIF 7223 kb) Additional file 5: Figure S5. Central putative monovalent cation-binding site and filter-like structure in SsTRIC compared to those of valinomycin and KcsA channel. a The ligands around the putative K + (purple sphere, likely with partial occupancy and mixed with water molecules) at pore center of SsTRIC. Zoom-in view shows the trigonal coordination mode of K + ; the side chain of Phe104 is located below K + . Dashed lines indicate the close interactions between K + and its ligands with the bond length (Å) labeled nearby. The interaction between K + and the backbone amide is not uncommon, as some similar interactions are found in the structure of gramicidin (PDB code 2IZQ, the K + -NH bond length is 3.1-3.4 Å, slightly weaker than the one we observed.). b The trigonal antiprism-type coordination of K + in valinomycin (Cambridge Crystallographic Data Centre/CCDC accession code: VALINK). DVA and V indicate D-and L-valines, respectively. Gray dashed lines are the coordination bonds absent in SsTRIC. c The square antiprism-type coordination of K + in KcsA. (PDB:1K4C). In panels a-c, K + ions and water molecules are shown as purple and red spheres, respectively. The amino acid residues coordinating K + ions are highlighted as ball-and-stick models. d and e The single layer filter-like structure sandwiched among M2 kink, M5 kink, Phe104, and Phe16 in SsTRIC. The putative K + ion and a nearby water molecule are shown as sphere models. Amino acid residues (Phe16, Thr41, Gly45, Phe104, Val129, and Gly133) involved in binding K + (or Na + ) and water binding are shown as stick models. View in e is rotated~90°with respect to view in a. f The canonical four-layer filter structure in KcsA channel. Purple funnels in e and f indicate the presumed permeation pathways for monovalent cations. (TIF 1214 kb) Additional file 6: Figure S6. In vivo assay on the function of SsTRIC and various mutants. a Side view of the SsTRIC monomer with the residues chosen for alanine-scanning mutagenesis highlighted in cyan or yellow. b Growth curves of E. coli cells expressing wild-type (WT) and various mutants of SsTRIC. c Control data for the cell growth assays. EcMscL: the large-conductance mechanosensitive channel (MscL) from E. coli. Overexpression of wild-type EcMscL does not inhibit the cell growth, as it remains tightly closed at the resting state, while its gain-of-function mutant G26H has a dramatic inhibitory effect on cell growth. They are included as controls for comparison with the phenotypes of SsTRIC and its mutants. Vector: empty pET21b vector used for expressing SsTRIC and mutants, included as a negative control; the R187A mutation site is relatively distant from the gating region and thus serves as an internal control. The error bars represent the standard errors of mean values (SEM, n = 3). Additional file 8: Figure S8. Asymmetric features of the internal repeats in SsTRIC monomer. a Structure of an SsTRIC monomer represented as a cartoon ribbon model. Convex quadrilaterals in light blue and pink highlight the regions of motifs A and B, respectively. The dark solid ellipse at the center shows the position of a pseudo-C2 axis relating the two motifs. The view is along the pseudo-C2 axis running nearly parallel to the membrane plane. b Superposition of motif B (golden) with motif A (blue). c Alignment of the amino acid sequences of motifs A and B. Identical residues are in dark boxes; similar ones are in gray backgrounds. The blue and golden arrows highlight the key residues around the pore region. d The intracellular vestibule in motif B is occluded. e The extracellular vestibule of the pore harbored in motif A is open. The views are approximately along the pore axis. The key amino acid residues surrounding the pore region are shown as stick models. The asterisk symbols indicate the approximate location of the pore. (TIF 3075 kb) Additional file 9: Figure S9. Comparison between SsTRIC and CeTRIC-B2 structures. a Structures of SsTRIC and CeTRIC-B2 monomers superposed and viewed from the same angle. The SsTRIC (left) and CeTRIC-B2 (right, PDB code: 5EIK) backbones are represented as cartoons, and the lipid cofactor (PIP 2 ) and key amino acid residues are shown as sticks. The plug motif/ molecule is highlighted in cyan, while those involved in binding the plug are colored in yellow. b The position of plug motif in SsTRIC (left) compared to that of PIP 2 in CeTRIC-B2 (right). The key residues involved in binding the plug motif/PIP 2 are represented as yellow sticks. (TIF 2656 kb) Additional file 10: Figure S10. Proposed model accounting for the gating mechanism of SsTRIC channel. At the closed state, the plug motif forms tight interactions with the pore lumen surface and blocks the cytoplasmic entrance. Upon activation by voltage, movement of the charged residues under electric field leads to destabilization of the interactions between plug motif and the potential voltage-sensing residues. The plug is thereby dislocated from the blocking position so that the cytoplasmic entrance of the hourglass-shaped pore is open for K + ions to permeate through the pore. Blue and red sticks represent the positively and negatively charged residues, respectively; purple spheres represent K + ions. For clarity, M3 and M7 helices are omitted in the model. For K v or Na v channels, they generally adopt a string of positively charged amino acid residues (Arg and Lys) on a transmembrane helix (S4 helix) of the voltage-sensing domain (VSD) to sense electrical signals [51][52][53]. While SsTRIC channel does not contain a canonical VSD in its structure, it does have nine positively charged residues (Arg and Lys) distributed asymmetrically on two solvent-exposed surfaces (Fig. 4b), reflecting the positive-inside rule for membrane proteins [54]. The eight positively charged residues on the luminal side are not randomly distributed. Among them, six cluster around the cytoplasmic gate area, and three of them (Lys 27, Lys 23 on M1 helix and Arg137 on M5 helix) form an array of positively charged regions lining the interfacial groove between M1 and M5 (Fig. 4b). This array of positively charged residues is stabilized by three acidic residues: Asp97, Asp138, and Glu143 (Fig. 4c). Furthermore, three key residues from the plug motif interact directly or indirectly with these potential voltage-sensing residues (Fig. 4c). Hence, the Velcro-like structure may serve as a basis for the potential voltage-dependent regulation of SsTRIC channel activity. (TIF 743 kb) Additional file 11: Supporting data for Fig. 2 and 5e-g. (XLSX 19 kb)