skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Collision-free poisson motion planning in ultra high-dimensional molecular conformation spaces

Abstract

The function of protein, RNA, and DNA is modulated by fast, dynamic exchanges between three-dimensional conformations. Conformational sampling of biomolecules with exact and nullspace inverse kinematics, using rotatable bonds as revolute joints and noncovalent interactions as holonomic constraints, can accurately characterize these native ensembles. However, sampling biomolecules remains challenging owing to their ultra-high dimensional configuration spaces, and the requirement to avoid (self-) collisions, which results in low acceptance rates. In this paper, we present two novel mechanisms to overcome these limitations. First, we introduce temporary constraints between near-colliding links. The resulting constraint varieties instantaneously redirect the search for collision-free conformations, and couple motions between distant parts of the linkage. Second, we adapt a randomized Poisson-disk motion planner, which prevents local oversampling and widens the search, to ultra-high dimensions. Tests on several model systems show that the sampling acceptance rate can increase from 16% to 70%, and that the conformational coverage in loop modeling measured as average closeness to existing loop conformations doubled. Finally, correlated protein motions identified with our algorithm agree with those from MD simulations.

Authors:
ORCiD logo [1]; ORCiD logo [2]; ORCiD logo [3]
  1. Stanford Univ., CA (United States). Molecular and Cellular Physiology; SLAC National Accelerator Lab., Menlo Park, CA (United States). Bioscience Division
  2. Univ. of Erlangen-Nuremberg, Erlangen (Germany). Chair of Applied Dynamics
  3. SLAC National Accelerator Lab., Menlo Park, CA (United States). Bioscience Division
Publication Date:
Research Org.:
SLAC National Accelerator Laboratory (SLAC), Menlo Park, CA (United States); Stanford Univ., CA (United States); Univ. of Erlangen-Nuremberg, Erlangen (Germany)
Sponsoring Org.:
USDOE; Deutsche Telekom Stiftung (Germany); Novo Nordisk Foundation (Denmark)
OSTI Identifier:
1471534
Grant/Contract Number:  
AC02-76SF00515; NNF15OC0015268
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
Journal of Computational Chemistry
Additional Journal Information:
Journal Volume: 39; Journal Issue: 12; Journal ID: ISSN 0192-8651
Publisher:
Wiley
Country of Publication:
United States
Language:
English
Subject:
59 BASIC BIOLOGICAL SCIENCES; 37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CHEMISTRY; 97 MATHEMATICS AND COMPUTING; high dimensional motion and path planning; collision-avoidance; poisson sampling; computational biology; inverse kinematics; molecular rigidity

Citation Formats

Fonseca, Rasmus, Budday, Dominik, and van den Bedem, Henry. Collision-free poisson motion planning in ultra high-dimensional molecular conformation spaces. United States: N. p., 2018. Web. doi:10.1002/jcc.25138.
Fonseca, Rasmus, Budday, Dominik, & van den Bedem, Henry. Collision-free poisson motion planning in ultra high-dimensional molecular conformation spaces. United States. https://doi.org/10.1002/jcc.25138
Fonseca, Rasmus, Budday, Dominik, and van den Bedem, Henry. 2018. "Collision-free poisson motion planning in ultra high-dimensional molecular conformation spaces". United States. https://doi.org/10.1002/jcc.25138. https://www.osti.gov/servlets/purl/1471534.
@article{osti_1471534,
title = {Collision-free poisson motion planning in ultra high-dimensional molecular conformation spaces},
author = {Fonseca, Rasmus and Budday, Dominik and van den Bedem, Henry},
abstractNote = {The function of protein, RNA, and DNA is modulated by fast, dynamic exchanges between three-dimensional conformations. Conformational sampling of biomolecules with exact and nullspace inverse kinematics, using rotatable bonds as revolute joints and noncovalent interactions as holonomic constraints, can accurately characterize these native ensembles. However, sampling biomolecules remains challenging owing to their ultra-high dimensional configuration spaces, and the requirement to avoid (self-) collisions, which results in low acceptance rates. In this paper, we present two novel mechanisms to overcome these limitations. First, we introduce temporary constraints between near-colliding links. The resulting constraint varieties instantaneously redirect the search for collision-free conformations, and couple motions between distant parts of the linkage. Second, we adapt a randomized Poisson-disk motion planner, which prevents local oversampling and widens the search, to ultra-high dimensions. Tests on several model systems show that the sampling acceptance rate can increase from 16% to 70%, and that the conformational coverage in loop modeling measured as average closeness to existing loop conformations doubled. Finally, correlated protein motions identified with our algorithm agree with those from MD simulations.},
doi = {10.1002/jcc.25138},
url = {https://www.osti.gov/biblio/1471534}, journal = {Journal of Computational Chemistry},
issn = {0192-8651},
number = 12,
volume = 39,
place = {United States},
year = {Fri Jan 05 00:00:00 EST 2018},
month = {Fri Jan 05 00:00:00 EST 2018}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 3 works
Citation information provided by
Web of Science

Figures / Tables:

Figure  1 Figure 1: Kinematic representation of molecules and constraints. a) A molecular graph with double bonds and partial double bonds highlighted. After edge-contractions these atoms become a single rigid body vertex. The remaining bonds can rotate around the bond axis and are described by dihedral angles. b) Small section of amore » protein molecule with hydrogen bonds marked as purple lines and a nitrogen-hydrogen clash marked with van der Waals spheres and the bisecting plane. c) A dynamic clash-avoiding constraint (dCC) allows individual motions of the near-clashing atoms p1 and p2 along directions tc in the plane with normal vector nc, but only a joint move along nc. This allows atoms to slide past each other, but prevents them from getting closer.« less

Save / Share:

Works referenced in this record:

Rapid Sampling of Molecular Motions with Prior Information Constraints
journal, February 2009


Random Coordinate Descent with Spinor-matrices and Geometric Filters for Efficient Loop Closure
journal, February 2013


Computational Models of Protein Kinematics and Dynamics: Beyond Simulation
journal, July 2012


Transition Networks for the Comprehensive Characterization of Complex Conformational Change in Proteins
journal, May 2006


Correlated Motions and Interactions at the Onset of the DNA-Induced Partial Unfolding of Ets-1
journal, February 2009


Algorithm and Data Structures for Efficient Energy Maintenance during Monte Carlo Simulation of Proteins
journal, October 2004


Stochastic sampling in computer graphics
journal, January 1986


Sampling-based exploration of folded state of a protein under kinematic and geometric constraints
journal, October 2011


Probing RNA Native Conformational Ensembles with Structural Constraints
journal, May 2016


Frustration-guided motion planning reveals conformational transitions in proteins: BUDDAY et al.
journal, July 2017


A Comparison of Methods for Generating Poisson Disk Distributions
journal, March 2008


An NMA-guided path planning approach for computing large-amplitude conformational changes in proteins
journal, July 2007


Integrative, dynamic structural biology at atomic resolution—it's about time
journal, March 2015


Bounding Volumes for Proteins: A Comparative Study
journal, October 2012


A Data-Driven Evolutionary Algorithm for Mapping Multibasin Protein Energy Landscapes
journal, September 2015


Coupled Motions in β 2 AR:Gαs Conformational Ensembles
journal, February 2016


A Distal Mutation Perturbs Dynamic Amino Acid Networks in Dihydrofolate Reductase
journal, June 2013


A stochastic roadmap method to model protein structural transitions
journal, December 2015


A path planning approach for computing large-amplitude motions of flexible molecules
journal, June 2005


Smooth continuous transition between tasks on a kinematic control level: Obstacle avoidance as a control problem
journal, September 2013


Simulating Protein Motions with Rigidity Analysis
journal, July 2007


Motion planning algorithms for molecular simulations: A survey
journal, July 2012


The experimental survey of protein-folding energy landscapes
journal, August 2005


Automated design of the surface positions of protein helices: Protein surface design
journal, June 1997


Automated identification of functional dynamic contact networks from X-ray crystallography
journal, August 2013


Resultants and loop closure
journal, January 2005


Probabilistic roadmaps for path planning in high-dimensional configuration spaces
journal, January 1996


A spatial data structure for fast Poisson-disk sample generation
journal, July 2006


Modeling protein conformational transitions by a combination of coarse-grained normal mode analysis and robotics-inspired methods
journal, January 2013


Biomolecular Simulation: A Computational Microscope for Molecular Biology
journal, June 2012


Characterizing RNA ensembles from NMR data with kinematic models
journal, August 2014


Geometric analysis characterizes molecular rigidity in generic and non-generic protein configurations
journal, October 2015


Spheres, molecules, and hidden surface removal
conference, January 1994


Exploring the gas access routes in a [NiFeSe] hydrogenase using crystals pressurized with krypton and oxygen
journal, August 2020


A spatial data structure for fast Poisson-disk sample generation
conference, January 2006


Structure and Dynamics of the G121V Dihydrofolate Reductase Mutant: Lessons from a Transition-State Inhibitor Complex
text, January 2012


Tracing conformational changes in proteins
conference, November 2009


Spheres, molecules, and hidden surface removal
journal, October 1998


Correlated Motions and Interactions at the Onset of the DNA-Induced Partial Unfolding of Ets-1
journal, February 2009


Spheres, molecules, and hidden surface removal
journal, October 1998


Automated identification of functional dynamic contact networks from X-ray crystallography
journal, August 2013


Characterizing RNA ensembles from NMR data with kinematic models
journal, August 2014


Tracing conformational changes in proteins
journal, January 2010


Rapid Sampling of Molecular Motions with Prior Information Constraints
journal, February 2009


Figures/Tables have been extracted from DOE-funded journal article accepted manuscripts.