Validation of Molecular Docking Calculations Involving FGF-1 and FGF-2

Synthesis and biological evaluation of direct thrombin inhibitors bearing 4-(piperidin-1-yl)pyridine at the P1 position with potent anticoagulant activity

The design and synthesis of a new class of nonpeptide direct thrombin inhibitors, built on the structure of 1-(pyridin-4-yl)piperidine-4-carboxamide, are described. Starting from a strongly basic 1-amidinopiperidine derivative (6) showing poor thrombin (fIIa) and factor Xa (fXa) inhibition activities, anti-fIIa activity and artificial membrane permeability were considerably improved by optimizing the basic P1 and the X-substituted phenyl P4 binding moieties. Structure-activity relationship studies, usefully complemented with molecular modeling results, led us to identify compound 13b, which showed excellent fIIa inhibition (Ki = 6 nM), weak anti-Xa activity (Ki = 5.64 μM), and remarkable selectivity over other serine proteases (e.g., trypsin). Compound 13b showed in vitro anticoagulant activity in the low micromolar range and significant membrane permeability. In mice (ex vivo), 13b demonstrated anticoagulant effects at 2 h after oral dosing (100 mg·kg(-1)), with a significant 43% prolongation of the activated partial thromboplastin time (aPTT), over controls (P < 0.05).

A leishmaniasis study: structure-based screening and molecular dynamics mechanistic analysis for discovering potent inhibitors of spermidine synthase

Protozoa Leishmania donovani (Ld) is the main cause of the endemic disease leishmaniasis. Spermidine synthase (SS), an important enzyme in the synthetic pathway of polyamines in Ld, is an essential element for the survival of this protozoan. Targeting SS may provide an important aid for the development of drugs against Ld. However, absence of tertiary structure of spermidine synthase of Leishmania donovani (LSS) limits the possibilities of structure based drug designing. Presence of the same enzyme in the host itself further challenges the drug development process. We modeled the tertiary structure of LSS using homology modeling approach making use of homologous X-ray crystallographic structure of spermidine synthase of Trypanosoma cruzi (TSS) (2.5Å resolution). The modeled structure was stabilized using Molecular Dynamics simulations. Based on active site structural differences between LSS and human spermidine synthase (HSS), we screened a large dataset of compounds against modeled protein using Glide virtual screen docking and selected two best inhibitors based on their docking scores (-10.04 and -13.11 respectively) with LSS and having least/no binding with the human enzyme. Finally Molecular Dynamics simulations were used to assess the dynamic stability of the ligand bound structures and to elaborate on the binding modes. This article is part of a Special Issue entitled: Computational Methods for Protein Interaction and Structural Prediction.

Substitution of aminomethyl at the meta-position enhances the inactivation of O6-alkylguanine-DNA alkyltransferase by O6-benzylguanine

O(6)-Benzylguanine is an irreversible inactivator of O(6)-alkylguanine-DNA alkyltransferase currently in clinical trials to overcome alkyltransferase-mediated resistance to certain cancer chemotherapeutic alkylating agents. In order to produce more soluble alkyltransferase inhibitors, we have synthesized three aminomethyl-substituted O(6)-benzylguanines and the three methyl analogs and found that the substitution of aminomethyl at the meta-position greatly enhances inactivation of alkyltransferase, whereas para-substitution has little effect and ortho-substitution virtually eliminates activity. Molecular modeling of their interactions with alkyltransferase provided a molecular explanation for these results. The square of the correlation coefficient (R(2)) obtained between E-model scores (obtained from GLIDE XP/QPLD docking calculations) vs log(ED(50)) values via a linear regression analysis was 0.96. The models indicate that the ortho-substitution causes a steric clash interfering with binding, whereas the meta-aminomethyl substitution allows an interaction of the amino group to generate an additional hydrogen bond with the protein.

Identification and characterization of coumestans as novel HCV NS5B polymerase inhibitors

The hepatitis C virus (HCV) NS5B is essential for viral RNA replication and is therefore a prime target for development of HCV replication inhibitors. Here, we report the identification of a new class of HCV NS5B inhibitors belonging to the coumestan family of phytoestrogens. Based on the in vitro NS5B RNA-dependent RNA polymerase (RdRp) inhibition in the low micromolar range by wedelolactone, a naturally occurring coumestan, we evaluated the anti-NS5B activity of four synthetic coumestan analogues bearing different patterns of substitutions in their A and D rings, and observed a good structure-activity correlation. Kinetic characterization of coumestans revealed a noncompetitive mode of inhibition with respect to nucleoside triphosphate (rNTP) substrate and a mixed mode of inhibition towards the nucleic acid template, with a major competitive component. The modified order of addition experiments with coumestans and nucleic acid substrates affected the potencies of the coumestan inhibitors. Coumestan interference at the step of NS5B-RNA binary complex formation was confirmed by cross-linking experiments. Molecular docking of coumestans within the allosteric site of NS5B yielded significant correlation between their calculated binding energies and IC(50) values. Coumestans thus add to the diversifying pool of anti-NS5B agents and provide a novel scaffold for structural refinement and development of potent NS5B inhibitors.

Can we use docking and scoring for hit-to-lead optimization?

Docking and scoring is currently one of the tools used for hit finding and hit-to-lead optimization when structural information about the target is known. Docking scores have been found useful for optimizing ligand binding to reproduce experimentally observed binding modes. The question is, can docking and scoring be used reliably for hit-to-lead optimization? To illustrate the challenges of scoring for hit-to-lead optimization, the relationship of docking scores with experimentally determined IC(50) values measured in-house were tested. The influences of the particular target, crystal structure, and the precision of the scoring function on the ability to differentiate between actives and inactives were analyzed by calculating the area under the curve of receiver operator characteristic curves for docking scores. It was found that for the test sets considered, MW and sometimes ClogP were as useful as GlideScores and no significant difference was observed between SP and XP scores for differentiating between actives and inactives. Interpretation by an expert is still required to successfully utilize docking and scoring in hit-to-lead optimization.

3D QSAR and molecular docking studies of benzimidazole derivatives as hepatitis C virus NS5B polymerase inhibitors

The urgent need for novel HCV antiviral agents has provided an impetus for understanding the structural requisites of NS5B polymerase inhibitors at the molecular level. Toward this objective, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) of 67 HCV NS5B polymerase inhibitors were performed using two methods. First, ligand-based 3D QSAR studies were performed based on the lowest energy conformations employing the atom fit alignment method. Second, receptor-based 3D QSAR models were derived from the predicted binding conformations obtained by docking all NS5B inhibitors at the allosteric binding site of NS5B (PDB ID: 2dxs). Results generated from the ligand-based model were found superior (r2cv values of 0.630 for CoMFA and 0.668 for CoMSIA) to those obtained by the receptor-based model (r2cv values of 0.536 and 0.561 for CoMFA and CoMSIA, respectively). The predictive ability of the models was validated using a structurally diversified test set of 22 compounds that had not been included in a preliminary training set of 45 compounds. The predictive r2 values for the ligand-based CoMFA and CoMSIA models were 0.734 and 0.800, respectively, while the corresponding predictive r2 values for the receptor-based CoMFA and CoMSIA models were 0.538 and 0.639, respectively. The greater potency of the tryptophan derivatives over that of the tyrosine derivatives was interpreted based on CoMFA steric and electrostatic contour maps. The CoMSIA results revealed that for a NS5B inhibitor to have appreciable inhibitory activity it requires hydrogen bond donor and acceptor groups at the 5-position of the indole ring and an R substituent at the chiral carbon, respectively. Interpretation of the CoMFA and CoMSIA contour maps in context of the topology of the allosteric binding site of NS5B provided insight into NS5B-inhibitor interactions. Taken together, the present 3D QSAR models were found to accurately predict the HCV NS5B polymerase inhibitory activity of structurally diverse test set compounds and to yield reliable clues for further optimization of the benzimidazole derivatives in the data set.

An Experimental and Molecular-Modeling Study of the Binding of Linked Sulfated Tetracyclitols to FGF-1 and FGF-2

The experimental binding affinities of a series of linked sulfated tetracyclitols [Cyc2N-R-NCyc2, where Cyc = C6H6(OSO3Na)3 and R = (CH2)n (n = 2-10), p-xylyl or (C2H4)2-Ncyc] for the fibroblast growth factors FGF-1 and FGF-2 have been measured by using a surface plasmon resonance assay. The KD values range from 7.0 nM to 1.1 microM for the alkyl-linked ligands. The binding affinity is independent of the flexibility of the linker, as replacement of the alkyl linker with a rigid p-xylyl group did not affect the KD. Calculations suggest that binding modes for the p-xylyl-linked ligand are similar to those calculated for the flexible alkyl-linked tetracyclitols. The possible formation of cross-linked FGF:cyclitol complexes was examined by determining KD values at increasing protein concentrations. No changes in KD were observed; this suggesting that only 1:1 complexes are formed under these assay conditions. Monte Carlo multiple-minima calculations of low-energy conformers of the FGF-bound ligands showed that all of the sulfated tetracyclitol ligands can bind effectively in the heparan sulfate-binding sites of FGF-1 and FGF-2. Binding affinities of these complexes were estimated by the Linear Interaction Energy (LIE) method to within a root-mean-square deviation of 1 kcal mol(-1) of the observed values. The effect of incorporating cations to balance the overall charge of the complexes during the LIE calculations was also explored.