Cytochrome P450BM-3 catalyzes NADPH-dependent metabolism of arachidonic acid to nearly enantiomerically pure 18(R)-hydroxyeicosatetraenoic acid and 14(S),15(R)-epoxyeicosatrienoic acid (80 and 20% of total products, respectively). P450BM-3 oxidizes arachidonic acid with a rate of 3.2 ± 0.4 μmol/min/nmol at 30°C, the fastest ever reported for an NADPH-dependent, P450-catalyzed reaction. Fatty acid, oxygen, and NADPH are utilized in an approximately 1:1:1 molar ratio, demonstrating efficient coupling of electron transport to monooxygenation. Eicosapentaenoic and eicosatrienoic acids, two arachidonic acid analogs that differ in the properties of the C-15-C-18 carbons, are also actively metabolized by P450BM-3 (1.4 ± 0.2 and 2.9 ± 0.1 μmol/min/nmol at 30°C, respectively). While the 17,18-olefinic bond of eicosapentaenoic acid is epoxidized with nearly absolute regio- and stereochemical selectivity to 17(S),18(R)-epoxyeicosatetraenoic acid (≥99% of total products, 97% optical purity), P450BM-3 is only moderately regioselective during hydroxylation of the eicosatrienoic acid ω-1, ω-2, and ω-3 sp³ carbons, with 17-, 18-, and 19-hydroxyeicosatrienoic acid formed in a ratio of 2.4:2.2:1, respectively. Based on the above and on a model of arachidonic acid-bound P450BM-3, we propose: 1) the formation by P450BM-3 of a single oxidant species capable of olefinic bond epoxidation and sp³ carbon hydroxylation and 2) that product chemistry and, thus, catalytic outcome are critically dependent on active site spatial coordinates responsible for substrate binding and productive orientation between heme-bound active oxygen and acceptor carbon bond(s).

Cytochrome P450 BM-3 catalyzes the high turnover regio- and stereoselective metabolism of arachidonic and eicosapentaenoic acids. To map structural determinants of productive active site fatty acid binding, we mutated two amino acid residues, arginine 47 and phenylalanine 87, which flank the surface and heme ends of the enzyme's substrate access channel, respectively. Replacement of arginine 47 with glutamic acid resulted in a catalytically inactive mutant. Replacement of arginine 47 with alanine yielded a protein with reduced substrate binding affinity and arachidonate sp³ carbon hydroxylation activity (72% of control wild type). On the other hand, arachidonic and eicosapentaenoic acid epoxidation was significantly enhanced (154 and 137%, of control wild type, respectively). As with wild type, the alanine 47 mutant generated (18R)-hydroxyeicosatetraenoic, (14S,15R)-epoxyeicosatrienoic, and (17S,18R)-epoxyeicosatetraenoic acids nearly enantiomerically pure. Replacement of phenylalanine 87 with valine converted cytochrome P450 BM-3 into a regio- and stereoselective arachidonic acid epoxygenase ((14S,15R)epoxyeicosatrienoic acid, 99% of total products). Conversely, metabolism of eicosapentaenoic acid by the valine 87 mutant yielded a mixture of (14S,15R)- and (17S,18R)-epoxyeicosatetraenoic acids (26 and 69% of total, 94 and 96% optical purity, respectively). Finally, replacement of phenylalanine 87 with tyrosine yielded an inactive protein. We propose that: (a) fatty acid oxidation by P450 BM-3 is incompatible with the presence of residues with negatively charged side chains at the surface opening of the substrate access channel or a polar aromatic side chain in the vicinity of the heme iron; (b) the high turnover regio- and stereoselective metabolism of arachidonic and eicosapentaenoic acids involves charge-dependent anchoring of the fatty acids at the mouth of the access channel by arginine 47, as well as steric gating of the heme-bound oxidant by phenylalanine 87; and (c) substrate binding coordinates, as opposed to oxygen chemistries, are the determining factors responsible for reaction rates, product chemistry, and, thus, catalytic outcome.

There is substantial interest in implementing a bioinformatics tool that allows the design of oligonucleotides to support the development of in vitro gene synthesis. Current protocols to make long synthetic DNA molecules rely on the in vitro assembly of a set of short oligonucleotides, either by ligase chain reaction (LCR) or by assembly PCR. Ideally, such oligonucleotides should represent both strands of the final DNA molecule. They should be adjacent on the same strand and overlap the complementary oligonucleotides from the second strand to ensure good hybridization during assembly. This implies that the thermodynamic properties of each oligonucleotide have to be consistent across the set. Furthermore, any given oligonucleotide has to be totally specific to its target to avoid the creation of incorrectly assembled sequences. We have developed Gene2Oligo (http://berry.engin.umich.edu/gene2oligo/), a web-based tool that divides a long input DNA sequence into a set of adjacent oligonucleotides representing both DNA strands. The length of the oligonucleotides is dynamically optimized to ensure both the specificity and the uniform melting temperatures necessary for in vitro gene synthesis. We have successfully designed and used a set of oligonucleotides to synthesize the Saccharomyces cerevisiae cytochrome b5 by using both LCR and assembly PCR.

The design of a family shuffling strategy (CLERY: Combinatorial Libraries Enhanced by Recombination in Yeast) associating PCR-based and in vivo recombination and expression in yeast is described. This strategy was tested using human cytochrome P450 CYP1A1 and CYP1A2 as templates, which share 74% nucleotide sequence identity. Construction of highly shuffled libraries of mosaic structures and reduction of parental gene contamination were two major goals. Library characterization involved multiprobe hybridization on DNA macro-arrays. The statistical analysis of randomly selected clones revealed a high proportion of chimeric genes (86%) and a homogeneous representation of the parental contribution among the sequences (55.8 ± 2.5% for parental sequence 1A2). A microtiter plate screening system was designed to achieve colorimetric detection of polycyclic hydrocarbon hydroxylation by transformed yeast cells. Full sequences of five randomly picked and five functionally selected clones were analyzed. Results confirmed the shuffling efficiency and allowed calculation of the average length of sequence exchange and mutation rates. The efficient and statistically representative generation of mosaic structures by this type of family shuffling in a yeast expression system constitutes a novel and promising tool for structure–function studies and tuning enzymatic activities of multicomponent eucaryote complexes involving non-soluble enzymes.

Quantitative structure-activity relationships may bring invaluable information on structural elements of both enzymes and substrates that, together, govern substrate specificity. Buried active sites in cytochrome P450 enzymes are connected to the solvent by a network of channels exiting at the distal surface of the protein. This review presents different in silico tools that were developed to uncover such channels in P450 crystal structures. It also lists some of the experimental evidence that actually suggest that these predicted channels might indeed play a critical role in modulating P450 functions. Amino acid residues at the entrance of the channels may participate to a first global ligand recognition of ligands by P450 enzymes before they reach the buried active site. Moreover, different P450 enzymes show different networks of predicted channels. The plasticity of P450 structures is also important to take into account when looking at how channels might play their role.