Kinetic deuterium isotope effects for 7-alkoxycoumarin O-dealkylation reactions catalyzed by human cytochromes P450 and in liver microsomes. Rate-limiting C-H bond breaking in cytochrome P450 1A2 substrate oxidation

Temperature-Dependent Kinetic Isotope Effects in R67 Dihydrofolate Reductase from Path-Integral Simulations

Calculation of temperature-dependent kinetic isotope effects (KIE) in enzymes presents a significant theoretical challenge. Additionally, it is not trivial to identify enzymes with available experimental accurate intrinsic KIEs in a range of temperatures. In the current work, we present a theoretical study of KIEs in the primitive R67 dihydrofolate reductase (DHFR) enzyme and compare with experimental work. The advantage of R67 DHFR is its significantly lower kinetic complexity compared to more evolved DHFR isoforms. We employ mass-perturbation-based path-integral simulations in conjunction with umbrella sampling and a hybrid quantum mechanics-molecular mechanics Hamiltonian. We obtain temperature-dependent KIEs in good agreement with experiments and ascribe the temperature-dependent KIEs primarily to zero-point energy effects. The active site in the primitive enzyme is found to be poorly preorganized, which allows excessive water access to the active site and results in loosely bound reacting ligands.

A Modified Arrhenius Approach to Thermodynamically Study Regioselectivity in Cytochrome P450-Catalyzed Substrate Conversion

The regio- (and stereo-)selectivity and specific activity of cytochrome P450s are determined by the accessibility of potential sites of metabolism (SOMs) of the bound substrate relative to the heme, and the activation barrier of the regioselective oxidation reaction(s). The accessibility of potential SOMs depends on the relative binding free energy (ΔΔGbind ) of the catalytically active substrate-binding poses, and the probability of the substrate to adopt a transition-state geometry. An established experimental method to measure activation energies of enzymatic reactions is the analysis of reaction rate constants at different temperatures and the construction of Arrhenius plots. This is a challenge for multistep P450-catalyzed processes that involve redox partners. We introduce a modified Arrhenius approach to overcome the limitations in studying P450 selectivity, which can be applied in multiproduct enzyme catalysis. Our approach gives combined information on relative activation energies, ΔΔGbind values, and collision entropies, yielding direct insight into the basis of selectivity in substrate conversion.

Heterologous expression and functional characterization of the NADPH-cytochrome P450 reductase from Capsicum annuum

Two NADPH-cytochrome P450 reductase (CPR) genes (CaCPR1 and CaCPR2) were isolated from hot pepper (Capsicum annuum L. cv. Bukang). At the red ripe stage, the expression level of CaCPR1 was more than 6-fold greater than that in leaves or flowers. It gradually increased during fruit ripening. The CaCPR2 gene seemed to be expressed constitutively in all of the tested tissues. To investigate the enzymatic properties of CaCPR1, the cDNA of CaCPR1 was heterologously expressed in Escherichia coli without any modification of amino acid sequences, and CaCPR1 was purified. The enzymatic properties of CaCPR1 were confirmed using cytochrome c and cytochrome b5 as protein substrates. The CaCPR1 could support human CYP1A2-catalyzed reaction. It also reduced tetrazolium salts and ferricyanide. These results show that CaCPR1 is the major CPR in most hot pepper tissues. It is suggested that the CaCPR1 can be used a prototype for studying biological functions and biotechnological applications of plant CPRs.

Characterization of the active site properties of CYP4F12

Cytochrome P450 4F12 is a drug-metabolizing enzyme that is primarily expressed in the liver, kidney, colon, small intestine, and heart. The properties of CYP4F12 that may impart an increased catalytic selectivity (decreased promiscuity) were explored through in vitro metabolite elucidation, kinetic isotope effect experiments, and computational modeling of the CYP4F12 active site. By using astemizole as a probe substrate for CYP4F12 and CYP3A4, it was observed that although CYP4F12 favored astemizole O-demethylation as the primary route of metabolism, CYP3A4 was capable of metabolizing astemizole at multiple sites on the molecule. Deuteration of astemizole at the site of O-demethylation resulted in an isotope effect of 7.1 as well as an 8.3-fold decrease in the rate of clearance for astemizole by CYP4F12. Conversely, although an isotope effect of 3.8 was observed for the formation of the O-desmethyl metabolite when deuterated astemizole was metabolized by CYP3A4, there was no decrease in the clearance of astemizole. Development of a homology model of CYP4F12 based on the crystal structure of cytochrome P450 BM3 predicted an active site volume for CYP4F12 that was approximately 76% of the active site volume of CYP3A4. As predicted, multiple favorable binding orientations were available for astemizole docked into the active site of CYP3A4, but only a single binding orientation with the site of O-demethylation oriented toward the heme was identified for CYP4F12. Overall, it appears that although CYP4F12 may be capable of binding similar ligands to other cytochrome P450 enzymes such as CYP3A4, the ability to achieve catalytically favorable orientations may be inherently more difficult because of the increased steric constraints of the CYP4F12 active site.

Application of a deuterium replacement strategy to modulate the pharmacokinetics of 7-(3,5-dimethyl-1H-1,2,4-triazol-1-yl)-3-(4-methoxy-2-methylphenyl)-2,6-dimethylpyrazolo[5,1-b]oxazole, a novel CRF1 antagonist

Deuterium isotope effects were evaluated as a strategy to optimize the pharmacokinetics of 7-(3,5-dimethyl-1H-1,2,4-triazol-1-yl)-3-(4-methoxy-2-methylphenyl)-2,6-dimethylpyrazolo[5,1-b]oxazole (NVS-CRF38), a novel corticotropin-releasing factor receptor 1 (CRF1) antagonist. In an attempt to suppress O-demethylation of NVS-CRF38 without losing activity against the CRF1 receptor, the protons at the site of metabolism were replaced with deuterium. For in vitro and in vivo studies, intrinsic primary isotope effects (KH/KD) were determined by the ratio of intrinsic clearance (CLint) obtained for NVS-CRF38 and deuterated NVS-CRF38. In vitro kinetic isotope effects (KH/KD) were more pronounced when CLint values were calculated based on the rate of formation of the O-desmethyl metabolite (KH/KD ∼7) compared with the substrate depletion method (KH/KD ∼2). In vivo isotope effects were measured in rats after intravenous (1 mg/kg) and oral (10 mg/kg) administration. For both administration routes, isotope effects calculated from in vivo CLint corresponding to all biotransformation pathways were lower (KH/KD ∼2) compared with CLint values calculated from the O-demethylation reaction alone (KH/KD ∼7). Comparative metabolite identification studies were undertaken using rat and human microsomes to explore the potential for metabolic switching. As expected, a marked reduction of the O-demethylated metabolite was observed for NVS-CRF38; however, levels of NVS-CRF38's other metabolites increased, compensating to some extent for the isotope effect.

The Flavin-Containing Reductase Domain of Cytochrome P450 BM3 Acts as a Surrogate for Mammalian NADPH-P450 Reductase

Cytochrome P450 BM3 (CYP102A1) from Bacillus megaterium is a self-sufficient monooxygenase that consists of a heme domain and FAD/FMN-containing reductase domain (BMR). In this report, the reduction of 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) and 5-cyano-2,3-ditolyl tetrazolium chloride (CTC) by BMR was evaluated as a method for monitoring BMR activity. The electron transfer proceeds from NADPH to BMR and then to BMR substrates, MTT and CTC. MTT and CTC are monotetrazolium salts that form formazans upon reduction. The reduction of MTT and CTC followed classical Michaelis-Menten kinetics (k_cat=4120 min⁻¹, K_m=77 μM for MTT and k_cat=6580 min⁻¹, K_m=51 μM for CTC). Our continuous assay using MTT and CTC allows the simple, rapid measurement of BMR activity. The BMR was able to metabolize mitomycin C and doxorubicin, which are anticancer drug substrates for CPR, producing the same metabolites as those produced by CPR. Moreover, the BMR was able to interact with CYP1A2 and transfer electrons to promote the oxidation reactions of substrates by CYP1A2 and CYP2E1 in humans. The results of this study suggest the possibility of the utilization of BMR as a surrogate for mammalian CPR.

Measuring kinetic isotope effects in enzyme reactions using time-resolved electrospray mass spectrometry

Kinetic isotope effect (KIE) measurements are a powerful tool for studying enzyme mechanisms; they can provide insights into microscopic catalytic processes and even structural constraints for transition states. However, KIEs have not come into widespread use in enzymology, due in large part to the requirement for prohibitively cumbersome experimental procedures and daunting analytical frameworks. In this work, we introduce time-resolved electrospray ionization mass spectrometry (TRESI-MS) as a straightforward, precise, and inexpensive method for measuring KIEs. Neither radioisotopes nor large amounts of material are needed and kinetic measurements for isotopically "labeled" and "unlabeled" species are acquired simultaneously in a single "competitive" assay. The approach is demonstrated first using a relatively large isotope effect associated with yeast alcohol dehydrogenase (YADH) catalyzed oxidation of ethanol. The measured macroscopic KIE of 2.19 ± 0.05 is consistent with comparable measurements in the literature but cannot be interpreted in a way that provides insights into isotope effects in individual microscopic steps. To demonstrate the ability of TRESI-MS to directly measure intrinsic KIEs and to characterize the precision of the technique, we measure a much smaller (12)C/(13)C KIE associated specifically with presteady state acylation of chymotrypsin during hydrolysis of an ester substrate.

Increase of human CYP1B1 activities by acidic phospholipids and kinetic deuterium isotope effects on CYP1B1 substrate oxidation

The effect of phospholipids on the kinetic parameters of three substrates, 7-ethoxy-4 -(trifluoromethyl)coumarin (7-EFC), 7-ethoxycoumarin (7-EC) and 17β-estradiol (E(2)), of human CYP1B1 was studied. In general, anionic phospholipids, phosphatidic acid and cardiolipin increased catalytic efficiency by increasing k(cat) values or decreasing K(m) values. The advantages of using the 7-EFC as a substrate over 7-EC and E(2) include high k(cat), low K(m) and high catalytic efficiency. Spectral binding titrations indicated that the binding affinity of 7-EFC to CYP1B1 in the presence or absence of phospholipids is higher than that of 7-EC or E(2). Furthermore, phosphatidylcholine increased the binding affinity of the substrates to the CYP1B1. High non-competitive intermolecular kinetic deuterium isotope effects (values 5.4-12) were observed for O-deethylation of 7-EFC and 7-EC with deuterium substitution at the ethoxy group, indicating that the C-H bond-breaking step makes a major contribution to the rate of these CYP1B1-catalyzed reactions. However, the intermolecular kinetic deuterium isotope effect is ~2 for the E(2) 4-hydroxylation reaction, indicating that the C-H bond-breaking step contributes only partially to the rate of this CYP1B1-catalyzed reaction. These results indicate that the reaction mechanism of CYP1B1-catalyzed reactions is distinct for each substrate.

Brain chemistry: how does P450 catalyze the O-demethylation reaction of 5-methoxytryptamine to yield serotonin?

Density functional theory has been applied to elucidate the mechanism of the O-demethylation reaction that generates serotonin from 5-methoxytryptamine (5-MT); a process that is efficiently catalyzed by P450 CYP2D6. Two substrates, the neutral 5-MT and the protonated 5-MTH(+), were used to probe the reactivity of CYP2D6 compound I. Notably, the H-abstraction process is found to be slightly more facile for 5-MT. However, our DFT augmented by docking results show that the amino acid Glu216 in the active site holds the NH(3)(+) tail of the 5-MTH(+) substrate in an upright conformation and thereby controls the regioselectivity of the bond activation. Thus, the substrate protonation serves an important function in maximizing the yield of serotonin. This finding is in accord with experimental conclusions that 5-MTH(+) serves as the substrate for the CYP2D6 enzyme. The study further shows that the H-abstraction follows two-state reactivity (TSR), whereas the rebound path may involve more states due to the appearance of both Fe(IV) and Fe(III) electromers during the reaction of 5-MTH(+).

Fast prediction of cytochrome P450 mediated drug metabolism

Cytochrome P450 mediated metabolism of drugs is one of the major determinants of their kinetic profile, and prediction of this metabolism is therefore highly relevant during the drug discovery and development process. A new rule-based method, based on results from density functional theory calculations, for predicting activation energies for aliphatic and aromatic oxidations by cytochromes P450 is developed and compared with several other methods. Although the applicability of the method is currently limited to a subset of P450 reactions, these reactions describe more than 90 % of the metabolites. The rules employed are relatively few and general, and when combined with solvent-accessible surface area calculations to account for steric accessibility, the method gives a major P450 metabolite as first-ranked position for 75 % of the substrates, and ranked in the top three for 90 % of substrates for a set of 20 substrates. In combination with docking, it can predict isoform-specific metabolism, and we apply this on CYP1A2 with very good results on 81 substrates, for which we find a major metabolite ranked in the top three for 90 % of the substrates (100 % in the training set and 87 % in the larger test set).

Heterologous expression and characterization of wild-type human cytochrome P450 1A2 without conventional N-terminal modification in Escherichia coli

In this study, wild-type human CYP1A2 without the conventional N-terminal modification (second codon GCT) or the truncation of the N-terminal hydrophobic region was functionally expressed in Escherichia coli. Its enzymatic properties were compared with N-terminally modified CYP1A2. Although modified CYP1A2 is almost all high-spin, some wild-type CYP1A2 shifted to low-spin. Spectral binding titrations with several ligands could be performed with wild-type enzyme, but not with modified enzyme. Kinetic parameters for several substrates were similar for the two CYP1A2 enzymes. However, the oxidation rates of phenacetin by modified enzyme were approximately 2-fold higher than those by wild-type enzyme. The intermolecular isotope effects were approximately 2 for phenacetin O-deethylation catalyzed by both enzymes. However, the wild-type enzyme, but not the modified enzyme, increased C-hydroxylation when O-deethylation rates were lowered by deuterium substitution. Molecular switching indicates that phenacetin rotates within the active site of wild-type enzyme and suggests a looser conformation in the active site of the wild-type enzyme than of the modified enzyme. These results reveal that the overall enzymatic properties of wild-type CYP1A2 enzyme are quite similar to those of modified CYP1A2, although its active site environment seems to differ from that of the modified enzyme.