Monoprotic and diprotic NH tautomerism in reduced oligoazaacenes, the pyrazinacenes, was studied by using first principles simulations. Stepwise reductions in the metadynamics-sampled free energy profile were observed during consecutive monoprotic tautomerizations, with energy barriers gradually reducing with increasing proton separation during monoprotic processes. This is accompanied by an increasing contribution from the quinoidal electronic structure, as evidenced by the computed highest occupied molecular orbital (HOMO) structure. An unusual odd-even effect in the free energy profiles is also observed upon changing the length of the pyrazinacene. Calculated HOMO structures reveal an increasing tendency for delocalization of pyrazine lone pairs with an increasing number of ring annelations. The influence of tautomerism on the pyrazine lone pair delocalization was also observed. Tautomers with protons situated centrally on the pyrazinacene backbone are predicted to be more stable due to a combination of (enamine) delocalization and a loss of Clar sextet resonance stabilization in tautomers with protons at terminal pyrazine rings. Experimental evidence suggesting the structure of pyrazinacene tautomers is included and discussed as a support to the calculation.