Summary
O Visual Dynamics é uma ferramenta de código aberto que acelera implementações e aprendizado em simulação de dinâmica molecular usando Gromacs. O protocolo apresentado irá guiá-lo através das etapas para realizar uma simulação proteína-ligante preparada no ACPYPE com facilidade e etapas gerais para outros modelos de simulação.
Abstract
A Dinâmica Visual (VD) é uma ferramenta web que visa facilitar o uso e a aplicação da Dinâmica Molecular (MD) executada no Gromacs, permitindo que usuários sem familiaridade computacional executem simulações de curta duração para fins de validação, demonstração e ensino. É verdade que os métodos quânticos são os mais precisos. No entanto, atualmente não há viabilidade computacional para realizar os experimentos que o MD realiza. A ferramenta descrita aqui tem recebido melhorias contínuas ao longo dos últimos dois anos. Este protocolo descreverá o que é necessário para executar uma simulação em VD com um complexo proteína-ligante previamente preparado em ACPYPE e algumas direções gerais sobre os outros modelos de simulação disponíveis. Para a simulação detalhada, será utilizada a proteína ligadora de FK506 de Plasmodium vivax complexada com o inibidor D5 (PDB ID: 4mgv), e todos os arquivos utilizados serão fornecidos. Observe que este protocolo dirá todas as opções a serem usadas para alcançar os mesmos resultados apresentados, mas essas opções não são necessariamente as únicas disponíveis.
Introduction
De acordo com a definição da IUPAC, MD é o procedimento de simulação que consiste em calcular o movimento de átomos em uma molécula ou de átomos individuais ou moléculas em sólidos, líquidos e gases, de acordo com as leis do movimento de Newton. As forças que atuam sobre os átomos, necessárias para simular seu movimento, são comumente calculadas usando campos de força da mecânica molecular1. Pode ser aplicado a qualquer fenômeno que busque extrair informações em nível molecular e muitas vezes atômico2.
A DM é uma das técnicas incorporadas à bioinformática, especificamente à bioinformática estrutural. Com ele, é possível obter características cinéticas e termodinâmicas de estruturas biomoleculares. Por exemplo, estabilidade macromolecular, identificação de sítios alostéricos, elucidação de mecanismos de atividade enzimática, reconhecimento molecular e propriedades de complexos com pequenas moléculas, associação entre proteínas, enovelamento de proteínas e sua hidratação3. Além disso, a DM permite uma ampla gama de estudos, incluindo o design molecular (amplamente utilizado no design de medicamentos), na determinação da estrutura e seu refinamento (raios-X, RMN e modelagem de proteínas)3. Os resultados obtidos ao final de um MD são os mais ricos e completos em termos de simulação não quântica4. A DM clássica é muito mais eficiente do que se poderia esperar de uma consideração completa da física dos sistemas biomoleculares devido ao número de aproximações substanciais. Notavelmente, os efeitos dinâmicos quânticos são geralmente ignorados3. No entanto, implementar um experimento de MD não é trivial5. Requer conhecimento de computação, especialmente o Terminal Linux, já que a maioria dos softwares de bioinformática estrutural é feita para isso. Mesmo com esse conhecimento, aprender os comandos e a parametrização do Gromacs é outra curva de aprendizado íngreme.
Desde sua primeira aplicação à biologia em 19776, muito evoluiu devido ao aumento do processamento computacional e à codificação aprimorada. Há mais de duas décadas, foi lançado o primeiro software MD destinado a problemas biológicos, ou seja, Gromacs7, AMBER8 e NAMD9.
Desde suas primeiras versões, esses softwares ainda continuam sendo os mais usados e citados. No entanto, eles continuam com as mesmas dificuldades comuns de implementação que afligem os pesquisadores que não são especialistas em computação5. Alguns têm etapas complexas de instalação e configuração, às vezes exigindo amplo conhecimento sobre o hardware em que serão executados para obter o máximo dele e documentação técnica altamente centrada no computador. É necessária uma maneira mais fácil de interagir com eles, além da linha de comando e dos parâmetros infinitos.
Uma interface atua como intermediária entre o processo lógico a ser realizado e o humano10. O paradigma de como o software é executado evoluiu à medida que os recursos de computação melhoraram. O primeiro paradigma digital foram as interfaces de linha de comando (CLI), seguidas pela evolução para as interfaces gráficas de usuário (GUI) conhecidas11. Seguindo o ciclo evolutivo, a interface produzida pela World Wide Web (ou simplesmente WEB) é considerada uma evolução das GUIs11. Esses três paradigmas atualmente coexistem dependendo dos desenvolvedores. Os aplicativos CLI usam comandos textuais no console do sistema operacional. Os aplicativos GUI, também chamados de áreas de trabalho gráficas, usam uma interface gráfica composta de janelas, botões e outros componentes. É específico e pré-programado para um sistema operacional. A principal diferença da CLI é o uso do mouse como elemento adicional na interação homem-máquina12. As aplicações WEB, apesar de serem confundidas com uma GUI, são mais complexas de desenvolver, mas são mais versáteis e de longe as mais ágeis em operação. Além disso, dependem apenas de um software intérprete denominado navegador, que possibilita que a aplicação cliente se comunique com o servidor por meio de uma rede independente do sistema operacional13.
O software de bioinformática estrutural geralmente usa paradigmas CLI e GUI. Alguns exemplos de software clássico que usam CLI são o Modeller14 para modelagem de similaridade, o Autodock15 para acoplamento molecular e o Gromacs16 para dinâmica molecular. Exemplos de softwares que adotam o tipo GUI são SwissPDBviewer17, Pymol18, VMD19, UCSF Chimera20, Autodock tools15, PyRx21, Biovia22, Maestro23 e Moe24, entre outros.
Com o surgimento das tecnologias Hypertext Markup Language versão 5 (HTML5)25, Cascading Style Sheets (CSS)26 e Javascript27 , entre outras, muitas aplicações de bioinformática estrutural puderam ser trazidas para a WEB, tornando-se assim mais acessíveis. Exemplos de modelagem de similaridade de servidores WEB são o MODWEB28, que usa o Modeller14 como back-end e o Swissmodel29. Exemplos de servidores de aplicativos da Web para docking molecular são Haddock30, Swissdock31, Cluspro32, Dockthor33 e outros.
Enquanto as metodologias de análise estrutural, modelagem e encaixe evoluíram dos paradigmas CLI para GUI e, finalmente, para WEB, o MD continua a ser suportado principalmente pela execução de linha de comando (tipo CLI). Algumas boas iniciativas têm surgido para melhorar esse panorama. Exemplos dessas iniciativas são a implementação de plug-ins em softwares existentes, como o plug-in QwikMD para VMD34, o plug-in GROMACS para PyMOL e a opção Molecular Dynamics Simulation no UCSF Chimera20, alguns aplicativos CLI novos e mais fáceis, como ASGARD35, Gmx_qk36 e CHAPERONg37, e uma plataforma web robusta, BioBB-Wfs38. Embora o uso desses plugins e aplicativos seja um avanço, sua implementação ainda é um desafio para a maioria dos pesquisadores não qualificados. Dificuldades comuns incluem problemas na instalação e configuração do software MD, que muitas vezes comprometem a execução completa da simulação5.
Em 2022, o software Visual Dynamics para simulação computacional baseada na web foi disponibilizado pelo Laboratório de Bioinformática e Química Medicinal da Fiocruz Rondônia39. Sua versão inicial foi construída em Python e Flask, permitindo simulações de sistemas com proteínas livres (apoenzimas) por apenas 2 ns. Posteriormente, foi aprimorado para incluir uma versão de simulação automatizada com ligantes preparados usando PRODRG40.
O VD foi construído para auxiliar todos os pesquisadores da área de biofísica estrutural, biotecnologia e áreas afins que possuem limitações no conhecimento computacional; a ferramenta permite que esses pesquisadores testem suas hipóteses envolvendo simulações de MD a partir de qualquer sistema operacional e sem acesso a um computador de alto desempenho (HPC). O objetivo deste trabalho é apresentar as novas funcionalidades do Visual Dynamics versão 3.0. Além disso, visa introduzir um protocolo de uso atualizado para a ferramenta e destacar as limitações a serem abordadas no futuro, juntamente com as estatísticas de uso até o presente momento (Figura 1).
Protocol
Representative Results
Discussion
Automatizar processos não é fácil, mas também é menos difícil do que reprogramar um sistema do zero. Gromacs é atualmente o software de simulação molecular mais popular e é constantemente atualizado. O Departamento de Química Biofísica da Universidade de Groningen inicialmente o desenvolveu, e agora é mantido pelo Laboratório de Ciências da Vida da Universidade de Estocolmo43.
Para qualquer novo usuário, aprender técnicas de simulação é uma longa jor…
Disclosures
The authors have nothing to disclose.
Acknowledgements
Este trabalho contou com o apoio da Fundação Oswaldo Cruz (Fiocruz), da Fundação para o Desenvolvimento Científico e Tecnológico em Saúde (Fiotec), do Instituto Nacional de Ciência e Tecnologia de Epidemiologia da Amazônia Ocidental – INCT-EpiAmO, da Fundação Rondônia de Amparo ao Desenvolvimento das Ações Científicas e Tecnológicas e à Pesquisa do Estado de Rondônia (FAPERO), a Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) e o 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
- . . The IUPAC Compendium of Chemical Terminology: The Gold Book. , (2019).
- Karplus, M., McCammon, J. A. Molecular dynamics simulations of biomolecules. Nat Str Biol. 9 (9), 646-652 (2002).
- Adcock, S. A., McCammon, J. A. Molecular dynamics: Survey of methods for simulating the activity of proteins. Chem Rev. 106 (5), 1589-1615 (2006).
- Stillinger, F. H., Rahman, A. Improved simulation of liquid water by molecular dynamics. J Chem Phys. 60 (4), 1545-1557 (1974).
- Zinovjev, K., Van Der Kamp, M. W. Enlighten2: molecular dynamics simulations of protein-ligand systems made accessible. Bioinformatics. 36 (20), 5104-5106 (2020).
- McCammon, J. A., Gelin, B. R., Karplus, M. Dynamics of folded proteins. Nature. 267 (5612), 585-590 (1977).
- 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).
- 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).
- Nelson, M. T. NAMD: A parallel, object-oriented molecular dynamics program. Int J Supercomp Appl High Performance Comp. 10 (4), 251-268 (1996).
- . . The Art of human-computer interface design. , (1990).
- Wigdor, D., Wixon, D. . Brave NUI World: Designing Natural User Interfaces for Touch and Gesture. , (2014).
- Unwin, A., Heike, H. . GUI and Command-line – Conflict or Synergy. , (2000).
- Clark, D. Developing Web Applications. Beginning C# Object-Oriented Programming. , 243-263 (2011).
- Eswar, N., et al. Comparative protein structure modeling using Modeller. Curr Protoc Bioinformatics. 5 (5.6), (2006).
- Morris, G. M. AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J Comp Chem. 30 (16), 2785-2791 (2009).
- 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).
- Guex, N., Peitsch, M. C. SWISS-MODEL and the Swiss-Pdb Viewer: An environment for comparative protein modeling. Electrophoresis. 18 (15), 2714-2723 (1997).
- . The PyMOL molecular graphics system Available from: https://pymol.org/ (2023)
- Humphrey, W., Dalke, A., Schulten, K. VMD: Visual molecular dynamics. J Mol Graph. 14 (1), 33-38 (1996).
- Pettersen, E. F., et al. UCSF Chimera-A visualization system for exploratory research and analysis. J Comp Chem. 25 (13), 1605-1612 (2004).
- . PyRx-python prescription v. 0.8 Available from: https://pyrx.sourceforge.io/ (2023)
- . PyRx-python prescription v. 0.8 Available from: https://pyrx.sourceforge.io/ (2023)
- . Maestro Available from: https://www.schrodinger.com/platform/products/maestro/ (2024)
- . Molecular operating environment (MOE) Available from: https://www.chemcomp.com/en/Products.htm (2023)
- . The Syntax, Vocabulary and APIs of HTML5 Available from: https://dev.w3.org/html5/html-author/ (2023)
- . Cascading Style Sheets Available from: https://www.w3.org/Style/CSS/ (2023)
- . ECMAScript 2020 language specification, 11th edition Available from: https://www.ecma-international.org/publications-and-standards/standards/ecma-262/ (2023)
- Pieper, U., et al. ModBase, a database of annotated comparative protein structure models, and associated resources. Nucleic Acids Res. 39 (Database), D465-D474 (2011).
- 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).
- 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).
- 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).
- Kozakov, D., et al. The ClusPro web server for protein-protein docking. Nat Prot. 12 (2), 255-278 (2017).
- 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).
- Ribeiro, J. V., et al. QwikMD – Integrative molecular dynamics toolkit for novices and experts. Sci Rep. 6 (1), 26536 (2016).
- Rodríguez Martínez, A., et al. . ASGARD. A simple and automatic GROMACS tool to analyze Molecular Dynamic simulations. , (2023).
- 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).
- 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).
- 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).
- Vieira, I. H. P., et al. Visual dynamics: a WEB application for molecular dynamics simulation using GROMACS. BMC Bioinfo. 24 (1), 107 (2023).
- 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).
- 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).
- Sousa Da Silva, A. W., Vranken, W. F. ACPYPE – AnteChamber PYthon Parser interfacE. BMC Res Notes. 5 (1), 367 (2012).
- . Gromacs developmet and developers Available from: https://www.gromacs.org/development.html (2023)
- Andrio, P., et al. BioExcel building blocks, a software library for interoperable biomolecular simulation workflows. Sci Data. 6 (1), 169 (2019).
- Páll, S., et al. Heterogeneous parallelization and acceleration of molecular dynamics simulations in GROMACS. J Chem Phys. 153 (13), 134110 (2020).
