Testing Kinetic Identities Involving Transition-Path Properties Using Single-Molecule Folding Trajectories

Observing the base-by-base search for native structure along transition paths during the folding of single nucleic acid hairpins

Biomolecular folding involves searching among myriad possibilities for the native conformation, but the elementary steps expected from theory for this search have never been detected directly. We probed the dynamics of folding at high resolution using optical tweezers, measuring individual trajectories as nucleic acid hairpins passed through the high-energy transition states that dominate kinetics and define folding mechanisms. We observed brief but ubiquitous pauses in the transition states, with a dwell time distribution that matched microscopic theories of folding quantitatively. The sequence dependence suggested that pauses were dominated by microbarriers from nonnative conformations during the search by each nucleotide residue for the native base-pairing conformation. Furthermore, the pauses were position dependent, revealing subtle local variations in energy-landscape roughness and allowing the diffusion coefficient describing the microscopic dynamics within the barrier to be found without reconstructing the shape of the energy landscape. These results show how high-resolution measurements can elucidate key microscopic events during folding to test fundamental theories of folding.

Pioglitazone Inhibits Metal Cluster Transfer of mitoNEET by Stabilizing the Labile Fe-N Bond Revealed at Single-Bond Level

Outer mitochondrial membrane protein mitoNEET (mNT) is a target of the type 2 diabetes drug pioglitazone. It contains a labile Fe₂S₂(His)₁(Cys)₃ metal cluster with a single Fe-N(His87) coordinating bond and can transfer its cluster to acceptor proteins. Previous ensemble studies showed that pioglitazone's binding inhibited the transfer by stabilizing the cluster, and histidine 87 may be the key mediator. Here we used atomic force microscopy-based single-molecule force spectroscopy (AFM-SMFS) to study the unfolding process of mNT dimer in the absence and presence of pioglitazone, which can distinguish the binding effect for different regions of a protein. By developing a two-step strategy using different mNT monomers with respective purification tags, we solve the problem that the classic polyprotein formation disables the mNT to dimerize. As a result, a polyprotein including a stable, naturally noncovalently bound mNT homodimer is obtained, which is required for reliable AFM measurement and pioglitazone binding. Then, the dissociation rate (k_off) of the metal cluster was measured, showing a 10-fold decrease upon pioglitazone binding, while the other parts decreased only 3-fold, verifying that pioglitazone mainly stabilizes the cluster. Moreover, when the Fe(III)-N(His87) bond was ruptured, this effect for the remaining Fe₂S₂(Cys)₃ intermediate largely disappeared. Consequently, AFM results revealed that pioglitazone inhibited the metal cluster transfer of mNT by stabilizing the labile Fe(III)-N(His87) bond. In addition, an alternative method to build a natural, noncovalently bound protein dimer or complex for reliable single-molecule measurement was developed.

Estimating transition path times and shapes from single-molecule photon trajectories: A simulation analysis

In a two-state molecular system, transition paths comprise the portions of trajectories during which the system transits from one stable state to the other. Because of their low population, it is essentially impossible to obtain information on transition paths from experiments on a large sample of molecules. However, single-molecule experiments such as laser optical tweezers or Förster resonance energy transfer (FRET) spectroscopy have allowed transition-path durations to be estimated. Here, we use molecular simulations to test the methodology for obtaining information on transition paths in single-molecule FRET by generating photon trajectories from the distance trajectories obtained in the simulation. Encouragingly, we find that this maximum likelihood analysis yields transition-path times within a factor of 2-4 of the values estimated using a good coordinate for folding, but tends to systematically underestimate them. The underestimation can be attributed partly to the fact that the large changes in the end-end distance occur mostly early in a folding trajectory. However, even if the transfer efficiency is a good reaction coordinate for folding, the assumption that the transition-path shape is a step function still leads to an underestimation of the transition-path time as defined here. We find that allowing more flexibility in the form of the transition path model allows more accurate transition-path times to be extracted and points the way toward further improvements in methods for estimating transition-path time and transition-path shape.

Cycle Completion Times Probe Interactions with Environment

Recent measurements of the durations of nonequilibrium processes provide valuable information on microscopic mechanisms and energetics. Theory for corresponding experiments to date is well-developed for single-particle systems only. Little is known for interacting systems in nonequilibrium environments. Here we introduce and study a basic model for cycle processes interacting with an environment that can exhibit a net particle flow. We find a surprising richness of cycle time variations with environmental conditions. This manifests itself in unequal cycle times τ⁺ and τ^- in forward and backward cycle directions with both asymmetries τ^- < τ⁺ and τ^- > τ⁺, speeding up of backward cycles by interactions, and dynamical phase transitions, where cycle times become multimodal functions of the bias. The model allows us to relate these effects to specific microscopic mechanisms, which can be helpful for interpreting experiments.

Probing microscopic conformational dynamics in folding reactions by measuring transition paths

Transition paths comprise those parts of a folding trajectory where the molecule passes through the high-energy transition states separating folded and unfolded conformations. The transition states determine the folding kinetics and mechanism but are difficult to observe because of their brief duration. Single-molecule experiments have in recent years begun to characterize transition paths in folding reactions, allowing the microscopic conformational dynamics that occur as a molecule traverses the energy barriers to be probed directly. Here we review single-molecule fluorescence and force spectroscopy measurements of transition-path properties, including the time taken to traverse the paths, the local velocity along them, the path shapes, and the variability within these measurements reflecting differences between individual barrier crossings. We discuss how these measurements have been related to theories of folding as diffusion over an energy landscape to deduce properties such as the diffusion coefficient, and how they are being combined with simulations to obtain enhanced atomistic understanding of folding. The richly detailed information available from transition path measurements holds great promise for improved understanding of microscopic mechanisms in folding.

Exact Solutions for Distributions of First-Passage, Direct-Transit, and Looping Times in Symmetric Cusp Potential Barriers and Wells

For a particle diffusing in one dimension, the distribution of its first-passage time from point a to point b is determined by the durations of the particle trajectories that start from point a and are terminated as soon as they touch point b for the first time. Any such trajectory consists of looping and direct-transit segments. The latter is the final part of the trajectory that leaves point a and goes to point b without returning to point a. The rest of the trajectory is the looping segment that makes numerous loops which begin and end at the same point a without touching point b. In this article we discuss general relations between the first-passage time distribution and those for the durations of the two segments. These general relations allow us to find exact solutions for the Laplace transforms of the distributions of the first-passage, direct-transit, and looping times for transitions between two points separated by a symmetric cusp potential barrier or well of arbitrary height and depth, respectively. The obtained Laplace transforms are inverted numerically, leading to nontrivial dependences of the resulting distributions on the barrier height and the well depth.

Measuring the average shape of transition paths during the folding of a single biological molecule

Transition paths represent the parts of a reaction where the energy barrier separating products and reactants is crossed. They are essential to understanding reaction mechanisms, yet many of their properties remain unstudied. Here, we report measurements of the average shape of transition paths, studying the folding of DNA hairpins as a model system for folding reactions. Individual transition paths were detected in the folding trajectories of hairpins with different sequences held under tension in optical tweezers, and path shapes were computed by averaging all transitions in the time domain, 〈t(x)〉, or by averaging transitions of a given duration in the extension domain, 〈x(t|τ)〉 _τ Whereas 〈t(x)〉 was close to straight, with only a subtle curvature, 〈x(t|τ)〉 _τ had more pronounced curvature that fit well to theoretical expectations for the dominant transition path, returning diffusion coefficients similar to values obtained previously from independent methods. Simulations suggested that 〈t(x)〉 provided a less reliable representation of the path shape than 〈x(t|τ)〉 _τ , because it was far more sensitive to the effects of coupling the molecule to the experimental force probe. Intriguingly, the path shape variance was larger for some hairpins than others, indicating sequence-dependent changes in the diversity of transition paths reflective of differences in the character of the energy barriers, such as the width of the barrier saddle-point or the presence of parallel paths through multiple barriers between the folded and unfolded states. These studies of average path shapes point the way forward for probing the rich information contained in path shape fluctuations.

Generalized Langevin Equation as a Model for Barrier Crossing Dynamics in Biomolecular Folding

Conformational memory in single-molecule dynamics has attracted recent attention and, in particular, has been invoked as a possible explanation of some of the intriguing properties of transition paths observed in single-molecule force spectroscopy (SMFS) studies. Here we study one candidate for a non-Markovian model that can account for conformational memory, the generalized Langevin equation with a friction force that depends not only on the instantaneous velocity but also on the velocities in the past. The memory in this model is determined by a time-dependent friction memory kernel. We propose a method for extracting this kernel directly from an experimental signal and illustrate its feasibility by applying it to a generalized Rouse model of a SMFS experiment, where the memory kernel is known exactly. Using the same model, we further study how memory affects various statistical properties of transition paths observed in SMFS experiments and evaluate the performance of recent approximate analytical theories of non-Markovian dynamics of barrier crossing. We argue that the same type of analysis can be applied to recent single-molecule observations of transition paths in protein and DNA folding.

Instrumental Effects in the Dynamics of an Ultrafast Folding Protein under Mechanical Force

The analysis and interpretation of single molecule force spectroscopy (smFS) experiments is often complicated by hidden effects from the measuring device. Here we investigate these effects in our recent smFS experiments on the ultrafast folding protein gpW, which has been previously shown to fold without crossing a free energy barrier in the absence of force (i.e., downhill folding). Using atomic force microscopy (AFM) smFS experiments, we found that a very small force of ∼5 pN brings gpW near its unfolding midpoint and results in two-state (un)folding patterns that indicate the emergence of a force-induced free energy barrier. The change in the folding regime is concomitant with a 30,000-fold slowdown of the folding and unfolding times, from a few microseconds that it takes gpW to (un)fold at the midpoint temperature to seconds in the AFM. These results are puzzling because the barrier induced by force in the folding free energy landscape of gpW is far too small to account for such a difference in time scales. Here we use recently developed theoretical methods to resolve the origin of the strikingly slow dynamics of gpW under mechanical force. We find that, while the AFM experiments correctly capture the equilibrium distance distribution, the measured dynamics are entirely controlled by the response of the cantilever and polyprotein linker, which is much slower than the protein conformational dynamics. This interpretation is likely applicable to the folding of other small biomolecules in smFS experiments, and becomes particularly important in the case of systems with fast folding dynamics and small free energy barriers, and for instruments with slow response times.