Estudio estructural del rol inhibitorio de derivados de naftaleno frente a la enzima Plpro de SARS-CoV

Abstract

La proteasa SARS-CoV PLpro está asociada a diferentes funciones que ocasionan el antagonismo de la respuesta inmune innata antiviral del huésped, como lo es la deubiquitinación y la desISGilación de las proteínas de las células de este. Dada la función que desempeña, la atención en la PLpro como diana terapéutica para el diseño de fármacos contra el SARS-CoV, MERS-CoV y actualmente el SARS-CoV-2, ha ido en aumento y se han reportado compuestos candidatos para detener la actividad de esta enzima. Entre estos, se reportó un conjunto de inhibidores no covalentes basados en naftaleno, que incluye compuestos que muestran una alta capacidad inhibitoria contra la PLpro de SARS-CoV. La comprensión de las características estructurales que facilitan la interacción entre la PLpro de coronavirus y sus ligandos es de gran utilidad para entender el funcionamiento de esta enzima y su capacidad para ser intervenida con inhibidores más efectivos. En este trabajo se revelaron características estructurales de los derivados del naftilo que explican su actividad contra la enzima PLpro de SARS-CoV, a través de un estudio detallado de las interacciones ligando-receptor. El protocolo empleado se inició con el método docking molecular, el cual permitió considerar las orientaciones de los ligandos en el sitio de unión de la proteasa. Posterior a ello, se realizaron diversos análisis quimiométricos: se utilizó el servidor LigRMSD para comparar las orientaciones de los inhibidores estudiados, se realizó un análisis de huellas de interacción para describir las interacciones recurrentes entre los residuos de la proteína y los grupo del ligando en los complejos y finalmente se construyeron modelos 3D-QSAR para explicar la relación entre la estructura de los compuestos y su actividad inhibitoria. Se describieron las interacciones entre diferentes porciones del ligando y varios residuos claves (Asp165, Tyr269, Gln270 y Tyr265), además de la ocupación de estos inhibidores en los subsitios S3 y S4 de la proteasa. Se encontró que la serie de compuestos derivados del naftaleno adoptan orientaciones similares en el sitio de unión de PLpro, con una alta semejanza con las orientaciones de los ligandos co-cristalizados. Sumado a esto, los modelos 3D-QSAR construidos en base a diferentes estrategias de alineamiento, no cumplieron en su totalidad con los requisitos estadísticos para ser considerados modelos con alto poder predictivo. En base a estos resultados, se propone el uso de una nueva metodología para explicar la actividad biológica de la serie de inhibidores. // ABSTRACT: The SARS-CoV protease PLpro is associated with different functions that cause antagonism of the host innate antiviral immune response, such as deubiquitination and deISGylation of host cell proteins. Given the role it plays, attention on PLpro as a therapeutic target for the design of drugs against SARS-CoV, MERS-CoV and currently SARS-CoV-2, has been increasing and candidate compounds have been reported to reduce the activity of this enzyme. Studying the structural features that facilitate the interaction between coronavirus PLpro and its ligands is of great utility in understanding the functioning of this enzyme and its ability to be targeted with more effective inhibitors. In this work, structural features of naphthyl derivatives against SARS-CoV PLpro enzyme were revealed through a detailed study of ligand-receptor interactions. The protocol started with the molecular docking method, which allowed to consider the ligand orientations at the protease binding site. Subsequently, several chemometric analyses were performed: the LigRMSD server was used to compare the orientations of the studied inhibitors, an interaction fingerprint analysis was performed to describe the recurrent interactions between the protein residues and the ligand groups in the complexes; finally, 3D-QSAR models were constructed to explain the relationship between the structure of the compounds and their inhibitory activity. The interactions between different portions of the ligand and key residues (Asp165, Tyr269, Gln270 and Tyr265) were described, in addition to the occupancy of these inhibitors at the S3 and S4 sub-sites of the protease. It was found that the serie of naphthalene-derived compounds adopted similar orientations at the PLpro binding site, with high similarity to the orientations of co-crystallized ligands. In addition, the 3D-QSAR models constructed based on different alignment strategies did not fully meet the statistical requirements to be considered models with high predictive power. Based on these results, a new methodology is proposed to explain the biological activity of the inhibitor series.