Solution NMR Structure of Selenium-binding Protein from Methanococcus vannielii*

Selenium is an important nutrient. The lack of selenium will suppress expression of various enzymes that will lead to cell abnormality and diseases. However, high concentrations of free selenium are toxic to the cell because it adversely affects numerous cell metabolic pathways. In Methanonoccus vannielii, selenium transport in the cell is established by the selenium-binding protein, SeBP. SeBP sequesters selenium during transport, thus regulating the level of free selenium in the cell, and delivers it specifically to the selenophosphate synthase enzyme. In solution, SeBP is an oligomer of 8.8-kDa subunits. It is a symmetric pentamer. The solution structure of SeBP was determined by NMR spectroscopy. Each subunit of SeBP is composed of an α-helix on top of a 4-stranded twisted β-sheet. The stability of the five subunits stems mainly from hydrophobic interactions and supplemented by hydrogen bond interactions. The loop containing Cys59 has been shown to be important for selenium binding, is flexible, and adopts multiple conformations. However, the cysteine accessibility is restricted in the structure, reducing the possibility of the binding of free selenium readily. Therefore, a different selenium precursor or other factors might be needed to facilitate opening of this loop to expose Cys59 for selenium binding.

Selenium is an important nutrient. The lack of selenium will suppress expression of various enzymes that will lead to cell abnormality and diseases. However, high concentrations of free selenium are toxic to the cell because it adversely affects numerous cell metabolic pathways. In Methanonoccus vannielii, selenium transport in the cell is established by the selenium-binding protein, SeBP. SeBP sequesters selenium during transport, thus regulating the level of free selenium in the cell, and delivers it specifically to the selenophosphate synthase enzyme. In solution, SeBP is an oligomer of 8.8-kDa subunits. It is a symmetric pentamer. The solution structure of SeBP was determined by NMR spectroscopy. Each subunit of SeBP is composed of an ␣-helix on top of a 4-stranded twisted ␤-sheet. The stability of the five subunits stems mainly from hydrophobic interactions and supplemented by hydrogen bond interactions. The loop containing Cys 59 has been shown to be important for selenium binding, is flexible, and adopts multiple conformations. However, the cysteine accessibility is restricted in the structure, reducing the possibility of the binding of free selenium readily. Therefore, a different selenium precursor or other factors might be needed to facilitate opening of this loop to expose Cys 59 for selenium binding.
Early studies examining selenium metabolism focused on the toxic potential of this unique element. Symptoms of selenium poisoning were observed in animals that ingested plants containing high levels of selenium (1). The interest in studying the biochemical nature of this trace element has increased dramatically, because it was discovered that low levels of selenium are essential (2). Defects such as liver necrosis, white muscle disease, and certain cardiac and skeletal muscle degenerations are linked to deficiency in selenium in mammals (for review please see Ref. 3). Recently, selenium has been shown to play a central role in redox processes in the cell (4,5). The amount of sele-nium that a cell intakes that is necessary for survival is in a very narrow range (6). Therefore, a sophisticated regulating mechanism is needed to control selenium levels in the cells. Despite its importance, the molecular details of selenium metabolism and transport in living cells are still not completely understood.
One pathway of selenium incorporation is associated with the biosynthesis of proteins containing selenocysteine. This pathway is best characterized in Escherichia coli. The insertion of the selenocysteine amino acid in a specific location in a protein is directed by the in-frame UGA codon and requires the products of four genes selA, selB, selC, and selD in E. coli (7)(8)(9). The activated form of selenium in this pathway is selenium phosphate and is produced by the selD gene product, selenophosphate synthetase (8,10).
Considering the toxicity of selenium, a transport carrier is needed to convert inorganic selenium to the appropriate precursor and deliver it to the proper biosynthetic pathways. The transport mechanism must be highly specific to selenium (11). Furthermore, due to the highly reactive and unstable nature of selenium, the transport mechanism has to preserve the physical state of the selenium during the transport.
Recently a novel new family of proteins that lacks selenocysteine but can bind selenium has been identified in the anaerobic methane-producing organisms. In Methanococcus vannielii, the gene corresponding to a family member of these proteins encodes an 8.8-kDa polypeptide. This protein was purified from the organism as a 42-kDa protein and referred to as the selenium-binding protein (SeBP) 3 (12,13). The recombinant protein was shown to bind selenium when the trace element was added to the growth media of the E. coli. However, the selenium incorporation into the protein was extremely low (13). Further characterization of SeBP shows that its oligomeric state is preserved under extreme conditions of pH, temperature, and denaturant (14). Under denaturing SDS-PAGE conditions, several protein bands corresponding to various oligomers such as dimer, trimer, and tetramer could still be found. This observation was attributed to the high content of hydrophobic amino acids in the SeBP sequence (13).
SeBP contains a single cysteine residue (Cys 59 ) per subunit that is important for selenium binding. Alkylation of this cysteine abolishes selenium incorporation (14). A number of studies on E. coli system have shown that selenodiglutathione (GS-Se-SG) can be used as the precursor from which selenium substrate can be derived (14). Incubation of glyceraldehyde-3phosphate dehydrogenase (GAPDH) from E. coli with this precursor in the presence of selenophosphate synthetase has resulted in selenium incorporation into GAPDH (15). Biochemical studies show that a significant incorporation of selenium can only be observed when SeBP was completely reduced and incubated with the selenodiglutathione as donor in the presence of mild denaturant (1 M guanidine-HCl). The cysteine residue in SeBP appears to be inaccessible. A total of 0.25 equivalents of selenium per SeBP subunit was observed (14). This stoichiometry suggests a 1:4 coordination ratio of the selenium by the cysteines. Despite all of these biochemical data, it is still not clear how SeBP can bind selenium with high selectivity and transport it to the proper biosynthetic pathways. Our study attempts to address, at the atomic level, the high stability of the SeBP oligomer, its mechanism for selenium binding, the accessibility of the cysteine residue, and possible transport of the selenium to the proper enzyme.
First, we determined that SeBP exists as a pentamer in solution. Subsequently, we solved the solution structure of SeBP by Nuclear Magnetic Resonance spectroscopy. The five subunits form a barrel-shaped assembly, where the reactive cysteine is located at the loop close to the core. Solving the SeBP structure is the first step in obtaining a full picture of how this potentially toxic element can be safely managed as an important nutrient for a healthy organism.

MATERIALS AND METHODS
Recombinant Protein-The protein was expressed using E. coli BL21(DE3) harboring pSBP2 (13). The homogeneity of the protein was accomplished using a 4-step purification procedure: anion exchange chromatography (DEAE-Sepharose), hydrophobic interaction chromatography (phenyl-Sepharose), affinity chromatography (anti-SeBP antibody beads), followed by size exclusion chromatography (TSK SW2000). Particularly, anti-SeBP antibody used here was obtained from rabbit by immunizing recombinant SeBP that was purified using the method described previously (13). The anti-SeBP was conjugated to Affi-Gel 15 (Bio-Rad) to prepare anti-SeBP antibody beads. The protein sample from previous steps was mixed with the beads, then the bound proteins were eluted using Immuno-Pure Gentle Ag/Ab elution buffer (Pierce).
Uniformly 15 N-labeled protein was obtained from the culture using minimal medium containing 15 (16). A mixture of 13 C-labeled subunit and unlabeled subunit was prepared as follows: the uniformly 15 N-, 13 C-labeled protein and unlabeled protein were mixed at a 1:1 ratio, denatured using urea, then refolded by removing urea through dialysis. All NMR samples contained 0. Multi-angle Light Scattering-A size exclusion chromatography column (BioSep-SEC-S 2000, Phenomenex) was coupled to a multi-angle detector (DAWN EOS, Wyatt Technology) equipped with a 30 milliwatt linearly polarized GaAs laser for light scattering photometry. The average molecular weight was determined by the Debye fitting method in the ASTRA software package (Wyatt Technology). A refractive index of buffer of 1.330 and a refractive index increment (dn/dc) of 0.184 were used. The samples studied were 0.5 mM SeBP in a volume of 200 l. The buffer was 20 mM Tris/HCl, pH 7.1, and the flow rate used was 0.5 ml/min.
Mass Spectrometry-The mass of the intact form of Sebp was determined by direct infusion at 1 l/min into an electrospray mass spectrometer (Agilent G1969A mass spectrometer with time of flight detector). The mass range from 145 to 200 amu was scanned at 10,000 transients/scan. Voltages were: capillary, 5000 V; fragmentor, 235 V; skimmer, 60 V; and octopole RF, 250 V. The ion source gas temperature was 150°C, drying gas flow was 10 ml/min, and nebulizer gas pressure was 10 psig. It is important to decrease the gas temperature from the usual 350°C to preserve the quaternary structure.
NMR Spectroscopy-Spectra were acquired at 45°C on Bruker 600 or 800 MHz NMR spectrometers, processed using NMRPipe (17), and analyzed with PIPP (18). For the sequential assignment, HN(CO)CACB (19) and HNCACB (20) experiments were carried out using the 15 N-, 13 C-, 2 H-enriched protein. The assignments of 1 H and 13 C resonances were done based on four-dimensional 15 N/ 13 C-edited NOESY (21) and four-dimensional 13 C/ 13 C-edited NOESY (22) experiments using the 15 N-, 13 C-labeled protein. Stereospecific assignments for the methyls of leucines and valines were achieved from 1 H-13 C CT-HSQC experiments (23) using the partially 13 Clabeled protein (16). The spin-lattice relaxation time constants in the rotating frame ( 15 N T 1 ) for backbone amides were calculated from the peak intensities measured using a conventional pulse program (24) with the 15 N-labeled protein.
Structure Calculation-Peak intensities from NOESY experiments were translated into proton-proton distance restraints. TALOS (29) predicted the backbone dihedral angles (, ) from 13 C ␣ , 13 C ␤ , 13 CЈ, 1 H ␣ , and 15 N H chemical shifts. Statistically significant angles were used as structural restraints with at least 10°margins. Generic hydrogen bond distance restraints were employed for ␣-helical regions that were determined based on secondary 13 C ␣ and 13 C ␤ chemical shifts and medium range NOE patterns. Hydrogen bond distance restraints were also employed for ␤-strand regions based on the observation of the 3h J NCЈ connectivity and NOE patterns in addition to the prediction from secondary 13 C ␣ and 13 C ␤ chemical shifts. Residual dipolar couplings for N-H and C-H were calculated from the difference in corresponding scalar couplings measured in the presence and absence of Pf1 phage. Residual dipolar coupling restraints were separated into mobile and rigid regions determined based on the backbone dynamic data.
Structures of SeBP were calculated by distance geometry and simulated annealing protocol (30) with the incorporation of dipolar coupling restraints (31) using the XPLOR-NIH (32). Structure calculations employed 811 proton-proton distance restraints, 40 hydrogen bond restraints, 36 and 36 angle restraints, 63 N-H, 42 C ␣ -H ␣ , and 54 side chain C-H dipolar couplings per subunit. Among 811 proton-proton distance restraints, 250 were intraresidue, 163 were sequential, 99 were mediumrange (2 Ͻ i Ͻ 4), and 299 were long-range (i Ͼ 5). The structure of SeBP is defined by 43 intersubunit proton-proton distance restraints. The final set of 10 structures was chosen from 100 calculated structures on the basis of their overall energies. No distance and dihedral restraints are consistently violated by more than 0.4 Å and 5°, respectively. All molecular illustrations were produced using the program MOLMOL (33).  15 N HSQC spectra shows that the number of cross peaks for backbone amides is roughly 75% of the number of residues of SeBP, indicating that the protein is composed of a symmetric assembly of the 81-resiude monomer. B, peak intensity of the 1 H-15 N HSQC cross peaks is plotted as a function of the residue number. The residues from Ala 54 to Leu 65 were not observed, and nearby residues exhibited low intensities. C, enlarged view of 1 H-15 N cross peaks corresponding to residues Gly 48 -Ile 53 and Gly 66 -Lys 72 . The cross peaks for residues next to the non-observable region, such as Leu 52 , Ile 53 , Gly 66 , and Tyr 67 , show a doublet or a triplet, indicating that this region undergoes a conformational exchange.  Denaturation Study of SeBP-The mutation of Ile-9 to Ala was made using the QuikChangeII site-directed mutagenesis kit (Stratagene), and the mutant protein I9A was produced and purified as a wild-type protein. The wild-type and mutant SeBP were incubated in 50 mM sodium phosphate, pH 7.0, 1 mM EDTA, 100 mM NaCl, 10 mM DTT for 30 min at room temperature in the anaerobic chamber. Reduced proteins were dialyzed against the same buffer without DTT. Differential scanning calorimetry (DSC) was carried out with a N-DSC II (Calorimetry Sciences) at a heating rate of 1°C/min. Sample proteins at 0.8 ϳ 1 mg/ml were used for each scanning. The data analysis was performed using CpCalc V.2.1 (Calorimetry Sciences).

RESULTS
Preliminary NMR Studies-A total of 64 resonances corresponding to backbone H-N could be observed in the 1 H N -15 N HSQC spectrum of the 81-residue polypeptide (Fig.  1A). Resonances of the N-terminal two residues could not be observed, along with the peaks corresponding to Asp 13 -Ile 15 and Phe 55 -Leu 65 . The fact that we could not observe these resonances suggests that some parts of SeBP are undergoing conformational exchange and contain some structural variations. We eliminated possible solvent exchange contribution to the missing resonances by pH titration experiments (pH 6.3-7.5) that showed no significant changes in the NMR spectrum as well as resonance line widths of SeBP (data not shown). The freshly prepared NMR sample did not contain any disulfide bonds. The sample that was left for several weeks after preparation showed the same HSQC spectrum as the freshly prepared sample. The sample in which DTT was added also showed the same HSQC spectrum. Together with the previous finding that Cys 59 is not easily modified (14), we can eliminate the possibility of missing resonances due to the oxidative state of Cys 59 . Upon closer inspection of the HSQC spectrum, multiple resonances with decreased intensity can be observed for residues preceding and trailing the region corresponding to Phe 55 -Gly 65 (Fig. 1, B and C). This indicates multiple con-  formations that are sampled at different rates for different residues. It should be noted, however, that we were able to assign some side-chain atoms for Ile 15 , Pro 16 , Phe 55 , Ile 57 , and Leu 65 , based on the analysis of NOESY spectra.
Oligomeric State of SeBP-Based on earlier studies, SeBP was proposed to be a homotetramer (13), (14). The 4:1 stoichiome-try of SeBP to selenium would suggest that a tetramer is indeed plausible. While the molecular mass of the subunit, as derived from MALDI-TOF, was 8,804 Da, the estimated size of the protein obtained from native PAGE as well as size exclusion chromatography was 42 kDa. Interestingly the denaturing SDS-PAGE showed the highest SeBP band migrated as a 33-kDa protein. Therefore it is natural to deduce from the available molecular size data that SeBP was tetrameric. However, for determining structures by NMR, it is crucial for the exact state of the oligomers to be determined.
The size of the SeBP oligomer will be directly manifested in its NMR relaxation. The 15 N transverse (T 2 ) relaxation time is a strong measure of the overall correlation time or the overall tumbling rate of the molecule, which in turn reflects its size. The 15 N T 2 values measured at 42°C and 800 MHz 1 H frequency are shown in Fig. 2. The average 15 N T 2 for the nonflexible portions of the protein is 47 ms. Based on this average value and using the Debye-Stokes-Einstein relationship, assuming a spherical oligomer, the size of SeBP is estimated to be 45 kDa. This is in the range of the previously estimated size of the protein. In addition, a light scattering experiment was performed on the SeBP. The result is shown in Fig. 3. A monodispersed distribution of molecular weights clustered around 43 kDa was found. This is consistent with the NMR relaxation data. What is unequivocally convincing data of the oligomeric size of SeBP was obtained from the mass spectrometry measurement. Under mildly ionizing conditions, the oligomeric state of SeBP could be preserved during the ion spray stage of the mass spectrometer, albeit as a fraction of the population of the full sample. The measured mass of SeBP under this condition is 44,116.9 Da. All of these data taken together indicate that SeBP is pentameric.
Structure of a Subunit of SeBP-NMR data were collected to provide a total of 811 NOEs per subunit for the determination of the SeBP structure. Dihedral restraints were obtained based on TALOS analysis of the chemical shifts of SeBP. The dihedral and NOE data indicated that SeBP contains a single ␣-helix and four ␤-strands (Fig. 4A). In addition, dipolar couplings were measured for 1 H-15 N and 1 H-13 C for the protein aligned in Pf1 phage. An NMR experiment to measure scalar coupling through the hydrogen bond was performed on SeBP to establish the hydrogen bond network present in the protein. Because of the size of SeBP; thus short T 2 values and limited signal to noise, only 4 hydrogen bonds could be directly observed. This however provides important information on how the ␤-strands are arranged. A twisted ␤-sheet is formed by two strands, ␤2 and ␤4, in an antiparallel fashion as the through-hydrogen bond 3h J NCЈ connectivities were observed between Tyr 21 -H N and Lys 72 -CЈ, between  Leu 23 -H N and Ile 70 -CЈ, between Ala 26 -H N and Gly 68 -CЈ, as well as between Ile 70 -H N and Gly 24 -CЈ, in addition to the NOE connectivities. In contrast, the ␤1and ␤3-strands form a short, parallel ␤-sheet that is defined mainly by NOEs. The structure of a single subunit of SeBP is illustrated in Fig.  4, B and C. Each subunit of SeBP is a twisted ␤-sheet containing four strands with the single ␣-helix capping one side of the sheet. Strand ␤1 consists of residues Ile 8 -Thr 10 and is followed by a loop connecting it to strand ␤2 (Leu 20 -Ser 27 ). Some resonances of the residues (Asp 13 and Glu 14 ) in the loop connecting two strands, ␤1 and ␤2, were not observed. Residues Thr 28 -Asn 32 form the loop connecting strand ␤2 to the single helix in SeBP. The helix is composed of residues Val 33 -Ala 46 . This helix is followed by a short loop (Lys 47 -Met 50 ) that connects it to strand ␤3 (Gly 51 -Ile 53 ) in the structure. A long loop (Ala 54 -Gly 66 ) containing the single cysteine residue (Cys 59 ) in SeBP connects strand ␤3 to strand ␤4 (Tyr 67 -Lys 72 ). Almost all 1 H N -15 N resonances (Ala 54 -Leu 65 ) belonging to residues in this long loop were not observed. As indicated by the increased 15 N transverse (T 2 ) relaxation times at both termini (Fig. 2), the N-terminal (residues Met 1 -Phe 7 ) and C-terminal loops (residues Glu 75 -Ala 81 ) are flexible. The other regions showed welldefined structures except the ␤3-␤4 loop that showed structural dispersion. This is in agreement with the HSQC spectra that suggested possible conformational variations.
Overall Structure of the SeBP-The SeBP structures were calculated from 5 sets of a completely extended starting structure to a pentamer of a barrel-shaped assembly. Several of the long range NOEs obtained from four-dimensional 13 C/ 13 C-edited NOESY experiment using uniformly 15 N-, 13 C-labeled protein sample, such as: Ile 8 -Ser 22 , Ile 9 -Ile 25 , Ile 25 -Phe 55 , and Ile 57 -Leu 65 , were also detected in three-dimensional 13 C-edited, 12 C-filtered NOESY experiments using the mixture of 13 C-labeled subunit and unlabeled subunit, and therefore included in the calculation as intersubunit contacts. The structures converged to a single conformation (Fig. 5) with the lowest overall energies and acceptable precision as shown in Table 1. The N and C termini of each subunit are located at one side of the pentamer (designated as top), while the loop containing Cys 59 is at the opposite side (designated as bottom). The top of the SeBP seems to have a slight opening, suggesting a possible solvent entry, while the bottom is closed. The center of SeBP consists primarily of residues of which NMR resonances were not observed, such as Asp 13 -Glu 14 and Phe 55 -Cys 59 .
NMR resonances from several residues in the loop containing the single cysteine in SeBP could not be observed. However, the conformation of this loop can still be defined, although not as well as the rest of the protein. Despite the ␤3-␤4 loop being long (12 residues), it is mostly positioned inside the barrel, which has a scaffold composed by the ␣-helix and ␤-strands. This seems to restrict the loop enough from being completely disordered. In addition, several NOE contacts involving side chains of Phe 55 , Ile 57 , and Leu 65 to other parts of the protein assist in defining the possible conformations of this loop. It is still important to point out that this loop does not adopt a single conformation. This is reflected in the high r.m.s.d. values in the calculated NMR structures for this loop.
One important question is to understand the high stability of SeBP. The primary forces that stabilize the subunit complex in SeBP are hydrophobic interactions between the subunits. For example, Ile 25 showed intersubunit NOEs to Ile 9 and Val 37 (Fig.  5C). NOEs between Ile 8 and Ser 22 were also observed. In the ␤1-strand, Ile 8 -O and Thr 10 -H N form hydrogen bonds with the H N and O of Leu 52 from the ␤3-strand to form a parallel ␤-sheet, while backbone H N and O atoms of Ile 9 are free, and facing toward Ser 22 . There is a high likelihood of a hydrogen bond formation between the backbone H N atom of Ile 9 and the side chain O ␥ atom of Ser 22 . Such a hydrogen bond would enhance the interactions between subunits and increases the stability of the SeBP.
The electrostatic surfaces of the SeBP are shown in Fig. 6. The outside of the barrel shows bands of positive and negative surfaces. The inner side from the top shows a concentrated negatively charged region. In the bottom of the SeBP, where the ␤3-␤4 loop has undefined conformations, a mix of positive and negative charges can be found.
Stability of SeBP-To support our NMR study that showed the stability of SeBP is due to hydrophobic interactions among subunits, we mutated a residue at the interface between the subunits, and investigated its thermal denaturation. We selected Ile 9 to be mutated because it resides at the interface, and showed intersubunit NOEs to Gly 24 and Ile 25 , being involved in hydrophobic interactions. The structural integrity of I9A mutant protein was confirmed to be the same as that of the wild-type protein by circular dichroism (data not shown). The thermal stability of the I9A mutant protein, however, was found to be very different from that of the wild-type protein. As shown in Fig. 7, the wild-type SeBP shows a thermal transition at 88.0°C while the I9A mutant shows a much earlier transition at 70.8°C in DSC. The mutant protein maintains less interaction forces compared with the wild-type protein, and thus is more heat-sensitive.
Structural Similarity to Other Proteins-The combination of four strands and a helix in a protein structure obviously is not uncommon. However, their arrangement in SeBP is quite unique. The twisted ␤-sheet is formed by aligning the ␤-strands 1, 3, 4, and 2 in order. Strands ␤1, ␤2, and ␤3 are parallel and strand ␤4 is antiparallel with the rest. A search of fold similarity was performed using the DALI program (34) against a representative set of structures in the PDB data base. The DALI search revealed that the pentameric assembly of SeBP is similar to that of a pentameric protein found in other bacteria. The list from the DALI search includes 1Y2I from Shigella flexneri and 2GTC from Bacillus cereus. These two proteins appear to be homologues of each other with 52% sequence identity, but do not show significant similarity to SeBP. The crystal structures 1Y2I and 2GTC were determined as part of a structural genomics effort, and their function is yet unknown. Their structural comparisons are illustrated in Fig. 7. When compared with SeBP, two-thirds of the 1Y2I and 2GTC structures exhibit similarity to the SeBP. The 1Y2I protein has one cysteine in its sequence, but it is located in the middle of strand ␤2 (indicated by a ball in Fig. 8). There is no cysteine in the sequence of 2GTC. The ␤3-strands of 1Y2I and 2GTC are elongated and contain a ␤-sheet along with ␤2and ␤4-strands in an antiparallel fashion. Instead, the corresponding region of SeBP is a loop, of which the backbone amide resonances are unobservable, and expected not to adopt a fixed ␤-sheet conformation. The overall scaffold of SeBP is similar to those of 1Y2I and 2GTC, suggesting perhaps a general architecture for highly stable proteins where the bottom part of it serves as a platform to define its unique biological function. The coordinates for the 10 lowest energy structures of SeBP have been deposited in the PDB with the accession number 2JZ7.

DISCUSSION
Selenium-binding protein of M. vannielii exists as a pentamer in solution. The interaction between the subunits is stabilized by hydrophobic interaction as well as hydrogen bonds. The strong subunit associations could be related to the protein function to transport selenium in the cell. A stable pentamer provides an excellent framework for safely sequestering the selenium from other factors in various cellular environments. Breakdown in the SeBP during the selenium transport will have adverse consequences to the organism because of the high toxicity of this element. The low number of dynamic and exposed loops connecting various secondary structure elements in SeBP contributes to its stability as well. This feature in conjunction with the amino acid compositions of the loops can reduce the potential of degradation of SeBP by proteolysis.
The single cysteine (Cys 59 ) located in the ␤3-␤4 loop is the key residue for selenium reactivity. The loop does not adopt a defined conformation. NMR data indicated that the loop undergoes conformational exchange and conformational variations. However, the SeBP structure revealed that the motion of the loop is limited within the barrel. At the outset, it was anticipated that the Cys residues would be accessible to selenium derivatives, but it was found in the previous study (14) that the Cys residues were not completely exposed without a denatur- ant. Our structure supports the fact that not all Cys residues are freely accessible. This partial accessibility may be related to the regulated reaction of SeBP with the selenium derivatives. The high specificity for selenium must come from the local environment around the Cys 59 . Several charged (Asp 31 , Asp 61 , Lys 63 ) as well as some hydrophobic residues (Thr 58 , Ala 60 , Phe 64 , Leu 65 ) are surrounding the Cys 59 . Even though the conformation of this loop is not well determined, therefore the relative placement of these residues are not fixed with respect to Cys 59 , they do offer a possibility to provide specificity for the substrate. In fact, the multiple conformations adopted by this loop may be necessary to provide selectivity in selenium binding.
Access to the Cys 59 can be achieved from either side of the SeBP core. Entrance from the top of the core, where a small opening is already present, has to involve interaction of the substrate with residues lining the core, which are mostly hydrophobic. Access from the bottom of the core will have to involve the opening of the Cys 59 loop to expose the cysteine. Either mechanisms will affect the on and off rate of substrate binding. This will also add to the stability of the bound selenium. The need for different selenium precursor or other factors to induce the opening of the Cys 59 loop to expose the cysteine to bind selenium cannot be ruled out. On the other hand, it is likely that only interaction with the selenophosphate synthase can induce a conformational change in the SeBP that exposes the bound selenium for the enzyme. This provides the selectivity in delivery of the selenium for this particular enzyme.
Selenium can be found in a variety of compounds, where the most common oxidation numbers of selenium are 6, 4, 2, -1, and -2, including selenophosphate (H 2 SePO 3 Ϫ ) and selenodiglutathione (GS-Se-SG). The possible coordination states of selenium as well as the measured stoichiometry between selenium and SeBP (14) do not seem to correspond to the five cysteines found in the pentameric SeBP. It is still not clear what important role one additional cysteine plays in the binding of selenium. Naturally a structure of SeBP with bound selenium will address this question and complete the picture of molecular interactions that define the function of this protein. As shown previously, the binding of selenium to SeBP in vitro is not trivial (14). This has to be carried out under partially denaturing and anaerobic conditions with the SeBP fully reduced, and it involves selenodiglutathione as the precursor of the selenium. So far the yield of bound selenium to SeBP under our sample conditions is only a few percent, eliminating the possibility of structurally studying this complex. It is clear that future structural studies are needed to address some of the questions regarding the selenium biosynthetic pathways in this organism. Nevertheless, our current SeBP structure affords the first glimpse of how this protein can manage the potentially toxic selenium in the cell and deliver it as an important nutrient for the organism.