Binding free energies of inhibitors to iron porphyrin complex as a model for Cytochrome P450

Bidirectional Charge Transfer at the Heme Iron in Reversible and Quasi-irreversible Cytochrome P450 Inhibition

The coordination bonding between inhibitor ligands and heme iron plays a critical role in disrupting the essential catalytic functions of cytochrome P450 enzymes (P450s). Despite its intrinsic importance and consequential implications for human health, our current understanding of coordination bonding in P450 inhibition remains limited. To address this knowledge gap, we conducted a systematic theoretical analysis of the complexes between a ferric or a ferrous heme model and representative inhibitor ligands. Specifically, we evaluated the charge-transfer (CT) effect within these complexes by employing a series of theoretical methods based on density functional theory (DFT). Through a comprehensive analysis, we unveiled the relative significance of ligand-to-heme forward CT in the ferric and ferrous complexes of reversible inhibitors. In contrast, backward CT dominates over forward CT in the ferrous heme complexes of quasi-irreversible inhibitors. Further analysis using the compact frontier orbital method underscores the elevated electron-accepting abilities of quasi-irreversible inhibitors for π backdonation, which greatly amplifies their binding affinity for the ferrous heme. This study sheds light on the intricate mechanisms underlying P450 inhibition and provides valuable insights for future inhibitor design and development.

Ni2+ and Cu2+ complexes of N-(2,6-dichlorophenyl)-N-mesityl formamidine dithiocarbamate structural and functional properties as CYP3A4 potential substrates

Metal compounds continued to attract diverse applications due to their malleability in several capacities. In this study, we present our findings on the crystal structures and functional properties of Ni2+ and Cu2+ complexes of N'-(2,6-dichlorophenyl)-N-mesitylformamidine dithiocarbamate (L) comprising [Ni-(L)2] (1) and [Cu-(L)2] (2) with a four-coordinate metal center. We established the two complex structures through 1H and 13C nuclear magnetic resonance (NMR), elemental, and single-crystal X-ray analysis. The analyses showed that the two complexes are isomorphous, having P21/c as a space group and a unit-cell similarity index (π) of 0.002. The two complexes conform to a distorted square planar geometry around the metal centers. The calculated and experimental data, including bond lengths, angles, and NMR values, are similar. Hirshfeld surface analysis revealed the variational contribution of the different types of intermolecular contacts driven by the crystal lattice of the two solvated complexes. Our knowledge of the potential biological implication of these structures enabled us to probe the compounds as prospective CYP3A4 inhibitors. This approach mimics current trends in pharmaceutical design and biomedicine by incorporating potentially active molecules into various media to predict their biological efficacies. The simulations show appreciable binding of compounds 1 and 2 to CYP3A4 with average interaction energies of -97 and -87 kcal/mol, respectively. The protein attains at least five conformational states in the three studied models using a Gaussian Mixture Model-based clustering and free energy prediction. Electric field analysis shows the crucial residues to substrate binding at the active site, enabling CYP3A4 structure to function prediction. The predicted inhibition with these Ni2+ and Cu2+ complexes indicates that CYP3A4 overexpression in a diseased state like cancer would reduce, thereby increasing the chemotherapeutic compounds' shelf-lives for adsorption. This multidimensional study addresses various aspects of molecular metal electronics, including their application as substrate-mimicking inhibitors. The outcome would enable further research on bio-metal compounds of critical potential.

Energy Decomposition Analysis of the Nature of Coordination Bonding at the Heme Iron Center in Cytochrome P450 Inhibition

Drug compounds or their metabolic intermediates (MIs) sometimes inhibit the function of cytochrome P450 enzymes (P450s) by forming a coordination bond with the Fe(III) heme or Fe(II) heme of P450s. Such inhibition is one of the major causes of drug-drug interactions (DDIs), a subject of longstanding academic and practical interest. However, such coordination bonding is not fully understood at the quantum mechanical level, thus hampering rational improvement of the accuracy of DDI-related predictions. In this work, we employed density functional theory (DFT) and the generalized Kohn-Sham energy decomposition analysis (GKS-EDA) scheme to investigate the nature of the coordination bonding formed in the reversible and quasi-irreversible inhibition of P450s. The GKS-EDA results highlighted a previously unrecognized role of the electron correlation effect in P450 inhibition. The correlation effect tends to be larger in Fe(II) complexes of MI-type inhibitors and is particularly prominent for the nitrosoalkane ligand. An additional natural bond orbital (NBO) analysis provided insight into the relative significance of the σ donation and π backdonation effects in various heme-inhibitor complexes.

Use of density functional theory in drug metabolism studies

INTRODUCTION

The cytochrome P450 enzymes (CYPs) metabolize many drug compounds. They catalyze a wide variety of reactions, and potentially, a large number of different metabolites can be generated. Density functional theory (DFT) has, over the past decade, been shown to be a powerful tool to rationalize and predict the possible metabolites generated by the CYPs as well as other drug-metabolizing enzymes.

AREAS COVERED

We review applications of DFT on reactions performed by the CYPs and other drug-metabolizing enzymes able to perform oxidation reactions, with an emphasis on predicting which metabolites are produced. We also cover calculations of binding energies for complexes in which the ligands interact directly with the heme iron atom.

EXPERT OPINION

DFT is a useful tool for prediction of the site of metabolism. The use of small models of the enzymes work surprisingly well for most CYP isoforms. This is probably due to the fact that the binding of the substrates is not the major determinant. When binding of the substrate plays a significant role, the well-known issue of determining the free energy of binding is the challenge. How approaches taking the protein environment into account, like docking, MD and QM/MM, can be used are discussed.

Understanding the mechanism of cytochrome P450 3A4: recent advances and remaining problems

Cytochromes P450 (CYPs) represent a diverse group of heme-thiolate proteins found in almost all organisms. CYPs share a common protein fold but differ in substrate selectivity and catalyze a wide variety of monooxygenation reactions via activation of molecular oxygen. Among 57 human P450s, the 3A4 isoform (CYP3A4) is the most abundant and the most important because it metabolizes the majority of administered drugs. A remarkable feature of CYP3A4 is its extreme promiscuity in substrate specificity and cooperative substrate binding, which often leads to undesirable drug-drug interactions and toxic side effects. Owing to its importance in drug development and therapy, CYP3A4 has been the most extensively studied mammalian P450. In this review we provide an overview on recent progress and remaining problems in the CYP3A4 research.