Revista de Salud VIVE. Revista de Investigación en Salud
https://doi.org/10.33996/revistavive.v1i2.10
Volumen 1 | No. 2 | Mayo - Agosto 2018
http://revistavive.org
ISSN: 2664 – 3243
pp. 103 – 111
Modelado de interacciones del ADN con fármacos de actividad anti-cancerígena mediante el método ONIOM
Modeling of DNA interactions with anti-cancer activity drugs using the ONIOM method
Modelando interações de DNA com drogas anticancerígenas usando o método ONIOM
Fabiana Salazar
Universidad de Carabobo/ Instituto Venezolano de Investigaciones Científicas (IVIC), Venezuela.
Recibido enero 2018 / Revisión febrero 2018 / Aceptado 1 de mayo 2018
RESUMEN
Introducción: Los complejos metálicos juegan un importante papel en la lucha contra el cáncer; estos fármacos interaccionan con el ADN mediante enlaces covalentes con las bases nitrogenadas, provocando distorsión de la hebra para su posterior reconocimiento por proteínas que inducen la muerte celular. Objetivo: Caracterizar a través de métodos computacionales basados en Mecánica Cuántica y Mecánica Molecular, la interacción entre diversos fármacos anti-cancerígenos y el ADN. Materiales y métodos: Se optimizaron todas las estructuras de dichos fármacos (complejos de Pt(II), Ag(I) y Pd(II)) y las bases nitrogenadas del ADN (A, C, G y T), empleando DFT al nivel B3LYP y las bases atómicas 6-31++G(d,p) para los átomos livianos, mientras que el dodecaedro del ADN se optimizó mediante mecánica molecular empleando UFF como campo de fuerza. Resultados: A través del análisis de los mapas de ESP y Energía de Ionización Local Promedio, se determinó que el sitio más nucleofílico en las bases nitrogenadas correspondió al N1 de la adenina. La especie más electrofilia fue el cis-diamindiaquoplatino (II) con aumento del carácter electrofílico tras la sustitución de grupos salientes por ligandos agua. Los sistemas fármaco/base nitrogenada se caracterizaron a través de QTAIM, observándose una interacción de capa cerrada en todos los sistemas. Los sistemas de Pt y Pd con guanina resultaron ser los más estables. Conclusión: Se determinó que el método ONIOM fue capaz de modelar la torcedura que sufre la molécula de ADN al coordinarse con complejos platinados, por lo tanto, para fines de observar geometrías moleculares, resultó confiable.
Palabras clave: ADN; fármacos; anticancerígenos; cáncer; molecular
ABSTRACT
Introduction: Introduction: Metal complexes play an important role in the fight against cancer; these drugs interact with DNA through covalent bonds with nitrogenous bases, causing distortion of the strand for subsequent recognition by proteins that induce cell death. Objective: To characterize, through computational methods based on Quantum Mechanics and Molecular Mechanics, the interaction between various anti-cancer drugs and DNA. Materials and methods: All the structures of these drugs (complexes of Pt (II), Ag (I) and Pd (II)) and the nitrogen bases of DNA (A, C, G and T) were optimized, using DFT at the level B3LYP and the atomic bases 6-31 ++ G (d, p) for light atoms, while the DNA dodecahedron was optimized by molecular mechanics using UFF as a force field. Results: Through the analysis of the ESP and Average Local Ionization Energy maps, it was determined that the most nucleophilic site in the nitrogenous bases corresponded to the N1 of the adenine. The most electrophilia species was cis-diamindiaquoplatin (II) with an increase in electrophilic character after substitution of leaving groups with water ligands. Nitrogenous drug / base systems were characterized through QTAIM, with a closed layer interaction observed in all systems. The Pt and Pd systems with guanine proved to be the most stable. Conclusion: It was determined that the ONIOM method was able to model the twist that the DNA molecule undergoes when coordinating with platinum complexes, therefore, for the purpose of observing molecular geometries, it was reliable.
Key words: DNA; drugs; anticancer; cancer; molecular
RESUMO
Introdução: Complexos metálicos desempenham um papel importante na luta contra o câncer; essas drogas interagem com o DNA através de ligações covalentes com as bases de nitrogênio, causando distorção do fio para posterior reconhecimento por proteínas que induzem a morte celular. Objetivo: Caracterizar através de métodos computacionais baseados em Mecânica Quântica e Mecânica Molecular, a interação entre várias drogas anticancerígenas e DNA. Materiais e métodos: Otimizado todas as estruturas dessas drogas (Pt(II), complexos Ag(I) e Pd(II) e bases de nitrogênio de DNA (A, C, G e T), usando DFT a nível b3LYP e bases atômicas 6-31++G(d,p) para átomos de luz , enquanto o dodecaedro de DNA foi otimizado pela mecânica molecular usando uff como campo de força. Resultados: Através da análise dos mapas esp e energia da ionização local média, foi determinado que o local mais nucleílico dos bases de nitrogênio correspondiam ao N1 da adenina. A espécie mais eletrofilalfoi cis-diamindiaquoplatin (II) com maior caráter eletrofilico após substituição de grupos de saída por ligas de água. Os sistemas de nitrogênio droga/base foram caracterizados através da QTAIM, com uma interação de camada fechada observada em todos os sistemas. Os sistemas pt e PD com guanine acabaram por ser os mais estáveis. Conclusão: Foi determinado que o método ONIOM foi capaz de modelar a torção da molécula de DNA coordenando com complexos banhados a platina, portanto, com o propósito de observar geometrias moleculares, mostrou-se confiável.
Palavras-chave: DNA; drogas; anticâncer; câncer; molecular
INTRODUCCIÓN
|
L |
os métodos basados en Mecánica Cuántica (QM por sus siglas en inglés) definen el estado de un sistema a partir de la ecuación de Schrödinger, permitiendo aproximar la energía de una molécula con un tratamiento explícito de su estructura electrónica. Esto implica un alto costo computacional (en términos de consumo de tiempo y memoria) asociado al cálculo de términos matemáticos que cuantifican la interacción de los electrones, proporcionando una descripción precisa del comportamiento de los mismos y por lo tanto, de la estructura de los sistemas químicos (1).
Debido a que el costo computacional resulta ser relativamente bajo en comparación con los métodos mecano-cuánticos, es posible con los computadores que se disponen en la actualidad modelar sistemas de gran tamaño (de hasta 104 átomos). Como herramienta, la química computacional puede ser aplicable a otras áreas como la física, la biología y la bioquímica. Una aplicación específica de la química computacional ha sido en el modelaje de las interacciones y el comportamiento de los sistemas biológicos (2).
Por si solos, los métodos MM y QM parecen ser exiguos al modelarlos y esta carencia llevó a pensar en nuevos métodos computacionales. De esta premisa surgen los métodos híbridos (entre ellos el método ONIOM) capaces de aprovechar las ventajas de la mecánica cuántica y la mecánica molecular (3). Este método permite estudiar los sistemas complejos al fragmentarlos en dos o más regiones a los cuales se les aplican diferentes niveles de teoría. Un ejemplo de ello es el estudio de las interacciones de los ácidos nucleicos con complejos metálicos, donde a la región del sitio químicamente activo, se aplican métodos ab initio y en el resto del sistema métodos de mecánica molecular (3,4).
En consecuencia a lo expuesto, se plantea estudiar las propiedades electrónicas y termodinámicas de fármacos anticancerígenos, con el objetivo de modelar las interacciones entre estos y el ADN y analizar la influencia de los complejos metálicos en la estructura helicoidal.
El cáncer es un grupo de enfermedades caracterizadas por el crecimiento celular no regulado y la invasión y diseminación de células desde el sitio de origen a otros sitios en el cuerpo (metástasis) (5). Es una enfermedad genética, es decir, es causada por cambios en los genes que controlan el funcionamiento de las células, específicamente aquellos que controlan el crecimiento y la división. Los cambios pueden heredarse o adquirirse durante la vida de una persona, ya sea como resultado de errores en el ADN (Ácido desoxirribonucleico) o por exposición a sustancias carcinógenas. Algunos cambios afectan solo una unidad de ADN, lo que se llama un nucleótido, mientras que otros cambios comprenden tramos más grandes (6).
Entre las terapias anti-cancerígenas implementadas con mayor frecuencia se encuentra la quimioterapia, que debe su nombre al uso de sustancias químicas para tal fin. Estas sustancias suelen ser fármacos a base de complejos de metales de transición, mayormente complejos platino (Pt) como el cis-diamindicloroplatino (II) (7). El cis-diamindicloroplatino (II) (o cisplatino como comúnmente se le conoce), es un complejo de platino (II) con geometría plano cuadrada y cuyos ligandos (−Cl y −NH3) se encuentran en posición cis. Su uso clínico generalizado ha provocado el aumento de su toxicidad y resistencia en el organismo, por lo que hoy en día solo se usa en combinación con otros compuestos (7). Estos complejos de platino tienen como blanco de acción al ADN, con el cual interaccionan covalentemente a través del N7 de la guanina y formando productos conocidos como aductos (7). Así, las interacciones entre el ADN y los complejos de platino (II) son objeto de investigaciones que buscan comprender el efecto que estos fármacos ejercen sobre la estructura del ADN (8). En este sentido, la química computacional se ha convertido en una herramienta útil para su estudio, permitiendo el modelado y análisis de estos sistemas, y la comprensión de las interacciones moleculares que tienen lugar en los mismos (1).
Esto conlleva a tratar el sistema ADN/fármaco con un método multiescala, es decir, un método que combine una descripción cuántica (QM) en la región de interés (donde el fármaco se enlaza a las bases del ADN) con una descripción clásica (MM) en el resto del entorno (1). A este tratamiento se le conoce como método QM: MM (en este caso particular ONIOM) y se ha evidenciado un aumento en su uso para estudiar los procesos biológicos (4).
El incremento del cáncer en la población mundial ha llevado a organizaciones como la Organización Mundial de la Salud (OMS) al desarrollo de programas informativos sobre las causas, riesgos, prevención y tratamiento de esta enfermedad. Una de las medidas para el control del cáncer está enfocada a comprender el mecanismo de acción de los fármacos (6). El blanco de acción es diverso, pueden actuar sobre el ADN rompiendo la hélice, interfiriendo con las proteínas relacionadas con el ADN o modificando la expresión de genes específicos (9). Los complejos de Platino (II), por ejemplo, actúan directamente sobre la hélice del ADN provocando su distorsión y posteriormente la apoptosis (9). Este mecanismo de acción se ha convertido un punto clave para comprender la actividad biológica de los complejos metálicos (8).
El modelaje de los fármacos y sus interacciones mediante técnicas computacionales permite dilucidar propiedades moleculares (distribución de la densidad electrónica, energía de ionización, basicidad y acidez, nucleofilicidad) (9) y termodinámicas responsables de la afinidad global y especificidad de los mismos (9). Estudios combinados pueden realizarse a través del método ONIOM, el cual fragmenta el sistema en varias regiones a fin de estudiar cada región mediante distintos niveles de teoría (1). Las interacciones de los ácidos nucleicos con complejos metálicos han sido estudiadas y aplicadas exitosamente mediante el método ONIOM (9), por lo que se hará uso del mismo para estudiar las interacciones del ADN con complejos de Pt (II).
MATERIALES Y MÉTODOS
|
L |
a investigación se llevó a cabo empleando paquetes de software libre para el cálculo de propiedades moleculares y la construcción y visualización de las geometrías moleculares. Los cálculos ONIOM, QM y MM se realizaron mediante el programa Gaussian. Para la construcción de los modelos moleculares se usó el programa Avogadro versión 1.2.0, el cual es un editor y visualizador de moléculas avanzado diseñado para el uso multiplataforma en química computacional (10).
Sistema de estudio
El estudio se realizó sobre distintos fármacos anti-cancerígenos, las bases nitrogenadas del ADN (A, C, G y T), sistemas fármaco/base nitrogenada y el sistema ADN/fármaco. Los fármacos modelados fueron tres complejos de Pt (II) (cisplatino, carboplatino y oxaliplatino) y dos (2) complejos metal-NHC de Pd (II) y Ag (I) respectivamente. La hebra de ADN presentó un total de doce nucleótidos (dodecaedro) en la secuencia d (CCTCTGGTCTCC) y los fármacos de interés interaccionaron con ésta mediante el N7 de la guanina a través de enlaces intracadena (al ser estos los de mayor interés biológico). La secuencia de ADN empleada se tomó de un estudio RMN realizado por Gelaso y Lippard (11). Donde determinaron la estructura de la macromolécula al interaccionar con el cis-diaminplatino (II) a través de un enlace intracadena d (GpG). Los sistemas ADN/fármaco se estudiaron independientemente uno del otro mediante el método ONIOM (QM: MM) mientras que los sistemas fármaco/base nitrogenada se estudiaron a través cálculos mecano-cuántica a nivel DFT.
Archivos de entrada en Gaussian
Las estructuras moleculares de los sistemas de estudio construidas se ejecutaron con el programa GaussView versión 5.0. El cuadro de diálogo Gaussian Calculación Setup (Configuración de cálculos de Gaussian) proporciona acceso a todas las principales características de Gaussian (12).
Cálculos mecano-cuánticos (QM): optimización de geometría
Gaussian09 permite realizar cálculos mecano-cuántica a diversos niveles de teoría. Las funciones de polarización se eligieron para modelar la distribución electrónica entre átomos enlazados y las funciones difusas para considerar la densidad electrónica lejos de los núcleos que poseen pares solitarios de electrones (nitrógeno, fósforo y oxígeno) (13).
Finalmente se especifica el nombre del archivo wfn que se quiere generar. Gen permite especificar un conjunto de bases dado por el usuario (14). Opt permite que se realice una optimización geométrica. El algoritmo predeterminado es el Algoritmo de Berny usando gediis. MaxCycle establece el número máximo de pasos de optimización (15). Input hace que el conjunto de base actual se imprima en una forma adecuada para usar como entrada de conjunto de base general y, por lo tanto, puede usarse para agregar o modificar conjuntos de bases estándar (16).
Cálculos de mecánica molecular (MM): optimización de la geometría del ADN
La estructura molecular del ADN se exportó desde Avogadro en formato .mol para ser leído con GaussView y se optimizó mediante el módulo de mecánica molecular de Gaussian09 utilizando UFF como campo de fuerza. Uff indica uno de los tres campos de fuerza que Gaussian tiene disponible para cálculos de MM y ONIOM (17).
Cálculo ONIOM (QM: MM): optimización de la geometría del sistema ADN fármaco
Las geometrías optimizadas de los fármacos y el ADN se usaron para modelar los sistemas ADN/fármaco, en los cuales el centro metálico se unió covalententemente al ADN mediante un enlace d (GpG). Las geometrías se optimizaron a través del método ONIOM.
Cálculo de propiedades moleculares
El cálculo del Potencial Electroestático (ESP) y la Energía de Ionización Local Promedio se realizó mediante el programa Multiwfn versión 3.4.1 (18). Multiwfn realiza análisis de funciones de onda electrónicas de manera clara. El archivo .wfn se ejecutó generando dos archivos para cada propiedad: .xyz y .cub, y mediante el programa wfn to binary (19).
Cálculo QTAIM
La caracterización computacional de las interacciones entre los fármacos y sistemas, fármaco/base nitrogenada se realizó mediante análisis topológico de la densidad electrónica a través del programa AIM-UC versión 1.6.4 (20). AIM-UC es un programa que realiza algunos gráficos relacionados con QTAIM de Richard Bader a partir de un mallado 3D de densidad de carga electrónica (20).
RESULTADOS
|
A |
través del análisis de los mapas de ESP y Energía de Ionización Local Promedio, se determinó que el sitio más nucleofílico en las bases nitrogenadas correspondió al N1 de la adenina. La especie más electrofilia fue el cis-diamindiaquoplatino (II) con aumento del carácter electrofílico tras la sustitución de grupos salientes por ligandos agua. Los sistemas fármaco/base nitrogenada se caracterizaron a través de QTAIM, observándose una interacción de capa cerrada en todos los sistemas. Los sistemas de Pt y Pd con guanina resultaron ser los más estables. Se detalla a continuación:
