Advertisement
JBC

HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Originally published In Press as doi:10.1074/jbc.M607946200 on January 23, 2007

J. Biol. Chem., Vol. 282, Issue 14, 10544-10552, April 6, 2007
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental Data
Right arrow All Versions of this Article:
282/14/10544    most recent
M607946200v1
Right arrow Submit a Letter to Editor
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowRequest Permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Mao, W.
Right arrow Articles by Schuler, M. A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Mao, W.
Right arrow Articles by Schuler, M. A.
Social Bookmarking
 Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Reddit   Add to Technorati  
What's this?

Allelic Variation in the Depressaria pastinacella CYP6AB3 Protein Enhances Metabolism of Plant Allelochemicals by Altering a Proximal Surface Residue and Potential Interactions with Cytochrome P450 Reductase*Formula

Wenfu Mao{ddagger}, Sanjeewa G. Rupasinghe§, Arthur R. Zangerl{ddagger}, May R. Berenbaum{ddagger}, and Mary A. Schuler§1

From the Departments of {ddagger}Entomology and §Cell and Developmental Biology, University of Illinois, Urbana, Illinois 61801

Received for publication, August 18, 2006 , and in revised form, January 5, 2007.


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 EXPERIMENTAL PROCEDURES
 RESULTS
 DISCUSSION
 REFERENCES
 
CYP6AB3v1, a cytochrome P450 monooxygenase in Depressaria pastinacella (parsnip webworm), is highly specialized for metabolizing imperatorin, a toxic furanocoumarin in the apiaceous host plants of this insect. Cloning and heterologous expression of CYP6AB3v2, an allelic variant identified in D. pastinacella, reveals that it metabolizes imperatorin at a rate (Vmax of 10.02 pmol/min/pmol of cytochrome P450 monooxygenase (P450)) significantly higher than CYP6AB3v1 (Vmax of 2.41 pmol/min/pmol) when supplemented with even low levels of cytochrome P450 reductase. Comparisons of the NADPH consumption rates for these variants indicate that CYP6AB3v2 utilizes this electron source at a faster rate than does CYP6AB3v1. Molecular modeling of the five amino acid differences between these variants and their potential interactions with P450 reductase suggests that replacement of Val92 on the proximal face of CYP6AB3v1 with Ala92 in CYP6AB3v2 affects interactions with P450 reductase so as to enhance its catalytic activity. Allelic variation at this locus potentially allows D. pastinacella to adapt to both intraspecific and interspecific variation in imperatorin concentrations in its host plants.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 EXPERIMENTAL PROCEDURES
 RESULTS
 DISCUSSION
 REFERENCES
 
Cytochrome P450 monooxygenases (P450s)2 are heme-containing enzymes that catalyze the NADPH-dependent reductive cleavage of molecular oxygen to produce functionalized organic products and a molecule of water, making them crucial for phase I metabolism in a wide range of organisms. Although considerable allelic variation in both the coding and the regulatory regions is known to occur in many vertebrate P450 genes (13), the extent of allelic variation in insect P450 loci is not clear.

Much of what is known of allelic variation in insect P450s results from studies of insecticide resistance. Here, most examples of P450-mediated resistance to insecticides involve regulatory changes, often transposon-mediated (4). Only a few examples exist where point mutations in coding regions affect reactivities toward insecticides. Perhaps the best evidence of such allelic variation is the existence of three amino acid variations in CYP6A2 (R335S, L336V, V476L) in the DDT-resistant RDDTR strain of Drosophila melanogaster; these substitutions occur near the CYP6A2 active site and increase catalytic efficiency of the CYP6A2 enzyme against DDT, 7-ethoxycoumarin, and 7-benzoyloxycoumarin (5). Allelic variation in the form of amino acid substitutions has also been documented in CYP6D1, CYP6D3, and CYP6X1 as well as CYP6A2 (68), but the functional contributions of these to resistance have not yet been determined.

In addition to contributing to insecticide resistance, P450s play an important role in mediating resistance to plant allelochemicals (defense compounds) in herbivorous insects (9). However, these allelochemicals present a toxicological challenge to insects that differs in several fundamental ways from that presented by insecticides. Although insecticides are generally single compounds designed to kill and applied in a broadcast fashion, plant defense chemicals tend to occur as mixtures of structurally similar compounds that are idiosyncratically distributed in the environment and expressed in a site-specific manner (10). Thus, although regulatory changes leading to overexpression may be common in P450-mediated insecticide resistance, allelic variations in coding sequences that affect the relative efficiency of substrate metabolism may be of adaptive value in herbivorous insects, particularly those with a narrow range of chemically similar host plants. Such insects are likely to encounter differences in the relative content and composition of toxins among their hosts, and allelic variation in P450-mediated metabolism may allow specialization at the level of insect populations and evolution of host plant responses to herbivory.

Depressaria pastinacella, the parsnip webworm, is restricted throughout its range on host plants in the genera Pastinaca and Heracleum (11, 12). Like many apioid umbellifers, Pastinaca sativa (wild parsnip) and Heracleum lanatum (cow parsnip) produce phototoxic furanocoumarins capable of intercalating between base pairs in DNA, forming covalent cross-links with the 5,6 double bond of pyrimidines (13). As a consequence of their ability to interfere with DNA replication, furanocoumarins are highly toxic to a wide range of organisms, including insects (14). As shown in Table 1, the amounts of five furanocoumarins present in fruits of both parsnip webworm hosts differ substantially both within species (12) and between host plant species (11). As is the case with several other insects, P450s are responsible for detoxifying furanocoumarins produced in their host plants (1521).


View this table:
[in this window]
[in a new window]

 
TABLE 1
Range of variation in furanocoumarin content of Depressaria pastinacella host plants in North America Data are from Zangerl and Berenbaum (12) with concentration measured in µg/mg seed.

 
The majority of furanocoumarin-metabolizing P450 genes identified to date have been members of the CYP6B subfamily, including CYP6B1 from Papilio polyxenes (15), the CYP6B4 and CYP6B17 groups from Papilio glaucus and Papilio canadensis (19, 22), and CYP6B8, CYP6B9, and CYP6B27 from Helicoverpa zea (23, 24). Furanocoumarin-metabolizing P450s of polyphagous (generalist) species, such as H. zea and P. glaucus, catalyze metabolism of a wide range of furanocoumarins with lower activities than the P450s of specialist species, such as P. polyxenes, which feed exclusively on furanocoumarin-containing host plants (17, 20, 2529). Molecular modelings of these P450s have indicated that at least some of the P450s in generalist insects are structurally more flexible and capable of binding a more diverse array of compounds when compared with some of the P450s in specialist insects (26, 27).

To date, the P450 associated with furanocoumarin metabolism in D. pastinacella CYP6AB3v1 has the narrowest range of substrates among xenobiotic-metabolizing P450s yet identified (21). Molecular analysis has indicated that this P450 transcript is strongly induced by the mixture of furanocoumarins extracted from the reproductive parts of P. sativa, allowing it to be cloned along with a few other P450 sequences (30). A molecular model of CYP6AB3v1 contains three-dimensional elements similar to those predicted in the catalytic sites of the furanocoumarin-metabolizing CYP6B proteins in P. polyxenes, P. glaucus, and H. zea (2527, 31) including residues in the B-helix and B'-C loop of SRS1, the I-helix of SRS4, and the beta-turn in beta-sheet 4 of SRS6 (30). Amino acids are also conserved in SRS1, SRS4, and SRS5 that contribute to form the catalytic site of CYP6B1 from the specialist P. polyxenes, CYP6B4 from the generalist P. glaucus, and CYP6B8 from the generalist H. zea. Specifically, Arg101, Val103, Glu111, Leu113, Asn116, Phe118, and Asp121 in SRS1 of D. pastinacella CYP6AB3 are absolutely conserved in all four of these CYP6B proteins (30); because of small deletions in these other proteins, these CYP6AB3 positions correspond to positions 102–119 in CYP6B1 and CYP6B4 and positions 104–121 in CYP6B8. As well, 9 of 14 positions in SRS4 and three of nine positions in SRS5 are conserved in these four proteins. Among these conserved positions, Phe116, Ala303, and Thr307 in CYP6B1 (corresponding to Phe118, Ala314, and Thr318 in CYP6AB3) have been identified as contacts for the furanocoumarin substrates of CYP6B1, and Val104 is predicted in the catalytic site just out of substrate contact range (25, 31, 32).

The CYP6AB3v1 protein heterologously expressed in Sf9 cells is highly specialized for metabolizing imperatorin with Vmax and Km of 2.41 pmol/min/pmol of P450 and 94.28 µM, respectively (21). As shown in Table 1, this furanocoumarin is at the same time the most abundant furanocoumarin in both North American hosts of D. pastinacella and the most variable in concentration. To determine whether CYP6AB3 allelic variants exist that metabolize imperatorin more effectively, a cDNA library constructed with RNA from the midguts of parsnip webworm larvae induced with parsnip seed furanocoumarins was probed with the original CYP6AB3v1 cDNA. Co-expression of the closely related CYP6AB3v2 cDNA isolated in this process with house fly NADPH-dependent P450 reductase has indicated that this second P450 variant metabolizes imperatorin significantly more rapidly than CYP6AB3v1. Biochemical analyses coupled with molecular modeling provide insight into the mechanistic reasons for the enhanced activity of CYP6AB3v2, which differs in only five positions from CYP6AB3v1. The fact that mutagenesis of Val92 to Ala92 in CYP6AB3v1 caused it to metabolize imperatorin nearly as effectively as CYP6AB3v2 indicates that this amino acid switch is the one primarily responsible for the enhanced activity of CYP6AB3v2 toward imperatorin.


    EXPERIMENTAL PROCEDURES
 TOP
 ABSTRACT
 INTRODUCTION
 EXPERIMENTAL PROCEDURES
 RESULTS
 DISCUSSION
 REFERENCES
 
Reagents—XbaI and BamHI restriction enzymes were bought from Invitrogen. Angelicin, bergapten, isopimpinellin, sphondin, and xanthotoxin were purchased from Indofine Chemical Co. (Belle Mead, NJ). Imperatorin was bought from Feinbiochemica (Heidelberg, Germany). NADPH was obtained from Sigma. Sf9 insect cells, SF-900 serum-free medium, pFast-Bac1 vector, and DH10BAC competent cells were purchased from Invitrogen. Penicillin/streptomycin was obtained from Bio-Whittaker (Walkersville, MD). HPLC solvents were bought from Fisher.

cDNA Library Selection—The full-length cDNA library used by Li et al. (30) was rescreened with the CYP6AB3v1 cDNA at medium stringency (65 °C, 5x SSC, 1% SDS, 5x Denhardt's solution). From 2000 clones, five positive clones over 1.6 kb in length were sequenced using T3 and T7 vector primers and CYP6AB3-specific internal primers. All sequences were defined by sequencing both cDNA strands.

Expression of CYP6AB3v2 and CYP6AB3v1 V92A Mutant in Sf9 Cells—The CYP6AB3v2 cDNA was inserted into the pFast-Bac vector using a PCR strategy with two primers synthesized against the N- and C-terminal ends of the coding sequence. The CYP6AB3v2 forward primer, C3-R (5'-GGCGGATCCGATGTATTTTTTAATTGCATTAG-3') contains the translation start codon and a BamHI site on its 5' end, and the reverse primer, C3-L (5'-CGCTCTAGATTAAACCTTTCTTGCAACCAATG-3'), contains the stop codon and an XbaI site on its 5' end. After PCR amplification of the CYP6AB3v2 clone, the PCR products were digested with BamHI and XbaI and ligated into BamHI-XbaI-cut pFastBac vector.

The CYP6AB3v1 V92A mutant was constructed using the CYP6AB3v1-containing pFastBac clone described in Mao et al. (21) as a template for single-stranded DNA site-directed mutagenesis (33). For this, the 5'-GAACTAATAAAGAGGGTCATGGCTACGGATTTTCAATACTTCCAC-3' primer (mutant nucleotides are underlined) was used instead of the pair of primers in the QuikChangeTM site-directed mutagenesis kit developed by Stratagene (La Jolla, CA). The PCR amplification was performed using the method of Humma et al. (33) except that annealing was done at 60 °C. The final mutants were confirmed by sequencing of the full length of the coding region with vector and internal primers.

All procedures for expression of the recombinant viruses in Sf9 cells were performed as described by Invitrogen protocols. As described previously for CYP6AB3v1 (21), one plate (100 x 20 mm) was seeded with 8 x 106 cells in 10 ml of SF-900 serum-free medium supplemented with 8–10% fetal bovine serum, 50 µg/ml streptomycin sulfate, and 50 units/ml penicillin and co-infected 1 h later with the recombinant CYP6AB3 viruses and house fly P450 reductase virus at varying multiplicities of infection (m.o.i.). Hemin was added to a final concentration of 2 µg/ml 24 h after the initial infection. Baculovirus-infected insect cells were harvested 72 h after infection by centrifugation at 1000 x g for 5 min and washed two times with 0.1 M ice-cold phosphate buffer (pH 7.8) (PB). The cell pellet was resuspended in one-fifth volume of ice-cold cell lysate buffer (0.1 M phosphate buffer (pH 7.8), 1.1 mM EDTA, 0.1 mM dithiothreitol, 0.5 mM phenylmethylsulfonyl fluoride, 5 µg/ml (w/v) leupeptin, 20% glycerol) and lysed by sonication on ice (3 x 15 s). After the sonicated extract was centrifuged at 3200 x g for 5 min at 4 °C, the supernatant was recentrifuged at 112,500 x g for 1 h at 4 °C to collect microsomes, and these were resuspended in one-tenth volume of cold lysate buffer. The concentration of expressed P450s was measured by carbon monoxide (CO) difference analysis (34), and aliquots of the expressed P450s were stored at -80 °C.

Analysis of Wild-type and Mutant Metabolic Activities—Reaction mixtures of 500 µl were set up with 4 µl of a 5 mM furanocoumarin stock in methanol, 50 µl of NADPH (5 mg/ml in 0.1 M PB, pH 7.8) for different treatments or 50 µl of PB for the control without NADPH, 30 pmol of CYP6AB3v2 expressed at an m.o.i. ratio of 1:0.1 (CYP6AB3v2/house fly P450 reductase), and a suitable volume of 0.1 M PB (pH 7.8). After incubating for 30 min in a shaking water bath at 30 °C, all the reactions were terminated by incubating for 5 min at 70 °C to avoid breakdown of imperatorin derivatives. After this inactivation, 5 µl of 1 mM psoralen was added as an internal control, and unmetabolized furanocoumarins were extracted with one 500-µl aliquot of ethyl acetate and analyzed by normal phase HPLC with a solvent system containing 55% cyclohexane, 42% isopropyl ether, 3% amyl alcohol.

CYP6AB3v2 and CYP6AB3v1 V92A mutant activities at varying substrate concentrations were determined using final concentrations of 10, 20, 30, 40, 60, and 80 µM imperatorin in reactions that were terminated at 15 min. Assays were replicated four times, and the data were analyzed using Prism 3.03 software.

Variations in m.o.i. Ratios—Five m.o.i. ratios of the recombinant CYP6AB3 virus to recombinant house fly P450 reductase virus (1:0.01, 1:0.05, 1:0.1, 1:0.5, 1:1) were used to infect Sf9 cells in six plates of 8 x 106 cells/plate to determine the effect of varying ratios on CYP6AB3v2 protein levels. Four ratios (1:0, 1:0.1, 1:0.25, 1:0.5) were used for experiments on CYP6AB3v2 and CYP6AB3v1 V92A mutant activities with each treatment replicated three times. After preparing microsomal samples, their P450 concentrations and activities were determined as above using imperatorin as the substrate with the rates of metabolism expressed as pmol of imperatorin metabolized/min/pmol of P450, and means were compared by analysis of variance.

Comparison of NADPH Consumption Rates—To compare the NADPH consumption rates of CYP6AB3v2, CYP6AB3v1, and its V92A variant co-expressed with house fly P450 reductase at m.o.i. ratios of 1:0.25, 5:0.25, and 1:0.25, respectively, 500-µl reaction mixtures containing 1 µl of 2.5 mM imperatorin in methanol, 20 pmol of a P450, and a suitable volume of 0.1 M PB (pH 7.8) were mixed and loaded into a 500-µl microcuvette. For all, the sample baseline was monitored at 340 nm on a Cary 100 spectrophotometer at room temperature until it stabilized and zeroed. After monitoring for 1 min at 340 nm, the reaction was initiated by the addition of 1 µl of 2.5 mg/ml NADPH (6.7 µM final concentration), and NADPH consumption was monitored for another 9 min.

Purification and Characterization of Imperatorin Metabolites—Imperatorin metabolism reactions by wild-type CYP6AB3v1 and CYP6AB3v2 proteins were carried out in 500-µl reactions as for substrate determinations except that CYP6AB3v1 was co-expressed at an m.o.i. ratio of 1:0.5 and CYP6AB3v2 was expressed at an m.o.i. ratio of 1:0.1. After 500 µl of the ethyl acetate-extracted products and remaining substrate was dried in a SpeedVac centrifuge, the residue was dissolved with 500 µl of methanol, and 50 µl of this solution was loaded on reverse phase 4.7 x 150-mm C-18 Novapak column (Waters Corp., Franklin, MA). The sample was gradient-eluted at a flow rate of 1 ml/min with solvent A composed of 5% MeOH, 1.9% glacial acetic acid, 93.1% water by volume and solvent B being 100% acetonitrile. The CYP6AB3v2 imperatorin metabolite peak was collected and MS-MS analyzed at the Mass Spectrometry Service Facility at the University of Illinois using solvent A alone for 5 min followed by a linear gradient changing from 0 to 100% solvent B over 50 min and then 100% solvent B for 5 min. The absorbance at 254 nm was monitored using a Waters 996 photodiode array detector.

Homology Modeling—The molecular model of CYP6AB3v2 was constructed using MOE programs (Chemical Computing Group Inc., Montreal, Canada) by introducing four amino acid replacements (V92A, E164Q, K214R, I303L listed with the residue in CYP6AB3v1 first) into the CYP6AB3v1 model described in Mao et al. (21) and reminimizing the model using the CHARMM22 force field.

The homology model of the FMN-binding domain of house fly P450 reductase was generated in MOE by aligning the house fly P450 reductase sequence with the FMN-binding domain of Bacillus megaterium P450 reductase (Protein Data Bank code 1BVY [PDB] ) (35) and rat P450 reductase (Protein Data Bank code 1AMO [PDB] ) (36). Ten models were generated using 1AMO as the primary template, the best model (as ranked by MOE) was selected, and the FMN coordinates from the rat P450 reductase were copied into it. The model was energy-minimized using the CHARMM22 force field, and the complex between CYP6AB3v2 and the FMN domain of house fly P450 reductase was constructed by overlaying the CYP6AB3v2 and the P450 reductase models with the crystal structure of the CYP102-P450 reductase FMN domain (1BVY) (35).


Figure 1
View larger version (71K):
[in this window]
[in a new window]

 
FIGURE 1.
Amino acid alignments. The derived sequences of the CYP6AB3v1 (AY295773), CYP6AB3v2 (EF111024), and CYP6AB7 (EF111023) proteins are shown with substrate recognition sequences (SRS1–6) underlined in red. The amino acids different between these three proteins are highlighted in magenta, green, and yellow.

 

    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 EXPERIMENTAL PROCEDURES
 RESULTS
 DISCUSSION
 REFERENCES
 
cDNA Cloning and Characterization—Of the 2000 clones screened from the furanocoumarin-induced parsnip webworm midgut cDNA library (31) with the existing CYP6AB3v1 cDNA, 14 positive clones were obtained including four CYP6AB3v1 (C4, 8, 11, 26), five CYP6AB3v2 (C3, 6, 7, 10, 23), and four CYP6AB7 clones (C1, 2, 12, 18). Alignments of these nucleotide sequences (Supplemental Fig. 1) indicate that individual cDNAs have the same lengths in their coding region and significantly different lengths to their polyadenylation sites. Alignments of the CYP6AB3v1, CYP6AB3v2, and CYP6AB7 protein sequences (Fig. 1) indicate that CYP6AB3v2 shares 99 and 87.9% amino acid identity with CYP6AB3v1 and CYP6AB7, respectively. Differences between the CYP6AB3v1 and CYP6AB3v2 variants are limited to five amino acids (H2Y, V92A, E164Q, K214R, I303L listed with the residue in CYP6AB3v1 first) with only one of these (K214R) occurring in SRS2.

Effect of Co-expressed House Fly P450 Reductase on Expression and Bioactivity of CYP6AB3v2— For comparison of wild-type CYP6AB3v2 activities with and without co-expressed insect P450 reductase, CYP6AB3v2 recombinant virus was expressed in Sf9 cells with varying m.o.i. between 1/0.01 and 1/1 (CYP6AB3v2/P450 reductase) (Fig. 2). As in the expression of the first CYP6AB3v1 variant (21) and other P450s in Sf9 cells, the expression level of CYP6AB3v2 as monitored by carbon monoxide (CO) difference analysis decreased as the m.o.i. of recombinant house fly P450 reductase increased; at an m.o.i. ratio of 1:1, the P450 concentration of CYP6AB3v2 was just 0.1 µM (Fig. 2a). However, in contrast to the many other P450s that we have analyzed, this variant metabolizes imperatorin at a rate of 8.0 pmol/min/pmol of P450 when supplemented with even the lowest levels of P450 reductase (m.o.i. ratio 1:0.1), and this activity remains constant as m.o.i. ratios increase to 1:0.5 (Fig. 2b). In the absence of co-expressed P450 reductase, CYP6AB3v2 metabolizes imperatorin at a rate of 2.39 pmol/min/pmol of P450.

Substrate Reactivities, Kinetics, and NADPH Consumption Rates of the Wild-type CYP6AB3v2 and CYP6AB3v1 Proteins— Analysis of other furanocoumarins (xanthotoxin, bergapten, isopimpinellin, angelicin, sphondin) indicated that none were metabolized by CYP6AB3v2 even when P450 reductase was co-expressed at an m.o.i. ratio of 1:0.1. Additional experiments performed to define the kinetics of CYP6AB3v2 enzyme indicated that it metabolizes imperatorin with Vmax and Km values of 14.18 pmol/min/pmol of P450 and 63.22 µM, respectively, at an m.o.i. ratio of 1:0.25 (Fig. 3). In comparison, wild-type CYP6AB3v1 metabolizes imperatorin with Vmax and Km values of 2.41 pmol/min/pmol of P450 and 94.28 µM, respectively, at non-limiting P450 reductase levels (21).


Figure 2
View larger version (40K):
[in this window]
[in a new window]

 
FIGURE 2.
Effect of co-expressed house fly P450 reductase on synthesis of CYP6AB3 protein and associated activities. Sf9 cells were co-infected with recombinant CYP6AB3v2 virus at a constant m.o.i. of 1, and recombinant house fly P450 reductase viruses at varying m.o.i. were cultured for 72 h after infection and monitored for P450 content by CO difference analyses (a) and imperatorin metabolism by CYP6AB3v2 and CYP6AB3v1 V92A mutant proteins (b). Results are the mean ± S.E. of two independent replicates each done in triplicate.

 


Figure 3
View larger version (24K):
[in this window]
[in a new window]

 
FIGURE 3.
CYP6AB3v2 and CYP6AB3v1 V92A mutant kinetics. Wild-type CYP6AB3v2 and CYP6AB3v1 V92A mutant proteins co-expressed with house fly P450 reductase at an m.o.i. ratio of 1:0.25 were assayed for their activities at varying imperatorin concentrations with 30 pmol of P450/500 µl of reaction. These results represent the mean of two independent replicates, each done in quadruplicate.

 
As the above experiments indicate, CYP6AB3v2 reaches its highest activity when co-expressed with very low levels of P450 reductase (m.o.i. ratio of 1:0.1), suggesting that CYP6AB3v2 couples more efficiently with P450 reductase than CYP6AB3v1 with P450 reductase. Imperatorin-dependent NADPH consumption rates were monitored in microsomes derived from Sf9 cells co-expressing CYP6AB3v1, CYP6AB3v2, or CYP6AB3v1 V92A mutant and P450 reductase at m.o.i. ratios of 5:0.25, 1:0.25, and 1:0.25, respectively. In these, the rate of NADPH consumption was significantly faster for CYP6AB3v2 and CYP6AB3v1 V92A reactions containing equal P450 concentrations (as monitored by CO difference analyses) than for CYP6AB3v1 (Fig. 4).


Figure 4
View larger version (37K):
[in this window]
[in a new window]

 
FIGURE 4.
Comparison of NADPH consumption rates. NADPH consumption curves obtained with equal concentrations of P450 (as defined by CO difference analysis) are graphed for wild-type CYP6AB3v1 and wild-type CYP6AB3v2 and CYP6AB3v1 V92A mutant proteins all co-expressed with house fly P450 reductase.

 
Imperatorin Metabolites Generated by CYP6AB3v1 and CYP6AB3v2—Considering that the K214R replacement in CYP6AB3v2 occurs in one of the six SRS regions that potentially affect substrate orientation and reactivities, the metabolites of imperatorin generated by CYP6AB3v2 were determined and compared with those generated by CYP6AB3v1. HPLC analyses of these metabolites (Fig. 5) provided the initial evidence that those produced by CYP6AB3v1 and CYP6AB3v2 were identical. Comparisons of the spectra for these metabolites with that of imperatorin indicated that the absorption peak at 222.6 nm was decreased and that the absorption peaks at 249.7 and 301.6 nm remained constant as typical for furanocoumarin spectra (Fig. 5). Subsequent tandem MS-MS analysis of the imperatorin metabolite at m/z 286.9 indicated that it fragmented to three peaks at m/z 269.2, 245.1, and 203.2 (Fig. 6a). Among these, the peak at m/z 203 contains the furanocoumarin core structure, the peak at m/z 269.2 corresponds to the loss of one oxygen atom from the core structure, and the peak at m/z 245.1 corresponds to the loss of a CH2(CH3)2 group (Fig. 6b). These data indicate that the metabolite arises as a result of cleavage of the double bond on the isoprenyl side chain of imperatorin and the formation of an epoxide intermediate as shown previously for CYP6AB3v1 (21). We have concluded that the metabolites generated by CYP6AB3v1 and CYP6AB3v2 are the same.

Homology Modeling—To compare these two P450 variants in the most effective manner, the three-dimensional model of CYP6AB3v2 was generated by computationally mutating the CYP6AB3v1 model (21) and minimizing the CYP6AB3v2 replacement model. Overlays of these models indicate that the V92A residue varying between CYP6AB3v1 and CYP6AB3v2 occurs in part of the B-helix on the proximal surface of these proteins. Relevant to the differences observed in NADPH consumption rates, this V92A variation occurs in a region predicted to interact with the FMN domain of the P450 reductase and lies close to the heme (Fig. 7). The model of the CYP6AB3v2-FMN domain complex that we have built indicates that several hydrogen bonds are possible between Ala92 in CYP6AB3v2 and P450 reductase in this region (not shown). In this position, substitution of Ala92 in CYP6AB3v2 for the larger Val92 in CYP6AB3v1 is likely to disrupt this hydrogen-bonding network and other interactions important for electron transfer, thereby reducing substrate turnover and NADPH consumption rates in CYP6B1v1. Of the three other variations, the E164Q variation is located at the C terminus of the D-helix on the outer surface of the protein with no possible contacts with P450 reductase, the K214R variation is located at the C-terminal of the F-helix on the outer surface with no possible contacts with P450 reductase, and I303L is located at the C terminus of the I-helix near predicted exit pathways 1 and 2C (designated as in Wade et al. (37)) with no possible contacts with P450 reductase. Of these last two residues, there is some possibility that residue 214 contacts the substrate and residue 303 affects product exit.


Figure 5
View larger version (25K):
[in this window]
[in a new window]

 
FIGURE 5.
Reaction products derived from wild-type CYP6AB3v1 and CYP6AB3v2 proteins. The products derived from imperatorin in reactions containing wild-type CYP6AB3v1 and CYP6AB3v2 co-expressed with house fly P450 reductase at m.o.i. ratios of 1:0.5 and 1:0.1, respectively, were fractionated on a reverse phase C18 HPLC column and compared with control reactions lacking NADPH. The spectra of imperatorin and its metabolites are shown in the insets above and below one another.

 
Substrate Reactivity of the CYP6AB3v1 V92A Mutant—To conclusively determine whether the V92A replacement in the CYP6AB3v1 variant is responsible for the higher activity of the CYP6AB3v2 variant, this single site change was introduced into the wild-type CYP6AB3v1 protein, and the rate of imperatorin metabolism by the CYP6AB3v1 V92A mutant was monitored at varying house fly P450 reductase virus levels. Like the CYP6AB3v2 protein expressed without P450 reductase supplementation, the CYP6AB3v1 V92A mutant protein also metabolizes imperatorin but at a lower velocity (0.5 pmol/min/pmol of P450) than for the CYP6AB3v2 protein (Fig. 2b). Supplementation with increasing amounts of P450 reductase virus enhances activity of the CYP6AB3v1 V92A mutant protein but not as significantly at the lowest P450 reductase viral titer (1:0.1) as seen for the CYP6AB3v2 protein. Its activity remains constant (6.0–6.4 pmol/min/pmol of P450) at m.o.i. ratios higher than 1:0.25 and slightly lower than the activity of CYP6AB3v2 (7.8–8.2 pmol/min/pmol of P450). Further analyses shown in Fig. 3 have determined that the CYP6AB3v1 V92A mutant protein metabolizes imperatorin with Vmax and Km values of 11.38 pmol/min/pmol of P450 and 32.10 µM, respectively. NADPH consumption assays have indicated that the CYP6AB3v1 V92A mutant protein consumes NADPH at a rate substantially more similar to the CYP6AB3v2 protein than to wild-type CYP6AB3v1 protein (Fig. 4).


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 EXPERIMENTAL PROCEDURES
 RESULTS
 DISCUSSION
 REFERENCES
 
Advances in x-ray crystallography, homology modeling, and site-directed mutagenesis have yielded crucial insights into the structural basis of P450 catalytic efficiencies. With all previous studies identifying residues in insect P450 catalytic sites affecting substrate recognition, turnover, and product release (25, 31, 38, 39), this study is the first to identify a residue affecting interactions with P450 reductase, the one electron transfer partner essential to endoplasmic reticulum-localized P450s. The natural evolution of residue 92 in CYP6AB3v1 from valine to the smaller alanine side chain enhances the NADPH consumption rate of the CYP6AB3v2 protein 2.4-fold when the initial slopes on these curves are compared.


Figure 6
View larger version (17K):
[in this window]
[in a new window]

 
FIGURE 6.
Determination of imperatorin metabolite structure using MS-MS methods. The MS-MS spectrum of the 286 dalton purified metabolite of imperatorin by wild-type CYP6AB3v2 and its corresponding fragmentation pattern are shown in a and b, respectively.

 
As early as 1988, Nadler and Strobel (40) found that chemical modification of P450 reductase affected electrostatic interactions between vertebrate P450s and P450 reductase that were necessary for docking and electron transfer step(s) between these proteins. Structural comparisons of the vertebrate P450s with crystal structures of bacterial CYP101 and CYP102 available in the earliest of these investigations suggested that residues in the B-, C-, G-, and L-helices might interact with residues in P450 reductase. Subsequent analyses of site-directed replacements on several vertebrate endoplasmic reticulum-localized P450s effectively demonstrated that residues in the B-, J'-, K-, and L-helices of human CYP17A1 (41, 42), in the C-helix, F-G loop and the heme-binding region of rabbit CYP2B4 (43, 44) and in the B- and G-helices, G-H loop, and beta2 (2) and beta1 (3) strands of CYP1A1 and CYP1A2 (45, 46) were involved in the formation of complexes with P450 reductase and/or cytochrome b5 (47). The V92A replacement in our newly defined CYP6AB3v2 variant exists in a region of the B-helix that is predicted to be positioned on the external surface of the CYP6AB3 protein in close proximity to residues in the beta2-{alpha}3 loop and the beta1-{alpha}2 loop and N terminus of {alpha}4.


Figure 7
View larger version (61K):
[in this window]
[in a new window]

 
FIGURE 7.
Model of the CYP6AB3-P450 reductase complex. The model of D. pastinacella wild-type CYP6AB3v2 is shown with the model of the FMN-binding domain of house fly P450 reductase with interactions in the B-helix region of the P450 shown in stick format and the FMN domain of the P450 reductase shown in ball-and-stick format. The position of the V92A variation between wild-type CYP6AB3v1 and CYP6AB3v2 is shown with Val92 in orange and Ala92 in aqua.

 
The electrostatic potential of the loop between strand 5 and helix F on the P450 reductase surface is largely negative (36) with clusters of negatively charged residues proposed to be involved in the recognition of specific positively charged residues on their associated P450 (48). Mutagenesis of specific residues within two negatively charged clusters in the rat P450 reductase (cluster I: Asp207-Asp208-Asp209; cluster II: Glu213-Glu214-Asp215) has now demonstrated that negatively charged residues within the FMN domain are important for the interaction of rat P450 reductase with P450s (49). Positively charged residues on the proximal P450 surface play an important role in this interaction cluster near the center of the proximal face and heme-binding domain and consist of a four-helix bundle (the D-, E-, I-, and L-helices), the J- and K-helices and the "meander" region (5052).

Site-directed mutagenesis and chemical modifications have identified basic residues important for the P450 reductase interactions on the proximal surface of the P450. These residues are located in three distinct secondary structural elements and include the B-helix, the C-helix, and the loop between the K''-helix and heme-binding region. Specifically, Arg122, Arg126, Arg133, and Lys139 in the C-helix and Lys433, Arg422, and Arg443 in the loop between the K''-helix and heme-binding region located on the proximal surface of CYP2B4 participate in the interaction with the P450 reductase (43). Lys99 and Lys105 in the B-helix, Lys440, Lys453, Lys463, and Arg455 in the loop between the K''-helix and heme-binding region, as well as an Arg cluster (Arg135-Arg136-Arg137) in the C-helix on the proximal surface of CYP1A2 might also be involved in the interaction with P450 reductase (45). Chemical modification of Lys384 (N-terminal of the beta1–3 strand), Arg422, Lys433, and Arg473 (loop between the K''-helix and heme-binding region) in CYP2B1 (corresponding to the same residues in CYP2B4) caused 95% loss of benzphetamine demethylation activity (46, 53). Basic arginine and lysine residues on the surface of CYP6AB3 including Lys126, Arg129, Arg131, Lys140, Lys142, Lys441, Arg452, and Arg458 correspond to some of these clusters in CYP2B4 and CYP1A2. Of these, Lys126, Arg129, Arg131, Lys140, and Lys142 locate to the C-helix, and Lys441 and Arg458 locate to the loop between the K''-helix and the heme-binding region. In the three-dimensional P450 structure, the B-helix and the K''-helix to the heme-binding region loop come close to each other on the proximal surface, bringing Val92 (B-helix) close to the basic residue cluster on loop between the K''-helix and heme-binding region and likely account for the dramatic effects that this particular residue has in defining the activity of CYP6AB3. The 2.4-fold increase in NADPH consumption rate by the CYP6AB3v2 protein is in the same range as the 4.3-fold increase in Vmax for this enzyme. Other factors contributing to the higher activity of this variant may be slight differences in substrate binding affinity but, because imperatorin does not induce strong spectral shifts upon binding to these proteins, substrate binding affinities cannot be reliably measured.

Allelic variation at the CYP6AB3 locus may be instrumental in allowing D. pastinacella to utilize host plants that differ quantitatively and qualitatively in furanocoumarin composition; although all of its host plants in North America contain high concentrations of furanocoumarins, the precise content and composition vary between species and among populations within species. The parsnip webworm displays a remarkable capacity for chemical phenotype matching, according to which the distribution and abundance of furanocoumarins in host plants is mirrored by the detoxification profile of the local webworm population (12, 54, 55). Given the extreme substrate specificity of CYP6AB3 relative to furanocoumarins, there are likely other P450s in this species that contribute to furanocoumarin metabolism; allelic variation at these other loci may allow for the rapid local adaptation that takes place in this herbivore.


    FOOTNOTES
 
* This project was funded by National Science Foundation Grants DEB 0235773 and 0612376 (to M. R. B. and A. R. Z.) and National Institutes of Health Grant R01 GM071826 (to M. A. S.). The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. Back

Formula The on-line version of this article (available at http://www.jbc.org) contains supplemental DNA sequence comparisons. Back

1 To whom correspondence should be addressed: Dept. of Cell and Developmental Biology, University of Illinois at Urbana-Champaign, 1201 W. Gregory Dr., 161 Edward R. Madigan Laboratory (ERML), Urbana, IL 61801. Tel.: 217-333-8784; Fax: 217-244-1336; E-mail: maryschu{at}uiuc.edu.

2 The abbreviations used are: P450, cytochrome P450 monooxygenase; CO, carbon monoxide; DDT, 1,1,1-trichloro 2,2-bis(p-chlorophenyl)ethane; m.o.i., multiplicity of infection; HPLC, high pressure liquid chromatography; MS-MS, tandem mass spectrometry; PB, phosphate buffer. Back



    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 EXPERIMENTAL PROCEDURES
 RESULTS
 DISCUSSION
 REFERENCES
 

  1. Ingleman-Sundberg, M. (2002) Toxicology 181–182, 447-452[CrossRef][Medline] [Order article via Infotrieve]
  2. Ingleman-Sundberg, M. (2005) Pharmacogenomics J. 5, 6-13[CrossRef][Medline] [Order article via Infotrieve]
  3. Jinno, H., Tanaka-Kagawa, T., Ohno, A., Makino, Y., Matsushima, E., Hanioka, N., and Ando, M. (2003) Drug Metab. Dispos. 31, 398-403[Abstract/Free Full Text]
  4. Li, X., Schuler, M. A., and Berenbaum, M. R. (2006) Annu. Rev. Entomol. 52, 231-253[CrossRef]
  5. Amichot, M., Tarès, S., Brun-Barale, A., Arthaud, L., Bride, J.-M., and Bergé, J.-B. (2004) Eur. J. Biochem. 271, 1250-1257[Medline] [Order article via Infotrieve]
  6. Feyereisen, R. (2005) in Comprehensive Molecular Science (Gilbert, L. I., Iatrou, K., and Gill, S. S., ed), Vol. 4, pp. 1-77, Elvesier Press, Oxford, UK[CrossRef]
  7. Kasai, S., and Scott, J. G. (2001) Insect Mol. Biol. 10, 191-196[CrossRef][Medline] [Order article via Infotrieve]
  8. Zhu, Y. C., and Snodgrass, G. L. (2003) Insect Mol. Biol. 12, 39-49[CrossRef][Medline] [Order article via Infotrieve]
  9. Berenbaum, M. R. (1999) in Molecular Biology of the Toxic Response, pp. 553-571, Taylor & Francis, Philadelphia, PA
  10. Berenbaum, M. R. (1995) Proc. Natl. Acad. Sci. U. S. A. 92, 2-8[Abstract/Free Full Text]
  11. Berenbaum, M. R., and Zangerl, A. R. (1991) Oikos 62, 153-159[CrossRef]
  12. Zangerl, A. R., and Berenbaum, M. R. (2003) Evolution 57, 806-815[CrossRef][Medline] [Order article via Infotrieve]
  13. Murray, R. D. H., Mendez, J., and Brown, S. A. (1982) The Natural Coumarins: Occurrence, Chemistry and Biochemistry, John Wiley & Sons, New York
  14. Berenbaum, M. R., and Zangerl, A. R. (1996) Phytochemistry 30, 1-24
  15. Cohen, M. B., Schuler, M. A., and Berenbaum, M. R. (1992) Proc. Natl. Acad. Sci. U. S. A. 89, 10920-10924[Abstract/Free Full Text]
  16. Hung, C.-F., Prapaipong, H., Berenbaum, M. R., and Schuler, M. A. (1995) Insect Biochem. Mol. Biol. 25, 89-99[CrossRef]
  17. Hung, C.-F., Berenbaum, M. R., and Schuler, M. A. (1997) Insect Biochem. Mol. Biol. 27, 377-385[CrossRef][Medline] [Order article via Infotrieve]
  18. Li, X., Zangerl, A. R., Schuler, M. A., and Berenbaum, M. R. (2000) J. Econ. Entomol. 93, 18-25[Medline] [Order article via Infotrieve]
  19. Li, W., Berenbaum, M. R., and Schuler, M. A. (2001) Insect Biochem. Mol. Biol. 31, 999-1011[CrossRef][Medline] [Order article via Infotrieve]
  20. Mao, W., Berhow, M., Zangerl, A. R., McGovern, J., and Berenbaum, M. R. (2006) J. Chem. Ecol. 32, 523-536[CrossRef][Medline] [Order article via Infotrieve]
  21. Mao, W., Rupasinghe, S., Zangerl, A., Schuler, M. A., and Berenbaum, M. R. (2006) Insect Mol. Biol. 15, 169-179[CrossRef][Medline] [Order article via Infotrieve]
  22. Li, W., Petersen, R. A., Schuler, M. A., and Berenbaum, M. R. (2002) Insect Mol. Biol. 11, 543-551[CrossRef][Medline] [Order article via Infotrieve]
  23. Li, X., Berenbaum, M. R., and Schuler, M. A. (2000) Insect Biochem. Mol. Biol. 30, 75-84[CrossRef][Medline] [Order article via Infotrieve]
  24. Li, X., Berenbaum, M. R., and Schuler, M. A. (2002) Insect Biochem. Mol. Biol. 32, 311-320[CrossRef][Medline] [Order article via Infotrieve]
  25. Chen, J.-S., Berenbaum, M. R., and Schuler, M. A. (2002) Insect Mol. Biol. 11, 175-186[CrossRef][Medline] [Order article via Infotrieve]
  26. Li, W., Schuler, M. A., and Berenbaum, M. R. (2003) Proc. Natl. Acad. Sci. U. S. A. 100, 14593-14598[Abstract/Free Full Text]
  27. Li, X., Baudry, J., Berenbaum, M. R., and Schuler, M. A. (2004) Proc. Natl. Acad. Sci. U. S. A. 101, 2939-2944[Abstract/Free Full Text]
  28. Wen, Z., Pan, L., Berenbaum, M. R., and Schuler, M. A. (2003) Insect Biochem. Mol. Biol. 33, 937-947[CrossRef][Medline] [Order article via Infotrieve]
  29. Sasabe, M., Wen, Z., Berenbaum, M. R., and Schuler, M. A. (2004) Gene (Amst.) 338, 163-175[CrossRef][Medline] [Order article via Infotrieve]
  30. Li, W., Zangerl, A. R., Schuler, M. A., and Berenbaum, M. R. (2004) Insect Mol. Biol. 13, 603-613[CrossRef][Medline] [Order article via Infotrieve]
  31. Baudry, J., Li, W., Pan, L., Berenbaum, M. R., and Schuler, M. A. (2003) Protein Eng. 16, 577-587[Abstract/Free Full Text]
  32. Wen, Z., Rupasinghe, S., Niu, G., Berenbaum, M. R., and Schuler, M. A. (2006) Mol. Biol. Evol. 23, 2434-2443[Abstract/Free Full Text]
  33. Humma, L. M., Farmerie, W. G., Wallace, M. R., and Johnson, J. A. (2000) BioTechniques 29, 970-972[Medline] [Order article via Infotrieve]
  34. Omura, T., and Sato, R. (1964) J. Biol. Chem. 239, 2379-2385[Free Full Text]
  35. Sevrioukova, I. F., Li, H., Zhang, H., Peterson, J. A., and Poulos, T. L. (1999) Proc. Natl. Acad. Sci. U. S. A. 96, 1863-1868[Abstract/Free Full Text]
  36. Wang, M., Roberts, D. L., Paschke, R., Shea, T. M., Masters, B. S. S., and Kim, J. J. P. (1997) Proc. Natl. Acad. Sci. U. S. A. 94, 8411-8416[Abstract/Free Full Text]
  37. Wade, R. C., Winn, P. J., Schlichting, I., and Sudarko (2004) J. Inorg. Biochem. 98, 1175-1182[CrossRef][Medline] [Order article via Infotrieve]
  38. Pan, L., Wen, Z., Baudry, J., Berenbaum, M. R., and Schuler, M. A. (2004) Arch. Biochem. Biophys. 422, 31-41[CrossRef][Medline] [Order article via Infotrieve]
  39. Wen, Z., Baudry, J., Berenbaum, M. R., and Schuler, M. A. (2005) Protein Eng. Des. Sel. 18, 191-199[Abstract/Free Full Text]
  40. Nadler, S. G., and Strobel, H. W. (1988) Arch. Biochem. Biophys. 261, 418-429[CrossRef][Medline] [Order article via Infotrieve]
  41. Lee-Robichaud, P., Akhtar, M. E., and Akhtar, M. (1998) Biochem. J. 332, 293-296[Medline] [Order article via Infotrieve]
  42. Lee-Robichaud, P., Akhtar, M. E., and Akhtar, M. (1999) Biochem. J. 342, 309-312[CrossRef][Medline] [Order article via Infotrieve]
  43. Bridges, A., Gruenke, L., Chang, Y. T., Vakser, I. A., Loew, G., and Waskell, L (1998) J. Biol. Chem. 273, 17036-17049[Abstract/Free Full Text]
  44. Lehnerer, M., Schulze, J., Achterhold, K., Lewis, D. F., and Hlavica, P. (2000) J. Biochem. (Tokyo) 127, 163-169[Abstract/Free Full Text]
  45. Shimizu, T., Tateishi, T., Hatano, M., and Fujii-Kuriyama, Y. (1991) J. Biol. Chem. 266, 3372-3375[Abstract/Free Full Text]
  46. Shen, S., and Strobel, H. W. (1992) Arch. Biochem. Biophys. 294, 83-90[CrossRef][Medline] [Order article via Infotrieve]
  47. Hlavica, P., Schulze, J., and Lewis, D. F. V. (2003) J. Inorg. Biochem. 96, 279-297[CrossRef][Medline] [Order article via Infotrieve]
  48. Nisimoto, Y., and Otsuka-Murakami, H. (1988) Biochemistry 27, 5869-5876[CrossRef][Medline] [Order article via Infotrieve]
  49. Shen, A. L., and Kasper, C. B. (1995) J. Biol. Chem. 270, 27475-27480[Abstract/Free Full Text]
  50. Graham-Lorence, S., and Peterson, J. A. (1996) FASEB J. 10, 206-214[Medline] [Order article via Infotrieve]
  51. Graham-Lorence, S., and Peterson, J. A. (1999) Arch. Biochem. Biophys. 369, 24-29[CrossRef][Medline] [Order article via Infotrieve]
  52. Peterson, J. A., and Graham-Lorence, S. (1998) Structure (Lond.) 6, 1079-1085[Medline] [Order article via Infotrieve]
  53. Shen, S., and Strobel, H. W. (1993) Arch. Biochem. Biophys. 304, 257-265[CrossRef][Medline] [Order article via Infotrieve]
  54. Berenbaum, M. R., and Zangerl, A. R. (1998) Proc. Natl. Acad. Sci. U. S. A. 95, 13743-13748[Abstract/Free Full Text]
  55. Zangerl, A. R., and Berenbaum, M. R. (2006) Ecology (Wash. D. C.) 87, 3070-3081

Add to CiteULike CiteULike   Add to Complore Complore   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us   Add to Digg Digg   Add to Reddit Reddit   Add to Technorati Technorati    What's this?


This article has been cited by other articles:


Home page
Plant Physiol.Home page
M. R. Berenbaum and A. R. Zangerl
Facing the Future of Plant-Insect Interaction Research: Le Retour a la "Raison d'Etre"
Plant Physiology, March 1, 2008; 146(3): 804 - 811.
[Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental Data
Right arrow All Versions of this Article:
282/14/10544    most recent
M607946200v1
Right arrow Submit a Letter to Editor
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowRequest Permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Mao, W.
Right arrow Articles by Schuler, M. A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Mao, W.
Right arrow Articles by Schuler, M. A.
Social Bookmarking
 Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Reddit   Add to Technorati  
What's this?


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 All ASBMB Journals   Molecular and Cellular Proteomics 
 Journal of Lipid Research   ASBMB Today 
Copyright © 2007 by the American Society for Biochemistry and Molecular Biology.
Advertisement
spacer
Advertisement
Advertisement