Barrier Crossing Dynamics from Single-Molecule Measurements

The significance of fuzzy boundaries of the barrier regions in single-molecule measurements of failed barrier crossing attempts

A recent ground-breaking experimental study [Lyons et al., Phys. Rev. X 14(1), 011017 (2024)] reports on measuring the temporal duration and the spatial extent of failed attempts to cross an activation barrier (i.e., "loops") for a folding transition in a single molecule and for a Brownian particle trapped within a bistable potential. Within the model of diffusive dynamics, however, both of these quantities are, on average, exactly zero because of the recrossings of the barrier region boundary. That is, an observer endowed with infinite spatial and temporal resolution would find that finite loops do not exist (or, more precisely, form a set of measure zero). Here we develop a description of the experiment that takes the "fuzziness" of the boundaries caused by finite experimental resolution into account and show how the experimental uncertainty of localizing the point, in time and space, where the barrier is crossed leads to observable distributions of loop times and sizes. Although these distributions generally depend on the experimental resolution, this dependence, in certain cases, may amount to a simple resolution-dependent factor and, therefore, the experiments do probe inherent properties of barrier crossing dynamics.

Arrhenius law for interacting diffusive systems

Finding the mean time it takes for a particle to escape from a metastable state due to thermal fluctuations is a fundamental problem in physics, chemistry, and biology. Here, we consider the escape rate of interacting diffusive particles, from a deep potential trap within the framework of the macroscopic fluctuation theory-a nonequilibrium hydrodynamic theory. For systems without excluded volume, our investigation reveals adherence to the well-established Arrhenius law. However, in the presence of excluded volume, a universality class emerges, fundamentally altering the escape rate. Remarkably, the modified escape rate within this universality class is independent of the interactions at play. The universality class, demonstrating the importance of excluded volume effects, may bring insights to the interpretation of escape processes in the realm of chemical physics.

Understanding the heterogeneity intrinsic to protein folding

Relating the native fold of a protein to its amino acid sequence remains a fundamental problem in biology. While computer algorithms have demonstrated recently their prowess in predicting what structure a particular amino acid sequence will fold to, an understanding of how and why a specific protein fold is achieved remains elusive. A major challenge is to define the role of conformational heterogeneity during protein folding. Recent experimental studies, utilizing time-resolved FRET, hydrogen-exchange coupled to mass spectrometry, and single-molecule force spectroscopy, often in conjunction with simulation, have begun to reveal how conformational heterogeneity evolves during folding, and whether an intermediate ensemble of defined free energy consists of different sub-populations of molecules that may differ significantly in conformation, energy and entropy.

Fick-Jacobs description and first passage dynamics for diffusion in a channel under stochastic resetting

The transport of particles through channels is of paramount importance in physics, chemistry, and surface science due to its broad real world applications. Much insight can be gained by observing the transition paths of a particle through a channel and collecting statistics on the lifetimes in the channel or the escape probabilities from the channel. In this paper, we consider the diffusive transport through a narrow conical channel of a Brownian particle subject to intermittent dynamics, namely, stochastic resetting. As such, resetting brings the particle back to a desired location from where it resumes its diffusive phase. To this end, we extend the Fick-Jacobs theory of channel-facilitated diffusive transport to resetting-induced transport. Exact expressions for the conditional mean first passage times, escape probabilities, and the total average lifetime in the channel are obtained, and their behavior as a function of the resetting rate is highlighted. It is shown that resetting can expedite the transport through the channel-rigorous constraints for such conditions are then illustrated. Furthermore, we observe that a carefully chosen resetting rate can render the average lifetime of the particle inside the channel minimal. Interestingly, the optimal rate undergoes continuous and discontinuous transitions as some relevant system parameters are varied. The validity of our one-dimensional analysis and the corresponding theoretical predictions is supported by three-dimensional Brownian dynamics simulations. We thus believe that resetting can be useful to facilitate particle transport across biological membranes-a phenomenon that can spearhead further theoretical and experimental studies.

Learning continuous potentials from smFRET

Potential energy landscapes are useful models in describing events such as protein folding and binding. While single-molecule fluorescence resonance energy transfer (smFRET) experiments encode information on continuous potentials for the system probed, including rarely visited barriers between putative potential minima, this information is rarely decoded from the data. This is because existing analysis methods often model smFRET output assuming, from the onset, that the system probed evolves in a discretized state space to be analyzed within a hidden Markov model (HMM) paradigm. By contrast, here, we infer continuous potentials from smFRET data without discretely approximating the state space. We do so by operating within a Bayesian nonparametric paradigm by placing priors on the family of all possible potential curves. As our inference accounts for a number of required experimental features raising computational cost (such as incorporating discrete photon shot noise), the framework leverages a structured-kernel-interpolation Gaussian process prior to help curtail computational cost. We show that our structured-kernel-interpolation priors for potential energy reconstruction from smFRET analysis accurately infers the potential energy landscape from a smFRET binding experiment. We then illustrate advantages of structured-kernel-interpolation priors for potential energy reconstruction from smFRET over standard HMM approaches by providing information, such as barrier heights and friction coefficients, that is otherwise inaccessible to HMMs.

The Effect of Time Resolution on Apparent Transition Path Times Observed in Single-Molecule Studies of Biomolecules

Single-molecule experiments have now achieved a time resolution allowing observation of transition paths, the brief trajectory segments where the molecule undergoing an unfolding or folding transition enters the energetically or entropically unfavorable barrier region from the folded/unfolded side and exits to the unfolded/folded side, thereby completing the transition. This resolution, however, is yet insufficient to identify the precise entrance/exit events that mark the beginning and the end of a transition path: the nature of the diffusive dynamics is such that a molecular trajectory will recross the boundary between the barrier region and the folded/unfolded state, multiple times, at a time scale much shorter than that of the typical experimental resolution. Here we use theory and Brownian dynamics simulations to show that, as a result of such recrossings, the apparent transition path times are generally longer than the true ones. We quantify this effect using a simple model where the observed dynamics is a moving average of the true dynamics and discuss experimental implications of our results.

On distributions of barrier crossing times as observed in single-molecule studies of biomolecules

Single-molecule experiments that monitor time evolution of molecular observables in real time have expanded beyond measuring transition rates toward measuring distributions of times of various molecular events. Of particular interest is the first-passage time for making a transition from one molecular configuration ( a ) to another ( b ) and conditional first-passage times such as the transition path time, which is the first-passage time from a to b conditional upon not leaving the transition region intervening between a and b . Another experimentally accessible (but not yet studied experimentally) observable is the conditional exit time, i.e., the time to leave the transition region through a specified boundary. The distributions of such times contain a wealth of mechanistic information about the transitions in question. Here, we use the first and the second (and, if desired, higher) moments of these distributions to characterize their relative width for the model in which the experimental observable undergoes Brownian motion in a potential of mean force. We show that although the distributions of transition path times are always narrower than exponential (in that the ratio of the standard deviation to the distribution's mean is always less than 1), distributions of first-passage times and of conditional exit times can be either narrow or broad, in some cases displaying long power-law tails. The conditional exit time studied here provides a generalization of the transition path time that also allows one to characterize the temporal scales of failed barrier crossing attempts.

The dynamics of subunit rotation in a eukaryotic ribosome

Protein synthesis by the ribosome is coordinated by an intricate series of large-scale conformational rearrangements. Structural studies can provide information about long-lived states, however biological kinetics are controlled by the intervening free-energy barriers. While there has been progress describing the energy landscapes of bacterial ribosomes, very little is known about the energetics of large-scale rearrangements in eukaryotic systems. To address this topic, we constructed an all-atom model with simplified energetics and performed simulations of subunit rotation in the yeast ribosome. In these simulations, the small subunit (SSU; ~1MDa) undergoes spontaneous and reversible rotations (~8°). By enabling the simulation of this rearrangement under equilibrium conditions, these calculations provide initial insights into the molecular factors that control dynamics in eukaryotic ribosomes. Through this, we are able to identify specific inter-subunit interactions that have a pronounced influence on the rate-limiting free-energy barrier. We also show that, as a result of changes in molecular flexibility, the thermodynamic balance between the rotated and unrotated states is temperature-dependent. This effect may be interpreted in terms of differential molecular flexibility within the rotated and unrotated states. Together, these calculations provide a foundation, upon which the field may begin to dissect the energetics of these complex molecular machines.

Hydrogen Bonds between Ions of Opposite and Like Charge in Hydroxyl-Functionalized Ionic Liquids: an Exhaustive Examination of the Interplay between Global and Local Motions and Intermolecular Hydrogen Bond Lifetimes and Kinetics

Hydroxyl-functionalized ionic liquids (ILs) represent a new interesting class of ILs where hydrogen bonds (HBs) play an important role: here, "typical" HBs between cations and anions (ca) are competing with "atypical" HBs connecting pairs of cations (cc). We study the equilibrium and kinetics of (cc) and (ca) HBs in 1-(n-hydroxyalkyl)-pyridinium bis(trifluoromethlysulfonyl)imide [HOC_nPy][NTf₂] ILs by means of molecular dynamics simulations. (cc) HBs are found to be between 0.96 and 3.76 kJ mol^-1 stronger than their (ca) counterparts, depending on the alkyl chain length. HB lifetimes and kinetics are analyzed by means of HB population and reactive flux correlation functions. Essentially, four different HB lifetimes have to be considered, spanning about 3 orders of magnitude, each valid in its own right and each associated with different aspects of HB breaking and HB reformation. The long-time limiting behavior of the HB population correlation function is controlled by diffusion of the ions and can be quantitatively described by analytical expressions. The short-time HB behavior is tied to the localized dynamics of the hydroxyl group exploring its local solvation environment. A minimalist kinetic two-domain model is introduced to realistically describe the time evolution of the HB population correlation function for both (ca) and (cc) HBs over 5 orders of magnitude. By employing the reactive flux method, we determine the kinetics of HB breaking, unaffected by diffusion processes. We determine both, the ultrafast upper boundary and the average rate of HB breaking, allowing recrossing-events during the transient relaxation time period. For sufficiently long alkyl chains, all those computed HB lifetimes indicate a higher kinetic stability of (cc) HBs over (ca) HBs; for short chains, it is vice-versa.