Identity of distributions of direct uphill and downhill translocation times for particles traversing membrane channels

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.

Counter-Intuitive Features of Particle Dynamics in Nanopores

Using the framework of a continuous diffusion model based on the Smoluchowski equation, we analyze particle dynamics in the confinement of a transmembrane nanopore. We briefly review existing analytical results to highlight consequences of interactions between the channel nanopore and the translocating particles. These interactions are described within a minimalistic approach by lumping together multiple physical forces acting on the particle in the pore into a one-dimensional potential of mean force. Such radical simplification allows us to obtain transparent analytical results, often in a simple algebraic form. While most of our findings are quite intuitive, some of them may seem unexpected and even surprising at first glance. The focus is on five examples: (i) attractive interactions between the particles and the nanopore create a potential well and thus cause the particles to spend more time in the pore but, nevertheless, increase their net flux; (ii) if the potential well-describing particle-pore interaction occupies only a part of the pore length, the mean translocation time is a non-monotonic function of the well length, first increasing and then decreasing with the length; (iii) when a rectangular potential well occupies the entire nanopore, the mean particle residence time in the pore is independent of the particle diffusivity inside the pore and depends only on its diffusivity in the bulk; (iv) although in the presence of a potential bias applied to the nanopore the "downhill" particle flux is higher than the "uphill" one, the mean translocation times and their distributions are identical, i.e., independent of the translocation direction; and (v) fast spontaneous gating affects nanopore selectivity when its characteristic time is comparable to that of the particle transport through the pore.

Time-Resolved Statistics of Snippets as General Framework for Model-Free Entropy Estimators

Irreversibility is commonly quantified by entropy production. An external observer can estimate it through measuring an observable that is antisymmetric under time reversal like a current. We introduce a general framework that allows us to infer a lower bound on entropy production through measuring the time-resolved statistics of events with any symmetry under time reversal, in particular, time-symmetric instantaneous events. We emphasize Markovianity as a property of certain events rather than of the full system and introduce an operationally accessible criterion for this weakened Markov property. Conceptually, the approach is based on snippets as particular sections of trajectories between two Markovian events, for which a generalized detailed balance relation is discussed.

Waiting Time Distributions in Hybrid Models of Motor-Bead Assays: A Concept and Tool for Inference

In single-molecule experiments, the dynamics of molecular motors are often observed indirectly by measuring the trajectory of an attached bead in a motor-bead assay. In this work, we propose a method to extract the step size and stalling force for a molecular motor without relying on external control parameters. We discuss this method for a generic hybrid model that describes bead and motor via continuous and discrete degrees of freedom, respectively. Our deductions are solely based on the observation of waiting times and transition statistics of the observable bead trajectory. Thus, the method is non-invasive, operationally accessible in experiments and can, in principle, be applied to any model describing the dynamics of molecular motors. We briefly discuss the relation of our results to recent advances in stochastic thermodynamics on inference from observable transitions. Our results are confirmed by extensive numerical simulations for parameters values of an experimentally realized F1-ATPase assay.

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.

First-passage times in conical varying-width channels biased by a transverse gravitational force: Comparison of analytical and numerical results

We study the crossing time statistic of diffusing point particles between the two ends of expanding and narrowing two-dimensional conical channels under a transverse external gravitational field. The theoretical expression for the mean first-passage time for such a system is derived under the assumption that the axial diffusion in a two-dimensional channel of smoothly varying geometry can be approximately described as a one-dimensional diffusion in an entropic potential with position-dependent effective diffusivity in terms of the modified Fick-Jacobs equation. We analyze the channel crossing dynamics in terms of the mean first-passage time, combining our analytical results with extensive two-dimensional Brownian dynamics simulations, allowing us to find the range of applicability of the one-dimensional approximation. We find that the effective particle diffusivity decreases with increasing amplitude of the external potential. Remarkably, the mean first-passage time for crossing the channel is shown to assume a minimum at finite values of the potential amplitude.

Inferring entropy production rate from partially observed Langevin dynamics under coarse-graining

The entropy production rate (EPR) measures time-irreversibility in systems operating far from equilibrium. The challenge in estimating the EPR for a continuous variable system is the finite spatiotemporal resolution and the limited accessibility to all of the nonequilibrium degrees of freedom. Here, we estimate the irreversibility in partially observed systems following oscillatory dynamics governed by coupled overdamped Langevin equations. We coarse-grain an observed variable of a nonequilibrium driven system into a few discrete states and estimate a lower bound on the total EPR. As a model system, we use hair-cell bundle oscillations driven by molecular motors, such that the bundle tip position is observed, but the positions of the motors are hidden. In the observed variable space, the underlying driven process exhibits second-order semi-Markov statistics. The waiting time distributions (WTD), associated with transitions among the coarse-grained states, are non-exponential and convey the information on the broken time-reversal symmetry. By invoking the underlying time-irreversibility, we calculate a lower bound on the total EPR from the Kullback-Leibler divergence (KLD) between WTD. We show that the mean dwell-time asymmetry factor - the ratio between the mean dwell-times along the forward direction and the backward direction, can qualitatively measure the degree of broken time reversal symmetry and increases with finer spatial resolution. Finally, we apply our methodology to a continuous-time discrete Markov chain model, coarse-grained into a linear system exhibiting second-order semi-Markovian statistics, and demonstrate the estimation of a lower bound on the total EPR from irreversibility manifested only in the WTD.

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.

Crowding breaks the forward/backward symmetry of transition times in biased random walks

Microscopic mechanisms of natural processes are frequently understood in terms of random walk models by analyzing local particle transitions. This is because these models properly account for dynamic processes at the molecular level and provide a clear physical picture. Recent theoretical studies made a surprising discovery that in complex systems, the symmetry of molecular forward/backward transition times with respect to local bias in the dynamics may be broken and it may take longer to go downhill than uphill. The physical origins of these phenomena remain not fully understood. Here, we explore in more detail the microscopic features of the symmetry breaking in the forward/backward transition times by analyzing exactly solvable discrete-state stochastic models. In particular, we consider a specific case of two random walkers on a four-site periodic lattice as the way to represent the general systems with multiple pathways. It is found that the asymmetry in transition times depends on several factors that include the degree of deviation from equilibrium, the particle crowding, and methods of measurements of dynamic properties. Our theoretical analysis suggests that the asymmetry in transition times can be explored experimentally for determining the important microscopic features of natural processes by quantitatively measuring the local deviations from equilibrium and the degrees of crowding.

Effect of Memory and Inertial Contribution on Transition-Time Distributions: Theory and Simulations

Transition paths refer to the time taken by molecules to cross a barrier separating two molecular conformations. In this work, we study how memory, as well as inertial contribution in the dynamics along a reaction coordinate, can affect the distribution of the transition-path time. We use a simple model of dynamics governed by a generalized Langevin equation with a power-law memory along with the inertial term, which was neglected in previous studies, where memory effects were explored only in the overdamped limit. We derive an approximate expression for the transit-time distribution and discuss our results for the short- and long-time limits and also compare it with known results in the high friction (overdamped) limit as well as in the Markovian limit. We have developed a numerical algorithm to test our theoretical results against extensive numerical simulations.

Barrier Crossing Dynamics from Single-Molecule Measurements

Chemists visualize chemical reactions as motion along one-dimensional "reaction coordinates" over free energy barriers. Various rate theories, such as transition state theory and the Kramers theory of diffusive barrier crossing, differ in their assumptions regarding the mathematical specifics of this motion. Direct experimental observation of the motion along reaction coordinates requires single-molecule experiments performed with unprecedented time resolution. Toward this goal, recent single-molecule studies achieved time resolution sufficient to catch biomolecules in the act of crossing free energy barriers as they fold, bind to their targets, or undergo other large structural changes, offering a window into the elusive reaction "mechanisms". This Perspective describes what we can learn (and what we have already learned) about barrier crossing dynamics through synergy of single-molecule experiments, theory, and molecular simulations. In particular, I will discuss how emerging experimental data can be used to answer several questions of principle. For example, is motion along the reaction coordinate diffusive, is there conformational memory, and is reduction to just one degree of freedom to represent the reaction mechanism justified? It turns out that these questions can be formulated as experimentally testable mathematical inequalities, and their application to experimental and simulated data has already led to a number of insights. I will also discuss open issues and current challenges in this fast evolving field of research.

Broad distributions of transition-path times are fingerprints of multidimensionality of the underlying free energy landscapes

Recent single-molecule experiments have observed transition paths, i.e., brief events where molecules (particularly biomolecules) are caught in the act of surmounting activation barriers. Such measurements offer unprecedented mechanistic insights into the dynamics of biomolecular folding and binding, molecular machines, and biological membrane channels. A key challenge to these studies is to infer the complex details of the multidimensional energy landscape traversed by the transition paths from inherently low-dimensional experimental signals. A common minimalist model attempting to do so is that of one-dimensional diffusion along a reaction coordinate, yet its validity has been called into question. Here, we show that the distribution of the transition path time, which is a common experimental observable, can be used to differentiate between the dynamics described by models of one-dimensional diffusion from the dynamics in which multidimensionality is essential. Specifically, we prove that the coefficient of variation obtained from this distribution cannot possibly exceed 1 for any one-dimensional diffusive model, no matter how rugged its underlying free energy landscape is: In other words, this distribution cannot be broader than the single-exponential one. Thus, a coefficient of variation exceeding 1 is a fingerprint of multidimensional dynamics. Analysis of transition paths in atomistic simulations of proteins shows that this coefficient often exceeds 1, signifying essential multidimensionality of those systems.

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.

Detailed balance for diffusion in a potential with trapping and forward-backward symmetry of trapping time distributions

For particles diffusing in a potential, detailed balance guarantees the absence of net fluxes at equilibrium. Here, we show that the conventional detailed balance condition is a special case of a more general relation that works when the diffusion occurs in the presence of a distributed sink that eventually traps the particle. We use this relation to study the lifetime distribution of particles that start and are trapped at specified initial and final points. It turns out that when the sink strength at the initial point is nonzero, the initial and final points are interchangeable, i.e., the distribution is independent of which of the two points is initial and which is final. In other words, this conditional trapping time distribution possesses forward–backward symmetry.

Biased Random Walk in Crowded Environment: Breaking Uphill/Downhill Symmetry of Transition Times

Various natural processes can be analyzed using the concept of random walks. For a single random walker, the mean waiting times for uphill and downhill transitions between neighboring sites are equal. Here we investigate the uphill/downhill symmetry of waiting times for transitions of a tracer in crowded environment using exactly solvable one-dimensional stochastic models. It is found that, unexpectedly, the time to move in the direction of the bias (downhill) is always longer than the time to move against the bias (uphill). The degree of asymmetry depends on the particle density, the strength of the bias, and the size of the system. The microscopic origin of the symmetry breaking is discussed.

Transition path dynamics in the binding of intrinsically disordered proteins: A simulation study

Association of proteins and other biopolymers is a ubiquitous process in living systems. Recent single-molecule measurements probe the dynamics of association in unprecedented detail by measuring the properties of association transition paths, i.e., short segments of molecular trajectories between the time the proteins are close enough to interact and the formation of the final complex. Interpretation of such measurements requires adequate models for describing the dynamics of experimental observables. In an effort to develop such models, here we report a simulation study of the association dynamics of two oppositely charged, disordered polymers. We mimic experimental measurements by monitoring intermonomer distances, which we treat as "experimental reaction coordinates." While the dynamics of the distance between the centers of mass of the molecules is found to be memoryless and diffusive, the dynamics of the experimental reaction coordinates displays significant memory and can be described by a generalized Langevin equation with a memory kernel. We compute the most commonly measured property of transition paths, the distribution of the transition path time, and show that, despite the non-Markovianity of the underlying dynamics, it is well approximated as one-dimensional diffusion in the potential of mean force provided that an apparent value of the diffusion coefficient is used. This apparent value is intermediate between the slow (low frequency) and fast (high frequency) limits of the memory kernel. We have further studied how the mean transition path time depends on the ionic strength and found only weak dependence despite strong electrostatic attraction between the polymers.

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.

Mechanism of the electroneutral sodium/proton antiporter PaNhaP from transition-path shooting

Na⁺/H⁺ antiporters exchange sodium ions and protons on opposite sides of lipid membranes. The electroneutral Na⁺/H⁺ antiporter NhaP from archaea Pyrococcus abyssi (PaNhaP) is a functional homolog of the human Na⁺/H⁺ exchanger NHE1, which is an important drug target. Here we resolve the Na⁺ and H⁺ transport cycle of PaNhaP by transition-path sampling. The resulting molecular dynamics trajectories of repeated ion transport events proceed without bias force, and overcome the enormous time-scale gap between seconds-scale ion exchange and microseconds simulations. The simulations reveal a hydrophobic gate to the extracellular side that opens and closes in response to the transporter domain motion. Weakening the gate by mutagenesis makes the transporter faster, suggesting that the gate balances competing demands of fidelity and efficiency. Transition-path sampling and a committor-based reaction coordinate optimization identify the essential motions and interactions that realize conformational alternation between the two access states in transporter function.

Cation-proton antiporters mediate selective ion exchange across cellular membranes to control pH, salt concentration and cell volume. Here the authors present a transition-path sampling method that overcomes the timescale gap between simulations (µs) and transport processes (s), which allows them to resolve the Na⁺ and H⁺ transport cycle of the Na⁺/H⁺ antiporter NhaP from Pyrococcus abyssi.

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.

Experimental evidence of symmetry breaking of transition-path times

While thermal rates of state transitions in classical systems have been studied for almost a century, associated transition-path times have only recently received attention. Uphill and downhill transition paths between states at different free energies should be statistically indistinguishable. Here, we systematically investigate transition-path-time symmetry and report evidence of its breakdown on the molecular- and meso-scale out of equilibrium. In automated Brownian dynamics experiments, we establish first-passage-time symmetries of colloids driven by femtoNewton forces in holographically-created optical landscapes confined within microchannels. Conversely, we show that transitions which couple in a path-dependent manner to fluctuating forces exhibit asymmetry. We reproduce this asymmetry in folding transitions of DNA-hairpins driven out of equilibrium and suggest a topological mechanism of symmetry breakdown. Our results are relevant to measurements that capture a single coordinate in a multidimensional free energy landscape, as encountered in electrophysiology and single-molecule fluorescence experiments.

Oscillations in the mean transition time of a particle scattered on a double slit potential

Scattering through a double slit potential is one of the most fundamental problems in quantum mechanics. It is well understood that due to the superposition of amplitudes, one observes a spatial interference pattern in the scattered wavefunction reflecting the superposition of amplitudes coming from both slits. However, the effect of the double slit on the mean time it takes to traverse the slit has not been considered previously. Using a transition path time formalism, we show that when a single Gaussian wavepacket is scattered through a double slit potential, one finds not only oscillations in the scattered density resulting from the spatial interference created by the splitting of the wavepacket but also an oscillatory pattern in the mean scattering time. Long times are associated with low values of a suitably defined momentum, and short times with higher values. The double slit thus serves as a momentum filtering device. We also find an interference pattern in the time averaged momentum weak value profile of the scattered particle implying that the double slit also acts as a weak momentum filter. These results not only demonstrate the value of considering transition path time distributions in their quantum mechanical context but also present a challenge to semiclassical approximations-can they account for temporal interference?

Transition path times reveal memory effects and anomalous diffusion in the dynamics of protein folding

Recent single-molecule experiments probed transition paths of biomolecular folding and, in particular, measured the time biomolecules spend while crossing their free energy barriers. A surprising finding from these studies is that the transition barriers crossed by transition paths, as inferred from experimentally observed transition path times, are often lower than the independently determined free energy barriers. Here we explore memory effects leading to anomalous diffusion as a possible origin of this discrepancy. Our analysis of several molecular dynamics trajectories shows that the dynamics of common reaction coordinates used to describe protein folding is subdiffusive, at least at sufficiently short times. We capture this effect using a one-dimensional fractional Brownian motion (FBM) model, in which the system undergoes a subdiffusive process in the presence of a potential of mean force, and show that this model yields much broader distributions of transition path times with stretched exponential long-time tails. Without any adjustable parameters, these distributions agree well with the transition path times computed directly from protein trajectories. We further discuss how the FBM model can be tested experimentally.

Communication: Coordinate-dependent diffusivity from single molecule trajectories

Single-molecule observations of biomolecular folding are commonly interpreted using the model of one-dimensional diffusion along a reaction coordinate, with a coordinate-independent diffusion coefficient. Recent analysis, however, suggests that more general models are required to account for single-molecule measurements performed with high temporal resolution. Here, we consider one such generalization: a model where the diffusion coefficient can be an arbitrary function of the reaction coordinate. Assuming Brownian dynamics along this coordinate, we derive an exact expression for the coordinate-dependent diffusivity in terms of the splitting probability within an arbitrarily chosen interval and the mean transition path time between the interval boundaries. This formula can be used to estimate the effective diffusion coefficient along a reaction coordinate directly from single-molecule trajectories.

Broken detailed balance and non-equilibrium dynamics in living systems: a review

Living systems operate far from thermodynamic equilibrium. Enzymatic activity can induce broken detailed balance at the molecular scale. This molecular scale breaking of detailed balance is crucial to achieve biological functions such as high-fidelity transcription and translation, sensing, adaptation, biochemical patterning, and force generation. While biological systems such as motor enzymes violate detailed balance at the molecular scale, it remains unclear how non-equilibrium dynamics manifests at the mesoscale in systems that are driven through the collective activity of many motors. Indeed, in several cellular systems the presence of non-equilibrium dynamics is not always evident at large scales. For example, in the cytoskeleton or in chromosomes one can observe stationary stochastic processes that appear at first glance thermally driven. This raises the question how non-equilibrium fluctuations can be discerned from thermal noise. We discuss approaches that have recently been developed to address this question, including methods based on measuring the extent to which the system violates the fluctuation-dissipation theorem. We also review applications of this approach to reconstituted cytoskeletal networks, the cytoplasm of living cells, and cell membranes. Furthermore, we discuss a more recent approach to detect actively driven dynamics, which is based on inferring broken detailed balance. This constitutes a non-invasive method that uses time-lapse microscopy data, and can be applied to a broad range of systems in cells and tissue. We discuss the ideas underlying this method and its application to several examples including flagella, primary cilia, and cytoskeletal networks. Finally, we briefly discuss recent developments in stochastic thermodynamics and non-equilibrium statistical mechanics, which offer new perspectives to understand the physics of living systems.

Instantaneous Tunneling Flight Time for Wavepacket Transmission through Asymmetric Barriers

The time it takes a particle to tunnel through the asymmetric Eckart barrier potential is investigated using Gaussian wavepackets, where the barrier serves as a model for the potential along a chemical reaction coordinate. We have previously shown that the, in principle experimentally measurable, tunneling flight time, which determines the time taken by the transmitted particle to traverse the barrier, vanishes for symmetric potentials like the Eckart and square barrier [ Petersen , J. ; Pollak , E. J. Phys. Chem. Lett. 2017 , 9 , 4017 ]. Here we show that the same result is obtained for the asymmetric Eckart barrier potential, and therefore, the zero tunneling flight time seems to be a general result for one-dimensional time-independent potentials. The wavepacket dynamics is simulated using both an exact quantum mechanical method and a classical Wigner prescription. The excellent agreement between the two methods shows that quantum coherences are not important in pure one-dimensional tunneling and reinforces the conclusion that the tunneling flight time vanishes.

First passage, looping, and direct transition in expanding and narrowing tubes: Effects of the entropy potential

We study transitions of diffusing particles between the left and right ends of expanding and narrowing conical tubes. In an expanding tube, such transitions occur faster than in the narrowing tube of the same length and radius variation rate. This happens because the entropy potential pushes the particle towards the wide tube end, thus accelerating the transitions in the expanding tube and slowing them down in the narrowing tube. To gain deeper insight into how the transitions occur, we divide each trajectory into the direct-transit and looping segments. The former is the final part of the trajectory, where the particle starting from the left tube end goes to the right end without returning to the left one. The rest of the trajectory is the looping segment, where the particle, starting from the left tube end, returns to this end again and again until the direct transition happens. Our focus is on the durations of the two segments and their sum, which is the duration of the particle first passage between the left and right ends of the tube. We approach the problem using the one-dimensional description of the particle diffusion along the tube axis in terms of the modified Fick-Jacobs equation. This allows us to derive analytical expressions for the Laplace transforms of the probability densities of the first-passage, direct-transit, and looping times, which we use to find the mean values of these random variables. Our results show that the direct transits are independent of the entropy potential and occur as in free diffusion. However, this "free diffusion" occurs with the effective diffusivity entering the modified Fick-Jacobs equation, which is smaller than the particle diffusivity in a cylindrical tube. This is the only way how the varying tube geometry manifests itself in the direct transits. Since direct-transit times are direction-independent, the difference in the first-passage times in the tubes of the two types is due to the difference in the durations of the looping segments in the expanding and narrowing tubes. Obtained analytical results are supported by three-dimensional Brownian dynamics simulations.

A new insight into diffusional escape from a biased cylindrical trap

Recent experiments with single biological nanopores, as well as single-molecule fluorescence spectroscopy and pulling studies of protein and nucleic acid folding raised a number of questions that stimulated theoretical and computational investigations of barrier crossing dynamics. The present paper addresses a closely related problem focusing on trajectories of Brownian particles that escape from a cylindrical trap in the presence of a force F parallel to the cylinder axis. To gain new insights into the escape dynamics, we analyze the "fine structure" of these trajectories. Specifically, we divide trajectories into two segments: a looping segment, when a particle unsuccessfully tries to escape returning to the trap bottom again and again, and a direct-transit segment, when it finally escapes moving without touching the bottom. Analytical expressions are derived for the Laplace transforms of the probability densities of the durations of the two segments. These expressions are used to find the mean looping and direct-transit times as functions of the biasing force F. It turns out that the force-dependences of the two mean times are qualitatively different. The mean looping time monotonically increases as F decreases, approaching exponential F-dependence at large negative forces pushing the particle towards the trap bottom. In contrast to this intuitively appealing behavior, the mean direct-transit time shows rather counterintuitive behavior: it decreases as the force magnitude, |F|, increases independently of whether the force pushes the particles to the trap bottom or to the exit from the trap, having a maximum at F = 0.

Transition Path Time Distribution, Tunneling Times, Friction, and Uncertainty

A quantum mechanical transition path time probability distribution is formulated and its properties are studied using a parabolic barrier potential model. The average transit time is well defined and readily calculated. It is smaller than the analogous classical mechanical average transit time, vanishing at the crossover temperature. It provides a direct route for determining tunneling times. The average time may be also used to define a coarse grained momentum of the system for the passage from one side of the barrier to the other. The product of the uncertainty in this coarse grained momentum with the uncertainty in the location of the particle is shown under certain conditions to be smaller than the ℏ/2 formal uncertainty limit. The model is generalized to include friction in the form of a bilinear interaction with a harmonic bath. Using an Ohmic friction model one finds that increasing the friction, increases the transition time. Only moderate values of the reduced friction coefficient are needed for the quantum transition time and coarse grained uncertainty to approach the classical limit which is smaller than ℏ/2 when the friction is not too small. These results show how one obtains classical dynamics from a pure quantum system without invoking any further assumptions, approximations, or postulates.

Transition path time distribution and the transition path free energy barrier

Free energy profile, showing why the transition path barrier is lower than the free energy of activation.

On the applicability of entropy potentials in transport problems

Transport in confined structures of varying geometry has become the subject of growing attention in recent years since such structures are ubiquitous in biology and technology. In analyzing transport in systems of this type, the notion of entropy potentials is widely used. Entropy potentials naturally arise in one-dimensional description of equilibrium distributions in multidimensional confined structures. However, their application to transport problems requires some caution. In this article we discuss such applications and summarize the results of recent studies exploring the limits of applicability. We also consider an example of a transport problem in a system of varying geometry, where the conventional approach is inapplicable since the geometry changes abruptly. In addition, we demonstrate how the entropy potential can be used to analyze optimal transport through a tree-dimensional cosine-shaped channel.

Computation of transit times using the milestoning method with applications to polymer translocation

Milestoning is an efficient approximation for computing long-time kinetics and thermodynamics of large molecular systems, which are inaccessible to brute-force molecular dynamics simulations. A common use of milestoning is to compute the mean first passage time (MFPT) for a conformational transition of interest. However, the MFPT is not always the experimentally observed timescale. In particular, the duration of the transition path, or the mean transit time, can be measured in single-molecule experiments, such as studies of polymers translocating through pores and fluorescence resonance energy transfer studies of protein folding. Here we show how to use milestoning to compute transit times and illustrate our approach by applying it to the translocation of a polymer through a narrow pore.

Parallel sub-micrometre channels with different dimensions for laser scattering detection

A novel and simple approach for the realization of polymer sub-micrometre channels is introduced by exploiting replica molding of Pt wires deposited by focused ion beam. We fabricate arrays of parallel channels with typical dimensions down to 600 nm and with variable height. We characterize the pressure-driven transport of polymer colloids through the channels in terms of the translocation frequency, amplitude and duration by implementing a laser scattering detection technique. We propose a prototype application of the presented platform such as the in situ sizing and sensing of populations of particles with different dimensions down to 50 nm.

Sequence dependence of the binding energy in chaperone-driven polymer translocation through a nanopore

We study the translocation of stiff polymers through a nanopore, driven by the chemical-potential gradient exerted by binding proteins (chaperones) on the trans side of the pore. Bound chaperones prevent backsliding through the pore and, therefore, partially rectify the polymer passage. We show that the sequence of chain monomers with different binding affinity for the chaperones significantly affects the translocation dynamics. In particular, we investigate the effect of the nearest-neighbor adjacency probability of the two monomer types. Depending on the magnitude of the involved binding energies, the translocation speed may either increase or decrease with the adjacency probability. We determine the mean first passage time and show that, by tuning the effective binding energy, the motion changes continuously from purely diffusive to ballistic translocation.

Facilitated translocation of polypeptides through a single nanopore

The transport of polypeptides through nanopores is a key process in biology and medical biotechnology. Despite its critical importance, the underlying kinetics of polypeptide translocation through protein nanopores is not yet comprehensively understood. Here, we present a simple two-barrier, one-well kinetic model for the translocation of short positively charged polypeptides through a single transmembrane protein nanopore that is equipped with negatively charged rings, simply called traps. We demonstrate that the presence of these traps within the interior of the nanopore dramatically alters the free energy landscape for the partitioning of the polypeptide into the nanopore interior, as revealed by significant modifications in the activation free energies required for the transitions of the polypeptide from one state to the other. Our kinetic model permits the calculation of the relative and absolute exit frequencies of the short cationic polypeptides through either opening of the nanopore. Moreover, this approach enabled quantitative assessment of the kinetics of translocation of the polypeptides through a protein nanopore, which is strongly dependent on several factors, including the nature of the translocating polypeptide, the position of the traps, the strength of the polypeptide-attractive trap interactions and the applied transmembrane voltage.

Bridging timescales and length scales: from macroscopic flux to the molecular mechanism of antibiotic diffusion through porins

Our aim in this study was to provide an atomic description of ampicillin translocation through OmpF, the major outer membrane channel in Escherichia coli and main entry point for beta-lactam antibiotics. By applying metadynamics simulations, we also obtained the energy barriers along the diffusion pathway. We then studied the effect of mutations that affect the charge and size at the channel constriction zone, and found that in comparison to the wild-type, much lower energy barriers are required for translocation. The expected higher translocation rates were confirmed on the macroscopic scale by liposome-swelling assays. A microscopic view on the millisecond timescale was obtained by analysis of temperature-dependent ion current fluctuations in the presence of ampicillin and provide the enthalpic part of the energy barrier. By studying antibiotic translocation over various timescales and length scales, we were able to discern its molecular mechanism and rate-limiting interactions, and draw biologically relevant conclusions that may help in the design of drugs with enhanced permeation rates.

Transition times in the low-noise limit of stochastic dynamics

We study the transition time distribution for a particle moving between two wells of a multidimensional potential in the low-noise limit of overdamped Langevin dynamics. Possible transition paths are restricted to a thin tube surrounding the most probable trajectory. We demonstrate that finding the transition time distribution reduces to a one-dimensional problem. The resulting transition time distribution has a universal and compact form. We suggest that transition barriers can be estimated from a single-temperature experiment if both the life times and the transition times are measured.

Drift and diffusion in periodic potentials: upstream and downstream step times are distributed identically

This note deals with particles diffusing in one dimension in the presence of a periodic potential and a uniform driving force. We show that (i) the probabilities for the particle to make a step of the length L, where L is the period, in the upstream and downstream directions are independent of the periodic potential, and (ii) the distributions of the step time are independent of the step direction. These two characteristics are used to derive expressions for the effective drift velocity and diffusion coefficient.

Polyelectrolyte entry and transport through an asymmetric alpha-hemolysin channel

We study the entry and transport of a polyelectrolyte, dextran sulfate (DS), through an asymmetric alpha-hemolysin protein channel inserted into a planar lipid bilayer. We compare the dynamics of the DS chains as they enter the channel at the opposite stem or vestibule sides. Experiments are performed at the single-molecule level by using an electrical method. The frequency of current blockades varies exponentially as a function of applied voltage. This frequency is smaller for the stem entrance than for the vestibule one, due to a smaller coupling with the electric field and a larger activation energy for entry. The value of the activation energy is quantitatively interpreted as an entropic effect of chain confinement. The translocation time decreases when the applied voltage increases and displays an exponential variation which is independent of the stem or vestibule sides.

Weak disorder strongly improves the selective enhancement of diffusion in a tilted periodic potential

The diffusion of an overdamped Brownian particle in a tilted periodic potential is known to exhibit a pronounced enhancement over the free thermal diffusion within a small interval of tilt values. Here we show that weak disorder in the form of small, time-independent deviations from a strictly spatially periodic potential may further boost this diffusion peak by orders of magnitude. Our general theoretical predictions are in excellent agreement with experimental observations.

Fluctuation theorem for channel-facilitated membrane transport of interacting and noninteracting solutes

In this paper, we discuss the fluctuation theorem for channel-facilitated transport of solutes through a membrane separating two reservoirs. The transport is characterized by the probability, P(n)(t), that n solute particles have been transported from one reservoir to the other in time t. The fluctuation theorem establishes a relation between P(n)(t) and P-(n)(t): The ratio P(n)(t)/P-(n)(t) is independent of time and equal to exp(nbetaA), where betaA is the affinity measured in the thermal energy units. We show that the same fluctuation theorem is true for both single- and multichannel transport of noninteracting particles and particles which strongly repel each other.

Counting translocations of strongly repelling particles through single channels: fluctuation theorem for membrane transport

Transport of strongly repelling particles through a single membrane channel is analyzed assuming that the channel cannot be occupied by more than one particle. An exact solution is found for the Laplace transform of the probability P_{n}(t) that n particles have been transported in time t. This transform is used to find the flux through the channel and to show that P_{n}(t) and P_{-n}(t) are related by the fluctuation theorem. The solution is obtained using an observation that P_{n}(t) is the propagator for a non-Markovian random walk, which can be found by solving a set of integral equations.