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.

Pollution caused by nanoplastics: adverse effects and mechanisms of interaction via molecular simulation

The continuous increase in the production of synthetic plastics for decades and the inadequate disposal of plastic waste have resulted in a considerable increase of these materials in aquatic environments, which has developed into a major environmental concern. In addition to conventional parameters, the relevance of the environmental monitoring of microplastics (MPs) and nanoplastics (NPs) has been highlighted by the scientific community due to the potential adverse effects these materials pose to the ecosystem as well as to human health. The literature has registered an increasing interest in understanding the mechanisms, at the molecular level, of the interaction between NPs and other compounds using molecular simulation techniques. The present review aims to: (i) summarize the force fields conventionally used to describe NPs by molecular simulations; (ii) discuss the effects of NPs in the structural and dynamical properties of biological membranes; (iii) evaluate how NPs affect the folding of proteins; (iv) discuss the mechanisms by which NPs adsorb contaminants from the environment. NPs can affect the secondary structure of proteins and change the lateral organization and diffusion of lipid membranes. As a result, they may alter the lipid digestion in the gastrointestinal system representing a risk to the assimilation of the nutrients by humans. The adsorption of contaminants on MPs and NPs can potentiate their harmful effects on human health, due to a possible synergism. Therefore, understanding the mechanisms involved in these interactions is crucial to predict dangerous combinations and outline action strategies that reduce negative impacts on ecosystems and human health. Depending on the chemical properties of contaminants and NPs, electrostatic and/or van der Waals interactions can be more relevant in explaining the adsorption process. Finally, we conclude by highlighting gaps in the literature and the critical aspects for future investigations.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

UPLC-MS-Based Serum Metabolic Profiling Reveals Potential Biomarkers for Predicting Propofol Responsiveness in Females

Although previous studies have shown that certain factors interfere with the sensitivity of propofol, the mechanisms for interindividual variability in response to propofol remain unclear. This study aimed to screen the metabolites to predict patients' sensitivity to propofol and to identify metabolic pathways to explore possible mechanisms associated with propofol resistance. Sera from 40 female patients undergoing elective hysteroscopic surgery in a prospective cohort propofol study were obtained before the administration of propofol. The patients' responsiveness to propofol was differentiated based on propofol effect-site concentration. Serum samples from two sets, a discovery set (n = 24) and an independent validation set (n = 16), were analyzed using ultraperformance liquid chromatography coupled with mass spectrometry based untargeted metabolomics. In the discovery set, 494 differential metabolites were screened out, and then 391 potential candidate biomarkers with the area under receiver operating characteristic curve >0.80 were selected. Pathway analysis showed that the pathway of glycerophospholipid metabolism was the most influential pathway. In the independent validation set, six potential biomarkers enabled the discrimination of poor responders from good and intermediate responders, which might be applied to predict propofol sensitivity. The mass spectrometry data are available via MetaboLights (http://www.ebi.ac.uk/metabolights/login) with the identifier MTBLS2311.

A relationship between membrane permeation and partitioning of nitroaromatic explosives and their functional groups. A computational study

Nitroaromatic explosives, such as 2,4,6-trinitrotoluene, are representative aromatic compounds, which are generally highly toxic.

Dynamic structure factor of a lipid bilayer in the presence of a high electric field

The influence of a high average electric field (∼1 V/nm) in the hydrophobic interior of a bilayer lipid membrane on short-wavelength in-plane phononic motions of lipid chains is considered. The average electric field is assumed to be nearly constant on a picosecond time scale and a nanometer length scale. This field may be induced, for instance, by externally applied subnanosecond electric pulses or the membrane dipole potential. Using a generalized hydrodynamic approach, we derive a corresponding electrohydrodynamic model generalized to high wave numbers. In the considered approximation, all electric field effects are reduced only to a constant contribution to the generalized isothermal compressibility modulus. The corresponding dynamic structure factor for a lipid bilayer is derived. We show that due to polarization effects, the high field can critically impact the dynamics of longitudinal acousticlike modes at wave numbers near the major peak of the static structure factor. We estimate quantitatively that for typical lipid bilayers, transverse high electric fields can cause strong phonon energy softening, enhancement of phonon population, and formation of a gap in the dispersion of excitation frequency. The results obtained agree with simulations of the initiation of lipid bilayer electropores, suggesting that the proposed model reproduces the essential features of the field's impact on atomic density fluctuations. The proposed mechanism may have significant implications for the understanding of electroporation, passive molecular transport, and spontaneous pore formation in lipid bilayers.

Can Nanoplastics Alter Cell Membranes?

Whilst the formation of plastic nanoparticles (nanoplastics) from plastic wastes has been unequivocally evidenced, little is known about the effects of these materials on living organisms at the subcellular or molecular levels. In the present contribution we show through molecular dynamics simulations that polyethylene nanoparticles dissolve in the hydrophobic core of lipid bilayers into a network of disentangled, single polymeric chains. The thereby induced structural and dynamic changes in the bilayer alter vital functions of the cell membrane, which if lacking a mechanism to decompose the polymer chains may result in the death of the cell.

Response of microbial membranes to butanol: interdigitation vs. disorder

Biobutanol production by fermentation is potentially a sustainable alternative to butanol production from fossil fuels. However, the toxicity of butanol to fermentative bacteria, resulting largely from cell membrane fluidization, limits production titers and is a major factor limiting the uptake of the technology. Here, studies were undertaken, in vitro and in silico, on the butanol effects on a representative bacterial (i.e. Escherichia coli) inner cell membrane. A critical butanol : lipid ratio for stability of 2 : 1 was observed, computationally, consistent with complete interdigitation. However, at this ratio the bilayer was ∼20% thicker than for full interdigitation. Furthermore, butanol intercalation induced acyl chain bending and increased disorder, measured as a 27% lateral diffusivity increase experimentally in a supported lipid bilayer. There was also a monophasic Tm reduction in butanol-treated large unilamellar vesicles. Both behaviours are inconsistent with an interdigitated gel. Butanol thus causes only partial interdigitation at physiological temperatures, due to butanol accumulating at the phospholipid headgroups. Acyl tail disordering (i.e. splaying and bending) fills the subsequent voids. Finally, butanol short-circuits the bilayer and creates a coupled system where interdigitated and splayed phospholipids coexist. These findings will inform the design of strategies targeting bilayer stability for increasing biobutanol production titers.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

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.

Propagation of a thermo-mechanical perturbation on a lipid membrane

The propagation of sound waves on lipid monolayers supported on water has been previously studied during the melting transition. Since changes in volume, area, and compressibility in lipid membranes have biological relevance, the observed sound propagation is of paramount importance. However, it is unknown what would occur on a lipid bilayer, which is a more approximate model of a cell membrane. With the aim to answer this relevant question, we built an experimental setup to assemble long artificial lipid membranes. We found that if these membranes are heated in order to force local melting, a thermo-mechanical perturbation propagates a long distance. Our findings may support the existence of solitary waves, postulated to explain the propagation of isentropic signals together with the action potential in neurons.

DPPC-cholesterol phase diagram using coarse-grained Molecular Dynamics simulations

Cholesterol-phospholipid bilayers continue to be the current state of the art in membrane models and serve as representative systems for studying the effect of cholesterol on the cell membrane. As the mixing of different lipid species requires long spatio-temporal scales, coarse-grained models have gained increasing popularity in modeling such membrane systems. In this paper, a systematic study of the MARTINI coarse-grained model for the DPPC-cholesterol binary system has been performed. We construct the phase diagram of DPPC lipid bilayers in the presence of different cholesterol concentrations and at different temperatures using coarse-grained Molecular Dynamics (MD) simulations with the MARTINI force field. The phase diagram based on the condensation effect is directly comparable to available experimental data and demonstrates qualitative agreement over all cholesterol concentrations. Self-assembled bilayers quantitatively reproduce experimental observables, such as lateral diffusion of lipids, electron density, area per lipid and lipid order parameters. The phase diagram of the DPPC-cholesterol binary system also reveals the profound effect of cholesterol on the physical properties of phospholipid bilayers such lipid order, diffusion, and fluidity. Cholesterol induces the liquid-ordered phase, which increases the fluidity of the phospholipid hydrocarbon chains above the gel to liquid-crystalline phase transition temperature and decreases it below the phase transition. The present study suggests that the MARTINI force field can be successfully used to obtain molecular level insights into cholesterol-DPPC model membranes.

In silico pharmacology: Drug membrane partitioning and crossing

Over the past decade, molecular dynamics (MD) simulations have become particularly powerful to rationalize drug insertion and partitioning in lipid bilayers. MD simulations efficiently support experimental evidences, with a comprehensive understanding of molecular interactions driving insertion and crossing. Prediction of drug partitioning is discussed with respect to drug families (anesthetics; β-blockers; non-steroidal anti-inflammatory drugs; antioxidants; antiviral drugs; antimicrobial peptides). To accurately evaluate passive permeation coefficients turned out to be a complex theoretical challenge; however the recent methodological developments based on biased MD simulations are particularly promising. Particular attention is paid to membrane composition (e.g., presence of cholesterol), which influences drug partitioning and permeation. Recent studies concerning in silico models of membrane proteins involved in drug transport (influx and efflux) are also reported here. These studies have allowed gaining insight in drug efflux by, e.g., ABC transporters at an atomic resolution, explicitly accounting for the mandatory forces induced by the surrounded lipid bilayer. Large-scale conformational changes were thoroughly analyzed.

Exploring the Effects on Lipid Bilayer Induced by Noble Gases via Molecular Dynamics Simulations

Noble gases seem to have no significant effect on the anesthetic targets due to their simple, spherical shape. However, xenon has strong narcotic efficacy and can be used clinically, while other noble gases cannot. The mechanism remains unclear. Here, we performed molecular dynamics simulations on phospholipid bilayers with four kinds of noble gases to elucidate the difference of their effects on the membrane. Our results showed that the sequence of effects on membrane exerted by noble gases from weak to strong was Ne, Ar, Kr and Xe, the same order as their relative narcotic potencies as well as their lipid/water partition percentages. Compared with the other three kinds of noble gases, more xenon molecules were distributed between the lipid tails and headgroups, resulting in membrane’s lateral expansion and lipid tail disorder. It may contribute to xenon’s strong anesthetic potency. The results are well consistent with the membrane mediated mechanism of general anesthesia.