Water-Polymer Coupling Induces a Dynamical Transition in Microgels

Hydration-induced dynamical changes in lyophilised and weakly hydrated apoferritin: insights from molecular dynamics simulation

The dynamics and functionality of proteins are significantly influenced by their interaction with water. For lyophilised (i.e. h ≤ 0.05 where h = g of H2O per g of protein) and weakly hydrated systems (i.e. h ≤ 0.38) hydration generally enhances protein mobility above the so-called 'dynamical transition' temperature (Td > 220 K). However, water-induced mobility hindrance at low temperatures (T < 175 K) has been reported in various proteins of varying secondary structure; namely green fluorescent protein (GFP), pig liver esterase, lysozyme, ribonuclease A (RNAse A) and apoferritin. By focussing on the dynamic behaviour of the apoferritin molecule, this study proposes mechanisms driving these hydration-induced mobility changes, particularly the less understood hindrance at low temperatures. Using atomistic molecular dynamics (MD) simulations of horse spleen apoferritin in the lyophilised (h = 0.05) and weakly hydrated (h = 0.31) states, we report here the impact of water on protein dynamics as a function of temperature. Through residue-specific mean squared displacement (MSD), radial distribution function (RDF), solvent accessible surface area (SASA), local hydration degree and hydrogen bonding analyses, we demonstrate that while water proximity directly correlates with mobility enhancement at high temperatures, the hydration-induced mobility reduction observed at temperatures below 175 K is primarily propagated through the protein backbone.

Water slowing down drives the occurrence of the low temperature dynamical transition in microgels

The protein dynamical transition marks an increase in atomic mobility and the onset of anharmonic motions at a critical temperature (T d), which is considered relevant for protein functionality. This phenomenon is ubiquitous, regardless of protein composition, structure and biological function and typically occurs at large protein content, to avoid water crystallization. Recently, a dynamical transition has also been reported in non-biological macromolecules, such as poly(N-isopropyl acrylamide) (PNIPAM) microgels, bearing many similarities to proteins. While the generality of this phenomenon is well-established, the role of water in the transition remains a subject of debate. In this study, we use atomistic molecular dynamics (MD) simulations and elastic incoherent neutron scattering (EINS) experiments with selective deuteration to investigate the microscopic origin of the dynamical transition and distinguish water and PNIPAM roles. While a standard analysis of EINS experiments would suggest that the dynamical transition occurs in PNIPAM and water at a similar temperature, simulations reveal a different perspective, also qualitatively supported by experiments. From room temperature down to about 180 K, PNIPAM exhibits only modest changes of dynamics, while water, being mainly hydration water under the probed extreme confinement, significantly slows down and undergoes a mode-coupling transition from diffusive to activated. Our findings therefore challenge the traditional view of the dynamical transition, demonstrating that it occurs in proximity of the water mode-coupling transition, shedding light on the intricate interplay between polymer and water dynamics.

Phenol release from pNIPAM hydrogels: scaling molecular dynamics simulations with dynamical density functional theory

We employed molecular dynamic simulations (MD) and the Bennett's acceptance ratio method to compute the free energy of transfer, ΔG_trans, of phenol, methane, and 5-fluorouracil (5-FU), between bulk water and water-pNIPAM mixtures of different polymer volume fractions, ϕ_p. For this purpose, we first calculate the solvation free energies in both media to obtain ΔG_trans. Phenol and 5-FU (a medication used to treat cancer) attach to the pNIPAM surface so that they show negative values of ΔG_trans irrespective of temperature (above or below the lower critical solution temperature of pNIPAM, T_c). Conversely, methane switches the ΔG_trans sign when considering temperatures below (positive) and above (negative) T_c. In all cases, and contrasting with some theoretical predictions, ΔG_trans maintains a linear behavior with the pNIPAM concentration up to large polymer densities. We have also employed MD to compute the diffusion coefficient, D, of phenol in water-pNIPAM mixtures as a function of ϕ_p in the diluted limit. Both ΔG_trans and D as a function of ϕ_p are required inputs to obtain the release halftime of hollow pNIPAM microgels through Dynamic Density Functional Theory (DDFT). Our scaling strategy captures the experimental value of 2200 s for 50 μm radius microgels with no cavity, for ϕ_p ≃ 0.83 at 315 K.

Adsorption Energetics of Amino Acid Analogs on Polymer/Water Interfaces Studied by All-Atom Molecular Dynamics Simulation and a Theory of Solutions

Energetics of adsorption was addressed with all-atom molecular dynamics simulation on the interfaces of poly(2-methoxyethyl acrylate) (PMEA), poly(methyl methacrylate) (PMMA), and poly(butyl acrylate) (PBA) with water. A wide variety of adsorbate solutes were examined, and the free energy of adsorption was computed with the method of energy representation. It was found that the adsorption free energy was favorable (negative) for all the combinations of solute and polymer, and among PMEA, PMMA, and PBA, the strongest adsorption was observed on PMMA for the hydrophobic solutes and on PMEA for the hydrophilic ones. According to the decomposition of the adsorption free energy into the contributions from polymer and water, it was seen that the polymer contribution is larger in magnitude with the solute size. The total free energy of adsorption was correlated well with the solvation free energy in bulk water only for hydrophobic solutes. The roles of the intermolecular interaction components such as electrostatic, van der Waals, and excluded-volume were further studied, and the electrostatic component was influential only in determining the polymer dependences of the adsorption propensities of hydrophilic solutes. The extent of adsorption was shown to be ranked by the van der Waals component in the solute-polymer interaction separately over the hydrophilic and hydrophobic solutes, with the excluded-volume effect from water pointed out to also drive the adsorption.

Modified Flory–Rehner Theory Describes Thermotropic Swelling Transition of Smart Copolymer Microgels

In the present article, we use an improved Flory–Rehner theory to describe the swelling behavior of copolymer microgels, where the interaction parameter is modeled by a Hill-like equation for a cooperative thermotropic transition. This description leads to very good fits of the swelling curves of the copolymer microgels at different comonomer contents (30 mol%, 50 mol% and 70 mol%) obtained by photon correlation spectroscopy. Fixed parameters, which are universally applicable for the respective monomers given in our previous work, are used to fit the swelling curves. The analysis of the swelling curves yields physically reasonable and meaningful results for the remaining adjustable parameters. The comonomer content of the statistical copolymer microgels poly(NNPAM-co-NIPAM), poly(NIPAM-co-NIPMAM) and poly(NIPMAM-co-NNPAM) is determined by nuclear magnetic resonance spectroscopy and is in agreement with the nominal comonomer feed used in the synthesis. To investigate the volume phase transition at a molecular level, swelling curves are also measured by Fourier transformation infrared spectroscopy. The obtained swelling curves are also fitted using the Hill-like model. The fits provide physically reasonable parameters too, consistent with the results from photon correlation spectroscopy.

Modeling Solution Behavior of Poly(N-isopropylacrylamide): A Comparison between Water Models

Water is known to play a fundamental role in determining the structure and functionality of macromolecules. The same crucial contribution is also found in the in silico description of polymer aqueous solutions. In this work, we exploit the widely investigated synthetic polymer poly(N-isopropylacrylamide) (PNIPAM) to understand the effect of the adopted water model on its solution behavior and to refine the computational setup. By means of atomistic molecular dynamics simulations, we perform a comparative study of PNIPAM aqueous solution using two advanced water models: TIP4P/2005 and TIP4P/Ice. The conformation and hydration features of an atactic 30-mer at infinite dilution are probed at a range of temperature and pressure suitable to detect the coil-to-globule transition and to map the P–T phase diagram. Although both water models can reproduce the temperature-induced coil-to-globule transition at atmospheric pressure and the polymer hydration enhancement that occurs with increasing pressure, the PNIPAM–TIP4P/Ice solution shows better agreement with experimental findings. This result can be attributed to a stronger interaction of TIP4P/Ice water with both hydrophilic and hydrophobic groups of PNIPAM, as well as to a less favorable contribution of the solvent entropy to the coil-to-globule transition.

Influence of hydration on segmental chain librations and dynamical transition in lipid bilayers

Continuous wave electron paramagnetic resonance spectroscopy of chain-labeled phospholipids is used to investigate the effects of hydration on the librational oscillations and the dynamical transition of phospholipid membranes in the low-temperature range 120-270 K. Bilayers of dipalmitoylphostatidiycholine (DPPC) spin-labeled at the first acyl chain segments and at the methyl ends and prepared at full, low, and very low hydration are considered. The segmental mean-square angular amplitudes of librations, 〈α²〉, are larger in the bilayer interior than at the polar/apolar interface and larger in the fully and low hydrated than in the very low hydrated membranes. For chain segments at the beginning of the hydrocarbon region, 〈α²〉-values are markedly restricted and temperature independent in DPPC with the lowest water content, whereas they increase with temperature in the low and fully hydrated bilayers, particularly at the highest temperatures. For chain segments at the chain termini, the librational amplitudes increase progressively, first slowly and then more rapidly with temperature in bilayers at any level of hydration. From the temperature dependence of the mean-square librational amplitude, the dynamical transition is detected around 240 K at the polar/apolar interface in fully and low hydrated DPPC and at around 225 K at the inner hydrocarbon region for bilayers at any hydration condition. At the dynamical transition the bilayers cross low energy barriers of activation energy in the range 10-20 kJ/mol. The results highlight biophysical properties of DPPC bilayers at low-temperature and provide evidence of the effects of the hydration on the dynamical transition in bilayers.

Dynamics of Water and Other Molecular Liquids Confined Within Voids and on Surface of Lignin Aggregates in Aging Bio Crude Oils

Neutron scattering methods were employed to study the microscopic structure and dynamics of Bio Crude Oils (BCOs) and their lignin fractions. The structure of the carbonaceous aggregates was investigated using Small Angle Neutron Scattering to reveal a fractal hierarchy as well as a growth of the aggregates as the aging of the BCO proceeds. Elastic Neutron Scattering measurements indicate that BCO liquid phase, comprised of water and other hydrogenated molecular liquids, is in a state of extreme confinement. Quasi-Elastic Neutron Scattering yields information on the molecular motions, indicating that long range translational diffusion is suppressed and only localized dynamics take place on the tens of picosecond time range. The obtained results provide quantitative information on the molecular activity, as aging proceed, in these reactive materials of relevance as potential renewable energy sources.

Thermoresponsivity of poly(N-isopropylacrylamide) microgels in water-trehalose solution and its relation to protein behavior

HYPOTHESES

Additives are commonly used to tune macromolecular conformational transitions. Among additives, trehalose is an excellent bioprotectant and among responsive polymers, PNIPAM is the most studied material. Nevertheless, their interaction mechanism so far has only been hinted without direct investigation, and, crucially, never elucidated in comparison to proteins. Detailed insights would help understand to what extent PNIPAM microgels can effectively be used as synthetic biomimetic materials, to reproduce and study, at the colloidal scale, isolated protein behavior and its sensitivity to interactions with specific cosolvents or cosolutes.

EXPERIMENTS

The effect of trehalose on the swelling behavior of PNIPAM microgels was monitored by dynamic light scattering; Raman spectroscopy and molecular dynamics simulations were used to explore changes of solvation and dynamics across the swelling-deswelling transition at the molecular scale.

FINDINGS

Strongly hydrated trehalose molecules develop water-mediated interactions with PNIPAM microgels, thereby preserving polymer hydration below and above the transition while drastically inhibiting local motions of the polymer and of its hydration shell. Our study, for the first time, demonstrates that slowdown of dynamics and preferential exclusion are the principal mechanisms governing trehalose effect on PNIPAM microgels, at odds with preferential adsorption of alcohols, but in full analogy with the behavior observed in trehalose-protein systems.

Mobility of bound water in PNIPAM microgels

Polymer-solvent interactions play a crucial role in the stimuli-responsive behaviour of polymer networks. They influence the swelling/deswelling behaviour as well as the dynamics of the polymer chains. Scattering experiments provide insight into the polymer-water interaction of poly(N-isopropylacrylamide) (PNIPAM) microgels cross-linked with N,N'-methylenebisacrylamide (BIS) in dried and humidified state. The water mobility is studied by means of neutron spin-echo spectroscopy and neutron backscattering spectroscopy. The residual water amount has been determined with Karl Fischer titration. For both degrees of humidification, the relaxation time of the water molecules is much larger than that of free water due to the strong interactions with the polymer network and is only weakly depending on temperature and length scale of observation. The possible influence of the water on methyl group rotations is discussed.

Controlling Growth of Poly (Triethylene Glycol Acrylate-Co-Spiropyran Acrylate) Copolymer Liquid Films on a Hydrophilic Surface by Light and Temperature

A quartz crystal microbalance with dissipation monitoring (QCM-D) was employed for in situ investigations of the effect of temperature and light on the conformational changes of a poly (triethylene glycol acrylate-co-spiropyran acrylate) (P (TEGA-co-SPA)) copolymer containing 12-14% of spiropyran at the silica-water interface. By monitoring shifts in resonance frequency and in acoustic dissipation as a function of temperature and illumination conditions, we investigated the evolution of viscoelastic properties of the P (TEGA-co-SPA)-rich wetting layer growing on the sensor, from which we deduced the characteristic coil-to-globule transition temperature, corresponding to the lower critical solution temperature (LCST) of the PTEGA part. We show that the coil-to-globule transition of the adsorbed copolymer being exposed to visible or UV light shifts to lower LCST as compared to the bulk solution: the transition temperature determined acoustically on the surface is 4 to 8 K lower than the cloud point temperature reported by UV/VIS spectroscopy in aqueous solution. We attribute our findings to non-equilibrium effects caused by confinement of the copolymer chains on the surface. Thermal stimuli and light can be used to manipulate the film formation process and the film's conformational state, which affects its subsequent response behavior.

Effect of trehalose on protein cryoprotection: Insights into the mechanism of slowing down of hydration water

We study, with molecular dynamics simulations, a lysozyme protein immersed in a water-trehalose solution upon cooling. The aim is to understand the cryoprotectant role played by this disaccharide through the modifications that it induces on the slow dynamics of protein hydration water with its presence. The α-relaxation shows a fragile to strong crossover about 20° higher than that in the bulk water phase and 15° higher than that in lysozyme hydration water without trehalose. The protein hydration water without trehalose was found to show a second slower relaxation exhibiting a strong to strong crossover coupled with the protein dynamical transition. This slower relaxation time importantly appears enormously slowed down in our cryoprotectant solution. On the other hand, this long-relaxation in the presence of trehalose is also connected with a stronger damping of the protein structural fluctuations than that found when the protein is in contact with the pure hydration water. Therefore, this appears to be the mechanism through which trehalose manifests its cryoprotecting function.

Molecular insights on poly(N-isopropylacrylamide) coil-to-globule transition induced by pressure

By using extensive all-atom molecular dynamics simulations of an atactic linear polymer chain, we provide microscopic insights into poly(N-isopropylacrylamide) (PNIPAM) coil-to-globule transition addressing the roles played by both temperature and pressure. We detect a coil-to-globule transition up to large pressures, showing a reentrant behavior of the critical temperature with increasing pressure in agreement with experimental observations. Furthermore, again confirming the experimental findings, we report the existence at high pressures of a new kind of globular state. It is characterized by a more structured hydration shell that is closer to PNIPAM hydrophobic domains, as compared to the globular state observed at atmospheric pressure. Our results highlight that temperature and pressure induce a PNIPAM coil-to-globule transition through different molecular mechanisms, opening the way for a systematic use of both thermodynamic variables to tune the location of the transition and the properties of the associated swollen/collapsed states.

Atomic scale investigation of the volume phase transition in concentrated PNIPAM microgels

Combining elastic incoherent neutron scattering and differential scanning calorimetry, we investigate the occurrence of the volume phase transition (VPT) in very concentrated poly-(N-isopropyl-acrylamide) (PNIPAM) microgel suspensions, from a polymer weight fraction of 30 wt. % up to dry conditions. Although samples are arrested at the macroscopic scale, atomic degrees of freedom are equilibrated and can be probed in a reproducible way. A clear signature of the VPT is present as a sharp drop in the mean square displacement of PNIPAM hydrogen atoms obtained by neutron scattering. As a function of concentration, the VPT gets smoother as dry conditions are approached, whereas the VPT temperature shows a minimum at about 43 wt. %. This behavior is qualitatively confirmed by calorimetry measurements. Molecular dynamics simulations are employed to complement experimental results and gain further insights into the nature of the VPT, confirming that it involves the formation of an attractive gel state between the microgels. Overall, these results provide evidence that the VPT in PNIPAM-based systems can be detected at different time- and length-scales as well as under overcrowded conditions.

Solvation dynamics of N-substituted acrylamide polymers and the importance for phase transition behavior

Functional groups present in thermo-responsive polymers are known to play an important role in aqueous solutions by manifesting their coil-to-globule conformational transition in a specific temperature range. Understanding the role of these functional groups and their interactions with water is of great interest as it may allow us to control both the nature and temperature of this coil-to-globule transition. In this work, polyacrylamide (PAAm), poly(N-isopropylacrylamide) (PNIPAm), and poly(N-isopropylmethacrylamide) (PNIPMAm) solvated in water are studied with the goal of discovering the structure of the solvent and its interaction with these polymers in determining the polymer conformations. Specifically, all-atom molecular dynamics (MD) simulations were performed on polymer chains with 30 monomer units (30-mers) at 295 K, 310 K and 320 K, which is below and above the lower critical solution temperature (LCST) of PNIPAm (LCST = 305 K) and PNIPMAm (LCST = 315 K), respectively. The MD simulation trajectories suggest that changes in the functional groups in the backbone and side-chains alter the water solvation shell around the polymer. This results in a change in the residence time probability and hydrogen bond characteristics of water at simulated temperatures. Specifically, water molecules reside for longer times near PAAm (no LCST) and PNIPMAm (LCST = 315 K) chains as compared to PNIPAm. This might be one of the possible causes for the higher LCST of PNIPMAm as compared to that of PNIPAm. These results can guide experimentalists and theoreticians to design new polymer structures with tailor-made LCST transitions while controlling the water solvation shell around the functional group.

Coil-Globule Transition Thermodynamics of Poly(N-isopropylacrylamide)

Thermosensitive polymers such as poly(N-isopropylacrylamide) (PNIPAM) undergo a phase transition in aqueous solution from a random-coil structural ensemble to a globule structural ensemble at the lower critical solution temperature (LCST). Above this temperature, PNIPAM agglomerates and becomes insoluble, whereas it is soluble below the temperature. Thus, thermosensitive polymers represent essential targets for several applications, e.g., in drug delivery. Although their ability to change structure in response to a temperature alteration is highly relevant for industrial processes, their thermodynamic properties are mostly qualitatively understood, and the quantitative thermodynamic picture is still elusive. In this study, we used a combined atomistic molecular dynamics and well-tempered metadynamics simulation approach to estimate coil-globule transition thermodynamics. An isotactic 30-mer of PNIPAM was investigated over a broad temperature range between 200 and 360 K. The transition from the globule to the random-coil structure was observed with well-tempered metadynamics. For the first time, the free energy surface of PNIPAM was estimated and it is shown that the simulation results are in line with the experimentally observed thermosensitive behavior. Below the LCST, the random-coil ensemble represents the global energy minimum and is thermodynamically favored by 21 ± 9 kJ/mol compared to the globule ensemble; both are separated by a barrier of 49 ± 14 kJ/mol. In contrast, above the LCST, the globule ensemble is thermodynamically favored by 21 ± 8 kJ/mol over the random-coil ensemble. The barrier from random-coil to globule is 17 ± 10 kJ/mol.

Unsaturated lipid bilayers at cryogenic temperature: librational dynamics of chain-labeled lipids from pulsed and CW-EPR

Fully hydrated bilayers of monounsaturated palmitoyloleoylphosphatidylcholine (POPC) and diunsaturated dioleoylphosphatidylcholine (DOPC) lipids have low main phase transition temperatures (271 K for POPC and 253 K for DOPC). Two-pulse echo detected spectra, combined with continuous wave electron paramagnetic resonance spectroscopy, are employed to study the low-temperature lamellar phases of the POPC and DOPC unsaturated bilayers that are usually studied in the fluid state. Phosphatidylcholine spin-labeled at C-5 and C-16 carbon atom positions along the acyl chain were used and the temperature varied over the range 77-270 K. Segmental chain librational oscillations of small amplitude and with correlation time in the subnanosecond to nanosecond range are found in both membranes. The mean-square angular amplitude, α2, of librations increases with temperature, is larger close to the bilayer midplane than close to the first acyl chain segments, and is larger in diunsaturated than in monounsaturated bilayers. In the inner hydrocarbon region of both lipid matrices, α2 increases first slowly and linearly with temperature and then more rapidly, and a dynamical transition is detected in the range 190-210 K. Compared to dipalmitoylphosphatidylcholine bilayers of fully saturated symmetric chain lipids, the presence of double bonds in the acyl chain enhances the intensity of librational motion which is characterized by larger angular variations at the terminal methyl ends. These findings highlight biophysical properties of unsaturated bilayers in the frozen state, including a detailed characterization of segmental chain dynamics and the evidence of a dynamical transition that appears to be a generic feature in hydrated macromolecular systems. These results can also be relevant in regulating membrane physical properties and function at higher physiological temperatures.

Characterization of hydration water in supercooled water-trehalose solutions: The role of the hydrogen bonds network

The structural and dynamical properties of hydration water in aqueous solutions of trehalose are studied with molecular dynamics simulation. We simulate the systems in the supercooled region to investigate how the interaction with the trehalose molecules modifies the hydrogen bond network, the structural relaxation, and the diffusion properties of hydration water. The analysis is performed by considering the radial distribution functions, the residence time of water molecules in the hydration shell, the two body excess entropy, and the hydrogen bond water-water and water-trehalose correlations of the hydration water. The study of the two body excess entropy shows the presence of a fragile to strong crossover in supercooled hydration water also found in the relaxation time of the water-water hydrogen bond correlation function, and this is in agreement with predictions of the mode coupling theory and of previous studies of the oxygen-oxygen density correlators [A. Iorio et al., J. Mol. Liq. 282, 617 (2019); Sci. China: Phys., Mech. Astron. 62, 107011 (2019)]. The water-trehalose hydrogen bond correlation function instead evidences a strong to strong crossover in the relaxation time, and this crossover is related to a trehalose dynamical transition. This signals the role that the strong interplay between the soluted molecules and the surrounding solvent has in determining the dynamical transition common to both components of the system that happens upon cooling and that is similar to the well known protein dynamical transition. We connect our results with the cryoprotecting role of trehalose molecules.