Theoretical study of the ligand-CYP2B4 complexes: effect of structure on binding free energies and heme spin state.

The molecular origins of temperature-dependent ligand-binding affinities and ligand-induced heme spin state conversion have been investigated using free energy analysis and DFT calculations for substrates and inhibitors of cytochrome P450 2B4 (CYP2B4), employing models of CYP2B4 based on CYP2C5(3LVdH)/CYP2C9 crystal structures, and the results compared with experiment. DFT calculations indicate that large heme-ligand interactions (ca. -15 kcal/mol) are required for inducing a high to low spin heme transition, which is correlated with large molecular electrostatic potentials (approximately -45 kcal/mol) at the ligand heteroatom. While type II ligands often contain oxygen and nitrogen heteroatoms that ligate heme iron, DFT results indicate that BP and MF heme complexes, with weak substrate-heme interactions (ca. -2 kcal/mol), and modest MEPS minima (>-35 kcal/mol) are high spin. In contrast, heme complexes of the CYP2B4 inhibitor, 4PI, the product of benzphetamine metabolism, DMBP, and water are low spin, have substantial heme-ligand interaction energies (<-15 kcal/mol) and deep MEPS minima (<-45 kcal/mol) near their heteroatoms. MMPBSA analysis of MD trajectories were made to estimate binding free energies of these ligands at the heme binding site of CYP2B4. In order to initially assess the realism of this approach, the binding free energy of 4PI inhibitor was computed and found to be a reasonable agreement with experiment: -7.7 kcal/mol [-7.2 kcal/mol (experiment)]. BP was determined to be a good substrate [-6.3 kcal/mol (with heme-ligand water), -7.3 kcal/mol (without ligand water)/-5.8 kcal/mol (experiment)], whereas the binding of MF was negligible, with only marginal binding binding free energy of -1.7 kcal/mol with 2-MF bound [-3.8 kcal/mol (experiment)], both with and without retained heme-ligand water. Analysis of the free energy components reveal that hydrophobic/nonpolar contributions account for approximately 90% of the total binding free energy of these substrates and are the source of their differential and temperature-dependent CYP2B4 binding. The results indicate the underlying origins of the experimentally observed differential binding affinities of BP and MF, and indicate the plausibility of the use of models derived from moderate sequence identity templates in conjunction with approximate free energy methods in the estimation of ligand-P450 binding affinities.

Determination of the rate of reduction of oxyferrous cytochrome P450 2B4 by 5-deazariboflavin adenine dinucleotide T491V cytochrome P450 reductase.

The use of 5-deazaFAD T491V cytochrome P450 reductase has made it possible to directly measure the rate of electron transfer to microsomal oxyferrous cytochrome (cyt) P450 2B4. In this reductase the FMN moiety can be reduced to the hydroquinone, FMNH(2), while the 5-deazaFAD moiety remains oxidized [Zhang, H., et al. (2003) Biochemistry 42, 6804-6813]. The rate of electron transfer from 5-deazaFAD cyt P450 reductase to oxyferrous cyt P450 was determined by rapidly mixing the ferrous cyt P450-2-electron-reduced 5-deazaFAD T491V reductase complex with oxygen in the presence of substrate. The 5-deazaFAD T491V reductase which can only donate a single electron reduces the oxyferrous cyt P450 and oxidizes to the air-stable semiquinone, with rate constants of 8.4 and 0.37 s(-1) at 15 degrees C. Surprisingly, oxyferrous cyt P450 turns over more slowly with a rate constant of 0.09 s(-1), which is the rate of catalysis under steady-state conditions at 15 degrees C (k(cat) = 0.08 s(-1)). In contrast, the rate constant for electron transfer from ferrous cyt b(5) to oxyferrous cyt P450 is 10 s(-1) with oxyferrous cyt P450 and cyt b(5) simultaneously undergoing spectral changes. Quantitative analyses by LC-MS/MS revealed that the product, norbenzphetamine, was formed with a coupling efficiency of 52% with cyt b(5) and 32% with 5-deazaFAD T491V reductase. Collectively, these results suggest that during catalysis a relatively stable reduced oxyferrous intermediate of cyt P450 is formed in the presence of cyt P450 reductase but not cyt b(5) and that the rate-limiting step in catalysis follows introduction of the second electron.

Preparation and characterization of a 5'-deazaFAD T491V NADPH-cytochrome P450 reductase.

NADPH-cytochrome P450 reductase is a flavoprotein which contains both an FAD and FMN cofactor. Since the distribution of electrons is governed solely by the redox potentials of the cofactors, there are nine different ways the electrons can be distributed and hence nine possible unique forms of the protein. More than one species of reductase will exist at a given level of oxidation except when the protein is either totally reduced or oxidized. In an attempt to unambiguously characterize the redox properties of the physiologically relevant FMNH(2) form of the reductase, the T491V mutant of NADPH-cytochrome P450 reductase has been reconstituted with 5'-deazaFAD which binds to the FAD-binding site of the reductase with a K(d) of 94 nM. The 5'-deazaFAD cofactor does not undergo oxidation or reduction under our experimental conditions. The molar ratio of FMN to 5'-deazaFAD in the reconstituted reductase was 1.1. Residual FAD accounted for less than 5% of the total flavins. Addition of 2 electron equivalents to the 5'-deazaFAD T491V reductase from dithionite generated a stoichiometric amount of the FMN hydroquinone form of the protein. The 5'-deazaFAD moiety remained oxidized under these conditions due to its low redox potential (-650 mV). The 2-electron-reduced 5'-deazaFAD reductase was capable of transferring only a single electron from its FMN domain to its redox partners, ferric cytochrome c and cytochrome b(5). Reduction of the cytochromes and oxidation of the reductase occurred simultaneously. The FMNH(2) in the 5'-deazaFAD reductase autoxidizes with a first-order rate constant of 0.007 s(-)(1). Availability of a stable NADPH-cytochrome P450 reductase capable of donating only a single electron to its redox partners provides a unique tool for investigating the electron-transfer properties of an intact reductase molecule.