How Does a Hydrophobic Macromolecule ... - ACS Publications

Oct 4, 2016 - Indrajit Tah and Jagannath Mondal*. Tata Institute of Fundamental Research, Center for Interdisciplinary Sciences, 21 Brundavan Colony, ...
0 downloads 0 Views 3MB Size
Subscriber access provided by CORNELL UNIVERSITY LIBRARY

Article

How Does a Hydrophobic Macromolecule Respond to Mixed Osmolyte Environment? Indrajit Tah, and Jagannath Mondal J. Phys. Chem. B, Just Accepted Manuscript • DOI: 10.1021/acs.jpcb.6b08378 • Publication Date (Web): 04 Oct 2016 Downloaded from http://pubs.acs.org on October 7, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry B is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

How Does a Hydrophobic Macromolecule Respond to Mixed Osmolyte Environment? Indrajit Tah and Jagannath Mondal∗ Tata Institute of Fundamental Research, Center for Interdisciplinary sciences, 21 Brundavan Colony, Narsingi, Hyderabad, India E-mail: [email protected],+914024195021 Abstract The role of the protecting osmolyte Trimethyl N-oxide (TMAO) in counteracting the denaturing effect of urea on a protein is quite well established. However, the mechanistic role of osmolytes on the hydrophobic interaction underlying protein folding is a topic of contention and is emerging as a key area of biophysical interest. Although recent experiment and computer simulation have established that individual aqueous solution of TMAO and urea respectively stabilizes and destabilizes the collapsed conformation of a hydrophobic polymer, it remains to be explored how a mixed aqueous solution of protecting and denaturing osmolytes influences the conformations of the polymer. In order to bridge the gap, we have simulated the conformational behavior of both a model hydrophobic polymer and a synthetic polymer polystyrene in an aqueous mixture of TMAO and urea. Intriguingly, our free energy based simulations on both the systems show that even though a pure aqueous solution of TMAO stabilizes the collapsed or globular conformation of the hydrophobic polymer, addition of TMAO to an aqueous solution of urea further destabilizes the collapsed conformation of the hydrophobic polymer. We also observe that the extent of destabilization in a mixed osmolyte solution is relatively higher than that in pure aqueous urea solution.

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The reinforcement of the denaturation effect of the hydrophobic macromolecule in a mixed osmolyte solution is in stark contrast to the well-known counteracting role of TMAO in proteins under denaturing condition of urea. In both model and realistic systems, our results show that in a mixed aqueous solution, greater number of cosolutes preferentially bind to the extended conformation of the polymer relative to that in the collapsed conformation, thereby complying with Tanford-Wyman preferential solvation theory disfavoring the collapsed conformation. The results are robust across a range of osmolyte concentrations and multiple cosolute forcefields. Our findings unequivocally imply that the action of mixed osmolyte solution on hydrophobic polymer is significantly distinct from that of proteins.

2 ACS Paragon Plus Environment

Page 2 of 28

Page 3 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Introduction Various organisms utilize a class of small aqueous cosolutes known as osmolytes 1–8 to cope up with the osmotic stress. Trimethyl N-oxide (TMAO) 9 is one of the natural osmolytes found in various marine animals. As a protecting osmolyte, TMAO is known to stabilize the conformational structure of proteins by counteracting the action of a denaturing osmolyte, such as urea. 1,10–13 However, the role of osmolytes on the hydrophobic interaction, the major driving force underlying protein conformation, have only started to garner significant attention recently. Initial works on action of osmolytes on hydrophobic interaction have focussed on understanding, via simulation and experiments, the effect of an aqueous osmolyte solution of a hydrophobic macromolecule. 14–18 The prevalent findings suggest that, the globular or so-called collapsed conformation of the macromolecule is stabilized at a low 19 TMAO concentration and destabilized at high urea concentration. However, the question remains how an aqueous mixture of protecting and denaturing osmolyte can influence the collapse behavior of a hydrophobic polymer. Here we show, using extensive free-energy-based computer simulation of a model polymer and a synthetic polymer polystyrene, that while the individual solution of TMAO and urea respectively stabilize and destabilize the collapsed conformation of a polymer, a mixed environment of both the osmolytes further destabilizes the collapsed conformation, even compared to that in an individual aqueous solution of urea. These observations are in line with the seminal theories of preferential binding by Tanford and Wyman and imply that the action of mixed osmolytes on hydrophobic macromolecules are distinctly different from that of a protein. Over the years, numerous experimental 20–29 and theoretical studies 14–16,30–32 have focussed majorly on proteins and sparingly on other macromolecules in single or mixedcomponent osmolyte solutions involving TMAO as a protecting osmolyte and urea as a denaturing osmolyte, resulting in hypothesis of different direct/indirect mechanisms and contradictory conclusions. In case of proteins, it is now established that TMAO counteracts the denaturing effect of urea towards protein but the mechanistic details have differed 3 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

in many studies. Both indirect osmolyte-water induced mechanism 33–36 and direct proteinosmolyte induced interaction 14,30–32,37–39 have been hypothesized by experiments and theory. One of the most prevalent hypotheses in the counteraction role of TMAO has been the preferential exclusion of TMAO from the surface of the protein 37,38,40 and the preferential binding of urea on the surface of protein. Overall, majority of the aforementioned works on proteins in a mixed solvent environment have suggested that the two opposing osmolytes’ actions are independent of each other and counter-active in nature. However, on the contrary, Ganguly et al 32 has shown that urea and TMAO are mutually excluded from protein surfaces, accompanied with an increase of osmolyte-osmolyte aggregation, which points to a non-additive, mixed-osmolyte effect. Proteins are generally comprised of heterogenous amino acids and hence isolating individual contributions to a solvent induced change in stability is a difficult task. Towards this end, homopolymers are better suited, which gives an incentive for exploring the action of osmolytes on hydrophobic interaction, one of the central driving forces for protein folding. Initial works on action of osmolytes on hydrophobic interaction had focussed on simulating an aqueous osmolyte solution of homopolymeric chain whose hydrophobicity can be controlled by tuning the dispersion interactions. 14–18 Recently, one of us 30 have systematically explored the collapse behavior of a model Lennard-Jones homopolymers with a controllable dispersion interaction in aqueous solutions of TMAO and urea. The work by Mondal et al 30 had shown that, at a high polymer dispersion interaction and low TMAO concentration, 19 aqueous TMAO acts to stabilize and aqueous urea acts to destabilize the globular structure of the model polymers, as observed in proteins. However, it was found that both urea and TMAO strongly bind to the polymer surface, in contrast to TMAO’s so-called exclusion mechanism from the protein surface. In compliance with the classic theories of Tanford and Wyman, 2,3 Mondal et al’s 30 study highlighted that it is actually the relative difference in preferential binding of osmolytes between collapsed and extended state, rather than the absolute magnitude of preferential binding in either state, which needs to be considered in

4 ACS Paragon Plus Environment

Page 4 of 28

Page 5 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

a proper treatment of osmolyte mediated polymer collapse. This work later stimulated a collaboration 31 of single molecule force spectroscopy experiment by Walker group and MD simulation by Berne group to investigate the effect of TMAO and urea on a synthetic hydrophobic polymer polystyrene, which lent credence to the predictions made by the Mondal et al 30 for the model polymer. However, most of the reported works investigating the role of osmolytes on hydrophobic macromolecules have either studied the action of pure aqueous solution of TMAO or pure aqueous solution of urea on the hydrophobic polymer. But, how would a hydrophobic macromolecule conformationally behave in a mixed aqueous solution of TMAO and urea? To bridge this gap, in the present work, we computationally investigate the behavior of a single hydrophobic polymer chain in a mixed aqueous solution of TMAO and urea. Specifically, we have simulated the conformational behavior of both the model hydrophobic polymer and synthetic hydrophobic polymer polystyrene,, 30,31 in a mixed aqueous environment of TMAO and urea and compared the results with that in an aqueous media of pure TMAO or pure urea. Our results are consistent with Wyman-Tanford theory 2,3 of preferential binding of cosolute over solvent towards macromolecules. The robustness of the results has been validated across multiple combination of osmolyte-mixtures of varying concentration, variable polymer-dispersion interactions and the variation of TMAO forcefield. The remainder of the manuscript is organized as follows: The methods section describes the details of molecular models and computational methodology in two independent subsections. The results section is divided into three subsections. The first subsection of results describes our simulated result on action of mixtures of osmolytes on a model polymer and tests the robustness of the result at multiple different combinations of the osmolyte mixtures and the polymer dispersion interactions. The second subsection provides a molecular-level mechanistic interpretation of the observed result using the theory of preferential solvation. The third subsection validates the underlying principle using a synthetic polymer namely polystyrene in a mixed solvent media. Finally, the conclusions and outlook section provides

5 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

an overview of our conclusions and possible biophysical implications.

Material and Methodology Simulation model: We have simulated two different systems in the current work: A model bead-on-a-spring polymer chain ( 32 beads) ( Figure 1a)) and a synthetic polymer namely 20-mer of polystyrene (Figure 1b)) in different aqueous media of osmolytes. Both of the systems have been previously studied by Mondal et al 30,31 extensively in the context of osmolyte-polymer interaction and the details of the models can be found else-where. 16,30,31 Here we provide a brief overview of the models.

a)

b)

Figure 1: A representative configuration of two systems of hydrophobic macromolecule studied in the manuscript: a) model polymer in extended configuration ( orange colored) and b) 20-mer of polystyrene in an aqueous mixture of TMAO (red colored) and urea (green colored).Water molecules are represented by blue dots. Briefly, the model polymer beads are uncharged and interact with the surrounding particles via dispersion interactions modelled using Lennard-Jones potential. The bead diameter σ is 0.4 nm. The hydrophobic character of the polymer is tuned via varying Lennard Jones interaction energy b . In this work, majority of the simulation has been performed by using b = 1.0 kJ/mol. However, we have also checked for the robustness of our result by using 6 ACS Paragon Plus Environment

Page 6 of 28

Page 7 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

an alternative b = 0.6 kJ/mol. The beads are connected by harmonic bond and harmonic angle, the values of which can be found in the previous works of Mondal et al 30 and Zangi et all 16 . Apart from the bond and the angle restraints, the 1-4 nonbonding interaction is also present among the beads. We have simulated water molecules via SPC/E model, 41 urea via Kirkwood-Buff potential, 42 and TMAO using the forcefield developed by shea and coworkers 32,43 which is considered an improvement over Kast forcefield. 44 We have nevertheless validated our result via repeating the simulation using two other forcefields of TMAO developed by Kast et al 44 and Netz and coworkers. 32,45 The geometric combining rules were used to model the intermediate interactions of both  and σ. We have simulated a total of ten different osmolyte compositions by inserting necessary number of TMAO, urea and water molecules in a box consisting of polymer: polymer solvated in a) neat water, b) 0.5 M aqueous TMAO solution, c) 1 M aqueous TMAO, d) 1M aqueous urea solution, e) 4.5 M aqueous urea solution f) 7M aqueous urea solution, g) aqueous mixture of 0.5 M TMAO + 4.5 M urea, h) aqueous mixture of 1 M TMAO + 1 M urea solution i) aqueous mixture of 1M TMAO + 4.5M urea, j) aqueous mixtures of 1M TMAO+7M urea and investigated the polymer conformational behavior in the respective mixtures. The approximate box dimension of the system, after energy minimization and NPT equilibration, ranged from 5.5 nm to 5.8 nm. In the second part of the work, We have also simulated a synthetic hydrophobic polymer namely a 20mer of polystyrene, previously simulated by Mondal et al. 31 To be consistent with previous work , we have used the identical parameter and concentration for polystyrene, TMAO, water and urea, as used previously. In this case, 4 different systems were simulated by inserting necessary number of TMAO, urea and water molecules : 20-mer of polystyrene solvated in a) neat water, b) 1 M aqueous TMAO, c) 7 M aqueous urea and d) aqueous mixture of 1 M TMAO + 7 M urea. The approximate box dimension of the system after energy minimization and NPT equilibration ranged from 6 to 6.5 nm. Simulation Method : The majority of the simulation efforts was invested in obtaining

7 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

free energy landscape of polymer conformation (both model Lennard Jones-polymer and polystyrene) along the radius of gyration of the polymer in each of the aforementioned osmolyte solutions. Towards this end, we have performed umbrella sampling simulations 46 using radius of gyration of the polymer as a collective variable. To generate the representative configurations required for each of the so-called ’umbrella sampling windows’, we first performed independent equilibrium Molecular Dynamics simulations for 100 ns at constant pressure of 1 bar and temperature of 300 K, starting with an extended all-trans configuration of the polymer. To avoid any bias and to supplement any missing configurations for subsequent umbrella sampling simulations, these simulations were also repeated starting with a collapsed configuration of the polymer. For umbrella sampling simulations, the values of the radius of gyration ranged from 0.4 nm to 1.2 nm at a spacing of 0.05 nm for the model Lennard-Jones polymer and 0.700 nm to 1.6 nm at a spacing of 0.05 nm for 20-mer of polystyrene. Each of the umbrella sampling simulations started with configuration corresponding to desired radius of gyration, as obtained from the aforementioned equilibrium simulations. Then each of the configurations were subjected to a harmonic restraint to ensure a gaussian distribution of the radius of gyration around each desired value of the radius of gyration. Each of the umbrella sampling windows were sampled for 20 ns using NPT ensemble. We used Nose-Hoover thermostat 47,48 for maintaining the average temperature of 300 K and the Parrinello-Rahman barostat 49 for maintaining the average pressure of 1 bar. All the water molecules were simulated as rigid molecule using SETTLE algorithm. 50 Finally, Weighted Histogram Analysis Method ( WHAM) 51,52 was used over the last 10 ns of each of the umbrella sampled trajectories to generate unbiased histograms and the corresponding potentials of mean force or free energies. The total simulation length for all umbrella sampling simulations for each osmolyte solution was 340 ns for the model Lennard-Jones polymer ( and hence a total of 3.4 microsecond simulation for 10 different systems) and 380 ns for polystyrene chain. All simulations were performed using Gromacs 5.0.6 53 patched with Plumed 54 plugin to enable umbrella sampling

8 ACS Paragon Plus Environment

Page 8 of 28

Page 9 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

simulation along the radius of gyration. We employed an experimentally relevant quantity called preferential binding coefficient (Γs ) 2,3   Nstot − ns · nw Γs = ns − tot Nw − nw

(1)

where ns is the number of cosolute ( urea or TMAO) called bound to the polymer and Nstot is the total number of cosolutes in the system. On the other hand, nw is the number of water molecules bound to the polymer and Nwtot is the total number of water molecules in the system. Γs quantifies the excess of cosolute molecules s in the polymer solvation shell as compared to its average concentration in the solution. We calculated the values of Kp and Γ as a function of distance from the polymer. Towards this end, we employed similar protocol as previously implemented by Mondal et al. 30,31 To compare the relative preferential binding of each of the cosolutes, we took the central atom (O for water, N for TMAO and C for Carbon) of the cosolutes, solvent and any polymer bead for which the distance is shortest. Towards this end, we carried out additional simulations which are as follows: To obtain better statistics, we froze the polymer either in the collapsed or in the extended conformation of the polymer and then simulated the systems for 15 ns and repeated for five independent velocity seeds (hence a total of 75 ns for each of the conformation). We have also verified the effect of the magnitude of force used to restrain the polymer in either collapsed or extended conformation by repeating the simulation in the presence of a weak harmonic force along the radius of gyration of polymer. As depicted in figure S1, we have found that the qualitative results hold good irrespective of the magnitude of restraints. The simulations were repeated for 15 ns in each of the collapsed and extended conformation of polymer for each osmolyte solutions. In each cases, Γs was averaged over these ensemble of trajectories, and the standard deviations were obtained via block averaging. The entire protocol was repeated for both the model 32-bead polymer and the 20-mer of polystyrene. We also have validated our result against another relevant quantity called local-bulk

9 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 28

partition coefficient Kp , which is defined as 7

Kp =

hns i Nwtot · . hnw i Nstot

(2)

Here hnX i is the average number of molecules of type X bound to polymer and NXtot is the total number of molecules of type X in the system (where X = s stands for the cosolvent (urea or TMAO) and X = w stands for water). Kp is intensive and reflects the affinity of the cosolvent for the polymer regardless of the exposed surface area of the polymer.

Results and Discussion Effect of osmolyte mixtures on a model polymer Even though we have used multiple combinations of osmolytes, we first mainly focus on the conformational behavior of the model polymer in an aqueous mixture of 1 M TMAO and 4.5 M urea. Both concentrations are widely used to study the effects of the respective osmolytes in vitro and in silico, majorly based on the fact that a very low concentration of TMAO is known to counteract the action of high concentration of urea in a protein. Nevertheless, we have also explored multiple other concentrations of the osmolytes which are described in Method section and the relevant results will also be discussed in the aforementioned paragraphs. To quantitatively explore the entire conformational degrees of freedom, we simulate the free energy landscape of the polymer conformations by performing a series of umbrella sampling simulations along the radius of gyration of the polymer, in individual solution namely neat water, an aqueous solution of 1M TMAO, aqueous solution of 4.5 M urea and finally a mixed osmolyte solution of 1M TMAO and 4.5 M urea. Figure 2a and 2b compares the free energy landscape and the corresponding probability distributions of the polymer conformations along the radius of gyration in each of the aforementioned solutions.

10 ACS Paragon Plus Environment

Page 11 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

b)

a)

εb = 1.0 kJ/mol

Figure 2: (a) Potentials of mean force and (b) corresponding probability distribution along radius of gyration for polymer chain with an b = 1.0 kJ/mol. The representative polymer configurations corresponding to collapsed ( Rg ∼ 0.45 nm), hairpin (Rg ∼ 0.60 nm) and extended ( Rg ∼ 1.0 nm) conformations which are the locations of major basins. As evident from figure 2a and 2b, in all aqueous media, the free energetically most favorable conformation of polymer corresponds to a collapsed conformation ( Rg ∼ 0.45 nm), which is representative of a hydrophobic polymer, along with a local basin corresponding to the extended conformation ( Rg ∼ 1.0 nm ) and a local intermediate representing a semicompact hairpin conformation ( Rg ∼ 0.6 nm). However, as gleaned from the computed free energy profiles in figure 2 a, the extent of relative free energy barrier to extended conformation is systematically influenced by different osmolyte solution. Relative to neat water, the extended state gets more destabilized in 1M TMAO solution and gets more stabilized in 4.5 M urea solution, thereby apparently behaving like a protecting and denaturing osmolyte respectively, as is general trend in case of a globular protein. However, in a mixture of these two osmolyte solution ( 1M TMAO + 4.5M urea), the polymer extended state gets further stabilized, even compared to that in pure urea solution. The relative stabilization of the extended state of polymer in a mixed environment of osmolyte, even compared to that in pure urea solution implies that, in a mixed osmolyte solution, TMAO actually promotes the

11 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

denaturation propensity of the polymer. This act of reinforcement of extended or unfolded state of a hydrophobic polymer by TMAO in a mixed osmolyte solution is distinctly different from well-known counteraction mechanism of TMAO and urea prevalent in proteins and thereby indicates a completely different mechanism that underlies a polymer in a mixed osmolyte solution. Our result is robust with respect to the change of osmolyte concentration in the mixture, different hydrophobicity of the polymer and variation of TMAO forcefield. Figure 3a delineates the concentration dependence of free energy profiles of polymer along radius of gyration, obtained using umbrella sampling techniques. As evident from the free energy profiles of the polymer in 0.5 M TMAO and 1 M TMAO in figure 3a, addition of TMAO molecules in neat water monotonically increases the free energy barrier to the extended conformation, thereby systematically destabilizing the polymer’s extended state relative to neat water. However, gradual addition of TMAO in an aqueous media of 4.5 M urea solution (thereby constituting mixed media of 0.5 M TMAO with 4.5 M urea and 1 M TMAO with 4.5 M urea) monotonically stabilizes the polymer’s extended state relative to that in pure urea solution. The overall result points to an interesting observation that addition of TMAO in neat water protects the hydrophobic polymer against unfolding while addition of TMAO in pure urea solution promotes the formation of extended state of the same polymer. The observed trend of further stabilization of the extended conformation in a mixed osmolyte media, relative to pure aqueous media holds qualitatively true with respect to the variation of urea concentration from 1M to 7M. As illustrated in figure S2 a) and b), both an aqueous mixture of 1M TMAO + 7M urea and an aqueous mixture of 1M TMAO + 1M urea stabilize the extended polymer conformation relative to individual aqueous osmolyte solution, thereby acting as a potentially stronger denaturing agent than the pure aqueous urea media. We note that we have not explored the TMAO concentration beyond 1M majorly because it has been recently observed 19 that pure aqueous TMAO itself acts as a denaturant at a higher concentration.

12 ACS Paragon Plus Environment

Page 12 of 28

Page 13 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

b)

a)

εb=0.6 kJ/mol

ε =1.0 kJ/mol b

d)

c)

ε b=1.0 kJ/mol

εb=1.0 kJ/mol

Figure 3: (a)Concentration dependence of Potentials of mean force along polymer ( = 1.0 kJ/mol) radius of gyration using Shea et al’s 32,43 forcefield of TMAO . (b) Comparison of Potentials of mean force along the radius of gyration with b = 0.6 kJ/mol. Validation of robustness of trend of Potentials of mean force along polymer radius of gyration c) with Kast forcefield 44 of TMAO and d) with Netz forcefield 32,45 of TMAO with polymer dispersion interaction b = 1.0 kJ/mol. We have also explored the effect of polymer dispersion interaction in the overall trend of the free energy profile. As depicted in figure 3b, decreasing the polymer dispersion interaction from  = 1 kJ/mol to 0.6 kJ/mol ( and thereby reducing the polymer-cosolvent interaction via combination rule) makes the extended conformation of the polymer freeenergetically most favorable configuration. However, the qualitative effect of mixed osmolyte media on the polymer free energy landscape remains the same, albeit weaker than that with  = 1.0 kJ/mol. The overall trend qualitatively holds true even when repeating the simulation using two independent TMAO forcefields developed by Kast et al 44 of TMAO, ( figure 3c) 13 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 28

and Netz and coworkers 32,45 (figure 3d).

Towards a preferential-solvation based molecular interpretation The aforementioned behavior of the hydrophobic polymer in different osmolyte solutions, although distinct from the precedent notion in protein, is in reasonable agreement with seminal Wyman-Tanford approach of preferential solvation. The Wyman-Tanford approach hypothesizes the preferential binding coefficient Γs (see equation 1) 2,3 as the key quantity to describe the dependence of the polymer conformational equilibria on cosolutes: Γs > 0 if the osmolyte preferentially accumulates next to the polymer and is < 0 if it is excluded from the polymer. Note that Γs is closely related to a frequently used parameter known as Kirkwood Buff integral. 32,42 The effect of preferential binding on a conformational equilibrium between the collapsed and the extended configurations C E (with an equilibrium constant K) is usually interpreted in terms of the thermodynamic analysis put forward by Wyman and Tanford, 2,3 which leads to ∂ ln K = ∆ΓC→E , ∂ ln as

(3)

where as is the activity of the cosolute in the binary solution. According to Eq. 3 and illustration of Figure 4a, an increase in the concentration of the cosolute would lead to the macromolecular unfolding if ∆Γ > 0, and in contrast would favor the folded state over the unfolded one if ∆Γ < 0. Figure 4b-d illustrates and compares the profiles of Γs as a function of distance from the polymer surface, for the collapsed and the extended state of polymer in different osmolyte solution. Irrespective of different osmolyte solution, we find Γs > 0 for both the collapsed and the extended conformation of the polymer at a relevant distance from the polymer close to first solvation shell of the polymer. In the context of the aforementioned approach of theory of preferential binding, Γs > 0 implies that both TMAO and urea preferentially binds to the polymer surface. However, the sign of ∆Γ i.e. the relative preferential binding of cosolutes in

14 ACS Paragon Plus Environment

Page 15 of 28

a)

b) ΔΓ < 0

k1 TMAO

k2 Collapsed (C)

Extended(E)

d)

ΔΓ > 0 (TMAO+Urea)

c) Urea

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

ΔΓ > 0 aq media of 1 M TMAO + 4.5 M urea

Figure 4: (a) A schematic representation Wyman-Tanford theory of preferential binding coefficient relating ∆Γ with conformational equilibrium constant. Comparison of Profiles of Preferential binding coefficient of b) TMAO in 1 M TMAO solution, (c) urea in 4.5 M urea solution. (d) TMAO+urea in mixed aqueous solution of 1 M TMAO and 4.5 M urea solution as a function of the distance from the polymer surface ( b = 1.0 kJ/mol) in a collapsed (red) or extended (black) configuration. Vertical dashed lines indicate the positions of the polymer first solvation shell in each case (5.5 ˚ A), and the black arrows in panels b, c and d represent the relative sign of ∆Γ. collapsed versus extended conformation, is distinctly different in different osmolyte solution, which determines the direction of conformational equilibrium. As depicted in figure 4b, in aqueous TMAO solution, ∆Γ < 0, implying relatively higher preferential binding of TMAO in collapsed conformation, which promotes the stabilization of collapsed state of polymer relative to neat water. On the other hand, as reflected in figure 4c and d in both aqueous urea solution and mixed TMAO-Urea solution, ∆Γ > 0, implying relatively lower preferential binding of cosolutes ( only urea in case of aqueous urea solution and all constituent cosolutes ( TMAO and urea) in a mixed solution) in collapsed conformation. Hence, in accordance with Wyman-Tanford approach, collapsed conformation of polymer is

15 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

destabilized in both pure urea and mixed TMAO-Urea solution, relative to neat water. In fact, higher value of ∆Γ in mixed osmolyte solution than that in pure urea solution suggests relatively higher extent of destabilization in mixed osmolyte solution than that in pure urea solution, which is consistent with our free-energy-based results. We note that, in the present work, in the case of mixed osmolyte solution, we have focussed on computing the preferential binding coefficient of the combination of TMAO + Urea. In other words, we have analyzed the cumulative preferential interaction of both TMAO and urea in order to discern their combined interaction on a polymer. We have never attempted to implement the WymanTanford relation to dissect the contribution arising due to only one of the two osmolytes in the mixture . Hence, in a mixture, we are using the Wyman-Tanford relation for effectively a two-component solvent mixtures , constituting of combination of TMAO + Urea as cosolute and water as solvent and hence Wyman-Tanford relation, in its current context, will hold fine. We also note that plausibility of extension of the Wyman-Tanford model to dissect the contribution of each of the cosolutes in a three-component mixture is worth exploring in a future work. We also quantified another relevant quantity local bulk partition coefficient (Kp ) which reflects the affinity of the cosolvent for the polymer regardless of the exposed surface area of the polymer. As shown in figure S3 in SI text, irrespective of osmolyte solution, Kp > 1, suggesting preferential binding of osmolytes on a polymer surface over water, complementing our finding based on Γ. More over, identical to the trend of Γ, relative to extended conformation, the local bulk partition coefficient of osmolyte in 1M TMAO solution is higher in collapsed conformation while lower in both 4.5M Urea and mixed solution of 1M TMAO and 4.5M Urea. Delving into the representative configuration of the collapsed and extended conformations of the polymer in its solvation shell provides us a better underlying picture of osmolyte’s action on polymer in different solution. As depicted in figure 5a, in aqueous TMAO solution, the first solvation shell of collapsed conformation has more TMAO molecules bound than

16 ACS Paragon Plus Environment

Page 16 of 28

Page 17 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

that in extended conformation, mainly because of prospect of more hydrophobic contacts between hydrophobic methyl groups of TMAO and collapsed state of polymer chain. a)

c)

b)

d) 1M TMAO

4.5M Urea

1M TMAO + 4.5M Urea

Figure 5: Representative snapshots of respective cosolutes bound first solvation shell of polymer in collapsed and extended conformation: a) TMAO molecules in 1M TMAO b) urea molecules in 4.5M urea, c) TMAO and urea molecules in mixture of 1M TMAO and 4.5M urea. the snapshots corresponding to collapsed configurations are kept larger so as to delineate the solvation shell clearly. d) Comparison of radial distribution function of TMAO and urea in different osmolyte solutions. On the other hand, in aqueous urea solution, as represented by figure 5b, more urea molecules accumulate next to extended conformation than collapsed conformation, mainly due to favorable access to higher surface area. However, in a mixed osmolyte solution, urea preferentially occupies the first solvation shell of the polymer, leaving very less number of sites for TMAO to occupy the polymer-surface. As a result, as shown in figure 5c, in the presence of urea molecules, in a mixture, TMAO molecules have very little rooms for making significant TMAO-polymer hydrophobic contacts in collapsed conformation, unlike that in pure aqueous TMAO solution. On the other hand, presence of higher surface area in extended conformation provides more sites to TMAO molecules in a mixed osmolyte environment. This is also quantified by polymer-cosolute radial distribution function as well ( figure 5d). Consequently, the cumulative number of cosolutes ( i.e. TMAO and urea 17 ACS Paragon Plus Environment

The Journal of Physical Chemistry

combined) in extended conformation is relatively higher than that in collapsed conformation in a mixed osmolyte solution.

How does a mixed osmolyte solution act towards a real polymer? Our aforementioned finding of effect of osmolyte solutions on model lennard Jones polymer is in reasonable agreement with simulation of a synthetic polymer namely polystyrene (degrees of polymerization= 20) in different osmolyte solution.

a)

ΔΓ < 0

polystyrene in osmolyte solution TMAO

b) polystyrene in aq. TMAO

polystyrene in mixture of TMAO+urea

polystyrene in aq. urea

(TMAO+Urea)

c) Urea

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 28

ΔΓ > 0

d) ΔΓ> 0

Figure 6: (a) Potentials of mean force along the radius of gyration for polystyrene in neat water, individual osmolyte solution and mixed osmolyte solution. Progression of preferential binding coefficient as a function of distance from polymer surface for b) TMAO in TMAO solution, c) urea in urea solution (d) (TMAO+urea) in mixed osmolyte solution. for b-d, both collapsed (red) or extended (black) configuration for polystyrene has been compared. Vertical dashed lines indicate the positions of the polymer first solvation shell in each case (5 ˚ A), and the black arrows in panels b, c and d represent the relevant sign of ∆Γ, thereby connecting the Wyman-Tanford theory with polymer conformational equilibrium.

18 ACS Paragon Plus Environment

Page 19 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

As depicted from the series of free energy profiles as a function of radius of gyration in figure 6a, the free energetically most stable conformation of the polystyrene chain corresponds to radius of gyration ∼ 0.750-0.800 nm, the so-called collapsed state of the polymer. However, the extent of stabilization of collapsed state of polystyrene relative to that of extended conformation is strongly dependent on the osmolyte composition of the solution. In commensurate with the aforementioned result of model polymer, the collapsed state of the polystyrene is relatively more stabilized in aqueous media of TMAO and less stabilized in aqueous solution of urea, relative to neat water. However, in a more striking agreement with the model polymer, addition of TMAO molecules to an aqueous solution of urea further destabilizes the collapsed state of the polystyrene, relative to aqueous media. Most importantly, the trend of preferential binding results in polystyrene are also in line with the Wyman-Tanford theory of preferential binding. As shown in figure 6b-d, the difference between the preferential binding coefficient of extended and collapsed configuration i.e. ∆Γ is negative in aqueous TMAO and positive in both aqueous urea and mixed solution of TMAO and urea, implying that, relative to neat water, collapsed conformation will be more stabilized in aqueous TMAO solution and more destabilized in aqueous urea and mixed solution of TMAO and urea solution.

Conclusions To conclude, our simulation results shed light on a unique role of mixed osmolyte solution on hydrophobic polymer, which is distinct from that in protein. The effect of osmolyte on macromolecule is usually investigated in protein and hence it is natural to compare the qualitative trend of our polymer-based result with the prevalent result of protein. In accordance with what is observed in protein, collapsed conformation of polymer is respectively stabilized and destabilized in aqueous solution of TMAO and aqueous solution of urea. However, the role of mixed solution of TMAO and urea is distinctly different in a hydrophobic polymer

19 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

compared to that of protein : while TMAO is known to counteract the effect of urea in a protein, our simulation result with a hydrophobic polymer suggests that TMAO, in a mixed environment, will reinforce the effect of urea by further destabilizing the collapsed conformation of polymer. The Wyman-Tanford theory of preferential binding coefficient sheds light on the key difference in the molecular mechanism of osmolyte-polymer interaction compared to that of osmolyte-protein interaction. For protein, preferential binding coefficient Γ is found to be < 0 for TMAO while > 0 for urea, implying that TMAO preferentially excludes from the surface of protein while urea prefers to bind to the surface of protein. However, for a hydrophobic polymer, the preferential binding coefficient Γ is > 0 for both TMAO and urea, which suggests that both TMAO and urea would bind to polymer surface. We believe that the preferential binding of TMAO to a hydrophobic polymer in contrast to preferential exclusion of the same from a protein surface has key implication in the end result of the polymer conformational behavior in a mixture of TMAO and urea, compared to that of protein. We find that urea preferentially coats the polymer surface in a mixed osmolyte solution, which provides very little prospect for TMAO to make hydrophobic contacts in collapsed state, relative to that in extended conformation, which can accommodate more TMAO molecules due to the larger surface area. Hence, cumulative number of cosolute molecules close to the polymer surface is higher in extended conformation of polymer than that in collapsed conformation, giving rise to positive ∆Γ and hence destabilization of collapsed conformation relative to aqueous media of urea solution. We hope that our current work on conformational behavior of hydrophobic polymer in a mixed osmolyte solution will stimulate novel experiments and theoretical efforts to dissect and understand the osmolyte-macromolecule interaction in the context of hydrophobic interaction. We also note that in a similar context, the mutually enhancing effect of TMAO and urea on denaturation` u has also recently been observed by Ganguly et al 39 where it is shown that TMAO at high concentration reduces hydrophobic aggregation and while a mixed osmolyte solution of TMAO and urea reduces the aggregation further. Moreover, recently polymer conformational behavior at

20 ACS Paragon Plus Environment

Page 20 of 28

Page 21 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

higher TMAO concentrations have been explored reporting intriguing observation. 19 Future direction will involve the effect of wide range of osmolyte concentrations 55 and other intrinsic properties of polymer namely charge and degree of polymerization and cononsolvency 56 on the hydrophobic interaction.

Acknowledgments This work was supported by computing resources obtained from shared facility of TIFR Center for Interdisciplinary Sciences, India. JM would like to acknowledge intramural research grants obtained from TIFR, India and Ramanujan Fellowship of DST, India. Part of the work was carried out in San diego supercomputing resources provided by XSEDE ( TG-CHE150024). JM would like to thank Dr. Isaac Li and Dr. Kalyaneswar Mandal for a critical reading of the manuscript and Dr. Guillaume Stirnemann for useful discussions.

Supporting Information Available Illustration of effect of restraining forces on preferential binding, free energy profiles at different combination of TMAO and urea concentrations, local-bulk partition coefficient(Kp ) for different osmolyte solution This material is available free of charge via the Internet at http://pubs.acs.org/.

References (1) Yancey, P. H.; Clark, M. E.; Hand, S. C.; Bowlus, R. D.; Somero, G. N. Living With Water-stress- Evolution of Osmolyte Systems. Science 1982, 217, 1214–1222. (2) Wyman, J. Linked Functions and Reciprocal Effects in Hemoglobin: A Second Look. Adv. Protein Chem. 1964, 19, 223–286.

21 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(3) Tanford, C. Extension of the Theory of Linked Functions to Incorporate the Effects of Protein Hydration. J. Mol. Biol. 1969, 39, 539–544. (4) Timasheff, S. N. Protein-solvent Preferential Interactions, Protein Hydration, and the Modulation of Biochemical Reactions By Solvent Components. Proc. Natl. Acad. Sci. USA 2002, 99, 9721–9726. (5) Parsegian, V. A.; Rand, R. P.; Rau, D. C. Osmotic Stress, Crowding, Preferential Hydration, and Binding: A Comparison of Perspectives. Proc. Natl. Acad. Sci. USA 2000, 97, 3987–3992. (6) Bolen, D.; Rose, G. Structure and Energetics of the Hydrogen-Bonded Backbone in Protein Folding. Annu. Rev. Biochem. 2008, 77, 339–362. (7) Courtenay, E. S.; Capp, M. W.; Anderson, C. F.; Record, M. T., Jr Vapor Pressure Osmometry Studies of Osmolyte-Protein Interactions: Implications For the Action Of Osmoprotectants in Vivo and For the Interpretation of ”Osmotic Stress” Experiments in Vitro. Biochemistry 2000, 39, 4455–4471. (8) Canchi, D. R.; Garc´ıa, A. E. Cosolvent Effects on Protein Stability. Annu. Rev. Phys. Chem. 2013, 64, 273–293. (9) Ufnal, M.; Zadlo, A.; Ostaszewski, R. TMAO: A Small Molecule of Great Expectations. Nutrition 2015, 31, 1317 – 1323. (10) Street, T. O.; Bolen, D. W.; Rose, G. D. A Molecular Mechanism for Osmolyte-Induced Protein Stability. Proc. Natl. Acad. Sci. USA 2006, 103, 13997–14002. (11) Wang, A. J.; Bolen, D. W. A Naturally Occurring Protective System in Urea-rich Cells: Mechanism of Osmolyte Protection of Proteins Against Urea Denaturation. Biochemistry 1997, 36, 9101–9108.

22 ACS Paragon Plus Environment

Page 22 of 28

Page 23 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(12) Hochachka, P. W.; Somero, G. N. Biochemical Adaptation. Mechanism and Process in Physiological Evolution; Oxford University Press, 2002. (13) Lin, T. Y.; Timasheff, S. N. Why Do Some Organisms Use a Urea-methylamine Mixture as Osmolyte- Thermodynamic Compensation of Urea and Trimethylamine N-oxide Interactions with protein. Biochemistry 1994, 33, 12695–12701. (14) Athawale, M. V.; Dordick, J. S.; Garde, S. Osmolyte trimethylamine N-oxide Does Not Affect the Strength of Hydrophobic Interactions: Origin of Osmolyte Compatibility. Biophys. J. 2005, 89, 858–866. (15) Paul, S.; Patey, G. N. The Influence of Urea and Trimethylamine N-oxide on Hydrophobic Interactions. J. Phys. Chem. B 2007, 111, 7932–7933. (16) Zangi, R.; Zhou, R.; Berne, B. J. Urea’s Action on Hydrophobic Interactions. J. Am. Chem. Soc. 2009, 131, 1535–1541. (17) Rodriguez-Ropero, F.; van der Vegt, N. F. A. Direct Osmolyte-Macromolecule Interactions Confer Entropic Stability to Folded States. J. Phys. Chem. B 2014, 118, 7327– 7334. (18) Reddy, P. M.; Taha, M.; Venkatesu, P.; Kumar, A.; Lee, M.-J. Destruction of Hydrogen Bonds of Poly(N-isopropylacrylamide) Aqueous Solution By Trimethylamine N-oxide. J. Chem. Phys. 2012, 136, 234904. (19) Rodriguez-Ropero, F.; Rotzscher, P.; van der Vegt, N. F. A. Comparison of Different TMAO Force Fields and Their Impact on the Folding Equilibrium of a Hydrophobic Polymer. J. Phys. Chem. B 2016, 120,8757–8767. (20) Tatzelt, J.; Prusiner, S.; Welch, W. Chemical Chaperones Interfere With The Formation of Scrapie Prion Protein. EMBO J. 1996, 15, 6363–6373.

23 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(21) Nandi, P. K.; Bera, A.; Sizaret, P.-Y. Osmolyte Trimethylamine N-oxide Converts Recombinant Alpha-helical Prion Protein to Its Soluble Beta-structured Form at High Temperature. J. Mol. Biol. 2006, 362, 810–820. (22) Rose, G. D.; Fleming, P. J.; Banavar, J. R.; Maritan, A. A backbone-based theory of protein folding. Proc. Natl. Acad. Sci. USA 2006, 103, 16623–16633. (23) Tseng, H. C.; Graves, D. J. Natural Methylamine Osmolytes, Trimethylamine N-oxide and Betaine, Increase Tau-induced Polymerization of Microtubules. Biochem. Biophys. Res. Comm. 1998, 250, 726–730. (24) Scaramozzino, F.; Peterson, D. W.; Farmer, P.; Gerig, J. T.; Graves, D. J.; Lew, J. TMAO Promotes Fibrillization and Microtubule Assembly Activity in the C-terminal Repeat Region of Tau. Biochemistry 2006, 45, 3684–3691. (25) Liu, Y.; Bolen, D. The Peptide Backbone Plays a Dominant Role in Protein Stabilization by Naturally-occuring Osmolytes. Biochemistry 1995, 34, 12884–12891. (26) Auton, M.; Roesgen, J.; Sinev, M.; Holthauzen, L. M. F.; Bolen, D. W. Osmolyte Effects on Protein Stability and Solubility: A Balancing Act Between Backbone and Side-chains. Biophys. Chem. 2011, 159, 90–99. (27) Auton, M.; Ferreon, A. C. M.; Bolen, D. W. Metrics that Differentiate the Origins of Osmolyte Effects on Protein Stability: A Test of the Surface Tension Proposal. J. Mol. Biol. 2006, 361, 983 – 992. (28) Pica, A.; Graziano, G. On Urea’s Ability To Stabilize the Globular State of Poly(Nisopropylacrylamide). Phys. Chem. Chem. Phys. 2016, 18, 14426–14433. (29) Zhang, Y.; Cremer, P. S. Chemistry of Hofmeister Anions and Osmolytes. Annu. Rev. Phys. Chem. 2010, 61, 63–83.

24 ACS Paragon Plus Environment

Page 24 of 28

Page 25 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(30) Mondal, J.; Stirnemann, G.; Berne, B. J. When Does Trimethylamine N-Oxide Fold a Polymer Chain and Urea Unfold It? J. Phys. Chem. B 2013, 117, 8723–8732. (31) Mondal, J.; Halverson, D.; Li, I. T.; Stirnemann, G.; Walker, G. C.; Berne, B. J. How Osmolytes Influence Hydrophobic Polymer Conformations: A Unified View From Experiment and Theory. Proc. Natl. Acad. Sci. USA 2015, 112, 9270–9275. (32) Ganguly, P.; Hajari, T.; Shea, J.-E.; van der Vegt, N. F. A. Mutual Exclusion of Urea and Trimethylamine N-Oxide from Amino Acids in Mixed Solvent Environment. J. Phys. Chem. Lett. 2015, 6, 581–585. (33) Bennion, B. J.; Daggett, V. Counteraction of Urea-induced Protein Denaturation By Trimethylamine N-oxide: A Chemical Chaperone at Atomic Resolution. Proc. Natl. Acad. Sci. USA 2004, 101, 6433–6438. (34) Sharp, K. A.; Madan, B.; Manas, E.; Vanderkooi, J. M. Water Structure Changes Induced By Hydrophobic And Polar Solutes Revealed By Simulations And Infrared Spectroscopy. J. Chem. Phys. 2001, 114, 1791–1796. (35) Hunger, J.; Tielrooij, K.-J.; Buchner, R.; Bonn, M.; Bakker, H. J. Complex Formation in Aqueous Trimethylamine-N-oxide (TMAO) Solutions. J. Phys. Chem. B 2012, 116, 4783–4795. (36) Sagle, L. B.; Zhang, Y.; Litosh, V. A.; Chen, X.; Cho, Y.; Cremer, P. S. Investigating the Hydrogen-Bonding Model of Urea Denaturation. J. Am. Chem. Soc. 2009, 131, 9304–9310, PMID: 19527028. (37) Cho, S. S.; Reddy, G.; Straub, J. E.; Thirumalai, D. Entropic Stabilization of Proteins by TMAO. J. Phys. Chem. B 2011, 115, 13401–13407. (38) Hu, C. Y.; Lynch, G. C.; Kokubo, H.; Pettitt, B. M. Trimethylamine N-oxide Influence On The Backbone of Proteins: An Oligoglycine Model. Proteins 2010, 78, 695–704. 25 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(39) Ganguly, P.; van der Vegt, N. F. A.; Shea, J.-E. Hydrophobic Association in Mixed ¨ ıTMAO Solutions. J. Phys. Chem. Lett. 2016, 7, 3052–3059. Urea A` (40) Canchi, D. R.; Jayasimha, P.; Rau, D. C.; Makhatadze, G. I.; Garcia, A. E. Molecular Mechanism for the Preferential Exclusion of TMAO from Protein Surfaces. J. Phys. Chem. B 2012, 116, 12095–12104. (41) Berendsen, H.; Grigera, J.; Straatsma, T. P. The Missing Term in Effective Pair Potential. J.Phys.Chem. 1987, 91, 6269-6271. (42) Kirkwood, J. G.; Buff, F. P. The Statistical Mechanical Theory of Solutions. I. J. Chem. Phys. 1951, 19, 774–777. (43) Larini, L.; Shea, J.-E. Double Resolution Model for Studying TMAO/Water Effective Interactions. J. Phys. Chem. B 2013, 117, 13268–13277. (44) Kast, K. M.; Brickman, J.; Kast, S. M.; Berry, R. S. Binary Phases of Aliphatic NOxides and Water: Force Field Development and Molecular Dynamics Simulation. J. Phys. Chem. A 2003, 107, 5342–5351. (45) Schneck, E.; Horinek, D.; Netz, R. R. Insight into the Molecular Mechanisms of Protein Stabilizing Osmolytes from Global Force-Field Variations. J. Phys. Chem. B 2013, 117, 8310–8321. (46) Torrie, G.; Valleau, J. Nonphysical Sampling Distributions In Monte Carlo Free-energy Estimation: Umbrella sampling. Journal of Computational Physics 1977, 23, 187 – 199. (47) Nos´e, S. A Molecular Dynamics Method For Simulations In The Canonical Ensemble. Mol. Phys. 1984, 52, 255. (48) Hoover, W. Canonical dynamics: Equilibrium Phase-space Distributions. Phys. Rev. A 1985, 31, 1695. 26 ACS Paragon Plus Environment

Page 26 of 28

Page 27 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(49) Parrinello, M.; Rahman, A. Polymorphic Transitions In Single Crystals: A New Molecular Dynamics Method. J. App. Phys. 1981, 52, 7182–7190. (50) Miyamoto, S.; Kollman, P. Settle: An Analytical Version of The SHAKE and RATTLE Algorithm For Rigid Water Models. J. Comput. Chem. 1992, 13, 952–962. (51) Kumar, S.; Bouzida, D.; Swendsen, R. H.; Kollman, P. A.; Rosenberg, J. M. The Weighted Histogram Analysis Method For Free-energy Calculations on Biomolecules. I. The method. J. Comput. Chem. 1992, 13, 1011–1021. (52) A.Grossfield, WHAM: The Weighted Histogram Analysis Method, version 2.0.9. http://membrane.urmc.rochester.edu/content/version-209 , 2014. (53) Hess, B.; Kutzner, C.; Van der Spoel, D.; Lindahl, E. GROMACS 4 : Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. J. Chem. Theory Comput. 2008, 4, 435–447. (54) Tribello, G. A.; Bonomi, M.; Branduardi, D.; Camilloni, C.; Bussi, G. PLUMED2: New Feathers For An Old Bird. Comput. Phys. Commun. 2014, 185, 604 – 613. (55) Rani, A.; Jayaraj, A.; Jayaram, B.; Pannuru, V. Trimethylamine-N-oxide Switches From Stabilizing Nature: A Mechanistic Outlook Through Experimental Techniques And Molecular Dynamics Simulation. Sci. Rep. 2016, 6, 23656. (56) Mochizuki, K.; Pattenaude, S. R.; Ben-Amotz, D. Influence of Cononsolvency on the Aggregation of Tertiary Butyl Alcohol in Methanol-Water Mixtures. J. Am. Chem. Soc. 2016, 138, 9045–9048, PMID: 27363494.

27 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Collapsed (C)

Extended(E)

28 ACS Paragon Plus Environment

Page 28 of 28