Método acoplado Autodock–PM6 para seleccionar la mejor pose en estudios de acoplamiento molecular

Autodock–PM6 method to choose the better pose in molecular docking studies

A Autodock–PM6 metodo pra seleccionar o melhor pose en estudios de acoplamiento molecular

Margarita Velásquez¹, Juan Drosos¹, Carlos Gueto¹, Johana Márquez¹, Ricardo Vivas–Reyes¹**

¹ Departamento de Química, Facultad de Ciencias, Universidad Nacional de Colombia, sede Bogotá, Crr. 30 No. 45–03, Bogotá DC, Colombia. Grupo de investigación Estado Sólido y Catálisis Ambiental (Esca).
*Autor de Correspondencia. Email: rvivasr@unicartagena.edu.co

Recibido: 15 de enero de 2013 • Aceptado: 29 de abril de 2013

Resumen

El acoplamiento molecular (conocido como docking) es una técnica de mecánica molecular ampliamente utilizada para predecir energías y modos de enlace entre ligandos y proteínas, información de gran utilidad en el estudio de nuevos compuestos con efectos terapéuticos. No obstante, los resultados obtenidos mediante esta técnica tienden a la subjetividad, debido a que los programas utilizados para llevarla a cabo proporcionan más de un criterio de selección de la mejor pose. En la presente investigación se aplicó el método semiempírico PM6 a los resultados del acoplamiento, obteniendo con ello mejorías en el proceso de selección de la mejor pose pues se obtuvieron poses con alta probabilidad de unión al sitio activo de su receptor y con energías de unión menores a las reportadas por los criterios de selección ofrecidos por el programa de docking.

Palabras Clave: acoplamiento molecular, método semiempírico PM6, autodock, pose.

Abstract

The acoplamiento molecular (conocido como docking) es una técnica de mecánica molecular ampliamente utilizada para predecir energías y modos de enlace entre ligandos y proteínas, lo que proporciona información de gran utilidad para el estudio de nuevos compuestos con efectos terapéuticos. No obstante, los resultados obtenidos mediante esta técnica tienden a la subjetividad, debido a que los programas utilizados para llevarla a cabo proporcionan más de un criterio de selección de la mejor pose. En la presente investigación, se aplicó el método semiempírico PM6 a los resultados del acoplamiento, obteniendo con ello mejorías en el proceso de selección de la mejor pose al observarse finalmente poses con alta probabilidad de unión al sitio activo de su receptor y con energías de unión menores a las reportadas por los criterios de selección ofrecidos por el programa de docking.

Keywords: Molecular Coupling, PM6 semi-empirical method, AutoDock, Pose.

Resumo

El acoplamiento molecular (conocido como docking) es una técnica de mecánica molecular ampliamente utilizada para predecir energías y modos de enlace entre ligandos y proteínas, lo que proporciona información de gran utilidad para el estudio de nuevos compuestos con efectos terapéuticos. No obstante, los resultados obtenidos mediante esta técnica tienden a la subjetividad, debido a que los programas utilizados para llevarla a cabo proporcionan más de un criterio de selección de la mejor pose. En la presente investigación, se aplicó el método semiempírico PM6 a los resultados del acoplamiento, obteniendo con ello mejorías en el proceso de selección de la mejor pose al observarse finalmente poses con alta probabilidad de unión al sitio activo de su receptor y con energías de unión menores a las reportadas por los criterios de selección ofrecidos por el programa de docking.

Palavras-chave: acoplamiento molecular, método semiempírico PM6, autodock, pose.

Introducción

Anteriormente, la búsqueda de nuevos compuestos con efectos terapéuticos implicaba la síntesis a través de varias rutas de los supuestos compuestos con la propiedad buscada o el cribado de productos naturales, esfuerzos que pueden ser categorizados como exploraciones aleatorias o al azar, en lugar de orientaciones racionales. En los últimos años se han dado a conocer una serie de grandes avances en los cálculos teóricos que han favorecido el uso y la popularidad de los estudios in silico (computacionales) para el análisis de moléculas de interés biológico; ello se ha convertido en parte integral de la investigación industrial y académica dirigida al diseño y descubrimiento de fármacos (1).

Las herramientas computacionales han evolucionando de tal forma que se han transformando en tecnologías cada vez más importantes para la búsqueda de moléculas candidatas a fármacos, mediante la se lección de moléculas denominadas 'cabezas de serie' o leads a partir de bases de datos. Entre estas herramientas se encuentra el modelado molecular, técnica que permite obtener moléculas –reales o virtuales– que poseen gran probabilidad de exhibir una acción específica y que hace posible, incluso, predecir su biodisponibilidad y toxicidad (1, 2).

A su vez, dentro de las herramientas de modelado molecular se encuentra el acoplamiento molecular (docking), técnica de gran utilidad para predecir la estructura de los complejos intermoleculares que se establecen entre dos o más moléculas. Uno de los casos más estudiados es la interacción proteína–ligando, cuyo objetivo principal es predecir energías y modos de enlace, y que suele hacerse previamente a un estudio experimental (3, 6).

En estudios donde se aplica el acoplamiento molecular, usualmente se considera la mejor pose aquella que reporta la menor energía de unión. Alternativamente, puede ser seleccionada del cluster más poblado entre los clusters obtenidos después del acoplamiento molecular y, en ocasiones, también se selecciona el menor valor de RMSD (Root Mean Square Deviation) con respecto a una estructura de referencia (4, 7).

Aunque los parámetros energéticos son muy importantes al momento de efectuar esta selección, es conveniente tener en cuenta que los campos de fuerza empleados por los métodos de mecánica molecular (MM) (entre éstos, el acoplamiento molecular) llevan consigo algunas aproximaciones que lógicamente conducirían a resultados no tan reales. Estas aproximaciones se relacionan principalmente con el hecho de considerar sólo la disposición nuclear de los átomos, prescindiendo de los electrones, al asumir que éstos se distribuirán óptimamente a su alrededor.

Persiguiendo este objetivo, en un estudio previo (8) realizado por nuestro grupo de investigación, se utilizó el método semiempírico de mecánica cuántica PM6 (9, 10) con el cual se recuperaría de los estudios de docking molecular la información perdida en términos de interacciones y energía de acoplamiento. Para tal estudio, la metodología arrojó resultados satisfactorios (9).

Este hecho se utilizó como punto de partida para el desarrollo del presente artículo, en el cual se evaluó la efectividad que tiene la aplicación del método semiempírico PM6 como complemento a los estudios de docking, utilizando para ello un grupo de complejos pertenecientes a la serie Astex Diverse (11) diseñada para mejorar el rendimiento y la exactitud de los programas de docking y que se distingue por contener estructuras modernas (archivos PDB depositados después del 11 de agosto del 2000) con resolución de 2.5 Å o menor y con ligandos similares a fármacos.

Finalmente, el uso de esta metodología acoplada permite obtener, poses con alta probabilidad de unión al sitio activo de su receptor para cada ligando y con energías de unión menores a las reportadas por los criterios de selección ofrecidos por el programa de docking.

Materiales y métodos

Selección de los complejos proteína–ligando

Para el desarrollo de este trabajo se utilizaron 22 complejos proteína–ligando (Tabla 1) tomados de la base de datos Protein Data Bank (PDB) (12) y seleccionados de la serie de prueba Astex Diverse (11). Esta serie está constituida por 85 complejos proteína–ligando que cumplen principalmente con los siguientes criterios: las dianas incluidas son estructuralmente diversas, aproximadamente el 90% de estas son utilizadas como blancos directos en proyectos de descubrimiento de fármacos y el resto están dentro de las llamadas moléculas off target, entre las cuales se encuentran los citocromos P450. Por otra parte, la mayoría de ligandos clasifican, ya sea como drogas aprobadas, en etapa de ensayos clínicos, o compuestos obtenidos en proyectos de descubrimiento de fármacos. Por último, se utilizaron estructuras de alta calidad para la cuales la densidad electrónica del ligando es compatible con el modo de unión (11).

Evaluación del programa de docking

Para el acoplamiento se utilizó el programa AutoDock v4.2 (13, 14). Se realizó un docking ciego entre cada uno de los 22 complejos proteína– ligando y su respectivo receptor, solicitándose para tal fin la generación de 100 poses en cada caso.

Posteriormente, se realizó una inspección visual de los resultados en el programa AutoDockTools v1.5.4 (15) (ADT), a fin de observar el grado de similitud existente entre las poses generadas y las poses nativas.

Preparación de receptores y ligandos

La preparación de los receptores y ligandos se llevó a cabo utilizando el programa Discovery Studio 2.5 (16).

El tratamiento de los receptores consistió en la extracción del ligando y la eliminación de moléculas de agua y cofactores con los cuales vienen resueltas sus estructuras cristalinas.

Por su parte, el tratamiento de los ligandos se basó en la revisión de las posiciones de sus átomos y de todos sus enlaces, tomando como referencia la estructura 2D mostrada en la base de datos PDB. En este punto también fueron adicionados los átomos de hidrógeno a las estructuras.

Acoplamiento molecular

El acoplamiento se realizó bajo los siguientes parámetros: libre rotación de todos los enlaces de los ligandos con libertad conformacional, adición de cargas de Gasteiger a la proteína y al ligando, adición de hidrógenos a la proteína, generación de 100 poses para cada ligando y utilización del Algoritmo Genético Lamarckiano (LGA). Este último se utilizó bajo las siguientes condiciones: 100 corridas del GA, 150 tamaño de población, 2 500 000 máximo número de evaluaciones, 27 000 máximo número de generaciones, 0.02 tasa de mutación del gen y 0.08 rata de cruza.

Selección de las mejores poses

Para seleccionar la mejor pose para cada ligando se utilizaron los siguientes criterios:

Criterio 1. Energía de unión del complejo: se toma como mejor pose aquella que reporta la menor energía de unión entre todas las poses obtenidas.

Criterio 2. Distribución de energía (cluster): se toma como mejor pose aquella que ocupa la primera posición dentro del clúster más poblado.

Recategorización de los resultados del docking utilizando el método semiempírico PM6

Acoplamiento molecular

Los resultados obtenidos de los cálculos de acoplamiento molecular fueron categorizados de acuerdo con los criterios 1 y 2, previamente mencionados (Tablas 2 y 3).

. Este porcentaje está representado por los ligandos DBQ, IH5, TQ3 y PFP, para los cuales la mejor pose seleccionada no coincidió con la pose obtenida mediante el criterio1 (Tabla 4).

Al observar en detalle y comparar –usando el programa AutoDockTools (ADT)–, la pose nativa vs. la pose obtenida para cada ligando, se observó gran similitud conformacional entre éstas para la mayor parte de los casos, lo cual fue respaldado por los valores de RMSD(Root Mean Square Deviation) arrojados para dichas poses (Tablas 2 y 3 y Figura 1).

No obstante, se encontraron casos particulares en los cuales no se observó dicha similitud (ligandos PFP y ZZ2) (Figura 2 y Tablas 2 y 3). En estos casos, las orientaciones espaciales de los átomos de la mejor pose difieren notablemente de los átomos de la pose nativa, hecho que lleva a pensar que este tipo de criterio en cierta medida no ofrece una selección 100% confiable.

Al aplicar el criterio de distribución de la energía (cluster) (Tabla 3) y comparar, se hizo evidente que aproximadamente en el 82% de los casos las poses coincidían con las obtenidas mediante el criterio 1 (Tabla 2), lo que era de esperarse ya que, generalmente, en el cluster más poblado se encuentran agrupadas las poses de menor energía. Sin embargo, es necesario fijarse en el 18% de los casos restantes en los cuales no se presentaron coincidencias