Abstract
Two new thermophilic branched chain amino acid transaminases have been identified within the genomes of different hyper-thermophilic archaea, Geoglobus acetivorans, and Archaeoglobus fulgidus. These enzymes belong to the class IV of transaminases as defined by their structural fold. The enzymes have been cloned and over-expressed in Escherichia coli and the recombinant enzymes have been characterized both biochemically and structurally. Both enzymes showed high thermostability with optimal temperature for activity at 80 and 85°C, respectively. They retain good activity after exposure to 50% of the organic solvents, ethanol, methanol, DMSO and acetonitrile. The enzymes show a low activity to (R)-methylbenzylamine but no activity to (S)-methylbenzylamine. Both enzymes have been crystallized and their structures solved in the internal aldimine form, to 1.9 Å resolution for the Geoglobus enzyme and 2.0 Å for the Archaeoglobus enzyme. Also the Geoglobus enzyme structure has been determined in complex with the amino acceptor α-ketoglutarate and the Archaeoglobus enzyme in complex with the inhibitor gabaculine. These two complexes have helped to determine the conformation of the enzymes during enzymatic turnover and have increased understanding of their substrate specificity. A comparison has been made with another (R) selective class IV transaminase from the fungus Nectria haematococca which was previously studied in complex with gabaculine. The subtle structural differences between these enzymes has provided insight regarding their different substrate specificities.
Introduction
Transaminases (TAms; aminotransferases; EC2.6.1.-) catalyse a transfer of an amino group between a donor substrate and an acceptor molecule (Mehta et al., ). They are involved in the metabolism of most natural amino acids and collectively have a broad substrate range. This makes them very useful for the biocatalytic synthesis of chiral amino acid analogs (Leuchtenberger et al., ; Tufvesson et al., 2011) such as tert-leucine (Taylor et al., 1998) and homophenylalanine (Cho et al., ).
The TAms belong to a large group of enzymes that utilize pyridoxal 5′-phosphate (PLP) (Hayashi, ; Jansonius, ). PLP is the biologically active form of vitamin B6 (Braunstein and Shemyakin, ; Metzler et al., ). The mechanism of these enzymes involves the PLP cofactor forming a covalent bond (a Schiff base) with the NZ atom of the active-site lysine (internal aldimine). In the first half of the TAm reaction, a Schiff base is formed between the donor substrate and PLP (external aldimine). After several reaction stages, the substrate amino group is transferred to the cofactor to form pyridoxamine 5′-phosphate (PMP) with the release of the ketone product. In the second half-reaction the PMP amino group is transferred to an acceptor ketone or aldehyde to form the PLP internal aldimine and a primary amine or a chiral amino product. Most TAms use an α-ketoglutarate (AKG) acceptor substrate to form a glutamate product. Some TAms, however utilize the pyruvate acceptor to form alanine and are usually referred to as substrate:pyruvate TAms (Ward and Wohlgemuth, 2010; Sayer et al., 2012).
The TAms are currently widely used in industrial biocatalysis for both the resolution and the asymmetric synthesis of chiral amines which are used as building blocks for drug intermediates. There is specific interest in TAms capable of the production of both (S) and (R)-chiral intermediates of non α-amino acid substrates. These enzymes catalyse amination of substrates with a ketone/aldehyde component distal to the carboxyl group [ω-TAms; (Malik et al., )] or with no carboxyl group [amine TAms; (Shin et al., 2003)]. Several Pfam class III (S)-selective ωTAms (Punta et al., 2012) from Vibrio fluvialis (Shin et al., 2003), Chromobacterium violaceum and Pseudomonas aeruginosa (Kaulmann et al., ), Arthrobacter sp. and Bacillus megaterium (van Oosterwijk et al., 2016) have been biochemically and structurally characterized (Humble et al., ; Midelfort et al., ; Sayer et al., 2013). They show activity toward the (S)-isoform of α-methylbenzylamine [(S)-MBA] and also toward other amines and ω-amino acids such as β-alanine. The amino acceptor can be pyruvate and a variety of other ketones and aldehydes (Kaulmann et al., ).
Most α-TAms and (S)-ωTAms, are specific toward the (S)-enantiomer of the substrate and have a PLP type I fold. Enzymes of the Pfam TAm class IV have a PLP type IV fold and catalyse transamination of branched chain L-amino acids (BCAAs), D-amino acids or ω-amino acids. The catalysis in class IV TAms occurs on the re-side of the cofactor, which differs from most of the other PLP enzymes where the reactions occur on the si-side. Since both (S) and (R)-enantiomers can be turned over within the same PLP type IV scaffold, enzymes capable of transamination of (R)-ω-substituted amines are found within this class of TAms. By comparing the sequence for class IV enzymes which are potentially capable of (R)-ωTAm activity a sequence fingerprint has been established (Höhne et al., ). Biochemical characterization of a number of these enzymes has demonstrated a high enantioselectivity toward the substrates (R)-α-methylbenzylamine [(R)-MBA], (2R)-amino-4-phenylbutane and (2R)-aminohexane.
The production of the (R)-ω-substituted chiral amines is a challenging process using conventional chemical synthesis. Therefore, enzymes capable of catalyzing the selective transamination of such compounds have received significant interest. An example of such an enzyme is a mutant (R)-specific transaminase from Arthrobacter sp. (ArRMut11), created in a collaboration between Codexis and Merck, which was able to catalyse specific amination of the sterically demanding 1,3-ketoamides to generate the (R)-chiral amine for the manufacture of Sitagliptin, a treatment for type II diabetes (Savile et al., 2010). This mutant ArRMut11 enzyme has been reported in reactions for amination of bicyclic ketones including tetralone (Mutti et al., ; Pressnitz et al., 2013). The crystal structures have been reported for several class IV TAm enzymes including the fungal (R)-amine:pyruvate ωTAms from Aspergillus fumigatus (Thomsen et al., 2014), Aspergillus terreus (Łyskowski et al., ) and Nectria haematococca (Sayer et al., 2014), and the bacterial (R)-selective ωTAm ArRMut11 from Arthrobacter sp (Guan et al., ).
There is an increasing demand for enzymes which are more robust to the demanding conditions used in industry. Enzymes found in thermophilic organisms have increased thermostability and are more tolerant to organic solvents and proteolytic cleavage. Solvent stability is advantageous since non-natural substrates used industrially often require the addition of organic solvents to the reaction mixture for substrate solubilisation (Littlechild et al., ). Also the biocatalytic process can be carried out at elevated temperatures where many non-natural substrates have improved solubility when using a thermostable enzyme which can be reused through several reaction cycles. This reduces the overall cost of the enzyme in the industrial process which is often a limitation in the development of a biocatalytic process.
Branched chain TAms (BCATs) catalyse reversible transamination of branched chain amino acids (shown in Scheme 1). Recently, archaeal thermophilic BCATs have been biochemically characterized from Thermococcus sp (Uchida et al., 2014) and biochemically and structurally studied from the thermophile Thermoproteus uzoniensis (Boyko et al., ).
Scheme 1
This paper reports the identification, biochemical and structural characterization of two new thermostable archaeal class IV TAms from Geoglobus acetivorans (Querellou et al., 2009) and Archaeoglobus fulgidus (Stetter, 1988). Both of these hyper-thermophilic archaea have been isolated from different deep sea hydrothermal vents and they share 79.2% sequence identity. The structures of the G. acetivorans enzyme have been determined in the internal aldimine form and in complex with the amino acceptor AKG and the structures of the A. fulgidus enzyme in the internal aldimine form and in complex with the inhibitor gabaculine. The different structural complexes of these related enzymes have given further insight into the overall mechanism of BCATs and their high stability for industrial application and their substrate specificity.
Results and Discussion
Enzyme Cloning, Expression, and Purification
The genes encoding two putative BCATs were identified in the genomes of G. acetivorans (Mardanov et al., ) and A. fulgidus (Klenk et al., ). Both proteins called GEO1900 and AF0933 have been cloned and over-expressed in a soluble form in Escherichia coli and have been purified to homogeneity using metal affinity and size exclusion chromatography. The recombinant BCATs GEO1900 (MW of subunit 32.6 kDa, 292 amino acids) and AF0933 (MW of subunit 32.4 kDa, 290 amino acids) are closely related with a sequence identity of 79.2% and 94.8% similarity. When purified by high resolution gel filtration chromatography the native molecular weight of the two enzymes varied with the GEO1900 approximately 70 kDa, indicating that the enzyme was a homodimer with small amount of tetramer in solution (Figure S1). However, the AF0933 enzyme had a native molecular mass of approximately 220 kDa as determined by size exclusion chromatography, indicating that it forms a homo-hexamer with only small amounts of a homo-dimeric enzyme observed (Figure S2). Both proteins showed absorption at 420 nm indicating that the cofactor PLP was bound in the aldimine form (data not shown). Both of the GEO1900 and AF0933 proteins run as a single band of 36 kDa on SDS-PAGE which is in line with the calculated molecular mass of the His-tagged proteins (35.4 kDa).
Biochemical Characterization of GEO1900 and AF0933
Both enzymes, GEO1900 and AF0933, showed significant activity toward keto analogs of BCAAs (Tables 1, 2). The apparent kinetic parameters for the AF0933 catalyzed transamination reaction between L-glutamate and the branched chain oxoacids were similar with a preference for 3-methyl-2-oxopentanoate. The extrapolated Vmax and Km values for L-glutamate were significantly higher (Table 2 and Figure S7). These values of Vmax and Km point to the low affinity binding of L-Glu in the active site. This was observed earlier for other BCATs from P. aeruginosa and Gluconobacter oxydans (Norton and Sokatch, ; Tachiki and Tochikura, 1973; Kanda et al., ). For GEO1900 the extrapolated Vmax and Km values for pyruvate and AKG were (1.9 ± 0.1 U/mg and 10.6 ± 1.9 mM) and (18.2 ±1.2 U/mg and 1.9 ± 0.6 mM), respectively. The calculated catalytic efficiency constant for GEO1900 toward AKG was found to be higher than toward pyruvate, 9.8 s−1 mM−1 compared to 0.1 s−1 mM−1. Both values fall within the range of the calculated values for the canonical BCATs (Bezsudnova et al., ).
Table 1
| Amino acceptor | Specific activity (U/mg) | |
|---|---|---|
| GEO1900a | AF0933b | |
| 3-methyl-2-oxobutyrate | n.a. | 5.14 |
| 4-methyl-2-oxopentanoate | 7.7 | 3.05 |
| 3-methyl-2-oxopentanoate | 7.1 | 3.23 |
| 2-oxopentanoate | n.a. | 6.51 |
| 2-oxohexanoate | n.a. | 4.26 |
| 2-oxooctanoate | n.a. | 2.4 |
| 2-oxobutyrate | 1.3 | n.a. |
| pyruvate | 0.56 | 0.73 |
| phenylglyoxylate | n.a. | 4.55 |
| β-phenylpyruvate | n.a. | 2.84 |
| indole-3-pyruvate | n.a. | 1.02 |
The specific activity of GEO1900 and AF0933 toward different keto acids.
5 mM L-glutamate was used as the amino donor and 5 mM keto acid was used at 65°C in 50 mM phosphate buffer, pH 8.0; supplemented with 60 μM PLP.
10 mM L-glutamate was used as the amino donor and 10 mM keto acid was used at 70°C in 100 mM Tris-HCl buffer pH 7.5 supplemented with 0.1 mM PLP.
Not all donor acceptor pairs were tested for both enzymes.
Table 2
| Substrate | Vmax (U/mg) | Km (mM) | kcat/Km (s−1 mM−1) |
|---|---|---|---|
| 3-methyl-2-oxobutyratea | 6.3 | 0.42 | 8.85 |
| 4-methyl-2-oxopentanoatea | 5.5 | 0.25 | 12.98 |
| 3-methyl-2-oxopentanoatea | 3.9 | 0.13 | 17.7 |
| L-glutamateb | 64.6* | 121 | 0.32 |
The apparent kinetic parameters of the transamination reaction catalyzed by AF0933.
At 10 mM L-glutamate.
At 5 mM 3-methyl-2-oxobutyrate.
The Vmax for L-glutamate is an extrapolated value calculated at saturated concentration of co-substrate.
As for the activity toward amino donors the profile of activity toward amino acids for GEO1900 (Table S1) was similar to that reported for related BCATs from E.coli, Bacillus brevis and Pseudomonas aeruginosa (Norton and Sokatch, ; Lee-Peng et al., ; Inoue et al., ; Kanda et al., ; Yu et al., 2014). GEO1900 was highly active toward branched-chain amino acids and their isomers. The AF0933 catalyzed the transamination between L-valine and AKG with a specific activity of 5.76 U/mg and between L-alanine and 3-methyl-2-oxobutanoic acid with a lower specific activity of 0.16 U/mg. When the enantio-preference of the GEO1900 and AF0933 in the reaction with the (R) and (S)-isomers of methylbenzylamine (MBA) was investigated using the acetophenone assay they both showed a small but significant activity toward (R)-MBA at pH 7.5, with a value of 2 and 10 mU/mg with 2-oxobutyrate as the keto substrate, respectively. Neither BCAT showed any activity toward (S)-MBA.
The temperature optimum for the transamination reaction catalyzed by GEO1900 and AF0933 was found to be 80 and 85°C, respectively (Figures S3, S5). The pH optimum for amination of keto acids with L-Glu was between pH 7.0 and 8.0 for both enzymes. The thermostability of GEO1900 and AF0933 was evaluated using the analysis of the residual activity after pre-incubation of the enzyme for 2 h at different temperatures. Significant loss of activity of GEO1900 was observed at temperatures higher than 50°C. The pre-incubation of the enzyme at 70°C resulted in 25% loss of activity after 2 h (Figure S4). At the same time more than 60% of the AF0933 residual activity remained after pre-incubation at 75°C for 2 h (Figure S6). The pre-incubation of AF0933 at 85°C for 60 min resulted in a drop of the residual activity to 30%.
The solvent stability of GEO1900 and AF0933 was tested in ethanol, methanol, DMSO and acetonitrile. After 1 h of incubation in up to 50% of these solvents GEO1900 retained close to 100% activity at both 20°C and 55°C. For AF0933 after 1 h incubation of the enzyme in up to 50% of the tested solvents at the similar temperatures 80% of the original activity remained.
Structural Studies
The GEO1900 holoform (internal aldimine form; GEO1900_holo) and the unproductive complex of holoenzyme with the keto substrate AKG-(GEO1900_AKG) have been determined to 1.9 and 2.2 Å resolution, respectively, in the trigonal space group P3221. There are three GEO1900 subunits per asymmetric unit, which form a tight hexamer with 32 point group symmetry with their symmetry mates. The AF0933 structure has been determined in the holoform (AF0933_holo) and in complex with the inhibitor gabaculine (GABC) (AF0933_GABC) to 2.1 and 2.0 Å resolution, respectively, in the orthorhombic P212121 space group. The asymmetric unit contains a hexameric AF0933 molecule which is very similar to GEO1900 hexamer observed in the crystal. All of the structures have been refined to low R-factors with good stereochemical parameters as shown in Table 3. In the GEO1900 structures two N-terminal amino acids which are not present in AF0933, were not defined in the electron density. The interdomain loop residues 121–129 could be fully built into the subunits A and B in GEO1900_AKG structure. The bound ligand AKG was found only in subunit A and was clearly defined in the electron density (Figure 1). The interdomain loop could only be fully built in one subunit (D) out of the six subunits in both the AF0933_holo and AF0933_GABC structures. The crystal contacts in AF0933 appear to be contributing to the lower disorder of the loop in the D subunit. The inhibitor, gabaculine was not present in the electron density in the subunit F of AF0933, therefore this subunit was built as an internal aldimine in the AF0933_GABC structure. In other subunits the gabaculine-PLP adduct was modeled with partial occupancy in the range of 0.6–0.7. There are no residues in a cis conformation in the well-defined regions of the GEO1900 and AF0933 structures. Residue Glu249 is a Ramachandran plot outlier (Ramakrishnan and Ramachandran, 1965) in all subunits of both AF0933 structures. The corresponding Glu251 has similar main chain angles in the GEO1900 structures, but lies in a generously allowed region of the Ramachandran plot (Laskowski et al., ). The fungal class IV TAms from Aspergillus and Nectria species (Łyskowski et al., ; Sayer et al., 2014; Thomsen et al., 2014) have a Gly residue at this position. Many residues in both of the GEO1900 and AF0933 structures have been modeled with multiple conformations of some amino acid side chains. The residues Gln107 and Asn108 were modeled with alternative conformations of the main chain in most subunits of the AF0933_holo structure and in the subunits E and F of the AF0933_GABC structure.
Table 3
| Holo enzyme GEO1900 | α-ketoglutarate complex GEO1900 | Holo enzyme AF0933 | Gabaculine complex AF0933 | |
|---|---|---|---|---|
| Diffraction source | BL41XU, Spring8 | BL41XU, Spring8 | I04-1, Diamond | I04-1, Diamond |
| Wavelength (Å) | 1.0 | 1.0 | 0.92 | 0.92 |
| Space group | P3221 | P3221 | P212121 | P212121 |
| a, b, c (Å) | 117.65, 117.65, 135.98 | 117.31, 117.31, 135.31 | 72.92, 139.12, 167.57 | 73.21, 140.23, 168.32 |
| α, β, γ (°) | 90, 90, 120 | 90, 90, 120 | 90, 90, 90 | 90, 90, 90 |
| Resolution range (Å) | 58.83–1.90 (2.00–1.90)a | 67.66–2.20 (2.26–2.20) | 69.60–2.10 (2.14–2.10) | 67.17–1.98 (2.01–1.98) |
| Number of unique reflections | 85638 (11949) | 55018 (4457) | 99985 (4906) | 119864 (5675) |
| Completeness (%) | 99.6 (98.6) | 99.9 (99.8) | 99.9 (99.6) | 99.0 (95.6) |
| Average redundancy | 19.8 (19.3) | 17.7 (18.4) | 6.3 (5.4) | 5.6 (4.7) |
| <I/σ(I)> | 31.75 (4.0) | 15.1 (4.9) | 10.0 (0.8) | 11.2 (0.8) |
| Rmeas (%)b | 6.8 (79.2) | 13.3 (67.4) | 14.7 (251.1) | 12.8 (233.5) |
| CC1/2 (%; Diederichs and Karplus, ) | 100 (92.9) | 99.8 (92.6) | 99.7 (29.2) | 99.7 (24.9) |
| Overall B factor from Wilson plot (Å2)c | 39.8 | 42.3 | 51.3 | 44.8 |
| Rfact (%) | 15.3 | 15.0 | 18.3 | 18.1 |
| Rfree (%) | 18.8 | 20.8 | 21.8 | 20.9 |
| Refined protein atoms | 6671 | 6836 | 14432 | 14440 |
| Refined solvent atoms | 626 | 508 | 754 | 828 |
| Refined ligand atoms | 0 | 10 | 0 | 125 |
| AVERAGE B FACTOR (Å2) | ||||
| Protein | 37.9 | 38.2 | 45.8 | 42.7 |
| Solvent | 45.3 | 40.9 | 50.6 | 47.9 |
| R.m.s.d. bond lengths (Å) | 0.021 | 0.019 | 0.013 | 0.013 |
| R.m.s.d. bond angles (°) | 1.99 | 2.11 | 1.78 | 1.69 |
| RAMACHANDRAN PLOT ANALYSIS, RESIDUES IN (%)d | ||||
| Most favored regions | 90.5 | 89.9 | 90.6 | 90.7 |
| Additional allowed regions | 8.7 | 8.8 | 8.6 | 8.5 |
| Generously allowed regions | 0.8 | 1.1 | 0.4 | 0.4 |
| Disallowed regions | 0 | 0.3 | 0.4 | 0.4 |
Summary of data processing and refinement statistics.
Values for the highest resolution shell are given in parentheses.
meas = Σh[m/(m – 1)]1/2 Σi|Ih, i| – < Ih>/ Σh ΣiIh, i (Karplus and Diederichs, ).
Wilson B-factor was estimated by SFCHECK (Vaguine et al., 1999).
Ramachandran plot analysis was performed by PROCHECK (Laskowski et al., ).
Figure 1
Subunit Structure
The subunit of the archaeal BCATs has a PLP type IV fold typical of other BCATs, D-amino acid TAms and (R)-selective ωTAms (Mehta et al., ; Höhne et al., ) and consists of two α/β domains, a small domain and a large domain connected by two inter-domain loops (Peisach et al., 1998; Okada et al., ; Goto et al., ; Hirotsu et al., ; Boyko et al., ) (Figure 2).
Figure 2
The small domain (residues 1–119 and 287–292 in GEO1900) is composed of three α-helices wrapped around a seven stranded antiparallel β-sheet with topology −1, +4x, +1, −2x, −1, +4x as defined by PROMOTIF (Richardson, 1981; Kabsch and Sander,
Quaternary Structure of BCATs
A dimer formation is required for the BCAT catalytic activity, since residues from both subunits in the dimer contribute to the active site (Figure 3). However, both GEO1900 and AF0933 form hexamers (trimers of dimers) in the crystal (Figure 4), with the GEO1900 hexamer located on a crystallographic dyad. Dimerization of GEO1900 results in 2400 Å2 or 18% of each subunit solvent accessible surface area being buried (Krissinel and Henrick,
Figure 3

A cartoon representation of the dimer of the AF0933 BCAT gabaculine complex. The gabaculine-PLP covalent adduct molecules are shown as spheres.
Figure 4

A cartoon representation of the AF0933 BCAT holo structure in the native hexameric form. The PLP cofactor molecules are shown as spheres.
Active Site
The active site of the archaeal BCATs is located at the bottom of the cleft formed by residues of both domains of one subunit as well as of the small domain of an adjacent subunit within the catalytic dimer. The same arrangement has also been observed in other BCATs (Okada et al.,
In the holo and complex structures of the archaeal BCATs the PLP moiety is clearly visible in the electron density map and has B factors similar to those of the neighboring protein residues, indicating a full occupancy of the cofactor. The PLP forms a number of contacts with the protein molecule as found in other BCATs (Hutson,
Ligand Binding and Enzyme Negative Cooperativity
A bound AKG molecule is clearly visible in the electron density (Figure 1) of one subunit of the GEO1900_AKG complex and has a B-factor comparable to the average for the protein. A contact analysis reveals that α-carboxylic as well as α-keto groups of the ligand are tightly bound in the active site via hydrogen bonds with the highly conserved residues Tyr89, Ala249, Thr248, Arg91 as well as Ala250 and Gly188 (via a conserved water molecule D238) (Figure S8). Although most of the residues around the γ-carboxylic group of the ligand are hydrophobic, this group is also connected to the protein via several hydrogen bonds with Arg91, Tyr29′ and Leu103′ (via the water molecule–D923).
Comparative analysis of the holo form of GEO1900 and its complex with AKG reveals that the binding of the acceptor substrate does not result in any significant conformational changes in the active site, except for the rotation of the side chains of Leu101′ and Leu103′, accompanied by a small main chain displacement (distance between corresponding Cα-atoms is about 0.6 Å). In addition the AKG binding residues 120–124 of the inter-domain loop shift toward the active site. This positions the side chains of residues Trp120 and Leu123 so that they shield the ligand from the solvent (Figures S9, S10). Interestingly, in the case of the two archaeal BCATs the residues Leu123 and Tyr124 are not part of the consensus sequence X-G-X-Y-L (where the first X is an aromatic residue) which is typical of other BCATs, rather this motif is X-G-X-L-Y in both GEO1900 and AF0933. This permutation (Tyr to Leu and vice versa) leads to a change in the shape of the active site cavity, making the A-pocket more extended in the direction of Leu123.
The inhibitor gabaculine makes an irreversible complex with the cofactor PLP due to the migration of the double bond of the external aldimine onto the gabaculine ring (Figure 5). However, the PLP-gabaculine adduct in the AF0933 complex structure has only partial occupancy of the gabaculine in the enzyme active site. Occupancy of the ligand in subunit F is below the level for an easy build into the model structure, therefore the internal aldimine was modeled into this subunit. The active site of the AF0933 BCAT is much more restrictive for binding of the gabaculine moiety than that of the fungal (R)-selective ωTAm from Nectria haematococca previously studied by us where it was found at 100% occupancy (Sayer et al., 2014). However, when 100:1 excess of gabaculine to PLP was added to the crystallization droplet it was expected to find full occupancy of the PLP-ligand adduct, unless AF0933 BCAT exhibits some negative cooperativity for substrate binding. Binding of the gabaculine inhibitor to one BCAT subunit appears to be sensed by neighboring subunits, with conformational adjustments preventing the binding of the inhibitor in the other subunits. Negative cooperativity has been reported earlier for other PLP enzymes (Stetefeld et al., 2006; Ruzicka and Frey, 2010) and is in agreement with the fact that the AKG binding in the GEO1900 complex structure in only found in one of the three subunits.
Figure 5

A view from outside onto the active site cavity entrance using a surface representation of the AF0933 BCAT complex with the inhibitor gabaculine (shown as spheres).
Inter-domain Loop Conformation
The conformation of the loop 121–129 is significantly different between the BCAT enzymes and other class IV enzymes such as (R)-selective Nectria TAm where this loop covers the active site, creating an enclosed environment which would facilitate the activation of the amine nitrogen. In the AF0933 BCAT enzyme this loop is poorly ordered and moves away from the active site. In other BCATs including the GEO1900 BCAT this loop is closer to the active site but does not cover it as tightly as seen in the Nectria TAm. However, the loop appears to become more ordered in the GEO1900_AKG complex.
Comparison of Structures of the Native Enzymes and the Complexes of GEO1900 and AF0933
From biochemical analysis it is known that AKG does not inhibit the activity of GEO1900 in the concentration range tested (Figure S11). It would appear that the crystal contacts play a role in increasing the affinity of GEO1900 toward AKG in one single subunit of the enzyme. Superposition of the subunits from the asymmetric unit gave a RMSD between the Cα atoms of about 0.2 and 0.5 Å in case of GEO1900_holo and GEO1900_AKG, respectively. In the case of the GEO1900_holo structure the conformations of the flexible loops have small differences between the different subunits. However, these are likely to be due to crystal packing effects. The structure of the complex GEO1900_AKG demonstrates that a major difference between the subunits can be assigned to a movement of the inter-domain loop which has a sufficiently different conformation in subunit A (as shown in Figure S9), where the ligand was bound, when compared to the other subunits. This loop was reported to shield the active site from the solvent upon the substrate binding (Hirotsu et al.,
Superposition of the GEO1900_holo subunits on the GEO1900_AKG subunits reveals a further major difference in the orientation of the inter-domain loop. This loop is usually disordered in the absence of a bound substrate (Okada et al.,
Although the AF0933 BCAT GABC complex was obtained by co-crystallization with gabaculine, differences in the main chain atom positions between the holoenzyme structure and the complex are minimal. No movement of the interdomain loop upon formation of the PLP-gabaculine complex was observed. The carboxyl group of the gabaculine-PLP adduct points away from both distal and proximal carboxyl group binding pockets of the BCAT. This is probably due to the orientation of this adduct being sterically hindered by the shape of the AF0933 BCAT active site cleft.
Comparison to Structures of Other BCATs and Other CLASS IV Transaminases
The BLAST search found several BCATs with known three-dimensional structure which have a moderate primary sequence similarity to GEO1900 and AF0933. The structure-based sequence alignment (Figure 6), shows the comparison of the GEO1900 and AF0933 sequences with other BCATs and another class IV TAm enzyme from the fungal species Nectria. This shows that the Nectria enzyme has an extended N-terminus and also a difference in the loop residues 99–108 (AF0933 numbering) between β5 and β6 of the N-terminal small domain. This loop contributes to the active site and is much larger in the GEO1900 and AF0933 BCATs (Figure 7).
Figure 6

Multiple sequence alignment of different BCATs, A. flugidus, G. acetivorans, T. thermophilus, E.coli, T. uzoniensis, N. haematococca. Arrows indicate β-strands, and helical curves denote α-helices of the structure of AF0933 above and Nectria TAm below. The active site lysine is highlighted in red. The figure was prepared with ESPript3 (Robert and Gouet, 2014).
Figure 7

An overlay of the Cα trace of the AF0933 BCAT (green and coral coil) gabaculine complex with the Nectria TAm (ice blue and gold coil) gabaculine complex (Sayer et al., 2014) to illustrate the differences between the ligand and loop conformations within the active site between the two enzymes. The PLP gabaculine adduct is shown as a stick model for the AF0933 BCAT structure and as thick lines for the Nectria amine TAm structure. Interdomain loops and the loop connecting β5 and β6 of the N-terminal small domain of the adjacent subunit in the catalytic dimer covering the active site cavity, have different conformation between the AF0933 BCAT and the Nectria TAm as shown.
The superposition of the GEO1900_holo subunit and the subunits of the AF0933 and other homologous enzymes using PDBeFOLD showed that the fold of the new thermostable BCATs is very similar to related enzymes. Small differences are observed within the flexible loops that connect the secondary structure elements. However, the inter-domain loop conformation is similar between GEO1900_AKG and the complexes of the E.coli BCAT with substrates bound (PDB codes−1I1M and 1I1L).
A comparison of the structures of GEO1900 and AF0933 has been made with other class IV TAms. The overall structures of these enzymes are similar for DATAs, BCATs and (R)-selective amine:pyruvate TAms. It appears that the differences are in the amino acid residues lining the active site cavity that determine substrate specificities. The BCAT catalyzed transamination between branched chain L-amino acids and the AKG acceptor produces the branched chain ketoacid and L-glutamate. The DATAs catalyze the exchange of the amino group between D-amino acids and the AKG. The (R)-selective amine:pyruvate TAms catalyze transfer of an amino group between a (R)-ω-substrate and pyruvate. Höhne et al. (
Structural Basis for GEO1900 and AF0933 Thermostability
Both BCATs GEO1900 and AF0933 demonstrate high thermal and solvent stability. The closest sequence homologs of these enzymes with a known 3-D structure are from another thermostable archaeal BCAT from T. uzoniensis (45% sequence identity to AF0933; Boyko et al.,
Features that are known to confer increased thermostability of proteins include an increase in ionic interactions or salt bridges which offer more thermostability when found in clusters, increased hydrophobicity and the shortening of surface loops (Littlechild et al.,
The AF0933 BCAT has much shorter loops connecting secondary structure elements in comparison with the E. coli BCAT. Particularly shortened are the loops connecting strands β3 to β4 and β5 to β6 of the N-terminal domain and β5 toβ6 of the large domain. The C-terminal inter-domain loop, which is only 2 residues long in the AF0933 BCAT, is 9 residues long in the E. coli enzyme. Shorter loops also mean that there is a higher percentage of residues belonging to the secondary structure, 55.9% in AF0933 BCAT compared to 49.2% in the E. coli BCAT. The inter-subunit interface in the catalytic dimer is mainly hydrophobic both in the AF0933 and in the E. coli BCAT, therefore it appears that hydrophobic interactions contribute equally to the stability of both proteins.
Thus, an increased number of ion pairs per subunit alongside shorter surface loops are the main contributing factors to the increased thermostability of AF0933 when compared to the E.coli BCAT.
Application to Biocatalysis
Both GEO1900 and AF0933 make a good thermostable and solvent tolerant scaffold which can be optimized for production of natural and synthetic chiral amines of interest to the pharmaceutical industries. The new BCAT enzymes have a chiral preference for the (R)-isomer of MBA but with low activity. However, they offer the possibility to be further optimized for this and other substrates of interest using the knowledge of the active site cavity obtained from the current structural studies. The active site cavities of GEO1900 and AF0933 can be extended to better accommodate MBA and other bulky substrates. The comparison of GEO1900 and AF0933 with a (R)-specific TAm from a Nectria sp. shows that this fungal enzyme has an additional long N-terminal α-helix that increases the depth of the active site of this enzyme. The presence of this N-terminal helix shown in Figure 8, a much tighter closure over the active site of the interdomain loop, and a different conformation of the loop between β5 and β6 of the N-terminal small domain, creates a difference in the substrate environment of the Nectria TAm in comparison with that in the GEO1900 and AF0933 BCATs. This would potentially alter the local dielectric constant within the fungal enzyme active site thereby facilitating the observed amination reactions of ω-keto acids and different aldehydes and ketones by this TAm (Sayer et al., 2014).
Figure 8

A cartoon representation of the superimposition of the catalytic dimers of the AF0933 gabaculine complex with the Nectria TAm gabaculine complex, illustrating the additional α-helix at the N-terminal region of Nectria TAm as shown in Figure 6. The PLP gabaculine adduct is shown in ball and stick mode.
Conclusions
The approach developed in the THERMOGENE project was to search for these new archaeal TAm enzymes within the untapped resource of “natures” biodiversity in order to find enzymes that have evolved to be active under the extreme growth conditions of their host organisms. This approach has been successful and has resulted in the discovery of new thermostable BCATs displaying new features which would have been missed in a directed evolution approach using already identified mesophilic enzymes. The thermophilic archaea G. acetivorans and A. fulgidus were isolated from different hydrothermal sources. The G. acetivorans has been isolated from samples of deep-sea hydrothermal vents of the Ashadze field, the deepest of the known World Ocean hydrothermal fields (Querellou et al., 2009) whereas the A. fulgidus was isolated in shallow hyperthermal sea water vents in Volcano Island in Sicily (Stetter, 1988) yet they have related BCAT enzymes with 79% sequence identity and similar properties. The GEO1900 and AF0933 enzymes have been characterized both biochemically and structurally. They both demonstrate significant thermostability and organic solvent stability and are promising candidates for industrial applications. Both the AF0933 and the GEO1900 are very thermostable as determined by the stability trials conducted. The AF0933 is a hexamer in solution as measured by gel filtration studies whereas the GEO1900 is mainly found as a dimer. The crystallographic studies have provided structures of both enzymes in the internal aldimine form, to 1.9 Å resolution for the Geoglobus enzyme and 2.0 Å for the Archaeoglobus enzyme. The structures of the enzymes with substrate and inhibitor have revealed the conformational changes occurring in the enzymes during catalytic turnover. A comparison of the Archaeoglobus inhibitor complex structure with the related (R)-selective fungal Nectria enzyme structure has produced further insights into differences in substrate specificity between the two enzymes. The new robust (R) selective BCAT enzymes will add to the enzyme “tool box” available for industrial biocatalysis and will have applications for production of new important chiral intermediates for the pharmaceutical industries using a sustainable route.
Materials and Methods
All reagents were purchased from Sigma–Aldrich, Buchs, Switzerland unless otherwise stated. The chromatography columns were obtained from GE Healthcare (Little Chalfont, UK).
Expression and Purification
One gene encoding a class IV transaminase was identified in the complete genome sequence of G. acetivorans (Mardanov et al.,
Another class IV transaminase gene was identified by searching known thermophilic archaeal genomes using the BCAT from N. haematococca as a query sequence. This identified a 290 amino acid protein in A. fulgidus (Klenk et al.,
Enzyme Activity Assays
The overall transamination reaction of GEO1900 was assayed by a discontinuous method by measuring the rate of formation of L-glutamic acid from AKG in the presence of BCAAs in the direct reaction or the formation of AKG from L-glutamate acid in the presence of keto acids in the reverse reaction. In the direct reaction assay GEO1900 (5–30 μg) was pre-incubated in 1 ml of the reaction mixture, containing 50 mM phosphate buffer, pH 8.0, 100 mM NaCl, 60 μM PLP, and 5 mM of the amino acid at 65°C for 3 min. The reaction was initiated by the addition of 5 mM AKG. Samples (100–300 μl) were taken at different time points and frozen to stop the reaction. Then the L-glutamic acid concentration was evaluated in the samples after dilution to 1 ml by spectrophotometry employing the glutamate dehydrogenase (GluDH) assay at 25°C in 50 mM Tris-HCl buffer, pH 9.0, supplemented with 1 mM NAD and 0.5 U GluDH. In the reverse reaction assay the enzyme (5–30 μg) was preincubated in the reaction mixture (1 ml) containing 50 mM phosphate buffer, pH 8.0, 100 mM NaCl, 60 μM PLP, and 5 mM L-glutamic acid at 65°C for 3 min. The reaction was initiated by the addition of 5 mM ketoacid. Samples (100–300 μl) were taken at several time points and frozen to stop the reaction. The α-ketoglutarate concentration was evaluated in the samples after dilution to 1 ml by spectrophotometry employing the GluDH assay at 25°C in 20 mM sodium phosphate buffer, pH 7.3, supplemented with 0.09 mM NADH, 30 mM ammonium chloride and 0.5 U GluDH. The specific activity of GEO1900 was calculated from the initial linear region of the progress curve of the L-glutamic acid/ α-ketoglutarate formation.
Steady-state kinetic parameters for pyruvate and AKG were determined from the substrate saturation curves for the overall reaction between the keto substrate and 5 mM L-Leu as a co-substrate. The saturation curves were analyzed using the Michaelis-Menten model. The GEO1900 activity toward (R) and (S)-MBA was measured spectrophotometrically at 65°C by the acetophenone assay (Schätzle et al., 2011; Boyko et al.,
The AF0933 BCAT activity was measured as 3-methyl-2-oxobutyrate dependent formation of AKG from L-glutamate. One unit of enzyme activity corresponded to the conversion of 1 μmol of substrate consumed or product formed per min. The standard assay mixture contained 0.1 M Tris-HCl, pH 8, 0.1 mM PLP, 5 mM 3-methyl-2-oxobutyrate, 10 mM L-glutamate and protein. Samples were incubated at 70°C up to 40 min, the reaction was stopped on ice and the amount of AKG formed was quantified by measuring NADH consumption in the reductive amination of AKG to L-glutamate at 340 nm. The quantification mixture (0.3 ml) was incubated at 20°C for 10 min and contained 0.1 M Tris-HCl, pH 7.5, 0.3 mM NADH, 40 mM NH4Cl, 5–100 μl of the samples and 0.32 U of L-glutamate dehydrogenase. Protein was determined by the Bradford method.
Kinetic constants for 3-methyl-2-oxobutyrate, 4-methyl-2-oxopentanoate and 3-methyl-2-oxopentanoate were determined with substrate concentrations up to 15 mM, 10 mM L-glutamate and 0.7–4 μg protein using standard conditions. Kinetic constants for L-glutamate were determined with substrate concentrations up to 80 mM, 5 mM 3-methyl-2-oxobutanoic acid and 2.7–9.1 μg protein using standard conditions. The substrate specificity was tested with 1 or 10 mM amino acceptors, 2.5 μg protein using standard conditions.
The reverse reaction was measured as L-valine dependent decrease of AKG. The assay mixture contained 0.1 M Tris-HCl, pH 8.0, 0.1 mM PLP, 10 mM L-valine, 10 mM AKG and 5.3 μg protein. Standard conditions were used for further incubation and quantification.
AF0933 activity was also measured as 3-methyl-2-oxobutyrate dependent formation of pyruvate from L-alanine. The assay mixture contained 0.1 M Tris-HCl, pH 8.0, 0.1 mM PLP, 5 mM 3-methyl-2-oxobutyrate, up to 20 mM L-alanine and 2.7 μg protein. Samples were incubated at 70°C up to 20 min, the reaction was stopped on ice and the amount of pyruvate formed was quantified by measuring the reduction of pyruvate to lactate at 340 nm. The quantification mixture (0.3 ml) was incubated at 20°C for 10 min and contained 0.1 M Tris-HCl, pH 7.5, 0.3 mM NADH, 15–100 μl of the samples and 1.8 U of L-lactate dehydrogenase (Roche, Burgess Hill, UK).
The AF0933 activity with (R) and (S)-MBA was determined by the acetophenone assay. The assay mixture of 1 ml contained 0.1 M Tris-HCl, pH 8.0, 0.1 mM PLP, 2 mM 3-methyl-2-oxobutyrate, 2-oxobutyrate or pyruvate as amino acceptor, 2 mM (R)-MBA or (S)-MBA as amino donor and 16 μg protein. Samples were incubated at 65°C up to 180 min, the reaction was stopped on ice and the absorbance was measured at 245 nm. The molar extinction coefficient of 11.6 mM−1 cm−1 was used for acetophenone.
Effects of pH, Temperature, and Solvents on the Transamination Reaction
The pH optimum of the reaction for GEO1900 with 5 mM L-Leu, and 5 mM AKG was determined at 65°C using 50 mM sodium phosphate, pH 6.0–8.0, 50 mM BICINE containing 100 mM NaCl.
The temperature optimum of the transamination reaction between 5 mM L-Leu, and 5 mM AKG was determined in the temperature range from 40 to 95°C in 50 mM phosphate buffer, pH 7.9, supplemented with 100 mM NaCl. The thermostability of GEO1900 was determined by incubating a 1.0 mg/mL enzyme samples in 100 mM phosphate buffer, pH 8.0 at different temperatures between 30 and 70°C for 2 h. After the thermal treatment, the remaining enzyme activity was determined at 65°C as described above.
The pH optimum of AF0933 was determined with 3.6 μg protein, standard conditions and the following buffers each at a concentration of 0.1 M: MES (pH 5.5, 6.0, 6.5), Bis-Tris (pH 6.5, 7.0), HEPES (pH 7.0, 7.5, 8.0), Tris-HCl (pH 7.5, 8.0, 8.5), Bicine (pH 8.5, 9.0) and Piperazine (pH 9.0, 9.5, 10.0).
The temperature optimum of the reaction of AF0933 was determined between 50°C and 95°C in 0.1 M HEPES, pH 7.0, containing 0.1 mM PLP, 5 mM 3-methyl-2-oxobutyrate, 10 mM L-glutamate and 2.7–5.4 μg protein.
The thermostability of AF0933 was measured between 55 and 95°C, therefore protein was incubated in 3 volumes of buffer (0.1M HEPES, pH 7.0 adjusted to the respective temperature) up to 120 min. Samples (2.8 μg protein) were chilled on ice and remaining BCAT activity was determined in 0.1 M HEPES, pH 7.0 using standard conditions.
The influence of solvents (ethanol, DMSO, methanol and acetonitrile) on activity was analyzed between 10 and 50% solvent. The respective solvents were incubated with 2.2 μg protein (for AF0933) and 50 μg protein (for GEO1900) for 1 h at 20°C and 55°C then remaining activity was measured with standard conditions.
Protein Crystallization
For crystallization fractions of GEO1900 after the size-exclusion chromatography step were concentrated up to 10 mg/ml using an Amicon Ultra 30 kDa MWCO (Merck Millipore, Darmstadt, Germany) centrifugal filter units and stored at −70°C in 30 mM Tris-HCl buffer, pH 8.0, supplemented with 300 mM NaCl and 100 μM PLP. High-throughput crystallization screening was setup using a robotic system (Rigaku Automation, USA) at the Resource Centers division of the NBICS Centre of the National Research Centre “Kurchatov Institute” (Boyko et al.,
Prior to crystallization the AF0933 enzyme was concentrated to ~10 mg ml−1 using a 10 kDa Vivaspin membrane (Vivaproducts, Littleton, Massachusetts, USA). Microbatch crystallization trials were set up using an Oryx 6 crystallization robot (Douglas Instruments, Hungerford, England) using the JCSG+ screen (Molecular Dimensions, Newmarket, England) (Newman et al.,
Crystals of the AF0933 holoenzyme were grown using several different PEG and methylpentanediol (MPD) conditions. The best crystals grew in 100 mM cacodylate buffer pH 6.5, 200 mM magnesium chloride and 50% v/v PEG200. To obtain the inhibitor bound complex AF0933 was co-crystallized with 5 mM gabaculine under the same conditions.
Data Collection and Processing
The GEO1900 data were collected at the BL41XU beamline of the Spring8 synchrotron (Harima, Japan) using a Pilatus detector. The data were indexed, integrated, and scaled using XDS (Kabsch,
Data for AF0933 were collected on beamline I04-1 at the Diamond Synchrotron light source (Oxford, UK) under cryo conditions (100 K in a stream of gaseous nitrogen). Data were processed using XDS (Kabsch,
Structure Solution and Refinement
The structure of GEO1900_AKG was solved at 2.2 Å resolution by the molecular replacement (MR) method using the BALBES MR pipeline (Long et al.,
The AF0933 holo-structure was phased by MR implemented in MOLREP (Vagin and Teplyakov, 2010) using a sequence modified model (Lebedev et al.,
Refinement of all structures was carried out with REFMAC5 (Murshudov et al.,
Structure Analysis and Validation
PROCHECK was used for assessment of the quality of the models (Laskowski et al.,
The atomic coordinates and structure factors for both GEO1900 and AF0933 have been deposited in the Protein Data Bank. The PDB codes for GEO1900 holo enzyme as 5CM0 and GEO1900 complexed with AKG 5E25. The PDB codes for AF0933 holo enzyme 5MQZ and AF0933 complexed with gabaculine 5MR0.
Statements
Author contributions
EB, NR, AM, PJ, and CS performed the discovery, cloning, and expression of the novel enzymes, while EB, PJ, CS, AN, and KB performed the protein purification, crystallization, and structure determination studies with MI. KB carrying out data analysis and structural refinement. AN, TS, MS, and J-MS performed the biochemical characterization of functional properties and the substrate specificity studies with PS. JL, VP, PS, MI, and EB coordinated the work and wrote the manuscript with input from all authors.
Funding
This work was funded by the ERA-IB project THERMOGENE. The crystallographic study of AF0933 was carried out at Exeter, UK with support from the University of Kiel, Germany for the enzyme kinetic analysis. The Exeter JAL group would like to thank the BBSRC funded grant BB/L002035/1 and the University of Exeter for support.
Acknowledgments
The PS group would like to thank the BMBF grant 031A22 and the University of Kiel for support. The Exeter group would also like to thank the Diamond Synchrotron Light Source for access to beamline I04-1 (proposals No. MX8889) and the beamline scientists for assistance. The biochemical and crystallographic studies of GEO1900 were carried out in the Research Center of Biotechnology of the Russian Academy of Science in Moscow who would like to thank support from the Russian Federal Space Agency, and the Russian Science Foundation (project No 14-24-00172, in part of protein purification and structural work).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2019.00007/full#supplementary-material
References
1
AltschulS.GishW.MillerW.MyersE.LipmanD. (1990). Basic local alignment search tool. J. Mol. Biol.215, 403–410. 10.1016/S0022-2836(05)80360-2
2
BattyeT. G.KontogiannisL.JohnsonO.PowellH. R.LeslieA. G. (2011). iMOSFLM: a new graphical interface for diffraction-image processing with MOSFLM. Acta Crystallogr. D Biol. Crystallogr.67, 271–281. 10.1107/S0907444910048675
3
BezsudnovaE. Y.BoykoK. M.PopovV. O. (2017). Properties of bacterial and archaeal branched-chain amino acid aminotransferases. Biochemistry13, 1572–1591. 10.1134/S0006297917130028
4
BoykoK. M.LipkinA. V.PopovV. O.KovalchukM. V. (2013). From gene to structure: the protein factory of the NBICS Centre of Kurchatov Institute. Crystallogr. Rep.58, 442–449. 10.1134/S106377451105004X
5
BoykoK. M.StekhanovaT. N.NikolaevaA. Y.MardanovA. V.RakitinA. L.RavinN. V.et al. (2016). First structure of archaeal branched-chain amino acid aminotransferase from Thermoproteus uzoniensis specific for L-amino acids and R-amines. Extremophiles20, 215–225. 10.1007/s00792-016-0816-z
6
BradfordM. M. (1976). Rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem.72, 248–254. 10.1016/0003-2697(76)90527-3
7
BraunsteinA. E.ShemyakinM. M. (1953). The theory of processes of amino acid metabolism catalysed by pyridoxal 5′-phosphate-dependent enzymes. Biokhimia18, 393–411.
8
BricogneG.BlancE.BrandlM.FlensburgC.KellerP.PaciorekP.et al. (2016). BUSTER Version 2.10.3.Cambridge: Global Phasing Ltd.
9
ChenC. D.LinC. H.ChuankhayanP.HuangY. C.HsiehY. C.HuangT. F.et al. (2012). Crystal structures of complexes of the branched-chain aminotransferase from Deinococcus radiodurans with alpha-ketoisocaproate and L-glutamate suggest the radiation resistance of this enzyme for catalysis. J. Bacteriol.194, 6206–6216. 10.1128/JB.01659-12
10
ChoB. K.SeoJ. H.KangT. W.KimB. G. (2003). Asymmetric synthesis of L-homophenylalanine by equilibrium-shift using recombinant aromatic L-amino acid transaminase. Biotechnol. Bioeng.83, 226–234. 10.1002/bit.10661
11
CowtanK. (2010). Recent developments in classical density modification. Acta Crystallogr. D Biol. Crystallogr.66, 470–478. 10.1107/S090744490903947X
12
DementievaI. S.ZakomirdinaL. N.SinitzinaN. I.AntsonA. A.WilsonK. S.IsupovM. N.et al. (1994). Crystallization and preliminary X-ray investigation of holotryptophanases from Escherichia coli and Proteus vulgaris. J. Mol. Biol.235, 783–786. 10.1006/jmbi.1994.1033
13
DiederichsK.KarplusP. A. (1997). Improved R-factors for diffraction data analysis in macromolecular crystallography. Nat. Struct. Biol.4, 269–275. 10.1038/nsb0497-269
14
EmsleyP.LohkampB.ScottW. G.CowtanK. (2010). Features and development of Coot. Acta Crystallogr D Biol Crystallogr66, 486–501. 10.1107/S0907444910007493
15
EvansP. (2006). Scaling and assessment of data quality. Acta Crystallogr. D Biol. Crystallogr.62, 72–82. 10.1107/S0907444905036693
16
EvansP. R.MurshudovG. N. (2013). How good are my data and what is the resolution?Acta Crystallogr. D Biol. Crystallogr.69, 1204–1214. 10.1107/S0907444913000061
17
GotoM.MiyaharaI.HayashiH.KagamiyamaH.HirotsuK. (2003). Crystal structures of branched-chain amino acid aminotransferase complexed with glutamate and glutarate: true reaction intermediate and double substrate recognition of the enzyme. Biochemistry42, 3725–3733. 10.1021/bi026722f
18
GuanL. J.OhtsukaJ.OkaiM.MiyakawaT.MaseT.ZhiY.et al. (2015). A new target region for changing the substrate specificity of amine transaminases. Sci. Rep.5:10753. 10.1038/srep10753
19
HayashiH. (1995). Pyridoxal enzymes mechanistic diversity and uniformity. J. Biochem.118, 463–473. 10.1093/oxfordjournals.jbchem.a124931
20
HirotsuK.GotoM.OkamotoA.MiyaharaI. (2005). Dual substrate recognition of aminotransferases. Chem. Rec.5, 160–172. 10.1002/tcr.20042
21
HöhneM.SchatzleS.JochensH.RobinsK.BornscheuerU. T. (2010). Rational assignment of key motifs for function guides in silico enzyme identification. Nat. Chem. Biol.6, 807–813. 10.1038/nchembio.447
22
HumbleM. S.CassimjeeK. E.HåkanssonM.KimbungY. R.WalseB.AbediV.et al (2012). Crystal structures of the Chromobacterium violaceum ω-transaminase reveal major structural rearrangements upon binding of coenzyme PLP. FEBS J.279, 779–792. 10.1111/j.1742-4658.2012.08468.x
23
HutchinsonE. G.ThorntonJ. M. (1996). PROMOTIF—a program to identify and analyze structural motifs in proteins. Protein Sci.5, 212–220. 10.1002/pro.5560050204
24
HutsonS. (2001). Structure and function of branched-chain aminotransferases. Prog. Nucl. Acid Res. Mol. Biol.70, 175–206. 10.1016/S0079-6603(01)70017-7
25
InoueK.KuramitsuS.AkiK.WatanabeY.TakagiT.NishigaiM.et al. (1988). Branched-chain amino acid aminotransferase of Escherichia coli: overproduction and properties. J. Biochem.104, 777–784. 10.1093/oxfordjournals.jbchem.a122549
26
JansoniusJ. N. (1998). Structure evolution and action of vitamin B6-dependent enzymes. Curr. Opin. Struct. Biol.8, 759–769. 10.1016/S0959-440X(98)80096-1
27
KabschW. (2010). XDS. Acta Crystallogr. D Biol. Crystallogr.66, 125–132. 10.1107/S0907444909047337
28
KabschW.SanderC. (1983). Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features. Biopolymers22, 2577–2637. 10.1002/bip.360221211
29
KandaM.HoriK.KurotsuT.OhgishiK.HanawaT.SaitoY. (1995). Purification and properties of branched chain amino acid aminotransferase from gramicidin S-producing Bacillus brevis. J. Nutr. Sci. Vitaminol.41, 51–60. 10.3177/jnsv.41.51
30
KarplusP. A.DiederichsK. (2012). Linking crystallographic model and data quality. Science336, 1030–1033. 10.1126/science.1218231
31
KaulmannU.SmithiesK.SmithM. E. B.HailesH. C.WardJ. M. (2007). Substrate spectrum of ω-transaminase from Chromobacterium violaceum DSM30191 and its potential for biocatalysis. Enzyme Microb. Technol.41, 628–637. 10.1016/j.enzmictec.2007.05.011
32
KlenkH. P.ClaytonR. A.TombJ. F.WhiteO.NelsonK. E.KetchumK. A.et al. (1997). The complete genome sequence of the hyper-thermophilic, sulphate-reducing archaeon Archaeoglobus fulgidus. Nature390, 364–370. 10.1038/37052
33
KrissinelE.HenrickK. (2004). Secondary-structure matching (SSM), a new tool for fast protein structure alignment in three dimensions. Acta Crystallogr. D Biol. Crystallogr.60, 2256–2268. 10.1107/S0907444904026460
34
KrissinelE.HenrickK. (2007). Inference of macromolecular assemblies from crystalline state. J. Mol. Biol.372, 774–797. 10.1016/j.jmb.2007.05.022
35
LaskowskiR. A.MacArthurM. W.MossD. S.ThorntonJ. M. (1993). PROCHECK: a program to check the stereochemical quality of protein structures. J. Appl. Cryst.26, 283–291. 10.1107/S0021889892009944
36
LaskowskiR. A.SwindellsM. B. (2011). LigPlot+: multiple ligand-protein interaction diagrams for drug discovery. J. Chem. Inf. Model.51, 2778–2786. 10.1021/ci200227u
37
LebedevA. A.VaginA. A.MurshudovG. N. (2008). Model preparation in MOLREP and examples of model improvement using X-ray data. Acta Crystallogr. D Biol. Crystallogr.64, 33–39. 10.1107/S0907444907049839
38
Lee-PengF.-C.HermodsonM. A.KohlhawG. B. (1979). Transaminase B from Escherichia coli: quaternary structure, amino-terminal sequence, substrate specificity, and absence of a separate valine-α-ketoglutarate activity. J. Bacteriol.139, 339–345.
39
LeuchtenbergerW.HuthmacherK.DrauzK. (2005). Biotechnological production of amino acids and derivatives: current status and prospects. Appl. Microbiol. Biotechnol.69, 1–8. 10.1007/s00253-005-0155-y
40
LittlechildJ.NovakH.JamesP.SayerC. (2013). Thermophilic microbes in environmental and industrial biotechnology, in Biotechnology of Thermophiles, eds SatyanarayanaT.KawarabayasiY.LittlechildJ. (Dordrecht: Springer), 481–509.
41
LittlechildJ. A.GuyJ.ConnellyS.MallettL.WaddellS.RyeC. A.et al. (2007). Natural methods of protein stabilization: thermostable biocatalysts. Biochem. Soc. Trans.35, 1558–1563. 10.1042/BST0351558
42
LongF.VaginA. A.YoungP.MurshudovG. N. (2007). BALBES: a molecular-replacement pipeline. Acta Crystallogr. D Biol. Crystallogr.64, 125–132. 10.1107/S0907444907050172
43
ŁyskowskiA.GruberC.SteinkellnerG.SchürmannM.SchwabH.GruberK.et al. (2014). Crystal structure of an (R)-selective ω-transaminase from Aspergillus terreus. PLoS ONE9:e87350. 10.1371/journal.pone.0087350
44
MalikM. S.ParkE. S.ShinJ. S. (2012). Features and technical applications of ω-transaminases. Appl. Microbiol. Biotechnol.94, 1163–1171. 10.1007/s00253-012-4103-3
45
MardanovA. V.GumerovV. M.BeletskyA. V.ProkofevaM. I.Bonch-OsmolovskayaE. A.RavinN. V.et al. (2011). Complete genome sequence of the thermoacidophilic crenarchaeon Thermoproteus uzoniensis 768-20. J. Bacteriol.193, 3156–3157. 10.1128/JB.00409-11
46
McNicholasS.PottertonE.WilsonK. S.NobleM. E. M. (2011). Presenting your structures: the CCP4mg molecular-graphics software. Acta Crystallogr. Sect. D Biol. Crystallogr.67, 386–394. 10.1107/S0907444911007281
47
MehtaP. K.HaleT. I.ChristenP. (1993). Aminotransferases: demonstration of homology and division into evolutionary subgroups. Eur. J. Biochem.214, 549–561. 10.1111/j.1432-1033.1993.tb17953.x
48
MetzlerD. E.IkawaM.SnellE. E. (1954). A general mechanism for vitamin B6-catalyzed reactions. J. Am. Chem. Soc.76, 648–652. 10.1021/ja01632a004
49
MidelfortK. S.KumarR.HanS.KarmilowiczM. J.McConnellK.GehlhaarD. K.et al. (2013). Redesigning and characterizing the substrate specificity and activity of Vibrio fluvialis aminotransferase for the synthesis of imagabalin. Protein Eng. Des. Sel.26, 25–33. 10.1093/protein/gzs065
50
MurshudovG. N.SkubákP.LebedevA. A.PannuN. S.SteinerR. A.NichollsR. A.et al. (2011). REFMAC5 for the refinement of macromolecular crystal structures. Acta Crystallogr. D Biol. Crystallogr.67, 355–367. 10.1107/S0907444911001314
51
MuttiF. G.FuchsC. S.PressnitzD.SattlerJ. H.KroutilW. (2011). Stereoselectivity of four (R)-selective transaminases for the asymmetric amination of ketones. Adv. Synth. Catal.353, 3227–3233. 10.1002/adsc.201100558
52
NewmanJ.EganD.WalterT. S.MegedR.BerryI.Ben JelloulM.et al. (2005). Towards rationalization of crystallisation screening for small- to medium-sized laboratories: the PACT/JCSG+ strategy. Acta Crystallogr. D Biol. Crystallogr.61, 1426–1431. 10.1107/S0907444905024984
53
NortonJ. E.SokatchJ. R. (1970). Purification and partial characterization of the branched chain amino acid transaminase of Pseudomonas aeruginosa. Biochim. Biophys. Acta206, 261–269. 10.1016/0005-2744(70)90109-9
54
OkadaK.HirotsuK.HayashiH.KagamiyamaH. (2001). Structures of Escherichia coli branched-chain amino acid aminotransferase and its complexes with 4-methylvalerate and 2-methylleucine: induced fit and substrate recognition of the enzyme. Biochemistry40, 7453–7463. 10.1021/bi010384l
55
PadillaJ. E.YeatesT. O. (2003). A statistic for local intensity differences: robustness to anisotropy and pseudo-centering and utility for detecting twinning. Acta Crystallogr. D Biol. Crystallogr.59, 1124–1130. 10.1107/S0907444903007947
56
PannuN. S.MurshudovG. N.DodsonE. J.ReadR. J. (1998). Incorporation of prior phase information strengthens maximum-likelihood structure refinement. Acta Crystallogr. D Biol. Crystallogr.54, 1285–1294. 10.1107/S0907444998004119
57
PeisachD.ChipmanD. M.Van OphemP. W.ManningJ. M.RingeD. (1998). Crystallographic study of steps along the reaction pathway of D-amino acid aminotransferase. Biochemistry37, 4958–4967. 10.1021/bi972884d
58
PressnitzD.FuchsC. S.SattlerJ. H.KnausT.MacherouxP.MuttiF. G.et al. (2013). Asymmetric amination of tetralone and chromanone derivatives employing ω-transaminases. ACS Catal.3, 555–559. 10.1021/cs400002d
59
PuntaM.CoggillP. C.EberhardtR. Y.MistryJ.TateJ.BoursnellC.et al. (2012). The Pfam protein families database. Nucleic Acids Res.14, D290–D301. 10.1093/nar/gkr1065
60
QuerellouJ.Bonch-OsmolovskayaE. A.SlobodkinA. I. (2009). Geoglobus acetivorans sp. nov., an iron (III)-reducing archaeon from a deep-sea hydrothermal vent. Int. J. Syst. Evol. Microbiol.59, 2880–2883. 10.1099/ijs.0.011080-0
61
RamakrishnanC.RamachandranG. N. (1965). Stereochemical criteria for polypeptide and protein chain conformations. II. Allowed conformations for a pair of peptide units. Biophys. J.5, 909–933. 10.1016/S0006-3495(65)86759-5
62
RichardsonJ. S. (1981). The anatomy and taxonomy of protein structure. Adv. Prot. Chem.34, 167–339. 10.1016/S0065-3233(08)60520-3
63
RobertX.GouetP. (2014). Deciphering key features in protein structures with the new ENDscript server. Nucleic Acids Res.42, 320–324. 10.1093/nar/gku316
64
RuzickaF. J.FreyP. A. (2010). Kinetic and spectroscopic evidence of negative cooperativity in the action of lysine 2,3-aminomutase. J. Phys. Chem. B114, 16118–16124. 10.1021/jp103856m
65
SavileC. K.JaneyJ. M.MundorffE. C.MooreJ. C.TamS.JarvisW. R.et al. (2010). Biocatalytic asymmetric synthesis of chiral amines from ketones applied to sitagliptin manufacture. Science329, 305–309. 10.1126/science.1188934
66
SayerC.BommerM.IsupovM.WardJ.LittlechildJ. (2012). Crystal structure and substrate specificity of the thermophilic serine:pyruvate aminotransferase from Sulfolobus solfataricus. Acta Crystallogr. D Biol. Crystallogr.68, 763–772. 10.1107/S0907444912011274
67
SayerC.IsupovM. N.WestlakeA.LittlechildJ. A. (2013). Structural studies of Pseudomonas and Chromobacterium ω-aminotransferases provide insights into their differing substrate specificity. Acta Crystallogr. D Biol. Crystallogr.69, 564–576. 10.1107/S0907444912051670
68
SayerC.Martinez-TorresR. J.RichterN.IsupovM. N.HailesH. C.LittlechildJ. A.et al. (2014). The substrate specificity, enantioselectivity and structure of the (R)-selective amine : pyruvate transaminase from Nectria haematococca. FEBS J.281, 2240–2253. 10.1111/febs.12778
69
SchätzleS.Steffen-MunsbergF.ThontowiA.HöhneM.RobinsK.et al. (2011). Enzymatic asymmetric synthesis of enantiomerically pure aliphatic, aromatic and arylaliphatic amines with (R)-selective amine transaminases. Adv. Synth. Catal.353, 2439–2445. 10.1002/adsc.201100435
70
ShinJ. S.YunH.JangJ. W.ParkI.KimB. G. (2003). Purification, characterization, and molecular cloning of a novel amine:pyruvate transaminase from Vibrio fluvialis JS17. Appl. Microbiol. Biotechnol.61, 463–471. 10.1007/s00253-003-1250-6
71
StetefeldJ.JennyM.BurkhardP. (2006). Intersubunit signaling in glutamate-1-semialdehyde-aminomutase. Proc. Natl. Acad. Sci. U.S.A.103, 13688–13693. 10.1073/pnas.0600306103
72
StetterK. O. (1988). Archaeoglobus fulgidus gen. nov., sp. nov. a new taxon of extremely thermophilic archaebacteria. Syst. Appl. Microbiol.10, 172–173. 10.1016/S0723-2020(88)80032-8
73
TachikiT.TochikuraT. (1973). Separation of L-leucine-pyruvate and L-leucine-alpha-ketoglutarate transaminases in Acetobacter suboxydans and identification of their reaction products. Agric. Biol. Chem.37, 1439–1448. 10.1080/00021369.1973.10860842
74
TaylorP. P.PantaleoneD. P.SenkpeilR. F.FotheringhamI. G. (1998). Novel biosynthetic approaches to the production of unnatural amino acids using transaminases. Trends Biotechnol.16, 412–418. 10.1016/S0167-7799(98)01240-2
75
ThomsenM.SkaldenL.PalmG. J.HöhneM.BornscheuerU. T.HinrichsW. (2014). Crystallographic characterization of the (R)-selective amine transaminase from Aspergillus fumigatus. Acta Crystallogr. D Biol. Crystallogr.70, 1086–1093. 10.1107/S1399004714001084
76
TufvessonP.Lima-RamosJ.JensenJ. S.Al-HaqueN.NetoW.WoodleyJ. M. (2011). Process considerations for the asymmetric synthesis of chiral amines using transaminases. Biotechnol. Bioeng.108, 1479–1493. 10.1002/bit.23154
77
UchidaY.HayashiH.WashioT.YamasakiR.KatoS.OikawaT. (2014). Cloning and characterization of a novel fold-type I branched-chain amino acid aminotransferase from the hyper-thermophilic archaeon Thermococcus sp. CKU-1 Extremophiles18, 589–602. 10.1007/s00792-014-0642-0
78
VaginA.TeplyakovA. (2010). Molecular replacement with MOLREP. Acta Crystallogr. D Biol. Crystallogr.66, 22–25. 10.1107/S0907444909042589
79
VaguineA. A.RichelleJ.WodakS. J. (1999). SFCHECK: a unified set of procedures for evaluating the quality of macromolecular structure-factor data and their agreement with the atomic model. Acta Crystallogr. D Biol. Crystallogr.55, 191–205. 10.1107/S0907444998006684
80
van OosterwijkN.WilliesS.HekelaarJ.Terwisscha van ScheltingaA. C.TurnerN. J.DijkstraB. W. (2016). Structural basis of the substrate range and enantioselectivity of two (S)-selective ω-transaminases. Biochemistry.55, 4422–4431. 10.1021/acs.biochem.6b00370
81
VriendG. (1990). WHAT IF: A molecular modeling and drug design program. J. Mol. Graph.8, 52–56. 10.1016/0263-7855(90)80070-V
82
WardJ.WohlgemuthR. (2010). High-yield biocatalytic amination reactions in organic synthesis. Curr. Org. Chem.14, 1914–1927. 10.2174/138527210792927546
83
WinnM. D.BallardC. C.CowtanK. D.DodsonE. J.EmsleyP.EvansP. R.et al. (2011). Overview of the CCP4 suite and current developments. Acta Crystallogr. D Biol. Crystallogr.67, 235–242. 10.1107/S0907444910045749
84
WinterG.LobleyC. M. C.PrinceS. M. (2013). Decision making in xia2. Acta Crystallogr. D Biol. Crystallogr.69, 1260–1273. 10.1107/S0907444913015308
85
YennawarN.DunbarJ.ConwayM.HutsonS.FarberG. (2001). The structure of human mitochondrial branched-chain aminotransferase. Acta Crystallogr. D Biol. Crystallogr.57, 506–515. 10.1107/S0907444901001925
86
YuX.WangX.EngelP. C. (2014). The specificity and kinetic mechanism of branched-chain amino acid aminotransferase from Escherichia coli studied with a new improved coupled assay procedure and the enzyme's potential for biocatalysis. FEBS J.281, 391–400. 10.1111/febs.12609
Summary
Keywords
thermophilic archaea, branched-chain aminotransferases, substrate specificity, X-ray structural analysis, biocatalysis
Citation
Isupov MN, Boyko KM, Sutter J-M, James P, Sayer C, Schmidt M, Schönheit P, Nikolaeva AY, Stekhanova TN, Mardanov AV, Ravin NV, Bezsudnova EY, Popov VO and Littlechild JA (2019) Thermostable Branched-Chain Amino Acid Transaminases From the Archaea Geoglobus acetivorans and Archaeoglobus fulgidus: Biochemical and Structural Characterization. Front. Bioeng. Biotechnol. 7:7. doi: 10.3389/fbioe.2019.00007
Received
26 October 2018
Accepted
08 January 2019
Published
24 January 2019
Volume
7 - 2019
Edited by
Roberto Contestabile, Sapienza University of Rome, Italy
Reviewed by
Stefano Bruno, University of Parma, Italy; Samanta Raboni, University of Parma, Italy
Updates

Check for updates
Copyright
© 2019 Isupov, Boyko, Sutter, James, Sayer, Schmidt, Schönheit, Nikolaeva, Stekhanova, Mardanov, Ravin, Bezsudnova, Popov and Littlechild.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Vladimir O. Popov vpopov@exchange.inbi.ras.ruJennifer A. Littlechild j.a.littlechild@exeter.ac.uk
This article was submitted to Bioprocess Engineering, a section of the journal Frontiers in Bioengineering and Biotechnology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.