Sequential assignment of the 1H and 31P resonances of the double stranded deoxynucleotide d (ATGCAT)2 by 2D-NMR correlation spectroscopy

The interactions of substituted pyrido[1,2-e]purines with oligonucleotides depend on the amphiphilic properties of their side chain

Three pyrido[1,2-e]purines of increasing hydrophilicity have been synthesized to evaluate as anticancer agents. These drugs interact quite differently with a synthetic oligodeoxynucleotide d(CGATCG)2. [1] is very hydrophobic due to a phenyl residue in its side chain. It only shows limited interactions with the minihelix without any evidence of intercalation. [2] and [3], on the other hand, have one ([2]) or two ([3]) hydroxyl groups in their acyl chain and present rather amphiphilic properties. The result is a similar intercalation of these derivatives between C and G base pairs as revealed by intermolecular nOe, 1H and 31P chemical shift variations. Models for the intercalation of [2] are proposed using energy minimizations and molecular dynamics (MD) calculations subject to restraints from nOe connectivities. Simulations and experiments indicate weak stability and thus fast exchange of [2] in its intercalation site.

AT selectivity and DNA minor groove binding: modelling, NMR and structural studies of the interactions of propamidine and pentamidine with d(CGCGAATTCGCG)2

A molecular modelling strategy has been developed to identify potential binding sites for bis(amidine) ligands in the minor groove of duplex DNA. Calculations of interaction energy for propamidine and pentamidine with d(CGCGAAT TCGCG)2 show that this duplex contains two symmetrically equivalent binding sites of identical affinity, each displaced by 0.3-0.4 bp from the centre of the AT segment. The ligands occupy groove sites spanning approximately 4 and 4-5 bp, respectively with asymmetric binding to the 5'-AATT sequence. The DNA-bis(amidine) interactions have been examined by high-resolution 1H-NMR. The patterns of induced changes in DNA proton chemical shift and the DNA-ligand NOEs confirm that both agents bind in the AT minor groove in a non-centrosymmetric fashion. Detailed structures were determined for each complex using a NOE-restrained simulated annealing procedure, showing that the B-type DNA conformation is not significantly altered upon complexation with either ligand. The free DNA duplex has previously been shown to be extensively hydrated in the minor groove [Kubinec, M.G. and Wemmer, D.E. (1992) J. Am, Chem. Soc. 114, 8739-8740 Liepinsh, E. Otting, G. and Wüthrich, K. (1992) Nucleic Acids Res. 20. 6549-6553]. We detect hydration water close to the A(H2) protons in the presence of propamidine, which may stabilise certain waters against exchange. This conclusion supports recent crystallographic analyses, suggesting that such ligands may use water molecules to bridge between amidinium protons and host DNA bases Details of the ligand interactions with AT-tract DNA duplexes can now be compared for the subsequences 5'-AAT, 5'-AATT and 5'-AAATTT.

Conformational flexibility in DNA duplexes containing single G.G mismatches

Purine-purine mismatches can base-pair in a variety of configurations depending on solution conditions. The G.G mismatch, which also occurs in the G-quartet structure, has been shown by both x-ray crystallography and NMR to adopt G(anti).G(syn) mispairs, with very different hydrogen bonding patterns [Skelly, J., Edwards, K., Jenkins, T. C. & Neidle, S. (1993) Proc. Natl Acad. Sci. USA 90, 804-808; Cognet, J. A. H., Gabarro-Arpa, J., Le Bret, M., van der Marel, G. A. van Boom, J. H. & Fazakerley, G. V. (1991) Nucleic Acids. Res. 19, 6771-6779] while we have recently suggested the presence of weakly hydrogen-bonded G(anti).G(anti) pairs in solution [Borden, K. L. B., Jenkins, T. C., Skelly, J. V., Brown, T. & Lane, A. N. (1992) Biochemistry 31, 5411-5422]. Spectral overlap and additional exchange processes have made detailed structural analysis difficult in these mismatched oligomers. We have used NMR to characterise the conformations of four duplexes containing single G.G mismatches, including a nonamer d(CATCGGATG), two undecamers d(GCATTGAATGC) and d(CATGTGACGTG) that can each form a self-complementary duplex with a single G.G mispair in the centre, and the non-self-complementary d(GTAACGACATG).d(CATGTGGTTAC). The three self-complementary duplexes have a single set of NMR resonances, and all four duplexes show evidence of conformational exchange at the mismatch site. The N1H resonances of the mismatched G residues each integrate to two protons, ruling out the enol tautomer. They resonate between 10.5-10.7 ppm, far upfield of the Watson-Crick hydrogen-bonded GN1H and exchange readily with water protons. Intraresidue GH8-H1' NOE intensities are two-threefold larger for the mismatched G residues than in G.C base pairs, indicating the presence of syn conformations. NOE time courses for the self-complementary duplexes were consistent with an equimolar mixture of G(syn).G(anti) and G (anti).G(syn) states. By symmetry, these states must be interconverting at a rate that is fast on the chemical shift timescale. In the non-self-complementary undecamer, the NOE data indicated that the distinguishable mismatched G residues also spend a significant, but different, fraction of the time in both the syn and anti conformations. The rate constant for the syn/anti transition in the non-self-complementary undecamer was determined as approximately 14,000 s-1 at 303 K from rotating frame T1 measurements, and the apparent frequency difference was > 250 Hz. Calculations based on NOEs and coupling constants showed that the duplexes are overall in the B form.(ABSTRACT TRUNCATED AT 250 WORDS)

Interaction of minor-groove-binding diamidine ligands with an asymmetric DNA duplex. NMR and molecular modelling studies

The aromatic diamidines berenil and propamidine bind reversibly to A+T-rich sites in the minor groove of B-form DNA duplexes. Based on extensive solution and crystallographic information we have designed a non-self-complementary double-stranded DNA sequence, d(GCAATGAGCG).d(CGCTCATTGC), that should contain a near-ideal binding site for berenil and a poorer site for the larger propamidine molecule, viz. d(AAT).d(ATT). 1H-NMR studies show that both ligands bind with 1:1 stoichiometry to the embedded 5'-AAT site, and induce numerous shifts of NMR resonances of DNA protons located in the minor groove. In addition, interactions with each strand can be distinguished by NOE spectroscopy due to the inherent asymmetry of the DNA. Detailed modelling based on experimental data show that no significant distortion of the B-DNA duplex is induced by either ligand. Sufficient NOE data were obtained to determine the position of the bound ligands in each complex. These conclusions are in agreement with predictions from molecular modelling calculations that provide a microscopic energy profile for interaction with the minor groove tract. Such calculations reveal an unexpected heterogeneous 5'-ATGA binding site that includes a spanned guanosine. This secondary binding site accounts for the extensive chemical shift perturbation induced by these ligands. The structures of the free DNA and the reversible complexes formed with each ligand molecule have been refined using an NOE-restrained isothermal annealing procedure. These structures confirm that the introduced ligands effect minimal perturbation of the helix, with binding to the 5'-AAT base sequence largely determined by specific non-bonded interactions.

Two mutant binding sites of the activating transcription factor region within the E2A promoter of adenovirus exist in a novel conformation

The solution conformations of two mutants of the ATF (activating transcription factor) binding site within the E2A promoter of adenovirus have been determined in vitro by NMR and CD methods. Both sequences have conformations which incorporate north-like sugar puckers in helices which are stacked in a B-like manner as seen with the parent binding site (Borden, K.L.B. (1993) Biochemistry 32, 6506-6514). The PATFm sequence has similar binding affinities and specificities to ATF while the PM2 oligonucleotide is recognized by a different subset of proteins within the ATF family. Both sequences contain unusual amounts of sugar puckers in north-like conformations but the specific distribution of north-like and south-like structures differs between them. These data indicate that the existence of this novel conformation is not characteristic of only the parent sequence. Further, there is a sequence dependent component as illustrated by the variation of the distribution of the north-like sugar puckers within the two mutant oligonucleotides. Differences in sugar pucker conformation can cause bending of the helix and will alter the phosphate backbone surface of the oligonucleotide. Both factors are important to the protein nucleic acid recognition process and thus to cellular control of transcription.

Properties of multiple G.A mismatches in stable oligonucleotide duplexes

The solution structure of the deoxydecanucleotide [sequence: see text] has been determined by NMR methods. This duplex, which contains six G.A mismatches and four Watson-Crick base pairs, is thermodynamically more stable than a decamer where T.A base pairs are substituted for the G.A mismatches, and is less stable than the duplex that contains G.C base pairs. Circular-dichroism spectroscopy indicates an overall B-like conformation for the decamer, but stronger than usual base stacking. 1H-NMR spectroscopy revealed that the N1H groups of the mismatched guanine residues are not hydrogen bonded, and 31P-NMR showed the presence of BII phosphate conformations for the GpA steps. Detailed analysis of the NMR data showed that all nucleotides have anti glycosidic torsion angles and S type sugar puckers. The G.A mismatches pair in the amino form as originally proposed by Li et al. [Li, Y., Zon, G. & Wilson, W. D. (1991) Proc. Natl Acad. Sci. USA 88, 26-30], which results in extensive base-base stacking between the tandem G.A base pairs and their nearest neighbours. The terminal G.A base pairs are less stable than the central base pairs and show evidence of an equilibrium between two conformations, one involving BII phosphate.

Thermodynamic stability and solution conformation of tandem G.A mismatches in RNA and RNA.DNA hybrid duplexes

G.A mismatches form a variety of hydrogen-bonded structures in DNA, most of which destabilise the duplex. Tandem G.A mismatches in the context YGAR (Y = pyrimidine, R = purine), however, form base pairs using the amino group of the guanine residue [Li, Y., Zon. G. & Wilson, W.D. (1991) Proc. Natl Acad. Sci. USA 88, 26-30], which permits extensive base-base stacking, leading to a slight stabilisation of the helix [Ebel, S., Lane, A. N. & Brown, T. (1992) Biochemistry 31, 12083-12086]. We have measured the thermodynamic stability of several RNA and RNA.DNA hybrid duplexes containing tandem G.A mismatches. The RNA duplexes are intrinsically much more stable than the corresponding DNA duplexes and the mutations are destabilising in all cases. NOE and coupling-constant data show that all of the sugars are in the C3'-endo range of conformations, and glycosidic torsion angles are in the range -160 degrees to -180 degrees in [sequence: see text]. Both sequential NOE intensities and circular-dichroism measurements indicate that the global conformation of the mismatched RNA is A-like. The N1H group of the mismatched guanine residue is not involved in hydrogen bonding with the adenine residue, indicating the presence of the amino-pairing scheme. Determination of the structure using 'loose' NMR-derived constraints shows that the potential energies of the imino-paired and amino-paired forms are similar, but substantially higher than energy-minimised RNA. Using tighter constraints derived from more extensive analysis of one-dimensional and two-dimensional NOE data showed that the amino-paired structure agrees with the constraint data better than the imino-paired structure, and also accounts for unusual chemical shifts and the lack of hydrogen bonding of the guanine N1H group. Resulting molecular models show that the amino-paired mismatches are not as extensively stacked on the neighbouring part of the duplex as in the B-DNA analogues, largely accounting for the lower thermodynamic stability in the RNA duplexes.

The activating transcription factor region within the E2A promoter exists in a novel conformation

The ATF (activating transcription factor) binding site within the E2A promoter region of adenovirus is shown to exist in a novel conformation in vitro via nuclear magnetic resonance methods. This novel conformation may be important to the protein DNA recognition process. This conformation has characteristics of both A- and B-form DNA. From circular dichroism and through-space-based NMR experiments, it is clear that the overall helical structure is B-like. However, the 1H-1H coupling constant information indicates that most of the sugar puckers of the individual nucleotides are in the C3'-endo/C4'-exo range which is more characteristic of A-form DNA. The sugar conformation can also be described as a mixture of two states, C3'-endo and C2'-endo, where many of the sugars exist mainly in the C3'-endo state. These data show that the conformation of the sugar puckers does not determine the nature of the overall base stacking on the DNA. Helical parameters were calculated from NOE build-up curves for half of the dinucleotide pairs. Severe spectral overlap on the nuclear Overhauser based spectra prevented determination of the helical parameters for all of the dinucleotide base-pairs. Energy minimization and molecular dynamic simulation methods using the sugar pucker and glycosidic torsion angles determined from the NMR data as constraints were carried out in order to demonstrate that such a conformation was energetically favorable.

NMR and molecular modeling studies of the interaction of berenil and pentamidine with d(CGCAAATTTGCG)2

The interaction of two anti-trypanosomal agents, berenil and pentamidine, with the A+T-rich dodecamer d(CGCAAATTTGCG)2 has been examined by high-resolution 1H-NMR, optical spectroscopy, and molecular modeling. Proton assignments for the free DNA and each DNA-ligand complex were obtained using nuclear Overhauser enhancement spectroscopy and total correlation spectroscopy. Complexation induces large changes in chemical shift for protons in the DNA minor groove for the A5-T9 segment, and intermolecular NOEs reveal contacts between the DNA bases and each ligand. The asymmetric binding site for berenil indicated by the NMR data suggests that at least two overlapping sites are involved. Rapid exchange between symmetrically-equivalent binding sites, via dissociative rearrangement, is consistent with retention of twofold degeneracy for both the ligand and the DNA host. Calculations of binding energy confirm that this DNA duplex contains overlapping sites of similar binding affinity. In contrast, the larger pentamidine molecule occupies a site that spans four or five bp, with asymmetric binding to the minor-groove 5'-ATTT sequence. The B-type conformation of the DNA is not altered substantially by either ligand.

Conformational properties of the G.G mismatch in d(CGCGAATTGGCG)2 determined by NMR

The conformational properties of the DNA duplex d(CGCGAATTGGCG)2, which contains two noncomplementary G.G base pairs, have been examined in aqueous solution by 1H and 31P NMR as a function of temperature. The G.G mismatch is highly destabilizing, with a Tm value 35 K below that observed for the native EcoRI dodecamer. The dodecamer appears symmetric in the NMR spectra and exists largely as an average B-type DNA conformation. However, the 1H and 31P NMR spectra give evidence of considerable conformational heterogeneity at the mismatched nucleotides and their nearest neighbors, which increases with increasing temperature. There is no evidence for a significant population of the syn purine conformation. The imino protons of the mispaired bases G4 and G9 are degenerate, resonate at high field, and exchange readily with solvent. These results indicate that the mispaired bases are only weakly hydrogen-bonded and are only partially stacked into the helix. On raising the temperature, the duplex shows increasing exchange between two or more conformations originating from the mismatch sites. However, these additional conformations maintain their Watson-Crick hydrogen bonding. The increase in chemical exchange is consistent with a quasimelting process for which the G.G sites provide local nuclei. Extensive modeling studies by dynamic annealing have confirmed that the G(anti).G(anti) conformation is favored and that the mispairs are poorly stacked within the helix. The results explain both the poor thermal stability and low hypochromicity of this duplex.

Sequence-specific NMR assignments of the trp repressor from Escherichia coli using three-dimensional 15N/1H heteronuclear techniques

Sequence-specific 15N and 1H assignments for the trp holorepressor from Escherichia coli are reported. The trp repressor consists of two identical 107-residue subunits which are highly helical in the crystal state [Schevitz, R., Otwinowski, Z., Joachimiak, A., Lawson, C. L. & Sigler, P. B. (1985) Nature 317, 782-786]. The high helical content and the relatively large size of the protein (Mr = 25,000) make it difficult to assign even the main-chain resonances by conventional homonuclear two-dimensional NMR methods. However, we have now assigned the main-chain resonances of 94% of the residues by using three-dimensional 15N/1H heteronuclear experiments on a sample of protein uniformly labelled with 15N. The additional resolution obtained by spreading out the signals into three dimensions proved indispensable in making these assignments. In particular, we have been able to resolve signals from residues in the N-terminal region of the A helix for the first time in solution. The observed NOE results confirm that the repressor is highly helical in solution, and contains no extended chain conformations.

Solution conformation and dynamics of the octadeoxy-nucleotide d(CACTAGTG)2: a multinuclear n.m.r. relaxation study

The conformation and internal dynamics of the octadeoxynucleotide d(CACTAGTG)2 have been examined by 1H- and 13C-n.m.r. relaxation. All the non-exchangeable protons, the seven phosphate resonances, and most of the 13C resonances of the proton-bearing carbons have been assigned by conventional two-dimensional n.m.r. methods. The average conformations of each nucleotide have been determined using time-dependent one-dimensional n.O.e.'s and 3JHH values derived from both NOESY and 2-quantum-filtered COSY experiments. All glycosidic torsion angles are anti, and in the range -95 to 125 degrees, in which the pyrimidines have a significantly larger angle than the purines. All sugars were found mainly (greater than 80%) in the conformation range C-2'endo to C-3'exo. The DNA fragment is within the B-family of conformations. The cytosine H-6-H-5 vectors move with an apparent correlation time of 3 ns at 25 degrees. Cross-relaxation rate constants for the H-1'-H-2b vectors and some H-2a-H-2b and H-2a-H-3' vectors were measured, from which order parameters were determined. The order parameters are all in the range 0.7-0.9, which is consistent with only moderate internal mobility on the sub-ns time scale. The (1H)-13C n.O.e. and the spin-lattice relaxation rate constant show that the terminal residues are relatively more mobile than the internal residues, and that the C-2'-H and C-3'-H vectors move with order parameters of 0.6-0.75.

Determination of the orientations of tryptophan analogues bound to the trp repressor and the relationship to activation

The antirepressor indole 3-propanoate has been shown by X-ray crystallography to bind in a different orientation compared with the natural corepressor for the tryp repressor, L-tryptophan (Lawson, C.L. & Sigler, P. B. (1988) Nature 333, 869-871). This suggests a simple difference between what constitutes a corepressor versus an antirepressor. We have used visible absorption and 1H-NMR spectroscopy to characterise the nature of several ligand-repressor complexes and DNA-binding assays to assess the relative operator binding affinities. 5-Fluorotryptophan binds with similar affinity and in the same orientation as L-tryptophan, and is an equally effective corepressor. In contrast, the tight-binding antirepressor indole 3-acrylate binds in the same orientation as indole 3-propanoate. Indole, also an antirepressor, also binds in the indole-3-propanoate orientation. 5-Methyltryptamine, a corepressor, shows spectroscopic characteristics of both tryptophan and indoleacrylate, though NOEs indicate that the tryptophan orientation is preferred. These results indicate that the ammonium group in the side chain is essential both for activation and binding in the L-tryptophan orientation. Antirepressors, lacking the ammonium group, bind in the more favourable indole-3-propanoate orientation. Differences in the NMR signatures of the different repressor-ligand complexes indicate that the details of the conformations depend on the nature of the ligands and their orientation within the binding site. Despite any conformational rearrangement of the protein on binding, dissociation of ligands is facile: 5-fluorotryptophan dissociates rapidly at 313 K. These findings complement and extend the X-ray and thermodynamic analyses of ligand binding.

N.m.r. determination of the solution conformation and dynamics of the A.G mismatch in the d(CGCAAATTGGCG)2 dodecamer

A.G base-paired mismatches that occur during replication are among the most difficult to detect by repair enzymes. Such purine.purine mispairs can exist in two conformations, one of which is stabilized by protons [Gao & Patel (1988) J. Am. Chem. Soc. 110, 5178-5182]. We have undertaken a 1H-n.m.r. and 31P-n.m.r. study of the mismatched dodecamer d(CGCAAATTGGCG)2 as a function of both temperature and pH to determine the conformational features of the A.G mismatch. At pH greater than 7 the mispaired bases are each in the anti conformation and are stacked in the B-like helix. As the pH is decreased, a second conformation becomes populated (apparent pKa approx. 5.9) with concomitant changes in the chemical shifts of protons of the mispaired bases and their nearest neighbours. Data from two-dimensional nuclear-Overhauser-enhancement spectroscopy show unequivocally that, at low pH, the dominant conformation is one in which the mismatched G residues are in the syn conformation and are hydrogen-bonded to the A residues that remain in the anti conformation. Residues not adjacent to the A.G sites are almost unaffected by the transition or the mispairing, suggesting considerable local flexibility of the unconstrained duplexes. Despite the bulging of the mispaired bases, the conformation of the A(anti).G(anti) duplex is very similar to the native dodecamer, whereas the AH+(anti).G(syn) duplex shows a greater variation in the backbone conformation at the mismatched site. According to the chemical shifts, the duplex retains twofold symmetry in solution. The equilibrium between the syn and anti conformations of G9/G21 is strongly dependent on pH, but only weakly dependent on temperature (delta H approx. 16 kJ.mol-1). The first-order rate constant for the transition is approx. 9 s-1 at 283 K and approx. 60 s-1 at 298 K, with an activation enthalpy of approx. 100 kJ.mol-1. The stabilization of the A(anti).G(syn) conformation by protons is consistent with models invoking N1 protonation of adenine. Using the derived glycosidic torsion angles we have used restrained molecular dynamics to build models of the neutral and protonated d(CGCAAATTGGCG)2 oligomers. The results confirm that the A(anti).G(anti) and AH+(anti).G(syn) conformations are favoured at high pH and low pH respectively, in accord with n.m.r. and single-crystal X-ray data.

1H- and 31P-NMR studies of ditercalinium binding to a d(GCGC)2 and d(CCTATAGG)2 minihelices: a sequence specificity study

The structures of the complexes formed in aqueous solution between ditercalinium, a bis-intercalating drug, and both the self-complementary tetranucleotide d(GCGC)2 and octanucleotide d(CCTATAGG)2, have been investigated by 400-MHz 1H-nmr and 162-MHz 31P-nmr. All the nonexchangeable protons, as well as the exchangeable imino protons and the phosphorus signals, have been assigned. Both oligonucleotides have been shown to adopt a right-handed B-DNA type structure. The addition of ditercalinium to the oligonucleotides lead to the formation of complexes in slow exchange at the nmr time scale with the free helices. At all drug-to-helix ratios studied, the ditercalinium was found in the bound form, whereas free and complexed oligonucleotides were in slow exchange, allowing resonance assignments through two-dimensional chemical exchange experiments. for d(GCGC)2 the strong upfield shifts induced on all aromatic protons of both the bases and the drug by complexation with ditercalinium suggest an interaction by intercalation of the two rings. However, the loss of twofold symmetry upon binding, as well as the chemical shift variation of the drug proton signals of one of the chromophores with temperature and concentration, favor a model in which the drug-nucleotide complexes have one ring of the drug intercalated and the other stacked on top of the external base pair. The intermolecular contacts between drug protons and nucleotide protons give a defined geometry for complexation that is consistent with the proposed model. In contrast, with d(CCTATAGG)2 several drug-nucleotide complexes were formed and a large increase in line broadening was observed at high drug-to-DNA ratios, precluding a detailed analysis of these complexes. However, the large upfield shift in the imino proton resonances together with the shielding of the ditercalinium ring protons favor a model with bis-intercalation of ditercalinium. This model is supported by the downfield shift of at least 4 out of 14 phosphorus signals. The results are compared with those obtained on ditercalinium binding to the homologous sequences d(CGCG)2 and d(TTCGCGAA)2, and discussed in terms of sequence specificity.

Interaction of berenil with the EcoRI dodecamer d(CGCGAATTCGCG)2 in solution studied by NMR

The conformation of the EcoRI dodecamer d(CGCGAATTCGCG)2 has been examined in solution by 1H and 31P NMR. Spin-spin coupling constants and nuclear Overhauser (NOE) enhancement spectroscopy show that all deoxyriboses lie in the south domain, with a small admixture of the north conformation (0-20%). The time dependence of the nuclear Overhauser enhancements also reveals a relatively uniform conformation at the glycosidic bonds (average angle, chi = -114 degrees). The average helical twist is 36.5 degrees (9.8 base pairs per turn). Tilt angles are small (in the range 0 to -10 degrees), and roll angles are poorly determined. Unlike single-crystal X-ray studies of the same sequence, there is no evidence for asymmetry in the structure. Both the NOE intensities and 31P relaxation data imply conformational anomalies at the C3-G4/C9-G10 and the A5-A6/T7-T8 steps. Berenil binds in 1:1 stoichiometry to the dodecamer with high affinity (Kd = 1 microM at 298 K) and causes substantial changes in chemical shifts of the sugar protons of nucleotides Ado 5-Cyt 9 and of the H2 resonances of the two Ado residues. No significant asymmetry appears to be induced in the DNA conformation on binding, and there is no evidence for intercalation, although the binding site is not centrosymmetric. NOEs are observed between the aromatic protons of berenil and the H1' of both Thy 7 and Thy 8, as well as to Ado 5 and Ado 6 H2. These results firmly establish that berenil binds via the minor groove and closely approaches the nucleotides Ado 6, Thy 7, and Thy 8. On the basis of quantitative NOE spectroscopy and measurements of spin-spin coupling constants, changes in the conformations of the nucleotides are found to be small. Using the observed NOEs between the ligand and the DNA together with the derived glycosidic torsion angles, we have built models that satisfy all of the available solution data. The berenil molecule binds at the 5'-AAT (identical to 5'-ATT on the complementary strand) site such that (i) favorable hydrogen bonds are formed between the charged amidinium groups and the N3 atoms of Ado 6 and Ado 18 and (ii) the ligand is closely isohelical with the floor of the minor groove.

The solution conformations of a mutant trp operator determined by n.m.r. spectroscopy

The principal conformational features of the mutant trp operator [d(CGTACTGATTAATCAGTACG)2] have been determined by n.m.r. at different temperatures. The sugar puckers were determined from J-resolved spectroscopy and high-resolution homonuclear Hartmann-Hahn spectroscopic (HOHAHA) experiments. Extensive one-dimensional nuclear-Overhauser-enhancement (NOE) data sets were acquired at 25 degrees C using irradiation times of 50, 100, 200, 300 and 500 ms to generate sufficient NOE information to determine the individual nucleotide conformations, and place limits on the local helical parameters, using multi-spin least-squares fitting and searching in conformation space with the program NUCFIT [Lane (1990) Biochim. Biophys. Acta 1049, 189-204]. The conformations of the nucleotides are well determined, and show significant sequence-dependent variation. Pyrimidine residues on average have a wider range of sugar conformations and smaller glycosidic torsion angles than purine residues. The helical parameters are in general less well determined, though clear evidence was obtained for sequence-dependent variation of the helical twist. The overall mean fractional deviation of the calculated from the observed NOEs was 0.108. The conformations of the base-pairs TAAT are temperature-dependent [Lane (1989) Biochem. J. 259, 715-724]. NOESY spectra were recorded at 10, 25 and 40 degrees C, using mixing times inversely proportional to the overall tumbling time to allow changes in the conformation to be described. A more detailed analysis was made using one-dimensional NOEs collected for nucleotides involved in the conformational transitions. There are significant temperature-dependent changes in the conformations of the central base-pairs from T9 to T13 with the largest changes in the glycosidic torsion angle occurring for A11 and A12 (up to 30 degrees). The orientation of the base-pairs T9-A12:T10-A11 also changes, with an increase in the base-pair roll and an unwinding of the helix as the temperature is increased. The conformational changes are qualitatively similar to those observed in a related sequence (Lefèvre, Lane & Jardetzky (1988) Biochemistry 27, 1086-1094]. The conformation is also similar to the wild-type sequence and to that observed in the crystal state in the complex with the trp holorepressor. The similarity suggests that the mutation produces a poorer operator by virtue of removal of essential functional groups within the major groove.

NMR studies of the interaction of chromomycin A3 with small DNA duplexes. Binding to GC-containing sequences

The interaction of chromomycin A3 with the oligodeoxyribonucleotides 1, d(ATGCAT), 2, d(ATCGAT), 3, d(TATGCATA), and 4, d(ATAGCTAT), has been investigated by 1H and 31P NMR. In the presence of Mg2+, chromomycin binds strongly to the three GC-containing oligomers 1, 3, and 4 but not to the CG-containing oligomer 2. The proton chemical shift changes for 1 and 3 are similar, and these DNA duplexes appear to bind with a stoichiometry of 2 drugs:1 Mg2+:1 duplex. The same stoichiometry of 2 drugs:1 duplex is confirmed with 4; however, proton chemical shift changes differ. An overall C2 symmetry is exhibited by the drug complex with 1, 3, and 4. At a molar ratio of 2.0 (drugs:duplex), no free DNA proton NMR signals remain. Two-dimensional nuclear Overhauser exchange spectroscopy (NOESY) of the saturated chromomycin complex with 1 and 3 positions both chromomycinone hydroxyls and the E carbohydrates in the minor groove and provides evidence suggesting that the B carbohydrates lie on the major-groove side. This is supported by several dipolar coupling cross-peaks between the drug and the DNA duplex. Drug-induced conformational changes in duplex 1 are evaluated over a range of NOESY mixing times and found to possess some characteristics of both B-DNA and A-DNA, where the minor groove is wider and shallower. A widening of the minor groove is essential for the DNA duplex to accommodate two drug molecules. This current minor-groove model is a substantial revision of our earlier major-groove model [Keniry, M.A., Brown, S.C., Berman, E., & Shafer, R.H. (1987) Biochemistry 26, 1058-1067] and is in agreement with the model recently proposed by Gao and Patel [Gao, X., & Patel, D. J. (1989a) Biochemistry 28, 751-762].

Solution conformations of d(CGTACG)2 determined from NMR experiments

The conformations of all the nucleotides in the hexamer d(CGTACG)2 have been determined using time-dependent one- and two-dimensional nuclear Overhauser enhancements (NOEs) and the program NUCFIT (see previous article). The glycosidic torsion angles are well determined, the fraction of the C2' endo state for the sugar puckers is less well determined, and the pseudorotation phase angle is poorly determined by the NOEs. The average glycosidic torsion angle is -107 +/- 9 degrees, and the deoxyriboses of the purine residues have a higher fraction of the C2' endo state than those of the pyrimidine residues. There is good agreement between the one- and two-dimensional NOE data. Of the helical parameters, the local rise and twist are moderately well determined, but the roll and tilt of the bases are not well described. The overall structure belongs to the B family of conformations, as previously described by Gronenborn et al. (Biochem. J. (1984) 221, 723-736), but there are significant differences which can be ascribed to the improved treatment of the spin-diffusion and motional averaging possible with the program NUCFIT. The results obtained using NUCFIT are compared with those from restrained energy minimisation calculations using distance restraints obtained from NUCFIT.

Reassessment of structural characteristics of the d(CGCG)2:actinomycin D complex from complete 1H and 31P NMR

Complexes formed between Actinomycin D (ActD) and the tetranucleotides d(AGCT)2 and d(CGCG)2 were studied in detail by one and two-dimensional 1H and 31P NMR. The 31P two dimensional chemical exchange experiment, at room temperature on saturated complexes (1:1), showed unambiguously that the asymmetrical phenoxazone ring binds to the unique GC site under the two possible orientations in the d(AGCT)2 tetranucleotide but adopts a single orientation in the d(CGCG)2 tetranucleotide. For the d(CGCG)2:Act D saturated complex, complete assignments of all protons and phosphorus signals of the two-nucleotide strands, as well as of the two cyclic pentapeptide chains has allowed us to study in details the conformational features of the complex from NOE and coupling constants analysis. The tetranucleotide remains in a right-handed duplex, but the sugar puckers are modified for residues at the intercalation site. A uniform C2' endo pucker is observed for residues on the strand facing the quinoid side of the phenoxazone ring while a C2' endo-C3-endo equilibrium about 60% of C2' endo is proposed for the two residues on the strand facing the benzenoid side of the phenoxazone ring. In contrast to previous studies on ActD-DNA interactions, we have been able to measure the 3J phosphorus-proton coupling constants at the intercalation site but also adjacent to it, showing that 31P chemical shifts are not simply related to the backbone conformation. Molecular mechanics calculations, using empirical distances deduced from NOE effects as restrained distances during minimizations, led to a model differing mainly from those previously published by orientation of the N methyl groups of both N-Methyl-Valines.

Solution structure of the parallel-stranded duplex oligonucleotide alpha-d(TCTAAAC)-beta-d(AGATTTG) via complete relaxation matrix analysis of the NOE effects and molecular mechanics calculations

The solution structure of the duplex formed by the association of the unnatural oligonucleotide alpha-d(TCTAAAC) with its natural and parallel complementary sequence beta-d(AGATTTG) was investigated by nuclear magnetic resonance spectroscopy and constrained molecular mechanics calculations. The structure was refined on the basis of interproton distances determined by NOE measurements for a series of mixing times. The NOE values were converted to distances by using the complete 134 x 134 relaxation matrix including all proton dipole-dipole interactions and spin diffusion. The computation of the relaxation matrix requires the Cartesian coordinates of the oligonucleotide, which are not known, a priori. To avoid this ambiguity, we used an iterative procedure in which the new distance constraints are obtained by using the complete relaxation matrix calculated from the previous structure. After three iterations, the process converged. The unnatural duplex alpha-d(TCTAAAC)-beta-d(AGATTTG) adopts in solution a right-helical structure with Watson-Crick base pairing, an anti conformation on the glycosyl linkage on the beta-strand, a syn conformation on the alpha-strand, and a 3'-exo conformation of the deoxyriboses for both sugar anomers. The three-dimensional structure obtained allowed us to describe the local heterogeneity of the duplex.

The influence of tryptophan on mobility of residues in the trp repressor of Escherichia coli

The relative mobility of residues in the trp repressor of Escherichia coli has been examined in the absence and presence of the corepressor L-tryptophan by one- and two-dimensional 1H NMR. A comparison of relative intensities of cross peaks in NOESY and COSY spectra allowed a rigid Tyr and a mobile Tyr residue, three mobile Ser residues and three mobile Lys residues to be detected. The two Tyr residues were assigned by selective nitration with tetranitromethane. The singly nitrated molecule (on Tyr7) binds the trp operator with an affinity close to that of the unmodified repressor. Measurements of the intraring cross-relaxation rate constant as a function of temperature for Tyr7 shows the presence of considerable internal motion on the subnanosecond time scale in the flexible N-terminal arm. The order parameter, S2, characterising the motion is 0.35, which increases to about 0.5 in the presence of Trp. Trp decreases both the amplitude of the motion and the rate of the motion. At least three of the six Ser residues of the trp repressor have greater mobility than expected for a rigid body, and two of the Ser residues are sensitive to the presence of Trp. The more mobile Ser residues are probably those on the N-terminal arm and the C-terminal sequence. These results complement the single-crystal X-ray diffraction studies for which the electron density of the first ten and last three amino acid residues is weak. The solution data are consistent with proposals that the flexible N-terminal arm of the trp repressor makes important contacts with the DNA.

N.m.r. assignments and temperature-dependent conformational transitions of a mutant trp operator-promoter in solution

A total of 145 protons in the mutant trp operator-promoter sequence CGTACTGATTAATCAGTACG were assigned by one-dimensional and two-dimensional n.m.r. methods. Except at the sites of mutation (underlined), the chemical shifts and other n.m.r. parameters are very similar to those observed in the symmetrized wild-type sequence [Lefèvre, Lane & Jardetzky (1987) Biochemistry 26, 5076-5090]. Spin-spin-relaxation rate constants of the resolved base protons and intra- and inter-nucleotide nuclear-Overhauser-enhancement intensities argue for a sequence-dependent structure similar to that of the wild-type, except at and close to the sites of the mutation. The overall tumbling time as a function of temperature was determined from cross-relaxation rate constants for the H-6-H-5 vectors of the four cytosine residues. The values are consistent with the oligonucleotide maintaining a double-helical conformation over the entire temperature range 5-45 degrees C, and that internal motions of the bases are of small amplitude on the subnanosecond time scale. The temperature-dependence of chemical shifts, spin-spin-relaxation rate constants and cross-relaxation rate constants show the occurrence of two conformational transitions localized to the TTAA sequence in the centre of the molecule. The thermodynamics of the transition at the lower temperature (tm = 16 degrees C) were analysed according to a two-state process. The mid-point temperature is about 6 degrees C higher than in the wild-type sequence. The conformational transition does not lead to rupture of the Watson-Crick hydrogen bonds, but probably involves changes in the propellor twists of T.A-9 and T.A-10. The second transition occurs at about 40 degrees C, but cannot be fully characterized. This conformational variability seems to be a property of the sequence TTAA, and may have functional significance in bacterial promoters.

Novel 1H-31P heteronuclear coherence transfer spectroscopy based on spin locking MLEV-17 and its application to signal assignment for a short DNA fragment

We have applied the MLEV-17 proton spin locking pulse to 1H-31P 2D heteronuclear correlation spectroscopy. Using this technique, long-range correlations between a proton and phosphorus that is up to 5 bonds distant are observed. Thus, the 3'-, 4'- and 5'-proton resonances can be traced 'sequentially' along the chain on the 2D NMR chart. Application of this technique to the complete assignment of the sugar proton and phosphorus resonances for d(ApGpA) is reported. This technique can also be used for convenient assignment of the phosphorus resonances of oligonucleotides as correlations between the 1'-proton and its 3'-side phosphorus are observed.

Why do all lariat RNA introns have adenosine as the branch-point nucleotide? Conformational study of naturally-occurring branched trinucleotides and its eleven analogues by 1H-, 31P-NMR and CD spectroscopy

1H-NMR conformational studies of six branched triribonucleotides where the branch-point nucleotide was either U, C or G (4-9) have been carried out by assigning 1H resonances through 2D NMR and then observing the temperature-dependent (i) chemical shifts of the aromatic and the anomeric protons, and (ii) shifts of the equilibrium of N and S pseudorotamer populations of each sugar moiety. The data have been compared with those of 2'----5' dimers (1-3) and other branched trimers (10-16). It emerged that all the branched trimers (4-16) adopt a conformational state closer to the corresponding 2'----5' dimers than the corresponding 3'----5' dimers. A temperature-dependent 31P chemical shift study confirmed that the conformational constraint is mainly associated with the 2'----5' phosphate linkage. Although, it appeared with the CD data that when C or especially when U is at the branch-point the overall constraint is weak. This suggests that even if these trimers adopt a 2'----5' dimer geometry, there is a lack of stabilization by strong stackings within the molecule. This is in sharp contrast with the results found for A (10-16) and to a smaller extent for G (8, 9) at the branch-point.

Practical aspects of 2D NMR for assigning the non-exchangeable protons in DNA-RNA fragments

Using the branched trimer A2'-5'A3'-5'A as an example, different 2D NMR experiments such as homonuclear correlations via single-quantum coherence, double-quantum filtration via double-quantum coherence, isotropic mixing, relayed connectivities, cross relaxation and 31P/1H correlations are presented together with the corresponding spectra. The discussion pointed out the advantages and the difficulties of 14 pulse sequences, used for assigning non-exchangeable protons in DNA or RNA fragments. The emphasis has been put on the methodological aspects of 2D NMR as a simple technique, allowing an easier assimilation of the new coming 2D NMR pulse sequences.

Alterations in the d(CpGpT) structure in solution as a result of [PtCl(diethylenetriamine)]+ binding

The trinucleotide d(CpGpT) reacts with [PtCl(dien)]Cl (dien = diethylenetriamine) to yield as a single adduct Pt(dien)[d(CpGpT)-N7(2)]. The structure of this adduct in solution has been analysed with the aid of NMR spectroscopy and compared with that of the unmodified trinucleotide. A change in the population of the S conformer of the guanosine deoxyribose ring and a syn preference of the guanine residue are the most important changes occurring upon platination. As a result the dC-dG stack disappears, whereas the dG-dT stack is hardly affected. The CD spectra of both platinated and free d(CpGpT) confirm the different nature of the two molecules.

Study of structure, base-pair opening kinetics and proton exchange mechanism of the d-(AATTGCAATT) self-complementary oligodeoxynucleotide in solution

Using proton magnetic resonance, we have investigated the structure and the base-pair opening kinetics of the d-(AATTGCAATT) self-complementary duplex. All the non-exchangeable (except H5',5") and most exchangeable proton resonances have been assigned. The structure belongs to the B family. Imino proton exchange, measured by line broadening, longitudinal relaxation and magnetization transfer from water, is catalyzed by proton acceptors. The base-pair lifetimes, obtained by extrapolation of the exchange times to infinite concentration of ammonia are 2 and 3 milliseconds for internal A.Ts and 18 ms for G.C at 15 degrees C. In the absence of added catalysts, the imino proton of the first A.T base pair exchanges faster than that of the unpaired thymidine of the duplex formed by the sequence d-(AATTGCAATTT). This gives strong evidence for intrinsic exchange catalysis. The exchange of adenine amino protons from the closed state has been observed. Hence amino proton exchange is ill-suited for the investigation of base-pair opening kinetics.

Resonance assignments of non-exchangeable protons in B type DNA oligomers, an overview

The chemical shifts of 1H resonances of non exchangeable protons (except H5', H5" and adenine H2) of over six hundred nucleotides have been collected. The influence which the base of the nucleotide itself as well as the bases on its 5' and 3' side exert on the chemical shifts of the various resonances has been investigated. Most of the resonances appear to be predominantly influenced by only one base. For H2', H2", H3', H4' and H6/H8 this is the base of the central nucleotide, for H5(C) and CH3(T) it is the one on the 5' side and for H1' it is the one on the 3' side. Chemical shift distribution profiles are presented which allow an estimation of the probability of finding a particular resonance at a particular position in the spectrum.

2D-NMR studies of the unnatural duplex alpha-d(TCTAAAC)-beta-d(AGATTTG)

The unnatural oligonucleotide alpha-d(TCTAAAC) was synthesized and was found more resistant towards endonucleases than its beta-analog. 2D-NMR experiments allowed the assignment of all non-exchangeable aromatic and sugar protons except for the overlapping 5' -5" resonances, as well as the exchangeable imino protons of the parallel hybrid duplex alpha-d (TCTAAAC)-beta-d(AGATTTG). NMR studies show that the strength of the association between the alpha-strand and the beta parallel strand is equivalent to that between their anti-parallel complementary beta-analogs beta-d(CAAATCT) and beta-d(AGATTTG). NOE data provide evidence that both duplexes form stable right-helical duplexes with an anti-conformation on the glycosyl linkages and a Watson-Crick pairing. NOESY and COSY spectra allowed us to determine that alpha and beta deoxyriboses adopt a 3' -exo conformation.

DNA and RNA: NMR studies of conformations and dynamics in solution

The early NMR research on nucleic acids was of a qualitative nature and was restricted to partial characterization of short oligonucleotides in aqueous solution. Major advances in magnet design, spectrometer electronics, pulse techniques, data analysis and computational capabilities coupled with the availability of pure and abundant supply of long oligonucleotides have extended these studies towards the determination of the 3-D structure of nucleic acids in solution.

An NMR study of the polymorphous behavior of the mismatched DNA octamer d(m5C-G-m5C-G-T-G-m5C-G) in solution. The B, Z, and hairpin forms

The polymorphism exhibited by the mismatched octamer d(m5C-G-m5C-G-T-G-m5C-G), as a function of the temperature, DNA concentration and ionic strength, was investigated by means of NMR spectroscopy. It is shown that this partly self-complementary DNA fragment, under conditions of low DNA concentration (0.4 mM) and low ionic strength, exclusively prefers to adopt a monomeric hairpin form, which consists of a stem of three Watson-Crick-type base pairs and a loop of only two residues. This in striking contrast with earlier intimations in literature, which postulated that in oligonucleotides loop formations containing only two residues are sterically impossible. Moreover, the hairpin form displays an unusual stability in comparison with previously reported hairpins. A Tm of 332 K and a delta H degree of -130 kJ.mol-1 were calculated for the hairpin to random coil transition. At high DNA concentration (8 mM) and/or upon the addition of sodium chloride the hairpin form occurs in slow exchange with a B-DNA dimer structure (approximately 20% at 270 K, no added salt), which comprises two central GxT-mismatched base pairs with the bases as major tautomers. At higher ionic strength (greater than 100 mM NaCl), or upon the addition of methanol, a third species appears, which is in slow exchange with both the B dimer and the hairpin form. This third species could be identified with a Z DNA form, comprising two GxT mismatches with the bases as major tautomers, with the guanine bases syn and the cytosine and thymine bases anti.

Alkyl phosphotriester modified oligodeoxyribonucleotides. VI. NMR and UV spectroscopic studies of ethyl phosphotriester (Et) modified Rp-Rp and Sp-Sp duplexes, (d[GGAA(Et)TTCC])2

1H NMR chemical shift assignments for the title compounds were made for all but a few H5' and H5" signals using two-dimensional nuclear Overhauser effect (2D-NOE) data, which was also used for the first time to assign absolute configuration at phosphorus. The chemical shifts were, in general, similar to those reported [Broido, M.S., et al. (1985) Eur. J. Biochem. 150, 117-128] for the B-like conformation of the unmodified, parent duplex, [d(GGAATTCC)]2. Differences in chemical shifts for corresponding protons were mostly localized to the AA(Et)TT region, and showed some stereochemical dependence. Unambiguous assignment of the phosphotriester 31P signals was achieved in a novel way using selective insensitive nucleus enhancement by polarization transfer (selective INEPT) NMR. The Rp-Rp duplex melted ca. 11 degrees C lower than either the Sp-Sp or parent duplexes, as evidenced by Tm and variable temperature 1H/31P NMR measurements. The 2D-NOE data for the Rp-Rp duplex suggested possible steric interactions between the ethyl group and the H3' of the flanking A residue. At low ionic strength, the Sp-Sp and parent duplexes had similar stability but at high ionic strength the Sp-Sp duplex was less stable.

31P NMR studies of the binding of the oligonucleotide (Ap)3A to an oligodeoxythymidylate covalently linked to an acridine derivative

31P NMR was used to study the specific interaction of an oligodeoxynucleotide containing four thymines and covalently attached to an acridine derivative through its 3'-phosphate [(Tp)4(CH2)5Acr] with a complementary oligoribonucleotide (Ap)3A. 31P-1H and 1H-1H chemical shift correlation spectroscopies were jointly used to provide the assignment of the phosphorus resonances. A downfield shift of two phosphorus resonances of (Tp)4(CH2)5Acr and of two phosphorus resonances of (Ap)3A was observed upon complex formation. The assignment of the phosphorus resonances which are downfield shifted allowed us to propose a model involving an equilibrium between several 1:1 complexes where the acridine ring is intercalated between different A.T base pairs.

The 31P-NMR spectrum of the dodecamer d(GACGATATCGTC)

The resonances in the 31P-NMR spectrum of the dodecamer d(GACGATATCGTC) have been assigned by regiospecific labelling with oxygen-17. All 11 resonances are clearly resolved at 26 degrees C. Most noticeably, individual resonances of the dinucleoside phosphates d(CpG), d(TpC), d(GpA) and d(ApT) which occur more than once can clearly be distinguished. This indicates that the position of the phosphate group in the oligomer influences its 31P-NMR shift. This observation is in agreement with what has been found for the 31P-NMR spectra of d(CGCGAATTCGCG) [Ott, J. and Eckstein, F. (1985) Biochemistry 24] and d(GGAATTCC) [Connolly, B.A. and Eckstein, F. (1984) Biochemistry 23, 5523-5527]. In general, the chemical shift appears the more at higher field the more central the dinucleoside phosphate is located in the oligomer. Exceptions are the resonances of dinucleoside phosphates of the type 5'-PyPu-3' which appear at lower field than expected from this rule. A reasonable correlation between 31P-NMR chemical shifts and the sum function of the base plane roll angles derived from Calladine's rule [Calladine, C.R. (1982) J. Mol. Biol. 161, 343-352] exists.