Rational engineering of cytochromes P450 2B6 and 2B11 for enhanced stability: Insights into structural importance of residue 334

Use of engineered cytochromes P450 for accelerating drug discovery and development

Numerous steps in drug development, including the generation of authentic metabolites and late-stage functionalization of candidates, necessitate the modification of often complex molecules, such as natural products. While it can be challenging to make the required regio- and stereoselective alterations to a molecule using purely chemical catalysis, enzymes can introduce changes to complex molecules with a high degree of stereo- and regioselectivity. Cytochrome P450 enzymes are biocatalysts of unequalled versatility, capable of regio- and stereoselective functionalization of unactivated CH bonds by monooxygenation. Collectively they catalyze over 60 different biotransformations on structurally and functionally diverse organic molecules, including natural products, drugs, steroids, organic acids and other lipophilic molecules. This catalytic versatility and substrate range makes them likely candidates for application as potential biocatalysts for industrial chemistry. However, several aspects of the P450 catalytic cycle and other characteristics have limited their implementation to date in industry, including: their lability at elevated temperature, in the presence of solvents, and over lengthy incubation times; the typically low efficiency with which they metabolize non-natural substrates; and their lack of specificity for a single metabolic pathway. Protein engineering by rational design or directed evolution provides a way to engineer P450s for industrial use. Here we review the progress made to date toward engineering the properties of P450s, especially eukaryotic forms, for industrial application, and including the recent expansion of their catalytic repertoire to include non-natural reactions.

Determinants of thermostability in the cytochrome P450 fold

Cytochromes P450 are found throughout the biosphere in a wide range of environments, serving a multitude of physiological functions. The ubiquity of the P450 fold suggests that it has been co-opted by evolution many times, and likely presents a useful compromise between structural stability and conformational flexibility. The diversity of substrates metabolized and reactions catalyzed by P450s makes them attractive starting materials for use as biocatalysts of commercially useful reactions. However, process conditions impose different requirements on enzymes to those in which they have evolved naturally. Most natural environments are relatively mild, and therefore most P450s have not been selected in Nature for the ability to withstand temperatures above ~40°C, yet industrial processes frequently require extended incubations at much higher temperatures. Thus, there has been considerable interest and effort invested in finding or engineering thermostable P450 systems. Numerous P450s have now been identified in thermophilic organisms and analysis of their structures provides information as to mechanisms by which the P450 fold can be stabilized. In addition, protein engineering, particularly by directed or artificial evolution, has revealed mutations that serve to stabilize particular mesophilic enzymes of interest. Here we review the current understanding of thermostability as it applies to the P450 fold, gleaned from the analysis of P450s characterized from thermophilic organisms and the parallel engineering of mesophilic forms for greater thermostability. We then present a perspective on how this information might be used to design stable P450 enzymes for industrial application. This article is part of a Special Issue entitled: Cytochrome P450 biodiversity and biotechnology, edited by Erika Plettner, Gianfranco Gilardi, Luet Wong, Vlada Urlacher, Jared Goldstone.

Functional importance of a peripheral pocket in mammalian cytochrome P450 2B enzymes

The functional importance of a peripheral pocket found in previously published X-ray crystal structures of CYP2B4 and CYP2B6 was probed using a biophysical approach. Introduction of tryptophan within the pocket of CYP2B4 at F202 or I241 leads to marked impairment of 7-ethoxy-4-(trifluoromethyl)coumarin (7-EFC) or 7-benzyloxyresorufin O-dealkylation efficiency; a similar substitution at F195, near the surface access to the pocket, does not affect these activities. The analogous CYP2B6 F202W mutant is inactive in the 7-EFC O-dealkylation assay. The stoichiometry of 7-EFC deethylation suggested that the decreased activity of F202W and I241W in CYP2B4 and lack of activity of F202W in CYP2B6 coincided with a sharp increase in the flux of reducing equivalents through the oxidase shunt to produce excess water. The results indicate that the chemical identity of residues within this peripheral pocket, but not at the mouth of the pocket, is important in substrate turnover and redox coupling, likely through effects on active site topology.

SNP genetic polymorphisms of MDR-1, CYP1A2 and CYPB11 genes in four canine breeds upon toxicological evaluation

The fields of pharmacogenetics and pharmacogenomics have become increasingly promising regarding the clinical application of genetic data to aid in prevention of adverse reactions. Specific screening tests can predict which animals express modified proteins or genetic sequences responsible for adverse effects associated with a drug. Among the genetic variations that have been investigated in dogs, the multidrug resistance gene (MDR) is the best studied. However, other genes such as CYP1A2 and CYP2B11 control the protein syntheses involved in the metabolism of many drugs. In the present study, the MDR-1, CYP1A2 and CYP2B11 genes were examined to identify SNP polymorphisms associated with these genes in the following four canine breeds: Uruguayan Cimarron, Border Collie, Labrador Retriever and German Shepherd. The results revealed that several SNPs of the CYP1A2 and CYP2B11 genes are potential targets for drug sensitivity investigations.

Functional characterization of cytochromes P450 2B from the desert woodrat Neotoma lepida

Mammalian detoxification processes have been the focus of intense research, but little is known about how wild herbivores process plant secondary compounds, many of which have medicinal value or are drugs. cDNA sequences that code for three enzymes of the cytochrome P450 (CYP) 2B subfamily, here termed 2B35, 2B36, and 2B37 have been recently identified from a wild rodent, the desert woodrat (Malenke et al., 2012). Two variant clones of each enzyme were engineered to increase protein solubility and to facilitate purification, as reported for CYP2B enzymes from multiple species. When expressed in Escherichia coli each of the woodrat proteins gave the characteristic maximum at 450nm in a reduced carbon monoxide difference spectrum but generally expressed at lower levels than rat CYP2B1. Two enzymes, 2B36 and 2B37, showed dealkylation activity with the model substrates 7-ethoxy-4-(trifluoromethyl)coumarin and 7-benzyloxyresorufin, whereas 2B35 was inactive. Binding of the monoterpene (+)-α-pinene produced a Type I shift in the absorbance spectrum of each enzyme. Mutation of 2B37 at residues 114, 262, or 480, key residues governing ligand interactions with other CYP2B enzymes, did not significantly change expression levels or produce the expected functional changes. In summary, two catalytic and one ligand-binding assay are sufficient to distinguish among CYP2B35, 2B36, and 2B37. Differences in functional profiles between 2B36 and 2B37 are partially explained by changes in substrate recognition site residue 114, but not 480. The results advance our understanding of the mechanisms of detoxification in wild mammalian herbivores and highlight the complexity of this system.

Multidimensional epistasis and fitness landscapes in enzyme evolution

The conventional analysis of enzyme evolution is to regard one single salient feature as a measure of fitness, expressed in a milieu exposing the possible selective advantage at a given time and location. Given that a single protein may serve more than one function, fitness should be assessed in several dimensions. In the present study we have explored individual mutational steps leading to a triple-point-mutated human GST (glutathione transferase) A2-2 displaying enhanced activity with azathioprine. A total of eight alternative substrates were used to monitor the diverse evolutionary trajectories. The epistatic effects of the mutations on catalytic activity were variable in sign and magnitude and depended on the substrate used, showing that epistasis is a multidimensional quality. Evidently, the multidimensional fitness landscape can lead to alternative trajectories resulting in enzymes optimized for features other than the selectable markers relevant at the origin of the evolutionary process. In this manner the evolutionary response is robust and can adapt to changing environmental conditions.

Plasticity of CYP2B enzymes: structural and solution biophysical methods

In the past three years, major advances in understanding cytochrome P450 2B (CYP2B) structure-function relationships have been made through determination of multiple ligand-bound and one ligand-free X-ray crystal structure of CYP2B4 and one ligand-bound X-ray crystal structure of CYP2B6. These structures have provided insight into the features that provide the high degree of plasticity of the enzymes. A combination of a phenylalanine cluster that allows for concerted movement of helices F through G and a conserved set of electrostatic interactions involving Arg(262) facilitates movement of this region to accommodate binding of ligands of various sizes without perturbing most of the P450 fold. Integrating solution based techniques such as NMR or deuterium exchange mass spectrometry (DXMS) with computational methods including molecular docking has provided further insight into enzyme behavior upon ligand binding. In addition, extended molecular dynamics simulations have provided a link between an open and a closed conformation of ligand-free CYP2B4 found in crystal structures. Other studies revealed the utility of rational engineering in improving stability of P450s to facilitate structural studies. The solution and computational results combined with the X-ray crystal structures yield a comprehensive picture of how these enzymes adopt different conformations to bind various ligands.

Investigation by site-directed mutagenesis of the role of cytochrome P450 2B4 non-active-site residues in protein-ligand interactions based on crystal structures of the ligand-bound enzyme

Residues located outside the active site of cytochromes P450 2B have exhibited importance in ligand binding, structural stability and drug metabolism. However, contributions of non-active-site residues to the plasticity of these enzymes are not known. Thus, a systematic investigation was undertaken of unique residue-residue interactions found in crystal structures of P450 2B4 in complex with 4-(4-chlorophenyl)imidazole (4-CPI), a closed conformation, or in complex with bifonazole, an expanded conformation. Nineteen mutants distributed over 11 sites were constructed, expressed in Escherichia coli and purified. Most mutants showed significantly decreased expression, especially in the case of interactions found in the 4-CPI structure. Six mutants (H172A, H172F, H172Q, L437A, E474D and E474Q) were chosen for detailed functional analysis. Among these, the K(s) of H172F for bifonazole was ∼ 20 times higher than for wild-type 2B4, and the K(s) of L437A for 4-CPI was ∼ 50 times higher than for wild-type, leading to significantly altered inhibitor selectivity. Enzyme function was tested with the substrates 7-ethoxy-4-(trifluoromethyl)coumarin, 7-methoxy-4-(trifluoromethyl)coumarin and 7-benzyloxyresorufin (7-BR). H172F was inactive with all three substrates, and L437A did not turn over 7-BR. Furthermore, H172A, H172Q, E474D and E474Q showed large changes in k(cat)/K(M) for each of the three substrates, in some cases up to 50-fold. Concurrent molecular dynamics simulations yielded distances between some of the residues in these putative interaction pairs that are not consistent with contact. The results indicate that small changes in the protein scaffold lead to large differences in solution behavior and enzyme function.

Analysis of Cytochrome P450 Conserved Sequence Motifs between Helices E and H: Prediction of Critical Motifs and Residues in Enzyme Functions

Rational approaches have been extensively used to investigate the role of active site residues in cytochrome P450 (CYP) functions. However, recent studies using random mutagenesis suggest an important role for non-active site residues in CYP functions. Meta-analysis of the random mutants showed that 75% of the functionally important non-active site residues are present in 20% of the entire protein between helices E and H (E-H) and conserved sequence motif (CSM) between 7 and 11. The CSM approach was developed recently to investigate the functional role of non-active site residues in CYP2B4. Furthermore, we identified and analyzed the CSM in multiple CYP families and subfamilies in the E-H region. Results from CSM analysis showed that CSM 7, 8, 10, and 11 are conserved in CYP1, CYP2, and CYP3 families, while CSM 9 is conserved only in CYP2 family. Analysis of different CYP2 subfamilies showed that CYP2B and CYP2C have similar characteristics in the CSM, while the characteristics of CYP2A and CYP2D subfamilies are different. Finally, we analyzed CSM 7, 8, 10, and 11, which are common in all the CYP families/subfamilies analyzed, in fifteen important drug-metabolizing CYPs. The results showed that while CSM 8 is most conserved among these CYPs, CSM 7, 9, and 10 have significant variations. We suggest that CSM8 has a common role in all the CYPs that have been analyzed, while CSM 7, 10, and 11 may have relatively specific role within the subfamily. We further suggest that these CSM play important role in opening and closing of the substrate access/egress channel by modulating the flexible/plastic region of the protein. Thus, site-directed mutagenesis of these CSM can be used to study structure-function and dynamic/plasticity-function relationships and to design CYP biocatalysts.

Comparison of in vitro metabolism of ticlopidine by human cytochrome P450 2B6 and rabbit cytochrome P450 2B4

A recent X-ray crystal structure of a rabbit cytochrome P450 2B4 (CYP2B4)-ticlopidine complex indicated that the compound could be modeled with either the thiophene or chlorophenyl group oriented toward the heme prosthetic group. Subsequent NMR relaxation and molecular docking studies suggested that orientation with the chlorophenyl ring closer to the heme was the preferred one. To evaluate the predictive value of these findings, the oxidation of ticlopidine by reconstituted CYP2B4 was studied and compared with CYP2B6, in which the thiophene portion of the molecule likely orients toward the heme. In vitro incubation of ticlopidine with both enzymes yielded the same set of metabolites: 7-hydroxyticlopidine (M1), 2-oxoticlopidine (M2), 5-(2-chlorobenzyl)thieno[3,2-c]pyridin-5-ium metabolite (M3), 5-(2-chlorobenzyl)thieno[3,2-c]pyridin-5-ium metabolite (M4), ticlopidine N-oxide (M5), and ticlopidine S-oxide dimer, a dimerization product of ticlopidine S-oxide (M6). The rates of metabolite formation deviated markedly from linearity with time, consistent with the known inactivation of CYP2B6 by ticlopidine. Fitting to a first-order equation yielded similar rate constants (k(obs)) for both enzymes. However, the amplitude (R(max)) of M1 and M6 formation was 4 to 5 times higher for CYP2B6 than CYP2B4, indicating a greater residence time of ticlopidine with its thiophene ring closer to heme in CYP2B6. In contrast, CYP2B4 formed M4 and M5 in more abundance than CYP2B6, indicating an alternate orientation. Overall, the results suggest that the preferential orientation of ticlopidine in the active site of CYP2B4 predicted by X-ray crystallography and NMR studies is unproductive and that ticlopidine likely reorients within CYP2B4 to a more productive mode.

Constrained water access to the active site of cytochrome P450 from the piezophilic bacterium Photobacterium profundum

Living species inhabiting ocean deeps must adapt to high hydrostatic pressure. This adaptation, which must enable functioning under conditions of promoted protein hydration, is especially important for proteins such as cytochromes P450 that exhibit functionally important hydration-dehydration dynamics. Here we study the interactions of substrates with cytochrome P450-SS9, a putative fatty acid hydroxylase from the piezophilic bacterium Photobacterium profundum SS9, and characterize the protein's barotropic properties. Comparison of P450-SS9 with cytochrome P450BM-3, a mesophilic fatty acid hydroxylase, suggests that P450-SS9 is characterized by severely confined accessibility and low water occupancy of the active site. This feature may reveal a mechanism of structural adaptation of the piezophilic enzyme. We also demonstrate that saturated and unsaturated fatty acids exert opposite effects on solvent accessibility and hydration of the active site. Modulation of the protein conformation by fatty acids is hypothesized to have an important physiological function in the piezophile.

Enhanced thermostability of methyl parathion hydrolase from Ochrobactrum sp. M231 by rational engineering of a glycine to proline mutation

Protein thermostability can be increased by some glycine to proline mutations in a target protein. However, not all glycine to proline mutations can improve protein thermostability, and this method is suitable only at carefully selected mutation sites that can accommodate structural stabilization. In this study, homology modeling and molecular dynamics simulations were used to select appropriate glycine to proline mutations to improve protein thermostability, and the effect of the selected mutations was proved by the experiments. The structure of methyl parathion hydrolase (MPH) from Ochrobactrum sp. M231 (Ochr-MPH) was constructed by homology modeling, and molecular dynamics simulations were performed on the modeled structure. A profile of the root mean square fluctuations of Ochr-MPH was calculated at the nanosecond timescale, and an eight-amino acid loop region (residues 186-193) was identified as having high conformational fluctuation. The two glycines nearest to this region were selected as mutation targets that might affect protein flexibility in the vicinity. The structures and conformational fluctuations of two single mutants (G194P and G198P) and one double mutant (G194P/G198P) were modeled and analyzed using molecular dynamics simulations. The results predicted that the mutant G194P had the decreased conformational fluctuation in the loop region and might increase the thermostability of Ochr-MPH. The thermostability and kinetic behavior of the wild-type and three mutant enzymes were measured. The results were consistent with the computational predictions, and the mutant G194P was found to have higher thermostability than the wild-type enzyme.