NMR study of general anesthetic interaction with nAChR beta2 subunit

Molecular Aspects Implicated in Dantrolene Selectivity with Respect to Ryanodine Receptor Isoforms

Dantrolene is an intra-cellularly acting skeletal muscle relaxant used for the treatment of the rare genetic disorder, malignant hyperthermia (MH). In most cases, MH susceptibility is caused by dysfunction of the skeletal ryanodine receptor (RyR1) harboring one of nearly 230 single-point MH mutations. The therapeutic effect of dantrolene is the result of a direct inhibitory action on the RyR1 channel, thus suppressing aberrant Ca2+ release from the sarcoplasmic reticulum. Despite the almost identical dantrolene-binding sequence exits in all three mammalian RyR isoforms, dantrolene appears to be an isoform-selective inhibitor. Whereas RyR1 and RyR3 channels are competent to bind dantrolene, the RyR2 channel, predominantly expressed in the heart, is unresponsive. However, a large body of evidence suggests that the RyR2 channel becomes sensitive to dantrolene-mediated inhibition under certain pathological conditions. Although a consistent picture of the dantrolene effect emerges from in vivo studies, in vitro results are often contradictory. Hence, our goal in this perspective is to provide the best possible clues to the molecular mechanism of dantrolene’s action on RyR isoforms by identifying and discussing potential sources of conflicting results, mainly coming from cell-free experiments. Moreover, we propose that, specifically in the case of the RyR2 channel, its phosphorylation could be implicated in acquiring the channel responsiveness to dantrolene inhibition, interpreting functional findings in the structural context.

X-Ray Crystallographic Studies for Revealing Binding Sites of General Anesthetics in Pentameric Ligand-Gated Ion Channels

X-ray crystallography is a powerful tool in structural biology and can offer insight into structured-based understanding of general anesthetic action on various relevant molecular targets, including pentameric ligand-gated ion channels (pLGICs). In this chapter, we outline the procedures for expression and purification of pLGICs. Optimization of crystallization conditions, especially to achieve high-resolution structures of pLGICs bound with general anesthetics, is also presented. Case studies of pLGICs bound with the volatile general anesthetic isoflurane, 2-bromoethanol, and the intravenous general anesthetic ketamine are revisited.

Solution NMR Studies of Anesthetic Interactions with Ion Channels

NMR spectroscopy is one of the major tools to provide atomic resolution protein structural information. It has been used to elucidate the molecular details of interactions between anesthetics and ion channels, to identify anesthetic binding sites, and to characterize channel dynamics and changes introduced by anesthetics. In this chapter, we present solution NMR methods essential for investigating interactions between ion channels and general anesthetics, including both volatile and intravenous anesthetics. Case studies are provided with a focus on pentameric ligand-gated ion channels and the voltage-gated sodium channel NaChBac.

Assessing anesthetic activity through modulation of the membrane dipole potential

There is great individual variation in response to general anesthetics (GAs) leading to difficulties in optimal dosing and sometimes even accidental awareness during general anesthesia (AAGA). AAGA is a rare, but potentially devastating, complication affecting between 0.1% and 2% of patients undergoing surgery. The development of novel personalized screening techniques to accurately predict a patient’s response to GAs and the risk of AAGA remains an unmet clinical need. In the present study, we demonstrate the principle of using a fluorescent reporter of the membrane dipole potential, di-8-ANEPPs, as a novel method to monitor anesthetic activity using a well-described inducer/noninducer pair. The membrane dipole potential has previously been suggested to contribute a novel mechanism of anesthetic action. We show that the fluorescence ratio of di-8-ANEPPs changed in response to physiological concentrations of the anesthetic, 1-chloro-1,2,2-trifluorocyclobutane (F₃), but not the structurally similar noninducer, 1,2-dichlorohexafluorocyclobutane (F₆), to artificial membranes and in vitro retinal cell systems. Modulation of the membrane dipole provides an explanation to overcome the limitations associated with the alternative membrane-mediated mechanisms of GA action. Furthermore, by combining this technique with noninvasive retinal imaging technologies, we propose that this technique could provide a novel and noninvasive technique to monitor GA susceptibility and identify patients at risk of AAGA.

Fluorine-19 NMR and computational quantification of isoflurane binding to the voltage-gated sodium channel NaChBac

Voltage-gated sodium channels (Na_V) play an important role in general anesthesia. Electrophysiology measurements suggest that volatile anesthetics such as isoflurane inhibit Na_V by stabilizing the inactivated state or altering the inactivation kinetics. Recent computational studies suggested the existence of multiple isoflurane binding sites in Na_V, but experimental binding data are lacking. Here we use site-directed placement of ¹⁹F probes in NMR experiments to quantify isoflurane binding to the bacterial voltage-gated sodium channel NaChBac. ¹⁹F probes were introduced individually to S129 and L150 near the S4-S5 linker, L179 and S208 at the extracellular surface, T189 in the ion selectivity filter, and all phenylalanine residues. Quantitative analyses of ¹⁹F NMR saturation transfer difference (STD) spectroscopy showed a strong interaction of isoflurane with S129, T189, and S208; relatively weakly with L150; and almost undetectable with L179 and phenylalanine residues. An orientation preference was observed for isoflurane bound to T189 and S208, but not to S129 and L150. We conclude that isoflurane inhibits NaChBac by two distinct mechanisms: (i) as a channel blocker at the base of the selectivity filter, and (ii) as a modulator to restrict the pivot motion at the S4-S5 linker and at a critical hinge that controls the gating and inactivation motion of S6.

Molecular mechanisms of drug action: an emerging view

Volatile anesthetics serve as useful probes of a conserved biological process that is essential to the proper functioning of the central nervous system. A kinetic and thermodynamic analysis of their unusual pharmacological and physiological characteristics has led to a general, predictive theory in which small molecules that adsorb to membranes modulate ion channel function by altering physical properties of membrane bilayers. A kinetic model that is both parsimonious and falsifiable has been developed to test this mechanism. This theory leads to predictions about the structure, function, origin, and evolution of synapses, the etiology of several diseases and disease symptoms affecting the brain, and the mechanism of action of several drugs that are used therapeutically. Neuronal membranes may offer an appealing drug target, given the large number of compounds that adsorb to interfaces and hence membranes.

Atomistic study of lipid membranes containing chloroform: looking for a lipid-mediated mechanism of anesthesia

The molecular mechanism of general anesthesia is still a controversial issue. Direct effect by linking of anesthetics to proteins and indirect action on the lipid membrane properties are the two hypotheses in conflict. Atomistic simulations of different lipid membranes subjected to the effect of small volatile organohalogen compounds are used to explore plausible lipid-mediated mechanisms. Simulations of homogeneous membranes reveal that electrostatic potential and lateral pressure transversal profiles are affected differently by chloroform (anesthetic) and carbon tetrachloride (non-anesthetic). Simulations of structured membranes that combine ordered and disordered regions show that chloroform molecules accumulate preferentially in highly disordered lipid domains, suggesting that the combination of both lateral and transversal partitioning of chloroform in the cell membrane could be responsible of its anesthetic action.

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.

The Role of Atomic Polarization in the Thermodynamics of Chloroform Partitioning to Lipid Bilayers

In spite of extensive research and use in medical practice, the precise molecular mechanism of volatile anesthetic action remains unknown. The distribution of anesthetics within lipid bilayers and potential targeting to membrane proteins is thought to be central to therapeutic function. Therefore, obtaining a molecular level understanding of volatile anesthetic partitioning into lipid bilayers is of vital importance to modern pharmacology. In this study we investigate the partitioning of the prototypical anesthetic, chloroform, into lipid bilayers and different organic solvents using molecular dynamics simulations with potential models ranging from simplified coarse-grained MARTINI to additive and polarizable CHARMM all-atom force fields. Many volatile anesthetics display significant inducible dipole moments, which correlate with their potency, yet the exact role of molecular polarizability in their stabilization within lipid bilayers remains unknown. We observe that explicit treatment of atomic polarizability makes it possible to accurately reproduce solvation free energies in solvents with different polarities, allowing for quantitative studies in heterogeneous molecular distributions, such as lipid bilayers. We calculate the free energy profiles for chloroform crossing lipid bilayers to reveal a role of polarizability in modulating chloroform partitioning thermodynamics via the chloroform-induced dipole moment and highlight competitive binding to the membrane core and toward the glycerol backbone that may have significant implications for understanding anesthetic action.

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.

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.

Computational studies on the interactions of inhalational anesthetics with proteins

Despite the widespread clinical use of anesthetics since the 19th century, a clear understanding of the mechanism of anesthetic action has yet to emerge. On the basis of early experiments by Meyer, Overton, and subsequent researchers, the cell's lipid membrane was generally concluded to be the primary site of action of anesthetics. However, later experiments with lipid-free globular proteins, such as luciferase and apoferritin, shifted the focus of anesthetic action to proteins. Recent experimental studies, such as photoaffinity labeling and mutagenesis on membrane proteins, have suggested specific binding sites for anesthetic molecules, further strengthening the proteocentric view of anesthetic mechanism. With the increased availability of high-resolution crystal structures of ion channels and other integral membrane proteins, as well as the availability of powerful computers, the structure-function relationship of anesthetic-protein interactions can now be investigated in atomic detail. In this Account, we review recent experiments and related computer simulation studies involving interactions of inhalational anesthetics and proteins, with a particular focus on membrane proteins. Globular proteins have long been used as models for understanding the role of protein-anesthetic interactions and are accordingly examined in this Account. Using selected examples of membrane proteins, such as nicotinic acetyl choline receptor (nAChR) and potassium channels, we address the issues of anesthetic binding pockets in proteins, the role of conformation in anesthetic effects, and the modulation of local as well as global dynamics of proteins by inhaled anesthetics. In the case of nicotinic receptors, inhalational anesthetic halothane binds to the hydrophobic cavity close to the M2-M3 loop. This binding modulates the dynamics of the M2-M3 loop, which is implicated in allosterically transmitting the effects to the channel gate, thus altering the function of the protein. In potassium channels, anesthetic molecules preferentially potentiate the open conformation by quenching the motion of the aromatic residues implicated in the gating of the channel. These simulations suggest that low-affinity drugs (such as inhalational anesthetics) modulate the protein function by influencing local as well as global dynamics of proteins. Because of intrinsic experimental limitations, computational approaches represent an important avenue for exploring the mode of action of anesthetics. Molecular dynamics simulations-a computational technique frequently used in the general study of proteins-offer particular insight in the study of the interaction of inhalational anesthetics with membrane proteins.

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

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

Anesthetic effects on the structure and dynamics of the second transmembrane domains of nAChR alpha4beta2

Channel functions of the neuronal alpha4beta2 nicotinic acetylcholine receptor (nAChR), one of the most widely expressed subtypes in the brain, can be inhibited by volatile anesthetics. Our Na(+) flux experiments confirmed that the second transmembrane domains (TM2) of alpha4 and beta2 in 2:3 stoichiometry, (alpha4)(2)(beta2)(3), could form pentameric channels, whereas the alpha4 TM2 alone could not. The structure, topology, and dynamics of the alpha4 TM2 and (alpha4)(2)(beta2)(3) TM2 in magnetically aligned phospholipid bicelles were investigated using solid-state NMR spectroscopy in the absence and presence of halothane and isoflurane, two clinically used volatile anesthetics. (2)H NMR demonstrated that anesthetics increased lipid conformational heterogeneity. Such anesthetic effects on lipids became more profound in the presence of transmembrane proteins. PISEMA experiments on the selectively (15)N-labeled alpha4 TM2 showed that the TM2 formed transmembrane helices with tilt angles of 12 degrees +/-1 degrees and 16 degrees +/-1 degrees relative to the bicelle normal for the alpha4 and (alpha4)(2)(beta2)(3) samples, respectively. Anesthetics changed the tilt angle of the alpha4 TM2 from 12 degrees +/-1 degrees to 14 degrees +/-1 degrees , but had only a subtle effect on the tilt angle of the (alpha4)(2)(beta2)(3) TM2. A small degree of wobbling motion of the helix axis occurred in the (alpha4)(2)(beta2)(3) TM2. In addition, a subset of the (alpha4)(2)(beta2)(3) TM2 exhibited counterclockwise rotational motion around the helix axis on a time scale slower than 10(-4) s in the presence of anesthetics. Both helical tilting and rotational motions have been identified computationally as critical elements for ion channel functions. This study suggested that anesthetics could alter these motions to modulate channel functions.

Loosening and reorganization of fluid phospholipid bilayers by chloroform

The mixing behavior of an exchangeable phospholipid (A) with an exchangeable sterol (B) in host bilayers made from 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) containing varying concentrations of cholesterol has been examined via the nearest-neighbor recognition method. At low sterol concentrations (i.e., 2.5 mol %), the mixing between A and B is close to ideal. Incremental increases in the sterol concentration to 40 mol % led to net increases in the affinity between A and B. Similar mixing experiments that were carried out in the presence of chloroform showed a leveling effect, where moderate sterol-phospholipid affinity was observed in all cases. These results, together with the fact that the number of chloroform molecules that are absorbed per phospholipid is essentially constant and independent of the sterol content, support a model in which chloroform favors solvation of the phospholipids and a common membrane state is produced. Fluorescence measurements and Raman spectra have also shown that chloroform significantly loosens both cholesterol-poor and cholesterol-rich membranes made from DPPC. In a broader context, these results suggest a fundamentally new mechanism of anesthesia, where the anesthetic, by solvating the lipid components, profoundly changes the lateral organization of the lipid framework.

Anesthetic modulation of protein dynamics: insight from an NMR study

Mistic (membrane integrating sequence for translation of integral membrane protein constructs) comprises the four-alpha-helix bundle scaffold found in the transmembrane domains of the Cys-loop receptors that are plausible targets for general anesthetics. Nuclear magnetic resonance (NMR) studies of anesthetic halothane interaction with Mistic in dodecyl phosphocholine (DPC) micelles provide an experimental basis for understanding molecular mechanisms of general anesthesia. Halothane was found to interact directly with Mistic, mostly in the interfacial loop regions. Although the presence of halothane had little effect on Mistic structure, (15)N NMR relaxation dispersion measurements revealed that halothane affected Mistic's motion on the microsecond-millisecond time scale. Halothane shifted the equilibrium of chemical exchange in some residues and made the exchange faster or slower in comparison to the original state in the absence of halothane. The motion on the microsecond-millisecond time scale in several residues disappeared in response to the addition of halothane. Most of the residues experiencing halothane-induced dynamics changes also exhibited profound halothane-induced changes in chemical shift, suggesting that dynamics modification of these residues might result from their direct interaction with halothane molecules. Allosteric modulation by halothane also contributed to dynamics changes, as reflected in residues I52 and Y82 where halothane introduction brought about dynamics changes but not chemical shift changes. The study suggests that inhaled general anesthetics could act on proteins via altering protein motion on the microsecond-millisecond time scale, especially motion in the flexible loops that link different alpha helices. The validation of anesthetic effect on protein dynamics that are potentially correlated with protein functions is a critical step in unraveling the mechanisms of anesthetic action on proteins.