Proton exchange and base-pair lifetimes in a deoxy-duplex containing a purine-pyrimidine step and in the duplex of inverse sequence

Homology recognition without double-stranded DNA-strand separation in D-loop formation by RecA

RecA protein and RecA/Rad51 orthologues are required for homologous recombination and DNA repair in all living creatures. RecA/Rad51 catalyzes formation of the D-loop, an obligatory recombination intermediate, through an ATP-dependent reaction consisting of two phases: homology recognition between double-stranded (ds)DNA and single-stranded (ss)DNA to form a hybrid-duplex core of 6-8 base pairs and subsequent hybrid-duplex/D-loop processing. How dsDNA recognizes homologous ssDNA is controversial. The aromatic residue at the tip of the β-hairpin loop (L2) was shown to stabilize dsDNA-strand separation. We tested a model in which dsDNA strands were separated by the aromatic residue before homology recognition and found that the aromatic residue was not essential to homology recognition, but was required for D-loop processing. Contrary to the model, we found that the double helix was not unwound even a single turn during search for sequence homology, but rather was unwound only after the homologous sequence was recognized. These results suggest that dsDNA recognizes its homologous ssDNA before strand separation. The search for homologous sequence with homologous ssDNA without dsDNA-strand separation does not generate stress within the dsDNA; this would be an advantage for dsDNA to express homology-dependent functions in vivo and also in vitro.

Base Dynamics in the HhaI Protein Binding Site

Protein-DNA interactions play an important role in numerous biological functions within the living cell. In many of these interactions, the DNA helix is significantly distorted upon protein-DNA complex formation. The HhaI restriction-modification system is one such system, where the methylation target is flipped out of the helix when bound to the methyltransferase. However, the base flipping mechanism is not well understood. The dynamics of the binding site of the HhaI methyltransferase and endonuclease (underlined) within the DNA oligomer [d(G1A2T3A4G5C6G7C8T9A10T11C12)]2 are studied using deuterium solid-state NMR (SSNMR). SSNMR spectra obtained from DNAs deuterated on the base of nucleotides within and flanking the [5'-GCGC-3']2 sequence indicate that all of these positions are structurally flexible. Previously, conformational flexibility within the phosphodiester backbone and furanose ring within the target sequence has been observed and hypothesized to play a role in the distortion mechanism. However, whether that distortion was occurring through an active or passive mechanism remained unclear. These NMR data demonstrate that although the [5'-GCGC-3']2 sequence is dynamic, the target cytosine is not passively flipping out of the double-helix on the millisecond-picosecond time scale. Additionally, although previous studies have shown that both the furanose ring and phosphodiester backbone experience a change in dynamics upon methylation, which may play a role in recognition and cleavage by the endonuclease, our observations here indicate that methylation has no effect on the dynamics of the base itself.

Dynamics Studies of DNA with Non-canonical Structure Using NMR Spectroscopy

The non-canonical structures of nucleic acids are essential for their diverse functions during various biological processes. These non-canonical structures can undergo conformational exchange among multiple structural states. Data on their dynamics can illustrate conformational transitions that play important roles in folding, stability, and biological function. Here, we discuss several examples of the non-canonical structures of DNA focusing on their dynamic characterization by NMR spectroscopy: (1) G-quadruplex structures and their complexes with target proteins; (2) i-motif structures and their complexes with proteins; (3) triplex structures; (4) left-handed Z-DNAs and their complexes with various Z-DNA binding proteins. This review provides insight into how the dynamic features of non-canonical DNA structures contribute to essential biological processes.

How does temperature impact the conformation of single DNA molecules below melting temperature?

The double stranded DNA molecule undergoes drastic structural changes during biological processes such as transcription during which it opens locally under the action of RNA polymerases. Local spontaneous denaturation could contribute to this mechanism by promoting it. Supporting this idea, different biophysical studies have found an unexpected increase in the flexibility of DNA molecules with various sequences as a function of the temperature, which would be consistent with the formation of a growing number of locally denatured sequences. Here, we take advantage of our capacity to detect subtle changes occurring on DNA by using high throughput tethered particle motion to question the existence of bubbles in double stranded DNA under physiological salt conditions through their conformational impact on DNA molecules ranging from several hundreds to thousands of base pairs. Our results strikingly differ from previously published ones, as we do not detect any unexpected change in DNA flexibility below melting temperature. Instead, we measure a bending modulus that remains stable with temperature as expected for intact double stranded DNA.

Parallel triplex structure formed between stretched single-stranded DNA and homologous duplex DNA

The interaction between the single-stranded DNA and the homologous duplex DNA is essential for DNA homologous repair. Here, we report that parallel triplex structure can form spontaneously between a mechanically extended ssDNA and a homologous dsDNA in protein-free condition. The triplex has a contour length close to that of a B-form DNA duplex and remains stable after force is released. The binding energy between the ssDNA and the homologous dsDNA in the triplex is estimated to be comparable to the basepairing energy in a B-form dsDNA. As ssDNA is in a similar extended conformation within recombinase-coated nucleoprotein filaments, we propose that the parallel triplex may form and serve as an intermediate during recombinase-catalyzed homologous joint formation.

Sequence dependency of canonical base pair opening in the DNA double helix

The flipping-out of a DNA base from the double helical structure is a key step of many cellular processes, such as DNA replication, modification and repair. Base pair opening is the first step of base flipping and the exact mechanism is still not well understood. We investigate sequence effects on base pair opening using extensive classical molecular dynamics simulations targeting the opening of 11 different canonical base pairs in two DNA sequences. Two popular biomolecular force fields are applied. To enhance sampling and calculate free energies, we bias the simulation along a simple distance coordinate using a newly developed adaptive sampling algorithm. The simulation is guided back and forth along the coordinate, allowing for multiple opening pathways. We compare the calculated free energies with those from an NMR study and check assumptions of the model used for interpreting the NMR data. Our results further show that the neighboring sequence is an important factor for the opening free energy, but also indicates that other sequence effects may play a role. All base pairs are observed to have a propensity for opening toward the major groove. The preferred opening base is cytosine for GC base pairs, while for AT there is sequence dependent competition between the two bases. For AT opening, we identify two non-canonical base pair interactions contributing to a local minimum in the free energy profile. For both AT and CG we observe long-lived interactions with water and with sodium ions at specific sites on the open base pair.

The DNA double helix, a molecule that stores biological information, has become an iconic image of biomedical research. In order to use or repair the information it carries, the bases that are stacked in the helix need to be chemically exposed. This can happen either by separating the two strands in the helix or by flipping out individual bases. Here, we focus on the latter process. Usually proteins are involved in interactions with bases, but it is still unclear if bases are pulled out actively by proteins or if they act on spontaneously flipped bases. Although experiments can detect base pair opening, it is difficult to detect which base moves in which direction. Here, we present results from molecular dynamics simulations using a recently developed sampling method which improves the statistics in the simulations by enhancing the probability of the base pair opening event. We observe differences in probability, modes and mechanism of opening that depend not only on the types of the bases in the pair, but also strongly on their neighbors. This provides essential information for understanding how DNA functions.

Enzyme-Substrate Binding Kinetics Indicate That Photolyase Recognizes an Extrahelical Cyclobutane Thymidine Dimer

Escherichia coli DNA photolyase is a DNA-repair enzyme that repairs cyclobutane pyrimidine dimers (CPDs) that are formed on DNA upon exposure of cells to ultraviolet light. The light-driven electron-transfer mechanism by which photolyase catalyzes the CPD monomerization after the enzyme-substrate complex has formed has been studied extensively. However, much less is understood about how photolyase recognizes CPDs on DNA. It has been clearly established that photolyase, like many other DNA-repair proteins, requires flipping of the CPD site into an extrahelical position. Photolyase is unique in that it requires the two dimerized pyrimidine bases to flip rather than just a single damaged base. In this paper, we perform direct measurements of photolyase binding to CPD-containing undecamer DNA that has been labeled with a fluorophore. We find that the association constant of ∼2 × 10(6) M(-1) is independent of the location of the CPD on the undecamer DNA. The binding kinetics of photolyase are best described by two rate constants. The slower rate constant is ∼10(4) M(-1) s(-1) and is most likely due to steric interference of the fluorophore during the binding process. The faster rate constant is on the order of 2.5 × 10(5) M(-1) s(-1) and reflects the binding of photolyase to the CPD on the DNA. This result indicates that photolyase finds and binds to a CPD lesion 100-4000 times slower than other DNA-repair proteins. In light of the existing literature, we propose a mechanism in which photolyase recognizes a CPD that is flipped into an extrahelical position via a three-dimensional search.

Mechanistic aspects of hydration of guanine radical cations in DNA

The mechanistic aspects of hydration of guanine radical cations, G(•+) in double- and single-stranded oligonucleotides were investigated by direct time-resolved spectroscopic monitoring methods. The G(•+) radical one-electron oxidation products were generated by SO4(•-) radical anions derived from the photolysis of S2O8(2-) anions by 308 nm laser pulses. In neutral aqueous solutions (pH 7.0), after the complete decay of SO4(•-) radicals (∼5 μs after the actinic laser flash) the transient absorbance of neutral guanine radicals, G(-H)(•) with maximum at 312 nm, is dominant. The kinetics of decay of G(-H)(•) radicals depend strongly on the DNA secondary structure. In double-stranded DNA, the G(-H)(•) decay is biphasic with one component decaying with a lifetime of ∼2.2 ms and the other with a lifetime of ∼0.18 s. By contrast, in single-stranded DNA the G(-H)(•) radicals decay monophasically with a ∼ 0.28 s lifetime. The ms decay component in double-stranded DNA is correlated with the enhancement of 8-oxo-7,8-dihydroguanine (8-oxoG) yields which are ∼7 greater than in single-stranded DNA. In double-stranded DNA, it is proposed that the G(-H)(•) radicals retain radical cation character by sharing the N1-proton with the N3-site of C in the [G(•+):C] base pair. This [G(-H)(•):H(+)C ⇆ G(•+):C] equilibrium allows for the hydration of G(•+) followed by formation of 8-oxoG. By contrast, in single-stranded DNA, deprotonation of G(•+) and the irreversible escape of the proton into the aqueous phase competes more effectively with the hydration mechanism, thus diminishing the yield of 8-oxoG, as observed experimentally.

Inhibitory effect of hydroxypropyl methylcellulose acetate succinate on drug recrystallization from a supersaturated solution assessed using nuclear magnetic resonance measurements

We examined the inhibitory effect of hydroxypropyl methylcellulose acetate succinate (HPMC-AS) on drug recrystallization from a supersaturated solution using carbamazepine (CBZ) and phenytoin (PHT) as model drugs. HPMC-AS HF grade (HF) inhibited the recrystallization of CBZ more strongly than that by HPMC-AS LF grade (LF). 1D-1H NMR measurements showed that the molecular mobility of CBZ was clearly suppressed in the HF solution compared to that in the LF solution. Interaction between CBZ and HF in a supersaturated solution was directly detected using nuclear Overhauser effect spectroscopy (NOESY). The cross-peak intensity obtained using NOESY of HF protons with CBZ aromatic protons was greater than that with the amide proton, which indicated that CBZ had hydrophobic interactions with HF in a supersaturated solution. In contrast, no interaction was observed between CBZ and LF in the LF solution. Saturation transfer difference NMR measurement was used to determine the interaction sites between CBZ and HF. Strong interaction with CBZ was observed with the acetyl substituent of HPMC-AS although the interaction with the succinoyl substituent was quite small. The acetyl groups played an important role in the hydrophobic interaction between HF and CBZ. In addition, HF appeared to be more hydrophobic than LF because of the smaller ratio of the succinoyl substituent. This might be responsible for the strong hydrophobic interaction between HF and CBZ. The intermolecular interactions between CBZ and HPMC-AS shown by using NMR spectroscopy clearly explained the strength of inhibition of HPMC-AS on drug recrystallization.

Synthesis and characterization of ruthenium(II)-oligopyridine-peptide conjugates. Interactions of the diasteromeres Δ- and Λ-[Ru(bpy)2(4-COY-4'-Mebpy)]Cl2 (Y=Gly-Lys(1)-Lys(2)CONH2, Lys(1)-Gly-Lys(2)CONH2, Lys(1)-Lys(2)-GlyCONH2) with the oligonucleotide d(5'-CGCGAATTCGCG-3')2

Diastereomeric complexes of the general formulae Λ- and Δ-[Ru(bpy)2(4-COY-4'-Mebpy)]Cl2 where bpy=2,2'-bipyridine and Y=Gly-Lys(1)-Lys(2)CONH2, Lys(1)-Gly-Lys(2)CONH2, Lys(1)-Lys(2)-GlyCONH2, were synthesized and characterized. The ability of these compounds to bind to the oligonucleotide duplex d(5'-CGCGAATTCGCG-3') was studied with NMR techniques. Complex Λ-2, Λ-[Ru(bpy)2(4-COLys(1)-Gly-Lys(2)CONH2),4'-Mebpy)]Cl2 (Mebpy=methyl-2,2'-bipyridine), interacts non-specifically causing changes for both complex and oligonucleotide (1)H NMR signals. Both Λ-1, Λ-[Ru(bpy)2(4-COGly-Lys(1)-Lys(2)CONH2),4'-Mebpy)]Cl2 and Λ-3, Λ-[Ru(bpy)2(4-COLys(1)-Lys(2)-GlyCONH2),4'-Mebpy)]Cl2, were bound to the oligonucleotide through both lysine aliphatic chains, indicating that the side chains of the sequential lysines create a kind of "clamp" to connect the complex with the oligonucleotide. Complex Δ-1, Δ-[Ru(bpy)2(4-COGly-Lys(1)-Lys(2)CONH2),4'-Mebpy)]Cl2, interacts with the oligonucleotide duplex with both lysine side chains in a manner similar to Λ-1. Δ-2, Δ-[Ru(bpy)2(4-COLys(1)-Gly-Lys(2)CONH2),4'-Mebpy)]Cl2, interacts with the oligonucleotide with the bipyridine ligands. In addition, the formation of a hydrogen bond between the Gly-NH and the carbonyl groups of the oligonucleotide bases was detected. A completely different binding mode was observed for Δ-3 Δ-[Ru(bpy)2(4-COLys(1)-Lys(2)-GlyCONH2),4'-Mebpy)]Cl2, which at a ratio of 1:1 ([Ru]/[nucleotide]) opens the oligonucleotide strands. In addition, participation of all three peptidic NH of Δ-3 in hydrogen bonds was observed.

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.

Temperature dependence of DNA persistence length

We have determined the temperature dependence of DNA persistence length, a, using two different methods. The first approach was based on measuring the j-factors of short DNA fragments at various temperatures. Fitting the measured j-factors by the theoretical equation allowed us to obtain the values of a for temperatures between 5°C and 42°C. The second approach was based on measuring the equilibrium distribution of the linking number between the strands of circular DNA at different temperatures. The major contribution into the distribution variance comes from the fluctuations of DNA writhe in the nicked circular molecules which are specified by the value of a. The computation-based analysis of the measured variances was used to obtain the values of a for temperatures up to 60°C. We found a good agreement between the results obtained by these two methods. Our data show that DNA persistence length strongly depends on temperature and accounting for this dependence is important in quantitative comparison between experimental results obtained at different temperatures.

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.

Heteroduplex joint formation free of net topological change by Mhr1, a mitochondrial recombinase

Homologous pairing, an essential process for homologous recombination, is the formation of a heteroduplex joint by an invading single-stranded DNA tail and a complementary sequence within double-stranded DNA (dsDNA). The base rotation of the parental dsDNA, to switch from parental base pairs to heteroduplex ones with the invading single-stranded DNA, sterically requires vertical extension between adjacent base pairs, which inevitably induces untwisting of the dsDNA. RecA is a prototype of the RecA/Rad51/Dmc1 family proteins, which promote ATP-dependent homologous pairing in homologous DNA recombination in vivo, except in mitochondria. As predicted by the requirement for the untwisting, dsDNA bound to RecA is extended and untwisted, and homologous pairing by RecA in vitro is extensively stimulated by the negative supercoils of dsDNA substrates. D-loop formation in negatively supercoiled dsDNA, which serves as an assay for homologous pairing, is also catalyzed in an ATP-independent manner by proteins structurally unrelated to RecA, such as Mhr1. Mhr1 is required for yeast mitochondrial DNA recombination instead of RecA family proteins. Inconsistent with the topological requirements, tests for the effects of negative supercoils revealed that Mhr1 catalyzes homologous pairing with relaxed closed circular dsDNA much more efficiently than with negatively supercoiled dsDNA. Topological analyses indicated that neither the process nor the products of homologous pairing by Mhr1 involve a net topological change of closed circular dsDNA. This would be favorable for homologous recombination in mitochondria, where dsDNA is unlikely to be under topological stress toward unwinding. We propose a novel topological mechanism wherein Mhr1 induces untwisting without net topological change.

2D NMR study of the DNA duplex d(CTCTC*A*ACTTCC).d(GGAAGTTGAGAG) cross-linked by the antitumor-active dirhodium(II,II) unit at the cytosine-adenine step

The 2D NMR analysis in solution of the DNA duplex d(CTCTC*A*ACTTCC).d(GGAAGTTGAGAG) binding to the dirhodium unit cis-[Rh2(mu-O2CCH3)2(eta1-O2CCH3)]+ showed that an unprecedented intrastrand adduct, dsII, is formed with the dirhodium unit cross-linking in the major groove residues C5 and A6 (indicated with asterisks), also corroborated by enzyme digestion studies. Formation of the dirhodium complex dsII destabilizes significantly the duplex as indicated by the substantial decrease in its melting temperature (DeltaTm = -22.9 degrees C). The reduced thermal stability of dsII is attributed to the decreased stacking of the bases and the complete disruption and/or weakening of the hydrogen bonds within the base pairs in the immediate vicinity of the metalation site (C5.G20 and A6.T19), but the effects due to the metal binding are more severe for the base pairs in the 5' direction to the lesion site. The NMR spectroscopic data indicate that Watson-Crick hydrogen bonding is completely disrupted for the C5.G20 site and considerably weakened for A6.T19. In dsII, the bases C5 and A6 bind to eq positions of the dirhodium unit cis-[Rh2(mu-O2CCH3)2(eta1-O2CCH3)]+, which retains one monodentate and two bridging acetate groups, presumably due to steric reasons. Binding of A6 takes place via N7, whereas binding of the C5 base takes place via the exocyclic N4 site, resulting in the anti-cytosine rotamer with respect to site N3 in its metal-stabilized rare iminooxo form.

NMR analysis of duplex d(CGCGATCGCG)2 modified by Λ- and Δ-[Ru(bpy)2(m-GHK)]Cl2 and DNA photocleavage study

The interaction of the diastereomeric complexes Lambda-[Ru(bpy)2(m-GHK)]Cl2 and Delta-[Ru(bpy)2(m-GHK)]Cl2 (bpy is 2,2'-bipyridine, GHK is glycine-L-histidine-L-lysine) with the deoxynucleotide duplex d(5'-CGCGATCGCG)2 was studied by means of 1H NMR spectroscopy. At a Delta-isomer to DNA ratio of 1:1, significant shifts for the metal complex are observed, whereas there is negligible effect on the oligonucleotide protons and only one intermolecular nuclear Overhauser effect (NOE) is present at the 2D nuclear Overhauser enhancement spectroscopy spectrum. The 1Eta NMR spectrum at ratio 2:1 is characterized by a slight shift for the Delta-isomer's bpy aromatic protons as well as significant shifts for the decanucleotide G4 H1' and Eta2'', A5 H2, G10 H1', T6 NH and G2 NH protons. Furthermore, at ratio 2:1, 11 intermolecular NOEs are observed. The majority of the NOEs involve the sugar Eta2' and Eta2'' protons sited in the major groove of the decanucleotide. Increasing the Delta-isomer to d(CGCGATCGCG)2 ratio to 5:1 results in noteworthy spectral changes. The Delta-isomer's proton shifts are reduced, whereas significant shifts are observed for the decanucleotide protons, especially the sugar protons, as well as for the exchangeable protons. Interaction is characterized by the presence of only one intermolecular NOE. Furthermore, there is significant broadening of the imino proton signals as the ratio of the Delta-isomer to DNuAlpha increases, which is attributed to the opening of the two strands of the duplex. The Lambda-isomer, on the other hand, approaches the minor groove of the oligonucleotide and interacts only weakly, possibly by electrostatic interactions. Photocleavage studies were also conducted with the plasmid pUC19 and a 158-bp restriction fragment, showing that both diastereomers cleave DNA with similar efficiency, attacking mainly the guanines of the sequence probably by generating active oxygen species.

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.

Solvation of nucleosides in aqueous mixtures of organic solvents: relevance to DNA open basepairs

Toward the goal of understanding how open basepairs in DNA interact with their heterogeneous environment, we have studied the steady-state intrinsic fluorescence properties of the purine and pyrimidine deoxynucleosides in organic solvents in the presence of small amounts of water. The organic solvents used in the present study were: n-butanol, acetonitrile, methanol, n-propanol, isopropanol, and isobutanol. For n-butanol and acetonitrile, which have a high degree of amphiphilicity and weak hydrogen bonding ability, respectively, the fluorescence spectral properties of the purines are found to depend on the sequence of steps in which the aqueous mixtures were formed. By contrast, no such dependence was observed in the mixtures with any of the other solvents used in the present study. Moreover, no such dependence was observed for the pyrimidines. These findings suggest that the final solvation network around the purines is dependent on the nature of the environment to which they were initially exposed. This would tend to present an impediment to the closing of AT or GC basepairs in DNA that become open as a result of structural fluctuations, DNA bending, or protein-DNA interactions.

Alteration of A.T base-pair opening kinetics by the ammonium cation in DNA A-tracts

We have investigated by NMR the effects of NH(4)(+) on the chemical shifts, on the structure, and on the imino proton exchange kinetics of two duplexes containing an A-tract, [d(CGCGAATTCGCG)](2) and [d(GCA(4)T(4)GC)](2), and of a B-DNA duplex,[d(CGCGATCGCG)](2). Upon NH(4)(+) addition to [d(CGCGAATTCGCG)](2), the adenosine H2 protons, the thymidine imino protons, and the guanosine imino proton of the adjacent G.C pair show unambiguous chemical shifts. Similar shifts are observed in the A-tract of [d(GCA(4)T(4)GC)](2) and for the A5(H2) proton of the B DNA duplex [d(CGCGATCGCG)](2). The localization of the shifted protons suggests an effect related to NH(4)(+) binding in the minor groove. The cross-peak intensities of the NOESY spectra collected at low and high NH(4)(+) concentrations are comparable, and the COSY spectra do not show any change of the sugar pucker. This indicates a modest effect of ammonium binding on the duplex structures. Nevertheless, the imino proton exchange catalysis by ammonia provides evidence for a substantial effect of NH(4)(+) binding on the A.T base-pair kinetics in the A-tracts. Proton exchange experiments performed at high and low NH(4)(+) concentrations show the occurrence of two native conformations in proportions depending on the NH(4)(+) concentration. The base-pair lifetimes and the open-state lifetimes of each conformation are distinct. Exchange from each conformation proceeds via a single open state. But if, and only if, the NH(4)(+) concentration is kept larger than 1 M, the A.T imino proton exchange times of A-tract sequences exhibit a linear dependence versus the inverse of the NH(3) proton acceptor concentration. This had been interpreted as an indication for two distinct base-pair opening modes (Wärmländer, S., Sen, A., and Leijon, M. (2000) Biochemistry 39, 607-615).

Unrestrained 5 ns molecular dynamics simulation of a cisplatin-DNA 1,2-GG adduct provides a rationale for the NMR features and reveals increased conformational flexibility at the platinum binding site

A 5 ns unrestrained molecular dynamics (MD) simulation of the DNA duplex d(GCCG*G*ATCGC)-d(GCGATCCGGC), bearing a cis-Pt(NH(3))(2)(2+) unit crosslinking the two G* guanine bases, is reported. The MD trajectory was a posteriori correlated with NMR data determined for the same adduct, and it is shown that interproton distances and the characteristic chemical shifts are accounted for by the simulation. The simulation and its confrontation with the NMR data have confirmed the finding derived early from static models that the cytosine complementary to the 5' G*, C17, is mobile with respect to its adjacent bases. However, in contrast to our previous description of this mobility, which included rupture of the Watson-Crick hydrogen bonds and formation of non-Watson-Crick hydrogen bonds, the MD simulation indicated that the G*4-C17 pair moves continuously along a trajectory roughly perpendicular to the local helix axis, with retention of all three Watson-Crick hydrogen bonds. The simulation indicated the reversible formation of a hydrogen bond between the 5' oriented NH(3) ligand of platinum and the C3pG*4 phosphate group, in accord with our former prediction. Furthermore, the simulation has disclosed previously undetected BI <=> BII transitions at the G*5pA6 and A6pT7 steps, connected to formation/rupture of a hydrogen bond between the 3' oriented NH(3) ligand of platinum and the N7 atom of A6. All these conformational equilibria affect the form of the minor groove and increase the conformational flexibility at the platination site, and are thus likely to facilitate recognition by cellular proteins.

Imino proton exchange and base-pair kinetics in RNA duplexes

Using NMR magnetization transfer from water and ammonia-catalyzed exchange of the imino proton, we have measured the base-pair lifetimes and the dissociation constants of six RNA duplexes: [r(CGCGAUCGCG)](2), [r(CGCGAAUUCGCG)](2), [r(CCUUUCGAAAGG)](2), [r(CGCACGUGCG)](2), [r(GGU(8)CC).r(GGA(8)CC)], and [poly(rA).poly(rU)], and we compare them with those of their DNA homologues. As predicted by a two-state (closed/open) model of the pair, the imino proton exchange times decrease linearly vs. the inverse of catalyst concentration. As in DNA duplexes, base pairs open one at a time, and the kinetics is in most cases insensitive to the nature of the adjacent residues. The lifetime of the r(G.C) pairs, 40 to 50 ms, is longer than that of the equivalent in the corresponding oligodeoxynucleotides, and the dissociation constants, about 10(-)(7), are slightly smaller. The r(A.U) opening and closing rates are much larger than those of the d(A.T) pairs, but the stabilities are comparable.

Threading bis-intercalation of a macrocyclic bisacridine at abasic sites in DNA: nuclear magnetic resonance and molecular modeling study

The macrocyclic bisacridine (CBA) has been reported previously to specifically recognize single-stranded nucleic acid structures, especially DNA hairpins. The binding of the drug with an abasic site-containing oligonucleotide, was investigated by (1)H NMR and molecular modeling. We have used a DNA undecamer, the d(C(1)G(2)C(3)A(4)C(5)X(6)C(7)A(8)C(9)G(10)C(11)) x d(G(12)C(13)G(14)T(15)G(16)T(17)G(18)T(19)G(2)(0)C(21)G(22)) duplex in which the X residue is a stable analogue of the abasic site [3-hydroxy-2-(hydroxymethyl) tetrahydrofuran]. Analysis of the NMR data reveals that the bisacridine molecule forms two different intercalation complexes in a 80/20 (+/- 10) ratio. For the major complex, a molecular modeling study was performed guided by nineteen intermolecular drug-DNA restraints, determined from NOESY spectra. In this model, the ligand interacts in the threading binding mode with an acridine ring intercalated between the C(7)-A(8) and T(15)-G(16) base pairs, while the other acridine ring resides in the abasic pocket. The two linker chains are positioned in the minor and in the major groove, respectively. A comparable study was performed to evaluate the interaction of CBA with the parent unmodified duplex in which X(6) was replaced by an adenine residue. No complex formation was observed when operating in identical conditions. This shows the selective binding of CBA to the abasic site and its potential interest to target the abasic site lesion.

Decreased imino proton exchange and base-pair opening in the IHF-DNA complex measured by NMR

Integration Host Factor, IHF, is an E. coli DNA binding protein that imposes a substantial bend on DNA. Previous footprinting studies and bending assays have characterized several recognition sequences in the bacterial and lambda phage genome as unique in the way they are bound by IHF. We have chosen one of the lambda phage sites, H1, for study because it presents a small yet sequence-specific substrate for NMR analysis of the complex. A 19 base-pair duplex, H19, corresponding to the recognition sequence at the H1 site was constructed by isotopically labeling one of the strands with 15N. (1H, 15N) heteronuclear NMR experiments aided in assigning the imino proton resonances of the DNA alone and in complex with IHF. The NMR results are consistent with a mode of binding observed in the recent crystal structure of IHF bound to another of its sites from the lambda phage genome. Additionally, the dramatic change that IHF imposes on the imino proton chemical shifts is indicative of a severe deviation from canonical B-DNA structure. In order to understand the dynamic properties of the DNA in the complex with IHF, the exchange rates of the imino protons with the solvent have been measured for H19 with and without IHF bound. A drastic reduction in exchange is observed for the imino protons in the IHF bound DNA. In the DNA-protein complex, groups of adjacent base-pair exchange at the same rate, and appear to close more slowly than the rate of imino proton exchange with bulk water, since their exchange rate is independent of catalyst concentration. We infer that segments of the double helix as large as 6 bp open in a cooperative process, and remain open much longer than is typical for opening fluctuations in naked duplex DNA. We discuss these results in terms of the specific protein-DNA contacts observed in the crystal structure.

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.

Dynamics of a bis-intercalator DNA complex by 1H-detected natural abundance 13C NMR spectroscopy

The dynamics of the DNA oligomer d(CGCTAGCG)2 (CTSYM) and its complex with the dye 1,1-(4,4,8,8-tetramethyl-4, 8-diazaundecamethylene)-bis-4-(3-methyl-2,3-dihydro-(benzo-1, 3-thiazole)-2-methylidene)-quinolinium tetraiodide (TOTO) (CTSYMTOTO) bis-intercalated at the 5'-CT-3' sequence steps have been determined from NMR relaxation parameters. Longitudinal and transverse 13C relaxation rates and heteronuclear NOE relaxation data were acquired and have been analyzed in the context of the Lipari and Szabo model-free formalism. The overall rotational correlation time for the CTSYM is 3.44 ns and the CTSYMTOTO is 3.48 ns. The generalized order parameters (S2) for methine carbons in the CTSYM and CTSYMTOTO are relatively high but nonuniform for the molecules and show sequence context and conformation-dependent variations. Average values of S2 = 0.79 +/- 0.02 for the CTSYM, S2 = 0.80 +/- 0.04 for the CTSYMTOTO aromatic spins, S2 = 0.76 +/- 0.02 for the CTSYM, and S2 = 0.83 +/- 0.05 for the CTSYMTOTO deoxyribose spins were found. The S2 values for the 5'-terminal deoxyribose are lower than for the other residues. The DNA backbone in CTSYMTOTO is distorted and elongated at the site of intercalation, and the C3' atom of the C3 deoxyribose residue has a very low S2 = 0.57 +/- 0.06. The low order for this spin is interpreted in terms of exchange between the C2'-endo and O1'-endo conformations of the C3 deoxyribose. Significant chemical exchange processes were found for most of the aromatic spins in CTSYM that are interpreted in terms of microsecond to millisecond time scale dynamics. The microsecond to millisecond dynamics of the bases in CTSYM are quenched upon TOTO complex formation due to unwinding of the helix and an increase in the surface area of the bases in mutual contact and the large surface area in contact with the intercalated dye. The derived order parameters combined with the solution structure provide motional models for conformational changes induced in the backbone in response to the ligand binding.

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.

Dynamics in psoralen-damaged DNA by 1H-detected natural abundance 13C NMR spectroscopy

The dynamics of the DNA oligomer d(GCGTACGC)2 and the 4'-(hydroxymethyl)-4,5',8-trimethylpsoralen-DNA furanside monoadduct (MAf) of this oligomer have been determined from NMR relaxation parameters. Longitudinal and transverse 13C relaxation rates and heteronuclear NOE relaxation data have been measured at natural abundance and have been analyzed in the context of the Lipari and Szabo model-free formalism. The generalized order parameters for methine carbons in the octamer sequence d(GCGTACGC)2 (UM) are relatively and uniformly high for the entire molecule. The generalized order parameters for methine carbons in the MAf are significantly lower for the deoxyribose bearing the damaged thymidine base and for the bases flanking the lesion on the undamaged strand, indicating additional conformational flexibility due to the lesion. The order parameters for the bases on the damaged strand flanking the lesion remain high. Analysis of the relaxation data indicates substantial chemical exchange for the adenosine residues in the UM TpA site, and this chemical exchange is quenched upon MAf formation. These data are discussed in terms of a model for DNA damage recognition by the nucleotide excision repair system.

Bending and circularization of site-specific and stereoisomeric carcinogen-DNA adducts

The potent tumorigen and mutagen (+)-7(R),8(S)-dihydroxy-9(S), 10(R)-epoxy-7,8,9,10-tetrahydrobenzo[a]pyrene ((+)-anti-BPDE) is a metabolite of benzo[a]pyrene that binds predominantly to the exocyclic amino group of guanine residues in DNA in vivo and in vitro. While the (-)-7S,8R,9R,10Senantiomer, (-)-anti-BPDE, also reacts with DNA to form similar covalent N2-deoxyguanosyl adducts, this diol epoxide is nontumorigenic and its mutagenic activities are different from those of (+)-anti-BPDE. In this work, T4 ligase-induced cyclization methods have been employed to demonstrate that the (+)-anti-[BP]-N2-dG lesions (G*) cause significantly greater amounts of bending and circularization of the one-base overhang undecamer duplex 5'-d(CACAT[G*]TACAC).d(TGTACATGTGG) than the stereoisomeric oligonucleotide duplex with G* = (-)-anti-[BP]-N2-dG. In the case of the (+)-anti-BPDE-modified oligonucleotides, the ratio of circular to linear DNA multimers reaches values of 8-9 for circle contour sizes of 99-121 base pairs, while for the (-)-anti-[BP]-N2-dG-modified DNA this ratio reaches a maximum value of only approximately 1 at 154-176 base pairs. Assuming a planar circle DNA model, the inferred bending angles for 90-92% of the observed circular ligation products range from 30 to 51 degrees per (+)-trans-anti-[BP]-N2-dG lesion and from 20 to 40 degrees per (-)-trans-anti-[BP]-N2-dG lesion. In the case of unmodified DNA, the probability of circular product formation is at least 1 order of magnitude less efficient than in the BPDE-modified sequences and about 90% of the circular products exhibit bending angles in the range of 14 -19 degrees . In the most abundant circular products observed experimentally, the bending angles are 40 degrees and 26 +/- 2 degrees per (+)-anti-[BP]- or (-)-anti-[BP]-modified 11-mer; these values correspond to a net contribution of 21-26 degrees and 5-19 degrees , respectively, to the observed overall bending per lesion. The coexistence of circular DNA molecules of different sizes and, therefore, different average bending angles per lesion, suggest that the lesions induce both torsional flexibility and flexible bends, which permit efficient cyclization, especially in the case of (+)-trans-[BP]-N2-dG adducts. The NMR characteristics of (+)-trans-[BP]-N2-dG lesion in the 11-mer duplex 5'-d(CACAT[G*]TACAC).d(GTGTACATGTG) indicate that all base pairs are intact, except at the underlined base pairs. This suggests a distortion in the normal conformation of the duplex on the 5'-side of the modified guanosine residue, which may be due to bending enhanced base pair opening and bending induced by the bulky carcinogen residue. The implications of base sequence-dependent flexibilities and conformational mobilities of anti-[BP]-N2-dG lesions on DNA replication and mutation are discussed.

Internal motions of nucleic acid structures and the determination of base-pair lifetimes

Over the years, imino proton exchange measurements have provided a description of base-pair opening and of properties of internal motions of nucleic acids such as the individual opening of base pairs and the insensitivity of base-pair lifetimes to the nature and stability of neighboring pairs. A recent determination of base-pair lifetimes in d(CGCGATCGCG) conflicts with the original measurements and their interpretation. This question is analyzed in the present work. We emphasize the importance of high concentrations of exchange catalyst (eg 1 M NH3) for the accurate determination of base-pair lifetimes. These concentrations entail a high ionic strength, which can lead to aggregation, enhanced magnetic relaxation and underestimation of base-pair lifetimes if exchange is measured by its effect on the proton relaxation rate. Magnetization transfer which provides a more direct method for the measurement of proton exchange rates is therefore preferred. We show that the lack of measurements at high catalyst concentration is responsible for the discrepancy mentioned above. Measurements by the magnetization transfer method validate the original interpretation and inversion recovery experiments illustrate the effect of the ionic strength on the relaxation rate of the imino protons.

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.

DNA structural forms

In the years that have passed since the publication of Wolfram Saenger's classic book on nucleic acid structure (Saenger, 1984), a considerable amount of new data has been accumulated on the range of conformations which can be adopted by DNA. Many unusual species have joined the DNA zoo, including new varieties of two, three and four stranded helices. Much has been learnt about intrinsic DNA curvature, dynamics and conformational transitions and many types of damaged or deformed DNA have been investigated. In this article, we will try to summarise this progress, pointing out the scope of the various experimental techniques used to study DNA structure, and, where possible, trying to discern the rules which govern the behaviour of this subtle macromolecule. The article is divided into six major sections which begin with a general discussion of DNA structure and then present successively, B-DNA, DNA deformations, A-DNA, Z-DNA and DNARNA hybrids. An extensive set of references is included and should serve the reader who wishes to delve into greater detai.

Base-catalysis of imino proton exchange in DNA: effects of catalyst upon DNA structure and dynamics

Characterization of the kinetics and energetics of base-pair opening in nucleic acids relies upon measurements of the rates of exchange of imino protons with water protons at high concentrations of the exchange catalyst. Under these conditions, the exchange catalyst may affect structural or dynamic properties of the nucleic acid molecule and thus, limit the significance of the exchange data. To address this problem, we have used NMR spectroscopy to characterize the effects of a catalyst of imino proton exchange, namely, ammonia upon the structure and dynamics of the self-complementary DNA dodecamer [d(CGCAGATCTGCG)]2. The changes in structure were monitored in proton NOESY and DQF-COSY experiments and in phosphorus spectra at 15 degrees C and at ammonia concentrations ranging from 0.002 to 0.5 M. The results indicate that ammonia induces subtle changes in the solution conformation of the dodecamer, but the overall structure is maintained close to the B-type DNA structure. However, the relaxation rates (i.e., transverse, longitudinal, and cross-relaxation rates) of several non-exchangeable protons were found to increase by approximately 50% upon changing ammonia concentration from 0.002 to 0.5 M. The increases were comparable for all protons investigated suggesting that they originate from an ammonia-induced increase in the overall correlation time of the DNA dodecamer. Numerical analysis revealed that the catalyst-induced enhancements in proton relaxation can alter significantly the calculated values of the exchange rates of imino protons, especially those obtained from measurements of the line widths of these proton resonances.

Rapid estimation of relative amide proton exchange rates of 15 N-labelled proteins by a straightforward water selective NOESY-HSQC experiment

A straightforward heteronuclear pseudo-3D NOESY-HSQC pulse sequence using radiation damping to selectively invert magnetization at the water frequency was developed to estimate the amide proton exchange rates in 15N-labelled proteins. The peak intensities in the resultant 2D spectrum allow a direct classification of amide proton exchange rates according to short (ms), intermediate (ms to s) or long (> or = s) residence times. This method was successfully used for the analysis of amide proton exchange rates in the 15N-labelled FruR DNA-binding domain and pertinent information about its dynamics was obtained.

Base-pair exchange kinetics of the imino and amino protons of the 3'-phenazinium tethered DNA-RNA duplex, r(5'GAUUGAA3'):d(5'TCAATC3'-Pzn), and their comparison with those of B-DNA duplex

The dynamics of the opening-closing of the constituent base-pairs as well as of the exchange kinetics of the base-paired imino and amino protons with water in a DNA-RNA hybrid, [5'r(G1A2U3U4G5A6A7)3']:5'p[d(T8C9A10A11T12C13)]3'-Pzn ] duplex (I), are reported here in details for the first time. The exchange kinetics of amino and imino protons in the DNA-RNA hybrid (duplex I) have been compared with identical studies on the following B-DNA duplexes: d(C1G2T3A4C5G6)2 (II), d[p(5'T1G2T3T4T5G6G7C8)3']:d[p(5'C9C10A11A12A13C14A15)3'] (III), d(C5G6C7G8A9A10T11T12C13G14C15G16)2 (IV) and d(C1G2C3G4C5G6C7G8A9A10T11T12C13G14C15G16C17G18C19G20)2 (V). This comparative study shows that the life-times tau o of various base-pairs in the DNA-RNA hybrid (I) varies in the range of approximately 1 ms, and they are quite comparable to those of the shorter B-DNA duplexes (II) and (III), but very different from the tau o of the larger duplexes (IV) and (V): the tau o for the base pair of T11 and T12 residues in the 20-mer (duplex V) are 2.9 +/- 2.3 ms and 23.2 +/- 8.9 ms, respectively, while the corresponding tau o in the 12-mer (duplex IV) are 2.8 +/- 2.2 ms and 17.4 +/- 5.4 ms. It has also been shown that the total energy of activation (Ea) assessed from the exchange rates of both imino and amino protons, representing energetic contributions from both base-pair and helix opening-closing as well as from the exchange process of the imino protons from the open state with the bound water, is close to the Ea of the short B-DNA duplex (Ea approximately 28-47 kcal/mol).

Effect of nucleoside modifications on the structure and thermal stability of Escherichia coli valine tRNA

Transfer RNA transcribed in vitro lacks the base modifications found in native tRNA. To understand the effect of base modifications on the structure of tRNA, the downfield region of the 1H NMR spectrum of in vitro transcribed E coli tRNAVal in aqueous phosphate buffer in the presence of excess Mg2+ was investigated. The resonances of all imino protons involved in hydrogen bonds in the helical stem regions and in tertiary interactions were assigned using two-dimensional nuclear Overhauser enhancement spectroscopy (NOESY) and one-dimensional difference nuclear Overhauser effect (NOE) methods. In addition, some aromatic C2 and C8 proton resonances as well as one amino proton resonance were assigned. The chemical shifts of the assigned resonances of unmodified E coli tRNAVal were compared with those of the native tRNA molecule under similar solution conditions. The similarity of the NMR data for unmodified and modified tRNA indicates that the in vitro transcribed tRNA has nearly the same solution structure as the native molecule in the presence of excess Mg2+. The only significant differences were the chemical shifts of resonances corresponding to protons in (or interacting with) bases, indicating the possibility of local structural perturbations. The thermal stability of E coli modified and unmodified tRNAVal in the presence of Mg2+ was also investigated by analyzing the temperature dependence of the imino proton spectra. Several tertiary interactions involving modified nucleosides in native E coli tRNAVal are less stable in the absence of base modifications.

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)

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.

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.