Computer simulation of small molecule permeation across a lipid bilayer: dependence on bilayer properties and solute volume, size, and cross-sectional area

Cell membrane permeation is required for most drugs to reach their biological target, and understanding this process is therefore crucial for rational drug design. Recent molecular dynamics simulations have studied the permeation of eight small molecules through a phospholipid bilayer. Unlike experiments, atomistic simulations allow the direct calculation of diffusion and partition coefficients of solutes at different depths inside a lipid membrane. Further analyses of the simulations suggest that solute diffusion is less size-dependent and solute partitioning more size-dependent than was commonly thought.     

Permeability of lipid bilayers to water and ionic solutes

The lipid bilayer moiety of biological membranes is considered to be the primary barrier to free diffusion of water and solutes. This conclusion arises from observations of lipid bilayer model membrane systems, which are generally less permeable than biological membranes. However, the nature of the permeability barrier remains unclear, particularly with respect to ionic solutes. For instance, anion permeability is significantly greater than cation permeability, and permeability to proton-hydroxide is orders of magnitude greater than other monovalent inorganic ions. In this review, we first consider bilayer permeability to water and discuss proposed permeation mechanisms which involve transient defects arising from thermal fluctuations. We next consider whether such defects can account for ion permeation, including proton-hydroxide flux. We conclude that at least two varieties of transient defects are required to explain permeation of water and ionic solutes.     

Does CO2 permeate through aquaporin-1?

Aquaporins facilitate water permeation across biological membranes. Additionally, glycerol and other small neutral solutes are permeated by related aquaglyceroporins. The role of aquaporins in gas permeation has been a long-standing and controversially discussed issue. We present an extensive set of atomistic molecular dynamics simulations that address the question of CO(2) permeation through human aquaporin-1. Free energy profiles derived from the simulations display a barrier of approximately 23 kJ/mol in the aromatic/arginine constriction region of the water pore, whereas a barrier of approximately 4 kJ/mol was observed for a palmitoyloleoylphosphatidylethanolamine lipid bilayer membrane. The results indicate that significant aquaporin-1-mediated CO(2) permeation is to be expected only in membranes with a low intrinsic CO(2) permeability.     

膜上相互作用对平板双分子层脂膜电性质的影响

以平板双分子层脂膜作为生物膜的简单模型,建立用平板双分子层脂膜电性质研究药物－生物膜相互作用的方法。研究以具有典型特征的物质－表面活性剂、自由基、金属手性配合物与平板膜的相互作用引起膜电性质的规律性改变;重组人Ｂ型血红细胞膜与溶液中抗Ｂ单克隆抗体发生特异相互作用时,膜电阻快速下降,下降的速率与加入的抗体量成正相关。在研究发生在平板膜上的典型反应的基础上,通过对膜电性质的监测和分析,从而确认平板双分     

Structure and dynamics of sphingomyelin bilayer: insight gained through systematic comparison to phosphatidylcholine

Sphingomyelin, one of the main lipid components of biological membranes, is actively involved in various cellular processes such as protein trafficking and signal transduction. In particular, specific lateral domains enriched in sphingomyelin and cholesterol have been proposed to play an important functional role in biomembranes, although their precise characteristics have remained unclear. A thorough understanding of the functional role of membranes requires detailed knowledge of their individual lipid components. Here, we employ molecular dynamics simulations to conduct a systematic comparison of a palmitoylsphingomyelin (PSM, 16:0-SM) bilayer with a membrane that comprises dipalmitoylphosphatidylcholine (DPPC) above the main phase transition temperature. We clarify atomic-scale properties that are specific to sphingomyelin due to its sphingosine moiety, and further discuss their implications for SM-rich membranes. We find that PSM bilayers, and in particular the dynamics of PSM systems, are distinctly different from those of a DPPC bilayer. When compared with DPPC, the strong hydrogen bonding properties characteristic to PSM are observed to lead to considerable structural changes in the polar headgroup and interface regions. The strong ordering of PSM acyl chains and specific ordering effects in the vicinity of a PSM-water interface reflect this issue clearly. The sphingosine moiety and related hydrogen bonding further play a crucial role in the dynamics of PSM bilayers, as most dynamic properties, such as lateral and rotational diffusion, are strongly suppressed. This is most evident in the rotational motion characterized by spin-lattice relaxation times and the decay of hydrogen bond autocorrelation functions that are expected to be important in complexation of SM with other lipids in many-component bilayers. A thorough understanding of SM bilayers would greatly benefit from nuclear magnetic resonance experiments for acyl chain ordering and dynamics, allowing full comparison of these simulations to experiments.     

Atomistic Simulations of Pore Formation and Closure in Lipid Bilayers

Cellular membranes separate distinct aqueous compartments, but can be breached by transient hydrophilic pores. A large energetic cost prevents pore formation, which is largely dependent on the composition and structure of the lipid bilayer. The softness of bilayers and the disordered structure of pores make their characterization difficult. We use molecular-dynamics simulations with atomistic detail to study the thermodynamics, kinetics, and mechanism of pore formation and closure in DLPC, DMPC, and DPPC bilayers, with pore formation free energies of 17, 45, and 78 kJ/mol, respectively. By using atomistic computer simulations, we are able to determine not only the free energy for pore formation, but also the enthalpy and entropy, which yields what is believed to be significant new insights in the molecular driving forces behind membrane defects. The free energy cost for pore formation is due to a large unfavorable entropic contribution and a favorable change in enthalpy. Changes in hydrogen bonding patterns occur, with increased lipid-water interactions, and fewer water-water hydrogen bonds, but the total number of overall hydrogen bonds is constant. Equilibrium pore formation is directly observed in the thin DLPC lipid bilayer. Multiple long timescale simulations of pore closure are used to predict pore lifetimes. Our results are important for biological applications, including the activity of antimicrobial peptides and a better understanding of membrane protein folding, and improve our understanding of the fundamental physicochemical nature of membranes.     

Atomistic Simulations of Pore Formation and Closure in Lipid Bilayers

Cellular membranes separate distinct aqueous compartments, but can be breached by transient hydrophilic pores. A large energetic cost prevents pore formation, which is largely dependent on the composition and structure of the lipid bilayer. The softness of bilayers and the disordered structure of pores make their characterization difficult. We use molecular-dynamics simulations with atomistic detail to study the thermodynamics, kinetics, and mechanism of pore formation and closure in DLPC, DMPC, and DPPC bilayers, with pore formation free energies of 17, 45, and 78 kJ/mol, respectively. By using atomistic computer simulations, we are able to determine not only the free energy for pore formation, but also the enthalpy and entropy, which yields what is believed to be significant new insights in the molecular driving forces behind membrane defects. The free energy cost for pore formation is due to a large unfavorable entropic contribution and a favorable change in enthalpy. Changes in hydrogen bonding patterns occur, with increased lipid-water interactions, and fewer water-water hydrogen bonds, but the total number of overall hydrogen bonds is constant. Equilibrium pore formation is directly observed in the thin DLPC lipid bilayer. Multiple long timescale simulations of pore closure are used to predict pore lifetimes. Our results are important for biological applications, including the activity of antimicrobial peptides and a better understanding of membrane protein folding, and improve our understanding of the fundamental physicochemical nature of membranes.     

Combining molecular modeling with experimental methodologies: mechanism of membrane permeation and accumulation of ofloxacin

The interaction between ofloxacin, as a model drug of the fluoroquinolone class, and biomembranes was examined as the possible initial step in a transmembrane diffusion process. Dipalmitoylphosphatidylcholine was used for the preparation of biomembrane models. The influence of environmental conditions and protonation on molecular physicochemical behavior, and hence on the membrane interaction, was investigated by differential scanning calorimetry (DSC). This technique has been shown to be very effective in the interpretation of interactions of drug microspeciations with biomembranes. These findings suggest that the interaction occurred owing to ionic and hydrophobic forces showing how the passage through the membrane is mainly favored in the pH interval 6–7.4. It was demonstrated that a pH gradient through model membranes may be responsible for a poorly homogeneous distribution of ofloxacin (or other related fluoroquinolones), which justifies the in vivo accumulation properties of this drug. DSC experiments, which are in agreement with computational data, also showed that the complexing capability of ofloxacin with regard to Mg⁺⁺ or Ca⁺⁺ may govern the drug entrance into bacterial cells before the DNA Girase inhibition and could ensure the formation of hydrophobic and more fluid phospholipid domains on the surface of the model membrane. These regions are more permeable with regard to various solutes, as well as ofloxacin, allowing a so-called ‘self-promoted entrance pathway’. The combination of experimental methodologies with computational data allowed a further rationalization of the results and opened new perspectives into the mechanism of action of ofloxacin, namely its interaction with lipid bilayers and drug–divalent cation complex formation, which might be extended to the entire fluoroquinolone class. Ofloxacin accumulation within Escherichia coli ATCC 25922 was measured as a function of time. Also in this example, the environmental conditions influenced ofloxacin penetration and accumulation. The in vitro experiments, reported here, show that a suitable balance of hydrophilic and hydrophobic fluoroquinolone properties needs to occur for there to be increased drug permeation.     

Lipid Gymnastics: Evidence of Complete Acyl Chain Reversal in Oxidized Phospholipids from Molecular Simulations

In oxidative environments, biomembranes contain oxidized lipids with short, polar acyl chains. Two stable lipid oxidation products are PoxnoPC and PazePC. PoxnoPC has a carbonyl group, and PazePC has an anionic carboxyl group pendant at the end of the short, oxidized acyl chain. We have used MD simulations to explore the possibility of complete chain reversal in OXPLs in POPC-OXPL mixtures. The polar AZ chain of PazePC undergoes chain reversal without compromising the lipid bilayer integrity at concentrations up to 25% OXPL, and the carboxyl group points into the aqueous phase. Counterintuitively, the perturbation of overall membrane structural and dynamic properties is stronger for PoxnoPC than for PazePC. This is because of the overall condensing and ordering effect of sodium ions bound strongly to the lipids in the PazePC simulations. The reorientation of AZ chain is similar for two different lipid force fields. This work provides the first molecular evidence of the “extended lipid conformation” in phospholipid membranes. The chain reversal of PazePC lipids decorates the membrane interface with reactive, negatively charged functional groups. Such chain reversal is likely to exert a profound influence on the structure and dynamics of biological membranes, and on membrane-associated biological processes.     

Effect of sodium bicarbonate as a pharmaceutical formulation excipient on the interaction of fluvastatin with membrane phospholipids

Excipients in the pharmaceutical formulation of oral drugs are notably employed to improve drug stability. However, they can affect drug absorption and bioavailability. Passive transport through intestinal cell walls is the main absorption mechanism of drugs and, thus, involves an interaction with the membrane lipids. Therefore in this work, the effect of the excipient NaHCO₃ on the interaction of the anticholesterolemic drug fluvastatin sodium (FS) with membrane phospholipids was investigated by ¹H NMR and FTIR spectroscopy. Sodium bicarbonate is often combined with fluvastatin for oral delivery to prevent its degradation. We have used model DMPC/DMPS membranes to mimic the phospholipid content of gut cell membranes. The results presented in this work show a 100% affinity of FS for the membrane phospholipids that is not modified by the presence of the excipient. However, NaHCO₃ is shown to change the interaction mechanism of the drug. According to our data, FS enters the DMPC/DMPS bilayer interface by interacting with the lipids’ polar headgroups and burying its aromatic moieties into the apolar core. Moreover, lipid segregation takes place between the anionic and zwitterionic lipids in the membranes due to a preferential interaction of FS with phosphatidylserines. The excipient counteracts this favored interaction without affecting the drug affinity and location in the bilayer. This work illustrates that preferential interactions with lipids can be involved in passive drug permeation mechanisms and gives evidence of a possible nonpassive role of certain excipients in the interaction of drugs with membrane lipids.     

Uptake and release protocol for assessing membrane binding and permeation by way of isothermal titration calorimetry

The activity of many biomolecules and drugs crucially depends on whether they bind to biological membranes and whether they translocate to the opposite lipid leaflet and trans aqueous compartment. A general strategy to measure membrane binding and permeation is the uptake and release assay, which compares two apparent equilibrium situations established either by the addition or by the extraction of the solute of interest. Only solutes that permeate the membrane sufficiently fast do not show any dependence on the history of sample preparation. This strategy can be pursued for virtually all membrane-binding solutes, using any method suitable for detecting binding. Here, we present in detail one example that is particularly well developed, namely the nonspecific membrane partitioning and flip-flop of small, nonionic solutes as characterized by isothermal titration calorimetry. A complete set of experiments, including all sample preparation procedures, can typically be accomplished within 2 days. Analogous protocols for studying charged solutes, virtually water-insoluble, hydrophobic compounds or specific ligands are also considered.     

Mycobacterial cell wall: Structure and role in natural resistance to antibiotics

Abstract Mycobacteria show a high degree of intrinsic resistance to most antibiotics and chemotherapeutic agents. The low permeability of the mycobacterial cell wall, with its unusual structure, is now known to be a major factor in this resistance. Thus hydrophilic agents cross the cell wall slowly because the myobacterial porin is inefficient in allowing the permeation of solutes and exists in low concentration. Lipophilic agents are presumably slowed down by the lipid bilayer which is of unusually low fluidity and abnormal thickness. Nevertheless, the cell wall barrier alone cannot produce significant levels of drug resistance, which requires synergistic contribution from a second factor, such as the enzymatic inactivation of drugs.     

Electrical oscillation and fluctuation in phospholipid membranes. Phospholipids can form a channel without protein

Fluctuations and/or step-wise changes in membrane potential and electrical current were observed in bilayer membranes of dioleoylphosphatidylcholine (DOPC) in the absence of any channel protein. The DOPC membranes consisted of three types: black lipid membranes, pipette-clamp membranes and lipid membranes transferred to porous filter paper by conventional Langmuir-Blodgett techniques. This finding is significant since phospholipids are the main constituents of biomembranes. Lipid molecules with a cis double bond in their carbon skeleton are suggested to be important in the gating or excitation of biomembranes.     

1-[2-Hydroxy-3-octadecan-1'-oate]propyl-2'',2'',5'',5''-tetramethyl pyrolidine-N-oxyl-3''-carboxylate as a potential spin probe for membrane structure studies

The synthesis of a new minimum steric perturbing proxyl nitroxide, which is a derivative of glycerol and contains a stearic acid moiety, has been carried out. Its localization in model membrane L-alpha-dipalmitoyl phosphatidyl choline (DPPC) was ascertained with the help of ESR, DSC, 1H and 31P NMR techniques. The nitroxide was used for detecting the changes in the phase transition temperature of the model membranes in the presence and absence of drugs. The permeation of the vasodilating drug epinephrine has also been studied using this spin label. The results prove the potential applicability of the new spin probe in the spin labeling of biomembranes.     

Effect of ethylene oxide and propylene oxide block copolymers on the permeability of bilayer lipid membranes to small solutes including doxorubicin

The effects of ethylene oxide and propylene oxide block copolymers (pluronics) on the permeability of several weak acids and bases through bilayer lipid membranes have been studied by the methods of monitoring (1) pH shifts near planar bilayers, (2) doxorubicin fluorescence quenching inside liposomes, and (3) current transients in the presence of hydrophobic anions. It has been shown that pluronics facilitate the permeation of comparatively large molecules (such as 2-n-undecylmalonic acid and doxorubicin) across lipid bilayers, while the permeation of small solutes (such as ammonium and acetic acid) remains unaffected. Pluronics also accelerate the translocation of large hydrophobic anions (tetraphenylborate). The effect of pluronics correlates with the content of propylene oxide units: it is enhanced when the portion of polypropylene oxide block in the copolymer is increased. The action of the pluronic on lipid membrane permeability differs from the effect of the conventional detergent Triton X-100, which does not affect doxorubicin transport if added at concentrations similar to those used for pluronics. It has been proposed that pluronics accelerate the processes of solute diffusion within lipid bilayers (in a structure-dependent manner) rather than influencing the rate of solute adsorption/desorption on the membrane surface. We suppose that the effect of pluronics on doxorubicin permeation across lipid bilayers along with the known effect on the multidrug resistance protein determines its influence on the therapeutic activity of anthracycline drugs.     

Molecular modeling and simulation of Mycobacterium tuberculosis cell wall permeability

The low permeability of the mycobacterial cell wall is thought to contribute to the intrinsic drug resistance of mycobacteria. In this study, the permeability of the Mycobacterium tuberculosis cell wall is studied by computer simulation. Thirteen known drugs with diverse chemical structures were modeled as solutes undergoing transport across a model for the M. tuberculosis cell wall. The properties of the solute-membrane complexes were investigated by means of molecular dynamics simulation, especially the diffusion coefficients of the solute molecules inside the cell wall. The molecular shape of the solute was found to be an important factor for permeation through the M. tuberculosis cell wall. Predominant lateral diffusion within, as opposed to transverse diffusion across, the membrane/cell wall system was observed for some solutes. The extent of lateral diffusion relative to transverse diffusion of a solute within a biological cell membrane may be an important finding with respect to absorption distribution, metabolism, elimination, and toxicity properties of drug candidates. Molecular similarity measures among the solutes were computed, and the results suggest that compounds having high molecular similarity will display similar transport behavior in a common membrane/cell wall environment. In addition, the diffusion coefficients of the solute molecules across the M. tuberculosis cell wall model were compared to those across the monolayers of dipalmitoylphosphatidylethanolamine and dimyristoylphosphatidylcholine, are two common phospholipids in bacterial and animal membranes. The differences among these three groups of diffusion coefficients were observed and analyzed.     

Lipid distribution and transport across cellular membranes

In eukaryotic cells, the membranes of different intracellular organelles have different lipid composition, and various biomembranes show an asymmetric distribution of lipid types across the membrane bilayer. Membrane lipid organization reflects a dynamic equilibrium of lipids moving across the bilayer in both directions. In this review, we summarize data supporting the role of specific membrane proteins in catalyzing transbilayer lipid movement, thereby controlling and regulating the distribution of lipids over the leaflets of biomembranes.     

The permeability of liposomes to nonelectrolytes. I. Activation energies for permeation.

The effect of temperature on the permeability of nonelectrolytes across liposomal membranes above and below their transition temperature has been studied by using an osmotic method. Below their transition temperature, liposomes are osmotically insensitive structures but, on addition of gramicidin A, the water permeability so increased that the permeability of solutes could be studied. The measured activation energies for permeation of a variety of nonelectrolytes has been found to increase when a) there is an increase in the capability of the solutes to form hydrogen bonds in water, b) the cholesterol concentration in the membranes increases and c) the membranes pass from a liquid-crystalline to a solid-crystalline state. The change in the activation energy for permeation per hydrogen bond is about 1.8 kcal/mole for all the different liposome systems investigated; the only solute tested that deviated from this correlation was urea, whose activation energy for permeation across a gramicidin-containing system was much lower than expected from its hydrogen-bonding capacity. This finding suggests that urea is permeating across the gramicidin pore. Although the literature contains only incomplete data relating the activation energies for permeation of nonelectrolytes across biological membranes to their hydrogen-bonding capacity, the available evidence suggests that there is a similar correlation to that found in liposomes. Thus, the average increase in the activation energy per hydrogen bond for permeation across ox red cell membranes (Jacobs, Glassman & Parpart, J. Cell. Comp. Physiol. 7:197, 1935) is 2.2 plus or minus 0.4 kcal/mole, a value that is similar to that obtained in liposomes. However, the activation energies for water and urea are - in such a system - very much lower than expected, suggesting that they, too, are permeating by some parallel route such as an aqueous pore.     

A molecular toolkit for highly insulating tethered bilayer lipid membranes on various substrates

Tethered bilayer lipid membranes (tBLMs) are promising model architectures that mimic the structure and function of natural biomembranes. They provide a fluid, stable, and electrically sealing platform for the study of membrane related processes, specifically, the function of incorporated membrane proteins. This paper presents a generic approach toward the synthesis of functional tBLMs adapted for application to various surfaces. The central element of a tethered membrane consists of a lipid bilayer. Its proximal layer is covalently attached via a spacer unit to a solid support, either gold or silicon oxide. The membranes are characterized optically by using surface plasmon resonance spectroscopy (SPR) or ellipsometry and electrically by using electrochemical impedance spectroscopy (EIS). The bilayer membranes obtained show high electrical barrier properties and can be used to incorporate and study small membrane proteins in a functional form.