Thermotoga neapolitana adenylate kinase is highly active at 30 degrees C

Protein tolerance to random circular permutation correlates with thermostability and local energetics of residue-residue contacts

Adenylate kinase (AK) orthologs with a range of thermostabilities were subjected to random circular permutation, and deep mutational scanning was used to evaluate where new protein termini were nondisruptive to activity. The fraction of circularly permuted variants that retained function in each library correlated with AK thermostability. In addition, analysis of the positional tolerance to new termini, which increase local conformational flexibility, showed that bonds were either functionally sensitive to cleavage across all homologs, differentially sensitive, or uniformly tolerant. The mobile AMP-binding domain, which displays the highest calculated contact energies, presented the greatest tolerance to new termini across all AKs. In contrast, retention of function in the lid and core domains was more dependent upon AK melting temperature. These results show that family permutation profiling identifies primary structure that has been selected by evolution for dynamics that are critical to activity within an enzyme family. These findings also illustrate how deep mutational scanning can be applied to protein homologs in parallel to differentiate how topology, stability, and local energetics govern mutational tolerance.

Engineering a Seven Enzyme Biotransformation using Mathematical Modelling and Characterized Enzyme Parts

Multi-step enzyme reactions offer considerable cost and productivity benefits. Process models offer a route to understanding the complexity of these reactions, and allow for their optimization. Despite the increasing prevalence of multi-step biotransformations, there are few examples of process models for enzyme reactions. From a toolbox of characterized enzyme parts, we demonstrate the construction of a process model for a seven enzyme, three step biotransformation using isolated enzymes. Enzymes for cofactor regeneration were employed to make this in vitro reaction economical. Good modelling practice was critical in evaluating the impact of approximations and experimental error. We show that the use and validation of process models was instrumental in realizing and removing process bottlenecks, identifying divergent behavior, and for the optimization of the entire reaction using a genetic algorithm. We validated the optimized reaction to demonstrate that complex multi-step reactions with cofactor recycling involving at least seven enzymes can be reliably modelled and optimized.

Circular permutation profiling by deep sequencing libraries created using transposon mutagenesis

Deep mutational scanning has been used to create high-resolution DNA sequence maps that illustrate the functional consequences of large numbers of point mutations. However, this approach has not yet been applied to libraries of genes created by random circular permutation, an engineering strategy that is used to create open reading frames that express proteins with altered contact order. We describe a new method, termed circular permutation profiling with DNA sequencing (CPP-seq), which combines a one-step transposon mutagenesis protocol for creating libraries with a functional selection, deep sequencing and computational analysis to obtain unbiased insight into a protein's tolerance to circular permutation. Application of this method to an adenylate kinase revealed that CPP-seq creates two types of vectors encoding each circularly permuted gene, which differ in their ability to express proteins. Functional selection of this library revealed that >65% of the sampled vectors that express proteins are enriched relative to those that cannot translate proteins. Mapping enriched sequences onto structure revealed that the mobile AMP binding and rigid core domains display greater tolerance to backbone fragmentation than the mobile lid domain, illustrating how CPP-seq can be used to relate a protein's biophysical characteristics to the retention of activity upon permutation.

The Structure of a Thermophilic Kinase Shapes Fitness upon Random Circular Permutation

Proteins can be engineered for synthetic biology through circular permutation, a sequence rearrangement in which native protein termini become linked and new termini are created elsewhere through backbone fission. However, it remains challenging to anticipate a protein's functional tolerance to circular permutation. Here, we describe new transposons for creating libraries of randomly circularly permuted proteins that minimize peptide additions at their termini, and we use transposase mutagenesis to study the tolerance of a thermophilic adenylate kinase (AK) to circular permutation. We find that libraries expressing permuted AKs with either short or long peptides amended to their N-terminus yield distinct sets of active variants and present evidence that this trend arises because permuted protein expression varies across libraries. Mapping all sites that tolerate backbone cleavage onto AK structure reveals that the largest contiguous regions of sequence that lack cleavage sites are proximal to the phosphotransfer site. A comparison of our results with a range of structure-derived parameters further showed that retention of function correlates to the strongest extent with the distance to the phosphotransfer site, amino acid variability in an AK family sequence alignment, and residue-level deviations in superimposed AK structures. Our work illustrates how permuted protein libraries can be created with minimal peptide additions using transposase mutagenesis, and it reveals a challenge of maintaining consistent expression across permuted variants in a library that minimizes peptide additions. Furthermore, these findings provide a basis for interpreting responses of thermophilic phosphotransferases to circular permutation by calibrating how different structure-derived parameters relate to retention of function in a cellular selection.

A transposase strategy for creating libraries of circularly permuted proteins

A simple approach for creating libraries of circularly permuted proteins is described that is called PERMutation Using Transposase Engineering (PERMUTE). In PERMUTE, the transposase MuA is used to randomly insert a minitransposon that can function as a protein expression vector into a plasmid that contains the open reading frame (ORF) being permuted. A library of vectors that express different permuted variants of the ORF-encoded protein is created by: (i) using bacteria to select for target vectors that acquire an integrated minitransposon; (ii) excising the ensemble of ORFs that contain an integrated minitransposon from the selected vectors; and (iii) circularizing the ensemble of ORFs containing integrated minitransposons using intramolecular ligation. Construction of a Thermotoga neapolitana adenylate kinase (AK) library using PERMUTE revealed that this approach produces vectors that express circularly permuted proteins with distinct sequence diversity from existing methods. In addition, selection of this library for variants that complement the growth of Escherichia coli with a temperature-sensitive AK identified functional proteins with novel architectures, suggesting that PERMUTE will be useful for the directed evolution of proteins with new functions.

Mesophilic and hyperthermophilic adenylate kinases differ in their tolerance to random fragmentation

The extent to which thermostability influences the location of protein fragmentation sites that allow retention of function is not known. To evaluate this, we used a novel transposase-based approach to create libraries of vectors that express structurally-related fragments of Bacillus subtilis adenylate kinase (BsAK) and Thermotoga neapolitana adenylate kinase (TnAK) with identical modifications at their termini, and we selected for variants in each library that complement the growth of Escherichia coli with a temperature-sensitive adenylate kinase (AK). Mutants created using the hyperthermophilic TnAK were found to support growth with a higher frequency (44%) than those generated from the mesophilic BsAK (6%), and selected TnAK mutants complemented E. coli growth more strongly than homologous BsAK variants. Sequencing of functional clones from each library also identified a greater dispersion of fragmentation sites within TnAK. Nondisruptive fission sites were observed within the AMP binding and core domains of both AK homologs. However, only TnAK contained sites within the lid domain, which undergoes dynamic fluctuations that are critical for catalysis. These findings implicate the flexible lid domain as having an increased sensitivity to fission events at physiological temperatures. In addition, they provide evidence that comparisons of nondisruptive fission sites in homologous proteins could be useful for finding dynamic regions whose conformational fluctuations are important for function, and they show that the discovery of protein fragments that cooperatively function in mesophiles can be aided by the use of thermophilic enzymes as starting points for protein design.

Crystal structure of the zinc-, cobalt-, and iron-containing adenylate kinase from Desulfovibrio gigas: a novel metal-containing adenylate kinase from Gram-negative bacteria

Adenylate kinases (AK) from Gram-negative bacteria are generally devoid of metal ions in their LID domain. However, three metal ions, zinc, cobalt, and iron, have been found in AK from Gram-negative bacteria. Crystal structures of substrate-free AK from Desulfovibrio gigas with three different metal ions (Zn(2+), Zn-AK; Co(2+), Co-AK; and Fe(2+), Fe-AK) bound in its LID domain have been determined by X-ray crystallography to resolutions 1.8, 2.0, and 3.0 Å, respectively. The zinc and iron forms of the enzyme were crystallized in space group I222, whereas the cobalt-form crystals were C2. The presence of the metals was confirmed by calculation of anomalous difference maps and by X-ray fluorescence scans. The work presented here is the first report of a structure of a metal-containing AK from a Gram-negative bacterium. The native enzyme was crystallized, and only zinc was detected in the LID domain. Co-AK and Fe-AK were obtained by overexpressing the protein in Escherichia coli. Zn-AK and Fe-AK crystallized as monomers in the asymmetric unit, whereas Co-AK crystallized as a dimer. Nevertheless, all three crystal structures are very similar to each other, with the same LID domain topology, the only change being the presence of the different metal atoms. In the absence of any substrate, the LID domain of all holoforms of AK was present in a fully open conformational state. Normal mode analysis was performed to predict fluctuations of the LID domain along the catalytic pathway.

A selection that reports on protein-protein interactions within a thermophilic bacterium

Many proteins can be split into fragments that exhibit enhanced function upon fusion to interacting proteins. While this strategy has been widely used to create protein-fragment complementation assays (PCAs) for discovering protein-protein interactions within mesophilic organisms, similar assays have not yet been developed for studying natural and engineered protein complexes at the temperatures where thermophilic microbes grow. We describe the development of a selection for protein-protein interactions within Thermus thermophilus that is based upon growth complementation by fragments of Thermotoga neapolitana adenylate kinase (AK(Tn)). Complementation studies with an engineered thermophile (PQN1) that is not viable above 75 degrees C because its adk gene has been replaced by a Geobacillus stearothermophilus ortholog revealed that growth could be restored at 78 degrees C by a vector that coexpresses polypeptides corresponding to residues 1-79 and 80-220 of AK(Tn). In contrast, PQN1 growth was not complemented by AK(Tn) fragments harboring a C156A mutation within the zinc-binding tetracysteine motif unless these fragments were fused to Thermotoga maritima chemotaxis proteins that heterodimerize (CheA and CheY) or homodimerize (CheX). This enhanced complementation is interpreted as arising from chemotaxis protein-protein interactions, since AK(Tn)-C156A fragments having only one polypeptide fused to a chemotaxis protein did not complement PQN1 to the same extent. This selection increases the maximum temperature where a PCA can be used to engineer thermostable protein complexes and to map protein-protein interactions.

Dynamics in Thermotoga neapolitana adenylate kinase: 15N relaxation and hydrogen-deuterium exchange studies of a hyperthermophilic enzyme highly active at 30 degrees C

Backbone conformational dynamics of Thermotoga neapolitana adenylate kinase in the free form (TNAK) and inhibitor-bound form (TNAK*Ap5A) were investigated at 30 degrees C using (15)N NMR relaxation measurements and NMR monitored hydrogen-deuterium exchange. With kinetic parameters identical to those of Escherichia coli AK (ECAK) at 30 degrees C, TNAK is a unique hyperthermophilic enzyme. These catalytic properties make TNAK an interesting and novel model to study the interplay between protein rigidity, stability, and activity. Comparison of fast time scale dynamics (picosecond to nanosecond) in the open and closed states of TNAK and ECAK at 30 degrees C reveals a uniformly higher rigidity across all domains of TNAK. Within this framework of a rigid TNAK structure, several residues located in the AMP-binding domain and in the core-lid hinge regions display high picosecond to nanosecond time scale flexibility. Together with the recent comparison of ECAK dynamics with those of hyperthermophilic Aquifex aeolicus AK (AAAK), our results provide strong evidence for the role of picosecond to nanosecond time scale fluctuations in both stability and activity. In the slow time scales, TNAK's increased rigidity is not uniform but localized in the AMP-binding and lid domains. The core domain amides of ECAK and TNAK in the open and closed states show comparable protection against exchange. Significantly, the hinges framing the lid domain show similar exchange data in ECAK and TNAK open and closed forms. Our NMR relaxation and hydrogen-deuterium exchange studies therefore suggest that TNAK maintains high activity at 30 degrees C by localizing flexibility to the hinge regions that are key to facilitating conformational changes.

Cobalt-, zinc- and iron-bound forms of adenylate kinase (AK) from the sulfate-reducing bacterium Desulfovibrio gigas: purification, crystallization and preliminary X-ray diffraction analysis

Adenylate kinase (AK; ATP:AMP phosphotransferase; EC 2.7.4.3) is involved in the reversible transfer of the terminal phosphate group from ATP to AMP. AKs contribute to the maintenance of a constant level of cellular adenine nucleotides, which is necessary for the energetic metabolism of the cell. Three metal ions, cobalt, zinc and iron(II), have been reported to be present in AKs from some Gram-negative bacteria. Native zinc-containing AK from Desulfovibrio gigas was purified to homogeneity and crystallized. The crystals diffracted to beyond 1.8 A resolution. Furthermore, cobalt- and iron-containing crystal forms of recombinant AK were also obtained and diffracted to 2.0 and 3.0 A resolution, respectively. Zn(2+)-AK and Fe(2+)-AK crystallized in space group I222 with similar unit-cell parameters, whereas Co(2+)-AK crystallized in space group C2; a monomer was present in the asymmetric unit for both the Zn(2+)-AK and Fe(2+)-AK forms and a dimer was present for the Co(2+)-AK form. The structures of the three metal-bound forms of AK will provide new insights into the role and selectivity of the metal in these enzymes.

Structure and biochemical characterization of an adenylate kinase originating from the psychrophilic organism Marinibacillus marinus

Adenylate kinases (AKs; EC 2.7.4.3) are essential members of the NMP kinase family that maintain cellular homeostasis by the interconversion of AMP, ADP and ATP. AKs play a critical role in adenylate homeostasis across all domains of life and have been used extensively as prototypes for the study of protein adaptation and the relationship of protein dynamics and stability to function. To date, kinetic studies of psychrophilic AKs have not been performed. In order to broaden understanding of extremophilic adaptation, the kinetic parameters of adenylate kinase from the psychrophile Marinibacillus marinus were examined and the crystal structure of this cold-adapted enzyme was determined at 2.0 A resolution. As expected, the overall structure and topology of the psychrophilic M. marinus AK are similar to those of mesophilic and thermophilic AKs. The thermal denaturation midpoint of M. marinus AK (321.1 K) is much closer to that of the mesophile Bacillus subtilis (320.7 K) than the more closely related psychrophile B. globisporus (316.4 K). In addition, the enzymatic properties of M. marinus AK are quite close to those of the mesophilic AK and suggests that M. marinus experiences temperature ranges in which excellent enzyme function over a broad temperature range (293-313 K) has been retained for the success of the organism. Even transient loss of AK function is lethal and as a consequence AK must be robust and be well adapted to the environment of the host organism.

A new type of metal-binding site in cobalt- and zinc-containing adenylate kinases isolated from sulfate-reducers Desulfovibrio gigas and Desulfovibrio desulfuricans ATCC 27774

Adenylate kinase (AK) mediates the reversible transfer of phosphate groups between the adenylate nucleotides and contributes to the maintenance of their constant cellular level, necessary for energy metabolism and nucleic acid synthesis. The AK were purified from crude extracts of two sulfate-reducing bacteria (SRB), Desulfovibrio (D.) gigas NCIB 9332 and Desulfovibrio desulfuricans ATCC 27774, and biochemically and spectroscopically characterised in the native and fully cobalt- or zinc-substituted forms. These are the first reported adenylate kinases that bind either zinc or cobalt and are related to the subgroup of metal-containing AK found, in most cases, in Gram-positive bacteria. The electronic absorption spectrum is consistent with tetrahedral coordinated cobalt, predominantly via sulfur ligands, and is supported by EPR. The involvement of three cysteines in cobalt or zinc coordination was confirmed by chemical methods. Extended X-ray absorption fine structure (EXAFS) indicate that cobalt or zinc are bound by three cysteine residues and one histidine in the metal-binding site of the "LID" domain. The sequence 129Cys-X5-His-X15-Cys-X2-Cys of the AK from D. gigas is involved in metal coordination and represents a new type of binding motif that differs from other known zinc-binding sites of AK. Cobalt and zinc play a structural role in stabilizing the LID domain.

Associative mechanism for phosphoryl transfer: a molecular dynamics simulation of Escherichia coli adenylate kinase complexed with its substrates

The ternary complex of Escherichia coli adenylate kinase (ECAK) with its substrates adenosine monophosphate (AMP) and Mg-ATP, which catalyzes the reversible transfer of a phosphoryl group between adenosine triphosphate (ATP) and AMP, was studied using molecular dynamics. The starting structure for the simulation was assembled from the crystal structures of ECAK complexed with the bisubstrate analog diadenosine pentaphosphate (AP(5)A) and of Bacillus stearothermophilus adenylate kinase complexed with AP(5)A, Mg(2+), and 4 coordinated water molecules, and by deleting 1 phosphate group from AP(5)A. The interactions of ECAK residues with the various moieties of ATP and AMP were compared to those inferred from NMR, X-ray crystallography, site-directed mutagenesis, and enzyme kinetic studies. The simulation supports the hypothesis that hydrogen bonds between AMP's adenine and the protein are at the origin of the high nucleoside monophosphate (NMP) specificity of AK. The ATP adenine and ribose moieties are only loosely bound to the protein, while the ATP phosphates are strongly bound to surrounding residues. The coordination sphere of Mg(2+), consisting of 4 waters and oxygens of the ATP beta- and gamma-phosphates, stays approximately octahedral during the simulation. The important role of the conserved Lys13 in the P loop in stabilizing the active site by bridging the ATP and AMP phosphates is evident. The influence of Mg(2+), of its coordination waters, and of surrounding charged residues in maintaining the geometry and distances of the AMP alpha-phosphate and ATP beta- and gamma-phosphates is sufficient to support an associative reaction mechanism for phosphoryl transfer.

A GTP:AMP phosphotransferase, Adk2p, in Saccharomyces cerevisiae. Role of the C terminus in protein folding/stabilization, thermal tolerance, and enzymatic activity

Adenylate kinases participate in maintaining the homeostasis of cellular nucleotides. Depending on the yeast strains, the GTP:AMP phosphotransferase is encoded by the nuclear gene ADK2 with or without a single base pair deletion/insertion near the 3' end of the open reading frame, and the corresponding protein exists as either Adk2p (short) or Adk2p (long) in the mitochondrial matrix. These two forms are identical except that the three C-terminal residues of Adk2p (short) are changed in Adk2p (long), and the latter contains an additional nine amino acids at the C terminus of the protein. The short form of Adk2p has so far been considered to be inactive (Schricker, R., Magdolen, V., Strobel, G., Bogengruber, E., Breitenbach, M., and Bandlow, W. (1995) J. Biol. Chem. 270, 31103-31110). Using purified proteins, we show that at the physiological temperature for yeast growth (30 degrees C), both short and long forms of Adk2p are enzymatically active. However, in contrast to the short form, Adk2p (long) is quite resistant to thermal inactivation, urea denaturation, and degradation by trypsin. Unfolding of the long form by high concentrations of urea greatly stimulated its import into isolated mitochondria. Using an integration-based gene-swapping approach, we found that regardless of the yeast strains used, the steady state levels of endogenous Adk2p (long) in mitochondria were 5-10-fold lower compared with those of Adk2p (short). Together, these results suggest that the modified C-terminal domain in Adk2p (long) is not essential for enzyme activity, but it contributes to and strengthens protein folding and/or stability and is particularly important for maintaining enzyme activity under stress conditions.