Gravity Persistent Signal 1 (GPS1) reveals novel cytochrome P450s involved in gravitropism.

PREMISE: Gravity is an important environmental factor that affects growth and development of plants. In response to changes in gravity, directional growth occurs along the major axes and lateral branches of both shoots and roots. The gravity persistent signal (gps) mutants of Arabidopsis thaliana were previously identified as having an altered response to gravity when reoriented relative to the gravity vector in the cold, with the gps1 mutant exhibiting a complete loss of tropic response under these conditions. METHODS: Thermal asymmetric interlaced (TAIL) PCR was used to identify the gene defective in gps1. Gene expression data, molecular modeling and computational substrate dockings, quantitative RT-PCR analyses, reporter gene fusions, and physiological analyses of knockout mutants were used to characterize the genes identified. RESULTS: Cloning of the gene defective in gps1 and genetic complementation revealed that GPS1 encodes CYP705A22, a cytochrome P450 monooxygenase (P450). CYP705A5, a closely related family member, was identified as expressed specifically in roots in response to gravistimulation, and a mutation affecting its expression resulted in a delayed gravity response, increased flavonol levels, and decreased basipetal auxin transport. Molecular modeling coupled with in silico substrate docking and diphenylboric acid 2-aminoethyl ester (DBPA) staining indicated that these P450s are involved in biosynthesis of flavonoids potentially involved in auxin transport. CONCLUSION: The characterization of two novel P450s (CYP705A22 and CYP705A5) and their role in the gravity response has offered new insights into the regulation of the genetic and physiological controls of plant gravitropism.

Rice P450 reductases differentially affect P450-mediated metabolism in bacterial expression systems.

We describe cloning and characterization of three rice (Oryza sativa) NADPH-cytochrome P450 reductases (OsCPRs; E.C.1.6.2.4) that are potential donors to plant P450s, including tryptamine 5-hydroxylase (T5H) in serotonin synthesis and cinnamate 4-hydroxylase (C4H) in phenylpropanoid synthesis. All three OsCPR transcripts are induced to varying degrees by stresses. Co-expression of full-length OsCPR1, OsCPR2 and OsCPR3 with either T5H or C4H in E. coli indicated that the OsCPR2/T5H and OsCPR2/C4H constructs displayed the highest T5H and C4H catalytic activities. The N-terminal residues of OsCPR2 were required for peak electron transfer activity to P450 even though deletion mutants with short N-terminal deletions were capable of reducing cytochrome c.

A gain-of-function polymorphism controlling complex traits and fitness in nature.

Identification of 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 quantitative trait locus 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 (branched-chain methionine allocation) loci controlling this variation and by two selectively favored amino acid changes in the glucosinolate-biosynthetic cytochrome P450 proteins 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.

VITAL NMR: using chemical shift derived secondary structure information for a limited set of amino acids to assess homology model accuracy.

Homology modeling is a powerful tool for predicting protein structures, whose success depends on obtaining a reasonable alignment between a given structural template and the protein sequence being analyzed. In order to leverage greater predictive power for proteins with few structural templates, we have developed a method to rank homology models based upon their compliance to secondary structure derived from experimental solid-state NMR (SSNMR) data. Such data is obtainable in a rapid manner by simple SSNMR experiments (e.g., (13)C-(13)C 2D correlation spectra). To test our homology model scoring procedure for various amino acid labeling schemes, we generated a library of 7,474 homology models for 22 protein targets culled from the TALOS+/SPARTA+ training set of protein structures. Using subsets of amino acids that are plausibly assigned by SSNMR, we discovered that pairs of the residues Val, Ile, Thr, Ala and Leu (VITAL) emulate an ideal dataset where all residues are site specifically assigned. Scoring the models with a predicted VITAL site-specific dataset and calculating secondary structure with the Chemical Shift Index resulted in a Pearson correlation coefficient (-0.75) commensurate to the control (-0.77), where secondary structure was scored site specifically for all amino acids (ALL 20) using STRIDE. This method promises to accelerate structure procurement by SSNMR for proteins with unknown folds through guiding the selection of remotely homologous protein templates and assessing model quality.

Crystal structure of a phenol-coupling P450 monooxygenase involved in teicoplanin biosynthesis.

The lipoglycopeptide antibiotic teicoplanin has proven efficacy against gram-positive pathogens. Teicoplanin is distinguished from the vancomycin-type glycopeptide antibiotics, by the presence of an additional cross-link between the aromatic amino acids 1 and 3 that is catalyzed by the cytochrome P450 monooxygenase Orf6* (CYP165D3). As a goal towards understanding the mechanism of this phenol-coupling reaction, we have characterized recombinant Orf6* and determined its crystal structure to 2.2-Å resolution. Although the structure of Orf6* reveals the core fold common to other P450 monooxygenases, there are subtle differences in the disposition of secondary structure elements near the active site cavity necessary to accommodate its complex heptapeptide substrate. Specifically, the orientation of the F and G helices in Orf6* results in a more closed active site than found in the vancomycin oxidative enzymes OxyB and OxyC. In addition, Met226 in the I helix replaces the more typical Gly/Ala residue that is positioned above the heme porphyrin ring, where it forms a hydrogen bond with a heme iron-bound water molecule. Sequence comparisons with other phenol-coupling P450 monooxygenases suggest that Met226 plays a role in determining the substrate regiospecificity of Orf6*. These features provide further insights into the mechanism of the cross-linking mechanisms that occur during glycopeptide antibiotics biosynthesis.

Citrus flavonoids luteolin, apigenin, and quercetin inhibit glycogen synthase kinase-3ß enzymatic activity by lowering the interaction energy within the binding cavity.

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.

A substrate-specific cytochrome P450 monooxygenase, CYP6AB11, from the polyphagous navel orangeworm (Amyelois transitella).

The navel orangeworm Amyelois transitella (Walker) (Lepidoptera: Pyralidae) is a serious pest of many tree crops in California orchards, including almonds, pistachios, walnuts and figs. To understand the molecular mechanisms underlying detoxification of phytochemicals, insecticides and mycotoxins by this species, full-length CYP6AB11 cDNA was isolated from larval midguts using RACE PCR. Phylogenetic analysis of this insect cytochrome P450 monooxygenase established its evolutionary relationship to a P450 that selectively metabolizes imperatorin (a linear furanocoumarin) and myristicin (a natural methylenedioxyphenyl compound) in another lepidopteran species. Metabolic assays conducted with baculovirus-expressed P450 protein, P450 reductase and cytochrome b(5) on 16 compounds, including phytochemicals, mycotoxins, and synthetic pesticides, indicated that CYP6AB11 efficiently metabolizes imperatorin (0.88 pmol/min/pmol P450) and slowly metabolizes piperonyl butoxide (0.11 pmol/min/pmol P450). LC-MS analysis indicated that the imperatorin metabolite is an epoxide generated by oxidation of the double bond in its extended isoprenyl side chain. Predictive structures for CYP6AB11 suggested that its catalytic site contains a doughnut-like constriction over the heme that excludes aromatic rings on substrates and allows only their extended side chains to access the catalytic site. CYP6AB11 can also metabolize the principal insecticide synergist piperonyl butoxide (PBO), a synthetic methylenedioxyphenyl compound, albeit slowly, which raises the possibility that resistance may evolve in this species after exposure to synergists under field conditions.

Arabidopsis CYP72C1 is an atypical cytochrome P450 that inactivates brassinosteroids.

Cytochrome P450 monooxygenases (P450s) are a diverse family of proteins that have specialized roles in secondary metabolism and in normal cell development. Two P450s in particular, CYP734A1 and CYP72C1, have been identified as brassinosteroid-inactivating enzymes important for steroid-mediated signal transduction in Arabidopsis thaliana. Genetic analyses have demonstrated that these P450s modulate growth throughout plant development. While members of the CYP734A subfamily inactivate brassinosteroids through C-26 hydroxylation, the biochemical activity of CYP72C1 is unknown. Because CYP734A1 and CYP72C1 in Arabidopsis diverge more than brassinosteroid inactivating P450s in other plants, this study examines the structure and biochemistry of each enzyme. Three-dimensional models were generated to examine the substrate binding site structures and determine how they might affect the function of each P450. These models have indicated that the active site of CYP72C1 does not contain several conserved amino acids typically needed for substrate hydroxylation. Heterologous expression of these P450s followed by substrate binding analyses have indicated that CYP734A1 binds active brassinosteroids, brassinolide and castasterone, as well as their upstream precursors whereas CYP72C1 binds precursors more effectively. Seedling growth assays have demonstrated that the genetic state of CYP734A1, but not CYP72C1, affected responsiveness to high levels of exogenous brassinolide supporting our observations that CYP72C1 acts on brassinolide precursors. Although there may be some overlap in their physiological function, the distinct biochemical functions of these proteins in Arabidopsis has significant potential to fine-tune the levels of different brassinosteroid hormones throughout plant growth and development.

Peptides from purified soybean beta-conglycinin inhibit fatty acid synthase by interaction with the thioesterase catalytic domain.

Fatty acid synthase (FAS) is uniquely expressed at high levels in cancer cells and adipose tissue. The objectives of this study were to identify, purify and validate soy FAS inhibitory peptides and to predict their binding modes. Soy peptides were isolated from hydrolysates of purified beta-conglycinin by co-immunoprecipitation and identified using LC-MS/MS. Three peptides, KNPQLR, EITPEKNPQLR and RKQEEDEDEEQQRE, inhibited FAS. The biological activity of these peptides was confirmed by their inhibitory activity against purified chicken FAS (IC(50) = 79, 27 and 16 mum, respectively) and a high correlation (r = -0.7) with lipid accumulation in 3T3-L1 adipocytes. The FAS inhibitory potency of soy peptides also correlated with their molecular mass, pI value and the number of negatively charged and hydrophilic residues. Molecular modeling predicted that the large FAS inhibitory peptides (EITPEKNPQLR and RKQEEDEDEEQQRE) bond to the thioesterase domain of human FAS with lower interaction energies (-442 and -353 kcal.mol(-1), respectively) than classical thioesterase inhibitors (Orlistat, -91 kcal.mol(-1) and C75, -51 kcal.mol(-1)). Docking studies suggested that soy peptides blocked the active site through interactions within the catalytic triad, the interface cavity and the hydrophobic groove in the human FAS thioesterase domain. FAS thioesterase inhibitory activities displayed by the synthetic soy peptides EITPEKNPQLR and RKQEEDEDEEQQRE (IC(50) = 10.1 +/- 1.6 and 10.7 +/- 4.4 mum, respectively) were higher than C75 (58.7 mum) but lower than Orlistat (0.9 mum). This is the first study to identify FAS inhibitory peptides from purified beta-conglycinin hydrolysates and predict their binding modes at the molecular level, leading to their possible use as nutraceuticals.

Quercetin-metabolizing CYP6AS enzymes of the pollinator Apis mellifera (Hymenoptera: Apidae).

Although the honey bee (Apis mellifera) genome contains far fewer cytochrome P450 genes associated with xenobiotic metabolism than other insect genomes sequenced to date, the CYP6AS subfamily, apparently unique to hymenopterans, has undergone an expansion relative to the genome of the jewel wasp (Nasonia vitripennis). The relative dominance of this family in the honey bee genome is suggestive of a role in processing phytochemicals encountered by honey bees in their relatively unusual diet of honey (comprising concentrated processed nectar of many plant species) and bee bread (a mixture of honey and pollen from many plant species). In this study, quercetin was initially suggested as a shared substrate for CYP6AS1, CYP6AS3, and CYP6AS4, by its presence in honey, extracts of which induce transcription of these three genes, and by in silico substrate predictions based on a molecular model of CYP6AS3. Biochemical assays with heterologously expressed CYP6AS1, CYP6AS3, CYP6AS4 and CYP6AS10 enzymes subsequently confirmed their activity toward this substrate. CYP6AS1, CYP6AS3, CYP6AS4 and CYP6AS10 metabolize quercetin at rates of 0.5+/-0.1, 0.5+/-0.1, 0.2+/-0.1, and 0.2+/-0.1 pmol quercetin/ pmol P450/min, respectively. Substrate dockings and sequence alignments revealed that the positively charged amino acids His107 and Lys217 and the carbonyl group of the backbone between Leu302 and Ala303 are essential for quercetin orientation in the CYP6AS3 catalytic site and its efficient metabolism. Multiple replacements in the catalytic site of CYP6AS4 and CYP6AS10 and repositioning of the quercetin molecule likely account for the lower metabolic activities of CYP6AS4 and CYP6AS10 compared to CYP6AS1 and CYP6AS3.

Aflatoxin B1 detoxification by CYP321A1 in Helicoverpa zea.

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.

Identification and characterization of topoisomerase II inhibitory peptides from soy protein hydrolysates.

Topoisomerases are targets of several anticancer agents because their inhibition impedes the processes of cell proliferation and differentiation in carcinogenesis. With very limited information available on the inhibitory activities of peptides derived from dietary proteins, the objectives of this study were to employ co-immunoprecipitation to identify inhibitory peptides in soy protein hydrolysates in a single step and to investigate their molecular interactions with topoisomerase II. For this, soy protein isolates were subjected to simulated gastrointestinal digestion with pepsin and pancreatin, and the human topoisomerase II inhibitory peptides were co-immunoprecipitated and identified on a CapLC- Micromass Q-TOF Ultima API system. The inhibitory activity of these peptides from soy isolates toward topoisomerase II was confirmed using three synthetic peptides, FEITPEKNPQ, IETWNPNNKP,and VFDGEL, which have IC 50 values of 2.4, 4.0, and 7.9 mM, respectively. The molecular interactions of these peptides evaluated by molecular docking revealed interaction energies with the topoisomerase II C-terminal domain (CTD) (-186 to -398 kcal/mol) that were smaller than for the ATPase domain (-169 to -357 kcal/mol) and that correlated well with our experimental IC 50 values ( R (2) = 0.99). In conclusion, three peptides released from in vitro gastrointestinal enzyme digestion of soy proteins inhibited human topoisomerase II activity through binding to the active site of the CTD domain.

Comparative molecular modeling of Anopheles gambiae CYP6Z1, a mosquito P450 capable of metabolizing DDT.

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.

Helicoverpa zea CYP6B8 and CYP321A1: different molecular solutions to the problem of metabolizing plant toxins and insecticides.

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.

High-yield expression and purification of isotopically labeled cytochrome P450 monooxygenases for solid-state NMR spectroscopy.

Cytochrome P450 monooxygenases (P450s), which represent the major group of drug metabolizing enzymes in humans, also catalyze important synthetic and detoxicative reactions in insects, plants and many microbes. Flexibilities in their catalytic sites and membrane associations are thought to play central roles in substrate binding and catalytic specificity. To date, Escherichia coli expression strategies for structural analysis of eukaryotic membrane-bound P450s by X-ray crystallography have necessitated full or partial removal of their N-terminal signal anchor domain and, often, replacement of residues more peripherally associated with the membrane (such as the F-G loop region). Even with these modifications, investigations of P450 structural flexibility remain challenging with multiple single crystal conditions needed to identify spatial variations between substrate-free and different substrate-bound forms. To overcome these limitations, we have developed methods for the efficient expression of 13C- and 15N-labeled P450s and analysis of their structures by magic-angle spinning solid-state NMR (SSNMR) spectroscopy. In the presence of co-expressed GroEL and GroES chaperones, full-length (53 kDa) Arabidopsis 13C,15N-labeled His4CYP98A3 is expressed at yields of 2-4 mg per liter of minimal media without the necessity of generating side chain modifications or N-terminal deletions. Precipitated His4CYP98A3 generates high quality SSNMR spectra consistent with a homogeneous, folded protein. These data highlight the potential of these methodologies to contribute to the structural analysis of membrane-bound proteins.

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.