Discovery and evaluation of potential sonic hedgehog signaling pathway inhibitors using pharmacophore modeling and molecular dynamics simulations

Sonic hedgehog (Shh) plays an important role in the activation of Shh signaling pathway that regulates preservation and rebirth of adult tissues. An abnormal activation of this pathway has been identified in hyperplasia and various tumorigenesis. Hence the inhibition of this pathway using a Shh inhibitor might be an efficient way to treat a wide range of malignancies. This study was done in order to develop a lead chemical candidate that has an inhibitory function in the Shh signaling pathway. We have generated common feature pharmacophore models using three-dimensional (3D) structural information of robotnikinin, an inhibitor of the Shh signaling pathway, and its analogs. These models have been validated with fit values of robotnikinin and its analogs, and the best model was used as a 3D structural query to screen chemical databases. The hit compounds resulted from the screening docked into a proposed binding site of the Shh named pseudo-active site. Molecular dynamics (MD) simulations were performed to investigate detailed binding modes and molecular interactions between the hit compounds and functional residues of the pseudo-active site. The results of the MD simulation analyses revealed that the hit compounds can bind the pseudo-active site with high affinity than robotnikinin. As a result of this study, a candidate inhibitor (GK 03795) was selected as a potential lead to be employed in future Shh inhibitor design.     

Theoretical studies on the interaction between the nitrile-based inhibitors and the catalytic triad of Cathepsin K

Computational studies on the interaction of novel inhibitor compounds with the Cathepsin K protease have been performed to study the inhibition properties of the inhibitor compounds. The quantum chemical calculations have been performed to analyze the molecular geometries, structural stability, reactivity, nature of interaction, and the charge transfer properties using B3LYP level of theory by implementing 6-311g(d,p) basis set. The calculated C–S and N–H…N bond lengths of the inhibitor-triad complexes are found to agree well with the previous literature results. The chemical reactivity of the inhibitors and catalytic triad are analyzed through frontier molecular orbital analysis and found that the inhibitors are subjected to nucleophilic attack by the catalytic triad. The nature of inhibition of the inhibitor compounds is examined using the quantum theory of Atoms in Molecules analysis and found to be partially covalent. The NBO stabilization energy for the Cys – inhibitor are found to be most stable than the other interactions. The molecular dynamic simulations were performed to study the influence of dynamic of the active site on the QM results. The many body decomposition interaction energy calculated for the final results of MD simulation reveals that the dynamic of the active site induces significant changes in the interaction energy and occupancy of H-bonds plays a major role in the stabilizing the active site inhibitor interactions. The present study reveals that the inhibitor compounds can inhibit the proteolytic activity of the proteases on binding with the catalytic active site.     

Small molecule inhibitors of IκB kinase β: A chip-based screening and molecular docking simulation

A chip-based screening system for IκB kinase β (IKKβ) has been developed by physically immobilizing the substrate IκBα on a glass matrix using a calixarene linker. Phosphorylation of IκBα by IKKβ and ATP was quantitated using a fluorescently labeled antibody. Using this efficient assay system a chemical library of 2000 bioactive compounds was screened against IKKβ and four were identified as good inhibitors, namely, aurintricarboxylic acid, diosmin, ellagic acid, and hematein. None of them have been reported to be an inhibitor of IKKβ although they were implicated in various NFκB-mediated biological processes. Our enzyme-based assay showed that IC50 of the four inhibitors is comparable with that of IKK-16, a previously known strong inhibitor. Molecular docking simulation shows that the hydrophobic moiety of an inhibitor interacts with the four hydrophobic residues (Leu21, Val29, Val152, and Ile165) of the active site. The MM-PBSA calculation suggests that these hydrophobic interactions appear to be the predominant contributor to the binding free energy. As IKKβ is ubiquitously expressed in various cell types and executes many biological functions, the enzyme and cell specificity of the four inhibitors need to be rigorously tested before accepted as a drug candidate.     

Targeting Outer Membrane Protein Component AdeC for the Discovery of Efflux Pump Inhibitor against AdeABC Efflux Pump of Multidrug Resistant <Emphasis Type="BoldItalic">Acinetobacter baumannii</Emphasis>

The structure and functioning of multidrug efflux systems provide us with a better understanding of the transport of various antibiotics, thus giving a path for the discovery of effective compounds for combating the multidrug resistance in Acinetobacter baumannii. In the present study, a number of computational techniques have been used to search for an inhibitor for the RND efflux pump, AdeABC, of A. baumannii targeting specifically its outermost component, i.e., AdeC. We have prepared the three-dimensional structure for AdeC using MODELLER v9.16 and identified its active binding site using SiteMap. Using high-throughput virtual screening, we identified compounds from a large library of biogenic compounds on the basis of their effective interaction at the binding site of AdeC. The validation of docking step was performed by plotting ROC curve (enrichment calculations). The docked complexes were further analyzed for their binding free energies by molecular mechanics using Generalized Born model and Solvent Accessibility (MMGBSA). The molecular dynamics simulation was performed for AdeC-ZINC77257599 complex using GROMACS. The present rational drug designing, molecular mechanics and molecular dynamics data provided an inhibitor, i.e, ZINC77257599 [(3R,4Z,6E,8E)-3-hydroxy-2,2,4-trimethyl-10-oxazol-5-yl-deca-4,6,8-trienamide], for the outer membrane protein component (AdeC) of efflux pump AdeABC of A. baumannii.     

Discovery of promising FtsZ inhibitors by E-pharmacophore, 3D-QSAR,molecular docking study,and molecular dynamics simulation

Filamentous temperature-sensitive protein Z (FtsZ), playing a key role in bacterial cell division, is regarded as a promising target for the design of antimicrobial agent. This study is looking for potential high-efficiency FtsZ inhibitors. Ligand-based pharmacophore and E-pharmacophore, virtual screening and molecular docking were used to detect promising FtsZ inhibitors, and molecular dynamics simulation was used to study the stability of protein-ligand complexes in this paper. Sixty-three inhibitors from published literatures with pIC₅₀ ranging from 2.483 to 5.678 were collected to develop ligand-based pharmacophore model. 4DXD bound with 9PC was selected to develop the E-pharmacophore model. The pharmacophore models validated by test set method and decoy set were employed for virtual screening to exclude inactive compounds against ZINC database. After molecular docking, ADME analysis, IFD docking and MM-GBSA, 8 hits were identified as potent FtsZ inhibitors. A 50?ns molecular dynamics simulation was implemented on the compounds to assess the stability between potent inhibitors and FtsZ. The results indicated that the candidate compounds had a high docking score and were strongly combined with FtsZ by forming hydrogen bonding interactions with key amino acid residues, and van der Waals forces and hydrophobic interactions had significant contribution to the stability of the binding. Molecular dynamics simulation results showed that the protein-ligand compounds performed well in both the stability and flexibility of the simulation process.     

Computational insights into the binding pattern of mitochondrial calcium uniporter inhibitor through homology modeling,molecular dynamics simulation,binding free energy prediction and density functional theory calculation

Abstract

The mitochondrial calcium uniporter (MCU) is the critical protein of the inner mitochondrial membrane that is the primary mediator for calcium uptake into the mitochondrial matrix. Herein we built the optimal homology model of human MCU which was refined through all-atom molecular dynamics simulation. Then, the binding mode of known inhibitor was predicted through molecular docking method, along with molecular dynamics simulation and binding free energy calculation to verify the docking result and stability of the protein-inhibitor complex. Finally, density functional theory (DFT) calculation enhanced our understanding of the molecular interaction of MCU inhibitor. Our research would provide a deeper insight into the interactions between human MCU and its inhibitor, which boosts to develop novel therapy against MCU related disease.

Communicated by Ramaswamy H. Sarma     

Intermolecular enzyme-ligand animation in the active site of porcine pancreatic elastase with acetyl-alanine-proline-alanine by means of molecular dynamics calculation

An interactive molecular display system has been developed to enable the visualization of the 540 conformations derived from a molecular dynamics calculation of the fluctuations of the enzyme-ligand complex formed between porcine pancreatic elastase (PPE) and acetyl-alanine-proline-alanine (APA). Dynamical interactions between the receptor and the inhibitor are observed at the active site, e.g. the pyrrolidone ring of the ligand 2-proline residue is observed to flex via restrained dihedral angle rotations; the terminal acetate moiety is seen to move between two adjacent binding loci. An animated molecular graphics display, linked to a molecular or stochastic dynamics method, is an instructive and predictive tool for investigating dynamical interactions of enzyme-ligand binding.     

Structure-based discovery of a novel non-peptidic small molecular inhibitor of caspase-3

Ac-DNLD-CHO is a novel caspase-3 specific peptide inhibitor that was rationally designed by our computational strategy. The specificity was shown to be due to the specific interaction of NLD moiety with the active site of caspase-3 on the basis of docking mode and site-directed mutagenesis analyses. Here, we computationally screened non-peptidic small molecular inhibitors of caspase-3 from our chemical library using a reliable pharmacophore derived from the specific binding mode of NLD. Through in vitro enzyme assay of the screened candidate compounds, we discovered a novel caspase-3 specific small molecular inhibitor, CS4566, which has a unique scaffold structure. The binding mode of CS4566 to caspase-3 mimics that of NLD, especially LD moiety. This represents a promising lead compound for creating non-peptidic pharmaceuticals for caspase-mediated diseases, such as neurodegenerative disorders.     

An explorative study on Staphylococcus aureus MurE inhibitor: induced fit docking,binding free energy calculation,and molecular dynamics

Staphylococcus aureus MurE enzyme catalyzes the addition of l-lysine as third residue of the peptidoglycan peptide moiety. Due to the high substrate specificity and its ubiquitous nature among bacteria, MurE enzyme is considered as one of the potential target for the development of new therapeutic agents. In the present work, induced fit docking (IFD), binding free energy calculation, and molecular dynamics (MD) simulation were carried out to elucidate the inhibition potential of 2-thioxothiazolidin-4-one based inhibitor 1 against S. aureus MurE enzyme. The inhibitor 1 formed majority of hydrogen bonds with the central domain residues Asn151, Thr152, Ser180, Arg187, and Lys219. Binding free-energy calculation by MM-GBSA approach showed that van der Waals (ΔG_vdW, ?57.30?kcal/mol) and electrostatic solvation (ΔG_solv, ?36.86?kcal/mol) energy terms are major contributors for the inhibitor binding. Further, 30-ns MD simulation was performed to validate the stability of ligand–protein complex and also to get structural insight into mode of binding. Based on the IFD and MD simulation results, we designed four new compounds D1–D4 with promising binding affinity for the S. aureus MurE enzyme. The designed compounds were subjected to the extra-precision docking and binding free energy was calculated for complexes. Further, a 30-ns MD simulation was performed for D1/4C13 complex.     

On the inhibition of HIV-1 protease by hydrazino-ureas displaying the N→CO interaction

To create novel HIV-1 protease (HIV PR) inhibitors, we have extended our investigations of the N→CO interaction as a moiety that reproduces electrostatic properties of the transition state of peptidolysis. Consequently, we prepared a series of compounds with an unusual hydrazino-urea core. In polar protic media, these adopt solely a cyclic constitution displaying the interaction on one side of the molecule while offering a urea moiety on the opposite side meant to hydrogen-bond with the enzyme flaps. Each inhibitor candidate was obtained via a key series of three synthetic steps employing carbonyl-di-imidazole (CDI). It was thus possible to efficiently fuse two independent building blocks, a hydrazine and a protected aminoaldehyde in a convergent manner. NMR and UV analysis proved that all compounds, when dissolved in polar protic media, existed exclusively in the cyclic constitution exhibiting the N→CO interaction. In total, five inhibitor candidates were tested with HIV PR for their potency. The one carrying the least bulk in peripheral substituents showed the highest activity. Its very low molecular weight (365 g/mol) holds great promise for future improvements in affinity without violating Lipinski’s rule of remaining within the limit of 500 g/mol.     

Molecular modeling based approach, synthesis and in vitro assay to new indole inhibitors of hepatitis C NS3/4A serine protease

In an attempt to identify potential HCV NS3 protease inhibitors lead compounds, a series of novel indoles (10a-g) was designed. Molecular modeling study, including fitting to a 3D-pharmacophore model of the designed molecules (10a-g), with HCV NS3 protease hypothesis using catalyst program was fulfilled. Also, the molecular docking into the NS3 active site was examined using Discovery Studio 2.5 software. Several compounds showed significant high simulation docking score and fit values. The designed compounds with high docking score and fit values were synthesized and biologically evaluated in vitro using an NS3 protease binding assay. It appears that most of the tested compounds reveal promising inhibitory activity against NS3 protease. Of these, compounds 10a and 10b demonstrated potent HCV NS3 protease inhibitors with IC₅₀ values of 9 and 12 ??g/mL, respectively. The experimental serine protease inhibitor activities of compounds 10a-g were consistent with their molecular modeling results. Inhibitors from this class have promising characteristics for further development as anti-HCV agents.     

Protein surface roughness accounts for binding free energy of Plasmepsin II‐ligand complexes

The calculation of absolute binding affinities for protein‐inhibitor complexes remains as one of the main challenges in computational structure‐based ligand design. The present work explored the calculations of surface fractal dimension (as a measure of surface roughness) and the relationship with experimental binding free energies of Plasmepsin II complexes. Plasmepsin II is an attractive target for novel therapeutic compounds to treat malaria. However, the structural flexibility of this enzyme is a drawback when searching for specific inhibitors. Concerning that, we performed separate explicitly solvated molecular dynamics simulations using the available high‐resolution crystal structures of different Plasmepsin II complexes. Molecular dynamics simulations allowed a better approximation to systems dynamics and, therefore, a more reliable estimation of surface roughness. This constitutes a novel approximation in order to obtain more realistic values of fractal dimension, because previous works considered only x‐ray structures. Binding site fractal dimension was calculated considering the ensemble of structures generated at different simulation times. A linear relationship between binding site fractal dimension and experimental binding free energies of the complexes was observed within 20 ns. Previous studies of the subject did not uncover this relationship. Regression model, coined FD model, was built to estimate binding free energies from binding site fractal dimension values. Leave‐one‐out cross‐validation showed that our model reproduced accurately the absolute binding free energies for our training set (R² = 0.76; <|error|> =0.55 kcal/mol; SD_error = 0.19 kcal/mol). The fact that such a simple model may be applied raises some questions that are addressed in the article.     

Automatic identification and representation of protein binding sites for molecular docking.

Molecular docking is a popular way to screen for novel drug compounds. The method involves aligning small molecules to a protein structure and estimating their binding affinity. To do this rapidly for tens of thousands of molecules requires an effective representation of the binding region of the target protein. This paper presents an algorithm for representing a protein's binding site in a way that is specifically suited to molecular docking applications. Initially the protein's surface is coated with a collection of molecular fragments that could potentially interact with the protein. Each fragment, or probe, serves as a potential alignment point for atoms in a ligand, and is scored to represent that probe's affinity for the protein. Probes are then clustered by accumulating their affinities, where high affinity clusters are identified as being the "stickiest" portions of the protein surface. The stickiest cluster is used as a computational binding "pocket" for docking. This method of site identification was tested on a number of ligand-protein complexes; in each case the pocket constructed by the algorithm coincided with the known ligand binding site. Successful docking experiments demonstrated the effectiveness of the probe representation.     

Integrating docking and molecular dynamics approaches for a series of proline-based 2,5-diketopiperazines as novel αβ-tubulin inhibitors

In this work, docking tools were utilized in order to study the binding properties of more than five hundred of proline-based 2,5-diketopiperazine in the binding site of αβ-tubulin. Results revealed that 20 compounds among them showed lower binding energies in comparison with Tryprostatin-A, a well known tubulin inhibitor and therefore could be potential inhibitors of tubulin. However, the precise evaluation of binding poses represents the similar binding modes for all of these compounds and Tryprostatin-A. Finally, the best docked complex was subjected to a 25 ns molecular dynamics simulation to further validate the proposed binding mode of this compound.     

Study on the binding mode of a pyrrolotriazin derivative with JAK2 by docking and MD simulation

The pyrrolotriazin derivative 2-(4-(4-((7-(3-(N-methylmethylsulfonamido)phenyl)pyrrolo [2,1-f][1,2,4]triazin-2-yl)amino)phenyl)piperidin-1-yl)acetamide (PPA) is a potential Janus kinase 2 (JAK2) inhibitor. The binding mode between PPA and JAK2 was investigated by using a combined method of docking, molecular dynamics (MD) simulation and binding free-energy calculation. The docking calculations preliminarily indicated that there were two possible binding modes 1 and 2; MD simulations and binding free-energy calculations identified that binding mode 1 was more stable and favourable, with the lower MM-PBSA binding free energy of ?34.00?±?0.17?kcal/mol. Moreover, some valuable binding information is revealed as follows: the inhibitor PPA is suitably located at the ATP-binding site of JAK2 and the hydrophobic interaction plays an essential role. PPA not only interacts with residues Leu855, Val863, Ala880, Tyr931, Leu932 and Leu983 via hydrophobic interaction but also interacts with Ser936 and Asp994 by hydrogen bonds. These two factors are advantageous for PPA to strongly bind to JAK2. These results help to understand the action mechanisms and designing new compounds with a higher affinity to JAK2.     

Molecular insights on analogs of HIV PR inhibitors toward HTLV‐1 PR through QM/MM interactions and molecular dynamics studies: comparative structure analysis of wild and mutant HTLV‐1 PR

Retroviruses HTLV‐1 and HIV‐1 are the primary causative agents of fatal adult T‐cell leukemia and acquired immune deficiency syndrome (AIDS) disease. Both retroviruses are similar in characteristics mechanism, and it encodes for protease that mainly involved in the viral replication process. On the basis of the therapeutic success of HIV‐1 PR inhibitors, the protease of HTLV‐1 is mainly considered as a potential target for chemotherapy. At the same time, structural similarities in both enzymes that originate HIV PR inhibitors can also be an HTLV‐1 PR inhibitor. But the expectations failed because of rejection of HIV PR inhibitors from the HTLV‐1 PR binding pocket. In this present study, the reason for the HIV PR inhibitor rejection from the HTLV‐1 binding site was identified through sequence analysis and molecular dynamics simulation method. Functional analysis of M37A mutation in HTLV PR clearly shows that the MET37 specificity and screening of potential inhibitors targeting MET37 is performed by using approved 90% similar HIV PR inhibitor compounds. From this approach, we report few compounds with a tendency to accept/donate electron specifically to an important site residue MET37 in HTLV‐1 PR binding pocket. Copyright © 2014 John Wiley & Sons, Ltd.     

Cyclopropane derivatives as potential human serine racemase inhibitors: unveiling novel insights into a difficult target

d-Serine is the co-agonist of NMDA receptors and binds to the so-called glycine site. d-Serine is synthesized by human serine racemase (SR). Over activation of NMDA receptors is involved in many neurodegenerative diseases and, therefore, the inhibition of SR might represent a novel strategy for the treatment of these pathologies. SR is a very difficult target, with only few compounds so far identified exhibiting weak inhibitory activity. This study was aimed at the identification of novel SR inhibitor by mimicking malonic acid, the best-known SR inhibitor, with a cyclopropane scaffold. We developed, synthesized, and tested a series of cyclopropane dicarboxylic acid derivatives, complementing the synthetic effort with molecular docking. We identified few compounds that bind SR in high micromolar range with a lack of significant correlation between experimental and predicted binding affinities. The thorough analysis of the results can be exploited for the development of more potent SR inhibitors.     

Computer-aided rational molecular design of argifin-derivatives with increased inhibitory activity against chitinase B from Serratia marcescens

Argifin, a novel pentapeptide chitinase inhibitor isolated from Gliocladium fungal culture, is a promising candidate for the development of new fungicides, insecticides, and anti-asthma medications. In this study, we undertook rational molecular design of argifin-derivatives and tested them against chitinase B from Serratia marcescens (SmChiB). The work involved molecular dynamics simulation with explicit water molecules, the molecular docking calculation, and free-energy analysis using the molecular mechanics Poisson–Boltzmann surface area method. The custom-designed derivatives were synthesized via effective solid phase synthesis, developed recently in our laboratory, and their inhibitory activities were measured against SmChiB. Finally, we identified and obtained a derivative which exhibited 28-fold more inhibition than argifin itself, a compound in which the d-Ala(5) of argifin was replaced with d-Leu and the 4-benzylpiperdine was attached to l-Asp(4).     

A simple screening method for detecting bindings between oligopeptides and HLA-DR molecules on filter papers: possible application for mapping of putative helper T-cell epitopes on MSP1 of Plasmodium falciparum

Binding capacities of synthetic peptides to HLA-DR molecules were tested on filter papers to identify putative helper T-cell epitopes on a malarial protein. The antigen tested was the merozoite surface glycoprotein 1 (MSP1) of Plasmodium falciparum, a vaccine candidate targeting the asexual erythrocytic stage. Bindings between synthetic oligopeptides and HLA-DR molecules were tested. Such bindings were not non-specific, and a known helper T-cell epitope peptide showed positive binding to the restricting HLA-DR molecule. By using this screening system, we observed the unequal distribution of HLA-DR-binding peptides in 10 out of 17 MSP1 blocks tested. Block #6 of MSP1 seemed to show the highest frequency in the positive binding; on the other hand, blocks #1 and #17, both of which were thought to be vaccine candidate regions, contained fewer HLA-DR binding peptides. This was not inconsistent with the results that block #17 was less stimulatory to peripheral T cells than block #6. The peptides with positive binding to HLA-DR showed actual epitope activities when we tested peptide-driven proliferation of human bulk T-cell lines, and association between the two parameters was statistically significant (P<0.001). For more detailed information for vaccine development, peptides with both IgG- and HLA-DR binding activities were mapped in block #17 of MSP1. Together with these results, we demonstrate that our simple screening system seems to provide essential information for vaccine development through uncovering locations of putative epitopes for human helper T cells.     

Combinatorial Clustering of Residue Position Subsets Predicts Inhibitor Affinity across the Human Kinome

The protein kinases are a large family of enzymes that play fundamental roles in propagating signals within the cell. Because of the high degree of binding site similarity shared among protein kinases, designing drug compounds with high specificity among the kinases has proven difficult. However, computational approaches to comparing the 3-dimensional geometry and physicochemical properties of key binding site residue positions have been shown to be informative of inhibitor selectivity. The Combinatorial Clustering Of Residue Position Subsets (ccorps) method, introduced here, provides a semi-supervised learning approach for identifying structural features that are correlated with a given set of annotation labels. Here, ccorps is applied to the problem of identifying structural features of the kinase atp binding site that are informative of inhibitor binding. ccorps is demonstrated to make perfect or near-perfect predictions for the binding affinity profile of 8 of the 38 kinase inhibitors studied, while only having overall poor predictive ability for 1 of the 38 compounds. Additionally, ccorps is shown to identify shared structural features across phylogenetically diverse groups of kinases that are correlated with binding affinity for particular inhibitors; such instances of structural similarity among phylogenetically diverse kinases are also shown to not be rare among kinases. Finally, these function-specific structural features may serve as potential starting points for the development of highly specific kinase inhibitors.