Interactions of beta-blockers with model lipid membranes: Molecular view of the interaction of acebutolol, oxprenolol, and propranolol with phosphatidylcholine vesicles by time-dependent fluorescence shift and molecular dynamics simulations

Time-resolved fluorescence of tryptophan characterizes membrane perturbation by cyclic lipopeptides

Viscosin is a membrane-permeabilizing, cyclic lipopeptide (CLiP) produced by Pseudomonas species. Here, we have studied four synthetic analogs (L1W, V4W, L5W, and L7W), each with one leucine (Leu; L) or valine residue exchanged for tryptophan (Trp; W) by means of time-resolved fluorescence spectroscopy of Trp. To this end, we recorded the average fluorescence lifetime, rotational correlation time and limiting anisotropy, dipolar relaxation time and limiting extent of relaxation, rate constant of acrylamide quenching, effect of H2O-D2O exchange, and time-resolved half-width of the spectrum in the absence and presence of POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) liposomes. Structure, localization, and hydration of the peptides were described by molecular dynamics simulations. The combination of the parameters provides a good description of the molecular environments of the Trp positions and the behavior of viscosin as a whole. Of particular value for characterizing the impact of viscosin on the membrane is the dipolar relaxation of Trp4 in V4W, which is deeply embedded in the hydrophobic core. The limiting relaxation level represents the membrane perturbation-unlike typical membrane probes-at the site of the perturbant. Fractions of Trp4 relax at different rates; the one not in contact with water upon excitation relaxes via recruitment of a water molecule on the 10-ns timescale. This rate is sensitive to the concerted membrane perturbation by more than one lipopeptide, which appears at high lipopeptide concentration and is assumed a prerequisite for the final formation of a membrane-permeabilizing defect. Trp7 relaxes primarily with respect to neighboring Ser residues. Trp5 flips between a membrane-inserted and surface-exposed orientation.

Exploring the pH-Responsive Interaction of β-Blocker Drug Propranolol with Biomimetic Micellar Media: Fluorescence and Electronic Absorption Studies

Interaction of neutral and charged lipophilic beta-blocker drug, propranolol (PPL) with biomimicking nanocavities formed by micelles bearing same and opposite charges namely, cationic cetyltrimethylammonium bromide (CTAB), a surface-active ionic liquid 1-hexadecyl-3-methylimidazolium chloride (HDMIC) and anionic sodium dodecyl sulphate (SDS) have been investigated using fluorescence and absorption spectroscopic techniques. Binding of PPL to SDS at pH < pKa is characterised by biphasic interactions with decrease in fluorescence intensity at lower concentrations and subsequent increase post micellization. All the surfactants show significant interactions with the neutral drug molecule at pH > pKa, which is evident from the strongest binding constant ( K b ) values at pH 10.4. Results of quenching studies indicate that the location of drug molecule is determined by its charge, which is influenced by both pH and charge on micelle surface. For PPL-CTAB and PPL-HDMIC systems, quenching was strongest at pH 10.4, moderate at pH 7.4 and was absent at pH 3.5. However, the PPL-SDS system displayed similar K SV values at all pH conditions, suggesting that the probe is at the same position regardless of pH. Non-covalent interactions, which play crucial role in biological systems, are similarly the primary driving force governing the interaction between PPL and surfactant micelles.

Potential of liposomes and lipid membranes for the separation of β-blockers by capillary electromigration and liquid chromatographic techniques

β-Blockers belong to a frequently used class of drugs primarily used to treat heart and circulatory conditions. Here we describe the use of lipid vesicles and liposomes as cell membrane biomimicking models in capillary electromigration (CE) and liquid chromatography (LC) techniques for the investigation of interactions between lipid membranes and β-blockers. In addition to liposomes, the use of commercial intravenous lipid emulsions, and their interactions with β-blockers are also discussed. Different CE and LC instrumental techniques designed for these purposes are introduced. Other methodologies for studying interactions between β-blockers and lipid membranes are also briefly discussed, and the different methodologies are compared. The aim is to give the reader a good overview on the status of the use of liposomes and lipids in CE and LC for studying β-blocker interactions.

Effect of liposome composition on β-blocker interactions studied by capillary electrokinetic chromatography

Liposome capillary electrokinetic chromatography was used to investigate the interactions between three β-blockers of different hydrophobicity and various liposome solutions. The studied β-blockers comprised alprenolol, propranolol, and carvedilol. The composition of the liposome solutions, containing 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine, 1,2-dioleoyl-sn-glycero-3-phos-phoethanolamine, 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine, and cholesterol in various molar ratios, was designed by a response surface methodology-central composite design approach. Subsequently, after conducting the liposome capillary electrokinetic chromatography experiments and determining the retention factors from the electrophoretic mobilities of the compounds, and further calculating the distribution coefficients, an analysis of variance was performed. After extracting the statistical models, optimal operational conditions were obtained based on the developed models. To further investigate the interactions between the β-blockers and the liposomes, nanoplasmonic sensing experiments were carried out on two different liposome systems. The overall results demonstrate the strong influence of cholesterol and 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine on the distribution coefficients.

Propranolol induces large-scale remodeling of lipid bilayers: tubules, patches, and holes

Herein, we report fluorescence microscopy analysis of the interaction between propranolol (PPN), a beta-adrenergic blocking agent, and planar supported lipid bilayers (SLBs), as model membranes. The results indicate that PPN can remarkably promote largescale remodeling in SLBs with various lipid compositions. It was found that PPN insertion induces the formation of long microtubules that can retract into hemispherical caps on the surface of the bilayer. These transformations are dynamic, partially reversible, and dependent upon the drug concentration. Quantitative analysis revealed a three-step model for PPN-lipid bilayer interaction, with the first step involving interfacial electrostatic adsorption, the second step centered on hydrophobic insertion, and the third step associated with membrane disruption and hole formation. By introducing cholesterol, phosphoethanolamine, phosphatidylglycerol, and phosphatidylserine lipids into the phosphocholine SLBs, it was illustrated that both the chemistry of the lipid headgroups and the packing of lipid acyl chains can substantially affect the particular steps in the interactions between PPN and lipid bilayers. Our findings may help to elucidate the possible mechanisms of PPN interaction with lipid membranes, the toxic behavior and overdosage scenarios of beta-blockers, and provide valuable information for drug development and modification.

Effect of Ester Moiety on Structural Properties of Binary Mixed Monolayers of Alpha-Tocopherol Derivatives with DPPC

Phospholipid membranes are ubiquitous components of cells involved in physiological processes; thus, knowledge regarding their interactions with other molecules, including tocopherol ester derivatives, is of great importance. The surface pressure–area isotherms of pure α-tocopherol (Toc) and its derivatives (oxalate (OT), malonate (MT), succinate (ST), and carbo analog (CT)) were studied in Langmuir monolayers in order to evaluate phase formation, compressibility, packing, and ordering. The isotherms and compressibility results indicate that, under pressure, the ester derivatives and CT are able to form two-dimensional liquid-condensed (LC) ordered structures with collapse pressures ranging from 27 mN/m for CT to 44 mN/m for OT. Next, the effect of length of ester moiety on the surface behavior of DPPC/Toc derivatives’ binary monolayers at air–water interface was investigated. The average molecular area, elastic modulus, compressibility, and miscibility were calculated as a function of molar fraction of derivatives. Increasing the presence of Toc derivatives in DPPC monolayer induces expansion of isotherms, increased monolayer elasticity, interrupted packing, and lowered ordering in monolayer, leading to its fluidization. Decreasing collapse pressure with increasing molar ratio of derivatives indicates on the miscibility of Toc esters in DPPC monolayer. The interactions between components were analyzed using additivity rule and thermodynamic calculations of excess and total Gibbs energy of mixing. Calculated excess area and Gibbs energy indicated repulsion between components, confirming their partial mixing. In summary, the mechanism of the observed phenomena is mainly connected with interactions of ionized carboxyl groups of ester moieties with DPPC headgroup moieties where formed conformations perturb alignment of acyl chains, resulting in increasing mean area per molecule, leading to disordering and fluidization of mixed monolayer.

What Does Time-Dependent Fluorescence Shift (TDFS) in Biomembranes (and Proteins) Report on?

The organization of biomolecules and bioassemblies is highly governed by the nature and extent of their interactions with water. These interactions are of high intricacy and a broad range of methods based on various principles have been introduced to characterize them. As these methods view the hydration phenomena differently (e.g., in terms of time and length scales), a detailed insight in each particular technique is to promote the overall understanding of the stunning “hydration world.” In this prospective mini-review we therefore critically examine time-dependent fluorescence shift (TDFS)—an experimental method with a high potential for studying the hydration in the biological systems. We demonstrate that TDFS is very useful especially for phospholipid bilayers for mapping the interfacial region formed by the hydrated lipid headgroups. TDFS, when properly applied, reports on the degree of hydration and mobility of the hydrated phospholipid segments in the close vicinity of the fluorophore embedded in the bilayer. Here, the interpretation of the recorded TDFS parameters are thoroughly discussed, also in the context of the findings obtained by other experimental techniques addressing the hydration phenomena (e.g., molecular dynamics simulations, NMR spectroscopy, scattering techniques, etc.). The differences in the interpretations of TDFS outputs between phospholipid biomembranes and proteins are also addressed. Additionally, prerequisites for the successful TDFS application are presented (i.e., the proper choice of fluorescence dye for TDFS studies, and TDFS instrumentation). Finally, the effects of ions and oxidized phospholipids on the bilayer organization and headgroup packing viewed from TDFS perspective are presented as application examples.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

Lipid Hydroperoxide Compromises the Membrane Structure Organization and Softens Bending Rigidity

Lipid hydroperoxides are key mediators of diseases and cell death. In this work, the structural and dynamic perturbations induced by the hydroperoxidized POPC lipid (POPC-OOH) in fluid POPC membranes, at both 23 and 37 °C, were addressed using advanced small-angle X-ray scattering (SAXS) and fluorescence methodologies. Notably, SAXS reveals that the hydroperoxide group decreases the lipid bilayer bending rigidity. This alteration disfavors the bilayer stacking and increases the swelling in-between stacked bilayers. We further investigated the changes in the apolar/polar interface of hydroperoxide-containing membranes through time-resolved fluorescence/anisotropy experiments of the probe TMA-DPH and time-dependent fluorescence shifts of Laurdan. A shorter mean fluorescence lifetime for TMA-DPH was obtained in enriched POPC-OOH membranes, revealing a higher degree of hydration near the membrane interface. Moreover, a higher microviscosity near TMA-DPH and lower order are predicted for these oxidized membranes, at variance with the usual trend of variation of these two parameters. Finally, the complex relaxation process of Laurdan in pure POPC-OOH membranes also indicates a higher membrane hydration and viscosity in the close vicinity of the -OOH moiety. Altogether, our combined approach reveals that the hydroperoxide group promotes alterations in the membrane structure organization, namely, at the level of membrane order, viscosity, and bending rigidity.

In-vitro prediction of the membranotropic action of emerging organic contaminants using a liposome-based multidisciplinary approach

According to ISO 17402:2008 more knowledge is needed on processes controlling bioavailability of organic species so as to close the still existing gap between chemical measurements and biological effects. The bioavailability concept encompasses the investigation of the degree of penetration of target species across biological membranes. In addition, REACH (Registration, Evaluation, Authorisation and restriction of Chemicals) guidelines promote the use of in-vitro methods against conventional ecotoxicological tests because of the ethical controversy of in-vivo tests. This work is aimed at filling the gap by proposing a multidisciplinary approach based on high-resolution and low-resolution empirical techniques, and theoretical quantum mechanics for the in-vitro investigation of the bioavailability and membranotropic effects of organic emerging contaminants, including bioaccumulation, via passive diffusion across lipid bilayers. Phosphatidylcholine (PC) liposomes are selected as biomembrane surrogates, and contaminant effects are explored by (i) fluorescence anisotropy and generalized polarization assays using membrane fluorescence probes (laurdan and prodan) and UV-Vis spectroscopy, (ii) ¹H NMR measurements to ascertain supramolecular interactions with PC and (iii) molecular dynamics simulations. In particular, un-regulated model compounds with distinct physico-chemical properties that are representative of three different classes of emerging contaminants in environmental compartments are chosen for validation of the holistic approach: (i) diclofenac as a model of anti-inflammatory drug; (ii) triclosan as an anti-microbial agent; and (iii) bisphenol A as a plastic-borne compound, and compared with chlorpyrifos as a legacy insecticide. Laurdan anisotropic measurements are in good agreement with ¹H NMR data and both approaches pinpoint that triclosan and chlorpyrifos are highly bioaccumulative in membranes. Molecular dynamic studies indicate that the lateral diffusion of the lipid bilayer is much lower with the incorporation of either triclosan or chlorpyrifos into the bilayer. The theoretical simulations also allowed estimating absolute bioavailability data under passive diffusion (<0.1%, 63%, 73% and 89% for diclofenac, bisphenol A, triclosan and chlorpyrifos, respectively) given as the percentage of time that a given species is located in the region of the fatty acyl chains. Our findings indicate that PC-based liposome assays serve as a fast and cost-effective in-vitro approach, notwithstanding its low resolution features, for environmental bioavailability studies of emerging contaminants for which insufficient or inconsistent ecotoxicological data are identified in the literature.

Electrostatic Environment of Proteorhodopsin Affects the pKa of Its Buried Primary Proton Acceptor

The protonation state of embedded charged residues in transmembrane proteins (TMPs) can control the onset of protein function. It is understood that interactions between an embedded charged residue and other charged or polar residues in the moiety would influence its pKa, but how the surrounding environment in which the TMP resides affects the pKa of these residues is unclear. Proteorhodopsin (PR), a light-responsive proton pump from marine bacteria, was used as a model to examine externally accessible factors that tune the pKa of its embedded charged residue, specifically its primary proton acceptor D97. The pKa of D97 was compared between PR reconstituted in liposomes with different net headgroup charges and equilibrated in buffer with different ion concentrations. For PR reconstituted in net positively charged compared to net negatively charged liposomes in low-salt buffer solutions, a drop of the apparent pKa from 7.6 to 5.6 was observed, whereas intrinsic pKa modeled with surface pH calculated from Gouy-Chapman predictions found an opposite trend for the pKa change, suggesting that surface pH does not account for the main changes observed in the apparent pKa. This difference in the pKa of D97 observed from PR reconstituted in oppositely charged liposome environments disappeared when the NaCl concentration was increased to 150 mM. We suggest that protein-intrinsic structural properties must play a role in adjusting the local microenvironment around D97 to affect its pKa, as corroborated with observations of changes in protein side-chain and hydration dynamics around the E-F loop of PR. Understanding the effect of externally controllable factors in tuning the pKa of TMP-embedded charged residues is important for bioengineering and biomedical applications relying on TMP systems, in which the onset of functions can be controlled by the protonation state of embedded residues.

Dynamic Protonation Dramatically Affects the Membrane Permeability of Drug-like Molecules

Permeability (P_m) across biological membranes is of fundamental importance and a key factor in drug absorption, distribution, and development. Although the majority of drugs will be charged at some point during oral delivery, our understanding of membrane permeation by charged species is limited. The canonical model assumes that only neutral molecules partition into and passively permeate across membranes, but there is mounting evidence that these processes are also facile for certain charged species. However, it is unknown whether such ionizable permeants dynamically neutralize at the membrane surface or permeate in their charged form. To probe protonation-coupled permeation in atomic detail, we herein apply continuous constant-pH molecular dynamics along with free energy sampling to study the permeation of a weak base propranolol (PPL), and evaluate the impact of including dynamic protonation on P_m. The simulations reveal that PPL dynamically neutralizes at the lipid-tail interface, which dramatically influences the permeation free energy landscape and explains why the conventional model overestimates the assigned intrinsic permeability. We demonstrate how fixed-charge-state simulations can account for this effect, and propose a revised model that better describes pH-coupled partitioning and permeation. Our results demonstrate how dynamic changes in protonation state may play a critical role in the permeation of ionizable molecules, including pharmaceuticals and drug-like molecules, thus requiring a revision of the standard picture.

Interaction of anions with lipid cubic phase membranes, an electrochemical impedance study

HYPOTHESIS

Electrochemical impedance spectroscopy is useful to monitor anionic interactions with a Lipid Cubic Phase, as previously demonstrated for cationic interaction (Khani Meynaq et al., 2016). It was expected that the smaller hydrophilic anions, acetate and chloride, would interact differently than the large tryptophan anion with its hydrophobic tail.

EXPERIMENT

The impedance measurements enabled estimation of resistances and capacitances of a freestanding lipid cubic phase membrane at exposure to 4 and 40 mM solutions of NaCl, NaOAc and NaTrp. Small-angle X-ray scattering was used for cubic phase identification and to track structural changes within the cubic phase when exposed to the different electrolytes.

FINDINGS

The membrane resistance increases at exposure to the electrolytes in the order Cl^- < OAc^- < Trp^-. The membrane resistance decreases with time at exposure to the hydrophilic anions and increases with time at Trp^- exposure. The membrane capacitances were lower for NaTrp compared to NaCl and NaOAc at the corresponding concentrations which is consistent with the results from SAXRD. It is concluded that Trp^- ions do not enter the aqueous channels of the cubic phase but are strongly adsorbed to the membrane/electrolyte interface leading to large alteration of the lipid phase structure and a high membrane resistance.

Micro-pharmacokinetics: Quantifying local drug concentration at live cell membranes

Fundamental equations for determining pharmacological parameters, such as the binding affinity of a ligand for its target receptor, assume a homogeneous distribution of ligand, with concentrations in the immediate vicinity of the receptor being the same as those in the bulk aqueous phase. It is, however, known that drugs are able to interact directly with the plasma membrane, potentially increasing local ligand concentrations around the receptor. We have previously reported an influence of ligand-phospholipid interactions on ligand binding kinetics at the β₂-adrenoceptor, which resulted in distinct "micro-pharmacokinetic" ligand profiles. Here, we directly quantified the local concentration of BODIPY630/650-PEG8-S-propranolol (BY-propranolol), a fluorescent derivative of the classical β-blocker propranolol, at various distances above membranes of single living cells using fluorescence correlation spectroscopy. We show for the first time a significantly increased ligand concentration immediately adjacent to the cell membrane compared to the bulk aqueous phase. We further show a clear role of both the cell membrane and the β₂-adrenoceptor in determining high local BY-propranolol concentrations at the cell surface. These data suggest that the true binding affinity of BY-propranolol for the β₂-adrenoceptor is likely far lower than previously reported and highlights the critical importance of understanding the "micro-pharmacokinetic" profiles of ligands for membrane-associated proteins.

Membrane/Water Partition Coefficients of Bile Salts Determined Using Laurdan as a Fluorescent Probe

The interaction of liposomal membranes composed of soybean phosphatidylcholine with the bile salts (BSs) cholate (Ch), glycocholate (GC), chenodeoxycholate (CDC), and glycochenodeoxycholate (GCDC) was studied. The BSs differed with regard to their lipophilicity, pKa values, and the size of their hydrophilic moiety. Their membrane interactions were investigated using Laurdan as a membrane-anchored fluorescent dye. The apparent membrane/water partition coefficient, D, at pH 7.4 was calculated from binding plots and compared with direct binding measurements using ultracentrifugation as a reference. The Laurdan-derived LogD values at pH 7.4 were found to be 2.10 and 2.25 for the trihydroxy BSs, i.e., Ch and GC, and 2.85 and 2.75 for the dihydroxy BSs, i.e., CDC and GCDC, respectively. For the membrane-associated glycine-conjugated GC and GCDC (pK_a values of ∼3.9), no differences in the Laurdan spectra of the respective BS were found at pH 6.8, 7.4, and 8.2. Unconjugated Ch and CDC (pK_a values of ∼5.0) showed pronounced differences at the three pH values. Furthermore, the kinetics of membrane adsorption and transbilayer movement differed between conjugated and unconjugated BSs as determined with Laurdan-labeled liposomes.

Fluorescent Lipid Nanoparticles as Biomembrane Models for Exploring Emerging Contaminant Bioavailability Supported by Density Functional Theory Calculations

Experimental sensing schemes and thermodynamic in-silico studies are combined holistically in this manuscript so as to give new insights into the bioavailability of environmental contaminants via permeation across lipid nanoparticles (liposomes) as a mimicry of biological membranes. Using Prodan and Laurdan as fluorescent membrane probes, phosphatidylcholine-based unilamellar liposomes are harnessed to investigate membranotropic effects of alkyl esters of p-hydroxybenzoic acid and triclosan in vitro on the basis of steady-state fluorescence anisotropy, light scattering, and generalized polarization measurements. The feasibility of the analytical responses to ascertain differences in temperature-dependent contaminant bioavailability is investigated in detail. High level density functional theory (DFT) calculations (RI-BP86-D3/def2-SVP) have been resorted to investigate noncovalent 1:1 complexes of the fluorescent probes and emerging contaminants with dipalmitoylphosphatidylcholine, as a minimalist model of a lipid nanoparticle, to evaluate both the interaction energies and the geometries of the complexes. This information can be related to the degree of penetration of the guest across the lipid bilayer. Our experimental results supported by in-silico DFT calculations and ecotoxicological data let us to conclude that simple analytical measurements of liposomal changes in lipid packaging, permeability, and fluidity are appropriate to foresee the potential bioavailability and toxicity of emerging contaminants.

In silico pharmacology: Drug membrane partitioning and crossing

Over the past decade, molecular dynamics (MD) simulations have become particularly powerful to rationalize drug insertion and partitioning in lipid bilayers. MD simulations efficiently support experimental evidences, with a comprehensive understanding of molecular interactions driving insertion and crossing. Prediction of drug partitioning is discussed with respect to drug families (anesthetics; β-blockers; non-steroidal anti-inflammatory drugs; antioxidants; antiviral drugs; antimicrobial peptides). To accurately evaluate passive permeation coefficients turned out to be a complex theoretical challenge; however the recent methodological developments based on biased MD simulations are particularly promising. Particular attention is paid to membrane composition (e.g., presence of cholesterol), which influences drug partitioning and permeation. Recent studies concerning in silico models of membrane proteins involved in drug transport (influx and efflux) are also reported here. These studies have allowed gaining insight in drug efflux by, e.g., ABC transporters at an atomic resolution, explicitly accounting for the mandatory forces induced by the surrounded lipid bilayer. Large-scale conformational changes were thoroughly analyzed.

Time-resolved fluorescence in lipid bilayers: selected applications and advantages over steady state

Fluorescence methods are versatile tools for obtaining dynamic and topological information about biomembranes because the molecular interactions taking place in lipid membranes frequently occur on the same timescale as fluorescence emission. The fluorescence intensity decay, in particular, is a powerful reporter of the molecular environment of a fluorophore. The fluorescence lifetime can be sensitive to the local polarity, hydration, viscosity, and/or presence of fluorescence quenchers/energy acceptors within several nanometers of the vicinity of a fluorophore. Illustrative examples of how time-resolved fluorescence measurements can provide more valuable and detailed information about a system than the time-integrated (steady-state) approach will be presented in this review: 1), determination of membrane polarity and mobility using time-dependent spectral shifts; 2), identification of submicroscopic domains by fluorescence lifetime imaging microscopy; 3), elucidation of membrane leakage mechanisms from dye self-quenching assays; and 4), evaluation of nanodomain sizes by time-resolved Förster resonance energy transfer measurements.

Sampling errors in free energy simulations of small molecules in lipid bilayers

Free energy simulations are a powerful tool for evaluating the interactions of molecular solutes with lipid bilayers as mimetics of cellular membranes. However, these simulations are frequently hindered by systematic sampling errors. This review highlights recent progress in computing free energy profiles for inserting molecular solutes into lipid bilayers. Particular emphasis is placed on a systematic analysis of the free energy profiles, identifying the sources of sampling errors that reduce computational efficiency, and highlighting methodological advances that may alleviate sampling deficiencies. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Computer modelling studies of the bilayer/water interface

This review summarises high resolution studies on the interface of lamellar lipid bilayers composed of the most typical lipid molecules which constitute the lipid matrix of biomembranes. The presented results were obtained predominantly by computer modelling methods. Whenever possible, the results were compared with experimental results obtained for similar systems. The first and main section of the review is concerned with the bilayer-water interface and is divided into four subsections. The first describes the simplest case, where the interface consists only of lipid head groups and water molecules and focuses on interactions between the lipid heads and water molecules; the second describes the interface containing also mono- and divalent ions and concentrates on lipid-ion interactions; the third describes direct inter-lipid interactions. These three subsections are followed by a discussion on the network of direct and indirect inter-lipid interactions at the bilayer interface. The second section summarises recent computer simulation studies on the interactions of antibacterial membrane active compounds with various models of the bacterial outer membrane. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Cholesterol, sphingolipids, and glycolipids: what do we know about their role in raft-like membranes?

Lipids rafts are considered to be functional nanoscale membrane domains enriched in cholesterol and sphingolipids, characteristic in particular of the external leaflet of cell membranes. Lipids, together with membrane-associated proteins, are therefore considered to form nanoscale units with potential specific functions. Although the understanding of the structure of rafts in living cells is quite limited, the possible functions of rafts are widely discussed in the literature, highlighting their importance in cellular functions. In this review, we discuss the understanding of rafts that has emerged based on recent atomistic and coarse-grained molecular dynamics simulation studies on the key lipid raft components, which include cholesterol, sphingolipids, glycolipids, and the proteins interacting with these classes of lipids. The simulation results are compared to experiments when possible.