Selenium derivatization of nucleic acids for crystallography.

The high-resolution structure of the DNA (5'-GTGTACA-C-3') with the selenium derivatization at the 2'-position of T2 was determined via MAD and SAD phasing. The selenium-derivatized structure (1.28 A resolution) with the 2'-Se modification in the minor groove is isomorphorous to the native structure (2.0 A). To directly compare with the conventional bromine derivatization, we incorporated bromine into the 5-postion of T4, determined the bromine-derivatized DNA structure at 1.5 A resolution, and found that the local backbone torsion angles and solvent hydration patterns were altered in the structure with the Br incorporation in the major groove. Furthermore, while the native and Br-derivatized DNAs needed over a week to form reasonable-size crystals, we observed that the Se-derivatized DNAs grew crystals overnight with high-diffraction quality, suggesting that the Se derivatization facilitated the crystal formation. In addition, the Se-derivatized DNA sequences crystallized under a broader range of buffer conditions, and generally had a faster crystal growth rate. Our experimental results indicate that the selenium derivatization of DNAs may facilitate the determination of nucleic acid X-ray crystal structures in phasing and high-quality crystal growth. In addition, our results suggest that the Se derivatization can be an alternative to the conventional Br derivatization.

INTRODUCTION

X-ray crystallography is a powerful tool for 3D structure determinations of DNA and RNA molecules and nucleic acid–protein complexes, which provides insight into the structure–function relationships of DNA–drug complexes, functional RNAs, (such as catalytic RNAs and snoRNAs), and RNA–protein and DNA–protein interactions (1–7). Besides the difficulties related to crystallization, however, heavy-atom derivatization for phase determination, a limiting factor in nucleic acid X-ray crystallography, has largely slowed down structural determination of new structures, including DNA–protein complexes. The conventional approaches (8,9) for DNA and RNA derivatization, such as heavy-atom soaking and co-crystallization, have proved to be much more difficult for nucleic acids than for proteins, probably because nucleic acids often lack specific binding sites for metal ions. In the conventional halogen derivatization, bromine (K-edge, 0.920 Å) is used to derivatize DNAs via 5-bromo-deoxyuridine (thymidine analog) and RNAs via 5-bromo-uridine (8,9) in multi-wavelength anomalous dispersion (MAD) phasing. The derivatives containing the bromine derivatization are relatively stable because of formation of vinyl halides. However, it is difficult to introduce bromine to other positions because of the chemical nature of the halogen. For instance, as bromide is a good leaving group, bromine placed on the ribose 2′-position can be displaced by the pyrimidine exo-2-oxygen, the purine 3-nitrogen, or the near-by nucleotide phosphate groups. Though the 8-bromo-purine nucleosides can been easily prepared (10), the 8-Br-substituted purine bases tend to adopt a syn-conformation (11), which perturbs the local structure. Therefore, in the conventional halogen derivatization, the halogens are primarily limited to the 5-position of uracil (sometimes cytosine), which places the halogen in the major groove. As the major groove of A-form dsDNA or dsRNA is narrow and deep, the halogen derivatization may perturb these structures. The halogen derivatization can also cause base-stacking disruption, hydration pattern change and other structural perturbations. In addition, halogen derivatives, such as the bromine derivatives, are light sensitive, and long-time exposure to X-ray or ultraviolet (UV) sources may cause decomposition (12,13). These factors may have caused failures of bromine derivative crystallization, diffraction, and MAD phasing in many cases.

The selenomethionine derivatization method (14–17) has revolutionized the protein X-ray crystallography, and over two-thirds of novel protein structures have been determined recently by the selenomethionine MAD or single-wavelength anomalous dispersion (SAD) phasing (18,19). Indirect derivatization method of RNA using selenomethionine-labeled U1A for phase and structure determination was also successfully demonstrated though it is labor-intensive (20). Recently, we were the first (Huang, Egli and the co-workers) to demonstrate the direct derivatization of DNAs and RNAs with selenium for nucleic acid X-ray crystallography via MAD phasing (21–27). This atom-specific derivatization has great potential in structural and functional studies of nucleic acids. So far, this novel derivatization methodology has been utilized in X-ray crystal structure studies of RNA and DNA molecules by several laboratories (23,24,28–30), including structure analysis and determination of ribozymes and riboswitches. To further facilitate the selenium derivatization of DNAs and RNAs, the enzymatic strategies have also been developed to prepare selenium-derivatized nucleic acids (SeNA) via DNA polymerization and RNA transcription (31,32). As the synthesis and purification of SeNAs are much easier than those of proteins, nucleic acid–protein complexes may also be derivatized by SeNAs instead of their protein counterparts. It appears that the selenium derivatization strategy has great potential in X-ray crystal structure determination of nucleic acids and their protein complexes.

Our UV thermal denaturation and NMR studies of the Se-modified duplexes have indicated the same stability as the non-modified DNAs and RNAs (25–27). In order to further investigate the selenium derivatization strategy, we determined the structure of the DNA (5′-GTGTACAC-3′) at high-resolution (1.28 Å) using the 2′-Se derivatization and also report here the structural comparison between the Se and Br derivatization. It is worth to mention that the Se derivatization is placed on the sugar ring and located in the minor groove of the duplex, while the Br derivatization occurs to the base and is located in the major groove. The structure study reveals that the DNA duplex structure with the Se derivation is identical to the native structure (33) while the Br derivatization causes local structure perturbations in all three bromine-derivatized DNA structures. This report demonstrates that the Se derivatization can be an alternative to the conventional Br derivatization. In addition, we also demonstrated that the Se derivatization may facilitate determination of nucleic acid crystal structures in phasing and high-quality crystal growth.

MATERIALS AND METHODS

Synthesis and incorporation of selenium into oligonucleotides

For the simplicity of the synthesis, 2′-Me-Se-dU (22) is used to replace T2 in the self-complementary DNA (GTGTACAC, Figure 1A). In order to directly compare the Se derivatization and the conventional Br derivatization strategy, we also determined the DNA structure with the Br derivatization at the T4 position (Figure 1B), in the presence and absence of the Se derivatization at the T2 position. The native DNA (5′-GTGTACAC-3′), the Se-DNA (5′-GdUSeGTACAC-3′), the Br-DNA (5′-GTGdUBrACAC-3′), and the Se-Br-DNA (5′-GdUSeGdUBrAC-A-C-3′) were synthesized on solid phase by following the developed strategy (22,26) and were purified by high-performance liquid chromatography (HPLC) using a 21.2 × 250 mm Zorbax, RX-C8 column at a flow rate of 10 ml/min [Buffer A: 20 mM triethylammonium acetate (TEAAc, pH 7.1); Buffer B: 20 mM TEAAc (pH 7.1) in 50% aqueous acetonitrile]. The purified oligonucleotides were analyzed by MS and HPLC to confirm their sequences and purity.

Figure 1

(A and B) Structures of the 2′-Se modification and 5-Br modification of DNA oligonucleotide.

Crystallization

These purified oligonucleotides (1 mM) were first heated to 90°C for 1 min, and the samples were then allowed to cool slowly to 25°C. Though the Se-DNA crystallized in the native buffer (33,34), crystallization conditions were still screened using the Nucleic Acid Mini Screen kit (Hampton Research, ). A total of l μl of the oligonucleotide solutions (1 mM) was typically used in each of these screens. Crystallization was carried out using the hanging drop method by vapor diffusion at 25°C. In the cases of the Se-DNA and the Se-Br-DNA, their crystals with diffraction quality were identified in many identical buffer conditions. Crystals of the Se-DNA (5′-GdUSeGTA-CAC-3′), grown in buffer #7 of the kit [10% v/v MPD, 40 mM sodium cacodylate (pH 6.0), 12 mM spermine tetra-HCI, 80 mM potassium chloride and 20 mM magnessium chloride], were identified to give the highest diffraction resolution, while crystals of the Se-Br-DNA (5′-GdUSeG-dUBrACAC-3′), grown in buffer #9 [10% v/v MPD, 40 mM sodium cacodylate (pH 6.0), 12 mM spermine tetra-HCI, 80 mM sodium chloride, 12 mM potassium chloride and 20 mM magnessium chloride], gave highest resolution. The Br-DNA crystal was grown under the native conditions (33, 34, 40% v/v MPD, 4 mM spermine and 8 mM magnessium chloride). The best Br-DNA crystal was obtained after a week under the lower temperature (15°C). The size of these crystals generating good diffraction is approximately 0.1 × 0.1 × 0.1 mm.

Data collection

30% glycerol or PEG 6000 was used as a cryoprotectant while data collection was taken under the liquid nitrogen stream at 99°K. The Se-DNA and Br-DNA data were collected at beam line X12C, and the Se-Br-DNA datasets were collected at X25, X29 and X12 in NSLS of Brookhaven National Laboratory. A number of crystals were scanned to find the one with strong anomalous scattering at the K-edge absorption of selenium and/or bromine. The distance of the detector to the crystals was set to 90 mm. The chosen wavelengths for selenium and bromine MAD are listed in Tables 1 and 2. The crystals were exposed for 10 or 15 s per image after one degree rotation, and a total of 180 images were taken for each dataset. The additional reference datasets were collected at 1.10 Å wavelength and were used for the final structure refinement. All data were processed using HKL2000 and DENZO/SCALEPACK (35).

Table 1

Summary of data collection and phasing statistics for Se-DNA (1Z7I)

	Selenium K-edge
	Peak	Inflection	Remote
Wavelength Å	0.9790	0.9794	0.9400
Resolution range Å (last shell)	50.0–1.4 (1.49–1.4)	50.0–1.4 (1.49–1.4)	50.0–1.4 (1.49–1.4)
Unique reflections	4478 (441)	4425 (420)	4396 (402)
Completeness %	97.6 (99.1)	96.3 (94.0)	95.2 (90.0)
Rmerge %	6.9 (40.2)	7.2 (44.0)	7.3 (45.0)
<I/σ(I)>	14.0 (3.8)	13.9 (3.9)	13.2 (3.5)
Redundancy	13.0 (10.5)	11.5 (9.7)	10.5 (9.7)
R-Cullis (Friedel)	0.513 (0.330)	0.371 (0.279)	−(0.428)
Phasing power λ3/λι (Friedel)	2.026 (4.345)	3.155 (5.021)	−(3.405)
Figure of merit (Friedel)	0.426 (0.508)	0.502 (0.539)	−(0.380)
Overall figure of merit	0.817

Rmerge = Σ|I − |/ΣI.

Table 2

Data collection and phasing statistics for Se-Br-DNA (2DLJ)

	Selenium K-edge			Bromine K-edge
	Peak	Inflection	Remote	Peak	Inflection	Remote
Wavelength Å	0.9790	0.9794	0.9400	0.9196	0.9199	0.9400
Resolution range Å (last shell)	50.0–1.80 (1.86–1.80)	50.0–1.80 (1.86–1.80)	50.0–1.80 (1.86–1.80)	50.0–1.80 (1.86–1.80)	50.0–1.80 (1.86–1.80)	50.0–1.80 (1.86–1.80)
Unique reflections	2283 (222)	2284 (219)	2288 (220)	2290 (220)	2292 (220)	2288 (220)
Completeness %	99.8 (100.0)	99.8 (100.0)	99.9 (100.0)	99.8 (100.0)	99.9 (100.0)	99.9 (100.0)
Rmerge %	7.8 (36.7)	7.2 (31.8)	6.7 (32.9)	7.9 (39.5)	7.9 (41.0)	7.8 (40.4)
I/σ(I)	11.4 (3.1)	12.2 (3.7)	12.1 (3.7)	11.4 (3.1)	11.1 (2.6)	11.1 (2.7)
Redundancy	12.8 (12.8)	12.9 (13.0)	12.9 (13.0)	12.9 (13.1)	12.9 (13.0)	12.9 (13.0)
R-Cullis (Friedel)	0.709 (0.455)	0.469 (0.359)	−(0.510)	0.686 (0.417)	0.587 (0.407)	−(0.482)
Phasing power (Friedel)	1.221 (3.223)	2.365 (4.043)	−(2.826)	1.205 (4.002)	1.646 (4.053)	−(3.395)
Figure of merit (Friedel)	0.321 (0.429)	0.425 (0.477)	−(0.311)	0.324 (0.465)	0.385 (0.478)	−(0.356)
Overall figure of merit	0.731			0.761

Rmerge = Σ|I − |/ΣI.

Summary of data collection and phasing statistics for Se-DNA (1Z7I)

Rmerge = Σ|I − |/ΣI.

Data collection and phasing statistics for Se-Br-DNA (2DLJ)

Rmerge = Σ|I − |/ΣI.

Structure determination and refinement

The positions of the selenium and bromine atoms were identified by the heavy-atom search scripts provided by CNS (36). Following the refinement of the selenium and bromine positions, the experimental phases were calculated and extended using the MAD or SAD procedures in CNS. The initial phased maps were then improved by the solvent flattening and density modification procedure. The model of single-strand of DNA was manually built into the electron density map using O (37), and then the refinement was carried out. The refinement protocol includes simulated annealing, positional refinement, restrained B-factor refinement and bulk solvent correction. The stereo-chemical topology and geometrical restrain parameters of DNA/RNA (38) have been applied. The topologies and parameters for modified dU with selenium (UMS) and dU with bromine (BRU) were constructed and applied. After several cycles of refinement and the model rebuilding, a number of highly ordered waters were added. Final, the occupancies of selenium and bromine were refined. Cross-validation (39) with a 5–10% test set was monitored during the refinement. The σA-weighted maps (40) of the (2m∣Fo∣ − D∣Fc∣) and the difference (m∣Fo∣ − D∣Fc∣) density maps were computed and used throughout the model building.

RESULTS AND DISCUSSION

Crystallization facilitation by the Se derivatization

We performed the crystallization screening to study the crystallizability of these modified oligonucleotides, and found that the Se-DNA formed quality crystals in many buffers of the Hampton kit overnight (approximate size: 0.1 × 0.1 × 0.1 mm). The Se-Br-DNA also formed quality crystals overnight in the same buffers, though the crystal sizes were smaller (approximately 0.05 × 0.05 × 0.05 mm). However, both the Br-DNA and the native DNA did not grow crystals in any buffers of the kit (24 buffers) over several months. While the native DNA and the Br-DNA crystallize in the native conditions (33,34), it is worth to mention that these two Se-derivatized DNAs can also crystallize in the native conditions.

To further demonstrate this interesting observation, we have synthesized and purified four other native DNAs [(41,42), PDB ID: 395D & 1QPH] and their selenium derivatives (Se-DNAs): 5′-GdUSeACGTAC-3′, 5′-GTAC-GCGdUSeAC-3′, 5′-CCGdUSeACGTACGG-3′, and 5′-GA-CCACGTGGdUSeC-3′. The Se derivatization was incorporated into different regions of these DNA molecules, including the 5′, middle and 3′ regions. These self-complementary DNA molecules were screened using the Hampton kit. High-quality single crystals (average size: 0.1 × 0.1 × 0.1 mm) of these Se-DNA sequences were formed in three days under many different buffer conditions: 5′-GdUSeACGTAC-3′ in buffer #9, #12, #13, #22 and #23, 5′-GTACGCGdUSeAC-3′ in buffer #8, #10, #12, #17, #19 and #24, 5′-CCGdUSeACGTACGG-3′ in buffer #11, #12 and #21, and 5′-GACCACGTGGdUSeC-3′ in buffer #18, #22 and #23. The corresponding native DNA sequences, however, did not crystallize over two weeks in any of these 24 buffers. It seems that it is quite challenging to crystallize the native molecules. In summary, these observations suggest that these Se-derivatized DNAs can crystallize in broader buffer conditions and have better crystallizability than these native and Br-derivatized DNAs.

Our observations indicate that the 2′-selenium derivatization may facilitate the crystal growth, which is probably assisted by the better packing of the Se-derivatized DNA. As the bulky 2′-Se modification may be able to lock the sugar pucker into the 2′-exo conformation, the created rigidity of the Se-derivatized DNAs may facilitate better molecular packing and faster crystal growth. Furthermore, the structure determination of the modified DNAs (the Se-DNA, the Se-Br-DNA and the Br-DNA), grown under different buffer conditions, indicated that the different buffers may not change the crystal packing and the space group of the crystal form.

MAD and SAD phasing

As one Se atom should enable the phase determination for DNAs and RNAs with up to 30 nt (23), one selenium or one bromine atom in the octamer DNA is sufficient for MAD or SAD phasing. As shown in Table 1, we have used the Se-DNA data only up to 1.4 Å for MAD and SAD phasing. The overall figures of merit (FOM) of the initial phases for both Se-DNA MAD and SAD datasets were 0.816 and 0.442, respectively. This produced an excellent electron density map showing a single-strand of DNA after solvent flattening and density modification. The reference data were used as the amplitudes of the structure factors and combined with the phase probabilities from MAD/SAD, and the phases were extended from 1.4 to 1.28 Å for the refinement of Se-DNA. The MAD datasets for the Se-Br-DNA were available to up to 1.8 Å resolutions. The data at the selenium K-edge and the bromine K-edge were processed separately. As shown in Table 2, the overall FOMs of the initial phases were 0.731 and 0.761, respectively from the selenium and bromine K-edge datasets. Both produced excellent electron density maps clearly showing the DNA strands. The experienced phases were extended to 1.5 Å for the refinement of the Se-Br-DNA structures. Likewise, the Br-DNA structure was determined. The data collection summary, phasing statistics and refinement statistics of the Se-DNA, Se-Br-DNA, and Br-DNA are shown in Tables 1–3.

Table 3

Data collection and refinement statistics in Se-DNA, Se-Br-DNA and Br-DNA structures

Structure (PDB ID)	Se-DNA (1Z7I)	Se-Br-DNA (2D LJ)	Se-Br-DNA (2GPX)	Br-DNA (2H05)
Data collection
Space group	P43212	P43212	P43212	P43212
Cell dimensions: a = b, c (Å),	42.076, 23.935	42.369, 23.744	42.252, 24.133	42.048, 24.270
Resolution range Å (last shell)	50.00–1.28 (1.36–1.28)	50.00–1.50 (1.55–1.50)	50.00–1.60 (1.66–1.60)	50.00–1.80 (1.86–1.80)
Unique reflections	6325 (552)	3777 (354)	3161 (297)	2248 (203)
Completeness %	95.3 (70.0)	99.7 (99.2)	99.7 (99.3)	99.4 (95.8)
Rmerge %	8.0 (32.8)	5.1 (27.8)	8.4 (49.9)	4.5 (32.3)
I/σ(I)	14.8 (4.3)	14.4 (4.9)	13.1 (4.2)	12.1 (2.7)
Redundancy	12.4 (10.7)	12.3 (9.1)	12.4 (9.0)	12.1 (7.1)
Refinement
Resolution range Å (last shell)	20.81–1.28 (1.36–1.28)	18.16–1.50 (1.59–1.50)	18.90–1.60 (1.70–1.60)	18.81–1.80 (1.91–1.80)
Rwork %	17.7 (23.7)	18.3 (18.9)	18.7 (21.3)	21.3 (27.3)
Rfree %	20.0 (23.1)	22.8 (22.7)	23.8 (24.4)	24.8 (41.3)
Number of reflections	5275 (552)	3597 (473)	2909 (289)	2036 (247)
Number of atoms
Nucleic Acid (single)	161	160	160	160
Heavy-atoms and ion	1 Se	1 Se, 1 Br	1 Se, 1 Br, 1 Ba2+	1 Br
Water	44	43	31	33
R.m.s.d
Bond length Å	0.012	0.011	0.010	0.009
Bond angle	1.8	1.6	1.8	1.9
Average B-factors Å2
All atoms	15.0	18.8	25.9	32.5
Wilson plot	14.9	18.2	23.6	29.8
Overall anisotropic B-values
B11/B22/B33	0.70/0.70/−1.41	−0.29/−0.29/0.58	0.31/0.31/−0.62	2.20/2.20/−4.41
Bulk solvent correction
Solvent density e/Å3	0.30	0.32	0.34	0.32
B-factors Å2	22.2	34.5	40.05	41.04
Coordinates error (c.-v.) 5Å
Esd. from Luzzatt plot Å	0.15	0.19	0.22	0.23
Esd. from SIGMAA Å	0.10	0.05	0.12	0.20

Rmerge = Σ|I − |/ΣI.

Data collection and refinement statistics in Se-DNA, Se-Br-DNA and Br-DNA structures

Rmerge = Σ|I − |/ΣI.

These Se and/or Br modified DNA structures have the same tetragonal space group P43212, as the native structure (33). These duplex models are generated by applying the symmetry operator (y, x, −z) of P43212. Figure 2 shows the electron density maps with the duplex models of the Se and Br-derivatized DNA structures. The resulting structures of the 2′-Se-DNA (1Z7I), the Se-Br-DNA (2DLJ) and the Br-DNA (2H05) were refined at 1.28, 1.5 and 1.8 Å resolution, respectively. The second Se-Br-DNA structure (2GPX), containing a barium ion, was refined at 1.6 Å resolution; its electron density map (identical to 2DLJ) is not shown. Figure 3 illustrates the details of the electron density maps of the modified dUs in the structures. Figure 3A, B and C represent the structures of the 2′-Se-dU in the Se-DNA, 5-Br-dU in the Se-Br-DNA and 5-Br-dU in the Br-DNA, respectively.

Figure 2

Electron density maps and models of the derivatized DNA duplexes (Stereo view). Red and green balls represent Se and Br, respectively. Contours are at the 1.2σ level. (A) The structure of the Se-DNA (GdUSeGTACAC) at 1.28 Å resolution (1Z7I). (B) the structure of the Se-Br-DNA (GdUSeGdUBrACAC) at 1.5 Å resolution (2DLJ). (C) The structure of the Br-DNA (GTGdUBrACAC) at 1.8 Å resolution (2H05).

Figure 3

Electron density maps and models of the Se and Br modified nucleotides in the derivatized DNAs (Stereo view). Red and green balls represent Se and Br, respectively. Contours are at the 1.2σ level. (A) The structure of the 2′-Se-dU in d(GdUSeGTACAC) at 1.28 Å resolution (1Z7I). (B) The structure of the 5-Br-dU in d(GdUSeGdUBrACAC) at 1.50 Å resolution (2DLJ). (C) The structure of the 5-Br-dU in d(GTGdUBrACAC) at 1.80 Å resolution (2H05).

Electron density maps and models of the derivatized DNA duplexes (Stereo view). Red and green balls represent Se and Br, respectively. Contours are at the 1.2σ level. (A) The structure of the Se-DNA (GdUSeGTACAC) at 1.28 Å resolution (1Z7I). (B) the structure of the Se-Br-DNA (GdUSeGdUBrACAC) at 1.5 Å resolution (2DLJ). (C) The structure of the Br-DNA (GTGdUBrACAC) at 1.8 Å resolution (2H05).

Electron density maps and models of the Se and Br modified nucleotides in the derivatized DNAs (Stereo view). Red and green balls represent Se and Br, respectively. Contours are at the 1.2σ level. (A) The structure of the 2′-Se-dU in d(GdUSeGTACAC) at 1.28 Å resolution (1Z7I). (B) The structure of the 5-Br-dU in d(GdUSeGdUBrACAC) at 1.50 Å resolution (2DLJ). (C) The structure of the 5-Br-dU in d(GTGdUBrACAC) at 1.80 Å resolution (2H05).

Structures of the derivatized DNAs

The Se-DNA structure (1Z7I, in blue, Figure 4A) is superimposable over the native A-form DNA structure (33, 1DNS, in pink), which has the same tetragonal space group (P43212). The overall r.m.s.d. of the Se-DNA over the native is low (0.328 Å) and mainly contributed by the first and the last bases (as shown in Figure 4A), and the main structure is considered as the same as the native, although the spermine molecule bound to the native was not observed in the derivatized structure. The superimposed local dUSe structure of the Se-DNA (1Z7I, in blue) over T2 of the native structure (1DNS, in pink) is shown in Figure 4B. The derivatized local structure of the Se-DNA is also superimposable over that of the second native (34, 1D78, in gray, Figure 4B), which has a different space group (P6122). Therefore, the selenium modification causes no local perturbations when it is compared with both native forms. In addition, superimposing the dUSe in this Se-DNA structure (space group P43212, Figure 4C) with the dUSe in the reported DNA structure (23, GCGTAdUSeA-CGC, 1MA8 in yellow, space group P212121) indicates no difference, suggesting that the Se-derivatized local structure is stable and independent from its sequence environment.

Figure 4

Superimposed comparison of the native and the derivatized DNA structures (red and green balls represent Se and Br, respectively). (A) The comparison of the Se-DNA (GdUSeGTACAC, 1Z7I in blue) and the native (GTGTACAC, 1DNS in pink). (B) The comparison of the Se-DNA dUse (1Z7I in blue) and the two native T2 (1DNS in pink and 1D78 in gray). (C) The comparison of the dUse of the Se-DNA (1Z7I in blue) and the reported DNA (GCGTAdUSeACGC, 1MA8 in yellow). (D) The comparison of the Br-DNA (GTGdUBrACAC, 2H05 in green) and the native (1DNS in pink). (E) The comparison of the Se-Br-DNA (GdUSeGdUBrACAC, 2DLJ in cyan) and the native (1DNS in pink).

Superimposed comparison of the native and the derivatized DNA structures (red and green balls represent Se and Br, respectively). (A) The comparison of the Se-DNA (GdUSeGTACAC, 1Z7I in blue) and the native (GTGTACAC, 1DNS in pink). (B) The comparison of the Se-DNA dUse (1Z7I in blue) and the two native T2 (1DNS in pink and 1D78 in gray). (C) The comparison of the dUse of the Se-DNA (1Z7I in blue) and the reported DNA (GCGTAdUSeACGC, 1MA8 in yellow). (D) The comparison of the Br-DNA (GTGdUBrACAC, 2H05 in green) and the native (1DNS in pink). (E) The comparison of the Se-Br-DNA (GdUSeGdUBrACAC, 2DLJ in cyan) and the native (1DNS in pink).

Figure 4D and E show the superimposed Br-DNA (2H05, in green) and Se-Br-DNA (2DLJ in cyan) over the native (1DNS, in pink), respectively. We observed the local perturbations in both Br-derivatized DNA structures, thought the global structures were very close to the native structure with the same crystal form (1DNS, 33). The overall r.m.s.d. values of the Se-Br-DNA and the Br-DNA over the native are 0.713 and 0.642 Å, respectively. In addition, the comparison of the Br-DNA dUBr with the T4 of the native DNA (1DNS) and the Se-DNA indicated that the Br incorporation could alter the hydration pattern (Figure 5A and B) and the water networking (Figure 6) at the derivatized site, though the sugar pucker of the dUBr was almost not changed from that of the native (Figure 5A). The alterations at the dG3 site of the Br-DNA (Figure 5C and Table 4), including the dG sugar pucker change and ∼107 degree rotation about its C4′–C5′ bond, were also observed. This rotation caused large local perturbations on the dG3, with its C5′ and P deviations by 1.30 and 0.63 Å, respectively. The C2′ of the Br-DNA T2 is also deviated by 0.60 Å (Figure 5D), which also created the flattened sugar pucker as the dG3 (Figure 5C). Unsurprisingly, the same perturbations were also observed with the Se-Br-derivatized DNA (2DLJ) while the perturbations were not observed with the Se-derivatized DNA (1Z7I). These perturbations caused by the Br derivatization may explain why some Br-derivatized DNAs have been difficult to crystallize even after successful crystallization of their natives. Our observation may serve as a reminder in the bromine derivatization.

Figure 5

The superimposed comparison of the native and Br-DNA local structures (green ball represents Br). The green and gray dash lines represent hydrogen and non-hydrogen bonds, respectively. (A) The comparison of the dUBr of the Br-DNA (GTGdUBrACAC, 2H05 in green, and water molecules in blue balls) and the T4 of the native DNA (GTGTACAC, 1DNS in pink). (B) The comparison of the dUBr of the Br-DNA (2H05 in green, and its water molecules in blue balls) and the T4 of the Se-DNA (GdUSeGTACAC, 1Z7I in blue, its water molecules in gray). (C) The comparison of the dG3 of the Br-DNA (in green) and the native (1DNS in pink). (D) The comparison of the T2 of the Br-DNA (in green) and the native (1DNS in pink).

Figure 6

The bromine hydration pattern and the local structures. Water molecules, Se, and Br are in blue, red and green balls, respectively. The green and gray dash lines represent hydrogen and non-hydrogen bonds, respectively. (A) The TGT backbone structure of the native DNA (GTGTACAC, 1D78). (B) The dUSeGT backbone structure of the Se-DNA (GdUSeGTACAC, 1Z7I). (C) The dUSeGdUBr backbone structure of the Se-Br-DNA structure (GdUSeGdUBrACAC, 2GPX). (D) The TGdUBr backbone structure of the Br-DNA (GTGdUBrACAC, 2H05).

Table 4

Summary of the atom deviations [the Br-DNA (2H05) versus the native structure (1DNS)]

Atom (Å)	C2′	C4′	C5′	O5′	P
dUBr versus T4	0.25	0.11	0.24	0.08	0.08
dG3 versus dG3	0.27	0.53	1.30	1.01	0.63
T2 versus T2	0.60	0.08	0.30	0.09	0.22

The superimposed comparison of the native and Br-DNA local structures (green ball represents Br). The green and gray dash lines represent hydrogen and non-hydrogen bonds, respectively. (A) The comparison of the dUBr of the Br-DNA (GTGdUBrACAC, 2H05 in green, and water molecules in blue balls) and the T4 of the native DNA (GTGTACAC, 1DNS in pink). (B) The comparison of the dUBr of the Br-DNA (2H05 in green, and its water molecules in blue balls) and the T4 of the Se-DNA (GdUSeGTACAC, 1Z7I in blue, its water molecules in gray). (C) The comparison of the dG3 of the Br-DNA (in green) and the native (1DNS in pink). (D) The comparison of the T2 of the Br-DNA (in green) and the native (1DNS in pink).

The bromine hydration pattern and the local structures. Water molecules, Se, and Br are in blue, red and green balls, respectively. The green and gray dash lines represent hydrogen and non-hydrogen bonds, respectively. (A) The TGT backbone structure of the native DNA (GTGTACAC, 1D78). (B) The dUSeGT backbone structure of the Se-DNA (GdUSeGTACAC, 1Z7I). (C) The dUSeGdUBr backbone structure of the Se-Br-DNA structure (GdUSeGdUBrACAC, 2GPX). (D) The TGdUBr backbone structure of the Br-DNA (GTGdUBrACAC, 2H05).

Summary of the atom deviations [the Br-DNA (2H05) versus the native structure (1DNS)]

Local perturbations and hydration patterns

The native d(GTGTACAC) structure in tetragonal space group [1DNS, (33)] did not include any water molecules in the refined structure, while the other native structure in hexagonal space group [1D78, (34)], determined at 1.4 Å resolution, has included 37 waters. The TGT local structures of the tetragonal native DNA (GTGTACAC), the hexagonal native, and the Se-DNA are identical (Figures 4A, 6A and B). Instead of comparing with those of the tetragonal native structure due to the absence of the water molecules, we compared the water networking and local backbone structure of the Se-DNA structure with those of the hexagonal native structure. The water molecules were found forming hydrogen bonds to the oxygen atoms of the 5′-phosphate groups of TGT in the native and Se-derivatized structures (Figure 6A and B). In the hexagonal native structure (1D78, Figure 6A), two water molecules (W1 and W3) formed hydrogen bonds (the green dash lines) to the oxygen atoms of the T2 and dG3 5′-phosphate groups. In the Se-DNA structure (1Z7I, Figure 6B), three water molecules (W1, W3 and W4) formed the hydrogen bonds to each other and to the oxygen atoms of the dUSe, dG3 and T4 5′-phosphate groups, forming the same water networking as the native. In the Se-Br-DNA structure (2GPX in Figure 6C), four water molecules (W1, W2, W3 and W4) were found in the region, connecting with each other and the oxygen atoms of the dUSe, dG3 and T4 5′-phosphate groups. These four water molecules formed a different water networking from those of the native and Se-DNA structures. Similar to the Se-Br-DNA structure (Figure 6C), the same water networking (W2, W3 and W4) was found in the Br-DNA (GTGdUBrACAC, 2H05 and Figure 6D). Interestingly, in both the Br-DNA and Se-Br-DNA structures, the W2 water insertion, which can probably facilitate the 5′-phosphate perturbation and the C4′-C5′ bond rotation of the dG3, was observed. This inserted water (W2) formed the hydrogen bonds with W1, the dG3 5′-phosphate oxygen, and W3, which was pushed closer to the dUBr bromine atom in both structures. In both the Br-DNA structure (Figure 6D) and the Se-Br-DNA structure (Figure 6C), we observed that ordered W3 molecules are very close to the bromine atom (3.02 Å and 2.95 Å, respectively), which suggests a hydrogen bond formation. This hydrogen bond may weaken the Br-C5 bond, which probably is one reason of the debromination reported previously (12).

Our experimental results indicate that though it does not cause the change in the sugar pucker of the dUBr, the Br derivatization placed in major groove caused significant conformation changes of the sugar puckers in the dG3 and T2 and the large changes of the hydration pattern and water molecule networking, which may lead to the local perturbation. Unlike the bromine derivatization, the 2′-selenium derivatization is placed in the minor groove, and the hydration alteration in the minor groove was not observed, probably due to disorder of the water molecules near the 2′-position of T2. No ordered water molecules near the 2′-positions in the Se-derivatized and native structures (33,34) have been observed. Probably, the hydration alteration caused by the 2′-selenium modification is insignificant.

CONCLUSIONS

To facilitate X-ray crystal structure studies of nucleic acids, we are in a process of developing selenium derivatization strategies by synthesizing SeNA via atom-specific replacement of oxygen with selenium and studying SeNA structures by comparing them with the structure-known natives and the conventional bromine derivatization. By using both MAD and SAD phasing, we have demonstrated the X-ray crystal structure determination of the 2′-Se-derivtized DNA at high-resolution. The Se-derivatized DNA structure (1.28 Å resolution) is identical to the native structure with the same crystal form (33, 2.0 Å resolution). Besides the demonstration of the formation of these Se-derivatized crystals with high-diffraction quality and in broader buffer conditions, we also demonstrated that these Se-derivatizated DNAs crystallized faster than these native and Br-derivatized DNAs. Furthermore, by comparing with the conventional bromine derivative, which retained the same global structure as the native, we found that the bromine derivatization caused the local perturbations, such as the backbone rotation, the water networking disruption, and the hydrate pattern alteration. Our experimental results indicate that the selenium derivatization may assist the studies of the nucleic acid X-ray crystal structures in two major ways: phasing and crystal growth with high-quality, which leads to structure determination with high-resolution. These experimental results suggest that the selenium derivatization can be an alternative to the conventional bromine derivatization. This Se derivatization strategy via the atom-specific substitution will significantly facilitate X-ray crystal structure studies of nucleic acids and their protein complexes.