Search

Smoke intoxication is a central event in mass burn incidents, and toxic smoke acts at different levels of the body, blocking breathing and oxygenation. The majority of these patients require early induction of anesthesia to preserve vital functions. We studied the influence of hemoglobin (HMG) and myoglobin (MGB) blockade by hydrochloric acid (HCl) in an interaction model with gaseous anesthetics using molecular docking techniques. In the next part of the study, molecular dynamics (MD) simulations were performed on the top-scoring ligand-receptor complexes to investigate the stability of the ligand-receptor complexes and the interactions between ligands and receptors in more detail. Through docking analysis, we observed that hemoglobin creates more stable complexes with anesthetic gases than myoglobin. Intoxication with gaseous hydrochloric acid produces conformational and binding energy changes of anesthetic gases to the substrate (both the pathway and the binding site), the most significant being recorded in the case of desflurane and sevoflurane, while for halothane and isoflurane, they remain unchanged. According to our theoretical model, the selection of anesthetic agents for patients affected by fire smoke containing hydrochloric acid is critical to ensure optimal anesthetic effects. In this regard, our model suggests that halothane and isoflurane are the most suitable choices for predicting the anesthetic effects in such patients when compared to sevoflurane and desflurane.

The Solubility of Halothane in Blood and Tissue Homogenates. By Larson CP, Eger EI, Severinghaus JW. Anesthesiology 1962; 23:349-55. Measured samples of human and bovine blood, human hemoglobin, and tissue homogenates from human fat and both human and bovine liver, kidney, muscle, whole brain, and separated gray and white cortex were added to stoppered 2,000-ml Erlenmeyer flasks. To each flask, 0.1 ml of liquid halothane was added under negative pressure using a calibrated micropipette. After the flask was agitated for 2 to 4 h to achieve equilibrium between the gas and blood or tissue contents, a calibrated infrared halothane analyzer was used to measure the concentration of halothane vapor. Calculated partition coefficients ranged from 0.7 for water to 2.3 for blood and from 3.5 for human or bovine kidney to 6 for human whole brain or liver and 8 for human muscle. Human peritoneal fat had a value of 138. The human blood-gas partition coefficient of 2.3 as determined by this equilibration method was well below the previously published value of 3.6., (Copyright © 2021, the American Society of Anesthesiologists. All Rights Reserved.)

Although instrumental for optimizing their pharmacological activity, a molecular understanding of the preferential interactions given by volatile anesthetics is quite poor. This paper confirms the ability of halothane to work as a hydrogen-bond (HB) donor and gives the first experimental proof that halothane also works as a halogen-bond (HaB) donor in the solid state and in solution. A halothane/hexamethylphosphortriamide co-crystal is described and its single-crystal X-ray structure shows short HaBs between bromine, or chlorine, and the phosphoryl oxygen. New UV/Vis absorption bands appear upon addition of diazabicyclooctane and tetra(n-butyl)ammonium iodide to halothane solutions, indicating that nitrogen atoms and anions may mediate the HaB-driven binding processes involving halothane as well. The ability of halothane to work as a bidentate/tridentate tecton by acting as a HaB and HB donor gives an atomic rationale for the eudismic ratio shown by this agent., (© 2019 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

We report results of a molecular dynamics simulation study of the effect of one general anesthetic, halothane, on some properties of mixed DPPC/DPPE phospholipid membranes. This is a suitable model for the study of simple, two-phospholipid membrane systems. From the simulation runs, we determined several membrane properties for five different molecular proportions of DPPC/DPPE. The effect of halothane on the studied membrane properties (area per lipid molecule, density of membrane, order parameter, etc.) was rather small. The distribution of halothane is not uniform through the bilayer thickness. Instead, there is a maximum of anesthetic concentration around 1.2 nm from the center of the membrane. The anesthetic molecule is located close to the phospholipid headgroups. The position of the halothane density maximum depends slightly on the DPPC/DPPE molar proportion. Snapshots taken over the plane of the membrane, as well as calculated two-dimensional radial distribution functions show that the anesthetic has no preference for either phospholipid (DPPC or DPPE). Our results indicate that this anesthetic molecule has only small effects on DPPC/DPPE mixed membranes. In addition, halothane displays no preferential location around DPPC or DPPE. This is probably due to the hydrophobic nature of halothane and to the fact that the chosen phospholipids have the same hydrophobic tails.

Anesthetics can interact with a wide variety of proteins in the body, including ion channels and alter their activity, but little is known about the molecular mechanisms of the interactions responsible for the functional activity. Characterization of the nature of anesthetic-protein interactions therefore is important and requires the complete analysis of the binding energetics. Isothermal titration calorimetry (ITC) is the only technique that allows quantitative determination of all thermodynamic parameters, including the equilibrium binding constant (K B ), the standard Gibbs free energy change (ΔG), the enthalpy change (ΔH), the entropy change (ΔS), heat capacity change (ΔC p ), and stoichiometry (n) of the reaction. ITC does not require any labeling or modification of the interacting partners analyzed and can be performed in solution with small amounts of reagents. In this chapter we describe the general properties of the ITC method, highlighting some critical aspects of experimental planning and data analysis, with practical application to anesthetic-protein interactions., (© 2018 Elsevier Inc. All rights reserved.)

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., (© 2018 Elsevier Inc. All rights reserved.)

The variation in phase-transition temperatures of dipalmitoylphosphatidylcholine (DPPC) bilayer membrane by adding two membrane-active ligands, a long-chain fatty acid (palmitic acid (PA)) and an inhalation anesthetic (halothane (HAL)), was investigated by light-transmittance measurements and fluorometry. By assuming the thermodynamic colligative property for the bilayer membrane at low ligand concentrations, the partitioning behavior of these ligands into the DPPC bilayer membrane was considered. It was proved from the differential partition coefficients between two phases that PA has strong affinity with the gel (lamellar gel) phase in a micro-molal concentration range and makes the bilayer membrane more ordered, while HAL has strong affinity with the liquid crystalline phase in a milli-molal concentration range and does the bilayer membrane more disordered. The transfer volumes of both ligands from the aqueous solution to each phase of the DPPC bilayer membrane showed that the preferential partitioning of the PA molecule into the gel (lamellar gel) produces about 20% decrease in transfer volume as compared with the liquid crystalline phase, whereas that of the HAL molecule into the liquid crystalline phase does about twice increase in transfer volume as compared with the gel (ripple gel) phase. Furthermore, changes in thermotropic and barotropic phase behavior of the DPPC bilayer membrane by adding the ligand was discussed from the viewpoint of the ligand partitioning. Reflecting the contrastive partitioning of PA and HAL into the pressure-induced interdigitated gel phase among the gel phases, it was revealed that PA suppresses the formation of the interdigitated gel phase under high pressure while HAL promotes it. These results clearly indicate that each phase of the DPPC bilayer membrane has a potential to recognize various ligand molecules., (Copyright © 2017 Elsevier B.V. All rights reserved.)

The effect of four general anesthetics, namely chloroform, halothane, diethyl ether, and enflurane on the free volume fraction and lateral pressure profiles in a fully hydrated dipalmitoylphosphatidylcholime (DPPC) membrane is investigated by means of computer simulation. In order to find changes that can be related to the molecular mechanism of anesthesia as well as its pressure reversal, the simulations are performed both at atmospheric and high (1000 bar) pressures. The obtained results show that the additional free volume occurring in the membrane is localized around the anesthetic molecules themselves. Correspondingly, the fraction of the free volume is increased in the outer of the two membrane regions (i.e., at the outer edge of the hydrocarbon phase) where anesthetic molecules prefer to stay in every case. As a consequence, the presence of anesthetics decreases the lateral pressure in the nearby region of the lipid chain ester groups, in which the anesthetic molecules themselves do not penetrate. Both of these changes, occurring upon introducing anesthetics in the membrane, are clearly reverted by the increase of the global pressure. These findings are in accordance both with the more than 60 years old "critical volume hypothesis" of Mullins, and with the more recent "lateral pressure hypothesis" of Cantor. Our results suggest that if anesthesia is indeed caused by conformational changes of certain membrane-bound proteins, induced by changes in the lateral pressure profile, as proposed by Cantor, the relevant conformational changes are expected to occur in the membrane region where the ester groups are located.

Neutron diffraction measurements demonstrate that hydrostatic pressure promotes liquid-ordered (Lo) domain formation in lipid membranes prepared as both oriented multilayers and unilamellar vesicles made of a canonical ternary lipid mixture for which demixing transitions have been extensively studied. The results demonstrate an unusually large dependence of the mixing transition on hydrostatic pressure. Additionally, data at 28 °C show that the magnitude of increase in Lo caused by 10 MPa pressure is much the same as the decrease in Lo produced by twice minimum alveolar concentrations (MAC) of general anesthetics such as halothane, nitrous oxide, and xenon. Therefore, the results may provide a plausible explanation for the reversal of general anesthesia by hydrostatic pressure.

We present the first set of ab initio calculations (vertical energies and oscillator strengths) of the valence and Rydberg transitions of the anaesthetic compound halothane (CF3CHBrCl). These results are complemented by high-resolution vacuum ultraviolet photoabsorption measurements over the wavelength range 115-310 nm (10.8-4.0 eV). The spectrum reveals several new features that were not previously reported in the literature. Spin-orbit effects have been considered in the calculations for the lowest-lying states, allowing us to explain the broad nature of the 6.1 and 7.5 eV absorption bands assigned to σ*(C-Br) ← nBr and σ*(C-Cl) ← n(Cl) transitions. Novel absolute photoabsorption cross sections from electron scattering data were derived in the 4.0-40.0 eV range. The measured absolute photoabsorption cross sections have been used to calculate the photolysis lifetime of halothane in the upper stratosphere (20-50 km).

Previous studies have demonstrated that volatile anesthetics could produce local anesthesia. Emulsified isoflurane at 8% has been reported to produce epidural anesthetic effect in rabbits. This study was designed to investigate the long-term epidural anesthetic effect of emulsified halothane in rabbits. In this study, 40 healthy adult rabbits (weighting 2.0-2.5 kg) with an epidural catheter were randomly divided into 4 groups (n=10/group), receiving epidural administration of 1% lidocaine (lido group), 8% emulsified isoflurane 1ml (8% E-iso group), 8% emulsified halothane (8% E-Halo group) and 12% emulsified halothane (12% E-Halo group). After administration, sensory and motor functions as well as consciousness state were assessed until 60 minutes after sensory and motor function returned to its baseline or at least for 180 min. After epidural anesthesia, all the rabbits were continuously observed for 7 days and sacrificed for pathological evaluations. As a result, all the four study solutions produced typical epidural anesthesia. Onset times of sensory and motor function blockade were similar among the four groups (P>0.05). Duration of sensory blockade in 12% E-Halo group (83±13 min) was significantly longer than other groups: 51±12 min in 8% E-Halo group (P<0.01), 57±8 min in 8% E-iso group (P<0.01) and 47±9 min in lido group (P<0.01). Duration of sensory blockade in 8% E-iso group is longer than lido group (P<0.05). Duration of motor blockade in 12% E-Halo group (81±12 min) was also significantly longer than other groups: 40±8 min in 8% E-Halo group (P<0.01), 37±3 min in 8% E-iso group (P<0.01), 37±6 min in lido group (P<0.01). Normal consciousness was found in the rabbits from 8% E-Halo, 8% E-iso and lido groups while there were four rabbits in 12% E-Halo group (4/10) showed a light sedation. For all the rabbits, no pathological injury was found. The present study demonstrates that emulsified halothane produces reversible concentration-dependent epidural anesthesia and at 12% (v/v), emulsified halothane could produce long-term anesthesia without pathological injury.

Background: The assessment of volatile agents' consumption can be performed by weighing vapourisers before and after use. This method is technically demanding and unavailable for retrospective analysis of anaesthesia records. Therefore, a method based on calculations from fresh gas flow and agent concentration is presented here., Methods: The presented calculation method herein enables a precise estimation of volatile agent consumption when average fresh gas flows and volatile agent concentrations are known. A pre-condition for these calculations is the knowledge of the vapour amount deriving from 1 ml fluid volatile agent. The necessary formulas for these calculations and an example for a sevoflurane anaesthesia are presented., Results: The amount of volatile agent vapour deriving from 1 ml of fluid agent are for halothane 229 ml, isoflurane 195 ml, sevoflurane 184 m, and desflurane 210 ml. The constant for sevoflurane is used in a fictitious clinical case to exemplify the calculation of its consumption in daily routine resulting in a total expenditure of 23.6 ml liquid agent., Conclusions: By application of the presented specific volatile agent constants and equations, it becomes easy to calculate volatile agent consumption if the fresh gas flows and the resulting inhaled concentration of the volatile agent are known. By this method, it is possible to extract data about volatile agent consumption both ways: (1) retrospectively from sufficiently detailed and accurate anaesthesia recordings, as well as (2) by application of this method in a prospective setting. Therefore, this method is a valuable contribution to perform pharmacoeconomical surveys., (© 2014 The Acta Anaesthesiologica Scandinavica Foundation. Published by John Wiley & Sons Ltd.)

The α7 nicotinic acetylcholine receptor (nAChR), assembled as homomeric pentameric ligand-gated ion channels, is one of the most abundant nAChR subtypes in the brain. Despite its importance in memory, learning and cognition, no structure has been determined for the α7 nAChR TM domain, a target for allosteric modulators. Using solution state NMR, we determined the structure of the human α7 nAChR TM domain (PDB ID: 2MAW) and demonstrated that the α7 TM domain formed functional channels in Xenopus oocytes. We identified the associated binding sites for the anesthetics halothane and ketamine; the former cannot sensitively inhibit α7 function, but the latter can. The α7 TM domain folds into the expected four-helical bundle motif, but the intra-subunit cavity at the extracellular end of the α7 TM domain is smaller than the equivalent cavity in the α4β2 nAChRs (PDB IDs: 2LLY; 2LM2). Neither drug binds to the extracellular end of the α7 TM domain, but two halothane molecules or one ketamine molecule binds to the intracellular end of the α7 TM domain. Halothane and ketamine binding sites are partially overlapped. Ketamine, but not halothane, perturbed the α7 channel-gate residue L9'. Furthermore, halothane did not induce profound dynamics changes in the α7 channel as observed in α4β2. The study offers a novel high-resolution structure for the human α7 nAChR TM domain that is invaluable for developing α7-specific therapeutics. It also provides evidence to support the hypothesis: only when anesthetic binding perturbs the channel pore or alters the channel motion, can binding generate functional consequences., (Copyright © 2013 Elsevier B.V. All rights reserved.)

The advent of the polymerase chain reaction and the availability of data from various global human genome projects should make it possible, using a DNA sample isolated from white blood cells, to diagnose rapidly and accurately almost any monogenic condition resulting from single nucleotide changes. DNA-based diagnosis for malignant hyperthermia (MH) is an attractive proposition, because it could replace the invasive and morbid caffeine-halothane/in vitro contracture tests of skeletal muscle biopsy tissue. Moreover, MH is preventable if an accurate diagnosis of susceptibility can be made before general anesthesia, the most common trigger of an MH episode. Diagnosis of MH using DNA was suggested as early as 1990 when the skeletal muscle ryanodine receptor gene (RYR1), and a single point mutation therein, was linked to MH susceptibility. In 1994, a single point mutation in the α 1 subunit of the dihydropyridine receptor gene (CACNA1S) was identified and also subsequently shown to be causative of MH. In the succeeding years, the number of identified mutations in RYR1 has grown, as has the number of potential susceptibility loci, although no other gene has yet been definitively associated with MH. In addition, it has become clear that MH is associated with either of these 2 genes (RYR1 and CACNA1S) in only 50% to 70% of affected families. While DNA testing for MH susceptibility has now become widespread, it still does not replace the in vitro contracture tests. Whole exome sequence analysis makes it potentially possible to identify all variants within human coding regions, but the complexity of the genome, the heterogeneity of MH, the limitations of bioinformatic tools, and the lack of precise genotype/phenotype correlations are all confounding factors. In addition, the requirement for demonstration of causality, by in vitro functional analysis, of any familial mutation currently precludes DNA-based diagnosis as the sole test for MH susceptibility. Nevertheless, familial DNA testing for MH susceptibility is now widespread although limited to a positive diagnosis and to those few mutations that have been functionally characterized. Identification of new susceptibility genes remains elusive. When new genes are identified, it will be the role of the biochemists, physiologists, and biophysicists to devise functional assays in appropriate systems. This will remain the bottleneck unless high throughput platforms can be designed for functional work. Analysis of entire genomes from several individuals simultaneously is a reality. DNA testing for MH, based on current criteria, remains the dream.

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., (Copyright © 2012 Elsevier B.V. All rights reserved.)

The effect of four general anesthetic molecules, i.e., chloroform, halothane, diethyl ether and enflurane, on the properties of a fully hydrated dipalmitoylphosphatidylcholine (DPPC) membrane is studied in detail by long molecular dynamics simulations. Furthermore, to address the problem of pressure reversal, the effect of pressure on the anesthetic containing membranes is also investigated. In order to ensure sufficient equilibration and adequate sampling, the simulations performed have been at least an order of magnitude longer than the studies reported previously in the literature on general anesthetics. The results obtained can help in resolving several long-standing contradictions concerning the effect of anesthetics, some of which were the consequence of too short simulation time used in several previous studies. More importantly, a number of seeming contradictions are found to originate from the fact that different anesthetic molecules affect the membrane structure differently in several respects. In particular, halothane, being able to weakly hydrogen bound to the ester group of the lipid tails, is found to behave in a markedly different way than the other three molecules considered. Besides, we also found that two changes, namely lateral expansion of the membrane and increasing local disorder in the lipid tails next to the anesthetic molecules, are clearly induced by all four anesthetic molecules tested here in the same way, and both of these effects are reverted by the increase in pressure.

The first quantitative insight has been obtained into the effects that volatile anesthetics have on the interactions and lateral organization of lipids in model membranes that mimic "lipid rafts". Specifically, nearest-neighbor recogntion measurements, in combination with Monte Carlo simulations, have been used to investigate the action of isoflurane, halothane, and chloroform on the compactness and lateral organization of cholesterol-rich bilayers of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) in the liquid-ordered (l(o)) phase. All three anesthetics induce a similar weakening of sterol-phospholipid association, corresponding to ca. 30 cal/mol of lipid at clinically relevant concentrations. Monte Carlo lattice simulations show that the lateral organization of the l(o) phase, under such conditions, remains virtually unchanged. In sharp contrast to their action on the l(o) phase, these anesthetics have been found to have a similar strengthening effect on sterol-phospholipid association in the liquid-disordered (l(d)) phase. The possibility of discrete complexes being formed between DPPC and these anesthetics and the biological relevance of these findings are discussed., (© 2011 American Chemical Society)

The structures and intermolecular interactions in the halogen bonded complexes of anaesthetics (chloroform, halothane, enflurane and isoflurane) with formaldehyde were studied by ab initio MP2 and CCSD(T) methods. The CCSD(T)/CBS calculated binding energies of these complexes are between -2.83 and -4.21 kcal mol(-1). The largest stabilization energy has been found for the C-Br···O bonded halothane···OCH(2) complex. In all complexes the C-X bond length (where X = Cl, Br) is slightly shortened, in comparison to a free compound, and an increase of the C-X stretching frequency is observed. The electrostatic interaction was excluded as being responsible for the C-X bond contraction. It is suggested that contraction of the C-X bond length can be explained in terms of the Pauli repulsion (the exchange overlap) between the electron pairs of oxygen and halogen atoms in the investigated complexes. This is supported by the DFT-SAPT results, which indicate that the repulsive exchange energy overcompensates the electrostatic one. Moreover, the dispersion and electrostatic contributions cover about 95% of the total attraction forces, in these complexes.

G protein-activated inwardly rectifying potassium (GIRK or Kir3) channels are directly gated by the βγ subunits of G proteins and contribute to inhibitory neurotransmitter signaling pathways. Paradoxically, volatile anesthetics such as halothane inhibit these channels. We find that neuronal Kir3 currents are highly sensitive to inhibition by halothane. Given that Kir3 currents result from increased Gβγ available to the channels, we asked whether reducing available Gβγ to the channel would adversely affect halothane inhibition. Remarkably, scavenging Gβγ using the C-terminal domain of β-adrenergic receptor kinase (cβARK) resulted in channel activation by halothane. Consistent with this effect, channel mutants that impair Gβγ activation were also activated by halothane. A single residue, phenylalanine 192, occupies the putative Gβγ gate of neuronal Kir3.2 channels. Mutation of Phe-192 at the gate to other residues rendered the channel non-responsive, either activated or inhibited by halothane. These data indicated that halothane predominantly interferes with Gβγ-mediated Kir3 currents, such as those functioning during inhibitory synaptic activity. Our report identifies the molecular correlate for anesthetic inhibition of Kir3 channels and highlights the significance of these effects in modulating neurotransmitter-mediated inhibitory signaling.

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

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., (Copyright © 2010 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

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.

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.

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., (Copyright 2009 Elsevier B.V. All rights reserved.)

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.

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.

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.

We find that the sensitivity of gramicidin A channels to the anesthetic halothane is highly lipid dependent. Specifically, exposure of membranes made of lamellar DOPC to halothane in concentrations close to clinically relevant reduces channel lifetimes by 1 order of magnitude. At the same time, gramicidin channels in membranes of nonlamellar DOPE are affected little, if at all, by halothane. We attribute this difference in channel behavior to a difference in the stress of lipid packing into a planar lipid bilayer, wherein the higher stress of DOPE packing reduces the degree of halothane partitioning into the hydrophobic interior.

A nitrile-derived amino acid, Phe(CN), has been used as an internal spectroscopic probe to study the binding of an inhalational anesthetic to a model membrane protein. The infrared spectra from experiment showed a blue-shift of the nitrile vibrational frequency in the presence of the anesthetic halothane. To interpret the infrared results and explore the nature of the interaction between halothane and the model protein, all-atom molecular dynamics (MD) simulations have been used to probe the structural and dynamic properties of the protein in the presence and absence of one halothane molecule. The frequency shift analyzed from MD simulations agrees well with the experimental infrared results. Decomposition of the forces acting on the nitrile probes demonstrates an indirect impact on the probes from halothane, namely a change of the protein's electrostatic local environment around the probes induced by halothane. Although the halothane remains localized within the designed hydrophobic binding cavity, it undergoes a significant amount of translational and rotational motion, modulated by the interaction of the trifluorine end of halothane with backbone hydrogens of the residues forming the cavity. This dominant interaction between halothane and backbone hydrogens outweighs the direct interaction between halothane and the nitrile groups, making it a good "spectator" probe of the halothane-protein interaction. These MD simulations provide insight into action of anesthetic molecules on the model membrane protein, and also support the further development of nitrile-labeled amino acids as spectroscopic probes within the designed binding cavity.

We previously reported the synthesis and structural characterization of a model membrane protein comprised of an amphiphilic 4-helix bundle peptide with a hydrophobic domain based on a synthetic ion channel and a hydrophilic domain with designed cavities for binding the general anesthetic halothane. In this work, we synthesized an improved version of this halothane-binding amphiphilic peptide with only a single cavity and an otherwise identical control peptide with no such cavity, and applied x-ray reflectivity to monolayers of these peptides to probe the distribution of halothane along the length of the core of the 4-helix bundle as a function of the concentration of halothane. At the moderate concentrations achieved in this study, approximately three molecules of halothane were found to be localized within a broad symmetric unimodal distribution centered about the designed cavity. At the lowest concentration achieved, of approximately one molecule per bundle, the halothane distribution became narrower and more peaked due to a component of approximately 19A width centered about the designed cavity. At higher concentrations, approximately six to seven molecules were found to be uniformly distributed along the length of the bundle, corresponding to approximately one molecule per heptad. Monolayers of the control peptide showed only the latter behavior, namely a uniform distribution along the length of the bundle irrespective of the halothane concentration over this range. The results provide insight into the nature of such weak binding when the dissociation constant is in the mM regime, relevant for clinical applications of anesthesia. They also demonstrate the suitability of both the model system and the experimental technique for additional work on the mechanism of general anesthesia, some of it presented in the companion parts II and III under this title.

We demonstrate that cyano-phenylalanine (Phe(CN)) can be utilized to probe the binding of the inhalational anesthetic halothane to an anesthetic-binding, model ion channel protein hbAP-Phe(CN). The Trp to Phe(CN) mutation alters neither the alpha-helical conformation nor the 4-helix bundle structure. The halothane binding properties of this Phe(CN) mutant hbAP-Phe(CN), based on fluorescence quenching, are consistent with those of the prototype, hbAP1. The dependence of fluorescence lifetime as a function of halothane concentration implies that the diffusion of halothane in the nonpolar core of the protein bundle is one-dimensional. As a consequence, at low halothane concentrations, the quenching of the fluorescence is dynamic, whereas at high concentrations the quenching becomes static. The 4-helix bundle structure present in aqueous detergent solution and at the air-water interface, is preserved in multilayer films of hbAP-Phe(CN), enabling vibrational spectroscopy of both the protein and its nitrile label (-CN). The nitrile groups' stretching vibration band shifts to higher frequency in the presence of halothane, and this blue-shift is largely reversible. Due to the complexity of this amphiphilic 4-helix bundle model membrane protein, where four Phe(CN) probes are present adjacent to the designed cavity forming the binding site within each bundle, all contributing to the infrared absorption, molecular dynamics (MD) simulation is required to interpret the infrared results. The MD simulations indicate that the blue-shift of -CN stretching vibration induced by halothane arises from an indirect effect, namely an induced change in the electrostatic protein environment averaged over the four probe oscillators, rather than a direct interaction with the oscillators. hbAP-Phe(CN) therefore provides a successful template for extending these investigations of the interactions of halothane with the model membrane protein via vibrational spectroscopy, using cyano-alanine residues to form the anesthetic binding cavity.

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.

Abeta peptide is the major component of senile plaques (SP) which accumulates in AD (Alzheimer's disease) brain. Reports from different laboratories indicate that anesthetics interact with Abeta peptide and induce Abeta oligomerization. The molecular mechanism of Abeta peptide interactions with these anesthetics was not determined. We report molecular details for the interactions of uniformly (15)N labeled Abeta40 with different anesthetics using 2D nuclear magnetic resonance (NMR) experiments. At high concentrations both isoflurane and propofol perturb critical amino acid residues (G29, A30 and I31) of Abeta peptide located in the hinge region leading to Abeta oligomerization. In contrast, these three specific residues do not interact with thiopental and subsequently no Abeta oligomerization was observed. However, studies with combined anesthetics (thiopental and halothane), showed perturbation of these residues (G29, A30 and I31) and subsequently Abeta oligomerization was found. Perturbation of these specific Abeta residues (G29, A30 and I31) by different anesthetics could play an important role to induce Abeta oligomerization.

A variety of experiments suggest that membrane proteins are important targets of anesthetic molecules, and that ion channels interact differently with anesthetics in their open and closed conformations. The availability of an open and a closed structural model for the KirBac1.1 potassium channel has made it possible to perform a comparative analysis of the interactions of anesthetics with the same channel in its open and closed states. To this end, all-atom molecular dynamics simulations supplemented by normal mode analysis have been employed to probe the interactions of the inhalational anesthetic halothane with both an open and closed conformer of KirBac1.1 embedded in a lipid bilayer. Normal mode analysis on the closed and open channel, in the presence and absence of halothane, reveals that the anesthetic modulates the global as well as the local dynamics of both conformations differently. In the case of the open channel, the observed reduction of flexibility of residues in the inner helices suggests a functional modification action of anesthetics on ion channels. In this context, preferential quenching of the aromatic residue motion and modulation of global dynamics by halothane may be seen as steps toward potentiating or favoring open state conformations. These molecular dynamics simulations provide the first insights into possible specific interactions between anesthetic molecules and ion channels in different conformations.

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.

The complex of halothane (CHClBrCF(3)) and dimethyl ether has been investigated experimentally in solutions of liquid krypton using infrared spectroscopy and theoretically using ab initio calculations at the MP2/6-311++G(d,p) level. The formation of a 1:1 complex was experimentally detected. The most stable ab initio geometry found is the one in which the C--H bond of halothane interacts with the oxygen atom of dimethyl ether. The complexes in which the chlorine or the bromine atom of halothane interacts with the oxygen atom of the ether were found to be local energy minima and were less stable by 14.5 and 9.3 kJ mol(-1), respectively, than the global minimum. The formation of a single complex species was observed in the infrared spectra; the standard complexation enthalpy of this complex was determined to be -12.3(8) kJ mol(-1). Analysis of the observed complexation shifts supports the identification of the complex as the hydrogen-bonded species. The C--H stretching vibration of halothane was found to show a redshift upon complexation of 19(2) cm(-1). The infrared intensity ratios epsilon(complex)/epsilon(monomer) for the fundamental and its first overtone were measured to be 6.5(1) and 0.31(1). The frequency shift was analyzed using Morokuma-type analysis, and the infrared intensity ratios were rationalized using a model including the mechanical and electric anharmonicity of the C--H stretching fundamental.

Background: Halothane, isoflurane, and sevoflurane exert negative inotropic side effects, generally mediated via a reduced availability of intracellular calcium. Other possible mechanisms include modified intracellular calcium handling, impaired actomyosin cross-bridge cycling, and/or alteration of calcium-induced conformational changes of the regulatory troponin complex., Methods: We investigated the effect of halothane, isoflurane, and sevoflurane on calcium-dependent kinetics of isolated human recombinant cardiac troponin C labeled with IAANS (HrcTnC(IAANS)) using stopped-flow and calcium titration techniques., Results: Calcium concentration at half-maximal fluorescence intensity (K(d)) in the control group was 2.1 +/- 0.1 mM. Volatile anesthetics increased calcium sensitivity in a concentration-dependent fashion sevoflurane (K(d) 1.5-1.7 mM, P = 0.001) > halothane (K(d) 1.7-1.9 mM, P < 0.01) > isoflurane (K(d) 1.8-1.9 mM, P < 0.05). The rate constant of conformational changes after rapid dissociation of calcium from HrcTnC(IAANS) (k(off(c))) was moderately prolonged at 4 degrees C by halothane and isoflurane > sevoflurane., Conclusion: These mechanisms may counteract the effects of lower calcium availability, and can be responsible for abbreviated, and possibly incomplete, relaxation of cardiac muscle fibers in the presence of volatile anesthetics.

Computational methods designed to predict and visualize ligand protein binding interactions were used to characterize volatile anesthetic (VA) binding sites and unoccupied pockets within the known structures of VAs bound to serum albumin, luciferase, and apoferritin. We found that both the number of protein atoms and methyl hydrogen, which are within approximately 8 A of a potential ligand binding site, are significantly greater in protein pockets where VAs bind. This computational approach was applied to structures of calmodulin (CaM), which have not been determined in complex with a VA. It predicted that VAs bind to [Ca(2+)](4)-CaM, but not to apo-CaM, which we confirmed with isothermal titration calorimetry. The VA binding sites predicted for the structures of [Ca(2+)](4)-CaM are located in hydrophobic pockets that form when the Ca(2+) binding sites in CaM are saturated. The binding of VAs to these hydrophobic pockets is supported by evidence that halothane predominantly makes contact with aliphatic resonances in [Ca(2+)](4)-CaM (nuclear Overhauser effect) and increases the Ca(2+) affinity of CaM (fluorescence spectroscopy). Our computational analysis and experiments indicate that binding of VA to proteins is consistent with the hydrophobic effect and the Meyer-Overton rule.

Molecular dynamics simulations have been performed to investigate the partitioning of the volatile anesthetic halothane from an aqueous phase into a coexisting hydrated bilayer, composed of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) lipids, with embedded alpha-helical peptide bundles based on the membrane-bound portions of the alpha- and delta-subunits, respectively, of nicotinic acetylcholine receptor. In the molecular dynamics simulations halothane molecules spontaneously partitioned into the DOPC bilayer and then preferentially occupied regions close to lipid headgroups. A single halothane molecule was observed to bind to tyrosine (Tyr-277) residue in the alpha-subunit, an experimentally identified specific binding site. The binding of halothane attenuated the local loop dynamics of alpha-subunit and significantly influenced global concerted motions suggesting anesthetic action in modulating protein function. Steered molecular dynamics calculations on a single halothane molecule partitioned into a DOPC lipid bilayer were performed to probe the free energy profile of halothane across the lipid-water interface and rationalize the observed spontaneous partitioning. Partitioned halothane molecules affect the hydrocarbon chains of the DOPC lipid, by lowering of the hydrocarbon tilt angles. The anesthetic molecules also caused a decrease in the number of peptide-lipid contacts. The observed local and global effects of anesthetic binding on protein motions demonstrated in this study may underlie the mechanism of action of anesthetics at a molecular level.

Inhaled anesthetic molecule occupancy of a protein internal cavity depends in part on the volumes of the guest molecule and the host site. Current algorithms to determine volume and surface area of cavities in proteins whose structures have been determined and cataloged make no allowance for shape or small degrees of shape adjustment to accommodate a guest. We developed an algorithm to determine spheroid dimensions matching cavity volume and surface area and applied it to screen the cavities of 6,658 nonredundant structures stored in the Protein Data Bank (PDB) for potential targets of halothane (2-bromo-2-chloro-1,1,1-trifluoroethane). Our algorithm determined sizes of prolate and oblate spheroids matching dimensions of each cavity found. If those spheroids could accommodate halothane (radius 2.91 A) as a guest, we determined the packing coefficient. 394,766 total cavities were identified. Of 58,681 cavities satisfying the fit criteria for halothane, 11,902 cavities had packing coefficients in the range of 0.46-0.64. This represents 20.3% of cavities large enough to hold halothane, 3.0% of all cavities processed, and found in 2,432 protein structures. Our algorithm incorporates shape dependence to screen guest-host relationships for potential small molecule occupancy of protein cavities. Proteins with large numbers of such cavities are more likely to be functionally altered by halothane.

In spite of numerous investigations, the molecular mechanism of general anesthetics action is still not well understood. It has been shown that the anesthetic potency is related to the ability of an anesthetic to partition into the membrane. We have investigated changes in structure, dynamics and forces of interaction in supported dipalmitoylphosphatidylcholine (DPPC) bilayers in the presence of the general anesthetic halothane. In the present study, we measured the forces of interaction between the probe and the bilayer using an atomic force microscope. The changes in force curves as a function of anesthetic incorporation were analyzed. Force measurements were in good agreement with AFM imaging data, and provided valuable information on bilayer thickness, structural transitions, and halothane-induced changes in electrostatic and adhesive properties.

Inhaled anesthetics are thought to alter the conformational states of Cys-loop ligand-gated ion channels (LGICs) by binding within discrete cavities that are lined by portions of four alpha-helical transmembrane domains. Because Cys-loop LGICs are complex molecules that are notoriously difficult to express and purify, scaled-down models have been used to better understand the basic molecular mechanisms of anesthetic action. In this study, stopped-flow fluorescence spectroscopy was used to define the kinetics with which inhaled anesthetics interact with (Aalpha(2)-L1M/L38M)(2), a four-alpha-helix bundle protein that was designed to model anesthetic binding sites on Cys-loop LGICs. Stopped-flow fluorescence traces obtained upon mixing (Aalpha(2)-L1M/L38M)(2) with halothane revealed immediate, fast, and slow components of quenching. The immediate component, which occurred within the mixing time of the spectrofluorimeter, was attributed to direct quenching of tryptophan fluorescence upon halothane binding to (Aalpha(2)-L1M/L38M)(2). This was followed by a biexponential fluorescence decay containing fast and slow components, reflecting anesthetic-induced conformational transitions. Fluorescence traces obtained in studies using sevoflurane, isoflurane, and desflurane, which poorly quench tryptophan fluorescence, did not contain the immediate component. However, these anesthetics did produce the fast and slow components, indicating that they also alter the conformation of (Aalpha(2)-L1M/L38M)(2). Cyclopropane, an anesthetic that acts with unusually low potency on Cys-loop LGICs, acted with low apparent potency on (Aalpha(2)-L1M/L38M)(2). These results suggest that four-alpha-helix bundle proteins may be useful models of in vivo sites of action that allow the use of a wide range of techniques to better understand how anesthetic binding leads to changes in protein structure and function.

The nature and the sites of interactions between anesthetic halothane and homodimeric Delta5-3-ketosteroid isomerase (KSI) are characterized by flexible ligand docking and confirmed by 1H-15N NMR. The dynamics consequence of halothane interaction and the implication of the dynamic changes to KSI function are studied by multiple 5-ns molecular dynamics simulations in the presence and absence of halothane. Both docking and MD simulations show that halothane prefer the amphiphilic dimeric interface to the hydrophobic active site of KSI. Halothane occupancy at the dimer interface disrupted the intersubunit hydrogen bonding formed either directly through side chains of polar residues or indirectly through the mediation of the interfacial water molecules. Moreover, in the presence of halothane, the exchange rate of the bound waters with bulk water was increased. Halothane perturbation to the dimer interface affected the overall flexibility of the active site. This action is likely to contribute to the halothane-induced reduction of the KSI activity. The allosteric halothane modulation of the dynamics-function relationship of KSI without direct competition at the enzymatic active sites may be generalized to offer a unifying explanation of anesthetic action on a diverse range of multidomain neuronal proteins that are potentially relevant to clinical general anesthesia.

Multiple binding sites for anaesthetics in HSA (human serum albumin) make solution studies difficult to interpret. In the present study, we expressed the wild-type HSA domain 3 (wtHSAd3), a peptide with two known anaesthetic binding sites in a yeast expression system. We also expressed a site-directed mutant of domain 3 (Y411Wd3). The stability and secondary structure of the constructed fragments were determined by HX (hydrogen-tritium exchange) and CD spectroscopy. The binding of two general anaesthetics, 2-bromo-2-chloro-1,1,1-trifluoroethane and propofol, to wtHSAd3 and Y411Wd3 was determined using isothermal titration calorimetry, HX and intrinsic tryptophan fluorescence quenching. Although the expressed fragments are less stable than intact wtHSA as indicated by both CD and HX, they retain the secondary structure and anaesthetic-binding characteristics of an intact HSA molecule, but with fewer binding sites. Y411Wd3 had decreased affinity for propofol but not for 2-bromo-2-chloro-1,1,1-trifluoroethane, consistent with steric hindrance. Retention of structural features and anaesthetic binding properties with fewer binding sites in this truncated protein provide feasibility for using scaled-down models of otherwise intractable systems to gain an understanding of anaesthetic binding requirements and binding-stability relationships.

Physiologic sites for inhaled anesthetics are presumed to be cavities within transmembrane 4-alpha-helix bundles of neurotransmitter receptors, but confirmation of binding and structural detail of such sites remains elusive. To provide such detail, we screened soluble proteins containing this structural motif, and found only one that exhibited evidence of strong anesthetic binding. Ferritin is a 24-mer of 4-alpha-helix bundles; both halothane and isoflurane bind with K(A) values of approximately 10(5) M(-1), higher than any previously reported inhaled anesthetic-protein interaction. The crystal structures of the halothane/apoferritin and isoflurane/apoferritin complexes were determined at 1.75 A resolution, revealing a common anesthetic binding pocket within an interhelical dimerization interface. The high affinity is explained by several weak polar contacts and an optimal host/guest packing relationship. Neither the acidic protons nor ether oxygen of the anesthetics contribute to the binding interaction. Compared with unliganded apoferritin, the anesthetic produced no detectable alteration of structure or B factors. The remarkably high affinity of the anesthetic/apoferritin complex implies greater selectivity of protein sites than previously thought, and suggests that direct protein actions may underlie effects at lower than surgical levels of anesthetic, including loss of awareness.

To gain insights into the molecular level mechanism of drug action at the membrane site, we have carried out extensive molecular dynamics simulations of a model membrane in the presence of a volatile anesthetic using a coarse-grain model. Six different anesthetic (halothane)/lipid (dimyristoylphosphatidylcholine) ratios have been investigated, going beyond the low doses typical of medical applications. The volatile anesthetics were introduced into a preassembled fully hydrated 512-molecule lipid bilayer and each of the molecular dynamics simulations were carried out at ambient conditions, using the NPT ensemble. The area per lipid increases monotonically with the halothane concentration and the lamellar spacing decreases, whereas the lipid bilayer thickness shows no appreciable differences and only a slight increase upon addition of halothane. The density profiles of the anesthetic molecules display a bimodal distribution along the membrane normal with maxima located close to the lipid-water interface region. We have studied how halothane molecules fluctuate between the two maxima of the bimodal distribution and we observed a different mechanism at low and high anesthetic concentrations. Through the investigation of the reorientational motions of the lipid tails, we found that the anesthetic molecules increase the segmental order of the lipids close to the membrane surface.

The ion chemistry of 1-bromo-1-chloro-2,2,2-trifluoroethane (the common anesthetic halothane) in air plasma at atmospheric pressure was investigated by atmospheric pressure chemical ionization mass spectrometry (APCI-MS). The major positive ion observed at low declustering (API interface) energies is the ionized dimer, M(+.)M, an unexpectedly abundant species which possibly is stabilized by two H-bonding interactions. At higher energies [M--HF](+.) and [M--Br](+) prevail; the former, corresponding to ionized olefin [ClBrC=CF(2)](+.), appears to originate from M(+.)M and is quite stable towards fragmentation. The latter fragment ion ([M--Br](+)) and its analogue, [M--Cl](+), which is also observed though at much lower abundance, are originally ethyl cations (+)CHX--CF(3) (X = Br, Cl) which, upon collisional activation, rearrange and fragment to CHFX(+) via elimination of CF(2). All of the above described ions are also observed in humid air: in addition, the oxygenated ion [ClBrC=CFOH](+.) also forms in humid air via water addition to [ClBrC=CF(2)](+.) and HF elimination, as observed earlier for ionized trichloroethene. In contrast with similar chloro- and fluoro-substituted ethanes, halothane does not react with H(3)O(+) in the APCI plasma, a result confirmed by selected ion APCI triple-quadrupole (TQ) experiments. Major negative ions formed from halothane in the air plasma are Br(-) and, to a lesser extent, Cl(-), and their complexes with neutral halothane. APCI-TQ experiments indicated that Br(-) and Cl(-) are formed via reaction of halothane with O(2) (-.), O(2) (-.)(H(2)O) and O(3) (-.), possibly via dissociative electron transfer or nucleophilic substitution. Competing proton transfer was also observed in the reaction with O(2) (-.) and, at high halothane pressure, also with O(2) (-.)(H(2)O); at lower pressures the molecular anion M(-.) was observed instead. The other minor anions of the air plasma, NO(2) (-), N(2)O(2) (-.) and NO(3) (-), were found to be unreactive towards halothane., (Copyright (c) 2005 John Wiley & Sons, Ltd.)

To determine if occupancy of interfacial pockets in oligomeric proteins by volatile anesthetic molecules can allosterically regulate oligomerization equilibria, variants of a three-helix bundle peptide able to form higher oligomers were studied with analytical ultracentrifugation, hydrogen exchange and modeling. Halothane shifted the oligomerization equilibria towards the oligomer only in a mutation predicted to create sufficient volume in the hexameric pocket. Other mutations at this residue, predicted to create a too small or too polar pocket, were unaffected by halothane. Inhaled anesthetic modulation of oligomerization interactions is a novel and potentially generalizable biophysical basis for some anesthetic actions.