Unifying characteristics of sites of anesthetic action revealed by combined use of anesthetics and non-anesthetics

Critical role of water in the binding of volatile anesthetics to proteins

Numerous small molecules exhibit drug-like properties by low-affinity binding to proteins. Such binding is known to be influenced by water, the detailed picture of which, however, remains unclear. One particular example is the controversial role of water in the binding of general anesthetics to proteins as an essential step in general anesthesia. Here we demonstrate that a critical amount of hydration water is a prerequisite for anesthetic-protein binding. Using nuclear magnetic resonance, the concurrent adsorption of hydration water and bound anesthetics on model proteins are simultaneously measured. Halothane binding on proteins can only take place after protein hydration reaches a threshold hydration level of ∼0.31 g of water/g of proteins at the relative water vapor pressure of ∼0.95. Similar dependence on hydration is also observed for several other anesthetics. The ratio of anesthetic partial pressures at which two different anesthetics reach the same fractional load is correlated with the anesthetic potency. The binding of nonimmobilizers, which are structurally similar to known anesthetics but unable to produce anesthesia, does not occur even after the proteins are fully hydrated. Our results provide the first unambiguous experimental evidence that water is absolutely required to enable anesthetic-protein interactions, shedding new light on the general mechanism of molecular recognition and binding.

Four-alpha-helix bundle with designed anesthetic binding pockets. Part I: structural and dynamical analyses

The four-alpha-helix bundle mimics the transmembrane domain of the Cys-loop receptor family believed to be the protein target for general anesthetics. Using high resolution NMR, we solved the structure (Protein Data Bank ID: 2I7U) of a prototypical dimeric four-alpha-helix bundle, (Aalpha(2)-L1M/L38M)(2,) with designed specific binding pockets for volatile anesthetics. Two monomers of the helix-turn-helix motif form an antiparallel dimer as originally designed, but the high-resolution structure exhibits an asymmetric quaternary arrangement of the four helices. The two helices from the N-terminus to the linker (helices 1 and 1') are associated with each other in the dimer by the side-chain ring stacking of F12 and W15 along the long hydrophobic core and by a nearly perfect stretch of hydrophobic interactions between the complementary pairs of L4, L11, L18, and L25, all of which are located at the heptad e position along the helix-helix dimer interface. In comparison, the axes of the two helices from the linker to the C-terminus (helices 2 and 2') are wider apart from each other, creating a lateral access pathway around K47 from the aqueous phase to the center of the designed hydrophobic core. The site of the L38M mutation, which was previously shown to increase the halothane binding affinity by approximately 3.5-fold, is not part of the hydrophobic core presumably involved in the anesthetic binding but shows an elevated transverse relaxation (R(2)) rate. Qualitative analysis of the protein dynamics by reduced spectral density mapping revealed exchange contributions to the relaxation at many residues in the helices. This observation was confirmed by the quantitative analysis using the Modelfree approach and by the NMR relaxation dispersion measurements. The NMR structures and Autodock analysis suggest that the pocket with the most favorable amphipathic property for anesthetic binding is located between the W15 side chains at the center of the dimeric hydrophobic core, with the possibility of two additional minor binding sites between the F12 and F52 ring stacks of each monomer. The high-resolution structure of the designed anesthetic-binding protein offers unprecedented atomistic details about possible sites for anesthetic-protein interactions that are essential to the understanding of molecular mechanisms of general anesthesia.

Four-alpha-helix bundle with designed anesthetic binding pockets. Part II: halothane effects on structure and dynamics

As a model of the protein targets for volatile anesthetics, the dimeric four-alpha-helix bundle, (Aalpha(2)-L1M/L38M)(2), was designed to contain a long hydrophobic core, enclosed by four amphipathic alpha-helices, for specific anesthetic binding. The structural and dynamical analyses of (Aalpha(2)-L1M/L38M)(2) in the absence of anesthetics (another study) showed a highly dynamic antiparallel dimer with an asymmetric arrangement of the four helices and a lateral accessing pathway from the aqueous phase to the hydrophobic core. In this study, we determined the high-resolution NMR structure of (Aalpha(2)-L1M/L38M)(2) in the presence of halothane, a clinically used volatile anesthetic. The high-solution NMR structure, with a backbone root mean-square deviation of 1.72 A (2JST), and the NMR binding measurements revealed that the primary halothane binding site is located between two side-chains of W15 from each monomer, different from the initially designed anesthetic binding sites. Hydrophobic interactions with residues A44 and L18 also contribute to stabilizing the bound halothane. Whereas halothane produces minor changes in the monomer structure, the quaternary arrangement of the dimer is shifted by about half a helical turn and twists relative to each other, which leads to the closure of the lateral access pathway to the hydrophobic core. Quantitative dynamics analyses, including Modelfree analysis of the relaxation data and the Carr-Purcell-Meiboom-Gill transverse relaxation dispersion measurements, suggest that the most profound anesthetic effect is the suppression of the conformational exchange both near and remote from the binding site. Our results revealed a novel mechanism of an induced fit between anesthetic molecule and its protein target, with the direct consequence of protein dynamics changing on a global rather than a local scale. This mechanism may be universal to anesthetic action on neuronal proteins.

Large-scale molecular dynamics simulations of general anesthetic effects on the ion channel in the fully hydrated membrane: the implication of molecular mechanisms of general anesthesia

Interactions of volatile anesthetics with the central nervous system are characterized by low yet specific binding affinities. Although neurotransmitter-gated ion channels are considered the primary anesthetic targets, the mechanism of action at the molecular level remains elusive. We consider here the theoretical implications of channel dynamics on anesthetic action in a simplified membrane-channel system. Large-scale 2.2-ns all-atom molecular dynamics simulations were performed to study the effects of halothane, a clinical anesthetic, on a gramicidin A (gA) channel in a fully hydrated dimyristoyl phosphatidylcholine membrane. In agreement with experimental results, anesthetics preferentially target the anchoring residues at the channel-lipid-water interface. Although the anesthetic effect on channel structure is minimal, the presence of halothane profoundly affects channel dynamics. For 2.2-ns simulation, the rms fluctuation of gA backbone in the lipid core increases from approximately equal 1 A in the absence of anesthetics to approximately equal 1.5 A in the presence of halothane. Autocorrelation analysis reveals that halothane (i) has no effect on the subpicosecond librational motion, (ii) prolongs the backbone autocorrelation time in the 10- to 100-ps time scale, and (iii) significantly decreases the asymptotic values of generalized order parameter and correlation time of nanosecond motions for the inner but not the outer residues. The simulation results discount the viewpoint of a structure-function paradigm that overrates the importance of structural fitting between general anesthetics and yet-unidentified hydrophobic protein pockets. Instead, the results underscore the global, as opposed to local, effects of anesthetics on protein dynamics as the underlying mechanisms for the action of general anesthetics and possibly of other low-affinity drugs.

NMR structures of the second transmembrane domain of the human glycine receptor alpha(1) subunit: model of pore architecture and channel gating

Glycine receptors (GlyR) are the primary inhibitory receptors in the spinal cord and belong to a superfamily of ligand-gated ion channels (LGICs) that are extremely sensitive to low-affinity neurological agents such as general anesthetics and alcohols. The high-resolution pore architecture and the gating mechanism of this superfamily, however, remain unclear. The pore-lining second transmembrane (TM2) segments of the GlyR alpha(1) subunit are unique in that they form functional homopentameric channels with conductance characteristics nearly identical to those of an authentic receptor (Opella, S. J., J. Gesell, A. R. Valente, F. M. Marassi, M. Oblatt-Montal, W. Sun, A. F. Montiel, and M. Montal. 1997. Chemtracts Biochem. Mol. Biol. 10:153-174). Using NMR and circular dichroism (CD), we determined the high-resolution structures of the TM2 segment of human alpha(1) GlyR and an anesthetic-insensitive mutant (S267Y) in dodecyl phosphocholine (DPC) and sodium dodecyl sulfate (SDS) micelles. The NMR structures showed right-handed alpha-helices without kinks. A well-defined hydrophilic path, composed of side chains of G2', T6', T10', Q14', and S18', runs along the helical surfaces at an angle approximately 10-20 degrees relative to the long axis of the helices. The side-chain arrangement of the NMR-derived structures and the energy minimization of a homopentameric TM2 channel in a fully hydrated DMPC membrane using large-scale computation suggest a model of pore architecture in which simultaneous tilting movements of entire TM2 helices by a mere 10 degrees may be sufficient to account for the channel gating. The model also suggests that additional residues accessible from within the pore include L3', T7', T13', and G17'. A similar pore architecture and gating mechanism may apply to other channels in the same superfamily, including GABA(A), nACh, and 5-HT(3) receptors.

Membrane structural perturbations caused by anesthetics and nonimmobilizers: a molecular dynamics investigation

The structural perturbations of the fully hydrated dimyristoyl-phosphatidylcholine bilayer induced by the presence of hexafluoroethane C(2F6), a "nonimmobilizer," have been examined by molecular dynamics simulations and compared with the effects produced by halothane CF3CHBrCl, an "anesthetic," on a similar bilayer (DPPC) (Koubi et al., Biophys. J. 2000. 78:800). We find that the overall structure of the lipid bilayer and the zwitterionic head-group dipole orientation undergo only a slight modification compared with the pure lipid bilayer, with virtually no change in the potential across the interface. This is in contrast to the anesthetic case in which the presence of the molecule led to a large perturbation of the electrostatic potential across to the membrane interface. Similarly, the analysis of the structural and dynamical properties of the lipid core are unchanged in the presence of the nonimmobilizer although there is a substantial increase in the microscopic viscosity for the system containing the anesthetic. These contrasting perturbations of the lipid membrane caused by those quite similarly sized molecules may explain the difference in their physiological effects as anesthetics and nonimmobilizers, respectively.

General anesthetic binding to gramicidin A: the structural requirements

There is a distinct possibility that general anesthetics exert their action on the postsynaptic receptor channels. The structural requirements for anesthetic binding in transmembrane channels, however, are largely unknown. High-resolution (1)H nuclear magnetic resonance and direct photoaffinity labeling were used in this study to characterize the volatile anesthetic binding sites in gramicidin A (gA) incorporated into sodium dodecyl sulfate (SDS) micelles and into dimyristoylphosphatidylcholine (DMPC) bilayers, respectively. To confirm that the structural arrangement of the peptide side chains can affect anesthetic binding, gA in nonchannel forms in methanol was also analyzed. The addition of volatile anesthetic halothane to gA in SDS with a channel conformation caused a concentration-dependent change in resonant frequencies of the indole amide protons of W9, W11, W13, and W15, with the most profound changes in W9. These frequency changes were observed only for gA carefully prepared to ensure a channel conformation and were absent for gA in methanol. For gA in DMPC bilayers, direct [(14)C]halothane photolabeling and microsequencing demonstrated dominant labeling of W9, less labeling of W11 and W13, and no significant labeling of W15. In methanol, gA showed much less labeling of any residues. Inspection of the 3-D structure of gA suggests that the spatial arrangements of the tryptophan residues in the channel form of gA, combined with the amphiphilic regions of lipid, create a favorable anesthetic binding motif.

A single glycine residue at the entrance to the first membrane-spanning domain of the gamma-aminobutyric acid type A receptor beta(2) subunit affects allosteric sensitivity to GABA and anesthetics

Site-directed mutagenesis of the gamma-aminobutyric acid type A (GABA(A)) receptor beta(2) subunit has demonstrated that conversion of a conserved glycine residue located at the entrance to the first transmembrane domain into the homologous rho(1) residue phenylalanine alters the modulating effects of four different i.v. anesthetics: pentobarbital, alphaxalone, etomidate, and propofol. Using the baculovirus expression system in Spodoptera frugiperda 9 cells, anesthetic-induced enhancement of [(3)H]muscimol and [(3)H]flunitrazepam binding in receptors containing the beta(2)(G219F) point mutation displayed a significantly reduced efficacy in modulation by all four i.v. anesthetics tested. Furthermore, GABA(A) receptors containing the alpha(1)(G223F) point mutation also significantly decreased the maximal effect of etomidate- and propofol-induced enhancement of ligand binding. Conversely, the homologous point mutation in rho(1) receptors (F261G) changed the i.v. anesthetic-insensitive receptor to confer anesthetic modulation of [(3)H]muscimol binding. Consistent with the binding, functional analysis of pentobarbital-enhanced GABA currents recorded with whole-cell patch clamp demonstrated the beta(2)(G219F) subunit mutation eliminated the potentiating effect of the anesthetic. Similarly, propofol-enhanced GABA currents were potentiated less in alpha(1)beta(2)(G219F)gamma(2) receptors than in alpha(1)beta(2)gamma(2) receptors. Although ligand binding displayed comparable K(D) values for muscimol among wild-type, alpha(1)beta(2)gamma(2), and mutant receptors, patch-clamp recordings showed that alpha(1)beta(2)(G219F)gamma(2) receptors had a significantly more potent response to GABA than did alpha(1)beta(2)gamma(2) or alpha(1)(G223F)beta(2)gamma(2). The alpha(1)beta(2)(G219F)gamma(2) receptors also were more sensitive to direct channel activation by pentobarbital and propofol in the absence of GABA. These results suggest that the first transmembrane glycine residue on the beta(2) subunit may be important for conformational or allosteric interactions of channel gating by both GABA and anesthetics.

NMR study of volatile anesthetic binding to nicotinic acetylcholine receptors

New lines of evidence suggest that volatile anesthetics interact specifically with proteins. Direct binding analysis, however, has been largely limited to soluble proteins. In this study, specific interaction was investigated between isoflurane, a clinically important volatile anesthetic, and membrane-bound nicotinic acetylcholine receptors (nAChRs) from Torpedo electroplax, using (19)F nuclear magnetic resonance spectroscopy and gas chromatography. The receptors were reconstituted into 1, 2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) lipid vesicles. After correcting for nonspecific partitioning into the lipid, the equilibrium dissociation constant, K(d), of isoflurane binding to nAChR at 15 degrees C was found to be 0.36 +/- 0.03 mM. This value is within the clinically relevant concentration range of the agent. Based on the receptor concentrations in the vesicle suspension assayed by the bicinchoninic acid method and the fraction of bound isoflurane, X(b), determined by gas chromatography, an estimate of an average of 9-10 specifically bound isoflurane molecules can be made for each receptor, or two for each subunit. Upon binding, the transverse relaxation time constant (T(2)) of (19)F resonance of isoflurane is decreased by nearly three orders of magnitude, indicating a dramatic reduction in the mobility of specifically bound isoflurane. Kinetic analysis reveals that the off rate of binding, k(-1), is 1.7 x 10(4) s(-1). The on rate, k(+1), can thus be calculated to be approximately 4.8 x 10(7) M(-1) s(-1), suggesting a nearly diffusion-limited association. This is in contrast to anesthetic binding to a soluble protein, bovine serum albumin (BSA), where k(+1) and k(-1) are at least an order of magnitude slower. It is concluded that the presence of lipids may be critical for the correct evaluation of binding kinetics between volatile anesthetics and neuronal receptors.

Distinctly different interactions of anesthetic and nonimmobilizer with transmembrane channel peptides

Although it plays no clinical role in general anesthesia, gramicidin A, a transmembrane channel peptide, provides an excellent model for studying the specific interaction between volatile anesthetics and membrane proteins at the molecular level. We show here that a pair of structurally similar volatile anesthetic and nonimmobilizer (nonanesthetic), 1-chloro-1,2,2-trifluorocyclobutane (F3) and 1, 2-dichlorohexafluorocyclobutane (F6), respectively, interacts differently with the transmembrane peptide. With 400 microM gramicidin A in a vesicle suspension of 60 mM phosphatidylcholine-phosphatidylglycerol (PC/PG), the intermolecular cross-relaxation rate constants between (19)F of F3 and (1)H in the chemical shift regions for the indole and backbone amide protons were 0.0106 +/- 0.0007 (n = 12) and 0.0105 +/- 0.0014 (n = 8) s(-1), respectively. No cross-relaxation was measurable between (19)F of F6 and protons in these regions. Sodium transport study showed that with 75 microM gramicidin A in a vesicle suspension of 66 mM PC/PG, F3 increased the (23)Na apparent efflux rate constant from 149.7 +/- 7.2 of control (n = 3) to 191.7 +/- 12.2 s(-1) (n = 3), and the apparent influx rate constant from 182.1 +/- 15.4 to 222.8 +/- 21.7 s(-1) (n = 3). In contrast, F6 had no effects on either influx or efflux rate. It is concluded that the ability of general anesthetics to interact with amphipathic residues near the peptide-lipid-water interface and the inability of nonimmobilizer to do the same may represent some characteristics of anesthetic-protein interaction that are of importance to general anesthesia.

Structural consequences of anesthetic and nonimmobilizer interaction with gramicidin A channels

Although interactions of general anesthetics with soluble proteins have been studied, the specific interactions with membrane bound-proteins that characterize general anesthesia are largely unknown. The structural modulations of anesthetic interactions with synaptic ion channels have not been elucidated. Using gramicidin A as a simplified model for transmembrane ion channels, we have recently demonstrated that a pair of structurally similar volatile anesthetic and nonimmobilizer, 1-chloro-1,2,2-trifluorocyclobutane (F3) and 1,2-dichlorohexafluorocyclobutane (F6), respectively, have distinctly different effects on the channel function. Using high-resolution NMR structural analysis, we show here that neither F3 nor F6 at pharmacologically relevant concentrations can significantly affect the secondary structure of the gramicidin A channel. Although both the anesthetic F3 and the nonimmobilizer F6 can perturb residues at the middle section of the channel deep inside the hydrophobic region in the sodium dodecyl sulfate micelles, only F3, but not F6, can significantly alter the chemical shifts of the tryptophan indole N-H protons near the channel entrances. The results are consistent with the notion that anesthetics cause functional change of the channel by interacting with the amphipathic domains at the peptide-lipid-water interface.