Biophysical insight into the heparin‐peptide interaction and its modulation by a small molecule

The heparin‐protein interaction plays a vital role in numerous physiological and pathological processes. Not only is the binding mechanism of these interactions poorly understood, studies concerning their therapeutic targeting are also limited. Here, we have studied the interaction of the heparin interacting peptide (HIP) from Tat (which plays important role in HIV infections) with heparin. Isothermal titration calorimetry binding exhibits distinct biphasic isotherm with two different affinities in the HIP‐heparin complex formation. Overall, the binding was mainly driven by the nonionic interactions with a small contribution from ionic interactions. The stoichiometric analysis suggested that the minimal site for a single HIP molecule is a chain of 4 to 5 saccharide molecules, also supported by docking studies. The investigation was also focused on exploiting the possibility of using a small molecule as an inhibitor of the HIP‐heparin complex. Quinacrine, because of its ability to mimic the HIP interactions with heparin, was shown to successfully modulate the HIP‐heparin interactions. This result demonstrates the feasibility of inhibiting the disease relevant heparin‐protein interactions by a small molecule, which could be an effective strategy for the development of future therapeutic agents.     

Crystal structure at 2.8 A of Huntingtin-interacting protein 1 (HIP1) coiled-coil domain reveals a charged surface suitable for HIP1 protein interactor (HIPPI)

Huntington's disease is a genetic neurological disorder that is triggered by the dissociation of the huntingtin protein (htt) from its obligate interaction partner Huntingtin-interacting protein 1 (HIP1). The release of the huntingtin protein permits HIP1 protein interactor (HIPPI) to bind to its recognition site on HIP1 to form a HIPPI/HIP1 complex that recruits procaspase-8 to begin the process of apoptosis. The interaction module between HIPPI and HIP1 was predicted to resemble a death-effector domain. Our 2.8-Å crystal structure of the HIP1 371-481 subfragment that includes F432 and K474, which is important for HIPPI binding, is not a death-effector domain but is a partially opened coiled coil. The HIP1 371-481 model reveals a basic surface that we hypothesize to be suitable for binding HIPPI. There is an opened region next to the putative HIPPI site that is highly negatively charged. The acidic residues in this region are highly conserved in HIP1 and a related protein, HIP1R, from different organisms but are not conserved in the yeast homologue of HIP1, sla2p. We have modeled ∼ 85% of the coiled-coil domain by joining our new HIP1 371-481 structure to the HIP1 482-586 model (Protein Data Bank code: 2NO2). Finally, the middle of this coiled-coil domain may be intrinsically flexible and suggests a new interaction model where HIPPI binds to a U-shaped HIP1 molecule.     

Tumor suppressor PTEN affects tau phosphorylation: deficiency in the phosphatase activity of PTEN increases aggregation of an FTDP-17 mutant Tau

Huntington disease (HD) is caused by expansion of a polyglutamine (polyQ) domain in the protein known as huntingtin (htt), and the disease is characterized by selective neurodegeneration. Expansion of the polyQ domain is not exclusive to HD, but occurs in eight other inherited neurodegenerative disorders that show distinct neuropathology. Yet in spite of the clear genetic defects and associated neurodegeneration seen with all the polyQ diseases, their pathogenesis remains elusive. The present review focuses on HD, outlining the effects of mutant htt in the nucleus and neuronal processes as well as the role of cell-cell interactions in HD pathology. The widespread expression and localization of mutant htt and its interactions with a variety of proteins suggest that mutant htt engages multiple pathogenic pathways. Understanding these pathways will help us to elucidate the pathogenesis of HD and to target therapies effectively.     

Integrative computational protocol for the discovery of inhibitors of the <Emphasis Type="Italic">Helicobacter pylori</Emphasis> nickel response regulator (NikR)

In order to identify novel inhibitors of the Helicobacter pylori nickel response regulator (HpNikR) an integrative protocol was performed for half a million compounds retrieved from the ZINC database. We firstly implement a structure-based virtual screening to build a library of potential inhibitors against the HpNikR using a docking analysis (AutoDock Vina). The library was then used to perform a hierarchical clustering of docking poses, based on protein-contact footprints calculation from the multiple conformations given by the AutoDock Vina software, and the drug-protein interaction analyses to identify and remove potential promiscuous compounds likely interacting with human proteins, hence causing drug side effects. 250 drug-like compounds were finally proposed as non-promicuous potential inhibitors for HpNikR. These compounds target the DNA-binding sites of HpNikR so that HpNikR-compound binding could be able to mimic key interactions in the DNA-protein recognition process. HpNikR inhibitors with promising potential against H. pylori could also act against other human bacterial pathogens due to the conservation of targeting motif of NikR involved in DNA-protein interaction.     

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.     

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.     

Huntingtin-interacting protein HIP14 is a palmitoyl transferase involved in palmitoylation and trafficking of multiple neuronal proteins

In neurons, posttranslational modification by palmitate regulates the trafficking and function of signaling molecules, neurotransmitter receptors, and associated synaptic scaffolding proteins. However, the enzymatic machinery involved in protein palmitoylation has remained elusive. Here, using biochemical assays, we show that huntingtin (htt) interacting protein, HIP14, is a neuronal palmitoyl transferase (PAT). HIP14 shows remarkable substrate specificity for neuronal proteins, including SNAP-25, PSD-95, GAD65, synaptotagmin I, and htt. Conversely, HIP14 is catalytically invariant toward paralemmin and synaptotagmin VII. Exogenous HIP14 enhances palmitoylation-dependent vesicular trafficking of several acylated proteins in both heterologous cells and neurons. Moreover, interference with endogenous expression of HIP14 reduces clustering of PSD-95 and GAD65 in neurons. These findings define HIP14 as a mammalian palmitoyl transferase involved in the palmitoylation and trafficking of multiple neuronal proteins.     

Rational design of novel diketoacid-containing ferrocene inhibitors of HIV-1 integrase

Molecular interaction field, density functional, and docking studies of novel potential ferrocene inhibitors of HIV-1 integrase (IN) are reported. The high docking scores, analysis of the ligand-receptor interactions in the active site as well as the molecular interaction potential calculations at the binding site of the receptor indicate important features for novel HIV-1 IN inhibitors. We also confirm in this work a novel binding trench in HIV-1 integrase, recently reported in a theoretical work by other authors. This observation may be interesting since the lack of detailed structural information about IN-ligand interactions has hampered the design of IN inhibitors. Our proposed ligands are open to experimental synthesis and testing.     

Drug screening of α-amylase inhibitors as candidates for treating diabetes

In the present study, the identification of potential α-amylase inhibitors is explored as a potential strategy for treating type-2 diabetes mellitus. A computationally driven approach using molecular docking was employed to search for new α-amylase inhibitors. The interactions of potential drugs with the enzyme's active site were investigated and compared with the contacts established by acarbose (a reference drug for α-amylase inhibition) in the crystallographic structure 1B2Y. For this active site characterization, both molecular docking and molecular dynamics simulations were performed, and the residues involved in the α-amylase–acarbose complex were considered to analyse the potential drug's interaction with the enzyme. Two potential α-amylase inhibitors (AN-153I105594 and AN-153I104845) have been selected following this computational strategy. Both compounds established a large number of interactions with key binding site α-amylase amino acids and obtained a comparable docking score concerning the reference drug (acarbose). Aiming to further analyse candidates' properties, their ADME (absorption, distribution, metabolism, excretion) parameters, druglikeness, organ toxicity, toxicological endpoints and median lethal dose (LD₅₀) were estimated. Overall estimations are promising for both candidates, and in silico toxicity predictions suggest that a low toxicity should be expected.     

Docking and molecular dynamics study on the inhibitory activity of <Emphasis Type="Italic">N,N</Emphasis>-disubstituted-trifluoro-3-amino-2-propanols-based inhibitors of cholesteryl ester transfer protein

Extensive studies suggest direct links between cholesteryl ester transfer protein (CETP), high-density lipoproteins-cholesterol level and cardiovascular diseases. Many therapeutic approaches are aimed at the CETP. A series of N, N-disubstituted-trifluoro-3-amino-2-propanol analogues are among the most highly potent and selective inhibitors of CETP described to date. For in-depth investigation into the structural and chemical features responsible for exploring the binding pocket of these compounds, as well as for the binding recognition mechanism concerned, we performed a series of automated molecular docking operations. Moreover, the docking results were quite robust as further validated by molecular dynamics. The docking results reveal that the binding site mainly consists of two hydrophobic regions (P1 and P2 site) which are able to accommodate the lipophilic arms of the compounds investigated. Val421 in P1 site and Met194 in P2 site could be considered to be two important residues in forming the two hydrophobic regions. The presence of residues Phe197 and Phe463 in P2 site may be responsible for the binding recognition through π-π stacking interactions. The hydrophobic 3-phenoxy substituent may be important in creating the preferable inhibitive capability for increasing the binding potency. The hydrophobic character of the tetrafluoroethoxybenzyl group at position 3 displays better hydrophobicity than a shorter hydrophobic substituent. An interaction model of CETP-inhibitors is derived that can be successfully used to explain the different biologic activities of these inhibitors. It is anticipated that the findings reported here may provide very useful information or clues for designing effective drugs for the therapeutic treatment of CETP-related cardiovascular diseases.     

Design of new chemotherapeutics against the deadly anthrax disease. Docking and molecular dynamics studies of inhibitors containing pyrrolidine and riboamidrazone rings on nucleoside hydrolase from Bacillus anthracis

Anthrax is a disease caused by Bacillus anthracis, a dangerous biological warfare agent already used for both military and terrorist purposes. An important selective target for chemotherapy against this disease is nucleoside hydrolase (NH), an enzyme still not found in mammals. Having this in mind we have performed molecular docking studies, aiming to analyze the three-dimensional positioning of six known inhibitors of Trypanosoma vivax NH (TvNH) in the active site of B. anthracis NH (BaNH). We also analyzed the main interactions of these compounds with the active site residues of BaNH and the relevant factors to biological activity. These results, together with further molecular dynamics (MD) simulations, pointed out to the most promising compounds as lead for the design of potential inhibitors of BaNH. Most of the docking and MD results obtained corroborated to each other. Additionally, the docking results also suggested a good correlation with experimental data.     

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.     

Expression and characterization of full-length human huntingtin, an elongated HEAT repeat protein

Huntington disease is an inherited neurodegenerative disorder that is caused by expanded CAG trinucleotide repeats, resulting in a polyglutamine stretch of >37 on the N terminus of the protein huntingtin (htt). htt is a large (347 kDa), ubiquitously expressed protein. The precise functions of htt are not clear, but its importance is underscored by the embryonic lethal phenotype in htt knock-out mice. Despite the fact that the htt gene was cloned 13 years ago, little is known about the properties of the full-length protein. Here we report the expression and preliminary characterization of recombinant full-length wild-type human htt. Our results support a model of htt composed entirely of HEAT repeats that stack to form an elongated superhelix.     

Phytosterols and triterpenes from Morinda lucida Benth (Rubiaceae) as potential inhibitors of anti-apoptotic BCL-XL,BCL-2, and MCL-1: an in-silico study

Deregulation of the normal cellular apoptotic function is a fundamental element in the etiology of most cancers and the anti-apoptotic B cell lymphoma 2 (BCL?2) protein family is known to play crucial role in the regulation of this function. Overexpression of this protein family has been implicated in some cancers, such that agents that could inhibit their over-activity are now being explored for anticancer drug development. A number of studies have revealed the anticancer potential of Morinda lucida-derived extracts and compounds. In search of more inhibitors of this anti-apoptotic protein family from plant resources, 47 compounds, identified in Morinda lucida Benth (Rubiaceae) were screened for their inhibitory activities against BCL-XL, BCL-2, and MCL-1 by molecular docking using BINDSURF, while binding interactions of the top compounds were viewed with PyMOL. Druglikeness and Absorption–Distribution–Metabolism–Excretion (ADME) parameters of the top 6 compounds from docking study were evaluated using SuperPred webserver. Results revealed that out of the 47 compounds, 2 triterpenes (ursolic acid and oleanolic acid) and 4 phytosterols (cycloartenol, campesterol, stigmasterol, and β-sitosterol) have higher binding affinities for the selected BCL-2 proteins, compared to known standard inhibitors; these compounds also fulfill oral drugability of Lipinski rule of five. Therefore, since these Morinda lucida-derived phytosterols and triterpenes show high binding affinity toward the selected anti-apoptotic proteins and exhibited good drugability characteristics, they qualify for further study on drug development against cancers characterized by overexpression of this family of protein.     

Discovery of novel SecA inhibitors of Candidatus Liberibacter asiaticus by structure based design

Candidatus Liberibacter asiaticus is the causal agent of Huanglongbing (HLB) disease of citrus. Current management practices have not been able to control HLB and stop the spread of HLB. The current study is focused on screening small molecule inhibitors against SecA protein of Ca. L. asiaticus. Homology modeling, structure based virtual screening and molecular docking methods have been used to find the novel inhibitory compounds against SecA activity at ATP binding region. At 20 μm 17 compounds showed >50% inhibition and four compounds had more than 65% inhibition. The most active compound has IC₅₀ value of 2.5 μM. The differences between the activities of the compounds are explained by their inter-molecular interactions at ATP binding site.     

Structure based de novo design of novel glycogen synthase kinase 3 inhibitors

Glycogen synthase kinase 3 (GSK-3) is a serine/threonine kinase that has captured great attention in drug discovery projects. Structure based design has been successfully carried out to find a novel class of GSK-3 inhibitors using the Ludi de novo ligand design program. A total of 15 potential leads are suggested from the study. The structures have been validated through detailed analysis of the Ludi score values and by molecular docking experiment using FlexX. The hits have been further verified through: (1) visual examination of how well the hits dock into the GSK-3beta binding site; (2) comparative analysis of their FlexX, G_Score, PMF_Score, ChemScore, and D_scores values; (3) a comparative investigation of the docking scores of the hits with those of the reported inhibitors after calibration of the docking procedure with 17 previously reported inhibitors; (4) determination of the binding mode of the hits and comparison with that of the so far known inhibitors. Hits retaining interactions with the common amino acids of GSK-3beta binding site were taken to represent potential leads. Structurally the hits designed are mainly flat nitrogen heterocycles. These hits are expected to be important additions to the search of GSK-3 inhibitors and may provide invaluable insights to further understand the structural basis of catalysis and inhibition of this kinase.     

Computational Insights into the Inhibitory Mechanism of Human AKT1 by an Orally Active Inhibitor,MK-2206

The AKT signaling pathway has been identified as an important target for cancer therapy. Among small-molecule inhibitors of AKT that have shown tremendous potential in inhibiting cancer, MK-2206 is a highly potent, selective and orally active allosteric inhibitor. Promising preclinical anticancer results have led to entry of MK-2206 into Phase I/II clinical trials. Despite such importance, the exact binding mechanism and the molecular interactions of MK-2206 with human AKT are not available. The current study investigated the exact binding mode and the molecular interactions of MK-2206 with human AKT isoforms using molecular docking and (un)binding simulation analyses. The study also involved the docking analyses of the structural analogs of MK-2206 to AKT1 and proposed one as better inhibitor. The Dock was used for docking simulations of MK-2206 into the allosteric site of AKT isoforms. The Ligplot+ was used for analyses of polar and hydrophobic interactions between AKT isoforms and the ligands. The MoMa-LigPath web server was used to simulate the ligand (un)binding from the binding site to the surface of the protein. In the docking and (un)binding simulation analyses of MK-2206 with human AKT1, the Trp-80 was the key residue and showed highest decrease in the solvent accessibility, highest number of hydrophobic interactions, and the most consistent involvement in all (un)binding simulation phases. The number of molecular interactions identified and calculated binding energies and dissociation constants from the co-complex structures of these isoforms, clearly explained the varying affinity of MK-2206 towards these isoforms. The (un)binding simulation analyses identified various additional residues which despite being away from the binding site, play important role in initial binding of the ligand. Thus, the docking and (un)binding simulation analyses of MK-2206 with AKT isoforms and its structure analogs will provide a suitable model for studying drug-protein interaction and will help in designing better drugs.     

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.     

Discovery of acetylcholinesterase peripheral anionic site ligands through computational refinement of a directed library

The formation of beta-amyloid plaques in the brain is a key neurodegenerative event in Alzheimer's disease. Small molecules capable of binding to the peripheral anionic site of acetylcholinesterase (AChE) have been shown to inhibit the AChE-induced aggregation of the beta-amyloid peptide. Using the combination of a computational docking model and experimental screening, five compounds that completely blocked the amyloidogenic effect of AChE were rapidly identified from an approximately 200-member library of compounds designed to disrupt protein-protein interactions. Critical to this docking model was the inclusion of two explicit water molecules that are tightly bound to the enzyme. Interestingly, none of the tested compounds inhibited the related enzyme butyrylcholinesterase (BuChE) up to their aqueous solubility limits. These compounds are among the most potent inhibitors of amyloid beta-peptide aggregation and are equivalent only to propidium, a well-characterized AChE peripheral anionic site binder and aggregation inhibitor.