Molecular dynamics (MD) simulations at normal and high temperature were used to study the flexibility and malleability of three microsomal cytochromes P450 (CYPs): CYP3A4, CYP2C9, and CYP2A6. Comparison of B-factors (describing the atomic fluctuations) between X-ray and MD data shows that the X-ray B-factors are significantly lower in the regions where the crystal contacts occur than for other regions. Consequently, the conclusions about CYP flexibility based solely on the X-ray data might be misleading. Comparison of flexibility patterns of the three CYPs enabled common features and variations in flexibility and malleability of the studied CYPs to be identified. The previously described pattern of flexibility in topological elements of microsomal CYPs (a rigid heme binding core, a malleable distal side and intermediately flexible proximal side) was confirmed. These topological features provide an important combination of high stereo- and regio-specificity (mediated by the relative rigidity in the neighborhood of the heme), together with high substrate promiscuity due to the more flexible active site and the malleability of the distal side. The data acquired here show that the malleability of the three studied CYPs correlates with their substrate specificity: CYP2A6 has a narrow substrate range and is the most rigid, CYP3A4 is the most promiscuous CYP known and is the most malleable, and CYP2C9 is intermediate in terms of both its substrate specificity and malleability. Thus, the malleability of CYPs is probably a major determinant of their substrate specificity.
Aryl Hydrocarbon Hydroxylases Binding Sites