Multi-predictor pathogenicity filtering and membrane-specific stability analysis identify and characterize candidate destabilizing missense variants of the GIP receptor (GIPR)

Resumo Técnico: Filtragem de Patogenicidade Multi-preditor e Análise de Estabilidade Específica de Membrana de Variantes de GIPR

Definição do Problema
O receptor do polipeptídeo insulinotrópico dependente de glicose (GIPR) é um alvo terapêutico validado para o diabetes tipo 2 e obesidade, mas as consequências funcionais de suas variantes missense permanecem amplamente não caracterizadas. Existe uma lacuna metodológica significativa na análise de variantes: a maioria dos preditores de patogenicidade e estabilidade foi desenvolvida para proteínas solúveis e globulares, apresentando desempenho insatisfatório quando aplicados a proteínas integrais de membrana, como o GIPR. Consequentemente, há uma carência de triagens abrangentes e conscientes da estrutura que possam distinguir variantes funcionalmente críticas de variantes benignas dentro do ambiente da bicamada lipídica.

Metodologia
Os autores empregaram um pipeline in silico rigoroso e de múltiplos estágios para priorizar variantes de GIPR missense que provavelmente sejam tanto patogênicas quanto estruturalmente desestabilizadoras. O fluxo de trabalho procedeu da seguinte forma:

1. Recuperação e Expansão de Variantes: 763 rsIDs únicos de missense foram recuperados do dbSNP. Após a expansão de posições multi-alélicas, isso resultou em 893 variantes de nucleotídeo único (SNVs).
2. Filtragem Sequencial de Patogenicidade: As variantes foram filtradas através de um painel de preditores complementares para garantir alta especificidade:
   • PolyPhen-2: Aplicado com um limiar de probabilidade estrito de ≥ 0,9 (modelo HumDiv).
   • Filtro de Raridade: Variantes com uma frequência alélica menor (MAF) global ≥ 0,005 no gnomAD v4.1.1 foram excluídas.
   • Conservação e Preditores de Ensemble: As variantes passaram pelo SIFT 4G (escore < 0,05), PROVEAN (escore ≤ −2,5), AlphaMissense (escore ≥ 0,564), MetaRNN (escore ≥ 0,5) e REVEL (escore ≥ 0,5).
   • Patogenicidade Específica de Membrana: O filtro de patogenicidade final foi o mCSM-membrane, um preditor treinado especificamente em proteínas transmembranares, que classificou as variantes como patogênicas ou benignas.
3. Análise de Estabilidade: Para as variantes que passaram pelos filtros de patogenicidade, os efeitos de estabilidade foram avaliados usando:
   • Módulo de Estabilidade mCSM-membrane: Utilizado como o principal filtro quantitativo com um limiar de ΔΔG ≤ −1,0 kcal/mol (indicando desestabilização).
   • DynaMut2: Utilizado para corroborar a direção do efeito (sinal de ΔΔG), reconhecendo seu treinamento em proteínas solúveis.
   • Missense3D-TM: Uma segunda ferramenta específica de membrana usada para identificar violações estereoquímicas específicas (ex: cargas enterradas, prolinas quebra-hélice).
4. Caracterização Estrutural: O conjunto final de variantes foi mapeado no modelo AlphaFold de comprimento total (AF-P48546-F1-v6) e na estrutura de criomicroscopia eletrônica (PDB: 7RBT). A conservação evolutiva foi avaliada via phyloP, GERP++ RS e ConSurf. O servidor HOPE foi utilizado para caracterizar as consequências estruturais, incluindo mudanças em tamanho, carga, hidrofobicidade e contatos moleculares.

Resultos Principais
O pipeline reduziu o conjunto inicial de 893 SNVs para um conjunto final de 14 variantes raras que atenderam a todos os critérios de patogenicidade e desestabilização estrutural.

• Distribuição: Todas as 14 variantes mapearam exclusivamente para as sete hélices transmembranares (TM), abrangendo seis das sete hélices (TM1, TM2, TM3, TM5, TM6, TM7). Nenhuma variante foi encontrada em TM4 ou nos loops extracelulares/intracelulares.
• Conservação: Todas as 14 posições mostraram forte restrição evolutiva, com escores phyloP elevados (4,55–9,55), escores GERP++ RS (>4,0) e graus ConSurf (7–9).
• Predições de Estabilidade:
  • mCSM-membrane previu todas as 14 variantes como desestabilizadoras (intervalo de ΔΔG: −1,000 a −1,716 kcal/mol).
  • DynaMut2 concordou com a direção da desestabilização para 10 variantes; foi discordante (prevendo neutro/estabilização) para quatro variantes, três das quais eram substituições de Glicina para Valina.
  • Missense3D-TM classificou 13 das 14 variantes como prejudiciais, com a única exceção sendo S179F.
• Mecanismos Estruturais: A análise HOPE identificou cinco mecanismos desestabilizadores recorrentes:
  1. Perda de Empacotamento do Núcleo/Formação de Cavidade: Y141N (a mais desestabilizadora, ΔΔG = −1,716 kcal/mol) remove uma tirosina enterrada, perdendo uma ligação de hidrogênio com Asp191 e criando uma cavidade.
  2. Substituição de Glicina Enterrada: G144V, G238V e G38V substituem glicinas pequenas e flexíveis por valinas volumosas em interfaces de hélice densamente compactadas, perturbando a geometria.
  3. Introdução de Cargas Enterradas ou Voltadas para o Lipídio: Cinco variantes (L152R, L241R, L318R, P302R, L344R) introduzem resíduos de arginina no núcleo hidrofóbico ou superfícies voltadas para o lipídio, incorrendo em penalidades eletrostáticas significativas.
  4. Disrupção da Geometria da Hélice: Q224P introduz uma prolina na TM3, prevista para quebrar a α-hélice e interromper a interface de ligação do peptídeo. P302L, P302R e P348H removem prolinas nativas, alterando a rigidez do esqueleto.
  5. Contato de Modulador: Três variantes (L344Q, L344R, P348H) contatam um modulador alostérico de pequena molécula ligado (LSN3556672), sugerindo potenciais implicações farmacogenéticas.

Significância e Alegações
O artigo afirma fornecer a primeira triagem abrangente e consciente da estrutura da variação missense de GIPR que considera especificamente o ambiente de membrana. Ao combinar preditores de patogenicidade ortogonais com ferramentas de estabilidade específicas de membrana, os autores geraram um "conjunto priorizado e mecanisticamente anotado" de 14 variantes candidatas.

Os autores enfatizam que esses resultados servem como uma base para futuros estudos funcionais e farmacogenéticos, em vez de uma classificação clínica definitiva. O estudo destaca que, embora as variantes sejam muito raras para serem determinantes comuns de doenças individualmente, elas revelam posições críticas essenciais para a integridade do receptor. Especificamente, a identificação de variantes próximas ao sítio de ligação do GIP (Q224P) e ao bolso do modulador alostérico sugere que esses candidatos são alvos de alta prioridade para validação experimental em relação à sinalização de incretina e resposta a drogas. O trabalho ressalta a necessidade de usar preditores específicos de membrana para evitar falsos negativos e imprecisões associadas à aplicação de ferramentas de proteínas solúveis em receptores transmembranares.