The hydroxylation of nonreactive CÀHb onds can be easily catalyzed by av ariety of metalloenzymes, especially cytochrome P450s (P450s). The mechanism of P450 mediated hydroxylation has been intensively studied, both experimentally and theoretically.H owever,u nderstandingt he regio-and stereoselectivities of substrates hydroxylated by P450s remains ag reat challenge. Herein, we use am ultiscale modeling approach to investigatet he selectivity of testosterone (TES) and dihydrotestosterone (DHT)h ydroxylation catalyzed by two important P450s,C YP3A4 and CYP19A1. For CYP3A4, two distinct binding modes for TES/DHTw ere predicted by dockings and molecular dynamics simulations, in whicht he experimentally identified sites of metabolism of TES/DHTc an access to the catalytic center.T he regio-and stereoselectivities of TES/DHTh ydroxylation were further evaluated by quantum mechanical andO NIOM calculations. For CYP19A1, we found that sites 1b,2 b and 19 can access the catalytic center, with the intrinsic reactivity 2b > 1b > 19. However, ourO NIOM calculations indicate that the hydroxylation is favored at site 19 for both TES and DHT,w hich is consistentw ith the experiments and reflects the importance of the catalytic environment in determining the selectivity. Our study unravels the mechanism underlyingt he selectivity of TES/DHT hydroxylationm ediated by CYP3A4 and CYP19A1a nd is helpful for understanding the selectivity of other substrates that are hydroxylatedb yP450s.[a] J.