Expression of human cytochromes P450 1A1 and P450 1A2 as fused enzymes with yeast NADPH-cytochrome P450 oxidoreductase in transgenic tobacco plants

Among 11 isoforms of the human cytochrome P450 enzymes metabolizing xenobiotics, CYP 1A1 and CYP 1A2 were major P450 species in the metabolism of the herbicides chlortoluron and atrazine in a yeast expression system. CYP1A2 was more active in the metabolism of both herbicides than CYP1A1. The fused enzymes of CYP1A1 and CYP1A2 with yeast NADPH-cytochrome P450 oxidoreductase were functionally active in the microsomal fraction of the yeast Saccharomyces cerevisiae and showed increased specific activity towards 7-ethoxyresorufin as compared to CYP1A1 and CYP1A2 alone. Then, both fused enzymes were each expressed in the microsomes of tobacco (Nicotiana tabacum cv. Samsun NN) plants. The transgenic plants expressing the CYP1A2 fusion enzyme had higher resistance to the herbicide chlortoluron than the plants expressing the CYP1A1 fusion enzyme did. The transgenic plants expressing the CYP1A2 fused enzyme metabolized chlortoluron to a larger extent to its non-phytotoxic metabolites through N-demethylation and ring-methyl hydroxylation as compared to the plants expressing the CYP1A1 fused enzyme. Thus, the possibility of increasing the herbicide resistance in the transgenic plants by the selection of P450 species and the fusion with P450 reductase is discussed.     

Exploring the structure characteristics and major channels of cytochrome P450 2A6, 2A13, and 2E1 with pilocarpine

The majority of cytochromes P450 play a critical role in metabolism of endogenous and exogenous substrates, some of its products are carcinogens. Therefore, inhibition of P450 enzymes activity can promote the detoxification and elimination of chemical carcinogens. In this study, molecular dynamics (MD) simulations and adaptive steered molecular dynamics (ASMD) simulations were performed to explore the structure features and channel dynamics of three P450 isoforms 2A6, 2A13, and 2E1 bound with the common inhibitor pilocarpine. The binding free energy results combined with the PMF calculations give a reasonable ranking of binding affinity, which are consistent with the experimental data. Our results uncover how a sequence divergence of different CYP2 enzymes causes individual variations in major channel selections. On the basis of channel bottleneck and energy decomposition analysis, we propose a gating mechanism of their respective major channels in three enzymes, which may be attributed to a reversal of Phe209 in CYP2A6/2A13, as well as the rotation of Phe116 and Phe298 in CYP2E1. The hydrophobic residues not only make strong hydrophobic interactions with inhibitor, but also act as gatekeeper to regulate the opening of channel. The present study provides important insights into the structure–function relationships of three cytochrome P450s and the molecular basis for development of potent inhibitors.     

细胞色素P450与除草剂抗性转基因植物

介绍了除草剂代谢有关的细胞色素P450基因及其应用,已从动植物体中分离具有除草剂代谢活性的细胞色素P450基因,通过转基因方法,成功培育出抗除草剂的转基因植物。     

Plant cytochrome P450-mediated herbicide metabolism

In the last two decades it has become apparent that enzymes of the P450 monooxygenase (P450) superfamily are responsible for the Phase I metabolism of numerous herbicides representing several classes of organic compounds. The majority of experimental evidence for P450 involvement in herbicide metabolism has been derived from in vitro studies in which the catalytic activity of plant microsomes towards herbicidal substrates was measured in the presence of various P450 inhibitors and activators. While the studies with microsomes elicited much appreciation for the pivotal roles of plant P450s in herbicide metabolism, detailed characterization of these enzymes only became possible after the isolation of genes encoding specific isoforms responsible for herbicide conversion. Several lines of evidence suggest that the development of herbicide resistance in weeds by enhanced detoxification is frequently associated with elevated levels of P450 activity. Enhanced detoxification-based herbicide resistance is particularly difficult to control, because it can involve resistance to multiple, chemically unrelated classes of herbicides. Continued research efforts are aimed at elucidating the role of P450s in the metabolic fates of herbicides in plants and the development of herbicide resistance in weeds. Recent advances made in the isolation and genetic manipulation of P450 enzymes have created new opportunities for their application in engineering herbicide tolerance and bioremediation.     

分子伴侣HdeA与底物蛋白SurA作用机制的模拟研究

分子伴侣HdeA与底物蛋白间的相互作用可帮助底物蛋白复性,这是肠道致病菌得以在酸性环境中幸存的重要原因之一.为探究HdeA发挥伴侣活性的作用机制,本研究采用分子对接和分子动力学的方法,模拟了HdeA与底物蛋白SurA间的相互作用,计算了二者的结合自由能.通过分析HdeA-SurA复合物体系的作用模式、氢键作用以及能量分解的结果,确定了HdeA与底物蛋白SurA结合时发挥重要作用的关键氨基酸残基.该研究结果为以后采用实验手段探究HdeA与底物蛋白之间的作用提供了重要的理论参考,同时为今后设计与开发HdeA的抑制剂提供了理论指导依据.     

Molecular basis of the recognition of arachidonic acid by cytochrome P450 2E1 along major access tunnel

Cytochrome P450 2E1 is widely known for its ability to oxidize both low molecular weight xenobiotics and endogenous fatty acids (e.g., arachidonic acid (AA)). In this study, we investigated the structural features of the AA‐bound CYP2E1 complex utilizing molecular dynamics (MD) and found that the distinct binding modes for both AA and fatty acid analog are conserved. Moreover, multiple random acceleration MD simulations and steered MD simulations uncovered the most possible tunnel for fatty acids. The main attractions are derived from three key residues, His107, Ala108, and His109, whose side chains reorient to keep ligands bound via hydrogen bonds during the initial unbinding process. More importantly, based on the calculated binding free energy results, we hypothesize that the hydrogen bonds between the receptor and the ligand are the most important contributors involved in the binding affinity. Thus, it is inferred that the hydrogen bonds between these three residues and the ligand may help offer insights into the structural basis of the different ligand egress mechanisms for fatty acids and small weight compounds. Our investigation provides detailed atomistic insights into the structural features of human CYP2E1–fatty acid complex structures. Furthermore, the ligand‐binding characteristics obtained in the present study are helpful for both experimental and computational studies of CYPs and may allow future researchers to achieve desirable changes in enzymatic activities. © 2014 Wiley Periodicals, Inc. Biopolymers 103: 53–66, 2015.     

Herbicide resistance of transgenic rice plants expressing human CYP1A1

Cytochrome P450 monooxygenases (P450s) metabolize herbicides to produce mainly non-phytotoxic metabolites. Although rice plants endogenously express multiple P450 enzymes, transgenic plants expressing other P450 isoforms might show improved herbicide resistance or reduce herbicide residues. Mammalian P450s metabolizing xenobiotics are reported to show a broad and overlapping substrate specificity towards lipophilic foreign chemicals, including herbicides. These P450s are ideal for enhancing xenobiotic metabolism in plants. A human P450, CYP1A1, metabolizes various herbicides with different structures and modes of herbicide action. We introduced human CYP1A1 into rice plants, and the transgenic rice plants showed broad cross-resistance towards various herbicides and metabolized them. The introduced CYP1A1 enhanced the metabolism of chlorotoluron and norflurazon. The herbicides were metabolized more rapidly in the transgenic rice plants than in non-transgenic controls. Transgenic rice plants expressing P450 might be useful for reducing concentrations of various chemicals in the environment.     

Resistance cost of a cytochrome P450 herbicide metabolism mechanism but not an ACCase target site mutation in a multiple resistant Lolium rigidum population

Costs of resistance are predicted to reduce plant productivity in herbicide-resistant weeds. Lolium rigidum herbicide-susceptible individuals (S), individuals possessing cytochrome P450-based herbicide metabolism (P450) and multiple resistant individuals possessing a resistant ACCase and enhanced cytochrome P450 metabolism (ACCase/P450) were grown in the absence of mutual plant interaction to estimate plant growth traits. Both P450 and ACCase/P450 resistant phenotypes produced less above-ground biomass than the S phenotype during the vegetative stage. Reduced biomass production in the resistant phenotypes corresponded to a reduced relative growth rate and a lower net assimilation rate and rate of carbon fixation. There were no significant differences between the two resistant phenotypes, suggesting that costs of resistance are associated with P450 metabolism-based resistance. There were no differences in reproductive output among the three phenotypes, indicating that the cost of P450 resistance during vegetative growth is compensated during the production of reproductive structures. The P450-based herbicide metabolism is shown to be associated with physiological resistance costs, which may be manipulated by agronomic management to reduce the evolution of herbicide resistance.     

Molecular dynamic modeling of CYP51B in complex with azole inhibitors

Cytochrome P450 14α-sterol demethylase (CYP51), the key enzyme in sterol biosynthesis pathway, is an important target protein of cholesterol-lowering agents, antifungal drugs, and herbicides. CYP51B enzyme is one of the CYP51 family members. In the present study, we have chosen four representative inhibitors of CYP51B: diniconazole (Din), fluconazole (Flu), tebuconazole (Teb), and voriconazole (Vor), and launched to investigate the binding features of CYP51B-inhibitor and gating characteristics via molecular docking and molecular dynamics (MD) simulations. The results show that the ranking of binding affinities among these inhibitors is Din > Teb > Vor > Flu. Din shows the strongest binding affinity, whereas Flu shows the weakest binding affinity. More importantly, based on the calculated binding free energy results, we hypothesize that the nonpolar interactions are the most important contributors, and three key residues (Thr77, Ala258, and Lys454) play crucial role in protein-inhibitor binding. Besides, residue Phe180 may play a molecular switch role in the process of the CYP51B-Teb and CYP51B-Vor binding. Additionally, Tunnel analysis results show that the major tunnel of Din, Flu, and Teb is located between helix K, FG loop, and β4 hairpin (Tunnel II).The top ranked possible tunnel (Tunnel II) corresponds to Vor exits through helix K, F and helix J. This study further revealed the CYP51B relevant structural characteristics at the atomic level as well as provided a basis for rational design of new and more efficacious antifungal agents.     

Investigation of ligand selectivity in CYP3A7 by molecular dynamics simulations

Cytochrome P450 (CYP) 3A7 plays a crucial role in the biotransformation of the metabolized endogenous and exogenous steroids. To compare the metabolic capabilities of CYP3A7–ligands complexes, three endogenous ligands were selected, namely dehydroepiandrosterone (DHEA), estrone, and estradiol. In this study, a three-dimensional model of CYP3A7 was constructed by homology modeling using the crystal structure of CYP3A4 as the template and refined by molecular dynamics simulation (MD). The docking method was adopted, combined with MD simulation and the molecular mechanics generalized born surface area method, to probe the ligand selectivity of CYP3A7. These results demonstrate that DHEA has the highest binding affinity, and the results of the binding free energy were in accordance with the experimental conclusion that estrone is better than estradiol. Moreover, several key residues responsible for substrate specificity were identified on the enzyme. Arg372 may be the most important residue due to the low interaction energies and the existence of hydrogen bond with DHEA throughout simulation. In addition, a cluster of Phe residues provides a hydrophobic environment to stabilize ligands. This study provides insights into the structural features of CYP3A7, which could contribute to further understanding of related protein structures and dynamics.     

Herbicide-resistant tobacco plants expressing the fused enzyme between rat cytochrome P4501A1 (CYP1A1) and yeast NADPH-cytochrome P450 oxidoreductase.

Transgenic tobacco (Nicotiana tabacum cv Xanthi) plants expressing a genetically engineered fused enzyme between rat cytochrome P4501A1 (CYP1A1) and yeast NADPH-cytochrome P450 oxidoreductase were produced. The expression plasmid pGFC2 for the fused enzyme was constructed by insertion of the corresponding cDNA into the expression vector pNG01 under the control of the cauliflower mosaic virus 35S promoter and nopaline synthase gene terminator. The fused enzyme cDNA was integrated into tobacco genomes by Agrobacterium infection techniques. In transgenic tobacco plants, the fused enzyme protein was localized primarily in the microsomal fraction. The microsomal monooxygenase activities were approximately 10 times higher toward both 7-ethoxycoumarin and benzo[a]pyrene than in the control plant. The transgenic plants also showed resistance to the herbicide chlortoluron.     

Genetic control of a cytochrome P450 metabolism-based herbicide resistance mechanism in Lolium rigidum

The dynamics of herbicide resistance evolution in plants are influenced by many factors, especially the biochemical and genetic basis of resistance. Herbicide resistance can be endowed by enhanced rates of herbicide metabolism because of the activity of cytochrome P450 enzymes, although in weedy plants the genetic control of cytochrome P450-endowed herbicide resistance is poorly understood. In this study we have examined the genetic control of P450 metabolism-based herbicide resistance in a well-characterized Lolium rigidum biotype. The phenotypic resistance segregation in herbicide resistant and susceptible parents, F1, F2 and backcross (BC) families was analyzed as plant survival following treatment with the chemically unrelated herbicides diclofop-methyl or chlorsulfuron. Dominance and nuclear gene inheritance was observed in F1 families when treated at the recommended field doses of both herbicides. The segregation values of P450 herbicide resistance phenotypic traits observed in F2 and BC families was consistent with resistance endowed by two additive genes in most cases. In obligate out-crossing species such as L. rigidum, herbicide selection can easily result in accumulation of resistance genes within individuals.     

Modelling and predicting the binding mechanics of HIV P1053-0.5β antibody complex

Simulating antigen–antibody interactions is crucial in understanding the mechanics of antigen–antibody binding in medical science. In this study, molecular dynamics simulations are performed to analyse the dissociation of the P1053-0.5β antibody complex structure. The two-dimensional free energy profiles of the complex structure are extracted using the weighted histogram analysis method, and the binding pathway is then predicted using a modified form of the MaxFlux-PRM method. The simulation results suggest that 10 amino residues (i.e. Leu11, Val13, Asp34, Arg112, Thr101, Gly127, Val229, Ser231, Ile235 and Arg236) play a key role in relaxing the antibody structure, thereby facilitating the binding of the 0.5β antibody-P1053 peptide system.     

Molecular modeling of cytochrome P450 1A1: enzyme-substrate interactions and substrate binding affinities

Human cytochrome P450 1A1, which is present in lungs, plays an important role in the metabolic activation of chemical carcinogens, and in particular, is thought to be linked to lung cancer. The mechanism of carcinogenesis is related to the enzyme's ability to oxidize highly toxic compounds, such as polycyclic aromatic hydrocarbons (PAHs), to their carcinogenic derivatives. In order to better understand P450 1A1 function, a homology model of this enzyme has been constructed. The model has been based on the structure of P450 2C5, the first mammalian P450 to be crystallized. The coordinates of the model have been calculated using a consensus strategy, and the resulting structure has been evaluated with the ProStat and Profiles-3D programs. P450 1A1 substrates, such as benzo[a]pyrene, ethoxyresorufin and methoxyresorufin, were then docked into the active site of the model, and key amino acid residues able to interact with the substrate, have been identified. The analysis of enzyme-substrate interactions indicated that hydrophobic interactions are mainly responsible for binding of these substrates in the active site. Moreover, the non-bond enzyme-substrate interaction energy for ethoxyresorufin was lower than that for methoxyresorufin, which is consistent with higher activity of 1A1 towards the former substrate. Key residue Val-382 may play an important role in these interactions. Additionally, we performed binding free energy calculations for the three substrates. The obtained values were similar to those observed experimentally, which suggests that this approach might be useful for prediction of binding constants.     

Insight into the binding modes and inhibition mechanisms of adamantyl-based 1,3-disubstituted urea inhibitors in the active site of the human soluble epoxide hydrolase

Soluble epoxide hydrolase (sEH) is a promising new target for treating hypertension and inflammation. Considerable efforts have been devoted to develop novel inhibitors. In this study, the binding modes and interaction mechanisms of a series of adamantyl-based 1,3-disubstituted urea inhibitors were investigated by molecular docking, molecular dynamics simulations, binding free energy calculations, and binding energy decomposition analysis. Based on binding affinity, the most favorable binding mode was determined for each inhibitor. The calculation results indicate that the total binding free energy (ΔG_TOT, the sum of enthalpy ΔG_MM-GB/SA, and entropy ?TΔS) presents a good correlation with the experimental inhibitory activity (IC₅₀, r²?=?.99). The van der Waals energy contributes most to the total binding free energy (ΔG_TOT). A detailed discussion on the interactions between inhibitors and those residues located in the active pocket is made based on hydrogen bond and binding modes analysis. According to binding energy decomposition, the residues Asp333 and Trp334 contribute the most to binding free energy in all systems. Furthermore, Hip523 plays a major role in determining this class of inhibitor-binding orientations. Combined with the results of hydrogen bond analysis and binding free energy, we believe that the conserved hydrogen bonds play a role only in anchoring the inhibitors to the exact site for binding and the number of hydrogen bonds may not directly relate to the binding free energy. The results we obtained will provide valuable information for the design of high potency sEH inhibitors.     

Probing the binding mechanism of mercaptoguanine derivatives as inhibitors of HPPK by docking and molecular dynamics simulations

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) is a promising antimicrobial target involved in the folate biosynthesis pathway. Although, the results from crystallographic studies of HPPK have attracted a great interest in the design of novel HPPK inhibitors, the mechanism of action of HPPK due to inhibitor binding remains questionable. Recently, mercaptoguanine derivatives were reported to inhibit the pyrophosphoryl transfer mechanism of Staphylococcus aureus HPPK (SaHPPK). The present study is an attempt to understand the SaHPPK-inhibitors binding mechanism and to highlight the key residues that possibly involve in the complex formation. To decipher these questions, we used the state-of-the-art advanced insilico approach such as molecular docking, molecular dynamics (MD), molecular mechanics-generalized Born surface area approach. Domain cross correlation and principle component analysis were applied to the snapshots obtained from MD revealed that the compounds with high binding affinity stabilize the conformational dynamics of SaHPPK. The binding free energy estimation showed that the van der Waals and electrostatic interactions played a vital role for the binding mechanism. Additionally, the predicted binding free energy was in good agreement with the experimental values (R² = .78). Moreover, the free energy decomposition on per-residue confirms the key residues that significantly contribute to the complex formation. These results are expected to be useful for rational design of novel SaHPPK inhibitors.     

pH‐sensitive residues in the p19 RNA silencing suppressor protein from carnation Italian ringspot virus affect siRNA binding stability

Tombusviruses, such as Carnation Italian ringspot virus (CIRV), encode a protein homodimer called p19 that is capable of suppressing RNA silencing in their infected hosts by binding to and sequestering short‐interfering RNA (siRNA) away from the RNA silencing pathway. P19 binding stability has been shown to be sensitive to changes in pH but the specific amino acid residues involved have remained unclear. Using constant pH molecular dynamics simulations, we have identified key pH‐dependent residues that affect CIRV p19–siRNA binding stability at various pH ranges based on calculated changes in the free energy contribution from each titratable residue. At high pH, the deprotonation of Lys60, Lys67, Lys71, and Cys134 has the largest effect on the binding stability. Similarly, deprotonation of several acidic residues (Asp9, Glu12, Asp20, Glu35, and/or Glu41) at low pH results in a decrease in binding stability. At neutral pH, residues Glu17 and His132 provide a small increase in the binding stability and we find that the optimal pH range for siRNA binding is between 7.0 and 10.0. Overall, our findings further inform recent experiments and are in excellent agreement with data on the pH‐dependent binding profile.     

Engineering herbicide metabolism in tobacco and Arabidopsis with CYP76B1, a cytochrome P450 enzyme from Jerusalem artichoke

The Jerusalem artichoke (Helianthus tuberosus) xenobiotic inducible cytochrome P450, CYP76B1, catalyzes rapid oxidative dealkylation of various phenylurea herbicides to yield nonphytotoxic metabolites. We have found that increased herbicide metabolism and tolerance can be achieved by ectopic constitutive expression of CYP76B1 in tobacco (Nicotiana tabacum) and Arabidopsis. Transformation with CYP76B1 conferred on tobacco and Arabidopsis a 20-fold increase in tolerance to linuron, a compound detoxified by a single dealkylation, and a 10-fold increase in tolerance to isoproturon or chlortoluron, which need successive catalytic steps for detoxification. Two constructs for expression of translational fusions of CYP76B1 with P450 reductase were prepared to test if they would yield even greater herbicide tolerance. Plants expressing these constructs had lower herbicide tolerance than CYP76B1 alone, which is apparently a consequence of reduced stability of the fusion proteins. In all cases, increased herbicide tolerance results from more extensive metabolism, as demonstrated with exogenously fed phenylurea. Beside increased herbicide tolerance, expression of CYP76B1 has no other visible phenotype in the transgenic plants. Our data indicate that CYP76B1 can function as a selectable marker for plant transformation, allowing efficient selection in vitro and in soil-grown plants. Plants expressing CYP76B1 may also be a potential tool for phytoremediation of contaminated sites.     

Inducible cross-tolerance to herbicides in transgenic potato plants with the rat CYP1A1 gene

A gene of the enzyme involved in xenobiotic metabolism in mammalian liver was introduced into potato to confer inducible herbicide tolerance. A rat cytochrome P450 monooxygenase, CYP1A1 cDNA, was kept under the control of the tobacco PR1a promoter in order to apply the system of chemical inducible expression using the plant activator Benzothiadiazole (BTH). Transgenic plants were obtained based on the kanamycin resistance test and PCR analysis. Northern-blot analysis revealed the accumulation of mRNA corresponding to rat CYP1A1 in the transgenic plants treated with BTH (3.0 μmol/pot), whereas no accumulation of the corresponding mRNA occurred without BTH treatment. These transgenic plants also produced a protein corresponding to CYP1A1 in the leaves by BTH treatment. The transgenic plants with BTH application showed a much-higher tolerance to the phenylurea herbicides chlortoluron and methabenzthiazuron than non-transgenic plants. These findings indicated that the ability of metabolizing the two herbicides to less-toxic derivatives was displayed in the transgenic plants after BTH treatment. Transgenic plants harboring the CYP1A1 cDNA fused with the yeast P450 reductase (YR) gene under the control of PR1a were also produced. Although the plants showed a lower expression level of the fused gene than transgenic plants with CYP1A1 cDNA alone, they were tolerant to herbicides. These facts suggested that the CYP1A1 enzyme fused with YR showed a higher specific activity than CYP1A1 alone. This study demonstrated that the mammalian cDNA for the de-toxification enzyme of herbicides under the control of the PR1a promoter conferred chemical-inducible herbicide tolerance on potato. Received: 15 March 2001 / Accepted: 14 June 2001     

Identification by site-directed mutagenesis of two lysine residues in cholesterol side chain cleavage cytochrome P450 that are essential for adrenodoxin binding.

Utilizing site-directed mutagenesis and an Escherichia coli expression system for bovine cholesterol side chain cleavage cytochrome P450, lysine residues at 377 and 381 are found to play crucial roles in binding bovine adrenodoxin, required for transfer of electrons to mitochondrial P450s. These lysine residues are conserved among mitochondrial P450s and have been implicated previously by chemical modification studies as being important for adrenodoxin binding. In the present study, site-directed mutagenesis producing either neutral or positive amino acids at 377 or 381 has no effect on the structure of side chain cleavage cytochrome P450 as determined spectrally or on the enzymatic conversion of cholesterol to pregnenolone. However, the estimated Ks of adrenodoxin binding is increased approximately 150-600-fold depending on the particular mutation. Therefore these conserved positively charged residues in mitochondrial P450s are the key sites for adrenodoxin binding which is electrostatic in nature.