Anesthetic sensitivity of the Gloeobacter violaceus proton-gated ion channel

Multisite binding of a general anesthetic to the prokaryotic pentameric Erwinia chrysanthemi ligand-gated ion channel (ELIC)

Pentameric ligand-gated ion channels (pLGICs), such as nicotinic acetylcholine, glycine, γ-aminobutyric acid GABA_A/C receptors, and the Gloeobacter violaceus ligand-gated ion channel (GLIC), are receptors that contain multiple allosteric binding sites for a variety of therapeutics, including general anesthetics. Here, we report the x-ray crystal structure of the Erwinia chrysanthemi ligand-gated ion channel (ELIC) in complex with a derivative of chloroform, which reveals important features of anesthetic recognition, involving multiple binding at three different sites. One site is located in the channel pore and equates with a noncompetitive inhibitor site found in many pLGICs. A second transmembrane site is novel and is located in the lower part of the transmembrane domain, at an interface formed between adjacent subunits. A third site is also novel and is located in the extracellular domain in a hydrophobic pocket between the β7–β10 strands. Together, these results extend our understanding of pLGIC modulation and reveal several specific binding interactions that may contribute to modulator recognition, further substantiating a multisite model of allosteric modulation in this family of ion channels.

Multiscale modeling of solvation in chemical and biological nanosystems and in nanoporous materials

Statistical–mechanical, 3D-RISM-KH molecular theory of solvation (3D reference interaction site model with the Kovalenko–Hirata closure) is promising as an essential part of multiscale methodology for chemical and biomolecular nanosystems in solution. 3D-RISM-KH explains the molecular mechanisms of self-assembly and conformational stability of synthetic organic rosette nanotubes (RNTs), aggregation of prion proteins and β-sheet amyloid oligomers, protein-ligand binding, and function-related solvation properties of complexes as large as the Gloeobacter violaceus pentameric ligand-gated ion channel (GLIC) and GroEL/ES chaperone. Molecular mechanics/Poisson–Boltzmann (generalized Born) surface area [MM/PB(GB)SA] post-processing of molecular dynamics (MD) trajectories involving SA empirical nonpolar terms is replaced with MM/3D-RISM-KH statistical–mechanical evaluation of the solvation thermodynamics. 3D-RISM-KH has been coupled with multiple time-step (MTS) MD of the solute biomolecule driven by effective solvation forces, which are obtained analytically by converging the 3D-RISM-KH integral equations at outer time-steps and are calculated in between by using solvation force coordinate extrapolation (SFCE) in the subspace of previous solutions to 3D-RISM-KH. The procedure is stabilized by the optimized isokinetic Nosé–Hoover (OIN) chain thermostatting, which enables gigantic outer time-steps up to picoseconds to accurately calculate equilibrium properties. The multiscale OIN/SFCE/3D-RISM-KH algorithm is implemented in the Amber package and illustrated on a fully flexible model of alanine dipeptide in aqueous solution, exhibiting the computational rate of solvent sampling 20 times faster than standard MD with explicit solvent. Further substantial acceleration can be achieved with 3D-RISM-KH efficiently sampling essential events with rare statistics such as exchange and localization of solvent, ions, and ligands at binding sites and pockets of the biomolecule. 3D-RISM-KH was coupled with ab initio complete active space self-consistent field (CASSCF) and orbital-free embedding (OFE) Kohn–Sham (KS) density functional theory (DFT) quantum chemistry methods in an SCF description of electronic structure, optimized geometry, and chemical reactions in solution. The (OFE)KS-DFT/3D-RISM-KH multi-scale method is implemented in the Amsterdam Density Functional (ADF) package and extensively validated against experiment for solvation thermochemistry, photochemistry, conformational equilibria, and activation barriers of various nanosystems in solvents and ionic liquids (ILs). Finally, the replica RISM-KH-VM molecular theory for the solvation structure, thermodynamics, and electrochemistry of electrolyte solutions sorbed in nanoporous materials reveals the molecular mechanisms of sorption and supercapacitance in nanoporous carbon electrodes, which is drastically different from a planar electrical double layer.

Reversal of ion-charge selectivity renders the pentameric ligand-gated ion channel GLIC insensitive to anaesthetics

pLGICs (pentameric ligand-gated ion channels) are a family of structurally homologous cation and anion channels involved in neurotransmission. Cation-selective members of the pLGIC family are typically inhibited by general anaesthetics, whereas anion-selective members are potentiated. GLIC is a prokaryotic cation pLGIC and can be inhibited by clinical concentrations of general anaesthetics. The introduction of three mutations, Y221A (Y-3'A), E222P (E-2'P) and N224R (N0'R), at the selectivity filter and one, A237T (A13'T), at the hydrophobic gate, converted GLIC into an anion channel. The mutated GLIC (GLIC4) became insensitive to the anaesthetics propofol and etomidate, as well as the channel blocker picrotoxin. MD (molecular dynamics) simulations revealed changes in the structure and dynamics of GLIC4 in comparison with GLIC, particularly in the tilting angles of the pore-lining helix [TM2 (transmembrane helix 2)] that consequently resulted in different pore radius and hydration profiles. Propofol binding to an intra-subunit site of GLIC shifted the tilting angles of TM2 towards closure at the hydrophobic gate region, consistent with propofol inhibition of GLIC. In contrast, the pore of GLIC4 was much more resilient to perturbation from propofol binding. The present study underscores the importance of pore dynamics and conformation to anaesthetic effects on channel functions.

Modeling complex biological systems: From solution chemistry to membranes and channels

Complex biological systems are intimately linked to their environment, a very crowded and equally complex solution compartmentalized by fluid membranes. Modeling such systems remains challenging and requires a suitable representation of these solutions and their interfaces. Here, we focus on particle-based modeling at an atomistic level using molecular dynamics (MD) simulations. As an example, we discuss important steps in modeling the solution chemistry of an ion channel of the ligand-gated ion channel receptor family, a major target of many drugs including anesthetics and addiction treatments. The bacterial pentameric ligand-gated ion channel (pLGIC) called GLIC provides clues about the functional importance of solvation, in particular for mechanisms such as permeation and gating. We present some current challenges along with promising novel modeling approaches.

Conscious processing: implications for general anesthesia

PURPOSE OF REVIEW

Can a general anesthetic binding to its membrane receptors alter global brain activity to cause loss of consciousness?

RECENT FINDINGS

The identification of the neurobiological mechanisms of conscious processing in awake patients that are altered by general anesthetics and the atomic structure of the general anesthetic's binding site.

SUMMARY

An important feature of general anesthesia is a preferential inhibition of global feedback connectivity when general anesthetics bind to allosteric sites of gamma-aminobutyric acid receptors present in the cerebral cortex.

Structure and pharmacology of pentameric receptor channels: from bacteria to brain

Orthologs of the pentameric receptor channels that mediate fast synaptic transmission in the central and peripheral nervous systems have been found in several bacterial species and in a single archaea genus. Recent X-ray structures of bacterial and invertebrate pentameric receptors point to a striking conservation of the structural features within the whole family, even between distant prokaryotic and eukaryotic members. These structural data reveal general principles of molecular organization that allow allosteric membrane proteins to mediate chemoelectric transduction. Notably, several conformations have been solved, including open and closed channels with distinct global tertiary and quaternary structure. The data reveal features of the ion channel architecture and of diverse categories of binding sites, such as those that bind orthosteric ligands, including neurotransmitters, and those that bind allosteric modulators, such as general anesthetics, ivermectin, or lipids. In this review, we summarize the most recent data, discuss insights into the mechanism of action in these systems, and elaborate on newly opened avenues for drug design.

Molecular mechanism for the dual alcohol modulation of Cys-loop receptors

Cys-loop receptors constitute a superfamily of pentameric ligand-gated ion channels (pLGICs), including receptors for acetylcholine, serotonin, glycine and γ-aminobutyric acid. Several bacterial homologues have been identified that are excellent models for understanding allosteric binding of alcohols and anesthetics in human Cys-loop receptors. Recently, we showed that a single point mutation on a prokaryotic homologue (GLIC) could transform it from a channel weakly potentiated by ethanol into a highly ethanol-sensitive channel. Here, we have employed molecular simulations to study ethanol binding to GLIC, and to elucidate the role of the ethanol-enhancing mutation in GLIC modulation. By performing 1-µs simulations with and without ethanol on wild-type and mutated GLIC, we observed spontaneous binding in both intra-subunit and inter-subunit transmembrane cavities. In contrast to the glycine receptor GlyR, in which we previously observed ethanol binding primarily in an inter-subunit cavity, ethanol primarily occupied an intra-subunit cavity in wild-type GLIC. However, the highly ethanol-sensitive GLIC mutation significantly enhanced ethanol binding in the inter-subunit cavity. These results demonstrate dramatic effects of the F(14′)A mutation on the distribution of ligands, and are consistent with a two-site model of pLGIC inhibition and potentiation.

Communication from one nerve cell to the next is an essential process for brain and muscle function. Nerve impulses result in release of transmitter molecules from one cell that bind to receptors on the next cell. Transmitter binding opens a pore in each receptor and ions flow across the membrane, leading to either enhancement or inhibition of new nerve impulses. These receptors are modulated by numerous drugs, including alcohols and anesthetics; identifying the precise location of modulator binding is critical for drug development. We have used computer simulation methods to model alcohol diffusion and binding to a receptor. By modifying a single residue in the receptor, we were able to move the location of the binding site and dramatically alter alcohol modulation, which supports a model with two separate binding sites for enhancement and inhibition in this family of receptors.

Conformational transitions underlying pore opening and desensitization in membrane-embedded Gloeobacter violaceus ligand-gated ion channel (GLIC)

Direct structural insight into the mechanisms underlying activation and desensitization remain unavailable for the pentameric ligand-gated channel family. Here, we report the structural rearrangements underlying gating transitions in membrane-embedded GLIC, a prokaryotic homologue, using site-directed spin labeling and electron paramagnetic resonance (EPR) spectroscopy. We particularly probed the conformation of pore-lining second transmembrane segment (M2) under conditions that favor the closed and the ligand-bound desensitized states. The spin label mobility, intersubunit spin-spin proximity, and the solvent-accessibility parameters in the two states clearly delineate the underlying protein motions within M2. Our results show that during activation the extracellular hydrophobic region undergoes major changes involving an outward translational movement, away from the pore axis, leading to an increase in the pore diameter, whereas the lower end of M2 remains relatively immobile. Most notably, during desensitization, the intervening polar residues in the middle of M2 move closer to form a solvent-occluded barrier and thereby reveal the location of a distinct desensitization gate. In comparison with the crystal structure of GLIC, the structural dynamics of the channel in a membrane environment suggest a more loosely packed conformation with water-accessible intrasubunit vestibules penetrating from the extracellular end all the way to the middle of M2 in the closed state. These regions have been implicated to play a major role in alcohol and drug modulation. Overall, these findings represent a key step toward understanding the fundamentals of gating mechanisms in this class of channels.

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.

Diffusive nature of xenon anesthetic changes properties of a lipid bilayer: molecular dynamics simulations

Effects of general anesthesia can be controllable by the ambient pressure. We perform molecular dynamics simulations for a 1-palmitoyl-2-oleoyl phosphatidylethanolamine lipid bilayer with or without xenon molecules by changing the pressure to elucidate the mechanism of the pressure reversal of general anesthesia. According to the diffusive nature of xenon molecules in the lipid bilayer, a decrease in the orientational order of the lipid tails, an increase in the area and volume per lipid molecule, and an increase in the diffusivity of lipid molecules are observed. We show that the properties of the lipid bilayer with xenon molecules at high pressure come close to those without xenon molecules at 0.1 MPa. Furthermore, we find that xenon molecules are concentrated in the middle of the lipid bilayer at high pressures by the pushing effect and that the diffusivity of xenon molecules is suppressed. These results suggest that the pressure reversal originates from a jamming and suppression of the diffusivity of xenon molecules in lipid bilayers.

General anesthetics predicted to block the GLIC pore with micromolar affinity

Although general anesthetics are known to modulate the activity of ligand-gated ion channels in the Cys-loop superfamily, there is at present neither consensus on the underlying mechanisms, nor predictive models of this modulation. Viable models need to offer quantitative assessment of the relative importance of several identified anesthetic binding sites. However, to date, precise affinity data for individual sites has been challenging to obtain by biophysical means. Here, the likely role of pore block inhibition by the general anesthetics isoflurane and propofol of the prokaryotic pentameric channel GLIC is investigated by molecular simulations. Microscopic affinities are calculated for both single and double occupancy binding of isoflurane and propofol to the GLIC pore. Computations are carried out for an open-pore conformation in which the pore is restrained to crystallographic radius, and a closed-pore conformation that results from unrestrained molecular dynamics equilibration of the structure. The GLIC pore is predicted to be blocked at the micromolar concentrations for which inhibition by isofluorane and propofol is observed experimentally. Calculated affinities suggest that pore block by propofol occurs at signifcantly lower concentrations than those for which inhibition is observed: we argue that this discrepancy may result from binding of propofol to an allosteric site recently identified by X-ray crystallography, which may cause a competing gain-of-function effect. Affinities of isoflurane and propofol to the allosteric site are also calculated, and shown to be 3 mM for isoflurane and for propofol; both anesthetics have a lower affinity for the allosteric site than for the unoccupied pore.

Although general anesthesia is performed every day on thousands of people, its detailed microscopic mechanisms are not known. What is known is that general anesthetic drugs modulate the activity of ion channels in the central nervous system. These channels are proteins that open in response to binding of neurotransmitter molecules, creating an electric current through the cell membrane and thus propagating nerve impulses between cells. One possible mechanism for ion channel inhibition by anesthetics is that the drugs bind inside the pore of the channels, blocking ion current. Here we investigate such a pore block mechanism by computing the strength of the drugs' interaction with the pore – and hence the likelihood of binding, in the case of GLIC, a bacterial channel protein. The results, obtained from numerical simulations of atomic models of GLIC, indicate that the anesthetics isoflurane and propofol have a tendency to bind in the pore that is strong enough to explain blocking of the channel, even at low concentration of the drugs.

Desensitization mechanism in prokaryotic ligand-gated ion channel

Crystal structures of Gloeobacter violaceus ligand-gated ion channel (GLIC), a proton-gated prokaryotic homologue of pentameric ligand-gated ion channel (LGIC) from G. violaceus, have provided high-resolution models of the channel architecture and its role in selective ion conduction and drug binding. However, it is still unclear which functional states of the LGIC gating scheme these crystal structures represent. Much of this uncertainty arises from a lack of thorough understanding of the functional properties of these prokaryotic channels. To elucidate the molecular events that constitute gating, we have carried out an extensive characterization of GLIC function and dynamics in reconstituted proteoliposomes by patch clamp measurements and EPR spectroscopy. We find that GLIC channels show rapid activation upon jumps to acidic pH followed by a time-dependent loss of conductance because of desensitization. GLIC desensitization is strongly coupled to activation and is modulated by voltage, permeant ions, pore-blocking drugs, and membrane cholesterol. Many of these properties are parallel to functions observed in members of eukaryotic LGIC. Conformational changes in loop C, measured by site-directed spin labeling and EPR spectroscopy, reveal immobilization during desensitization analogous to changes in LGIC and acetylcholine binding protein. Together, our studies suggest conservation of mechanistic aspects of desensitization among LGICs of prokaryotic and eukaryotic origin.

Binding site and affinity prediction of general anesthetics to protein targets using docking

BACKGROUND

The protein targets for general anesthetics remain unclear. A tool to predict anesthetic binding for potential binding targets is needed. In this study, we explored whether a computational method, AutoDock, could serve as such a tool.

METHODS

High-resolution crystal data of water-soluble proteins (cytochrome C, apoferritin, and human serum albumin), and a membrane protein (a pentameric ligand-gated ion channel from Gloeobacter violaceus [GLIC]) were used. Isothermal titration calorimetry (ITC) experiments were performed to determine anesthetic affinity in solution conditions for apoferritin. Docking calculations were performed using DockingServer with the Lamarckian genetic algorithm and the Solis and Wets local search method (http://www.dockingserver.com/web). Twenty general anesthetics were docked into apoferritin. The predicted binding constants were compared with those obtained from ITC experiments for potential correlations. In the case of apoferritin, details of the binding site and their interactions were compared with recent cocrystallization data. Docking calculations for 6 general anesthetics currently used in clinical settings (isoflurane, sevoflurane, desflurane, halothane, propofol, and etomidate) with known 50% effective concentration (EC(50)) values were also performed in all tested proteins. The binding constants derived from docking experiments were compared with known EC(50) values and octanol/water partition coefficients for the 6 general anesthetics.

RESULTS

All 20 general anesthetics docked unambiguously into the anesthetic binding site identified in the crystal structure of apoferritin. The binding constants for 20 anesthetics obtained from the docking calculations correlate significantly with those obtained from ITC experiments (P = 0.04). In the case of GLIC, the identified anesthetic binding sites in the crystal structure are among the docking predicted binding sites, but not the top ranked site. Docking calculations suggest a most probable binding site located in the extracellular domain of GLIC. The predicted affinities correlated significantly with the known EC(50) values for the 6 frequently used anesthetics in GLIC for the site identified in the experimental crystal data (P = 0.006). However, predicted affinities in apoferritin, human serum albumin, and cytochrome C did not correlate with these 6 anesthetics' known experimental EC(50) values. A weak correlation between the predicted affinities and the octanol/water partition coefficients was observed for the sites in GLIC.

CONCLUSION

We demonstrated that anesthetic binding sites and relative affinities can be predicted using docking calculations in an automatic docking server (AutoDock) for both water-soluble and membrane proteins. Correlation of predicted affinity and EC(50) for 6 frequently used general anesthetics was only observed in GLIC, a member of a protein family relevant to anesthetic mechanism.

Recognition of anesthetic barbiturates by a protein binding site: a high resolution structural analysis

Barbiturates potentiate GABA actions at the GABA_A receptor and act as central nervous system depressants that can induce effects ranging from sedation to general anesthesia. No structural information has been available about how barbiturates are recognized by their protein targets. For this reason, we tested whether these drugs were able to bind specifically to horse spleen apoferritin, a model protein that has previously been shown to bind many anesthetic agents with affinities that are closely correlated with anesthetic potency. Thiopental, pentobarbital, and phenobarbital were all found to bind to apoferritin with affinities ranging from 10–500 µM, approximately matching the concentrations required to produce anesthetic and GABAergic responses. X-ray crystal structures were determined for the complexes of apoferritin with thiopental and pentobarbital at resolutions of 1.9 and 2.0 Å, respectively. These structures reveal that the barbiturates bind to a cavity in the apoferritin shell that also binds haloalkanes, halogenated ethers, and propofol. Unlike these other general anesthetics, however, which rely entirely upon van der Waals interactions and the hydrophobic effect for recognition, the barbiturates are recognized in the apoferritin site using a mixture of both polar and nonpolar interactions. These results suggest that any protein binding site that is able to recognize and respond to the chemically and structurally diverse set of compounds used as general anesthetics is likely to include a versatile mixture of both polar and hydrophobic elements.

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.

Cys-loop receptor channel blockers also block GLIC

The Gloeobacter ligand-gated ion channel (GLIC) is a bacterial homolog of vertebrate Cys-loop ligand-gated ion channels. Its pore-lining region in particular has a high sequence homology to these related proteins. Here we use electrophysiology to examine a range of compounds that block the channels of Cys-loop receptors to probe their pharmacological similarity with GLIC. The data reveal that a number of these compounds also block GLIC, although the pharmacological profile is distinct from these other proteins. The most potent compound was lindane, a GABA(A) receptor antagonist, with an IC₅₀ of 0.2 μM. Docking studies indicated two potential binding sites for this ligand in the pore, at the 9' or between the 0' and 2' residues. Similar experiments with picrotoxinin (IC₅₀ = 2.6 μM) and rimantadine (IC₅₀ = 2.6 μM) reveal interactions with 2'Thr residues in the GLIC pore. These locations are strongly supported by mutagenesis data for picrotoxinin and lindane, which are less potent in a T2'S version of GLIC. Overall, our data show that the inhibitory profile of the GLIC pore has considerable overlap with those of Cys-loop receptors, but the GLIC pore has a unique pharmacology.

Expression and characterization of the bacterial mechanosensitive channel MscS in Xenopus laevis oocytes

We have successfully expressed and characterized mechanosensitive channel of small conductance (MscS) from Escherichia coli in oocytes of the African clawed frog, Xenopus laevis. MscS expressed in oocytes has the same single-channel conductance and voltage dependence as the channel in its native environment. Two hallmarks of MscS activity, the presence of conducting substates at high potentials and reversible adaptation to a sustained stimulus, are also exhibited by oocyte-expressed MscS. In addition to its ease of use, the oocyte system allows the user to work with relatively large patches, which could be an advantage for the visualization of membrane deformation. Furthermore, MscS can now be compared directly to its eukaryotic homologues or to other mechanosensitive channels that are not easily studied in E. coli.

Structural basis for alcohol modulation of a pentameric ligand-gated ion channel

Despite its long history of use and abuse in human culture, the molecular basis for alcohol action in the brain is poorly understood. The recent determination of the atomic-scale structure of GLIC, a prokaryotic member of the pentameric ligand-gated ion channel (pLGIC) family, provides a unique opportunity to characterize the structural basis for modulation of these channels, many of which are alcohol targets in brain. We observed that GLIC recapitulates bimodal modulation by n-alcohols, similar to some eukaryotic pLGICs: methanol and ethanol weakly potentiated proton-activated currents in GLIC, whereas n-alcohols larger than ethanol inhibited them. Mapping of residues important to alcohol modulation of ionotropic receptors for glycine, γ-aminobutyric acid, and acetylcholine onto GLIC revealed their proximity to transmembrane cavities that may accommodate one or more alcohol molecules. Site-directed mutations in the pore-lining M2 helix allowed the identification of four residues that influence alcohol potentiation, with the direction of their effects reflecting α-helical structure. At one of the potentiation-enhancing residues, decreased side chain volume converted GLIC into a highly ethanol-sensitive channel, comparable to its eukaryotic relatives. Covalent labeling of M2 positions with an alcohol analog, a methanethiosulfonate reagent, further implicated residues at the extracellular end of the helix in alcohol binding. Molecular dynamics simulations elucidated the structural consequences of a potentiation-enhancing mutation and suggested a structural mechanism for alcohol potentiation via interaction with a transmembrane cavity previously termed the "linking tunnel." These results provide a unique structural model for independent potentiating and inhibitory interactions of n-alcohols with a pLGIC family member.

Functional prokaryotic-eukaryotic chimera from the pentameric ligand-gated ion channel family

Pentameric ligand-gated ion channels (pLGICs), which mediate chemo-electric signal transduction in animals, have been recently found in bacteria. Despite clear sequence and 3D structure homology, the phylogenetic distance between prokaryotic and eukaryotic homologs suggests significant structural divergences, especially at the interface between the extracellular (ECD) and the transmembrane (TMD) domains. To challenge this possibility, we constructed a chimera in which the ECD of the bacterial protein GLIC is fused to the TMD of the human α1 glycine receptor (α1GlyR). Electrophysiology in Xenopus oocytes shows that it functions as a proton-gated ion channel, thereby locating the proton activation site(s) of GLIC in its ECD. Patch-clamp experiments in BHK cells show that the ion channel displays an anionic selectivity with a unitary conductance identical to that of the α1GlyR. In addition, pharmacological investigations result in transmembrane allosteric modulation similar to the one observed on α1GlyR. Indeed, the clinically active drugs propofol, four volatile general anesthetics, alcohols, and ivermectin all potentiate the chimera while they inhibit GLIC. Collectively, this work shows the compatibility between GLIC and α1GlyR domains and points to conservation of the ion channel and transmembrane allosteric regulatory sites in the chimera. This provides evidence that GLIC and α1GlyR share a highly homologous 3D structure. GLIC is thus a relevant model of eukaryotic pLGICs, at least from the anionic type. In addition, the chimera is a good candidate for mass production in Escherichia coli, opening the way for investigations of "druggable" eukaryotic allosteric sites by X-ray crystallography.

X-ray structures of general anaesthetics bound to a pentameric ligand-gated ion channel

General anaesthetics have enjoyed long and widespread use but their molecular mechanism of action remains poorly understood. There is good evidence that their principal targets are pentameric ligand-gated ion channels (pLGICs) such as inhibitory GABA(A) (γ-aminobutyric acid) receptors and excitatory nicotinic acetylcholine receptors, which are respectively potentiated and inhibited by general anaesthetics. The bacterial homologue from Gloeobacter violaceus (GLIC), whose X-ray structure was recently solved, is also sensitive to clinical concentrations of general anaesthetics. Here we describe the crystal structures of the complexes propofol/GLIC and desflurane/GLIC. These reveal a common general-anaesthetic binding site, which pre-exists in the apo-structure in the upper part of the transmembrane domain of each protomer. Both molecules establish van der Waals interactions with the protein; propofol binds at the entrance of the cavity whereas the smaller, more flexible, desflurane binds deeper inside. Mutations of some amino acids lining the binding site profoundly alter the ionic response of GLIC to protons, and affect its general-anaesthetic pharmacology. Molecular dynamics simulations, performed on the wild type (WT) and two GLIC mutants, highlight differences in mobility of propofol in its binding site and help to explain these effects. These data provide a novel structural framework for the design of general anaesthetics and of allosteric modulators of brain pLGICs.

Anesthetic sites and allosteric mechanisms of action on Cys-loop ligand-gated ion channels

PURPOSE

The Cys-loop ligand-gated ion channel superfamily is a major group of neurotransmitter-activated receptors in the central and peripheral nervous system. The superfamily includes inhibitory receptors stimulated by γ-aminobutyric acid (GABA) and glycine and excitatory receptors stimulated by acetylcholine and serotonin. The first part of this review presents current evidence on the location of the anesthetic binding sites on these channels and the mechanism by which binding to these sites alters their function. The second part of the review addresses the basis for this selectivity, and the third part describes the predictive power of a quantitative allosteric model showing the actions of etomidate on γ-aminobutyric acid type A receptors (GABA(A)Rs).

PRINCIPAL FINDINGS

General anesthetics at clinical concentrations inhibit the excitatory receptors and enhance the inhibitory receptors. The location of general anesthetic binding sites on these receptors is being defined by photoactivable analogues of general anesthetics. The receptor studied most extensively is the muscle-type nicotinic acetylcholine receptor (nAChR), and progress is now being made with GABA(A)Rs. There are three categories of sites that are all in the transmembrane domain: 1) within a single subunit's four-helix bundle (intrasubunit site; halothane and etomidate on the δ subunit of AChRs); 2) between five subunits in the transmembrane conduction pore (channel lumen sites; etomidate and alcohols on nAChR); and 3) between two subunits (subunit interface sites; etomidate between the α1 and β2/3 subunits of the GABA(A)R).

CONCLUSIONS

These binding sites function allosterically. Certain conformations of a receptor bind the anesthetic with greater affinity than others. Time-resolved photolabelling of some sites occurs within milliseconds of channel opening on the nAChR but not before. In GABA(A)Rs, electrophysiological data fit an allosteric model in which etomidate binds to and stabilizes the open state, increasing both the fraction of open channels and their lifetime. As predicted by the model, the channel-stabilizing action of etomidate is so strong that higher concentrations open the channel in the absence of agonist. The formal functional paradigm presented for etomidate may apply to other potent general anesthetic drugs. Combining photolabelling with structure-function mutational studies in the context of allosteric mechanisms should lead us to a more detailed understanding of how and where these important drugs act.

GABA(A) receptors as molecular targets of general anesthetics: identification of binding sites provides clues to allosteric modulation

Purpose

The purpose of this review is to summarize current knowledge of detailed biochemical evidence for the role of γ-aminobutyric acid type A receptors (GABA_A–Rs) in the mechanisms of general anesthesia.

Principal findings

With the knowledge that all general anesthetics positively modulate GABA_A-R-mediated inhibitory transmission, site-directed mutagenesis comparing sequences of GABA_A-R subunits of varying sensitivity led to identification of amino acid residues in the transmembrane domain that are critical for the drug actions in vitro. Using a photo incorporable analogue of the general anesthetic, R(+)etomidate, we identified two transmembrane amino acids that were affinity labelled in purified bovine brain GABA_A-R. Homology protein structural modelling positions these two residues, αM1-11’ and βM3-4’, close to each other in a single type of intersubunit etomidate binding pocket at the β/α interface. This position would be appropriate for modulation of agonist channel gating. Overall, available information suggests that these two etomidate binding residues are allosterically coupled to sites of action of steroids, barbiturates, volatile agents, and propofol, but not alcohols. Residue α/βM2-15’ is probably not a binding site but allosterically coupled to action of volatile agents, alcohols, and intravenous agents, and α/βM1-(-2’) is coupled to action of intravenous agents.

Conclusions

Establishment of a coherent and consistent structural model of the GABA_A-R lends support to the conclusion that general anesthetics can modulate function by binding to appropriate domains on the protein. Genetic engineering of mice with mutation in some of these GABA_A-R residues are insensitive to general anesthetics in vivo, suggesting that further analysis of these domains could lead to development of more potent and specific drugs.

Molecular dynamics and brownian dynamics investigation of ion permeation and anesthetic halothane effects on a proton-gated ion channel

Bacterial Gloeobacter violaceus pentameric ligand-gated ion channel (GLIC) is activated to cation permeation upon lowering the solution pH. Its function can be modulated by anesthetic halothane. In the present work, we integrate molecular dynamics (MD) and Brownian dynamics (BD) simulations to elucidate the ion conduction, charge selectivity, and halothane modulation mechanisms in GLIC, based on recently resolved X-ray crystal structures of the open-channel GLIC. MD calculations of the potential of mean force (PMF) for a Na(+) revealed two energy barriers in the extracellular domain (R109 and K38) and at the hydrophobic gate of transmembrane domain (I233), respectively. An energy well for Na(+) was near the intracellular entrance: the depth of this energy well was modulated strongly by the protonation state of E222. The energy barrier for Cl(-) was found to be 3-4 times higher than that for Na(+). Ion permeation characteristics were determined through BD simulations using a hybrid MD/continuum electrostatics approach to evaluate the energy profiles governing the ion movement. The resultant channel conductance and a near-zero permeability ratio (P(Cl)/P(Na)) were comparable to experimental data. On the basis of these calculations, we suggest that a ring of five E222 residues may act as an electrostatic gate. In addition, the hydrophobic gate region may play a role in charge selectivity due to a higher dehydration energy barrier for Cl(-) ions. The effect of halothane on the Na(+) PMF was also evaluated. Halothane was found to perturb salt bridges in GLIC that may be crucial for channel gating and open-channel stability, but had no significant impact on the single ion PMF profiles.

Anesthetic binding in a pentameric ligand-gated ion channel: GLIC

Cys-loop receptors are molecular targets of general anesthetics, but the knowledge of anesthetic binding to these proteins remains limited. Here we investigate anesthetic binding to the bacterial Gloeobacter violaceus pentameric ligand-gated ion channel (GLIC), a structural homolog of cys-loop receptors, using an experimental and computational hybrid approach. Tryptophan fluorescence quenching experiments showed halothane and thiopental binding at three tryptophan-associated sites in the extracellular (EC) domain, transmembrane (TM) domain, and EC-TM interface of GLIC. An additional binding site at the EC-TM interface was predicted by docking analysis and validated by quenching experiments on the N200W GLIC mutant. The binding affinities (K(D)) of 2.3 ± 0.1 mM and 0.10 ± 0.01 mM were derived from the fluorescence quenching data of halothane and thiopental, respectively. Docking these anesthetics to the original GLIC crystal structure and the structures relaxed by molecular dynamics simulations revealed intrasubunit sites for most halothane binding and intersubunit sites for thiopental binding. Tryptophans were within reach of both intra- and intersubunit binding sites. Multiple molecular dynamics simulations on GLIC in the presence of halothane at different sites suggested that anesthetic binding at the EC-TM interface disrupted the critical interactions for channel gating, altered motion of the TM23 linker, and destabilized the open-channel conformation that can lead to inhibition of GLIC channel current. The study has not only provided insights into anesthetic binding in GLIC, but also demonstrated a successful fusion of experiments and computations for understanding anesthetic actions in complex proteins.

Multiple binding sites for the general anesthetic isoflurane identified in the nicotinic acetylcholine receptor transmembrane domain

An extensive search for isoflurane binding sites in the nicotinic acetylcholine receptor (nAChR) and the proton gated ion channel from Gloebacter violaceus (GLIC) has been carried out based on molecular dynamics (MD) simulations in fully hydrated lipid membrane environments. Isoflurane introduced into the aqueous phase readily partitions into the lipid membrane and the membrane-bound protein. Specifically, isoflurane binds persistently to three classes of sites in the nAChR transmembrane domain: (i) An isoflurane dimer occludes the pore, contacting residues identified by previous mutagenesis studies; analogous behavior is observed in GLIC. (ii) Several nAChR subunit interfaces are also occupied, in a site suggested by photoaffinity labeling and thought to positively modulate the receptor; these sites are not occupied in GLIC. (iii) Isoflurane binds to the subunit centers of both nAChR alpha chains and one of the GLIC chains, in a site that has had little experimental targeting. Interpreted in the context of existing structural and physiological data, the present MD results support a multisite model for the mechanism of receptor-channel modulation by anesthetics.