Structure-based discovery of novel cruzain inhibitors with distinct trypanocidal activity profiles.

Over 110 years after the first formal description of Chagas disease, the trypanocidal drugs thus far available have limited efficacy and several side effects. This encourages the search for novel treatments that inhibit T. cruzi targets. One of the most studied anti-T. cruzi targets is the cysteine protease cruzain; it is associated with metacyclogenesis, replication, and invasion of the host cells. We used computational techniques to identify novel molecular scaffolds that act as cruzain inhibitors. First, with a docking-based virtual screening, we identified compound 8, a competitive cruzain inhibitor with a Ki of 4.6 µM. Then, aided by molecular dynamics simulations, cheminformatics, and docking, we identified the analog compound 22 with a Ki of 27 µM. Surprisingly, despite sharing the same isoquinoline scaffold, compound 8 presented higher trypanocidal activity against the epimastigote forms, while compound 22, against the trypomastigotes and amastigotes. Taken together, compounds 8 and 22 represent a promising scaffold for further development of trypanocidal compounds as drug candidates for treating Chagas disease.

Pose, duplicate, then elaborate: Steps towards increased affinity for inhibitors targeting the specificity surface of the Pim-1 kinase.

In this study, fragment-sized hits binding to Pim-1 kinase with initially modest affinity were further optimized by combining computational, synthetic and crystallographic expertise, eventually resulting in potent ligands with affinities in the nanomolar range that address rarely-targeted regions of Pim-1 kinase. Starting from a set of crystallographically validated, chemically distinct fragments that bind to Pim-1 kinase but lack typical nucleotide mimetic structures, a library of extended fragments was built by exhaustive in silico reactions. After docking, minimization, clustering, visual inspection of the top-ranked compounds, and evaluation of ease of synthetic accessibility, either the original compound or a close derivative was synthesized and tested against Pim-1. For compounds showing the highest degree of Pim-1 inhibition the binding mode was determined crystallographically. Following a structure-guided approach, these were further optimized in a subsequent design cycle improving the compound's initial affinity by several orders of magnitude while synthesizing only a comparatively modest number of derivatives. The combination of computational and experimental approaches resulted in the development of a reasonably potent, novel molecular scaffold for inhibition of Pim-1 that targets specific surface regions, such as the interaction with R122 and P123 of the hinge region, which has been less frequently investigated in similar studies.

Isoforms of GPR35 have distinct extracellular N-termini that allosterically modify receptor-transducer coupling and mediate intracellular pathway bias.

Within the intestine, the human G protein-coupled receptor (GPCR) GPR35 is involved in oncogenic signaling, bacterial infections, and inflammatory bowel disease. GPR35 is known to be expressed as two distinct isoforms that differ only in the length of their extracellular N-termini by 31 amino acids, but detailed insights into their functional differences are lacking. Through gene expression analysis in immune and gastrointestinal cells, we show that these isoforms emerge from distinct promoter usage and alternative splicing. Additionally, we employed optical assays in living cells to thoroughly profile both GPR35 isoforms for constitutive and ligand-induced activation and signaling of 10 different heterotrimeric G proteins, ligand-induced arrestin recruitment, and receptor internalization. Our results reveal that the extended N-terminus of the long isoform limits G protein activation yet elevates receptor-ß-arrestin interaction. To better understand the structural basis for this bias, we examined structural models of GPR35 and conducted experiments with mutants of both isoforms. We found that a proposed disulfide bridge between the N-terminus and extracellular loop 3, present in both isoforms, is crucial for constitutive G13 activation, while an additional cysteine contributed by the extended N-terminus of the long GPR35 isoform limits the extent of agonist-induced receptor-ß-arrestin2 interaction. The pharmacological profiles and mechanistic insights of our study provide clues for the future design of isoform-specific GPR35 ligands that selectively modulate GPR35-transducer interactions and allow for mechanism-based therapies against, for example, inflammatory bowel disease or bacterial infections of the gastrointestinal system.

Systematic investigation of PRMT6 substrate recognition reveals broad specificity with a preference for an RG motif or basic and bulky residues.

Protein arginine methyltransferase 6 (PRMT6) catalyses the asymmetric dimethylation of arginines on numerous substrate proteins within the human cell. In particular, PRMT6 methylates histone H3 arginine 2 (H3R2) which affects both gene repression and activation. However, the substrate specificity of PRMT6 has not been comprehensively analysed. Here, we systematically characterise the substrate recognition motif of PRMT6, finding that it has broad specificity and recognises the RG motif. Working with a H3 tail peptide as a template, on which we made 204 amino acid substitutions, we use targeted mass spectrometry to measure their effect on PRMT6 in vitro activity. We first show that PRMT6 methylates R2 and R8 in the H3 peptide, although H3R8 is methylated with lower efficiency and is not an in vivo PRMT6 substrate. We then quantify the effect of 194 of these amino acid substitutions on methylation at both H3R2 and H3R8. In both cases, we find that PRMT6 tolerates essentially any amino acid substitution in the H3 peptide, but that positively charged and bulky residues are preferred near the target arginine. We show that PRMT6 also has preference for glycine, but only in the position immediately following the target arginine. This indicates that PRMT6 recognises the RG motif rather than the RGG motif. We further confirm this preference for the RG motif on another PRMT6 substrate, histone H4R3. This broad specificity and recognition of RG rather than RGG are distinctive among the PRMT family and has implications for the development of drugs to selectively target PRMT6. DATABASES: Panorama Public (https://panoramaweb.org/PRMT6motif.url); ProteomeXchange (PXD016711).

Analyzing Kinase Similarity in Small Molecule and Protein Structural Space to Explore the Limits of Multi-Target Screening.

While selective inhibition is one of the key assets for a small molecule drug, many diseases can only be tackled by simultaneous inhibition of several proteins. An example where achieving selectivity is especially challenging are ligands targeting human kinases. This difficulty arises from the high structural conservation of the kinase ATP binding sites, the area targeted by most inhibitors. We investigated the possibility to identify novel small molecule ligands with pre-defined binding profiles for a series of kinase targets and anti-targets by in silico docking. The candidate ligands originating from these calculations were assayed to determine their experimental binding profiles. Compared to previous studies, the acquired hit rates were low in this specific setup, which aimed at not only selecting multi-target kinase ligands, but also designing out binding to anti-targets. Specifically, only a single profiled substance could be verified as a sub-micromolar, dual-specific EGFR/ErbB2 ligand that indeed avoided its selected anti-target BRAF. We subsequently re-analyzed our target choice and in silico strategy based on these findings, with a particular emphasis on the hit rates that can be expected from a given target combination. To that end, we supplemented the structure-based docking calculations with bioinformatic considerations of binding pocket sequence and structure similarity as well as ligand-centric comparisons of kinases. Taken together, our results provide a multi-faceted picture of how pocket space can determine the success of docking in multi-target drug discovery efforts.

A Taxicab geometry quantification system to evaluate the performance of in silico methods: a case study on adenosine receptors ligands.

Among still comparatively few G protein-coupled receptors, the adenosine A2A receptor has been co-crystallized with several ligands, agonists as well as antagonists. It can thus serve as a template with a well-described orthosteric ligand binding region for adenosine receptors. As not all subtypes have been crystallized yet, and in order to investigate the usability of homology models in this context, multiple adenosine A1 receptor (A1AR) homology models had been previously obtained and a library of lead-like compounds had been docked. As a result, a number of potent and one selective ligand toward the intended target have been identified. However, in in vitro experimental verification studies, many ligands also bound to the A2AAR and the A3AR subtypes. In this work we asked the question whether a classification of the ligands according to their selectivity was possible based on docking scores. Therefore, we built an A3AR homology model and docked all previously found ligands to all three receptor subtypes. As a metric, we employed an in vitro/in silico selectivity ranking system based on taxicab geometry and obtained a classification model with reasonable separation. In the next step, the method was validated with an external library of, selective ligands with similarly good performance. This classification system might also be useful in further screens.

Direct targeting of Gαq and Gα11 oncoproteins in cancer cells.

Somatic gain-of-function mutations of GNAQ and GNA11, which encode α subunits of heterotrimeric Gαq/11 proteins, occur in about 85% of cases of uveal melanoma (UM), the most common cancer of the adult eye. Molecular therapies to directly target these oncoproteins are lacking, and current treatment options rely on radiation, surgery, or inhibition of effector molecules downstream of these G proteins. A hallmark feature of oncogenic Gαq/11 proteins is their reduced intrinsic rate of hydrolysis of guanosine triphosphate (GTP), which results in their accumulation in the GTP-bound, active state. Here, we report that the cyclic depsipeptide FR900359 (FR) directly interacted with GTPase-deficient Gαq/11 proteins and preferentially inhibited mitogenic ERK signaling rather than canonical phospholipase Cß (PLCß) signaling driven by these oncogenes. Thereby, FR suppressed the proliferation of melanoma cells in culture and inhibited the growth of Gαq-driven UM mouse xenografts in vivo. In contrast, FR did not affect tumor growth when xenografts carried mutated B-RafV600E as the oncogenic driver. Because FR enabled suppression of malignant traits in cancer cells that are driven by activating mutations at codon 209 in Gαq/11 proteins, we envision that similar approaches could be taken to blunt the signaling of non-Gαq/11 G proteins.

Three stories on Eph kinase inhibitors: From in silico discovery to in vivo validation.

Several selective and potent EphB4 inhibitors have been discovered, optimized and biophysically characterized by our groups over the past years. On the outset of these discoveries high throughput docking techniques were applied. Herein, we review the optimization campaigns started from three of these hits (Xan-A1, Pyr-A1 and Qui-A1) with emphasis on their in depth in vitro and in vivo characterization, together with previously unpublished angiogenesis and fluorescence based assays.

Identifying modulators of CXC receptors 3 and 4 with tailored selectivity using multi-target docking.

The G protein-coupled receptors of the C-X-C subfamily form a group among the chemokine receptors whose endogenous ligands are peptides with a common Cys-X-Cys motif. The CXC chemokine receptors 3 and 4 (CXCR3, CXCR4), which are investigated in this study, are linked to severe diseases such as cancer, multiple sclerosis, and HIV infections. Of particular interest, this receptor pair potentially forms a target for a polypharmacological drug treatment. Considering known ligands from public databases, such dual binders have not been identified yet. We therefore applied large-scale docking to the structure of CXCR4 and a homology model of CXCR3 with the goal to predict such dual binders, as well as compounds selective for either one of the receptors. Using signaling and biochemical assays, we showed that more than 50% of these predictions were correct in each category, yielding ligands with excellent binding efficiencies. These results highlight that docking is a suitable tool for the identification of ligands with tailored binding profiles to GPCRs, even when using homology models. More importantly, we present novel CXCR3-CXCR4 dual modulators that might pave the road to understanding the mechanisms of polypharmacological inhibition of these receptors.

Limits of ligand selectivity from docking to models: in silico screening for A(1) adenosine receptor antagonists.

G protein-coupled receptors (GPCRs) are attractive targets for pharmaceutical research. With the recent determination of several GPCR X-ray structures, the applicability of structure-based computational methods for ligand identification, such as docking, has increased. Yet, as only about 1% of GPCRs have a known structure, receptor homology modeling remains necessary. In order to investigate the usability of homology models and the inherent selectivity of a particular model in relation to close homologs, we constructed multiple homology models for the A(1) adenosine receptor (A(1)AR) and docked ∼2.2 M lead-like compounds. High-ranking molecules were tested on the A(1)AR as well as the close homologs A(2A)AR and A(3)AR. While the screen yielded numerous potent and novel ligands (hit rate 21% and highest affinity of 400 nM), it delivered few selective compounds. Moreover, most compounds appeared in the top ranks of only one model. These findings have implications for future screens.

The hunt for 8-oxoguanine deaminase.

An enzyme from Pseudomonas aeruginosa, Pa0142 (gi|9945972), that is able to catalyze the deamination of 8-oxoguanine (8-oxoG) to uric acid has been identified for the first time. 8-Oxoguanine is formed by the oxidation of guanine residues within DNA by reactive oxygen species, and this lesion results in G:C to T:A transversions. The value of k(cat)/K(m) for the deamination of 8-oxoG by Pa0142 at pH 8.0 and 30 degrees C is 2.0 x 10(4) M(-1) s(-1). This enzyme can also catalyze the deamination of isocystosine and guanine at rates that are approximately an order of magnitude lower. The three-dimensional structure of a homologous enzyme (gi|44264246) from the Sargasso Sea has been determined by X-ray diffraction methods to a resolution of 2.2 A (PDB entry). The enzyme folds as a (beta/alpha)(8) barrel and is a member of the amidohydrolase superfamily with a single zinc in the active site. This enzyme catalyzes the deamination of 8-oxoG with a k(cat)/K(m) value of 2.7 x 10(5) M(-1) s(-1). Computational docking of potential high-energy intermediates for the deamination reaction to the X-ray crystal structure suggests that active-site binding of 8-oxoG is facilitated by hydrogen-bond interactions from a conserved glutamine that follows beta-strand 1 with the carbonyl group at C6, a conserved tyrosine that follows beta-strand 2 with N7, and a conserved cysteine residue that follows beta-strand 4 with the carbonyl group at C8. A bioinformatic analysis of available protein sequences suggests that approximately 200 other bacteria possess an enzyme capable of catalyzing the deamination of 8-oxoG.

A double-headed cathepsin B inhibitor devoid of warhead.

Most synthetic inhibitors of peptidases have been targeted to the active site for inhibiting catalysis through reversible competition with the substrate or by covalent modification of catalytic groups. Cathepsin B is unique among the cysteine peptidase for the presence of a flexible segment, known as the occluding loop, which can block the primed subsites of the substrate binding cleft. With the occluding loop in the open conformation cathepsin B acts as an endopeptidase, and it acts as an exopeptidase when the loop is closed. We have targeted the occluding loop of human cathepsin B at its surface, outside the catalytic center, using a high-throughput docking procedure. The aim was to identify inhibitors that would interact with the occluding loop thereby modulating enzyme activity without the help of chemical warheads against catalytic residues. From a large library of compounds, the in silico approach identified [2-[2-(2,4-dioxo-1,3-thiazolidin-3-yl)ethylamino]-2-oxoethyl] 2-(furan-2-carbonylamino) acetate, which fulfills the working hypothesis. This molecule possesses two distinct binding moieties and behaves as a reversible, double-headed competitive inhibitor of cathepsin B by excluding synthetic and protein substrates from the active center. The kinetic mechanism of inhibition suggests that the occluding loop is stabilized in its closed conformation, mainly by hydrogen bonds with the inhibitor, thus decreasing endoproteolytic activity of the enzyme. Furthermore, the dioxothiazolidine head of the compound sterically hinders binding of the C-terminal residue of substrates resulting in inhibition of the exopeptidase activity of cathepsin B in a physiopathologically relevant pH range.

Structure-based tailoring of compound libraries for high-throughput screening: discovery of novel EphB4 kinase inhibitors.

High-throughput docking is a computational tool frequently used to discover small-molecule inhibitors of enzymes or receptors of known three-dimensional structure. Because of the large number of molecules in chemical libraries, automatic procedures to prune multimillion compound collections are useful for high-throughput docking and necessary for in vitro screening. Here, we propose an anchor-based library tailoring approach (termed ALTA) to focus a chemical library by docking and prioritizing molecular fragments according to their binding energy which includes continuum electrostatics solvation. In principle, ALTA does not require prior knowledge of known inhibitors, but receptor-based pharmacophore information (hydrogen bonds with the hinge region) is additionally used here to identify molecules with optimal anchor fragments for the ATP-binding site of the EphB4 receptor tyrosine kinase. The 21,418 molecules of the focused library (from an initial collection of about 730,000) are docked into EphB4 and ranked by force-field-based energy including electrostatic solvation. Among the 43 compounds tested in vitro, eight molecules originating from two different anchors show low-micromolar activity in a fluorescence-based enzymatic assay. Four of them are active in a cell-based assay and are potential anti-angiogenic compounds.

Automatic and efficient decomposition of two-dimensional structures of small molecules for fragment-based high-throughput docking.

The computer program DAIM (Decomposition and Identification of Molecules) has been developed to automatically break up compounds in small-molecule libraries for fragment-based docking as well as database analysis. Here, DAIM is evaluated on 130 ligands derived from known crystal structures of ligand-protein complexes. The decomposition and a new fingerprint-based identification technique are used to select anchor fragments for docking. The docking results show that the DAIM selection is superior to size-based or random selection of fragments. To evaluate the usefulness for analyzing the fragment composition of a large library, DAIM is applied to a collection of about 1.85 million commercially available compounds. Interestingly, it is found that the set of most frequent cyclic and acyclic fragments originating from the decomposition of the 1.85 million molecules shows a large overlap with the most frequent fragments in a library of 5120 known drugs. DAIM has been successfully used in the in silico screening for inhibitors of beta-secretase and EphB4 kinase by fragment-based high-throughput docking. Possible future applications for de novo ligand design are briefly discussed.

Obsessive-compulsive disorder and rectal prolapse.

A 47-year-old woman with a long-standing history of obsessive-compulsive disorder relating to dirt and germs is presented. Her fear of developing bowel cancer led her to manually evacuate faeces from her rectum five times a day and to a resultant rectal prolapse. Treatment involved prolonged graduated exposure to the patient's feared contaminants and ritual avoidance. After 5 months of inpatient therapy, the patient reported a subjective 70% improvement in her symptoms. As her obsessive-compulsive disorder symptoms improved, the patient's rectal prolapse disappeared.

Automated docking of highly flexible ligands by genetic algorithms: a critical assessment.

An improved version of the fragment-based flexible ligand docking approach SEED-FFLD is tested on inhibitors of human immunodeficiency virus type 1 protease, human alpha-thrombin and the estrogen receptor beta. The docking results indicate that it is possible to correctly reproduce the binding mode of inhibitors with more than ten rotatable bonds if the strain in their covalent geometry upon binding is not large. A high degree of convergence towards a unique binding mode in multiple runs of the genetic algorithm is proposed as a necessary condition for successful docking.