Computer-aided design of selective COX-2 inhibitors: comparative molecular field analysis, comparative molecular similarity indices analysis, and docking studies of some 1,2-diarylimidazole derivatives.

Comparative molecular field analysis and comparative molecular similarity indices analysis were performed on 114 analogues of 1,2-diarylimidazole to optimize their cyclooxygenase-2 (COX-2) selective antiinflammatory activities. These studies produced models with high correlation coefficients and good predictive abilities. Docking studies were also carried out wherein these analogues were docked into the active sites of both COX-1 and COX-2 to analyze the receptor ligand interactions that confer selectivity for COX-2. The most active molecule in the series (53) adopts an orientation similar to that of SC-558 (4-[5-(4-bromophenyl)-3-trifluoromethyl-1H-1-pyrozolyl]-1-benzenesulfonamide) inside the COX-2 active site while the least active molecule (101) optimizes in a different orientation. In the active site, there are some strong hydrogen-bonding interactions observed between residues His90, Arg513, and Phe518 and the ligands. Additionally, a correlation of the quantitative structure-activity relationship data and the docking results is found to validate each other and suggests the importance of the binding step in overall drug action.

NMR study of cyclic peptides with renin inhibitor activity.

Several cyclic analogues of renin inhibitors, based on Glu-D-Phe-Lys motif have been investigated by NMR spectroscopy and molecular dynamics calculations (MD). The 15 membered macrocycle, resulting from Glu and Lys side-chain cyclization, exhibits conformational preference. The structural evidence from NMR shows the presence of hydrogen bond between Lys NH and Glu side-chain carbonyl, resulting in a 10 membered pseudo beta-turn-like structure. The structure of the cyclic moiety is similar in all the peptides, which takes at least two conformations around Calpha-Cbeta in Glu side chain. The restrained MD calculations further support such observations and show that the macrocycle is fairly rigid, with two conformations about the Glu Calpha-Cbeta bond. The linear peptide appendages, which are essential for activity in cyclic peptides, show an extended structure in the beta-region of Ramchandran plot. These calculations also demonstrate that for the most active peptide, two major conformers each exist about the Calpha-CO bond of the Lys, D-Trp and Leu residues. In this peptide, the cyclic moiety presents a negatively charged surface formed due to the carbonyl oxygens, which are thus available to form hydrogen bonds with the receptor. The linear fragment presents further binding sites with a surface which has the hydrophobic side chains of D-Trp, Leu and D-Met on one side and carbonyls on the other side.

Anthrax lethal factor inhibitors as potential countermeasure of the infection.

Anthrax Lethal Factor (LF) is a zinc-dependent metalloprotease, one of the virulence factor of anthrax infection. Three forms of the anthrax infection have been identified: cutaneous (through skin), gastrointestinal (through alimentary tract), and pulmonary (by inhalation of spores). Anthrax toxin is composed of protective antigen (PA), lethal factor (LF), and edema factor (EF). Protective antigen mediates the entry of Lethal Factor/Edema Factor into the cytosol of host cells. Lethal factor (LF) inactivates mitogen-activated protein kinase kinase inducing cell death, and EF is an adenylyl cyclase impairing host defenses. In the past few years, extensive studies are undertaken to design inhibitors targeting LF. The current review focuses on the small molecule inhibitors targeting LF activity and its structure activity relationships (SAR).

Design, synthesis and screening studies of potent thiazol-2-amine derivatives as fibroblast growth factor receptor 1 inhibitors.

Fibroblast growth factor receptor 1 (FGFR1) a tyrosine kinase receptor, plays important roles in angiogenesis, embryonic development, cell proliferation, cell differentiation, and wound healing. The FGFR isoforms and their receptors (FGFRs) considered as a potential targets and under intense research to design potential anticancer agents. Fibroblast growth factors (FGF's) and its growth factor receptors (FGFR) plays vital role in one of the critical pathway in monitoring angiogenesis. In the current study, quantitative pharmacophore models were generated and validated using known FGFR1 inhibitors. The pharmacophore models were generated using a set of 28 compounds (training). The top pharmacophore model was selected and validated using a set of 126 compounds (test set) and also using external validation. The validated pharmacophore was considered as a virtual screening query to screen a database of 400,000 virtual molecules and pharmacophore model retrieved 2800 hits. The retrieved hits were subsequently filtered based on the fit value. The selected hits were subjected for docking studies to observe the binding modes of the retrieved hits and also to reduce the false positives. One of the potential hits (thiazole-2-amine derivative) was selected based the pharmacophore fit value, dock score, and synthetic feasibility. A few analogues of the thiazole-2-amine derivative were synthesized. These compounds were screened for FGFR1 activity and anti-proliferative studies. The top active compound showed 56.87% inhibition of FGFR1 activity at 50 µM and also showed good cellular activity. Further optimization of thiazole-2-amine derivatives is in progress.

Ligand coordinate analysis of SC-558 from the active site to the surface of COX-2: a molecular dynamics study.

We have performed a ligand coordinate analysis to monitor the movement of the inhibitor SC-558 from the active site of the COX-2 protein to the exterior using molecular dynamics techniques. This study provides an insight into the intermolecular interactions formed by the ligand during this journey. The published crystal structure of COX-2 with SC-558 in the active site (1cx2) was taken, and the ligand was moved incrementally in 13 steps. At each of these points on the path, exhaustive minimization and dynamics calculations were performed. The role of water was found to be important in these computations. An average structure was obtained from 250 conformations at each point and minimized. At each point on the path, the 10 lowest-energy conformations were also selected; a consideration of the average and lowest conformations provides fine details on the consistency of specific and strong interactions, and also on the geometry of the ligand. The movement of the ligand through the protein may be divided into three stages that are distinguished from each other because of energy and geometry discontinuities in both the ligand and the protein. The first of these covers the region between the active site and the point at 5.8 A from it. The second, which covers the distance between 8.2 and 10.0 A and is associated with maximum energetic and structural instability, is of critical importance. The third stage covers the distance between 10.5 A and the exterior and represents a stage of increasing hydration and expulsion of the ligand from the protein. Our results provide a confirmation for the existence of a shallow cavity near the protein surface in which the ligand is bound reversibly. By examining the residues that show maximum mobility, one obtains an idea of the gating mechanism that governs the entry and exit of the protein into or from the deep pocket that contains the active site. We note, however, that the variation of the root-mean-square deviation of all residues begins to increase almost as soon as the ligand leaves the active site, and even before there are any changes in the gate inter-residue distances. This loosening of the protein even before the gate opens might be a part of the enthalpy-entropy balance that accompanies the ligand's passage through the protein. Our results provide an energy profile of the ligand during its entry/exit into/from the protein and can, in principle, enable one to assess the residence time, which in turn may be associated or indirectly correlated with adverse cardiovascular side effects of nonsteroidal anti-inflammatory drugs. We believe that similar analyses for other selected COX-2-specific inhibitors will provide a measure (or prediction) of possible toxicity effects.

Virtual screening of 4-anilinoquinazoline analogues as EGFR kinase inhibitors: importance of hydrogen bonds in the evaluation of poses and scoring functions.

Virtual Screening (VS) is a computational technique that allows selection and ranking of possible hits from a library of compounds. We have carried out VS on 128 selected EGFR kinase inhibitors with GOLD and LigandFit. From the experimental crystal structure of the erlotinib-EGFR complex, three key hydrogen bonds were identified as responsible for anchoring the ligand in the active site. These are of the N-H...N, O(w)-H...N, and C-H...O types. Failure to include the hydrogen-bonded water molecule that forms the O(w)-H...N bond leads to incorrect results. Of the three interactions, the C-H...O formed by an activated C-H group is the best conserved. On the basis of the efficacy of these hydrogen bonds, the poses were classified into one of three categories: close, shifted, and misoriented. In the VS context, all three interactions need to be modeled correctly so that correct poses and affinities are obtained, and this happens in ligands of the close variety. Cross scoring wherein the poses from one software are input into another for scoring and consensus scoring wherein the scores from various software packages are weighted are also helpful in obtaining better agreements.

Analogue based design of MMP-13 (Collagenase-3) inhibitors.

3D-QSAR studies using MFA and RSA methods were performed on a series of 39MMP-13 inhibitors. Model developed by MFA method has a r(2)(cv) (cross-validated) of 0.616 while its r(2) (conventional) value is 0.822. For the RSA model r(2)(cv) and r(2) are 0.681 and 0.847, respectively. Both the models indicate good internal as well as external predictive abilities. These models provide crucial information about the field descriptors for the design of potential inhibitors of MMP-13.

Design of EGFR kinase inhibitors: a ligand-based approach and its confirmation with structure-based studies.

Three-dimensional quantitative structure-activity relationship (3D-QSAR) models were developed for 100 anilinoquinazolines, inhibiting epidermal growth factor receptor (EGFR) kinase. The studies included molecular field analysis (MFA) and receptor surface analysis (RSA). The cross-validated r2 (r2cv) values are 0.81 and 0.79 for MFA and RSA, respectively. The predictive ability of these models was validated by 28 test set molecules. The results of the best QSAR model were further compared with structure-based investigations using docking studies with the crystal structure of EGFR kinase domain. The results helped to understand the nature of substituents at the 6- and 7-positions, thereby providing new guidelines for the design of novel inhibitors.

3-D QSAR studies on new dibenzyltin(IV) anticancer agents by comparative molecular field analysis (CoMFA).

Dibenzyltin(IV)dichloride and dibenzyltin(IV)diisothiocyanate derivatives with N,S-donor ligands show significant cytotoxic activities against human cancer cell lines, and are well compared to analogous dialkyltin(IV) derivatives. CoMFA models were generated for the first time for these organotin derivatives using the cytotoxic activities (against two human tumor cell lines, MCF-7, a mammary carcinoma and WiDr, a colon carcinoma) of 21 complexes. High r(2) and r(2)(cv) values for both CoMFA models indicate good predictive power for the models.