Dissecting the multistep reaction pathway of an RNA enzyme by single-molecule kinetic “fingerprinting”

  1. Shixin Liu*,
  2. Gregory Bokinsky*,
  3. Nils G. Walter, and
  4. Xiaowei Zhuang*,,§,
  1. Departments of *Chemistry and Chemical Biology and
  2. Physics, and
  3. §Howard Hughes Medical Institute, Harvard University, Cambridge, MA 02138; and
  4. Department of Chemistry, Single-Molecule Analysis Group, University of Michigan, Ann Arbor, MI 48109
  1. Edited by Robert J. Silbey, Massachusetts Institute of Technology, Cambridge, MA, and approved April 4, 2007 (received for review November 30, 2006)

  1. Fig. 1.

    Multistep catalytic reaction by the hairpin ribozyme. (A) The 2WJ hairpin ribozyme used in this study. An AC5 linker was added to the ribozyme to keep substrate or the 5′ product linked to the ribozyme. It has been shown that the ribozyme with AC5 linker behaves similarly to the nicked 2WJ ribozyme (11, 17). Substrate is colored in orange. The cleavage site is indicated by a red arrow. Biotin and fluorophores Cy3 and Cy5 were attached as indicated. (B) Multistep cleavage reaction scheme. (C) Representative time trace of Cy3 and Cy5 fluorescence and the corresponding FRET values of a ribozyme molecule undergoing cleavage. Mg2+ (12 mM) was added at 100 s. The trace was recorded at 1 frame per second (fps), and the HMM analysis was applied to the FRET trace (pink line). (D) Histogram of FRET values obtained from many time traces as shown in C, but counting only the part after Mg2+ addition. The three peaks at FRET = 0.18, 0.30, and 0.76 represent product-released (PR), undocked (UL and UC), and docked (DL and DC) states, respectively. (E) FRET histogram for ribozymes complexed with noncleavable substrate (with the native 2′-OH of A-1 substituted with 2′-OMe) in the presence of 12 mM Mg2+.


  2. Fig. 2.

    Binding and dissociation of the cleavage product 3′P. (A) Schematic of experimental principle. 3′P with fluorescence quencher dabcyl (3′P-D) was added to the immobilized and cleaved ribozyme in the presence of Mg2+. The green, red, and black dots indicate Cy3, Cy5, and dabcyl, respectively. (B) Representative time traces of Cy3 and Cy5 fluorescence of a single ribozyme molecule in the presence of 500 nM 3′P-D. The trace was recorded at the frame rate of 33 fps. In the undocked state, 3′P-D stochastically bound to and dissociated from the ribozyme, leading to repetitive quenching and dequenching of the Cy5 signal, whereas the signal from the more distant Cy3 stayed large. In the docked state, 3′P-D remained bound, quenching signals from Cy3 and Cy5, now both in proximity to the quencher. (C) Cumulative dwell-time histogram of the 3′P unbound state in the presence of 500 nM 3′P-D. The histogram (circles) is fit to a single exponential (red line), yielding an observed binding rate constant of 2.8 s−1. (D) Cumulative dwell-time histogram of the 3′P bound state (circles), which is fit to a single exponential (red line), yielding a dissociation rate constant of 2.6 s−1.


  3. Fig. 3.

    Docking kinetics of the ligated and cleaved forms of the ribozyme. (A) Reaction scheme to determine k dock L with Mg2+ added at 100 s to induce docking, shown together with a representative single-molecule FRET time trace (1 fps) and HHM analysis (pink line). k dock L is determined from the time (T L) between Mg2+ addition and docking. (B) Reaction scheme to determine k dock C with 3′P added at 100 s to induce docking, shown with a representative FRET time trace (1 fps) and HHM analysis (pink line). k dock C is determined from the time (T C) between 3′P binding and the first docking event. (C) The cumulative histogram of T L (circles), which is fit to a single exponential (red line), yielding a docking rate constant k dock L = 0.013 s−1. (D) The cumulative histogram of T C (circles), which is fit to a single exponential (red line), yielding a docking rate constant k dock C = 0.012 s−1.


  4. Fig. 4.

    Equilibrium constants for docking and the internal chemistry. (A) Schematic of the sequential buffer exchange experiment. Equilibrium among state UL, DL, DC, and UC was first reached by placing molecules in a solution containing 12 mM Mg2+ and 30 μM 3′P for >1 h. At time t 0, an unfolding buffer containing no Mg2+ or 3′P but 10 mM EDTA was added to force undocking and product release. Then a folding buffer containing 12 mM Mg2+ but no 3′P was added at t 1 to trigger docking. Molecules in state UL, DL, DC, and UC at t 0 showed distinct kinetic fingerprints, allowing them to be clearly distinguished. (B) Representative FRET time traces (1 fps) for each of the four scenarios. The first and second buffer exchanges were at 0 and 100 s, respectively. Scenarios 1, 2, 3, or 4 correspond to a molecule at state UL, DL, DC, or UC at t 0, respectively.


  5. Fig. 5.

    Dwell-time analysis of the last docking event before 3′P release. (A) Reaction scheme with Mg2+ added at 100 s to induce docking and cleavage, shown with a representative FRET time trace (1 fps). The HMM analysis is shown in a pink line. The dwell time T of the last FRET = 0.8 event before releasing 3′P was analyzed. (B) The cumulative histogram of T (circles) is fit to a single exponential (red line) with a decay constant of m = 0.0032 s−1, where m is a function of k undock L, k undock C, k lig, and k cleav.


  6. Fig. 6.

    The overall cleavage kinetics catalyzed by the 2WJ hairpin ribozyme. (A) Summary of the rate constants along the reaction pathway for the major (stably docked) population of the enzyme determined in this work. (B) (Upper) Number of molecules transitioned to the PR state per 30-s interval (N) as a function of time with Mg2+ being added at t = 0. Red curve indicates data simulated (N sim,i) using the kinetic rate constants determined above without any adjustable parameter and black curve indicates data obtained experimentally (N exp,i). (Inset) Integrated number of molecules accumulated in the PR state as a function of time (simulated curve: red; experimental curve: black). (Lower) The residual plot showing the difference between the simulated and experimental data. A standardized residual is calculated as ( Formula.


Footnotes

  • To whom correspondence should be addressed. E-mail: zhuang{at}chemistry.harvard.edu
« Previous | Next Article »Table of Contents