Refine
Document Type
- Article (4) (remove)
Has Fulltext
- no (4)
Is part of the Bibliography
- no (4)
Keywords
- MSM (2)
- Convex Optimization (1)
- Girsanov (1)
- Reversible Markov Chain (1)
- Reweighting (1)
Institute
In recent years Markov State Models (MSMs) have attracted a consid-
erable amount of attention with regard to modelling conformation changes and associated function of biomolecular systems. They have been used successfully, e.g., for peptides including time-resolved spectroscopic experiments, protein function and protein folding , DNA and RNA, and ligand-receptor interaction in drug design and more complicated multivalent scenarios. In this article a novel reweighting scheme is introduced that allows to construct an MSM for certain molecular system out of an MSM for a similar system. This permits studying how molecular properties on long timescales differ between similar molecular systems without performing full molecular dynamics simulations for each system under con-
sideration. The performance of the reweighting scheme is illustrated for simple test cases including one where the main wells of the respective energy landscapes are located differently and an alchemical transformation of butane to pentane where the dimension of the state space is changed.
In many applications one is interested to compute transition probabilities of a Markov chain.
This can be achieved by using Monte Carlo methods with local or global sampling points.
In this article, we analyze the error by the difference in the $L^2$ norm between the true transition probabilities and the approximation
achieved through a Monte Carlo method.
We give a formula for the error for Markov chains with locally computed sampling points. Further, in the case of reversible Markov chains, we will deduce a formula for the error when sampling points are computed globally.
We will see that in both cases the error itself can be approximated with Monte Carlo methods.
As a consequence of the result, we will derive surprising properties of reversible Markov chains.
Reversible Markov chains are the basis of many applications. However, computing transition probabilities by a finite sampling of a Markov chain can lead to truncation errors. Even if the original Markov chain is reversible, the approximated Markov chain might be non-reversible and will lose important properties, like the real valued spectrum. In this paper, we show how to find the closest reversible Markov chain to a given transition matrix. It turns out that this matrix can be computed by solving a convex minimization problem.