Solution structure of the CpG containing d(CTTCGAAG)2 oligonucleotide: NMR data and energy calculations are compatible with a BI/BII equilibrium at CpG

Mechanical codes of chemical-scale specificity in DNA motifs

In gene transcription, certain sequences of double-stranded (ds)DNA play a vital role in nucleosome positioning and expression initiation. That dsDNA is deformed to various extents in these processes leads us to ask: Could the genomic DNA also have sequence specificity in its chemical-scale mechanical properties? We approach this question using statistical machine learning to determine the rigidity between DNA chemical moieties. What emerges for the polyA, polyG, TpA, and CpG sequences studied here is a unique trigram that contains the quantitative mechanical strengths between bases and along the backbone. In a way, such a sequence-dependent trigram could be viewed as a DNA mechanical code. Interestingly, we discover a compensatory competition between the axial base-stacking interaction and the transverse base-pairing interaction, and such a reciprocal relationship constitutes the most discriminating feature of the mechanical code. Our results also provide chemical-scale understanding for experimental observables. For example, the long polyA persistence length is shown to have strong base stacking while its complement (polyAc) exhibits high backbone rigidity. The mechanical code concept enables a direct reading of the physical interactions encoded in the sequence which, with further development, is expected to shed new light on DNA allostery and DNA-binding drugs.

First Structure of a Designed Minor Groove Binding Heterocyclic Cation that Specifically Recognizes Mixed DNA Base Pair Sequences

The high-resolution NMR structure of the first heterocyclic, non-amide, organic cation that strongly and selectively recognizes mixed AT/GC bp (bp=base pair) sequences of DNA in a 1:1 complex is described. Compound designs of this type provide essential methods for control of functional, non-genomic DNA sequences and have broad cell uptake capability, based on studies from animals to humans. The high-resolution structural studies described in this report are essential for understanding the molecular basis for the sequence-specific binding as well as for new ideas for additional compound designs for sequence-specific recognition. The molecular features, in this report, explain the mechanism of recognition of both A⋅T and G⋅C bps and are an interesting molecular recognition story. Examination of the experimental structure and the NMR restrained molecular dynamics model suggests that recognition of the G⋅C base pair involves two specific H-bonds. The structure illustrates a wealth of information on different DNA interactions and illustrates an interfacial water molecule that is a key component of the complex.

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.

Automatic workflow for the classification of local DNA conformations

Background

A growing number of crystal and NMR structures reveals a considerable structural polymorphism of DNA architecture going well beyond the usual image of a double helical molecule. DNA is highly variable with dinucleotide steps exhibiting a substantial flexibility in a sequence-dependent manner. An analysis of the conformational space of the DNA backbone and the enhancement of our understanding of the conformational dependencies in DNA are therefore important for full comprehension of DNA structural polymorphism.

Results

A detailed classification of local DNA conformations based on the technique of Fourier averaging was published in our previous work. However, this procedure requires a considerable amount of manual work. To overcome this limitation we developed an automatic classification method consisting of the combination of supervised and unsupervised approaches. A proposed workflow is composed of k-NN method followed by a non-hierarchical single-pass clustering algorithm. We applied this workflow to analyze 816 X-ray and 664 NMR DNA structures released till February 2013. We identified and annotated six new conformers, and we assigned four of these conformers to two structurally important DNA families: guanine quadruplexes and Holliday (four-way) junctions. We also compared populations of the assigned conformers in the dataset of X-ray and NMR structures.

Conclusions

In the present work we developed a machine learning workflow for the automatic classification of dinucleotide conformations. Dinucleotides with unassigned conformations can be either classified into one of already known 24 classes or they can be flagged as unclassifiable. The proposed machine learning workflow permits identification of new classes among so far unclassifiable data, and we identified and annotated six new conformations in the X-ray structures released since our previous analysis. The results illustrate the utility of machine learning approaches in the classification of local DNA conformations.

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.

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.

DNA conformations and their sequence preferences

The geometry of the phosphodiester backbone was analyzed for 7739 dinucleotides from 447 selected crystal structures of naked and complexed DNA. Ten torsion angles of a near-dinucleotide unit have been studied by combining Fourier averaging and clustering. Besides the known variants of the A-, B- and Z-DNA forms, we have also identified combined A + B backbone-deformed conformers, e.g. with alpha/gamma switches, and a few conformers with a syn orientation of bases occurring e.g. in G-quadruplex structures. A plethora of A- and B-like conformers show a close relationship between the A- and B-form double helices. A comparison of the populations of the conformers occurring in naked and complexed DNA has revealed a significant broadening of the DNA conformational space in the complexes, but the conformers still remain within the limits defined by the A- and B- forms. Possible sequence preferences, important for sequence-dependent recognition, have been assessed for the main A and B conformers by means of statistical goodness-of-fit tests. The structural properties of the backbone in quadruplexes, junctions and histone-core particles are discussed in further detail.

Genome-scale relationships between cytosine methylation and dinucleotide abundances in animals

In mammalian genomes CpGs occur at one-fifth their expected frequency. This is accepted as resulting from cytosine methylation and deamination of 5-methylcytosine leading to TpG and CpA dinucleotides. The corollary that a CpG deficit should correlate with TpG excess has not hitherto been systematically tested at a genomic level. I analyzed genome sequences (human, chimpanzee, mouse, pufferfish, zebrafish, sea squirt, fruitfly, mosquito, and nematode) to do this and generally to assess the hypothesis that CpG deficit, TpG excess, and other data are accountable in terms of 5-methylcytosine mutation. In all methylated genomes local CpG deficit decreases with higher G + C content. Local TpG surplus, while positively associated with G + C level in mammalian genomes but negatively associated with G + C in nonmammalian methylated genomes, is always explicable in terms of the CpG trend under the methylation model. Covariance of dinucleotide abundances with G + C demonstrates that correlation analyses should control for G + C. Doing this reveals a strong negative correlation between local CpG and TpG abundances in methylated genomes, in accord with the methylation hypothesis. CpG deficit also correlates with CpT excess in mammals, which may reflect enhanced cytosine mutation in the context 5'-YCG-3'. Analyses with repeat-masked sequences show that the results are not attributable to repetitive elements.

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.

Solution structure of a purine rich hexaloop hairpin belonging to PGY/MDR1 mRNA and targeted by antisense oligonucleotides

A preferential target of antisense oligonucleotides directed against human PGY/MDR1 mRNA is a hairpin containing a stem with a G*U wobble pair, capped by the purine-rich 5'r(GGGAUG)3' hexaloop. This hairpin is studied by multidimensional NMR and restrained molecular dynamics, with special emphasis on the conformation of south sugars and non-standard phosphate linkages evidenced in both the stem and the loop. The hairpin is found to be highly structured. The G*U wobble pair, a strong counterion binding site, displays structural particularities that are characteristic of this type of mismatch. The upper part of the stem undergoes distortions that optimize its interactions with the beginning of the loop. The loop adopts a new fold in which the single-stranded GGGA purine tract is structured in A-like conformation stacked in continuity of the stem and displays an extensive hydrogen bonding surface for recognition. The remarkable hairpin stability results from classical inter- and intra-strand interactions reinforced by numerous hydrogen bonds involving unusual backbone conformations and ribose 2'-hydroxyl groups. Overall, this work emphasizes numerous features that account for the well-ordered structure of the whole hairpin and highlights the loop properties that facilitate interaction with antisense oligonucleotides.

M.TaqI facilitates the base flipping via an unusual DNA backbone conformation

MD simulations have been carried out to understand the dynamical behavior of the DNA substrate of the Thermus aquaticus DNA methyltransferase (M.TaqI) in the methylation process at N6 of adenine. As starting structures, an x-ray structure of M.TaqI in complex with DNA and cofactor analogue (PDB code: 1G 38) and free decamer d(GTTCGATGTC)(2) were taken. The x-ray structure shows two consecutive BII substates that are not observed in the free decamer. These consecutive BII substates are also observed during our simulation. Additionally, their facing backbones adopt the same conformations. These double facing BII substates are stable during the last 9 ns of the trajectories and result in a stretched DNA structure. On the other hand, protein-DNA contacts on 5' and 3' phosphodiester groups of the partner thymine of flipped adenine have changed. The sugar and phosphate parts of thymine have moved further into the empty space left by the flipping base without the influence of protein. Furthermore, readily high populated BII substates at the GpA step of palindromic tetrad TCGA rather than CpG step are observed in the free decamer. On the contrary, the BI substate at the GpA step is observed on the flipped adenine strand. A restrained MD simulation, reproducing the BI/BII pattern in the complex, demonstrated the influence of the unusual backbone conformation on the dynamical behavior of the target base. This finding along with the increased nearby interstrand phosphate distance is supportive to the N6-methylation mechanism.

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.

The role of DNA structure and dynamics in the recognition of bovine papillomavirus E2 protein target sequences

The papillomavirus E2 transcription and replication factors bind to the DNA consensus ACCGN(4)CGGT sequence (E2-BS), through both direct and indirect readout mechanisms. The two symmetric half-sites ACCG.CGGT are highly conserved in the genomes and are hydrogen bound with E2. Although E2 does not contact the N4 spacer, the affinities are modulated by the base composition of this DNA part. Nevertheless, the origin of either the global recognition mechanism or the spacer effect remains unclear, particularly in the case of the bovine papillomavirus type 1 E2 (BPV-1-E2) system, used as model to study the papillomaviruses. We present, herein, studies carried out on oligomers differently recognized by the BPV-1-E2 protein and based on molecular dynamic simulations including counterions and water. The sequences contain the conserved half-sites but three different spacers (CCAT, ACGT and AAAC), resulting in very high, high and low affinity targets for BPV-1-E2. In order to estimate how much the free DNAs resemble the bound conformations, comparisons are made with two DNAs extracted from E2-BS-BPV-1 crystallographic complexes, representative of high and moderate affinity structures. The analysis of 15 ns trajectories reveals that the ACCG/CGGT half-sites, whatever the spacer, have the same behavior and adopt average stable base-pair parameters very close to those of the bound conformations. In contrast, the three different free spacers strongly differ in their BI <--> BII backbone dynamics. The low affinity AAAC spacer exhibits stable BI backbone conformations, the high affinity ACGT spacer is characterized by a dramatic instability of the CpG phosphate groups, and the CpA and GpG backbones in the very high affinity CCAT.ATGG spacer are trapped in BII conformations. All resemble more of the moderate affinity complex DNA than the high affinity one. Nevertheless, the particular behavior of the CCAT and ACGT backbones allows the emergence of BII-rich spacers, a configuration reproducing both local and global helical features of the bound DNA conformation of the high affinity complex and favoring the minor groove curvature required in the complex. In particular, the CCAT-containing site spends almost half of the time in this form that well mimics the bound one. Thus, we propose that the E2 protein could take advantage of the invariant favorable structures of the half-sites to form a pre-complex, but would require a specific spacer intrinsic malleability to lock the interaction. Finally, the backbone conformational states, by their ability to translate information coded in the sequence into structural properties, provide insight into the mechanisms that contribute to fine binding site selection and specific nucleic acid ligand recognition.

Molecular dynamics simulations of DNA curvature and flexibility: helix phasing and premelting

Recent studies of DNA axis curvature and flexibility based on molecular dynamics (MD) simulations on DNA are reviewed. The MD simulations are on DNA sequences up to 25 base pairs in length, including explicit consideration of counterions and waters in the computational model. MD studies are described for ApA steps, A-tracts, for sequences of A-tracts with helix phasing. In MD modeling, ApA steps and A-tracts in aqueous solution are essentially straight, relatively rigid, and exhibit the characteristic features associated with the B'-form of DNA. The results of MD modeling of A-tract oligonucleotides are validated by close accord with corresponding crystal structure results and nuclear magnetic resonance (NMR) nuclear Overhauser effect (NOE) and residual dipolar coupling (RDC) structures of d(CGCGAATTCGCG) and d(GGCAAAAAACGG). MD simulation successfully accounts for enhanced axis curvature in a set of three sequences with phased A-tracts studied to date. The primary origin of the axis curvature in the MD model is found at those pyrimidine/purine YpR "flexible hinge points" in a high roll, open hinge conformational substate. In the MD model of axis curvature in a DNA sequence with both phased A-tracts and YpR steps, the A-tracts appear to act as positioning elements that make the helix phasing more precise, and key YpR steps in the open hinge state serve as curvature elements. Our simulations on a phased A-tract sequence as a function of temperature show that the MD simulations exhibit a premelting transition in close accord with experiment, and predict that the mechanism involves a B'-to-B transition within A-tracts coupled with the prediction of a transition in key YpR steps from the high roll, open hinge, to a low roll, closed hinge substate. Diverse experimental observations on DNA curvature phenomena are examined in light of the MD model with no serious discrepancies. The collected MD results provide independent support for the "non-A-tract model" of DNA curvature. The "junction model" is indicated to be a special case of the non-A-tract model when there is a Y base at the 5' end of an A-tract. In accord with crystallography, the "ApA wedge model" is not supported by MD.

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.

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.

Structural effects of cytosine methylation on DNA sugar pucker studied by FTIR

This FTIR investigation concerns structural consequences of 5-methylation of cytosine in a DNA decamer in solution. Methylation of DNA is an important functional signal in transcription, but its effect on DNA structure is variable and not fully understood. Here, single and multiple 5-methylcytosine substitutions are introduced into the self-complementary sequence d(CCGGCGCCGG)(2). No major structural effect of methylation on the DNA duplex in solution is seen in the IR spectra: The overall B-form character of the backbone and S-type of sugar puckering are maintained in all the studied sequences, in agreement with previous literature. However, certain significant effects are detected in the IR regions sensitive to sugar pucker and glycosidic torsional angle. A single or multiple 5-methylcytosine substitution in d(CCGGCGCCGG)(2) leads to a doublet splitting of the S-type 840-820 cm(-1) sugar conformational band. The results suggest the coexistence of two different major sugar puckers within the S-conformational family, with an increased relative contribution of the C2'-endo type of sugar in the methylated sequences. In addition, a partial or full downshift of the guanosine/anti marker band at 1,375 cm(-1) in the methylated sequences reflects a change in the value of the dihedral angle chi of guanosine upon methylation. The IR spectra are interpreted in terms of localized transitions between the BI and BII subconformational states of the B-DNA backbone caused by the methylation. An increased amount of the BII subconformer in the methylated sequences should give rise to a structurally more rigid conformation, in agreement with earlier observations on DNA backbone dynamics and bending flexibility in methylated DNA.

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

The conformational exchange of the phosphate and deoxyribose groups of the DNA oligomers d(GCGTACGC)(2) and d(CGCTAGCG)(2) have been investigated using a combination of homonuclear and heteronuclear NMR techniques. Two-state exchange between phosphate B(I) and B(II) conformations and deoxyribose N and S conformations was expressed as percent population of the major conformer, %B(I) or %S. Sequence context-dependent variations in %B(I) and %S were observed. The positions of the phosphate and deoxyribose equilibria provide a quantitative measure of the ps to ns timescale dynamic exchange processes in the DNA backbone. Linear correlations between %B(I), %S, and previously calculated model free (13)C order parameters (S(2)) were observed. The %B(I) of the phosphates were found to be correlated to the S(2) of the flanking C3' and C4' atoms. The %B(I) was also found to be correlated with the %S and C1' S(2) of the deoxyribose ring 5' of the phosphates. The %B(I) of opposing phosphates is correlated, while the %B(I) of sequential phosphates is anti-correlated. These correlations suggest that conformational exchange processes in DNA are coupled to each other and are modulated by DNA base sequence, which may have important implications for DNA-protein interactions.

Impact of CpG methylation on structure, dynamics and solvation of cAMP DNA responsive element

Methylation of CpG motifs in DNA is involved in the control of gene expression and in several other epigenic effects. It suppresses also the immuno-stimulation properties of bacterial or viral DNAs that contain CPGS: However, effects of methylation on the DNA structure and dynamics are not clear. Here we carried out a 10 ns MD simulation, confronted to an NMR analysis, of a hexadecanucleotide with the cAMP responsive element (CRE) DNA methylated at its center: d(GAGATGAmCGTCATCTC)(2) (CREmet). Methylation does not introduce significant structure modification but reduces the dynamics. Molecular mechanics and generalized Born solvation energy calculations showed that the stiffness of CREmet arises from both a restriction of the conformational space by the bulky methyl groups and a folding of DNA around the hydrophobic methyls. The latter effect is favored when the GpA steps belonging to the TGA binding half-sites adopt the BII conformation. The inability of the methylated DNAs to interact with their protein partners-either transcription factors for gene regulation or a Toll-like receptor for immunostimulation-could result from both the obstacle created by methyls, preventing crucial interactions, and the loss of DNA flexibility, reducing its adaptability. Results are discussed in the light of NMR and crystallographic data.

Solution conformation of an RNA--DNA hybrid duplex containing a pyrimidine RNA strand and a purine DNA strand

RNA--DNA hybrid duplexes are involved in transcription, replication and reverse transcription of nucleic acids. Information on such duplexes may shed some light on the mechanism of these processes. For this purpose, the influence of base composition on the structure of a polypyrimidine--polypurine RNA--DNA duplex r(cucuccuucucuu). d(GAGAGGAAGAGAA) has been studied using 1H, 31P and 13C NMR experiments, molecular modeling (JUMNA program) and NOE back-calculation methods. The resulting structure of the 13-mer hybrid duplex shows that the RNA strand is in the expected A-type conformation while the DNA strand is in a very flexible conformation. In the DNA strand, the desoxyribose sugars retain the C2'-endo B-type conformation. The duplex helical parameters (such as inclination, twist and displacement of the bases) are close to the A-type conformation. No bending was observed for the global axis curvature. The major groove width is close to the B-form value and the minor groove width is intermediate between standard values for A and B-forms. These results are in favour of the independence of minor groove size (where RNase H interacts) and the base composition of the hybrid duplexes.

Bending and adaptability to proteins of the cAMP DNA-responsive element: molecular dynamics contrasted with NMR

DNA bending is assumed to play a crucial role during recognition of the cAMP-responsive element (CRE) by transcription factors. However, diverging results have been obtained for the bending direction of the unbound double helix. The refined NMR structures present a bend directed toward the minor groove, while biochemical methods conclude that there is a bend toward the major groove. The present 10-ns molecular dynamics (MD) simulation of d(GAGATGACGTCATCTC)(2), which contains the octamer CRE in its center, was carried out with AMBER in explicit water and counterions. It shows that CRE is a flexible segment, although it is bent slightly toward the major groove (7 degrees -8 degrees ) on the average. The MD structure agrees with both the biochemical results and unrefined NMR data. The divergence with the NMR refined structures suggests an improper electrostatic parameterization in the refinement software. The malleability of the central CpG is certainly the major contribution to the curving of the whole CRE segment in both the unbound and bound states. Comparison with the crystal structure of CRE bound to GCN4 shows that the deformation induced by the protein is concentrated mainly on the CpG step, rendering the bound structure of CRE closer to the structure of the 12-0 tetradecanoylphorbol-beta-acetate-responsive element.

An A-type double helix of DNA having B-type puckering of the deoxyribose rings

DNA usually adopts structure B in aqueous solution, while structure A is preferred in mixtures of trifluoroethanol (TFE) with water. However, the octamer d(CCCCGGGG) and other d(C(n)G(n)) fragments of DNA provide CD spectra that suggest that the base-pairs are stacked in an A-like fashion even in aqueous solution. Yet, d(CCCCGGGG) undergoes a cooperative TFE-induced transition into structure A, indicating that an important part of the aqueous duplex retains structure B. NMR spectroscopy shows that puckering of the deoxyribose rings is of the B-type. Hence, combination of the information provided by CD spectroscopy and NMR spectroscopy suggests an unprecedented double helix of DNA in which A-like base stacking is combined with B-type puckering of the deoxyribose rings. In order to determine whether this combination is possible, we used molecular dynamics to simulate the duplex of d(CCCCGGGG). Remarkably, the simulations, completely unrestrained by the experimental data, provided a very stable double helix of DNA, exhibiting just the intermediate B/A features described above. The double helix contained well-stacked guanine bases but almost unstacked cytosine bases. This generated a hole in the double helix center, which is a property characteristic for A-DNA, but absent from B-DNA. The minor groove was narrow at the double helix ends but wide at the central CG step where the Watson-Crick base-pairs were buckled in opposite directions. The base-pairs stacked tightly at the ends but stacking was loose in the duplex center. The present double helix, in which A-like base stacking is combined with B-type sugar puckering, is relevant to replication and transcription because both of these phenomena involve a local B-to-A transition.

31P NMR analysis of the DNA conformation induced by protein binding SRY/DNA complexes

Complexes of the HMG box protein SRY with two duplexes of 8 and 14 base pairs have been studied by 31P NMR and complete assignment of all phosphorus signals of the bound DNA duplexes are presented. While for the free DNA, all 31P signals display limited spectral dispersion (< 0.8 p.p.m.) for the bound duplexes, 31P resonances are spread over 2 p.p.m. Based on the previously published 3D structure of hSRY-HMG, with the 8 mer it is demonstrated that the upfield shifted resonances correspond to the site of partial intercalation of an isoleucine side chain into the DNA. Moreover, the observation of significant difference in linewidths between the two duplexes allows to estimate lifetime of the complexes from 31P-31P 2D exchange experiments.

Sequence-dependent DNA structure: tetranucleotide conformational maps

A database of X-ray crystal structures of double helical DNA oligomers has been used to analyse the role of the sugar-phosphate backbone in coupling the conformational properties of neighbouring dinucleotide steps. The base step parameters which are most strongly coupled to the backbone degrees of freedom are slide and shift, and these are the two dinucleotide step parameters which show strong correlations along a sequence: the value of slide follows the values in the neighbouring steps, whereas shift tends to alternate. This conformational coupling is mediated by the shared furanose rings at the step junctions: a change in the value of slide causes a change in the mean value of the same strand 3' and 5'-chi torsion angle, and a change in the mean value of the 3' and 5' sugar pseudo-rotation phase angle, P; a change in the value of shift causes a difference between the same strand 3' and 5'-chi in A-DNA and a difference between the 3' and 5'-P in B-DNA. We have used a database of tetranucleotide X-ray crystal structures to parameterise a simple model for the coupling of slide and shift. Using this junction model together with our dinucleotide step potential energy maps described previously, we can in principle calculate the structure of any DNA oligomer. The parameterisation indicates that the rotational step parameters are accurate to within 5 degrees, and the translational step parameters are accurate to within 0.5 A. The model has been used to study the potential energy surfaces of all possible tetranucleotide sequences, and the calculations agree well with the experimental data from X-ray crystal structures. Some dinucleotide steps are context independent (AA/TT, AT and TA), because the conformational properties of all possible neighbouring steps are compatible. When the conformational properties of the neighbours are not compatible, the behaviour of a step cannot be understood at the dinucleotide level. Thus the conformations of CG, GC and GG/CC are all strongly context dependent. The remaining mixed sequence steps show weakly context-dependent behaviour. The approach allows the calculation of the relative stability and flexibility of tetranucleotide sequences, and the results indicate why TATA is used as an origin of replication. Clear predictions are made about sequences which have not yet been characterised crystallographically. In particular, poly(CCA).poly(TGG) is predicted to have an unusual structure which lies between the C and D-DNA polymorphs.

A DNA hairpin with a single residue loop closed by a strongly distorted Watson-Crick G x C base-pair

Our previous NMR and modeling studies have shown that the single-stranded 19mer oligonucleotides d(AGCTTATC-ATC-GATAA GCT) -ATC- and d(AGCTTATC-GAT-GATAAGCT) -GAT- encompassing the strongest topoisomerase II cleavage site in pBR322 DNA could form stable hairpin structures. A new sheared base-pair, the pyrimidine-purine C x A, was found to close the single base -ATC- loop, while -GAT- displayed a flexible loop of three/five residues with no stabilizing interactions. Now we report a structural study on -GAC-, an analog of -GAT-, derived through the substitution of the loop residue T by C. The results obtained from NMR, non-denaturing PAGE, UV-melting, circular dichroism experiments and restrained molecular dynamics indicate that -GAC- adopts a hairpin structure folded through a single residue loop. In the -GAC- hairpin the direction of the G9 sugar is reversed relative to the C8 sugar, thus pushing the backbone of the loop into the major groove. The G9 x C11 base-pair closing the loop is thus neither a sheared base-pair nor a regular Watson-Crick one. Although G9 and C11 are paired through hydrogen bonds of Watson-Crick type, the base-pair is not planar but rather adopts a wedge-shaped geometry with the two bases stacked on top of each other in the minor groove. The distortion decreases the sugar C1'-C1' distance between the paired G9 and C11, to 8 A versus 11 A in the standard B-DNA. The A10 residue at the center of the loop interacts with the G9 x C11 base-pair, and seems to contribute to the extra thermal stability displayed by -GAC- compared to -GAT-. Test calculations allowed us to identify the experimental NOEs critical for inducing the distorted G.C Watson-Crick base-pair. The preference of -GAC- for a hairpin structure rather than a duplex is confirmed by the diffusion constant values obtained from pulse-field gradient NMR experiments. All together, the results illustrate the high degree of plasticity of single-stranded DNAs which can accommodate a variety of turn-loops to fold up on themselves.

Impact of C5-cytosine methylation on the solution structure of d(GAAAACGTTTTC)2. An NMR and molecular modelling investigation

The solution structures of d(GAAAACGTTTTC)2 and of its methylated derivative d(GAAAAMe5CGTTTTC)2 have been determined by NMR and molecular modelling in order to examine the impact of cytosine methylation on the central CpG conformation. Detailed 1H NMR and 31P NMR investigation of the two oligomers includes quantitative NOESY, 2D homonuclear Hartmann-Hahn spectroscopy, double-quantum-filtered COSY and heteronuclear 1H-31P correlation. Back-calculations of NOESY spectra and simulations of double-quantum-filtered COSY patterns were performed to gain accurate information on interproton distances and sugar phase angles. Molecular models under experimental constraints were generated by energy minimization by means of the molecular mechanics program JUMNA. The MORASS software was used to iteratively refine the structures obtained. After methylation, the oligomer still has a B-DNA conformation. However, there are differences in the structural parameters and the thermal stability as compared to the unmethylated molecule. Careful structural analysis shows that after methylation CpG departs from the usual conformation observed in other ACGT tetramers with different surroundings. Subtle displacements of bases, sugars and backbone imposed by the steric interaction of the two methyl groups inside the major groove are accompanied by severe pinching of the minor groove at the C-G residues.

Helix morphology changes in B-DNA induced by spontaneous B(I)<==>B(II) substrate interconversion

Investigations of spontaneous, i.e. not forced, B-DNA's B(I)<==>B(II) substate transitions are carried out on the d(CGCGAATTCGCG)2 EcoRI dodecamer sequence using Molecular Dynamics Simulations. Analysis of the resulting transition processes with respect to the backbone angles reveals concerted changes not only for backbone angles epsilon, zeta, and beta, but also for the 5'-delta and 5'-chi angles. For alpha and delta inside the interconverting base step, a change is seen in short lived B(II) conformers. With respect to base morphology distinct changes are observed for buckle, propeller twist, shift, roll and twist, as well as x-displacement and tip. The base mainly involved in the changes is identified as the base preceding the interconverting phosphate. Altogether single B(I)<==>B(II) interconversions result only in local distortions represented by the larger spread of most parameters. Comparison of the atomic positional fluctuations derived from the simulation with those obtained from the static X-ray structure results in striking similarities.

Sequence effects on energetic and structural properties of phosphorothioate DNA: a molecular modelling study

Phosphorothioate (PS) oligonucleotides constitute a new class of potent drugs, resulting from the replacement of one anionic oxygen of the phosphodiester backbone by one sulphur atom. This replacement confers chirality to the phosphorus atom (PSS or PSR) and alters the energetic, structural and biological properties of B-DNA. These properties were assessed by molecular mechanics calculations on a set of regular sequences, d(YR)8.d(YR)8 and d(RR)8.d(YY)8 (R, purine; Y, pyrimidine). Results indicated: (i) destabilisation of both the PS(R)and the PSS oligomers, the loss of total energy being mainly due to a variation in the electrostatic term; (ii) an additional chirality effect, due to van der Waals and backbone angle energies, larger for PSS oligomers than for PSR oligomers; (iii) a clear sequence effect on stability, particularly from the base immediately preceding the PS group. Even though the PS group alters the stability of oligomers, it does not significantly modify the conformation. Altogether, our molecular modelling data parallel the available experimental data. Our results reveal that sequence effects on the energetic properties of PS oligomers are local and additive. Therefore, studies of the set of the 10 unique double-stranded modified dinucleotide steps included in regular oligomers could be used to predict the behaviour of any double-stranded PS-DNA.

Conformational variation of the central CG site in d(ATGACGTCAT)2 and d(GAAAACGTTTTC)2. An NMR, molecular modelling and 3D-homology investigation

The determination of the solution structure of two self-complementary oligomers d(ATGACGTCAT)2 (CG10) and d(GAAAACGTTTTC)2 (CG12), both containing the 5'-pur-ACGT-pyr-3' sequence, is reported. The impact of the base context on the conformation of the central CpG site has been examined by a combined approach of: (a) 2D 1H-NMR and 31P-NMR; (b) molecular mechanics under experimental constraints; (c) back-calculations of NOESY spectra and iterative refinements of distances; and (d) 3D-homology search of the central tetrad ACGT within the complete oligonucleotides. A full NMR study of each fragment is achieved by means of standard 2D experiments: NOESY, 2D homonuclear Hartmann-Hahn spectroscopy, double-quantum-filtered COSY and heteronuclear 1H-31P correlation. Sugar phase angle, epsilon-zeta difference angle and NOE-derived distances are input as experimental constraints to generate molecular models by energy minimization with the help of jumna. The morass program is used to iteratively refine the structures obtained. The similarity of the two ACGTs within the whole oligonucleotides is investigated. Both the decamer and the dodecamer adopt a B-like DNA conformation. However, the helical parameters within this conformational type are significantly different in CG12 and CG10. The central CpG step conformation is not locked by its nearest environment (5'A and 3'T) as seen from the structural analysis of ACGT in the two molecules. In CG12, despite the presence of runs of A-T pairs, CpG presents a high twist of 43 degrees and a sugar phase at the guanine of about 180 degrees, previously observed in other ACGT-containing-oligomers. Conversely, ACGT in CG10 exhibits strong inclinations, positive rolls, a flat profile of sugar phase, twist and glycosidic angles, as a result of the nucleotide sequence extending beyond the tetrad. The structural specificity of CG10 and its flexibility (as reflected by its energy) are tentatively related to the process of recognition of the cyclic AMP response element by its cognate protein.

Molecular dynamic simulations of environment and sequence dependent DNA conformations: the development of the BMS nucleic acid force field and comparison with experimental results

Molecular dynamic (MD) simulations using the BMS nucleic acid force field produce environment and sequence dependent DNA conformations that closely mimic experimentally derived structures. The parameters were initially developed to reproduce the potential energy surface, as defined by quantum mechanics, for a set of small molecules that can be used as the building blocks for nucleic acid macromolecules (dimethyl phosphate, cyclopentane, tetrahydrofuran, etc.). Then the dihedral parameters were fine tuned using a series of condensed phase MD simulations of DNA and RNA (in zero added salt, 4M NaCl, and 75% ethanol solutions). In the tuning process the free energy surface for each dihedral was derived from the MD ensemble and fitted to the conformational distributions and populations observed in 87 A- and B-DNA x-ray and 17 B-DNA NMR structures. Over 41 nanoseconds of MD simulations are presented which demonstrate that the force field is capable of producing stable trajectories, in the correct environments, of A-DNA, double stranded A-form RNA, B-DNA, Z-DNA, and a netropsin-DNA complex that closely reproduce the experimentally determined and/or canonical DNA conformations. Frequently the MD averaged structure is closer to the experimentally determined structure than to the canonical DNA conformation. MD simulations of A- to B- and B- to A-DNA transitions are also shown. A-DNA simulations in a low salt environment cleanly convert into the B-DNA conformation and converge into the RMS space sampled by a low salt simulation of the same sequence starting from B-DNA. In MD simulations using the BMS force field the B-form of d(GGGCCC)2 in a 75% ethanol solution converts into the A-form. Using the same methodology, parameters, and conditions the A-form of d(AAATTT)2 correctly converts into the B-DNA conformation. These studies demonstrate that the force field is capable of reproducing both environment and sequence dependent DNA structures. The 41 nanoseconds (nsec) of MD simulations presented in this paper paint a global picture which suggests that the DNA structures observed in low salt solutions are largely due to the favorable internal energy brought about by the nearly uniform screening of the DNA electrostatics. While the conformations sampled in high salt or mixed solvent environments occur from selective and asymmetric screening of the phosphate groups and DNA grooves, respectively, brought about by sequence induced ion and solvent packing.

Solution structure of a non-palindromic 16 base-pair DNA related to the HIV-1 kappa B site: evidence for BI-BII equilibrium inducing a global dynamic curvature of the duplex

1H and 31P NMR spectroscopy have been used together with molecular modelling to determine the fine structure of a non-palindromic 16 bp DNA containing the NF-kappa B binding site. Much emphasis has been placed upon NMR optimization of both two-dimensional 31P NMR techniques to extract structural information defining the phosphodiester backbone conformation and selective homonuclear 2D COSY experiments to determine sugar conformations. NMR data show evidence for a dynamic behaviour of steps flanking the ten base-pairs of the NF-kappa B binding site. A BI-BII equilibrium at these steps is demonstrated and two models for each extreme conformation are proposed in agreement with NMR data. In the refined BII structures, the NF-kappa B binding site exhibits an intrinsic curvature towards the major groove that is magnified by the four flanking steps in the BII conformation. Furthermore, the base-pairs are translated into the major groove. Thus, we present a novel mode of dynamic intrinsic curvature compatible with the DNA curvature observed in the X-ray structure of the p50-DNA complex.

DNA binding discrimination of the murine DNA cytosine-C5 methyltransferase

Mammalian DNA cytosine-C5 methyltransferase modifies the CpG dinucleotide in the context of many different genomic sequences. A rigorous DNA binding assay was developed for the murine enzyme and used to define how sequences flanking the CpG dinucleotide affect the stability of the enzyme:DNA complex. Oligonucleotides containing a single CpG site form reversible 1:1 complexes with the enzyme that are sequence-specific. A guanine/cytosine-rich 30 base-pair sequence, a mimic of the GC-box cis-element, bound threefold more tightly than an adenine/thymine-rich sequence, a mimic of the cyclic AMP responsive element. However, the binding discrimination between hemi- and unmethylated forms of these DNA substrates was small, as we previously observed at the K(m)DNA level (Biochemistry, 35, 7308-7315 (1996)). Single-stranded substrates are bound much more weakly than double-stranded DNA forms. An in vitro screening method was used to select for CpG flanking sequence preferences of the DNA methyltransferase from a large, divergent population of DNA substrates. After five iterative rounds of increasing selective pressure, guanosine/cytosine-rich sequences were abundant and contributed to binding stabilization for at least 12 base-pairs on either side of a central CpG. Our results suggest a read-out of sequence-dependent conformational features, such as helical flexibility, minor groove dimensions and critical phosphate orientation and mobility, rather than interactions with specific bases over the course of two complete helical turns. Thus, both studies reveal a preference for guanosine/cytosine deoxynucleotides flanking the cognate CpG. The enzyme specificity for similar sequences in the genome may contribute to the in vivo functions of this vital enzyme.

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

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