Residual Dipolar Coupling Restraints Research Articles

Molecular structure refinement based on residual dipolar couplings using magnetic-field rotational sampling.

A method for structure refinement of molecules based on residual dipolar coupling (RDC) data is proposed. It calculates RDC values using magnetic-field rotational sampling of the rotational degrees of freedom of a molecule in conjunction with molecule-internal configurational sampling. By applying rotational sampling, as is occurring in the experiment, leading to observable RDCs, the method stays close to the experiment. It avoids the use of an alignment tensor and, therefore, the assumptions that the overall rotation of the molecule is decoupled from its internal motions and that the molecule is rigid. Two simple molecules, a relatively rigid and a very flexible cyclo-octane molecule with eight aliphatic side chains containing 24 united atoms, serve as so-called "toy model" test systems. The method demonstrates the influence of molecular flexibility, force-field dominance, and the number of RDC restraints available on the outcome of structure refinement based on RDCs. Magnetic-field rotational sampling is basically equivalent but more efficient than explicitly sampling the rotational degrees of freedom of the molecule. In addition, the performance of the method is less dependent on the number NRDC of measured RDC-values available. The restraining forces bias the overall orientation distribution of the molecule correctly. This study suggests that the information content of RDCs with respect to molecular structure is limited.

A three-dimensional potential of mean force to improve backbone and sidechain hydrogen bond geometry in Xplor-NIH protein structure determination.

We introduce a new hydrogen bonding potential of mean force generated from high-quality crystal structures for use in Xplor-NIH structure calculations. This term applies to hydrogen bonds involving both backbone and sidechain atoms. When used in structure refinement calculations of 10 example protein systems with experimental distance, dihedral and residual dipolar coupling restraints, we demonstrate that the new term has superior performance to the previously developed hydrogen bonding potential of mean force used in Xplor-NIH.

Protein structure prediction using sparse NOE and RDC restraints with Rosetta in CASP13.

Computational methods that produce accurate protein structure models from limited experimental data, for example, from nuclear magnetic resonance (NMR) spectroscopy, hold great potential for biomedical research. The NMR-assisted modeling challenge in CASP13 provided a blind test to explore the capabilities and limitations of current modeling techniques in leveraging NMR data which had high sparsity, ambiguity, and error rate for protein structure prediction. We describe our approach to predict the structure of these proteins leveraging the Rosetta software suite. Protein structure models were predicted de novo using a two-stage protocol. First, low-resolution models were generated with the Rosetta de novo method guided by nonambiguous nuclear Overhauser effect (NOE) contacts and residual dipolar coupling (RDC) restraints. Second, iterative model hybridization and fragment insertion with the Rosetta comparative modeling method was used to refine and regularize models guided by all ambiguous and nonambiguous NOE contacts and RDCs. Nine out of 16 of the Rosetta de novo models had the correct fold (global distance test total score > 45) and in three cases high-resolution models were achieved (root-mean-square deviation < 3.5 å). We also show that a meta-approach applying iterative Rosetta + NMR refinement on server-predicted models which employed non-NMR-contacts and structural templates leads to substantial improvement in model quality. Integrating these data-assisted refinement strategies with innovative non-data-assisted approaches which became possible in CASP13 such as high precision contact prediction will in the near future enable structure determination for large proteins that are outside of the realm of conventional NMR.

Accuracy of MD solvent models in RNA structure refinement assessed via liquid-crystal NMR and spin relaxation data

Molecular dynamics (MD) simulations play an important role in characterizing Ribonucleic Acid (RNA) structure, augmenting information from experimental techniques such as Nuclear Magnetic Resonance (NMR). In this work, we examine the accuracy of structural representation resulting from application of a number of explicit and implicit solvent models and refinement protocols against experimental data ranging from high density of residual dipolar coupling (RDC) restraints to completely unrestrained simulations. For a prototype A-form RNA helix, our results indicate that AMBER RNA force field with either implicit or explicit solvent can produce a realistic dynamic representation of RNA helical structure, accurately cross-validating with respect to a diverse array of NMR observables. In refinement against NMR distance restraints, modern MD force fields are found to be equally adequate, with high fidelity cross-validation to the residual dipolar couplings (RDCs) and residual chemical shift anisotropies (RCSAs), while slightly over-estimating structural order as monitored via NMR relaxation data. With restraints trimmed to encode only for base pairing information, cross-validation quality significantly deteriorates, now exhibiting a pronounced dependence on the choice of the solvent model. This deterioration is found to be partially reversible by increasing planarity restraints on the nucleobase geometry. For completely unrestrained MD simulations, the choice of water model becomes very important, with the best-performing TIP4P-Ew accurately reproducing both the RDC and RCSA data, while closely matching the NMR-derived order parameters. The information provided here will serve as a foundation for MD-based refinement of solution state NMR structures of RNA.

Configuration determination by residual dipolar couplings: accessing the full conformational space by molecular dynamics with tensorial constraints.

Residual dipolar couplings (RDCs) and other residual anisotropic NMR parameters provide valuable structural information of high quality and quantity, bringing detailed structural models of flexible molecules in solution in reach. The corresponding data interpretation so far is directly or indirectly based on the concept of a molecular alignment tensor, which, however, is ill-defined for flexible molecules. The concept is typically applied to a single or a small set of lowest energy structures, ignoring the effect of vibrational averaging. Here, we introduce an entirely different approach based on time averaged molecular dynamics with dipolar couplings as tensorial orientational restraints that can be used to solve structural problems in molecules of any size without the need of introducing an explicit molecular alignment tensor into the computation. RDC restraints are represented by their full 3D interaction tensor in the laboratory frame, for which pseudo forces are calculated using a secular dipolar Hamiltonian as the target. The resulting rotational averaging of each individual tensorial restraint leads to structural ensembles that best fulfil the experimental data. Using one-bond RDCs, the approach has been implemented in the force field procedures of the program COSMOS and extensively tested. A concise theoretical introduction, including the special treatment of force fields for stable and fast MD simulations, as well as applications regarding configurational analyses of small to medium-sized organic molecules with different degrees of flexibility, is given. The observed results are discussed in detail.

Enhanced Sampling of Interdomain Motion Using Map-Restrained Langevin Dynamics and NMR: Application to Pin1

Many signaling proteins consist of globular domains connected by flexible linkers that allow for substantial domain motion. Because these domains often serve as complementary functional modules, the possibility of functionally important domain motions arises. To explore this possibility, we require knowledge of the ensemble of protein conformations sampled by interdomain motion.Measurements of NMR residual dipolar couplings (RDCs) of backbone HN bonds offer a per-residue characterization of interdomain dynamics, as the couplings are sensitive to domain orientation. A challenge in reaching this potential is the need to interpret the RDCs as averages over dynamic ensembles of domain conformations.Here, we address this challenge by introducing an efficient protocol for generating conformational ensembles appropriate for flexible, multi-domain proteins. The protocol uses map-restrained self-guided Langevin dynamics simulations to promote collective, interdomain motion while restraining the internal domain motion to near rigidity. Critically, the simulations retain an all-atom description for facile inclusion of site-specific NMR RDC restraints. The result is the rapid generation of conformational ensembles consistent with the RDC data.We illustrate this protocol on human Pin1, a two-domain peptidyl-prolyl isomerase relevant for cancer and Alzheimer's disease. The results include the ensemble of domain orientations sampled by Pin1, as well as those of a dysfunctional variant, I28A-Pin1. The differences between the ensembles corroborate our previous spin relaxation results that showed weakened interdomain contact in the I28A variant relative to wild type. Our protocol extends our abilities to explore the functional significance of protein domain motions.

Progressive Stereo Locking (PSL): A Residual Dipolar Coupling Based Force Field Method for Determining the Relative Configuration of Natural Products and Other Small Molecules.

Establishing the relative configuration of a bioactive natural product represents the most challenging part in determining its structure. Residual dipolar couplings (RDCs) are sensitive probes of the relative spatial orientation of internuclear vectors. We adapted a force field structure calculation methodology to allow free sampling of both R and S configurations of the stereocenters of interest. The algorithm uses a floating alignment tensor in a simulated annealing protocol to identify the conformations and configurations that best fit experimental RDC and distance restraints (from NOE and J-coupling data). A unique configuration (for rigid molecules) or a very small number of configurations (for less rigid molecules) of the structural models having the lowest chiral angle energies and reasonable magnitudes of the alignment tensor are provided as the best predictions of the unknown configuration. For highly flexible molecules, the progressive locking of their stereocenters into their statistically dominant R or S state dramatically reduces the number of possible relative configurations. The result is verified by checking that the same configuration is obtained by initiating the locking from different regions of the molecule. For all molecules tested having known configurations (with conformations ranging from mostly rigid to highly flexible), the method accurately determined the correct configuration.

Structural Analysis of Multi-Helical RNAs by NMR–SAXS/WAXS: Application to the U4/U6 di-snRNA

NMR and SAXS (small-angle X-ray scattering)/WAXS (wide-angle X-ray scattering) are highly complementary approaches for the analysis of RNA structure in solution. Here we describe an efficient NMR–SAXS/WAXS approach for structural investigation of multi-helical RNAs. We illustrate this approach by determining the overall fold of a 92-nt 3-helix junction from the U4/U6 di-snRNA. The U4/U6 di-snRNA is conserved in eukaryotes and is part of the U4/U6.U5 tri-snRNP, a large ribonucleoprotein complex that comprises a major subunit of the assembled spliceosome. Helical orientations can be determined by X-ray scattering data alone, but the addition of NMR RDC (residual dipolar coupling) restraints improves the structure models. RDCs were measured in two different external alignment media and also by magnetic susceptibility anisotropy. The resulting alignment tensors are collinear, which is a previously noted problem for nucleic acids. Including WAXS data in the calculations produces models with significantly better fits to the scattering data. In solution, the U4/U6 di-snRNA forms a 3-helix junction with a planar Y-shaped structure and has no detectable tertiary interactions. Single-molecule Förster resonance energy transfer data support the observed topology. A comparison with the recently determined cryo-electron microscopy structure of the U4/U6.U5 tri-snRNP illustrates how proteins scaffold the RNA and dramatically alter the geometry of the U4/U6 3-helix junction.

Pressure-induced structural transition of mature HIV-1 protease from a combined NMR/MD simulation approach.

We investigate the pressure-induced structural changes in the mature human immunodeficiency virus type 1 protease dimer, using residual dipolar coupling (RDC) measurements in a weakly oriented solution. (1)DNH RDCs were measured under high-pressure conditions for an inhibitor-free PR and an inhibitor-bound complex, as well as for an inhibitor-free multidrug resistant protease bearing 20 mutations (PR20). While PR20 and the inhibitor-bound PR were little affected by pressure, inhibitor-free PR showed significant differences in the RDCs measured at 600 bar compared with 1 bar. The structural basis of such changes was investigated by MD simulations using the experimental RDC restraints, revealing substantial conformational perturbations, specifically a partial opening of the flaps and the penetration of water molecules into the hydrophobic core of the subunits at high pressure. This study highlights the exquisite sensitivity of RDCs to pressure-induced conformational changes and illustrates how RDCs combined with MD simulations can be used to determine the structural properties of metastable intermediate states on the folding energy landscape.

High-definition NMR structure of PED/PEA-15 death effector domain reveals details of key polar side chain interactions

Death effector domain (DED) proteins constitute a subfamily of the large death domain superfamily that is primarily involved in apoptosis pathways. DED structures have characteristic side chain–side chain interactions among polar residues on the protein surface, forming a network of hydrogen bonds and salt bridges. The polar interaction network is functionally important in promoting protein–protein interactions by maintaining optimal side chain orientations. We have refined the solution DED structure of the PED/PEA-15 protein, a representative member of DED subfamily, using traditional NMR restraints with the addition of residual dipolar coupling (RDC) restraints from two independent alignment media, and employed the explicit solvent refinement protocol. The newly refined DED structure of PED/PEA-15 possesses higher structural quality as indicated by WHAT IF Z-scores, with most significant improvement in the backbone conformation normality quality factor. This higher quality DED structure of PED/PEA-15 leads to the identification of a number of key polar side chain interactions, which are not typically observed in NMR protein structures. The elucidation of polar side chain interactions is a key step towards the understanding of protein–protein interactions involving the death domain superfamily. The NMR structures with extensive details of protein structural features are thereby termed high-definition (HD) NMR structures.

Model of a Six Immunoglobulin-Like Domain Fragment of Filamin A (16–21) Built Using Residual Dipolar Couplings

Filamins are actin-binding proteins that participate in a wide range of cell functions, including cell morphology, locomotion, membrane protein localization, and intracellular signaling. The three filamin isoforms found in humans, filamins A, B, and C, are highly homologous, and their roles are partly complementary. In addition to actin, filamins interact with dozens of other proteins that have roles as membrane receptors and channels, enzymes, signaling intermediates, and transcription factors. Filamins are composed of an N-terminal actin-binding domain and 24 filamin-type immunoglobulin-like domains (FLN) that form tail-to-tail dimers with their C-terminal FLN domain. Many of the filamin interactions including those for glycoprotein Ibα and integrins have been mapped to the region comprising FLN domains 16-21. Traditionally, FLN domains have been viewed as independent folding units, arranged in a linear chain joined with flexible linkers. Recent structural findings have shown that consecutive FLNs form more intricate superstructures. The crystal structure of filamin A domains 19-21 (FLNa19-21) revealed that domains 20 and 21 fold together and that the domain interaction can be autoregulatory. The solution structure of domains 18-19 showed a similar domain interaction, whereas domain pair 16-17 has a completely different domain packing mode. In this study, we characterize the domain organization of the FLNa domain sextet 16-21 using NMR spectroscopy. A structure model of this 60-kDa protein has been built using residual dipolar coupling restraints. RDCs and (15)N relaxation data have been used to characterize interdomain motions.

Paramagnetic relaxation enhancement for the characterization of the conformational heterogeneity in two-domain proteins

Multidomain proteins are often composed of rigid domains that can reorient in solution more or less freely. Calmodulin (CaM) is a two domain protein which can experience a large degree of conformational freedom thanks to a mobile linker connecting the N-terminal and C-terminal domains of the protein. The maximum occurrences (MOs) of the possible protein conformations have been analyzed using the paramagnetic relaxation enhancements (PREs) induced by a gadolinium(III) ion together with the paramagnetic pseudocontact shift and residual dipolar coupling restraints measured in the presence of terbium(III), thulium(III) or dysprosium(III) ions. The results suggest that the PREs provide complementary information useful for improving the description of the conformational heterogeneity of the protein. The data, acquired at 298 K and at 278 K, suggest that compact conformations are disfavoured by decreasing the temperature.

Interaction with Polyglutamine-expanded Huntingtin Alters Cellular Distribution and RNA Processing of Huntingtin Yeast Two-hybrid Protein A (HYPA)

Huntington disease (HD) is an autosomal inherited disorder that causes the deterioration of brain cells. The polyglutamine (polyQ) expansion of huntingtin (Htt) is implicated in the pathogenesis of HD via interaction with an RNA splicing factor, Htt yeast two-hybrid protein A/forming-binding protein 11 (HYPA/FBP11). Besides the pathogenic polyQ expansion, Htt also contains a proline-rich region (PRR) located exactly in the C terminus to the polyQ tract. However, how the polyQ expansion influences the PRR-mediated protein interaction and how this abnormal interaction leads to the biological consequence remain elusive. Our NMR structural analysis indicates that the PRR motif of Htt cooperatively interacts with the tandem WW domains of HYPA through domain chaperoning effect of WW1 on WW2. The polyQ-expanded Htt sequesters HYPA to the cytosolic location and then significantly reduces the efficiency of pre-mRNA splicing. We propose that the toxic gain-of-function of the polyQ-expanded Htt that causes dysfunction of cellular RNA processing contributes to the pathogenesis of HD.

Nuclear magnetic resonance structure of calcium‐binding protein 1 in a Ca2+‐bound closed state: Implications for target recognition

Calcium-binding protein 1 (CaBP1), a neuron-specific member of the calmodulin (CaM) superfamily, regulates the Ca(2+) -dependent activity of inositol 1,4,5-triphosphate receptors (InsP3Rs) and various voltage-gated Ca(2+) channels. Here, we present the NMR structure of full-length CaBP1 with Ca(2+) bound at the first, third, and fourth EF-hands. A total of 1250 nuclear Overhauser effect distance measurements and 70 residual dipolar coupling restraints define the overall main chain structure with a root-mean-squared deviation of 0.54 Å (N-domain) and 0.48 Å (C-domain). The first 18 residues from the N-terminus in CaBP1 (located upstream of the first EF-hand) are structurally disordered and solvent exposed. The Ca(2+) -saturated CaBP1 structure contains two independent domains separated by a flexible central linker similar to that in calmodulin and troponin C. The N-domain structure of CaBP1 contains two EF-hands (EF1 and EF2), both in a closed conformation [interhelical angles = 129° (EF1) and 142° (EF2)]. The C-domain contains EF3 and EF4 in the familiar Ca(2+) -bound open conformation [interhelical angles = 105° (EF3) and 91° (EF4)]. Surprisingly, the N-domain adopts the same closed conformation in the presence or absence of Ca(2+) bound at EF1. The Ca(2+) -bound closed conformation of EF1 is reminiscent of Ca(2+) -bound EF-hands in a closed conformation found in cardiac troponin C and calpain. We propose that the Ca(2+) -bound closed conformation of EF1 in CaBP1 might undergo an induced-fit opening only in the presence of a specific target protein, and thus may help explain the highly specialized target binding by CaBP1.

Solution Structures of Ca2+-CIB1 and Mg2+-CIB1 and Their Interactions with the Platelet Integrin αIIb Cytoplasmic Domain

The calcium- and integrin-binding protein 1 (CIB1) is a ubiquitous Ca(2+)-binding protein and a specific binding partner for the platelet integrin αIIb cytoplasmic domain, which confers the key role of CIB1 in hemostasis. CIB1 is also known to be involved in apoptosis, embryogenesis, and the DNA damage response. In this study, the solution structures of both Ca(2+)-CIB1 and Mg(2+)-CIB1 were determined using solution-state NMR spectroscopy. The methyl groups of Ile, Leu, and Val were selectively protonated to compensate for the loss of protons due to deuteration. The solution structure of Ca(2+)-CIB1 possesses smaller opened EF-hands in its C-domain compared with available crystal structures. Ca(2+)-CIB1 and Mg(2+)-CIB1 have similar structures, but the N-lobe of Mg(2+)-CIB1 is slightly more opened than that of Ca(2+)-CIB1. Additional NMR experiments, such as chemical shift perturbation and methyl group solvent accessibility as measured by a nitroxide surface probe, were carried out to further characterize the structures of Ca(2+)-CIB1 and Mg(2+)-CIB1 as well as their interactions with the integrin αIIb cytoplasmic domain. NMR measurements of backbone amide proton slow motion (microsecond to millisecond) dynamics confirmed that the C-terminal helix of Ca(2+)-CIB1 is displaced upon αIIb binding. The EF-hand III of both Ca(2+)-CIB1 and Mg(2+)-CIB1 was identified to be directly involved in the interaction of CIB1 with αIIb. Together, these data illustrate that CIB1 behaves quite differently from related EF-hand regulatory calcium-binding proteins, such as calmodulin or neuronal calcium sensor proteins.

Solution structure of an alternate conformation of helix27 fromEscherichia coli16S rRNA

Helix (H)27 of 16S ribosomal (r)RNA from Escherichia coli was dubbed the "switch helix" when mutagenesis suggested that two alternative base pair registers may have distinct functional roles in the bacterial ribosome. Although more recent genetic analyses suggest that H27 conformational switching is not required for translation, previous solution studies demonstrated that the isolated E. coli H27 can dynamically convert between the 885 and 888 conformations. Here, we have solved the nuclear magnetic resonance solution structure of a locked 888 conformation. NOE and residual dipolar coupling restraints reveal an architecture that markedly differs from that of the 885 conformation found in crystal structures of the bacterial ribosome. In place of the loop E motif that characterizes the 885 conformer and that the 888 conformer cannot adopt, we find evidence for an asymmetrical A-rich internal loop stabilized by stacking interactions among the unpaired A's. Comparison of the isolated H27 888 solution structure with the 885 crystal structure within the context of the ribosome suggests a difference in overall length of H27 that presents one plausible reason for the absence of H27 conformational switching within the sterically confining ribosome.

Experimentally Restrained Molecular Dynamics Simulations for Characterizing the Open States of Cytochrome P450cam

Residual dipolar couplings (RDCs) were used as restraints in fully solvated molecular dynamics simulations of reduced substrate- and carbonmonoxy-bound cytochrome P450(cam) (CYP101A1), a 414-residue soluble monomeric heme-containing camphor monooxygenase from the soil bacterium Pseudomonas putida. The (1)D(NH) residual dipolar couplings used as restraints were measured in two independent alignment media. A soft annealing protocol was used to heat the starting structures while incorporating the RDC restraints. After production dynamics, structures with the lowest total violation energies for RDC restraints were extracted to identify ensembles of conformers accessible to the enzyme in solution. The simulations result in substrate orientations different from that seen in crystallographic structures and a more open and accessible enzyme active site and largely support previously reported differences between the open and closed states of CYP101A1.

The solution structure of the Mg2+ form of soybean calmodulin isoform 4 reveals unique features of plant calmodulins in resting cells

Soybean calmodulin isoform 4 (sCaM4) is a plant calcium-binding protein, regulating cellular responses to the second messenger Ca(2+). We have found that the metal ion free (apo-) form of sCaM4 possesses a half unfolded structure, with the N-terminal domain unfolded and the C-terminal domain folded. This result was unexpected as the apo-forms of both soybean calmodulin isoform 1 (sCaM1) and mammalian CaM (mCaM) are fully folded. Because of the fact that free Mg(2+) ions are always present at high concentrations in cells (0.5-2 mM), we suggest that Mg(2+) should be bound to sCaM4 in nonactivated cells. CD studies revealed that in the presence of Mg(2+) the initially unfolded N-terminal domain of sCaM4 folds into an alpha-helix-rich structure, similar to the Ca(2+) form. We have used the NMR backbone residual dipolar coupling restraints (1)D(NH), (1)D(C alpha H alpha), and (1)D(C'C alpha) to determine the solution structure of the N-terminal domain of Mg(2+)-sCaM4 (Mg(2+)-sCaM4-NT). Compared with the known structure of Ca(2+)-sCaM4, the structure of the Mg(2+)-sCaM4-NT does not fully open the hydrophobic pocket, which was further confirmed by the use of the fluorescent probe ANS. Tryptophan fluorescence experiments were used to study the interactions between Mg(2+)-sCaM4 and CaM-binding peptides derived from smooth muscle myosin light chain kinase and plant glutamate decarboxylase. These results suggest that Mg(2+)-sCaM4 does not bind to Ca(2+)-CaM target peptides and therefore is functionally similar to apo-mCaM. The Mg(2+)- and apo-structures of the sCaM4-NT provide unique insights into the structure and function of some plant calmodulins in resting cells.

Accurate Solution Structures of Proteins from X-ray Data and a Minimal Set of NMR Data: Calmodulin−Peptide Complexes As Examples

A strategy for the accurate determination of protein solution structures starting from X-ray data and a minimal set of NMR data is proposed and successfully applied to two complexes of calmodulin (CaM) with target peptides not previously described. Its implementation in the present case is based on the use of lanthanide ions as substitutes for calcium in one of the four calcium binding sites of CaM and the collection of pseudocontact shift (pcs) and residual dipolar coupling (rdc) restraints induced by the paramagnetic metals. Starting from the crystal structures, new structural models are calculated that are in excellent agreement with the paramagnetic restraints and differ significantly from the starting crystal structures. In particular, in both complexes, a change in orientation of the first helix of the N-terminal CaM domain and of the whole C-terminal domain is observed. The simultaneous use of paramagnetic pcs and rdc restraints has the following crucial advantages: (i) it allows one to assess the possible presence of interdomain conformational freedom, which cannot be detected if the rdc values are derived from external orienting media; (ii) in the absence of significant conformational freedom, the global orientation tensor can be independently and precisely determined from pcs values, which are less sensitive than rdc values to the presence of local structural inaccuracies, and therefore (iii) the relative rearrangement of a domain or a secondary structure element with respect to the metal-bearing domain can be detected.

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 Methanococcus 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 alpha-helix on top of a 4-stranded twisted beta-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.