Quantum Mechanics/Molecular Mechanics Modeling of Drug Metabolism: Mexiletine N-Hydroxylation by Cytochrome P450 1A2

The involvement of CYP1A2 in biodegradation of dioxins in pigs

2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) is one of the most harmful chemicals showing resistance to biodegradation. The majority of TCDD effects is mediated by the aryl hydrocarbon receptor (AhR) pathway. TCDD binding to AhR results in the activation of cytochrome P450 enzymes (CYP1A1, CYP1A2, CYP1B1) involved in dioxin biodegradation. The goal of the study was to explore the potentialrole of CYP1A2 in the metabolism of TCDD. We investigated a molecular structure of CYP1A2 and the binding selectivity and affinity between the pig CYP1A2 and: 1/ DiCDD or TCDD (dioxins differing in toxicity and biodegradability) or 2/ their selected metabolites. pCYP1A2 demonstrated higher affinity towards DiCDD and TCDD than other pCYP1 enzymes. All dioxin-pCYP1A2 complexes were found to be stabilized by hydrophobic interactions. The calculated distances between the heme oxygen and the dioxin carbon nearest to the oxygen, reflecting the hydroxylating potential of CYP1A2, were higher than in other pCYP1 enzymes. The distances between the heme iron and the nearest dioxin carbon exceeded 5 Å, a distance sufficient to allow the metabolites to leave the active site. However, the molecular dynamics simulations revealed that two access channels of CYP1A2 were closed upon binding the majority of the examined dioxins. Moreover, the binding of dioxin metabolites did not promote opening of channel S–an exit for hydroxylated products. It appears that the undesired changes in the behavior of access channels prevail over the hydroxylating potential of CYP1A2 towards TCDD and the favorable distances, ultimately trapping the metabolites at the enzyme’s active site.

Hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) Simulation: A Tool for Structure-based Drug Design and Discovery

Quantum mechanics (QM) is physics based theory which explains the physical properties of nature at the level of atoms and sub-atoms. Molecular mechanics (MM) construct molecular systems through the use of classical mechanics. So, hybrid quantum mechanics and molecular mechanics (QM/MM) when combined together can act as computer-based methods which can be used to calculate structure and property data of molecular structures. Hybrid QM/MM combines the strengths of QM with accuracy and MM with speed. QM/MM simulation can also be applied for the study of chemical process in solutions as well as in the proteins, and has a great scope in structure-based drug design (CADD) and discovery. Hybrid QM/MM also applied to HTS, to derive QSAR models and due to availability of many protein crystal structures; it has a great role in computational chemistry, especially in structure- and fragment-based drug design. Fused QM/MM simulations have been developed as a widespread method to explore chemical reactions in condensed phases. In QM/MM simulations, the quantum chemistry theory is used to treat the space in which the chemical reactions occur; however the rest is defined through molecular mechanics force field (MMFF). In this review, we have extensively reviewed recent literature pertaining to the use and applications of hybrid QM/MM simulations for ligand and structure-based computational methods for the design and discovery of therapeutic agents.

Reactive Metabolites from Thiazole-Containing Drugs: Quantum Chemical Insights into Biotransformation and Toxicity

Drugs containing thiazole and aminothiazole groups are known to generate reactive metabolites (RMs) catalyzed by cytochrome P450s (CYPs). These RMs can covalently modify essential cellular macromolecules and lead to toxicity and induce idiosyncratic adverse drug reactions. Molecular docking and quantum chemical hybrid DFT study were carried out to explore the molecular mechanisms involved in the biotransformation of thiazole (TZ) and aminothiazole (ATZ) groups leading to RM epoxide, S-oxide, N-oxide, and oxaziridine. The energy barrier required for the epoxidation is 13.63 kcal/mol, that is lower than that of S-oxidation, N-oxidation, and oxaziridine formation (14.56, 17.90, and 20.20, kcal/mol respectively). The presence of the amino group in ATZ further facilitates all the metabolic pathways, for example, the barrier for the epoxidation reaction is reduced by ∼2.5 kcal/mol. Some of the RMs/their isomers are highly electrophilic and tend to form covalent bonds with nucleophilic amino acids, finally leading to the formation of metabolic intermediate complexes (MICs). The energy profiles of these competitive pathways have also been explored.

Histidine versus Cysteine-Bearing Heme-Dependent Halogen Peroxidases: Parallels and Differences for Cl- Oxidation

The homodimeric myeloperoxidase (MPO) features a histidine as a proximal ligand and a sulfonium linkage covalently attaching the heme porphyrin ring to the protein. MPO is able to catalyze Cl^- oxidation with about the same efficiency as chloroperoxidase at pH 7.0. In this study, we seek to explore the parallels and differences between the histidine and cysteine heme-dependent halogen peroxidases. Transition states, reaction barriers, and relevant thermodynamic properties are calculated on protein models. Together with electronic structure calculations, it gives an overview of the reaction mechanisms and of the factors that determine the selectivity between one- and two-electron paths. Conclusions point to the innate oxidizing nature of MPO with the ester and sulfonium linkages hiking up the reactivity to enable chloride oxidation. The installation of a deprotonated imidazolate as a proximal ligand does not shift the equilibrium from one- to two-electron events without influencing the chemistry of the oxidation reaction.

Multiscale simulation approaches to modeling drug-protein binding

Simulations can provide detailed insight into the molecular processes involved in drug action, such as protein-ligand binding, and can therefore be a valuable tool for drug design and development. Processes with a large range of length and timescales may be involved, and understanding these different scales typically requires different types of simulation methodology. Ideally, simulations should be able to connect across scales, to analyze and predict how changes at one scale can influence another. Multiscale simulation methods, which combine different levels of treatment, are an emerging frontier with great potential in this area. Here we review multiscale frameworks of various types, and selected applications to biomolecular systems with a focus on drug-ligand binding.

Coupled electron and proton transfer in the piperidine drug metabolism pathway by the active species of cytochromes P450

KS-DFT and MSDFT studies reveal a novel CEPT step that triggers ring contraction of piperidines by P450.

Molecular Dynamics Simulation Framework to Probe the Binding Hypothesis of CYP3A4 Inhibitors

The Cytochrome P450 family of heme-containing proteins plays a major role in catalyzing phase I metabolic reactions, and the CYP3A4 subtype is responsible for the metabolism of many currently marketed drugs. Additionally, CYP3A4 has an inherent affinity for a broad spectrum of structurally diverse chemical entities, often leading to drug-drug interactions mediated by the inhibition or induction of the metabolic enzyme. The current study explores the binding of selected highly efficient CYP3A4 inhibitors by docking and molecular dynamics (MD) simulation protocols and their binding free energy calculated using the WaterSwap method. The results indicate the importance of binding pocket residues including Phe57, Arg105, Arg106, Ser119, Arg212, Phe213, Thr309, Ser312, Ala370, Arg372, Glu374, Gly481 and Leu483 for interaction with CYP3A4 inhibitors. The residue-wise decomposition of the binding free energy from the WaterSwap method revealed the importance of binding site residues Arg106 and Arg372 in the stabilization of all the selected CYP3A4-inhibitor complexes. The WaterSwap binding energies were further complemented with the MM(GB/PB)SA results and it was observed that the binding energies calculated by both methods do not differ significantly. Overall, our results could guide towards the use of multiple computational approaches to achieve a better understanding of CYP3A4 inhibition, subsequently leading to the design of highly specific and efficient new chemical entities with suitable ADMETox properties and reduced side effects.

Phe-125 and Phe-226 of pig cytochrome P450 1A2 stabilize the binding of aflatoxin B1 and 7-ethoxyresorufin through the key CH/π interactions

Cytochrome P450 1A2 (CYP1A2) plays important roles in the metabolism of many planar and aromatic drugs and also contributes to the bioactivation of aflatoxin B₁ (AFB₁) in vivo. To date, the structural basis for CYP1A2's preference to the planar substrates remains unclear. Herein, we investigated the structure-activity relationships for pig CYP1A2 catalyzing AFB₁ and 7-ethoxyresorufin (7-ER). A molecular docking study was performed based on a constructed model of pig CYP1A2, which predicted the contributions of Thr-118, Thr-124, Phe-125, Phe-226, Leu-260, and Asp-313 to the substrate catalysis. Site-directed mutagenesis and kinetic analyses exhibited the common grounds: Phe-125, Phe-226 and Asp-313 were vital to AFB₁ oxidation (including exo-epoxidation and 9A-hydroxylation) and ethoxyresorufin O-deethylation. Meanwhile, Phe-125 and Phe-226 formed CH/π interactions with AFB₁/7-ER, and Asp-313 formed hydrogen bonds with them. Based on other published reports, this study further emphasizes the critical roles of Phe-125 and Phe-226 in recognizing the planar substrates. Our findings highlight the structural basis of pig CYP1A2 specifically catalyzing AFB₁ and 7-ER, and may help to elucidate the underlying mechanism of CYP1A2's metabolic preference to the planar and aromatic substrates.

Multiple CH/π Interactions Maintain the Binding of Aflatoxin B₁ in the Active Cavity of Human Cytochrome P450 1A2

Human cytochrome P450 1A2 (CYP1A2) is one of the key CYPs that activate aflatoxin B₁ (AFB₁), a notorious mycotoxin, into carcinogenic exo-8,9-epoxides (AFBO) in the liver. Although the structure of CYP1A2 is available, the mechanism of CYP1A2-specific binding to AFB₁ has not been fully clarified. In this study, we used calculation biology to predict a model of CYP1A2 with AFB₁, where Thr-124, Phe-125, Phe-226, and Phe-260 possibly participate in the specific binding. Site-directed mutagenesis was performed to construct mutants T124A, F125A, F226A, and F260A. Escherichia coli-expressed recombinant proteins T124A, F226A, and F260A had active structures, while F125A did not. This was evidenced by Fe²⁺∙Carbon monoxide (CO)-reduced difference spectra and circular dichroism spectroscopy. Mutant F125A was expressed in HEK293T cells. Steady kinetic assays showed that T124A had enhanced activity towards AFB₁, while F125A, F226A, and F260A were significantly reduced in their ability to activate AFB₁, implying that hydrogen bonds between Thr-124 and AFB₁ were not important for substrate-specific binding, whereas Phe-125, Phe-226, and Phe-260 were essential for the process. The computation simulation and experimental results showed that the three key CH/π interactions between Phe-125, Phe-226, or Phe-260 and AFB₁ collectively maintained the stable binding of AFB₁ in the active cavity of CYP1A2.

Is CYP1B1 involved in the metabolism of dioxins in the pig?

2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) is the most difficult to biodegradate and the most toxic dioxin congener. Previously, we demonstrated in silico the ability of pig CYP1A1 to hydroxylate 2,7-dichlorodibenzo-p-dioxin (DiCDD), but not TCDD. To increase our knowledge concerning the low effectiveness of TCDD biodegradability, we analyzed in silico the binding selectivity and affinity between pig CYP1B1 and the two dioxins by means of molecular modeling. We also compared the effects of TCDD and DiCDD on CYP1B1 gene expression (qRT-PCR) and catalytic (EROD) activity in porcine granulosa cells. It was found that DiCDD and TCDD were stabilized within the pig CYP1B1 active site by hydrophobic interactions. The analysis of substrate channel availability revealed that both dioxins opened the exit channel S, allowing metabolites to leave the enzyme active site. Moreover, DiCDD and TCDD increased the CYP1B1 gene expression and catalytic activity in porcine granulosa cells. On the other hand, TCDD demonstrated higher than DiCDD calculated affinity to pig CYP1B1, hindering TCDD exit from the active site. The great distance between CYP1B1's heme and TCDD also might contribute to the lower hydroxylation effectiveness of TCDD compared to that of DiCDD. Moreover, the narrow active site of pig CYP1B1 may immobilize TCDD molecule, inhibiting its hydroxylation. The results of the access channel analysis and the distance from pig CYP1B1's heme to TCDD suggest that the metabolizing potential of pig CYP1B1 is higher than that of pig CYP1A1. However, this potential is probably not sufficiently high to considerably improve the slow TCDD biodegradation.

QM/MM Simulations with the Gaussian Electrostatic Model: A Density-based Polarizable Potential

The use of advanced polarizable potentials in quantum mechanical/molecular mechanical (QM/MM) simulations has been shown to improve the overall accuracy of the calculation. We have developed a density-based potential called the Gaussian electrostatic model (GEM), which has been shown to provide very accurate environments for QM wave functions in QM/MM. In this contribution we present a new implementation of QM/GEM that extends our implementation to include all components (Coulomb, exchange-repulsion, polarization, and dispersion) for the total intermolecular interaction energy in QM/MM calculations, except for the charge-transfer term. The accuracy of the method is tested using a subset of water dimers from the water dimer potential energy surface reported by Babin et al. ( J. Chem. Theory Comput. 2013 9, 5395-5403). Additionally, results of the new implementation are contrasted with results obtained with the classical AMOEBA potential. Our results indicate that GEM provides an accurate MM environment with average root-mean-square error <0.15 kcal/mol for every intermolecular interaction energy component compared with SAPT2+3/aug-cc-pVTZ reference calculations.

Oxygen Activation and Radical Transformations in Heme Proteins and Metalloporphyrins

As a result of the adaptation of life to an aerobic environment, nature has evolved a panoply of metalloproteins for oxidative metabolism and protection against reactive oxygen species. Despite the diverse structures and functions of these proteins, they share common mechanistic grounds. An open-shell transition metal like iron or copper is employed to interact with O2 and its derived intermediates such as hydrogen peroxide to afford a variety of metal-oxygen intermediates. These reactive intermediates, including metal-superoxo, -(hydro)peroxo, and high-valent metal-oxo species, are the basis for the various biological functions of O2-utilizing metalloproteins. Collectively, these processes are called oxygen activation. Much of our understanding of the reactivity of these reactive intermediates has come from the study of heme-containing proteins and related metalloporphyrin compounds. These studies not only have deepened our understanding of various functions of heme proteins, such as O2 storage and transport, degradation of reactive oxygen species, redox signaling, and biological oxygenation, etc., but also have driven the development of bioinorganic chemistry and biomimetic catalysis. In this review, we survey the range of O2 activation processes mediated by heme proteins and model compounds with a focus on recent progress in the characterization and reactivity of important iron-oxygen intermediates. Representative reactions initiated by these reactive intermediates as well as some context from prior decades will also be presented. We will discuss the fundamental mechanistic features of these transformations and delineate the underlying structural and electronic factors that contribute to the spectrum of reactivities that has been observed in nature as well as those that have been invented using these paradigms. Given the recent developments in biocatalysis for non-natural chemistries and the renaissance of radical chemistry in organic synthesis, we envision that new enzymatic and synthetic transformations will emerge based on the radical processes mediated by metalloproteins and their synthetic analogs.

Ab Initio QM/MM Modeling of the Rate-Limiting Proton Transfer Step in the Deamination of Tryptamine by Aromatic Amine Dehydrogenase

Aromatic amine dehydrogenase (AADH) and related enzymes are at the heart of debates on the roles of quantum tunneling and protein dynamics in catalysis. The reaction of tryptamine in AADH involves significant quantum tunneling in the rate-limiting proton transfer step, shown by large H/D primary kinetic isotope effects (KIEs), with unusual temperature dependence. We apply correlated ab initio combined quantum mechanics/molecular mechanics (QM/MM) methods, at levels up to local coupled cluster theory (LCCSD(T)/(aug)-cc-pVTZ), to calculate accurate potential energy surfaces for this reaction, which are necessary for quantitative analysis of tunneling contributions and reaction dynamics. Different levels of QM/MM treatment are tested. Multiple pathways are calculated with fully flexible transition state optimization by the climbing-image nudged elastic band method at the density functional QM/MM level. The average LCCSD(T) potential energy barriers to proton transfer are 16.7 and 14.0 kcal/mol for proton transfer to the two carboxylate atoms of the catalytic base, Asp128β. The results show that two similar, but distinct pathways are energetically accessible. These two pathways have different barriers, exothermicity and curvature, and should be considered in analyses of the temperature dependence of reaction and KIEs in AADH and other enzymes. These results provide a benchmark for this prototypical enzyme reaction and will be useful for developing empirical models, and analyzing experimental data, to distinguish between different conceptual models of enzyme catalysis.

Structural-functional adaptations of porcine CYP1A1 to metabolize polychlorinated dibenzo-p-dioxins

Polychlorinated dibenzo-p-dioxins (PCDDs) are widespread by-products of human industrial activity. They accumulate in tissues of animals and humans, exerting numerous adverse effects on different systems. In living organisms, dioxins are metabolized by enzymes of the cytochrome P450 family, including CYP1A1. Particular dioxin congeners differ in their toxicity level and ability to undergo biodegradation. Since the molecular mechanisms underlying dioxin susceptibility or resistance to biodegradation are unknown, in the present study the molecular interactions between five selected dioxins and porcine CYP1A1 protein were investigated. It was found that the ability of a dioxin to undergo CYP1A1-mediated degradation is associated mainly with the number and position of chlorine atoms in the dioxin molecule. Among all examined congeners, 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) demonstrated the highest affinity to CYP1A1 and, at the same time, the greatest distance to the active site of the enzyme. Interestingly, in contrast to other dioxins, the binding of the TCDD molecule to the porcine CYP1A1 active site resulted in a rapid and continuous closure of substrate channels. All the information may help to explain the extended half-life of TCDD in living organisms as well as its high toxicity.

Tryptophan lyase (NosL): mechanistic insights into amine dehydrogenation and carboxyl fragment migration by QM/MM calculations

QM/MM calculations suggest two feasible pathways for the breaking of the C–C bond of the substrate. The breaking of the Cα–Cβ bond leads to the final product, whereas the cleavage of the Cα–C bond will terminate in the EPR-trapped radical intermediate.