Anesthetic modulation of protein dynamics: insight from an NMR study

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

Ethanol and halothane differently modulate HLA class I and class II oligomerization. A new look at the mode of action of anesthetic agents through fluorescence spectroscopy

The field of research considering the working mechanism of anesthetic agents is a complex one and the site or sites of action of general anesthetics are yet to be elucidated. Through the years, on the molecular level, the discussion has shifted from the lipid theories to the more specific interaction with the proteins responsible for the signal transduction. While this approach led to several models, they offer, at best, partial explanations for the observed phenomena. Anesthetic agents interact with many systems, of which the neuronal is best studied, leaving interaction with the immune defense system relatively unexplored. In this study we focus on the interaction of ethanol and halothane with the co-localization on the membrane of HLA I and II molecules. We show that ethanol tends to randomize the distribution of HLA I and II molecules, while halothane increases the clustering of HLA I proteins. The notion that anesthetics modulate cell function by disrupting clustering and thereby promoting a random distribution is a novel approach that may explain the general involvement of many systems during exposition to anesthetic drugs. In this study we show the disturbance of co-localization of molecules that may form a functional network. The relevance of this finding depends on the importance of these networks for extracellular and intracellular processes.     

A model for binding of inhalation anesthetics to membranes: two dose-dependent distinctly different binding modes

Inhalation anesthetics currently in clinical use, such as halothane, methoxyflurane, enflurane, isoflurane, etc., are polar hydrophobic molecules, except nitrous oxide, which is an apolar and weak anesthetic, incapable of inducing surgical stage anesthesia. Experimental data are accumulating that these potent amphipathic inhalation anesthetics preferentially bind membranes and macromolecules on the surface at clinical concentrations. The anesthetic binding to lipid membranes in the low concentration range is characterized by a saturable curve approaching to a limiting value. When the anesthetic concentration s greatly increased above the clinical range, the binding starts to exceed the limiting saturation value. Our model for anesthetic binding to membranes consists of two parts: Langmuir-type adsorption to the membrane surface at the low concentration range and penetration into the hydrophobic core at the high concentration range. The present communication provides a statistical-thermodynamic basis to analyze this twostep interaction. An expression is derived for membrane capacitance as a function of anesthetic concentration, which explains the experimental data well. Binding parameters of anesthetics are estimated according to the theory.This study was supported by NIH grants GM 25716 and GM 26950, and by the Medical Research Service of the Veterans Administration.     

Temperature dependence of thermodynamic activity in volatile anesthetics: correlation between anesthetic potency and activity

Temperature dependence of the saturated concentration and the activity coefficient of anesthetics (1-propanol, diethyl ether, chloroform, and halothane) in water were evaluated using vapor pressure and H NMR measurement. We found that these physical values (quantities) correlate with anesthetic potencies estimated according to the thermodynamic equilibrium model. The anesthetic potency for hydrophilic anesthetic (diethyl ether) decreased with decreasing temperature because of the temperature specificity of this saturated concentration. In contrast, potencies of hydrophobic anesthetics (chloroform and halothane) increased with decreasing temperature because of the temperature specificity of those activity coefficients. By assuming that anesthetics interact with hydrated water of cell membranes, the temperature dependence of anesthetic potencies in vivo is qualitatively explicable.     

Halothane changes the domain structure of a binary lipid membrane

X-ray and neutron diffraction studies of a binary lipid membrane demonstrate that halothane at physiological concentrations produces a pronounced redistribution of lipids between domains of different lipid types identified by different lamellar d-spacings and isotope composition. In contrast, dichlorohexafluorocyclobutane (F6), a halogenated nonanesthetic, does not produce such significant effects. These findings demonstrate a specific effect of inhalational anesthetics on mixing phase equilibria of a lipid mixture.     

Influence of ethanol on lipid membranes: from lateral pressure profiles to dynamics and partitioning

We have combined experiments with atomic-scale molecular dynamics simulations to consider the influence of ethanol on a variety of lipid membrane properties. We first employed isothermal titration calorimetry together with the solvent-null method to study the partitioning of ethanol molecules into saturated and unsaturated membrane systems. The results show that ethanol partitioning is considerably more favorable in unsaturated bilayers, which are characterized by their more disordered nature compared to their saturated counterparts. Simulation studies at varying ethanol concentrations propose that the partitioning of ethanol depends on its concentration, implying that the partitioning is a nonideal process. To gain further insight into the permeation of alcohols and their influence on lipid dynamics, we also employed molecular dynamics simulations to quantify kinetic events associated with the permeation of alcohols across a membrane, and to characterize the rotational and lateral diffusion of lipids and alcohols in these systems. The simulation results are in agreement with available experimental data and further show that alcohols have a small but non-vanishing effect on the dynamics of lipids in a membrane. The influence of ethanol on the lateral pressure profile of a lipid bilayer is found to be prominent: ethanol reduces the tension at the membrane-water interface and reduces the peaks in the lateral pressure profile close to the membrane-water interface. The changes in the lateral pressure profile are several hundred atmospheres. This supports the hypothesis that anesthetics may act by changing the lateral pressure profile exerted on proteins embedded in membranes.     

Distribution of local anesthetics in lipid membranes

Absorption of local anesthetics into lipid membranes and adsorption onto their surfaces were studied as a function of the pH of aqueous bulk solutions by measuring lipid vesicle electrophoretic mobility, the partition of the anesthetics between the aqueous and membrane phases by the use of fluorescence and radioactive tracer methods, and the effect of the anesthetics on interfacial tension of lipid monolayers formed at the oil/aqueous interface.

At a pH much lower than the pK_a value of the local anesthetic, the charged form of the local anesthetic was only adsorbed onto the membrane surface, as determined from vesicle electrophoretic mobility, radioisotope tracer and the monolayer surface tension studies. Surface partition coefficients of the charged form of the local anesthetics on phosphatidylcholine and phosphatidylserine membranes were obtained from the data of electrophoretic mobilities for lipid vesicles. The surface partition coefficients of various local anesthetics paralleled those of the bulk partition coefficients.

As the pH of the solutions increased, the adsorbed amount of the charged form of the anesthetic at the membrane interface decreased, while the absorption of the uncharged form of the local anesthetic into the membrane increased. The total amount of local anesthetic adsorbed per unit area of the membrane generally increased as the pH of the solution increased. This was also observed from the measurements of the fluorescence of local anesthetics adsorbed into the membranes. At lower pH than that corresponding to the pK_a value of the local anesthetic, the amount of anesthetic adsorbed depended greatly upon the membrane surface charge. At a higher pH than its pK_a, it did not depend appreciably on the surface charge density of the membrane but did depend on the bulk partition coefficients between the aqueous and oil phases.     

Binding of the general anesthetics chloroform and 2,2,2-trichloroethanol to the hydrophobic core of a four-alpha-helix bundle proteins

The structural features of general anesthetic binding sites on proteins are being examined using a defined model system consisting of a four-alpha-helix bundle scaffold with a hydrophobic core. Previous work suggested that halothane binding to the four-alpha-helix bundle was improved by (1) introducing a cavity into the hydrophobic core and (2) substituting a methionine side-chain in place of an alpha-helical heptad e position leucine. In this study, the ability of the general anesthetics chloroform and 2,2,2-trichloroethanol to bind to the hydrophobic core of the four-alpha-helix bundle (Aalpha2-L38M)2 is explored. The halogenated alkane chloroform binds with a dissociation constant (Kd) = 1.4 +/- 0.2 mM, whereas 2,2,2-trichloroethanol binds with a Kd = 19.5 +/- 1.2 mM. The affinity of both general anesthetics for the hydrophobic core of the four-alpha-helix bundle approximates their whole animal effective concentration in 50% of test subjects' (EC50) values, as shown previously for halothane. Tryptophan phosphorescence decay rates at 77 K are accelerated by a factor of 4.5 by both bound halothane and chloroform, indicating that the heavy-atom effect is responsible for a portion of the observed fluorescence quenching. Because heavy-atom effects are operative only at short distances, the findings indicate that these general anesthetics are binding in the vicinity of the indole rings of W15 in the hydrophobic core of the four-alpha-helix bundle scaffold. The results indicate that chloroform, halothane and 2,2,2-trichloroethanol may occupy the same sites on protein targets.     

Effects of isoflurane, halothane, and chloroform on the interactions and lateral organization of lipids in the liquid-ordered phase

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.     

Competitive adsorption of Ca and local anesthetics on lipid membrane surfaces

Adsorption fo tertriary amine local anesthetics and Ca²⁺ onto lipid membranes having various negative surface charge densities was studied by measuring lipid vesicle electrophoretic mobility.

As the surface charge density of the membrane was reduced, the adsorption of the local anesthetics dominated that of the divalent cation. For a relatively high negatively charged membrane, the adsorption of both local anesthetic and Ca²⁺ became comparable and competitive.

It is deduced that the major factor for the adsorption of local anesthetic onto lipid membranes is due to simple physical partitioning between aqueous and membrane phases, and not due to ionic type of binding as seen for divalent cations with membranes. However, the adsorption of anesthetics is influenced by the surface potential of membranes which is in turn related to the surface concentration of local anesthetics near the membrane.

The amounts of competitive adsorption of divalent cations and local anesthetics are analyzed with respect to their bulk concentrations and various surface charge densities of the membranes. With the results of the above studies, a possible interpretation for the interaction site as well as the mode of adsorption of local anesthetics onto axon membranes is made in relation to divalent cation concentrations in the bulk phases.     

Water solvent and local anesthetics: A computational study

There are various experimental studies regarding the toxicity and the time of action of local anesthetics, which contain general insights about their pharmacological and physicochemical properties. Although a detailed microscopic analysis of the local anesthetics would contribute to understanding these properties, there are relatively few theoretical studies about these molecules. In this article, we present the results from calculations performed for three local anesthetics: tetracaine, procaine, and lidocaine, both in their charged and uncharged forms, in aqueous environment. We have used the density functional theory and molecular dynamics simulations to study the structural characteristics of these compounds. The radial distribution function g(r) was used to examine the structure of water molecules surrounding different regions of the local anesthetics. We demonstrated the nonhomogeneous character of the anesthetics with respect to their affinity to water solvent molecules as well as the modifications in their affinity to water caused by changes in their charge state. We also observed that the biological potency of the anesthetics is more related to the behavior of specific groups within the molecule, which are responsible for the interaction with the lipid phase of membranes, rather than the general properties of the molecule as a whole. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Steered molecular dynamics simulations of Na+ permeation across the gramicidin A channel

The potential of mean forces (PMF) governing Na+ permeation through gramicidin A (gA) channels with explicit water and membrane was characterized using steered molecular dynamics (SMD) simulations. Constant-force SMD with a steering force parallel to the channel axis revealed at least seven energy wells in each monomer of the channel dimer. Except at the channel dimer interface, each energy well is associated with at least three and at most four backbone carbonyl oxygens and two water oxygens in a pseudo-hexahedral or pseudo-octahedral coordination with the Na+ ion. Repeated constant-velocity SMD by dragging a Na+ ion from each energy well in opposite directions parallel to the channel axis allowed the computation of the PMF across the gA channel, revealing a global minimum corresponding to Na+ binding sites near the entrance of gA at +/-9.3 A from the geometric center of the channel. The effect of volatile anesthetics on the PMF was also analyzed in the presence of halothane molecules. Although the accuracy of the current PMF calculation from SMD simulations is not yet sufficient to quantify the PMF difference with and without anesthetics, the comparison of the overall PMF profiles nevertheless confirms that the anesthetics cause insignificant changes to the structural makeup of the free energy wells along the channel and the overall permeation barrier. On average, the PMF appears less rugged in the outer part of the channel in the presence of anesthetics, consistent with our earlier finding that halothane interaction with anchoring residues makes the gA channel more dynamic. A causal relationship was observed between the reorientation of the coordinating backbone carbonyl oxygen and Na+ transit from one energy well to another, suggesting the possibility that even minute changes in the conformation of pore-lining residues due to dynamic motion could be sufficient to trigger the ion permeation. Because some of the carbonyl oxygens contribute to Na+ coordination in two adjacent energy wells, our SMD results reveal that the atomic picture of ion "hopping" through a gA channel actually involves a Na+ ion being carried in a relay by the coordinating oxygens from one energy well to the next. Steered molecular dynamics complements other computational approaches as an attractive means for the atomistic interpretation of experimental permeation studies.     

Morphological transitions in model membrane systems by the addition of anesthetics

The mechanism of anesthetic action on membranes is still an open question, regardless of their extensive use in medical practice. It has been proposed that anesthetics may have the effect of promoting pore formation across membranes or at least switching transmembrane channels. In both cases this may be the result of changes in the interfacial curvature of the membrane due to the presence of anesthetic molecules. Aqueous solutions of surfactants display phases that mimic, in a simplified manner, real biological membranes. Therefore, in this study, two nonionic surfactant systems C16E6/H2O in concentrated solution and C10E3/H2O in dilute solution have been used as model membranes for the investigation of the effects of six common anesthetics (halothane, sodium thiopental, lidocaine base form and hydrochloride, prilocaine hydrochloride, and ketamine hydrochloride). Both binary surfactant-water systems exhibit phase transitions from the lamellar phase, Lalpha, that has zero spontaneous curvature and zero monolayer curvature to phases with more local interfacial curvature. These are the random mesh phase, Mh1(0), which consists of lamellae pierced by water-filled pores with local areas of positive interfacial curvature and the sponge phase, L3, that consists of the lamellar phase with interlamellae attachments, often referred to as a "melted" cubic phase, possessing negative monolayer curvature. Small-angle X-ray scattering and 2H NMR experiments upon the C16E6/2H2O system and optical observations of the C10E3/H2O system showed that all anesthetics employed in this study cause a shift in the Mh1(0) to Lalpha phase transition temperature and in the Lalpha to L3 transition temperature, respectively. All of the anesthetics studied bind to the interfacial region of the surfactant systems. Two types of behavior were observed on anesthetic addition: type I anesthetics, which decreased interfacial curvature, and type II, which increased it. However, at physiological pH both types of anesthetics decreased interfacial curvature.     

Ab Initio Calculation of structures and properties of halogenated general anesthetics: halothane and sevoflurane

To correctly analyze the effects of general anesthetics on their potential targets by large‐scale molecular simulation, the structural parameters and partial atomic charges of the anesthetics are of determinant importance. Geometric optimizations using the Hartree–Fock and the B3LYP density functional theory methods with the large 6‐311+G(2d,p) basis set were performed to determine the structures and charge distributions of two halogenated anesthetics, 2‐bromo‐2‐chloro‐1,1,1‐trifluoroethane (halothane) and fluoromethyl‐2,2,2,‐trifluoro‐1‐(trifluoromethyl) ethyl ether (sevoflurane). The calculated bond lengths and angles are within 3% of the corresponding experimental values reported for the similar molecular groups. Charges are assigned using the Mulliken population analysis and the electrostatic potential (ESP) based on the Merz–Kollman–Singh scheme. The atoms‐in‐molecules (AIM) theory is also used to assign the charges in halothane. The dipole moments calculated with the Mulliken population analysis and ESP for the structures optimized by B3LYP/6‐311+(2d,p) were respectively 1.355 and 1.430 D for halothane and 2.255 and 2.315 D for sevoflurane. These are in excellent agreement with the experimental values of 1.41 and 2.33 D for halothane and sevoflurane, respectively. The calculated structures and partial charge distributions can be readily parameterized for molecular mechanics and molecular dynamics simulations involving these halogenated agents. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 436–444, 2001     

Effects of lidocaine on the expansion of lipid monolayer at air/water interface in relation to the local anesthesia

Lidocaine compounds have widely been used as local anesthetics. Regarding the molecular mechanism for anesthesia by local anesthetics, two hypotheses have been proposed. The first one is that molecules of local anesthetics penetrate into the hydrophobic region of cell membrane and expand the membrane volume, resulting in a change of protein conformation that blocks sodium permeability. The second hypothesis is that molecules of local anesthetics are directly adsorbed into the receptors of anesthetics in the protein channel without expanding the cell membrane. However, these proposals have never been examined systematically. In this study, the expansion of cell membrane by lidocaine compounds was investigated by employing lipid monolayer at the air/water interface as the mimetic system for cell membrane. It was found that oil-soluble lidocaine contracted the area/molecule of lipids in the monolayer of phosphatidyl choline, sphingomyelin, DS-PL95E and lipoid, but expand the monolayer of phosphatidyl ethanolamine only in a certain range of mixing ratios. Thus, this study can provide an evidence that lidocaine yields anesthesia effect by adsorbing into receptors in the protein channel rather than expanding the cell membrane.     

Diffusion of lipid-like single-molecule fluorophores in the cell membrane

The dicyanomethylenedihydrofuran (DCDHF) class of single-molecule fluorophores contains an amine donor and a dicyanomethylenedihydrofuran acceptor linked by a conjugated unit (benzene, naphthalene, or styrene). Molecules in this class have a number of useful properties in addition to those usually required for single-molecule studies (such as high fluorescence quantum yield and photostability), including second-order optical nonlinearity, large ground-state dipole moment, and sensitivity to local environment. Moreover, most DCDHF molecules have amphiphilic structures, with a polar dicyanomethylenedihydrofuran headgroup and nonpolar hydrocarbon tails on the amine or furan ring, and can be used as fluorescent lipid analogues for live cell imaging. Here we demonstrate that individual molecules of several different DCDHF lipid analogues can be observed diffusing in the plasma membrane of Chinese hamster ovary cells. The photophysical and diffusive behaviors of the DCDHF lipid analogues in membranes are described and are found to be competitive with the well-known lipid probe N-(6-tetramethylrhodaminethiocarbamoyl)-1,2-dihexadecanoyl-sn-glycero-3-phosphoethanolamine.     

Interactions of anesthetics with the membrane-water interface

Although the potency of conventional anesthetics correlates with lipophilicity, an affinity to water also is essential. It was recently found that compounds with very low affinities to water do not produce anesthesia regardless of their lipophilicity. This finding implies that clinical anesthesia might arise because of interactions at molecular sites near the interface of neuronal membranes with the aqueous environment and, therefore, might require increased concentrations of anesthetic molecules at membrane interfaces. As an initial test of this hypothesis, we calculated in molecular dynamics simulations the free energy profiles for the transfer of anesthetic 1,1,2-trifluoroethane and nonanesthetic perfluoroethane across water-membrane and water-hexane interfaces. Consistent with the hypothesis, it was found that trifluoroethane, but not perfluoroethane, exhibits a free energy minimum and, therefore, increased concentrations at both interfaces. The transfer of trifluoroethane from water to the nonpolar hexane or interior of the membrane is accompanied by a considerable, solvent-induced shift in the conformational equilibrium around the C-C bond.     

α-生育酚在模型生物膜中的分子动力学模拟

用分子动力学方法模拟了280, 310和350 K下α-生育酚在二豆蔻酰磷脂酰胆碱、二豆蔻酰磷脂酰乙醇胺、二硬脂酰磷脂酰胆碱和二硬脂酰磷脂酰乙醇胺双层膜中的性质, 包括了空间位置、氢键、取向和动力学性质, 取得了如下的结论. 第一, 生育酚头部的羟基一般位于脂双层亲疏水界面的下方, 升高温度将促进羟基向膜双层的中心移动, 在350 K时观察到了在上下两个单层间的翻转. 第二, 生育酚主要与磷脂的酯基形成氢键, 几乎不与磷脂酰乙醇胺的氨基形成氢键; 比较生育酚与磷脂酰胆碱和乙醇胺形成的氢键后发现, 后者更稳定. 第三, 生育酚的头部在膜中取向多变, 与膜的法线夹角不固定, 尾部的构象也很复杂. 第四, 在温度较低时, 生育酚的侧向扩散系数与磷脂的相当, 但在350 K时其扩散速度明显加快; 在垂直方向生育酚的扩散速度很慢.     

Increase in water permeability of negatively charged liposomal membrane by local anesthetics.

Effect of the local anesthetics dibucaine, tetracaine, lidocaine and procaine on the water permeability of phospholipid membrane was examined using liposomes composed of bovine heart cardiolipin and egg yolk phosphatidylcholine in a molar ratio of 2/98 by monitoring the osmotic shrinkage of liposomes in hypertonic glucose solution at pH 7.3 and 30 degrees C. These local anesthetics greatly accelerated the water permeability by destabilizing the membrane structure. The effect was found to be governed by the hydrophobicity of the anesthetics. There was also a significant correlation between the membrane destabilizing actions and the anesthetic activities.