Collinearity of Decomposed Energy Terms in MM-GBSA Binding Free Energy Calculations

Resumo Técnico: Colinearidade de Termos de Energia Decompostos em Cálculos de Energia Livre de Ligação MM-GBSA

Definição do Problema
O método de Mecânica Molecular-Generalized Born Surface Area (MM-GBSA) é uma abordagem de estado final padrão para estimar energias livres de ligação em design de fármacos computacional. Uma prática comum envolve a decomposição da energia livre total em componentes distintos — van der Waals ($E_{vdW}$), eletrostática ($E_{elec}$), solvatação polar ($G_{GB}$) e solvatação não polar ($G_{SA}$) — para correlacionar esses termos com dados experimentais ou atribuir significados físicos específicos a eles. Este estudo contesta a suposição fundamental de que esses termos decompostos possuem coeficientes de ajuste independentes. Os autores postulam que tratar esses termos como preditores estatisticamente independentes em modelos de regressão pode ser inválido devido às redundâncias matemáticas e físicas inerentes ao formalismo MM-GBSA.

Metodologia
Os autores empregaram uma abordagem dupla combinando derivações analíticas e simulações de dinâmica molecular (MD) em larga escala para investigar as relações estatísticas entre os termos de energia.

• Derivação Analítica: O estudo examinou matematicamente as dependências estruturais entre os componentes de energia, focando especificamente em como o modelo Generalized Born (GB) é projetado para realizar a blindagem (screening) da eletrostática no vácuo.
• Dados de Simulação: A análise utilizou extensos dados de trajetórias de MD. Os autores avaliaram tanto as flutuações instantâneas dentro das trajetórias quanto valores médios através de um conjunto diversificado de 139 complexos proteína-proteína.
• Comparação de Protocolos: A robustez dos achados foi testada através de diferentes protocolos MM-GBSA, incluindo abordagens de trajetória única e de três trajetórias.
• Análise de Entropia: O estudo também incorporou a Entropia de Interação e correções de entropia $C_2$ para avaliar sua dependência em relação às flutuações eletrostáticas subjacentes.

Principais Contribuições e Resultados
A principal contribuição deste trabalho é o diagnóstico sistemático da severa colinearidade entre os termos de energia decompostos padrão no MM-GBSA. Os resultados demonstram:

1. Redundância entre Eletrostática e Solvatação Polar: A energia de interação Coulombiana proteína-ligante ($E_{elec}$) e a correção de solvatação GB ($G_{GB}$) exibem uma colinearidade quase perfeita ($R^2 \geq 0,99$). Isso reflete sua relação projetada, onde o termo GB serve como um mecanismo de blindagem do solvente para a eletrostática no vácuo.
2. Correlação entre Van der Waals e Solvatação Não Polar: Da mesma forma, a energia de interação de van der Waals e o termo de solvatação não polar ($G_{SA}$) mostram forte correlação, pois ambos são primariamente impulsionados pela área de superfície enterrada.
3. Dependência de Entropia: As correções de entropia de interação e $C_2$ foram encontradas fortemente dependentes das flutuações eletrostáticas subjacentes, reforçando ainda mais a redundância dentro da decomposição de energia.
4. Universalidade: Essas correlações persistem independentemente de os dados serem analisados ao nível de flutuações instantâneas de trajetória ou como médias de conjunto. Os achados são válidos para o conjunto diversificado de 139 complexos proteína-proteína e permanecem consistentes tanto em protocolos de trajetória única quanto de três trajetórias.

Significância e Alegações
O artigo afirma que a suposição de coeficientes de ajuste independentes para os termos decompostos do MM-GBSA é estatisticamente infundada. Consequentemente, os autores alertam contra o uso desses termos individuais como preditores independentes em modelos de regressão, pois fazê-lo pode levar ao sobreajuste (overfitting) ou à interpretação errônea das contribuições físicas.

Em vez disso, o estudo sugere uma estratégia de modelagem mais robusta: combinar os termos altamente correlacionados em contribuições efetivas e agregadas. Especificamente, os autores propõem agrupar os termos em contribuições polares, não polares e entrópicas efetivas. Ao reconhecer e abordar essa colinearidade, o estudo visa facilitar uma modelagem preditiva mais interpretável e estatisticamente robusta em cálculos de energia livre de ligação.