A real-time proximity querying algorithm for haptic-based molecular docking

Intermolecular binding underlies every metabolic and regulatory processes of the cell, and the therapeutic and pharmacological properties of drugs. Molecular docking systems model and simulate these interactions in silico and allow us to study the binding process. Haptic-based docking provides an immersive virtual docking environment where the user can interact with and guide the molecules to their binding pose. Moreover, it allows human perception, intuition and knowledge to assist and accelerate the docking process, and reduces incorrect binding poses. Crucial for interactive docking is the real-time calculation of interaction forces. For smooth and accurate haptic exploration and manipulation, force-feedback cues have to be updated at a rate of 1 kHz. Hence, force calculations must be performed within 1ms. To achieve this, modern haptic-based docking approaches often utilize pre-computed force grids and linear interpolation. However, such grids are time-consuming to pre-compute (especially for large molecules), memory hungry, can induce rough force transitions at cell boundaries and cannot be applied to flexible docking. Here we propose an efficient proximity querying method for computing intermolecular forces in real time. Our motivation is the eventual development of a haptic-based docking solution that can model molecular flexibility. Uniquely in a haptics application we use octrees to decompose the 3D search space in order to identify the set of interacting atoms within a cut-off distance. Force calculations are then performed on this set in real time. The implementation constructs the trees dynamically, and computes the interaction forces of large molecular structures (i.e. consisting of thousands of atoms) within haptic refresh rates. We have implemented this method in an immersive, haptic-based, rigid-body, molecular docking application called Haptimol_RD. The user can use the haptic device to orientate the molecules in space, sense the interaction forces on the device, and guide the molecules to their binding pose. Haptimol_RD is designed to run on consumer level hardware, i.e. there is no need for specialized/proprietary hardware

The Journal of Computer-Aided Molecular Design: a bibliometric note

Summarizes the articles in, and the citations to, volumes 2-24 of the Journal of Computer-Aided Molecular Design. The citations to the journal come from almost 2000 different sources that span a very wide range of academic subjects, with the most heavily cited articles being descriptions of software systems and of computational methods

The Predicted Binding Site and Dynamics of Peptide Inhibitors to the Methuselah GPCR from Drosophila melanogaster

Peptide inhibitors of Methuselah (Mth), a G protein-coupled receptor (GPCR), were reported that can extend the life span of Drosophila melanogaster. Mth is a class B GPCR, which is characterized by a large, N-terminal ectodomain that is often involved with ligand recognition. The crystal structure of the Mth ectodomain, which binds to the peptide inhibitors with high affinity, was previously determined. Here we report the predicted structures for RWR motif peptides in complex with the Mth ectodomain. We studied representatives of both Pro-class and Arg-class RWR motif peptides and identified ectodomain residues Asp139, Phe130, Asp127, and Asp78 as critical in ligand binding. To validate these structures, we predicted the effects of various ligand mutations on the structure and binding to Mth. The binding of five mutant peptides to Mth was characterized experimentally by surface plasmon resonance, revealing measured affinities that are consistent with predictions. The electron density map calculated from our MD structure compares well with the experimental map of a previously determined peptide/Mth crystal structure and could be useful in refining the current low-resolution data. The elucidation of the ligand binding site may be useful in analyzing likely binding sites in other class B GPCRs

LightDock: a new multi-scale approach to protein–protein docking

Computational prediction of protein–protein complex structure by docking can provide structural and mechanistic insights for protein interactions of biomedical interest. However, current methods struggle with difficult cases, such as those involving flexible proteins, low-affinity complexes or transient interactions. A major challenge is how to efficiently sample the structural and energetic landscape of the association at different resolution levels, given that each scoring function is often highly coupled to a specific type of search method. Thus, new methodologies capable of accommodating multi-scale conformational flexibility and scoring are strongly needed. We describe here a new multi-scale protein–protein docking methodology, LightDock, capable of accommodating conformational flexibility and a variety of scoring functions at different resolution levels. Implicit use of normal modes during the search and atomic/coarse-grained combined scoring functions yielded improved predictive results with respect to state-of-the-art rigid-body docking, especially in flexible cases.B.J-G was supported by a FPI fellowship from the Spanish Ministry of Economy and Competitiveness. This work was supported by I+D+I Research Project grants BIO2013-48213-R and BIO2016-79930-R from the Spanish Ministry of Economy and Competitiveness. This work is partially supported by the European Union H2020 program through HiPEAC (GA 687698), by the Spanish Government through Programa Severo Ochoa (SEV-2015-0493), by the Spanish Ministry of Science and Technology (TIN2015-65316-P) and the Departament d’Innovació, Universitats i Empresa de la Generalitat de Catalunya, under project MPEXPAR: Models de Programaciói Entorns d’Execució Paral·lels (2014-SGR-1051).Peer ReviewedPostprint (author's final draft