Applying forces to elastic network models of large biomolecules using a haptic feedback device

Elastic network models of biomolecules have proved to be relatively good at predicting global conformational changes particularly in large systems. Software that facilitates rapid and intuitive exploration of conformational change in elastic network models of large biomolecules in response to externally applied forces would therefore be of considerable use, particularly if the forces mimic those that arise in the interaction with a functional ligand. We have developed software that enables a user to apply forces to individual atoms of an elastic network model of a biomolecule through a haptic feedback device or a mouse. With a haptic feedback device the user feels the response to the applied force whilst seeing the biomolecule deform on the screen. Prior to the interactive session normal mode analysis is performed, or pre-calculated normal mode eigenvalues and eigenvectors are loaded. For large molecules this allows the memory and number of calculations to be reduced by employing the idea of the important subspace, a relatively small space of the first M lowest frequency normal mode eigenvectors within which a large proportion of the total fluctuation occurs. Using this approach it was possible to study GroEL on a standard PC as even though only 2.3% of the total number of eigenvectors could be used, they accounted for 50% of the total fluctuation. User testing has shown that the haptic version allows for much more rapid and intuitive exploration of the molecule than the mouse version

Advances in Human-Protein Interaction - Interactive and Immersive Molecular Simulations

International audienc

Software Introduction: Methodological advances for interacting with biomolecules using haptics

Over the past 15 years we have been developing tools for interacting with biomolecules using haptics. Interactions with biomolecules in the virtual world are made via a haptic-feedback device that is able to resist inputs from the user or even act to move the user’s hand in response to molecular forces. Here we highlight the key methodological advances made in the development of these tools including Haptimol ISAS, a tool for interacting with a molecule’s solvent accessible surface, Haptimol ENM, a tool for applying forces to an elastic network model of a biomolecule, DockIT (formerly Haptimol RD), for interactive rigid docking, and Haptimol FlexiDock, for interactive docking that models flexibility in the receptor molecule

Exploration of Reaction Pathways and Chemical Transformation Networks

For the investigation of chemical reaction networks, the identification of all relevant intermediates and elementary reactions is mandatory. Many algorithmic approaches exist that perform explorations efficiently and automatedly. These approaches differ in their application range, the level of completeness of the exploration, as well as the amount of heuristics and human intervention required. Here, we describe and compare the different approaches based on these criteria. Future directions leveraging the strengths of chemical heuristics, human interaction, and physical rigor are discussed.Comment: 48 pages, 4 figure

Adaptive GPU-accelerated force calculation for interactive rigid molecular docking using haptics

Molecular docking systems model and simulate in silico the interactions of intermolecular binding. Haptics-assisted docking enables the user to interact with the simulation via their sense of touch but a stringent time constraint on the computation of forces is imposed due to the sensitivity of the human haptic system. To simulate high fidelity smooth and stable feedback the haptic feedback loop should run at rates of 500 Hz to 1 kHz. We present an adaptive force calculation approach that can be executed in parallel on a wide range of Graphics Processing Units (GPUs) for interactive haptics-assisted docking with wider applicability to molecular simulations. Prior to the interactive session either a regular grid or an octree is selected according to the available GPU memory to determine the set of interatomic interactions within a cutoff distance. The total force is then calculated from this set. The approach can achieve force updates in less than 2 ms for molecular structures comprising hundreds of thousands of atoms each, with performance improvements of up to 90 times the speed of current CPU-based force calculation approaches used in interactive docking. Furthermore, it overcomes several computational limitations of previous approaches such as pre-computed force grids, and could potentially be used to model receptor flexibility at haptic refresh rates

Simulating molecular docking with haptics

Intermolecular binding underlies various 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 the study of the binding process. In molecular docking, haptics enables the user to sense the interaction forces and intervene cognitively in the docking process. Haptics-assisted docking systems provide an immersive virtual docking environment where the user can interact with the molecules, feel the interaction forces using their sense of touch, identify visually the binding site, and guide the molecules to their binding pose. Despite a forty-year research e�ort however, the docking community has been slow to adopt this technology. Proprietary, unreleased software, expensive haptic hardware and limits on processing power are the main reasons for this. Another signi�cant factor is the size of the molecules simulated, limited to small molecules. The focus of the research described in this thesis is the development of an interactive haptics-assisted docking application that addresses the above issues, and enables the rigid docking of very large biomolecules and the study of the underlying interactions. Novel methods for computing the interaction forces of binding on the CPU and GPU, in real-time, have been developed. The force calculation methods proposed here overcome several computational limitations of previous approaches, such as precomputed force grids, and could potentially be used to model molecular exibility at haptic refresh rates. Methods for force scaling, multipoint collision response, and haptic navigation are also reported that address newfound issues, particular to the interactive docking of large systems, e.g. force stability at molecular collision. The i ii result is a haptics-assisted docking application, Haptimol RD, that runs on relatively inexpensive consumer level hardware, (i.e. there is no need for specialized/proprietary hardware)

Virtual environment for studying the docking interactions of rigid biomolecules with haptics

Haptic technology facilitates user interaction with the virtual world via the sense of touch. In molecular docking, haptics enables the user to sense the interaction forces during the docking process. Here we describe a haptics-assisted interactive software tool, called Haptimol RD, for the study of docking interactions. By utilising GPU-accelerated proximity querying methods very large systems can now be studied. Methods for force scaling, multipoint collision response and haptic navigation are described that address force stability issues that are particular to the interactive docking of large systems. Thus Haptimol RD expands, for the first time, the use of interactive biomolecular haptics to the study of protein-protein interactions. Unlike existing approaches, Haptimol RD is designed to run on relatively inexpensive consumer-level hardware and is freely available to the community

Haptic-assisted interactive molecular docking incorporating receptor flexibility

Haptic-assisted interactive docking tools immerse the user in an environment where intuition and knowledge can be used to help guide the docking process. Here we present such a tool where the user “holds” a rigid ligand via a haptic device through which they feel interaction forces with a flexible receptor biomolecule. To ensure forces transmitted through the haptic device are smooth and stable, they must be updated at a rate greater than 500 Hz. Due to this time constraint, the majority of haptic docking tools do not attempt to model the conformational changes that would occur when molecules interact during binding. Our haptic-assisted docking tool, “Haptimol Flexidock”, models a receptor’s conformational response to forces of interaction with a ligand whilst maintaining the required haptic refresh rate. In order to model receptor flexibility we use the method of linear response for which we determine the variance-covariance matrix of atomic fluctuations from the trajectory of an explicit-solvent Molecular Dynamics simulation of the ligand-free receptor molecule. Key to satisfying the time constraint is an eigenvector decomposition of the variance-covariance matrix which enables a good approximation to the conformational response of the receptor to be calculated rapidly. This exploits a feature of protein dynamics whereby most fluctuation occurs within a relatively small subspace. The method is demonstrated on Glutamine Binding Protein in interaction with glutamine, and Maltose Binding Protein in interaction with maltose. For both proteins, the movement that occurs when the ligand is docked near to its binding site matches the experimentally determined movement well. It is thought that this tool will be particularly useful for structure-based drug design

Interactive molecular docking with haptics and advanced graphics

Biomolecular interactions underpin many of the processes that make up life. Molecular docking is the study of these interactions in silico. Interactive docking applications put the user in control of the docking process, allowing them to use their knowledge and intuition to determine how molecules bind together. Interactive molecular docking applications often use haptic devices as a method of controlling the docking process. These devices allow the user to easily manipulate the structures in 3D space, whilst feeling the forces that occur in response to their manipulations. As a result of the force refresh rate requirements of haptic devices, haptic assisted docking applications are often limited, in that they model the interacting proteins as rigid, use low fidelity visualisations or require expensive propriety equipment to use. The research in this thesis aims to address some of these limitations. Firstly, the development of a visualisation algorithm capable of rendering a depiction of a deforming protein at an interactive refresh rate, with per-pixel shadows and ambient occlusion, is discussed. Then, a novel approach to modelling molecular flexibility whilst maintaining a stable haptic refresh rate is developed. Together these algorithms are presented within Haptimol FlexiDock, the first haptic-assisted molecular docking application to support receptor flexibility with high fidelity graphics, whilst also maintaining interactive refresh rates on both the haptic device and visual display. Using Haptimol FlexiDock, docking experiments were performed between two protein-ligand pairs: Maltodextrin Binding Protein and Maltose, and glutamine Binding Protein and Glucose. When the ligand was placed in its approximate binding site, the direction of over 80% of the intra-molecular movement aligned with that seen in the experimental structures. Furthermore, over 50% of the expected backbone motion was present in the structures generated with FlexiDock. Calculating the deformation of a biomolecule in real time, whilst maintaining an interactive refresh rate on the haptic device (> 500Hz) is a breakthrough in the field of interactive molecular docking, as, previous approaches either model protein flexibility, but fail to achieve the required haptic refresh rate, or do not consider biomolecular flexibility at all