Molecular Dynamic Investigations of the Mutational Effects on Structural Characteristics and Tunnel Geometry in CYP17A1

Cytochrome P450 (CYP) 17A1 is a dual-function monooxygenase with a critical role in the synthesis of many human steroid hormones. The enzyme is an important target for the treatment of breast and prostate cancers that proliferate in response to estrogens and androgens. Despite the ample experimental mutagenesis data, the molecular origin and the structural motifs for the enzymatic activities deficiencies have not been rationalized at the atomic resolution. To this end, we have investigated the effects on structural characteristics and tunnel geometry upon single point mutations in CYP17A1. The MD simulation results combined with PMF calculations and MM-GBSA calculations render an access mechanism which encapsulates the effects of mutations on the changes in both structural flexibility and tunnel dynamics, bridging the gap between the theory and the experimentally observed results of enzymatic activity decrease. The underlying molecular mechanism of the heterogeneities in open/closed conformational changes, as well as the wider opening of their respective major tunnels between wt17A1 and two mutants, may be attributed to the closer distances of hydrophobic residues or the disruption of a hydrophobic core. The knowledge of ligand binding characteristics and key residues contributions could guide future experimental and computational work on CYPs so that desirable changes in their enzymatic activities may be achieved. The present study provides important insights into the structure-function relationships of CYP17A1 protein, which could contribute to further understanding about 17-hydroxylase deficiencies and may also improve the understanding of polycystic ovary disease.