Abstract
This work proposes a computational procedure for structure-based lead generation and optimization, which relies on the complementarity of the protein-ligand interactions. This procedure takes as input the known structure of a protein-ligand complex. Retaining the positions of the ligand heavy atoms in the protein binding site it designs structurally similar compounds considering all possible combinations of atomic species (N, C, O, CH(3), NH, etc). Compounds are ranked based on a score which incorporates energetic contributions evaluated using molecular mechanics force fields. This procedure was used to design new inhibitor molecules for three serine/threonine protein kinases (p38 MAP kinase, p42 MAP kinase (ERK2), and c-Jun N-terminal kinase 3 (JNK3)). For each enzyme, the calculations produce a set of potential inhibitors whose scores are in agreement with IC50 data and Ki values. Furthermore, the native ligands for each protein target, scored within the five top-ranking compounds predicted by our method, one of the top-ranking compounds predicted to inhibit JNK3 was synthesized and his inhibitory activity confirmed against ATP hydrolysis. Our computational procedure is therefore deemed to be a useful tool for generating chemically diverse molecules active against known target proteins.