Zachary Lee Johnson1, Cheom-Gil Cheong, Seok-Yong Lee. 1. Department of Biochemistry and Ion Channel Research Unit, Duke University Medical Center, 2 Genome Court, Durham, North Carolina 27710, USA.
Abstract
Nucleosides are required for DNA and RNA synthesis, and the nucleoside adenosine has a function in a variety of signalling processes. Transport of nucleosides across cell membranes provides the major source of nucleosides in many cell types and is also responsible for the termination of adenosine signalling. As a result of their hydrophilic nature, nucleosides require a specialized class of integral membrane proteins, known as nucleoside transporters (NTs), for specific transport across cell membranes. In addition to nucleosides, NTs are important determinants for the transport of nucleoside-derived drugs across cell membranes. A wide range of nucleoside-derived drugs, including anticancer drugs (such as Ara-C and gemcitabine) and antiviral drugs (such as zidovudine and ribavirin), have been shown to depend, at least in part, on NTs for transport across cell membranes. Concentrative nucleoside transporters, members of the solute carrier transporter superfamily SLC28, use an ion gradient in the active transport of both nucleosides and nucleoside-derived drugs against their chemical gradients. The structural basis for selective ion-coupled nucleoside transport by concentrative nucleoside transporters is unknown. Here we present the crystal structure of a concentrative nucleoside transporter from Vibrio cholerae in complex with uridine at 2.4 Å. Our functional data show that, like its human orthologues, the transporter uses a sodium-ion gradient for nucleoside transport. The structure reveals the overall architecture of this class of transporter, unravels the molecular determinants for nucleoside and sodium binding, and provides a framework for understanding the mechanism of nucleoside and nucleoside drug transport across cell membranes.
Nucleosides are required for DNA and RNA synthesis, and the nucleoside adenosine has a function in a variety of signalling processes. Transport of nucleosides across cell membranes provides the major source of nucleosides in many cell types and is also responsible for the termination of adenosine signalling. As a result of their hydrophilic nature, nucleosides require a specialized class of integral membrane proteins, known as nucleoside transporters (NTs), for specific transport across cell membranes. In addition to nucleosides, NTs are important determinants for the transport of nucleoside-derived drugs across cell membranes. A wide range of nucleoside-derived drugs, including anticancer drugs (such as Ara-C and gemcitabine) and antiviral drugs (such as zidovudine and ribavirin), have been shown to depend, at least in part, on NTs for transport across cell membranes. Concentrative nucleoside transporters, members of the solute carrier transporter superfamily SLC28, use an ion gradient in the active transport of both nucleosides and nucleoside-derived drugs against their chemical gradients. The structural basis for selective ion-coupled nucleoside transport by concentrative nucleoside transporters is unknown. Here we present the crystal structure of a concentrative nucleoside transporter from Vibrio cholerae in complex with uridine at 2.4 Å. Our functional data show that, like its human orthologues, the transporter uses a sodium-ion gradient for nucleoside transport. The structure reveals the overall architecture of this class of transporter, unravels the molecular determinants for nucleoside and sodium binding, and provides a framework for understanding the mechanism of nucleoside and nucleoside drug transport across cell membranes.
Humans have three isoforms of CNT and the substrate specificity and tissue distribution are different among the isoforms [14-16]. Knowledge of the mechanism of these transporters would help us not only to understand physiological processes associated with nucleosides but also to provide a framework for future drug design to improve nucleoside drug delivery. The major barrier towards achieving a mechanistic understanding is the lack of atomic structures that reveal the origins of nucleoside specificity and the principles of function of CNTs. For structure determination and functional studies we chose a CNT homolog from Vibrio cholerae (vcCNT) due to its high sequence homology to hCNTs (39% identical to hCNT3, Supplementary Figure 1) and optimal biochemical stability.To test whether vcCNT can transport nucleosides, and if so what ion gradient vcCNT utilizes, we carried out a radioactive nucleoside uptake assay using recombinant vcCNT-reconstituted liposomes. Studies have shown that the human CNTs use Na+ and the E. coli CNT uses H+ for nucleoside transport [17-19]. Transport activity was measured by monitoring uptake of [5,6-3H]-uridine both in the presence and absence of a Na+ gradient (Figure 1a). We chose uridine because all of the CNTs that have been characterized to date transport uridine [14,17-19]. When a Na+ gradient was present, vesicles containing vcCNT accumulated radioactive uridine significantly more than the control empty vesicles. When a Na+ gradient was not present (Na+ is replaced with choline in the solution), uridine uptake was significantly reduced. Further characterization showed that uridine uptake is electrogenic and independent of pH gradients (Supplementary Figure 2). These data clearly demonstrate that vcCNT, like human CNTs, utilizes a Na+ gradient to transport nucleosides.
Figure 1
vcCNT is a Na+-coupled nucleoside transporter with a trimeric architecture. a, Time course of the uptake of 2.4 μM [5,6-3H]-uridine into vesicles containing vcCNT in the presence of a sodium gradient (squares, NaCl), in the absence of a sodium gradient (triangles, ChCl), and the control empty vesicles in the presence of sodium (diamonds, empty). b, Cartoon representation of the vcCNT trimer viewed from the cytoplasm. Individual protomers are colored blue, red, and green. c, Cartoon representation for the vcCNT trimer viewed parallel to the membrane. The putative membrane bilayer is indicated in c with grey boxes. The dimensions of the putative membrane bilayer, extracellular region, and intracellular region are shown.
We crystallized and solved the structure of vcCNT at 2.4 Å (Supplementary Table 1). Experimental phases to 3.5 Å were obtained by single anomalous dispersion (SAD) from a platinum-soaked crystal (see Methods). The final model contains a single protomer in the asymmetric unit and is of good quality with a free R factor of 22.8% (Supplementary Table 1). vcCNT crystallizes as a homotrimer that is shaped like an inverted triangular basin with its mouth facing the intracellular side and a knob-like structure facing the extracellular side (Figure 1b and 1c). The three-fold axis coincides with the crystallographic 6-fold axis and is perpendicular to the membrane. When viewed from the intracellular side each side of the triangle formed by the trimer is ~92 Å, and when viewed parallel to the membrane the trimer is ~57 Å in height (Figure 1). The membrane-embedded region lies approximately in the middle of the transporter judging from the positions of three amphipathic helices on each protomer. Based on the predicted location of the membrane bilayer, the mouth of the basin penetrates into the membrane plane (Figure 1c and Supplementary Figure 3). The presence of several polar amino acids on the basin surface likely allows bulk aqueous solution to reach deep into the membrane bilayer and access this surface of the transporter.The structure reveals that each protomer contains 8 transmembrane helices (TM1-8), two reentrant helix-turn-helix hairpins (HP1 and HP2) with opposite orientations in the membrane, and three interfacial helices (IH1-3) that run parallel to the membrane (Figure 2a). The orientation of the structure relative to the membrane is consistent with previous accessibility studies [20-22] and conforms to the positive-inside rule (Supplementary Figure 4). hCNTs are predicted to contain three additional N-terminal TMs compared to prokaryotic CNTs, suggesting an 11-TM topology for hCNTs [20].
Figure 2
Topology and fold of the vcCNT protomer. a, Schematic representation of vcCNT topology. The group of helices under the pink triangular background is related to the group of helices under the cyan triangular background by two-fold pseudo-symmetry with the symmetry axis parallel to the membrane. b, Cartoon representation of the vcCNT protomer fold. Only helices comprising the scaffold domain and the group of helices under the pink triangular background are labeled. c, Same with b but rotated 120°. Only TM1, TM6, and helices comprising the group of helices under the cyan triangular background are labeled. d, Same with c but rotated 120°. The scaffold domain is removed to show the pseudo-two-fold relationship between the two groups of helices colored as pink and cyan triangular backgrounds. The two-fold symmetry axis runs through the nucleoside at the center of the transporter.
Each protomer can be grouped into two subdomains based on their locations (at the outer and inner regions) relative to the center of the protomer. TM1, TM2, IH1, EH (extracellular helices), TM3, and TM6 are located at the outer part of the transporter and they appear to be important for maintaining the overall architecture of the transporter. These helices in the outer region form a scaffold for the transporter (henceforth termed the scaffold domain) (Figure 2). Trimerization contacts are mediated by part of the scaffold domain - IH1, EH, TM3, and TM6 (Figure 2b). IH1 is a 40 Å-long amphipathic helix that is most likely situated at the water-membrane interface, and TM6 is ~60 Å long and tilted almost 60° with respect to the membrane normal. EH protrudes into the extracellular solution and is ~16 Å in length. Since IH1 is involved in trimerization, is constrained to be at the membrane-water interface, and interacts with many TMs (TM2-TM5), it most likely plays an important role in building and maintaining the overall structure of vcCNT. Surrounded by the scaffold domain, many conserved amino acids implicated in nucleoside transport are localized at the inner domain (henceforth termed the transport domain). The transport domain is composed of two structural groups that are related by an internal two-fold pseudo-symmetry. The first group is composed of IH2, HP1, TM4a/b, and TM5 (pink background in Figure 2a), and the second group is composed of IH3, HP2, TM7a/b, and TM8 (cyan background in Figure 2a). These two groups, separated by TM6 in sequence, can be superimposed with an r.m.s.d. of 3.3 Å with the two-fold symmetry operator running though the middle of TM6 and parallel to the membrane (Figure 2d). There is no significant amino acid sequence homology between these two groups (~10% sequence identity). This two-fold relationship positions the tips of HP1, HP2, and unwound regions of TM4 and TM7 at the center of the transport domain, which is located slightly below the middle of the membrane plane (Figure 2d). Sequence alignment of hCNTs and vcCNT reveals high sequence conservation around HP1, HP2, and the two unwound helices (TM4 and TM7), indicating the functional importance of this region (Supplementary Figure 1).To the best of our knowledge, the overall fold of vcCNT is novel although local structural elements such as helical hairpins and unwound helices have been observed previously [23,24].The crystal structure suggests that vcCNT adopts a trimeric configuration. To test whether the stoichiometry of vcCNT is trimeric and further validate the physiological relevancy of our crystal structure, we performed structure-guided disulfide bridge cross-linking experiments (Supplementary Figure 5). Our cysteine mutants readily form disulfide cross-linked trimers in both detergent micelles and cell membranes under oxidizing conditions, therefore our crystal structure reflects a physiologically relevant oligomerization state. Given the sequence conservation of most of the amino acids involved in trimerization, we propose that the stoichiometry of both eukaryotic and prokaryotic CNT family members is trimeric (Supplementary Figure 1).The structure contains three deep clefts (one per subunit) at the intracellular side facing the center of the trimer (Figure 3a). A simulated-annealing OMIT map clearly shows that the electron density in the cleft is that of uridine (Figure 3b). The bound uridine in the cleft faces the intracellular basin of the transporter; however, it is not free to be released into the intracellular solution because TM6 and TM7b partly cover the binding site (Figure 3a and 3b and Supplementary Figure 6).
Figure 3
Nucleoside binding site and sodium binding site. a, Cut away surface representation of vcCNT viewed parallel to the membrane. The putative membrane bilayer is shown with lines. Uridine is shown as spheres. b, The vcCNT protomer viewed from the center of the trimer. The scaffold domain is shown in ribbon representation and the transport domain is shown in cartoon representation. The blue mesh, covered by TM6 (green), is an Fo-Fc simulated annealing OMIT map, contoured at 4 σ, showing density for uridine. The inset is a zoomed-in image of the simulated annealing OMIT electron density. c, The nucleoside-binding site showing HP1 (red), TM4b (orange), HP2 (blue), and TM7b (cyan). Hydrogen bonds are shown as dashed lines. d, The sodium-binding site is located between HP1 (red) and the unwound region of TM4 (orange). The blue mesh is an Fo-Fc simulated annealing OMIT map, contoured at 6 σ, showing density for sodium. Coordination of the sodium ion is shown as dashed lines. e, The sodium-binding site is near the nucleoside-binding site.
The nucleoside-binding site is located at the center of the internal two-fold symmetry and lined by the tips of HP1 and HP2 and the unwound regions of TM4 and TM7 (Figure 3b and 3c). Inspection of the interactions between uridine and the binding site shows that many polar or charged amino acids from HP1, HP2, TM4, and TM7 interact with the uracil base and ribose. HP1 and TM4b are responsible for interacting with the uracil base. The side chains of Gln 154, Thr 155, and Glu 156 from HP1 interact with the uracil base either directly (Gln 154) or indirectly through a water molecule (Thr 155 and Glu 156) (Figure 3c). Val 188 from TM4b interacts with the uracil base via van der Waals interactions. The involvement of these amino acid residues in the interactions with the nucleoside is consistent with mutational studies by Young and colleagues: the residue corresponding to Gln 154, together with that corresponding to Val 188, is important for nucleoside specificity of hCNTs; Glu 156 is critical for Na+-nucleoside coupled transport [25-27].HP2 and TM7 are responsible for the interactions with ribose. The side chains of Glu 332 (HP2), Asn 368 (TM7), and Ser 371 (TM7) interact with the ribose either directly (Glu 332 and Ser 371) or indirectly through a water molecule (Asn 368). Mutation of the residue corresponding to Glu 332 in hCNTs was shown to have significant functional effects on both nucleoside binding and the rate of transport [21,26].Since vcCNT is a sodium-coupled transporter, it must contain at least one sodium-binding site. Initial hints regarding the location of the sodium-binding site came from an Fo-Fc map that shows a strong peak (~8σ) between the tip of HP1 and the unwound region of TM4. Further inspection reveals that this site is octahedrally coordinated by three backbone carbonyls, two side chain hydroxyl groups, and a water molecule (Figure 3d) with the typical distances (~2.4 Å) for Na+ coordination by protein [28]. Since it is difficult to distinguish Na+ from H2O based on electron density alone, we performed the valence test that was used to identify the Na+ sites in the crystal structure of LeuT [28,29]. The result (ν = 1.02) is consistent with the bound ion in this site being Na+. Therefore we suggest that vcCNT has at least one Na+-binding site at this position between HP1 and the unwound region of TM4.The sodium-binding site is near the nucleoside-binding site, but the Na+ does not directly interact with the nucleoside (Figure 3e). What is the functional role of Na+ binding at this site? Since key amino acid residues for the interactions with the nucleoside base are localized on HP1 and TM4b, we hypothesize that the role of Na+-binding is to bring HP1 close to TM4 to complete the formation of the nucleoside-binding site.In the field of sodium-coupled secondary transporters, the current paradigm revolves around the alternating access mechanism [30]. Assuming the alternating access model is utilized, our crystal structure probably represents an inward-facing occluded conformation (Figure 3a and Supplementary Figure 6). On the bases of our crystal structure, previous accessibility studies of hCNTs, and studies of the sodium-coupled aspartate transporter GltPh that shares certain structural features with vcCNT (Supplementary Figure 7), we propose a hypothetical model to explain sodium-coupled nucleoside transport by vcCNT (Figure 4). A rigid-body movement of the transport domain across TM6 while the scaffold domain is held in place could allow for the transition from an inward-facing to an outward-facing conformation. Our reasoning in proposing this mechanism is provided in the Supplementary Information.
Figure 4
Hypothetical mechanism of nucleoside transport. Only a single protomer is shown for simplicity, viewed from the center of the trimer (same orientation as 3b). On the left, the transporter adopts an inward-facing conformation. Uridine is bound between the tips of HP1 and HP2, and TM6 (grey), serving as a hydrophobic barrier, partially covers uridine. The blue dotted lines demarcate the mobile part of the transport domain and the grey region denotes the scaffold domain. A rigid-body motion of the nucleoside-binding site across TM6 can expose the binding site to the extracellular side.
The crystal structure provides an important step toward understanding the mechanism of nucleoside transport by the CNT family. However, many important issues remain to be elucidated including the structural bases for nucleoside and nucleoside-drug specificity and the conformational changes associated with each step along the transport cycle.
Methods Summary
vcCNT was expressed in E. coliC41 cells. Cells were lysed and protein was solubilized with dodecyl maltoside. Protein was purified by Co2+ affinity chromatography followed by gel filtration chromatography in the presence of decyl maltoside and uridine.vcCNT was crystallized in the presence of 100 mM CaCl2, 37–42% PEG400, and 100 mM buffer (pH 5.6–9.5). Platinum (Pt) derivatives were prepared by soaking in K2Pt(CNS)6. A partial poly ALA model was built with phases to 5.0 Å from single isomorphous replacement with anomalous scattering (SIRAS). A complete model was built by combining phases from single anomalous dispersion using a new platinum derivative with molecular replacement phases from the partial model (MR-SAD) to 3.5 Å. Structure refinement was then performed against the 2.4 Å native data. The final model is of good quality with good Ramachandran statistics (98% favored and 2% allowed). The model contains residues 2 to 416 (residues 230 to 240 missing), a uridine, 5 decylmaltosides (1 full and 4 partial), a sodium ion and 43 bound waters. Disordered residues Lys 226 and Glu 345 were modeled as Ala.For the flux assay, protein was reconstituted into lipid vesicles and then vesicles were diluted into buffer containing either choline chloride or NaCl. Flux was initiated by the addition of 2.4 μM [5,6-3H]-uridine and 1 μM valinomycin.
Full methods
Expression and Purification
The vcCNT gene was cloned from Vibrio cholerae genomic DNA into a modified pET26 vector that contains a pelB leader sequence and a PreScission Protease cleavable His10-MBP fusion. The original vector was a gift from Raimund Dutzler (University of Zurich) and was further modified. The vector containing vcCNT was expressed in E. coliC41 (DE3) cells. Cells were lysed by homogenizer (AVESTIN) and protein was solubilized with 30 mM dodecyl maltoside (DDM) for two hours at 4 °C. Solubilized lysates were spun down to remove the insoluble fraction and supernatants were applied to Talon Co2+ affinity resin. After binding, protein was eluted with imidazole and digested overnight with PreScission Protease. The digestion mixture was concentrated and applied to a Superdex 200 size exclusion column in the presence of 5 mM decyl maltoside (DM), 1 mM uridine, 150 mM NaCl, 20 mM Tris pH 8.0, and 2 mM dithiothreitol. Peak fractions corresponding to vcCNT were collected for further experiments.
Crystallization
Initial crystallization conditions were obtained at the high-throughput crystallization-screening laboratory at the Hauptman-Woodward Institute [31]. Protein was concentrated to ~10 mg/mL and mixed 1:1 with crystallization solution containing 100 mM CaCl2, 37–42% PEG400, and 100 mM buffer. Crystals grew over a wide pH range (5.6–9.5), but data for structure solution was collected on crystals grown at pH 9.0 (100 mM Tris) and pH 9.5 (100 mM glycine). Crystals were grown using the microbatch-under-oil method. After 10–14 days, crystals were harvested, transferred to cryo solution containing 32.5% PEG400, and flash frozen in liquid nitrogen. Platinum derivatives were prepared by soaking in 2.5 mM K2Pt(CNS)6 for 2–4 h then transferring to cryo solution and flash freezing.
Structure Determination
The data were collected on beam lines 22ID, 22BM, and 24ID-C at the Advanced Photon Source. The data were processed using HKL2000[32]. Crystals of vcCNT diffract to 2.4 Å Bragg spacings and belong to the space group P63. Extensive screening of crystals was performed because approximately 70% of crystals are merohedrally twinned with the twinning operator perpendicular to the crystallographic 6-fold axis leading to the apparent space group of P6322. We found that SeMet-substituted crystals are almost always twinned with significant twinning fractions (20–45%), which necessitated heavy atom-soaked crystals as the choice for a de novo phasing method for the vcCNT structure. Initial phases to 5.0 Å were obtained by single isomorphous replacement with anomalous scattering (SIRAS) from a platinum-soaked crystal. Pt sites were found by SHELXD [33], phasing was calculated to 5.0 Å Bragg spacings using SOLVE [34] with the figure of merit 0.34 and the density was modified by solvent flattening using RESOLVE [35]. A partial model was built by manually placing idealized poly-Ala helices into the solvent-flattened electron density map using COOT [36]. After further extensive screening of platinum derivatives, we found a derivative that diffracted to 3.1 Å with significant anomalous signal but was non-isomorphous with any of our native data. Platinum sites were found from an anomalous difference Fourier map using the partial model phases. Combined phases of single anomalous dispersion from the platinum derivative and molecular replacement from the partial model (MR-SAD) were calculated at 4.0 Å with the figure of merit 0.35 and extended to 3.5 Å by solvent flattening with PHASER [37] and RESOLVE [35] using the PHENIX interface [38]. After iterative cycles of manual adjustment of poly-Ala helices and calculation of combined MR-SAD phases, the electron density map was of excellent quality and allowed us to place side chains into the partial model. We also collected native data to 2.4 Å at a long wavelength (1.6 Å) and identified sulfur sites using an anomalous difference Fourier map, which helped guide the model building. After about 70% of manual model building was complete, molecular replacement was done using the partial model against the 2.4 Å native data for further model building and refinement. Structure refinement was done using PHENIX [38]. An anomalous difference Fourier map with the native data to 2.4 Å collected at a long wavelength (1.6 Å) identified six sulfur sites and helped guide the model building. The final model is of good quality with Rwork/Rfree of 19.6/22.8 % and contains residues 2 to 416 (residues 230 to 240 missing), a uridine, 5 DMs (1 full and 4 partial), a sodium ion, and 43 bound waters.
Vesicle Reconstitution and Flux Assay
Protein was reconstituted into lipid vesicles containing 10 mg/mL of 3:1 POPE:POPG at a mass ratio of 1:500 protein:lipid as previously described [39]. Vesicles were reconstituted in the presence of 200 mM KCl, 20 mM HEPES pH 7.4, and 100 mM choline chloride. Reconstituted vesicles were then flash frozen and thawed 3 times then extruded through a 1.0 μm filter using the Avanti Mini-Extruder.For the flux assay, vesicles were diluted 20-fold into buffer containing 200 mM KCl, 20 mM HEPES pH 7.4, and either 100 mM choline chloride or 100 mM NaCl. Flux was initiated by the addition of 2.4 μM [5,6-3H]-uridine and 1 μM valinomycin. All experiments were performed in triplicate at 30 °C. Vesicles were harvested on GF/B glass microfiber filters (Whatman) and counted by scintillation the following day.
Authors: John R Mackey; Lori L Jennings; Marilyn L Clarke; Cheryl L Santos; Laith Dabbagh; Michaela Vsianska; Sheryl L Koski; Robert W Coupland; Stephen A Baldwin; James D Young; Carol E Cass Journal: Clin Cancer Res Date: 2002-01 Impact factor: 12.531
Authors: S R Hamilton; S Y Yao; J C Ingram; D A Hadden; M W Ritzel; M P Gallagher; P J Henderson; C E Cass; J D Young; S A Baldwin Journal: J Biol Chem Date: 2001-05-25 Impact factor: 5.157
Authors: Sylvia Y M Yao; Amy M L Ng; Melissa D Slugoski; Kyla M Smith; Ras Mulinta; Edward Karpinski; Carol E Cass; Stephen A Baldwin; James D Young Journal: J Biol Chem Date: 2007-08-17 Impact factor: 5.157
Authors: J R Mackey; R S Mani; M Selner; D Mowles; J D Young; J A Belt; C R Crawford; C E Cass Journal: Cancer Res Date: 1998-10-01 Impact factor: 12.701
Authors: Airlie J McCoy; Ralf W Grosse-Kunstleve; Paul D Adams; Martyn D Winn; Laurent C Storoni; Randy J Read Journal: J Appl Crystallogr Date: 2007-07-13 Impact factor: 3.304
Authors: Mathieu Coincon; Povilas Uzdavinys; Emmanuel Nji; David L Dotson; Iven Winkelmann; Saba Abdul-Hussein; Alexander D Cameron; Oliver Beckstein; David Drew Journal: Nat Struct Mol Biol Date: 2016-02-01 Impact factor: 15.369
Authors: Natalia Grañé-Boladeras; Christopher M Spring; W J Brad Hanna; Marçal Pastor-Anglada; Imogen R Coe Journal: Cell Mol Life Sci Date: 2016-06-06 Impact factor: 9.261
Authors: Nathan H Joh; Gevorg Grigoryan; Yibing Wu; William F DeGrado Journal: Philos Trans R Soc Lond B Biol Sci Date: 2017-08-05 Impact factor: 6.237