Propofol modulates the lipid phase transition and localizes near the headgroup of membranes

Micro-Solvation of Propofol in Propylene Glycol-Water Binary Mixtures: Molecular Dynamics Simulation Studies

The water microstructure around propofol plays a crucial role in controlling their solubility in the binary mixture. The unusual nature of such a water microstructure can influence both translational and reorientational dynamics, as well as the water hydrogen bond network near propofol. We have carried out all-atom molecular dynamics simulations of five different compositions of the propylene glycol (PG)/water binary mixture containing propofol (PFL) molecules to investigate the differential behavior of water microsolvation shells around propofol, which is likely to control the propofol solubility. It is evident from the simulation snapshots for various compositions that the PG at high molecular ratio favors the water cluster and extended chainlike network that percolates within the PG matrix, where the propofol is in the dispersed state. We estimated that the radial distribution function indicates higher ordered water microstructure around propofol for high PG content, as compared to the lower PG content in the PG/water mixture. So, the hydrophilic PG regulates the stability of the water micronetwork around propofol and its solubility in the binary mixture. We observed that the translational and rotational mobility of water belonging to the propofol microsolvation shell is hindered for high PG content and relaxed toward the low PG molecular ratio in the PG/water mixture. It has been noticed that the structural relaxation of the hydrogen bond formed between the propofol and the water molecules present in the propofol microsolvation shell for all five compositions is found to be slower for high PG content and becomes faster on the way to low PG content in the mixture. Simultaneously, we calculated the intermittent residence time correlation function of the water molecules belonging to the microsolvation shell around the propofol for five different compositions and found a faster short time decay followed up with long time components. Again, the origin of such long time decay is primarily from the structural relaxation of the microsolvation shell around the propofol, where the high PG content shows the slower structural relaxation that turns faster as the PG content approaches to the other end of the compositions. So, our studies showed that the slower structural relaxation of the microsolvation shell around propofol for a high PG molecular ratio in the PG/water mixture correlate well with the extensive ordering of the water microstructure and restricted water mobility and facilitates the dissolution process of propofol in the binary mixture.

Magic mushroom extracts in lipid membranes

The active hallucinogen of magic mushrooms, psilocin, is being repurposed to treat nicotine addiction and treatment-resistant depression. Psilocin belongs to the tryptamine class of psychedelic compounds which include the hormone serotonin. It is believed that psilocin exerts its effect by binding to the serotonin 5-HT_2A receptor. However, recent in-vivo evidence suggests that psilocin may employ a different mechanism to exert its effects. Membrane-mediated receptor desensitization of neurotransmitter receptors is one such mechanism. We compare the impact of the neutral and charged versions of psilocin and serotonin on the properties of zwitterionic and anionic lipid membranes using molecular dynamics simulations and calorimetry. Both compounds partition to the lipid interface and induce membrane thinning. The tertiary amine in psilocin, as opposed to the primary amine in serotonin, limits psilocin's impact on the membrane although more psilocin partitions into the membrane than serotonin. Calorimetry corroborates that both compounds induce a classical melting point depression like anesthetics do. Our results also lend support to a membrane-mediated receptor-binding mechanism for both psilocin and serotonin and provide physical insights into subtle chemical changes that can alter the membrane-binding of psychedelic compounds.

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.

Partitioning of Seven Different Classes of Antibiotics into LPS Monolayers Supports Three Different Permeation Mechanisms through the Outer Bacterial Membrane

The outer membrane (OM) of Gram-negative (G-) bacteria presents a barrier for many classes of antibacterial agents. Lipopolysaccharide (LPS), present in the outer leaflet of the OM, is stabilized by divalent cations and is considered to be the major impediment for antibacterial agent permeation. However, the actual affinities of major antibiotic classes toward LPS have not yet been determined. In the present work, we use Langmuir monolayers formed from E. coli Re and Rd types of LPS to record pressure-area isotherms in the presence of antimicrobial agents. Our observations suggest three general types of interactions. First, some antimicrobials demonstrated no measurable interactions with LPS. This lack of interaction in the case of cefsulodin, a third-generation cephalosporin antibiotic, correlates with its low efficacy against G- bacteria. Ampicillin and ciprofloxacin also show no interactions with LPS, but in contrast to cefsulodin, both exhibit good efficacy against G- bacteria, indicating permeation through common porins. Second, we observe substantial intercalation of the more hydrophobic antibiotics, novobiocin, rifampicin, azithromycin, and telithromycin, into relaxed LPS monolayers. These largely repartition back to the subphase with monolayer compression. We find that the hydrophobic area, charge, and dipole all show correlations with both the mole fraction of antibiotic retained in the monolayer at the monolayer-bilayer equivalence pressure and the efficacies of these antibiotics against G- bacteria. Third, amine-rich gentamicin and the cationic antimicrobial peptides polymyxin B and colistin show no hydrophobic insertion but are instead strongly driven into the polar LPS layer by electrostatic interactions in a pressure-independent manner. Their intercalation stably increases the area per molecule (by up to 20%), which indicates massive formation of defects in the LPS layer. These defects support a self-promoted permeation mechanism of these antibiotics through the OM, which explains the high efficacy and specificity of these antimicrobials against G- bacteria.

In vitro alteration on erythrocytes mechanical properties by propofol, remifentanil and vecuronium

Several studies report flow disturbance and microcirculation disorders upon anesthesia treatment. These alterations are often related to blood rheology changes. In this work, it was attempted to make a detailed description of the alterations in erythrocyte mechanical properties by the action of propofol, remifentanil, and vecuronium. For this, an in vitro study was performed on red blood cell samples from healthy donors incubated with solutions of propofol (4 μg/mL whole blood), remifentanil (10 ng/mL plasma), and vecuronium (0.15 μg/mL plasma). Erythrocyte viscoelastic parameters were determined by octuplicate using a Reómetro Eritrocitario. Also, a Wilcoxon signed rank-test with Yates correction for continuity was performed to analyze the overall alteration in the mechanical properties of erythrocytes. Statistical analysis showed that the three studied anesthetics changed the erythrocyte mechanical properties at different parts of the membrane. These results would imply an interaction of these anesthetics with the erythrocyte membrane. Finally, this could conduce to alterations in microcirculation.

Effects of gabergic phenols on the dynamic and structure of lipid bilayers: A molecular dynamic simulation approach

γ-Aminobutyric acid (GABA) is the major inhibitory neurotransmitter in the vertebrate and invertebrate nervous system. GABAA receptors are activated by GABA and their agonists, and modulated by a wide variety of recognized drugs, including barbiturates, anesthetics, and benzodiazepines. The phenols propofol, thymol, chlorothymol, carvacrol and eugenol act as positive allosteric modulators on GABAA-R receptor. These GABAergic phenols interact with the lipid membrane, therefore, their anesthetic activity could be the combined result of their specific activity (with receptor proteins) as well as nonspecific interactions (with surrounding lipid molecules) modulating the supramolecular organization of the receptor environment. Therefore, we aimed to contribute to a description of the molecular events that occur at the membrane level as part of the mechanism of general anesthesia, using a molecular dynamic simulation approach. Equilibrium molecular dynamics simulations indicate that the presence of GABAergic phenols in a DPPC bilayer orders lipid acyl chains for carbons near the interface and their effect is not significant at the bilayer center. Phenols interacts with the polar interface of phospholipid bilayer, particularly forming hydrogen bonds with the glycerol and phosphate group. Also, potential of mean force calculations using umbrella sampling show that propofol partition is mainly enthalpic driven at the polar region and entropic driven at the hydrocarbon chains. Finally, potential of mean force indicates that propofol partition into a gel DPPC phase is not favorable. Our in silico results were positively contrasted with previous experimental data.

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.

Interactions between model cell membranes and the neuroactive drug propofol

Vibrational sum frequency spectroscopy (VSFS) complemented by surface pressure isotherm and neutron reflectometry (NR) experiments were employed to investigate the interactions between propofol, a small amphiphilic molecule that currently is the most common general anaesthetic drug, and phospholipid monolayers. A series of biologically relevant saturated phospholipids of varying chain length from C₁₈ to C₁₄ were spread on either pure water or propofol (2,6-bis(1-methylethyl)phenol) solution in a Langmuir trough, and the change in the molecular structure of the film, induced by the interaction with propofol, was studied with respect to the surface pressure. The results from the surface pressure isotherm experiments revealed that propofol, as long as it remains at the interface, enhances the fluidity of the phospholipid monolayer. The VSF spectra demonstrate that for each phospholipid the amount of propofol in the monolayer region decreases with increasing surface pressure. Such squeeze out is in contrast to the enhanced interactions that can be exhibited by more complex amphiphilic molecules such as peptides. At surface pressures of 22-25 mN m^-1, which are relevant for biological cell membranes, most of the propofol has been expelled from the monolayer, especially in the case of the C₁₆ and C₁₈ phospholipids that adopt a liquid condensed phase packing of its alkyl tails. At lower surface pressures of 5 mN m^-1, the effect of propofol on the structure of the alkyl tails is enhanced when the phospholipids are present in a liquid expanded phase. Specifically, for the C₁₆ phospholipid, NR data reveal that propofol is located exclusively in the head group region, which is rationalized in the context of previous studies. The results imply a non-homogeneous distribution of propofol in the plane of real cell membranes, which is an inference that requires urgent testing and may help to explain why such low concentration of the drug are required to induce general anaesthesia.

Selective effect of cell membrane on synaptic neurotransmission

Atomistic molecular dynamics simulations were performed with 13 non-peptidic neurotransmitters (NTs) in three different membrane environments. The results provide compelling evidence that NTs are divided into membrane-binding and membrane-nonbinding molecules. NTs adhere to the postsynaptic membrane surface whenever the ligand-binding sites of their synaptic receptors are buried in the lipid bilayer. In contrast, NTs that have extracellular ligand-binding sites do not have a similar tendency to adhere to the membrane surface. This finding is a seemingly simple yet important addition to the paradigm of neurotransmission, essentially dividing it into membrane-independent and membrane-dependent mechanisms. Moreover, the simulations also indicate that the lipid composition especially in terms of charged lipids can affect the membrane partitioning of NTs. The revised paradigm, highlighting the importance of cell membrane and specific lipids for neurotransmission, should to be of interest to neuroscientists, drug industry and the general public alike.

Membrane Affinity of Platensimycin and Its Dialkylamine Analogs

Membrane permeability is a desired property in drug design, but there have been difficulties in quantifying the direct drug partitioning into native membranes. Platensimycin (PL) is a new promising antibiotic whose biosynthetic production is costly. Six dialkylamine analogs of PL were synthesized with identical pharmacophores but different side chains; five of them were found inactive. To address the possibility that their activity is limited by the permeation step, we calculated polarity, measured surface activity and the ability to insert into the phospholipid monolayers. The partitioning of PL and the analogs into the cytoplasmic membrane of E. coli was assessed by activation curve shifts of a re-engineered mechanosensitive channel, MscS, in patch-clamp experiments. Despite predicted differences in polarity, the affinities to lipid monolayers and native membranes were comparable for most of the analogs. For PL and the di-myrtenyl analog QD-11, both carrying bulky sidechains, the affinity for the native membrane was lower than for monolayers (half-membranes), signifying that intercalation must overcome the lateral pressure of the bilayer. We conclude that the biological activity among the studied PL analogs is unlikely to be limited by their membrane permeability. We also discuss the capacity of endogenous tension-activated channels to detect asymmetric partitioning of exogenous substances into the native bacterial membrane and the different contributions to the thermodynamic force which drives permeation.

Cholesterol expels ibuprofen from the hydrophobic membrane core and stabilizes lamellar phases in lipid membranes containing ibuprofen

There is increasing evidence that common drugs, such as aspirin and ibuprofen, interact with lipid membranes. Ibuprofen is one of the most common over the counter drugs in the world, and is used for relief of pain and fever. It interacts with the cyclooxygenase pathway leading to inhibition of prostaglandin synthesis. From X-ray diffraction of highly oriented model membranes containing between 0 and 20 mol% ibuprofen, 20 mol% cholesterol, and dimyristoylphosphatidylcholine (DMPC), we present evidence for a non-specific interaction between ibuprofen and cholesterol in lipid bilayers. At a low ibuprofen concentrations of 2 mol%, three different populations of ibuprofen molecules were found: two in the lipid head group region and one in the hydrophobic membrane core. At higher ibuprofen concentrations of 10 and 20 mol%, the lamellar bilayer structure is disrupted and a lamellar to cubic phase transition was observed. In the presence of 20 mol% cholesterol, ibuprofen (at 5 mol%) was found to be expelled from the membrane core and reside solely in the head group region of the bilayers. 20 mol% cholesterol was found to stabilize lamellar membrane structure and the formation of a cubic phase at 10 and 20 mol% ibuprofen was suppressed. The results demonstrate that ibuprofen interacts with lipid membranes and that the interaction is strongly dependent on the presence of cholesterol.

Macroscopic and macromolecular specificity of alkylphenol anesthetics for neuronal substrates

We used a photoactive general anesthetic called meta-azi-propofol (AziPm) to test the selectivity and specificity of alkylphenol anesthetic binding in mammalian brain. Photolabeling of rat brain sections with [³H]AziPm revealed widespread but heterogeneous ligand distribution, with [³H]AziPm preferentially binding to synapse-dense areas compared to areas composed largely of cell bodies or myelin. With [³H]AziPm and propofol, we determined that alkylphenol general anesthetics bind selectively and specifically to multiple synaptic protein targets. In contrast, the alkylphenol anesthetics do not bind to specific sites on abundant phospholipids or cholesterol, although [³H]AziPm shows selectivity for photolabeling phosphatidylethanolamines. Together, our experiments suggest that alkylphenol anesthetic substrates are widespread in number and distribution, similar to those of volatile general anesthetics, and that multi-target mechanisms likely underlie their pharmacology.

Sites and functional consequence of VDAC-alkylphenol anesthetic interactions

General anesthetics have previously been shown to bind mitochondrial VDAC. Here, using a photoactive analog of the anesthetic propofol, we determined that alkylphenol anesthetics bind to Gly56 and Val184 on rat VDAC1. By reconstituting rat VDAC into planar bilayers, we determined that propofol potentiates VDAC gating with asymmetry at the voltage polarities; in contrast, propofol does not affect the conductance of open VDAC. Additional experiments showed that propofol also does not affect gramicidin A properties that are sensitive to lipid bilayer mechanics. Together, this suggests propofol affects VDAC function through direct protein binding, likely at the lipid-exposed channel surface, and that gating can be modulated by ligand binding to the distal ends of VDAC β-strands where Gly56 and Val184 are located.