Abstract:Computational methods for creating in silico libraries of molecular descriptors (e.g., collision cross sections) are becoming increasingly prevalent due to the limited number of authentic reference materials available for traditional library building. These so-called "reference-free metabolomics" methods require sampling sets of molecular conformers in order to produce high accuracy property predictions. Due to the computational cost of the subsequent calculations for each conformer, there is a need to sample the most relevant subset and avoid repeating calculations on conformers that are nearly identical. The goal of this study is to introduce a heuristic method of finding the most dissimilar conformers from a larger population in order to help speed up reference-free calculation methods and maintain a high property prediction accuracy. Finding the set of the n items most dissimilar from each other out of a larger population becomes increasingly difficult and computationally expensive as either n or the population size grows large. Because there exists a pairwise relationship between each item and all other items in the population, finding the set of the n most dissimilar items is different than simply sorting an array of numbers. For instance, if you have a set of the most dissimilar n = 4 items, one or more of the items from n = 4 might not be in the set n = 5. An exact solution would have to search all possible combinations of size n in the population exhaustively. We present an open-source software called similarity downselection (SDS), written in Python and freely available on GitHub. SDS implements a heuristic algorithm for quickly finding the approximate set(s) of the n most dissimilar items. We benchmark SDS against a Monte Carlo method, which attempts to find the exact solution through repeated random sampling. We show that for SDS to find the set of n most dissimilar conformers, our method is not only orders of magnitude faster, but it is also more accurate than running Monte Carlo for 1,000,000 iterations, each searching for set sizes n = 3-7 out of a population of 50,000. We also benchmark SDS against the exact solution for example small populations, showing that SDS produces a solution close to the exact solution in these instances. Using theoretical approaches, we also demonstrate the constraints of the greedy algorithm and its efficacy as a ratio to the exact solution.

Conformational Pruning via the Permutation Invariant Root Mean Square Deviation of Atomic Positions

Canonized then Minimized RMSD for Three-Dimensional Structures

A new method on reconstructing protein structure from NOESY distances

Scaled Alternating Steepest Descent Algorithm Applied for Protein Structure Determination from Nuclear Magnetic Resonance Data.

Quantifying Near-Symmetric Molecular Distortion Using Symmetry-Coordinate Structural Decomposition

Molecular Gas Phase Conformational Ensembles

Exact solution for minimization of root mean square deviation with G-RMSD to determine molecular similarity

Measuring Intrinsic Disorder and Tracking Conformational Transitions Using Rosetta ResidueDisorder

Efficient Conformational Search Based on Structural Dissimilarity Sampling: Applications for Reproducing Structural Transitions of Proteins.

Protein Retrieval via Integrative Molecular Ensembles (PRIME) through extended similarity indices

Equivalent Distance Geometry Error for Molecular Conformation Comparison

Structure optimization in an off-lattice protein model

Understanding and Quantifying Molecular Flexibility: Torsion Angular Bin Strings

Eurecon: Equidistant uniform rigid-body ensemble constructor

Experimental determination of optimal root-mean-square deviations of macromolecular bond lengths and angles from their restrained ideal values

First-Principles Molecular Structure Search with a Genetic Algorithm

Efficient molecular conformation generation with quantum-inspired algorithm

Similarity Downselection: Finding the n Most Dissimilar Molecular Conformers for Reference-Free Metabolomics

Regularized Molecular Conformation Fields

Quantum refinement in real and reciprocal space using the Phenix and ORCA software