Skip to content

Science1 publisher2 min readPublished

A network trained on crystal hydrogens disputes the tautomers of 126 PDBbind ligands

Tautomer assignment has been inference wherever protein crystallography cannot see hydrogen atoms. An NYU group borrowed a million experimental answers from small-molecule crystals and trained a network to make the call.

The Scientist · Science desk

Illustration accompanying A network trained on crystal hydrogens disputes the tautomers of 126 PDBbind ligands

What happened

  • Protein Data Bank structures are usually too low in resolution to locate hydrogen atoms, so the tautomeric state of a bound ligand is inferred by the scientists who model it.
  • Xiaolin Pan mined the Cambridge Structural Database, where high-resolution small-molecule crystals do show hydrogens, and assembled more than 1.1 million experimental tautomeric states.
  • The NYU group trained a graph neural network on that set to pick the stable tautomer straight from a 2D molecular structure, with no 3D coordinates and no quantum-mechanical calculation.
  • Run on 5,075 PDBbind ligands with more than one possible tautomer, the model judged 126 of the deposited assignments likely wrong and proposed alternatives with better hydrogen bonding.

Compiled by The ScientistSomething wrong?How this is made

Why it matters

  • decision Anyone using PDBbind ligands as ground truth for training or benchmarking a scoring function now chooses between re-curating a flagged subset and knowingly keeping chemistry the model disputes.
  • capability Quantum-mechanical treatment is too expensive to run across a full library; a network that reads 2D structures puts tautomer curation inside the preprocessing budget.
  • constraint The model's experimental grounding is hydrogen positions in crystals, so it cannot settle which tautomer dominates in water or inside a protein pocket.
  • precedent If depositors accept Zhang's framing, individual database entries become open to re-annotation requests on hydrogen placement alone, without touching the underlying protein structure.

The percentage needs its denominator. The 5,075 PDBbind ligands in the test were selected because each admits more than one tautomer, so the 126 disagreements sit inside the population where disagreement was possible at all [13]. That is about one in 40 of the ambiguous cases [17].

In those 126 cases the disagreement is about chemical plausibility. "Although this may seem like a small change, different tautomers of the same molecule can alter how a molecule interacts with a protein target," said Yingkai Zhang, the NYU chemistry professor and senior author of the study, which is published in Chemical Science [4][2]. The model moved a hydrogen and the resulting hydrogen-bonding pattern improved; in one ligand the reassignment let the molecule form additional hydrogen bonds with nearby protein residues [13][15]. "Reassigning these tautomers generally produced more chemically reasonable interactions. This does not mean that the experimentally determined protein structures themselves are incorrect; rather, our results suggest that the previously assigned chemical representation may warrant revision," Zhang said [14].

The training labels come out of crystals [11]. Molecules in the Cambridge Structural Database sit packed in a lattice; a drug candidate binds in a solvated pocket. The other source of experimental truth is solution-phase measurement, and those datasets hold a few hundred molecules [8]. Take 300 as a working size: the crystal-derived set is then roughly 3,700 times larger [18].

phys.org's account puts the throughput at 4.6 million compounds in hours [16]. The evidence offered for the corrections is the hydrogen-bonding pattern, and the account does not report a docking score or an affinity estimate recomputed before and after the reassignment [19]. For a screening campaign, the test worth running is whether any of the 126 corrections changes a rank order.

On the evidence reported, the sensible use is a bulk pass whose disagreements go into a review queue. Two conditions attach to that. The labels are solid-state, and the check on the reassignments is internal chemical reasonableness rather than an independent measurement of where the hydrogen sits in the bound complex [13]. Calling library curation finished would take a solution-phase test set far larger than the few hundred molecules now available [8], or a downstream result showing that corrected tautomers predict binding affinities better than the deposited ones.

What to watch

  • Whether the model is tested against solution-phase tautomer measurements, the data type it was not trained on.
  • Whether PDBbind or Protein Data Bank curators act on any of the 126 proposed ligand reassignments.
  • Whether re-docking or free-energy calculations on the corrected ligands predict known affinities better than the deposited assignments.
Loading claim ledger
Loading source directory links
Loading share composer
Loading topic controls
Loading related stories