Pentalenene formation mechanisms redux.

Quantum chemical calculations are used to assess various means of lowering the barrier for the dyotropic rearrangement previously proposed to occur during the carbocation rearrangement process promoted by pentalenene synthase. Several means of lowering this barrier, including a stepwise pathway for dyotropic rearrangement, are uncovered.

Asymmetric ligand binding facilitates conformational transitions in pentameric ligand-gated ion channels.

The anesthetic propofol inhibits the currents of the homopentameric ligand-gated ion channel GLIC, yet the crystal structure of GLIC with five propofol molecules bound symmetrically shows an open-channel conformation. To address this dilemma and determine if the symmetry of propofol binding sites affects the channel conformational transition, we performed a total of 1.5 µs of molecular dynamics simulations for different GLIC systems with propofol occupancies of 0, 1, 2, 3, and 5. GLIC without propofol binding or with five propofol molecules bound symmetrically, showed similar channel conformation and hydration status over multiple replicates of 100-ns simulations. In contrast, asymmetric binding to one, two or three equivalent sites in different subunits accelerated the channel dehydration, increased the conformational heterogeneity of the pore-lining TM2 helices, and shifted the lateral and radial tilting angles of TM2 toward a closed-channel conformation. The results differentiate two groups of systems based on the propofol binding symmetry. The difference between symmetric and asymmetric groups is correlated with the variance in the propofol-binding cavity adjacent to the hydrophobic gate and the force imposed by the bound propofol. Asymmetrically bound propofol produced greater variance in the cavity size that could further elevate the conformation heterogeneity. The force trajectory generated by propofol in each subunit over the course of a simulation exhibits an ellipsoidal shape, which has the larger component tangential to the pore. Asymmetric propofol binding creates an unbalanced force that expedites the channel conformation transitions. The findings from this study not only suggest that asymmetric binding underlies the propofol functional inhibition of GLIC, but also advocate for the role of symmetry breaking in facilitating channel conformational transitions.

Isoflurane alters the structure and dynamics of GLIC.

Pentameric ligand-gated ion channels are targets of general anesthetics. Although the search for discrete anesthetic binding sites has achieved some degree of success, little is known regarding how anesthetics work after the events of binding. Using the crystal structures of the bacterial Gloeobacter violaceus pentameric ligand-gated ion channel (GLIC), which is sensitive to a variety of general anesthetics, we performed multiple molecular dynamics simulations in the presence and absence of the general anesthetic isoflurane. Isoflurane bound to several locations within GLIC, including the transmembrane pocket identified crystallographically, the extracellular (EC) domain, and the interface of the EC and transmembrane domains. Isoflurane also entered the channel after the pore was dehydrated in one of the simulations. Isoflurane disrupted the quaternary structure of GLIC, as evidenced in a striking association between the binding and breakage of intersubunit salt bridges in the EC domain. The pore-lining helix experienced lateral and inward radial tilting motion that contributed to the channel closure. Isoflurane binding introduced strong anticorrelated motions between different subunits of GLIC. The demonstrated structural and dynamical modulations by isoflurane aid in the understanding of the underlying mechanism of anesthetic inhibition of GLIC and possibly other homologous pentameric ligand-gated ion channels.

The role of structured water in mediating general anesthetic action on alpha4beta2 nAChR.

Water is an essential component for many biological processes. Pauling proposed that water might play a critical role in general anesthesia by forming water clathrates around anesthetic molecules. To examine potential involvement of water in general anesthesia, we analyzed water within alpha4beta2 nAChR, a putative protein target hypersensitive to volatile anesthetics. Experimental structure-derived closed- and open-channel nAChR systems in a fully hydrated lipid bilayer were examined using all-atom molecular dynamics simulations. At the majority of binding sites in alpha4beta2 nAChR, halothane replaced the slow-exchanging water molecules and caused a regional water population decrease. Only two binding sites had an increased quantity of water in the presence of halothane, where water arrangements resemble clathrate-like structures. The small number of such clathrate-like water clusters suggests that the formation of water clathrates is unlikely to be a primary cause for anesthesia. Despite the decrease in water population at most of the halothane binding sites, the number of sites that can be occupied transiently by water is actually increased in the presence of halothane. Many of these water sites were located between two subunits or in regions containing agonist binding sites or critical structural elements for transducing agonist binding to channel gating. Changes in water sites in the presence of halothane affected water-mediated protein-protein interactions and the protein dynamics, which can have direct impact on protein function. Collectively, water contributes to the action of anesthetics in proteins by mediating interactions between protein subunits and altering protein dynamics, instead of forming water clathrates around anesthetics.

General anesthetic binding to neuronal alpha4beta2 nicotinic acetylcholine receptor and its effects on global dynamics.

The neuronal alpha4beta2 nicotinic acetylcholine receptor (nAChR) is a target for general anesthetics. Currently available experimental structural information is inadequate to understand where anesthetics bind and how they modulate the receptor motions essential to function. Using our published open-channel structure model of alpha4beta2 nAChR, we identified and evaluated six amphiphilic interaction sites for the volatile anesthetic halothane via flexible ligand docking and subsequent 20-ns molecular dynamics simulations. Halothane binding energies ranged from -6.8 to -2.4 kcal/mol. The primary binding sites were located at the interface of extracellular and transmembrane domains, where halothane perturbed conformations of, and widened the gap among, the Cys loop, the beta1-beta2 loop, and the TM2-TM3 linker. The halothane with the highest binding affinity at the interface between the alpha4 and beta2 subunits altered interactions between the protein and nearby lipids by competing for hydrogen bonds. Gaussian network model analyses of the alpha4beta2 nAChR structures at the end of 20-ns simulations in the absence or presence of halothane revealed profound changes in protein residue mobility. The concerted motions critical to protein function were also perturbed considerably. Halothane's effect on protein dynamics was not confined to the residues adjacent to the binding sites, indicating an action on a more global scale.

Design and synthesis of propeller-shaped dispiroisoxazolinopiperidinochromanones.

A library of novel, propeller-shaped dispirotriheterocyclic isoxazolinopiperidinochromanones is reported. Each rigid dispirotriheterocycle was prepared in five linear steps from commercially available tert-butyl 4-oxopiperidine-1-carboxylate and various derivatives of 1-(2-hydroxyphenyl)ethanone, benzaldehyde oxime, and carboxylic acids. Computational chemistry was employed to analyze the three-dimensional geometries of these dispirotriheterocycles, as well as to generate chemoinformatic bioavailability data. X-ray crystallographic structure determination verified the regioselectivity of the nitrile oxide 1,3-dipolar cycloaddition reaction. The resulting library of compounds has been added to the National Institutes of Health repository (approximately 10 mg of each with > or =90% purity) for pilot-scale biomedical studies with bioassay data available at the National Center for Biotechnology Information PubChem database.

Structure of the pentameric ligand-gated ion channel ELIC cocrystallized with its competitive antagonist acetylcholine.

ELIC, the pentameric ligand-gated ion channel from Erwinia chrysanthemi, is a prototype for Cys-loop receptors. Here we show that acetylcholine is a competitive antagonist for ELIC. We determine the acetylcholine-ELIC cocrystal structure to a 2.9-Å resolution and find that acetylcholine binding to an aromatic cage at the subunit interface induces a significant contraction of loop C and other structural rearrangements in the extracellular domain. The side chain of the pore-lining residue F247 reorients and the pore size consequently enlarges, but the channel remains closed. We attribute the inability of acetylcholine to activate ELIC primarily to weak cation-π and electrostatic interactions in the pocket, because an acetylcholine derivative with a simple quaternary-to-tertiary ammonium substitution activates the channel. This study presents a compelling case for understanding the structural underpinning of the functional relationship between agonism and competitive antagonism in the Cys-loop receptors, providing a new framework for developing novel therapeutic drugs.

Structure of the pentameric ligand-gated ion channel GLIC bound with anesthetic ketamine.

Pentameric ligand-gated ion channels (pLGICs) are targets of general anesthetics, but a structural understanding of anesthetic action on pLGICs remains elusive. GLIC, a prokaryotic pLGIC, can be inhibited by anesthetics, including ketamine. The ketamine concentration leading to half-maximal inhibition of GLIC (58 µM) is comparable to that on neuronal nicotinic acetylcholine receptors. A 2.99 Å resolution X-ray structure of GLIC bound with ketamine revealed ketamine binding to an intersubunit cavity that partially overlaps with the homologous antagonist-binding site in pLGICs. The functional relevance of the identified ketamine site was highlighted by profound changes in GLIC activation upon cysteine substitution of the cavity-lining residue N152. The relevance is also evidenced by changes in ketamine inhibition upon the subsequent chemical labeling of N152C. The results provide structural insight into the molecular recognition of ketamine and are valuable for understanding the actions of anesthetics and other allosteric modulators on pLGICs.

Cannabinoids suppress inflammatory and neuropathic pain by targeting α3 glycine receptors.

Certain types of nonpsychoactive cannabinoids can potentiate glycine receptors (GlyRs), an important target for nociceptive regulation at the spinal level. However, little is known about the potential and mechanism of glycinergic cannabinoids for chronic pain treatment. We report that systemic and intrathecal administration of cannabidiol (CBD), a major nonpsychoactive component of marijuana, and its modified derivatives significantly suppress chronic inflammatory and neuropathic pain without causing apparent analgesic tolerance in rodents. The cannabinoids significantly potentiate glycine currents in dorsal horn neurons in rat spinal cord slices. The analgesic potency of 11 structurally similar cannabinoids is positively correlated with cannabinoid potentiation of the α3 GlyRs. In contrast, the cannabinoid analgesia is neither correlated with their binding affinity for CB1 and CB2 receptors nor with their psychoactive side effects. NMR analysis reveals a direct interaction between CBD and S296 in the third transmembrane domain of purified α3 GlyR. The cannabinoid-induced analgesic effect is absent in mice lacking the α3 GlyRs. Our findings suggest that the α3 GlyRs mediate glycinergic cannabinoid-induced suppression of chronic pain. These cannabinoids may represent a novel class of therapeutic agents for the treatment of chronic pain and other diseases involving GlyR dysfunction.

Higher susceptibility to halothane modulation in open- than in closed-channel alpha4beta2 nAChR revealed by molecular dynamics simulations.

The neuronal alpha4beta2 nicotinic acetylcholine receptor (nAChR) is a potential molecular target for general anesthetics. It is unclear, however, whether anesthetic action produces the same effect on the open and closed channels. Computations parallel to our previous open channel study (J. Phys. Chem. B 2009, 113, 12581) were performed on the closed-channel alpha4beta2 nAChR to investigate the conformation-dependent anesthetic effects on channel structures and dynamics. Flexible ligand docking and over 20 ns molecular dynamics simulations revealed similar halothane-binding sites in the closed and open channels. The sites with relatively high binding affinities (approximately -6.0 kcal/mol) were identified at the interface of extracellular (EC) and transmembrane (TM) domains or at the interface between alpha4 and beta2 subunits. Despite similar sites for halothane binding, the closed-channel conformation showed much less sensitivity than the open channel to the structural and dynamical perturbations from halothane. Compared to the systems without anesthetics, the amount of water inside the pore decreased by 22% in the presence of halothane in the open channel but only by 6% in the closed channel. Comparison of the nonbonded interactions at the EC/TM interfaces suggested that the beta2 subunits were more prone than the alpha4 subunits to halothane binding. In addition, our data support the notion that halothane exerts its effect by disturbing the quaternary structure and dynamics of the channel. The study concludes that sensitivity and global dynamics responsiveness of alpha4beta2 nAChR to halothane are conformation dependent. The effect of halothane on the global dynamics of the open-channel conformation might also account for the action of other inhaled general anesthetics.

Unresponsive correlated motion in alpha7 nAChR to halothane binding explains its functional insensitivity to volatile anesthetics.

Neuronal nicotinic acetylcholine receptors (nAChRs) have been implicated as targets for general anesthetics, but the functional responses to anesthetic modulation vary considerably among different subtypes of nAChRs. Inhaled general anesthetics, such as halothane, could effectively inhibit the channel activity of the alpha4beta2 nAChR but not the homologous alpha7 nAChR. To understand why alpha7 is insensitive to inhaled general anesthetics, we performed multiple sets of 20 ns molecular dynamics (MD) simulations on the closed- and open-channel alpha7 in the absence and presence of halothane and critically compared the results with those from our studies on the alpha4beta2 nAChR (Liu et al. J. Phys. Chem. B 2009, 113, 12581 and Liu et al. J. Phys. Chem. B 2010, 114, 626). Several halothane binding sites with fairly high binding affinities were identified in both closed- and open-channel alpha7, suggesting that a lack of sensitive functional responses of the alpha7 nAChR to halothane in the previous experiments was unlikely due to a lack of halothane interaction with alpha7. The binding affinities of halothane in alpha7 seemed to be protein conformation-dependent. Overall, halothane affinity was higher in the closed-channel alpha7. Halothane binding to alpha7 did not induce profound changes in alpha7 structure and dynamics that could be related to the channel function. In contrast, correlated motion of the open-channel alpha4beta2 was reduced substantially in the presence of halothane, primarily due to the more susceptible nature of beta2 to anesthetic modulation. The amphiphilic extracellular and transmembrane domain interface of the beta2 subunit is attractive to halothane and is susceptible to halothane perturbation, which may be why alpha4beta2 is functionally more sensitive to halothane than alpha7.

Mechanistic studies on the stereoselective formation of beta-mannosides from mannosyl iodides using alpha-deuterium kinetic isotope effects.

Stereoselective synthesis of beta-mannosides is one of the most challenging linkages to achieve in carbohydrate chemistry. Both the anomeric effect and the C2 axial substituent favor the formation of the axial glycoside (alpha-product). Herein, we describe mechanistic studies on the beta-selective glycosidation of trimethylene oxide (TMO) using mannosyl iodides. Density functional calculations (at the B3LYP/6-31+G(d,p):LANL2DZ level) suggest that formation of both alpha- and beta-mannosides involve loose S(N)2-like transition-state structures with significant oxacarbenium character, although the transition structure for formation of the alpha-mannoside is significantly looser. alpha-Deuterium kinetic isotope effects (alpha-DKIEs) based upon these computed transition state geometries match reasonably well with the experimentally measured values: 1.16 +/- 0.02 for the beta-linkage (computed to be 1.15) and 1.19 +/- 0.05, see table 2 for the alpha-analogue (computed to be 1.26). Since it was unclear if beta-selectivity resulted from a conformational constraint induced by the anomeric iodide, a 4,6-O-benzylidine acetal was used to lock the iodide into a chairlike conformation. Both experiments and calculations on this analogue suggest that it does not mirror the behavior of mannosyl iodides lacking bridging acetal protecting groups.

Phenylglycine and sulfonamide correctors of defective delta F508 and G551D cystic fibrosis transmembrane conductance regulator chloride-channel gating.

Mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) chloride channel cause cystic fibrosis. The delta F508 mutation produces defects in channel gating and cellular processing, whereas the G551D mutation produces primarily a gating defect. To identify correctors of gating, 50,000 diverse small molecules were screened at 2.5 microM (with forskolin, 20 microM) by an iodide uptake assay in epithelial cells coexpressing delta F508-CFTR and a fluorescent halide indicator (yellow fluorescent protein-H148Q/I152L) after delta F508-CFTR rescue by 24-h culture at 27 degrees C. Secondary analysis and testing of >1000 structural analogs yielded two novel classes of correctors of defective delta F508-CFTR gating ("potentiators") with nanomolar potency that were active in human delta F508 and G551D cells. The most potent compound of the phenylglycine class, 2-[(2-1H-indol-3-yl-acetyl)-methylamino]-N-(4-isopropylphenyl)-2-phenylacetamide, reversibly activated delta F508-CFTR in the presence of forskolin with K(a) approximately 70 nM and also activated the CFTR gating mutants G551D and G1349D with K(a) values of approximately 1100 and 40 nM, respectively. The most potent sulfonamide, 6-(ethylphenylsulfamoyl)-4-oxo-1,4-dihydroquinoline-3-carboxylic acid cycloheptylamide, had K(a) approximately 20 nM for activation of delta F508-CFTR. In cell-attached patch-clamp experiments, phenylglycine-01 (PG-01) and sulfonamide-01 (SF-01) increased channel open probability >5-fold by the reduction of interburst closed time. An interesting property of these compounds was their ability to act in synergy with cAMP agonists. Microsome metabolism studies and rat pharmacokinetic analysis suggested significantly more rapid metabolism of PG-01 than SF-03. Phenylglycine and sulfonamide compounds may be useful for monotherapy of cystic fibrosis caused by gating mutants and possibly for a subset of delta F508 subjects with significant delta F508-CFTR plasma-membrane expression.