Identification of operons in the hyperthermophilic archaeon Pyrococcus furiosus represents an important step to understanding the regulatory mechanisms that enable the organism to adapt and thrive in extreme environments. We have predicted operons in P.furiosus by combining the results from three existing algorithms using a neural network (NN). These algorithms use intergenic distances, phylogenetic profiles, functional categories and gene-order conservation in their operon prediction. Our method takes as inputs the confidence scores of the three programs, and outputs a prediction of whether adjacent genes on the same strand belong to the same operon. In addition, we have applied Gene Ontology (GO) and KEGG pathway information to improve the accuracy of our algorithm. The parameters of this NN predictor are trained on a subset of all experimentally verified operon gene pairs of Bacillus subtilis. It subsequently achieved 86.5% prediction accuracy when applied to a subset of gene pairs for Escherichia coli, which is substantially better than any of the three prediction programs. Using this new algorithm, we predicted 470 operons in the P.furiosus genome. Of these, 349 were validated using DNA microarray data.

Porcine epidemic diarrhea virus (PEDV), an enteric coronavirus, induces severe vomiting and acute watery diarrhea in unweaned piglets. The pig industry has suffered tremendous financial losses due to the high mortality rate of piglets caused by PEDV. Consequently, a simple and rapid on-site diagnostic technology is crucial for preventing and controlling PEDV. This study established a detection method for PEDV using recombinase-aided amplification (RAA) and Pyrococcus furiosus Argonaute (PfAgo), which can detect 100 copies of PEDV without cross-reactivity with other pathogens. The entire reaction of RAA and PfAgo to detect PEDV does not require sophisticated instruments, and the reaction results can be observed with the naked eye. Overall, this integrated RAA-PfAgo cleavage assay is a practical tool for accurately and quickly detecting PEDV. KEY POINTS: • PfAgo has the potential to serve as a viable molecular diagnostic tool for the detection and diagnosis of viral genomes • The RAA-PfAgo detection technique has a remarkable level of sensitivity and specificity • The RAA-PfAgo detection system can identify PEDV without needing advanced equipment.

Synthesis and assembly of ribosomal components are fundamental cellular processes and generally well-conserved within the main groups of organisms. Yet, provocative variations to the general schemes exist. We have discovered an unusual processing pathway of pre-rRNA in extreme thermophilic archaea exemplified by Pyrococcus furiosus. The large subunit (LSU) rRNA is produced as a circularly permuted form through circularization followed by excision of Helix 98. As a consequence, the terminal domain VII that comprise the binding site for the signal recognition particle is appended to the 5´ end of the LSU rRNA that instead terminates in Domain VI carrying the Sarcin-Ricin Loop, the primary interaction site with the translational GTPases. To our knowledge, this is the first example of a true post-transcriptional circular permutation of a main functional molecule and the first example of rRNA fragmentation in archaea.

Dehydroquinate synthase (DHQS) catalyses the second step of the shikimate pathway to aromatic compounds. DHQS from the archaeal hyperthermophile Pyrococcus furiosus was insoluble when expressed in Escherichia coli but was partially solubilised when KCl was included in the cell lysis buffer. A purification procedure was developed, involving lysis by sonication at 30°C followed by a heat treatment at 70°C and anion exchange chromatography. Purified recombinant P. furiosus DHQS is a dimer with a subunit Mr of 37,397 (determined by electrospray ionisation mass spectrometry) and is active over broad pH and temperature ranges. The kinetic parameters are K(M) (3-deoxy-D-arabino-heptulosonate 7-phosphate) 3.7 μM and k(cat) 3.0 sec(-1) at 60°C and pH 6.8. EDTA inactivates the enzyme, and enzyme activity is restored by several divalent metal ions including (in order of decreasing effectiveness) Cd(2+), Co(2+), Zn(2+), and Mn(2+). High activity of a DHQS in the presence of Cd(2+) has not been reported for enzymes from other sources, and may be related to the bioavailability of Cd(2+) for P. furiosus. This study is the first biochemical characterisation of a DHQS from a thermophilic source. Furthermore, the characterisation of this hyperthermophilic enzyme was carried out at elevated temperatures using an enzyme-coupled assay.

The methyltransferase fibrillarin is the catalytic component of ribonucleoprotein complexes that direct site-specific methylation of precursor ribosomal RNA and are critical for ribosome biogenesis in eukaryotes and archaea. Here we report the crystal structure of a fibrillarin ortholog from the hyperthermophilic archaeon Pyrococcus furiosus at 1.97A resolution. Comparisons of the X-ray structures of fibrillarin orthologs from Methanococcus jannashii and Archaeoglobus fulgidus reveal nearly identical backbone configurations for the catalytic C-terminal domain with the exception of a unique loop conformation at the S-adenosyl-l-methionine (AdoMet) binding pocket in P. furiosus. In contrast, the N-terminal domains are divergent which may explain why some forms of fibrillarin apparently homodimerize (M. jannashii) while others are monomeric (P. furiosus and A. fulgidus). Three positively charged amino acids surround the AdoMet-binding site and sequence analysis indicates that this is a conserved feature of both eukaryotic and archaeal fibrillarins. We discuss the possibility that these basic residues of fibrillarin are important for RNA-guided rRNA methylation.

The archaellum is the macromolecular machinery that Archaea use for propulsion or surface adhesion, enabling them to proliferate and invade new territories. The molecular composition of the archaellum and of the motor that drives it appears to be entirely distinct from that of the functionally equivalent bacterial flagellum and flagellar motor. Yet, the structure of the archaellum machinery is scarcely known. Using combined modes of electron cryo-microscopy (cryoEM), we have solved the structure of the Pyrococcus furiosus archaellum filament at 4.2 Å resolution and visualise the architecture and organisation of its motor complex in situ. This allows us to build a structural model combining the archaellum and its motor complex, paving the way to a molecular understanding of archaeal swimming motion.

The crystal structure of the ferritin from the archaeon, hyperthermophile and anaerobe Pyrococcus furiosus (PfFtn) is presented. While many ferritin structures from bacteria to mammals have been reported, until now only one was available from archaea, the ferritin from Archaeoglobus fulgidus (AfFtn). The PfFtn 24-mer exhibits the 432 point-group symmetry that is characteristic of most ferritins, which suggests that the 23 symmetry found in the previously reported AfFtn is not a common feature of archaeal ferritins. Consequently, the four large pores that were found in AfFtn are not present in PfFtn. The structure has been solved by molecular replacement and refined at 2.75-Angstrom resolution to R = 0.195 and R(free) = 0.247. The ferroxidase center of the aerobically crystallized ferritin contains one iron at site A and shows sites B and C only upon iron or zinc soaking. Electron paramagnetic resonance studies suggest this iron depletion of the native ferroxidase center to be a result of a complexation of iron by the crystallization salt. The extreme thermostability of PfFtn is compared with that of eight structurally similar ferritins and is proposed to originate mostly from the observed high number of intrasubunit hydrogen bonds. A preservation of the monomer fold, rather than the 24-mer assembly, appears to be the most important factor that protects the ferritin from inactivation by heat.

Because base deaminations, which are promoted by high temperature, ionizing radiation, aerobic respiration and nitrosative stress, produce mutations during replication, deaminated bases must be repaired quickly to maintain genome integrity. Recently, we identified a novel lesion-specific endonuclease, PfuEndoQ, from Pyrococcus furiosus, and PfuEndoQ may be involved in the DNA repair pathway in Thermococcales of Archaea. PfuEndoQ recognizes a deaminated base and cleaves the phosphodiester bond 5' of the lesion site. To elucidate the structural basis of the substrate recognition and DNA cleavage mechanisms of PfuEndoQ, we determined the structure of PfuEndoQ using X-ray crystallography. The PfuEndoQ structure and the accompanying biochemical data suggest that PfuEndoQ recognizes a deaminated base using a highly conserved pocket adjacent to a Zn2+-binding site and hydrolyses a phosphodiester bond using two Zn2+ ions. The PfuEndoQ-DNA complex is stabilized by a Zn-binding domain and a C-terminal helical domain, and the complex may recruit downstream proteins in the DNA repair pathway.

The bacterial acquired immune system consists of clustered regularly interspaced short palindromic repeats (CRISPRs) and CRIPSR-associated (Cas) genes, which include Cas-module repeat-associated mysterious proteins (Cmr). The six Cmr proteins of Pyrococcus furiosus (pfCmr1-pfCmr6) form a Cmr effector complex that functions against exogenous nucleic acid. Among the Cmr proteins, the role of pfCmr5 and its involvement in the complex's cleavage activity have been obscure. The elucidated pfCmr5 structure has two inserted α-helices compared with the other trimeric Cmr5 structure. However, pfCmr5 exists as a monomeric protein both in the crystalline state and in solution. In vitro assays indicate that pfCmr5 interacts with pfCmr4. These structural and biophysical data might help in understanding the complicated and ill-characterized Cmr effector complex.

Although copper is in many cases an essential micronutrient for cellular life, higher concentrations are toxic. Therefore, all living cells have developed strategies to maintain copper homeostasis. In this manuscript, we have analyzed the transcriptome-wide response of Pyrococcus furiosus to increased copper concentrations and described the essential role of the putative copper-sensing metalloregulator CopR in the detoxification process. To this end, we employed biochemical and biophysical methods to characterize the role of CopR. Additionally, a copR knockout strain revealed an amplified sensitivity in comparison to the parental strain towards increased copper levels, which designates an essential role of CopR for copper homeostasis. To learn more about the CopR-regulated gene network, we performed differential gene expression and ChIP-seq analysis under normal and 20 μM copper-shock conditions. By integrating the transcriptome and genome-wide binding data, we found that CopR binds to the upstream regions of many copper-induced genes. Negative-stain transmission electron microscopy and 2D class averaging revealed an octameric assembly formed from a tetramer of dimers for CopR, similar to published crystal structures from the Lrp family. In conclusion, we propose a model for CopR-regulated transcription and highlight the regulatory network that enables Pyrococcus to respond to increased copper concentrations.

The hyperthermophilic archaeon Pyrococcus furiosus grows by fermenting peptides and carbohydrates to organic acids. In the terminal step, acyl-CoA synthetase (ACS) isoenzymes convert acyl-CoA derivatives to the corresponding acid and conserve energy in the form of ATP. ACS1 and ACS2 were previously purified from P. furiosus and have α 2 β 2 structures but the genome contains genes encoding three additional α-subunits. The ten possible combinations of α and β genes were expressed in E. coli and each resulted in stable and active α 2 β 2 isoenzymes. The α-subunit of each isoenzyme determined CoA-based substrate specificity and between them they accounted for the CoA derivatives of fourteen amino acids. The β-subunit determined preference for adenine or guanine nucleotides. The GTP-generating isoenzymes are proposed to play a role in gluconeogenesis by producing GTP for GTP-dependent phosphoenolpyruvate carboxykinase and for other GTP-dependent processes. Transcriptional and proteomic data showed that all ten isoenzymes are constitutively expressed indicating that both ATP and GTP are generated from the metabolism of most of the amino acids. A phylogenetic analysis showed that the ACSs of P. furiosus and other members of the Thermococcales are evolutionarily distinct from those found throughout the rest of biology, including those of other hyperthermophilic archaea.

Several in vitro studies document the function of the transcriptional regulator TrmBL1 of Pyrococcus furiosus. These data indicate that the protein can act as repressor or activator and is mainly involved in transcriptional control of sugar uptake and in the switch between glycolysis and gluconeogenesis. The aim of this study was to complement the in vitro data with an in vivo analysis using ChIP-seq to explore the genome-wide binding profile of TrmBL1 under glycolytic and gluconeogenic growth conditions.

The need to find more stable catalysts has encouraged the study of naturally resilient enzymes, such as those found in extremophile organisms. In the present work, the influence of rare codons on the expression in Escherichia coli of the lipase Pf2001Δ60 from Pyrococcus furiosus was evaluated. Expression was carried out in two E. coli strains, BL21(DE3)pLysS and the rare tRNA supplemented Rosetta(DE3)pLysS. 3(2) factorial design was used to appraise the influence of temperature and inducer concentration on enzyme expression every hour for the 4-h expression period. Four response surfaces were constructed for each time, and the statistical parameters were evaluated. Lipase production was twice as high for Rosetta(DE3)pLysS than for BL21(DE3)pLysS. The factorial design indicated that optimal expression occurred at 30 °C after 4h, with lipase production of 240 U/L. The analysis of statistical parameters during the expression time showed that the velocity at which the enzyme was produced affected cell growth and maximum activity, with a higher speed leading to lower expression and cell growth. The presence of rare tRNAs prevented bottlenecks in lipase expression, and the experimental design was shown to be important for maximizing the production strategies and minimizing the metabolic load to which the host is subjected.

DNA is constantly damaged by endogenous and environmental influences. Deaminated adenine (hypoxanthine) tends to pair with cytosine and leads to the A:T→G:C transition mutation during DNA replication. Endonuclease V (EndoV) hydrolyzes the second phosphodiester bond 3' from deoxyinosine in the DNA strand, and was considered to be responsible for hypoxanthine excision repair. However, the downstream pathway after EndoV cleavage remained unclear. The activity to cleave the phosphodiester bond 5' from deoxyinosine was detected in a Pyrococcus furiosus cell extract. The protein encoded by PF1551, obtained from the mass spectrometry analysis of the purified fraction, exhibited the corresponding cleavage activity. A putative homolog from Thermococcus kodakarensis (TK0887) showed the same activity. Further biochemical analyses revealed that the purified PF1551 and TK0887 proteins recognize uracil, xanthine and the AP site, in addition to hypoxanthine. We named this endonuclease Endonuclease Q (EndoQ), as it may be involved in damaged base repair in the Thermococcals of Archaea.

Transcription factor B recruiting factor 1 (TFB-RF1; PF1088) is a transcription regulator which activates transcription on archaeal promoters containing weak TFB recognition elements (BRE) by recruiting TFB to the promoter. The mechanism of activation is described in detail, but nothing is known about the biological function of this protein in Pyrococcus furiosus. The protein is located in an operon structure together with the hypothetical gene pf1089 and western blot as well as end-point RT-PCR experiments revealed an extremely low expression rate of both proteins. Furthermore, conditions to induce the expression of the operon are not known. By introducing an additional copy of tfb-RF1 using a Pyrococcus shuttle vector we could circumvent the lacking expression of both proteins under standard growth conditions as indicated by western blot as well as end-point RT-PCR experiments. A ChIP-seq experiment revealed an additional binding site of TFB-RF1 in the upstream region of the pf1011/1012 operon, beside the expected target of the pf1089/tfb-RF1 region. This operon codes for a putative ABC transporter which is most-related to a multidrug export system and in vitro analysis using gel shift assays, DNase I footprinting and in vitro transcription confirmed the activator function of TFB-RF1 on the corresponding promoter. These findings are also in agreement with in vivo data, as RT-qPCR experiments also indicate transcriptional activation of both operons. Taken together, the overexpression strategy of tfb-RF1 enabled the identification of an additional operon of the TFB-RF1 regulon which indicates a transport-related function and provides a promising starting position to decipher the physiological function of the TFB-RF1 gene regulatory network in P. furiosus.

Alpha amylase hydrolyzes α-bonds of polysaccharides such as starch and produces malto-oligosaccharides. Its starch saccharification applications make it an essential enzyme in the textile, food and brewing industries. Commercially available α-amylase is mostly produced from Bacillus or Aspergillus. A hyper-thermostable and Ca 2++ independent α-amylase from Pyrococcus furiosus (PFA) expressed in E.coli forms insoluble inclusion bodies and thus is not feasible for industrial applications.

Although Pyrococcus furiosus is one of the best studied hyperthermophilic archaea, to date no experimental investigation of the extent of protein secretion has been performed. We describe experimental verification of the extracellular proteome of P. furiosus grown on starch. LC-MS/MS-based analysis of culture supernatants led to the identification of 58 proteins. Fifteen of these proteins had a putative N-terminal signal peptide (SP), tagging the proteins for translocation across the membrane. The detected proteins with predicted SPs and known function were almost exclusively involved in important extracellular functions, like substrate degradation or transport. Most of the 43 proteins without predicted N-terminal signal sequences are known to have intracellular functions, mainly (70 %) related to intracellular metabolism. In silico analyses indicated that the genome of P. furiosus encodes 145 proteins with N-terminal SPs, including 21 putative lipoproteins and 17 with a class III peptide. From these we identified 15 (10 %; 7 SPI, 3 SPIII and 5 lipoproteins) under the specific growth conditions of this study. The putative lipoprotein signal peptides have a unique sequence motif, distinct from the motifs in bacteria and other archaeal orders.

Ethanol is an important target for the renewable production of liquid transportation fuels. It can be produced biologically from pyruvate, via pyruvate decarboxylase, or from acetyl-CoA, by alcohol dehydrogenase E (AdhE). Thermophilic bacteria utilize AdhE, which is a bifunctional enzyme that contains both acetaldehyde dehydrogenase and alcohol dehydrogenase activities. Many of these organisms also contain a separate alcohol dehydrogenase (AdhA) that generates ethanol from acetaldehyde, although the role of AdhA in ethanol production is typically not clear. As acetyl-CoA is a key central metabolite that can be generated from a wide range of substrates, AdhE can serve as a single gene fuel module to produce ethanol through primary metabolic pathways. The focus here is on the hyperthermophilic archaeon Pyrococcus furiosus, which grows by fermenting sugar to acetate, CO2 and H2 . Previously, by the heterologous expression of adhA from a thermophilic bacterium, P. furiosus was shown to produce ethanol by a novel mechanism from acetate, mediated by AdhA and the native enzyme aldehyde oxidoreductase (AOR). In this study, the AOR gene was deleted from P. furiosus to evaluate ethanol production directly from acetyl-CoA by heterologous expression of the adhE gene from eight thermophilic bacteria. Only AdhEs from two Thermoanaerobacter strains showed significant activity in cell-free extracts of recombinant P. furiosus and supported ethanol production in vivo. In the AOR deletion background, the highest amount of ethanol (estimated 61% theoretical yield) was produced when adhE and adhA from Thermoanaerobacter were co-expressed.

Nucleases play important roles in all DNA transactions, including replication, repair, and recombination. Many different nucleases from bacterial and eukaryotic organisms have been identified and functionally characterized. However, our knowledge about the nucleases from Archaea, the third domain of life, is still limited. We searched for 3'-5' exonuclease activity in the hyperthermophilic archaeon, Pyrococcus furiosus, and identified a protein with the target activity. The purified protein, encoded by PF2046, is composed of 229 amino acids with a molecular weight of 25,596, and displayed single-strand specific 3'-5' exonuclease activity. The protein, designated as PfuExo I, forms a stable trimeric complex in solution and excises the DNA at every two nucleotides from the 3' to 5' direction. The amino acid sequence of this protein is conserved only in Thermococci, one of the hyperthermophilic classes in the Euryarchaeota subdomain in Archaea. The newly discovered exonuclease lacks similarity to any other proteins with known function, including hitherto reported 3'-5' exonucleases. This novel nuclease may be involved in a DNA repair pathway conserved in the living organisms as a specific member for some hyperthermophilic archaea.

Pyrococcus furiosus PF2050 is an uncharacterized putative protein that contains two DUF2666 domains. Functional and structural studies of PF2050 have not previously been performed. In this study, we determined the crystal structure of PF2050. The structure of PF2050 showed that the two DUF2666 domains interact tightly, forming a globular structure. Each DUF2666 domain comprises an antiparallel β-sheet and an α-helical bundle. One side of the PF2050 structure has a positively charged basic cleft, which may have a DNA-binding function. Furthermore, we confirmed that PF2050 interacts with circular and linear dsDNA.