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

Micro-Solvation of Propofol in Propylene Glycol-Water Binary Mixtures: Molecular Dynamics Simulation Studies

The water microstructure around propofol plays a crucial role in controlling their solubility in the binary mixture. The unusual nature of such a water microstructure can influence both translational and reorientational dynamics, as well as the water hydrogen bond network near propofol. We have carried out all-atom molecular dynamics simulations of five different compositions of the propylene glycol (PG)/water binary mixture containing propofol (PFL) molecules to investigate the differential behavior of water microsolvation shells around propofol, which is likely to control the propofol solubility. It is evident from the simulation snapshots for various compositions that the PG at high molecular ratio favors the water cluster and extended chainlike network that percolates within the PG matrix, where the propofol is in the dispersed state. We estimated that the radial distribution function indicates higher ordered water microstructure around propofol for high PG content, as compared to the lower PG content in the PG/water mixture. So, the hydrophilic PG regulates the stability of the water micronetwork around propofol and its solubility in the binary mixture. We observed that the translational and rotational mobility of water belonging to the propofol microsolvation shell is hindered for high PG content and relaxed toward the low PG molecular ratio in the PG/water mixture. It has been noticed that the structural relaxation of the hydrogen bond formed between the propofol and the water molecules present in the propofol microsolvation shell for all five compositions is found to be slower for high PG content and becomes faster on the way to low PG content in the mixture. Simultaneously, we calculated the intermittent residence time correlation function of the water molecules belonging to the microsolvation shell around the propofol for five different compositions and found a faster short time decay followed up with long time components. Again, the origin of such long time decay is primarily from the structural relaxation of the microsolvation shell around the propofol, where the high PG content shows the slower structural relaxation that turns faster as the PG content approaches to the other end of the compositions. So, our studies showed that the slower structural relaxation of the microsolvation shell around propofol for a high PG molecular ratio in the PG/water mixture correlate well with the extensive ordering of the water microstructure and restricted water mobility and facilitates the dissolution process of propofol in the binary mixture.

Alchemical Free Energy Calculations on Membrane-Associated Proteins

Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.

Effects of general anesthetics on synaptic transmission and plasticity

General anesthetics depress excitatory and/or enhance inhibitory synaptic transmission principally by modulating the function of glutamatergic or GABAergic synapses, respectively, with relative anesthetic agent-specific mechanisms. Synaptic signaling proteins, including ligand- and voltage-gated ion channels, are targeted by general anesthetics to modulate various synaptic mechanisms, including presynaptic neurotransmitter release, postsynaptic receptor signaling, and dendritic spine dynamics to produce their characteristic acute neurophysiological effects. As synaptic structure and plasticity mediate higher-order functions such as learning and memory, long-term synaptic dysfunction following anesthesia may lead to undesirable neurocognitive consequences depending on the specific anesthetic agent and the vulnerability of the population. Here we review the cellular and molecular mechanisms of transient and persistent general anesthetic alterations of synaptic transmission and plasticity.

Open-channel blocking action of volatile anaesthetics desflurane and sevoflurane on human voltage-gated Kv 1.5 channel

BACKGROUND AND PURPOSE

Volatile anaesthetics have been shown to differentially modulate mammalian Shaker-related voltage-gated potassium (K_v 1.x) channels. This study was designed to investigate molecular and cellular mechanisms underlying the modulatory effects of desflurane or sevoflurane on human K_v 1.5 (hK_v 1.5) channels.

EXPERIMENTAL APPROACH

Thirteen single-point mutations were constructed within pore domain of hK_v 1.5 channel using site-directed mutagenesis. The effects of desflurane or sevoflurane on heterologously expressed wild-type and mutant hK_v 1.5 channels were examined by whole-cell patch-clamp technique. A computer simulation was conducted to predict the docking pose of desflurane or sevoflurane within hK_v 1.5 channel.

KEY RESULTS

Both desflurane and sevoflurane increased hK_v 1.5 current at mild depolarizations but decreased it at strong depolarizations, indicating that these anaesthetics produce both stimulatory and inhibitory actions on hK_v 1.5 channels. The inhibitory effect of desflurane or sevoflurane on hK_v 1.5 channels arose primarily from its open-channel blocking action. The inhibitory action of desflurane or sevoflurane on hK_v 1.5 channels was significantly attenuated in T480A, V505A, and I508A mutant channels, compared with wild-type channel. Computational docking simulation predicted that desflurane or sevoflurane resides within the inner cavity of channel pore and has contact with Thr479, Thr480, Val505, and Ile508.

CONCLUSION AND IMPLICATIONS

Desflurane and sevoflurane exert an open-channel blocking action on hK_v 1.5 channels by functionally interacting with specific amino acids located within the channel pore. This study thus identifies a novel molecular basis mediating inhibitory modulation of hK_v 1.5 channels by desflurane and sevoflurane.

Identification of the Initial Steps in Signal Transduction in the α4β2 Nicotinic Receptor: Insights from Equilibrium and Nonequilibrium Simulations

Nicotinic acetylcholine receptors (nAChRs) modulate synaptic transmission in the nervous system. These receptors have emerged as therapeutic targets in drug discovery for treating several conditions, including Alzheimer's disease, pain, and nicotine addiction. In this in silico study, we use a combination of equilibrium and nonequilibrium molecular dynamics simulations to map dynamic and structural changes induced by nicotine in the human α4β2 nAChR. They reveal a striking pattern of communication between the extracellular binding pockets and the transmembrane domains (TMDs) and show the sequence of conformational changes associated with the initial steps in this process. We propose a general mechanism for signal transduction for Cys-loop receptors: the mechanistic steps for communication proceed firstly through loop C in the principal subunit, and are subsequently transmitted, gradually and cumulatively, to loop F of the complementary subunit, and then to the TMDs through the M2-M3 linker.

Anesthetics: from modes of action to unconsciousness and neurotoxicity

Modern anesthetic compounds and advanced monitoring tools have revolutionized the field of medicine, allowing for complex surgical procedures to occur safely and effectively. Faster induction times and quicker recovery periods of current anesthetic agents have also helped reduce health care costs significantly. Moreover, extensive research has allowed for a better understanding of anesthetic modes of action, thus facilitating the development of more effective and safer compounds. Notwithstanding the realization that anesthetics are a prerequisite to all surgical procedures, evidence is emerging to support the notion that exposure of the developing brain to certain anesthetics may impact future brain development and function. Whereas the data in support of this postulate from human studies is equivocal, the vast majority of animal research strongly suggests that anesthetics are indeed cytotoxic at multiple brain structure and function levels. In this review, we first highlight various modes of anesthetic action and then debate the evidence of harm from both basic science and clinical studies perspectives. We present evidence from animal and human studies vis-à-vis the possible detrimental effects of anesthetic agents on both the young developing and the elderly aging brain while discussing potential ways to mitigate these effects. We hope that this review will, on the one hand, invoke debate vis-à-vis the evidence of anesthetic harm in young children and the elderly, and on the other hand, incentivize the search for better and less toxic anesthetic compounds.

Propofol inhibits the voltage-gated sodium channel NaChBac at multiple sites

Voltage-gated sodium (NaV) channels are important targets of general anesthetics, including the intravenous anesthetic propofol. Electrophysiology studies on the prokaryotic NaV channel NaChBac have demonstrated that propofol promotes channel activation and accelerates activation-coupled inactivation, but the molecular mechanisms of these effects are unclear. Here, guided by computational docking and molecular dynamics simulations, we predict several propofol-binding sites in NaChBac. We then strategically place small fluorinated probes at these putative binding sites and experimentally quantify the interaction strengths with a fluorinated propofol analogue, 4-fluoropropofol. In vitro and in vivo measurements show that 4-fluoropropofol and propofol have similar effects on NaChBac function and nearly identical anesthetizing effects on tadpole mobility. Using quantitative analysis by 19F-NMR saturation transfer difference spectroscopy, we reveal strong intermolecular cross-relaxation rate constants between 4-fluoropropofol and four different regions of NaChBac, including the activation gate and selectivity filter in the pore, the voltage sensing domain, and the S4-S5 linker. Unlike volatile anesthetics, 4-fluoropropofol does not bind to the extracellular interface of the pore domain. Collectively, our results show that propofol inhibits NaChBac at multiple sites, likely with distinct modes of action. This study provides a molecular basis for understanding the net inhibitory action of propofol on NaV channels.

Ketamine Inhibition of the Pentameric Ligand-Gated Ion Channel GLIC

Ketamine inhibits pentameric ligand-gated ion channels (pLGICs), including the bacterial pLGIC from Gloeobacter violaceus (GLIC). The crystal structure of GLIC shows R-ketamine bound to an extracellular intersubunit cavity. Here, we performed molecular dynamics simulations of GLIC in the absence and presence of R- or S-ketamine. No stable binding of S-ketamine in the original cavity was observed in the simulations, largely due to its unfavorable access to residue D154, which provides important electrostatic interactions to stabilize R-ketamine binding. Contrary to the symmetric binding shown in the crystal structure, R-ketamine moved away from some of the binding sites and was bound to GLIC asymmetrically at the end of simulations. The asymmetric binding is consistent with the experimentally measured negative cooperativity of ketamine binding to GLIC. In the presence of R-ketamine, all subunits showed changes in structure and dynamics, irrespective of binding stability; the extracellular intersubunit cavity expanded and intersubunit electrostatic interactions involved in channel activation were altered. R-ketamine binding promoted a conformational shift toward closed GLIC. Conformational changes near the ketamine-binding site were propagated to the interface between the extracellular and transmembrane domains, and further to the pore-lining TM2 through two pathways: pre-TM1 and the β1-β2 loop. Both signaling pathways have been predicted previously using the perturbation-based Markovian transmission model. The study provides a structural and dynamics basis for the inhibitory modulation of ketamine on pLGICs.

IPPA08 allosterically enhances the action of imidacloprid on nicotinic acetylcholine receptors

Our previous study showed that IPPA08, a cis-configuration neonicotinoid compound with unique oxabridged substructure, acted as a specific synergist to neonicotinoid insecticides targeting nicotinic acetylcholine receptors (nAChRs). Heteropentamer nAChRs have diverse characteristics and can form canonical and noncanonical subunit interfaces. While canonical interfaces have been exploited as targets of many drugs, noncanonical interfaces have received less attention. In this study, the mechanism of IPPA08 synergism was evaluated on hybrid nAChRs consisting of three α1 subunits from the brown planthopper and two rat β1 subunits (Nlα1/rβ2) expressed in Xenopus oocytes. IPPA08 alone evoked inward currents, but only at very high concentrations, greater than 1 mM. However, at concentrations below 200 μM, IPPA08 slowed the decay of inward currents evoked by imidacloprid, but not by acetylcholine, and also increased the sensitivity of Nlα1/rβ2 to imidacloprid. Both modulations by IPPA08 were concentration-dependent in the same concentration range of 10-150 μM. Experimentally induced mutations in canonical (α+/β-) and noncanonical (β+/α-) interfaces of Nlα1/rβ2 receptors were also examined to evaluate the presence of possible binding sites for IPPA08 on the receptors. Our results showed that mutations in the canonical interfaces affected only the potency of IPPA08 as an agonist, while mutations in the noncanonical interfaces affected only the synergistic action of IPPA08. Based on these results, we propose that at low concentrations IPPA08 can act as a positive allosteric modulator of noncanonical interfaces, and likely slow the decay of currents through stabilizing the open-channel state caused by the action of imidacloprid on canonical interfaces.

Heteromeric α7β2 Nicotinic Acetylcholine Receptors in the Brain

The α7 nicotinic acetylcholine receptor (α7 nAChR) is highly expressed in the brain, where it maintains various neuronal functions including (but not limited to) learning and memory. In addition, the protein expression levels of α7 nAChRs are altered in various brain disorders. The classic rule governing α7 nAChR assembly in the mammalian brain was that it was assembled from five α7 subunits to form a homomeric receptor pentamer. However, emerging evidence demonstrates the presence of heteromeric α7 nAChRs in heterologously expressed systems and naturally in brain neurons, where α7 subunits are co-assembled with β2 subunits to form a novel type of α7β2 nAChR. Interestingly, the α7β2 nAChR exhibits distinctive function and pharmacology from traditional homomeric α7 nAChRs. We review recent advances in probing the distribution, function, pharmacology, pathophysiology, and stoichiometry of the heteromeric α7β2 nAChR, which have provided new insights into the understanding of a novel target of cholinergic signaling.

Direct Pore Binding as a Mechanism for Isoflurane Inhibition of the Pentameric Ligand-gated Ion Channel ELIC

Pentameric ligand-gated ion channels (pLGICs) are targets of general anesthetics, but molecular mechanisms underlying anesthetic action remain debatable. We found that ELIC, a pLGIC from Erwinia chrysanthemi, can be functionally inhibited by isoflurane and other anesthetics. Structures of ELIC co-crystallized with isoflurane in the absence or presence of an agonist revealed double isoflurane occupancies inside the pore near T237(6′) and A244(13′). A pore-radius contraction near the extracellular entrance was observed upon isoflurane binding. Electrophysiology measurements with a single-point mutation at position 6′ or 13′ support the notion that binding at these sites renders isoflurane inhibition. Molecular dynamics simulations suggested that isoflurane binding was more stable in the resting than in a desensitized pore conformation. This study presents compelling evidence for a direct pore-binding mechanism of isoflurane inhibition, which has a general implication for inhibitory action of general anesthetics on pLGICs.

Pharmacological and molecular studies on the interaction of varenicline with different nicotinic acetylcholine receptor subtypes. Potential mechanism underlying partial agonism at human α4β2 and α3β4 subtypes

To determine the structural components underlying differences in affinity, potency, and selectivity of varenicline for several human (h) nicotinic acetylcholine receptors (nAChRs), functional and structural experiments were performed. The Ca2+ influx results established that: (a) varenicline activates (μM range) nAChR subtypes with the following rank sequence: hα7>hα4β4>hα4β2>hα3β4>>>hα1β1γδ; (b) varenicline binds to nAChR subtypes with the following affinity order (nM range): hα4β2~hα4β4>hα3β4>hα7>>>Torpedo α1β1γδ. The molecular docking results indicating that more hydrogen bond interactions are apparent for α4-containing nAChRs in comparison to other nAChRs may explain the observed higher affinity; and that (c) varenicline is a full agonist at hα7 (101%) and hα4β4 (93%), and a partial agonist at hα4β2 (20%) and hα3β4 (45%), relative to (±)-epibatidine. The allosteric sites found at the extracellular domain (EXD) of hα3β4 and hα4β2 nAChRs could explain the partial agonistic activity of varenicline on these nAChR subtypes. Molecular dynamics simulations show that the interaction of varenicline to each allosteric site decreases the capping of Loop C at the hα4β2 nAChR, suggesting that these allosteric interactions limit the initial step in the gating process. In conclusion, we propose that in addition to hα4β2 nAChRs, hα4β4 nAChRs can be considered as potential targets for the clinical activity of varenicline, and that the allosteric interactions at the hα3β4- and hα4β2-EXDs are alternative mechanisms underlying partial agonism at these nAChRs.

Atomistic models of general anesthetics for use in in silico biological studies

While small molecules have been used to induce anesthesia in a clinical setting for well over a century, a detailed understanding of the molecular mechanism remains elusive. In this study, we utilize ab initio calculations to develop a novel set of CHARMM-compatible parameters for the ubiquitous modern anesthetics desflurane, isoflurane, sevoflurane, and propofol for use in molecular dynamics (MD) simulations. The parameters generated were rigorously tested against known experimental physicochemical properties including dipole moment, density, enthalpy of vaporization, and free energy of solvation. In all cases, the anesthetic parameters were able to reproduce experimental measurements, signifying the robustness and accuracy of the atomistic models developed. The models were then used to study the interaction of anesthetics with the membrane. Calculation of the potential of mean force for inserting the molecules into a POPC bilayer revealed a distinct energetic minimum of 4–5 kcal/mol relative to aqueous solution at the level of the glycerol backbone in the membrane. The location of this minimum within the membrane suggests that anesthetics partition to the membrane prior to binding their ion channel targets, giving context to the Meyer–Overton correlation. Moreover, MD simulations of these drugs in the membrane give rise to computed membrane structural parameters, including atomic distribution, deuterium order parameters, dipole potential, and lateral stress profile, that indicate partitioning of anesthetics into the membrane at the concentration range studied here, which does not appear to perturb the structural integrity of the lipid bilayer. These results signify that an indirect, membrane-mediated mechanism of channel modulation is unlikely.

Paying attention to smell: cholinergic signaling in the olfactory bulb

The tractable, layered architecture of the olfactory bulb (OB), and its function as a relay between odor input and higher cortical processing, makes it an attractive model to study how sensory information is processed at a synaptic and circuit level. The OB is also the recipient of strong neuromodulatory inputs, chief among them being the central cholinergic system. Cholinergic axons from the basal forebrain modulate the activity of various cells and synapses within the OB, particularly the numerous dendrodendritic synapses, resulting in highly variable responses of OB neurons to odor input that is dependent upon the behavioral state of the animal. Behavioral, electrophysiological, anatomical, and computational studies examining the function of muscarinic and nicotinic cholinergic receptors expressed in the OB have provided valuable insights into the role of acetylcholine (ACh) in regulating its function. We here review various studies examining the modulation of OB function by cholinergic fibers and their target receptors, and provide putative models describing the role that cholinergic receptor activation might play in the encoding of odor information.

Purification and functional characterization of a truncated human α4β2 nicotinic acetylcholine receptor

Nicotinic acetylcholine receptors (nAChR) are abundant in the brain and are essential in cognitive function, learning and memory. Previous efforts on α4β2 nAChR had been focused on functional and pharmacological characterization, where high expression yield is not essential. For structural studies though, large amounts of pure protein is important but heterologous overexpression of membrane proteins can be a burdensome task, especially if high amounts are required. In the current study, a truncated mutant of the human α4β2 nAChR was designed in order to improve expression and solubility and to obtain material suitable for high resolution structural studies. We showed that the wild type α4β2 nAChR presented low expression and solubilization yield both of which were improved with the truncated construct. The truncated nAChR showed similar binding profile to the wild type, was purified by a two-step chromatography and isolated in high purity and adequate quantity.

Pentameric Ligand-gated Ion Channels : Insights from Computation

Pentameric ligand-gated ion channels (pLGICs) conduct upon the binding of an agonist and are fundamental to neurotransmission. New insights into the complex mechanisms underlying pLGIC gating, ion selectivity, and modulation have recently been gained via a series of crystal structures in prokaryotes and C .elegans, as well as computational studies relying on these structures. Here we review contributions from a variety of computational approaches, including normal mode analysis, automated docking, and fully atomistic molecular dynamics simulation. Examples from our own research, particularly concerning interactions with general anesthetics and lipids, are used to illustrate predictive results complementary to crystallographic studies.

Insights into distinct modulation of α7 and α7β2 nicotinic acetylcholine receptors by the volatile anesthetic isoflurane

Nicotinic acetylcholine receptors (nAChRs) are targets of general anesthetics, but functional sensitivity to anesthetic inhibition varies dramatically among different subtypes of nAChRs. Potential causes underlying different functional responses to anesthetics remain elusive. Here we show that in contrast to the α7 nAChR, the α7β2 nAChR is highly susceptible to inhibition by the volatile anesthetic isoflurane in electrophysiology measurements. Isoflurane-binding sites in β2 and α7 were found at the extracellular and intracellular end of their respective transmembrane domains using NMR. Functional relevance of the identified β2 site was validated via point mutations and subsequent functional measurements. Consistent with their functional responses to isoflurane, β2 but not α7 showed pronounced dynamics changes, particularly for the channel gate residue Leu-249(9'). These results suggest that anesthetic binding alone is not sufficient to generate functional impact; only those sites that can modulate channel dynamics upon anesthetic binding will produce functional effects.

GABA binding to an insect GABA receptor: a molecular dynamics and mutagenesis study

RDL receptors are GABA-activated inhibitory Cys-loop receptors found throughout the insect CNS. They are a key target for insecticides. Here, we characterize the GABA binding site in RDL receptors using computational and electrophysiological techniques. A homology model of the extracellular domain of RDL was generated and GABA docked into the binding site. Molecular dynamics simulations predicted critical GABA binding interactions with aromatic residues F206, Y254, and Y109 and hydrophilic residues E204, S176, R111, R166, S176, and T251. These residues were mutated, expressed in Xenopus oocytes, and their functions assessed using electrophysiology. The data support the binding mechanism provided by the simulations, which predict that GABA forms many interactions with binding site residues, the most significant of which are cation-π interactions with F206 and Y254, H-bonds with E204, S205, R111, S176, T251, and ionic interactions with R111 and E204. These findings clarify the roles of a range of residues in binding GABA in the RDL receptor, and also show that molecular dynamics simulations are a useful tool to identify specific interactions in Cys-loop receptors.

NMR resolved multiple anesthetic binding sites in the TM domains of the α4β2 nAChR

The α4β2 nicotinic acetylcholine receptor (nAChR) has significant roles in nervous system function and disease. It is also a molecular target of general anesthetics. Anesthetics inhibit the α4β2 nAChR at clinically relevant concentrations, but their binding sites in α4β2 remain unclear. The recently determined NMR structures of the α4β2 nAChR transmembrane (TM) domains provide valuable frameworks for identifying the binding sites. In this study, we performed solution NMR experiments on the α4β2 TM domains in the absence and presence of halothane and ketamine. Both anesthetics were found in an intra-subunit cavity near the extracellular end of the β2 transmembrane helices, homologous to a common anesthetic binding site observed in X-ray structures of anesthetic-bound GLIC (Nury et al., [32]). Halothane, but not ketamine, was also found in cavities adjacent to the common anesthetic site at the interface of α4 and β2. In addition, both anesthetics bound to cavities near the ion selectivity filter at the intracellular end of the TM domains. Anesthetic binding induced profound changes in protein conformational exchanges. A number of residues, close to or remote from the binding sites, showed resonance signal splitting from single to double peaks, signifying that anesthetics decreased conformation exchange rates. It was also evident that anesthetics shifted population of two conformations. Altogether, the study comprehensively resolved anesthetic binding sites in the α4β2 nAChR. Furthermore, the study provided compelling experimental evidence of anesthetic-induced changes in protein dynamics, especially near regions of the hydrophobic gate and ion selectivity filter that directly regulate channel functions.

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.

NMR structures of the transmembrane domains of the α4β2 nAChR

The α4β2 nicotinic acetylcholine receptor (nAChR) is the predominant heteromeric subtype of nAChRs in the brain, which has been implicated in numerous neurological conditions. The structural information specifically for the α4β2 and other neuronal nAChRs is presently limited. In this study, we determined structures of the transmembrane (TM) domains of the α4 and β2 subunits in lauryldimethylamine-oxide (LDAO) micelles using solution NMR spectroscopy. NMR experiments and size exclusion chromatography-multi-angle light scattering (SEC-MALS) analysis demonstrated that the TM domains of α4 and β2 interacted with each other and spontaneously formed pentameric assemblies in the LDAO micelles. The Na(+) flux assay revealed that α4β2 formed Na(+) permeable channels in lipid vesicles. Efflux of Na(+) through the α4β2 channels reduced intra-vesicle Sodium Green™ fluorescence in a time-dependent manner that was not observed in vesicles without incorporating α4β2. The study provides structural insight into the TM domains of the α4β2 nAChR. It offers a valuable structural framework for rationalizing extensive biochemical data collected previously on the α4β2 nAChR and for designing new therapeutic modulators.

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.

NMR structure and dynamics of a designed water-soluble transmembrane domain of nicotinic acetylcholine receptor

The nicotinic acetylcholine receptor (nAChR) is an important therapeutic target for a wide range of pathophysiological conditions, for which rational drug designs often require receptor structures at atomic resolution. Recent proof-of-concept studies demonstrated a water-solubilization approach to structure determination of membrane proteins by NMR (Slovic et al., PNAS, 101: 1828-1833, 2004; Ma et al., PNAS, 105: 16537-42, 2008). We report here the computational design and experimental characterization of WSA, a water-soluble protein with ~83% sequence identity to the transmembrane (TM) domain of the nAChR α1 subunit. Although the design was based on a low-resolution structural template, the resulting high-resolution NMR structure agrees remarkably well with the recent crystal structure of the TM domains of the bacterial Gloeobacter violaceus pentameric ligand-gated ion channel (GLIC), demonstrating the robustness and general applicability of the approach. NMR T(2) dispersion measurements showed that the TM2 domain of the designed protein was dynamic, undergoing conformational exchange on the NMR timescale. Photoaffinity labeling with isoflurane and propofol photolabels identified a common binding site in the immediate proximity of the anesthetic binding site found in the crystal structure of the anesthetic-GLIC complex. Our results illustrate the usefulness of high-resolution NMR analyses of water-solubilized channel proteins for the discovery of potential drug binding sites.

Alcohol-binding sites in distinct brain proteins: the quest for atomic level resolution

Defining the sites of action of ethanol on brain proteins is a major prerequisite to understanding the molecular pharmacology of this drug. The main barrier to reaching an atomic-level understanding of alcohol action is the low potency of alcohols, ethanol in particular, which is a reflection of transient, low-affinity interactions with their targets. These mechanisms are difficult or impossible to study with traditional techniques such as radioligand binding or spectroscopy. However, there has been considerable recent progress in combining X-ray crystallography, structural modeling, and site-directed mutagenesis to define the sites and mechanisms of action of ethanol and related alcohols on key brain proteins. We review such insights for several diverse classes of proteins including inwardly rectifying potassium, transient receptor potential, and neurotransmitter-gated ion channels, as well as protein kinase C epsilon. Some common themes are beginning to emerge from these proteins, including hydrogen bonding of the hydroxyl group and van der Waals interactions of the methylene groups of ethanol with specific amino acid residues. The resulting binding energy is proposed to facilitate or stabilize low-energy state transitions in the bound proteins, allowing ethanol to act as a "molecular lubricant" for protein function. We discuss evidence for characteristic, discrete alcohol-binding sites on protein targets, as well as evidence that binding to some proteins is better characterized by an interaction region that can accommodate multiple molecules of ethanol.

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.

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.

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.

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.

Halothane solvation in water and organic solvents from molecular simulations with new polarizable potential function

The partitioning of a substrate from one phase into another is a complex process with widespread applications: from chemical technology to the pharmaceutical industry. One particularly well-known and well-studied example is 2-bromo-2-chloro-1,1,1-trifluoroethane (halothane) trafficking through the lipid bilayer. Halothane is a model volatile anesthetic known to impact functions of model lipid bilayers, altering the structure and thickness upon its partitioning from the bulk phase. A number of theoretical and experimental investigations suggest the importance of electronic polarizability, determining a preference for halothane to partition in the interfacial systems as in lipid bilayers or binary solvents. The recently published protocol for the development of polarizable force fields based on the classical Drude model has provided fresh impetus to efforts directed at understanding the molecular principles governing complex thermodynamics of the hydrophobic hydration. Here, molecular simulations were combined with free energy simulations to study solvation of halothane in polarizable water and methanol. The absolute free energy of halothane solvation in different solvents (water, methanol, and n-hexane) has been evaluated for additive and polarizable models. It was found that both additive and polarizable models provide an adequate description of the halothane solvation in high-dielectric (polar) solvents such as water, but explicit accounting for electronic polarization is imperative for a correct description of the solvation thermodynamics in nonpolar systems. To study halothane dynamics in binary mixtures, all-atom molecular dynamics (MD) simulations for halothane-methanol mixtures in a wide range of concentrations were performed alongside an analysis of structural organization, dynamics, and thermodynamic properties to dissect the molecular determinants of the halothane solvation in polar and amphiphilic liquids such as methanol. Additionally, a theoretical test of the hypothesis on the weak hydrogen bonding of halothane and methanol in the condensed phase is provided, which was presented on the basis of spectroscopic analysis of the C-H vibrations in different gas-phase complexes. The simulations performed in the condensed phase suggest that hydrophobic interactions between halothane and methanol play a dominant role in preferential solvation.

NMR structure of the transmembrane domain of the n-acetylcholine receptor beta2 subunit

Nicotinic acetylcholine receptors (nAChRs) are involved in fast synaptic transmission in the central and peripheral nervous system. Among the many different types of subunits in nAChRs, the beta2 subunit often combines with the alpha4 subunit to form alpha4beta2 pentameric channels, the most abundant subtype of nAChRs in the brain. Besides computational predictions, there is limited experimental data available on the structure of the beta2 subunit. Using high-resolution NMR spectroscopy, we solved the structure of the entire transmembrane domain (TM1234) of the beta2 subunit. We found that TM1234 formed a four-helix bundle in the absence of the extracellular and intracellular domains. The structure exhibited many similarities to those previously determined for the Torpedo nAChR and the bacterial ion channel GLIC. We also assessed the influence of the fourth transmembrane helix (TM4) on the rest of the domain. Although secondary structures and tertiary arrangements were similar, the addition of TM4 caused dramatic changes in TM3 dynamics and subtle changes in TM1 and TM2. Taken together, this study suggests that the structures of the transmembrane domains of these proteins are largely shaped by determinants inherent in their sequence, but their dynamics may be sensitive to modulation by tertiary and quaternary contacts.

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.