Shear Viscosity of Model Mixtures by Nonequilibrium Molecular Dynamics III. Effect of Quadrupolar Interactions

Abstract

We present new results for thermodynamic properties and viscosities of pure quadrupolar fluids, a pure dipolar quadrupolar fluid, nonquadrupolar/quadrupolar mixtures, and quadrupolar/quadrupolar mixtures. It is evident that, the addition of quadrupolar interactions to the pure Ar and the addition of quadrupolar interactions to the pure dipolar Ar, leads to higher viscosities as was observed in the addition of dipolar interactions to the pure Ar [Lee and Cummings, J. Chem. Phys., 105, 2044 (1996)]. The total energies and the mixture densities show a linear dependence for both nonquadrupolar Ar/quadrupolar Kr (case B) and quadrupolar Ar/quadrupolar Kr (case C), and the linearity of case C is better than that of case B. This is not consistent with the idea that in the cases of the dipolar mixtures, dipolar/dipolar and nondipolar/nondipolar mixture are likely to be more ideal than nondipolar/dipolar mixtures. This is mainly due to the weaker interaction of quadrupole-quadrupole than that of dipoledipole.     

Non Equilibrium Molecular Dynamics Simulation of Coupled Heat- and Mass Transport Across a Liquid/Vapor Interface

Abstract

The transport properties of bulk liquid, gas and at the gas/liquid interface were studied for two binary Lennard-Jones/spline mixtures by use of nonequilibrium molecular dynamics. One of the mixtures was an ideal isotope mixture, the other a non-ideal mixture. The simulations gave the thermal conductivity, mutual diffusion coefficient, heat flux, mass flux, and the changes in these quantities across the interface. The local entropy production was expressed in terms of fluxes and thermodynamic forces, and numercial estimates are given. It was shown that the largest contribution to the total entropy production occurs in the vapor phase under the chosen conditions. We expect, however that if the mass flux were larger, the major contribution to the entropy production would come from the liquid phase.     

The antitumor molecular mechanism of Alisma orientalis with c-myc DNA: multi-spectroscopic analysis and molecular simulation

Abstract

We prepared extracts of Alisma orientalis from Sichuan and Fujian Province, China. Based on the ratio of alisol B 23-acetate (23B) to alisol A 24-acetate (24A) in two Alisma orientalis extracts, we prepared two mixtures of 24A and 23B (24A:23B?=?1:3 or 1:10). The antitumor molecular mechanism of the monomers 24A and 23B, the two mixtures and the effective components of Alisma orientalis from different habitats were studied. The MTT assay suggested that the difference in the antitumor activity of Alisma orientalis from different habitats was correlated to the ratio of 24A to 23B. The multi-spectroscopic analysis suggested that the effective components, the monomers and mixtures interacted with c-myc DNA in a partial intercalation manner. The binding strength of the alisol acetates to c-myc DNA was consistent with the anticancer activity, indicating that c-myc DNA was the anticancer target. The molecular simulation indicated that the mixtures were all directly bound to different base pairs of c-myc DNA for a superimposed effect, which led to the binding strength of the mixtures to c-myc DNA was stronger than that of the monomers. The molecules in the 1:3 mixture were all bound to different base pairs of c-myc DNA. However, for the 1:10 mixture, seven molecules of 23B bound to the side chain of 24A, resulting in the mixture with a long chain structure which increased the steric hindrance of 24A. As a result, affinity between 24A and c-myc DNA in the 1:10 mixture was weaker than that in the 1:3 mixture.

The antitumor molecular mechanism of the alisol monomers 24A and 23B, the mixtures with different proportions and the effective components of Alisma orientalis from different habitats were studied. The order of the antitumor activity was as follows: Sichuan?>?Fujian, 24A-23B (1:3) > 24A-23B (1:10) > 23B?>?24A. The antitumor activity of Alisma orientalis from different habitats was consistent with the mixtures which were designed according to the contents of the active ingredients of the medicinal materials, indicating that the antitumor activity of Alisma orientalis from Sichuan is better than that from Fujian which is related to the contents of 24A and 23B and the proportion of 1:3 is better than 1:10. The binding strength of the mixtures to c-myc DNA was consistent with the anticancer activity. The mixtures were all directly bound to different base pairs of c-myc DNA for a superimposed effect, which led to the strength of the interaction of the mixtures to c-myc DNA was stronger than that of the monomers. For the 24A-23B (1:3) mixture, the four small molecules bound to c-myc DNA directly and interacted with different base pairs of c-myc DNA. While for the 24A-23B (1:10) mixture, 24A and three 23B molecules interacted with c-myc DNA, the remaining seven 23B molecules bound to the side chain of 24A, which increased the steric hindrance. The binding of the mixture to c-myc DNA was decreased.

Communicated by Ramaswamy H. Sarma     

A Study of Conformational Equilibria in Chain Molecules. I. Liquid n-Butane

Abstract

We have performed molecular dynamics simulations for liquid n-butane in order to understand liquid structures in terms of both inter- and intra-molecular interactions. Each n-butane molecule consists of four sites interacting with LJ potential and only a dihedral angle is taken into account as the internal degree of freedom. The population of gauche conformations with respect to the ideal gas state is found to increase in the liquid state. To investigate how the intermolecular interaction affects the dihedral angle distribution, we also adopt the repulsive LJ potential (RLJ) model. It is found that the nearest neighbor packing of the methyl and/or methylene groups can be approximately represented by using only the repulsive interaction. From the dihedral angle distribution, however, the rate of the shift of RLJ model to gauche is larger than that of LJ model and the attractive force also plays a significant role in the conformational equilibrium.     

Effect of temperature on adsorption of mixtures in porous materials

The effect of temperature on the adsorption of a simple mixture (Ar/Kr) in disordered porous materials is investigated by means of molecular simulation. In the larger mesopores of porous silica glasses, capillary condensation occurs upon decreasing the temperature. At temperatures above the capillary condensation temperature, Kr is preferentially adsorbed at the pore surface and Ar adsorption occurs in regions of low Kr density. For temperatures below the capillary condensation temperature, Ar density surprisingly increases as temperature increases, the behaviour that is consistent with an over-solubility effect. In contrast, in the disordered sub-nanoporous carbon, filling of the pores occurs in a reversible and continuous way upon decreasing the temperature, owing to the small size and amorphous shape of the pores. These results show that the crossover between capillary condensation and continuous reversible filling observed for pure fluids in pores also exists for mixtures. We also show that the Kr selectivity exhibits a minimum in the disordered porous silica that is located at the capillary condensation temperature. In contrast, in the disordered porous carbon where no capillary condensation occurs, the selectivity decreases monotonically with increasing the temperature. These results shed light on low-temperature adsorption of mixtures confined in porous materials and provide a guide to design efficient phase separation processes.     

The Rheology of Liquid Pentane Isomers by Non-Equilibrium Molecular Dynamics Simulations

Abstract

In this paper, non-equilibrium molecular dynamics (NEMD) simulations of planar Couette flow are reported for an expanded collapsed atom model for liquid pentane isomers at 273.15 K. The strain rate dependent viscosity for liquid pentane isomers exhibits shear-thinning and a linear dependence on γ^1/2. Newtonian viscosities for liquid pentane isomers obtained by a linear extrapolation to zero strain rate are: 0.256cP for normal pentane, 0.219cP for isopentane, and 0.168cP for neopentane. The strain rate dependent pressure difference and normal stress difference vary nearly linearly with the γ^3/2 law and the γ law, respectively, for all three liquid pentane isomers. The overall trend of the square of radius of gyration and end-to-end distance for normal pentane is a linear increase with strain rate. For isopentane, the trend hardly changes for the range of shear rate in this study. The alignment angle decreases with increasing strain rate and the alignment angle of the straight chain alkane is less than that of the branched chain alkane. The average percentage of C?C?C?C trans for normal pentane as a function of strain rate is in excellent correlation with the square of the radius of gyration and the average end-to-end distance. Applying the strain rate in the x-direction, the alignment angle is forced to decrease and the percentage of C?C?C?C trans increases with increasing strain rate.     

Simulation of Polymer Solutions by Dissipative Particle Dynamics

Abstract

Dissipative Particle Dynamics (DPD) is employed to model the dynamics and rheology of polymer solutions, and suspensions of spherical particles with adsorbed polymers. Static and dynamic scaling relationships for the variation of radius of gyration and relaxation time with polymer chain length are reviewed, demonstrating that the DPD polymer solution model correctly represents the effects of hydrodynamic interaction and excluded volume. Rheological simulations for both polymer solutions and polymer-sphere suspensions predict Newtonian viscosities at low shear rate followed by shear-thinning behavior as a reduced shear rate of unity is approached. Both the Newtonian viscosity and the extent of shear-thinning are greatly enhanced in the case of good solvents, compared to the viscosity curves for polymers and polymer-spheres structures dissolved in theta solvents and poor solvents.     

On the nature of bonding in HCOOH...Ar and HCOOH...Kr complexes

The chemical interaction in HCOOH...Ng (Ng=Ar, Kr) complex was analyzed by topological analysis of the electron density based on Atoms-In-Molecules theory. For all computationally stable equilibrium structures of 1:1 HCOOH...Ng complexes, an intermolecular bond path with a bond critical point was found and perturbation of formic acid (FA) atomic basins and electron density was observed. The intermolecular interaction between the two complex subunits can be classified, according to its nature, as a closed-shell van der Waals type of interaction. However, one of the computed structures (complex II), pictures a noble gas atom attached linearly to the acidic O–H tail of FA. In this particular case, the electron density at the intermolecular bond critical point was found to resemble a hydrogen-bonded system and thus, may be termed a hydrogen-bond-like interaction. This change in the nature of the interaction is also shown by large perturbations of the FA properties found for this complex structure. The structural and vibrational perturbations are larger than for the other three structures and they increase for the Kr complexes compared to the Ar complex.Figure Electron density analysis of HCOOH...Ng (Ng=Ar,Kr) complex.     

The Rheology of 6-Propyl Duodecane and 5-Dibutyl Nonane by Non-Equilibrium Molecular Dynamics Simulations

Abstract

In recent papers, we reported non-equilibrium molecular dynamics (NEMD) simulations of planar Couette flow for liquid n- and i-butane, and liquid n-decane and 4-propyl heptane, using two collapsed atom models and an atomistically detailed model. It was found that the collapsed atom models predict the viscosities of the n-butane and n-decane quite well, and that the atomistically detailed model does not yield quantitative agreement with the viscosity of the n-alkanes or the branched alkanes, but it does have the one positive feature that the calculated viscosities of the branched alkanes are higher than these of the n-alkanes. In the present paper, we report results of NEMD simulations of planar Couette flow for liquid 6-propyl duodecane and 5-dibutyl nonane at 296 K and 0.782 g/cc, using an expanded collapsed atom model for simplicity. The strain rate dependent viscosity shows three different regions—two shear thinning ones and a Newtonian one. The slopes of the log-log plot for the branched-chain alkanes at high strain rate are different from those at intermediate strain rate, which is characterized as a rheological behavior of branched-chain alkanes. The Newtonian viscosity of the branched-chain alkanes can be extrapolated from the plateau value of the strain rate dependent viscosity at low strain rate as for straight-chain alkanes [J. Chem. Phys., 105, 1214 (1996)]. The results indicate that more-branched alkanes have a larger viscosity than less-branched C₁₇ alkanes.     

Studies of DNA dumbbells. V. A DNA triplex formed between a 28 base-pair DNA dumbbell substrate and a 16 base linear single strand

CD spectra and melting curves were collected for a 28 base-pair DNA fragment in the form of a DNA dumbbell (linked on both ends by T₄ single-strand loops) and the same DNA sequence in the linear form (without end loops). The central 16 base pairs (bp) of the 28-bp duplex region is the poly(pu) sequence: 5′-AGGAAGGAGGAAAGAG-3′. Mixtures of the dumbbell and linear DNAs with the 16-base single-strand sequence 5′-TCCTTCCTCCTTTCTC-3′ were also prepared and studied. At 22°C, CD measurements of the mixtures in 950 mM NaCl, 10 mM sodium acetate, 1 mM EDTA, pH 5.5, at a duplex concentration of 1.8 μM, provided evidence for triplex formation. Spectroscopic features of the triplexes formed with either a dumbbell or linear substrate were quite similar. Melting curves of the duplex molecules alone and in mixtures with the third strand were collected as a function of duplex concentration from 0.16 to 2.15 μM. Melting curves of the dumbbell alone and mixtures with the third strand were entirely independent of DNA concentration. In contrast, melting curves of the linear duplex alone or mixed with the third strand were concentration dependent. At identical duplex concentrations, the dumbbell alone melts ～20°C higher than the linear duplex. The curve of the linear duplex displayed a significant pretransition probably due to end fraying. On melting curves of mixtures of the dumbbell or linear duplex with the third strand, a low temperature transition with much lower relative hyperchromicity change (～ 5%) was observed. This transition was attributed to the melting of a new molecular species, e.g., the triplex formed between the duplex and single-strand DNA molecules. In the case of the dumbbell/single-strand mixture, these melting transitions of the triplex and the dumbbell were entirely resolvable. In contrast, the melting transitions of the linear duplex and the triplex overlapped, thereby preventing their clear distinction. To analyze the data, a three-state equilibrium model is presented. The analysis utilizes differences in relative absorbance vs temperature curves of dumbbells (or linear molecules) alone and in mixtures with the third strand. From the model analysis a straightforward derivation of f_T(T), the fraction of triplex as a function of temperature, was obtained. Analysis of f_T vs temperature curves, in effect melting curves of the triplexes, provided evaluation of thermodynamic parameters of the melting transition. For the triplex formed with the dumbbell substrate, the total transition enthalpy is ΔH_T = 118.4 ± 12.8 kcal/mol (7.4 ± 0.8 kcal/mol per triplet unit) and the total transition entropy is ΔS_T = 344 ± 36.8 cal/K · mol (eu) (21.5 ± 2.3 eu per triple unit). The transition curves of the triplex formed with the linear duplex substrate displayed two distinct regions. A broad pretransition region from f_T = 0 to 0.55 and a higher, sharper transition above f_T = 0.55. The transition parameters derived from the lower temperature region of the curve are ΔH′_T = 44.8 ± 9.6 kcal/mol and ΔS′_T = 112 ± 33.6 eu (or ΔH′ = 2.8 ± 0.6 kcal/mol and ΔS′ = 7.0 ± 2.1 eu per triplet). These values are probably too small to correspond to actual melting of the triplex but instead likely reveal effects of end fraying of the duplex substrate on triplex stability. Transition parameters of the upper transition are ΔH′_T = 128.0 ± 2.3 kcal/mol and ΔS′_T = 379.2 ± 6.4 eu (ΔH′ = 8.0 ± 0.2 kcal/mol and ΔS′ = 23.7 ± 0.4 eu per triplet) in good agreement (within experimental error) with the transition parameters of the triplex formed with the dumbbell substrate. Supposing this upper transition reflects actual dissociation of the third strand from the linear duplex substrate this triplex is comparable in thermodynamic stability to the triplex formed with a dumbbell substrate. Even so, the biphasic melting character of the linear triplex obscures the whole analysis, casting doubt on its absolute reliability. Apparently triplexes formed with a dumbbell substrate offer technical advantages over triplexes formed from linear or hairpin duplex substrates for studies of DNA triplex stability. © 1993 John Wiley & Sons, Inc.     

Chemical Potential,Partial Enthalpy and Partial Volume of Mixtures by NPT Molecular Dynamics

Abstract

Equilibrium NPT molecular dynamics computer simulations have been used to determine the chemical potential, partial enthalpy and partial volume of model Ar-Kr mixtures using newly devised non-intrusive particle insertion and particle swap techniques [P. Sindzingre et al. Chemical Physics, 129 (1989) 213]. In this report we examine, for the first time, in some detail the relative convergence statistics of the particle swap and particle insertion methods for these properties for binary Lennard-Jones (LJ) mixtures. Both species are represented by single-site Lennard-Jones pair potentials with Lorentz-Berthelot rules for the cross-species interactions. We show that, over the whole phase diagram and especially in the vicinity of the fluid-solid coexistence line, the particle swap method gives significantly better statistics than the particle insertion method for the difference in chemical potential of the two species, partial enthalpy and partial volume of each species. Also, we find that, using the particle swap method, the difference in the chemical potential converges more rapidly than the differences in the partial enthalpy and volume.     

Thermodynamics of phosphatidylcholine-cholesterol mixed model membranes in the liquid crystalline state studied by the orientational order parameter.

It is shown that good estimates of the activity of cholesterol in phosphatidylcholine-cholesterol mixed model membranes are obtained by examining the orientational order parameter S of cholestane spin probe (CSL) that is obtained from electron spin resonance by spectral simulation. By introducing thermodynamic stability conditions of liquid mixtures, the variation of activity (or S) as a function of cholesterol mole fraction is utilized to predict the concentration at which the phase separation occurs. These results for DMPC and cholesterol binary mixtures agree very well with those of Tempo-partitioning experiments. The comparison of activity coefficients and the phase boundary in DMPC/cholesterol mixtures with those of POPC/cholesterol mixtures suggests that acyl chain unsaturation leads to poorer mixing of cholesterol in phosphatidylcholine model membranes at higher temperatures (i.e., greater than 35 degrees C). In ternary solutions of DMPC, POPC, and cholesterol, it is found that cholesterol shows less deviation from ideality than in either of the two binary mixtures, and this implies that the phase separation occurs at higher cholesterol concentration than in either of the two binary mixtures. The present analysis suggests that there may not be a critical point in DMPC/cholesterol mixtures, even though phase separation does occur.     

Micellization of fatty acyl-CoA mixtures and its relevance to the fatty acyl selectivity of acyltransferases

The critical micelle concentrations (CMCs) of palmitoyl-CoA/stearoyl-CoA and palmitoyl-CoA/oleoyl-CoA mixtures in 0.050 M KPi, pH 7.4, a buffer used in enzymatic studies, were determined by fluorescence. Mixed micelle solution theory, analogous to the thermodynamic treatment of vapor pressure, was applied to calculate monomer and micelle compositions. The behavior of the palmitoyl-CoA/stearoyl-CoA mixture is ideal, while the palmitoyl-CoA/oleoyl-CoA mixture, although not exhibiting ideal behavior, can be fitted reasonably well by nonideal theory. In both mixtures, selective micellization takes place and, unlike the case of pure fatty acyl-CoAs, above the CMC of the mixtures the concentration of molecules free in solution is strongly dependent upon total concentration. The information derived from the present physical studies becomes important in enzymatic studies with membrane-bound acyltransferases, where selectivity toward various fatty acyl donors, presented as binary mixtures, is frequently observed.     

Enhancement of the viscosity of mucin by serum albumin.

The interaction of serum albumin with a model epithelial mucin from pig stomach was explored by rotary viscometry. During 30 min of incubation of human serum albumin(20mg/ml) and pig gastric mucin (8mg/ml) in iso-osmotic buffers at 37 degrees C, the solution became markedly viscous. Viscosity enhancement was proportional to albumin concentration (2-40mg/ml), was most pronounced under conditions of low shear rate (less than 45S-1), and was considerably greater than the additive or multiplicative viscosity values calculated from albumin or mucin solutions measured separately. The viscous mucin-albumin complex was destroyed by high shear rates (greater than 90S-1), but slowly re-formed under zero shear conditions. Elevation of pH (7 to 9), ionic strength (0.1 to 1.0), and addition of disodium EDTA (5mM) did not cause marked or specific alterations in the viscosity of the mixture, suggesting that electrostatic interactions probably do not stabilize mucin-albumin complexes. Urea (7M) and heating (35 to 55 degrees C) caused a major increase in the viscosity of mucin and mucin-albumin mixtures, suggesting that rupture of hydrogen bonds, unfolding and partial denaturation of mucin promotes greater intertangling (possibly hydrophobic interactions) between mucin and albumin molecules. The implications of mucin-albumin interaction in diseases associated with mucus obstruction are briefly discussed.     

Molecular dynamics simulations of noble gases encapsulated in C60 Fullerene

Molecular dynamics simulations of Helium (He), Neon (Ne), Argon (Ar), Krpton (Kr) and Xenon (Xe) encapsulated in C₆₀ are discussed, as well simulations of Fullerenes containing anywhere from two to four He atoms. Even for single atom encapsulates, no species resides at the geometric center of the Fullerene cage. Smaller atoms sit more off-center than larger ones, and He appears to be a special case in both centering and dynamics. Some encapsulated species stabilize the cage by stifling radial fluctuations and others disrupt it; adding Ne seems to have the most stabilizing effect, while Kr and Xe cause the largest radial atomic excursions. Multiple He encapsulates tend to stabilize the cage; such systems are very stressed and show structure over a wide temperature range. Based on dynamical information quadruple He seems to be close to the packing limit for C₆₀.     

Molecular dynamics simulation study of probe diffusion in liquid n-alkanes

We performed molecular dynamics simulations for the probe diffusion and friction dynamics of Lennard-Jones (LJ) particles modelled for methyl yellow (MY) in liquid n-alkanes of C₁₂–C₂₀₀ at temperatures of 318, 418, 518 and 618 K. Two LJ particles are chosen: MY1 with a mass of 114 g/mol, LJ parameters of σ = 4.0 Å and ? = 0.4 kJ/mol, and MY2 with a mass of 225 g/mol, σ = 6.0 Å and ? = 0.6 kJ/mol. We observed a clear transition in the power law dependence of MY2 diffusion on the molecular weight of n-alkanes at lower temperatures of 318 and 418 K. The sharp transitions occur near n-dotriacontane (C₃₂). However, no such transition is found for MY1 at all the temperatures and for MY2 at higher temperatures of 518 and 618 K. We also calculated the friction constants of both MY probe molecules in liquid n-alkanes. For the larger probe molecule (MY2), at lower temperatures, a large deviation of slope from the linear dependence of the friction of MY2 on the chain length of n-alkane is observed, which indicates a large reduction of friction in longer chains when compared with the shorter chains, enhancing the diffusion of the probe molecules (MY2).     

Variants of a cloned synthetic lactose operator: I. A palindromic dimer lactose operator derived from one strand of the cloned 40-base pair operator

Starting with one strand of the 40-bp synthetic operator (Sadler et al., 1978), we have constructed and cloned a 66-bp, palindromic DNA segment with the following sequence

where the horizontal arrows indicate the locations of the two 21-bp “core? operator sequences in this segment and the vertical arrow designates the dyad axis of symmetry. Upon denaturation and rapid renaturation, each strand of this fragment forms a hairpin molecule still retaining an EcoRI cohesive end. Two hairpin molecules can be joined with T4 DNA ligase to form a duplex DNA molecule having no ends (dumbbell form A). Denaturation and rapid renaturation of dumbbell A yields a mixture of two dumbbell forms: dumbbell A which is a substrate for EcoKL, and a new form, dumbbell B, which is not a substrate. Each of the conformations of this DNA fragment have been purified and all are active in binding lactose repressor in vitro.     

Effects of temperature and composition on the viscosity of respiratory gases

The steady-state sensitivity of resistance pneumotachographs is proportional to viscosity. Dynamic characteristics of pneumotachographs, pressure transducers, and mass spectrometers are also viscosity dependent. We derive linear equations to approximate the viscosities of O2, N2, CO2, H2O, He, N2O, and Ar for temperatures between 20 and 40 degrees C by using published viscosity data and a nonlinear extrapolation equation. We verify the accuracy of the extrapolation equation by comparison with published data. Our linear equations for pure gas viscosities yield standard errors less than 0.35 microP. We also compare a nonlinear equation for calculating the viscosities of mixtures of gases with published measured viscosities of dry air, humid air, and He-O2 and N2-CO2 mixtures. The maximum difference between published and calculated values is 1.3% for 10% CO2 in N2. All other differences are less than 0.38%. For saturated humid air at 35 degrees C, a linear concentration-weighted combination of viscosities differs from our nonlinear equation by 4.9, 2.1, and 1.7% at barometric pressures of 32, 83, and 100 kPa, respectively. By use of our method, the viscosity of normal respiratory gases can be calculated to within 1% of measured values.     

Properties and Anti-HIV Activity of Nicked Dumbbell Oligonucleotides

Abstract

We have designed a new type of oligodeoxyribonucleotide. These oligodeoxyribonucleotides form two hairpin loop structures with base pairs (sense and antisense) in the double helical stem at the 3′ and 5′-ends (nicked dumbbell oligonucleotides). The nicked dumbbell oligonucleotides are molecules with free ends that are more resistant to exonuclease attack. Furthermore, the nicked dumbbell oligonucleotide containing phosphorothioate (P=S) bonds in the hairpin loops has increased nuclease resistance, as compared to the unmodified nicked oligonucleotide. The binding of the nicked dumbbell oligonucleotide to RNA is lower than that of a single-stranded DNA. We also describe the anti-HIV activity of nicked dumbbell oligonucleotides.

     

Two photon fluorescence microscopy of coexisting lipid domains in giant unilamellar vesicles of binary phospholipid mixtures

Images of giant unilamellar vesicles (GUVs) formed by different phospholipid mixtures (1,2-dipalmitoyl-sn-glycero-3-phosphocholine/1, 2-dilauroyl-sn-glycero-3-phosphocholine (DPPC/DLPC) 1:1 (mol/mol), and 1,2-dipalmitoyl-sn-glycero-3-phosphoethanolamine/1, 2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPE/DPPC), 7:3 and 3:7 (mol/mol) at different temperatures were obtained by exploiting the sectioning capability of a two-photon excitation fluorescence microscope. 6-Dodecanoyl-2-dimethylamino-naphthalene (LAURDAN), 6-propionyl-2-dimethylamino-naphthalene (PRODAN), and Lissamine rhodamine B 1,2-dihexadecanoyl-sn-glycero-3-phosphoethanolamine (N-Rh-DPPE) were used as fluorescent probes to reveal domain coexistence in the GUVs. We report the first characterization of the morphology of lipid domains in unsupported lipid bilayers. From the LAURDAN intensity images the excitation generalized polarization function (GP) was calculated at different temperatures to characterize the phase state of the lipid domain. On the basis of the phase diagram of each lipid mixture, we found a homogeneous fluorescence distribution in the GUV images at temperatures corresponding to the fluid region in all lipid mixtures. At temperatures corresponding to the phase coexistence region we observed lipid domains of different sizes and shapes, depending on the lipid sample composition. In the case of GUVs formed by DPPE/DPPC mixture, the gel DPPE domains present different shapes, such as hexagonal, rhombic, six-cornered star, dumbbell, or dendritic. At the phase coexistence region, the gel DPPE domains are moving and growing as the temperature decreases. Separated domains remain in the GUVs at temperatures corresponding to the solid region, showing solid-solid immiscibility. A different morphology was found in GUVs composed of DLPC/DPPC 1:1 (mol/mol) mixtures. At temperatures corresponding to the phase coexistence, we observed the gel domains as line defects in the GUV surface. These lines move and become thicker as the temperature decreases. As judged by the LAURDAN GP histogram, we concluded that the lipid phase characteristics at the phase coexistence region are different between the DPPE/DPPC and DLPC/DPPC mixtures. In the DPPE/DPPC mixture the coexistence is between pure gel and pure liquid domains, while in the DLPC/DPPC 1:1 (mol/mol) mixture we observed a strong influence of one phase on the other. In all cases the domains span the inner and outer leaflets of the membrane, suggesting a strong coupling between the inner and outer monolayers of the lipid membrane. This observation is also novel for unsupported lipid bilayers.