Andrea
Pica
a and
Giuseppe
Graziano
*b
aEuropean Molecular Biology Laboratory, Grenoble Outstation, 71 Avenue des Martyrs, Grenoble, France
bDipartimento di Scienze e Tecnologie, Università del Sannio, Via Francesco de Sanctis snc, 82100 Benevento, Italy. E-mail: graziano@unisannio.it
First published on 29th November 2019
The T(collapse) of poly(N-isopropylacrylamide), PNIPAM, shows a nonlinear dependence on the concentration of NaSCN or NaClO4; in the case of NaClO4, for example, at very low concentrations of the salt, T(collapse) increases with the concentration, while it has an opposite trend at higher NaClO4 concentrations [J. Am. Chem. Soc., 2005, 127, 14505]. These puzzling experimental data can be rationalized by considering that low charge density and poorly hydrated ions, such as thiocyanate and perchlorate, interact preferentially with the surface of the polymer, and cause an increase of the magnitude of the energetic term that stabilizes swollen conformations at low salt concentrations. However, as both swollen and collapsed PNIPAM conformations are accessible to such ions in view of their large conformational freedom, the difference in the number of ions bound to PNIPAM surface upon collapse changes little on increasing the salt concentration. Thus, the energetic term that favors swollen conformations increases with salt concentration to a lesser extent than the solvent-excluded volume term (linked to the density increase caused by salt addition to water), that favors collapsed conformations, leading to a nonlinear trend of T(collapse).
In previous papers, we have already proposed that the increase in entropy is related to an increase in translational entropy of water molecules caused by the decrease in solvent-excluded volume associated with PNIPAM collapse (i.e., an expansion of the available configurational space).6,10 The solvent-excluded volume can reliably be measured by the solvent-accessible surface area of the molecule,11,12 that in the case of water is named WASA. Computer simulations on all-atom models have shown that the average WASA of the conformations populating the C-state is significantly larger than the average WASA of the conformations populating the G-state.13,14 This means that water molecules gain translational entropy on pushing PNIPAM chains to collapse. This entropic driving force has a geometric ground and can be captured by means of a theoretical approach that uses simple geometric models,6,10 and analytic relationships from classic scaled particle theory,15–17 SPT. Indeed, the solvent-excluded volume effect can be measured by calculating the reversible work associated with cavity creation in a liquid18 (including water and aqueous solutions), ΔGc. A fundamental point is that, on keeping fixed the van der Waals volume, VvdW, of the cavity and changing its shape, the ΔGc magnitude is markedly affected.11 Specifically, on increasing the cavity WASA (i.e., passing from a spherical cavity to prolate spherocylindrical cavities of increasing length), the ΔGc magnitude increases in an almost linear manner.11 Therefore, the solvent-excluded volume effect is not measured by the cavity VvdW, but by the cavity WASA. In addition: (a) the ΔGc magnitude is larger in water than in other liquids as a consequence of its larger number density, originating in the small size of water molecules and the strength of H-bonds;18–21 (b) the ΔGc magnitude in water increases with temperature because the density decreases to a very small extent over the 0–100 °C temperature range, due to the strength of H-bonds in comparison to the random thermal energy.18–21 According to this theoretical approach, the density of water and aqueous solutions and its temperature dependence play a fundamental role for the energetics of PNIPAM collapse. Its application has provided a coherent rationalization of several features of such transition: (a) the effect of several sodium salts on T(collapse);10 (b) the effect of urea, tetramethylurea and TMAO on T(collapse);6,22,23 (c) the co-non-solvency phenomenon in water–methanol solutions around room temperature.24,25
However, other experimental results have not been yet addressed in their entire complexity. Interesting results by Cremer and co-workers26,27 showed that: (a) NaSCN causes a small increase in T(collapse) at low concentration and, on further increasing the NaSCN concentration, the T(collapse) value remains practically constant; in other words, NaSCN addition causes a small stabilization of the C-state at low salt concentration, but it has no effect at 1 M concentration; (b) NaClO4 causes a small increase in T(collapse) at very low concentration and, on further increasing the NaClO4 concentration, the T(collapse) value decreases significantly; in other words, NaClO4 stabilizes the C-state at very low concentration, but stabilizes the G-state at high concentration. This means that, for both salts, the dependence of T(collapse) on salt concentration is non-linear, in contrast with the linear decrease caused by the addition of salts like NaCl, consisting of well hydrated ions (see Fig. 1). It is important to underscore that: (a) the NaSCN effect on PNIPAM T(collapse) was first detected by Schild and Tirrell;3 (b) Cremer and co-workers obtained qualitatively similar results also for elastin-like polypeptides,28 regardless of the chemical heterogeneity of the latter in comparison to PNIPAM chemical homogeneity.
![]() | ||
Fig. 1 Trend of PNIPAM T(collapse) as a function of salt concentration for NaCl, NaSCN and NaClO4, for polymer samples with a weight-average molecular weight of 170 kDa; data are from ref. 27. |
These findings are very interesting because: (a) both NaSCN and NaClO4 are denaturing agents of globular proteins29–31 (i.e., thiocyanate and perchlorate are at the far-right end of the Hofmeister anion series26); (b) PNIPAM is considered to be a simplified model of globular proteins;4,5 (c) the two salts, at high concentration, do not stabilize the C-state of PNIPAM, that should resemble the denatured state of globular proteins. In the present study, we analyze the effect of NaSCN and NaClO4 on PNIPAM collapse over a large salt concentration range by means of our theoretical approach, providing an interesting and reliable explanation.
ΔGtr = [ΔGc(C) − ΔGc(G)] − T·ΔSconf + [Ea(C) − Ea(G) + ΔE(intra)] = ΔΔGc − T·ΔSconf + ΔEa | (1) |
It is necessary to devise a reliable procedure to calculate the three terms in eqn (1). The ΔΔGc = ΔGc(C) − ΔGc(G) contribution is calculated by means of a simple geometric model:6,10 the G-state is a sphere, and the C-state is a prolate spherocylinder, having the same VvdW of the sphere representing the G-state, but a significantly larger WASA. As in previous applications, a PNIPAM chain of 40 monomers in the G-state is modelled as a sphere of radius a = 10 Å, VvdW = 4189 Å3 and WASA = 1633 Å2, whereas the C-state is modelled as a prolate spherocylinder of radius a = 5 Å, cylindrical length l = 46.7 Å, VvdW = 4189 Å3 and WASA = 2393 Å2 (experimental measurements34 show that the volume change associated with PNIPAM collapse in water is negligibly small). The geometric models for the G-state and the C-state (rough and rigid representations of the average geometric properties possessed by chains populating the two ensembles) are assumed to be not affected by salt addition to water (this assumption should work because, in water, the interior volume of G-state contains a lot of water molecules, which will gradually be replaced by anions on adding salt). Note that the ΔΔGc contribution: (a) is always positive because ΔGc increases with cavity WASA, even though the cavity VvdW is kept fixed;11 (b) is calculated by means of the classic SPT formulae for spherical and spherocylindrical cavities in a hard sphere fluid mixture16,17 (the pressure–volume term is neglected for its smallness at P = 1 atm). A fundamental variable is the volume packing density of the hard sphere fluid mixture (i.e., aqueous solutions) ξ3 = (π/6)·∑ρj·σj3, where ρj is the number density, in molecules per Å3, of the species j and σj is the corresponding hard sphere diameter; ξ3 represents the fraction of the total liquid volume actually occupied by ions and molecules. The physical reliability of classic SPT formulae is grounded in its geometric basis.11,35
Experimental values of the density of water,36 and the considered aqueous salt solutions,37 have been used to perform calculations over the 5–40 °C temperature range. The use of experimental density is important because it allows one to take into account, even though in an indirect way, the actual interactions existing among molecules and ions that determine the liquid density.10,20–24 The following effective hard sphere diameters have been used and considered to be temperature-independent: (a) σ(H2O) = 2.80 Å,38 corresponding to the position of the first maximum in the oxygen–oxygen radial distribution function of water,39 at room temperature and 1 atm; (b) σ(Na+) = 2.02 Å, which is two times the ionic radius first proposed by Wasastjerna;40 (c) σ(SCN−) = 3.94 Å, corresponding to the effective hard sphere diameter of carbon dioxide;41 (d) σ(ClO4−) = 4.80 Å, which is the Pauling-type diameter reported by Marcus.42
The T·ΔSconf term is calculated by considering that each monomer gains a temperature-independent conformational entropy upon swelling (i.e., the coupling between the conformational degrees of freedom of the backbone and those of the side chains is neglected):
T·ΔSconf = T·Nres·ΔSconf(res) | (2) |
In our previous studies,6,10,22–24 the ΔEa quantity was fixed to −166.5 kJ mol−1 in pure water, and was considered to be temperature-independent in view of the limited temperature range important for PNIPAM collapse. The ΔEa quantity was considered to be not affected by the addition to water of high charge density anions,10 because the latter have strong and preferential electrostatic attractions with water molecules. In contrast, the ΔEa quantity is expected to be larger in magnitude in aqueous solutions containing NaSCN or NaClO4, because the low charge density and poorly hydrated44,45 thiocyanate and perchlorate anions should have preferential attractions with PNIPAM isopropyl groups, as already pointed out by means of both experimental and computational studies.46 However, no theoretical relationship is used to arrive at the ΔEa values in NaSCN or NaClO4 aqueous solutions; the adopted procedure will become clear in the following.
![]() | ||
Fig. 2 Experimental density of water and aqueous 0.2, 0.4, 0.6, 0.8 and 1 M NaSCN solutions over the 5–40 °C temperature range and 1 atm (panel A); experimental density of water and aqueous 0.2, 0.4, 0.6, 0.8 and 1 M NaClO4 solutions over the 5–40 °C temperature range and 1 atm (panel B); data are from ref. 36 and 37. |
[NaSCN]/M | d(20 °C)/g ml−1 | T(collapse)/°C | ΔΔGc/kJ mol−1 | T·ΔSconf/kJ mol−1 | ΔEa/kJ mol−1 |
---|---|---|---|---|---|
0 | 0.9982 | 31.0 | 215.2 | 48.7 | −166.5 |
0.1 | 1.0025 | 31.7 | 215.7 | 48.8 | −166.9 |
0.2 | 1.0067 | 31.9 | 216.2 | 48.8 | −167.4 |
0.4 | 1.0150 | 32.0 | 217.0 | 48.8 | −168.2 |
0.6 | 1.0233 | 31.8 | 217.4 | 48.8 | −168.6 |
0.8 | 1.0315 | 31.5 | 218.1 | 48.7 | −169.4 |
1.0 | 1.0395 | 31.0 | 218.9 | 48.7 | −170.2 |
[NaClO4]/M | d(20 °C)/g ml−1 | T(collapse)/°C | ΔΔGc/kJ mol−1 | T·ΔSconf/kJ mol−1 | ΔEa/kJ mol−1 |
---|---|---|---|---|---|
0 | 0.9982 | 31.0 | 215.2 | 48.7 | −166.5 |
0.1 | 1.0060 | 31.0 | 216.1 | 48.7 | −167.4 |
0.2 | 1.0138 | 30.5 | 217.2 | 48.6 | −168.6 |
0.4 | 1.0293 | 29.0 | 219.1 | 48.3 | −170.8 |
0.6 | 1.0447 | 27.0 | 220.6 | 48.0 | −172.6 |
0.8 | 1.0601 | 24.5 | 222.0 | 47.6 | −174.4 |
1.0 | 1.0754 | 22.0 | 223.4 | 47.2 | −176.2 |
DSC measurements have shown that the enthalpy change associated with PNIPAM collapse transition, ΔHtr, decreases slightly on increasing the NaSCN concentration, and the process becomes less cooperative, as indicated by DSC peak broadening.53,54 It is worth noting that the macroscopic ΔHtr quantity cannot be considered an exact measure of the ΔEa quantity with a changed sign because the former accounts also for the energy change due to the water–water and water–ion structural reorganization associated with PNIPAM collapse.
Taking into account both the experimental T(collapse) values in the presence of the two salts and the grounds of the theoretical approach [i.e., the T·ΔSconf contribution is not affected by salt addition to water, and ΔGtr = 0 at T(collapse) because PNIPAM collapse is a first-order phase transition], it is possible to arrive at the ΔEa estimates that are listed in the last column of Tables 1 and 2, in the assumption that the classic SPT-ΔΔGc values are correct. Even though the ΔEa values are just estimates, not to be used for quantitative analysis, they can help to provide a qualitative but reliable picture. Looking at the numbers and keeping in mind that the positive ΔΔGc quantity rises significantly with temperature (see Fig. 4), one can state that: (a) on increasing the NaSCN concentration, the increase in the ΔEa magnitude is almost entirely counterbalanced by the ΔΔGc increase (note that the ΔHtr decrease measured via DSC53,54 does not contrast with the increase in the ΔEa magnitude); (b) on increasing the NaClO4 concentration, the increase in the ΔEa magnitude is overwhelmed by the ΔΔGc increase. The difference between SCN− and ClO4− ions mainly comes from the larger size of the latter and its better attractions with water molecules (i.e., the NaClO4 solubility in water is larger than that of NaSCN); these two factors lead to larger ξ3 values (see Fig. 3), and so to larger ΔΔGc values (see Fig. 4) that stabilize the G-state. In general, it should be recognized that such differences are not large because the PNIPAM conformational stability is ruled by a delicate balance of contrasting quantities [see eqn (1)].
At first sight, the emerged scenario does appear strange because PNIPAM is considered to be a simplified model of globular proteins, and NaSCN and NaClO4 are denaturing agents of the native state of globular proteins29–31 (i.e., in the latter case the destabilizing ΔEa quantity overwhelms the stabilizing ΔΔGc one on increasing the salt concentration). The point to recognize and underscore is that the G-state of PNIPAM does not resemble the native state of globular proteins. Quasi-elastic neutron scattering measurements,55 and MD simulations56,57 have shown that PNIPAM G-state is characterized by an ensemble of different globular conformations, containing in their interior volume a lot of water and cosolute molecules or ions. In other words, PNIPAM G-state does not resemble a “crystal molecule”58 in view of its large conformational freedom. This means that: (a) low charge density and highly polarizable ions, not having strong electrostatic attractions with water molecules, such as SCN− and ClO4−, are present inside globular conformations, where they have attractive interactions with PNIPAM moieties; (b) in the case of PNIPAM, the ΔEa magnitude is not large and does not increase significantly on increasing the salt concentration because the difference in accessibility of the two PNIPAM states to SCN− or ClO4− anions is smaller than expected; (c) this interpretation is supported by the experimental finding that T(collapse) nonlinearity versus NaSCN and NaClO4 concentration depends on PNIPAM molecular weight (see Fig. 7e and f in ref. 59); the structural organization of PNIPAM chains in the G-state should depend on their length and the accessibility to SCN− or ClO4− anions is expected to increase on increasing the molecular weight, in line with the measured, stronger decrease of T(collapse);59 (d) in the case of globular proteins, the ΔEa magnitude is large because there is a marked difference in accessibility between the native state and the denatured one to SCN− or ClO4− anions (i.e., the interior volume of native state proves to be inaccessible to water and cosolute molecules or ions since it is very well packed60). In other words, the important structural differences existing between the G-state of PNIPAM and the native state of globular proteins do not manifest themselves in the collapse thermodynamics in water, but in the response of T(collapse) to the addition of cosolutes or salts to water. In fact, the addition of urea, a well-known denaturing agent of globular proteins, stabilizes the G-state of PNIPAM,22 such as the addition of a significant NaClO4 quantity. The reliability of this analysis is also supported by the experimental finding that NaClO4 addition to water causes the stabilization of the molten globule state of cytochrome c,61 and staphylococcal nuclease,62 by recognizing that the molten globule state of globular proteins should somewhat resemble the G-state of PNIPAM.
The effect of NaSCN and NaClO4 addition on the collapse transition of PNIPAM has already been analyzed by two groups. On one side, Cremer and co-workers26–28,59 fitted their own T(collapse) experimental data by means of a heuristic relationship consisting of a constant, a linear, and a nonlinear term. The linear term should describe the G-state stabilization due to depletion effects of the added ions, whereas the nonlinear term has the shape of a Langmuir binding isotherm to describe the binding of low charge density ions to PNIPAM surface. The linear term was considered to account for effects related to the increase in surface tension caused by salt addition to water,63 and to the negative hydration entropy of the anions.64 Even though there is no clear theoretical ground behind their heuristic T(collapse) expression, Cremer and co-workers26 found that SCN− or ClO4− anions possess the largest values of binding constant to PNIPAM surface, in line with the present analysis. On the other side, Heyda and Dzubiella13 devised an approximate expression for T(collapse), resembling the one used by Cremer and co-workers, starting from a second order Taylor expansion of ΔGtr, considered to be a function of temperature and salt concentration, and used it to fit experimental data. The approximate T(collapse) expression of Heyda and Dzubiella is indeed very interesting (something similar has recently been used in a different context65), but unfortunately it cannot provide a molecular rationalization of the phenomenon because it is not grounded in statistical mechanics. For instance: (a) it makes direct use of the experimental value of the entropy change associated with PNIPAM collapse transition,13 that is a macroscopic quantity coming from diverse molecular contributions; (b) its linear term in salt concentration, the m-value, is positive for both NaSCN and NaClO4, the derivative (∂m/∂Csalt) proves to be negative for both salts, but markedly larger in magnitude for NaClO4, without any explanation proposed (look at Table 3 in ref. 13). This fitting exercise cannot be considered satisfactory to arrive at a complete understanding.
In conclusion, coherent application of our theoretical approach leads to the following rationalization. For very low concentrations of NaSCN or NaClO4, T(collapse) increases because the ΔEa magnitude increases to a greater extent than the ΔΔGc one. However, poorly hydrated ions, such as SCN− or ClO4−, interact so well with PNIPAM surface that they are bound in a significant manner to both the C-state and the G-state, and the difference in the number of ions bound to PNIPAM surface changes little on increasing the salt concentration. This implies that the ΔEa magnitude increases with salt concentration to a lesser extent than the ΔΔGc quantity, leading to a nonlinear trend of T(collapse). In other words, beyond a given salt concentration, ΔΔGc becomes dominant in the Gibbs energy balance. At this particular concentration point, the salt changes its qualitative action: it does not stabilize the C-state, but starts to stabilize the G-state, as in the case of sodium perchlorate.
This journal is © the Owner Societies 2020 |