Variational transition state theory calculations including Multidimensional tunneling (VTST/MT) for excited-state tautomerization in the 1:1 7-azaindole:H2O complex were performed. Electronic Structures and energies for reactant, product, transition state, and potential energy curves along the reaction coordinate were computed at the CASSCF(10,9)/6-31G(d,p) level of theory. The potential energies were corrected by second-order multireference perturbation theory to take the dynamic electron correlation into consideration. The final potential energy Curves along, the reaction coordinate were generated at the MRPT2//CASSCF(10,9)/6-31G(d,p) level. Two protons in the excited-state tautomerization are transferred concertedly, albeit asynchronously. The Position of the variational transition state is very different from the conventional transition state, and is highly dependent on isotopic Substitution. Rate constants were calculated using VTST/MT, and were oil the order of 10(-6) s(-1) at room temperature. The HH/DD kinetic isotope effects are consistent with experimental observations; consideration of both tunneling and variational effects was essential to predict the experimental values correctly.