Using Physical Chemistry To Differentiate Nicotinic from Cholinergic Agonists at the Nicotinic Acetylcholine Receptor

Electrostatic Contributions to the Binding Free Energy of Nicotine to the Acetylcholine Binding Protein

Biomolecular binding relies on specific attractive interactions between two partner molecules, including electrostatics, dispersion, hydrophobicity, and solvation. Assessing the contributions of electrostatic interactions to binding is key to the understanding of ligand binding mechanisms and the design of improved biomolecular binders. For example, nicotine is a well-known agonist of nicotinic acetylcholine receptors (nAChRs), but the molecular mechanisms for the differential action of nicotine on brain and muscle nAChRs remain elusive. In this work, we have chosen the acetylcholine binding protein (AChBP) in complex with nicotine as a model system to interrogate the electrostatic contributions to nicotine binding. Our absolute binding free energy simulations confirm that nicotine binds AChBP predominantly in its protonated (charged) form. By comparing energetic contributions from decomposed interactions for either neutral or charged nicotine, our calculations shed light on the nature of the binding of nicotine to the AChBP. The preferred binding of charged nicotine over neutral nicotine originates from its stronger electrostatic interactions with AChBP, a cation-π interaction to a tryptophan residue and a hydrogen bond between nicotine and the backbone carbonyl of the tryptophan, whereas the major force driving the binding process appears to be van der Waals interactions. The various nonelectrostatic terms can also indirectly modulate the electrostatic interactions through fine-tuning the binding pose of the ligand in the binding site, providing an explanation of why the binding specificity of nicotine to the brain versus muscle nAChRs is driven by electrostatic interaction, given that the immediate binding site residues, including the key tryptophan residue, are identical in the two receptors.

Cobalt-catalyzed highly diastereoselective [3 + 2] carboannulation reactions: facile access to substituted indane derivatives

Efficient oxidative [3 + 2] annulation reaction involving aryl hydrazones and heterobicyclic alkenes has been realized with inexpensive and earth-abundant cobalt salts under aerobic conditions. The reaction proceeds via directing-group-assisted C-H activation and exo-selective migratory insertion, followed by the intramolecular nucleophilic attack of the alkylcobalt(III) species onto the imine with high anti-diastereoselectivity to provide complex indane derivatives. The generation of three contiguous stereogenic centers within the indanyl unit and the avoidance of the use of stoichiometric amounts of metal oxidants make this transformation more valuable and appealing.

Mechanism of biomolecular recognition of trimethyllysine by the fluorinated aromatic cage of KDM5A PHD3 finger

The understanding of biomolecular recognition of posttranslationally modified histone proteins is centrally important to the histone code hypothesis. Despite extensive binding and structural studies on the readout of histones, the molecular language by which posttranslational modifications on histone proteins are read remains poorly understood. Here we report physical-organic chemistry studies on the recognition of the positively charged trimethyllysine by the electron-rich aromatic cage containing PHD3 finger of KDM5A. The aromatic character of two tryptophan residues that solely constitute the aromatic cage of KDM5A was fine-tuned by the incorporation of fluorine substituents. Our thermodynamic analyses reveal that the wild-type and fluorinated KDM5A PHD3 fingers associate equally well with trimethyllysine. This work demonstrates that the biomolecular recognition of trimethyllysine by fluorinated aromatic cages is associated with weaker cation-π interactions that are compensated by the energetically more favourable trimethyllysine-mediated release of high-energy water molecules that occupy the aromatic cage.

A Methyl Scan of the Pyrrolidinium Ring of Nicotine Reveals Significant Differences in Its Interactions with α7 and α4β2 Nicotinic Acetylcholine Receptors

The two major nicotinic acetylcholine receptors (nAChRs) in the brain are the α4β2 and α7 subtypes. A "methyl scan" of the pyrrolidinium ring was used to detect differences in nicotine's interactions with these two receptors. Each methylnicotine was investigated using voltage-clamp and radioligand binding techniques. Methylation at each ring carbon elicited unique changes in nicotine's receptor interactions. Replacing the 1'-N-methyl with an ethyl group or adding a second 1'-N-methyl group significantly reduced interaction with α4β2 but not α7 receptors. The 2'-methylation uniquely enhanced binding and agonist potency at α7 receptors. Although 3'- and 5'-trans-methylations were much better tolerated by α7 receptors than α4β2 receptors, 4'-methylation decreased potency and efficacy at α7 receptors much more than at α4β2 receptors. Whereas cis-5'-methylnicotine lacked agonist activity and displayed a low affinity at both receptors, trans-5'-methylnicotine retained considerable α7 receptor activity. Differences between the two 5'-methylated analogs of the potent pyridyl oxymethylene-bridged nicotine analog A84543 were consistent with what was found for the 5'-methylnicotines. Computer docking of the methylnicotines to the Lymnaea acetylcholine binding protein crystal structure containing two persistent waters predicted most of the changes in receptor affinity that were observed with methylation, particularly the lower affinities of the cis-methylnicotines. The much smaller effects of 1'-, 3'-, and 5'-methylations and the greater effects of 2'- and 4'-methylations on nicotine α7 nAChR interaction might be exploited for the design of new drugs based on the nicotine scaffold. SIGNIFICANCE STATEMENT: Using a comprehensive "methyl scan" approach, we show that the orthosteric binding sites for acetylcholine and nicotine in the two major brain nicotinic acetylcholine receptors interact differently with the pyrrolidinium ring of nicotine, and we suggest reasons for the higher affinity of nicotine for the heteromeric receptor. Potential sites for nicotine structure modification were identified that may be useful in the design of new drugs targeting these receptors.

Specific Encapsulation of Acetylcholine Chloride by a Self-Assembled Molecular Capsule with Sulfonamido Moiety

New molecular capsule 1₂, which encapsulates only acetylcholine chloride in the ion-pair form, has been developed. Cavitand 1 with four sulfonamido moieties on the upper rim of tetraimino-cavitand self-assembled to form a stable molecular capsule in the presence of an acetylcholine chloride guest through eight intermolecular -NH···O═S hydrogen bonds, two from each of the four paired sulfonamide units.

Probing Backbone Hydrogen Bonds in Proteins by Amide-to-Ester Mutations

All proteins contain characteristic backbones formed of consecutive amide bonds, which can engage in hydrogen bonds. However, the importance of these is not easily addressed by conventional technologies that only allow for side-chain substitutions. By contrast, technologies such as nonsense suppression mutagenesis and protein ligation allow for manipulation of the protein backbone. In particular, replacing the backbone amide groups with ester groups, that is, amide-to-ester mutations, is a powerful tool to examine backbone-mediated hydrogen bonds. In this minireview, we showcase examples of how amide-to-ester mutations can be used to uncover pivotal roles of backbone-mediated hydrogen bonds in protein recognition, folding, function, and structure.

Absolute Configuration and Pharmacology of the Poison Frog Alkaloid Phantasmidine

Phantasmidine, a rigid congener of the well-known nicotinic acetylcholine receptor agonist epibatidine, is found in the same species of poison frog ( Epipedobates anthonyi). Natural phantasmidine was found to be a 4:1 scalemic mixture, enriched in the (2a R,4a S,9a S) enantiomer by chiral-phase LC-MS comparison to the synthetic enantiomers whose absolute configurations were previously established by Mosher's amide analysis. The major enantiomer has the opposite S configuration at the benzylic carbon to natural epibatidine, whose benzylic carbon is R. Pharmacological characterization of the synthetic racemate and separated enantiomers established that phantasmidine is ∼10-fold less potent than epibatidine, but ∼100-fold more potent than nicotine in most receptors tested. Unlike epibatidine, phantasmidine is sharply enantioselective in its activity and the major natural enantiomer whose benzylic carbon has the 4a S configuration is more active. The stereoselective pharmacology of phantasmidine is ascribed to its rigid and asymmetric shape as compared to the nearly symmetric conformations previously suggested for epibatidine enantiomers. While phantasmidine itself is too toxic for direct therapeutic use, we believe it is a useful platform for the development of potent and selective nicotinic agonists, which may have value as pharmacological tools.

Heterobicyclic Core Retained Hydroarylations through C-H Activation: Synthesis of Epibatidine Analogues

The heterobicyclic core retained hydroarylation of oxa/azabenzonorbornadienes with quinoline N-oxides has been achieved under rhodium catalysis, giving quinoline N-oxide substituted heterobicyclic structures with excellent regioselectivity and in good yields. As the first example of the direct introduction of quinoline N-oxides onto heterobicyclic structures, the strained heterobicyclic core was well retained in the reaction. The products could be successfully transformed into a series of useful compounds, including epibatidine analogues.

Probing for and Quantifying Agonist Hydrogen Bonds in α6β2 Nicotinic Acetylcholine Receptors

Designing subtype-selective agonists for neuronal nicotinic acetylcholine receptors is a challenging and significant goal aided by intricate knowledge of each subtype's binding patterns. We previously reported that in α6β2 receptors, acetylcholine makes a functional cation-π interaction with Trp149, but nicotine and TC299423 do not, suggesting a distinctive binding site. This work explores hydrogen binding at the backbone carbonyl associated with α6β2 Trp149. Substituting residue i + 1, Thr150, with its α-hydroxy analogue (Tah) attenuates the carbonyl's hydrogen bond accepting ability. At α6(T150Tah)β2, nicotine shows a 24-fold loss of function, TC299423 shows a modest loss, and acetylcholine shows no effect. Nicotine was further analyzed via a double-mutant cycle analysis utilizing N'-methylnicotinium, which indicated a hydrogen bond in α6β2 with a ΔΔG of 2.6 kcal/mol. Thus, even though nicotine does not make the conserved cation-π interaction with Trp149, it still makes a functional hydrogen bond to its associated backbone carbonyl.

Activation and Desensitization of Peripheral Muscle and Neuronal Nicotinic Acetylcholine Receptors by Selected, Naturally-Occurring Pyridine Alkaloids

Teratogenic alkaloids can cause developmental defects due to the inhibition of fetal movement that results from desensitization of fetal muscle-type nicotinic acetylcholine receptors (nAChRs). We investigated the ability of two known teratogens, the piperidinyl-pyridine anabasine and its 1,2-dehydropiperidinyl analog anabaseine, to activate and desensitize peripheral nAChRs expressed in TE-671 and SH-SY5Y cells. Activation-concentration response curves for each alkaloid were obtained in the same multi-well plate. To measure rapid desensitization, cells were first exposed to five potentially-desensitizing concentrations of each alkaloid in log₁₀ molar increments from 10 nM to 100 µM and then to a fixed concentration of acetylcholine (ACh), which alone produces near-maximal activation. The fifty percent desensitization concentration (DC₅₀) was calculated from the alkaloid concentration-ACh response curve. Agonist fast desensitization potency was predicted by the agonist potency measured in the initial response. Anabaseine was a more potent desensitizer than anabasine. Relative to anabaseine, nicotine was more potent to autonomic nAChRs, but less potent to the fetal neuromuscular nAChRs. Our experiments have demonstrated that anabaseine is more effective at desensitizing fetal muscle-type nAChRs than anabasine or nicotine and, thus, it is predicted to be more teratogenic.

Plant alkaloids that cause developmental defects through the disruption of cholinergic neurotransmission

The exposure of a developing embryo or fetus to alkaloids from plants, plant products, or plant extracts has the potential to cause developmental defects in humans and animals. These defects may have multiple causes, but those induced by piperidine and quinolizidine alkaloids arise from the inhibition of fetal movement and are generally referred to as multiple congenital contracture-type deformities. These skeletal deformities include arthrogyrposis, kyposis, lordosis, scoliosis, and torticollis, associated secondary defects, and cleft palate. Structure-function studies have shown that plant alkaloids with a piperidine ring and a minimum of a three-carbon side-chain α to the piperidine nitrogen are teratogenic. Further studies determined that an unsaturation in the piperidine ring, as occurs in gamma coniceine, or anabaseine, enhances the toxic and teratogenic activity, whereas the N-methyl derivatives are less potent. Enantiomers of the piperidine teratogens, coniine, ammodendrine, and anabasine, also exhibit differences in biological activity, as shown in cell culture studies, suggesting variability in the activity due to the optical rotation at the chiral center of these stereoisomers. In this article, we review the molecular mechanism at the nicotinic pharmacophore and biological activities, as it is currently understood, of a group of piperidine and quinolizidine alkaloid teratogens that impart a series of flexure-type skeletal defects and cleft palate in animals.

Spectroscopic, crystallographic and theoretical studies of lasalocid complex with ammonia and benzylamine

A natural antibiotic--Lasalocid is able to form stable complexes with ammonia and organic amines. New complexes of lasalocid with benzylamine and ammonia were obtained in the crystal forms and studied using X-ray, FT-IR, (1)H NMR, (13)C NMR and DFT methods. These studies have shown that in both complexes the proton is transferred from the carboxylic group to the amine group with the formation of a pseudo-cyclic structure of lasalocid anion complexing the protonated amine or NH4(+) cation. The spectroscopic and DFT studies demonstrated that the structure of the complex formed between Lasalocid and benzylamine in the solid is also conserved in the solution and gas phase. In contrast, the structure of the complex formed between lasalocid and ammonium cation found in the solid state undergoes dissociation in chloroform solution accompanied with a change in the coordination form of the NH4(+) cation.

Functional probes of drug-receptor interactions implicated by structural studies: Cys-loop receptors provide a fertile testing ground

Structures of integral membrane receptors provide valuable models for drug–receptor interactions across many important classes of drug targets and have become much more widely available in recent years. However, it remains to be determined to what extent these images are relevant to human receptors in their biological context and how subtle issues such as subtype selectivity can be informed by them. The high precision structural modifications enabled by unnatural amino acid mutagenesis on mammalian receptors expressed in vertebrate cells allow detailed tests of predictions from structural studies. Using the Cys-loop superfamily of ligand-gated ion channels, we show that functional studies lead to detailed binding models that, at times, are significantly at odds with the structural studies on related invertebrate proteins. Importantly, broad variations in binding interactions are seen for very closely related receptor subtypes and for varying drugs at a given binding site. These studies highlight the essential interplay between structural studies and functional studies that can guide efforts to develop new pharmaceuticals.

Probing the non-canonical interface for agonist interaction with an α5 containing nicotinic acetylcholine receptor

Nicotinic acetylcholine receptors (nAChRs) containing the α5 subunit are of interest because genome-wide association studies and candidate gene studies have identified polymorphisms in the α5 gene that are linked to an increased risk for nicotine dependence, lung cancer, and/or alcohol addiction. To probe the functional impact of an α5 subunit on nAChRs, a method to prepare a homogeneous population of α5-containing receptors must be developed. Here we use a gain of function (9') mutation to isolate populations of α5-containing nAChRs for characterization by electrophysiology. We find that the α5 subunit modulates nAChR rectification when co-assembled with α4 and β2 subunits. We also probe the α5-α4 interface for possible ligand-binding interactions. We find that mutations expected to ablate an agonist-binding site involving the α5 subunit have no impact on receptor function. The most straightforward interpretation of this observation is that agonists do not bind at the α5-α4 interface, in contrast to what has recently been demonstrated for the α4-α4 interface in related receptors. In addition, our mutational results suggest that the α5 subunit does not replace the α4 or β2 subunits and is relegated to occupying only the auxiliary position of the pentameric receptor.

An unusual pattern of ligand-receptor interactions for the α7 nicotinic acetylcholine receptor, with implications for the binding of varenicline

The α7 nicotinic acetylcholine receptor shows broad pharmacology, complicating the development of subtype-specific nicotinic receptor agonists. Here we use unnatural amino acid mutagenesis to characterize binding to α7 by the smoking cessation drug varenicline (Chantix; Pfizer, Groton, CT), an α4β2-targeted agonist that shows full efficacy and modest potency at the α7 receptor. We find that unlike binding to its target receptor, varenicline does not form a cation-π interaction with TrpB, further supporting a unique binding mode for the cationic amine of nicotinic agonists at the α7 receptor. We also evaluate binding to the complementary face of the receptor's binding site by varenicline, the endogenous agonist acetylcholine, and the potent nicotine analog epibatidine. Interestingly, we find no evidence for functionally important interactions involving backbone NH and CO groups thought to bind the canonical agonist hydrogen bond acceptor of the nicotinic pharmacophore, perhaps reflecting a lesser importance of this pharmacophore element for α7 binding. We also show that the Trp55 and Leu119 side chains of the binding site's complementary face are important for the binding of the larger agonists epibatidine and varenicline, but dispensable for binding of the smaller, endogenous agonist acetylcholine.

Binding interactions with the complementary subunit of nicotinic receptors

The agonist-binding site of nicotinic acetylcholine receptors (nAChRs) spans an interface between two subunits of the pentameric receptor. The principal component of this binding site is contributed by an α subunit, and it binds the cationic moiety of the nicotinic pharmacophore. The other part of the pharmacophore, a hydrogen bond acceptor, has recently been shown to bind to the complementary non-α subunit via the backbone NH of a conserved Leu. This interaction was predicted by studies of ACh-binding proteins and confirmed by functional studies of the neuronal (CNS) nAChR, α4β2. The ACh-binding protein structures further suggested that the hydrogen bond to the backbone NH is mediated by a water molecule and that a second hydrogen bonding interaction occurs between the water molecule and the backbone CO of a conserved Asn, also on the non-α subunit. Here, we provide new insights into the nature of the interactions between the hydrogen bond acceptor of nicotinic agonists and the complementary subunit backbone. We studied both the nAChR of the neuromuscular junction (muscle-type) and a neuronal subtype, (α4)2(β4)3. In the muscle-type receptor, both ACh and nicotine showed a strong interaction with the Leu NH, but the potent nicotine analog epibatidine did not. This interaction was much attenuated in the α4β4 receptor. Surprisingly, we found no evidence for a functionally significant interaction with the backbone carbonyl of the relevant Asn in either receptor with an array of agonists.

Variations in binding among several agonists at two stoichiometries of the neuronal, α4β2 nicotinic receptor

Drug-receptor binding interactions of four agonists, ACh, nicotine, and the smoking cessation compounds varenicline (Chantix) and cytisine (Tabex), have been evaluated at both the 2:3 and 3:2 stoichiometries of the α4β2 nicotinic acetylcholine receptor (nAChR). Previous studies have established that unnatural amino acid mutagenesis can probe three key binding interactions at the nAChR: a cation-π interaction, and two hydrogen-bonding interactions to the protein backbone of the receptor. We find that all drugs make a cation-π interaction to TrpB of the receptor. All drugs except ACh, which lacks an N(+)H group, make a hydrogen bond to a backbone carbonyl, and ACh and nicotine behave similarly in acting as a hydrogen-bond acceptor. However, varenicline is not a hydrogen-bond acceptor to the backbone NH that interacts strongly with the other three compounds considered. In addition, we see interesting variations in hydrogen bonding interactions with cytisine that provide a rationalization for the stoichiometry selectivity seen with this compound.

Probing the effects of residues located outside the agonist binding site on drug-receptor selectivity in the nicotinic receptor

The nicotinic acetylcholine receptors (nAChRs) are a family of closely related but pharmacologically distinct neurotransmitter-gated ion channels. They are therapeutic targets for a wide range of neurological disorders, and a key issue in drug development is selective targeting among the more than 20 subtypes of nAChRs that are known. The present work evaluates a proposed hydrogen bonding interaction involving a residue known as the "loop B glycine" that distinguishes receptors that are highly responsive to ACh and nicotine from those that are much less so. We have performed structure-function studies on the loop B site, including unnatural amino acid mutagenesis, in three different nAChR subtypes and found that the correlation between agonist potency and this residue is strong. Low potency receptor subtypes have a glycine at this key site, and mutation to a residue with a side chain converts a low potency receptor to a high potency receptor. Innately high potency receptors have a lysine at the loop B site and show a decrease in potency for the reverse mutation (i.e., introducing a glycine). This residue lies outside of the agonist binding site, and studies of other residues at the agonist binding site show that the details of how changes at the loop B glycine site impact agonist potency vary for differing receptor subtypes. This suggests a model in which the loop B residue influences the global shape of the agonist binding site rather than modulating any specific interaction.

On the interaction of aliphatic amines and ammonium ions with carboxylic acids in solution and in receptor pockets

Association energies of the acetate ion with cationic amines bearing one to three methyl groups were calculated in the range of -14 to -17 kcal/mol in aqueous solution by means of the IEF-PCM method at the CCSD(T)/CBS//MP2/aug-cc-pvdz and DFT/B97D/CBS//B97D/aug-cc-pvtz levels. The main stabilization factor for the association is the possibility for the formation of an ionic intermolecular hydrogen bond between the elements of the complex. For a quaternary ammonium ion, the favorable electrostatic interaction energy is the only driving force, and the stabilization energy for the complex is reduced to -4 kcal/mol. The internal free energies of the ion-pair tautomers of the studied species are higher by 10-15 kcal/mol in water than those for the neutral, hydrogen-bonded forms. Monte Carlo free energy perturbation calculations at T = 298 K and p = 1 atm predict -11 to -16 kcal/mol relative solvation free energy in favor of the corresponding ionic form. As a result, the ion-pair tautomer is the prevailing form in aqueous solution and on the extracellular surface of a receptor. Modeling the complex of a protonated ligand interacting with an Asp/Glu carboxylate side-chain in the binding cavity of a receptor, two strongly bound water molecules were considered so as to form hydrogen-bonded water bridges between the elements of the ion-pair. Nonetheless, the low polarity environment mimicked by a chloroform solvent cannot stabilize the ionic tautomer. A proton jump was predicted, which suggests that acetylcholine, an inherent cation by structure, might have evolved as the natural agonist for muscarinic receptors because a quaternary ammonium system assures the maintenance of the ion-pair form with a carboxylate side-chain in a protein cavity, the latter perhaps then being needed for receptor activation.

Two neuronal nicotinic acetylcholine receptors, alpha4beta4 and alpha7, show differential agonist binding modes

Nicotinic acetylcholine receptors (nAChRs) are pentameric, neurotransmitter-gated ion channels responsible for rapid excitatory neurotransmission in the central and peripheral nervous systems, resulting in skeletal muscle tone and various cognitive effects in the brain. These complex proteins are activated by the endogenous neurotransmitter ACh as well as by nicotine and structurally related agonists. Activation and modulation of nAChRs has been implicated in the pathology of multiple neurological disorders, and as such, these proteins are established therapeutic targets. Here we use unnatural amino acid mutagenesis to examine the ligand binding mechanisms of two homologous neuronal nAChRs: the α4β4 and α7 receptors. Despite sequence identity among the residues that form the core of the agonist-binding site, we find that the α4β4 and α7 nAChRs employ different agonist-receptor binding interactions in this region. The α4β4 receptor utilizes a strong cation-π interaction to a conserved tryptophan (TrpB) of the receptor for both ACh and nicotine, and nicotine participates in a strong hydrogen bond with a backbone carbonyl contributed by TrpB. Interestingly, we find that the α7 receptor also employs a cation-π interaction for ligand recognition, but the site has moved to a different aromatic amino acid of the agonist-binding site depending on the agonist. ACh participates in a cation-π interaction with TyrA, whereas epibatidine participates in a cation-π interaction with TyrC2.

The structural basis of function in Cys-loop receptors

Cys-loop receptors are membrane-spanning neurotransmitter-gated ion channels that are responsible for fast excitatory and inhibitory transmission in the peripheral and central nervous systems. The best studied members of the Cys-loop family are nACh, 5-HT3, GABAA and glycine receptors. All these receptors share a common structure of five subunits, pseudo-symmetrically arranged to form a rosette with a central ion-conducting pore. Some are cation selective (e.g. nACh and 5-HT3) and some are anion selective (e.g. GABAA and glycine). Each receptor has an extracellular domain (ECD) that contains the ligand-binding sites, a transmembrane domain (TMD) that allows ions to pass across the membrane, and an intracellular domain (ICD) that plays a role in channel conductance and receptor modulation. Cys-loop receptors are the targets for many currently used clinically relevant drugs (e.g. benzodiazepines and anaesthetics). Understanding the molecular mechanisms of these receptors could therefore provide the catalyst for further development in this field, as well as promoting the development of experimental techniques for other areas of neuroscience.In this review, we present our current understanding of Cys-loop receptor structure and function. The ECD has been extensively studied. Research in this area has been stimulated in recent years by the publication of high-resolution structures of nACh receptors and related proteins, which have permitted the creation of many Cys loop receptor homology models of this region. Here, using the 5-HT3 receptor as a typical member of the family, we describe how homology modelling and ligand docking can provide useful but not definitive information about ligand interactions. We briefly consider some of the many Cys-loop receptors modulators. We discuss the current understanding of the structure of the TMD, and how this links to the ECD to allow channel gating, and consider the roles of the ICD, whose structure is poorly understood. We also describe some of the current methods that are beginning to reveal the differences between different receptor states, and may ultimately show structural details of transitions between them.

Theoretical study of the structures and hydrogen-bond properties of new alternated heterocyclic compounds

The conformational preferences of a new bis-pyrrole derivative and its bis-pyridazine precursor have been investigated through quantum chemistry calculations (HF, DFT(MPWB1K), LMP2) and observations in the solid state. The global energetic minima are planar for both structures, with the conformational preferences being explained by pi-electronic conjugation between the aromatic systems and the occurrence of intramolecular hydrogen bonds (HB). For the bis-pyridazine derivative, the all-anti preferred conformation results from CH...Nsp(2) HB whereas the all-syn conformation of the bis-pyrrole is partly due to NH...Nsp(2) HB. For both systems, the validity of the theoretical conformational features is confirmed through the excellent agreement with the experimental data available. Calculations of electrostatic potential computed on the molecular surface of the various structures and their building blocks allow the variations to be rationalized in terms of molecular structure and are used to analyze the HB donor and acceptor sites of the compounds. The HB interaction sites predicted from the quantum chemical calculations are confirmed through the HB interactions observed in relevant crystal structures.

Nicotinic pharmacophore: the pyridine N of nicotine and carbonyl of acetylcholine hydrogen bond across a subunit interface to a backbone NH

Pharmacophore models for nicotinic agonists have been proposed for four decades. Central to these models is the presence of a cationic nitrogen and a hydrogen bond acceptor. It is now well-established that the cationic center makes an important cation-pi interaction to a conserved tryptophan, but the donor to the proposed hydrogen bond acceptor has been more challenging to identify. A structure of nicotine bound to the acetylcholine binding protein predicted that the binding partner of the pharmacophore's second component was a water molecule, which also hydrogen bonds to the backbone of the complementary subunit of the receptors. Here we use unnatural amino acid mutagenesis coupled with agonist analogs to examine whether such a hydrogen bond is functionally significant in the alpha4beta2 neuronal nAChR, the receptor most associated with nicotine addiction. We find evidence for the hydrogen bond with the agonists nicotine, acetylcholine, carbamylcholine, and epibatidine. These data represent a completed nicotinic pharmacophore and offer insight into the design of new therapeutic agents that selectively target these receptors.

Chemical scale studies of the Phe-Pro conserved motif in the cys loop of Cys loop receptors

The functions of two conserved residues, Phe(135) and Pro(136), located at the apex of the Cys loop of the nicotinic acetylcholine receptor are investigated. Both residues were substituted with natural and unnatural amino acids, focusing on the role of aromaticity at Phe(135), backbone conformation at Pro(136), side chain polarity and volume, and the specific interaction between the aromatic side chain and the proline. NMR spectroscopy studies of model peptides containing proline and unnatural proline analogues following a Phe show a consistent increase in the population of the cis conformer relative to peptides lacking the Phe. In the receptor, a strong interaction between the Phe and Pro residues is evident, as is a strong preference for aromaticity and hydrophobicity at the Phe site. A similar influence of hydrophobicity is observed at the proline site. In addition, across a simple homologous series of proline analogues, the results reveal a correlation between receptor function and cis bias at the proline backbone. This could suggest a significant role for the cis proline conformer at this site in receptor function.

Locating GABA in GABA receptor binding sites

The Cys-loop family of ligand-gated ion channels contains both vertebrate and invertebrate members that are activated by GABA (gamma-aminobutyric acid). Many of the residues that are critical for ligand binding have been identified in vertebrate GABA(A) and GABA(C) receptors, and specific interactions between GABA and some of these residues have been determined. In the present paper, I show how a cation-pi interaction for one of the binding site residues has allowed the production of models of GABA docked into the binding site, and these orientations are supported by mutagenesis and functional data. Surprisingly, however, the residue that forms the cation-pi interaction is not conserved, suggesting that GABA occupies subtly different locations even in such closely related receptors.

The effect of substituent rotation on aromatic quadrupole moments

We recently reported substituent constants for the accurate prediction of molecular quadrupole moments of mono-, di-, tri- and tetra-substituted aromatics. Four of the substituents in the study, –OH, –NO2, –NH2, and –CH3, were polyatomic and for these groups the substituent constants only hold for the lowest energy, or near-lowest energy, geometries. Herein, we report a computational investigation of the effect of rotation of –OH, –NO2, –NH2, and –CH3 groups on the aromatic quadrupole moment, Θzz. As expected, rotation of these substituents significantly affects the aromatic Θzz value; however, the affects are clearly periodic. Additionally, we have modified the methods to best employ our substituent constants.

Insights into docking and scoring neuronal alpha4beta2 nicotinic receptor agonists using molecular dynamics simulations and QM/MM calculations

A combined quantum mechanical (QM)-polarized docking and molecular dynamics approach to study the binding mode and to predict the binding affinity of ligands acting at the alpha4beta2-nAChR is presented. The results obtained in this study indicate that the quantum mechanical/molecular mechanics docking protocol well describes the charge-driven interactions occurring in the binding of nicotinic agonists, and it is able to represent the polarization effects on the ligand exerted by the surrounding atoms of the receptor at the binding site. This makes it possible to properly score agonists of alpha4beta2-nAChR and to reproduce the experimental binding affinity data with good accuracy, within a mean error of 2.2 kcal/mol. Moreover, applying the QM-polarized docking to an ensemble of nAChR conformations obtained from MD simulations enabled us to accurately capture nAChR-ligand induced-fit effects.