Molecular simulation study of CYP2B6 polymorphism with and without psoralen

CYP2B6 is a polymorphic enzyme with a large number of variants which may lead to functional changes in enzyme activity and substrate selectivity. In this study, CYP2B6 and its three variants with and without psoralen, a mechanism-based inactivator, were investigated using molecular simulation method. The obtained docking orientation of psoralen was in agreement with previously identified site of metabolism. Stability analysis showed that the three variants displayed more flexibility than CYP2B6.1, and CYP2B6.34 was the most flexible one without psoralen binding. However, in the presence of psoralen, CYP2B6.34 became more rigidity. Tunnel analysis indicates that the bottleneck change of tunnels may be correlated to the increased or decreased activity of variants. Binding free energy analysis shows that van der Waals interaction dominates the binding of psoralen. CYP2B6.34 has the highest affinity to psoralen with lowest binding free energy. Ile114, Phe115 and heme contribute largely to the binding of psoralen with CYP2B6.6, while Phe206 and Leu363 play important roles for CYP2B6.1 and CYP2B6.4. These computational observations suggest that the increased activity of CYP2B6.4 and reduced activity of CYP2B6.6 may be due to changes in regional structures.