Insight into the redox partner interaction mechanism in cytochrome P450BM-3 using molecular dynamics simulations

Mechanisms of Electron Transfer Rate Modulations in Cytochrome P450 BM3

Bacterial cytochromes P450 BM3 (CYP450 BM3) catalyze reactions of industrial importance. Despite many successful biotransformations, robust (re)design for novel applications remains challenging. Rational design and evolutionary approaches are not always successful, highlighting a lack of complete understanding of the mechanisms of electron transfer (ET) modulations. Thus, the full potential of CYP450 reactions remains under-exploited. In this work, we report the first molecular dynamics (MD)-based explicit prediction of BM3 ET parameters (reorganization energies; λ and ET free energies; ΔG°), and log ET rates (log k_ET) using the Marcus theory. Overall, the calculated ET rates for the BM3 wild-type (WT), mutants (F393 and L86), ligand-bound state, and ion concentrations agree well with experimental data. In ligand-free (LF) BM3, mutations modulate k_ET via ET ΔG°. Simulations show that the experimental ET rate enhancement is due to increased driving force (more negative ΔG°) upon ligation. This increase is related to the protein reorganization required to accommodate the ligand in the binding pocket rather than binding interactions with the ligand. Our methodology (CYPWare 1.0) automates all the stages of the MD simulation step-up, energy calculations, and estimation of ET parameters. CYPWare 1.0 and this work thus represent an important advancement in the CYP450 ET rate predictions, which has the potential to guide the redesign of ET enzymes. This program and a Web tool are available on GitHub for academic research.

Advanced 96-microtiter plate based bioelectrochemical platform reveals molecular short cut of electron flow in cytochrome P450 enzyme

In bioelectrocatalysis, immobilised redox enzymes are activated in a bioelectronic interface without redox equivalents such as NADPH, thus enabling heterogeneous flow chemistry. The functional contact between enzyme and electrode requires a high degree of optimisation regarding choice of electrode material, electrode pre-treatment, enzyme immobilisation and reaction conditions. So far, however, there are no systems that can easily enable an optimisation procedure at a higher throughput. Here, we present an advanced platform with a vertical divided cell architecture in conjunction with a developed 96-multipotentiostat to be able to drive redox enzymes in 96 well microtiter plate based multielectrode arrays. This platform controls 96 independent three-electrode setups with arbitrary working electrode materials. We demonstrate its applicability in a mutation study of cytochrome P450 BM3 using indium tin oxide as electrode material and the 7-ethoxycoumarin product quantification assay. We show that the bioelectrocatalytic activity of P450 BM3 can be amplified when the cofactor FAD is erased from the enzyme by a single point mutation, so that FMN becomes the first electron entry point. Bioelectrocatalysis thus offers an approach to enzyme simplification as a remedy for the inherent instability of self-sufficient cytochrome P450 enzymes. In addition, we examined native and artificial enzyme activation with respect to ionic strength and buffer composition. The optimal conditions of the activation types differ substantially from each other and exhibit a new molecular facet in enzyme characteristics. In a proof-of-principle we demonstrate that the platform is also compatible with raw cell extracts, thus opening the door for random mutagenesis screenings.

Characterization of the structure and interactions of P450 BM3 using hybrid mass spectrometry approaches

The cytochrome P450 monooxygenase P450 BM3 (BM3) is a biotechnologically important and versatile enzyme capable of producing important compounds such as the medical drugs pravastatin and artemether, and the steroid hormone testosterone. BM3 is a natural fusion enzyme comprising two major domains: a cytochrome P450 (heme-binding) catalytic domain and a NADPH-cytochrome P450 reductase (CPR) domain containing FAD and FMN cofactors in distinct domains of the CPR. A crystal structure of full-length BM3 enzyme is not available in its monomeric or catalytically active dimeric state. In this study, we provide detailed insights into the protein-protein interactions that occur between domains in the BM3 enzyme and characterize molecular interactions within the BM3 dimer by using several hybrid mass spectrometry (MS) techniques, namely native ion mobility MS (IM-MS), collision-induced unfolding (CIU), and hydrogen-deuterium exchange MS (HDX-MS). These methods enable us to probe the structure, stoichiometry, and domain interactions in the ∼240 kDa BM3 dimeric complex. We obtained high-sequence coverage (88–99%) in the HDX-MS experiments for full-length BM3 and its component domains in both the ligand-free and ligand-bound states. We identified important protein interaction sites, in addition to sites corresponding to heme-CPR domain interactions at the dimeric interface. These findings bring us closer to understanding the structure and catalytic mechanism of P450 BM3.

Engineering Chemoselectivity in Hemoprotein-Catalyzed Indole Amidation

Here we report a cytochrome P450 variant that catalyzes C₂-amidation of 1-methylindoles with tosyl azide via nitrene transfer. Before evolutionary optimization the enzyme exhibited two undesired side reactivities resulting in reduction of the putative iron-nitrenoid intermediate or cycloaddition between the two substrates to form triazole products. We speculated that triazole formation was a promiscuous cycloaddition activity of the P450 heme domain, while sulfonamide formation likely arose from surplus electron transfer from the reductase domain. Directed evolution involving mutagenesis of both the heme and reductase domains delivered an enzyme providing the desired indole amidation products with up to 8400 turnovers, 90% yield, and a shift in chemoselectivity from 2:19:1 to 110:12:1 in favor of nitrene transfer over reduction or triazole formation. This work expands the substrate scope of hemoprotein nitrene transferases to heterocycles and highlights the adaptability of the P450 scaffold to solve challenging chemoselectivity problems in non-natural enzymatic catalysis.

Choreography of the Reductase and P450BM3 Domains Toward Electron Transfer Is Instigated by the Substrate

The driving force for the electron transfer (ET) step in the catalytic cycle of cytochrome P450_BM3 is investigated using three sets of 1 μs molecular dynamic simulations for the resting state of P450 in complex with the flavin (FMN) in the reductase domain. These sets involve the following: (i) substrate-free (SF), (ii) substrate (N-palmitoyl glycine, i.e., NPG)-bound (SB), and (iii) SB with the semiquinone radical anion (SQ^-) of FMN. Starting from the X-ray structure of the SF heme domain, we observe that the α1-helix (of the reductase) and the C-helix (of the heme) undergo reorientation to a parallel orientation, which is the thermochemically stable form. The reorientation of the helices pushes away the FMN to a distance of 18.4 Å from the heme's center. When the substrate binds it causes the I-helix of the heme domain to kink and push the C-helix toward the α1-helix, thereby locking the latter two into a stabilized perpendicular conformation, wherein the FMN-heme distance is 12 Å. The distance drops further in the SQ^- form, and upon QM/MM geometry optimization the two moieties approach 8.8 Å, which enhances the ET rate (by 10⁴-10⁶ fold) to the heme's Fe³⁺ ion. These motions are driven by hydrogen bond strengthening between the C- and the α1-helices. Finally, substrate binding leads to formation of an organized water chain connecting the FMN and heme moieties. The water channel assists the ET and couples it to the proton transfer steps that should activate O₂ and create the oxo-iron active species.

Modulating proposed electron transfer pathways in P450BM3 led to improved activity and coupling efficiency

Electrochemical in vitro reduction of P450 enzymes is a promising alternative to in vivo applications. Previously we presented three engineered P450_BM3 variants for aniline hydroxylation, equipped with a carbon nanotube binding-peptide (CNT-tag) for self-assembly on CNT electrodes. Compared to wildtype P450_BM3 the NADPH-dependent activity was enhanced, but the coupling efficiency remained low. For P450_BM3 Verma, Schwaneberg and Roccatano (2014, Biopolymers 101, 197-209) calculated putative electron transfer pathways (eTPs) by MD simulations. We hypothesised that knockouts of these transfer pathways would alter the coupling efficiency of the system. The results revealed no improved system for the electrically-driven P450s. For the NADPH-driven P450s, however, the most active eTP-mutant showed a 13-fold increased activity and a 32-fold elevated coupling efficiency using NADPH as reducing equivalent. This suggests an alternative principle of electron transport for the reduction by NADPH and an electrode, respectively. The work presents moreover a tool to improve the coupling and activity of P450s with non-natural substrates.

Crystallographic insights into a cobalt (III) sepulchrate based alternative cofactor system of P450 BM3 monooxygenase

P450 BM3 is a multi-domain heme-containing soluble bacterial monooxygenase. P450 BM3 and variants are known to oxidize structurally diverse substrates. Crystal structures of individual domains of P450 BM3 are available. However, the spatial organization of the full-length protein is unknown. In this study, crystal structures of the P450 BM3 M7 heme domain variant with and without cobalt (III) sepulchrate are reported. Cobalt (III) sepulchrate acts as an electron shuttle in an alternative cofactor system employing zinc dust as the electron source. The crystal structure shows a binding site for the mediator cobalt (III) sepulchrate at the entrance of the substrate access channel. The mediator occupies an unusual position which is far from the active site and distinct from the binding of the natural redox partner (FAD/NADPH binding domain).

In Silico Studies of Small Molecule Interactions with Enzymes Reveal Aspects of Catalytic Function

Small molecules, such as solvent, substrate, and cofactor molecules, are key players in enzyme catalysis. Computational methods are powerful tools for exploring the dynamics and thermodynamics of these small molecules as they participate in or contribute to enzymatic processes. In-depth knowledge of how small molecule interactions and dynamics influence protein conformational dynamics and function is critical for progress in the field of enzyme catalysis. Although numerous computational studies have focused on enzyme-substrate complexes to gain insight into catalytic mechanisms, transition states and reaction rates, the dynamics of solvents, substrates, and cofactors are generally less well studied. Also, solvent dynamics within the biomolecular solvation layer play an important part in enzyme catalysis, but a full understanding of its role is hampered by its complexity. Moreover, passive substrate transport has been identified in certain enzymes, and the underlying principles of molecular recognition are an area of active investigation. Enzymes are highly dynamic entities that undergo different conformational changes, which range from side chain rearrangement of a residue to larger-scale conformational dynamics involving domains. These events may happen nearby or far away from the catalytic site, and may occur on different time scales, yet many are related to biological and catalytic function. Computational studies, primarily molecular dynamics (MD) simulations, provide atomistic-level insight and site-specific information on small molecule interactions, and their role in conformational pre-reorganization and dynamics in enzyme catalysis. The review is focused on MD simulation studies of small molecule interactions and dynamics to characterize and comprehend protein dynamics and function in catalyzed reactions. Experimental and theoretical methods available to complement and expand insight from MD simulations are discussed briefly.

Insights into an efficient light-driven hybrid P450 BM3 enzyme from crystallographic, spectroscopic and biochemical studies

BACKGROUND

In order to perform selective CH functionalization upon visible light irradiation, Ru(II)-diimine functionalized P450 heme enzymes have been developed. The sL407C-1 enzyme containing the Ru(bpy)₂PhenA (bpy=2,2'-bipyridine and PhenA=5-acetamido-1,10-phenanthroline) photosensitizer (1) covalently attached to the non-native single cysteine L407C of the P450BM3 heme domain mutant, displays high photocatalytic activity in the selective CH bond hydroxylation of several substrates.

METHODS

A combination of X-ray crystallography, site-directed mutagenesis, transient absorption measurements and enzymatic assays was used to gain insights into its photocatalytic activity and electron transfer pathway.

RESULTS

The crystal structure of the sL407C-1 enzyme was solved in the open and closed conformations revealing a through-space electron transfer pathway involving highly conserved, F393 and Q403, residues. Several mutations of these residues (F393A, F393W or Q403W) were introduced to probe their roles in the overall reaction. Transient absorption measurements confirm rapid electron transfer as heme reduction is observed in all four hybrid enzymes. Compared to the parent sL407C-1, photocatalytic activity was negligible in the dF393A-1 enzyme while 60% increase in activity with total turnover numbers of 420 and 90% product conversion was observed with the dQ403W-1 mutant.

CONCLUSIONS

In the sL407C-1 enzyme, the photosensitizer is ideally located to rapidly deliver electrons, using the naturally occurring electron transfer pathway, to the heme center in order to activate molecular dioxygen and sustain photocatalytic activity.

GENERAL SIGNIFICANCE

The results shed light on the design of efficient light-driven biocatalysts and the approach can be generalized to other members of the P450 superfamily.

Versatile biocatalysis of fungal cytochrome P450 monooxygenases

Cytochrome P450 (CYP) monooxygenases, the nature’s most versatile biological catalysts have unique ability to catalyse regio-, chemo-, and stereospecific oxidation of a wide range of substrates under mild reaction conditions, thereby addressing a significant challenge in chemocatalysis. Though CYP enzymes are ubiquitous in all biological kingdoms, the divergence of CYPs in fungal kingdom is manifold. The CYP enzymes play pivotal roles in various fungal metabolisms starting from housekeeping biochemical reactions, detoxification of chemicals, and adaptation to hostile surroundings. Considering the versatile catalytic potentials, fungal CYPs has gained wide range of attraction among researchers and various remarkable strategies have been accomplished to enhance their biocatalytic properties. Numerous fungal CYPs with multispecialty features have been identified and the number of characterized fungal CYPs is constantly increasing. Literature reveals ample reviews on mammalian, plant and bacterial CYPs, however, modest reports on fungal CYPs urges a comprehensive review highlighting their novel catalytic potentials and functional significances. In this review, we focus on the diversification and functional diversity of fungal CYPs and recapitulate their unique and versatile biocatalytic properties. As such, this review emphasizes the crucial issues of fungal CYP systems, and the factors influencing efficient biocatalysis.

Structure, dynamics, and function of the monooxygenase P450 BM-3: insights from computer simulations studies

The monooxygenase P450 BM-3 is a NADPH-dependent fatty acid hydroxylase enzyme isolated from soil bacterium Bacillus megaterium. As a pivotal member of cytochrome P450 superfamily, it has been intensely studied for the comprehension of structure-dynamics-function relationships in this class of enzymes. In addition, due to its peculiar properties, it is also a promising enzyme for biochemical and biomedical applications. However, despite the efforts, the full understanding of the enzyme structure and dynamics is not yet achieved. Computational studies, particularly molecular dynamics (MD) simulations, have importantly contributed to this endeavor by providing new insights at an atomic level regarding the correlations between structure, dynamics, and function of the protein. This topical review summarizes computational studies based on MD simulations of the cytochrome P450 BM-3 and gives an outlook on future directions.

Cytochromes P450 as promising catalysts for biotechnological application: chances and limitations

Cytochromes P450 (CYPs) belong to the superfamily of heme b containing monooxygenases with currently more than 21,000 members. These enzymes accept a vast range of organic molecules and catalyze diverse reactions. These extraordinary capabilities of CYP systems that are unmet by other enzymes make them attractive for biotechnology. However, the complexity of these systems due to the need of electron transfer from nicotinamide adenine dinucleotide (phosphate) (NAD(P)H) via redox partner proteins for the initial hydroxylation step limits a broader technical implementation of CYP enzymes. There have been several reviews during the past years tackling the potential CYPs for synthetic application. The aim of this review is to give a critical overview about possibilities and chances for application of these interesting catalysts as well as to discuss drawbacks and problems related to their use. Solutions to overcome these limitations will be demonstrated, and several selected examples of successful CYP applications under industrial conditions will be reviewed.