A Comprehensive Computational Study on the Thermodynamics and Kinetics of Tetrahydrobiopterin Regeneration Process

One of the most crucial enzymatic cofactors in the human body is tetrahydrobiopterin, which is acquired through biological synthesis and self-regeneration. During this regenerative process, it undergoes oxidation, deprotonation, further oxidation, and subsequent deprotonation, resulting in the formation of quinonoid-dihydrobiopterin, which then undergoes tautomerization to yield dihydrobiopterin. This study presents the thermodynamic and kinetic properties associated with each stage of the regeneration process using theoretical calculations. The redox potentials for oxidation steps and the pKa values for deprotonation steps are determined employing the Born-Haber cycle and the direct change of free energy in implicit solvent models. The redox metabolites are characterized and confirmed from their calculated absorption spectra using the time-dependent density functional theory method. For the tautomerization steps, an IRC calculation is executed, and rate constants are computed using Eyring's Transition State Theory (TST). The tunnelling probability of the H atom during the tautomerization process is incorporated using Wigner's tunnelling correction in the calculation of the rate constant. Notably, we identify the N3 atom as the most probable deprotonation site for H₃B⁺ and predict its geometry based on our calculations. Furthermore, we elucidate the spectral properties of intermediates involved in the regeneration process, highlighting key electronic transitions responsible for their excitations. Our results indicate that each step of tautomerization occurs along vibrational bending modes. We have observed that these tautomerization processes have high activation energies by optimising transition states. Additionally, considering tunnelling correction can significantly affect the reaction rates associated with these processes. These results provide a comprehensive understanding of the thermodynamics and kinetics of the regeneration process of tetrahydrobiopterin, which will help in the modulation of its biological activity.