Method Article

Neue Funktionen in Visual Dynamics 3.0

DOI:

10.3791/66964

August 9th, 2024

In This Article

Summary

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

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

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

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

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

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

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

1. Zugriff auf die Software und Registrierung neuer Benutzer

  1. 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.
  2. Klicken Sie oben rechts auf Anmelden , um auf den Anmeldebildschirm des Systems zuzugreifen. Füllen Sie die Felder Benutzername/E-Mail und Passwort aus und klicken Sie auf Anmelden. Nach der Anmeldung hat der Benutzer Zugriff auf seinen Einreichungsbereich für Simulationen. Sie können auch Tutorials und Nutzungsstatistiken für VD anzeigen.

2. Einreichung der Apoenzym-Simulation

  1. Klicken Sie in der linken Seitenleiste auf Neue Simulation . Klicken Sie auf dem angezeigten Bildschirm auf die Schaltfläche APO, die sich auf das Apoenzym bezieht.
  2. Laden Sie die kostenlose Proteindatei 4mvg.pdb hoch. Wählen Sie das AMBER94 Kraftfeld (oder eine andere geeignete Option).
    HINWEIS: 4mvg.pdb kann aus der Ergänzungsdatei 1 bezogen oder direkt von der Proteindatenbank (PDB) heruntergeladen werden.
  3. Wählen Sie das TIP3P-Wassermodell aus. Wählen Sie den kubischen Quader aus. Wählen Sie einen Abstand von 0,5 nm zwischen Protein und Box-Rand.
    HINWEIS: Bei den ausgewählten Optionen handelt es sich nur um Vorschläge. Alle anderen Optionen funktionieren in VD.
  4. Aktivieren Sie die Option Auf unseren Servern ausführen , um die Simulation mit den übermittelten Dateien und Parametern auszuführen. Nach dem Klicken auf Simulation ausführen sieht der Benutzer die Entwicklung der Simulationsschritte auf dem Bildschirm und erhält nach Abschluss eine E-Mail-Benachrichtigung über den Simulationsstatus.

3. Einreichung der Simulation eines Enzyms, das mit einem in ACPYPE hergestellten Liganden komplexiert ist

  1. Öffnen Sie mit UCSF Chimera20 den Protein-Liganden-Komplex 4mgv.pdb41, klicken Sie unter Select auf Residue, und legen Sie den Code auf D5I fest. Klicken Sie dann unter Datei auf PDB speichern, wählen Sie Nur ausgewählte Atome speichern, legen Sie den Dateinamen auf ligand.pdb fest, und klicken Sie auf Speichern.
  2. Übermitteln Sie die im vorherigen Schritt generierte Datei ligand.pdb an den Bio2Byte ACPYPE-Server42, wobei aus den Ausgabedateien ligand_NEW.itp und ligand_NEW.pdb in diesem Experiment verwendet und bereitgestellt werden.
    HINWEIS: Die Dateien ligand_NEW.itp und ligand_NEW.pdb können aus den Zusatzdateien 2 und 3 abgerufen werden.
  3. Klicken Sie in der linken Seitenleiste auf Neue Simulation . Klicken Sie auf die Schaltfläche Protein + Ligand (in ACPYP hergestellt).
  4. Laden Sie die kostenlose Proteindatei 4mvg.pdb hoch. Wählen Sie die in ACPYPE vorbereiteten Ligandendateien aus: ligand_NEW.itp und ligand_NEW.pdb. Wählen Sie das AMBER94 Kraftfeld aus.
  5. Wählen Sie das TIP3P-Wassermodell aus. Wählen Sie den kubischen Quader aus. Wählen Sie einen Abstand von 0,5 nm zwischen Protein und Box-Rand. Bei den ausgewählten Optionen handelt es sich lediglich um Vorschläge. Alle anderen Optionen funktionieren in VD.
  6. Aktivieren Sie die Option Auf unseren Servern ausführen , um die Simulation mit den übermittelten Dateien und Parametern auszuführen. Nach dem Klicken auf Simulation ausführen sieht der Benutzer die Entwicklung der Simulationsschritte auf dem Bildschirm und erhält nach Abschluss eine E-Mail-Benachrichtigung über den Simulationsstatus.

4. Zugriff auf die Simulationsergebnisse

  1. Klicken Sie in der linken Seitenleiste auf Meine Simulationen .
  2. Klicken Sie auf MDP-Dateien herunterladen , um die von der Plattform verwendeten Simulationskonfigurationsdateien auf den Computer des Benutzers herunterzuladen.
  3. Laden Sie Simulationselemente wie unten beschrieben herunter.
    1. Klicken Sie auf Befehle , um die Liste der von der Plattform ausgeführten Befehle auf den Computer des Benutzers herunterzuladen.
    2. Klicken Sie auf GROMACS Log , um die Protokolldatei mit den nacheinander ausgegebenen .gmx-Befehlen auf den Computer des Benutzers herunterzuladen.
    3. Klicken Sie auf Ergebnisse , um Dateien, die mit .gmx-Befehlen generiert wurden, wie z. B. _npt.gro, _pr.edr, _pr.tpr, _pr_PBC.gro und pr_PBC.xtc, auf den Computer des Benutzers herunterzuladen.
    4. Klicken Sie auf Figurengrafiken , um Diagramme zur Analyse der einzelnen Simulationsschritte im Bild- und .xvg-Format auf den Computer des Benutzers herunterzuladen.

Results

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

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 automatisierten Modus über das WEB wurde der Test für einen Protein-Liganden-Komplex durchgeführt, der in ACPYP unter Verwendung der Struktur des FK506-bindenden Proteins aus Plasmodium vivax hergestellt wurde, das mit dem Inhibitor D5 (PDB-ID: 4mgv)41 komplexiert wurde. Die Vorbereitung erfolgte nach dem beschriebenen Protokoll, und die Ergebnisse der Analyse sind in Abbildung 2A-D dargestellt.

Abbildung 2A zeigt die mittlere quadratische Abweichung (RMSD) zwischen der anfänglichen Proteinstruktur und ihrer Simulation über einen Zeitraum von 5 ns (wie im System festgelegt). Das Proteinrückgrat wies während der gesamten Simulation einen RMSD von weniger als 2,5 Å auf. Abbildung 2B zeigt den Gyrationsradius (Rg), der die Kompaktheit des Proteins während der 5-ns-Simulation beschreibt. Dieses Diagramm zeigt Rg in den drei Koordinaten x, y und z sowie den Gesamtwert. Abbildung 2C veranschaulicht die Root Mean Square Fluctuation (RMSF), die den durchschnittlichen Fluktuationsabstand jeder Aminosäure in der Proteinstruktur während der 5-ns-Simulation darstellt. Abbildung 2D zeigt die Energievariation in kJ/mol für das System während des Energieminimierungsprozesses unter Verwendung der Methode des steilsten Abstiegs. Aus dieser Grafik geht hervor, dass sich das System mit einer maximalen Kraft von weniger als 1000 kJ/mol/nm stabilisiert hat.

Ein weiterer Anwendungsfall ist, wenn der Benutzer es auf seinem eigenen Server ausführen möchte, auf dem Gromacs installiert ist. Diese Form der Nutzung wurde in diesem Protokoll nicht behandelt, da sie ein moderates Maß an Kenntnissen über das Linux-Terminal und CLI-Anwendungen erfordert. Diese Methode verwendet VD als Generator für änderbare Gromacs-Befehle, die das gleiche Ergebnis wie VD erzielen, wenn sie wie generiert ausgeführt werden. Sie laden die MDP-Dateien und die Datei mit den generierten Befehlen herunter (Abbildung 3A). Alles ist in VD verfügbar. Wenn sie es im Linux-Terminal ausführen, erhalten sie Ergebnisse wie in Abbildung 3B.

figure-results-1
Abbildung 1: Verwenden von VD zum Generieren von Befehlsskripten. (A) Liste der Befehle, die im VD generiert werden, wenn die Option zum Herunterladen der Befehlsliste ausgewählt wird. (B) Beispiel für die Ausgabe der Ausführung eines Befehls. Kopieren Sie einfach, fügen Sie es ein und führen Sie es in der lokalen Linux-Eingabeaufforderung aus. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

figure-results-2
Abbildung 2: Zugriffs- und Nutzungsstatistiken von Visual Dynamics. Mehr als 4 Tausend Single-IP-Nutzer aus rund 63 Ländern haben bereits auf VD zugegriffen. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

figure-results-3
Abbildung 3: Beispiele für Ausgabediagramme im Zusammenhang mit der automatischen VD-Analyse. (A) Mittlere quadratische Abweichung (RMSD) über einen Verlauf von 5 ns. (B) Gyrationsradius (Rg) über einen Verlauf von 5 ns. (c) Mittlere quadratische Fluktuation (RMSF) über einen Zeitraum von 5 ns. Die x-Achse des RMSF-Diagramms stellt die Anzahl der Aminosäuren im Enzym dar, und die Bildgebungssoftware gibt eine minimale Größe aus. Wenn das Protein größer ist, vergrößert sich die x-Achse. (D) Energievariation während des Energieminimierungsprozesses. Die Energievariation endet, wenn der Prozess einen angemessenen Wert erreicht. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Ergänzende Abbildung 1: Vereinfachtes Flussdiagramm des Prozesses, den Visual Dynamics zum Einrichten einer neuen Simulation verwendet. Bitte klicken Sie hier, um diese Datei herunterzuladen.

Ergänzende Abbildung 2: Vereinfachtes Flussdiagramm des Prozesses, den Visual Dynamics verwendet, um die Ausführung von Simulationen in der Warteschlange zu verwalten. Bitte klicken Sie hier, um diese Datei herunterzuladen.

Ergänzende Datei 1: Protein-pdb-Datei. Bitte klicken Sie hier, um diese Datei herunterzuladen.

Ergänzende Datei 2: Liganden-ITP-Datei. Bitte klicken Sie hier, um diese Datei herunterzuladen.

Ergänzende Datei 3: Liganden-pdb-Datei. Bitte klicken Sie hier, um diese Datei herunterzuladen.

Discussion

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

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 Simulationstechniken eine langwierige Reise. VD bietet sich als Alternative an, um diesen Lernprozess zu erleichtern und Experimente zu validieren, bevor Ressourcen in die Durchführung von Simulationen für längere Zeiträume investiert werden. Der Benutzer kann eine Simulation von 5 ns ausführen und so den Nutzen einer weiteren Erweiterung der Simulation beurteilen.

Die Automatisierung von Verkettungen und Parameterübergaben an Gromacs-Befehle ist nur durch die Art und Weise möglich, wie Linux sein Terminal bereitstellt und verwaltet. Eine der Hauptfunktionen von VD, die eine vollständige Automatisierung ermöglichte, ist in Abbildung 3B zu sehen. Bei einem manuellen Ansatz ist es notwendig, dass der Benutzer zu Beginn der Ausführung interagiert, um die zu analysierenden Gruppen auszuwählen16. In diesem Fall wurde die Gruppe 4 (Backbone) als Standard ausgewählt. Um diesen Schritt zu umgehen und die Benutzerinteraktion zu automatisieren, wird das Befehlsecho "4 4" im Hintergrund an den GMX-Befehl übergeben, zusammen mit einem Verkettungstool, das als Pipe (|) in der Befehlszeile selbst bezeichnet wird. Dies ist nützlich für Programmierer, die ein ähnliches Tool wie VD entwickeln möchten und ermöglicht die Automatisierung und Emulation der Benutzerinteraktion.

Derzeit ist VD auf die Ausführung von Simulationen mit einer Dauer von 5 ns auf dem WEB-Server beschränkt. In der vorherigen Version lag dieser Grenzwert nur bei 2 ns39. Wenn eine längere Simulationszeit benötigt wird, wird empfohlen, das Skript und die MDP-Dateien herunterzuladen, um sie auf dem eigenen Computer auszuführen, auf dem Gromacs installiert ist. Die Optionen System neutralisieren und Wasserstoffatome ignorieren sind immer aktiviert und können nicht deaktiviert werden. Ebenso ist die Option Doppelte Genauigkeit verwenden immer deaktiviert. Diese Optionen werden angezeigt, um den Benutzer zu informieren. Aus Performance-Gründen sind sie fixiert. Die aktuelle Version führt keine Simulationen mit Liganden durch, die in OPLS- oder CHARMM-Kraftfeldern präpariert wurden. Diese Einschränkung ist rein technischer und vorübergehender Natur.

Strukturell verfügt VD über zwei Anwendungen: ein Web-Frontend, d. h. die Anwendung, die der Benutzer sieht, und ein Server-Backend, d. h. die Anwendung, die die harte Arbeit beim Einrichten und Ausführen des MD übernimmt. Das Web-Frontend ist eine einfache Schnittstelle mit Formularen, um zu erfassen, was die Benutzer über Auswahlen, Umschalter und Dateiauswahlen tun möchten. Das Flussdiagramm des Web-Frontends entfällt vorerst, da es Schritt für Schritt im Protokoll beschrieben wird.

Das Server-Backend besteht aus zwei Diensten: einer Anwendung zum Verwalten und Einrichten von MDs und einer weiteren zum gleichzeitigen Verwalten und Ausführen von MDs. Die ergänzende Abbildung 1 zeigt den Ablauf zum Einrichten eines MD, und die ergänzende Abbildung 2 zeigt den Ablauf der Anwendung in der Ausführungswarteschlange.

Neben Visual Dynamics gibt es auch andere automatisierte/Web-Tools zum Ausführen von MD, jedes mit seinen eigenen Vor- und Nachteilen. ASGARD35, Gmx_qk36 und CHAPERONg37 sind hervorragende Werkzeuge, die die Verwendung von GROMACS14 erleichtern, aber immer noch CLI-Anwendungen sind, die ein Mindestmaß an Kenntnissen über die Shell erfordern. BioBB-Wfs38 ist die ähnlichste Anwendung zu VD39, sie bietet auch viele Dinge, die VD nicht bietet, da sie BioBB44 verwendet, eine Bibliothek, die aus Wrappern um beliebte MD-Tools besteht. Die Bildschirme können für jemanden, der gerade erst anfängt, etwas zu komplex sein. Der Prozess ist im Vergleich zu VD recht komplex. Nach vielen Schritten wird eine Art universelle Komponente namens Common Workflow Language (CWL) für die Ausführung auf HPCs generiert. Der Benutzer muss jedoch Zugriff auf diese Server haben, was natürlich nicht kostenlos ist.

In der aktuellen öffentlichen Instanz von VD ist Gromacs installiert, um den Prozessor zu verwenden, aber es kann auch mit Anweisungen zur Verwendung der Grafikprozessor (GPU) des Computers installiert werden, wodurch die Simulationsleistung erheblich verbessertwird 45. Darüber hinaus ist das in VD generierte Befehlsskript GPU-transparent, was bedeutet, dass es auf einem Computer mit installierter GPU ausgeführt werden kann, ohne dass etwas an den Befehlen geändert werden muss.

Zu diesem Zeitpunkt ist VD ein großartiges Werkzeug, um MD jemandem vorzustellen, der neu darin ist, aber es sollte am besten von denjenigen verwendet werden, die mit MD und Gromacs vertraut sind. Es kann in Simulationsklassen verwendet werden, in denen der Professor eine schnelle Herangehensweise ohne Installations- oder Parameterkonfigurationsprobleme benötigt, um eine Übung auszuführen. Darüber hinaus kann VD verwendet werden, um eine Simulation zu validieren. Wenn es 5 ns läuft, ist es vertretbar, in die Erweiterung der Simulation mit den in VD.In generierten Dateien zu investieren. Im Allgemeinen arbeitet die Gruppe daran, stärker automatisierte Analysen wie h-Bindungen und die Bindungsenergie der molekularen Mechanik Poisson-Boltzmann-Oberfläche (MMPBSA) einzubeziehen und gleichzeitig kontinuierlich die Benutzerfreundlichkeit zu verbessern.

Disclosures

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

Die Autoren haben nichts offenzulegen.

Acknowledgements

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

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

List of materials used in this article
NameCompanyCatalog NumberComments
ACPYPE ServerBio2ByteErhältlich bei https://www.bio2byte.be/acpype/
Erhältlich bei https://plasma-gate.weizmann.ac.il/Grace/
GROMACS SoftwareGROMACS TeamInstallationsanleitung bei https://manual.gromacs.org/current/install-guide/index.html
Die Struktur des FK506-bindenden Proteins
Aus Plasmodium vivax, das mit dem
-Inhibitor D5
komplexiert ist, RCSB-ProteindatenbankVerfügbar unter https://www.rcsb.org/structure/4mgv
Enthält bereits den Liganden, der mit dem Makromolekül komplexiert ist.
GRACE Software Plasma Laboratory am Weizmann Institute of Science

References

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,
  1. The IUPAC Compendium of Chemical Terminology: The Gold Book. , International Union of Pure and Applied Chemistry (IUPAC). Research Triangle Park, NC. (2019).">International Union of Pure and Applied Chemistry (IUPAC). The IUPAC Compendium of Chemical Terminology: The Gold Book. , International Union of Pure and Applied Chemistry (IUPAC). Research Triangle Park, NC. (2019).
  2. Molecular dynamics simulations of biomolecules. Nat Str Biol. 9 (9), 646-652 (2002).">Karplus, M., McCammon, J. A. Molecular dynamics simulations of biomolecules. Nat Str Biol. 9 (9), 646-652 (2002).
  3. Molecular dynamics: Survey of methods for simulating the activity of proteins. Chem Rev. 106 (5), 1589-1615 (2006).">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. Improved simulation of liquid water by molecular dynamics. J Chem Phys. 60 (4), 1545-1557 (1974).">Stillinger, F. H., Rahman, A. Improved simulation of liquid water by molecular dynamics. J Chem Phys. 60 (4), 1545-1557 (1974).
  5. Enlighten2: molecular dynamics simulations of protein-ligand systems made accessible. Bioinformatics. 36 (20), 5104-5106 (2020).">Zinovjev, K., Van Der Kamp, M. W. Enlighten2: molecular dynamics simulations of protein-ligand systems made accessible. Bioinformatics. 36 (20), 5104-5106 (2020).
  6. Dynamics of folded proteins. Nature. 267 (5612), 585-590 (1977).">McCammon, J. A., Gelin, B. R., Karplus, M. Dynamics of folded proteins. Nature. 267 (5612), 585-590 (1977).
  7. GROMACS: A message-passing parallel molecular dynamics implementation. Comp Phys Comm. 91 (1-3), 43-56 (1995).">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. 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).">Pearlman, D. A., et al. 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. NAMD: A parallel, object-oriented molecular dynamics program. Int J Supercomp Appl High Performance Comp. 10 (4), 251-268 (1996).">Nelson, M. T., et al. 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. , Addison-Wesley Pub. Co. Reading, Mass. (1990).">Addison-Wesley Pub. Co. The Art of human-computer interface design. , Addison-Wesley Pub. Co. Reading, Mass. (1990).
  11. Brave NUI World: Designing Natural User Interfaces for Touch and Gesture. , Elsevier Science. Saint Louis. (2014).">Wigdor, D., Wixon, D. Brave NUI World: Designing Natural User Interfaces for Touch and Gesture. , Elsevier Science. Saint Louis. (2014).
  12. GUI and Command-line - Conflict or Synergy. , (2000).">Unwin, A., Heike, H. GUI and Command-line - Conflict or Synergy. , (2000).
  13. Beginning C# Object-Oriented Programming. , 243-263 (2011).">Clark, D. Developing Web Applications. Beginning C# Object-Oriented Programming. , 243-263 (2011).
  14. Comparative protein structure modeling using Modeller. Curr Protoc Bioinformatics. 5, 10(2006).">Eswar, N., et al. Comparative protein structure modeling using Modeller. Curr Protoc Bioinformatics. 5, 10(2006).
  15. AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J Comp Chem. 30 (16), 2785-2791 (2009).">Morris, G. M., et al. AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J Comp Chem. 30 (16), 2785-2791 (2009).
  16. GROMACS: Fast, flexible, and free. J Comp Chem. 26 (16), 1701-1718 (2005).">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. SWISS-MODEL and the Swiss-Pdb Viewer: An environment for comparative protein modeling. Electrophoresis. 18 (15), 2714-2723 (1997).">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. https://pymol.org/ (2023).">The PyMOL molecular graphics system. , Available from: https://pymol.org/ (2023).
  19. VMD: Visual molecular dynamics. J Mol Graph. 14 (1), 33-38 (1996).">Humphrey, W., Dalke, A., Schulten, K. VMD: Visual molecular dynamics. J Mol Graph. 14 (1), 33-38 (1996).
  20. UCSF Chimera-A visualization system for exploratory research and analysis. J Comp Chem. 25 (13), 1605-1612 (2004).">Pettersen, E. F., et al. UCSF Chimera-A visualization system for exploratory research and analysis. J Comp Chem. 25 (13), 1605-1612 (2004).
  21. https://pyrx.sourceforge.io/ (2023).">PyRx-python prescription v. 0.8. , Available from: https://pyrx.sourceforge.io/ (2023).
  22. https://www.3ds.com/products/biovia/discovery-studio (2023).">BIOVIA Discovery Studio 2017 R2: A comprehensive predictive science application for the Life Sciences. , Available from: https://www.3ds.com/products/biovia/discovery-studio (2023).
  23. https://www.schrodinger.com/platform/products/maestro/ (2024).">Maestro. , Available from: https://www.schrodinger.com/platform/products/maestro/ (2024).
  24. https://www.chemcomp.com/en/Products.htm (2023).">Molecular operating environment (MOE). , Available from: https://www.chemcomp.com/en/Products.htm (2023).
  25. https://dev.w3.org/html5/html-author/ (2023).">The Syntax, Vocabulary and APIs of HTML5. , Available from: https://dev.w3.org/html5/html-author/ (2023).
  26. https://www.w3.org/Style/CSS/ (2023).">Cascading Style Sheets. , Available from: https://www.w3.org/Style/CSS/ (2023).
  27. https://www.ecma-international.org/publications-and-standards/standards/ecma-262/ (2023).">ECMAScript 2020 language specification, 11th edition. , Available from: https://www.ecma-international.org/publications-and-standards/standards/ecma-262/ (2023).
  28. ModBase, a database of annotated comparative protein structure models, and associated resources. Nucleic Acids Res. 39, Database D465-D474 (2011).">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. The SWISS-MODEL Repository and associated resources. Nucleic Acids Res. 37, Database D387-D392 (2009).">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. The HADDOCK web server for data-driven biomolecular docking. Nat Prot. 5 (5), 883-897 (2010).">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. SwissDock, a protein-small molecule docking web service based on EADock DSS. Nucleic Acids Res. 39, suppl W270-W277 (2011).">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. The ClusPro web server for protein-protein docking. Nat Prot. 12 (2), 255-278 (2017).">Kozakov, D., et al. The ClusPro web server for protein-protein docking. Nat Prot. 12 (2), 255-278 (2017).
  33. 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).">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. QwikMD - Integrative molecular dynamics toolkit for novices and experts. Sci Rep. 6 (1), 26536(2016).">Ribeiro, J. V., et al. QwikMD - Integrative molecular dynamics toolkit for novices and experts. Sci Rep. 6 (1), 26536(2016).
  35. ASGARD. A simple and automatic GROMACS tool to analyze Molecular Dynamic simulations. , (2023).">Rodríguez Martínez, A., et al. ASGARD. A simple and automatic GROMACS tool to analyze Molecular Dynamic simulations. , (2023).
  36. 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).">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. CHAPERONg: A tool for automated GROMACS-based molecular dynamics simulations and trajectory analyses. Comp Str Biotechnol J. 21, 4849-4858 (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).
  38. BioExcel building blocks workflows (BioBB-Wfs), an integrated web-based platform for biomolecular simulations. Nucleic Acids Res. 50 (W1), W99-W107 (2022).">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), W99-W107 (2022).
  39. Visual dynamics: a WEB application for molecular dynamics simulation using GROMACS. BMC Bioinfo. 24 (1), 107(2023).">Vieira, I. H. P., et al. Visual dynamics: a WEB application for molecular dynamics simulation using GROMACS. BMC Bioinfo. 24 (1), 107(2023).
  40. PRODRG a tool for high-throughput crystallography of protein-ligand complexes. Acta Crystallographica Sect D Biol Crystallography. 60 (8), 1355-1363 (2004).">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. Adamantyl derivative as a potent inhibitor of Plasmodium FK506 binding protein 35. ACS Med Chem Lett. 4 (11), 1097-1101 (2013).">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. ACPYPE - AnteChamber PYthon Parser interfacE. BMC Res Notes. 5 (1), 367(2012).">Sousa Da Silva, A. W., Vranken, W. F. ACPYPE - AnteChamber PYthon Parser interfacE. BMC Res Notes. 5 (1), 367(2012).
  43. https://www.gromacs.org/development.html (2023).">Gromacs developmet and developers. , Available from: https://www.gromacs.org/development.html (2023).
  44. BioExcel building blocks, a software library for interoperable biomolecular simulation workflows. Sci Data. 6 (1), 169(2019).">Andrio, P., et al. BioExcel building blocks, a software library for interoperable biomolecular simulation workflows. Sci Data. 6 (1), 169(2019).
  45. Heterogeneous parallelization and acceleration of molecular dynamics simulations in GROMACS. J Chem Phys. 153 (13), 134110(2020).">Páll, S., et al. Heterogeneous parallelization and acceleration of molecular dynamics simulations in GROMACS. J Chem Phys. 153 (13), 134110(2020).

Reprints and Permissions

Request permission to reuse the text or figures of this JoVE article

Request Permission

Tags

Molecular DynamicsVisual DynamicsGromacs SimulationProtein Ligand ComplexSimulation PipelineForce Field SelectionACPYPE PreparationTIP3P Water ModelRoot Mean SquareEnergy Minimization

Related Articles