Atomistic Molecular Dynamics Simulations of Propofol and Fentanyl in Phosphatidylcholine Lipid Bilayers

Anesthetic Mechanisms: Synergistic Interactions With Lipid Rafts and Voltage-Gated Sodium Channels

Despite successfully utilizing anesthetics for over 150 years, the mechanism of action remains relatively unknown. Recent studies have shown promising results, but due to the complex interactions between anesthetics and their targets, there remains a clear need for further mechanistic research. We know that lipophilicity is directly connected to anesthetic potency since lipid solubility relates to anesthetic partition into the membrane. However, clinically relevant concentrations of anesthetics do not significantly affect lipid bilayers but continue to influence various molecular targets. Lipid rafts are derived from liquid-ordered phases of the plasma membrane that contain increased concentrations of cholesterol and sphingomyelin and act as staging platforms for membrane proteins, including ion channels. Although anesthetics do not perturb membranes at clinically relevant concentrations, they have recently been shown to target lipid rafts. In this review, we summarize current research on how different types of anesthetics-local, inhalational, and intravenous-bind and affect both lipid rafts and voltage-gated sodium channels, one of their major targets, and how those effects synergize to cause anesthesia and analgesia. Local anesthetics block voltage-gated sodium channel pores while also disrupting lipid packing in ordered membranes. Inhalational anesthetics bind to the channel pore and the voltage-sensing domain while causing an increase in the number, size, and diameter of lipid rafts. Intravenous anesthetics bind to the channel primarily at the voltage-sensing domain and the selectivity filter, while causing lipid raft perturbation. These changes in lipid nanodomain structure possibly give proteins access to substrates that have translocated as a result of these structural alterations, resulting in lipid-driven anesthesia. Overall, anesthetics can impact channel activity either through direct interaction with the channel, indirectly through the lipid raft, or both. Together, these result in decreased sodium ion flux into the cell, disrupting action potentials and producing anesthetic effects. However, more research is needed to elucidate the indirect mechanisms associated with channel disruption through the lipid raft, as not much is known about anionic lipid products and their influence over voltage-gated sodium channels. Anesthetics' effect on S-palmitoylation, a promising mechanism for direct and indirect influence over voltage-gated sodium channels, is another auspicious avenue of research. Understanding the mechanisms of different types of anesthetics will allow anesthesiologists greater flexibility and more specificity when treating patients.

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.

Investigating the Permeation Mechanism of Typical Phthalic Acid Esters (PAEs) and Membrane Response Using Molecular Dynamics Simulations

Phthalic acid esters (PAEs) are typical environmental endocrine disrupters, interfering with the endocrine system of organisms at very low concentrations. The plasma membrane is the first barrier for organic pollutants to enter the organism, so membrane permeability is a key factor affecting their biological toxicity. In this study, based on computational approaches, we investigated the permeation and intramembrane aggregation of typical PAEs (dimethyl phthalate, DMP; dibutyl phthalate, DBP; di-2-ethyl hexyl phthalate, DEHP), as well as their effects on membrane properties, and related molecular mechanisms were uncovered. Our results suggested that PAEs could enter the membrane spontaneously, preferring the headgroup-acyl chain interface of the bilayer, and the longer the side chain (DEHP > DBP > DMP), the deeper the insertion. Compared with the shortest DMP, DEHP apparently increased membrane thickness, order, and rigidity, which might be due to its stronger hydrophobicity. Potential of means force (PMF) analysis revealed the presence of an energy barrier located at the water-membrane interface, with a maximum value of 2.14 kcal mol−1 obtained in the DEHP-system. Therefore, the difficulty of membrane insertion is also positively correlated with the side-chain length or hydrophobicity of PAE molecules. These findings will inspire our understanding of structure-activity relationship between PAEs and their effects on membrane properties, and provide a scientific basis for the formulation of environmental pollution standards and the prevention and control of small molecule pollutants.

Interaction With the Lipid Membrane Influences Fentanyl Pharmacology

Overdose deaths from fentanyl have reached epidemic proportions in the USA and are increasing worldwide. Fentanyl is a potent opioid agonist that is less well reversed by naloxone than morphine. Due to fentanyl’s high lipophilicity and elongated structure we hypothesised that its unusual pharmacology may be explained by its interactions with the lipid membrane on route to binding to the μ-opioid receptor (MOPr). Through coarse-grained molecular dynamics simulations, electrophysiological recordings and cell signalling assays, we determined how fentanyl and morphine access the orthosteric pocket of MOPr. Morphine accesses MOPr via the aqueous pathway; first binding to an extracellular vestibule, then diffusing into the orthosteric pocket. In contrast, fentanyl may take a novel route; first partitioning into the membrane, before accessing the orthosteric site by diffusing through a ligand-induced gap between the transmembrane helices. In electrophysiological recordings fentanyl-induced currents returned after washout, suggesting fentanyl deposits in the lipid membrane. However, mutation of residues forming the potential MOPr transmembrane access site did not alter fentanyl’s pharmacological profile in vitro. A high local concentration of fentanyl in the lipid membrane, possibly in combination with a novel lipophilic binding route, may explain the high potency and lower susceptibility of fentanyl to reversal by naloxone.

Effect of Ionic Strength on Ibuprofenate Adsorption on a Lipid Bilayer of Dipalmitoylphosphatidylcholine from Molecular Dynamics Simulations

In this work, the free energy change in the process of transferring ibuprofenate from the bulk solution to the center of a model of the dipalmitoylphosphatidylcholine bilayer at different NaCl concentrations was calculated. Two minima were found in the free energy profile: a local minimum, located in the vicinity of the membrane, and the global free energy minimum, found near the headgroup region. The downward shift of free energy minima with increasing NaCl concentration is consistent with the results of previous works. Conversely, the upward shift of the free energy maximum with increasing ionic strength is due to the competition of sodium ions and lipids molecules to coordinate with ibuprofenate and neutralize its charge. In addition, normal molecular dynamics simulations were performed to study the effects of the ibuprofenate on the lipid bilayer and in the presence of a high ibuprofenate concentration. The effect of ionic strength on the properties of the lipid bilayer and on lipid-drug interactions was analyzed. The area per lipid shrinking with increasing ionic strength, volume of lipids, and thickness of the bilayer is consistent with the experimental results. At a very high ibuprofenate concentration, the lipid bilayer dehydrates, and it consequently transforms into the gel phase, thus blocking the permeation.

Towards Quantum-Chemical Modeling of the Activity of Anesthetic Compounds

The modeling of the activity of anesthetics is a real challenge because of their unique electronic and structural characteristics. Microscopic approaches relevant to the typical features of these systems have been developed based on the advancements in the theory of intermolecular interactions. By stressing the quantum chemical point of view, here, we review the advances in the field highlighting differences and similarities among the chemicals within this group. The binding of the anesthetics to their partners has been analyzed by Symmetry-Adapted Perturbation Theory to provide insight into the nature of the interaction and the modeling of the adducts/complexes allows us to rationalize their anesthetic properties. A new approach in the frame of microtubule concept and the importance of lipid rafts and channels in membranes is also discussed.

Modulation of Phospholipid Bilayer Properties by Simvastatin

Simvastatin (Zocor) is one of the most prescribed drugs for reducing high cholesterol. Although simvastatin is ingested in its inactive lactone form, it is converted to its active dihydroxyheptanoate form by carboxylesterases in the liver. The dihydroxyheptanoate form can also be converted back to its original lactone form. Unfortunately, some of the side effects associated with the intake of simvastatin and other lipophilic statins at higher doses include statin-associated myopathy (SAM) and, in more severe cases, kidney failure. While the cause of SAM is unknown, it is hypothesized that these side effects are dependent on the localization of statins in lipid bilayers and their impact on bilayer properties. In this work, we carry out all-atom molecular dynamics simulations on both the lactone and dihydroxyheptanoate forms of simvastatin (termed "SN" and "SA", respectively) with a pure 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid bilayer and a POPC/cholesterol (30 mol %) binary mixture as membrane models. Additional simulations were carried out with multiple simvastatin molecules to mimic in vitro conditions that produced pleiotropic effects. Both SN and SA spontaneously diffused into the lipid bilayer, and a longer simulation time of 4 μs was needed for the complete incorporation of multiple SAs into the bilayer. By constructing potential mean force and electron density profiles, we find that SN localizes deeper within the hydrophobic interior of the bilayer and that SA has a greater tendency to form hydrogen-bonding interactions with neighboring water molecules and lipid headgroups. For the pure POPC bilayer, both SN and SA increase membrane order, while membrane fluidity increases for the POPC/cholesterol bilayer.

A novel small-molecule fatty acid synthase inhibitor with antitumor activity by cell cycle arrest and cell division inhibition

Fatty acid synthase (FASN), the key enzyme in de novo lipogenesis, is an attractive therapeutic target for diseases characterized by excessive lipid accumulation. Many FASN inhibitors have failed in the clinical trial phase, largely because of poor solubility and safety. In this study, we generated a novel small-molecule FASN inhibitor by structure-based virtual screening. PFI09, the lead compound, is easy to synthesize, and inhibits the lipid synthesis in OP9 mammalian cell line and Caenorhabditis elegans as well as the proliferation of several cancer cell lines via the blockade of FASN. Mechanistic investigations show that PFI09 induces S-phase arrest, cell division reduction and apoptosis. We also develop a chemically stable analog of PFI09, MFI03, which reduces the proliferation of PC3 tumor cells both in vitro and in vivo, without toxicity to mice. In summary, our data suggest that MFI03 is an effective FASN inhibitor and a promising antineoplastic drug candidate.

Unassisted N-acetyl-phenylalanine-amide transport across membrane with varying lipid size and composition: kinetic measurements and atomistic molecular dynamics simulation

Biological membranes are essential to preserve structural integrity and regulate functional properties through the permeability of nutrients, pharmaceutical drugs, and neurotransmitters of a living cell. The movement of acetylated and amidated phenylalanine (NAFA) across 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) membrane bilayers is investigated to probe physical transport. The rate of transport is measured experimentally applying parallel artificial membrane permeation assay (PAMPA). At the physiological temperature, 310 K, the measured time constants in the neutral pH were ∼6 h in DOPC and ∼3 h in POPC, while in a more acidic condition, at a pH 4.8, the time constants were ∼8 h in both lipids. Computationally, we have expanded our transport study of three aromatic dipeptides across a bilayer composed of DOPC¹⁸. In this study, we have examined the effects of lipid composition and bilayer size on the passive transport of NAFA by simulating the dipeptide in three different bilayers, with 50 DOPC lipids, 50 POPC lipids, and 40 POPC molecules. Specifically, atomistic molecular dynamics simulations with umbrella sampling were used to calculate the potential of mean force for the passive permeation of NAFA across the bilayers. Diffusion constants were then calculated by numerically solving the Smoluchowski equation. Permeability coefficients and mean first passage times were then calculated. Structural properties - Ramachandran plots, sidechain torsions, peptide insertion angles, radial distribution functions, and proximal peptide water molecules - were also examined to determine the effect of system size and lipid type. In terms of systems size, we observed a small decrease in the highest barrier of the potential of mean force and fewer sampled sidechain dihedral angle conformations with 40 versus 50 POPC lipids due to weaker membrane deformations within a smaller lipid bilayer. In terms of lipid type, DOPC contains two monounsaturated acyl chains compared to only one such acyl chain in POPC; therefore, DOPC bilayers are less ordered and more easily deformed, as seen by a much broader potential of mean force profile. The NAFA in DOPC lipid also transitioned to an internally hydrogen-bonded backbone conformation at lower membrane depths than in POPC. Similarly, as for other aromatic dipeptides, NAFA tends to insert into the membrane sidechain-first, remains mostly desolvated in the membrane center, and exhibits slow reorientations within the bilayer in both DOPC and POPC. With a joint experimental and computational study we have gained a new insight into the rate of transport and the underlying microscopic mechanism in different lipid bilayer conditions of the simplest hydrophobic aromatic dipeptide.Communicated by Ramaswamy H. Sarma.

In silico studies of the interactions between propofol and fentanyl using Gaussian accelerated molecular dynamics

Fentanyl is a potent opioid analgesic, which for decades has been used routinely in surgical and therapeutic applications. In addition to its analgesic properties, fentanyl also possesses anesthetic properties, which are not well understood. Fentanyl is used in the general anesthesia process to induce and maintain anesthesia in combination with the general anesthetic propofol, which fentanyl is known to potentiate. As the atomic-level mechanism behind the potentiation of propofol is unclear, we have used classical molecular dynamics simulations to study the interactions of these drugs with the Gloeobacter violaceus ion channel (GLIC). This ion channel has been identified as a target for many anesthetic drugs. We identified multiple binding sites using flooding style and Gaussian accelerated molecular dynamics (GaMD) simulations, showing fentanyl acting as a stabiliser that holds propofol within binding sites. Our extensive GaMD simulations were also able to show the pathway by which propofol blocks the channel pore, which has previously been suggested as a mechanism for ion channel modulation. General anesthesia is a multi-drug process and this study provides the first insight into the interactions between two different drugs in the anesthesia process in a relevant biological environment.Communicated by Ramaswamy H. Sarma.