How homogeneous is the cytoplasm? There is increasing evidence that proteins in concentrated solution form relatively long-lived clusters. Wilson Poon and his colleagues showed in 2004 that this effect, previously rather anecdotal, is real and general, applying to colloidal particles as well as proteins (Stradner et al., Nature 432, 492; 2004). They showed using small-angle X-ray and neutron scattering that lysozyme forms clusters of about 3-10 molecules at volume fractions of between 0.05 and 0.2, which they say are equilibrium structures resulting from the interplay of short-ranged attractive (van de Waals) and electrostatic repulsive forces. Now Peter Vekilov at Houston and his coworkers broadly support that notion by looking at concentrated solution of bacterial lumazine synthase using light scattering (Gliko et al., J. Phys. Chem. B 111, 3106; 2007). They see clusters with lifetimes of around 10 s and a mean radius of about 350 nm (individual molecules are bout 15.6 nm in diameter). Changing the protein concentration changes the cluster concentration (which can reach a volume fraction of 0.001), but not the cluster size. But Velikov et al. say that cluster formation and size is dominated by kinetics, not thermodynamics: these clusters are metastable with respect to both protein crystals (i.e. in supersaturated solution) and to well-dispersed solution.
All of this recalls the flurry of interest in clustering excited by the work of Geckeler and Samal on C60 (Chem. Comm. 2001, 2224), which was rather breathlessly touted as a possible mechanism for homeopathy. That work was truly odd, as it seemed to suggest that the cluster size increased with increasing dilution. I’m not aware that the result has been reproduced. Needless to say, the homeopathy connection makes no sense (at best, you get a few ‘active’ bottles and the vast majority containing just water); but at the very least, there’s no reason to regard clustering per se as perplexing or odd.
The question of a vapour gap at the interface of water and a hydrophobic surface, and its relation to the long-ranged hydrophobic attraction, seems to be resolving itself. It seems now that a very thin depletion layer exists. But the role of dissolved gases in forming nanobubbles remains to be fully resolved (see ‘Why does water do that’ below). They have been seen by various methods, but with the proviso that they could possibly be an artefact of the probe technique, and that they haven’t been shown definitively to be gaseous anyway. Also, small nanobubbles should have a small radius of curvature and thus a large Laplace pressure, promoting their dissolution. William Ducker and colleagues at the University of Melbourne have now shown that flat gas bubbles, about 5-80 nm thick and 4 microns across, can exist at such a hydrophobic interface for over an hour (Phys. Rev. Lett. 98, 136101; 2007). This size means that the internal pressure is barely above atmospheric. But the bubbles form only when a particular protocol is followed for introducing the gas layer (carbon dioxide): in other words, “the presence of the gas phase depends on the previous history of the interface.”
Ivan Brovchenko and colleagues recently linked the low-hydration polymorphic transitions of B-DNA to the presence (or not) of a fully connected (percolating) network of water molecules in the hydration sphere (Brovchenko et al., Phys. Rev. Lett. 97, 137801; 2006). They’ve now extended that work by looking at the percolation transition for both B- and A-DNA (Brovchenko et al., J. Phys. Chem. B 111, 3258; 2007). Although the percolation thresholds (that is, the surface coverage of water on the DNA molecules) are virtually identical, the mechanisms are quite different in each case: the threshold corresponds to the appearance of a spanning water network in the major groove of B-DNA, but the minor groove of A-DNA. It isn’t clear, then, whether the near-coincidence of the two thresholds is indeed just coincidence or has some deeper physical cause. In any event, there are also insights here into how ions can alter the hydrogen-bonding patterns and thus shift the thresholds.
Water seems to play an important role in electron transfer between some protein redox centres – this was discussed nicely by Gray & Winkler recently (PNAS 102, 3534; 2005). They enumerated the various ways, direct and indirect, that water might facilitate electron hopping. Agostino Migliore and colleagues in Modena have now used ab initio calculations for the copper active sites of azurin to figure out how water-mediated pathways are functioning in this case (Migliore et al., J. Phys. Chem. B, advance online publication, doi:10.1021/jp068773i). But it’s a complicated picture that emerges, in which no one physical pathway seems to be responsible for what is observed. Sorry, but it’s hard to put this one into more of a nutshell than that.
Finally, and perhaps even more cryptically, I want to flag up a paper by Bruce Berne and colleagues in JACS ASAP (doi:10.1021/ja068305m) on the effect of ions on the hydrophobic interaction between two plates. This complements Berne’s recent study of much the same thing for hydrophobic particles (J. Phys. Chem. B 110, 22736; 2006). The phenomenon is of course intimately related to salting-in/out and Hofmeister effects, and as such, contains a lot of important information on the effect of electrolytes on protein aggregation and folding. So there’s a lot of good stuff here to digest, and I’m not going to manage that in a hurry without risking indigestion (or more probably, misapprehension). Worth spending time on.
Tuesday, March 27, 2007
Tuesday, March 20, 2007
How thick-skinned is hydration?
How much water do you need to fully solvate a protein? There are many studies of protein behaviour at low water coverage, going back to the suggestion by Rupley and Careri (Adv. Protein Chem. 41, 37-172; 1991) that proteins seem to require about 0.4 g of water per gram of protein to achieve their normal functionality. Roland Winter and coworkers have investigated the notion of a percolation transition in water coverage that brings the protein dynamics to life (Oleinikova et al., J. Phys. Chem. B 109, 1988-1998; 2005; Smolin et al., J Phys. Chem. B 109, 10995-11005; 2005). Mehdi Bagheri Hamaneh and Matthias Buck at Case Western have looked at the question in a rather different light: how much water do you need to put around a protein in order to be able to simulate it realistically? They find (Biophys. J. 92, L49; 2007) that you don't need to fill up your simulation box with explicit water – a shell just two or three layers thick (using the CHARMM22/CAMP potential function) will do the job well enough. That's computationally cheap, and I suppose implies that there's not really much excuse for failing to model hydration explicitly. It also implies that the celebrated cooperativity of water dynamics in the hydration shell does not appear to extend very far – at least, perhaps one should say, for the case of lysozyme considered here.
More on the 'glass transition' at around 220 K, seemingly shown now (Chen et al., PNAS 103, 901; 2006) to be a fragile-to-strong crossover. This applies also to DNA (Chen et al., J. Chem. Phys. 125, 171103; 2006), and now Sow-Hsin Chen at MIT and coworkers have found similar behaviour for RNA, again at 220 K (http://xxx.arxiv.org/abs/physics/0703166). So this seems to be pretty universal behaviour for biomacromolecules, reinforcing the idea that the change in dynamics is imposed by the hydration water.
Masahiro Kinoshita has sent me a couple of his papers from Chem. Phys. Lett. (see them here and here) which explore his idea that "the major driving force in protein folding is a gain in water entropy". In a nutshell, they say that "a protein is designed to fold into the structure that maximizes the entropy of water under the requirement that sufficiently many intramolecular hydrogen bonds be formed to compensate the dehydration penalty." In other words, as I understand it, the enthalpies balance and what's left (governing stability) is the water entropy change. That's intriguing; I'm still struggling to see how this ties up with Jack Dunitz's suggestion (Science 264, 670; 1994; Chem. Biol. 2, 709-712; 1995) that transferring a water molecule from an ordered binding site where it is bound by an ‘average’ hydrogen bond to the bulk involves an overall free-energy change that is close to zero, and with ideas about the importance of interactions at specific residues – the 'hotspots' discussed in the last post, and Ariel Fernandez's notion of dehydrons, for instance. One day, perhaps, it will all make sense to me.
More on the 'glass transition' at around 220 K, seemingly shown now (Chen et al., PNAS 103, 901; 2006) to be a fragile-to-strong crossover. This applies also to DNA (Chen et al., J. Chem. Phys. 125, 171103; 2006), and now Sow-Hsin Chen at MIT and coworkers have found similar behaviour for RNA, again at 220 K (http://xxx.arxiv.org/abs/physics/0703166). So this seems to be pretty universal behaviour for biomacromolecules, reinforcing the idea that the change in dynamics is imposed by the hydration water.
Masahiro Kinoshita has sent me a couple of his papers from Chem. Phys. Lett. (see them here and here) which explore his idea that "the major driving force in protein folding is a gain in water entropy". In a nutshell, they say that "a protein is designed to fold into the structure that maximizes the entropy of water under the requirement that sufficiently many intramolecular hydrogen bonds be formed to compensate the dehydration penalty." In other words, as I understand it, the enthalpies balance and what's left (governing stability) is the water entropy change. That's intriguing; I'm still struggling to see how this ties up with Jack Dunitz's suggestion (Science 264, 670; 1994; Chem. Biol. 2, 709-712; 1995) that transferring a water molecule from an ordered binding site where it is bound by an ‘average’ hydrogen bond to the bulk involves an overall free-energy change that is close to zero, and with ideas about the importance of interactions at specific residues – the 'hotspots' discussed in the last post, and Ariel Fernandez's notion of dehydrons, for instance. One day, perhaps, it will all make sense to me.
Tuesday, March 13, 2007
Why are hot spots hot?
In protein-protein interactions, some residues do nearly all the work. These 'hot spots' are apparently responsible for most of the binding energy, something that becomes apparent if they are mutated to alanine. It has been suggested that hot spots rely on being dry – sheltered from water by a surrounding ring of protective residues, called the O ring. But that idea hasn't been well tested. Now Maria Ramos and coworkers at the University of Porto in Portugal have studied it using MD simulations (see paper here). They find that indeed hot (and 'warm') spots in the interaction of an immunoglobulin and a lysozyme do seem to be relatively inaccessible to water. Moreover, the water molecules that do get past the O ring have relatively short residence times, much the same as those of bulk water – they don't particularly want to be there.
Greg Voth and his coworkers have a paper in JACS looking at proton transport in cytochrome c oxidase. It has been suggested that a proton is held in a 'trap' in the bovine form of this enzyme before being transported once a residue elsewhere is deprotonated. Voth and colleagues show that this process happens for a bacterial cytochrome c oxidase (Rhodobacter sphaeroides) too, and that it depends on a hydrogen-bonded network of about five water molecules that is somewhat comparable to the proton-release and transport complex found in bacteriorhodopsin (see Garczarek & Gerwert, Nature 439, 109; 2006).
Greg Voth and his coworkers have a paper in JACS looking at proton transport in cytochrome c oxidase. It has been suggested that a proton is held in a 'trap' in the bovine form of this enzyme before being transported once a residue elsewhere is deprotonated. Voth and colleagues show that this process happens for a bacterial cytochrome c oxidase (Rhodobacter sphaeroides) too, and that it depends on a hydrogen-bonded network of about five water molecules that is somewhat comparable to the proton-release and transport complex found in bacteriorhodopsin (see Garczarek & Gerwert, Nature 439, 109; 2006).
Sunday, March 4, 2007
Water from first principles?
A paper in Science this week claims to present a new, improved effective pair potential for water calculated from first principles, which works well both for water dimers and for the bulk liquid. I was a little puzzled by this, since there has been plenty of previous work in this area and it wasn’t immediately clear what the new trick was here that had enabled the claimed improvements. Having been asked to write about the work for Chemistry World, I contacted two experts in this field. Both were strongly critical of the paper.
First, the general context, which is nicely explained by one of my advisers:
“Modeling intermolecular interactions in water in terms of simple pair potentials is difficult because these potentials miss the cooperative effects of the H bonds. These are important: it is because of these effects that the “dipole” moment of a molecule in condensed phase is greatly enhanced compared to that of a molecule in gas phase (here I put dipole in quotes because this quantity cannot be unambiguously defined in condensed phase). Compared to simple liquids water has a much more open structure dictated by the underlying H-bond network.
Popular empirical potentials for water try to model the basic physics underlying H bonds interactions, which are largely of electrostatic origin in terms of interactions between point charges spatially located in a way that mimics, albeit very approximately, the charge distribution in the water molecule. These potentials are tuned to reproduce a number of properties of the liquid but have limited transferability to different environments as they miss the cooperative effects mentioned above.”
Existing classical potentials for water describe well the structure of the liquid (even very well I would say!). The difficulty with this approach, apart from the issue of transferability, is that we also need to describe the interaction between water and other molecules, material systems etc. in order to model solvation, interfacial water, confined water etc. In addition, there are properties such as the dielectric properties that depend on the actual electronic structure of the molecules in a given environment. Furthermore water molecules can dissociate. A good potential for the intermolecular interaction does not tell us everything.”
But in terms of the Science paper by Bukowski et al., he says:
“I am not too impressed by the contribution of Bukowski et al. They show that an extremely good potential for the dimer (a pair potential) is not sufficient for a good description of the liquid, as one could have expected. Including non-additive many body effects they obtain a more decent agreement with the experimental liquid structure. Interestingly, and to some extent unexpectedly (at least to me), simple polarization effects seem to be doing most of the job. However, in the end they have just another rigid water potential which, judging from their pair correlation function, appears to be almost as good as the best existing rigid empirical potentials for water. Of course conceptually it is not the same thing: the potential of Bukowski et al. is not tuned to reproduce experiment but is derived from accurate quantum mechanical calculations on the dimer (and in the most accurate case also on the trimer). This is an important achievement but difficult to generalize to a wide range of possible contexts, including e.g. flexible molecules, solvation effects, hydrophobic and hydrophilic conditions, confinement, systems other than water, etc. Keeping the attention on the liquid structure, I do not think that this potential gives us a better insight than what we already know on the structure of the H bond network in water.
I do not think that the work of Bukowski et al. goes anywhere beyond ab-initio molecular dynamics on the flight. [This is the kind of approach used by Dominik Marx, Roberto Car and others.] The latter models water in terms of nuclei and electrons, the former in terms of rigid intermolecular potentials. The latter produces flexible molecules by construction and describes in detail the interplay between electronic structure and nuclear dynamics. As such ab-initio molecular dynamics on the flight is applicable to the most general range of situations, including for instance proton transfer effects (the Grotthus mechanism) that are not allowed by any rigid classical intermolecular potential (whether ab-initio or empirical!). The main limitation of ab-initio molecular dynamics on the flight, apart from numerical cost, is due to the limited accuracy of existing approximations of density functional theory, but as these improve or if highly efficient and more accurate electronic structure methods are established, immediately this progress could be transferred to the accuracy of ab-initio molecular dynamics on the flight. It is not so with the method of Bukowski et al. which cannot go beyond the accuracy provided by a simple polarization approximation of the non additive many-body effects. The only system that Bukowski et al. describe better than ab-initio MD on the flight is the potential of the water dimer, for which they rely on the most accurate available quantum chemical methods. Some of this accuracy is lost when they go to condensed phase. Judging from their current results, they do not have a better description of liquid water than that provided by existing empirical potential. Overall, the insight on a number of properties, structural, electronic etc. provided by ab-initio MD (on the flight) is vastly superior.”
That squares with my second adviser, who says:
“I'm very disappointed such a paper is appearing in Science. First, it gives an incredibly misleading account of the existing ab-initio literature on water, almost entirely based on DFT, and it omits a very recent quantum chemistry work (by S. Xantheas, JCP). Second, it "sells" a quantum chemistry based ab-initio potential as giving excellent agreement with experiment, when the agreement the authors find is not much better than what is already out there, in my opinion, using DFT based methods. This Science work just has different kind of disagreement with experiments, w.r.t. those found with DFT, but overall it is not much better (actually it is a little worse); most importantly: this work does NOT solve any of the open problems out there on the structure and/or properties of water.
Let me elaborate a bit on these points.
1) Account of the existing ab-initio literature. In the last ~ 15 years, most of the ab-initio simulations (no fit to or input from experiment) of water have been performed using Density Functional Theory (DFT) based methods, and they have been carried out mainly with two functionals (so call gradient corrected energy functionals, BLYP and PBE). Both of these functionals have been believed to give good agreement with experiment for several years (specifically until ~ 2004). However the good agreement with experiment did NOT came from a good performance of the theory but was somehow fortuitous, due to numerical inaccuracies in the solution of the Kohn-Sham equations (DFT equations) for the electrons. In 2004, Schwegler et al. (JCP 2004) and Grossman et al. (JCP 2004) pointed out these inaccuracies and showed that water correlation functions (g(r)--also discussed in the Science report) are over-structured and diffusion too slow (w.r.t. to experiment), when the numerics is done right and numerical inaccuracies are removed. The fortuitous agreement with experiment there had nothing to do with a fortuitous choice of functionals, as stated by the present authors. It had to do with numerics adopted in integration techniques. Those results were later confirmed by Sit and Marzari (JCP 2005), Serra and Artacho (several papers appeared in JCP) and others (none of these papers are mentioned in the Science report). So the introduction of this Science paper carelessly dismisses an approach (ab-initio simulations based on DFT) that, although not in full quantitative agreement with experiment, can describe well the salient, qualitative features of liquid water. On top of this, at the end of the paper the reader realizes that such a dismissed approach gives an agreement with experiment which is similar to the one found by the authors (Actually, I'd like to claim that DFT gives a more consistent agreement with experiment than the one presented in the Science paper--see below).
I also note that in the introduction of this paper dispersion forces in water are declared "non negligible", with no reference and no discussion. There is no clear experimental and/or theoretical evidence to support this statement. It may well be so, but nobody knows right now.
2) Agreement with experiment found here, wrt to existing ab-initio simulations based on DFT. I believe that overall the g(r)s found here are not in better agreement with experiment than those described by DFT ( see original papers). Note that they get right the first peak of g(r) but not the second, implying that their angular distribution (if they had computed it!) is quite inaccurate. Their diffusion coefficient is in agreement with experiment, but admittedly (see their own statement), this agreement may be fortuitous as they neglected proton quantum effects of the monomer, and flexibility of the monomer. Most importantly: the coordination number they find is an overestimate of what is accepted in the field as a reasonable number extracted from experiment. They find 5.6 and they compared it with 4.8, when the accepted value in the field is more like 4.3/4.5. Whatever number you consider, their value is a big overestimate, giving a liquid over-coordinated with respect to experiment and with too many hydrogen bonds. Their computed internal energies are worst that those obtained with empirical potentials.”
Not encouraging, then. And in any event, in terms of understanding hydration, the messages were as follows:
Adviser 1: “I cannot predict what would be the result of applying the scheme of Bukowski et al. to studying hydrophobic hydration. In hydrophobic hydration there is not just water but also the hydrophobic substance (unless this is the vacuum). Since their potential should be more transferable than standard empirical potentials, it should describe better the water close to a hydrophobic solute. In principle, however, all the interactions need to be included and treated with comparable accuracy. In our study of solvated methane both the water and the methane molecules have deformable electronic clouds that play an important role in the outcome of the calculation. These effects are described to some extent in terms of polarization effects by Bukowski et al., but they should also construct a methane-methane and a methane-water potential and include polarization effects beyond pure water in order to tackle the solvation problem.”
Adviser 2: “Even if the authors found the best ab-initio potential for liquid water fitted to ab-initio gas phase data, if they wanted to describe solvation they would have to start all over again, as they would have to do the fitting to gas phase water containing the solvated molecule or they would have to add other pieces to the potential. This is why fitted potential (whatever they are fitted to: ab-initio data, experiment, etc.) will always have serious drawbacks. If you want to study systems that have not been fitted, containing other entities, well... here you go, you have to start again with your fit.”
So there it is. I don’t like dumping on a paper, but this one has come out in a very high-profile journal where it will get a lot of recognition. So I thought it is only right to put the record straight.
First, the general context, which is nicely explained by one of my advisers:
“Modeling intermolecular interactions in water in terms of simple pair potentials is difficult because these potentials miss the cooperative effects of the H bonds. These are important: it is because of these effects that the “dipole” moment of a molecule in condensed phase is greatly enhanced compared to that of a molecule in gas phase (here I put dipole in quotes because this quantity cannot be unambiguously defined in condensed phase). Compared to simple liquids water has a much more open structure dictated by the underlying H-bond network.
Popular empirical potentials for water try to model the basic physics underlying H bonds interactions, which are largely of electrostatic origin in terms of interactions between point charges spatially located in a way that mimics, albeit very approximately, the charge distribution in the water molecule. These potentials are tuned to reproduce a number of properties of the liquid but have limited transferability to different environments as they miss the cooperative effects mentioned above.”
Existing classical potentials for water describe well the structure of the liquid (even very well I would say!). The difficulty with this approach, apart from the issue of transferability, is that we also need to describe the interaction between water and other molecules, material systems etc. in order to model solvation, interfacial water, confined water etc. In addition, there are properties such as the dielectric properties that depend on the actual electronic structure of the molecules in a given environment. Furthermore water molecules can dissociate. A good potential for the intermolecular interaction does not tell us everything.”
But in terms of the Science paper by Bukowski et al., he says:
“I am not too impressed by the contribution of Bukowski et al. They show that an extremely good potential for the dimer (a pair potential) is not sufficient for a good description of the liquid, as one could have expected. Including non-additive many body effects they obtain a more decent agreement with the experimental liquid structure. Interestingly, and to some extent unexpectedly (at least to me), simple polarization effects seem to be doing most of the job. However, in the end they have just another rigid water potential which, judging from their pair correlation function, appears to be almost as good as the best existing rigid empirical potentials for water. Of course conceptually it is not the same thing: the potential of Bukowski et al. is not tuned to reproduce experiment but is derived from accurate quantum mechanical calculations on the dimer (and in the most accurate case also on the trimer). This is an important achievement but difficult to generalize to a wide range of possible contexts, including e.g. flexible molecules, solvation effects, hydrophobic and hydrophilic conditions, confinement, systems other than water, etc. Keeping the attention on the liquid structure, I do not think that this potential gives us a better insight than what we already know on the structure of the H bond network in water.
I do not think that the work of Bukowski et al. goes anywhere beyond ab-initio molecular dynamics on the flight. [This is the kind of approach used by Dominik Marx, Roberto Car and others.] The latter models water in terms of nuclei and electrons, the former in terms of rigid intermolecular potentials. The latter produces flexible molecules by construction and describes in detail the interplay between electronic structure and nuclear dynamics. As such ab-initio molecular dynamics on the flight is applicable to the most general range of situations, including for instance proton transfer effects (the Grotthus mechanism) that are not allowed by any rigid classical intermolecular potential (whether ab-initio or empirical!). The main limitation of ab-initio molecular dynamics on the flight, apart from numerical cost, is due to the limited accuracy of existing approximations of density functional theory, but as these improve or if highly efficient and more accurate electronic structure methods are established, immediately this progress could be transferred to the accuracy of ab-initio molecular dynamics on the flight. It is not so with the method of Bukowski et al. which cannot go beyond the accuracy provided by a simple polarization approximation of the non additive many-body effects. The only system that Bukowski et al. describe better than ab-initio MD on the flight is the potential of the water dimer, for which they rely on the most accurate available quantum chemical methods. Some of this accuracy is lost when they go to condensed phase. Judging from their current results, they do not have a better description of liquid water than that provided by existing empirical potential. Overall, the insight on a number of properties, structural, electronic etc. provided by ab-initio MD (on the flight) is vastly superior.”
That squares with my second adviser, who says:
“I'm very disappointed such a paper is appearing in Science. First, it gives an incredibly misleading account of the existing ab-initio literature on water, almost entirely based on DFT, and it omits a very recent quantum chemistry work (by S. Xantheas, JCP). Second, it "sells" a quantum chemistry based ab-initio potential as giving excellent agreement with experiment, when the agreement the authors find is not much better than what is already out there, in my opinion, using DFT based methods. This Science work just has different kind of disagreement with experiments, w.r.t. those found with DFT, but overall it is not much better (actually it is a little worse); most importantly: this work does NOT solve any of the open problems out there on the structure and/or properties of water.
Let me elaborate a bit on these points.
1) Account of the existing ab-initio literature. In the last ~ 15 years, most of the ab-initio simulations (no fit to or input from experiment) of water have been performed using Density Functional Theory (DFT) based methods, and they have been carried out mainly with two functionals (so call gradient corrected energy functionals, BLYP and PBE). Both of these functionals have been believed to give good agreement with experiment for several years (specifically until ~ 2004). However the good agreement with experiment did NOT came from a good performance of the theory but was somehow fortuitous, due to numerical inaccuracies in the solution of the Kohn-Sham equations (DFT equations) for the electrons. In 2004, Schwegler et al. (JCP 2004) and Grossman et al. (JCP 2004) pointed out these inaccuracies and showed that water correlation functions (g(r)--also discussed in the Science report) are over-structured and diffusion too slow (w.r.t. to experiment), when the numerics is done right and numerical inaccuracies are removed. The fortuitous agreement with experiment there had nothing to do with a fortuitous choice of functionals, as stated by the present authors. It had to do with numerics adopted in integration techniques. Those results were later confirmed by Sit and Marzari (JCP 2005), Serra and Artacho (several papers appeared in JCP) and others (none of these papers are mentioned in the Science report). So the introduction of this Science paper carelessly dismisses an approach (ab-initio simulations based on DFT) that, although not in full quantitative agreement with experiment, can describe well the salient, qualitative features of liquid water. On top of this, at the end of the paper the reader realizes that such a dismissed approach gives an agreement with experiment which is similar to the one found by the authors (Actually, I'd like to claim that DFT gives a more consistent agreement with experiment than the one presented in the Science paper--see below).
I also note that in the introduction of this paper dispersion forces in water are declared "non negligible", with no reference and no discussion. There is no clear experimental and/or theoretical evidence to support this statement. It may well be so, but nobody knows right now.
2) Agreement with experiment found here, wrt to existing ab-initio simulations based on DFT. I believe that overall the g(r)s found here are not in better agreement with experiment than those described by DFT ( see original papers). Note that they get right the first peak of g(r) but not the second, implying that their angular distribution (if they had computed it!) is quite inaccurate. Their diffusion coefficient is in agreement with experiment, but admittedly (see their own statement), this agreement may be fortuitous as they neglected proton quantum effects of the monomer, and flexibility of the monomer. Most importantly: the coordination number they find is an overestimate of what is accepted in the field as a reasonable number extracted from experiment. They find 5.6 and they compared it with 4.8, when the accepted value in the field is more like 4.3/4.5. Whatever number you consider, their value is a big overestimate, giving a liquid over-coordinated with respect to experiment and with too many hydrogen bonds. Their computed internal energies are worst that those obtained with empirical potentials.”
Not encouraging, then. And in any event, in terms of understanding hydration, the messages were as follows:
Adviser 1: “I cannot predict what would be the result of applying the scheme of Bukowski et al. to studying hydrophobic hydration. In hydrophobic hydration there is not just water but also the hydrophobic substance (unless this is the vacuum). Since their potential should be more transferable than standard empirical potentials, it should describe better the water close to a hydrophobic solute. In principle, however, all the interactions need to be included and treated with comparable accuracy. In our study of solvated methane both the water and the methane molecules have deformable electronic clouds that play an important role in the outcome of the calculation. These effects are described to some extent in terms of polarization effects by Bukowski et al., but they should also construct a methane-methane and a methane-water potential and include polarization effects beyond pure water in order to tackle the solvation problem.”
Adviser 2: “Even if the authors found the best ab-initio potential for liquid water fitted to ab-initio gas phase data, if they wanted to describe solvation they would have to start all over again, as they would have to do the fitting to gas phase water containing the solvated molecule or they would have to add other pieces to the potential. This is why fitted potential (whatever they are fitted to: ab-initio data, experiment, etc.) will always have serious drawbacks. If you want to study systems that have not been fitted, containing other entities, well... here you go, you have to start again with your fit.”
So there it is. I don’t like dumping on a paper, but this one has come out in a very high-profile journal where it will get a lot of recognition. So I thought it is only right to put the record straight.
Tuesday, February 27, 2007
Dynamics fast and slow
How do motions of a peptide chain depend on those of water molecules in the hydration shell? It's clear enough that the solvent dynamics can span a wide range of timescales, depending on how the water molecules interact with the protein. But in unfolded and molten-globule states, it has been suggested that there are rapid fluctuations between various helical and beta-sheet-like states that are 'lubricated' by picosecond rearrangements of the hydrogen-bonding network in water. This lubrication, enabled by the rapidity of H-bond making and breaking, is presumed to enable protein folding. Neil Hunt and coworkers at the University of Strathclyde now say that they've seen such motions, with timescales of a few tenths of a picosecond, in the alpha-helix-to-random-coil transition of a homo-polypeptide (poly-L-lysine) using optical Kerr-effect spectroscopy. The paper is here.
Motions in the hydration shell that are one or two orders of magnitude slower have been studied by Dongping Zhong and colleagues at Ohio State, in a paper here. They're looking at apomyoglobin, and specifically at the water dynamics around a tryptophan group (Trp7), which serves as a convenient fluorescent chromophore. Through both experiment and simulation, they find slow relaxation on timescales of 5-87 ps. Previous studies of such slow dynamics have offered divergent interpretations. Ahmed Zewail and coworkers have suggested that these dynamics are due to the effect of the protein's potential field on the hydration water (J. Phys. Chem. B 107, 13218; 2003). Bertil Halle thinks that these water dynamics close to a protein aren't so different from those in the bulk (PNAS 102, 13867; 2005). Zhong and colleagues say that the slow relaxation is due to strongly coupled water-protein motions. If either the water or the protein is frozen in the simulations, the slow component disappears. I guess that supports the contention of Bizzarri and Cannistraro that the dynamics of the protein and solvent are so strongly coupled that they ‘should be conceived as a single entity' (J. Phys. Chem. B 106, 6617-6633; 2002).
Motions in the hydration shell that are one or two orders of magnitude slower have been studied by Dongping Zhong and colleagues at Ohio State, in a paper here. They're looking at apomyoglobin, and specifically at the water dynamics around a tryptophan group (Trp7), which serves as a convenient fluorescent chromophore. Through both experiment and simulation, they find slow relaxation on timescales of 5-87 ps. Previous studies of such slow dynamics have offered divergent interpretations. Ahmed Zewail and coworkers have suggested that these dynamics are due to the effect of the protein's potential field on the hydration water (J. Phys. Chem. B 107, 13218; 2003). Bertil Halle thinks that these water dynamics close to a protein aren't so different from those in the bulk (PNAS 102, 13867; 2005). Zhong and colleagues say that the slow relaxation is due to strongly coupled water-protein motions. If either the water or the protein is frozen in the simulations, the slow component disappears. I guess that supports the contention of Bizzarri and Cannistraro that the dynamics of the protein and solvent are so strongly coupled that they ‘should be conceived as a single entity' (J. Phys. Chem. B 106, 6617-6633; 2002).
Friday, February 16, 2007
Quantum hydrophobicity, and water in drug design
A paper in PNAS by Ned Wingreen at Princeton uses first-principles quantum molecular dynamics to look at the hydrophobic interaction of two hydrated methane molecules. I’d say this is primarily a methodological paper – that’s pretty much the angle Guilia Galli takes in her accompanying commentary – in that it aims to establish how well classical simulation approaches do in capturing the nature of the interaction. The answer is: not particularly. Classical force fields predict two free-energy minima, one at methane-methane contact and a second, rather shallow, at a small separation corresponding to one intervening solvent layer. But the relative stability of these minima is rather sensitive to the precise force-field parameters and can be reversed for some values. The quantum simulations reveal a deeper potential well at contact – this is always the stable configuration. But there is a succession of shallow minima at several other methane-methane separations, implying the existence of several relatively stable hydration cages. The authors talk about these configurations in terms of ‘clathrate-like cages’, but in fact it seems that they don’t have any well-defined hydrogen-bonding arrangements: two independent simulations at a specific methane-methane separation gave water structures that could not be superimposed. (There was, however, apparently some consistency in the ‘hydrogen-bonded rings’ in between the two methanes.) Wingreen and colleagues suggest that the shallow minima are the result of relatively well packed configurations for water in the hydration shells. But I don’t really know what this means. Normally, considerations of molecular packing in liquids are governed by the short-ranged repulsion between molecules. But ‘well packed’ is an ambiguous term for water, where optimal hydrogen bonding means that the waters sit rather further apart than equivalent spherical molecules would do. I’m assuming ‘well packed’ here refers to unbroken, unstrained H-bonding configurations…?
**********
“Until a few years ago it was common practice to ignore water molecules in protein binding sites”, say Jonathan Essex and colleagues at Southampton University in an interesting paper in JACS. But now, they point out, there is increasing interest in designing ligands that will displace particular water molecules in drug binding. Conceivably, this might make the ligands more selective and the binding energy more favourable.
But it’s not clear whether that will necessarily be so. Despite the entropic advantage of expelling bound water from a binding cleft, one can’t generalize about the consequent free energy change. Whether or not it is advantageous to incorporate a water molecule at the binding interface hinges on a delicate balance. Confining a water molecule clearly has an entropic penalty, but this might be repaid by the enthalpic gains of hydrogen-bond formation – an issue that must itself be weighed against the average number of hydrogen bonds that a bulk water molecule engages in. Jack Dunitz (Science 264, 670 (1994); Chem. Biol. 2, 709-712 (1995)) has estimated that transferring a water molecule from an ordered binding site where it is bound by an ‘average’ hydrogen bond to the bulk involves an overall free-energy change that is close to zero. So it is not obvious which way the scales will tip in any instance. John Ladbury and his coworkers (D. A. Renzoni, M. J. J. M. Zvelebil, T. Lundbäck & J. E. Ladbury, in J. E. Ladbury & P. R. Connelly (eds), Structure-Based Drug Design: Thermodynamics, Modeling and Strategy 161-180, Landes Bioscience (1997)) have thought about the implications for drug design.
The message is illustrated in the binding of various inhibitors of HIV-1 protease, one of the key targets in AIDS therapies. Crystal structures show that some of these, such as KNI-272, bind to the enzyme via a bridging water molecule. Other inhibitors, such as DMP450, have been designed specifically to exclude this water molecule while mimicking its hydrogen-bonding capacity, and have found to bind more strongly. Li and Lazaridis (JACS 125, 6636-6637 (2003)) have calculated that displacement of the bound water by DMP450 is in itself unfavourable relative to KNI-272, but that this cost is outweighed by the lower cost of desolvating DMP450 to form the bound complex. So the consequences of eliminating the water molecule are both highly specific and not obvious.
With all this in mind, Essex and colleagues have sought a way of classifying water molecules in protein binding sites according to how easily displaced they are by ligands. By studying the thermodynamics of six proteins complexed with a variety of ligands, they say that the water molecules in the binding sites seem to come in two classes: those that are readily displaced (by at least some ligands), and those that never are. The latter, unsurprisingly, turn out to be more tightly bound according to MC simulations. All the same, the authors say that “no linear correlation exists between the binding free energies of water molecules and the change in binding affinity of ligands displacing the water molecules.” Yet they conclude that if we can identify the ‘conserved’ water molecules – those that do not get displaced come what may – then these can be usefully used in the design of drug docking: in effect, they serve as ‘part of the protein’, available for hydrogen bonding to the ligand. This paper supplies some heuristics for deciding which water molecules are conserved or displacable.
**********
“Until a few years ago it was common practice to ignore water molecules in protein binding sites”, say Jonathan Essex and colleagues at Southampton University in an interesting paper in JACS. But now, they point out, there is increasing interest in designing ligands that will displace particular water molecules in drug binding. Conceivably, this might make the ligands more selective and the binding energy more favourable.
But it’s not clear whether that will necessarily be so. Despite the entropic advantage of expelling bound water from a binding cleft, one can’t generalize about the consequent free energy change. Whether or not it is advantageous to incorporate a water molecule at the binding interface hinges on a delicate balance. Confining a water molecule clearly has an entropic penalty, but this might be repaid by the enthalpic gains of hydrogen-bond formation – an issue that must itself be weighed against the average number of hydrogen bonds that a bulk water molecule engages in. Jack Dunitz (Science 264, 670 (1994); Chem. Biol. 2, 709-712 (1995)) has estimated that transferring a water molecule from an ordered binding site where it is bound by an ‘average’ hydrogen bond to the bulk involves an overall free-energy change that is close to zero. So it is not obvious which way the scales will tip in any instance. John Ladbury and his coworkers (D. A. Renzoni, M. J. J. M. Zvelebil, T. Lundbäck & J. E. Ladbury, in J. E. Ladbury & P. R. Connelly (eds), Structure-Based Drug Design: Thermodynamics, Modeling and Strategy 161-180, Landes Bioscience (1997)) have thought about the implications for drug design.
The message is illustrated in the binding of various inhibitors of HIV-1 protease, one of the key targets in AIDS therapies. Crystal structures show that some of these, such as KNI-272, bind to the enzyme via a bridging water molecule. Other inhibitors, such as DMP450, have been designed specifically to exclude this water molecule while mimicking its hydrogen-bonding capacity, and have found to bind more strongly. Li and Lazaridis (JACS 125, 6636-6637 (2003)) have calculated that displacement of the bound water by DMP450 is in itself unfavourable relative to KNI-272, but that this cost is outweighed by the lower cost of desolvating DMP450 to form the bound complex. So the consequences of eliminating the water molecule are both highly specific and not obvious.
With all this in mind, Essex and colleagues have sought a way of classifying water molecules in protein binding sites according to how easily displaced they are by ligands. By studying the thermodynamics of six proteins complexed with a variety of ligands, they say that the water molecules in the binding sites seem to come in two classes: those that are readily displaced (by at least some ligands), and those that never are. The latter, unsurprisingly, turn out to be more tightly bound according to MC simulations. All the same, the authors say that “no linear correlation exists between the binding free energies of water molecules and the change in binding affinity of ligands displacing the water molecules.” Yet they conclude that if we can identify the ‘conserved’ water molecules – those that do not get displaced come what may – then these can be usefully used in the design of drug docking: in effect, they serve as ‘part of the protein’, available for hydrogen bonding to the ligand. This paper supplies some heuristics for deciding which water molecules are conserved or displacable.
Thursday, February 8, 2007
Why does water do that?
Following on from the suggestion that the ‘glass transition’ of hydrated proteins at around 220 K is in fact a change in the dynamical state of the water of hydration, related to the second critical point of water in the supercooled, high-pressure state (see the previous post below), Gene Stanley and his coworkers have now released a preprint that aims to relate this behaviour to the molecular structure of the liquid. To recap, the idea is that the ‘transition’ corresponds to the crossing of the so-called Widom line, the locus of the maximum in the correlation length that extends like a ‘ghost’ of the liquid-liquid phase transition beyond the critical point at which this transition vanishes. (At the critical point itself, this correlation length diverges.) What does this mean for the nature of the hydrogen-bonded network? Well, it’s subtle. In thermodynamic terms, the crossover corresponds to a change from non-Arrhenius dynamics at high temperature (the activation energy for water diffusion depends on temperature) to Arrhenius dynamics at low temperature (temperature-independent activation energy). The simulations by Stanley and colleagues now suggest that this crossover shows up in terms of the tendency to form non-bifurcated hydrogen bonds. Bifurcated H-bonds are fairly common in the liquid state at ambient conditions, and seem to be the defects that make diffusion facile. The Widom line corresponds to the point at which the derivative of the probability of forming non-bifurcated bonds with respect to temperature is maximal. Hmm, not an easy quantity to visualize. But it does mean that above the Widom line, supercooled water has fewer non-bifurcated hydrogen bonds, and so is less ‘tetrahedral’ and denser (like the high-density liquid phase), than below this line. That’s intuitive enough – crossing the Widom line as temperature decreases corresponds to the supercooled liquid becoming less dense, more structured and changing from a high-density-liquid-like to a low-density-liquid-like state: a ghost of the liquid-liquid phase transition itself, but with no abrupt change of thermodynamic variables.
Lawrence Pratt at Los Alamos and his coworkers have posted a preprint entitled ‘What is special about water as a matrix of life?’. This, as I recall, is basically the paper that Lawrence presented at the Varenna meeting, convened to discuss that very question in early 2005. The paper aims to address the title question in general, and hydrophobic effects in particular, by moving away from an emphasis on structure (for example, the classic Kauzmann model of entropically driven hydrophobic attraction due to the release of ‘structured water’ in the space between hydrophobes) and focusing instead on what the authors call the ‘engineering characteristics of the liquid’, such as its equation of state. The key message seems to be that, as a solvent for life, water represents a safe bet: the liquid state exists over a wide temperature range (compared with other simple small-molecule solvents), and within that range there is rather little variation in thermodynamic variables: response functions such as the compressibility and thermal expansion coefficient, as well as the nature of the solvophobic effect, vary little. It seems to me that this raises several questions (which are not objections), such as: does life require a wide temperature range, or does it just fill up whatever niches are available? (Thermophiles seem very ancient; would life have happily persisted if all Earth’s water was warm?) Does the ‘structured’ character of water play any necessary role in this scheme? (It does seem to be biologically important that water forms directional H-bonds.) And what does underlie the hydrophobic attraction, and how general is it?
Speaking of which, more evidence for the role of nanobubbles in the long-range hydrophobic attraction is provided in a study of nanoparticle adsorption onto a quartz microbalance in gassed and degassed water by Sangmin Jeon and colleagues in Korea (Langmuir 23, 1623-1625 (2007)).
Lawrence Pratt at Los Alamos and his coworkers have posted a preprint entitled ‘What is special about water as a matrix of life?’. This, as I recall, is basically the paper that Lawrence presented at the Varenna meeting, convened to discuss that very question in early 2005. The paper aims to address the title question in general, and hydrophobic effects in particular, by moving away from an emphasis on structure (for example, the classic Kauzmann model of entropically driven hydrophobic attraction due to the release of ‘structured water’ in the space between hydrophobes) and focusing instead on what the authors call the ‘engineering characteristics of the liquid’, such as its equation of state. The key message seems to be that, as a solvent for life, water represents a safe bet: the liquid state exists over a wide temperature range (compared with other simple small-molecule solvents), and within that range there is rather little variation in thermodynamic variables: response functions such as the compressibility and thermal expansion coefficient, as well as the nature of the solvophobic effect, vary little. It seems to me that this raises several questions (which are not objections), such as: does life require a wide temperature range, or does it just fill up whatever niches are available? (Thermophiles seem very ancient; would life have happily persisted if all Earth’s water was warm?) Does the ‘structured’ character of water play any necessary role in this scheme? (It does seem to be biologically important that water forms directional H-bonds.) And what does underlie the hydrophobic attraction, and how general is it?
Speaking of which, more evidence for the role of nanobubbles in the long-range hydrophobic attraction is provided in a study of nanoparticle adsorption onto a quartz microbalance in gassed and degassed water by Sangmin Jeon and colleagues in Korea (Langmuir 23, 1623-1625 (2007)).
Subscribe to:
Posts (Atom)