DIME: Dynamics Inferred from Monte Carlo Ensembles via Continuous-Time Markov Chains
TL;DR
A framework for constructing a continuous-time Markov chain (CTMC) from a static, Boltzmann-weighted ensemble is presented, enabling the generation of physically plausible kinetic trajectories without requiring extensive molecular dynamics simulation.
Abstract
Monte Carlo (MC) sampling and related enhanced sampling methods can efficiently explore the conformational space of complex biomolecular systems and accurately estimate equilibrium populations of conformational states. However, the absence of temporal information prevents these methods from directly describing molecular kinetics. Here we present a framework for constructing a continuous-time Markov chain (CTMC) from a static, Boltzmann-weighted ensemble, enabling the generation of physically plausible kinetic trajectories without requiring extensive molecular dynamics (MD) simulation. The method includes (i) clustering the MC ensemble into metastable discrete states whose populations define the stationary distribution, (ii) building a rate matrix over these states that satisfies detailed balance by design, with transition rates set by the free energy barriers separating the states, estimated directly from the free energy surface of the ensemble, and (iii) calibrating the absolute time scale by matching the slowest relaxation mode to a short MD reference trajectory. We validate the approach across six molecular systems spanning a broad range of complexity: n-pentane and alanine dipeptide as low-dimensional benchmarks; a WLALL pentapeptide; the achiral peptide AIB9, whose left- and right-handed helices are related by an exact symmetry that provides a reference-free test of the inferred kinetics; wild-type chignolin as a β-hairpin folder; and the millisecond native-state dynamics of bovine pancreatic trypsin inhibitor (BPTI) from a long reference trajectory. In each case the inferred model reproduces MD-derived rate matrices, mean first-passage times, and implied time scales following a short calibration step. Since the method relies on the free energy surface, it can equally well be applied to surfaces obtained from enhanced sampling simulations, as demonstrated for alanine dipeptide using well-tempered metadynamics. The method provides a principled route to kinetic insight when MD convergence is prohibitively expensive but thermodynamic sampling is tractable and serves as a complement to existing Markov state model approaches that build kinetic models from MD data.