Summary

Summary

Visual Dynamics ist ein Open-Source-Tool, das die Implementierung und das Lernen in der Molekulardynamik-Simulation mit Gromacs beschleunigt. Das vorgestellte Protokoll führt Sie mit Leichtigkeit durch die Schritte zur Durchführung einer Protein-Liganden-Simulation, die in ACPYPE vorbereitet wurde, und durch allgemeine Schritte zu anderen Simulationsmodellen.

Abstract

Visual Dynamics (VD) ist ein Web-Tool, das darauf abzielt, die Verwendung und Anwendung von Molecular Dynamics (MD) zu erleichtern, das in Gromacs ausgeführt wird, und es Benutzern ohne rechnerische Erfahrung ermöglicht, Kurzzeitsimulationen zu Validierungs-, Demonstrations- und Lehrzwecken durchzuführen. Es stimmt, dass Quantenmethoden am genauesten sind. Derzeit gibt es jedoch keine rechnerische Machbarkeit, um die von MD durchgeführten Experimente durchzuführen. Das hier beschriebene Tool wurde in den letzten Jahren kontinuierlich verbessert. Dieses Protokoll beschreibt, was erforderlich ist, um eine Simulation in VD mit einem Protein-Liganden-Komplex durchzuführen, der zuvor in ACPYPE hergestellt wurde, und einige allgemeine Anweisungen zu den anderen verfügbaren Simulationsmodellen. Für die detaillierte Simulation wird das FK506-bindende Protein aus Plasmodium vivax verwendet, das mit dem Inhibitor D5 (PDB-ID: 4mgv) komplexiert ist, und alle verwendeten Dateien werden zur Verfügung gestellt. Beachten Sie, dass dieses Protokoll angibt, dass jede Option verwendet werden soll, um die gleichen Ergebnisse zu erzielen, aber diese Optionen sind nicht unbedingt die einzigen, die verfügbar sind.

Introduction

Nach der IUPAC-Definition ist MD das Simulationsverfahren, das darin besteht, die Bewegung von Atomen in einem Molekül oder von einzelnen Atomen oder Molekülen in Festkörpern, Flüssigkeiten und Gasen nach den Newtonschen Bewegungsgesetzen zu berechnen. Die auf Atome wirkenden Kräfte, die notwendig sind, um ihre Bewegung zu simulieren, werden üblicherweise mit Hilfe von Kraftfeldern aus der Molekularmechanikberechnet 1. Es kann auf jedes Phänomen angewendet werden, das versucht, Informationen auf molekularer und oft atomarer Ebene zu extrahieren2.

MD ist eine der Techniken, die in die Bioinformatik integriert sind, insbesondere in die strukturelle Bioinformatik. Damit ist es möglich, kinetische und thermodynamische Eigenschaften biomolekularer Strukturen zu erhalten. Zum Beispiel die makromolekulare Stabilität, die Identifizierung allosterischer Stellen, die Aufklärung der Mechanismen der enzymatischen Aktivität, die molekulare Erkennung und die Eigenschaften von Komplexen mit kleinen Molekülen, die Assoziation zwischen Proteinen, die Proteinfaltung und ihre Hydratation3. Darüber hinaus ermöglicht MD eine breite Palette von Studien, einschließlich des molekularen Designs (das im Wirkstoffdesign weit verbreitet ist), zur Bestimmung der Struktur und ihrer Verfeinerung (Röntgen-, NMR- und Proteinmodellierung)3. Die Ergebnisse, die am Ende einer MD erzielt werden, sind die reichhaltigsten und vollständigsten in Bezug auf die Nicht-Quantensimulation4. Die klassische MD ist aufgrund der Vielzahl substantieller Approximationen viel effizienter, als man bei einer vollständigen Berücksichtigung der Physik biomolekularer Systeme erwarten könnte. Bemerkenswert ist, dass quantendynamische Effekte in der Regel ignoriert werden3. Die Implementierung eines MD-Experiments ist jedoch nicht trivial5. Es erfordert Kenntnisse in der Informatik, insbesondere im Linux-Terminal, da die meisten strukturellen Bioinformatik-Software dafür gemacht sind. Selbst mit diesem Wissen ist das Erlernen von Gromacs-Befehlen und Parametrisierung eine weitere steile Lernkurve.

Seit seiner ersten Anwendung in der Biologie im Jahr 19776 hat sich viel weiterentwickelt, was auf eine verstärkte Rechenverarbeitung und eine verbesserte Kodierung zurückzuführen ist. Vor mehr als zwei Jahrzehnten wurde die erste MD-Software für biologische Probleme auf den Markt gebracht, nämlich Gromacs7, AMBER8 und NAMD9.

Seit ihren ersten Versionen sind diese Software immer noch die am häufigsten verwendeten und zitierten. Sie haben jedoch weiterhin die gleichen allgemeinen Implementierungsschwierigkeiten, die Forscher plagen, die keine Computerspezialisten sind5. Einige haben komplexe Installations- und Konfigurationsschritte, die manchmal umfangreiche Kenntnisse über die Hardware erfordern, auf der sie ausgeführt werden, um das Beste daraus zu machen, und eine stark computerzentrierte technische Dokumentation. Eine einfachere Möglichkeit, mit ihnen zu kommunizieren, abgesehen von der Befehlszeile und unendlichen Parametern, ist erforderlich.

Eine Schnittstelle fungiert als Vermittler zwischen dem auszuführenden logischen Prozess und dem Menschen10. Das Paradigma, wie Software ausgeführt wird, hat sich mit der Verbesserung der Rechenressourcen weiterentwickelt. Das erste digitale Paradigma waren die Command Line Interfaces (CLI), gefolgt von der Entwicklung zu den bekannten grafischen Benutzeroberflächen (GUI)11. Dem Evolutionszyklus folgend, wird die vom World Wide Web (oder einfach WEB) erzeugte Schnittstelle als eine Evolution der GUIsbetrachtet 11. Diese drei Paradigmen existieren derzeit je nach Entwickler nebeneinander. CLI-Anwendungen verwenden Textbefehle auf der Betriebssystemkonsole. GUI-Anwendungen, auch grafische Desktops genannt, verwenden eine grafische Oberfläche, die aus Fenstern, Schaltflächen und anderen Komponenten besteht. Es ist spezifisch und für ein Betriebssystem vorprogrammiert. Der Hauptunterschied zur CLI besteht in der Verwendung der Maus als zusätzliches Element in der Mensch-Maschine-Interaktion12. WEB-Anwendungen sind zwar mit einer GUI verwechselt, aber komplexer in der Entwicklung, aber vielseitiger und bei weitem am agilsten im Betrieb. Darüber hinaus hängen sie nur von einer Interpretersoftware ab, die als Browser bezeichnet wird und die es der Client-Anwendung ermöglicht, mit dem Server über ein vom Betriebssystemunabhängiges Netzwerk 13 zu kommunizieren.

Strukturelle Bioinformatik-Software verwendet am häufigsten CLI- und GUI-Paradigmen. Einige Beispiele für klassische Software, die CLI verwendet, sind Modeller14 für die Ähnlichkeitsmodellierung, Autodock15 für molekulares Docking und Gromacs16 für die Molekulardynamik. Beispiele für Software, die den GUI-Typ verwenden, sind unter anderem SwissPDBviewer17, Pymol18, VMD19, UCSF Chimera20, Autodock Tools15, PyRx21, Biovia22, Maestro23 und Moe24.

Mit dem Aufkommen der Hypertext Markup Language Version 5 (HTML5)25, Cascading Style Sheets (CSS)26 und Javascript27 konnten unter anderem viele strukturelle Bioinformatik-Anwendungen ins WEB gebracht und so leichter zugänglich gemacht werden. Beispiele für Similarity Modeling WEB-Server sind MODWEB28, das Modeller14 als Backend verwendet, und Swissmodel29. Beispiele für Webanwendungsserver für molekulares Docking sind Haddock30, Swissdock31, Cluspro32, Dockthor33 und andere.

Während sich die Methoden der Strukturanalyse, Modellierung und des Andockens von CLI-Paradigmen zu GUI und schließlich zu WEB entwickelt haben, wird MD weiterhin hauptsächlich durch Befehlszeilenausführung (CLI-Typ) unterstützt. Es sind einige gute Initiativen entstanden, um dieses Panorama zu verbessern. Beispiele für diese Initiativen sind die Implementierung von Plugins in bestehende Software, wie z. B. das QwikMD-Plugin für VMD34, das GROMACS-Plugin für PyMOL und die Option für die Simulation der Molekulardynamik in UCSF Chimera20, einige neue und einfachere CLI-Anwendungen wie ASGARD35, Gmx_qk36 und CHAPERONg37 sowie eine robuste Webplattform, BioBB-Wfs38. Obwohl die Verwendung dieser Plugins und Anwendungen ein Fortschritt ist, ist ihre Implementierung für die meisten ungelernten Forscher immer noch eine Herausforderung. Häufige Schwierigkeiten sind Probleme bei der Installation und Konfiguration der MD-Software, die oft die vollständige Ausführung der Simulation beeinträchtigen5.

Im Jahr 2022 wurde die Visual Dynamics-Software für webbasierte Computersimulationen vom Laboratório de Bioinformática e Química Medicinal in Fiocruz Rondônia39 zur Verfügung gestellt. Die erste Version wurde in Python und Flask erstellt und ermöglichte Simulationen von Systemen mit freien Proteinen (Apoenzymen) für nur 2 ns. Anschließend wurde es um eine automatisierte Simulationsversion mit Liganden erweitert, die mit PRODRG40 präpariert wurden.

VD wurde entwickelt, um alle Forscher auf dem Gebiet der strukturellen Biophysik, Biotechnologie und verwandten Bereichen zu unterstützen, die nur begrenzte rechnerische Kenntnisse haben. Das Tool ermöglicht es diesen Forschern, ihre Hypothesen mit MD-Simulationen von jedem Betriebssystem aus und ohne Zugang zu einem Hochleistungsrechner (HPC) zu testen. Der Zweck dieser Arbeit besteht darin, die neuen Funktionen von Visual Dynamics Version 3.0 vorzustellen. Darüber hinaus soll ein aktualisiertes Nutzungsprotokoll für das Tool eingeführt und die in Zukunft zu behebenden Einschränkungen sowie Nutzungsstatistiken bis zum heutigen Zeitpunkt hervorgehoben werden (Abbildung 1).

Protocol

1. Zugriff auf die Software und Registrierung neuer Benutzer Besuchen Sie die Webseite von Visual Dynamics (VD). Klicken Sie oben rechts auf das Symbol +Registrieren , um ein Konto zu erstellen. Registrieren Sie sich, um die Software zu nutzen.HINWEIS: Es sind nur institutionelle E-Mail-Adressen zulässig. Der Benutzer erhält eine E-Mail-Benachrichtigung, sobald seine Registrierung genehmigt wurde. Klicken Sie oben rechts auf Anmelden , um auf de…

Representative Results

VD bietet eine vollständig autonome Simulationsausführung, die keine Benutzereingriffe oder vom Benutzer bereitgestellte Rechenressourcen erfordert. Nachdem der Benutzer eine Simulation zur Ausführung eingereicht hat, kann er sie verlassen, seine Computer ausschalten und die Simulation wird weiter ausgeführt. Es ermöglicht Benutzern auch, von jedem Gerät aus auf die Ergebnisse zuzugreifen, sei es ein Laptop oder ein mobiles Gerät. Als Beispiel für die Verwendung von VD im automatisiert…

Discussion

Die Automatisierung von Prozessen ist nicht einfach, aber auch weniger schwierig, als ein System von Grund auf neu zu programmieren. Gromacs ist derzeit die beliebteste molekulare Simulationssoftware und wird ständig aktualisiert. Das Institut für Biophysikalische Chemie der Universität Groningen hat es ursprünglich entwickelt und wird heute vom Labor für Biowissenschaften der Universität Stockholmbetreut 43.

Für jeden neuen Benutzer ist das Erlernen von Simulati…

Offenlegungen

The authors have nothing to disclose.

Acknowledgements

Diese Arbeit wurde unterstützt von der Fundação Oswaldo Cruz (Fiocruz), der Fundação para o Desenvolvimento Científico e Tecnológico em Saúde (Fiotec), dem Instituto Nacional de Ciência e Tecnologia de Epidemiologia da Amazônia Ocidental – INCT-EpiAmO, der Fundação Rondônia de Amparo ao Desenvolvimento das Ações Científicas e Tecnológicas e à Pesquisa do Estado de Rondônia (FAPERO), die Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) und der 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.

Referenzen

  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).

Tags

Play Video

Diesen Artikel zitieren
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