Kinetic analyses of divalent cation-dependent EcoRV digestions on a DNA-immobilized quartz crystal microbalance

Enzymatic digestion with a type IIP restriction endonuclease EcoRV was investigated on a DNA-immobilized 27-MHz quartz crystal microbalance (QCM). Real-time observations of both the enzyme binding process and the DNA cleavage process of EcoRV were followed by frequency (mass) changes on the QCM, which were dependent on divalent cations such as Ca(2+) or Mg(2+). In the presence of Ca(2+), the site-specific binding of EcoRV to DNA could be observed, without the catalytic process. On the other hand, in the presence of Mg(2+), both the binding of the enzyme to the specific DNA (mass increase) and the site-specific cleavage reaction (mass decrease) could be observed continuously from QCM frequency changes. From time courses of frequency (mass) changes, each kinetic parameter, namely binding rate constants (k(on)), dissociation rate constants (k(off)), dissociation constants (K(d)) of EcoRV to DNA, and catalytic rate constant (k(cat)) of the cleavage reaction, could be determined. The binding kinetic parameters of EcoRV in the presence of Ca(2+) were consistent with those of the binding process followed by the cleavage process in the presence of Mg(2+). The k(cat) value obtained by the QCM method was also consistent with that obtained by other methods. This study is the first to simultaneously determine k(on), k(off), and k(cat) for a type IIP restriction endonuclease on one device.     

Osmotic pressure effects on EcoRV cleavage and binding

Investigations of DNA binding proteins frequently measure pH and salt dependence, but relatively few studies measure protein binding in high concentrations of small molecules often found in vivo. We have measured kinetics of the restriction enzyme EcoRV in concentrated solutions of three small cosolvents that produce osmotic pressures from 0.25 to 2.5 mol/kg (6 to 62 atm or water activity of 0.995 to 0.956). We have correlated DNA cleavage and binding parameters with four solution parameters (dielectric constant, viscosity, water concentration, and water activity). We found that the responses of maximum velocity (Vmax) and the dissociation constant for nonspecific binding (Kd,ns) best correlate with water activity. The Michaelis constant (Km) correlates with both water activity and solution viscosity, the latter due to the highly dilute reactant concentrations, which make enzyme-substrate combination diffusion limited. Dielectric constant does not influence any of the kinetic parameters, which is consistent with a view that protein and DNA are preferentially hydrated, and excluded solutes cannot affect the local dielectric constant.     

Investigation of the inhibitory role of phosphorothioate internucleotidic linkages on the catalytic activity of the restriction endonuclease EcoRV

The inhibitory effect of phosphorothioate residues, located within one strand of double-stranded DNA, on the hydrolytic activity of the restriction endonuclease EcoRV was investigated. Specific incorporation of a phosphorothioate group at the site of cleavage yielded the sequence 5'-GATsATC-3'. This modified sequence was cleaved at a relative rate of 0.1 compared to the unmodified substrate. Substrates 5'-GATsAsTC-3' and 5'-GsATsATC-3', both containing one additional phosphorothioate substitution, were linearized at a rate of 0.04 relative to unmodified DNA. However, under the same conditions, fully dAMPS-substituted DNA was found to be virtually resistant to the hydrolytic activity of EcoRV. Further experiments showed that double-stranded DNA fragments generated by PCR containing phosphorothioate groups within both strands are potent inhibitors of EcoRV catalysis. The inhibition was independent of whether the inhibitor fragment contained an EcoRV recognition site. We concluded that substitution of the phosphate group at the site of cleavage by a phosphorothioate residue decreases the rate of EcoRV-catalyzed hydrolysis most significantly. Substitution of other phosphate groups within the recognition sequence plays a limited role in enzyme inhibition. The presence of multiple dNMPS residues at regions of the DNA removed from the EcoRV recognition site may decrease the amount of enzyme available for catalysis by nonspecific binding to EcoRV.     

Crowding effects on EcoRV kinetics and binding

The cytosol of the cell contains high concentrations of small and large macromolecules, but experimental data are often obtained in dilute solutions that do not reflect in vivo conditions. We have studied the crowding effect that large macromolecules have on EcoRV cleavage by adding high-molecular-weight Ficoll 70 to reaction solutions. Results indicate that Ficoll has surprisingly little effect on overall EcoRV reaction velocity because of offsetting increases in V(max) and K(m), and stronger nonspecific binding. The changes in measured parameters can largely be attributed to the excluded volume effects on reactant activities and the slowing of protein diffusion. Covolume reduction upon binding appears to reinforce nonspecific binding strength, and k(cat) appears to be slowed by stronger nonspecific binding, which slows product release. The data also suggest that effective Ficoll particle volume decreases as its concentration increases above a few weight percent, which may be due to Ficoll interpenetration or compression.     

Solution parameters modulating DNA binding specificity of the restriction endonuclease EcoRV

The DNA binding stringency of restriction endonucleases is crucial for their proper function. The X-ray structures of the specific and non-cognate complexes of the restriction nuclease EcoRV are considerably different suggesting significant differences in the hydration and binding free energies. Nonetheless, the majority of studies performed at pH 7.5, optimal for enzymatic activity, have found a < 10-fold difference between EcoRV binding constants to the specific and nonspecific sequences in the absence of divalent ions. We used a recently developed self-cleavage assay to measure EcoRV-DNA competitive binding and to evaluate the influence of water activity, pH and salt concentration on the binding stringency of the enzyme in the absence of divalent ions. We find the enzyme can readily distinguish specific and nonspecific sequences. The relative specific-nonspecific binding constant increases strongly with increasing neutral solute concentration and with decreasing pH. The difference in number of associated waters between specific and nonspecific DNA-EcoRV complexes is consistent with the differences in the crystal structures. Despite the large pH dependence of the sequence specificity, the osmotic pressure dependence indicates little change in structure with pH. The large osmotic pressure dependence means that measurement of protein-DNA specificity in dilute solution cannot be directly applied to binding in the crowded environment of the cell. In addition to divalent ions, water activity and pH are key parameters that strongly modulate binding specificity of EcoRV.     

Calcium buffer injections delay cleavage in Xenopus laevis blastomeres

Microinjection of calcium buffers into the two-cell Xenopus laevis embryo delays cell division in a dose-dependent manner. Four calcium buffers in the BAPTA series with different affinities for calcium were used to distinguish between a localized calcium gradient regulating cleavage and the global calcium concentration regulating this event. DibromoBAPTA (Kd = 1.5 microM) was found to delay cleavage at the lowest intracellular concentration (1.3 mM) of the four buffers tested. The effectiveness of the calcium buffers was dependent upon the buffer dissociation constant but not in a linear fashion. The concentration of buffer required to delay cleavage increased as the buffer's dissociation constant shifted above or below that of the optimum buffer, dibromoBAPTA. This relationship between a calcium buffer's effectiveness at delaying cleavage and its calcium affinity provides support for the hypothesis that a calcium concentration gradient is required for normal cell cycle progression (Speksnijder, J. E., A. L. Miller, M. H. Weisenseel, T.-H. Chen, and L. F. Jaffe. 1989. Proc. Natl. Acad. Sci. USA. 86:6607-6611). DibromoBAPTA was also injected with two different amounts of coinjected calcium to test the possibility that the free calcium concentration of the buffer solution is the important parameter for delaying cleavage. However, we found that changes in buffer concentration have a much stronger effect than changes in the free calcium concentration. This observation supports the hypothesis that BAPTA-type buffers exert their effect by shuttling calcium from regions of high concentration to those of lower concentration, reducing any calcium concentration gradients present in the Xenopus embryo.     

EcoRV restriction endonuclease binds all DNA sequences with equal affinity

In the presence of MgCl2, the EcoRV restriction endonuclease cleaves its recognition sequence on DNA at least a million times more readily than any other sequence. In this study, the binding of the EcoRV restriction enzyme to DNA was examined in the absence of Mg2+. With each DNA fragment tested, several DNA-protein complexes were detected by electrophoresis through polyacrylamide. No differences were observed between isogenic DNA molecules that either contained or lacked the EcoRV recognition site. The number of complexes with each fragment varied with the length of the DNA. Three complexes were formed with a DNA molecule of 55 base pairs, corresponding to the DNA bound to 1, 2, or 3 molecules of the protein, while greater than 15 complexes were formed with a DNA of 381 base pairs. A new method was developed to analyze the binding of a protein to multiple sites on DNA. The method showed that the EcoRV enzyme binds to all DNA sequences, including the EcoRV recognition site, with the same equilibrium constant, though two molecules of the protein bind preferentially to adjacent sites on the DNA in a cooperative fashion. All of the complexes with a substrate that contained the EcoRV site dissociated upon addition of competitor DNA, but when the competitor was mixed with MgCl2, a fraction of the substrate was cleaved at the EcoRV site. The fraction cleaved was due mainly to the translocation of the enzyme from nonspecific sites on the DNA to the specific site.     

DNA-tension dependence of restriction enzyme activity reveals mechanochemical properties of the reaction pathway

Type II restriction endonucleases protect bacteria against phage infections by cleaving recognition sites on foreign double-stranded DNA (dsDNA) with extraordinary specificity. This capability arises primarily from large conformational changes in enzyme and/or DNA upon target sequence recognition. In order to elucidate the connection between the mechanics and the chemistry of DNA recognition and cleavage, we used a single-molecule approach to measure rate changes in the reaction pathway of EcoRV and BamHI as a function of DNA tension. We show that the induced-fit rate of EcoRV is strongly reduced by such tension. In contrast, BamHI is found to be insensitive, providing evidence that both substrate binding and hydrolysis are not influenced by this force. Based on these results, we propose a mechanochemical model of induced-fit reactions on DNA, allowing determination of induced-fit rates and DNA bend angles. Finally, for both enzymes a strongly decreased association rate is obtained on stretched DNA, presumably due to the absence of intradomain dissociation/re-association between non-specific sites (jumping). The obtained results should apply to many other DNA-associated proteins.     

Positively charged C-terminal subdomains of EcoRV endonuclease: contributions to DNA binding, bending, and cleavage

The carboxy-terminal subdomains of the homodimeric EcoRV restriction endonuclease each bear a net charge of +4 and are positioned on the inner concave surface of the 50 degree DNA bend that is induced by the enzyme. A complete kinetic and structural analysis of a truncated EcoRV mutant lacking these domains was performed to assess the importance of this diffuse charge in facilitating DNA binding, bending, and cleavage. At the level of formation of an enzyme-DNA complex, the association rate for the dimeric mutant enzyme was sharply decreased by 10(3)-fold, while the equilibrium dissociation constant was weakened by nearly 10(6)-fold compared with that of wild-type EcoRV. Thus, the C-terminal subdomains strongly stabilize the enzyme-DNA ground-state complex in which the DNA is known to be bent. Further, the extent of DNA bending as observed by fluorescence resonance energy transfer was also significantly decreased. The crystal structure of the truncated enzyme bound to DNA and calcium ions at 2.4 A resolution reveals that the global fold is preserved and suggests that a divalent metal ion crucial to catalysis is destabilized in the active site. This may explain the 100-fold decrease in the rate of metal-dependent phosphoryl transfer observed for the mutant. These results show that diffuse positive charge associated with the C-terminal subdomains of EcoRV plays a key role in DNA association, bending, and cleavage.     

Reactions of BglI and other type II restriction endonucleases with discontinuous recognition sites

Type II restriction enzymes generally recognize continuous sequences of 4-8 consecutive base pairs on DNA, but some recognize discontinuous sites where the specified sequence is interrupted by a defined length of nonspecific DNA. To date, a mechanism has been established for only one type II endonuclease with a discontinuous site, SfiI at GGCCNNNNNGGCC (where N is any base). In contrast to orthodox enzymes such as EcoRV, dimeric proteins that act at a single site, SfiI is a tetramer that interacts with two sites before cleaving DNA. BglI has a similar recognition sequence (GCCNNNNNGGC) to SfiI but a crystal structure like EcoRV. BglI and several other endonucleases with discontinuous sites were examined to see if they need two sites for their DNA cleavage reactions. The enzymes included some with sites containing lengthy segments of nonspecific DNA, such as XcmI (CCANNNNNNNNNTGG). In all cases, they acted at individual sites. Elongated recognition sites do not necessitate unusual reaction mechanisms. Other experiments on BglI showed that it bound to and cleaved DNA in the same manner as EcoRV, thus further delineating a distinct group of restriction enzymes with similar structures and a common reaction mechanism.     

Divalent metal dependence of site-specific DNA binding by EcoRV endonuclease.

Measurements of binding equilibria of EcoRV endonuclease to DNA, for a series of base-analogue substrates, demonstrate that expression of sequence selectivity is strongly enhanced by the presence of Ca2+ ions. Binding constants were determined for short duplex oligodeoxynucleotides containing the cognate DNA site, three cleavable noncognate sites, and a fully nonspecific site. At pH 7.5 and 100 mM NaCl, the full range of specificity from the specific (tightest binding) to nonspecific (weakest binding) sites is 0.9 kcal/mol in the absence of metal ions and 5.8 kcal/mol in the presence of Ca2+. Precise determination of binding affinities in the presence of the active Mg2+ cofactor was found to be possible for substrates retaining up to 1.6% of wild-type activity, as determined by the rate of phosphoryl transfer. These measurements show that Ca2+ is a near-perfect analogue for Mg2+ in binding reactions of the wild-type enzyme with DNA base-analogue substrates, as it provides identical DeltaDeltaG degrees bind values among the cleavable noncognate sites. Equilibrium dissociation constants of wild-type and base-analogue sites were also measured for the weakly active EcoRV mutant K38A, in the presence of either Mg2+ or Ca2+. In this case, Ca2+ allows expression of a greater degree of specificity than does Mg2+. DeltaDeltaG degrees bind values of K38A toward specific versus nonspecific sites are 6.1 kcal/mol with Ca2+ and 3.9 kcal/mol with Mg2+, perhaps reflecting metal-specific conformational changes in the ground-state ternary complexes. The enhancement of binding specificity provided by divalent metal ions is likely to be general to many restriction endonucleases and other metal-dependent nucleic acid-modifying enzymes. These results strongly suggest that measurements of DNA binding affinities for EcoRV, and likely for many other restriction endonucleases, should be performed in the presence of divalent metal ions.     

Biotin-avidin microplate assay for the quantitative analysis of enzymatic methylation of DNA by DNA methyltransferases

An assay is described to measure methylation of biotinylated oligonucleotide substrates by DNA methyltransferases using [methyl-3H]-AdoMet. After the methylation reaction the oligonucleotides are immobilized on an avidin-coated microplate. The incorporation of [3H] into the DNA is quenched by addition of unlabeled AdoMet to the binding buffer. Unreacted AdoMet and enzyme are removed by washing. To release the radioactivity incorporated into the DNA, the wells are incubated with a non-specific endonuclease and the radioactivity determined by liquid scintillation counting. As an example, we have studied methylation of DNA by the EcoRV DNA methyltransferase. The reaction progress curves measured with this assay are linear with respect to time. Methylation rates linearly increase with enzyme concentration. The rates are comparable to results obtained with the same enzyme using a different assay. The biotin-avidin assay is inexpensive, convenient, quantitative, fast and well suited to process many samples in parallel. The accuracy of the assay is high, allowing to reproduce results within +/- 10%. The assay is very sensitive as demonstrated by the detection of incorporation of 0.8 fmol methyl groups into the DNA. Under the experimental conditions, this corresponds to methylation of only 0.03% of all target sites of the substrate. Using this assay, the DNA methylation activity of some M.EcoRV variants could be detected that was not visible by other in vitro methylation assays.     

Mode of action of topoisomerase II-targeting agents at a specific DNA sequence. Uncoupling the DNA binding, cleavage and religation events.

Methods of uncoupling the DNA binding, cleavage and religation reactions of topoisomerase II were employed to investigate the influence of topoisomerase II-directed drugs on the individual steps in the enzyme's catalytic cycle. A special DNA substrate containing a major topoisomerase II interaction site, which can be cleaved by the enzyme in the absence of any concomitant religation, was used to examine the effect of topoisomerase II-directed agents upon the DNA cleavage reaction. The experiment demonstrated that the topoisomerase II targeting agent Ro 15-0216 stimulates the DNA cleavage reaction extensively, whereas the traditional topoisomerase II inhibitor, mAMSA, has only a minor effect on this reaction. Topoisomerase II trapped in the cleavage complexes can religate to the 3' hydroxyl end of another DNA strand. Using this religation assay, it was demonstrated that the major effect of mAMSA is an inhibition of the enzyme's religation reaction, whereas Ro 15-0216 has no effect on this reaction. Recently, considerable attention has been given to drugs preventing topoisomerase II from introducing DNA cleavages. In the present paper the initial non-covalent DNA binding reaction of topoisomerase II was investigated under conditions excluding enzyme-mediated DNA cleavage. This demonstrated that the anthracycline, aclarubicin, prevents topoisomerase II from performing its initial non-covalent DNA binding reaction and thereby abolishes the DNA cleavage reaction of the enzyme. The results presented here demonstrate that profound differences exist in the mode of action of different agents targeting topoisomerase II, and that the enzyme can be affected by such agents at both its DNA binding, cleavage and religation subreactions.     

Sliding and jumping of single EcoRV restriction enzymes on non-cognate DNA

The restriction endonuclease EcoRV can rapidly locate a short recognition site within long non-cognate DNA using 'facilitated diffusion'. This process has long been attributed to a sliding mechanism, in which the enzyme first binds to the DNA via nonspecific interaction and then moves along the DNA by 1D diffusion. Recent studies, however, provided evidence that 3D translocations (hopping/jumping) also help EcoRV to locate its target site. Here we report the first direct observation of sliding and jumping of individual EcoRV molecules along nonspecific DNA. Using fluorescence microscopy, we could distinguish between a slow 1D diffusion of the enzyme and a fast translocation mechanism that was demonstrated to stem from 3D jumps. Salt effects on both sliding and jumping were investigated, and we developed numerical simulations to account for both the jump frequency and the jump length distribution. We deduced from our study the 1D diffusion coefficient of EcoRV, and we estimated the number of jumps occurring during an interaction event with nonspecific DNA. Our results substantiate that sliding alternates with hopping/jumping during the facilitated diffusion of EcoRV and, furthermore, set up a framework for the investigation of target site location by other DNA-binding proteins.     

Analysis of the recognition mechanism involved in the EcoRV catalyzed cleavage of DNA using modified oligodeoxynucleotides.

We have prepared a series of undecadeoxynucleotides that contain changes in the functional group pattern present within the EcoRV recognition site - GATATC-. Oligonucleotides were synthesized on solid phase using normal and modified beta-cyanoethylphosphoramidites and analyzed in steady state cleavage experiments with the EcoRV restriction endonuclease. The following groups appear to interact strongly with the enzyme, since their modification or substitution renders the oligonucleotides refractory to cleavage: the exocyclic NH2-groups of both A residues, the N7 of the first A residue, the exocyclic NH2-group of the C residue and the CH3-groups of both T residues. The exocyclic NH-group of the G residue supports effective recognition, since its absence lowers the kcat of the cleavage reaction. The N7 of the second A residue and the C5 position of the C residue apparently are not recognized by EcoRV; their substitution by -CH- or modification with -Br or -CH3, resp., does not considerably change the rate of cleavage. All oligonucleotides investigated compete with the unmodified substrate for binding to the enzyme. We conclude that EcoRV recognizes its substrate presumably through hydrogen bonds to the exocyclic NH2-group and the N7 of the first A residue, the exocyclic NH2-groups of the second A and the C residue, as well as through hydrophobic interactions with both T residues.     

Kinetics of target site localization of a protein on DNA: a stochastic approach

It is widely recognized that the cleaving rate of a restriction enzyme on target DNA sequences is several orders-of-magnitude faster than the maximal one calculated from the diffusion-limited theory. It was therefore commonly assumed that the target site interaction of a restriction enzyme with DNA has to occur via two steps: one-dimensional diffusion along a DNA segment, and long-range jumps coming from association-dissociation events. We propose here a stochastic model for this reaction which comprises a series of one-dimensional diffusions of a restriction enzyme on nonspecific DNA sequences interrupted by three-dimensional excursions in the solution until the target sequence is reached. This model provides an optimal finding strategy which explains the fast association rate. Modeling the excursions by uncorrelated random jumps, we recover the expression of the mean time required for target site association to occur given by Berg et al. in 1981, and we explicitly give several physical quantities describing the stochastic pathway of the enzyme. For competitive target sites we calculate two quantities: processivity and preference. By comparing these theoretical expressions to recent experimental data obtained for EcoRV-DNA interaction, we quantify: 1), the mean residence time per binding event of EcoRV on DNA for a representative one-dimensional diffusion coefficient; 2), the average lengths of DNA scanned during the one-dimensional diffusion (during one binding event and during the overall process); and 3), the mean time and the mean number of visits needed to go from one target site to the other. Further, we evaluate the dynamics of DNA cleavage with regard to the probability for the restriction enzyme to perform another one-dimensional diffusion on the same DNA substrate following a three-dimensional excursion.     

Linear diffusion of the restriction endonuclease EcoRV on DNA is essential for the in vivo function of the enzyme.

Linear diffusion along DNA is a mechanism of enhancing the association rates of proteins to their specific recognition sites on DNA. It has been demonstrated for several proteins in vitro, but to date in no case in vivo. Here we show that the restriction endonuclease EcoRV slides along the DNA, scanning approximately 1000 bp in one binding event. This process is critically dependent on contacts between amino acid residues of the protein and the backbone of the DNA. The disruption of single hydrogen bonds and, in particular, the alteration of electrostatic interactions between amino acid side chains of the protein and phosphate groups of the DNA interfere with or abolish effective sliding. The efficiency of linear diffusion is dependent on salt concentration, having a maximum at 50 mM NaCl. These results suggest that a nonspecific and mobile binding mode capable of linear diffusion is dependent on a subtle balance of forces governing the interaction of the enzyme and the DNA. A strong correlation between the ability of EcoRV mutants to slide along the DNA in vitro and to protect Escherichia coli cells from phage infection demonstrates that linear diffusion occurs in vivo and is essential for effective phage restriction.     

Effects of sodium on cell surface and intracellular 3H-naloxone binding sites

The binding of the opiate antagonist 3H-naloxone was examined in rat whole brain homogenates and in crude subcellular fractions of these homogenates (nuclear, synaptosomal, and mitochondrial fractions) using buffers that approximated intra- (low sodium concentration) and extracellular (high sodium concentration) fluids. Saturation studies showed a two-fold decrease in the dissociation constant (Kd) in all subcellular fractions examined in extracellular buffer compared to intracellular buffer. In contrast, there was no significant effect of the buffers on the Bmax. Thus, 3H-naloxone did not distinguish between binding sites present on cell surface and intracellular tissues in these two buffers. These results show that the "sodium effect" of opiate antagonist binding is probably not a function of altered selection of intra- and extracellular binding sites.     

Selective inhibition of DNA gyrase in vitro by a GC specific eight-ring hairpin polyamide at nanomolar concentration

The influence of an eight-ring hairpin DNA minor groove binder on the gyrase mediated DNA supercoiling and cleavage reaction step of the enzyme was investigated. The results demonstrate that supercoiling is affected by the hairpin polyamide in the millimolar concentration range while the enzyme catalyzed cleavage of a 162 bp fragment of pBR322 containing a single strong gyrase site is effectively inhibited at nanomolar concentration. As demonstrated by footprint analysis the latter effect is caused by a specific binding of the hairpin forming polyamide to the enzyme recognition site (GGCC), which indicates that the gyrase activity to produce a double strand break is blocked at this site. The pyrrole-imidazole hairpin polyamide is the most potent inhibitor of the gyrase mediated cleavage reaction compared to other known anti-gyrase active DNA binding agents.     

Binding kinetics and footprinting of TaqI endonuclease: effects of metal cofactors on sequence-specific interactions.

Restriction endonucleases achieve sequence-specific recognition and strand cleavage through the interplay of base, phosphate backbone, and metal cofactor interactions. In this study, we investigate the binding kinetics of TaqI endonuclease using the wild-type enzyme and a binding proficient, catalysis deficient mutant TaqI-D137A both in the absence of a metal cofactor and in the presence of Mg2+ or Ca2+. As demonstrated by gel mobility shift analyses, TaqI endonuclease requires a metal cofactor for achieving high-affinity specific binding to its cognate sequence, TCGA. In the absence of a metal cofactor, the enzyme binds all DNA sequences (TaqI cognate site, star site, and nonspecific site) with essentially equal affinity, thereby exhibiting little discrimination. The dissociation constant of the cognate sequence in the presence of Mg2+ at 60 degrees C is 0. 26 nM, a value comparable to our previously reported Km of 0.5 nM measured under steady-state conditions. The TaqI-TCGA-Mg2+ complex is stable, with a half-life of 21 min at 60 degrees C. The boundary of the protein-DNA interface is approximated to be about 18 bp as determined by DNase I footprinting. Data from this study support the notion that a metal cofactor plays a critical role for achieving sequence-specific discrimination in a subset of nucleases, including TaqI, EcoRV, and others.