A high-resolution structure of a DNA-chromomycin-Co(II) complex determined from pseudocontact shifts in nuclear magnetic resonance.

BACKGROUND: The drug chromomycin-A(3) binds to the minor groove of DNA and requires a divalent metal ion for complex formation. (1)H, (31)P and (13)C pseudocontact shifts occurring in the presence of a tightly bound divalent cobalt ion in the complex between d(TTGGCCAA)(2) and chromomycin-A(3) have been used to determine the structure of the complex. The accuracy of the structure was verified by validation with nuclear Overhauser enhancements (NOEs) and J-coupling constants not used in the structure calculation. RESULTS: The final structure was determined to 0.7 A resolution. The structure was compared with a structure obtained in an earlier study using NOEs, in order to assess the accuracy of NOEs in giving global structural information for a DNA complex. Although some basic features of the structures agreed, they differed substantially in the fine structural details and in the DNA axis curvature generated by the drug. The distortion of base-pair planarity that was observed in the NOE structure was not seen in our structure. Differences in drug orientation and hydrogen bonding also occurred. The curvature and elongation of the DNA that was obtained previously was not found to occur in our study. CONCLUSIONS: The use of pseudocontact shifts has enabled us to obtain a high-precision global structure of the chromomycin-DNA complex, which provides an accurate template on which to consider targeting minor groove binding drugs. The effect of such binding is not propagated far along the helix but is restricted to a local kink in the axis that reverts to its original direction within four base pairs.

Interactions of human nucleotide excision repair protein XPA with DNA and RPA70 Delta C327: chemical shift mapping and 15N NMR relaxation studies.

Human XPA is an essential component in the multienzyme nucleotide excision repair (NER) pathway. The solution structure of the minimal DNA binding domain of XPA (XPA-MBD: M98-F219) was recently determined [Buchko et al. (1998) Nucleic Acids Res. 26, 2779-2788, Ikegami et al. (1998) Nat. Struct. Biol. 5, 701-706] and shown to consist of a compact zinc-binding core and a loop-rich C-terminal subdomain connected by a linker sequence. Here, the solution structure of XPA-MBD was further refined using an entirely new class of restraints based on pseudocontact shifts measured in cobalt-substituted XPA-MBD. Using this structure, the surface of XPA-MBD which interacts with DNA and a fragment of the largest subunit of replication protein A (RPA70 Delta C327: M1-Y326) was determined using chemical shift mapping. DNA binding in XPA-MBD was highly localized in the loop-rich subdomain for DNA with or without a lesion [dihydrothymidine (dhT) or 6-4-thymidine-cytidine (64TC)], or with DNA in single- or double-stranded form, indicating that the character of the lesion itself is not the driving force for XPA binding DNA. RPA70 Delta C327 was found to contact regions in both the zinc-binding and loop-rich subdomains. Some overlap of the DNA and RPA70 Delta C327 binding regions was observed in the loop-rich subdomain, indicating a possible cooperative DNA-binding mode between XPA and RPA70 Delta C327. To complement the chemical shift mapping data, the backbone dynamics of free XPA-MBD and XPA-MBD bound to DNA oligomers containing dhT or 64TC lesions were investigated using 15N NMR relaxation data. The dynamic analyses for the XPA-MBD complexes with DNA revealed localized increases and decreases in S2 and an increase in the global correlation time. Regions of XPA-MBD with the largest increases in S2 overlapped regions having the largest chemical shifts changes upon binding DNA, indicating that the loop-rich subdomain becomes more rigid upon binding DNA. Interestingly, S2 decreased for some residues in the zinc-binding core upon DNA association, indicating a possible concerted structural rearrangement on binding DNA.

Structure determination by restrained molecular dynamics using NMR pseudocontact shifts as experimentally determined constraints.

The structure of a DNA octamer d(TTGGCCAA)(2) complexed to chromomycin-A(3) and a single divalent cobalt ion has been solved by using the pseudocontact shifts due to the unpaired electrons on the cobalt. A protocol was developed and critically evaluated for using the pseudocontact shifts in structure determination. The pseudocontact shifts were input as experimental restraints in molecular dynamics simulations with or without NOE constraints. Both the magnitude and orientation of the susceptibility anisotropy tensor required for the shift calculations were determined during the simulations by iterative refinement. The pseudocontact shifts could be used to define the structure to a very high precision and accuracy compared with a corresponding NOE-determined structure. Convergence was obtained from different starting structures and tensors. A structure determination using both NOE's and pseudocontact shifts revealed a general agreement between the two data sets. However, some evidence for a discrepancy between NOE's and pseudocontact shifts was observed in the backbone and terminal base pairs of the DNA. Violations in shift or NOE restraints remaining in the final structures were examined and may be a reflection of motional averaging of the constraints and evidence for flexibility. This work demonstrates that pseudocontact shifts are a powerful tool for NMR structure determination.

Nuclear magnetic resonance characterization of a paramagnetic DNA-drug complex with high spin cobalt; assignment of the 1H and 31P NMR spectra, and determination of electronic, spectroscopic and molecular properties.

The proton NMR spectrum of the ternary complex between the octamer duplex d(TTGGCCAA)2, two molecules of the drug chromomycin-A3, and a divalent cobalt ion has been assigned. Assignment procedures used standard two-dimensional techniques and relied upon the expected NOE contacts observed in the equivalent diamagnetic complex containing zinc. The magnetic susceptibility tensor for the cobalt was determined and used to calculate shifts for all nuclei, aiding in the assignment process and verification. Relaxation, susceptibility, temperature and field dependence studies of the paramagnetic spectrum enabled determination of electronic properties of the octahedral cobalt complex. The electronic relaxation tau(s) was determined to be 2.5 +/- 1.5 ps; the effective isotropic g value was found to be 2.6 +/- 0.2, indicating strong spin-orbit coupling. The magnetic susceptibility tensor was determined to be chi(xx) = 8.9 x 10(-3) cm3/mol, chi(yy) = 9.5 x 10(-3) cm3/mol, chi(zz) = 12.8 * 10(-3) cm3/mol. A tentative rotational correlation time of 8 ns was obtained for the complex. Both macroscopic and microscopic susceptibility measurements revealed deviations from Curie behavior over the temperature range accessible in the study. Non-selective relaxation rates were found to be inaccurate for defining distances from the metal center. However, pseudocontact shifts could be calculated with high accuracy using the dipolar shift equation. Isotropic hyperfine shifts were factored into contact and dipolar terms, revealing that the dipolar shift predominates and that contact shifts are relatively small.

Ligand-induced perturbations in Urtica dioica agglutinin.

The binding of the trisaccharide, N,N',N"-triacetylchitotriose, to Urtica dioica agglutinin (UDA) was investigated using 1H NMR spectroscopy. UDA is a small antiviral plant lectin containing two homologous 43-amino acid domains. Carbohydrate-induced pertubations occur in one domain of UDA at trisaccharide concentrations below equimolar. Residues in the second domain are shifted at higher carbohydrate concentrations. This data confirms the presence of two binding sites of non-identical affinities per UDA monomer. Qualitative analysis of the 2D NOESY spectra indicates that UDA contains two short stretches of antiparallel beta-sheet. The 1H resonance assignments for both antiparallel beta-sheet sequences have been completed and there is one beta-stretch per domain. A number of these beta-sheet residues are perturbed in the presence of carbohydrate.

Protein structure refinement based on paramagnetic NMR shifts: applications to wild-type and mutant forms of cytochrome c.

A new approach to NMR solution structure refinement is introduced that uses paramagnetic effects on nuclear chemical shifts as constraints in energy minimization or molecular dynamics calculations. Chemical shift differences between oxidized and reduced forms of horse cytochrome c for more than 300 protons were used as constraints to refine the structure of the wild-type protein in solution and to define the structural changes induced by a Leu 94 to Val mutation. A single round of constrained minimization, using the crystal structure as the starting point, converged to a low-energy structure with an RMS deviation between calculated and observed pseudo-contact shifts of 0.045 ppm, 7.5-fold lower than the starting structure. At the same time, the procedure provided stereospecific assignments for more than 45 pairs of methylene protons and methyl groups. Structural changes caused by the mutation were determined to a precision of better than 0.3 A. Structure determination based on dipolar paramagnetic (pseudocontact) shifts is applicable to molecules containing anisotropic paramagnetic centers with short electronic relaxation times, including numerous naturally occurring metalloproteins, as well as proteins or nucleic acids to which a paramagnetic metal ion or ligand may be attached. The long range of paramagnetic shift effects (up to 20 A from the iron in the case of cytochrome c) provides global structural constraints, which, in conjunction with conventional NMR distance and dihedral angle constraints, will enhance the precision of NMR solution structure determination.

Two-dimensional COSY and two-dimensional NOE spectroscopy of d(AC)4.d(GT)4: extraction of structural constraints.

Pure absorption phase, proton two-dimensional nuclear Overhauser effect (2D NOE) and double-quantum-filtered COSY (DQF-COSY) spectra were recorded for d(AC)4.d(GT)4. A full proton resonance assignment was made, except for the 5' and 5" protons. A new semiautomatic method for improved quantitation of 2D NOE peak intensities was developed, and its limitations and usefulness were examined. With this new method, 2D NOE intensity sets at several mixing times were obtained. Simulations of the 1'2', 1'2", and 2'3' DQF-COSY cross-peaks were compared with experimental data, establishing an alternating sugar pucker for the alternating purine-pyrimidine sequence. Scalar coupling constants for the sugar ring protons, derived from the fitting of the simulated spectra, are reported. Complete relaxation matrix analysis of the 2D NOE spectrum verified this alternating structure for all NOE interactions between nonexchangeable protons. Both the DQF-COSY and the 2D NOE results qualitatively indicate that the structure of d(AC)4.d(GT)4 resembles wrinkled D-DNA in aqueous solution.

Solution structure studies of d(AC)4.d(GT)4 via restrained molecular dynamics simulations with NMR constraints derived from two-dimensional NOE and double-quantum-filtered COSY experiments.

The structure of d(AC)4.d(GT)4 is investigated by constrained molecular dynamics simulations. The constraints include proton pair distances derived from 2D NOE intensities by using the iterative relaxation matrix analysis algorithm MARDIGRAS and sugar pucker phases and amplitudes derived from double-quantum-filtered COSY spectra. Molecular dynamics runs on simulated intensity and distance sets as well as the experimental data were carried out to determine the effects of starting structure, distance constraint derivation, energy functions, and experimental errors on the end result. It was found that structural details could not be elucidated within about 1.5-A overall atomic deviation. This limitation is due in part to the accuracy of the experimental data but, more importantly, is attributable to the quantity of experimental constraints available and to imperfections in the force field utilized in the molecular dynamics calculations. Within the limits of the method, some structural characteristics of d(AC)4.d(GT)4 could be elucidated.

In vivo 19F-NMR of 5-fluorouracil incorporation into RNA and metabolites in Escherichia coli cells.

19F resonances from RNA with 5-fluorouracil incorporated could be observed in intact Escherichia coli cells, as well as in tRNA isolated from the cells. 19F-NMR signals from the metabolic breakdown products of the fluorinated RNA were also detected in vivo. By observing the 19F-NMR spectrum, variations in the metabolic disposition of administered 5-fluorouracil could be monitored as a function of time and be compared when the cells were deprived of oxygen and other nutrients, subjected to ethidium bromide treatment, or grown in the presence of mitomycin C.