Interactive molecular dynamics in virtual reality for accurate flexible protein-ligand docking

Molecular Docking Improved with Human Spatial Perception Using Virtual Reality

Adaptive steered molecular dynamics (ASMD) is a computational biophysics method in which an external force is applied to a selected set of atoms or a specific reaction coordinate to induce a particular molecular motion. Virtual reality (VR) based methods for protein-ligand docking are beneficial for visualizing on-the-fly interactive molecular dynamics and performing promising docking trajectories. In this paper, we propose a novel method to guide ASMD with optimal trajectories collected from human experiences using interactive molecular dynamics in virtual reality (iMD-VR). We also explain the benefits of using VR as a tool for expediting the process of ligand binding, outlining an experimental protocol that enables iMD-VR users to guide Amprenavir into and out of the binding pockets of HIV-1 protease and recreate their respective crystallographic binding poses within 5 minutes. Later, we discuss our analysis of the results from iMD-VR-assisted ASMD simulation and assess its performance compared to a standard ASMD simulation. From the accuracy point of view, our proposed method calculates higher Potential Mean Force (PMF) values consistently relative to a standard ASMD simulation with an almost twofold increase in all the experiments. Finally, we describe the novelty of the research and discuss results showcasing a faster and more effective convergence of the ligand to the protein's binding site as compared to a standard molecular dynamics simulation, proving the effectiveness of VR in the field of drug discovery. Future work includes the development of an artificial intelligence algorithm capable of predicting optimal binding trajectories for many protein-ligand pairs, as well as the required force needed to steer the ligand to follow the said trajectory.

Interactive molecular dynamics in virtual reality for modelling materials and catalysts

Interactive molecular dynamics simulation in virtual reality (iMD-VR) is emerging as a promising technique in molecular science. Here, we demonstrate its use in a range of fifteen applications in materials science and heterogeneous catalysis. In this work, the iMD-VR package Narupa is used with the MD package, DL_POLY [1]. We show how iMD-VR can be used to: (i) investigate the mechanism of lithium fast ion conduction by directing the formation of defects showing that vacancy transport is favoured over interstitialcy mechanisms, and (ii) guide a molecule through a zeolite pore to explore diffusion within zeolites, examining in detail the motion of methyl n-hexanoate in H-ZSM-5 zeolite and identifying bottlenecks restricting diffusion. iMD-VR allows users to manipulate these systems intuitively, to drive changes in them and observe the resulting changes in structure and dynamics. We make these simulations available, as a resource for both teaching and research. All simulation files, with videos, can be found online (https://doi.org/10.5281/zenodo.8252314) and are provided as open-source material.

Free energy along drug-protein binding pathways interactively sampled in virtual reality

We describe a two-step approach for combining interactive molecular dynamics in virtual reality (iMD-VR) with free energy (FE) calculation to explore the dynamics of biological processes at the molecular level. We refer to this combined approach as iMD-VR-FE. Stage one involves using a state-of-the-art 'human-in-the-loop' iMD-VR framework to generate a diverse range of protein-ligand unbinding pathways, benefitting from the sophistication of human spatial and chemical intuition. Stage two involves using the iMD-VR-sampled pathways as initial guesses for defining a path-based reaction coordinate from which we can obtain a corresponding free energy profile using FE methods. To investigate the performance of the method, we apply iMD-VR-FE to investigate the unbinding of a benzamidine ligand from a trypsin protein. The binding free energy calculated using iMD-VR-FE is similar for each pathway, indicating internal consistency. Moreover, the resulting free energy profiles can distinguish energetic differences between pathways corresponding to various protein-ligand conformations (e.g., helping to identify pathways that are more favourable) and enable identification of metastable states along the pathways. The two-step iMD-VR-FE approach offers an intuitive way for researchers to test hypotheses for candidate pathways in biomolecular systems, quickly obtaining both qualitative and quantitative insight.

W254 in furin functions as a molecular gate promoting anti-viral drug binding: Elucidation of putative drug tunneling and docking by non-equilibrium molecular dynamics

Furins are serine endoproteases that process precursor proteins into their biologically active forms, and they play essential roles in normal metabolism and disease presentation, including promoting expression of bacterial virulence factors and viral pathogenesis. Thus, furins represent vital targets for development of antimicrobial and antiviral therapeutics. Recent experimental evidence indicated that dichlorophenyl (DCP)-pyridine "BOS" drugs (e.g., BOS-318) competitively inhibit human furin by an induced-fit mechanism in which tryptophan W254 in the furin catalytic cleft (FCC) functions as a molecular gate, rotating nearly 180o through a steep energy barrier about its chi-1 dihedral to an "open" orientation, exposing a buried (i.e., cryptic) hydrophobic pocket 1. Once exposed, the non-polar DCP group of BOS-318, and similar halo-phenyl groups of analogs, enter the cryptic pocket, stabilizing drug binding. Here, we demonstrate flexible-receptor docking of BOS-318 (and various analogs) was unable to emulate the induced-fit motif, even when tryptophan was replaced with less bulky phenylalanine or glycine. While either substitution allowed access to the hydrophobic pocket for most ligands tested, optimal binding was observed only for W254, inferring a stabilizing effect of the indole sidechain. Furthermore, non-equilibrium steered molecular dynamics (sMD) in which the bound drugs (or their fragments) were extracted from the FCC did not cause closure of the open W254 gate, consistent with the thermodynamic stability of the open or closed W254 orientations. Finally, interactive molecular dynamics (iMD) revealed two putative conduits of drug entry and binding into the FCC, each coupled with W254 dihedral rotation and opening of the cryptic pocket. The iMD simulations further revealed ligand entry and binding in the FCC is likely driven in part by energy fluxes stemming from disruption and re-formation of ligand and protein solvation shells during drug migration from the solution phase into the FCC.

Co-creation environment with cloud virtual reality and real-time artificial intelligence toward the design of molecular robots

This paper describes the design philosophy for our cloud-based virtual reality (VR) co-creation environment (CCE) for molecular modeling. Using interactive VR simulation can provide enhanced perspectives in molecular modeling for intuitive live demonstration and experimentation in the CCE. Then the use of the CCE can enhance knowledge creation by bringing people together to share and create ideas or knowledge that may not emerge otherwise. Our prototype CCE discussed here, which was developed to demonstrate our design philosophy, has already enabled multiple members to log in and touch virtual molecules running on a cloud server with no noticeable network latency via real-time artificial intelligence techniques. The CCE plays an essential role in the rational design of molecular robot parts, which consist of bio-molecules such as DNA and protein molecules.

Open-Source Browser-Based Tools for Structure-Based Computer-Aided Drug Discovery

We here outline the importance of open-source, accessible tools for computer-aided drug discovery (CADD). We begin with a discussion of drug discovery in general to provide context for a subsequent discussion of structure-based CADD applied to small-molecule ligand discovery. Next, we identify usability challenges common to many open-source CADD tools. To address these challenges, we propose a browser-based approach to CADD tool deployment in which CADD calculations run in modern web browsers on users’ local computers. The browser app approach eliminates the need for user-initiated download and installation, ensures broad operating system compatibility, enables easy updates, and provides a user-friendly graphical user interface. Unlike server apps—which run calculations “in the cloud” rather than on users’ local computers—browser apps do not require users to upload proprietary information to a third-party (remote) server. They also eliminate the need for the difficult-to-maintain computer infrastructure required to run user-initiated calculations remotely. We conclude by describing some CADD browser apps developed in our lab, which illustrate the utility of this approach. Aside from introducing readers to these specific tools, we are hopeful that this review highlights the need for additional browser-compatible, user-friendly CADD software.

The emerging potential of interactive virtual reality in drug discovery

INTRODUCTION

The potential of virtual reality (VR) to contribute to drug design and development has been recognized for many years. A recent advance is to use VR not only to visualize and interact with molecules, but also to interact with molecular dynamics simulations 'on the fly' (interactive molecular dynamics in VR, IMD-VR), which is useful for flexible docking and examining binding processes and conformational changes.

AREAS COVERED

The authors use the term 'interactive VR' to refer to software where interactivity is an inherent part of the user VR experience e.g. in making structural modifications or interacting with a physically rigorous molecular dynamics (MD) simulation, as opposed to using VR controllers to rotate and translate the molecule for enhanced visualization. Here, they describe these methods and their application to problems relevant to drug discovery, highlighting the possibilities that they offer in this arena.

EXPERT OPINION

The ease of viewing and manipulating molecular structures and dynamics, using accessible VR hardware, and the ability to modify structures on the fly (e.g. adding or deleting atoms) - and for groups of researchers to work together in the same virtual environment - makes modern interactive VR a valuable tool to add to the armory of drug design and development methods.

Group VR experiences can produce ego attenuation and connectedness comparable to psychedelics

With a growing body of research highlighting the therapeutic potential of experiential phenomenology which diminishes egoic identity and increases one’s sense of connectedness, there is significant interest in how to elicit such ‘self-transcendent experiences’ (STEs) in laboratory contexts. Psychedelic drugs (YDs) have proven particularly effective in this respect, producing subjective phenomenology which reliably elicits intense STEs. With virtual reality (VR) emerging as a powerful tool for constructing new perceptual environments, we describe a VR framework called ‘Isness-distributed’ (Isness-D) which harnesses the unique affordances of distributed multi-person VR to blur conventional self-other boundaries. Within Isness-D, groups of participants co-habit a shared virtual space, collectively experiencing their bodies as luminous energetic essences with diffuse spatial boundaries. It enables moments of ‘energetic coalescence’, a new class of embodied intersubjective experience where bodies can fluidly merge, enabling participants to include multiple others within their self-representation. To evaluate Isness-D, we adopted a citizen science approach, coordinating an international network of Isness-D 'nodes'. We analyzed the results (N = 58) using 4 different self-report scales previously applied to analyze subjective YD phenomenology (the inclusion of community in self scale, ego-dissolution inventory, communitas scale, and the MEQ30 mystical experience questionnaire). Despite the complexities associated with a distributed experiment like this, the Isness-D scores on all 4 scales were statistically indistinguishable from recently published YD studies, demonstrating that distributed VR can be used to design intersubjective STEs where people dissolve their sense of self in the connection to others.

Online tools to easily build virtual molecular models for display in augmented and virtual reality on the web

Several groups developed in the last years augmented and virtual reality (AR/VR) software to visualize 3D molecules, most rather static, limited in content, and requiring software installs, some even requiring expensive hardware. We launched in 2020 moleculARweb (https://molecularweb.epfl.ch), a website that offers interactive content for chemistry and structural biology education through commodity web-based AR that works on consumer devices like smartphones, tablets and laptops. Among thousands of users, teachers increasingly request more biological macromolecules to be available, a demand that we cannot address individually. Therefore, to allow users to build their own material, we built a web interface where they can create online AR experiences in few steps starting from Protein Data Bank, AlphaFold or custom uploaded structures, or from virtual objects/scenes exported from the Visual Molecular Dynamics program, without any programming knowledge. The web tool also returns WebXR sessions for viewing and manipulating the models in WebXR-compatible devices including smartphones, tablets, and also immersive VR headsets with WebXR-capable browsers, where models can be manipulated even with bare hands when supported by the device. The tool is accessible for free at https://molecularweb.epfl.ch/pages/pdb2ar.html.

UDock2: interactive real-time multi-body protein-protein docking software

MOTIVATION

Protein-protein docking aims at predicting the geometry of protein interactions to gain insights into the mechanisms underlying these processes and develop new strategies for drug discovery. Interactive and user-oriented manipulation tools can support this task complementary to automated software.

RESULTS

This article presents an interactive multi-body protein-protein docking software, UDock2, designed for research but also usable for teaching and popularization of science purposes due to its high usability. In UDock2, the users tackle the conformational space of protein interfaces using an intuitive real-time docking procedure with on-the-fly scoring. UDock2 integrates traditional computer graphics methods to facilitate the visualization and to provide better insight into protein surfaces, interfaces, and properties.

AVAILABILITY AND IMPLEMENTATION

UDock2 is open-source, cross-platform (Windows and Linux), and available at http://udock.fr. The code can be accessed at https://gitlab.com/Udock/Udock2.

XR game development as a tool for authentic, experiential, and collaborative learning

This brief review provides insights into the recent advances made by the Unity game engine into education and the ways that learning resources can be used to teach novice bioscience students to use the VR/AR affordances of the platform. We present our own practises for implementation of a novel gamification task for final-year bioscience students and reflect on the changes we have had to make in order to overcome challenges presented by the COVID-19 pandemic.

Discovery of SARS-CoV-2 Mpro peptide inhibitors from modelling substrate and ligand binding

The main protease (Mpro) of SARS-CoV-2 is central to viral maturation and is a promising drug target, but little is known about structural aspects of how it binds to its 11 natural cleavage sites. We used biophysical and crystallographic data and an array of biomolecular simulation techniques, including automated docking, molecular dynamics (MD) and interactive MD in virtual reality, QM/MM, and linear-scaling DFT, to investigate the molecular features underlying recognition of the natural Mpro substrates. We extensively analysed the subsite interactions of modelled 11-residue cleavage site peptides, crystallographic ligands, and docked COVID Moonshot-designed covalent inhibitors. Our modelling studies reveal remarkable consistency in the hydrogen bonding patterns of the natural Mpro substrates, particularly on the N-terminal side of the scissile bond. They highlight the critical role of interactions beyond the immediate active site in recognition and catalysis, in particular plasticity at the S2 site. Building on our initial Mpro-substrate models, we used predictive saturation variation scanning (PreSaVS) to design peptides with improved affinity. Non-denaturing mass spectrometry and other biophysical analyses confirm these new and effective 'peptibitors' inhibit Mpro competitively. Our combined results provide new insights and highlight opportunities for the development of Mpro inhibitors as anti-COVID-19 drugs.

Exploring human-guided strategies for reaction network exploration: Interactive molecular dynamics in virtual reality as a tool for citizen scientists

The emerging fields of citizen science and gamification reformulate scientific problems as games or puzzles to be solved. Through engaging the wider non-scientific community, significant breakthroughs may be made by analyzing citizen-gathered data. In parallel, recent advances in virtual reality (VR) technology are increasingly being used within a scientific context and the burgeoning field of interactive molecular dynamics in VR (iMD-VR) allows users to interact with dynamical chemistry simulations in real time. Here, we demonstrate the utility of iMD-VR as a medium for gamification of chemistry research tasks. An iMD-VR "game" was designed to encourage users to explore the reactivity of a particular chemical system, and a cohort of 18 participants was recruited to playtest this game as part of a user study. The reaction game encouraged users to experiment with making chemical reactions between a propyne molecule and an OH radical, and "molecular snapshots" from each game session were then compiled and used to map out reaction pathways. The reaction network generated by users was compared to existing literature networks demonstrating that users in VR capture almost all the important reaction pathways. Further comparisons between humans and an algorithmic method for guiding molecular dynamics show that through using citizen science to explore these kinds of chemical problems, new approaches and strategies start to emerge.

Protein-protein interactions at a glance: Protocols for the visualization of biomolecular interactions

Protein-protein interactions (PPIs) play a key role in many biological processes and are intriguing targets for drug discovery campaigns. Advancements in experimental and computational techniques are leading to a growth of data accessibility, and, with it, an increased need for the analysis of PPIs. In this respect, visualization tools are essential instruments to represent and analyze biomolecular interactions. In this chapter, we reviewed some of the available tools, highlighting their features, and describing their functions with practical information on their usage.

Current trends in computer aided drug design and a highlight of drugs discovered via computational techniques: A review

Computer-aided drug design (CADD) is one of the pivotal approaches to contemporary pre-clinical drug discovery, and various computational techniques and software programs are typically used in combination, in a bid to achieve the desired outcome. Several approved drugs have been developed with the aid of CADD. On SciFinder®, we evaluated more than 600 publications through systematic searching and refining, using the terms, virtual screening; software methods; computational studies and publication year, in order to obtain data concerning particular aspects of CADD. The primary focus of this review was on the databases screened, virtual screening and/or molecular docking software program used. Furthermore, we evaluated the studies that subsequently performed molecular dynamics (MD) simulations and we reviewed the software programs applied, the application of density functional theory (DFT) calculations and experimental assays. To represent the latest trends, the most recent data obtained was between 2015 and 2020, consequently the most frequently employed techniques and software programs were recorded. Among these, the ZINC database was the most widely preferred with an average use of 31.2%. Structure-based virtual screening (SBVS) was the most prominently used type of virtual screening and it accounted for an average of 57.6%, with AutoDock being the preferred virtual screening/molecular docking program with 41.8% usage. Following the screening process, 38.5% of the studies performed MD simulations to complement the virtual screening and GROMACS with 39.3% usage, was the popular MD software program. Among the computational techniques, DFT was the least applied whereby it only accounts for 0.02% average use. An average of 36.5% of the studies included reports on experimental evaluations following virtual screening. Ultimately, since the inception and application of CADD in pre-clinical drug discovery, more than 70 approved drugs have been discovered, and this number is steadily increasing over time.

VRmol: an integrative web-based virtual reality system to explore macromolecular structure

SUMMARY

Structural visualization and analysis are fundamental to explore macromolecular functions. Here, we present a novel integrative web-based virtual reality (VR) system-VRmol, to visualize and study molecular structures in an immersive virtual environment. Importantly, it is integrated with multiple online databases and is able to couple structure studies with associated genomic variations and drug information in a visual interface by cloud-based drug docking. VRmol thus can serve as an integrative platform to aid structure-based translational research and drug design.

AVAILABILITY AND IMPLEMENTATION

VRmol is freely available (https://VRmol.net), with detailed manual and tutorial (https://VRmol.net/docs). The code of VRmol is available as open source under the MIT license at http://github.com/kuixu/VRmol.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Docking Paradigm in Drug Design

Docking is in demand for the rational computer aided structure based drug design. A review of docking methods and programs is presented. Different types of docking programs are described. They include docking of non-covalent small ligands, protein-protein docking, supercomputer docking, quantum docking, the new generation of docking programs and the application of docking for covalent inhibitors discovery. Taking into account the threat of COVID-19, we present here a short review of docking applications to the discovery of inhibitors of SARS-CoV and SARS-CoV-2 target proteins, including our own result of the search for inhibitors of SARS-CoV-2 main protease using docking and quantum chemical post-processing. The conclusion is made that docking is extremely important in the fight against COVID-19 during the process of development of antivirus drugs having a direct action on SARS-CoV-2 target proteins.

Interactive Molecular Dynamics in Virtual Reality Is an Effective Tool for Flexible Substrate and Inhibitor Docking to the SARS-CoV-2 Main Protease

The main protease (Mpro) of the SARS-CoV-2 virus is one focus of drug development efforts for COVID-19. Here, we show that interactive molecular dynamics in virtual reality (iMD-VR) is a useful and effective tool for creating Mpro complexes. We make these tools and models freely available. iMD-VR provides an immersive environment in which users can interact with MD simulations and so build protein complexes in a physically rigorous and flexible way. Recently, we have demonstrated that iMD-VR is an effective method for interactive, flexible docking of small molecule drugs into their protein targets (Deeks et al. PLoS One 2020, 15, e0228461). Here, we apply this approach to both an Mpro inhibitor and an oligopeptide substrate, using experimentally determined crystal structures. For the oligopeptide, we test against a crystallographic structure of the original SARS Mpro. Docking with iMD-VR gives models in agreement with experimentally observed (crystal) structures. The docked structures are also tested in MD simulations and found to be stable. Different protocols for iMD-VR docking are explored, e.g., with and without restraints on protein backbone, and we provide recommendations for its use. We find that it is important for the user to focus on forming binding interactions, such as hydrogen bonds, and not to rely on using simple metrics (such as RMSD), in order to create realistic, stable complexes. We also test the use of apo (uncomplexed) crystal structures for docking and find that they can give good results. This is because of the flexibility and dynamic response allowed by the physically rigorous, atomically detailed simulation approach of iMD-VR. We make our models (and interactive simulations) freely available. The software framework that we use, Narupa, is open source, and uses commodity VR hardware, so these tools are readily accessible to the wider research community working on Mpro (and other COVID-19 targets). These should be widely useful in drug development, in education applications, e.g., on viral enzyme structure and function, and in scientific communication more generally.

Combining Virtual Reality Visualization with Ensemble Molecular Dynamics to Study Complex Protein Conformational Changes

Molecular dynamics (MD) simulations are increasingly used to elucidate relationships between protein structure, dynamics, and their biological function. Currently, it is extremely challenging to perform MD simulations of large-scale structural rearrangements in proteins that occur on millisecond timescales or beyond, as this requires very significant computational resources, or the use of cumbersome "collective variable" enhanced sampling protocols. Here, we describe a framework that combines ensemble MD simulations and virtual reality visualization (eMD-VR) to enable users to interactively generate realistic descriptions of large amplitude, millisecond timescale protein conformational changes in proteins. Detailed tests demonstrate that eMD-VR substantially decreases the computational cost of folding simulations of a WW domain, without the need to define collective variables a priori. We further show that eMD-VR generated pathways can be combined with Markov state models to describe the thermodynamics and kinetics of large-scale loop motions in the enzyme cyclophilin A. Our results suggest eMD-VR is a powerful tool for exploring protein energy landscapes in bioengineering efforts.