Dioleoylphosphatidylcholine Lipid Bilayers Research Articles

Structure and Dynamics of Dioleoyl-Phosphatidylcholine Bilayers under the Influence of Quercetin and Rutin.

Quercetin and rutin, two widely studied flavonoids with applications foreseen in the sectors of pharmaceutical and cosmetic industries, have been chosen as model compounds for a detailed structural and dynamical investigation onto their influence on fluid lipid bilayers. Combining global small angle X-ray scattering analysis with molecular dynamics, various changes in the properties of dioleoyl-phosphatidylcholine (DOPC) bilayers have been determined. The solubility of quercetin in DOPC membranes is assured up to 12 mol %, whereas rutin, with additional glucose and rhamnose groups, are fully soluble only up to 6 mol %. Both flavonoids induce an increase in membrane undulations and thin the bilayers slightly (<1 Å) in a concentration dependent manner, wherein quercetin shows a stronger effect. Concomitantly, in the order of 2-4%, the adjacent bilayer distance increases with the flavonoid's concentration. Partial molecular areas of quercetin and rutin are determined to be 26 and 51 Å2, respectively. Simulated averaged areas per molecule confirm these estimates. A 60° tilted orientation of quercetin is observed with respect to the bilayer normal, whereas the flavonoid moiety of rutin is oriented more perpendicular (α-angle 30°) to the membrane surface. Both flavonoid moieties are located at a depth of 12 and 16 Å for quercetin and rutin, respectively, while their anionic forms display a location closer to the polar interface. Finally, at both simulated concentrations (1.5 and 12 mol %), DOPC-rutin systems induce a stronger packing of the pure DOPC lipid bilayer, mainly due to stronger attractive electrostatic interactions in the polar lipid head region.

Kinetic evolution of DOPC lipid bilayers exposed to α-cyclodextrins.

Cyclodextrins are cyclic oligosaccharides capable of forming inclusion complexes with a variety of molecules, and as such have been recognized as a pharmaceutical and biotechnological asset. Cyclodextrins are known to interact with the components of cell membranes, and this correlates with a significant degree of cytotoxicity. In this work, we report on the mechanism of degradation of a model dioleoyl-phosphatidylcholine (DOPC) bilayer exposed to a solution with increasing concentrations of α-cyclodextrins. By combining optical fluorescence microscopy and quartz-crystal microbalance experiments, we study the evolution of supported lipid bilayers (SLBs) and giant unilamellar vesicles (GUVs). The rate of lipid removal is found to display a strong nonlinear dependence on the cyclodextrin concentration. A mechanism involving lipid aggregates is proposed.

Effects of cholesterol concentration on the interaction of cytarabine with lipid membranes: a molecular dynamics simulation study

Liposomal cytarabine, DepoCyt, is a chemotherapy agent which is used in cancer treatment. This form of cytarabine has more efficacy and fewer side effects relative to the other forms. Since DepoCyt contains the cytarabine encapsulated within phosphatidylcholine and the sterol molecules, we modeled dioleoylphosphatidylcholine (DOPC)/cholesterol bilayer membrane as a carrier for cytarabine to study drug–bilayer interactions. For this purpose, we performed a series of united-atom molecular dynamics (MD) simulations for 25 ns to investigate the interactions between cytarabine and cholesterol-containing DOPC lipid bilayers. Only the uncharged form of cytarabine molecule was investigated. In this study, different levels of the cholesterol content (0, 20, and 40%) were used. MD simulations allowed us to determine dynamical and structural properties of the bilayer membrane and to estimate the preferred location and orientation of the cytarabine molecule inside the bilayer membrane. Properties such as membrane thickness, area per lipid, diffusion coefficient, mass density, bilayer packing, order parameters, and intermolecular interactions were examined. The results show that by increasing the cholesterol concentration in the lipid bilayers, the bilayer thickness increases and area per lipid decreases. Moreover, in accordance with the experiments, our calculations show that cholesterol molecules have ordering effect on the hydrocarbon acyl chains. Furthermore, the cytarabine molecule preferentially occupies the polar region of the lipid head groups to form specific interactions (hydrogen bonds). Our results fully support the experimental data. Our finding about drug–bilayer interaction is crucial for the liposomal drug design.

Molecular dynamics studies of the antimicrobial peptides piscidin 1 and its mutants with a DOPC lipid bilayer

Piscidin 1 (Pis-1) has a high broad-spectrum activity against bacteria, fungi, and viruses but it also has a moderate hemolytic activities. To improve the antibacterial activity and to reduce toxicity, mutants Pis-1AA (G8A/G13A double mutant) and Pis-1PG (G8P mutant) have been designed based on the crystal structure of Pis-1. Eighteen independent molecular dynamics (MD) simulations of Pis-1 and its mutants with membranes are conducted in this article. Furthermore, 60 independent MD simulations of three peptides in water box have also been discussed for comparison. The results indicate that the unfolding process starts at the middle of the peptide. Pis-1 disrupts easily in the region of Val10-Lys14. Pis-1PG has a flexible N-terminal region, and the interaction between N-terminal and C-terminal is very weak. Pis-1AA has the most stable helical structure. In addition, percentage of native contacts and hydrogen bonds analysis are also performed. Lipid-peptide interaction analysis suggests that Pis-1 and Pis-1AA has a stronger interaction with the zwitterionic dioleoylphosphatidylcholine (DOPC) lipid bilayer than Pis-1PG. When compared with the results of peptide with membrane, peptides are unstable and unfolding quickly in water solution. Our results are applicable in examining diversities on hemolytic, antibacterial, and selectivity of antimicrobial peptides.

Role of Membrane Tension on Physical Properties of DOPC Lipid Bilayer

Mechanical stimuli trigger a wide range of cellular and biochemical processes in the cell. However the underlying molecular mechanisms of this triggering process are not well established. Numerous studies have demonstrated that in many cases mechanochemical signal conversion originates at the cell membrane. In the present study, a series of systematic molecular dynamics simulations were performed to characterize the changes in dioleoylphosphatidylcholine (DOPC) lipid bilayer properties in response to the membrane tension. All simulations were carried out using NAMD program and CHARMM36 force field. Our studies show that membrane tension causes changes in membrane properties such as area per lipid, volume, membrane thickness, lateral diffusion and dipole potential. The results also indicate that the lateral diffusion of lipid molecules is anomalous in nature due to the non-exponential distribution of waiting times. The simulation results were compared with experimental data obtained from fluorescence spectroscopy experiments using fluorescent lipid tracers and dipole-potential sensitive fluorescent probes.

A molecular view of the cholesterol condensing effect in DOPC lipid bilayers.

The condensing effect of cholesterol in dioleoylphosphatidylcholine (DOPC) lipid bilayers was systematically investigated via atomistic molecular dynamics (MD) simulation. Fourteen independent 200 ns simulations, spanning the entire range of cholesterol mole fraction (x(c)) in DOPC bilayers (i.e., from x(c) = 0 to 0.66), were performed at 323 K. The molecular areas occupied by DOPC and cholesterol at different distances from the bilayer center were analyzed using a slicing method based on the VDW radii of atoms. Curiously, while the average area per lipid and the cholesterol tilt angle, with respect to the bilayer normal, both show monotonic decreases as x(c) increases, the average bilayer height shows a significant decrease for x(c) > 0.35, following an initial increase. The calculated partial-specific areas of lipids clearly show the condensing effect of cholesterol. The VDW areal analysis showed that the condensing effect is limited only to the cholesterol sterol ring region, where the acyl chains of DOPC are severely compressed by cholesterol. As x(c) increases, the headgroups of DOPC gradually expand along the bilayer-aqueous interface to occupy more lateral area. Thus, it confirmed a key prediction of the umbrella model. At high cholesterol mole fractions, the calculated area per DOPC and area per cholesterol using some existing methods showed an inconsistent result: both increase while the overall area per lipid decreases. The inconsistency stems from the problematic assumption that cholesterol and DOPC have cylindrical shape and the same height. Our results showed that the total area of a PC/cholesterol bilayer is primarily determined by the molecular packing in the cholesterol sterol ring region. An alternative analysis of area per molecule within this region is proposed, which takes into account the cholesterol tilt angle and the practical incompressibility of cholesterol sterol rings. The new calculation shows that the majority of the area lost due to the cholesterol condensing effect is taken from PC molecules.

The MARTINI Coarse-Grained Force Field: Extension to Proteins.

Many biologically interesting phenomena occur on a time scale that is too long to be studied by atomistic simulations. These phenomena include the dynamics of large proteins and self-assembly of biological materials. Coarse-grained (CG) molecular modeling allows computer simulations to be run on length and time scales that are 2-3 orders of magnitude larger compared to atomistic simulations, providing a bridge between the atomistic and the mesoscopic scale. We developed a new CG model for proteins as an extension of the MARTINI force field. Here, we validate the model for its use in peptide-bilayer systems. In order to validate the model, we calculated the potential of mean force for each amino acid as a function of its distance from the center of a dioleoylphosphatidylcholine (DOPC) lipid bilayer. We then compared amino acid association constants, the partitioning of a series of model pentapeptides, the partitioning and orientation of WALP23 in DOPC lipid bilayers and a series of KALP peptides in dimyristoylphosphatidylcholine and dipalmitoylphosphatidylcholine (DPPC) bilayers. A comparison with results obtained from atomistic models shows good agreement in all of the tests performed. We also performed a systematic investigation of the partitioning of five series of polyalanine-leucine peptides (with different lengths and compositions) in DPPC bilayers. As expected, the fraction of peptides partitioned at the interface increased with decreasing peptide length and decreasing leucine content, demonstrating that the CG model is capable of discriminating partitioning behavior arising from subtle differences in the amino acid composition. Finally, we simulated the concentration-dependent formation of transmembrane pores by magainin, an antimicrobial peptide. In line with atomistic simulation studies, disordered toroidal pores are formed. In conclusion, the model is computationally efficient and effectively reproduces peptide-lipid interactions and the partitioning of amino acids and peptides in lipid bilayers.

Lateral Distribution of Cholesterol in Dioleoylphosphatidylcholine Lipid Bilayers: Cholesterol-Phospholipid Interactions at High Cholesterol Limit

Lateral organization of cholesterol in dioleoyl-phosphatidylcholine (DOPC) lipid bilayers at high cholesterol concentration (>45 mol%) was investigated using steady-state fluorescence anisotropy and fluorescent resonance energy transfer techniques. The recently devised Low Temperature Trap method was used to prepare compositionally uniform cholesterol/DOPC liposomes to avoid the problem of lipid demixing. The fluorescence anisotropy of diphenylhexatrience chain-labeled phosphatidylcholine (DPH-PC) in these liposomes exhibited local maxima at cholesterol mol fractions of 0.50 and 0.57, and a sharp drop at 0.67. For the liposomes labeled with both dehydroergosterol and DPH-PC, the fluorescent resonance energy transfer efficiency from dehydroergosterol to DPH-PC displayed a steep jump at cholesterol mol fraction of 0.5, and dips at 0.57 and 0.68. These results indicate the presence of highly ordered cholesterol regular distribution domains at those observed critical compositions. The observed critical mol fraction at 0.67 agreed favorably with the solubility limit of cholesterol in DOPC bilayers as independently measured by light scattering and optical microscopy. The regular distribution at 0.57 was previously predicted from a Monte Carlo simulation based on the Umbrella model. The results strongly support the hypothesis that the primary requirement for cholesterol-phospholipid mixing is that the polar phospholipid headgroups need to cover the nonpolar body of cholesterol to avoid the exposure of cholesterol to water.

Depolymerization of phospholamban in the presence of calcium pump: a fluorescence energy transfer study.

Phospholamban (PLB), a 52-amino acid protein, regulates the Ca-ATPase (calcium pump) in cardiac sarcoplasmic reticulum (SR) through PLB phosphorylation mediated by beta-adrenergic stimulation. The mobility of PLB on SDS-PAGE indicates a homopentamer, and it has been proposed that the pentameric structure of PLB is important for its regulatory function. However, the oligomeric structure of PLB must be determined in its native milieu, a lipid bilayer containing the Ca-ATPase. Here we have used fluorescence energy transfer (FET) to study the oligomeric structure of PLB in SDS and dioleoylphosphatidylcholine (DOPC) lipid bilayers reconstituted in the absence and presence of Ca-ATPase. PLB was labeled, specifically at Lys 3 in the cytoplasmic domain, with amine-reactive fluorescent donor/acceptor pairs. FET between donor- and acceptor-labeled subunits of PLB in SDS solution and DOPC lipid bilayers indicated the presence of PLB oligomers. The dependence of FET efficiency on the fraction of acceptor-labeled PLB in DOPC bilayers indicated that it is predominantly an oligomer having 9-11 subunits, with approximately 10% of the PLB as monomer, and the distance between dyes on adjacent PLB subunits is 0.9 +/- 0.1 nm. When labeled PLB was reconstituted with purified Ca-ATPase, FET indicated the depolymerization of PLB into smaller oligomers having an average of 5 subunits, with a concomitant increase in the fraction of monomer to 30-40% and a doubling of the intersubunit distance. We conclude that PLB exists primarily as an oligomer in membranes, and the Ca-ATPase affects the structure of this oligomer, but the Ca-ATPase binds preferentially to the monomer and/or small oligomers. These results suggest that the active inhibitory species of PLB is a monomer or an oligomer having fewer than 5 subunits.