Solution structure of the EcoRI DNA sequence: refinement of NMR-derived distance geometry structures by NOESY spectrum back-calculations

Single-Base Lesions and Mismatches Alter the Backbone Conformational Dynamics in DNA

DNA damage has been implicated in numerous human diseases, particularly cancer, and the aging process. Single-base lesions and mismatches in DNA can be cytotoxic or mutagenic and are recognized by a DNA glycosylase during the process of base excision repair. Altered local dynamics and conformational properties in damaged DNAs have previously been suggested to assist in recognition and specificity. Herein, we use solution nuclear magnetic resonance to quantify changes in BI-BII backbone conformational dynamics due to the presence of single-base lesions in DNA, including uracil, dihydrouracil, 1,N6-ethenoadenine, and T:G mismatches. Stepwise changes to the %BII and ΔG of the BI-BII dynamic equilibrium compared to those of unmodified sequences were observed. Additionally, the equilibrium skews toward endothermicity for the phosphates nearest the lesion/mismatched base pair. Finally, the phosphates with the greatest alterations correlate with those most relevant to the repair of enzyme binding. All of these results suggest local conformational rearrangement of the DNA backbone may play a role in lesion recognition by repair enzymes.

Pillar[5]arene-based [1]rotaxane: high-yield synthesis, characterization and application in Knoevenagel reaction

Quantitative synthesis of a pillararene-based [1]rotaxane has been achieved via a “self-threading-stoppering” approach, followed by its first application in organic catalysis.

An FTIR investigation of flanking sequence effects on the structure and flexibility of DNA binding sites

Fourier transform infrared (FTIR) spectroscopy and a library of FTIR marker bands have been used to examine the structure and relative flexibilities conferred by different flanking sequences on the EcoRI binding site. This approach allowed us to examine unique peaks and subtle changes in the spectra of d(AAAGAATTCTTT)(2), d(TTCGAATTCGAA)(2), and d(CGCGAATTCGCG)(2) and thereby identify local changes in base pairing, base stacking, backbone conformation, glycosidic bond rotation, and sugar puckering in the studied sequences. The changes in flanking sequences induce differences in the sugar puckers, glycosidic bond rotation, and backbone conformations. Varying levels of local flexibility are observed within the sequences in agreement with previous biological activity assays. The results also provide supporting evidence for the presence of a splay in the G(4)-C(9) base pair of the EcoRI binding site and a potential pocket of flexibility at the G(4) cleavage site that have been proposed in the literature. In sum, we have demonstrated that FTIR is a powerful methodology for studying the effect of flanking sequences on DNA structure and flexibility, for it can provide information about the local structure of the nucleic acid and the overall relative flexibilities conferred by different flanking sequences.

Impact of pyrophosphate and O-ethyl-substituted pyrophosphate groups on DNA structure

Design of novel DNA probes to inhibit specific repair pathways is important for basic science applications and for use as therapeutic agents. As shown previously, single pyrophosphate (PP) and O-ethyl-substituted pyrophosphate (SPP) modifications can inhibit the DNA glycosylase activities on damaged DNA. To understand the structural basis of this inhibition, the influence of the PP and SPP internucleotide groups on the helical parameters and geometry of a double-stranded DNA was studied by using molecular modeling tools including molecular dynamics and quantum mechanical-molecular mechanical (QM/MM) approaches. Native and locally modified PP- and SPP-containing DNA duplexes of dodecanucleotide d(C1G2C3G4A5A6T7T8C9G10C11G12) were simulated in aqueous solution. The energies and forces were computed by using the PBE0/6-31+G** approach in the QM part and the AMBER force-field parameters in the MM part. Analysis of the local base-pair helical parameters, internucleotide distances, and overall global structure at the located stationary points revealed a close similarity of the initial and modified duplexes, with only torsion angles of the main chain being altered in the vicinity of introduced chemical modification. Results show that the PP and SPP groups are built into a helix structure without elongation of the internucleotide distance due to flipping-out of phosphate group from the sugar-phosphate backbone. The mechanism of such embedding has only a minor impact on the base pairs stacking and Watson-Crick interactions. Biochemical studies revealed that the PP and SPP groups immediately 5', but not 3', to the 8-oxoguanosine (8oxodG) inhibit translesion synthesis by a DNA polymerase in vitro. These results suggest that subtle perturbations of the DNA backbone conformation influence processing of base lesions.

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.

Direct observation of dipolar couplings between distant protons in weakly aligned nucleic acids

Under conditions where macromolecules are aligned very weakly with respect to an external magnetic field, Brownian diffusion no longer averages internuclear dipole-dipole interactions to zero. The resulting residual dipolar coupling, although typically 3 orders of magnitude weaker than in a fully aligned sample, can readily be measured by solution NMR methods. To date, application of this idea has focused primarily on pairs of nuclei separated by one or two covalent bonds, where the internuclear separation is known and the measured dipolar coupling provides direct information on the orientation of the internuclear vector. A method is described that allows observation of dipolar interactions over much larger distances. By decoupling nearest-neighbor interactions, it is readily possible to observe direct dipolar interactions between protons separated by up to 12 A. The approach is demonstrated for the DNA dodecamer d(CGCGAATTCGCG)2, where direct interactions are observed between protons up to three base pairs apart.

Assessment of the molecular dynamics structure of DNA in solution based on calculated and observed NMR NOESY volumes and dihedral angles from scalar coupling constants

To assess the accuracy of the molecular dynamics (MD) models of nucleic acids, a detailed comparison between MD-calculated and NMR-observed indices of the dynamical structure of DNA in solution has been carried out. The specific focus of our comparison is the oligonucleotide duplex, d(CGCGAATTCGCG)(2), for which considerable structural data have been obtained from crystallography and NMR spectroscopy. An MD model for the structure of d(CGCGAATTCGCG)(2) in solution, based on the AMBER force field, has been extended with a 14 ns trajectory. New NMR data for this sequence have been obtained in order to allow a detailed and critical comparison between the calculated and observed parameters. Observable two-dimensional (2D) nuclear Overhauser effect spectroscopy (NOESY) volumes and scalar coupling constants were back-calculated from the MD trajectory and compared with the corresponding NMR data. The comparison of these results indicate that the MD model is in generally good agreement with the NMR data, and shows closer accord with experiment than back-calculations based on the crystal structure of d(CGCGAATTCGCG)(2) or the canonical A or B forms of the sequence. The NMR parameters are not particularly sensitive to the known deficiency in the AMBER MD model, which is a tendency toward undertwisting of the double helix when the parm.94 force field is used. The MD results are also compared with a new determination of the solution structure of d(CGCGAATTCGCG)(2) using NMR dipolar coupling data.

Modeling furanose ring dynamics in DNA

Determination of the conformational flexibility of the furanose ring is of vital importance in understanding the structure of DNA. In this work we have applied a model of furanose ring motion to the analysis of deuterium line shape data obtained from sugar rings in solid hydrated DNA. The model describes the angular trajectories of the atoms in the furanose ring in terms of pseudorotation puckering amplitude (q) and the pseudorotation puckering phase phi. Fixing q, the motion is thus treated as Brownian diffusion through an angular-dependent potential U(phi). We have simulated numerous line shapes varying the adjustable parameters, including the diffusion coefficient D, pseudorotation puckering amplitude q, and the form of the potential U(phi). We have used several forms of the potential, including equal double-well potentials, unequal double-well potentials, and a potential truncated to "second order" in the Fourier series. To date, we have obtained best simulations for both equilibrium and nonequilibrium (partially relaxed) solid-state deuterium NMR line shapes for the sample [2' '-2H]-2'-deoxycytidine at the position C3 (underlined) in the DNA sequence [d(CGCGAATTCGCG)]2, using a double-well potential with an equal barrier height of U(0) = 5.5k(B)T ( approximately 3.3 kcal/mol), a puckering amplitude of q = 0.4 A, and a diffusion coefficient characterizing the underlying stochastic jump rate D = 9.9 x 10(8) Hz. Then the rate of flux for the C-D bond over the barrier, i.e., the escape velocity or the overall rate of puckering between modes, was found to be 0.7 x 10(7) Hz.

The first example of a Hoogsteen base-paired DNA duplex in dynamic equilibrium with a Watson-Crick base-paired duplex--a structural (NMR), kinetic and thermodynamic study

A single-point substitution of the O4' oxygen by a CH2 group at the sugar residue of A6 (i.e. 2'-deoxyaristeromycin moiety) in a self-complementary DNA duplex, 5'-d(C1G2C3G4A5A6T7T8C9G10C11G12)2(-3), has been shown to steer the fully Watson-Crick basepaired DNA duplex (1A), akin to the native counterpart, to a doubly A6:T7 Hoogsteen basepaired (1B) B-type DNA duplex, resulting in a dynamic equilibrium of (1A)<==>(1B): Keq = k1/k(-1) = 0.56+/-0.08. The dynamic conversion of the fully Watson-Crick basepaired (1A) to the partly Hoogsteen basepaired (1B) structure is marginally kinetically and thermodynamically disfavoured [k1 (298K) = 3.9 0.8 sec(-1); deltaHdegrees++ = 164+/-14 kJ/mol; -TdeltaS degrees++ (298K) = -92 kJ/mol giving a deltaG degrees++ 298 of 72 kJ/mol. Ea (k1) = 167 14 kJ/mol] compared to the reverse conversion of the Hoogsteen (1B) to the Watson-Crick (1A) structure [k-1 (298K) = 7.0 0.6 sec-1, deltaH degrees++ = 153 13 kJ/mol; -TdeltaSdegrees++ (298K) = -82 kJ/mol giving a deltaGdegrees++(298) of 71 kJ/mol. Ea (k-1) = 155 13 kJ/mol]. Acomparison of deltaGdegrees++(298) of the forward (k1) and backward (k-1) conversions, (1A)<==>(1B), shows that there is ca 1 kJ/mol preference for the Watson-Crick (1A) over the double Hoogsteen basepaired (1B) DNA duplex, thus giving an equilibrium ratio of almost 2:1 in favour of the fully Watson-Crick basepaired duplex. The chemical environments of the two interconverting DNA duplexes are very different as evident from their widely separated sets of chemical shifts connected by temperature-dependent exchange peaks in the NOESY and ROESY spectra. The fully Watson-Crick basepaired structure (1A) is based on a total of 127 intra, 97 inter and 17 cross-strand distance constraints per strand, whereas the double A6:T7 Hoogsteen basepaired (1B) structure is based on 114 intra, 92 inter and 15 cross-strand distance constraints, giving an average of 22 and 20 NOE distance constraints per residue and strand, respectively. In addition, 55 NMR-derived backbone dihedral constraints per strand were used for both structures. The main effect of the Hoogsteen basepairs in (1B) on the overall structure is a narrowing of the minor groove and a corresponding widening of the major groove. The Hoogsteen basepairing at the central A6:T7 basepairs in (1B) has enforced a syn conformation on the glycosyl torsion of the 2'-deoxyaristeromycin moiety, A6, as a result of substitution of the endocyclic 4'-oxygen in the natural sugar with a methylene group in A6. A comparison of the Watson-Crick basepaired duplex (1A) to the Hoogsteen basepaired duplex (1B) shows that only a few changes, mainly in alpha, sigma and gamma torsions, in the sugar-phosphate backbone seem to be necessary to accommodate the Hoogsteen basepair.

Absence of minor groove monovalent cations in the crosslinked dodecamer C-G-C-G-A-A-T-T-C-G-C-G

The structure of a crosslinked B -DNA dodecamer of sequence C-G-C-G-A-A-T-T-C-G-C-G has been solved to a resolution of 1.43 A. The dithiobis-propane crosslink, -CH2-CH2-CH2-S-S-CH2-CH2-CH2-, bridges N7 atoms of adenine bases 6 and 18 in the two central base-pairs within the major groove. The crosslink is sufficiently long that no bending is induced in the helix, which is essentially isostructural with the native unlinked dodecamer at 1.9 A. A constellation of solvent peaks tentatively fitted as a spermine molecule in that earlier analysis is now seen at higher resolution to be a well-defined octahedral magnesium hexahydrate complex in the major groove. One end of the duplex curves around that complex to produce a roll-bend near base-pairs 3-5, and an overall bend in helix axis, as has long been noted. Two other magnesium complexes connect the helices and help to knit the crystal lattice together. No evidence exists for partial sodium or potassium ion substitution for solvent water molecules within the minor groove spine of hydration, as had been suggested previously: not coordination geometry and environment, nor B values, nor calculated valence values, nor difference map analyses. Indeed, the very numbers that have been claimed in support of partial substitution by sodium or potassium ions are reproduced with the present crystals, which by chemical analysis contains only one trace sodium ion per 160 bp, and one potassium ion per 41 bp. In contrast, our crystals contain one Mg2+ per base-pair, meaning that phosphate group charge neutrality is accomplished by divalent cations, not monovalent ions. Three of these magnesium cations per duplex are localized and visible in the X-ray analysis, and nine are disordered and invisible. Hence although binding of monovalent cations within the minor groove of A -tracts on occasion may be a consequence of groove narrowing, it cannot be the cause of that narrowing. Cations, contrary to what has been claimed, are not in charge.

The solution conformation of a carbocyclic analog of the Dickerson-Drew dodecamer: comparison with its own X-ray structure and that of the NMR structure of the native counterpart

The NMR conformation of a carbocyclic analog of the Dickerson-Drew dodecamer [d(CGC-GAAT*T*CGCG)]2 containing 6'-alpha-Me carbocyclic thymidines (T*) has been determined and compared with that of its X-ray structure. The solution structure of the 6'-alpha-Me carbocyclic thymidine modified duplex has also been compared with the solution structure of the corresponding unmodified Dickerson-Drew duplex solved by us under the same experimental conditions. The NMR structures have been based on 24 experimental distance and torsion constraints per residue for [d(CGCGAAT*T*CGCG)]2 (1) and on 21 constraints per residue for the natural counterpart. In general, both final NMR structures are more close to the B-type DNA. The cyclopentane moieties of the carbocyclic thymidine residues adopt C1'-exo B-DNA type puckers (the phase angles P = 136-139 degrees and the puckering amplitudes psi = 36-37 degrees) that are close to their previously published crystal C1'-exo or C2'-endo puckers. The main differences between the two NMR structures are for beta(T*8) and epsilon, xi(T*7) backbone torsions (27-50 degrees ), for basepair twist for the 7-8 and 8-9 basepair steps (5-6 degrees), tilt for the 8-9 step (7 degrees), roll for the 7-8 step (7 degrees), shift for the 7-8 step (0.9A) and slide for the 9-10 step (0.6A). The relatively small deviations of helical structure parameters lead to structural isomorphism of these duplexes in aqueous solutions (atomic RMSD = 1.0A). The difference of the minor groove widths (less than 0.7A) in the core part of the modified duplex in comparison with the native one is much smaller than the difference between the X-ray structures of these duplexes. A detailed comparison of NMR and X-ray structure parameters showed significant monotonic differences (0.9-2.5A) for all basepair slides in both duplexes. Deviations between NMR and X-ray structure parameters for the modified duplex were also found for basepair tilt of the 4-5 step (13 degrees), rolls for the 8-9 and 10-11 steps (16 and 19 degrees), twist of the 3-4 step (8 degrees) and shift of the 9-10 step (0.9A).

The residence time of the bound water in the hydrophobic minor groove of the carbocyclic-nucleoside analogs of Dickerson-Drew dodecamers

The residence time of the bound water molecules in the antisense oligodeoxyribonucleotides containing 7'-alpha-methyl (TMe) carbocyclic thymidines in duplex (I), d5'(1C2G3C4G5A6A7TMc8TMe9C10G11C12G)23', and 6'-alpha-hydroxy (TOH) carbocyclic thymidines in duplex (II), d5'(1C2G3C4G5AOH6AOH7TOH8 TOH9C10G11C12G)2(3), have been investigated using a combination of NOESY and ROESY experiments. Because of the presence of 7'-alpha-methyl groups of TMe in the centre of the minor groove in duplex (I), the residence time of the bound water molecule is shorter than 0.3 ns. The dramatic reduction of the residence time of the water molecule in the minor groove in duplex (I) compared with the natural counterpart has been attributed to the replacement of second shell of hydration and disruption of hydrogen-bonding with 04' in the minor groove by hydrophobic alpha-methyl groups, as originally observed in the X-ray study. This effect could not be attributed to the change of the width of the minor groove because a comparative NMR study of the duplex (I) and its natural counterpart showed that the widths of their minor grooves are more or less unchanged (r.m.s.d change in the core part is <0.63A). For duplex (II) with polar 6'-alpha-hydroxyl groups pointed to the minor groove, the correlation time is much longer than 0.36 ns as a result of the stabilising hydrogen-bonding interaction with N3 or 02 of the neighbouring nucleotides.

Exploratory studies on azole carboxamides as nucleobase analogs: thermal denaturation studies on oligodeoxyribonucleotide duplexes containing pyrrole-3-carboxamide

In order to study base pairing properties of the amide group in DNA duplexes, a nucleoside analog, 1-(2'-deoxy-beta-D-ribofuranosyl)pyrrole-3-carboxamide, was synthesized by a new route from the ester, methyl 1-(2'-deoxy-3',5'-di-O-p -toluoyl-beta-D-erythro-pentofuranosyl)pyrrole-3-carboxylate, obtained from the coupling reaction between 1-chloro-2-deoxy-3,5-di-O -toluoyl-d-erythropentofuranose and methyl pyrrole-3-carboxylate by treatment with dimethylaluminum amide. 1-(2'-Deoxy-beta-D-ribofuranosyl)pyrrole-3-carboxamide was incorporated into a series of oligodeoxyribonucleotides by solid-phase phosphoramidite technology. The corresponding oligodeoxyribonucleotides with 3-nitropyrrole in the same position in the sequence were synthesized for UV comparison of helix-coil transitions. The thermal melting studies indicate that pyrrole-3-carboxamide, which could conceptually adopt either a dA-like or a dI-like hydrogen bond conformation, pairs with significantly higher affinity to T than to dC. Pyrrole-3-carboxamide further resembles dA in the relative order of its base pairing preferences (T >dG >dA >dC). Theoretical calculations on the model compound N-methylpyrrole-3-carboxamide using density functional theory show little difference in the preference for a syntau versus anti conformation about the bond from pyrrole C3 to the amide carbonyl. The amide groups in both the minimized antitau and syntau conformations are twisted out of the plane of the pyrrole ring by 6-14 degrees. This twist may be one source of destabilization when the amide group is placed in the helix. Another contribution to the difference in stability between the base pairs of pyrrole-3-carboxamide with T and pyrrole-3-carboxamide with C may be the presence of a hydrogen bond in the former involving an acidic proton (N3-H of T).

A 5-nanosecond molecular dynamics trajectory for B-DNA: analysis of structure, motions, and solvation

We report the results of four new molecular dynamics (MD) simulations on the DNA duplex of sequence d(CGCGAATTCGCG)2, including explicit consideration of solvent water, and a sufficient number of Na+ counterions to provide electroneutrality to the system. Our simulations are configured particularly to characterize the latest MD models of DNA, and to provide a basis for examining the sensitivity of MD results to the treatment of boundary conditions, electrostatics, initial placement of solvent, and run lengths. The trajectories employ the AMBER 4.1 force field. The simulations use particle mesh Ewald summation for boundary conditions, and range in length from 500 ps to 5.0 ns. Analysis of the results is carried out by means of time series for conformationalm, helicoidal parameters, newly developed indices of DNA axis bending, and groove widths. The results support a dynamically stable model of B-DNA for d(CGCGAATTCGCG)2 over the entire length of the trajectory. The MD results are compared with corresponding crystallographic and NMR studies on the d(CGCGAATTCGCG)2 duplex, and placed in the context of observed behavior of B-DNA by comparisons with the complete crystallographic data base of B-form structures. The calculated distributions of mobile solvent molecules, both water and counterions, are displayed. The calculated solvent structure of the primary solvation shell is compared with the location of ordered solvent positions in the corresponding crystal structure. The results indicate that ordered solvent positions in crystals are roughly twice as structured as bulk water. Detailed analysis of the solvent dynamics reveals evidence of the incorporation of ions in the primary solvation of the minor groove B-form DNA. The idea of localized complexation of otherwise mobile counterions in electronegative pockets in the grooves of DNA helices introduces an additional source of sequence-dependent effects on local conformational, helicoidal, and morphological structure, and may have important implications for understanding the functional energetics and specificity of the interactions of DNA and RNA with regulatory proteins, pharmaceutical agents, and other ligands.

Kinetics of DNA hydration

The hydration of the d(CGCGAATTCGCG) B-DNA duplex in solution was studied by nuclear magnetic relaxation dispersion (NMRD) of the water nuclei 1H, 2H, and 17O, and by nuclear Overhauser effects (NOEs) in high-resolution two-dimensional 1H NMR spectra. By comparing results from the free duplex with those from its complex with netropsin, water molecules in the "spine of hydration" in the AATT region of the minor groove could be distinguished from hydration water elsewhere in the duplex. The 2H and 17O relaxation dispersions yield a model-independent residence time of 0.9(+/-0.1) ns at 4 degrees C for five highly ordered water molecules in the spine. When corrected for frequency offset effects, the NOE data yield the same residence time as the NMRD data, giving credence to both methods. At 27 degrees C, the residence time is estimated to 0.2 ns, a factor of 40 shorter than the tumbling time of the duplex. The NMRD data show that all water molecules associated with the duplex, except the five molecules in the spine, have residence times significantly shorter than 1 ns at 4 degrees C. There is thus no long-lived hydration structure associated with the phosphate backbone. In contrast to 2H and 17O, the 1H relaxation dispersion is dominated by labile DNA protons and therefore provides little information about DNA hydration.

Site-specific dynamics in DNA: experiments

This chapter reviews the dynamics information obtained from experimental magnetic resonance studies of site-specifically labeled duplex DNA. A previous review (43) discusses the dynamics of duplex DNA; it develops a theory that shows how magnetic resonance experiments are used to detect those dynamics. The methods for obtaining information about dynamics as well as a summary of what is now known about the site-specific dynamics of DNA are presented. This review contains two methods sections which present results using electron paramagnetic resonance and nuclear magnetic resonance active probes.

Modifying the helical structure of DNA by design: recruitment of an architecture-specific protein to an enforced DNA bend

BACKGROUND

Proteins can force DNA to adopt distorted helical structures that are rarely if ever observed in naked DNA. The ability to synthesize DNA that contains defined helical aberrations would offer a new avenue for exploring the structural and energetic plasticity of DNA. Here we report a strategy for the enforcement of non-canonical helical structures through disulfide cross-linking; this approach is exemplified by the design and synthesis of an oligonucleotide containing a pronounced bend.

RESULTS

A localized bend was site-specifically introduced into DNA by the formation of a disulfide cross-link between the 5' adenines of a 5'-AATT-3' region in complementary strands of DNA. The DNA bend was characterized by high-resolution NMR structure determination of a cross-linked dodecamer and electrophoretic mobility assays on phased multimers, which together indicate that the cross-linked tetranucleotide induces a helical bend of approximately 30 degrees and a modest degree of unwinding. The enforced bend was found to stimulate dramatically the binding of an architecture-specific protein, HMG-D, to the DNA. DNase I foot-printing analysis revealed that the protein is recruited to the section of DNA that is bent.

CONCLUSIONS

The present study reports a novel approach for the investigation of non-canonical DNA structures and their recognition by architecture-specific proteins. The mode of DNA bending induced by disulfide cross-linking resembles that observed in structures of protein-DNA complexes. The results reveal common elements in the DNA-binding mode employed by sequence-specific and architecture-specific HMG proteins.

Molecular dynamics simulations of the effects of ring-saturated thymine lesions on DNA structure

The effect of thymine lesions produced by radiation or oxidative damage on DNA structure was studied by molecular dynamics simulations of native and damaged DNA. Thymine in position 7 of native dodecamer d(CGCGAATTCGCG)2 was replaced by one of the four thymine lesions 5-hydroxy-5,6-dihydrothymine, 6-hydroxy-5,6-dihydrothymine (thymine photohydrate), 5,6-dihydroxy-5,6-dihydro-thymine (thymine glycol), and 5,6-dihydrothymine. Simulations were performed with Assisted Model Building with Energy Refinement force field. Solvent was represented by a rectangular box of water with periodic boundary conditions applied. A constant temperature and constant volume protocol was used. The observed level of distortions of DNA structure depends on the specific nature of the lesion. The 5,6-dihydrothymine does not cause distinguishable perturbations to DNA. Other lesions produce a dramatic increase in the rise parameter between the lesion and the 5' adjacent adenine. These changes are accompanied by weakening of Watson-Crick hydrogen bonds in the A6-T19 base pair on the 5' side of the lesion. The lesioned bases also show negative values of inclination relative to the helical axis. No changes in the pattern of backbone torsional angles are observed with any of the lesions incorporated into DNA. The structural distortions in DNA correlate well with known biological effects of 5,6-dihydrothymine and thymine glycol on such processes as polymerase action or recognition by repair enzymes.

The effects of N7-methylguanine on duplex DNA structure

BACKGROUND

Non-enzymatic methylation of DNA by endogenous and exogenous agents produces a variety of adducts, of which the predominant one is N7-methyl-2'-deoxyguanosine (m7dG). Although it is known that living organisms counter the deleterious effects of m7dG by producing adduct-specific DNA repair proteins, the molecular basis for specific recognition and catalysis by these proteins is poorly understood. In addition to its role as an endogenous DNA adduct, m7dG is also widely used as an in vitro probe of protein-DNA interactions. We set out to examine whether incorporation of m7dG into DNA affects duplex DNA structure.

RESULTS

We carried out a large-scale synthesis of a dodecamer containing the m7dG adduct at a single, defined position. Because the instability of m7dG precludes its incorporation into oligonucleotides by standard solid-phase methods, a novel strategy employing chemical and enzymatic synthesis was used. Characterization of the m7dG-containing dodecamer by NMR reveals no structural distortion; indeed, m7dG appears to encourage a modest shift toward a more characteristically B-form duplex.

CONCLUSIONS

These results argue strongly against induced DNA distortion as a mechanism for specific recognition of m7dG by adduct-specific repair proteins. The broad substrate specificity of these repair proteins disfavors a model involving direct recognition of aberrantly placed methyl groups; hence, it may be that m7dG is recognized indirectly, perhaps by its effects on the dynamics of DNA. On the other hand, the evidence presented here suggests that m7dG interferes directly with sequence-specific recognition by DNA-binding proteins by steric blockage or by masking of required contact functionalities. The synthetic methodology used here should be generally applicable to high-resolution structural studies of oligonucleotides bearing adducts that are unstable to the conditions of solid-phase DNA synthesis.

Correlations of crystal structures of DNA oligonucleotides with enantioselective recognition by Rh(phen)2phi3+: probes of DNA propeller twisting in solution

Rh(phen)2phi3+ (phen = 1,10-phenanthroline; phi = 9,10-phenanthrenequinone diimine) avidly binds to DNA via intercalation from the major groove and upon photoactivation produces strand scission with single base 5' asymmetry. Enantiomers of Rh(phen)2phi3+, which lack hydrogen-bonding substituents in ancillary positions, distinguish a DNA site through shape-selection; site recognition depends upon the local variations in structure at a binding site. Here, we examine the application of delta-Rh(phen)2phi3+ as a sequence-dependent structural probe and, in particular, as a probe of DNA propeller twisting in solution, by comparing directly cleavage results using delta- and lambda-Rh(phen)2phi3+ on crystallographically characterized oligonucleotides with several sequence-dependent crystallographic parameters. The three oligonucleotides examined in this study are the Dickerson-Drew dodecamer, 5'-CGCGAATTCGCG-3', the NarI dodecamer, 5'-ACCGGCGCCACA-3', and the CG decamer, 5'-CCAACGTTGG-3', all of which have been crystallized in the B-form. Enantioselective cleavage and reaction favored by the delta-isomer is found to be governed locally by the opening of the site in the major groove. A correlation is demonstrated between cleavage by delta-Rh(phen)2phi3+ and the opening in the major groove that results from the change in propeller twist (differential propeller twist) at a base step. When the major groove is closed as a result of a change in propeller twist, there is little cleavage evident by either enantiomer; at sites which are indicated crystallographically to be open in the major groove, a direct correlation is observed between enantioselective cleavage and the degree of opening. A trend of higher enantioselectivity at sites possessing higher twist angles is also observed.(ABSTRACT TRUNCATED AT 250 WORDS)

The amplitude of local angular motion of purines in DNA in solution

Nuclear magnetic resonance and optical experiments are combined to determine the rms amplitude of local angular motion of purines in DNA in solution. A 12 base-pair duplex DNA with the sequence d(CGCGAATTCGCG)2 is deuterated at the H8 positions of adenine and guanine by exchange with solvent at 55 degrees C. The deuterium nmr spectrum of this DNA is measured at 30 mg/mL at 30 degrees C in an 11.76 Tesla magnet (76.75 MHz). The time-resolved fluorescence polarization anisotropies (FPA) of this same sample and also a greatly diluted sample (0.215 mg/mL) were measured after addition of ethidium. FPA measurements of the dilute sample yield the hydrodynamic radius, RH = 9.94 +/- 0.2 A, while those at the nmr concentration are employed to characterize the collective motions in terms of either an enhanced viscosity or dimer formation. The rms amplitude of local angular motion was determined by analyzing the 2H-nmr spectrum, in particular the line width, using recently developed theory for the transverse relaxation rate (RQ2) together with essential information about the collective motions from these and other optical studies. When the principal-axis frame of the electric field gradient tensor is assumed to undergo overdamped libration around each of its three body-fixed axes in an isotropic deflection potential, then the rms amplitude of local angular motion around any single axis is found to lie in the range 10 degrees-11 degrees, provided the high DNA concentration acts to enhance the viscosity, and is about 9 degrees-11 degrees, if it acts to produce end-to-end dimers. The proton nmr relaxation data of Eimer et al. are reanalyzed and shown to yield an rms amplitude of angular motion of the cytosine H5-H6 internuclear vector of 9 degrees-10 degrees, depending upon its orientation with respect to the helix axis. In all of these analyses, full account is taken of the collective twisting and bending deformations, which have a small but significant effect on the results. It is shown that the rms amplitudes of local angular motion do not depend strongly on the model (potential), provided that isotropic rotation around the same number of axes is allowed and that one compares rms angles of the same dimensionality. The rms amplitudes of local angular motion in solution are comparable to those observed for the same sequence at low levels of hydration in the solid state.

Determination of distances involving exchangeable protons from NOESY spectra of two DNA dodecamers

Two-dimensional proton nuclear Overhauser effect (NOESY) spectra were obtained as a function of mixing time for two DNA dodecamers, 5'-CGCGAATTCGCG-3' and 5'-CGCAAATTTGCG-3', in aqueous solutions. The time evolution of cross-peak volumes was quantitatively analyzed based on the approximate solution of the relaxation/exchange equation over a range of mixing times from 30 to 150 ms. Inter-proton distances, involving exchangeable protons in key locations of the structures, were calculated and compared to the distances predicted by the crystal structures. These NMR-derived distances enhance the number of constraints, and their accuracy, for determination of solution structures of the two DNA dodecamers.

Molecular dynamics simulations of B-DNA: an analysis of the role of initial molecular configuration, randomly assigned velocity distribution, long integration times, and nonconstrained termini

Molecular dynamics simulations of three DNA sequences using the AMBER 3.0 force field were performed with implicit inclusion of water through a distance-dependent dielectric constant and solvated counterions. Simulations of the self-complementary DNA dodecamer d(CGCGAATTCGCG) were started from a regular B-DNA structure and the x-ray single crystal B-DNA structure. Although mean convergence during the 89-ps calculation was confirmed, localized differences in backbone torsionals and base-pair helicoidals were observed. A nanosecond simulation of the nonself-complementary 14 base-pair DNA d(GGCGGAATTGGCGG) indicates that most structural parameters stabilize within the first 100-200 ps, while isolated features show low-frequency oscillations throughout the calculation. The lack of harmonic constraints on the ends of the molecules was shown not to perturb the structural dynamics of the internal oligonucleotide beyond the external 2 base pairs. Comparison of three simulations of the nonself-complementary 14 base-pair DNA d(GGCGAAATTCGCGG), identical in all respects other than the assignment of initial Maxwellian atomic velocity distributions, revealed the inherent systematic variability. The three calculations result in nearly superimposable global structures, with localized variations in torsionals and helicoidals. Our results provide a basis for performing a comparative analysis of the effect of DNA sequence on localized structure.

Investigation of solution structure of d(GAATTTAAATTC)2 by 1H NMR, molecular dynamics, mechanics, refinement by back-calculation of the NOESY spectrum and analysis of this structure using X-ray data

1H NMR spectroscopy, restrained molecular mechanics and dynamics and refinements after back-calculation of the NOESY spectrum have been performed to study the structure of the d(GAATTTAAATTC)2 duplex and to determine whether it is bent or not. It is found that the duplex adopts a B-type conformation; all sugar conformations belong to the C2' endo region and purines have a larger pseudorotation angle as compared to pyrimidines. The cross-strand AH2(n)-AH1' (m + 1) distance (where (n) and (m) are complementary residues), crucial for an anomalous A/T tract structure, is large on the TA step and gradually decreases at the 3' and 5' ends of the TTTAAA tract and follows the rules proposed previously [Chuprina, V.P., Lipanov, A.A., Fedoroff, O.Yu, Kim, S.G., Kintanar, A., and Reid, B.R., Proc. Natl. Acad. Sci. U.S.A. 88, 9087 (1991)]. The changes in this distance correlate with those in the T1 value for AH2 protons which we measured for several oligonucleotide sequences. A total number of about 250 interproton distance constraints were determined from NOESY spectra and were used for structure determination by molecular mechanics, dynamics and refinement by back-calculations. It is shown that these data are not enough to determine whether the duplex is bent or not. The whole family of B-type conformations including bent and straight structures fit well with the available NMR data. In principle, additional non-NMR data could be used in order to reduce the number of the allowable structures. The refinement of the structures with additional different non-NMR constraints (used as a driving force) on P-P or H1'-H1' minor groove width distances in the TA region shows a very good correlation between these distances and the angle of bending of the dodecamer. The more the minor groove width increases in the TA region the more the duplex is bent at the major groove of this region. On the other hand, there is also a very good correlation between P-P, H1'-H1' and AH2-H1' cross-strand distances as follows from analysis of X-ray B-type structures. These two correlations, together with the increased AH2-H1' cross-strand NMR distance in the TA region of the dodecamer indicate that the duplex should be characterized by a wider minor groove in the TA region and be bent in the major groove in this region.

Computational challenges for macromolecular structure determination by X-ray crystallography and solution NMR-spectroscopy

Macromolecular structure determination by X-ray crystallography and solution NMR spectroscopy has experienced unprecedented growth during the past decade.

Solution structure of the TnAn DNA duplex GCCGTTAACGCG containing the HpaI restriction site

The solution structure of the self-complementary DNA duplex [d(GCCGTTAACGGC)]2, which contains the HpaI restriction site GTTAAC, has been elucidated by two-dimensional NMR, distance geometry (DG), and NOE back-calculation methods. Initial distance constraints were determined by polynomial fitting the two-spin initial NOE rates; backbone constraints from NOE and J-coupling observations (Kim et al., 1992) were included. RMSDs between initial-distance-refined structures derived from random-embedded DG, A-DNA, and B-DNA starting structures were all in the range 0.5-1.0 A, indicating good convergence properties of the algorithm, regardless of the starting structure. A semiautomatic back-calculation refinement procedure was developed and used to generate more refined structures for which the BKCALC-simulated NOE volumes matched the experimental data. The six final structures refined from various starting structures exhibit very good agreement with the experimental data (R values = 0.18) and converge well to within 0.8-A RMSD differences for the central 8 base pairs. The torsion and pseudorotation phase angles were found to be well determined by the data, and the local helical parameters for each base step converged quite well. The final structures show that the central T6-A7 step is somewhat underwound (twist angle ca. 29 degrees), with a large negative cup and a normal (wide) minor groove width, while the T5-T6 and A7-A8 steps have a partially narrowed minor groove.

Solution structure of [d(GCGTATACGC)]2

The solution structure of the alternating pyrimidine-purine DNA duplex [d(GCGTATACGC)]2 has been determined using two-dimensional nuclear magnetic resonance techniques and distance geometry methods. Backbone distance constraints derived from experimental nuclear Overhauser enhancement and J-coupling torsion angle constraints were required to adequately define the conformation of the inter-residue backbone linkages and to avoid underwinding of the duplex. The distance geometry structures were further refined by back-calculation of the two-dimensional nuclear Overhauser enhancement spectra to correct spin-diffusion distance errors. Fifteen final structures for [d(GCGTATACGC)]2 were generated from the refined experimental distance bounds. These structures all exhibit fully wound B-form geometry with small penalty values (< 1.5 A) against the distance bounds and small pair-wise root-mean-square deviation values (typically 0.6 A to 1.5 A). The final structures exhibit positive base-pair inclination with respect to the helix axis, a marked alternation in rise and twist, and are shorter and wider than classical fiber B-form DNA. The purines were found to adopt a sugar pucker close to the C-2'-endo conformation while pyrimidine sugars exhibited significantly lower pseudorotation phase angles in the C-1'-exo to C-2'-endo range. The minor groove cross-strand steric clashes at pyrimidine-purine steps that would exist in pure B-DNA are attenuated by an increased rise at these steps (and an increased roll angle at TpA steps). Concomitantly the backbone torsion angles of the pyrimidine moieties have larger gamma values, larger epsilon values, and smaller zeta values than the purines. The structures generated by distance geometry methods were also compared with those obtained from restrained molecular dynamics with empirical force-field potentials. The results indicate that the nuclear magnetic resonance/distance geometry approach alone is capable of elucidating most of the salient structural features of double-stranded helical nucleic acids in solution without resorting to empirical energy potentials and without using any structural assumptions from crystallographic data.

Solution structure of [d(ATGAGCGAATA)]2. Adjacent G:A mismatches stabilized by cross-strand base-stacking and BII phosphate groups

The solution structure of a rather unusual B-form duplex [d(ATGAGCGAATA)]2 has been determined using two-dimensional nuclear magnetic resonance (2D-NMR) and distance geometry methods. This sequence forms a stable ten base-pair B-form duplex with 3' overhangs and two pairs of adjacent G:A mismatches paired via a sheared hydrogen-bonding scheme. All non-exchangeable protons, including the stereo-specific H-5'S/H-5'R of the 3G and 7G residues, were assigned by 2D-NMR. The phosphorus spectrum was assigned using heteronuclear correlation with H-3' and H-4' reasonances. The complete assignments reveal several unusual nuclear Overhauser enhancements (NOEs) and unusual chemical shifts for the neighboring G:A mismatch pairs and their adjacent nucleotides. Inter-proton distances were derived from time-dependent NOEs and used to generate initial structures, which were further refined by iterative back-calculation of the two-dimensional nuclear Overhauser enhancement spectra; 22 final structures were calculated from the refined distance bounds. All these final structures exhibit fully wound helical structures with small penalty values against the refined distance bounds and small pair-wise root-mean-square deviation values (typically 0.5 A to 0.9 A). The two helical strands exchange base stacking at both of the two G:A mismatch sites, resulting in base stacking down each side rather than down each strand of the twisted duplex. Very large twist angles (77 degrees) were found at the G:A mismatch steps. All the final structures were found to have BII phosphate conformations at the adjacent G:A mismatch sites, consistent with observed downfield 31P chemical shifts and Monte-Carlo conformational search results. Our results support the hypothesis that 31P chemical shifts are related to backbone torsion angles. These BII phosphate conformations in the adjacent G:A mismatch step suggest that hydrogen bonding of the G:A pair G-NH2 to a nearby phosphate oxygen atom is unlikely. The unusual structure of the duplex may be stabilized by strong interstrand base stacking as well as intrastrand stacking, as indicated by excellent base overlap within the mismatch stacks.