Discovery of Novel Adenosine Receptor Agonists That Exhibit Subtype Selectivity

A series of N6-bicyclic and N6-(2-hydroxy)cyclopentyl derivatives of adenosine were synthesized as novel A1R agonists and their A1R/A2R selectivity assessed using a simple yeast screening platform. We observed that the most selective, high potency ligands were achieved through N6-adamantyl substitution in combination with 5′-N-ethylcarboxamido or 5′-hydroxymethyl groups. In addition, we determined that 5′-(2-fluoro)thiophenyl derivatives all failed to generate a signaling response despite showing an interaction with the A1R. Some selected compounds were also tested on A1R and A3R in mammalian cells revealing that four of them are entirely A1R-selective agonists. By using in silico homology modeling and ligand docking, we provide insight into their mechanisms of recognition and activation of the A1R. We believe that given the broad tissue distribution, but contrasting signaling profiles, of adenosine receptor subtypes, these compounds might have therapeutic potential.

Cell-cell fusion induced by the Ig3 domain of receptor FGFRL1 in CHO cells.

FGFRL1 is a single-pass transmembrane protein with three extracellular Ig domains. When overexpressed in CHO cells or related cell types, it induces cell-cell fusion and formation of large, multinucleated syncytia. For this fusion-promoting activity, only the membrane-proximal Ig domain (Ig3) and the transmembrane domain are required. It does not matter whether the transmembrane domain is derived from FGFRL1 or from another receptor, but the distance of the Ig3 domain to the membrane is crucial. Fusion can be inhibited with soluble recombinant proteins comprising the Ig1-Ig2-Ig3 or the Ig2-Ig3 domains as well as with monoclonal antibodies directed against Ig3. Mutational analysis reveals a hydrophobic site in Ig3 that is required for fusion. If a single amino acid from this site is mutated, fusion is abolished. The site is located on a β-sheet, which is part of a larger β-barrel, as predicted by computer modeling of the 3D structure of FGFRL1. It is possible that this site interacts with a target protein of neighboring cells to trigger cell-cell fusion.

Molecular mechanisms of prostacyclin receptor signaling: The intracellular domains in G protein coupling

To better understand the mechanisms of how the human prostacyclin receptor (1P) mediates vasodilation and platelet anti-aggregation through Gs protein coupling, a strategy integrating multiple approaches including high resolution NMR experiments, synthetic peptide, fluorescence spectroscopy, molecular modeling, and recombinant protein was developed and used to characterize the structure/function relationship of important segments and residues of the IP receptor and the α-subunit of the Gs protein (Gαs). The first (iLP1) and third (iLP3) intracellular loops of the IP receptor, as well as the Gαs C-terminal domain, relevant to the Gs-mediated IP receptor signaling, were first identified by observation of the effects of the mini gene-expressed corresponding protein segments in HEK293 cells which co-expressed the receptor and Gαs. Evidence of the IP iLP1 domain interacted with the Gαs C-terminal domain was observed by fluorescence and NMR spectroscopic studies using a constrained synthetic peptide, which mimicked the IP iLP1 domain, and the synthetic peptide, which mimicked Gαs C-terminal domain. The solution structural models and the peptide-peptide interaction of the two synthetic protein segments were determined by high resolution NMR spectroscopy. The important residues in the corresponding domains of the IP receptor and the Gαs predicted by NMR chemical shift mapping were used to guide the identification of their protein-protein interaction in cells. A profile of the residues Arg42 - Ala48 of the IP iLP1 domain and the three residues Glu392 ∼ Leu394 of the Gαs C-terminal domain involved in the IP/Gs protein coupling were confirmed by recombinant proteins. The data revealed an intriguing speculation on the mechanisms of how the signal of the ligand-activated IP receptor is transmitted to the Gs protein in regulating vascular functions and homeostasis, and also provided substantial insights into other prostanoid receptor signaling. ^

Advanced Protein Modeling Method: Benchmarking and Applications in Computer-Aided Drug Discovery

Development of homology modeling methods will remain an area of active research. These methods aim to develop and model increasingly accurate three-dimensional structures of yet uncrystallized therapeutically relevant proteins e.g. Class A G-Protein Coupled Receptors. Incorporating protein flexibility is one way to achieve this goal. Here, I will discuss the enhancement and validation of the ligand-steered modeling, originally developed by Dr. Claudio Cavasotto, via cross modeling of the newly crystallized GPCR structures. This method uses known ligands and known experimental information to optimize relevant protein binding sites by incorporating protein flexibility. The ligand-steered models were able to model, reasonably reproduce binding sites and the co-crystallized native ligand poses of the β2 adrenergic and Adenosine 2A receptors using a single template structure. They also performed better than the choice of template, and crude models in a small scale high-throughput docking experiments and compound selectivity studies. Next, the application of this method to develop high-quality homology models of Cannabinoid Receptor 2, an emerging non-psychotic pain management target, is discussed. These models were validated by their ability to rationalize structure activity relationship data of two, inverse agonist and agonist, series of compounds. The method was also applied to improve the virtual screening performance of the β2 adrenergic crystal structure by optimizing the binding site using β2 specific compounds. These results show the feasibility of optimizing only the pharmacologically relevant protein binding sites and applicability to structure-based drug design projects.

Selectivity of diacylhydrazine insecticides to the predatory bug Orius laevigatus: In vivo and modeling/docking experiments

BACKGROUND: Knowledge of pesticide selectivity to natural enemies is necessary for a successful implementation of biological and chemical control methods in integrated pest management (IPM) programs. Diacylhydrazine (DAH)-based ecdysone agonists also known as molting-accelerating compounds (MACs) are considered a selective group of insecticides, and their compatibility with predatory Heteroptera, which are used as biological control agents, is known. However, their molecular mode of action has not been explored in beneficial insects such as Orius laevigatus (Fieber) (Hemiptera: Anthocoridae). RESULTS: In this project in vivo toxicity assays demonstrated that the DAH-based RH-5849, tebufenozide and methoxyfenozide have no toxic effect against O. laevigatus. The ligand-binding domain (LBD) of the ecdysone receptor (EcR) of O. laevigatus was sequenced and a homology protein model was constructed which confirmed a cavity structure with 12 ?-helixes, harboring the natural insect molting hormone 20-hydroxyecdysone. However, docking studies showed that a steric clash occurred for the DAH-based insecticides due to a restricted extent of the ligand-binding cavity of the EcR of O. laevigatus. CONCLUSIONS: The insect toxicity assays demonstrated that MACs are selective for O. laevigatus. The modeling/docking experiments are indications that these pesticides do not bind with the LBD-EcR of O. laevigatus and support that they show no biological effects in the predatory bug. These data help in explaining the compatible use of MACs together with predatory bugs in IPM programs. Keywords: Orius laevigatus, selectivity, diacylhydrazine insecticides, ecdysone receptor, homology modelling, docking studies.

Activating mutations in the extracellular domain of the fibroblast growth factor receptor 2 function by disruption of the disulfide bond in the third immunoglobulin-like domain

Multiple human skeletal and craniosynostosis disorders, including Crouzon, Pfeiffer, Jackson–Weiss, and Apert syndromes, result from numerous point mutations in the extracellular region of fibroblast growth factor receptor 2 (FGFR2). Many of these mutations create a free cysteine residue that potentially leads to abnormal disulfide bond formation and receptor activation; however, for noncysteine mutations, the mechanism of receptor activation remains unclear. We examined the effect of two of these mutations, W290G and T341P, on receptor dimerization and activation. These mutations resulted in cellular transformation when expressed as FGFR2/Neu chimeric receptors. Additionally, in full-length FGFR2, the mutations induced receptor dimerization and elevated levels of tyrosine kinase activity. Interestingly, transformation by the chimeric receptors, dimerization, and enhanced kinase activity were all abolished if either the W290G or the T341P mutation was expressed in conjunction with mutations that eliminate the disulfide bond in the third immunoglobulin-like domain (Ig-3). These results demonstrate a requirement for the Ig-3 cysteine residues in the activation of FGFR2 by noncysteine mutations. Molecular modeling also reveals that noncysteine mutations may activate FGFR2 by altering the conformation of the Ig-3 domain near the disulfide bond, preventing the formation of an intramolecular bond. This allows the unbonded cysteine residues to participate in intermolecular disulfide bonding, resulting in constitutive activation of the receptor.

Molecular basis for the recognition of two structurally different major histocompatibility complex/peptide complexes by a single T-cell receptor

2C is a typical alloreactive cytotoxic T lymphocyte clone that recognizes two different ligands. These ligands are adducts of the allo-major histocompatibility complex (MHC) molecule H-2Ld and an endogenous octapeptide, and of the self-MHC molecule H-2Kb and another peptide. MHC-binding and T-cell assays with synthetic peptides in combination with molecular modeling studies were employed to analyze the structural basis for this crossreactivity. The molecular surfaces of the two complexes differ greatly in densities and distributions of positive and negative charges. However, modifications of the peptides that increase similarity decrease the capacities of the resulting MHC peptide complexes to induce T-cell responses. Moreover, the roles of the peptides in ligand recognition are different for self- and allo-MHC-restricted T-cell responses. The self-MHC-restricted T-cell responses were finely tuned to recognition of the peptide. The allo-MHC-restricted responses, on the other hand, largely ignore modifications of the peptide. The results strongly suggest that adaptation of the T-cell receptor to the different ligand structures, rather than molecular mimicry by the ligands, is the basis for the crossreactivity of 2C. This conclusion has important implications for T-cell immunology and for the understanding of immunological disorders.

Genetically altered AMPA-type glutamate receptor kinetics in interneurons disrupt long-range synchrony of gamma oscillation

Gamma oscillations synchronized between distant neuronal populations may be critical for binding together brain regions devoted to common processing tasks. Network modeling predicts that such synchrony depends in part on the fast time course of excitatory postsynaptic potentials (EPSPs) in interneurons, and that even moderate slowing of this time course will disrupt synchrony. We generated mice with slowed interneuron EPSPs by gene targeting, in which the gene encoding the 67-kDa form of glutamic acid decarboxylase (GAD67) was altered to drive expression of the α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) glutamate receptor subunit GluR-B. GluR-B is a determinant of the relatively slow EPSPs in excitatory neurons and is normally expressed at low levels in γ-aminobutyric acid (GABA)ergic interneurons, but at high levels in the GAD-GluR-B mice. In both wild-type and GAD-GluR-B mice, tetanic stimuli evoked gamma oscillations that were indistinguishable in local field potential recordings. Remarkably, however, oscillation synchrony between spatially separated sites was severely disrupted in the mutant, in association with changes in interneuron firing patterns. The congruence between mouse and model suggests that the rapid time course of AMPA receptor-mediated EPSPs in interneurons might serve to allow gamma oscillations to synchronize over distance.

Structural analysis of p185c-neu and epidermal growth factor receptor tyrosine kinases: oligomerization of kinase domains.

The epidermal growth factor receptor (EGFR) and p185c-neu proteins associate as dimers to create an efficient signaling assembly. Overexpression of these receptors together enhances their intrinsic kinase activity and concomitantly results in oncogenic cellular transformation. The ectodomain is able to stabilize the dimer, whereas the kinase domain mediates biological activity. Here we analyze potential interactions of the cytoplasmic kinase domains of the EGFR and p185c-neu tyrosine kinases by homology molecular modeling. This analysis indicates that kinase domains can associate as dimers and, based on intermolecular interaction calculations, that heterodimer formation is favored over homodimers. The study also predicts that the self-autophosphorylation sites located within the kinase domains are not likely to interfere with tyrosine kinase activity, but may regulate the selection of substrates, thereby modulating signal transduction. In addition, the models suggest that the kinase domains of EGFR and p185c-neu can undergo higher order aggregation such as the formation of tetramers. Formation of tetrameric complexes may explain some of the experimentally observed features of their ligand affinity and hetero-receptor internalization.

D1 cap region involved in the receptor recognition and neural cell survival activity of human ciliary neurotrophic factor.

Human ciliary neurotrophic factor (hCNTF), which promotes the cell survival and differentiation of motor and other neurons, is a protein belonging structurally to the alpha-helical cytokine family. hCNTF was subjected to three-dimensional structure modeling and site-directed mutagenesis to analyze its structure-function relationship. The replacement of Lys-155 with any other amino acid residue resulted in abolishment of neural cell survival activity, and some of the Glu-153 mutant proteins had 5- to 10-fold higher biological activity. The D1 cap region (around the boundary between the CD loop and helix D) of hCNTF, including both Glu-153 and Lys-155, was shown to play a key role in the biological activity of hCNTF as one of the putative receptor-recognition sites. In this article, the D1 cap region of the 4-helix-bundle proteins is proposed to be important in receptor recognition and biological activity common to alpha-helical cytokine proteins reactive with gp130, a component protein of the receptors.

Compensatory T-Cell Regulation in Unaffected Relatives of SLE Patients, and Opposite IL-2/CD25-Mediated Effects Suggested by Coreferentiality Modeling

In human systemic lupus erythematosus (SLE), diverse autoantibodies accumulate over years before disease manifestation. Unaffected relatives of SLE patients frequently share a sustained production of autoantibodies with indiscriminable specificity, usually without ever acquiring the disease. We studied relations of IgG autoantibody profiles and peripheral blood activated regulatory T-cells (aTregs), represented by CD4(+)CD25(bright) T-cells that were regularly 70-90% Foxp3(+). We found consistent positive correlations of broad-range as well as specific SLE-associated IgG with aTreg frequencies within unaffected relatives, but not patients or unrelated controls. Our interpretation: unaffected relatives with shared genetic factors compensated pathogenic effects by aTregs engaged in parallel with the individual autoantibody production. To study this further, we applied a novel analytic approach named coreferentiality that tests the indirect relatedness of parameters in respect to multivariate phenotype data. Results show that independently of their direct correlation, aTreg frequencies and specific SLE-associated IgG were likely functionally related in unaffected relatives: they significantly parallelled each other in their relations to broad-range immunoblot autoantibody profiles. In unaffected relatives, we also found coreferential effects of genetic variation in the loci encoding IL-2 and CD25. A model of CD25 functional genetic effects constructed by coreferentiality maximization suggests that IL-2-CD25 interaction, likely stimulating aTregs in unaffected relatives, had an opposed effect in SLE patients, presumably triggering primarily T-effector cells in this group. Coreferentiality modeling as we do it here could also be useful in other contexts, particularly to explore combined functional genetic effects.

Potent cyclic antagonists of the complement C5a receptor on human polymorphonulcear leukocytes. Relationships between structures and activity

Human C5a is a plasma protein with potent chemoattractant and pro-inflammatory properties, and its overexpression correlates with severity of inflammatory diseases. C5a binds to its G protein-coupled receptor (C5aR) on polymorphonuclear leukocytes (PMNLs) through a high-affinity helical bundle and a low-affinity C terminus, the latter being solely responsible for receptor activation. Potent and selective C5a antagonists are predicted to be effective anti-inflammatory drugs, but no pharmacophore for small molecule antagonists has yet been developed, and it would significantly aid drug design. We have hypothesized that a turn conformation is important for activity of the C terminus of C5a and herein report small cyclic peptides that are stable turn mimics with potent antagonism at C5aR on human PMNLs. A comparison of solution structures for the C terminus of C5a, small acyclic peptide ligands, and cyclic antagonists supports the importance of a turn for receptor binding. Competition between a cyclic antagonist and either C5a or an acyclic agonist for C5aR on PMNLs supports a common or overlapping binding site on the C5aR. Structure-activity relationships for 60 cyclic analogs were evaluated by competitive radioligand binding with C5a (affinity) and myeloperoxidase release (antagonist potency) from human PMNLs, with 20 compounds having high antagonist potencies (IC50, 20 nM(-1) muM). Computer modeling comparisons reveal that potent antagonists share a common cyclic backbone shape, with affinity-determining side chains of defined volume projecting from the cyclic scaffold. These results define a new pharmacophore for C5a antagonist development and advance our understanding of ligand recognition and receptor activation of this G protein-coupled receptor.

Characterization of the structure of RAMP1 by mutagenesis and molecular modeling

Receptor activity modifying proteins (RAMPs) are a family of single-pass transmembrane proteins that dimerize with G-protein-coupled receptors. They may alter the ligand recognition properties of the receptors (particularly for the calcitonin receptor-like receptor, CLR). Very little structural information is available about RAMPs. Here, an ab initio model has been generated for the extracellular domain of RAMP1. The disulfide bond arrangement (Cys 27-Cys82, Cys40-Cys72, and Cys 57-Cys104) was determined by site-directed mutagenesis. The secondary structure (a-helices from residues 29-51, 60-80, and 87-100) was established from a consensus of predictive routines. Using these constraints, an assemblage of 25,000 structures was constructed and these were ranked using an all-atom statistical potential. The best 1000 conformations were energy minimized. The lowest scoring model was refined by molecular dynamics simulation. To validate our strategy, the same methods were applied to three proteins of known structure; PDB:1HP8, PDB:1V54 chain H (residues 21-85), and PDB:1T0P. When compared to the crystal structures, the models had root mean-square deviations of 3.8 Å, 4.1 Å, and 4.0 Å, respectively. The model of RAMP1 suggested that Phe93, Tyr 100, and Phe101 form a binding interface for CLR, whereas Trp74 and Phe92 may interact with ligands that bind to the CLR/RAMP1 heterodimer. © 2006 by the Biophysical Society.

A key role for transmembrane prolines in calcitonin receptor-like receptor agonist binding and signalling:implications for family B G-protein-coupled receptors

Calcitonin receptor like-receptor is a family B G-protein coupled receptor (GPCR). It requires receptor activity modifying protein (RAMP) 1 to give a calcitonin gene-related peptide (CGRP) receptor. Little is known of how members of this receptor family function. Proline residues often form important kinks in alpha-helices. Therefore, all proline residues within the transmembrane helices of the receptor (Pro241, Pro244 in helix 4, Pro275 in helix 5, Pro321 and Pro331 in helix 6) were mutated to alanine. Pro241 Pro275, and Pro321 are highly conserved throughout all family B GPCRs. The binding of CGRP and its ability to stimulate cAMP production were investigated in mutant and wild-type receptors after transient transfection into COS-7 cells with RAMP1. The P321A mutation significantly decreased the pEC(50) for CGRP and reduced its affinity but did not change cell-surface expression. Antagonist binding [CGRP(8-37) and 1-piperidinecarboxamide N-[2-[[5amino-1-[[4-(4-pyridinyl)-1-piperazinyl]carbonyl]pentyl]amino]-1-[(3 5-dibromo-4-hydroxyphenyl)methyl]-2-oxoethyl]-4-(1,4-dihydro-2-oxo-3(2H)-quina zolinyl) (BIBN4096BS)] was little altered by the mutation. Adrenomedullin-mediated signaling was disrupted when P321A was coexpressed with RAMP1, RAMP2, or RAMP3. The P331A mutant produced a moderate reduction in CGRP binding and receptor activation. Mutation of the other residues had no effect on receptor function. Thus, Pro321 and Pro331 are required for agonist binding and receptor activation. Modeling suggested that Pro321 induces a bend in helix 6, bringing its C terminus near that of helix 3, as seen in many family A GPCRs. This is abolished in P321A. P321A-I325P predicted to restore this conformation, showed wild-type activation. Modeling can also rationalize the effects of transmembrane proline mutants previously reported for another family B GPCR, the VPAC(1) receptor.

Diverse functional motifs within the three intracellular loops of the CGRP1 receptor

The CGRP1 receptor exists as a heterodimeric complex between a single-pass transmembrane accessory protein (RAMP1) and a family B G-protein-coupled receptor (GPCR) called the calcitonin receptor-like receptor (CLR). This study investigated the structural motifs found in the intracellular loops (ICLs) of this receptor. Molecular modeling was used to predict active and inactive conformations of each ICL. Conserved residues were altered to alanine by site-directed mutagenesis. cAMP accumulation, cell-surface expression, agonist affinity, and CGRP-stimulated receptor internalization were characterized. Within ICL1, L147 and particularly R151 were important for coupling to Gs. R151 may interact directly with the G-protein, accessing it following conformational changes involving ICL2 and ICL3. At the proximal end of ICL3, I290 and L294, probably lying on the same face of an α helix, formed a G-protein coupling motif. The largest effects on coupling were observed with I290A; additionally, it reduced CGRP affinity and impaired internalization. 1290 may interact with TM6 to stabilize the conformation of ICL3, but it could also interact directly with Gs. R314, at the distal end of ICL3, impaired G-protein coupling and to a lesser extent reduced CGRP affinity; it may stabilize the TM6-ICL3 junction by interacting with the polar headgroups of membrane phospholipids. Y215 and L214 in ICL2 are required for cell-surface expression; they form a microdomain with H216 which has the same function. This study reveals similarities between the activation of CLR and other GPCRs in the role of TM6 and ICL3 but shows that other conserved motifs differ in their function. © 2006 American Chemical Society.