Electronic Structure of Compound I in Human Isoforms of Cytochrome P450 from QM/MM Modeling
- 26 August 2005
- journal article
- Published by American Chemical Society (ACS) in Journal of the American Chemical Society
- Vol. 127 (37) , 12900-12908
- https://doi.org/10.1021/ja0520924
Abstract
Human cytochromes P450 play a vital role in drug metabolism. The key step in substrate oxidation involves hydrogen atom abstraction or C=C bond addition by the oxygen atom of the Compound I intermediate. The latter has three unpaired electrons, two on the Fe-O center and one shared between the porphyrin ring and the proximal cysteinyl sulfur atom. Changes in its electronic structure have been suggested to affect reactivity. The electronic and geometric structure of Compound I in three important human subfamilies of cytochrome P450 (P450, 2C, 2B, and 3A) that are major contributors to drug metabolism is characterized here using combined quantum mechanical/molecular mechanical (QM/MM) calculations at the B3LYP:CHARMM27 level. Compound I is remarkably similar in all isoforms, with the third unpaired electron located mainly on the porphyrin ring, and this prediction is not very sensitive to details of the QM/MM methodology, such as the DFT functional, the basis set, or the size of the QM region. The presence of substrate also has no effect. The main source of variability in spin density on the cysteinyl sulfur (from 26 to 50%) is the details of the system setup, such as the starting protein geometry used for QM/MM minimization. This conformational effect is larger than the differences between human isoforms, which are therefore not distinguishable on electronic grounds, so it is unlikely that observed large differences in substrate selectivity can be explained to a large extent in these terms.Keywords
This publication has 50 references indexed in Scilit:
- Prediction of binding modes for ligands in the cytochromes P450 and other heme‐containing proteinsProteins-Structure Function and Bioinformatics, 2005
- Aromatic Hydroxylation by Cytochrome P450: Model Calculations of Mechanism and Substituent EffectsJournal of the American Chemical Society, 2003
- Regioselectivity of CYP2B6: homology modeling, molecular dynamics simulation, dockingJournal of Molecular Modeling, 2002
- The High-Valent Compound of Cytochrome P450: The Nature of the Fe−S Bond and the Role of the Thiolate Ligand as an Internal Electron DonorPublished by Wiley ,2000
- All-Atom Empirical Potential for Molecular Modeling and Dynamics Studies of ProteinsThe Journal of Physical Chemistry B, 1998
- Delocalization over the heme and the axial ligands of one of the two oxidizing equivalents stored above the ferric state in the peroxidase and catalase Compound-I intermediates: indirect participation of the proximal axial ligand of iron in the oxidation reactions catalyzed by heme-based peroxidases and catalases?JBIC Journal of Biological Inorganic Chemistry, 1996
- Optimization of the Lennard-Jones parameters for a combinedab initio quantum mechanical and molecular mechanical potential using the 3-21G basis setJournal of Computational Chemistry, 1996
- Heme-Containing OxygenasesChemical Reviews, 1996
- Evidence for Compound I Formation in the Reaction of Cytochrome-P450cam with m-Chloroperbenzoic AcidBiochemical and Biophysical Research Communications, 1994
- CHARMM: A program for macromolecular energy, minimization, and dynamics calculationsJournal of Computational Chemistry, 1983