Evaluación de modelos por homología del receptor RDL de insectos para screening virtual: influencia del estado conformacional del templado en pLGICs

FELSZTYNA, IVÁN; VILLARREAL, MARCOS A.; GARCÍA, DANIEL A.; MIGUEL, VIRGINIA

Virtual

Simposio; 6to Simposio Argentino de Jóvenes Investigadores en Bioinformática; 2021

ICSB Regional Student Group Argentina

Los canales iónicos pentaméricos activados por ligando (pLGICs) constituyen una gran familia de receptores transmembrana ampliamente expresados en animales y en algunas bacterias. Esta familia incluye a los receptores del ácido γ-aminobutírico (GABA), entre otros. El homopentámero de la subunidad RDL es el principal receptor de GABA en el sistema nervioso de los insectos. Este receptor presenta diferencias estructurales con los receptores GABAA de vertebrados que le otorgan un comportamiento farmacológico particular. Por ello, el receptor RDL (R-RDL) es uno de los blancos más relevantes para la unión de insecticidas. Debido a las dificultades que implica la cristalización de pLGICs, se ha utilizado en numerosos estudios el modelado por homología para obtener estructuras de estas proteínas y realizar estudios computacionales de su unión a ligandos. Sin embargo, no ha sido estudiado en profundidad el impacto que podría tener el estado conformacional del templado (o estructura molde) en el desempeño del modelo obtenido para un screening virtual. En este trabajo nos proponemos obtener modelos por homología del R-RDL en distintos estados conformacionales y evaluar su desempeño en un screening virtual retrospectivo de insecticidas bloqueantes del canal. Para ello, se obtuvieron quince modelos del R-RDL basados en distintos templados de la familia pLGICs, cuyas estructuras representan los tres estados conformacionales conocidos para estos canales iónicos: cerrado, abierto y desensibilizado. Sobre estos modelos, se realizaron ensayos de docking molecular con un conjunto de ligandos activos e inactivos. Para evaluar el desempeño, se calcularon la curva ROC y los factores de enriquecimiento para cada modelo. Además, se realizaron simulaciones de dinámica molecular (SDM) de los mejores modelos para cada estado conformacional. Las estructuras iniciales se tomaron de las poses obtenidas para el insecticida fipronil, un compuesto canónico para el sitio de unión. Los parámetros de desempeño mostraron variaciones de acuerdo a los estados conformacionales de los modelos, particularmente en cuanto al enriquecimiento temprano de ligandos activos. Se evaluó la correlación de estos parámetros con distintas variables, con el objetivo de analizar cuáles fueron los factores determinantes para la correcta identificación de ligandos. Se observó que la identidad de secuencia no resultó relevante, mientras que algunas propiedades estructurales del canal, como el área y el volumen accesibles al solvente y el diámetro del poro a la altura de ciertos residuos, pueden dar cuenta de las variaciones en el desempeño. Los mejores resultados se obtuvieron con un modelo basado en un templado en estado conformacional cerrado. El análisis de las SDM confirma que las interacciones entre los residuos del sitio de unión y el insecticida fipronil coinciden con los resultados experimentales disponibles en el modelo cerrado, mientras que no se presentan de manera estable en los otros dos estados conformacionales. Tomando en conjunto estos resultados, puede destacarse que es necesario explorar diversos templados para obtener modelos por homología confiables en este receptor de la familia pLGICs, siendo particularmente relevante el estado conformacional del canal. El modelo que presentó las mejores métricas de desempeño podría ser utilizado en un screening virtual prospectivo de nuevos compuestos insecticidas.