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.

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.

Supramolecular structures of enzyme clusters

The structural characterization of subtilisin mesoscale clusters, which were previously shown to induce supramolecular order in biocatalytic self-assembly of Fmocdipeptides, was carried out by synchrotron small-angle X-ray, dynamic, and static light scattering measurements. Subtilisin molecules self-assemble to form supramolecular structures in phosphate buffer solutions. Structural arrangement of subtilisin clusters at 55 degrees Centigrade was found to vary systematically with increasing enzyme concentration. Static light scattering measurements showed the cluster structure to be consistent with a fractal-like arrangement, with fractal dimension varying from 1.8 to 2.6 with increasing concentration for low to moderate enzyme concentrations. This was followed by a structural transition around the enzyme concentration of 0.5 mg mL-1 to more compact structures with significantly slower relaxation dynamics, as evidenced by dynamic light scattering measurements. These concentration-dependent supramolecular enzyme clusters provide tunable templates for biocatalytic self-assembly.

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.

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.

Molecular dynamics simulation insight into the mechanism of phenol adsorption at low coverages from aqueous solutions on microporous carbons

MD simulation studies showing the influence of porosity and carbon surface oxidation on phenol adsorption from aqueous solutions on carbons are reported. Based on a realistic model of activated carbon, three carbon structures with gradually changed microporosity were created. Next, a different number of surface oxygen groups was introduced. The pores with diameters around 0.6 nm are optimal for phenol adsorption and after the introduction of surface oxygen functionalities, adsorption of phenol decreases (in accordance with experimental data) for all studied models. This decrease is caused by a pore blocking effect due to the saturation of surface oxygen groups by highly hydrogen-bounded water molecules.

Adsorption from aqueous solutions on opened carbon nanotubes: organic compounds speed up delivery of water from inside

We report the results of first systematic studies of organic adsorption from aqueous solutions onto relatively long single walled carbon nanotubes (four tubes, in initial and oxidised forms). Using molecular dynamics simulations (GROMACS package) we discuss the behaviour of tube-water as well as tube-adsorbate systems, for three different adsorbates (benzene, phenol and paracetamol).

Radiative efficiencies for fluorinated esters: indirect global warming potentials of hydrofluoroethers

Density Functional Theory (DFT) has been used with an empirically-derived correction for the wavenumbers of vibrational band positions to predict the infrared spectra of several fluorinated esters (FESs). Radiative efficiencies (REs) were then determined using the method of Pinnock et al. and these were used with atmospheric lifetimes from the literature to determine the direct global warming potentials of FESs. FESs, in particular fluoroalkylacetates, alkylfluoroacetates and fluoroalkylformates, are potential greenhouse gases and their likely long atmospheric lifetimes and relatively large REs, compared to their parent HFEs, make them active contributors to global warming. Here, we use the concept of indirect global warming potential (indirect GWP) to assess the contribution to the warming of several commonly used HFEs emitted from the Earth's surface, explicitly taking into account that these HFEs will be converted into the corresponding FESs in the troposphere. The indirect GWP can be calculated using the radiative efficiencies and lifetimes of the HFE and its degradation FES products. We found that the GWPs of those studied HFEs which have the smallest direct GWP can be increased by 100-1600% when taking account of the cumulative effect due to the secondary FESs formed during HFE atmospheric oxidation. This effect may be particularly important for non-segregated HFEs and some segregated HFEs, which may contribute significantly more to global warming than can be concluded from examination of their direct GWPs.