Phospholipid complexation of general anesthetics in fluid bilayers

How to Determine Lipid Interactions in Membranes from Experiment Through the Ising Model

The determination and the meaning of interactions in lipid bilayers are discussed and interpreted through the Ising model. Originally developed to understand phase transitions in ferromagnetic systems, the Ising model applies equally well to lipid bilayers. In the case of a membrane, the essence of the Ising model is that each lipid is represented by a site on a lattice and that the interaction of each site with its nearest neighbors is represented by an energy parameter ω. To calculate the thermodynamic properties of the system, such as the enthalpy, the Gibbs energy, and the heat capacity, the partition function is derived. The calculation of the configurational entropy factor in the partition function, however, requires approximations or the use of Monte Carlo (MC) simulations. Those approximations are described. Ultimately, MC simulations are used in combination with experiment to determine the interaction parameters ω in lipid bilayers. Several experimental approaches are described, which can be used to obtain interaction parameters. They include nearest-neighbor recognition, differential scanning calorimetry, and Förster resonance energy transfer. Those approaches are most powerful when used in combination of MC simulations of Ising models. Lipid membranes of different compositions are discussed, which have been studied with these approaches. They include mixtures of cholesterol, saturated (ordered) phospholipids, and unsaturated (disordered) phospholipids. The interactions between those lipid species are examined as a function of molecular properties such as the degree of unsaturation and the acyl chain length. The general rule that emerges is that interactions between different lipids are usually unfavorable. The exception is that cholesterol interacts favorably with saturated (ordered) phospholipids. However, the interaction of cholesterol with unsaturated phospholipids becomes extremely unfavorable as the degree of unsaturation increases.

Heat Capacity of DPPC/Cholesterol Mixtures: Comparison of Single Bilayers with Multibilayers and Simulations

The excess heat capacity (Δ C _p) of mixtures of dipalmitoylphosphatidylcholine (DPPC) and cholesterol (Chol) is examined in detail in large unilamellar vesicles (LUVs), both experimentally, using differential scanning calorimetry (DSC), and theoretically, using a three-state Ising model. The model postulates that DPPC can access three conformational states: gel, liquid-disordered (L_d), and liquid-ordered (L_o). The L_o state, however, is only available if coupled with interaction with an adjacent Chol. Δ C _p was calculated using Monte Carlo simulations on a lattice and compared to experiment. The DSC results in LUVs are compared with literature data on multilamellar vesicles (MLVs). The enthalpy change of the complete phase transition from gel to L_d is identical in LUVs and MLVs, and the melting temperatures ( T_m) are similar. However, the DSC curves in LUVs are significantly broader, and the maxima of Δ C _p are accordingly smaller. The parameters in the Ising model were chosen to match the DSC curves in LUVs and the nearest-neighbor recognition (NNR) data. The model reproduces the NNR data very well. It also reproduces the phase transition in DPPC, the freezing point depression induced by Chol, and the broad component of Δ C _p in DPPC/Chol LUVs. However, there is a sharp component, between 5 and 15 mol % Chol, that the model does not reproduce. The broad component of Δ C _p becomes dominant as Chol concentration increases, indicating that it involves melting of the L_o phase. Because the simulations reproduce this component, the conclusions regarding the nature of the phase transition at high Chol concentrations and the structure of the L_o phase are important: there is no true phase separation in DPPC/Chol LUVs. There are large domains of gel and L_o phase coexisting below T_m of DPPC, but above T_m the three states of DPPC are mixed with Chol, although clusters persist.

Binding of serotonin to lipid membranes

Serotonin (5-hydroxytryptamine, 5-HT) is a prevalent neurotransmitter throughout the animal kingdom. It exerts its effect through the specific binding to the serotonin receptor, but recent research has suggested that neural transmission may also be affected by its nonspecific interactions with the lipid matrix of the synaptic membrane. However, membrane-5-HT interactions remain controversial and superficially investigated. Fundamental knowledge of this interaction appears vital in discussions of putative roles of 5-HT, and we have addressed this by thermodynamic measurements and molecular dynamics (MD) simulations. 5-HT was found to interact strongly with lipid bilayers (partitioning coefficient ~1200 in mole fraction units), and this is highly unusual for a hydrophilic solute like 5-HT which has a bulk, oil-water partitioning coefficient well below unity. It follows that membrane affinity must rely on specific interactions, and the MD simulations identified the salt-bridge between the primary amine of 5-HT and the lipid phosphate group as the most important interaction. This interaction anchored cationic 5-HT in the membrane interface with the aromatic ring system pointing inward and a prevailing residence between the phosphate and the carbonyl groups of the lipid. The unprotonated form of 5-HT shows the opposite orientation, with the primary amine pointing toward the membrane core. Partitioning of 5-HT was found to decrease lipid chain order. These distinctive interactions of 5-HT and model membranes could be related to nonspecific effects of this neurotransmitter.

Atomistic study of lipid membranes containing chloroform: looking for a lipid-mediated mechanism of anesthesia

The molecular mechanism of general anesthesia is still a controversial issue. Direct effect by linking of anesthetics to proteins and indirect action on the lipid membrane properties are the two hypotheses in conflict. Atomistic simulations of different lipid membranes subjected to the effect of small volatile organohalogen compounds are used to explore plausible lipid-mediated mechanisms. Simulations of homogeneous membranes reveal that electrostatic potential and lateral pressure transversal profiles are affected differently by chloroform (anesthetic) and carbon tetrachloride (non-anesthetic). Simulations of structured membranes that combine ordered and disordered regions show that chloroform molecules accumulate preferentially in highly disordered lipid domains, suggesting that the combination of both lateral and transversal partitioning of chloroform in the cell membrane could be responsible of its anesthetic action.

The Role of Atomic Polarization in the Thermodynamics of Chloroform Partitioning to Lipid Bilayers

In spite of extensive research and use in medical practice, the precise molecular mechanism of volatile anesthetic action remains unknown. The distribution of anesthetics within lipid bilayers and potential targeting to membrane proteins is thought to be central to therapeutic function. Therefore, obtaining a molecular level understanding of volatile anesthetic partitioning into lipid bilayers is of vital importance to modern pharmacology. In this study we investigate the partitioning of the prototypical anesthetic, chloroform, into lipid bilayers and different organic solvents using molecular dynamics simulations with potential models ranging from simplified coarse-grained MARTINI to additive and polarizable CHARMM all-atom force fields. Many volatile anesthetics display significant inducible dipole moments, which correlate with their potency, yet the exact role of molecular polarizability in their stabilization within lipid bilayers remains unknown. We observe that explicit treatment of atomic polarizability makes it possible to accurately reproduce solvation free energies in solvents with different polarities, allowing for quantitative studies in heterogeneous molecular distributions, such as lipid bilayers. We calculate the free energy profiles for chloroform crossing lipid bilayers to reveal a role of polarizability in modulating chloroform partitioning thermodynamics via the chloroform-induced dipole moment and highlight competitive binding to the membrane core and toward the glycerol backbone that may have significant implications for understanding anesthetic action.

Physical and physicochemical factors effecting transport of chlorohydrocarbon gases from lung alveolar air to blood as measured by the causation of narcosis

This systematic investigation examines gas transport in the lung for two sets of chlorohydrocarbons (CHCs): the chloromethanes (C1) and chloroethanes (C2). The C1 series includes chloromethane, methylene chloride, chloroform, and carbon tetrachloride, and the C2 series includes chloroethane, 1,2-dichloroethane, 1, 1, 2-trichloroethane, and 1, 1, 2, 2-tetrachloroethane. Most CHC gases cause narcosis. The comprehensive narcosis work of Lehmann and colleagues on CHCs was used as a basis for the narcosis endpoint in the present examination. The sites for narcosis are located in the brain (midline cortex and posterior parietal area), the spine, and at many peripheral nerve sites. Central nervous system (CNS) exposure executes a multisite, neural transmission set of inhibitions that promotes rapid loss of consciousness, sensory feeling, and current and stored memory while providing temporary amnesia. Absorption into the system requires dissolution into many lipid membranes and binding to lipoproteins. Lipophilicity is a CHC property shared with many anesthetics according to the Meyer-Overton Rule. Many structurally different lipid chemicals produce the narcosis response when the lipid concentration exceeds -67 mM. This suggests narcotic or anesthetic dissolution into CNS membranes until the lipid organization is disrupted or perturbed. This perturbation includes loading of Na(+)- and K(+)-channel transmembrane lipoprotein complexes and disrupting their respective channel functional organizations. The channel functions become attenuated or abrogated until the CHC exposure ceases and CHC loading reverses. This investigation demonstrates how the CHC physical and chemical properties influence the absorption of these CHCs via the lung and the alveolar system on route to the blood. Narcosis in test animals was used here as an objective biological endpoint to study the effects of the physical factors Bp, Vp, Kd (oil: gas) partition, Henry's constant (HK), and water solubility (S%) on gas transport. Narcosis is immediate after gas exposure and requires no chemical activation only absorption into the blood and circulation to CNS narcotic sites. The three physical factors Bp, K(d) (oil: air), and S% vary directly with unitary narcosis (UN) whereas Vp and HK vary inversely with UN in linear log-log relationships for the C2 series but not for the C1 series. Physicochemical properties of C1 series gases indicate why they depart from what is usually assumed to be an Ideal Gas. An essential discriminating process in the distal lung is the limiting alveolar film layer (AFL) and the membrane layer of the alveolar acini. The AFL step influences gas uptake by physically limiting the absorption process. Interaction with and dissolution into aqueous solvent of the AFL is required for transport and narcotic activity. Narcotics or anesthetics must engage the aqueous AFL with sufficient strength to allow transport and absorption for downstream CNS binding. CHCs that do not engage well with the AFL are not narcotic. Lipophilicity and amphipathicity are also essential solvency properties driving narcotics' transport through the alveolar layer, delivery to the blood fats and lipoproteins, and into critical CNS lipids, lipoproteins, and receptor sites that actuate narcosis. AFL disruption is thought to be strongly related to a number of serious pulmonary diseases such acute respiratory distress syndrome, infant respiratory distress syndrome, emphysema, chronic obstructive pulmonary disease, asthma, chronic bronchitis, pneumonia, pulmonary infections, and idiopathic pulmonary fibrosis. The physical factors (Bp, Vp, Kd [oil: gas] partition, Henry's constant, and water solubility [S%]) combine to affect a specific transport through the AFL if lung C > C(0) (threshold concentration for narcosis). The degree of blood CHC absorption depends on dose, lipophilicity, and lung residence time. AFL passage can be manipulated by physical factors of increased pressure (kPa) or increased gas exposure (moles). Molecular lipophilicity facilitates narcosis but lipophilicity alone does not explain narcosis. Vapor pressure is also required for narcosis. Narcotic activity apparently requires stereospecific processing in the AFL and/or down-stream inhibition at stereospecific lipoproteins at CNS inhibitory sites. It is proposed that CHCs likely cannot proceed through the AFL without perturbation or disruption of the integrity of the AFL at the alveoli. CHC physicochemical properties are not expected to allow their transport through the AFL as physiological CO(2) and O(2) naturally do in respiration. This work considers CHC inspiration and systemic absorption into the blood with special emphasis on the CHC potential perturbation effects on the lipid, protein liquid layer supra to the alveolar membrane (AFL). A heuristic gas transport model for the CHCs is presented as guidance for this examination. The gas transport model can be used to study absorption for other gas delivery endpoints of environmental concern such as carcinogens.

Effects of isoflurane, halothane, and chloroform on the interactions and lateral organization of lipids in the liquid-ordered phase

The first quantitative insight has been obtained into the effects that volatile anesthetics have on the interactions and lateral organization of lipids in model membranes that mimic "lipid rafts". Specifically, nearest-neighbor recogntion measurements, in combination with Monte Carlo simulations, have been used to investigate the action of isoflurane, halothane, and chloroform on the compactness and lateral organization of cholesterol-rich bilayers of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) in the liquid-ordered (l(o)) phase. All three anesthetics induce a similar weakening of sterol-phospholipid association, corresponding to ca. 30 cal/mol of lipid at clinically relevant concentrations. Monte Carlo lattice simulations show that the lateral organization of the l(o) phase, under such conditions, remains virtually unchanged. In sharp contrast to their action on the l(o) phase, these anesthetics have been found to have a similar strengthening effect on sterol-phospholipid association in the liquid-disordered (l(d)) phase. The possibility of discrete complexes being formed between DPPC and these anesthetics and the biological relevance of these findings are discussed.