Search

Biased agonists, which selectively stimulate certain signaling pathways controlled by a G protein-coupled receptor (GPCR), hold great promise as drugs that maximize efficacy while minimizing dangerous side effects. Biased agonists of the μ-opioid receptor (μOR) are of particular interest as a means to achieve analgesia through G protein signaling without dose-limiting side effects such as respiratory depression and constipation. Rational structure-based design of biased agonists remains highly challenging, however, because the ligand-mediated interactions that are key to activation of each signaling pathway remain unclear. We identify several compounds for which the R- and S-enantiomers have distinct bias profiles at the μOR. These compounds serve as excellent comparative tools to study bias because the identical physicochemical properties of enantiomer pairs ensure that differences in bias profiles are due to differences in interactions with the μOR binding pocket. Atomic-level simulations of compounds at μOR indicate that R- and S-enantiomers adopt different poses that form distinct interactions with the binding pocket. A handful of specific interactions with highly conserved binding pocket residues appear to be responsible for substantial differences in arrestin recruitment between enantiomers. Our results offer guidance for rational design of biased agonists at μOR and possibly at related GPCRs.

The development of G protein-biased agonists for the μ-opioid receptor (MOR) offers a clear drug discovery rationale for improved analgesia and reduced side-effects of opiate pharmacotherapy. However, our understanding of the molecular mechanisms governing ligand bias is limited, which hinders our ability to rationally design biased compounds. We have investigated the role of MOR binding site residues W320 and Y328 in controlling bias, by receptor mutagenesis. The pharmacology of a panel of ligands in a cAMP and a β-arrestin2 assay were compared between the wildtype and mutated receptors, with bias factors calculated by operational analysis using ΔΔlog(τ/KA) values. [3H]diprenorphine competition binding was used to estimate affinity changes. Introducing the mutations W320A and Y328F caused changes in pathway bias, with different patterns of change between ligands. For example, DAMGO increased relative β-arrestin2 activity at the W320A mutant, whilst its β-arrestin2 response was completely lost at Y328F. In contrast, endomorphin-1 gained activity with Y328F but lost activity at W320A, in both pathways. For endomorphin-2 there was a directional shift from cAMP bias at the wildtype towards more β-arrestin2 bias at W320A. We also observe clear uncoupling between mutation-driven changes in function and binding affinity. These findings suggest that the mutations influenced the balance of pathway activation in a ligand-specific manner, thus identifying residues in the MOR binding pocket that govern ligand bias. This increases our understanding of how ligand/receptor binding interactions can be translated into agonist-specific pathway activation.

Valproic acid (VPA) is an anticonvulsant drug that is also used to treat migraines and bipolar disorder. Its proposed biological targets include human voltage-gated sodium channels, among other membrane proteins. We used the prokaryotic NavMs sodium channel, which has been shown to be a good exemplar for drug binding to human sodium channels, to examine the structural and functional interactions of VPA. Thermal melt synchrotron radiation circular dichroism spectroscopic binding studies of the full-length NavMs channel (which includes both pore and voltage sensor domains), and a pore-only construct, undertaken in the presence and absence of VPA, indicated that the drug binds to and destabilizes the channel, but not the pore-only construct. This is in contrast to other antiepileptic compounds that have previously been shown to bind in the central hydrophobic core of the pore region of the channel, and that tend to increase the thermal stability of both pore-only constructs and full-length channels. Molecular docking studies also indicated that the VPA binding site is associated with the voltage sensor, rather than the hydrophobic cavity of the pore domain. Electrophysiological studies show that VPA influences the block and inactivation rates of the NavMs channel, although with lower efficacy than classical channel-blocking compounds. It thus appears that, while VPA is capable of binding to these voltage-gated sodium channels, it has a very different mode and site of action than other anticonvulsant compounds.

Diseases such as chronic pain with complex etiologies are unlikely to respond to single, target-specific therapeutics but rather require intervention at multiple points within a perturbed disease system. Such approaches are being enabled by the rise of computational methods to identify key points of intervention and by new screening techniques that focus on a relevant condition or phenotype, rather than a specific target. Here we apply an in silico network pharmacology approach to identify small-molecule compounds with the potential to selectively disrupt the structure of a chronic-pain specific disease network, which we validate using a novel phenotypic screen that recapitulates key aspects of neuronal and pain biology by measuring changes in neuronal excitability in native sensory neurons. The combination of network pharmacology with a phenotypic screen is a powerful approach; we show that hit rates increase from 26% to 42%. This represents a rational approach to the discovery of compounds with a poly-pharmacology based therapeutic value, which will be vital for the discovery of treatments for complex disease.

Summary Small conductance potassium (SK) ion channels define neuronal firing rates by conducting the after-hyperpolarization current. They are key targets in developing therapies where neuronal firing rates are dysfunctional, such as in epilepsy, Parkinson's, and amyotrophic lateral sclerosis (ALS). Here, we characterize a binding pocket situated at the intracellular interface of SK2 and calmodulin, which we show to be shared by multiple small-molecule chemotypes. Crystallization of this complex revealed that riluzole (approved for ALS) and an analog of the anti-ataxic agent (4-chloro-phenyl)-[2-(3,5-dimethyl-pyrazol-1-yl)-pyrimidin-4-yl]-amine (CyPPA) bind to and allosterically modulate via this site. Solution-state nuclear magnetic resonance demonstrates that riluzole, NS309, and CyPPA analogs bind at this bipartite pocket. We demonstrate, by patch-clamp electrophysiology, that both classes of ligand interact with overlapping but distinct residues within this pocket. These data define a clinically important site, laying the foundations for further studies of the mechanism of action of riluzole and related molecules.

Cholinergic hypofunction is associated with decreased attention and cognitive deficits in the central nervous system in addition to compromised motor function. Consequently, stimulation of cholinergic neurotransmission is a rational therapeutic approach for the potential treatment of a variety of neurological conditions. High affinity choline uptake (HACU) into acetylcholine (ACh)-synthesizing neurons is critically mediated by the sodium- and pH-dependent high-affinity choline transporter (CHT, encoded by the SLC5A7 gene). This transporter is comparatively well-characterized but otherwise unexplored as a potential drug target. We therefore sought to identify small molecules that would enable testing of the hypothesis that positive modulation of CHT mediated transport would enhance activity-dependent cholinergic signaling. We utilized existing and novel screening techniques for their ability to reveal both positive and negative modulation of CHT using literature tools. A screening campaign was initiated with a bespoke compound library comprising both the Pfizer Chemogenomic Library (CGL) of 2,753 molecules designed specifically to help enable the elucidation of new mechanisms in phenotypic screens and 887 compounds from a virtual screening campaign to select molecules with field-based similarities to reported negative and positive allosteric modulators. We identified a number of previously unknown active and structurally distinct molecules that could be used as tools to further explore CHT biology or as a starting point for further medicinal chemistry.

Selective modulators of the γ-amino butyric acid (GABAA) family of receptors have the potential to treat a range of disease states related to cognition, pain, and anxiety. While the development of various α subunit-selective modulators is currently underway for the treatment of anxiety disorders, a mechanistic understanding of the correlation between their bioactivity and efficacy, based on ligand–target interactions, is currently still lacking. In order to alleviate this situation, in the current study we have analyzed, using ligand- and structure-based methods, a data set of 5440 GABAA modulators. The Spearman correlation (ρ) between binding activity and efficacy of compounds was calculated to be 0.008 and 0.31 against the α1 and α2 subunits of GABA receptor, respectively; in other words, the compounds had little diversity in structure and bioactivity, but they differed significantly in efficacy. Two compounds were selected as a case study for detailed interaction analysis due to the small difference in their structures and affinities (ΔpKi(comp1_α1 – comp2_α1) = 0.45 log units, ΔpKi(comp1_α2 – comp2_α2) = 0 log units) as compared to larger relative efficacies (ΔRE(comp1_α1 – comp2_α1) = 1.03, ΔRE(comp1_α2 – comp2_α2) = 0.21). Docking analysis suggested that His-101 is involved in a characteristic interaction of the α1 receptor with both compounds 1 and 2. Residues such as Phe-77, Thr-142, Asn-60, and Arg-144 of the γ chain of the α1γ2 complex also showed interactions with heterocyclic rings of both compounds 1 and 2, but these interactions were disturbed in the case of α2γ2 complex docking results. Binding pocket stability analysis based on molecular dynamics identified three substitutions in the loop C region of the α2 subunit, namely, G200E, I201T, and V202I, causing a reduction in the flexibility of α2 compared to α1. These amino acids in α2, as compared to α1, were also observed to decrease the vibrational and dihedral entropy and to increase the hydrogen bond content in α2 in the apo state. However, freezing of both α1 and α2 was observed in the ligand-bound state, with an increased number of internal hydrogen bonds and increased entropy. Therefore, we hypothesize that the amino acid differences in the loop C region of α2 are responsible for conformational changes in the protein structure compared to α1, as well as for the binding modes of compounds and hence their functional signaling.

The DNA Polymerase α (Pol α)/primase complex initiates DNA synthesis in eukaryotic replication. In the complex, Pol α and primase cooperate in the production of RNA-DNA oligonucleotides that prime synthesis of new DNA. Here we report crystal structures of the catalytic core of yeast Pol α in unliganded form, bound to an RNA primer/DNA template and extending an RNA primer with deoxynucleotides. We combine the structural analysis with biochemical and computational data to demonstrate that Pol α specifically recognizes the A-form RNA/DNA helix and that the ensuing synthesis of B-form DNA terminates primer synthesis. The spontaneous release of the completed RNA-DNA primer by the Pol α/primase complex simplifies current models of primer transfer to leading- and lagging strand polymerases. The proposed mechanism of nucleotide polymerization by Pol α might contribute to genomic stability by limiting the amount of inaccurate DNA to be corrected at the start of each Okazaki fragment. DOI: http://dx.doi.org/10.7554/eLife.00482.001, eLife digest During mitosis, a cell duplicates its DNA and then divides, ultimately generating two genetically identical daughter cells. In eukaryotes, the process of DNA duplication occurs at multiple sites throughout the genome: at each site, the antiparallel strands of the parental DNA separate and provide a template for DNA polymerase (Pol), the enzyme that synthesizes the two new DNA strands. Duplication of the DNA proceeds in both directions from each site through the polymerization of nucleotides to form new strands of DNA that are complementary to the template strands. However, since DNA polymerases can only polymerize nucleotides in one direction, the 5′ to 3′ direction, synthesis of the so-called leading strand proceeds continuously, whereas the other, lagging strand is synthesized in fragments. The task of duplicating the bulk of the DNA is shared between Pol δ, which is primarily responsible for synthesis of the lagging strand, and Pol ε, which fulfils the same role for the leading strand. However, Pols δ and ε cannot initiate DNA synthesis by themselves; short RNA-DNA chains called primers must also be paired to each template strand. Production of the primers requires the concerted action of two more enzymes: an RNA polymerase known as primase, and another DNA polymerase called Pol α. It is known that completion of the RNA-DNA primer requires Pol α to increase the length of the RNA segment by adding extra nucleotides, but the details of this process are poorly understood. Perera et al. combined crystallographic, biochemical and computational evidence to describe how Pol α first recognizes and then extends the RNA strand in the primer. They found that Pol α recognizes the particular shape of double helix—an A-form helix—that is formed by the DNA template and the RNA primer. The geometry of this helix prompts the Pol α enzyme to start adding nucleotides to the RNA in the primer. Perera et al. determined that once a full turn of double-helix DNA has been synthesized, Pol α is no longer in direct contact with the A-form helix, which causes the enzyme to disengage and terminate polymerization, leaving behind the now complete RNA-DNA primer. Perera et al. offer a new paradigm for understanding the initiation of DNA synthesis in eukaryotic replication. Their work suggests that Pol α has the ability to discriminate between different shapes of the primer-template helix, thus providing a mechanistic understanding of primer release. The spontaneous release of the primer offers a simple and elegant way to limit DNA synthesis by Pol α, a polymerase that is prone to error, and to make the RNA-DNA primer directly available for extension by Pol δ and Pol ε. DOI: http://dx.doi.org/10.7554/eLife.00482.002

Fibroblast growth factors (FGFs) are recognized targets for the development of therapies against angiogenesis-driven diseases, including cancer. The formation of a ternary complex with the transmembrane tyrosine kinase receptors (FGFRs), and heparan sulphate proteoglycans (HSPGs) is required for FGF2 pro-angiogenic activity. Here by using a combination of techniques including Nuclear Magnetic Resonance, Molecular Dynamics, Surface Plasmon Resonance and cell-based binding assays we clarify the molecular mechanism of inhibition of an angiostatic small molecule, sm27, mimicking the endogenous inhibitor of angiogenesis, thrombospondin-1. NMR and MD data demonstrate that sm27 engages the heparin-binding site of FGF2 and induces long-range dynamics perturbations along FGF2/FGFR1 interface regions. The functional consequence of the inhibitor binding is an impaired FGF2 interaction with both its receptors, as demonstrated by SPR and cell-based binding assays. We propose that sm27 antiangiogenic activity is based on a twofold-direct and allosteric-mechanism, inhibiting FGF2 binding to both its receptors.

G protein-coupled receptors (GPCRs) are a major family of membrane receptors in eukaryotic cells and play a crucial role in various biological processes. They represent a family of protein targets with significant therapeutic value, and accordingly more than 30% of prescription drugs are GPCR ligands [1]. Extending previous attempts to map the pharmacological space solely based on ligand chemical similarity,[2,3] we in this work relate GPCRs pharmacological space by combining structure-activity data from ChEMBL and WOMBAT that covers 167 human GPCRs and 67k ligands. By including more information from the ligand side in our analysis than previous studies, we hence attempted to construct a more detailed map of the pharmacological space. A statistical approach similar to the "Similarity Ensemble Approach" (SEA)[2] was implemented to relate proteins based on the chemical similarity of their ligands, and to rank the significance of the resulting similarity scores. A prospective external validation dataset was then employed to confirm new relationship between ligands and different GPCRs, providing mechanistic evidence for observed side effects of drugs in the dataset. The results of the study aim to contribute to a better understanding of the overlap of GPCRs in chemical space, and to the cross-reactivity observed even among distant biological targets, as defined by their sequence similarities[4] . Relevant applications range from understanding drug side effects to the design of drugs with a desired polypharmacological profile.