Introducción
El virus SARS-CoV-2 sorprendió a la humanidad por su facilidad de transmisión, cuyo mecanismo involucraba la difusión en el aire de pequeñas gotas expulsadas por las personas infectadas al respirar o toser. Este virus altamente infeccioso fue la causa principal de la enfermedad respiratoria denominada COVID-19 y en enero de 2020 fue declarada emergencia de salud pública por la Organización Mundial de la Salud (OMS) [1], [2]. La estructura macromolecular del SARS-CoV-2 es compleja debido a sus diferentes componentes, que incluyen una envoltura lipídica, ARN y ciertas proteínas [3], [4]. Estas proteínas son las responsables de unirse al receptor ACE2 de las células humanas, lo cual provoca que el virus entre e infecte al huésped. La enzima que desempeña un papel crucial en la replicación del virus SARS-CoV-2 es la proteasa 3CL [5]. Esta proteasa facilita la hidrólisis de los enlaces peptídicos en más de 11 sitios de la poliproteína viral; se destacan los residuos HIS41 y CYS145 localizados en el sitio activo de la proteasa 3CL [5], [6]. Ambos residuos son críticos para la actividad catalítica de la proteasa, por lo que pueden ser utilizados como la diana catalítica en el tratamiento del SARS-CoV-2 [7].
El acoplamiento molecular (molecular docking en inglés) es una técnica de simulación molecular que ha permitido realizar estudios de la interacción molecular de compuestos químicos con sitios activos de proteínas específicas para determinar sus propiedades farmacéuticas [8], [9]. La interacción molecular efectiva entre un fármaco y un sitio activo tiene un impacto directo en la disminución, la replicación y proliferación de virus como el SARS-CoV-2. En el desarrollo de la pandemia, el tratamiento médico con productos naturales fue utilizado ampliamente como mecanismo preventivo del SARS-CoV-2. Sin embargo, en muchos casos, la identificación de los compuestos activos que producían la inhibición de la actividad de la proteína proteasa 3CL no fue realizada. En Latinoamérica, la planta de tomillo común (Thymus vulgaris) fue utilizada para el tratamiento del COVID-19. Esta planta contiene un aceite esencial formado por compuestos orgánicos mezclados con propiedades antibacterianas, antifúngicas y antivíricas [10]. Para este aceite esencial obtenido del tomillo, el timol (2-isopropil-5-metilfenol) corresponde a una de las especies más importantes presentes en este aceite, con el cual se logró una recuperación favorable de pacientes con el padecimiento del COVID-19 (ver Figura 1). En otros estudios, la actividad de este compuesto fue evaluada y comparada con otros productos naturales como el carvacrol y la bis-demetoxi-curcumina [11]-[13]. En estos estudios, el timol mostró una buena actividad inhibidora, la cual estuvo asociada a su efectivo enlazamiento con la transmembrana proteasa serina 2 [12].
En el contexto planteado, para este trabajo se determinaron las energías de interacción molecular entre el timol y los residuos de aminoácidos presentes en el sitio activo de la proteasa 3CL mediante acoplamiento molecular. A la vez, usando la teoría de átomos en moléculas (QTAIM por sus siglas en inglés) en combinación con la teoría de interacciones no covalentes (NCI por sus siglas en inglés), se analizaron los tipos de interacciones moleculares presentes entre los residuos de aminoácidos identificados y el compuesto orgánico 2-isopropil-5-metilfenol. Así mismo, otros estudios de acoplamiento molecular han sugerido que este compuesto químico puede interactuar con el sitio activo de la proteasa 3CL, especialmente con los residuos HIS41 y CYS145 [14]-[16]. Por lo tanto, la evaluación de los tipos de interacciones moleculares presentes entre estos residuos de aminoácidos y el timol es de gran importancia farmacéutica para comprender a nivel molecular el mecanismo de inhibición de este sitio activo presente en la proteína proteasa 3CL.
Materiales y métodos
Preparación de la proteína proteasa 3CLy el compuesto timol (2-isopropil-5-metilfenol)
La estructura molecular de la proteasa 3CL fue obtenida de la base de datos RSCB (Research Collaboratory for Structural Bioinformatics, PBD ID: 6LU7) en formato Protein Data Bank (PDB) [4]. Para preparar la proteasa 3CL para el acoplamiento molecular, se eliminaron las moléculas de agua y el inhibidor N3 utilizando el programa Chimera-1.16 [17]. La estructura resultante se guardó en formato PDB y se utilizó para realizar el acoplamiento molecular con el programa AutoDock Tools-1.5.6 [18]. Por otro lado, la estructura molecular del timol se obtuvo de la base de datos PubChem (ht-tps://pubchem.ncbi.nlm.nih.gov/) [19]-[21] en formato SDF. Para optimizar la estructura molecular del timol, se utilizó el programa Avogradro-1.2 [22] empleando el campo de fuerza MMFF94 [23]. Finalmente, las geometrías optimizadas se guardaron en formato *.mol2 para su uso en el acoplamiento molecular.
Acoplamiento molecular
Para el acoplamiento molecular se utilizaron los programas Auto-Dock-4.2.6 y AutoDock Tools-1.5.6 [18]. Se agregaron los hidrógenos esenciales y las cargas de Kollman a la proteína, y las cargas de Gasteiger al ligando utilizando AutoDock Tools-1.5.6 [18]. El sitio utilizado como objetivo para el acoplamiento molecular fue el que contiene los residuos HIS41 y CYS145 de la proteasa 3CL. Se colocó una celda de 40 Å de longitud usando las coordenadas XYZ: -10,977; 16,402; 67,709. Además, se utilizaron los parámetros del algoritmo genético [24] para realizar el acoplamiento molecular, utilizando 50 conformaciones y repitiendo el acoplamiento por triplicado. La energía de interacción molecular y la constante de inhibición se determinaron a partir de la conformación con la menor energía del complejo proteína-ligando.
Cálculos cuánticos
Todos los cálculos cuánticos fueron realizados utilizando el programa Orca-5.0.3 [25], [26]. Los cálculos de energía se hicieron con el funcional M06-2X//aug-cc-pVDZ en fase gaseosa para la molécula de timol con los residuos de la proteasa 3CLpro (HIS41 y CYS145). Además, se evaluaron las interacciones moleculares entre el compuesto y los residuos utilizando la Teoría Cuántica de Átomos en Moléculas (QTAIM) [27], [28], la función de localización de electrones (ELF) [29], el método de gradiente reducido de la densidad (RDG por sus siglas en inglés) [30], [31] y el modelo de gradiente independiente basado en la partición de Hirshfeld (IGMH) [32]. Por último, los cálculos de QTAIM, ELF, RDG e IGMH se realizaron con el programa Multiwfn-3.8 [33] empleando los archivos *.wfn generados con el módulo orca_2aim.
Resultados y discusión
Acoplamiento molecular
El acoplamiento molecular se ha convertido en una herramienta valiosa para investigar la interacción entre una proteína y un compuesto denominado ligando, lo que la hace útil en el descubrimiento y diseño de fármacos [8], [9]. En este estudio se utilizó el método de acoplamiento molecular para evaluar la interacción entre la proteasa 3CL del SARS-CoV-2 y el timol, con el objetivo de determinar la energía de unión del complejo proteína-ligando y su constante de inhibición. Los resultados del estudio muestran que el timol se une con la región activa de la proteasa con una energía de unión de -15,833 (0,017) kJ/mol y una constante de inhibición de 1,29 (0,06) mM. Además, se observa la interacción molecular del timol con seis residuos de aminoácidos de la proteína: HIS41, LEU141, GLY143, SER144, CYS145 y HIS163 (ver Figura 3) mediante interacciones tipo enlaces de hidrógeno e interacciones hidrofóbicas.
Desde el punto de vista energético, la interacción molecular del ti-mol en el sitio activo es más favorable comparado con compuestos sintéticos como el Paritaprevir y el Simeprevir que mostraron valor de energía de interacción de -8,80 kJ/mol y 8,78 kJ/mol, respectivamente [16]. Del mismo modo, el timol muestra un valor de constante de inhibición mayor a los reportados para estos compuestos (0,36 |UM para el Paritaprevir y 0,37 ßM). Sin embargo, estos valores sugieren una mayor concentración de timol para inhibir el sitio activo.
Los residuos HIS41, LEU141, GLY143, SER144, CYS145 y HIS163 corresponden a los aminoácidos que conforman uno de los sitios activos de la proteasa 3CL del SARS-CoV-2. La presencia de los residuos HIS41, GLY143 y HIS163 es esencial para la actividad enzimática de la proteasa, mientras que el residuo LEU141 contribuye a la especificidad de los sustratos y la eficacia de los inhibidores [2], [34], [35]. Por su parte, el residuo SER144 actúa como nucleó-filo en la reacción de hidrólisis de los sustratos. Por último, el residuo CYS145 forma un enlace disulfuro con otro residuo de cisteína para estabilizar la estructura de la proteasa [36]. En general, estos residuos desempeñan un papel fundamental en la replicación del virus y son importantes para el desarrollo de terapias antivirales. Al observar la figura 3 podemos ver la formación de un enlace de hidrógeno entre el oxígeno del timol y el residuo CYS145 (distancia igual a 1,913 Å), lo que sugiere que este ligando es específico para la inhibición de la proteasa 3CLpro del SARS-CoV-2.
Análisis QTAIM
La Teoría Cuántica de Átomos en Moléculas (QTAIM) desarrollada por Bader [27], [28] es una herramienta muy útil para comprender la estructura química y la reactividad a partir de la densidad electrónica (p(r)), y también es utilizada para el análisis de interacciones no covalentes. La existencia de puntos críticos de enlace (BCP) entre los residuos HIS41 y CYS145 de la proteasa 3CLpro con el timol sugiere la presencia de varias interacciones moleculares (ver Figura 4). Entre estas interacciones, las más destacadas son los enlaces de hidrógenos y las interacciones dispersivas tipo Van der Waals. La Tabla 1 muestra los parámetros topológicos evaluados. El signo del valor del Laplaciano de la densidad electrónica (ï72 p(r)) permite evaluar el comportamiento de la densidad electrónica en el punto crítico de enlace, y la densidad de energía total (H(r)) permite evaluar la estabilidad de la interacción molecular [37], [38]. Estos valores son reportados en unidades atómicas (u.a.).
El análisis de los puntos críticos de enlace (BCP) permite identificar la naturaleza de la interacción molecular. Inicialmente, entre el timol y el residuo HIS41 se observa la formación de un punto crítico de enlace (BCP: 98) con un valor del Laplaciano electrónico de 0,0601 u.a. y una densidad de energía total igual a -0,0004 u.a. El Laplaciano electrónico al ser positivo indica una acumulación de densidad electrónica en el punto crítico, la cual es típica de enlaces covalentes. A su vez, este valor de H(r) indica una interacción molecular débil con una densidad electrónica a lo largo de la ruta de enlace baja (ver Figura 5) [28], [39], [40]. Esto sugiere que la interacción molecular, aunque presenta características covalentes, se asemeja más a una interacción molecular hidrógeno-hidrógeno muy débil entre estas moléculas [41].
Por su parte, en la interacción molecular entre el timol y el residuo CYS145 se observan cuatro puntos críticos de enlace (BCP: 76, 111, 113 y 123). Todos estos puntos críticos presentan valores positivos de Laplaciano (ÍVp(r) > 0), lo cual es indicativo de una dispersión de la densidad electrónica a lo largo de la ruta de enlace. Adicionalmente, los valores de energía total (H(r) > 0) positivos son característicos de interacciones de capa cerrada (closed-shell en inglés), es decir, interacciones tipo enlaces iónicos, enlaces de hidrógenos y Van der Waals. Los valores del Laplaciano electrónico y la energía total en el punto crítico de enlace 111 confirman una interacción molecular efectiva entre el timol y el residuo CYS145 del tipo de enlace de hidrógeno. Como se puede observar en la Figura 6, la densidad electrónica del oxígeno asociada a los pares libres de electrones se orienta hacia el átomo de hidrógeno. Además, se observa la formación de un punto crítico de enlace entre el oxígeno del timol y el grupo -CH de la CYS145 (ver representación en la Figura 6). La magnitud de la energía total sugiere una interacción débil, del tipo Van der Waals. En este caso, el átomo de hidrógeno del grupo -CH actúa como un aceptor de electrones provenientes del átomo de oxígeno. Sin embargo, la naturaleza débil de este enlace se debe a la falta de un enlace covalente directo entre el hidrógeno y el oxígeno [42], [43]
Adicionalmente, se observa la formación de un punto crítico de enlace entre el átomo de azufre del residuo CYS145 y un carbono del anillo aromático del timol (ver representación en la Figura 7). Los valores del Laplaciano electrónico y de la energía total sugieren una interacción covalente o parcialmente covalente de capa cerrada, debido a que los anillos aromáticos tienen densidades electrónicas altas como consecuencia de la deslocalización electrónica. Este hecho puede afectar la naturaleza y la fuerza de la interacción entre el átomo de azufre S19 y el átomo de carbono C24. Por otra parte, se observa la formación de un punto crítico de enlace entre un grupo -CH del timol y otro -CH del residuo CYS145 (ver Figura 8); este tipo de interacción es conocida como enlace de dihidrógeno y es muy similar al enlace de hidrógeno convencional, pero existe una fuerte interacción electrostática entre los átomos de hidrógenos con cargas parciales opuestas. Además, existe una contribución de fuerza de atracción de Van der Waals [44], [45]. El valor del Laplaciano electrónico obtenido de 0,0142 u.a. se encuentra dentro del rango reportado para este tipo de interacciones intermoleculares (0,0140,139 u.a.) [46].
Función de localización de electrones
La función de localización de electrones (ELF) es una herramienta útil para determinar la probabilidad de localizar pares de electrones. Además, ayuda a comprender el comportamiento de los electrones en un sistema molecular. En esta parte, el método ELF fue utilizado para el estudio de las interacciones entre el timol y los residuos HIS41 yCYS145 de la proteasa 3CL del virus SARS-CoV-2.
En la Figura 9 se muestra el mapa ELF para la interacción molecular entre el timol y el residuo HIS41 en el plano para los átomos H54, H12 y N14. Este mapa revela zonas de color rojo asociadas a un carácter covalente, lo cual indica una alta localización de electrones, característica de un enlace covalente, lo que se correlaciona con lo obtenido por el análisis QTAIM.
Por su parte, en la Figura 10 se presentan los mapas ELF para la interacción entre el timol (átomos C24 y O27) y el residuo CYS145 (átomos S19 y H14). En ambos casos se puede observar que los pares libres de electrones participan en la interacción molecular, tanto para el átomo de azufre (S19) como para el átomo de oxígeno (O27). Los pares de electrones libres del oxígeno están orientados hacia el átomo de hidrógeno del grupo amino de la CYS145. Este comportamiento es característico para la formación de un enlace de hidrógeno. Finalmente, para el átomo de carbono se observa que la distribución electrónica está influenciada por la naturaleza del anillo aromático.
En la Figura 11 se observan los mapas ELF para el timol (átomos O27 y H52) con el residuo CYS145 (átomo H43). En ambos casos se muestra que el átomo de hidrógeno del residuo CYS145 interviene en ambas interacciones moleculares. Adicionalmente, se observan zonas de alta concentración electrónica (color rojo) entre el hidrógeno H52 del timol y el H43 del residuo CYS145 característico de una interacción estérica no covalente hidrógeno-hidrógeno, lo cual corrobora lo observado para el Laplaciano electrónico en la Figura 8 con la presencia de un punto crítico de enlace entre el H52 del timol y el H43 del residuo CYS145 [41].
Método de gradiente reducido de la densidad (RDG)
En el análisis de interacción no covalente, uno de los métodos más utilizados es el gradiente reducido de la densidad (RDG por sus siglas en inglés); este método fue desarrollado por Johnson et al. [30], [31] y permite obtener la representación de las interacciones de enlaces de hidrógenos (H-B), Van der Waals (VdW) y efectos estéricos (E.E.). Adicionalmente, este método distingue las interacciones débiles como áreas de baja densidad electrónica y bajos gradientes de densidad reducida, por lo cual a partir de la isosuperficie RDG podemos identificar las diferentes regiones de interacción intermolecular simplemente examinando sus colores [47]. En el gráfico 2D de la Figura 12, la zona azul indica la presencia de interacciones no covalentes fuertes y atractivas como los enlaces de hidrógenos, la zona roja señala la presencia de interacciones repulsivas como efectos estéricos y la zona verde muestra la presencia de interacciones débiles como interacciones débiles de Van der Waals.
En la Figura 12 se observa una isosuperficie de color azul ubicada entre el átomo de oxígeno del timol y el grupo N-H del residuo CYS145, lo cual indica la presencia de un enlace de hidrógeno. Por su parte, los valores de signo (A2)p permiten determinar la naturaleza de las interacciones moleculares. Aquí se observó que el enlace de hidrógeno (H-B) está comprendido entre -0,030 u.a. < s¡gno(A2) p < -0,020 u.a., mientras que las interacciones moleculares de Van der Waals (isosuperficie de color verde), en un rango entre -0,010 u.a. < signo (A2)p < 0,010 u.a. Además, las isosuperficies de color rojo indican el efecto estérico que se encuentra en el anillo aromático del timol, lo que resulta en una fuerte repulsión. Estas interacciones repulsivas se observaron dentro de 0,020 u.a. < signo (A2)p < 0,030 u.a.
Modelo de gradiente independiente basado en la partición de Hirshfeld (IGMH)
El modelo de gradiente independiente basado en la partición de Hirshfeld (IGMH) permite calcular la densidad electrónica y el gradiente de la densidad electrónica de una molécula, basado en la partición de Hirshfeld, el cual consiste en dividir la densidad electrónica de la molécula en fragmentos atómicos [32], [47]. Adicionalmente, con este modelo se pueden visualizar y cuantificar las interacciones intermoleculares e intramoleculares. De allí, la importancia para el estudio de las interacciones entre el timol y los residuos HIS41 y CYS145. Uno de los parámetros por obtener a través de este modelo es el índice intrínseco de la fuerza de enlace para interacciones débiles (IBSIW) que permite cuantificar las fuerzas de las interacciones formadas [48].
Como se puede observar en la Figura 13, el análisis IGMH revela cuatro interacciones entre el timol y los residuos de la proteasa 3CL. Las zonas de color azul revelan una alta concentración de densidad electrónica característica de enlaces de hidrógeno; por su parte, las zonas de color rojo señalan una baja densidad electrónica. Los valores de IBSIW permiten ordenar en función de la fuerza de la interacción; el enlace de hidrógeno formado entre la CYS145 y el timol es el de mayor magnitud, lo cual indica una fuerte interacción.
Conclusiones
Se emplearon métodos computacionales para analizar las interacciones moleculares entre el timol y los residuos clave en el sitio activo de la proteasa 3CL del SARS-CoV-2. Los resultados obtenidos a través del acoplamiento molecular revelaron que el timol se une al sitio activo de la 3CLpro, con una energía de interacción de -3,784 kcal/mol y una constante de inhibición de 1,29 mM. El mapeo de las interacciones mostró la formación de enlaces de hidrógeno e interacciones hidrofóbicas con los residuos de aminoácidos HIS41, LEU141, GLY143, SER144, CYS145 e HIS163.
Asimismo, se caracterizaron las interacciones moleculares entre el timol y los residuos HIS41 y CYS145 utilizando las teorías QTAIM, ELF, RDG e IGMH. Estas teorías permitieron identificar la presencia de enlaces de hidrógeno y otras interacciones no covalentes. El análisis QTAIM mostró puntos críticos de enlace y la existencia de un enlace de hidrógeno entre el grupo OH del timol y el residuo CYS145 con valores de Laplaciano y densidad de energía total de 0,01036 y 0,0027 u.a. Además, ELF reveló cómo los pares de electrones libres en el oxígeno se orientan hacia el hidrógeno en CYS145, lo que favorece la formación del enlace de hidrógeno.
Finalmente, los análisis de RDG e IGMH para las interacciones no covalentes confirmaron la presencia de la formación del enlace de hidrógeno entre el timol y la CYS145 y la presencia de interacciones del tipo Van der Waals con la HIS41. Adicionalmente, en la formación de estas diversas interacciones el timol puede unirse de forma efectiva con los residuos CYS145 y HIS41 localizados en la región activa de la proteasa 3CLpro.