Email updates

Keep up to date with the latest news and content from BMC Bioinformatics and BioMed Central.

This article is part of the supplement: Highlights from the Ninth International Society for Computational Biology (ISCB) Student Council Symposium 2013

Open Access Meeting abstract

Towards breaking the curse of dimensionality in computational methods for the conformational analysis of molecules

Han Cheng Lie

Author Affiliations

Institute of Mathematics – Freie Universität Berlin, Berlin, Germany

BMC Bioinformatics 2014, 15(Suppl 3):A2  doi:10.1186/1471-2105-15-S3-A2

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2105/15/S3/A2


Published:11 February 2014

© 2014 Lie; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Background

In computational molecular biology, the conformational analysis of a molecule refers to the identification of conformations - the clusters of molecular configurations sharing a large-scale geometric structure - and the description of conformation dynamics. These data help to describe the structures and functions of molecules which are relevant in biology. In applications such as computational drug design, they help to identify suitable compounds for disrupting processes involved in disease mechanisms. Today, many computational methods for conformational analysis are mesh-based: they involve discretizing configuration space using regular lattices, and simulating trajectories on energy landscapes. These methods suffer from the ‘curse of dimensionality’ in that the cost increases exponentially with the number of atoms in the molecule. Given that many important molecules in biology - such as proteins - have thousands of atoms, the curse of dimensionality severely limits the applicability of these methods. Finding an alternative approach to conformational analysis that does not suffer from the curse of dimensionality would enlarge the set of molecules and biological processes which may be studied using computational methods.

Results

We construct a mesh-free method [1] using random sampling of energy landscapes, Monte Carlo quadrature, and coarse-graining [2]. We provide proof of concept by testing the method on a toy model of a molecule from the literature and find that the results of our method are comparable to those obtained by an exact, mesh-based method. Using results from linear programming theory and discrete geometry, we establish that the computational cost of our method increases only polynomially with the size of the molecule – an improvement over the exponential increase in the case of mesh-based methods.

Conclusions

Our findings demonstrate that one does not need meshes for conformational analysis. For the specific problem we study in our paper – that of finding conformations and computing conformation transition rates – our findings suggest that mesh-free methods may be just as accurate as mesh-based methods, while costing less computationally. Our findings emphasize the utility of using ideas from different mathematical fields when devising new methods to tackle computational challenges.

References

  1. Lie HC, Fackeldey K, Weber M: A Square Root Approximation of Transition Rates for a Markov State Model.

    SIAM Journal of Matrix Analysis and Applications 2013, 34(2):738-756. Publisher Full Text OpenURL

  2. Kube S, Weber M: A Coarse-Graining Method for the Identification of Transition Rates between Molecular Conformations.

    Journal of Chemical Physics 2007, 126(2):024103. PubMed Abstract | Publisher Full Text OpenURL