Ile115Leu mutation in the SRS1 region of an insect cytochrome P450 (CYP6B1) compromises substrate turnover via changes in a predicted product release channel

The Chromosome-Level Genome of Hestina assimilis (Lepidoptera: Nymphalidae) Reveals the Evolution of Saprophagy-Related Genes in Brush-Footed Butterflies

Most butterflies feed on nectar, while some saprophagous butterflies forage on various non-nectar foods. To date, little is known about the genomic and molecular shifts associated with the evolution of the saprophagous feeding strategy. Here, we assembled the high-quality chromosome-level genome of Hestina assimilis to explore its saprophagous molecular and genetic mechanisms. This chromosome-level genome of H. assimilis is 412.82 Mb, with a scaffold N50 of 15.70 Mb. In total, 98.11% of contigs were anchored to 30 chromosomes. Compared with H. assimilis and other Nymphalidae butterflies, the genes of metabolism and detoxification experienced expansions. We annotated 80 cytochrome P450 (CYP) genes in the H. assimilis genome, among which genes belonging to the CYP4 subfamily were significantly expanded (p < 0.01). These P450 genes were unevenly distributed and mainly concentrated on chromosomes 6-9. We identified 33 olfactory receptor (OR), 20 odorant-binding protein (OBP), and six gustatory receptor (GR) genes in the H. assimilis genome, which were fewer than in the nectarivorous Danaus plexippus. A decreased number of OBP, OR, and GR genes implied that H. assimilis should resort less to olfaction and gustation than their nectarivorous counterparts, which need highly specialized olfactory and gustatory functions. Moreover, we found one site under positive selection occurred in residue 996 (phenylalanine) of GR genes exclusive to H. assimilis, which is conservative in most lineages. Our study provides support for the adaptive evolution of feeding habits in butterflies.

Single amino acid variations drive functional divergence of cytochrome P450s in Helicoverpa species

Divergence of gene function is a hallmark of evolution, but assessing such divergence in one species or between species requires information on functional alterations of the alleles and homologs. Here, we explore the functional divergence of two paralogs, CYP6AE19 and CYP6AE20, from Helicoverpa armigera, and two close orthologs, CYP6B8 and CYP6B7, from two related species (Helicoverpa zea and H. armigera); although there is high sequence identity within each pair of enzymes, the latter P450 of each pair has lost metabolic competence towards the plant allelochemical xanthotoxin. Multiple chimeric and single/double site mutants were created by exchanging the diverse substrate recognition sites (SRSs) and amino acids within each pair of P450s. Heterologous expression in Sf9 cells and in vitro metabolism studies showed that the exchange of SRS4 swapped the activity of CYP6AE19 and CYP6AE20, and subsequent site-directed mutagenesis demonstrated that the CYP6AE20 V318M substitution causes a gain-of-function towards xanthotoxin. Meanwhile, a single amino acid substitution (L489P) in SRS6 was found to swap activity between the CYP6B orthologs. Sequence alignments of CYP6AE paralogs and all reported insect xanthotoxin-metabolizing P450s suggest M318 and P489 are essential for the catalytic activities of CYP6AE paralogs and CYP6B orthologs, respectively, but P450s in different subfamilies may have different mechanisms towards the same substrate. Our findings demonstrate that a single amino acid substitution can suffice to alter substrate metabolism and this functional divergence resulting from natural mutations will help to further our understanding of the process of natural selection of P450 genes and their role in insect-host plant interactions.

Xenobiotic responses in insects

Insects have evolved a powerful detoxification system to protect themselves against environmental and anthropogenic xenobiotics including pesticides and nanoparticles. The resulting tolerance to insecticides is an immense problem in agriculture. In this study, we summarize advances in our understanding of insect xenobiotic responses: the detoxification strategies and the regulation mechanisms against xenobiotics including nanoparticles, the problem of response specificity and the potential usefulness of this study field for an elaborate pest management. In particular, we highlight that versatility of the detoxification system relies on the relatively unspecific recognition of a broad range of potential toxic substances that trigger either of various canonical xenobiotic responses signaling pathways, including CncC/Keap1, HR96, AHR/ARNT, GPCR, and MAPK/CREB. However, it has emerged that the actual response to an inducer may nevertheless be specific. There are two nonexclusive possibilities that may explain response specificity: (1) differential cross-talk between the known pathways and (2) additional, yet unidentified regulators and pathways of detoxification. Hence, a deeper and broader understanding of the regulation mechanisms of xenobiotic response in insects in the future might facilitate the development and application of highly efficient and environmentally friendly pest control methods, allowing us to face the challenge of the world population growth.

The overexpression and variant of CYP6G4 associated with propoxur resistance in the housefly, Musca domestica L

BACKGROUND

The control of the housefly, Musca domestica, heavily relies on the application of insecticides. Propoxur, a carbamate, has been widely used for vector control. The housefly populations with high propoxur resistance display point mutations and overexpression of acetylcholinesterase. However, the roles of cytochrome P450 monoxygenases (P450s), as important detoxification enzymes, remain poorly understand in the housefly resistant to propoxur.

RESULTS

P450s were implied to contribute to propoxur resistance based on the synergism of piperonyl butoxide (PBO) and the increase of P450 enzyme activity in the near-isogenic line propoxur resistant strain (N-PRS). Five P450 genes (CYP6G4, CYP6A25, CYP304A1, CYP6D3, and CYP6A1) by RNA-sequencing comparison were significantly up-regulated in the N-PRS strain with >1035-fold resistance to propoxur. A total of 13 non-synonymous mutations of three P450 genes (CYP6G4, CYP6D3, and CYP6D8) were found in the N-PRS strain. The amino acid substitutions of CYP6D3 and CYP6D8 were probably not resistance-associated single nucleotide polymorphisms (SNPs) because they were also found in the aabys susceptible strain. However, CYP6G4 variant in the N-PRS strain was not found in the aabys strain. The conjoint analysis of mutations and a series of genetic crosses exhibited that the housefly propoxur resistance was strongly associated with the mutations of CYP6G4 gene.

CONCLUSION

Our results suggested that a combination of up-regulated transcript levels and mutations of CYP6G4 contributed to propoxur resistance in the housefly. © 2021 Society of Chemical Industry.

Saltational evolution of a pesticide-metabolizing cytochrome P450 in a global crop pest

BACKGROUND

The cotton bollworm, Helicoverpa armigera (Hübner), is a damaging insect pest threatening agricultural crops worldwide as a result of its resistance to insecticides. Metabolic resistance to pyrethroid insecticides is conferred by the chimeric P450 enzyme CYP337B3, produced by unequal crossing-over between CYP337B1 and CYP337B2. CYP337B3 is 99.7% similar to CYP337B1 except for the 177 N-terminal amino acids (AAs) containing the substrate recognition site 1 from CYP337B2. Here, we studied the structure-function relationship of CYP337B3 and CYP337B1 to determine the AAs that enable CYP337B3 to efficiently hydroxylate the 4'-carbon position of fenvalerate, which neither CYP337B1 nor CYP337B2 can do.

RESULTS

Site-directed mutagenesis showed that the L114F substitution in CYP337B3 reduced its 4'-hydroxylation activity by 89%, but the reciprocal F114L substitution in CYP337B1 increased its 4'-hydroxylation activity to only 49% of the level of CYP337B3. Docking models showed that AA 114 seems to have different functions in CYP337B1 and CYP337B3. Antibodies detected two- to three-fold more CYP337B1 than CYP337B3 in larval cuticle, which along with a 49% 4'-hydroxylation activity increase due to a F114L substitution in vivo might be expected to provide as much protection for the larva against exposure to fenvalerate as CYP337B3. However, CYP337B3 is present at much higher frequencies than CYP337B1-CYP337B2 in most populations, including those recently invading South America.

CONCLUSION

The metabolic resistance to pyrethroids in H. armigera has evolved by saltational evolution - by a single mutation, an unequal crossing-over, producing a larger selective advantage than could be attained gradually by stepwise improvement of the parental enzyme. © 2021 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Dual Substrate Specificity of the Rutinosidase from Aspergillus niger and the Role of Its Substrate Tunnel

Rutinosidases (α-l-rhamnopyranosyl-(1-6)-β-d-glucopyranosidases, EC 3.2.1.168, CAZy GH5) are diglycosidases that cleave the glycosidic bond between the disaccharide rutinose and the respective aglycone. Similar to many retaining glycosidases, rutinosidases can also transfer the rutinosyl moiety onto acceptors with a free -OH group (so-called transglycosylation). The recombinant rutinosidase from Aspergillus niger (AnRut) is selectively produced in Pichia pastoris. It can catalyze transglycosylation reactions as an unpurified preparation directly from cultivation. This enzyme exhibits catalytic activity towards two substrates; in addition to rutinosidase activity, it also exhibits β-d-glucopyranosidase activity. As a result, new compounds are formed by β-glucosylation or rutinosylation of acceptors such as alcohols or strong inorganic nucleophiles (NaN₃). Transglycosylation products with aliphatic aglycones are resistant towards cleavage by rutinosidase, therefore, their side hydrolysis does not occur, allowing higher transglycosylation yields. Fourteen compounds were synthesized by glucosylation or rutinosylation of selected acceptors. The products were isolated and structurally characterized. Interactions between the transglycosylation products and the recombinant AnRut were analyzed by molecular modeling. We revealed the role of a substrate tunnel in the structure of AnRut, which explained the unusual catalytic properties of this glycosidase and its specific transglycosylation potential. AnRut is attractive for biosynthetic applications, especially for the use of inexpensive substrates (rutin and isoquercitrin).

Cytochrome P450s Cyp4p1 and Cyp4p2 associated with the DDT tolerance in the Drosophila melanogaster strain 91-R

Cytochrome P450s are part of a super-gene family that has undergone gene duplication, divergence, over-expression and, in some cases, loss of function. One such case is the 91-R and 91-C strains of common origin, in Drosophila melanogaster, whereby 91-R (DDT resistant strain) overexpresses Cyp4p1 and Cyp4p2 and both genes are lost in 91-C (DDT susceptible strain). In this study, we used a comparative approach to demonstrate that transcription of Cyp4p1 and Cyp4p2 were constitutively up-regulated in the Drosophila melanogaster strain 91-R as compared to another DDT susceptible strain Canton-S which does not have a loss of function of these genes. Furthermore, significantly increased expression of Cyp4p1 and Cyp4p2 was induced in 91-R in response to sublethal DDT exposure, however, such induction did not occur in the DDT treated Canton-S. Additionally, fixed nucleotide variation within putative transcription factor binding sites of Cyp4p1 and Cyp4p2 promoters were observed between 91-R and Canton-S, however, their impact on transcription remains to be determined. Two GAL4/UAS transgenic strains with integrated heat shock-inducible Cyp4p1- or Cyp4p2-RNAi constructs within wild-type genetic backgrounds were developed. Following heat shock induction of Cyp4p1 and Cyp4p2 knockdown, these transgenic lines showed increased DDT mortality as compared to their corresponding non-heat shock controls. These results provide a functional link of Cyp4p1 and Cyp4p2 in conferring tolerance to DDT exposure.

Engineering of lysine cyclodeaminase conformational dynamics for relieving substrate and product inhibitions in the biosynthesis of l-pipecolic acid

Efficient biocatalytic process construction by relieving substrate and product inhibitions via identification and engineering of enzyme conformational dynamics.

Comparative CYP-omic analysis between the DDT-susceptible and -resistant Drosophila melanogaster strains 91-C and 91-R

BACKGROUND

Cytochrome P450 monooxygenases (P450s) are involved in the biosynthesis of endogenous intracellular compounds and the metabolism of xenobiotics, including chemical insecticides. We investigated the structural and expression level variance across all P450 genes with respect to the evolution of insecticide resistance under multigenerational dichlorodiphenyltrichloroethane (DDT) selection.

RESULTS

RNA-sequencing (RNA-seq) and reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) indicated that the transcript levels of seven P450 genes were significantly up-regulated and three P450 genes were down-regulated in the DDT-resistant strain 91-R, as compared to the control strain 91-C. The overexpression of Cyp6g1 was associated with the presence of an Accord and an HMS-Beagle element insertion in the 5' upstream region in conjunction with copy number variation in the 91-R strain, but not in the 91-C strain. A total of 122 (50.2%) fixed nonsynonymous (amino acid-changing) mutations were found between 91-C and 91-R, and 20 (8.2%) resulted in amino acid changes within functional domains. Three P450 proteins were truncated as a result of premature stop codons and fixed between strains.

CONCLUSION

Our results demonstrate that a combination of changes in P450 protein-coding regions and transcript levels are possibly associated with DDT resistance, and thereby suggest that selection for variant function may occur within this gene family in response to chronic DDT exposure. © 2018 Society of Chemical Industry.

Ligand Access Channels in Cytochrome P450 Enzymes: A Review

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

Phylogenetic and functional characterization of ten P450 genes from the CYP6AE subfamily of Helicoverpa armigera involved in xenobiotic metabolism

The cotton bollworm, Helicoverpa armigera, is a generalist herbivore widely distributed over the world and is a major lepidopteran pest on cotton. Studies, especially from Asia, show that it relies on cytochrome P450 monooxygenases with broad substrate specificities to protect itself from pesticides. The number of P450s may have expanded in the processes of coping with the wide diversity of phytochemicals that the insect encounters among its numerous host plants. In order to examine the metabolic capabilities of these P450s, we focused here on all ten P450s of the Helicoverpa armigera CYP6AE subfamily, which can be easily induced by plant toxins and pyrethroids. These P450s, along with cytochrome P450 reductase, were heterologously expressed in insect cells and compared functionally. In vitro metabolism showed that all CYP6AE subfamily members can convert esfenvalerate to 4'-hydroxyesfenvalerate efficiently except CYP6AE20. In contrast, none of the recombinant CYP6AE enzymes could metabolise gossypol under our experimental conditions. Epoxidation capabilities were observed in the CYP6AE subfamily, aldrin can be converted to dieldrin at rates up to 0.45 ± 0.04 pmol/min/pmol P450. Seven P450s in this subfamily can metabolise imidacloprid, but with lower efficiency than Bemisia tabaci CYP6CM1vQ. CYP6AE20 had virtually no metabolic competence to these four compounds but could metabolise several model fluorogenic substrates. These results showed the broad substrate spectrum of H. armigera CYP6AE P450s and suggest a limited role of gossypol upon the evolution of H. armigera CYP6AE genes.

Genomic exaptation enables Lasius niger adaptation to urban environments

Background

The world is rapidly urbanizing, and only a subset of species are able to succeed in stressful city environments. Efficient genome-enabled stress response appears to be a likely prerequisite for urban adaptation. Despite the important role ants play in the ecosytem, only the genomes of ~13 have been sequenced so far. Here, we present the draft genome assembly of the black garden ant Lasius niger – the most successful urban inhabitant of all ants – and we compare it with the genomes of other ant species, including the closely related Camponotus floridanus.

Results

Sequences from 272 M Illumina reads were assembled into 41,406 contigs with total length of 245 MB, and N50 of 16,382 bp, similar to other ant genome assemblies enabling comparative genomic analysis. Remarkably, the predicted proteome of L. niger is significantly enriched relative to other ant genomes in terms of abundance of domains involved in nucleic acid binding, DNA repair, and nucleotidyl transferase activity, reflecting transposable element proliferation and a likely genomic response. With respect to environmental stress, we note a proliferation of various detoxification genes, including glutatione-S-transferases and those in the cytochrome P450 families. Notably, the CYP9 family is highly expanded with 19 complete and 21 nearly complete members - over twice as many compared to other ants. This family exhibits the signatures of strong directional selection, with eleven positively selected positions in ligand-binding pockets of enzymes. Gene family contraction was detected for several components of the olfactory system, accompanied by instances of both directional selection and relaxation.

Conclusions

Our results suggest that the success of L. niger in urbanized areas may be the result of fortuitous coincidence of several factors, including the expansion of the CYP9 cytochrome family due to coevolution with parasitic fungi, the diversification of DNA repair systems as an answer to proliferation of retroelements, and the reduction of olfactory system and behavioral preadaptations from non-territorial subdominant life strategies found in natural environments. Diversification of cytochromes and DNA repair systems along with reduced odorant communication are the basis of L. niger pollutant resistance and polyphagy, while non-territorial and mobilization strategies allows more efficient exploitation of large but patchy food sources.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0867-x) contains supplementary material, which is available to authorized users.

Enzyme Tunnels and Gates As Relevant Targets in Drug Design

Many enzymes contain tunnels and gates that are essential to their function. Gates reversibly switch between open and closed conformations and thereby control the traffic of small molecules-substrates, products, ions, and solvent molecules-into and out of the enzyme's structure via molecular tunnels. Many transient tunnels and gates undoubtedly remain to be identified, and their functional roles and utility as potential drug targets have received comparatively little attention. Here, we describe a set of general concepts relating to the structural properties, function, and classification of these interesting structural features. In addition, we highlight the potential of enzyme tunnels and gates as targets for the binding of small molecules. The different types of binding that are possible and the potential pharmacological benefits of such targeting are discussed. Twelve examples of ligands bound to the tunnels and/or gates of clinically relevant enzymes are used to illustrate the different binding modes and to explain some new strategies for drug design. Such strategies could potentially help to overcome some of the problems facing medicinal chemists and lead to the discovery of more effective drugs.

Engineering of an epoxide hydrolase for efficient bioresolution of bulky pharmaco substrates

Optically pure epoxides are essential chiral precursors for the production of (S)-propranolol, (S)-alprenolol, and other β-adrenergic receptor blocking drugs. Although the enzymatic production of these bulky epoxides has proven difficult, here we report a method to effectively improve the activity of BmEH, an epoxide hydrolase from Bacillus megaterium ECU1001 toward α-naphthyl glycidyl ether, the precursor of (S)-propranolol, by eliminating the steric hindrance near the potential product-release site. Using X-ray crystallography, mass spectrum, and molecular dynamics calculations, we have identified an active tunnel for substrate access and product release of this enzyme. The crystal structures revealed that there is an independent product-release site in BmEH that was not included in other reported epoxide hydrolase structures. By alanine scanning, two mutants, F128A and M145A, targeted to expand the potential product-release site displayed 42 and 25 times higher activities toward α-naphthyl glycidyl ether than the wild-type enzyme, respectively. These results show great promise for structure-based rational design in improving the catalytic efficiency of industrial enzymes for bulky substrates.

CAVER 3.0: a tool for the analysis of transport pathways in dynamic protein structures

Tunnels and channels facilitate the transport of small molecules, ions and water solvent in a large variety of proteins. Characteristics of individual transport pathways, including their geometry, physico-chemical properties and dynamics are instrumental for understanding of structure-function relationships of these proteins, for the design of new inhibitors and construction of improved biocatalysts. CAVER is a software tool widely used for the identification and characterization of transport pathways in static macromolecular structures. Herein we present a new version of CAVER enabling automatic analysis of tunnels and channels in large ensembles of protein conformations. CAVER 3.0 implements new algorithms for the calculation and clustering of pathways. A trajectory from a molecular dynamics simulation serves as the typical input, while detailed characteristics and summary statistics of the time evolution of individual pathways are provided in the outputs. To illustrate the capabilities of CAVER 3.0, the tool was applied for the analysis of molecular dynamics simulation of the microbial enzyme haloalkane dehalogenase DhaA. CAVER 3.0 safely identified and reliably estimated the importance of all previously published DhaA tunnels, including the tunnels closed in DhaA crystal structures. Obtained results clearly demonstrate that analysis of molecular dynamics simulation is essential for the estimation of pathway characteristics and elucidation of the structural basis of the tunnel gating. CAVER 3.0 paves the way for the study of important biochemical phenomena in the area of molecular transport, molecular recognition and enzymatic catalysis. The software is freely available as a multiplatform command-line application at http://www.caver.cz.

Cytochrome P450 2B diversity and dietary novelty in the herbivorous, desert woodrat (Neotoma lepida)

Detoxification enzymes play a key role in plant-herbivore interactions, contributing to the on-going evolution of ecosystem functional diversity. Mammalian detoxification systems have been well studied by the medical and pharmacological industries to understand human drug metabolism; however, little is known of the mechanisms employed by wild herbivores to metabolize toxic plant secondary compounds. Using a wild rodent herbivore, the desert woodrat (Neotoma lepida), we investigated genomic structural variation, sequence variability, and expression patterns in a multigene subfamily involved in xenobiotic metabolism, cytochrome P450 2B (CYP2B). We hypothesized that differences in CYP2B expression and sequence diversity could explain differential abilities of woodrat populations to consume native plant toxins. Woodrats from two distinct populations were fed diets supplemented with either juniper (Juniperus osteosperma) or creosote bush (Larrea tridentata), plants consumed by woodrats in their respective desert habitats. We used Southern blot and quantitative PCR to determine that the genomic copy number of CYP2B in both populations was equivalent, and similar in number to known rodent copy number. We compared CYP2B expression patterns and sequence diversity using cloned hepatic CYP2B cDNA. The resulting sequences were very diverse, and clustered into four major clades by amino acid similarity. Sequences from the experimental treatments were distributed non-randomly across a CYP2B tree, indicating unique expression patterns from woodrats on different diets and from different habitats. Furthermore, within each major CYP2B clade, sequences shared a unique combination of amino acid residues at 13 sites throughout the protein known to be important for CYP2B enzyme function, implying differences in the function of each major CYP2B variant. This work is the most comprehensive investigation of the genetic diversity of a detoxification enzyme subfamily in a wild mammalian herbivore, and contributes an initial genetic framework to our understanding of how a wild herbivore responds to critical changes in its diet.

Active-site hydration and water diffusion in cytochrome P450cam: a highly dynamic process

Long-timescale molecular dynamics simulations (300 ns) are performed on both the apo- (i.e., camphor-free) and camphor-bound cytochrome P450cam (CYP101). Water diffusion into and out of the protein active site is observed without biased sampling methods. During the course of the molecular dynamics simulation, an average of 6.4 water molecules is observed in the camphor-binding site of the apo form, compared to zero water molecules in the binding site of the substrate-bound form, in agreement with the number of water molecules observed in crystal structures of the same species. However, as many as 12 water molecules can be present at a given time in the camphor-binding region of the active site in the case of apo-P450cam, revealing a highly dynamic process for hydration of the protein active site, with water molecules exchanging rapidly with the bulk solvent. Water molecules are also found to exchange locations frequently inside the active site, preferentially clustering in regions surrounding the water molecules observed in the crystal structure. Potential-of-mean-force calculations identify thermodynamically favored trans-protein pathways for the diffusion of water molecules between the protein active site and the bulk solvent. Binding of camphor in the active site modifies the free-energy landscape of P450cam channels toward favoring the diffusion of water molecules out of the protein active site.

Homology modelling of Drosophila cytochrome P450 enzymes associated with insecticide resistance

BACKGROUND

Overexpression of the cytochrome P450 gene Cyp6g1 confers resistance against DDT and a broad range of other insecticides in Drosophila melanogaster Meig. In the absence of crystal structures of CYP6G1 or complexes with its substrates, structural studies rely on homology modelling and ligand docking to understand P450-substrate interactions.

RESULTS

Homology models are presented for CYP6G1, a P450 associated with resistance to DDT and neonicotinoids, and two other enzymes associated with insecticide resistance in D. melanogaster, CYP12D1 and CYP6A2. The models are based on a template of the X-ray structure of the phylogenetically related human CYP3A4, which is known for its broad substrate specificity. The model of CYP6G1 has a much smaller active site cavity than the template. The cavity is also 'V'-shaped and is lined with hydrophobic residues, showing high shape and chemical complementarity with the molecular characteristics of DDT. Comparison of the DDT-CYP6G1 complex and a non-resistant CYP6A2 homology model implies that tight-fit recognition of this insecticide is important in CYP6G1. The active site can accommodate differently shaped substrates ranging from imidacloprid to malathion but not the pyrethroids permethrin and cyfluthrin.

CONCLUSION

The CYP6G1, CYP12D1 and CYP6A2 homology models can provide a structural insight into insecticide resistance in flies overexpressing P450 enzymes with broad substrate specificities.

Decreased detoxification genes and genome size make the human body louse an efficient model to study xenobiotic metabolism

The human body louse, Pediculus humanus humanus, has one of the smallest insect genomes, containing ∼10 775 annotated genes. Annotation of detoxification [cytochrome P450 monooxygenase (P450), glutathione-S-transferase (GST), esterase (Est) and ATP-binding cassette transporter (ABC transporter)] genes revealed that they are dramatically reduced in P. h. humanus compared to other insects except for Apis mellifera. There are 37 P450, 13 GST and 17 Est genes present in P. h. humanus, approximately half the number found in Drosophila melanogaster and Anopheles gambiae. The number of putatively functional ABC transporter genes in P. h. humanus and Ap. mellifera are the same (36) but both have fewer than An. gambiae (44) or Dr. melanogaster (65). The reduction of detoxification genes in P. h. humanus may be a result of this louse's simple life history, in which it does not encounter a wide variety of xenobiotics. Neuronal component genes are highly conserved across different insect species as expected because of their critical function. Although reduced in number, P. h. humanus still retains at least a minimum repertoire of genes known to confer metabolic or toxicokinetic resistance to xenobiotics (eg Cyp3 clade P450s, Delta GSTs, B clade Ests and B/C subfamily ABC transporters), suggestive of its high potential for resistance development.

Metabolic enzymes associated with xenobiotic and chemosensory responses in Nasonia vitripennis

The numbers of glutathione S-transferase, cytochrome P450 and esterase genes in the genome of the hymenopteran parasitoid Nasonia vitripennis are about twice those found in the genome of another hymenopteran, the honeybee Apis mellifera. Some of the difference is associated with clades of these families implicated in xenobiotic resistance in other insects and some is in clades implicated in hormone and pheromone metabolism. The data support the hypothesis that the eusocial behaviour of the honeybee and the concomitant homeostasis of the nest environment may obviate the need for as many gene/enzyme systems associated with xenobiotic metabolism as are found in other species, including N. vitripennis, that are thought to encounter a wider range of potentially toxic xenobiotics in their diet and habitat.

Pathways and mechanisms for product release in the engineered haloalkane dehalogenases explored using classical and random acceleration molecular dynamics simulations

Eight mutants of the DhaA haloalkane dehalogenase carrying mutations at the residues lining two tunnels, previously observed by protein X-ray crystallography, were constructed and biochemically characterized. The mutants showed distinct catalytic efficiencies with the halogenated substrate 1,2,3-trichloropropane. Release pathways for the two dehalogenation products, 2,3-dichloropropane-1-ol and the chloride ion, and exchange pathways for water molecules, were studied using classical and random acceleration molecular dynamics simulations. Five different pathways, denoted p1, p2a, p2b, p2c, and p3, were identified. The individual pathways showed differing selectivity for the products: the chloride ion releases solely through p1, whereas the alcohol releases through all five pathways. Water molecules play a crucial role for release of both products by breakage of their hydrogen-bonding interactions with the active-site residues and shielding the charged chloride ion during its passage through a hydrophobic tunnel. Exchange of the chloride ions, the alcohol product, and the waters between the buried active site and the bulk solvent can be realized by three different mechanisms: (i) passage through a permanent tunnel, (ii) passage through a transient tunnel, and (iii) migration through a protein matrix. We demonstrate that the accessibility of the pathways and the mechanisms of ligand exchange were modified by mutations. Insertion of bulky aromatic residues in the tunnel corresponding to pathway p1 leads to reduced accessibility to the ligands and a change in mechanism of opening from permanent to transient. We propose that engineering the accessibility of tunnels and the mechanisms of ligand exchange is a powerful strategy for modification of the functional properties of enzymes with buried active sites.

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

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 (V(max) of 10.02 pmol/min/pmol of cytochrome P450 monooxygenase (P450)) significantly higher than CYP6AB3v1 (V(max) 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 Val(92) on the proximal face of CYP6AB3v1 with Ala(92) 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.

The evolutionary ecology of insect resistance to plant chemicals

Understanding the diversity of insect responses to chemical pressures (e.g. plant allelochemicals and pesticides) in their local ecological context represents a key challenge in developing durable pest control strategies. To what extent do the resistance mechanisms evolved by insects to deal with the chemical defences of plants differ from those that have evolved to resist insecticides? Here, we review recent advances in our understanding of insect resistance to plant chemicals, with a special emphasis on their underlying molecular basis, evaluate costs associated with each resistance trait, and discuss the ecological and evolutionary significance of these findings.

A deficit of detoxification enzymes: pesticide sensitivity and environmental response in the honeybee

The honeybee genome has substantially fewer protein coding genes ( approximately 11 000 genes) than Drosophila melanogaster ( approximately 13 500) and Anopheles gambiae ( approximately 14 000). Some of the most marked differences occur in three superfamilies encoding xenobiotic detoxifying enzymes. Specifically there are only about half as many glutathione-S-transferases (GSTs), cytochrome P450 monooxygenases (P450s) and carboxyl/cholinesterases (CCEs) in the honeybee. This includes 10-fold or greater shortfalls in the numbers of Delta and Epsilon GSTs and CYP4 P450s, members of which clades have been recurrently associated with insecticide resistance in other species. These shortfalls may contribute to the sensitivity of the honeybee to insecticides. On the other hand there are some recent radiations in CYP6, CYP9 and certain CCE clades in A. mellifera that could be associated with the evolution of the hormonal and chemosensory processes underpinning its highly organized eusociality.

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.

Inhibition of CYP6B1-mediated detoxification of xanthotoxin by plant allelochemicals in the black swallowtail (Papilio polyxenes)

The structural and biosynthetic diversity of allelochemicals in plants is thought to arise from selection for additive toxicity as a consequence of toxin mixture or for enhanced toxicity as a result of synergism. In order to understand how insects cope with this type of plant defense, we tested the effects of some allelochemicals in host plants of the black swallowtail Papilio polyxenes on the xanthotoxin-metabolic activity of CYP6B1, the principal enzyme responsible for the detoxification of furanocoumarins in this caterpillar. Additionally, the effects of some synthetic compounds not normally encountered by P. polyxenes on CYP6B1 were tested. These studies demonstrate that the integrity of furanocoumarin structure is important for competitive binding to the active site of CYP6B1, even though the carbonyl group on the pyranone ring apparently does not affect its inhibitory capacity, as in the case of furanochromones. Angular furanocoumarins are generally less phototoxic to many organisms than linear furanocoumarins due to their reduced capacity for cross-linking DNA strands, yet they are more toxic than linear furanocoumarins to black swallowtail larvae. This enhanced toxicity in vivo may be due to the ability of angular furanocoumarins to bind to the active site of CYP6B1 without being rapidly metabolized. This binding reduces the availability of CYP6B1 to metabolize other linear furanocoumarins. The structure-activity relationships for methylenedioxyphenyl compounds, flavonoids, imidazole, and imidazole derivatives are also discussed in light of their capacity to inhibit the xanthotoxin-metabolic activity of CYP6B1.

Multiple molecular recognition mechanisms. Cytochrome P450--a case study

Biomolecular recognition is complex. The balance between the different molecular properties that contribute to molecular recognition, such as shape, electrostatics, dynamics and entropy, varies from case to case. This, along with the extent of experimental characterization, influences the choice of appropriate computational approaches to study biomolecular interactions. Here, we present computational studies of cytochrome P450 enzymes and their interactions with small molecules and with other proteins. These interactions exemplify some of the diversity of molecular determinants of binding affinity and specificity observed for proteins and we discuss some of the challenges that they pose for molecular modelling and simulation.

Regulation of an insect cytochrome P450 monooxygenase gene (CYP6B1) by aryl hydrocarbon and xanthotoxin response cascades

Many organisms respond to toxic compounds in their environment by inducing regulatory networks controlling the expression and activity of cytochrome P450 monooxygenase (P450s) detoxificative enzymes. In particular, black swallowtail (Papilio polyxenes) caterpillars respond to xanthotoxin, a toxic phytochemical in their hostplants, by activating transcription of the CYP6B1 promoter via several regions located within 150 nt of the transcription initiation site. One such element is the xenobiotic response element to xanthotoxin (XRE-Xan) that lies upstream of consensus XRE-AhR (xenobiotic response element to the aryl hydrocarbon receptor) and OCT-1 (octamer-1 binding site) element known to be utilized in mammalian aryl hydrocarbon response cascades. Two-plasmid transfections conducted in Sf9 cells have indicated that XRE-Xan, XRE-AhR and a number of other proximal elements, but not OCT-1, are critical for basal as well as xanthotoxin- and benzo[alpha]pyrene-induced transcription of the CYP6B1 promoter. Four-plasmid transfections with vectors co-expressing the spineless (Ss) and tango (Tgo) proteins, the Drosophila melanogaster homologues of mammalian AhR and ARNT, have indicated that these proteins enhance basal expression of the CYP6B1 promoter but not the magnitude of its xanthotoxin and benzo[alpha]pyrene induction. Based on these results, we propose that these Drosophila transcription factors modulate basal expression of this promoter in a ligand-independent manner and attenuate its subsequent responses to planar aryl hydrocarbons (benzo[alpha]pyrene) and allelochemicals (xanthotoxin).