Time-resolved fluorescence studies of nucleotide flipping by restriction enzymes

Restriction enzymes Ecl18kI, PspGI and EcoRII-C, specific for interrupted 5-bp target sequences, flip the central base pair of these sequences into their protein pockets to facilitate sequence recognition and adjust the DNA cleavage pattern. We have used time-resolved fluorescence spectroscopy of 2-aminopurine-labelled DNA in complex with each of these enzymes in solution to explore the nucleotide flipping mechanism and to obtain a detailed picture of the molecular environment of the extrahelical bases. We also report the first study of the 7-bp cutter, PfoI, whose recognition sequence (T/CCNGGA) overlaps with that of the Ecl18kI-type enzymes, and for which the crystal structure is unknown. The time-resolved fluorescence experiments reveal that PfoI also uses base flipping as part of its DNA recognition mechanism and that the extrahelical bases are captured by PfoI in binding pockets whose structures are quite different to those of the structurally characterized enzymes Ecl18kI, PspGI and EcoRII-C. The fluorescence decay parameters of all the enzyme-DNA complexes are interpreted to provide insight into the mechanisms used by these four restriction enzymes to flip and recognize bases and the relationship between nucleotide flipping and DNA cleavage.     

How PspGI, catalytic domain of EcoRII and Ecl18kI acquire specificities for different DNA targets

Restriction endonucleases Ecl18kI and PspGI/catalytic domain of EcoRII recognize CCNGG and CCWGG sequences (W stands for A or T), respectively. The enzymes are structurally similar, interact identically with the palindromic CC:GG parts of their recognition sequences and flip the nucleotides at their centers. Specificity for the central nucleotides could be influenced by the strength/stability of the base pair to be disrupted and/or by direct interactions of the enzymes with the flipped bases. Here, we address the importance of these contributions. We demonstrate that wt Ecl18kI cleaves oligoduplexes containing canonical, mismatched and abasic sites in the central position of its target sequence CCNGG with equal efficiencies. In contrast, substitutions in the binding pocket for the extrahelical base alter the Ecl18kI preference for the target site: the W61Y mutant prefers only certain mismatched substrates, and the W61A variant cuts exclusively at abasic sites, suggesting that pocket interactions play a major role in base discrimination. PspGI and catalytic domain of EcoRII probe the stability of the central base pair and the identity of the flipped bases in the pockets. This ‘double check’ mechanism explains their extraordinary specificity for an A/T pair in the flipping position.     

DNA base flipping by both members of the PspGI restriction-modification system

The PspGI restriction–modification system recognizes the sequence CCWGG. R.PspGI cuts DNA before the first C in the cognate sequence and M.PspGI is thought to methylate N4 of one of the cytosines in the sequence. M.PspGI enhances fluorescence of 2-aminopurine in DNA if it replaces the second C in the sequence, while R.PspGI enhances fluorescence when the fluorophore replaces adenine in the central base pair. This strongly suggests that the methyltransferase flips the second C in the recognition sequence, while the endonuclease flips both bases in the central base pair out of the duplex. M.PspGI is the first N4-cytosine MTase for which biochemical evidence for base flipping has been presented. It is also the first type IIP methyltransferase whose catalytic activity is strongly stimulated by divalent metal ions. However, divalent metal ions are not required for its base-flipping activity. In contrast, these ions are required for both base flipping and catalysis by the endonuclease. The two enzymes have similar temperature profiles for base flipping and optimal flipping occurs at temperatures substantially below the growth temperature of the source organism for PspGI and for the catalytic activity of endonuclease. We discuss the implications of these results for DNA binding by these enzymes and their evolutionary origin.     

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.     

The Effect of a G:T Mispair on the Dynamics of DNA

Distortions in the DNA sequence such as damages or mispairs are specifically recognized and processed by DNA repair enzymes. A particular challenge for the enzymatic specificity is the recognition of a wrongly-placed native nucleotide such as thymine in T:G mispairs. An important step of substrate binding which is observed in many repair proteins is the flipping of the target base out of the DNA helix into the enzyme’s active site. In this work we investigate how much the intrinsic dynamics of mispaired DNA is changed compared to canonical DNA. Our molecular dynamics simulations of DNA with and without T:G mispairs show significant differences in the conformation of paired and mispaired DNA. The wobble pair T:G shows local distortions such as twist, shear and stretch which deviate from canonical B form values. Moreover, the T:G mispair is found to be kinetically less stable, exhibiting two states with respect to base opening: a closed state comparable to the canonical base pairs, and a more open state, indicating a proneness for base flip. In addition, we observe that the thymine base in a T:G mispair is significantly more probable to be flipped than thymine in a T:A pair or cytosine in a C:G pair. Such local deformations and in particular the existence of a second, more-open state can be speculated to help the target-site recognition by repair enzymes.     

Nucleotide flipping by restriction enzymes analyzed by 2-aminopurine steady-state fluorescence

Many DNA modification and repair enzymes require access to DNA bases and therefore flip nucleotides. Restriction endonucleases (REases) hydrolyze the phosphodiester backbone within or in the vicinity of the target recognition site and do not require base extrusion for the sequence readout and catalysis. Therefore, the observation of extrahelical nucleotides in a co-crystal of REase Ecl18kI with the cognate sequence, CCNGG, was unexpected. It turned out that Ecl18kI reads directly only the CCGG sequence and skips the unspecified N nucleotides, flipping them out from the helix. Sequence and structure conservation predict nucleotide flipping also for the complexes of PspGI and EcoRII with their target DNAs (/CCWGG), but data in solution are limited and indirect. Here, we demonstrate that Ecl18kI, the C-terminal domain of EcoRII (EcoRII-C) and PspGI enhance the fluorescence of 2-aminopurines (2-AP) placed at the centers of their recognition sequences. The fluorescence increase is largest for PspGI, intermediate for EcoRII-C and smallest for Ecl18kI, probably reflecting the differences in the hydrophobicity of the binding pockets within the protein. Omitting divalent metal cations and mutation of the binding pocket tryptophan to alanine strongly increase the 2-AP signal in the Ecl18kI–DNA complex. Together, our data provide the first direct evidence that Ecl18kI, EcoRII-C and PspGI flip nucleotides in solution.     

Uncoupling of nucleotide flipping and DNA bending by the t4 pyrimidine dimer DNA glycosylase

Bacteriophage T4 pyrimidine dimer glycosylase (T4-Pdg) is a base excision repair protein that incises DNA at cyclobutane pyrimidine dimers that are formed as a consequence of exposure to ultraviolet light. Cocrystallization of T4-Pdg with substrate DNA has shown that the adenosine opposite the 5'-thymine of a thymine-thymine (TT) dimer is flipped into an extrahelical conformation and that the DNA backbone is kinked 60 degrees in the enzyme-substrate (ES) complex. To examine the kinetic details of the precatalytic events in the T4-Pdg reaction mechanism, investigations were designed to separately assess nucleotide flipping and DNA bending. The fluorescent adenine base analogue, 2-aminopurine (2-AP), placed opposite an abasic site analogue, tetrahydrofuran, exhibited a 2.8-fold increase in emission intensity when flipped in the ES complex. Using the 2-AP fluorescence signal for nucleotide flipping, kon and koff pre-steady-state kinetic measurements were determined. DNA bending was assessed by fluorescence resonance energy transfer using fluorescent donor-acceptor pairs located at the 5'-ends of oligonucleotides in duplex DNA. The fluorescence intensity of the donor fluorophore was quenched by 15% in the ES complex as a result of an increased efficiency of energy transfer between the labeled ends of the DNA in the bent conformation. Kinetic analyses of the bending signal revealed an off rate that was 2.5-fold faster than the off rate for nucleotide flipping. These results demonstrate that the nucleotide flipping step can be uncoupled from the bending of DNA in the formation of an ES complex.     

Structural studies of DNA fragments: the G.T wobble base pair in A, B and Z DNA; the G.A base pair in B-DNA

The crystal structures of five double helical DNA fragments containing non-Watson-Crick complementary base pairs are reviewed. They comprise four fragments containing G.T base pairs: two deoxyoctamers d(GGGGCTCC) and d(GGGGTCCC) which crystallise as A type helices; a deoxydodecamer d(CGCGAATTTGCG) which crystallises in the B-DNA conformation; and the deoxyhexamer d(TGCGCG), which crystallises as a Z-DNA helix. In all four duplexes the G and T bases form wobble base pairs, with bases in the major tautomer forms and hydrogen bonds linking N1 of G with O2 of T and O6 of G with N3 of T. The X-ray analyses establish that the G.T wobble base pair can be accommodated in the A, B or Z double helix with minimal distortion of the global conformation. There are, however, changes in base stacking in the neighbourhood of the mismatched bases. The fifth structure, d(CGCGAATTAGCG), contains the purine purine mismatch G.A where G is in the anti and A in the syn conformation. The results represent the first direct structure determinations of base pair mismatches in DNA fragments and are discussed in relation to the fidelity of replication and mismatch recognition.     

Binding of EcoP15I DNA methyltransferase to DNA reveals a large structural distortion within the recognition sequence

EcoP15I DNA methyltransferase, a member of the type III restriction-modification system, binds to the sequence 5'-CAGCAG-3' transferring a methyl group from S-adenosyl-l-methionine to the second adenine base. We have investigated protein-DNA interactions in the methylase-DNA complex by three methods. Determination of equilibrium dissociation constants indicated that the enzyme had higher affinity for DNA containing mismatches at the target base within the recognition sequence. Potassium permanganate footprinting studies revealed that there was a hyper-reactive permanganate cleavage site coincident with adenine that is the target base for methylation. More importantly, to detect DNA conformational alterations within the enzyme-DNA complexes, we have used a fluorescence-based assay. When EcoP15I DNA methyltransferase bound to DNA containing 2-aminopurine substitutions within the cognate sequence, an eight to tenfold fluorescent enhancement resulting from enzymatic flipping of the target adenine base was observed. Furthermore, fluorescence spectroscopy analysis showed that the changes attributable to structural distortion were specific for only the bases within the recognition sequence. More importantly, we observed that both the adenine bases in the recognition site appear to be structurally distorted to the same extent. While the target adenine base is probably flipped out of the DNA duplex, our results also suggest that fluorescent enhancements could be derived from protein-DNA interactions other than base flipping. Taken together, our results support the proposed base flipping mechanism for adenine methyltransferases.     

Stopped-flow and mutational analysis of base flipping by the Escherichia coli Dam DNA-(adenine-N6)-methyltransferase

By stopped-flow kinetics using 2-aminopurine as a probe to detect base flipping, we show here that base flipping by the Escherichia coli Dam DNA-(adenine-N6)-methyltransferase (MTase) is a biphasic process: target base flipping is very fast (k(flip)>240 s(-1)), but binding of the flipped base into the active site pocket of the enzyme is slow (k=0.1-2 s(-1)). Whereas base flipping occurs in the absence of S-adenosyl-l-methionine (AdoMet), binding of the target base in the active site pocket requires AdoMet. Our data suggest that the tyrosine residue in the DPPY motif conserved in the active site of DNA-(adenine-N6)-MTases stacks to the flipped target base. Substitution of the aspartic acid residue of the DPPY motif by alanine abolished base flipping, suggesting that this residue contacts and stabilizes the flipped base. The exchange of Ser188 located in a loop next to the active center by alanine led to a seven- to eightfold reduction of k(flip), which was also reduced with substrates having altered GATC recognition sites and in the absence of AdoMet. These findings provide evidence that the enzyme actively initiates base flipping by stabilizing the transition state of the process. Reduced rates of base flipping in substrates containing the target base in a non-canonical sequence demonstrate that DNA recognition by the MTase starts before base flipping. DNA recognition, cofactor binding and base flipping are correlated and efficient base flipping takes place only if the enzyme has bound to a cognate target site and AdoMet is available.     

Effects of bulge composition and flanking sequence on the kinking of DNA by bulged bases

We recently showed that bulged bases kink duplex DNA, with the degree of kinking increasing in roughly equal increments as the number of bases in the bulge increases from one to four [Hsieh, C.-H., & Griffith, J.D. (1989) Proc. Natl. Acad. Sci. U.S.A. 86, 4833-4837]. Here we have examined the kinking of DNA by single A, C, G, or T bulges with different neighboring base pairs. Synthetic 30 base pair (bp) duplex DNAs containing 2 single-base bulges spaced by 10 bp were ligated head to tail, and their electrophoretic behavior in highly cross-linked gels was examined. All bulge-containing DNAs showed marked electrophoretic retardations as compared to non-bulge-containing DNA. Regardless of the sequence of the flanking base pairs, purine bulges produced greater retardations than pyrimidine bulges. Furthermore, C and T bulges produced the same retardations as did G and A bulges. Bulged DNA containing different flanking base pairs showed marked differences in electrophoretic mobility. For C-bulged DNA, the greatest retardations were observed with G.C neighbors, the least with T.A neighbors, and an intermediate amount with a mixture of neighboring base pairs. For A-bulged DNA, the retardations were greatest with G.C neighbors, less with T.A neighbors, even less with a mixture of neighboring base pairs, and finally least with C.G neighbors. Thus flanking base pairs affect C-bulged DNA and A-bulged DNA differently, and G.C and C.G flanking base pairs were seen to have very different effects. These results imply an important role of base stacking in determining how neighboring base pairs influence the kinking of DNA by a single-base bulge.     

Energy barriers and rates of tautomeric transitions in DNA bases: ab initio quantum chemical study

Tautomeric transitions of DNA bases are proton transfer reactions, which are important in biology. These reactions are involved in spontaneous point mutations of the genetic material. In the present study, intrinsic reaction coordinates (IRC) analyses through ab initio quantum chemical calculations have been carried out for the individual DNA bases A, T, G, C and also A:T and G:C base pairs to estimate the kinetic and thermodynamic barriers using MP2/6-31G** method for tautomeric transitions. Relatively higher values of kinetic barriers (about 50-60 kcal/mol) have been observed for the single bases, indicating that tautomeric alterations of isolated single bases are quite unlikely. On the other hand, relatively lower values of the kinetic barriers (about 20-25 kcal/mol) for the DNA base pairs A:T and G:C clearly suggest that the tautomeric shifts are much more favorable in DNA base pairs than in isolated single bases. The unusual base pairing A':C, T':G, C':A or G':T in the daughter DNA molecule, resulting from a parent DNA molecule with tautomeric shifts, is found to be stable enough to result in a mutation. The transition rate constants for the single DNA bases in addition to the base pairs are also calculated by computing the free energy differences between the transition states and the reactants.     

The structure of guanosine-thymidine mismatches in B-DNA at 2.5-A resolution

The structure of the deoxyoligomer d(C-G-C-G-A-A-T-T-T-G-C-G) was determined at 2.5-A resolution by single crystal x-ray diffraction techniques. The final R factor is 18% with the location of 71 water molecules. The oligomer crystallizes in a B-DNA-type conformation, with two strands interacting to form a dodecamer duplex. The double helix consists of four A X T and six G X C Watson-Crick base pairs and two G X T mismatches. The G X T pairs adopt a "wobble" structure with the thymine projecting into the major groove and the guanine into the minor groove. The mispairs are accommodated in the normal double helix by small adjustments in the conformation of the sugar phosphate backbone. A comparison with the isomorphous parent compound containing only Watson-Crick base pairs shows that any changes in the structure induced by the presence of G X T mispairs are highly localized. The global conformation of the duplex is conserved. The G X T mismatch has already been studied by x-ray techniques in A and Z helices where similar results were found. The geometry of the mispair is essentially identical in all structures so far examined, irrespective of the DNA conformation. The hydration is also similar with solvent molecules bridging the functional groups of the bases via hydrogen bonds. Hydration may be an important factor in stabilizing G X T mismatches. A characteristic of Watson-Crick paired A X T and G X C bases is the pseudo 2-fold symmetry axis in the plane of the base pairs. The G X T wobble base pair is pronouncedly asymmetric. This asymmetry, coupled with the disposition of functional groups in the major and minor grooves, provides a number of features which may contribute to the recognition of the mismatch by repair enzymes.     

DNA oligonucleotides with A, T, G or C opposite an abasic site: structure and dynamics

Abasic sites are common DNA lesions resulting from spontaneous depurination and excision of damaged nucleobases by DNA repair enzymes. However, the influence of the local sequence context on the structure of the abasic site and ultimately, its recognition and repair, remains elusive. In the present study, duplex DNAs with three different bases (G, C or T) opposite an abasic site have been synthesized in the same sequence context (5′-CCA AAG₆ XA₈C CGG G-3′, where X denotes the abasic site) and characterized by 2D NMR spectroscopy. Studies on a duplex DNA with an A opposite the abasic site in the same sequence has recently been reported [Chen,J., Dupradeau,F.-Y., Case,D.A., Turner,C.J. and Stubbe,J. (2007) Nuclear magnetic resonance structural studies and molecular modeling of duplex DNA containing normal and 4′-oxidized abasic sites. Biochemistry, 46, 3096–3107]. Molecular modeling based on NMR-derived distance and dihedral angle restraints and molecular dynamics calculations have been applied to determine structural models and conformational flexibility of each duplex. The results indicate that all four duplexes adopt an overall B-form conformation with each unpaired base stacked between adjacent bases intrahelically. The conformation around the abasic site is more perturbed when the base opposite to the lesion is a pyrimidine (C or T) than a purine (G or A). In both the former cases, the neighboring base pairs (G6-C21 and A8-T19) are closer to each other than those in B-form DNA. Molecular dynamics simulations reveal that transient H-bond interactions between the unpaired pyrimidine (C20 or T20) and the base 3′ to the abasic site play an important role in perturbing the local conformation. These results provide structural insight into the dynamics of abasic sites that are intrinsically modulated by the bases opposite the abasic site.     

Structural Studies of DNA Fragments: The G·T Wobble Base Pair in A,B and Z DNA; The G·A Base Pair in B-DNA

Abstract

The crystal structures of five double helical DNA fragments containing non-Watson-Crick complementary base pairs are reviewed. They comprise four fragments containing G·T base pairs: two deoxyoctamers d(GGGGCTCC) and d(GGGGTCCC) which crystallise as A type helices; a deoxydodecamer d(CGCGAATTTGCG) which crystallises in the B-DNA conformation; and the deoxyhexamer d(TGCGCG), which crystallises as a Z-DNA helix. In all four duplexes the G and T bases form wobble base pairs, with bases in the major tautomer forms and hydrogen bonds linking N1 of G with 02 of T and 06 of G with N3 of T. The X-ray analyses establish that the G·T wobble base pair can be accommodated in the A, B or Z double helix with minimal distortion of the global conformation. There are, however, changes in base stacking in the neighbourhood of the mismatched bases. The fifth structure, d(CGCGAATTAGCG), contains the purine purine mismatch G·A where G is in the anti and A in the syn conformation. The results represent the first direct structure determinations of base pair mismatches in DNA fragments and are discussed in relation to the fidelity of replication and mismatch recognition.     

The P1 phage replication protein RepA contacts an otherwise inaccessible thymine N3 proton by DNA distortion or base flipping

The RepA protein from bacteriophage P1 binds DNA to initiate replication. RepA covers one face of the DNA and the binding site has a completely conserved T that directly faces RepA from the minor groove at position +7. Although all four bases can be distinguished through contacts in the major groove of B-form DNA, contacts in the minor groove cannot easily distinguish between A and T bases. Therefore the 100% conservation at this position cannot be accounted for by direct contacts approaching into the minor groove of B-form DNA. RepA binding sites with modified base pairs at position +7 were used to investigate contacts with RepA. The data show that RepA contacts the N3 proton of T at position +7 and that the T=A hydrogen bonds are already broken in the DNA before RepA binds. To accommodate the N3 proton contact the T₊₇ /A_+7^′ base pair must be distorted. One possibility is that T₊₇ is flipped out of the helix. The energetics of the contact allows RepA to distinguish between all four bases, accounting for the observed high sequence conservation. After protein binding, base pair distortion or base flipping could initiate DNA melting as the second step in DNA replication.     

Spectroscopic and theoretical insights into sequence effects of aminofluorene-induced conformational heterogeneity and nucleotide excision repair

A systematic spectroscopic and computational study was conducted in order to probe the influence of base sequences on stacked (S) versus B-type (B) conformational heterogeneity induced by the major dG adduct derived from the model carcinogen 7-fluoro-2-aminofluorene (FAF). We prepared and characterized eight 12-mer DNA duplexes (-AG*N- series, d[CTTCTAG*NCCTC]; -CG*N- series, d[CTTCTCG*NCCTC]), in which the central guanines (G*) were site-specifically modified with FAF with varying flanking bases (N = G, A, C, T). S/B heterogeneity was examined by CD, UV, and dynamic 19F NMR spectroscopy. All the modified duplexes studied followed a typical dynamic exchange between the S and B conformers in a sequence dependent manner. Specifically, purine bases at the 3'-flanking site promoted the S conformation (G > A > C > T). Simulation analysis showed that the S/B energy barriers were in the 14-16 kcal/mol range. The correlation times (tau = 1/kappa) were found to be in the millisecond range at 20 degrees C. The van der Waals energy force field calculations indicated the importance of the stacking interaction between the carcinogen and neighboring base pairs. Quantum mechanics calculations showed the existence of correlations between the total interaction energies (including electrostatic and solvation effects) and the S/B population ratios. The S/B equilibrium seems to modulate the efficiency of Escherichia coli UvrABC-based nucleotide excision repair in a conformation-specific manner: i.e., greater repair susceptibility for the S over B conformation and for the -AG*N- over the -CG*N- series. The results indicate a novel structure-function relationship, which provides insights into how bulky DNA adducts are accommodated by UvrABC proteins.     

HhaI DNA methyltransferase uses the protruding Gln237 for active flipping of its target cytosine

Access to a nucleotide by its rotation out of the DNA helix (base flipping) is used by numerous DNA modification and repair enzymes. Despite extensive studies of the paradigm HhaI methyltransferase, initial events leading to base flipping remained elusive. Here we demonstrate that the replacement of the target C:G pair with the 2-aminopurine:T pair in the DNA or shortening of the side chain of Gln237 in the protein severely perturb base flipping, but retain specific DNA binding. Kinetic analyses and molecular modeling suggest that a steric interaction between the protruding side chain of Gln237 and the target cytosine in B-DNA reduces the energy barrier for flipping by 3 kcal/mol. Subsequent stabilization of an open state by further 4 kcal/mol is achieved through specific hydrogen bonding of the side chain to the orphan guanine. Gln237 thus plays a key role in actively opening the target C:G pair by a "push-and-bind" mechanism.     

M.HhaI binds tightly to substrates containing mismatches at the target base.

The (cytosine-5) DNA methyltransferase M.HhaI causes its target cytosine base to be flipped completely out of the DNA helix upon binding. We have investigated the effects of replacing the target cytosine by other, mismatched bases, including adenine, guanine, thymine and uracil. We find that M.HhaI binds more tightly to such mismatched substrates and can even transfer a methyl group to uracil if a G:U mismatch is present. Other mismatched substrates in which the orphan guanine is changed exhibit similar behavior. Overall, the affinity of DNA binding correlates inversely with the stability of the target base pair, while the nature of the target base appears irrelevant for complex formation. The presence of a cofactor analog. S-adenosyl-L-homocysteine, greatly enhances the selectivity of the methyltransferase for cytosine at the target site. We propose that the DNA methyltransferases have evolved from mismatch binding proteins and that base flipping was, and still is, a key element in many DNA-enzyme interactions.