Statistically optimal analysis of state-discretized trajectory data from multiple thermodynamic states

Hao Wu, Antonia S J S Mey, Edina Rosta, Frank Noé

Research output: Contribution to journalArticlepeer-review

Abstract

We propose a discrete transition-based reweighting analysis method (dTRAM) for analyzing configuration-space-discretized simulation trajectories produced at different thermodynamic states (temperatures, Hamiltonians, etc.) dTRAM provides maximum-likelihood estimates of stationary quantities (probabilities, free energies, expectation values) at any thermodynamic state. In contrast to the weighted histogram analysis method (WHAM), dTRAM does not require data to be sampled from global equilibrium, and can thus produce superior estimates for enhanced sampling data such as parallel/simulated tempering, replica exchange, umbrella sampling, or metadynamics. In addition, dTRAM provides optimal estimates of Markov state models (MSMs) from the discretized state-space trajectories at all thermodynamic states. Under suitable conditions, these MSMs can be used to calculate kinetic quantities (e.g., rates, timescales). In the limit of a single thermodynamic state, dTRAM estimates a maximum likelihood reversible MSM, while in the limit of uncorrelated sampling data, dTRAM is identical to WHAM. dTRAM is thus a generalization to both estimators.

Original languageEnglish
Pages (from-to)214106
JournalThe Journal of Chemical Physics
Volume141
Issue number21
DOIs
Publication statusPublished - 7 Dec 2014

Keywords / Materials (for Non-textual outputs)

  • Ion Channels/chemistry
  • Likelihood Functions
  • Markov Chains
  • Molecular Dynamics Simulation
  • Proteins/chemistry
  • Thermodynamics

Fingerprint

Dive into the research topics of 'Statistically optimal analysis of state-discretized trajectory data from multiple thermodynamic states'. Together they form a unique fingerprint.

Cite this