Bioinformatics can boost metabolomics research.

Metabolomics is the modern term for the field of small molecule research in biology and biochemistry. Currently, metabolomics is undergoing a transition where the classic analytical chemistry is combined with modern cheminformatics and bioinformatics methods, paving the way for large-scale data analysis. We give some background on past developments, highlight current state-of-the-art approaches, and give a perspective on future requirements.

SQM/COSMO Scoring Function at the DFTB3-D3H4 Level: Unique Identification of Native Protein-Ligand Poses.

We have recently introduced the "SQM/COSMO" scoring function which combines a semiempirical quantum mechanical description of noncovalent interactions at the PM6-D3H4X level and the COSMO implicit model of solvation. This approach outperformed standard scoring functions but faced challenges with a metalloprotein featuring a Zn2+···S- interaction. Here, we invoke SCC-DFTB3-D3H4, a higher-level SQM method, and observe improved behavior for the metalloprotein and high-quality results for the other systems. This method holds promise for diverse protein-ligand complexes including metalloproteins.

A Deep Hydrophobic Binding Cavity is the Main Interaction for Different Y2 R Antagonists.

The neuropeptide Y2 receptor (Y2 R) is involved in various pathophysiological processes such as epilepsy, mood disorders, angiogenesis, and tumor growth. Therefore, the Y2 R is an interesting target for drug development. A detailed understanding of the binding pocket could facilitate the development of highly selective antagonists to study the role of Y2 R in vitro and in vivo. In this study, several residues crucial to the interaction of BIIE0246 and SF-11 derivatives with Y2 R were investigated by signal transduction assays. Using the experimental results as constraints, the antagonists were docked into a comparative structural model of the Y2 R. Despite differences in size and structure, all three antagonists display a similar binding site, including a deep hydrophobic cavity formed by transmembrane helices (TM) 4, 5, and 6, as well as a hydrophobic patch at the top of TM2 and 7. Additionally, we suggest that the antagonists block Q3.32 , a position that has been shown to be crucial for binding of the amidated C terminus of NPY and thus for receptor activation.

On-resin Diels-Alder reaction with inverse electron demand: an efficient ligation method for complex peptides with a varying spacer to optimize cell adhesion.

Solid phase peptide synthesis (SPPS) is the method of choice to produce peptides. Several protecting groups enable specific modifications. However, complex peptide conjugates usually require a rather demanding conjugation strategy, which is mostly performed in solution. Herein, an efficient strategy is described using an on-resin Diels-Alder reaction with inverse electron demand (DARinv). This method is compatible with the standard Fmoc/tBu strategy and is easy to monitor. As a proof of concept a titanium binding peptide was modified with a cyclic cell binding peptide (RGD) by DARinv on a solid support applying different tetrazines and alkenes. The generated bulky DARinv linkers were employed to act as the required spacer for RGD mediated cell adhesion on titanium. In vitro studies demonstrated improved cell spreading on DARinv-conjugated peptides and revealed, in combination with molecular dynamics-simulation, new insights into the design of spacers between the RGD peptide and the surface. Performing the DARinv on resin expands the toolbox of SPPS to produce complex peptide conjugates under mild, catalyst free conditions with reduced purification steps. The resulting conjugate can be effectively exploited to promote cell adhesion on biomaterials.

The SQM/COSMO filter: reliable native pose identification based on the quantum-mechanical description of protein-ligand interactions and implicit COSMO solvation.

Current virtual screening tools are fast, but reliable scoring is elusive. Here, we present the 'SQM/COSMO filter', a novel scoring function featuring a quantitative semiempirical quantum mechanical (SQM) description of all types of noncovalent interactions coupled with implicit COSMO solvation. We show unequivocally that it outperforms eight widely used scoring functions. The accuracy and chemical generality of the SQM/COSMO filter make it a perfect tool for late stages of virtual screening.

Engineered bacterial polyester hydrolases efficiently degrade polyethylene terephthalate due to relieved product inhibition.

Recent studies on the enzymatic degradation of synthetic polyesters have shown the potential of polyester hydrolases from thermophilic actinomycetes for modifying or degrading polyethylene terephthalate (PET). TfCut2 from Thermobifida fusca KW3 and LC-cutinase (LCC) isolated from a compost metagenome are remarkably active polyester hydrolases with high sequence and structural similarity. Both enzymes exhibit an exposed active site in a substrate binding groove located at the protein surface. By exchanging selected amino acid residues of TfCut2 involved in substrate binding with those present in LCC, enzyme variants with increased PET hydrolytic activity at 65°C were obtained. The highest activity in hydrolyzing PET films and fibers were detected with the single variant G62A and the double variant G62A/I213S. Both variants caused a weight loss of PET films of more than 42% after 50 h of hydrolysis, corresponding to a 2.7-fold increase compared to the wild type enzyme. Kinetic analysis based on the released PET hydrolysis products confirmed the superior hydrolytic activity of G62A with a fourfold higher hydrolysis rate constant and a 1.5-fold lower substrate binding constant than those of the wild type enzyme. Mono-(2-hydroxyethyl) terephthalate is a strong inhibitor of TfCut2. A determination of the Rosetta binding energy suggested a reduced interaction of G62A with 2PET, a dimer of the PET monomer ethylene terephthalate. Indeed, G62A revealed a 5.5-fold lower binding constant to the inhibitor than the wild type enzyme indicating that its increased PET hydrolysis activity is the result of a relieved product inhibition by mono-(2-hydroxyethyl) terephthalate. Biotechnol. Bioeng. 2016;113: 1658-1665. © 2016 Wiley Periodicals, Inc.

A unique serpin P1' glutamate and a conserved ß-sheet C arginine are key residues for activity, protease recognition and stability of serpinA12 (vaspin).

SerpinA12 (vaspin) is thought to be mainly expressed in adipose tissue and has multiple beneficial effects on metabolic, inflammatory and atherogenic processes related to obesity. KLK7 (kallikrein 7) is the only known protease target of vaspin to date and is inhibited with a moderate inhibition rate. In the crystal structure, the cleavage site (P1-P1') of the vaspin reactive centre loop is fairly rigid compared with the flexible residues before P2, possibly supported by an ionic interaction of P1' glutamate (Glu(379)) with an arginine residue (Arg(302)) of the ß-sheet C. A P1' glutamate seems highly unusual and unfavourable for the protease KLK7. We characterized vaspin mutants to investigate the roles of these two residues in protease inhibition and recognition by vaspin. Reactive centre loop mutations changing the P1' residue or altering the reactive centre loop conformation significantly increased inhibition parameters, whereas removal of the positive charge within ß-sheet C impeded the serpin-protease interaction. Arg(302) is a crucial contact to enable vaspin recognition by KLK7 and it supports moderate inhibition of the serpin despite the presence of the detrimental P1' Glu(379), which clearly represents a major limiting factor for vaspin-inhibitory activity. We also show that the vaspin-inhibition rate for KLK7 can be modestly increased by heparin and demonstrate that vaspin is a heparin-binding serpin. Noteworthily, we observed vaspin as a remarkably thermostable serpin and found that Glu(379) and Arg(302) influence heat-induced polymerization. These structural and functional results reveal the mechanistic basis of how reactive centre loop sequence and exosite interaction in vaspin enable KLK7 recognition and regulate protease inhibition as well as stability of this adipose tissue-derived serpin.

Unwinding of the C-Terminal Residues of Neuropeptide Y is critical for Y2 Receptor Binding and Activation.

Despite recent breakthroughs in the structural characterization of G-protein-coupled receptors (GPCRs), there is only sparse data on how GPCRs recognize larger peptide ligands. NMR spectroscopy, molecular modeling, and double-cycle mutagenesis studies were integrated to obtain a structural model of the peptide hormone neuropeptide Y (NPY) bound to its human G-protein-coupled Y2 receptor (Y2R). Solid-state NMR measurements of specific isotope-labeled NPY in complex with in vitro folded Y2R reconstituted into phospholipid bicelles provided the bioactive structure of the peptide. Guided by solution NMR experiments, it could be shown that the ligand is tethered to the second extracellular loop by hydrophobic contacts. The C-terminal α-helix of NPY, which is formed in a membrane environment in the absence of the receptor, is unwound starting at T(32) to provide optimal contacts in a deep binding pocket within the transmembrane bundle of the Y2R.

Influence of the hinge region and its adjacent domains on binding and signaling patterns of the thyrotropin and follitropin receptor.

Glycoprotein hormone receptors (GPHR) have a large extracellular domain (ECD) divided into the leucine rich repeat (LRR) domain for binding of the glycoprotein hormones and the hinge region (HinR), which connects the LRR domain with the transmembrane domain (TMD). Understanding of the activation mechanism of GPHRs is hindered by the unknown interaction of the ECD with the TMD and the structural changes upon ligand binding responsible for receptor activation. Recently, our group showed that the HinR of the thyrotropin receptor (TSHR) can be replaced by those of the follitropin (FSHR) and lutropin receptor (LHCGR) without effects on surface expression and hTSH signaling. However, differences in binding characteristics for bovine TSH at the various HinRs were obvious. To gain further insights into the interplay between LRR domain, HinR and TMD we generated chimeras between the TSHR and FSHR. Our results obtained by the determination of cell surface expression, ligand binding and G protein activation confirm the similar characteristics of GPHR HinRs but they also demonstrate an involvement of the HinR in ligand selectivity indicated by the observed promiscuity of some chimeras. While the TSHR HinR contributes to specific binding of TSH and its variants, no such contribution is observed for FSH and its analog TR4401 at the HinR of the FSHR. Furthermore, the charge distribution at the poorly characterized LRR domain/HinR transition affected ligand binding and signaling even though this area is not in direct contact with the ligand. In addition our results also demonstrate the importance of the TMD/HinR interface. Especially the combination of the TSHR HinR with the FSHR-TMD resulted in a loss of cell surface expression of the respective chimeras. In conclusion, the HinRs of GPHRs do not only share similar characteristics but also behave as ligand specific structural and functional entities.

Proteolytic activation of prochemerin by kallikrein 7 breaks an ionic linkage and results in C-terminal rearrangement.

The excessive accumulation of adipose tissue in obesity is associated with multiple inflammatory dermatological diseases. Chemerin, a chemoattractant adipokine, dependent on proteolytical activation, is highly expressed in skin. Different proteases have been reported to activate prochemerin, but none is inherently expressed in human skin. In the present study, we identified a tissue-specific protease and investigated the underlying mechanism of activation at the molecular level. We characterized human KLK7 (kallikrein 7) as a prochemerin processing protease in vitro converting prochemerin into active chemerinF(156). The activating truncation by the protease might trigger a structural rearrangement leading to an increased affinity of chemerin to CMKLR1 (chemokine-like receptor 1). Molecular modelling and experimental data suggest an underlying ionic interaction in prochemerin C-terminal domains. These findings provide a general molecular basis for the necessity of C-terminal processing of prochemerin. Moreover, immunohistochemistry was used to investigate prochemerin, KLK7 and the recently identified KLK7 inhibitor vaspin expression in human skin biopsies, and distinct co-localization in psoriatic biopsies was observed. On the basis of these results, it is hypothesized that KLK7 activity may contribute to the development of psoriatic lesions as a consequence of excessive chemerin activation and impaired protease activity regulation by vaspin. Therefore this interaction represents an interesting target for psoriasis therapy and treatment of other obesity-related diseases.

Selective mode of action of guanidine-containing non-peptides at human NPFF receptors.

The binding pocket of both NPFF receptors was investigated, focusing on subtype-selective behavior. By use of four nonpeptidic compounds and the peptide mimetics RF9 and BIBP3226, agonistic and antagonistic properties were characterized. A set of Ala receptor mutants was generated. The binding pocket was narrowed down to the upper part of transmembrane helices V, VI, VII and the extracellular loop 2. Positions 5.27 and 6.59 have been shown to have a strong impact on receptor activation and were suggested to form an acidic, negatively charged binding pocket in both NPFF receptor subtypes. Additionally, position 7.35 was identified to play an important role in functional selectivity. According to docking experiments, the aryl group of AC-216 interacts with position 7.35 in the NPFF(1) but not in the NPFF(2) receptor. These results provide distinct insights into the receptor specific binding pockets, which is necessary for the development of drugs to address the NPFF system.

ParaDockS: a framework for molecular docking with population-based metaheuristics.

Molecular docking is a simulation technique that aims to predict the binding pose between a ligand and a receptor. The resulting multidimensional continuous optimization problem is practically unsolvable in an exact way. One possible approach is the combination of an optimization algorithm and an objective function that describes the interaction. The software ParaDockS is designed to hold different optimization algorithms and objective functions. At the current stage, an adapted particle-swarm optimizer (PSO) is implemented. Available objective functions are (i) the empirical objective function p-Score and (ii) an adapted version of the knowledge-based potential PMF04. We tested the docking accuracy in terms of reproducing known crystal structures from the PDBbind core set. For 73% of the test instances the native binding mode was found with an rmsd below 2 A. The virtual screening efficiency was tested with a subset of 13 targets and the respective ligands and decoys from the directory of useful decoys (DUD). ParaDockS with PMF04 shows a superior early enrichment. The here presented approach can be employed for molecular docking experiments and virtual screenings of large compound libraries in academia as well as in industrial research and development. The performance in terms of accuracy and enrichment is close to the results of commercial software solutions.

Virtual screening and biological characterization of novel histone arginine methyltransferase PRMT1 inhibitors.

Lysine and arginine methyltransferases participate in the posttranslational modification of histones and regulate key cellular functions. Protein arginine methyltransferase 1 (PRMT1) has been identified as an essential component of mixed lineage leukemia (MLL) oncogenic complexes, revealing its potential as a novel therapeutic target in human cancer. The first potent arginine methyltransferase inhibitors were recently discovered by random- and target-based screening approaches. Herein we report virtual and biological screening for novel inhibitors of PRMT1. Structure-based virtual screening (VS) of the Chembridge database composed of 328 000 molecules was performed with a combination of ligand- and target-based in silico approaches. Nine inhibitors were identified from the top-scored docking solutions; these were experimentally tested using human PRMT1 and an antibody-based assay with a time-resolved fluorescence readout. Among several aromatic amines, an aliphatic amine and an amide were also found to be active in the micromolar range.

Thiobarbiturates as sirtuin inhibitors: virtual screening, free-energy calculations, and biological testing.

NAD+-dependent histone deacetylases (sirtuins) are enzymes that cleave acetyl groups from lysine residues in histones and other proteins. Potent selective sirtuin inhibitors are interesting tools for the investigation of the biological functions of these enzymes and may be future drugs for the treatment of cancer or neurodegenerative diseases. Herein we present the results from a protein-based virtual screen of a commercial database with subsequent biological testing of the most promising compounds. The combination of docking and in vitro experimental testing resulted in the identification of novel sirtuin inhibitors with thiobarbiturate structure. To rationalize the experimental results, free-energy calculations were carried out by molecular mechanics Poisson-Boltzmann/surface area (MM-PBSA) calculations. A significant correlation between calculated binding free energies and measured Sirt2 inhibitory activities was observed. The analyses suggested a molecular basis for the interaction of the identified thiobarbiturate derivatives with human Sirt2. Based on the docking and MM-PBSA calculations we synthesized and tested five further thiobarbiturates. The MM-PBSA method correctly predicted the activity of the novel thiobarbiturates. The identified compounds will be used to further explore the therapeutic potential of sirtuin inhibitors.

Structure-activity studies on splitomicin derivatives as sirtuin inhibitors and computational prediction of binding mode.

NAD (+)-dependent histone deacetylases (sirtuins) are enzymes that cleave acetyl groups from lysines in histones and other proteins. Potent selective sirtuin inhibitors are interesting tools for the investigation of the biological functions of those enzymes and may be future drugs for the treatment of cancer. Splitomicin was among the first two inhibitors that were discovered for yeast sirtuins but showed rather weak inhibition on human enzymes. We present detailed structure-activity relationships on splitomicin derivatives and their inhibition of recombinant Sirt2. To rationalize our experimental results, ligand docking followed by molecular mechanics Poisson-Boltzmann/surface area (MM-PBSA) calculations were carried out. These analyses suggested a molecular basis for the interaction of the beta-phenylsplitomicins with human Sirt2. Protein-based virtual screening resulted in the identification of a novel Sirt2 inhibitor chemotype. Selected inhibitors showed antiproliferative properties and tubulin hyperacetylation in MCF7 breast cancer cells and are promising candidates for further optimization as potential anticancer drugs.

Generation of a homology model of the human histamine H(3) receptor for ligand docking and pharmacophore-based screening.

The human histamine H(3) receptor (hH(3)R) is a G-protein coupled receptor (GPCR), which modulates the release of various neurotransmitters in the central and peripheral nervous system and therefore is a potential target in the therapy of numerous diseases. Although ligands addressing this receptor are already known, the discovery of alternative lead structures represents an important goal in drug design. The goal of this work was to study the hH(3)R and its antagonists by means of molecular modelling tools. For this purpose, a strategy was pursued in which a homology model of the hH(3)R based on the crystal structure of bovine rhodopsin was generated and refined by molecular dynamics simulations in a dipalmitoylphosphatidylcholine (DPPC)/water membrane mimic before the resulting binding pocket was used for high-throughput docking using the program GOLD. Alternatively, a pharmacophore-based procedure was carried out where the alleged bioactive conformations of three different potent hH(3)R antagonists were used as templates for the generation of pharmacophore models. A pharmacophore-based screening was then carried out using the program Catalyst. Based upon a database of 418 validated hH(3)R antagonists both strategies could be validated in respect of their performance. Seven hits obtained during this screening procedure were commercially purchased, and experimentally tested in a [(3)H]N(alpha)-methylhistamine binding assay. The compounds tested showed affinities at hH(3)R with K ( i ) values ranging from 0.079 to 6.3 muM.

Structure-activity studies on suramin analogues as inhibitors of NAD+-dependent histone deacetylases (sirtuins).

Suramin is a symmetric polyanionic naphthylurea originally used for the treatment of trypanosomiasis and onchocerciasis. Suramin and diverse analogues exhibit a broad range of biological actions in vitro and in vivo, including, among others, antiproliferative and antiviral activity. Suramin derivatives usually target purinergic binding sites. Class III histone deacetylases (sirtuins) are amidohydrolases that require nicotinamide adenine dinucleotide (NAD(+)) as a cofactor for their catalytic mechanism(.) Deacetylation of the target proteins leads to a change in conformation and alters the activity of the proteins in question. Suramin was reported to inhibit human sirtuin 1 (SIRT1). We tested a diverse set of suramin analogues to elucidate the inhibition of the NAD(+)-dependent histone deacetylases SIRT1 and SIRT2 and discovered selective inhibitors of human sirtuins with potency in the two-digit nanomolar range. In addition, the structural requirements for the binding of suramin derivatives to sirtuins were investigated by molecular docking. The recently published X-ray crystal structure of human SIRT5 in complex with suramin and the human SIRT2 structure were used to analyze the interaction mode of the novel suramin derivatives.

Adenosine mimetics as inhibitors of NAD+-dependent histone deacetylases, from kinase to sirtuin inhibition.

NAD+-dependent histone deacetylases, sirtuins, cleave acetyl groups from lysines of histones and other proteins to regulate their activity. Identification of potent selective inhibitors would help to elucidate sirtuin biology and could lead to useful therapeutic agents. NAD+ has an adenosine moiety that is also present in the kinase cofactor ATP. Kinase inhibitors based upon adenosine mimesis may thus also target NAD+-dependent enzymes. We present a systematic approach using adenosine mimics from one cofactor class (kinase inhibitors) as a viable method to generate new lead structures in another cofactor class (sirtuin inhibitors). Our findings have broad implications for medicinal chemistry and specifically for sirtuin inhibitor design. Our results also raise a question as to whether selectivity profiling for kinase inhibitors should be limited to ATP-dependent targets.