Abstract:Motivation: Backbone resonance assignment is a critical bottleneck in studies of protein structure, dynamics and interactions by nuclear magnetic resonance (NMR) spectroscopy. A minimalist approach to assignment, which we call 'contact-based', seeks to dramatically reduce experimental time and expense by replacing the standard suite of through-bond experiments with the through-space (nuclear Overhauser enhancement spectroscopy, NOESY) experiment. In the contact-based approach, spectral data are represented in a graph with vertices for putative residues (of unknown relation to the primary sequence) and edges for hypothesized NOESY interactions, such that observed spectral peaks could be explained if the residues were 'close enough'. Due to experimental ambiguity, several incorrect edges can be hypothesized for each spectral peak. An assignment is derived by identifying consistent patterns of edges (e.g. for alpha-helices and beta-sheets) within a graph and by mapping the vertices to the primary sequence. The key algorithmic challenge is to be able to uncover these patterns even when they are obscured by significant noise. Results: This paper develops, analyzes and applies a novel algorithm for the identification of polytopes representing consistent patterns of edges in a corrupted NOESY graph. Our randomized algorithm aggregates simplices into polytopes and fixes inconsistencies with simple local modifications, called rotations, that maintain most of the structure already uncovered. In characterizing the effects of experimental noise, we employ an NMR-specific random graph model in proving that our algorithm gives optimal performance in expected polynomial time, even when the input graph is significantly corrupted. We confirm this analysis in simulation studies with graphs corrupted by up to 500% noise. Finally, we demonstrate the practical application of the algorithm on several experimental beta-sheet datasets. Our approach is able to eliminate a large majority of noise edges and to uncover large consistent sets of interactions. Availability: Our algorithm has been implemented in the platform-independent Python code. The software can be freely obtained for academic use by request from the authors.

Voronota-LT: efficient, flexible and solvent-aware tessellation-based analysis of atomic interactions

An extension to VORO++ for multithreaded computation of Voronoi cells

Voronoi topology in two dimensions: theory, algorithms, and applications

ParVoro++: A Scalable Parallel Algorithm for Constructing 3D Voronoi Tessellations Based on Kd-Tree Decomposition

High-order finite element method for atomic structure calculations

ezAlign: A Tool for Converting Coarse-Grained Molecular Dynamics Structures to Atomistic Resolution for Multiscale Modeling

votess: A multi-target, GPU-capable, parallel Voronoi tessellator

Efficient Molecular Approach to Quantifying Solvent-Mediated Interactions

On Open and Strong-Scaling Tools for Atom Probe Crystallography: High-Throughput Methods for Indexing Crystal Structure and Orientation

Voronoi Graph -- Improved raycasting and integration schemes for high dimensional Voronoi diagrams

Linear-scaling local natural orbital CCSD(T) approach for open-shell systems: algorithm, benchmarks, and large-scale applications

Fast, scalable, and interactive software for Landau-de Gennes numerical modeling of nematic topological defects

Characterizing structural features of two-dimensional particle systems through Voronoi topology

Robust structural identification via polyhedral template matching

From Peptides to Nanostructures: A Euclidean Transformer for Fast and Stable Machine Learned Force Fields

Advancing non-atom-centered basis methods for more accurate interaction energies: benchmarks and large-scale applications

Level Set Restricted Voronoi Tessellation for Large scale Spatial Statistical Analysis

Driving atomic structures of molecules, crystals, and complex systems with local similarity kernels

Martinize2 and Vermouth: Unified Framework for Topology Generation

Toward transferable interatomic van der Waals interactions without electrons: The role of multipole electrostatics and many-body dispersion

An efficient randomized algorithm for contact-based NMR backbone resonance assignment