An SDP-Based Divide-and-Conquer Algorithm for Large-Scale Noisy Anchor-Free Graph Realization

• Authors: Ngai-Hang Z. Leung;Kim-Chuan Toh
• Venue: SIAM Journal on Scientific Computing
• Year: 2009

We propose the DISCO algorithm for graph realization in $\mathbb{R}^d$, given sparse and noisy short-range intervertex distances as inputs. Our divide-and-conquer algorithm works as follows. When a group has a sufficiently small number of vertices, the basis step is to form a graph realization by solving a semidefinite program. The recursive step is to break a large group of vertices into two smaller groups with overlapping vertices. These two groups are solved recursively, and the subconfigurations are stitched together, using the overlapping atoms, to form a configuration for the larger group. At intermediate stages, the configurations are improved by gradient descent refinement. The algorithm is applied to the problem of determining protein moleculer structure. Tests are performed on molecules taken from the Protein Data Bank database. For each molecule, given 20-30% of the inter-atom distances less than 6Å that are corrupted by a high level of noise, DISCO is able to reliably and efficiently reconstruct the conformation of large molecules. In particular, given 30% of distances with 20% multiplicative noise, a 13000-atom conformation problem is solved within an hour with a root mean square deviation of 1.6Å.