Structure determination and analysis of helix parameters in the DNA decamer d(CATGGCCATG)2 comparison of results from NMR and crystallography

Stiffer double-stranded DNA in two-dimensional confinement due to bending anisotropy

Using analytical approach and Monte Carlo (MC) simulations, we study the elastic behavior of the intrinsically twisted elastic ribbons with bending anisotropy, such as double-stranded DNA (dsDNA), in two-dimensional (2D) confinement. We show that, due to the bending anisotropy, the persistence length of dsDNA in 2D conformations is always greater than three-dimensional (3D) conformations. This result is in consistence with the measured values for DNA persistence length in 2D and 3D in equal biological conditions. We also show that in two dimensions, an anisotropic, intrinsically twisted polymer exhibits an implicit twist-bend coupling, which leads to the transient curvature increasing with a half helical turn periodicity along the bent polymer.

DNA Recognition by a Novel Bis-Intercalator, Potent Anticancer Drug XR5944

XR5944 is a potent anticancer drug with a novel DNA binding mode: DNA bisintercalationg with major groove binding. XR5944 can bind the estrogen response element (ERE) sequence to block ER-ERE binding and inhibit ERα activities, which may be useful for overcoming drug resistance to currently available antiestrogen treatments. This review discusses the progress relating to the structure and function studies of specific DNA recognition of XR5944. The sites of intercalation within a native promoter sequence appear to be different from the ideal binding site and are context- and sequence- dependent. The structural information may provide insights for rational design of improved EREspecific XR5944 derivatives, as well as of DNA bis-intercalators in general.

Solution structure of a 2:1 complex of anticancer drug XR5944 with TFF1 estrogen response element: insights into DNA recognition by a bis-intercalator

XR5944, a deoxyribonucleic acid (DNA) bis-intercalator with potent anticancer activity, can bind the estrogen response element (ERE) sequence to inhibit estrogen receptor-α activities. This novel mechanism of action may be useful for overcoming drug resistance to currently available antiestrogen treatments, all of which target the hormone-receptor complex. Here we report the nuclear magnetic resonance solution structure of the 2:1 complex of XR5944 with the naturally occurring TFF1-ERE, which exhibits important and unexpected features. In both drug–DNA complexes, XR5944 binds strongly at one intercalation site but weakly at the second site. The sites of intercalation within a native promoter sequence appear to be context and sequence dependent. The binding of one drug molecule influences the binding site of the second. Our structures underscore the fact that the DNA binding of a bis-intercalator is directional and different from the simple addition of two single intercalation sites. Our study suggests that improved XR5944 bis-intercalators targeting ERE may be designed through optimization of aminoalkyl linker and intercalation moieties at the weak binding sites.

Evidence from thermodynamics that DNA photolyase recognizes a solvent-exposed CPD lesion

Binding of a cis,syn-cyclobutane pyrimidine dimer (CPD) to Escherichia coli DNA photolyase was examined as a function of temperature, enzyme oxidation state, salt, and substrate conformation using isothermal titration calorimetry. While the overall ΔG° of binding was relatively insensitive to most of the conditions examined, the enthalpic and entropic terms that make up the free energy of binding are sensitive to the conditions of the experiment. Substrate binding to DNA photolyase is generally driven by a negative change in enthalpy. Electrostatic interactions and protonation are affected by the oxidation state of the required FAD cofactor and substrate conformation. The fully reduced enzyme appears to bind approximately two additional water molecules as part of substrate binding. More significantly, the experimental change in heat capacity strongly suggests that the CPD lesion must be flipped out of the intrahelical base stacking prior to binding to the protein; the DNA repair enzyme appears to recognize a solvent-exposed CPD as part of its damage recognition mechanism.

Extended weak bonding interactions in DNA: pi-stacking (base-base), base-backbone, and backbone-backbone interactions

We report on several weak interactions in nucleic acids, which, collectively, can make a nonnegligible contribution to the structure and stability of these molecules. Fragments of DNA were obtained from previously determined accurate experimental geometries and their electron density distributions calculated using density functional theory (DFT). The electron densities were analyzed topologically according to the quantum theory of atoms in molecules (AIM). A web of closed-shell bonding interactions is shown to connect neighboring base pairs in base-pair duplexes and in dinuleotide steps. This bonding underlies the well-known pi-stacking interaction between adjacent nucleic acid bases and is characterized topologically for the first time. Two less widely appreciated modes of weak closed-shell interactions in nucleic acids are also described: (i) interactions between atoms in the bases and atoms belonging to the backbone (base-backbone) and (ii) interactions among atoms within the backbone itself (backbone-backbone). These interactions include hydrogen bonding, dihydrogen bonding, hydrogen-hydrogen bonding, and several other weak closed-shell X-Y interactions (X, Y = O, N, C). While each individual interaction is very weak and typically accompanied by perhaps 0.5-3 kcal/mol, the sum total of these interactions is postulated to play a role in stabilizing the structure of nucleic acids. The Watson-and-Crick hydrogen bonding is also characterized in detail at the experimental geometries as a prelude to the discussion of the modes of interactions listed in the title.

DNA structural similarity in the 2:1 complexes of the antitumor drugs trabectedin (Yondelis) and chromomycin A3 with an oligonucleotide sequence containing two adjacent TGG binding sites on opposing strands

Yondelis (trabectedin) is an antitumor ecteinascidin that binds covalently to the 2-amino group of the central guanine in the minor groove of selected DNA pyrimidine-G-G and purine-G-C triplets. Chromomycin A3 is an aureolic acid derivative that binds noncovalently to the DNA minor groove in G/C-rich triplet sites as a metal-chelated dimer. Despite their different binding modes, the cytotoxicity profiles of these two drugs, as assessed in the COMPARE analysis carried out by the National Cancer Institute on data from 60 human tumor cell lines, are highly correlated (Pearson's correlation coefficient of 0.96). We now report that in an oligonucleotide containing the "natural bending element" TGGCCA, the structural distortions inflicted by the tail-to-tail bonding of two trabectedin molecules to adjacent target sites on opposing strands are strikingly similar to those observed in a crystal containing d(TTGGCCAA)2 and two bound chromomycin A3 molecules arranged in a head-to-tail orientation in the minor groove. In both complexes, the double helix is characterized by being considerably unwound and possessing a notably widened minor groove. Binding of the drugs to this sequence could be favored by the distinct bends at each of the TpG steps that are already present in the free oligonucleotide. Simultaneous drug binding to the two strands in the manner described here is proposed to stabilize the helical structure of duplex DNA to prevent or hamper strand separation and stall replication and transcription forks.

Structure of DNA sequence d-TGATCA by two-dimensional nuclear magnetic resonance spectroscopy and restrained molecular dynamics

The 5' d-TpG 3' element is a part of DNA sequences involved in regulation of gene expression and is also a site for intercalation of several anticancer drugs. Solution conformation of DNA duplex d-TGATCA containing this element has been investigated by two-dimensional NMR spectroscopy. Using a total of 12 torsional angles and 121 distance constraints, structural refinement has been carried out by restrained molecular dynamics (rMDs) in vacuum up to 100 ps. The structure is characterized by a large positive roll at TpG/CpA base pair step and large negative propeller twist for AT and TA base pairs. The backbone torsional angle, gamma(O5'-C5'-C4'-C3'), of T1 residue adopts a trans-conformation which is corroborated by short intra nucleotide T1H6-T1H5' (3.7A) distance in nuclear overhauser effect spectroscopy (NOESY) spectra while the backbone torsional angle, beta(P-O5'-C5'-C4'), exists in trans as well as gauche state for T1 and C5 residues. There is evidence of significant flexibility of the sugar-phosphate backbone with rapid inter-conversion between two different conformers at TpG/CpA base pair step. The base sequence dependent variations and local structural heterogeneity have important implications in specific recognition of DNA by ligands.

Novel DNA bis-intercalation by MLN944, a potent clinical bisphenazine anticancer drug

The new bisphenazine anticancer drug MLN944 is a novel cytotoxic agent with exceptional anti-tumor activity against a range of human and murine tumor models both in vitro and in vivo. MLN944 has recently entered Phase I clinical trials. Despite the structural similarity with its parent monophenazine carboxamide and acridine carboxamide anticancer compounds, MLN944 appears to work by a distinct mechanism of inhibiting DNA transcription rather than the expected mechanism of topoisomerase I and II inhibition. Here we present the first NMR structure of MLN944 complexed with d(ATGCAT)(2) DNA duplex, demonstrating a novel binding mode in which the two phenazine rings bis-intercalate at the 5'-TpG site, with the carboxamide amino linker lying in the major groove of DNA. The MLN944 molecule adopts a significantly unexpected conformation and side chain orientation in the DNA complex, with the N10 on the phenazine ring protonated at pH 7. The phenazine chromophore of MLN944 is very well stacked with the flanking DNA base pairs using the parallel base-stacking intercalation binding mode. The DNA sequence specificity and the groove recognition of MLN944 binding is determined by several site-specific hydrogen bond interactions with the central G:C base pair as well as the favorable stacking interactions with the 5'-flanking thymine. The specific binding site of MLN944 is known to be recognized by a number of important transcription factors. Our electrophoretic gel mobility shift assay results demonstrated that the c-Jun DNA binding to the AP-1 site is significantly inhibited by MLN944 in a dose-dependent manner. Thus, the exceptional biological activity of MLN944 may be due to its novel DNA binding mode leading to a unique mechanism of action.

The structural basis of DNA flexibility

Although the average physico-chemical properties of a long DNA molecule may approximate to those of a thin isotropic homogeneous rod, DNA behaves more locally as an anisotropic heterogeneous rod. This bending anisotropy is sequence dependent and to a first approximation reflects both the geometry and stability of individual base steps. The biological manipulation and packaging of the molecule often depend crucially on local variations in both bending and torsional flexibility. However, whereas the probability of DNA untwisting can be strongly correlated with a high bending flexibility, DNA bending, especially when the molecule is tightly wrapped on a protein surface, may be energetically favoured by a less flexible sequence whose preferred configuration conforms more closely to that of the complementary protein surface. In the latter situation the lower bending flexibility may be more than compensated for on binding by a reduced required deformation energy relative to a fully isotropic DNA molecule.

Structure of DNA Hexamer Sequence d-CGATCG by Two-dimensional Nuclear Magnetic Resonance Spectroscopy and Restrained Molecular Dynamics

Solution conformation of self-complementary DNA duplex d-CGATCG, containing 5' d-CpG 3' site for intercalation of anticancer drug, daunomycin and adriamycin, has been investigated by nuclear magnetic resonance (NMR) spectroscopy. Complete resonance assignments of all the protons (except some H5'/H5" protons) have been obtained following standard procedures based on double quantum filtered correlation spectroscopy (dQF COSY) and two-dimensional nuclear Overhauser effect (NOE) spectra. Analysis of sums of coupling constants in one-dimensional NMR spectra, cross peak patterns in dQF COSY spectra and inter proton distances shows that the DNA sequence assumes a conformation close to the B-DNA family. The deoxyribose sugar conformation is in dynamic equilibrium with predominantly S-type conformer and a minor N-type conformer with N<-->S equilibrium varying with temperature. At 325 K, the mole fraction of the N-conformer increases for some of the residues by approximately 9%. Using a total of 10 spin-spin coupling constants and 112 NOE intensities, structural refinement has been carried out using Restrained Molecular Dynamics (rMD) with different starting structures, potential functions and rMD protocols. It is observed that pseudorotation phase angle of deoxyribose sugar for A3 and T4 residues is approximately 180 degrees and approximately 120 degrees, respectively while all other residues are close to C2'endo-conformation. A large propeller twist (approximately -18 degrees) and smallest twist angle (approximately 31 degrees) at A3pT4 step, in the middle of the sequence, a wider (12 A) and shallower (3.0 A) major groove with glycosidic bond rotation as high anti at both the ends of hexanucleotide are observed. The structure shows base-sequence dependent variations and hence strong local structural heterogeneity, which may have implications in ligand binding.

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 fluorescence 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.

The S. cerevisiae architectural HMGB protein NHP6A complexed with DNA: DNA and protein conformational changes upon binding

NHP6A is a non-sequence-specific DNA-binding protein from Saccharomyces cerevisiae which belongs to the HMGB protein family. Previously, we have solved the structure of NHP6A in the absence of DNA and modeled its interaction with DNA. Here, we present the refined solution structures of the NHP6A-DNA complex as well as the free 15bp DNA. Both the free and bound forms of the protein adopt the typical L-shaped HMGB domain fold. The DNA in the complex undergoes significant structural rearrangement from its free form while the protein shows smaller but significant conformational changes in the complex. Structural and mutational analysis as well as comparison of the complex with the free DNA provides insight into the factors that contribute to binding site selection and DNA deformations in the complex. Further insight into the amino acid determinants of DNA binding by HMGB domain proteins is given by a correlation study of NHP6A and 32 other HMGB domains belonging to both the DNA-sequence-specific and non-sequence-specific families of HMGB proteins. The resulting correlations can be rationalized by comparison of solved structures of HMGB proteins.

An FTIR and CD study of the structural effects of G-tract length and sequence context on DNA conformation in solution

Fourier transform infrared (FTIR) and CD spectroscopy have been used to investigate the structural effects of G-tract length and flanking sequence on the conformation of DNA G-tracts in aqueous solution. Particularly, a possible predisposition for A-form features has been probed, since this may be important for protein-DNA interactions. Five different G-tract-containing DNA duplexes have been studied: d[CATGGCCATG](2), d[CATGGGCCCATG](2), d[CATGGGGCCCCATG](2,) d[AGGGGCCCCT](2), and d[TGGGGCCCCA](2). In addition, a DNA duplex lacking a G-tract center was probed (d[CATATGCATATG](2)). The CD and FTIR results show that the G-tract-containing sequences are all in a dominating B-DNA conformation in solution. However, certain spectral variations reflect structural effects of sequence context and G-tract length. CD spectra and FTIR results in the 1800-1500 cm(-1) region show that the base-stacking pattern is greatly affected by the sequence context. The FTIR backbone 1250-1000 cm(-1) region shows the antisymmetric non-bridging phosphate vibration around 1225 cm(-1) in all sequences, demonstrating the overall B-conformation of the backbone. The FTIR sugar 900-800 cm(-1) region shows variable contributions of two bands around 865 cm(-1) and 840 cm(-1), reflecting the N and S-type of sugar pucker. The relative intensities of the 865 cm(-1) and 840 cm(-1) bands have been proposed in the literature to quantitatively yield the contribution of N and S-type of sugar pucker, respectively. This correlation is supported by the present study. Furthermore, the contributions of N-type sugar in the DNA sequences studied indicate structural propensities that agree with trends in reported crystal structures of the same sequences: (1) d[CATGGCCATG](2), for which FTIR shows the lowest contribution of N-type sugar puckering in solution, crystallizes in a B-like conformation; (2) d[AGGGGCCCCT](2), with the highest degree of N-type sugar puckering of all the sequences studied, crystallizes in an A-like conformation; (3) d[CATGGGCCCATG](2), with an N-type contribution intermediate between that of d[CATGGCCATG](2) and d[AGGGGCCCCT](2), crystallizes in an A/B intermediate conformation.

Localization of Fe(2+) at an RTGR sequence within a DNA duplex explains preferential cleavage by Fe(2+) and H2O2

Nicking of duplex DNA by the iron-mediated Fenton reaction occurs preferentially at a limited number of sequences. Of these, purine-T-G-purine (RTGR) is of particular interest because it is a required element in the upstream regulatory regions of many genes involved in iron and oxidative-stress responses. In order to study the basis of this preferential nicking, NMR studies were undertaken on the RTGR-containing duplex oligonucleotide, d(CGCGATATGACACTAG)/d(CTAGTGTCATATCGCG). One-dimensional and two-dimensional 1H NMR measurements show that Fe(2+) interacts preferentially and reversibly at the ATGA site within the duplex at a rate that is rapid relative to the chemical-shift timescale, while selective paramagnetic NMR line-broadening of the ATGA guanine H8 suggests that Fe(2+) interacts with the guanine N7 moiety. Localization at this site is supported by Fe(2+) titrations of a duplex containing a 7-deazaguanine substitution in place of the guanine in the ATGA sequence. The addition of a 100-fold excess of Mg(2+) over Fe(2+) does not affect the Fe(2+)-dependent broadening. When the ATGA site in the duplex is replaced by ATGT, an RTGR site (GTGA) is created on the opposite strand. Preferential iron localization then takes place at the 3' guanine in GTGA but no longer at the guanine in ATGT, consistent with the lack of preferential cleavage of ATGT sites relative to ATGA sites.

Solution structure of the dodecamer d-(CATGGGCC-CATG)2 is B-DNA. Experimental and molecular dynamics study

The DNA duplex d-(CATGGGCCCATG)2 has been studied in solution by FTIR, NMR and CD. The experimental approaches have been complemented by series of large-scale unrestrained molecular dynamics simulation with explicit inclusion of solvent and counterions. Typical proton-proton distances extracted from the NMR spectra and the CD spectra are completely in agreement with slightly modified B-DNA. By molecular dynamics simulation, starting from A-type sugar pucker, a spontaneous repuckering to B-type sugar pucker was observed. Both experimental and theoretical approaches suggest for the dodecamer d-(CATGGGCCCATG)2 under solution conditions puckering of all 2'-deoxyribose residues in the south conformation (mostly C2'-endo) and can exclude significant population of sugars in the north conformation (C3'-endo). NMR, FTIR and CD data are in agreement with a B-form of the dodecamer in solution. Furthermore, the duplex shows a cooperative B-A transition in solution induced by addition of trifluorethanol. This contrasts a recently published crystal structure of the same oligonucleotide found as an intermediate between B- and A-DNA where 23 out of 24 sugar residues were reported to adopt the north (N-type) conformation (C3'-endo) like in A-DNA (Ng, H. L., Kopka, M. L. and Dickerson, R. E., Proc. Natl. Acad. Sci. U S A 97, 2035-2039 (2000)). The simulated structures resemble standard B-DNA. They nevertheless show a moderate shift towards A-type stacking similar to that seen in the crystal, despite the striking difference in sugar puckers between the MD and X-ray structures. This is in agreement with preceding MD reports noticing special stacking features of G-tracts exhibiting a tendency towards the A-type stacking supported by the CD spectra also reflecting the G-tract stacking. MD simulations reveal several noticeable local conformational variations, such as redistribution of helical twist and base pair roll between the central GpC steps and the adjacent G-tract segments, as well as a substantial helical twist variability in the CpA(TpG) steps combined with a large positive base pair roll. These local variations are rather different from those seen in the crystal.

Conformations of an adenine bulge in a DNA octamer and its influence on DNA structure from molecular dynamics simulations

Molecular dynamics simulations have been applied to the DNA octamer d(GCGCA-GAAC). d(GTTCGCGC), which has an adenine bulge at the center to determine the pathway for interconversion between the stacked and extended forms. These forms are known to be important in the molecular recognition of bulges. From a total of ~35 ns of simulation time with the most recent CHARMM27 force field a variety of distinct conformations and subconformations are found. Stacked and fully looped-out forms are in excellent agreement with experimental data from NMR and x-ray crystallography. Furthermore, in a number of conformations the bulge base associates with the minor groove to varying degrees. Transitions between many of the conformations are observed in the simulations and used to propose a complete transition pathway between the stacked and fully extended conformations. The effect on the surrounding DNA sequence is investigated and biological implications of the accessible conformational space and the suggested transition pathway are discussed, in particular for the interaction of the MS2 replicase operator RNA with its coat protein.

Kinetic analysis of the stepwise formation of a long-range DNA interstrand cross-link by a dinuclear platinum antitumor complex: evidence for aquated intermediates and formation of both kinetically and thermodynamically controlled conformers

Reported here is a detailed study of the kinetics and mechanism of formation of a 1,4 GG interstrand cross-link by [(trans-PtCl(NH(3))(2))(2)(mu-NH(2)(CH(2))(n)NH(2))](2+) (1,1/t,t (n = 6), 1), the prototype of a novel class of platinum antitumor complexes. The reaction of the self-complementary 12-mer duplex 5'-[d(ATATGTACATAT)(2)] with (15)N-1 has been studied at 298 K, pH 5.4, by [(1)H,(15)N] HSQC 2D NMR spectroscopy. Initial electrostatic interactions with the duplex are observed for 1 and the monoaqua monochloro species (2). Aquation of 1 to yield 2 occurs with a pseudo-first-order rate constant of (4.15 +/- 0.04) x 10(-5) s(-1). 2 then undergoes monofunctional binding to the guanine N7 of the duplex to form 3 (G/Cl) with a rate constant of 0.47 +/- 0.06 M(-(1) s(-1). There is an electrostatic interaction between the unbound [PtN(3)Cl] group of 3 and the duplex, which is consistent with H-bonding interactions observed in the molecular model of the monofunctional (G/Cl) adduct. Closure of 3 to form the 1,4 GG interstrand cross-link (5) most likely proceeds via the aquated (G/H(2)O) intermediate (4) (pseudo-first-order rate constant = (3.62 +/- 0.04) x 10(-5) s(-1)) followed by closure of 4 to form 5 (rate constant = (2.7 +/- 1.5) x 10(-3) s(-1)). When closure is treated as direct from 3 (G/Cl) the rate constant is (3.39 +/- 0.04) x 10(-5) s(-1). Closure is ca. 10-55-fold faster than that found for 1,2 GG intrastrand cross-link formation by the diaqua form of cisplatin. Changes in the (1)H and (15)N shifts of the interstrand cross-link 5 indicate that the initially formed conformer (5(i)) converts irreversibly into other product conformer(s) 5(f). The NMR data for 5(i) are consistent with a molecular model of the 1,4 GG interstrand cross-link on B-form DNA, which shows that the NH(2) protons have no contacts except with solvent. The NMR data for 5(f) show several distinct NH(2) environments indicative of interactions between the NH(2) protons and the DNA. HPLC characterization of the final product showed only one major product peak that was confirmed by ESI-FTICR mass spectroscopy to be a cross-linked adduct of (15)N-1 and the duplex. The potential significance of these findings to the antitumor activity of dinuclear platinum complexes is discussed.

Local conformational variations observed in B-DNA crystals do not improve base stacking: computational analysis of base stacking in a d(CATGGGCCCATG)(2) B<-->A intermediate crystal structure

The crystal structure of d(CATGGGCCCATG)(2) shows unique stacking patterns of a stable B<-->A-DNA intermediate. We evaluated intrinsic base stacking energies in this crystal structure using an ab initio quantum mechanical method. We found that all crystal base pair steps have stacking energies close to their values in the standard and crystal B-DNA geometries. Thus, naturally occurring stacking geometries were essentially isoenergetic while individual base pair steps differed substantially in the balance of intra-strand and inter-strand stacking terms. Also, relative dispersion, electrostatic and polarization contributions to the stability of different base pair steps were very sensitive to base composition and sequence context. A large stacking flexibility is most apparent for the CpA step, while the GpG step is characterized by weak intra-strand stacking. Hydration effects were estimated using the Langevin dipoles solvation model. These calculations showed that an aqueous environment efficiently compensates for electrostatic stacking contributions. Finally, we have carried out explicit solvent molecular dynamics simulation of the d(CATGGGCCCATG)(2) duplex in water. Here the DNA conformation did not retain the initial crystal geometry, but moved from the B<-->A intermediate towards the B-DNA structure. The base stacking energy improved in the course of this simulation. Our findings indicate that intrinsic base stacking interactions are not sufficient to stabilize the local conformational variations in crystals.

DNA-interactions and nuclear localisation of the chromosomal HMG domain protein SSRP1 from maize

The structure-specific recognition protein 1 (SSRP1) is a member of the protein family containing a high mobility group (HMG) domain DNA-binding motif. We have functionally characterised the 71.4 kDa Zm-SSRP1 protein from maize. The chromatin-associated Zm-SSRP1 is detected by immunoblot analysis in maize leaves, kernels and suspension culture cells, but not in roots. Mediated by its HMG domain, recombinant Zm-SSRP1 interacts structure-specifically with supercoiled DNA and DNA minicircles when compared with linear DNA. In linear duplex DNA, the protein does not recognise a specific sequence, but it binds preferentially to sequences containing the deformable dinucleotide TG, as demonstrated by a random oligonucleotide selection experiment. Zm-SSRP1 modulates DNA structure by bending the target sequence, since it promotes the circularisation of short DNA fragments in the presence of DNA ligase. Moreover, Zm-SSRP1 facilitates the formation of nucleoprotein structures, as measured using the bacterial site-specific beta-mediated recombination reaction. Analysis of the subcellular localisation of various SSRP1-GFP fusions revealed that, in contrast to HMG domain transcription factors, the nuclear localisation sequence of Zm-SSRP1 is situated within a 20-amino acid residue region adjacent to the HMG domain rather than within the DNA-binding domain. The results are discussed in the context of the likely function of SSRP1 proteins in transcription and replication.

Effect of neighboring bases on base-pair stacking orientation: a molecular dynamics study

It is generally believed that base-pair stacking interaction in DNA double helix is one of the strongest interactions that governs sequence directed structural variability. However, X-ray crystal structures of some base-paired doublet sequences have been seen to adopt different structures when flanked by different base-pairs. DNA crystal database, however, is still too small to make good statistical inference about effect of such flanking residues. Influence of neighboring residue on the local helical geometry of a base-paired doublet in B-DNA has been investigated here using molecular dynamics simulation. We have generated ensembles of structures for d(CA).d(TG) and d(AA).d(TT) base-paired doublets located at the centers of d(CGCGCAAAGCG).d(CGCTTTGCGCG) and d(CGCGAAAACGCG).d(CGCGTTTTCGCG) sequences along with their analogs by varying the bases either at 5'- or 3'- position to the central doublet. Comparison of base paired doublet parameters for the ensembles of structures show that stacking geometry of d(CA).d(TG) doublet depends on some of the flanking base-pairs. On the other hand d(AA).d(TT) doublet remains nearly unperturbed when the flanking residues are altered.

Sequence-dependent variation in DNA minor groove width dictates orientational preference of Hoechst 33258 in A-tract recognition: solution NMR structure of the 2:1 complex with d(CTTTTGCAAAAG)(2)

The solution structure of the dodecamer duplex d(CTTTTGCAAAAG)(2)and its 2:1 complex with the bis -benzimidazole Hoechst 33258 has been investigated by NMR and NOE-restrained molecular dynamics (rMD) simulations. Drug molecules are bound in each of the two A-tracts with the bulky N-methylpiperazine ring of each drug located close to the central TG (CA) step, binding essentially to the narrow minor groove of each A-tract. MD simulations over 1 ns, using an explicit solvation model, reveal time-averaged sequence-dependent narrowing of the minor groove from the 3'-end towards the 5'-end of each TTTT sequence. Distinct junctions at the TpG (CpA) steps, characterised by large positive roll, low helical and propeller twists and rapid AT base pair opening rates, add to the widening of the groove at these sites and appear to account for the bound orientation of the two drug molecules with the N-methylpiperazine ring binding in the wider part of the groove close to the junctions. Comparisons between the free DNA structure and the 2:1 complex (heavy atom RMSD 1.55 A) reveal that these sequence-dependent features persist in both structures. NMR studies of the sequence d(GAAAAGCTTTTC)(2), in which the A-tracts have been inverted with the elimination of the TpG junctions, results in loss of orientational specificity of Hoechst 33258 and formation of multiple bound species in solution, consistent with the drug binding in a number of different orientations.

Solution structure and dynamics of a complex between DNA and the antitumor bisnaphthalimide LU-79553: intercalated ring flipping on the millisecond time scale

Using a combination of nuclear magnetic resonance (NMR) spectroscopy experiments and molecular dynamics, we have analyzed the structure and dynamics of a complex between the bisnaphthalimide drug LU-79553 and the DNA duplex d(ATGCAT)(2). LU-79553 is a DNA-binding topoisomerase II inhibitor that is particularly effective against human solid tumors that are refractory to other drugs. We have found that the two naphthalimide chromophores of the drug bisintercalate at the TpG and CpA steps of the DNA hexanucleotide, stacking mainly with the purine G and A bases from opposite strands. The 3, 7-diazanonylene linker lies in the major groove of the DNA molecule, with its two amino groups hydrogen-bonded to the symmetry-related guanine bases. Unexpectedly, we have detected an unprecedented exchange process between two equivalent and intercalated states of the naphthalimide rings in the drug-DNA complex. The interconversion process takes place by rotational ring flipping, has an activation energy of 22 kcal mol(-)(1) for the two rings, and does not affect the aminoalkyl linker region of the drug. The exchange rate is intermediate to fast on the chemical shift time scale at 36 degrees C (1800 s(-)(1)) but slow at 2 degrees C (20 s(-)(1)). We have also observed limited flexibility for the drug linker on the picosecond time scale on the basis of NMR data and a time-averaged restrained molecular dynamics simulation. The implications of the structural and dynamic features of the DNA-LU-79553 complex on the binding specificity and on the antitumor activity of bisnaphthalimide agents are discussed.

Uniquely altered DNA replication fidelity conferred by an amino acid change in the nucleotide binding pocket of human immunodeficiency virus type 1 reverse transcriptase

Arginine 72 in human immunodeficiency virus type 1 reverse transcriptase (RT), a highly conserved residue among retroviral polymerases and telomerases, forms part of the binding pocket for the nascent base pair. We show here that replacement of Arg(72) by alanine strongly alters fidelity in a highly unusual manner. R72A reverse transcriptase is a frameshift and base substitution antimutator polymerase whose increased fidelity results both from increased nucleotide selectivity and from a decreased ability to extend mismatched primer termini. Thus, Arg(72)-substrate interactions in wild-type human immunodeficiency virus type 1 RT can stabilize incorrect nucleotides allowing misinsertion and promoting extension of mismatched and perhaps misaligned template-primers. In contrast to the higher fidelity at most sites, R72A RT is highly error-prone for misincorporations opposite template T in the sequence context: 5'-CTGG. Surprisingly, this results mostly from a 1200-fold increase in the apparent K(m) for correct dAMP incorporation. Thus, Arg(72) interactions with substrate are critical for the stability of the correct T.dAMP base pair when the 5'-CTGG sequence is present in the binding pocket for the nascent base pair. Collectively, the data show that a mutant polymerase may yield higher than normal average replication fidelity, yet paradoxically place specific sequences at very high risk of mutation.

Modeling helix-turn-helix protein-induced DNA bending with knowledge-based distance restraints

A crucial element of many gene functions is protein-induced DNA bending. Computer-generated models of such bending have generally been derived by using a presumed bending angle for DNA. Here we describe a knowledge-based docking strategy for modeling the structure of bent DNA recognized by a major groove-inserting alpha-helix of proteins with a helix-turn-helix (HTH) motif. The method encompasses a series of molecular mechanics and dynamics simulations and incorporates two experimentally derived distance restraints: one between the recognition helix and DNA, the other between respective sites of protein and DNA involved in chemical modification-enabled nuclease scissions. During simulation, a DNA initially placed at a distance was "steered" by these restraints to dock with the binding protein and bends. Three prototype systems of dimerized HTH DNA binding were examined: the catabolite gene activator protein (CAP), the phage 434 repressor (Rep), and the factor for inversion stimulation (Fis). For CAP-DNA and Rep-DNA, the root mean square differences between model and x-ray structures in nonhydrogen atoms of the DNA core domain were 2.5 A and 1.6 A, respectively. An experimental structure of Fis-DNA is not yet available, but the predicted asymmetrical bending and the bending angle agree with results from a recent biochemical analysis.

Solution structure of the HMG protein NHP6A and its interaction with DNA reveals the structural determinants for non-sequence-specific binding

NHP6A is a chromatin-associated protein from Saccharomyces cerevisiae belonging to the HMG1/2 family of non-specific DNA binding proteins. NHP6A has only one HMG DNA binding domain and forms relatively stable complexes with DNA. We have determined the solution structure of NHP6A and constructed an NMR-based model structure of the DNA complex. The free NHP6A folds into an L-shaped three alpha-helix structure, and contains an unstructured 17 amino acid basic tail N-terminal to the HMG box. Intermolecular NOEs assigned between NHP6A and a 15 bp 13C,15N-labeled DNA duplex containing the SRY recognition sequence have positioned the NHP6A HMG domain onto the minor groove of the DNA at a site that is shifted by 1 bp and in reverse orientation from that found in the SRY-DNA complex. In the model structure of the NHP6A-DNA complex, the N-terminal basic tail is wrapped around the major groove in a manner mimicking the C-terminal tail of LEF1. The DNA in the complex is severely distorted and contains two adjacent kinks where side chains of methionine and phenylalanine that are important for bending are inserted. The NHP6A-DNA model structure provides insight into how this class of architectural DNA binding proteins may select preferential binding sites.

High base pair opening rates in tracts of GC base pairs

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