Spontaneous base flipping in DNA and its possible role in methyltransferase binding

Adenine Methylation Enhances the Conformational Flexibility of an RNA Hairpin Tetraloop

The N6-methyladenosine modification is one of the most abundant post-transcriptional modifications in ribonucleic acid (RNA) molecules. Using molecular dynamics simulations and alchemical free-energy calculations, we studied the structural and energetic implications of incorporating this modification in an adenine mononucleotide and an RNA hairpin structure. At the mononucleotide level, we found that the syn configuration is more favorable than the anti configuration by 2.05 ± 0.15 kcal/mol. The unfavorable effect of methylation was due to the steric overlap between the methyl group and a nitrogen atom in the purine ring. We then probed the effect of methylation in an RNA hairpin structure containing an AUCG tetraloop, which is recognized by a "reader" protein (YTHDC1) to promote transcriptional silencing of long noncoding RNAs. While methylation had no significant conformational effect on the hairpin stem, the methylated tetraloop showed enhanced conformational flexibility compared to the unmethylated tetraloop. The increased flexibility was associated with the outward flipping of two bases (A6 and U7) which formed stacking interactions with each other and with the C8 and G9 bases in the tetraloop, leading to a conformation similar to that in the RNA/reader protein complex. Therefore, methylation-induced conformational flexibility likely facilitates RNA recognition by the reader protein.

Base Dynamics in the HhaI Protein Binding Site

Protein-DNA interactions play an important role in numerous biological functions within the living cell. In many of these interactions, the DNA helix is significantly distorted upon protein-DNA complex formation. The HhaI restriction-modification system is one such system, where the methylation target is flipped out of the helix when bound to the methyltransferase. However, the base flipping mechanism is not well understood. The dynamics of the binding site of the HhaI methyltransferase and endonuclease (underlined) within the DNA oligomer [d(G1A2T3A4G5C6G7C8T9A10T11C12)]2 are studied using deuterium solid-state NMR (SSNMR). SSNMR spectra obtained from DNAs deuterated on the base of nucleotides within and flanking the [5'-GCGC-3']2 sequence indicate that all of these positions are structurally flexible. Previously, conformational flexibility within the phosphodiester backbone and furanose ring within the target sequence has been observed and hypothesized to play a role in the distortion mechanism. However, whether that distortion was occurring through an active or passive mechanism remained unclear. These NMR data demonstrate that although the [5'-GCGC-3']2 sequence is dynamic, the target cytosine is not passively flipping out of the double-helix on the millisecond-picosecond time scale. Additionally, although previous studies have shown that both the furanose ring and phosphodiester backbone experience a change in dynamics upon methylation, which may play a role in recognition and cleavage by the endonuclease, our observations here indicate that methylation has no effect on the dynamics of the base itself.

Energy Landscapes for Base-Flipping in a Model DNA Duplex

We explore the process of base-flipping for four central bases, adenine, guanine, cytosine, and thymine, in a deoxyribonucleic acid (DNA) duplex using the energy landscape perspective. NMR imino-proton exchange and fluorescence correlation spectroscopy studies have been used in previous experiments to obtain lifetimes for bases in paired and extrahelical states. However, the difference of almost 4 orders of magnitude in the base-flipping rates obtained by the two methods implies that they are exploring different pathways and possibly different open states. Our results support the previous suggestion that minor groove opening may be favored by distortions in the DNA backbone and reveal links between sequence effects and the direction of opening, i.e., whether the base flips toward the major or the minor groove side. In particular, base flipping along the minor groove pathway was found to align toward the 5' side of the backbone. We find that bases align toward the 3' side of the backbone when flipping along the major groove pathway. However, in some cases for cytosine and thymine, the base flipping along the major groove pathway also aligns toward the 5' side. The sequence effect may be caused by the polar interactions between the flipping-base and its neighboring bases on either of the strands. For guanine flipping toward the minor groove side, we find that the equilibrium constant for opening is large compared to flipping via the major groove. We find that the estimated rates of base opening, and hence the lifetimes of the closed state, obtained for thymine flipping through small and large angles along the major groove differ by 6 orders of magnitude, whereas for thymine flipping through small angles along the minor groove and large angles along the major groove, the rates differ by 3 orders of magnitude.

Dynamics of 5R-Tg Base Flipping in DNA Duplexes Based on Simulations─Agreement with Experiments and Beyond

Damaged or mismatched DNA bases are normally thought to be able to flip out of the helical stack, providing enzymes with access to the faulty genetic information otherwise hidden inside the helix. Thymine glycol (Tg) is one of the most common products of nucleic acid damage. However, the static and dynamic structures of DNA duplexes affected by 5R-Tg epimers are still not clearly understood, including the ability of these to undergo spontaneous base flipping. Structural effects of the 5R-Tg epimers on the duplex DNA are herein studied using molecular dynamics together with reliable DFT based calculations. In comparison with the corresponding intact DNA, the cis-5R,6S-Tg epimer base causes little perturbation to the duplex DNA, and a barrier of 4.9 kcal mol^–1 is obtained by meta-eABF for cis-5R,6S-Tg base flipping out of the duplex DNA, comparable to the 5.4 kcal mol^–1 obtained for the corresponding thymine flipping in intact DNA. For the trans-5R,6R-Tg epimer, three stable local structures were identified, of which the most stable disrupts the Watson–Crick hydrogen-bonded G5/C20 base pair, leading to conformational distortion of the duplex. Interestingly, the relative barrier height of the 5R-Tg flipping is only 1.0 kcal mol^–1 for one of these trans-5R,6R-Tg epimers. Water bridge interactions were identified to be essential for 5R-Tg flipping. The study clearly demonstrates the occurrence of partial trans-5R,6R-Tg epimer flipping in solution.

Reaction Coordinate and Thermodynamics of Base Flipping in RNA

Base flipping is a key biophysical event involved in recognition of various ligands by ribonucleic acid (RNA) molecules. However, the mechanism of base flipping in RNA remains poorly understood, in part due to the lack of atomistic details on complex rearrangements in neighboring bases. In this work, we applied transition path sampling (TPS) methods to study base flipping in a double-stranded RNA (dsRNA) molecule that is known to interact with RNA-editing enzymes through this mechanism. We obtained an ensemble of 1000 transition trajectories to describe the base-flipping process. We used the likelihood maximization method to determine the refined reaction coordinate (RC) consisting of two collective variables (CVs), a distance and a dihedral angle between nucleotides that form stacking interactions with the flipping base. The free energy profile projected along the refined RC revealed three minima, two corresponding to the initial and final states and one for a metastable state. We suggest that the metastable state likely represents a wobbled conformation of nucleobases observed in NMR studies that is often characterized as the flipped state. The analyses of reactive trajectories further revealed that the base flipping is coupled to a global conformational change in a stem-loop of dsRNA.

Dynamics of spontaneous flipping of a mismatched base in DNA duplex

DNA base flipping is a fundamental theme in DNA biophysics. The dynamics for a B-DNA base to spontaneously flip out of the double helix has significant implications in various DNA-protein interactions but are still poorly understood. The spontaneous base-flipping rate obtained previously via the imino proton exchange assay is most likely the rate of base wobbling instead of flipping. Using the diffusion-decelerated fluorescence correlation spectroscopy together with molecular dynamics simulations, we show that a base of a single mismatched base pair (T-G, T-T, or T-C) in a double-stranded DNA can spontaneously flip out of the DNA duplex. The extrahelical lifetimes are on the order of 10 ms, whereas the intrahelical lifetimes range from 0.3 to 20 s depending on the stability of the base pairs. These findings provide detailed understanding on the dynamics of DNA base flipping and lay down foundation to fully understand how exactly the repair proteins search and locate the target mismatched base among a vast excess of matched DNA bases.

AN OVERVIEW OF THE STRETCHED INTERMEDIATE MODEL OF B–Z DNA TRANSITION

Recently, the stretched intermediate model was proposed for the B–Z deoxyribonucleic acid (DNA) transition based on simulation results carried out using the Stochastic Difference Equation (SDE) that showed unwinding and elongation of the oligomer during the transition. This model has proven to be successful in resolving the steric dilemma in short oligomers. However, extending the simulation method to larger DNA strands may prove to be computationally challenging. Such difficulty has led us to adopt a mean field approach using phenomenological interaction potentials to simulate the transition. Like the atomistic approach, the SDE simulations based on the mean field approach, also suggest the presence of a stretched intermediate during the transition.

Conformational diversity of short DNA duplex

Two 25 base-pair cDNA strands are encapsulated within an optically trapped nanodroplet, and we observe the kinetics of their hybridization in dynamic equilibrium via single-molecule fluorescence resonance energy transfer (FRET) as a function of temperature and of the solution's NaCl concentration. We have observed the duplex unfolding and refolding, and we have observed quasistable partially unfolded states under low salinity conditions. Furthermore, our measurements reveal that, even in conditions under which the duplex is stable, it undergoes conformational fluctuations in solution.

Opening of DNA double helix at room temperature: Application of alpha-cyclodextrin self-aggregates

Self-aggregation of alpha-cyclodextrin (alpha-CD) can induce DNA opening at room temperature (25 degrees C) owing to the hydroxyl groups on the surface of the spherical aggregates of alpha-CD, which promote hydrogen bonding with the flipped-out bases in DNA duplex prohibiting them from reverting back.

Backbone dynamics in the DNA HhaI protein binding site

The dynamics of the phosphodiester backbone in the [5'-GCGC-3'] 2 moiety of the DNA oligomer [d(G 1A 2T 3A 4 G 5 C 6 G 7 C 8T 9A 10T 11C 12)] 2 are studied using deuterium solid-state NMR (SSNMR). SSNMR spectra obtained from DNAs nonstereospecifically deuterated on the 5' methylene group of nucleotides within the [5'-GCGC-3'] 2 moiety indicated that all of these positions are structurally flexible. Previous work has shown that methylation reduces the amplitude of motion in the phosphodiester backbone and furanose ring of the same DNA, and our observations indicate that methylation perturbs backbone dynamics through not only a loss of mobility but also a change of direction of motion. These NMR data indicate that the [5'-GCGC-3'] 2 moiety is dynamic, with the largest amplitude motions occurring nearest the methylation site. The change of orientation of this moiety in DNA upon methylation may make the molecule less amenable to binding to the HhaI endonuclease.

The role of Arg165 towards base flipping, base stabilization and catalysis in M.HhaI

Arg165 forms part of a previously identified base flipping motif in the bacterial DNA cytosine methyltransferase, M.HhaI. Replacement of Arg165 with Ala has no detectable effect on either DNA or AdoMet affinity, yet causes the base flipping and restacking transitions to be decreased approximately 16 and 190-fold respectively, thus confirming the importance of this motif. However, these kinetic changes cannot account for the mutant's observed 10(5)-fold decreased catalytic rate. The mutant enzyme/cognate DNA cocrystal structure (2.79 A resolution) shows the target cytosine to be positioned approximately 30 degrees into the major groove, which is consistent with a major groove pathway for nucleotide flipping. The pyrimidine-sugar chi angle is rotated to approximately +171 degrees, from a range of -95 degrees to -120 degrees in B DNA, and -77 degrees in the WT M.HhaI complex. Thus, Arg165 is important for maintaining the cytosine positioned for nucleophilic attack by Cys81. The cytosine sugar pucker is in the C2'-endo-C3'-exo (South conformation), in contrast to the previously reported C3'-endo (North conformation) described for the original 2.70 A resolution cocrystal structure of the WT M.HhaI/DNA complex. We determined a high resolution structure of the WT M.HhaI/DNA complex (1.96 A) to better determine the sugar pucker. This new structure is similar to the original, lower resolution WT M.HhaI complex, but shows that the sugar pucker is O4'-endo (East conformation), intermediate between the South and North conformers. In summary, Arg165 plays significant roles in base flipping, cytosine positioning, and catalysis. Furthermore, the previously proposed M.HhaI-mediated changes in sugar pucker may not be an important contributor to the base flipping mechanism. These results provide insights into the base flipping and catalytic mechanisms for bacterial and eukaryotic DNA methyltransferases.

DNA bending by M.EcoKI methyltransferase is coupled to nucleotide flipping

The maintenance methyltransferase M.EcoKI recognizes the bipartite DNA sequence 5′-AACNNNNNNGTGC-3′, where N is any nucleotide. M.EcoKI preferentially methylates a sequence already containing a methylated adenine at or complementary to the underlined bases in the sequence. We find that the introduction of a single-stranded gap in the middle of the non-specific spacer, of up to 4 nt in length, does not reduce the binding affinity of M.EcoKI despite the removal of non-sequence-specific contacts between the protein and the DNA phosphate backbone. Surprisingly, binding affinity is enhanced in a manner predicted by simple polymer models of DNA flexibility. However, the activity of the enzyme declines to zero once the single-stranded region reaches 4 nt in length. This indicates that the recognition of methylation of the DNA is communicated between the two methylation targets not only through the protein structure but also through the DNA structure. Furthermore, methylation recognition requires base flipping in which the bases targeted for methylation are swung out of the DNA helix into the enzyme. By using 2-aminopurine fluorescence as the base flipping probe we find that, although flipping occurs for the intact duplex, no flipping is observed upon introduction of a gap. Our data and polymer model indicate that M.EcoKI bends the non-specific spacer and that the energy stored in a double-stranded bend is utilized to force or flip out the bases. This energy is not stored in gapped duplexes. In this way, M.EcoKI can determine the methylation status of two adenine bases separated by a considerable distance in double-stranded DNA and select the required enzymatic response.

The stretched intermediate model of B-Z DNA transition

There have been numerous attempts to describe the mechanism of B-Z transition. Our simulations based on the stochastic difference equation with length algorithm show that a short DNA oligomer will tend to unwind and overstretch during the transition. The overstretching of DNA is then understood from the Zhou, Zhang, and Ou-Yang model. Unlike the Harvey model, the stretched intermediate model does not pose any steric dilemma; we are able to show that the chain sense reversal progresses spontaneously using the stretched intermediate model. A nonlinear DNA model is used to describe the origins and mechanism of base rotation in the stretched intermediate state of DNA. We also propose an experiment that can verify the existence of a stretched intermediate state during B-Z transition, thus opening up fresh grounds for experimentation in this field.

Atomistic view of base flipping in DNA

Base flipping is essential for the enzyme-catalysed methylation of DNA. In our previous studies, the flipping of bases out of duplex DNA in DNA alone and when bound to the (cytosine-C5)-methyltransferase from HhaI (M.HhaI) were investigated via potential of mean force calculations. Insights into various experimental observations were obtained. In the present paper we present an overview of previous computational studies of base flipping along with new detailed structural and energetic analysis on atomic events that contribute to the free energy surfaces. The contributions from different intrinsic and environmental effects to the base-flipping process are explored, and experimental data derived from a variety of methods are reconciled. A detailed protein-facilitated base-flipping mechanism is proposed. Ground-state destabilization is achieved via disruption of the target base Watson-Crick interactions by substitution with favourable DNA-protein interactions. In addition, specific DNA-protein interactions and favourable solvation effects further promote target base flipping along the major groove through the protein matrix, and maximal interactions occur between the DNA and the protein upon reaching the fully flipped state. Other DNA binding proteins that involve distortion of DNA's conformation may use a similar mechanism to that by which M.HhaI facilitates base flipping.

A Hoogsteen base pair embedded in undistorted B-DNA

Hoogsteen base pairs within duplex DNA typically are only observed in regions containing significant distortion or near sites of drug intercalation. We report here the observation of a Hoogsteen base pair embedded within undistorted, unmodified B-DNA. The Hoogsteen base pair, consisting of a syn adenine base paired with an anti thymine base, is found in the 2.1 A resolution structure of the MATalpha2 homeodomain bound to DNA in a region where a specifically and a non-specifically bound homeodomain contact overlapping sites. NMR studies of the free DNA show no evidence of Hoogsteen base pair formation, suggesting that protein binding favors the transition from a Watson-Crick to a Hoogsteen base pair. Molecular dynamics simulations of the homeodomain-DNA complex support a role for the non-specifically bound protein in favoring Hoogsteen base pair formation. The presence of a Hoogsteen base pair in the crystal structure of a protein-DNA complex raises the possibility that Hoogsteen base pairs could occur within duplex DNA and play a hitherto unrecognized role in transcription, replication and other cellular processes.

Free energy and structural pathways of base flipping in a DNA GCGC containing sequence

Structural distortions of DNA are essential for its biological function due to the genetic information of DNA not being physically accessible in the duplex state. Base flipping is one of the simplest structural distortions of DNA and may represent an initial event in strand separation required to access the genetic code. Flipping is also utilized by DNA-modifying and repair enzymes to access specific bases. It is typically thought that base flipping (or base-pair opening) occurs via the major groove whereas minor groove flipping is only possible when mediated by DNA-binding proteins. Here, umbrella sampling with a novel center-of-mass pseudodihedral reaction coordinate was used to calculate the individual potentials of mean force (PMF) for flipping of the Watson-Crick (WC) paired C and G bases in the CCATGCGCTGAC DNA dodecamer. The novel reaction coordinate allowed explicit investigation of the complete flipping process via both the minor and major groove pathways. The minor and major groove barriers to flipping are similar for C base flipping while the major groove barrier is slightly lower for G base flipping. Minor groove flipping requires distortion of the WC partner while the flipping base pulls away from its partner during major groove flipping. The flipped states are represented by relatively flat free energy surfaces, with a small, local minimum observed for the flipped G base. Conserved patterns of phosphodiester backbone dihedral distortions during flipping indicate their essential role in the flipping process. During flipping, the target base tracks along the respective grooves, leading to hydrogen-bonding interactions with neighboring base-pairs. Such hydrogen-bonding interactions with the neighboring sequence suggest a novel mechanism of sequence dependence in DNA dynamics.