Temperature-Dependent Implicit-Solvent Model of Polyethylene Glycol in Aqueous Solution

Concentration and temperature dependent interactions and state diagram of dispersions of copolymer microgels

We investigate by means of small angle neutron scattering experiments and numerical simulations the interactions and inter-particle arrangements of concentrated dispersions of copolymer poly(N-isopropylacrylamide)-poly(ethylene glycol methyl ether methacrylate) (PNIPAM-PEGMA) microgels across the volume phase transition (VPT). The scattering data of moderately concentrated dispersions are accurately modeled at all temperatures by using a star polymer form factor and static structure factors calculated from the effective potential obtained from simulations. Interestingly, for temperatures below the VPT temperature (VPTT), the radius of gyration and blob size of the particles significantly decrease with increasing the effective packing fraction in the non-overlapping regime. This is attributed to the presence of charges in the system associated with the use of an ionic initiator in the synthesis. Simulations using the experimentally corroborated interaction potential are used to explore the state diagram in a wide range of effective packing fractions. Below and slightly above the VPTT, the system undergoes an arrest transition mainly driven by the soft repulsion between the particles. Only well above the VPTT the system is found to phase separate before arresting. Our results highlight the versatility and potential of copolymer PNIPAM-PEGMA microgels to explore different kinds of arrested states balancing attraction and repulsion by changing temperature and packing fraction.

Coarse-Grained Model Incorporating Short- and Long-Range Effective Potentials for the Fast Simulation of Micelle Formation in Solutions of Ionic Surfactants

A coarse-grained model comprising short- and long-range effective potentials, parametrized with the iterative Boltzmann inversion (IBI) method, is presented for capturing micelle formation in aqueous solutions of ionic surfactants using as a model system sodium dodecyl sulfate (SDS). In the coarse-grained (CG) model, each SDS molecule is represented as a sequence of four beads while each water molecule is modeled as a single bead. The proposed CG scheme involves ten potential energy functions: four of them describe bonded interactions and control the distribution functions of intramolecular degrees of freedom (bond lengths, valence angles, and dihedrals) along an SDS molecule while the other six account for intermolecular interactions between pairs of SDS and water beads and control the radial distribution functions. The nonbonded effective potentials between coarse-grained SDS molecules extend up to about 12 nm and capture structural and morphological features of the micellar solution both at short and long distances. The long-range component of these potentials, in particular, captures correlations between surfactant molecules belonging to different micelles and is essential to describe ordering associated with micelle formation. A new strategy is introduced for determining the effective potentials through IBI by using information (target distribution functions) extracted from independent atomistic simulations of a micellar reference system (a salt-free SDS solution at total surfactant concentration c_T equal to 103 mM, temperature T equal to 300 K, and pressure P equal to 1 atm) obtained through a multiscale approach described in an earlier study. It employs several optimization steps for bonded and nonbonded interactions and a gradual parametrization of the short- and long-range components of the latter, followed by reparametrization of the bonded ones. The proposed CG model can reproduce remarkably accurately the microstructure and morphology of the reference system within only a few hours of computational time. It is therefore very promising for future studies of structural and morphological behavior of various liquid surfactant formulations.

Weakly hydrated anions bind to polymers but not monomers in aqueous solutions

Weakly hydrated anions help to solubilize hydrophobic macromolecules in aqueous solutions, but small molecules comprising the same chemical constituents precipitate out when exposed to these ions. Here, this apparent contradiction is resolved by systematically investigating the interactions of NaSCN with polyethylene oxide oligomers and polymers of varying molecular weight. A combination of spectroscopic and computational results reveals that SCN^- accumulates near the surface of polymers, but is excluded from monomers. This occurs because SCN^- preferentially binds to the centre of macromolecular chains, where the local water hydrogen-bonding network is disrupted. These findings suggest a link between ion-specific effects and theories addressing how hydrophobic hydration is modulated by the size and shape of a hydrophobic entity.

Solutions and Condensed Phases of PEG2000 from All-Atom Molecular Dynamics

Extensive all-atom molecular dynamics studies of polyethylene glycol (PEG₂₀₀₀) when solvated and in the polymer bulk condensed phases were performed across a wide temperature range. We proposed two modified all-atom force field and observed the fate of the PEG₂₀₀₀ macromolecule when solvated in water, water with 4% ethanol, and ethyl acetate. In aqueous solutions, the macromolecule collapsed into a prolate spheroidal ball-like structure while adopting a rather elongated coiled structure in ethyl acetate. Inspection of the polymer-condensed phases across the 150-340 K temperature range enabled the atomistic view of the solid glass below the glass transition temperature of 230 K < T_g < 250 K and the rubber behavior above T_g. Predicted properties include the enthalpy, density, and cohesive energy temperature behavior, the specific heat, thermal expansivity, thermal compressibility, bulk modulus, and Hildebrand solubility parameter both below and above T_g. Within the polymer matrix, the PEG₂₀₀₀ macromolecules were entangled displaying a wide distribution of sizes that persisted when transitioning from the glass to the rubbery phases. Calculated properties agree very well with experiments when available or stand as crucial predictions while awaiting experimental measurement. Understanding the thermodynamics and structure of this useful polymer enables the efficient prediction of its behavior when building novel composite materials for nanomedicine and nanotherapeutics.

Polyethylene Glycol Crowder's Effect on Enzyme Aggregation, Thermal Stability, and Residual Catalytic Activity

Protein stability and performance in various natural and artificial systems incorporating many other macromolecules for therapeutic, diagnostic, sensor, and biotechnological applications attract increasing interest with the expansion of these technologies. Here we address the catalytic activity of lysozyme protein (LYZ) in the presence of a polyethylene glycol (PEG) crowder in a broad range of concentrations and temperatures in aqueous solutions of two different molecular mass PEG samples (M_w = 3350 and 10000 g/mol). The phase behavior of PEG-protein solutions is examined by using dynamic light scattering (DLS) and small-angle X-ray scattering (SAXS), while the enzyme denaturing is monitored by using an activity assay (AS) and circular dichroism (CD) spectroscopy. Molecular dynamic (MD) simulations are used to illustrate the effect of PEG concentration on protein stability at high temperatures. The results demonstrate that LYZ residual activity after 1 h incubation at 80 °C is improved from 15% up to 55% with the addition of PEG. The improvement is attributed to two underlying mechanisms. (i) Primarily, the stabilizing effect is due to the suppression of the enzyme aggregation because of the stronger PEG-protein interactions caused by the increased hydrophobicity of PEG and lysozyme at elevated temperatures. (ii) The MD simulations showed that the addition of PEG to some degree stabilizes the secondary structures of the enzyme by delaying unfolding at elevated temperatures. The more pronounced effect is observed with an increase in PEG concentration. This trend is consistent with CD and AS experimental results, where the thermal stability is strengthened with increasing of PEG concentration and molecular mass. The results show that the highest stabilizing effect is approached at the critical overlap concentration of PEG.

Insights into the Thermal Response of a Poly(ethylene oxide)-poly(propylene oxide)-poly(ethylene oxide) Triblock Polymer in Water

A thermal responsive block copolymer made up of ethylene oxide (EO) and propylene oxide (PO) blocks was simulated with optimized atomistic potentials and enhanced sampling methods over a range of temperatures. The results for the L42 pluronic polymer (EO)₄(PO)₂₂(EO)₄, which is known to undergo a transition in this temperature range, and the similarly sized (EO)₃₀ polymer, which does not, are compared. The thermal responsive L42 polymers in a dilute solution tend to aggregate, and this tendency gets stronger as temperature increases. The poly(ethylene oxide) polymer shows no such tendency. The aggregation is stabilized by the hydrophobic contact of the propylene oxide methyl groups, which outweighs a small loss in hydrogen bonds between the ether oxygens and water.

Analysis of Atomistic Potentials for Poly(ethylene glycol) Ethers

Two different potentials, the modified TraPPE-UA model of Fischer, J. [ J. Phys. Chem. B 2008, 112, 2388-2398] and the modified general AMBER force-field (GAFF) model of Barbosa, N. S. V. [ J. Mol. Model. 2017, 23, 194]are tested for a variety of temperature-dependent properties for neat and aqueous solutions for small poly(oxyethylene) (PEO) oligomers and larger polymers. A set of charges for PEO of arbitrary size is adapted for the modified GAFF model. Both models accurately reproduce experimental properties, but we find that the modified TraPPE-UA model is more accurate.

On the origin of the extremely different solubilities of polyethers in water

The solubilities of polyethers are surprisingly counter-intuitive. The best-known example is the difference between polyethylene glycol ([-CH₂-CH₂-O-]_n) which is infinitely soluble, and polyoxymethylene ([-CH₂-O-]_n) which is completely insoluble in water, exactly the opposite of what one expects from the C/O ratios of these molecules. Similar anomalies exist for oligomeric and cyclic polyethers. To solve this apparent mystery, we use femtosecond vibrational and GHz dielectric spectroscopy with complementary ab initio calculations and molecular dynamics simulations. We find that the dynamics of water molecules solvating polyethers is fundamentally different depending on their C/O composition. The ab initio calculations and simulations show that this is not because of steric effects (as is commonly believed), but because the partial charge on the O atoms depends on the number of C atoms by which they are separated. Our results thus show that inductive effects can have a major impact on aqueous solubilities.

Tuning the collapse transition of weakly charged polymers by ion-specific screening and adsorption

The experimentally observed swelling and collapse response of weakly charged polymers to the addition of specific salts displays quite convoluted behavior that is not easy to categorize. Here we use a minimalistic implicit-solvent/explicit-salt simulation model with a focus on ion-specific interactions between ions and a single weakly charged polyelectrolyte to qualitatively explain the observed effects. In particular, we demonstrate ion-specific screening and bridging effects cause collapse at low salt concentrations whereas the same strong ion-specific direct interactions drive re-entrant swelling at high concentrations. Consistently with experiments, a distinct salt concentration at which the salting-out power of anions inverts from the reverse to direct Hofmeister series is observed. At this so called isospheric point, the ion-specific effects vanish. Furthermore, with additional simplifying assumptions, an ion-specific mean-field model is developed for the collapse transition which quantitatively agrees with the simulations. Our work demonstrates the sensitivity of the structural behavior of charged polymers to the addition of specific salt beyond simple screening and shall be useful for further guidance of experiments.

Machine-Learning Enabled New Insights into the Coil-to-Globule Transition of Thermosensitive Polymers Using a Coarse-Grained Model

We present a computational framework that integrates coarse-grained (CG) molecular dynamics (MD) simulations and a data-driven machine-learning (ML) method to gain insights into the conformations of polymers in solutions. We employ this framework to study conformational transition of a model thermosensitive polymer, poly( N-isopropylacrylamide) (PNIPAM). Here, we have developed the first of its kind, a temperature-independent CG model of PNIPAM that can accurately predict its experimental lower critical solution temperature (LCST) while retaining the tacticity in the presence of an explicit water model. The CG model was extensively validated by performing CG MD simulations with different initial conformations, varying the radius of gyration of chain, the chain length, and the angle between the adjacent monomers of the initial configuration of PNIPAM (total simulation time = 90 μs). Moreover, for the first time, we utilize the nonmetric multidimensional scaling (NMDS) method, a data-driven ML approach, to gain further insights into the mechanisms and pathways of this coil-to-globule transition by analyzing CG MD simulation trajectories. NMDS analysis provides entirely new insights and shows multiple metastable states of PNIPAM during its coil-to-globule transition above the LCST.

Coarse-Grained Simulations of Aqueous Thermoresponsive Polyethers

Thermoresponsive polymers can change structure or solubility as a function of temperature. Block co-polymers of polyethers have a response that depends on polymer molecular weight and co-polymer composition. A coarse-grained model for aqueous polyethers is developed and applied to polyethylene oxide and polyethylene oxide-polypropylene oxide-polyethylene oxide triblock co-polymers. In this model, no interaction sites on hydrogen atoms are included, no Coulombic interactions are present, and all interactions are short-ranged, treated with a combination of two- and three-body terms. Our simulations find that The triblock co-polymers tend to associate at temperatures above 350 K. The aggregation is stabilized by contact between The hydrophobic methyl groups on The propylene oxide monomers and involves a large, favorable change in entropy.