Rearrangement of the water network hydrating a protein can provide a mechanism for allostery, according to a study by Peter Hamm and colleagues (B. Buchli et al., PNAS 10.1073/pnas.1306323110; paper here). They insert an azobenzene photoswitch in the binding groove of a PDZ domain protein, a common system for studying allostery, such that photo-induced isomerisation induces a conformational change similar to that which occurs on ligand binding. This enables them to use fast IR spectroscopy to look at the opening of the binding groove in a precisely controlled fashion, having first characterized the equilibrium structures using NMR. Their MD simulations show that a change in water density in the vicinity of the photoswitch right after switching propagates slowly through the water network over about 100 ns until it reaches the back of the protein. They suggest that this change in hydration structure could then induce, either dynamically or structurally, a remote allosteric change in protein conformation. In such a case, this would be a particularly dramatic example of how the hydration network is really a part of the functional apparatus of the protein.
How aquaporins transport water across membranes – and specifically how they do so without also transporting protons – has been a topic of much debate. The consensus has come to focus on the so-called NPA motif in the centre of the channel, a bottleneck through which water kolecules pass in single file, which seems to prohibit proton transport via electrostatic repulsion. Urszula Kosinska Eriksson at the University of Gothenburg and colleagues have recently reported a new high-resolution crystal structure of yeast aquaporin 1 which sheds new light on the issue (Science 340, 1346; 2013 – paper here). As Jeff Abramson and Armand Vartanian of ULCA explain in an accompanying perspective (Science 340, 1294 – paper here), the structure shows that proton transport isn’t (as some have suggested) blocked by hydrogen-bonding of a single water molecule in the NPA region to two asparagines. Rather, there are two independent waters here, but the interactions with the asparagines constrain their dynamics in such a way as to effectively break the ‘water wire’ threading through the channel. The authors say that water transport then seems to happen in pairwise fashion, similar to the ion transport in potassium channels – a possible example of convergent evolution to solve related problems.
A carefully structured water cluster also seems to play an important role in oxygen evolution during the photocycle of photosystem II, as Brandon Polander and Bridgette Barry of Georgia Tech deduce using laser flashes to induce the process and following it by FTIR spectroscopy (PNAS 10.1073/pnas1306532110 – paper here). The hydrogen-bonded cluster of five water molecules, bound to the catalytic Mn4CaO5 cluster, seems to become protonated during the S1→S2 part of the cycle, and it stores the proton until a later stage of the reaction. Ammonia can poison the reaction by disrupting the water cluster.
Water molecules that gain access to the interiors of globular proteins can act as probes of the intrinsic conformational dynamics that enable proteins to function. These waters in the ‘dry’ protein interior have been studied by magnetic relaxation dispersion spectroscopy, but atomistic models are needed to interpret those results. The problem is that such deep water penetration tends to be a rare event, demanding very long (millisecond) run times for simulations. A technique has recently been developed that enables this (D. E. Shaw et al., Science 330, 341; 2010), and now Filip Persson and Bertil Halle at Lund have used the method to compare MD with the MRD experiments for bovine pancreatic trypsin inhibitor (JACS 135, 8735; 2013 – paper here). They find that some of these internal hydration sites have water residence times of several microseconds, and that the water molecules gain access along single-file hydrogen-bonded chains.
Dewetting transtions are known to be important for at least some instances of hydrophobic assembly, and Bruce Berne and colleagues at Columbia now extend their earlier studies of dewetting-induced protein collapse to look at the potential role of such transitions in the docking of hydrophobic ligands in their binding pockets (J. Mondal et al., preprint http://www.arxiv.org/abs/1305.7505). In this way the solvent dynamics, which are retarded in the concave cavity, are explicitly included in the kinetics of the binding process – the process can be parametrized through a state variable that describes whether the pocket is ‘wet’ or ‘dry’, while the ligand diffuses across a potential-energy surface that can switch between these two states.
Daniel Sindhikara and Fumio Hirata of the Ritsumeikan University in Japan present a fast algorithm, based on the three-dimensional reference interaction site model (3D-RISM), for calculating the solvent distribution around solutes (J. Phys. Chem. B 117, 6718; 2013 – paper here). This is a wholly theoretical approach derived from the Ornstein-Zernicke equation. They say that it gives water positions and orientations that agree well with available experimental data, and demonstrate its use on HIV-1 protease.
Paul Ben Ishai at the Hebrew University of Jerusalem and colleagues have used quasielastic neutron scattering to look at salt effects on water dynamics (J. Phys. Chem. B 117, 7724; 2013 – paper here). They find that water diffusion is slower than in pure water, on average, in NaCl solution, but faster in KCl. They interpret the result in terms of structure-making and –breaking, saying that the disruption of the hydrogen-bonding network by potassium ions accounts for its apparent ‘lubricating’ effect.
The hydration state of arginine side chains can be deduced from its UV resonance-enhanced Raman spectrum, according to Sanford Asher at Pittsburgh and colleagues (Z. Hong et al., J. Phys. Chem. B 117, 7145; 2013 – paper here). Their density-functional calculations show that a particular vibration of this residue is sensitive to hydration. They use this signal to characterize differing degrees of hydration of Arg in two polyAla model peptides.
Monday, July 15, 2013
Friday, June 28, 2013
Complexity of osmolytes
More on the mechanisms of denaturation and osmolyte stabilization. Trimethylamine-N-oxide (TMAO) is known to confer protection against protein denaturation by urea. Rahul Sarma and Sandip Paul at the Indian Institute of Technology have carried out MD simulations to try to figure out why (J. Phys. Chem. B 117, 5691; 2013 – paper here). They say that TMAO interacts with a model peptide, N-methylacetamide (NMA), so as to cause some dehydration by replacing solvation water. However, this interaction with the peptide is relatively inefficient, especially because TMAO cannot donate its hydrogen to the backbone carbonyls. As a result, TMAO is not efficient at stabilizing the unfolded peptide chain. If the protein folds, then TMAO is more available for forming very strong hydrogen bonds with water, and indeed with urea. It’s a subtle argument in which the direct interactions between all four components – TMAO, urea, peptide and water – are implicated.
Much the same issues are explored in a new preprint by Bruce Berne and coworkers at Columbia (J. Mondal et al., http://www.arxiv.org/abs/1306.4642). They look at the interactions of urea and TMAO separately with a model 32-mer made of Lennard-Jones beads whose hydrophobicity can be varied. Although both small solutes interact strongly with the polymer, TMAO destabilizes extended conformations while urea stabilizes them. This seems to be because TMAO molecules are more stable (via van der Waals interactions) next to the collapsed rather than the extended polymer. The results provide further support for an interpretation of osmolyte effects that invokes direct interactions rather than indirect effects on “water structure”.
The subtlety of mixed-osmolyte solutions is revealed in simulations of a hydrophobic polymer in urea and guanidinium chloride by Payel Das at IBM Yorktown Heights and colleagues (P. Das et al., Langmuir 29, 4877; 2013 – paper here). They find that, while both of these molecules act as denaturants on their own, in combination they actually promote collapse of the polymer. It is hard to tell a simple story about why this happens, although apparently Paul Flory predicted back in 1955 that two good solvents could combine to induce polymer chain collapse. Although Gdm has the stronger interaction with the polymer, its enhanced concentration in the vicinity of the polymer attracts urea due to the favourable urea-Gdm interaction, with the result that urea is in fact preferentially adsorbed onto the polymer. This sets up a long-ranged interaction between the monomers mediated by their clouds of urea molecules, ultimately driving collapse.
The effect of progressive dehydration of a protein (lysozyme) on its structure and dynamics is studied using Raman spectroscopy by Gediminas Niaura of Vilnius University of colleagues (J. Phys. Chem. B 117, 4981; 2013 – paper here). They find that there is a structural change in the dry protein crystal that begins at about 7-10 wt% water, and that the native state is reached at about 35 wt% water (which amounts to appreciably more than monolayer coverage). In the dry state the protein is dynamically glassy, with a diminished proportion of alpha helices and an enhanced content of beta-sheet contacts.
And on the question of ‘dry’ proteins, Eric Gloaguen at the CNRS Laboratoire Francis Perrin in Gif-sur-Yvette and colleagues report that small peptides with aromatic residues will fold into hydrophobic domains even in the solvent-free gas phase, showing that this is a favourable conformation even in the earliest stages of protein folding (E. Gloaguen et al., J. Phys. Chem. B 117, 4945; 2013 – paper here).
Somedatta Pal and Sanjoy Bandyopadhyay at the Indian Institute of Technology in Kharagpur have used MD simulations to look at the complementary issue of how protein (here barstar) dynamics affect the dynamics of hydration water (J. Phys. Chem. B 117, 5848; 2013 – paper here). They compare the normal situation of the hydrated protein with that in which the protein is kept frozen, looking specifically at the effect of the change on the low-frequency vibrations of the water molecules. Freezing the protein results in stronger confinement of water bound to the protein surface, with a corresponding blue shift of the vibrational frequency for transverse water oscillations – but much less effect on longitudinal oscillations.
An intriguing paper by Nicholas Spencer and colleagues at ETH investigates hydration forces for glycoproteins using the surface force apparatus (R. M. Espinosa-Marzel et al., Biophys. J. 104, 2686; 2013 – paper here). They attach these onto hydrophobic and hydrophobic surfaces, and find that there is a rather long-ranged repulsive force (several tens of nanometres) between the surfaces. It is good to see this neglected class of biological molecules get some attention in terms of their fundamental hydration characteristics. However, I don’t buy the interpretation that the repulsion is due to some “long-ranged structuring of water”. The whole discussion around the results seems to be something of a throwback, starting with the whole notion of “vicinal water” and involving two-state water theories, kosmotropes and all the rest of the paraphernalia that tended to surround discussions of ‘water structure’ 20 years ago. “Ordering” of water over 30 nm or so just isn’t any longer consistent with what is known of other systems, and it’s been suggested to me that steric repulsion of the surfaces due to a few of the macromolecules protruding from the monolayer is a much more likely explanation of the effects observed.
But could it be that long-ranged water structuring is going to try to stage a little comeback? I ask this because it is invoked in another recent paper too, by Jan Christer Eriksson and Ulf Henriksson of the Royal Institute of Technology in Sweden (Langmuir 29, 4789; 2013 – paper here). They develop their earlier argument (Langmuir 23, 1126; 2007) that the long-ranged hydrophobic attraction might be accounted for by the formation of roughly cylindrical bridging water clusters that are “slightly more organized than the rest of the film”. Their analysis suggests that such a proposal can explain some recent measurements on water thin films with the SFA (Wang et al., J. Coll. Int. Sci. 364, 257; 2011). But I think I will stick with the bridging-nanobubbles idea.
Much the same issues are explored in a new preprint by Bruce Berne and coworkers at Columbia (J. Mondal et al., http://www.arxiv.org/abs/1306.4642). They look at the interactions of urea and TMAO separately with a model 32-mer made of Lennard-Jones beads whose hydrophobicity can be varied. Although both small solutes interact strongly with the polymer, TMAO destabilizes extended conformations while urea stabilizes them. This seems to be because TMAO molecules are more stable (via van der Waals interactions) next to the collapsed rather than the extended polymer. The results provide further support for an interpretation of osmolyte effects that invokes direct interactions rather than indirect effects on “water structure”.
The subtlety of mixed-osmolyte solutions is revealed in simulations of a hydrophobic polymer in urea and guanidinium chloride by Payel Das at IBM Yorktown Heights and colleagues (P. Das et al., Langmuir 29, 4877; 2013 – paper here). They find that, while both of these molecules act as denaturants on their own, in combination they actually promote collapse of the polymer. It is hard to tell a simple story about why this happens, although apparently Paul Flory predicted back in 1955 that two good solvents could combine to induce polymer chain collapse. Although Gdm has the stronger interaction with the polymer, its enhanced concentration in the vicinity of the polymer attracts urea due to the favourable urea-Gdm interaction, with the result that urea is in fact preferentially adsorbed onto the polymer. This sets up a long-ranged interaction between the monomers mediated by their clouds of urea molecules, ultimately driving collapse.
The effect of progressive dehydration of a protein (lysozyme) on its structure and dynamics is studied using Raman spectroscopy by Gediminas Niaura of Vilnius University of colleagues (J. Phys. Chem. B 117, 4981; 2013 – paper here). They find that there is a structural change in the dry protein crystal that begins at about 7-10 wt% water, and that the native state is reached at about 35 wt% water (which amounts to appreciably more than monolayer coverage). In the dry state the protein is dynamically glassy, with a diminished proportion of alpha helices and an enhanced content of beta-sheet contacts.
And on the question of ‘dry’ proteins, Eric Gloaguen at the CNRS Laboratoire Francis Perrin in Gif-sur-Yvette and colleagues report that small peptides with aromatic residues will fold into hydrophobic domains even in the solvent-free gas phase, showing that this is a favourable conformation even in the earliest stages of protein folding (E. Gloaguen et al., J. Phys. Chem. B 117, 4945; 2013 – paper here).
Somedatta Pal and Sanjoy Bandyopadhyay at the Indian Institute of Technology in Kharagpur have used MD simulations to look at the complementary issue of how protein (here barstar) dynamics affect the dynamics of hydration water (J. Phys. Chem. B 117, 5848; 2013 – paper here). They compare the normal situation of the hydrated protein with that in which the protein is kept frozen, looking specifically at the effect of the change on the low-frequency vibrations of the water molecules. Freezing the protein results in stronger confinement of water bound to the protein surface, with a corresponding blue shift of the vibrational frequency for transverse water oscillations – but much less effect on longitudinal oscillations.
An intriguing paper by Nicholas Spencer and colleagues at ETH investigates hydration forces for glycoproteins using the surface force apparatus (R. M. Espinosa-Marzel et al., Biophys. J. 104, 2686; 2013 – paper here). They attach these onto hydrophobic and hydrophobic surfaces, and find that there is a rather long-ranged repulsive force (several tens of nanometres) between the surfaces. It is good to see this neglected class of biological molecules get some attention in terms of their fundamental hydration characteristics. However, I don’t buy the interpretation that the repulsion is due to some “long-ranged structuring of water”. The whole discussion around the results seems to be something of a throwback, starting with the whole notion of “vicinal water” and involving two-state water theories, kosmotropes and all the rest of the paraphernalia that tended to surround discussions of ‘water structure’ 20 years ago. “Ordering” of water over 30 nm or so just isn’t any longer consistent with what is known of other systems, and it’s been suggested to me that steric repulsion of the surfaces due to a few of the macromolecules protruding from the monolayer is a much more likely explanation of the effects observed.
But could it be that long-ranged water structuring is going to try to stage a little comeback? I ask this because it is invoked in another recent paper too, by Jan Christer Eriksson and Ulf Henriksson of the Royal Institute of Technology in Sweden (Langmuir 29, 4789; 2013 – paper here). They develop their earlier argument (Langmuir 23, 1126; 2007) that the long-ranged hydrophobic attraction might be accounted for by the formation of roughly cylindrical bridging water clusters that are “slightly more organized than the rest of the film”. Their analysis suggests that such a proposal can explain some recent measurements on water thin films with the SFA (Wang et al., J. Coll. Int. Sci. 364, 257; 2011). But I think I will stick with the bridging-nanobubbles idea.
Friday, May 10, 2013
No ice please
I’m not a fan of the Frank-Evans iceberg picture of hydrophobic hydration. But I will happily admit that, on the question of whether water hydrating small nonpolar solutes is more/less ordered and mobile than that in the bulk, the evidence is mixed. Nuno Galamba of the University of Lisbon now offers some more support for a weak version of the iceberg picture (J. Phys. Chem. B 117, 2153; 2013 – paper here). His MD simulations show that a subset of water molecules in the first hydration shell of methane have “significantly enhanced tetrahedrality and a slightly larger number of hydrogen bonds”, as well as slower reorientational dynamics, relative to the bulk. He adds that these characteristics should not be visible in the rdfs deduced from neutron scattering, explaining why they have not been seen experimentally. Whether this view extends to large hydrophobes is another matter, but these results at least argue for some small degree of water ‘ordering’, even if this is very far from ice-like.
The role of hydration water in the dynamics of hydrophobic sidechains of peptides is explored by Daniela Russo at the ILL in Grenoble and colleagues, using inelastic neutron scattering and simulations (J. Phys. Chem. B 117, 2829; 2013 – paper here). They say that the activation energy for methyl group motions increases with increasing level of hydration but eventually reaches a plateau when an extended hydrogen-bonding network is established around the group, which happens when it is essentially surrounded by a single layer of water. These sidechain dynamics seem to have a critical impact on the flexibility of the peptide as a whole, and so the degree of hydration appears to determine the onset of conformational freedom for the entire polypeptide.
Neeraj Sinha and colleagues at the Centre of Biomedical Magnetic Resonance in Lucknow have taken on the challenge of exploring water-protein interactions in a rather complex system, the helical coat protein of the filamentous virus Pf1 (R. N. Purusottam et al., J. Phys. Chem. B 117, 2837; 2013 – paper here). Using 2D proton-N15 NMR, they deduce that the filamentous assembly has a highly hydrated core which not only acts as a ‘glue’ but might also mediate the interaction of the arg44 residue with DNA.
Membrane transporters are membrane protein channels that pump directional transport of small molecules across the membrane, coupled to chemical energy sources in the cell. Their operation has generally been considered to involve a carefully orchestrated sequence of conformational changes to ensure one-way and selective transport. But it has been found that sometimes water and ions get through too, and the question arises of whether this is passive, osmotically driven ‘leakage’ or stoichiometric co-transport of these species. Emad Tajkhorshid and colleages at Illinois at Urbana investigate this question using MD simulations, for several classes of membrane transporters (J. Li et al., PNAS 110, 7696; 2013 – paper here). They find that these generally support states in which there are channels that permit passive water flux – but that this does not interfere with the coordinated vectorial transport of the primary substrate. Thus the transporters are imperfect, but not problematically so: as the authors put it, “Given the soft mechanical properties of transporter proteins, it comes as no surprise to observe harmless imperfections in the overall gating motions, which manifest themselves in the formation of water-conducting states. It would, of course, be a concern if these channels were large enough to leak the substrate, and/or long-lived to allow very large amounts of smaller species to permeate across the membrane. Neither of these aspects appears to be the observed in our results, because the leaky states are only large enough for small species such as water. Furthermore, it appears that these states only transiently rise during the transport cycle, an attribute that might make them difficult to capture experimentally.”
Solvation of small molecules in water has often been described using point charges affixed to the solutes to represent polarization effects. David Cerutti at Rutgers and colleagues present a quantum chemical approach for fitting such partial charges to solutes like amino acids (D. S. Cerutti et al., J. Phys. Chem. B 117, 2328; 2013 – paper here). They say that their approach, which represents an evolution from the simplest point-charge models of previous decades, predicts substantially more polarization of amino acids than earlier efforts using AMBER force fields.
More on urea-induced protein denaturation, this time from Michela Candotti of the Institute of Research in Biomedicine in Barcelona and colleagues, who use MD, SAXS and NMR data to develop an atomistic picture of the unfolded states and the energetics of unfolding and refolding (M. Candotti et al., PNAS 110, 5933; 2013 – paper here). They conclude that urea’s denaturing influence is a combination of kinetic (disrupting stabilizing intramolecular contacts) and thermodynamic (stabilizing the extended conformation) effects. As I understand it, the results support models based on direct interactions of urea and protein, while also revealing that the urea-unfolded state is rather different to the denatured conformation that exists in pure water.
Hydration is evidently important to DNA conformation, influencing the A-B transition, interactions with DNA-binding proteins, and perhaps even affecting shape in a sequence-dependent manner. There is some evidence that the hydration water of DNA has collective dynamics of a glassy nature, and this notion is offered further support in inelastic neutron scattering experiments by Alessandro Paciaroni at the Università degli Studi of Perugia and colleagues (A. Paciaroni et al., J. Phys. Chem. B 117, 2026; 2013 – paper here). They find that at 100 K the large-wavevector scattering from water, related to coherent excitations, seems to imply a character related to amorphous ice. In other words, the interactions with DNA significantly alter the structure and dynamics of the interfacial water.
Melittin, one of the multi-subunit peptides that seems to aggregate by dewetting, is a component of been venom that acts as an antimicrobial. The melittin tetramer forms a pore that inserts into lipid membranes, and Max Berkowitz and colleagues at UNC now suggest that its effect is to create transient water-permeable channels, making the membranes leaky (K. P. Santo et al., J. Phys. Chem. B 117, 5031; 2013 – paper here). The technique they use for MD simulations can handle timescales of up to microseconds, and they show that the melittin peptides gradually aggregate into a kind of wedge that punctures the membranes and allows water to pass through, before falling apart again.
[FeFe] hydrogenase catalyses hydrogen-ion reduction to gaseous hydrogen, and could therefore be a useful biocatalyst. Simulations by Martin McCullagh and Greg Voth at Chicago now show that it seems to work by coupling electron transfer to proton transfer along a previously unknown water channel accessing the active site (J. Phys. Chem. B 117, 4062; 2013 – paper here).
There are seemingly strange rumours that fully miscible liquids are inhomogeneous on length scales of hundreds of nanometres, a conclusion suggested by some light-scattering and small-angle neutron-scattering studies. For example, Marián Sedlák of the Slovak Academy of Sciences reported such a claim in a series of papers in 2006 (e.g. J. Phys. Chem. B 110, 4329). Could these be solute clusters, or perhaps nanobubbles stabilized by adsorbed solute? Sedlák and Dmytro Rak now investigate that latter possibility (J. Phys. Chem. B 117, 2495; 2013 – paper here). They look at a range of solutes: magnesium sulphate, citric acid, urea, and t-butyl alcohol, and find no significant differences in the scattering from normal and degassed solutions, apparently ruling out the nanobubble interpretation.
The role of hydration water in the dynamics of hydrophobic sidechains of peptides is explored by Daniela Russo at the ILL in Grenoble and colleagues, using inelastic neutron scattering and simulations (J. Phys. Chem. B 117, 2829; 2013 – paper here). They say that the activation energy for methyl group motions increases with increasing level of hydration but eventually reaches a plateau when an extended hydrogen-bonding network is established around the group, which happens when it is essentially surrounded by a single layer of water. These sidechain dynamics seem to have a critical impact on the flexibility of the peptide as a whole, and so the degree of hydration appears to determine the onset of conformational freedom for the entire polypeptide.
Neeraj Sinha and colleagues at the Centre of Biomedical Magnetic Resonance in Lucknow have taken on the challenge of exploring water-protein interactions in a rather complex system, the helical coat protein of the filamentous virus Pf1 (R. N. Purusottam et al., J. Phys. Chem. B 117, 2837; 2013 – paper here). Using 2D proton-N15 NMR, they deduce that the filamentous assembly has a highly hydrated core which not only acts as a ‘glue’ but might also mediate the interaction of the arg44 residue with DNA.
Membrane transporters are membrane protein channels that pump directional transport of small molecules across the membrane, coupled to chemical energy sources in the cell. Their operation has generally been considered to involve a carefully orchestrated sequence of conformational changes to ensure one-way and selective transport. But it has been found that sometimes water and ions get through too, and the question arises of whether this is passive, osmotically driven ‘leakage’ or stoichiometric co-transport of these species. Emad Tajkhorshid and colleages at Illinois at Urbana investigate this question using MD simulations, for several classes of membrane transporters (J. Li et al., PNAS 110, 7696; 2013 – paper here). They find that these generally support states in which there are channels that permit passive water flux – but that this does not interfere with the coordinated vectorial transport of the primary substrate. Thus the transporters are imperfect, but not problematically so: as the authors put it, “Given the soft mechanical properties of transporter proteins, it comes as no surprise to observe harmless imperfections in the overall gating motions, which manifest themselves in the formation of water-conducting states. It would, of course, be a concern if these channels were large enough to leak the substrate, and/or long-lived to allow very large amounts of smaller species to permeate across the membrane. Neither of these aspects appears to be the observed in our results, because the leaky states are only large enough for small species such as water. Furthermore, it appears that these states only transiently rise during the transport cycle, an attribute that might make them difficult to capture experimentally.”
Solvation of small molecules in water has often been described using point charges affixed to the solutes to represent polarization effects. David Cerutti at Rutgers and colleagues present a quantum chemical approach for fitting such partial charges to solutes like amino acids (D. S. Cerutti et al., J. Phys. Chem. B 117, 2328; 2013 – paper here). They say that their approach, which represents an evolution from the simplest point-charge models of previous decades, predicts substantially more polarization of amino acids than earlier efforts using AMBER force fields.
More on urea-induced protein denaturation, this time from Michela Candotti of the Institute of Research in Biomedicine in Barcelona and colleagues, who use MD, SAXS and NMR data to develop an atomistic picture of the unfolded states and the energetics of unfolding and refolding (M. Candotti et al., PNAS 110, 5933; 2013 – paper here). They conclude that urea’s denaturing influence is a combination of kinetic (disrupting stabilizing intramolecular contacts) and thermodynamic (stabilizing the extended conformation) effects. As I understand it, the results support models based on direct interactions of urea and protein, while also revealing that the urea-unfolded state is rather different to the denatured conformation that exists in pure water.
Hydration is evidently important to DNA conformation, influencing the A-B transition, interactions with DNA-binding proteins, and perhaps even affecting shape in a sequence-dependent manner. There is some evidence that the hydration water of DNA has collective dynamics of a glassy nature, and this notion is offered further support in inelastic neutron scattering experiments by Alessandro Paciaroni at the Università degli Studi of Perugia and colleagues (A. Paciaroni et al., J. Phys. Chem. B 117, 2026; 2013 – paper here). They find that at 100 K the large-wavevector scattering from water, related to coherent excitations, seems to imply a character related to amorphous ice. In other words, the interactions with DNA significantly alter the structure and dynamics of the interfacial water.
Melittin, one of the multi-subunit peptides that seems to aggregate by dewetting, is a component of been venom that acts as an antimicrobial. The melittin tetramer forms a pore that inserts into lipid membranes, and Max Berkowitz and colleagues at UNC now suggest that its effect is to create transient water-permeable channels, making the membranes leaky (K. P. Santo et al., J. Phys. Chem. B 117, 5031; 2013 – paper here). The technique they use for MD simulations can handle timescales of up to microseconds, and they show that the melittin peptides gradually aggregate into a kind of wedge that punctures the membranes and allows water to pass through, before falling apart again.
[FeFe] hydrogenase catalyses hydrogen-ion reduction to gaseous hydrogen, and could therefore be a useful biocatalyst. Simulations by Martin McCullagh and Greg Voth at Chicago now show that it seems to work by coupling electron transfer to proton transfer along a previously unknown water channel accessing the active site (J. Phys. Chem. B 117, 4062; 2013 – paper here).
There are seemingly strange rumours that fully miscible liquids are inhomogeneous on length scales of hundreds of nanometres, a conclusion suggested by some light-scattering and small-angle neutron-scattering studies. For example, Marián Sedlák of the Slovak Academy of Sciences reported such a claim in a series of papers in 2006 (e.g. J. Phys. Chem. B 110, 4329). Could these be solute clusters, or perhaps nanobubbles stabilized by adsorbed solute? Sedlák and Dmytro Rak now investigate that latter possibility (J. Phys. Chem. B 117, 2495; 2013 – paper here). They look at a range of solutes: magnesium sulphate, citric acid, urea, and t-butyl alcohol, and find no significant differences in the scattering from normal and degassed solutions, apparently ruling out the nanobubble interpretation.
Thursday, April 25, 2013
Hydration and thermostability
What makes thermophilic proteins stable? And in particular, does hydration play a role? Those questions are examined by Fabio Sterpone of the Université Paris Diderot and colleagues (O. Rahaman et al., Phys. Chem. Chem. Phys. 15, 3570; 2013 – paper here). They study two homologous proteins, one from E. coli and the other from the thermophile Sulfolobus solfataricus, using MD simulations. The average water dynamics at the protein surface is the same in both cases, despite their different amino acid compositions, being slowed by a factor of 3-5 relative to the bulk. The authors conclude that this slowdown is primarily a geometric effect due to excluded volume, which explains why the protein sequence has little influence. This doesn’t exactly answer the initial question about thermostability – but it does suggest that there is nothing special about thermophilic proteins in terms of their hydration.
More on dewetting as a mechanism for protein folding and stability. Ruhong Zhou, Robert Matthews and their colleagues at Columbia have simulated a TIM barrel protein and found drying inside clusters of hydrophobic residues, signified by strong fluctuations in water density (P. Das et al., JACS 135, 1882; 2013 – paper here). In particular, a hydrophobic cluster (ILV) at the N-terminus shows drying that is weakened or suppressed by substituting some of the hydrophobic amino acids for less hydrophobic ones. A cluster at the C-terminus, meanwhile, seems to experience only partial drying that is unaffected by such substitutions. Experiments on the structure and stability of wild-type and mutant versions seem to back up these conclusions. The authors note that ILV clusters are common in several other protein motifs too.
Huib Bakker and colleagues at the FOM Institute in Amsterdam have used Förster resonant energy transfer between water OH stretches to find out where the water hydrating lipid membranes actually is (L. Piatkowski et al., J. Phys. Chem. B 117, 1367; 2013 – paper here). Somewhat surprisingly, they find that, even at rather low hydration levels, the water molecules are not evenly distributed among the lipid head groups, but form nanoscale clusters with an average intermolecular distance of 3.4 Å.
Our picture of the water surface continues to evolve. Is the water here structurally and dynamically different from that in the bulk? That idea is supported by experiments by James Skinner and colleagues at Wisconsin in simulations using a new three-body water potential (Y. Ni et al., PNAS 110, 1992; 2013 – paper here). They calculate that the hydrogen-bond switching dynamics are retarded at the water surface by a factor of about 3, although the rotational dynamics are actually a little faster. They say that vibrational 2D sum-frequency generation spectroscopy should be a good experimental method for investigating these dynamics, and calculate what the spectra should look like.
To make optimal use of plant cellulose as a feedstock for biofuels such as ethanol, it’s necessary to get it into aqueous solution. But cellulose is hard to solvate, which is why it is important to understand its hydration structure. Sylvia McLean at Oxford and colleagues have used neutron diffraction to look at the aqueous solvation of the disaccharide cellobiose, and find that there is (as has been suggested) a hydrogen bond across the glycosidic bond linking the two sugars (W. B. O’Dell et al., PLoS ONE 7, e45311; 2012 – paper here). There is competition from water molecules for the oxygen acceptor in this bond, however, with average occupancy of 50% for both water and the intramolecular OH donor.
More evidence that a hydrogen-bonded cluster of water molecules at the catalytic site plays a crucial role on photosynthesis: Bridgette Barry and colleagues at Georgia Tech use EPR to look at the rate of reaction of an intermediate neutral radical in the proton-coupled electron transfer that leads to oxygen evolution at the reaction centre of photosystem II, which contains a catalytic Mn4CaO5 cluster (J. M. Keough et al., J. Phys. Chem. B 117, 1296; 2013 – paper here). They find that this network is rearranged during the transition between the S0 and S2 states of the catalytic cycle, and that ammonia slows oxygen evolution because it disrupts the network by displacing water.
Suzi Jarvis and colleagues at University College Dublin have been using scanning probe microscopy for some time to probe hydration forces, and in their latest contribution they use a technique called frequency modulation AFM to look at hydration forces at the interface of mica and an electrolyte (J. I. Kilpatrick et al., JACS 135, 2628; 2013 – paper here). By ‘hydration force” here they mean the monotonically decaying force with a characteristic length of a few Å, without regard to its precise origin. They find that, relative to pure water, ions introduce or accentuate oscillations in the force as a function of distance due to the formation of distinct hydration layers. They point out that there are implications for obtaining atomic-resolution AFM images in aqueous saline solution.
Uzi Landman, Gary Schuster and colleagues at Georgia Tech offer a fascinating insight into the role of hydration water around DNA in the oxidation and consequent mutation of A/T-rich regions, which my be significant in the early stages of carcinogenesis due to stalling of replication (R. N. Barnett et al., JACS 135, 3904; 2013 – paper here). Their experiments and simulations indicate that oxidation of adenine leads to proton-coupled electron transfer to thymine, mediated by a water wire. This process can explain why nearly all the mutation due to such an oxidation event happens at thymine.
It’s sobering to realise that even now the reasons for the anomalous behaviour of water’s thermodynamic response functions, such as the divergence in the heat capacity and compressibility at low temperature, are still unknown. Francesco Mallamarce, Carmelo Corsaro and Gene Stanley discuss this issue with reference to experimental data on the power spectrum of sound velocity (F. Mallamarce et al., PNAS 110, 4899; 2013 – paper here). They conclude that the anomalies are associated with a structural transformation due to the appearance of an extended hydrogen-bonding network, which gives rise to viscoelastic behaviour in the liquid.
Also on bulk water, what happens to water’s structure above its critical point has been studied using X-ray Raman spectroscopy by Christoph Sahle of the Technical University of Dortmund and colleagues (C. Sahle et al., PNAS 110, 6301; 2013 – paper here). They find that, as one might expect, distortions of the hydrogen bonds are significant above the critical point, and the average coordination of each molecule decreases to just 0.6 at 600 oC and 134 MPa.
More on dewetting as a mechanism for protein folding and stability. Ruhong Zhou, Robert Matthews and their colleagues at Columbia have simulated a TIM barrel protein and found drying inside clusters of hydrophobic residues, signified by strong fluctuations in water density (P. Das et al., JACS 135, 1882; 2013 – paper here). In particular, a hydrophobic cluster (ILV) at the N-terminus shows drying that is weakened or suppressed by substituting some of the hydrophobic amino acids for less hydrophobic ones. A cluster at the C-terminus, meanwhile, seems to experience only partial drying that is unaffected by such substitutions. Experiments on the structure and stability of wild-type and mutant versions seem to back up these conclusions. The authors note that ILV clusters are common in several other protein motifs too.
Huib Bakker and colleagues at the FOM Institute in Amsterdam have used Förster resonant energy transfer between water OH stretches to find out where the water hydrating lipid membranes actually is (L. Piatkowski et al., J. Phys. Chem. B 117, 1367; 2013 – paper here). Somewhat surprisingly, they find that, even at rather low hydration levels, the water molecules are not evenly distributed among the lipid head groups, but form nanoscale clusters with an average intermolecular distance of 3.4 Å.
Our picture of the water surface continues to evolve. Is the water here structurally and dynamically different from that in the bulk? That idea is supported by experiments by James Skinner and colleagues at Wisconsin in simulations using a new three-body water potential (Y. Ni et al., PNAS 110, 1992; 2013 – paper here). They calculate that the hydrogen-bond switching dynamics are retarded at the water surface by a factor of about 3, although the rotational dynamics are actually a little faster. They say that vibrational 2D sum-frequency generation spectroscopy should be a good experimental method for investigating these dynamics, and calculate what the spectra should look like.
To make optimal use of plant cellulose as a feedstock for biofuels such as ethanol, it’s necessary to get it into aqueous solution. But cellulose is hard to solvate, which is why it is important to understand its hydration structure. Sylvia McLean at Oxford and colleagues have used neutron diffraction to look at the aqueous solvation of the disaccharide cellobiose, and find that there is (as has been suggested) a hydrogen bond across the glycosidic bond linking the two sugars (W. B. O’Dell et al., PLoS ONE 7, e45311; 2012 – paper here). There is competition from water molecules for the oxygen acceptor in this bond, however, with average occupancy of 50% for both water and the intramolecular OH donor.
More evidence that a hydrogen-bonded cluster of water molecules at the catalytic site plays a crucial role on photosynthesis: Bridgette Barry and colleagues at Georgia Tech use EPR to look at the rate of reaction of an intermediate neutral radical in the proton-coupled electron transfer that leads to oxygen evolution at the reaction centre of photosystem II, which contains a catalytic Mn4CaO5 cluster (J. M. Keough et al., J. Phys. Chem. B 117, 1296; 2013 – paper here). They find that this network is rearranged during the transition between the S0 and S2 states of the catalytic cycle, and that ammonia slows oxygen evolution because it disrupts the network by displacing water.
Suzi Jarvis and colleagues at University College Dublin have been using scanning probe microscopy for some time to probe hydration forces, and in their latest contribution they use a technique called frequency modulation AFM to look at hydration forces at the interface of mica and an electrolyte (J. I. Kilpatrick et al., JACS 135, 2628; 2013 – paper here). By ‘hydration force” here they mean the monotonically decaying force with a characteristic length of a few Å, without regard to its precise origin. They find that, relative to pure water, ions introduce or accentuate oscillations in the force as a function of distance due to the formation of distinct hydration layers. They point out that there are implications for obtaining atomic-resolution AFM images in aqueous saline solution.
Uzi Landman, Gary Schuster and colleagues at Georgia Tech offer a fascinating insight into the role of hydration water around DNA in the oxidation and consequent mutation of A/T-rich regions, which my be significant in the early stages of carcinogenesis due to stalling of replication (R. N. Barnett et al., JACS 135, 3904; 2013 – paper here). Their experiments and simulations indicate that oxidation of adenine leads to proton-coupled electron transfer to thymine, mediated by a water wire. This process can explain why nearly all the mutation due to such an oxidation event happens at thymine.
It’s sobering to realise that even now the reasons for the anomalous behaviour of water’s thermodynamic response functions, such as the divergence in the heat capacity and compressibility at low temperature, are still unknown. Francesco Mallamarce, Carmelo Corsaro and Gene Stanley discuss this issue with reference to experimental data on the power spectrum of sound velocity (F. Mallamarce et al., PNAS 110, 4899; 2013 – paper here). They conclude that the anomalies are associated with a structural transformation due to the appearance of an extended hydrogen-bonding network, which gives rise to viscoelastic behaviour in the liquid.
Also on bulk water, what happens to water’s structure above its critical point has been studied using X-ray Raman spectroscopy by Christoph Sahle of the Technical University of Dortmund and colleagues (C. Sahle et al., PNAS 110, 6301; 2013 – paper here). They find that, as one might expect, distortions of the hydrogen bonds are significant above the critical point, and the average coordination of each molecule decreases to just 0.6 at 600 oC and 134 MPa.
Friday, April 12, 2013
How antifreeze proteins work
A relatively short catch-up, with much more to come…
Martina Havenith at Bochum and her colleagues have extended their previous work on the mechanism of antifreeze proteins. Previously they looked at an antifreeze glycoprotein and found, using terahertz spectroscopy, that ice-binding by the protein seems to involve a long-range retardation of water H-bond dynamics, extending up to 2 nm from the molecular surface (S. Ebbinghaus et al., JACS 132, 12210; 2010). Now they find a similar effect operating for the antifreeze protein (AFP) of the fire-coloured beetle D. canadensis, a hyperactive insect (K. Meister et al., PNAS 110, 1617; 2013 – paper here). This contrasts, however, with the mechanism of another class of AFP, called wfAFP-1, which seems to operate only by short-ranged water binding to surface OH groups (S. Ebbinghaus et al., Biophys. J. 103, L20; 2012). As the authors say, “Nature is probably more inventive than initially thought and makes use of short- and long-range water perturbation to varying degrees in different classes of AFPs”.
Meanwhile, Ido Braslavsky of Ohio University and colleagues have used microfluidic methods to study the effects of AFPs on ice nucleation, and find that the growth of ice crystals is inhibited by the irreversible surface binding of the proteins (Y. Celik et al., PNAS 110, 1309; 2013 – paper here). These results help to rule out suggestions that direct binding of the AFPs to ice is not necessary to their mode of action.
More cases of water assisting receptor-substrate recognition and catalytic activity. First, Stephen Neidle at UCL and colleagues find that a cluster of 11 water molecules in an AT region of the minor groove of DNA seems to support the binding of three different small-molecule ligands (D.-G. Wei et al., JACS 135, 1369; 2013 – paper here). This cluster stabilizes the ligand by hydrogen bonding, and is also linked to (but distinct from) the well-known spine of hydration in B-DNA. Slight differences in binding mode with this cluster seem to account for the differences in binding affinity of the ligands: in other words, it is the water network that is calling the shots.
Second, Xiaoqing Wang and Hajime Hirao at Nanyang Technological University in Singapore say that the catalysis of myo-inositol monophosphatase (IMPase, a potential target for lithium treatment of bipolar disorder) is dependent on two bound water molecules (J. Phys. Chem. B 117, 833; 2012 – paper here). One provides the hydroxide ion that attacks the bound substrate. The other, coordinated to a magnesium ion, facilitates proton transfer leading to the product.
There seems to be an increasing perception that understanding the collective vibrations of proteins – their softness or rigidity – could offer insights into their enzymatic activity. Sow-Hsin Chen at MIT and coworkers support that view with an X-ray scattering study of the collective modes of hydrated lysozyme (Z. Wang et al., J. Phys. Chem. B 117, 1186; 2013 – paper here). They find that at low hydration levels both the collective ‘soft’ phonon modes and the enzymatic activity are much weaker or absent, suggesting a causal relationship: an indicator of the now familiar plasticizing effect of hydration.
C. Cametti at “La Sapienza” University of Rome and colleagues consider another aspect of protein hydration: how high concentrations of the protein (again lysozyme) can lead to clustering and a consequent decrease in average hydration number (C. Cametti et al., J. Phys. Chem. B 117, 104; 2013 – paper here). Their measurements of the dielectric spectra from 500 MHz to 50 GHz, which probe orientational relaxation, are consistent with this hypothesis of clustering into small aggregates at concentrations above about 100 mg/mL, which was first proposed by Stradner et al. (Nature 432, 492; 2004).
There is now some debate about whether the proposed two metastable liquid phases and associated critical point claimed on the basis of some simulations is real or not. That has been disputed by David Limmer and David Chandler at Berkeley, who extend their previous negative results using several different water potentials in a new preprint. I have written a commentary on the issue here.
Martina Havenith at Bochum and her colleagues have extended their previous work on the mechanism of antifreeze proteins. Previously they looked at an antifreeze glycoprotein and found, using terahertz spectroscopy, that ice-binding by the protein seems to involve a long-range retardation of water H-bond dynamics, extending up to 2 nm from the molecular surface (S. Ebbinghaus et al., JACS 132, 12210; 2010). Now they find a similar effect operating for the antifreeze protein (AFP) of the fire-coloured beetle D. canadensis, a hyperactive insect (K. Meister et al., PNAS 110, 1617; 2013 – paper here). This contrasts, however, with the mechanism of another class of AFP, called wfAFP-1, which seems to operate only by short-ranged water binding to surface OH groups (S. Ebbinghaus et al., Biophys. J. 103, L20; 2012). As the authors say, “Nature is probably more inventive than initially thought and makes use of short- and long-range water perturbation to varying degrees in different classes of AFPs”.
Meanwhile, Ido Braslavsky of Ohio University and colleagues have used microfluidic methods to study the effects of AFPs on ice nucleation, and find that the growth of ice crystals is inhibited by the irreversible surface binding of the proteins (Y. Celik et al., PNAS 110, 1309; 2013 – paper here). These results help to rule out suggestions that direct binding of the AFPs to ice is not necessary to their mode of action.
More cases of water assisting receptor-substrate recognition and catalytic activity. First, Stephen Neidle at UCL and colleagues find that a cluster of 11 water molecules in an AT region of the minor groove of DNA seems to support the binding of three different small-molecule ligands (D.-G. Wei et al., JACS 135, 1369; 2013 – paper here). This cluster stabilizes the ligand by hydrogen bonding, and is also linked to (but distinct from) the well-known spine of hydration in B-DNA. Slight differences in binding mode with this cluster seem to account for the differences in binding affinity of the ligands: in other words, it is the water network that is calling the shots.
Second, Xiaoqing Wang and Hajime Hirao at Nanyang Technological University in Singapore say that the catalysis of myo-inositol monophosphatase (IMPase, a potential target for lithium treatment of bipolar disorder) is dependent on two bound water molecules (J. Phys. Chem. B 117, 833; 2012 – paper here). One provides the hydroxide ion that attacks the bound substrate. The other, coordinated to a magnesium ion, facilitates proton transfer leading to the product.
There seems to be an increasing perception that understanding the collective vibrations of proteins – their softness or rigidity – could offer insights into their enzymatic activity. Sow-Hsin Chen at MIT and coworkers support that view with an X-ray scattering study of the collective modes of hydrated lysozyme (Z. Wang et al., J. Phys. Chem. B 117, 1186; 2013 – paper here). They find that at low hydration levels both the collective ‘soft’ phonon modes and the enzymatic activity are much weaker or absent, suggesting a causal relationship: an indicator of the now familiar plasticizing effect of hydration.
C. Cametti at “La Sapienza” University of Rome and colleagues consider another aspect of protein hydration: how high concentrations of the protein (again lysozyme) can lead to clustering and a consequent decrease in average hydration number (C. Cametti et al., J. Phys. Chem. B 117, 104; 2013 – paper here). Their measurements of the dielectric spectra from 500 MHz to 50 GHz, which probe orientational relaxation, are consistent with this hypothesis of clustering into small aggregates at concentrations above about 100 mg/mL, which was first proposed by Stradner et al. (Nature 432, 492; 2004).
There is now some debate about whether the proposed two metastable liquid phases and associated critical point claimed on the basis of some simulations is real or not. That has been disputed by David Limmer and David Chandler at Berkeley, who extend their previous negative results using several different water potentials in a new preprint. I have written a commentary on the issue here.
Friday, February 22, 2013
On icebergs
Still we ponder the nature of hydrophobic hydration: is it more ‘ice-like’ in some sense? Apparently so, say Dor Ben-Amotz and colleagues at Purdue (J. G. Davis et al., Nature 491, 582; 2012 – paper here). They use Raman spectroscopy to monitor vibrational spectra of water hydrating linear alcohols ranging from methanol to heptanol, and see more tetrahedral ordering of water molecules, and fewer weak hydrogen bonds, at low temperatures. But for hydrophobic chains longer than about 1 nm this hydration structure gives way to one in which the water is more disordered and has weaker H-bonds at higher temperatures. This scale-dependent crossover is reminiscent of that proposed by Lum, Chandler & Weeks above about 1 nm (J. Phys. Chem. B 103, 4570; 1999).
The nature of the hydration shells of methanol, ethanol and propanol have been probed using THz spectroscopy by Vladimir Matvejev of the Free University of Brussels and colleagues (J. Phys. Chem. B 116, 14071; 2012 – paper here). They estimate that the shells comprise about 14, 23, 23 and 31 molecules for methanol, ethanol, 1- and 2-propanol, respectively, and the water molecules are retarded by a factor of around 1.4.
In a somewhat similar vein, L. Luca of the University of Perugia and colleagues use GHz-THz light scattering to probe the hydration of mono- and disaccharides (J. Phys. Chem. B 116, 14760; 2012 – paper here). They find that slowing of water collective reorientation (by relatively large factors of 5-6) occurs only over relatively short distances (3-4 Å or essentially the first hydration layer), regardless of the size of the sugar molecules. This retardation involves considerably more water molecules than those few instantaneously hydrogen-bonded to the sugar.
A more generic approach to small-molecule solvation is described by Alla Oleinikova and Ivan Brovchenko of the Dortmund University of Technology (J. Phys. Chem. B jp306781y – paper here). They use MC simulations to study water structure around spherical solute particles 3-10 Å in size that vary from strongly hydrophobic to strongly hydrophilic. In all cases there is a density depletion relative to the bulk due simply to the missing-neighbour effect. For strongly hydrophobic particles there is a drying transition at the surface. Similar effects are seen for other fluids, but the directional hydrogen bonding of water enhances them.
Water in useful places: there is likely to be water molecules forging a hydrogen-bonded assembly in the mammalian photoreceptor melanopsin in the retina, which triggers the biological clock, according to Sivakumar Sekharan and colleagues at Cornell (JACS 134, 19536; 2012 – paper here). Their first-principles calculations of the active site suggest that two water molecules bridge the Schiff base and residues on the peptide, accounting for the blue shift of the optical absorption relative to the closely related rhodopsin.
Cytochrome c might control water access to its heme centre to tune the reduction potential via hydration changes, according to quantum MD simulations of Isabella Daidone of the University of L’Aquila and colleagues (C. A. Bortolotti et al., JACS 134, 13670; 2012 – paper here). They find that a yeast cytochrome has two channels that it can open to admit water, altering the reduction potential and thus refining processes of electron transfer.
How water evaporates from the air-water interface has been studied in detail by Patrick Varilly and David Chandler at Berkeley (J. Phys. Chem. B jp310070y – paper here). They find that the escape trajectory of a water molecule can be described in terms of two parameters: the distance from the instantaneous interface and the velocity along the surface normal. The results seem to imply that evaporation has zero activation energy, as some but not all experiments have suggested.
Fresh variety in the structures of ice confined to nanopores is reported by Jaeil Bai and Xiao Cheng Zeng of the University of Nebraska (PNAS 109, 21240; 2012 – paper here). Their simulations indicate that a bilayer of ice-like water molecules in slit-like pores about 7-9Å apart can be transformed under pressure to an amorphous phase and then to a very-high-density amorphous phase at 250K and 3 GPa. For rapid compression to 6 GPa an entirely new VHD ordered phase is found in which the water molecules are linked into square nanotubes.
Andreas Barth and colleagues at Stockholm University propose that studying changes in the water absorption bands in infrared spectroscopy can offer a way of quickly and remotely monitoring for binding of a ligand to a target protein in drug development (S. Kumar et al., J. Phys. Chem. B 116, 13968; 2012 – paper here). Specifically, they look at the bands diagnostic of ‘bound’ water in the binding cavity, and see this decline as ligand binding expels the water. They have so far looked at rather high protein concentrations, but think that the approach should work at lower concentrations with brighter light sources.
Hydrophobic drugs can be rendered more water-soluble by adding certain solutes known as hydrotropes. Although this has been long known, the mechanism is not clear. Seishi Shimizu of the University of York and colleagues have studied the question using thermodynamic theory (the so-called fluctuation theory of solution) and measurements of thermodynamic quantities (J. J. Booth et al., J. Phys. Chem. B 116, 14915; 2012 – paper here).
Junrong Zheng at Rice University and coworkers report evidence for strong segregation of ions and water in strong electrolyte solutions of thiocyanate ions (J. Phys. Chem. B jp310153n – paper here). The clustering is enhanced by the addition of strongly hydrated ions such as fluoride, whereas iodide ions tend to associate with the SCN clusters. In any event, the aqueous systems are strongly inhomogeneous.
Finally for now – though lots more still to catch up on – I have a brief article on the Chemistry World site (here) on a recent paper by Thomas Kühne and Rustam Khaliullin at Mainz (Nat. Commun. 4, 1450; 2013 – paper here) that seems to shed light on the arguments over XAS studies of liquid water by Anders Nilsson and colleagues (Wernet et al., Science 304, 995; 2004). Views on whether there was a controversy here in the first place will doubtless vary, but it seems that this new work provides a useful perspective on what the XAS work was showing. Anders gave me some rather extensive remarks on his view of the matter, which I thought might be usefully reported here. Forgive me for presenting them here unmediated – they will hopefully at least clarify Anders’ current position on this. Everything that follows is from him:
I find this paper to represent a very important step forward. It is very close to our original suggestion in the Wernet et al. Science paper that many molecules will be in instantaneous configurations with one strong and one weak hydrogen bond, i.e. asymmetric configurations. XAS measures the electronic structure on a timescale of a few attoseconds which means that the molecules have no time to move so we are only detecting snapshots of space-averaged frozen configurations (also stated in the original Wernet et al. paper).
However, there is a difference. In real water we expect the asymmetry between the strong and weak bonds to be much larger than based on the structures from the current simulation. It is mentioned in the methods section but not directly in the main text. Here is the underlying experimental evidence that requires no spectroscopy interpretation.
I attach a recent study (JCP in press, Skinner et al.) where we together with Benmore’s group have undertaken a major effort to determine a more accurate O-O pair distribution function (PDF) of ambient water. This is based on 4 new x-ray diffraction data sets (only 3 are consistent at high Q), all with a much higher Q cut-off than any previous measurements where we have also taken serious care to remove any OH contributions. Figure 9a shows the O-O PDF with very small error bars with a first peak height of 2.57. If you comparethis with the O-O PDF presented in the supplementary material it is clear that the simulations have an overstructured main peak. For the PBE functional, used in the main text, the peak height is 3.3. The TPSS –D3-FF simulation shown in the supplementary information has a peak height of 3.1, still overstructured but less so than the PBE. The consequence in the asymmetry is clearly visible in terms of the weak bond energy distribution from fig. 3a and S3a where the latter has more contributions towards lower energies. If we have to further dramatically understructure the liquid down to a peak height of 2.57 in the O-O PDF we expect the asymmetry to become much larger. In my opinion their asymmetry is only a lower limit.
Most likely the underlying reason for the asymmetry is in the many-body cooperativity effects that only an electronic structure simulation can capture as demonstrated in the current paper, and which is not represented by classical force fields. It has been known that the cooperativity effect is strongest when you have one strong donor and one strong acceptor bond. Thereby for 2 hydrogen bonded structures the energy per hydrogen bond is higher than for 4 bonded or tetrahedral coordinated structures. We discussed this in a few sentences in the Wernet et al. paper with references that I have underlined with yellow in the attachment. I also attach a recent review by Lars, Congcong and myself on water (A. Nilsson et al., J. Mol. Liq. 176, 2; 2012) where we discuss the importance to further develop simulations and in particular the importance of many-body effects in an electronic structure description where it is also essential that the latter includes van der Waals interactions.
Another interesting aspect in their study is the dynamics indicating that the strong and weak bonds switch place on ultrafast time scales. This is also in line with our previous discussions. I underline in the Leetma et al. paper (J. Chem. Phys. 129, 084502; 2008) how we discussed ultrafast measurements where the strongly H-bonded OH group in one water molecule could switch places with the weak one. The difference was only that we assumed that the switching occurred via librations whereas in their simulations it is mostly via translational motions of other molecules.
We had a last sentence in our Wernet et al. paper of a more speculative nature that with mostly only one strong donor and one strong acceptor bond per molecule, one-dimensional structures should appear which could be chains or rings. I noticed that they have a sentence regarding molecular chains at each instant connecting the strong bonds. Very nice.
Another point not related to the current paper that could also be of interest for you is the recent paper by Overduin and Patey (J. Phys. Chem. B 116, 12014; 2012) that discusses our PNAS Huang et al. paper based on simulations coming to similar conclusions regarding inhomogeneities in the liquid but only using a different language in terms of concentration fluctuations of two different classes of structures.
I think Kühne and Khaliullin put the picture in terms of a symmetry breaking in order to be not too far away from the most accepted tetrahedral picture of water. But if you really look at what they claim, which is also aligned with us, is that most molecules will be in an asymmetric position at all instances when taking a snapshot. It is simply that an OH group that is either weak or strong switches on a rather fast timescale. This would mean that the symmetric position is never really seen. It is only when you average over a long time that it looks like it is mostly in the tetrahedral position. Like me make analog. Take a pendulum that swings back and forth. The speed at the end points are close to zero and at the middle position it is highest. It means that the pendulum spends most time at the end points and extremely little time at the middle position. If you take the average position it will be the middle position but it is hardly ever visited. You could in such a picture claim that the equilibrium position is in the middle if you average over a single period but the pendulum will spend very little time there. The question is how the surrounding will be affected by such a motion. If I understand correctly the paper it is the surrounding that indeed infer the asymmetry. They claim that it is the translational motion of the other molecules that provides the mechanism of the switching. This should mean in my opinion that it is the end points of the asymmetric motion that makes the interaction with the surrounding and not the time averaged position. If it would have been dominated by a more internal motion such as librations it could have had less effect on the surrounding.
In water my own belief is that it is somewhat more complicated. Let me come to my current picture of water at the end of these comments.
I think it is currently not possible to observe the asymmetry based on the current status of pump-probe IR spectroscopy. Please don't quote me on this but I believe that we are still missing some major understanding about IR spectroscopy. The vibrational life time as a local oscillator in H20 is too short. Nearly all experiments on ultrafast dynamics are based on isotope substitution with HDO impurities in either H2O and D2O and measurements of the time resolved development of the OH or OD frequency. It is then assumed that HDO has equal probabilities for all molecular positions in the liquid. This assumption has never been proven. Based on our understanding of water in terms of fluctuations of two local structures, tetrahedral (low density water) and asymmetric (high density water) the latter has more contributions to the free energy through stronger bond energy whereas the latter has more entropy. I attached a slide from one of my presentations so you can see what I mean (hit a return in slide show mode to see the motion). There you also see the switching of the two bonds. Since the asymmetric configuration will provide more entropy and the HDO is already asymmetric it is likely that HDO will reside more in asymmetric configurations. Furthermore, we can anticipate through quantum effects that also the OD and OH hydrogen bonds are different. I would assume that OD will more likely be situated with the strong bond and OH more in the weak bond. This was observed in a recent PRL by Misha Bonn's group (attached) where they showed based on sum frequency generation spectroscopy and simulations that there is a strong preference for OH to point to the vacuum and OD to the liquid side for HDO at the water-air interface. At this point I think caution needs to be exercised regarding pump probe IR since these measurements might only be probing some fraction of all possible molecular motions. We (Lars) are currently further investigating the HDO dynamics based on quantum simulations. There are also recent indications that the pump pulse disturbs the dynamics and these measurements might not always represent the equilibrium motions.
We are currently using the new x-ray laser at Stanford, the Linac Coherent Light Source (LCLS) to various problems related to water. The machine provides completely new opportunities. We are planning a future step to open up dynamics. In particular x-ray correlation spectroscopy with x-ray lasers could provide new answers. This is probably 2 years away. Since the x-ray laser is fully coherent there will be a speckle pattern due to diffraction also from a disordered material. The plan is to make probe-probe measurement to thereby study equilibrium dynamics. You split the x-ray beam into two pulses with a controllable delay and the change in the speckle pattern will give information about how molecules have moved between the two pulses. These will be most challenging experiments but hopefully can be done and would open for completely new avenue's to probe dynamics in liquids. Naturally one challenge will be to not allow the first pulse to disturb the system but water is quite forgiving for hard x-rays since it is a low Z liquid.
We have demonstrated in a number of papers that this asymmetry can't been seen in diffraction since the data is not completely perfect. The difficulty is to have a technique that is only sensitive to the hydrogen bonding asymmetry around individual molecules. You can have a total hydrogen bonding average through a linear combination of water molecules with nearly no hydrogen bonds and fully tetrahedral water that on average would be like an asymmetric species. It is the rehybridization of the molecular orbital structure leading to local OH orbitals instead of delocalized H2O orbitals that makes XAS sensitive. I attach a paper where we discuss this (Nilsson section 3.5). It has to do with O2p and O2s hybridization and XAS only provides intensity for the O2p part in the orbitals. It is quite involved. Maybe it is time for a simpler review so many can understand these effects.
Here is our current picture of water. At high temperatures (close to boiling and above) water behaves as a simple liquid where most interactions are isotropic (dominated by van der Waals interactions). This is a structure where many molecules are in very disordered shells and without having a well defined first and second shell (this structure you get from ab initio MD simulations with van der Waals functionals even at ambient temperature, see discussion in my previous mail in paper "fluctuations in ambient water"). This type of non directional orientations gives high flexibility for various motions and high entropy. As the liquid cools down the water molecules starts to stick to each other through directional hydrogen bonds. This will appear in two classes of configurations, tetrahedral and asymmetric. In the tetrahedral structures each molecules are in 4 hydrogen bonds, this provides the lowest enthalpic energy. Since cooperatively effects makes the bonds stronger if water is bonded to other waters that are also in tetrahedral structures, these start to clump together in small local regions. Since the molecules are stuck with four bonds the motion is very restricted (see the attached slide) and thereby low entropy.
The other alternative is to form asymmetric structures. Here the hydrogen bond energy per molecule is higher than in the tetrahedral structures but since it is fewer bonds the enthalpic energy is lower. With less directional hydrogen bonds you have more flexibility for motions and thereby higher entropy (see slides). In this structure there are also non isotropic molecular interactions causing interstitial molecules. These are therefore called by us and others, preferable in the supercooled community, as high density water. Another way of viewing such species is that we start not with hexagonal ice but with high density amorphous or very high density amorphous ices where we have a large collapse of the second shell. In these ices we have also a local tetrahedral bonding but with other angles towards the second shell (often call interstitials). Here we induce asymmetric distortions around these tetrahedral bonded molecules in the first shell but simultaneous keeping the interstitials ( see fluctuations in ambient water). We have also seen this in water at an interface recently published in the Nature journal Scientific Reports (see attached Kaya et al.). There are fluctuations between the tetrahedral structures and asymmetric or high density structures. As we cool the liquid down the molecules in the asymmetric structures converts more and more into the tetrahedral structures which grows in size since the enthalpic energy contributes more to the free energy with decreasing temperature. The timescales in the fluctuations between these two classes should also slow down (not yet determined). There is also a continuous change in the asymmetric structures with temperature. The switching time is expected to slow down and thereby also the amplitude in the motion with decreasing temperature. We will be approaching more and more a local tetrahedral arrangement. However, not as in hexagonal ice but more towards high density ices with still many interstitials. That the density is decreasing below 4C is simply due to that we are at the same time converting many molecules into tetrahedral structures which have a more similar structure as hexagonal ice but disordered more towards low density amorphous ice (low density water) with the second shell at the tetrahedral angle causing an open network with low density. In my opinion we have to consider both these two classes of fluctuations (between tetrahedral and asymmetric and within asymmetric) separate but most likely there should also be some coupling.
It is these fluctuations between tetrahedral and asymmetric structures plus the varying fluctuations in the asymmetric structures that depends on temperature, pressure and interactions with solutes, interfaces, biomolecules etc. We can imagine that for instance an interface disfavors tetrahedral structures (such as in Kaya et al.) creating a dominance of asymmetric high density structures but through the interaction with the interface the switching time and amplitude will be affected between the strong and weak bonds.
The nature of the hydration shells of methanol, ethanol and propanol have been probed using THz spectroscopy by Vladimir Matvejev of the Free University of Brussels and colleagues (J. Phys. Chem. B 116, 14071; 2012 – paper here). They estimate that the shells comprise about 14, 23, 23 and 31 molecules for methanol, ethanol, 1- and 2-propanol, respectively, and the water molecules are retarded by a factor of around 1.4.
In a somewhat similar vein, L. Luca of the University of Perugia and colleagues use GHz-THz light scattering to probe the hydration of mono- and disaccharides (J. Phys. Chem. B 116, 14760; 2012 – paper here). They find that slowing of water collective reorientation (by relatively large factors of 5-6) occurs only over relatively short distances (3-4 Å or essentially the first hydration layer), regardless of the size of the sugar molecules. This retardation involves considerably more water molecules than those few instantaneously hydrogen-bonded to the sugar.
A more generic approach to small-molecule solvation is described by Alla Oleinikova and Ivan Brovchenko of the Dortmund University of Technology (J. Phys. Chem. B jp306781y – paper here). They use MC simulations to study water structure around spherical solute particles 3-10 Å in size that vary from strongly hydrophobic to strongly hydrophilic. In all cases there is a density depletion relative to the bulk due simply to the missing-neighbour effect. For strongly hydrophobic particles there is a drying transition at the surface. Similar effects are seen for other fluids, but the directional hydrogen bonding of water enhances them.
Water in useful places: there is likely to be water molecules forging a hydrogen-bonded assembly in the mammalian photoreceptor melanopsin in the retina, which triggers the biological clock, according to Sivakumar Sekharan and colleagues at Cornell (JACS 134, 19536; 2012 – paper here). Their first-principles calculations of the active site suggest that two water molecules bridge the Schiff base and residues on the peptide, accounting for the blue shift of the optical absorption relative to the closely related rhodopsin.
Cytochrome c might control water access to its heme centre to tune the reduction potential via hydration changes, according to quantum MD simulations of Isabella Daidone of the University of L’Aquila and colleagues (C. A. Bortolotti et al., JACS 134, 13670; 2012 – paper here). They find that a yeast cytochrome has two channels that it can open to admit water, altering the reduction potential and thus refining processes of electron transfer.
How water evaporates from the air-water interface has been studied in detail by Patrick Varilly and David Chandler at Berkeley (J. Phys. Chem. B jp310070y – paper here). They find that the escape trajectory of a water molecule can be described in terms of two parameters: the distance from the instantaneous interface and the velocity along the surface normal. The results seem to imply that evaporation has zero activation energy, as some but not all experiments have suggested.
Fresh variety in the structures of ice confined to nanopores is reported by Jaeil Bai and Xiao Cheng Zeng of the University of Nebraska (PNAS 109, 21240; 2012 – paper here). Their simulations indicate that a bilayer of ice-like water molecules in slit-like pores about 7-9Å apart can be transformed under pressure to an amorphous phase and then to a very-high-density amorphous phase at 250K and 3 GPa. For rapid compression to 6 GPa an entirely new VHD ordered phase is found in which the water molecules are linked into square nanotubes.
Andreas Barth and colleagues at Stockholm University propose that studying changes in the water absorption bands in infrared spectroscopy can offer a way of quickly and remotely monitoring for binding of a ligand to a target protein in drug development (S. Kumar et al., J. Phys. Chem. B 116, 13968; 2012 – paper here). Specifically, they look at the bands diagnostic of ‘bound’ water in the binding cavity, and see this decline as ligand binding expels the water. They have so far looked at rather high protein concentrations, but think that the approach should work at lower concentrations with brighter light sources.
Hydrophobic drugs can be rendered more water-soluble by adding certain solutes known as hydrotropes. Although this has been long known, the mechanism is not clear. Seishi Shimizu of the University of York and colleagues have studied the question using thermodynamic theory (the so-called fluctuation theory of solution) and measurements of thermodynamic quantities (J. J. Booth et al., J. Phys. Chem. B 116, 14915; 2012 – paper here).
Junrong Zheng at Rice University and coworkers report evidence for strong segregation of ions and water in strong electrolyte solutions of thiocyanate ions (J. Phys. Chem. B jp310153n – paper here). The clustering is enhanced by the addition of strongly hydrated ions such as fluoride, whereas iodide ions tend to associate with the SCN clusters. In any event, the aqueous systems are strongly inhomogeneous.
Finally for now – though lots more still to catch up on – I have a brief article on the Chemistry World site (here) on a recent paper by Thomas Kühne and Rustam Khaliullin at Mainz (Nat. Commun. 4, 1450; 2013 – paper here) that seems to shed light on the arguments over XAS studies of liquid water by Anders Nilsson and colleagues (Wernet et al., Science 304, 995; 2004). Views on whether there was a controversy here in the first place will doubtless vary, but it seems that this new work provides a useful perspective on what the XAS work was showing. Anders gave me some rather extensive remarks on his view of the matter, which I thought might be usefully reported here. Forgive me for presenting them here unmediated – they will hopefully at least clarify Anders’ current position on this. Everything that follows is from him:
I find this paper to represent a very important step forward. It is very close to our original suggestion in the Wernet et al. Science paper that many molecules will be in instantaneous configurations with one strong and one weak hydrogen bond, i.e. asymmetric configurations. XAS measures the electronic structure on a timescale of a few attoseconds which means that the molecules have no time to move so we are only detecting snapshots of space-averaged frozen configurations (also stated in the original Wernet et al. paper).
However, there is a difference. In real water we expect the asymmetry between the strong and weak bonds to be much larger than based on the structures from the current simulation. It is mentioned in the methods section but not directly in the main text. Here is the underlying experimental evidence that requires no spectroscopy interpretation.
I attach a recent study (JCP in press, Skinner et al.) where we together with Benmore’s group have undertaken a major effort to determine a more accurate O-O pair distribution function (PDF) of ambient water. This is based on 4 new x-ray diffraction data sets (only 3 are consistent at high Q), all with a much higher Q cut-off than any previous measurements where we have also taken serious care to remove any OH contributions. Figure 9a shows the O-O PDF with very small error bars with a first peak height of 2.57. If you comparethis with the O-O PDF presented in the supplementary material it is clear that the simulations have an overstructured main peak. For the PBE functional, used in the main text, the peak height is 3.3. The TPSS –D3-FF simulation shown in the supplementary information has a peak height of 3.1, still overstructured but less so than the PBE. The consequence in the asymmetry is clearly visible in terms of the weak bond energy distribution from fig. 3a and S3a where the latter has more contributions towards lower energies. If we have to further dramatically understructure the liquid down to a peak height of 2.57 in the O-O PDF we expect the asymmetry to become much larger. In my opinion their asymmetry is only a lower limit.
Most likely the underlying reason for the asymmetry is in the many-body cooperativity effects that only an electronic structure simulation can capture as demonstrated in the current paper, and which is not represented by classical force fields. It has been known that the cooperativity effect is strongest when you have one strong donor and one strong acceptor bond. Thereby for 2 hydrogen bonded structures the energy per hydrogen bond is higher than for 4 bonded or tetrahedral coordinated structures. We discussed this in a few sentences in the Wernet et al. paper with references that I have underlined with yellow in the attachment. I also attach a recent review by Lars, Congcong and myself on water (A. Nilsson et al., J. Mol. Liq. 176, 2; 2012) where we discuss the importance to further develop simulations and in particular the importance of many-body effects in an electronic structure description where it is also essential that the latter includes van der Waals interactions.
Another interesting aspect in their study is the dynamics indicating that the strong and weak bonds switch place on ultrafast time scales. This is also in line with our previous discussions. I underline in the Leetma et al. paper (J. Chem. Phys. 129, 084502; 2008) how we discussed ultrafast measurements where the strongly H-bonded OH group in one water molecule could switch places with the weak one. The difference was only that we assumed that the switching occurred via librations whereas in their simulations it is mostly via translational motions of other molecules.
We had a last sentence in our Wernet et al. paper of a more speculative nature that with mostly only one strong donor and one strong acceptor bond per molecule, one-dimensional structures should appear which could be chains or rings. I noticed that they have a sentence regarding molecular chains at each instant connecting the strong bonds. Very nice.
Another point not related to the current paper that could also be of interest for you is the recent paper by Overduin and Patey (J. Phys. Chem. B 116, 12014; 2012) that discusses our PNAS Huang et al. paper based on simulations coming to similar conclusions regarding inhomogeneities in the liquid but only using a different language in terms of concentration fluctuations of two different classes of structures.
I think Kühne and Khaliullin put the picture in terms of a symmetry breaking in order to be not too far away from the most accepted tetrahedral picture of water. But if you really look at what they claim, which is also aligned with us, is that most molecules will be in an asymmetric position at all instances when taking a snapshot. It is simply that an OH group that is either weak or strong switches on a rather fast timescale. This would mean that the symmetric position is never really seen. It is only when you average over a long time that it looks like it is mostly in the tetrahedral position. Like me make analog. Take a pendulum that swings back and forth. The speed at the end points are close to zero and at the middle position it is highest. It means that the pendulum spends most time at the end points and extremely little time at the middle position. If you take the average position it will be the middle position but it is hardly ever visited. You could in such a picture claim that the equilibrium position is in the middle if you average over a single period but the pendulum will spend very little time there. The question is how the surrounding will be affected by such a motion. If I understand correctly the paper it is the surrounding that indeed infer the asymmetry. They claim that it is the translational motion of the other molecules that provides the mechanism of the switching. This should mean in my opinion that it is the end points of the asymmetric motion that makes the interaction with the surrounding and not the time averaged position. If it would have been dominated by a more internal motion such as librations it could have had less effect on the surrounding.
In water my own belief is that it is somewhat more complicated. Let me come to my current picture of water at the end of these comments.
I think it is currently not possible to observe the asymmetry based on the current status of pump-probe IR spectroscopy. Please don't quote me on this but I believe that we are still missing some major understanding about IR spectroscopy. The vibrational life time as a local oscillator in H20 is too short. Nearly all experiments on ultrafast dynamics are based on isotope substitution with HDO impurities in either H2O and D2O and measurements of the time resolved development of the OH or OD frequency. It is then assumed that HDO has equal probabilities for all molecular positions in the liquid. This assumption has never been proven. Based on our understanding of water in terms of fluctuations of two local structures, tetrahedral (low density water) and asymmetric (high density water) the latter has more contributions to the free energy through stronger bond energy whereas the latter has more entropy. I attached a slide from one of my presentations so you can see what I mean (hit a return in slide show mode to see the motion). There you also see the switching of the two bonds. Since the asymmetric configuration will provide more entropy and the HDO is already asymmetric it is likely that HDO will reside more in asymmetric configurations. Furthermore, we can anticipate through quantum effects that also the OD and OH hydrogen bonds are different. I would assume that OD will more likely be situated with the strong bond and OH more in the weak bond. This was observed in a recent PRL by Misha Bonn's group (attached) where they showed based on sum frequency generation spectroscopy and simulations that there is a strong preference for OH to point to the vacuum and OD to the liquid side for HDO at the water-air interface. At this point I think caution needs to be exercised regarding pump probe IR since these measurements might only be probing some fraction of all possible molecular motions. We (Lars) are currently further investigating the HDO dynamics based on quantum simulations. There are also recent indications that the pump pulse disturbs the dynamics and these measurements might not always represent the equilibrium motions.
We are currently using the new x-ray laser at Stanford, the Linac Coherent Light Source (LCLS) to various problems related to water. The machine provides completely new opportunities. We are planning a future step to open up dynamics. In particular x-ray correlation spectroscopy with x-ray lasers could provide new answers. This is probably 2 years away. Since the x-ray laser is fully coherent there will be a speckle pattern due to diffraction also from a disordered material. The plan is to make probe-probe measurement to thereby study equilibrium dynamics. You split the x-ray beam into two pulses with a controllable delay and the change in the speckle pattern will give information about how molecules have moved between the two pulses. These will be most challenging experiments but hopefully can be done and would open for completely new avenue's to probe dynamics in liquids. Naturally one challenge will be to not allow the first pulse to disturb the system but water is quite forgiving for hard x-rays since it is a low Z liquid.
We have demonstrated in a number of papers that this asymmetry can't been seen in diffraction since the data is not completely perfect. The difficulty is to have a technique that is only sensitive to the hydrogen bonding asymmetry around individual molecules. You can have a total hydrogen bonding average through a linear combination of water molecules with nearly no hydrogen bonds and fully tetrahedral water that on average would be like an asymmetric species. It is the rehybridization of the molecular orbital structure leading to local OH orbitals instead of delocalized H2O orbitals that makes XAS sensitive. I attach a paper where we discuss this (Nilsson section 3.5). It has to do with O2p and O2s hybridization and XAS only provides intensity for the O2p part in the orbitals. It is quite involved. Maybe it is time for a simpler review so many can understand these effects.
Here is our current picture of water. At high temperatures (close to boiling and above) water behaves as a simple liquid where most interactions are isotropic (dominated by van der Waals interactions). This is a structure where many molecules are in very disordered shells and without having a well defined first and second shell (this structure you get from ab initio MD simulations with van der Waals functionals even at ambient temperature, see discussion in my previous mail in paper "fluctuations in ambient water"). This type of non directional orientations gives high flexibility for various motions and high entropy. As the liquid cools down the water molecules starts to stick to each other through directional hydrogen bonds. This will appear in two classes of configurations, tetrahedral and asymmetric. In the tetrahedral structures each molecules are in 4 hydrogen bonds, this provides the lowest enthalpic energy. Since cooperatively effects makes the bonds stronger if water is bonded to other waters that are also in tetrahedral structures, these start to clump together in small local regions. Since the molecules are stuck with four bonds the motion is very restricted (see the attached slide) and thereby low entropy.
The other alternative is to form asymmetric structures. Here the hydrogen bond energy per molecule is higher than in the tetrahedral structures but since it is fewer bonds the enthalpic energy is lower. With less directional hydrogen bonds you have more flexibility for motions and thereby higher entropy (see slides). In this structure there are also non isotropic molecular interactions causing interstitial molecules. These are therefore called by us and others, preferable in the supercooled community, as high density water. Another way of viewing such species is that we start not with hexagonal ice but with high density amorphous or very high density amorphous ices where we have a large collapse of the second shell. In these ices we have also a local tetrahedral bonding but with other angles towards the second shell (often call interstitials). Here we induce asymmetric distortions around these tetrahedral bonded molecules in the first shell but simultaneous keeping the interstitials ( see fluctuations in ambient water). We have also seen this in water at an interface recently published in the Nature journal Scientific Reports (see attached Kaya et al.). There are fluctuations between the tetrahedral structures and asymmetric or high density structures. As we cool the liquid down the molecules in the asymmetric structures converts more and more into the tetrahedral structures which grows in size since the enthalpic energy contributes more to the free energy with decreasing temperature. The timescales in the fluctuations between these two classes should also slow down (not yet determined). There is also a continuous change in the asymmetric structures with temperature. The switching time is expected to slow down and thereby also the amplitude in the motion with decreasing temperature. We will be approaching more and more a local tetrahedral arrangement. However, not as in hexagonal ice but more towards high density ices with still many interstitials. That the density is decreasing below 4C is simply due to that we are at the same time converting many molecules into tetrahedral structures which have a more similar structure as hexagonal ice but disordered more towards low density amorphous ice (low density water) with the second shell at the tetrahedral angle causing an open network with low density. In my opinion we have to consider both these two classes of fluctuations (between tetrahedral and asymmetric and within asymmetric) separate but most likely there should also be some coupling.
It is these fluctuations between tetrahedral and asymmetric structures plus the varying fluctuations in the asymmetric structures that depends on temperature, pressure and interactions with solutes, interfaces, biomolecules etc. We can imagine that for instance an interface disfavors tetrahedral structures (such as in Kaya et al.) creating a dominance of asymmetric high density structures but through the interaction with the interface the switching time and amplitude will be affected between the strong and weak bonds.
Thursday, December 6, 2012
Salts: how far do they go?
In my previous post I commented on a recent paper by Jeremy Smith and colleagues (L. Hong et al., Phys. Rev. Lett. 108, 238102; 2012) in which it was shown that the coupling between the dynamics of a protein and its hydration shell can propagate from the surface to the core to soften up the entire molecule. Martin Weik at Grenoble, Frans Mulder at Aarhus and their coworkers have now investigated this idea experimentally via neutron scattering from the small protein calbindin D9k P43G deuterated in the one case on the “outside” and in the other case on the “inside” (K. Wood et al., Angew. Chem. Int. Ed. 10.1002/anie.201205898 – paper here). They find that both the exterior and the core of the protein are sensitive to hydration: the dynamics of both regions undergo a hydration-dependent dynamical transition in the same temperature range of around 250 K. In other words, this transition does appear to be a “global” one.
Ruhong Zhou at IBM Yorktown Heights is continuing to investigate how denaturants affect protein conformation and folding. In a paper with Eugene Shaknovich and coworkers, he reports the surprising finding that in a mixture of the denaturants urea and guanidinium chloride, hen egg-white lysozyme and protein L will both collapse, albeit with an increase in non-native hydrophobic contacts (Z. Xia et al., JACS 134, 18266; 2012 – paper here). The collapse (it is evidently not ‘folding’ in the proper sense) is induced by the specific interactions of the denaturants with the protein surface: GdmCl is absorbed there due to electrostatic interactions, while urea also accumulates near the first hydration shell and introduces crowding. This is a reminder of how specific, rather than generic, interactions of proteins with denaturants and osmolytes will dominate the behaviour, sometimes in non-intuitive ways.
Another counter-intuitive result is reported by Haiping Fang of the Shanghai Institute of Applied Physics and coworkers, relating to the evaporation of water molecules from solid surfaces. They find that the evaporation rate doesn’t change monotonically as one progresses from hydrophilic to hydrophobic surfaces, as one might expect, but has a maximum (S. Wang et al., J. Phys. Chem. B 116, 13863; 2012 – paper here). This follows from the fact that, while interactions between water and the surface dominate in situations where the surface is well wetted, for highly hydrophobic surfaces the water gathers into isolated ‘surface droplets’, in which case water evaporation is controlled by water-water interactions. The maximum results from a crossover between these two competing mechanisms. The findings might be important not only for, e.g. water retention in soils but also for the possibility of drying transitions in hydrophobic aggregation.
It’s become clear that even in relatively dilute solution some soluble small organic molecules may not be homogeneously dispersed. That notion is backed up by neutron-scattering experiments of Lorna Dougan at the University of Leeds and colleagues, in which they investigate glutamine solutions (N. H. Rhys et al., J. Phys. Chem. B 116, 13308; 2012 – paper here). Polyglutamine stretches of proteins are quite common and apparently important – they feature, for example, in some proteins that aggregate to induce neurodegeneration. Such aggregates seem to have collapsed polyglutamine regions, even though glutamate might be expected to form favourable hydrogen-bonding interactions with water. This motivated the present study, to investigate the interactions of glutamine monomers in water. It appears that glutamine forms some dimers via hydrogen bonds in both the ‘backbone’ and ‘side chain’ of the molecules, even at concentrations of 30 mg per mL, revealing a strong propensity to self-associate. Later work will address glutamine oligomers and polymers.
Lorna, in collaboration again with Alan Soper, has also used neutron scattering to investigate the clustering/microsegregation of the cryoprotectant glycerol in water (J. J. Towey et al., J. Phys. Chem. B 116, 13898; 2012 – paper here). They find that glycerol’s cryobiological effects don’t seem to stem from any effect it has on the hydrogen-bonding ability of water. Rather, glycerol molecules simply replace water molecules to allow the waters to retain their hydrogen-bonding capacity over a wide range of glycerol concentration. However, the presence of glycerol does segregate water into clusters at higher concentrations. So the researchers propose that the cryoprotectant activity stems not from any disruption of ‘water structure’ per se, but from glycerol’s very ability to substitute for water, while suppressing ice formation because of water segregation.
It is sobering to see that the dissociation of an ion pair such as NaCl –that’s to say, how and why salt dissolves – is still not fully understood. Andrew Ballard at the University of Maryland and Christoph Dellago of the University of Vienna present simulations of this process, using TIP4P water (J. Phys. Chem. B 116, 13490; 2012 – paper here). They say that the ion dissociation is favoured energetically but opposed entropically because of the water molecules entering the hydration shells of the ions. As in Dellago’s earlier work with Philip Geissler and David Chandler (J. Phys. Chem. B 103, 3706 (1999)), they find that the inter-ion separation is not a good reaction coordinate, as structures with different relaxation behaviour can occur for the same separation. The solvent fluctuations continue to play a role in the process over relatively long ranges, even into the third hydration shell.
It would be nice to know how that finding of a relatively long-ranged perturbation of water sits with the apparently surprising results of Sheeba Jem Irudayam of the UNC at Chapel Hill and Richard Henchman at Manchester. Their MD simulations of the hydration of alkali metal and halide ions show that this perturbation extends over remarkably long distances (J. Chem. Phys. 137, 034508; 2012 – paper here). They find non-bulk-like structural characteristics over virtually the whole simulation box (about 2.2 nm on a side for most ions, 3.3 nm for lithium and fluoride, corresponding to around 375 and 1200 TIP4P water molecules respectively). Specifically, there is a very slight but detectable excess of H-bond acceptors around halide ions, and of donors around alkali metal ions, which compensates for the shorter-ranged perturbations to water structure in the first and second hydration spheres. They also find long-ranged deviations for noble gas solutes: slightly enhanced tetrahedrality and an oscillating excess and deficiency of donors and acceptors. Similar perturbations are found for the air-water interface. The authors say that such small effects are likely to be invisible to standard spectroscopic techniques, although there have been some hints of them in neutron-scattering studies by Alan Soper and coworkers. To what extent they might have important thermodynamic consequences remain to be seen, but these results are highly intriguing.
In concentrated salt solution, ion-dependent effects on the O-D stretch are seen for both the anion and cation in pump-probe IR spectroscopic experiments by Michaal Fayer and colleagues at Stanford (C. H. Giammanco et al., J. Phys. Chem. B 116, 13781; 2012 – paper here). They find that in concentrated solutions it no longer suffices to divide water’s relaxational modes into ion-associated and water-associated fractions, because there is no longer any bulk-like component of the solvent. In this situation, water reorientational motions are highly cooperative.
Ionic hydration in confined spaces is important for a number of reasons not connected to water in biology, not least for understanding the interfacial behaviour in supercapacitors and batteries with nanoporous carbon electrodes. Tomonori Ohba of Chiba University in Japan and colleagues have explored this issue using synchrotron XRD of electrolytes inside carbon nanotubes (JACS 134, 17850; 2012 – paper here). They conclude that hydration is significantly different inside nanotubes with an average internal pore diameter of 2 nm, relative to the bulk. They can’t evaluate hydration numbers or detailed hydration structures, but say that the hydration structuring is stronger under confinement and that the hydrogen-bonded network of the solvent is correspondingly stretched and weakened. I must confess that I struggle to find an intuitive picture of what is happening here – but this perturbation is at least consistent with experiments showing a pore-size dependence of capacitance in double-layer capacitors (e.g. Chmiola et al., Science 313, 1760; 2006).
Transport of water through nanopores such as carbon nanotubes and aquaporin has important implications both for biology (e.g. functioning of ion channels) and technology (water purification). Kuiwen Zhao and Huiying Wu of Shanghai Jiao Tong University have used MD simulations to study water and ion transport through arrays of carbon nanotubes, driven by osmotic pressure (J. Phys. Chem. B 116, 13459; 2012 – paper here). In particular, they have looked at the effect of the packing density of the pores, and find that at high packing densities there can be steric interference of ions (and their hydration spheres) entering the pore mouths. There may therefore be an optimal packing density for efficient transport through the pore array, rather than simply trying to pack pores as densely as possible.
Acid or base? Yes, it’s the air-water interface again, and this time Agustín Colussi at Caltech and colleagues report experiments which seem to indicate that the interface is Brønsted basic (H. Mishra et al., PNAS 109, 18679; 2012 – paper here). They say that the controversies and discrepancies that have previously plagues this question stem at least partly from a failure to recognize acidity as a relative concept referring not so much to proton or hydroxide concentrations as to the extent of proton sharing between conjugate acid/base pairs. The researchers use electrospray ionization mass spectrometry of interfacial layers to measure the degree of dissociation of carboxylic acids both in the dissolved aqueous phase and when it collides in the gas phase with a water jet – probing, respectively, the ‘inner’ and ‘outer’ side of the surface. The detection of carboxylate ions indicates the presence of hydroxide at the surface, for all pH>2. Moreover, this surface excess of hydroxide can account for the observed negative charge at the air-water surface.
Mischa Bonn and coworkers at Mainz say that the bond orientational behaviour of water at the air-water interface is sensitive to nuclear quantum effects (Y. Nagata et al., Phys. Rev. Lett. 109, 226101; 2012 – paper here). They report quantum MD simulations showing that, while H2O and D2O have indistinguishable structures at the interface, HDO is quite distinct, with OD bonds oriented into the liquid and OH bonds oriented towards the gas phase. This would be because OD groups are able to form relatively stronger hydrogen bonds owing to a quantum isotope effect. They say this finding is in good quantitative agreement with SFG spectroscopic studies, e.g. by Geri Richmond.
The use of fluorescence microscopy to study freeze-dried biological samples can reveal details of water and ionic content of cells at the sub-cellular, nanoscale level, according to a paper by Jean Michel at the Université de Reims Champagne Ardenne in France and coworkers (F. Nolin et al., J. Struct. Biol. 180, 352; 2012 – paper here). Elemental (e.g. ion) distributions can be deduced from EDXS analysis, while specific protein densities can be studied by GFP-labelling. It looks like a nice method for deducing larger-scale patterns of hydration and ion distribution, for example that occasioned by chromatin compaction.
Water can undergo capillary evaporation from between two hydrophobic surfaces at small separations, but exactly how this happens hasn’t been fully elucidated. Sumit Sharma and Pablo Debenedetti at Princeton have investigated the process using MC simulations, and find that the outcome depends on the size of the surfaces (J. Phys. Chem. B 116, 13282; 2012 – paper here). If they are sufficiently large (3 nm square), evaporation involves the formation of a tubular cavity spanning the gap of the slit-like pore – an activated event as described by classical nucleation theory. But for 1 nm square surfaces there is too little space to accommodate such a vapour cavity, and the gap simply empties entirely.
It’s well documented that water in such confined spaces may also show reduced diffusional mobility. Hiroki Matsubara at Tohoku University in Japan and colleagues have attempted to figure out why this happens for liquids in general using MD simulations (Phys. Rev. Lett. 109, 197801; 2012 – paper here). They find that OMCTS (an alkylsiloxane) between two mica surfaces has a diffusion coefficient that depends on the surface separation, at least in the range 23-64 Å (3-7 molecular layers). This is because diffusion is an activated process, and the activation energy increases both because of the entropic constraint on some molecular configurations under confinement and the lowering of the average molecular potential energy in the gap. To what extent these effects operate and are modified for water remains to be seen.
How water mediates the association of a protein with its ligand is one of the most interesting issues concerning water’s role in molecular biology. Francesco Paesani of the University of California at San Diego and his colleagues present MD simulations which further support the notion that the water is involved as an active participant (R. Baron et al., J. Phys. Chem. B 116, 13774; 2012 – paper here). They are interested in extracting signatures of time-dependent changes in water structure and dynamics accessible to ultrafast vibrational spectroscopy. They model the binding between an apolar cavity and a spherical hydrophobic ligand, and find that the expulsion of disordered water from the cavity and suppression of slow water density fluctuations on binding result in an unfavourable entropic contribution to the binding free energy. Meanwhile, the reorientational dynamics of the water hydrating the ligand speed up as it approaches the cavity, and this water becomes less tetrahedral. This should be accompanied by the appearance of a shoulder on the O-D stretch mode for mixtures of HOD/H2O, as used in some recent spectroscopic studies, owing to the increasing concentration of dangling O-D bonds.
Francesco and his coworkers also present a new ab initio model for water, called HBB2-pol (V. Babin et al., J. Phys. Chem. Lett. 3, 3765; 2012 - paper here). They say it provides agreement with experiment from water dimers to the liquid state, capturing both the observed structural and dynamic properties while a avoiding the computational complexity that has previously rendered models such as WHBB intractable beyond small clusters.
Hydrophobicity may be the driving force behind the formation of some knots in protein structures. This idea is supported by work by Jeremy Sanders at Cambridge and colleagues, who show that a synthetic molecule composed of hydrophobic polyaromatics linked by hydrophilic amino acids will trimerize into a trefoil knot as the most effective way to ‘bury’ the hydrophobic surfaces (N. Ponnuswamy et al., Science 338, 783; 2012 – paper here).
Ruhong Zhou at IBM Yorktown Heights is continuing to investigate how denaturants affect protein conformation and folding. In a paper with Eugene Shaknovich and coworkers, he reports the surprising finding that in a mixture of the denaturants urea and guanidinium chloride, hen egg-white lysozyme and protein L will both collapse, albeit with an increase in non-native hydrophobic contacts (Z. Xia et al., JACS 134, 18266; 2012 – paper here). The collapse (it is evidently not ‘folding’ in the proper sense) is induced by the specific interactions of the denaturants with the protein surface: GdmCl is absorbed there due to electrostatic interactions, while urea also accumulates near the first hydration shell and introduces crowding. This is a reminder of how specific, rather than generic, interactions of proteins with denaturants and osmolytes will dominate the behaviour, sometimes in non-intuitive ways.
Another counter-intuitive result is reported by Haiping Fang of the Shanghai Institute of Applied Physics and coworkers, relating to the evaporation of water molecules from solid surfaces. They find that the evaporation rate doesn’t change monotonically as one progresses from hydrophilic to hydrophobic surfaces, as one might expect, but has a maximum (S. Wang et al., J. Phys. Chem. B 116, 13863; 2012 – paper here). This follows from the fact that, while interactions between water and the surface dominate in situations where the surface is well wetted, for highly hydrophobic surfaces the water gathers into isolated ‘surface droplets’, in which case water evaporation is controlled by water-water interactions. The maximum results from a crossover between these two competing mechanisms. The findings might be important not only for, e.g. water retention in soils but also for the possibility of drying transitions in hydrophobic aggregation.
It’s become clear that even in relatively dilute solution some soluble small organic molecules may not be homogeneously dispersed. That notion is backed up by neutron-scattering experiments of Lorna Dougan at the University of Leeds and colleagues, in which they investigate glutamine solutions (N. H. Rhys et al., J. Phys. Chem. B 116, 13308; 2012 – paper here). Polyglutamine stretches of proteins are quite common and apparently important – they feature, for example, in some proteins that aggregate to induce neurodegeneration. Such aggregates seem to have collapsed polyglutamine regions, even though glutamate might be expected to form favourable hydrogen-bonding interactions with water. This motivated the present study, to investigate the interactions of glutamine monomers in water. It appears that glutamine forms some dimers via hydrogen bonds in both the ‘backbone’ and ‘side chain’ of the molecules, even at concentrations of 30 mg per mL, revealing a strong propensity to self-associate. Later work will address glutamine oligomers and polymers.
Lorna, in collaboration again with Alan Soper, has also used neutron scattering to investigate the clustering/microsegregation of the cryoprotectant glycerol in water (J. J. Towey et al., J. Phys. Chem. B 116, 13898; 2012 – paper here). They find that glycerol’s cryobiological effects don’t seem to stem from any effect it has on the hydrogen-bonding ability of water. Rather, glycerol molecules simply replace water molecules to allow the waters to retain their hydrogen-bonding capacity over a wide range of glycerol concentration. However, the presence of glycerol does segregate water into clusters at higher concentrations. So the researchers propose that the cryoprotectant activity stems not from any disruption of ‘water structure’ per se, but from glycerol’s very ability to substitute for water, while suppressing ice formation because of water segregation.
It is sobering to see that the dissociation of an ion pair such as NaCl –that’s to say, how and why salt dissolves – is still not fully understood. Andrew Ballard at the University of Maryland and Christoph Dellago of the University of Vienna present simulations of this process, using TIP4P water (J. Phys. Chem. B 116, 13490; 2012 – paper here). They say that the ion dissociation is favoured energetically but opposed entropically because of the water molecules entering the hydration shells of the ions. As in Dellago’s earlier work with Philip Geissler and David Chandler (J. Phys. Chem. B 103, 3706 (1999)), they find that the inter-ion separation is not a good reaction coordinate, as structures with different relaxation behaviour can occur for the same separation. The solvent fluctuations continue to play a role in the process over relatively long ranges, even into the third hydration shell.
It would be nice to know how that finding of a relatively long-ranged perturbation of water sits with the apparently surprising results of Sheeba Jem Irudayam of the UNC at Chapel Hill and Richard Henchman at Manchester. Their MD simulations of the hydration of alkali metal and halide ions show that this perturbation extends over remarkably long distances (J. Chem. Phys. 137, 034508; 2012 – paper here). They find non-bulk-like structural characteristics over virtually the whole simulation box (about 2.2 nm on a side for most ions, 3.3 nm for lithium and fluoride, corresponding to around 375 and 1200 TIP4P water molecules respectively). Specifically, there is a very slight but detectable excess of H-bond acceptors around halide ions, and of donors around alkali metal ions, which compensates for the shorter-ranged perturbations to water structure in the first and second hydration spheres. They also find long-ranged deviations for noble gas solutes: slightly enhanced tetrahedrality and an oscillating excess and deficiency of donors and acceptors. Similar perturbations are found for the air-water interface. The authors say that such small effects are likely to be invisible to standard spectroscopic techniques, although there have been some hints of them in neutron-scattering studies by Alan Soper and coworkers. To what extent they might have important thermodynamic consequences remain to be seen, but these results are highly intriguing.
In concentrated salt solution, ion-dependent effects on the O-D stretch are seen for both the anion and cation in pump-probe IR spectroscopic experiments by Michaal Fayer and colleagues at Stanford (C. H. Giammanco et al., J. Phys. Chem. B 116, 13781; 2012 – paper here). They find that in concentrated solutions it no longer suffices to divide water’s relaxational modes into ion-associated and water-associated fractions, because there is no longer any bulk-like component of the solvent. In this situation, water reorientational motions are highly cooperative.
Ionic hydration in confined spaces is important for a number of reasons not connected to water in biology, not least for understanding the interfacial behaviour in supercapacitors and batteries with nanoporous carbon electrodes. Tomonori Ohba of Chiba University in Japan and colleagues have explored this issue using synchrotron XRD of electrolytes inside carbon nanotubes (JACS 134, 17850; 2012 – paper here). They conclude that hydration is significantly different inside nanotubes with an average internal pore diameter of 2 nm, relative to the bulk. They can’t evaluate hydration numbers or detailed hydration structures, but say that the hydration structuring is stronger under confinement and that the hydrogen-bonded network of the solvent is correspondingly stretched and weakened. I must confess that I struggle to find an intuitive picture of what is happening here – but this perturbation is at least consistent with experiments showing a pore-size dependence of capacitance in double-layer capacitors (e.g. Chmiola et al., Science 313, 1760; 2006).
Transport of water through nanopores such as carbon nanotubes and aquaporin has important implications both for biology (e.g. functioning of ion channels) and technology (water purification). Kuiwen Zhao and Huiying Wu of Shanghai Jiao Tong University have used MD simulations to study water and ion transport through arrays of carbon nanotubes, driven by osmotic pressure (J. Phys. Chem. B 116, 13459; 2012 – paper here). In particular, they have looked at the effect of the packing density of the pores, and find that at high packing densities there can be steric interference of ions (and their hydration spheres) entering the pore mouths. There may therefore be an optimal packing density for efficient transport through the pore array, rather than simply trying to pack pores as densely as possible.
Acid or base? Yes, it’s the air-water interface again, and this time Agustín Colussi at Caltech and colleagues report experiments which seem to indicate that the interface is Brønsted basic (H. Mishra et al., PNAS 109, 18679; 2012 – paper here). They say that the controversies and discrepancies that have previously plagues this question stem at least partly from a failure to recognize acidity as a relative concept referring not so much to proton or hydroxide concentrations as to the extent of proton sharing between conjugate acid/base pairs. The researchers use electrospray ionization mass spectrometry of interfacial layers to measure the degree of dissociation of carboxylic acids both in the dissolved aqueous phase and when it collides in the gas phase with a water jet – probing, respectively, the ‘inner’ and ‘outer’ side of the surface. The detection of carboxylate ions indicates the presence of hydroxide at the surface, for all pH>2. Moreover, this surface excess of hydroxide can account for the observed negative charge at the air-water surface.
Mischa Bonn and coworkers at Mainz say that the bond orientational behaviour of water at the air-water interface is sensitive to nuclear quantum effects (Y. Nagata et al., Phys. Rev. Lett. 109, 226101; 2012 – paper here). They report quantum MD simulations showing that, while H2O and D2O have indistinguishable structures at the interface, HDO is quite distinct, with OD bonds oriented into the liquid and OH bonds oriented towards the gas phase. This would be because OD groups are able to form relatively stronger hydrogen bonds owing to a quantum isotope effect. They say this finding is in good quantitative agreement with SFG spectroscopic studies, e.g. by Geri Richmond.
The use of fluorescence microscopy to study freeze-dried biological samples can reveal details of water and ionic content of cells at the sub-cellular, nanoscale level, according to a paper by Jean Michel at the Université de Reims Champagne Ardenne in France and coworkers (F. Nolin et al., J. Struct. Biol. 180, 352; 2012 – paper here). Elemental (e.g. ion) distributions can be deduced from EDXS analysis, while specific protein densities can be studied by GFP-labelling. It looks like a nice method for deducing larger-scale patterns of hydration and ion distribution, for example that occasioned by chromatin compaction.
Water can undergo capillary evaporation from between two hydrophobic surfaces at small separations, but exactly how this happens hasn’t been fully elucidated. Sumit Sharma and Pablo Debenedetti at Princeton have investigated the process using MC simulations, and find that the outcome depends on the size of the surfaces (J. Phys. Chem. B 116, 13282; 2012 – paper here). If they are sufficiently large (3 nm square), evaporation involves the formation of a tubular cavity spanning the gap of the slit-like pore – an activated event as described by classical nucleation theory. But for 1 nm square surfaces there is too little space to accommodate such a vapour cavity, and the gap simply empties entirely.
It’s well documented that water in such confined spaces may also show reduced diffusional mobility. Hiroki Matsubara at Tohoku University in Japan and colleagues have attempted to figure out why this happens for liquids in general using MD simulations (Phys. Rev. Lett. 109, 197801; 2012 – paper here). They find that OMCTS (an alkylsiloxane) between two mica surfaces has a diffusion coefficient that depends on the surface separation, at least in the range 23-64 Å (3-7 molecular layers). This is because diffusion is an activated process, and the activation energy increases both because of the entropic constraint on some molecular configurations under confinement and the lowering of the average molecular potential energy in the gap. To what extent these effects operate and are modified for water remains to be seen.
How water mediates the association of a protein with its ligand is one of the most interesting issues concerning water’s role in molecular biology. Francesco Paesani of the University of California at San Diego and his colleagues present MD simulations which further support the notion that the water is involved as an active participant (R. Baron et al., J. Phys. Chem. B 116, 13774; 2012 – paper here). They are interested in extracting signatures of time-dependent changes in water structure and dynamics accessible to ultrafast vibrational spectroscopy. They model the binding between an apolar cavity and a spherical hydrophobic ligand, and find that the expulsion of disordered water from the cavity and suppression of slow water density fluctuations on binding result in an unfavourable entropic contribution to the binding free energy. Meanwhile, the reorientational dynamics of the water hydrating the ligand speed up as it approaches the cavity, and this water becomes less tetrahedral. This should be accompanied by the appearance of a shoulder on the O-D stretch mode for mixtures of HOD/H2O, as used in some recent spectroscopic studies, owing to the increasing concentration of dangling O-D bonds.
Francesco and his coworkers also present a new ab initio model for water, called HBB2-pol (V. Babin et al., J. Phys. Chem. Lett. 3, 3765; 2012 - paper here). They say it provides agreement with experiment from water dimers to the liquid state, capturing both the observed structural and dynamic properties while a avoiding the computational complexity that has previously rendered models such as WHBB intractable beyond small clusters.
Hydrophobicity may be the driving force behind the formation of some knots in protein structures. This idea is supported by work by Jeremy Sanders at Cambridge and colleagues, who show that a synthetic molecule composed of hydrophobic polyaromatics linked by hydrophilic amino acids will trimerize into a trefoil knot as the most effective way to ‘bury’ the hydrophobic surfaces (N. Ponnuswamy et al., Science 338, 783; 2012 – paper here).
Subscribe to:
Posts (Atom)