A pair of bases in the DNA, say, A-T, have a tautomerized form A*-T* (resulting from switching the sides of both protons along the hydrogen bridges).
I have studied how, by means of DFT computations, one can compute the energy difference between the normal and the anomalous forms. But, another point of interest is the stability of the tautomer.
I have read several times that those anomalous tautomers "decay" to the normal A-T form very quickly (much quicker than the time it lasts a round of DNA replication), which is interesting because it implies that tautomerization is not likely to be a source of spontaneous mutations.
What I would like to learn is how can one compute this expected life-time of the anomalous tautomer
I have an idea (which I don't know if it could be rigth or if maybe is completely worthless), which involves making a Born-Oppenheimer separation for the two protons (''after'' the DFT computation of the electronic ground state) and then study the non-adiabatic transitions between the different vibronic states of the couple of protons. I am happy to provide details, if someone whises so.
But perhaps "my method" is too messy, or naive (or both) and / or there are other standar approaches to compute approximately the expected mean life of anomalous tautomers such as the mentioned A*-T*.
I would be very grateful if someone could apport some insigth and/or references into how are these expected life-times calculated typically in computational chemistry!