Effects of volatile anesthetic on channel structure of gramicidin A

Orthosteric- versus allosteric-dependent activation of the GABAA receptor requires numerically distinct subunit level rearrangements

Anaesthetic molecules act on synaptic transmission via the allosteric modulation of ligand-gated chloride channels, such as hetero-oligomeric α₁β₂γ₂ GABA_A receptors. To elucidate the overall activation paradigm via allosteric versus orthosteric sites, we used highly homologous, but homo-oligomeric, ρ₁ receptors that are contrastingly insensitive to anaesthetics and respond partially to several full GABA α₁β₂γ₂ receptor agonists. Here, we coexpressed varying ratios of RNAs encoding the wild-type and the mutated ρ₁ subunits, which are anaesthetic-sensitive and respond with full efficacy to partial GABA agonists, to generate distinct ensembles of receptors containing five, four, three, two, one, or zero mutated subunits. Using these experiments, we then demonstrate that, in the pentamer, three anaesthetic-sensitive ρ₁ subunits are needed to impart full efficacy to the partial GABA agonists. By contrast, five anaesthetic-sensitive subunits are required for direct activation by anaesthetics alone, and only one anaesthetic-sensitive subunit is sufficient to confer the anaesthetic-dependent potentiation to the GABA current. In conclusion, our data indicate that GABA and anaesthetics holistically activate the GABA_A ρ₁ receptor through distinct subunit level rearrangements and suggest that in contrast to the global impact of GABA via orthosteric sites, the force of anaesthetics through allosteric sites may not propagate to the neighbouring subunits and, thus, may have only a local and limited effect on the ρ₁ GABA_A receptor model system.

Physical and physicochemical factors effecting transport of chlorohydrocarbon gases from lung alveolar air to blood as measured by the causation of narcosis

This systematic investigation examines gas transport in the lung for two sets of chlorohydrocarbons (CHCs): the chloromethanes (C1) and chloroethanes (C2). The C1 series includes chloromethane, methylene chloride, chloroform, and carbon tetrachloride, and the C2 series includes chloroethane, 1,2-dichloroethane, 1, 1, 2-trichloroethane, and 1, 1, 2, 2-tetrachloroethane. Most CHC gases cause narcosis. The comprehensive narcosis work of Lehmann and colleagues on CHCs was used as a basis for the narcosis endpoint in the present examination. The sites for narcosis are located in the brain (midline cortex and posterior parietal area), the spine, and at many peripheral nerve sites. Central nervous system (CNS) exposure executes a multisite, neural transmission set of inhibitions that promotes rapid loss of consciousness, sensory feeling, and current and stored memory while providing temporary amnesia. Absorption into the system requires dissolution into many lipid membranes and binding to lipoproteins. Lipophilicity is a CHC property shared with many anesthetics according to the Meyer-Overton Rule. Many structurally different lipid chemicals produce the narcosis response when the lipid concentration exceeds -67 mM. This suggests narcotic or anesthetic dissolution into CNS membranes until the lipid organization is disrupted or perturbed. This perturbation includes loading of Na(+)- and K(+)-channel transmembrane lipoprotein complexes and disrupting their respective channel functional organizations. The channel functions become attenuated or abrogated until the CHC exposure ceases and CHC loading reverses. This investigation demonstrates how the CHC physical and chemical properties influence the absorption of these CHCs via the lung and the alveolar system on route to the blood. Narcosis in test animals was used here as an objective biological endpoint to study the effects of the physical factors Bp, Vp, Kd (oil: gas) partition, Henry's constant (HK), and water solubility (S%) on gas transport. Narcosis is immediate after gas exposure and requires no chemical activation only absorption into the blood and circulation to CNS narcotic sites. The three physical factors Bp, K(d) (oil: air), and S% vary directly with unitary narcosis (UN) whereas Vp and HK vary inversely with UN in linear log-log relationships for the C2 series but not for the C1 series. Physicochemical properties of C1 series gases indicate why they depart from what is usually assumed to be an Ideal Gas. An essential discriminating process in the distal lung is the limiting alveolar film layer (AFL) and the membrane layer of the alveolar acini. The AFL step influences gas uptake by physically limiting the absorption process. Interaction with and dissolution into aqueous solvent of the AFL is required for transport and narcotic activity. Narcotics or anesthetics must engage the aqueous AFL with sufficient strength to allow transport and absorption for downstream CNS binding. CHCs that do not engage well with the AFL are not narcotic. Lipophilicity and amphipathicity are also essential solvency properties driving narcotics' transport through the alveolar layer, delivery to the blood fats and lipoproteins, and into critical CNS lipids, lipoproteins, and receptor sites that actuate narcosis. AFL disruption is thought to be strongly related to a number of serious pulmonary diseases such acute respiratory distress syndrome, infant respiratory distress syndrome, emphysema, chronic obstructive pulmonary disease, asthma, chronic bronchitis, pneumonia, pulmonary infections, and idiopathic pulmonary fibrosis. The physical factors (Bp, Vp, Kd [oil: gas] partition, Henry's constant, and water solubility [S%]) combine to affect a specific transport through the AFL if lung C > C(0) (threshold concentration for narcosis). The degree of blood CHC absorption depends on dose, lipophilicity, and lung residence time. AFL passage can be manipulated by physical factors of increased pressure (kPa) or increased gas exposure (moles). Molecular lipophilicity facilitates narcosis but lipophilicity alone does not explain narcosis. Vapor pressure is also required for narcosis. Narcotic activity apparently requires stereospecific processing in the AFL and/or down-stream inhibition at stereospecific lipoproteins at CNS inhibitory sites. It is proposed that CHCs likely cannot proceed through the AFL without perturbation or disruption of the integrity of the AFL at the alveoli. CHC physicochemical properties are not expected to allow their transport through the AFL as physiological CO(2) and O(2) naturally do in respiration. This work considers CHC inspiration and systemic absorption into the blood with special emphasis on the CHC potential perturbation effects on the lipid, protein liquid layer supra to the alveolar membrane (AFL). A heuristic gas transport model for the CHCs is presented as guidance for this examination. The gas transport model can be used to study absorption for other gas delivery endpoints of environmental concern such as carcinogens.

An atomistic model for simulations of the general anesthetic isoflurane

An atomistic model of isoflurane is constructed and calibrated to describe its conformational preferences and intermolecular interactions. The model, which is compatible with the CHARMM force field for biomolecules, is based on target quantities including bulk liquid properties, molecular conformations, and local interactions with isolated water molecules. Reference data is obtained from tabulated thermodynamic properties and high-resolution structural information from gas-phase electron diffraction, as well as DFT calculations at the B3LYP level. The model is tested against experimentally known solvation properties in water and oil, and shows quantitative agreement. In particular, isoflurane is faithfully described as lipophilic, yet nonhydrophobic, a combination of properties critical to its pharmacological activity. Intermolecular interactions of the model are further probed through simulations of the binding of isoflurane to a binding site in horse spleen apoferritin (HSAF). The observed binding mode compares well with crystallographic data, and the calculated binding affinities are compatible with experimental results, although both computational and experimental measurements are challenging and provide results with limited precision. The model is expected to be useful for detailed simulations of the elementary molecular processes associated with anesthesia. Full parameters are provided as Supporting Information.

Computational studies on the interactions of inhalational anesthetics with proteins

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

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.

Fluorescence quenching of gramicidin D in model membranes by halothane

Inhaled anesthetics were introduced in surgery over a century ago. To this day, the molecular mechanism of anesthetic action remains largely unknown. However, ion-channels of neuronal membranes are believed to be the most- likely molecular targets of inhaled anesthetics. In the study presented here, we investigated the interaction of a simplified ion-channel system, gramicidin, with halothane, a small haloalkane inhaled anesthetic in various environments. Fluorescence-quenching experiments of gramicidin D in dioleoylphosphatidylcholine (DOPC) large unilamellar vesicles (LUVS) have shown that halothane can directly interact with the ion channel (KSV = 66 M–1). Halothane quenched the fluorescence from tryptophan residues located at the lipid bilayer – aqueous interfaces as well as those tryptophans located deeper in the bilayer. Quenching data from gramicidin D in sodium dodecyl sulfide (SDS) micelles revealed that the tryptophan residues located at the micelle–solvent interface were preferentially quenched by halothane (KSV = 22 M–1). In 1-octanol, fluorescence quenching was observed, but with a lower KSV value (KSV = 6 M–1) than in DOPC LUVS and SDS micelles. Taken together, these results indicate that halothane interactions with gramicidin, mediated by a lipid bilayer, are the strongest, and that the mechanism of anesthetic action may also be lipid-mediated.

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.

Interactions of anesthetics with their targets: non-specific, specific or both?

What makes a general anesthetic a general anesthetic? We shall review first what general anesthesia is all about and which drugs are being used as anesthetics. There is neither a unique definition of general anesthesia nor any consensus on how to measure it. Diverse drugs and combinations of drugs generate general anesthetic states of sometimes very different clinical quality. Yet the principal drugs are still considered to belong to the same class of 'general anesthetics'. Effective concentrations of inhalation anesthetics are in the high micromolar range and above, and even for intravenous anesthetics they do not go below the micromolar range. At these concentrations, many molecular and higher level targets are affected by inhalation anesthetics, fewer probably by intravenous anesthetics. The only physicochemical characteristic shared by anesthetics is the correlation of their anesthetic potencies with hydrophobicity. These correlations depend on the group of general anesthetics considered. In this review, anesthetic potencies for many different targets are plotted against octanol/water partition coefficients as measure of hydrophobicity. Qualitatively, similar correlations result, suggesting several but weak interactions with proteins as being characteristic of anesthetic actions. The polar interactions involved are weak, being roughly equal in magnitude to hydrophobic interactions. Generally, intravenous anesthetics are noticeably more potent than inhalation anesthetics. They differ considerably more between each other in their interactions with various targets than inhalation anesthetics do, making it difficult to come to a decision which of these should be used in future studies as representative 'prototypical general anesthetics'.

The orientation effect of gramicidin A on bicelles and Eu3+-doped bicelles as studied by solid-state NMR and FT-IR spectroscopy

We have explored the effect of gramicidin A (gA) on bicelle (Bic) orientation in the absence and presence of Eu(3+) by (31)P and (2)H NMR at different DMPC/gA ratios. FT-IR spectroscopy was used to assess the lipid chain ordering and verify the transmembrane peptide conformation. Our results show a time-dependent flipping of the bilayer normal alignment at high temperatures and high proportion of gA. The results are explained by both the diamagnetic susceptibility anisotropy of the beta(6.3) helical peptides and viscosity of the lipid mixture. The concentration effect of gramicidin on Bic/Eu(3+) is compared to that on Eu(3+)-doped DMPC liposomes. The Bic/Eu(3+) system is no longer oriented in the presence of gA and adopts a vesicular morphology while the peptide incorporation induces the formation of ellipsoidal DMPC/Eu(3+) assemblies aligned with their normal parallel to the magnetic field. The difference is explained in terms of lipid chain disorder and size of the bilayers.

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.

Molecular dynamics simulations of C2F6 effects on gramicidin A: implications of the mechanisms of general anesthesia

It was recently postulated that the effects of general anesthetics on protein global dynamics might underlie a unitary molecular mechanism of general anesthesia. To verify that the specific dynamics effects caused by general anesthetics were not shared by nonanesthetic molecules, two parallel 8-ns all-atom molecular dynamics simulations were performed on a gramicidin A (gA) channel in a fully hydrated dimyristoylphosphatidylcholine membrane in the presence and absence of hexafluoroethane (HFE), which structurally resembles the potent anesthetic molecule halothane but produces no anesthesia. Similar to halothane, HFE had no measurable effects on the gA channel structure. In contrast to halothane, HFE produced no significant changes in the gA channel dynamics. The difference between halothane and HFE on channel dynamics can be attributed to their distinctly different distributions within the lipid bilayer and consequently to the different interactions of the anesthetic and the nonanesthetic molecules with the channel-anchoring tryptophan residues. The study further supports the notion that anesthetic-induced changes in protein global dynamics may play an important role in mediating anesthetic actions on proteins.

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

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