Summary

Новые возможности Visual Dynamics 3.0

Published: August 09, 2024
doi:

Summary

Visual Dynamics is an open-source tool that accelerates implementations and learning in molecular dynamics simulation using Gromacs. The presented protocol will guide you through the steps to perform a protein-ligand simulation prepared in ACPYPE with ease and general steps to other simulation models.

Abstract

Visual Dynamics (VD) is a web tool that aims to facilitate the use and application of Molecular Dynamics (MD) executed in Gromacs, allowing users without computational familiarity to run short-time simulations for validation, demonstration, and teaching purposes. It is true that quantum methods are the most accurate. However, there is currently no computational feasibility to carry out the experiments that MD performs. The tool described here has continuously received improvements over the course of the last couple of years. This protocol will describe what is needed to run a simulation in VD with a protein-ligand complex previously prepared in ACPYPE and some general directions on the other simulation models available. For the detailed simulation, the FK506-binding protein from Plasmodium vivax complexed with the inhibitor D5 (PDB ID: 4mgv) will be used, and all files used will be provided. Note that this protocol will tell every option to be used to achieve the same results presented, but these options are not necessarily the only ones available.

Introduction

According to the IUPAC definition, MD is the simulation procedure that consists of computing the motion of atoms in a molecule or of individual atoms or molecules in solids, liquids, and gases, according to Newton's laws of motion. The forces acting on atoms, necessary to simulate their motion, are commonly calculated using force fields from molecular mechanics1. It can be applied to any phenomenon that seeks to extract information at a molecular and often atomic level2.

MD is one of the techniques incorporated into bioinformatics, specifically structural bioinformatics. With it, it is possible to obtain kinetic and thermodynamic characteristics of biomolecular structures. For example, macromolecular stability, identification of allosteric sites, elucidation of mechanisms of enzymatic activity, molecular recognition and properties of complexes with small molecules, association between proteins, protein folding, and its hydration3. Furthermore, MD enables a wide range of studies, including molecular design (widely used in drug design), in determining the structure and its refinement (X-ray, NMR, and protein modeling)3. The results obtained at the end of an MD are the richest and most complete in terms of non-quantum simulation4.Classical MD is much more efficient than might be expected from a full consideration of the physics of biomolecular systems due to the number of substantial approximations. Notably, quantum dynamical effects are usually ignored3. However, implementing an MD experiment is not trivial5. It requires knowledge of computing, especially the Linux Terminal, as most structural bioinformatics software is made for it. Even with that knowledge, learning Gromacs commands and parametrization is another steep learning curve.

Since its first application to biology in 19776, much has evolved due to increased computational processing and improved coding. More than two decades ago, the first MD software intended for biological problems was launched, namely Gromacs7, AMBER8, and NAMD9.

Since their first versions, these software still remain the most used and cited. However, they continue with the same common implementation difficulties that plague researchers who are not computer specialists5. Some have complex installation and configuration steps, sometimes requiring extensive knowledge about the hardware it will run on to get the most out of it and highly computer-centric technical documentation. An easier way to interface with them, aside from the command line and infinite parameters, is needed.

An interface acts as an intermediary between the logical process to be performed and the human10. The paradigm of how software is executed has evolved as computing resources have improved. The first digital paradigm was the command line interfaces (CLI) followed by the evolution to the known graphical user interfaces (GUI)11. Following the evolutionary cycle, the interface produced by the World Wide Web (or simply WEB) is considered an evolution of GUIs11. These three paradigms currently co-exist depending on developers. CLI applications use textual commands on the operating system console. GUI applications, also called graphical desktops, use a graphical interface made up of windows, buttons, and other components. It is specific and pre-programmed for an operating system. The main difference from the CLI is the use of the mouse as an additional element in human-machine interaction12. WEB applications, despite being confused with a GUI, are more complex to develop but are more versatile and by far the most agile in operation. Furthermore, they only depend on an interpreter software called a browser, which makes it possible for the client application to communicate with the server through a network independent of the operating system13.

Structural bioinformatics software most commonly use CLI and GUI paradigms. Some examples of classic software that use CLI are Modeller14 for similarity modeling, Autodock15 for molecular docking, and Gromacs16 for molecular dynamics. Examples of software that adopt the GUI type are SwissPDBviewer17, Pymol18, VMD19, UCSF Chimera20, Autodock tools15, PyRx21, Biovia22, Maestro23, and Moe24, among others.

With the emergence of Hypertext Markup Language version 5 (HTML5)25, Cascading Style Sheets (CSS)26, and Javascript27 technologies, among others, many structural bioinformatics applications could be brought to the WEB, thus becoming more accessible. Examples of similarity modeling WEB servers are MODWEB28, which uses Modeller14 as a back-end and Swissmodel29. Examples of web application servers for molecular docking are Haddock30, Swissdock31, Cluspro32, Dockthor33, and others.

While structural analysis, modeling, and docking methodologies evolved from CLI paradigms to GUI and finally to WEB, MD continues to be mostly supported by command line execution (CLI type). Some good initiatives have emerged to improve this panorama. Examples of these initiatives are the implementation of plugins in existing software, such as QwikMD plugin to VMD34, GROMACS Plugin to PyMOL, and the Molecular Dynamics Simulation option in UCSF Chimera20, some new and easier CLI applications, such as ASGARD35, Gmx_qk36, and CHAPERONg37, and a robust web platform, BioBB-Wfs38. Although the use of these plugins and applications is an advance, their implementation is still a challenge for most unskilled researchers. Common difficulties include problems installing and configuring the MD software, which often compromise the full execution of the simulation5.

In 2022, the Visual Dynamics software for web-based computational simulation was made available by the Laboratório de Bioinformática e Química Medicinal at Fiocruz Rondônia39. Its initial version was built in Python and Flask, allowing simulations of systems with free proteins (apoenzymes) for only 2 ns. Subsequently, it was enhanced to include an automated simulation version with ligands prepared using PRODRG40.

VD was built to assist all researchers in the field of structural biophysics, biotechnology, and related areas who have limitations in computational knowledge; the tool allows these researchers to test their hypotheses involving MD simulations from any operational system and without access to a high-performance computer (HPC). The purpose of this work is to present the new features of Visual Dynamics version 3.0. Additionally, it aims to introduce an updated usage protocol for the tool and highlight the limitations to be addressed in the future, along with usage statistics up to the present moment (Figure 1).

Protocol

1. Accessing the software and new user registration Visit the Visual Dynamics (VD) web page. Click on the +Register icon at the top right to create an account. Register to use the software. NOTE: Only institutional email addresses are allowed. The user will receive an email notification once their registration is approved. Click on Login at the top right to access the system login screen. Fill in the username/email and password fields and cli…

Representative Results

VD provides a fully autonomous simulation execution that does not require user intervention or user-provided computational resources. After submitting a simulation to execution, the user can leave it, turn off their machines, and the simulation will continue running. It also allows users to access the results from any device, be it a laptop or mobile device. As an example of using VD in automated mode through the WEB, the test was made for a protein-ligand complex prepared in ACPYPE using the …

Discussion

Automating processes is not easy, but it is also less difficult than reprogramming a system from scratch. Gromacs is currently the most popular molecular simulation software, and it is constantly updated. The Department of Biophysical Chemistry at Groningen University initially developed it, and it is now maintained by the Life Sciences Laboratory at the University of Stockholm43.

For any new user, learning simulation techniques is a lengthy journey. VD emerges as an al…

Disclosures

The authors have nothing to disclose.

Acknowledgements

This work has been supported by the The Fundação Oswaldo Cruz (Fiocruz), the Fundação para o Desenvolvimento Científico e Tecnológico em Saúde (Fiotec), the Instituto Nacional de Ciência e Tecnologia de Epidemiologia da Amazônia Ocidental – INCT-EpiAmO, the Fundação Rondônia de Amparo ao Desenvolvimento das Ações Científicas e Tecnológicas e à Pesquisa do Estado de Rondônia (FAPERO), the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) and the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq).

Materials

ACPYPE Server Bio2Byte Available at https://www.bio2byte.be/acpype/
GRACE software Plasma Laboratory at the Weizmann Institute of Science Available at https://plasma-gate.weizmann.ac.il/Grace/
GROMACS software GROMACS Team Installation instructions at https://manual.gromacs.org/current/install-guide/index.html
The structure of the FK506-binding protein
From Plasmodium vivax complexed with the
inhibitor D5
RCSB Protein Data Bank Available at https://www.rcsb.org/structure/4mgv
Already contains the ligand complexed to the macromolecule.

References

  1. . . The IUPAC Compendium of Chemical Terminology: The Gold Book. , (2019).
  2. Karplus, M., McCammon, J. A. Molecular dynamics simulations of biomolecules. Nat Str Biol. 9 (9), 646-652 (2002).
  3. Adcock, S. A., McCammon, J. A. Molecular dynamics: Survey of methods for simulating the activity of proteins. Chem Rev. 106 (5), 1589-1615 (2006).
  4. Stillinger, F. H., Rahman, A. Improved simulation of liquid water by molecular dynamics. J Chem Phys. 60 (4), 1545-1557 (1974).
  5. Zinovjev, K., Van Der Kamp, M. W. Enlighten2: molecular dynamics simulations of protein-ligand systems made accessible. Bioinformatics. 36 (20), 5104-5106 (2020).
  6. McCammon, J. A., Gelin, B. R., Karplus, M. Dynamics of folded proteins. Nature. 267 (5612), 585-590 (1977).
  7. Berendsen, H. J. C., Van Der Spoel, D., Van Drunen, R. GROMACS: A message-passing parallel molecular dynamics implementation. Comp Phys Comm. 91 (1-3), 43-56 (1995).
  8. Pearlman, D. A., et al. AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of molecules. Comp Phys Comm. 91 (1-3), 1-41 (1995).
  9. Nelson, M. T. NAMD: A parallel, object-oriented molecular dynamics program. Int J Supercomp Appl High Performance Comp. 10 (4), 251-268 (1996).
  10. . . The Art of human-computer interface design. , (1990).
  11. Wigdor, D., Wixon, D. . Brave NUI World: Designing Natural User Interfaces for Touch and Gesture. , (2014).
  12. Unwin, A., Heike, H. . GUI and Command-line – Conflict or Synergy. , (2000).
  13. Clark, D. Developing Web Applications. Beginning C# Object-Oriented Programming. , 243-263 (2011).
  14. Eswar, N., et al. Comparative protein structure modeling using Modeller. Curr Protoc Bioinformatics. 5 (5.6), (2006).
  15. Morris, G. M. AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J Comp Chem. 30 (16), 2785-2791 (2009).
  16. Van Der Spoel, D., Lindahl, E., Hess, B., Groenhof, G., Mark, A. E., Berendsen, H. J. C. GROMACS: Fast, flexible, and free. J Comp Chem. 26 (16), 1701-1718 (2005).
  17. Guex, N., Peitsch, M. C. SWISS-MODEL and the Swiss-Pdb Viewer: An environment for comparative protein modeling. Electrophoresis. 18 (15), 2714-2723 (1997).
  18. . The PyMOL molecular graphics system Available from: https://pymol.org/ (2023)
  19. Humphrey, W., Dalke, A., Schulten, K. VMD: Visual molecular dynamics. J Mol Graph. 14 (1), 33-38 (1996).
  20. Pettersen, E. F., et al. UCSF Chimera-A visualization system for exploratory research and analysis. J Comp Chem. 25 (13), 1605-1612 (2004).
  21. . PyRx-python prescription v. 0.8 Available from: https://pyrx.sourceforge.io/ (2023)
  22. . PyRx-python prescription v. 0.8 Available from: https://pyrx.sourceforge.io/ (2023)
  23. . Maestro Available from: https://www.schrodinger.com/platform/products/maestro/ (2024)
  24. . Molecular operating environment (MOE) Available from: https://www.chemcomp.com/en/Products.htm (2023)
  25. . The Syntax, Vocabulary and APIs of HTML5 Available from: https://dev.w3.org/html5/html-author/ (2023)
  26. . Cascading Style Sheets Available from: https://www.w3.org/Style/CSS/ (2023)
  27. . ECMAScript 2020 language specification, 11th edition Available from: https://www.ecma-international.org/publications-and-standards/standards/ecma-262/ (2023)
  28. Pieper, U., et al. ModBase, a database of annotated comparative protein structure models, and associated resources. Nucleic Acids Res. 39 (Database), D465-D474 (2011).
  29. Kiefer, F., Arnold, K., Kunzli, M., Bordoli, L., Schwede, T. The SWISS-MODEL Repository and associated resources. Nucleic Acids Res. 37 (Database), D387-D392 (2009).
  30. De Vries, S. J., Van Dijk, M., Bonvin, A. M. J. J. The HADDOCK web server for data-driven biomolecular docking. Nat Prot. 5 (5), 883-897 (2010).
  31. Grosdidier, A., Zoete, V., Michielin, O. SwissDock, a protein-small molecule docking web service based on EADock DSS. Nucleic Acids Res. 39 (suppl), W270-W277 (2011).
  32. Kozakov, D., et al. The ClusPro web server for protein-protein docking. Nat Prot. 12 (2), 255-278 (2017).
  33. Santos, K. B., Guedes, I. A., Karl, A. L. M., Dardenne, L. E. Highly flexible ligand docking: Benchmarking of the DockThor program on the LEADS-PEP protein-peptide data set. J Chem Info Modeling. 60 (2), 667-683 (2020).
  34. Ribeiro, J. V., et al. QwikMD – Integrative molecular dynamics toolkit for novices and experts. Sci Rep. 6 (1), 26536 (2016).
  35. Rodríguez Martínez, A., et al. . ASGARD. A simple and automatic GROMACS tool to analyze Molecular Dynamic simulations. , (2023).
  36. Singh, H., Raja, A., Prakash, A., Medhi, B. Gmx_qk: An automated protein/protein-ligand complex simulation workflow bridged to MM/PBSA, based on Gromacs and Zenity-Dependent GUI for beginners in MD simulation study. J Chem Info Modeling. 63 (9), 2603-2608 (2023).
  37. Yekeen, A. A., Durojaye, O. A., Idris, M. O., Muritala, H. F., Arise, R. O. CHAPERONg: A tool for automated GROMACS-based molecular dynamics simulations and trajectory analyses. Comp Str Biotechnol J. 21, 4849-4858 (2023).
  38. Bayarri, G., Andrio, P., Hospital, A., Orozco, M., Gelpí, J. L. BioExcel building blocks workflows (BioBB-Wfs), an integrated web-based platform for biomolecular simulations. Nucleic Acids Res. 50 (W1), (2022).
  39. Vieira, I. H. P., et al. Visual dynamics: a WEB application for molecular dynamics simulation using GROMACS. BMC Bioinfo. 24 (1), 107 (2023).
  40. Schüttelkopf, A. W., Van Aalten, D. M. F. PRODRG: a tool for high-throughput crystallography of protein-ligand complexes. Acta Crystallographica Sect D Biol Crystallography. 60 (8), 1355-1363 (2004).
  41. Harikishore, A., et al. Adamantyl derivative as a potent inhibitor of Plasmodium FK506 binding protein 35. ACS Med Chem Lett. 4 (11), 1097-1101 (2013).
  42. Sousa Da Silva, A. W., Vranken, W. F. ACPYPE – AnteChamber PYthon Parser interfacE. BMC Res Notes. 5 (1), 367 (2012).
  43. . Gromacs developmet and developers Available from: https://www.gromacs.org/development.html (2023)
  44. Andrio, P., et al. BioExcel building blocks, a software library for interoperable biomolecular simulation workflows. Sci Data. 6 (1), 169 (2019).
  45. Páll, S., et al. Heterogeneous parallelization and acceleration of molecular dynamics simulations in GROMACS. J Chem Phys. 153 (13), 134110 (2020).
This article has been published
Video Coming Soon
Keep me updated:

.

Cite This Article
Henrique Provensi Vieira, I., Mendonça, E. A. M., Guariero, F. L., Guimarães, R. M. d. S., Zanchi, F. B. New Features in Visual Dynamics 3.0. J. Vis. Exp. (210), e66964, doi:10.3791/66964 (2024).

View Video