Alternative structures for alternating poly(dA-dT) tracts: the structure of the B-DNA decamer C-G-A-T-A-T-A-T-C-G

In-vitro and in-silico analysis and antitumor studies of novel Cu(II) and V(V) complexes of N-p-Tolylbenzohydroxamic acid

Copper(L2Cu) and vanadium(L2VOCl) complexes of N-p-tolylbenzohydroxamic acid (LH) ligand have been investigated for DNA binding efficacy by multiple analytical, spectral, and computational techniques. The results revealed that complexes as groove binders as evidenced by UV absorption. Fluorescence studies including displacement assay using classical intercalator ethidium bromide as fluorescent probe also confirmed as groove binders. The viscometric analysis too supports the inferences as strong groove binders for both the complexes. Molecular docking too exposed DNA as a target to the complexes which precisely binds L2Cu, in the minor groove region while L2VOCl in major groove region. Molecular dynamic simulation performed on L2Cu complex revealing the interaction of complex with DNA within 20 ns time. The complex stacked into the nitrogen bases of oligonucleotides and the bonding features were intrinsically preserved for longer simulation times. In-vitro cytotoxicity study was undertaken employing MTT assay against the breast cancer cell line (MCF-7). Potential cytotoxic activities were observed for L2Cu and L2VOCl complexes with IC50 values of showing 71 % and 74 % of inhibition respectively.

X-ray Structure Characterization of the Selective Recognition of AT Base Pair Sequences

The rational design of small molecules that target specific DNA sequences is a promising strategy to modulate gene expression. This report focuses on a diamidinobenzimidazole compound, whose selective binding to the minor groove of AT DNA sequences holds broad significance in the molecular recognition of AT-rich human promoter sequences. The objective of this study is to provide a more detailed and systematized understanding, at an atomic level, of the molecular recognition mechanism of different AT-specific sequences by a rationally designed minor groove binder. The specialized method of X-ray crystallography was utilized to investigate how the sequence-dependent recognition properties in general, A-tract, and alternating AT sequences affect the binding of diamidinobenzimidazole in the DNA minor groove. While general and A-tract AT sequences give a narrower minor groove, the alternating AT sequences intrinsically have a wider minor groove which typically constricts upon binding. A strong and direct hydrogen bond between the N-H of the benzimidazole and an H-bond acceptor atom in the minor groove is essential for DNA recognition in all sequences described. In addition, the diamidine compound specifically utilizes an interfacial water molecule for its DNA binding. DNA complexes of AATT and AAAAAA recognition sites show that the diamidine compound can bind in two possible orientations with a preference for water-assisted hydrogen bonding at either cationic end. The complex structures of AAATTT, ATAT, ATATAT, and AAAA are bound in a singular orientation. Analysis of the helical parameters shows a minor groove expansion of about 1 Å across all the nonalternating DNA complexes. The results from this systematic approach will convey a greater understanding of the specific recognition of a diverse array of AT-rich sequences by small molecules and more insight into the design of small molecules with enhanced specificity to AT and mixed DNA sequences.

Modulation of the helical properties of DNA: next-to-nearest neighbour effects and beyond

We used extensive molecular dynamics simulations to study the structural and dynamic properties of the central d(TpA) step in the highly polymorphic d(CpTpApG) tetranucleotide. Contrary to the assumption of the dinucleotide-model and its nearest neighbours (tetranucleotide-model), the properties of the central d(TpA) step change quite significantly dependent on the next-to-nearest (hexanucleotide) sequence context and in a few cases are modulated by even remote neighbours (beyond next-to-nearest from the central TpA). Our results highlight the existence of previously undescribed dynamical mechanisms for the transmission of structural information into the DNA and demonstrate the existence of certain sequences with special physical properties that can impact on the global DNA structure and dynamics.

6S RNA Mimics B-Form DNA to Regulate Escherichia coli RNA Polymerase

Noncoding RNAs (ncRNAs) regulate gene expression in all organisms. Bacterial 6S RNAs globally regulate transcription by binding RNA polymerase (RNAP) holoenzyme and competing with promoter DNA. Escherichia coli (Eco) 6S RNA interacts specifically with the housekeeping σ⁷⁰-holoenzyme (Eσ⁷⁰) and plays a key role in the transcriptional reprogramming upon shifts between exponential and stationary phase. Inhibition is relieved upon 6S RNA-templated RNA synthesis. We report here the 3.8 Å resolution structure of a complex between 6S RNA and Eσ⁷⁰ determined by single-particle cryo-electron microscopy and validation of the structure using footprinting and crosslinking approaches. Duplex RNA segments have A-form C3' endo sugar puckers but widened major groove widths, giving the RNA an overall architecture that mimics B-form promoter DNA. Our results help explain the specificity of Eco 6S RNA for Eσ⁷⁰ and show how an ncRNA can mimic B-form DNA to directly regulate transcription by the DNA-dependent RNAP.

Amphotericin B and anidulafungin directly interact with DNA and induce oxidative damage in the mammalian genome

Amphotericin B and anidulafungin are widely used antifungal drugs for the treatment of systemic and serious mycoses. Amphotericin B is a relatively toxic drug which has long been established. This study is first of its kind to systematically investigate the nature of binding to DNA, and to evaluate intercalation of AMP-B or ANIDULA with the aid of UV-Vis, ITC, and CD spectroscopy. The binding affinity of AMP-B with exclusion sites of 4.68 base pairs (1.2 × 10(5) M(-1)) was found to be higher than that of ANIDULA with exclusion sites of 6.67 base pairs (3.78 × 10(4) M(-1)); consistent with the binding affinity values obtained for AMP-B (10(5) M(-1)) and ANIDULA (10(4) M(-1)). The binding of two drugs with double-stranded DNA was favoured by negative enthalpy as well as negative entropy changes. The intercalation of drugs to duplex polynucleotide induced changes in the intrinsic CD spectra and revealed comparatively higher affinity towards AMP-B than ANIDULA. Molecular docking studies revealed that the negative binding energy was higher in the case of AMP-B reflecting more affinity towards single-stranded DNA. The results of the cytotoxicity, immunoblotting, and gene specific LA-QPCR assay have indicated that ANIDULA is less genotoxic than AMP-B. Hence, the superiority of ANIDULA over AMP-B as a systemic antifungal drug has been established beyond doubt.

The coumarin–pterocarpan conjugate – a natural product inspired hybrid molecular probe for DNA recognition

Thermally induced cascade sigmatropic rearrangement of a butynyloxycoumarin derivative has led to a quick access to the coumarin–pterocarpan hybrid molecule. Biophysical studies together with molecular modeling show that this nature-inspired hybrid molecule is capable of binding to the minor groove of DNA as a non-conventional entity.

Molecular basis for sequence-dependent induced DNA bending

With a growing understanding of the microstructural variations of DNA, it has become apparent that subtle conformational features are essential for specific DNA molecular recognition and function. DNA containing an A-tract has a narrow minor groove and a globally bent conformation but the structural features of alternating AT DNA are less well understood. Several studies indicate that alternating AT sequences are polymorphic with different global and local properties from A-tracts. The mobility of alternating AT DNA in gel electrophoresis is extensively reduced upon binding with minor-groove binding agents such as netropsin. Although this suggests that such complexes are bent, similarly to A-tract DNA, direct evidence and structural information on AT DNA and the induced conformational change is lacking. We have used NMR spectroscopy and residual dipolar coupling together with restrained molecular-dynamics simulations to determine the solution structures of an alternating AT DNA segment, with and without netropsin, in order to evaluate the molecular basis of the binding-induced effects. Complex formation causes a significant narrowing of the minor groove and a pronounced change in bending, from a slight bend towards the major groove for the free DNA to a pronounced bend towards the minor groove in the complex. This observation demonstrates that conformational features and the inherent malleability of AT sequences are essential for specific molecular recognition and function. These results take the field of DNA structures into new areas while opening up avenues to target novel DNA sequences.

Binding of Cu(II) complexes of oxicam NSAIDs to alternating AT and homopolymeric AT sequences: differential response to variation in backbone structure

Besides their principal functions as painkillers and anti-inflammatory agents, drugs belonging to the nonsteroidal anti-inflammatory drug (NSAID) group also have anticancer properties. Cu(II) complexes of these drugs enhance the anticancer effect. How they exert this effect is not clear. As a possible molecular mechanism, our group has already shown that the Cu(II) complexes of two oxicam NSAIDs with anticancer properties, viz. piroxicam and meloxicam, can directly bind to the DNA backbone. AT stretches are abundant in the eukaryotic genome. These stretches are more accessible to binding of different ligands, resulting in expression of different functions. AT stretches containing both alternating base pairs and homopolymeric bases in individual strands show subtle differences in backbone structures. It is therefore of interest to see how the Cu(II)-NSAID complexes respond to such differences in backbone structure. Binding studies of these complexes with alternating polydA-dT and homopolymeric polydA-polydT have been conducted using UV-vis absorption titration studies, UV melting studies and circular dichroism spectroscopy. Competitive binding with the standard intercalator ethidium bromide and the minor groove binder 4',6-diamidino-2-phenylindole was monitored using fluorescence to identify the possible binding mode. Our results show that Cu(II)-NSAID complexes are highly sensitive to the subtle differences in backbone structures of polydA-dT and polydA-polydT and respond to them by exhibiting different binding properties, such as binding constants, effect on duplex stability and binding modes. Both complexes have a similar binding mode with polydA-dT, which is intercalative, but for polydA-polydT, the results point to a mixed mode of binding.

Nucleic acid-metal ion interactions in the solid state

Metal ions play a key role in nucleic acid structure and activity. Elucidation of the rules that govern the binding of metal ions is therefore an essential step for better understanding of the nucleic acid functions. This review is as an update to a preceding one (Metal Ions Biol. Syst., 1996, 32, 91-134), in which we offered a general view of metal ion interactions with mono-, di-, tri-, and oligonucleotides in the solid state, based on their crystal structures reported before 1994. In this chapter, we survey all the crystal structures of metal ion complexes with nucleotides involving oligonucleotides reported after 1994 and we have tried to uncover new characteristic metal bonding patterns for mononucleotides and oligonucleotides with A-RNA and A/B/Z-DNA fragments that form duplexes. We do not cover quadruplexes, duplexes with metal-mediated base-pairs, tRNAs, rRNAs in ribosome, ribozymes, and nucleic acid-drug and -protein complexes. Factors that affect metal binding to mononucleotides and oligonucleotide duplexes are also dealt with.

Crystal structure of the chromodomain helicase DNA-binding protein 1 (Chd1) DNA-binding domain in complex with DNA

Chromatin remodelers are ATP-dependent machines that dynamically alter the chromatin packaging of eukaryotic genomes by assembling, sliding, and displacing nucleosomes. The Chd1 chromatin remodeler possesses a C-terminal DNA-binding domain that is required for efficient nucleosome sliding and believed to be essential for sensing the length of DNA flanking the nucleosome core. The structure of the Chd1 DNA-binding domain was recently shown to consist of a SANT and SLIDE domain, analogous to the DNA-binding domain of the ISWI family, yet the details of how Chd1 recognized DNA were not known. Here we present the crystal structure of the Saccharomyces cerevisiae Chd1 DNA-binding domain in complex with a DNA duplex. The bound DNA duplex is straight, consistent with the preference exhibited by the Chd1 DNA-binding domain for extranucleosomal DNA. Comparison of this structure with the recently solved ISW1a DNA-binding domain bound to DNA reveals that DNA lays across each protein at a distinct angle, yet contacts similar surfaces on the SANT and SLIDE domains. In contrast to the minor groove binding seen for Isw1 and predicted for Chd1, the SLIDE domain of the Chd1 DNA-binding domain contacts the DNA major groove. The majority of direct contacts with the phosphate backbone occur only on one DNA strand, suggesting that Chd1 may not strongly discriminate between major and minor grooves.

Coiled-coil conformation of a pentamidine–DNA complex

The coiled-coil structure formed by the complex of the DNA duplex d(ATATATATAT)2with pentamidine is presented. The duplex was found to have a mixed structure containing Watson–Crick and Hoogsteen base pairs. The drug stabilizes the coiled coil through the formation of cross-links between neighbouring duplexes. The central part of the drug is found in the minor groove as expected, whereas the charged terminal amidine groups protrude and interact with phosphates from neighbouring molecules. The formation of cross-links may be related to the biological effects of pentamidine, which is used as an antiprotozoal agent in trypanosomiasis, leishmaniasis and pneumonias associated with AIDS. The DNA sequence that was used is highly abundant in most eukaryotic genomes. However, very few data are available on DNA sequences which only contain A·T base pairs.

Conformations of the sugar-phosphate backbone in helical DNA crystal structures

The DNA backbone geometry was analyzed for 96 crystal structures of oligodeoxynucleotides. The ranges and mean values of the torsion angles for the observed subclasses of the A-, B-, and Z-DNA conformational types were determined by analyzing distributions of the torsion angles and scattergrams relating pairs of angles.

DNA polymorphism as an origin of adenine-thymine tract length-dependent threading intercalation rate

Binuclear ruthenium complexes that bind DNA by threading intercalation have recently been found to exhibit an exceptional kinetic selectivity for long polymeric adenine-thymine (AT) DNA. A series of oligonucleotide hairpin duplexes containing a central tract of 6-44 alternating AT base pairs have here been used to investigate the nature of the recognition mechanism. We find that, above a threshold AT tract length corresponding to one helix turn of B-DNA, a dramatic increase in threading intercalation rate occurs. In contrast, such length dependence is not observed for rates of unthreading. Intercalation by any mechanism that depends on the open end of the hairpin was found not to be important in the series of oligonucleotides used, as verified by including in the study a hairpin duplex cyclized by a copper-catalyzed "click" reaction. Our observations are interpreted in terms of a conformational pre-equilibrium, determined by the length of the AT tract. We finally find that mismatches or loops in the oligonucleotide facilitate the threading process, of interest for the development of mismatch-recognizing probes.

Strain softening in stretched DNA

The microscopic mechanics of DNA stretching was characterized using extensive molecular dynamics simulations. By employing an anisotropic pressure-control method, realistic force-extension dependences of effectively infinite DNA molecules were obtained. A coexistence of B and S DNA domains was observed during the overstretching transition. The simulations revealed that strain softening may occur in the process of stretching torsionally constrained DNA. The latter observation was qualitatively reconciled with available experimental data using a random-field Ising model.

C-->G base mutations in the CArG box of c-fos serum response element alter its bending flexibility. Consequences for core-SRF recognition

By binding to the CArG box sequence, the serum response factor (SRF) activates several muscle-specific genes, as well as genes that respond to mitogens. The core domain of the SRF (core-SRF) binds as a dimer to the CArG box C-5C-4A-3T-2A-1T+1T+2A+3G+4G+5 of the c-fos serum response element (SREfos). However, previous studies using 20-mer DNAs have shown that the binding stoichiometry of core-SRF is significantly altered by mutations C-5-->G (SREGfos) and C-5C-4-->GG (SREGGfos) of the CArG box [A Huet, A Parlakian, M-C Arnaud, J-M Glandières, P Valat, S Fermandjian, D Paulin, B Alpert & C Zentz (2005) FEBS J272, 3105-3119]. To understand these effects, we carried out a comparative analysis of the three 20-mer DNAs SREfos, SREGfos and SREGGfos in aqueous solution. Their CD spectra were of the B-DNA type with small differences generated by variations in the mutual arrangement of the base pairs. Analysis by singular value decomposition of a set of Raman spectra recorded as a function of temperature, revealed a premelting transition associated with a conformational shift in the DNA double helices from a bent to a linear form. Time-resolved fluorescence anisotropy shows that the fluorescein reporter linked to the oligonucleotide 5'-ends experiences twisting motions of the double helices related to the interconversion between bent and linear conformers. The three SREs present various bent populations submitted, however, to particular internal dynamics, decisive for the mutual adjustment of binding partners and therefore specific complex formation.

Minor groove deformability of DNA: a molecular dynamics free energy simulation study

The conformational deformability of nucleic acids can influence their function and recognition by proteins. A class of DNA binding proteins including the TATA box binding protein binds to the DNA minor groove, resulting in an opening of the minor groove and DNA bending toward the major groove. Explicit solvent molecular dynamics simulations in combination with the umbrella sampling approach have been performed to investigate the molecular mechanism of DNA minor groove deformations and the indirect energetic contribution to protein binding. As a reaction coordinate, the distance between backbone segments on opposite strands was used. The resulting deformed structures showed close agreement with experimental DNA structures in complex with minor groove-binding proteins. The calculated free energy of minor groove deformation was approximately 4-6 kcal mol(-1) in the case of a central TATATA sequence. A smaller equilibrium minor groove width and more restricted minor groove mobility was found for the central AAATTT and also a significantly ( approximately 2 times) larger free energy change for opening the minor groove. The helical parameter analysis of trajectories indicates that an easier partial unstacking of a central TA versus AT basepair step is a likely reason for the larger groove flexibility of the central TATATA case.

Effect of DNA structural flexibility on promoter strength--molecular dynamics studies of E. coli promoter sequences

To study possible correlations between promoter activity and the structural flexibility of the DNA helix, we have carried out unrestrained molecular dynamics simulations of the -10 consensus region sequence and five variants forming the -10 region of various Escherichia coli promoter sequences. Analyses of the trajectories obtained from the simulations show that the consensus sequence has a pattern of two structurally flexible nucleotide steps sandwiched between two stiff steps. In the other sequences, this pattern varies in consonance with the change in the sequence. The variations in the patterns show correlation with the promoter strength.

Predicting indirect readout effects in protein-DNA interactions

Recognition of DNA by proteins relies on direct interactions with specific DNA-functional groups, along with indirect effects that reflect variable energetics in the response of DNA sequences to twisting and bending distortions induced by proteins. Predicting indirect readout requires knowledge of the variations in DNA curvature and flexibility in the affected region, which we have determined for a series of DNA-binding sites for the E2 regulatory protein by using the cyclization kinetics method. We examined 16 sites containing different noncontacted spacer sequences, which vary by more than three orders of magnitude in binding affinity. For 15 of these sites, the variation in affinity was predicted within a factor of 3, by using experimental curvature and flexibility values and a statistical mechanical theory. The sole exception was traced to differential magnesium ion binding.

Crystal structure of an antiparallel DNA fragment with Hoogsteen base pairing

We report here an alternative double-helical structure of the DNA molecule. It has been found in the d(ATA(Br)UAT) and d(ATATAT) sequences by single-crystal x-ray crystallography. This sequence is found not only in TATA boxes, but also in other regulatory regions of DNA. Bases of the two antiparallel strands form Hoogsteen pairs, with adenines in the syn conformation. The structure is related neither to those found in triple helices nor to parallel DNA duplexes. Its conformational parameters are very similar to those of duplex DNA in the B form. Both forms may coexist under physiological conditions, although the Hoogsteen pairing greatly influences the recognition sites on DNA. Our results demonstrate that an alternative to the classical B-DNA double helix is possible.

Intrinsic bending and deformability at the T-A step of CCTTTAAAGG: a comparative analysis of T-A and A-T steps within A-tracts

Introduction of a T-A or pyrimidine-purine step into a straight and rigid A-tract can cause a positive roll deformation that kinks the DNA helix at that step. In CCTTTAAAGG, the central T-A step has an 8.6 degrees bend toward the major groove. We report the structural analysis of CCTTTAAAGG and a comparison with 25 other representative crystal structures from the NDB containing at least four consecutive A or T bases. On average, more local bending occurs at the disruptive T-A step (8.21 degrees ) than at an A-T step (5.71 degrees ). In addition, A-tracts containing an A-T step are more bent than are pure A-tracts, and hence A-A and A-T steps are not equivalent. All T-A steps examined exhibit positive roll, bending towards the major groove, while A-T steps display negative roll and bend slightly towards the minor groove. This illustrates how inherent negative and positive roll are, respectively, at A-T and T-A steps within A-tracts. T-A steps are more deformable, showing larger and more variable deformations of minor groove width, rise, cup, twist, and buckle. Standard deviations of twist, rise, and cup for T-A steps are 6.66 degrees, 0.55 A, and 15.90 degrees, versus 2.28 degrees, 0.21 A, and 2.99 degrees for A-T steps. Packing constraints determine which local values of these helical parameters an individual T-A step will adopt. For instance, with CCTTTAAAGG and three isomorphous structures, CGATTAATCG, CGATATATCG, and CGATCGATCG, crystal packing forces lead to a series of correlated changes: widened minor groove, large slide, low twist, and large rise. The difference in helical parameters between A-T steps lying within A-tracts, versus A-T steps within alternating AT sequences, demonstrates the importance of neighboring steps on the conformation of a given dinucleotide step.

1 A crystal structures of B-DNA reveal sequence-specific binding and groove-specific bending of DNA by magnesium and calcium

The 1 A resolution X-ray crystal structures of Mg(2+) and Ca(2+) salts of the B-DNA decamers CCAACGTTGG and CCAGCGCTGG reveal sequence-specific binding of Mg(2+) and Ca(2+) to the major and minor grooves of DNA, as well as non-specific binding to backbone phosphate oxygen atoms. Minor groove binding involves H-bond interactions between cross-strand DNA base atoms of adjacent base-pairs and the cations' water ligands. In the major groove the cations' water ligands can interact through H-bonds with O and N atoms from either one base or adjacent bases, and in addition the softer Ca(2+) can form polar covalent bonds bridging adjacent N7 and O6 atoms at GG bases. For reasons outlined earlier, localized monovalent cations are neither expected nor found.Ultra-high atomic resolution gives an unprecedented view of hydration in both grooves of DNA, permits an analysis of individual anisotropic displacement parameters, and reveals up to 22 divalent cations per DNA duplex. Each DNA helix is quite anisotropic, and alternate conformations, with motion in the direction of opening and closing the minor groove, are observed for the sugar-phosphate backbone. Taking into consideration the variability of experimental parameters and crystal packing environments among these four helices, and 24 other Mg(2+) and Ca(2+) bound B-DNA structures, we conclude that sequence-specific and strand-specific binding of Mg(2+) and Ca(2+) to the major groove causes DNA bending by base-roll compression towards the major groove, while sequence-specific binding of Mg(2+) and Ca(2+) in the minor groove has a negligible effect on helix curvature. The minor groove opens and closes to accommodate Mg(2+) and Ca(2+) without the necessity for significant bending of the overall helix. The program Shelxdna was written to facilitate refinement and analysis of X-ray crystal structures by Shelxl-97 and to plot and analyze one or more Curves and Freehelix output files.

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.

Crystal studies of B-DNA: the answers and the questions

The history of DNA crystallography is reviewed and is followed by discussion of the methods used for structure determination. The features of B-DNA molecular and crystal structures are described.

Signals for TBP/TATA box recognition

The TATA box-binding protein (TBP) recognizes its target sites (TATA boxes) by indirectly reading the DNA sequence through its conformation effects (indirect readout). Here, we explore the molecular mechanisms underlying indirect readout of TATA boxes by TBP by studying the binding of TBP to adenovirus major late promoter (AdMLP) sequence variants, including alterations inside as well as in the sequences flanking the TATA box. We measure here the dissociation kinetics of complexes of TBP with AdMLP targets and, by phase-sensitive assay, the intrinsic bending in the TATA box sequences as well as the bending of the same sequence induced by TBP binding. In these experiments we observe a correlation of the kinetic stability to sequence changes within the TATA recognition elements. Comparison of the kinetic data with structural properties of TATA boxes in known crystalline TBP/TATA box complexes reveals several "signals" for TATA box recognition, which are both on the single base-pair level, as well as larger DNA tracts within the TATA recognition element. The DNA bending induced by TBP on its binding sites is not correlated to the stability of TBP/TATA box complexes. Moreover, we observe a significant influence on the kinetic stability of alteration in the region flanking the TATA box. This effect is limited however to target sites with alternating TA sequences, whereas the AdMLP target, containing an A tract, is not influenced by these changes.

Base opening in RNA and DNA duplexes: implication for RNA stability

The energetics of a low-energy single base opening in several RNA duplex crystal structures has been calculated and compared to DNA duplexes. Base opening in RNA appears to have an overall preference towards the major groove, similar to results previously reported for B-DNA. Movement of each of the adenine, uracil, and cytosine bases into the minor groove is blocked by a high-energy barrier due to severe close contact with neighboring bases. Guanine bases are able to open towards both grooves because of the unique orientation of the base that avoids steric clash along the opening pathway. RNA bases are found to have a substantially smaller major groove opening extent than that of their B-DNA counterparts. A comparison with base opening behavior of A-DNA duplexes suggests that this difference results from helix constraint associated with A-form backbone conformation. The reduced opening extent correlates with the RNA duplex stability and is consistent with observed slower imino proton exchange rates in RNA duplexes.

Conformational flexibility of B-DNA at 0.74 A resolution: d(CCAGTACTGG)(2)

The affinity and specificity of a ligand for its DNA site is a function of the conformational changes between the isolated and complexed states. Although the structures of a hydroxypyrrole-imidazole-pyrrole polyamide dimer with 5'-CCAGTACTGG-3' and the trp repressor recognizing the sequence 5'-GTACT-3' are known, the baseline conformation of the DNA site would contribute to our understanding of DNA recognition by these ligands. The 0.74 A resolution structure of a B-DNA double helix, 5'-CCAGTACTGG-3', has been determined by X-ray crystallography. Six of the nine phosphates, two of four bound calcium ions and networks of water molecules hydrating the oligonucleotide have alternate conformations. By contrast, nine of the ten bases have a single, unique conformation with hydrogen atoms visible in most cases. The polyamide molecules alter the geometry of the phosphodiester backbone, and the water molecules mediating contacts in the trp repressor/operator complex are conserved in the unliganded DNA. Furthermore, the multiple conformational states, ions and hydration revealed by this ultrahigh resolution structure of a B-form oligonucleotide are potentially general considerations for understanding DNA-binding affinity and specificity by ligands.

The impact of abasic sites on DNA flexibility

We use internal coordinate molecular mechanics calculations to study the impact of abasic sites on the conformation and the mechanics of the DNA double helix. Abasic sites, which are common mutagenic lesions, are shown to locally modify both the groove geometry and the curvature of DNA in a sequence dependent manner. By controlled twisting and bending, it is also shown that these lesions modify the deformability of the duplex, generally increasing its flexibility, but again to an extent which depends on the nature of the abasic site and on the surrounding base sequence. Both the conformational and mechanical influence of this type of DNA damage may be significant for recognition and repair mechanisms.

Structural effects of DNA sequence on T.A recognition by hydroxypyrrole/pyrrole pairs in the minor groove

Synthetic polyamides composed of three types of aromatic amino acids, N-methylimidazole (Im), N-methylpyrrole (Py) and N-methyl-3-hydroxypyrrole (Hp) bind specific DNA sequences as antiparallel dimers in the minor groove. The side-by-side pairings of aromatic rings in the dimer afford a general recognition code that allows all four base-pairs to be distinguished. To examine the structural consequences of changing the DNA sequence context on T.A recognition by Hp/Py pairs in the minor groove, crystal structures of polyamide dimers (ImPyHpPy)(2) and the pyrrole counterpart (ImPyPyPy)(2) bound to the six base-pair target site 5'-AGATCT-3' in a ten base-pair oligonucleotide have been determined to a resolution of 2.27 and 2.15 A, respectively. The structures demonstrate that the principles of Hp/Py recognition of T.A are consistent between different sequence contexts. However, a general structural explanation for the non-additive reduction in binding affinity due to introduction of the hydroxyl group is less clear. Comparison with other polyamide-DNA cocrystal structures reveals structural themes and differences that may relate to sequence preference.

A structural snapshot of base-pair opening in DNA

The response of double-helical DNA to torsional stress may be a driving force for many processes acting on DNA. The 1.55-A crystal structure of a duplex DNA oligonucleotide d(CCAGGCCTGG)(2) with an engineered crosslink in the minor groove between the central guanine bases depicts how the duplex can accommodate such torsional stress. We have captured in the same crystal two rather different conformational states. One duplex contains a strained crosslink that is stabilized by calcium ion binding in the major groove, directly opposite the crosslink. For the other duplex, the strain in the crosslink is relieved through partial rupture of a base pair and partial extrusion of a cytosine accompanied by helix bending. The sequence used is the target sequence for the HaeIII methylase, and this partially flipped cytosine is the same nucleotide targeted for extrusion by the enzyme. Molecular dynamics simulations of these structures show an increased mobility for the partially flipped-out cytosine.

Water and ions in a high resolution structure of B-DNA

A detailed picture of hydration and counterion location in the B-DNA duplex d(GCGAATTCG) is presented. Detailed data have been obtained by single crystal x-ray diffraction at atomic resolution (0.89 A) in the presence of Mg(2+). The latter is the highest resolution ever obtained for a B-DNA oligonucleotide. Minor groove hydration is compared with that found in the Na(+) and Ca(2+) crystal forms of the related dodecamer d(CGCGAATTCGCG). High resolution data (1.45 A) of the Ca(2+) form obtained in our laboratory are used for that purpose. The central GAATTC has a very stable hydration spine identical in all cases, independent of duplex length and crystallization conditions (counterions, space group). However, the organization of the water molecules (tertiary and quaternary layers) associated with the central spine vary in each case.

DNA minor groove recognition of a non-self-complementary AT-rich sequence by a tris-benzimidazole ligand

The crystal structure of the non-self-complementary dodecamer DNA duplex formed by d(CG[5BrC]ATAT-TTGCG) and d(CGCAAATATGCG) has been solved to 2.3 A resolution, together with that of its complex with the tris-benzimidazole minor groove binding ligand TRIBIZ. The inclusion of a bromine atom on one strand in each structure enabled the possibility of disorder to be discounted. The native structure has an exceptional narrow minor groove, of 2.5-2.6 A in the central part of the A/T region, which is increased in width by approximately 0.8 A on drug binding. The ligand molecule binds in the central part of the sequence. The benzimidazole subunits of the ligand participate in six bifurcated hydrogen bonds with A:T base pair edges, three to each DNA strand. The presence of a pair of C-H...O hydrogen bonds has been deduced from the close proximity of the pyrrolidine group of the ligand to the TpA step in the sequence.

Modeling high-resolution hydration patterns in correlation with DNA sequence and conformation

Hydration around the DNA fragment d(C5T5).(A5G5) is presented from two molecular dynamics simulations of 10 and 12 ns total simulation time. The DNA has been simulated as a flexible molecule with both the CHARMM and AMBER force fields in explicit solvent including counterions and 0.8 M additional NaCl salt. From the previous analysis of the DNA structure B-DNA conformations were found with the AMBER force-field and A-DNA conformations with CHARMM parameters. High-resolution hydration patterns are compared between the two conformations and between C.G and T.A base-pairs from the homopolymeric parts of the simulated sequence. Crystallographic results from a statistical analysis of hydration sites around DNA crystal structures compare very well with the simulation results. Differences between the crystal sites and our data are explained by variations in conformation, sequence, and limitations in the resolution of water sites by crystal diffraction. Hydration layers are defined from radial distribution functions and compared with experimental results. Excellent agreement is found when the measured experimental quantities are compared with the equivalent distribution of water molecules in the first hydration shell. The number of water molecules bound to DNA was found smaller around T.A base-pairs and around A-DNA as compared to B-DNA. This is partially offset by a larger number of water molecules in hydrophobic contact with DNA around T.A base-pairs and around A-DNA. The numbers of water molecules in minor and major grooves have been correlated with helical roll, twist, and inclination angles. The data more fully explain the observed B-->A transition at low humidity.

TIT for TAT: the properties of inosine and adenosine in TATA box DNA

The sequence dependent conformation, flexibility and hydration properties of DNA molecules constitute selectivity determinants in the formation of protein-DNA complexes. TATA boxes in which AT basepairs (bp) have been substituted by IC bp (TITI box) allow for probing these selectivity determinants for the complexation with the TATA box-binding protein (TBP) with different sequences but identical chemical surfaces. The reference promoter Adenovirus 2 Major Late Promoter (mlp) is formed by the apposition of two sequences with very different dynamic properties: an alternating TATA sequence and an A-tract. For a comparative study, we carried out molecular dynamics simulations of two DNA oligomers, one containing the mlp sequence (2 ns), and the other an analog where AT basepairs were substituted by IC basepairs (1 ns). The simulations, carried out with explicit solvent and counterinons, yield straight purine tracts, the A-tract being stiffer than the I-tract, an alternating structure for the YRYR tracts, and hydration patterns that differ between the purine tracts and the alternating sequence tracts. A detailed analysis of the proposed interactions responsible for the stiffness of the purine tracts indicates that the stacking between the bases bears the strongest correlation to stiffness. The hydration properties of the minor groove in the two oligomers are distinctly different. Such differences are likely to be responsible for the stronger binding of TBP to mlp over the inosine-substituted variant. The calculations were made possible by the development, described here, of a new set of forcefield parameters for inosine that complement the published CHARMM all-hydrogen nucleic acid parametrization.

The structure of the most studied DNA fragment changes under the influence of ions: a new packing of d(CGCGAATTCGCG)

The title oligonucleotide and many related dodecamers have been extensively studied alone and as DNA-drug complexes. In practically all cases they were found to crystallize in the same space group, stabilized by interactions among the terminal guanine bases. Here we report new packing interactions (R3) in the presence of Ca2+. The oligonucleotides interact by placing their terminal guanines in the narrow groove of a neighbor molecule, an interaction which had never been found in dodecamers.

Hydration of the phosphate group in double-helical DNA

Water distributions around phosphate groups in 59 B-, A-, and Z-DNA crystal structures were analyzed. It is shown that the waters are concentrated in six hydration sites per phosphate and that the positions and occupancies of these sites are dependent on the conformation and type of nucleotide. The patterns of hydration that are characteristic of the backbone of the three DNA helical types can be attributed in part to the interactions of these hydration sites.

DNA-drug interaction measurements using surface plasmon resonance

The interactions of the drugs 2,7-bis[(diethylamino)-ethoxy]-fluoren-9-one dihydrochloride (Tilorone), 2,7-bis[(dipropylamino)-acetamido]-fluoren-9-one dihydrochloride (FA-2), 2'-(4-hydroxyphenyl)-5-(4-methyl-1-piperazinyl)-2,5'-bi-1H-benzimidazole trihydrochloride (Hoechst 33258), and hematoporphyrin IX derivative (HPD) with synthetic self-complementary DNA (36-b.p.; 5'-biotin-spacer-[d(CGCTATATAGCG)]3-3') were studied by SPR (Surface Plasmon Resonance). Monolayers of biotinylated DNA were immobilized on a streptavidin-dextran-gold triple-layer. Small portions of the drugs (approximately 30 pmol/ml) were injected in continuous flow. The mass corresponded to the amount of the bound molecules. Injections of 50 mM sodium hydroxide pulses separated the DNA double strands, releasing the effector molecules. Subsequent treatments with the effectors gave reproducible results. The maximum interaction between drug and DNA was observed in the case of Tilorone. 41 molecules could bind to the 36-b.p. DNA duplex. To investigate the microscopic behavior in condensed nucleic acid phases, SFM (Scanning Force Microscopy)-imaging and polarizing microscopic observations of DNA-effector complexes were carried out. Supplementary UV-absorption thermal denaturation curves of DNA with the above-mentioned effectors in dilute solutions were measured. As an additional aid to understand the geometries of DNA-drug interactions, computer simulations were performed and compared with the experimental data.

On the potential role of the amino nitrogen atom as a hydrogen bond acceptor in macromolecules

Crystallographic studies of duplex DNA have indicated that opposing exocyclic amino groups may form close NH⋯:N contacts. To study the nature of such interactions, we have examined the database of small molecule, high-resolution crystal structures for more accurate examples of this type of unconventional interaction. We have found cases where the amino groups in guanine and adenine bases accept hydrogen bonds from conventional donors, such as amino or hydroxyl groups. More frequently, the purine amino group was found to contact closely electropositive C-H groups. Searches of the nucleic acid structural databases also yielded several examples where the purine amino group is contacted by hydrogen bond donors in macromolecules. Ab initio calculations indicate that the hydrogen-amino contact is improved energetically when the amino group moves from the conventional geometry, where all atoms are co-planar with the base, to one in which the hydrogen atoms lie out of the plane and the nitrogen is at the apex of a pyramid, resulting in polarization of the amino group. The combined structural and theoretical data suggest that the amino group is flexible, and can accommodate close contacts, because the resulting polarization permits electropositive atoms to approach the amino group nitrogen more closely than expected for their conventional van der Waals radii. The flexibility of the amino group may permit particular DNA conformations that enforce hydrogen-amino contacts to optimize favorable stacking interactions, and it may play a role in the recognition of nucleosides. We speculate that the amino group can accept hydrogen bonds under special circumstances in macromolecules, and that this ability might play a mechanistic role in catalytic processes such as deamination or amino transfer.

Interresidue quiet NOEs for DNA structural studies

The potential utility of long-range NOEs in DNA has not been exploited since the observed signals have contributions both from the direct magnetization route and from multiple diffusion pathways. The Quiet NOE approach can be used to select for the direct magnetization transfer pathway by suppressing spin diffusion. A single-band Quiet NOE, which allows detection of the direct NOEs between protons in a selected chemical shift window, has been demonstrated on two duplex DNAs, and the NOEs observed can contain important structural information.

Conformational properties of DNA strands containing guanine-adenine and thymine-adenine repeats

CD spectroscopy and PAGE were used to cooperatively analyze melting conformers of DNA strands containing GA and TA dinucleotide repeats. The 20mer (GA)10 formed a homoduplex in neutral solutions containing physiological concentrations of salts and this homoduplex was not destabilized even in the terminal (GA)3 hexamers of (GA)3(TA)4(GA)3, although the central (TA)4 portion of this oligonucleotide preserved the conformation adopted by (TA)10. This observation demonstrates that homoduplexes of alternating GA and TA sequences can co-exist in a single DNA molecule. Another 20mer, (GATA)5, adopted as a whole either the AT duplex, like (TA)10, or the GA duplex, like (GA)10, and switched between them reversibly. The concentration of salt controlled the conformational switching. Hence, guanine and thymine share significant properties regarding complementarity to adenine, while the TA and GA sequences can stack in at least two mutually compatible ways within the DNA duplexes analyzed here. These properties extend our knowledge of non-canonical structures of DNA.

Metal-nucleotide interactions: crystal structures of alkali (Li+, Na+, K+) and alkaline earth (Ca2+, Mg2+) metal complexes of adenosine 2'-monophosphate

Crystal structures of lithium, sodium, potassium, calcium and magnesium salts of adenosine 2'-monophosphate (2'-AMP) have been obtained at atomic resolution by X-ray crystallographic methods. 2'-AMP.Li belongs to the monoclinic space group P2(1) with a = 7.472(3)A, b = 26.853(6) A, c = 9.184(1)A, b = 113.36(1)A and Z= 4. 2'-AMP.Na and 2'-AMP.K crystallize in the trigonal space groups P3(1) and P3(1)21 with a = 8.762(1)A, c = 34.630(5)A, Z= 6 and a = 8.931(4), Ac = 34.852(9)A and Z= 6 respectively while 2'-AMP.Ca and 2'-AMP.Mg belong to space groups P6(5)22 and P2(1) with cell parameters a = 9.487(2), c = 74.622(13), Z = 12 and a = 4.973(1), b = 10.023(2), c = 16.506(2), beta = 91.1(0) and Z = 2 respectively. All the structures were solved by direct methods and refined by full matrix least-squares to final R factors of 0.033, 0.028, 0.075, 0.069 and 0.030 for 2'-AMP.Li, 2'-AMP.Na, 2'- AMP.K, 2'-AMP.Ca and 2'-AMP.Mg, respectively. The neutral adenine bases in all the structures are in syn conformation stabilized by the O5'-N3 intramolecular hydrogen bond as in free acid and ammonium complex reported earlier. In striking contrast, the adenine base is in the anti geometry (chiCN = -156.4(2)degrees) in 2'-AMP.Mg. Ribose moieties adopt C2'-endo puckering in 2'-AMP.Li and 2'-AMP.Ca, C2'-endo-C3'-exo twist puckering in 2'-AMP.Na and 2'-AMP.K and a C3'-endo-C2'-exo twist puckering in 2'-AMP.Mg structure. The conformation about the exocyclic C4'-C5' bond is the commonly observed gauche-gauche (g+) in all the structures except the gauche- trans (g-) conformation observed in 2'-AMP.Mg structure. Lithium ions coordinate with water, ribose and phosphate oxygens at distances 1.88 to 1.99A. Na+ ions and K+ ions interact with phosphate and ribose oxygens directly and with N7 indirectly through a water oxygen. A distinct feature of 2'-AMP.Na and 2'-AMP.K structures is the involvement of ribose 04' in metal coordination. The calcium ion situated on a two-fold axis coordinates directly with three oxygens OW1, OW2 and O2 and their symmetry mates at distances 2.18 to 2.42A forming an octahedron. A classic example of an exception to the existence of the O5'-N3 intramolecular hydorgen bond is the 2'-AMP.Mg strucure. Magnesium ion forms an octahedral coordination with three water and three phosphate oxygens at distances ranging from 2.02 to 2.11 A. A noteworthy feature of its coordination is the indirect link with N3 through OW3 oxygen resulting in macrochelation between the base and the phosphate group. Greater affnity of metal clays towards 5' compared to 2' and 3' nucleotides (J. Lawless, E. Edelson, and L. Manning, Am. Chem. Soc. Northwest Region Meeting, Seattle. 1978) due to macrochelation infered from solution studies (S. S. Massoud, H. Sigel, Eur J. Biochem. 179, 451-458 (1989)) and interligand hydrogen bonding induced by metals postulated from metal-nucleotide structures in solid state (V. Swaminathan and M. Sundaralingam, CRC. Crit. Rev. Biochem. 6, 245-336 (1979)) are borne out by our structures also. The stacking patterns of adenine bases of both 2'-AMP.Na and 2'-AMP.K structures resemble the 2'-AMP.NH4 structure reported in the previous article. 2'-AMP.Li, 2'-AMP.Ca and 2'-AMP.Mg structures display base-ribose 04' stacking. An overview of interaction of monovalent and divalent cations with 2' and 5'-nucleotides has been presented.

A crystallographic approach to DNA bending: prediction of nucleosome formation by DNA triple repeats and other repetitive sequences

DNA bending is due to two main factors: the inherent curvature of the sequence and its flexibility. Most methods of analysis do not allow a differentiation between these two factors. In this paper I show that the flexibility of DNA sequences can be estimated from the standard deviation of roll values determined by X-ray crystallography for each base step. As an application of this approach, the nucleosome formation ability of triple repeat sequences has been determined and shown to be in agreement with the experimental results. Local variations in twist do not appear to have any influence on nucleosome formation.

Sensitivity of NMR internucleotide distances to B-DNA conformation: underlying mechanics

Nuclear magnetic resonance (NMR) spectroscopy, combining correlated spectroscopy (COSY) coupling constant measurements with nuclear Overhauser effect spectroscopy (NOESY) interatomic distances, should make it possible to determine an averaged solution structure for DNA oligomers. However, even if such data could be obtained with high accuracy, it is not clear which structural parameters of DNA would be determined. Here, the relationships between measurable internucleotide distances and helical parameters are systematically studied through molecular modelling. Investigations are carried out using four representative sequences, (ACGT)n, (TCGA)n, (AGCT)n and (TGCA)n, composed of repeated tetranucleotides belonging to oligomers previously studied by NMR. Correlations between interatomic distances become evident and strong connections between distances and inter-base helical parameters are observed. Results imply that twist, roll, shift and slide values can be accurately determined from NMR data. Sequence independent mechanical coupling which link backbone and sugar conformations to helical twist are also described.

Structure of a DNA analog of the primer for HIV-1 RT second strand synthesis

The non-self-complementary DNA decamer C-A-A-A-G-A-A-A-A-G/C-T-T-T-T-C-T-T-T-G is a DNA/DNA analogue of a portion of the polypurine tract or PPT, which is a RNA/DNA hybrid that serves as a primer for synthesis of the (+) DNA strand by HIV reverse transcriptase (RT), and which is not digested by the RNase H domain of reverse transcriptase following (-) strand synthesis. The same unusual conformation that eludes RNase H, thought to be a change in width of minor groove, may also be responsible for the inhibition of HIV RT by minor groove binding drugs such as distamycin and their bis-linked derivatives. The present X-ray crystal structure of this DNA decamer exhibits the usual properties of A-tract B-DNA under biologically relevant conditions: large propeller twist of base-pairs, narrowed minor groove, and a straight helix axis. Groove narrowing is fully developed in the A-A-A-A region, but not in the A-A-A region, which previous investigators have proposed as being too short to exhibit typical A-tract properties. The RNA/DNA hybrid produced by HIV reverse transcriptase during (-) strand synthesis presumably forms a "heteromerous" or H-helix with narrower minor groove than an A-helical RNA/RNA duplex. If the narrowing of minor groove in A-tract H-helices is comparable to that seen in A-tract B-helices, then the narrowed minor groove of the polypurine tract could make the second primer site both (1) impervious to RNase H digestion, and (2) susceptible to inhibition by minor groove binding drugs.

The structure of an RNA "kissing" hairpin complex of the HIV TAR hairpin loop and its complement

We have used nuclear magnetic resonance (NMR) to obtain the structure of an RNA "kissing" hairpin complex formed between the HIV-2 TAR hairpin loop and a hairpin with a complementary loop sequence. Kissing hairpins are important in natural antisense reactions; their complex is a specific target for protein binding. The complex has all six nucleotides of each loop paired to form a bent quasicontinuous helix of three coaxially stacked helices: two stems plus a loop-loop interaction helix. Experimental constraints derived from heteronuclear and homonuclear NMR data on 13C and 15N-labeled RNA led to a structure for the loop-loop helix with an average root-mean-square deviation of 0.83 (+/-0.10) A for 33 converged structures relative to the average structure. The loop-loop helix of the kissing complex is distorted compared to A-form RNA. Its major groove is blocked by the phosphodiester bonds that connect the first loop residue of each hairpin with its own stem, and it is flanked by two negatively charged phosphate clusters. The loop-loop helix has alternating helical twists between adjacent base-pairs. The base-pairs at the helix junctions are overwound and three base-pairs near the helix junctions adopt high propeller twists. All these changes reduce the distance needed for the bridging phosphodiester bonds connecting each stem and loop to cross the major groove of the loop-loop helix, and result in a deformed RNA helix with localized perturbations in the minor groove surface. The alternating helical twist pattern, plus other distortions in the loop-loop helix may be important for Rom protein recognition of the kissing hairpin complex.