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

Single-Molecule FRET at 10 MHz Count Rates

A bottleneck in many studies utilizing single-molecule Förster resonance energy transfer is the attainable photon count rate, as it determines the temporal resolution of the experiment. As many biologically relevant processes occur on time scales that are hardly accessible with currently achievable photon count rates, there has been considerable effort to find strategies to increase the stability and brightness of fluorescent dyes. Here, we use DNA nanoantennas to drastically increase the achievable photon count rates and observe fast biomolecular dynamics in the small volume between two plasmonic nanoparticles. As a proof of concept, we observe the coupled folding and binding of two intrinsically disordered proteins, which form transient encounter complexes with lifetimes on the order of 100 μs. To test the limits of our approach, we also investigated the hybridization of a short single-stranded DNA to its complementary counterpart, revealing a transition path time of 17 μs at photon count rates of around 10 MHz, which is an order-of-magnitude improvement compared to the state of the art. Concomitantly, the photostability was increased, enabling many seconds long megahertz fluorescence time traces. Due to the modular nature of the DNA origami method, this platform can be adapted to a broad range of biomolecules, providing a promising approach to study previously unobservable ultrafast biophysical processes.

Protein Folding Dynamics as Diffusion on a Free Energy Surface: Rate Equation Terms, Transition Paths, and Analysis of Single-Molecule Photon Trajectories

The rates of protein (un)folding are often described as diffusion on the projection of a hyperdimensional energy landscape onto a few (ideally one) order parameters. Testing such an approximation by experiment requires resolving the reactive transition paths of individual molecules, which is now becoming feasible with advanced single-molecule spectroscopic techniques. This has also sparked the interest of theorists in better understanding reactive transition paths. Here we focus on these issues aiming to establish (i) practical guidelines for the mechanistic interpretation of transition path times (TPT) and (ii) methods to extract the free energy surface and protein dynamics from the maximum likelihood analysis of photon trajectories (MLA-PT). We represent the (un)folding rates as diffusion on a 1D free energy surface with the FRET efficiency as a reaction coordinate proxy. We then perform diffusive kinetic simulations on surfaces with two minima and a barrier, but with different shapes (curvatures, barrier height, and symmetry), coupled to stochastic simulations of photon emissions that reproduce current SM-FRET experiments. From the analysis of transition paths, we find that the TPT is inversely proportional to the barrier height (difference in free energy between minimum and barrier top) for any given surface shape, and that dividing the TPT into climb and descent segments provides key information about the barrier's symmetry. We also find that the original MLA-PT procedure used to determine the TPT from experiments underestimates its value, particularly for the cases with smaller barriers (e.g., fast folders), and we suggest a simple strategy to correct for this bias. Importantly, we also demonstrate that photon trajectories contain enough information to extract the 1D free energy surface's shape and dynamics (if TPT is >4-5-fold longer than the interphoton time) using the MLA-PT directly implemented with a diffusive free energy surface model. When dealing with real (unknown) experimental data, the comparison between the likelihoods of the free energy surface and discrete kinetic three-state models can be used to evaluate the statistical significance of the estimated free energy surface.