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

Resumen Técnico: Colinealidad de los Términos de Energía Descompuestos en los Cálculos de la Energía Libre de Unión MM-GBSA

Planteamiento del Problema
El método de Áreas de Superficie de Born Generalizada de Mecánica Molecular (MM-GBSA) es un enfoque de estado final estándar para estimar las energías libres de unión en el diseño computacional de fármacos. Una práctica común consiste en descomponer la energía libre total en componentes distintos —van der Waals ($E_{vdW}$), electrostática ($E_{elec}$), solvatación polar ($G_{GB}$) y solvatación no polar ($G_{SA}$)— para correlacionar estos términos con datos experimentales o para asignarles significados físicos específicos. Este estudio cuestiona la suposición fundamental de que estos términos descompuestos poseen coeficientes de ajuste independientes. Los autores postulan que tratar estos términos como predictores estadísticamente independientes en modelos de regresión puede ser inválido debido a las redundancias matemáticas y físicas inherentes al formalismo de MM-GBSA.

Metodología
Los autores emplearon un enfoque dual que combina derivaciones analíticas y simulaciones de dinámica molecular (MD) a gran escala para investigar las relaciones estadísticas entre los términos de energía.

• Derivación Analítica: El estudio examinó matemáticamente las dependencias estructurales entre los componentes de energía, centrándose específicamente en cómo el modelo de Born Generalizada (GB) está diseñado para el apantallamiento de la electrostática en vacío.
• Datos de Simulación: El análisis utilizó extensos datos de trayectorias de MD. Los autores evaluaron tanto las fluctuaciones instantáneas dentro de las trayectorias como los valores promedio a través de un conjunto diverso de 139 complejos proteína-proteína.
• Comparación de Protocolos: La robustez de los hallazgos se probó en diferentes protocolos de MM-GBSA, incluyendo enfoques de trayectoria única y de tres trayectorias.
• Análisis de Entropía: El estudio también incorporó la Entropía de Interacción y las correcciones de entropía $C_2$ para evaluar su dependencia de las fluctuaciones electrostáticas subyacentes.

Contribuciones Clave y Resultados
La contribución principal de este trabajo es el diagnóstico sistemático de una severa colinealidad entre los términos de energía descompuestos estándar en MM-GBSA. Los resultados demuestran:

1. Redundancia entre Electrostática y Solvatación Polar: La energía de interacción Coulomb entre proteína y ligando ($E_{elec}$) y la corrección de solvatación GB ($G_{GB}$) exhiben una colinealidad casi perfecta ($R^2 \geq 0.99$). Esto refleja su relación diseñada, donde el término GB sirve como un mecanismo de apantallamiento de solvente para la electrostática en vacío.
2. Correlación entre Van der Waals y Solvatación No Polar: Del mismo modo, la energía de interacción de van der Waals y el término de solvatación no polar ($G_{SA}$) muestran una fuerte correlación, ya que ambos son impulsados primordialmente por el área de superficie enterrada.
3. Dependencia de la Entropía: Se encontró que la entropía de interacción y las correcciones de entropía $C_2$ dependen fuertemente de las fluctuaciones electrostáticas subyacentes, reforzando aún más la redundancia dentro de la descomposición de la energía.
4. Universalidad: Estas correlaciones persisten independientemente de si los datos se analizan al nivel de las fluctuaciones instantáneas de la trayectoria o como promedios de conjunto. Los hallazgos son válidos para el diverso conjunto de 139 complejos proteína-proteína y se mantienen consistentes tanto en los protocolos de una como de tres trayectorias.

Significancia y Reivindicaciones
El artículo sostiene que la suposición de coeficientes de ajuste independientes para los términos descompuestos de MM-GBSA es estadísticamente errónea. En consecuencia, los autores advierten contra el uso de estos términos individuales como predictores independientes en modelos de regresión, ya que hacerlo puede conducir al sobreajuste o a la mala interpretación de las contribuciones físicas.

En su lugar, el estudio sugiere una estrategia de modelado más robusta: combinar los términos altamente correlacionados en contribuciones efectivas y agregadas. Específicamente, los autores proponen agrupar los términos en contribuciones efectivas polares, no polares y entrópicas. Al reconocer y abordar esta colinealidad, el estudio pretende facilitar un modelado predictivo más interpretable y estadísticamente robusto en los cálculos de la energía libre de unión.