Identifying the causal genes that control complex trait variation remains challenging, limiting our appreciation of the evolutionary processes that influence polymorphisms in nature. We cloned a QTL that controls plant defensive chemistry, damage by insect herbivores, survival, and reproduction in the natural environments where this polymorphism evolved. These ecological effects are driven by duplications in the BCMA loci controlling this QTL and by two selectively favored amino acid changes in the glucosinolate-biosynthetic P450s that they encode. These changes cause a gain of novel enzyme function, modulated by allelic differences in catalytic rate and gene copy number. Ecological interactions in diverse environments likely contribute to the widespread polymorphism of this biochemical function.

One of the challenges faced in malarial control is the acquisition of insecticide resistance that has developed in mosquitoes that are vectors for this disease. Anopheles gambiae, which has been the major mosquito vector of the malaria parasite Plasmodium falciparum in Africa, has over the years developed resistance to insecticides including dieldrin, 1,1-bis(p-chlorophenyl)-2,2,2-trichloroethane (DDT), and pyrethroids. Previous microarray studies using fragments of 230 An. gambiae genes identified five P450 loci, including CYP4C27, CYP4H15, CYP6Z1, CYP6Z2, and CYP12F1, that showed significantly higher expression in the DDT-resistant ZAN/U strain compared with the DDT-susceptible Kisumu strain. To predict whether either of the CYP6Z1 and CYP6Z2 proteins might potentially metabolize DDT, we generated and compared molecular models of these two proteins with and without DDT docked in their catalytic sites. This comparison indicated that, although these two CYP6Z proteins share high sequence identity, their metabolic profiles were likely to differ dramatically from the larger catalytic site of CYP6Z1, potentially involved in DDT metabolism, and the more constrained catalytic site of CYP6Z2, not likely to metabolize DDT. Heterologous expressions of these proteins have corroborated these predictions: only CYP6Z1 is capable of metabolizing DDT. Overlays of these models indicate that slight differences in the backbone of SRS1 and variations of side chains in SRS2 and SRS4 account for the significant differences in their catalytic site volumes and DDT-metabolic capacities. These data identify CYP6Z1 as one important target for inhibitor design aimed at inactivating insecticide-metabolizing P450s in natural populations of this malarial mosquito.cytochrome P450 monooxygenases ͉ insecticides ͉ plant allelochemicals I n 2004, the World Health Organization reported that up to 2.7 million people die of malaria every year with 80-90% of these deaths occurring in Africa (ref. 1 and www.africanfront.com/ AIDS1.php). Many prevention and treatment strategies have been developed to tackle this life-threatening disease from the side of the mosquito vector and that of the human host (2, 3). These range from antimalarial drugs to indoor spraying of insecticides and use of bed nets treated with pyrethroid insecticides. Although these practices have helped reduce human mortality, various issues have emerged with mosquito vectors developing insecticide resistance and parasites developing drug resistance. Both of these significantly reduce the efficacy of current malarial prevention and treatment practices (2, 4-8).The development of insecticide resistance in mosquito vectors is illustrated by Anopheles gambiae, which has been the major mosquito vector of the malaria parasite Plasmodium falciparum in Africa (4, 6). Throughout the years, people have reported Anopheles resistance (albeit low-level resistance) to various insecticides including dieldrin (a cyclodiene-type insecticide), 1,1-bis(pchlorophenyl)-2,2,2-trichloroethane (DDT), and pe...

No abstract

Under continual exposure to naturally occurring plant toxins and synthetic insecticides, insects have evolved cytochrome P450 monooxygenases (P450s) capable of metabolizing a wide range of structurally different compounds. Two such P450s, CYP6B8 and CYP321A1, expressed in Helicoverpa zea (a lepidopteran) in response to plant allelochemicals and plant signaling molecules metabolize these compounds with varying efficiencies. While sequence alignments of these proteins indicate highly divergent substrate recognition sites (SRSs), homology models developed for them indicate that the two active site cavities have essentially the same volume with distinct shapes dictated by side-chain differences in SRS1 and SRS5. CYP6B8 has a narrower active site cavity extending from substrate access channel pw2a with a very narrow access to the ferryl oxygen atom. This predicted shape suggests that bulkier molecules bind further from the ferryl oxygen at positions that are not as effectively metabolized. In contrast, CYP321A1 is predicted to have a more spacious cavity allowing larger molecules to access the heme-bound oxygen. The metabolic profiles for several plant toxins (xanthotoxin, angelicin) and insecticides (cypermethrin, aldrin and diazinon) correlate well with these predictive models. The absence of Thr in the I helix of CYP321A1 and hydroxyl groups on many of its substrates suggests that this insect P450 mediates oxygen activation by a mechanism different from that employed by CYP107A1 and CYP158A1, which are two bacterial P450s also lacking Thr in their I helix, and most other P450s that contain Thr in their I helix.

Despite extensive primary sequence diversity, crystal structures of several bacterial cytochrome P450 monooxygenases (P450s) and a single eukaryotic P450 indicate that these enzymes share a structural core of alpha-helices and beta-sheets and vary in the loop regions contacting individual substrates. To determine the extent to which individual structural features are conserved among divergent P450s existing in a single biosynthetic pathway, we have modeled the structures of four highly divergent P450s (CYP73A5, CYP84A1, CYP75B1, CYP98A3) in the Arabidopsis phenylpropanoid pathway synthesizing lignins, flavonoids and anthocyanins. Analysis of these models has indicated that, despite primary sequence identities as low as 13%, the structural cores and several loop regions of these P450s are highly conserved. Substrate docking indicated that all four enzymes employ a common strategy to identify their substrates in that their cinnamate-derived substrates align along helix I with their aromatic ring positioned towards the C-terminus of this helix and their aliphatic tails positioned towards the N-terminus. Further similarity was observed in the way the substrates contact the consensus P450 substrate recognition sites (SRS). Residues predicted to contact the aromatic ring region exist in SRS5, SRS6 and the C-terminal portion of SRS4 and residues contacting the distal end of each substrate exist in SRS1, SRS2 and the N-terminal portion of SRS4. Alignments of the regions contacting the aromatic ring region indicate that SRS4, SRS5 and SRS6 share higher degrees of sequence conservation than found in SRS1, SRS2 or the full-length protein.

Pancreatic cancer studies have shown that inhibition of glycogen synthase kinase-3β (GSK-3β) leads to decreased cancer cell proliferation and survival by abrogating nuclear factor κB (NFκB) activity. In this investigation, various citrus compounds, including flavonoids, phenolic acids, and limonoids, were individually investigated for their inhibitory effects on GSK-3β by using a luminescence assay. Of the 22 citrus compounds tested, the flavonoids luteolin, apigenin, and quercetin had the highest inhibitory effects on GSK-3β, with 50% inhibitory values of 1.5, 1.9, and 2.0 μM, respectively. Molecular dockings were then performed to determine the potential interactions of each citrus flavonoid with GSK-3β. Luteolin, apigenin, and quercetin were predicted to fit within the binding pocket of GSK-3β with low interaction energies (-76.4, -76.1, and -84.6 kcal·mol(-1), respectively) and low complex energies (-718.1, -688.1, and -719.7 kcal·mol(-1), respectively). Our results indicate that several citrus flavonoids inhibit GSK-3β activity and suggest that these have potential to suppress the growth of pancreatic tumors.

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.

The polyphagous corn earworm Helicoverpa zea frequently encounters aflatoxins, mycotoxins produced by the pathogens Aspergillus flavus and A. parasiticus, which infect many of this herbivore's host plants. While aflatoxin B1 metabolism by midgut enzymes isolated from fifth instars feeding on control diets was not detected, this compound was metabolized by midgut enzymes isolated from larvae consuming diets supplemented with xanthotoxin, coumarin, or indole-3-carbinol, phytochemicals that are likely to co-occur with aflatoxin in infected host plants. Of the two metabolites generated, the main derivative identified in midguts induced with these chemicals and in reactions containing heterologously expressed CYP321A1 was aflatoxin P1 (AFP1), an O-demethylated product of AFB1. RT-PCR gel blots indicated that the magnitude of CYP321A1 transcript induction by these chemicals is associated with the magnitude of increase in the metabolic activities of induced midgut enzymes (coumarin>xanthotoxin>indole 3-carbinol). These results indicate that induction of P450s, such as CYP321A1, plays an important role in reducing AFB1 toxicity to H. zea. Docking of AFB1 in the molecular models of CYP321A1 and CYP6B8 highlights differences in their proximal catalytic site volumes that allow only CYP321A1 to generate the AFP1 metabolite.