G. Mathematics of Computing
Refine
Document Type
- ZIB-Report (3)
Language
- English (3)
Has Fulltext
- yes (3)
Is part of the Bibliography
- no (3)
Keywords
- Markov State Models (1)
- Markov chains (1)
- Proteins, Conformation Space, Meshfree Methods (1)
- Voronoi (1)
- linear programming (1)
- meshfree methods (1)
- metastability (1)
Institute
- Computational Molecular Design (3) (remove)
Trajectory- or mesh-based methods for analyzing the dynamical behavior of large molecules tend to be impractical due to the curse of dimensionality - their computational cost increases exponentially with the size of the molecule. We propose a method to break the curse by a novel square root approximation of transition rates, Monte Carlo quadrature and a discretization approach based on solving linear programs. With randomly sampled points on the molecular energy landscape and randomly generated discretizations of the molecular configuration space as our initial data, we construct a matrix describing the transition rates between adjacent discretization regions. This transition rate matrix yields a Markov State Model of the molecular dynamics. We use Perron cluster analysis and coarse-graining techniques in order to identify metastable sets in configuration space and approximate the transition rates between the metastable sets. Application of our method to a simple energy landscape on a two-dimensional configuration space provides proof of concept and an example for which we compare the performance of different discretizations. We show that the computational cost of our method grows only polynomially with the size of the molecule. However, finding discretizations of higher-dimensional configuration spaces in which metastable sets can be identified remains a challenge.
This work aims to develop a new algorithm to calculate the free energy of water molecules by using a deterministic way.
For this purpose, we assume a closed system confined to a physical volume, having water molecules in gas phase.
To calculate the free energy of this sytem we utilized Mayer cluster expansion and the fluctuation free integration
method.
Particle methods have become indispensible in conformation dynamics to
compute transition rates in protein folding, binding processes and
molecular design, to mention a few.
Conformation dynamics requires at a decomposition of a molecule's position
space into metastable conformations.
In this paper, we show how this decomposition
can be obtained via the design of either ``soft'' or ``hard''
molecular conformations.
We show, that the soft approach results in a larger metastabilitiy of
the decomposition and is thus more advantegous. This is illustrated
by a simulation of Alanine Dipeptide.