A Ligand predication tool based on modeling and reasoning with imprecise probabilistic knowledge

Ligand prediction has been driven by a fundamental desire to understand more about how biomolecules recognize their ligands and by the commercial imperative to develop new drugs. Most of the current available software systems are very complex and time-consuming to use. Therefore, developing simple and efficient tools to perform initial screening of interesting compounds is an appealing idea. In this paper, we introduce our tool for very rapid screening for likely ligands (either substrates or inhibitors) based on reasoning with imprecise probabilistic knowledge elicited from past experiments. Probabilistic knowledge is input to the system via a user-friendly interface showing a base compound structure. A prediction of whether a particular compound is a substrate is queried against the acquired probabilistic knowledge base and a probability is returned as an indication of the prediction. This tool will be particularly useful in situations where a number of similar compounds have been screened experimentally, but information is not available for all possible members of that group of compounds. We use two case studies to demonstrate how to use the tool.

Identification of Determinants of Glucose-Dependent Insulinotropic Polypeptide Receptor That Interact with N-Terminal Biologically Active Region of the Natural Ligand

Glucose-dependent insulinotropic polypeptide receptor (GIPR), a member of family B of the G-protein coupled receptors, is a potential therapeutic target for which discovery of nonpeptide ligands is highly desirable. Structure-activity relationship studies indicated that the N-terminal part of glucose-dependent insulinotropic polypeptide (GIP) is crucial for biological activity. Here, we aimed at identification of residues in the GIPR involved in functional interaction with N-terminal moiety of GIP. A homology model of the transmembrane core of GIPR was constructed, whereas a three-dimensional model of the complex formed between GIP and the N-terminal extracellular domain of GIPR was taken from the crystal structure. The latter complex was docked to the transmembrane domains of GIPR, allowing in silico identification of putative residues of the agonist binding/activation site. All mutants were expressed at the surface of human embryonic kidney 293 cells as indicated by flow cytometry and confocal microscopy analysis of fluorescent GIP binding. Mutation of residues Arg183, Arg190, Arg300, and Phe357 caused shifts of 76-, 71-, 42-, and 16-fold in the potency to induce cAMP formation, respectively. Further characterization of these mutants, including tests with alanine-substituted GIP analogs, were in agreement with interaction of Glu3 in GIP with Arg183 in GIPR. Furthermore, they strongly supported a binding mode of GIP to GIPR in which the N-terminal moiety of GIP was sited within transmembrane helices (TMH) 2, 3, 5, and 6 with biologically crucial Tyr1 interacting with Gln224 (TMH3), Arg300 (TMH5), and Phe357 (TMH6). These data represent an important step toward understanding activation of GIPR by GIP, which should facilitate the rational design of therapeutic agents.

Ligand and structure-based methodologies for the prediction of the activity of G protein-coupled receptor ligands

Accurate in silico models for the quantitative prediction of the activity of G protein-coupled receptor (GPCR) ligands would greatly facilitate the process of drug discovery and development. Several methodologies have been developed based on the properties of the ligands, the direct study of the receptor-ligand interactions, or a combination of both approaches. Ligand-based three-dimensional quantitative structure-activity relationships (3D-QSAR) techniques, not requiring knowledge of the receptor structure, have been historically the first to be applied to the prediction of the activity of GPCR ligands. They are generally endowed with robustness and good ranking ability; however they are highly dependent on training sets. Structure-based techniques generally do not provide the level of accuracy necessary to yield meaningful rankings when applied to GPCR homology models. However, they are essentially independent from training sets and have a sufficient level of accuracy to allow an effective discrimination between binders and nonbinders, thus qualifying as viable lead discovery tools. The combination of ligand and structure-based methodologies in the form of receptor-based 3D-QSAR and ligand and structure-based consensus models results in robust and accurate quantitative predictions. The contribution of the structure-based component to these combined approaches is expected to become more substantial and effective in the future, as more sophisticated scoring functions are developed and more detailed structural information on GPCRs is gathered.

Discovery of substituted maleimides as liver X receptor agonists and determination of a ligand-bound crystal structure

Substituted 3-(phenylamino)-1H-pyrrole-2,5-diones were identified from a high throughput screen as inducers of human ATP binding cassette transporter A1 expression. Mechanism of action studies led to the identification of GSK3987 (4) as an LXR ligand. 4 recruits the steroid receptor coactivator-1 to human LXR alpha and LXRP with EC(50)s of 40 nM, profiles as an LXR agonist in functional assays, and activates LXR though a mechanism that is similar to first generation LXR agonists.

Identification of high-affinity binding sites for the insulin sensitizer rosiglitazone (BRL-49653) in rodent and human adipocytes using a radioiodinated ligand for peroxisomal proliferator-activated receptor gamma.

A radioiodinated ligand, [125I]SB-236636 [(S)-(-)3-[4-[2-[N-(2-benzoxazolyl)-N-methylamino]ethoxy]3-[125I]iodophenyl]2-ethoxy propanoic acid], which is specific for the ? isoform of the peroxisomal proliferator activated receptor (PPAR?), was developed. [125I]SB-236636 binds with high affinity to full-length human recombinant PPAR?1 and to a GST (glutathione S-transferase) fusion protein contg. the ligand binding domain of human PPAR?1 (KD = 70 nM). Using this ligand, the authors characterized binding sites in adipose-derived cells from rat, mouse and humans. In competition expts., rosiglitazone (BRL-49653), a potent antihyperglycemic agent, binds with high affinity to sites in intact adipocytes (IC50 = 12, 4 and 9 nM for rat, 3T3-L1 and human adipocytes, resp.). Binding affinities (IC50) of other thiazolidinediones for the ligand binding domain of PPAR?1 were comparable with those detd. in adipocytes and reflected the rank order of potencies of these agents as stimulants of glucose transport in 3T3-L1 adipocytes and antihyperglycemic agents in vivo: rosiglitazone > pioglitazone > troglitazone. Competition of [125I]SB-236636 binding was stereoselective in that the IC50 value of SB-219994, the (S)-enantiomer of an ?-trifluoroethoxy propanoic acid insulin sensitizer, was 770-fold lower than that of SB-219993 [(R)-enantiomer] at recombinant human PPAR?1. The higher binding affinity of SB-219994 also was evident in intact adipocytes and reflected its 100-fold greater potency as an antidiabetic agent. The results strongly suggest that the high-affinity binding site for [125I]SB-236636 in intact adipocytes is PPAR? and that the pharmacol. of insulin-sensitizer binding in rodent and human adipocytes is very similar and, moreover, predictive of antihyperglycemic activity in vivo.

Thrombopoietin, flt3-ligand and c-kit-ligand modulate HOX gene expression in expanding cord blood CD133(+) cells

Haemopoietic stem/progenitor cell (HSPC) development is regulated by extrinsic and intrinsic stimuli. Extrinsic modulators include growth factors and cell adhesion molecules, whereas intrinsic regulation is achieved with many transcription factor families, of which the HOX gene products are known to be important in haemopoiesis. Umbilical cord blood CD133(+) HSPC proliferation potential was tested in liquid culture with 'TPOFLK' (thrombopoietin, flt-3 ligand and c-kit ligand, promoting HSPC survival and self-renewal), in comparison to 'K36EG' (c-kit-ligand, interleukins-3 and -6, erythropoietin and granulocyte colony-stimulating factor, inducing haemopoietic differentiation). TPOFLK induced a higher CD133(+) HSPC proliferation (up to 60-fold more, at week 8) and maintained a higher frequency of the primitive colony-forming cells than K36EG. Quantitative polymerase chain reaction analysis revealed opposite expression patterns for specific HOX genes in expanding cord blood CD133(+) HSPC. After 8 weeks in liquid culture, TPOFLK increased the expression of HOX B3, B4 and A9 (associated with uncommitted HSPC) and reduced the expression of HOX B8 and A10 (expressed in committed myeloid cells) when compared to K36EG. These results suggest that TPOFLK induces CD133(+) HSPC proliferation, self-renewal and maintenance, up-regulation of HOX B3, B4 and A9 and down-regulation of HOX B8 and A10 gene expression.

ADP ribose is an endogenous ligand for the purinergic P2Y1 receptor

The mechanism by which extracellular ADP ribose (ADPr) increases intracellular free Ca2+ concentration ([Ca2+](i)) remains unknown. We measured [Ca2+](i) changes in fura-2 loaded rat insulinoma INS-1E cells, and in primary beta-cells from rat and human. A phosphonate analogue of ADPr (PADPr) and 8-Bromo-ADPr (8Br-ADPr) were synthesized. ADPr increased [Ca2+](i) in the form of a peak followed by a plateau dependent on extracellular Ca2+. NAD(+), cADPr, PADPr, 8Br-ADPr or breakdown products of ADPr did not increase [Ca2+](i). The ADPr-induced [Ca2+](i) increase was not affected by inhibitors of TRPM2, but was abolished by thapsigargin and inhibited when phospholipase C and IP3 receptors were inhibited. MRS 2179 and MRS 2279, specific inhibitors of the purinergic receptor P2Y1, completely blocked the ADPrinduced [Ca2+](i) increase. ADPr increased [Ca2+](i) in transfected human astrocytoma cells (1321N1) that express human P2Y1 receptors, but not in untransfected astrocytoma cells. We conclude that ADPr is a specific agonist of P2Y1 receptors. (c) 2010 Elsevier Ireland Ltd. All rights reserved.

A structural analysis of G-quadruplex/ligand interactions

This focused review article discusses in detail, all available high-resolution small molecule ligand/G-quadruplex structural data derived from crystallographic and NMR based techniques, in an attempt to understand key factors in ligand binding and to highlight the biological importance of these complexes. In contrast to duplex DNA, G-quadruplexes are four-stranded nucleic acid structures folded from guanine rich repeat sequences stabilized by the stacking of guanine G-quartets and extensive Watson-Crick/Hoogsteen hydrogen bonding. Thermally stable, these topologies can play a role in telomere regulation and gene expression. The core structures of G-quadruplexes form stable scaffolds while the loops have been shown, by the addition of small molecule ligands, to be sufficiently adaptable to generate new and extended binding platforms for ligands to associate, either by extending G-quartet surfaces or by forming additional planar dinucleotide pairings. Many of these structurally characterised loop rearrangements were totally unexpected opening up new opportunities for the design of selective ligands. However these rearrangements do significantly complicate attempts to rationally design ligands against well defined but unbound topologies, as seen for the series of napthalene diimides complexes. Drawing together previous findings and with the introduction of two new crystallographic quadruplex/ligand structures we aim to expand the understanding of possible structural adaptations available to quadruplexes in the presence of ligands, thereby aiding in the design of new selective entities. (C) 2011 Elsevier Masson SAS. All rights reserved.

Sequential C-H and C-Ru bond formation and cleavage during the thermally induced rearrangement of aryl ruthenium(II) complexes with [C6H3(CH2NMe2)(2)-2,6](-) as a bidentate eta(2)-C,N coordinated ligand. The crystal structures of the isomeric pairs [RuCl{eta(6)-C10H14}{eta(2)-C,N-C6H3(CH2NMe2)(2)-2,n}] (n = 4 or 6) and [Ru(eta(5)-C5H5){(eta(2)-C,N-C6H3(CH2NMe2)(2)-2,n} (PPH3)] (n = 4 or 6)

New air-stable ruthenium(II) complexes that contain the aryldiamine ligand [C6H3(CH2-NMe2)(2)-2,6](-) (NCN) are described. These complexes are [RuCl{eta(2)-C,N-C6H3(CH2NMe2)(2)-2,6}(eta(6)-C10H14)] (2; C10H14 = p-cymene = C6H4Me-Pr-i-4), [Ru{eta(2)-C,N-C6H3(CH2NMe2)(2)-2,6}(eta(5)-C5H5)(PPh3)] (5), and their isomeric forms [RuCl{eta(2)-C,N-C6H3(CH2NMe2)(2)-2,4}(eta(6)-C10H14)] (3) and [Ru{eta(2)-C,N-C6H3(CH2NMe2)(2)-2,4}(eta(5)-C5H5)(PPh3)] (6), respectively. Complex 2 has been prepared from the reaction of [Li(NCN)](2) with [RuCl2(eta(6)-C10H14)](2), whereas complex 5 has been prepared by the treatment of [RuCl{eta(3)-N,C,N-C6H3(CH2NMe2)(2)-2,6}(PPh3)] (4) with [Na(C5H5)](n). Both 2 and 5 are formally 18-electron ruthenium(II) complexes in which the monoanionic potentially tridentate coordinating ligand NCN is eta(2)-C,N-bonded, In solution (halocarbon solvent at room temperature or in aromatic solvents at elevated temperature), the intramolecular rearrangements of 2 and 5 afford complexes 3 and 6, respectively. This is a result of a shift of the metal-C-aryl bond from position-1 to position-3 on the aromatic ring of the NCN ligand. The mechanism of the isomerization is proposed to involve a sequence of intramolecular oxidative addition and reductive elimination reactions of both aromatic and aliphatic C-H bonds. This is based on results from deuterium labeling, spectroscopic studies, and some kinetic experiments. The mechanism is proposed to contain fully reversible steps in the case of 5, but a nonreversible step involving oxidative addition of a methyl NCH2-H bond in the case of 2. The solid-state structures of complexes 2, 3, 5, and 6 have been determined by single-crystal X-ray diffraction. A new dinuclear 1,4-phenylene-bridged bisruthenium(II) complex, [1,4-{RuCl(eta(6)-C10H14)}(2){C-6(CH2NMe2)(4)-2,3,5,6-C,N,C',N'}] (9) has also been prepared from the dianionic ligand [C-6(CH2NMe2)(4)-2,3,5,6](2-) (C2N4). The C2N4 ligand is in an eta(2)-C,N-eta(2)-C',N'-bis(bidentate) bonding mode. Compound 9 does not isomerize in solution (halocarbon solvent), presumably because of the absence of an accessible C-aryl-H bond. Complex 9 could not be isolated in an analytically pure form, probably because of its high sensitivity to air and very low solubility, which precludes recrystallization.