NMR evidence for base dynamics at all TpA steps in DNA

NMR evidence is presented indicating that the exceptional conformational dynamics found at TpA steps in DNA is general to all immediate sequence contexts. One easily tractable NMR parameter that is sensitive to TpA base dynamics is the resonance linewidth of the TpA adenine H2 proton. This resonance experiences a temperature-dependent broadening due to conformational dynamics. Unusual dynamics at TpA steps were originally observed in the sequence context (T)pTpTpApAp(A). We have since shown that the evidence for TpA dynamics persists when either the thymine preceding the TpA step or the adenine following the TpA step is preserved [McAteer et al., Nucleic Acids Res. 23, 3962-3966 (1995)]. Here, in order establish whether or not exceptional TpA dynamics occurs in all DNA sequence contexts, we investigated a series of DNA sequences of the form GCNaTANbNbTANaGC, where N=A,T,C,G. In this family of sequences, all 16 possible immediate sequence context environments of the form NaTANb were examined using 10 DNA sequences. Our NMR results show that the TpA adenine H2 resonance contains a temperature dependent excess linewidth indicative of dynamics in all 16 sequence context environments. By studying a complete set of sequence contexts, it was possible to recognize trends relating resonance parameters and sequence environment. For example, the magnitude of the maximum linewidth is largely determined by the identity of the nucleotide following the TpA step and the magnitude of the linewidth maximum is moderately correlated (r=0.56) with the temperature of the linewidth maximum. The physical basis for these correlations is discussed.     

Flexibility of the B-DNA backbone: effects of local and neighbouring sequences on pyrimidine-purine steps.

The structurally correlated dihedral angles epsilon and zeta are known for their large variability within the B-DNA backbone. We have used molecular modelling to study both energetic and mechanical features of these variations which can produce BI/BII transitions. Calculations were carried out on DNA oligomers containing either YpR or RpY dinucleotides steps within various sequence environments. The results indicate that CpA and CpG steps favour the BI/BII transition more than TpA or any RpY step. The stacking energy and its intra- and inter-strand components explain these effects. Analysis of neighbouring base pairs reveals that BI/BII transitions of CpG and CpA are easiest within (Y)n(R)n sequences. These can also induce a large vibrational amplitude for TpA steps within the BI conformation.     

Calculated distortions of duplex DNA by a cis, syn cyclobutane thymine dimer are unaffected by a 3' TpA step.

Molecular dynamics simulations were performed on the duplex DNA dodecamers d(CGCGAA TT CGCG): d(CGCGAATTCGCG) and d(GCACGAA TT AAG): d(CTTAATTCGTGC), where TT denotes a cis, syn cyclobutane thymine dimer. The constant temperature and pressure algorithm of the AMBER 4.1 molecular-modeling package was used with explicit water and counterions, periodic boundary conditions and electrostatic interactions evaluated by the particle-mesh Ewald method. Results were analyzed by the CURVES algorithm and its implementation in DIALS and WINDOWS. Calculated distortions of DNA structure by the thymine dimer were qualitatively and quantitatively similar for the two sequences. Despite the enhanced flexibility of the native TpA dinucleotide step, major deviations from the B-DNA values of helicoidal parameters were found only at the Ap and p dinucleotide steps in both sequences. Only the AT base pairs of the two sequences that contain the 5' thymine of the dimers exhibited weakened Watson-Crick hydrogen bonds and anomalous stretching. Hence, we conclude that the pattern of structural perturbations responsible for recognition of cis, syn thymine dimers by repair enzymes is not sensitive to their sequence context.     

Sequence Recognition of DNA by Protein-Induced Conformational Transitions

The binding of proteins to specific sequences of DNA is an important feature of virtually all DNA transactions. Proteins recognize specific DNA sequences using both direct readout (sensing types and positions of DNA functional groups) and indirect readout (sensing DNA conformation and deformability). Previously we showed that the P22 c2 repressor N-terminal domain (P22R NTD) forces the central non-contacted 5'-ATAT-3' sequence of the DNA operator into the B′ state, a state known to affect DNA hydration, rigidity and bending. Usually the B′ state, with a narrow minor groove and a spine of hydration, is reserved for A-tract DNA (TpA steps disrupt A-tracts). Here, we have co-crystallized P22R NTD with an operator containing a central 5′-ACGT-3′ sequence in the non-contacted region. C·G base pairs have not previously been observed in the B′ state and are thought to prevent it. However, P22R NTD induces a narrow minor groove and a spine of hydration to 5'-ACGT-3'. We observe that C·G base pairs have distinctive destabilizing and disordering effects on the spine of hydration. It appears that the reduced stability of the spine results in a higher energy cost for the B to B′ transition. The differential effect of DNA sequence on the barrier to this transition allows the protein to sense the non-contacted DNA sequence.     

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.     

Formation of cis-diamminedichloroplatinum(II) 1,2-intrastrand cross-links on DNA is flanking-sequence independent

Mapping of cis-diamminedichloroplatinum(II) (cis-DDP, cisplatin) DNA adducts over >3000 nucleotides was carried out using a replication blockage assay. The sites of inhibition of modified T4 DNA polymerase, also referred to as stop sites, were analyzed to determine the effects of local sequence context on the distribution of intrastrand cisplatin cross-links. In a 3120 base fragment from replicative form M13mp18 DNA containing 24.6% guanine, 25.5% thymine, 26.9% adenine and 23.0% cytosine, 166 individual stop sites were observed at a bound platinum/nucleotide ratio of 1-2 per thousand. The majority of stop sites (90%) occurred at G(n>2) sequences and the remainder were located at sites containing an AG dinucleotide. For all of the GG sites present in the mapped sequences, including those with Gn(>)2, 89% blocked replication, whereas for the AG sites only 17% blocked replication. These blockage sites were independent of flanking nucleotides in a sequence of N(1)G*G*N(2) where N(1), N(2) = A, C, G, T and G*G* indicates a 1,2-intrastrand platinum cross-link. The absence of long-range sequence dependence was confirmed by monitoring the reaction of cisplatin with a plasmid containing an 800 bp insert of the human telomere repeat sequence (TTAGGG)(n). Platination reactions monitored at several formal platinum/nucleotide ratios or as a function of time reveal that the telomere insert was not preferentially damaged by cisplatin. Both replication blockage and telomere-insert plasmid platination experiments indicate that cisplatin 1,2-intrastrand adducts do not form preferentially at G-rich sequences in vitro.     

Compensating bends in a 16-base-pair DNA oligomer containing a T(3)A(3) segment: A NMR study of global DNA curvature

In-phase ligated DNA containing T(n)A(n) segments fail to exhibit the retarded polyacrylamide gel electrophoresis (PAGE) migration observed for in-phase ligated A(n)T(n) segments, a behavior thought to be correlated with macroscopic DNA curvature. The lack of macroscopic curvature in ligated T(n)A(n) segments is thought to be due to cancellation of bending in regions flanking the TpA steps. To address this issue, solution-state NMR, including residual dipolar coupling (RDC) restraints, was used to determine a high-resolution structure of [d(CGAGGTTTAAACCTCG)2], a DNA oligomer containing a T3A3 tract. The overall magnitude and direction of bending, including the regions flanking the central TpA step, was measured using a radius of curvature, Rc, analysis. The Rc for the overall molecule indicated a small magnitude of global bending (Rc = 138 +/- 23 nm) towards the major groove, whereas the Rc for the two halves (72 +/- 33 nm and 69 +/- 14 nm) indicated greater localized bending into the minor groove. The direction of bending in the regions flanking the TpA step is in partial opposition (109 degrees), contributing to cancellation of bending. The cancellation of bending did not correlate with a pattern of roll values at the TpA step, or at the 5' and 3' junctions, of the T3A3 segment, suggesting a simple junction/roll model is insufficient to predict cancellation of DNA bending in all T(n)A(n) junction sequence contexts. Importantly, Rc analysis of structures refined without RDC restraints lacked the precision and accuracy needed to reliably measure bending.     

Radiation chemistry of a dinucleoside monophosphate and its sequence isomer

A combination of high-performance liquid chromatography (HPLC) and nuclear magnetic resonance (NMR) spectroscopy was used to analyze the products of X-irradiated aqueous solutions of the dinucleoside monophosphate thymidylyl(3'-5')-2'-deoxyadenosine, d(TpA), and its sequence isomer 2'-deoxyadenylyl(3'-5')thymidine, d(ApT). The products of d(TpA) include both bases and nucleotides and a variety of thymine modifications of d(TpA) including the two cis and two trans glycol stereoisomers, two cis monohydroxy derivatives, an N-formamide derivative, and the hydroxymethyl derivative. Attention is focused on using NMR spectral features to distinguish among the various stereoisomers. The radiation chemistry of d(ApT) is also explored and differences in product formation compared with d(TpA) are described, particularly the formation of two products involving modification of adenine base. The potential of the HPLC-NMR approach to the study of radiation chemistry in DNA model compounds is discussed.     

NMR investigation of Hoogsteen base pairing in quinoxaline antibiotic--DNA complexes: comparison of 2:1 echinomycin, triostin A and [N-MeCys3,N-MeCys7] TANDEM complexes with DNA oligonucleotides.

Hoogsteen base pairs have been demonstrated to occur in base pairs adjacent to the CpG binding sites in complexes of triostin A and echinomycin with a variety of DNA oligonucleotides. To understand the relationship of these unusual base pairs to the sequence specificity of these quinoxaline antibiotics, the conformation of the base pairs flanking the YpR binding sites of the 2:1 drug-DNA complexes of triostin A with [d(ACGTACGT)]2 and of the TpA specific [N-MeCys3, N-MeCys7] TANDEM with [d(ATACGTAT)]2 have been studied by 1H NMR spectroscopy. In both the 2:1 triostin A-DNA complex and the 2:1 [N-MeCys3, N-MeCys7] TANDEM-DNA complex, the terminal A.T base pairs are Hoogsteen base paired with the 5' adenine in the syn conformation. This indicates that both TpA specific and CpG specific quinoxaline antibiotics are capable of inducing Hoogsteen base pairs in DNA. However, in both 2:1 complexes, Hoogsteen base pairing is limited to the terminal base pairs. In the 2:1 triostin A complex, the internal adenines are anti and in the 2:1 [N-MeCys3, N-MeCys7] TANDEM-DNA complex, the internal guanines are anti regardless of pH, which indicates that the central base pairs of both complexes form Watson-Crick base pairs. This indicates that the sequence dependent nature of Hoogsteen base pairing is the same in TpA specific and CpG specific quinoxaline antibiotic-DNA complexes. We have calculated a low resolution three-dimensional structure of the 2triostin A-[d(ACGTACGT)]2 complex and compared it with other CpG specific quinoxaline antibiotic-DNA complexes. The role of stacking in the formation of Hoogsteen base pairs in these complexes is discussed.     

Molecular structure of the netropsin-d(CGCGATATCGCG) complex: DNA conformation in an alternating AT segment

The molecular structure of the complex between a minor groove binding drug (netropsin) and the DNA dodecamer d(CGCGATATCGCG) has been solved and refined by single-crystal X-ray diffraction analysis to a final R factor of 20.0% to 2.4-A resolution. The crystal is similar to that of the other related dodecamers with unit cell dimensions of a = 25.48 A, b = 41.26 A, and c = 66.88 A in the space group P2(1)2(1)2(1). In the complex, netropsin binds to the central ATAT tetranucleotide segment in the narrow minor groove of the dodecamer B-DNA double helix as expected. However, in the structural refinement the drug is found to fit the electron density in two orientations equally well, suggesting the disordered model. This agrees with the results from solution studies (chemical footprinting and NMR) of the interactions between minor groove binding drugs (e.g., netropsin and distamycin A) and DNA. The stabilizing forces between drug and DNA are provided by a combination of ionic, van der Waals, and hydrogen-bonding interactions. No bifurcated hydrogen bond is found between netropsin and DNA in this complex due to the unique dispositions of the hydrogen-bond acceptors (N3 of adenine and O2 of thymine) on the floor of the DNA minor groove. Two of the four AT base pairs in the ATAT stretch have low propeller twist angles, even though the DNA has a narrow minor groove. Alternating helical twist angles are observed in the ATAT stretch with lower twist in the ApT steps than in the TpA step.     

Role of DNA flexibility in sequence-dependent activity of uracil DNA glycosylase

Uracil DNA glycosylase (UDG) is a base excision repair enzyme that specifically recognizes and removes uracil from double- or single-stranded DNA. The efficiency of the enzyme depends on the DNA sequence surrounding the uracil. Crystal structures of UDG in complex with DNA reveal that the DNA is severely bent and distorted in the region of the uracil. This suggests that the sequence-dependent efficiency of the enzyme may be related to the energetic cost of DNA distortion in the process of specific damage recognition. To test this hypothesis, molecular dynamics simulations were performed on two sequences representing extreme cases of UDG efficiency, AUA/TAT (high efficiency) and GUG/CAC (low efficiency). Analysis of the simulations shows that the effective bending force constants are lower for the AUA/TAT sequence, indicating that this sequence is more flexible than the GUG/CAC sequence. Fluorescence lifetimes of the adenine analogue 2-aminopurine (2AP), replacing adenine opposite the uracil, are shorter in the context of the AUA/TAT sequence, indicating more dynamic base-base interaction and greater local flexibility than in the GUG/CAC sequence. Furthermore, the K(M) of Escherichia coli UDG for the AUA/TAT sequence is 10-fold smaller than that for the GUG/CAC sequence, while the k(cat) is only 2-fold smaller. This indicates that differences in UDG efficiency largely arise from differences in binding and not catalysis. These results link directly flexibility near the damaged DNA site with the efficiency of DNA repair.     

Identification of intrinsic dynamics in a DNA sequence preferentially cleaved by topoisomerase II enzyme

Topoisomerase II enzymes are essential enzymes that modulate DNA topology and play a role in chromatin compaction. While these enzymes appear to recognize and cleave the DNA in a nonrandom fashion, factors that underlie enzyme specificity remain an enigma. To gain new insights on these topics, we undertake, using NMR and molecular dynamics methods, studies of the structural and dynamic features of a 21 bp DNA segment preferentially cleaved by topoisomerases II. The large size of the oligonucleotide did not hamper the determination of structures of sufficient quality, and numerous interesting correlations between helicoidal parameters already depicted in crystals and molecular dynamics simulations are recovered here. The main feature of the sequence is the occurrence of a large opening of the base pairs in a four-residue AT-rich region located immediately at the 5′ end of one of the cleaved sites. This opening seems to be largely dependent on sequence context, since a similar opening is not found in the other AT base pairs of the sequence. Furthermore, two adenine nucleotides of the same portion of the oligonucleotide present slow internal motions at the NMR timescale, revealing particular base dynamics. In conclusion, this AT-rich region presents the most salient character in the sequence and could be involved in the preferential cleavage by topoisomerase II. The examination of preferred sites in the literature pointed out the frequent occurrence of AT-rich sequences, namely matrix attachment region and scaffold attachment region sequences, at the sites cleaved by topoisomerase II. We could infer that the particular flexibility of these sequences plays an important role in enabling the formation of a competent cleavage complex. The sequences could then be selected based on their facility to undertake conformational change during the complex formation, rather than purely based on binding affinity.     

DNA sequence recognition by the antitumor antibiotic CC-1065 and analogs of CC-1065

The DNA base pair preferences of the antitumor antibiotic CC-1065 and two analogs of CC-1065 were studied by following the rate of covalent bond formation (N-3 adenine adduct) with DNA oligomers containing the 5'NNTTA* and 5'NNAAA* sequences (N = nucleotide, A* = alkylated adenine). The rate of adduct formation of CC-1065 is greatly affected by DNA base changes at the fourth and fifth positions of the bonding site for the 5'NNAAA sequences, but not the 5'NNTTA sequences. However, an analog of CC-1065 containing the same alkylating moiety as CC-1065, but not the third fused ring system or additional methylene and oxygen substituents, shows similar rates of adduct formation for all sequences. A second analog of CC-1065 containing three fused ring systems, but not the methylene and oxygen substituents of CC-1065, shows rates of adduct formation with the same sequence dependence as CC-1065, but does not distinguish between the sequences to the degree shown by CC-1065. Adduct formation of CC-1065, but not the analogs, competes with a reversibly bound species. Thymine bases to the 3' side of a potentially reactive adenine or a cytosine base at the fifth position from the bonding adenine create reversible binding sites which decrease the rate of adduct formation of CC-1065. The sequence 5'GCGAATT binds CC-1065 only reversibly. This sequence can compete for CC-1065 with covalent bonding sequences if the sites are located in different oligomers, or if the sites are located (overlapped or not overlapped) in the same oligomer. The results of these competitive binding experiments suggest that the transfer of CC-1065 from the reversible binding site to the covalent bonding site with both sites located on a single DNA duplex, not overlapped, occurs through an equilibrium of CC-1065 in solution, not by migration of CC-1065 in the minor groove.     

The Role of Context-Dependent Mutations in Generating Compositional and Codon Usage Bias in Grass Chloroplast DNA

Abstract The influence of local base composition on mutations in chloroplast DNA (cpDNA) is studied in detail and the resulting, empirically derived, mutation dynamics are used to analyze both base composition and codon usage bias. A 4 × 4 substitution matrix is generated for each of the 16 possible flanking base combinations (contexts) using 17,253 noncoding sites, 1309 of which are variable, from an alignment of three complete grass chloroplast genome sequences. It is shown that substitution bias at these sites is correlated with flanking base composition and that the A+T content of these flanking sites as well as the number of flanking pyrimidines on the same strand appears to have general influences on substitution properties. The context-dependent equilibrium base frequencies predicted from these matrices are then applied to two analyses. The first examines whether or not context dependency of mutations is sufficient to generate average compositional differences between noncoding cpDNA and silent sites of coding sequences. It is found that these two classes of sites exist, on average, in very different contexts and that the observed mutation dynamics are expected to generate significant differences in overall composition bias that are similar to the differences observed in cpDNA. Context dependency, however, cannot account for all of the observed differences: although silent sites in coding regions appear to be at the equilibrium predicted, noncoding cpDNA has a significantly lower A+T content than expected from its own substitution dynamics, possibly due to the influence of indels. The second study examines the codon usage of low-expression chloroplast genes. When context is accounted for, codon usage is very similar to what is predicted by the substitution dynamics of noncoding cpDNA. However, certain codon groups show significant deviation when followed by a purine in a manner suggesting some form of weak selection other than translation efficiency. Overall, the findings indicate that a full understanding of mutational dynamics is critical to understanding the role selection plays in generating composition bias and sequence structure.     

Structural bioinformatics of DNA: a web-based tool for the analysis of molecular dynamics results and structure prediction

We report here the release of a web-based tool (MDDNA) to study and model the fine structural details of DNA on the basis of data extracted from a set of molecular dynamics (MD) trajectories of DNA sequences involving all the unique tetranucleotides. The dynamic web interface can be employed to analyze the first neighbor sequence context effects on the 10 unique dinucleotide steps of DNA. Functionality is included to build all atom models of any user-defined sequence based on the MD results. The backend of this interface is a relational database storing the conformational details of DNA obtained in 39 different MD simulation trajectories comprising all the 136 unique tetranucleotide steps. Examples of the use of this data to predict DNA structures are included. Availability: http://humphry.chem.wesleyan.edu:8080/MDDNA. Supplementary information: Supplementary data including color figures are available at Bioinformatics online.     

Invariant and variable base stacking geometries in B-DNA and A-DNA

We calculated the interatomic distances between all couples of non-hydrogen atoms belonging to the neighboring Watson-Crick base pairs in the available crystal structures of DNA. Their standard deviations revealed remarkably large differences in the variability of the base stacking geometries of the particular steps. In line with experimental studies in solution, (CpA)-(TpG) and (TpA).(TpA) were identified as the most variable or flexible steps in the crystal structures of B-DNA. On the other hand, base stacking geometries of the (ApT).(ApT) steps were the most invariant, which was very surprising because all three steps composed only of C and G were much more flexible. This finding suggests that conformational stability of DNA and the rigidity have different origins. Furthermore, the nucleotide sequence dependence of the flexibility was almost reversed in A-DNA because the most flexible steps in B-DNA were the least flexible in A-DNA. The most invariant steps of B-DNA were variable in A-DNA. The (ApT).(ApT) step was a notable exception to this rule because it belonged to the most rigid steps in both B-DNA and A-DNA. The present results are fully consistent with the properties that poly(dA-dT).poly(dA-dT), poly(dA).poly(dT), poly(dAdC).poly(dG-dT) and poly(dA-dG).poly(dC-dT) exhibit in solution.     

Water and ion binding around r(UpA)12 and d(TpA)12 oligomers--comparison with RNA and DNA (CpG)12 duplexes

The structural and dynamic properties of the water and ion first coordination shell of the r(A-U) and d(A-T) base-pairs embedded within the r(UpA)12 and d(TpA)12 duplexes are described on the basis of two 2.4 ns molecular dynamics simulations performed in a neutralizing aqueous environment with 0.25 M added KCl. The results are compared to previous molecular dynamics simulations of the r(CpG)12 and d(CpG)12 structures performed under similar conditions. It can be concluded that: (i) RNA helices are more rigid than DNA helices of identical sequence, as reflected by the fact that RNA duplexes keep their initial A-form shape while DNA duplexes adopt more sequence-specific shapes. (ii) Around these base-pairs, the water molecules occupy 21 to 22 well-defined hydration sites, some of which are partially occupied by potassium ions. (iii) These hydration sites are occupied by an average of 21.9, 21.0, 20.1, and 19.8 solvent molecules (water and ions) around the r(G=C), r(A-U), d(G=C), and d(A-T) pairs, respectively. (iv) From a dynamic point of view, the stability of the hydration shell is the strongest for the r(G=C) pairs and the weakest for the d(A-T) pairs. (v) For RNA, the observed long-lived hydration patterns are essentially non-sequence dependent and involve water bridges located in the deep groove and linking OR atoms of adjacent phosphate groups. Maximum lifetimes are close to 400 ps. (vi) In contrast, for DNA, long-lived hydration patterns are sequence dependent and located in the minor groove. For d(CpG)12, water bridges linking the (G)N3 and (C)O2 with the O4' atoms of adjacent nucleotides with 400 ps maximum lifetimes are characterized while no such bridges are observed for d(TpA)12. (vii) Potassium ions are observed to bind preferentially to deep/major groove atoms at RpY steps, essentially d(GpC), r(GpC), and r(ApU), by forming ion-bridges between electronegative atoms of adjacent base-pairs. On average, about half an ion is observed per base-pair. Positive ion-binding determinants are related to the proximity of two or more electronegative atoms. Negative binding determinants are associated with the electrostatic and steric hindrance due to the proximity of electropositive amino groups and neutral methyl groups. Potassium ions form only transient contacts with phosphate groups.     

Induced fit DNA recognition by a minor groove binding analogue of Hoechst 33258: fluctuations in DNA A tract structure investigated by NMR and molecular dynamics simulations

NMR analysis and molecular dynamics simulations of d(GGTAATTACC)₂ and its complex with a tetrahydropyrimidinium analogue of Hoechst 33258 suggest that DNA minor groove recognition in solution involves a combination of conformational selection and induced fit, rather than binding to a preorganised site. Analysis of structural fluctuations in the bound and unbound states suggests that the degree of induced fit observed is primarily a consequence of optimising van der Waals contacts with the walls of the minor groove resulting in groove narrowing through: (i) changes in base step parameters, including increased helical twist and propeller twist; (ii) changes to the sugar–phosphate backbone conformation to engulf the bound ligand; (iii) suppression of bending modes at the TpA steps. In contrast, the geometrical arrangement of hydrogen bond acceptors on the groove floor appears to be relatively insensitive to DNA conformation (helical twist and propeller twist). We suggest that effective recognition of DNA sequences (in this case an A tract structure) appears to depend to a significant extent on the sequence being flexible enough to be able to adopt the geometrically optimal conformation compatible with the various binding interactions, rather than involving ‘lock and key’ recognition.     

Slow motion in the CAA*TTG sequence of a DNA decamer duplex studied by NMR

In DNA duplexes, pyrimidine-purine steps are believed to be flexible or conformationally unstable. Indeed, several DNA crystal structures exhibit a multitude of conformations for CpA*TpG steps. The question arises of whether this structural flexibility is accompanied by dynamical flexibility, i.e., a question pertaining to the energy barrier between conformations. Except for TpA steps, slow motions on the microsecond-to-millisecond time scale have not been detected in duplexes until now. In the present study, such slow motion was investigated by 1H, 13C, and 15N NMR relaxation measurements on a DNA decamer d(CATTTGCATC)*d(GATGCAAATG). The DNA decamer was enriched with 15% 13C and 98% 15N isotopes for each adenosine and guanosine residue. Three lines of evidence support the notion of slow motion in the CAA*TTG moiety. Analysis of (15)N relaxation showed that the order parameter, S2, of guanosine imino NH groups was about 0.8, similar to that of CH groups for this oligomer. The strong temperature dependence of guanosine NH S2 in the CAA*TTG sequence indicated the presence of a large-amplitude motion. Signals of adenosine H8 protons in the CAA*TTG sequence were broadened in 2D 1H NOESY spectra, which also suggested the existence of slow motion. As well as being smaller than for other adenine residues, the 1H T2 values exhibited a magnetic field strength dependence for all adenosine H8 signals in the ATTTG*CAAAT region, suggesting slow motions more pronounced at the first adenosine in the CAA*TTG sequence but extending over the CAAAT*ATTTG region. This phenomenon was further examined by the pulse field strength dependence of the 1H, 13C, and 15N T1rho values. 1H and 13C T1rho values showed a pulse field strength dependence, but 15N T1rho did not. Assuming a two-site exchange process, an exchange time constant of 20-300 micros was estimated for the first adenosine in the CAA sequence. The exact nature of this motion remains unknown.     

Role of stacking interactions in the binding sequence preferences of DNA bis-intercalators: insight from thermodynamic integration free energy simulations

The major structural determinant of the preference to bind to CpG binding sites on DNA exhibited by the natural quinoxaline bis-intercalators echinomycin and triostin A, or the quinoline echinomycin derivative, 2QN, is the 2-amino group of guanine (G). However, relocation of this group by means of introduction into the DNA molecule of the 2-aminoadenine (=2,6-diaminopurine, D) base in place of adenine (A) has been shown to lead to a drastic redistribution of binding sites, together with ultratight binding of 2QN to the sequence DTDT. Also, the demethylated triostin analogs, TANDEM and CysMeTANDEM, which bind with high affinity to TpA steps in natural DNA, bind much less tightly to CpI steps, despite the fact that both adenosine and the hypoxanthine-containing nucleoside, inosine (I), provide the same hydrogen bonding possibilities in the minor groove. To study both the increased binding affinity of 2QN for DTDT relative to GCGC sites and the remarkable loss of binding energy between CysMeTANDEM and ICIC compared with ATAT, a series of thermodynamic integration free energy simulations involving conversions between DNA base pairs have been performed. Our results demonstrate that the electrostatic component of the stacking interactions between the heteroaromatic rings of these compounds and the bases that make up the intercalation sites plays a very important role in the modulation of their binding affinities.