Simulating the effect of alcohol on the structure of a membrane

Advances on nanoformulation approaches for delivering plant-derived antioxidants: A case of quercetin

Oxidative stress has been implicated in tumorigenic, cardiovascular, neuro-, and age-related degenerative changes. Antioxidants minimize the oxidative damage through neutralization of reactive oxygen species (ROS) and other causative agents. Ever since the emergence of COVID-19, plant-derived antioxidants have received enormous attention, particularly in the Indian subcontinent. Quercetin (QCT), a bio-flavonoid, exists in the glycosylated form in fruits, berries and vegetables. The antioxidant potential of QCT analogs relates to the number of free hydroxyl groups in their structure. Despite presence of these groups, QCT exhibits substantial hydrophobicity. Formulation scientists have tested nanotechnology-based approaches for its improved solubilization and delivery to the intended site of action. By the virtue of its hydrophobicity, QCT gets encapsulated in nanocarriers carrying hydrophobic domains. Apart from passive accumulation, active uptake of such formulations into the target cells can be facilitated through well-studied functionalization strategies. In this review, we have discussed the approaches of improving solubilization and bioavailability of QCT with the use of nanoformulations.

Structural properties of cationic surfactant-fatty alcohol bilayers: insights from dissipative particle dynamics

Bilayers, self-assembled by cationic surfactants and fatty alcohols in water, are the basic units of lamellar gel networks - creamy formulations extensively used in cosmetics and pharmaceutics. Mesoscopic modelling and study of the bilayers formed by single- or double-tail cationic surfactants (CTAC or DHDAC), and fatty alcohols (FAs) in the lamellar fluid and gel phases were employed. Fatty alcohols with alkyl tail equal to or greater than the surfactant alkyl tail, i.e., C₁₆FA or C₁₈FA and C₂₂FA, were considered. A model formulation was explored with the FA concentration greater than that of the surfactant and the structure of the fluid and gel bilayers in tensionless state characterised via the density profiles across the bilayers, orientational order parameters of the surfactant and FA chains, intrinsic analysis of the bilayer interfaces, and bending rigidity. The intrinsic analysis allows identification and quantification of the coexistence of the interdigitated and non-interdigitated phases present within the gel bilayers. The FA chains were found to conform the primary scaffolding of the bilayers while the surfactant chains tessellate bilayer monolayers from their water-hydrophobic interface. Further, the overlap of the FA chains from the apposed monolayers of the fluid bilayers rises with increasing FA length. Finally, the prevalence of the non-interdigitated phase over the interdigitated phase within the gel bilayers becomes enhanced upon the FA length increase with a preference of the surfactant chains to reside in the non-interdigitated phase rather than the interdigitated phase.

Design of Lipid-Based Nanocarriers via Cation Modulation of Ethanol-Interdigitated Lipid Membranes

Short-chain alcohols (i.e., ethanol) can induce membrane interdigitation in saturated-chain phosphatidylcholines (PCs). In this process, alcohol molecules intercalate between phosphate heads, increasing lateral separation and favoring hydrophobic interactions between opposing acyl chains, which interpenetrate forming an interdigitated phase. Unraveling mechanisms underlying the interactions between ethanol and model lipid membranes has implications for cell biology, biochemistry, and for the formulation of lipid-based nanocarriers. However, investigations of ethanol-lipid membrane systems have been carried out in deionized water, which limits their applicability. Here, using a combination of small- and wide-angle X-ray scattering, small-angle neutron scattering, and all-atom molecular dynamics simulations, we analyzed the effect of varying CaCl₂ and NaCl concentrations on ethanol-induced interdigitation. We observed that while ethanol addition leads to the interdigitation of bulk phase 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) bilayers in the presence of CaCl₂ and NaCl regardless of the salt concentration, the ethanol-induced interdigitation of vesicular DPPC depends on the choice of cation and its concentration. These findings unravel a key role for cations in the ethanol-induced interdigitation of lipid membranes in either bulk phase or vesicular form.

Increasing Solvent Tolerance to Improve Microbial Production of Alcohols, Terpenoids and Aromatics

Fuels and polymer precursors are widely used in daily life and in many industrial processes. Although these compounds are mainly derived from petrol, bacteria and yeast can produce them in an environment-friendly way. However, these molecules exhibit toxic solvent properties and reduce cell viability of the microbial producer which inevitably impedes high product titers. Hence, studying how product accumulation affects microbes and understanding how microbial adaptive responses counteract these harmful defects helps to maximize yields. Here, we specifically focus on the mode of toxicity of industry-relevant alcohols, terpenoids and aromatics and the associated stress-response mechanisms, encountered in several relevant bacterial and yeast producers. In practice, integrating heterologous defense mechanisms, overexpressing native stress responses or triggering multiple protection pathways by modifying the transcription machinery or small RNAs (sRNAs) are suitable strategies to improve solvent tolerance. Therefore, tolerance engineering, in combination with metabolic pathway optimization, shows high potential in developing superior microbial producers.

Response of microbial membranes to butanol: interdigitation vs. disorder

Biobutanol production by fermentation is potentially a sustainable alternative to butanol production from fossil fuels. However, the toxicity of butanol to fermentative bacteria, resulting largely from cell membrane fluidization, limits production titers and is a major factor limiting the uptake of the technology. Here, studies were undertaken, in vitro and in silico, on the butanol effects on a representative bacterial (i.e. Escherichia coli) inner cell membrane. A critical butanol : lipid ratio for stability of 2 : 1 was observed, computationally, consistent with complete interdigitation. However, at this ratio the bilayer was ∼20% thicker than for full interdigitation. Furthermore, butanol intercalation induced acyl chain bending and increased disorder, measured as a 27% lateral diffusivity increase experimentally in a supported lipid bilayer. There was also a monophasic Tm reduction in butanol-treated large unilamellar vesicles. Both behaviours are inconsistent with an interdigitated gel. Butanol thus causes only partial interdigitation at physiological temperatures, due to butanol accumulating at the phospholipid headgroups. Acyl tail disordering (i.e. splaying and bending) fills the subsequent voids. Finally, butanol short-circuits the bilayer and creates a coupled system where interdigitated and splayed phospholipids coexist. These findings will inform the design of strategies targeting bilayer stability for increasing biobutanol production titers.

Physiological Genomics of Multistress Resistance in the Yeast Cell Model and Factory: Focus on MDR/MXR Transporters

The contemporary approach of physiological genomics is vital in providing the indispensable holistic understanding of the complexity of the molecular targets, signalling pathways and molecular mechanisms underlying the responses and tolerance to stress, a topic of paramount importance in biology and biotechnology. This chapter focuses on the toxicity and tolerance to relevant stresses in the cell factory and eukaryotic model yeast Saccharomyces cerevisiae. Emphasis is given to the function and regulation of multidrug/multixenobiotic resistance (MDR/MXR) transporters. Although these transporters have been considered drug/xenobiotic efflux pumps, the exact mechanism of their involvement in multistress resistance is still open to debate, as highlighted in this chapter. Given the conservation of transport mechanisms from S. cerevisiae to less accessible eukaryotes such as plants, this chapter also provides a proof of concept that validates the relevance of the exploitation of the experimental yeast model to uncover the function of novel MDR/MXR transporters in the plant model Arabidopsis thaliana. This knowledge can be explored for guiding the rational design of more robust yeast strains with improved performance for industrial biotechnology, for overcoming and controlling the deleterious activities of spoiling yeasts in the food industry, for developing efficient strategies to improve crop productivity in agricultural biotechnology.

Partition of common anesthetic molecules in the liquid disordered phase domain of a composite multicomponent membrane

Despite a vast clinical application of anesthetics, the molecular level of understanding of general anesthesia is far from our reach. Using atomistic molecular dynamics simulation, we study the effects of common anesthetics: ethanol, chloroform, and methanol in the fully hydrated symmetric multicomponent lipid bilayer membrane comprised of an unsaturated palmitoyl-oleoyl-phosphatidyl-choline (POPC), a saturated palmitoyl-sphingomyelin, and cholesterol, which exhibits phase coexistence of liquid-ordered (l_{o})-liquid-disordered (l_{d}) phase domains. We find that the mechanical and physical properties such as the thickness and rigidity of the membrane are reduced while the lateral expansion of the membrane is exhibited in the presence of anesthetic molecules. Our simulation shows both lateral and transverse heterogeneity of the anesthetics in the composite multicomponent lipid membrane. Both ethanol and chloroform partition in the POPC-rich l_{d} phase domain, while methanol is distributed in both l_{o}-l_{d} phase domains. Chloroform can penetrate deep into the membrane, while methanol partitions mostly at the water layer closed to the head group and ethanol at the neck of the lipids in the membrane.

Farnesol-Induced Disruption of the Staphylococcus aureus Cytoplasmic Membrane

Farnesol, a sesquiterpene alcohol with an aliphatic carbon chain, inhibited the growth of Staphylococcus aureus and induced the leakage of potassium ions. We investigated the action of farnesol on the cytoplasmic membrane of S. aureus. No ion channels that would account for the loss of potassium ions were detected. Electron paramagnetic resonance measurements suggested that farnesol proceeds into the cytoplasmic membrane of S. aureus cells, where it induces the disordering and eventual disruption of the cytoplasmic membrane. This was supported by the result that the effects of farnesol decreased by the addition of carotenoid which was the stabilizing reagent for the lipid bilayer.

Structural changes of a sodium dodecyl sulfate (SDS) micelle induced by alcohol molecules

Coarse-grained dynamical simulations have been performed to investigate the behavior of a surfactant micelle in the presence of six different alcohols: hexanol, octanol, decanol, dodecanol, tetradecanol, and hexadecanol. The self-assembly of sodium dodecyl sulfate (SDS) is modified by the alcohol molecules into cylindrical and bilayer micelles as a function of the alcohol/SDS mass ratio. Therefore, in order to understand, from a molecular point of view, how SDS and alcohol molecules self-organize to form the new micelles, different studies were carried out. Analysis of micelle structures, density profiles, and parameters of order were conducted to characterize the shape and size of those micelles. The density profiles revealed that the alcohol molecules were located at the water-micelle interface next to the SDS molecules at low alcohol/SDS mass ratio. At high alcohol/SDS mass ratios, alcohol molecules moved to the middle of the micelle by increasing their size and by producing a structural change. Moreover, micelle structures and sizes were influenced not only by the alcohol/SDS mass ratio but also by the order of the SDS and alcohol tails. Finally, the size of the micelles and enthalpy calculations were used as order parameters to determine a structural phase diagram of alcohol/SDS mixtures in water. Graphical Abstract Structural transition of SDS/alcohol mixtures.

Ethanol induces the formation of water-permeable defects in model bilayers of skin lipids

We observe that ethanol can induce the formation of water-permeable defects in model bilayers of skin lipids and propose this as a new mechanism of action of ethanol as a membrane modulator.

Lipid diffusion in alcoholic environment

We have studied the effects of a high concentration of butanol and octanol on the phase behavior and on the lateral mobility of 1,2-palmitoyl-sn-glycero-3-phosphocholine (DPPC) by means of differential scanning calorimetry and pulsed-gradient stimulated-echo (PGSTE) NMR spectroscopy. A lowering of the lipid transition from the gel to the liquid-crystalline state for the membrane-alcohol systems has been observed. NMR measurements reveal three distinct diffusions in the DPPC-alcohol systems, characterized by a high, intermediate, and slow diffusivity, ascribed to the water, the alcohol, and the lipid, respectively. The lipid diffusion process is promoted in the liquid phase while it is hindered in the interdigitated phase due to the presence of alcohols. Furthermore, in the interdigitated phase, lipid lateral diffusion coefficients show a slight temperature dependence. To the best of our knowledge, this is the first time that lateral diffusion coefficients on alcohol with so a long chain, and at low temperatures, are reported. By the Arrhenius plots of the temperature dependence of the diffusion coefficients, we have evaluated the apparent activation energy in both the liquid and in the interdigitated phase. The presence of alcohol increases this value in both phases. An explanation in terms of a free volume model that takes into account also for energy factors is proposed.

Chronic alcoholism-mediated impairment in the medulla oblongata: a mechanism of alcohol-related mortality in traumatic brain injury?

Alcohol-related traumatic brain injury (TBI) is a common condition in medical and forensic practice, and results in high prehospital mortality. We investigated the mechanism of chronic alcoholism-related mortality by examining the effects of alcohol on the synapses of the medulla oblongata in a rat model of TBI. Seventy adult male Sprague-Dawley rats were randomly assigned to either ethanol (EtOH) group, EtOH-TBI group, or control groups (water group, water-TBI group). To establish chronic alcoholism model, rats in the EtOH group were given EtOH twice daily (4 g/kg for 2 weeks and 6 g/kg for another 2 weeks). The rats also received a minor strike on the occipital tuberosity with an iron pendulum. Histopathologic and ultrastructure changes and the numerical density of the synapses in the medulla oblongata were examined. Expression of postsynaptic density-95 (PSD-95) in the medulla oblongata was measured by ELISA. Compared with rats in the control group, rats in the chronic alcoholism group showed: (1) minor axonal degeneration; (2) a significant decrease in the numerical density of synapses (p < 0.01); and (3) compensatory increase in PSD-95 expression (p < 0.01). Rats in the EtOH-TBI group showed: (1) high mortality (50%, p < 0.01); (2) inhibited respiration before death; (3) severe axonal injury; and (4) decrease in PSD-95 expression (p < 0.05). Chronic alcoholism induces significant synapse loss and axonal impairment in the medulla oblongata and renders the brain more susceptible to TBI. The combined effects of chronic alcoholism and TBI induce significant synapse and axon impairment and result in high mortality.

Investigations on the melting and bending modulus of polymer grafted bilayers using dissipative particle dynamics

Understanding the influence of polymer grafted bilayers on the physicomechanical properties of lipid membranes is important while developing liposomal based drug delivery systems. The melting characteristics and bending moduli of polymer grafted bilayers are investigated using dissipative particle dynamics simulations as a function of the amount of grafted polymer and lipid tail length. Simulations are carried out using a modified Andersen barostat, whereby the membrane is maintained in a tensionless state. For lipids made up of four to six tail beads, the transition from the low temperature L(β) phase to the L(α) phase is lowered only above a grafting fraction of G(f)=0.12 for polymers made up of 20 beads. Below G(f)=0.12 small changes are observed only for the HT(4) bilayer. The bending modulus of the bilayers is obtained as a function of G(f) from a Fourier analysis of the height fluctuations. Using the theory developed by Marsh et al. [Biochim. Biophys. Acta 1615, 33 (2003)] for polymer grafted membranes, the contributions to the bending modulus due to changes arising from the grafted polymer and bilayer thinning are partitioned. The contributions to the changes in κ from bilayer thinning were found to lie within 11% for the lipids with four to six tail beads, increasing to 15% for the lipids containing nine tail beads. The changes in the area stretch modulus were also assessed and were found to have a small influence on the overall contribution from membrane thinning. The increase in the area per head group of the lipids was found to be consistent with the scalings predicted by self-consistent mean field results.

Effect of polymer grafting on the bilayer gel to liquid-crystalline transition

Grafted polymers on the surface of lipid membranes have potential applications in liposome-based drug delivery and supported membrane systems. The effect of polymer grafting on the phase behavior of bilayers made up of single-tail lipids is investigated using dissipative particle dynamics. The bilayer is maintained in a tensionless state using a barostat. Simulations are carried out by varying the grafting fraction, G(f), defined as the ratio of the number of polymer molecules to the number of lipid molecules, and the length of the lipid tails. At low G(f), the bilayer shows a sharp transition from the gel (L(beta)) to the liquid-crystalline (L(alpha)) phase. This main melting transition temperature is lowered as G(f) is increased, and above a critical value of G(f), the interdigitated L(betaI) phase is observed prior to the main transition. The temperature range over which the intermediate phases are observed is a function of the lipid tail length and G(f). At higher grafting fractions, the presence of the L(betaI) phase is attributed to the increase in the area per head group due to the lateral pressure exerted by the polymer brush. The areal expansion and decrease in the melting temperatures as a function of G(f) were found to follow the scalings predicted by the self-consistent mean field theories for grafted polymer membranes. Our study shows that the grafted polymer density can be used to effectively control the temperature range and occurrence of a given bilayer phase.

Using ergosterol to mitigate the deleterious effects of ethanol on bilayer structure

In wine fermentations, yeast is exposed to concentrated ethanol solutions. Ergosterol, a sterol that is found in lower eukaryotic membranes, helps preserve the structural integrity of yeast membranes in stressful environmental conditions. A premature arrest in ethanol production due to unknown metabolic changes in yeasts results in undesirably large concentrations of residual sugar and may be caused by the formation of an ethanol-induced interdigitated phase. We use atomistic molecular dynamics simulations to examine the induction of the interdigitated phase in model yeast membranes that contain either 0, 10, 20, 25 mol % ergosterol in ethanol concentrations of 0, 10, 15 vol %. The 25 mol % ergosterol system shows a similar level of interdigitation for the 0 and 10 vol % ethanol solutions, indicating that ergosterol molecules in this system are able to effectively counteract the disruptive behavior of ethanol molecules. However, at a 15 vol % ethanol solution, the amount of interdigitation triples and this ethanol concentration is similar to the concentrations found in stuck fermentations. The other three ergosterol concentrations studied (0, 10, 20 mol %) show larger quantities of interdigitation in the 10 vol % ethanol solution than the 0 vol % solution. Thus, the 25 mol % ergosterol bilayer, which is representative of the ergosterol concentrations seen in yeast membranes, is unique in the systems examined in its ability to delay the onset of ethanol-induced interdigitation. The concentration of ergosterol affects the permeability of a fluid-phase bilayer, where the 10 mol % ergosterol bilayer is more permeable to ethanol than either a bilayer containing no ergosterol molecules or a bilayer containing 20/25 mol % ergosterol. This lipid permeability appears to be correlated with the existence of a lipid region whose lipids neither have direct contact with ergosterol molecules nor exhibit bulk lipid/lipid interactions.

Effect of cholesterol on the structure of a phospholipid bilayer

Cholesterol plays an important role in regulating the properties of phospholipid membranes. To obtain a detailed understanding of the lipid-cholesterol interactions, we have developed a mesoscopic water-lipid-cholesterol model. In this model, we take into account the hydrophobic-hydrophilic interactions and the structure of the molecules. We compute the phase diagram of dimyristoylphosphatidylcholine-cholesterol by using dissipative particle dynamics and show that our model predicts many of the different phases that have been observed experimentally. In quantitative agreement with experimental data our model also shows the condensation effect; upon the addition of cholesterol, the area per lipid decreases more than one would expect from ideal mixing. Our calculations show that this effect is maximal close to the main-phase transition temperature, the lowest temperature for which the membrane is in the liquid phase, and is directly related to the increase of this main-phase transition temperature upon addition of cholesterol. We demonstrate that no condensation is observed if we slightly change the structure of the cholesterol molecule by adding an extra hydrophilic head group or if we decrease the size of the hydrophobic part of cholesterol.

The influence of 1-alkanols and external pressure on the lateral pressure profiles of lipid bilayers

The suggestion by Robert Cantor, that drug-induced pressure changes in lipid bilayers can change the conformational equilibrium between open and closed states of membrane proteins and thereby cause anesthesia, attracted much attention lately. Here, we studied the effect of both large external pressure and of 1-alkanols of different chain lengths--some of them anesthetics, others not--on the lateral pressure profiles across dimyristoylphosphatidylcholine (DMPC) bilayers by molecular dynamics simulations. For a pure DMPC bilayer, high pressure both reduced and broadened the tension at the interface hydrophobic/hydrophilic and diminished the repulsion between the phospholipid headgroups. Whereas the effect of ethanol on the lateral pressure profile was similar to the effect of a large external pressure on a DMPC bilayer, long-chain 1-alkanols significantly amplified local maxima and minima in the lateral pressure profile. For most 1-alkanols, external pressure had moderate effects and did not reverse the changes 1-alkanols exerted on the pressure profile. Nevertheless, assuming the bent helix model as a simple geometric model for the transmembrane region of a membrane protein, protein conformational equilibria were shifted in opposite directions by addition of 1-alkanols and additional application of external pressure.

1-Alkanols and membranes: a story of attraction

Although 1-alkanols have long been known to act as penetration enhancers and anesthetics, the mode of operation is not yet understood. In this study, long-time molecular dynamics simulations have been performed to investigate the effect of 1-alkanols of various carbon chain lengths onto the structure and dynamics of dimyristoylphosphatidylcholine bilayers. The simulations were complemented by microcalorimetry, continuous bleaching and film balance experiments. In the simulations, all investigated 1-alkanols assembled inside the lipid bilayer within tens of nanoseconds. Their hydroxyl groups bound preferentially to the lipid carbonyl group and the hydrocarbon chains stretched into the hydrophobic core of the bilayer. Both molecular dynamics simulations and experiments showed that all 1-alkanols drastically affected the bilayer properties. Insertion of long-chain 1-alkanols decreased the area per lipid while increasing the thickness of the bilayer and the order of the lipids. The bilayer elasticity was reduced and the diffusive motion of the lipids within the bilayer plane was suppressed. On the other hand, integration of ethanol into the bilayer enlarged the area per lipid. The bilayer became softer and lipid diffusion was enhanced.

Phase behavior of model lipid bilayers

We investigated the phase behavior of double-tail lipids, as a function of temperature, headgroup interaction and tail length. At low values of the head-head repulsion parameter a(hh), the bilayer undergoes with increasing temperature the transitions from the subgel phase L(c) via the flat gel phase L(beta) to the fluid phase L(alpha). For higher values of a(hh), the transition from the L(c) to the L(alpha) phase occurs via the tilted gel phase L(beta)(') and the rippled phase P(beta)('). The occurrence of the L(beta)(') phase depends on tail length. We find that the rippled structure (P(beta)(')) occurs if the headgroups are sufficiently surrounded by water and that the ripple is a coexistence between the L(c) or L(beta)(') phase and the L(alpha) phase. The anomalous swelling, observed at the P(beta)(') --> L(alpha) transition, is not directly related to the rippled phase, but a consequence of conformational changes of the tails.

How alcohol chain-length and concentration modulate hydrogen bond formation in a lipid bilayer

Molecular dynamics simulations are used to measure the change in properties of a hydrated dipalmitoylphosphatidylcholine bilayer when solvated with ethanol, propanol, and butanol solutions. There are eight oxygen atoms in dipalmitoylphosphatidylcholine that serve as hydrogen bond acceptors, and two of the oxygen atoms participate in hydrogen bonds that exist for significantly longer time spans than the hydrogen bonds at the other six oxygen atoms for the ethanol and propanol simulations. We conclude that this is caused by the lipid head group conformation, where the two favored hydrogen-bonding sites are partially protected between the head group choline and the sn-2 carbonyl oxygen. We find that the concentration of the alcohol in the ethanol and propanol simulations does not have a significant influence on the locations of the alcohol/lipid hydrogen bonds, whereas the concentration does impact the locations of the butanol/lipid hydrogen bonds. The concentration is important for all three alcohol types when the lipid chain order is examined, where, with the exception of the high-concentration butanol simulation, the alcohol molecules having the longest hydrogen-bonding relaxation times at the favored carbonyl oxygen acceptor sites also have the largest order in the upper chain region. The lipid behavior in the high-concentration butanol simulation differs significantly from that of the other alcohol concentrations in the order parameter, head group rotational relaxation time, and alcohol/lipid hydrogen-bonding location and relaxation time. This appears to be the result of the system being very near to a phase transition, and one occurrence of lipid flip-flop is seen at this concentration.

Artifacts in dynamical simulations of coarse-grained model lipid bilayers

With special focus on dissipative particle dynamics simulations of anisotropic and complex soft matter, such as lipid bilayers in water, we have investigated the occurrence of artifacts in the results obtained from dynamical simulations of coarse-grained particle-based models. The particles are modeled by beads that interact via soft repulsive conservative forces (as defined in dissipative particle dynamics simulations), harmonic bond potentials, as well as bending potentials imparting stiffness to the lipid tails. Two different update schemes are investigated: dissipative particle dynamics with a velocity-Verlet-like integration scheme [G. Besold, I. Vattulainen, M. Karttunen, and J. M. Polson, Phys. Rev. E 63, R7611 (2000)] and Lowe-Andersen thermostatting [C. P. Lowe, Europhys. Lett. 47, 145 (1999)] with the standard velocity-Verlet integration algorithm. By varying the integration time step, we examine various physical quantities, in particular pressure profiles and kinetic bead temperatures, for their sensitivity to artifacts caused by the specific combination of integration technique and the thermostat. We then propose a simple fingerprint method that allows monitoring the presence of simulation artifacts.

Alcohols reduce lateral membrane pressures: predictions from molecular theory

We explore the effects of alcohols on fluid lipid bilayers using a molecular theory with a coarse-grained model. We show that the trends predicted from the theory in the changes in area per lipid, alcohol concentration in the bilayer, and area compressibility modulus, as a function of alcohol chain length and of the alcohol concentration in the solvent far from the bilayer, follow those found experimentally. We then use the theory to study the effect of added alcohol on the lateral pressure profile across the membrane, and find that added alcohol reduces the surface tensions at both the headgroup/solvent and headgroup/tailgroup interfaces, as well as the lateral pressures in the headgroup and tailgroup regions. These changes in lateral pressures could affect the conformations of membrane proteins, providing a nonspecific mechanism for the biological effects of alcohols on cells.

Modelling of proteins in membranes

This review describes some recent theories and simulations of mesoscopic and microscopic models of lipid membranes with embedded or attached proteins. We summarize results supporting our understanding of phenomena for which the activities of proteins in membranes are expected to be significantly affected by the lipid environment. Theoretical predictions are pointed out, and compared to experimental findings, if available. Among others, the following phenomena are discussed: interactions of interfacially adsorbed peptides, pore-forming amphipathic peptides, adsorption of charged proteins onto oppositely charged lipid membranes, lipid-induced tilting of proteins embedded in lipid bilayers, protein-induced bilayer deformations, protein insertion and assembly, and lipid-controlled functioning of membrane proteins.

Can purely repulsive soft potentials predict micelle formation correctly?

We determine the free energy of micelle formation for surfactants in a solvent by employing a hybrid semi-grand Monte Carlo simulation scheme in combination with umbrella sampling and configurational bias techniques. We compare results of two surfactant models: one based on Lennard-Jones interactions and one based on the soft repulsive potential that is often used in dissipative particle dynamics (DPD). The free energies of micellization in both models show similar behavior. However, although the critical micelle concentration for the Lennard-Jones systems lies within the experimental range, it is 13 orders of magnitude lower for the soft repulsive model. We discuss the implication for the applicability of soft repulsive potentials for the study of micelle formation.

Configurational Temperature in Membrane Simulations Using Dissipative Particle Dynamics

The use of excessively long time steps in dissipative particle dynamics simulations may produce simulation artifacts due to the generation of configurations which are not representative of the desired canonical ensemble. The configurational temperature, among other quantities, may be used to assess the extent of the deviation from equilibrium. This paper presents results for simulations of models of water and lipid bilayer membranes to illustrate the nature of the problems.

Simulation studies of protein-induced bilayer deformations, and lipid-induced protein tilting, on a mesoscopic model for lipid bilayers with embedded proteins

Biological membranes are complex and highly cooperative structures. To relate biomembrane structure to their biological function it is often necessary to consider simpler systems. Lipid bilayers composed of one or two lipid species, and with embedded proteins, provide a model system for biological membranes. Here we present a mesoscopic model for lipid bilayers with embedded proteins, which we have studied with the help of the dissipative particle dynamics simulation technique. Because hydrophobic matching is believed to be one of the main physical mechanisms regulating lipid-protein interactions in membranes, we considered proteins of different hydrophobic length (as well as different sizes). We studied the cooperative behavior of the lipid-protein system at mesoscopic time- and lengthscales. In particular, we correlated in a systematic way the protein-induced bilayer perturbation, and the lipid-induced protein tilt, with the hydrophobic mismatch (positive and negative) between the protein hydrophobic length and the pure lipid bilayer hydrophobic thickness. The protein-induced bilayer perturbation was quantified in terms of a coherence length, xi(P), of the lipid bilayer hydrophobic thickness profile around the protein. The dependence on temperature of xi(P), and the protein tilt-angle, were studied above the main-transition temperature of the pure system, i.e., in the fluid phase. We found that xi(P) depends on mismatch, i.e., the higher the mismatch is, the longer xi(P) becomes, at least for positive values of mismatch; a dependence on the protein size appears as well. In the case of large model proteins experiencing extreme mismatch conditions, in the region next to the so-called lipid annulus, there appears an undershooting (or overshooting) region where the bilayer hydrophobic thickness is locally lower (or higher) than in the unperturbed bilayer, depending on whether the protein hydrophobic length is longer (or shorter) than the pure lipid bilayer hydrophobic thickness. Proteins may tilt when embedded in a too-thin bilayer. Our simulation data suggest that, when the embedded protein has a small size, the main mechanism to compensate for a large hydrophobic mismatch is the tilt, whereas large proteins react to negative mismatch by causing an increase of the hydrophobic thickness of the nearby bilayer. Furthermore, for the case of small, peptidelike proteins, we found the same type of functional dependence of the protein tilt-angle on mismatch, as was recently detected by fluorescence spectroscopy measurements.