Employing a hybrid QM/MM simulation we explored the reaction dynamics of the hydroxyl radical scavenging activity of daidzein, a soy isoflavone. Our simulations illustrate that the highly reactive hydroxyl radical can participate in hydrogen abstraction reaction with both OH functional groups of daidzein and can form stable daidzein radicals. We found that the reaction involving the 4'-OH site of daidzein is energetically favorable over the other reaction pathway involving the 7-OH site of daidzein by similar to 29 kcal/mol. The high enthalpic stabilization involved in daidzein radical formation at the 4'-OH site can be partly attributed to better solvation through hydrogen-bonding interactions with water and higher electron density delocalization of radical over the adjacent aromatic ring. As evident from the QM/MM dynamics, both HAT pathways led to formation of ketones at the 7-OH and 4'-OH sites of daidzein, respectively, and the adjacent aromatic rings appear in a p-quinonoid form, a highly stable resonating structure. The suitability of the QM/MM methodology to study the reaction mechanism, identification of intermediate states, and pathways of flavonoid radical stabilization reported here opens up a new possibility to study a similar reaction mechanism in other systems.