Thermodynamic study of benzocaine insertion into different lipid bilayers

Influence of Lipid Composition on the Insertion Process of Glyphosate into Membranes: A Thermodynamic Study

In this work, molecular dynamics simulations were applied to investigate the influence of lipid composition of the model membrane on the insertion of glyphosate (in its charged state, GLYP2-). The profiles of free energy, entropy and enthalpy were obtained through umbrella sampling calculations, for lipid bilayers composed by only 1,2-dipalmitoyl-sn-glycerol-3-phosphocholine (DPPC), only 1,2-dipalmitoyl-sn-glycerol-3-phosphoserine (DPPS) or a symmetric binary mixture of DPPC and DPPS. In general, the location, the values of minima and maxima of the free energy, and the trend of free energy profiles are influenced by the lipid composition of the lipid bilayer. The driving force in the glyphosate insertion process depends on the lipid composition of the membrane model. If the lipid bilayer is composed solely of DPPS or DPPC, GLYP2- insertion is driven by a favorable enthalpic change. However, if the membrane is composed of a mixture of both lipids, this process is driven by a favorable entropic change. In the lipid bilayer containing DPPS, the glyphosate was found to penetrate hydrated and coordinated with Na+ ions, in contrast to the pure zwitterionic lipid bilayer which penetrated only hydrated. This effect is independent of the concentration of sodium ions present in the bulk solution.

Antibacterial Effect of Chitosan-Gold Nanoparticles and Computational Modeling of the Interaction between Chitosan and a Lipid Bilayer Model

Pathogenic bacteria have the ability to develop antibiotic resistance mechanisms. Their action consists mainly in the production of bacterial enzymes that inactivate antibiotics or the appearance of modifications that prevent the arrival of the drug at the target point or the alteration of the target point itself, becoming a growing problem for health systems. Chitosan-gold nanoparticles (Cs-AuNPs) have been shown as effective bactericidal materials avoiding damage to human cells. In this work, Cs-AuNPs were synthesized using chitosan as the reducing agent, and a systematic analysis of the influence of the synthesis parameters on the size and zeta potential of the Cs-AuNPs and their UV-vis spectra was carried out. We used a simulation model to characterize the interaction of chitosan with bacterial membranes, using a symmetric charged bilayer and two different chitosan models with different degrees of the chitosan amine protonation as a function of pH, with the aim to elucidate the antibacterial mechanism involving the cell wall disruption. The Cs-AuNP antibacterial activity was evaluated to check the simulation model.

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.

Molecular dynamics simulations of glyphosate in a DPPC lipid bilayer

Extensive molecular dynamics simulations have been performed to study the effect of glyphosate (in their neutral and charged forms, GLYP and GLYP2-, respectively) on fully hydrated DiPalmitoylPhosphatidylCholine (DPPC) lipid bilayer. First, we calculated the free energy profile (using the Umbrella Sampling technique) for both states of charge of glyphosate. The minimum value for the free energy for GLYP is ∼-60 kJ mol-1 located at z = ±1.7 nm (from the lipid bilayer center), and there is almost no maximum at the center of the lipid bilayer. By contrast, the minimum for GLYP2- is ∼-35 kJ mol-1 located at z = ± 1.4 nm (from the lipid bilayer center), and the maximum reaches ∼35 kJ mol-1 at the center of the lipid bilayer. Then, different lipid bilayer properties were analyzed for different glyphosate:lipid (G:L) ratios. The mean area per lipid was slightly affected, increasing only 5% (in the presence of glyphosate at high concentrations), which is in agreement with the slight decrease in deuterium order parameters. As for the thickness of the bilayer, it is observed that the state of charge produces opposite effects. On one hand, the neutral state produces an increase in the thickness of the lipid bilayer; on the other, the charged form produces a decrease in the thickness, which not depend linearly on the G:L ratios, either. The orientation of the DPPC head groups is practically unaffected throughout the range of the G:L ratios studied. Finally, the mobility of the lipids of the bilayer is strongly affected by the presence of glyphosate, considerably increasing its lateral diffusion coefficient noteworthy (one order of magnitude), with increasing G:L ratio.

Protonation state of inhibitors determines interaction sites within voltage-gated sodium channels

Voltage-gated sodium channels are essential for carrying electrical signals throughout the body, and mutations in these proteins are responsible for a variety of disorders, including epilepsy and pain syndromes. As such, they are the target of a number of drugs used for reducing pain or combatting arrhythmias and seizures. However, these drugs affect all sodium channel subtypes found in the body. Designing compounds to target select sodium channel subtypes will provide a new therapeutic pathway and would maximize treatment efficacy while minimizing side effects. Here, we examine the binding preferences of nine compounds known to be sodium channel pore blockers in molecular dynamics simulations. We use the approach of replica exchange solute tempering (REST) to gain a more complete understanding of the inhibitors' behavior inside the pore of NavMs, a bacterial sodium channel, and NavPas, a eukaryotic sodium channel. Using these simulations, we are able to show that both charged and neutral compounds partition into the bilayer, but neutral forms more readily cross it. We show that there are two possible binding sites for the compounds: (i) a site on helix 6, which has been previously determined by many experimental and computational studies, and (ii) an additional site, occupied by protonated compounds in which the positively charged part of the drug is attracted into the selectivity filter. Distinguishing distinct binding poses for neutral and charged compounds is essential for understanding the nature of pore block and will aid the design of subtype-selective sodium channel inhibitors.

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.

Role of charged lipids in membrane structures - Insight given by simulations

Lipids and proteins are the main components of cell membranes. It is becoming increasingly clear that lipids, in addition to providing an environment for proteins to work in, are in many cases also able to modulate the structure and function of those proteins. Particularly charged lipids such as phosphatidylinositols and phosphatidylserines are involved in several examples of such effects. Molecular dynamics simulations have proved an invaluable tool in exploring these aspects. This so-called computational microscope can provide both complementing explanations for the experimental results and guide experiments to fruitful directions. In this paper, we review studies that have utilized molecular dynamics simulations to unravel the roles of charged lipids in membrane structures. We focus on lipids as active constituents of the membranes, affecting both general membrane properties as well as non-lipid membrane components, mainly proteins. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Composition effect on the aggregate/solution interface of a nematic lyotropic liquid crystal

We study how the concentration of phospholipids, decanol and NaCl affects the bilayer and the liquid crystal/solution interface structures of lyotropic nematic liquid crystals.

From dioxin to dioxin congeners: understanding the differences in hydrophobic aggregation in water and absorption into lipid membranes by means of atomistic simulations

Molecular dynamics simulations indicate that dioxin likely accumulates in cell membranes more than its congeners (anthracene, tetrahydrodioxin, 3,3′,5,5′-tetrachlorobiphenyl, and 1,2-dihydroxytetrahydrodibenzo-p-dioxin).

Hydrophobic aggregation and collective absorption of dioxin into lipid membranes: insights from atomistic simulations

Dioxins are a highly toxic class of chlorinated aromatic chemicals. They have been extensively studied, but several molecular-level details of their action are still missing. Here we present molecular dynamics simulations of their absorption and diffusion through cell membranes. We show that, due to their hydrophobic character, dioxins can quickly penetrate into a lipid membrane, both as single molecules and as aggregates. We find clear evidence for their ability to accumulate in cell membranes. Our free energy calculations indicate that subsequent transport into the cell is unlikely to be a simple diffusive process.

Volatile anesthetics inhibit sodium channels without altering bulk lipid bilayer properties

Although general anesthetics are clinically important and widely used, their molecular mechanisms of action remain poorly understood. Volatile anesthetics such as isoflurane (ISO) are thought to alter neuronal function by depressing excitatory and facilitating inhibitory neurotransmission through direct interactions with specific protein targets, including voltage-gated sodium channels (Na(v)). Many anesthetics alter lipid bilayer properties, suggesting that ion channel function might also be altered indirectly through effects on the lipid bilayer. We compared the effects of ISO and of a series of fluorobenzene (FB) model volatile anesthetics on Na(v) function and lipid bilayer properties. We examined the effects of these agents on Na(v) in neuronal cells using whole-cell electrophysiology, and on lipid bilayer properties using a gramicidin-based fluorescence assay, which is a functional assay for detecting changes in lipid bilayer properties sensed by a bilayer-spanning ion channel. At clinically relevant concentrations (defined by the minimum alveolar concentration), both the FBs and ISO produced prepulse-dependent inhibition of Na(v) and shifted the voltage dependence of inactivation toward more hyperpolarized potentials without affecting lipid bilayer properties, as sensed by gramicidin channels. Only at supra-anesthetic (toxic) concentrations did ISO alter lipid bilayer properties. These results suggest that clinically relevant concentrations of volatile anesthetics alter Na(v) function through direct interactions with the channel protein with little, if any, contribution from changes in bulk lipid bilayer properties. Our findings further suggest that changes in lipid bilayer properties are not involved in clinical anesthesia.

Molecular dynamics simulation of the partitioning of benzocaine and phenytoin into a lipid bilayer

Molecular dynamics simulations were used to examine the partitioning behaviour of the local anaesthetic benzocaine and the anti-epileptic phenytoin into lipid bilayers, a factor that is critical to their mode of action. Free energy methods are used to quantify the thermodynamics of drug movement between water and octanol as well as for permeation across a POPC membrane. Both drugs are shown to favourably partition into the lipid bilayer from water and are likely to accumulate just inside the lipid headgroups where they may alter bilayer properties or interact with target proteins. Phenytoin experiences a large barrier to cross the centre of the bilayer due to less favourable energetic interactions in this less dense region of the bilayer. Remarkably, in our simulations both drugs are able to pull water into the bilayer, creating water chains that extend back to bulk, and which may modify the local bilayer properties. We find that the choice of atomic partial charges can have a significant impact on the quantitative results, meaning that careful validation of parameters for new drugs, such as performed here, should be performed prior to their use in biomolecular simulations.