Proton exchange and internal motions in two chromomycin dimer-DNA oligomer complexes

Individual basepair stability of DNA and RNA studied by NMR-detected solvent exchange

In this study, we have optimized NMR methodology to determine the thermodynamic parameters of basepair opening in DNA and RNA duplexes by characterizing the temperature dependence of imino proton exchange rates of individual basepairs. Contributions of the nuclear Overhauser effect to exchange rates measured with inversion recovery experiments are quantified, and the influence of intrinsic and external catalysis exchange mechanisms on the imino proton exchange rates is analyzed. Basepairs in DNA and RNA have an approximately equal stability, and the enthalpy and entropy values of their basepair dissociation are correlated linearly. Furthermore, the compensation temperature, T(c), which is derived from the slope of the correlation, coincides with the melting temperature, and duplex unfolding occurs at that temperature where all basepairs are equally thermodynamically stable. The impact of protium-deuterium exchange of the imino hydrogen on the free energy of RNA basepair opening is investigated, and it is found that two A·U basepairs show distinct fractionation factors.

HEXIM1 targets a repeated GAUC motif in the riboregulator of transcription 7SK and promotes base pair rearrangements

7SK snRNA, an abundant RNA discovered in human nucleus, regulates transcription by RNA polymerase II (RNAPII). It sequesters and inhibits the transcription elongation factor P-TEFb which, by phosphorylation of RNAPII, switches transcription from initiation to processive elongation and relieves pauses of transcription. This regulation process depends on the association between 7SK and a HEXIM protein, neither isolated partner being able to inhibit P-TEFb alone. In this work, we used a combined NMR and biochemical approach to determine 7SK and HEXIM1 elements that define their binding properties. Our results demonstrate that a repeated GAUC motif located in the upper part of a hairpin on the 5'-end of 7SK is essential for specific HEXIM1 recognition. Binding of a peptide comprising the HEXIM Arginine Rich Motif (ARM) induces an opening of the GAUC motif and stabilization of an internal loop. A conserved proline-serine sequence in the middle of the ARM is shown to be essential for the binding specificity and the conformational change of the RNA. This work provides evidences for a recognition mechanism involving a first event of induced fit, suggesting that 7SK plasticity is involved in the transcription regulation.

Direct observation of the temperature-induced melting process of the Salmonella fourU RNA thermometer at base-pair resolution

In prokaryotes, RNA thermometers regulate a number of heat shock and virulence genes. These temperature sensitive RNA elements are usually located in the 5'-untranslated regions of the regulated genes. They repress translation initiation by base pairing to the Shine-Dalgarno sequence at low temperatures. We investigated the thermodynamic stability of the temperature labile hairpin 2 of the Salmonella fourU RNA thermometer over a broad temperature range and determined free energy, enthalpy and entropy values for the base-pair opening of individual nucleobases by measuring the temperature dependence of the imino proton exchange rates via NMR spectroscopy. Exchange rates were analyzed for the wild-type (wt) RNA and the A8C mutant. The wt RNA was found to be stabilized by the extraordinarily stable G14-C25 base pair. The mismatch base pair in the wt RNA thermometer (A8-G31) is responsible for the smaller cooperativity of the unfolding transition in the wt RNA. Enthalpy and entropy values for the base-pair opening events exhibit linear correlation for both RNAs. The slopes of these correlations coincide with the melting points of the RNAs determined by CD spectroscopy. RNA unfolding occurs at a temperature where all nucleobases have equal thermodynamic stabilities. Our results are in agreement with a consecutive zipper-type unfolding mechanism in which the stacking interaction is responsible for the observed cooperativity. Furthermore, remote effects of the A8C mutation affecting the stability of nucleobase G14 could be identified. According to our analysis we deduce that this effect is most probably transduced via the hydration shell of the RNA.

Sequence-dependent base pair opening in DNA double helix

Preservation of genetic information in DNA relies on shielding the nucleobases from damage within the double helix. Thermal fluctuations lead to infrequent events of the Watson-Crick basepair opening, or DNA "breathing", thus making normally buried groups available for modification and interaction with proteins. Fluctuational basepair opening implies the disruption of hydrogen bonds between the complementary bases and flipping of the base out of the helical stack. Prediction of sequence-dependent basepair opening probabilities in DNA is based on separation of the two major contributions to the stability of the double helix: lateral pairing between the complementary bases and stacking of the pairs along the helical axis. The partition function calculates the basepair opening probability at every position based on the loss of two stacking interactions and one base-pairing. Our model also includes a term accounting for the unfavorable positioning of the exposed base, which proceeds through a formation of a highly constrained small loop, or a ring. Quantitatively, the ring factor is found as an adjustable parameter from the comparison of the theoretical basepair opening probabilities and the experimental data on short DNA duplexes measured by NMR spectroscopy. We find that these thermodynamic parameters suggest nonobvious sequence dependent basepair opening probabilities.

A nuclear magnetic resonance investigation of the energetics of basepair opening pathways in DNA

The opening of basepairs plays a key role in DNA replication and transcription, and in the action of DNA repair and modification enzymes. In this article, we have used proton exchange to define the energetics of the pathways for basepair opening in two DNA 17-mer duplexes. The rates of exchange of imino protons with solvent protons were measured by NMR spectroscopy for each DNA duplex, as a function of the concentration of exchange catalyst and of temperature. The measurements provided the rates and the equilibrium constants of the opening reactions for individual basepairs at different temperatures. These temperature dependences were used to calculate the enthalpies and the free energies of the barrier to opening and of the open state for each basepair. The results reveal the existence of three distinct patterns of enthalpy changes in the opening reactions. The patterns differ from each other in the location of the kinetic opening barrier relative to the open state. Neighboring bases, which are one or more positions removed from the opening basepair, influence the enthalpic pattern of the opening pathway. The free energies of the opening barriers are found to be linearly related to the free energies of the open state. This correlation is analyzed in terms of rate-equilibrium free energy relationships previously observed in other systems, and suggests that the transition state in the opening reaction is closer to the native closed state of the basepair than to its open state.

1H NMR determination of base-pair lifetimes in oligonucleotides containing single base mismatches

Proton nuclear magnetic resonance (NMR) spectroscopy is employed to characterize the kinetics of base-pair opening in a series of 9mer duplexes containing different single base mismatches. The imino protons from the different mismatched, as well as fully matched, duplexes are assigned from the imino-imino region in the WATERGATE NOESY spectra. The exchange kinetics of the imino protons are measured from selective longitudinal relaxation times. In the limit of infinite exchange catalyst concentration, the exchange times of the mismatch imino protons extrapolate to much shorter lifetimes than are commonly observed for an isolated GC base pair. Different mismatches exhibit different orders of base-pair lifetimes, e.g. a TT mismatch has a shorter base-pair lifetime than a GG mismatch. The effect of the mismatch was observed up to a distance of two neighboring base pairs. This indicates that disruption in the duplex caused by the mismatch is quite localized. The overall order of base-pair lifetimes in the selected sequence context of the base pair is GC > GG > AA > CC > AT > TT. Interestingly, the fully matched AT base pair has a shorter base-pair lifetime relative to many of the mismatches. Thus, in any given base pair, the exchange lifetime can exhibit a strong dependence on sequence context. These findings may be relevant to the way mismatch recognition is accomplished by proteins and small molecules.

HIV-1 nucleocapsid protein activates transient melting of least stable parts of the secondary structure of TAR and its complementary sequence

The nucleocapsid protein NCp7 of HIV-1 possesses a nucleic acid chaperone activity that is critical in minus and plus strand transfer during reverse transcription. The minus strand transfer notably relies on the ability of NCp7 to destabilize the stable stem with five contiguous, double-stranded segments of both the TAR sequence at the 3' end of the viral genome and the complementary sequence, cTAR, in minus strong-stop DNA. In order to examine the nature and the extent of NCp7 destabilizing activity, we investigated, by absorbance and fluorescence spectroscopy, the interaction of TAR and cTAR with a (12-55)NCp7 peptide containing the zinc-finger motifs but lacking the ability to aggregate the oligonucleotides. The absorbance changes in the UV band of cTAR show that seven to eight base-pairs, on average, are melted per oligonucleotide at a ratio of one peptide to 7.5 nucleotides. In contrast, the melting of TAR does not exceed an average of one base-pair per oligonucleotide. This may be linked to the greater stability of TAR, since a strong correlation between NCp7 destabilizing effect and oligonucleotide stability was observed. The effect of (12-55)NCp7 on the stem terminus was investigated by using a cTAR molecule doubly labeled at the 3' and 5' ends by a donor/acceptor couple. In the absence of the peptide, about 80 % of the oligonucleotides are in a dark non-fluorescent state, having a close proximity of the two dyes. The remaining 20 % are distributed between three fluorescent species, having either the terminal segment, the two terminal segments or all segments of the stem melted. This is in line with a fraying mechanism wherein the stem terminus fluctuates rapidly between open and closed states. Addition of (12-55)NCp7 shifts the equilibrium toward the open species, suggesting that NC enhances fraying of the stem terminus. Taken together, our data suggest that NCp7 activates the transient opening of base-pairs in the least stable parts of the stem. Also, this activity of NCp7 was found to be dependent on the zinc-finger motifs, since no melting was observed with a fingerless NCp7 peptide.

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.

High base pair opening rates in tracts of GC base pairs

Sequence-dependent structural features of the DNA double helix have a strong influence on the base pair opening dynamics. Here we report a detailed study of the kinetics of base pair breathing in tracts of GC base pairs in DNA duplexes derived from 1H NMR measurements of the imino proton exchange rates upon titration with the exchange catalyst ammonia. In the limit of infinite exchange catalyst concentration, the exchange times of the guanine imino protons of the GC tracts extrapolate to much shorter base pair lifetimes than commonly observed for isolated GC base pairs. The base pair lifetimes in the GC tracts are below 5 ms for almost all of the base pairs. The unusually rapid base pair opening dynamics of GC tracts are in striking contrast to the behavior of AT tracts, where very long base pair lifetimes are observed. The implication of these findings for the structural principles governing spontaneous helix opening as well as the DNA-binding specificity of the cytosine-5-methyltransferases, where flipping of the cytosine base has been observed, are discussed.

Determination of the residence time of water molecules hydrating B'- DNA and B-DNA, by one-dimensional zero-enhancement nuclear Overhauser effect spectroscopy

The residence time of water in the minor groove of the d(CGCGAATTCGCG) duplex has been determined by a recent measurement combining nuclear Overhauser enhancements (NOE, ROE) and 17O relaxation dispersion. The time is in the range of nanoseconds, so that it may be measured by a rather simple method proposed here, namely the choice of conditions such that the NOE between the observed DNA proton and a nearby water proton is zero. This condition is realized when the residence time of the water molecule is 0.178 times the nuclear magnetic resonance period (e.g. 0.297 ns at 600 MHz). It may be achieved by varying the magnetic field and/or the temperature. The zero-NOE measurement may be performed by one-dimensional NMR, and has therefore good sensitivity. We have developed excitation sequences which suppress two spurious contributions to the NOE: from neighboring exchangeable protons and from H3' protons whose chemical shift is close to that of water. The method is applied here to the comparison of residence times of water next to B-DNA and next to B'-DNA, the latter corresponding to better stacked, propeller-twisted base-pairs and a correspondingly narrower minor groove. In the minor groove of [d(CGCGAATTCGCG)]2, a B'-DNA duplex, the residence time of the water molecule next to H2 of adenine(6) (underlined), is 0.6 ns at 10 degreesC, in good agreement with the value obtained previously. The residence time is slightly but distinctly shorter for the water next to A5, suggesting non-cooperative departure of these two molecules which are presumed to be part of the hydration spine. Near A5 and A4 of [d(AAAAATTTTT)]2, another B'-DNA duplex, the residence times are approximately twice as long, but the activation enthalpies are about the same, ca. 38 kJ/mol. The residence time in the minor groove of the regular B-DNA sequence d(CGCGATCGCG) was 0.3 ns at 10 degreesC, shorter than in the case of the B'-DNA sequences by factors of 2 and 4, respectively. The temperature dependence is less, with an activation enthalpy of 27 kJ/mol. The major groove residence times are comparable for the three sequences, and a few times shorter than those of minor groove water. A value of 0.36 ns, or even more in case of rotation of water, is obtained around -8 degreesC. The most striking aspect of these results is the relatively small difference in the residence times of reputedly fast and slow-exchanging water molecules bound to DNA in biological conditions. This suggests that the spine of hydration is perhaps not a major stabilizer of the B'-DNA structure as compared with B-DNA.

Formation of an amino-acid-binding pocket through adaptive zippering-up of a large DNA hairpin loop

BACKGROUND

In vitro selection has identified DNA aptamers that target cofactors, amino acids, peptides and proteins. Structure determination of such ligand-DNA aptamer complexes should elucidate the details of adaptive DNA structural transitions, binding-pocket architectures and ligand recognition. We have determined the solution structure of the complex of a DNA aptamer containing a guanine-rich 18-residue hairpin loop that binds L-argininamide with approximately 100 microM affinity.

RESULTS

The DNA aptamer generates its L-argininamide-binding pocket by adaptive zippering up the 18-residue loop through formation of Watson-Crick pairs, mismatch pairs and base triples, while maximizing stacking interactions. Three of the four base triples involve minor-groove recognition through sheared G.A mismatch formation. The unique fold is also achieved through positioning of an adenine residue deep within the minor groove and through nestling of a smaller loop within the larger loop on complex formation. The accessibility to the unique L-argininamide-binding pocket is restricted by a base pair that bridges across one side of the major-groove-binding site. The guanidinium group of the bound L-argininamide aligns through intermolecular hydrogen-bond formation with the base edges of nonadjacent guanine and cytosine residues while being sandwiched between the planes of nonadjacent guanine residues.

CONCLUSIONS

The available structures of L-arginine/L-argininamide bound to their DNA and RNA targets define the common principles and patterns associated with molecular recognition, as well as the diversity of intermolecular hydrogen-bonding alignments associated with the distinct binding pockets.

Proton exchange kinetics in [d(ACGTATACGT)]2-echinomycin and [d(ACGTTAACGT)]2-echinomycin complexes

Based on imino proton exchange catalysis, base-pair lifetimes and apparent dissociation constants are reported on the complexes formed by bisintercalation of echinomycin at the CpG steps of the d(ACGTATACGT)2 and d(ACGTTAACGT)2 duplexes. The lifetimes of the four central A x T base pairs between two echinomycin binding sites are much shorter than in the free duplexes. The destabilization of base pairs adjacent to the binding sites is propagated one additional base pair away from the binding site.

Unique base-pair breathing dynamics in PNA-DNA hybrids

Kinetic and thermodynamic parameters, derived from 1H-NMR measurements of the imino proton exchange rates upon titration with the exchange catalyst ammonia, are reported for two mixed-sequence peptide nucleic acid (PNA)-DNA hybrids and their counterpart DNA duplex. The exchange times of the imino protons in the PNA strands extrapolate to very short base-pair lifetimes in the limit of infinite exchange catalyst concentration. This is not due to generally less stable base-pairs in PNA-DNA hybrids, since the lifetimes, apparent dissociation constants and thermodynamic stability (DeltaG degrees ) of the innermost DNA guanine imino protons are similar in the hybrid duplexes and in the DNA duplex. In addition, the apparent dissociation constants determined for PNA bases of the hybrids are of the same order as those of the corresponding bases in the DNA duplex. An exchange process from the closed state was found to be inconsistent with the experimental data. From these results, we conclude that opening and closing rates of the PNA guanine and thymine bases are at least two orders of magnitude higher than those of the corresponding bases in the DNA duplex. Unusual kinetics in the hybrids is also evident from the destabilization of the complementary DNA strand thymine bases, which exhibit base-pair dissociation constants increased by approximately two orders of magnitude compared to what is observed in the DNA duplex, while the DNA strand guanine bases are largely unaffected. The general pattern of the base-pair dynamics in the hybrids obtained when using trimethylamine as an exchange catalyst is the same as when using ammonia. However, the long base-pair lifetimes i. e. those of the DNA duplex and the guanine bases of the DNA strands in the hybrids, are approximately three to five times longer than when using ammonia. Thus, all opening events sensed by ammonia are not accessible to trimethylamine. These observations are discussed in regard to the mechanism of base-pair opening and the nature of the open state.

Imino proton exchange and base-pair kinetics in the AMP-RNA aptamer complex

We report on the dynamics of base-pair opening in the ATP-binding asymmetric internal loop and flanking base-pairs of the AMP-RNA aptamer complex by monitoring the exchange characteristics of the extremely well resolved imino protons in the NMR spectrum of the complex. The kinetics of imino proton exchange as a function of basic pH or added ammonia catalyst are used to measure the apparent base-pair dissociation constants and lifetimes of Watson-Crick and mismatched base-pairs, as well as the solvent accessibility of the unpaired imino protons in the complex. The exchange characteristics of the imino protons identify the existence of four additional hydrogen bonds stabilizing the conformation of the asymmetric ATP-binding internal loop that were not detected by NOEs and coupling constants alone, but are readily accommodated in the previously reported solution structure of the AMP-RNA aptamer complex published from our laboratory. The hydrogen exchange kinetics of the non-Watson-Crick pairs in the asymmetric internal loop of the AMP-RNA aptamer complex have been characterized and yield apparent dissociation constants (alphaKd) that range from 10(-2) to 10(-7). Surprisingly, three of these alphaKd values are amongst the lowest measured for all base-pairs in the AMP-RNA aptamer complex. Comparative studies of hydrogen exchange of the imino protons in the free RNA aptamer and the AMP-RNA aptamer complex establish that complexation stabilizes not only the bases within the ATP-binding asymmetric internal loop, but also the flanking stem base-pairs (two pairs on either side) of the binding site. We also outline some preliminary results related to the exchange properties of a sugar 2'-hydroxyl proton of a guanosine residue involved in a novel hydrogen bond that has been shown to contribute to the immobilization of the bound AMP by the RNA aptamer, and whose resonance is narrow and downfield shifted in the spectrum.

The effect of cross-links on the conformational dynamics of duplex DNA

The base pair lifetimes and apparent dissociation constants of a 21 base DNA hairpin and an analog possessing a disulfide cross-link bridging the 3'- and 5'-terminal bases were determined by measuring imino proton exchange rates as a function of exchange catalyst concentration and temperature. A comparison of the lifetimes and apparent dissociation constants for corresponding base pairs of the two hairpins indicates that the cross-link neither increases the number of base pairs involved in fraying nor alters the lifetime, dissociation constant, or the opened structure from which exchange occurs for the base pairs that are not frayed. The cross-link does, however, stabilize the frayed penultimate base pair of the stem duplex. Significantly, it appears that the disulfide cross-link is more effective at preventing fraying of the penultimate base pair than is the 5 base hairpin loop. Because this disulfide cross-link can be incorporated site specifically, and does not adversely affect static or dynamic properties of DNA, it should prove very useful in studies of nucleic acid structure and function.

Acid-induced exchange of the imino proton in G.C pairs

Acid-induced catalysis of imino proton exchange in G.C pairs of DNA duplexes is surprisingly fast, being nearly as fast as for the isolated nucleoside, despite base-pair dissociation constants in the range of 10(-5) at neutral or basic pH. It is also observed in terminal G.C pairs of duplexes and in base pairs of drug-DNA complexes. We have measured imino proton exchange in deoxyguanosine and in the duplex (ATATAGATCTATAT) as a function of pH. We show that acid-induced exchange can be assigned to proton transfer from N7-protonated guanosine to cytidine in the open state of the pair. This is faster than transfer from neutral guanosine (the process of intrinsic catalysis previously characterized at neutral ph) due to the lower imino proton pK of the protonated form, 7.2 instead of 9.4. Other interpretations are excluded by a study of exchange catalysis by formiate and cytidine as exchange catalysts. The cross-over pH between the regimes of pH-independent and acid-induced exchange rates is more basic in the case of base pairs than in the mononucleoside, suggestive of an increase by one to two decades in the dissociation constant of the base pair upon N7 protonation of G. Acid-induced catalysis is much weaker in A.T base pairs, as expected in view of the low pK for protonation of thymidine.

NMR studies and restrained-molecular-dynamics calculations of a long A+T-rich stretch in DNA. Effects of phosphate charge and solvent approximations

The nonamer duplex d(GCAAAAACG).d(CGTTTTTGC) was studied by 1H-NMR at 500 MHz. With the exception of the H5' and H5" sugar protons, all protons were assigned by two-dimensional NMR experiments [NOE spectroscopy (NOESY), double-quantum-filtered J-correlated spectroscopy (DQF-COSY) and total correlation spectroscopy (TOCSY)]. The exchange kinetics of the imino protons of the Watson-Crick base pairing were studied at 15 degrees C by measuring inversion-recovery rates under conditions of extensive ammonia base catalysis. Extrapolation to infinite base concentration gave anomalous long lifetimes for the A-tract in accordance with previous results [Leroy, J.-L., Charettier, E., Kochoyan, M. & Guéron, M. (1988) Biochemistry 27, 8894-8898]. On average, 11 NOESY distance constraints/nucleotide were evaluated using the complete relaxation matrix approach. Deoxyribose coupling constants were obtained from simulations of the DQF-COSY cross-peaks, assuming a rapid two-state equilibrium between a C2'-endo and C3'-endo conformer. The sugars were found to be predominantly in the C2'-endo conformation. The NMR-derived distance and torsion constraints were implemented into three different restrained-molecular-dynamics (rMD) protocols, two in vacuo, with different charges on the phosphate group and the third with the solvent explicitly included. All protocols displayed good convergence from different starting structures. The structures derived from the three protocols satisfied experimental restraints equally well and had similar final energies. Although the overall pattern of sequence dependence of helical parameters shows some resemblance in all structures, we find that the absolute amplitudes of the parameters are largely dependent on the rMD protocols, particularly the twist parameters. The minor groove distance P(n + 2)-P(m + 2) varies from 0.7 nm to 1.2 nm in the three protocols. Still the NOESY-derived anomalously short distances AH2(n)-H1'(m + 1) and AH2(n)-H1'(n + 1), n and m denote complementary residues, which are assumed to be indicative of a compressed minor groove, are kept in all calculated structures.

Base-pair opening and spermine binding--B-DNA features displayed in the crystal structure of a gal operon fragment: implications for protein-DNA recognition

A sequence that is represented frequently in functionally important sites involving protein-DNA interactions is GTG/CAC, suggesting that the trimer may play a role in regulatory processes. The 2.5 A resolution structure of d(CGGTGG)/d(CCACCG), a part of the interior operator (OI, nucleotides +44 to +49) of the gal operon, co-crystallized with spermine, is described herein. The crystal packing arrangement in this structure is unprecedented in a crystal of B-DNA, revealing a close packing of columns of stacked DNA resembling a 5-stranded twisted wire cable. The final structure contains one hexamer duplex, 17 water molecules and 1.5 spermine molecules per crystallographic asymmetric unit. The hexamer exhibits base-pair opening and shearing at T.A resulting in a novel non-Watson-Crick hydrogen-bonding scheme between adenine and thymine in the GTG region. The ability of this sequence to adopt unusual conformations in its GTG region may be a critical factor conferring sequence selectivity on the binding of Gal repressor. In addition, this is the first conclusive example of a crystal structure of spermine with native B-DNA, providing insight into the mechanics of polyamine-DNA binding, as well as possible explanations for the biological action of spermine.

Solution structure of the mithramycin dimer-DNA complex

We have characterized the NMR parameters for the complexes formed by the Mg(2+)-coordinated mithramycin dimer with self-complementary d(T-G-G-C-C-A) and d(T-C-G-C-G-A) duplexes. The solution structure of the latter complex has been determined using a combined NMR-molecular dynamics study including relaxation matrix refinement. The Mg(2+)-coordinated mithramycin dimer-d(T-C-G-C-G-A) complex exhibits a 2-fold center of symmetry with the divalent cation coordinated aglycons positioned opposite the central (G3-C4).(G3-C4) segment such that the aglycon C8 hydroxyl oxygens form symmetrical sequence-specific hydrogen bonds to guanine amino protons in the complex. The C-D-E trisaccharide segments of each monomer in the mithramycin dimer adopt extended conformations, are positioned inside the minor groove, and are directed toward either end of the duplex. The C-D saccharide component of one monomer and the aglycon of the other monomer in the mithramycin dimer share a widened minor groove with the hydrophobic edges of the C and D sugars interacting with individual strands of the duplex. The E-sugar ring is positioned in the floor of the minor groove, and its hydroxyl-bearing face interacts with both strands of the duplex through hydrogen-bonding and hydrophobic intermolecular interactions. The A-B disaccharide and the hydrophilic side chain form intermolecular contacts with the sugar-phosphate backbone in the complex. The antiparallel alignment of divalent cation coordinated monomers in the mithramycin dimer results in the two outwardly directed C-D-E trisaccharide segments generating a right-handed continuous hexasaccharide domain that spans six base pairs in the minor groove of the duplex. The solution structure of the mithramycin dimer-DNA complex reported in this study and the solution structure of the chromomycin dimer-DNA complex reported previously [Gao, X., Mirau, P., & Patel, D. J. (1992) J. Mol. Biol. 223, 259-279] show global similarities, as well as local differences that are of interest. All four nucleotides in the tetranucleotide segment of the duplex centered about the sequence-specific (G-C).(G-C) step adopt A-DNA sugar puckers and glycosidic torsion angles in the chromomycin dimer-DNA complex, while only the central cytidine adopts an A-DNA sugar pucker and glycosidic torsion angle in the mithramycin dimer-DNA complex.(ABSTRACT TRUNCATED AT 400 WORDS)

Acid multimers of oligodeoxycytidine strands: stoichiometry, base-pair characterization, and proton exchange properties

The structure recently proposed for the acid form of the oligonucleotide 5'-d(TC5) is a four-strand "tetrad" in which two parallel-stranded, base-paired duplexes are intimately associated, with their hemiprotonated C-C+ base pairs face-to-face and fully intercalated, in a so-called "i-motif" (Gehring et al., 1993). We use the amino and imino proton spectra to establish the structure and symmetry of the base pairs, properties which are a primary element in the resolution of the acid form describe above. The amino proton spectrum gives the best lower limit (8 x 10(4) s-1) on the rate of the imino proton jumping process which is responsible for the base-pair symmetry. The stoichiometry of the acid form of other deoxycytidine sequences is studied by gel filtration chromatography and in one case by an NMR equilibrium titration. In all cases, i.e., d(C12), d(T2C8T2), d(C4TC4), d(TC5), d(C5), d(C4), d(TC4), d(TC3T), and d(TC3), the acid form elutes as a tetramer. A single-strand component is also present in some cases. But no dimer is observed, except for some samples prepared by quenching from high temperatures. The characteristic H1'-H1' interresidue NOESY cross-peaks of the d(TC5) structure (Gehring et al., 1993) are also found in all the tetramers where they have been searched for, i.e., those of d(T2C8T2), d(C4TC4), d(TC3T), and d(TC3) (not shown), suggesting that these tetramers also are built on the i-motif and that such structures may be formed generally by strands containing a stretch of as little as three deoxycytidines. From the NMR titration of d(TC3), we derive a free energy of -7.6 kJ/mol per cytidine base pair for the formation of the tetramer from single strands. The free energy released by packing a base pair into the i-motif is comparable to that released in forming the base pair itself. Imino proton exchange is limited by base-pair opening, thanks to efficient intrinsic exchange catalysis: this explains the lack of effect of added catalysts. The base-pair lifetime is hundreds of times longer than in any DNA duplex, presumably due to the base-pair intercalation geometry. The variation of the lifetime along the sequence of the d(TC5) tetramer provides support for the recently proposed structure. The internal amino proton exchanges from the open state of the C-C+ pair, at a rate compatible with a pK of 9 appropriate for C+. But the external proton exchanges from the closed state, as with a pK of 17!(ABSTRACT TRUNCATED AT 400 WORDS)

Effects of sequence and length on imino proton exchange and base pair opening kinetics in DNA oligonucleotide duplexes

The base catalysed imino proton exchange in DNA oligonucleotides of different sequences and lengths was studied by 1H-NMR saturation recovery experiments. The self-complementary sequences studied were GCGCGAATTCGCGC (I), CGCGAATTCGCG (II), GCGAATTCGC (III), and CGCGATCGCG (IV). The evaluation of base pair lifetimes was made after correction for the measured 'absence of added catalyst' effect which was found to be characterized by recovery times of 400-500 ms for the AT base pairs and 250-300 ms for the GC base pairs at 15 degrees C. End effects with rapid exchange is noticeable up to 3 base pairs from either end of the duplexes. The inner hexamer cores GAATTC of sequences I-II show similar base pair lifetime patterns, around 30 ms for the innermost AT, 5-10 ms for the outer AT and 20-50 ms for the GC base pairs at 15 degrees C. The shorter sequences III and particularly IV show much shorter lifetimes in their central AT base pairs (11 ms and 1 ms, respectively).

Quantitative footprinting analysis of the chromomycin A3--DNA interaction

Chromomycin A3 (CHR) binding to the duplex d(CAAGTCTGGCCATCAGTC).d(GACTGATGGCCAGACTTG) has been studied using quantitative footprinting methods. Previous NMR studies indicated CHR binds as a dimer in the minor groove. Analysis of autoradiographic spot intensities derived from DNase I cleavage of the 18-mer in the presence of various amounts of CHR revealed that the drug binds as a dimer to the sequence 5'-TGGCCA-3',3'-ACCGGT-5' in the 18-mer with a binding constant of (2.7 +/- 1.4) x 10(7) M-1. Footprinting and fluorescence data indicate that the dimerization constant for the drug in solution is approximately 10(5) M-1. Since it has been suggested that CHR binding alters DNA to the A configuration, quantitative footprinting studies using dimethyl sulfate, which alkylates at N-7 of guanine in the major groove, were also carried out. Apparently, any drug-induced alteration in DNA structure does not affect cleavage by DMS enough to be observed by these experiments.

Proton exchange in DNA-luzopeptin and DNA-echinomycin bisintercalation complexes: rates and processes of base-pair opening

Imino proton exchange studies are reported on the complexes formed by bisintercalation of luzopeptin around the two central A.T pairs of the d(CCCATGGG) and d(AGCATGCT) duplexes and of echinomycin around the two central C.G pairs of the d(AAACGTTT) and d(CCAAACGTTTGG) duplexes. The depsipeptide backbone of the drugs occupies the minor groove of the complexes at the bisintercalation site. The exchange time of the amide protons of the depsipeptide rings provides a lower estimate of the complex lifetime: 20 min at 15 degrees C for the echinomycin complexes and 4 days at 45 degrees C for the luzopeptin complexes. The exchange time of imino protons is always shorter than the complex lifetime. Hence, base pairs open even within the complexed oligomers. For the two base pairs sandwiched between the aromatic rings of the drug, the base-pair lifetime is strongly increased, and the dissociation constant is correspondingly reduced. Hence, the lifetime of the open state is unchanged. This suggests similar open states in the free duplex and in the complex. In contrast to the sandwiched base pairs, the base pairs flanking the intercalation site are not stabilized in the complex. Thus, the action of the bisintercalating drug may be compared to a vise clamping the inner base pairs. Analysis suggests that base-pair opening may require prior unwinding or bending of the DNA duplex.

Structure refinement of the chromomycin dimer-DNA oligomer complex in solution

We have refined the initial docking model of the Mg(II)-co-ordinated chromomycin-d(A2G2C2T2) complex (2 drug equivalents per duplex) by a complete relaxation matrix analysis simulation of the two-dimensional nuclear Overhauser effect (NOESY) spectrum of the complex in 2H2O solution. This relaxation matrix refined structure of the complex exhibits the following characteristics. (1) We observe an unwound and elongated duplex that exhibits characteristics distinct from the A and B-DNA family of helices at the central (G-G-C-C).(G-G-C-C) chromomycin dimer binding and flanking sites. On the other hand sugar puckers, glycosidic torsion angles, displacement of the base-pairs from the helix axis and the minor groove width for this central tetranucleotide segment all fall within the A-family of helical parameters. (2) The chromomycin monomers are aligned in a head-to-tail orientation in the Mg(II)-co-ordinated dimer in the complex. The chromophores are aligned with a slight tilt relative to each other and make an angle of 75 degrees between their planes. The C-D-E trisaccharide segments from individual monomers adopt an extended conformation that projects in opposite directions in the dimer. The divalent metal cation is co-ordinated to the O(1) carbonyl and O(9) enolate atoms of the chromophores and aligns them such that the O(9)-Mg-O(9) angle is 170 degrees while all other O-Mg-O angles are in the 95(+/- 15)degrees range. (3) The sequence specificity of the chromomycin dimer for the widened and shallower (G3-G4-C5-C6).(G3-G4-C5-C6) minor groove binding site is associated with intermolecular hydrogen bonds formed between the OH group at C(8) of the chromophore and the minor groove NH2 group at position 2 and N(3) groups of G4 and between the O(1) oxygen of the E-sugar and the minor groove NH2 group at position 2 of G3 in the complex. (4) Additional intermolecular interactions are primarily van der Waals contacts between anomeric and adjacent CH2 protons on each sugar in the C-D-E trisaccharide segments of the chromomycin dimer and the minor groove surface of the DNA. These results provide insights into the induced conformational transitions required to generate a complementary match between the drug dimer and its DNA binding site on complex formation.