Peroxidase activity stabilization of cytochrome P450(BM3) by rational analysis of intramolecular electron transfer

Combined quantum mechanical and molecular mechanical (QM/MM) calculations were used to explore the electron pathway involved in the suicide inactivation of cytochrome P450BM3 from Bacillus megaterium. The suicide inactivation is a common phenomenon observed for heme peroxidases, in which the enzyme is inactivated as a result of self-oxidation mediated by highly oxidizing enzyme intermediates formed during the catalytic cycle. The selected model was a mutant comprising only the heme domain (CYPBM3 21B3) that had been previously evolved to efficiently catalyze hydroxylation reactions with hydrogen peroxide (H2O2) as electron acceptor. An extensive mapping of residues involved in electron transfer routes was obtained from density functional calculations on activated heme (i.e. Compound I) and selected amino acid residues. Identification of oxidizable residues (electron donors) was performed by selectively activating/deactivating different quantum regions. This method allowed a rational identification of key oxidizable targets in order to replace them for less oxidizable residues by site-directed mutagenesis. The residues W96 and F405 were consistently predicted by the QM/MM electron pathway to hold high spin density; single and double mutants of P450BM3 on these positions (W96A, F405L, W96A/F405L) resulted in a more stable variants in the presence of hydrogen peroxide, displaying a similar reaction rate than P450BM3 21B3. Furthermore, mass spectrometry confirmed these oxidation sites and corroborated the possible routes described by QM/MM electron transfer (ET) pathways.