Hypothesis of a proton switch in QM/MM modelling of interaction of dUMP analogues with thymidylate synthase

Quantum mechanical, molecular mechanics and molecular dynamics (MD) methods were used to investigate initial steps of 2'-deoxyuridine-5'-monophosphate (dUMP) methylation catalysed by the thymidylate synthase (TS) enzyme. The amino acid residues surrounding the active site within a 10 radius sphere were modelled with the combined quantum mechanical (B3LYP/LANL2DZ) and molecular mechanics ONIOM double-layer method. The results indicated the initial nucleophilic attack of Cys146 on dUMP to be concerted with formation of a hydrogen bond to the oxygen O4 of dUMP. Moreover, the proton in the vicinity of the O4 atom appears to act as a 'proton switch': if a proton is present near O4, it stabilises the S(Cys146)-C6(dUMP) sulphur-carbon bond, but if it is absent, the sulphur-carbon bond does not form. If the O4 oxygen is replaced by sulphur atom, the 'switch effect' does not occur. The suggested correlation between the strength of hydrogen bond involving O4 oxygen and the ability of dUMP to form bonds at C6 corresponds well to the crystal structures of TS complexes available in the Protein Data Bank. In the vast majority of crystal structures, the presence of the S(Cys146)-C6(dUMP) bond was coupled with the presence of hydrogen bond between the dUMP O4 atom and the conserved Asn177. The 'proton switch' hypothesis is supported also by the results of MD studies of TS binary complexes, suggesting that average distance separating S(Cys146) and C6(dUMP) becomes distinctly shorter in the presence of hydrogen bonding between Asn177 and O4.