An official website of the United States government
Single-molecule fluorescence resonance energy transfer (FRET) experiments are commonly used to study the dynamics of molecular machines. While in vivo molecular processes often break time-reversal symmetry, the temporal directionality of cyclically operating molecular machines is often not evident from single-molecule FRET trajectories, especially in the most common two-color FRET studies. Solving a more quantitative problem of estimating the energy dissipation/entropy production by a molecular machine from single-molecule data is even more challenging. Here, we present a critical assessment of several practical methods of doing so, including Markov-model-based methods and a model-free approach based on an information-theoretical measure of entropy production that quantifies how (statistically) dissimilar observed photon sequences are from their time reverses. The Markov model approach is computationally feasible and may outperform model free approaches, but its performance strongly depends on how well the assumed model approximates the true microscopic dynamics. Markov models are also not guaranteed to give a lower bound on dissipation. Meanwhile, model-free, information-theoretical methods systematically underestimate entropy production at low photoemission rates, and long memory effects in the photon sequences make these methods demanding computationally. There is no clear winner among the approaches studied here, and all methods deserve to belong to a comprehensive data analysis toolkit.
more »
« less
Non-Markov models of single-molecule dynamics from information-theoretical analysis of trajectories
Whether single-molecule trajectories, observed experimentally or in molecular simulations, can be described using simple models such as biased diffusion is a subject of considerable debate. Memory effects and anomalous diffusion have been reported in a number of studies, but directly inferring such effects from trajectories, especially given limited temporal and/or spatial resolution, has been a challenge. Recently, we proposed that this can be achieved with information-theoretical analysis of trajectories, which is based on the general observation that non-Markov effects make trajectories more predictable and, thus, more “compressible” by lossless compression algorithms. Toy models where discrete molecular states evolve in time were shown to be amenable to such analysis, but its application to continuous trajectories presents a challenge: the trajectories need to be digitized first, and digitization itself introduces non-Markov effects that depend on the specifics of how trajectories are sampled. Here we develop a milestoning-based method for information-theoretical analysis of continuous trajectories and show its utility in application to Markov and non-Markov models and to trajectories obtained from molecular simulations.
more »
« less
- Award ID(s):
- 1955552
- PAR ID:
- 10499263
- Publisher / Repository:
- AIP
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 159
- Issue:
- 6
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Small integration time steps limit molecular dynamics (MD) simulations to millisecond time scales. Markov state models (MSMs) and equation-free approaches learn low-dimensional kinetic models from MD simulation data by performing configurational or dynamical coarse-graining of the state space. The learned kinetic models enable the efficient generation of dynamical trajectories over vastly longer time scales than are accessible by MD, but the discretization of configurational space and/or absence of a means to reconstruct molecular configurations precludes the generation of continuous all-atom molecular trajectories. We propose latent space simulators (LSS) to learn kinetic models for continuous all-atom simulation trajectories by training three deep learning networks to (i) learn the slow collective variables of the molecular system, (ii) propagate the system dynamics within this slow latent space, and (iii) generatively reconstruct molecular configurations. We demonstrate the approach in an application to Trp-cage miniprotein to produce novel ultra-long synthetic folding trajectories that accurately reproduce all-atom molecular structure, thermodynamics, and kinetics at six orders of magnitude lower cost than MD. The dramatically lower cost of trajectory generation enables greatly improved sampling and greatly reduced statistical uncertainties in estimated thermodynamic averages and kinetic rates.more » « less
-
Rotating disc electrode (RDE) voltammetry has been widely adopted for the study of heterogenized molecular electrocatalysts for multi-step fuel-forming reactions but this tool has never been comprehensively applied to their homogeneous analogues. Here, the utility and limitations of RDE techniques for mechanistic and kinetic analysis of homogeneous molecular catalysts that mediate multi-electron, multi-substrate redox transformations are explored. Using the ECEC′ reaction mechanism as a case study, two theoretical models are derived based on the Nernst diffusion layer model and the Hale transformation. Current–potential curves generated by these computational strategies are compared under a variety of limiting conditions to identify conditions under which the more minimalist Nernst Diffusion Layer approach can be applied. Based on this theoretical treatment, strategies for extracting kinetic information from the plateau current and the foot of the catalytic wave are derived. RDEV is applied to a cobaloxime hydrogen evolution reaction (HER) catalyst under non-aqueous conditions in order to experimentally validate this theoretical framework and explore the feasibility of RDE as a tool for studying homogeneous catalysts. Crucially, analysis of the foot-of-the-wave via this theoretical framework provides rate constants for elementary reaction steps that agree with those extracted from stationary voltammetric methods, supporting the application of RDE to study homogeneous fuel-forming catalysts. Finally, obstacles encountered during the kinetic analysis of cobaloxime, along with the voltammetric signatures used to diagnose this reactivity, are discussed with the goal of guiding groups working to improve RDE set-ups and help researchers avoid misinterpretation of RDE data.more » « less
-
Abstract Understanding phase transition between the liquid and gaseous states has gained significant interest, and has been ubiquitously observed in many places ranging from natural systems to water–energy nexus and thermal management applications. Phase transition phenomena at liquid–vapor interfaces are greatly governed by intermolecular‐level kinetics, which requires the use of empirical parameters in continuum‐level relations to explain the discrete nature of molecular particles. Despite its significance, it has been a great challenge to find detailed expressions of empirical parameters such as accommodation coefficients, which represent the probabilities for phase transition of liquid or vapor molecules at the interface. Here, direct statistical measurements of accommodation coefficients are reported by tracking the trajectories of liquid and vapor molecules in molecular simulations. The measurements reveal that evaporation and condensation coefficients are different by ≈50%, whereas they have been assumed to be equal in most previous studies. Then, the indirect measurement method is studied from a perspective of theoretical genetics based on the diffusion approximation. A good agreement between two approaches suggests that diffusion approximation can contribute to provide empirical parameters with a cost‐effective method.more » « less
-
An issue for molecular dynamics simulations is that events of interest often involve timescales that are much longer than the simulation time step, which is set by the fastest timescales of the model. Because of this timescale separation, direct simulation of many events is prohibitively computationally costly. This issue can be overcome by aggregating information from many relatively short simulations that sample segments of trajectories involving events of interest. This is the strategy of Markov state models (MSMs) and related approaches, but such methods suffer from approximation error because the variables defining the states generally do not capture the dynamics fully. By contrast, once converged, the weighted ensemble (WE) method aggregates information from trajectory segments so as to yield unbiased estimates of both thermodynamic and kinetic statistics. Unfortunately, errors decay no faster than unbiased simulation in WE as originally formulated and commonly deployed. Here, we introduce a theoretical framework for describing WE that shows that the introduction of an approximate stationary distribution on top of the stratification, as in nonequilibrium umbrella sampling (NEUS), accelerates convergence. Building on ideas from MSMs and related methods, we generalize the NEUS approach in such a way that the approximation error can be reduced systematically. We show that the improved algorithm can decrease the simulation time required to achieve the desired precision by orders of magnitude.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1063/5.0158930
