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

Machine learning shows torsion angle preferences in left-handed and right-handed quadruplex DNAs

Left-handed G quadruplexes (LHG4) have been recently discovered as a new class of G quadruplexes. The biological functions of LHG4s are still unknown, but they share a striking resemblance to Z-DNA in their helicity and jagged phosphate backbone. To further understand structural features of the LHG4s that define their left handedness, we have employed human-interpretable machine-learning methods to classify right- and left-handed G4s purely based on torsional angle analysis. Our results reveal the importance of the α, β, δ, and χ angles in left-handed structuring across both Z-DNAs and LHG4s. Our analysis may serve as the first step to understanding the conditions of formation for LHG4s and their potential biological relevance.

RNA polymerase II trapped on a molecular treadmill: Structural basis of persistent transcriptional arrest by a minor groove DNA binder

Elongating RNA polymerase II (Pol II) can be paused or arrested by a variety of obstacles. These obstacles include DNA lesions, DNA-binding proteins, and small molecules. Hairpin pyrrole-imidazole (Py-Im) polyamides bind to the minor groove of DNA in a sequence-specific manner and induce strong transcriptional arrest. Remarkably, this Py-Im-induced Pol II transcriptional arrest is persistent and cannot be rescued by transcription factor TFIIS. In contrast, TFIIS can effectively rescue the transcriptional arrest induced by a nucleosome barrier. The structural basis of Py-Im-induced transcriptional arrest and why TFIIS cannot rescue this arrest remain elusive. Here we determined the X-ray crystal structures of four distinct Pol II elongation complexes (Pol II ECs) in complex with hairpin Py-Im polyamides as well as of the hairpin Py-Im polyamides-dsDNA complex. We observed that the Py-Im oligomer directly interacts with RNA Pol II residues, introduces compression of the downstream DNA duplex, prevents Pol II forward translocation, and induces Pol II backtracking. These results, together with biochemical studies, provide structural insight into the molecular mechanism by which Py-Im blocks transcription. Our structural study reveals why TFIIS fails to promote Pol II bypass of Py-Im-induced transcriptional arrest.

The Impact of the HydroxyMethylCytosine epigenetic signature on DNA structure and function

We present a comprehensive, experimental and theoretical study of the impact of 5-hydroxymethylation of DNA cytosine. Using molecular dynamics, biophysical experiments and NMR spectroscopy, we found that Ten-Eleven translocation (TET) dioxygenases generate an epigenetic variant with structural and physical properties similar to those of 5-methylcytosine. Experiments and simulations demonstrate that 5-methylcytosine (mC) and 5-hydroxymethylcytosine (hmC) generally lead to stiffer DNA than normal cytosine, with poorer circularization efficiencies and lower ability to form nucleosomes. In particular, we can rule out the hypothesis that hydroxymethylation reverts to unmodified cytosine physical properties, as hmC is even more rigid than mC. Thus, we do not expect dramatic changes in the chromatin structure induced by differences in physical properties between d(mCpG) and d(hmCpG). Conversely, our simulations suggest that methylated-DNA binding domains (MBDs), associated with repression activities, are sensitive to the substitution d(mCpG) ➔ d(hmCpG), while MBD3 which has a dual activation/repression activity is not sensitive to the d(mCpG) d(hmCpG) change. Overall, while gene activity changes due to cytosine methylation are the result of the combination of stiffness-related chromatin reorganization and MBD binding, those associated to 5-hydroxylation of methylcytosine could be explained by a change in the balance of repression/activation pathways related to differential MBD binding.

In Eukaryotic cells, DNA epigenetic modifications play an important role in gene expression and regulation, and protein recognition. In this work we investigate the physical implications of cytosine 5-hydroxymethylation on DNA, its structural and flexibility differences with methylated and unmodified cytosine using molecular dynamics, biophysical experiments and NMR spectroscopy.

In particular the effect of hydroxyl group on free energy of nucleosome and Methyl binding Protein (MBD) binding, comparing in silico and experimental data to shed light on the effect of the reduced flexibility and the direct protein-DNA recognition.

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.

Accurate geometrical restraints for Watson-Crick base pairs

Geometrical restraints provide key structural information for the determination of biomolecular structures at lower resolution by experimental methods such as crystallography or cryo-electron microscopy. In this work, restraint targets for nucleic acids bases are derived from three different sources and compared: small-molecule crystal structures in the Cambridge Structural Database (CSD), ultrahigh-resolution structures in the Protein Data Bank (PDB) and quantum-mechanical (QM) calculations. The best parameters are those based on CSD structures. After over two decades, the standard library of Parkinson et al. [(1996), Acta Cryst. D52, 57-64] is still valid, but improvements are possible with the use of the current CSD database. The CSD-derived geometry is fully compatible with Watson-Crick base pairs, as comparisons with QM results for isolated and paired bases clearly show that the CSD targets closely correspond to proper base pairing. While the QM results are capable of distinguishing between single and paired bases, their level of accuracy is, on average, nearly two times lower than for the CSD-derived targets when gauged by root-mean-square deviations from ultrahigh-resolution structures in the PDB. Nevertheless, the accuracy of QM results appears sufficient to provide stereochemical targets for synthetic base pairs where no reliable experimental structural information is available. To enable future tests for this approach, QM calculations are provided for isocytosine, isoguanine and the iCiG base pair.

DNA conformational transitions inferred from re-evaluation of m|Fo| - D|Fc| electron-density maps

Conformational flexibility of DNA plays important roles in biological processes such as transcriptional regulation and DNA packaging etc. To understand the mechanisms of these processes, it is important to analyse when, where and how DNA shows conformational variations. Recent analyses have indicated that conventional refinement methods do not always provide accurate models of crystallographic heterogeneities and that some information on polymorphism has been overlooked in previous crystallographic studies. In the present study, the m|Fo| - D|Fc| electron-density maps of double-helical DNA crystal structures were calculated at a resolution equal to or better than 1.5 Å and potential conformational transitions were found in 27% of DNA phosphates. Detailed analyses of the m|Fo| - D|Fc| peaks indicated that some of these unassigned densities correspond to ZI ↔ ZII or A/B → BI conformational transitions. A relationship was also found between ZI/ZII transitions and metal coordination in Z-DNA from the detected peaks. The present study highlights that frequent transitions of phosphate backbones occur even in crystals and that some of these transitions are affected by the local molecular environment.

Influence of BII Backbone Substates on DNA Twist: A Unified View and Comparison of Simulation and Experiment for All 136 Distinct Tetranucleotide Sequences

Reliable representation of the B-DNA base-pair step twist is one of the crucial requirements for theoretical modeling of DNA supercoiling and other biologically relevant phenomena in B-DNA. It has long been suspected that the twist is inaccurately described by current empirical force fields. Unfortunately, comparison of simulation results with experiments is not straightforward because of the presence of BII backbone substates, whose populations may differ in experimental and simulation ensembles. In this work, we provide a comprehensive view of the effect of BII substates on the overall B-DNA helix twist and show how to reliably compare twist values from experiment and simulation in two scenarios. First, for longer DNA segments freely moving in solution, we show that sequence-averaged twists of different BI/BII ensembles can be compared directly because of approximate cancellation of the opposing BII effects. Second, for sequence-specific data, such as a particular base-pair step or tetranucleotide twist, can be compared only for a clearly defined BI/BII backbone conformation. For the purpose of force field testing, we designed a compact set of fourteen 22-base-pair B-DNA duplexes (Set 14) containing all 136 distinct tetranucleotide sequences and carried out a total of 84 μs of molecular dynamics simulations, primarily with the OL15 force field. Our results show that the ff99bsc0εζ_OL1χ_OL4, parmbsc1, and OL15 force fields model the B-DNA helical twist in good agreement with X-ray and minicircle ligation experiments. The comprehensive understanding obtained regarding the effect of BII substates on the base-pair step geometry should aid meaningful comparisons of various conformational ensembles in future research.

A targeted DNA substrate mechanism for the inhibition of HIV-1 integrase by inhibitors with antiretroviral activity

We recently reported that viral DNA could be the primary target of raltegravir (RAL), an efficient anti‐HIV‐1 drug, which acts by inhibiting integrase. To elucidate this mechanism, we conducted a comparative analysis of RAL and TB11, a diketoacid abandoned as an anti‐HIV‐1 drug for its weak efficiency and marked toxicity, and tested the effects of the catalytic cofactor Mg²⁺ (5 mm) on drug‐binding properties. We used circular dichroism and fluorescence to determine drug affinities for viral DNA long terminal repeats (LTRs) and peptides derived from the integrase active site and DNA retardation assays to assess drug intercalation into DNA base pairs. We found that RAL bound more tightly to LTR ends than did TB11 (a diketo acid bearing an azido group) and that Mg²⁺ significantly increased the affinity of both RAL and TB11. We also observed a good relationship between drug binding with processed LTR and strand transfer inhibition. This unusual type of inhibition was caused by Mg²⁺‐assisted binding of drugs to DNA substrate, rather than to enzyme. Notably, while RAL bound exclusively to the cleavable/cleaved site, TB11 further intercalated into DNA base pairs and interacted with the integrase‐derived peptides. These unwanted binding sites explain the weaker bioavailability and higher toxicity of TB11 compared with the more effective RAL.

DNA polymorphism in crystals: three stable conformations for the decadeoxynucleotide d(GCATGCATGC)

High-resolution structures of DNA fragments determined using X-ray crystallography or NMR have provided descriptions of a veritable alphabet of conformations. They have also shown that DNA is a flexible molecule, with some sequences capable of adopting two different structures. Here, the first example is presented of a DNA fragment that can assume three different and distinct conformations in crystals. The decanucleotide d(GCATGCATGC) was previously reported to assume a single-stranded double-fold structure. In one of the two crystal structures described here the decamer assumes both the double-fold conformation and, simultaneously, the more conventional B-type double-helical structure. In the other crystal the sequence assumes the A-type double-helical conformation. These results, taken together with CD spectra, which were recorded as the decamer was titrated against four metal ions and spermine, indicate that the molecule may exist as a mixed population of structures in solution. Small differences in the environmental conditions, such as the concentration of metal ion, may decide which of these crystallizes out. The results also support the idea that it may be possible for DNA to change its structure to suit the binding requirements of proteins or drugs.

Unraveling the sequence-dependent polymorphic behavior of d(CpG) steps in B-DNA

We have made a detailed study of one of the most surprising sources of polymorphism in B-DNA: the high twist/low twist (HT/LT) conformational change in the d(CpG) base pair step. Using extensive computations, complemented with database analysis, we were able to characterize the twist polymorphism in the d(CpG) step in all the possible tetranucleotide environment. We found that twist polymorphism is coupled with BI/BII transitions, and, quite surprisingly, with slide polymorphism in the neighboring step. Unexpectedly, the penetration of cations into the minor groove of the d(CpG) step seems to be the key element in promoting twist transitions. The tetranucleotide environment also plays an important role in the sequence-dependent d(CpG) polymorphism. In this connection, we have detected a previously unexplored intramolecular C-H···O hydrogen bond interaction that stabilizes the low twist state when 3'-purines flank the d(CpG) step. This work explains a coupled mechanism involving several apparently uncorrelated conformational transitions that has only been partially inferred by earlier experimental or theoretical studies. Our results provide a complete description of twist polymorphism in d(CpG) steps and a detailed picture of the molecular choreography associated with this conformational change.

B-DNA characteristics are preserved in double stranded d(A)3·d(T)3 and d(G)3·d(C)3 mini-helixes: conclusions from DFT/M06-2X study

We report the results of the first comprehensive DFT study on the d(A)3·d(T)3 and d(G)3·d(C)3 nucleic acid duplexes. The ability of mini-helixes to preserve the conformation of B-DNA in the gas phase and under the influence of such factors as: solvent, uncompensated charge, and counter-ions was evaluated using M06-2X functional with 6-31G(d,p) basis set. The accuracy of the models was ascertained based on their ability to reproduce key structural features of natural B-DNA. Analysis of the helicity suggests that the helical conformations adopt geometrical parameters which are close to those of the B-DNA form. The torsion angles fall somewhere between the values observed for BI/BII conformational classes. The comparative analysis of parameters of isolated Watson-Crick base pairs versus B-DNA-like conformations indicates the same tendency of base-pair polarization and hydration. Specifically, effects of polarization of nucleobases in continuum type dielectric medium mimicking water are stronger than those caused by the presence of backbone. Polar environment as well as the presence of counterions stabilizes duplexes, facilitating helix formation. Substantial conformational changes of nucleotides upon duplex formation decrease the binding energy. In spite of structural and energetic changes, the placement of a mini-helix into the gas phase does not lead to significant disruption of the structure. On the contrary, the duplex preserves its helicity and the strands remain bound.

DNA conformational flexibility study using phosphate backbone neutralization model

Due to the critical role of DNA in the processes of the cell cycle, the structural and physicochemical properties of DNA have long been of concern. In the present work, the effect of interplay between the DNA duplex and metal ions in solution on the DNA structure and conformational flexibility is studied by comparing the structure and dynamic conformational behavior of a duplex in a normal form and its “null isomer” using molecular dynamics methods. It was found that the phosphate neutralization changes the cation atmosphere around the DNA duplex greatly, increases the major groove width, decreases the minor groove width, and reduces the global bending direction preference. We also noted that the probability of BI phosphate linkages increases significantly because of the charge reduction in the backbone phosphate groups. More importantly, we found that the electrostatic effect on the DNA conformational flexibility is dependent on the sequence; that is, the phosphate backbone neutralization induces the global dynamic bending to be less flexible for GC-rich sequences but more flexible for AT-rich sequences.

Site-specific stabilization of DNA by a tethered major groove amine, 7-aminomethyl-7-deaza-2'-deoxyguanosine

A cationic 7-aminomethyl-7-deaza-2′-deoxyguanosine (7amG) was incorporated site-specifically into the self-complementary duplex d(G¹A²G³A⁴X⁵C⁶G⁷C⁸T⁹C¹⁰T¹¹C¹²)₂ (X = 7amG). This construct placed two positively charged amines adjacent to the major groove edges of two symmetry-related guanines, providing a model for probing how cation binding in the major groove modulates the structure and stability of DNA. Molecular dynamics calculations restrained by nuclear magnetic resonance (NMR) data revealed that the tethered cationic amines were in plane with the modified base pairs. The tethered amines did not form salt bridges to the phosphodiester backbone. There was also no indication of the amines being capable of hydrogen bonding to flanking DNA bases. NMR spectroscopy as a function of temperature revealed that the X⁵ imino resonance remained sharp at 55 °C. Additionally, two 5′-neighboring base pairs, A⁴:T⁹ and G³:C¹⁰, were stabilized with respect to the exchange of their imino protons with solvent. The equilibrium constant for base pair opening at the A⁴:T⁹ base pair determined by magnetization transfer from water in the absence and presence of added ammonia base catalyst decreased for the modified duplex compared to that of the A⁴:T⁹ base pair in the unmodified duplex, which confirmed that the overall fraction of the A⁴:T⁹ base pair in the open state of the modified duplex decreased. This was also observed for the G³:C¹⁰ base pair, where αK_op for the G³:C¹⁰ base pair in the modified duplex was 3.0 × 10⁶ versus 4.1 × 10⁶ for the same base pair in the unmodified duplex. In contrast, equilibrium constants for base pair opening at the X⁵:C⁸ and C⁶:G⁷ base pairs did not change at 15 °C. These results argue against the notion that electrostatic interactions with DNA are entirely entropic and suggest that major groove cations can stabilize DNA via enthalpic contributions to the free energy of duplex formation.

Effect of initial ion positions on the interactions of monovalent and divalent ions with a DNA duplex as revealed with atomistic molecular dynamics simulations

Monovalent (Na(+)) and divalent (Mg(2+)) ion distributions around the Dickerson-Drew dodecamer were studied by atomistic molecular dynamics (MD) simulations with AMBER molecular modeling software. Different initial placements of ions were tried and the resulting effects on the ion distributions around DNA were investigated. For monovalent ions, results were found to be nearly independent of initial cation coordinates. However, Mg(2+) ions demonstrated a strong initial coordinate dependent behavior. While some divalent ions initially placed near the DNA formed essentially permanent direct coordination complexes with electronegative DNA atoms, Mg(2+) ions initially placed further away from the duplex formed a full, nonexchanging, octahedral first solvation shell. These fully solvated cations were still capable of binding with DNA with events lasting up to 20 ns, and in comparison were bound much longer than Na(+) ions. Force field parameters were also investigated with modest and little differences arising from ion (ions94 and ions08) and nucleic acid description (ff99, ff99bsc0, and ff10), respectively. Based on known Mg(2+) ion solvation structure, we conclude that in most cases Mg(2+) ions retain their first solvation shell, making only solvent-mediated contacts with DNA duplex. The proper way to simulate Mg(2+) ions around DNA duplex, therefore, should begin with ions placed in the bulk water.

Pyrrole-imidazole polyamides distinguish between double-helical DNA and RNA

Groove specificity: pyrrole-imidazole polyamides are well-known for their specific interactions with the minor groove of DNA. However, polyamides do not show similar binding to duplex RNA, and a structural rationale for the molecular-level discrimination of nucleic acid duplexes by minor-groove-binding ligands is presented.

Mg2+ in the major groove modulates B-DNA structure and dynamics

This study investigates the effect of Mg(2+) bound to the DNA major groove on DNA structure and dynamics. The analysis of a comprehensive dataset of B-DNA crystallographic structures shows that divalent cations are preferentially located in the DNA major groove where they interact with successive bases of (A/G)pG and the phosphate group of 5'-CpA or TpG. Based on this knowledge, molecular dynamics simulations were carried out on a DNA oligomer without or with Mg(2+) close to an ApG step. These simulations showed that the hydrated Mg(2+) forms a stable intra-strand cross-link between the two purines in solution. ApG generates an electrostatic potential in the major groove that is particularly attractive for cations; its intrinsic conformation is well-adapted to the formation of water-mediated hydrogen bonds with Mg(2+). The binding of Mg(2+) modulates the behavior of the 5'-neighboring step by increasing the BII (ε-ζ>0°) population of its phosphate group. Additional electrostatic interactions between the 5'-phosphate group and Mg(2+) strengthen both the DNA-cation binding and the BII character of the 5'-step. Cation binding in the major groove may therefore locally influence the DNA conformational landscape, suggesting a possible avenue for better understanding how strong DNA distortions can be stabilized in protein-DNA complexes.

Reference simulations of noncanonical nucleic acids with different χ variants of the AMBER force field: quadruplex DNA, quadruplex RNA and Z-DNA

Refinement of empirical force fields for nucleic acids requires their extensive testing using as wide range of systems as possible. However, finding unambiguous reference data is not easy. In this paper, we analyze four systems which we suggest should be included in standard portfolio of molecules to test nucleic acids force fields, namely, parallel and antiparallel stranded DNA guanine quadruplex stems, RNA quadruplex stem, and Z-DNA. We highlight parameters that should be monitored to assess the force field performance. The work is primarily based on 8.4 μs of 100-250 ns trajectories analyzed in detail followed by 9.6 μs of additional selected back up trajectories that were monitored to verify that the results of the initial analyses are correct. Four versions of the Cornell et al. AMBER force field are tested, including an entirely new parmχ(OL4) variant with χ dihedral specifically reparametrized for DNA molecules containing syn nucleotides. We test also different water models and ion conditions. While improvement for DNA quadruplexes is visible, the force fields still do not fully represent the intricate Z-DNA backbone conformation.

Structure of the tetradecanucleotide d(CCCCGGTACCGGGG)2 as an A-DNA duplex

The crystal structure of the tetradecanucleotide sequence d(CCCCGGTACCGGGG)(2) has been determined at 2.5 Å resolution in the tetragonal space group P4(1). This sequence was designed with the expectation of a four-way junction. However, the sequence crystallized as an A-DNA duplex and represents more than one full turn of the A-helix. The crystallographic asymmetric unit consists of one tetradecanucleotide duplex. The structural parameters of the A-type DNA duplex structure and the crystal-packing arrangement are described. One Mn(2+) ion was identified with direct coordination to the N7 position of G(13) and a water molecule at the major-groove side of the C(2)·G(13) base pair.

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.

B-DNA structure is intrinsically polymorphic: even at the level of base pair positions

Increasingly exact measurement of single crystal X-ray diffraction data offers detailed characterization of DNA conformation, hydration and electrostatics. However, instead of providing a more clear and unambiguous image of DNA, highly accurate diffraction data reveal polymorphism of the DNA atomic positions and conformation and hydration. Here we describe an accurate X-ray structure of B-DNA, painstakingly fit to a multistate model that contains multiple competing positions of most of the backbone and of entire base pairs. Two of ten base-pairs of CCAGGCCTGG are in multiple states distinguished primarily by differences in slide. Similarly, all the surrounding ions are seen to fractionally occupy discrete competing and overlapping sites. And finally, the vast majority of water molecules show strong evidence of multiple competing sites. Conventional resolution appears to give a false sense of homogeneity in conformation and interactions of DNA. In addition, conventional resolution yields an average structure that is not accurate, in that it is different from any of the multiple discrete structures observed at high resolution. Because base pair positional heterogeneity has not always been incorporated into model-building, even some high and ultrahigh-resolution structures of DNA do not indicate the full extent of conformational polymorphism.

Thrombin-aptamer recognition: a revealed ambiguity

Aptamers are structured oligonucleotides that recognize molecular targets and can function as direct protein inhibitors. The best-known example is the thrombin-binding aptamer, TBA, a single-stranded 15-mer DNA that inhibits the activity of thrombin, the key enzyme of coagulation cascade. TBA folds as a G-quadruplex structure, as proved by its NMR structure. The X-ray structure of the complex between TBA and human α-thrombin was solved at 2.9-Å resolution, but did not provide details of the aptamer conformation and the interactions with the protein molecule. TBA is rapidly processed by nucleases. To improve the properties of TBA, a number of modified analogs have been produced. In particular, a modified TBA containing a 5′-5′ polarity inversion site, mTBA, has higher stability and higher affinity toward thrombin with respect to TBA, although it has a lower inhibitory activity. We present the crystal structure of the thrombin–mTBA complex at 2.15-Å resolution; the resulting model eventually provides a clear picture of thrombin–aptamers interaction, and also highlights the structural bases of the different properties of TBA and mTBA. Our findings open the way for a rational design of modified aptamers with improved potency as anticoagulant drugs.

Bifacial nucleoside as a surrogate for both T and A in duplex DNA

A new functional bifacial nucleoside derived from 7-aminopyrimido[4,5-d]pyrimidine-2,4(1H,3H)-dione, a Janus-type nucleobase, has been synthesized and incorporated into DNA oligonucleotides. The nucleobase, having self-complementary H-bonding faces, mimics both T and A and engages in the corresponding Watson-Crick-like base pairs, forming stable duplexes.

High regularity of Z-DNA revealed by ultra high-resolution crystal structure at 0.55 A

The crystal structure of a Z-DNA hexamer duplex d(CGCGCG)(2) determined at ultra high resolution of 0.55 Å and refined without restraints, displays a high degree of regularity and rigidity in its stereochemistry, in contrast to the more flexible B-DNA duplexes. The estimations of standard uncertainties of all individually refined parameters, obtained by full-matrix least-squares optimization, are comparable with values that are typical for small-molecule crystallography. The Z-DNA model generated with ultra high-resolution diffraction data can be used to revise the stereochemical restraints applied in lower resolution refinements. Detailed comparisons of the stereochemical library values with the present accurate Z-DNA parameters, shows in general a good agreement, but also reveals significant discrepancies in the description of guanine-sugar valence angles and in the geometry of the phosphate groups.

Structural basis for cyclic Py-Im polyamide allosteric inhibition of nuclear receptor binding

Pyrrole-imidazole polyamides are a class of small molecules that can be programmed to bind a broad repertoire of DNA sequences, disrupt transcription factor−DNA interfaces, and modulate gene expression pathways in cell culture experiments. In this paper we describe a high-resolution X-ray crystal structure of a β-amino turn-linked eight-ring cyclic Py-Im polyamide bound to the central six base pairs of the sequence d(5′-CCAGTACTGG-3′)₂, revealing significant modulation of DNA shape. We compare the DNA structural perturbations induced by DNA-binding transcripton factors, androgen receptor and glucocorticoid receptor, in the major groove to those induced by cyclic polyamide binding in the minor groove. The cyclic polyamide is an allosteric modulator that perturbs the DNA structure in such a way that nuclear receptor protein binding is no longer compatible. This allosteric perturbation of the DNA helix provides a molecular basis for disruption of transcription factor−DNA interfaces by small molecules, a minimum step in chemical control of gene networks.

Sequence-specific ultrasonic cleavage of DNA

We investigated the phenomenon of ultrasonic cleavage of DNA by analyzing a large set of cleavage patterns of DNA restriction fragments using polyacrylamide gel electrophoresis. The cleavage intensity of individual phosphodiester bonds was found to depend on the nucleotide sequence and the position of the bond with respect to the ends of the fragment. The relative intensities of cleavage of the central phosphodiester bond in 16 dinucleotides and 256 tetranucleotides were determined by multivariate statistical analysis. We observed a remarkable enhancement of the mean values of the relative intensities of cleavage (cleavage rates) in phosphodiester bonds following deoxycytidine, which diminished in the row of dinucleotides: d(CpG) > d(CpA) > d(CpT) >> d(CpC). The cleavage rates for all pairs of complementary dinucleotides were significantly different from each other. The effect of flanking nucleotides in tetranucleotides on cleavage rates of all 16 types of central dinucleotides was also statistically significant. The sequence-dependent ultrasonic cleavage rates of dinucleotides are consistent with reported data on the intensity of the conformational motion of their 5'-deoxyribose. As a measure of local conformational dynamics, cleavage rates may be useful for characterizing functional regions of the genome.

Using soft X-rays for a detailed picture of divalent metal binding in the nucleosome

Divalent metals associate with DNA in a site-selective manner, which can influence nucleosome positioning, mobility, compaction, and recognition by nuclear factors. We previously characterized divalent metal binding in the nucleosome core using hard (short-wavelength) X-rays allowing high-resolution crystallographic determination of the strongest affinity sites, which revealed that Mn(2+) associates with the DNA major groove in a sequence- and conformation-dependent manner. In this study, we obtained diffraction data with soft X-rays at the Mn(2+) absorption edge for a core particle crystal in the presence of 10 mM MnSO(4), mimicking prevailing Mg(2+) concentration in the nucleus. This provides an exceptional view of counterion binding in the nucleosome through identification of 45 divalent metal binding sites. In addition to that at the well-characterized major interparticle interface, only one other histone-divalent metal binding site is found, which corresponds to a symmetry-related counterpart on the 'free' H2B alpha1 helix C-terminus. This emphasizes the importance of the alpha-helix dipole in ion binding and suggests that the H2B motif may serve as a nucleation site in nucleosome compaction. The 43 sites associated with the DNA are characterized by (1) high-affinity direct coordination at the most electrostatically favorable major groove locations, (2) metal hydrate binding to the major groove, (3) direct coordination to phosphate groups at sites of high charge density, (4) metal hydrate binding in the minor groove, or (5) metal hydrate-divalent anion pairing. Metal hydrates are found within the minor groove only at locations displaying a narrow range of high-intermediate width and to which histone N-terminal tails are not associated or proximal. This indicates that divalent metals and histone tails can both collaborate and compete in minor groove association, which sheds light on nucleosome solubility and chromatin compaction behavior.

Specificity of LTR DNA recognition by a peptide mimicking the HIV-1 integrase {alpha}4 helix

HIV-1 integrase integrates retroviral DNA through 3′-processing and strand transfer reactions in the presence of a divalent cation (Mg²⁺ or Mn²⁺). The α4 helix exposed at the catalytic core surface is essential to the specific recognition of viral DNA. To define group determinants of recognition, we used a model composed of a peptide analogue of the α4 helix, oligonucleotides mimicking processed and unprocessed U5 LTR end and 5 mM Mg²⁺. Circular dichroism, fluorescence and NMR experiments confirmed the implication of the α4 helix polar/charged face in specific and non-specific bindings to LTR ends. The specific binding requires unprocessed LTR ends—i.e. an unaltered 3′-processing site CA↓GT3′—and is reinforced by Mg²⁺ (K_d decreases from 2 to 0.8 nM). The latter likely interacts with the ApG and GpT3′ steps of the 3′-processing site. With deletion of GT3′, only persists non-specific binding (K_d of 100 μM). Proton chemical shift deviations showed that specific binding need conserved amino acids in the α4 helix and conserved nucleotide bases and backbone groups at LTR ends. We suggest a conserved recognition mechanism based on both direct and indirect readout and which is subject to evolutionary pressure.

Evaluation of parameters critical for observing nucleic acids inside living Xenopus laevis oocytes by in-cell NMR spectroscopy

In-cell NMR spectroscopy of proteins in different cellular environments is a well-established technique that, however, has not been applied to nucleic acids so far. Here, we show that isotopically labeled DNA and RNA can be observed inside the eukaryotic environment of Xenopus laevis oocytes by in-cell NMR spectroscopy. One limiting factor for the observation of nucleic acids in Xenopus oocytes is their reduced stability. We demonstrate that chemical modification of DNA and RNA can protect them from degradation and can significantly enhance their lifetime. Finally, we show that the imino region of the NMR spectrum is devoid of any oocyte background signals enabling the detection even of isotopically nonlabeled molecules.

Polarizable atomic multipole X-ray refinement: application to peptide crystals

Recent advances in computational chemistry have produced force fields based on a polarizable atomic multipole description of biomolecular electrostatics. In this work, the Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field is applied to restrained refinement of molecular models against X-ray diffraction data from peptide crystals. A new formalism is also developed to compute anisotropic and aspherical structure factors using fast Fourier transformation (FFT) of Cartesian Gaussian multipoles. Relative to direct summation, the FFT approach can give a speedup of more than an order of magnitude for aspherical refinement of ultrahigh-resolution data sets. Use of a sublattice formalism makes the method highly parallelizable. Application of the Cartesian Gaussian multipole scattering model to a series of four peptide crystals using multipole coefficients from the AMOEBA force field demonstrates that AMOEBA systematically underestimates electron density at bond centers. For the trigonal and tetrahedral bonding geometries common in organic chemistry, an atomic multipole expansion through hexadecapole order is required to explain bond electron density. Alternatively, the addition of interatomic scattering (IAS) sites to the AMOEBA-based density captured bonding effects with fewer parameters. For a series of four peptide crystals, the AMOEBA-IAS model lowered R(free) by 20-40% relative to the original spherically symmetric scattering model.

A rare nucleotide base tautomer in the structure of an asymmetric DNA junction

The single-crystal structure of a DNA Holliday junction assembled from four unique sequences shows a structure that conforms to the general features of models derived from similar constructs in solution. The structure is a compact stacked-X form junction with two sets of stacked B-DNA-type arms that coaxially stack to form semicontinuous duplexes interrupted only by the crossing of the junction. These semicontinuous helices are related by a right-handed rotation angle of 56.5 degrees, which is nearly identical to the 60 degree angle in the solution model but differs from the more shallow value of approximately 40 degrees for previous crystal structures of symmetric junctions that self-assemble from single identical inverted-repeat sequences. This supports the model in which the unique set of intramolecular interactions at the trinucleotide core of the crossing strands, which are not present in the current asymmetric junction, affects both the stability and geometry of the symmetric junctions. An unexpected result, however, is that a highly wobbled A.T base pair, which is ascribed here to a rare enol tautomer form of the thymine, was observed at the end of a CCCC/GGGG sequence within the stacked B-DNA arms of this 1.9 A resolution structure. We suggest that the junction itself is not responsible for this unusual conformation but served as a vehicle for the study of this CG-rich sequence as a B-DNA duplex, mimicking the form that would be present in a replication complex. The existence of this unusual base lends credence to and defines a sequence context for the "rare tautomer hypothesis" as a mechanism for inducing transition mutations during DNA replication.

Small local variations in B-form DNA lead to a large variety of global geometries which can accommodate most DNA-binding protein motifs

BACKGROUND

An important question of biological relevance is the polymorphism of the double-helical DNA structure in its free form, and the changes that it undergoes upon protein-binding. We have analysed a database of free DNA crystal structures to assess the inherent variability of the free DNA structure and have compared it with a database of protein-bound DNA crystal structures to ascertain the protein-induced variations.

RESULTS

Most of the dinucleotide steps in free DNA display high flexibility, assuming different conformations in a sequence-dependent fashion. With the exception of the AA/TT and GA/TC steps, which are 'A-phobic', and the GG/CC step, which is 'A-philic', the dinucleotide steps show no preference for A or B forms of DNA. Protein-bound DNA adopts the B-conformation most often. However, in certain cases, protein-binding causes the DNA backbone to take up energetically unfavourable conformations. At the gross structural level, several protein-bound DNA duplexes are observed to assume a curved conformation in the absence of any large distortions, indicating that a series of normal structural parameters at the dinucleotide and trinucleotide level, similar to the ones in free B-DNA, can give rise to curvature at the overall level.

CONCLUSION

The results illustrate that the free DNA molecule, even in the crystalline state, samples a large amount of conformational space, encompassing both the A and the B-forms, in the absence of any large ligands. A-form as well as some non-A, non-B, distorted geometries are observed for a small number of dinucleotide steps in DNA structures bound to the proteins belonging to a few specific families. However, for most of the bound DNA structures, across a wide variety of protein families, the average step parameters for various dinucleotide sequences as well as backbone torsion angles are observed to be quite close to the free 'B-like' DNA oligomer values, highlighting the flexibility and biological significance of this structural form.

Sheng Ding: using chemical methods to control cell fate

Profiles provide insights into the lives, backgrounds, career paths, and futures of scientists who serve as Experts on ACS Chemical Biology's online Ask the Expert feature. Ding will begin answering your questions in mid-March, 2009. Readers are encouraged to submit their questions to the Experts at http://community.acs.org/chembiol/ . The editors will post the most interesting exchanges on the web site.

Functionalized 2'-amino-alpha-L-LNA: directed positioning of intercalators for DNA targeting

Chemically modified oligonucleotides are increasingly applied in nucleic acid based therapeutics and diagnostics. LNA (locked nucleic acid) and its diastereomer alpha-L-LNA are two promising examples thereof that exhibit increased thermal and enzymatic stability. Herein, the synthesis, biophysical characterization, and molecular modeling of N2'-functionalized 2'-amino-alpha-L-LNA is described. Chemoselective N2'-functionalization of protected amino alcohol 1 followed by phosphitylation afforded a structurally varied set of target phosphoramidites, which were incorporated into oligodeoxyribonucleotides. Incorporation of pyrene-functionalized building blocks such as 2'-N-(pyren-1-yl)carbonyl-2'-amino-alpha-L-LNA (monomer X) led to extraordinary increases in thermal affinity of up to +19.5 degrees C per modification against DNA targets in particular. In contrast, incorporation of building blocks with small nonaromatic N2'-functionalities such as 2'-N-acetyl-2'-amino-alpha-L-LNA (monomer V) had detrimental effects on thermal affinity toward DNA/RNA complements with decreases of as much as -16.5 degrees C per modification. Extensive thermal DNA selectivity, favorable entropic contributions upon duplex formation, hybridization-induced bathochromic shifts of pyrene absorption maxima and increases in circular dichroism signal intensity, and molecular modeling studies suggest that pyrene-functionalized 2'-amino-alpha-L-LNA monomers W-Y having short linkers between the bicyclic skeleton and the pyrene moiety allow high-affinity hybridization with DNA complements and precise positioning of intercalators in nucleic acid duplexes. This rigorous positional control has been utilized for the development of probes for emerging therapeutic and diagnostic applications focusing on DNA targeting.

Crystallographic analysis of a sex-specific enhancer element: sequence-dependent DNA structure, hydration, and dynamics

The crystal structure of a sex-specific enhancer element is described at a resolution of 1.6 A. This 16-bp site, designated Dsx(A), functions in the regulation of a genetic switch between male and female patterns of gene expression in Drosophila melanogaster. Related sites are broadly conserved in metazoans, including in the human genome. This enhancer element is unusually rich in general regulatory sequences related to DNA recognition by multiple classes of eukaryotic transcription factors, including the DM motifs, homeodomain, and high mobility group box. Whereas free DNA is often crystallized as an A-form double helix, Dsx(A) was crystallized as B-DNA and thus provides a model for the prebound conformation of diverse regulatory DNA complexes. Sequence-dependent conformational properties that extend features of shorter B-DNA fragments with respect to double helical parameters, groove widths, hydration, and binding of divalent metal ions are observed. The structure also exhibits a sequence-dependent pattern of isotropic thermal B-factors, suggesting possible variation in the local flexibility of the DNA backbone. Such fluctuations are in accord with structural variability observed in prior B-DNA structures. We speculate that sites of intrinsic flexibility within a DNA control element provide hinges for its protein-directed reorganization in a transcriptional preinitiation complex.

Transition-state analysis of the DNA repair enzyme MutY

The transition state (TS) structure of MutY-catalyzed DNA hydrolysis was solved using multiple kinetic isotope effect (KIE) measurements. MutY is a base excision repair enzyme which cleaves adenine from 8-oxo-G:A mismatches in vivo, and also from G:A mismatches in vitro. TS analysis of G:A-DNA hydrolysis revealed a stepwise S(N)1 (D(N)*A(N)(double dagger)) mechanism proceeding through a highly reactive oxacarbenium ion intermediate which would have a lifetime in solution of <10(-10) s. C-N bond cleavage is reversible; the N-glycoside bond breaks and reforms repeatedly before irreversible water attack on the oxacarbenium ion. KIEs demonstrated that MutY uses general acid catalysis by protonating N7. It enforces a 3'-exo sugar ring conformation and other sugar ring distortions to stabilize the oxacarbenium ion. Combining the experimental TS structure with the previously reported crystal structure of an abortive Michaelis complex elucidates the step-by-step catalytic sequence.

Electrostatic versus nonelectrostatic effects in DNA sequence discrimination by divalent ions Mg2+ and Mn2+

Mg2+ and Mn2+ ions are critical to the functioning of phosphoryl transfer enzymes, such as restriction endonucleases. Although these ions play similar roles in the chemical steps, they govern substrate specificity via modulating sequence discrimination by up to a factor of 10(5) with Mg2+ and only up to a factor of 10 with Mn2+. To explain whether such diversity originates in fundamental differences in the electronic structures of the nucleobase-hydrated-metal ion complexes, structures and interaction energies were determined at the density functional (DFT) and second-order Møller-Plesset (MP2) levels of theory. Although both metal ions favor identical binding sites, Mn2+ complexes exhibit greater distortions from the ideal octahedral geometry and larger variability than the corresponding Mg2+ systems. In inner-shell complexes, with direct contact between the metal and the nucleobase, Mg2+ is preferred over Mn2+ in the gas phase, due primarily to nonelectrostatic effects. The interaction energies of the two metal ions are more similar in the outer-shell complexes, likely due to reduced charge transfer between the hydrated metal ion and the base moieties. Inclusion of solvation effects can amplify the relative nucleobase preferences of Mg2+ and Mn2+, indicating that bulk hydration modulates the balance between electrostatic and nonelectrostatic terms. In most cases, the base substitutions in solution are facilitated more by Mn2+ than by Mg2+. Electrostatic properties of the environment were demonstrated to have a major influence on the nucleobase preferences of the two metal ions. Overall, quantum chemical calculations suggest that the contrasting selectivity of Mg2+ and Mn2+ cofactors toward nucleobases derives from the larger flexibility of the Mn2+ complexes accompanied by the excessive polarization and charge-transfer effects as well as less favorable solvation.

Specific Interactions of Divalent Metal Ions with a DNA Duplex Containing the d(CA)n/(GT)nTandem Repeat

Divalent metal ions are essential for maintaining functional states of the DNA molecule. Their participation in DNA structure is modulated by the base sequence and varies depending on the nature of the ion. The present investigation addresses the interaction of Ca2+ ions with a tandem repeat of two CA dinucleotides, (CA)2/(TG)2. The binding of Ca2+ to the repeat is monitored by nuclear magnetic resonance (NMR) spectroscopy using chemical shift mapping. Parallel experiments monitor binding of Mg2+ ions to the repeat as well as binding of each ion to a DNA duplex in which the (CA)2/(TG)2 repeat is eliminated. The results reveal that the direction and the magnitude of chemical shift changes induced by Ca2+ ions in the NMR spectra of the repeat are different from those induced by Mg2+ ions. The differences between the two cations are significantly diminished by the elimination of the (CA)2/(TG)2 repeat. These findings suggest a specific interaction of Ca2+ ions with the (CA)2/(TG)2 motif. The specificity of the interaction resides in the two A-T base pairs of the repeat, and it involves the major groove of the first A-T base pair and both grooves of the second A-T base pair.

Complicated water orientations in the minor groove of the B-DNA decamer d(CCATTAATGG)2 observed by neutron diffraction measurements

It has long been suspected that the structure and function of a DNA duplex can be strongly dependent on its degree of hydration. By neutron diffraction experiments, we have succeeded in determining most of the hydrogen (H) and deuterium (D) atomic positions in the decameric d(CCATTAATGG)₂ duplex. Moreover, the D positions in 27 D₂O molecules have been determined. In particular, the complex water network in the minor groove has been observed in detail. By a combined structural analysis using 2.0 Å resolution X-ray and 3.0 Å resolution neutron data, it is clear that the spine of hydration is built up, not only by a simple hexagonal hydration pattern (as reported in earlier X-ray studies), but also by many other water bridges hydrogen-bonded to the DNA strands. The complexity of the hydration pattern in the minor groove is derived from an extraordinary variety of orientations displayed by the water molecules.

Influence of minor groove substituents on the structure of DNA Holliday junctions

The inosine-containing sequence d(CCIGTACm(5)CGG) is shown to crystallize as a four-stranded DNA junction. This structure is nearly identical to the antiparallel junction formed by the parent d(CCGGTACm(5)()CGG) sequence [Vargason, J. M., and Ho, P. S. (2002) J. Biol. Chem. 277, 21041-21049] in terms of its conformational geometry, and inter- and intramolecular interactions within the DNA and between the DNA and solvent, even though the 2-amino group in the minor groove of the important G(3).m(5)C(8) base pair of the junction core trinucleotide (italicized) has been removed. In contrast, the analogous 2,6-diaminopurine sequence d(CCDGTACTGG) crystallizes as resolved duplex DNAs, just like its parent sequence d(CCAGTACTGG) [Hays, F. A., Vargason, J. M., and Ho, P. S. (2003) Biochemistry 42, 9586-9597]. These results demonstrate that it is not the presence or absence of the 2-amino group in the minor groove of the R(3).Y(8) base pair that specifies whether a sequence forms a junction, but the positions of the extracyclic amino and keto groups in the major groove. Finally, the study shows that the arms of the junction can accommodate perturbations to the B-DNA conformation of the stacked duplex arms associated with the loss of the 2-amino substituent, and that two hydrogen bonding interactions from the C(7) and Y(8) pyrimidine nucleotides to phosphate oxygens of the junction crossover specify the geometry of the Holliday junction.

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.