Does Cytosine Methylation Stabilize the BI Substate of DNA?

Abstract

The local conformation of DNA plays a critical role in the recognition-and-binding process of transcription factors. In particular, the most common form of the DNA double helix, B-DNA, exists as a conformational equilibrium between BI and BII substates that are associated with changes in major and minor groove dimensions. It has been hypothesized that cytosine methylation, the most common epigenetic modification and one that is often found in so-called CpG islands, modulates transcription factor binding affinity by altering the BI-BII equilibrium in the neighborhood of methylation sites. Indeed, previous molecular dynamics studies have suggested that cytosine methylation leads to a uniform stabilization of the BI substate in CpG islands. In this work we first benchmark the ability of the latest Amber DNA force fields and a variety of three-point water models to reproduce BI-BII equilibrium in the Dickerson dodecamer and reduce terminal base pair fraying. Using the best-performing combinations, we then re-examine the conformational effects of cytosine methylation in a prototypical CpG island, (GC)5. We find that the BI substate is stabilized for GpC steps in such sequences, but that the BII substate is stabilized for the CpG steps. More interestingly, the Amber ff99bsc0 DNA force field used in previous studies yields opposite conclusions. These findings highlight the need for continued fixed-charge force field development for nucleic acids, as well as solution state experimental datasets involving epigenetic modification that can serve as benchmarks for simulation accuracy.