NMR evidence for mechanical coupling of phosphate B(I)-B(II) transitions with deoxyribose conformational exchange in DNA

Experimental detection of conformational transitions between forms of DNA: problems and prospects

Under different conditions, the DNA double helix can take different geometric forms. Of the large number of its conformations, in addition to the "canonical" B form, the A, C, and Z forms are widely known, and the D, Hoogsteen, and X forms are less known. DNA locally takes the A, C, and Z forms in the cell, in complexes with proteins. We compare different methods for detecting non-canonical DNA conformations: X-ray, IR, and Raman spectroscopy, linear and circular dichroism in both the infrared and ultraviolet regions, as well as NMR (measurement of chemical shifts and their anisotropy, scalar and residual dipolar couplings and inter-proton distances from NOESY (nuclear Overhauser effect spectroscopy) data). We discuss the difficulties in applying these methods, the problems of theoretical interpretation of the experimental results, and the prospects for reliable identification of non-canonical DNA conformations.

Kinetics and thermodynamics of BI-BII interconversion altered by T:G mismatches in DNA

T:G mismatches in DNA result in humans primarily from deamination of methylated CpG sites. They are repaired by redundant systems, such as thymine DNA glycosylase (TDG) and methyl-binding domain enzyme (MBD4), and maintenance of these sites has been implicated in epigenetic processes. The process by which these enzymes identify a canonical DNA base in the incorrect basepairing context remains a mystery. However, the conserved contacts of the repair enzymes with the DNA backbone suggests a role for protein-phosphate interaction in the recognition and repair processes. We have used 31P NMR to investigate the energetics of DNA backbone BI-BII interconversion, and for this work have focused on alterations to the activation barriers to interconversion and the effect of a mismatch compared with canonical DNA. We have found that alterations to the ΔG of interconversion for T:G basepairs are remarkably similar to U:G basepairs in the form of stepwise differences in ΔG of 1-2 kcal/mol greater than equivalent steps in unmodified DNA, suggesting a universality of this result for TDG substrates. Likewise, we see perturbations to the free energy (∼1 kcal/mol) and enthalpy (2-5 kcal/mol) of activation for the BI-BII interconversion localized to the phosphates flanking the mismatch. Overall our results strongly suggest that the perturbed backbone energetics in T:G basepairs play a significant role in the recognition process of DNA repair enzymes.

Variability of the DNA Backbone Geometry in DNA-Protein Complexes: Experimental Data Analysis

We have analyzed and compared the available experimental data (PDB) on the backbone geometry of the DNA in solution (NMR), in crystals (X-rays), and in complexes with proteins (X-rays and cryo-electron microscopy). The deoxyribose (pseudorotational angle τ₀) and ε/ζ (BI-BII transition in phosphates) flexibilities are practically the same in the four samples. The α/γ mobility is minimal in crystalline DNA: on the histograms, there is one canonical and one noncanonical t/t peak. The α/γ mobility increases in DNA solutions (three more noncanonical peaks) and is maximal in DNA-protein complexes (another additional peak). On a large amount of data, we have confirmed that the three main degrees of freedom of the sugar-phosphate backbone are "orthogonal": changes in any of the angles τ₀, (ζ-ε), and (γ-α) occur, as a rule, at a constant (usually canonical) value of any other. In the DNA-protein complexes, none of the geometrical parameters commonly used to distinguish the A and B forms of DNA, except for Zp and its simpler analog Zp', show an unambiguous correlation with τ₀. Proteins, binding to DNA, in 59% of cases change the local shape of the helix up to the characteristic of the A-form without switching the deoxyribose conformation from south to north. However, we have found simple local characteristics of one nucleotide that correlate with the angles τ₀ and (ζ-ε). These are the angles C3'C1'N* and C4'C3'P(2), respectively. They are orthogonal in DNA-protein complexes exactly as the pair τ₀ and (ζ-ε). Most characteristics of DNA in complexes with proteins are the same in X-ray and in cryo-EM data, except for the histogram for the angle τ₀. We offer a possible explanation for this difference. We also discuss the artifacts on the ε/ζ histogram for DNA in solutions caused by the currently used NMR refinement protocols.

A dynamic view of DNA structure within the nucleosome: Biological implications

Most of eukaryotic cellular DNA is packed in nucleosome core particles (NCPs), in which the DNA (DNA_NCP) is wrapped around histones. The influence of this organization on the intrinsic local dynamics of DNA is largely unknown, in particular because capturing such information from experiments remains notoriously challenging. Given the importance of dynamical properties in DNA functions, we addressed this issue using CHARMM36 MD simulations of a nucleosome containing the NCP positioning 601 sequence and four related free dodecamers. Comparison between DNA_NCP and free DNA reveals a limited impact of the dense DNA-histone interface on correlated motions of dinucleotide constituents and on fluctuations of inter base pair parameters. A characteristic feature intimately associated with the DNA_NCP super-helical path is a set of structural periodicities that includes a marked alternation of regions enriched in backbone BI and BII conformers. This observation led to uncover a convincing correspondence between the sequence effect on BI/BII propensities in both DNA_NCP and free DNA, strengthening the idea that the histone preference for particular DNA sequences relies on those intrinsic structural properties. These results offer for the first time a detailed view of the DNA dynamical behavior within NCP. They show in particular that the DNA_NCP dynamics is substantial enough to preserve the ability to structurally adjust to external proteins, for instance remodelers. Also, fresh structural arguments highlight the relevance of relationships between DNA sequence and structural properties for NCP formation. Overall, our work offers a more rational framework to approach the functional, biological roles of NCP.

Structural features of DNA that determine RNA polymerase II core promoter

Background

The general structure and action of all eukaryotic and archaeal RNA polymerases machinery have an astonishing similarity despite the diversity of core promoter sequences in different species. The goal of our work is to find common characteristics of DNA region that define it as a promoter for the RNA polymerase II (Pol II).

Results

The profiles of a large number of physical and structural characteristics, averaged over representative sets of the Pol II minimal core promoters of the evolutionary divergent species from animals, plants and unicellular fungi were analysed. In addition to the characteristics defined at the base-pair steps, we, for the first time, use profiles of the ultrasonic cleavage and DNase I cleavage indexes, informative for internal properties of each complementary strand.

Conclusions

DNA of the core promoters of metazoans and Schizosaccharomyces pombe has similar structural organization. Its mechanical and 3D structural characteristics have singular properties at the positions of TATA-box. The minor groove is broadened and conformational motion is decreased in that region. Special characteristics of conformational behavior are revealed in metazoans at the region, which connects the end of TATA-box and the transcription start site (TSS). The intensities of conformational motions in the complementary strands are periodically changed in opposite phases. They are noticeable, best of all, in mammals. Such conformational features are lacking in the core promoters of S. pombe. The profiles of Saccharomyces cerevisiae core promoters significantly differ: their singular region is shifted down thus pointing to the uniqueness of their structural organization. Obtained results may be useful in genetic engineering for artificial modulation of the promoter strength.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3292-z) contains supplementary material, which is available to authorized users.

The intrinsic mechanics of B-DNA in solution characterized by NMR

Experimental characterization of the structural couplings in free B-DNA in solution has been elusive, because of subtle effects that are challenging to tackle. Here, the exploitation of the NMR measurements collected on four dodecamers containing a substantial set of dinucleotide sequences provides new, consistent correlations revealing the DNA intrinsic mechanics. The difference between two successive residual dipolar couplings (ΔRDCs) involving C6/8-H6/8, C3′-H3′ and C4′-H4′ vectors are correlated to the ³¹P chemical shifts (δP), which reflect the populations of the BI and BII backbone states. The δPs are also correlated to the internucleotide distances (D_inter) involving H6/8, H2′ and H2″ protons. Calculations of NMR quantities on high resolution X-ray structures and controlled models of DNA enable to interpret these couplings: the studied ΔRDCs depend mostly on roll, while D_inter are mainly sensitive to twist or slide. Overall, these relations demonstrate how δP measurements inform on key inter base parameters, in addition to probe the BI↔BII backbone equilibrium, and shed new light into coordinated motions of phosphate groups and bases in free B-DNA in solution. Inspection of the 5′ and 3′ ends of the dodecamers also supplies new information on the fraying events, otherwise neglected.

Simulations Meet Experiment to Reveal New Insights into DNA Intrinsic Mechanics

The accurate prediction of the structure and dynamics of DNA remains a major challenge in computational biology due to the dearth of precise experimental information on DNA free in solution and limitations in the DNA force-fields underpinning the simulations. A new generation of force-fields has been developed to better represent the sequence-dependent B-DNA intrinsic mechanics, in particular with respect to the BI ↔ BII backbone equilibrium, which is essential to understand the B-DNA properties. Here, the performance of MD simulations with the newly updated force-fields Parmbsc0εζOLI and CHARMM36 was tested against a large ensemble of recent NMR data collected on four DNA dodecamers involved in nucleosome positioning. We find impressive progress towards a coherent, realistic representation of B-DNA in solution, despite residual shortcomings. This improved representation allows new and deeper interpretation of the experimental observables, including regarding the behavior of facing phosphate groups in complementary dinucleotides, and their modulation by the sequence. It also provides the opportunity to extensively revisit and refine the coupling between backbone states and inter base pair parameters, which emerges as a common theme across all the complementary dinucleotides. In sum, the global agreement between simulations and experiment reveals new aspects of intrinsic DNA mechanics, a key component of DNA-protein recognition.

NMR proton chemical shift prediction of C·C mismatches in B-DNA

Accurate prediction of DNA chemical shifts facilitates resonance assignment and allows recognition of different conformational features. Based on the nearest neighbor model and base pair replacement approach, we have determined a set of triplet chemical shift values and correction factors for predicting the proton chemical shifts of B-DNA containing an internal C·C mismatch. Our results provide a reliable chemical shift prediction with an accuracy of 0.07 ppm for non-labile protons and 0.09 ppm for labile protons. In addition, we have also shown that the correction factors for C·C mismatches can be used interchangeably with those for T·T mismatches. As a result, we have generalized a set of correction factors for predicting the flanking residue chemical shifts of pyrimidine·pyrimidine mismatches.

Toward Improved Description of DNA Backbone: Revisiting Epsilon and Zeta Torsion Force Field Parameters

We present a refinement of the backbone torsion parameters ε and ζ of the Cornell et al. AMBER force field for DNA simulations. The new parameters, denoted as εζOL1, were derived from quantum-mechanical calculations with inclusion of conformation-dependent solvation effects according to the recently reported methodology (J. Chem. Theory Comput. 2012, 7(9), 2886-2902). The performance of the refined parameters was analyzed by means of extended molecular dynamics (MD) simulations for several representative systems. The results showed that the εζOL1 refinement improves the backbone description of B-DNA double helices and G-DNA stem. In B-DNA simulations, we observed an average increase of the helical twist and narrowing of the major groove, thus achieving better agreement with X-ray and solution NMR data. The balance between populations of BI and BII backbone substates was shifted towards the BII state, in better agreement with ensemble-refined solution experimental results. Furthermore, the refined parameters decreased the backbone RMS deviations in B-DNA MD simulations. In the antiparallel guanine quadruplex (G-DNA) the εζOL1 modification improved the description of non-canonical α/γ backbone substates, which were shown to be coupled to the ε/ζ torsion potential. Thus, the refinement is suggested as a possible alternative to the current ε/ζ torsion potential, which may enable more accurate modeling of nucleic acids. However, long-term testing is recommended before its routine application in DNA simulations.

Simulating DNA by molecular dynamics: aims, methods, and validation

The structure and dynamics of the B-DNA double helix involves subtle sequence-dependent effects which are decisive for its function, but difficult to characterize. These structural and dynamic effects can be addressed by simulations of DNA sequences in explicit solvent. Here, we present and discuss the state-of-art of B-DNA molecular dynamics simulations with the major force fields in use today. We explain why a critical analysis of the MD trajectories is required to assess their reliability, and estimate the value and limitations of these models. Overall, simulations of DNA bear great promise towards deciphering the structural and physical subtleties of this biopolymer, where much remains to be understood.

Structure, Stiffness and Substates of the Dickerson-Drew Dodecamer

The Dickerson-Drew dodecamer (DD) d-[CGCGAATTCGCG]₂ is a prototypic B-DNA molecule whose sequence-specific structure and dynamics have been investigated by many experimental and computational studies. Here, we present an analysis of DD properties based on extensive atomistic molecular dynamics (MD) simulations using different ionic conditions and water models. The 0.6-2.4-µs-long MD trajectories are compared to modern crystallographic and NMR data. In the simulations, the duplex ends can adopt an alternative base-pairing, which influences the oligomer structure. A clear relationship between the BI/BII backbone substates and the basepair step conformation has been identified, extending previous findings and exposing an interesting structural polymorphism in the helix. For a given end pairing, distributions of the basepair step coordinates can be decomposed into Gaussian-like components associated with the BI/BII backbone states. The nonlocal stiffness matrices for a rigid-base mechanical model of DD are reported for the first time, suggesting salient stiffness features of the central A-tract. The Riemann distance and Kullback-Leibler divergence are used for stiffness matrix comparison. The basic structural parameters converge very well within 300 ns, convergence of the BI/BII populations and stiffness matrices is less sharp. Our work presents new findings about the DD structural dynamics, mechanical properties, and the coupling between basepair and backbone configurations, including their statistical reliability. The results may also be useful for optimizing future force fields for DNA.

Probing sequence-specific DNA flexibility in a-tracts and pyrimidine-purine steps by nuclear magnetic resonance (13)C relaxation and molecular dynamics simulations

Sequence-specific DNA flexibility plays a key role in a variety of cellular interactions that are critical for gene packaging, expression, and regulation, yet few studies have experimentally explored the sequence dependence of DNA dynamics that occur on biologically relevant time scales. Here, we use nuclear magnetic resonance (NMR) carbon spin relaxation combined with molecular dynamics (MD) simulations to examine the picosecond to nanosecond dynamics in a variety of dinucleotide steps as well as in varying length homopolymeric A(n)·T(n) repeats (A(n)-tracts, where n = 2, 4, or 6) that exhibit unusual structural and mechanical properties. We extend the NMR spin relaxation time scale sensitivity deeper into the nanosecond regime by using glycerol and a longer DNA duplex to slow overall tumbling. Our studies reveal a structurally unique A-tract core (for n > 3) that is uniformly rigid, flanked by junction steps that show increasing sugar flexibility with A-tract length. High sugar mobility is observed at pyrimidine residues at the A-tract junctions, which is encoded at the dinucleotide level (CA, TG, and CG steps) and increases with A-tract length. The MD simulations reproduce many of these trends, particularly the overall rigidity of A-tract base and sugar sites, and suggest that the sugar-backbone dynamics could involve transitions in sugar pucker and phosphate backbone BI ↔ BII equilibria. Our results reinforce an emerging view that sequence-specific DNA flexibility can be imprinted in dynamics occurring deep within the nanosecond time regime that is difficult to characterize experimentally at the atomic level. Such large-amplitude sequence-dependent backbone fluctuations might flag the genome for specific DNA recognition.

Differential gene expression activity among species-specific polypyrimidine/polypurine motifs in mu opioid receptor gene promoters

The mu opioid receptor (MOR) is the principle molecular target of opioid analgesics. An appropriate understanding of MOR gene expression across species is critical for understanding its analgesic functions in humans. Here, we undertake a cross-species analysis of the polymorphic polypyrimidine/polypurine (PPy/u) motif, a key enhancer of MOR gene expression. The mouse PPy/u motif is highly homologous to those of rat (67%) and human (83%), but drives reporter gene expression tenfold and fivefold more effectively than those of rat and human, respectively. Circular dichroism profiles of PPy/u oligonucleotides from different species showed that they are primarily different in structure. Conformational studies of reporter plasmids using confocal Raman spectra, S1 nuclease and restriction enzymes demonstrated that the structural difference is the result of changes in the phosphodiester backbone. Furthermore, these conformational disparities produce differences in torsional stress, as shown by topoisomerase II relaxation and activation of different levels of gene expression under hypertonic conditions. This study demonstrates that homologous PPy/u motifs adopt unique species-specific conformations with different mechanisms and activities for gene expression. We further discuss how structural aspects of transcription regulatory elements, rather than the sequence itself, are significant when studying functional gene expression regulatory elements.

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.

Understanding the sequence-dependence of DNA groove dimensions: implications for DNA interactions

Background

The B-DNA major and minor groove dimensions are crucial for DNA-protein interactions. It has long been thought that the groove dimensions depend on the DNA sequence, however this relationship has remained elusive. Here, our aim is to elucidate how the DNA sequence intrinsically shapes the grooves.

Methodology/Principal Findings

The present study is based on the analysis of datasets of free and protein-bound DNA crystal structures, and from a compilation of NMR ³¹P chemical shifts measured on free DNA in solution on a broad range of representative sequences. The ³¹P chemical shifts can be interpreted in terms of the BI↔BII backbone conformations and dynamics. The grooves width and depth of free and protein-bound DNA are found to be clearly related to the BI/BII backbone conformational states. The DNA propensity to undergo BI↔BII backbone transitions is highly sequence-dependent and can be quantified at the dinucleotide level. This dual relationship, between DNA sequence and backbone behavior on one hand, and backbone behavior and groove dimensions on the other hand, allows to decipher the link between DNA sequence and groove dimensions. It also firmly establishes that proteins take advantage of the intrinsic DNA groove properties.

Conclusions/Significance

The study provides a general framework explaining how the DNA sequence shapes the groove dimensions in free and protein-bound DNA, with far-reaching implications for DNA-protein indirect readout in both specific and non specific interactions.

DNA structures from phosphate chemical shifts

For B-DNA, the strong linear correlation observed by nuclear magnetic resonance (NMR) between the (31)P chemical shifts (deltaP) and three recurrent internucleotide distances demonstrates the tight coupling between phosphate motions and helicoidal parameters. It allows to translate deltaP into distance restraints directly exploitable in structural refinement. It even provides a new method for refining DNA oligomers with restraints exclusively inferred from deltaP. Combined with molecular dynamics in explicit solvent, these restraints lead to a structural and dynamical view of the DNA as detailed as that obtained with conventional and more extensive restraints. Tests with the Jun-Fos oligomer show that this deltaP-based strategy can provide a simple and straightforward method to capture DNA properties in solution, from routine NMR experiments on unlabeled samples.

Intrinsic flexibility of B-DNA: the experimental TRX scale

B-DNA flexibility, crucial for DNA–protein recognition, is sequence dependent. Free DNA in solution would in principle be the best reference state to uncover the relation between base sequences and their intrinsic flexibility; however, this has long been hampered by a lack of suitable experimental data. We investigated this relationship by compiling and analyzing a large dataset of NMR ³¹P chemical shifts in solution. These measurements reflect the BI ↔ BII equilibrium in DNA, intimately correlated to helicoidal descriptors of the curvature, winding and groove dimensions. Comparing the ten complementary DNA dinucleotide steps indicates that some steps are much more flexible than others. This malleability is primarily controlled at the dinucleotide level, modulated by the tetranucleotide environment. Our analyses provide an experimental scale called TRX that quantifies the intrinsic flexibility of the ten dinucleotide steps in terms of Twist, Roll, and X-disp (base pair displacement). Applying the TRX scale to DNA sequences optimized for nucleosome formation reveals a 10 base-pair periodic alternation of stiff and flexible regions. Thus, DNA flexibility captured by the TRX scale is relevant to nucleosome formation, suggesting that this scale may be of general interest to better understand protein-DNA recognition.

Investigation of transcription factor Ndt80 affinity differences for wild type and mutant DNA: a molecular dynamics study

Molecular dynamics simulations and free energy calculations have been performed on the transcription factor Ndt80 either in complex with the native DNA sequence or with a mutant DNA with a switched central base pair, C5-G5' to G5-C5'. This mutant has been shown to have a 100-fold decrease in binding affinity of Ndt80, and in this study we explain this both structurally and energetically. The major interactions between the protein and the DNA were maintained in the simulations, apart from around the mutation site. The crystal structure of the Ndt80-DNA complex revealed that R177 makes a base specific bidentate interaction with G5' which is part of a conserved 5'-YpG-3' step. In the simulation with the mutant DNA, the side chain of R177 changes conformation and makes three new stable hydrogen bonds to the DNA backbone. This in turn induces a conformational change in the DNA backbone of the T6'-G5' step from the unusual BII state to the canonical BI state. The affinity difference for the protein-DNA complex with the native DNA compared with the mutant DNA is only about 3 kcal/mol. The free energy calculations of the base pair switch indicated a larger difference than what was found experimentally, about 7.7 kcal/mol, but this is explained in structural terms using the 10 ns simulations of the solvated complexes and the rearrangement of the R177 side chain.

31P NMR investigation of backbone dynamics in DNA binding sites

The backbone conformation of DNA plays an important role in the indirect readout mechanisms for protein--DNA recognition events. Thus, investigating the backbone dynamics of each step in DNA binding sequences provides useful information necessary for the characterization of these interactions. Here, we use 31P dynamic NMR to characterize the backbone conformation and dynamics in the Dickerson dodecamer, a sequence containing the EcoRI binding site, and confirm solid-state 2H NMR results showing that the C3pG4 and C9pG10 steps experience unique dynamics and that these dynamics are quenched upon cytosine methylation. In addition, we show that cytosine methylation affects the conformation and dynamics of neighboring nucleotide steps, but this effect is localized to only near neighbors and base-pairing partners. Last, we have been able to characterize the percent BII in each backbone step and illustrate that the C3pG4 and C9pG10 favor the noncanonical BII conformation, even at low temperatures. Our results demonstrate that 31P dynamic NMR provides a robust and efficient method for characterizing the backbone dynamics in DNA. This allows simple, rapid determination of sequence-dependent dynamical information, providing a useful method for studying trends in protein-DNA recognition events.

Cytosine ribose flexibility in DNA: a combined NMR 13C spin relaxation and molecular dynamics simulation study

Using (13)C spin relaxation NMR in combination with molecular dynamic (MD) simulations, we characterized internal motions within double-stranded DNA on the pico- to nano-second time scale. We found that the C-H vectors in all cytosine ribose moieties within the Dickerson-Drew dodecamer (5'-CGCGAATTCGCG-3') are subject to high amplitude motions, while the other nucleotides are essentially rigid. MD simulations showed that repuckering is a likely motional model for the cytosine ribose moiety. Repuckering occurs with a time constant of around 100 ps. Knowledge of DNA dynamics will contribute to our understanding of the recognition specificity of DNA-binding proteins such as cytosine methyltransferase.

Quantification of DNA BI/BII backbone states in solution. Implications for DNA overall structure and recognition

The backbone states of B-DNA influence its helical parameters, groove dimensions, and overall curvature. Therefore, detection and fine characterization of these conformational states are desirable. Using routine NMR experiments on a nonlabeled B-DNA oligomer and analyzing high-resolution X-ray structures, we investigated the relationship between interproton distances and backbone conformational states. The three H2'i-H6/8i+1, H2' 'i-H6/8i+1, and H6/8i-H6/8i+1 sequential distances were found cross-correlated and linearly coupled to epsilon-zeta values in X-ray structures and 31P chemical shifts (deltaP) in NMR that reflect the interconversion between the backbone BI (epsilon-zeta < 0 degrees ) and BII (epsilon-zeta > 0 degrees) states. These relationships provide a detailed check of the NMR data consistency and the possibility to extend the set of restraints for structural refinement through various extrapolations. Furthermore, they allow translation of deltaP in terms of BI/BII ratios. Also, comparison of many published deltaP in solution to crystal data shows that the impact of sequence on the BI/BII propensities is similar in both environments and is therefore an intrinsic and general property of B-DNA. This quantification of the populations of BI and BII is of general interest because these sequence-dependent backbone states act on DNA overall structure, a key feature for DNA-protein-specific recognition.

Refinement of the AMBER force field for nucleic acids: improving the description of alpha/gamma conformers

We present here the parmbsc0 force field, a refinement of the AMBER parm99 force field, where emphasis has been made on the correct representation of the alpha/gamma concerted rotation in nucleic acids (NAs). The modified force field corrects overpopulations of the alpha/gamma = (g+,t) backbone that were seen in long (more than 10 ns) simulations with previous AMBER parameter sets (parm94-99). The force field has been derived by fitting to high-level quantum mechanical data and verified by comparison with very high-level quantum mechanical calculations and by a very extensive comparison between simulations and experimental data. The set of validation simulations includes two of the longest trajectories published to date for the DNA duplex (200 ns each) and the largest variety of NA structures studied to date (15 different NA families and 97 individual structures). The total simulation time used to validate the force field includes near 1 mus of state-of-the-art molecular dynamics simulations in aqueous solution.

Conformations and Dynamics of the Phosphodiester Backbone of a DNA Fragment That Bears a Strong Topoisomerase II Cleavage Site

The dynamics of the DNA phosphodiester backbone conformations have been studied for a strong topoisomerase II cleavage site (site 22) using molecular dynamics simulations in explicit water and in the presence of sodium ions. We investigated the backbone motions and more particularly the BI/BII transitions involving the epsilon and zeta angles. The consensus cleavage site is adjacent to the phosphate which shows the most important phosphodiester backbone flexibility in the sequence. We infer that these latter properties could be responsible for the preferential cleavage at this site possibly through the perturbation of the cleavage/ligation activities of the topoisomerase II. More generally, the steps pur-pur and pyr-pur are those presenting the highest BII contents. Relations are observed between the backbone phosphodiester BI/BII transitions and the flexibility of the deoxyribose sugar and the helical parameters such as roll. The roll is sequence dependent when the related phosphate is in the BI form, whereas this appears not to be true when it is in the BII form. The BI/BII transitions are associated with water migration, and new relations are observed with counterions. Indeed, it is observed that a strong coupling exists between the BII form and the presence of sodium ions near the adjacent sugar deoxyribose. The presence of sodium ions in the O4' surroundings or their binding could assist the BI to BII transition by furnishing energy. The implications of these new findings and, namely, their importance in the context of the sequence-dependent behavior of BI/BII transitions will be investigated in future studies.

Sequence preference for BI/BII conformations in DNA: MD and crystal structure data analysis

Deciphering sequence information from sugar-phosphate backbone is finely tuned through the conformational substates of DNA. BII conformation, one of the conformational substates of B-DNA, is known to play a key role in DNA-protein recognition. BI and BII are identified by the epsilon-zeta difference, which is negative in BI and positive in BII. Our analysis of MD and crystal structures shows that BII conformation is sequence specific and dinucleotides GC, CG, CA, TG, TA show high preference to take up BII conformation, while TT, TC, CT, CC dinucleotides rarely take up this conformation. Significant changes were observed in the dinucleotide parameters viz. twist, roll, and slide for the steps having BII conformation. Interestingly, the magnitude of variation in the dinucleotide parameters is seen to depend mainly on two factors, the magnitude of epsilon-zeta difference and the presence or absence of BII conformation in the second strand, across the WC base-paired dinucleotide step. Based on these two factors, the conformational substate of a dinucleotide step can be further classified as BI.BI (BI conformation in both strands), BI.BII (BI conformation in one strand and BII conformation in the other), and BII.BII (BII conformation in both strands). The occurrence of BII in both strands was found to be quite rare and thus, it can be concluded that BI.BI and BI.BII hybrid steps are more favorable than a BII.BII step. In conformity with the sequence preference seen for dinucleotides in each strand, BII.BII combination of backbone conformation was observed only for GC, CG, CA, and TG containing dinucleotide steps. We further classified BII.BII step as strong BII and weak BII depending on the magnitude of the average epsilon-zeta difference. The dinucleotide steps which belong to the category of strong BII, have large twist, high positive slide and negative roll values, while those in the weak BII group have roll, twist, and slide values similar to that of hybrid BI.BII steps. This conformational property could be contributing to the groove opening/closing and thus can modulate protein-DNA interaction.

Discrimination between BI and BII conformational substates of B-DNA based on sugar-base interproton distances

Molecular dynamics (MD) simulations of four water-solvated DNA duplexes were used to generate a database of approximately 27000 dinucleotide conformations. Analyzing this database, we investigated the relationship between so-called BI-BII transitions and short-range interproton distances. Four H-H distances were found particularly sensitive to BI-BII transitions: internucleotide H1'(n)-H68(n+1), H2'(n)-H68(n+1), and H2"(n)-H68(n+1), and intranucleotide H2"(n)-H68(n). Determination of these distances using classical NOESY spectroscopy can thus provide valuable indications on the existence of BII substates, complementing the existing method based on (31)P chemical shifts and (31)P-(1)H spin-spin coupling constants.

DNA fine structure and dynamics in crystals and in solution: The impact of BI/BII backbone conformations

Sugar-phosphate backbone conformations are an important structural element for a complete understanding of specific recognition in nucleic acid-protein interactions. They can be involved both in early stages of target discrimination and in structural adaptation upon binding. In the first part of this study, we have analyzed high-resolution structures of double-stranded B-DNA either isolated or bound to proteins, and explored the impact of both the standard BI and the unusual BII phosphate backbone conformations on neighboring sugar puckers and on selected helical parameters. Correlations are found to be similar for free and bound DNA, and in both categories, the possible facing backbone conformations (BI.BI, BI.BII, and BII.BII) define well-characterized substates in the B-DNA conformational space. Notably, BII.BII steps are characterized by specific, and sequence-independent, structural effects involving reduced standard deviations for almost all conformational parameters. In the second part of this work, we analyze four 10 ns molecular dynamics simulations in explicit solvent on the DNA targets of NF-kappaB and bovine papillomavirus E2 proteins, highlighting the multiplicity of backbone dynamical behavior. These results show sequence effects on the percentages of BI and BII conformers, the preferential state of facing backbones, the occurrence of coupled transitions. The backbone states can consequently be seen as a mechanism for transmitting information from the bases to the phosphate groups and thus for modulating the overall structural properties of the target DNA.

NMR structure of an alpha-L-LNA:RNA hybrid: structural implications for RNase H recognition

Alpha-L-LNA (alpha-L-ribo configured locked nucleic acid) is a nucleotide analogue that raises the thermostability of nucleic acid duplexes by up to approximately 4 degrees C per inclusion. We have determined the NMR structure of a nonamer alpha-L-LNA:RNA hybrid with three alpha-L-LNA modifications. The geometry of this hybrid is intermediate between A- and B-type, all nucleobases partake in Watson-Crick base pairing and base stacking, and the global structure is very similar to that of the corresponding unmodified hybrid. The sugar-phosphate backbone is rearranged in the vicinity of the modified nucleotides. As a consequence, the phosphate groups following the modified nucleotides are rotated into the minor groove. It is interesting that the alpha-L-LNA:RNA hybrid, which has an elevation in melting temperature of 17 degrees C relative to the corresponding DNA:RNA hybrid, retains the global structure of this hybrid. To our knowledge, this is the first example of such a substantial increase in melting temperature of a nucleic acid analogue that does not act as an N-type (RNA) mimic. alpha-L-LNA:RNA hybrids are recognised by RNase H with subsequent cleavage of the RNA strand, albeit with slow rates. We attempt to rationalise this impaired enzyme activity from the rearrangement of the sugar-phosphate backbone of the alpha-L-LNA:RNA hybrid.

Conformational characteristics and correlations in crystal structures of nucleic acid oligonucleotides: evidence for sub-states

Sugar phosphate backbone conformations are a structural element inextricably involved in a complete understanding of specific recognition nucleic acid ligand interactions, from early stage discrimination of the correct target to complexation per se, including any structural adaptation on binding. The collective results of high resolution DNA, RNA and protein/DNA crystal structures provide an opportunity for an improved and enhanced statistical analysis of standard and unusual sugar-phosphate backbone conformations together with corresponding dinucleotide sequence effects as a basis for further exploration of conformational effects on binding. In this study, we have analyzed the conformations of all relevant crystal structures in the nucleic acids data base, determined the frequency distribution of all possible epsilon, zeta, alpha, beta and gamma backbone angle arrangements within four nucleic acid categories (A-RNA and A-DNA, free and bound B-DNA) and explored the relationships between backbone angles, sugar puckers and selected helical parameters. The trends in the correlations are found to be similar regardless of the nucleic acid category. It is interesting that specific structural effects exhibited by the different unusual backbone sub-states are in some cases contravariant. Certain alpha/gamma changes are accompanied by C3' endo (north) sugars, small twist angles and positive values of base pair roll, and favor a displacement of nucleotide bases towards the minor groove compared to that of canonical B form structures. Unusual epsilon/zeta combinations occur with C2' (south) sugars, high twist angles, negative values of base pair roll, and base displacements towards the major groove. Furthermore, any unusual backbone correlates with a reduced dispersion of equilibrium structural parameters of the whole double helix, as evidenced by the reduced standard deviations of almost all conformational parameters. Finally, a strong sequence effect is displayed in the free oligomers, but reduced somewhat in the ligand bound forms. The most variable steps are GpA and CpA, and, to a lesser extent, their partners TpC and TpG. The results provide a basis for considering if the variable and non-variable steps within a biological active sequence precisely determine morphological structural features as the curvature direction, the groove depth, and the accessibility of base pair for non covalent associations.

Operator recognition by the phage 434 cI repressor: MD simulations of free and bound 50-bp DNA reveal important differences between the OR1 and OR2 sites

Using molecular dynamics simulations in explicit solvent, we investigated the behavior of a 50-bp DNA sequence containing the 434 bacteriophage operators OR1 and OR2 separated by an 8-bp spacer. Two simulations of 1 ns each were carried out, with DNA alone and with DNA complexed to dimers of the R1-69 DNA binding domain of the phage 434 cI repressor protein at the OR1 and OR2 sites. Strong correlations among average structural parameters are observed between our simulations and available experimental data for the bound OR1/OR2 subsites. In the free state, some differences appear between the three relevant fragments (OR1, the spacer, and OR2). Unbound OR1 exhibits a large, shallow major groove into which the base atoms protrude and is also bent toward the major groove. This structure is maintained because structural fluctuations are weak. Unbound OR2 resembles canonical B-DNA although the structural parameters show greater fluctuations, essentially due to a malleable step (the innermost CpA/TpG), absent in OR1. Complexation with the proteins slightly alters the base positions but strongly modifies the sugar and backbone motions. The most crucial repressor effects are changes in the flexibility of the OR1/OR2 sites. Structural fluctuations are enhanced for OR1, conferring a favorable energetic contribution to the OR1 binding, whereas they are reduced for OR2. Therefore, both structural and dynamic properties of DNA suggest OR1 is the most attractive site for the repressor, which may explain the different binding association constants observed for the OR1 and OR2 sites. Finally, we also investigated the impact of the protein on the DNA backbone dynamics and find that direct or indirect interactions facilitate the DNA structural variations required for achieving complementarity with the protein.

Alpha/gamma transitions in the B-DNA backbone

In the crystal structures of protein complexes with B-DNA, alpha and gamma DNA backbone torsion angles often exhibit non-canonical values. It is not known if these alternative backbone conformations are easily accessible in solution and can contribute to the specific recognition of DNA by proteins. We have analysed the coupled transition of the alpha and gamma torsion angles within the central GpC step of a B-DNA dodecamer by computer simulations. Five stable or metastable non-canonical alpha/gamma sub-states are found. The most favourable pathway from the canonical alpha/gamma structure to any unusual form involves a counter-rotation of alpha and gamma, via the trans conformation. However, the corresponding free energy indicates that spontaneous flipping of the torsions is improbable in free B-DNA. This is supported by an analysis of the available high resolution crystallographic structures showing that unusual alpha/gamma states are only encountered in B-DNA complexed to proteins. An analysis of the structural consequences of alpha/gamma transitions shows that the non-canonical backbone geometry influences essentially the roll and twist values and reduces the equilibrium dispersion of structural parameters. Our results support the hypothesis that unusual alpha/gamma backbones arise during protein-DNA complexation, assisting the fine structural adjustments between the two partners and playing a role in the overall complexation free energy.

Sequence-dependent folding of DNA three-way junctions

Three-way DNA junctions can adopt several different conformers, which differ in the coaxial stacking of the arms. These structural variants are often dominated by one conformer, which is determined by the DNA sequence. In this study we have compared several three-way DNA junctions in order to assess how the arrangement of bases around the branch point affects the conformer distribution. The results show that rearranging the different arms, while retaining their base sequences, can affect the conformer distribution. In some instances this generates a structure that appears to contain parallel coaxially stacked helices rather than the usual anti-parallel arrangement. Although the conformer equilibrium can be affected by the order of purines and pyrimidines around the branch point, this is not sufficient to predict the conformer distribution. We find that the folding of three-way junctions can be separated into two groups of dinucleotide steps. These two groups show distinctive stacking properties in B-DNA, suggesting there is a correlation between B-DNA stacking and coaxial stacking in DNA junctions.

Structural differences in the NOE-derived structure of G-T mismatched DNA relative to normal DNA are correlated with differences in (13)C relaxation-based internal dynamics

Detailed description of the characteristics of mismatched DNA that are distinct from normal DNA is vital to the understanding of how mismatch repair proteins are able to recognize and repair these DNA lesions. To this end, we have used nuclear Overhauser effect spectroscopy (NOESY)-based distance restraints and (13)C relaxation measurements to solve the solution structures and measure some of the internal dynamics of the G-T mismatched DNA oligomer d(CCATGCGTGG)(2) (GT) and its parent DNA sequence d(CCACGCGTGG)(2) (GC). In GT, the mismatched G7 is structurally perturbed much more than the mismatched T4 relative to their corresponding bases in GC. The degree of G7 displacement differs from previous high-resolution structures of G-T mismatch-containing B-DNA, suggesting a dependence of G-T mismatch-induced structural perturbation on sequence context. The internal dynamics of GC and GT differ on multiple timescales. The mismatched G7 of GT contains spins that decrease significantly in order in GT compared to GC, while spins in C6, T8, and A3 have significantly higher order in GT compared to GC. Linear correlations between helical parameters of GC and GT and the order of C-1' and aromatic methine carbon atoms relate differences in internal dynamics to the structures quantitatively. The dynamic differences between the normal and mismatched DNA signify changes in local flexibility that may be exploited by the mismatch repair system to bind mismatched DNA preferentially while ignoring normal DNA.