Common and Distinct Interactions of Chemical Inhibitors with Cytochrome P450 CYP1A2, CYP2A6 and CYP2B6 Enzymes

BACKGROUND

Tobacco smoking is a leading cause of preventable disease and death globally. Nicotine is the main addictive component in tobacco. Nicotine is eliminated from the body by biotransformation in the liver to inactive metabolites. This reaction is catalyzed by the cytochrome P450 2A6 (CYP2A6) enzyme. Administering chemical inhibitors of CYP2A6 has been shown to slow down the elimination of nicotine with consequent reduction in number of cigarettes smoked. We have systematically developed small molecule CYP2A6 inhibitors with good balance between potency and CYP selectivity.

OBJECTIVE

During this process we have noticed that many potent CYP2A6 inhibitors also inhibit other human liver CYP forms, most notably CYP1A2 and CYP2B6. This study aimed at defining common and distinct features of ligand binding to CYP1A2, CYP2A6 and CYP2B6 active sites.

METHODS

We used our previous chemical inhibitor databases to construct improved 3-dimensional quantitative structureactivity relationship (3D-QSAR) models for CYP1A2, CYP2A6 and CYP2B6.

RESULTS

Combined 3D-QSAR and docking procedures yielded precise information about the common and distinct interactions of inhibitors and the enzyme active sites. Positioning of hydrogen bond donor/acceptor atoms and the shape and volume of the compound defined the potency and specificity of inhibition. A novel potent and selective CYP1A2 inhibitor was found.

CONCLUSION

This in silico approach will provide a means for very rapid and high throughput prediction of cross-inhibition of these three CYP enzymes.

Structural and biophysical characterization of human cytochromes P450 2B6 and 2A6 bound to volatile hydrocarbons: analysis and comparison

X-ray crystal structures of complexes of cytochromes CYP2B6 and CYP2A6 with the monoterpene sabinene revealed two distinct binding modes in the active sites. In CYP2B6, sabinene positioned itself with the putative oxidation site located closer to the heme iron. In contrast, sabinene was found in an alternate conformation in the more compact CYP2A6, where the larger hydrophobic side chains resulted in a significantly reduced active-site cavity. Furthermore, results from isothermal titration calorimetry indicated a much more substantial contribution of favorable enthalpy to sabinene binding to CYP2B6 as opposed to CYP2A6, consistent with the previous observations with (+)-α-pinene. Structural analysis of CYP2B6 complexes with sabinene and the structurally similar (3)-carene and comparison with previously solved structures revealed how the movement of the F206 side chain influences the volume of the binding pocket. In addition, retrospective analysis of prior structures revealed that ligands containing -Cl and -NH functional groups adopted a distinct orientation in the CYP2B active site compared with other ligands. This binding mode may reflect the formation of Cl-π or NH-π bonds with aromatic rings in the active site, which serve as important contributors to protein-ligand binding affinity and specificity. Overall, the findings from multiple techniques illustrate how drugs metabolizing CYP2B6 and CYP2A6 handle a common hydrocarbon found in the environment. The study also provides insight into the role of specific functional groups of the ligand that may influence the binding to CYP2B6.

Understanding a substrate's product regioselectivity in a family of enzymes: a case study of acetaminophen binding in cytochrome P450s

Product regioselectivity as influenced by molecular recognition is a key aspect of enzyme catalysis. We applied large-scale two-dimensional (2D) umbrella sampling (USP) simulations to characterize acetaminophen (APAP) binding in the active sites of the family of Cytochrome P450 (CYP) enzymes as a case study to show the different regioselectivity exhibited by a single substrate in comparative enzymes. Our results successfully explain the experimentally observed product regioselectivity for all five human CYPs included in this study, demonstrating that binding events play an important role in determining regioselectivity. In CYP2C9 and CYP3A4, weak interactions in an overall large active site cavity result in a fairly small binding free energy difference between APAP reactive binding states, consistent with experimental results that show little preference for resulting metabolites. In contrast, in CYP1A2 and CYP2E1, APAP is strongly restrained by a compact binding pocket, leading to a preferred binding conformation. The calculated binding equilibrium of APAP within the compact active site of CYP2A6 is able to predict the experimentally documented product ratios and is also applied to explain APAP regioselectivity in CYP1A2 and CYP2C9. APAP regioselectivity seems to be related to the selectivity for one binding conformation over another binding conformation as dictated by the size and shape of the active site. Additionally, unlike docking and molecular dynamics (MD), our free energy calculations successfully reproduced a unique APAP pose in CYP3A4 that had been reported experimentally, suggesting this approach is well suited to find the realistic binding pose and the lowest-energy starting structure for studying the chemical reaction step in the future.