Novos recursos no Visual Dynamics 3.0

Resumo

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.

Resumo

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.

Introdução

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 molecular¹. Pode ser aplicado a qualquer fenômeno que busque extrair informações em nível molecular e muitas vezes atômico².

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ção³. 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)³. Os resultados obtidos ao final de um MD são os mais ricos e completos em termos de simulação não quântica⁴. 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 ignorados³. No entanto, implementar um experimento de MD não é trivial⁵. 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 1977⁶, 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, Gromacs⁷, AMBER⁸ e NAMD⁹.

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ção⁵. 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 humano¹⁰. 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 GUIs¹¹. 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áquina¹². 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 operacional¹³.

O software de bioinformática estrutural geralmente usa paradigmas CLI e GUI. Alguns exemplos de software clássico que usam CLI são o Modeller¹⁴ para modelagem de similaridade, o Autodock¹⁵ para acoplamento molecular e o Gromacs¹⁶ para dinâmica molecular. Exemplos de softwares que adotam o tipo GUI são SwissPDBviewer¹⁷, Pymol¹⁸, VMD¹⁹, UCSF Chimera²⁰, Autodock tools¹⁵, PyRx²¹, Biovia²², Maestro²³ e Moe²⁴, entre outros.

Com o surgimento das tecnologias Hypertext Markup Language versão 5 (HTML5)²⁵, Cascading Style Sheets (CSS)²⁶ e Javascript²⁷ , 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 MODWEB²⁸, que usa o Modeller¹⁴ como back-end e o Swissmodel²⁹. Exemplos de servidores de aplicativos da Web para docking molecular são Haddock³⁰, Swissdock³¹, Cluspro³², Dockthor³³ 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 VMD³⁴, o plug-in GROMACS para PyMOL e a opção Molecular Dynamics Simulation no UCSF Chimera²⁰, alguns aplicativos CLI novos e mais fáceis, como ASGARD³⁵, Gmx_qk³⁶ e CHAPERONg³⁷, e uma plataforma web robusta, BioBB-Wfs³⁸. 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ção⁵.

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ônia³⁹. 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 PRODRG⁴⁰.

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

Protocolo

1. Acesso ao software e registro de novo usuário

1. Visite a página da Web do Visual Dynamics (VD). Clique no ícone +Registrar no canto superior direito para criar uma conta. Registre-se para usar o software.
   NOTA: Somente endereços de e-mail institucionais são permitidos. O usuário receberá uma notificação por e-mail assim que seu registro for aprovado.
2. Clique em Login no canto superior direito para acessar a tela de login do sistema. Preencha os campos de nome de usuário/e-mail e senha e clique em Entrar. Após o login, o usuário terá acesso à sua área de envio de simulação. Eles também podem visualizar tutoriais e estatísticas de uso do VD.

2. Submissão de simulação de apoenzima

1. Clique em Nova Simulação na barra lateral esquerda. Na tela que aparece, clique no botão APO, que se refere à apoenzima.
2. Carregue o arquivo 4mvg.pdb de proteína grátis. Selecione o campo de força AMBER94 (ou qualquer outra opção apropriada).
   NOTA: 4mvg.pdb pode ser obtido no Arquivo Suplementar 1 ou baixado diretamente do Protein Data Bank (PDB).
3. Selecione o modelo de água TIP3P. Selecione a caixa cúbica. Selecione a distância de 0,5 nm entre a proteína e a borda da caixa.
   NOTA: As opções selecionadas são apenas sugestões. Todas as outras opções funcionam no VD.
4. Marque a opção Executar em Nossos Servidores para executar a simulação com os arquivos e parâmetros enviados. Após clicar em Executar simulação, o usuário verá a evolução das etapas da simulação na tela e receberá uma notificação por e-mail sobre o status da simulação após a conclusão.

3. Submissão de Simulação de Enzima Complexada com Ligante Preparado em ACPYPE

1. Usando o UCSF Chimera²⁰, abra o complexo proteína-ligante 4mgv.pdb⁴¹, em Selecionar, clique em Resíduo e defina o código como D5I. Em seguida, em Arquivo, clique em Salvar PDB, selecione Salvar Somente Átomos Selecionados, defina o nome do arquivo como ligand.pdb e clique em Salvar.
2. Envie o arquivo ligand.pdb gerado na etapa anterior para o servidor Bio2Byte ACPYPE⁴², a partir dos arquivos de saída, ligand_NEW.itp e ligand_NEW.pdb serão os utilizados e fornecidos neste experimento.
   NOTA: Os arquivos ligand_NEW.itp e ligand_NEW.pdb podem ser obtidos no Arquivo Suplementar 2 e no Arquivo Suplementar 3.
3. Clique em Nova Simulação na barra lateral esquerda. Clique no botão Proteína + Ligante (preparado em ACPYPE).
4. Carregue o arquivo 4mvg.pdb de proteína grátis. Selecione os Arquivos Ligantes preparados no ACPYPE: ligand_NEW.itp e ligand_NEW.pdb. Selecione o campo de força AMBER94.
5. Selecione o modelo de água TIP3P. Selecione a caixa cúbica. Selecione a distância de 0,5 nm entre a proteína e a borda da caixa. As opções selecionadas são apenas sugestões. Todas as outras opções funcionam no VD.
6. Marque a opção Executar em Nossos Servidores para executar a simulação com os arquivos e parâmetros enviados. Após clicar em Executar simulação, o usuário verá a evolução das etapas da simulação na tela e receberá uma notificação por e-mail sobre o status da simulação após a conclusão.

4. Acessando os resultados da simulação

1. Clique em Minhas simulações na barra lateral esquerda.
2. Clique em Baixar arquivos MDP para baixar os arquivos de configuração de simulação usados pela plataforma para o computador do usuário.
3. Faça o download dos elementos de simulação conforme descrito abaixo.
   1. Clique em Comandos para baixar a lista de comandos executados pela plataforma para o computador do usuário.
   2. Clique em GROMACS Log para baixar o arquivo de log que contém as saídas do comando .gmx em sequência para o computador do usuário.
   3. Clique em Resultados para baixar arquivos gerados por comandos .gmx, como _npt.gro, _pr.edr, _pr.tpr, _pr_PBC.gro e pr_PBC.xtc para o computador do usuário.
   4. Clique em Figura Gráfica para baixar gráficos para analisar cada etapa da simulação em formato de imagem e .xvg para o computador do usuário.

Resultados

O VD fornece uma execução de simulação totalmente autônoma que não requer intervenção do usuário ou recursos computacionais fornecidos pelo usuário. Depois de enviar uma simulação para execução, o usuário pode deixá-la, desligar suas máquinas e a simulação continuará em execução. Ele também permite que os usuários acessem os resultados de qualquer dispositivo, seja um laptop ou dispositivo móvel.

Como exemplo de uso de VD em modo automatizado através da WEB, o teste f...

Discussão

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 Estocolmo⁴³.

Para qualquer novo usuário, aprender técnicas de simulação é uma longa jor...

Materiais

Name	Company	Catalog Number	Comments
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.