|
Advertisement | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 282, Issue 14, 10544-10552, April 6, 2007
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*
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
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).
|
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
-turn in
-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 102119 in CYP6B1 and CYP6B4 and positions 104121 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 |
|---|
|
|
|---|
cDNA Library SelectionThe 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 CellsThe 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 810% 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 ActivitiesReaction 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. RatiosFive 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 RatesTo 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 MetabolitesImperatorin 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 ModelingThe 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).
|
| RESULTS |
|---|
|
|
|---|
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).
|
|
|
Homology ModelingTo 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.
|
| DISCUSSION |
|---|
|
|
|---|
|
2 (2) and
1 (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
2-
3 loop and the
1-
2 loop and N terminus of
4.
|
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
13 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 |
|---|
The on-line version of this article (available at http://www.jbc.org) contains supplemental DNA sequence comparisons. ![]()
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. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
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] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| All ASBMB Journals | Molecular and Cellular Proteomics |
| Journal of Lipid Research | ASBMB Today |