Phospholipid surface density determines the partitioning and permeability of acetic acid in DMPC:cholesterol bilayers

Recent Experimental Developments in Studying Passive Membrane Transport of Drug Molecules

The ability to measure the passive membrane permeation of drug-like molecules is of fundamental biological and pharmaceutical importance. Of significance, passive diffusion across the cellular membranes plays an effective role in the delivery of many pharmaceutical agents to intracellular targets. Hence, approaches for quantitative measurement of membrane permeability have been the topics of research for decades, resulting in sophisticated biomimetic systems coupled with advanced techniques. In this review, recent developments in experimental approaches along with theoretical models for quantitative and real-time analysis of membrane transport of drug-like molecules through mimetic and living cell membranes are discussed. The focus is on time-resolved fluorescence-based, surface plasmon resonance, and second-harmonic light scattering approaches. The current understanding of how properties of the membrane and permeant affect the permeation process is discussed.

Tail-Oxidized Cholesterol Enhances Membrane Permeability for Small Solutes

Cholesterol renders mammalian cell membranes more compact by reducing the amount of voids in the membrane structure. Because of this, cholesterol is known to regulate the ability of cell membranes to prevent the permeation of water and water-soluble molecules through the membranes. Meanwhile, it is also known that even seemingly tiny modifications in the chemical structure of cholesterol can lead to notable changes in membrane properties. The question is, how significantly do these small changes in cholesterol structure affect the permeability barrier function of cell membranes? In this work, we applied fluorescence methods as well as atomistic molecular dynamics simulations to characterize changes in lipid membrane permeability induced by cholesterol oxidation. The studied 7β-hydroxycholesterol (7β-OH-chol) and 27-hydroxycholesterol (27-OH-chol) represent two distinct groups of oxysterols, namely, ring- and tail-oxidized cholesterols, respectively. Our previous research showed that the oxidation of the cholesterol tail has only a marginal effect on the structure of a lipid bilayer; however, oxidation was found to disturb membrane dynamics by introducing a mechanism that allows sterol molecules to move rapidly back and forth across the membrane-bobbing. Herein, we show that bobbing of 27-OH-chol accelerates fluorescence quenching of NBD-lipid probes in the inner leaflet of liposomes by dithionite added to the liposomal suspension. Systematic experiments using fluorescence quenching spectroscopy and microscopy led to the conclusion that the presence of 27-OH-chol increases membrane permeability to the dithionite anion. Atomistic molecular dynamics simulations demonstrated that 27-OH-chol also facilitates water transport across the membrane. The results support the view that oxysterol bobbing gives rise to successive perturbations to the hydrophobic core of the membrane, and these perturbations promote the permeation of water and small water-soluble molecules through a lipid bilayer. The observed impairment of permeability can have important consequences for eukaryotic organisms. The effects described for 27-OH-chol were not observed for 7β-OH-chol which represents ring-oxidized sterols.

Predicting the Time of Entry of Nanoparticles in Lipid Membranes

The number of engineered nanoparticles for applications in the biomedical arena has grown tremendously over the last years due to advances in the science of synthesis and characterization. For most applications, the crucial step is the transport through a physiological cellular membrane. However, the behavior of nanoparticles in a biological matrix is a very complex problem that depends not only on the type of nanoparticle but also on its size, shape, phase, surface charge, chemical composition, and agglomeration state. In this paper, we introduce a streamlined theoretical model that predicts the average time of entry of nanoparticles in lipid membranes, using a combination of molecular dynamics simulations and statistical approaches. The model identifies four parameters that separate the contributions of nanoparticle characteristics (i.e., size, shape, solubility) from the membrane properties (density distribution). This factorization allows the inclusion of data obtained from both experimental and computational sources, as well as a rapid estimation of large sets of permutations in membranes. The robustness of the model is supported by experimental data carried out in lipid vesicles encapsulating graphene quantum dots as nanoparticles. Given the high level of interest across multiple areas of study in modulating intracellular targets, and the need to understand and improve the applications of nanoparticles and to assess their effect on human health (i.e., cytotoxicity, bioavailability), this work contributes to the understanding and prediction of interactions between nanoparticles and lipid membranes.

Bobbing of Oxysterols: Molecular Mechanism for Translocation of Tail-Oxidized Sterols through Biological Membranes

Translocation of sterols between cellular membrane leaflets is of key importance in membrane organization, dynamics, and signaling. We present a novel translocation mechanism that differs in a unique manner from the established ones. The bobbing mechanism identified here is demonstrated for tail-oxidized sterols, but is expected to be viable for any molecule containing two polar centers at the opposite sides of the molecule. The mechanism renders translocation across a lipid membrane possible without a change in molecular orientation. For tail-oxidized sterols, the bobbing mechanism provides an exceptionally facile means to translocate these signaling molecules across membrane structures and may thus represent an important pathway in the course of their biological action.

A Universal Correlation Predicts Permeability Coefficients of Fluid- and Gel-Phase Phospholipid and Phospholipid-Cholesterol Bilayers for Arbitrary Solutes

The permeability of gel-phase phospholipids is typically about an order of magnitude lower than that of the same compositions in the fluid phase, yet a quantitative description of the ordering factors leading to this difference has been elusive. The present analysis examines these factors with particular focus on the area per phospholipid chain, Ac, and its relationship to the minimum area per molecule in the crystalline state, A0. It is shown that fluid- and gel-phase phospholipid permeabilities can be reconciled by postulating a minimum area per chain Ac,0 = 17.1 Å(2), substantially less than one would estimate by dividing the accepted value A0 = 40.8 Å(2) by 2. An extended data set of phospholipid and phospholipid-cholesterol bilayer permeability data extending over 9 orders of magnitude is analyzed and correlated according to the developed relationship (N = 85, s = 0.3024, r(2) = 0.9332). Individual permeability values are consequently predicted to within an average deviation of 10(0.3024) or about a factor of 2. The analysis is broadly applicable in the fluid phase but is restricted to gel-phase phospholipid compositions that do not contain cholesterol. Guidance for the latter scenario is provided.

Physical insights into the blood-brain barrier translocation mechanisms

The number of individuals suffering from diseases of the central nervous system (CNS) is growing with an aging population. While candidate drugs for many of these diseases are available, most of these pharmaceutical agents cannot reach the brain rendering most of the drug therapies that target the CNS inefficient. The reason is the blood-brain barrier (BBB), a complex and dynamic interface that controls the influx and efflux of substances through a number of different translocation mechanisms. Here, we present these mechanisms providing, also, the necessary background related to the morphology and various characteristics of the BBB. Moreover, we discuss various numerical and simulation approaches used to study the BBB, and possible future directions based on multi-scale methods. We anticipate that this review will motivate multi-disciplinary research on the BBB aiming at the design of effective drug therapies.

Design of Hydrated Porphyrin-Phospholipid Bilayers with Enhanced Magnetic Resonance Contrast

Computer simulations are used to design more hydrated bilayers, formed from amine-modified porphyrin-phospholipids (PoPs). Experiments confirm that the new constructs give rise to bilayers with greater water content. When chelated with manganese, amine-modified PoPs provide improved contrast for magnetic resonance and are safely used for imaging in vivo.

Cholesterol affects C₆₀ translocation across lipid bilayers

Cholesterol plays an important role in regulating the structural properties of phospholipid membranes and further influences the permeability of molecules and nanoparticles. However, nanoparticles' translocation across phospholipid membranes in the presence of cholesterol on the molecular scale is rarely studied. Here, we performed coarse-grained molecular dynamics simulations to probe the translocation of C60, one of the most popular nanoparticles, across dipalmitoylphosphatidylcholine bilayers with different concentrations of cholesterol molecules (0-50 mol%). The results reveal that the presence of cholesterol molecules induces lower area per lipid, larger bilayer thickness, and more ordered orientation of lipid tails. The higher the concentration of cholesterol molecules, the more significant is the condensing effect of lipid bilayer as just mentioned. Besides, dynamic processes, free energy profiles and permeability coefficients further indicate that the permeability of C60 decreases with increasing cholesterol concentration, which can be explained by the condensation effect and reduced free volume. Our researches provide an explicit description of the impact of cholesterol on C60 translocation across lipid bilayers.

A microscopic multiphase diffusion model of viable epidermis permeability

A microscopic model of passive transverse mass transport of small solutes in the viable epidermal layer of human skin is formulated on the basis of a hexagonal array of cells (i.e., keratinocytes) bounded by 4-nm-thick, anisotropic lipid bilayers and separated by 1-μm layers of extracellular fluid. Gap junctions and tight junctions with adjustable permeabilities are included to modulate the transport of solutes with low membrane permeabilities. Two keratinocyte aspect ratios are considered to represent basal and spinous cells (longer) and granular cells (more flattened). The diffusion problem is solved in a unit cell using a coordinate system conforming to the hexagonal cross section, and an efficient two-dimensional treatment is applied to describe transport in both the cell membranes and intercellular spaces, given their thinness. Results are presented in terms of an effective diffusion coefficient, D¯(epi), and partition coefficient, K¯(epi/w), for a homogenized representation of the microtransport problem. Representative calculations are carried out for three small solutes-water, L-glucose, and hydrocortisone-covering a wide range of membrane permeability. The effective transport parameters and their microscopic interpretation can be employed within the context of existing three-layer models of skin transport to provide more realistic estimates of the epidermal concentrations of topically applied solutes.

Permeability of fluid-phase phospholipid bilayers: assessment and useful correlations for permeability screening and other applications

Permeability data (P(lip/w) ) for liquid crystalline phospholipid bilayers composed of egg lecithin and dimyristoylphosphatidylcholine (DMPC) are analyzed in terms of a mathematical model that accounts for free surface area and chain-ordering effects in the bilayer as well as size and lipophilicity of the permeating species. Free surface area and chain ordering are largely determined by temperature and cholesterol content of the membrane, molecular size is represented by molecular weight, and lipophilicity of the barrier region is represented by the 1,9-decadiene/water partition coefficient, following earlier work by Xiang, Anderson, and coworkers. A correlating variable χ = MW(n) σ/(1 -σ) is used to link the results from different membrane systems, where different values of n are tried, and σ denotes a reduced phospholipid density. The group (1 -σ)/σ is a measure of free surface area, but can also be interpreted in terms of free volume. A single exponential function of χ is developed that is able to correlate 39 observations of P(lip/w) for different compounds in egg lecithin at low density, and 22 observations for acetic acid in DMPC at higher densities, spanning nine orders of magnitude to within an rms error for log 10 P(lip/w) of 0.20. The best fit found for n = 0.87 ultimately makes χ much closer to the ratio of molecular to free volumes than surface areas. The results serve as a starting point for estimating passive permeability of cell membranes to nonionized solutes as a function of temperature and cholesterol content of the membrane.

A correlation for 1,9-decadiene/water partition coefficients

An important series of papers by Xiang, Anderson, and coworkers has established the strong correlation between phospholipid bilayer membrane permeability and the 1,9-decadiene/water partition coefficient over a wide range of compounds, elevating the importance of K(decadiene/w) as a predictor of molecular bioavailability. On the basis of a 58-point dataset developed by these authors, this research note develops an optimal correlation predicting log(10) K(decadiene/w) in terms of the octanol/water partition coefficient and four of the Abraham solvation parameters, namely A (hydrogen bond acidity), S (polarity/polarizability), E (excess molar refraction), and V (McGowan characteristic volume). The fitted dataset is described to within a root-mean-square error of 0.42, and the probable error in making a prediction for a compound not present therein is 0.49. It is shown that this correlation error for K(decadiene/w) is the dominant source of uncertainty in applying a comprehensive new model of phospholipid bilayer membrane permeability developed in a companion paper (Nitsche and Kasting, submitted for publication), which superposes the effects of molecular size and lipid density upon the decadiene lipophilicity scale. Thus, more experimental studies to augment the limited existing database on K(decadiene/w) are called for.

Testing physical models of passive membrane permeation

The biophysical basis of passive membrane permeability is well-understood, but most methods for predicting membrane permeability in the context of drug design are based on statistical relationships that indirectly capture the key physical aspects. Here, we investigate molecular mechanics-based models of passive membrane permeability and evaluate their performance against different types of experimental data, including parallel artificial membrane permeability assays (PAMPA), cell-based assays, in vivo measurements, and other in silico predictions. The experimental data sets we use in these tests are diverse, including peptidomimetics, congeneric series, and diverse FDA approved drugs. The physical models are not specifically trained for any of these data sets; rather, input parameters are based on standard molecular mechanics force fields, such as partial charges, and an implicit solvent model. A systematic approach is taken to analyze the contribution from each component in the physics-based permeability model. A primary factor in determining rates of passive membrane permeation is the conformation-dependent free energy of desolvating the molecule, and this measure alone provides good agreement with experimental permeability measurements in many cases. Other factors that improve agreement with experimental data include deionization and estimates of entropy losses of the ligand and the membrane, which lead to size-dependence of the permeation rate.

Large influence of cholesterol on solute partitioning into lipid membranes

Cholesterol plays an important role in maintaining the correct fluidity and rigidity of the plasma membrane of all animal cells, and hence, it is present in concentrations ranging from 20 to 50 mol %. Whereas the effect of cholesterol on such mechanical properties has been studied exhaustively over the last decades, the structural basis for cholesterol effects on membrane permeability is still unclear. Here we apply systematic molecular dynamics simulations to study the partitioning of solutes between water and membranes. We derive potentials of mean force for six different solutes permeating across 20 different lipid membranes containing one out of four types of phospholipids plus a cholesterol content varying from 0 to 50 mol %. Surprisingly, cholesterol decreases solute partitioning into the lipid tail region of the membranes much more strongly than expected from experiments on macroscopic membranes, suggesting that a laterally inhomogeneous cholesterol concentration and permeability may be required to explain experimental findings. The simulations indicate that the cost of breaking van der Waals interactions between the lipid tails of cholesterol-containing membranes account for the reduced partitioning rather than the surface area per phospholipid, which has been frequently suggested as a determinant for solute partitioning. The simulations further show that the partitioning is more sensitive to cholesterol (i) for larger solutes, (ii) in membranes with saturated as compared to membranes with unsaturated lipid tails, and (iii) in membranes with smaller lipid head groups.

Dermal clearance model for epidermal bioavailability calculations

A computational model for estimating dermal clearance in humans of arbitrary, nonmetabolized solutes is presented. The blood capillary component employs slit theory with contributions from both small (10 nm) and large (50 nm) slits. The lymphatic component is derived from previously reported clearance measurements of dermal and subcutaneous injections of (131)I-albumin in humans. Model parameters were fitted to both blood capillary permeability data and lymphatic clearance data. Small molecules are cleared largely by the blood and large molecules by the lymph. The combined model shows a crossover behavior at approximately 29 kDa, in acceptable agreement with the reported value of 16 kDa. When combined with existing models for stratum corneum permeability and appropriate measures of tissue binding, the developed model has the potential to significantly improve tissue concentration estimates for large or highly protein-bound permeants following dermal exposure.

Optimizing PK properties of cyclic peptides: the effect of side chain substitutions on permeability and clearance()

A series of cyclic peptides were designed and prepared to investigate the physicochemical properties that affect oral bioavailabilty of this chemotype in rats. In particular, the ionization state of the peptide was examined by the incorporation of naturally occurring amino acid residues that are charged in differing regions of the gut. In addition, data was generated in a variety of in vitro assays and the usefulness of this data in predicting the subsequent oral bioavailability observed in the rat is discussed.

Human skin permeation of 3-O-alkyl carbamate prodrugs of naltrexone

N-Monoalkyl and N,N-dialkyl carbamate prodrugs of naltrexone (NTX), an opioid antagonist, were synthesized and their in vitro permeation across human skin was determined. Relevant physicochemical properties were also determined. Most prodrugs exhibited lower melting points, lower aqueous solubilities, and higher oil solubilities than NTX. The flux values from N-monoalkyl carbamate prodrugs were significantly higher than those from NTX and N,N-dialkyl carbamates. The melting points of N-monoalkyl carbamate prodrugs were quite low compared to the N,N-dialkyl carbamate prodrugs and NTX. Heats of fusion for the N,N-dialkyl carbamate prodrugs were higher than that for NTX. N-Monoalkyl carbamate prodrugs had higher stratum corneum/vehicle partition coefficients than their N,N-dialkyl counterparts. Higher percent prodrug bioconversion to NTX in skin appeared to be related to increased skin flux. N,N-Dialkyl carbamate prodrugs were more stable in buffer and in plasma than N-monoalkyl carbamate prodrugs. In conclusion, N-monoalkyl carbamate prodrugs of NTX improved the systemic delivery of NTX across human skin in vitro. N,N-Dialkyl substitution in the prodrug moiety decreased skin permeation and plasma hydrolysis to the parent drug. The cross-sectional area of the carbamate head group was the major determinant of flux of the N-monoalkyl and N,N-dialkyl carbamate prodrugs of NTX.

Alternative measures of lipophilicity: from octanol-water partitioning to IAM retention

This review describes lipophilicity parameters currently used in drug design and QSAR studies. After a short historical overview, the complex nature of lipophilicity as the outcome of polar/nonpolar inter- and intramolecular interactions is analysed and considered as the background for the discussion of the different lipophilicity descriptors. The first part focuses on octanol-water partitioning of neutral and ionisable compounds, evaluates the efficiency of predictions and provides a short description of the experimental methods for the determination of distribution coefficients. A next part is dedicated to reversed-phase chromatographic techniques, HPLC and TLC in lipophilicity assessment. The two methods are evaluated for their efficiency to simulate octanol-water and the progress achieved in the refinement of suitable chromatographic conditions, in particular in the field of HPLC, is outlined. Liposomes as direct models of biological membranes are examined and phospolipophilicity is compared to the traditional lipophilicity concept. Difficulties associated with liposome-water partitioning are discussed. The last part focuses on Immobilised Artificial Membrane (IAM) chromatography as an alternative which combines membrane simulation with rapid measurements. IAM chromatographic retention is compared to octanol-water and liposome-water partitioning as well as to reversed-phase retention and its potential to predict biopartitioning and biological activities is discussed.

Quantitative visualization of passive transport across bilayer lipid membranes

The ability to predict and interpret membrane permeation coefficients is of critical importance, particularly because passive transport is crucial for the effective delivery of many pharmaceutical agents to intracellular targets. We present a method for the quantitative measurement of the permeation coefficients of protonophores by using laser confocal scanning microscopy coupled to microelectrochemistry, which is amenable to precise modeling with the finite element method. The technique delivers well defined and high mass transport rates and allows rapid visualization of the entire pH distribution on both the cis and trans side of model bilayer lipid membranes (BLMs). A homologous series of carboxylic acids was investigated as probe molecules for BLMs composed of soybean phosphatidylcholine. Significantly, the permeation coefficient decreased with acyl tail length contrary to previous work and to Overton's rule. The reasons for this difference are considered, and we suggest that the applicability of Overton's rule requires re-evaluation.

Structural determinants of water permeability through the lipid membrane

Despite intense study over many years, the mechanisms by which water and small nonelectrolytes cross lipid bilayers remain unclear. While prior studies of permeability through membranes have focused on solute characteristics, such as size, polarity, and partition coefficient in hydrophobic solvent, we focus here on water permeability in seven single component bilayers composed of different lipids, five with phosphatidylcholine headgroups and different chain lengths and unsaturation, one with a phosphatidylserine headgroup, and one with a phosphatidylethanolamine headgroup. We find that water permeability correlates most strongly with the area/lipid and is poorly correlated with bilayer thickness and other previously determined structural and mechanical properties of these single component bilayers. These results suggest a new model for permeability that is developed in the accompanying theoretical paper in which the area occupied by the lipid is the major determinant and the hydrocarbon thickness is a secondary determinant. Cholesterol was also incorporated into DOPC bilayers and X-ray diffuse scattering was used to determine quantitative structure with the result that the area occupied by DOPC in the membrane decreases while bilayer thickness increases in a correlated way because lipid volume does not change. The water permeability decreases with added cholesterol and it correlates in a different way from pure lipids with area per lipid, bilayer thickness, and also with area compressibility.

Theory of passive permeability through lipid bilayers

Recently measured water permeability through bilayers of different lipids is most strongly correlated with the area per lipid A rather than with other structural quantities such as the thickness. This paper presents a simple three-layer theory that incorporates the area dependence in a physically realistic way and also includes the thickness as a secondary modulating parameter. The theory also includes the well-known strong correlation of permeability upon the partition coefficients of general solutes in hydrocarbon environments (Overton's rule). Two mathematical treatments of the theory are given; one model uses discrete chemical kinetics and one model uses the Nernst-Planck continuum equation. The theory is fit to the recent experiments on water permeability in the accompanying paper.

Interaction of piroxicam and meloxicam with DMPG/DMPC mixed vesicles: Anomalous partitioning behavior

Small unilamellar vesicles (SUVs) formed from a mixture of dimyristoylphosphatidylcholine (zwitterionic lipid with bulkier headgroup) and dimyristoylphosphatidylglycerol (anionic lipid with relatively smaller headgroup) allows better modulation of the physical properties of lipid bilayers compared to SUVs formed by a single type of lipid, providing us with a better model system to study the effect of membrane parameters on the partitioning of small molecules. Membrane parameter like packing of the vesicles is more pronounced in the gel phase and hence the study was carried out in the gel phase. Mixed vesicles formed from DMPG and DMPC with the mole percent ratio of 100:0, 90:10 and 80:20 were used for this study. As examples of polar solutes, piroxicam and meloxicam, two Non Steroidal Anti-inflammatory Drugs (NSAIDs) were chosen. The pH was adjusted to 2.8 in order to eliminate the presence of anionic forms of the drugs that would not approach the vesicles containing negatively charged DMPG (50% deprotonated at pH 2.8). Surface potential measured by using TNS (2,6-p-toluidinonaphthalene sulfonate, sodium salt) as surface charge sensitive probe showed no significant changes in the surface electrostatics in increasing DMPC content from 0 to 20%. Transmission electron microscopy (TEM) was used to characterize SUVs of different composition at pH 2.8. The average diameter of the mixed vesicles was found to be smaller than that formed by DMPG and DMPC alone. Partition coefficient (K(P)) of piroxicam and meloxicam was measured using intrinsic fluorescence of these molecules. K(P) value of piroxicam decreases with increase in DMPC content whereas it increases with DMPC content in case of meloxicam. This anomalous behavior of partitioning is unexpected since there was no significant change in surface pH of the vesicles and has been explained in terms of lipid packing and water penetration in the lipid bilayer.

Liposomal drug transport: a molecular perspective from molecular dynamics simulations in lipid bilayers

Computational methods to predict drug permeability across biomembranes prior to synthesis are increasingly desirable to minimize the investment in drug design and development. Significant progress in molecular dynamics (MD) simulation methodologies applied to lipid bilayer membranes, for example, is making it possible to move beyond characterization of the membranes themselves to explore various thermodynamic and kinetic processes governing membrane binding and transport. Such methods are also likely to be directly applicable to the design and optimization of liposomal delivery systems. MD simulations are particularly valuable in addressing issues that are difficult to explore in laboratory experiments due to the heterogeneity of lipid bilayer membranes at the molecular level. Insights emerging from MD simulations are contributing to an understanding of which regions within bilayers are most and least favored by solutes at equilibrium as the solute structure is varied, local diffusivities of permeants, and the origin of the amplified selectivity to permeant size imposed by lipid bilayer membranes, particularly as changes in composition increase acyl chain ordering.

Behaviour of small solutes and large drugs in a lipid bilayer from computer simulations

To reach their biological target, drugs have to cross cell membranes, and understanding passive membrane permeation is therefore crucial for rational drug design. Molecular dynamics simulations offer a powerful way of studying permeation at the single molecule level. Starting from a computer model proven to be able to reproduce the physical properties of a biological membrane, the behaviour of small solutes and large drugs in a lipid bilayer has been studied. Analysis of dihedral angles shows that a few nano seconds are sufficient for the simulations to converge towards common values for those angles, even if the starting structures belong to different conformations. Results clearly show that, despite their difference in size, small solutes and large drugs tend to lie parallel to the bilayer normal and that, when moving from water solution into biomembranes, permeants lose degrees of freedom. This explains the experimental observation that partitioning and permeation are highly affected by entropic effects and are size-dependent. Tilted orientations, however, occur when they make possible the formation of hydrogen bonds. This helps to understand the reason why hydrogen bonding possibilities are an important parameter in cruder approaches which predict drug absorption after administration. Interestingly, hydration is found to occur even in the membrane core, which is usually considered an almost hydrophobic region. Simulations suggest the possibility for highly polar compounds like acetic acid to cross biological membranes while hydrated. These simulations prove useful for drug design in rationalising experimental observations and predicting solute behaviour in biomembranes.

Fatty acid flip-flop and proton transport determined by short-circuit current in planar bilayers

The effect of palmitic acid (PA) and oleic acid (OA) on electrical parameters of planar membranes was studied. We found a substantial difference between the effects of PA and OA on proton transfer. PA induced a small increase in conductance, requiring a new technique for estimating proton-mediated currents across low-conductance planar bilayers in which an electrometer is used to measure the transmembrane current under virtual short circuit (SCC). Open-circuit voltage and SCC were used to determine proton and leak conductances. OA caused a marked increase in membrane conductance, allowing the use of a voltage-clamp technique. From SCC data, we were able to estimate the flip-flop rate constants for palmitate (1 x 10(-6) s(-1)) and oleate (49 x 10(-6) s(-1)) anions. Cholesterol, included in the membrane-forming solution, decreased importantly the leak conductance both in membranes unmodified by FA and in membranes modified by PA added to the bath.

Effects of palmitic acid and cholesterol on proton transport across black lipid membranes

We studied the effect of palmitic acid (PA) and cholesterol (approximately 17 wt.%) on proton translocation across asolectin (charged) and diphytanoylphosphatidylcholine (DPhPC, neutral) black lipid membranes (BLMs). Potential difference (PD), short circuit current (SCC), and conductance (G(total)) were measured with a digital electrometer. Membranes were exposed to pH gradients (0.4-2.0 units), followed by PA addition to bath (symmetrically, 40-65 microM). The membrane conductive pathway was subdivided into an unspecific and a proton-related routes. A computer program estimated the conductances (G(un) and G(H)) of the two pathways from the measured parameters. No significant differences in proton selectivity were found between DPhPC membranes and DPhPC/cholesterol membranes. By contrast, cholesterol incorporation into asolectin increases membranes selectivity to proton. Cholesterol dramatically reduced G(un) reflecting, probably, its ability of inducing order in lipid chains. In asolectin membranes, PA increases proton selectivity, probably by acting as a proton shuttle according to the model proposed by Kamp and et al. [Biochemistry 34 (1995) 11928]. Cholesterol incorporation into asolectin membranes eliminates the PA-induced increase in proton selectivity. In DPhPC and DPhPC/cholesterol membranes, PA does not affect proton selectivity. These results are discussed in terms of the presence of cardiolipin (CL) in asolectin, cholesterol/PA interactions, and cholesterol order-inducing effects on acyl-chains.

Correlation of membrane order and dynamics derived from time-resolved fluorescence measurements with solute permeability

The relevance of order and dynamics of phospholipid bilayer membranes as detected with fluorescent probe molecules embedded in the membranes for describing their permeability properties was studied. Order parameters (S) and rotational diffusion coefficients (Dperpendicular) of 1,6-diphenyl-1,3,5-hexatriene (DPH) and 1-(4-trimethylammoniumphenyl)-6-phenyl-1,3,5-hexatriene (TMA-DPH) in unilamellar vesicles were determined by time-resolved fluorescence spectroscopy. Vesicles consisting of combinations of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine, 1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine, egg sphingomyelin and cholesterol were studied at 288, 298, and 308 degrees K. Permeability coefficients (P) of the model permeant D-[14C]mannitol were determined. A model was proposed for correlating P with both S and Dperpendicular, where S is linked to the average free surface area per lipid molecule and Dperpendicular reflects lipid thermal motion and, thus, redistribution rate of free surface area of the bilayer. P values ranging from 0.9 to 12.4 x 10(-11) cm/s were well described by the model. This supports the notion that permeation depends on membrane structural and dynamic properties. While changes in both S and Dperpendicular, at relative significance varying with the situation, appeared responsible for the effect of lipid composition on permeability, the effect of temperature on P was related primarily to Dperpendicular. P correlated better with S and Dperpendicular obtained with TMA-DPH rather than DPH. The location of the fluorescent probe molecules within the membranes is discussed as the cause for this difference.

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.

Cholesterol effects on the phospholipid condensation and packing in the bilayer: a molecular simulation study

A 15-ns molecular dynamics simulation of the fully hydrated liquid-crystalline dimyristoylphosphatidylcholine-cholesterol (DMPC-Chol) bilayer containing approximately 22 mol% Chol was carried out. The generated trajectory was analysed to investigate the mechanism of the Chol condensing effect on DMPC hydrocarbon chains and the influence of Chol on the chain packing in the membrane. Chol was found to induce stronger van der Waals interactions among the chains, whereas its interactions with the chains were weak. In the DMPC-Chol bilayer, as in the DMPC bilayer, DMPC chains were regularly packed around a chosen chain but around a Chol molecule they were not. DMPC gamma chains made closer contacts with Chol than the beta chains.

Effect of liposome type and membrane fluidity on drug–membrane partitioning analyzed by immobilized liposome chromatography

Immobilized liposome chromatography (ILC) has been proven to be a useful method for the study or rapid screening of drug-membrane interactions. To obtain an adequate liposomal membrane phase for ILC, unilamellar liposomes were immobilized in gel beads by avidin-biotin binding. The retardation of 15 basic drugs on the liposome column could be converted to membrane partitioning coefficients, K(LM). The effects of small or large unilamellar liposomes and multilamellar liposomes on the drug-membrane partitioning were compared. The K(LM) values for both small and large liposomes were similar, but higher than those for the multilamellar liposomes. The basic drugs showed stronger partitioning into negatively charged liposomes than into either neutral liposomes or positively charged liposomes. The membrane fluidity of the immobilized liposomes was modulated by incorporating cholesterol into the liposomal membranes, by changing the acyl chain length and degree of unsaturation of the phospholipids, and by changing the temperature for ILC runs. Our data show that K(LM) obtained using ILC correlated well with those reported by batch studies using free liposomes. It is concluded that negatively charged or cholesterol-containing large unilamellar liposomes are suitable models for the ILC analysis of drug-membrane interactions.

Effect of therapeutic ultrasound on partition and diffusion coefficients in human stratum corneum

Application of various enhancers including ultrasound and chemicals has been shown to enhance transdermal drug transport. Most of these enhancers increase transdermal transport by increasing either partition or diffusion coefficients in lipid bilayers. Although the effect of such enhancers on skin permeability has been measured in many cases, the effect of the same enhancers on solute partition and diffusion coefficients has not been measured since such measurements are difficult. In this paper, we describe application of a new method that can be used to determine the effect of enhancers on partition and diffusion coefficients. This method is based on two independent measurements of the transport properties of the SC and a theoretical model. This method was tested for five solutes including estradiol, napthol, aldosterone, lidocaine, and testosterone. Measurements based on our method showed that the primary effect of therapeutic ultrasound (1 MHz) is on diffusion coefficient and not partition coefficient. Specifically, ultrasound enhanced diffusion coefficients of these solutes by up to 15-fold. However, it did not significantly enhance partition coefficients. The method described in this paper is quite general and can be used to measure the effect of various enhancers on partition and diffusion coefficients.

Effect of cholesterol on molecular transport of organic cations across liposome bilayers probed by second harmonic generation

The effect of cholesterol on the molecular transport of an organic cation, malachite green (MG), across large unilamellar dioleolyphosphatidylglycerol (DOPG) liposome bilayers with 0-50 mol% cholesterol was studied by second harmonic generation (SHG). Because SHG is a surface-specific technique, it requires no labeled molecule, quencher, or shifting agent to distinguish the location of the solute molecules. An additional important feature of SHG is that it is sensitive only to the probe molecules bound to the liposome, whereas other methods can only differentiate between molecules that are outside and those inside the liposome. The transport kinetics of MG across the liposome bilayers was observed in real time, and the results show that cholesterol retards the rate of transport of MG across liposome bilayers. The rate was found to decrease by six times for 50 mol% cholesterol content compared with cholesterol-free liposomes. This demonstrates the applicability of SHG to investigation of the effect of liposome composition on the transport kinetics across the liposome bilayers.

A 1H-NMR study of the permeation of glycolic acid through phospholipid membranes

The transmembrane permeability coefficient of the alpha-hydroxyacid, glycolic acid, has been measured for egg phosphatidylcholine large unilamellar vesicles. The determination of the vesicle concentration independent first-order rate constant for membrane transport and the permeability coefficient were made using an NMR technique employing shift agents. Both the temperature dependence and the dependence on cholesterol content were investigated. The activation energy and the Arrhenius pre-exponential factor were found to be dependent on the cholesterol content. A marked increase in both parameters was observed up to 20 mol% cholesterol, with a further, small increase up to 50%. The pH dependence of permeability was also investigated. An increase in permeability is observed with a decrease in pH, providing a possible explanation for the effectiveness of glycolic acid in skin treatment.

Influence of chain ordering on the selectivity of dipalmitoylphosphatidylcholine bilayer membranes for permeant size and shape

The effects of lipid chain packing and permeant size and shape on permeability across lipid bilayers have been investigated in gel and liquid crystalline dipalmitoylphosphatidylcholine (DPPC) bilayers by a combined NMR line-broadening/dynamic light scattering method using seven short-chain monocarboxylic acids (formic acid, acetic acid, propionic acid, butyric acid, valeric acid, isovaleric acid, and trimethylacetic acid) as permeants. The experimental permeability coefficients are compared with the predictions of a bulk solubility diffusion model in which the bilayer membrane is represented as a slab of bulk hexadecane. Deviations of the observed permeability coefficients (Pm) from the values predicted from solubility diffusion theory (Po) lead to the determination of a correction factor, the permeability decrement f (= Pm/Po), to account for the effects of chain ordering. The natural logarithm of f has been found to correlate linearly with the inverse of the bilayer free surface area with slopes of 25 +/- 2, 36 +/- 3, 45 +/- 8, 32 +/- 12, 33 +/- 4, 49 +/- 12, and 75 +/- 6 A2 for formic acid, acetic acid, propionic acid, butyric acid, valeric acid, isovaleric acid, and trimethylacetic acid, respectively. The slope, which measures the sensitivity of the permeability coefficient of a given permeant to bilayer chain packing, exhibits an excellent linear correlation (r = 0.94) with the minimum cross-sectional area of the permeant and a poor correlation (r = 0.59) with molecular volume, suggesting that in the bilayer interior the permeants prefer to move with their long principal axis along the bilayer normal. Based on these studies, a permeability model combining the effects of bilayer chain packing and permeant size and shape on permeability across lipid membranes is developed.

Avidin-biotin immobilization of unilamellar liposomes in gel beads for chromatographic analysis of drug-membrane partitioning

To construct a homogeneous lipid membrane chromatographic phase, biotinylated unilamellar liposomes of small and large sizes (SUVs and LUVs, respectively) were immobilized in avidin- or streptavidin-derived gel beads in amounts up to 55 micromol phospholipid/ml gel bed at yields above 50%. The immobilized liposomes exhibited excellent stability due to avidin-biotin multiple-site binding. The trapped volume and size distribution of the immobilized liposomes (0.33-0.42 microl/micromol lipid and 20-30 nm diameter for SUVs, 1.7-1.9 microl/micromol lipid and 80-120 nm for LUVs) indicated the unilamellarity and integrity of the immobilized liposomes. Partitioning of 15 pharmaceutical drugs into the bilayers of LUVs immobilized in different gel matrices correlated very well, as shown by chromatographic drug retention analysis. The partitioning of several beta-blockers into the immobilized LUVs showed a close correlation with their partitioning, reported in the literature, into free liposomes. The avidin-biotin-immobilized unilamellar liposomes can thus be used for chromatographic analysis and screening of solute-membrane interactions.

Phase structures of binary lipid bilayers as revealed by permeability of small molecules

The effects of changes in bilayer phase structure on the permeability of acetic acid and trimethylacetic acid were studied in large unilamellar vesicles (LUVs) composed of dipalmitoylphosphatidylcholine (DPPC)/cholesterol (CHOL), dihexadecylphosphatidylcholine (DHPC)/CHOL, or DPPC/dimyristoylphosphatidylcholine (DMPC) using an NMR line-broadening method. Phase transitions were induced by changes in temperature and lipid composition (i.e., XCHOL was varied from 0.0 to 0.5 and XDMPC from 0.0 to 1.0). In DPPC/CHOL and DHPC/CHOL bilayers, the addition of CHOL induces only a modest change in the permeability coefficient (Pm) of acetic acid in the gel-phase (Pbeta') but significantly reduces Pm in ordered and disordered liquid-crystalline phases (Lo and Lalpha). Abrupt changes in slopes in semi-logarithmic plots of Pm vs. XCHOL occur at specific values of XCHOL and temperature corresponding to the boundaries between Pbeta' and Lo or between Lalpha and Lo phases. In most respects, phase diagrams generated from the break points in plots of Pm vs. XCHOL obtained at various temperatures in DHPC/CHOL and DPPC/CHOL bilayers closely resemble those constructed previously for DPPC/CHOL bilayers using NMR and DSC methods. Above Tm, the phase diagrams generated from permeability data reveal the presence of both the disordered (Lalpha) and the ordered (Lo) liquid-crystalline phases, as well as the two-phase coexistence region. In DPPC/DMPC bilayers, the addition of DMPC increases Pm dramatically in the gel phase but only slightly in the liquid-crystalline phase. Abrupt changes in slopes in semi-logarithmic plots of Pm vs. XDMPC also occur at specific values of XDMPC and temperature, from which a phase diagram can be constructed which closely resembles diagrams obtained previously by other methods. These correlations indicate that trans-bilayer permeability measurements can be used to construct lipid bilayer phase diagrams. Positive deviations of Pm from predicted values based on the phase lever rule are observed in the two-phase coexistence regions with the degree of the deviation depending on bilayer chemical composition and temperature. These results may reflect a specific contribution of the interfacial region between two phases to higher solute permeability or may be due to the higher lateral compressibility of lipid bilayers in the two-phase coexistence region.

Association and release of prostaglandin E1 from liposomes

PGE1-lipid interactions were studied in several liposome systems. Data from both circular dichroic (CD) measurements and differential scanning calorimetry (DSC) indicated that PGE1 in the protonated form seeks the less polar environment of the lipid bilayer. CD measurements made on PGE1 in solution showed that the wavelength of maximum absorbance red shifted approximately 8 nm with decreasing solvent polarity. The CD spectrum of liposomal PGE1 prepared in pH 4.5 but not pH 7.2 buffer was also red shifted. There was no red shift in the CD spectrum of PGE1 detected at pH 4.5 in the absence of phospholipid. DSC measurements on DSPC bilayers prepared with 5 mol% PGE1 at pH 4.5 but not pH 7.2 revealed an almost complete loss of the pre-transition as well as broadening of the main phase transition. The amount of 3H-PGE1 initially associated with EPC, POPC or DSPC liposomes was determined using size exclusion filters and centrifugation. This amount was found to be dependent on the pH of the buffer (pH 4.5 >> pH 7.2) and fluidity of the bilayer (EPC = POPC > DSPC), but independent of the lamellarity of the liposome. In all cases, addition of cholesterol reduced the amount of PGE1 associated with the liposome. The time-dependent release of PGE1 from the liposomes was determined by rapidly diluting the sample 100-fold into pH 7.2 buffer. Lipid saturation was a key factor influencing this release. Gel-phase liposomes of DSPC showed a rapid initial release (t(1/2) < 2 min) of PGE1, corresponding to the amount in the outer monolayer, followed by a very slow, almost negligible release of the remaining PGE1. A rapid initial release also occurred in fluid-phase membranes, followed by a more gradual release of the remaining PGE1 over several hours. This release rate could be slowed by increasing the lamellarity of these liposomes, or adding cholesterol to decrease the fluidity of the membrane.

Permeability of acetic acid across gel and liquid-crystalline lipid bilayers conforms to free-surface-area theory

Solubility-diffusion theory, which treats the lipid bilayer membrane as a bulk lipid solvent into which permeants must partition and diffuse across, fails to account for the effects of lipid bilayer chain order on the permeability coefficient of any given permeant. This study addresses the scaling factor that must be applied to predictions from solubility-diffusion theory to correct for chain ordering. The effects of bilayer chemical composition, temperature, and phase structure on the permeability coefficient (Pm) of acetic acid were investigated in large unilamellar vesicles by a combined method of NMR line broadening and dynamic light scattering. Permeability values were obtained in distearoylphosphatidylcholine, dipalmitoylphosphatidylcholine, dimyristoylphosphatidylcholine, and dilauroylphosphatidylcholine bilayers, and their mixtures with cholesterol, at various temperatures both above and below the gel-->liquid-crystalline phase transition temperatures (Tm). A new scaling factor, the permeability decrement f, is introduced to account for the decrease in permeability coefficient from that predicted by solubility-diffusion theory owing to chain ordering in lipid bilayers. Values of f were obtained by division of the observed Pm by the permeability coefficient predicted from a bulk solubility-diffusion model. In liquid-crystalline phases, a strong correlation (r = 0.94) between f and the normalized surface density sigma was obtained: in f = 5.3 - 10.6 sigma. Activation energies (Ea) for the permeability of acetic acid decreased with decreasing phospholipid chain length and correlated with the sensitivity of chain ordering to temperature, [symbol: see text] sigma/[symbol: see text](1/T), as chain length was varied. Pm values decreased abruptly at temperatures below the main phase transition temperatures in pure dipalmitoylphosphatidylcholine and dimyristoylphosphatidylcholine bilayers (30-60-fold) and below the pretransition in dipalmitoylphosphatidylcholine bilayers (8-fold), and the linear relationship between in f and sigma established for liquid-crystalline bilayers was no longer followed. However, in both gel and liquid-crystalline phases in f was found to exhibit an inverse correlation with free surface area (in f = -0.31 - 29.1/af, where af is the average free area (in square angstroms) per lipid molecule). Thus, the lipid bilayer permeability of acetic acid can be predicted from the relevant chain-packing properties in the bilayer (free surface area), regardless of whether chain ordering is varied by changes in temperature, lipid chain length, cholesterol concentration, or bilayer phase structure, provided that temperature effects on permeant dehydration and diffusion and the chain-length effects on bilayer barrier thickness are properly taken into account.

Mechanism of unassisted ion transport across membrane bilayers

To establish how charged species move from water to the nonpolar membrane interior and to determine the energetic and structural effects accompanying this process, we performed molecular dynamics simulations of the transport of Na+ and Cl- across a lipid bilayer located between two water lamellae. The total length of molecular dynamics trajectories generated for each ion was 10 ns. Our simulations demonstrate that permeation of ions into the membrane is accompanied by the formation of deep, asymmetric thinning defects in the bilayer, whereby polar lipid head groups and water penetrate the nonpolar membrane interior. Once the ion crosses the midplane of the bilayer the deformation "switches sides"; the initial defect slowly relaxes, and a defect forms in the outgoing side of the bilayer. As a result, the ion remains well solvated during the process; the total number of oxygen atoms from water and lipid head groups in the first solvation shell remains constant. A similar membrane deformation is formed when the ion is instantaneously inserted into the interior of the bilayer. The formation of defects considerably lowers the free energy barrier to transfer of the ion across the bilayer and, consequently, increases the permeabilities of the membrane to ions, compared to the rigid, planar structure, by approximately 14 orders of magnitude. Our results have implications for drug delivery using liposomes and peptide insertion into membranes.