Abstract: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.

Extraction of kinetics from equilibrium distributions of states using the Metropolis Monte Carlo method

Biomolecular conformational changes and ligand binding: from kinetics to thermodynamics.

Frequency Adaptive Metadynamics for the Calculation of Rare-Event Kinetics.

Monte Carlo Sampling of Protein Folding by Combining an All-Atom Physics-Based Model with a Native State Bias

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

Replica exchange molecular dynamics method for protein folding simulation.

Free Energy Landscapes, Diffusion Coefficients, and Kinetic Rates from Transition Paths

Molecular free energies, rates, and mechanisms from data-efficient path sampling simulations

Constrained Hybrid Monte Carlo Sampling Made Simple for Chemical Reaction Simulations

Computing Long Timescale Biomolecular Dynamics using Quasi-Stationary Distribution Kinetic Monte Carlo (QSD-KMC)

Experiment-friendly kinetic analysis of single molecule data in and out of equilibrium

Deconvolution of dynamic mechanical networks

Nonequilibrium candidate Monte Carlo: A new tool for efficient equilibrium simulation

Identifying a Feasible Transition Pathway between Two Conformational States for a Protein

Atomic-Level Features for Kinetic Monte Carlo Models of Complex Chemistry from Molecular Dynamics Simulations

A Fast Species Redistribution Approach to Accelerate the Kinetic Monte Carlo Simulation for Heterogeneous Catalysis

The Folding Thermodynamics and Kinetics of Crambin Using an All-Atom Monte Carlo Simulation

Monte Carlo Protein Folding: Simulations of Met-Enkephalin with Solvent-Accessible Area Parameterizations

Multistate Method to Efficiently Account for Tautomerism and Protonation in Alchemical Free-Energy Calculations

Dominant Pathways in Protein Folding

Kinetic and thermodynamic analysis of proteinlike heteropolymers: Monte Carlo histogram technique