Molecular Dynamics Simulations of Concentrated and Dilute Aqueous Solutions of Poly(N-Isopropylacrylamide) Using a Modified OPLS-AA Model

We perform molecular dynamics (MD) simulations of aqueous poly(N-isopropylacrylamide) (PNIPAM) solutions using the OPLS-AA model and its modified version in combination with the TIP4P/2005 water model. The original OPLS-AA model underestimates the demixing temperature of concentrated solutions and the coil-globule transition temperature in dilute solutions, whereas the modified model reproduces both. It was reported that the transition from the globule state to the coil state never occurred in MD simulations with the original OPLS-AA model. This problem is avoided by using the modified model because of the fast dynamics near the coil-globule transition temperature, which is higher for the modified model than for the original model by ∼60 K. The distribution functions of two structural order parameters calculated from well equilibrated trajectories of the modified model clearly show that the coil-globule transition is not a discontinuous first-order transition, as suggested in previous MD studies, but a continuous transition as predicted from mean field theory. We also find that the effect of polymer-polymer hydrogen bonding on the dynamics of aggregated polymer chains is less significant than assumed in experimental studies.

A generic model for pH-sensitive collapse of hydrophobic polymers

The hydrophobic effect is an important contributor to the stability of proteins and may be influenced by many factors including the pH of the solution. To simplify the study of pH effects on proteins, we parameterize biologically motivated titratable monomers which we insert into the sequence of a hydrophobic polymer and study via constant pH molecular dynamics (MD) simulations. We calculate the potential of mean force of the polymer as a function of its radius of gyration at different pH values and observe that the collapsed state of the polymer is destabilized when the titratable monomer is more charged (high pH for an acid and low pH for a base). Further, the extent of the destabilization is influenced by the position of the titratable monomer along the polymer sequence. The pKa value of the titratable monomer is also observed to be sensitive to polymer conformation, in agreement with protein studies. We further study a zwitterionic polymer with an acidic and a basic monomer in the same sequence which presents a pH-dependent hairpin formation. Our model provides a simplified yet powerful framework to study pH effects on the hydrophobic effect, providing insights into mechanisms governing the behavior of intrinsically disordered proteins (IDPs) and pH-sensitive drug delivery, among other applications.

Anion Competition at Positively Charged Surfactant Monolayers

Interactions of anions with hydrophobic surfaces of proteins and water-soluble polymers depend on the ability of the ions to shed their hydration shells. At positively charged surfactant monolayers, the interactions of anions are less well understood. Due to the interplay of electrostatic surface forces, hydration effects, and ion-ion interactions in the electrostatic double layer, a comprehensive microscopic picture remains elusive. Herein, we study the interactions of chloride, bromide, and a mixture of these two anions at the aqueous interface of dihexadecyldimethylammonium (DHDA+) and dioctadecyldimethylammonium (DODA+) cationic monolayers. Using molecular dynamics simulations and three surface-sensitive X-ray scattering techniques, we demonstrate that bromide interacts preferentially over chloride with both monolayers. The structure of the two monolayers and their interfacial electron density profiles obtained from the simulations quantitatively reproduce the experimental data. We observe that chloride and bromide form contact ion pairs with the quaternary ammonium groups on both monolayers. However, ion pairing with bromide leads to a greater reduction in the number of water molecules hydrating the anion, resulting in more energetically stable ion pairs. This leads to long-range (>3 nm) lateral correlations between bromide ions on the structured DODA+ monolayer. These observations indicate that ion hydration is the dominant factor determining the interfacial electrolyte structure.

Interactions between Water and a Hydrophobic Polymer

To investigate the mechanisms of interactions between a hydrophobic polymer and water, molecular dynamics calculations and Raman spectroscopic measurements of cis-1,4-polyisoprene (PI)-water systems were performed. The results show that PI in water undergoes a coil-globule transition at around 248 K. The transition is attributed to changes in the density and diffusivity of water. The volume expansion of the supercooled liquid water induces the coil structure of PI. The phase separation of PI from water with an increase in the self-diffusion coefficient of water molecules results in the globule structure of PI. The self-diffusion coefficient of free water with PI is larger than that of pure water because PI has an effect to decrease the hydrogen-bonding strength of water. The result suggests that the effects of the coexisting water are important factors governing the physical and chemical properties of hydrophobic polymers.

Solvation Behavior of Elastin-like Polypeptides in Divalent Metal Salt Solutions

The effects of CaCl2 and MgCl2 on the cloud point temperature of two different elastin-like polypeptides (ELPs) were studied using a combination of cloud point measurements, molecular dynamics simulations, and infrared spectroscopy. Changes in the cloud point for the ELPs in aqueous divalent metal cation solutions were primarily governed by two competing interactions: the cation-amide oxygen electrostatic interaction and the hydration of the cation. In particular, Ca2+ cations can more readily shed their hydration shells and directly contact two amide oxygens by the formation of ion bridges. By contrast, Mg2+ cations were more strongly hydrated and preferred to partition toward the amide oxygens along with their hydration shells. In fact, although hydrophilic ELP V5A2G3 was salted-out at low concentrations of MgCl2, it was salted-in at higher salt concentrations. By contrast, CaCl2 salted the ELP sharply out of solution at higher salt concentrations because of the bridging effect.

Ions Speciation at the Water-Air Interface

In typical aqueous systems, including naturally occurring sweet and salt water and tap water, multiple ion species are co-solvated. At the water-air interface, these ions are known to affect the chemical reactivity, aerosol formation, climate, and water odor. Yet, the composition of ions at the water interface has remained enigmatic. Here, using surface-specific heterodyne-detected sum-frequency generation spectroscopy, we quantify the relative surface activity of two co-solvated ions in solution. We find that more hydrophobic ions are speciated to the interface due to the hydrophilic ions. Quantitative analysis shows that the interfacial hydrophobic ion population increases with decreasing interfacial hydrophilic ion population at the interface. Simulations show that the solvation energy difference between the ions and the intrinsic surface propensity of ions determine the extent of an ion's speciation by other ions. This mechanism provides a unified view of the speciation of monatomic and polyatomic ions at electrolyte solution interfaces.

Product from sessile droplet evaporation of PNIPAM/water system above LCST: A block or micro/nano-particles?

PNIPAM as a stimuli-responsive polymer has generated extreme interests due to its versatile applications. However, there is no research report on whether PNIPAM micro/nano-particles can be extracted from its suspension after phase separation. In the present work, LCST-type phase separation in self-synthesized PNIPAM/water system was investigated in depth by dividing the DLS testing process into four stages. In addition to quenching duration, temperature rise process, quenching temperature and PNIPAM concentration all have a great influence on particle size of the suspension. Meanwhile, the steady-state rheology and dynamic viscoelasticity results show that PNIPAM micro/nano-particles in the suspension are "soft" that can deform. Finally, FE-SEM was used to observe the morphology of dehydrated PNIPAM extracted by sessile droplet evaporation under different conditions. The results indicate that these "soft" particles are easier to fuse together, do not have sufficient mechanical strength to maintain their spherical morphology after dehydration. But the above fusion can be suppressed by adjusting evaporation conditions to acquire smaller PNIPAM particles which have sufficient mechanical properties to keep their basic particle morphology. Further, by changing evaporation pressure to positive or negative pressure, dehydrated PNIPAM micro/nano-particles with excellent uniformity and separation can be obtained. This work will provide guidance for extracting micro/nano-particles from polymer/diluent systems with LCST.

A simple theory for interfacial properties of dilute solutions

Recent studies suggest that cosolute mixtures may exert significant non-additive effects upon protein stability. The corresponding liquid–vapor interfaces may provide useful insight into these non-additive effects. Accordingly, in this work, we relate the interfacial properties of dilute multicomponent solutions to the interactions between solutes. We first derive a simple model for the surface excess of solutes in terms of thermodynamic observables. We then develop a lattice-based statistical mechanical perturbation theory to derive these observables from microscopic interactions. Rather than adopting a random mixing approximation, this dilute solution theory (DST) exactly treats solute–solute interactions to lowest order in perturbation theory. Although it cannot treat concentrated solutions, Monte Carlo (MC) simulations demonstrate that DST describes the interactions in dilute solutions with much greater accuracy than regular solution theory. Importantly, DST emphasizes a fundamental distinction between the “intrinsic” and “effective” preferences of solutes for interfaces. DST predicts that three classes of solutes can be distinguished by their intrinsic preference for interfaces. While the surface preference of strong depletants is relatively insensitive to interactions, the surface preference of strong surfactants can be modulated by interactions at the interface. Moreover, DST predicts that the surface preference of weak depletants and weak surfactants can be qualitatively inverted by interactions in the bulk. We also demonstrate that DST can be extended to treat surface polarization effects and to model experimental data. MC simulations validate the accuracy of DST predictions for lattice systems that correspond to molar concentrations.