Molecular analysis of CYP321A1, a novel cytochrome P450 involved in metabolism of plant allelochemicals (furanocoumarins) and insecticides (cypermethrin) in Helicoverpa zea

Cytochrome P450 monooxygenases play a significant role in the detoxification of hostplant allelochemicals and synthetic insecticides in Lepidoptera. In the corn earworm Helicoverpa zea, a noctuid of considerable economic importance, metabolisms of xanthotoxin, a toxic furanocoumarin, and alpha-cypermethrin, an insecticide, are mediated by at least one P450 with a catalytic site capable of accepting both substrates. To further the characterization of P450s in this species, we have cloned three full-length cDNAs encoding two CYP4M subfamily members and a novel CYP321A subfamily member. RNA analyses have demonstrated that the CYP321A1 gene is highly induced (51-fold) in larval midguts in response to xanthotoxin but not cypermethrin. Both CYP4M genes are expressed at negligible levels that are not increased by xanthotoxin or cypermethrin. Baculovirus-mediated expression of the full-length CYP321A1 cDNA has demonstrated that the CYP321A1 protein metabolizes xanthotoxin and angelicin, like the CYP6B1 protein in the furanocoumarin specialist Papilio polyxenes, and alpha-cypermethrin, like the CYP6B8 protein previously characterized in H. zea. In contrast, the CYP4M7 protein does not metabolize xanthotoxin at any detectable level. We conclude that at least two xanthotoxin-inducible P450s from highly divergent subfamilies (CYP6B and CYP321A) contribute to the resistance of H. zea larvae to toxic furanocoumarins and insecticides. Genomic PCR analysis indicates that the CYP321A1 gene has evolved independently from the CYP6B genes known to be present in this insect.     

Cross-resistance to alpha-cypermethrin after xanthotoxin ingestion in Helicoverpa zea (Lepidoptera: Noctuidae)

Cytochrome P450 monooxygenases (P450) are membrane-bound hemoproteins that play important roles in conferring protection against both naturally occurring phytochemicals and synthetic organic insecticides. Despite the potential for common modes of detoxification, cross-resistance between phytochemicals and synthetic organic insecticides has rarely been documented. In this study, we examined the responses of a susceptible strain of corn earworm, Helicoverpa zea (Boddie), a polyphagous noctuid, to exposure by an allelochemical infrequently encountered in its host plants and by an insecticide widely used for control purposes. Within a single generation, survivors of xanthotoxin exposure displayed higher levels of tolerance to alpha-cypermethrin than did unexposed control larvae. The F1 offspring of xanthotoxin-exposed survivors also displayed higher alpha-cypermethrin tolerance than did offspring of unexposed control larvae, suggesting that increased alpha-cypermethrin tolerance after xanthotoxin exposure represents, at least in part, heritable resistance. Administration of piperonyl butoxide, a P450 synergist, demonstrated that resistance to both xanthotoxin and alpha-cypermethrin is P450-mediated. Alpha-cypermethrin-exposed survivors, however, failed to show superior growth on xanthotoxin diets. Assays with control larvae, larvae induced by both xanthotoxin and alpha-cypermethrin, and survivors of LD50 doses of both compounds indicated that H. zea midgut P450s are capable of metabolizing both xanthotoxin and alpha-cypermethrin. Metabolism of each compound is significantly inhibited by the presence of the other compound, suggesting that at least one form of P450 in H. zea midguts degrades both compounds and may constitute the biochemical basis for possible cross-resistance. Compared with control larvae, xanthotoxin- and alpha-cypermethrin-induced larvae displayed 2- to 4-fold higher P450-mediated metabolism of both compounds. However, xanthotoxin- and alpha-cypermethrin-exposed survivors exhibited much higher (2.5- to 11-fold) metabolism of both compounds than did the induced larvae. The metabolism results, like the bioassay results, are consistent with the interpretation that increased alpha-cypermethrin tolerance after xanthotoxin exposure is attributable mainly to heritable resistance.     

Molecular docking of substrates and inhibitors in the catalytic site of CYP6B1, an insect cytochrome p450 monooxygenase

Furanocoumarins represent plant toxins that are used in the treatment of a variety of skin diseases and are metabolized by cytochrome p450 monooxygenases (p450s) existing in insects such as Papilio polyxenes (the black swallowtail). To elucidate the active site in the CYP6B1 protein that is the principal p450 existing in this species, we have constructed a homology model of it based on sequence and structure alignments with the bacterial CYP102 protein whose crystal structure has been defined and with the insect CYP6B4 protein that also metabolizes furanocoumarins. In the derived CYP6B1 model, Phe116 and His117 in SRS1, Phe371 in SRS5 and Phe484 in SRS6 contribute to the formation of a resonant network that stabilizes the p450's catalytic site and allows for interactions with its furanocoumarin substrates. The first two of these residues are absolutely conserved in all members of the insect CYP6B subfamily and the last two are variable in different members of the CYP6B subfamily. A combination of theoretical and experimental docking analyses of two substrates (xanthotoxin and bergapten) and two inhibitors (coumarin and pilocarpine) of this p450 provide significant information on the positioning of furanocoumarins within this catalytic pocket. Molecular replacement models based on the results of variations at two of these critical amino acids provide support for our furanocoumarin-docked model and begin to rationalize the altered substrate reactivities observed in experimental analyses.     

Metabolism of linear and angular furanocoumarins by Papilio polyxenes CYP6B1 co-expressed with NADPH cytochrome P450 reductase

One challenge in the heterologous expression of microsomal cytochrome P450 monooxygenases (P450s) is fulfilling their obligatory requirement for electrons transferred from NADPH P450 reductase. We have established co-expression parameters for Papilio polyxenes CYP6B1 and house fly P450 reductase in baculovirus-infected Sf9 cells that allow for efficient expression of both components and significantly enhance metabolic turnover of this insect P450's substrates. These expression conditions have allowed us to reexamine the turnover capacities of CYP6B1 toward linear and angular furanocoumarins present in the host plants for the specialist caterpillar P. polyxenes. Coexpression of CYP6B1 and P450 reductase at equivalent viral concentrations [MOI (multiplicity of infection) ratio of 1] results in turnover rates for the linear furanocoumarins xanthotoxin and psoralen, which are increased 32-33 fold over the turnover rates obtained with CYP6B1 expressed alone. The turnover rate for the angular furanocoumarin angelicin is also significantly increased to 4.76 nmol/min/nmol P450 compared to its barely detectable level obtained with CYP6B1 expressed alone. Substrate binding analyses indicate that all three of these compounds elicit typical type I binding spectra but with varying magnitudes and affinities that are indicative of each substrate's effectiveness at coordinating with the heme iron. The relative proportions of high spin state generated with these substrates are consistent with CYP6B1 metabolizing these furanocoumarins in the rank order xanthotoxin>psoralen>angelicin. These differential activities for CYP6B1 suggest that it may have been an ancient participant in the coevolutionary arms race between papilionid butterflies and their apiaceous host plants. Due to its ability to handle a range of furanocoumarin structures, CYP6B1 may have contributed to P. polyxenes' early colonization of linear furanocoumarin-containing plants and to its subsequent colonization of angular furanocoumarin-containing plants.     

Cytochrome P450 and actin genes expressed in Helicoverpa zea and Helicoverpa armigera: paralogy/orthology identification, gene conversion and evolution.

Molecular phylogenetic analysis was conducted using conserved cytoplasmic actin and diversified cytochrome P450 (P450) sequences isolated from Helicoverpa zea and Helicoverpa armigera, two species thought to be closely related based on allozyme analyses. These sequences were compared in turn with published sequences from other insects to gain insight into how different gene families evolve. In Bombyx mori and these Helicoverpa species, cytoplasmic actin genes are present as a pair of tandemly duplicated paralogs with coding sequence identities as high as 95.5% (B. mori), 98.9% (H. zea) and 98.5% (H. armigera) due to recent 5'-polar gene conversions. Phylogeny and interspecies comparisons assign the six actin genes into two orthologous groups: HaA3a/HzA3a/BmA3 and HaA3b/HzA3b/BmA4, which exhibit more similarities between H. zea and H. armigera than between Helicoverpa species and B. mori. Like the actin genes in H. zea, four CYP6B genes exist as two pairs of duplicated paralogs with recent 5'-polar gene conversions. Interspecific comparisons and phylogeny analysis identified three groups of orthologous CYP6B genes: H. zea CYP6B8 or CYP6B28/H. armigera CYP6B7, H. zea CYP6B27/H. armigera CYP6B6, and H. zea CYP6B9/H. armigera CYP6B2/Heliothis virescens CYP6B10. The low degree of divergence in the first two of these groups is comparable to allelic variation within a single species. These orthologous relationships and the high degrees of similarity in both actin and P450 genes strongly indicate that these Helicoverpa species are extremely closely related.     

CYP6B1 and CYP6B3 of the black swallowtail (Papilio polyxenes): adaptive evolution through subfunctionalization

Gene duplication provides essential material for functional divergence of proteins and hence allows organisms to adapt to changing environments. Following duplication events, redundant paralogs may undergo different evolutionary paths via processes known as nonfunctionalization, neofunctionalization, or subfunctionalization. Studies of adaptive evolution at the molecular level have progressed rapidly by computationally analyzing nucleotide substitution patterns but such studies are limited by the absence of information relating to alterations of function of the encoded enzymes. In this respect, evolution of the Papilio polyxenes cytochrome P450 monooxygenases (P450s) responsible for the adaptation of this insect to furanocoumarin-containing host plants provides an excellent model for elucidating the evolutionary fate of duplicated genes. Evidence from sequence and functional analysis in combination with molecular modeling indicates that the paralogous CYP6B1 and CYP6B3 genes in P. polyxenes have probably evolved via subfunctionalization after the duplication event by which they arose. Both enzymes have been under independent purifying selection as evidenced by the low dN/dS ratio in both the coding region and substrate recognition sites. Both enzymes have maintained their ability to metabolize linear and angular furanocoumarins albeit at different efficiencies. Comparisons of molecular models developed for the CYP6B3 and CYP6B1 proteins highlight differences in their binding modes that account for their different activities toward linear and angular furanocoumarins. That P. polyxenes maintains these 2 furanocoumarin-metabolizing loci with somewhat different activities and expression patterns provides this species with the potential to acquire P450s with novel functions while maintaining those most critical to its exclusive feeding on its current range of host plants.     

Cloning and characterization of two cytochrome P450 CYP6AX1 and CYP6AY1 cDNAs from Nilaparvata lugens Stål (Homoptera: Delphacidae)

Two full‐length P450 cDNAs, CYP6AX1 and CYP6AY1, were cloned from the brown planthopper Nilaparvata lugens Stål (Homoptera: Delphacidae). Both CYP6AX1 and CYP6AY1 are typical microsomal P450s and their deduced amino acid sequences share common characteristics with other members of the insect P450 CYP6 family. CYP6AX1 and CYP6AY1 show the highest percent identity (36%) of amino acid to each other; both of them have 31–33% amino acid identity with CYP6B1 from Papilio polyxenes (Lepidoptera: Papilionidae), CYP6B4 from Papilio glaucus (Lepidoptera: Papilionidae), and CYP6B8 from Helicoverpa zea (Lepidoptera: Noctuidae). Phylogenetic analysis showed the clustering of CYP6AX1 and CYP6AY1 was in the clade including CYP6AE1 from Depressaria pastinacella (Lepidoptera: Oecophoridae) and the CYP6B family members from Helicoverpa and Papilio species. Northern blot analysis revealed that both of the P450s were induced by the resistant rice variety B5 (Oryza sativa L), and CYP6AY1 was expressed at a higher level than CYP6AX1. The results suggest that more than one P450s are likely involved in metabolism of rice allelochemicals and that they are possibly important components in adaptation of Nilaparvata lugens to host rice. Arch. Insect Biochem. Physiol. 64:88–99, 2007. © 2007 Wiley‐Liss, Inc.     

Identification of variable amino acids in the SRS1 region of CYP6B1 modulating furanocoumarin metabolism

The homology model of Papilio polyxenes CYP6B1 places Ile115, one of two variable amino acids, in the SRS1 of various CYP6B subfamily proteins in close proximity to the heme and Ala113, another variable amino acid, in a more distal position. We have constructed mutant CYP6B1 proteins altered at either of these positions and homology models of each based on multiple alignments with crystallized P450 proteins. The homology models suggest the existence of significant structural diversity in the hydrogen bond network surrounding the heme as a result of single point mutations in SRS1. Mutagenesis of Ile115 or Ala113 to other residues present in the insect CYP6B subfamily indicates that these amino acids control the spin state of the heme and, as a result, the catalytic activity of this monooxygenase. In particular, the I115L mutation significantly increases the spin state of the heme coordinately with 2- to 4-fold increases in its turnover of linear furanocoumarins. Other A113V, A113L, A113Q, and A113E mutations display more variation in their effects but, in each case, strong correlations exist between furanocoumarin turnover and heme spin state. These data demonstrate that variable amino acids in SRS1 of the insect CYP6B subfamily exert dramatic effects on the range of furanocoumarins metabolized, even when they occur in positions potentially distal from the substrate. These effects are possibly mediated through rearrangement of the local hydrogen bond network.     

Permethrin Induction of Multiple Cytochrome P450 Genes in Insecticide Resistant Mosquitoes,Culex quinquefasciatus

The expression of some insect P450 genes can be induced by both exogenous and endogenous compounds and there is evidence to suggest that multiple constitutively overexpressed P450 genes are co-responsible for the development of resistance to permethrin in resistant mosquitoes. This study characterized the permethrin induction profiles of P450 genes known to be constitutively overexpressed in resistant mosquitoes, Culex quinquefasciatus. The gene expression in 7 of the 19 P450 genes CYP325K3v1, CYP4D42v2, CYP9J45, (CYP) CPIJ000926, CYP325G4, CYP4C38, CYP4H40 in the HAmCqG8 strain, increased more than 2-fold after exposure to permethrin at an LC50 concentration (10 ppm) compared to their acetone treated counterpart; no significant differences in the expression of these P450 genes in susceptible S-Lab mosquitoes were observed after permethrin treatment. Eleven of the fourteen P450 genes overexpressed in the MAmCqG6 strain, CYP9M10, CYP6Z12, CYP9J33, CYP9J43, CYP9J34, CYP306A1, CYP6Z15, CYP9J45, CYPPAL1, CYP4C52v1, CYP9J39, were also induced more than doubled after exposure to an LC50 (0.7 ppm) dose of permethrin. No significant induction in P450 gene expression was observed in the susceptible S-Lab mosquitoes after permethrin treatment except for CYP6Z15 and CYP9J39, suggesting that permethrin induction of these two P450 genes are common to both susceptible and resistant mosquitoes while the induction of the others are specific to insecticide resistant mosquitoes. These results demonstrate that multiple P450 genes are co-up-regulated in insecticide resistant mosquitoes through both constitutive overexpression and induction mechanisms, providing additional support for their involvement in the detoxification of insecticides and the development of insecticide resistance.     

Modifications in the N-terminus of an insect cytochrome P450 enhance production of catalytically active protein in baculovirus-Sf9 cell expression systems

Although baculovirus vectors are powerful tools for the heterologous expression of proteins in insect cell cultures, some insect and plant microsomal P450 proteins are not effectively expressed in this system. Hypothesizing that their expression failures might result from collisions between their N-terminal sequences and adjacent cytosolic sequences, we compared and mutated the N-terminus of Papilio multicaudatus CYP6B33, which is inappropriately folded in Sf9 cells, to sequences present in its Papilio polyxenes CYP6B1 counterpart, which is efficiently expressed and appropriately folded. Molecular modeling of the three differences in the linker separating the signal anchor domain (SAD) and the cytosolic domain identified Val32 in CYP6B33 as a residue potentially important for folding and/or positioning of the cytosolic domain. Mutation of Val32 to Ala32 in the CYP6B33 linker (CYP6B33 V32A mutant) or replacement of the CYP6B33 SAD with that of CYP6B1 (CYP6B1 1-20/CYP6B33 21-500 mutant) allowed for significant P450 expression, indicating that complex interactions involving both the signal anchor and membrane linker affect folding and activity of P450s in this heterologous expression system.     

MOLECULAR CLONING AND SEQUENCE ANALYSIS OF THREE NEW FULL‐LENGTH cDNAs OF CYTOCHROME P450 (CYP6N3v1‐v3)*FROM AEDES ALBOPICTUS**

Abstract The three new full‐length cDNA sequences including the complete 5′‐and 3′‐ untranslated regions (UTR) coding for cytochrome P450s from Aedes albopictus have been obtained. The P450 proteins deduced from the nucleotide sequences shared 58.6% ‐ 62.4% amino acid identity with CYP6N1 and CYP6N2 from Anopheles gambiae, and 99% with each other. The three new complete sequences have been submitted and named as CYP6N3v1, CYP6N3v2 and CYP6N3v3 by the P450 Nomenclature Committee. The original cDNAs were obtained by rapid amplification of cDNA ends (RACE) approach with several pairs of gene specific primers based on the cDNA fragment previously obtained from deltamethrin‐resistant strain of Ae. albopictus. Further analysis showed that the three new sequences are present in both resistant strain and susceptible strain and might be effectively translated. In addition, the 5′‐ and 3′‐UTRs were compared between the CYP6N3vl‐v3 and other known insect P450s. The multiplicity of trans‐lational control of insect P450 genes was discussed.