The dynamic impact of CpG methylation in DNA

DNA Methylation in the Adaptive Response to Exercise

Emerging evidence published over the past decade has highlighted the role of DNA methylation in skeletal muscle function and health, including as an epigenetic transducer of the adaptive response to exercise. In this review, we aim to synthesize the latest findings in this field to highlight: (1) the shifting understanding of the genomic localization of altered DNA methylation in response to acute and chronic aerobic and resistance exercise in skeletal muscle (e.g., promoter, gene bodies, enhancers, intergenic regions, un-annotated regions, and genome-wide methylation); (2) how these global/regional methylation changes relate to transcriptional activity following exercise; and (3) the factors (e.g., individual demographic or genetic features, dietary, training history, exercise parameters, local epigenetic characteristics, circulating hormones) demonstrated to alter both the pattern of DNA methylation after exercise, and the relationship between DNA methylation and gene expression. Finally, we discuss the changes in non-CpG methylation and 5-hydroxymethylation after exercise, as well as the importance of emerging single-cell analyses to future studies-areas of increasing focus in the field of epigenetics. We anticipate that this review will help generate a framework for clinicians and researchers to begin developing and testing exercise interventions designed to generate targeted changes in DNA methylation as part of a personalized exercise regimen.

The Dynamic Regulation of G-Quadruplex DNA Structures by Cytosine Methylation

It is well known that certain non B-DNA structures, including G-quadruplexes, are key elements that can regulate gene expression. Here, we explore the theory that DNA modifications, such as methylation of cytosine, could act as a dynamic switch by promoting or alleviating the structural formation of G-quadruplex structures in DNA or RNA. The interaction between epigenetic DNA modifications, G4 formation, and the 3D architecture of the genome is a complex and developing area of research. Although there is growing evidence for such interactions, a great deal still remains to be discovered. In vivo, the potential effect that cytosine methylation may have on the formation of DNA structures has remained largely unresearched, despite this being a potential mechanism through which epigenetic factors could regulate gene activity. Such interactions could represent novel mechanisms for important biological functions, including altering nucleosome positioning or regulation of gene expression. Furthermore, promotion of strand-specific G-quadruplex formation in differentially methylated genes could have a dynamic role in directing X-inactivation or the control of imprinting, and would be a worthwhile focus for future research.

Detecting DNA Methylations in the Hyperthermoacidophilic Crenarchaeon Sulfolobus acidocaldarius Using SMRT Sequencing

DNA methylations are one of the most well-known epigenetic modifications along with histone modifications and noncoding RNAs. They are found at specific sites along the DNA in all domains of life, with 5-mC and 6-mA/4-mC being well-characterized in eukaryotes and bacteria respectively, and they have not only been described as contributing to the structure of the double helix itself but also as regulators of DNA-based processes such as replication, transcription, and recombination. Different methods have been developed to accurately identify and/or map methylated motifs to decipher the involvement of DNA methylations in regulatory networks that affect the cellular state.Although DNA methylations have been detected along archaeal genomes, their involvement as regulators of DNA-based processes remains the least known. To highlight the importance of DNA methylations in the control of key cellular mechanisms and their dynamics in archaea cells, we have used single-molecule real-time (SMRT) sequencing. This sequencing technology allows the identification and direct mapping of the methylated motifs along the genome of an organism. In this chapter, we present a step-by-step protocol for detecting DNA methylations in the hyperthermophilic crenarchaeon Sulfolobus acidocaldarius using SMRT sequencing. This protocol can easily be adapted to other prokaryotes.

Single molecule characterization of the binding kinetics of a transcription factor and its modulation by DNA sequence and methylation

The interaction of transcription factors with their response elements in DNA is emerging as a highly complex process, whose characterization requires measuring the full distribution of binding and dissociation times in a well-controlled assay. Here, we present a single-molecule assay that exploits the thermal fluctuations of a DNA hairpin to detect the association and dissociation of individual, unlabeled transcription factors. We demonstrate this new approach by following the binding of Egr1 to its consensus motif and the three binding sites found in the promoter of the Lhb gene, and find that both association and dissociation are modulated by the 9 bp core motif and the sequences around it. In addition, CpG methylation modulates the dissociation kinetics in a sequence and position-dependent manner, which can both stabilize or destabilize the complex. Together, our findings show how variations in sequence and methylation patterns synergistically extend the spectrum of a protein's binding properties, and demonstrate how the proposed approach can provide new insights on the function of transcription factors.

Small extracellular vesicles in cancer

Extracellular vesicles (EV) are lipid-bilayer enclosed vesicles in submicron size that are released from cells. A variety of molecules, including proteins, DNA fragments, RNAs, lipids, and metabolites can be selectively encapsulated into EVs and delivered to nearby and distant recipient cells. In tumors, through such intercellular communication, EVs can regulate initiation, growth, metastasis and invasion of tumors. Recent studies have found that EVs exhibit specific expression patterns which mimic the parental cell, providing a fingerprint for early cancer diagnosis and prognosis as well as monitoring responses to treatment. Accordingly, various EV isolation and detection technologies have been developed for research and diagnostic purposes. Moreover, natural and engineered EVs have also been used as drug delivery nanocarriers, cancer vaccines, cell surface modulators, therapeutic agents and therapeutic targets. Overall, EVs are under intense investigation as they hold promise for pathophysiological and translational discoveries. This comprehensive review examines the latest EV research trends over the last five years, encompassing their roles in cancer pathophysiology, diagnostics and therapeutics. This review aims to examine the full spectrum of tumor-EV studies and provide a comprehensive foundation to enhance the field. The topics which are discussed and scrutinized in this review encompass isolation techniques and how these issues need to be overcome for EV-based diagnostics, EVs and their roles in cancer biology, biomarkers for diagnosis and monitoring, EVs as vaccines, therapeutic targets, and EVs as drug delivery systems. We will also examine the challenges involved in EV research and promote a framework for catalyzing scientific discovery and innovation for tumor-EV-focused research.

Combined Effects of Methylated Cytosine and Molecular Crowding on the Thermodynamic Stability of DNA Duplexes

Methylated cytosine within CpG dinucleotides is a key factor for epigenetic gene regulation. It has been revealed that methylated cytosine decreases DNA backbone flexibility and increases the thermal stability of DNA. Although the molecular environment is an important factor for the structure, thermodynamics, and function of biomolecules, there are few reports on the effects of methylated cytosine under a cell-mimicking molecular environment. Here, we systematically investigated the effects of methylated cytosine on the thermodynamics of DNA duplexes under molecular crowding conditions, which is a critical difference between the molecular environment in cells and test tubes. Thermodynamic parameters quantitatively demonstrated that the methylation effect and molecular crowding effect on DNA duplexes are independent and additive, in which the degree of the stabilization is the sum of the methylation effect and molecular crowding effect. Furthermore, the effects of methylation and molecular crowding correlate with the hydration states of DNA duplexes. The stabilization effect of methylation was due to the favorable enthalpic contribution, suggesting that direct interactions of the methyl group with adjacent bases and adjacent methyl groups play a role in determining the flexibility and thermodynamics of DNA duplexes. These results are useful to predict the properties of DNA duplexes with methylation in cell-mimicking conditions.

Solid-state NMR spectroscopy

Solid-state nuclear magnetic resonance (NMR) spectroscopy is an atomic-level method used to determine the chemical structure, three-dimensional structure, and dynamics of solids and semi-solids. This Primer summarizes the basic principles of NMR as applied to the wide range of solid systems. The fundamental nuclear spin interactions and the effects of magnetic fields and radiofrequency pulses on nuclear spins are the same as in liquid-state NMR. However, because of the anisotropy of the interactions in the solid state, the majority of high-resolution solid-state NMR spectra is measured under magic-angle spinning (MAS), which has profound effects on the types of radiofrequency pulse sequences required to extract structural and dynamical information. We describe the most common MAS NMR experiments and data analysis approaches for investigating biological macromolecules, organic materials, and inorganic solids. Continuing development of sensitivity-enhancement approaches, including 1H-detected fast MAS experiments, dynamic nuclear polarization, and experiments tailored to ultrahigh magnetic fields, is described. We highlight recent applications of solid-state NMR to biological and materials chemistry. The Primer ends with a discussion of current limitations of NMR to study solids, and points to future avenues of development to further enhance the capabilities of this sophisticated spectroscopy for new applications.

Methylation dependent gold adsorption behaviour identifies cancer derived extracellular vesicular DNA

Extracellular vesicles (EV) play a major role in intercellular communication by transmitting cellular materials (e.g. protein, RNA) among distant cells. Recent evidence suggests that they could also contribute to carrying DNA which could inform on the mutational status of the parent tumour DNA. Thus, the fundamental analysis of evDNA could open a better understanding of tumour metastasis and provide new pathways for noninvasive detection and monitoring of cancer. To explore the potential of evDNA for diagnostics, the isolation of pure evDNA from body fluids free of cfDNA contamination is crucial. Herein, we use a liposome based model system to develop an improved evDNA isolation protocol free from cfDNA contamination and evaluate the methylation dependent physicochemical properties of evDNA to develop a simple test for detecting cancer evDNA. Using a highly sensitive multiplex microelectrode device, we demonstrate that serum-evDNA derived from cancer patients show different solution and surface based properties than normal evDNA due to their different methylation landscape (i.e. methylscape). This microdevice allows simultaneous analysis of multiple samples in a single platform from as low as 500 pg μL^-1 of evDNA.

Oxidative damage to epigenetically methylated sites affects DNA stability, dynamics and enzymatic demethylation

DNA damage can affect various regulatory elements of the genome, with the consequences for DNA structure, dynamics, and interaction with proteins remaining largely unexplored. We used solution NMR spectroscopy, restrained and free molecular dynamics to obtain the structures and investigate dominant motions for a set of DNA duplexes containing CpG sites permuted with combinations of 5-methylcytosine (mC), the primary epigenetic base, and 8-oxoguanine (oxoG), an abundant DNA lesion. Guanine oxidation significantly changed the motion in both hemimethylated and fully methylated DNA, increased base pair breathing, induced BI→BII transition in the backbone 3′ to the oxoG and reduced the variability of shift and tilt helical parameters. UV melting experiments corroborated the NMR and molecular dynamics results, showing significant destabilization of all methylated contexts by oxoG. Notably, some dynamic and thermodynamic effects were not additive in the fully methylated oxidized CpG, indicating that the introduced modifications interact with each other. Finally, we show that the presence of oxoG biases the recognition of methylated CpG dinucleotides by ROS1, a plant enzyme involved in epigenetic DNA demethylation, in favor of the oxidized DNA strand. Thus, the conformational and dynamic effects of spurious DNA oxidation in the regulatory CpG dinucleotide can have far-reaching biological consequences.

Effect of Methylation on Local Mechanics and Hydration Structure of DNA

Cytosine methylation affects mechanical properties of DNA and potentially alters the hydration fingerprint for recognition by proteins. The atomistic origin for these effects is not well understood, and we address this via all-atom molecular dynamics simulations. We find that the stiffness of the methylated dinucleotide step changes marginally, whereas the neighboring steps become stiffer. Stiffening is further enhanced for consecutively methylated steps, providing a mechanistic origin for the effect of hypermethylation. Steric interactions between the added methyl groups and the nonpolar groups of the neighboring nucleotides are responsible for the stiffening in most cases. By constructing hydration maps, we found that methylation also alters the surface hydration structure in distinct ways. Its resistance to deformation may contribute to the stiffening of DNA for deformational modes lacking steric interactions. These results highlight the sequence- and deformational-mode-dependent effects of cytosine methylation.

Epigenetically reprogrammed methylation landscape drives the DNA self-assembly and serves as a universal cancer biomarker

Epigenetic reprogramming in cancer genomes creates a distinct methylation landscape encompassing clustered methylation at regulatory regions separated by large intergenic tracks of hypomethylated regions. This methylation landscape that we referred to as Methylscape is displayed by most cancer types, thus may serve as a universal cancer biomarker. To-date most research has focused on the biological consequences of DNA Methylscape changes whereas its impact on DNA physicochemical properties remains unexplored. Herein, we examine the effect of levels and genomic distribution of methylcytosines on the physicochemical properties of DNA to detect the Methylscape biomarker. We find that DNA polymeric behaviour is strongly affected by differential patterning of methylcytosine, leading to fundamental differences in DNA solvation and DNA-gold affinity between cancerous and normal genomes. We exploit these Methylscape differences to develop simple, highly sensitive and selective electrochemical or colorimetric one-step assays for the detection of cancer. These assays are quick, i.e., analysis time ≤10 minutes, and require minimal sample preparation and small DNA input.

DNA methylation is an epigenetic modification that control genetic programs. Here, the authors found that the methylation landscape influences the physicochemical properties of DNA and that it can serve as a universal cancer biomarker, and developed a one-step assay for the detection of cancer DNA.

Of Men and Mice: Modeling the Fragile X Syndrome

The Fragile X Syndrome (FXS) is one of the most common forms of inherited intellectual disability in all human societies. Caused by the transcriptional silencing of a single gene, the fragile x mental retardation gene FMR1, FXS is characterized by a variety of symptoms, which range from mental disabilities to autism and epilepsy. More than 20 years ago, a first animal model was described, the Fmr1 knock-out mouse. Several other models have been developed since then, including conditional knock-out mice, knock-out rats, a zebrafish and a drosophila model. Using these model systems, various targets for potential pharmaceutical treatments have been identified and many treatments have been shown to be efficient in preclinical studies. However, all attempts to turn these findings into a therapy for patients have failed thus far. In this review, I will discuss underlying difficulties and address potential alternatives for our future research.

A sequence environment modulates the impact of methylation on the torsional rigidity of DNA

We explain, with atomic-level detail, how DNA methylation may contribute to the regulation of biological processes. We show that DNA sequence environment modulates the impact of methylation on the physical properties of the double helix.

Computational Modeling Approach in Probing the Effects of Cytosine Methylation on the Transcription Factor Binding to DNA

INTRODUCTION

Cytosine methylation at CpG dinucleotides is a chief mechanism in epigenetic modification of gene expression patterns. Previous studies demonstrated that increased CpG methylation of Sp1 sites at -268 and -346 of protein kinase C ε promoter repressed the gene expression.

MATERIALS & METHODS

The present study investigated the impact of CpG methylation on the Sp1 binding via molecular modeling and electrophoretic mobility shift assay. Each of the Sp1 sites contain two CpGs. Methylation of either CpG lowered the binding affinity of Sp1, whereas methylation of both CpGs produced a greater decrease in the binding affinity. Computation of van der Waals (VDW) energy of Sp1 in complex with the Sp1 sites demonstrated increased VDW values from one to two sites of CpG methylation. Molecular modeling indicated that single CpG methylation caused underwinding of the DNA fragment, with the phosphate groups at C1, C4 and C5 reoriented from their original positions. Methylation of both CpGs pinched the minor groove and increased the helical twist concomitant with a shallow, hydrophobic major groove. Additionally, double methylation eliminated hydrogen bonds on recognition helix residues located at positions -1 and 1, which were essential for interaction with O6/N7 of G-bases. Bonding from linker residues Arg565, Lys595 and Lys596 were also reduced. Methylation of single or both CpGs significantly affected hydrogen bonding from all three Sp1 DNA binding domains, demonstrating that the consequences of cytosine modification extend beyond the neighboring nucleotides.

RESULTS

The results indicate that cytosine methylation causes subtle structural alterations in Sp1 binding sites consequently resulting in inhibition of side chain interactions critical for specific base recognition and reduction of the binding affinity of Sp1.

Unmethylated and methylated CpG dinucleotides distinctively regulate the physical properties of DNA

In eukaryotic cells, DNA has to bend significantly to pack inside the nucleus. Physical properties of DNA such as bending flexibility and curvature are expected to affect DNA packaging and partially determine the nucleosome positioning patterns inside a cell. DNA CpG methylation, the most common epigenetic modification found in DNA, is known to affect the physical properties of DNA. However, its detailed role in nucleosome formation is less well-established. In this study, we evaluated the effect of defined CpG patterns (unmethylated and methylated) on DNA structure and their respective nucleosome-forming ability. Our results suggest that the addition of CpG dinucleotides, either as a (CG)n stretch or (CGX8 )n repeats at 10 bp intervals, lead to reduced hydrodynamic radius and decreased nucleosome-forming ability of DNA. This effect is more predominant for a DNA stretch ((CG)5) located in the middle of a DNA fragment. Methylation of CpG sites, surprisingly, seems to reduce the difference in DNA structure and nucleosome-forming ability among DNA constructs with different CpG patterns. Our results suggest that unmethylated and methylated CpG patterns can play very different roles in regulating the physical properties of DNA. CpG methylation seems to reduce the DNA conformational variations affiliated with defined CpG patterns. Our results can have significant bearings in understanding the nucleosome positioning pattern in living organisms modulated by DNA sequences and epigenetic features.

Effects of cytosine modifications on DNA flexibility and nucleosome mechanical stability

Cytosine can undergo modifications, forming 5-methylcytosine (5-mC) and its oxidized products 5-hydroxymethylcytosine (5-hmC), 5-formylcytosine (5-fC) and 5-carboxylcytosine (5-caC). Despite their importance as epigenetic markers and as central players in cellular processes, it is not well understood how these modifications influence physical properties of DNA and chromatin. Here we report a comprehensive survey of the effect of cytosine modifications on DNA flexibility. We find that even a single copy of 5-fC increases DNA flexibility markedly. 5-mC reduces and 5-hmC enhances flexibility, and 5-caC does not have a measurable effect. Molecular dynamics simulations show that these modifications promote or dampen structural fluctuations, likely through competing effects of base polarity and steric hindrance, without changing the average structure. The increase in DNA flexibility increases the mechanical stability of the nucleosome and vice versa, suggesting a gene regulation mechanism where cytosine modifications change the accessibility of nucleosomal DNA through their effects on DNA flexibility.

Visualizing translocation dynamics and nascent transcript errors in paused RNA polymerases in vivo

BACKGROUND

Transcription elongation is frequently interrupted by pausing signals in DNA, with downstream effects on gene expression. Transcription errors also induce prolonged pausing, which can lead to a destabilized genome by interfering with DNA replication. Mechanisms of pausing associated with translocation blocks and misincorporation have been characterized in vitro, but not in vivo.

RESULTS

We investigate the pausing pattern of RNA polymerase (RNAP) in Escherichia coli by a novel approach, combining native elongating transcript sequencing (NET-seq) with RNase footprinting of the transcripts (RNET-seq). We reveal that the G-dC base pair at the 5' end of the RNA-DNA hybrid interferes with RNAP translocation. The distance between the 5' G-dC base pair and the 3' end of RNA fluctuates over a three-nucleotide width. Thus, the G-dC base pair can induce pausing in post-translocated, pre-translocated, and backtracked states of RNAP. Additionally, a CpG sequence of the template DNA strand spanning the active site of RNAP inhibits elongation and induces G-to-A errors, which leads to backtracking of RNAP. Gre factors efficiently proofread the errors and rescue the backtracked complexes. We also find that pausing events are enriched in the 5' untranslated region and antisense transcription of mRNA genes and are reduced in rRNA genes.

CONCLUSIONS

In E. coli, robust transcriptional pausing involves RNAP interaction with G-dC at the upstream end of the RNA-DNA hybrid, which interferes with translocation. CpG DNA sequences induce transcriptional pausing and G-to-A errors.

Solid-state NMR studies of nucleic acid components

Recent applications of solid-state NMR spectroscopy to studies of nucleic acids and their components.

Dynamics of nucleosome assembly and effects of DNA methylation

The nucleosome is the fundamental packing unit of the eukaryotic genome, and CpG methylation is an epigenetic modification associated with gene repression and silencing. We investigated nucleosome assembly mediated by histone chaperone Nap1 and the effects of CpG methylation based on three-color single molecule FRET measurements, which enabled direct monitoring of histone binding in the context of DNA wrapping. According to our observation, (H3-H4)2 tetramer incorporation must precede H2A-H2B dimer binding, which is independent of DNA termini wrapping. Upon CpG methylation, (H3-H4)2 tetramer incorporation and DNA termini wrapping are facilitated, whereas proper incorporation of H2A-H2B dimers is inhibited. We suggest that these changes are due to rigidified DNA and increased random binding of histones to DNA. According to the results, CpG methylation expedites nucleosome assembly in the presence of abundant DNA and histones, which may help facilitate gene packaging in chromatin. The results also indicate that the slowest steps in nucleosome assembly are DNA termini wrapping and tetramer positioning, both of which are affected heavily by changes in the physical properties of DNA.

Enhanced sampling simulations of DNA step parameters

A novel approach for the selection of step parameters as reaction coordinates in enhanced sampling simulations of DNA is presented. The method uses three atoms per base and does not require coordinate overlays or idealized base pairs. This allowed for a highly efficient implementation of the calculation of all step parameters and their Cartesian derivatives in molecular dynamics simulations. Good correlation between the calculated and actual twist, roll, tilt, shift, and slide parameters is obtained, while the correlation with rise is modest. The method is illustrated by its application to the methylated and unmethylated 5'-CATGTGACGTCACATG-3' double stranded DNA sequence. One-dimensional umbrella simulations indicate that the flexibility of the central CG step is only marginally affected by methylation.

DNA methylation reduces binding and cleavage by bleomycin

In a recent study, we described the enhanced double-strand cleavage of hairpin DNAs by Fe·bleomycin (Fe·BLM) that accompanies increasingly strong binding of this antitumor agent and suggested that this effect may be relevant to the mechanism by which BLM mediates its antitumor effects. Because the DNA in tumor cells is known to be hypomethylated on cytidine relative to that in normal cells, it seemed of interest to study the possible effects of methylation status on BLM-induced double-strand DNA cleavage. Three hairpin DNAs found to bind strongly to bleomycin, and their methylated counterparts, were used to study the effect of methylation on bleomycin-induced DNA degradation. Under conditions of limited DNA cleavage, there was a significant overall decrease in the cleavage of methylated hairpin DNAs. Cytidine methylation was found to result in decreased BLM-induced cleavage at the site of methylation and to result in enhanced cleavage at adjacent nonmethylated sites. For two of the three hairpin DNAs studied, methylation was accompanied by a dramatic decrease in the binding affinity for Fe·BLM, suggesting the likelihood of diminished double-strand cleavage. The source of the persistent binding of BLM by the third hairpin DNA was identified. Also identified was the probable molecular mechanism for diminished binding and cleavage of the methylated DNAs by BLM. The possible implications of these findings for the antitumor selectivity of bleomycin are discussed.

Coupled binding-bending-folding: The complex conformational dynamics of protein-DNA binding studied by atomistic molecular dynamics simulations

BACKGROUND

Protein-DNA binding often involves dramatic conformational changes such as protein folding and DNA bending. While thermodynamic aspects of this behavior are understood, and its biological function is often known, the mechanism by which the conformational changes occur is generally unclear. By providing detailed structural and energetic data, molecular dynamics simulations have been helpful in elucidating and rationalizing protein-DNA binding.

SCOPE OF REVIEW

This review will summarize recent atomistic molecular dynamics simulations of the conformational dynamics of DNA and protein-DNA binding. A brief overview of recent developments in DNA force fields is given as well.

MAJOR CONCLUSIONS

Simulations have been crucial in rationalizing the intrinsic flexibility of DNA, and have been instrumental in identifying the sequence of binding events, the triggers for the conformational motion, and the mechanism of binding for a number of important DNA-binding proteins.

GENERAL SIGNIFICANCE

Molecular dynamics simulations are an important tool for understanding the complex binding behavior of DNA-binding proteins. With recent advances in force fields and rapid increases in simulation time scales, simulations will become even more important for future studies. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Modeling biominerals formed by apatites and DNA

Different aspects of biominerals formed by apatite and DNA have been investigated using computer modeling tools. Firstly, the structure and stability of biominerals in which DNA molecules are embedded into hydroxyapatite and fluoroapatite nanopores have been examined by combining different molecular mechanics methods. After this, the early processes in the nucleation of hydroxyapatite at a DNA template have been investigated using molecular dynamics simulations. Results indicate that duplexes of DNA adopting a B double helix can be encapsulated inside nanopores of hydroxyapatite without undergoing significant distortions in the inter-strand hydrogen bonds and the intra-strand stacking. This ability of hydroxyapatite is practically independent of the DNA sequence, which has been attributed to the stabilizing role of the interactions between the calcium atoms of the mineral and the phosphate groups of the biomolecule. In contrast, the fluorine atoms of fluoroapatite induce pronounced structural distortions in the double helix when embedded in a pore of the same dimensions, resulting in the loss of its most relevant characteristics. On the other hand, molecular dynamics simulations have allowed us to observe the formation of calcium phosphate clusters at the surface of the B-DNA template. Electrostatic interactions between the phosphate groups of DNA and Ca(2+) have been found to essential for the formation of stable ion complexes, which were the starting point of calcium phosphate clusters by incorporating PO3(4) from the solution.

Effects of 2'-O-methyl nucleotide substitution on EcoRI endonuclease cleavage activities

To investigate the effect of sugar pucker conformation on DNA-protein interactions, we used 2′-O-methyl nucleotide (2′-OMeN) to modify the EcoRI recognition sequence -TGAATTCT-, and monitored the enzymatic cleavage process using FRET method. The 2′-O-methyl nucleotide has a C3′-endo sugar pucker conformation different from the C2′-endo sugar pucker conformation of native DNA nucleotides. The initial reaction velocities were measured and the kinetic parameters, K_m and V_max were derived using Michaelis-Menten equation. Experimental results showed that 2′-OMeN substitutions for the EcoRI recognition sequence decreased the cleavage efficiency for A2, A3 and T4 substitutions significantly, and 2′-OMeN substitution for T5 residue inhibited the enzymatic activity completely. In contrast, substitutions for G1 and C6 could maintain the original activity. 2′-fluoro nucleic acid (2′-FNA) and locked nucleic acid (LNA) having similar C3′-endo sugar pucker conformation also demonstrated similar enzymatic results. This position-dependent enzymatic cleavage property might be attributed to the phosphate backbone distortion caused by the switch from C2′-endo to C3′-endo sugar pucker conformation, and was interpreted on the basis of the DNA-EcoRI structure. These 2′-modified nucleotides could behave as a regulatory element to modulate the enzymatic activity in vitro, and this property will have potential applications in genetic engineering and biomedicine.

NMR structure of 2'-O-(2-methoxyethyl) modified and C5-methylated RNA dodecamer duplex

The solution-state structure of 2'-O-(2-methoxyethly) substituted dodecamer r(*CG*CGAA*U*U*CG*C)d(G), 2'-MOE RNA, with all cytosines and uracils methylated at the C5-position has been determined by NMR spectroscopy. The chemical modifications were used to improve the oligonucleotide's drug-like properties. The 2'-MOE group drives pseudorotational equilibrium of the ribofuranose moiety to the N-type conformation and supposedly results in structural preorganization leading to high affinity of a modified oligonucleotide towards its complementary biological target, improved pharmacokinetic and toxicological properties. The high melting temperature of the antiparallel duplex structure adopted by 2'-MOE RNA was explained through the formation of a stable A-form RNA consistent with effective base-pairing and stacking interactions. The comparison of the solution-state structure with the crystal structure of a non-methylated analogue shows an increase in the stacking at the base pair steps for the C5-methylated 2'-MOE RNA duplex. The MOE substituents adopt a well-defined structure in the minor groove with the predominant gauche conformations around the ethylene bond.

Most methylation-susceptible DNA sequences in human embryonic stem cells undergo a change in conformation or flexibility upon methylation

DNA methylation in eukaryotes occurs on the cytosine bases in CG, CHG, and CHH (where H indicates non-G nucleotides) contexts and provides an important epigenetic mark in various biological processes. However, the structural and physical properties of methylated DNA are poorly understood. Using nondenaturing polyacrylamide gel electrophoresis, we performed a systematic study of the influence of DNA methylation on the conformation and physical properties of DNA for all CG, CHG, and CHH contexts. In the CG context, methylated multimers of the CG/CG-containing unit fragment migrated in gels slightly faster than their unmethylated counterparts. In the CHG context, both homo- and hemimethylation caused retarded migration of multimers of the CAG/CTG-containing fragment. In the CHH context, methylation caused or enhanced retarded migration of the multimers of CAA/TTG-, CAT/ATG-, CAC/GTG-, CTA/TAG-, or CTT/AAG-containing fragments. These results suggest that methylation increases DNA rigidity in the CG context and introduces distortions into several CHG and CHH sequences. More interestingly, we found that nearly all of the methylation repertoires in the CHG context and 98% of those in the CHH context in human embryonic stem cells were species that undergo conformational changes upon methylation. Similarly, most of the methylation repertoires in the Arabidopsis CHG and CHH contexts were sequences with methylation-induced distortion. We hypothesize that the methylation-induced properties or conformational changes in DNA may facilitate nucleosome formation, which provides the essential mechanism for alterations of chromatin density.

Conformation effects of CpG methylation on single-stranded DNA oligonucleotides: analysis of the opioid peptide dynorphin-coding sequences

Single-stranded DNA (ssDNA) is characterized by high conformational flexibility that allows these molecules to adopt a variety of conformations. Here we used native polyacrylamide gel electrophoresis (PAGE), circular dichroism (CD) spectroscopy and nuclear magnetic resonance (NMR) spectroscopy to show that cytosine methylation at CpG sites affects the conformational flexibility of short ssDNA molecules. The CpG containing 37-nucleotide PDYN (prodynorphin) fragments were used as model molecules. The presence of secondary DNA structures was evident from differences in oligonucleotide mobilities on PAGE, from CD spectra, and from formation of A-T, G-C, and non-canonical G-T base pairs observed by NMR spectroscopy. The oligonucleotides displayed secondary structures at 4°C, and some also at 37°C. Methylation at CpG sites prompted sequence-dependent formation of novel conformations, or shifted the equilibrium between different existing ssDNA conformations. The effects of methylation on gel mobility and base pairing were comparable in strength to the effects induced by point mutations in the DNA sequences. The conformational effects of methylation may be relevant for epigenetic regulatory events in a chromatin context, including DNA-protein or DNA-DNA recognition in the course of gene transcription, and DNA replication and recombination when double-stranded DNA is unwinded to ssDNA.

Effects of histone acetylation and CpG methylation on the structure of nucleosomes

Nucleosomes are the fundamental packing units of the eukaryotic genome. A nucleosome core particle comprises an octameric histone core wrapped around by ~147bp DNA. Histones and DNA are targets for covalent modifications mediated by various chromatin modification enzymes. These modifications play crucial roles in various gene regulation activities. A group of common hypotheses for the mechanisms of gene regulation involves changes in the structure and structural dynamics of chromatin induced by chromatin modifications. We employed single molecule fluorescence methods to test these hypotheses by monitoring the structure and structural dynamics of nucleosomes before and after histone acetylation and DNA methylation, two of the best-conserved chromatin modifications throughout eukaryotes. Our studies revealed that these modifications induce changes in the structure and structural dynamics of nucleosomes that may contribute directly to the formation of open or repressive chromatin conformation.

Effects of DNA methylation on the structure of nucleosomes

Nucleosomes are the fundamental packing units of the eukaryotic genome. Understanding the dynamic structure of a nucleosome is a key to the elucidation of genome packaging in eukaryotes, which is tied to the mechanisms of gene regulation. CpG methylation of DNA is an epigenetic modification associated with the inactivation of transcription and the formation of a repressive chromatin structure. Unraveling the changes in the structure of nucleosomes upon CpG methylation is an essential step toward the understanding of the mechanisms of gene repression and silencing by CpG methylation. Here we report single-molecule and ensemble fluorescence studies showing how the structure of a nucleosome is affected by CpG methylation. The results indicate that CpG methylation induces tighter wrapping of DNA around the histone core accompanied by a topology change. These findings suggest that changes in the physical properties of nucleosomes induced upon CpG methylation may contribute directly to the formation of a repressive chromatin structure.

Cytosine methylation alters DNA mechanical properties

DNA methylation plays an essential role in transcriptional control of organismal development in epigenetics, from turning off a specific gene to inactivation of entire chromosomes. While the biological function of DNA methylation is becoming increasingly clear, the mechanism of methylation-induced gene regulation is still poorly understood. Through single-molecule force experiments and simulation we investigated the effects of methylation on strand separation of DNA, a crucial step in gene expression. Molecular force assay and single-molecule force spectroscopy revealed a strong methylation dependence of strand separation. Methylation is observed to either inhibit or facilitate strand separation, depending on methylation level and sequence context. Molecular dynamics simulations provided a detailed view of methylation effects on strand separation, suggesting the underlying physical mechanism. According to our study, methylation in epigenetics may regulate gene expression not only through mechanisms already known but also through changing mechanical properties of DNA.

Solid state 2H NMR analysis of furanose ring dynamics in DNA containing uracil

DNA damage has been implicated in numerous human diseases, particularly cancer, and the aging process. Single-base lesions, such as uracil, in DNA can be cytotoxic or mutagenic and are recognized by a DNA glycosylase during the process of base excision repair. Increased dynamic properties in lesion-containing DNAs have been suggested to assist recognition and specificity. Deuterium solid-state nuclear magnetic resonance (SSNMR) has been used to directly observe local dynamics of the furanose ring within a uracil:adenine (U:A) base pair and compared to a normal thymine:adenine (T:A) base pair. Quadrupole echo lineshapes, <T(1Z)>, and <T(2e)> relaxation data were collected, and computer modeling was performed. The results indicate that the relaxation times are identical within the experimental error, the solid lineshapes are essentially indistinguishable above the noise level, and our lineshapes are best fit with a model that does not have significant local motions. Therefore, U:A base pair furanose rings appear to have essentially identical dynamic properties as a normal T:A base pair, and the local dynamics of the furanose ring are unlikely to be the sole arbiter for uracil recognition and specificity in U:A base pairs.

Site-specific DNA structural and dynamic features revealed by nucleotide-independent nitroxide probes

In site-directed spin labeling, a covalently attached nitroxide probe containing a chemically inert unpaired electron is utilized to obtain information on the local environment of the parent macromolecule. Studies presented here examine the feasibility of probing local DNA structural and dynamic features using a class of nitroxide probes that are linked to chemically substituted phosphorothioate positions at the DNA backbone. Two members of this family, designated as R5 and R5a, were attached to eight different sites of a dodecameric DNA duplex without severely perturbing the native B-form conformation. Measured X-band electron paramagnetic resonance (EPR) spectra, which report on nitroxide rotational motions, were found to vary depending on the location of the label (e.g., duplex center vs termini) and the surrounding DNA sequence. This indicates that R5 and R5a can provide information on the DNA local environment at the level of an individual nucleotide. As these probes can be attached to arbitrary nucleotides within a nucleic acid sequence, they may provide a means to "scan" a given DNA molecule in order to interrogate its local structural and dynamic features.

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.

Conformational dynamics of an intact virus: order parameters for the coat protein of Pf1 bacteriophage

This study has examined the atomic-level dynamics of the protein in the capsid of filamentous phage Pf1. This capsid consists of approximately 7,300 small subunits of only 46 aa in a helical array around a highly extended, circular single-stranded DNA molecule of 7,349 nt. Measurements were made of site-specific, solid-state NMR order parameters, S, the values which are dimensionless quantities between 0 (mobile) and 1 (static) that characterize the amplitudes of molecular bond angular motions that are faster than microseconds. It was found that the protein subunit backbone is very static, and of particular interest, it appears to be static at residues glycine 15 and glutamine 16 where it had been previously thought to be mobile. In contrast to the backbone, several side chains display large-amplitude angular motions. Side chains on the virion exterior that interact with solvent are highly mobile, but surprisingly, the side chains of residues arginine 44 and lysine 45 near the DNA deep in the interior of the virion are also highly mobile. The large-amplitude dynamic motion of these positively charged side chains in their interactions with the DNA were not previously expected. The results reveal a highly dynamic aspect of a DNA-protein interface within a virus.

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.

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.

The N6-Methyl Group of Adenine Further Increases the BI Stability of DNA Compared to C5-Methyl Groups

Methylated DNA bases are natural modifications which play an important role in protein-DNA interactions. Recent experimental and theoretical results have shown an influence of the base modification on the conformational behavior of the DNA backbone. MD simulations of four different B-DNA dodecamers (d(GC)(6), d(AT)(6), d(G(5mCG)(5)C), and d(A(T6mA)(5)T)) have been performed with the aim to examine the influence of methyl groups on the B-DNA backbone behavior. An additional control simulation of d(AU)(6) has also been performed to examine the further influence of the C5-methyl group in thymine. Methyl groups in the major groove (as in C5-methylcytosine, thymine, or N6-methyladenine) decrease the BII substate population of RpY steps. Due to methylation a clearer distinction of the BI substate stability between YpR and RpY (CpG/GpC or TpA/ApT) steps arises. A positive correlation between the BII substate population and base stacking distances is seen only for poly(GC). A methyl group added into the major groove increases mean water residence times around the purine N7 atom, which may stabilize the BI substate by improving the hydration network between the DNA backbone and the major groove. The N6-methyl group also forms a water molecule bridge between the N6 and O4 atoms, and thus further stabilizes the BI substate.

Recent developments in solid-state magic-angle spinning, nuclear magnetic resonance of fully and significantly isotopically labelled peptides and proteins

In recent years, a large number of solid-state nuclear magnetic resonance (NMR) techniques have been developed and applied to the study of fully or significantly isotopically labelled ((13)C, (15)N or (13)C/(15)N) biomolecules. In the past few years, the first structures of (13)C/(15)N-labelled peptides, Gly-Ile and Met-Leu-Phe, and a protein, Src-homology 3 domain, were solved using magic-angle spinning NMR, without recourse to any structural information obtained from other methods. This progress has been made possible by the development of NMR experiments to assign solid-state spectra and experiments to extract distance and orientational information. Another key aspect to the success of solid-state NMR is the advances made in sample preparation. These improvements will be reviewed in this contribution. Future prospects for the application of solid-state NMR to interesting biological questions will also briefly be discussed.

A determining influence for CpG dinucleotides on nucleosome positioning in vitro

DNA sequence information that directs the translational positioning of nucleosomes can be attenuated by cytosine methylation when a short run of CpG dinucleotides is located close to the dyad axis of the nucleosome. Here, we show that point mutations introduced to re-pattern methylation at the (CpG)3 element in the chicken betaA-globin promoter sequence themselves strongly influenced nucleosome formation in reconstituted chromatin. The disruptive effect of cytosine methylation on nucleosome formation was found to be determined by the sequence context of CpG dinucleotides, not just their location in the positioning sequence. Additional mutations indicated that methylation can also promote the occupation of certain nucleosome positions. DNase I analysis demonstrated that these genetic and epigenetic modifications altered the structural characteristics of the (CpG)3 element. Our findings support a proposal that the intrinsic structural properties of the DNA at the -1.5 site, as occupied by (CpG)3 in the nucleosome studied, can be decisive for nucleosome formation and stability, and that changes in anisotropic DNA bending or flexibility at this site explain why nucleosome positioning can be exquisitely sensitive to genetic and epigenetic modification of the DNA sequence.

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.

C5-Methylation of Cytosine in B-DNA Thermodynamically and Kinetically Stabilizes BI

We have examined the backbone dynamics of two alternating purine-pyrimidine dodecamers. One sequence consists of "pure" GC bases; the other one contains 5-methylcytosines. The effect of the methyl groups on the backbone substates BI/BII was investigated by means of molecular dynamics. The methylation influences, on one hand, the transition barrier between BI and BII and, on the other hand, the state of equilibrium. The kinetic consequences are an increase of the DeltaG of Gp5mC steps by 1.5 kcal/mol and a decrease of the DeltaG of 5mCpG steps by 0.8 kcal/mol (compared with the nonmethylated DNA). Thus, the additive group differentiates between the two occurring dinucleotide steps and renders the phosphate of the 5-methylcytosine more rigid, as proposed by experimental studies. The thermodynamic consequences are an increase of the DeltaG of Gp5mC steps by 1.1 kcal/mol and a decrease of the DeltaG of 5mCpG steps by 0.8 kcal/mol. The reason for this shift in equilibrium is still not completely clear on a molecular basis. But we can conclude that the indirect readout of DNA is influenced by methylation.

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.

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.

Sequence-dependent dynamics of duplex DNA: the applicability of a dinucleotide model

The short-time (submicrosecond) bending dynamics of duplex DNA were measured to determine the effect of sequence on dynamics. All measurements were obtained from a single site on duplex DNA, using a single, site-specific modified base containing a rigidly tethered, electron paramagnetic resonance active spin probe. The observed dynamics are interpreted in terms of single-step sequence-dependent bending force constants, determined from the mean squared amplitude of bending relative to the end-to-end vector using the modified weakly bending rod model. The bending dynamics at a single site are a function of the sequence of the nucleotides constituting the duplex DNA. We developed and examined several dinucleotide-based models for flexibility. The models indicate that the dominant feature of the dynamics is best explained in terms of purine- and pyrimidine-type steps, although distinction is made among all 10 unique steps: It was found that purine-purine steps (which are the same as pyrimidine-pyrimidine steps) were near average in flexibility, but the pyrimidine-purine steps (5' to 3') were nearly twice as flexible, whereas purine-pyrimidine steps were more than half as flexible as average DNA. Therefore, the range of stepwise flexibility is approximately fourfold and is characterized by both the type of base pair step (pyrimidine/purine combination) and the identity of the bases within the pair (G, A, T, or C). All of the four models considered here underscore the complexity of the dependence of dynamics on DNA sequence with certain sequences not satisfactorily explainable in terms of any dinucleotide model. These findings provide a quantitative basis for interpreting the dynamics and kinetics of DNA-sequence-dependent biological processes, including protein recognition and chromatin packaging.

Bending and flexibility of methylated and unmethylated EcoRI DNA

We used cyclization kinetics experiments and Monte Carlo simulations to determine a structural model for a DNA decamer containing the EcoRI restriction site. Our findings agree well with recent crystal and NMR structures of the EcoRI dodecamer, where an overall bend of seven degrees is distributed symmetrically over the molecule. Monte Carlo simulations indicate that the sequence has a higher flexibility, assumed to be isotropic, compared to that of a "generic" DNA sequence. This model was used as a starting point for the investigation of the effect of cytosine methylation on DNA bending and flexibility. While methylation did not affect bend magnitude or direction, it resulted in a reduction in bending flexibility and under-winding of the methylated nucleotides. We demonstrate that our approach can augment the understanding of DNA structure and dynamics by adding information about the global structure and flexibility of the sequence. We also show that cyclization kinetics can be used to study the properties of modified nucleotides.

Modeling furanose ring dynamics in DNA

Determination of the conformational flexibility of the furanose ring is of vital importance in understanding the structure of DNA. In this work we have applied a model of furanose ring motion to the analysis of deuterium line shape data obtained from sugar rings in solid hydrated DNA. The model describes the angular trajectories of the atoms in the furanose ring in terms of pseudorotation puckering amplitude (q) and the pseudorotation puckering phase phi. Fixing q, the motion is thus treated as Brownian diffusion through an angular-dependent potential U(phi). We have simulated numerous line shapes varying the adjustable parameters, including the diffusion coefficient D, pseudorotation puckering amplitude q, and the form of the potential U(phi). We have used several forms of the potential, including equal double-well potentials, unequal double-well potentials, and a potential truncated to "second order" in the Fourier series. To date, we have obtained best simulations for both equilibrium and nonequilibrium (partially relaxed) solid-state deuterium NMR line shapes for the sample [2' '-2H]-2'-deoxycytidine at the position C3 (underlined) in the DNA sequence [d(CGCGAATTCGCG)]2, using a double-well potential with an equal barrier height of U(0) = 5.5k(B)T ( approximately 3.3 kcal/mol), a puckering amplitude of q = 0.4 A, and a diffusion coefficient characterizing the underlying stochastic jump rate D = 9.9 x 10(8) Hz. Then the rate of flux for the C-D bond over the barrier, i.e., the escape velocity or the overall rate of puckering between modes, was found to be 0.7 x 10(7) Hz.