Identification and characterization of candidates involved in production of OMEGAs in microalgae: a gene mining and phylogenomic approach

Abstract

Optimizing the production of the high-value renewables such as OMEGAs through pathway engineering requires an in-depth understanding of the structure–function relationship of genes involved in the OMEGA biosynthetic pathways. In this preliminary study, our rationale is to identify and characterize the ～221 putative genes involved in production of OMEGAs using bioinformatic analysis from the Streptophyte (plants), Chlorophyte (green algae), Rhodophyta (red algae), and Bacillariophyta (diatoms) lineages based on their phylogenomic profiling, conserved motif/domain organization and physico-chemical properties. The MEME suite predicted 12 distinct protein domains, which are conserved among these putative genes. The phylogenomic analysis of the putative candidate genes [such as FAD2 (delta-12 desaturase); ECR (enoyl-CoA reductase); FAD2 (delta-12 desaturase); ACOT (acyl CoA thioesterase); ECH (enoyl-CoA hydratase); and ACAT (acetyl-CoA acyltransferase)] with similar domains and motif patterns were remarkably well conserved. Furthermore, the subcellular network prediction of OMEGA biosynthetic pathway genes revealed a unique interaction between the light-dependent chlorophyll biosynthesis and glycerol-3-phosphate dehydrogenase, which predicts a major cross-talk between the key essential pathways. Such bioinformatic analysis will provide insights in finding the key regulatory genes to optimize the productivity of OMEGAs in microalgal cell factories.     

Human brown fat inducible thioesterase variant 2 cellular localization and catalytic function

The mammalian brown fat inducible thioesterase variant 2 (BFIT2), also known as ACOT11, is a multimodular protein containing two consecutive hotdog-fold domains and a C-terminal steroidogenic acute regulatory protein-related lipid transfer domain (StarD14). In this study, we demonstrate that the N-terminal region of human BFIT2 (hBFIT2) constitutes a mitochondrial location signal sequence, which undergoes mitochondrion-dependent posttranslational cleavage. The mature hBFIT2 is shown to be located in the mitochondrial matrix, whereas the paralog "cytoplasmic acetyl-CoA hydrolase" (CACH, also known as ACOT12) was found in the cytoplasm. In vitro activity analysis of full-length hBFIT2 isolated from stably transfected HEK293 cells demonstrates selective thioesterase activity directed toward long chain fatty acyl-CoA thioesters, thus distinguishing the catalytic function of BFIT2 from that of CACH. The results from a protein-lipid overlay test indicate that the hBFIT2 StarD14 domain binds phosphatidylinositol 4-phosphate.     

Functional characterization of thioesterase superfamily member 1/Acyl-CoA thioesterase 11: implications for metabolic regulation

Thioesterase superfamily member 1 (Them1; synonyms acyl-CoA thioesterase 11 and StarD14) is highly expressed in brown adipose tissue and limits energy expenditure in mice. Them1 is a putative fatty acyl-CoA thioesterase that comprises tandem hot dog-fold thioesterase domains and a lipid-binding C-terminal steroidogenic acute regulatory protein-related lipid transfer (START) domain. To better define its role in metabolic regulation, this study examined the biochemical and enzymatic properties of Them1. Purified recombinant Them1 dimerized in solution to form an active fatty acyl-CoA thioesterase. Dimerization was induced by fatty acyl-CoAs, coenzyme A (CoASH), ATP, and ADP. Them1 hydrolyzed a range of fatty acyl-CoAs but exhibited a relative preference for long-chain molecular species. Thioesterase activity varied inversely with temperature, was stimulated by ATP, and was inhibited by ADP and CoASH. Whereas the thioesterase domains of Them1 alone were sufficient to yield active recombinant protein, the START domain was required for optimal enzyme activity. An analysis of subcellular fractions from mouse brown adipose tissue and liver revealed that Them1 contributes principally to the fatty acyl-CoA thioesterase activity of microsomes and nuclei. These findings suggest that under biological conditions, Them1 functions as a lipid-regulated fatty acyl-CoA thioesterase that could be targeted for the management of metabolic disorders.     

Direct observation of chemo-mechanical coupling in DnaK by single-molecule force experiments

Protein allostery requires a communication channel for functional regulation between distal sites within a protein. In the molecular chaperone Hsp70, a two-domain enzyme, the ATP/ADP status of an N-terminal nucleotide-binding domain regulates the substrate affinity of a C-terminal substrate-binding domain. Recently available three-dimensional structures of Hsp70 in ATP/ADP states have provided deep insights into molecular pathways of allosteric signals. However, direct mechanical probing of long-range allosteric coupling between the ATP hydrolysis step and domain states is missing. Using laser optical tweezers, we examined the mechanical properties of a truncated two-domain DnaK(1–552ye) in apo/ADP/ATP- and peptide-bound states. We find that in the apo and ADP states, DnaK domains are mechanically stable and rigid. However, in the ATP state, substrate-binding domain (SBD)¹ye is mechanically destabilized as the result of interdomain docking followed by the unfolding of the α-helical lid. By observing the folding state of the SBD, we could observe the continuous ATP/ADP cycling of the enzyme in real time with a single molecule. The SBD lid closure is strictly coupled to the chemical steps of the ATP hydrolysis cycle even in the presence of peptide substrate.     

Dissecting the N-Ethylmaleimide-sensitive Factor: REQUIRED ELEMENTS OF THE N AND D1 DOMAINS*

N-Ethylmaleimide-sensitive factor (NSF) is a homo-hexameric member of the AAA⁺ (ATPases associated with various cellular activities plus) family. It plays an essential role in most intracellular membrane trafficking through its binding to and disassembly of soluble NSF attachment protein (SNAP) receptor (SNARE) complexes. Each NSF protomer contains an N-terminal domain (NSF-N) and two AAA domains, a catalytic NSF-D1 and a structural NSF-D2. This study presents detailed mutagenesis analyses of NSF-N and NSF-D1, dissecting their roles in ATP hydrolysis, SNAP·SNARE binding, and complex disassembly. Our results show that a positively charged surface on NSF-N, bounded by Arg⁶⁷ and Lys¹⁰⁵, and the conserved residues in the central pore of NSF-D1 (Tyr²⁹⁶ and Gly²⁹⁸) are involved in SNAP·SNARE binding but not basal ATP hydrolysis. Mutagenesis of Sensor 1 (Thr³⁷³–Arg³⁷⁵), Sensor 2 (Glu⁴⁴⁰–Glu⁴⁴²), and Arginine Fingers (Arg³⁸⁵ and Arg³⁸⁸) in NSF-D1 shows that each region plays a discrete role. Sensor 1 is important for basal ATPase activity and nucleotide binding. Sensor 2 plays a role in ATP- and SNAP-dependent SNARE complex binding and disassembly but does so in cis and not through inter-protomer interactions. Arginine Fingers are important for SNAP·SNARE complex-stimulated ATPase activity and complex disassembly. Mutants at these residues have a dominant-negative phenotype in cells, suggesting that Arginine Fingers function in trans via inter-protomer interactions. Taken together, these data establish functional roles for many of the structural elements of the N domain and of the D1 ATP-binding site of NSF.     

Structure,function, and regulation of thioesterases

Thioesterases are present in all living cells and perform a wide range of important biological functions by catalysing the cleavage of thioester bonds present in a diverse array of cellular substrates. Thioesterases are organised into 25 families based on their sequence conservation, tertiary and quaternary structure, active site configuration, and substrate specificity. Recent structural and functional characterisation of thioesterases has led to significant changes in our understanding of the regulatory mechanisms that govern enzyme activity and their respective cellular roles. The resulting dogma changes in thioesterase regulation include mechanistic insights into ATP and GDP-mediated regulation by oligomerisation, the role of new key regulatory regions, and new insights into a conserved quaternary structure within TE4 family members. Here we provide a current and comparative snapshot of our understanding of thioesterase structure, function, and regulation across the different thioesterase families.     

Inflammatory stimuli induce acyl-CoA thioesterase 7 and remodeling of phospholipids containing unsaturated long (≥C20)-acyl chains in macrophages

Acyl-CoA thioesterase 7 (ACOT7) is an intracellular enzyme that converts acyl-CoAs to FFAs. ACOT7 is induced by lipopolysaccharide (LPS); thus, we investigated downstream effects of LPS-induced induction of ACOT7 and its role in inflammatory settings in myeloid cells. Enzymatic thioesterase activity assays in WT and ACOT7-deficient macrophage lysates indicated that endogenous ACOT7 contributes a significant fraction of total acyl-CoA thioesterase activity toward C20:4-, C20:5-, and C22:6-CoA, but contributes little activity toward shorter acyl-CoA species. Lipidomic analyses revealed that LPS causes a dramatic increase, primarily in bis(monoacylglycero)phosphate species containing long (≥C20) polyunsaturated acyl-chains in macrophages, and that the limited effect observed by ACOT7 deficiency is restricted to glycerophospholipids containing 20-carbon unsaturated acyl-chains. Furthermore, ACOT7 deficiency did not detectably alter the ability of LPS to induce cytokines or prostaglandin E₂ production in macrophages. Consistently, although ACOT7 was induced in macrophages from diabetic mice, hematopoietic ACOT7 deficiency did not alter the stimulatory effect of diabetes on systemic inflammation or atherosclerosis in LDL receptor-deficient mice. Thus, inflammatory stimuli induce ACOT7 and remodeling of phospholipids containing unsaturated long (≥C20)-acyl chains in macrophages, and, although ACOT7 has preferential thioesterase activity toward these lipid species, loss of ACOT7 has no major detrimental effect on macrophage inflammatory phenotypes.≥     

Structural insight into activation mechanism of Toxoplasma gondii nucleoside triphosphate diphosphohydrolases by disulfide reduction

The intracellular parasite Toxoplasma gondii produces two nucleoside triphosphate diphosphohydrolases (NTPDase1 and -3). These tetrameric, cysteine-rich enzymes require activation by reductive cleavage of a hitherto unknown disulfide bond. Despite a 97% sequence identity, both isozymes differ largely in their ability to hydrolyze ATP and ADP. Here, we present crystal structures of inactive NTPDase3 as an apo form and in complex with the product AMP to resolutions of 2.0 and 2.2 Å, respectively. We find that the enzyme is present in an open conformation that precludes productive substrate binding and catalysis. The cysteine bridge 258–268 is identified to be responsible for locking of activity. Crystal structures of constitutively active variants of NTPDase1 and -3 generated by mutation of Cys²⁵⁸–Cys²⁶⁸ show that opening of the regulatory cysteine bridge induces a pronounced contraction of the whole tetramer. This is accompanied by a 12° domain closure motion resulting in the correct arrangement of all active site residues. A complex structure of activated NTPDase3 with a non-hydrolyzable ATP analog and the cofactor Mg²⁺ to a resolution of 2.85 Å indicates that catalytic differences between the NTPDases are primarily dictated by differences in positioning of the adenine base caused by substitution of Arg⁴⁹² and Glu⁴⁹³ in NTPDase1 by glycines in NTPDase3.     

Biochemical and Kinetic Characterization of the Recombinant ADP-Forming Acetyl Coenzyme A Synthetase from the Amitochondriate Protozoan Entamoeba histolytica

Entamoeba histolytica, an amitochondriate protozoan parasite that relies on glycolysis as a key pathway for ATP generation, has developed a unique extended PP_i-dependent glycolytic pathway in which ADP-forming acetyl-coenzyme A (CoA) synthetase (ACD; acetate:CoA ligase [ADP-forming]; EC 6.2.1.13) converts acetyl-CoA to acetate to produce additional ATP and recycle CoA. We characterized the recombinant E. histolytica ACD and found that the enzyme is bidirectional, allowing it to potentially play a role in ATP production or in utilization of acetate. In the acetate-forming direction, acetyl-CoA was the preferred substrate and propionyl-CoA was used with lower efficiency. In the acetyl-CoA-forming direction, acetate was the preferred substrate, with a lower efficiency observed with propionate. The enzyme can utilize both ADP/ATP and GDP/GTP in the respective directions of the reaction. ATP and PP_i were found to inhibit the acetate-forming direction of the reaction, with 50% inhibitory concentrations of 0.81 ± 0.17 mM (mean ± standard deviation) and 0.75 ± 0.20 mM, respectively, which are both in the range of their physiological concentrations. ATP and PP_i displayed mixed inhibition versus each of the three substrates, acetyl-CoA, ADP, and phosphate. This is the first example of regulation of ACD enzymatic activity, and possible roles for this regulation are discussed.     

Biochemical and Kinetic Characterization of the Eukaryotic Phosphotransacetylase Class IIa Enzyme from Phytophthora ramorum

Phosphotransacetylase (Pta), a key enzyme in bacterial metabolism, catalyzes the reversible transfer of an acetyl group from acetyl phosphate to coenzyme A (CoA) to produce acetyl-CoA and P_i. Two classes of Pta have been identified based on the absence (Pta^I) or presence (Pta^II) of an N-terminal regulatory domain. Pta^I has been fairly well studied in bacteria and one genus of archaea; however, only the Escherichia coli and Salmonella enterica Pta^II enzymes have been biochemically characterized, and they are allosterically regulated. Here, we describe the first biochemical and kinetic characterization of a eukaryotic Pta from the oomycete Phytophthora ramorum. The two Ptas from P. ramorum, designated PrPta^II1 and PrPta^II2, both belong to class II. PrPta^II1 displayed positive cooperativity for both acetyl phosphate and CoA and is allosterically regulated. We compared the effects of different metabolites on PrPta^II1 and the S. enterica Pta^II and found that, although the N-terminal regulatory domains share only 19% identity, both enzymes are inhibited by ATP, NADP, NADH, phosphoenolpyruvate (PEP), and pyruvate in the acetyl-CoA/P_i-forming direction but are differentially regulated by AMP. Phylogenetic analysis of bacterial, archaeal, and eukaryotic sequences identified four subtypes of Pta^II based on the presence or absence of the P-loop and DRTGG subdomains within the N-terminal regulatory domain. Although the E. coli, S. enterica, and P. ramorum enzymes all belong to the IIa subclass, our kinetic analysis has indicated that enzymes within a subclass can still display differences in their allosteric regulation.     

Characterization of the molecular properties of KtrC,a second RCK domain that regulates a Ktr channel in Bacillus subtilis

RCK (regulating conductance of K⁺) domains are common regulatory domains that control the activity of eukaryotic and prokaryotic K⁺ channels and transporters. In bacteria these domains play roles in osmoregulation, regulation of turgor and membrane potential and in pH homeostasis. Whole-genome sequencing unveiled RCK gene redundancy, however the biological role of this redundancy is not well understood. In Bacillus subtilis, there are two closely related RCK domain proteins (KtrA and KtrC) that regulate the activity of the Ktr cation channels. KtrA has been well characterized but little is known about KtrC. We have characterized the structural and biochemical proprieties of KtrC and conclude that KtrC binds ATP and ADP, just like KtrA. However, in solution KtrC exist in a dynamic equilibrium between octamers and non-octameric species that is dependent on the bound ligand, with ATP destabilizing the octameric ring relative to ADP. Accordingly, KtrC-ADP crystal structures reveal closed octameric rings similar to those in KtrA, while KtrC-ATP adopts an open assembly with RCK domains forming a super-helix. In addition, both KtrC-ATP and -ADP octamers are stabilized by the signaling molecule cyclic-di-AMP, which binds to KtrC with high affinity. In contrast, c-di-AMP binds with 100-fold lower affinity to KtrA. Despite these differences we show with an E. coli complementation assay that KtrC and KtrA are interchangeable and able to form functional transporters with both KtrB and KtrD. The distinctive properties of KtrC, in particular ligand-dependent assembly/disassembly, suggest that this protein has a specific physiological role that is distinct from KtrA.     

Effects of adenosine phosphates and nicotinamide nucleotides on pyruvate carboxylase from baker''s yeast

1. Pyruvate carboxylase from baker's yeast is inhibited by ADP, AMP and adenosine at pH8.0 in the presence of magnesium chloride concentrations equal to or higher than the ATP concentration. The adenine moiety is essential for the inhibitory effect. 2. In the absence of acetyl-CoA (an allosteric activator) ADP, AMP and adenosine are competitive inhibitors with respect to ATP. In the presence of acetyl-CoA, besides the effect with respect to ATP, AMP competes with acetyl-CoA, whereas ADP and adenosine are non-competitive inhibitors with respect to the activator. 3. Pyruvate carboxylase is inhibited by NADH. The inhibition is competitive with respect to acetyl-CoA and specific with respect to NADH, since NAD(+), NADP(+) and NADPH do not affect the enzyme activity. In the absence of acetyl-CoA, NAD(+), NADH, NADP(+) and NADPH do not inhibit pyruvate carboxylase. 4. Pyruvate carboxylase is inhibited by ADP, AMP and NADH at pH6.5, in the presence of 12mm-Mg(2+), 0.75mm-Mn(2+) and 0.5mm-ATP, medium conditions similar to those existing inside the yeast cell. The ADP and NADH effects are consistent with a regulation of enzyme activity by the intracellular [ATP]/[ADP] ratio and secondarily by NADH concentration. These mechanisms would supplement the already known control of yeast pyruvate carboxylase by acetyl-CoA and l-aspartate. Inhibition by AMP is less marked and its physiological role is perhaps limited.     

海洋链霉菌酰基辅酶A硫酯酶基因的克隆及分析

酰基辅酶A硫酯酶(ACOT)是催化脂酰辅酶A水解成自由脂肪酸和辅酶A的一类酶,而II型ACOT对底物具有较高的特异性,在脂肪酸合成途径中起着重要的作用,II型ACOT的缺乏或失调会导致机体脂肪酸代谢的紊乱,从而引起一系列疾病。从海洋链霉菌L1基因组DNA中克隆到786碱基的ACOT II基因,生物信息学分析表明其拟编码262个氨基酸,与Acyl-CoA thioesterase II (sequence ID:WP 069742521.1)相似性达到100%。BLAST比对发现,其属于热狗超家族成员,并包含一个II型硫酯酶催化反应的特有结构域。并构建了表达载体pET32a-ACOT II,0.5 mmol/L IPTG诱导目的基因表达,SDS-PAGE显示其分子量约为37.0 kD。通过对链霉菌中II型ACOT基因的研究,可为进一步深入研究II型硫酯酶的分子结构和生物学功能,进而指导药物研发提供参考。     

Redox-linked Gating of Nucleotide Binding by the N-terminal Domain of Adenosine 5′-Phosphosulfate Kinase

Adenosine 5′-phosphosulfate kinase (APSK) catalyzes the phosphorylation of adenosine 5′-phosphosulfate (APS) to 3′-phosphoadenosine-5′-phosphosulfate (PAPS). Crystallographic studies of APSK from Arabidopsis thaliana revealed the presence of a regulatory intersubunit disulfide bond (Cys⁸⁶–Cys¹¹⁹). The reduced enzyme displayed improved catalytic efficiency and decreased effectiveness of substrate inhibition by APS compared with the oxidized form. Here we examine the effect of disulfide formation and the role of the N-terminal domain on nucleotide binding using isothermal titration calorimetry (ITC) and steady-state kinetics. Formation of the disulfide bond in A. thaliana APSK (AtAPSK) inverts the binding affinities at the ATP/ADP and APS/PAPS sites from those observed in the reduced enzyme, consistent with initial binding of APS as inhibitory, and suggests a role for the N-terminal domain in guiding nucleotide binding order. To test this, an N-terminal truncation variant (AtAPSKΔ96) was generated. The resulting protein was completely insensitive to substrate inhibition by APS. ITC analysis of AtAPSKΔ96 showed decreased affinity for APS binding, although the N-terminal domain does not directly interact with this ligand. Moreover, AtAPSKΔ96 displayed reduced affinity for ADP, which corresponds to a loss of substrate inhibition by formation of an E·ADP·APS dead end complex. Examination of the AtAPSK crystal structure suggested Arg⁹³ as important for positioning of the N-terminal domain. ITC and kinetic analysis of the R93A mutant also showed a complete loss of substrate inhibition and altered nucleotide binding affinities, which mimics the effect of the N-terminal deletion. These results show how thiol-linked changes in AtAPSK alter the energetics of binding equilibria to control its activity.     

Regulation of Energy Transduction and Electron Transfer in Cytochrome c Oxidase by Adenine Nucleotides

Cytochrome c oxidase from bovine heart contains seven high-affinity binding sites for ATP or ADP and three additional only for ADP. One binding site for ATP or ADP, located at the matrix-oriented domain of the heart-type subunit VIaH, increases the H⁺/e^– stoichiometry of the enzyme from heart or skeletal muscle from 0.5 to 1.0 when bound ATP is exchanged by ADP. Two further binding sites for ATP or ADP, located at the cytosolic and the matrix domain of subunit IV, increases the K _M for Cytochrome c and inhibit the respiratory activity at high ATP/ADP ratios, respectively. We propose that thermogenesis in mammals is related to subunit VIaL of cytochrome c oxidase with a H⁺/e^– stoichiometry of 0.5 compared to 1.0 in the enzyme from bacteria or ectotherm animals. This hypothesis is supported by the lack of subunit VIa isoforms in cytochrome c oxidase from fish.     

Probing l-Pyruvate Kinase Regulatory Phosphorylation Site by Mutagenesis

The activity of L-type pyruvate kinase (L-PK, ATP:pyruvate 2-O-phosphotransferase, EC 2.7.1.40) is regulated by phosphorylation of serine residue 12 of the N-terminal regulatory domain MEGPAGYLRR¹⁰AS ¹² VAQLTQEL²⁰GTAFF of the protein. In this report we studied the effect of the point mutations around this phosphorylation site on the catalytic properties of this enzyme, by introducing amino acids A, L, K, Q and E into positions 9, 10 and 13 of this peptide sequence. It was found that some of these mutations in positions 9 and 10, although occurring at great distances from the enzyme??s active site, affected the enzyme??s activity by decreasing the effectiveness of phosphoenolpyruvate binding (PEP) with the enzyme, but had practically no influence on the binding effectiveness of the second substrate ADP. A similar asymmetric effect on the binding of these substrates was previously observed after phosphorylation of the enzyme regulatory N-domain peptide, and also after proteolytic truncation of the same N-terminal part of L-PK. All these results could be explained by the internal complex formation between the N-domain peptide and the enzyme??s main body. The present study delineated the specificity of the internal binding site and revealed the possibility that the regulatory effect could be modulated by selecting mutation sites and amino acids introduced into the N-terminal domain structure.     

Active creatine kinase is present in matrix vesicles isolated from femurs of chicken embryo: Implications for bone mineralization

Proteomic analysis of matrix vesicles (MVs) isolated from 17-day-old chicken embryo femurs revealed the presence of creatine kinase. In this report we identified the enzyme functionally and suggest that the enzyme may participate in the synthesis of ATP from ADP and phosphocreatine within the lumen of these organelles. Then, ATP is converted by nucleotide hydrolyzing enzymes such as Na⁺, K⁺-ATPase, protein kinase C, or alkaline phosphatase to yield inorganic phosphate (P_i), a substrate for mineralization. Alternatively, ATP can be hydrolyzed by a nucleoside triphosphate pyrophosphatase phosphodiesterase 1 producing inorganic pyrophosphate (PP_i), a mineralization inhibitor. In addition, immunochemical evidence indicated that VDAC 2 is present in MVs that may serve as a transporter of nucleotides from the extracellular matrix. We discussed the implications of ATP production and hydrolysis by MVs as regulatory mechanisms for mineralization.