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.     

X-ray crystal structures,density functional theory and docking on deacetylase enzyme for antiproliferative activity of hispolon derivatives on HCT116 colon cancer

The antiproliferative action of hispolon derivatives is stronger than that of related curcumin against several tumor cell lines. Hispolon size, smaller than curcumin, fits better than curcumin into the active site of HDAC6, an enzyme involved in deacetylation of lysine residues. HDACs are considered potential targets for tumor drug discovery and hydroxamates are known inhibitors of HDACs. One of them, SAHA (Vorinostat) is used in clinical studies. Investigations into possible mechanisms for hispolon derivatives active against the HCT116 colon tumor cell line are done after examining the structural results obtained from hispolon X-ray crystal structures as well as performing associated computational docking and Density Functional Theory techniques on HDAC6. These studies show preference for the HDAC6 active site by chelating the Zn center, in contrast with other ineffective hispolon derivatives, that establish only a single bond to the metal center. Structure activity relationships make clear that hydrogenation of the hispolon bridge also leads to single bond (non chelate) hispolon-Zn binding, and consistently nullifies the antiproliferative action against HCT116 tumor.     

Computer-aided molecular design of 1H-imidazole-2,4-diamine derivatives as potential inhibitors of <Emphasis Type="Italic">Plasmodium falciparum</Emphasis> DHFR enzyme

Design and discovery of new potential inhibitors of Plasmodium falciparum dihydrofolate reductase (PfDHFR), equally active against both the wild-type and mutant strains, is urgently needed. In this study, a computer-aided molecular design approach that involved ab initio molecular orbital and density functional theory calculations, along with molecular electrostatic potential analysis, and molecular docking studies was employed to design 15 1H-imidazole-2,4-diamine derivatives as potential inhibitors of PfDHFR enzyme. Visual inspection of the binding modes of the compounds demonstrated that they all interact, via H-bond interactions, with key amino acid residues (Asp54, Ileu/Leu164, Asn/Ser108 and Ile14) similar to those of WR99210 (3) in the active site of the enzymes used in the study. These interactions are known to be essential for enzyme inhibition. These compounds showed better or comparable binding affinities to that of the bound ligand (WR99210). In silico toxicity predictions, carried out using TOPKAT software, also indicated that the compounds are non-toxic.     

Efficient docking of peptides to proteins without prior knowledge of the binding site

Reliability in docking of ligand molecules to proteins or other targets is an important challenge for molecular modeling. Applications of the docking technique include not only prediction of the binding mode of novel drugs, but also other problems like the study of protein-protein interactions. Here we present a study on the reliability of the results obtained with the popular AutoDock program. We have performed systematical studies to test the ability of AutoDock to reproduce eight different protein/ligand complexes for which the structure was known, without prior knowledge of the binding site. More specifically, we look at factors influencing the accuracy of the final structure, such as the number of torsional degrees of freedom in the ligand. We conclude that the Autodock program package is able to select the correct complexes based on the energy without prior knowledge of the binding site. We named this application blind docking, as the docking algorithm is not able to "see" the binding site but can still find it. The success of blind docking represents an important finding in the era of structural genomics.     

Protein-Ligand interaction studies on 2, 4, 6- trisubstituted triazine derivatives as anti-malarial DHFR agents using AutoDock

The dihydrofolate reductase (DHFR) domain of P. falciparum is one of the few well defined targets in malarial chemotherapy. The enzyme catalyzes the nicotinamide adenine dinucleotide phosphate (NADPH) dependent reduction of dihydrofolate to tetrahydrofolate. Protein-ligand interactions were studied using DHFR protein 2BL9, extracted from PDB to evaluate the strength of affinity of various molecules towards ligand binding site and to study the extent of correlation between experimental values and computational dock scores. AutoDock runs resulted in binding energy scores from -7.14 to -10.72 kcal/mol. Among the five inhibitors (Bioorganic and Medicinal Chemistry Letters 15 2005 531-533) selected for docking studies, an excellent correlation was observed in all cases, for instance, experimentally reported most active molecule 2a (MIC: 1μg/ml) showed a high dock score (-10.72 kcal/mol) than the remaining inhibitors. Therefore, molecular docking using AutoDock suggests the importance of evaluating the prediction accuracy of various molecules as evidenced by a correlation coefficient of 0.961 between experimental activities and AutoDock binding energies.     

Toward the virtual screening of Cdc25A phosphatase inhibitors with the homology modeled protein structure

Cdc25 phosphatases have been considered as attractive drug targets for anticancer therapy due to the correlation of their overexpression with a wide variety of cancers. As a method for the discovery of novel inhibitors of Cdc25 phosphatases, we have evaluated the computer-aided drug design protocol involving the homology modeling of Cdc25A and virtual screening with the two docking tools: FlexX and the modified AutoDock program implementing the effects of ligand solvation in the scoring function. The homology modeling with the X-ray crystal structure of Cdc25B as a template provides a high-quality structure of Cdc25A that enables the structure-based inhibitor design. Of the two docking programs under consideration, AutoDock is found to be more accurate than FlexX in terms of scoring putative ligands. A detailed binding mode analysis of the known inhibitors shows that they can be stabilized in the active site of Cdc25A through the simultaneous establishment of the multiple hydrogen bonds and the hydrophobic interactions. The present study demonstrates the usefulness of the modified AutoDock program as a docking tool for virtual screening of new Cdc25 phosphatase inhibitors as well as for binding mode analysis to elucidate the activities of known inhibitors. Figure Structures and available IC50 values (in μM) of the twenty Cdc25 phosphatase inhibitors seeded in docking library     

The molecular basis of urokinase inhibition: from the nonempirical analysis of intermolecular interactions to the prediction of binding affinity

Urokinase-type plasminogen activator (uPA) is a trypsin-like serine protease that plays a crucial role in angiogenesis process. In addition to its physiological role in healthy organisms, angiogenesis is extremely important in cancer growth and metastasis, resulting in numerous attempts to understand its control and to develop new approaches to anticancer therapy. The α-aminoalkylphosphonate diphenyl esters are well known as highly efficient serine protease inhibitors. However, their mode of binding has not been verified experimentally in details. For a group of average and potent phosphonic inhibitors of urokinase, flexible docking calculations were performed to gain an insight into the active site interactions responsible for observed enzyme inhibition. The docking results are consistent with the previously suggested mode of inhibitors binding. Subsequently, rigorous ab initio study of binding energy was carried out, followed by its decomposition according to the variation–perturbation procedure to reveal stabilization energy constituents with clear physical meaning. Availability of the experimental inhibitory activities and comparison with theoretical binding energy allows for the validation of theoretical models of inhibition, as well as estimation of the possible potential for binding affinity prediction. Since the docking results accompanied by molecular mechanics optimization suggested that several crucial active site contacts were too short, the optimal distances corresponding to the minimum ab initio interaction energy were also evaluated. Despite the deficiencies of force field-optimized enzyme-inhibitor structures, satisfactory agreement with experimental inhibitory activity was obtained for the electrostatic interaction energy, suggesting its possible application in the binding affinity prediction. Figure The comparison of an arrangement of inhibitors within uPA active. Amino acid residues form S1 pocket that binds the variable part of inhibitor molecules     

Effect of hesperetin on tyrosinase: inhibition kinetics integrated computational simulation study

Tyrosinase inhibitors have potential applications in medicine, cosmetics and agriculture to prevent hyperpigmentation or browning effects. Some of the flavonoids mostly found in herbal plants and fruits are revealed as tyrosinase inhibitors. We studied the inhibitory effects of one such flavonoid, hesperetin, on mushroom tyrosinase using inhibition kinetics and computational simulation. Hesperetin reversibly inhibited tyrosinase in a competitive manner with K_i = 4.03 ± 0.26 mM. Measurements of ANS-binding fluorescence showed that hesperetin induced the hydrophobic disruption of tyrosinase. For further insight, we used the docking algorithms to simulate binding between tyrosinase and hesperetin. Simulation was successful (binding energies for Dock6.3: −34.41 kcal/mol and for AutoDock4.2: −5.67 kcal/mol) and showed that a copper ion coordinating with 3 histidine residues (HIS61, HIS85, and HIS259) within the active site pocket was chelated via hesperetin binding. Our study provides insight into the inhibition of tyrosinase in response to flavonoids. A combination of inhibition kinetics and computational prediction may facilitate the identification of potential natural tyrosinase inhibitors such as flavonoids and the prediction of their inhibitory mechanisms.     

Docking studies of nickel-peptide deformylase (PDF) inhibitors: exploring the new binding pockets

The binding modes of a series of known activity inhibitors docking to Peptide deformylase (PDF) have been studied using molecular docking software AutoDock3.0.5. In this study, good correlation (R(2)=0.894) between calculated binding energies and experimental inhibitory activities is obtained. We find that some shallow pockets near the known active pocket are very important which can accommodate the side-chains of the inhibitor. Moreover, a new binding pocket is also explored. All these may provide something useful for designing the potent inhibitors.     

Structure modeling of novel DNA glycosylase enzyme from oral pathogen Streptococcus sanguinis

The novel 3-methyladenine DNA glycosylase enzyme from oral pathogen Streptococcus sanguinisin involves in DNA repair mechanisms and participates in base excision repair. Its 3D structure is still unknown which may be a potential drug target, therefore here we proposed its putative 3D structure by homology modeling approach. EsyPred3d software produced more precise modeled structure as compare to Swiss model software. The modeled structure was further verified by PROCHECK analysis and subjected to functional site prediction servers for active site residues prediction. The functional site was further validated by molecular docking approach with ligand EDA (3- [2- Deoxyribofuranosyl] - 3H- 1, 3, 4, 5A, 8-Pentaaza- Asindacene-5- monophosphate) from 1F4R. The EDR docked at the cavity of modeled structure of 3-methyladenine DNA glycosylase enzyme with highest Patchdock score of 3966 and lowest Autodock 4 docking energy of -10.30 Kcal/mol. The YA51, LA105, RA107 residues are surrounding the EDA and matching with ligand binding residues predicted by PROFUNC server.     

Homology modeling of human serum paraoxonase1 and its molecular interaction studies with aspirin and cefazolin

Human serum paraoxonase1 (HuPON1) belongs to the family of A-esterases (EC.3.1.8.1). It is associated with HDL particle and prevents atherosclerosis by cleaving lipid hydroperoxides and other proatherogenic molecules of oxidized low density lipoproteins (LDL). Since the precise structure of HuPON1 is not yet available, the structure-function relationship between HuPON1 and activators/inhibitors is still unknown. Therefore, a theoretical model of HuPON1 was generated using homology modelling and precise molecular interactions of an activator aspirin and an inhibitor cefazolin with PON1 were studied using Autodock software. The ligand binding residues were found to be similar to the predicted active site residues. Both cefazolin and aspirin were found to dock in the vicinity of the predicted active sites of PON1; cefazolin bound at residues N166, S193 and Y71, while aspirin at residues N309, I310 and L311. Binding region in the PON1 by prediction (3D2GO server) and docking studies provide useful insight into mechanism of substrate and inhibitor binding to the enzyme active site.     

An Integrated Study of Tyrosinase Inhibition by Rutin: Progress using a Computational Simulation

Abstract

Tyrosinase inhibition studies have recently gained the attention of researchers due to their potential application values. We simulated docking (binding energies for AutoDock Vina: ?9.1 kcal/mol) and performed a molecular dynamics simulation to verify docking results between tyrosinase and rutin. The docking results suggest that rutin mostly interacts with histidine residues located in the active site. A 10 ns molecular dynamics simulation showed that one copper ion at the tyrosinase active site was responsible for the interaction with rutin. Kinetic analyses showed that rutin-mediated inactivation followed a first-order reaction and mono- and biphasic rate constants occurred with rutin. The inhibition was a typical competitive type with K_i = 1.10 ± 0.25 mM. Measurements of intrinsic and ANS-binding fluorescences showed that rutin showed a relatively strong binding affinity for tyrosinase and one possible binding site that could be a copper was detected accompanying with a hydrophobic exposure of tyrosinase. Cell viability testing with rutin in HaCaT keratinocytes showed that no toxic effects were produced. Taken together, rutin has the potential to be a potent antipigment agent. The strategy of predicting tyrosinase inhibition based on hydroxyl group number and computational simulation may prove useful for the screening of potential tyrosinase inhibitors.     

An integrated study of tyrosinase inhibition by rutin: progress using a computational simulation

Tyrosinase inhibition studies have recently gained the attention of researchers due to their potential application values. We simulated docking (binding energies for AutoDock Vina: -9.1 kcal/mol) and performed a molecular dynamics simulation to verify docking results between tyrosinase and rutin. The docking results suggest that rutin mostly interacts with histidine residues located in the active site. A 10 ns molecular dynamics simulation showed that one copper ion at the tyrosinase active site was responsible for the interaction with rutin. Kinetic analyses showed that rutin-mediated inactivation followed a first-order reaction and mono- and biphasic rate constants occurred with rutin. The inhibition was a typical competitive type with K(i) = 1.10±0.25 mM. Measurements of intrinsic and ANS-binding fluorescences showed that rutin showed a relatively strong binding affinity for tyrosinase and one possible binding site that could be a copper was detected accompanying with a hydrophobic exposure of tyrosinase. Cell viability testing with rutin in HaCaT keratinocytes showed that no toxic effects were produced. Taken together, rutin has the potential to be a potent anti-pigment agent. The strategy of predicting tyrosinase inhibition based on hydroxyl group number and computational simulation may prove useful for the screening of potential tyrosinase inhibitors.     

Modeling of human M1 aminopeptidases for in silico screening of potential Plasmodium falciparum alanine aminopeptidase (PfA-M1) specific inhibitors

Plasmodium falciparum alanine M1-aminopeptidase (PfA-M1) is a validated target for anti-malarial drug development. Presence of significant similarity between PfA-M1 and human M1-aminopeptidases, particularly within regions of enzyme active site leads to problem of non-specificity and off-target binding for known aminopeptidase inhibitors. Molecular docking based in silico screening approach for off-target binding has high potential but requires 3D-structure of all human M1-aminopeptidaes. Therefore, in the present study 3D structural models of seven human M1-aminopeptidases were developed. The robustness of docking parameters and quality of predicted human M1-aminopeptidases structural models was evaluated by stereochemical analysis and docking of their respective known inhibitors. The docking scores were in agreement with the inhibitory concentrations elucidated in enzyme assays of respective inhibitor enzyme combinations (r2≈0.70). Further docking analysis of fifteen potential PfA-M1 inhibitors (virtual screening identified) showed that three compounds had less docking affinity for human M1-aminopeptidases as compared to PfA-M1. These three identified potential lead compounds can be validated with enzyme assays and used as a scaffold for designing of new compounds with increased specificity towards PfA-M1.     

Theoretical studies on the interaction of biphenyl inhibitors with Mycobacterium tuberculosis protein tyrosine phosphatase MptpB

MptpB is an essential secreted virulence factor for M. tuberculosis. Inhibition of MptpB impairs mycobacterial survival in host macrophages and thus helps reduce tuberculosis infections. However, the binding mode of the biphenyl inhibitors, which are known as some of the most potent MptpB inhibitors, remains unclear. In this study, to understand the interactions between biphenyl inhibitors and MptpB, docking and molecular dynamics simulations were carried out using AutoDock and GROMACS softwares. Calculation results show that all the biphenyl inhibitors can be docked to the binding site of MptpB, with the acid warheads forming a hydrogen bond network at the active site. But the binding modes of other terminals of these inhibitors are different. The cyclohexyl and trifluoromethyl substituents at R1 and R2 sites are necessary for the inhibitors to adopt their double-site binding mechanism. The estimated binding affinities are basically consistent with the experimental results. MD simulations show that these binding complexes display different stability.     

Computational analyses of curcuminoid analogs against kinase domain of HER2

Background

Human epidermal growth factor receptor 2 (HER2) has an important role in cancer aggressiveness and poor prognosis. HER2 has been used as a drug target for cancers. In particular, to effectively treat HER2-positive cancer, small molecule inhibitors were developed to target HER2 kinase. Knowing that curcumin has been used as food to inhibit cancer activity, this study evaluated the efficacy of natural curcumins and curcumin analogs as HER2 inhibitors using in vitro and in silico studies. The curcumin analogs considered in this study composed of 4 groups classified by their core structure, β-diketone, monoketone, pyrazole, and isoxazole.

Results

In the present study, both computational and experimental studies were performed. The specificity of curcumin analogs selected from the docked results was examined against human breast cancer cell lines. The screened curcumin compounds were then subjected to molecular dynamics simulation study. By modifying curcumin analogs, we found that protein-ligand affinity increases. The benzene ring with a hydroxyl group could enhance affinity by forming hydrophobic interactions and the hydrogen bond with the hydrophobic pocket. Hydroxyl, carbonyl or methoxy group also formed hydrogen bonds with residues in the adenine pocket and sugar pocket of HER2-TK. These modifications could suggest the new drug design for potentially effective HER2-TK inhibitors. Two outstanding compounds, bisdemethylcurcumin (AS-KTC006) and 3,5-bis((E)-3,4-dimethoxystyryl)isoxazole (AS-KTC021 ),were well oriented in the binding pocket almost in the simulation time, 30 ns. This evidence confirmed the results of cell-based assays and the docking studies. They possessed more distinguished interactions than known HER2-TK inhibitors, considering them as a promising drug in the near future.

Conclusions

The series of curcumin compounds were screened using a computational molecular docking and followed by human breast cancer cell lines assay. Both AS-KTC006 and AS-KTC021 could inhibit breast cancer cell lines though inhibiting of HER2-TK. The intermolecular interactions were confirmed by molecular dynamics simulation studies. This information would explore more understanding of curcuminoid structures and HER2-TK.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-261) contains supplementary material, which is available to authorized users.     

Kinetic and docking studies of phenol-based inhibitors of carbonic anhydrase isoforms I, II, IX and XII evidence a new binding mode within the enzyme active site

Carbonic anhydrases (CAs, EC 4.2.1.1) are inhibited by sulfonamides, inorganic anions, phenols, coumarins (acting as prodrugs) and polyamines. A novel class of CA inhibitors (CAIs), interacting with the CA isozymes I, II (cytosolic) and IX, XII (transmembrane, tumor-associated) in a different manner, is reported here. Kinetic measurements allowed us to identify hydroxy-/methoxy-substituted benzoic acids as well as di-/tri-methoxy benzenes as submicromolar-low micromolar inhibitors of the four CA isozymes. Molecular docking studies of a set of such inhibitors within CA I and II allowed us to understand the inhibition mechanism. This new class of inhibitors binds differently compared to all other classes of inhibitors known to date: they were found between the phenol-binding site and the coumarin-binding site, filling thus the middle of the enzyme cavity. They exploit different interactions with amino acid residues and water molecules from the CA active site compared to other classes of inhibitors, offering the possibility to design CAIs with an interesting inhibition profile compared to the clinically used sulfonamides/sulfamates.     

Exploration of Virtual Candidates for Human HMG-CoA Reductase Inhibitors Using Pharmacophore Modeling and Molecular Dynamics Simulations

3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGR) is a rate-controlling enzyme in the mevalonate pathway which involved in biosynthesis of cholesterol and other isoprenoids. This enzyme catalyzes the conversion of HMG-CoA to mevalonate and is regarded as a drug target to treat hypercholesterolemia. In this study, ten qualitative pharmacophore models were generated based on chemical features in active inhibitors of HMGR. The generated models were validated using a test set. In a validation process, the best hypothesis was selected based on the statistical parameters and used for virtual screening of chemical databases to find novel lead candidates. The screened compounds were sorted by applying drug-like properties. The compounds that satisfied all drug-like properties were used for molecular docking study to identify their binding conformations at active site of HMGR. The final hit compounds were selected based on docking score and binding orientation. The HMGR structures in complex with the hit compounds were subjected to 10 ns molecular dynamics simulations to refine the binding orientation as well as to check the stability of the hits. After simulation, binding modes including hydrogen bonding patterns and molecular interactions with the active site residues were analyzed. In conclusion, four hit compounds with new structural scaffold were suggested as novel and potent HMGR inhibitors.     

Understanding the conformational flexibility and electrostatic properties of curcumin in the active site of rhAChE via molecular docking,molecular dynamics,and charge density analysis

Acetylcholinesterase (AChE) is an important enzyme responsible for Alzheimer’s disease, as per report, keto-enol form of curcumin inhibits this enzyme. The present study aims to understand the binding mechanism of keto-enol curcumin with the recombinant human Acetylcholinesterase (rhAChE) from its conformational flexibility, intermolecular interactions, charge density distribution, and the electrostatic properties at the active site of rhAChE. To accomplish this, a molecular docking analysis of curcumin with the rhAChE was performed, which gives the structure and conformation of curcumin in the active site of rhAChE. Further, the charge density distribution and the electrostatic properties of curcumin molecule (lifted from the active site of rhAChE) were determined from the high level density functional theory (DFT) calculations coupled with the charge density analysis. On the other hand, the curcumin molecule was optimized (gas phase) using DFT method and further, the structure and charge density analysis were also carried out. On comparing the conformation, charge density distribution and the electrostatic potential of the active site form of curcumin with the corresponding gas phase form reveals that the above said properties are significantly altered when curcumin is present in the active site of rhAChE. The conformational stability and the interaction of curcumin in the active site are also studied using molecular dynamics simulation, which shows a large variation in the conformational geometry of curcumin as well as the intermolecular interactions.