Structural Characterization of Photopolymerizable Binary Liposomes Containing Diacetylenic and Saturated Phospholipids

The use of liposomes to encapsulate materials has received widespread attention for drug delivery, transfection, diagnostic reagent, and as immunoadjuvants. Phospholipid polymers form a new class of biomaterials with many potential applications in medicine and research. Of interest are polymeric phospholipids containing a diacetylene moiety along their acyl chain since these kinds of lipids can be polymerized by Ultra-Violet (UV) irradiation to form chains of covalently linked lipids in the bilayer. In particular the diacetylenic phosphatidylcholine 1,2-bis(10,12-tricosadiynoyl)- sn-glycero-3-phosphocholine (DC8,9PC) can form intermolecular cross-linking through the diacetylenic group to produce a conjugated polymer within the hydrocarbon region of the bilayer. As knowledge of liposome structures is certainly fundamental for system design improvement for new and better applications, this work focuses on the structural properties of polymerized DC8,9PC:1,2-dimyristoyl-sn-glycero-3-phusphocholine (DMPC) liposomes. Liposomes containing mixtures of DC8,9PC and DMPC, at different molar ratios, and exposed to different polymerization cycles, were studied through the analysis of the electron spin resonance (ESR) spectra of a spin label incorporated into the bilayer, and the calorimetric data obtained from differential scanning calorimetry (DSC) studies. Upon irradiation, if all lipids had been polymerized, no gel-fluid transition would be expected. However, even samples that went through 20 cycles of UV irradiation presented a DSC band, showing that around 80% of the DC8,9PC molecules were not polymerized. Both DSC and ESR indicated that the two different lipids scarcely mix at low temperatures, however few molecules of DMPC are present in DC8,9PC rich domains and vice versa. UV irradiation was found to affect the gel fluid transition of both DMPC and DC8,9PC rich regions, indicating the presence of polymeric units of DC8,9PC in both areas, A model explaining lipids rearrangement is proposed for this partially polymerized system.

Many-body energy decomposition of hydrogen-bonded glycine clusters in gas-phase

A detailed analysis of the many-body contribution to the interaction energies of the gas-phase hydrogen-bonded glycine clusters, (Gly)(N), N = 1-4 is presented. The energetics of the hydrogen-bonded dimer, trimer and tetramer complexes have been analyzed using density-functional theory. The magnitude of the two-through four-body energy terms have been calculated and compared. The relaxation energy and the two-body energy terms are the principal contributors to the total binding energy. Four-body contribution is negligible. However, the three-body contribution is found to be sizable and the formation of the cyclic glycine trimer presents geometric strains that make it less favorable. (C) 2010 Elsevier B.V. All rights reserved.

Thermal Phase Behavior of DMPG Bilayers in Aqueous Dispersions as Revealed by (2)H- and (31)P-NMR

The synthetic lipid 1,2-dimyristoyl-sn-3-phosphoglycerol (DMPG), when dispersed in water/NaCl exhibits a complex phase behavior caused by its almost unlimited swelling in excess water. Using deuterium ((2)H)- and phosphorus ((31)P)-NMR we have studied the molecular properties of DMPG/water/NaCl dispersions as a function of lipid and NaCl concentration. We have measured the order profile of the hydrophobic part of the lipid bilayer with deuterated DMPG while the orientation of the phosphoglycerol headgroup was deduced from the (31)P NMR chemical shielding anisotropy. At temperatures > 30 degrees C we observe well-resolved (2)H- and (31)P NMR spectra not much different from other liquid crystalline bilayers. From the order profiles it is possible to deduce the average length of the flexible fatty acyl chain. Unusual spectra are obtained in the temperature interval of 20-25 degrees C, indicating one or several phase transitions. The most dramatic changes are seen at low lipid concentration and low ionic strength. Under these conditions and at 25 degrees C, the phosphoglycerol headgroup rotates into the hydrocarbon layer and the hydrocarbon chains show larger flexing motions than at higher temperatures. The orientation of the phosphoglycerol headgroup depends on the bilayer surface charge and correlates with the degree of dissociation of DMPG-Na(+). The larger the negative surface charge, the more the headgroup rotates toward the nonpolar region.

Solvent effects on the electronic absorption spectrum of camphor using continuum, discrete or explicit approaches

We address the effect of solvation on the lowest electronic excitation energy of camphor. The solvents considered represent a large variation in-solvent polarity. We consider three conceptually different ways of accounting for the solvent using either an implicit, a discrete or an explicit solvation model. The solvatochromic shifts in polar solvents are found to be in good agreement with the experimental data for all three solvent models. However, both the implicit and discrete solvation models are less successful in predicting solvatochromic shifts for solvents of low polarity. The results presented suggest the importance of using explicit solvent molecules in the case of nonpolar solvents. (C) 2009 Elsevier B.V. All rights reserved.

NMR Chemical Shielding and Spin-Spin Coupling Constants of Liquid NH(3): A Systematic Investigation using the Sequential QM/MM Method

The NMR spin coupling parameters, (1)J(N,H) and (2)J(H,H), and the chemical shielding, sigma((15)N), of liquid ammonia are studied from a combined and sequential QM/MM methodology. Monte Carlo simulations are performed to generate statistically uncorrelated configurations that are submitted to density functional theory calculations. Two different Lennard-Jones potentials are used in the liquid simulations. Electronic polarization is included in these two potentials via an iterative procedure with and without geometry relaxation, and the influence on the calculated properties are analyzed. B3LYP/aug-cc-pVTZ-J calculations were used to compute the V(N,H) constants in the interval of -67.8 to -63.9 Hz, depending on the theoretical model used. These can be compared with the experimental results of -61.6 Hz. For the (2)J(H,H) coupling the theoretical results vary between -10.6 to -13.01 Hz. The indirect experimental result derived from partially deuterated liquid is -11.1 Hz. Inclusion of explicit hydrogen bonded molecules gives a small but important contribution. The vapor-to-liquid shifts are also considered. This shift is calculated to be negligible for (1)J(N,H) in agreement with experiment. This is rationalized as a cancellation of the geometry relaxation and pure solvent effects. For the chemical shielding, U(15 N) Calculations at the B3LYP/aug-pcS-3 show that the vapor-to-liquid chemical shift requires the explicit use of solvent molecules. Considering only one ammonia molecule in an electrostatic embedding gives a wrong sign for the chemical shift that is corrected only with the use of explicit additional molecules. The best result calculated for the vapor to liquid chemical shift Delta sigma((15)N) is -25.2 ppm, in good agreement with the experimental value of -22.6 ppm.

Solvent Effects in Chemical Processes. Water-Assisted Proton Transfer Reaction of Pterin in Aqueous Environment

Pterins are members of a family of heterocyclic compounds present in a wide variety of biological systems and may exist in two forms, corresponding to an acid and a basic tautomer. In this work, the proton transfer reaction between these tautomeric forms was investigated in the gas phase and in aqueous solution. In gas phase, the intramolecular mechanism was carried out for die isolated pterin by quantum mechanical second-order Moller-Plesset Perturbation theory (MP2/aug-cc-pVDZ) calculations and it indicates that the acid form is more stable than the basic form by -1.4 kcal/mol with a barrier of 34.2 kcal/mol with respect to the basic form. In aqueous solution, the role of the water molecules in the proton transfer reaction was analyzed in two separated parts, the direct participation of one water molecule in the reaction path, called water-assisted mechanism, and the complementary participation of the aqueous solvation. The water-assisted mechanism was carried out for one pterin-water cluster by quantum mechanical calculations and it indicates that the acid form is still more stable by -3.3 kcal/mol with a drastic reduction of 70% of the barrier, The bulk solution effect on the intramolecular and water-assisted mechanisms was included by free energy perturbation implemented on Monte Carlo simulations. The bulk water effect is found to be substantial and decisive when the reaction path involves the water-assisted mechanism. In this case, the free energy barrier is only 6.7 kcal/mol and the calculated relative Gibbs free energy for the two tautomers is -11.2 kcal/mol. This value is used to calculate the pK(a) value of 8.2 +/- 0.6 that is in excellent agreement with the experimental result of 7.9.

Dipole polarizability and Rayleigh light scattering by the hydrated electron

Assuming the existence of a confined state of the electron in bulk water the polarizability of the hydrated electron is analyzed. Statistically uncorrelated supermolecular structures composed of seven water molecules (first solvation shell) with an extra electron were extracted from classical Monte Carlo simulation and used in quantum mechanical second-order Moller-Plesset calculations. It is found that the bound excess electron contributes with 274 a.u. to the total dipole polarizability of 345 a.u. for (H(2)O)(7)(-). From the calculated polarizabilities the Rayleigh elastic light scattering properties are inferred and found to considerably enhance activity and light depolarization. (C) 2009 Elsevier B.V. All rights reserved.

The degree of polydispersity in micellar dispersions: A sum rule approach

Statistical properties of a two-dimensional ideal dispersion of polydisperse micelles are derived by analyzing the convergence properties of a sum rule set by mass conservation. Internal micellar degrees of freedom are accounted for by a microscopic model describing small displacements of the constituting amphiphiles with respect to their equilibrium positions. The transfer matrix (TM) method is employed to compute internal micelle partition function. We show that the conditions under which the sum rule is saturated by the largest eigenvalue of the TM determine the value of amphiphile concentration above which the dispersion becomes highly polydisperse and micelle sizes approach a Schultz distribution. (C) 2011 Elsevier B.V. All rights reserved.

Polarization and Spectral Shift of Benzophenone in Supercritical Water

Monte Carlo simulation and quantum mechanics calculations based on the INDO/CIS and TD-DFT methods were utilized to study the solvatochromic shift of benzophenone when changing the environment from normal water to supercritical (P = 340.2 atm and T = 673 K) condition. Solute polarization increases the dipole moment of benzophenone, compared to gas phase, by 88 and 35% in normal and supercritical conditions, giving the in-solvent dipole value of 5.8 and 4.2 D, respectively. The average number of solute-solvent hydrogen bonds was analyzed, and a large decrease of 2.3 in normal water to only 0.8 in the supercritical environment was found. By using these polarized models of benzophenone in the two different conditions of water, we performed MC simulations to generate statistically uncorrelated configurations of the solute surrounded by the solvent molecules and subsequent quantum mechanics calculations on these configurations. When changing from normal to supercritical water environment, INDO/CIS calculations explicitly considering all valence electrons of the 235 solvent water molecules resulted in a solvatochromic shift of 1425 cm(-1) for the most intense transition of benzophenone, that is, slightly underestimated in comparison with the experimentally inferred result of 1700 cm(-1). TD-B3LYP/6-311+G(2d,p) calculations on the same configurations but with benzophenone electrostatically embedded in the 320 water molecules resulted in a solvatochromic shift of 1715 cm(-1) for this transition, in very good agreement with the experimental result. When using the unpolarized model of the benzophenone, this calculated solvatochromic shift was only 640 cm(-1). Additional calculations were also made by using BHandHLYP/6-311+G(2d,p) to analyze the effect of the asymptotic decay of the exchange functional. This study indicates that, contrary to the general expectation, there is a sizable solute polarization even in the low-density regime of supercritical condition and that the inclusion of this polarization is important for a reliable description of the spectral shifts considered here.

Polarization and solvatochromic shift of ortho-betaine in water

Betaine dyes are known to show very large transition energy shifts in different solvents. The ortho-betaine molecule - a simple two-ring prototype of the E-T(30) Reichardt dye - has been investigated theoretically from a combined statistical and quantum mechanics approach. Using sequential Monte Carlo (MC) simulations and MP2/cc-pVDZ calculations the in-water dipole moment of ortho-betaine is obtained as 12.30 +/- 0.05 D. This result shows a considerable increase of 75% compared to the in-vacuum dipole moment. For comparison, the use of a polarizable continuum model using the same MP2/cc-pVDZ leads to an in-water dipole moment of 11.6 D, in good agreement. This large polarization is incorporated in the classical potential for another MC simulation to generate solute-solvent configurations and to obtain the contribution of the polarization effect in the solvatochromic shift. Using statistically uncorrelated configurations and supermolecular INDO/CIS calculations, including the solute and, explicitly, 230 solvent water molecules, the statistically converged calculated shift is obtained here as 6360 cm(-1), in good agreement with the experimental result of 7550 cm(-1). (c) 2007 Elsevier B.V. All rights reserved.

On the interaction of the anthraquinone barbaloin with negatively charged DMPG bilayers

Barbaloin is a bioactive glycosilated 1,8-dihydroxyanthraquinone present in several exudates from plants, Such as Aloe vera, which are used for cosmetic or food purposes. It has been shown that barbaloin interacts with DMPG (dimyristoylphosphatidylglycerol) model membranes, altering the bilayer structure (Alves, D. S.; Perez-Fons, L.; Estepa, A.; Micol, V. Biochem. Pharm. 2004, 68, 549). Considering that ESR (electron spin resonance) of spin labels is one of the best techniques to monitor structural properties at the molecular level, the alterations caused by the anthraquinone barbaloin on phospholipid bilayers will be discussed here via the ESR signal of phospholipid spin probes intercalated into the membranes. In DMPG at high ionic strength (10 mM Hepes pH 7.4 + 100 mM NaCl), a system that presents a gel-fluid transition around 23 degrees C, 20 mol % barbaloin turns the gel phase more rigid, does not alter much the fluid phase packing, but makes the lipid thermal transition less sharp. However, in a low-salt DMPG dispersion (10 mM Hepes pH 7.4 + 2 mM NaCl), which presents a rather complex gel-fluid thermal transition (Lamy-Freund, M. T.; Riske, K. A. Chem. Phys. Lipids 2003, 122, 19), barbaloin strongly affects bilayer structural properties, both in the gel and fluid phases, extending the transition region to much higher temperature values. The position of barbaloin in DMPG bilayers will be discussed on the basis of ESR results, in parallel with data from sample viscosity, DSC (differential scanning calorimetry), and SAXS (small-angle X-ray scattering).

Self-assembling of phenothiazine compounds investigated by small-angle X-ray scattering and electron paramagnetic resonance spectroscopy

Small-angle X-ray scattering (SAXS) and electron paramagnetic resonance (EPR) have been carried out to investigate the structure of the self-aggregates of two phenothiazine drugs, chlorpromazine (CPZ) and trifluoperazine (TFP), in aqueous solution. In the SAXS studies, drug solutions of 20 and 60 mM, at pH 4.0 and 7.0, were investigated and the best data fittings were achieved assuming several different particle form factors with a homogeneous electron density distribution in respect to the water environment. Because of the limitation of scattering intensity in the q range above 0.15 angstrom(-1), precise determination of the aggregate shape was not possible and all of the tested models for ellipsoids, cylinders, or parallelepipeds fitted the experimental data equally well. The SAXS data allows inferring, however, that CPZ molecules might self-assemble in a basis set of an orthorhombic cell, remaining as nanocrystallites in solution. Such nanocrystals are composed of a small number of unit cells (up to 10, in c-direction), with CPZ aggregation numbers of 60-80. EPR spectra of 5- and 16-doxyl stearic acids bound to the aggregates were analyzed through simulation, and the dynamic and magnetic parameters were obtained. The phenothiazine concentration in EPR experiments was in the range of 5-60 mM. Critical aggregation concentration of TFP is lower than that for CPZ, consistent with a higher hydrophobicity of TFP. At acidic pH 4.0 a significant residual motion of the nitroxide relative to the aggregate is observed, and the EPR spectra and corresponding parameters are similar to those reported for aqueous surfactant micelles. However, at pH 6.5 a significant motional restriction is observed, and the nitroxide rotational correlation times correlate very well with those estimated for the whole aggregated particle from SAXS data. This implies that the aggregate is densely packed at this pH and that the nitroxide is tightly bound to it producing a strongly immobilized EPR spectrum. Besides that, at pH 6.5 the differences in motional restriction observed between 5- and 16-DSA are small, which is different from that observed for aqueous surfactant micelles.

Calculations of vibrational frequencies, Raman activities and degrees of depolarization for complexes involving water, methanol and ethanol

Raman activities and degrees of depolarization are reported for 14 complexes involving methanol, ethanol and water using the MP2/aug-cc-pVDZ model. For ethanol both trans and gauche isomers are considered. The red-shifts of the OH stretching and the blue shifts of the bending tau(CO-OH) mode were analyzed for the proton-donor molecules upon hydrogen bond. The shift of the nu(CO) stretching mode of the alcohol molecules are also analyzed and found to be specific giving characterization of the amphoteric relation, being positive for the proton-acceptor and negative for the proton-donor molecule. (c) 2008 Elsevier B.V. All rights reserved.

Conformational behavior of different possible ways of oligoglycine formation in a solvent-free environment

The possible ways for glycine oligopeptide formation in gas phase, both in the extended P-strand like conformation and folded 2(7)-ribbon like conformations are analyzed using quantum chemical calculations. We focus on the sequential formation of peptide bond through upgradation of the immediate lower order molecule and observe the consequences in other related processes like oligoglycine formation through simultaneous peptide linkage of n glycine monomers and interchange of molecular conformation through peptide linkage. A comparison is made between the structures and binding energies obtained for both conformers. All binding energies are increased by the zero-point energy contribution. The role of electron correlation effects is briefly analyzed. The folded 2(7)-ribbon-like conformations in vacuo are found to be more stable in comparison to the extended structure. (c) 2007 Elsevier B.V. All rights reserved.

Electronic and magnetic properties of Ti and Fe on graphene

The structural, electronic and magnetic properties of Fe and Ti atomic wires and the complete covering when adsorbed on graphene are presented through ab initio calculations based on density functional theory. The most stable configurations are investigated for Fe and Ti in different concentrations adsorbed on the graphene surface, and the corresponding binding energies are calculated. The results show a tendency of the Ti atoms to cover uniformly the graphene surface, whereas the Fe atoms form clusters. The adsorption of the transition metal on the graphene surface changes significantly the electronic density of states near the graphene Fermi region. In all arrangements studied, a charge transfer is observed from the adsorbed species to the graphene surface due to the high hybridizations between the systems.