Skip to main content

Drug Assays

Down Amongst the Water Molecules

This new paper shows one reason why it’s so tricky to calculate compound binding in an active site (which is what we’d want to do in order to do effective in silico virtual screening). The authors (a multicenter team from York, Demuris, Vernalis, and St. Jude) are looking at a virulence protein from H. influenzae called SiaP. It binds sialic acid, but it also binds a whole set of water molecules that can be seen in the X-ray data. These are arranged into a hydrogen-bonded network that bridges that protein surface and the ligand, and while this is an unusually extensive example, this sort of thing is quite common.

Even single water molecules can have a big influence, and to make things more amusing, they can show up with either an enthalpic or entropic effect (or both, why not) and with either positive or negative influence on the overall free energy of binding. You’d think that something the size of a water molecule would be only a pawn-level piece on the binding-pocket chess board, and sometimes that’s all they are. But they can suddenly turn into rooks and queens depending on the situation (not a bad analogy, considering that their hydrogen bonds are so directional) and they can at the same time turn into either white pieces or black ones. Which situation would make it rather hard to write a useful chess-playing program, and it makes it hard to write a useful binding-prediction one, too.

At right are two X-ray structures of N-acetylneuraminic acid bound to SiaP – the left one is the native protein, and you can count ten water molecules involved in the binding. On the right, though is a point mutant (A11N, alanine to asparagine, down at the bottom of the figure). Most of the water molecules are still present, and in the same place. But the network around water #3 has been disrupted. And now is the time to mention that the Kd of Neu5Ac with the wild-type SiaP is 30 nM, while the Kd with the mutant form is about 42,000 nM. You lose 1400x binding affinity with that one distortion of the water hydrogen-bonding network. Otherwise, the positions of the protein residues and of the ligand itself are within experimental error between the two structures.

This brings to mind another effect that will be familiar to fragment-based drug discovery types. If you look at a series of X-ray structures of different fragments against the same target and ask someone to rank-order their binding affinities just from those, they can’t do it. In some cases you might be able to pick out the best one, if it makes a notably different set of hydrogen bonds or other interactions, but if they’re all roughly the same binding mode you can generally forget it. The differences are too subtle, and they’re too subtly distributed between entropic and enthalpic effects on the ligand, the associated water molecules, and the protein structure itself.

If you’d shown me those two structures above with no other information, I probably would have guessed that the one on the left would be better, because most of the time another hydrogen bonded water interaction is a good thing. But not always: that bound water #4 very likely helped the entropic term when it left its structured state, for example, and sometimes that wins out over the enthalpic effects. And I certainly wouldn’t have guessed a 1400x drop in affinity!

One interesting aspect of this paper is that the crystal structures of both species were collected at cryogenic temperature and at room temperature, and the latter revealed more disorder in the mutant Asn residue on warming. That’s another factor that not everyone appreciates when they look at an X-ray-structure: you’re generally seeing the lowest frozen well of the energy landscape, and that may or may not be as informative as you’d want about the situation at RT (which is where you’re probably running your binding assays!) We have a hard enough time understanding static structures, and adding in dynamic effects just makes things even messier, but things do tend to move around out here in the real world. The movement in the Asn side chain complicated the authors’ attempts to work out a thermodynamic model, and this is actually a pretty well-behaved system compared to many others that you’ll see.

At any rate, the more detailed thermodynamic study laid out in the paper indicates that the majority of the energy change is indeed due to disruption of the water network, and there’s a strong contribution of the enthalpy term. As the paper shows, if you go through the literature looking at the estimates of enthalpic and entropic energy changes for desolvation of a single water molecule, you get numbers all over the place, and the differences between them are all of a magnitude to show up very strongly as differences in binding affinity. It looks like this case agrees more with the high side of those estimates, but anyone who’s in love with virtual screening and computational drug design should consider the facts that (A) these values can be so high and so meaningful to binding, and (B) that the best experts in the field are still disagreeing about what those numbers are. It’s quite likely that trying to assign a single average value to such things isn’t too useful anyway, given the wide variety of water molecule environments, but that’s not helping to simplify anyone’s life, is it?

17 comments on “Down Amongst the Water Molecules”

  1. Peter Kenny says:

    One question would be whether the changes to the water network are a cause of the affinity reduction or merely a response to to the introduction of a side chain that is larger and higher maintenance (in the context of hydrogen bonding) than alanine. There could also have implications for how we interpret about alanine scanning results. In case it’s of interest, I’ve linked a blog post on hydrogen bonding asymmetry as the URL for this comment.

  2. John Wayne says:

    I’d like to give the authors credit for doing something hard; keep up the good work.

  3. Yajun Zheng says:

    “The movement in the Asp side chain complicated the authors’ attempts to work out a thermodynamic model, and this is actually a pretty well-behaved system compared to many others that you’ll see.” Here Asp should be Asn.

  4. luysii says:

    Old friend from grad school days Don Voet in his Biochemistry text notes “Water is a chemically reactive liquid with such extraordinary physical properties that, if chemists had discovered it in recent times, it would undoubtedly been classified as an exotic substance.”

    1. metaphysician says:

      Silly question time: do any of the other chalcogen hydrides ( H2S, H2Se, etc ) show similarly exotic properties to H2O, at least when they are at suitable temperature to be liquids?

      1. DRP says:

        Not really, as far as I am aware, as the interaction between the chalcogen and hydrogen isn’t very strong (too much of a miss match in orbital energies) – hence H2S is a gas despite being a much heavier molecule.

      2. TWS says:

        Yeah, I discovered H2S has very different properties to H2O during my undergrad years. Very much lower boiling / melting points, much smaller range over which the gas is a liquid and (only really of relevance to the project I was working on) an infrared absorption cross-section about the size of a small pea compared to water which is of course an extremely strong absorber of IR light. During our very early work, this did lead to us opening up a gas cell to check if the process was indeed evolving H2S, as we couldn’t see a thing in the spectrum. Let’s just say it didn’t take us long to realise the fault was with the poor sensitivity of IR spectroscopy towards that particular gas.

        Later on we used a multi-pass gas cell with silver (!) mirrors in direct contact with the gas (this time we knew the risk, but it’s all we had). They lasted a surprisingly long time (about a year and a half with more or less continuous use) due to the protective coating, but eventually the H2S found a way through, blackening the edges of the mirrors. Switched over to gold mirrors from then on which have near enough identical performance in the wavelength ranges we were interested in.

    2. eub says:

      We’d call it a self-assembling mesh of nanobots except really ångströbots.

  5. Julien says:

    There is a bunch of molecular modelling methodologies that could have a go at this sort of problem in principle (FEP, GCMC, 3D-RISM, IFST, cell theory …to name just a few). Will be interesting to see which of these methodologies can pick up the effect on the ligand binding energetics.

  6. PR says:

    It’s all focused on binding of the ligand in the active site, but isn’t it possible that the shift in Kd can be explained by more stable inactive binding site in the case of mutated enzyme? I.e. there could even be some stabilizing water molecule bound to Asn which makes the ligand-free pocket more compact and less prone to opening and accepting ligand. No?

    1. XJ says:

      Excellent point. It’s the difference between the free energies of the apo state and inhibitor-bound state that determines the binding affinity in each individual cases. WT and mutant proteins have different apo state free energies. Only comparing the bound state of the two proteins can therefore over or underestimate the differences in binding affinities. Water molecules only make things more complicated on top of that.

      1. PR says:

        Exactly. I didn’t read the paywalled paper, but I would love to see comparison of structures of the two apoenzymes, WT and A11N. I could imagine the amide from Asn reaching somewhere accross the ligand-free binding pocket and thus stabilizing it. Water or not.

  7. KN says:

    Water is a perfect fragment, as they say in FBDD.

  8. eub says:

    If I remember my P-Chem, that’s a difference of ln(1500) * 273 K * 8.3145 J/mol/K = 16.6 kJ/mol. And the more-structured water variant has to be winning on enthalpy because it’s losing on entropy, right?

    I see a loss of four hydrogen bonds, which implies 4 kJ/mol for a hydrogen bond, which is out the high side of the normal fuzzy range for it. Were these super-happy hydrogen bonds or is there something else going on?

  9. Tc says:

    Wasn’t it just a week ago that structure didn’t matter to Derek and now we are talking about water networks?

    1. Derek Lowe says:

      Hah! It’s not that structure *can’t* matter, though – in fact eventually it’s going to have to. And for that to happen, we need to understand this sort of thing.

Leave a Reply

Your email address will not be published. Required fields are marked *

Time limit is exhausted. Please reload CAPTCHA.

This site uses Akismet to reduce spam. Learn how your comment data is processed.