Partition equilibrium of inhalation anesthetics and alcohols between water and membranes of phospholipids with varying acyl chain-lengths

Ligand binding at the protein-lipid interface: strategic considerations for drug design

Many drug targets are embedded within the phospholipid bilayer of cellular membranes, including G protein-coupled receptors, ion channels, transporters and membrane-bound enzymes. Increasing evidence from biophysical and structural studies suggests that many small-molecule drugs commonly associate with these targets at binding sites at the protein-phospholipid interface. Without a direct path from bulk solvent to a binding site, a drug must first partition in the phospholipid membrane before interacting with the protein target. This membrane access mechanism necessarily affects the interpretation of potency data, structure-activity relationships, pharmacokinetics and physicochemical properties for drugs that target these sites. With an increasing number of small-molecule intramembrane binding sites revealed through X-ray crystallography and cryogenic electron microscopy, we suggest that ligand-lipid interactions likely play a larger role in small-molecule drug action than commonly appreciated. This Perspective introduces key concepts and drug design considerations to aid discovery teams operating within this target space, and discusses challenges and future opportunities in the field.

Prediction of Partition Coefficients of Environmental Toxins Using Computational Chemistry Methods

The partitioning of compounds between aqueous and other phases is important for predicting toxicity. Although thousands of octanol-water partition coefficients have been measured, these represent only a small fraction of the anthropogenic compounds present in the environment. The octanol phase is often taken to be a mimic of the inner parts of phospholipid membranes. However, the core of such membranes is typically more hydrophobic than octanol, and other partition coefficients with other compounds may give complementary information. Although a number of (cheap) empirical methods exist to compute octanol-water (log k _OW) and hexadecane-water (log k _HW) partition coefficients, it would be interesting to know whether physics-based models can predict these crucial values more accurately. Here, we have computed log k _OW and log k _HW for 133 compounds from seven different pollutant categories as well as a control group using the solvation model based on electronic density (SMD) protocol based on Hartree-Fock (HF) or density functional theory (DFT) and the COSMO-RS method. For comparison, XlogP3 (log k _OW) values were retrieved from the PubChem database, and KowWin log k _OW values were determined as well. For 24 of these compounds, log k _OW was computed using potential of mean force (PMF) calculations based on classical molecular dynamics simulations. A comparison of the accuracy of the methods shows that COSMO-RS, KowWin, and XlogP3 all have a root-mean-square deviation (rmsd) from the experimental data of ≈0.4 log units, whereas the SMD protocol has an rmsd of 1.0 log units using HF and 0.9 using DFT. PMF calculations yield the poorest accuracy (rmsd = 1.1 log units). Thirty-six out of 133 calculations are for compounds without known log k _OW, and for these, we provide what we consider a robust prediction, in the sense that there are few outliers, by averaging over the methods. The results supplied may be instrumental when developing new methods in computational ecotoxicity. The log k _HW values are found to be strongly correlated to log k _OW for most compounds.

Modeling the impact of biota on polychlorinated biphenyls (PCBs) fate and transport in Lake Ontario using a population-based multi-compartment fugacity approach

Organisms have long been treated as receptors in exposure studies of polychlorinated biphenyls (PCBs) and other persistent organic pollutants (POPs). The influences of environmental pollution on organisms are well recognized. However, the impact of biota on PCB transport in an environmental system has not been considered in sufficient detail. In this study, a population-based multi-compartment fugacity model is developed by reconfiguring the organisms as populated compartments and reconstructing all the exchange processes between the organism compartments and environmental compartments, especially the previously ignored feedback routes from biota to the environment. We evaluate the model performance by simulating the PCB concentration distribution in Lake Ontario using published loading records. The lake system is divided into three environment compartments (air, water, and sediment) and several organism groups according to the dominant local biotic species. The comparison indicates that the simulated results are well-matched by a list of published field measurements from different years. We identify a new process, called Facilitated Biotic Intermedia Transport (FBIT), to describe the enhanced pollution transport that occurs between environmental media and organisms. As the hydrophobicity of PCB congener increases, the organism population exerts greater influence on PCB mass flows. In a high biomass scenario, the model simulation indicates significant FBIT effects and biotic storage effects with hydrophobic PCB congeners, which also lead to significant shifts in systemic contaminant exchange rates between organisms and the environment.

Chloroform alters interleaflet coupling in lipid bilayers: an entropic mechanism

The interaction of the two leaflets of the plasmatic cell membrane is conjectured to play an important role in many cell processes. Experimental and computational studies have investigated the mechanisms that modulate the interaction between the two membrane leaflets. Here, by means of coarse-grained molecular dynamics simulations, we show that the addition of a small and polar compound such as chloroform alters interleaflet coupling by promoting domain registration. This is interpreted in terms of an entropic gain that would favour frequent chloroform commuting between the two leaflets. The implication of this effect is discussed in relation to the general anaesthetic action.

Separation and enhanced detection of anesthetic compounds using solid phase micro-extraction (SPME)-Raman spectroscopy

Polydimethylsiloxane (PDMS)-based solid-phase micro-extraction (SPME) was used along with Raman spectroscopy (RS) to separate and enhance the detection of five anesthetic compounds (halothane, propofol, isoflurane, enflurane, and etomidate) from aqueous and serum phases. Raman signals in the spectral ranges 250-450 cm(-1) and 950-1050 cm(-1) allowed the unique characterization of all five compounds when extracted into the PDMS phase. The SPME-RS detection of clinically relevant concentrations of aqueous propofol (6.5 μM) and halothane (200 μM) is shown. We quantify the partition coefficient for aqueous halothane in PDMS as log K = 1.9 ± 0.2. Solid-phase micro-extraction of the anesthetics makes their detection possible without the strong autofluorescent interference of serum proteins. Because of low solubility and/or weak Raman scattering, we found it challenging to detect enflurane, isoflurane, and etomidate directly from the aqueous phase, but could we do so with SPME enhancement. These studies show the potential of SPME-RS as a method for the direct detection of anesthetics in blood.

Structure-based prediction of drug distribution across the headgroup and core strata of a phospholipid bilayer using surrogate phases

Solvation of drugs in the core (C) and headgroup (H) strata of phospholipid bilayers affects their physiological transport rates and accumulation. These characteristics, especially a complete drug distribution profile across the bilayer strata, are tedious to obtain experimentally, to the point that even simplified preferred locations are only available for a few dozen compounds. Recently, we showed that the partition coefficient (P) values in the system of hydrated diacetyl phosphatidylcholine (DAcPC) and n-hexadecane (C16), as surrogates of the H- and C-strata of the bilayer composed of the most abundant mammalian phospholipid, PC, agree well with the preferred bilayer location of compounds. High P values are typical for lipophiles accumulating in the core, and low P values are characteristic of cephalophiles preferring the headgroups. This simple pattern does not hold for most compounds, which usually have more even distribution and may also accumulate at the H/C interface. To model complete distribution, the correlates of solvation energies are needed for each drug state in the bilayer: (1) for the H-stratum it is the DAcPC/W P value, calculated as the ratio of the C16/W and C16/DAcPC (W for water) P values; (2) for the C-stratum, the C16/W P value; (3) for the H/C interface, the P values for all plausible molecular poses are characterized using the fragment DAcPC/W and C16/W solvation parameters for the parts of the molecule embedded in the H- and C-strata, respectively. The correlates, each scaled by two Collander coefficients, were used in a nonlinear, mass-balance based model of intrabilayer distribution, which was applied to the easily measurable overall P values of compounds in the DMPC (M = myristoyl) bilayers and monolayers as the dependent variables. The calibrated model for 107 neutral compounds explains 94% of experimental variance, achieves similar cross-validation levels, and agrees well with the nontrivial, experimentally determined bilayer locations for 27 compounds. The resulting structure-based prediction system for intrabilayer distribution will facilitate more realistic modeling of passive transport and drug interactions with those integral membrane proteins, which have the binding sites located in the bilayer, such as some enzymes, influx and efflux transporters, and receptors. If only overall bilayer accumulation is of interest, the 1-octanol/W P values suffice to model the studied set.

Improved spatial resolution of matrix-assisted laser desorption/ionization imaging of lipids in the brain by alkylated derivatives of 2,5-dihydroxybenzoic acid

RATIONALE

Matrix-assisted laser desorption/ionization (MALDI) is one of the major techniques for mass spectrometry imaging (MSI) of biological systems along with secondary-ion mass spectrometry (SIMS) and desorption electrospray mass spectrometry (DESI). The inherent variability of MALDI-MSI signals within intact tissues is related to the heterogeneity of both the sample surface and the matrix crystallization. To circumvent some of these limitations of MALDI-MSI, we have developed improved matrices for lipid analysis based on structural modification of the commonly used matrix 2,5-dihydroxybenzoic acid (DHB).

METHODS

We have synthesized DHB containing -C6H13 and -C12H25 alkyl chains and applied these matrices to rat brain using a capillary sprayer. We utilized a Bruker Ultraflex II MALDI-TOF/TOF mass spectrometer to analyze lipid extracts and tissue sections, and examined these sections with polarized light microscopy and differential interference contrast microscopy.

RESULTS

O-alkylation of DHB yields matrices, which, when applied to brain sections, follow a trend of phase transition from crystals to an oily layer in the sequence DHB → DHB-C6H13 → DHB-C12H25 . MALDI-MSI images acquired with DHB-C12H25 exhibited a considerably higher density of lipids than DHB.

CONCLUSIONS

Comparative experiments with DHB and DHB-C12H25 are presented, which indicate that the latter matrix affords higher lateral resolution than the former.

Disruption of phosphatidylcholine monolayers and bilayers by perfluorobutane sulfonate

Perfluoroalkyl acids (PFAAs) are persistent environmental contaminants resistant to biological and chemical degradation due to the presence of carbon-fluorine bonds. These compounds exhibit developmental toxicity in vitro and in vivo. The mechanisms of toxicity may involve partitioning into lipid bilayers. We investigated the interaction between perfluorobutane sulfonate (PFBS), an emerging PFAA, and model phosphatidylcholine (PC) lipid assemblies (i.e., dimyristoyl-, dipalmitoyl- and distearoylphosphatidylcholine) using fluorescence anisotropy and Langmuir monolayer techniques. PFBS decreased the transition temperature and transition width of PC bilayers. The apparent membrane partition coefficients ranged from 4.9 × 10(2) to 8.2 × 10(2). The effects on each PC were comparable. The limiting molecular area of PC monolayers increased, and the surface pressure at collapse decreased in a concentration-dependent manner. The compressibility of all three PCs was decreased by PFBS. In summary, PFBS disrupted different model lipid assemblies, indicating potential for PFBS to be a human toxicant. However, the effects of PFBS are not as pronounced as those seen with longer chain PFAAs.

Alcohol's effects on lipid bilayer properties

Alcohols are known modulators of lipid bilayer properties. Their biological effects have long been attributed to their bilayer-modifying effects, but alcohols can also alter protein function through direct protein interactions. This raises the question: Do alcohol's biological actions result predominantly from direct protein-alcohol interactions or from general changes in the membrane properties? The efficacy of alcohols of various chain lengths tends to exhibit a so-called cutoff effect (i.e., increasing potency with increased chain length, which that eventually levels off). The cutoff varies depending on the assay, and numerous mechanisms have been proposed such as: limited size of the alcohol-protein interaction site, limited alcohol solubility, and a chain-length-dependent lipid bilayer-alcohol interaction. To address these issues, we determined the bilayer-modifying potency of 27 aliphatic alcohols using a gramicidin-based fluorescence assay. All of the alcohols tested (with chain lengths of 1-16 carbons) alter the bilayer properties, as sensed by a bilayer-spanning channel. The bilayer-modifying potency of the short-chain alcohols scales linearly with their bilayer partitioning; the potency tapers off at higher chain lengths, and eventually changes sign for the longest-chain alcohols, demonstrating an alcohol cutoff effect in a system that has no alcohol-binding pocket.

Model and cell membrane partitioning of perfluorooctanesulfonate is independent of the lipid chain length

Perfluorooctanesulfonic acid (PFOS) is a persistent environmental pollutant that may cause adverse health effects in humans and animals by interacting with and disturbing of the normal properties of biological lipid assemblies. To gain further insights into these interactions, we investigated the effect of PFOS potassium salt on dimyristoyl- (DMPC), dipalmitoyl- (DPPC) and distearoylphosphatidylcholine (DSPC) model membranes using fluorescence anisotropy measurements and differential scanning calorimetry (DSC) and on the cell membrane of HL-60 human leukemia cells and freshly isolated rat alveolar macrophages using fluorescence anisotropy measurements. PFOS produced a concentration-dependent decrease of the main phase transition temperature (T(m)) and an increased peak width (DeltaT(w)) in both the fluorescence anisotropy and the DSC experiments, with a rank order DMPC>DPPC>DSPC. PFOS caused a fluidization of the gel phase of all phosphatidylcholines investigated, but had the opposite effect on the liquid-crystalline phase. The apparent partition coefficients of PFOS between the phosphatidylcholine bilayer and the bulk aqueous phase were largely independent of the phosphatidylcholine chain length and ranged from 4.4x10(4) to 8.8x10(4). PFOS also significantly increased the fluidity of membranes of cells. These findings suggest that PFOS readily partitions into lipid assemblies, independent of their composition, and may cause adverse biological effects by altering their fluidity in a manner that depends on the membrane cooperativity and state (e.g., gel versus liquid-crystalline phase) of the lipid assembly.

Molecular packing in 1-hexanol–DMPC bilayers studied by molecular dynamics simulation

The structure and molecular packing density of a "mismatched" solute, 1-hexanol, in lipid membranes of dimyristoyl phosphatidylcholine (DMPC) was studied by molecular dynamics simulations. We found that the average location and orientation of the hexanol molecules matched earlier experimental data on comparable systems. The local density or molecular packing in DMPC-hexanol was elucidated through the average Voronoi volumes of all heavy (non-hydrogen) atoms. Analogous analysis was conducted on trajectories from simulations of pure 1-hexanol and pure (hydrated) DMPC bilayers. The results suggested a positive volume change, DeltaV(m), of 4 cm(3) mol(-1) hexanol partitioned at 310 K in good accordance with experimental values. Analysis of the apparent volumes of each component in the pure and mixed states further showed that DeltaV(m) reflects a balance between a substantial increase in the packing density of the alcohol upon partitioning and an even stronger loosening in the packing of the lipid. Furthermore, analysis of Voronoi volumes along the membrane normal identifies a distinctive depth dependence of the changes in molecular packing. The outer (interfacial) part of the lipid acyl chains (up to C8) is stretched by about 4%. Concomitantly, the average lateral area per chain decreases and these two effects compensate so that the overall packing density in the outer region, where the hexanol molecules are located, remains practically constant. The core of the bilayer (C9-C13) is slightly thinned. The average lateral area per chain in this region expands, resulting in a looser packing density. The net effect in the core is a 2-3% decrease in density corresponding to a total volume increase of approximately 14 cm(3) mol(-1) hexanol partitioned.

Packing properties of 1-alkanols and alkanes in a phospholipid membrane

We have used vibrating tube densitometry to investigate the packing properties of four alkanes and a homologous series of ten alcohols in fluid-phase membranes of dimyristoyl phosphatidylcholine (DMPC). It was found that the volume change of transferring these compounds from their pure states into the membrane, DeltaV(m)(pure-->mem), was positive for small (C4-C6) 1-alkanols while it was negative for larger alcohols and all alkanes. The magnitude of DeltaV(m)(pure-->mem) ranged from about +4 cm3/mol for alcohols with an alkyl chain about half the length of the fatty acids of DMPC, to -10 to -15 cm3/mol for the alkanes and long chain alcohols. On the basis of these observations, previously published information on the structure of the membrane-solute complexes and the free volume properties of (pure) phospholipid membranes, we suggest that two effects dominate the packing properties of hydrophobic solutes in DMPC. First, perturbation of the tightly packed interfacial zone around the ester bonds and first few methylene groups of DMPC brings about a positive contribution to DeltaV(m)(pure-->mem). This effect dominates the volume behavior for alcohols like 1-butanol, 1-pentanol and 1-hexanol. More hydrophobic solutes penetrate into the membrane core, which is loosely packed. In this region, they partially occupy interstitial (or free-) volume, which bring about a denser molecular packing and generate a negative contribution to DeltaV(m)(pure-->mem).

A kinetic and thermodynamic study on hydrolysis of sodium laurate in aqueous phase accompanied by transfer into oil phases containing different organic additives (I)

The kinetic and thermodynamic behavior at the interface between an aqueous solution of sodium laurate (NaLA) and various oil phases comprised primarily of benzene (Bz) and/or different organic compounds including amphiphiles has been investigated in regard to the hydrolysis of NaLA accelerated at the interface, transfer of lauric acid (LA) into oil phase and reverse transfer of Bz into aqueous phase in addition to interface tension. The contact of aqueous NaLA solution with the oil phase was found to accompany the mass transfer of LA and simultaneously promote the hydrolysis of NaLA in water phase. Analysis of the change of OH- ion concentration ([OH-]) over time allowed us to treat the events as a first order reaction. From the rate constant data the activation parameters such as the activation enthalpy and entropy, both of which control the transfer of LA molecules, were determined. The parameters were found to depend greatly on varied situations of the oil phase, being clearly able to explain the physicochemical behavior of the interface. Comparing the cases where the oil phase is one of the respective single systems such as Bz, dodecane (C12) and dodecylbenzene (C12Bz), C12Bz resulted in the lowest rate constant. The transfer (or hydrolysis) rate was measured for the amphiphile-added oil systems as a function of amphiphile concentration. When 0.206 M C16OH-Bz came in contact with aqueous phase, emulsion formation at the interface layer was brought about with approximately zero activation enthalpy, leading to facile or spontaneous transfer of LA. In addition, UV absorbance representing the transfer of Bz from the oil phase to the aqueous phase also demonstrated the effects of added amphiphiles on the action of the interface.

Liposome fluidization and melting point depression by pressurized CO2 determined by fluorescence anisotropy

The influence of CO2 on the bilayer fluidity of liposomes, which are representative of model cellular membranes, was examined for the first time at the elevated pressures (up to 13.9 MPa) associated with CO2-based processing of liposomes and microbial sterilization. Fluidization and melting point depression of aqueous dipalmitoylphosphatidylcholine (DPPC) liposomes by pressurized CO2 (present as an excess phase) were studied by steady-state fluorescence anisotropy using the membrane probe 1,6-diphenyl-1,3,5-hexatriene (DPH). Isothermal experiments revealed reversible, pressure-dependent fluidization of DPPC bilayers at temperatures corresponding to near-gel (295 K) and fluid (333 K) phases at atmospheric pressure, where the gel-to-fluid phase transition (Tm) occurs at approximately 315 K. Isobaric measurements (PCO2 =1.8, 7.0, and 13.9 MPa) of DPH anisotropy demonstrate substantial melting point depression (DeltaTm = -4.8 to -18.5 K) and a large broadening of the gel-fluid phase transition region, which were interpreted using conventional theories of melting point depression. Liposome fluidity is influenced by CO2 accumulation in the hydrocarbon core and polar headgroup region, as well as the formation of carbonic acid and/or the presence of buffering species under elevated CO2 pressure.

Partitioning of ABA into bilayers of Di-saturated phosphatidylcholines as measured by DSC

Using differential scanning calorimetry, we have investigated partitioning of the plant hormone abscisic acid into a homologous series of di-saturated phosphatidylcholines increasing in chain length from C(14) to C(19). Partition coefficients calculated from the shift in T(m) range from 1280 for DiC(14)PC to 480 for DiC(19)PC. The free energy of transfer is chain-length independent with a value of DeltaG = -17.4 kJ/mol and an enthalpic contribution of DeltaH = -22.6 kJ/mol. The low net entropic contribution of -TDeltaS = -5.2 J/mol agrees with the concept of the bilayer effect, but differs from that of the entropy-driven classic hydrophobic effect valid for partitioning between bulk solvents. Preferential location of the hormone in the outer region of the membrane is indicated by characteristic changes in the transition profiles and by comparison with partitioning into organic solvents whose dielectric constants model the interior and exterior regions of the bilayer. Differences in partitioning and surface pKa between the biologically active ct-ABA and the inactive tt-isomer are discussed for biological relevance.

An index of lipid phase diagrams

There is a growing awareness of the utility of lipid phase behavior data in studies of membrane-related phenomena. Such miscibility information is commonly reported in the form of temperature-composition (T-C) phase diagrams. The current index is a conduit to the relevant literature. It lists lipid phase diagrams, their components and conditions of measurement, and complete bibliographic information. The main focus of the index is on lipids of membrane origin where water is the dispersing medium. However, it also includes records on acylglycerols, fatty acids, cationic lipids, and detergent-containing systems. The miscibility of synthetic and natural lipids with other lipids, with water, and with biomolecules (proteins, nucleic acids, carbohydrates, etc.) and non-biological materials (drugs, anesthetics, organic solvents, etc.) is within the purview of the index. There are 2188 phase diagram records in the index, the bulk (81%) of which refers to binary (two-component) T-C phase diagrams. The remainder is made up of more complex (ternary, quaternary) systems, pressure-T phase diagrams, and other more exotic miscibility studies. The index covers the period from 1965 through to July, 2001.

Solubilization study of local anesthetics into sodium dodecyl sulfate micelle using anesthetic cation selective electrodes

The free concentrations of local anesthetic cations in equilibrium with sodium dodecyl sulfate (SDS) micelle which solubilized the anesthetic were determined by using ion-selective electrodes sensitive to local anesthetics, procaine (PC), lidocaine (LC), and mepivacaine (MC). Solubilizate distribution between water and SDS micelle was analyzed by means of the stepwise mass-action model. Association constant, K(1), was found to depend upon the anesthetic concentration, which decreased exponentially as the concentration of free anesthetic increased. Therefore, K(1) should include the interaction function &phi;(A) as K(1)=K(int)exp{-&phi;(A)} where K(int) is an intrinsic association constant. &phi;(A) is an increasing function of the anesthetic concentration, which means that occupation of a solubilization site by a local anesthetic cation makes sequential solubilization more difficult. The binding affinity of an anesthetic with SDS micelle increased in the following order PC<LC<MC.The critical micelle concentration (CMC) of mixed micelle was determined as a function of the concentration of free anesthetic. The CMC decreased with an increasing amount of anesthetics solubilized. All the anesthetic compositions in the micelle calculated thermodynamically from the CMC data were larger than the corresponding ones in the aqueous phase. Although the local anesthetics used here do not form micelles by themselves, the CMC vs composition curve can be regarded as a part of a micellar phase diagram showing the negative azeotropic behavior, which reflects the attractive interaction between the anionic surfactant micelle and the local anesthetic cation.

Alkane derivative-bacteriorhodopsin interaction: proton transport and protein structure

The effects of alkane derivatives, 1-alcohols (ROH), aliphatic amine hydrochlorides (RNH(2).HCl) and sodium aliphatic carboxylates (ROONa), on the proton pumping activity of bacteriorhodopsin (bR) in a purple membrane have been examined. Photocurrents in bR in the purple membrane adsorbed onto polyester thin film were recorded before and after exposure to these test substances. The peak photocurrent in bR was reversibly suppressed by each substance. From the dose-response curve, the concentrations required to reduce the peak capacitive current by 50% were determined for each homolog and then the standard free energies per CH(2) for the adsorption of the alkane derivatives to the site of action were estimated: -3.13 kJ mol(-1) for ROH, -3.05 kJ mol(-1) for RNH(3)(+), and -2.95 kJ mol(-1) for ROO(-). The proton pumping activity of bR was mainly suppressed by the hydrophobic interaction with the additive. The relative potencies of the functional groups of the alkane derivatives were almost comparable between 1-octanol (C(8)OH) and octylamine hydrochloride (C(8)NH(3)(+)) and about 10 times less effective for sodium octanoate (C(8)OO(-)) than for others. The addition of C(8)OH or C(8)OO(-) changed the absorption spectra of bR with a maximum at 560 nm to the spectra of the intermediate state with a maximum at 480 nm, while the C(8)NH(3)(+) decreased the intensity of the 560 nm band only with no blue-shift by the 480 nm band. We conclude that the action of the alkane derivatives is nonspecific and directed to all organized purple membrane structures and that the binding sites of the ROH and ROO(-) are different from that of RNH(3)(+).

Partition coefficients of charged and uncharged local anesthetics into dipalmitoylphosphatidylcholine bilayer membrane: estimation from pH dependence on the depression of phase transition temperatures

Effects of the local anesthetics, dibucaine, bupivacaine and lidocaine on the phase transition temperatures of dipalmitoylphosphatidylcholine (DPPC) bilayer membrane were studied by the optical method. We focus our attention on pH dependence of the depression of main transition and pretransition temperatures. The temperatures of both transitions of DPPC bilayer membrane were depressed by the addition of anesthetics; the higher the value of pH, the larger the depression of main transition temperature and/or pretransition temperature by anesthetics. By extending the colligative thermodynamic framework to the depression of main transition temperature by an anesthetic, we can estimate the differential partition coefficient, which is defined by the difference in partition coefficients of an anesthetic into the ripple gel and liquid crystal phases. The difference in partition coefficient between the lamellar and ripple gel phases can also be estimated from the depression of pretransition temperature. Since the differential partition coefficients include both contributions of the charged and uncharged anesthetics, we could estimate the partition coefficients of the charged and uncharged anesthetic into the membranes from the pH dependence of differential partition coefficients. The liquid crystalline membrane of DPPC bilayer was more receptive to the uncharged local anesthetics than the charged species. The partition coefficients of the charged and uncharged anesthetics into the liquid crystalline phase of DPPC bilayer membrane were 3540 and 249000 (for dibucaine), 1120 and 83900 (for bupivacaine), 256 and 11700 (for lidocaine), respectively. The transfer free energy of uncharged anesthetics from the aqueous phase to the liquid crystalline membrane was well correlated to the local anesthetic potency.

Binding of small alcohols to a lipid bilayer membrane: does the partitioning coefficient express the net affinity?

The total vapor pressures at 26 degreesC of binary (water-alcohol) and ternary (water-alcohol-vesicle) systems were measured for six short chain alcohols. The vesicles were unilamellar dipalmitoyl phosphatidylcholine (DMPC). The data was used to evaluate the effect of vesicles on the chemical potential of alcohols expressed as the preferential binding parameter of the alcohol-lipid interaction, gamma23. This quantity is a thermodynamic (model-free) measure of the net strength of membrane-alcohol interactions. For the smaller investigated alcohols (methanol, ethanol and 1-propanol) gamma23 was negative. This is indicative of so-called preferential hydration, a condition where the affinity of the membrane for water is higher than the affinity for the alcohol. For the longer alcohols (1-butanol, 1-pentanol, 1-hexanol) gamma23 was positive and increasing with increasing chain length. This demonstrates preferential binding, i.e. enrichment of alcohol in the membrane and a concomitant depletion of the solute in the aqueous bulk. The measured values of gamma23 were compared to the number of alcohol-membrane contacts specified by partitioning coefficients from the literature. It was found that for the small alcohols the number of alcohol-membrane contacts is much larger than the number of preferentially bound solutes. This discrepancy, which is theoretically expected in cases of very weak binding, becomes less pronounced with increasing alcohol chain length, and when the partitioning coefficient exceeds approximately 3 on the molal scale (10(2) in mole fraction units) it vanishes. Based on this, relationships between structural and thermodynamic interpretations of membrane partitioning are discussed.

Effect of local anesthetics on the bilayer membrane of dipalmitoylphosphatidylcholine: interdigitation of lipid bilayer and vesicle-micelle transition

The phase transitions of dipalmitoylphosphatidylcholine (DPPC) bilayer membrane were observed by means of differential scanning calorimetry (DSC) as a function of the concentration of local anesthetics, dibucaine (DC x HCl), tetracaine (TC x HCl), lidocaine (LC x HCl) and procaine hydrochlorides (PC x HCl). LC x HCl and PC x HCl depressed monotonously the temperatures of the main- and pre-transition of DPPC bilayer membrane. The enthalpy changes of both transitions decreased slightly with an increase in anesthetic concentration up to 160 mmol kg(-1). In contrast, the addition of TC x HCl or DC x HCl, having the ability to form a micelle by itself, induced the complex phase behavior of DPPC bilayer membrane including the vesicle-to-micelle transition. The depression of both temperatures of the main- and pre-transition, which is accompanied with a decrease in enthalpy, was observed by the addition of TC x HCl up to 21 mmol kg(-1) or DC x HCl up to 11 mmol kg(-1). The pretransition disappeared when these concentrations of anesthetic were added, and the interdigitated gel phase appeared above these concentrations. The appearance of the interdigitated gel phase, instead of the ripple gel phase, brings about the stabilization of the gel phase by 1.8-2.4 kcal mol(-1). In the concentration range of 70-120 mmol kg(-1) TC x HCl (or 40-60 mmol kg(-1) DC x HCl), the enthalpy of the main transition exhibited a drastic decrease, resulting in the virtual disappearance of the main transition. This process includes the decrease in vesicle size with increasing anesthetic concentration, resulting in the mixed micelle of DPPC and anesthetics. Therefore, in this range of anesthetic concentration, the DPPC vesicle solubilized an anesthetic which coexists with the DPPC-anesthetic mixed micelle. Above the concentration of 120 mmol kg(-1) TC x HCl (or 60 mmol kg(-1) DC x HCl), there exists the DPPC-anesthetic mixed micelle. Two types of new transitions concerned with the mixed micelle of DPPC and micelle-forming anesthetics were observed by DSC.

Unifying characteristics of sites of anesthetic action revealed by combined use of anesthetics and non-anesthetics

1. The usefulness of nonanesthetics in the study of mechanisms of general anesthesia lies in the possibility to identify the unifying characteristics of molecular sites that are shared by the anesthetics but not by the structurally similar nonanesthetics. 2. In model membranes, pairs of structurally similar anesthetics and nonanesthetics showed distinctly different submolecular distributions. 3. This difference may be the underlying cause for the different anesthetic and nonanesthetic interaction with gramicidin A, a model transmembrane cation channel. 4. Generalization of our findings suggests that the nature of the sites, whether in lipids or proteins, must be neither extremely hydrophilic nor extremely lipophilic, but amphiphilic.

Membrane-buffer partition coefficients of a local anesthetic tetracaine monitored by an anesthetic sensor; effects of temperature and pH

Binding of a local anesthetic tetracaine (TC) to dimyristoylphosphatidylcholine (DMPC) bilayer membrane was studied by the potentiomerty with an ion-selective electrode sensitive to TC cation. DMPC membrane-buffer partition coefficient (K(app)) was determined in mole fraction unit as a function of pH for the lamellar gel (at 12 degrees C), ripple gel (at 20 degrees C), and liquid crystal (at 30 degrees C) phases. The partition coefficients of charged (K+) and uncharged TC (K0) into the DMPC membranes were estimated from the pH-dependence of K(app). The three states of DMPC membranes were more receptive to the uncharged TC than the charged species.

Phases and phase transitions of the phosphatidylcholines

LIPIDAT (http://www.lipidat.chemistry.ohio-state.edu) is an Internet accessible, computerized relational database providing access to the wealth of information scattered throughout the literature concerning synthetic and biologically derived polar lipid polymorphic and mesomorphic phase behavior and molecular structures. Here, a review of the data subset referring to phosphatidylcholines is presented together with an analysis of these data. This subset represents ca. 60% of all LIPIDAT records. It includes data collected over a 43-year period and consists of 12,208 records obtained from 1573 articles in 106 different journals. An analysis of the data in the subset identifies trends in phosphatidylcholine phase behavior reflecting changes in lipid chain length, unsaturation (number, isomeric type and position of double bonds), asymmetry and branching, type of chain-glycerol linkage (ester, ether, amide), position of chain attachment to the glycerol backbone (1,2- vs. 1,3-) and head group modification. Also included is a summary of the data concerning the effect of pressure, pH, stereochemical purity, and different additives such as salts, saccharides, amino acids and alcohols, on phosphatidylcholine phase behavior. Information on the phase behavior of biologically derived phosphatidylcholines is also presented. This review includes 651 references.

The effect of ethanol on the structure of phosphatidylserine bilayers

Thermotrophic and structural effects of ethanol on phosphatidylserine (PS) membranes were investigated by differential scanning calorimetry (DSC) and X-ray diffraction. It was found that up to 15% (v/v) added ethanol, there is little change in the melting temperature of the phospholipid and no change in the interbilayer (d) spacing in the gel phase, indicating that there is no interdigitation of the hydrocarbon chains. Above the melting temperature of the phospholipid, a large decrease of the d spacing, due primarily to a decrease in the thickness of the bilayer, was found. Ethanol molecules located in the headgroup region apparently expand the area available to the headgroups with concomitant coiling of the acyl chains, resulting in marked thinning of the lipid layer.

Different distribution of fluorinated anesthetics and nonanesthetics in model membrane: a 19F NMR study

Despite their structural resemblance, a pair of cyclic halogenated compounds, 1-chloro-1,2,2-trifluorocyclobutane (F3) and 1,2-dichlorohexafluorocyclobutane (F6), exhibit completely different anesthetic properties. Whereas the former is a potent general anesthetic, the latter produces no anesthesia. Two linear compounds, isoflurane and 2,3-dichlorooctofluorobutane (F8), although not a structural pair, also show the same anesthetic discrepancy. Using 19F nuclear magnetic spectroscopy, we investigated the time-averaged submolecular distribution of these compounds in a vesicle suspension of phosphatidylcholine lipids. A two-site exchange model was used to interpret the observed changes in resonance frequencies as a function of the solubilization of these compounds in membrane and in water. At clinically relevant concentrations, the anesthetics F3 and isoflurane distributed preferentially to regions of the membrane that permit easy contact with water. The frequency changes of these two anesthetics can be well characterized by the two-site exchange model. In contrast, the nonanesthetics F6 and F8 solubilized deeply into the lipid core, and their frequency change significantly deviated from the prediction of the model. It is concluded that although anesthetics and nonanesthetics may show similar hydrophobicity in bulk solvents such as olive oil, their distributions in various regions in biomembranes, and hence their effective concentrations at different submolecular sites, may differ significantly.

Influence of membrane lipid packing on T2-weighted magnetic resonance images: study of relaxation parameters in model membrane systems

The detailed mechanisms leading to soft tissue contrast in MRI are not known. To explore the physical basis for this phenomena, the effect of acyl chain packing of phospholipid model membranes on water proton relaxation were investigated at 20 MHz at 40 degrees C. Three variables affecting lipid packing were examined: chain composition, toluene concentration, and pH/ionic strength. It is demonstrated that the T2 of water protons is strongly influenced by changes in lipid packing, while T1 is generally insensitive to these variables. It is also shown that small changes in water content (7 wt%) altered both T1 and T2 in all systems studied. The results suggest that the physical properties of biological membranes in tissues are an important factor in determining T2.

Alcohol binding to liposomes by 2H NMR and radiolabel binding assays: does partitioning describe binding?

Implicit within the concept of membrane-buffer partition coefficients of solutes is a nonspecific solvation mechanism of solute binding. However, (2)H NMR studies of the binding of (2)H(6)-ethanol and [1-(2)H(2)] n-hexanol to phosphatidylcholine vesicles have been interpreted as evidence for two distinct alcohol binding modes. One binding mode was reported to be at the membrane surface. The second mode was reported to be within the bilayer interior. An examination of the (2)H NMR binding studies, together with direct radiolabel binding assays, shows that other interpretations of the data are more plausible. The results are entirely consistent with partitioning (nonspecific binding) as the sole mode of alcohol binding to liposomes, in accord with our previous thermodynamic interpretation of alcohol action in phosphatidylcholine liposomes.

Primary alcohols and phosphatidylcholine metabolism in rat brain synaptosomal membranes via phospholipase D

Phospholipase D of rat brain synaptosomal membranes was tested with phosphatidylcholine as the substrate for its specificity in the use of primary alcohols as transphosphatidylation co-substrates. The efficiency of the reaction was related to the hydrophobicity and the membrane penetrating capacity of the alcohol molecule. Phosphatidylalcohol formation could be detected up to 1-octanol but not for alcohols with longer hydrocarbon chains (C(9), C(10)). With increasing alcohol concentration the transphosphatidylation activity of the phospholipase D reached an optimum and then declined abruptly. Alcohol concentrations required for maximal transphosphatidylation reaction generally decreased with increasing hydrophobicities of the alcohols. Nevertheless 1-butanol and 4-chloro-1-butanol were the most efficient cosubstrates, sharing identical optimal conditions. Transphosphatidylation works at the cost of phosphatidic acid formation. Phosphatidic acid itself was transformed to diacylglycerol, probably by a contaminating phosphatidic acid phosphohydrolase.

Concentration-dependent binding of the chiral beta-blocker oxprenolol to isoelectric or negatively charged unilamellar vesicles

Large unilamellar vesicles (LUVs) of different lipid compositions were used to study the type of binding of the beta-blocking cationic agent oxprenolol to the lipid matrix of biological membranes at a physiologic pH value of 7.4. When isoelectric membranes of pure egg lecithin or egg lecithin/cholesterol (7:3 mol/mol) were used, a linear relationship between membrane-bound and free oxprenolol indicated a constant molar partition coefficient of 54 or 44 between the liposomal and the aqueous phase over a wide concentration range of the drug up to 25 mM. This pointed to deep insertion of the drug molecules into the hydrophobic membrane interior. Drug binding to membranes of negatively charged phosphatidylserine from bovine brain was cooperative with a Hill coefficient h of 3.4 at concentrations below 0.5 mM and a molar ratio Re of bound drug per lipid of 1:10. Above drug concentrations of 2.5 mM and Re = 1:5, a constant molar partition coefficient of 33 could be estimated. R-oxprenolol or S-oxprenolol, as well as the racemic drug, showed no differences in membrane binding, even with egg lecithin LUVs containing 20 mol% of the negatively charged (2S, 4R)-N-(hexadecanoyl)-4-hydroxyproline, which has a pronounced chiral headgroup. Our results suggest that enantioselective interactions of the chiral oxprenolol with the chiral lipids of biological membranes can be excluded. Furthermore, surface adsorption of the drug is probable only on the negatively charged cytosolic side of biological plasma membranes, whereas on the isoelectric exterior the cationic drug is inserted deeply into the membrane.

Multiphasic modelling of ligand/acceptor interactions. The hydrophobicity-dependent binding of relatively small amphiphilic substances to acceptor proteins and the nature and facedness of acceptor sites

The modelling of multiphasic ligand/acceptor equilibrium binding systems proceeds at three logically distinct levels: (1) A suitable response quantity, e.g. the amount of acceptor-bound ligand nEL, is expressed as a function of the ligand concentrations [Li] (L = A,B,...) in the compartment i that contains the acceptor sites. One thus obtains a response function nEL = f1([Li]). In general, the equilibrium constants KL contained in such mathematical models are physically ill-defined. (2) Each local concentration [Li] is further expressed as a function of [Laq], the corresponding concentration in the aqueous phase, leading to nEL = f2([Laq]). In this way, the constants KL are transformed into effective constants K'L which (i) can be assessed experimentally and (ii) depend on ligand hydrophobicity in a way that is characteristic of the binding site. Formulation of the functions f1 and F2 only requires knowledge of the reactions in which the acceptor sites participate directly. (3) For each ligand, the experimentally accessible total ligand concentration Lt is expressed as a function of [Laq], leading to concentration balance equations Lt = Lt([Laq]). The latter transformation takes account of any reactions, distinct from ligand/acceptor interaction, in which the ligands are involved, e.g. binding to additional protein sites. As a result of steps 2 and 3, each binding system is described by a set of simultaneous equations dependent on the auxiliary variable [Laq]: (i) the response function f2([Laq]) and (ii) a concentration balance for each ligand Lt = Lt([Laq]). The formulae are rendered more conscise and their discussion and application to data fitting are simplified by introducing, for each ligand L, a function FL characterising the distribution of unbound monomeric ligand over the various partition compartments. When the acceptor acts on unbound ligand, the formulae are further expressed in terms of a new auxiliary variable i.e. the total concentration of unbound monomeric ligand microL. In contrast to data analysis as a function of local concentrations, analysis in terms of total ligand concentrations avoids losing sight of alternate hypotheses about the nature of the binding sites. The present formulation has also permitted clarification of several consequences of the multiphasic nature of the binding systems that, as yet, have been poorly recognised.(ABSTRACT TRUNCATED AT 400 WORDS)

An approximate model and empirical energy function for solute interactions with a water-phosphatidylcholine interface

An empirical model of a liquid crystalline (L alpha phase) phosphatidylcholine (PC) bilayer interface is presented along with a function which calculates the position-dependent energy of associated solutes. The model approximates the interface as a gradual two-step transition, the first step being from an aqueous phase to a phase of reduced polarity, but which maintains a high enough concentration of water and/or polar head group moieties to satisfy the hydrogen bond-forming potential of the solute. The second transition is from the hydrogen bonding/low polarity region to an effectively anhydrous hydrocarbon phase. The "interfacial energies" of solutes within this variable medium are calculated based upon atomic positions and atomic parameters describing general polarity and hydrogen bond donor/acceptor propensities. This function was tested for its ability to reproduce experimental water-solvent partitioning energies and water-bilayer partitioning data. In both cases, the experimental data was reproduced fairly well. Energy minimizations carried out on beta-hexyl glucopyranoside led to identification of a global minimum for the interface-associated glycolipid which exhibited glycosidic torsion angles in agreement with prior results (Hare, B.J., K.P. Howard, and J.H. Prestegard. 1993. Biophys. J. 64:392-398). Molecular dynamics simulations carried out upon this same molecule within the simulated interface led to results which were consistent with a number of experimentally based conclusions from previous work, but failed to quantitatively reproduce an available NMR quadrupolar/dipolar coupling data set (Sanders, C.R., and J.H. Prestegard. 1991. J. Am. Chem. Soc. 113:1987-1996). The proposed model and functions are readily incorporated into computational energy modeling algorithms and may prove useful in future studies of membrane-associated molecules.