A parallel version of the Branch & Prune algorithm for the Molecular Distance Geometry Problem

We consider the Molecular Distance Geometry Problem (MDGP), which is the problem of finding the conformation of a molecule from some known distances between its atoms. Such distances can be estimated by performing experiments of Nuclear Magnetic Resonance (NMR). Unfortunately, data obtained during these experiments are usually noisy and affected by errors. In particular, some of the estimated distances can be wrong, typically because assigned to the wrong pair of atoms. When particular assumptions are satisfied, the problem can be discretized, and solved by employing an ad-hoc algorithm called Branch & Prune (BP). However, this algorithm has been proved to be less efficient than a meta-heuristic algorithm when the percentage of wrong distances is large. We propose a parallel version of the BP algorithm which is able to handle this kind of instances. The scalability of the proposed algorithm allows for solving very large instances containing wrong distances. Implementation details of the algorithm in C/MPI are discussed, and computational experiments, performed on the nation-wide grid infrastructure Grid5000, are presented.

[1]  Gordon M. Crippen,et al.  Distance Geometry and Molecular Conformation , 1988 .

[2]  I. Kuntz,et al.  Distance geometry. , 1989, Methods in enzymology.

[3]  Kim-Chuan Toh,et al.  A Distributed SDP Approach for Large-Scale Noisy Anchor-Free Graph Realization with Applications to Molecular Conformation , 2008, SIAM J. Sci. Comput..

[4]  W. Gronwald,et al.  Automated structure determination of proteins by NMR spectroscopy , 2004 .

[5]  William Gropp,et al.  Users guide for mpich, a portable implementation of MPI , 1996 .

[6]  Leo Liberti,et al.  A Branch-and-Prune algorithm for the Molecular Distance Geometry Problem , 2008, Int. Trans. Oper. Res..

[7]  Carlile Lavor,et al.  On generating Instances for the Molecular Distance Geometry Problem , 2006 .

[8]  Franck Cappello,et al.  Grid'5000: A Large Scale And Highly Reconfigurable Experimental Grid Testbed , 2006, Int. J. High Perform. Comput. Appl..

[9]  Leo Liberti,et al.  Comparisons between an exact and a metaheuristic algorithm for the molecular distance geometry problem , 2009, GECCO.

[10]  Anthony Skjellum,et al.  A High-Performance, Portable Implementation of the MPI Message Passing Interface Standard , 1996, Parallel Comput..

[11]  T. Schlick Molecular modeling and simulation , 2002 .

[12]  Leo Liberti,et al.  Molecular distance geometry methods: from continuous to discrete , 2010, Int. Trans. Oper. Res..

[13]  Robin K. Harris,et al.  Encyclopedia of nuclear magnetic resonance , 1996 .

[14]  Jorge J. Moré,et al.  Global Continuation for Distance Geometry Problems , 1995, SIAM J. Optim..

[15]  Leo Liberti,et al.  The discretizable molecular distance geometry problem , 2006, Computational Optimization and Applications.

[16]  Tamar Schlick,et al.  Molecular Modeling and Simulation: An Interdisciplinary Guide , 2010 .

[17]  Jon M. Kleinberg,et al.  Reconstructing a three-dimensional model with arbitrary errors , 1996, STOC '96.

[18]  Leo Liberti,et al.  Molecular Distance Geometry Problem , 2009, Encyclopedia of Optimization.

[19]  Andrea Grosso,et al.  Solving molecular distance geometry problems by global optimization algorithms , 2009, Comput. Optim. Appl..