Effect of headgroup size, charge, and solvent structure on polymer−micelle interactions, studied by molecular dynamics simulations

We performed atomistic molecular dynamics simulations of anionic and cationic micelles in the presence of poly(ethylene oxide) (PEO) to understand why nonionic water-soluble polymers such as PEO interact strongly with anionic micelles but only weakly with cationic micelles. Our micelles include sodium n-dodecyl sulfate (SDS), n-dodecyl trimethylammonium chloride (DTAC), n-dodecyl ammonium chloride (DAC), and micelles in which we artificially reverse the sign of partial charges in SDS and DTAC. We observe that the polymer interacts hydrophobically with anionic SDS but only weakly with cationic DTAC and DAC, in agreement with experiment. However, the polymer also interacts with the artificial anionic DTAC but fails to interact hydrophobically with the artificial cationic SDS, illustrating that large headgroup size does not explain the weak polymer interaction with cationic micelles. In addition, we observe through simulation that this preference for interaction with anionic micelles still exists in a dipolar "dumbbell" solvent, indicating that water structure and hydrogen bonding alone cannot explain this preferential interaction. Our simulations suggest that direct electrostatic interactions between the micelle and polymer explain the preference for interaction with anionic micelles, even though the polymer overall carries no net charge. This is possible given the asymmetric distribution of negative charges on smaller atoms and positive charges oil larger units in the polymer chain.

Molecular dynamics simulations of threadlike cetyltrimethylammonium chloride micelles: effects of sodium chloride and sodium salicylate salts

We use atomistic molecular dynamics simulations to probe the effects of added sodium chloride (NaCl) and sodium salicylate (NaSal) salts on the spherical-to-threadlike micelle shape transition in aqueous solutions of cetyltrimethylammonium chloride (CTAC) surfactants. Long threadlike micelles are found to be unstable and break into spherical micelles at low concentrations or NaCl, but remain stable for 20 ns above a threshold value of [NaCl] approximate to 3.0 M, which is about 2.5 times larger than the experimental salt concentration at which the transition between spherical and rodlike micelles occurs. The chloride counterions associate weakly oil the surface of the CTAC micelles with the degree of counterion dissociation decreasing slightly with increasing [NaCl] on spherical micelles, but dropping significantly on the threadlike micelles tit high [NaCl]. This effect indicates that the electrolyte ions drive the micellar shape transition by screening the electrostatic repulsions between the micellar headgroups, The aromatic salicylate counterions, on the other hand, penetrate inside the micelle with their hydrophilic groups staying in the surfactant headgroup region and the hydrophobic groups partially embedded into the hydrophobic core of the micelle. The strong association of the salicylate ions with the surfactant headgroups leads to dense packing of the surfactant molecules, which effectively reduces the surface area per surfactant, and increases intramicellar ordering of the surfactant headgroups, favoring the formation of long threadlike micelles. Simulation predictions of the geometric and electrostatic properties of the spherical and threadlike micelles are in good agreement with experiments.

Molecular dynamics simulation of interactions between a sodium dodecyl sulfate micelle and a poly(ethylene oxide) polymer

We have performed atomistic molecular dynamics simulations of an anionic sodium dodecyl sulfate (SDS) micelle and a nonionic poly(ethylene oxide) (PEO) polymer in aqueous solution. The micelle consisted of 60 surfactant molecules, and the polymer chain lengths varied from 20 to 40 monomers. The force field parameters for PEO were adjusted by using 1,2-dimethoxymethane (DME) as a model compound and matching its hydration enthalpy and conformational behavior to experiment. Excellent agreement with previous experimental and simulation work was obtained through these modifications. The simulated scaling behavior of the PEO radius of gyration was also in close agreement with experimental results. The SDS-PEO simulations show that the polymer resides on the micelle surface and at the hydrocarbon-water interface, leading to a selective reduction in the hydrophobic contribution to the solvent-accessible surface area of the micelle. The association is mainly driven by hydrophobic interactions between the polymer and surfactant tails, while the interaction between the polymer and sulfate headgroups on the micelle surface is weak. The 40-monomer chain is mostly wrapped around the micelle, and nearly 90% of the monomers are adsorbed at low PEO concentration. Simulations were also performed with multiple 20-monomer chains, and gradual addition of polymer indicates that about 120 monomers are required to saturate the micelle surface. The stoichiometry of the resulting complex is in close agreement with experimental results, and the commonly accepted "beaded necklace" structure of the SDS-PEO complex is recovered by our simulations.

Time-resolved gas-phase kinetic study of the germylene addition reaction, CH3C CCH3

Time-resolved kinetic studies of the reaction of germylene, GeH2, generated by laser. ash photolysis of 3,4-dimethyl-1-germacyclopent-3-ene, have been carried out to obtain rate constants for its bimolecular reaction with 2-butyne, CH3C CCH3. The reaction was studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas, at five temperatures in the range 300-556 K. The second order rate constants obtained by extrapolation to the high pressure limits at each temperature, fitted the Arrhenius equation: log(k(infinity)/cm(3) molecule(-1) s(-1)) = (-10.46 +/- 10.06) + (5.16 +/- 10.47) kJ mol(-1)/ RT ln 10 Calculations of the energy surface of the GeC4H8 reaction system were carried out employing the additivity principle, by combining previous quantum chemical calculations of related reaction systems. These support formation of 1,2-dimethylvinylgermylene (rather than 2,3-dimethylgermirene) as the end product. RRKM calculations of the pressure dependence of the reaction are in reasonable agreement with this finding. The reactions of GeH2 with C2H2 and with CH3CRCCH3 are compared and contrasted.

Gas-phase kinetics of chlorosilylene reactions. I. ClSiH + Me3SiH: absolute rate measurements and theoretical calculations for prototype Si−H insertion reactions

Time-resolved studies of chlorosilylene, CISiH, generated by the 193 nm laser flash photolysis of 1-chloro-1-silacyclopent-3-ene, have been carried out to obtain rate constants for its bimolecular reaction with trimethylsilane, Me3SiH, in the gas phase. The reaction was studied at total pressures up to 100 torr (with and without added SF6) over the temperature range 297-407 K. The rate constants were found to be pressure independent and gave the following Arrhenius equation: log(k/cm(3) molecule(-1) s(-1)) = (-13.97 +/- 0.25) + (12.57 +/- 1.64) kJ mol(-1)/RT In 10. The Arrhenius parameters are consistent with a mechanism involving an intermediate complex, whose rearrangement is the rate-determining step. Quantum chemical calculations of the potential energy surface for this reaction and also the reactions of CISiH with SiH4 and the other methylsilanes support this conclusion. Comparisons of both experiment and theory with the analogous Si-H insertion processes of SiH2 and SiMe2 show that the main factor causing the lower reactivity of ClSiH is the secondary energy barrier. The calculations also show the existence of a novel intramolecular H-atom exchange process in the complex of ClSiH with MeSiH3.

The gas-phase reaction between silylene and 2-butyne: kinetics, isotope studies, pressure dependence studies and quantum chemical calculations

Time-resolved kinetic studies of the reactions of silylene, SiH2, and dideutero-silylene, SiD2, generated by laser. ash photolysis of phenylsilane and phenylsilane-d(3), respectively, have been carried out to obtain rate coefficients for their bimolecular reactions with 2-butyne, CH3C CCH3. The reactions were studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas at five temperatures in the range 294-612 K. The second-order rate coefficients, obtained by extrapolation to the high pressure limits at each temperature, fitted the Arrhenius equations where the error limits are single standard deviations: log(k(H)(infinity)/cm(3) molecule(-1) s(-1)) = (-9.67 +/- 0.04) + (1.71 +/- 0.33) kJ mol(-1)/RTln10 log(k(D)(infinity)/cm(3) molecule(-1) s(-1)) = (-9.65 +/- 0.01) + (1.92 +/- 0.13) kJ mol(-1)/RTln10 Additionally, pressure-dependent rate coefficients for the reaction of SiH2 with 2-butyne in the presence of He (1-100 Torr) were obtained at 301, 429 and 613 K. Quantum chemical (ab initio) calculations of the SiC4H8 reaction system at the G3 level support the formation of 2,3-dimethylsilirene [cyclo-SiH2C(CH3)=C(CH3)-] as the sole end product. However, reversible formation of 2,3-dimethylvinylsilylene [CH3CH=C(CH3)SiH] is also an important process. The calculations also indicate the probable involvement of several other intermediates, and possible products. RRKM calculations are in reasonable agreement with the pressure dependences at an enthalpy value for 2,3-dimethylsilirene fairly close to that suggested by the ab initio calculations. The experimental isotope effects deviate significantly from those predicted by RRKM theory. The differences can be explained by an isotopic scrambling mechanism, involving H - D exchange between the hydrogens of the methyl groups and the D-atoms in the ring in 2,3-dimethylsilirene-1,1-d(2). A detailed mechanism involving several intermediate species, which is consistent with the G3 energy surface, is proposed to account for this.

Influence of end-capping on the self-assembly of model amyloid peptide fragments

The influence of charge and aromatic stacking interactions on the self-assembly of a series of four model amyloid peptides has been examined. The four model peptides are based on the KLVFF motif from the amyloid Beta peptide, ABeta(16-20) extended at the N terminus with two Beta-alanine residues. We have studied NH2-BetaABetaAKLVFF-COOH (FF), NH2-BetaABetaAKLVFCOOH (F), CH3CONH-BetaABetaAKLVFF-CONH2 (CapF), and CH3CONH-BetaABetaAKLVFFCONH2 (CapFF). The former two are uncapped (net charge plus 2) and differ by one hydrophobic phenylalanine residue; the latter two are the analogous capped peptides (net charge plus 1). The self-assembly characteristics of these peptides are remarkably different and strongly dependent on concentration. NMR shows a shift from carboxylate to carboxylic acid forms upon increasing concentration. Saturation transfer measurements of solvent molecules indicate selective involvement of phenylalanine residues in driving the self-assembly process of CapFF due presumably to the effect of aromatic stacking interactions. FTIR spectroscopy reveals beta-sheet features for the two peptides containing two phenylalanine residues but not the single phenylalanine residue, pointing again to the driving force for self-assembly. Circular dichroism (CD) in dilute solution reveals the polyproline II conformation, except for F which is disordered. We discuss the relationship of this observation to the significant pH shift observed for this peptide when compared the calculated value. Atomic force microscopy and cryogenic-TEM reveals the formation of twisted fibrils for CapFF, as previously also observed for FF. The influence of salt on the self-assembly of the model beta-sheet forming capped peptide CapFF was investigated by FTIR. Cryo-TEM reveals that the extent of twisting decreases with increased salt concentration, leading to the formation of flat ribbon structures. These results highlight the important role of aggregation-induced pKa shifts in the self-assembly of model beta-sheet peptides.

Epitaxial growth of ultrathin palladium films on Re{0001}

Ultrathin bimetallic layers create unusual magnetic and surface chemical effects through the modification of electronic structure brought on by low dimensionality, polymorphism, reduced screening, and epitaxial strain. Previous studies have related valence and core-level shifts to surface reactivity through the d-band model of Hammer and Nørskov, and in heteroepitaxial films this band position is determined by competing effects of coordination, strain, and hybridization of substrate and overlayer states. In this study we employ the epitaxially matched Pd on Re{0001} system to grow films with no lateral strain. We use a recent advancement in low-energy electron diffraction to expand the data range sufficiently for a reliable determination of the growth sequence and out-of-plane surface relaxation as a function of film thickness. The results are supported by scanning tunneling microscopy and X-ray photoemission spectroscopy, which show that the growth is layer-by-layer with significant core-level shifts due to changes in film structure, morphology, and bonding.

Puroindoline-a, a lipid binding protein from common wheat, spontaneously forms prolate protein micelles in solution

The self-assembly in solution of puroindoline-a (Pin-a), an amphiphilic lipid binding protein from common wheat, was investigated by small angle neutron scattering, dynamic light scattering and size exclusion chromatography. Pin-a was found to form monodisperse prolate ellipsoidal micelles with a major axial radius of 112 +/- 4.5 A ˚ and minor axial radius of 40.4 +/- 0.18 A ˚ . These protein micelles were formed by the spontaneous self-assembly of 38 Pin-a molecules in solution and were stable over a wide pH range (3.5–11) and at elevated temperatures (20–65 degC). Pin-a micelles could be disrupted upon addition of the non-ionic surfactant dodecyl-b-maltoside, suggesting that the protein self-assembly is driven by hydrophobic forces, consisting of intermolecular interactions between Trp residues located within a well-defined Trp-rich domain of Pin-a.

Highly proton-ordered water structures on oxygen precovered Ru{0001}

We present helium scattering measurements of a water ad-layer grown on a O(2 1)/Ru(0001) surface. The adsorbed water layer results in a well ordered helium diffraction pattern with systematic extinctions of diffraction spots due to glide line symmetries. The data reflects a well-defined surface structure that maintains proton order even at surprisingly high temperatures of 140 K. The diffraction data we measure is consistent with a structure recently derived from STM measurements performed at 6 K. Comparison with recent DFT calculation is in partial agreement, suggesting that these calculations might be underestimating the contribution of relative water molecule orientations to the binding energy.

A study of the ethene-ozone reaction with photoelectron spectroscopy: measurement of product branching ratios and atmospheric implications

The ozone-ethene reaction has been investigated at low pressure in a flow-tube interfaced to a u.v. photoelectron spectrometer. Photoelectron spectra recorded as a function of reaction time have been used to estimate partial pressures of the reagents and products, using photoionization cross-sections for selected photoelectron bands of the reagents and products, which have been measured separately. Product yields compare favourably with results of other studies, and the production of oxygen and acetaldehyde have been measured as a function of time for the first time. A reaction scheme developed for the ozone-ethene reaction has been used to simulate the reagents and products as a function of time. The results obtained are in good agreement with the experimental measurements. For each of the observed products, the simulations allow the main reaction (or reactions) for production of that product to be established. The product yields have been used in a global model to estimate their global annual emissions in the atmosphere. Of particular interest are the calculated global annual emissions of formaldehyde (0.96 ± 0.10 Tg) and formic acid, (0.05 ± 0.01 Tg) which are estimated as 0.04% and 0.7% of the total annual emission respectively.

A kinetics and mechanistic study of the atmospherically relevant reaction between molecular chlorine and dimethyl sulfide (DMS)

A gas-phase kinetics study of the atmospherically important reaction between Cl2 and dimethyl sulfide (DMS)Cl2 + CH3SCH3 → productshas been made using a flow-tube interfaced to a photoelectron spectrometer. The rate constant for this reaction has been measured at 1.6 and 3.0 torr at T = (294 ± 2) K as (3.4 ± 0.7) × 10−14 cm3 molecule−1 s−1. Reaction (1) has been found to proceed via an intermediate, (CH3)2SCl2, to give CH3SCH2Cl and HCl as the products. The mechanism of this reaction and the structure of the intermediate were investigated using electronic structure calculations. A comparison of the mechanisms of the reactions between Cl atoms and DMS, and Cl2 and DMS has been made and the relevance of the results to atmospheric chemistry is discussed.

Modeling the tetraphenylalanine-PEG hybrid amphiphile: from DFT calculations on the peptide to molecular dynamics simulations on the conjugate

The conformational properties of the hybrid amphiphile formed by the conjugation of a hydrophobic peptide with four phenylalanine (Phe) residues and hydrophilic poly(ethylene glycol), have been investigated using quantum mechanical calculations and atomistic molecular dynamics simulations. The intrinsic conformational preferences of the peptide were examined using the building-up search procedure combined with B3LYP/ 6-31G(d) geometry optimizations, which led to the identification of 78, 78, and 92 minimum energy structures for the peptides containing one, two, and four Phe residues. These peptides tend to adopt regular organizations involving turn-like motifs that define ribbon or helicallike arrangements. Furthermore, calculations indicate that backbone ... side chain interactions involving the N-H of the amide groups and the pi clouds of the aromatic rings play a crucial role in Phe-containing peptides. On the other hand,MD simulations on the complete amphiphile in aqueous solution showed that the polymer fragment rapidly unfolds maximizing the contacts with the polar solvent, even though the hydrophobic peptide reduce the number of waters of hydration with respect to an individual polymer chain of equivalent molecular weight. In spite of the small effect of the peptide in the hydrodynamic properties of the polymer, we conclude that the two counterparts of the amphiphile tend to organize as independent modules.

Lipid binding interactions of antimicrobial plant seed defence proteins: puroindoline-a and β-purothionin

The indolines and thionins are basic, amphiphilic and cysteine-rich proteins found in cereals; puroindoline-a (Pin-a) and β-purothionin (β-Pth) are members of these families in wheat (Triticum aestivum). Pin-a and β-Pth have been suggested to play a significant role in seed defence against microbial pathogens, making the interaction of these proteins with model bacterial membranes an area of potential interest. We have examined the binding of these proteins to lipid monolayers composed of 1,2-dipalmitoyl-sn-glycero-3-phospho-(1'-rac-glycerol) (DPPG) using a combination of neutron reflectometry, Brewster angle microscopy, and infrared spectroscopy. Results showed that both Pin-a and β-Pth interact strongly with condensed phase DPPG monolayers, but the degree of penetration was different. β-Pth was shown to penetrate the lipid acyl chain region of the monolayer and remove lipids from the air/liquid interface during the adsorption process, suggesting this protein may be able to both form membrane spanning ion channels and remove membrane phospholipids in its lytic activity. Conversely, Pin-a was shown to interact mainly with the head-group region of the condensed phase DPPG monolayer and form a 33 Å thick layer below the lipid film. The differences between the interfacial structures formed by these two proteins may be related to the differing composition of the Pin-a and β-Pth hydrophobic regions.