Premium
XModeScore : a novel method for accurate protonation/tautomer‐state determination using quantum‐mechanically driven macromolecular X‐ray crystallographic refinement
Author(s) -
Borbulevych Oleg,
Martin Roger I.,
Tickle Ian J.,
Westerhoff Lance M.
Publication year - 2016
Publication title -
acta crystallographica section d
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 7.374
H-Index - 138
ISSN - 2059-7983
DOI - 10.1107/s2059798316002837
Subject(s) - tautomer , protonation , macromolecule , crystallography , quantum , state (computer science) , chemistry , computational chemistry , stereochemistry , physics , quantum mechanics , computer science , algorithm , organic chemistry , biochemistry , ion
Gaining an understanding of the protein–ligand complex structure along with the proper protonation and explicit solvent effects can be important in obtaining meaningful results in structure‐guided drug discovery and structure‐based drug discovery. Unfortunately, protonation and tautomerism are difficult to establish with conventional methods because of difficulties in the experimental detection of H atoms owing to the well known limitations of X‐ray crystallography. In the present work, it is demonstrated that semiempirical, quantum‐mechanics‐based macromolecular crystallographic refinement is sensitive to the choice of a protonation‐state/tautomer form of ligands and residues, and can therefore be used to explore potential states. A novel scoring method, called XModeScore , is described which enumerates the possible protomeric/tautomeric modes, refines each mode against X‐ray diffraction data with the semiempirical quantum‐mechanics (PM6) Hamiltonian and scores each mode using a combination of energetic strain (or ligand strain) and rigorous statistical analysis of the difference electron‐density distribution. It is shown that using XModeScore it is possible to consistently distinguish the correct bound protomeric/tautomeric modes based on routine X‐ray data, even at lower resolutions of around 3 Å. These X‐ray results are compared with the results obtained from much more expensive and laborious neutron diffraction studies for three different examples: tautomerism in the acetazolamide ligand of human carbonic anhydrase II (PDB entries 3hs4 and 4k0s), tautomerism in the 8HX ligand of urate oxidase (PDB entries 4n9s and 4n9m) and the protonation states of the catalytic aspartic acid found within the active site of an aspartic protease (PDB entry 2jjj). In each case, XModeScore applied to the X‐ray diffraction data is able to determine the correct protonation state as defined by the neutron diffraction data. The impact of QM‐based refinement versus conventional refinement on XModeScore is also discussed.