Sequence-dependent dynamics of duplex DNA: the applicability of a dinucleotide model

The short-time (submicrosecond) bending dynamics of duplex DNA were measured to determine the effect of sequence on dynamics. All measurements were obtained from a single site on duplex DNA, using a single, site-specific modified base containing a rigidly tethered, electron paramagnetic resonance active spin probe. The observed dynamics are interpreted in terms of single-step sequence-dependent bending force constants, determined from the mean squared amplitude of bending relative to the end-to-end vector using the modified weakly bending rod model. The bending dynamics at a single site are a function of the sequence of the nucleotides constituting the duplex DNA. We developed and examined several dinucleotide-based models for flexibility. The models indicate that the dominant feature of the dynamics is best explained in terms of purine- and pyrimidine-type steps, although distinction is made among all 10 unique steps: It was found that purine-purine steps (which are the same as pyrimidine-pyrimidine steps) were near average in flexibility, but the pyrimidine-purine steps (5' to 3') were nearly twice as flexible, whereas purine-pyrimidine steps were more than half as flexible as average DNA. Therefore, the range of stepwise flexibility is approximately fourfold and is characterized by both the type of base pair step (pyrimidine/purine combination) and the identity of the bases within the pair (G, A, T, or C). All of the four models considered here underscore the complexity of the dependence of dynamics on DNA sequence with certain sequences not satisfactorily explainable in terms of any dinucleotide model. These findings provide a quantitative basis for interpreting the dynamics and kinetics of DNA-sequence-dependent biological processes, including protein recognition and chromatin packaging.     

Modeling DNA Thermodynamics under Torsional Stress

Negatively twisted DNA is essential to many biological functions. Due to torsional stress, duplex DNA can have local, sequence-dependent structural defects. In this work, a thermodynamic model of DNA was built to qualitatively predict the local sequence-dependent mechanical instabilities under torsional stress. The results were compared to both simulation of a coarse-grained model and experiment results. By using the Kirkwood superposition approximation, we built an analytical model to represent the free energy difference ΔW of a hydrogen-bonded basepair between the B-form helical state and the basepair opened (or locally melted) state, within a given sequence under torsional stress. We showed that ΔW can be well approximated by two-body interactions with its nearest-sequence-neighbor basepairs plus a free energy correction due to long-range correlations. This model is capable of rapidly predicting the position and thermodynamics of local defects in a given sequence. The result qualitatively matches with an in vitro experiment for a long DNA sequence (>4000 basepairs). The 12 parameters used in this model can be further quantitatively refined when more experimental data are available.     

Modeling DNA Thermodynamics under Torsional Stress

Negatively twisted DNA is essential to many biological functions. Due to torsional stress, duplex DNA can have local, sequence-dependent structural defects. In this work, a thermodynamic model of DNA was built to qualitatively predict the local sequence-dependent mechanical instabilities under torsional stress. The results were compared to both simulation of a coarse-grained model and experiment results. By using the Kirkwood superposition approximation, we built an analytical model to represent the free energy difference ΔW of a hydrogen-bonded basepair between the B-form helical state and the basepair opened (or locally melted) state, within a given sequence under torsional stress. We showed that ΔW can be well approximated by two-body interactions with its nearest-sequence-neighbor basepairs plus a free energy correction due to long-range correlations. This model is capable of rapidly predicting the position and thermodynamics of local defects in a given sequence. The result qualitatively matches with an in vitro experiment for a long DNA sequence (>4000 basepairs). The 12 parameters used in this model can be further quantitatively refined when more experimental data are available.     

Flexibility of duplex DNA on the submicrosecond timescale

Using a site-specific, Electron Paramagnetic Resonance (EPR)-active spin probe that is more rigidly locked to the DNA than any previously reported, the internal dynamics of duplex DNAs in solution were studied. EPR spectra of linear duplex DNAs containing 14-100 base pairs were acquired and simulated by the stochastic Liouville equation for anisotropic rotational diffusion using the diffusion tensor for a right circular cylinder. Internal motions have previously been assumed to be on a rapid enough time scale that they caused an averaging of the spin interactions. This assumption, however, was found to be inconsistent with the experimental data. The weakly bending rod model is modified to take into account the finite relaxation times of the internal modes and applied to analyze the EPR spectra. With this modification, the dependence of the oscillation amplitude of the probe on position along the DNA was in good agreement with the predictions of the weakly bending rod theory. From the length and position dependence of the internal flexibility of the DNA, a submicrosecond dynamic bending persistence length of around 1500 to 1700 A was found. Schellman and Harvey (Biophys. Chem. 55:95-114, 1995) have estimated that, out of the total persistence length of duplex DNA, believed to be about 500 A, approximately 1500 A is accounted for by static bends and 750 A by fluctuating bends. A measured dynamic persistence length of around 1500 A leads to the suggestion that there are additional conformations of the DNA that relax on a longer time scale than that accessible by linear CW-EPR. These measurements are the first direct determination of the dynamic flexibility of duplex DNA in 0.1 M salt.     

Experimental mapping of DNA duplex shape enabled by global lineshape analyses of a nucleotide-independent nitroxide probe

Sequence-dependent variation in structure and dynamics of a DNA duplex, collectively referred to as ‘DNA shape’, critically impacts interactions between DNA and proteins. Here, a method based on the technique of site-directed spin labeling was developed to experimentally map shapes of two DNA duplexes that contain response elements of the p53 tumor suppressor. An R5a nitroxide spin label, which was covalently attached at a specific phosphate group, was scanned consecutively through the DNA duplex. X-band continuous-wave electron paramagnetic resonance spectroscopy was used to monitor rotational motions of R5a, which report on DNA structure and dynamics at the labeling site. An approach based on Pearson''s coefficient analysis was developed to collectively examine the degree of similarity among the ensemble of R5a spectra. The resulting Pearson''s coefficients were used to generate maps representing variation of R5a mobility along the DNA duplex. The R5a mobility maps were found to correlate with maps of certain DNA helical parameters, and were capable of revealing similarity and deviation in the shape of the two closely related DNA duplexes. Collectively, the R5a probe and the Pearson''s coefficient-based lineshape analysis scheme yielded a generalizable method for examining sequence-dependent DNA shapes.     

Stretched DNA Investigated Using Molecular-Dynamics and Quantum-Mechanical Calculations

We combined atomistic molecular-dynamics simulations with quantum-mechanical calculations to investigate the sequence dependence of the stretching behavior of duplex DNA. Our combined quantum-mechanical/molecular-mechanical approach demonstrates that molecular-mechanical force fields are able to describe both the backbone and base-base interactions within the highly distorted nucleic acid structures produced by stretching the DNA from the 5′ ends, which include conformations containing disassociated basepairs, just as well as these force fields describe relaxed DNA conformations. The molecular-dynamics simulations indicate that the force-induced melting pathway is sequence-dependent and is influenced by the availability of noncanonical hydrogen-bond interactions that can assist the disassociation of the DNA basepairs. The biological implications of these results are discussed.     

Modeling translocation dynamics of strand displacement DNA synthesis by DNA polymerase I

A model is presented for the translocation dynamics of the strand displacement DNA synthesis by DNA polymerases such as polymerase I family. (i) The model gives an explanation to the experimental results which showed that the rate of strand displacement DNA synthesis is nearly consistent with that of single stranded primer extension synthesis, although the two are expected to have substantial differences in their energetics. (ii) During strand displacement DNA synthesis, the pausing at the specific sequence is considered to be due to an affinity of the fingers subdomain for the specific sequence of dsDNA downstream of the single strand. The theoretical results on the sequence-dependent pausing dynamics such as the mean pausing lifetimes and the distribution of the pausing lifetime are consistent with the experimental data. Moreover, predicted results are presented for the binding affinity of the fingers subdomain for the specific sequence of dsDNA and the dependence of the mean sequence-dependent pausing lifetime on the external force acting on the polymerase.     

Sequence-dependent twist-stretch coupling in DNA

Recent single-molecule micromanipulation experiments on DNA subject to small distortion revealed positive coupling between DNA stretching and twisting—for instance, DNA elongates when overtwisted. Here we propose a method to calculate the twist-stretch coupling constant specific to a DNA fragment of a given sequence. The method employs a sequence-dependent dinucleotide force field and is based on constrained minimization of the fragment's deformation energy. Using a force field inferred from atomistic molecular dynamics simulations, we obtain the twist-stretch coupling for random sequence to be 0.30 nm/turn, close to experimental values. An exhaustive calculation for all oligomers of nine basepairs yields values between 0.14 and 0.45 nm/turn, positively correlated with the contents of pyrimidine-purine steps in the sequence. Our method is simple to use and allows one to explore the hypothesis that some sequences may be optimized for twist-stretch coupling.     

Effect of a neutralized phosphate backbone on the minor groove of B-DNA: molecular dynamics simulation studies

Alternative models have been presented to provide explanations for the sequence-dependent variation of the DNA minor groove width. In a structural model groove narrowing in A-tracts results from direct, short-range interactions among DNA bases. In an electrostatic model, the narrow minor groove of A-tracts is proposed to respond to sequence-dependent localization of water and cations. Molecular dynamics simulations on partially methylphosphonate substituted helical chains of d(TATAGGCCTATA) and d(CGCGAATTCGCG) duplexes have been carried out to help evaluate the effects of neutralizing DNA phosphate groups on the minor groove width. The results show that the time-average minor groove width of the GGCC duplex becomes significantly more narrow on neutralizing the phosphate backbone with methylphosphonates. The minor groove of the AATT sequence is normally narrow and the methylphosphonate substitutions have a smaller but measurable affect on this sequence. These results and models provide a system that can be tested by experiment and they support the hypothesis that the electrostatic environment around the minor groove affects the groove width in a sequence-dependent dynamic and time-average manner.     

PELDOR analysis of enzyme-induced structural changes in damaged DNA duplexes

PELDOR (pulsed electron-electron double resonance) spectroscopy was applied to determine spin-spin distances in spin-labeled DNA duplexes (13-mer and 17-mer) containing the damaged sites 8-oxoguanine or uncleavable abasic site analogue tetrahydrofuran. The lesions were located in one strand of the DNA, and two nitroxyl spin labels were attached at the 5'- and 3'-ends of the complementary strand. PELDOR data allow us to obtain distances between the two spin labels in DNAs, which turned out to be around 5 nm for the 13-mer DNA and around 6 nm for 17-mer DNA. Results of PELDOR measurements were supported by molecular dynamics calculations. Study of the interaction of DNA fragments with DNA repair enzyme 8-oxoguanine-DNA glycosylase from E. coli (Fpg protein) showed that this interaction leads to a noticeable decrease of the distance between spin labels, which indicates the enzyme-induced bending of the DNA duplex. This bending may be important for the mechanisms of recognition of damaged sites by DNA repair enzymes.     

Sulfur signaling: is the agent sulfide or sulfane?

Triplex-forming oligonucleotides (TFOs) are sequence-dependent DNA binders that may be useful for DNA targeting and detection. A sensitive and convenient method to monitor triplex formation by a TFO and its target DNA duplex is required for the application of TFO probes. Here we describe a novel design by which triplex formation can be monitored homogeneously without prelabeling the target duplex. The design uses a TFO probe tagged with a fluorophore that undergoes fluorescence resonance energy transfer with fluorescent dyes that intercalate into the target duplex. Through color compensation analysis, the specific emission of the TFO probe reveals the status of the triple helices. We used this method to show that triple helix formation with TFOs is magnesium dependent. We also demonstrated that the TFO probe can be used for detection of sequence variation in melting analysis and for DNA quantitation in real-time polymerase chain reaction.     

Anthracycline antibiotic arugomycin binds in both grooves of the DNA helix simultaneously: an NMR and molecular modelling study.

Perturbations to the 1H and 31P chemical shifts of DNA resonances together with twenty-four intermolecular nuclear Overhauser effects show that the anthracycline antibiotic arugomycin intercalates between the basepairs of the hexamer duplex d(5'-GCATGC)2 at the 5'-CpA and 5'-TpG binding sites. In the complex two drug molecules are bound per duplex with full retention of the dyad symmetry. Arugomycin adopts a threaded binding orientation with chains of sugars positioned in both the major and minor groove of the helix simultaneously. The complex is stabilized by hydrogen bonding, electrostatic and van der Waals interactions principally in the major groove and involving substituents on the rigidly oriented bicycloamino-glucose sugar of the antibiotic. A specific hydrogen bond is identified between the C2'-hydroxyl and the guanine N7 at the intercalation site. Together, interactions in the major groove appear to account for the intercalation specificity of arugomycin that requires both a guanine and thymine at the intercalation site. We are unable to identify any sequence specific interactions between the minor groove and the arugarose sugar (S1) which binds only weakly, through van der Walls contacts, over the d(GCA).d(TGC) trinucleotide sequence. The data indicate that the sugar chains of arugomycin are flexible and play little part in the interaction of the antibiotic with DNA. The intensity of sequential internucleotide NOEs identifies the intercalation site as being assymmetric. A family of conformers computed using restrained energy minimisation and molecular dynamics indicate that basepair buckling is a feature of the anthracycline intercalation site that may serve to maximise intermolecular van der Waals interactions by wrapping the basepairs around the antibiotic chromophore.     

Steady-state and time-resolved fluorescence studies indicate an unusual conformation of 2-aminopurine within ATAT and TATA duplex DNA sequences

2-Aminopurine (2-AP), a fluorescent analog of adenine, has been widely used as a probe for local DNA conformation, since excitation and emission characteristics and fluoresence lifetimes of 2-AP vary in a sequence-dependent manner within DNA. Using steady-state and time-resolved fluorescence techniques, we report that 2-AP appears to be unusually stacked in the internal positions of ATAT and TATA in duplex DNA. The excitation wavelength maxima for 2-AP within these contexts were red shifted, indicating reduced solvent exposure for the fluorophore. Furthermore, in these contexts, 2-AP fluorescence was resistant to acrylamide-dependent collisional quenching, suggesting that the fluorophore is protected by its stacked position within the duplex. This conclusion was further reinforced by the presence of a secondary peak at 275 nm in the fluorescence excitation spectra that is indicative of efficient excitation energy transfer from nearby non-fluorescent DNA bases. Fluorescence anisotropy decay and internal angular ‘wobbling’ motion measurements of 2-AP within these alternating AT contexts were also consistent with the fluorophore being highly constrained and immobile within the base stack. When these fluorescence characteristics are compared with those of 2-AP within other duplex DNA sequence contexts, they are unique.     

A Statistical Thermodynamic Model for Investigating the Stability of DNA Sequences from Oligonucleotides to Genomes

We describe the development and testing of a simple statistical mechanics methodology for duplex DNA applicable to sequences of any composition and extensible to genomes. The microstates of a DNA sequence are modeled in terms of blocks of basepairs that are assumed to be fully closed (paired) or open. This approach generates an ensemble of bubblelike microstates that are used to calculate the corresponding partition function. The energies of the microstates are calculated as additive contributions from hydrogen bonding, basepair stacking, and solvation terms parameterized from a comprehensive series of molecular dynamics simulations including solvent and ions. Thermodynamic properties and nucleotide stability constants for DNA sequences follow directly from the partition function. The methodology was tested by comparing computed free energies per basepair with the experimental melting temperatures of 60 oligonucleotides, yielding a correlation coefficient of −0.96. The thermodynamic stability of genic/nongenic regions was tested in terms of nucleotide stability constants versus sequence for the Escherichia coli K-12 genome. It showed clear differentiation of the genes from promoters and captures genic regions with a sensitivity of 0.94. The statistical thermodynamic model presented here provides a seemingly new handle on the challenging problem of interpreting genomic sequences.     

Sequence-dependence of the energetics of opening of at basepairs in DNA

Proton exchange and nuclear magnetic resonance spectroscopy are being used to characterize the energetics of opening of AT/TA basepairs in the DNA dodecamer 5'-d(GCTATAAAAGGG)-3'/5'-d(CCCTTTTATAGC)-3'. The dodecamer contains the TATA box of the adenovirus major late promoter. The equilibrium constants for opening of each basepair are measured from the dependence of the exchange rates of imino protons on ammonia concentration. The enthalpy, entropy, and free energy changes in the opening reaction of each basepair are determined from the temperature dependence of the exchange rates. The results reveal that the opening enthalpy changes encompass a wide range of values, namely, from 17 to 29 kcal/mol. The largest values are observed for the AT basepairs in 7th and 8th positions. These values, and the exchange rates of the corresponding imino protons, suggest that these two basepairs open in a single concerted reaction. The enthalpy changes for opening of the central six basepairs are correlated to the opening entropy changes. This enthalpy-entropy compensation minimizes the variations in the opening free energies among these central basepairs. Deviations from the enthalpy-entropy compensation pattern are observed for basepairs located close to the ends of the duplex structure, suggesting a different mode of opening for these basepairs.     

A Statistical Thermodynamic Model for Investigating the Stability of DNA Sequences from Oligonucleotides to Genomes

We describe the development and testing of a simple statistical mechanics methodology for duplex DNA applicable to sequences of any composition and extensible to genomes. The microstates of a DNA sequence are modeled in terms of blocks of basepairs that are assumed to be fully closed (paired) or open. This approach generates an ensemble of bubblelike microstates that are used to calculate the corresponding partition function. The energies of the microstates are calculated as additive contributions from hydrogen bonding, basepair stacking, and solvation terms parameterized from a comprehensive series of molecular dynamics simulations including solvent and ions. Thermodynamic properties and nucleotide stability constants for DNA sequences follow directly from the partition function. The methodology was tested by comparing computed free energies per basepair with the experimental melting temperatures of 60 oligonucleotides, yielding a correlation coefficient of −0.96. The thermodynamic stability of genic/nongenic regions was tested in terms of nucleotide stability constants versus sequence for the Escherichia coli K-12 genome. It showed clear differentiation of the genes from promoters and captures genic regions with a sensitivity of 0.94. The statistical thermodynamic model presented here provides a seemingly new handle on the challenging problem of interpreting genomic sequences.     

Predicting sequence-dependent melting stability of short duplex DNA oligomers

Many important applications of DNA sequence-dependent hybridization reactions have recently emerged. This has sparked a renewed interest in analytical calculations of sequence-dependent melting stability of duplex DNA. In particular, for many applications it is often desirable to accurately predict the transition temperature, or t_m, of short duplex DNA oligomers (∼ 20 base pairs or less) from their sequence and concentration. The thermodynamic analytical method underlying these predictive calculations is based on the nearest-neighbor model. At least 11 sets of nearest-neighbor sequence-dependent thermodynamic parameters for DNA have been published. These sets are compared. Use of the nearest-neighbor sets in predicting t_m from the DNA sequence is demonstrated, and the ability of the nearest-neighbor parameters to provide accurate predictions of experimental t_m's of short duplex DNA oligomers is assessed. © 1998 John Wiley & Sons, Inc. Biopoly 44: 217–239, 1997     

Statistical mechanics of sequence-dependent circular DNA and its application for DNA cyclization

DNA cyclization is potentially the most powerful approach for systematic quantitation of sequence-dependent DNA bending and flexibility. We extend the statistical mechanics of the homogeneous DNA circle to a model that considers discrete basepairs, thus allowing for inhomogeneity, and apply the model to analysis of DNA cyclization. The theory starts from an iterative search for the minimum energy configuration of circular DNA. Thermodynamic quantities such as the J factor, which is essentially the ratio of the partition functions of circular and linear forms, are evaluated by integrating the thermal fluctuations around the configuration under harmonic approximation. Accurate analytic expressions are obtained for equilibrium configurations of homogeneous circular DNA with and without bending anisotropy. J factors for both homogeneous and inhomogeneous DNA are evaluated. Effects of curvature, helical repeat, and bending and torsional flexibility in DNA cyclization are analyzed in detail, revealing that DNA cyclization can detect as little as one degree of curvature and a few percent change in flexibility. J factors calculated by our new approach are well consistent with Monte Carlo simulations, whereas the new theory has much greater efficiency in computations. Simulation of experimental results has been demonstrated.