DNA Decamer Research Articles

Direct energy functional minimization under orthogonality constraints

The direct energy functional minimization problem in electronic structure theory, where the single-particle orbitals are optimized under the constraint of orthogonality, is explored. We present an orbital transformation based on an efficient expansion of the inverse factorization of the overlap matrix that keeps orbitals orthonormal. The orbital transformation maps the orthogonality constrained energy functional to an approximate unconstrained functional, which is correct to some order in a neighborhood of an orthogonal but approximate solution. A conjugate gradient scheme can then be used to find the ground state orbitals from the minimization of a sequence of transformed unconstrained electronic energy functionals. The technique provides an efficient, robust, and numerically stable approach to direct total energy minimization in first principles electronic structure theory based on tight-binding, Hartree-Fock, or density functional theory. For sparse problems, where both the orbitals and the effective single-particle Hamiltonians have sparse matrix representations, the effort scales linearly with the number of basis functions N in each iteration. For problems where only the overlap and Hamiltonian matrices are sparse the computational cost scales as O(M2N), where M is the number of occupied orbitals. We report a single point density functional energy calculation of a DNA decamer hydrated with 4003 water molecules under periodic boundary conditions. The DNA fragment containing a cis-syn thymine dimer is composed of 634 atoms and the whole system contains a total of 12,661 atoms and 103,333 spherical Gaussian basis functions.

Rational Design of Hetero-ring-Expanded Guanine Analogs with Enhanced Properties for Modified DNA Building Blocks

The properties and modes of recognition of physiological DNAs associated with the four natural nucleobases might be extended, in principle, by the design of non-natural nucleobase derivatives. The goal is an expansion of the genetic alphabet, with the possible outcome of producing new DNAs with improved physical or biological properties. In this work, a new series of hetero-ring-expanded guanine analogs are proposed, and their relevant structural characteristics and electronic properties are determined by density functional theory. The stabilities of the decamer DNA duplexes (dn.dC)10 (where n represents the corresponding expanded guanine analog designed here) are also examined, using molecular dynamics. The simulations show that the designed motifs can form stable DNA-like structures. We determined the pairing energies for the Watson-Crick (WC) hydrogen-bonded dimers between the expanded G-analogs and the natural C, and found that the pairing energies are close to those of the natural GC pair. The calculated adiabatic ionization potentials (IPs) of the size-expanded guanine analogs and their base pairs, and the corresponding vertical ionization potentials, show that some are distinctly smaller than the corresponding natural versions. The HOMO-LUMO energy gaps for most of the size-expanded guanine analogs and their WC base pairs are considerably lower than those of the corresponding natural base and base pairs. Thus, the expanded G bases may be considered as DNA genetic motifs, and they may serve as building blocks for potential biological applications and the development of molecular electronic devices.

1H-NMR Studies of Duplex DNA Decamer Containing a Uracil Cyclobutane Dimer: Implications Regarding the High UV Mutagenecity of CC Photolesions¶

To determine the origin of the UV-specific CC to TT tandem mutation at the CC site, we made a duplex DNA decamer containing a uracil cis–syn cyclobutane dimer (CBD) as the deaminated model of a cytosine dimer. Two-dimensional 1H-NMR spectroscopy studies were performed on this sequence where two adenines (Ade) were opposite to the uracil dimer. Two imino protons of the uracil dimer were found to retain Watson–Crick hydrogen bonding with the opposite Ade, although the 5′-U(NH) of the dimer site showed unusual upfield shift like that of the 5′-T(NH) of the TT dimer, which seemed to be associated with deshielding by the flanking base rather than with reduced hydrogen bonding. (McAteer et al. 1998, J. Mol. Biol. 282:1013–1032). Hydrogen bondings at the dimer site were also supported by detecting typical strong nuclear Overhauser effects (NOE) between two imino protons and the opposite Ade H2 or NH2. But sequential NOE interactions of base protons with sugar protons were absent at the two flanking nucleotides of the 5′ side of the uracil dimer and at the intradimer site, contrasting with its thymine analog where sequential NOE was absent only at the A4–T5 step. In addition, NOE cross peak for U5(NH) ↔ A4(H2) was detected, although the NOE interactions of U6(NH) with A7(H2) and A17(H2) were not observed in contrast to the thymine dimer duplex. This different local structural alteration may be affected by the induced right-hand twisted puckering mode of cis–syn cyclobutane ring of the uracil dimer in the B-DNA duplex, even though the isolated uracil dimer had left-hand twisted puckering rigidly. In parallel, these observations may be correlated with observed differences in mutagenic properties between cis–syn UU dimer and cis–syn TT dimer.

A Pyrazolato‐Bridged Dinuclear Platinum(II) Complex Induces Only Minor Distortions upon DNA‐Binding

The cytotoxic, pyrazolato-bridged dinuclear platinum(II) complex [(cis-{Pt(NH3)2})2(mu-OH)(mu-pz)]2+ (pz=pyrazolate) has been found to cross-link two adjacent guanines of a double-stranded DNA decamer without destabilizing the duplex and without changing the directionality of the helix axis. A 1H NMR study of the oligonucleotide d(CTCTG*G*TCTC)-d(GAGACCAGAG), cross-linked at the two G* guanines by [(cis-{Pt(NH3)2})2(mu-pz)]3+, and molecular dynamics simulations of the explicitly solvated duplex were performed to characterize the structural details of the adduct. The dinuclear platinum cross-link unwinds the helix by approximately 15 degrees , that is, to a similar extent as the widely used antitumor drug cisplatin, but, in contrast to the latter, induces no significant bend in the helix axis. The Watson-Crick base-pairing remains intact, and the melting temperature of the duplex is unaffected by the cross-link. The helical twist is considerably reduced between the two platinated bases, as becomes manifest in an unusually short sequential H1'-H1' distance. This unwinding also affects the sugar ring of the guanosine in the 3'-position to the cross-link, which presents an N<-->S equilibrium. This is the first cytotoxic platinum complex that has been successfully designed by envisioning the structural consequences of its binding to DNA.

Similarities and differences in interaction of K + and Na + with condensed ordered DNA. A molecular dynamics computer simulation study

Four 20 ns molecular dynamics simulations have been performed with two counterions, K+ or Na+, at two water contents, 15 or 20 H2O per nucleotide. A hexagonal simulation cell comprised of three identical DNA decamers [d(5′-ATGCAGTCAG) × d(5′-TGACTGCATC)] with periodic boundary condition along the DNA helix was used. The simulation setup mimics the DNA state in oriented DNA fibers or in crystals of DNA oligomers. Variation of counterion nature and water content do not alter averaged DNA structure. K+ and Na+ binding to DNA are different. K+ binds to the electronegative sites of DNA bases in the major and the minor grooves, while Na+ interacts preferentially with the phosphate groups. Increase of water causes a shift of both K+ and Na+ from the first hydration shell of O1P/O2P and of the DNA bases in the minor groove with lesser influence for the cation binding to the bases in the major groove. Mobility of both water and cations in the K–DNA systems is faster than in the Na–DNA systems: Na+ organizes and immobilizes water structure around itself and near DNA while for K+ water is less organized and more dynamic.

Design and characterization of a multimeric DNA binding protein using Sac7d and GCN4 as templates

The protein Sac7d belongs to a class of small chromosomal proteins from the hyperthermophilic archaeon Sulfolobus acidocaldarius. Sac7d is extremely stable to heat, acid, and chemical agents. This protein is a monomer and it binds DNA without any particular sequence preference, while inducing a sharp kink in the DNA. By appending a leucine-zipper-like helical peptide derived from the yeast transcriptional activator GCN4 to the C-terminal end of Sac7d, the modified monomers (denoted S7dLZ) are expected to interact with each other via hydrophobic force to form a parallel dimer. The recombinant S7dLZ was expressed in Escherichia coli and purified by heating and ion-exchange chromatography. The formation of dimer was detected by gel-filtration chromatography and chemical cross-link. The results of surface plasmon resonance and circular dichroism experiments showed that the DNA-binding capacity was retained. Furthermore, X-ray diffraction analysis of single crystals of S7dLZ in complex with DNA decamer CCTATATAGG showed that the leucine-zipper segments of S7dLZ were associated into an antiparallel four-helix bundle. There are two DNA fragments bound to each S7dLZ tetramer in the crystal. This model works as a successful template that endows protein a new function without losing original properties.

Probing the Influence of Stereoelectronic Effects on the Biophysical Properties of Oligonucleotides: Comprehensive Analysis of the RNA Affinity, Nuclease Resistance, and Crystal Structure of Ten 2‘-O-Ribonucleic Acid Modifications,

The syntheses of 10 new RNA 2'-O-modifications, their incorporation into oligonucleotides, and an evaluation of their properties such as RNA affinity and nuclease resistance relevant to antisense activity are presented. All modifications combined with the natural phosphate backbone lead to significant gains in terms of the stability of hybridization to RNA relative to the first-generation DNA phosphorothioates (PS-DNA). The nuclease resistance afforded in particular by the 2'-O-modifications carrying a positive charge surpasses that of PS-DNA. However, small electronegative 2'-O-substituents, while enhancing the RNA affinity, do not sufficiently protect against degradation by nucleases. Similarly, oligonucleotides containing 3'-terminal residues modified with the relatively large 2'-O-[2-(benzyloxy)ethyl] substituent are rapidly degraded by exonucleases, proving wrong the assumption that steric bulk will generally improve protection against nuclease digestion. To analyze the factors that contribute to the enhanced RNA affinity and nuclease resistance we determined crystal structures of self-complementary A-form DNA decamer duplexes containing single 2'-O-modified thymidines per strand. Conformational preorganization of substituents, favorable electrostatic interactions between substituent and sugar-phosphate backbone, and a stable water structure in the vicinity of the 2'-O-modification all appear to contribute to the improved RNA affinity. Close association of positively charged substituents and phosphate groups was observed in the structures with modifications that protect most effectively against nucleases. The promising properties exhibited by some of the analyzed 2'-O-modifications may warrant a more detailed evaluation of their potential for in vivo antisense applications. Chemical modification of RNA can also be expected to significantly improve the efficacy of small interfering RNAs (siRNA). Therefore, the 2'-O-modifications introduced here may benefit the development of RNAi therapeutics.

Structure of human POT1 bound to telomeric single-stranded DNA provides a model for chromosome end-protection

The POT1 (protection of telomeres 1) protein binds the single-stranded overhang at the ends of chromosomes in diverse eukaryotes. It is essential for chromosome end-protection in the fission yeast Schizosaccharomyces pombe, and it is involved in regulation of telomere length in human cells. Here, we report the crystal structure at a resolution of 1.73 A of the N-terminal half of human POT1 (hPOT1) protein bound to a telomeric single-stranded DNA (ssDNA) decamer, TTAGGGTTAG, the minimum tight-binding sequence indicated by in vitro binding assays. The structure reveals that hPOT1 contains two oligonucleotide/ oligosaccharide-binding (OB) folds; the N-terminal OB fold binds the first six nucleotides, resembling the structure of the S. pombe Pot1pN-ssDNA complex, whereas the second OB fold binds and protects the 3' end of the ssDNA. These results provide an atomic-resolution model for chromosome end-capping.

Short lexitropsin that recognizes the DNA minor groove at 5'-ACTAGT-3': understanding the role of isopropyl-thiazole.

Isopropyl-thiazole ((iPr)Th) represents a new addition to the building blocks of nucleic acid minor groove-binding molecules. The DNA decamer duplex d(CGACTAGTCG)(2) is bound by a short lexitropsin of sequence formyl-PyPy(iPr)Th-Dp (where Py represents N-methyl pyrrole, (iPr)Th represents thiazole with an isopropyl group attached, and Dp represents dimethylaminopropyl). NMR data indicate ligand binding in the minor groove of DNA to the sequence 5'-ACT(5)AG(7)T-3' at a 2:1 ratio of ligand to DNA duplex. Ligand binding, assisted by the enhanced hydrophobicity of the (iPr)Th group, occurs in a head-to-tail fashion, the formyl headgroups being located toward the 5'-ends of the DNA sequence. Sequence reading is augmented through hydrogen bond formation between the exocyclic amine protons of G(7) and the (iPr)Th nitrogen, which lies on the minor groove floor. The B(I)/B(II) DNA backbone equilibrium is altered at the T(5) 3'-phosphate position to accommodate a B(II) configuration. The ligands bind in a staggered mode with respect to one another creating a six base pair DNA reading frame. The introduction of a new DNA sequence-reading element into the recognition jigsaw, combined with an extended reading frame for a small lexitropsin with enhanced hydrophobicity, holds great promise in the development of new, potentially commercially viable drug lead candidates for gene targeting.

Structures of the hyperthermophilic chromosomal protein Sac7d in complex with DNA decamers.

The protein Sac7d belongs to a class of small chromosomal proteins from the hyperthermophilic archaeon Sulfolobus acidocaldarius. Two new crystal forms of Sac7d in complex with the DNA decamers CCTATATAGG and CCTACGTAGG were obtained and their structures were determined by molecular replacement. The refined models yielded R/Rfree values of 0.221/0.257 and 0.248/0.290 at 1.9 and 2.2 A resolution, respectively. The protein structures are similar to the previously determined structure of Sac7d-GCGATCGC (PDB code 1azp), but the DNA molecules are more bent overall, by 14-20 degrees. The relative positions of the Sac7d protein and the bound DNA also differ by rotations of 6-10 degrees and translations of 1.0-2.4 A. In addition to the water molecules in the central cavity, three additional conserved water molecules are found that mediate the protein-DNA interactions. The decamer DNA fragments form virtual double helices in the crystal, with a unit length of eight base pairs. The molecular packing of the new crystal forms differs from that of 1azp. The terminal nucleotides are opened up and form triple base pairs with other DNA molecules. Through lattice contacts, the Sac7d molecule also makes additional interactions with DNA, whereas only limited protein-protein interactions are seen.

NMR structure of the DNA decamer duplex containing double T*G mismatches of cis-syn cyclobutane pyrimidine dimer: implications for DNA damage recognition by the XPC-hHR23B complex.

The cis-syn cyclobutane pyrimidine dimer (CPD) is a cytotoxic, mutagenic and carcinogenic DNA photoproduct and is repaired by the nucleotide excision repair (NER) pathway in mammalian cells. The XPC-hHR23B complex as the initiator of global genomic NER binds to sites of certain kinds of DNA damage. Although CPDs are rarely recognized by the XPC-hHR23B complex, the presence of mismatched bases opposite a CPD significantly increased the binding affinity of the XPC-hHR23B complex to the CPD. In order to decipher the properties of the DNA structures that determine the binding affinity for XPC-hHR23B to DNA, we carried out structural analyses of the various types of CPDs by NMR spectroscopy. The DNA duplex which contains a single 3' T*G wobble pair in a CPD (CPD/GA duplex) induces little conformational distortion. However, severe distortion of the helical conformation occurs when a CPD contains double T*G wobble pairs (CPD/GG duplex) even though the T residues of the CPD form stable hydrogen bonds with the opposite G residues. The helical bending angle of the CPD/GG duplex was larger than those of the CPD/GA duplex and properly matched CPD/AA duplex. The fluctuation of the backbone conformation and significant changes in the widths of the major and minor grooves at the double T*G wobble paired site were also observed in the CPD/GG duplex. These structural features were also found in a duplex that contains the (6-4) adduct, which is efficiently recognized by the XPC-hHR23B complex. Thus, we suggest that the unique structural features of the DNA double helix (that is, helical bending, flexible backbone conformation, and significant changes of the major and/or minor grooves) might be important factors in determining the binding affinity of the XPC-hHR23B complex to DNA.

Effect of bulky lesions on DNA: solution structure of a DNA duplex containing a cholesterol adduct.

The three-dimensional solution structure of two DNA decamers of sequence d(CCACXGGAAC)-(GTTCCGGTGG) with a modified nucleotide containing a cholesterol derivative (X) in its C1 '(chol)alpha or C1 '(chol)beta diastereoisomer form has been determined by using NMR and restrained molecular dynamics. This DNA derivative is recognized with high efficiency by the UvrB protein, which is part of the bacterial nucleotide excision repair, and the alpha anomer is repaired more efficiently than the beta one. The structures of the two decamers have been determined from accurate distance constraints obtained from a complete relaxation matrix analysis of the NOE intensities and torsion angle constraints derived from J-coupling constants. The structures have been refined with molecular dynamics methods, including explicit solvent and applying the particle mesh Ewald method to properly evaluate the long range electrostatic interactions. These calculations converge to well defined structures whose conformation is intermediate between the A- and B-DNA families as judged by the root mean square deviation but with sugar puckerings and groove shapes corresponding to a distorted B-conformation. Both duplex adducts exhibit intercalation of the cholesterol group from the major groove of the helix and displacement of the guanine base opposite the modified nucleotide. Based on these structures and molecular dynamics calculations, we propose a tentative model for the recognition of damaged DNA substrates by the UvrB protein.

Solution Structure of a DNA Duplex Containing an α-Anomeric Adenosine: Insights into Substrate Recognition by Endonuclease IV

The cytotoxic α anomer of adenosine, generated in situ by radicals, must be recognized and repaired to maintain genomic stability. Endonuclease IV (Endo IV), a member of the base excision repair (BER) enzyme family, in addition to acting on abasic sites, has the auxiliary function of removing this mutagenic nucleotide in Escherichia coli. We have employed enzymatic, thermodynamic, and structural studies on DNA duplexes containing a central α-anomeric adenosine residue to characterize the role of DNA structure on recognition and catalysis by Endo IV. The enzyme recognizes and cleaves our αA-containing DNA duplexes at the site of the modification. The NMR solution structure of the DNA decamer duplex establishes that the single α-anomeric adenosine residue is intrahelical and stacks in a reverse Watson–Crick fashion consistent with the slight decrease in thermostability. However, the presence of this lesion confers significant changes to the global duplex conformation, resulting from a kink of the helical axis into the major groove and an opening of the minor groove emanating from the α-anomeric site. Interestingly, the conformation of the flanking base-paired segments is not greatly altered from a B-type conformation. The global structural changes caused by this lesion place the DNA along the conformational path leading to the DNA structure observed in the complex. Thus, it appears that the α-anomeric lesion facilitates recognition by Endo IV.

Molecular dynamics simulation study of oriented polyamine- and Na-DNA: sequence specific interactions and effects on DNA structure.

Molecular dynamics (MD) computer simulations have been carried out on four systems that correspond to an infinite array of parallel ordered B-DNA, mimicking the state in oriented DNA fibers and also being relevant for crystals of B-DNA oligonucleotides. The systems were all comprised of a periodical hexagonal cell with three identical DNA decamers, 15 water molecules per nucleotide, and counterions balancing the DNA charges. The sequence of the double helical DNA decamer was d(5'-ATGCAGTCAG)xd(5'-TGACTGCATC). The counterions were the two natural polyamines spermidine(3+) (Spd(3+)) and putrescine(2+) (Put(2+)), the synthetic polyamine diaminopropane(2+) (DAP(2+)), and the simple monovalent cation Na(+). This work compares the specific structures of the polyamine- and Na-DNA systems and how they are affected by counterion interactions. It also describes sequence-specific hydration and interaction of the cations with DNA. The local DNA structure is dependent on the nature of the counterion. Even the very similar polyamines, Put(2+) and DAP(2+), show clear differences in binding to DNA and in effect on hydration and local structure. Generally, the polyamines disorder the hydration of the DNA around their binding sites whereas Na(+) being bound to DNA attracts and organizes water in its vicinity. Cation binding at the selected sites in the minor and in the major groove is compared for the different polyamines and Na(+). We conclude that the synthetic polyamine (DAP(2+)) binds specifically to several structural and sequence-specific motifs on B-DNA, unlike the natural polyamines, Spd(3+) and Put(2+). This specificity of DAP(2+) compared to the more dynamic behavior of Spd(3+) and Put(2+) may explain why the latter polyamines are naturally occurring in cells.

Crystallization and preliminary X-ray diffraction studies of DNA decamer d(CCAGGCCTGG) complexed with cobalt(III)-pepleomycin.

Crystals of the self-complementary oligonucleotide d(CCAGGCCTGG) [Heinemann & Alings (1989), J. Mol. Biol. 210, 369-381] complexed with the hydroperoxide of cobalt-pepleomycin (CoPEP) were obtained by the hanging-drop vapour-diffusion method at 298 K. An X-ray diffraction data set was collected to 2.8 A at 100 K. The crystal belonged to space group P2(1)2(1)2(1), with unit-cell parameters a = 34.52, b = 59.88, c = 72.93 A.

Insight into G[bond]T mismatch recognition using molecular dynamics with time-averaged restraints derived from NMR spectroscopy.

Molecular dynamics (MD) simulations were conducted for a G[bond]T mismatch-containing DNA decamer, d(CCATGCGTGG)(2), and its Watson-Crick parent sequence, d(CCACGCGTGG)(2). Dynamics in unrestrained MD trajectories were in poor agreement with prior (13)C NMR studies. However, the accuracy of the trajectories was improved by the use of time-averaged interatomic distance restraints derived from (1)H NMR. Postprocess smoothing of the trajectories further improved accuracy. Comparison of restrained and smoothed trajectories of the two DNA molecules revealed distinct differences in dynamics. The major groove width of the mismatched oligomer was more variable over the course of the simulation compared to its parent sequence. Greater variability in helical parameters stretch and opening for the mismatches indicated less kinetically stable base pairing. Interbase helical parameters rise, roll, and tilt were also more variable in certain base steps involving mismatched bases. These dynamic differences between normal and G[bond]T mismatched DNA reflect differences in local flexibility that may play a role in mismatch recognition by the MutS. A potential alternate G[bond]T mismatch binding mode for MutS is also proposed.

A molecular dynamics simulation study of oriented DNA with polyamine and sodium counterions: diffusion and averaged binding of water and cations.

Four different molecular dynamics (MD) simulations have been performed for ordered DNA decamers, d(5'-ATGCAGTCAG)*d(5'-TGACTGCATC). The counterions were the two natural polyamines spermidine3+ (Spd3+) and putrescine2+ (Put2+), the synthetic polyamine diaminopropane2+ (DAP2+) and Na+. The simulation set-up corresponds to an infinite array of parallel DNA mimicking the state in oriented DNA fibers or crystals. This work describes general properties of polyamine and Na+ binding to DNA. Simulated diffusion coefficients show satisfactory agreement with experimental NMR diffusion data of comparable systems. The interaction of the polyamines with DNA is dynamic in character and the cations mostly form short-lived contacts with the electronegative binding sites of DNA. Polyamines, Na+ and water interact most frequently with the charged phosphate atoms with preference for association from the minor groove side with O1P over O2P. There is a strong anti-correlation in the cation binding to the electronegative groups of DNA, i.e. the presence of a cation near one of the DNA sites repels other cations from binding to this and to the other sites separated by <7.5 A from each other. In contrast to the other polyamines, DAP2+ is able to form 'bridges' connecting neighboring phosphate groups along the DNA strand. A small fraction of DAP2+ and Put2+ can be found in the major grooves, while Spd3+ is absent there. The results of the MD simulations reveal principal differences in the polyamine-DNA interactions between the natural (Spd3+, Put2+ and spermine4+) and synthetic (DAP2+) polyamines.

NMR study on the interaction between RPA and DNA decamer containing cis-syn cyclobutane pyrimidine dimer in the presence of XPA: implication for damage verification and strand-specific dual incision in nucleotide excision repair.

In mammalian cells, nucleotide excision repair (NER) is the major pathway for the removal of bulky DNA adducts. Many of the key NER proteins are members of the XP family (XPA, XPB, etc.), which was named on the basis of its association with the disorder xerodoma pigmentosum. Human replication protein A (RPA), the ubiquitous single-stranded DNA-binding protein, is another of the essential proteins for NER. RPA stimulates the interaction of XPA with damaged DNA by forming an RPA-XPA complex on damaged DNA sites. Binding of RPA to the undamaged DNA strand is most important during NER, because XPA, which directs the excision nucleases XPG and XPF, must bind to the damaged strand. In this study, nuclear magnetic resonance (NMR) spectroscopy was used to assess the binding of the tandem high affinity DNA-binding domains, RPA-AB, and of the isolated domain RPA-A, to normal DNA and damaged DNA containing the cyclobutane pyrimidine dimer (CPD) lesion. Both RPA-A and RPA-AB were found to bind non- specifically to both strands of normal and CPD- containing DNA duplexes. There were no differences observed when binding to normal DNA duplex was examined in the presence of the minimal DNA-binding domain of XPA (XPA-MBD). However, there is a drastic difference for CPD-damaged DNA duplex as both RPA-A and RPA-AB bind specifically to the undamaged strand. The strand-specific binding of RPA and XPA to the damaged duplex DNA shows that RPA and XPA play crucial roles in damage verification and guiding cleavage of damaged DNA during NER.

Structural origins of adenine-tract bending.

DNA sequences containing short adenine tracts are intrinsically curved and play a role in transcriptional regulation. Despite many high-resolution NMR and x-ray studies, the origins of curvature remain disputed. Long-range restraints provided by 85 residual dipolar couplings were measured for a DNA decamer containing an adenine (A)(4)-tract and used to refine the structure. The overall bend in the molecule is a result of in-phase negative roll in the A-tract and positive roll at its 5' junction, as well as positive and negative tilt inside the A-tract and near its junctions. The bend magnitude and direction obtained from NMR structures is 9.0 degrees into the minor groove in a coordinate frame located at the third AT base pair. We evaluated long-range and wedge models for DNA curvature and concluded that our data for A-tract curvature are best explained by a "delocalized bend" model. The global bend magnitude and direction of the NMR structure are in excellent agreement with the junction model parameters used to rationalize gel electrophoretic data and with preliminary results of a cyclization kinetics assay from our laboratory.

New complex of post-activated neocarzinostatin chromophore with DNA: bulge DNA binding from the minor groove.

Neocarzinostatin (NCS-chrom), a natural enediyne antitumor antibiotic, undergoes either thiol-dependent or thiol-independent activation, resulting in distinctly different DNA cleavage patterns. Structures of two different post-activated NCS-chrom complexes with DNA have been reported, revealing strikingly different binding modes that can be directly related to the specificity of DNA chain cleavage caused by NCS-chrom. The third structure described herein is based on recent studies demonstrating that glutathione (GSH) activated NCS-chrom efficiently cleaves DNA at specific single-base sites in sequences containing a putative single-base bulge. In this structure, the GSH post-activated NCS-chrom (NCSi-glu) binds to a decamer DNA, d(GCCAGAGAGC), from the minor groove. This binding triggers a conformational switch in DNA from a loose duplex in the free form to a single-strand, tightly folded hairpin containing a bulge adenosine embedded between a three base pair stem. The naphthoate aromatic moiety of NCSi-glu intercalates into a GG step flanked by the bulge site, and its substituent groups, the 2-N-methylfucosamine carbohydrate ring and the tetrahydroindacene, form a complementary minor groove binding surface, mostly interacting with the GCC strand in the duplex stem of DNA. The bulge site is stabilized by the interactions involving NCSi-glu naphthoate and GSH tripeptide. The positioning of NCSi-glu is such that only single-chain cleavage via hydrogen abstraction at the 5'-position of the third base C (which is opposite to the putative bulge base) in GCC is possible, explaining the observed single-base cleavage specificity. The reported structure of the NCSi-glu-bulge DNA complex reveals a third binding mode of the antibiotic and represents a new family of minor groove bulge DNA recognition structures. We predict analogue structures of NCSi-R (R = glu or other substituent groups) may be versatile probes for detecting the existence of various structures of nucleic acids. The NMR structure of this complex, in combination with the previously reported NCSi-gb-bulge DNA complex, offers models for specific recognition of DNA bulges of various sizes through binding to either the minor or the major groove and for single-chain cleavage of bulge DNA sequences.