Distribution of halothane in a dipalmitoylphosphatidylcholine bilayer from molecular dynamics calculations

Understanding the Molecular Mechanism of Anesthesia: Effect of General Anesthetics and Structurally Similar Non-Anesthetics on the Properties of Lipid Membranes

General anesthesia can be caused by various, chemically very different molecules, while several other molecules, many of which are structurally rather similar to them, do not exhibit anesthetic effects at all. To understand the origin of this difference and shed some light on the molecular mechanism of general anesthesia, we report here molecular dynamics simulations of the neat dipalmitoylphosphatidylcholine (DPPC) membrane as well as DPPC membranes containing the anesthetics diethyl ether and chloroform and the structurally similar non-anesthetics n-pentane and carbon tetrachloride, respectively. To also account for the pressure reversal of anesthesia, these simulations are performed both at 1 bar and at 600 bar. Our results indicate that all solutes considered prefer to stay both in the middle of the membrane and close to the boundary of the hydrocarbon domain, at the vicinity of the crowded region of the polar headgroups. However, this latter preference is considerably stronger for the (weakly polar) anesthetics than for the (apolar) non-anesthetics. Anesthetics staying in this outer preferred position increase the lateral separation between the lipid molecules, giving rise to a decrease of the lateral density. The lower lateral density leads to an increased mobility of the DPPC molecules, a decreased order of their tails, an increase of the free volume around this outer preferred position, and a decrease of the lateral pressure at the hydrocarbon side of the apolar/polar interface, a change that might well be in a causal relation with the occurrence of the anesthetic effect. All these changes are clearly reverted by the increase of pressure. Furthermore, non-anesthetics occur in this outer preferred position in a considerably smaller concentration and hence either induce such changes in a much weaker form or do not induce them at all.

Alteration of lipid bilayer mechanics by volatile anesthetics: Insights from μs-long molecular dynamics simulations

Very few drugs in clinical practice feature the chemical diversity, narrow therapeutic window, unique route of administration, and reversible cognitive effects of volatile anesthetics. The correlation between their hydrophobicity and their potency and the increasing amount of evidence suggesting that anesthetics exert their action on transmembrane proteins, justifies the investigation of their effects on phospholipid bilayers at the molecular level, given the strong functional and structural link between transmembrane proteins and the surrounding lipid matrix. Molecular dynamics simulations of a model lipid bilayer in the presence of ethylene, desflurane, methoxyflurane, and the nonimmobilizer 1,2-dichlorohexafluorocyclobutane (also called F6 or 2N) at different concentrations highlight the structural consequences of VA partitioning in the lipid phase, with a decrease of lipid order and bilayer thickness, an increase in overall lipid lateral mobility and area-per-lipid, and a marked reduction in the mechanical stiffness of the membrane, that strongly correlates with the compounds' hydrophobicity.

•

Molecular simulations of lipid bilayer interaction with volatile anesthetics

•

Comparison of volatile anesthetics' and nonimmobilizers' effects on lipid bilayers

•

Ligand-dependent partitioning of the compounds in the lipid phase

•

Effects on bilayer thickness, stiffness, and lipid order upon ligand partitioning

Molecular dynamics simulation study of adsorption of anionic–nonionic surfactants at oil/water interfaces

Four anionic–nonionic surfactants with the same headgroups and different units of oxygen ethyl (EO) and oxygen propyl (PO) were adopted to investigate the influence on oil/water interfacial tensions in this article. Molecular dynamics (MD) simulations were conducted to study the interfacial property of the four surfactants. Four parameters were proposed to reveal the effecting mechanism of molecular structure on interfacial tension, which included the interfacial thickness, order parameter of the hydrophobic chain, radial distribution function, and the solvent accessible surface area. In addition, the electrostatic potential of the four surfactants was calculated. The research results indicated that the interface facial mask formed by the surfactants, which contained three EO or three PO units was more stable, and it was easier for the surfactants of six EO or six PO units to form a microemulsion at higher concentrations. The adsorption mechanism of the anionic–nonionic surfactant systems at the oil/water interfaces was supplemented at a molecular level, which provided fundamental guidance for an in-depth understanding of the optimal selection of the surfactants in enhancing oil recovery.

Four anionic–nonionic surfactants with the same headgroups and different units of oxygen ethyl (EO) and oxygen propyl (PO) were adopted to investigate the influence on oil/water interfacial tensions in this article.

Effect of general anesthetics on the properties of lipid membranes of various compositions

Computer simulations of four lipid membranes of different compositions, namely neat DPPC and PSM, and equimolar DPPC-cholesterol and PSM-cholesterol mixtures, are performed in the presence and absence of the general anesthetics diethylether and sevoflurane both at 1 and 600 bar. The results are analyzed in order to identify membrane properties that are potentially related to the molecular mechanism of anesthesia, namely that change in the same way in any membrane with any anesthetics, and change oppositely with increasing pressure. We find that the lateral lipid density satisfies both criteria: it is decreased by anesthetics and increased by pressure. This anesthetic-induced swelling is attributed to only those anesthetic molecules that are located close to the boundary of the apolar phase. This lateral expansion is found to lead to increased lateral mobility of the lipids, an effect often thought to be related to general anesthesia; to an increased fraction of the free volume around the outer preferred position of anesthetics; and to the decrease of the lateral pressure in the nearby range of the ester and amide groups, a region into which anesthetic molecules already cannot penetrate. All these changes are reverted by the increase of pressure. Another important finding of this study is that cholesterol has an opposite effect on the membrane properties than anesthetics, and, correspondingly, these changes are less marked in the presence of cholesterol. Therefore, changes in the membrane that can lead to general anesthesia are expected to occur in the membrane domains of low cholesterol content.

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.

Ligand partitioning into lipid bilayer membranes under high pressure: Implication of variation in phase-transition temperatures

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

Physical insights into the blood-brain barrier translocation mechanisms

The number of individuals suffering from diseases of the central nervous system (CNS) is growing with an aging population. While candidate drugs for many of these diseases are available, most of these pharmaceutical agents cannot reach the brain rendering most of the drug therapies that target the CNS inefficient. The reason is the blood-brain barrier (BBB), a complex and dynamic interface that controls the influx and efflux of substances through a number of different translocation mechanisms. Here, we present these mechanisms providing, also, the necessary background related to the morphology and various characteristics of the BBB. Moreover, we discuss various numerical and simulation approaches used to study the BBB, and possible future directions based on multi-scale methods. We anticipate that this review will motivate multi-disciplinary research on the BBB aiming at the design of effective drug therapies.

Analysis of pulmonary surfactant by Fourier transform infrared spectroscopy after exposure to sevoflurane and isoflurane

Pulmonary surfactant, consisting primarily of phospholipids and four surfactant-specific proteins, is among the first structures that is exposed to inhalation anesthetics. Consequently, changes of pulmonary surfactant due to this exposure could cause respiratory complications after long anesthetic procedures. Fourier transform infrared (FTIR) spectroscopy was used to explore the effects of two inhalation anesthetics, sevoflurane and isoflurane, on a commercially available pulmonary surfactant. The research was primarily focused on the effect of anesthetics on the lipid component of the surfactant. Four different concentrations of anesthetics were added, and the doses were higher from the low clinical doses typically used. Recorded spectra were analyzed using principal component analysis, and the Student's t-test was performed to confirm the results. The exposure to both anesthetics induced similar changes, consistent with the increase of the anesthetic concentration. The most pronounced effect was on the hydrophilic head group of phospholipids, which is in agreement with the disruption of the hydrogen bond, caused by the anesthetics. A change in the band intensities of CH2 stretching vibrations, indicative of a disordering effect of anesthetics on the hydrophobic tails of phospholipids, was also observed. Changes induced by isoflurane appear to be more pronounced than those induced by sevoflurane. Furthermore, our results suggest that FTIR spectroscopy is a promising tool in studying anesthetic effects on pulmonary surfactant.

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.

Study of procaine and tetracaine in the lipid bilayer using molecular dynamics simulation

Despite available experimental results, the molecular mechanism of action of local anesthetics upon the nervous system and contribution of the cell membrane to the process are still controversial. In this work, molecular dynamics simulations were performed to investigate the effect of two clinically used local anesthetics, procaine and tetracaine, on the structure and dynamics of a fully hydrated dimyristoylphosphatidylcholine lipid bilayer. We focused on comparing the main effects of uncharged and charged drugs on various properties of the lipid membrane: mass density distribution, diffusion coefficient, order parameter, radial distribution function, hydrogen bonding, electrostatic potential, headgroup angle, and water dipole orientation. To compare the diffusive nature of anesthetic through the lipid membrane quantitatively, we investigated the hexadecane/water partition coefficient using expanded ensemble simulation. We predicted the permeability coefficient of anesthetics in the following order: uncharged tetracaine > uncharged procaine > charged tetracaine > charged procaine. We also shown that the charged forms of drugs are more potent in hydrogen bonding, disturbing the lipid headgroups, changing the orientation of water dipoles, and increasing the headgroup electrostatic potential more than uncharged drugs, while the uncharged drugs make the lipid diffusion faster and increase the tail order parameter. The results of these simulation studies suggest that the different forms of anesthetics induce different structural modifications in the lipid bilayer, which provides new insights into their molecular mechanism.

Naratriptan aggregation in lipid bilayers: perspectives from molecular dynamics simulations

In order to understand the interaction between naratriptan and a fully hydrated bilayer of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidyl-choline (POPC), we carried out molecular dynamics simulations. The simulations were performed considering neutral and protonated ionization states, starting from different initial conditions. At physiological pH, the protonated state of naratriptan is predominant. It is expected that neutral compounds could have larger membrane partition than charged compounds. However, for the specific case of triptans, it is difficult to study neutral species in membranes experimentally, making computer simulations an interesting tool. When the naratriptan molecules were originally placed in water, they partitioned between the bilayer/water interface and water phase, as has been described for similar compounds. From this condition, the drugs displayed low access to the hydrophobic environment, with no significant effects on bilayer organization. The molecules anchored in the interface, due mainly to the barrier function of the polar and oriented lipid heads. On the other hand, when placed inside the bilayer, both neutral and protonated naratriptan showed self-aggregation in the lipid tail environment. In particular, the protonated species exhibited a pore-like structure, dragging water through this environment. Graphical Abstract Different behaviour of Naratriptan and Sumatriptan, when the drugs were originally placed in the lipid core.

The effect of anaesthetics on the properties of a lipid membrane in the biologically relevant phase: a computer simulation study

Phospholipid membranes containing four different general anaesthetic molecules are simulated in the biologically relevant L_α phase at atmospheric and high pressures.

Incorporation of ubiquinone in supported lipid bilayers on ITO

Ubiquinone (UQ) is one of the main electron and proton shuttle molecules in biological systems, and dipalmitoylphosphatidylcholine (DPPC) is one of the most used model lipids. Supported planar bilayers (SPBs) are extensively accepted as biological model membranes. In this study, SPBs have been deposited on ITO, which is a semiconductor with good electrical and optical features. Specifically, topographic atomic force microscopy (AFM) images and force curves have been performed on SPBs with several DPPC:UQ ratios to study the location and the interaction of UQ in the SPB. Additionally, cyclic voltammetry has been used to understand the electrochemical behavior of DPPC:UQ SPBs. Obtained results show that, in our case, UQ is placed in two main different positions in SPBs. First, between the DPPC hydrophobic chains, fact that originates a decrease in the breakthrough force of the bilayer, and the second between the two leaflets that form the SPBs. This second position occurs when increasing the UQ content, fact that eventually forms UQ aggregates at high concentrations. The formation of aggregates produces an expansion of the SPB average height and a bimodal distribution of the breakthrough force. The voltammetric response of UQ depends on its position on the bilayer.

Exploring the Free Energy Landscape of Solutes Embedded in Lipid Bilayers

Free energy calculations are vital for our understanding of biological processes on an atomistic scale and can offer insight to various mechanisms. However, in some cases, degrees of freedom (DOFs) orthogonal to the reaction coordinate have high energy barriers and/or long equilibration times, which prohibit proper sampling. Here we identify these orthogonal DOFs when studying the transfer of a solute from water to a model membrane. Important DOFs are identified in bulk liquids of different dielectric nature with metadynamics simulations and are used as reaction coordinates for the translocation process, resulting in two- and three-dimensional space of reaction coordinates. The results are in good agreement with experiments and elucidate the pitfalls of using one-dimensional reaction coordinates. The calculations performed here offer the most detailed free energy landscape of solutes embedded in lipid bilayers to date and show that free energy calculations can be used to study complex membrane translocation phenomena.

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.

Diffusive nature of xenon anesthetic changes properties of a lipid bilayer: molecular dynamics simulations

Effects of general anesthesia can be controllable by the ambient pressure. We perform molecular dynamics simulations for a 1-palmitoyl-2-oleoyl phosphatidylethanolamine lipid bilayer with or without xenon molecules by changing the pressure to elucidate the mechanism of the pressure reversal of general anesthesia. According to the diffusive nature of xenon molecules in the lipid bilayer, a decrease in the orientational order of the lipid tails, an increase in the area and volume per lipid molecule, and an increase in the diffusivity of lipid molecules are observed. We show that the properties of the lipid bilayer with xenon molecules at high pressure come close to those without xenon molecules at 0.1 MPa. Furthermore, we find that xenon molecules are concentrated in the middle of the lipid bilayer at high pressures by the pushing effect and that the diffusivity of xenon molecules is suppressed. These results suggest that the pressure reversal originates from a jamming and suppression of the diffusivity of xenon molecules in lipid bilayers.

A novel approach to the investigation of passive molecular permeation through lipid bilayers from atomistic simulations

Predicting the permeability coefficient (P) of drugs permeating through the cell membrane is of paramount importance in drug discovery. We here propose an approach for calculating P based on bias-exchange metadynamics. The approach allows constructing from atomistic simulations a model of permeation taking explicitly into account not only the "trivial" reaction coordinate, the position of the drug along the direction normal to the lipid membrane plane, but also other degrees of freedom, for example, the torsional angles of the permeating molecule, or variables describing its solvation/desolvation. This allows deriving an accurate picture of the permeation process, and constructing a detailed molecular model of the transition state, making a rational control of permeation properties possible. We benchmarked this approach on the permeation of ethanol molecules through a POPC membrane, showing that the value of P calculated with our model agrees with the one calculated by a long unbiased molecular dynamics of the same system.

Fluorescence anisotropy of diphenylhexatriene and its cationic Trimethylamino derivative in liquid dipalmitoylphosphatidylcholine liposomes: opposing responses to isoflurane

Background

The mechanism of action of volatile general anesthetics has not yet been resolved. In order to identify the effects of isoflurane on the membrane, we measured the steady-state anisotropy of two fluorescent probes that reside at different depths. Incorporation of anesthetic was confirmed by shifting of the main phase transition temperature.

Results

In liquid crystalline dipalmitoylphosphatidylcholine liposomes, isoflurane (7-25 mM in the bath) increases trimethylammonium-diphenylhexatriene fluorescence anisotropy by ~0.02 units and decreases diphenylhexatriene anisotropy by the same amount.

Conclusions

The anisotropy data suggest that isoflurane decreases non-axial dye mobility in the headgroup region, while increasing it in the tail region. We propose that these results reflect changes in the lateral pressure profile of the membrane.

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.

In vivo optical virtual biopsy of human oral mucosa with harmonic generation microscopy

Recent clinical studies on human skin indicated that in vivo multi-harmonic generation microscopy (HGM) can achieve sub-micron resolution for histopathological analysis with a high penetration depth and leave no energy or photodamages in the interacted tissues. It is thus highly desired to apply HGM for in vivo mucosa histopathological diagnosis. In this paper, the first in vivo optical virtual biopsy of human oral mucosa by using epi-HGM is demonstrated. We modified an upright microscope to rotate the angle of objective for in vivo observation. Our clinical study reveals the capability of HGM to in vivo image cell distributions in human oral mucosa, including epithelium and lamina propria with a high penetration depth greater than 280 μm and a high spatial resolution better than 500 nm. We also found that the third-harmonic-generation (THG) contrast on nucleus depends strongly on its thicknesses, in agreement with a numerical simulation. Besides, 4% acetic acid was found to be able to enhance the THG contrast of nucleus in oral mucosa, while such enhancement was found to decay due to the metabolic clearance of the contrast enhancer by the oral mucosa. Our clinical study indicated that, the combined epi-THG and epi-second-harmonic-generation (SHG) microscopy is a promising imaging tool for in vivo noninvasive optical virtual biopsy and disease diagnosis in human mucosa.

Membrane fusion induced by small molecules and ions

Membrane fusion is a key event in many biological processes. These processes are controlled by various fusogenic agents of which proteins and peptides from the principal group. The fusion process is characterized by three major steps, namely, inter membrane contact, lipid mixing forming the intermediate step, pore opening and finally mixing of inner contents of the cells/vesicles. These steps are governed by energy barriers, which need to be overcome to complete fusion. Structural reorganization of big molecules like proteins/peptides, supplies the required driving force to overcome the energy barrier of the different intermediate steps. Small molecules/ions do not share this advantage. Hence fusion induced by small molecules/ions is expected to be different from that induced by proteins/peptides. Although several reviews exist on membrane fusion, no recent review is devoted solely to small moleculs/ions induced membrane fusion. Here we intend to present, how a variety of small molecules/ions act as independent fusogens. The detailed mechanism of some are well understood but for many it is still an unanswered question. Clearer understanding of how a particular small molecule can control fusion will open up a vista to use these moleucles instead of proteins/peptides to induce fusion both in vivo and in vitro fusion processes.

Partitioning of 2,6-Bis(1H-Benzimidazol-2-yl)pyridine fluorophore into a phospholipid bilayer: complementary use of fluorescence quenching studies and molecular dynamics simulations

Successful use of fluorescence sensing in elucidating the biophysical properties of lipid membranes requires knowledge of the distribution and location of an emitting molecule in the bilayer. We report here that 2,6-bis(1H-benzimidazol-2-yl)pyridine (BBP), which is almost non-fluorescent in aqueous solutions, reveals a strong emission enhancement in a hydrophobic environment of a phospholipid bilayer, making it interesting for fluorescence probing of water content in a lipid membrane. Comparing the fluorescence behavior of BBP in a wide variety of solvents with those in phospholipid vesicles, we suggest that the hydrogen bonding interactions between a BBP fluorophore and water molecules play a crucial role in the observed "light switch effect". Therefore, the loss of water-induced fluorescence quenching inside a membrane are thought to be due to deep penetration of BBP into the hydrophobic, water-free region of a bilayer. Characterized by strong quenching by transition metal ions in solution, BBP also demonstrated significant shielding from the action of the quencher in the presence of phospholipid vesicles. We used the increase in fluorescence intensity, measured upon titration of probe molecules with lipid vesicles, to estimate the partition constant and the Gibbs free energy (ΔG) of transfer of BBP from aqueous buffer into a membrane. Partitioning BBP revealed strongly favorable ΔG, which depends only slightly on the lipid composition of a bilayer, varying in a range from -6.5 to -7.0kcal/mol. To elucidate the binding interactions of the probe with a membrane on the molecular level, a distribution and favorable location of BBP in a POPC bilayer were modeled via atomistic molecular dynamics (MD) simulations using two different approaches: (i) free, diffusion-driven partitioning of the probe molecules into a bilayer and (ii) constrained umbrella sampling of a penetration profile of the dye molecule across a bilayer. Both of these MD approaches agreed with regard to the preferred location of a BBP fluorophore within the interfacial region of a bilayer, located between the hydrocarbon acyl tails and the initial portion of the lipid headgroups. MD simulations also revealed restricted permeability of water molecules into this region of a POPC bilayer, determining the strong fluorescence enhancement observed experimentally for the membrane-partitioned form of BBP.

Partitioning and localization of environment-sensitive 2-(2'-pyridyl)- and 2-(2'-pyrimidyl)-indoles in lipid membranes: a joint refinement using fluorescence measurements and molecular dynamics simulations

Fluorescence of environment-sensitive dyes is widely applied to monitor local structure and solvation dynamics of biomolecules. It has been shown that, in comparison with a parent indole fluorophore, fluorescence of 2-(2'-pyridyl)-5-methylindole (5M-PyIn-0) and 2-[2'-(4',6'-dimethylpyrimidyl)]-indole (DMPmIn-0) is remarkably sensitive to hydrogen bonding with protic partners. Strong fluorescence, observed for these compounds in nonpolar and polar aprotic solvents, is efficiently quenched in aqueous solution. This study demonstrates that 5M-PyIn-0 and DMPmIn-0, which are almost nonemitting in aqueous solution, become highly fluorescent upon titrating with phospholipid vesicles. The fluorescence enhancement is accompanied by a significant blue shift of emission maximum. The Gibbs free energy of membrane partitioning, measured by the increase in the steady-state fluorescence intensities during transfer from an aqueous environment to a lipid bilayer, is very favorable for both compounds, being in a range from -7.1 to -8.0 kcal/mol and depending only slightly on lipid composition of the membrane. The fluorescence enhancement upon membrane partitioning is indicative of the loss of the specific hydrogen-bonding interactions between the excited fluorophore and water molecules, causing efficient fluorescence quenching in bulk water. This conclusion is supported by atomistic molecular dynamics (MD) simulations, demonstrating that both 5M-PyIn-0 and DMPmIn-0 bind rapidly and partition deeply into a lipid bilayer. MD simulations also show a rapid, nanosecond-scale decrease in the probability of solute-solvent hydrogen bonding during passive diffusion of the probe molecules from bulk water into a lipid bilayer. At equilibrium conditions, both 5M-PyIn-0 and DMPmIn-0 prefer deep localization within the hydrophobic, water-free region of the bilayer. A free energy profile of penetration across a bilayer estimated using MD umbrella sampling shows that both indole derivatives favor residence in a rather wide potential energy well located 10-15 Å from the bilayer center.

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.

Anaesthetic mechanisms: update on the challenge of unravelling the mystery of anaesthesia

General anaesthesia is administered each day to thousands of patients worldwide. Although more than 160 years have passed since the first successful public demonstration of anaesthesia, a detailed understanding of the anaesthetic mechanism of action of these drugs is still lacking. An important early observation was the Meyer-Overton correlation, which associated the potency of an anaesthetic with its lipid solubility. This work focuses attention on the lipid membrane as a likely location for anaesthetic action. With the advent of cellular electrophysiology and molecular biology techniques, tools to dissect the components of the lipid membrane have led, in recent years, to the widespread acceptance of proteins, namely receptors and ion channels, as more likely targets for the anaesthetic effect. Yet these accumulated data have not produced a comprehensive explanation for how these drugs produce central nervous system depression. In this review, we follow the story of anaesthesia mechanisms research from its historical roots to the intensely neurophysiological research regarding it today. We will also describe recent findings that identify specific neuroanatomical locations mediating the actions of some anaesthetic agents.

Isovaleric, methylmalonic, and propionic acid decrease anesthetic EC50 in tadpoles, modulate glycine receptor function, and interact with the lipid 1,2-dipalmitoyl-Sn-glycero-3-phosphocholine

INTRODUCTION

Elevated concentrations of isovaleric (IVA), methylmalonic (MMA), and propionic acid are associated with impaired consciousness in genetic diseases (organic acidemias). We conjectured that part of the central nervous system depression observed in these disorders was due to anesthetic effects of these metabolites. We tested three hypotheses. First, that these metabolites would have anesthetic-sparing effects, possibly being anesthetics by themselves. Second, that these compounds would modulate glycine and gamma-aminobutyric acid (GABA(A)) receptor function, increasing chloride currents through these channels as potent clinical inhaled anesthetics do. Third, that these compounds would affect physical properties of lipids.

METHODS

Anesthetic EC(50)s were measured in Xenopus laevis tadpoles. Glycine and GABA(A) receptors were expressed in Xenopus laevis oocytes and studied using two-electrode voltage clamping. Pressure-area isotherms of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) monolayers were measured with and without added organic acids.

RESULTS

IVA acid was an anesthetic in tadpoles, whereas MMA and propionic acid decreased isoflurane's EC(50) by half. All three organic acids concentration-dependently increased current through alpha(1) glycine receptors. There were minimal effects on alpha(1)beta(2)gamma(2s) GABA(A) receptors. The organic acids increased total lateral pressure (surface pressure) of DPPC monolayers, including at mean molecular areas typical of bilayers.

CONCLUSION

IVA, MMA, and propionic acid have anesthetic effects in tadpoles, positively modulate glycine receptor function and affect physical properties of DPPC monolayers.

A hypothesis on the origin and evolution of the response to inhaled anesthetics

In this article, I present an evolutionary explanation for why organisms respond to inhaled anesthetics. It is conjectured that organisms today respond to inhaled anesthetics owing to the sensitivity of ion channels to inhaled anesthetics, which in turn has arisen by common descent from ancestral, anesthetic-sensitive ion channels in one-celled organisms (i.e., that the response to anesthetics did not arise as an adaptation of the nervous system, but rather of ion channels that preceded the origin of multicellularity). This sensitivity may have been refined by continuing selection at synapses in multicellular organisms. In particular, it is hypothesized that 1) the beneficial trait that was selected for in one-celled organisms was the coordinated response of ion channels to compounds that were present in the environment, which influenced the conformational equilibrium of ion channels; 2) this coordinated response prevented the deleterious consequences of entry of positive charges into the cell, thereby increasing the fitness of the organism; and 3) these compounds (which may have included organic anions, cations, and zwitterions as well as uncharged compounds) mimicked inhaled anesthetics in that they were interfacially active, and modulated ion channel function by altering bilayer properties coupled to channel function. The proposed hypothesis is consistent with known properties of inhaled anesthetics. In addition, it leads to testable experimental predictions of nonvolatile compounds having anesthetic-like modulatory effects on ion channels and in animals, including endogenous compounds that may modulate ion channel function in health and disease. The latter included metabolites that are increased in some types of end-stage organ failure, and genetic metabolic diseases. Several of these predictions have been tested and proved to be correct.

Distribution and favorable binding sites of pyrroloquinoline and its analogues in a lipid bilayer studied by molecular dynamics simulations

The distribution of 1H-pyrrolo[3,2-h]quinoline (PQ), 11H-dipyrido[2,3-a]carbazole (PC) and 7-azaindole (7AI) at a water/membrane interface has been investigated by molecular dynamics (MD) simulations. The MD study focused on favorable binding sites of the azaaromatic probes across a dipalmitoylphosphatidylcholine (DPPC) bilayer. Our simulations show that PQ and PC are preferably accommodated at the hydrocarbon core of the bilayer below the glycerol moiety. In addition, it is found that the hydrophobic aromatic parts of the probes are located inside a more ordered region of DPPC, consisting of hydrophobic lipid chains. In contrast to PQ and PC, 7AI is characterized by a broad distribution between a DPPC interface and water, so that the three preferable binding sites are found across a water/membrane interface. It is found that in the sequence 7AI-PQ-PC, due to the increase of the number of aromatic rings and, hence, the hydrophobic character of the probes, the depth of the probe localization is gradually shifted deeper inside the hydrocarbon core of the bilayer. We found that the probe-lipid hydrogen-bonding contributes weakly to the favorable localizations of the azaaromatic probes inside the DPPC bilayer, so that the probe localization is mainly driven by electrostatic dipole-dipole and van der Waals interactions.

Interaction of anesthetics with open and closed conformations of a potassium channel studied via molecular dynamics and normal mode analysis

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