Determination of the Structure of the Catabolic N-Succinylornithine Transaminase (AstC) from Escherichia coli

Escherichia coli possesses two acyl ornithine aminotransferases, one catabolic (AstC) and the other anabolic (ArgD), that participate in L-arginine metabolism. Although only 58% identical, the enzymes have been shown to be functionally interchangeable. Here we have purified AstC and have obtained X-ray crystal structures of apo and holo-AstC and of the enzyme complexed with its physiological substrate, succinylornithine. We compare the structures obtained in this study with those of ArgD from Salmonella typhimurium obtained elsewhere, finding several notable differences. Docking studies were used to explore the docking modes of several substrates (ornithine, succinylornithine and acetylornithine) and the co-substrate glutamate/α-ketogluterate. The docking studies support our observations that AstC has a strong preference for acylated ornithine species over ornithine itself, and suggest that the increase in specificity associated with acylation is caused by steric and desolvation effects rather than specific interactions between the substrate and enzyme.


Introduction
The catabolism of L-arginine by bacteria has been an area of extensive study. In part, these studies have been driven by the role of L-arginine catabolism in various potential applications. For example, disruption of L-arginine catabolism can lead to strains that over-produce L-arginine and may be useful for the multibillion dollar amino acid production industry [1] and L-arginine can be used as a feedstock for the synthesis of polyamines such as 1,4-diaminobutane [2] and putrescine [3], which have utility in polymer production. Enzymes used in the catabolism of arginine, and related enzymes, have other potential biotechnological applications, particularly in the development of efficient routes for the production of chiral amines, used in the production of pharmaceutical products [4].
The arginine-succinyl transferase pathway of E. coli provides nitrogen for growth via a series of succinylated intermediates ( Fig. 1) [16]. The pathway is encoded by the astCADBE operon, which is highly regulated in response to nitrogen and general starvation through the coordinated actions of the transcription factors RpoN, RpoS, NtrC, CRP and ArgR [17].
The second E. coli SOAT is the anabolic ArgD, which is 58.6% identical to AstC. ArgD is required for the biosynthesis of arginine via a series of N-acetyl intermediates, including N-acetylornithine. In 1998, Ledwidge and Blanchard showed that ArgD is identical to the long-sought dapC gene product (DapATase) which catalyzes the glutamate-dependent transamination of N-succinyl-L-2-amino-6-oxopimelate to generate N-succinyl-L,L-diaminopimelate, required for the biosynthesis of lysine [19]. X-ray structures for ArgD from a number of different bacteria including Salmonella typhimurium (e.g. PDB code 2PB2 and 2PB0), Thermus thermophilus (PDB code 1VEF), Aquifex aeolicus VF5 (PDB code 2EH6), Campylobacter jejuni (PDB code 3NX3) and Thermotoga maritima (PDB code 2ORD) have been determined, demonstrating that this enzyme is also a classical dimeric PLP-dependent, ab-fold protein.
ArgD from E. coli and S. typhimurium have been characterized biochemically and shown to have k cat values of 0.22-0.82 s 21 and 0.61-1.6 s 21 for ornithine (ORN) and acetylornithine (AO), respectively, and K M values of 4,500-640 mM and 150-37 mM for ORN and AO, respectively [19,20] . The equilibrium constant (K eq ) for this reaction is close to 1, suggesting that the enzyme can operate efficiently in either direction depending upon the relative concentrations of substrate and product. Indeed, AstC expressed from a plasmid is able to complement an argD deletion [16] [21].
Here we describe the structural and biochemical characterization of the catabolic SOAT from E. coli, AstC, finding differences from the structures of the related anabolic SOAT enzyme (ArgD). We also propose a plausible explanation for the greater specificity of AstC for AO compared with ORN itself.

Purification of AstC
The native, untagged form (i.e. no plasmid or clone used) of AstC from E. coli (BL21 DE3) was purified from a bacterial growth in a rich medium (2xYT) using affinity (Ni-NTA) and size exclusion chromatography. Although AstC is not one of the 10 E. coli proteins found to bind tightly to Ni-NTA columns by Tiwari et al, it is likely that under the growth conditions used in their study AstC was not highly expressed [22]. We estimate the purity of the protein after these two steps to be .95% (from Coomassie stained SDS PAGE analysis, data not shown). The protein was subsequently concentrated in a 10 kDa Amicon spin concentrator to approximately 5 mg/mL for crystallization trials. The concentrated protein solution was clear and colorless. Holo-enzyme was obtained by adding an excess of PLP (,10 fold) to the apo-enzyme and incubating this on ice for several hours.

Crystallography
A long, needle-like crystal of the PLP-free enzyme grew from a reservoir containing 1.4 M sodium malonate (pH 7), 10% (v/v) malate-MES-Tris buffer [23] pH 5 at 20u C in droplets set up in a ratio of 50:25:25 protein:reservoir:silver bullet #62 (HR2-996-62, Hampton Research, USA) with a total volume of 400 nL -see Figure S1-A. The Schiff base-enzyme complex (PLP bound) crystals were grown at 8u C in 1.17 M sodium malonate (pH 7), 0.09 M Tris chloride pH 8.0, 0.02 M di-ethylammonium formate in drops that were 300 nL initial volume, consisting of equal volumes of protein and reservoir -see Figure S1-B. The crystals with PLP bound in the protein, but the Schiff base broken were grown under similar conditions to those with the Schiff base intact, but ORN as a dry solid was added to the crystallization drop 2 days prior to harvesting the crystals and collecting diffraction data. This same procedure was used to add aKG to holo-AstC crystals. 10 mM glutamic acid in crystallization cryo-solution (2.2 M sodium malonate pH 7) was added to crystals 3 days before data collection. Co-crystals of AstC in the presence of PLP and SO were grown at 20u C as a 50:50 ratio of the protein complex with 1-4-1.6 M ammonium sulfate, and pH 7-9 (300 nL initial drop size), see Figure S1-C. All of the crystallizations were set up in SD-2 (IDEX Corp, USA) 96-well sitting drop plates using a Phoenix robot (Art Robbins Industries, USA). Crystals grown with the PLP-enzyme complex were chunkier and grew more reproducibly than the crystals of the apo enzyme.
All X-ray data were collected at the MX-2 beamline of the Australian Synchrotron (see Table 1). The non-PLP crystals adopted the R32:h rhombohedral spacegroup, whereas all crystals grown with added PLP were found the C2 spacegroup. There were 2 molecules in the asymmetric unit of the H32 crystals and 4 molecules in the asymmetric unit of the C2 crystals. Multiple data sets were collected for almost all soaks and co-crystals, and generally 360 degrees of data were collected in 0.5 degree oscillation increments for a total of 720 frames at 13,000 eV X-ray energy. The crystals grown in sodium malonate had no additional cryo-protectant added and were flash-cooled directly in the nitrogen cryo-stream. Crystals grown in ammonium sulfate were cryo-protected by adding 1.5 mL of 2.5 M sodium malonate to the crystallization drop before harvesting the crystals and cryo-cooling in the cold nitrogen stream. The data were indexed using XDS [24] and scaled using SCALA [25](CCP4 software suite [26]). The initial non-PLP bound structure was solved by molecular replacement using Phaser [27] with PDB 2PB2. Subsequent structures were solved using the A protomer from the non-PLP bound structure as the starting model by molecular replacement with Phaser. All structures were built using Coot [28] and refined with Refmac [29].

Synthesis of Succinylornithine
To a solution of Boc-L-Orn-OH (200 mg, 0.86 mmol) in a 1:1 mixture of water and dimethylformamide (8 mL) was added triethylamine (480 mL, 3.44 mmol) and succinic anhydride (103 mg, 1.03 mmol). The colorless solution was left to stir at room temperature under nitrogen atmosphere for 14 h and was concentrated under reduced pressure to give a white colored solid. This material was subjected to flash chromatography (silica, 5:95 acetic acid: 10 % methanol in dichloromethane) to yield Boc protected SO as a white colored foam (246 mg, 86%).  Figure S2 for a schematic representation of the reaction.

Biochemical Assays
Transaminase activity was measured by glutamate dehydrogenase (GDH) assay [30,31]. Briefly, the rate of glutamate formation via the amination of aKG was monitored by the inclusion of glutamate dehydrogenase and NAD + in the reaction, effectively recycling the aKG and reducing one molecule of NAD + to NADH per glutamate formed. The concentration of NADH was monitored at 340 nm. Reactions were conducted in 10 mM Tris-HCl (pH 8.0) containing 200 mM aKG, 500 mM NAD + , 5 U GDH and 0-2,000 mM ORN, SO, AO. The reactions were carried out in 100 mL volumes in a 96-well, UV-transparent microtiter plate and the change of absorbance at 340 nm was followed using a SpectroMax 190 spectrophotometer reading each well once every 10 seconds. Concentrations of aKG above 200 mM inhibited GDH activity ( Figure S3), preventing use of aKG at higher concentrations. However, the rate of turnover by GDH far exceeded the rate of turnover by AstC (.1,000-fold), as confirmed using ammonia as a positive control (data not shown). The concentration of aKG was therefore assumed to be non-limiting throughout the assays.

Thermal Denaturation Studies
The AstC protein was used at 110 mM (5 mg/mL) and four duplications of the apo-AstC enzyme and four duplications the holo-enzyme (PLP/enzyme complex) were set up. Each experimental well contained 0.3 mL protein +0.3 mL 1/10 SYPRO dye (Sigma S5692) +19.4 mL Tris buffered saline solution.

Modelling
All substrates were prepared by assigning protonation states using a neutral pH model as well as AM1BCC partial atomic charges [32] with QuacPac from OpenEye Scientific Software (Santa Fe, NM). Conformations for each ligand were enumerated using OpenEye's Omega software. Hydrogens were added to the crystal structures using Reduce [33]. All docking calculations were performed with FRED 2.2.5 [34]. Two constraints were used during the docking process. The constraints were added such that an amine, ketone, or carboxylic acid had to be positioned within 3 Å of the Schiff base and that the ligand must form a hydrogen bond with Arg141 (i.e. the interaction needed to position the ligand for catalysis as revealed in our structure with bound SO). These constraints assure that all successfully docked poses will be consistent with our knowledge of this enzyme's catalytic mechanism. Poses were scored and optimized with the Chemgauss3 scoring function [34] and the top 100 poses were ranked by a consensus of the following scoring functions: Chemgaus3, Piecewise Linear Potential (PLP), Shapegauss, Chemscore, OEChemscore, and Chemical Gaussian Overlay (CGO). The electrostatic potential for the protein active site was calculated using the Poisson-Boltzmann equation as implemented in the Zap toolkit from Openeye Scientific Software.

Analysis of Purified AstC
The native, non-tagged AstC that was purified from E. coli was the apo-protein, with no bound PLP, as suggested by the lack of color of the purified protein, and confirmed by a UV-Vis absorbance spectrum, which showed no absorbance in the region of 320 to 500 nm (data not shown) which is characteristic of PLP containing proteins [35]. This was unexpected, as PLP is not limiting under the normal growth conditions used here (2x YT media) and it is generally bound covalently to a lysine as an aldimine for most PLP enzymes (for example [20]).
It was possible to reconstitute the holo-enzyme by the addition of PLP (Sigma P9255) to the apo-protein. The reconstituted enzyme was more heat stable than the apo-protein; holo-AstC and apo-AstC had apparent T m values of 46uC and 43.5uC, respectively ( Figure 2).
The activity of the reconstituted holo-AstC was tested against ORN, SO and AO ( Table 2). The K M data for ORN and AO are in reasonable agreement with those recorded by Ledwidge and Blanchard [19], who found K M values of 4,500 and 150 mM for ORN and AO respectively (compared with 4,391 and 338 mM in this study). The k cat values from the Ledwidge study were significantly lower than those described in this study; 4.7 s 21 for AO and 3.5 s 21 for ORN from this study compared to 0.61 s 21 for AO and 0.22 s 21 for ORN [19]. One possible explanation for this discrepancy is that the ArgD used in the earlier study was not fully occupied with PLP.
The more recent study from Rajaram and co-workers with S. typhimurium ArgD (St-ArgD) [20] also reported steady state kinetic values for ORN (k cat = 0.82 s 21 , K M = 640 mM) and AO (k cat = 1.6 s 21 , K M = 37 mM). Although the values for K M obtained by Rajaram et al are somewhat lower than our own or those of Ledwidge and Blanchard, the order of substrate preference is the same in all three cases.
We also tested the activity of purified holo-AstC against its physiological substrate, SO. The

Structures of Apo-and Holo-AstC
We were able to obtain crystals of the apo-protein, although they diffracted to only medium resolution (2.75 Å ). The structure of apo-AstC is shown in Figures 3A and 3B. The A and B protomers in the asymmetric unit differ substantially, with the A chain having electron density for almost the entire length of the molecule (Gln3 to Ser402) and the B chain starting at residue Ala21 and ending at Ser402. The B chain also has no electron density corresponding to two loops in the A chain from His142 to Ala163 and Val274 to Gly283, respectively. The two protomers do not form a dimer, rather dimers are formed through crystallographic symmetry for both, the A with another A protomer and the B with another B protomer. The loops missing from the B protomer form part of the dimer interface; analysis of the PISA [36] output shows that the A-A interface buries about 4763 Å 2 whereas the B-B interface buries only 1906 Å 2 . Consistent with our spectrophotometric data, there was no electron density for the PLP co-factor in the crystal structures of protein obtained directly from purification.
Crystals of holo-AstC diffracted to a higher resolution than the apo-protein (2.20 Å ). Figures 4A and 4B show a dimer of holo-AstC, with PLP clearly bound in the cofactor site at the dimer interface. There is unambiguous density for the cofactor, and there is continuous density between Lys252 and the PLP, indicating that the expected aldimine had formed between the lysine residue and the PLP (see Figure 5A). There are four independent protomers in the asymmetric unit in the C2 crystal form, two full dimers, each protomer being almost equally resolved. The PLP makes interactions with other residues in the binding site: the N1 of PLP is 2.6 to 2.8 Å from the Asp223 sidechain carboxylate, the O3 is 2.9 to 3.0 Å to the Gln226 sidechain, the PO 3 oxygens are 2.9 to 3.1 Å to the Gly105 backbone nitrogen, 2.8 to 3.0 Å from the Ala106 backbone nitrogen, 2.9 to 3.1 Å from the Thr281 backbone nitrogen and 2.5 to 2.8 Å from the Thr281 sidechain oxygen (the distances vary slightly between the four protomers in the asymmetric unit). In this case, the A-B and C-D dimer interfaces are also effectively equal, with a PISA analysis showing that the average surface area buried is 4752 Å 2 (4764 Å 2 for A-B and 4740 Å 2 for C-D). Dynamic light scattering indicates a hydrodynamic radius mean value 4.5 nm, indicative of a molecular weight of approximately 100 kDa (data not shown), suggesting that AstC is also a dimer in solution.

Substrate/co-substrate Binding
To better understand the intricacies of the reaction, ORN, aKG and GLU were soaked into crystals of holo-AstC. The crystals soaked with GLU diffracted poorly (i.e. worse than 3 Å ), these crystals were not of sufficient quality for further analyses. We found very little change in the structures obtained from aKG soaked crystals and were unable to detect electron density corresponding to the co-substrate.
We were also unable to find extra electron density for the ORN soaked crystals, but there was a distinct change in the active site, with the PLP-Schiff base bond to Lys252 clearly broken (see Figure 5B). This finding is consistent with the completion of the first half reaction (i.e. deamination of the substrate), producing an aldehyde and pyridoxamine. The lack of electron density for the substrate may be due to the high K M of AstC for ORN ( Table 2).
Electron density corresponding to the external aldimine enzyme:substrate complex was present in the structures obtained from crystals co-crystallized with SO. In these crystals the Lys252-PLP Schiff base was absent and continuous electron density between PLP and SO was observed (see Figure 5C). The SO makes remarkably few direct interactions with the protein and it seems that it is held in place mainly by the covalent bond to the cofactor. In addition to this, one of the carboxylates makes a salt bridge to Arg141, which varies in length from 2.8 to 3.1 Å . The other interactions are through water molecules to the protein (0 to 2 per protomer). There were also some minor changes outside of the cofactor binding site. The Phe138 side chain moves slightly, there is some movement in the loop consisting of residues 151-154 and the side chain of Tyr18 moves to accommodate the SO moiety. Apart from these differences, all three structures are remarkably similar.
Using molecular docking we were able to predict plausible conformations for those ligands that were not observed in the Xray structures. SO, ORN and AO were docked into the active site of the structure with PLP covalently bound to Lys252, while aKG was docked into the non-covalently bound pyridoxamine phosphate structure. By constraining the docking process to produce poses that form hydrogen bonds with Arg141 and place either an amine, a ketone, or a carboxylate functional group within 3 Å of the catalytic nitrogen, we were able to obtain reasonable docked conformations for all molecules. Furthermore, using this method, the crystallographic ligand (SO) was docked within 1.8 Å of the crystallographic conformation ( Figure 6C). Additionally, for all but AO, the top ranking docked conformation was in an orientation consistent with the known catalytic mechanism. In the case of AO, the highest-ranking conformation consistent with the known mechanism was ranked 6 th (out of 100).
We see that aKG is positioned with the ketone carbon ready for attack by the amino nitrogen of pyridoxamine phosphate. Additionally, the carbonyl group is stabilized by interactions with Arg141 ( Figure 6B). For the ornithine-derived substrates we see the same trend. The side chain amine is positioned proximally to the Schiff nitrogen of the covalently bound PLP. The backbone portion of ORN and its acylated derivatives are stabilized by interactions with Arg141 ( Figure 6C). We also see for many docked conformations, Glu195 and Arg194 together serve to stabilize different charged groups near the Schiff base.
Although AstC can catalyze the reverse reaction (glutamate to aketoglutarate), in order to do so, GLU must bind to the PLP-bound state and interact with the aldimine. Upon docking GLU into this structure we find that the backbone carboxylate group interacts with Arg141 positioning the amino nitrogen for catalysis ( figure 6D). It may at first be surprising to see that GLU and aKG do not interact in the same fashion, but this seems reasonable considering that they must bind to two different states of the enzyme for catalysis.

AstC vs. ArgD
A comparison of the AstC structure with the ArgD structures from Samonella typhimurium (PDB codes 2PB0 and 2PB2) shows that the overall fold is the same but with several discrete differences  Figure 3A shows the superposition of apo-AstC protomer B (burnt orange) on protomer A (cyan) without PLP present in the structure. The N-terminus at the top right and the first 20 residues of protomer B are not modelled due to a lack of electron density for these residues. The other clear difference shown in 3A is loop 142-163, which is missing in protomer B (lower right). An approximate 90 degree rotation in 3B highlights the missing loop of residues 274-283 (lower left side of 3B). doi:10.1371/journal.pone.0058298.g003  Figures 4A and 4B show the native dimer structure as a ribbon Ca trace with the PLP cofactor shown in a magenta colored stick representation. Figure 4B  ( Figure 7A). The start of the structural alignment (as performed in Coot [28]) is around Ala21/Phe23 (A or B chains, respectively), although there are several conserved residues prior to this in the sequence. 244 out of 405/406 residues in the sequence are identical, with Ser402 being the last residue seen in the AstC structure corresponding to residue 405 in 2PB2 (the last 4 residues of AstC, Arg-Gly-Ser-Ser, are not conserved, do not have electron density and are not modelled). Although the sequence identity is less than 50% for the 100 residues in the C-termini, the backbones of the structures essentially overlap in this region.
The St-ArgD structures are most similar to the B protomer of the apo-enzyme structure of AstC reported here, in that the Nterminal ,23 amino acids of the structures are missing as are six residues in a loop from 277-282. These six residues are not modelled in three out of four independent St-ArgD monomers. The electron density for these residues is unambiguous in all of the  Figure 5A shows the 2Fo-Fc density (blue mesh) for the PLP cofactor bound to Lys252 in a Schiff base. The B protomer is in green and the A protomer loop of 277-284 is shown in magenta. The electron density is set at 1.5 sigma. Figure 5B shows the density for PLP and Lys252 after the crystals were soaked with ORN. It is clear that the covalent bond between Lys252 and the PLP is broken and we have now modelled the PLP with an additional oxygen atom. Figure 5C shows the electron density for Lys252, PLP and the SO found in the SO soaked crystals. There is clean electron density that shows the SO covalently bound to the PLP and a break from Lys252. All density figures are 2Fo-Fc maps and were made with the ligands in the model. doi:10.1371/journal.pone.0058298.g005 PLP-bound structures of AstC and the A protomer of the PLP-free structure. In the AstC structures, the loop 277-282 ventures into the catalytic site of the neighboring protomer and makes important interactions to the cofactor (the distance from the Thr281 backbone nitrogen to an oxygen of PLP varies from 2.75 to 3.10 Å and the hydroxyl of Thr281 has a hydrogen bond of 2.5 to 3.0 Å to a phosphate oxygen). The equivalent St-ArgD Thr284 (in the one loop that exists in 2PB0) is also within hydrogen bonding distance to an oxygen of the PLP phosphate (2.54 Å ). As this is part of the dimer interface, it is no surprise that the inferred dimer interfaces for St-ArgD are smaller than for AstC, with PISA giving 3600 Å 2 for 2PB0 and 2744 Å 2 for 2PB2 compared with 4752 Å 2 for holo-AstC.
The structure factors deposited in the RCSB for 2PB0/2PB2 give electron density maps that clearly show that the PLP is not covalently bound to Lys255 via a Schiff base and that the PLP cofactors are not fully occupied in the pocket of ArgD (see Figure  S4), in contrast to what is stated in the paper describing the structure [20]. This partial occupancy suggests that the recombinant St-ArgD was not fully saturated with PLP, and is consistent with the observation that AstC purified as an apo-protein, suggesting a low affinity of both ArgD and AstC for their PLP co-factor. These two observations lend credence to the suggestion that the earlier activity assay work of Ledwidge and Blanchard may have been affected by having limiting PLP in the protein sample. The lack of internal aldimine in the ArgD structures is similar to the structure of AstC after soaking with ORN, suggesting that the majority of St-ArgD had undergone a half reaction.
At the N-termini of the AstC and St-ArgD structures there are a number of differences ( Figure 7B). The AstC structures start at residue Gln3 and come into alignment with 2PB2 at residue Ala21 (Ala24 in the A chain of 2PB2). In the B chain of 2PB2, the model of St-ArgD starts at residue Leu18 and comes into alignment with AstC at residue Phe26/Ile27 (23/24 in the AstC numbering). An interesting feature is that the A chain of AstC at residue Trp13 This is the highest scoring pose that was correctly oriented for catalysis (6 out of 100). 6B: Shown here is the top scoring conformation of aketoglutarate (grey) docked into SO transaminase with pyridoxamine phosphate (magenta). Hydrogen bonds with Arg 141 are shown in black. 6C: Shown here is the top-scoring pose of SO (grey, SO-m) docked into the active site of SO transaminase. Pyridoxyl phosphate is shown in magenta covalently bound to Lys 252. Hydrogen bonds between SO and Arg 141 are shown in black. In cyan is the product succinyl-glutamic semialdehyde (SO-c, derived from the crystal structure by replacing the covalently bound amine group with the appropriate aldehyde). Hydrogen bonds from this ligand to Arg 141 are shown in cyan. In this image the docked structure (grey) has a root-mean-square distance of 1.8Å relative to the crystal structure (cyan). 6D: Shown here is glutamate (grey) docked into the active site of SO transaminase. Pyridoxyl phosphate is shown in magenta covalently bound to Lys 252. Hydrogen bonds between Arg 141 and glutamate are shown in black. This is the highest scoring pose that was correctly oriented for catalysis (2 out of 100). doi:10.1371/journal.pone.0058298.g006 starts to overlap with the model for the 2PB2 structure at Ser283B. That is to say, Ser283-Thr284 of the B chain of 2PB2 roughly occupies the same space as Trp13-Met14 of the A chain of AstC, before the 2PB2 structure moves off between two beta-strands and then rejoins the AstC model for chain B at residue Gly287 (Gly284 in AstC). The reciprocal relationship is also true but with an extra residue, where part of the N-terminus of the AstC B chain is overlapping with the middle of the A chain of 2PB2. We believe this is due to the missing loop in the 2PB2/2PB0 structures, residues 277-282, which are clearly defined in the AstC structures but have inadequate electron density in the 2PB2/2PB0 structures for accurate structural modelling. In the one monomer of 2PB0 where there is sufficient electron density to model these residues, it follows the same path as in the AstC structures. But where these residues are missing, the 2PB2/2PB0 models deviate from this path and start the model from a completely different place (the Nterminus of the dimer partner in the AstC structures).
In the analysis of the electron density for structures 2PB0/2PB2, we found that there is clear density for a valine at position 298 in the sequence (see Figure S5), which is defined as an alanine in the sequence given by Rajaram et al and is an alanine in the Uniprot database. We expect that one of two scenarios is possible: there was a PCR error during the cloning of this gene from the genomic DNA used, or that there was a mutation at this position in the S. typhimurium genomic DNA. As this is a single base mutation (GCX for Ala, GUX for Val) either possibility is reasonable.
There are five ArgD structures currently available in the PDB from various organisms (S. typhimurium, T. thermophilus, T. maritima, C. jejuni and A. aeolicus, see Table 3). Superposing these with the AstC structure (using the SSM algorithm as implemented in Coot) shows that the N-termini do not align well until about position 20 (AstC numbering) and the first helix varies somewhat, but past the first helix, the alignment is very good with only small differences in the majority of the structures (there is some variation in 7 different turn/loop regions: Ala93 to Ala96, Asp122 to Ser129, Gly149 to Pro160, Thr266 to Gly284, Lys298 to Asn306, Thr321 to Glu331 and Asn344 to Lys351). There is also some variation in the Ctermini. Of note is the region Thr266 to Gly284, which is the region that complements the binding site of the second protomer in the dimer. One structure, 3NX3, has this in an 'open' form and the 2PB2/2PB0 structures have this as disordered. Around the PLP site the residues are highly conserved, the only changes being Val225 is sometimes isoleucine and Ser104 is threonine in one instance (2ORD). The other changes nearby are: Asn109 to valine or isoleucine, His139 to serine and Tyr282 to phenylalanine, with these all being four or more Å ngstroms away.
From the biochemical data (Table 2), it appears that AstC and ArgD are biochemically similar, consistent with the ability of astC to complement an argD deletion [16]. Both have comparable steady state kinetic values, with a strong preference for the acylated ornithine derivatives AO (AstC and ArgD) and SO (for AstC, at least) over ORN itself.
When we look at the distribution of scores for the docked conformations, we find that we are able to reproduce the correct rank ordering of the ORN substrates relative to their experimental  Figure 7A shows the structure of 2PB0 (in burnt orange) with a dimer of AstC overlayed (cyan). The structures are very similar with the major differences being in the N-terminus (highlighted in 7B) and the loop of 277-286. Figure 7B again has the 2PB0 structure in orange and AstC in cyan, but the differences in 2PB0 are highlighted in red. On the top right hand side we have the start of the 2PB0 B chain, Leu18, which does not come into alignment with AstC until residue Phe26. At the top we have the two structures aligned at residue Phe276 where the density falls off for 2PB0 and there is a gap to residue Gly282 (hidden in the helix of the N-terminus of AstC). The 2PB0 model for this part of the chain comes back into alignment with AstC just before the helix at residue Gly287 (bottom of figure). doi:10.1371/journal.pone.0058298.g007  Figure S6, middle). Examination of the various terms in our scoring functions reveal that the increased affinity of the various ORN substrates results mostly from steric and desolvation effects (within the relative accuracy of the scoring function; Figure S6, bottom). This hypothesis is supported by the fact that the acetyl or succinyl moieties do not make any specific interactions with the active site, but rather are stabilized by positive electrostatic potential observed in this portion of the cavity ( Figure S7). This hypothesis is sensible when considering the fact that a promiscuous enzyme binding different substrates in similar orientations would likely need to rely on non-specific interactions to accommodate all of its substrates. Having numerous specific interactions such hydrogen bond donors/acceptors in the active site would incur an energetic penalty when a given ligand cannot satisfy these interactions. Thus the most reasonable strategy for positioning different ligands in the same orientation is to rely on non-specific interactions to stabilize the enzyme-substrate complex. This is consistent with what we observe. It seems that E. coli possesses two enzymes (AstC and ArgD) with interchangeable functionality and that the two have likely evolved from a common ancestral gene. Although the proteins are biochemically similar, they are components of differentially regulated transcriptional units, i.e. astCADBE and argD, and in these distinct regulatory contexts they perform different physiological roles. The I-A-D (innovation-amplification-divergence) model of evolution as proposed by Nä svall et al [37] proposes that a gene encoding a protein that has developed a promiscuous functionality may be duplicated, allowing both the physiological and promiscuous functionalities to become more efficient in their separate tasks [38][39][40]. Here, however, it is through divergence of the regulation of duplicate genes that each gene product has become adapted for specific physiological roles (anabolism and catabolism).
In conclusion, the E. coli AstC protein, involved in arginine catabolism, is a PLP containing transaminase that is structurally similar to the related protein ArgD, which is the transaminase enzyme used in arginine synthesis. The enzymes show similar kinetic values, although the current work does raise the possibility that previously published values of the ArgD kinetic parameters may be underestimated, due to a lack of saturating PLP in the enzyme samples. Docking studies show that AstC has a strong preference for acylated ornithine species over ORN, and suggest that the increase in specificity associated with acylation is caused by steric and desolvation effects rather than specific interactions between the substrate and enzyme.  Figure S4 The PLP density found for 2PB0. Shows the negative density associated with the PLP-Lys255 bond in the 2PB0 structure as deposited in the Protein Data Bank. The orientation of 2PB0 is roughly equivalent to the AstC structure in Figure 4A. 2Fo-Fc density is shown in the blue chicken wire contoured at 1.0 sigma, the Fo-Fc density map is shown in the red chicken wire contoured at 3.0 sigma. (TIF) Figure S5 Mutation found in the structure of 2PB0. This figure shows the 2Fo-Fc density seen for Ala298 from 2PB0 and from 2PB2 (all four monomers show equivalent density). It is clear that this residue should be modelled as a valine, as the peaks for the Cc1 and Cc2 are 10-11 sigma in the difference maps. (TIF) Figure S6 Plots of docked conformations using different scoring algorithms. Top: Box-whisker plot of relative scores for correct (sidechain amine near catalytic nitrogen) versus backwards (back-bone amine near catalytic nitrogen) docked conformations using three different scoring functions: Shapegauss, PLP, and Chemgauss3. Middle: Box-whisker plots for the relative scores of docked conformations separated by substrate and by scoring function. Bottom: Box-whisker plot for the relative scores of docked conformations for different sub-terms of the Chemgauss3 scoring function. (TIF) Figure S7 The electrostatic potential of the AstC binding pocket. Shown here is the binding pocket of succinylornithine transaminase with the electrostatic potential of the apo pocket mapped to the molecular surface (colored from 25k B T in red to 5k B T in blue). The crystallographic ligand succinalglutamic semialdehyde is shown in cyan. The presence of large patches with highly positive potential can readily be seen.