Evolutionary drivers of thermoadaptation in enzyme catalysis

Parallel molecular mechanisms for enzyme temperature adaptation

The mechanisms that underly the adaptation of enzyme activities and stabilities to temperature are fundamental to our understanding of molecular evolution and how enzymes work. Here, we investigate the molecular and evolutionary mechanisms of enzyme temperature adaption, combining deep mechanistic studies with comprehensive sequence analyses of thousands of enzymes. We show that temperature adaptation in ketosteroid isomerase (KSI) arises primarily from one residue change with limited, local epistasis, and we establish the underlying physical mechanisms. This residue change occurs in diverse KSI backgrounds, suggesting parallel adaptation to temperature. We identify residues associated with organismal growth temperature across 1005 diverse bacterial enzyme families, suggesting widespread parallel adaptation to temperature. We assess the residue properties, molecular interactions, and interaction networks that appear to underly temperature adaptation.

Hinge-shift mechanism as a protein design principle for the evolution of β-lactamases from substrate promiscuity to specificity

TEM-1 β-lactamase degrades β-lactam antibiotics with a strong preference for penicillins. Sequence reconstruction studies indicate that it evolved from ancestral enzymes that degraded a variety of β-lactam antibiotics with moderate efficiency. This generalist to specialist conversion involved more than 100 mutational changes, but conserved fold and catalytic residues, suggesting a role for dynamics in enzyme evolution. Here, we develop a conformational dynamics computational approach to rationally mold a protein flexibility profile on the basis of a hinge-shift mechanism. By deliberately weighting and altering the conformational dynamics of a putative Precambrian β-lactamase, we engineer enzyme specificity that mimics the modern TEM-1 β-lactamase with only 21 amino acid replacements. Our conformational dynamics design thus re-enacts the evolutionary process and provides a rational allosteric approach for manipulating function while conserving the enzyme active site.

Ancestral lysosomal enzymes with increased activity harbor therapeutic potential for treatment of Hunter syndrome

We show the successful application of ancestral sequence reconstruction to enhance the activity of iduronate-2-sulfatase (IDS), thereby increasing its therapeutic potential for the treatment of Hunter syndrome—a lysosomal storage disease caused by impaired function of IDS. Current treatment, enzyme replacement therapy with recombinant human IDS, does not alleviate all symptoms, and an unmet medical need remains. We reconstructed putative ancestral sequences of mammalian IDS and compared them with extant IDS. Some ancestral variants displayed up to 2-fold higher activity than human IDS in in vitro assays and cleared more substrate in ex vivo experiments in patient fibroblasts. This could potentially allow for lower dosage or enhanced therapeutic effect in enzyme replacement therapy, thereby improving treatment outcomes and cost efficiency, as well as reducing treatment burden. In summary, we showed that ancestral sequence reconstruction can be applied to lysosomal enzymes that function in concert with modern enzymes and receptors in cells.

•

Reconstruction of ancestral lysosomal enzymes that function in complex cellular context

•

Ancestral iduronate-2-sulfatases with increased activity compared with the human enzyme

•

Increased clearance of substrate in patient fibroblasts indicates therapeutic potential

Biomolecular QM/MM Simulations: What Are Some of the "Burning Issues"?

QM/MM simulations have become an indispensable tool in many chemical and biochemical investigations. Considering the tremendous degree of success, including recognition by a 2013 Nobel Prize in Chemistry, are there still "burning challenges" in QM/MM methods, especially for biomolecular systems? In this short Perspective, we discuss several issues that we believe greatly impact the robustness and quantitative applicability of QM/MM simulations to many, if not all, biomolecules. We highlight these issues with observations and relevant advances from recent studies in our group and others in the field. Despite such limited scope, we hope the discussions are of general interest and will stimulate additional developments that help push the field forward in meaningful directions.

Computational design of enzymes for biotechnological applications

Enzymes are the natural catalysts that execute biochemical reactions upholding life. Their natural effectiveness has been fine-tuned as a result of millions of years of natural evolution. Such catalytic effectiveness has prompted the use of biocatalysts from multiple sources on different applications, including the industrial production of goods (food and beverages, detergents, textile, and pharmaceutics), environmental protection, and biomedical applications. Natural enzymes often need to be improved by protein engineering to optimize their function in non-native environments. Recent technological advances have greatly facilitated this process by providing the experimental approaches of directed evolution or by enabling computer-assisted applications. Directed evolution mimics the natural selection process in a highly accelerated fashion at the expense of arduous laboratory work and economic resources. Theoretical methods provide predictions and represent an attractive complement to such experiments by waiving their inherent costs. Computational techniques can be used to engineer enzymatic reactivity, substrate specificity and ligand binding, access pathways and ligand transport, and global properties like protein stability, solubility, and flexibility. Theoretical approaches can also identify hotspots on the protein sequence for mutagenesis and predict suitable alternatives for selected positions with expected outcomes. This review covers the latest advances in computational methods for enzyme engineering and presents many successful case studies.

Structure, Function, and Thermal Adaptation of the Biotin Carboxylase Domain Dimer from Hydrogenobacter thermophilus 2-Oxoglutarate Carboxylase

2-Oxoglutarate carboxylase (OGC), a unique member of the biotin-dependent carboxylase family from the order Aquificales, captures dissolved CO₂ via the reductive tricarboxylic acid (rTCA) cycle. Structure and function studies of OGC may facilitate adaptation of the rTCA cycle to increase the level of carbon fixation for biofuel production. Here we compare the biotin carboxylase (BC) domain of Hydrogenobacter thermophilus OGC with the well-studied mesophilic homologues to identify features that may contribute to thermal stability and activity. We report three OGC BC X-ray structures, each bound to bicarbonate, ADP, or ADP-Mg²⁺, and propose that substrate binding at high temperatures is facilitated by interactions that stabilize the flexible subdomain B in a partially closed conformation. Kinetic measurements with varying ATP and biotin concentrations distinguish two temperature-dependent steps, consistent with biotin's rate-limiting role in organizing the active site. Transition state thermodynamic values derived from the Eyring equation indicate a larger positive ΔH^⧧ and a less negative ΔS^⧧ compared to those of a previously reported mesophilic homologue. These thermodynamic values are explained by partially rate limiting product release. Phylogenetic analysis of BC domains suggests that OGC diverged prior to Aquificales evolution. The phylogenetic tree identifies mis-annotations of the Aquificales BC sequences, including the Aquifex aeolicus pyruvate carboxylase structure. Notably, our structural data reveal that the OGC BC dimer comprises a "wet" dimerization interface that is dominated by hydrophilic interactions and structural water molecules common to all BC domains and likely facilitates the conformational changes associated with the catalytic cycle. Mutations in the dimerization domain demonstrate that dimerization contributes to thermal stability.

Bayesian genome scale modelling identifies thermal determinants of yeast metabolism

The molecular basis of how temperature affects cell metabolism has been a long-standing question in biology, where the main obstacles are the lack of high-quality data and methods to associate temperature effects on the function of individual proteins as well as to combine them at a systems level. Here we develop and apply a Bayesian modeling approach to resolve the temperature effects in genome scale metabolic models (GEM). The approach minimizes uncertainties in enzymatic thermal parameters and greatly improves the predictive strength of the GEMs. The resulting temperature constrained yeast GEM uncovers enzymes that limit growth at superoptimal temperatures, and squalene epoxidase (ERG1) is predicted to be the most rate limiting. By replacing this single key enzyme with an ortholog from a thermotolerant yeast strain, we obtain a thermotolerant strain that outgrows the wild type, demonstrating the critical role of sterol metabolism in yeast thermosensitivity. Therefore, apart from identifying thermal determinants of cell metabolism and enabling the design of thermotolerant strains, our Bayesian GEM approach facilitates modelling of complex biological systems in the absence of high-quality data and therefore shows promise for becoming a standard tool for genome scale modeling.

While temperature impacts the function of all cellular components, it’s hard to rule out how the temperature dependence of cell phenotypes emerged from the dependence of individual components. Here, the authors develop a Bayesian genome scale modelling approach to identify thermal determinants of yeast metabolism.

The high-energy transition state of the glutamate transporter homologue GltPh

Membrane transporters mediate cellular uptake of nutrients, signaling molecules, and drugs. Their overall mechanisms are often well understood, but the structural features setting their rates are mostly unknown. Earlier single-molecule fluorescence imaging of the archaeal model glutamate transporter homologue GltPh from Pyrococcus horikoshii suggested that the slow conformational transition from the outward- to the inward-facing state, when the bound substrate is translocated from the extracellular to the cytoplasmic side of the membrane, is rate limiting to transport. Here, we provide insight into the structure of the high-energy transition state of GltPh that limits the rate of the substrate translocation process. Using bioinformatics, we identified GltPh gain-of-function mutations in the flexible helical hairpin domain HP2 and applied linear free energy relationship analysis to infer that the transition state structurally resembles the inward-facing conformation. Based on these analyses, we propose an approach to search for allosteric modulators for transporters.

Evolution of Protein Structure and Stability in Global Warming

This review focuses on the molecular signatures of protein structures in relation to evolution and survival in global warming. It is based on the premise that the power of evolutionary selection may lead to thermotolerant organisms that will repopulate the planet and continue life in general, but perhaps with different kinds of flora and fauna. Our focus is on molecular mechanisms, whereby known examples of thermoresistance and their physicochemical characteristics were noted. A comparison of interactions of diverse residues in proteins from thermophilic and mesophilic organisms, as well as reverse genetic studies, revealed a set of imprecise molecular signatures that pointed to major roles of hydrophobicity, solvent accessibility, disulfide bonds, hydrogen bonds, ionic and π-electron interactions, and an overall condensed packing of the higher-order structure, especially in the hydrophobic regions. Regardless of mutations, specialized protein chaperones may play a cardinal role. In evolutionary terms, thermoresistance to global warming will likely occur in stepwise mutational changes, conforming to the molecular signatures, such that each "intermediate" fits a temporary niche through punctuated equilibrium, while maintaining protein functionality. Finally, the population response of different species to global warming may vary substantially, and, as such, some may evolve while others will undergo catastrophic mass extinction.

Heat induces end to end repetitive association in P. furiosus L-asparaginase which enables its thermophilic property

It remains undeciphered how thermophilic enzymes display enhanced stability at elevated temperatures. Taking l-asparaginase from P. furiosus (PfA) as an example, we combined scattering shapes deduced from small-angle X-ray scattering (SAXS) data at increased temperatures with symmetry mates from crystallographic structures to find that heating caused end-to-end association. The small contact point of self-binding appeared to be enabled by a terminal short β-strand in N-terminal domain, Leu¹⁷⁹-Val-Val-Asn¹⁸² (LVVN). Interestingly, deletion of this strand led to a defunct enzyme, whereas suplementation of the peptide LVVN to the defunct enzyme restored structural frameworkwith mesophile-type functionality. Crystal structure of the peptide-bound defunct enzyme showed that one peptide ispresent in the same coordinates as in original enzyme, explaining gain-of lost function. A second peptide was seen bound to the protein at a different location suggesting its possible role in substrate-free molecular-association. Overall, we show that the heating induced self-assembly of native shapes of PfA led to an apparent super-stable assembly.

Hidden Conformational States and Strange Temperature Optima in Enzyme Catalysis

The existence of temperature optima in enzyme catalysis that occur before protein melting sets in can be described by different types of kinetic models. Such optima cause distinctly curved Arrhenius plots and have, for example, been observed in several cold-adapted enzymes from psychrophilic species. The two main explanations proposed for this behavior either invoke conformational equilibria with inactive substrate-bound states or postulate differences in heat capacity between the reactant and transition states. Herein, we analyze the implications of the different types of kinetic models in terms of apparent activation enthalpies, entropies, and heat capacities, using the catalytic reaction of a cold-adapted α-amylase as a prototypic example. We show that the behavior of these thermodynamic activation parameters is fundamentally different between equilibrium and heat capacity models, and in the α-amylase case, computer simulations have shown the former model to be correct. A few other enzyme-catalyzed reactions are also discussed in this context.

Ancestral sequence reconstruction produces thermally stable enzymes with mesophilic enzyme-like catalytic properties

Enzymes have high catalytic efficiency and low environmental impact, and are therefore potentially useful tools for various industrial processes. Crucially, however, natural enzymes do not always have the properties required for specific processes. It may be necessary, therefore, to design, engineer, and evolve enzymes with properties that are not found in natural enzymes. In particular, the creation of enzymes that are thermally stable and catalytically active at low temperature is desirable for processes involving both high and low temperatures. In the current study, we designed two ancestral sequences of 3-isopropylmalate dehydrogenase by an ancestral sequence reconstruction technique based on a phylogenetic analysis of extant homologous amino acid sequences. Genes encoding the designed sequences were artificially synthesized and expressed in Escherichia coli. The reconstructed enzymes were found to be slightly more thermally stable than the extant thermophilic homologue from Thermus thermophilus. Moreover, they had considerably higher low-temperature catalytic activity as compared with the T. thermophilus enzyme. Detailed analyses of their temperature-dependent specific activities and kinetic properties showed that the reconstructed enzymes have catalytic properties similar to those of mesophilic homologues. Collectively, our study demonstrates that ancestral sequence reconstruction can produce a thermally stable enzyme with catalytic properties adapted to low-temperature reactions.

Quantifying the Relationship between Conformational Dynamics and Enzymatic Activity in Ribonuclease HI Homologues

Ribonuclease HI (RNHI), a ubiquitous, non-sequence-specific endonuclease, cleaves the RNA strand in RNA/DNA hybrids. RNHI functions in replication and genome maintenance, and retroviral reverse transcriptases contain an essential ribonuclease H domain. Nuclear magnetic resonance (NMR) spectroscopy combined with molecular dynamics (MD) simulations suggests a model in which the extended handle region domain of Escherichia coli RNHI populates (substrate-binding-competent) "open" and (substrate-binding-incompetent) "closed" states, while the thermophilic Thermus thermophilus RNHI mainly populates the closed state at 300 K [Stafford, K. A., Robustelli, P., and Palmer, A. G., III (2013) PLoS Comput. Biol. 9, 1-10]. In addition, an in silico-designed mutant E. coli Val98Ala RNHI was predicted to populate primarily the closed state. The work presented here validates this model and confirms the predicted properties of the designed mutant. MD simulations suggest that the conformational preferences of the handle region correlate with the conformations of Trp85, Thr92, and Val101. NMR residual dipolar coupling constants, three-bond scalar coupling constants, and chemical shifts experimentally define the conformational states of these residues and hence of the handle domain. These NMR parameters correlate with the Michaelis constants for RNHI homologues, confirming the important role of the handle region in the modulation of substrate recognition and illustrating the power of NMR spectroscopy in dissecting the conformational preferences underlying enzyme function.

Empowering Global Chemical Biology at the Dawn of the New Decade

On January 22-24, 2020, scientific luminaries across the far-flung corners of chemical biology gathered in Geneva, Switzerland, to deliver their latest and greatest discoveries in the field. Generously supported by the Swiss National Science Foundation (SNSF), our academic partners, and industrial and journal sponsors, this chemical biology symposium in our opinion will remain memorable for several years to come, not only because of the diversity in scientific topics delivered by our invited eminent speakers as detailed herein, but it is also one-of-a-kind conference which reflected multidimensional balance-balance in age and gender, across these speakers. Such a remarkable speaker line-up doubtless attracted >200 attendees from academia and industry in and around Switzerland and beyond, representing a huge swathe of subfields of science interfacing chemistry and biology. Poster presentations from students and postdocs further spotlighted the exciting diversity in the field: spanning biosynthesis, optochemical genetics, genetic code expansion, lipid chemical biology, redox perturbation, microfluidics screening, membrane signaling, immune modulation, DNA circuits, and synthetic and computational biology. This notable heterogeneity in scientific topics also went hand-in-hand with the diverse representations of student/postdoc trainees from 56 institutions covering 14 countries worldwide, allowing us to witness science as a truly global enterprise.

Predicting the optimal growth temperatures of prokaryotes using only genome derived features

MOTIVATION

Optimal growth temperature is a fundamental characteristic of all living organisms. Knowledge of this temperature is central to the study of a prokaryote, the thermal stability and temperature dependent activity of its genes, and the bioprospecting of its genome for thermally adapted proteins. While high throughput sequencing methods have dramatically increased the availability of genomic information, the growth temperatures of the source organisms are often unknown. This limits the study and technological application of these species and their genomes. Here, we present a novel method for the prediction of growth temperatures of prokaryotes using only genomic sequences.

RESULTS

By applying the reverse ecology principle that an organism's genome includes identifiable adaptations to its native environment, we can predict a species' optimal growth temperature with an accuracy of 5.17°C root-mean-square error and a coefficient of determination of 0.835. The accuracy can be further improved for specific taxonomic clades or by excluding psychrophiles. This method provides a valuable tool for the rapid calculation of organism growth temperature when only the genome sequence is known.

AVAILABILITY AND IMPLEMENTATION

Source code, genomes analyzed and features calculated are available at: https://github.com/DavidBSauer/OGT_prediction.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Embedding Functional Biomacromolecules within Peptide-Directed Metal-Organic Framework (MOF) Nanoarchitectures Enables Activity Enhancement

Rationally tailoring a robust artificial coating can enhance the life-time of fragile biomacromolecules. However, the coating also can restrain the activity of the guest because of the decreased substrate accessibility. Herein, we report a peptide-directed strategy that enables in situ tailoring of the MOF-shrouded biohybrids into controllable nanoarchitectures. The MOF biohybrid can be shaped from different 3D microporous architectures into a 2D mesoporous layer by a peptide modulator. Using this mild strategy, we show that the nanoarchitectures of the MOF coatings significantly affect the biological functions of the contained biomacromolecules. The biomacromolecules entrapped within the novel 2D mesoporous spindle-shaped MOFs (2D MSMOFs) have significantly increased bioactivity compared to when encased within the hitherto explored 3D microporous MOFs. The improvement results from the shortened diffusion path and enlarged pore channel in 2D MSMOFs. Meanwhile, the thin 2D MSMOF layer also can provide excellent protection of the hosted biomacromolecules or protein-scaffolded biominerals through structural confinement.

Hybrid Thermophilic/Mesophilic Enzymes Reveal a Role for Conformational Disorder in Regulation of Bacterial Enzyme I

Conformational disorder is emerging as an important feature of biopolymers, regulating a vast array of cellular functions, including signaling, phase separation, and enzyme catalysis. Here we combine NMR, crystallography, computer simulations, protein engineering, and functional assays to investigate the role played by conformational heterogeneity in determining the activity of the C-terminal domain of bacterial Enzyme I (EIC). In particular, we design chimeric proteins by hybridizing EIC from thermophilic and mesophilic organisms, and we characterize the resulting constructs for structure, dynamics, and biological function. We show that EIC exists as a mixture of active and inactive conformations and that functional regulation is achieved by tuning the thermodynamic balance between active and inactive states. Interestingly, we also present a hybrid thermophilic/mesophilic enzyme that is thermostable and more active than the wild-type thermophilic enzyme, suggesting that hybridizing thermophilic and mesophilic proteins is a valid strategy to engineer thermostable enzymes with significant low-temperature activity.

Computer simulations explain the anomalous temperature optimum in a cold-adapted enzyme

Cold-adapted enzymes from psychrophilic species show the general characteristics of being more heat labile, and having a different balance between enthalpic and entropic contributions to free energy barrier of the catalyzed reaction compared to mesophilic orthologs. Among cold-adapted enzymes, there are also examples that show an enigmatic inactivation at higher temperatures before unfolding of the protein occurs. Here, we analyze these phenomena by extensive computer simulations of the catalytic reactions of psychrophilic and mesophilic α-amylases. The calculations yield temperature dependent reaction rates in good agreement with experiment, and also elicit the anomalous rate optimum for the cold-adapted enzyme, which occurs about 15 °C below the melting point. This result allows us to examine the structural basis of thermal inactivation, which turns out to be caused by breaking of a specific enzyme-substrate interaction. This type of behaviour is also likely to be relevant for other enzymes displaying such anomalous temperature optima.

Temperature, Dynamics, and Enzyme-Catalyzed Reaction Rates

We review the adaptations of enzyme activity to different temperatures. Psychrophilic (cold-adapted) enzymes show significantly different activation parameters (lower activation enthalpies and entropies) from their mesophilic counterparts. Furthermore, there is increasing evidence that the temperature dependence of many enzyme-catalyzed reactions is more complex than is widely believed. Many enzymes show curvature in plots of activity versus temperature that is not accounted for by denaturation or unfolding. This is explained by macromolecular rate theory: A negative activation heat capacity for the rate-limiting chemical step leads directly to predictions of temperature optima; both entropy and enthalpy are temperature dependent. Fluctuations in the transition state ensemble are reduced compared to the ground state. We show how investigations combining experiment with molecular simulation are revealing fundamental details of enzyme thermoadaptation that are relevant for understanding aspects of enzyme evolution. Simulations can calculate relevant thermodynamic properties (such as activation enthalpies, entropies, and heat capacities) and reveal the molecular mechanisms underlying experimentally observed behavior.

The Soybean Lipoxygenase-Substrate Complex: Correlation between the Properties of Tunneling-Ready States and ENDOR-Detected Structures of Ground States

Hydrogen tunneling in enzymatic C-H activation requires a dynamical sampling among ground-state enzyme-substrate (E-S) conformations, which transiently generates a tunneling-ready state (TRS). The TRS is characterized by a hydrogen donor-acceptor distance (DAD) of 2.7 Å, ∼0.5 Å shorter than the dominant DAD of optimized ground states. Recently, a high-resolution, 13C electron-nuclear double-resonance (ENDOR) approach was developed to characterize the ground-state structure of the complex of the linoleic acid (LA) substrate with soybean lipoxygenase (SLO). The resulting enzyme-substrate model revealed two ground-state conformers with different distances between the target C11 of LA and the catalytically active cofactor [Fe(III)-OH]: the active conformer "a", with a van der Waals DAD of 3.1 Å between C11 and metal-bound hydroxide, and an inactive conformer "b", with a distance that is almost 1 Å longer. Herein, the structure of the E-S complex is examined for a series of six variants in which subtle structural modifications of SLO have been introduced either at a hydrophobic side chain near the bound substrate or at a remote residue within a protein network whose flexibility influences hydrogen transfer. A remarkable correlation is found between the ENDOR-derived population of the active ground-state conformer a and the kinetically derived differential enthalpic barrier for D versus H transfer, ΔEa, with the latter increasing as the fraction of conformer a decreases. As proposed, ΔEa provides a "ruler" for the DAD within the TRS. ENDOR measurements further corroborate the previous identification of a dynamical network coupling the buried active site of SLO to the surface. This study shows that subtle imperfections within the initial ground-state structures of E-S complexes are accompanied by compromised geometries at the TRS.

Protein thermal stability

Proteins, in general, fold to a well-organized three-dimensional structure in order to function. The stability of this functional shape can be perturbed by external environmental conditions, such as temperature. Understanding the molecular factors underlying the resistance of proteins to the thermal stress has important consequences. First of all, it can aid the design of thermostable enzymes able to perform efficient catalysis in the high-temperature regime. Second, it is an essential brick of knowledge required to decipher the evolutionary pathways of life adaptation on Earth. Thanks to the development of atomistic simulations and ad hoc enhanced sampling techniques, it is now possible to investigate this problem in silico, and therefore provide support to experiments. After having described the methodological aspects, the chapter proposes an extended discussion on two problems. First, we focus on thermophilic proteins, a perfect model to address the issue of thermal stability and molecular evolution. Second, we discuss the issue of how protein thermal stability is affected by crowded in vivo-like conditions.

Thermostabilization of VPR, a kinetically stable cold adapted subtilase, via multiple proline substitutions into surface loops

Protein stability is a widely studied topic, there are still aspects however that need addressing. In this paper we examined the effects of multiple proline substitutions into loop regions of the kinetically stable proteinase K-like serine protease VPR, using the thermostable structural homologue AQUI as a template. Four locations for proline substitutions were chosen to imitate the structure of AQUI. Variants were produced and characterized using differential scanning calorimetry (DSC), circular dichroism (CD), steady state fluorescence, acrylamide fluorescence quenching and thermal inactivation experiments. The final product VPR_ΔC_N3P/I5P/N238P/T265P was greatly stabilized which was achieved without any noticeable detrimental effects to the catalytic efficiency of the enzyme. This stabilization seems to be derived from the conformation restrictive properties of the proline residue in its ability to act as an anchor point and strengthen pre-existing interactions within the protein and allowing for these interactions to prevail when thermal energy is applied to the system. In addition, the results underline the importance of the synergy between distant local protein motions needed to result in stabilizing effects and thus giving an insight into the nature of the stability of VPR, its unfolding landscape and how proline residues can infer kinetic stability onto protein structures.

Differences in thermal structural changes and melting between mesophilic and thermophilic dihydrofolate reductase enzymes

The thermal resistance of two homolog enzymes is investigated, with an emphasis on their local stability and flexibility, and on the possible implications regarding their reactivity.

Highly thermostable carboxylic acid reductases generated by ancestral sequence reconstruction

Carboxylic acid reductases (CARs) are biocatalysts of industrial importance. Their properties, especially their poor stability, render them sub-optimal for use in a bioindustrial pipeline. Here, we employed ancestral sequence reconstruction (ASR) - a burgeoning engineering tool that can identify stabilizing but enzymatically neutral mutations throughout a protein. We used a three-algorithm approach to reconstruct functional ancestors of the Mycobacterial and Nocardial CAR1 orthologues. Ancestral CARs (AncCARs) were confirmed to be CAR enzymes with a preference for aromatic carboxylic acids. Ancestors also showed varied tolerances to solvents, pH and in vivo-like salt concentrations. Compared to well-studied extant CARs, AncCARs had a Tm up to 35 °C higher, with half-lives up to nine times longer than the greatest previously observed. Using ancestral reconstruction we have expanded the existing CAR toolbox with three new thermostable CAR enzymes, providing access to the high temperature biosynthesis of aldehydes to drive new applications in biocatalysis.

Enhancing the thermostability and activity of uronate dehydrogenase from Agrobacterium tumefaciens LBA4404 by semi-rational engineering

Background

Glucaric acid, one of the aldaric acids, has been declared a “top value-added chemical from biomass”, and is especially important in the food and pharmaceutical industries. Biocatalytic production of glucaric acid from glucuronic acid is more environmentally friendly, efficient and economical than chemical synthesis. Uronate dehydrogenases (UDHs) are the key enzymes for the preparation of glucaric acid in this way, but the poor thermostability and low activity of UDH limit its industrial application. Therefore, improving the thermostability and activity of UDH, for example by semi-rational design, is a major research goal.

Results

In the present work, three UDHs were obtained from different Agrobacterium tumefaciens strains. The three UDHs have an approximate molecular weight of 32 kDa and all contain typically conserved UDH motifs. All three UDHs showed optimal activity within a pH range of 6.0–8.5 and at a temperature of 30 °C, but the UDH from A. tumefaciens (At) LBA4404 had a better catalytic efficiency than the other two UDHs (800 vs 600 and 530 s−1 mM−1). To further boost the catalytic performance of the UDH from AtLBA4404, site-directed mutagenesis based on semi-rational design was carried out. An A39P/H99Y/H234K triple mutant showed a 400-fold improvement in half-life at 59 °C, a 5 °C improvement in $$ {\text{T}}_{ 5 0}^{ 1 0} $$ T 50 10 value and a 2.5-fold improvement in specific activity at 30 °C compared to wild-type UDH.

Conclusions

In this study, we successfully obtained a triple mutant (A39P/H99Y/H234K) with simultaneously enhanced activity and thermostability, which provides a novel alternative for the industrial production of glucaric acid from glucuronic acid.

Targeted Heating of Enzyme Systems Based on Photothermal Materials

This study demonstrates that the enzymatic reaction rate can be increased significantly by targeted heating of the microenvironment around the enzyme, while maintaining the reaction system at environmental temperature. Enzyme molecules are covalently attached to the surface of Fe₃ O₄ @reduced graphite oxide (rGO). Under visible-light irradiation, the reaction rate catalyzed by the enzyme-Fe₃ O₄ @rGO system is clearly enhanced relative to that of the free enzyme and a mixture of free enzyme and Fe₃ O₄ @rGO. This local heating mechanism contributes to promotion of the enzymatic reactions of the targeted heating of the enzyme (THE) system, which has been validated by using different enzymes, including lipase, glucose oxidase, and organophosphorus hydrolase. These results indicate that targeted heating of the catalytic centers has the same effect on speeding up reactions as that of traditional heating methods, which treat the whole reaction system. As an example, it is shown that the THE system promotes the sensitivity of an enzyme screen-printed electrode by 14 times at room temperature, which implies that the THE system can be advantageous in improving enzyme efficiency, especially if heating the entire system is impossible or could lead to degradation of substrates or damage of components, such as in vitro bioanalysis of frangible molecules or in vivo diagnosis.

How to resurrect ancestral proteins as proxies for ancient biogeochemistry

Throughout the history of life, enzymes have served as the primary molecular mediators of biogeochemical cycles by catalyzing the metabolic pathways that interact with geochemical substrates. The byproducts of enzymatic activities have been preserved as chemical and isotopic signatures in the geologic record. However, interpretations of these signatures are limited by the assumption that such enzymes have remained functionally conserved over billions of years of molecular evolution. By reconstructing ancient genetic sequences in conjunction with laboratory enzyme resurrection, preserved biogeochemical signatures can instead be related to experimentally constrained, ancestral enzymatic properties. We may thereby investigate instances within molecular evolutionary trajectories potentially tied to significant biogeochemical transitions evidenced in the geologic record. Here, we survey recent enzyme resurrection studies to provide a reasoned assessment of areas of success and common pitfalls relevant to ancient biogeochemical applications. We conclude by considering the Great Oxidation Event, which provides a constructive example of a significant biogeochemical transition that warrants investigation with ancestral enzyme resurrection. This event also serves to highlight the pitfalls of facile interpretation of paleophenotype models and data, as applied to two examples of enzymes that likely both influenced and were influenced by the rise of atmospheric oxygen - RuBisCO and nitrogenase.

Theory of Nonequilibrium Free Energy Transduction by Molecular Machines

Biomolecular machines are protein complexes that convert between different forms of free energy. They are utilized in nature to accomplish many cellular tasks. As isothermal nonequilibrium stochastic objects at low Reynolds number, they face a distinct set of challenges compared with more familiar human-engineered macroscopic machines. Here we review central questions in their performance as free energy transducers, outline theoretical and modeling approaches to understand these questions, identify both physical limits on their operational characteristics and design principles for improving performance, and discuss emerging areas of research.

Emergence of a Negative Activation Heat Capacity during Evolution of a Designed Enzyme

Temperature influences the reaction kinetics and evolvability of all enzymes. To understand how evolution shapes the thermodynamic drivers of catalysis, we optimized the modest activity of a computationally designed enzyme for an elementary proton-transfer reaction by nearly 4 orders of magnitude over 9 rounds of mutagenesis and screening. As theorized for primordial enzymes, the catalytic effects of the original design were almost entirely enthalpic in origin, as were the rate enhancements achieved by laboratory evolution. However, the large reductions in ΔH^⧧ were partially offset by a decrease in TΔS^⧧ and unexpectedly accompanied by a negative activation heat capacity, signaling strong adaptation to the operating temperature. These findings echo reports of temperature-dependent activation parameters for highly evolved natural enzymes and are relevant to explanations of enzymatic catalysis and adaptation to changing thermal environments.

Establishment of mesophilic-like catalytic properties in a thermophilic enzyme without affecting its thermal stability

Thermophilic enzymes are generally more thermally stable but are less active at moderate temperatures than are their mesophilic counterparts. Thermophilic enzymes with improved low-temperature activity that retain their high stability would serve as useful tools for industrial processes especially when robust biocatalysts are required. Here we show an effective way to explore amino acid substitutions that enhance the low-temperature catalytic activity of a thermophilic enzyme, based on a pairwise sequence comparison of thermophilic/mesophilic enzymes. One or a combination of amino acid(s) in 3-isopropylmalate dehydrogenase from the extreme thermophile Thermus thermophilus was/were substituted by a residue(s) found in the Escherichia coli enzyme at the same position(s). The best mutant, which contained three amino acid substitutions, showed a 17-fold higher specific activity at 25 °C compared to the original wild-type enzyme while retaining high thermal stability. The kinetic and thermodynamic parameters of the mutant showed similar patterns along the reaction coordinate to those of the mesophilic enzyme. We also analyzed the residues at the substitution sites from a structural and phylogenetic point of view.

Following the Evolutionary Track of a Highly Specific l-Arginine Oxidase by Reconstruction and Biochemical Analysis of Ancestral and Native Enzymes

Following the evolutionary track of enzymes can help elucidate how enzymes attain their characteristic functions, such as thermal adaptation and substrate selectivity, during the evolutionary process. Ancestral sequence reconstruction (ASR) is effective for following evolutionary processes if sufficient sequence data are available. Selecting sequences from the data to generate a curated sequence library is necessary for the successful design of artificial proteins by ASR. In this study, we tried to follow the evolutionary track of l-arginine oxidase (AROD), a flavin adenine dinucleotide (FAD)-dependent amino acid oxidase (LAAO) that exhibits high specificity for l-arginine. The library was generated by selecting sequences in which the 15th, 50th, 332nd, and 580th residues are Gly, Ser, Trp, and Thr, respectively. We excluded sequences that are either extremely short or long and those with a low degree of sequence identity. Three ancestral ARODs (AncARODn0, AncARODn1, and AncARODn2) were designed using the library. Subsequently, we expressed the ancestral ARODs as well as native Oceanobacter kriegii AROD (OkAROD) in bacteria. AncARODn0 is phylogenetically most remote from OkAROD, whereas AncARODn2 is most similar to OkAROD. Thermal stability was gradually increased by extending AROD sequences back to the progenitor, while the temperature at which the residual activity is half of the maximum measured activity (T _1/2) of AncARODn0 was >20°C higher than that of OkAROD. Remarkably, only AncARODn0 exhibited broad substrate selectivity similar to that of conventional promiscuous LAAO. Taken together, our findings led us to infer that AROD may have evolved from a highly thermostable and promiscuous LAAO.IMPORTANCE In this study, we attempted to infer the molecular evolution of a recently isolated FAD-dependent l-arginine oxidase (AROD) that oxidizes l-arginine to 2-ketoarginine. Utilizing 10 candidate AROD sequences, we obtained a total of three ancestral ARODs. In addition, one native AROD was obtained by cloning one of the candidate ARODs. The candidate sequences were selected utilizing a curation method defined in this study. All the ARODs were successfully expressed in Escherichia coli for analysis of their biochemical functions. The catalytic activity of our bacterially expressed ancestral ARODs suggests that our ASR was successful. The ancestral AROD that is phylogenetically most remote from a native AROD has the highest thermal stability and substrate promiscuity. Our findings led us to infer that AROD evolved from a highly thermostable and promiscuous LAAO. As an application, we can design artificial ARODs with improved functions compared with those of native ones.

Substrate inhibition imposes fitness penalty at high protein stability

Proteins are only moderately stable. It has long been debated whether this narrow range of stabilities is solely a result of neutral drift toward lower stability or purifying selection against excess stability-for which no experimental evidence was found so far-is also at work. Here, we show that mutations outside the active site in the essential Escherichia coli enzyme adenylate kinase (Adk) result in a stability-dependent increase in substrate inhibition by AMP, thereby impairing overall enzyme activity at high stability. Such inhibition caused substantial fitness defects not only in the presence of excess substrate but also under physiological conditions. In the latter case, substrate inhibition caused differential accumulation of AMP in the stationary phase for the inhibition-prone mutants. Furthermore, we show that changes in flux through Adk could accurately describe the variation in fitness effects. Taken together, these data suggest that selection against substrate inhibition and hence excess stability may be an important factor determining stability observed for modern-day Adk.

Protein engineers turned evolutionists-the quest for the optimal starting point

The advent of laboratory directed evolution yielded a fruitful crosstalk between the disciplines of molecular evolution and bio-engineering. Here, we outline recent developments in both disciplines with respect to how one can identify the best starting points for directed evolution, such that highly efficient and robust tailor-made enzymes can be obtained with minimal optimization. Directed evolution studies have highlighted essential features of engineer-able enzymes: highly stable, mutationally robust enzymes with the capacity to accept a broad range of substrates. Robust, evolvable enzymes can be inferred from the natural sequence record. Broad substrate spectrum relates to conformational plasticity and can also be predicted by phylogenetic analyses and/or by computational design. Overall, an increasingly powerful toolkit is becoming available for identifying optimal starting points including network analyses of enzyme superfamilies and other bioinformatics methods.

Enhancing the Thermostability of Highly Active and Glucose-Tolerant β-Glucosidase Ks5A7 by Directed Evolution for Good Performance of Three Properties

A high-performance β-glucosidase for efficient cellulose hydrolysis needs to excel in thermostability, catalytic efficiency, and resistance to glucose inhibition. However, it is challenging to achieve superb properties in all three aspects in a single enzyme. In this study, a hyperactive and glucose-tolerant β-glucosidase Ks5A7 was employed as the starting point. Four rounds of random mutagenesis were then performed, giving rise to a thermostable mutant 4R1 with five amino acid substitutions. The half-life of 4R1 at 50 °C is 8640-fold that of Ks5A7 (144 h vs 1 min). Meanwhile, 4R1 had a higher specific activity (374.26 vs 243.18 units·mg^-1) than the wild type with a similar glucose tolerance. When supplemented to Celluclast 1.5L, the mutant significantly enhanced the hydrolysis of pretreated sugar cane bagasse, improving the released glucose concentration by 44%. With excellent performance in thermostability, activity, and glucose tolerance, 4R1 will serve as an exceptional catalyst for industrial applications.

Biophysical Spandrels form a Hot-Spot for Kosmotropic Mutations in Bacteriophage Thermal Adaptation

Temperature plays a dominating role in protein structure and function, and life has evolved myriad strategies to adapt proteins to environmental thermal stress. Cellular systems can utilize kosmotropic osmolytes, the products of complex biochemical pathways, to act as chemical chaperones. These extrinsic molecules, e.g., trehalose, alter local water structure to modulate the strength of the hydrophobic effect and increase protein stability. In contrast, simpler genetic systems must rely on intrinsic mutation to affect protein stability. In naturally occurring microvirid bacteriophages of the subfamily Bullavirinae, capsid stability is randomly distributed across the phylogeny, suggesting it is not phylogenetically linked and could be altered through adaptive mutation. We hypothesized that these phages could utilize an adaptive mechanism that mimics the stabilizing effects of the kosmotrope trehalose through mutation. Kinetic stability of wild-type ID8, a relative of ΦX174, displays a saturable response to trehalose. Thermal adaptation mutations in ID8 improve capsid stability and reduce responsiveness to trehalose suggesting the mutations move stability closer to the kosmotropic saturation point, mimicking the kosmotropic effect of trehalose. These mutations localize to and modulate the hydrophobicity of a cavern formation at the interface of phage coat and spike proteins-an evolutionary spandrel. Across a series of genetically distinct phages, responsiveness to trehalose correlates positively with cavern hydrophobicity suggesting that the level of hydrophobicity of the cavern may provide a biophysical gating mechanism constraining or permitting adaptation in a lineage-specific manner. Our results demonstrate that a single mutation can exploit pre-existing, non-adaptive structural features to mimic the adaptive effects of complex biochemical pathways.

Nonequivalence of Second Sphere "Noncatalytic" Residues in Pentaerythritol Tetranitrate Reductase in Relation to Local Dynamics Linked to H-Transfer in Reactions with NADH and NADPH Coenzymes

Many enzymes that catalyze hydride transfer reactions work via a mechanism dominated by quantum mechanical tunneling. The involvement of fast vibrational modes of the reactive complex is often inferred in these reactions, as in the case of the NAD(P)H-dependent pentaerythritol tetranitrate reductase (PETNR). Herein, we interrogated the H-transfer mechanism in PETNR by designing conservative (L25I and I107L) and side chain shortening (L25A and I107A) PETNR variants and using a combination of experimental approaches (stopped-flow rapid kinetics, X-ray crystallography, isotope/temperature dependence studies of H-transfer and NMR spectroscopy). X-ray data show subtle changes in the local environment of the targeted side chains but no major structural perturbation caused by mutagenesis of these two second sphere active site residues. However, temperature dependence studies of H-transfer revealed a coenzyme-specific and complex thermodynamic equilibrium between different reactive configurations in PETNR–coenzyme complexes. We find that mutagenesis of these second sphere “noncatalytic” residues affects differently the reactivity of PETNR with NADPH and NADH coenzymes. We attribute this to subtle, dynamic structural changes in the PETNR active site, the effects of which impact differently in the nonequivalent reactive geometries of PETNR−NADH and PETNR−NADPH complexes. This inference is confirmed through changes observed in the NMR chemical shift data for PETNR complexes with unreactive 1,4,5,6-tetrahydro-NAD(P) analogues. We show that H-transfer rates can (to some extent) be buffered through entropy–enthalpy compensation, but that use of integrated experimental tools reveals hidden complexities that implicate a role for dynamics in this relatively simple H-transfer reaction. Similar approaches are likely to be informative in other enzymes to understand the relative importance of (distal) hydrophobic side chains and dynamics in controlling the rates of enzymatic H-transfer.

Reconstructing the evolutionary history of F420-dependent dehydrogenases

During the last decade the number of characterized F₄₂₀-dependent enzymes has significantly increased. Many of these deazaflavoproteins share a TIM-barrel fold and are structurally related to FMN-dependent luciferases and monooxygenases. In this work, we traced the origin and evolutionary history of the F₄₂₀-dependent enzymes within the luciferase-like superfamily. By a thorough phylogenetic analysis we inferred that the F₄₂₀-dependent enzymes emerged from a FMN-dependent common ancestor. Furthermore, the data show that during evolution, the family of deazaflavoproteins split into two well-defined groups of enzymes: the F₄₂₀-dependent dehydrogenases and the F₄₂₀-dependent reductases. By such event, the dehydrogenases specialized in generating the reduced deazaflavin cofactor, while the reductases employ the reduced F₄₂₀ for catalysis. Particularly, we focused on investigating the dehydrogenase subfamily and demonstrated that this group diversified into three types of dehydrogenases: the already known F₄₂₀-dependent glucose-6-phosphate dehydrogenases, the F₄₂₀-dependent alcohol dehydrogenases, and the sugar-6-phosphate dehydrogenases that were identified in this study. By reconstructing and experimentally characterizing ancestral and extant representatives of F₄₂₀-dependent dehydrogenases, their biochemical properties were investigated and compared. We propose an evolutionary path for the emergence and diversification of the TIM-barrel fold F₄₂₀-dependent dehydrogenases subfamily.

Linear Eyring Plots Conceal a Change in the Rate-Limiting Step in an Enzyme Reaction

The temperature dependence of psychrophilic and mesophilic (R)-3-hydroxybutyrate dehydrogenase steady-state rates yields nonlinear and linear Eyring plots, respectively. Solvent viscosity effects and multiple- and single-turnover pre-steady-state kinetics demonstrate that while product release is rate-limiting at high temperatures for the psychrophilic enzyme, either interconversion between enzyme–substrate and enzyme–product complexes or a step prior to it limits the rate at low temperatures. Unexpectedly, a similar change in the rate-limiting step is observed with the mesophilic enzyme, where a step prior to chemistry becomes rate-limiting at low temperatures. This observation may have implications for past and future interpretations of temperature–rate profiles.

Ancestral diterpene cyclases show increased thermostability and substrate acceptance

Bacterial diterpene cyclases are receiving increasing attention in biocatalysis and synthetic biology for the sustainable generation of complex multicyclic building blocks. Herein, we explore the potential of ancestral sequence reconstruction (ASR) to generate remodeled cyclases with enhanced stability, activity, and promiscuity. Putative ancestors of spiroviolene synthase, a bacterial class I diterpene cyclase, display an increased yield of soluble protein of up to fourfold upon expression in the model organism Escherichia coli. Two of the resurrected enzymes, with an estimated age of approximately 1.7 million years, display an upward shift in thermostability of 7-13 °C. Ancestral spiroviolene synthases catalyze cyclization of the natural C₂₀ -substrate geranylgeranyl diphosphate (GGPP) and also accept C₁₅ farnesyl diphosphate (FPP), which is not converted by the extant enzyme. In contrast, the consensus sequence generated from the corresponding multiple sequence alignment was found to be inactive toward both substrates. Mutation of a nonconserved position within the aspartate-rich motif of the reconstructed ancestral cyclases was associated with modest effects on activity and relative substrate specificity (i.e., k_cat /K_M for GGPP over k_cat /K_M for FPP). Kinetic analyses performed at different temperatures reveal a loss of substrate saturation, when going from the ancestor with highest thermostability to the modern enzyme. The kinetics data also illustrate how an increase in temperature optimum of biocatalysis is reflected in altered entropy and enthalpy of activation. Our findings further highlight the potential and limitations of applying ASR to biosynthetic machineries in secondary metabolism.

Ancient mechanisms for the evolution of the bicoid homeodomain's function in fly development

The ancient mechanisms that caused developmental gene regulatory networks to diversify among distantly related taxa are not well understood. Here we use ancestral protein reconstruction, biochemical experiments, and developmental assays of transgenic animals carrying reconstructed ancestral genes to investigate how the transcription factor Bicoid (Bcd) evolved its central role in anterior-posterior patterning in flies. We show that most of Bcd's derived functions are attributable to evolutionary changes within its homeodomain (HD) during a phylogenetic interval >140 million years ago. A single substitution from this period (Q50K) accounts almost entirely for the evolution of Bcd's derived DNA specificity in vitro. In transgenic embryos expressing the reconstructed ancestral HD, however, Q50K confers activation of only a few of Bcd's transcriptional targets and yields a very partial rescue of anterior development. Adding a second historical substitution (M54R) confers regulation of additional Bcd targets and further rescues anterior development. These results indicate that two epistatically interacting mutations played a major role in the evolution of Bcd's controlling regulatory role in early development. They also show how ancestral sequence reconstruction can be combined with in vivo characterization of transgenic animals to illuminate the historical mechanisms of developmental evolution.

Structure and dynamics conspire in the evolution of affinity between intrinsically disordered proteins

In every established species, protein-protein interactions have evolved such that they are fit for purpose. However, the molecular details of the evolution of new protein-protein interactions are poorly understood. We have used nuclear magnetic resonance spectroscopy to investigate the changes in structure and dynamics during the evolution of a protein-protein interaction involving the intrinsically disordered CREBBP (CREB-binding protein) interaction domain (CID) and nuclear coactivator binding domain (NCBD) from the transcriptional coregulators NCOA (nuclear receptor coactivator) and CREBBP/p300, respectively. The most ancient low-affinity "Cambrian-like" [540 to 600 million years (Ma) ago] CID/NCBD complex contained less secondary structure and was more dynamic than the complexes from an evolutionarily younger "Ordovician-Silurian" fish ancestor (ca. 440 Ma ago) and extant human. The most ancient Cambrian-like CID/NCBD complex lacked one helix and several interdomain interactions, resulting in a larger solvent-accessible surface area. Furthermore, the most ancient complex had a high degree of millisecond-to-microsecond dynamics distributed along the entire sequences of both CID and NCBD. These motions were reduced in the Ordovician-Silurian CID/NCBD complex and further redistributed in the extant human CID/NCBD complex. Isothermal calorimetry experiments show that complex formation is enthalpically favorable and that affinity is modulated by a largely unfavorable entropic contribution to binding. Our data demonstrate how changes in structure and motion conspire to shape affinity during the evolution of a protein-protein complex and provide direct evidence for the role of structural, dynamic, and frustrational plasticity in the evolution of interactions between intrinsically disordered proteins.

Coevolution of both Thermostability and Activity of Polyphosphate Glucokinase from Thermobifida fusca YX

Thermostability and specific activity of enzymes are two of the most important properties for industrial biocatalysts. Here, we developed a petri dish-based double-layer high-throughput screening (HTS) strategy for rapid identification of desired mutants of polyphosphate glucokinase (PPGK) from a thermophilic actinobacterium, Thermobifida fusca YX, with both enhanced thermostability and activity. Escherichia coli colonies representing a PPGK mutant library were grown on the first-layer Phytagel-based plates, which can remain solid for 1 h, even at heat treatment temperatures of more than 100°C. The second layer that was poured on the first layer contained agarose, substrates, glucose 6-phosphate dehydrogenase (G6PDH), the redox dye tetranitroblue tetrazolium (TNBT), and phenazine methosulfate. G6PDH was able to oxidize the product from the PPGK-catalyzed reaction and generate NADH, which can be easily examined by a TNBT-based colorimetric assay. The best mutant obtained after four rounds of directed evolution had a 7,200-fold longer half-life at 55°C, 19.8°C higher midpoint of unfolding temperature (T_m ), and a nearly 3-fold enhancement in specific activities compared to those of the wild-type PPGK. The best mutant was used to produce 9.98 g/liter myo-inositol from 10 g/liter glucose, with a theoretical yield of 99.8%, along with two other hyperthermophilic enzymes at 70°C. This PPGK mutant featuring both great thermostability and high activity would be useful for ATP-free production of glucose 6-phosphate or its derived products.IMPORTANCE Polyphosphate glucokinase (PPGK) is an enzyme that transfers a terminal phosphate group from polyphosphate to glucose, producing glucose 6-phosphate. A petri dish-based double-layer high-throughput screening strategy was developed by using ultrathermostable Phytagel as the first layer instead of agar or agarose, followed by a redox dye-based assay for rapid identification of ultrathermostable PPGK mutants. The best mutant featuring both great thermostability and high activity could produce glucose 6-phosphate from glucose and polyphosphate without in vitro ATP regeneration.