Integral membrane glutamate transporters couple the concentrative substrate transport to ion gradients. There is a wealth of structural and mechanistic information about this protein family. Recent studies of an archaeal homologue, GltPh, revealed transport rate heterogeneity, which is inconsistent with simple kinetic models; however, its structural and mechanistic determinants remain undefined. Here, we demonstrate that in a mutant GltPh, which exclusively populates the outward-facing state, at least two substates coexist in slow equilibrium, binding the substrate with different apparent affinities. Wild type GltPh shows similar binding properties, and modulation of the substate equilibrium correlates with transport rates. The low-affinity substate of the mutant is transient following substrate binding. Consistently, cryo-EM on samples frozen within seconds after substrate addition reveals the presence of structural classes with perturbed helical packing of the extracellular half of the transport domain in regions adjacent to the binding site. By contrast, an equilibrated structure does not show such classes. The structure at 2.2-Å resolution details a pattern of waters in the intracellular half of the domain and resolves classes with subtle differences in the substrate-binding site. We hypothesize that the rigid cytoplasmic half of the domain mediates substrate and ion recognition and coupling, whereas the extracellular labile half sets the affinity and dynamic properties.
Introduction
Membrane glutamate transporters pump substrates against their concentration gradients, serving critical biological functions across all kingdoms of life. In mammals, excitatory amino acid transporters recycle glutamate from the synaptic cleft into the glia (Freidman et al., 2020). In prokaryotes, orthologs take up nutrients, including glutamate, aspartate, neutral amino acids, or dicarboxylic acids (Kim et al., 2002; Burguiere et al., 2004; Youn et al., 2009). Transporters utilize energy from downhill ionic electrochemical gradients to carry concentrative substrate uptake. Excitatory amino acid transporters rely on inward Na+ and proton gradients and an outward K+ gradient (Zerangue and Kavanaugh, 1996). Prokaryotes couple transport to either proton or Na+ gradients (Tolner et al., 1995; Ryan et al., 2009).
These transporters are homotrimers with each protomer composed of a rigid scaffold trimerization domain and a mobile transport domain containing the ligand-binding sites. Protomers functions independently (Erkens et al., 2013; Ruan et al., 2017; Riederer et al., 2018; Georgieva et al., 2013; Koch et al., 2007; Grewer et al., 2005; Koch and Larsson, 2005). Transport domains translocate ligands across membranes by moving ∼15 Å from the outward-facing state (OFS) to the inward-facing state (IFS), in an elevator motion (Reyes et al, 2009; Garaeva et al, 2019; Arkhipova et al, 2020; Qiu et al, 2021). Studies in archaeal Na+-coupled transporters GltPh and GltTk led to a simple kinetic model of transport (Boudker et al., 2007; Reyes et al., 2013; Verdon et al., 2014; Oh and Boudker, 2018; Guskov et al., 2016; Riederer and Valiyaveetil, 2019; Wang and Boudker, 2020; Arkhipova et al., 2020; Alleva et al., 2020). Briefly, in the OFS, ion binding to Na1 and Na3 sites reveals the substrate and an additional sodium (Na2)-binding site through an opening of helical hairpin 2 (HP2), also preventing the translocation of Na+-only bound transport domain. Subsequent binding of the substrate and Na2 closes HP2, allowing translocation to the IFS and ligand release into the cytoplasm.
Recently, high-speed atomic force microscopy, single-molecule FRET (smFRET) TIRF microscopy, and 19F-NMR revealed a more complex picture of GltPh transport and dynamics (Huang et al., 2020; Ciftci et al., 2020; Akyuz et al., 2015; Matin et al., 2020; Erkens et al., 2013; Akyuz et al., 2013; Huysmans et al., 2021; Ciftci et al., 2021). These studies established the existence of additional conformational substates in OFS and IFS, of which some translocate and transport at different rates. Although it is expected that cryo-EM studies would resolve these proposed conformational substates from heterogeneous datasets, this so far does not appear to be the case (Wang and Boudker, 2020; Arkhipova et al., 2020). Therefore, the structural and mechanistic determinants of kinetic heterogeneity remain unclear.
Using a GltPh mutant that exclusively occupies the OFS, we show that substrate binding in the OFS is heterogeneous, consistent with at least two substates with different affinities. Salt composition and temperature modulate the substate populations, suggesting that they are in equilibrium. However, the conformational exchange must be slow, in the order of at least tens of seconds, for the substates to manifest in binding isotherms. Notably, a similar conformational equilibrium also exists in WT protein and potentially contributes to heterogeneous transport kinetics. The substrate-bound high-affinity state of the mutant transporter is expected to predominate after equilibration. Thus, to gain insights into the structure of the transient low-affinity substate, we conducted an extensive analysis of the cryo-EM imaging data collected on the transporter frozen within seconds after substrate addition. The identified structural classes reveal subtle differences. Specifically, we observed a subset of structural classes with differently packed helices in the extracellular half of the transport domain adjacent to the binding site, suggesting that the region is labile. We hypothesize that the ensemble of transient, more dynamic, less uniquely packed conformations comprises the low-affinity substate. Images of an equilibrated protein produced reconstructions at a uniquely high 2.2 Å resolution, revealing a complement of structured waters in the cytoplasmic side of the transport domain that may contribute to its conformational rigidity. We did not observe a conformational heterogeneity of the extracellular half of the transport domain in these data, which we attribute to the relaxation of the protein to the higher affinity state. Classifications instead revealed subtle differences in the substrate-binding site and the global orientations of the transport domains, which could also contribute to kinetic heterogeneity. Our results provide a framework in which modal kinetic behavior demonstrated by GltPh may be a result of subtle but long-lived structural differences.
Materials and methods
DNA manipulations and protein preparation
Mutations were introduced to the previously described GltPh CAT7 construct (Yernool et al., 2004) using PfuUltra II, and sequences were verified using Sanger sequencing (Macrogen). Proteins were expressed as C-terminal (His)8 fusions, separated by a thrombin cleavage site. Proteins were purified as previously described (Yernool et al., 2004). Briefly, crude membranes of DH10B Escherichia coli cells overexpressing GltPh were solubilized for 2 h in 20 mM HEPES/NaOH, pH 7.4, 200 mM NaCl, 5 mM monopotassium L-aspartate (L-Asp), and 40 mM n-dodecyl-β-D-maltopyranoside (DDM; Anatrace). After solubilization, the DDM was diluted to ∼8–10 mM, and after high-speed ultracentrifugation (40,000 rpm, Ti45 rotor), the supernatant was applied to pre-equilibrated Ni-NTA affinity resin (Qiagen) for 1 h. The resin was washed with seven column volumes of 20 mM HEPES/NaOH, pH 7.4, 200 mM NaCl, 5 mM L-Asp, and 40 mM imidazole. Subsequently, the resin was eluted with the same buffer with increased imidazole (250 mM). The (His)8 tag was cleaved by thrombin digestion overnight at 4°C, and the proteins were further purified by size-exclusion chromatography (SEC) in the appropriate buffer for subsequent experiments. The concentration of GltPh protomers was determined in a UV cuvette with a 10-mm pathlength (Starna Cells), using protein diluted 1:40, and an experimentally determined extinction coefficient of 57,400 M−1 cm−1 (Reyes et al., 2013).
Reconstitution and uptake assays
Liposomes able to maintain proton gradients were prepared using a 3:1 mixture of POPE and POPG (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine and 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-[1′-rac-glycerol]). Lipids were dried on the rotary evaporator for 2 h and under vacuum overnight. The resulting lipid film was hydrated by 10 freeze–thaw cycles at a concentration of 5 mg/ml in 50 mM potassium phosphate buffer and 100 mM potassium acetate, pH 7. The suspensions were extruded using a Mini-Extruder (Avanti) through 400-nm membranes (Whatman) 10 times to form unilamellar liposomes, and Triton X-100 was added to liposomes at a ratio of 1:2 (wt/wt).
P-GltPh for reconstitution was affinity purified, thrombin cleaved, and purified by SEC in 20 mM HEPES/Tris, pH 7.4, 200 mM NaCl, 1 mM L-Asp, and 7 mM DDM. Purified protein was added to destabilized liposomes at a ratio of 1:1,000 (wt/wt) and incubated for 30 min at 23°C. Detergent was removed with four rounds of incubation with SM-2 beads (Bio-Rad) at 80 mg beads per 1 ml of liposome suspension (2 h at 23°C twice, overnight at 23°C once, and 2 h at 23°C once). Before use, SM-2 beads were washed in methanol, rinsed thoroughly with distilled water, and equilibrated in the liposome internal buffer. After detergent removal, proteoliposomes were concentrated to 50 mg/ml by ultracentrifugation at 86,000 g for 40 min at 4°C, freeze–thawed three times, and extruded through 400-nm membranes 10 times.
Uptakes were initiated by diluting reconstituted proteoliposomes 1:100 in the appropriate reaction buffer preincubated at 30°C. At the indicated time points, 200-μl reaction aliquots were removed and diluted in 2 ml of ice-cold quenching buffer (20 mM HEPES/Tris, pH 7, and 200 mM LiCl). The quenched reaction was immediately filtered through a 0.22-µm filter membrane (Millipore Sigma) and washed three times with 2 ml quenching buffer. Washed membranes were inserted into scintillation vials, and the membranes were soaked overnight in 5 ml Econo-Safe Counting Cocktail. Radioactivity in liposomes was measured using an LS6500 scintillation counter (Beckman Coulter).
Isothermal titration calorimetry (ITC)
For the TFB-TBOA competition experiments, protein was first titrated with TBOA to saturation. Then, the concentration of TFB-TBOA was increased to a final concentration of 150 μM of TFB-TBOA, so that the overflow protein was also saturated. The appropriate amount of L-Asp was subsequently titrated into the saturated protein to yield a binding curve.
Cryo-EM imaging data collection
In both data sets, 3.5 μl protein at ∼4.5 mg/ml was applied to a glow-discharged QuantiFoil R 1.2/1.3, 300-mesh, gold grid (Electron Microscopy Sciences). Grids were blotted at room temperature and 100% humidity for 3 s at 0 blot force and plunge frozen in liquid ethane using a VitroBot Mark IV (FEI). Data S1 was collected on P-GltPh, which was affinity-purified, concentrated, and buffer-exchanged into 20 mM HEPES/Tris, pH 7.4, 99 mM K-gluconate, 1 mM Na-gluconate, and 0.8 mM DDM to remove the substrate. The protein was then SEC purified in 20 mM HEPES/Tris pH 7.4, 250 mM NaNO3, and 0.8 mM DDM. 2.5 μl of substrate-free protein was applied to the grid, then 1 μl of L-Asp was added and mixed just prior to freezing so that the final substrate concentration was 1 mM, the final protein concentration was ∼4.5 mg/ml, and protein was exposed to the substrate for ∼5 s prior to freezing (including blot time). Data S2 was collected on P-GltPh SEC purified in 20 mM HEPES/Tris, pH 7.4, 250 mM NaNO3, 1 mM L-Asp, and 0.8 mM DDM. The final buffer conditions of Data S2 were identical to Data S1.
All imaging data were collected on Titan Krios microscopes (FEI) operated at 300 kV. Data S1 was collected on a K2 Summit direct electron detector (Gatan). Automated data collection was performed in counting mode using Leginon software (Suloway et al., 2005), with a magnification of 22,500×, electron exposure of 70.23 e−/Å2, 50 frames/s, a defocus range of −1.0 to −2.0 μm, and pixel size of 1.073 Å. Data S2 was collected with a K3 Summit direct electron detector (Gatan). Automated data collection was performed in superresolution counting mode using SerialEM software (Mastronarde, 2005) with a magnification of 81,000×, electron exposure of 47.91 e−/Å2, 30 frames/s, defocus range of −0.5 to −2.5 μm, and pixel size of 0.53 Å.
Image processing
The frame stacks were motion-corrected using MotionCor2 (Zheng et al., 2017), with 2× binning in the case of Data S1, and contrast transfer function (CTF) estimation was performed using CTFFIND 4.1 (Rohou and Grigorieff, 2015). All datasets were processed using cryoSPARC 3.0 and Relion 3.0.8 simultaneously with default parameters unless otherwise stated (Su et al., 2020; Punjani et al., 2017; Zivanov et al., 2018). Specific information on the processing of each dataset is in Figs. S5 and S7. In brief, particles were nonspecifically picked from micrographs using the Laplacian of Gaussian (LoG) picker, aiming for ∼2,000 picks per micrograph. These particles were extracted at a box size of 240 pixels with 4× binning and then imported to cryoSPARC. The particles underwent one round of 2-D classification to remove artifacts, and then multiple rounds of heterogeneous refinement (C1) using eight total classes, seven of which were noisy volumes (created by one iteration of ab initio) and one was an unmasked 3-D model obtained from a previous processing pipeline. Once >95% of particles converged on a single class, the particles were converted back to Relion format via PyEM (Asarnow et al, 2019) and re-extracted at full box size. These particles were reimported to cryoSPARC, underwent three more rounds of heterogeneous refinement, then nonuniform (NU) refinement using C3 symmetry (dynamic mask threshold, dynamic mask near, and dynamic mask far were always set to 0.8, 12, and 28, respectively; Punjani et al., 2020). These particles were converted back to Relion format and underwent Bayesian polishing, using parameters obtained using 5,000 random particles within the set. After three more rounds of heterogeneous refinement in cryoSPARC and one round of NU-refinement, we performed local CTF refinement (minimum fit res 7 Å). After three more rounds of heterogeneous refinement and one round of NU-refinement, we performed global CTF refinement (three iterations, minimum fit res 7 Å; fitting trefoil, spherical aberration, and tetrafoil), three more rounds of heterogeneous refinement, and one round of NU-refinement.
Data S2 was also further processed to obtain the high-resolution structure by three additional rounds of polishing, local CTF refinement, and global CTF refinement as described above. During the polishing rounds, the box size and pixel size were rescaled as indicated in the supplement. After the second round of polishing, we classified single protomers by employing the relion_particle_symmetry_expand function in Relion to artificially expand the particle set three times (C3) so that each protomer rotated to the same position (Scheres, 2016). The expanded particle set was subjected to 3-D classification without alignment with T = 400 and 10 classes, using the refined C3 structure as a reference map. The exceptionally high T value was chosen to separate out subtle structural differences, and lower T values were also tested during processing. The local mask was created using a 20 Å map of the transport domain of Chain A of PDB accession no. 2NWX, with an initial binarization threshold of 0.01, extended by 3 pixels, and a soft-edge of 10 pixels. Of these, particle stacks from subsets of interest were separately used in cryoSPARC’s local refinement (C1), using the mask and map obtained from the most recent NU refinement. Single protomers in Data S1 were classified similarly, except we performed symmetry expansion after the first round of polishing due to the limited resolution of the dataset. After processing, the resulting half-maps were used as inputs for density modification implemented in PHENIX 1.19.1–4122 (Terwilliger et al., 2020; Adams et al., 2010), using a mask created from the NU-refinement job (threshold 0.1, dilation radius 15, soft padding width 5). All density maps were displayed using ChimeraX (Pettersen et al., 2021).
Model building and refinement
For atomic model building, the crystal structure of WT GltPh in the OFS (PDB accession no. 2NWX) was docked into the densities using UCSF Chimera (Pettersen et al., 2004). The model was first real-space refined in PHENIX (Adams et al., 2010). Then, chain A was adjusted manually, and ions were added in COOT (Emsley and Cowtan, 2004). Waters were initially added using phenix.douse, and subsequently manually inspected and adjusted. The resulting model underwent additional rounds of real-space refinement and validated using MolProbity (Chen et al., 2010). All structural models were displayed using ChimeraX (Pettersen et al., 2021). Per-residue Cα RMSDs were generated with Chimera (Meng et al., 2006). To cross-validate models, refined models (FSCsum) were randomly displaced an average of 0.3 Å using phenix.pdbtools. The displaced model was real-space refined against half-map 1 obtained through density modification to obtain FSCwork. The resulting model was validated against half-map 2 to obtain FSCfree.
Single-molecule dynamics assay
Thrombin-cleaved P-GltPh, containing C321A and N378C mutations, was labeled and reconstituted as described previously (Akyuz et al., 2013). Protein was SEC purified in 20 mM HEPES/Tris, 200 mM NaCl, 1 mM L-Asp, and 1 mM DDM. Purified protein was labeled using maleimide-activated LD555P-MAL and LD655-MAL dyes and biotin-PEG11 at a molar ratio of 4:5:10:2.5. Excess dyes were removed on a PDMiniTrap Sephadex G-25 desalting column (GE Healthcare).
All experiments were performed on a previously described home-built prism-based TIRF microscope constructed around a Nikon Eclipse Ti inverted microscope body (Juette et al., 2016). Microfluidic imaging chambers were passivated with biotin-PEG, as previously described (Akyuz et al., 2015). After passivation, the microfluidic channel was incubated with 0.8 μM streptavidin (Invitrogen) in T50 buffer (50 mM NaCl and 10 mM Tris, pH 7.5) for 7 min, then thoroughly rinsed with T50 buffer. Detergent-solubilized protein was immobilized by slowly flowing over the channel, and excess protein was removed by washing with 1 ml of 25 mM HEPES/Tris, pH 7.4, 200 mM KCl, and 1 mM DDM.
Buffers were supplemented with an oxygen-scavenging system composed of 2 mM protocatechuic acid and 50 nM protocatechuate-3,4-dioxygenase, as described previously (Aitken et al., 2008). The smFRET movies were recorded with 100-ms integration time using 80–100 mW laser power. All conditions tested were in the presence of 25 mM HEPES/Tris, 500 mM sodium salt, and 1 mM DDM, in the presence or absence of 1 mM L-Asp. Slides were washed with 25 mM HEPES/Tris, 200 mM KCl, and 1 mM DDM between experiments. Traces were analyzed using Spartan software (Juette et al., 2016). Trajectories were corrected for spectral cross talk and preprocessed automatically to exclude trajectories that lasted ≤15 frames and had a signal-to-noise ratio of ≤8. Traces with multiple photobleaching events (indicative of multiple sensors in a protein) or inconsistent total fluorescence intensity were also discarded.
Single-molecule transport assay and analysis
Thrombin-cleaved GltPh (C321A/N378C) was SEC purified in 20 mM HEPES/Tris, pH 7.4, 200 mM NaCl, and 0.1 mM L-Asp. Protein was labeled with maleimide-activated biotin-PEG11 (EZ-Link; Thermo Fisher Scientific) in the presence of N-ethylmaleimide at a molar ratio of 1:2:4 as previously described (Ciftci et al., 2020). Liposomes were prepared from a 3:1 (wt/wt) mixture of E. coli polar lipid extract (Avanti Polar Lipids) and egg phosphatidylcholine in SM-KCl buffer (50 mM HEPES/Tris, pH 7.4, and 200 mM KCl). Liposomes were extruded through 400-nm filters (Whatman Nucleopore) using a syringe extruder (Avanti), and destabilized by the addition of Triton X-100 at 1:2 (wt/wt) detergent-to-lipid ratio. Labeled GltPh was added to the liposome suspension at 1:1,000 (wt/wt) protein-to-lipid ratio at room temperature for 30 min. Detergent was removed by six rounds incubation of Bio-Beads (two rounds at 23°C for 2 h each, one round at 4°C overnight, and three rounds at 4°C for 2 h each). The excess substrate was removed by three rounds of centrifugation for 1 h at 49,192 g at 4°C, removal of the supernatant, addition of 1 ml fresh SM buffer, and three freeze/thaw cycles. Liposomes were concentrated to 50 mg/ml, and ccPEB1a-Y198F labeled with activated LD555P-MAL and LD655-MAL dyes as described (Ciftci et al., 2020) was added at a final concentration of 0.6 μM and encapsulated by two freeze–thaw cycles. To remove unencapsulated ccPEB1a-Y198F, 1 ml of SM-KCl buffer was added, and liposomes were centrifuged for 1 h at 49.192 g at 4°C. Supernatants were discarded, and liposomes were suspended at 50 mg/ml and extruded 12 times through 100-nm filters.
Single-transporter smFRET assays were performed on the same microscope described above, and the microfluidic imaging chambers were prepared in the same way. After coating with streptavidin, the channel was rinsed thoroughly with SM-K(X) buffer (50 mM HEPES/Tris, pH 7.4, containing 200 mM potassium salt buffer, where the anion (X) was changed based on the condition tested). Extruded liposomes were immobilized by slowly flowing over the channel, and excess liposomes were removed by washing with 1 ml SM-K(X) buffer.
Buffers were supplemented with an oxygen-scavenging system as above, and the smFRET movies were recorded with 400-ms integration time using 20–40-mW laser power. To confirm that liposomes were not leaky, a video was taken in SM-K(X) buffer containing 1 μM L-Asp and 1 μM valinomycin. No L-Asp uptake was observed under these conditions lacking Na+ gradient. After completion of this video, another video was initiated to record transport events. At ∼3 s into this video, SM-Na(X) buffer containing 1 μM L-Asp and 1 μM valinomycin was perfused into the channel.
Traces were analyzed using Spartan software (Juette et al., 2016). Trajectories were corrected for spectral cross talk and preprocessed automatically to exclude trajectories that lasted ≤15 frames, had a signal-to-noise ratio of ≤8, and had an initial FRET efficiency <0.4 or >0.7. Traces with multiple photobleaching events (indicative of multiple sensors in a liposome) or inconsistent total fluorescence intensity were also discarded. Remaining traces were sorted by either the presence or absence of observable transport events (responding and nonresponding traces, respectively), determined by an increase in FRET efficiency from ∼0.55–0.6 to ∼0.75–0.8.
Responding traces from each video were plotted as time-dependent mean FRET efficiency. Buffer replacement time was determined as described (Ciftci et al., 2020). Within each dataset, the data were normalized so that 0% was the first time point and 100% was the first time point +0.2 (the change in FRET efficiency upon saturation of ccPEB1a-Y198F). The resulting normalized data were multiplied by the fraction of the responding traces relative to the total traces. Data from three independent reconstitutions were merged and fitted to a triexponential function in GraphPad Prism 8.4.2, where Y0 = 0, plateau = 1, percentages were set between 0 and 100, and all rate constants were set to be shared between all data sets and >0.
Online supplemental material
Fig. S1 shows the location of P-GltPh mutations and their effect on function. Fig. S2 illustrates how complex ITC isotherms change depending on parameters. Fig. S3 shows binding of TFB-TBOA to WT GltPh. Fig. S4 shows the fitted parameters of WT GltPh uptake in the presence of different anions. Figs. S5 and S6 are processing flowcharts and model validations of structures from Data S1. Figs. S7, S9, and S10 are processing flowcharts and map validations of structures from Data S2. Fig. S8 shows the effect of P-GltPh mutations on structure. Fig. S11 shows the angles of protomer tilts between OFSout, OFSmid, and OFSin. Fig. S12 shows the changes in adjacent protomers in OFSout, OFSmid, and OFSin. Fig. S13 shows the conformation of extracellular helices in OFSout, OFSmid, and OFSin. Video 1 shows different viewing angles of Fig. 3 a. Videos 2 and 3 show 3-D variability (3DVA) components of OFSout/OFSmid and OFSmid/OFSin, respectively. Video 4 shows structural transitions between OFSout, OFSmid, and OFSin. Tables S1 and S2 show fitted parameters of L-Asp binding to P-GltPh at 10°C and 15°C, respectively. Tables S3 and S4 show model validations of Data S1 and S2, respectively. Table S5 shows structural differences between structures with packing heterogeneity and OFSout, OFSmid, and OFSin. Table S6 show the tilt states arising from 3-D classification of Data S2. Data S1 provides maps and models from P-GltPh where substrate was added just before freezing. Data S2 provides maps and models from P-GltPh where substrate was present throughout the purification.
Results
P-GltPh (S279E/D405N) reveals two outward-facing substrate-binding conformations modulated by temperature and salts
We generated a mutant, P-GltPh, that eliminates Na+ binding to Na1 (D405N) and introduces a protonatable residue at the tip of HP1 (S279E), mimicking amino acid sequence features of proton-coupled orthologues (Fig. S1, a and b). Our original intention was to test the pH dependence of this mutant. While proton gradients stimulated P-GltPh activity in the presence of Na+ (Fig. S1 c), this is not the focus of this study. We observed that P-GltPh had greatly diminished transport compared to the WT transporter (Ryan et al., 2009), prompting us to test its ability to translocate substrate using smFRET. We introduced a single cysteine mutation into a cysteine-free background (P-GltPh C321A/N378C), purified the protein in DDM, labeled with donor and acceptor fluorophores, and analyzed by smFRET as in earlier studies (Akyuz et al., 2013; Akyuz et al., 2015; Huysmans et al., 2021). Conformations of protomer pairs within individual GltPh trimers can be distinguished by FRET efficiency (EFRET) as either both in OFS (∼0.4), a mixture of OFS and IFS (∼0.6), or both in IFS (∼0.8). Most P-GltPh molecules occupy a low EFRET state in the presence of 500 mM sodium salts, regardless of anion or substrate presence (Fig. 1, a and b). Thus, this mutant is mainly in OFS or the intermediate iOFS (Huang et al., 2020; Verdon and Boudker, 2012), which our smFRET setup cannot distinguish.
Further characterization of substrate binding to the mutant yielded results incompatible with our current understanding of the glutamate transporter mechanism. We expected to observe simple 1:1 binding of aspartate to a single P-GltPh binding site in ITC experiments. Instead, we observed bimodal binding isotherms in 500 mM NaCl at 15°C (Fig. 1 d). Because this mutant exclusively occupies OFS, this result suggests heterogeneity in substrate binding to this state. A two-state model assuming two independent, nonidentical binding states (OFS1 and OFS2; Freire et al., 2009) is the simplest model that fits this data reliably. Several other binding models for complex equilibria, including cooperative and sequential binding, cannot fit our data. Notably, lack of coupling between protomers has been well established (Georgieva et al., 2013; Riederer et al., 2018; Erkens et al., 2013; Ruan et al., 2017). Taken together, the data suggest there are two dominant binding states, although there could be additional underlying complexity. Furthermore, it is likely that these two states represent two different conformations of the same site in the transporter population rather than two distinct binding sites within a protomer, since the sums of the apparent stoichiometries (n1 and n2 values) averaged 0.97 ± 0.26 (range of 0.74–1.53, n = 6; Tables S1 and S2). The bimodal isotherms reflect the presence of a conformation with a higher affinity and lower exothermic binding enthalpy (OFS1) and a conformation with a lower affinity and higher exothermic enthalpy (OFS2; Brautigam, 2015; Le et al., 2013). The two conformations must interconvert only slowly, or not at all, during the ITC experiment to manifest two distinct binding states. We do not observe bimodal isotherms in 50 mM NaCl (Fig. 1 c), likely because the lower substrate affinity at lower Na+ concentrations (Reyes et al., 2013; Boudker et al., 2007) blurs the distinctions between the states.
It is possible, in principle, that heterogeneous binding reflects incomplete occupancy of the sodium-binding sites. P-GltPh does not have Na1 due to D405N mutation, and Na2 requires aspartate binding and closure of HP2. Thus, Na3 is the only site that could have incomplete occupancy before substrate binding. However, the measured Na+KD-s for WT GltPh, 99–170 mM (Riederer and Valiyaveetil, 2019; Reyes et al., 2013; Hanelt et al., 2015), suggests that the sites are saturated at 500 mM NaCl. Furthermore, if the lower-affinity OFS2 were due to incomplete occupancy of the sodium sites, increasing NaCl concentration would eliminate it. Instead, in 1 M NaCl, we observed qualitative differences in the isotherms, consistent with an increased fraction of OFS2 (Fig. 1 e). Thus, the relative populations are unlikely to depend on the occupancy of the Na3 binding site, and NaCl-dependent changes might reflect general salting effects.
To test this, we determined OFS2 fraction in the presence of 500 mM sodium salts containing anions on different ends of the Hofmeister lyotropic series (gluconate− < Cl− < NO3−; Zhang and Cremer, 2006). Gluconate and nitrate have the opposite effects on protein structure; respectively, they decrease and increase the solubility of nonpolar moieties: “salting out” and “salting in” effects. The biphasic shape of the binding isotherms is less pronounced in gluconate than chloride and nitrate (Fig. 1, d, f, and g). The fitted OFS2 fraction increases from ∼28% in NaCl to ∼38% in NaNO3 (Fig. 1 h). It decreases in Na+-gluconate, although the binding parameters were difficult to model reproducibly. Decreasing temperature to 10°C resulted in the OFS2 fraction falling to ∼20% in NaCl (Tables S1 and S2). Increasing temperatures above 15°C resulted in protein aggregation.
The observation that temperature and chaotropic salts increase the OFS2 fraction suggests that OFS1 and OFS2 are in a slow equilibrium. Because later ions in the Hofmeister series favor OFS2, the state might feature greater solvent exposure of hydrophobic regions than OFS1. It is possible that looser helical packing leads to more extensive water accessibility, with energy costs offset by increased conformational entropy. We estimated the free energy differences between OFS1 and OFS2 based on the measured populations and binding free energies at 10°C and 15°C under the assumption that the states do not exchange significantly during the ITC experiment (Fig. 1 i). OFS1 and OFS2 are nearly isoenergetic before substrate binding, but OFS1 predominates when bound to L-Asp, reflecting higher affinity. Transiency of OFS2 may explain why alternate substrate-binding conformations have not been visible in structures. Nevertheless, it must be kinetically stable over the course of ITC titrations. Notably, nanomolar apparent binding affinities measured for P-GltPh suggest that substrate dissociation contributes little to the binding process observed in ITC. Thus, the two states may differ primarily in the binding on-rates, with the high-affinity state binding substrate faster than the low-affinity state. Regardless, these results strongly suggest that there is conformational heterogeneity in the transporter, manifesting in different binding mechanisms.
Heterogeneous substrate binding in WT GltPh
We also performed ITC experiments on the WT protein. When we used high protein concentrations to increase experimental sensitivity, we observed unusual features in binding isotherms (Fig. 2, a–c). Specifically, the initial L-Asp injections do not have constant heats as expected for a high-affinity single-site binding process. Instead, injection heats steadily decrease until an abrupt drop occurs when the ligand saturates the protein. The two-state model, where the two affinities are close but not identical and the higher-affinity state has a higher exothermic binding enthalpy, fits data well, though the binding parameters are not uniquely determined (Fig. S2 a). As in P-GltPh, we observed qualitative differences in the isotherms in different salts and temperatures. The isotherm in 500 mM Na-gluconate at 15°C has a particularly unusual shape (Fig. 2 a) but becomes more reminiscent of a single-site binding in more chaotropic salts (Fig. 2, b and c) or at higher temperatures (Fig. 2, d and e). As in P-GltPh, changes in the state populations can account for the changing isotherm shapes (Fig. S2 b), though faster exchange between states or altered enthalpies could also contribute. Collectively, our data suggest that WT GltPh has multiple binding states in temperature- and salt-modulated equilibrium.
WT transporter in saturating Na+ concentrations is predominantly in the OFS (Akyuz et al., 2015; Akyuz et al., 2013), but we cannot exclude that inward-facing protomers contribute to heterogeneous L-Asp binding. However, binding of the transport blocker TFB-TBOA, which has a 100-fold preference for OFS (Mcilwain et al., 2016; Boudker et al., 2007; Wang and Boudker, 2020), also produced bimodal ITC isotherms reminiscent of P-GltPh (Fig. S3, a–c). Due to the lower TFB-TBOA affinity, we could not determine precise binding parameters and quantify the salt effects. When GltPh saturated with TFB-TBOA was competed with L-Asp, we again observed bimodal isotherms (Fig. S3, d–f). Thus, the inhibitor binding is heterogeneous, and some level of conformational heterogeneity persists after binding.
We used a recently developed smFRET-based single-transporter assay to test if salt-modulated state populations correlated with transport rates (Ciftci et al., 2020). P-GltPh C321A/N378C mutant was labeled with PEG11-biotin and N-ethylmaleimide and reconstituted into liposomes; low protein-to-lipid ratios enriched vesicles containing at most one GltPh trimer. The proteoliposomes were then loaded with periplasmic glutamate/aspartate binding protein (PEB1a) Y198F/N73C/K149C mutant with reduced aspartate affinity labeled with maleimide-activated donor (LD555P) and acceptor (LD655) fluorophores (referred to altogether as ccPEB1a-Y198F). The proteoliposomes were immobilized in microscope chambers via biotinylated transporter and assayed for transport after perfusion with saturating Na+ and L-Asp concentrations. An increase in mean EFRET from ∼0.6 to ∼0.8 reflects saturation of the ccPEB1a-Y198F sensor by L-Asp molecules transported into vesicles. This assay previously established kinetic heterogeneity in WT GltPh transport (Ciftci et al., 2020), where at least three observable populations (fast, intermediate, and slow) transport at vastly different rates, and all contribute to mean uptake measured in bulk. Most WT transporters are slow, with turnover times of tens to hundreds of seconds.
Because GltPh mediates an uncoupled anion conductance, which dissipates the buildup of membrane potential due to electrogenic transport, it shows faster uptake in the presence of more permeant anions (gluconate− < Cl− < NO3−; Ryan and Mindell, 2007). Thus, we measured transport in K+-loaded proteoliposomes in the presence of ionophore valinomycin clamping the potential. Triexponential fits of the uptake kinetics suggest that most molecules are in the slow transporting population regardless of the salt, as expected. The fractions of the slow transporters were similar in Na-gluconate (81.1 ± 3.1) and NaCl (79.5 ± 3.4%) conditions but increased in NaNO3 conditions (87.3 ± 2.2%; Fig. 2 f and Fig. S4, a and b). Therefore NaNO3-favored OFS2 conformation might correlate with a slower transporter population.
Transient transport domain structures following substrate binding
smFRET has shown that P-GltPh is exclusively outward facing, making this mutant an excellent model to dissect differences between OFS1 and OFS2 using cryo-EM, where structural heterogeneity should reflect the binding heterogeneity. Because OFS2 is transient and OFS1 predominates at equilibrium, we optimized conditions to increase the probability of imaging OFS2. ITC analysis showed that elevated temperatures and chaotropic salts increased the OFS2 fraction (Fig. 1). Thus, we pre-equilibrated P-GltPh in 250 mM NaNO3 at 25°C and froze grids within ∼5 s after adding 1 mM L-Asp (Data S1).
When we refined particles with imposed C3 symmetry, we obtained density maps with an overall resolution of 3.0 Å (Fig. S5). To maximally retain heterogeneity, we used a data processing approach designed to pick the highest-quality particles regardless of conformation (Su et al., 2020). We then used symmetry expansion and focused classification of single protomers into 10 classes followed by local refinement, a processing approach that previously revealed OFS and iOFS in the WT GltPh ensemble (Huang et al., 2020). We did not find any iOFS classes and observed only OFS classes with similar overall structures. We refined models for four classes with the highest resolution, from 3.15 to 3.85 Å (Figs. S5 and S6, and Table S3). When we superimposed their isolated transport domains on the intracellular regions (HP1, TM8b, and TM7a), below the substrate-binding site, they aligned well (Fig. 3 a). In contrast, we observed displacements of helices in the extracellular halves above the substrate-binding site, most noticeable in HP2, TM8a, and TM7b (Fig. 3 a and Video 1). These observations suggest heterogeneity in packing of the extracellular half of the transport domain immediately after substrate binding. Notably, we did not observe any differences between the Na3 sites of these structures, further suggesting that heterogeneous substrate binding does not result from incomplete Na3 occupancy. Similar superpositions of transport domains from the previously reported structures of substrate-, inhibitor-, and Na+-only bound GltPh (Alleva et al., 2020; Boudker et al., 2007; Yernool et al., 2004) also picture differences in positions of TM7b and TM8a in addition to the expected differences in gating HP2 (Fig. 3 b). These observations further support conformational lability of these helices and suggest that ligand binding entails the restructuring of the entire region.
High-resolution equilibrated structures of P-GltPh
ITC experiments suggest that the low-affinity OFS2 should be transient in the presence of substrate, and the high-affinity OFS1 should predominate at equilibrium (Fig. 1 i). To visualize OFS1, we imaged P-GltPh in equilibrium conditions, purified in the presence of 250 mM NaNO3 and 1 mM L-Asp (Data S2). We refined the maps to 2.2 Å after imposing C3 symmetry (Fig. S7). The increased resolution compared with Data S1 may reflect reduced structural heterogeneity but can also be due to different microscopes, imaging parameters, or variations in grid preparations. The D405N mutation abolishes Na+ binding at Na1 (Boudker et al., 2007; Riederer et al., 2021 Preprint); in its place, we observed an excess density, suggesting that a water molecule replaces the ion (Fig S8 a). The S279E side chain points into the extracellular milieu, away from the transport domain (Fig. S8 b).
We found that the equilibrated transport domain structure is nearly identical to class A3 (with overall RMSD of 0.30 Å). Thus, we propose that this conformation corresponds to the higher-affinity OFS1 substate. The ensemble of other conformations observed in Data S1 (A1, A6, and A7), showing different helical packing of the transport domain, together make up the low-affinity OFS2 substate. Previously described V366A and A345V mutations in HP2 are on the interface between the hairpin and TM8a and TM7b helices, where they can disrupt the helical packing. Indeed, HDX-MS measurements suggested increased local dynamics in this region in Y204L/A345V/V366A GltPh mutant. Consistent with our hypothesis, these mutations also decreased substrate affinity, even though the crystal structure of the mutant pictured preserved substrate coordination and binding site details (Huysmans et al., 2021; Ciftci et al., 2021).
Interestingly, we observed several water densities within the transport domain of the equilibrated structure contributing to the hydrogen bond network between substrate- and ion-coordinating residues. All six resolved buried water molecules are below the substrate-binding site. In contrast, the extracellular half of the domain, corresponding to the labile helices in Data S1, appears “dry” (Fig. 3 c). We speculate that the extensive hydrogen bond network in the cytoplasmic half of the transport domain ensures its rigidity. In contrast, the extracellular half, less constrained by polar interactions, can sample multiple conformations with altered packing. Notably, chaotropic salts and elevated temperature favor OFS2 consistent with less well-packed, more dynamic, water-accessible structures.
We also looked for structural heterogeneity in Data S2. 3DVA analysis on a single protomer using the symmetry-expanded particle stack (Punjani and Fleet, 2021) revealed small movements of the transport domain (Videos 2 and 3). Focused 3-D classification of ∼1.6 million symmetry-expanded particles showed only protomers in OFS and yielded structural classes corresponding to the density variations seen in 3DVA analysis. We refined these classes to 2.36–2.65 Å resolutions (Figs. S9 and S10; and Table S4). Superimpositions of the refined trimers on trimerization regions (residues 150–195) showed three subtly different tilts of the classified protomers, consisting of movements of the transport domain and the peripheral parts of the scaffold (OFSout, OFSmid, and OFSin; Video 4). The largest tilt difference of 2.1° is between OFSout and OFSin transport domains. The tilt differences for OFSout/OFSmid and OFSmid/OFSin were ∼1.1° each (Fig. S11). The adjacent protomers are unaffected, suggesting that the movements occur independently in individual protomers (Video 4 and Fig. S12). Notably, we observed no rearrangements of the extracellular helices regardless of the tilts, and all tilt states most closely resembled class A3 in Data S1, consistent with the high-affinity OFS1 predominating at equilibrium (Fig. S13 and Table S5). The mechanistic basis of the tilts and their role in transport remain unclear. Similar transport domain tilts might also be present in Data S1; however, the moderate resolution of the maps prevents their visualization.
Analysis of the substrate-binding sites of Data S2 classes revealed that Asp-390 sampled multiple rotameric states. The highly conserved Asp-390 in TM8 does not coordinate the substrate but is critical for high-affinity binding—D390A mutant has a 1,000-fold lower affinity (Riederer and Valiyaveetil, 2019). Arg-397 in TM8 is the principal substrate-coordinating residue. Its guanidinium group forms hydrogen bonds with the L-Asp sidechain carboxylate and cation-π interactions with Tyr-317 in TM7. We found that Asp-390 can be in down or up rotamers, hydrogen-bonding to Arg-397 or Tyr-317, respectively (Fig. 4, a and b). Classifications also revealed a middle rotamer, perhaps representing an average of the two rotamers or a unique state (Fig. 4 c). Different Asp-390 rotamers do not result in an observable change of Arg-397 conformation but might alter the local electrostatics. Furthermore, tyrosine hydrogen bonding through the OH group potentiates cation-π interactions compared with phenylalanine (Gallivan and Dougherty, 1999), and Y317F mutation leads to a 10-fold loss of L-Asp affinity in GltPh (Riederer and Valiyaveetil, 2019). Thus, the up and down rotamers might alter substrate affinity. After additional rounds of sorting, we found that only subpopulations of OFSout, comprising ∼10% of all particles, featured up or middle rotamers (Fig. 4 c, Table S6, and Fig. S9). Thus, the preference of Asp-390 to hydrogen-bond with Arg-397 or Tyr-317 might be allosterically coupled to the position of the transport domain. Whether these structural heterogeneities contribute to the elevator dynamic or substrate binding heterogeneities is yet unclear.
Discussion
Serendipitously, we found that P-GltPh mutant has two slowly exchanging outward-facing conformational substates, OFS1 and 2, binding aspartate with different affinities and enthalpies (Fig. 1). smFRET and cryo-EM showed that P-GltPh is predominantly outward facing (Fig. 1, a and b). Thus, P-GltPh is an excellent model to consider the mechanism of heterogeneous binding in OFS. WT GltPh is less suitable because its conformational ensemble includes iOFS and IFS (Wang and Boudker, 2020; Reyes et al., 2013; Akyuz et al., 2013). Nevertheless, WT binding isotherms also suggest multiple binding conformations modulated by salts and temperature (Fig. 2). Interestingly, a recent saturation transfer difference NMR study reported an unusually low Hill coefficient of 0.69 for aspartate binding to liposome-reconstituted GltPh (Hall et al., 2020). A Hill coefficient below one may reflect negative cooperativity or result from multiple binding states with distinct affinities (Cattoni et al., 2015; Sevlever et al., 2020; Wang and Pan, 1996). Distinguishing these possibilities requires a kinetic approach (Cattoni et al., 2009).
Recent studies showed that GltPh, originating from a hyperthermophilic archaeon, exhibits activity modes at ambient temperature, where transporter subpopulations function with rates differing by orders of magnitude (Ciftci et al., 2020). Switching between the modes is rare, occurring on a timescale of hundreds of seconds. These modes were attributed to subpopulations with different elevator dynamics and intracellular substrate-release rates (Matin et al., 2020; Huysmans et al., 2021; Ciftci et al., 2021). Here, we observed that NaNO3 modulated the populations of the binding substates and increased the fraction of the slow transporters (Fig. 2 f). Therefore, the substrate-binding heterogeneity might contribute to the heterogeneous uptake rates or reflect different substrate-binding properties of fast and slow subpopulations. Some gain-of-function mutations, including R276S/M395R, A345V, and V366A, both increase GltPh elevator dynamics and reduce substrate affinity (Huysmans et al., 2021; Ciftci et al., 2021). Thus, structural perturbations can affect binding and translocation in concert, suggesting that the two processes share some of the determinants. Therefore, insights into the structural bases of heterogeneous binding might also illuminate the origins of heterogeneous transport kinetics.
Our cryo-EM structures of P-GltPh, imaged in conditions nearly identical to the ITC experiments, demonstrate that the extracellular half of the transport domain has a continuum of packing states upon binding (Fig. 3, a and b). Yet, after equilibration, we observed no such heterogeneity. We found no other pronounced structural heterogeneities that exist immediately after binding but are gone after equilibration. Furthermore, our structural classes do not display any evidence of altered substrate coordination that could explain affinity differences. Thus, heterogeneous binding observed in ITC correlates with different helical packing observed in cryo-EM structures. We propose that structural flexibility in the extracellular regions, including the gating HP2 and adjacent helices, alter the energetics of substrate binding. We further suggest that the structural class with the optimal packing corresponds to the high-affinity substate, and the classes with the continuum of helical packing arrangements together correspond to the low-affinity substate. In this model, following substrate binding, the transporter relaxes to the optimally packed state on a very slow time scale. The equilibration process involves helical repacking and might be slow because it requires rearrangement of a large transport domain or substrate release and rebinding. Consistently with our hypothesis, packing mutations A345V and V366A in the extracellular helices lead to decreased affinity and increased binding and dissociation kinetics (Huysmans et al., 2021; Ciftci et al., 2021). Notably, the closure of the HP2 tip over the substrate-binding site is unlikely to contribute to binding heterogeneity because it is a rapid process (Riederer and Valiyaveetil, 2019), and achieving tight binding requires a kinetically slow step (Ewers et al., 2013; Hanelt et al., 2015). Thus, we see helical repacking as the most, if not the only, feasible explanation of the experimentally observed heterogeneous binding.
Our structural analysis suggests that there is a structural specialization within the transport domain of glutamate transporters. A network of buried waters and polar residues within the cytoplasmic half of the domain (TM8b, TM7a, and HP1) might ensure exquisitely specific, evolutionarily conserved structure responsible for sodium selectivity and allosteric coupling between ion and substrate binding (Fig. 3 c). In contrast, the extracellular half (TN7b, HP2, and TM8a) is more hydrophobic and contains no resolved waters. Fewer constraining polar interactions may permit variable helical packing in this region, setting the dynamic properties—substrate affinities and elevator dynamics—of GltPh substates and, perhaps, different glutamate transporter homologues. Consistently, extensive mutagenesis in HP1 and adjacent helices did not identify gain-of-function mutants analogous to those in HP2 (Huysmans et al., 2021). These observations of structural specializations are reminiscent of “protein sectors,” where proteins can have discrete, independently evolving functional units in tertiary structure (Halabi et al., 2009).
The structural differences between the high- and low-affinity states are small and require analysis of high-resolution cryo-EM datasets. Similarly, modal gating in KcsA was attributed to subtle sidechain rearrangements (Chakrapani et al., 2011). Functional heterogeneity has been observed in ion channels, transporters, and enzymes. For example, nicotinic acetylcholine receptors feature opening bursts interspersed with short or long closed periods (Colquhoun and Sakmann, 1985), and ionotropic glutamate receptors show complex kinetics with multiple gating modes (Popescu, 2012). P-type ATPase also displayed periods of rapid transport interspersed with prolonged pauses (Veshaguri et al., 2016). As in GltPh, these distinct modes may be due to concerted subtle restructuring of protein regions occurring on long timescales.
Acknowledgments
Joseph A. Mindell served as editor.
We thank Drs. Eva Fortea, Maria Falzone, Philipp Schmidpeter, Biao Qiu, and Xiaoyu Wang for helpful discussions on cryo-EM data processing. We especially thank Navid Paknejad for in-depth assistance with optimizing cryo-EM processing workflows. We thank the Scientific Computing Unit at Weill Cornell Medical College for maintenance and support of computational resources. Also, we thank Bryce Delgado for preliminary ITC experiments, Will Eng for protein expression, and Vishnu Ghani for preparation of the PEB1a protein. We thank Dr. Scott Blanchard and members of the Blanchard lab for support and resources for smFRET experiments. Finally, we thank Drs. Erika Riederer, Francis Valiyaveetil, and Yun Huang for helpful discussions and exploratory experiments.
The work was supported by National Institutes of Health (NIH) grant F32 NS102325 (to K.D. Reddy), American Heart Association grant 19PRE34380215 (to D. Ciftci), and NIH grants R01NS064357 and R37NS085318 (to O. Boudker). Data S1 was collected with assistance from Carolina Hernandez at the Simons Electron Microscopy Center and National Resource for Automated Molecular Microscopy located at the New York Structural Biology Center, supported by grants from the Simons Foundation (349247), NYSTAR, and the NIH National Institute of General Medical Sciences (GM103310) with additional support from Agouron Institute (F00316) and NIH S10 OD019994-01. Data S2 was collected at the UMass cryo-EM facility with help from Dr. Kangkang Song and Dr. Chen Xu.
The authors declare no competing financial interests.
Author contributions: K.D. Reddy and O. Boudker designed the experiments, analyzed the data, refined the molecular models, and wrote the manuscript with input from all authors. K.D. Reddy and A.J. Scopelliti performed the cloning and radioactive transport assays. K.D. Reddy and D. Ciftci performed the smFRET dynamics and transport assays. K.D. Reddy performed the ITC and cryo-EM sample preparation and processing.