Single-molecule Taq DNA polymerase dynamics

Taq DNA polymerase functions at elevated temperatures with fast conformational dynamics—regimes previously inaccessible to mechanistic, single-molecule studies. Here, single-walled carbon nanotube transistors recorded the motions of Taq molecules processing matched or mismatched template–deoxynucleotide triphosphate pairs from 22° to 85°C. By using four enzyme orientations, the whole-enzyme closures of nucleotide incorporations were distinguished from more rapid, 20-μs closures of Taq’s fingers domain testing complementarity and orientation. On average, one transient closure was observed for every nucleotide binding event; even complementary substrate pairs averaged five transient closures between each catalytic incorporation at 72°C. The rate and duration of the transient closures and the catalytic events had almost no temperature dependence, leaving all of Taq’s temperature sensitivity to its rate-determining open state.


INTRODUCTION
DNA polymerases catalyze reactions essential to life-DNA replication and repair (1). They incorporate incoming deoxynucleotide triphosphates (dNTPs) into a nascent DNA strand that is complementary to a single-stranded DNA (ssDNA) template. This capability makes DNA polymerases workhorses for molecular biology and biotechnology (2). In particular, the thermostable DNA polymerase from Thermus aquaticus (Taq) is widely used in the polymerase chain reaction (PCR) (3).
The simplest model for DNA polymerase function is a catalytic cycle between open and closed conformations (1,4). In the open conformation, dNTP and ssDNA template bind to sites on the enzyme's "fingers" and "thumb" domains, respectively (5). After substrate binding, the two domains come together to form a closed conformation that either rejects the dNTP or catalytically incorporates it into the nascent ssDNA (4,6). While this simple model is broadly accepted, fundamental questions about the fidelity-checking mechanism have motivated searches for additional intermediate conformations (7,8). How does fidelity checking occur in the open conformation before the dNTP and DNA template can coordinate? Alternatively, if fidelity checking is postponed until the catalytically active, closed conformation, then how are misincorporation rates of 10 −6 errors per dNTP achieved? Early research hypothesized the necessity of partially closed, intermediate conformations where dNTP identification and mismatch rejection could occur before the closed conformation (6,9).
Over the past decade, multiple breakthroughs began to identify and characterize these intermediate conformations in Taq and homologous A-family DNA polymerases. In 2011, Wu and Beese (10) crystallized an "ajar" structure using a noncomplementary pairing and mutagenic modifications to stabilize the conformation. Independently, advances in single-molecule Förster resonance energy transfer (smFRET) inferred partially closed conformations from fluorescence emission populations (8,11,12) and developed a dynamic model of rapidly interconverting conformational populations influenced by complementarity (13,14). On the theoretical side, all-atom molecular dynamics (MD) advanced to produce the first detailed trajectories of closing and reopening transitions (15,16), followed by observations of interconversion dynamics among a diverse ensemble of atomic configurations (17,18). Together, this progress has revealed partially closed intermediates and dynamics between the open and closed conformations, but it has not directly observed the fidelitychecking mechanism nor resolved the debate over whether "conformational selection," "induced fitting," or a combination of the two is the best reductionist model (8,(19)(20)(21)(22)(23)(24)(25)(26). A broad diversity in the MD trajectories suggests that fidelity checking is more complex than a binary choice between two deterministic, concerted pathways.
Here, we contribute long-duration, single-molecule trajectories of Taq opening and closing upon either matched or mismatched dNTP-template pairs. Using a solid-state electronic technique with microsecond resolution, we monitored individual Taq molecules over minutes of activity, directly recording processive catalytic dNTP incorporations interspersed among 20-s, noncatalytic closures. The results agree with the glimpses provided by smFRET measurements (11,12,14) and MD trajectories (17), detailing the timing with which Taq tests for complementarity by closing upon matched or mismatched dNTPs. A further advantage of the electronic technique is its compatibility with a wide range of temperatures; here, Taq's singlemolecule activity was studied from 22° to 85°C to cover the range relevant for PCR.
To record these trajectories, we attached Taq molecules to fieldeffect transistors (FETs) composed of individual, single-walled carbon nanotubes (SWNTs; Fig. 1A). The technique has proven useful with other enzymes (27,28), including the Klenow fragment (KF) of DNA polymerase (29)(30)(31). Following methods similar to our previous work (29) and adopting the single-cysteine variants of Taq studied by Xu et al. (32), we used a pyrene-maleimide linker to bioconjugate individual Taq molecules to SWNT devices in four distinct orientations designed to transduce signals from the backside of Taq's intervening (R411C), palm (A814C), fingers (R695C), or thumb (E524C) domains (figs. S1 to S5). Here, the resultant enzyme-SWNT devices are each referred to as Taq domain (e.g., Taq-Intervening). Atomic force microscopy of an example Taq-Intervening device is shown in Fig. 1B and fig. S6. In each orientation, the movement of charged amino acid side chains closest to the SWNT electrostatically modulated the source-drain current I(t), whereas more distant domains were screened by mobile ions in the Taq activity buffer [40 mM Hepes, 50 mM KCl, 5 mM MgCl 2 (pH 8.5)], which had a Debye screening radius of 1.2 nm ( fig. S2). The electrical current as a function of time, I(t), was recorded while Taq interacted with a solution of homopolymeric "polyT" DNA primer-template [4 nM poly(dT) 42 fused to an M13 priming site and annealed with M13 primer; table S1] and/or different dNTPs (typically 10 M). The I(t) signals presented here were high-pass filtered at 15 Hz to remove DC components of 10 to 500 nA and the lowest-frequency fluctuations from the raw I(t) data. Materials and Methods provides a fuller description of the experimental methods for production and purification of each fulllength Taq variant and SWNT FET fabrication, bioconjugation, and measurement.

RESULTS
Aside from being acquired at 72°C, the signals generated by Taq-Intervening devices were essentially identical to previous electrical measurements of KF (33). In activity buffer (Fig. 1C) or with mismatched dNTPs (Fig. 1D), I(t) was a featureless band of 1/f noise. In the presence of matched dNTPs, additional two-level switching appeared ( Fig. 1E and fig. S7). At 72°C, the average rate of Taq's two-level excursions below the baseline was 61 s −1 , effectively equal to Taq's catalytic rate k cat at that temperature (34). Led by this agreement and the similarity to previous KF results (29), where event counting correctly enumerated DNA primer-template lengths (31), we assigned I(t) to a two-state model of dNTP incorporations. The I(t) = 0 baseline corresponded to Taq's open conformation, and each excursion recorded one dNTP incorporation occurring in Taq's closed conformation. For brevity, we henceforth refer to these excursions as "catalytic closures." Temperature dependence provided additional confirmation of the assignment of I(t) excursions to Taq's catalytic closures. From 22° to 85°C, the average rate of catalytic closures increased exponentially from 4 to 96 s −1 (Fig. 2). While exact values varied from one molecule to another, the range of rates and the inferred activation energy of 11 kcal/mol (shown in fig. S8) was comparable to previous reports of k cat (36)(37)(38)(39). As in ensemble measurements (36,38), we observed a softening of the Arrhenius rate law above 50°C, where thermal fluctuations decrease Taq-dNTP and primer-template binding stabilities. At 95°C, I(t) became disordered and not enumerable, but organized activity returned upon cooling devices below 85°C. Figures 1 and 2 are short samples taken from continuous I(t) records lasting no less than 5 min and, in some cases, as long as 45 min. At PCR temperatures, the full recordings contained thousands of catalytic closures per minute, each of which could be resolved and individually characterized in the manner depicted in Fig. 2B. For quantitative statistical analysis, each catalytic closure was assigned a duration  cat , a waiting time before the closure  open , and an average amplitude <I> measured from the event's instantaneous baseline. Transition times at the front and back step edges of each catalytic closure were also inspected but determined to be limited by the amplifier rise time (0.5 s). Unlike previous work with T4 lysozyme, which directly observed finite transition times averaging 37 s (40), most Taq open-to-closed transitions had sharp edges and rise times faster than 1 s. Figure 3 shows probability distribution functions for large populations of  open and  cat from Taq-Intervening activity at four temperatures.  open exhibited strong temperature dependence, consistent with previous ensemble work (6), whereas  cat was largely temperature independent. k cat was essentially determined by the former, because  open was 100 to 1000 times longer than  cat . Consequently, Fig. 3A is a single-molecule representation of Taq's k cat and its temperature dependence, illustrating the single-molecule, event-by-event composition of k cat and its instantaneous variability. Solid lines in Fig The  open distributions exceeded the exponential fits at longer durations, but this deviation occurred over many minutes of acquisition. Shorter subpopulations of  open fit single exponentials with minute-by-minute variations in slope. Similar variability was observed in  cat , causing its cumulative distributions to resemble double-exponential curves (Fig. 3B). This artificial stretching of single-exponential distributions can indicate external factors like temperature variation or dynamic disorder intrinsic to an enzyme (41)(42)(43)(44), but the latter explanation seems more likely because  cat had no other temperature sensitivity. Table 1 summarizes numeric values for the fits in Fig. 3 with the range (±) of single-exponential  cat fits that were observed experimentally at each temperature during extended recordings. Data from the ends of these ranges are responsible for the long tails in the cumulative distributions shown in Fig. 3.
Taq-Palm devices corroborated the findings from Taq-Intervening. The Taq-Palm orientation suffered from sparse activity, perhaps inhibited by its attachment to the SWNT FET, but otherwise minutelong bursts of I(t) excursions had the same characteristics as Taq-Intervening ( fig. S9). Even studied in detail, however, the singlemolecule trajectories from Taq-Palm and Taq-Intervening devices only observed k cat activity already well established by ensemble and smFRET studies. The palm and intervening domains participated in the enzyme's global transformations during catalysis (32), but these attachment sites were located too far from the incoming dNTPs to transduce the dynamic activity associated with complementarity checking ( fig. S2). Studying those signals required new orientations that attached the SWNT FET directly to the fingers or thumb domains.
Taq-Fingers devices produced I(t) signals with two distinct populations of events. The first population was identical to the Taq-Intervening catalytic closures, having the same  cat probability distribution and thermally activated rates k cat , and only occurring in polyT solutions with matched dNTPs [e.g., 2′-deoxyadenosine 5′-triphosphate (dATP)]. The second population of events, on the other hand, appeared in I(t) signals collected in both matched (Fig. 4A) and mismatched dNTPs (Fig. 4B). This population of events, which we term "transient closures," was also distinguished by an average duration < transient > = 18 ± 6 s, approximately eight times shorter than < cat >. In mismatched dNTPs (where transient closures were easiest to analyze because catalytic closures were absent), transient closures occurred at a mean rate k transient = 10 ± 8 s −1 with a singleexponential distribution of durations (Fig. 4C). Datasets acquired at 45°, 60°, and 75°C could not resolve any temperature dependence in the characteristics of these transient closures (fig. S10). As with catalytic closures, observing transient closures required the presence of both DNA primer-template and dNTP; these events disappeared in control measurements using buffer, polyT alone, or dNTPs alone (fig. S11). Similar signals were generated by Taq-Thumb devices (fig. S12).
In matched dNTPs, Taq-Fingers signals contained a mixture of the new transient closures along with the catalytic closures observed with Taq-Intervening. The two types of events overlapped in a biexponential distribution (Fig. 4C), making classification of individual events impossible in the 75-to 125-s region where the two subpopulations had comparable probabilities. Overall, though, the subpopulations were distinguishable and easily identified with help from the single-population datasets. The longer-duration subpopulation of   Fig. 4C shows the matched and mismatched population distributions in units of events per second. As noted above, single-exponential fits to the  cat and  transient distributions disguised the actual dynamics and variability in k cat and k transient . When calculated on 1-s intervals from single-molecule trajectories, both rates varied within a tenfold range. Figure 5 shows instantaneous k transient values (Fig. 5A) and the moment-to-moment accumulation of those values into log-normal probability distributions, with well-separated distributions for matched and mismatched dNTPs (Fig. 5B). The log-normal fitting indicates that k transient obeyed a thermodynamic rate equation ln(k transient /k 0 ) = −G ‡ /k B T, where G ‡ is the energy barrier for transient-closing motions, k 0 is an attempt rate, and k B is Boltzmann's constant. Plotting k transient on a relative energy axis G ‡ illustrates two important features. First, the distribution widths full width at half maximum = 0.65 ± 0.05 kcal/mol equaled 1 k B T, showing that thermal fluctuations in G ‡ governed the k transient variability. Figure S13 expands upon Fig. 5 with a comparison of k transient variability at 45°, 60°, and 72°C, demonstrating that the ~k B T width exceeded all other temperature dependence.
The second important feature of the G ‡ axis in Fig. 5 is the energetic interpretation it provides for k transient . The 30-fold difference between matched and mismatched k transient values is energetically summarized as a peak-to-peak energy difference of 2.1 ± 0.2 kcal/mol. This separation exactly matches the difference in binding affinities revealed by ensemble assays two decades ago (6,9,45,46). Thus, the preferential binding of matched dNTPs, well-accepted and established by ensemble measurements, reappeared with identical energetics when thousands of single-molecule events were fit to log-normal curves in Fig. 5.
Despite this one-to-one correspondence between dNTP binding and transient closures, the electronic technique did not sense dNTP binding directly. No I(t) excursions were produced in the absence of template, although Taq binds dNTP in the absence of template. Therefore, we conclude that any I(t) excursion observed electronically was not a binding event per se, but rather a closure of the Taqtemplate-dNTP complex induced subsequent to binding. Matched dNTPs bind at rates 30 times higher than unmatched dNTPs, and that difference manifested itself in the single-molecule data as a 30-fold increase of transient closures. This linear correspondence further suggests that mismatched dNTPs dissociate after a single transient closure. An efficient mismatch-rejection mechanism might operate in that manner to avoid multiple closure attempts on the same dNTP.
Other than the differences in rates and durations already noted, no other characteristics of the I(t) excursions distinguished matched from mismatched dNTPs or transient closures from catalytic ones. Three hypotheses were given special attention in the data analysis. First, the  transient and  cat events generated by Taq-Fingers were uniformly intermixed in dATP without evidence of bunching or sequencing, as if both were independent, memoryless Poisson processes. For instance,  cat events observed with Taq-Fingers had identical timing and statistics to those observed with Taq-Intervening. Second, the  transient durations were the same for matched dATP and mismatched 2′-deoxyguanosine 5′-triphosphate (dGTP), ruling out event duration as a possible direct measurement of any fidelity-checking step. Third, the transient and catalytic closures had identical distributions of the signal amplitude observed for Taq-Fingers (Fig. 6). This latter observation is important because <I> is electrostatically transduced by motions of charged amino acids near the attachment site. Because <I> did not distinguish matched and mismatched dNTPs or catalytic closures from transient closures, the SWNT FET was registering the same degree of amino acid motion in all three cases. On the other hand, <I> differed between Taq-Fingers, Taq-Intervening, and the other two orientations, reflecting the different compositions of signalgenerating amino acids at each attachment site.

DISCUSSION
Historically, the preferential binding of matched dNTPs has generated many questions about fidelity-checking mechanisms in polymerases. Before binding dNTP, Taq sits in an open conformation with the templating base far from the dNTP binding site on the fingers domain. Researchers have questioned how a preference for one dNTP over another arises and transmits to the fingers domain (6,9,45). Ultrashort finger-closing transients or partially closed intermediate conformations have both been proposed as possible keys to such a mechanism (12). However, the dynamical time scales of these transients have been proven too brief for direct observation by smFRET (13, 14, 47). Modern MD simulations have reproduced the 2 kcal/mol preference for binding matched dNTPs (15,17,18), but only in closed or partially closed conformations and with submicrosecond time scales. Nevertheless, transients into and through intermediates have grown to be widely anticipated feature of Taq's dynamics.
In the present work, the observations described above lead us to conclude that I(t) was not sensitive to Taq's partially closed conformations, especially any conformation key to fidelity checking (12). Instead, fidelity checking completely preceded the motions detected by I(t), leading to the one-to-one correspondence between binding and closure and the lack of any <I> difference between catalytic closures and transient ones. Close inspection of Taq's structure in the Taq-Fingers orientation suggests that one charged amino acid (D655) generated the I(t) excursions but only when Taq closed fully and not in the partially closed conformation reported by Wu and Beese (10). Therefore, we conclude that I(t) excursions reported the fully closed conformation and that I(t) was not sensitive to partially closed conformations. Distinguishing those may require carefully designed mutations in future experiments, as well as higher bandwidths if the event durations are less than 1 s.
Even  = 5). No model has suggested that polymerase might pass a fidelity checkpoint with a matched dNTP and yet still proceed through multiple noncatalytic closures. When the temperature was lowered, transient closures continued to occur at k transient = 280 s −1 , while the number of catalytic closures 1/k cat decreased exponentially, filling the entire  open duration with transient closures. This free-running process requires revisions to the linear, multistep diagrams often used to represent Taq's catalytic cycle (9,48).
Second, transient closures were observed upon binding mismatched dNTPs. Again, the premise of a fidelity checkpoint is to reduce misincorporations by reducing closures upon mismatched dNTPs. Here, Taq appears to have fully closed once upon every dNTP binding event, whether matched or mismatched. Closing upon mismatched dNTPs after fidelity checking may be counterintuitive, but that kind of motion has been inferred from the analysis of smFRET transition probabilities (14) and directly seen in MD trajectories (17). The MD research, in particular, has suggested that the mismatch-rejection mechanism might need to follow a trajectory in which Taq briefly accesses the closed conformation after the fidelity check and on the way to reopening and dNTP dissociation. In this type of cycle, each mismatched dNTP binding event is followed by one transient closure with a oneto-one correspondence as reported here.
Last, the observed dynamics lead us to propose a previously undescribed sequence in which matched dNTPs also dissociate after every transient closure. While expected for mismatched dNTPs, frequent dissociation of matched dNTPs is contrary to the optimization of catalytic rate. Nevertheless, our proposal is based on two observations. First, if Taq executes one transient closure per binding and dissociation cycle with mismatched dNTP, then it must perform a similar cycle with matched dNTPs to produce a k transient ratio that matches the ratio of ensemble binding rates. Otherwise, if matched dNTPs remained bound through multiple transient closures, then the k transient ratio would be multiple times the 30:1 ratio of binding events. Second, the I(t) recordings with matched dNTPs show that transient closures filled  open and occurred at a temperature-insensitive, constant rate k transient . These data suggest a rapid and continuous cycle of dNTP binding, closure, and dissociation that repeats indefinitely while awaiting the completion of some other, thermally activated and rate-limiting step.
This type of continuous cycling does not fit the common schemes used to represent Taq's catalytic cycle with matched dNTPs. In most of these schemes (9,48), dNTP binding is one of a linear series of steps that occurs during  open and that proceeds directly to catalysis after successful fidelity checking. If  open is instead filled with a freerunning, parallel cycle of dNTP binding, closure, and dissociation, then that cycle could require rethinking the schemes for Taq's steps and revisiting the interpretation of ensemble experiments where rate constants were extracted by assuming a linear scheme.
Transient closures with matched dNTPs may reflect misalignment of the molecular orbitals required for bond formation and breakage. The observation of multiple transient closures per catalytic closure may reflect some of Taq's dynamical interconversions as it previews the closed conformation and cycles back to open or partially open conformations to achieve correct template positioning, dNTP orientation, or alignment of active site and substrate functionalities, including Mg 2+ . On the other hand, the  transient distributions with matched and mismatched dNTPs had identical durations and they both fit single-step Poisson distributions, so any rearrangement or repositioning was not occurring during the transient closure. Furthermore, closures lasting less than 5 s were not counted here, leaving many additional opportunities for another population of the submicrosecond dynamics predicted by MD simulations (16,18). Once the molecular orbital alignment is correct, a catalytic closure can proceed (7) with the global conformational changes that extend to the palm and intervening domains (32). Figure S14 depicts a likely sequence of transient closures and catalytic closures among the mechanistic steps of the catalytic sequence.
This work opens frontiers in protein engineering and biophysics. The conceptual separation of Taq's catalytic cycle into a thermally activated step and temperature-independent motions offers a successful paradigm for the evolution of high-temperature enzymes. Directed evolution of industrially important enzymes could benefit from exploiting this paradigm to achieve Taq's catalytic rate acceleration at elevated temperatures. Such capabilities are notoriously hard to engineer by standard experimental and computational approaches (49,50). Furthermore, the dynamics of DNA replication is merely one case of an important but complex chemical process challenging state-of-the-art biophysical techniques. More broadly, enzyme evolution, optimization, and engineering could all benefit from the ability to track single-molecule chemical trajectories indefinitely with microsecond resolution. The added benefits of operating over a wide temperature range or with solvent compatibility make the solid-state technique of single-molecule electronics a promising approach to reveal timing, sequence, and dynamic kinetics in a wide range of systems.

Mutagenesis of Taq
A pET28c plasmid for ligation-independent cloning containing an open reading frame (ORF) encoding Taq fused to an N-terminal His 6 peptide epitope was mutated to abolish the polymerase's exonuclease activity (G46D) and introduce a single cysteine (R411C, E524C, R695C, or A814C) (table S2). Q5 site-directed mutagenesis was applied according to the manufacturer's instructions using Q5 Hot Start High-Fidelity DNA Polymerase and the appropriate oligos, Oligo1 to Oligo10 (table S1). The Kinase, Ligase, and DpnI mix (KLD mix) (5 l) was transformed into TOP10 Escherichia coli-competent cells, and transformants were plated on a kanamycin-supplemented (50 g/ml) agar plate before incubation at 37°C overnight. Five single colonies were selected to inoculate 4 ml of 2YT medium in a 15-ml culture tube supplemented with kanamycin (50 g/ml). The seed cultures were incubated at 37°C with shaking at 225 rpm for 8 to 12 hours. Plasmid DNA was isolated using the QIAprep spin miniprep kit according to the manufacturer's instructions. The successful subcloning of the ORF of the Taq variants was checked via DNA sequencing using the appropriate oligos, Oligo11 to Oligo14 (table S1).

Protein expression and purification
Protein expression and purification was conducted using procedures previously described (30,51), with modifications as follows. The plasmid encoding each Taq variant was transformed via heat shock into competent BL21 Star (DE3) E. coli. A single colony was transferred to LB medium (25 ml) supplemented with kanamycin (50 g/ml) and incubated at 37°C for 10 to 16 hours. An aliquot of the starter culture (10 ml) was transferred to LB medium (1 liter of LB in a 4-liter baffled flask). After reaching an optical density at 600 nm of ~0.4, the culture was induced through addition of IPTG (1 mM final concentration) before incubation, with shaking (225 rpm) at 37°C for 2 hours.  S4). For single-molecule measurements on SWNT devices, the purified recombinant protein was dialyzed using 3.5K MWCO 0.1-ml Slide-A-Lyzer mini dialysis device into attachment buffer [40 mM Hepes, 50 mM KCl, 5 mM MgCl 2 (pH 6.5)]. The protein concentration was determined by a BCA assay using the enzyme's estimated MW (www.expasy.org). For storage, the purified recombinant protein was dialyzed into storage buffer [20 mM tris-HCl, 100 mM KCl, 0.1 mM EDTA, 1 mM dithiothreitol, 0.5% Tween 20, 10% glycerol (pH 7.4)] and maintained to −80°C. Ensemble activity assays Before bioconjugation to SWNTs, Taq activity was examined using previously described ensemble activity assays (30,51), using Oligo15 and polyT or Oligo16 (table S1). The high-performance liquid chromatography-purified oligos were solubilized in water to 100 M and hybridized to the template through annealing from 45° to 5°C with a 5°C temperature decrease every 5 min. The positive control reactions contained the Taq variant (2 M), dATP (100 M), and DNA template-primer (5 M) in 1× Taq Reaction Buffer [20 mM tris-HCl, 50 mM KCl, 5 mM MgCl 2 (pH 8.3)]. The negative control reaction omitted Taq. Reactions were incubated at 72°C for 1 hour, before the addition of SYBR Green Nucleic Acid I stain, and then visualized by electrophoresis using 5% high-resolution agarose gel. Gels were imaged using a Typhoon scanner at 256 nm. A representative ensemble activity assay is shown in fig. S5.

Fabrication and measurement of SWNT FETs
SWNT FET devices were fabricated using techniques described previously (29). In brief, SWNTs were grown via chemical vapor deposition (CVD) from Fe 30 Mo 84 catalyst seeds (52). The Fe 30 Mo 84 catalyst nanoparticles were confined to Al 2 O 3 islands photolithographically defined on a 4″ p++ silicon wafer with a 250-nm thermal oxide (Silicon Quest International). A fast-heating CVD recipe facilitated growth of approximately one straight, 100-m-long SWNT per island. Immediately after growth, the wafer was baked in air (315°C, 30 min) to remove excess amorphous carbon, annealed (95%/5% flowing Ar/H 2 at 600°C, 60 min), and given a 15-nm protective Al 2 O 3 coating by atomic layer deposition (Savannah S200, Cambridge NanoTech). Electrodes were defined by lithographically patterning a bilayer resist (S1808 on LOR-A1), selectively etching the Al 2 O 3 (Transetch-N, Transene), e-beam evaporating 10-nm Pt and 40-nm Ni, and lifting off the resist (Remover PG, MicroChemicals). Devices were then passivated with poly(methyl methacrylate) (PMMA), with a small opening defined over each 2-m SWNT channel by e-beam lithography. In the preferred protocol, devices were individually etched in Transetch-N to remove the protective Al 2 O 3 in this windowed region immediately before biofunctionalization.
Biofunctionalization of the newly exposed channel proceeded in two stages. First, the device was soaked in a solution of 100 pM pyrene-maleimide diluted in 1 M of pyrene in ethanol for 2 min to coat the SWNT surface with a monolayer of pyrene and pyrene-maleimide linkers. Excess pyrene was then gently rinsed away with ethanol for 2 min in a drop-by-drop fashion. Next, the device was submerged in flowing attachment buffer [40 mM Hepes, 50 mM KCl, 5 mM MgCl 2 (pH 6.5)] for 10 min to equilibrate the surface. For bioconjugation, the device was incubated in a solution of the same buffer with 4 nM Taq for 5 min, followed by 10 min of continuous rinsing with buffer to remove nonspecific protein. This protocol was optimized to yield approximately one Taq attachment per 2 m of SWNT, which in our geometry equaled one active, single-molecule SWNT FET device for every pair of devices prepared ( fig. S6). Figure S15 summarizes the measurements done with all four Taq variants and the types of signals observed in each case.
Before measurements, devices were mounted in a holder, wirebonded, and compressed with a polydimethylsiloxane (PDMS) gasket. To achieve a proper seal between the PDMS and PMMA surface of the device, the PDMS gasket was sonicated in 80% ethanol and 20% isopropanol for 10 min before being clamped in place. A 100-m-wide microfluidic channel in the PDMS confined the liquid solutions to flow over the exposed SWNT FET device. The potential of the liquid was controlled with Pt counter electrodes and Ag/AgCl reference electrodes mounted in the fluid inlets and outlets.
During electrical measurements, electrodes were controlled and monitored using National Instruments data acquisition hardware (PCI-6281) and LabVIEW 2018 software. A transimpedance current preamplifier with 10 7 V/A gain (DHPCA-100 with 0.55-s rise time, or DLPCA-200 with 2.5-s rise time, FEMTO) was used to bias the drain electrode (20 to 100 mV) and measure the source-drain current I(t). The Pt counter electrode was set to electrostatically gate the SWNT near the middle of its linear operating range (typically −200 mV with respect to Pt), and then, signals were continuously digitized and stored at 2 MHz.

Statistical analysis
To isolate and analyze enzyme-induced I(t) signals, a denoising filter was implemented with the LabVIEW 2018 wavelet analysis library. Two parallel, nonadaptive bandpass filters used Haar wavelets with garrote thresholding from 66 ms to 256 s for a low-noise output channel or from 66 ms to 8 s for a sharper-edged, high-frequency output channel. Figure S7 shows an example of raw I(t) and the resulting filtered, high-frequency output channel. After denoising, an adaptive thresholding filter was applied to both channels to identify individual events, excluding fluctuations with durations ≤4 s to avoid counting external noise sources.
After identification, each event was individually characterized. In addition to the event duration and mean event amplitude reported here, tabulated parameters included total variation, kurtosis, and frequency components calculated from both the raw signal and the denoised outputs. No features distinguished transient closures from catalytic closures except for their durations, and the duration distribution was smooth and continuous when both present in a dataset. Categorization into transient and catalytic subsets was therefore guided by biexponential fitting and comparison against single-exponential distributions, as depicted in Fig. 4C. For quantitative statistical analysis, nonoverlapping subsets were artificially divided into events lasting <75 s (brief events) or >125 s (long events), ignoring the indeterminate events in between. This definition of subsets introduced systematic undercounting of event frequencies. It also limited our pursuit of quantitative comparisons among the different Taq mutants.
We note that the 4-s noise cutoff caused the analysis to ignore very short fluctuations and to introduce systematic error into the enumeration of rates. The cutoff truncates no less than 20% of events in a Poisson process having a mean duration <> = 18 s, which was the experimental value for transient closures. The same cutoff truncates only 3% of events in a Poisson process with <> = 150 s, which was a typical value for catalytic closures. All k transient and k cat rates stated here are based on events enumerated by the analysis, without any adjustment to account for this systematic undercounting.