Identification of essential operons with a rhamnose-inducible promoter in Burkholderia cenocepacia

Scanning of bacterial genomes to identify essential genes is of biological interest, for understanding the basic functions required for life, and of practical interest, for the identification of novel targets for new antimicrobial therapies. In particular, the lack of efficacious antimicrobial treatments for infections caused by the Burkholderia cepacia complex is causing high morbidity and mortality of cystic fibrosis patients and of patients with nosocomial infections. Here, we present a method based on delivery of the tightly regulated rhamnose-inducible promoter P(rhaB) for identifying essential genes and operons in Burkholderia cenocepacia. We demonstrate that different levels of gene expression can be achieved by using two vectors that deliver P(rhaB) at two different distances from the site of insertion. One of these vectors places P(rhaB) at the site of transposon insertion, while the other incorporates the enhanced green fluorescent protein gene (e-gfp) downstream from P(rhaB). This system allows us to identify essential genes and operons in B. cenocepacia and provides a new tool for systematically identifying and functionally characterizing essential genes at the genomic level.

Chemically engineering ligand selectivity at the free fatty acid receptor 2 based on pharmacological variation between species orthologs

When it is difficult to develop selective ligands within a family of related G-protein-coupled receptors (GPCRs), chemically engineered receptors activated solely by synthetic ligands (RASSLs) are useful alternatives for probing receptor function. In the present work, we explored whether a RASSL of the free fatty acid receptor 2 (FFA2) could be developed on the basis of pharmacological variation between species orthologs. For this, bovine FFA2 was characterized, revealing distinct ligand selectivity compared with human FFA2. Homology modeling and mutational analysis demonstrated a single mutation in human FFA2 of C4.57G resulted in a human FFA2 receptor with ligand selectivity similar to the bovine receptor. This was exploited to generate human FFA2-RASSL by the addition of a second mutation at a known orthosteric ligand interaction site, H6.55Q. The resulting FFA2-RASSL displayed a >100-fold loss of activity to endogenous ligands, while responding to the distinct ligand sorbic acid with pEC(50) values for inhibition of cAMP, 5.83 ± 0.11; Ca(2+) mobilization, 4.63 ± 0.05; ERK phosphorylation, 5.61 ± 0.06; and dynamic mass redistribution, 5.35 ± 0.06. This FFA2-RASSL will be useful in future studies on this receptor and demonstrates that exploitation of pharmacological variation between species orthologs is a powerful method to generate novel chemically engineered GPCRs.-Hudson, B. D., Christiansen, E., Tikhonova, I. G., Grundmann, M., Kostenis, E., Adams, D. R., Ulven, T., Milligan, G. Chemically engineering ligand selectivity at the free fatty acid receptor 2 based on pharmacological variation between species orthologs.

Extracellular Loop 2 of the Free Fatty Acid Receptor 2 Mediates Allosterism of a Phenylacetamide Ago-Allosteric Modulator

Allosteric agonists are powerful tools for exploring the pharmacology of closely related G protein-coupled receptors that have nonselective endogenous ligands, such as the short chain fatty acids at free fatty acid receptors 2 and 3 (FFA2/GPR43 and FFA3/GPR41, respectively). We explored the molecular mechanisms mediating the activity of 4-chloro-alpha-(1-methylethyl)-N-2-thiazolylbenzeneacetamide (4-CMTB), a recently described phenylacetamide allosteric agonist and allosteric modulator of endogenous ligand function at human FFA2, by combining our previous knowledge of the orthosteric binding site with targeted examination of 4-CMTB structure-activity relationships and mutagenesis and chimeric receptor generation. Here we show that 4-CMTB is a selective agonist for FFA2 that binds to a site distinct from the orthosteric site of the receptor. Ligand structure-activity relationship studies indicated that the N-thiazolyl amide is likely to provide hydrogen bond donor/acceptor interactions with the receptor. Substitution at Leu(173) or the exchange of the entire extracellular loop 2 of FFA2 with that of FFA3 was sufficient to reduce or ablate, respectively, allosteric communication between the endogenous and allosteric agonists. Thus, we conclude that extracellular loop 2 of human FFA2 is required for transduction of cooperative signaling between the orthosteric and an as-yet-undefined allosteric binding site of the FFA2 receptor that is occupied by 4-CMTB.

Selective Orthosteric Free Fatty Acid Receptor 2 (FFA2) Agonists IDENTIFICATION OF THE STRUCTURAL AND CHEMICAL REQUIREMENTS FOR SELECTIVE ACTIVATION OF FFA2 VERSUS FFA3

Free fatty acid receptor 2 (FFA2; GPR43) is a G protein-coupled seven-transmembrane receptor for short-chain fatty acids (SCFAs) that is implicated in inflammatory and metabolic disorders. The SCFA propionate has close to optimal ligand efficiency for FFA2 and can hence be considered as highly potent given its size. Propionate, however, does not discriminate between FFA2 and the closely related receptor FFA3 (GPR41). To identify FFA2-selective ligands and understand the molecular basis for FFA2 selectivity, a targeted library of small carboxylic acids was examined using holistic, label-free dynamic mass redistribution technology for primary screening and the receptor-proximal G protein [S-35] guanosine 5'-(3-O-thio) triphosphate activation, inositol phosphate, and cAMP accumulation assays for hit confirmation. Structure-activity relationship analysis allowed formulation of a general rule to predict selectivity for small carboxylic acids at the orthosteric binding site where ligands with substituted sp(3)-hybridized alpha-carbons preferentially activate FFA3, whereas ligands with sp(2)- or sp-hybridized alpha-carbons prefer FFA2. The orthosteric binding mode was verified by site-directed mutagenesis: replacement of orthosteric site arginine residues by alanine in FFA2 prevented ligand binding, and molecular modeling predicted the detailed mode of binding. Based on this, selective mutation of three residues to their non-conserved counterparts in FFA3 was sufficient to transfer FFA3 selectivity to FFA2. Thus, selective activation of FFA2 via the orthosteric site is achievable with rather small ligands, a finding with significant implications for the rational design of therapeutic compounds selectively targeting the SCFA receptors.

Bidirectional, iterative approach to the structural delineation of the functional "Chemoprint" in GPR40 for agonist recognition

GPR40, free fatty acid receptor 1 (FFAR1), is a member of the GPCR superfamily and a possible target for the treatment of type 2 diabetes. In this work, we conducted a bidirectional iterative investigation, including computational modeling and site-directed mutagenesis, aimed at delineating amino acid residues forming the functional "chemoprint" of GPR40 for agonist recognition. The computational and experimental studies revolved around the recognition of the potent synthetic agonist GW9508. Our experimentally supported model suggested that H137(4.56), R183(5.39), N244(6.55), and R258(7.35) are directly involved in interactions with the ligand. We have proposed a polarized NH-pi interaction between H137(4.56) and GW9508 as one of the contributing forces leading to the high potency of GW9508. The modeling approach presented in this work provides a general strategy for the exploration of receptor-ligand interactions in G-protein coupled receptors beginning prior to acquisition of experimental data.

P2Y(1) antagonists: Combining receptor-based modeling and QSAR for a quantitative prediction of the biological activity based on consensus scoring

P2Y(1) is an ADP-activated G protein-coupled receptor (GPCR). Its antagonists impede platelet aggregation in vivo and are potential antithrombotic agents. Combining ligand and structure-based modeling we generated a consensus model (LIST-CM) correlating antagonist structures with their potencies. We docked 45 antagonists into our rhodopsin-based human P2Y(1) homology model and calculated docking scores and free binding energies with the Linear Interaction Energy (LIE) method in continuum-solvent. The resulting alignment was also used to build QSAR based on CoMFA, CoMSIA, and molecular descriptors. To benefit from the strength of each technique and compensate for their limitations, we generated our LIST-CM with a PLS regression based on the predictions of each methodology. A test set featuring untested substituents was synthesized and assayed in inhibition of 2-MeSADP-stimulated PLC activity and in radioligand binding. LIST-CM outperformed internal and external predictivity of any individual model to predict accurately the potency of 75% of the test set.

Validated ligand binding sites in CCK receptors. Next step: Computer aided design of novel CCK ligands

Computer-aided drug design becomes an important part of G-protein coupled receptors (GPCR) drug discovery process that is applied for improving the efficiency of derivation and optimization of novel ligands. It represents the combination of methods that-use-structural information of a receptor binding site of known ligands to design new ligands. In this report, we give a brief description of ligand binding sites in cholecystokinin and gastrin receptors (CK1R and CCK2R) which were delineated using experimental and computational methods, and then, we show how the validated ligand binding sites can be used to design and improve novel ligands. The translation of the knowledge of ligand-binding sites of different GPCRs to computer-aided design of novel ligands is summarized.

Extracellular ionic locks determine variation in constitutive activity and ligand potency between species orthologs of the free fatty acid receptors FFA2 and FFA3

Free fatty acid receptors 2 and 3 (FFA2 and FFA3) are G protein-coupled receptors for short chain free fatty acids (SCFAs). They respond to the same set of endogenous ligands but with distinct rank-order of potency, such that acetate (C2) has been described as FFA2 selective while propionate (C3) is non-selective. Although C2 was confirmed to be selective for human FFA2 over FFA3, this ligand was not selective between the mouse orthologs. Moreover, although C3 was indeed not selective between the human orthologs it displayed clear selectivity for mouse FFA3 over mouse FFA2. This altered selectivity to C2 and C3 resulted from broad differences in SCFAs potency at the mouse orthologs. In studies to define the molecular basis for these observations marked variation in ligand-independent, constitutive activity was identified. The orthologs with higher potency for the SCFAs, human FFA2 and mouse FFA3, displayed high constitutive activity while the orthologs with lower potency for the agonist ligands, mouse FFA2 and human FFA3, did not. Sequence alignments of the 2nd extracellular loop identified single negatively charged residues in FFA2 and FFA3 not conserved between species and predicted to form ionic lock interactions with arginine residues within the FFA2 or FFA3 agonist binding pocket to regulate constitutive activity and SCFA potency. Reciprocal mutation of these residues between species orthologs resulted in the induction (or repression) of constitutive activity, and in most cases also yielded corresponding changes in SCFA potency.

flp-32 Ligand/Receptor Silencing Phenocopy Faster Plant Pathogenic Nematodes

Restrictions on nematicide usage underscore the need for novel control strategies for plant pathogenic nematodes such as Globodera pallida (potato cyst nematode) that impose a significant economic burden on plant cultivation activities. The nematode neuropeptide signalling system is an attractive resource for novel control targets as it plays a critical role in sensory and motor functions. The FMRFamide-like peptides (FLPs) form the largest and most diverse family of neuropeptides in invertebrates, and are structurally conserved across nematode species, highlighting the utility of the FLPergic system as a broad-spectrum control target. flp-32 is expressed widely across nematode species. This study investigates the role of flp-32 in G. pallida and shows that: (i) Gp-flp-32 encodes the peptide AMRNALVRFamide; (ii) Gp-flp-32 is expressed in the brain and ventral nerve cord of G. pallida; (iii) migration rate increases in Gp-flp-32-silenced worms; (iv) the ability of G. pallida to infect potato plant root systems is enhanced in Gp-flp-32-silenced worms; (v) a novel putative Gp-flp-32 receptor (Gp-flp-32R) is expressed in G. pallida; and, (vi) Gp-flp-32R silenced worms also display an increase in migration rate. This work demonstrates that Gp30 flp-32 plays an intrinsic role in the modulation of locomotory behaviour in G. pallida, and putatively interacts with at least one novel G-protein coupled receptor (Gp-flp-32R). This is the first functional characterisation of a parasitic nematode FLP-GPCR. © 2013 Atkinson et al.

Distinct CCK-2 Receptor Conformations Associated with β-Arrestin-2 Recruitment or Phospholipase-C Activation Revealed by a Biased Antagonist

Seven-transmembrane receptors (7TMRs), also termed G protein-coupled receptors (GPCRs), form the largest class of cell surface membrane receptors, involving several hundred members in the human genome. Near 30% of marketed pharmacological agents target 7TMRs. 7TMRs adopt multiple conformations upon agonist binding. Biased agonists, in contrast to non-biased agonists, are believed to stabilize conformations preferentially activating either G-protein- or ß-arrestin-dependent signalling pathways. However, proof that cognate conformations of receptors display structural differences within their binding site where biased agonism initiates, are still lacking. Here, we show that a non-biased agonist, cholecystokinin (CCK) induces conformational states of the CCK2R activating Gq-protein-dependent pathway (CCK2RG) or recruiting ß-arrestin2 (CCK2Rß) that are pharmacologically and structurally distinct. Two structurally unrelated antagonists competitively inhibited both pathways. A third ligand (GV150,013X), acted as a high affinity competitive antagonist on CCK2RG but was nearly inefficient as inhibitor of CCK2Rß. Several structural elements on both GV150,013X and in CCK2R binding cavity, which hinder binding of GV150,013X only to the CCK2Rß were identified. At last, proximity between two conserved amino acids from transmembrane helices 3 and 7 interacting through sulphur-aromatic interaction was shown to be crucial for selective stabilization of the CCK2Rß state. These data establish structural evidences for distinct conformations of a 7TMR associated with ß-arrestin-2 recruitment or G-protein coupling and validate relevance of the design of biased ligands able to selectively target each functional conformation of 7TMRs.

Defining the molecular basis for the first potent and selective orthosteric agonists of the FFA2 free fatty acid receptor

FFA2 is a G protein-coupled receptor that responds to short chain fatty acids (SCFAs) and has generated interest as a therapeutic target for metabolic and inflammatory conditions. However, definition of its functions has been slowed by a dearth of selective ligands that can distinguish it from the closely related FFA3. At present, the only selective ligands described for FFA2 suffer from either poor potency, altered signaling due to allosteric modes of action, or a lack of function at non-human orthologs of the receptor. To address the need for novel selective ligands, we synthesized two compounds potentially having FFA2 activity and examined the molecular basis of their function. These compounds were confirmed to be potent and selective FFA2 agonists that interact with the orthosteric binding site. A combination of ligand structure-activity relationship, pharmacological analysis, homology modeling, species ortholog comparisons and mutagenesis studies were then employed to define the molecular basis of selectivity and function of these ligands. From this, we identified key residues within both extracellular loop 2 (ECL2) and the transmembrane domain (TM) regions of FFA2 critical for ligand function. One of these ligands was active with reasonable potency at rodent orthologs of FFA2 and demonstrated the role of FFA2 in the regulation of lipolysis in murine 3T3-L1 adipocytes. Together, these findings describe the first potent and selective FFA2 orthosteric agonists and demonstrate key aspects of ligand interaction within the orthosteric binding site of FFA2 that will be invaluable in future ligand development at this receptor.

Addressing Selective Polypharmacology of Antipsychotic Drugs Targeting the Bioaminergic Receptors through the Receptor Dynamic Conformational Ensembles

Selective polypharmacology, where a drug acts on multiple rather than single molecular targets involved in a disease, emerges to develop a structure-based system biology approach to design drugs selectively targeting a disease-active protein network. We focus on the bioaminergic receptors that belong to the group of integral membrane signalling proteins coupled to the G protein and represent targets for therapeutic agents against schizophrenia and depression. Among them, it has been shown that the serotonin (5-HT2A and 5-HT6), dopamine (D2 and D3) receptors induce a cognition-enhancing effect (group 1), while the histamine (H1) and serotonin (5-HT2C) receptors lead to metabolic side effects and the 5-HT2B serotonin receptor causes pulmonary hypertension (group 2). Thus, the problem arises to develop an approach that allows identifying drugs targeting only the disease-active receptors, i.e. group 1. The recent release of several crystal structures of the bioaminergic receptors, involving the D3 and H1 receptors provides the possibility to model the structures of all receptors and initiate a study of the structural and dynamic context of selective polypharmacology. In this work, we use molecular dynamics simulations to generate a conformational space of the receptors and subsequently characterize its binding properties applying molecular probe mapping. All-against-all comparison of the generated probe maps of the selected diverse conformations of all receptors with the Tanimoto similarity coefficient (Tc) enable to separate the receptors of group 1 from group 2. The pharmacophore built based on the Tc-selected receptor conformations, using the multiple probe maps discovers structural features that can be used to design molecules selective towards the receptors of group 1. The importance of several predicted residues to ligand selectivity is supported by the available mutagenesis and ligand structure-activity relationships studies. In addition, the Tc-selected conformations of the receptors for group 1 show good performance in isolation of known ligands from a random decoy. Our computational structure-based protocol to tackle selective polypharmacology of antipsychotic drugs could be applied for other diseases involving multiple drug targets, such as oncologic and infectious disorders.

Simulations of Biased Agonists in the β2 Adrenergic Receptor with Accelerated Molecular Dynamics

The biased agonism of the G protein-coupled receptors (GPCRs), where in addition to a traditional G protein-signalling pathway a GPCR promotes intracellular signals though ß-arrestin, is a novel paradigm in pharmacology. Biochemical and biophysical studies have suggested that a GPCR forms a distinct ensemble of conformations signalling through the G protein and ß-arrestin. Here we report on the dynamics of the ß2 adrenergic receptor bound to the ß-arrestin and G protein biased agonists and the empty receptor to further characterize the receptor conformational changes caused by biased agonists. We use conventional and accelerated molecular dynamics (aMD) simulations to explore the conformational transitions of the GPCR from the active state to the inactive state. We found that aMD simulations enable monitoring the transition within the nanosecond timescale while capturing the known microscopic characteristics of the inactive states, such as the ionic lock, the inward position of F6.44, and water clusters. Distinct conformational states are shown to be stabilized by each biased agonist. In particular, in simulations of the receptor with the ß-arrestin biased agonist, N-cyclopentylbutanepherine we observe a different pattern of motions in helix 7 when compared to simulations with the G protein biased agonist, Salbutamol that involves perturbations of the network of interactions within the NPxxY motif. Understanding the network of interactions induced by biased ligands and the subsequent receptor conformational shifts will lead to development of more efficient drugs. © 2013 American Chemical Society