Computing the equilibrium properties of complex systems, such as free energy differences, is often hampered by rare events in the dynamics. Enhanced sampling methods may be used in order to speed up sampling by, for example, using high temperatures,
as in parallel tempering, or simulating with a biasing potential such as in the case of umbrella sampling. The equilibrium properties of the thermodynamic state of interest (e.g., lowest temperature or unbiased potential) can be computed using reweighting estimators such as the weighted histogram analysis method or the multistate Bennett acceptance ratio (MBAR). weighted histogram analysis method and MBAR produce unbiased estimates, the simulation samples from the global equilibria at their respective thermodynamic state--a requirement that can be prohibitively expensive for some simulations such as a large parallel tempering ensemble of an explicitly solvated biomolecule. Here, we introduce the transition-based reweighting analysis method (TRAM)--a class of estimators that exploit ideas from Markov modeling and only require the simulation data to be in local equilibrium within subsets of the configuration space. We formulate the expanded TRAM (xTRAM) estimator that is shown to be asymptotically unbiased and a generalization of MBAR. Using four exemplary systems of varying complexity, we demonstrate the improved convergence (ranging from a twofold improvement to several orders of magnitude) of xTRAM in comparison to a direct counting estimator and MBAR, with respect to the invested simulation effort. Lastly, we introduce a random-swapping simulation protocol that can be used with xTRAM, gaining orders-of-magnitude advantages over simulation protocols that require the constraint of sampling from a global equilibrium.
We explore the dynamical large-deviations of a lattice heteropolymer model of a protein by means of path sampling of trajectories. We uncover the existence of non-equilibrium dynamical phase-transitions in ensembles of trajectories between active and
inactive dynamical phases, whose nature depends on properties of the interaction potential. When the full heterogeneity of interactions due to the amino-acid sequence is preserved, as in a fully interacting model or in a heterogeneous version of the G={o} model where only native interactions are considered, the transition is between the equilibrium native state and a highly native but kinetically trapped state. In contrast, for the homogeneous G={o} model, where there is a single native energy and the sequence plays no role, the dynamical transition is a direct consequence of the static bi-stability between unfolded and native states. In the heterogeneous case the native-active and native-inactive states, despite their static similarity, have widely varying dynamical properties, and the transition between them occurs even in lattice proteins whose sequences are designed to make them optimal folders.
