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

Absolute Differential Cross-Sections for Elastic Electron Scattering from Sevoflurane Molecule in the Energy Range from 50-300 eV

We report the results of the measurements and calculations of the absolute differential elastic electron scattering cross-sections (DCSs) from sevoflurane molecule (C4H3F7O). The experimental absolute DCSs for elastic electron scattering were obtained for the incident electron energies from 50 eV to 300 eV, and for scattering angles from 25° to 125° using a crossed electron/target beams setup and the relative flow technique for calibration to the absolute scale. For the calculations, we have used the IAM-SCAR+I method (independent atom model (IAM) applying the screened additivity rule (SCAR) with interference terms included (I)). The molecular cross-sections were obtained from the atomic data by using the SCAR procedure, incorporating interference term corrections, by summing all the relevant atomic amplitudes, including the phase coefficients. In this approach, we obtain the molecular differential scattering cross-section (DCS), which, integrated over the scattered electron angular range, gives the integral scattering cross-section (ICS). Calculated cross-sections agree very well with experimental results, in the whole energy and angular range.

Modern Anesthetic Ethers Demonstrate Quantum Interactions with Entangled Photons

Despite decades of research, the mechanism of anesthetic-induced unconsciousness remains incompletely understood, with some advocating for a quantum mechanical basis. Despite associations between general anesthesia and changes in physical properties such as electron spin, there has been no empirical demonstration that general anesthetics are capable of functional quantum interactions. In this work, we studied the linear and non-linear optical properties of the halogenated ethers sevoflurane (SEVO) and isoflurane (ISO), using UV-Vis spectroscopy, time dependent-density functional theory (TD-DFT) calculations, classical two-photon spectroscopy, and entangled two-photon spectroscopy. We show that both of these halogenated ethers interact with pairs of 800 nm entangled photons while neither interact with 800 nm classical photons. By contrast, nonhalogenated diethyl ether does not interact with entangled photons. This is the first experimental evidence that halogenated anesthetics can directly undergo quantum interaction mechanisms, offering a new approach to understanding their physicochemical properties.

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

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

Mapping the sevoflurane-binding sites of calmodulin

General anesthetics, with sevoflurane (SF) being the first choice inhalational anesthetic agent, provide reversible, broad depressor effects on the nervous system yet have a narrow margin of safety. As characterization of low-affinity binding interactions of volatile substances is exceptionally challenging with the existing methods, none of the numerous cellular targets proposed as chief protagonists in anesthesia could yet be confirmed. The recognition that most critical functions modulated by volatile anesthetics are under the control of intracellular Ca²⁺ concentration, which in turn is primarily regulated by calmodulin (CaM), motivated us for characterization of the SF–CaM interaction. Solution NMR (Nuclear Magnetic Resonance) spectroscopy was used to identify SF-binding sites using chemical shift displacement, NOESY and heteronuclear Overhauser enhancement spectroscopy (HOESY) experiments. Binding affinities were measured using ITC (isothermal titration calorimetry). SF binds to both lobes of (Ca²⁺)₄-CaM with low mmol/L affinity whereas no interaction was observed in the absence of Ca²⁺. SF does not affect the calcium binding of CaM. The structurally closely related SF and isoflurane are shown to bind to the same clefts. The SF-binding clefts overlap with the binding sites of physiologically relevant ion channels and bioactive small molecules, but the binding affinity suggests it could only interfere with very weak CaM targets.

Theoretical investigation of atmospheric chemistry of volatile anaesthetic sevoflurane: reactions with the OH radicals and atmospheric fate of the alkoxy radical (CF3)2CHOCHFO: thermal decomposition vs. oxidation

Reaction profile (kcal mol⁻¹) for (CF₃)₂CHOCHFO radical at the M06-2X/6-311++G(d,p) level.

George M, Freedman H, Barakat KH, Damaraju S, Hameroff S, Tuszynski JA. Computational predictions of volatile anesthetic interactions with the microtubule cytoskeleton: implications for side effects of general anesthesia

The cytoskeleton is essential to cell morphology, cargo trafficking, and cell division. As the neuronal cytoskeleton is extremely complex, it is no wonder that a startling number of neurodegenerative disorders (including but not limited to Alzheimer's disease, Parkinson's disease and Huntington's disease) share the common feature of a dysfunctional neuronal cytoskeleton. Recently, concern has been raised about a possible link between anesthesia, post-operative cognitive dysfunction, and the exacerbation of neurodegenerative disorders. Experimental investigations suggest that anesthetics bind to and affect cytoskeletal microtubules, and that anesthesia-related cognitive dysfunction involves microtubule instability, hyper-phosphorylation of the microtubule-associated protein tau, and tau separation from microtubules. However, exact mechanisms are yet to be identified. In this paper the interaction of anesthetics with the microtubule subunit protein tubulin is investigated using computer-modeling methods. Homology modeling, molecular dynamics simulations and surface geometry techniques were used to determine putative binding sites for volatile anesthetics on tubulin. This was followed by free energy based docking calculations for halothane (2-bromo-2-chloro-1,1,1-trifluoroethane) on the tubulin body, and C-terminal regions for specific tubulin isotypes. Locations of the putative binding sites, halothane binding energies and the relation to cytoskeleton function are reported in this paper.

The conformational landscape of the volatile anesthetic sevoflurane

Most of the volatile haloorganic compounds used as anesthetics exhibit a heavy-atom frame large enough to allow for conformational changes. Even in the absence of directed intermolecular interactions, only some or just one of the possible conformations might have an appreciable abundance. In this realm, the structure of the anesthetic haloether sevoflurane (CH(2)F-O-CH(CF(3))(2)) has been resolved using Fourier-transform microwave (FT-MW) spectroscopy in a supersonic-jet expansion. In isolated conditions sevoflurane adopts a single conformation characterized by a gauche fluoromethoxy group and a near-symmetric orientation of the isopropyl group with respect to the ether plane (cis H-C(ipr)-O-C(F)). Substitution and effective structures have been calculated from the rotational spectra of all (13)C and (18)O monosubstituted isotopic species observed in natural abundance. The electric dipole moment components were determined from additional Stark effect measurements. The experimental structures and rotational data are compared with those obtained from supporting ab initio predictions using MP2 calculations and the B3LYP hybrid functional.

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

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

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

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

Volatile anesthetic binding to proteins is influenced by solvent and aliphatic residues

The main objective of this work was to characterize VA binding sites in multiple anesthetic target proteins. A computational algorithm was used to quantify the solvent exclusion and aliphatic character of amphiphilic pockets in the structures of VA binding proteins. VA binding sites in the protein structures were defined as the pockets with solvent exclusion and aliphatic character that exceeded minimum values observed in the VA binding sites of serum albumin, firefly luciferase, and apoferritin. We found that the structures of VA binding proteins are enriched in these pockets and that the predicted binding sites were consistent with experimental determined binding locations in several proteins. Autodock3 was used to dock the simulated molecules of 1,1,1,2,2-pentafluoroethane, difluoromethyl 1,1,1,2-tetrafluoroethyl ether, and sevoflurane and the isomers of halothane and isoflurane into these potential binding sites. We found that the binding of the various VA molecules to the amphiphilic pockets is driven primarily by VDW interactions and to a lesser extent by weak hydrogen bonding and electrostatic interactions. In addition, the trend in Delta G binding values follows the Meyer-Overton rule. These results suggest that VA potencies are related to the VDW interactions between the VA ligand and protein target. It is likely that VA bind to sites with a high degree of solvent exclusion and aliphatic character because aliphatic residues provide favorable VDW contacts and weak hydrogen bond donors. Water molecules occupying these sites maintain pocket integrity, associate with the VA ligand, and diminish the unfavorable solvation enthalpy of the VA. Water molecules displaced into the bulk by the VA ligand may provide an additional favorable enthalpic contribution to VA binding. Anesthesia is a component of many health related procedures, the outcomes of which could be improved with a better understanding of the molecular targets and mechanisms of anesthetic action.

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

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

Dynamics of firefly luciferase inhibition by general anesthetics: Gaussian and anisotropic network analyses

The new crystal structures of the product-bound firefly luciferase combined with the previously determined substrate-free structures allow for a detailed analysis of the dynamics basis for the luciferase enzymatic activities. Using the Gaussian network model and the anisotropic network model, we show here that the superposition of the three slowest anisotropic network model modes, consisting of the bending, rotating, and rocking motions of the C-domain, accounts for large rearrangement of domains from the substrate-free (open) to product-bound (closed) conformation and thus constitutes a critical component of the enzyme's functions. The analysis also offers a unique platform to reexamine the molecular mechanism of the anesthetic inhibition of the firefly luciferase. Through perturbing the protein backbone network by introducing additional nodes to represent anesthetics, we found that the presence of two representative anesthetics, halothane and n-decanol, in different regions of luciferase had distinctively different effects on the protein's global motion. Only at the interface of the C- and N-domains did the anesthetics cause the most profound reduction in the overall flexibility of the C-domain and the concomitant increase in the flexibility of the loop, where the substitution of a conserved lysine residue was found experimentally to lead to >2-3 orders of magnitude reduction in activity. These anesthetic-induced dynamics changes can alter the normal function of the protein, appearing as an epiphenomenon of an "inhibition". The implication of the study is that a leading element for general anesthetic action on proteins is to disrupt the modes of motion essential to protein functions.

Molecular surface electrostatic potentials and anesthetic activity

General anesthetics apparently act through weak, noncovalent and reversible interactions with certain sites in appropriate brain proteins. As a means of gaining insight into the factors underlying anesthetic potency, we have analyzed the computed electrostatic potentials V (S)(r) on the surfaces of 20 molecules with activities that vary between zero and high. Our results are fully consistent with, and help to interpret, what has been observed experimentally. We find that an intermediate level of internal charge separation is required; this is measured by Pi, the average absolute deviation of V (S)(r), and the approximate window is 7 < Pi < 13 kcal mol(-1). This fits in well with the fact that anesthetics need to be lipid soluble, but also to have some degree of hydrophilicity. We further show that polyhalogenated alkanes and ethers, which include the most powerful known anesthetics, have strong positive potentials, V (S,max), associated with their hydrogens, chlorines and bromines (but not fluorines). These positive sites may impede the functioning of key brain proteins, for example by disrupting their normal hydrogen-bond patterns. It has indeed been recognized for some time that the most active polyhalogenated alkanes and ethers contain hydrogens usually in combination with chlorines and/or bromines.

Anesthetic interaction with ketosteroid isomerase: insights from molecular dynamics simulations

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.

Structural basis for high‐affinity volatile anesthetic binding in a natural 4‐helix bundle protein

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.

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

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

The role of electrostatic interactions in human serum albumin binding and stabilization by halothane

Electrostatic interactions have been proposed as a potentially important force for anesthetics and protein binding but have not yet been tested directly. In the present study, we used wild-type human serum albumin (HSA) and specific site-directed mutants as a native protein model to investigate the role of electrostatic interactions in halothane binding. Structural geometry analysis of the HSA-halothane complex predicted an absence of significant electrostatic interactions, and direct binding (tryptophan fluorescence and zonal elution chromatography) and stability experiments (hydrogen exchange) confirmed that loss of charge in the binding sites, by charged to uncharged mutations and by changing ionic strength of the buffer, generally increased both regional (tryptophan region) and global halothane/HSA affinity. The results indicate that electrostatic interactions (full charges) either do not contribute or diminish halothane binding to HSA, leaving only the more general hydrophobic and van der Waals forces as the major contributors to the binding interaction.

Halogenated diazirines as photolabel mimics of the inhaled haloalkane anesthetics

The inhaled anesthetics are low affinity volatile compounds whose mechanism of action remains unclear, in part due to the difficulty of determining their binding targets. Photolabeling may help resolve this difficulty, and thus we have synthesized six compounds (four previously unreported) with structural and physical similarity to halothane (1-bromo-1-chloro-2,2,2-trifluoroethane), a commonly used clinical anesthetic. These compounds incorporate either a diazo, diazirine, or azido group to provide photolability in the long-UV range and to provide a highly reactive photolysis product. While several of the compounds have immobilizing activity in tadpoles, it is complicated by either toxicity or very low potency. One compound however, a halogenated three-carbon diazirine 4, is a potent anesthetic, is apparently nontoxic, potentiates GABA(A) Cl(-) currents, and stabilizes serum albumin, all of which are features of halothane. When tagged with radioactivity, this compound should serve as a reasonable probe of haloalkane anesthetic binding targets and sites.