Full cyclic coordinate descent: solving the protein loop closure problem in C<it>α </it>space

• Main Authors: Hamelryck Thomas, Boomsma Wouter
• Format: Article
• Language: English
• Published: BMC 2005-06-01
• Series: BMC Bioinformatics
• Online Access: http://www.biomedcentral.com/1471-2105/6/159

<p>Abstract</p> <p>Background</p> <p>Various forms of the so-called <it>loop closure problem </it>are crucial to protein structure prediction methods. Given an N- and a C-terminal end, the problem consists of finding a suitable segment of a certain length that bridges the ends seamlessly.</p> <p>In homology modelling, the problem arises in predicting loop regions. In <it>de novo </it>protein structure prediction, the problem is encountered when implementing local moves for Markov Chain Monte Carlo simulations.</p> <p>Most loop closure algorithms keep the bond angles fixed or semi-fixed, and only vary the dihedral angles. This is appropriate for a full-atom protein backbone, since the bond angles can be considered as fixed, while the (<it>φ</it>, <it>ψ</it>) dihedral angles are variable. However, many <it>de novo </it>structure prediction methods use protein models that only consist of C<it>α </it>atoms, or otherwise do not make use of all backbone atoms. These methods require a method that alters both bond and dihedral angles, since the pseudo bond angle between three consecutive C<it>α </it>atoms also varies considerably.</p> <p>Results</p> <p>Here we present a method that solves the loop closure problem for C<it>α </it>only protein models. We developed a variant of Cyclic Coordinate Descent (CCD), an inverse kinematics method from the field of robotics, which was recently applied to the loop closure problem. Since the method alters both bond and dihedral angles, which is equivalent to applying a full rotation matrix, we call our method Full CCD (FCDD). FCCD replaces CCD's vector-based optimization of a rotation around an axis with a singular value decomposition-based optimization of a general rotation matrix. The method is easy to implement and numerically stable.</p> <p>Conclusion</p> <p>We tested the method's performance on sets of random protein C<it>α </it>segments between 5 and 30 amino acids long, and a number of loops of length 4, 8 and 12. FCCD is fast, has a high success rate and readily generates conformations close to those of real loops. The presence of constraints on the angles only has a small effect on the performance. A reference implementation of FCCD in Python is available as supplementary information.</p>