Search

In this work, a new method for exploring conformational energy landscapes is described. The method, called transition-rapidly exploring random tree (T-RRT), combines ideas from statistical physics and robot path planning algorithms. A search tree is constructed on the conformational space starting from a given state. The tree expansion is driven by a double strategy: on the one hand, it is naturally biased toward yet unexplored regions of the space; on the other, a Monte Carlo-like transition test guides the expansion toward energetically favorable regions. The balance between these two strategies is automatically achieved due to a self-tuning mechanism. The method is able to efficiently find both energy minima and transition paths between them. As a proof of concept, the method is applied to two academic benchmarks and the alanine dipeptide. Copyright © 2011 Wiley Periodicals, Inc. more...

This paper presents a numerical method to compute all possible conformations of distance-constrained molecular loops, i.e., loops where some interatomic distances are held fixed, while others can vary. The method is general (it can be applied to single or multiple intermingled loops of arbitrary topology) and complete (it isolates all solutions, even if they form positive-dimensional sets). Generality is achieved by reducing the problem to finding all embeddings of a set of points constrained by pairwise distances, which can be formulated as computing the roots of a system of Cayley-Menger determinants. Completeness is achieved by expressing these determinants in Bernstein form and using a numerical algorithm that exploits such form to bound all root locations at any desired precision. The method is readily parallelizable, and the current implementation can be run on single- or multiprocessor machines. Experiments are included that show the method's performance on rigid loops, mobile loops, and multiloop molecules. In all cases, complete maps including all possible conformations are obtained, thus allowing an exhaustive analysis and visualization of all pseudo-rotation paths between different conformations satisfying loop closure. more...