Mechanistic and structural insights into the regioselectivity of an acyl-CoA fatty acid desaturase via directed molecular evolution

Membrane-bound fatty acid desaturases and related enzymes play a pivotal role in the biosynthesis of unsaturated and various unusual fatty acids. Structural insights into the remarkable catalytic diversity and wide range of substrate specificities of this class of enzymes remain limited due to the lack of a crystal structure. To investigate the structural basis of the double bond positioning (regioselectivity) of the desaturation reaction in more detail, we relied on a combination of directed evolution in vitro and a powerful yeast complementation assay to screen for Δx regioselectivity. After two selection rounds, variants of the bifunctional Δ12/Δ9-desaturase from the house cricket (Acheta domesticus) exhibited increased Δ9-desaturation activity on shorter chain fatty acids. This change in specificity was the result of as few as three mutations, some of them near the putative active site. Subsequent analysis of individual substitutions revealed an important role of residue Phe-52 in facilitating Δ9-desaturation of shorter chain acyl substrates and allowed for the redesign of the cricket Δ12/Δ9-desaturase into a 16:0-specific Δ9-desaturase. Our results demonstrate that a minimal number of mutations can have a profound impact on the regioselectivity of acyl-CoA fatty acid desaturases and include the first biochemical data supporting the acyl-CoA acyl carrier specificity of a desaturase able to carry out Δ12-desaturation.     

Metabolic engineering of Saccharomyces cerevisiae for production of Eicosapentaenoic Acid, using a novel {Delta}5-Desaturase from Paramecium tetraurelia

Very-long-chain polyunsaturated fatty acids, such as arachidonic acid (ARA), eicosapentaenoic acid (EPA), and docosahexaenoic acid (DHA), have well-documented importance in human health and nutrition. Sustainable production in robust host organisms that do not synthesize them naturally requires the coordinated expression of several heterologous desaturases and elongases. In the present study we show production of EPA in Saccharomyces cerevisiae using glucose as the sole carbon source through expression of five heterologous fatty acid desaturases and an elongase. Novel Δ5-desaturases from the ciliate protozoan Paramecium tetraurelia and from the microalgae Ostreococcus tauri and Ostreococcus lucimarinus were identified via a BLAST search, and their substrate preferences and desaturation efficiencies were assayed in a yeast strain producing the ω6 and ω3 fatty acid substrates for Δ5-desaturation. The Δ5-desaturase from P. tetraurelia was up-to-2-fold more efficient than the microalgal desaturases and was also more efficient than Δ5-desaturases from Mortierella alpina and Leishmania major. In vivo investigation of acyl carrier substrate specificities showed that the Δ5-desaturases from P. tetraurelia, O. lucimarinus, O. tauri, and M. alpina are promiscuous toward the acyl carrier substrate but prefer phospholipid-bound substrates. In contrast, the Δ5-desaturase from L. major showed no activity on phospholipid-bound substrate and thus appears to be an exclusively acyl coenzyme A-dependent desaturase.     

Histidine-41 of the Cytochrome b5 Domain of the Borage Δ6 Fatty Acid Desaturase Is Essential for Enzyme Activity

Unlike most other plant microsomal desaturases, the Delta6-fatty acid desaturase from borage (Borago officinalis) contains an N-terminal extension that shows homology to the small hemoprotein cytochrome (Cyt) b5. To determine if this domain serves as a functional electron donor for the Delta6-fatty acid desaturase, mutagenesis and functional analysis by expression in transgenic Arabidopsis was carried out. Although expression of the wild-type borage Delta6-fatty acid desaturase resulted in the synthesis and accumulation of Delta6-unsaturated fatty acids, this was not observed in plants transformed with N-terminally deleted forms of the desaturase. Site-directed mutagenesis was used to disrupt one of the axial heme-binding residues (histidine-41) of the Cyt b5 domain; expression of this mutant form of the Delta6-desaturase in transgenic plants failed to produce Delta6-unsaturated fatty acids. These data indicate that the Cyt b5 domain of the borage Delta6-fatty acid desaturase is essential for enzymatic activity.     

Distinct roles of endoplasmic reticulum cytochrome b5 and fused cytochrome b5-like domain for rat Delta6-desaturase activity

The Delta6-desaturase catalyzes key steps in long-chain polyunsaturated fatty acid biosynthesis. Although the gene coding for this enzyme has been isolated in diverse animal species, the protein structure remains poorly characterized. In this work, rat Delta6-desaturase expressed in COS-7 cells was shown to localize in the endoplasmic reticulum. As the enzyme contains an N-terminal cytochrome b5-like domain, we investigated by site-directed mutagenesis the role of this domain in the enzyme activity. The typical HPGG motif of the cytochrome b5-like domain, and particularly histidine in this motif, is required for the activity of the enzyme, whatever the substrate. Neither endogenous COS-7 cytochrome b5 nor coexpressed rat endoplasmic reticulum cytochrome b5 could rescue the activity of mutated forms of Delta6-desaturase. Moreover, when rat endoplasmic reticulum cytochrome b5 was coexpressed with wild-type desaturase, both proteins interacted and Delta6-desaturase activity was significantly increased. The identified interaction between these proteins is not dependent on the desaturase HPGG motif. These data suggest distinct and essential roles for both the desaturase cytochrome b5-like domain and free endoplasmic reticulum cytochrome b5 for Delta6-desaturase activity.     

The role of Δ6-desaturase acyl-carrier specificity in the efficient synthesis of long-chain polyunsaturated fatty acids in transgenic plants

The role of acyl‐CoA‐dependent Δ6‐desaturation in the heterologous synthesis of omega‐3 long‐chain polyunsaturated fatty acids was systematically evaluated in transgenic yeast and Arabidopsis thaliana. The acyl‐CoA Δ6‐desaturase from the picoalga Ostreococcus tauri and orthologous activities from mouse (Mus musculus) and salmon (Salmo salar) were shown to generate substantial levels of Δ6‐desaturated acyl‐CoAs, in contrast to the phospholipid‐dependent Δ6‐desaturases from higher plants that failed to modify this metabolic pool. Transgenic plants expressing the acyl‐CoA Δ6‐desaturases from either O. tauri or salmon, in conjunction with the two additional activities required for the synthesis of C20 polyunsaturated fatty acids, contained higher levels of eicosapentaenoic acid compared with plants expressing the borage phospholipid‐dependent Δ6‐desaturase. The use of acyl‐CoA‐dependent Δ6‐desaturases almost completely abolished the accumulation of unwanted biosynthetic intermediates such as γ‐linolenic acid in total seed lipids. Expression of acyl‐CoA Δ6‐desaturases resulted in increased distribution of long‐chain polyunsaturated fatty acids in the polar lipids of transgenic plants, reflecting the larger substrate pool available for acylation by enzymes of the Kennedy pathway. Expression of the O. tauriΔ6‐desaturase in transgenic Camelina sativa plants also resulted in the accumulation of high levels of Δ6‐desaturated fatty acids. This study provides evidence for the efficacy of using acyl‐CoA‐dependent Δ6‐desaturases in the efficient metabolic engineering of transgenic plants with high value traits such as the synthesis of omega‐3 LC‐PUFAs.     

In Silico Structural Studies and Molecular Docking Analysis of Delta6-desaturase in HUFA Biosynthetic Pathway

Fish are an important source of highly unsaturated fatty acids (HUFA) such as eicosapentaenoic acid EPA (20:5 n-3) and docosahexaenoic acid DHA (22:6 n-3) and play a significant role in human nutrition. The fatty acyl delta6-desaturase (Δ6 desaturase) is a rate-limiting enzyme in the biosynthetic pathway of highly unsaturated fatty acids (HUFA) that converts polyunsaturated fatty acids (PUFA) such as linoleic (18:2n-6) and α-linolenic (18:3n-3) acids into HUFA. In this study, fatty acyl Δ6 desaturase was identified from pangasius (Pangasianodon hypophthalmus) and further analyzed for sequenced-based characterization and 3D structural conformation. Sequenced-based analysis revealed some important secondary information such as physicochemical property. e.g., isoelectric point, extinction coefficient, aliphatic index, and grand average hydropathy, among others, and also post-translational modification sites were identified. An evolutionary-conserved stretch of amino acid residue and a functionally significant conserved structural ancestor, N-terminal cytochrome b5 and membrane FADS-like superfamily, were identified. Protein association analysis showed a high confidence score with acyl-CoA synthetase, elovl5, elovl2, and phospholipase A2. Herein, we report, for the first time, a 3D native structure of Δ6 desaturase protein by homology modeling approach; molecular docking analysis was performed with linoleic (18:2n-6) and α-linolenic (18:3n-3) acids, which are the two key substrates in the HUFA biosynthetic pathway. This work provides insight into the structural and functional characterization of Δ6 desaturase, which is involved in HUFA biosynthesis as a rate-limiting enzyme.     

The alternative pathway C20 Delta8-desaturase from the non-photosynthetic organism Acanthamoeba castellanii is an atypical cytochrome b5-fusion desaturase

A cDNA encoding a C20 Delta8-desaturase was isolated from the free-living soil amoeba, Acanthamoeba castellanii and functionally characterised by heterologous expression. The open reading frame of the A. castellanii C20 Delta8-desaturase showed similarity to other microsomal front-end desaturases, but the N-terminal domain contained a variant form of the conserved heme-binding motif in which H-P-G-G is replaced by H-P-A-G. Co-expression of the A. castellani Delta8-desaturase with the Isochrysis galbana Delta9-elongase in transgenic Arabidopsis plants confirmed the activity observed in yeast and its role in the alternative pathway for C20 polyunsaturated fatty acid synthesis. Acyl-CoA profiles of these transgenic plants revealed an unexpected accumulation of C20 fatty acids in the acyl-CoA pool. This is the first report of an alternative pathway C20 Delta8-desaturase from a non-photosynthetic organism, and also the first report of a front-end desaturase lacking the canonical cytochrome b5 domain.     

深黄被孢霉△~6-脂肪酸脱氢酶基因的克隆及序列分析

γ－亚麻酸（ＧＬＡ,Ｃ１８：３△６,９,１２）是由△６－脂肪酸脱氢酶以亚油酸（ＬＡ,Ｃ１８：２△９,１２）为底物,在Ｃ６位脱氢形成的。由于在人体中,γ－亚麻酸是花生四烯酸、前列腺素类和白三烯类等生理活性物质的前体物,而深黄被孢霉是目前用于微生物发酵生产γ－亚麻酸的主要菌株。本文根据脂肪酸脱氢酶的保守区设计引物,利用反转录聚合酶链式反应从丝状真菌深黄被孢霉中克隆了编码△６－脂肪酸脱氢酶的ｃＤＮＡ,全长为１３７４个核苷酸,编码４５７ 个氨基酸,但与其他位点的脂肪酸脱氢酶不同的是, △６－脂肪酸脱氢酶在其序列的 Ｎ 端特有细胞色素 ｂ５（Ｃｙｔｂ５）区。这是国际上对深黄被孢霉△６－脂肪酸脱氢酶基因的首次报道。     

The Evolution of Fatty Acid Desaturases and Cytochrome b5 in Eukaryotes

Desaturases that introduce double bonds into the fatty acids are involved in the adaptation of membrane fluidity to changes in the environment. Besides, polyunsaturated fatty acids (PUFAs) are increasingly recognized as important pharmaceutical and nutraceutical compounds. To successfully engineer organisms with increased stress tolerance or the ability to synthesize valuable PUFAs, detailed knowledge about the complexity of the desaturase family as well as understanding of the coevolution of desaturases and their cytochrome b5 electron donors is needed. We have constructed phylogenies of several hundred desaturase sequences from animals, plants, fungi and bacteria and of the cytochrome b5 domains that are fused to some of these enzymes. The analysis demonstrates the existence of three major desaturase acyl-CoA groups that share few similarities. Our results indicate that the fusion of Δ⁶-desaturase-like enzymes with their cytochrome b5 electron donor was a single event that took place in the common ancestor of all eukaryotes. We also propose the Δ⁶-desaturase-like enzymes as the most probable donor of the cytochrome b5 domain found in fungal Δ⁹-desaturases and argue that the recombination most likely happened soon after the separation of the animal and fungal ancestors. These findings answer some of the previously unresolved questions and contribute to the quickly expanding field of research on desaturases.     

Identification of a fatty acid Delta11-desaturase from the microalga Thalassiosira pseudonana

A set of genomic DNA sequences putatively encoding front-end desaturases were identified by in silico analysis of the draft genome of the marine microalga Thalassiosira pseudonana. Among these candidate genes, an open reading frame named TpdesN was found to be full-length, intronless, and constitutively expressed during cell cultivation. The predicted amino acid sequence of the corresponding protein, TpDESN, exhibited typical features of desaturases involved in the production of polyunsaturated fatty acids (PUFAs) in algae, i.e. a cytochrome b5-like domain at the N-terminus and three conserved histidine-rich motifs in the desaturase domain. Expression of TpDESN in Saccharomyces cerevisiae revealed that this enzyme was not involved in PUFA synthesis, but specifically desaturated palmitic acid 16:0 to 16:1Delta11. To our knowledge, until this report, Delta11-desaturase activity had only been detected in insect cells.     

A front-end desaturase from Chlamydomonas reinhardtii produces pinolenic and coniferonic acids by omega13 desaturation in methylotrophic yeast and tobacco

Pinolenic acid (PA; 18:3Delta(5,9,12)) and coniferonic acid (CA; 18:4Delta(5,9,12,15)) are Delta(5)-unsaturated bis-methylene-interrupted fatty acids (Delta(5)-UBIFAs) commonly found in pine seed oil. They are assumed to be synthesized from linoleic acid (LA; 18:2Delta(9,12)) and alpha-linolenic acid (ALA; 18:3Delta(9,12,15)), respectively, by Delta(5)-desaturation. A unicellular green microalga Chlamydomonas reinhardtii also accumulates PA and CA in a betain lipid. The expressed sequence tag (EST) resource of C. reinhardtii led to the isolation of a cDNA clone that encoded a putative fatty acid desaturase named as CrDES containing a cytochrome b5 domain at the N-terminus. When the coding sequence was expressed heterologously in the methylotrophic yeast Pichia pastoris, PA and CA were newly detected and comparable amounts of LA and ALA were reduced, demonstrating that CrDES has Delta(5)-desaturase activity for both LA and ALA. CrDES expressed in the yeast showed Delta(5)-desaturase activity on 18:1Delta(9) but not 18:1Delta(11). Unexpectedly, CrDES also showed Delta(7)-desaturase activity on 20:2Delta(11,14) and 20:3Delta(11,14,17) to produce 20:3Delta(7,11,14) and 20:4Delta(7,11,14,17), respectively. Since both the Delta(5) bond in C18 and the Delta(7) bond in C20 fatty acids are 'omega13' double bonds, these results indicate that CrDES has omega13 desaturase activity for omega9 unsaturated C18/C20 fatty acids, in contrast to the previously reported front-end desaturases. In order to evaluate the activity of CrDES in higher plants, transgenic tobacco plants expressing CrDES were created. PA and CA accumulated in the leaves of transgenic plants. The highest combined yield of PA and CA was 44.7% of total fatty acids, suggesting that PA and CA can be produced in higher plants on a large scale.     

Delta-8 desaturation activity varies among fatty acyl desaturases of teleost fish: high activity in delta-6 desaturases of marine species

The benefits of dietary fish and fish oil are derived from n-3 long-chain polyunsaturated fatty acids (LC-PUFA) that have beneficial effects in a range of human diseases and pathologies such as cardiovascular and other inflammatory disorders, neural development and neurological pathologies. The precursor of n-3 LC-PUFA, 18:3n-3 does not have the same beneficial effects prompting interest in the pathways of endogenous synthesis of LC-PUFA in vertebrates. The LC-PUFA biosynthesis pathway classically involves Δ6 and Δ5 fatty acyl desaturases (Fad), but it was recently shown that Δ6 Fad in mammals also displayed Δ8 activity demonstrating a possible alternative "Δ8-pathway" for the synthesis of LC-PUFA. Our primary hypothesis was that Δ8 desaturase activity would be a common feature of vertebrate Δ6 Fads, and so the aim of the present study was to determine the ability of teleostei Fads for Δ8 desaturation activity. To this end, cDNAs for Fads from a range of freshwater, diadromous and marine teleost fish species were assayed for Δ8 activity in the heterologous yeast expression system. In summary, the present study has demonstrated that Δ8 desaturation activity was also a characteristic of fish orthologs, although the activity varied notably between freshwater/diadromous and marine fish species, with the latter possessing Fads2-like proteins with Δ8 activity far higher than mammalian FADS2. The data showed that, generally, the fish Fad are technically υ-3 desaturases, with new double bonds introduced 3C beyond a pre-existing double bond. However, the ability of zebrafish and rabbitfish Fads, previously characterised as Δ6/Δ5 bifunctional desaturases, to introduce non-methylene interrupted double bonds in 20:3n-3 and 20:2n-6 suggested that a novel combination of regioselectivity modes operates within these enzymes.     

A second functional delta5 fatty acid desaturase in the cellular slime mould Dictyostelium discoideum.

A cDNA with homology to fatty acid desaturases was selected by searching the cDNA data bank of Dictyostelium discoideum (http://www. csm.biol.tsukuba.ac.jp/cDNAproject.html) with conserved histidine box motifs. Using this sequence, genomic DNA encoding the Delta5 desaturase was amplified from the genomic DNA of D. discoideum, and its desaturase activity was confirmed by the overexpression mutation in D. discoideum and the gain-of-function mutation in yeast. The cloned cDNA is 1565 nucleotides in length, and the deduced amino-acid sequence comprised 467 amino-acid residues containing an N-terminal cytochrome b5 domain that shared 43% identity with cytochrome b5 of Oryza sativa. The whole sequence was 42% identical to the Delta5 desaturase of Mortierella alpina. This desaturase is a novel member of the cytochrome b5-containing Delta5 fatty acid desaturase. As we have already reported one other Delta5 desaturase in Dictyostelium, this organism is the first to be confirmed as having two functional Delta5 fatty acid desaturase genes. The substrate specificities of the two functional Delta5 desaturases of D. discoideum were also examined.     

球等鞭金藻中Δ5去饱和酶基因的克隆与功能鉴定

球等鞭金藻(Isochrysis galbana)是一类单细胞海洋微藻,富含二十二碳六烯酸(DHA,22:6Δ4,7,10,13,16,19)。我们利用RACE的方法从球等鞭金藻cDNA文库中同源克隆到一个大小为1329 bp的cDNA片段,编码442个氨基酸的多肽,分子量约49.9 kD。生物信息学分析表明,其编码产物N端具有细胞色素b5结构域,以及与电子传递有关的三个富含组氨酸的结构域,与Pavlova salinaΔ5去饱和酶同源性最高,达56%,故将该基因命名为IgD5。酿酒酵母功能鉴定实验表明,其编码的蛋白质具有Δ5去饱和酶活性,能够将二高-γ-亚麻酸(DGLA,20:3Δ8,11,14)转化成花生四烯酸(AA,20:4Δ5,8,11,14),转化效率平均为34.6%,最高可达40.3%。     

Fatty acid desaturases from the microalga Thalassiosira pseudonana

Analysis of a draft nuclear genome sequence of the diatom Thalassiosira pseudonana revealed the presence of 11 open reading frames showing significant similarity to functionally characterized fatty acid front-end desaturases. The corresponding genes occupy discrete chromosomal locations as determined by comparison with the recently published genome sequence. Phylogenetic analysis showed that two of the T. pseudonana desaturase (Tpdes) sequences grouped with proteobacterial desaturases that lack a fused cytochrome b5 domain. Among the nine remaining gene sequences, temporal expression analysis revealed that seven were expressed in T. pseudonana cells. One of these, TpdesN, was previously characterized as encoding a Delta11-desaturase active on palmitic acid. From the six remaining putative desaturase genes, we report here that three, TpdesI, TpdesO and TpdesK, respectively encode Delta6-, Delta5- and Delta4-desaturases involved in production of the health beneficial polyunsaturated fatty acid DHA (docosahexaenoic acid). Furthermore, we show that one of the remaining genes, TpdesB, encodes a Delta8-sphingolipid desaturase with strong preference for dihydroxylated substrates.     

Newly identified essential amino acid residues affecting Δ8-sphingolipid desaturase activity revealed by site-directed mutagenesis

In order to identify amino acid residues crucial for the enzymatic activity of Δ(8)-sphingolipid desaturases, a sequence comparison was performed among Δ(8)-sphingolipid desaturases and Δ(6)-fatty acid desaturases from various plants. In addition to the known conserved cytb(5) (cytochrome b(5)) HPGG motif and three conserved histidine boxes, they share additional 15 completely conserved residues. A series of site-directed mutants were generated using our previously isolated Δ(8)-sphingolipid desaturase gene from Brassica rapa to evaluate the importance of these residues to the enzyme function. The mutants were functionally characterized by heterologous expression in yeast, allowing the identification of the products of the enzymes. The results revealed that residues H63, N203, D208, D210, and G368 were obligatorily required for the enzymatic activity, and substitution of the residues F59, W190, W345, L369 and Q372 markedly decreased the enzyme activity. Among them, replacement of the residues W190, L369 and Q372 also has significant influence on the ratio of the two enzyme products. Information obtained in this work provides the molecular basis for the Δ(8)-sphingolipid desaturase activity and aids in our understanding of the structure-function relationships of the membrane-bound desaturases.     

The role of cytochrome b5 fusion desaturases in the synthesis of polyunsaturated fatty acids

The biosynthetic pathway of polyunsaturated fatty acids (PUFAs) has been the subject of much interest over the last few years. Significant progress has been made in the identification of the enzymes required for PUFA synthesis; in particular, the fatty acid desaturases which are central to this pathway have now all been identified. These "front-end" desaturases are all members of the cytochrome b(5) fusion desaturase superfamily, since they contain an N-terminal domain that is orthologous to the microsomal cytochrome b(5). Examination of the primary sequence relationships between the various PUFA-specific cytochrome b(5) fusion desaturases and related fusion enzymes allows inferences regarding the evolution of this important enzyme class. More importantly, this knowledge helps underpin our understanding of polyunsaturated fatty acid biosynthesis.     

Identification of the substrate recognition region in the Δ6-fatty acid and Δ8-sphingolipid desaturase by fusion mutagenesis

Δ⁸-sphingolipid desaturase and Δ⁶-fatty acid desaturase share high protein sequence identity. Thus, it has been hypothesized that Δ⁶-fatty acid desaturase is derived from Δ⁸-sphingolipid desaturase; however, there is no direct proof. The substrate recognition regions of Δ⁶-fatty acid desaturase and Δ⁸-sphingolipid desaturase, which aid in understanding the evolution of these two enzymes, have not been reported. A blackcurrant Δ⁶-fatty acid desaturase and a Δ⁸-sphingolipid desaturase gene, RnD6C and RnD8A, respectively, share more than 80 % identity in their coding protein sequences. In this study, a set of fusion genes of RnD6C and RnD8A were constructed and expressed in yeast. The Δ⁶- and Δ⁸-desaturase activities of the fusion proteins were characterized. Our results indicated that (1) the exchange of the C-terminal 172 amino acid residues can lead to a significant decrease in both desaturase activities; (2) amino acid residues 114–174, 206–257, and 258–276 played important roles in Δ⁶-substrate recognition, and the last two regions were crucial for Δ⁸-substrate recognition; and (3) amino acid residues 114–276 of Δ⁶-fatty acid desaturase contained the substrate recognition site(s) responsible for discrimination between ceramide (a substrate of Δ⁸-sphingolipid desaturase) and acyl-PC (a substrate of Δ⁶-fatty acid desaturase). Substituting the amino acid residues 114-276 of RnD8A with those of RnD6C resulted in a gain of Δ⁶-desaturase activity in the fusion protein but a loss in Δ⁸-sphingolipid desaturase activity. In conclusion, several regions important for the substrate recognition of Δ8-sphingolipid desaturase and Δ⁶-fatty acid desaturase were identified, which provide clues in understanding the relationship between the structure and function in desaturases.