Molecular analysis of some camel cytochrome P450 enzymes reveals lower evolution and drug-binding properties

Camels bear unique genotypes and phenotypes for adaptation of their harsh environment. They have unique visual systems, sniffing, water metabolism, and heat-control mechanisms that are different from other creatures. The recent announcement for the complete sequence of camel genome will allow for the discovery of many secrets of camel life. In this context, the genetic bases of camel drug-metabolizing enzymes are still unknown. Furthermore, the genomic content of camel that rendered it highly susceptible to some drugs (as monensin and salinomycin) and became easily intoxicated needs to be investigated. The objectives of this work are the annotation of camel genome and retrieval of camel for cytochrome P450 (CYP) 1A1, 2C, and 3A enzymes. This is followed by comprehensive phylogenetic, evolution, molecular modeling, and docking studies. In comparison with the human enzymes, camel CYPs showed lower evolution rate, especially CYP1A1. Furthermore, the binding of monensin, salinomycin, alfa-naphthoflavone, felodepine, and ritonavir was weaker in camel enzymes. Interestingly, rerank score indicated instable binding of monensin and salinomycin with camel CYP1A1 as well as salinomycin with camel CYP2C. The results of this work suggest that camels are more susceptible to toxicity with compounds undergoing metabolic oxidation. This conclusion was based on lower evolution rate and lower binding potency of camels compared with the human enzymes.     

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.     

Histone acetyltransferase general control non‐repressed protein 5 (GCN5) affects the fatty acid composition of Arabidopsis thaliana seeds by acetylating fatty acid desaturase3 (FAD3)

Seed oils are important natural resources used in the processing and preparation of food. Histone modifications represent key epigenetic mechanisms that regulate gene expression, plant growth and development. However, histone modification events during fatty acid (FA) biosynthesis are not well understood. Here, we demonstrate that a mutation of the histone acetyltransferase GCN5 can decrease the ratio of α‐linolenic acid (ALA) to linoleic acid (LA) in seed oil. Using RNA‐Seq and ChIP assays, we identified FAD3, LACS2, LPP3 and PLAIIIβ as the targets of GCN5. Notably, the GCN5‐dependent H3K9/14 acetylation of FAD3 determined the expression levels of FAD3 in Arabidopsis thaliana seeds, and the ratio of ALA/LA in the gcn5 mutant was rescued to the wild‐type levels through the overexpression of FAD3. The results of this study indicated that GCN5 modulated FA biosynthesis by affecting the acetylation levels of FAD3. We provide evidence that histone acetylation is involved in FA biosynthesis in Arabidopsis seeds and might contribute to the optimization of the nutritional structure of edible oils through epigenetic engineering.     

cytochrome p450 superfamily in Arabidopsis thaliana: isolation of cDNAs,Differential Expression,and RFLP mapping of multiple cytochromes P450

We have isolated multiple cDNAs encoding cytochromes P450 (P450s) from Arabidopsis thaliana employing a PCR strategy. Degenerate oligonucleotide primers were designed from amino acid sequences conserved between two plant P450s, CYP71A1 and CYP73A2, including the heme-binding site and the proline-rich motif found in the N-terminal region, and 11 putative P450 fragments were amplified from first-strand cDNA from 7-day-old Arabidopsis as a template. With these PCR fragments as hybridization probes, 13 full-length and 3 partial cDNAs encoding different P450s have been isolated from an Arabidopsis cDNA library. These P450s have been assigned to either one of the established subfamilies: CYP71B, CYP73A, and CYP83A; or novel subfamilies: CYP76C, CYP83B, and CYP91A. The primary protein structures predicted from the cDNA sequences revealed that the regions around both the heme-binding site and the proline-rich motif were highly conserved among all these P450s. The N-terminal structures of the predicted P450 proteins suggested that these Arabidopsis P450s were located at the endoplasmic reticulum membrane. The loci of four P450 genes were determined by RFLP mapping. One of the clones, CYP71B2, was located at a position very close to the ga4 and gai mutations. RNA blot analysis showed expression patterns unique to each of the P450s in terms of tissue specificity and responsiveness to wounding and light/dark cycle, implicating involvement of these P450s in diverse metabolic processes.     

Probing ligand binding modes of human cytochrome P450 2J2 by homology modeling, molecular dynamics simulation, and flexible molecular docking

Cytochrome P450 (P450) 2J2 catalyzes epoxidation of arachidonic acid to eicosatrienoic acids, which are related to a variety of diseases such as coronary artery disease, hypertension, and carcinogenesis. Recent experimental data also suggest that P450 2J2 could be a novel biomarker and a potential target for cancer therapy. However, the active site topology and substrate specificity of this enzyme remain unclear. In this study, a three-dimensional model of human P450 2J2 was first constructed on the basis of the crystal structure of human P450 2C9 in complex with a substrate using homology modeling method, and refined by molecular dynamics simulation. Flexible docking approaches were then employed to dock four ligands into the active site of P450 2J2 in order to probe the ligand-binding modes. By analyzing the results, active site architecture and certain key residues responsible for substrate specificity were identified on the enzyme, which might be very helpful for understanding the enzyme's biological role and providing insights for designing novel inhibitors of P450 2J2.     

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.     

Molecular Modeling of Cytochrome P450 2B1: Mode of Membrane Insertion and Substrate Specificity

A molecular model of a mammalian membrane-bound cytochrome P450, rat P450 2B1, was constructed in order to elucidate its mode of attachment to the endoplasmic reticulum and the structural basis of substrate specificity. The model was primarily derived from the structure of P450BM-3, which as a class II P450 is the most functionally similar P450 of known structure. However, model development was also guided by the conserved core regions of P450cam and P450terp. To optimally align the P450 2B1 and P450BM-3 sequences, multiple alignment was performed using sequences of five P450s in the II family, followed by minor adjustments on the basis of secondary structure predictions. The resulting P450 2B1 homology model structure was refined by molecular dynamics heating, equilibration, simulation, and energy minimization. The model suggests that the F–G loop serves as both a hydrophobic membrane anchor and entrance channel for hydrophobic substrates from the membrane to the P450 active site. To assess the mode of substrate binding, benzphetamine, testosterone, and benzo[a]pyrene were docked into the active site. The hydrophobic substrate-binding pocket is consistent with the preferences of this P450 toward hydrophobic substrates, while the presence of an acidic Glu-105 in this pocket is consistent with the preference of this P450 for the cationic substrate benzphetamine. This model is thus consistent with several known experimental properties of this P450, such as membrane attachment and substrate selectivity.     

CYP707A3, a major ABA 8'-hydroxylase involved in dehydration and rehydration response in Arabidopsis thaliana

Abscisic acid (ABA) catabolism is one of the determinants of endogenous ABA levels affecting numerous aspects of plant growth and abiotic stress responses. The major ABA catabolic pathway is triggered by ABA 8'-hydroxylation catalysed by the cytochrome P450 CYP707A family. Among four members of Arabidopsis CYP707As, the expression of CYP707A3 was most highly induced in response to both dehydration and subsequent rehydration. A T-DNA insertional cyp707a3-1 mutant contained higher ABA levels in turgid plants, which showed a reduced transpiration rate and hypersensitivity to exogenous ABA during early seedling growth. On dehydration, the cyp707a3-1 mutant accumulated a higher amount of stress-induced ABA than the wild type, an event that occurred relatively later and was coincident with slow drought induction of CYP707A3. The cyp707a3 mutant plants exhibited both exaggerated ABA-inducible gene expression and enhanced drought tolerance. Conversely, constitutive expression of CYP707A3 relieved growth retardation by ABA, increased transpiration, and a reduction of endogenous ABA in both turgid and dehydrated plants. Taken together, our results indicate that CYP707A3 plays an important role in determining threshold levels of ABA during dehydration and after rehydration.     

Computational screening of fatty acid synthase inhibitors against thioesterase domain

Thioesterase (TE) domain of fatty acid synthase (FAS) is an attractive therapeutic target for design and development of anticancer drugs. In this present work, we search for the potential FAS inhibitors of TE domain from the ZINC database based on similarity search using three natural compounds as templates, including flavonoids, terpenoids, and phenylpropanoids. Molecular docking was used to predict the interaction energy of each screened ligand compared to the reference compound, which is methyl γ-linolenylfluorophosphonate (MGLFP). Based on this computational technique, rosmarinic acid and its eight analogs were observed as a new series of potential FAS inhibitors, which showed a stronger binding affinity than MGLFP. Afterward, nine docked complexes were studied by molecular dynamics simulations for investigating protein–ligand interactions and binding free energies using MM-PB(GB)SA, MM-3DRISM-KH, and QM/MM-GBSA methods. The binding free energy calculation indicated that the ZINC85948835 (R34) displayed the strongest binding efficiency against the TE domain of FAS. There are eight residues (S2308, I2250, E2251, Y2347, Y2351, F2370, L2427, and E2431) mainly contributed for the R34 binding. Moreover, R34 could directly form hydrogen bonds with S2308, which is one of the catalytic triad of TE domain. Therefore, our finding suggested that R34 could be a potential candidate as a novel FAS-TE inhibitor for further drug design.     

Long-chain acyl-CoA synthetase isoforms differ in preferences for eicosanoid species and long-chain fatty acids

Because the signaling eicosanoids, epoxyeicosatrienoic acids (EETs) and HETEs, are esterified to membrane phospholipids, we asked which long-chain acyl-CoA synthetase (ACSL) isoforms would activate these molecules and whether the apparent FA substrate preferences of each ACSL isoform might differ depending on whether it was assayed in mammalian cell membranes or as a purified bacterial recombinant protein. We found that all five ACSL isoforms were able to use EETs and HETEs as substrates and showed by LC-MS/MS that ACSLs produce EET-CoAs. We found differences in substrate preference between ACS assays performed in COS7 cell membranes and recombinant purified proteins. Similarly, preferences and Michaelis-Menten kinetics for long-chain FAs were distinctive. Substrate preferences identified for the purified ACSLs did not correspond to those observed in ACSL-deficient mouse models. Taken together, these data support the concept that each ACSL isoform exhibits a distinct substrate preference, but apparent substrate specificities depend upon multiple factors including membrane character, coactivators, inhibitors, protein interactions, and posttranslational modification.     

Molecular analysis of the myosin gene family in Arabidopsis thaliana

Myosin is believed to act as the molecular motor for many actin-based motility processes in eukaryotes. It is becoming apparent that a single species may possess multiple myosin isoforms, and at least seven distinct classes of myosin have been identified from studies of animals, fungi, and protozoans. The complexity of the myosin heavy-chain gene family in higher plants was investigated by isolating and characterizing myosin genomic and cDNA clones from Arabidopsis thaliana. Six myosin-like genes were identified from three polymerase chain reaction (PCR) products (PCR1, PCR11, PCR43) and three cDNA clones (ATM2, MYA2, MYA3). Sequence comparisons of the deduced head domains suggest that these myosins are members of two major classes. Analysis of the overall structure of the ATM2 and MYA2 myosins shows that they are similar to the previously-identified ATM1 and MYA1 myosins, respectively. The MYA3 appears to possess a novel tail domain, with five IQ repeats, a six-member imperfect repeat, and a segment of unique sequence. Northern blot analyses indicate that some of the Arabidopsis myosin genes are preferentially expressed in different plant organs. Combined with previous studies, these results show that the Arabidopsis genome contains at least eight myosin-like genes representing two distinct classes.     

On the specificity of lipid hydroperoxide fragmentation by fatty acid hydroperoxide lyase from Arabidopsis thaliana

Fatty acid hydroperoxide lyase (HPL) is a membrane associated P450 enzyme that cleaves fatty acid hydroperoxides into aldehydes and omega-oxo fatty acids. One of the major products of this reaction is (3Z)-hexenal. It is a constituent of many fresh smelling fruit aromas. For its biotechnological production and because of the lack of structural data on the HPL enzyme family, we investigated the mechanistic reasons for the substrate specificity of HPL by using various structural analogues of HPL substrates. To approach this 13-HPL from Arabidopsis thaliana was cloned and expressed in E. coli utilising a His-Tag expression vector. The fusion protein was purified by affinity chromatography from the E. coli membrane fractions and its pH optimum was detected to be pH 7.2. Then, HPL activity against the respective (9S)- and (13S)-hydroperoxides derived either from linoleic, alpha-linolenic or gamma-linolenic acid, respectively, as well as that against the corresponding methyl esters was analysed. Highest enzyme activity was observed with the (13S)-hydroperoxide of alpha-linolenic acid (13alpha-HPOT) followed by that with its methyl ester. Most interestingly, when the hydroperoxy isomers of gamma-linolenic acid were tested as substrates, 9gamma-HPOT and not 13gamma-HPOT was found to be a better substrate of the enzyme. Taken together from these studies on the substrate specificity it is concluded that At13HPL may not recognise the absolute position of the hydroperoxy group within the substrate, but shows highest activities against substrates with a (1Z4S,5E,7Z)-4-hydroperoxy-1,5,7-triene motif. Thus, At13HPL may not only be used for the production of C6-derived volatiles, but depending on the substrate may be further used for the production of Cg-derived volatiles as well.     

Indole-3-butyric acid in Arabidopsis thaliana

Indole-3-butyric acid (IBA) was identified by HPLC and GC-MS as an endogenous compound in plantlets of the crucifer Arabidopsis thaliana (L.) Heynh. A. thaliana was cultivated under sterile conditions as shaking culture in different liquid media with and without supply of hormones. Free and total IBA and indole-3-acetic acid (IAA) were determined at different stages of development during the culture period as well as in culture media of different initial pH values. The results showed that IAA was present in higher concentrations than IBA, but both hormones seemed to show the same behaviour under the different experimental conditions. Differences were found in the mode of conjugation of the two hormones. While IAA was mostly conjugated via amide bonds, the main IBA conjugates were ester bound. The ethylene concentration derived from the seedlings, when they were grown in flasks of different size, seemed not to influence the auxin content in the same cultures.     

Arabidopsis CYP85A2 Catalyzes Lactonization Reactions in the Biosynthesis of 2-Deoxy-7-oxalactone Brassinosteroids

Brassinolide (BL), a plant 7-oxalactone-type steroid hormone, is one of the active brassinosteroids (BRs) that regulates plant growth and development. BL is biosynthesized from castasterone by the cytochrome P450 monooxygenase, CYP85A2. We showed that a Pichia pastoris transformant that synchronously expresses Arabidopsis P450 reductase gene ATR1 and P450 gene CYP85A2 converts teasterone and typhasterol to 7-oxateasterone and 7-oxatyphasterol, respectively. Thus, CYP85A2 catalyzes the lactonization reactions of not only castasterone but also teasterone and typhasterol. The two 2-deoxy-7-oxalactone-type BRs were identified in Arabidopsis plants. Although the reversible conversion between 7-oxateasterone and 7-oxatyphasterol was observed in vivo, no conversion of 7-oxatyphasterol to BL was observed. The biological activity of 7-oxatyphasterol toward Arabidopsis hypocotyl elongation was nearly the same as that of castasterone. These results suggest that a new BR biosynthetic pathway, a BR lactonization pathway, functions in Arabidopsis and plays an important role in regulating the concentration of active BRs, even though the metabolism of 7-oxatyphasterol to BL is still unknown.     

π–π Stacking mediated drug–drug interactions in human CYP2E1

Because of having many low molecular mass substrates, CYP2E1 is of particular interests to the pharmaceutical industry. Many evidences showed that this enzyme can adopt multiple substrates to significantly reduce the oxidation rate of the substrates. The detailed mechanism for this observation is still unclear. In the current study, we employed GPU‐accelerated molecular dynamics simulations to study the multiple‐binding mode of human CYP2E1, with an aim of offering a mechanistic explanation for the unexplained multiple‐substrate binding. Our results showed that Thr303 and Phe478 were key factors for the substrate recognition and multiple‐substrate binding. The former can form a significant hydrogen bond to recognize and position the substrate in the productive binding orientation in the active site. The latter acted as a mediator for the substrate communications via π–π stacking interactions. In the multiple‐binding mode, the aforementioned π–π stacking interactions formed by the aromatic rings of both substrates and Phe478 drove the first substrate far away from the catalytic center, orienting in an additional binding position and going against the substrate metabolism. All these findings could give atomic insights into the detailed mechanism for the multiple‐substrate binding in human CYP2E1, providing useful information for the drug metabolism mechanism and personalized use of clinical drugs. Proteins 2013; © 2012 Wiley Periodicals, Inc.     

Comparative modeling and molecular docking of orphan human CYP4V2 protein with fatty acid substrates: Insights into substrate specificity

Cytochromes P450 (CYPs) are a super family of heme-containing enzymes well-known for their monooxgenase reaction. There are 57 CYP isoenzymes found in human which exhibit specific physiological functions. Thirteen members of this super family are classified as "orphan" CYP because of their unknown enzymatic functions. CYP4V2 is found to be a potential drug target for Bietti crystalline corneoretinal dystrophy (BCD). However, three-dimensional structure, the active site topology and substrate binding modes of CYP4V2 remain unclear. In this study, the three-dimensional model of CYP4V2 was constructed using the homology modeling method. Four possible fatty acid substrates namely, caprylic, lauric, myrisitc and palmitic acids were optimized and evaluated for drug likeness using Lipinski's rule of five. Further, these substrates were docked into active sites of CYP4V2 and several key residues responsible for substrate binding were identified. These findings will be helpful for the structure-based drug design and detailed characterization of the biological roles of CYP4V2.