A Cytoplasmic NAD(P)H-Dependent Polysulfide Reductase with Thiosulfate Reductase Activity from the Hyperthermophilic Bacterium Thermotoga maritima

Thermotoga maritima is an anaerobic hyperthermophilic bacterium that efficiently produces H2 by fermenting carbohydrates. High concentration of H2 inhibits the growth of T. maritima, and S0 could eliminate the inhibition and stimulate the growth through its reduction. The mechanism of T. maritima sulfur reduction, however, has not been fully understood. Herein, based on its similarity with archaeal NAD(P)H-dependent sulfur reductases (NSR), the ORF THEMA_RS02810 was identified and expressed in Escherichia coli, and the recombinant protein was characterized. The purified flavoprotein possessed NAD(P)H-dependent S0 reductase activity (1.3 U/mg for NADH and 0.8 U/mg for NADPH), polysulfide reductase activity (0.32 U/mg for NADH and 0.35 U/mg for NADPH), and thiosulfate reductase activity (2.3 U/mg for NADH and 2.5 U/mg for NADPH), which increased 3~4-folds by coenzyme A stimulation. Quantitative RT-PCR analysis showed that nsr was upregulated together with the mbx, yeeE, and rnf genes when the strain grew in S0- or thiosulfate-containing medium. The mechanism for sulfur reduction in T. maritima was discussed, which may affect the redox balance and energy metabolism of T. maritima. Genome search revealed that NSR homolog is widely distributed in thermophilic bacteria and archaea, implying its important role in the sulfur cycle of geothermal environments. IMPORTANCE The reduction of S0 and thiosulfate is essential in the sulfur cycle of geothermal environments, in which thermophiles play an important role. Despite previous research on some sulfur reductases of thermophilic archaea, the mechanism of sulfur reduction in thermophilic bacteria is still not clearly understood. Herein, we confirmed the presence of a cytoplasmic NAD(P)H-dependent polysulfide reductase (NSR) from the hyperthermophile T. maritima, with S0, polysulfide, and thiosulfate reduction activities, in contrast to other sulfur reductases. When grown in S0- or thiosulfate-containing medium, its expression was upregulated. And the putative membrane-bound MBX and Rnf may also play a role in the metabolism, which might influence the redox balance and energy metabolism of T. maritima. This is distinct from the mechanism of sulfur reduction in mesophiles such as Wolinella succinogenes. NSR homologs are widely distributed among heterotrophic thermophiles, suggesting that they may be vital in the sulfur cycle in geothermal environments.

High-temperature behavior of hyperthermostable Thermotoga maritima xylanase XYN10B after designed and evolved mutations

A hyperthermostable xylanase XYN10B from Thermotoga maritima (PDB code 1VBR, GenBank accession number KR078269) was subjected to site-directed and error-prone PCR mutagenesis. From the selected five mutants, the two site-directed mutants (F806H and F806V) showed a 3.3-3.5-fold improved enzyme half-life at 100 °C. The mutant XYNA generated by error-prone PCR showed slightly improved stability at 100 °C and a lower K_m. In XYNB and XYNC, the additional mutations over XYNA decreased the thermostability and temperature optimum, while elevating the K_m. In XYNC, two large side-chains were introduced into the protein's interior. Micro-differential scanning calorimetry (DSC) showed that the melting temperature (T_m) dropped in XYNB and XYNC from 104.9 °C to 93.7 °C and 78.6 °C, respectively. The detrimental mutations showed that extremely thermostable enzymes can tolerate quite radical mutations in the protein's interior and still retain high thermostability. The analysis of mutations (F806H and F806V) in a hydrophobic area lining the substrate-binding region indicated that active site hydrophobicity is important for high activity at extreme temperatures. Although polar His at 806 provided higher stability, the hydrophobic Phe at 806 provided higher activity than His. This study generates an understanding of how extreme thermostability and high activity are formed in GH10 xylanases. KEY POINTS: • Characterization and molecular dynamics simulations of TmXYN10B and its mutants • Explanation of structural stability of GH10 xylanase.

The encapsulin from Thermotoga maritima is a flavoprotein with a symmetry matched ferritin-like cargo protein

Bacterial nanocompartments, also known as encapsulins, are an emerging class of protein-based 'organelles' found in bacteria and archaea. Encapsulins are virus-like icosahedral particles comprising a ~ 25-50 nm shell surrounding a specific cargo enzyme. Compartmentalization is thought to create a unique chemical environment to facilitate catalysis and isolate toxic intermediates. Many questions regarding nanocompartment structure-function remain unanswered, including how shell symmetry dictates cargo loading and to what extent the shell facilitates enzymatic activity. Here, we explore these questions using the model Thermotoga maritima nanocompartment known to encapsulate a redox-active ferritin-like protein. Biochemical analysis revealed the encapsulin shell to possess a flavin binding site located at the interface between capsomere subunits, suggesting the shell may play a direct and active role in the function of the encapsulated cargo. Furthermore, we used cryo-EM to show that cargo proteins use a form of symmetry-matching to facilitate encapsulation and define stoichiometry. In the case of the Thermotoga maritima encapsulin, the decameric cargo protein with fivefold symmetry preferentially binds to the pentameric-axis of the icosahedral shell. Taken together, these observations suggest the shell is not simply a passive barrier-it also plays a significant role in the structure and function of the cargo enzyme.

Role of a high centrality residue in protein dynamics and thermal stability

Centralities determined from Residue Interaction Networks (RIN) in proteins have been used to predict aspects of their structure and dynamics. Here, we correlate the Eigenvector Centrality (Ec) with the rate constant for thermal denaturation (kden) of the HisF protein from Thermotoga maritima based on 12 single alanine substitution mutants. The molecular basis for this correlation was further explored by studying a mutant containing a replacement of a high Ec residue, Y182A, which displayed increased kden at 80 °C. The crystallographic structure of this mutant showed few changes, mostly in two flexible loops. The 1H-15N -HSQC showed only subtle changes of cross peak positions for residues located near the mutation site and scattered throughout the structure. However, the comparison of the RIN showed that Y182 is the vertex of a set of high centrality residues that spreads throughout the HisF structure, which is lacking in the mutant. Cross-correlation displacements of Cα calculated from a molecular dynamics simulation at different temperatures showed that the Y182A mutation reduced the correlated movements in the HisF structure above 70 °C. 1H-15N NMR chemical shift covariance using temperature as perturbation were consistent with these results. In conclusion the increase in temperature drives the structure of the mutant HisF-Y182A into a less connected state, richer in non-concerted motions, located predominantly in the C-terminal half of the protein where Y182 is placed. Conversely, wild-type HisF responds to increased temperature as a single unit. Hence the replacement of a high Ec residue alters the distribution of thermal energy through HisF structure.

Thermostable D-amino acid decarboxylases derived from Thermotoga maritima diaminopimelate decarboxylase

Diaminopimelate decarboxylases (DAPDCs) are highly selective enzymes that catalyze the common final step in different lysine biosynthetic pathways, i.e. the conversion of meso-diaminopimelate (DAP) to L-lysine. We examined the modification of the substrate specificity of the thermostable decarboxylase from Thermotoga maritima with the aim to introduce activity with 2-aminopimelic acid (2-APA) since its decarboxylation leads to 6-aminocaproic acid (6-ACA), a building block for the synthesis of nylon-6. Structure-based mutagenesis of the distal carboxylate binding site resulted in a set of enzyme variants with new activities toward different D-amino acids. One of the mutants (E315T) had lost most of its activity toward DAP and primarily acted as a 2-APA decarboxylase. We next used computational modeling to explain the observed shift in catalytic activities of the mutants. The results suggest that predictive computational protocols can support the redesign of the catalytic properties of this class of decarboxylating PLP-dependent enzymes.

Biocatalysts from cyanobacterial hapalindole pathway afford antivirulent isonitriles against MRSA

The emergence of resistance to frontline antibiotics has called for novel strategies to combat serious pathogenic infections. Methicillin-resistant Staphylococcus aureus [MRSA] is one such pathogen. As opposed to traditional antibiotics, bacteriostatic anti-virulent agents disarm MRSA, without exerting pressure, that cause resistance. Herein, we employed a thermophilic Thermotoga maritima tryptophan synthase (TmTrpB1) enzyme followed by an isonitrile synthase and Fe(II)-α-ketoglutarate-dependent oxygenase, in sequence as biocatalysts to produce antivirulent indole vinyl isonitriles. We report on conversion of simple derivatives of indoles to their C3-vinyl isonitriles, as the enzymes employed here demonstrated broader substrate tolerance. In toto, eight distinct L-Tryptophan derived α-amino acids (7) were converted to their bioactive vinyl isonitriles (3) by action of an isonitrile synthase (WelI1) and an Fe(II)-α-ketoglutarate-dependent oxygenase (WelI3) yielding structural variants possessing antivirulence against MRSA. These indole vinyl isonitriles at 10 μg/mL are effective as antivirulent compounds against MRSA, as evidenced through analysis of rabbit blood hemolysis assay. Based on a homology modelling exercise, of enzyme-substrate complexes, we deduced potential three dimensional alignments of active sites and glean mechanistic insights into the substrate tolerance of the Fe(II)-α-ketoglutarate-dependent oxygenase.

Scalable continuous evolution for the generation of diverse enzyme variants encompassing promiscuous activities

Enzyme orthologs sharing identical primary functions can have different promiscuous activities. While it is possible to mine this natural diversity to obtain useful biocatalysts, generating comparably rich ortholog diversity is difficult, as it is the product of deep evolutionary processes occurring in a multitude of separate species and populations. Here, we take a first step in recapitulating the depth and scale of natural ortholog evolution on laboratory timescales. Using a continuous directed evolution platform called OrthoRep, we rapidly evolve the Thermotoga maritima tryptophan synthase β-subunit (TmTrpB) through multi-mutation pathways in many independent replicates, selecting only on TmTrpB's primary activity of synthesizing L-tryptophan from indole and L-serine. We find that the resulting sequence-diverse TmTrpB variants span a range of substrate profiles useful in industrial biocatalysis and suggest that the depth and scale of evolution that OrthoRep affords will be generally valuable in enzyme engineering and the evolution of biomolecular functions.

Development of a Protein Scaffold for Arginine Sensing Generated through the Dissection of the Arginine-Binding Protein from Thermotoga maritima

Arginine is one of the most important nutrients of living organisms as it plays a major role in important biological pathways. However, the accumulation of arginine as consequence of metabolic defects causes hyperargininemia, an autosomal recessive disorder. Therefore, the efficient detection of the arginine is a field of relevant biomedical/biotechnological interest. Here, we developed protein variants suitable for arginine sensing by mutating and dissecting the multimeric and multidomain structure of Thermotoga maritima arginine-binding protein (TmArgBP). Indeed, previous studies have shown that TmArgBP domain-swapped structure can be manipulated to generate simplified monomeric and single domain scaffolds. On both these stable scaffolds, to measure tryptophan fluorescence variations associated with the arginine binding, a Phe residue of the ligand binding pocket was mutated to Trp. Upon arginine binding, both mutants displayed a clear variation of the Trp fluorescence. Notably, the single domain scaffold variant exhibited a good affinity (~3 µM) for the ligand. Moreover, the arginine binding to this variant could be easily reverted under very mild conditions. Atomic-level data on the recognition process between the scaffold and the arginine were obtained through the determination of the crystal structure of the adduct. Collectively, present data indicate that TmArgBP scaffolds represent promising candidates for developing arginine biosensors.

Specificity in the biosynthesis of the universal tRNA nucleoside N6-threonylcarbamoyl adenosine (t6A)-TsaD is the gatekeeper

N⁶-threonylcarbamoyl adenosine (t⁶A) is a nucleoside modification found in all kingdoms of life at position 37 of tRNAs decoding ANN codons, which functions in part to restrict translation initiation to AUG and suppress frameshifting at tandem ANN codons. In Bacteria the proteins TsaB, TsaC (or C2), TsaD, and TsaE, comprise the biosynthetic apparatus responsible for t⁶A formation. TsaC(C2) and TsaD harbor the relevant active sites, with TsaC(C2) catalyzing the formation of the intermediate threonylcarbamoyladenosine monophosphate (TC-AMP) from ATP, threonine, and CO₂, and TsaD catalyzing the transfer of the threonylcarbamoyl moiety from TC-AMP to A₃₇ of substrate tRNAs. Several related modified nucleosides, including hydroxynorvalylcarbamoyl adenosine (hn⁶A), have been identified in select organisms, but nothing is known about their biosynthesis. To better understand the mechanism and structural constraints on t⁶A formation, and to determine if related modified nucleosides are formed via parallel biosynthetic pathways or the t⁶A pathway, we carried out biochemical and biophysical investigations of the t⁶A systems from E. coli and T. maritima to address these questions. Using kinetic assays of TsaC(C2), tRNA modification assays, and NMR, our data demonstrate that TsaC(C2) exhibit relaxed substrate specificity, producing a variety of TC-AMP analogs that can differ in both the identity of the amino acid and nucleotide component, whereas TsaD displays more stringent specificity, but efficiently produces hn⁶A in E. coli and T. maritima tRNA. Thus, in organisms that contain modifications such as hn⁶A in their tRNA, we conclude that their origin is due to formation via the t⁶A pathway.

Enhanced Biomass Yield of and Saccharification in Transgenic Tobacco Over-Expressing β-Glucosidase

Here, we report an increase in biomass yield and saccharification in transgenic tobacco plants (Nicotiana tabacumL.) overexpressing thermostable β-glucosidase from Thermotoga maritima, BglB, targeted to the chloroplasts and vacuoles. The transgenic tobacco plants showed phenotypic characteristics that were significantly different from those of the wild-type plants. The biomass yield and life cycle (from germination to flowering and harvest) of the transgenic tobacco plants overexpressing BglB were 52% higher and 36% shorter than those of the wild-type tobacco plants, respectively, indicating a change in the genome transcription levels in the transgenic tobacco plants. Saccharification in biomass samples from the transgenic tobacco plants was 92% higher than that in biomass samples from the wild-type tobacco plants. The transgenic tobacco plants required a total investment (US$/year) corresponding to 52.9% of that required for the wild-type tobacco plants, but the total biomass yield (kg/year) of the transgenic tobacco plants was 43% higher than that of the wild-type tobacco plants. This approach could be applied to other plants to increase biomass yields and overproduce β-glucosidase for lignocellulose conversion.

Non-chromatographic purification of thermostable endoglucanase from Thermotoga maritima by fusion with a hydrophobic elastin-like polypeptide

Endoglucanase EG12B from Thermotoga maritima is a thermophilic cellulase that has great potential for industrial applications. Here, to enable the selective purification of EG12B in a simple and efficient manner, an elastin-like polypeptide (ELP), which acts as a thermally responsive polypeptide, was fused with EG12B to enable its inverse phase transition cycling (ITC). A small gene library comprising ELPs from ELP₅ to ELP₅₀ was constructed using recursive directional ligation by plasmid reconstruction. ELP₅₀ was added to the C-terminus of EG12B as a fusion tag to obtain the expression vector pET28-EG12B-ELP₅₀, which was transformed into Escherichia coli BL21 (DE3) to enable the expression of fusion protein via IPTG induction. Gray scanning analysis revealed that the EG12B-ELP₅₀ expression level was up to about 35% of the total cellular proteins. After three rounds of ITC, 8.14 mg of EG12B-ELP₅₀ was obtained from 500-mL lysogeny broth culture medium. The recovery rate and purification fold of EG12B-ELP₅₀ purified by ITC reached 78.1% and 11.8, respectively. The cellulase activity assay showed that EG12B-ELP₅₀ had a better thermostability, higher optimal temperature, and longer half-life than those of free EG12B. Overall, our results suggested that ELP₅₀ could be used as a favorable fusion tag, providing a rapid, simple, and inexpensive strategy for non-chromatographic target-protein purification.

A thermostable 5-enolpyruvylshikimate-3-phosphate synthase from Thermotoga maritima enhances glyphosate tolerance in Escherichia coli and transgenic Arabidopsis

5-Enolpyruvylshikimate-3-phosphate synthase (EPSPS) overexpression, attempting to provide excess EPSPS to combine with glyphosate, is one way to improve glyphosate resistance of plants. The EPSPS in extremophiles which is selected by nature to withstand the evolutionary pressure may possess some potential-specific biological functions. In this study, we reported the cloning, expression and enzymatic characterization of a novel Class II EPSPS AroAT. maritima from Thermotoga maritima MSB8. The enzyme showed low sequence identities with other EPSPSs, and was one of the most thermostable EPSPSs so far, which showed the optimum enzyme activity at 80 °C. The enzyme maintains the activity below 50 °C and in a wide range of pH 4.0-10, which indicated its stability under rough environment, especially in tropical regions and alkaline soil. Excellent Ki/Km value of AroAT. maritima suggested that the enzyme showed powerful competitive binding capacity of PEP over glyphosate and high glyphosate tolerance. Furthermore, aroAT. maritima gene was transformed into Arabidopsis thaliana. The transgenic lines were resistant to 15 mM glyphosate, which proved the application value in the cultivation of glyphosate-tolerant plants.

Differential Substrate Recognition by Maltose Binding Proteins Influenced by Structure and Dynamics

The genome of the hyperthermophile Thermotoga maritima contains three isoforms of maltose binding protein (MBP) that are high-affinity receptors for di-, tri-, and tetrasaccharides. Two of these proteins (tmMBP1 and tmMBP2) share significant sequence identity, approximately 90%, while the third (tmMBP3) shares less than 40% identity. MBP from Escherichia coli (ecMBP) shares 35% sequence identity with the tmMBPs. This subset of MBP isoforms offers an interesting opportunity to investigate the mechanisms underlying the evolution of substrate specificity and affinity profiles in a genome where redundant MBP genes are present. In this study, the X-ray crystal structures of tmMBP1, tmMBP2, and tmMBP3 are reported in the absence and presence of oligosaccharides. tmMBP1 and tmMBP2 have binding pockets that are larger than that of tmMBP3, enabling them to bind to larger substrates, while tmMBP1 and tmMBP2 also undergo substrate-induced hinge bending motions (∼52°) that are larger than that of tmMBP3 (∼35°). Small-angle X-ray scattering was used to compare protein behavior in solution, and computer simulations provided insights into dynamics of these proteins. Comparing quantitative protein-substrate interactions and dynamical properties of tmMBPs with those of the promiscuous ecMBP and disaccharide selective Thermococcus litoralis MBP provides insights into the features that enable selective binding. Collectively, the results provide insights into how the structure and dynamics of tmMBP homologues enable them to differentiate between a myriad of chemical entities while maintaining their common fold.

Structure and mechanism of a bacterial t6A biosynthesis system

The universal N(6)-threonylcarbamoyladenosine (t6A) modification at position 37 of ANN-decoding tRNAs is central to translational fidelity. In bacteria, t6A biosynthesis is catalyzed by the proteins TsaB, TsaC/TsaC2, TsaD and TsaE. Despite intense research, the molecular mechanisms underlying t6A biosynthesis are poorly understood. Here, we report biochemical and biophysical studies of the t6A biosynthesis system from Thermotoga maritima. Small angle X-ray scattering analysis reveals a symmetric 2:2 stoichiometric complex of TsaB and TsaD (TsaB2D2), as well as 2:2:2 complex (TsaB2D2E2), in which TsaB acts as a dimerization module, similar to the role of Pcc1 in the archaeal system. The TsaB2D2 complex is the minimal platform for the binding of one tRNA molecule, which can then accommodate a single TsaE subunit. Kinetic data demonstrate that TsaB2D2 alone, and a TsaB2D2E1 complex with TsaE mutants deficient in adenosine triphosphatase (ATPase) activity, can catalyze only a single cycle of t6A synthesis, while gel shift experiments provide evidence that the role of TsaE-catalyzed ATP hydrolysis occurs after the release of product tRNA. Based on these results, we propose a model for t6A biosynthesis in bacteria.

S-adenosyl-L-homocysteine hydrolase from a hyperthermophile (Thermotoga maritima) is expressed in Escherichia coli in inactive form - Biochemical and structural studies

Thermotoga maritima is a hyperthermophilic bacterium but its genome encodes a number of archaeal proteins including S-adenosyl-L-homocysteine hydrolase (SAHase), which regulates cellular methylation reactions. The question of proper folding and activity of proteins of extremophilic origin is an intriguing problem. When expressed in E.coli and purified (as a homotetramer) at room temperature, the hyperthermophilic SAHase from T.maritima was inactive. ITC study indicated that the protein undergoes heat-induced conformational changes, and enzymatic activity assays demonstrated that these changes are required to attain enzymatic activity. To explain the mechanism of thermal activation, two crystal structures of the inactive form of T. maritima SAHase (iTmSAHase) were determined for an incomplete binary complex with the reduced cofactor (NADH), and in a mixture of binary complexes with NADH and with adenosine. In contrast to active SAHases, in iTmSAHase only two of the four subunits contain a bound cofactor, predominantly in its non-reactive, reduced state. Moreover, the closed-like conformation of the cofactor-containing subunits precludes substrate delivery to the active site. The two other subunits cannot be involved in the enzymatic reaction either; although they have an open-like conformation, they do not contain the cofactor, whose binding site may be occupied by an adenosine molecule. The results suggest that this enzyme, when expressed in mesophilic cells, is arrested in the activity-incompatible conformation revealed by its crystal structures.

Impact of orientation and flexibility of peptide linkers on T. maritima lipase Tm1350 displayed on Bacillus subtilis spores surface using CotB as fusion partner

Fusion protein construction often requires peptide linkers for prolonged conformation, extended stability and enzyme activity. In this study a series of fusion between Thermotoga maritima lipase Tm1350 and Bacillus subtillis coat protein CotB, comprising of several peptide linkers, with different length, flexibility and orientations were constructed. Effects of temperature, pH and chemicals were examined, on the activity of displayed enzyme. The fusion protein with longer flexible linkers L9 [(GGGGS)₄] and L7 (GGGGS-GGGGS-EAAAK-EAAAK-GGGGS-GGGGS) possess 1.29 and 1.16-fold higher activity than the original, under optimum temperature and pH respectively. Moreover, spore surface displaying Tm1350 with L3 (EAAAK-GGGGS) and L9 ((GGGGS)₄) showed extended thermostably, maintaining 1.40 and 1.35-fold higher activity than the original respectively, at 80 °C after 5 h of incubation. The enzyme activity of linkers with different orientation, including L5, L6 and L7 was determined, where L7 maintained 1.05 and 1.27-fold higher activity than L5 and L6. Effect of 0.1% proteinase K, bromelain, 20% ethanol and 30% methanol was investigated. Linkers with appropriate Glycine residues (flexible) showed higher activity than Alanine residues (rigid). The activity of the displayed enzyme can be improved by maintaining orientation and flexibility of peptide linkers, to evaluate high activity and stability in industrial processes.

Semi-rational engineering of a thermostable aldo-keto reductase from Thermotoga maritima for synthesis of enantiopure ethyl-2-hydroxy-4-phenylbutyrate (EHPB)

A novel aldo-keto reductase Tm1743 characterized from Thermotoga maritima was explored as an effective biocatalyst in chiral alcohol production. Natural Tm1743 catalyzes asymmetric reduction of ethyl 2-oxo-4-phenylbutyrate (EOPB) at high efficiency, but the production of, ethyl (S)-2-hydroxy-4-phenylbutyrate ((S)-EHPB), which is less desirable, is preferred with an enantiomeric excess (ee) value of 76.5%. Thus, altering the enantioselectivity of Tm1743 to obtain the more valuable product (R)-EHPB for angiotensin drug synthesis is highly desired. In this work, we determined the crystal structure of Tm1743 in complex with its cofactor NADP+ at 2.0 Å resolution, and investigated the enantioselectivity of Tm1743 through semi-rational enzyme design. Molecular simulations based on the crystal structure obtained two binding models representing the pro-S and pro-R conformations of EOPB. Saturation mutagenesis studies revealed that Trp21 and Trp86 play important roles in determining the enantioselectivity of Tm1743. The best (R)- and (S)-EHPB preferring Tm1743 mutants, denoted as W21S/W86E and W21L/W118H, were identified; their ee values are 99.4% and 99.6% and the catalytic efficiencies are 0.81 and 0.12 mM-1s-1, respectively. Our work presents an efficient strategy to improve the enantioselectivity of a natural biocatalyst, which will serve as a guide for further exploration of new green catalysts for asymmetric reactions.

Substrate adaptabilities of Thermotogae mannan binding proteins as a function of their evolutionary histories

The Thermotogae possess a large number of ATP-binding cassette (ABC) transporters, including two mannan binding proteins, ManD and CelE (previously called ManE). We show that a gene encoding an ancestor of these was acquired by the Thermotogae from the archaea followed by gene duplication. To address the functional evolution of these proteins as a consequence of their evolutionary histories, we measured the binding affinities of ManD and CelE orthologs from representative Thermotogae. Both proteins bind cellobiose, cellotriose, cellotetraose, β-1,4-mannotriose, and β-1,4-mannotetraose. The CelE orthologs additionally bind β-1,4-mannobiose, laminaribiose, laminaritriose and sophorose while the ManD orthologs additionally only weakly bind β-1,4-mannobiose. The CelE orthologs have higher unfolding temperatures than the ManD orthologs. An examination of codon sites under positive selection revealed that many of these encode residues located near or in the binding site, suggesting that the proteins experienced selective pressures in regions that might have changed their functions. The gene arrangement, phylogeny, binding properties, and putative regulatory networks suggest that the ancestral mannan binding protein was a CelE ortholog which gave rise to the ManD orthologs. This study provides a window on how one class of proteins adapted to new functions and temperatures to fit the physiologies of their new hosts.

Hyperthermostable Thermotoga maritima xylanase XYN10B shows high activity at high temperatures in the presence of biomass-dissolving hydrophilic ionic liquids

The gene of Thermotoga maritima GH10 xylanase (TmXYN10B) was synthesised to study the extreme limits of this hyperthermostable enzyme at high temperatures in the presence of biomass-dissolving hydrophilic ionic liquids (ILs). TmXYN10B expressed from Pichia pastoris showed maximal activity at 100 °C and retained 92 % of maximal activity at 105 °C in a 30-min assay. Although the temperature optimum of activity was lowered by 1-ethyl-3-methylimidazolium acetate ([EMIM]OAc), TmXYN10B retained partial activity in 15-35 % hydrophilic ILs, even at 75-90 °C. TmXYN10B retained over 80 % of its activity at 90 °C in 15 % [EMIM]OAc and 15-25 % 1-ethyl-3-methylimidazolium dimethylphosphate ([EMIM]DMP) during 22-h reactions. [EMIM]OAc may rigidify the enzyme and lower V max. However, only minor changes in kinetic parameter K m showed that competitive inhibition by [EMIM]OAc of TmXYN10B is minimal. In conclusion, when extended enzymatic reactions under extreme conditions are required, TmXYN10B shows extraordinary potential.

Structural Dissection of the Active Site of Thermotoga maritima β-Galactosidase Identifies Key Residues for Transglycosylating Activity

Glycoside hydrolases, specifically β-galactosidases, can be used to synthesize galacto-oligosaccharides (GOS) due to the transglycosylating (secondary) activity of these enzymes. Site-directed mutagenesis of a thermoresistant β-galactosidase from Thermotoga maritima has been carried out to study the structural basis of transgalactosylation and to obtain enzymatic variants with better performance for GOS biosynthesis. Rational design of mutations was based on homologous sequence analysis and structural modeling. Analysis of mutant enzymes indicated that residue W959, or an alternative aromatic residue at this position, is critical for the synthesis of β-3'-galactosyl-lactose, the major GOS obtained with the wild-type enzyme. Mutants W959A and W959C, but not W959F, showed an 80% reduced synthesis of this GOS. Other substitutions, N574S, N574A, and F571L, increased the synthesis of β-3'-galactosyl-lactose about 40%. Double mutants F571L/N574S and F571L/N574A showed an increase of about 2-fold.

Tuning protein mechanics through an ionic cluster graft from an extremophilic protein

Proteins from extremophilic organisms provide excellent model systems to determine the role of non-covalent interactions in defining protein stability and dynamics as well as being attractive targets for the development of robust biomaterials. Hyperthermophilic proteins have a prevalence of salt bridges, relative to their mesophilic homologues, which are thought to be important for enhanced thermal stability. However, the impact of salt bridges on the mechanical properties of proteins is far from understood. Here, a combination of protein engineering, biophysical characterisation, single molecule force spectroscopy (SMFS) and molecular dynamics (MD) simulations directly investigates the role of salt bridges in the mechanical stability of two cold shock proteins; BsCSP from the mesophilic organism Bacillus subtilis and TmCSP from the hyperthermophilic organism Thermotoga maritima. Single molecule force spectroscopy shows that at ambient temperatures TmCSP is mechanically stronger yet, counter-intuitively, its native state can withstand greater deformation before unfolding (i.e. it is mechanically soft) compared with BsCSP. MD simulations were used to identify the location and quantify the population of salt bridges, and reveal that TmCSP contains a larger number of highly occupied salt bridges than BsCSP. To test the hypothesis that salt-bridges endow these mechanical properties on the hyperthermophilic CSP, a charged triple mutant (CTM) variant of BsCSP was generated by grafting an ionic cluster from TmCSP into the BsCSP scaffold. As expected CTM is thermodynamically more stable and mechanically softer than BsCSP. We show that a grafted ionic cluster can increase the mechanical softness of a protein and speculate that it could provide a mechanical recovery mechanism and that it may be a design feature applicable to other proteins.

Systems biology of the structural proteome

BACKGROUND

The success of genome-scale models (GEMs) can be attributed to the high-quality, bottom-up reconstructions of metabolic, protein synthesis, and transcriptional regulatory networks on an organism-specific basis. Such reconstructions are biochemically, genetically, and genomically structured knowledge bases that can be converted into a mathematical format to enable a myriad of computational biological studies. In recent years, genome-scale reconstructions have been extended to include protein structural information, which has opened up new vistas in systems biology research and empowered applications in structural systems biology and systems pharmacology.

RESULTS

Here, we present the generation, application, and dissemination of genome-scale models with protein structures (GEM-PRO) for Escherichia coli and Thermotoga maritima. We show the utility of integrating molecular scale analyses with systems biology approaches by discussing several comparative analyses on the temperature dependence of growth, the distribution of protein fold families, substrate specificity, and characteristic features of whole cell proteomes. Finally, to aid in the grand challenge of big data to knowledge, we provide several explicit tutorials of how protein-related information can be linked to genome-scale models in a public GitHub repository ( https://github.com/SBRG/GEMPro/tree/master/GEMPro_recon/).

CONCLUSIONS

Translating genome-scale, protein-related information to structured data in the format of a GEM provides a direct mapping of gene to gene-product to protein structure to biochemical reaction to network states to phenotypic function. Integration of molecular-level details of individual proteins, such as their physical, chemical, and structural properties, further expands the description of biochemical network-level properties, and can ultimately influence how to model and predict whole cell phenotypes as well as perform comparative systems biology approaches to study differences between organisms. GEM-PRO offers insight into the physical embodiment of an organism's genotype, and its use in this comparative framework enables exploration of adaptive strategies for these organisms, opening the door to many new lines of research. With these provided tools, tutorials, and background, the reader will be in a position to run GEM-PRO for their own purposes.

A Robust and Versatile Method of Combinatorial Chemical Synthesis of Gene Libraries via Hierarchical Assembly of Partially Randomized Modules

A major challenge in gene library generation is to guarantee a large functional size and diversity that significantly increases the chances of selecting different functional protein variants. The use of trinucleotides mixtures for controlled randomization results in superior library diversity and offers the ability to specify the type and distribution of the amino acids at each position. Here we describe the generation of a high diversity gene library using tHisF of the hyperthermophile Thermotoga maritima as a scaffold. Combining various rational criteria with contingency, we targeted 26 selected codons of the thisF gene sequence for randomization at a controlled level. We have developed a novel method of creating full-length gene libraries by combinatorial assembly of smaller sub-libraries. Full-length libraries of high diversity can easily be assembled on demand from smaller and much less diverse sub-libraries, which circumvent the notoriously troublesome long-term archivation and repeated proliferation of high diversity ensembles of phages or plasmids. We developed a generally applicable software tool for sequence analysis of mutated gene sequences that provides efficient assistance for analysis of library diversity. Finally, practical utility of the library was demonstrated in principle by assessment of the conformational stability of library members and isolating protein variants with HisF activity from it. Our approach integrates a number of features of nucleic acids synthetic chemistry, biochemistry and molecular genetics to a coherent, flexible and robust method of combinatorial gene synthesis.

Site-Directed Mutagenesis of a Hyperthermophilic Endoglucanase Cel12B from Thermotoga maritima Based on Rational Design

To meet the demand for the application of high activity and thermostable cellulases in the production of new-generation bioethanol from nongrain-cellulose sources, a hyperthermostable β-1,4-endoglucase Cel12B from Thermotoga maritima was selected for further modification by gene site-directed mutagenesis method in the present study, based on homology modeling and rational design. As a result, two recombinant enzymes showed significant improvement in enzyme activity by 77% and 87%, respectively, higher than the parental enzyme TmCel12B. Furthermore, the two mutants could retain 80% and 90.5% of their initial activity after incubation at 80°C for 8 h, while only 45% for 5 h to TmCel12B. The Km and Vmax of the two recombinant enzymes were 1.97±0.05 mM, 4.23±0.15 μmol·mg(-1)·min(-1) of TmCel12B-E225H-K207G-D37V, and 2.97±0.12 mM, 3.15±0.21 μmol·mg(-1)·min(-1) of TmCel12B-E225H-K207G, respectively, when using CMC-Na as the substrate. The roles of the mutation sites were also analyzed and evaluated in terms of electron density, hydrophobicity of the modeled protein structures. The recombinant enzymes may be used in the hydrolysis of cellulose at higher temperature in the future. It was concluded that the gene mutagenesis approach of a certain active residues may effectively improve the performance of cellulases for the industrial applications and contribute to the study the thermostable mechanism of thermophilic enzymes.

Differential contributions of the latch in Thermotoga maritima reverse gyrase to the binding of single-stranded DNA before and after ATP hydrolysis

Reverse gyrase catalyzes the ATP-dependent introduction of positive supercoils into DNA. Supercoiling requires the functional cooperation of its N-terminal helicase domain with the C-terminal topoisomerase domain. The helicase domain contains a superfamily 2 helicase core formed by two RecA domains, H1 and H2. We show here that a helicase domain lacking the latch, an insertion in H2, fails to close the cleft in the helicase core in response to nucleotide and DNA binding at the beginning of the catalytic cycle. In the presence of the pre-hydrolysis ATP analog ADP·BeFx, however, the closed conformer can still be formed in the absence of the latch. The helicase domain lacking the latch exhibits reduced DNA affinities. The energetic difference between the two nucleotide states involved in duplex separation is diminished, rationalizing the unwinding deficiency of reverse gyrase lacking the latch. The latch most strongly contributes to binding of single-stranded DNA in the post-hydrolysis state, before phosphate release. Our results are in line with contributions of the latch in determining the direction of strand passage, and in orienting the cleaved single-stranded DNA for re-ligation. At the same time, the latch may coordinate the re-ligation reaction with strand passage and with the nucleotide cycle.

Characterization of Glycosynthase Mutants Derived from Glycoside Hydrolase Family 10 Xylanases

Four xylanases belonging to glycoside hydrolase family 10-Thermotoga maritima XylB (TM), Clostridium stercorarium XynB (CS), Bacillus halodurans XynA (BH), and Cellulomonas fimi Cex (CF)-were converted to glycosynthases by substituting the nucleophilic glutamic acid residues with glycine, alanine, and serine. The glycine mutants exhibited the highest levels of glycosynthase activity with all four enzymes. All the glycine mutants formed polymeric beta-1,4-linked xylopyranose as a precipitate during reaction with alpha-xylobiosyl fluoride. Two glycine mutants (TM and CF) recognized X(2) as an effective acceptor molecule to prohibit the formation of the polymer, while the other two (CS and BH) did not. The difference in acceptor specificity is considered to reflect the difference in substrate affinity at their +2 subsites. The results agreed with the structural predictions of the subsite, where TM and CF exhibit high affinity at subsite 2, suggesting that the glycosynthase technique is useful for investigating the affinity of +subsites.

Thermodynamic properties of the β-glucosidase from Thermotoga maritima extend the upper limit of thermophilicity

Enzymes from thermophilic organisms are believed to be strong candidates for industrial applications due to their ability to withstand temperature-induced enzyme inactivation. The present study demonstrated molecular cloning, over-expression, purification and characterization of β-glucosidase from Thermotoga maritima. The bglA gene with a capacity to encode a 51 kDa enzyme was heterologously expressed in E. coli M15. The enzyme was produced @130 mgL(-1) in LB media and @440 mgL(-1) in Dubos salt medium accounting 40-47 % of total cellular soluble proteins when lactose was used as an inducer. The enzyme showed a peak activity between pH and temperature range of 5.0-7.0 and 80-100 °C, respectively. The activity was fairly stable up to 140 °C. The turnover rate (kcat) of the enzyme was 187.1±20 s(-1), whereas the Km and Vmax values were 0.56 mM and 238±2.4 IU mg(-1) protein, respectively. The enzyme was shown to have half-life of 136, 71 and 12.6 h at 80, 90 and 100 °C, respectively. Thermodynamics parameters including melting temperature (130 °C), activation energy for inactivation (36.92 kJmole(-1)), enthalpy (33.73 kJmole(-1)), Gibb's free energy (127.96 kJmole(-1)) and entropy (-246.46 Jmole(-1)K(-1)) have shown that the enzyme have enhanced hydrophobic interactions to prevent its thermal unfolding. These features endorse the industrial applications of the enzyme.

Improved activity of a thermophilic cellulase, Cel5A, from Thermotoga maritima on ionic liquid pretreated switchgrass

Ionic liquid pretreatment of biomass has been shown to greatly reduce the recalcitrance of lignocellulosic biomass, resulting in improved sugar yields after enzymatic saccharification. However, even under these improved saccharification conditions the cost of enzymes still represents a significant proportion of the total cost of producing sugars and ultimately fuels from lignocellulosic biomass. Much of the high cost of enzymes is due to the low catalytic efficiency and stability of lignocellulolytic enzymes, especially cellulases, under conditions that include high temperatures and the presence of residual pretreatment chemicals, such as acids, organic solvents, bases, or ionic liquids. Improving the efficiency of the saccharification process on ionic liquid pretreated biomass will facilitate reduced enzyme loading and cost. Thermophilic cellulases have been shown to be stable and active in ionic liquids but their activity is typically at lower levels. Cel5A_Tma, a thermophilic endoglucanase from Thermotoga maritima, is highly active on cellulosic substrates and is stable in ionic liquid environments. Here, our motivation was to engineer mutants of Cel5A_Tma with higher activity on 1-ethyl-3-methylimidazolium acetate ([C₂mim][OAc]) pretreated biomass. We developed a robotic platform to screen a random mutagenesis library of Cel5A_Tma. Twelve mutants with 25–42% improvement in specific activity on carboxymethyl cellulose and up to 30% improvement on ionic-liquid pretreated switchgrass were successfully isolated and characterized from a library of twenty thousand variants. Interestingly, most of the mutations in the improved variants are located distally to the active site on the protein surface and are not directly involved with substrate binding.

Improved recombinant cellulase expression in chloroplast of tobacco through promoter engineering and 5' amplification promoting sequence

Economical production of bioethanol from lignocellulosic biomass still faces many technical limitations. Cost-effective production of fermentable sugars is still not practical for large-scale production of bioethanol due to high costs of lignocellulolytic enzymes. Therefore, plant molecular farming, where plants are used as bioreactors, was developed for the mass production of cell wall degrading enzymes that will help reduce costs. Subcellular targeting is also potentially more suitable for the accumulation of recombinant cellulases. Herein, we generated transgenic tobacco plants (Nicotiana tabacum cv. SR1) that accumulated Thermotoga maritima BglB cellulase, which was driven by the alfalfa RbcsK-1A promoter and contained a small subunit of the rubisco complex transit peptide. The generated transformants possessed high specific BglB activity and did not show any abnormal phenotypes. Furthermore, we genetically engineered the RbcsK-1A promoter (MRbcsK-1A) and fused the amplification promoting sequence (aps) to MRbcsK-1A promoter to obtain high expression of BglB in transgenic plants. AMRsB plant lines with aps-MRbcsK-1A promoter showed the highest specific activity of BglB, and the accumulated BglB protein represented up to 9.3 % of total soluble protein. When BglB was expressed in Arabidopsis and tobacco plants, the maximal production capacity of recombinant BglB was 0.59 and 1.42 mg/g wet weight, respectively. These results suggests that suitable recombinant expression of cellulases in subcellular compartments such as chloroplasts will contribute to the cost-effective production of enzymes, and will serve as the solid foundation for the future commercialization of bioethanol production via plant molecular farming.

[Residues affecting hydrolysis of soy isoflavone glycosides, stability and catalytic properties of Thermotoga maritime β-glucosidase]

The thermostable β-glucosidase A (TmBglA) from Thermotoga maritime is a promising biocatalyst for production of isoflavone aglycones. Use of enzymes with high specificity for soy isoflavone conjugates is however essential for efficient hydrolysis. The effect of the amino acids located in the aglycone binding pocket with non-conserved residues between specificity groups in family 1 glycoside hydrolase (GH1) was studied using wild-type TmBglA and 3 exchange mutants (MI-TmBglA, M2-TmBglA, M1M2-TmBglA). Three mutants were expressed in Escherichia coli, purified and characterized. They had shifts in both optimum tem- perature and thermal stability, and their narrowing pH-activity curve caused by removing the ionized side chain in mutation. All mutants demonstrated the decreased catalytic efficiency more effectively revealed with natural glycoside, salicin, than with artificial substrate, p-nitrophenyl-β-D-glucopyranoside, suggesting that' these amino acids are the key residues to determine aglycone specificity. A lower hydrolysis of genistin and daidzin for M2-TmBglA than M1-TmBglA indicated that L400, A407 and E408 being preferable to V170, A171, V173, G 174 and H180 residues of Tm-BglA could be essential for soy isoflavone glycoside binding and catalysis.

Homology modeling of the CheW coupling protein of the chemotaxis signaling complex

Homology models of the E. coli and T. maritima chemotaxis protein CheW were constructed to assess the quality of structural predictions and their applicability in chemotaxis research: i) a model of E. coli CheW was constructed using the T. maritima CheW NMR structure as a template, and ii) a model of T. maritima CheW was constructed using the E. coli CheW NMR structure as a template. The conformational space accessible to the homology models and to the NMR structures was investigated using molecular dynamics and Monte Carlo simulations. The results show that even though static homology models of CheW may be partially structurally different from their corresponding experimentally determined structures, the conformational space they can access through their dynamic variations can be similar, for specific regions of the protein, to that of the experimental NMR structures. When CheW homology models are allowed to explore their local accessible conformational space, modeling can provide a rational path to predicting CheW interactions with the MCP and CheA proteins of the chemotaxis complex. Homology models of CheW (and potentially, of other chemotaxis proteins) should be seen as snapshots of an otherwise larger ensemble of accessible conformational space.

Characterization of an invertase with pH tolerance and truncation of its N-terminal to shift optimum activity toward neutral pH

Most invertases identified to date have optimal activity at acidic pH, and are intolerant to neutral or alkaline environments. Here, an acid invertase named uninv2 is described. Uninv2 contained 586 amino acids, with a 100 amino acids N-terminal domain, a catalytic domain and a C-terminal domain. With sucrose as the substrate, uninv2 activity was optimal at pH 4.5 and at 45°C. Removal of N-terminal domain of uninv2 has shifted the optimum pH to 6.0 while retaining its optimum temperaure at 45°C. Both uninv2 and the truncated enzyme retained highly stable at neutral pH at 37°C, and they were stable at their optimum pH at 4°C for as long as 30 days. These characteristics make them far superior to invertase from Saccharomyces cerevisiae, which is mostly used as industrial enzyme.

Discovery and characterization of a novel ATP/polyphosphate xylulokinase from a hyperthermophilic bacterium Thermotoga maritima

Xylulokinase (XK, E.C. 2.7.1.17) is one of the key enzymes in xylose metabolism and it is essential for the activation of pentoses for the sustainable production of biocommodities from biomass sugars. The open reading frame (TM0116) from the hyperthermophilic bacterium Thermotoga maritima MSB8 encoding a putative xylulokinase were cloned and expressed in Escherichia coli BL21 Star (DE3) in the Luria-Bertani and auto-inducing high-cell-density media. The basic biochemical properties of this thermophilic XK were characterized. This XK has the optimal temperature of 85 °C. Under a suboptimal condition of 60 °C, the k cat was 83 s⁻¹, and the K(m) values for xylulose and ATP were 1.24 and 0.71 mM, respectively. We hypothesized that this XK could work on polyphosphate possibly because this ancestral thermophilic microorganism utilizes polyphosphate to regulate the Embden-Meyerhof pathway and its substrate-binding residues are somewhat similar to those of other ATP/polyphosphate-dependent kinases. This XK was found to work on low-cost polyphosphate, exhibiting 41 % of its specific activity on ATP. This first ATP/polyphosphate XK could have a great potential for xylose utilization in thermophilic ethanol-producing microorganisms and cell-free biosystems for low-cost biomanufacturing without the use of ATP.

Structural and functional characterization of a noncanonical nucleoside triphosphate pyrophosphatase from Thermotoga maritima

The hyperthermophilic bacterium Thermotoga maritima has a noncanonical nucleoside triphosphatase that catalyzes the conversion of inosine triphosphate (ITP), deoxyinosine triphosphate (dITP) and xanthosine triphosphate (XTP) into inosine monophosphate (IMP), deoxyinosine monophosphate (IMP) and xanthosine monophosphate (XMP), respectively. The k(cat)/K(m) values determined at 323 and 353 K fall between 1.31 × 10(4) and 7.80 × 10(4) M(-1) s(-1). ITP and dITP are slightly preferred over XTP. Activity towards canonical nucleoside triphosphates (ATP and GTP) was not detected. The enzyme has an absolute requirement for Mg(2+) as a cofactor and has a preference for alkaline conditions. A protein X-ray structure of the enzyme with bound IMP was obtained at 2.15 Å resolution. The active site houses a well conserved network of residues that are critical for substrate recognition and catalysis. The crystal structure shows a tetramer with two possible dimer interfaces. One of these interfaces strongly resembles the dimer interface that is found in the structures of other noncanonical nucleoside pyrophosphatases from human (human ITPase) and archaea (Mj0226 and PhNTPase).

Genetic and biochemical dissection of a HisKA domain identifies residues required exclusively for kinase and phosphatase activities

Two-component signal transduction systems, composed of histidine kinases (HK) and response regulators (RR), allow bacteria to respond to diverse environmental stimuli. The HK can control both phosphorylation and subsequent dephosphorylation of its cognate RR. The majority of HKs utilize the HisKA subfamily of dimerization and histidine phosphotransfer (DHp) domains, which contain the phospho-accepting histidine and directly contact the RR. Extensive genetics, biochemistry, and structural biology on several prototypical TCS systems including NtrB-NtrC and EnvZ-OmpR have provided a solid basis for understanding the function of HK–RR signaling. Recently, work on NarX, a HisKA_3 subfamily protein, indicated that two residues in the highly conserved region of the DHp domain are responsible for phosphatase activity. In this study we have carried out both genetic and biochemical analyses on Myxococcus xanthus CrdS, a member of the HisKA subfamily of bacterial HKs. CrdS is required for the regulation of spore formation in response to environmental stress. Following alanine-scanning mutagenesis of the α1 helix of the DHp domain of CrdS, we determined the role for each mutant protein for both kinase and phosphatase activity. Our results indicate that the conserved acidic residue (E372) immediately adjacent to the site of autophosphorylation (H371) is specifically required for kinase activity but not for phosphatase activity. Conversely, we found that the conserved Thr/Asn residue (N375) was required for phosphatase activity but not for kinase activity. We extended our biochemical analyses to two CrdS homologs from M. xanthus, HK1190 and HK4262, as well as Thermotoga maritima HK853. The results were similar for each HisKA family protein where the conserved acidic residue is required for kinase activity while the conserved Thr/Asn residue is required for phosphatase activity. These data are consistent with conserved mechanisms for kinase and phosphatase activities in the broadly occurring HisKA family of sensor kinases in bacteria.

Bacterial histidine kinases (HK) serve as bifunctional enzymes capable of both phosphorylation and dephosphorylation of their cognate response regulators (RR). The majority of HKs (77%) belong to the HisKA subfamily. While both kinase and phosphatase functions have been assayed for HisKA proteins, relatively few examples have been studied to determine which residues are required for kinase and phosphatase activity. Recent studies on NarX, a HisKA_3 family protein, and the dedicated phosphatases CheZ and CheX illustrate requirements for two amino acids for phosphatase function. In this study, we undertook saturating mutagenesis of the proposed interaction surface between the HK and its cognate RR and conclude that only one residue (T/N) is required exclusively for phosphatase activity for HisKA family proteins in evolutionarily distant organisms Myxococcus xanthus and Thermotoga maritima. In addition, we identified only one residue (E/D), adjacent to the conserved site of phosphorylation, required exclusively for kinase activity within the highly conserved motif H-E/D-x-x-T/N. Because similar sequences are found in nearly all HisKA kinases, these residues provide excellent targets for dissection of kinase and phosphatase activities within this broadly occurring family of bacterial kinases.

Genes for the major structural components of Thermotogales species' togas revealed by proteomic and evolutionary analyses of OmpA and OmpB homologs

The unifying structural characteristic of members of the bacterial order Thermotogales is their toga, an unusual cell envelope that includes a loose-fitting sheath around each cell. Only two toga-associated structural proteins have been purified and characterized in Thermotoga maritima: the anchor protein OmpA1 (or Ompα) and the porin OmpB (or Ompβ). The gene encoding OmpA1 (ompA1) was cloned and sequenced and later assigned to TM0477 in the genome sequence, but because no peptide sequence was available for OmpB, its gene (ompB) was not annotated. We identified six porin candidates in the genome sequence of T. maritima. Of these candidates, only one, encoded by TM0476, has all the characteristics reported for OmpB and characteristics expected of a porin including predominant β-sheet structure, a carboxy terminus porin anchoring motif, and a porin-specific amino acid composition. We highly enriched a toga fraction of cells for OmpB by sucrose gradient centrifugation and hydroxyapatite chromatography and analyzed it by LC/MS/MS. We found that the only porin candidate that it contained was the TM0476 product. This cell fraction also had β-sheet character as determined by circular dichroism, consistent with its enrichment for OmpB. We conclude that TM0476 encodes OmpB. A phylogenetic analysis of OmpB found orthologs encoded in syntenic locations in the genomes of all but two Thermotogales species. Those without orthologs have putative isofunctional genes in their place. Phylogenetic analyses of OmpA1 revealed that each species of the Thermotogales has one or two OmpA homologs. T. maritima has two OmpA homologs, encoded by ompA1 (TM0477) and ompA2 (TM1729), both of which were found in the toga protein-enriched cell extracts. These annotations of the genes encoding toga structural proteins will guide future examinations of the structure and function of this unusual lineage-defining cell sheath.

ATPase cycle and DNA unwinding kinetics of RecG helicase

The superfamily 2 bacterial helicase, RecG, is a monomeric enzyme with a role in DNA repair by reversing stalled replication forks. The helicase must act specifically and rapidly to prevent replication fork collapse. We have shown that RecG binds tightly and rapidly to four-strand oligonucleotide junctions, which mimic a stalled replication fork. The helicase unwinds such DNA junctions with a step-size of approximately four bases per ATP hydrolyzed. To gain an insight into this mechanism, we used fluorescent stopped-flow and quenched-flow to measure individual steps within the ATPase cycle of RecG, when bound to a DNA junction. The fluorescent ATP analogue, mantATP, was used throughout to determine the rate limiting steps, effects due to DNA and the main states in the cycle. Measurements, when possible, were also performed with unlabeled ATP to confirm the mechanism. The data show that the chemical step of hydrolysis is the rate limiting step in the cycle and that this step is greatly accelerated by bound DNA. The ADP release rate is similar to the cleavage rate, so that bound ATP and ADP would be the main states during the ATP cycle. Evidence is provided that the main structural rearrangements, which bring about DNA unwinding, are linked to these states.

Characterization of multi-functional properties and conformational analysis of MutS2 from Thermotoga maritima MSB8

The MutS2 homologues have received attention because of their unusual activities that differ from those of MutS. In this work, we report on the functional characteristics and conformational diversities of Thermotoga maritima MutS2 (TmMutS2). Various biochemical features of the protein were demonstrated via diverse techniques such as scanning probe microscopy (SPM), ATPase assays, analytical ultracentrifugation, DNA binding assays, size chromatography, and limited proteolytic analysis. Dimeric TmMutS2 showed the temperature-dependent ATPase activity. The non-specific nicking endonuclease activities of TmMutS2 were inactivated in the presence of nonhydrolytic ATP (ADPnP) and enhanced by the addition of TmMutL. In addition, TmMutS2 suppressed the TmRecA-mediated DNA strand exchange reaction in a TmMutL-dependent manner. We also demonstrated that small-angle X-ray scattering (SAXS) analysis of dimeric TmMutS2 exhibited nucleotide- and DNA-dependent conformational transitions. Particularly, TmMutS2-ADPnP showed the most compressed form rather than apo-TmMutS2 and the TmMutS2-ADP complex, in accordance with the results of biochemical assays. In the case of the DNA-binding complexes, the stretched conformation appeared in the TmMutS2-four-way junction (FWJ)-DNA complex. Convergences of biochemical- and SAXS analysis provided abundant information for TmMutS2 and clarified ambiguous experimental results.

The structural basis of mode of activation and functional diversity: a case study with HtrA family of serine proteases

HtrA (High temperature requirement protease A) proteins that are primarily involved in protein quality control belong to a family of serine proteases conserved from bacteria to humans. HtrAs are oligomeric proteins that share a common trimeric pyramidal architecture where each monomer comprises a serine protease domain and one or two PDZ domains. Although the overall structural integrity is well maintained and they exhibit similar mechanism of activation, subtle conformational changes and structural plasticity especially in the flexible loop regions and domain interfaces lead to differences in their active site conformation and hence in their specificity and functions.

RNase P: at last, the key finds its lock

Apart from the ribosome, the crystal structure of the bacterial RNase P in complex with a tRNA, reported by Reiter and colleagues recently, constitutes the first example of a multiple turnover RNA enzyme. Except in rare exceptions, RNase P is ubiquitous and, like the ribosome, is older than the initial branch point of the phylogenetic tree. Importantly, the structure shows how the RNA and the protein moieties cooperate to process the pre-tRNA substrates. The catalytic site comprises some critical RNA residues spread over the secondary structure but gathered in a compact volume next to the protein, which helps recognize and orient the substrate. The discussion here outlines some important aspects of that crystal structure, some of which could apply to RNA molecules in general.

In planta differential targeting analysis of Thermotoga maritima Cel5A and CBM6-engineered Cel5A for autohydrolysis

The heterologous expression of glycosyl hydrolases in bioenergy crops can improve the lignocellulosic conversion process for ethanol production. We attempted to obtain high-level expression of an intact Thermotoga maritima endoglucanase, Cel5A, and CBM6-engineered Cel5A in transgenic tobacco plants for the mass production and autohydrolysis of endoglucanase. Cel5A expression was targeted to different subcellular compartments, namely, the cytosol, apoplast, and chloroplast, using the native form of the pathogenesis-related protein 1a (PR1a) and Rubisco activase (RA) transit peptides. Cel5A transgenic tobacco plants with the chloroplast transit peptide showed the highest average endoglucanase activity and protein accumulation up to 4.5% total soluble protein. Cel5A-CBM6 was targeted to the chloroplast and accumulated up to 5.2% total soluble protein. In terms of the direct conversion of plant tissue into free sugar, the Cel5A-CBM6 transgenic plant was 33% more efficient than the Cel5A transgenic plant. The protein stability of Cel5A and Cel5A-CBM6 in lyophilized leaf material is an additional advantage in the bioconversion process.

Probing for DNA damage with β-hairpins: similarities in incision efficiencies of bulky DNA adducts by prokaryotic and human nucleotide excision repair systems in vitro

Nucleotide excision repair (NER) is an important prokaryotic and eukaryotic defense mechanism that removes a large variety of structurally distinct lesions in cellular DNA. While the proteins involved are completely different, the mode of action of these two repair systems is similar, involving a cut-and-patch mechanism in which an oligonucleotide sequence containing the lesion is excised. The prokaryotic and eukaryotic NER damage-recognition factors have common structural features of β-hairpin intrusion between the two DNA strands at the site of the lesion. In the present study, we explored the hypothesis that this common β-hairpin intrusion motif is mirrored in parallel NER incision efficiencies in the two systems. We have utilized human HeLa cell extracts and the prokaryotic UvrABC proteins to determine their relative NER incision efficiencies. We report here comparisons of relative NER efficiencies with a set of stereoisomeric DNA lesions derived from metabolites of benzo[a]pyrene and equine estrogens in different sequence contexts, utilizing 21 samples. We found a general qualitative trend toward similar relative NER incision efficiencies for ∼65% of these substrates; the other cases deviate mostly by ∼30% or less from a perfect correlation, although several more distant outliers are also evident. This resemblance is consistent with the hypothesis that lesion recognition through β-hairpin insertion, a common feature of the two systems, is facilitated by local thermodynamic destabilization induced by the lesions in both cases. In the case of the UvrABC system, varying the nature of the UvrC endonuclease, while maintaining the same UvrA/B proteins, can markedly affect the relative incision efficiencies. These observations suggest that, in addition to recognition involving the initial modified duplexes, downstream events involving UvrC can also play a role in distinguishing and processing different lesions in prokaryotic NER.

Production of hyperthermostable GH10 xylanase Xyl10B from Thermotoga maritima in transplastomic plants enables complete hydrolysis of methylglucuronoxylan to fermentable sugars for biofuel production

Overcoming the recalcitrance in lignocellulosic biomass for efficient hydrolysis of the polysaccharides cellulose and hemicellulose to fermentable sugars is a research priority for the transition from a fossilfuel-based economy to a renewable carbohydrate economy. Methylglucuronoxylans (MeGXn) are the major components of hemicellulose in woody biofuel crops. Here, we describe efficient production of the GH10 xylanase Xyl10B from Thermotoga maritima in transplastomic plants and demonstrate exceptional stability and catalytic activities of the in planta produced enzyme. Fully expanded leaves from homotransplastomic plants contained enzymatically active Xyl10B at a level of 11-15% of their total soluble protein. Transplastomic plants and their seed progeny were morphologically indistinguishable from non-transgenic plants. Catalytic activity of in planta produced Xyl10B was detected with poplar, sweetgum and birchwood xylan substrates following incubation between 40 and 90 °C and was also stable in dry and stored leaves. Optimal yields of Xyl10B were obtained from dry leaves if crude protein extraction was performed at 85 °C. The transplastomic plant derived Xyl10B showed exceptional catalytic activity and enabled the complete hydrolysis of MeGXn to fermentable sugars with the help of a single accessory enzyme (α-glucuronidase) as revealed by the sugar release assay. Even without this accessory enzyme, the majority of MeGXn was hydrolyzed by the transplastomic plant-derived Xyl10B to fermentable xylose and xylobiose.

Co2+ selectivity of Thermotoga maritima CorA and its inability to regulate Mg2+ homeostasis present a new class of CorA proteins

CorA is a family of divalent cation transporters ubiquitously present in bacteria and archaea. Although CorA can transport both Mg(2+) and Co(2+) almost equally well, its main role has been suggested to be that of primary Mg(2+) transporter of prokaryotes and hence the regulator of Mg(2+) homeostasis. The reason is that the affinity of CorA for Co(2+) is relatively low and thus considered non-physiological. Here, we show that Thermotoga maritima CorA (TmCorA) is incapable of regulating the Mg(2+) homeostasis and therefore cannot be the primary Mg(2+) transporter of T. maritima. Further, our in vivo experiments confirm that TmCorA is a highly selective Co(2+) transporter, as it selects Co(2+) over Mg(2+) at >100 times lower concentrations. In addition, we present data that show TmCorA to be extremely thermostable in the presence of Co(2+). Mg(2+) could not stabilize the protein to the same extent, even at high concentrations. We also show that addition of Co(2+), but not Mg(2+), specifically induces structural changes to the protein. Altogether, these data show that TmCorA has the role of being the transporter of Co(2+) but not Mg(2+). The physiological relevance and requirements of Co(2+) in T. maritima is discussed and highlighted. We suggest that CorA may have different roles in different organisms. Such functional diversity is presumably a reflection of minor, but important structural differences within the CorA family that regulate the gating, substrate selection, and transport.

The genus Thermotoga: recent developments

The genus Thermotoga comprises extremely thermophilic (Topt > or = 70 degrees C) and hyperthermophilic (Topt > or = 80 degrees C) bacteria, which have been extensively studied for insights into the basis for life at elevated temperatures and for biotechnological opportunities (e.g. biohydrogen production, biocatalysis). Over the past decade, genome sequences have become available for a number of Thermotoga species, leading to functional genomics efforts to understand growth physiology as well as genomics-based identification and characterization of novel high-temperature biocatalysts. Discussed here are recent developments along these lines for this group of microorganisms.

Plant cell calcium-rich environment enhances thermostability of recombinantly produced alpha-amylase from the hyperthermophilic bacterium Thermotoga maritime

In the industrial processing of starch for sugar syrup and ethanol production, a liquefaction step is involved where starch is initially solubilized at high temperature and partially hydrolyzed with a thermostable and thermoactive alpha-amylase. Most amylases require calcium as a cofactor for their activity and stability, therefore calcium, along with the thermostable enzyme, are typically added to the starch mixture during enzymatic liquefaction, thereby increasing process costs. An attractive alternative would be to produce the enzyme directly in the tissue to be treated. In a proof of concept study, tobacco cell cultures were used as model system to test in planta production of a hyperthermophilic alpha-amylase from Thermotoga maritima. While comparable biochemical properties to recombinant production in Escherichia coli were observed, thermostability of the plant-produced alpha-amylase benefited significantly from high intrinsic calcium levels in the tobacco cells. The plant-made enzyme retained 85% of its initial activity after 3 h incubation at 100 degrees C, whereas the E. coli-produced enzyme was completely inactivated after 30 min under the same conditions. The addition of Ca(2+) or plant cell extracts from tobacco and sweetpotato to the E. coli-produced enzyme resulted in a similar stabilization, demonstrating the importance of a calcium-rich environment for thermostability, as well as the advantage of producing this enzyme directly in plant cells where calcium is readily available.

Structure of a putative beta-phosphoglucomutase (TM1254) from Thermotoga maritima

The structure of TM1254, a putative beta-phosphoglucomutase from T. maritima, was determined to 1.74 A resolution in a high-throughput structural genomics programme. Diffraction data were obtained from crystals belonging to space group P22(1)2(1), with unit-cell parameters a = 48.16, b = 66.70, c = 83.80 A, and were refined to an R factor of 19.2%. The asymmetric unit contained one protein molecule which is comprised of two domains. Structural homologues were found from protein databases that confirmed a strong resemblance between TM1254 and members of the haloacid dehalogenase (HAD) hydrolase family.

Structural analysis of semi-specific oligosaccharide recognition by a cellulose-binding protein of thermotoga maritima reveals adaptations for functional diversification of the oligopeptide periplasmic binding protein fold

Periplasmic binding proteins (PBPs) constitute a protein superfamily that binds a wide variety of ligands. In prokaryotes, PBPs function as receptors for ATP-binding cassette or tripartite ATP-independent transporters and chemotaxis systems. In many instances, PBPs bind their cognate ligands with exquisite specificity, distinguishing, for example, between sugar epimers or structurally similar anions. By contrast, oligopeptide-binding proteins bind their ligands through interactions with the peptide backbone but do not distinguish between different side chains. The extremophile Thermotoga maritima possesses a remarkable array of carbohydrate-processing metabolic systems, including the hydrolysis of cellulosic polymers. Here, we present the crystal structure of a T. maritima cellobiose-binding protein (tm0031) that is homologous to oligopeptide-binding proteins. T. maritima cellobiose-binding protein binds a variety of lengths of beta(1-->4)-linked glucose oligomers, ranging from two rings (cellobiose) to five (cellopentaose). The structure reveals that binding is semi-specific. The disaccharide at the nonreducing end binds specifically; the other rings are located in a large solvent-filled groove, where the reducing end makes several contacts with the protein, thereby imposing an upper limit of the oligosaccharides that are recognized. Semi-specific recognition, in which a molecular class rather than individual species is selected, provides an efficient solution for the uptake of complex mixtures.

Hydrolysis of soy isoflavone glycosides by recombinant beta-glucosidase from hyperthermophile Thermotoga maritima

A recombinant Thermotoga maritima beta-glucosidase A (BglA) was purified to homogeneity for performing enzymatic hydrolysis of isoflavone glycosides from soy flour. The kinetic properties K(m), k(cat), and k(cat)/K(m) of BglA towards isoflavone glycosides, determined using high-performance liquid chromatography, confirmed the higher efficiency of BglA in hydrolyzing malonylglycosides than non-conjugated glycosides (daidzin and genistin). During hydrolysis of soy flour by BglA at 80 degrees C, the isoflavone glycosides (soluble form) were extracted from soy flour (solid state) into the solution (liquid state) in thermal condition and converted to their aglycones (insoluble form), which mostly existed in the pellet to be separated from BglA in the reaction solution. The enzymatic hydrolysis in one-step and two-step approaches yielded 0.38 and 0.35 mg genistein and daidzein per gram of soy flour, respectively. The optimum conditions for conversion of isoflavone aglycones were 100 U per gram of soy flour, substrate concentration 25% (w/v), and incubation time 3 h for 80 degrees C.

Three-dimensional structural view of the central metabolic network of Thermotoga maritima

Metabolic pathways have traditionally been described in terms of biochemical reactions and metabolites. With the use of structural genomics and systems biology, we generated a three-dimensional reconstruction of the central metabolic network of the bacterium Thermotoga maritima. The network encompassed 478 proteins, of which 120 were determined by experiment and 358 were modeled. Structural analysis revealed that proteins forming the network are dominated by a small number (only 182) of basic shapes (folds) performing diverse but mostly related functions. Most of these folds are already present in the essential core (approximately 30%) of the network, and its expansion by nonessential proteins is achieved with relatively few additional folds. Thus, integration of structural data with networks analysis generates insight into the function, mechanism, and evolution of biological networks.