Identification of selective MMP-9 inhibitors through multiple e-pharmacophore,ligand-based pharmacophore,molecular docking,and density functional theory approaches

Matrix metalloproteinase-9 (MMP-9) is a significant target for the development of drugs for the treatment of arthritis, CNS disorders, and cancer metastasis. The structure-based and ligand-based methods were used for the virtual screening (VS) of database compounds to obtain potent and selective MMP-9 inhibitors. Experimentally known MMP-9 inhibitors were used to grow up ligand-based three pharmacophore models utilizing Schrodinger suite. The X-ray crystallographic structures of MMP-9 with different inhibitors were used to develop five energy-optimized structure-based (e-pharmacophore) models. All developed pharmacophores were validated and applied to screen the Zinc database. Pharmacophore matched compounds were subjected to molecular docking to retrieve hits with novel scaffolds. The molecules with diverse structures, high docking scores and low binding energies for various crystal structures of MMP-9, were selected as final hits. The Induced fit docking (IFD) analysis provided significant information about the driving of inhibitor to approve a suitable bioactive conformational position in the active site of protein. Since charge transfer reaction occurs during receptor–ligand interaction, therefore, electronic features of hits (ligands) are interesting parameters to explain the binding interactions. Density functional theory (DFT) at B3LYP/6-31G* level was utilized to explore electronic features of hits. The docking study of hits using AutoDock was helpful to establish the binding interactions. The study illustrates that the combined pharmacophore approach is advantageous to identify diverse hits which have better binding affinity to the active site of the enzyme for all possible bioactive conformations. The approach used in the study is worthy to design drugs for other targets.     

3D-QSAR and docking studies of aldehyde inhibitors of human cathepsin K

In order to better understand the structural and chemical features of human cathepsin K (CatK), which is an important cysteine protease in the pathogenesis of osteoporosis, the 3D-QSAR (CoMFA) studies were conducted on recently explored aldehyde compounds with known CatK inhibitory activities. The genetic algorithm of GOLD2.2 has been employed to position 59 aldehyde compounds into the active sites of CatK to determine the probable binding conformation. Good correlations between the predicted binding free energies and the experimental inhibitory activities suggested that the identified binding conformations of these potential inhibitors are reliable. The docking results also provided a reliable conformational alignment scheme for 3D-QSAR model. Based on the docking conformations, highly predictive comparative molecular field analysis (CoMFA) was performed with q2 value of 0.723. The predictive ability was validated by some compounds that were not included in the training set. Furthermore, the CoMFA model was mapped back to the binding sites of CatK, to get a better understanding of vital interactions between the aldehyde compounds and the protease. The CoMFA field distributions are in good agreement with the structural characteristics of the binding groove of the CatK, which suggested that the n-Bu in R4 position is the favor group substitute at P1 and moderate groups in R2 group are required on P2 substitute. In addition, 3D-QSAR results also demonstrated that aldehyde is an important pharmacophore because of electrostatic effect. These results, together with the good correlations between the inhibitory activities and the binding free energies predicted by GOLD2.2, demonstrated the power of combining docking/QSAR approach to explore the probable binding conformations of compounds at the active sites of the protein target, and further provided useful information in understanding the structural and chemical features of CatK in designing and finding new potential inhibitors.     

Docking and hydropathic scoring of polysubstituted pyrrole compounds with antitubulin activity

Compounds that bind at the colchicine site of tubulin have drawn considerable attention with studies indicating that these agents suppress microtubule dynamics and inhibit tubulin polymerization. Data for 18 polysubstituted pyrrole compounds are reported, including antiproliferative activity against human MDA-MB-435 cells and calculated free energies of binding following docking the compounds into models of alphabeta-tubulin. These docking calculations coupled with HINT interaction analyses are able to represent the complex structures and the binding modes of inhibitors such that calculated and measured free energies of binding correlate with an r(2) of 0.76. Structural analysis of the binding pocket identifies important intermolecular contacts that mediate binding. As seen experimentally, the complex with JG-03-14 (3,5-dibromo-4-(3,4-dimethoxyphenyl)-1H-pyrrole-2-carboxylic acid ethyl ester) is the most stable. These results illuminate the binding process and should be valuable in the design of new pyrrole-based colchicine site inhibitors as these compounds have very accessible syntheses.     

Curcumin binds in silico to anti-cancer drug target enzyme MMP-3 (human stromelysin-1) with affinity comparable to two known inhibitors of the enzyme

In silico interaction of curcumin with the enzyme MMP-3 (human stromelysin-1) was studied by molecular docking using AutoDock 4.2 as the docking software application. AutoDock 4.2 software serves as a valid and acceptable docking application to study the interactions of small compounds with proteins. Interactions of curcumin with MMP-3 were compared to those of two known inhibitors of the enzyme, PBSA and MPPT. The calculated free energy of binding (ΔG binding) shows that curcumin binds with affinity comparable to or better than the two known inhibitors. Binding interactions of curcumin with active site residues of the enzyme are also predicted. Curcumin appears to bind in an extendended conformation making extensive VDW contacts in the active site of the enzyme. Hydrogen bonding and pi-pi interactions with key active site residues is also observed. Thus, curcumin can be considered as a good lead compound in the development of new inhibitors of MMP-3 which is a potential target of anticancer drugs. The results of these studies can serve as a starting point for further computational and experimental studies.     

Screening for the selective inhibitors of MMP-9 from natural products based on pharmacophore modeling and molecular docking in combination with bioassay experiment,hybrid QM/MM calculation,and MD simulation

Matrix metalloproteinase-9 (MMP-9) has been considered as an attractive target involving cancer therapy. In this study, the 3D QSAR pharmacophore model of MMP-9 inhibitors is built, and its reliability is subsequently validated based on different methods. The built pharmacophore model consists of the four chemical features, including two hydrogen bond acceptors (HBA), one hydrophobic (HY), and one ring aromatic (RA). Among them, both HY and RA are found to be especially important features because they involve the interactions of inhibitors with the S1′ pocket of MMP-9, which determines the selectivity of MMP-9 inhibitors. By combining pharmacophore model with molecular docking, the virtual screening is carried out to identify the selective MMP-9 inhibitors from natural products. The four potential selective MMP-9 inhibitors of natural products are found. One of them was used to carry out the bioassay experiment inhibiting MMP-9, and the estimated IC₅₀ value of only 26.94 µM clearly shows its strongly inhibitory activity; besides, both the hybrid quantum mechanics/molecular mechanics (QM/MM) calculation and the molecular dynamics simulation are performed to examine the reliability regarding the binding mode of this inhibitor with MMP-9 active sites predicted by molecular docking. All the screened four natural products are found to well bind with the MMP-9 active sites by different kinds of interactions. Finally, the ADMET properties of screened four natural products are assessed. These screened MMP-9 inhibitors of natural products could be used as the lead compounds to perform structural modifications and optimizations in the future work.

Communicated by Ramaswamy H. Sarma     

Affinity of aporphines for the human 5-HT2A receptor: Insights from homology modeling and molecular docking studies

Analogs of nantenine were docked into a modeled structure of the human 5-HT_2A receptor using ICM Pro, GLIDE, and GOLD docking methods. The resultant docking scores were used to correlate with observed in vitro apparent affinity (K_e) data. The GOLD docking algorithm when used with a homology model of 5-HT_2A, based on a bovine rhodopsin template and built by the program MODELLER, gives results which are most in agreement with the in vitro results. Further analysis of the docking poses among members of a C1 alkyl series of nantenine analogs, indicate that they bind to the receptor in a similar orientation, but differently than nantenine. Besides an important interaction between the protonated nitrogen of the C1 alkyl analogs and residue Asp155, we identified Ser242, Phe234, and Gly238 as key residues responsible for the affinity of these compounds for the 5-HT_2A receptor. Specifically, the ability of some of these analogs to establish a H-bond with Ser242 and hydrophobic interactions with Phe234 and Gly238 appears to explain their enhanced affinity as compared to nantenine.     

Mechanistic insights into the inhibition of prostate specific antigen by beta-lactam class compounds

Prostate Specific Antigen (PSA) is a biomarker used in the diagnosis of prostate cancer and to monitor therapeutic response. However, its precise role in prostate carcinogenesis and metastasis remains largely unknown. A number of studies arguing in the favor of an active role of PSA in prostate cancer development and progression have implicated this serine protease in the release and activation of growth factors such as insulin-like growth factor 1 (IGF1) through cleavage of insulin like growth factor binding protein 3 and Transforming Growth Factor beta (TGF-beta) through cleavage of Latent TGF-beta. In contrast, other studies suggest that PSA activity might hinder tumor development and progression. In light of these contradictory findings, efficient inhibitors of PSA are needed for exploring its biological role in tumor development and metastasis. Towards the goal of developing potent inhibitors of PSA, we have explored the molecular mechanism of a series of beta-lactam based compounds on binding to PSA using activity assays, matrix assisted laser desorption ionization with a time-of-flight mass spectrometry, and GOLD docking methodology. The mass spectrometry experiments and the activity assays confirmed the time-dependent and covalent nature of beta-lactam binding. To gain insights on the reaction intermediates at the molecular level, we docked beta-lactam inhibitors to a homology modeled PSA using the GOLD docking program in noncovalent and covalent binding modes. The docking studies elucidated the molecular details of the early noncovalent Michaelis complex, the acyl-enzyme covalent complex, and the nature of conformational reorganization required for the long term stability of the covalent complex. Additionally, the molecular basis for the effect of stereochemistry of the lactam ring on the inhibitory potency was elucidated through docking of beta-lactam enantiomers. As a validation of our docking methodology, two novel enantiomers were synthesized and evaluated for their inhibitory potency using fluorogenic substrate based activity assays. Additionally, cis enantiomers of eight beta-lactam compounds reported in a previous study were docked and their GOLD scores and binding modes were analyzed in order to assess the general applicability of our docking results. The close agreement of our docking results with the experimental data validates the mechanistic insights revealed through the docking studies and paves the way for the design and development of potent and specific inhibitors of PSA.     

A computational docking study for prediction of binding mode of diospyrin and derivatives: Inhibitors of human and leishmanial DNA topoisomerase-I

A computational approach was utilized to study the relative binding modes of diospyrin (bisnaphthoquinonoid) with the crystal structure of human DNA-TopoI and the recently reported Leishmania donavani DNA-TopoI. Additionally, the binding site interactions of amino derivatives of diospyrin with human TopoI were studied extensively. Based on the docking results, binding modes of diospyrin with the human and leishmanial TopoI catalytic core were predicted. The parallel use of two efficient and predictive docking programs, GOLD and Ligandfit, allowed mutual validation of the predicted binding poses. A reasonably good correlation coefficient between the calculated docking scores and the experimentally determined cytotoxicity helped in validating the docking method. Furthermore, a structure-based pharmacophore model was developed for L. donavani DNA-TopoI inhibition which helped in elucidating the topological and spatial requirements of the ligand-receptor interactions. This study provides an understanding of the structural basis of ligand binding to the topoisomerase receptor, which may be used for the structure-based design of potent and novel ligands for anticancer and antileishmanial therapy. To our knowledge, this is the first report of a binding mode exploration study for diospyrin and its derivatives as inhibitors of the leishmanial and human TopoI enzymes.     

Type I collagen stabilization of matrix metalloproteinase-2

The activity of matrix metalloproteinase-2 (MMP-2) is regulated stringently on the posttranslational level. MMP-2 efficiently undergoes autolysis into inactive polypeptides in vitro, prompting the hypothesis that MMP-2 autolysis may function as an alternative mechanism for posttranslational control of MMP-2 in vivo. Moreover, MMP-2 binds to intact type I collagen fibrils; however, the functional consequences of this interaction have not been fully elucidated. To test the hypothesis that MMP-2 binding to type I collagen functions as a positive regulator of MMP-2 proteolytic potential, the effect of type I collagen on MMP-2 activity, inhibition by tissue inhibitor of metalloproteinase-2 (TIMP-2), and enzyme stability was examined. Here, we report that purified MMP-2 binds but does not cleave intact type I collagen. The presence of type I collagen affects neither enzymatic activity against a quenched fluorescent peptide substrate nor the kinetics of inhibition by TIMP-2. However, MMP-2 is stabilized from autolysis in the presence of type I collagen, but not by elastin, fibrinogen, or laminin. These data provide biochemical evidence that MMP-2 exosite interactions with type I collagen may function in the posttranslational control of MMP-2 activity by reducing the rate of autolytic inactivation.     

Molecular modeling of non-covalent binding of homochiral (3S,3'S)-astaxanthin to matrix metalloproteinase-13 (MMP-13)

Inhibitors for matrix metalloproteinases (MMPs) are under investigation for the treatment of various important chronic illnesses, including cancer, arthritis, and cardiovascular disease (CVD). In particular, MMP-13 is currently being probed as a potential key target in CVD and malignant disease due to its documented effects on extracellular matrix (ECM) remodeling, important in the pathophysiology of these diseases. Within the family of related mammalian MMP enzymes, MMP-13 possesses a large hydrophobic binding pocket relative to that of other MMPs. Homochiral astaxanthin (3S,3'S-AST; 3S,3'S-dihydroxy-beta,beta-carotene-4,4'-dione), an important antioxidant and anti-inflammatory xanthophyll carotenoid, is an active metabolite of several novel soft drugs in clinical development; it is also extensively used and tested as a human nutraceutical. In the current study, the prediction of the geometry and energetics of its binding to human MMP-13 was conducted with molecular modeling. The method used was found to predict the energy of binding of known ligands of MMP-13 with great precision. Blind docking using the whole protein target was then used in order to identify the possible binding site(s) of AST. AST was predicted to bind at several sites in close proximity to the active center. Subsequent analyses focused on the binding site at the atomic (i.e., amino acid sequence) level suggested that AST can bind to MMP-13 with high affinity and favorable energetics. Therefore, the modeling study predicts potential direct enzyme-inhibitory activity of AST against MMP-13, a behavior that may be exploited in mammalian systems in which pathological upregulation of MMP activity is paramount.     

Synthesis,antibacterial studies,and molecular modeling studies of 3,4-dihydropyrimidinone compounds

The syntheses of dihydropyrimidinones (DHPMs) using solvent-free grindstone chemistry method. All the synthesized compounds exhibited significant activity against pathogenic bacteria. The current effort has been developed to obtain new DHPM derivatives that focus on the bacterial ribosomal A site RNA as a drug target. Molecular docking simulation analysis was applied to confirm the target specificity of DHPMs. The crystal structure of bacterial 16S rRNA and human 40S rRNA was taken as receptors for docking. Finally, the docking score, binding site interaction analysis revealed that DHPMs exhibit more specificity towards 16S rRNA than known antibiotic amikacin. Accordingly, targeting the bacterial ribosomal A site RNA with potential drug leads promises to overcome the bacterial drug resistance. Even though, anti-neoplastic activities of DHPMs were also predicted through prediction of activity spectra for substances (PASS) tool. Further, the results establish that the DHPMs can serve as perfect leads against bacterial drug resistance.     

In silico screening and biological evaluation of inhibitors of Src-SH3 domain interaction with a proline-rich ligand

Src signalling and transduction are directly involved in cell growth, cell cycle, malignant transformation and cell migration, providing therapeutic opportunities through inhibition of Src. Here we report virtual screening for novel compounds that inhibit the Src-SH3 protein-protein interaction with a proline-rich peptide ligand. Computational docking of the ZINC compound database was performed using GOLD. Top-scoring compounds were assayed using a fluorescence polarization-based assay. A benzoquinoline derivative showed micromolar inhibition of binding between Src-SH3 and the proline-rich peptide. Several analogues were subsequently assayed showing the requirement of a linker between the benzoquinoline and phenyl rings, and electron donating substituents on the phenyl ring.     

New azoles with antifungal activity: Design, synthesis, and molecular docking

In order to search for many target compounds with excellent activities, a series of 1-(1H-1,2,4-triazol-1-yl)-2-(2,4-difluoro-phenyl)-3-[(4-substituted phenyl)-piperazin-1-yl]-propan-2-ols were designed, synthesized, and evaluated as antifungal agents. Results of preliminary antifungal tests against eight human pathogenic fungi in vitro showed that all the title compounds exhibited excellent activities with broad spectrum. Moreover, a molecular model for the binding between 5a and the active site of CACYP51 was provided based on the computational docking results.     

Discovery of potential pancreatic cholesterol esterase inhibitors using pharmacophore modelling, virtual screening, and optimization studies

Pancreatic cholesterol esterase (CEase) is a serine hydrolase involved in the hydrolysis of variety of lipids and transport of free cholesterol. In this study, pharmacophore hypotheses based on known inhibitors were generated using common feature pharmacophore generation protocol available in Discovery Studio program. The best pharmacophore model containing two hydrogen bond acceptor and three hydrophobic features was selected and validated. It was further used in screening three diverse chemical databases. Hit compounds were subjected to drug-likeness and molecular docking studies. Four hits, namely SEW00846, NCI0040784, GK03167, and CD10645, were selected based on the GOLD fitness score and interaction with active site amino acids. All hit compounds were further optimized to improve their binding in the active site. The optimized compounds were found to have improved binding at the active site. Strongly binding optimized hits at the active site can act as virtual leads in potent CEase inhibitor designing.     

3D structure generation,virtual screening and docking of human Ras-associated binding (Rab3A) protein involved in tumourigenesis

Rab3A is expressed predominantly in brain and synaptic vesicles. Rab3A is involved specifically in tethering and docking of synaptic vesicles prior to fusion which is a critical step in regulated release of neurotransmitters. The precise function of Rab3A is still not known. However, up-regulation of Rab3A has been reported in malignant neuroendocrine and breast cancer cells. In the present study, the structure of Rab3A protein was generated using MODELLER 9v8 software. The modeled protein structure was validated and subjected to molecular docking analyses. Docking with GTP was carried out on the binding site of Rab3A using GOLD software. The Rab3A-GTP complex has best GOLD fitness value of 77.73. Ligplot shows hydrogen bondings (S16, S17, V18, G19, K20, T21, S22, S31, T33, A35, S38, T39 and G65) and hydrophobic interacting residues (F25, F32, P34, F36, V37, D62 and A64) with the GTP ligands in the binding site of Rab3A protein. Here, the ligand molecules of NCI diversity set II from the ZINC database against the active site of the Rab3A protein were screened. For this purpose, the incremental construction algorithm of GLIDE and the genetic algorithm of GOLD were used. Docking results were analyzed for top ranking compounds using a consensus scoring function of X-Score to calculate the binding affinity and Ligplot was used to measure protein–ligand interactions. Five compounds which possess good inhibitory activity and may act as potential high affinity inhibitors against Rab3A active site were identified. The top ranking molecule (ZINC13152284) has a Glide score of ?6.65 kcal/mol, X-Score of ?3.02 kcal/mol and GOLD score of 64.54 with 03 hydrogen bonds and 09 hydrophobic contacts. This compound is thus a good starting point for further development of strong inhibitors.     

Macrocyclic inhibitors of the malarial aspartic proteases plasmepsin I, II, and IV

The first macrocyclic inhibitor of the Plasmodium falciparum aspartic proteases plasmepsin I, II, and IV with considerable selectivity over the human aspartic protease cathepsin D has been identified. A series of macrocyclic compounds were designed and synthesized. Cyclizations were accomplished using ring-closing metathesis with the second generation Grubbs catalyst. These compounds contain either a 13-membered or a 16-membered macrocycle and incorporate a 1,2-dihydroxyethylene as transition state mimicking unit. The binding mode of this new class of compounds was predicted with automated docking and molecular dynamics simulations, with an estimation of the binding affinities through the linear interaction energy (LIE) method.