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.

Constraint release in entangled binary blends of linear polymers: a molecular dynamics study

We present extensive molecular dynamics simulations of the dynamics of diluted long probe chains entangled with a matrix of shorter chains. The chain lengths of both components are above the entanglement strand length, and the ratio of their lengths is varied over a wide range to cover the crossover from the chain reptation regime to tube Rouse motion regime of the long probe chains. Reducing the matrix chain length results in a faster decay of the dynamic structure factor of the probe chains, in good agreement with recent neutron spin echo experiments. The diffusion of the long chains, measured by the mean square displacements of the monomers and the centers of mass of the chains, demonstrates a systematic speed-up relative to the pure reptation behavior expected for monodisperse melts of sufficiently long polymers. On the other hand, the diffusion of the matrix chains is only weakly perturbed by the diluted long probe chains. The simulation results are qualitatively consistent with the theoretical predictions based on constraint release Rouse model, but a detailed comparison reveals the existence of a broad distribution of the disentanglement rates, which is partly confirmed by an analysis of the packing and diffusion of the matrix chains in the tube region of the probe chains. A coarse-grained simulation model based on the tube Rouse motion model with incorporation of the probability distribution of the tube segment jump rates is developed and shows results qualitatively consistent with the fine scale molecular dynamics simulations. However, we observe a breakdown in the tube Rouse model when the short chain length is decreased to around N-S = 80, which is roughly 3.5 times the entanglement spacing N-e(P) = 23. The location of this transition may be sensitive to the chain bending potential used in our simulations.

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.

Reaction of a phospholipid monolayer with gas-phase ozone at the air-water interface: measurement of surface excess and surface pressure in real time

The reaction between gas-phase ozone and monolayers of the unsaturated lipid 1-palmitoy1-2-oleoyl-sn-glycero-3-phosphocholine, POPC, on aqueous solutions has been studied in real time using neutron reflection and surface pressure measurements. The reaction between ozone and lung surfactant, which contains POPC, leads to decreased pulmonary function, but little is known shout the changes that occur to the interfacial material as a result of oxidation. The results reveal that the initial reaction of ozone with POPC leads to a rapid increase in surface pressure followed by a slow decrease to very low values. The neutron reflection measurements, performed on an isotopologue of POPC with a selectively deuterated palmitoyl strand, reveal that the reaction leads to loss of this strand from the air-water interface. suggesting either solubilization of the product lipid or degradation of the palmitoyl strand by a reactive species. Reactions of H-1-POPC on D2O reveal that the headgroup region of the lipids in aqueous solution is not dramatically perturbed by the reaction of POPC monolayers with ozone supporting degradation of the palmitoyl strand rather than solubilization. The results are consistent with the reaction of ozone with the oleoyl strand of POPC at the air water interface leading to the formation of OH radicals. the highly reactive OH radicals produced can then go on to react with the saturated palmitoyl strands leading to the formation or oxidized lipids with shorter alkyl tails.

Temperature dependence of protein dynamics simulated with three different water models

The effect of variation of the water model on the temperature dependence of protein and hydration water dynamics is examined by performing molecular dynamics simulations of myoglobin with the TIP3P, TIP4P, and TIP5P water models and the CHARMM protein force field at temperatures between 20 and 300 K. The atomic mean-square displacements, solvent reorientational relaxation times, pair angular correlations between surface water molecules, and time-averaged structures of the protein are all found to be similar, and the protein dynamical transition is described almost indistinguishably for the three water potentials. The results provide evidence that for some purposes changing the water model in protein simulations without a loss of accuracy may be possible.

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.

Direct time-resolved study of the kinetics of the reaction of silylene with phenylsilane in the gas phase. Does SiH2 react with the aromatic ring?

Laser flash photolysis studies of silylene, SiH2, generated by the 193 nm laser flash photolysis phenylsilane, PhSiH3, have been carried out to obtain rate constants for its bimolecular reaction with PhSiH3 itself, in the gas phase. The reaction was studied in SF6 (mostly at 10 Torr total pressure) over the temperature range 298-595 K. The rate constants (also found to be pressure independent) gave the following Arrhenius equation: log(k/cm(3) molecule(-1) s(-1)) = (-9.92 +/- 0.04) + (3.31 +/- 0.27) kJ mol(-1)/RT ln 10 Similar investigations of the reaction of silylene with benzene, C6H6, (295-410 K) gave data suggestive of the fact that SiH2 might be reacting with photochemical products of C6H6 as well as with C6H6 itself. However, in the latter system, apparent rate constants were sufficiently low to indicate that in the reaction of SiH2 with PhSiH3 addition to the aromatic ring was unlikely to be in excess of 3% of the total. Quantum chemical calculations of the energy surface for SiH2 + C6H6 indicate that 7-silanorcaradiene and 7-silacycloheptatriene are possible products but that PhSiH3 formation is unlikely. RRKM calculations suggest that 7-silanorcaradiene should be the initial product but that it cannot be collisionally stabilized under experimental conditions

A unified titration formula

In a number of standard titrations, a volume of a monoprotic base MOH at a specific concentration is added to a volume of a monoprotic acid HA at a specific concentrations. Four different types of titration are possible, depending on whether the acid and base are strong or weak. A single unifying formula covering all four cases has been determined.

Construction of supramolecular helices and breaking the helicity by forming supramolecular β-sheet structures using suitable self-assembling pseudopeptide building blocks

Bis-valine derivatives or malonamide (Guha,S.; Drew, M.G.B. Small 2008, 4, 1993-2005) and a bis-valine derivative of 1,1-cyclopropone dicarboxamide were used as building blocks for the construction of supramolecular helical structures. The six-membered intramolecular hydrogen-bonded scaffold is formed, and this acts as a unique supramolecular synthon for the construction of a pseudopeptide-based supramolecular helical structure. However, in absence of this intramolecular hydrogen bond. intermolecular hydrogen bonds are formed among the peptide strands. This leads to a supramolecular beta-sheet structure. Proper selection of the supramolecular synthon (six-membered intramolecular hydrogenbonded scaffold) promotes supramolecular helix formation, and a deviation from this molecular structure dictates the disruption of supramolecular helicity. In this study, six crystal structures have been used to demonstrate that a change in the central angle and/or the central core structure of dicarboxamides can be used to design either a supramolecular helix or a beta-sheet.

Use of a reduced Schiff-Base ligand to prepare novel chloro-bridged chains of rare Cu(II) trinuclear complexes with mixed azido/oxo and chloro/oxo bridges

Two mixed bridged one-dimensional (1D) polynuclear complexes, [Cu3L2(mu(1,1)-N-3)(2)(mu-Cl)Cl](n) (1) and {[Cu3L2(mu-Cl)(3)Cl]center dot 0.46CH(3)OH}(n), (2), have been synthesized using the tridentate reduced Schiff-base ligand HL (2-[(2-dimethylamino-ethylamino)-methyl]-phenol). The complexes have been characterized by X-ray structural analyses and variable-temperature magnetic susceptibility measurements. In both complexes the basic trinuclear angular units are joined together by weak chloro bridges to form a 1D chain. The trinuclear structure of 1 is composed of two terminal square planar [Cu(L)(mu(1,1)-N-3)] units connected by a central Cu(II) atom through bridging nitrogen atoms of end-on azido ligands and the phenoxo oxygen atom of the tridentate ligand. These four coordinating atoms along with a chloride ion form a distorted trigonal bipyramidal geometry around the central Cu(II). The structure of 2 is similar; the only difference being a Cl bridge replacing the mu(1,1)-N-3 bridge in the trinuclear unit. The magnetic properties of both trinuclear complexes can be very well reproduced with a simple linear symmetrical trimer model (H = JS(i)S(i+1)) with only one intracluster exchange coupling (J) including a weak intertrimer interaction (.j) reproduced with the molecular field approximation. This model provides very satisfactory fits for both complexes in the whole temperature range with the following parameters: g = 2.136(3), J = 93.9(3) cm(-1) and zj= -0.90(3) cm(-1) (z = 2) for 1 and g = 2.073(7), J = -44.9(4) cm(-1) and zJ = -1.26(6) cm(-1) (z = 2) for 2.

A unique example of structural and magnetic diversity in four interconvertible copper(II)−azide complexes with the same Schiff base ligand: a monomer, a dimer, a chain, and a layer

Four new Cu(II)-azido complexes of formula [CuL(N-3)] (1), [CuL(N-3)](2) (2), [Cu7L2(N-3)(12)](n) (3), and [Cu2L(dmen)-(N-3)(3)](n) (4) (dmen = N,N-dimethylethylenediamine) have been synthesized using the same tridentate Schiff base ligand HL (2-[1-(2-dimethylaminoethylimino)ethyl]phenol, the condensation product of dmen and 2-hydroxyacetophenone). The four compounds have been characterized by X-ray structural analyses and variable-temperature magnetic susceptibility measurements. Complex 1 is mononuclear, whereas 2 is a single mu-1,1 azido-bridged dinuclear compound. The polymeric compound 3 possesses a 2D structure in which the Cu(II) ions are linked by phenoxo oxygen atoms and two different azide bridges (mu-1,1 and mu-1,1,3). The structure of complex 4 is a double helix in which two mu-1,3-azido-bridged alternating one-dimensional helical chains of CuL(N-3) and Cu(dmen)(N-3)(2) are joined together by weak mu-1,1 azido bridges and H-bonds. The complexes interconvert in solution and can be obtained in pure form by carefully controlling the conditions. The magnetic properties of compounds 1 and 2 show the presence of very weak antiferromagnetic exchange interactions mediated by a ligand pi overlap (J = -1.77) and by an asymmetric 1,1-N-3 bridge (J = -1.97 cm(-1)), respectively. Compound 3 presents, from the magnetic point of view, a decorated chain structure with both ferro- and antiferromagnetic interactions. Compound 4 is an alternating helicoidal chain with two weak antiferromagnetic exchange interactions (J -1.35 and -2.64 cm(-1)).

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.

Modeling the plant uptake of organic chemicals, including the soil-air-plant pathway

The soil−air−plant pathway is potentially important in the vegetative accumulation of organic pollutants from contaminated soils. While a number of qualitative frameworks exist for the prediction of plant accumulation of organic chemicals by this pathway, there are few quantitative models that incorporate this pathway. The aim of the present study was to produce a model that included this pathway and could quantify its contribution to the total plant contamination for a range of organic pollutants. A new model was developed from three submodels for the processes controlling plant contamination via this pathway: aerial deposition, soil volatilization, and systemic translocation. Using the combined model, the soil−air−plant pathway was predicted to account for a significant proportion of the total shoot contamination for those compounds with log KOA > 9 and log KAW < −3. For those pollutants with log KOA < 9 and log KAW > −3 there was a higher deposition of pollutant via the soil−air−plant pathway than for those chemicals with log KOA > 9 and log KAW < −3, but this was an insignificant proportion of the total shoot contamination because of the higher mobility of these compounds via the soil−root−shoot pathway. The incorporation of the soil−air−plant pathway into the plant uptake model did not significantly improve the prediction of the contamination of vegetation from polluted soils when compared across a range of studies. This was a result of the high variability between the experimental studies where the bioconcentration factors varied by 2 orders of magnitude at an equivalent log KOA. One potential reason for this is the background air concentration of the pollutants under study. It was found background air concentrations would dominate those from soil volatilization in many situations unless there was a soil hot spot of contamination, i.e., >100 mg kg−1.