Aplicação de modelagem molecular e de formalismo do CAMD (Computer-Aided Molecular Design) na elucidação do mecanismo de ação de inibidores de metalopropteinases de matriz

Publication year: 2015
Theses and dissertations in Portugués presented to the Universidade de São Paulo to obtain the academic title of Doutor. Leader: Barros, Silvia Berlanga de Moraes

As metaloproteinases de matriz (MMP) são enzimas superexpressas em quase todos os tumores humanos, sendo que os subtipos MMP-2 e MMP-9 têm sido associados ao potencial metastático e prognóstico desfavorável em neoplasias malignas como, por exemplo, melanoma metastático e glioma. Compostos capazes de inibir a atividade destas enzimas podem representar potenciais agentes terapêuticos. O composto 4-nerolidilcatecol (4-NC), isolado de plantas do gênero Pothomorphe, apresentou resultados promissores para o tratamento do melanoma e glioma e foi capaz de atuar em várias etapas bioquímicas importantes envolvidas na progressão dessas patologias, inclusive inibindo MMP-2 e MMP-9. No entanto, o mecanismo de ação do 4-NC não está completamente elucidado. O presente estudo envolveu a aplicação de métodos de modelagem molecular e de formalismos do planejamento de novas moléculas auxiliado por computador, CAMD (Computer-Aided Molecular Design) a fim de explorar a interação entre esta molécula e as enzimas MMP-2 e MMP-9, além de planejar novos inibidores para estes alvos. Análise exploratória de dados, que compreende a análise de agrupamentos hierárquicos e de componentes principais. foi desenvolvida para um conjunto de hidroxamatos (N=64) descritos como inibidores de MMP-2 e MMP-9, a fim de identificar as propriedades moleculares que mais influenciavam o processo de discriminação dos compostos. As propriedades termodinâmicas, eletrônicas e estéricas foram importantes para descrever os compostos mais ativos no conjunto de dados da MMP-2. Para a MMP-9, o coeficiente de distribuição (ClogD) em pH 1,5 foi relevante no processo de discriminação do conjunto. A presença de substituintes volumosos na porção R3 parece ser crucial para o conjunto de inibidores investigados. Esta região está envolvida em interações moleculares com a cavidade S1 de ambas as enzimas, mas há um limite de volume a ser considerado para estes substituintes. O formalismo QSAR-4D independente do receptor (IR) foi aplicado ao mesmo conjunto de dados e permitiu estabelecer o mapeamento do farmacóforo, além de explorar diferentes alinhamentos para a obtenção da hipótese de conformação bioativa prevista pelo melhor modelo de QSAR. OS modelos QSAR apresentaram boa capacidade de previsão, auxiliaram na proposição de novos inibidores e estimaram a atividade do 4-NC. Com o melhor modelo QSAR para MMP-9 (N=64), a atividade prevista para o 4-NC foi classificada na faixa dos inibidores com atividade moderada. Entretanto, o melhor modelo QSAR obtido para MMP-2 (N=38) não foi capaz de prever, de forma adequada, a atividade de compostos com arcabouço químico diferente daqueles utilizados na construção dos modelos. Estudos de ancoramento molecular foram desenvolvidos para investigar a orientação do 4-NC no sitio catalítico das duas enzimas e as interações que poderiam ser estabelecidas nestes complexos. Duas conformações favoráveis foram encontradas. Simulações computacionais de dinâmica molecular foram desenvolvidas com os complexos mais promissores selecionados nos estudos de ancoramento, a fim de obter informações mais detalhadas e de maior confiabilidade. sobre suas interações intermoleculares. O 4-NC tende a se orientar no sítio de forma a acomodar sua cadeia lateral no bolso S1 adjacente ao sítio catalítico em ambas as enzimas. Ensaios de zimografia também foram realizados com o objetivo de elucidar possíveis contribuições da cadeia lateral e do núcleo catecólico do 4-NC na atividade inibitória frente às enzimas em estudo. O núcleo catecólico parece ser o responsável por sua atividade, pois o composto 1,2dimetoxibenzeno, que possui as hidroxilas bloqueadas por grupos metil, não foi capaz de exercer atividade inibitória significante frente à MMP-2 e MMP-9. Estudos de voltametria reforçaram a hipótese de que o 4-NC tem a capacidade de quelar os íons zinco presentes no tampão de incubação
Matrix metalloproteinases (MMP) enzymes are overexpressed in almost all human tumors, and MMP-2 and MMP-9 subtypes have been associated with metastatic potential and poor prognosis in malignant tumors, such as metastatic melanoma and glioma. Compounds capable of inhibiting the activity of theses enzymes would be considered as potential therapeutic agents. The 4-nerolidylcatechol compound (4-NC), isolated from plants of genus Pothomorphe, has showed promising results in the treatment of melanoma and glioma, and was able to act in several important biochemical steps involved in the progression of these diseases, as well as inhibiting MMP-2 and MMP-9. However, the 4-NC mechanism of action is not completely understood. This study has involved the application of molecular modeling methods and formalisms of computer-aided molecular design (CAMD) in order to explore the interaction between 4-NC and MMP-2/MMP-9, and to design new inhibitors for these targets. Exploratory data analysis, which comprises hierarchical cluster analysis and principal components analysis, was performed to a set of hydroxamates (N=64). previously reported as MMP-2 and MMP-9 inhibitors, in order lo identify the molecular properties that is most critical for the discrimination process regarding the investigated compounds. The thermodynamic, electronic, and steric properties were: quite important to describe the highly active compounds in the data set of MMP-2, whereas the apparent partition coefficient (ClogD) at pH 1.5 was the property more relevant for MMP-9 data set. The presence of bulky substituents on the R3 moiety seems to be crucial for this set of inhibitors due to the molecular interaction with the S1 subsite of both enzymes. However, there is a limit regarding the substituents volume in this region. Receptor independent (RI) 4D-QSAR analysis was applied lo the same data set and it was possible to establish the pharmacophore mapping, besides to explore different alignments in order to generate the hypothesized bioactive conformation through the best QSAR model. The QSAR models have presented good predictability, assisted in proposing new inhibitors, and estimated the activity of 4-NC. Regarding the best QSAR model for MMP-9 (N=64), the 4-NC predicted activity was classified in the range of the moderate active inhibitors. The best QSAR model obtained for MMP-2 (N=38), however was not able to properly predict the activity for compounds with different chemical scaffold from those used to build up the QSAR model. Molecular docking studies have been developed to investigate the 4-NC binding mode into the catalytic site of the two enzymes and the interactions that could be established in those complexes. The results have shown two favorable conformers regarding the MMP inhibition. Molecular dynamics computational simulation were combined to molecular docking studies in order to obtain more detailed and reliable information regarding the intermolecular interactions of each complex. The 4-NC molecule tends to accommodate the side chain in the S1 pocket adjacent to the catalytic site in both enzymes. Experimental zymography assays were also performed to elucidate the possible contribution of the side chain and the catechol core in the 4-NC inhibitory activity against the MMP-2 and MMP-9 enzymes. The catechol core seems to be responsible for its activity, since the 1,2 dimethoxybenzene compound, which has the hydroxyl blocked by a methyl group, was not able to exert any significant inhibition on enzymes. Voltametric assays confirmed the hypothesis that 4-NC chelates zinc ions present in the incubation buffer

More related