Thomas Balle

Unraveling the High- and Low-Sensitivity Agonist Responses of Nicotinic Acetylcholine Receptors

The neuronal α4β2 nicotinic acetylcholine receptors exist as two distinct subtypes, (α4)2(β2)3 and (α4)3(β2)2, and biphasic responses to acetylcholine and other agonists have been ascribed previously to coexistence of these two receptor subtypes. We offer a novel and radical explanation for the observation of two distinct agonist sensitivities. Using different expression ratios of mammalian α4 and β2 subunits and concatenated constructs, we demonstrate that a biphasic response is an intrinsic functional property of the (α4)3(β2)2 receptor. In addition to two high-sensitivity sites at α4β2 interfaces, the (α4)3(β2)2 receptor contains a third low-sensitivity agonist binding site in the α4α4 interface. Occupation of this site is required for full activation and is responsible for the widened dynamic response range of this receptor subtype. By site-directed mutagenesis, we show that three residues, which differ between the α4β2 and α4α4 sites, control agonist sensitivity. The results presented here provide a basic insight into the function of pentameric ligand-gated ion channels, which enables modulation of the receptors with hitherto unseen precision; it becomes possible to rationally design therapeutics targeting subpopulations of specific receptor subtypes.

Open3DQSAR: a new open-source software aimed at high-throughput chemometric analysis of molecular interaction fields

AbstractOpen3DQSAR is a freely available open-source program aimed at chemometric analysis of molecular interaction fields. MIFs can be imported from different sources (GRID, CoMFA/CoMSIA, quantum-mechanical electrostatic potential or electron density grids) or generated by Open3DQSAR itself. Much focus has been put on automation through the implementation of a scriptable interface, as well as on high computational performance achieved by algorithm parallelization. Flexibility and interoperability with existing molecular modeling software make Open3DQSAR a powerful tool in pharmacophore assessment and ligand-based drug design. FigureOpen3DQSAR’s logo displaying PLS coefficient isocontours

A Unified Model of the GABAA Receptor Comprising Agonist and Benzodiazepine Binding Sites

We present a full-length α1β2γ2 GABA receptor model optimized for agonists and benzodiazepine (BZD) allosteric modulators. We propose binding hypotheses for the agonists GABA, muscimol and THIP and for the allosteric modulator diazepam (DZP). The receptor model is primarily based on the glutamate-gated chloride channel (GluCl) from C. elegans and includes additional structural information from the prokaryotic ligand-gated ion channel ELIC in a few regions. Available mutational data of the binding sites are well explained by the model and the proposed ligand binding poses. We suggest a GABA binding mode similar to the binding mode of glutamate in the GluCl X-ray structure. Key interactions are predicted with residues α1R66, β2T202, α1T129, β2E155, β2Y205 and the backbone of β2S156. Muscimol is predicted to bind similarly, however, with minor differences rationalized with quantum mechanical energy calculations. Muscimol key interactions are predicted to be α1R66, β2T202, α1T129, β2E155, β2Y205 and β2F200. Furthermore, we argue that a water molecule could mediate further interactions between muscimol and the backbone of β2S156 and β2Y157. DZP is predicted to bind with interactions comparable to those of the agonists in the orthosteric site. The carbonyl group of DZP is predicted to interact with two threonines α1T206 and γ2T142, similar to the acidic moiety of GABA. The chlorine atom of DZP is placed near the important α1H101 and the N-methyl group near α1Y159, α1T206, and α1Y209. We present a binding mode of DZP in which the pending phenyl moiety of DZP is buried in the binding pocket and thus shielded from solvent exposure. Our full length GABAA receptor is made available as Model S1.

Open3DALIGN: an open-source software aimed at unsupervised ligand alignment

An open-source, cross-platform software aimed at conformer generation and unsupervised rigid-body molecular alignment is presented. Different algorithms have been implemented to perform single and multi-conformation superimpositions on one or more templates. Alignments can be accomplished by matching pharmacophores, heavy atoms or a combination of the two. All methods have been successfully validated on eight comprehensive datasets previously gathered by Sutherland and co-workers. High computational performance has been attained through efficient parallelization of the code. The unsupervised nature of the alignment algorithms, together with its scriptable interface, make Open3DALIGN an ideal component of high-throughput, automated cheminformatics workflows.

Intersubunit Bridge Formation Governs Agonist Efficacy at Nicotinic Acetylcholine α4β2 Receptors: UNIQUE ROLE OF HALOGEN BONDING REVEALED

Background: Molecular features governing α4β2 nAChRs efficacy have remained elusive. Results: Binding studies, electrophysiology, and structural data from co-crystallization with Ls-AChBP are reported for a series of α4β2 agonists. Conclusion: Direct halogen bonds and an invariant Loop-C suggest that intersubunit bridge formation governs efficacy. Significance: The data provide a structural basis for understanding of efficacy levels at nAChRs. The α4β2 subtype of the nicotinic acetylcholine receptor has been pursued as a drug target for treatment of psychiatric and neurodegenerative disorders and smoking cessation aids for decades. Still, a thorough understanding of structure-function relationships of α4β2 agonists is lacking. Using binding experiments, electrophysiology and x-ray crystallography we have investigated a consecutive series of five prototypical pyridine-containing agonists derived from 1-(pyridin-3-yl)-1,4-diazepane. A correlation between binding affinities at α4β2 and the acetylcholine-binding protein from Lymnaea stagnalis (Ls-AChBP) confirms Ls-AChBP as structural surrogate for α4β2 receptors. Crystal structures of five agonists with efficacies at α4β2 from 21–76% were determined in complex with Ls-AChBP. No variation in closure of loop C is observed despite large efficacy variations. Instead, the efficacy of a compound appears tightly coupled to its ability to form a strong intersubunit bridge linking the primary and complementary binding interfaces. For the tested agonists, a specific halogen bond was observed to play a large role in establishing such strong intersubunit anchoring.

Rational design of alpha-conotoxin analogues targeting alpha7 nicotinic acetylcholine receptors: improved antagonistic activity by incorporation of proline derivatives.

Nicotinic acetylcholine receptors (nAChRs) are ligand-gated ion channels that belong to the superfamily of Cys loop receptors. Valuable insight into the orthosteric ligand binding to nAChRs in recent years has been obtained from the crystal structures of acetylcholine-binding proteins (AChBPs) that share significant sequence homology with the amino-terminal domains of the nAChRs. α-Conotoxins, which are isolated from the venom of carnivorous marine snails, selectively inhibit the signaling of neuronal nAChR subtypes. Co-crystal structures of α-conotoxins in complex with AChBP show that the side chain of a highly conserved proline residue in these toxins is oriented toward the hydrophobic binding pocket in the AChBP but does not have direct interactions with this pocket. In this study, we have designed and synthesized analogues of α-conotoxins ImI and PnIA[A10L], by introducing a range of substituents on the Pro6 residue in these toxins to probe the importance of this residue for their binding to the nAChRs. Pharmacological characterization of the toxin analogues at the α7 nAChR shows that although polar and charged groups on Pro6 result in analogues with significantly reduced antagonistic activities, analogues with aromatic and hydrophobic substituents in the Pro6 position exhibit moderate activity at the receptor. Interestingly, introduction of a 5-(R)-phenyl substituent at Pro6 in α-conotoxin ImI gives rise to a conotoxin analogue with a significantly higher binding affinity and antagonistic activity at the α7 nAChR than those exhibited by the native conotoxin.Nicotinic acetylcholine receptors (nAChRs) are ligand-gated ion channels that belong to the superfamily of Cys loop receptors. Valuable insight into the orthosteric ligand binding to nAChRs in recent years has been obtained from the crystal structures of acetylcholine-binding proteins (AChBPs) that share significant sequence homology with the amino-terminal domains of the nAChRs. alpha-Conotoxins, which are isolated from the venom of carnivorous marine snails, selectively inhibit the signaling of neuronal nAChR subtypes. Co-crystal structures of alpha-conotoxins in complex with AChBP show that the side chain of a highly conserved proline residue in these toxins is oriented toward the hydrophobic binding pocket in the AChBP but does not have direct interactions with this pocket. In this study, we have designed and synthesized analogues of alpha-conotoxins ImI and PnIA[A10L], by introducing a range of substituents on the Pro(6) residue in these toxins to probe the importance of this residue for their binding to the nAChRs. Pharmacological characterization of the toxin analogues at the alpha(7) nAChR shows that although polar and charged groups on Pro(6) result in analogues with significantly reduced antagonistic activities, analogues with aromatic and hydrophobic substituents in the Pro(6) position exhibit moderate activity at the receptor. Interestingly, introduction of a 5-(R)-phenyl substituent at Pro(6) in alpha-conotoxin ImI gives rise to a conotoxin analogue with a significantly higher binding affinity and antagonistic activity at the alpha(7) nAChR than those exhibited by the native conotoxin.

New insights into the GABAA receptor structure and orthosteric ligand binding: Receptor modeling guided by experimental data

GABAA receptors (GABAARs) are ligand gated chloride ion channels that mediate overall inhibitory signaling in the CNS. A detailed understanding of their structure is important to gain insights in, e.g., ligand binding and functional properties of this pharmaceutically important target. Homology modeling is a necessary tool in this regard because experimentally determined structures are lacking. Here we present an exhaustive approach for creating a high quality model of the α1β2γ2 subtype of the GABAAR ligand binding domain, and we demonstrate its usefulness in understanding details of orthosteric ligand binding. The model was constructed by using multiple templates and by incorporation of knowledge from biochemical/pharmacological experiments. It was validated on the basis of objective energy functions, its ability to account for available residue specific information, and its stability in molecular dynamics (MD) compared with that of the two homologous crystal structures. We then combined the model with extensive structure‐activity relationships available from two homologous series of orthosteric GABAAR antagonists to create a detailed hypothesis for their binding modes. Excellent agreement with key experimental data was found, including the ability of the model to accommodate and explain a previously developed pharmacophore model. A coupling to agonist binding was thereby established and discussed in relation to activation mechanisms. Our results highlight the importance of critical evaluation and optimization of each step in the homology modeling process. The approach taken here can greatly aid in increasing the understanding of GABAARs and related receptors where structural insight is limited and reliable models are difficult to obtain. Proteins 2011;

G protein- and agonist-bound serotonin 5-HT2A receptor model activated by steered molecular dynamics simulations.

A 5-HT(2A) receptor model was constructed by homology modeling based on the β(2)-adrenergic receptor and the G protein-bound opsin crystal structures. The 5-HT(2A) receptor model was transferred into an active conformation by an agonist ligand and a G(αq) peptide in four subsequent steered molecular dynamics (MD) simulations. The driving force for the transformation was the addition of several known intermolecular and receptor interhelical hydrogen bonds enforcing the necessary helical and rotameric movements. Subsquent MD simulations without constraints confirmed the stability of the activated receptor model as well as revealed new information about stabilizing residues and bonds. The active 5-HT(2A) receptor model was further validated by retrospective ligand screening of more than 9400 compounds, whereof 182 were known ligands. The results show that the model can be used in drug discovery for virtual screening and structure-based ligand design as well as in GPCR activation studies.

Ginkgolide X Is a Potent Antagonist of Anionic Cys-loop Receptors with a Unique Selectivity Profile at Glycine Receptors

The novel ginkgolide analog ginkgolide X was characterized functionally at human glycine and γ-aminobutyric acid type A receptors (GlyRs and GABAARs, respectively) in the fluorescence-based FLIPRTM Membrane Potential assay. The compound inhibited the signaling of all GABAAR subtypes included in the study with high nanomolar/low micromolar IC50 values, except the ρ1 receptor at which it was a significantly weaker antagonist. Ginkgolide X also displayed high nanomolar/low micromolar IC50 values at the homomeric α1 and α2 GlyRs, whereas it was inactive at the heteromeric α1β and α2β subtypes at concentrations up to 300 μm. Thus, the functional properties of the compound were significantly different from those of the naturally occurring ginkgolides A, B, C, J, and M but similar to those of picrotoxin. In a mutagenesis study the 6′ M2 residues in the GlyR ion channel were identified as the primary molecular determinant of the selectivity profile of ginkgolide X, and a 6′ M2 ring consisting of five Thr residues was found to be of key importance for its activity at the GABAAR. Conformational analysis and docking of low-energy conformations of the native ginkgolide A and ginkgolide X into a α1 GlyR homology model revealed two distinct putative binding sites formed by the 6′ M2 residues together with the 2′ residues and the 10′ and 13′ residues, respectively. Thus, we propose that the distinct functionalities of ginkgolide X compared with the other ginkgolides could arise from different flexibility and thus different binding modes to the ion channel of the anionic Cys-loop receptor.

Probing the Role of Backbone Hydrogen Bonds in Protein-Peptide Interactions by Amide-to-Ester Mutations

One of the most frequent protein-protein interaction modules in mammalian cells is the postsynaptic density 95/discs large/zonula occludens 1 (PDZ) domain, involved in scaffolding and signaling and emerging as an important drug target for several diseases. Like many other protein-protein interactions, those of the PDZ domain family involve formation of intermolecular hydrogen bonds: C-termini or internal linear motifs of proteins bind as β-strands to form an extended antiparallel β-sheet with the PDZ domain. Whereas extensive work has focused on the importance of the amino acid side chains of the protein ligand, the role of the backbone hydrogen bonds in the binding reaction is not known. Using amide-to-ester substitutions to perturb the backbone hydrogen-bonding pattern, we have systematically probed putative backbone hydrogen bonds between four different PDZ domains and peptides corresponding to natural protein ligands. Amide-to-ester mutations of the three C-terminal amides of the peptide ligand severely affected the affinity with the PDZ domain, demonstrating that hydrogen bonds contribute significantly to ligand binding (apparent changes in binding energy, ΔΔG = 1.3 to >3.8 kcal mol(-1)). This decrease in affinity was mainly due to an increase in the dissociation rate constant, but a significant decrease in the association rate constant was found for some amide-to-ester mutations suggesting that native hydrogen bonds have begun to form in the transition state of the binding reaction. This study provides a general framework for studying the role of backbone hydrogen bonds in protein-peptide interactions and for the first time specifically addresses these for PDZ domain-peptide interactions.