Non-additivity of sequence-specific enzyme-DNA interactions in the EcoRI DNA methyltransferase

Isospecific adenine DNA methyltransferases show distinct preferences towards DNA substrates

Here, we report results on systematic analysis of DNA substrate preferences of three N6-adenine β-class DNA methyltransferases that are part of the type II restriction-modification systems. The studied enzymes were: M.EcoVIII, M.HindIII and M.LlaCI, which although found in phylogenetically distant bacteria (γ-proteobacteria and low-GC Gram-positive bacteria), recognize the same palindromic specific sequence 5′-AAGCTT-3′ and catalyze formation of N6-methyladenine at the first A-residue. As expected overall the enzymes share the most analyzed features, but they show also some distinct differences in substrate recognition. Therefore DNA methylation reactions were carried out not only under standard, but also under relaxed conditions using DMSO or glycerol. We found that all of these enzymes preferred DNA containing a hemimethylated target site, but differ in modification of ssDNA, especially more pronounced for M.EcoVIII under relaxed conditions. In these conditions they also have shown varied preferences toward secondary sites, which differ by one nucleotide from specific sequence. They preferred sequences with substitutions at the 1^st (A¹ → G/C) and at the 2^nd position (A² → C), while sites with substitutions at the 3^rd position (G³ → A/C) were modified less efficiently. Kinetic parameters of the methylation reaction carried out by M.EcoVIII were determined. Methylation efficiency (k_cat/K_m) of secondary sites was 4.5–10 times lower when compared to the unmethylated specific sequences, whilst efficiency observed for the hemimethylated substrate was almost 4.5 times greater. We also observed a distinct effect of analyzed enzymes on unspecific interaction with DNA phosphate backbone. We concluded that for all three enzymes the most critical is the phosphodiester bond between G³-C⁴ nucleotides at the center of the target site.

Pattern matching through Chaos Game Representation: bridging numerical and discrete data structures for biological sequence analysis

Background

Chaos Game Representation (CGR) is an iterated function that bijectively maps discrete sequences into a continuous domain. As a result, discrete sequences can be object of statistical and topological analyses otherwise reserved to numerical systems. Characteristically, CGR coordinates of substrings sharing an L-long suffix will be located within 2^-L distance of each other. In the two decades since its original proposal, CGR has been generalized beyond its original focus on genomic sequences and has been successfully applied to a wide range of problems in bioinformatics. This report explores the possibility that it can be further extended to approach algorithms that rely on discrete, graph-based representations.

Results

The exploratory analysis described here consisted of selecting foundational string problems and refactoring them using CGR-based algorithms. We found that CGR can take the role of suffix trees and emulate sophisticated string algorithms, efficiently solving exact and approximate string matching problems such as finding all palindromes and tandem repeats, and matching with mismatches. The common feature of these problems is that they use longest common extension (LCE) queries as subtasks of their procedures, which we show to have a constant time solution with CGR. Additionally, we show that CGR can be used as a rolling hash function within the Rabin-Karp algorithm.

Conclusions

The analysis of biological sequences relies on algorithmic foundations facing mounting challenges, both logistic (performance) and analytical (lack of unifying mathematical framework). CGR is found to provide the latter and to promise the former: graph-based data structures for sequence analysis operations are entailed by numerical-based data structures produced by CGR maps, providing a unifying analytical framework for a diversity of pattern matching problems.

Structure, function and mechanism of exocyclic DNA methyltransferases

DNA MTases (methyltransferases) catalyse the transfer of methyl groups to DNA from AdoMet (S-adenosyl-L-methionine) producing AdoHcy (S-adenosyl-L-homocysteine) and methylated DNA. The C5 and N4 positions of cytosine and N6 position of adenine are the target sites for methylation. All three methylation patterns are found in prokaryotes, whereas cytosine at the C5 position is the only methylation reaction that is known to occur in eukaryotes. In general, MTases are two-domain proteins comprising one large and one small domain with the DNA-binding cleft located at the domain interface. The striking feature of all the structurally characterized DNA MTases is that they share a common core structure referred to as an 'AdoMet-dependent MTase fold'. DNA methylation has been reported to be essential for bacterial virulence, and it has been suggested that DNA adenine MTases (Dams) could be potential targets for both vaccines and antimicrobials. Drugs that block Dam could slow down bacterial growth and therefore drug-design initiatives could result in a whole new generation of antibiotics. The transfer of larger chemical entities in a MTase-catalysed reaction has been reported and this represents an interesting challenge for bio-organic chemists. In general, amino MTases could therefore be used as delivery systems for fluorescent or other reporter groups on to DNA. This is one of the potential applications of DNA MTases towards developing non-radioactive DNA probes and these could have interesting applications in molecular biology. Being nucleotide-sequence-specific, DNA MTases provide excellent model systems for studies on protein-DNA interactions. The focus of this review is on the chemistry, enzymology and structural aspects of exocyclic amino MTases.

Bacteriophage T4Dam (DNA-(adenine-N6)-methyltransferase): evidence for two distinct stages of methylation under single turnover conditions

We compared the (pre)steady-state and single turnover methylation kinetics of bacteriophage T4Dam (DNA-(adenine-N6)-methyltransferase)-mediated methyl group transfer from S-adenosyl-l-methionine (AdoMet) to oligodeoxynucleotide duplexes containing a single recognition site (palindrome 5'-GATC/5'-GATC) or some modified variant. T4Dam-AdoMet functions as a monomer under steady-state conditions (enzyme/DNA << 1), whereas under single turnover conditions (enzyme/DNA > 1), a catalytically active complex containing two Dam-AdoMet molecules is formed initially, and two methyl groups are transferred per duplex (to produce a methylated duplex and S-adenosyl-l-homocysteine (AdoHcy)). We propose that the single turnover reaction proceeds in two stages. First, two preformed T4Dam-AdoMet complexes bind opposite strands of the unmodified target site, and one enzyme molecule catalyzes the rapid transfer of the AdoMet-methyl group (kmeth1 = 0.21 s-1); this is 2.5-fold slower than the rate observed with monomeric T4Dam-AdoMet bound under pre-steady-state conditions for burst determination. In the second stage, methyl transfer to adenine in GATC on the complementary strand occurs at a rate that is 1 order of magnitude slower (kmeth2 = 0.023 s-1). We suggest that under single turnover conditions, methylation of the second strand is rate-limited by Dam-AdoHcy dissociation or its clearance from the methylated complementary strand. The hemimethylated duplex 5'-GATC/5'-GMTC also interacts with T4Dam-AdoMet complexes in two stages under single turnover reaction conditions. The first stage (kmeth1) reflects methylation by dimeric T4Dam-AdoMet productively oriented to the strand with the adenine residue capable of methylation. The slower second stage (kmeth2) reflects methylation by enzyme molecules non-productively oriented to the GMTC chain, which then have to re-orient to the opposite productive chain. Substitutions of bases and deletions in the recognition site affect the kinetic parameters in different fashions. When the GAT portion of GATC was disrupted, the proportion of the initial productive enzyme-substrate complexes was sharply reduced.

Preferential methylation of unmethylated DNA by Mammalian de novo DNA methyltransferase Dnmt3a

DNA methylation is an epigenetic modification of DNA. There are currently three catalytically active mammalian DNA methyltransferases, DNMT1, -3a, and -3b. DNMT1 has been shown to have a preference for hemimethylated DNA and has therefore been termed the maintenance methyltransferase. Although previous studies on DNMT3a and -3b revealed that they act as functional enzymes during development, there is little biochemical evidence about how new methylation patterns are established and maintained. To study this mechanism we have cloned and expressed Dnmt3a using a baculovirus expression system. The substrate specificity of Dnmt3a and molecular mechanism of its methylation reaction were then analyzed using a novel and highly reproducible assay. We report here that Dnmt3a is a true de novo methyltransferase that prefers unmethylated DNA substrates more than 3-fold to hemimethylated DNA. Furthermore, Dnmt3a binds DNA nonspecifically, regardless of the presence of CpG dinucleotides in the DNA substrate. Kinetic analysis supports an Ordered Bi Bi mechanism for Dnmt3a, where DNA binds first, followed by S-adenosyl-l-methionine.

Thermodynamic and kinetic analyses for understanding sequence-specific DNA recognition

Thermodynamic and kinetic analyses of biomolecular interactions reveal details of the energetic and dynamic features of molecular recognition processes, and complement structural analyses of the free and complexed conformations. The recent improvements in both isothermal titration calorimetry and surface plasmon resonance sensoring provide powerful tools for analysing biomolecular interactions in thermodynamic and kinetic approaches. The thermodynamic and kinetic parameters obtained for binding between protein and DNA indicate the mechanism of specific DNA recognition, in the high-resolution structures of the protein-DNA complexes. The effects of temperature and ionic strength reflect the conformational changes of the protein and DNA molecules upon complex formation, including important contributions of water and solutes. When combined with mutational studies, the interactions can be reduced to several energetic contributions from individual contacts. These studies should be useful to determine general features of protein functions in genetic regulation.

Interaction of the phage T4 Dam DNA-[N6-adenine] methyltransferase with oligonucleotides containing native or modified (defective) recognition sites

The DNA-[N 6-adenine]-methyltransferase (Dam MTase) of phage T4 catalyzes methyl group transfer from S-adenosyl-l-methionine (AdoMet) to the N6-position of adenine in the palindromic sequence, GATC. We have used a gel shift assay to monitor complex formation between T4 Dam and various synthetic duplex oligonucleotides, either native or modified/defective. The results are summarized as follows. (i) T4 Dam bound with approximately 100-fold higher affinity to a 20mer specific (GATC-containing) duplex containing the canonical palindromic methylation sequence, GATC, than to a non-specific duplex containing another palindrome, GTAC. (ii) Compared with the unmethylated duplex, the hemimethylated 20mer specific duplex had a slightly increased ( approximately 2-fold) ability to form complexes with T4 Dam. (iii) No stable complex was formed with a synthetic 12mer specific (GATC-containing) duplex, although T4 Dam can methylate it. This indicates that there is no relation between formation of a catalytically competent 12mer-Dam complex and one stable to gel electrophoresis. (iv) Formation of a stable complex did not require that both strands be contiguous or completely complementary. Absence of a single internucleotide phosphate strongly reduced complex formation only when missing between the T and C residues. This suggests that if T4 Dam makes critical contact(s) with a backbone phosphate(s), then the one between T and C is the only likely candidate. Having only one half of the recognition site intact on one strand was sufficient for stable complex formation provided that the 5'G.C base-pairs be present at both ends of the palindromic, GATC. Since absence of either a G or C abolished T4 Dam binding, we conclude that both strands are recognized by T4 Dam.

Dam methyltransferase from Escherichia coli: kinetic studies using modified DNA oligomers: nonmethylated substrates

Steady-state kinetics of the N6-adenine Dam methyltransferase have been measured using as substrates non-self-complementary tetradecanucleotide duplexes that contain the GATC target sequence. Modifications in the GATC target sequence of one or both of the strands included substitution of guanine by hypoxanthine, thymine by uracil or 5-ethyl-uracil and adenine by diamino-purine (2-amino-adenine). Thermodynamic parameters for the 14-mer duplexes were also determined. DNA methylation of duplexes containing single dl for dG substitution of the Dam recognition site was little perturbed compared with the canonical substrate. Replacement of dG residues by dl in both strands resulted in a decrease of the specificity constant. Substitution in both strands appears to be cumulative. Substitution of the methyl-accepting adenine residues by 2-amino-adenine resulted in surprisingly little perturbation. Dam methyltransferase is rather tolerant to different substitutions. The results show much less spread than those for the analogous hemimethylated substrates studied previously (Marzabal et al., 1995). The absence of the methylation marker appears to be deleterious to the specificity of the transition state of the active complex, while the binding of the DNA substrate to the enzyme appears to be mostly determined by the thermodynamic stability of the DNA duplex.

Dam methylase from Escherichia coli: kinetic studies using modified DNA oligomers: hemimethylated substrates

We have measured steady-state kinetics of the N6-adenine methyltransferase Dam Mtase using as substrates non-selfcomplementary tetradecamer duplexs (d[GCCGGATCTAGACG]-d[CGTCTAGATCC-GGC]) containing the hemimethylated GATC target sequence in one or the other strand and modifications in the GATC target sequence of the complementary strands. Modifications included substitution of guanine by hypoxanthine (I), thymine by uracil (U) or 5-ethyl-uracil (E) and adenine by 2,6-diamino-purine (D). Thermodynamic parameters were obtained from the concentration dependence of the melting temperature (Tm) of the duplexes. Large differences in DNA methylation of duplexes containing single dI for dG substitution of the Dam recognition site were observed compared with the canonical substrate, if the substitution involved the top strand (on the G.C rich side). Substitution in either strand by uracil (dU) or 5-ethyluracil (dE) resulted in small perturbation of the methylation patterns. When 2,6-diamino-purine (dD) replaced the adenine to be methylated, small, but significant methylation was observed. The kinetic parameters of the methylation reaction were compared with the thermodynamic free energies and significant correlation was observed.

Sequence-dependent effects on DNA stability resulting from guanosine replacements by inosine

The effect of replacing a G.C base-pair with an I.C base-pair on DNA stability was investigated for a related set of 14-mers. DNA melting analysis of the 14-mers was used to determine delta Hzero, delta Szero and delta G(zero)37 of the double to single stranded transition. All 14mers were shown to have B-DNA character by circular dichroism analysis. 14mers substituted with a single inosine in place of guanosine at different positions showed that consequences on DNA stability are sequence-dependent. Large changes in delta Hzero and delta Szero result when inosine is substituted within the trinucleotide sequence d(TCG).d(CGA) while substitution within d(TCC).d(GGA) causes minor changes in enthalpy and entropy. Moreover, some 14-mers with two inosine substitutions five base-pairs apart showed non-additive free energy changes for the double to single stranded transition. These results clearly indicate that the structural consequences of replacing a single guanosine with an inosine are transmitted over a significant distance.

Kinetic characterization of the EcaI methyltransferase

A kinetic analysis of the EcaI adenine-N6-specific methyltransferase (MTase) is presented. The enzyme catalyzes the transfer of a methyl group from S-adenosyl-L-methionine (AdoMet) to the adenine of the GGTNACC sequence with a random rapid-equilibrium mechanism. Experiments with a synthetic, 14-bp DNA substrate suggest that recognition of the specific site of DNA occurs after the binding of AdoMet. Proton concentration does not affect the dissociation constant of AdoMet while Vm and the dissociation constant of DNA show a maximum around pH 8. Increasing the amount of S-adenosyl-L-homocysteine decreases the inhibitory effect of methylated DNA which proves the active role of AdoMet in site recognition. Experiments with hemimethylated DNA show that the methylase binds the double-stranded DNA asymmetrically.

EcoRI DNA methyltransferase-DNA interactions

We present a novel strategy with synthetic hemimethylated DNA substrates containing uracil for thymine and inosine for guanosine replacements and EcoRI DNA methyltransferase to characterize the importance of major groove hydrophobic groups to the sequence-specific modification of DNA. The bacterial Mtase uses S-adenosyl-L-methionine to methylate the double-stranded DNA site 5'GAATTC3' at the N6 position of the central adenosine of each strand. Uracil substitution in either strand at the outer thymine (5'GAATUC3') causes 2.2- and 1.7-fold improvements in specificity (kcat/KmDNA). The fact that the specificity constant for the substrate containing uracil in both strands is identical to the value expected for noninteracting substitutions suggests that no significant methyltransferase-DNA interactions are altered beyond the site of either substitution. Similar analysis of the internal thymine (5'GAAUTC3') also shows these methyl groups to make a negative contribution to specificity, although the observed nonadditivity with the doubly modified substrate clearly shows methyltransferase-DNA interactions beyond the site of substitution to be affected in this case. To further probe the effect of analogue incorporation on methyltransferase-DNA interactions beyond the site of substitution, the relatively "silent" and additive uracil changes (5'GAATUC3') were combined with inosine for guanosine substitutions (e.g., 5'IAATTC3') known to have significant negative effects on specificity. In contrast to the additivity observed with the outer thymines, these studies show significant changes in methyltransferase-DNA interactions caused by the removal of the thymine methyls. Our results implicate a complex and flexible methyltransferase-DNA interface in which subtle structural changes in the substrate are transmitted over the entire canonical site.(ABSTRACT TRUNCATED AT 250 WORDS)