Table of Contents Author Guidelines Submit a Manuscript
Enzyme Research
Volume 2010, Article ID 125429, 7 pages
http://dx.doi.org/10.4061/2010/125429
Research Article

Cytochrome P450 102A2 Catalyzes Efficient Oxidation of Sodium Dodecyl Sulphate: A Molecular Tool for Remediation

Laboratory of Enzyme Technology, Department of Agricultural Biotechnology, Agricultural University of Athens, Iera Odos 75, 11855 Athens, Greece

Received 5 January 2010; Revised 11 April 2010; Accepted 21 May 2010

Academic Editor: Munishwar Nath Gupta

Copyright © 2010 Irene Axarli et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Bacterial cytochrome P450s (CYPs) constitute an important family of monooxygenase enzymes that carry out essential roles in the metabolism of endogenous compounds and foreign chemicals. In the present work we report the characterization of CYP102A2 from B. subtilis with a focus on its substrate specificity. CYP102A2 is more active in oxidation of sodium dodecyl sulphate (SDS) than any other characterized CYP. The effect of SDS and NADPH concentration on reaction rate showed nonhyperbolic and hyperbolic dependence, respectively. The enzyme was found to exhibit a bell-shaped curve for plots of activity versus pH, over pH values 5.9–8.5. The rate of SDS oxidation reached the maximum value approximately at pH 7.2 and the pH transition observed controlled by two ps in the acidic and basic pH range. The results are discussed in relation to the future biotechnology applications of CYPs.

1. Introduction

Cytochrome P450 monooxygenases (CYPs) play a key role in primary and secondary metabolic pathways and in drug detoxification. CYPs are heme-thiolate proteins which are widely distributed in animals, plants, and microorganisms [1, 2]. They catalyze the oxidation of nonactivated C–H bonds, often in a regio- and stereoselective manner, according to the equation: Among the members of the cytochrome P450 family, the monooxygenase from Bacillus megaterium (CYP102A1 or P450 BM3), it is well characterised. This protein is self-sufficient monooxygenase composed of an N-terminal heme monooxygenase linked to the C-terminal diflavin reductase domain on a single polypeptide chain [1, 36]. The X-ray structure of the monooxygenase domain of CYP102A1 was first solved in 1993 [5]. Since then, several structures of CYP102A1 mutants with and without bound substrates [6] and one including the FMN binding domain [7] have been solved. More recently, a fragment containing the FAD binding site was crystallized and its structure was resolved by X-ray analysis [8].

The genome sequence of the gram-positive model organism Bacillus subtilis 168 encodes eight cytochromes P450 [9]. Of these B. subtilis P450s, four have been wellcharacterized at the enzyme level. For example, BioI (CYP107H1) is involved in the early stages of synthesis of the vitamin biotin [10]. CYP152A1 catalyzes hydrogen peroxide dependent hydroxylation of long-chain fatty acids [11]. The two other P450 monooxygenases CYP102A2 (accession number O08394) and CYP102A3 (accession number O08336) within the Bacillus subtilis genome that shows high similarity to the Bacillus megaterium have been recently identified and characterised [1214].

Because of their broad substrate specificity and catalytic diversity, there is an increasing interest to use P450s in biotechnology [1421]. The products of these enzymes (e.g., enantiomerically pure hydroxylated or epoxidised fatty acids and their derivatives) might find applications in the preparation of fine chemicals for the synthesis of polymers, flavours, fragrances, for the production of pharmaceuticals, or the optimization of lead compounds and existing drugs.

In the present work we report the characterization of CYP102A2 from Bacillus subtilis with a focus on its ability to catalyze the efficient oxidation of sodium dodecyl sulphate.

2. Materials and Methods

2.1. Materials

The pCR T7/CT-TOPO TA Expression Kit was purchased from Invitrogen (USA). NADPH (tetrasodium salt, ca. 95 %), sodium dodecyl sulfate (SDS), lauric acid, all other organic substrates and crystalline bovine serum albumin (BSA) (fraction V) were obtained from Sigma-Aldrich Co. (USA).

2.2. Methods
2.2.1. Cloning, Expression, and Purification of CYP102A2 from E. coli BL21 (DE3) Cells

B. subtilis (ATCC 168) was grown at 30 in a medium containing 1% peptone, 0.5% yeast extract, and 1% NaCl. After 24 h, cells were pelleted by centrifugation and genomic DNA was isolated according to a standard procedure [22]. PCR was used to amplify the full-length gene CYP102A2 from genomic DNA using the oligo primers synthesised to the region of the gene from the ATG start codon (-ATGAAGGAAACAAGCCCGATTCCT CAGCCG-) and to the end of the gene finishing at the CTA stop codon (-TTTAGA TCTCTATATCCCTGCCCAGACATG-). The PCR reaction was carried out in a total volume of 50 L contained 0.1 M of each primer, 5 ng template genomic DNA, 200 μM of each dNTP, 5 Pfu buffer, and 2.1 units of Pfu DNA polymerase (Promega). The PCR procedure comprised 33 cycles of 2 minutes at C, 2 minutes at C, and 6 minutes at C. A final extension time at C for 20 minutes was performed after the 33 cycles. The resulting PCR amplicon was TOPO ligated into a T7 expression vector (pCR T7/CT-TOPO). The resulting expression construct (pCYP102A2) was sequenced along both strands and was used to transform competent BL21(DE3) E. coli cells. E. coli cells, harboring plasmid pCYP102A2, were grown at C in 1 L LB medium containing 100 g/mL ampicillin. The synthesis of CYP102A2 was induced by the addition of 1 mM IPTG when the absorbance at 600 nm was 0.5–0.6. Four hours after induction, cells (approx. 3 g) were harvested by centrifugation at 8000 r.p.m. and C for 20 min. CYP102A2 was purified by a method similar to that described elsewhere [23]. Protein purity was judged by SDS-PAGE. The enzyme was also cloned as a 6His-tagged protein using the same expression vector (pCR T7/CT-TOPO) and employing as and end specific primers the sequences: (-ATGAAGGAAACAAGCCCGATTCCTCAGCCG-) and (-TTTAGATCTCTATATCCCTGCCCAGAC-). The 6His-tagged enzyme was purified using NTA-Sepharose column according to the standard protocol. Protein purity was judged by SDS-PAGE. Kinetic comparison (results not shown) of the tagged and untagged enzyme showed that the extra 6His residues on the C-terminus did not interfere with the activity and function of the enzyme.

2.2.2. Assay of Enzyme Activity, Protein and Steady-State Kinetic Analysis

Enzyme activities were measured by determining the rate of NADPH conversion to and following the decrease of absorbance at 340 nm. The final assay volume of 1 mL contained 0.1 M potassium phosphate buffer, pH 7.2; 0.1 mM NADPH; 1.735 mM SDS and sample containing enzyme activity. One unit of enzyme activity is defined as the amount of enzyme that catalyses the conversion of 1 mol NADPH to NADP+ per minute at C. Protein concentration was determined at C by the method of Bradford using bovine serum albumin (fraction V) as standard [24].

Steady-state kinetic measurements were performed at C in 0.1 M potassium phosphate buffer, pH 7.2 by varying the concentration of the substrates (NADPH, SDS). Initial velocities were determined in the presence of 1.735 mM SDS, while the NADPH concentration range was 6.6–100 M. Michaelis-Menten kinetics was observed under these investigated conditions. The kinetic parameters and were calculated by nonlinear regression analysis of experimental steady-state data. Turnover numbers were calculated on the basis of active site per 119 kDa. Kinetic data were analyzed using the computer program GraFit (Erithacus Software Ltd.) [27]. When NADPH was used at a fixed concentration (0.1 mM), the SDS varied in the range of 0.34–1.73 mM. In this case the data are best fitted to the Hill function since the curves are nonhyperbolic (sigmoidal curves). For each experimental velocity curve, the value, ( is the substrate concentration at which v 0.5), and the Hill coefficient, , were determined by fitting the plotted versus substrate concentration to the Hill equation: Curve-fits were obtained using the GraFit (Erithacus Software, Ltd.) computer program.

2.2.3. Difference Spectroscopy

Difference spectral titrations were performed in a Perkin-Elmer Lamda16 double beam double monochromator UV-VIS spectrophotometer. Enzyme solution (0.8 mL; 0.077 mg enzyme in 20 mM potassium phosphate, pH 7.4) and enzyme solvent (0.8 mL; 20 mM potassium phosphate, pH 7.4) were placed in the sample and reference black-wall silica cuvettes (10 mm pathlength), respectively, and the baseline difference spectrum was recorded in the range 500–360 nm. Identical volumes of sodium dodecyl sulfate (SDS) solution were added to both cuvettes and the difference spectra were recorded after each addition.

2.2.4. pH Dependence of

The pH profile of CYP102A2 was conducted at C in 0.1 M potassium phosphate buffer (pH 5.9–8.5). values were estimated by fitting the experimental data to the equation reported by Blanchard and Cleland 1980 [28], using the computer program GraFit (Erithacus Software Ltd.).

2.2.5. Bioinformatics Analysis and Molecular Modelling

A molecular model of the heme domain of CYP102A2 was constructed using SWISS-MODEL (http://www.expasy.org/swissmod) [29], as described in Axarli et al., 2005 [12]. Atomic contacts in the complex CYP102A2-palmitoleic acid were analyzed by iMolTalk (http://i.moltalk.org/) [30] using a distance threshold set to 4 Angstroms.

3. Results and Discussion

3.1. Effect of pH on Vmax for the SDS Oxidation Reaction

P450 monooxygenases catalyze a broad range of reactions, with different members of the family exhibiting quite varied substrate specificity [12, 23]. The enzyme is more active in oxidation of SDS than any other characterized P450 monooxygenase and therefore detailed kinetic analysis was carried out using SDS and NADPH as substrates.

The effect of pH on the for the SDS oxidation reaction was investigated and the results are shown in Figure 1. CYP102A2 remained catalytically active in a broad pH range. The rate of SDS oxidation reached the maximum value approximately at pH 7.2 (0.1 M potassium phosphate buffer). The pH transition observed seems to be controlled by two in the acidic () and basic pH range: (). The -pH profile yields the , which reflects the ionization of the enzyme in complex with substrate [29]. As there is no ionizable groups on SDS in the pH range 5.9–8.5, one can conclude that these must be for an enzyme side chains.

125429.fig.001
Figure 1: The effect of pH on of the CYP102A2 SDS/NADPH reaction. Steady-state kinetic measurements were performed in 0.1 M potassium phosphate buffers adjusted to different pH values (5.9–8.5).
3.2. Effect of SDS and NADPH Concentration on the Enzyme Activity

The steady-state kinetics for CYP102A2 was determined by monitoring the SDS dependent oxidation of NADPH. The effect of SDS and NADPH concentration on the reaction rate was studied at 37 as shown in Figure 2 and the results are listed in Table 1. The dependence of reaction rate on NADPH and SDS concentration was hyperbolic and nonhyperbolic, respectively. Nonhyperbolic data sets showed sigmoidal characteristics and were accurately fitted by the Hill function. Similar nonhyperbolic kinetics was also observed using lauric acid as substrate (Table 1).

tab1
Table 1: Kinetic parameters of CYP102A2 using SDS and NADPH as substrates. Steady-state kinetic measurements were performed at 37 in 0.1 M potassium phosphate buffer, pH 7.2. All initial velocities were determined in triplicate. The kinetic parameters and for NADPH were calculated by nonlinear regression analysis of experimental steady-state data using the GraFit (Erithacus Software Ltd.) program. The value, , and the Hill coefficient, , for SDS were determined by fitting the plotted versus substrate concentration to the Hill equation using the GraFit (Erithacus Software Ltd.) program.
fig2
Figure 2: Kinetic analysis of CYP102A2 in 0.1 M potassium phosphate buffer, pH 7.2. (a) Initial velocity analysis with NADPH as variable substrate (6.6–100 M) and SDS at saturation concentration. (b) Initial velocity analysis with SDS as variable substrate (0.34–1.73 mM) and NADPH at saturation concentration. The plot of rate versus SDS concentration is nonhyperbolic and is fitted to the Hill function.

The actual role of nonhyperbolic dependence is not understood. However, it appears more likely that several structural changes throughout the extended hydrophobic active sites in the heme domain of CYP102A2 combine to alter the substrate binding modes in this enzyme and to produce the apparent nonhyperbolic dependence of kinetic on the substrate concentration [5]. Alternatively, the nonhyperbolic dependence probably indicates that CYP102A2 is able to bind simultaneously more molecules of SDS substrates. This which probably shows that the binding process reflects cooperative binding of more than one SDS molecule to the active sites of this enzyme. Atypical (non-Michaelis-Menten) kinetic features have also been observed for several drug metabolising enzyme such as CYP3A4 [5]. In the case of CYP3A4 it has been suggested that this atypical kinetics is due to the existence of, and interaction between, several binding sites on the enzyme. The large CYP3A4 active site may allow the simultaneous presence of multiple molecules and the exact binding conformations appear to depend on the substrates involved, their relative concentration, and affinity for the enzyme.

CYP102A2 exhibits higher values compared to CYP102A1 towards identical substrates. [23]. Comparing both the amino acid sequences of CYP102A1 and CYP102A2 (Figure 3) and their structures showed significant differences within the substrate access channel, which might explain the different kinetic constants found for CYP102A2. For example, the fatty acid substrate in the crystal structure of CYP102A1 [3] exhibits total 29 interactions. On the other hand the fatty acid substrate in the modeled structure of CYP102A1 exhibits fewer interactions (total 22 interactions). Arg47 and Tyr51 (Figure 3) in CYP102A1 is of particular importance since it is located at the entrance of the active site and interacts with the carboxylate group of the fatty acid substrate [3]. These residues are considered crucial for the proper positioning of the substrate [31, 32]. In CYP102A2, residue 48 (homologous to Arg47 in CYP102A1) is Gly (Figure 3), which rules out a direct effect on substrate binding. In addition, residue 51 (homologous to Tyr51) is Val in our homology model of CYP102A2. Its side-chain is ~6.5 Å away from the substrate which suggests that Val51 might not interact directly with the substrates. Instead of Arg47, Arg353 might interact directly with the carboxylate group of the fatty acid substrate, although the distance between these groups is rather long (7 Å).

125429.fig.003
Figure 3: Amino acid sequence alignments. Sequence alignments of heme domain of CYP102A1 (residues 1–472) of B. megaterium flavocytochrome P450 BM3 (NCBI accession number A34286) with the respective domain of CYP102A2 from B. subtilis (NCBI accession number O08394). The alignments were produced using Clustal W [25] and visualied using ESPript [26]. The secondary structure of CYP102A1 (pdb 1FAG) and numbering are shown above the alignment. Alpha helices and beta strands are represented as helices and arrows, respectively, and beta turns are marked with TT. Conserved areas are shown shaded. A column is framed, if more than 70% of its residues are similar according to physicochemical properties.

CYP102A2 undergoes inhibition by NADPH at concentrations above 0.1 mM. The apparent Michaelis constant is shown in Table 1. When NADPH was the variable substrate with several fixed concentrations of SDS, an intersecting pattern of Lineweaver-Burk plot was obtained (Figure 4). When the slopes and the intercepts of the primary Lineweaver-Burk plots were replotted against the reciprocals of the second fixed substrate, a linear relationship was obtained. Product inhibition studies were conducted with both NADP+ and NAD+. With NADPH as the variable substrate, NADP+ and NAD+ gave linear competitive inhibition with equal to 0.1 mM and 149.5 mM, respectively. The intersecting initial velocity pattern obtained and the competitive inhibition obtained using NADP+ and NAD+ are all consisted with a rapid-equilibrium ordered bi-bi kinetic mechanism.

125429.fig.004
Figure 4: Kinetic analysis of CYP102A2. Initial velocity analysis of CYP102A2 with NADPH as the variable substrate for several fixed concentrations of SDS (mM): 0.7, (○); 1 mM, (•); 2 mM, (□).
3.3. Spectral Binding Titrations

The tight binding of various long chain saturated and unsaturated fatty acids to P450 BM3 has been demonstrated by several groups [13, 23, 31, 32]. Binding of fatty acids to P450 BM3 induces a shift in equilibrium of the heme iron spin state toward the high-spin form, leading to changes in the absorption spectrum in the Soret region [23]. To analyse the binding of SDS to CYP102A2, difference spectral titrations were performed. Figure 5(a) illustrates the spectral perturbations observed for CYP102A2 on titration with SDS. SDS binding showed perturbation of the heme spectrum of CYP102A2, with shift of the resting (substrate free, low-spin) Soret band from 418 nm toward a new position (at 390 nm), typical of the high-spin forms. Gustafsson et al. [23], have reported that comparison of the extent of spin-state perturbation induced between SDS and other fatty acids shows a significant variation. Far more extensive shifts in spin-state equilibrium are observed for the unsaturated and branched-chain fatty acids than with the SDS.

fig5
Figure 5: SDS binding to CYP102A2. (a) The spectral changes induced on titration of CYP102A2 with SDS are shown. Arrows indicate directions of change of spectra induced by successive additions of SDS. Difference spectral titrations of CYP102A2 (curve a) and the CYP102A2-SDS 0.8653 mM (curve b) and CYP102A2-SDS 1.726 mM (curve c) and CYP102A2-SDS 2.1553 mM (curve d) complexes. (b) The difference absorbance at 420 nm as a function of the total SDS concentration.

Plot of the absorption change versus SDS concentration (Figure 5(b)) was fitted accurately to the Hill function, which probably indicates that the binding process reflects cooperative binding of more than one SDS molecule to the active sites of the enzymes. This finding further supports the results observed by kinetic analysis.

In conclusion, in this report we addressed questions regarding the kinetic properties of CYP102A2, a poorly characterized enzyme among the CYP. Detailed studies of bacterial CYPs are justified because of the considerable biotechnological potential of these enzymes. Further investigations of this interesting P450 monooxygenase are now in progress in this laboratory. This includes the isolation and characterization of mutants, obtained by in vitro directed evolution.

Conflict of Interest

None is declared.

5. Abbreviations Used

BMP, the heme domain of P450 BM3; CPR, NADPH: cytochrome P450 reductase; CYP, cytochrome P450; CYP102A1 or CYP102 or P450 BM3, cytochrome P450 monooxygenase from Bacillus megaterium; CYP102A2, cytochrome P450 monooxygenase from Bacillus subtilis; FAD, flavin adenine dinucleotide; FMN, flavin mononucleotide; NADP+, -nicotinamide-adenine dinucleotide phosphate; NADPH, -nicotinamide-adenine dinucleotide phosphate, reduced form; pCYP102A2, the expression vector pCR T7/CT-TOPO which contains the gene of CYP102A2; SDS, sodium dodecyl sulphate.

References

  1. L. O. Narhi and A. J. Fulco, “Characterization of a catalytically self-sufficient 119,000-dalton cytochrome P-450 monooxygenase induced by barbiturates in Bacillus megaterium,” The Journal of Biological Chemistry, vol. 261, no. 16, pp. 7160–7169, 1986. View at Google Scholar · View at Scopus
  2. D. R. Nelson, L. Koymans, and L. Koymans, “P450 superfamily: update on new sequences, gene mapping, accession numbers and nomenclature,” Pharmacogenetics, vol. 6, no. 1, pp. 1–42, 1996. View at Publisher · View at Google Scholar · View at Scopus
  3. H. Li and T. L. Poulos, “The structure of the cytochrome p450BM-3 haem domain complexed with the fatty acid substrate, palmitoleic acid,” Nature Structural Biology, vol. 4, no. 2, pp. 140–146, 1997. View at Publisher · View at Google Scholar · View at Scopus
  4. L. P. Wen and A. J. Fulco, “Cloning of the gene encoding a catalytically self-sufficient cytochrome P-450 fatty acid monooxygenase induced by barbiturates in Bacillus megaterium and its functional expression and regulation in heterologous (Escherichia coli) and homologous (Bacillus megaterium) hosts,” The Journal of Biological Chemistry, vol. 262, no. 14, pp. 6676–6682, 1987. View at Google Scholar
  5. K. G. Ravichandran, S. S. Boddupalli, C. A. Hasemann, J. A. Peterson, and J. Deisenhofer, “Crystal structure of hemoprotein domain of P450BM-3, a prototype for microsomal P450's,” Science, vol. 261, no. 5122, pp. 731–736, 1993. View at Google Scholar · View at Scopus
  6. D. C. Haines, D. R. Tomchick, M. Machius, and J. A. Peterson, “Pivotal role of water in the mechanism of P450BM-3,” Biochemistry, vol. 40, no. 45, pp. 13456–13465, 2001. View at Publisher · View at Google Scholar · View at Scopus
  7. I. F. Sevrioukova, H. Li, H. Zhang, J. A. Peterson, and T. L. Poulos, “Structure of a cytochrome P450-redox partner electron-transfer complex,” Proceedings of the National Academy of Sciences of the United States of America, vol. 96, no. 5, pp. 1863–1868, 1999. View at Google Scholar
  8. A. J. Warman, O. Roitel, and O. Roitel, “Flavocytochrome P450 BM3: an update on structure and mechanism of a biotechnologically important enzyme,” Biochemical Society Transactions, vol. 33, no. 4, pp. 747–753, 2005. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  9. F. Kunst, N. Ogasawara, and N. Ogasawara, “The complete genome sequence of the gram-positive bacterium Bacillus subtilis,” Nature, vol. 390, no. 6657, pp. 249–256, 1997. View at Google Scholar
  10. M. J. Cryle, N. J. Matovic, and J. J. De Voss, “Products of cytochrome P450Biol (CYP107H1)-catalyzed oxidation of fatty acids,” Organic Letters, vol. 5, no. 18, pp. 3341–3344, 2003. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  11. I. Matsunaga, A. Yamada, and A. Yamada, “Enzymatic reaction of hydrogen peroxide-dependent peroxygenase cytochrome P450s: kinetic deuterium isotope effects and analyses by resonance Raman spectroscopy,” Biochemistry, vol. 41, no. 6, pp. 1886–1892, 2002. View at Publisher · View at Google Scholar · View at Scopus
  12. I. Axarli, A. Prigipaki, and N. E. Labrou, “Engineering the substrate specificity of cytochrome P450 CYP102A2 by directed evolution: production of an efficient enzyme for bioconversion of fine chemicals,” Biomolecular Engineering, vol. 22, no. 1–3, pp. 81–88, 2005. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  13. M. C. Gustafsson, C. N. Palmer, R. C. Wolf, and C. von Wachenfeldt, “Fatty-acid-displaced transcriptional repressor, a conserved regulator of cytochrome P450 102 transcription in Bacillus species,” Archives of Microbiology, vol. 176, no. 6, pp. 459–464, 2001. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  14. O. Lentz, V. Urlacher, and R. D. Schmid, “Substrate specificity of native and mutated cytochrome P450 (CYP102A3) from Bacillus subtilis,” Journal of Biotechnology, vol. 108, no. 1, pp. 41–49, 2004. View at Publisher · View at Google Scholar · View at Scopus
  15. A. Andreadeli, D. Platis, V. Tishkov, V. Popov, and N. E. Labrou, “Structure-guided alteration of coenzyme specificity of formate dehydrogenase by saturation mutagenesis to enable efficient utilization of NADP+,” FEBS Journal, vol. 275, no. 15, pp. 3859–3869, 2008. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  16. M. Budde, M. Morr, R. D. Schmid, and V. B. Urlacher, “Selective hydroxylation of highly branched fatty acids and their derivatives by CYP102A1 from Bacillus megaterium,” ChemBioChem, vol. 7, no. 5, pp. 789–794, 2006. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  17. M. C. Damsten, B. M. A. van Vugt-Lussenburg, T. Zeldenthuis, J. S. B. de Vlieger, J. N. M. Commandeur, and N. P. E. Vermeulen, “Application of drug metabolising mutants of cytochrome P450 BM3 (CYP102A1) as biocatalysts for the generation of reactive metabolites,” Chemico-Biological Interactions, vol. 171, no. 1, pp. 96–107, 2008. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  18. F. P. Guengerich, “Cytochrome P450 enzymes in the generation of commercial products,” Nature Reviews Drug Discovery, vol. 1, no. 5, pp. 359–366, 2002. View at Google Scholar · View at Scopus
  19. C. R. Otey, M. Landwehr, J. B. Endelman, K. Hiraga, J. D. Bloom, and F. H. Arnold, “Structure-guided recombination creates an artificial family of cytochromes P450,” PLoS Biology, vol. 4, no. 5, article e112, 2006. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  20. V. B. Urlacher and S. Eiben, “Cytochrome P450 monooxygenases: perspectives for synthetic application,” Trends in Biotechnology, vol. 24, no. 7, pp. 324–330, 2006. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  21. V. B. Urlacher, S. Lutz-Wahl, and R. D. Schmid, “Microbial P450 enzymes in biotechnology,” Applied Microbiology and Biotechnology, vol. 64, no. 3, pp. 317–325, 2004. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  22. J. Sambrook, E. F. Fritsch, and T. Maniatis, Molecular Cloning: A Laboratory Manual, Vols 1, 2 and 3, Cold Spring Harbor Laboratory Press, New York, NY, USA, 1989.
  23. M. C. U. Gustafsson, O. Roitel, and O. Roitel, “Expression, purification, and characterization of Bacillus subtilis cytochromes P450 CYP102A2 and CYP102A3: flavocytochrome homologues of P450 BM3 from Bacillus megaterium,” Biochemistry, vol. 43, no. 18, pp. 5474–5487, 2004. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  24. M. A. Bradford, “A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein dye binding,” Analytical Biochemistry, vol. 72, no. 1-2, pp. 248–254, 1976. View at Google Scholar · View at Scopus
  25. J. D. Thompson, D. G. Higgins, and T. J. Gibson, “CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice,” Nucleic Acids Research, vol. 22, no. 22, pp. 4673–4680, 1994. View at Google Scholar
  26. P. Gouet, E. Courcelle, D. I. Stuart, and F. Métoz, “ESPript: analysis of multiple sequence alignments in PostScript,” Bioinformatics, vol. 15, no. 4, pp. 305–308, 1999. View at Publisher · View at Google Scholar
  27. R. J. Leatherbarrow, GraFit, Version 3, Erythacus Software, Staines, UK, 1998.
  28. J. S. Blanchard and W. W. Cleland, “Kinetic and chemical mechanisms of yeast formate dehydrogenase,” Biochemistry, vol. 19, no. 15, pp. 3543–3550, 1980. View at Google Scholar
  29. T. Schwede, J. Kopp, N. Guex, and M. C. Peitsch, “SWISS-MODEL: an automated protein homology-modeling server,” Nucleic Acids Research, vol. 31, no. 13, pp. 3381–3385, 2003. View at Publisher · View at Google Scholar
  30. A. V. Diemand and H. Scheib, “iMolTalk: an interactive, internet-based protein structure analysis server,” Nucleic Acids Research, vol. 32, pp. W512–W516, 2004. View at Publisher · View at Google Scholar · View at PubMed
  31. T. W. B. Ost, C. S. Miles, J. Murdoch, Y.-F. Cheung, G. A. Reid, S. K. Chapman, and A. W. Munro, “Rational re-design of the substrate binding site of flavocytochrome P450 BM3,” FEBS Letters, vol. 486, no. 2, pp. 173–177, 2000. View at Publisher · View at Google Scholar
  32. M. G. Joyce, H. M. Girvan, A. W. Munro, and D. Leys, “A single mutation in cytochrome P450 BM3 induces the conformational rearrangement seen upon substrate binding in the wild-type enzyme,” The Journal of Biological Chemistry, vol. 279, no. 22, pp. 23287–23293, 2004. View at Publisher · View at Google Scholar · View at PubMed