The cytochromes CYP3A4 and CYP3A5 share 84% sequence identity, but they exhibit different catalytic activities toward some substrates. Schisantherin E (SE) was recently identified as a selective substrate of CYP3A5, which exhibited catalytic efficiency that was more than 23 times higher than CYP3A4. At present, however, the structural determinants responsible for the different catalytic activities of the two enzymes toward SE have not been fully understood. In this study, a combination of molecular docking, molecular dynamic simulations, and binding free energy calculation was performed on the CYP3A4/CYP3A5-SE systems to investigate the issue. The results demonstrate that Ser119 in CYP3A4 and Glu374 in CYP3A5 formed direct hydrogen bonding with SE, respectively. Additionally, one water molecule located between the B-C loop and the I helix mediated different hydrogen-bonding networks between CYP3A4/3A5 and SE. The residue differences (Phe/Leu108 and Leu/Phe210) triggered the distinct conformational changes of the Phe-cluster residues, especially Phe213 and Phe215, which formed stronger hydrophobic interactions with SE in CYP3A5. The calculated binding free energies were consistent with the experimental results.