26 de setembro de 2025
Este artigo apresenta um protocolo para a construção de um modelo de farmacoforo de consenso integrando características moleculares de vários ligantes. Este método é aplicável aos esforços de descoberta de medicamentos direcionados a qualquer alvo biológico com conformações conhecidas de ligação ao ligante, permitindo a identificação dos principais recursos de interação para triagem virtual e design racional de medicamentos.
Nossa pesquisa explora a modelagem de farmacoforos com grandes bibliotecas de ligantes. Neste estudo, nos concentramos no SARS-CoV-2 Mpro para identificar padrões estruturais que orientam a descoberta de medicamentos. Avanços recentes na modelagem de farmacoforos usam conjuntos de dados maiores e mais diversificados, o que melhora a precisão na identificação dos principais alimentadores moleculares para a descoberta de medicamentos.
Simulações dinâmicas moleculares e aprendizado de máquina são usados atualmente para avançar a modelagem de farmacoforos e melhorar os resultados da descoberta de medicamentos. Um dos principais desafios é gerenciar conjuntos de dados de ligantes grandes e diversificados, melhorando a precisão da extração futura e reduzindo os custos computacionais. Mostramos que bibliotecas de ligantes caras melhoram a precisão futura do farmacoforo, permitindo previsões mais confiáveis e melhor identificação de potenciais candidatos a medicamentos.
Para começar, inicie o software PyMOL. Alinhe todos os complexos de ligantes de proteínas no software. Para cada complexo alinhado, extraia o conformador de ligantes e salve-o como um arquivo separado no formato SDF.
Carregue cada arquivo de ligante individualmente para farmá-lo, usando a opção Carregar recursos. Após o carregamento, use a opção Salvar sessão para baixar o arquivo farmacofore no formato JSON. Armazene todos os arquivos JSON do farmacophore baixados juntos em uma única pasta para uso posterior no ambiente do Google Colab.
Inicie o Google Colab em um navegador da Web e crie um novo bloco de anotações. Instale o Conda e o PyMOL e, em seguida, execute a célula clicando no ícone de reprodução ou pressionando shift e Enter para verificar a execução bem-sucedida. Uma barra verde aparecerá acima da célula se concluída corretamente.
Agora, instale o pacote CANFAR Python e importe os módulos necessários. Execute a célula clicando no ícone de reprodução para confirmar a instalação bem-sucedida. Uma mensagem de confirmação aparecerá após a execução bem-sucedida.
Em seguida, crie uma nova pasta para armazenar arquivos JSON do farmacoforo. Abra a pasta, clique com o botão direito dentro dela e selecione carregar para carregar os arquivos JSON. Para analisar e consolidar os recursos do farmacoforo, extraia os recursos de todos os arquivos carregados.
Clique no ícone de reprodução para executar a célula e gerar um quadro de dados consolidado contendo os recursos de cada ligante individual. Exiba todos os descritores de farmacoforos e, em seguida, execute a célula para exibir padrões de agrupamento. Agora, salve o modelo de farmacoforo no formato PyMOL.
Execute a célula para gerar um arquivo pse e, em seguida, salve o modelo como um arquivo JSON executando a célula de código correspondente. Gere saídas e dendogramas agrupados de recursos usando o código fornecido. Clique em reproduzir e exporte os dados do farmacoforo de consenso para um arquivo CSV executando o comando de exportação.
Use o modelo de farmacoforo de consenso gerado para triagem virtual. 100 proteases principais ou complexos Mpro cocristalizados com diferentes inibidores não covalentes foram alinhados para visualizar a diversidade de ligantes dentro da bolsa de ligação. Cada ligante foi extraído individualmente dos complexos alinhados para gerar arquivos separados para triagem virtual.
Os ligantes extraídos foram carregados no servidor pharmit e usados para gerar arquivos JSON estruturados para modelagem de farmacoforos. Os arquivos JSON produziram 1.450 características de farmacóforo, agrupadas em 110 clusters, incluindo 23 aromáticos, 30 aceitadores de ligações de hidrogênio, 16 doadores de ligações de hidrogênio, 36 clusters hidrofóbicos e cinco de ânions. Para construir o modelo de consenso, apenas os clusters mais densamente povoados foram retidos, enquanto os clusters de ânions foram excluídos, devido à insuficiência de membros.
O modelo farmacoforal de consenso final consistiu em 11 características, compreendendo três aromáticas, para aceitadores de ligações de hidrogênio, dois doadores de ligações de hidrogênio e duas características hidrofóbicas. Após a remoção da característica aromática menos representativa, aro um, dois compostos candidatos foram identificados. A estrutura química bidimensional do composto 101 267 741 foi visualizada para confirmar sua compatibilidade com o modelo farmacofórico.
O composto 101 267 741 se encaixa profundamente no bolso de ligação Mpro, ocupando os subbolsos S1 e S2 de forma mais eficaz do que o ligante de referência 38a, que se espalhou em direção ao bolso S1. O composto 101 267 741 formou 11 interações intermoleculares, sete ligações de hidrogênio e quatro contatos hidrofóbicos.
Este artigo apresenta um protocolo para a construção de um modelo de farmacóforo consenso mediante a integração de características moleculares de múltiplos ligantes, com foco específico na Mpro do SARS-CoV-2. Este método aprimora os esforços de descoberta de fármacos ao identificar características-chave de interação para triagem virtual e design racional de medicamentos.
A modelagem de farmacóforo consenso utilizando extensas bibliotecas de ligantes permite maior confiança preditiva na descoberta precoce de fármacos, particularmente para alvos como a Mpro do SARS-CoV-2. A integração de conformações ligante-associadas diversas reduz o viés do modelo e apoia triagem virtual robusta, impactando diretamente a identificação de compostos iniciais e a triagem de portfólios. Essa abordagem agiliza a seleção racional de candidatos para alvos com dados estruturais abundantes, aumentando a eficiência da P&D empresarial.
Este protocolo de modelagem de farmacoforo por consenso é adequado desde as fases iniciais da descoberta até a identificação de compostos líderes, aproveitando a biologia estrutural e a bioinformática para orientar triagens virtuais e a priorização de candidatos.