The structure of Pyrococcus furiosus glutamate dehydrogenase reveals a key role for ion-pair networks in maintaining enzyme stability at extreme temperatures

Molecular dynamics simulations at high temperatures of the Aeropyrum pernix L7Ae thermostable protein: Insight into the unfolding pathway

Most life forms on earth live at temperatures below 50 °C. Within these organisms are proteins that form the three-dimensional structures essential to their biological activity and function. However, some thermophilic life forms can resist higher temperatures and have corresponding adaptations to preserve protein function at these high temperatures. Among the structural factors responsible for this resistance of thermophilic proteins to high temperatures is the presence of additional hydrogen bonds in the thermophilic proteins, which means that the structure of the protein is more resistant to unfolding. Similarly, thermostable proteins are rich in structure-stabilizing salt bridges and/or disulfide bridges. In this context, we perform multiple replica molecular dynamics simulations at different temperatures on the Aeropyrum pernix (L7Ae) protein (from the crenarchaeal species A. pernix), known for its high melting temperature, and this in the aim to elucidate the structural factors responsible for its high thermostability. The results reveal that between the most sensitive regions of the protein to the increase of temperature are the loops L1, and L5, which surround the hydrophobic core region of the protein, besides the loop L9, and the C-terminal α5 region. This latter is the longer alpha helix of the protein secondary structure motifs and it is the first to be denaturated at 450 K, while the rest of the protein secondary structure motifs at this temperature were intact. The mechanism of unfolding that follows this protein at 550 K is similar to other thermophile proteins found in literature, with the opening of the loops that surround the hydrophobic core of the protein. So, the latter is completely exposed to the solvent, and partially denatured. The total denaturation process of the protein takes an average time of 40 ns to be achieved. Our investigation also shows that all the calculated salt bridges, with distances less than or equal to 6 A°, are on the periphery part of the protein, exposed to the solvent. However, the hydrophobic core of the protein is not involved in the formation of salt bridges, but rather with formation of some important hydrogen bondings that still persist even at 450 K. So, optimizing hydrogen bonding, near or within the core region, at high temperatures is a strategy that follows this thermostable protein to protect its hydrophobic core from denaturation, and ensure the thermal stability of the protein.

Structural Evolution of Primate Glutamate Dehydrogenase 2 as Revealed by In Silico Predictions and Experimentally Determined Structures

Glutamate dehydrogenase (GDH) interconverts glutamate to a-ketoglutarate and ammonia, interconnecting amino acid and carbohydrate metabolism. In humans, two functional GDH genes, GLUD1 and GLUD2, encode for hGDH1 and hGDH2, respectively. GLUD2 evolved from retrotransposition of the GLUD1 gene in the common ancestor of modern apes. These two isoenzymes are involved in the pathophysiology of human metabolic, neoplastic, and neurodegenerative disorders. The 3D structures of hGDH1 and hGDH2 have been experimentally determined; however, no information is available about the path of GDH2 structure changes during primate evolution. Here, we compare the structures predicted by the AlphaFold Colab method for the GDH2 enzyme of modern apes and their extinct primate ancestors. Also, we analyze the individual effect of amino acid substitutions emerging during primate evolution. Our most important finding is that the predicted structure of GDH2 in the common ancestor of apes was the steppingstone for the structural evolution of primate GDH2s. Two changes with a strong functional impact occurring at the first evolutionary step, Arg443Ser and Gly456Ala, had a destabilizing and stabilizing effect, respectively, making this step the most important one. Subsequently, GDH2 underwent additional modifications that fine-tuned its enzymatic properties to adapt to the functional needs of modern-day primate tissues.

A newly introduced salt bridge cluster improves structural and biophysical properties of de novo TIM barrels

Protein stability can be fine‐tuned by modifying different structural features such as hydrogen‐bond networks, salt bridges, hydrophobic cores, or disulfide bridges. Among these, stabilization by salt bridges is a major challenge in protein design and engineering since their stabilizing effects show a high dependence on the structural environment in the protein, and therefore are difficult to predict and model. In this work, we explore the effects on structure and stability of an introduced salt bridge cluster in the context of three different de novo TIM barrels. The salt bridge variants exhibit similar thermostability in comparison with their parental designs but important differences in the conformational stability at 25°C can be observed such as a highly stabilizing effect for two of the proteins but a destabilizing effect to the third. Analysis of the formed geometries of the salt bridge cluster in the crystal structures show either highly ordered salt bridge clusters or only single salt bridges. Rosetta modeling of the salt bridge clusters results in a good prediction of the tendency on stability changes but not the geometries observed in the three‐dimensional structures. The results show that despite the similarities in protein fold, the salt bridge clusters differently influence the structural and stability properties of the de novo TIM barrel variants depending on the structural background where they are introduced.

PDB Code(s): 7OSU, 7OT7, 7OSV, 7OT8 and 7P12;

Enzymatic analysis of truncation mutants of a type II pullulanase from Bifidobacterium adolescentis P2P3, a resistant starch-degrading gut bacterium

A putative type II pullulanase gene, pulP, was identified in Bifidobacterium adolescentis P2P3. PulP possesses an α-amylase domain at the N-terminus and a pullulanase type I domain at the C-terminus, as well as three carbohydrate-binding modules (one CBM25 and two CBM41s) between them. The native PulP and four truncated mutant recombinant proteins (PulPΔCΔP, PulPΔP, PulPΔAΔC, and PulPΔA), in which each of the two catalytic domains and/or the CBMs were deleted, were produced in Escherichia coli and their specific properties were characterized. The removal of either catalytic domain abolished the corresponding catalytic activity of the wild-type enzyme. Deletion of the C-terminal domain resulted in a drastic decrease in the optimal temperature and thermostability, indicating that the pullulanase domain might be related to the temperature dependency of the enzyme. In addition, the elimination of the CBMs in the mutant proteins led to a loss of binding affinity toward raw substrates as well as the loss of their hydrolysis activities compared to the wild-type enzyme. HPAEC and TLC analyses proved that PulP and its mutants could hydrolyze α-glucans into maltotriose as their main product. These results suggest that PulP may play an important role in α-glucan metabolism in B. adolescentis P2P3.

Thermostable lipases and their dynamics of improved enzymatic properties

Thermal stability is one of the most desirable characteristics in the search for novel lipases. The search for thermophilic microorganisms for synthesising functional enzyme biocatalysts with the ability to withstand high temperature, and capacity to maintain their native state in extreme conditions opens up new opportunities for their biotechnological applications. Thermophilic organisms are one of the most favoured organisms, whose distinctive characteristics are extremely related to their cellular constituent particularly biologically active proteins. Modifications on the enzyme structure are critical in optimizing the stability of enzyme to thermophilic conditions. Thermostable lipases are one of the most favourable enzymes used in food industries, pharmaceutical field, and actively been studied as potential biocatalyst in biodiesel production and other biotechnology application. Particularly, there is a trade-off between the use of enzymes in high concentration of organic solvents and product generation. Enhancement of the enzyme stability needs to be achieved for them to maintain their enzymatic activity regardless the environment. Various approaches on protein modification applied since decades ago conveyed a better understanding on how to improve the enzymatic properties in thermophilic bacteria. In fact, preliminary approach using advanced computational analysis is practically conducted before any modification is being performed experimentally. Apart from that, isolation of novel extremozymes from various microorganisms are offering great frontier in explaining the crucial native interaction within the molecules which could help in protein engineering. In this review, the thermostability prospect of lipases and the utility of protein engineering insights into achieving functional industrial usefulness at their high temperature habitat are highlighted. Similarly, the underlying thermodynamic and structural basis that defines the forces that stabilize these thermostable lipase is discussed. KEY POINTS: • The dynamics of lipases contributes to their non-covalent interactions and structural stability. • Thermostability can be enhanced by well-established genetic tools for improved kinetic efficiency. • Molecular dynamics greatly provides structure-function insights on thermodynamics of lipase.

Structure elucidation and docking analysis of 5M mutant of T1 lipase Geobacillus zalihae

5M mutant lipase was derived through cumulative mutagenesis of amino acid residues (D43E/T118N/E226D/E250L/N304E) of T1 lipase from Geobacillus zalihae. A previous study revealed that cumulative mutations in 5M mutant lipase resulted in decreased thermostability compared to wild-type T1 lipase. Multiple amino acids substitution might cause structural destabilization due to negative cooperation. Hence, the three-dimensional structure of 5M mutant lipase was elucidated to determine the evolution in structural elements caused by amino acids substitution. A suitable crystal for X-ray diffraction was obtained from an optimized formulation containing 0.5 M sodium cacodylate trihydrate, 0.4 M sodium citrate tribasic pH 6.4 and 0.2 M sodium chloride with 2.5 mg/mL protein concentration. The three-dimensional structure of 5M mutant lipase was solved at 2.64 Å with two molecules per asymmetric unit. The detailed analysis of the structure revealed that there was a decrease in the number of molecular interactions, including hydrogen bonds and ion interactions, which are important in maintaining the stability of lipase. This study facilitates understanding of and highlights the importance of hydrogen bonds and ion interactions towards protein stability. Substrate specificity and docking analysis on the open structure of 5M mutant lipase revealed changes in substrate preference. The molecular dynamics simulation of 5M-substrates complexes validated the substrate preference of 5M lipase towards long-chain p-nitrophenyl-esters.

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.

A second shell residue modulates a conserved ATP-binding site with radically different affinities for ATP

BACKGROUND

Prediction of ligand binding and design of new function in enzymes is a time-consuming and expensive process. Crystallography gives the impression that proteins adopt a fixed shape, yet enzymes are functionally dynamic. Molecular dynamics offers the possibility of probing protein movement while predicting ligand binding. Accordingly, we choose the bacterial F₁F_o ATP synthase ε subunit to unravel why ATP affinity by ε subunits from Bacillus subtilis and Bacillus PS3 differs ~500-fold, despite sharing identical sequences at the ATP-binding site.

METHODS

We first used the Bacillus PS3 ε subunit structure to model the B. subtilis ε subunit structure and used this to explore the utility of molecular dynamics (MD) simulations to predict the influence of residues outside the ATP binding site. To verify the MD predictions, point mutants were made and ATP binding studies were employed.

RESULTS

MD simulations predicted that E102 in the B. subtilis ε subunit, outside of the ATP binding site, influences ATP binding affinity. Engineering E102 to alanine or arginine revealed a ~10 or ~54 fold increase in ATP binding, respectively, confirming the MD prediction that E102 drastically influences ATP binding affinity.

CONCLUSIONS

These findings reveal how MD can predict how changes in the "second shell" residues around substrate binding sites influence affinity in simple protein structures. Our results reveal why seemingly identical ε subunits in different ATP synthases have radically different ATP binding affinities.

GENERAL SIGNIFICANCE

This study may lead to greater utility of molecular dynamics as a tool for protein design and exploration of protein design and function.

Biochemical characterization of a thermophilic hyaluronate lyase TcHly8C from Thermasporomyces composti DSM22891

Hyaluronic acid (HA) is an anionic linear polysaccharide abundantly distributed in the extracellular matrix of mammalian connective, growing, and tumor tissues. Hyaluronidase is used as an important drug diffusion promoter and a tool enzyme to produce HA oligosaccharides. However, there is no thermostable hyaluronidase suitable for application to date. In this study, a thermophilic hyaluronate lyase, TcHly8C, from Thermasporomyces composti DSM22891 was expressed in Escherichia coli. The recombinant TcHly8C was most active at 70 °C, and it retained about 30% of initial activity after incubation at 60 °C for 28 days. The half-lives of TcHly8C at 60 °C and 70 °C were 16.1 d and 2.3 h, respectively. The optimum pH of TcHly8C is 5.93, and it was stable at pH 6.15-10.90. The presence of Mg²⁺ could enhance its enzymatic activity significantly. K_m, k_cat, and k_cat/K_m of TcHly8C towards HA were 3.69 mg∙ml^-1, 17.82 s^-1, and 4.82 ml∙mg^-1∙s^-1, respectively. TcHly8C degraded HA in an exolytic mode, and the end product was unsaturated HA disaccharide (ΔUA-GlcNAc). Overall, our results show that TcHly8C is the first reported PL8 exo-type hyaluronate lyase with high thermostability, which provides a potential enzyme used in medicine and production of HA oligosaccharides.

Finding the generalized molecular principles of protein thermal stability

Are there any generalized molecular principles of thermal adaptation? Here, integrating the concepts of structural bioinformatics, sequence analysis, and classical knot theory, we develop a robust computational framework that seeks for mechanisms of thermal adaptation by comparing orthologous mesophilic-thermophilic and mesophilic-hyperthermophilic proteins of remarkable structural and topological similarities, and still leads us to context-independent results. A comprehensive analysis of 4741 high-resolution, non-redundant X-ray crystallographic structures collected from 11 hyperthermophilic, 32 thermophilic and 53 mesophilic prokaryotes unravels at least five "nearly universal" signatures of thermal adaptation, irrespective of the enormous sequence, structure, and functional diversity of the proteins compared. A careful investigation further extracts a set of amino acid changes that can potentially enhance protein thermal stability, and remarkably, these mutations are overrepresented in protein crystallization experiments, in disorder-to-order transitions and in engineered thermostable variants of existing mesophilic proteins. These results could be helpful to find a precise, global picture of thermal adaptation.

Thermostable iron oxide nanoparticle synthesis within recombinant ferritins from the hyperthermophile Pyrococcus yayanosii CH1

Thermostable nanoparticles have numerous applications in catalysis and in the oil/gas industry. However, synthesizing these nanoparticles requires expensive polymers. Here, a novel thermostable ferritin named PcFn, originally from the hyperthermophilic archaeon Pyrococcus yayanosii CH1, was overexpressed in Escherichia coli, purified and characterized, which could successfully direct the synthesis of thermostable magnetoferritins (M-PcFn) with monodispersed iron oxide nanoparticles in one step. Transmission electron microscopy and magnetic measurements show that the cores of the M-PcFn have an average diameter of 4.7 nm, are well-crystalline and superparamagnetic. Both the PcFn and M-PcFn can resist temperatures up to 110 °C, which is significantly higher than for human H-chain ferritin (HFn) and M-HFn, and comparable to temperatures previously reported for Pyrococcus furiosus ferritin (PfFn) and M-PfFn. After heating at 110 °C for 30 minutes, PcFn and M-PcFn maintained their secondary structures and PcFn retained 87.4% of its iron uptake activity. This remarkable thermostability of PcFn and M-PcFn suggests potential applications in elevated temperature environments.

Dissecting the Antenna in Human Glutamate Dehydrogenase: Understanding Its Role in Subunit Communication and Allosteric Regulation

Glutamate dehydrogenase (GDH) is a homohexameric enzyme that catalyzes the reversible oxidative deamination of l-glutamate. While GDH is found in all living organisms, only that from animals is highly allosterically regulated by a wide array of metabolites. Because only animal GDH has a 50-residue antenna domain, we hypothesized that it was critical for allostery. To this end, we previously replaced the antenna with the loop found in bacteria, and the resulting chimera was no longer regulated by purine nucleotides. Hence, it seemed logical that the purpose of the antenna is to exert the subunit communication necessary for heterotrophic allosteric regulation. Here, we revisit the antenna deletion studies by retaining 10 more of the human GDH (hGDH) residues without adding the bacterial loop. Unexpectedly, the results were profoundly different than before. The basal activity of the mutant is only ∼13% of that of the wild type but ∼100 times more sensitive to all allosteric activators. In contrast, the mutant is still affected by all of the tested inhibitors to approximately the same degree. The resulting antenna-less mutant retained its negative cooperativity with respect to the coenzyme, again suggesting that intersubunit communication is intact. Finally, the mutant still exhibits substrate inhibition, albeit there are differences in the details. We present a model in which the majority of the antenna is not directly involved in allosteric regulation per se but rather may be responsible for improving enzymatic efficiency by acting as a conduit for substrate binding energy between subunits.

Engineering lipases for temperature adaptation: Structure function correlation

Bacillus lipases are industrially attractive enzymes due to their broad substrate specificity and optimum alkaline pH. However, narrow temperature range of action and low thermostability restrain their optimal use and thus, necessitate attention. Several laboratories are engaged in protein engineering of Bacillus lipases to generate variants with improved attributes for decades using techniques such as directed evolution or rational design. This review summarizes the effect of mutations on the conformational changes through in silico modeling and their manifestation with respect to various biochemical parameters. Various studies have been put together to develop a perspective on the molecular basis of biocatalysis of lipases holding industrial importance.

Aggregation-promoting conditions necessary to create the complexes by acylphosphatase from the hyperthermophile Sulfolobus solfataricus

The structural transition from the globular to the amyloid form of proteins requires aggregation-promoting conditions. The protein example of this category is acylphosphatase from the hyperthermophile Sulfolobus solfataricus. This protein represents a structure with a well-defined hydrophobic core. This is why the complexation (including oligomerization) of this protein is of low probability. The chain fragment participating in aggregation in comparison to the status with respect to the fuzzy oil drop model is discussed in this paper.

Artificial Thermostable D-Amino Acid Dehydrogenase: Creation and Application

Many kinds of NAD(P)⁺-dependent L-amino acid dehydrogenases have been so far found and effectively used for synthesis of L-amino acids and their analogs, and for their sensing. By contrast, similar biotechnological use of D-amino acid dehydrogenase (D-AADH) has not been achieved because useful D-AADH has not been found from natural resources. Recently, using protein engineering methods, an NADP⁺-dependent D-AADH was created from meso-diaminopimelate dehydrogenase (meso-DAPDH). The artificially created D-AADH catalyzed the reversible NADP⁺-dependent oxidative deamination of D-amino acids to 2-oxo acids. The enzyme, especially thermostable one from thermophiles, was efficiently applicable to synthesis of D-branched-chain amino acids (D-BCAAs), with high yields and optical purity, and was useful for the practical synthesis of ¹³C- and/or ¹⁵N-labeled D-BCAAs. The enzyme also made it possible to assay D-isoleucine selectively in a mixture of isoleucine isomers. Analyses of the three-dimensional structures of meso-DAPDH and D-AADH, and designed mutations based on the information obtained made it possible to markedly enhance enzyme activity and to create D-AADH homologs with desired reactivity profiles. The methods described here may be an effective approach to artificial creation of biotechnologically useful enzymes.

Structure and function of the type III pullulan hydrolase from Thermococcus kodakarensis

Pullulan-hydrolysing enzymes, more commonly known as debranching enzymes for starch and other polysaccharides, are of great interest and have been widely used in the starch-saccharification industry. Type III pullulan hydrolase from Thermococcus kodakarensis (TK-PUL) possesses both pullulanase and α-amylase activities. Until now, only two enzymes in this class, which are capable of hydrolysing both α-1,4- and α-1,6-glycosidic bonds in pullulan to produce a mixture of maltose, panose and maltotriose, have been described. TK-PUL shows highest activity in the temperature range 95-100°C and has a pH optimum in the range 3.5-4.2. Its unique ability to hydrolyse maltotriose into maltose and glucose has not been reported for other homologous enzymes. The crystal structure of TK-PUL has been determined at a resolution of 2.8 Å and represents the first analysis of a type III pullulan hydrolyse. The structure reveals that the last part of the N-terminal domain and the C-terminal domain are significantly different from homologous structures. In addition, the loop regions at the active-site end of the central catalytic domain are quite different. The enzyme has a well defined calcium-binding site and possesses a rare vicinal disulfide bridge. The thermostability of TK-PUL and its homologues may be attributable to several factors, including the increased content of salt bridges, helical segments, Pro, Arg and Tyr residues and the decreased content of serine.

Marginal stability drives irreversible unfolding of large multi-domain family 3 glycosylhydrolases from thermo-tolerant yeast

Protein folding is an extremely complex and fast, yet perfectly defined process, involving interplay of many intra and inter-molecular forces. In vitro, these molecular interactions are reversible for many proteins e.g., smaller and monomeric, organized into single domains. However, refolding of larger multi-domain/multimeric proteins is much more complicated, proceeds in a hierarchal way and is often irreversible. In a comparative study on two large, multi-domain and multimeric isozymes, β-glucosidase I (BGLI) and β-glucosidase II (BGLII) from Pichia etchellsii, we studied spontaneous and assisted refolding under three denaturing conditions viz. GdnHCl, alkaline pH and heat. During refolding, higher refolding yields were obtained for BGLII in case of pH induced unfolding (13.89%±0.25) than BGLI (6%±0.85) while for GdnHCl induced unfolding, refolding was marginal (BGLI=5%±0.5; BGLII=6%±0.69). Thermal unfolding was irreversible while assisted refolding also showed little structural gain for both proteins. When the apparent free energies of unfolding (ΔG_U^app) were calculated from GdnHCl unfolding data, their values were strikingly found to be lower (BGLI ΔG_U^app=3.02kcal/mol; BGLII ΔG_U^app=2.99kcal/mol) than reported for globular (ΔG_U=5-15kcal/mol)/multimeric proteins (ΔG_U=23-29kcal/mol) indicating marginal stability results in low refolding.

Structure and function of the thermostableL-asparaginase fromThermococcus kodakarensis

L-Asparaginases catalyse the hydrolysis of asparagine to aspartic acid and ammonia. In addition, L-asparaginase is involved in the biosynthesis of amino acids such as lysine, methionine and threonine. These enzymes have been used as chemotherapeutic agents for the treatment of acute lymphoblastic leukaemia and other haematopoietic malignancies since the tumour cells cannot synthesize sufficient L-asparagine and are thus killed by deprivation of this amino acid. L-Asparaginases are also used in the food industry and have potential in the development of biosensors, for example for asparagine levels in leukaemia. The thermostable type I L-asparaginase fromThermococcus kodakarensis(TkA) is composed of 328 amino acids and forms homodimers in solution, with the highest catalytic activity being observed at pH 9.5 and 85°C. It has aKmvalue of 5.5 mMfor L-asparagine, with no glutaminase activity being observed. The crystal structure of TkA has been determined at 2.18 Å resolution, confirming the presence of two α/β domains connected by a short linker region. The N-terminal domain contains a highly flexible β-hairpin which adopts `open' and `closed' conformations in different subunits of the solved TkA structure. In previously solved L-asparaginase structures this β-hairpin was only visible when in the `closed' conformation, whilst it is characterized with good electron density in all of the subunits of the TkA structure. A phosphate anion resides at the active site, which is formed by residues from both of the neighbouring monomers in the dimer. The high thermostability of TkA is attributed to the high arginine and salt-bridge content when compared with related mesophilic enzymes.

Density discriminates between thermophilic and mesophilic proteins

Despite an intense interest and a remarkable number of studies on the subject, the relationships between thermostability and (primary, secondary and tertiary) structure of proteins are still not fully understood. Here, comparing the protein density - defined by the ratio between the residue number and protein excluded volume - for a set of thermophilic/mesophilic pairs, we provide evidence that this property is connected to the optimal growth temperature. In particular, our results indicate that thermophilic proteins have - in general - a lower density with respect to the mesophilic counterparts, being such a correlation more pronounced for optimal growth temperature differences greater than 40°C. The effect of the protein thermostability changes on the molecular shape is also presented.

Complex network analysis of thermostable mutants of Bacillus subtilis Lipase A

Three-dimensional structures of proteins that regulate their functions can be modelled using complex network based approaches for understanding the structure-function relationship. The six mutants of the protein Lipase A from Bacillus subtilis, harbouring 2 to 12 mutations, retain their function at higher temperatures with negligible variation in their overall three-dimensional crystallographic structures. This enhanced thermostability of the mutants questions the structure-function paradigm. In this paper, a coarse-grained complex network approach is used to elucidate the structural basis of enhanced thermostability in the mutant proteins, by uncovering small but significant local changes distributed throughout the structure, rendering stability to the mutants at higher temperatures. Community structure analysis of the six mutant protein networks uncovers the specific reorganisations among the nodes/residues that occur, in absence of overall structural variations, which induce enhanced rigidity underlying the increased thermostability. This study offers a novel and significant application of complex network analysis that proposes to be useful in the understanding and designing of thermostable proteins.

Crystal structure of the 2-iminoglutarate-bound complex of glutamate dehydrogenase from Corynebacterium glutamicum

The NADP⁺ -dependent glutamate dehydrogenase from Corynebacterium glutamicum (CgGDH) is considered to be one of the key enzymes in the industrial fermentation of glutamate due to its high glutamate-producing activity. We determined the crystal structure of CgGDH complexed with NADP⁺ and 2-iminoglutarate. Among six subunits of hexameric CgGDH-binding NADP⁺ , only four subunits bind 2-iminoglutarate in a closed form, while the other two are in an open form. In the closed form, 2-iminoglutarate is bound to the substrate-binding site with the 2-imino group stacked by the nicotinamide ring of the coenzyme, suggesting a prehydride transfer state in a hypothesized reaction scheme with the imino intermediate. We also conducted MD simulations and provide insights into the extreme preference for the glutamate-producing reaction of CgGDH.

DATABASE

The atomic coordinate and structure factors have been deposited in the RCSB PDB database under the accession number 5GUD.

Physical and molecular bases of protein thermal stability and cold adaptation

The molecular bases of thermal and cold stability and adaptation, which allow proteins to remain folded and functional in the temperature ranges in which their host organisms live and grow, are still only partially elucidated. Indeed, both experimental and computational studies fail to yield a fully precise and global physical picture, essentially because all effects are context-dependent and thus quite intricate to unravel. We present a snapshot of the current state of knowledge of this highly complex and challenging issue, whose resolution would enable large-scale rational protein design.

Crystal structure of a chimaeric bacterial glutamate dehydrogenase

Glutamate dehydrogenases (EC 1.4.1.2-4) catalyse the oxidative deamination of L-glutamate to α-ketoglutarate using NAD(P)(+) as a cofactor. The bacterial enzymes are hexameric, arranged with 32 symmetry, and each polypeptide consists of an N-terminal substrate-binding segment (domain I) followed by a C-terminal cofactor-binding segment (domain II). The catalytic reaction takes place in the cleft formed at the junction of the two domains. Distinct signature sequences in the nucleotide-binding domain have been linked to the binding of NAD(+) versus NADP(+), but they are not unambiguous predictors of cofactor preference. In the absence of substrate, the two domains move apart as rigid bodies, as shown by the apo structure of glutamate dehydrogenase from Clostridium symbiosum. Here, the crystal structure of a chimaeric clostridial/Escherichia coli enzyme has been determined in the apo state. The enzyme is fully functional and reveals possible determinants of interdomain flexibility at a hinge region following the pivot helix. The enzyme retains the preference for NADP(+) cofactor from the parent E. coli domain II, although there are subtle differences in catalytic activity.

Structural and functional insights into the stationary-phase survival protein SurE, an important virulence factor of Brucella abortus

The stationary-phase survival protein SurE from Brucella abortus (BaSurE) is a metal-dependent phosphatase that is essential for the survival of this bacterium in the stationary phase of its life cycle. Here, BaSurE has been biochemically characterized and its crystal structure has been determined to a resolution of 1.9 Å. BaSurE was found to be a robust enzyme, showing activity over wide ranges of temperature and pH and with various phosphoester substrates. The active biomolecule is a tetramer and each monomer was found to consist of two domains: an N-terminal domain, which forms an approximate α + β fold, and a C-terminal domain that belongs to the α/β class. The active site lies at the junction of these two domains and was identified by the presence of conserved negatively charged residues and a bound Mg(2+) ion. Comparisons of BaSurE with its homologues have revealed both common features and differences in this class of enzymes. The number and arrangement of some of the equivalent secondary structures, which are seen to differ between BaSurE and its homologues, are responsible for a difference in the size of the active-site area and the overall oligomeric state of this enzyme in other organisms. As it is absent in mammals, it has the potential to be a drug target.

Analysis of Comparative Sequence and Genomic Data to Verify Phylogenetic Relationship and Explore a New Subfamily of Bacterial Lipases

Thermostable and organic solvent-tolerant enzymes have significant potential in a wide range of synthetic reactions in industry due to their inherent stability at high temperatures and their ability to endure harsh organic solvents. In this study, a novel gene encoding a true lipase was isolated by construction of a genomic DNA library of thermophilic Aneurinibacillus thermoaerophilus strain HZ into Escherichia coli plasmid vector. Sequence analysis revealed that HZ lipase had 62% identity to putative lipase from Bacillus pseudomycoides. The closely characterized lipases to the HZ lipase gene are from thermostable Bacillus and Geobacillus lipases belonging to the subfamily I.5 with ≤ 57% identity. The amino acid sequence analysis of HZ lipase determined a conserved pentapeptide containing the active serine, GHSMG and a Ca²⁺-binding motif, GCYGSD in the enzyme. Protein structure modeling showed that HZ lipase consisted of an α/β hydrolase fold and a lid domain. Protein sequence alignment, conserved regions analysis, clustal distance matrix and amino acid composition illustrated differences between HZ lipase and other thermostable lipases. Phylogenetic analysis revealed that this lipase represented a new subfamily of family I of bacterial true lipases, classified as family I.9. The HZ lipase was expressed under promoter P_lac using IPTG and was characterized. The recombinant enzyme showed optimal activity at 65°C and retained ≥ 97% activity after incubation at 50°C for 1h. The HZ lipase was stable in various polar and non-polar organic solvents.

Rigidity versus flexibility: the dilemma of understanding protein thermal stability

The role of fluctuations in protein thermostability has recently received considerable attention. In the current literature a dualistic picture can be found: thermostability seems to be associated with enhanced rigidity of the protein scaffold in parallel with the reduction of flexible parts of the structure. In contradiction to such arguments it has been shown by experimental studies and computer simulation that thermal tolerance of a protein is not necessarily correlated with the suppression of internal fluctuations and mobility. Both concepts, rigidity and flexibility, are derived from mechanical engineering and represent temporally insensitive features describing static properties, neglecting that relative motion at certain time scales is possible in structurally stable regions of a protein. This suggests that a strict separation of rigid and flexible parts of a protein molecule does not describe the reality correctly. In this work the concepts of mobility/flexibility versus rigidity will be critically reconsidered by taking into account molecular dynamics calculations of heat capacity and conformational entropy, salt bridge networks, electrostatic interactions in folded and unfolded states, and the emerging picture of protein thermostability in view of recently developed network theories. Last, but not least, the influence of high temperature on the active site and activity of enzymes will be considered.

Insights into the unfolding pathway and identification of thermally sensitive regions of phytase from Aspergillus niger by molecular dynamics simulations

Thermal stability is of great importance in the application of commercial phytases. Phytase A (PhyA) is a monomeric protein comprising twelve α-helices and ten β-sheets. Comparative molecular dynamics (MD) simulations (at 310, 350, 400, and 500 K) revealed that the thermal stability of PhyA from Aspergillus niger (A. niger) is associated with its conformational rigidity. The most thermally sensitive regions were identified as loops 8 (residues 83-106), 10 (161-174), 14 (224-230), 17 (306-331), and 24 (442-444), which are present on the surface of the protein. It was observed that solvent-exposed loops denature before or show higher flexibility than buried residues. We observed that PhyA begins to unfold at loops 8 and 14, which further extends to loop 24 at the C-terminus. The intense movement of loop 8 causes the helix H2 and beta-sheet B3 to fluctuate at high temperature. The high flexibility of the H2, H10, and H12 helices at high temperature resulted in complete denaturation. The high mobility of loop 14 easily transfers to the adjacent helices H7, H8, and H9, which fluctuate and partially unfold at high temperature (500 K). It was also observed that the salt bridges Asp110-Lys149, Asp205-Lys277, Asp335-Arg136, Asp416-Arg420, and Glu387-Arg400 are important influences on the structural stability but not the thermostability, as the lengths of these salt bridges did not increase with rising temperature. The salt bridges Glu125-Arg163, Asp299-Arg136, Asp266-Arg219, Asp339-Lys278, Asp335-Arg136, and Asp424-Arg428 are all important for thermostability, as the lengths of these bridges increased dramatically with increasing temperature. Here, for the first time, we have computationally identified the thermolabile regions of PhyA, and this information could be used to engineer novel thermostable phytases. Numerous homologous phytases of fungal as well as bacterial origin are known, and these homologs show high sequence similarity. Our findings could prove useful in attempts to increase the thermostability of homologous phytases via protein engineering.

Synergistic function of four novel thermostable glycoside hydrolases from a long-term enriched thermophilic methanogenic digester

In biofuel production from lignocellulose, low thermostability and product inhibition strongly restrict the enzyme activities and production process. Application of multiple thermostable glycoside hydrolases, forming an enzyme "cocktail", can result in a synergistic action and therefore improve production efficiency and reduce operational costs. Therefore, increasing enzyme thermostabilities and compatibility are important for the biofuel industry. In this study, we reported the screening, cloning and biochemical characterization of four novel thermostable lignocellulose hydrolases from a metagenomic library of a long-term dry thermophilic methanogenic digester community, which were highly compatible with optimal conditions and specific activities. The optimal temperatures of the four enzymes, β-xylosidase, xylanase, β-glucosidase, and cellulase ranged from 60 to 75°C, and over 80% residual activities were observed after 2 h incubation at 50°C. Mixtures of these hydrolases retained high residual synergistic activities after incubation with cellulose, xylan, and steam-exploded corncob at 50°C for 72 h. In addition, about 55% dry weight of steam-exploded corncob was hydrolyzed to glucose and xylose by the synergistic action of the four enzymes at 50°C for 48 h. This work suggested that since different enzymes from a same ecosystem could be more compatible, screening enzymes from a long-term enriching community could be a favorable strategy.

Structural insight into the thermostable NADP(+)-dependent meso-diaminopimelate dehydrogenase from Ureibacillus thermosphaericus

Crystal structures of the thermostable meso-diaminopimelate dehydrogenase (DAPDH) from Ureibacillus thermosphaericus were determined for the enzyme in the apo form and in complex with NADP(+) and N-tris(hydroxymethyl)methyl-2-aminoethanesulfonic acid. The main-chain coordinates of the enzyme showed notable similarity to those of Symbiobacterium thermophilum DAPDH. However, the subunit arrangement of U. thermosphaericus DAPDH (a dimer) was totally different from that of the S. thermophilum enzyme (a hexamer). Structural comparison with the dimeric enzyme from the mesophile Corynebacterium glutamicum revealed that the presence of large numbers of intrasubunit and intersubunit hydrophobic interactions, as well as the extensive formation of intersubunit ion-pair networks, were likely to be the main factors contributing to the higher thermostability of U. thermosphaericus DAPDH. This differs from S. thermophilum DAPDH, within which the unique hexameric assembly is likely to be responsible for its high thermostability. Analysis of the active site of U. thermosphaericus DAPDH revealed the key factors responsible for the marked difference in substrate specificity between DAPDH and the D-amino acid dehydrogenase recently created from DAPDH by introducing five point mutations [Akita et al. (2012). Biotechnol. Lett. 34, 1693-1699; 1701-1702].

Prediction of the determinants of thermal stability by linear discriminant analysis: the case of the glutamate dehydrogenase protein family

Little is known about the determinants of thermal stability in individual protein families. Most of the knowledge on thermostability comes, in fact, from comparative analyses between large, and heterogeneous, sets of thermo- and mesophilic proteins. Here, we present a multivariate statistical approach aimed to detect signature sequences for thermostability in a single protein family. It was applied to the glutamate dehydrogenase (GDH) family, which is a good model for investigating this peculiar process. The structure of GDH consists of six subunits, each of them organized into two domains. Formation of ion-pair networks on the surface of the protein subunits, or increase in the inter-subunit hydrophobic interactions, have been suggested as important factors for explaining stability at high temperatures. However, identification of the amino acid changes that are involved in this process still remains elusive. Our approach consisted of a linear discriminant analysis on a set of GDH sequences from Archaea and Bacteria (33 thermo- and 36 mesophilic GDHs). It led to detection of 3 amino acid clusters as the putative determinants of thermal stability. They were localized at the subunit interface or in close proximity to the binding site of the NAD(P)(+) coenzyme. Analysis within the clusters led to prediction of 8 critical amino acid sites. This approach could have a wide utility, in the ligth of the notion that each protein family seems to adopt its own strategy for achieving thermostability.

Thermostability of in vitro evolved Bacillus subtilis lipase A: a network and dynamics perspective

Proteins in thermophilic organisms remain stable and function optimally at high temperatures. Owing to their important applicability in many industrial processes, such thermostable proteins have been studied extensively, and several structural factors attributed to their enhanced stability. How these factors render the emergent property of thermostability to proteins, even in situations where no significant changes occur in their three-dimensional structures in comparison to their mesophilic counter-parts, has remained an intriguing question. In this study we treat Lipase A from Bacillus subtilis and its six thermostable mutants in a unified manner and address the problem with a combined complex network-based analysis and molecular dynamic studies to find commonality in their properties. The Protein Contact Networks (PCN) of the wild-type and six mutant Lipase A structures developed at a mesoscopic scale were analyzed at global network and local node (residue) level using network parameters and community structure analysis. The comparative PCN analysis of all proteins pointed towards important role of specific residues in the enhanced thermostability. Network analysis results were corroborated with finer-scale molecular dynamics simulations at both room and high temperatures. Our results show that this combined approach at two scales can uncover small but important changes in the local conformations that add up to stabilize the protein structure in thermostable mutants, even when overall conformation differences among them are negligible. Our analysis not only supports the experimentally determined stabilizing factors, but also unveils the important role of contacts, distributed throughout the protein, that lead to thermostability. We propose that this combined mesoscopic-network and fine-grained molecular dynamics approach is a convenient and useful scheme not only to study allosteric changes leading to protein stability in the face of negligible over-all conformational changes due to mutations, but also in other molecular networks where change in function does not accompany significant change in the network structure.

Effects of metal ions on stability and activity of hyperthermophilic pyrolysin and further stabilization of this enzyme by modification of a Ca2+-binding site

Pyrolysin is an extracellular subtilase produced by the marine hyperthermophilic archaeon Pyrococcus furiosus. This enzyme functions at high temperatures in seawater, but little is known about the effects of metal ions on the properties of pyrolysin. Here, we report that the supplementation of Na(+), Ca(2+), or Mg(2+) salts at concentrations similar to those in seawater destabilizes recombinant pyrolysin but leads to an increase in enzyme activity. The destabilizing effect of metal ions on pyrolysin appears to be related to the disturbance of surface electrostatic interactions of the enzyme. In addition, mutational analysis of two predicted high-affinity Ca(2+)-binding sites (Ca1 and Ca2) revealed that the binding of Ca(2+) is important for the stabilization of this enzyme. Interestingly, Asn substitutions at residues Asp818 and Asp820 of the Ca2 site, which is located in the C-terminal extension of pyrolysin, resulted in improvements in both enzyme thermostability and activity without affecting Ca(2+)-binding affinity. These effects were most likely due to the elimination of unfavorable electrostatic repulsion at the Ca2 site. Together, these results suggest that metal ions play important roles in modulating the stability and activity of pyrolysin.

Thermal adaptation of dihydrofolate reductase from the moderate thermophile Geobacillus stearothermophilus

The thermal melting temperature of dihydrofolate reductase from Geobacillus stearothermophilus (BsDHFR) is ~30 °C higher than that of its homologue from the psychrophile Moritella profunda. Additional proline residues in the loop regions of BsDHFR have been proposed to enhance the thermostability of BsDHFR, but site-directed mutagenesis studies reveal that these proline residues contribute only minimally. Instead, the high thermal stability of BsDHFR is partly due to removal of water-accessible thermolabile residues such as glutamine and methionine, which are prone to hydrolysis or oxidation at high temperatures. The extra thermostability of BsDHFR can be obtained by ligand binding, or in the presence of salts or cosolvents such as glycerol and sucrose. The sum of all these incremental factors allows BsDHFR to function efficiently in the natural habitat of G. stearothermophilus, which is characterized by temperatures that can reach 75 °C.

Structural basis of thermal stability of the tungsten cofactor synthesis protein MoaB from Pyrococcus furiosus

Molybdenum and tungsten cofactors share a similar pterin-based scaffold, which hosts an ene-dithiolate function being essential for the coordination of either molybdenum or tungsten. The biosynthesis of both cofactors involves a multistep pathway, which ends with the activation of the metal binding pterin (MPT) by adenylylation before the respective metal is incorporated. In the hyperthermophilic organism Pyrococcus furiosus, the hexameric protein MoaB (PfuMoaB) has been shown to catalyse MPT-adenylylation. Here we determined the crystal structure of PfuMoaB at 2.5 Å resolution and identified key residues of α3-helix mediating hexamer formation. Given that PfuMoaB homologues from mesophilic organisms form trimers, we investigated the impact on PfuMoaB hexamerization on thermal stability and activity. Using structure-guided mutagenesis, we successfully disrupted the hexamer interface in PfuMoaB. The resulting PfuMoaB-H3 variant formed monomers, dimers and trimers as determined by size exclusion chromatography. Circular dichroism spectroscopy as well as chemical cross-linking coupled to mass spectrometry confirmed a wild-type-like fold of the protomers as well as inter-subunits contacts. The melting temperature of PfuMoaB-H3 was found to be reduced by more than 15°C as determined by differential scanning calorimetry, thus demonstrating hexamerization as key determinant for PfuMoaB thermal stability. Remarkably, while a loss of activity at temperatures higher than 50°C was observed in the PfuMoaB-H3 variant, at lower temperatures, we determined a significantly increased catalytic activity. The latter suggests a gain in conformational flexibility caused by the disruption of the hexamerization interface.

A role for A-to-I RNA editing in temperature adaptation

A-to-I RNA editing can recode mRNAs, giving organisms the option to express diverse, functionally distinct protein isoforms. Here, we propose that RNA editing is inherently geared for temperature adaptation because it tends to recode to smaller, less stabilizing amino acids. Studies on how editing affects protein function support this idea.

Analysis of protein thermostability enhancing factors in industrially important thermus bacteria species

Elucidation of evolutionary factors that enhance protein thermostability is a critical problem and was the focus of this work on Thermus species. Pairs of orthologous sequences of T. scotoductus SA-01 and T. thermophilus HB27, with the largest negative minimum folding energy (MFE) as predicted by the UNAFold algorithm, were statistically analyzed. Favored substitutions of amino acids residues and their properties were determined. Substitutions were analyzed in modeled protein structures to determine their locations and contribution to energy differences using PyMOL and FoldX programs respectively. Dominant trends in amino acid substitutions consistent with differences in thermostability between orthologous sequences were observed. T. thermophilus thermophilic proteins showed an increase in non-polar, tiny, and charged amino acids. An abundance of alanine substituted by serine and threonine, as well as arginine substituted by glutamine and lysine was observed in T. thermophilus HB27. Structural comparison showed that stabilizing mutations occurred on surfaces and loops in protein structures.

The thermostability of two kinds of recombinant ∆6-fatty acid desaturase with different N-terminal sequence lengths in low temperature

Two recombinant Rhizopus stolonifer ∆6-fatty acid desaturase enzymes with different-length N-termini were cloned and expressed in Saccharomyces cerevisiae strain INVScl: LRsD6D begins with the sequence of the N-terminal of the R. stolonifer ∆6-fatty acid desaturase native, encoding a deduced polypeptide of 459 amino acids (M-S-T-L-D-R-Q-S-I-F-T-I-K-E-L-E-S-I-S-Q-R-I-H-D-G-D-E-E-A-M-K-F), whereas SRsD6D begins with the amino acid sequence of the predicted ORF, encoding a deduced polypeptide of 430 amino acids (M-K-F) and LRsD6D is longer than SRsD6D by 29 amino acids (M-S-T-L-D-R-Q-S-I-F-T-I-K-E-L-E-S-I-S-Q-R-I-H-D-G-D-E-E-A). Bioinformatic analysis characterized the two recombinant ∆6-fatty acid desaturase enzymes with different-length N-termini, including three conserved histidine-rich motifs, hydropathy profile, and a cytochrome b5-like domain in the N-terminus. When the coding sequence was expressed in S. cerevisiae strain INVScl, the coding produced ∆6-fatty acid desaturase activity exhibited by RsD6D, leading to a novel peak corresponding to γ-linolenic acid methyl ester standards, which was detected with the same retention time. The residual activity of LRsD6D was 74 % at 15 °C for 4 h and that of SRsD6D was 43 %. Purified recombinant LRsD6D was more stable than SRsD6D, indicating that the N-terminal extension, containing mostly hydrophobic residues, affected the overall stability of recombinant LRsD6D.

Glutamate dehydrogenases: the why and how of coenzyme specificity

NAD(+) and NADP(+), chemically similar and with almost identical standard oxidation-reduction potentials, nevertheless have distinct roles, NAD(+) serving catabolism and ATP generation whereas NADPH is the biosynthetic reductant. Separating these roles requires strict specificity for one or the other coenzyme for most dehydrogenases. In many organisms this holds also for glutamate dehydrogenases (GDH), NAD(+)-dependent for glutamate oxidation, NADP(+)-dependent for fixing ammonia. In higher animals, however, GDH has dual specificity. It has been suggested that GDH in mitochondria reacts only with NADP(H), the NAD(+) reaction being an in vitro artefact. However, contrary evidence suggests mitochondrial GDH not only reacts with NAD(+) but maintains equilibrium using the same pool as accessed by β-hydroxybutyrate dehydrogenase. Another complication is the presence of an energy-linked dehydrogenase driving NADP(+) reduction by NADH, maintaining the coenzyme pools at different oxidation-reduction potentials. Its coexistence with GDH makes possible a futile cycle, control of which is not yet properly explained. Structural studies show NAD(+)-dependent, NADP(+)-dependent and dual-specificity GDHs are closely related and a few site-directed mutations can reverse specificity. Specificity for NAD(+) or for NADP(+) has probably emerged repeatedly during evolution, using different structural solutions on different occasions. In various GDHs the P7 position in the coenzyme-binding domain plays a key role. However, whereas in other dehydrogenases an acidic P7 residue usually hydrogen bonds to the 2'- and 3'-hydroxyls, dictating NAD(+) specificity, among GDHs, depending on detailed conformation of surrounding residues, an acidic P7 may permit binding of NAD(+) only, NADP(+) only, or in higher animals both.

A multi-factors rational design strategy for enhancing the thermostability of Escherichia coli AppA phytase

Despite recent advances in our understanding of the importance of protein surface properties for protein thermostability,there are seldom studies on multi-factors rational design strategy, so a more scientific, simple and effective rational strategy is urgent for protein engineering. Here, we first attempted to use a three-factors rational design strategy combining three common structural features, protein flexibility, protein surface, and salt bridges. Escherichia coli AppA phytase was used as a model enzyme to improve its thermostability. Moreover, the structure and enzyme features of the thermostable mutants designed by our strategy were analyzed roundly. For the single mutants, two (Q206E and Y311K), in five exhibited thermostable property with a higher success rate of prediction (40 %). For the multiple mutants, the themostable sites were combined with another site, I427L, we obtained by directed evolution, Q206E/I427L, Y311K/I427L, and Q206E/Y311K/I427L, all exhibited thermostable property. The Y311K/I427L doubled thermostability (61.7 %, and was compared to 30.97 % after being heated at 80 °C for 10 min) and catalytic efficiency (4.46 was compared to 2.37) improved more than the wild-type AppA phytase almost without hampering catalytic activity. These multi-factors of rational design strategy can be applied practically as a thermostabilization strategy instead of the conventional single-factor approach.