Partitioning of anesthetics into a lipid bilayer and their interaction with membrane-bound peptide bundles

Dynamic framework for large-scale modeling of membranes and peripheral proteins

In this chapter, we present a novel computational framework to study the dynamic behavior of extensive membrane systems, potentially in interaction with peripheral proteins, as an alternative to conventional simulation methods. The framework effectively describes the complex dynamics in protein-membrane systems in a mesoscopic particle-based setup. Furthermore, leveraging the hydrodynamic coupling between the membrane and its surrounding solvent, the coarse-grained model grounds its dynamics in macroscopic kinetic properties such as viscosity and diffusion coefficients, marrying the advantages of continuum- and particle-based approaches. We introduce the theoretical background and the parameter-space optimization method in a step-by-step fashion, present the hydrodynamic coupling method in detail, and demonstrate the application of the model at each stage through illuminating examples. We believe this modeling framework to hold great potential for simulating membrane and protein systems at biological spatiotemporal scales, and offer substantial flexibility for further development and parametrization.

Role of Disulphide Bonds in Membrane Partitioning of a Viral Peptide

The importance of disulphide bond in mediating viral peptide entry into host cells is well known. In the present work, we elucidate the role of disulphide (SS) bond in partitioning mechanism of membrane-active Hepatitis A Virus-2B (HAV-2B) peptide, which harbours three cysteine residues promoting formation of multiple SS-bonded states. The inclusion of SS-bond not only results in a compact conformation but also induces distorted α-helical hairpin geometry in comparison to SS-free state. Owing to these, the hydrophobic residues get buried, restricting the insertion of SS-bonded HAV-2B peptide into lipid packing defects and thus the partitioning of the peptide is completely or partly abolished. In this way, the disulphide bond can potentially regulate the partitioning of HAV-2B peptide such that the membrane remodelling effects of this viral peptide are significantly reduced. The current findings may have potential implications in drug designing, targeting the HAV-2B protein by promoting disulphide bond formation within its membrane-active region.

General anaesthetics as 'awakening agents'? Re-appraising the evidence for suggested 'pressure reversal' of anaesthesia

Increasing ambient pressure has been suggested to reverse general anaesthesia and provides support for the 'lipid theory'. Anaesthetic dissolution into cell membranes is said to cause their expansion to a critical volume. This triggers a sequence of events as basis of a unitary theory of anaesthestic mechanism. Pressure is argued to restore membrane volume to below critical level, reversing this process. We wished to review the original literature to assess internal consistency within and across papers, and to consider if alternative interpretations were possible. A literature search yielded 31 relevant 'pressure reversal' papers for narrative review, and 8 papers that allowed us to re-plot original data more consistently as 'dose-response' curves for the anaesthetics examined. Original studies were heterogenous for end-points, pressure ranges, species, and agents. Pressure effects were inconsistent, with narcosis at certain pressures and excitation at others, influenced by carrier gas (e.g., nitrogen vs helium). Pressure reversal (a right- or downward-shift on the re-plotted dose-response curves) was evident, but only in some species and at certain pressures and anaesthetic concentrations. However, even more striking was a novel 'awakening' effect of anaesthetics: i.e., anaesthetics reversed the narcotic effect of pressure, but this was limited to certain pressures at generally low anaesthetic concentrations. Contrary to the established view, 'pressure reversal' is not a universal phenomenon. The awakening effect of anaesthetics - described here for the first time - has equal evidence to support it, within the same literature, and is something that cannot be fully explained. Pressure cannot meaningfully be used to gain insight into anaesthetic mechanisms because of its heterogenous, non-specific and unpredictable effects on biological systems.

Effect of cholesterol on the membrane partitioning dynamics of hepatitis A virus-2B peptide

Understanding viral peptide detection and partitioning and the subsequent host membrane composition-based response is essential for gaining insights into the viral mechanism. Here, we probe the crucial role of the presence of membrane lipid packing defects, depending on the membrane composition, in allowing the viral peptide belonging to C-terminal Hepatitis A Virus-2B (HAV-2B) to detect, attach and subsequently partition into host cell membrane mimics. Using molecular dynamics simulations, we conclusively show that the hydrophobic residues in the viral peptide detect transiently present lipid packing defects, insert themselves into such defects, form anchor points and facilitate the partitioning of the peptide, thereby inducing membrane disruption. We also show that the presence of cholesterol significantly alters such lipid packing defects, both in size and in number, thus mitigating the partitioning of the membrane active viral peptide into cholesterol-rich membranes. Our results are in excellent agreement with previously published experimental data and further explain the role of lipid defects in understanding such data. These results show differential ways in which the presence and absence of cholesterol can alter the permeability of the host membranes to the membrane active peptide component of HAV-2B virus, via lipid packing defects, and can possibly be a part of the general membrane detection mechanism for viroporins.

Investigation on penetration of saffron components through lipid bilayer bound to spike protein of SARS-CoV-2 using steered molecular dynamics simulation

A coronavirus identified as COVID-19 is the reason for an infection outbreak which is started in December 2019. NO completely effective drugs and treatments are not recognized for this virus. Recently, saffron and its compounds were used to treat different viral diseases. Saffron extract and its major ingredients have shown antiviral effects. In this study, the steered molecular dynamics simulation was used for investigating the effect of four main components of saffron that include: crocin, crocetin, safranal, and picrocrocin as candidate for drug molecules, on COVID-19. The binding energies between drug molecules and spike protein and the main protease of the virus were evaluated. The obtained results based on Lennard-Jones and electrostatic potentials demonstrated that crocetin has a high affinity towards spike protein and also the main protease of the virus. Also, the quantum mechanics calculations elucidated that the crocetin could overcome energy barrier of lipid bilayer with strong dipole moment and polarizability. The pharmacokinetic and ADMET properties proved that crocetin could be a suitable drug candidate. So, crocetin could be a promising drug for treatment of COVID-19.

Competitive Interactions between Halothane and Isoflurane at the Carotid Body and TASK Channels

Background

The degree to which different volatile anesthetics depress carotid body hypoxic response relates to their ability to activate TASK potassium channels. Most commonly, volatile anesthetic pairs act additively at their molecular targets. We examined whether this applied to carotid body TASK channels.

Methods

We studied halothane and isoflurane effects on hypoxia-evoked rise in intracellular calcium (Ca2+i, using the indicator Indo-1) in isolated neonatal rat glomus cells, and TASK single-channel activity (patch clamping) in native glomus cells and HEK293 cell line cells transiently expressing TASK-1.

Results

Halothane (5%) depressed glomus cell Ca2+i hypoxic response (mean ± SD, 94 ± 4% depression; P < 0.001 vs. control). Isoflurane (5%) had a less pronounced effect (53 ± 10% depression; P < 0.001 vs. halothane). A mix of 3% isoflurane/1.5% halothane depressed cell Ca2+i response (51 ± 17% depression) to a lesser degree than 1.5% halothane alone (79 ± 15%; P = 0.001), but similar to 3% isoflurane alone (44 ± 22%; P = 0.224), indicating subadditivity. Halothane and isoflurane increased glomus cell TASK-1/TASK-3 activity, but mixes had a lesser effect than that seen with halothane alone: 4% halothane/4% isoflurane yielded channel open probabilities 127 ± 55% above control, versus 226 ± 12% for 4% halothane alone (P = 0.009). Finally, in HEK293 cell line cells, progressively adding isoflurane (1.5 to 5%) to halothane (2.5%) reduced TASK-1 channel activity from 120 ± 38% above control, to 88 ± 48% (P = 0.034).

Conclusions

In all three experimental models, the effects of isoflurane and halothane combinations were quantitatively consistent with the modeling of weak and strong agonists competing at a common receptor on the TASK channel.

Editor’s Perspective

What We Already Know about This Topic

What This Article Tells Us That Is New

Menthol Binding to the Human α4β2 Nicotinic Acetylcholine Receptor Facilitated by Its Strong Partitioning in the Membrane

We utilize various computational methodologies to study menthol's interaction with multiple organic phases, a lipid bilayer, and the human α4β2 nicotinic acetylcholine receptor (nAChR), the most abundant nAChR in the brain. First, force field parameters developed for menthol are validated in alchemical free energy perturbation simulations to calculate solvation free energies of menthol in water, dodecane, and octanol and compare the results against experimental data. Next, umbrella sampling is used to construct the free energy profile of menthol permeation across a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer. The results from a flooding simulation designed to study the water-membrane partitioning of menthol in a POPC lipid bilayer are used to determine the penetration depth and the preferred orientation of menthol in the bilayer. Finally, employing both docking and flooding simulations, menthol is shown to bind to different sites on the human α4β2 nAChR. The most likely binding mode of menthol to a desensitized membrane-embedded α4β2 nAChR is identified to be via a membrane-mediated pathway in which menthol binds to the sites at the lipid-protein interface after partitioning in the membrane. A rare but distinct binding mode in which menthol binds to the extracellular opening of receptor's ion permeation pore is also reported.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Molecular dynamics simulation study of the effect of halothane on mixed DPPC/DPPE phospholipid membranes

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.

Location and Character of Volatile General Anesthetics Binding Sites in the Transmembrane Domain of TRPV1

It has been proposed that general anesthesia results from direct multisite interactions with multiple and diverse ion channels in the brain. An understanding of the mechanisms by which general anesthetics modulate ion channels is essential to clarify their underlying behavior and their role in reversible immobilization and amnesia. Despite the fact that volatile general anesthetics are drugs that primarily induce insensitivity to pain, they have been reported to sensitize and active the vanilloid-1 receptor, TRPV1, which is known to mediate the response of the nervous system to certain harmful stimuli and which plays a crucial role in the pain pathway. Currently, the mechanism of action of anesthetics is unknown and the precise molecular sites of interaction have not been identified. Here, using ∼2.5 μs of classical molecular dynamics simulations and metadynamics, we explore these enigmas. Binding sites are identified and the strength of the association is further characterized using alchemical free-energy calculations. Anesthetic binding/unbinding proceeds primarily through a membrane-embedded pathway, and subsequently, a complex scenario is established involving multiple binding sites featuring single or multiple occupancy states of two small volatile drugs. One of the five anesthetic binding sites reported was previously identified experimentally, and another one, importantly, is identical to that of capsaicin, one of the chemical stimuli that activate TRPV1. However, in contrast to capsaicin, isoflurane and chloroform binding free-energies render modest to no association compared to capsaicin, suggesting a different activation mechanism. Uncovering chloroform and isoflurane modulatory sites will further our understanding of the TRPV1 molecular machinery and open the possibility of developing site-specific drugs.

A membrane-embedded pathway delivers general anesthetics to two interacting binding sites in the Gloeobacter violaceus ion channel

General anesthetics exert their effects on the central nervous system by acting on ion channels, most notably pentameric ligand-gated ion channels. Although numerous studies have focused on pentameric ligand-gated ion channels, the details of anesthetic binding and channel modulation are still debated. A better understanding of the anesthetic mechanism of action is necessary for the development of safer and more efficacious drugs. Herein, we present a computational study identifying two anesthetic binding sites in the transmembrane domain of the Gloeobacter violaceus ligand-gated ion channel (GLIC) channel, characterize the putative binding pathway, and observe structural changes associated with channel function. Molecular simulations of desflurane reveal a binding pathway to GLIC via a membrane-embedded tunnel using an intrasubunit protein lumen as the conduit, an observation that explains the Meyer-Overton hypothesis, or why the lipophilicity of an anesthetic and its potency are generally proportional. Moreover, employing high concentrations of ligand led to the identification of a second transmembrane site (TM2) that inhibits dissociation of anesthetic from the TM1 site and is consistent with the high concentrations of anesthetics required to achieve clinical effects. Finally, asymmetric binding patterns of anesthetic to the channel were found to promote an iris-like conformational change that constricts and dehydrates the ion pore, creating a 13.5 kcal/mol barrier to ion translocation. Together with previous studies, the simulations presented herein demonstrate a novel anesthetic binding site in GLIC that is accessed through a membrane-embedded tunnel and interacts with a previously known site, resulting in conformational changes that produce a non-conductive state of the channel.

The cellular membrane as a mediator for small molecule interaction with membrane proteins

The cellular membrane constitutes the first element that encounters a wide variety of molecular species to which a cell might be exposed. Hosting a large number of structurally and functionally diverse proteins associated with this key metabolic compartment, the membrane not only directly controls the traffic of various molecules in and out of the cell, it also participates in such diverse and important processes as signal transduction and chemical processing of incoming molecular species. In this article, we present a number of cases where details of interaction of small molecular species such as drugs with the membrane, which are often experimentally inaccessible, have been studied using advanced molecular simulation techniques. We have selected systems in which partitioning of the small molecule with the membrane constitutes a key step for its final biological function, often binding to and interacting with a protein associated with the membrane. These examples demonstrate that membrane partitioning is not only important for the overall distribution of drugs and other small molecules into different compartments of the body, it may also play a key role in determining the efficiency and the mode of interaction of the drug with its target protein. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Structural perturbation of a dipalmitoylphosphatidylcholine (DPPC) bilayer by warfarin and its bolaamphiphilic analogue: A molecular dynamics study

Compounds with nominally similar biological activity may exhibit differential toxicity due to differences in their interactions with cell membranes. Many pharmaceutical compounds are amphiphilic and can be taken up by phospholipid bilayers, interacting strongly with the lipid-aqueous interface whether or not subsequent permeation through the bilayer is possible. Bolaamphiphilic compounds, which possess two hydrophilic ends and a hydrophobic linker, can likewise undergo spontaneous uptake by bilayers. While membrane-spanning bolaamphiphiles can stabilize membranes, small molecules with this characteristic have the potential to create membrane defects via disruption of bilayer structure and dynamics. When compared to the amphiphilic therapeutic anticoagulant, warfarin, the bolaamphiphilic analogue, brodifacoum, exhibits heightened toxicity that goes beyond superior inhibition of the pharmacological target enzyme. We explore, herein, the consequences of anticoagulant accumulation in a dipalmitoylphosphatidylcholine (DPPC) bilayer. Coarse-grained molecular dynamics simulations reveal that permeation of phospholipid bilayers by brodifacoum causes a disruption of membrane barrier function that is driven by the bolaamphiphilic nature and size of this molecule. We find that brodifacoum partitioning into bilayers causes membrane thinning and permeabilization and promotes lipid flip-flop - phenomena that are suspected to play a role in triggering cell death. These phenomena are either absent or less pronounced in the case of the less toxic, amphiphilic compound, warfarin.

Modelling of noble anaesthetic gases and high hydrostatic pressure effects in lipid bilayers

Our objective was to study molecular processes that might be responsible for inert gas narcosis and high-pressure nervous syndrome. The classical molecular dynamics trajectories (200 ns) of dioleoylphosphatidylcholine (DOPC) bilayers simulated by the Berger force field were evaluated for water and the atomic distribution of noble gases around DOPC molecules in the pressure range of 1-1000 bar and at a temperature of 310 K. Xenon and argon have been tested as model gases for general anaesthetics, and neon has been investigated for distortions that are potentially responsible for neurological tremors in hyperbaric conditions. The analysis of stacked radial pair distribution functions of DOPC headgroup atoms revealed the explicit solvation potential of the gas molecules, which correlates with their dimensions. The orientational dynamics of water molecules at the biomolecular interface should be considered as an influential factor, while excessive solvation effects appearing in the lumen of membrane-embedded ion channels could be a possible cause of inert gas narcosis. All the noble gases tested exhibit similar order parameter patterns for both DOPC acyl chains, which are opposite of the patterns found for the order parameter curve at high hydrostatic pressures in intact bilayers. This finding supports the 'critical volume' hypothesis of anaesthesia pressure reversal. The irregular lipid headgroup-water boundary observed in DOPC bilayers saturated with neon in the pressure range of 1-100 bar could be associated with the possible manifestation of neurological tremors at the atomic scale. The non-immobiliser neon also demonstrated the highest momentum impact on the normal component of the DOPC diffusion coefficient representing the monolayer undulation rate, which indicates that enhanced diffusivity rather than atomic size is the key factor.

Spontaneous adsorption of coiled-coil model peptides K and E to a mixed lipid bilayer

A molecular description of the lipid-protein interactions underlying the adsorption of proteins to membranes is crucial for understanding, for example, the specificity of adsorption or the binding strength of a protein to a bilayer, or for characterizing protein-induced changes of membrane properties. In this paper, we extend an automated in silico assay (DAFT) for binding studies and apply it to characterize the adsorption of the model fusion peptides E and K to a mixed phospholipid/cholesterol membrane using coarse-grained molecular dynamics simulations. In addition, we couple the coarse-grained protocol to reverse transformation to atomistic resolution, thereby allowing to study molecular interactions with high detail. The experimentally observed differential binding of the peptides E and K to membranes, as well as the increased binding affinity of helical over unstructered peptides, could be well reproduced using the polarizable Martini coarse-grained (CG) force field. Binding to neutral membranes is shown to be dominated by initial binding of the positively charged N-terminus to the phospholipid headgroup region, followed by membrane surface-aligned insertion of the peptide at the interface between the hydrophobic core of the membrane and its polar headgroup region. Both coarse-grained and atomistic simulations confirm a before hypothesized snorkeling of lysine side chains for the membrane-bound state of the peptide K. Cholesterol was found to be enriched in peptide vicinity, which is probably of importance for the mechanism of membrane fusion. The applied sequential multiscale method, using coarse-grained simulations for the slow adsorption process of peptides to membranes followed by backward transformation to atomistic detail and subsequent atomistic simulations of the preformed peptide-lipid complexes, is shown to be a versatile approach to study the interactions of peptides or proteins with biomembranes.

Hydrocarbon molar water solubility predicts NMDA vs. GABAA receptor modulation

Background

Many anesthetics modulate 3-transmembrane (such as NMDA) and 4-transmembrane (such as GABA_A) receptors. Clinical and experimental anesthetics exhibiting receptor family specificity often have low water solubility. We hypothesized that the molar water solubility of a hydrocarbon could be used to predict receptor modulation in vitro.

Methods

GABA_A (α₁β₂γ_2s) or NMDA (NR1/NR2A) receptors were expressed in oocytes and studied using standard two-electrode voltage clamp techniques. Hydrocarbons from 14 different organic functional groups were studied at saturated concentrations, and compounds within each group differed only by the carbon number at the ω-position or within a saturated ring. An effect on GABA_A or NMDA receptors was defined as a 10% or greater reversible current change from baseline that was statistically different from zero.

Results

Hydrocarbon moieties potentiated GABA_A and inhibited NMDA receptor currents with at least some members from each functional group modulating both receptor types. A water solubility cut-off for NMDA receptors occurred at 1.1 mM with a 95% CI = 0.45 to 2.8 mM. NMDA receptor cut-off effects were not well correlated with hydrocarbon chain length or molecular volume. No cut-off was observed for GABA_A receptors within the solubility range of hydrocarbons studied.

Conclusions

Hydrocarbon modulation of NMDA receptor function exhibits a molar water solubility cut-off. Differences between unrelated receptor cut-off values suggest that the number, affinity, or efficacy of protein-hydrocarbon interactions at these sites likely differ.

Inhibition of voltage-gated Na⁺ channels by the synthetic cannabinoid ajulemic acid

BACKGROUND

The synthetic cannabinoid ajulemic acid has been demonstrated to alleviate pain in patients suffering from chronic neuropathic pain. Cannabinoids interact with several molecules within the pain circuit, including a potent inhibition of voltage-gated sodium channels. In this study, we closely characterized this property on neuronal and nonneuronal sodium channels.

METHODS

The inhibition of sodium inward currents by ajulemic acid was studied in vitro. Human embryonic kidney 293t cells were used as the expression system for Nav1.2, 1.3, 1.4, 1.5, 1.5N406K, 1.5F1760A, and 1.7; Nav1.8 was transiently expressed in ND7/23 cells. Nav1.2, Nav1.3, and Nav 1.8 were from rats, and Nav1.4, Nav1.5, and Nav1.7 were of human origin. Sodium currents were analyzed by means of the whole cell patch-clamp technique. The investigated concentrations of ajulemic acid were 0.1, 0.3, 1, 3, 10, and 30 μmol/L.

RESULTS

Ajulemic acid reversibly and concentration-dependently inhibited all voltage-gated sodium channel (Nav) isoforms investigated in this study, including Nav1.2, 1.3, 1.4, 1.5, 1.7, and 1.8. Tonic block of resting channels yielded half-maximal inhibitory concentration values between 2 and 9 μmol/L and was strongly enhanced on inactivated channels, suggesting state-dependent inhibition by ajulemic acid. Tonic block did not differ significantly when comparing Nav1.2 and Nav1.3, Nav1.4 and Nav1.5, and Nav1.7 and Nav1.8. Statistical analysis of other combinations of subunits (e.g., Nav1.2 and Nav1.4) by analysis of variance yielded a significant difference in block. Although we did not observe any relevant use-dependent block, ajulemic acid induced a strong hyperpolarizing shift of the voltage dependency of fast inactivation and modest shift of slow inactivation. The local anesthetic-insensitive Nav1.5 constructs N406K and F1760A displayed a preserved sensitivity to block by ajulemic acid. Finally, we found that low concentrations of ajulemic acid efficiently inhibited Navβ4 peptide-mediated resurgent currents in Nav1.5.

CONCLUSIONS

Our data suggest that block of sodium channels can be a relevant mechanism by which ajulemic acid alleviates neuropathic pain. The potent inhibition of resurgent currents and the preserved block on local anesthetic-insensitive channels indicates that ajulemic acid interacts with a conserved but yet unknown site of sodium channels.

Molecular Simulation of the Transport of Drugs across Model Membranes

This Perspective describes recent progress in the area of the molecular simulation of the interactions of hydrophobic and hydrophilic solutes with membranes. The ability to predict drug solubility prior to synthesis is an extremely desirable goal for pharmaceutical research. A major advantage of molecular dynamics is the ability to computationally probe both the drug solubility as well as the pathway for the transport of drugs across membranes. Computational methods to predict the interaction free energy of solutes with membranes have advanced significantly in recent years and can characterize the intra- and intermolecular state of the drug as well as perturbations of the drug to the membrane environment itself. In addition to a brief review on computational methods to characterize the transport of drugs across membranes, we will highlight recent discoveries and discuss the challenges and opportunities in the field.

Pentameric Ligand-gated Ion Channels : Insights from Computation

Pentameric ligand-gated ion channels (pLGICs) conduct upon the binding of an agonist and are fundamental to neurotransmission. New insights into the complex mechanisms underlying pLGIC gating, ion selectivity, and modulation have recently been gained via a series of crystal structures in prokaryotes and C .elegans, as well as computational studies relying on these structures. Here we review contributions from a variety of computational approaches, including normal mode analysis, automated docking, and fully atomistic molecular dynamics simulation. Examples from our own research, particularly concerning interactions with general anesthetics and lipids, are used to illustrate predictive results complementary to crystallographic studies.

Seeking structural specificity: direct modulation of pentameric ligand-gated ion channels by alcohols and general anesthetics

Alcohols and other anesthetic agents dramatically alter neurologic function in a wide range of organisms, yet their molecular sites of action remain poorly characterized. Pentameric ligand-gated ion channels, long implicated in important direct effects of alcohol and anesthetic binding, have recently been illuminated in renewed detail thanks to the determination of atomic-resolution structures of several family members from lower organisms. These structures provide valuable models for understanding and developing anesthetic agents and for allosteric modulation in general. This review surveys progress in this field from function to structure and back again, outlining early evidence for relevant modulation of pentameric ligand-gated ion channels and the development of early structural models for ion channel function and modulation. We highlight insights and challenges provided by recent crystal structures and resulting simulations, as well as opportunities for translation of these newly detailed models back to behavior and therapy.

Biophysical changes induced by xenon on phospholipid bilayers

Structural and dynamic changes in cell membrane properties induced by xenon, a volatile anesthetic molecule, may affect the function of membrane-mediated proteins, providing a hypothesis for the mechanism of general anesthetic action. Here, we use molecular dynamics simulation and differential scanning calorimetry to examine the biophysical and thermodynamic effects of xenon on model lipid membranes. Our results indicate that xenon atoms preferentially localize in the hydrophobic core of the lipid bilayer, inducing substantial increases in the area per lipid and bilayer thickness. Xenon depresses the membrane gel-liquid crystalline phase transition temperature, increasing membrane fluidity and lipid head group spacing, while inducing net local ordering effects in a small region of the lipid carbon tails and modulating the bilayer lateral pressure profile. Our results are consistent with a role for nonspecific, lipid bilayer-mediated mechanisms in producing xenon's general anesthetic action.

General anesthetics predicted to block the GLIC pore with micromolar affinity

Although general anesthetics are known to modulate the activity of ligand-gated ion channels in the Cys-loop superfamily, there is at present neither consensus on the underlying mechanisms, nor predictive models of this modulation. Viable models need to offer quantitative assessment of the relative importance of several identified anesthetic binding sites. However, to date, precise affinity data for individual sites has been challenging to obtain by biophysical means. Here, the likely role of pore block inhibition by the general anesthetics isoflurane and propofol of the prokaryotic pentameric channel GLIC is investigated by molecular simulations. Microscopic affinities are calculated for both single and double occupancy binding of isoflurane and propofol to the GLIC pore. Computations are carried out for an open-pore conformation in which the pore is restrained to crystallographic radius, and a closed-pore conformation that results from unrestrained molecular dynamics equilibration of the structure. The GLIC pore is predicted to be blocked at the micromolar concentrations for which inhibition by isofluorane and propofol is observed experimentally. Calculated affinities suggest that pore block by propofol occurs at signifcantly lower concentrations than those for which inhibition is observed: we argue that this discrepancy may result from binding of propofol to an allosteric site recently identified by X-ray crystallography, which may cause a competing gain-of-function effect. Affinities of isoflurane and propofol to the allosteric site are also calculated, and shown to be 3 mM for isoflurane and for propofol; both anesthetics have a lower affinity for the allosteric site than for the unoccupied pore.

Although general anesthesia is performed every day on thousands of people, its detailed microscopic mechanisms are not known. What is known is that general anesthetic drugs modulate the activity of ion channels in the central nervous system. These channels are proteins that open in response to binding of neurotransmitter molecules, creating an electric current through the cell membrane and thus propagating nerve impulses between cells. One possible mechanism for ion channel inhibition by anesthetics is that the drugs bind inside the pore of the channels, blocking ion current. Here we investigate such a pore block mechanism by computing the strength of the drugs' interaction with the pore – and hence the likelihood of binding, in the case of GLIC, a bacterial channel protein. The results, obtained from numerical simulations of atomic models of GLIC, indicate that the anesthetics isoflurane and propofol have a tendency to bind in the pore that is strong enough to explain blocking of the channel, even at low concentration of the drugs.

The Role of Atomic Polarization in the Thermodynamics of Chloroform Partitioning to Lipid Bilayers

In spite of extensive research and use in medical practice, the precise molecular mechanism of volatile anesthetic action remains unknown. The distribution of anesthetics within lipid bilayers and potential targeting to membrane proteins is thought to be central to therapeutic function. Therefore, obtaining a molecular level understanding of volatile anesthetic partitioning into lipid bilayers is of vital importance to modern pharmacology. In this study we investigate the partitioning of the prototypical anesthetic, chloroform, into lipid bilayers and different organic solvents using molecular dynamics simulations with potential models ranging from simplified coarse-grained MARTINI to additive and polarizable CHARMM all-atom force fields. Many volatile anesthetics display significant inducible dipole moments, which correlate with their potency, yet the exact role of molecular polarizability in their stabilization within lipid bilayers remains unknown. We observe that explicit treatment of atomic polarizability makes it possible to accurately reproduce solvation free energies in solvents with different polarities, allowing for quantitative studies in heterogeneous molecular distributions, such as lipid bilayers. We calculate the free energy profiles for chloroform crossing lipid bilayers to reveal a role of polarizability in modulating chloroform partitioning thermodynamics via the chloroform-induced dipole moment and highlight competitive binding to the membrane core and toward the glycerol backbone that may have significant implications for understanding anesthetic action.

Extension of a three-helix bundle domain of myosin VI and key role of calmodulins

The molecular motor protein myosin VI moves toward the minus-end of actin filaments with a step size of 30-36 nm. Such large step size either drastically limits the degree of complex formation between dimer subunits to leave enough length for the lever arms, or requires an extension of the lever arms' crystallographically observed structure. Recent experimental work proposed that myosin VI dimerization triggers the unfolding of the protein's proximal tail domain which could drive the needed lever-arm extension. Here, we demonstrate through steered molecular dynamics simulation the feasibility of sufficient extension arising from turning a three-helix bundle into a long α-helix. A key role is played by the known calmodulin binding that facilitates the extension by altering the strain path in myosin VI. Sequence analysis of the proximal tail domain suggests that further calmodulin binding sites open up when the domain's three-helix bundle is unfolded and that subsequent calmodulin binding stabilizes the extended lever arms.

Biophysical and computational studies of membrane penetration by the GRP1 pleckstrin homology domain

The pleckstrin homology (PH) domain of the general receptor for phosphoinositides 1 (GRP1) exhibits specific, high-affinity, reversible binding to phosphatidylinositol (3,4,5)-trisphosphate (PI(3,4,5)P(3)) at the plasma membrane, but the nature and extent of the interaction between this bound complex and the surrounding membrane environment remains unclear. Combining equilibrium and nonequilibrium molecular dynamics (MD) simulations, NMR spectroscopy, and monolayer penetration experiments, we characterize the membrane-associated state of GRP1-PH. MD simulations show loops flanking the binding site supplement the interaction with PI(3,4,5)P(3) through multiple contacts with the lipid bilayer. NMR data show large perturbations in chemical shift for these loop regions on binding to PI(3,4,5)P(3)-containing DPC micelles. Monolayer penetration experiments and further MD simulations demonstrate that mutating hydrophobic residues to polar residues in the flanking loops reduces membrane penetration. This supports a "dual-recognition" model of binding, with specific GRP1-PH-PI(3,4,5)P(3) interactions supplemented by interactions of loop regions with the lipid bilayer.

Generalized detailed fluctuation theorem under nonequilibrium feedback control

It has been shown recently that the Jarzynski equality is generalized under nonequilibrium feedback control [T. Sagawa and M. Ueda, Phys. Rev. Lett. 104, 090602 (2010)]. The presence of feedback control in physical systems should modify both the Jarzynski equality and the detailed fluctuation theorem [K. H. Kim and H. Qian, Phys. Rev. E 75, 022102 (2007)]. However, the generalized Jarzynski equality under forward feedback control has been proved by considering that the physical systems under feedback control should locally satisfy the detailed fluctuation theorem. We use the same formalism and derive the generalized detailed fluctuation theorem for nonequilibrium driven systems under feedback control. We find that the feedback control in a physical system should preserve the detailed fluctuation theorem if the system has the same feedback information measure in forward and reverse directions.

New insights into the molecular mechanisms of general anaesthetics

This paper provides new insights of how general anaesthetic research should be carried out in the future by an analysis of what we know, what we do not know and what we would like to know. I describe previous hypotheses on the mechanism of action of general anaesthetics (GAs) involving membranes and protein receptors. I provide the reasons why the GABA type A receptor, the NMDA receptor and the glycine receptor are strong candidates for the sites of action of GAs. I follow with a review on attempts to provide a mechanism of action, and how future research should be conducted with the help of physical and chemical methods.

Multiple binding sites for the general anesthetic isoflurane identified in the nicotinic acetylcholine receptor transmembrane domain

An extensive search for isoflurane binding sites in the nicotinic acetylcholine receptor (nAChR) and the proton gated ion channel from Gloebacter violaceus (GLIC) has been carried out based on molecular dynamics (MD) simulations in fully hydrated lipid membrane environments. Isoflurane introduced into the aqueous phase readily partitions into the lipid membrane and the membrane-bound protein. Specifically, isoflurane binds persistently to three classes of sites in the nAChR transmembrane domain: (i) An isoflurane dimer occludes the pore, contacting residues identified by previous mutagenesis studies; analogous behavior is observed in GLIC. (ii) Several nAChR subunit interfaces are also occupied, in a site suggested by photoaffinity labeling and thought to positively modulate the receptor; these sites are not occupied in GLIC. (iii) Isoflurane binds to the subunit centers of both nAChR alpha chains and one of the GLIC chains, in a site that has had little experimental targeting. Interpreted in the context of existing structural and physiological data, the present MD results support a multisite model for the mechanism of receptor-channel modulation by anesthetics.

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.

Modulation of KvAP unitary conductance and gating by 1-alkanols and other surface active agents

The actions of alcohols and anesthetics on ion channels are poorly understood. Controversy continues about whether bilayer restructuring is relevant to the modulatory effects of these surface active agents (SAAs). Some voltage-gated K channels (Kv), but not KvAP, have putative low affinity alcohol-binding sites, and because KvAP structures have been determined in bilayers, KvAP could offer insights into the contribution of bilayer mechanics to SAA actions. We monitored KvAP unitary conductance and macroscopic activation and inactivation kinetics in PE:PG/decane bilayers with and without exposure to classic SAAs (short-chain 1-alkanols, cholesterol, and selected anesthetics: halothane, isoflurane, chloroform). At levels that did not measurably alter membrane specific capacitance, alkanols caused functional changes in KvAP behavior including lowered unitary conductance, modified kinetics, and shifted voltage dependence for activation. A simple explanation is that the site of SAA action on KvAP is its entire lateral interface with the PE:PG/decane bilayer, with SAA-induced changes in surface tension and bilayer packing order combining to modulate the shape and stability of various conformations. The KvAP structural adjustment to diverse bilayer pressure profiles has implications for understanding desirable and undesirable actions of SAA-like drugs and, broadly, predicts that channel gating, conductance and pharmacology may differ when membrane packing order differs, as in raft versus nonraft domains.

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.

Distribution and dynamics of adamantanes in a lipid bilayer

The adamantanes are a class of compounds that have found use in the treatment of influenza A and Parkinson's disease, among others. The mode of action for influenza A is based on the adamantanes' interaction with the transmembrane M2 channel, whereas the treatment of Parkinson's disease is thought to relate to a channel block of N-methyl-D-aspartate receptors. An understanding of how these compounds interact with the lipid bilayer is thus of great interest. We used molecular-dynamics simulations to calculate the potential of mean force of adamantanes in a lipid bilayer. Our results demonstrate a preference for the interfacial region of the lipid bilayer for both protonated and deprotonated species, with the protonated species proving significantly more favorable. However, the protonated species have a large free-energy barrier in the center of the membrane. In contrast, there is no barrier (compared with aqueous solution) at the center of the bilayer for deprotonated species, suggesting that the permeant species is indeed the neutral form, as commonly assumed. We discuss the results with respect to proposed mechanisms of action and implications for drug-delivery in general.

Behavioral differences between phosphatidic acid and phosphatidylcholine in the presence of the nicotinic acetylcholine receptor

It has been found experimentally that negatively charged phosphatidic acid (PA) lipids and cholesterol molecules stabilize the nicotinic acetylcholine receptor (nAChR) in a functional resting state that can participate in an agonist-induced conformational change. In this study, we compare phosphatidylcholine (PC) and PA lipid behavior in the presence of the nAChR to determine why PC lipids do not support a functional nAChR. For lipids that are located within 1.0 nm of the protein, both PC and PA lipids have very similar order-parameter and bilayer-thickness values, which indicate that the annular lipid properties are protein-dependent. The most significant difference between the PC and PA bilayers is the formation of a lipid domain around the protein, which is visible in the PA bilayer but not the PC bilayer. This suggests that the PA domain may help stabilize the nAChR resting state. The PA lipids in the microdomain have a decreased order compared to a homogeneous PA bilayer and the lipids near the protein attempt to increase the free space in their vicinity by residing in multiple lateral planes.

The volatile anesthetic isoflurane perturbs conformational activation of integrin LFA-1 by binding to the allosteric regulatory cavity

The molecular and structural basis of anesthetic interactions with conformations and functionalities of cell surface receptors remains to be elucidated. We have demonstrated that the widely used volatile anesthetic isoflurane blocks the activation-dependent conformational conversion of integrin lymphocyte function associated antigen-1 (LFA-1), the major leukocyte cell adhesion molecule, to a high-affinity configuration. Perturbation of LFA-1 activation by isoflurane at clinically relevant concentrations leads to the inhibition of T-cell interactions with target cells as well as ligand-triggered intracellular signaling. Nuclear magnetic resonance spectroscopy reveals that isoflurane binds within a cavity in the LFA-1 ligand-binding domain, which is a previously identified drug-binding site for allosteric small-molecule antagonists that stabilize LFA-1 in a low-affinity conformation. These results provide a potential mechanism for the immunomodulatory properties of isoflurane.

Abeta peptide interactions with isoflurane, propofol, thiopental and combined thiopental with halothane: a NMR study

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.

Ecotoxicology of narcosis: stereoselectivity and potential target sites

The stereoselectivity of certain anesthetics is currently thought to be inconsistent with lipid theories of narcosis. The EC50-values of etomidate enantiomers for tadpole narcosis are now examined as a function of octanol/water partition coefficients, and enhancement factors for predicted over experimental EC50 values are determined from a calibration curve for non-selective narcosis. The unfavored S-(-)-enantiomers of etomidate and two analogues surprisingly still fulfill the Meyer-Overton rule. The R+-enantiomers of etomidate and four structural analogues are up to 34-fold more active than expected. The non-chiral anesthetic, propofol, is 8-fold more active than expected. It is concluded that there may be two pathways to tadpole narcosis: enhanced narcosis involving specific receptor binding sites and non-selective narcosis corresponding to the Meyer-Overton rule and operating on the lipid/protein interface.