Molecular study of the diffusional process of DMSO in double lipid bilayers

As a way to quantify the diffusion process of molecular compounds through biological membranes, we investigated in this study the dynamics of DMSO through an 1,2-Dipalmitoyl-sn-Glycero-3-Phosphocholine (DPPC) bilayer system. To properly account for the diffusion of DMSO due to a concentration gradient, a double DPPC bilayer was setup for our simulations. In such configuration, the aqueous phases can be explicitly associated with the extra and intracellular domains of the membrane, which is seldom the case in studies of single lipid bilayer due to the periodicity imposed by the simulations. DMSO molecules were initially contained in one of the aqueous phases (extracellular region) at a concentration of 5 wt.%. Molecular dynamics simulation was performed in this system for 95 ns at 350 K and 1 bar. The simulations showed that although many DMSO molecules penetrated the lipid bilayer, only about 10% of them crossed the bilayer to reach the other aqueous phase corresponding to the intracellular region of the membrane. The simulation time considered was insufficient to reach equilibrium of the DMSO concentration between the aqueous phases. However, the simulations provided sufficient information to estimate parameters to apply Fick's Law to model the diffusion process of the system. Using this model, we predicted that for the time considered in our simulation, the concentration of DMSO in the intracellular domain should have been about half of the actual value obtained. The model also predicted that equilibrium of the DMSO concentration in the system would be reached after about 2000 ns, approximately 20 times longer than the performed simulation.     

Molecular dynamic simulation study of cholesterol and conjugated double bonds in lipid bilayers

Conjugated linoleic acids (CLA) are found naturally in dairy products. Two isomers of CLA, that differ only in the location of cis and trans double bonds, are found to have distinct and different biological effects. The cis 9 trans 11 (C9T11) isomer is believed to have anti-carcinogenic effects, while the trans 10 cis 12 (T10C12) isomer is believed to be associated with anti-obesity effects. In this paper we extend earlier molecular dynamics (MD) simulations of pure CLA–phosphatidylcholine bilayers to investigate the comparative effects of cholesterol on bilayers composed of the two respective isomers. Simulations of phosphatidylcholine lipid bilayers in which the sn-2 chains contained one of the two isomers of CLA were performed in which, for each isomer, the simulated bilayers contained 10% and 30% cholesterol (Chol). From MD trajectories we calculate and compare structural properties of the bilayers, including areas per molecule, thickness of bilayers, tilt angle of cholesterols, order parameter profiles, and one and two-dimensional radial distribution function (RDF), as functions of Chol concentration. While the structural effect of cholesterol is approximately the same for both isomers, we find differences at an atomistic level in order parameter profiles and in two-dimensional radial distribution functions.     

Conjugated double bonds in lipid bilayers: a molecular dynamics simulation study

Conjugated linoleic acids (CLA) are found naturally in dairy products. Two isomers of CLA, that differ only in the location of cis and trans double bonds, are found to have distinct and different biological effects. The cis 9 trans 11 (C9T11) isomer is attributed to have the anti-carcinogenic effects, while the trans 10 cis 12 (T10C12) isomer is believed to be responsible for the anti-obesity effects. Since dietary CLA are incorporated into membrane phospholipids, we have used Molecular Dynamics (MD) simulations to investigate the comparative effects of the two isomers on lipid bilayer structure. Specifically, simulations of phosphatidylcholine lipid bilayers in which the sn-2 chains contained one of the two isomers of CLA were performed. Force field parameters for the torsional potential of double bonds were obtained from ab initio calculations. From the MD trajectories we calculated and compared structural properties of the two lipid bilayers, including areas per molecule, density profiles, thickness of bilayers, tilt angle of tail chains, order parameters profiles, radial distribution function (RDF) and lateral pressure profiles. The main differences found between bilayers of the two CLA isomers, are (1) the order parameter profile for C9T11 has a dip in the middle of sn-2 chain while the profile for T10C12 has a deeper dip close to terminal of sn-2 chain, and (2) the lateral pressure profiles show differences between the two isomers. Our simulation results reveal localized physical structural differences between bilayers of the two CLA isomers that may contribute to different biological effects through differential interactions with membrane proteins or cholesterol.     

Characterization of domain instabilities in lipid bilayers by Karhunen-Loeve analysis

We study the stability of lipid bilayers with artificial domains. In investigating different domain structures, we identify scenarios of stable and unstable arrangements of patches of mixed phospholipids. These are then characterized using Karhunen-Loeve Expansion (KLE), a special form of Principal Components Analysis (PCA). The simulation data are interrogated using KLE to reveal spatiotemporal patterns that explain relevant motions in the bilayer system. By projecting the high-dimensional dataset onto a small number of key modes, KLE reveals specific dynamic signatures that can help distinguish and characterize various domain instability mechanisms. We find that typically very few modes are responsible for describing a mechanism of instability to a reasonable extent and can clearly distinguish between stable and unstable arrangements. Different instability modes are characterized as they exhibit unique features like global deformation or local mixing modes.     

The modulation of thermal properties of vinblastine by cholesterol in membrane bilayers

It has been shown that the partitioning of vinblastine in 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine (DPPC) single and multiple bilayer dispersions induces partial interdigitation of the lipid alkyl chains. Similar behavior has been observed for abietic and ursodeoxycholic acids and may well be generalized for the partitioning of bulky amphoteric molecules, which tend to localize in the vicinity of the polar heads. For the present study, differential scanning calorimetry (DSC) has been employed to investigate the role of lipid molecular characteristics such as the alkyl chain length and the polarity of the head-group, as well as the impact of cholesterol upon vinblastine-induced interdigitation. It is found that vinblastine does not induce interdigitation in lipids with either shorter or longer alkyl chains than DPPC, or having head-groups of different polarity. In addition, it is shown that the presence of cholesterol in the lipid bilayer tends to modulate the phase behavior of the lipid/vinblastine bilayer system. Preliminary studies show that such properties directly affect the encapsulation efficiency and the pharmacokinetics of liposomes.     

The resealing process of lipid bilayers after reversible electrical breakdown

The resealing process of lipid bilayer membranes after reversible electrical breakdown was investigated using two voltage pulses switched on together. Electrical breakdown of the membranes was induced with a voltage pulse of high intensity and short duration. The time course of the change in membrane conductance after the application of the high (short) voltage pulse was measured with a longer voltage pulse of low amplitude. The decrease in membrane conductance during the resealing process could be fitted to a single exponential curve with a time constant of 10-2 μs in the temperature range between 2 and 20°C. The activation energy for this exponential decay process was found to be about 50 kJ/mol, which might indicate a diffusion process. Above 25°C the resealing process is controlled by two exponential processes.The data obtained for the time course of the resealing process can be explained in terms of pore formation in the membranes in response to the high electrical field strength. A radius of about 4 nm is calculated for the initial pore size. From the assumed exponential change of the pore area with progressive resealing time a diffusion constant of 10^?8 cm²/s for lateral lipid diffusion can be estimated.     

Temperature dependence of water self-diffusion through lipid bilayers assessed by NMR

The temperature dependence of the coefficient of water self-diffusion across plane-parallel multib-ilayers of dioleoylphosphatidylcholine oriented on a glass support was studied in the 20–60°C range by pulsed field gradient NMR. The coefficient for transbilayer diffusion of water proved almost four orders of magnitude smaller than for bulk water, and 10 times smaller than that for lateral diffusion of lipid under the same conditions. The temperature dependence obeyed the Arrhenius law with apparent activation energy of 41 kJ/mol, much higher than that for bulk water (18 kJ/mol). The experimental data were analyzed using the “dissolution-diffusion” model, by simulating water passage through membrane channels, and by examining water exchange in states with different modes of translational mobility, including pore channels and bilayer defects. Each approach could take into account the role of bilayer permeability and assess the apparent activation energy for water diffusion in the hydrophobic part of the bilayer, which proved close to the value for bulk water. Estimates were obtained for water diffusion coefficients in the system, coefficients of bilayer permeability for water, and the influence of bilayer defects on the lateral and transverse diffusion coefficients.     

Molecular dynamics simulation of dipalmitoylphosphatidylcholine modified with a MTSL nitroxide spin label in a lipid membrane

We investigate the interaction between dipalmitoylphosphatidylcholine (DPPC) and a nitroxide spin label in order to understand its influences on lipid structure and dynamics using molecular dynamics simulations. The system was modified by covalently attaching nitroxide spin labels to the headgroups of two DPPC molecules. (S-(2,2,5,5-tetramethyl-2,5-dihydro-1H-pyrrol-3-yl)methyl methanesulfonothioate) (MTSL) was used as the spin label. The label position and dynamics were analyzed as was the impact of the modified DPPC on the structure of the surrounding lipids. The modified DPPC molecules locate closer to the center of the membrane than unmodified DPPC molecules. The rotation of the spin label is unrestricted, but there are favored orientations. MTSL depresses the deuterium order parameters of the carbon atoms close to the headgroup in surrounding DPPC molecules. The spin label has no impact on order parameters of carbon atoms at the end of the lipid tails. The lateral diffusion constant of the modified DPPC is indistinguishable from unmodified DPPC molecules. These novel computational results suggest an experimental validation.     

Outer membrane proteins: comparing X-ray and NMR structures by MD simulations in lipid bilayers

The structures of three bacterial outer membrane proteins (OmpA, OmpX and PagP) have been determined by both X-ray diffraction and NMR. We have used multiple (7 × 15 ns) MD simulations to compare the conformational dynamics resulting from the X-ray versus the NMR structures, each protein being simulated in a lipid (DMPC) bilayer. Conformational drift was assessed via calculation of the root mean square deviation as a function of time. On this basis the ‘quality’ of the starting structure seems mainly to influence the simulation stability of the transmembrane β-barrel domain. Root mean square fluctuations were used to compare simulation mobility as a function of residue number. The resultant residue mobility profiles were qualitatively similar for the corresponding X-ray and NMR structure-based simulations. However, all three proteins were generally more mobile in the NMR-based than in the X-ray simulations. Principal components analysis was used to identify the dominant motions within each simulation. The first two eigenvectors (which account for >50% of the protein motion) reveal that such motions are concentrated in the extracellular loops and, in the case of PagP, in the N-terminal α-helix. Residue profiles of the magnitude of motions corresponding to the first two eigenvectors are similar for the corresponding X-ray and NMR simulations, but the directions of these motions correlate poorly reflecting incomplete sampling on a ∼10 ns timescale.     

Capturing the nanoscale complexity of cellular membranes in supported lipid bilayers

The lateral mobility of cell membranes plays an important role in cell signaling, governing the rate at which embedded proteins can interact with other biomolecules. The past two decades have seen a dramatic transformation in understanding of this environment, as the mechanisms and potential implications of nanoscale structure of these systems has become accessible to theoretical and experimental investigation. In particular, emerging micro- and nano-scale fabrication techniques have made possible the direct manipulation of model membranes at the scales relevant to these biological processes. This review focuses on recent advances in nanopatterning of supported lipid bilayers, capturing the impact of membrane nanostructure on molecular diffusion and providing a powerful platform for further investigation of the role of this spatial complexity on cell signaling.     

A low-temperature structural phase transition of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine bilayers in the gel phase

A new thermotropic phase transition, at ?30°C and atmospheric pressure, was found to occur in the gel phase of aqueous DPPC dispersions. The Raman spectral changes at this phase transition are similar to those observed in the gel phase of DMPC dispersions at ?60°C. The thermotropic phase transition at ?30°C is equivalent to the barotropic GII to GIII phase transition observed in DPPC at 1.7 kbar and 30°C. It is shown that the rate of the large angle interchain reorientational fluctuations decreases gradually with decreasing temperature, and that the orientationally disordered acyl chain structure of the GII phase is extended into the GIII phase of DPPC. The interchain interaction, arising from the damping of the reorientational fluctuations, increases with decreasing temperature in the GII gel phase as well as in the GIII gel phase.     

Effect of retinol and retinoic acid on permeability,electrical resistance and phase transition of lipid bilayers

Retinol and retinoic acid have been incorporated into the artificial membrane systems, planar bimolecular lipid membranes and liposomes, and their effects on several membrane parameters have been measured. 1. Retinol and retinoic acid increased the permeability of egg lecithin liposomes to K⁺, I^? and glucose when incorporated into the membranes at levels as low as 0.5 membrane mol%. Retinoic acid influenced permeability more than did retinol for each of the solutes tested. 2. Retinol and retinoic acid both decreased the electrical resistance of egg lecithin-planar bimolecular lipid membranes from 0.5 to 8 membrane mol%. Retinoic acid effected a larger change than did retinol. 3. Retinol and retinoic acid increased the permeability of dimyristoylphosphatidylcholine and dipalmitoylphosphatidylcholine liposomes to water at 1.0 and 3.0 membrane mol%. A larger effect on water permeability was measured for retinoic acid than for retinol. 4. Retinol and retinoic acid at 1.0 and 3.0 membrane mol% were shown to lower the phase-transition temperature of liposomes composed of dimyristoylphosphatidylcholine or dipalmitoylphosphatidylcholine. Phase-transition temperatures were monitored by abrupt changes in water permeability and liposome size associated with the transition. Retinoic acid lowered the phase-transition temperature of dimyristoylphosphatidylcholine liposomes more than did retinol, while both retinoids had almost the same effect on dipalmitoylphosphatidylcholine liposomes.     

Simulation of gel phase formation and melting in lipid bilayers using a coarse grained model

The transformation between a gel and a fluid phase in dipalmitoyl-phosphatidylcholine (DPPC) bilayers has been simulated using a coarse grained (CG) model by cooling bilayer patches composed of up to 8000 lipids. The critical step in the transformation process is the nucleation of a gel cluster consisting of 20-80 lipids, spanning both monolayers. After the formation of the critical cluster, a fast growth regime is entered. Growth slows when multiple gel domains start interacting, forming a percolating network. Long-lived fluid domains remain trapped and can be metastable on a microsecond time scale. From the temperature dependence of the rate of cluster growth, the line tension of the fluid-gel interface was estimated to be 3+/-2 pN. The reverse process is observed when heating the gel phase. No evidence is found for a hexatic phase as an intermediate stage of melting. The hysteresis observed in the freezing and melting transformation is found to depend both on the system size and on the time scale of the simulation. Extrapolating to macroscopic length and time scales, the transition temperature for heating and cooling converges to 295+/-5 K, in semi-quantitative agreement with the experimental value for DPPC (315 K). The phase transformation is associated with a drop in lateral mobility of the lipids by two orders of magnitude, and an increase in the rotational correlation time of the same order of magnitude. The lipid headgroups, however, remain fluid. These observations are in agreement with experimental findings, and show that the nature of the ordered phase obtained with the CG model is indeed a gel rather than a crystalline phase. Simulations performed at different levels of hydration furthermore show that the gel phase is stabilized at low hydration. A simulation of a small DPPC vesicle reveals that curvature has the opposite effect.     

Studies on the lack of cooperativity in the melting of lipid bilayers

The gel-to-fluid first-order melting transition of lipid bilayers is simulated by the use of a microscopic interaction model which includes a variable number of lipid-chain conformational states. The results suggest that the experimental observation of ‘continuous melting’ in pure wet lipid bilayers, rather than being ascribed to the presence of impurities, may be explained as a result of kinetically caused metastability of intermediate lipid-chain conformations.     

Validation of all-atom phosphatidylcholine lipid force fields in the tensionless NPT ensemble

A recently defined charge set, to be used in conjunction with the all-atom CHARMM27r force field, has been validated for a series of phosphatidylcholine lipids. The work of Sonne et al. successfully replicated experimental bulk membrane behaviour for dipalmitoylphosphatidylcholine (DPPC) under the isothermal-isobaric (NPT) ensemble. Previous studies using the defined CHARMM27r charge set have resulted in lateral membrane contraction when used in the tensionless NPT ensemble, forcing the lipids to adopt a more ordered conformation than predicted experimentally. The current study has extended the newly defined charge set to 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC) and 1-palmitoyl-2-docosahexaenoyl-sn-glycero-3-phosphatidylcholine (PDPC). Molecular dynamics simulations were run for each of the lipids (including DPPC) using both the CHARMM27r charge set and the newly defined modified charge set. In all three cases a significant improvement was seen in both bulk membrane properties and individual atomistic effects. Membrane width, area per lipid and the depth of water penetration were all seen to converge to experimental values. Deuterium order parameters generated with the new charge set showed increased disorder across the width of the bilayer and reflected both results from experiment and similar simulations run with united atom models. These newly validated models can now find use in mixed biological simulations under the tensionless ensemble without concern for lateral contraction.     

Implicit membrane treatment of buried charged groups: Application to peptide translocation across lipid bilayers

The energetic cost of burying charged groups in the hydrophobic core of lipid bilayers has been controversial, with simulations giving higher estimates than certain experiments. Implicit membrane approaches are usually deemed too simplistic for this problem. Here we challenge this view. The free energy of transfer of amino acid side chains from water to the membrane center predicted by IMM1 is reasonably close to all-atom free energy calculations. The shape of the free energy profile, however, for the charged side chains needs to be modified to reflect the all-atom simulation findings (IMM1-LF). Membrane thinning is treated by combining simulations at different membrane widths with an estimate of membrane deformation free energy from elasticity theory. This approach is first tested on the voltage sensor and the isolated S4 helix of potassium channels. The voltage sensor is stably inserted in a transmembrane orientation for both the original and the modified model. The transmembrane orientation of the isolated S4 helix is unstable in the original model, but a stable local minimum in IMM1-LF, slightly higher in energy than the interfacial orientation. Peptide translocation is addressed by mapping the effective energy of the peptide as a function of vertical position and tilt angle, which allows identification of minimum energy pathways and transition states. The barriers computed for the S4 helix and other experimentally studied peptides are low enough for an observable rate. Thus, computational results and experimental studies on the membrane burial of peptide charged groups appear to be consistent. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.     

A theoretical model for the association of amphiphilic transmembrane peptides in lipid bilayers

A theoretical model is proposed for the association of trans-bilayer peptides in lipid bilayers. The model is based on a lattice model for the pure lipid bilayer, which accounts accurately for the most important conformational states of the lipids and their mutual interactions and statistics. Within the lattice formulation the bilayer is formed by two independent monolayers, each represented by a triangular lattice, on which sites the lipid chains are arrayed. The peptides are represented by regular objects, with no internal flexibility, and with a projected area on the bilayer plane corresponding to a hexagon with seven lattice sites. In addition, it is assumed that each peptide surface at the interface with the lipid chains is partially hydrophilic, and therefore interacts with the surrounding lipid matrix via selective anisotropic forces. The peptides would therefore assemble in order to shield their hydrophilic residues from the hydrophobic surroundings. The model describes the self-association of peptides in lipid bilayers via lateral and rotational diffusion, anisotropic lipid-peptide interactions, and peptide-peptide interactions involving the peptide hydrophilic regions. The intent of this model study is to analyse the conditions under which the association of trans-bilayer and partially hydrophilic peptides (or their dispersion in the lipid matrix) is lipid-mediated, and to what extent it is induced by direct interactions between the hydrophilic regions of the peptides. The model properties are calculated by a Monte Carlo computer simulation technique within the canonical ensemble. The results from the model study indicate that direct interactions between the hydrophilic regions of the peptides are necessary to induce peptide association in the lipid bilayer in the fluid phase. Furthermore, peptides within each aggregate are oriented in such a way as to shield their hydrophilic regions from the hydrophobic environment. The average number of peptides present in the aggregates formed depends on the degree of mismatch between the peptide hydrophobic length and the lipid bilayer hydrophobic thickness: The lower the degree of mismatch is the higher this number is. Received: 30 December 1996 / Accepted: 9 May 1997     

Relaxation process after the cooling jump across the pretransition of dipalmitoylphosphatidylcholine bilayers

The relaxation processes involved in the pretransition of dipalmitoylphosphatidylcholine (DPPC) multilamellar vesicles were investigated by the measurements of transmitted light intensity under crossed polars and freeze fracture electron microscopy. Temperature jump experiments from the ripple (P_β′) phase into the gel (L_β′) phase suggest that there must be at least three relaxation stages, one of which has a relaxation time much shorter than those described previously. Freeze fracture electron microscopy observations of the surface structure during the relaxation process reveals that the defects in the ripple pattern, i.e. the open terminal ends of each ripple ridge line, may be the primary origin of the growth of the L_β′ phase, and that the phase conversion may take place exclusively at the pointed end of a ripple ridge line. Thus the L_β′ phase grows almost one-dimensionally between ripple structures.     

Molecular dynamics simulation of the dissolution process of a cellulose triacetate-II nano-sized crystal in DMSO

An understanding of the dissolution process of cellulose derivatives is important not only for basic research but also for industrial purposes. We investigated the dissolution process of cellulose triacetate II (CTA II) nano-sized crystal in DMSO solvent using molecular dynamics simulations. The nano-sized crystal consists of 18 CTA chains. During the 9 ns simulation, it was observed that one chain (C01) located at the corner of the lozenge crystal was solvated by the DMSO molecules and moved away from the remaining cluster into the DMSO solvent. The analysis showed that the breakage of the interaction between the H1, H3, and H5 hydrogens of the pyranose ring and the acetyl carbonyl oxygen in the C01 and C02 adjacent chains would be crucial for the dissolution of CTA. The DMSO molecules solvating around these atoms would prevent the re-crystallization of the CTA molecules and facilitate further dissolution.