We have performed hybrid density functional theory (DFT) calculations to investigate how chemical equilibria can be described in the S-3 state of the oxygen-evolving complex in photosystem II. For a chosen 340-atom model, 1 stable and 11 metastable intermediates have been identified within the range of 13 kcal mol(1) that differ in protonation, charge, spin, and conformational states. The results imply that reversible interconversion of these intermediates gives rise to dynamic equilibria that involve processes with relocations of protons and electrons residing in the Mn4CaO5 cluster, as well as bound water ligands, with concomitant large changes in the cluster geometry. Such proton tautomerism and redox isomerism are responsible for reversible activation/deactivation processes of substrate oxygen species, through which MnO and OO bonds are transiently ruptured and formed. These results may allow for a tentative interpretation of kinetic data on substrate water exchange on the order of seconds at room temperature, as measured by time-resolved mass spectrometry. The reliability of the hybrid DFT method for the multielectron redox reaction in such an intricate system is also addressed.