Visual Analysis of Molecular Conformations by Means of a Dynamic Density Mixture Model

We propose an approach for transforming the sampling of a molecular conformation distribution into an analytical model based on Hidden Markov Models. The model describes the sampled shape density as a mixture of multivariate unimodal densities. Thus, it delivers an interpretation of the sampled density as a set of typical shapes that appear with different probabilities and are characterized by their geometry, their variability and transition probabilities between the shapes. The gained model is used to identify atom groups of constant shape that are connected by metastable torsion angles. Based on this description an alignment for the original sampling is computed. As it takes into account the different shapes contained in the sampled set, this alignment allows to compute reasonable average shapes and meaningful shape density plots. Furthermore, it enables us to visualize typical conformations.

[1]  Charles Hansen,et al.  The Visualization Handbook , 2011 .

[2]  C. Schütte Conformational Dynamics: Modelling, Theory, Algorithm, and Application to Biomolecules , 1999 .

[3]  Heermann,et al.  Hybrid Monte Carlo method for condensed-matter systems. , 1992, Physical review. B, Condensed matter.

[4]  P. Deuflhard,et al.  Robust Perron cluster analysis in conformation dynamics , 2005 .

[5]  Maryam Shayegan Hastings,et al.  She Does Math!: Mathematics and Computer Science , 1995 .

[6]  Alexander Fischer,et al.  Identification of biomolecular conformations from incomplete torsion angle observations by hidden markov models , 2007, J. Comput. Chem..

[7]  Hans-Christian Hege,et al.  amira: A Highly Interactive System for Visual Data Analysis , 2005, The Visualization Handbook.

[8]  Daniel Baum,et al.  Visualizing dynamic molecular conformations , 2002, IEEE Visualization, 2002. VIS 2002..

[9]  Ross Moore,et al.  Applied Mathematics Entering the 21st Century , 2004 .

[10]  W. Kabsch A solution for the best rotation to relate two sets of vectors , 1976 .

[11]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[12]  P. Deuflhard,et al.  Molecular Conformation Dynamics and Computational Drug Design 1 , 2003 .

[13]  Robert van Liere,et al.  Interactive visualization of protein dynamics , 2000 .

[14]  Lawrence R. Rabiner,et al.  A tutorial on hidden Markov models and selected applications in speech recognition , 1989, Proc. IEEE.

[15]  P. Deuflhard,et al.  A Direct Approach to Conformational Dynamics Based on Hybrid Monte Carlo , 1999 .

[16]  P. Deuflharda,et al.  Identification of almost invariant aggregates in reversible nearly uncoupled Markov chains , 2000 .

[17]  W. Kabsch A discussion of the solution for the best rotation to relate two sets of vectors , 1978 .

[18]  H. Berendsen,et al.  Essential dynamics of proteins , 1993, Proteins.

[19]  Marcus Weber,et al.  Clustering by using a simplex structure , 2003 .

[20]  William H. Press,et al.  Numerical recipes in C , 2002 .

[21]  Xavier Pennec,et al.  Multiple Registration and Mean Rigid Shape - Application to the 3D case , 1996 .