Molecular modeling studies on nucleoside hydrolase from the biological warfare agent Brucella suis

Brucella suis is a dangerous biological warfare agent already used for military purposes. This bacteria cause brucellosis, a zoonosis highly infective and difficult to fight. An important selective target for chemotherapy against this disease is nucleoside hydrolase (NH), an enzyme still not found in mammals. We present here the first three-dimensional structure of B. suis NH (BsNH) and propose this enzyme as a molecular target to the drug design in the fight against brucellosis. In addition, we performed molecular docking studies, aiming to analyze the three-dimensional positioning of nine known inhibitors of Chritidia fasciculata NH (CfNH) in the active sites of BsNH and CfNH. We also analyzed the main interactions of some of these compounds inside the active site of BsNH and the relevant factors to biological activity. These results, together with further molecular dynamics (MD) simulations, pointed out to the most promising compound as lead for the design of potential inhibitors of BsNH. Most of the docking and MD results corroborated to each other and the docking results also suggested a good correlation with experimental data.     

Molecular modeling studies on nucleoside hydrolase from the biological warfare agent Brucella suis

Brucella suis is a dangerous biological warfare agent already used for military purposes. This bacteria cause brucellosis, a zoonosis highly infective and difficult to fight. An important selective target for chemotherapy against this disease is nucleoside hydrolase (NH), an enzyme still not found in mammals. We present here the first three-dimensional structure of B. suis NH (BsNH) and propose this enzyme as a molecular target to the drug design in the fight against brucellosis. In addition, we performed molecular docking studies, aiming to analyze the three-dimensional positioning of nine known inhibitors of Chritidia fasciculata NH (CfNH) in the active sites of BsNH and CfNH. We also analyzed the main interactions of some of these compounds inside the active site of BsNH and the relevant factors to biological activity. These results, together with further molecular dynamics (MD) simulations, pointed out to the most promising compound as lead for the design of potential inhibitors of BsNH. Most of the docking and MD results corroborated to each other and the docking results also suggested a good correlation with experimental data.     

Analysis of Bacillus anthracis nucleoside hydrolase via in silico docking with inhibitors and molecular dynamics simulation

As the enzyme nucleoside hydrolase (NH) is widely found in nature but has not yet been detected in mammals, it is considered an ideal target in the development of chemotherapy against parasitic diseases and bacterial infections like anthrax. Considering the risk that this biological warfare agent represents nowadays, the search for new drugs and new molecular targets in the development of chemotherapy against anthrax is imperative. On this basis, we performed docking studies of six known NH inhibitors at the active site of NH from Bacillus anthracis (BaNH). Subsequently, molecular dynamics (MD) simulations of these compounds inside BaNH were carried out in order to complement the docking studies and select the most promising compounds as leads for the design of potential BaNH inhibitors. Most of the docking and MD results obtained agreed well with each other and showed good correlation with experimental data.     

3D-QSAR,molecular docking,and molecular dynamic simulations for prediction of new Hsp90 inhibitors based on isoxazole scaffold

Heat shock protein 90(Hsp90), as a molecular chaperone, play a crucial role in folding and proper function of many proteins. Hsp90 inhibitors containing isoxazole scaffold are currently being used in the treatment of cancer as tumor suppressers. Here in the present studies, new compounds based on isoxazole scaffold were predicted using a combination of molecular modeling techniques including three-dimensional quantitative structure–activity relationship (3D-QSAR), molecular docking and molecular dynamic (MD) simulations. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were also done. The steric and electrostatic contour map of CoMFA and CoMSIA were created. Hydrophobic, hydrogen bond donor and acceptor of CoMSIA model also were generated, and new compounds were predicted by CoMFA and CoMSIA contour maps. To investigate the binding modes of the predicted compounds in the active site of Hsp90, a molecular docking simulation was carried out. MD simulations were also conducted to evaluate the obtained results on the best predicted compound and the best reported Hsp90 inhibitors in the 3D-QSAR model. Findings indicate that the predicted ligands were stable in the active site of Hsp90.     

Understanding microscopic binding of macrophage migration inhibitory factor with phenolic hydrazones by molecular docking, molecular dynamics simulations and free energy calculations

Macrophage migration inhibitory factor (MIF), an immunoregulatory protein, is a potential target for a number of inflammatory diseases. In the current work, the interactions between MIF and a series of phenolic hydrazones were studied by molecular docking, molecular dynamics (MD) simulations, binding free energy calculations, and binding energy decomposition analysis to determine the structural requirement for achieving favorable biological activity of phenolic hydrazones. First, molecular docking was used to predict the binding modes of inhibitors in the binding site of MIF. The good correlation between the predicted docking scores and the experimental activities shows that the binding conformations of the inhibitors in the active site of MIF are well predicted. Moreover, our results suggest that the flexibility of MIF is essential in ligand binding process. Then, MD simulations and MM/GBSA free energy calculations were employed to determine the dynamic binding process and compare the binding modes of the inhibitors with different activities. The predicted binding free energies given by MM/GBSA are not well correlated with the experimental activities for the two subsets of the inhibitors; however, for each subset, a good correlation between the predicted binding free energies and the experimental activities is achieved. The MM/GBSA free energy decomposition analysis highlights the importance of hydrophobic residues for the MIF binding of the studied inhibitors. Based on the essential factors for MIF-inhibitor interactions derived from the theoretical predictions, some derivatives were designed and the higher inhibitory activities of several candidates were confirmed by molecular docking studies. The structural insights obtained from our study are useful for designing potent inhibitors of MIF.     

Homology modeling of human CCR5 and analysis of its binding properties through molecular docking and molecular dynamics simulation

In this study, homology modeling, molecular docking and molecular dynamics simulation were performed to explore structural features and binding mechanism of some inhibitors of chemokine receptor type 5 (CCR5), and to construct a model for designing new CCR5 inhibitors for preventing HIV attachment to the host cell. A homology modeling procedure was employed to construct a 3D model of CCR5. For this procedure, the X-ray crystal structure of bovine rhodopsin (1F88A) at 2.80? resolution was used as template. After inserting the constructed model into a hydrated lipid bilayer, a 20ns molecular dynamics (MD) simulation was performed on the whole system. After reaching the equilibrium, twenty-four CCR5 inhibitors were docked in the active site of the obtained model. The binding models of the investigated antagonists indicate the mechanism of binding of the studied compounds to the CCR5 obviously. Moreover, 3D pictures of inhibitor-protein complex provided precious data regarding the binding orientation of each antagonist into the active site of this protein. One additional 20 ns MD simulation was performed on the initial structure of the CCR5-ligand 21 complex, resulted from the previous docking calculations, embedded in a hydrated POPE bilayer to explore the effects of the presence of lipid bilayer in the vicinity of CCR5-ligand complex. This article is part of a Special Issue entitled Protein translocation across or insertion into membranes.     

Viper venom hyaluronidase and its potential inhibitor analysis: a multipronged computational investigation

Viper venom hyaluronidase (VV-HYA) inhibitors have long been used as therapeutic agents for arresting the local and systemic effects caused during its envenomation. Henceforth, to understand its structural features and also to identify the best potential inhibitor against it the present computational study was undertaken. Structure-based homology modeling of VV-HYA followed by its docking and free energy-based ranking analysis of ligand, the MD simulations of the lead complex was also performed. The sequence analysis and homology modeling of VV-HYA revealed a distorted (β/α)₈ folding as in the case of hydrolases family of proteins. Molecular docking of the resultant 3D structure of VV-HYA with known inhibitors (compounds 1–25) revealed the importance of molecular recognition of hotspot residues (Tyr 75, Arg 288, and Trp 321) other than that of the active site residues. It also revealed that Trp 321 of VV-HYA is highly important for mediating π–π interactions with ligands. In addition, the molecular docking and comparative free energy binding analysis was investigated for the VV-HYA inhibitors (compounds 1–25). Both molecular docking and relative free energy binding analysis clearly confirmed the identification of sodium chromoglycate (compound 1) as the best potential inhibitor against VV-HYA. Molecular dynamics simulations additionally confirmed the stability of their binding interactions. Further, the information obtained from this work is believed to serve as an impetus for future rational designing of new novel VV-HYA inhibitors with improved activity and selectivity.     

Exploring docking methods for virtual screening: application to the identification of neuraminidase and Ftsz potential inhibitors

Virtual screenings based on molecular docking play a major role in medicinal chemistry for the identification of new bioactive molecules. For this purpose, several docking methods can be used. Here, using Arguslab as software and a Gold Platinum subset library of commercially available compounds from Asinex, two docking methods associated to the scoring function Ascore were employed to investigate virtual screenings. One method is based on a genetic algorithm and the other based on a shape-based method. As case studies, both docking techniques were explored by targeting the PC190723 binding site of FtsZ protein from Staphylococcus aureus and the active site of N8 neuraminidase from Influenza virus. Following four docking sequences for each docking engine, the genetic algorithm led to multiple docking results, whereas the shape-based method gave reproducible results. The present study shows that the stochastic nature of the genetic algorithm will require the biological evaluation of more compounds than the shape-based method. This study showed that both methods are complementary and also led to the identification of neuraminidase and FtsZ potential inhibitors.     

Novel virtual lead identification in the discovery of hematopoietic cell kinase (HCK) inhibitors: application of 3D QSAR and molecular dynamics simulation

High level of hematopoietic cell kinase (Hck) is associated with drug resistance in chronic myeloid leukemia. Additionally, Hck activity has also been connected with the pathogenesis of HIV-1 and chronic obstructive pulmonary disease. In this study, three-dimensional (3D) QSAR pharmacophore models were generated for Hck based on experimentally known inhibitors. A best pharmacophore model, Hypo1, was developed with high correlation coefficient (0.975), Low RMS deviation (0.60) and large cost difference (49.31), containing three ring aromatic and one hydrophobic aliphatic feature. It was further validated by the test set (r?=?0.96) and Fisher’s randomization method (95%). Hypo 1 was used as a 3D query for screening the chemical databases, and the hits were further screened by applying Lipinski’s rule of five and ADMET properties. Selected hit compounds were subjected to molecular docking to identify binding conformations in the active site. Finally, the appropriate binding modes of final hit compounds were revealed by molecular dynamics (MD) simulations and free energy calculation studies. Hence, we propose the final three hit compounds as virtual candidates for Hck inhibitors.     

Synthesis and biological evaluation of novel quinazoline-triazole hybrid compounds with potential use in Alzheimer’s disease

A library of twelve quinazoline-triazole hybrid compounds were designed, synthesized and evaluated as a novel class of acetylcholinesterase inhibitors to treat Alzheimer’s disease (AD). The biological assay results demonstrated the ability of several hybrid compounds to inhibit AChE enzyme (IC₅₀ range = 0.2–83.9 µM). To understand the high potential activity of these compounds, molecular docking simulations were performed to get better insights into the mechanism of binding of quinazoline-triazole hybrid compounds. As expected, compounds 8a and 9a-b bind to both catalytic anionic site (CAS) and peripheral anionic site (PAS) in the active site of AChE enzyme, which implicates that these compounds could act as dual binding site inhibitors. These compounds were not cytotoxic and they also displayed appropriated physicochemical as well as pharmacokinetic profile to be developed as novel anti-AD drug candidates.     

Novel inhibitors of anthrax edema factor

Several pathogenic bacteria produce adenylyl cyclase toxins, such as the edema factor (EF) of Bacillus anthracis. These disturb cellular metabolism by catalyzing production of excessive amounts of the regulatory molecule cAMP. Here, a structure-based method, where a 3D-pharmacophore that fit the active site of EF was constructed from fragments, was used to identify non-nucleotide inhibitors of EF. A library of small molecule fragments was docked to the EF-active site in existing crystal structures, and those with the highest HINT scores were assembled into a 3D-pharmacophore. About 10,000 compounds, from over 2.7 million compounds in the ZINC database, had a similar molecular framework. These were ranked according to their docking scores, using methodology that was shown to achieve maximum accuracy (i.e., how well the docked position matched the experimentally determined site for ATP analogues in crystal structures of the complex). Finally, 19 diverse compounds with the best AutoDock binding/docking scores were assayed in a cell-based assay for their ability to reduce cAMP secretion induced by EF. Four of the test compounds, from different structural groups, inhibited in the low micromolar range. One of these has a core structure common to phosphatase inhibitors previously identified by high-throughput assays of a diversity library. Thus, the fragment-based pharmacophore identified a small number of diverse compounds for assay, and greatly enhanced the selection process of advanced lead compounds for combinatorial design.     

Searching the conformational complexity and binding properties of HDAC6 through docking and molecular dynamic simulations

Histone deacetylases (HDACs) are a family of proteins involved in the deacetylation of histones and other non-histones substrates. HDAC6 belongs to class II and shares similar biological functions with others of its class. Nevertheless, its three-dimensional structure that involves the catalytic site remains unknown for exploring the ligand recognition properties. Therefore, in this contribution, homology modeling, 100-ns-long Molecular Dynamics (MD) simulation and docking calculations were combined to explore the conformational complexity and binding properties of the catalytic domain 2 from HDAC6 (DD2-HDAC6), for which activity and affinity toward five different ligands have been reported. Clustering analysis allowed identifying the most populated conformers present during the MD simulation, which were used as starting models to perform docking calculations with five DD2-HDAC6 inhibitors: Cay10603 (CAY), Rocilinostat (RCT), Tubastatin A (TBA), Tubacin (TBC), and Nexturastat (NXT), and then were also submitted to 100-ns-long MD simulations. Docking calculations revealed that the five inhibitors bind at the DD2-HDAC6 binding site with the lowest binding free energy, the same binding mode is maintained along the 100-ns-long MD simulations. Overall, our results provide structural information about the molecular flexibility of apo and holo DD2-HDAC6 states as well as insight of the map of interactions between DD2-HDAC6 and five well-known DD2-HDAC6 inhibitors allowing structural details to guide the drug design. Finally, we highlight the importance of combining different theoretical approaches to provide suitable structural models for structure-based drug design.     

Discovery of new inhibitors of aldose reductase from molecular docking and database screening

Aldose reductase (ALR2) is a target enzyme for the treatment of diabetic complications. Owing to the limited number of currently available drugs for the treatment of diabetic complications, the discovery of new inhibitors of ALR2 that can potentially be optimized as drugs appears highly desirable. In this study, a molecular docking analysis of the structures of more than 127,000 organic compounds contained in the National Cancer Institute database was performed to find and score molecules that are complementary to ALR2. Besides retrieving several carboxylic acid derivatives, which are known to generally inhibit aldose reductase, docking proposed other families of putative inhibitors such as sulfonic acids, nitro-derivatives, sulfonamides and carbonyl derivatives. Twenty-five compounds, chosen as the highest-scoring representatives of each of these families, were tested as aldose reductase inhibitors. Five of them were found to inhibit aldose reductase in the micromolar range. For these active compounds, selectivity with respect to the closely-related aldehyde reductase was determined by measuring the corresponding inhibitory activities. The structures of the complexes between the new lead inhibitors and aldose reductase, here refined with molecular mechanics and molecular dynamics calculations, suggest that new pharmacophoric groups can bind aldose reductase very efficiently. In the case of the family of the nitro-derivative inhibitors, a class of particularly interesting compounds, a round of optimizations was performed with the synthesis and biological evaluation of a series of derivatives aimed at testing the proposed binding mode and at improving interaction with active site residues. Starting from a hit compound having an IC(50) of 42 microM, the most potent compound synthesized showed a 10-fold increase in inhibitory activity and 10-fold selectivity with respect to ALR1, and structure--activity relationships of the designed compounds were in agreement with the proposed mode of binding at the active site.     

Correlation between biological activity and binding energy in systems of integrin with cyclic RGD-containing binders: a QM/MM molecular dynamics study

We here report a combined quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) study on the binding interactions between the α(V)β(3) integrin and eight cyclic arginine-glycine-aspartate (RGD) containing peptides. The initial conformation of each peptide within the binding site of the integrin was determined by docking the ligand to the reactive site of the integrin crystal structure with the aid of docking software FRED. The subsequent QM/MM MD simulations of the complex structures show that these eight cyclic RGD-peptides have a generally similar interaction mode with the binding site of the integrin to the cyclo(RGDf-N[M]V) analog found in the crystal structure. Still, there are subtle differences in the interactions of peptide ligands with the integrin, which contribute to the different inhibition activities. The averaged QM/MM protein-ligand interaction energy (IE) is remarkably correlated to the biological activity of the ligand. The IE, as well as a three-variable model which is somewhat interpretable, thus can be used to predict the bioactivity of a new ligand quantitatively, at least within a family of analogs. The present study establishes a helpful protocol for advancing lead compounds to potent inhibitors.     

Docking studies on the binding of quinoline derivatives and hematin to Plasmodium falciparum lactate dehydrogenase

The literature has reported that ferriprotoporphyrin IX (hematin) intoxicates the malarial parasite through competition with NADH for the active site of the enzyme lactate dehydrogenase (LDH). In order to avoid this, the parasite polymerizes hematin to hemozoin. The quinoline derivatives are believed to form complexes with dimeric hematin, avoiding the formation of hemozoin and still inhibiting LDH. In order to investigate this hypothesis we calculated the docking energies of NADH and some quinoline derivatives (in the free forms and in complex with dimeric hematin) in the active site of the Plasmodium falciparum LDH (PfLDH). Ours results showed better docking score values to the complexes when compared to the free compounds, pointing them as more efficient inhibitors of Pf_LDH. Further we performed Molecular Dynamics (MD) simulations studies on the best docking conformation of the complex chloroquine-dimeric hematin with PfLDH. Our in silico results corroborate experimental data suggesting a possible action route for the quinoline derivatives in the inhibition of PfLDH.     

Structural insights into the binding mode of flavonols with the active site of matrix metalloproteinase-9 through molecular docking and molecular dynamic simulations studies

Cartilage degradation in rheumatoid arthritis is mediated principally by the collagenases and gelatinases. Gelatinase B (also called matrix metalloproteinase 9 – MMP-9), is a valid target molecule which is known to participate in cartilage degradation as well as angiogenesis associated with the disease and inhibition of its activity shall prevent cartilage damage and angiogenesis. The focus of this study is to investigate the possibilities of MMP-9 inhibition by flavonol class of bioflavonoids by studying their crucial binding interactions at the active site of MMP 9 using molecular docking (Glide XP and QPLD) and further improvisation by post-docking MM-GBSA and molecular dynamic (MD) simulations. The results show that flavonols can convincingly bind to active site of MMP-9 as demonstrated by their stable interactions at the S1′ specificity pocket and favourable binding energies. Gossypin has emerged as a promising candidate with a docking score of ?14.618 kcal/mol, binding energy of ?79.97 kcal/mol and a stable MD pattern over 15 ns. In addition, interaction mechanisms with respect to catalytic site zinc are also discussed. Further, the drug-like characters of the ligands were also analysed using ADME analysis.     

Molecular dynamics simulation study and hybrid pharmacophore model development in human LTA4H inhibitor design

Human leukotriene A4 hydrolase (hLTA4H) is a bi-functional enzyme catalyzes the hydrolase and aminopeptidase functions upon the fatty acid and peptide substrates, respectively, utilizing the same but overlapping binding site. Particularly the hydrolase function of this enzyme catalyzes the rate-limiting step of the leukotriene (LT) cascade that converts the LTA4 to LTB4. This product is a potent pro-inflammatory activator of inflammatory responses and thus blocking this conversion provides a valuable means to design anti-inflammatory agents. Four structurally very similar chemical compounds with highly different inhibitory profile towards the hydrolase function of hLTA4H were selected from the literature. Molecular dynamics (MD) simulations of the complexes of hLTA4H with these inhibitors were performed and the results have provided valuable information explaining the reasons for the differences in their biological activities. Binding mode analysis revealed that the additional thiophene moiety of most active inhibitor helps the pyrrolidine moiety to interact the most important R563 and K565 residues. The hLTA4H complexes with the most active compound and substrate were utilized in the development of hybrid pharmacophore models. These developed pharmacophore models were used in screening chemical databases in order to identify lead candidates to design potent hLTA4H inhibitors. Final evaluation based on molecular docking and electronic parameters has identified three compounds of diverse chemical scaffolds as potential leads to be used in novel and potent hLTA4H inhibitor design.     

Molecular docking and 3D-QSAR studies of Yersinia protein tyrosine phosphatase YopH inhibitors

Three-dimensional quantitative structure-activity relationship (QSAR) studies were conducted on two classes of recently explored compounds with known YopH inhibitory activities. Docking studies were employed to position the inhibitors into the YopH active site to determine the probable binding conformation. Good correlations between the predicated binding free energies and the inhibitory activities were found for two subsets of phosphate mimetics: alpha-ketocarboxylic acid and squaric acid (R2=0.70 and 0.68, respectively). The docking results also provided a reliable conformational alignment scheme for 3D-QSAR modeling. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were performed based on the docking conformations, giving q2 of 0.734 and 0.754 for CoMFA and CoMSIA models, respectively. The 3D-QSAR models were significantly improved after removal of an outlier (q2=0.829 for CoMFA and q2=0.837 for CoMSIA). The predictive ability of the models was validated using a set of compounds that were not included in the training set. Mapping the 3D-QSAR models to the active site of YopH provides new insight into the protein-inhibitor interactions for this enzyme. These results should be applicable to the prediction of the activities of new YopH inhibitors, as well as providing structural implications for designing potent and selective YopH inhibitors as antiplague agents.     

Accounting for ligand-bound metal ions in docking small molecules on adenylyl cyclase toxins

The adenylyl cyclase toxins produced by bacteria (such as the edema factor (EF) of Bacillus anthracis and CyaA of Bordetella pertussis) are important virulence factors in anthrax and whooping cough. Co-crystal structures of these proteins differ in the number and positioning of metal ions in the active site. Metal ions bound only to the ligands in the crystal structures are not included during the docking. To determine what effect these "missing" metals have on docking results, the AutoDock, LigandFit/Cerius2, and FlexX programs were compared for their ability to correctly place substrate analogues and inhibitors into the active sites of the crystal structures of EF, CyaA, and mammalian adenylate cyclase. Protonating the phosphates of substrate analogues improved the accuracy of docking into the active site of CyaA, where the grid did not account for one of the three Mg2+ ions in the crystal structure. The AutoDock ranking (based on docking energies) of a test group of compounds was relatively unaffected by protonation of carboxyl groups. However, the ranking by FlexX-ChemScore varied significantly, especially for docking to CyaA, suggesting that alternate protonation states should be tested when screening compound libraries with this program. When the charges on the bound metal were set correctly, AutoDock was the most reliable program of the three tested with respect to positioning substrate analogues and ranking compounds according to their experimentally determined ability to inhibit EF.