Formation of left-handed helices by C2'-fluorinated nucleic acids under physiological salt conditions.

Recent findings in cell biology have rekindled interest in Z-DNA, the left-handed helical form of DNA. We report here that two minimally modified nucleosides, 2'F-araC and 2'F-riboG, induce the formation of the Z-form under low ionic strength. We show that oligomers entirely made of these two nucleosides exclusively produce left-handed duplexes that bind to the Zα domain of ADAR1. The effect of the two nucleotides is so dramatic that Z-form duplexes are the only species observed in 10 mM sodium phosphate buffer and neutral pH, and no B-form is observed at any temperature. Hence, in contrast to other studies reporting formation of Z/B-form equilibria by a preference for purine glycosidic angles in syn, our NMR and computational work revealed that sequential 2'F…H2N and intramolecular 3'H…N3' interactions stabilize the left-handed helix. The equilibrium between B- and Z- forms is slow in the 19F NMR time scale (≥ms), and each conformation exhibited unprecedented chemical shift differences in the 19F signals. This observation led to a reliable estimation of the relative population of B and Z species and enabled us to monitor B-Z transitions under different conditions. The unique features of 2'F-modified DNA should thus be a valuable addition to existing techniques for specific detection of new Z-binding proteins and ligands.

CGeNArate: a sequence-dependent coarse-grained model of DNA for accurate atomistic MD simulations of kb-long duplexes.

We present CGeNArate, a new model for molecular dynamics simulations of very long segments of B-DNA in the context of biotechnological or chromatin studies. The developed method uses a coarse-grained Hamiltonian with trajectories that are back-mapped to the atomistic resolution level with extreme accuracy by means of Machine Learning Approaches. The method is sequence-dependent and reproduces very well not only local, but also global physical properties of DNA. The efficiency of the method allows us to recover with a reduced computational effort high-quality atomic-resolution ensembles of segments containing many kilobases of DNA, entering into the gene range or even the entire DNA of certain cellular organelles.

Zeta potential of Z-DNA: A new signature to study B-Z transition in linear and branched DNA.

Zeta potential is commonly referred as surface charge density and is a key factor in modulating the structural and functional properties of nucleic acids. Although the negative charge density of B-DNA is well understood, there is no prior description of the zeta potential measurement of Z-DNA. In this study, for the first time we discover the zeta potential difference between B-DNA and lanthanum chloride-induced Z-DNA. A series of linear repeat i.e. (CG)n and (GC)n DNA as well as branched DNA (bDNA) structures was used for the B-to-Z DNA transition. Herein, the positive zeta potential of Z-DNA has been demonstrated as a powerful tool to discriminate between B-form and Z-form of DNA. The generality of the approach has been validated both in linear and bDNA nanostructures. Thus, we suggest zeta potential can be used as an ideal signature for the left-handed Z-DNA.

Z-DNA as a Tool for Nuclease-Free DNA Methyltransferase Assay.

Methylcytosines in mammalian genomes are the main epigenetic molecular codes that switch off the repertoire of genes in cell-type and cell-stage dependent manners. DNA methyltransferases (DMT) are dedicated to managing the status of cytosine methylation. DNA methylation is not only critical in normal development, but it is also implicated in cancers, degeneration, and senescence. Thus, the chemicals to control DMT have been suggested as anticancer drugs by reprogramming the gene expression profile in malignant cells. Here, we report a new optical technique to characterize the activity of DMT and the effect of inhibitors, utilizing the methylation-sensitive B-Z transition of DNA without bisulfite conversion, methylation-sensing proteins, and polymerase chain reaction amplification. With the high sensitivity of single-molecule FRET, this method detects the event of DNA methylation in a single DNA molecule and circumvents the need for amplification steps, permitting direct interpretation. This method also responds to hemi-methylated DNA. Dispensing with methylation-sensitive nucleases, this method preserves the molecular integrity and methylation state of target molecules. Sparing methylation-sensing nucleases and antibodies helps to avoid errors introduced by the antibody's incomplete specificity or variable activity of nucleases. With this new method, we demonstrated the inhibitory effect of several natural bio-active compounds on DMT. All taken together, our method offers quantitative assays for DMT and DMT-related anticancer drugs.

Z-DNA as a Touchstone for Additive Empirical Force Fields and a Refinement of the Alpha/Gamma DNA Torsions for AMBER.

Although current AMBER force fields are relatively accurate for canonical B-DNA, many noncanonical structures are still described incorrectly. As noncanonical motifs are attracting increasing attention due to the role they play in living organisms, further improvement is desirable. Here, we have chosen the Z-DNA molecule, which can be considered a touchstone of the universality of empirical force fields, since the noncanonical α and γ backbone conformations native to Z-DNA are also found in protein-DNA complexes, i-motif DNA, and other noncanonical DNAs. We show that spurious α/γ conformations occurring in simulations with current AMBER force fields, OL15 and bsc1, are largely due to inaccurate α/γ parametrization. Moreover, stabilization of native Z-DNA substates involving γ = trans conformations appears to be in conflict with the correct description of the canonical B-DNA structure. Because the balance of the native and spurious conformations is influenced by nonadditive effects, this is a difficult case for an additive dihedral energy scheme such as AMBER. We propose new α/γ parameters, denoted OL21, and show that they improve the stability of native α/γ Z-DNA substates while keeping the canonical DNA description virtually unchanged, thus representing a reasonable compromise within the additive force field framework. Although further extensive testing is needed, the new modification appears to be a promising step toward a more reliable description of noncanonical DNA motifs and provides the best performance for Z-DNA molecules among current AMBER force fields.

Copper chloride complexes with substituted 4'-phenyl-terpyridine ligands: synthesis, characterization, antiproliferative activities and DNA interactions.

Eleven copper chloride coordination compounds (1-11) with 4'-(4'-substituted-phenyl)-2,2':6',2''-terpyridine ligands bearing hydrogen (L1), cyano (L2), p-hydroxyl (L3), m-hydroxyl (L4), o-hydroxyl (L5), methoxyl (L6), iodo (L7), bromo (L8), chloro (L9), fluoro (L10) or methylsulfonyl (L11) were prepared and characterized by IR spectroscopy, elemental analysis and single crystal X-ray diffraction. Antiproliferative activities against tumor cells were investigated and DNA interactions were studied by circular dichroism spectroscopy and molecular modeling methods. In vitro data demonstrate that all the compounds exhibit higher antiproliferative activities as compared to cisplatin against five human carcinoma cell lines: A549, Bel-7402, Eca-109, HeLa and MCF-7. Compound 6 with methoxyl shows the best anti-proliferation activity. Spectrophotometric results reveal the strong affinity of the compounds for binding with DNA as intercalators and induce DNA conformational transitions. The results of molecular docking studies show that the compounds interact with DNA through π-π stacking, van der Waals forces, hydrophobic interactions and hydrogen bonds. The binding energies between compound 11 and three macromolecules, including DNA duplex, oligonucleotide and DNA-Topo I complex, are the lowest. The binding stability of compounds containing hydroxyl, methoxy and methylsulfonyl groups with biological macromolecules mainly relies on the hydrogen bonds. The ability of a compound to form hydrogen bonds can promote its binding to biological targets, thereby exhibiting high antiproliferative activity.

A novel ß-hairpin peptide derived from the ARC repressor selectively interacts with the major groove of B-DNA.

Transcription factors (TFs) have a remarkable role in the homeostasis of the organisms and there is a growing interest in how they recognize and interact with specific DNA sequences. TFs recognize DNA using a variety of structural motifs. Among those, the ribbon-helix-helix (RHH) proteins, exemplified by the MetJ and ARC repressors, form dimers that insert antiparallel ß-sheets into the major groove of DNA. A great chemical challenge consists of using the principles of DNA recognition by TFs to design minimized peptides that maintain the DNA affinity and specificity characteristics of the natural counterparts. In this context, a peptide mimic of an antiparallel ß-sheet is very attractive since it can be obtained by a single peptide chain folding in a ß-hairpin structure and can be as short as 14 amino acids or less. Herein, we designed eight linear and two cyclic dodeca-peptides endowed with ß-hairpins. Their DNA binding properties have been investigated using fluorescence spectroscopy together with the conformational analysis through circular dichroism and solution NMR. We found that one of our peptides, peptide 6, is able to bind DNA, albeit without sequence selectivity. Notably, it shows a topological selectivity for the major groove of the DNA which is the interaction site of ARC and many other DNA-binding proteins. Moreover, we found that a type I' ß-hairpin folding pattern is a favorite peptide structure for interaction with the B-DNA major groove. Peptide 6 is a valuable lead compound for the development of novel analogs with sequence selectivity.

Systematic Approach to Find the Global Minimum of Relaxation Dispersion Data for Protein-Induced B-Z Transition of DNA.

Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion spectroscopy is commonly used for quantifying conformational changes of protein in µs-to-ms timescale transitions. To elucidate the dynamics and mechanism of protein binding, parameters implementing CPMG relaxation dispersion results must be appropriately determined. Building an analytical model for multi-state transitions is particularly complex. In this study, we developed a new global search algorithm that incorporates a random search approach combined with a field-dependent global parameterization method. The robust inter-dependence of the parameters carrying out the global search for individual residues (GSIR) or the global search for total residues (GSTR) provides information on the global minimum of the conformational transition process of the Zα domain of human ADAR1 (hZαADAR1)-DNA complex. The global search results indicated that a α-helical segment of hZαADAR1 provided the main contribution to the three-state conformational changes of a hZαADAR1-DNA complex with a slow B-Z exchange process. The two global exchange rate constants, kex and kZB, were found to be 844 and 9.8 s-1, respectively, in agreement with two regimes of residue-dependent chemical shift differences-the "dominant oscillatory regime" and "semi-oscillatory regime". We anticipate that our global search approach will lead to the development of quantification methods for conformational changes not only in Z-DNA binding protein (ZBP) binding interactions but also in various protein binding processes.

Opposite roles of transcription elongation factors Spt4/5 and Elf1 in RNA polymerase II transcription through B-form versus non-B DNA structures.

Transcription elongation can be affected by numerous types of obstacles, such as nucleosome, pausing sequences, DNA lesions and non-B-form DNA structures. Spt4/5 and Elf1 are conserved transcription elongation factors that promote RNA polymerase II (Pol II) bypass of nucleosome and pausing sequences. Importantly, genetic studies have shown that Spt4/5 plays essential roles in the transcription of expanded nucleotide repeat genes associated with inherited neurological diseases. Here, we investigate the function of Spt4/5 and Elf1 in the transcription elongation of CTG•CAG repeat using an in vitro reconstituted yeast transcription system. We found that Spt4/5 helps Pol II transcribe through the CTG•CAG tract duplex DNA, which is in good agreement with its canonical roles in stimulating transcription elongation. In sharp contrast, surprisingly, we revealed that Spt4/5 greatly inhibits Pol II transcriptional bypass of CTG and CAG slip-out structures. Furthermore, we demonstrated that transcription elongation factor Elf1 individually and cooperatively with Spt4/5 inhibits Pol II bypass of the slip-out structures. This study uncovers the important functional interplays between template DNA structures and the function of transcription elongation factors. This study also expands our understanding of the functions of Spt4/5 and Elf1 in transcriptional processing of trinucleotide repeat DNA.

Sequence-dependent cost for Z-form shapes the torsion-driven B-Z transition via close interplay of Z-DNA and DNA bubble.

Despite recent genome-wide investigations of functional DNA elements, the mechanistic details about their actions remain elusive. One intriguing possibility is that DNA sequences with special patterns play biological roles, adopting non-B-DNA conformations. Here we investigated dynamics of thymine-guanine (TG) repeats, microsatellite sequences and recurrently found in promoters, as well as cytosine-guanine (CG) repeats, best-known Z-DNA forming sequence, in the aspect of Z-DNA formation. We measured the energy barriers of the B-Z transition with those repeats and discovered the sequence-dependent penalty for Z-DNA generates distinctive thermodynamic and kinetic features in the torque-induced transition. Due to the higher torsional stress required for Z-form in TG repeats, a bubble could be induced more easily, suppressing Z-DNA induction, but facilitate the B-Z interconversion kinetically at the transition midpoint. Thus, the Z-form by TG repeats has advantages as a torsion buffer and bubble selector while the Z-form by CG repeats likely behaves as torsion absorber. Our statistical physics model supports quantitatively the populations of Z-DNA and reveals the pivotal roles of bubbles in state dynamics. All taken together, a quantitative picture for the transition was deduced within the close interplay among bubbles, plectonemes and Z-DNA.

Biophysical interaction between self-assembled branched DNA nanostructures with bovine serum albumin and bovine liver catalase.

Branched DNA (bDNA) nanostructures have emerged as self-assembled biomaterials and are being considered for biomedical applications. Herein, we report the biophysical interaction between self-assembled bDNA nanostructure with circulating protein bovine serum albumin (BSA) and cellular enzyme bovine liver catalase (BLC). The binding between bDNA and BSA or BLC was confirmed through the decrease in fluorescence spectra. The Stern-Volmer data supports for non-covalent bonding with ~1 binding site in case of BSA and BLC thus advocating a static binding. Furthermore, FTIR and ITC study confirmed the binding of bDNAs with proteins through hydrogen bonding and van der Waals interaction. The negative free energy observed in ITC represent spontaneous reaction for BLC-bDNA interaction. The biophysical interaction between bDNA nanostructures and proteins was also supported by DLS and zeta potential measurement. With an increase in bDNA concentrations up to 100 nM, no significant change in absorbance and CD spectra was observed for both BLC and BSA which suggests structural stability and unaffected secondary conformation of proteins in presence of bDNA. Furthermore, the catalytic activity of BLC was unaltered in presence of bDNAscr even with increasing the incubation period from 1 h to 24 h. Interestingly, the time-dependent decrease in activity of BLC was protected by bDNAmix. The thermal melting study suggests a higher Tm value for proteins in presence of bDNAmix which demonstrates that interaction with bDNAmix increases the thermal stability of proteins. Collectively these data suggest that self-assembled DNA nanostructure may bind to BSA for facilitating circulation in plasma or binding to intracellular proteins like BLC for stabilization, however the secondary conformation of protein or catalytic activity of enzyme is unaltered in presence of bDNA nanostructure. Thus, the newly established genomic sequence-driven self-assembled DNA nanostructure can be explored for in vitro or in vivo experimental work in recent future.

Spin-projected QM/MM Free Energy Simulations for Oxidation Reaction of Guanine in B-DNA by Singlet Oxygen.

Guanine is the most susceptible base to oxidation damage induced by reactive oxygen species including singlet oxygen (1 O2 , 1 Δg ). We clarify whether the first step of guanine oxidation in B-DNA proceeds via either a zwitterionic or a diradical intermediate. The free energy profiles are calculated by means of a combined quantum mechanical and molecular mechanical (QM/MM) method coupled with the adaptive biasing force (ABF) method. To describe the open-shell electronic structure of 1 O2 correctly, the broken-symmetry spin-unrestricted density functional theory (BS-UDFT) with an approximate spin projection (AP) correction is applied to the QM region. We find that the effect of spin contamination on the activation and reaction free energies is up to ∼8 kcal mol-1 , which is too large to be neglected. The QM(AP-ULC-BLYP)/MM-based free energy calculations also reveal that the reaction proceeds through a diradical transition state, followed by a conversion to a zwitterionic intermediate. Our computed activation energy of 5.2 kcal mol-1 matches experimentally observed range (0∼6 kcal mol-1 ).

Protein-induced B-Z transition of DNA duplex containing a 2'-OMe guanosine.

Structural transformation of the canonical right-handed helix, B-DNA, to the non-canonical left-handed helix, Z-DNA, can be induced by the Zα domain of the human RNA editing enzyme ADAR1 (hZαADAR1). To characterize the site-specific preferences of binding and structural changes in DNA containing the 2'-O-methyl guanosine derivative (mG), titration of the imino proton spectra and chemical shift perturbations were performed on hZαADAR1 upon binding to Z-DNA. The structural transition between B-Z conformation as the changing ratio between DNA and protein showed a binding affinity of the modified DNA onto the Z-DNA binding protein similar to wild-type DNA or RNA. The chemical shift perturbation results showed that the overall structure and environment of the modified DNA revealed DNA-like properties rather than RNA-like characteristics. Moreover, we found evidence for two distinct regimes, "Z-DNA Sensing" and "Modification Sensing", based on the site-specific chemical shift perturbation between the DNA (or RNA) binding complex and the modified DNA-hZαADAR1 complex. Thus, we propose that modification of the sugar backbone of DNA with 2'-O-methyl guanosine promotes the changes in the surrounding α3 helical structural segment as well as the non-perturbed feature of the ß-hairpin region.

Environmental Effects on Guanine-Thymine Mispair Tautomerization Explored with Quantum Mechanical/Molecular Mechanical Free Energy Simulations.

DNA bases can adopt energetically unfavorable tautomeric forms that enable the formation of Watson-Crick-like (WC-like) mispairs, which have been proposed to give rise to spontaneous mutations in DNA and misincorporation errors in DNA replication and translation. Previous NMR and computational studies have indicated that the population of WC-like guanine-thymine (G-T) mispairs depends on the environment, such as the local nucleic acid sequence and solvation. To investigate these environmental effects, herein G-T mispair tautomerization processes are studied computationally in aqueous solution, in A-form and B-form DNA duplexes, and within the active site of a DNA polymerase λ variant. The wobble G-T (wG-T), WC-like G-T*, and WC-like G*-T forms are considered, where * indicates the enol tautomer of the base. The minimum free energy paths for the tautomerization from the wG-T to the WC-like G-T* and from the WC-like G-T* to the WC-like G*-T are computed with mixed quantum mechanical/molecular mechanical (QM/MM) free energy simulations. The reaction free energies and free energy barriers are found to be significantly influenced by the environment. The wG-T→G-T* tautomerization is predicted to be endoergic in aqueous solution and the DNA duplexes but slightly exoergic in the polymerase, with Arg517 and Asn513 providing electrostatic stabilization of G-T*. The G-T*→G*-T tautomerization is also predicted to be slightly more thermodynamically favorable in the polymerase relative to these DNA duplexes. These simulations are consistent with an experimentally driven kinetic misincorporation model suggesting that G-T mispair tautomerization occurs in the ajar polymerase conformation or concertedly with the transition from the ajar to the closed polymerase conformation. Furthermore, the order of the associated two proton transfer reactions is predicted to be different in the polymerase than in aqueous solution and the DNA duplexes. These studies highlight the impact of the environment on the thermodynamics, kinetics, and fundamental mechanisms of G-T mispair tautomerization, which plays a role in a wide range of biochemically important processes.

The complex formed between a synthetic RNA aptamer and the transcription repressor TetR is a structural and functional twin of the operator DNA-TetR regulator complex.

RNAs play major roles in the regulation of gene expression. Hence, designer RNA molecules are increasingly explored as regulatory switches in synthetic biology. Among these, the TetR-binding RNA aptamer was selected by its ability to compete with operator DNA for binding to the bacterial repressor TetR. A fortuitous finding was that induction of TetR by tetracycline abolishes both RNA aptamer and operator DNA binding in TetR. This enabled numerous applications exploiting both the specificity of the RNA aptamer and the efficient gene repressor properties of TetR. Here, we present the crystal structure of the TetR-RNA aptamer complex at 2.7 Å resolution together with a comprehensive characterization of the TetR-RNA aptamer versus TetR-operator DNA interaction using site-directed mutagenesis, size exclusion chromatography, electrophoretic mobility shift assays and isothermal titration calorimetry. The fold of the RNA aptamer bears no resemblance to regular B-DNA, and neither does the thermodynamic characterization of the complex formation reaction. Nevertheless, the functional aptamer-binding epitope of TetR is fully contained within its DNA-binding epitope. In the RNA aptamer complex, TetR adopts the well-characterized DNA-binding-competent conformation of TetR, thus revealing how the synthetic TetR-binding aptamer strikes the chords of the bimodal allosteric behaviour of TetR to function as a synthetic regulator.

A multi-modal coarse grained model of DNA flexibility mappable to the atomistic level.

We present a new coarse grained method for the simulation of duplex DNA. The algorithm uses a generalized multi-harmonic model that can represent any multi-normal distribution of helical parameters, thus avoiding caveats of current mesoscopic models for DNA simulation and representing a breakthrough in the field. The method has been parameterized from accurate parmbsc1 atomistic molecular dynamics simulations of all unique tetranucleotide sequences of DNA embedded in long duplexes and takes advantage of the correlation between helical states and backbone configurations to derive atomistic representations of DNA. The algorithm, which is implemented in a simple web interface and in a standalone package reproduces with high computational efficiency the structural landscape of long segments of DNA untreatable by atomistic molecular dynamics simulations.

Interaction of camel Lactoferrin derived peptides with DNA: a molecular dynamics study.

BACKGROUND: Lactoferrampin (LFampin), Lactoferricin (LFcin), and LFchimera are three well-known antimicrobial peptides derived from Lactoferrin and proposed as alternatives for antibiotics. Although the intracellular activity of these peptides has been previously demonstrated, their mode of action is not yet fully understood. Here, we performed a molecular dynamics simulation study to understand the molecular interactions between camel Lactoferrin derived peptides, including CLFampin, CLFcin, and CLFchimera, and DNA as an important intracellular target. RESULTS: Our results indicate that all three peptides bind to DNA, albeit with different propensities, with CLFchimera showing the highest binding affinity. The secondary structures of the peptides, modeled on Lactoferrin, did not undergo significant changes during simulation, supporting their functional relevance. Main residues involved in the peptide-DNA interaction were identified based on binding free energy estimates calculated over 200 ns, which, as expected, confirmed strong electrostatic interactions between DNA phosphate groups and positively charged peptide side chains. Interaction between the different concentrations of CLFchimera and DNA revealed that after binding of four copies of CLFchimera to DNA, hydrogen bonds between the two strands of DNA start to break from one of the termini. CONCLUSIONS: Importantly, our results revealed that there is no DNA-sequence preference for peptide binding, in line with a broad antimicrobial activity. Moreover, the results showed that the strength of the interaction between DNA and CLFchimera is concentration dependent. The insight provided by these results can be used for the rational redesign of natural antimicrobial peptides targeting the bacterial DNA.

Modeling DNA Flexibility: Comparison of Force Fields from Atomistic to Multiscale Levels.

Accurate parametrization of force fields (FFs) is of ultimate importance for computer simulations to be reliable and to possess a predictive power. In this work, we analyzed, in multi-microsecond simulations of a 40-base-pair DNA fragment, the performance of four force fields, namely, the two recent major updates of CHARMM and two from the AMBER family. We focused on a description of double-helix DNA flexibility and dynamics both at atomistic and at mesoscale level in coarse-grained (CG) simulations. In addition to the traditional analysis of different base-pair and base-step parameters, we extended our analysis to investigate the ability of the force field to parametrize a CG DNA model by structure-based bottom-up coarse-graining, computing DNA persistence length as a function of ionic strength. Our simulations unambiguously showed that the CHARMM36 force field is unable to preserve DNA's structural stability at over-microsecond time scale. Both versions of the AMBER FF, parmbsc0 and parmbsc1, showed good agreement with experiment, with some bias of parmbsc0 parameters for intermediate A/B form DNA structures. The CHARMM27 force field provides stable atomistic trajectories and overall (among the considered force fields) the best fit to experimentally determined DNA flexibility parameters both at atomistic and at mesoscale level.

Structures of GapR reveal a central channel which could accommodate B-DNA.

GapR is a nucleoid-associated protein required for the cell cycle of Caulobacter cresentus. We have determined new crystal structures of GapR to high resolution. As in a recently published structure, a GapR monomer folds into one long N-terminal α helix and two shorter α helices, and assembles into a tetrameric ring with a closed, positively charged, central channel. In contrast to the conclusions drawn from the published structures, we observe that the central channel of the tetramer presented here could freely accommodate B-DNA. Mutation of six conserved lysine residues lining the cavity and electrophoretic mobility gel shift experiments confirmed their role in DNA binding and the channel as the site of DNA binding. Although present in our crystals, DNA could not be observed in the electron density maps, suggesting that DNA binding is non-specific, which could be important for tetramer-ring translocation along the chromosome. In conjunction with previous GapR structures we propose a model for DNA binding and translocation that explains key published observations on GapR and its biological functions.

The static and dynamic structural heterogeneities of B-DNA: extending Calladine-Dickerson rules.

We present a multi-laboratory effort to describe the structural and dynamical properties of duplex B-DNA under physiological conditions. By processing a large amount of atomistic molecular dynamics simulations, we determine the sequence-dependent structural properties of DNA as expressed in the equilibrium distribution of its stochastic dynamics. Our analysis includes a study of first and second moments of the equilibrium distribution, which can be accurately captured by a harmonic model, but with nonlocal sequence-dependence. We characterize the sequence-dependent choreography of backbone and base movements modulating the non-Gaussian or anharmonic effects manifested in the higher moments of the dynamics of the duplex when sampling the equilibrium distribution. Contrary to prior assumptions, such anharmonic deformations are not rare in DNA and can play a significant role in determining DNA conformation within complexes. Polymorphisms in helical geometries are particularly prevalent for certain tetranucleotide sequence contexts and are always coupled to a complex network of coordinated changes in the backbone. The analysis of our simulations, which contain instances of all tetranucleotide sequences, allow us to extend Calladine-Dickerson rules used for decades to interpret the average geometry of DNA, leading to a set of rules with quantitative predictive power that encompass nonlocal sequence-dependence and anharmonic fluctuations.