RESUMO
The p38-mitogen-activated protein kinases (p38-MAPKs) belong to a family of serine-threonine kinases activated by pro-inflammatory or stressful stimuli that are known to be involved in several diseases. Their biological importance, related to the release of inflammatory pro-cytokines such as tumor necrosis factor-alpha (TNF-alpha) and interleukin-1 (IL-1), has generated many studies aiming at the development of selective inhibitors for the treatment of inflammatory diseases. In this work, we developed receptor-based three dimensional (3D) quantitative structure-activity relationship (QSAR) models for a series of 33 pyridinyl imidazole compounds [Liverton et al. (1999) 42:2180], using a methodology named free-energy force-field (FEFF) [Tokarski and Hopfinger (1997) 37:792], in which scaled intra- and intermolecular energy terms of the Assisted Model Building Energy Refinement (AMBER) force field combined with a hydration-shell solvation model are the independent variables used in the QSAR studies. Multiple temperature molecular-dynamics simulations (MDS) of ligand-protein complexes and genetic-function approximation (GFA) were employed using partial least squares (PLS) as the fitting functions to develop FEFF-3D-QSAR models for the binding process. The best model obtained in the FEFF-3D-QSAR receptor-dependent (RD) method shows the importance of the van der Waals energy change upon binding and the electrostatic energy in the interaction of ligands with the receptor. The QSAR equations described here show good predictability and may be regarded as representatives of the binding process of ligands to p38-MAPK. Additionally, we have compared the top FEFF-3D-QSAR model with receptor independent (RI) 4D-QSAR models developed in a recent study [Romeiro et al. (2005) 19:385].
Assuntos
Imidazóis/química , Inibidores de Proteínas Quinases/química , Relação Quantitativa Estrutura-Atividade , Proteínas Quinases p38 Ativadas por Mitógeno/antagonistas & inibidores , Interações Hidrofóbicas e Hidrofílicas , Imidazóis/farmacologia , Ligantes , Ligação Proteica , Conformação Proteica , Eletricidade Estática , Temperatura , TermodinâmicaRESUMO
The p38-mitogen-activated protein kinase (p38-MAPK) plays a key role in lipopolysaccharide-induced tumor necrosis factor-alpha (TNF-alpha) and interleukin-1 (IL-1) release during the inflammatory process, emerging as an attractive target for new anti-inflammatory agents. Four-dimensional quantitative structure-activity relationship (4D-QSAR) analysis [Hopfinger et al., J. Am. Chem. Soc., 119 (1997) 10509] was applied to a series of 33 (a training set of 28 and a test set of 5) pyridinyl-imidazole and pyrimidinyl-imidazole inhibitors of p38-MAPK, with IC50 ranging from 0.11 to 2100 nM [Liverton et al., J. Med. Chem., 42 (1999) 2180]. Five thousand conformations of each analogue were sampled from a molecular dynamics simulation (MDS) during 50 ps at a constant temperature of 303 K. Each conformation was placed in a 2 angstroms grid cell lattice for each of three trial alignments. 4D-QSAR models were constructed by genetic algorithm (GA) optimization and partial least squares (PLS) fitting, and evaluated by leave-one-out cross-validation technique. In the best models, with three to six terms, the adjusted cross-validated squared correlation coefficients, Q2adj, ranged from 0.67 to 0.85. Model D (Q2adj = 0.84) was identified as the most robust model from alignment 1, and it is representative of the other best models. This model encompasses new molecular regions as containing pharmacophore sites, such as the amino-benzyl moiety of pyrimidine analogs and the N1-substituent in the imidazole ring. These regions of the ligands should be further explored to identify better anti-inflammatory inhibitors of p38-MAPK.