Summary

Nuove funzionalità di Visual Dynamics 3.0

Published: August 09, 2024
doi:

Summary

Visual Dynamics è uno strumento open source che accelera le implementazioni e l’apprendimento nella simulazione della dinamica molecolare utilizzando Gromacs. Il protocollo presentato vi guiderà attraverso i passaggi per eseguire una simulazione proteina-ligando preparata in ACPYPE con facilità e i passaggi generali verso altri modelli di simulazione.

Abstract

Visual Dynamics (VD) è uno strumento web che mira a facilitare l’uso e l’applicazione della Molecular Dynamics (MD) eseguita in Gromacs, consentendo agli utenti senza familiarità computazionale di eseguire simulazioni di breve durata per scopi di convalida, dimostrazione e insegnamento. È vero che i metodi quantistici sono i più accurati. Tuttavia, attualmente non esiste una fattibilità computazionale per eseguire gli esperimenti eseguiti da MD. Lo strumento qui descritto ha ricevuto continui miglioramenti nel corso degli ultimi due anni. Questo protocollo descriverà ciò che è necessario per eseguire una simulazione in VD con un complesso proteina-ligando precedentemente preparato in ACPYPE e alcune indicazioni generali sugli altri modelli di simulazione disponibili. Per la simulazione dettagliata, verrà utilizzata la proteina legante FK506 di Plasmodium vivax complessata con l’inibitore D5 (PDB ID: 4mgv) e verranno forniti tutti i file utilizzati. Si noti che questo protocollo indicherà ogni opzione da utilizzare per ottenere gli stessi risultati presentati, ma queste opzioni non sono necessariamente le uniche disponibili.

Introduction

Secondo la definizione IUPAC, la MD è la procedura di simulazione che consiste nel calcolare il moto degli atomi in una molecola o di singoli atomi o molecole in solidi, liquidi e gas, secondo le leggi del moto di Newton. Le forze che agiscono sugli atomi, necessarie per simulare il loro movimento, sono comunemente calcolate utilizzando i campi di forza della meccanica molecolare1. Può essere applicato a qualsiasi fenomeno che cerchi di estrarre informazioni a livello molecolare e spesso atomico2.

La MD è una delle tecniche incorporate nella bioinformatica, in particolare nella bioinformatica strutturale. Con esso è possibile ottenere caratteristiche cinetiche e termodinamiche di strutture biomolecolari. Ad esempio, stabilità macromolecolare, identificazione di siti allosterici, chiarimento dei meccanismi di attività enzimatica, riconoscimento molecolare e proprietà di complessi con piccole molecole, associazione tra proteine, ripiegamento delle proteine e sua idratazione3. Inoltre, la MD consente un’ampia gamma di studi, tra cui la progettazione molecolare (ampiamente utilizzata nella progettazione di farmaci), per determinare la struttura e il suo raffinamento (raggi X, NMR e modellazione proteica)3. I risultati ottenuti al termine di un MD sono i più ricchi e completi in termini di simulazione non quantistica4. La MD classica è molto più efficiente di quanto ci si potrebbe aspettare da una considerazione completa della fisica dei sistemi biomolecolari a causa del numero di approssimazioni sostanziali. In particolare, gli effetti della dinamica quantistica vengono solitamente ignorati3. Tuttavia, l’implementazione di un esperimento MD non è banale5. Richiede la conoscenza dell’informatica, in particolare del terminale Linux, poiché la maggior parte dei software di bioinformatica strutturale è fatta per questo. Anche con questa conoscenza, l’apprendimento dei comandi e della parametrizzazione di Gromacs è un’altra curva di apprendimento ripida.

Dalla sua prima applicazione alla biologia nel 19776, molto si è evoluto grazie all’aumento dell’elaborazione computazionale e al miglioramento della codifica. Più di due decenni fa, è stato lanciato il primo software MD destinato ai problemi biologici, ovvero Gromacs7, AMBER8 e NAMD9.

Sin dalle loro prime versioni, questi software rimangono ancora oggi i più utilizzati e citati. Tuttavia, continuano con le stesse comuni difficoltà di implementazione che affliggono i ricercatori che non sono specialisti di computer5. Alcuni hanno passaggi di installazione e configurazione complessi, che a volte richiedono una conoscenza approfondita dell’hardware su cui verranno eseguiti per ottenere il massimo da esso e una documentazione tecnica altamente incentrata sul computer. È necessario un modo più semplice per interfacciarsi con loro, oltre alla riga di comando e ai parametri infiniti.

Un’interfaccia funge da intermediario tra il processo logico da eseguire e l’essere umano10. Il paradigma di esecuzione del software si è evoluto con il miglioramento delle risorse informatiche. Il primo paradigma digitale sono state le interfacce a riga di comando (CLI), seguite dall’evoluzione verso le note interfacce grafiche utente (GUI)11. Seguendo il ciclo evolutivo, l’interfaccia prodotta dal World Wide Web (o semplicemente WEB) è considerata un’evoluzione delle GUI11. Questi tre paradigmi attualmente coesistono a seconda degli sviluppatori. Le applicazioni CLI utilizzano comandi testuali sulla console del sistema operativo. Le applicazioni GUI, chiamate anche desktop grafici, utilizzano un’interfaccia grafica composta da finestre, pulsanti e altri componenti. È specifico e pre-programmato per un sistema operativo. La principale differenza rispetto alla CLI è l’uso del mouse come elemento aggiuntivo nell’interazione uomo-macchina12. Le applicazioni WEB, nonostante siano confuse con un’interfaccia grafica, sono più complesse da sviluppare ma sono più versatili e di gran lunga le più agili nel funzionamento. Inoltre, essi dipendono solo da un software di interpretazione, chiamato browser, che consente all’applicazione client di comunicare con il server attraverso una rete indipendente dal sistema operativo13.

I software di bioinformatica strutturale utilizzano più comunemente i paradigmi CLI e GUI. Alcuni esempi di software classici che utilizzano CLI sono Modeller14 per la modellazione di similarità, Autodock15 per il docking molecolare e Gromacs16 per la dinamica molecolare. Esempi di software che adottano il tipo GUI sono SwissPDBviewer17, Pymol18, VMD19, UCSF Chimera20, Autodock tools15, PyRx21, Biovia22, Maestro23 e Moe24, tra gli altri.

Con l’emergere delle tecnologie Hypertext Markup Language versione 5 (HTML5)25, Cascading Style Sheets (CSS)26 e Javascript27 , tra le altre, molte applicazioni di bioinformatica strutturale potrebbero essere portate sul WEB, diventando così più accessibili. Esempi di similarity modeling WEB server sono MODWEB28, che utilizza Modeller14 come back-end e Swissmodel29. Esempi di server di applicazioni web per il docking molecolare sono Haddock30, Swissdock31, Cluspro32, Dockthor33 e altri.

Mentre l’analisi strutturale, la modellazione e le metodologie di docking si sono evolute dai paradigmi CLI alla GUI e infine al WEB, MD continua ad essere per lo più supportato dall’esecuzione della riga di comando (tipo CLI). Sono emerse alcune buone iniziative per migliorare questo panorama. Esempi di queste iniziative sono l’implementazione di plug-in in software esistenti, come il plug-in QwikMD per VMD34, il plug-in GROMACS per PyMOL e l’opzione di simulazione dinamica molecolare in UCSF Chimera20, alcune nuove e più semplici applicazioni CLI, come ASGARD35, Gmx_qk36 e CHAPERONg37, e una robusta piattaforma web, BioBB-Wfs38. Sebbene l’uso di questi plugin e applicazioni sia un progresso, la loro implementazione è ancora una sfida per la maggior parte dei ricercatori non qualificati. Le difficoltà più comuni includono problemi di installazione e configurazione del software MD, che spesso compromettono la piena esecuzione della simulazione5.

Nel 2022, il software Visual Dynamics per la simulazione computazionale basata sul web è stato reso disponibile dal Laboratório de Bioinformática e Química Medicinal presso Fiocruz Rondônia39. La sua versione iniziale è stata costruita in Python e Flask, consentendo simulazioni di sistemi con proteine libere (apoenzimi) per soli 2 ns. Successivamente, è stato migliorato per includere una versione di simulazione automatizzata con ligandi preparati utilizzando PRODRG40.

VD è stato creato per assistere tutti i ricercatori nel campo della biofisica strutturale, della biotecnologia e delle aree correlate che hanno limitazioni nella conoscenza computazionale; lo strumento consente a questi ricercatori di testare le loro ipotesi che coinvolgono simulazioni MD da qualsiasi sistema operativo e senza l’accesso a un computer ad alte prestazioni (HPC). Lo scopo di questo lavoro è quello di presentare le nuove funzionalità di Visual Dynamics versione 3.0. Inoltre, mira a introdurre un protocollo di utilizzo aggiornato per lo strumento e ad evidenziare le limitazioni da affrontare in futuro, insieme alle statistiche di utilizzo fino al momento attuale (Figura 1).

Protocol

1. Accesso al software e registrazione di un nuovo utente Visita la pagina Web di Visual Dynamics (VD). Clicca sull’icona +Registrati in alto a destra per creare un account. Registrati per utilizzare il software.NOTA: Sono consentiti solo indirizzi e-mail istituzionali. L’utente riceverà una notifica via e-mail una volta approvata la registrazione. Clicca su Accedi in alto a destra per accedere alla schermata di accesso al sistema. Compila i camp…

Representative Results

VD fornisce un’esecuzione di simulazione completamente autonoma che non richiede l’intervento dell’utente o risorse computazionali fornite dall’utente. Dopo aver inviato una simulazione all’esecuzione, l’utente può uscirne, spegnere i computer e la simulazione continuerà a essere eseguita. Consente inoltre agli utenti di accedere ai risultati da qualsiasi dispositivo, sia esso un laptop o un dispositivo mobile. Come esempio di utilizzo della VD in modalità automatizzata attraverso il WEB, i…

Discussion

Automatizzare i processi non è facile, ma è anche meno difficile che riprogrammare un sistema da zero. Gromacs è attualmente il software di simulazione molecolare più popolare ed è costantemente aggiornato. Il Dipartimento di Chimica Biofisica dell’Università di Groningen lo ha inizialmente sviluppato, ed è ora gestito dal Laboratorio di Scienze della Vita dell’Università di Stoccolma43.

Per ogni nuovo utente, l’apprendimento delle tecniche di simulazione è un …

Disclosures

The authors have nothing to disclose.

Acknowledgements

Questo lavoro è stato sostenuto dalla Fundação Oswaldo Cruz (Fiocruz), dalla Fundação para o Desenvolvimento Científico e Tecnológico em Saúde (Fiotec), dall’Instituto Nacional de Ciência e Tecnologia de Epidemiologia da Amazônia Ocidental – INCT-EpiAmO, dalla Fundação Rondônia de Amparo ao Desenvolvimento das Ações Científicas e Tecnológicas e à Pesquisa do Estado de Rondônia (FAPERO), il Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) e il 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).

Play Video

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