Semiclassical molecular dynamics simulations have been combined with quantum chemistry calculations to provide detailed modeling of the methane and ethane hydroxylation reactions catalyzed by the hydroxylase enzymes of the soluble methane monooxygenase system. The experimental distribution of enantiomeric alcohols in the reaction of ethanes made chiral by the use of hydrogen isotopes is quantitatively reproduced and explained. The reaction dynamics involve a mixture of concerted and bound radical trajectories, and we characterize each of these reactive channels in detail. Diffusion of the bound radical intermediate at the active site core determines the global rate constant. The results also provide a qualitative rationale for the lack of ring-opened products derived from certain radical clock substrate probes and for the relative race constants and kinetic isotope effects exhibited by a variety of substrates.