Predicting reaction barriers of hydrogen atom transfer in proteins

Hydrogen atom transfer (HAT) reactions are important in many biological systems. As these reactions are hard to observe experimentally, it is of high interest to shed light on them using simulations. Here, we present a machine learning model based on graph neural networks for the prediction of activation energies of HAT reactions in proteins. It is trained on more than 17,000 energy barriers calculated using hybrid density functional theory. We built and evaluated the model in the context of HAT in collagen, but the same workflow can easily be applied to HAT reactions in other biological or synthetic polymers. We obtain for relevant reactions (small reaction distances) a model with good predictive power (R2 <0.9 and mean absolute error of <3 kcal/mol). As the inference speed is high, this model enables evaluations of many chemical situations in rapid succession. When combined with molecular dynamics in a kinetic Monte-Carlo scheme, the model paves the way toward reactive simulations.