Molecular dynamics of substrate complexes with hamster cytochrome P450c17 (CYP17): mechanistic approach to understanding substrate binding and activities

The cytochrome P450c17 isoforms from various animal species have different substrate selectivity, especially for 17,20-lyase activity. In particular, the human P450c17 selectively produces dehydroepiandrosterone with little androstenedione (AD). Hamster P450c17, on the other hand, produces both of these steroids at comparable rates. We thus investigated if computational analysis could explain the difference in activity profiles. Therefore, we inserted the four P450c17 substrates-pregnenolone, progesterone, and their 17alpha-hydroxylated forms-inside our hamster P450c17 model, which we derived from our human P450c17 model based on the crystal structure of P450BMP. We performed molecular dynamics (MD) simulations on the complexes and analyzed the resultant trajectories to identify amino acids that interact with substrates. Starting with substrates in two different orientations, we obtained two sets of binding trajectories in each case. The first set of trajectories reveal structural rearrangements that occur during binding, whereas the second set of trajectories reflects substrate orientations during catalysis. Our modeling suggests that three distinct steps are required for substrate selectivity and binding to the hamster P450c17: (1) recognition of the substrate at the putative substrate entrance, characterized by a pocket at the surface of the hamster P450c17 containing charged residues R96 and D116; (2) entry of the substrate into the active site, in an intermediate position directed by possible hydrogen bonding of the substrates with the heme D-ring propionate group, R96, R440, and T306; followed by (3) 90 degrees counterclockwise rotation of the substrates, positioning them in optimal position for reactivity, a process that may be directed by hydrogen bonding to the 110-112 region of the hamster P450c17. With some substrates, we obtained trajectories which suggest that major distortions in the I-helix and opening of the H-I loop occur during substrate binding. In conclusion, these modeling exercises provide insight to possible structural reorganizations that occur during substrate binding and suggest that amino acids that participate in three distinct steps of this process may all contribute to substrate binding and activity.