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

The Correlated States Theory of the Hydrophobic Effect

Starting from Frank and Evans' "iceberg" model of hydrophobic hydration of small molecules (the "microscopic" hydrophobic effect, HE) published in 1945, much has been done with respect to understanding the nature of HE and elaborating a quantitative theory able to describe the thermodynamic profile (the "signature" of HE) for the large volume of experimental data accumulated to date. Generally, three sets of approaches addressing this issue were suggested, ranging from the approval of the central role of the water shell to the complete denial of its role, with the focus placed on the solute and its interactions with surrounding water. For this reason, some controversy is still present in understanding the fundamental nature of HE, even at the "microscopic" scale. Nevertheless, the general tendency of the past decade seems to have shifted toward a greater role of the water environment in determining the thermodynamic profile of HE, with a designated place for solute-water interactions as a fine-tuning of thermodynamic observables. In the present work, we developed a novel view on HE at the microscopic scale, appearing as a consequence of solute-water correlated translational and orientational vibration motion, emerging as a new property of hydrophobic hydration/solvation (the Correlated States Theory of HE). We built a fully analytical description of this process, which has enabled us to quantify the "signature" of HE for extended thermodynamic data sets without employing molecular simulations or any numerical functions in the core of the theory. As a consequence, our approach provides a self-consistent view on the known major experimental manifestations of HE across an extended temperature range, addresses some controversial issues existing to date, and creates a new augmentation to current knowledge. Most importantly, the suggested approach offers a paradigm shift from the currently dominating views on HE as a consequence of water-water interactions and the "excluded volume effect" toward the central role of solute-water interactions, and provides the first nonempirical proof of the validity of SASA-based computational models of hydrophobic hydration/solvation, which have been utilized on an empirical basis for more than 40 years.

Role of electrostatics in cold adaptation: A comparative study of eury- and stenopsychrophilic triose phosphate isomerase

Psychrophilic (cold-active) organisms have developed enzymes that facilitate sufficient metabolic activity at low temperatures to sustain life. This occurs through molecular adaptations that tend to increase protein flexibility at the expense of stability. However, psychrophiles also vary in their growth conditions. Eurypsychrophiles thrive over a wide temperature range and often prefer temperatures above 20 °C, while stenopsychrophiles grow optimally below 15 °C and are more narrowly adapted to cold temperatures. To elucidate differences between these two classes of enzymes, we here compare the stability and unfolding kinetics of two orthologues of the basal household enzyme triose phosphate isomerase, one from the stenopsychrophilic Antarctic permafrost bacterium Rhodonellum psychrophilum (sTPI) and the other from the eurypsychrophilic Greenland ikaite column bacterium Rhodococcus sp. JG-3 (eTPI). Remarkably, sTPI proved significantly more thermostable and resistant to chemical denaturation than its eurypsychrophilic counterpart, eTPI, in the absence of ionic components in solution, whereas inclusion of electrostatic screening agents in the form of sodium chloride or the charged denaturant guanidinium chloride largely cancelled out this difference. Thus, electrostatics play a prominent role in stabilizing the stenopsychrophilic sTPI, and a mandatory low-temperature growth environment does not preclude the development of considerable thermotolerance for individual enzymes. We were able to increase the thermostability of sTPI using an evolutionary machine learning model, which transferred several sTPI residues into the eTPI active site. While the stabilizing effect was modest, the combination of individual mutations was additive, underscoring the potential of combining multiple beneficial mutations to achieve enhanced enzyme properties.

Physics-based modeling in the new era of enzyme engineering

Enzyme engineering is entering a new era characterized by the integration of computational strategies. While bioinformatics and artificial intelligence methods have been extensively applied to accelerate the screening of function-enhancing mutants, physics-based modeling methods, such as molecular mechanics and quantum mechanics, are essential complements in many objectives. In this Perspective, we highlight how physics-based modeling will help the field of computational enzyme engineering reach its full potential by exploring current developments, unmet challenges and emerging opportunities for tool development.

Computer Simulations of the Temperature Dependence of Enzyme Reactions

In this review we discuss the development of methodology for calculating the temperature dependence and thermodynamic activation parameters for chemical reactions in solution and in enzymes, from computer simulations. We outline how this is done by combining the empirical valence bond method with molecular dynamics free energy simulations. In favorable cases it turns out that such simulations can even capture temperature optima for the catalytic rate. The approach turns out be very useful both for addressing questions regarding the roles of enthalpic and entropic effects in specific enzymes and also for attacking evolutionary problems regarding enzyme adaptation to different temperature regimes. In the latter case, we focus on cold-adaptation of enzymes from psychrophilic species and show how computer simulations have revealed the basic mechanisms behind such adaptation. Understanding these mechanisms also opens up the possibility of designing the temperature dependence, and we highlight a recent example of this.

High-Spin States of Manganese(III) Enable Robust Cold-Adapted Activity of MnO2 Nanozymes

Developing novel cold-adapted nanozymes and elucidating their mechanisms of action remains a great challenge. Inspired by natural oxidases that utilize high-spin and high-valent metal-oxygen intermediates to achieve high efficiency at low temperatures, in this study, a series of MnOx nanomaterials with varied valence and spin states are synthesized. The activity assay revealed that the oxygen vacancy-engineered ε-MnO2 nanozyme displayed excellent cold-adapted oxidase-like properties, and no observable activity loss is observed in the temperature range of -20 to 45 °C. The superior performance is attributed to the high-spin Mn(III)-O species coupled with its induced Jahn-Teller effect, which facilitates the dissociation and activation of oxygen at low temperatures. As a proof of concept, an excellent cold-adapted δ-MnO2 nanozyme can be obtained using Mn3O4 as the precursor by regulating the spin state of Mn(III). Moreover, a novel and effective degradation strategy for corn stalk at low temperature is built based on the robust cold-adapted oxidase-like activity of ε-MnO2. These results not only provide new insights for the rational design of cold-adapted nanozymes but also broaden the application of nanozymes in low-temperature industrial processes.

Entropy in catalyst dynamics under confinement

Entropy during the dynamic structural evolution of catalysts has a non-trivial influence on chemical reactions. Confinement significantly affects the catalyst dynamics and thus impacts the reactivity. However, a full understanding has not been clearly established. To investigate catalyst dynamics under confinement, we utilize the active learning scheme to effectively train machine learning potentials for computing free energies of catalytic reactions. The scheme enables us to compute the reaction free energies and entropies of O2 dissociation on Pt clusters with different sizes confined inside a carbon nanotube (CNT) at the timescale of tens of nanoseconds while keeping ab initio accuracy. We observe an entropic effect owing to liquid-to-solid phase transitions of clusters at finite temperatures. More importantly, the confinement effect enhances the structural dynamics of the cluster and leads to a lower melting temperature than those of the bare cluster and cluster outside the CNT, consequently facilitating the reaction to occur at lower temperatures and preventing the catalyst from forming unfavorable oxides. Our work reveals the important influence of confinement on structural dynamics, providing useful insight into entropy in dynamic catalysis.

Computational Analysis of Heat Capacity Effects in Protein-Ligand Binding

Heat capacity effects in protein-ligand binding as measured by calorimetric experiments have recently attracted considerable attention, particularly in the field of enzyme inhibitor design. A significant negative heat capacity change upon ligand binding implies a marked temperature dependence of the binding enthalpy, which is of high relevance for attempts to optimize protein-ligand interactions. In this work, we address the question of how well such heat capacity changes can be predicted by computer simulations. We examine a series of human thrombin inhibitors that all bind with ΔCp values of about -0.4 kcal/mol/K and calculate heat capacity changes from plain molecular dynamics simulations of the bound and free states of the enzyme and ligand. The results show that accurate ΔCp estimates within a few tenths of a kcal/mol/K of the experimental values can be obtained with this approach. This allows us to address the structural and energetic origin of the negative heat capacity changes for the thrombin inhibitors, and it is found that conformational equilibria of the free ligands in solution make a major contribution to the observed effect.

Computational Insights into SARS-CoV-2 Main Protease Mutations and Nirmatrelvir Efficacy: The Effects of P132H and P132H-A173V

Nirmatrelvir, a pivotal component of the oral antiviral Paxlovid for COVID-19, targets the SARS-CoV-2 main protease (Mpro) as a covalent inhibitor. Here, we employed combined computational methods to explore how the prevalent Omicron variant mutation P132H, alone and in combination with A173V (P132H-A173V), affects nirmatrelvir's efficacy. Our findings suggest that P132H enhances the noncovalent binding affinity of Mpro for nirmatrelvir, whereas P132H-A173V diminishes it. Although both mutants catalyze the rate-limiting step more efficiently than the wild-type (WT) Mpro, P132H slows the overall rate of covalent bond formation, whereas P132H-A173V accelerates it. Comprehensive analysis of noncovalent and covalent contributions to the overall binding free energy of the covalent complex suggests that P132H likely enhances Mpro sensitivity to nirmatrelvir, while P132H-A173V may confer resistance. Per-residue decompositions of the binding and activation free energies pinpoint key residues that significantly affect the binding affinity and reaction rates, revealing how the mutations modulate these effects. The mutation-induced conformational perturbations alter drug-protein local contact intensities and the electrostatic preorganization of the protein, affecting noncovalent binding affinity and the stability of key reaction states, respectively. Our findings inform the mechanisms of nirmatrelvir resistance and sensitivity, facilitating improved drug design and the detection of resistant strains.

Interplay of structural preorganization and conformational sampling in UDP-glucuronic acid 4-epimerase catalysis

Understanding enzyme catalysis as connected to protein motions is a major challenge. Here, based on temperature kinetic studies combined with isotope effect measurements, we obtain energetic description of C-H activation in NAD-dependent UDP-glucuronic acid C4 epimerase. Approach from the ensemble-averaged ground state (GS) to the transition state-like reactive conformation (TSRC) involves, alongside uptake of heat (Δ H ‡  = 54 kJ mol-1), significant loss in entropy ( - T Δ S ‡  = 20 kJ mol-1; 298 K) and negative activation heat capacity (Δ C p ‡  = -0.64 kJ mol-1 K-1). Thermodynamic changes suggest the requirement for restricting configurational freedom at the GS to populate the TSRC. Enzyme variants affecting the electrostatic GS preorganization reveal active-site interactions important for precise TSRC sampling and H-transfer. Collectively, our study captures thermodynamic effects associated with TSRC sampling and establishes rigid positioning for C-H activation in an enzyme active site that requires conformational flexibility in fulfillment of its natural epimerase function.

Exploring the Antarctic aminopeptidase P from Pseudomonas sp. strain AMS3 through structural analysis and molecular dynamics simulation

Aminopeptidase P (APPro) is a crucial metalloaminopeptidase involved in amino acid cleavage from peptide N-termini, playing essential roles as versatile biocatalysts with applications ranging from pharmaceuticals to industrial processes. Despite acknowledging its potential for catalysis in lower temperatures, detailed molecular basis and biotechnological implications in cold environments are yet to be explored. Therefore, this research aims to investigate the molecular mechanisms underlying the cold-adapted characteristics of APPro from Pseudomonas sp. strain AMS3 (AMS3-APPro) through a detailed analysis of its structure and dynamics. In this study, structure analysis and molecular dynamics (MD) simulation of a predicted model of AMS3-APPro has been performed at different temperatures to assess structural flexibility and thermostability across a temperature range of 0-60 °C over 100 ns. The MD simulation results revealed that the structure were able to remain stable at low temperatures. Increased temperatures present a potential threat to the overall stability of AMS3-APPro by disrupting the intricate hydrogen bond networks crucial for maintaining structural integrity, thereby increasing the likelihood of protein unfolding. While the metal binding site at the catalytic core exhibits resilience at higher temperatures, highlighting its local structural integrity, the overall enzyme structure undergoes fluctuations and potential denaturation. This extensive structural instability surpasses the localized stability observed at the metal binding site. Consequently, these assessments offer in-depth understanding of the cold-adapted characteristics of AMS3-APPro, highlighting its capability to uphold its native conformation and stability in low-temperature environments. In summary, this research provides valuable insights into the cold-adapted features of AMS3-APPro, suggesting its efficient operation in low thermal conditions, particularly relevant for potential biotechnological applications in cold environments.

Why Do Empirical Valence Bond Simulations Yield Accurate Arrhenius Plots?

Computer simulations of the temperature dependence of enzyme reactions using the empirical valence bond (EVB) method have proven to give very accurate results in terms of the thermodynamic activation parameters. Here, we analyze the reasons for why such simulations are able to correctly capture activation enthalpies and entropies and how sensitive these quantities are to parametrization of the reactive potential energy function. We examine first the solution reference reaction for the enzyme ketosteroid isomerase, which corresponds to the acetate catalyzed deprotonation of the steroid in water. The experimentally determined activation parameters for this reaction turn out to be remarkably well reproduced by the calculations. By modifying the EVB potential so that the activation and reaction free energies become significantly shifted, we show that the activation entropy is basically invariant to such changes and that ΔS⧧ is instead determined by the specific mixture of the underlying force fields in the transition state region. The coefficients of this mixture do not change appreciably when the EVB potential is modified within reasonable limits, and hence, the estimate of ΔS⧧ becomes very robust. This is further verified by examining a more complex concerted hydride and proton transfer reaction in the enzyme hydroxybutyrate dehydrogenase.

Enthalpies and entropies of hydration from Monte Carlo simulations

The changes in free energy, enthalpy, and entropy for transfer of a solute from the gas phase into solution are the fundamental thermodynamic quantities that characterize the solvation process. Owing to the development of methods based on free-energy perturbation theory, computation of free energies of solvation has become routine in conjunction with Monte Carlo (MC) statistical mechanics and molecular dynamics (MD) simulations. Computation of the enthalpy change and by inference the entropy change is more challenging. Two methods are considered in this work corresponding to direct averaging for the solvent and solution and to computing the temperature derivative of the free energy in the van't Hoff approach. The application is for neutral organic solutes in TIP4P water using long MC simulations to improve precision. Definitive results are also provided for pure TIP4P water. While the uncertainty in computed free energies of hydration is ca. 0.05 kcal mol-1, it is ca. 0.4 kcal mol-1 for the enthalpy changes from either van't Hoff plots or the direct method with sampling for 5 billion MC configurations. Partial molar volumes of hydration are also computed by the direct method; they agree well with experimental data with an average deviation of 3 cm3 mol-1. In addition, the results permit breakdown of the errors in the free energy changes from the OPLS-AA force field into their enthalpic and entropic components. The excess hydrophobicity of organic solutes is enthalpic in origin.

A Manganese-Based Metal-Organic Framework as a Cold-Adapted Nanozyme

The development of cold-adapted enzymes with high efficiency and good stability is an advanced strategy to overcome the limitations of catalytic medicine in low and cryogenic temperatures. In this work, inspired by natural enzymes, a novel cold-adapted nanozyme based on a manganese-based nanosized metal-organic framework (nMnBTC) is designed and synthesized. The nMnBTC as an oxidase mimetic not only exhibits excellent activity at 0 °C, but also presents almost no observable activity loss as the temperature is increased to 45 °C. This breaks the traditional recognition that enzymes show maximum activity only under specific psychrophilic or thermophilic condition. The superior performance of nMnBTC as a cold-adapted nanozyme can be attributed to its high-catalytic efficiency at low temperature, good substrate affinity, and flexible conformation. Based on the robust performance of nMnBTC, a low-temperature antiviral strategy is developed to inactivate influenza virus H1N1 even at -20 °C. These results not only provide an important guide for the rational design of highly efficient artificial cold-adapted enzymes, but also pave a novel way for biomedical application in cryogenic fields.

Snapshots of the Reaction Coordinate of a Thermophilic 2'-Deoxyribonucleoside/ribonucleoside Transferase

Nucleosides are ubiquitous to life and are required for the synthesis of DNA, RNA, and other molecules crucial for cell survival. Despite the notoriously difficult organic synthesis of nucleosides, 2'-deoxynucleoside analogues can interfere with natural DNA replication and repair and are successfully employed as anticancer, antiviral, and antimicrobial compounds. Nucleoside 2'-deoxyribosyltransferase (dNDT) enzymes catalyze transglycosylation via a covalent 2'-deoxyribosylated enzyme intermediate with retention of configuration, having applications in the biocatalytic synthesis of 2'-deoxynucleoside analogues in a single step. Here, we characterize the structure and function of a thermophilic dNDT, the protein from Chroococcidiopsis thermalis (CtNDT). We combined enzyme kinetics with structural and biophysical studies to dissect mechanistic features in the reaction coordinate, leading to product formation. Bell-shaped pH-rate profiles demonstrate activity in a broad pH range of 5.5-9.5, with two very distinct pKa values. A pronounced viscosity effect on the turnover rate indicates a diffusional step, likely product (nucleobase1) release, to be rate-limiting. Temperature studies revealed an extremely curved profile, suggesting a large negative activation heat capacity. We trapped a 2'-fluoro-2'-deoxyarabinosyl-enzyme intermediate by mass spectrometry and determined high-resolution structures of the protein in its unliganded, substrate-bound, ribosylated, 2'-difluoro-2'-deoxyribosylated, and in complex with probable transition-state analogues. We reveal key features underlying (2'-deoxy)ribonucleoside selection, as CtNDT can also use ribonucleosides as substrates, albeit with a lower efficiency. Ribonucleosides are the building blocks of RNA and other key intracellular metabolites participating in energy and metabolism, expanding the scope of use of CtNDT in biocatalysis.

Exploring the molecular mechanism of cold-adaption of an alkaline protease mutant by molecular dynamics simulations and residue interaction network

Psychrophilic proteases have attracted enormous attention in past decades, due to their high catalytic activity at low temperatures in a wide range of industrial processes, especially in the detergent and leather industries. Among them, H5 is an alkaline protease mutant, which featuring psychrophilic-like behavior, but the reasons that H5 with higher activity at low temperatures are still poorly understood. Herein, the molecular dynamics (MD) simulations combined with residue interaction network (RIN) were utilized to investigate the mechanisms of the cold-adaption of mutant H5. The results demonstrated that two loops involved in the substrate binding G100-S104 and S125-S129 in H5 had higher mobility, and the distance enlargement between the two loops modulated the substrate's accessibility compared with wild type counterpart. Besides, H5 enhanced conformational flexibility by weakening salt bridges and increasing interaction with the solvent. In particular, the absence of Lys251-Asp197-Arg247 salt bridge network may contribute to the structural mobility. Based on the free energy landscape and molecular mechanics Poisson-Boltzmann surface area of the wild type and H5, it was elucidated that H5 possessed a large population of interconvertible conformations, resulting in the weaker substrate binding free energy. The calculated RIN topology parameters such as the average degree, average cluster coefficient, and average path length further verified that the mutant H5 attenuated residue-to-residue interactions. The investigation of the mechanisms by which how the residue mutation affects the stability and activity of enzymes provides a theoretical basis for the development of cold-adapted protease.

Helping proteins come in from the cold: 5 burning questions about cold-active enzymes

Enzymes from psychrophilic (cold-loving) organisms have attracted considerable interest over the past decades for their potential in various low-temperature industrial processes. However, we still lack large-scale commercialization of their activities. Here, we review their properties, limitations and potential. Our review is structured around answers to 5 central questions: 1. How do cold-active enzymes achieve high catalytic rates at low temperatures? 2. How is protein flexibility connected to cold-activity? 3. What are the sequence-based and structural determinants for cold-activity? 4. How does the thermodynamic stability of psychrophilic enzymes reflect their cold-active capabilities? 5. How do we effectively identify novel cold-active enzymes, and can we apply them in an industrial context? We conclude that emerging screening technologies combined with big-data handling and analysis make it reasonable to expect a bright future for our understanding and exploitation of cold-active enzymes.

Functional activity of E. coli RNase R in the Antarctic Pseudomonas syringae Lz4W

BACKGROUND

In Antarctic P. syringae RNase R play an essential role in the processing of 16S and 5S rRNA, thereby playing an important role in cold-adapted growth of the bacterium. This study is focused on deciphering the in vivo functional activity of mesophilic exoribonuclease R and its catalytic domain (RNB) in an evolutionary distant psychrophilic bacterium Pseudomonas syringae Lz4W.

RESULTS

Our results confirm that E. coli RNase R complemented the physiological functions of the psychrophilic bacterium P. syringae RNase R and rescued the cold-sensitive phenotype of Pseudomonas syringae ∆rnr mutant. More importantly, the catalytic domain (RNB) of the E. coli RNase R is also capable of alleviating the cold-sensitive growth defects of ∆rnr mutant as seen with the catalytic domain (RNB) of the P. syringae enzyme. The Catalytic domain of E. coli RNase R was less efficient than the Catalytic domain of P. syringae RNase R in rescuing the cold-sensitive growth of ∆rnr mutant at 4°C, as the ∆rnr expressing the RNBEc (catalytic domain of E. coli RNase R) displayed longer lag phase than the RNBPs (Catalytic domain of P. syringae RNase R) complemented ∆rnr mutant at 4°C. Altogether it appears that the E. coli RNase R and P. syringae RNase R are functionally exchangeable for the growth requirements of P. syringae at low temperature (4°C). Our results also confirm that in P. syringae the requirement of RNase R for supporting the growth at 4°C is independent of the degradosomal complex.

CONCLUSION

E. coli RNase R (RNase REc) rescues the cold-sensitive phenotype of the P. syringae Δrnr mutant. Similarly, the catalytic domain of E. coli RNase R (RNBEc) is also capable of supporting the growth of Δrnr mutant at low temperatures. These findings have a vast scope in the design and development of low-temperature-based expression systems.

Psychrophilic enzymes: strategies for cold-adaptation

Psychrophilic organisms thriving at near-zero temperatures synthesize cold-adapted enzymes to sustain cell metabolism. These enzymes have overcome the reduced molecular kinetic energy and increased viscosity inherent to their environment and maintained high catalytic rates by development of a diverse range of structural solutions. Most commonly, they are characterized by a high flexibility coupled with an intrinsic structural instability and reduced substrate affinity. However, this paradigm for cold-adaptation is not universal as some cold-active enzymes with high stability and/or high substrate affinity and/or even an unaltered flexibility have been reported, pointing to alternative adaptation strategies. Indeed, cold-adaptation can involve any of a number of a diverse range of structural modifications, or combinations of modifications, depending on the enzyme involved, its function, structure, stability, and evolutionary history. This paper presents the challenges, properties, and adaptation strategies of these enzymes.

Computational design of the temperature optimum of an enzyme reaction

Cold-adapted enzymes are characterized both by a higher catalytic activity at low temperatures and by having their temperature optimum down-shifted, compared to mesophilic orthologs. In several cases, the optimum does not coincide with the onset of protein melting but reflects some other type of inactivation. In the psychrophilic α-amylase from an Antarctic bacterium, the inactivation is thought to originate from a specific enzyme-substrate interaction that breaks around room temperature. Here, we report a computational redesign of this enzyme aimed at shifting its temperature optimum upward. A set of mutations designed to stabilize the enzyme-substrate interaction were predicted by computer simulations of the catalytic reaction at different temperatures. The predictions were verified by kinetic experiments and crystal structures of the redesigned α-amylase, showing that the temperature optimum is indeed markedly shifted upward and that the critical surface loop controlling the temperature dependence approaches the target conformation observed in a mesophilic ortholog.

Calculation of Heat Capacity Changes in Enzyme Catalysis and Ligand Binding

It has been suggested that heat capacity changes in enzyme catalysis may be the underlying reason for temperature optima that are not related to unfolding of the enzyme. If this were to be a common phenomenon, it would have major implications for our interpretation of enzyme kinetics. In most cases, the support for the possible existence of a nonzero (negative) activation heat capacity, however, only relies on fitting such a kinetic model to experimental data. It is therefore of fundamental interest to try to use computer simulations to address this issue. One way is simply to calculate the temperature dependence of the activation free energy and determine whether the relationship is linear or not. An alternative approach is to calculate the absolute heat capacities of the reactant and transition states from plain molecular dynamics simulations using either the temperature derivative or fluctuation formula for the enthalpy. Here, we examine these different approaches for a designer enzyme with a temperature optimum that is not caused by unfolding. Benchmark calculations for the heat capacity of liquid water are first carried out using different thermostats. It is shown that the derivative formula for the heat capacity is generally the most robust and insensitive to the thermostat used and its parameters. The enzyme calculations using this method give results in agreement with direct calculations of activation free energies and show no sign of a negative activation heat capacity. We also provide a simple scheme for the calculation of binding heat capacity changes, which is of clear interest in ligand design, and demonstrate it for substrate binding to the designer enzyme. Neither in that case do the simulations predict any negative heat capacity change.

Structure and Mechanism of a Cold-Adapted Bacterial Lipase

The structural origin of enzyme cold-adaptation has been the subject of considerable research efforts in recent years. Comparative studies of orthologous mesophilic-psychrophilic enzyme pairs found in nature are an obvious strategy for solving this problem, but they often suffer from relatively low sequence identity of the enzyme pairs. Small bacterial lipases adapted to distinctly different temperatures appear to provide an excellent model system for these types of studies, as they may show a very high degree of sequence conservation. Here, we report the first crystal structures of lipase A from the psychrophilic bacterium Bacillus pumilus, which confirm the high structural similarity to the mesophilic Bacillus subtilis enzyme, as indicated by their 81% sequence identity. We further employ extensive QM/MM calculations to delineate the catalytic reaction path and its energetics. The computational prediction of a rate-limiting deacylation step of the enzymatic ester hydrolysis reaction is verified by stopped-flow experiments, and steady-state kinetics confirms the psychrophilic nature of the B. pumilus enzyme. These results provide a useful benchmark for examining the structural basis of cold-adaptation and should now make it possible to disentangle the effects of the 34 mutations between the two enzymes on catalytic properties and thermal stability.

The Activation Parameters of a Cold-Adapted Short Chain Dehydrogenase Are Insensitive to Enzyme Oligomerization

The structural principles of enzyme cold adaptation are of fundamental interest both for understanding protein evolution and for biotechnological applications. It has become clear in recent years that structural flexibility plays a major role in tuning enzyme activity at low temperatures, which is reflected by characteristic changes in the thermodynamic activation parameters for psychrophilic enzymes, compared to those of mesophilic and thermophilic ones. Hence, increased flexibility of the enzyme surface has been shown to lead to a lower enthalpy and a more negative entropy of activation, which leads to higher activity in the cold. This immediately raises the question of how enzyme oligomerization affects the temperature dependence of catalysis. Here, we address this issue by computer simulations of the catalytic reaction of a cold-adapted bacterial short chain dehydrogenase in different oligomeric states. Reaction free energy profiles are calculated at different temperatures for the tetrameric, dimeric, and monomeric states of the enzyme, and activation parameters are obtained from the corresponding computational Arrhenius plots. The results show that the activation free energy, enthalpy, and entropy are remarkably insensitive to the oligomeric state, leading to the conclusion that assembly of the subunit interfaces does not compromise cold adaptation, even though the mobilities of interfacial residues are indeed affected.

Double mutations far from the active site affect cold activity in an Antarctic halophilic β-galactosidase

The Antarctic haloarchaeon, Halorubrum lacusprofundi, contains a polyextremophilic family 42 β-galactosidase, which we are using as a model for cold-active enzymes. Divergent amino acid residues in this 78 kDa protein were identified through comparative genomics and hypothesized to be important for cold activity. Six amino acid residues were previously mutated and five were shown by steady-state kinetic analysis to have altered temperature-dependent catalytic activity profiles via effects on K_m and/or k_cat compared to the wild-type enzyme. In this follow-up study, double-mutated enzymes were constructed and tested for temperature effects, including two new tandem residue pairs (N180T/A181T and T383A/S384A), and pairwise combination of the single residue mutations (N251D, F387L, I476V, and V482L). All double-mutated enzymes were found to be more catalytically active at moderate and/or less active at colder temperatures than wild-type, with both K_m and k_cat effects observed for the two tandem mutations. For pairwise combinations, a K_m effect was seen when the surface exposed F387L mutation located in a domain A TIM barrel α helix 19 Å from the active site was combined with two internal residues, N251D or V482L. When another surface exposed mutation I476V located in a coiled region of domain B 25 Å from the active site was paired with N251D or V482L, a k_cat effect was observed. These results indicate that temperature-dependent kinetic effects may be complex and subtle and are mediated by a combination of a small number of residues distant from the active site via changes to the hydration shell and/or perturbation of internal packing.

Structural and functional adaptation in extremophilic microbial α-amylases

Maintaining stable native conformation of a protein under a given ecological condition is the prerequisite for survival of organisms. Extremophilic bacteria and archaea have evolved to adapt under extreme conditions of temperature, pH, salt, and pressure. Molecular adaptations of proteins under these conditions are essential for their survival. These organisms have the capability to maintain stable, native conformations of proteins under extreme conditions. The enzymes produced by the extremophiles are also known as extremozyme, which are used in several industries. Stability and functionality of extremozymes under varying temperature, pH, and solvent conditions are the most desirable requirement of industry. α-Amylase is one of the most important enzymes used in food, pharmaceutical, textile, and detergent industries. This enzyme is produced by diverse microorganisms including various extremophiles. Therefore, understanding its stability is important from fundamental as well as an applied point of view. Each class of extremophiles has a distinctive set of dominant non-covalent interactions which are important for their stability. Static information obtained by comparative analysis of amino acid sequence and atomic resolution structure provides information on the prevalence of particular amino acids or a group of non-covalent interactions. Protein folding studies give the information about thermodynamic and kinetic stability in order to understand dynamic aspect of molecular adaptations. In this review, we have summarized information on amino acid sequence, structure, stability, and adaptability of α-amylases from different classes of extremophiles.

Chemical Mapping Exposes the Importance of Active Site Interactions in Governing the Temperature Dependence of Enzyme Turnover

Uncovering the role of global protein dynamics in enzyme turnover is needed to fully understand enzyme catalysis. Recently, we have demonstrated that the heat capacity of catalysis, ΔC _P ^‡, can reveal links between the protein free energy landscape, global protein dynamics, and enzyme turnover, suggesting that subtle changes in molecular interactions at the active site can affect long-range protein dynamics and link to enzyme temperature activity. Here, we use a model promiscuous enzyme (glucose dehydrogenase from Sulfolobus solfataricus) to chemically map how individual substrate interactions affect the temperature dependence of enzyme activity and the network of motions throughout the protein. Utilizing a combination of kinetics, red edge excitation shift (REES) spectroscopy, and computational simulation, we explore the complex relationship between enzyme-substrate interactions and the global dynamics of the protein. We find that changes in ΔC _P ^‡ and protein dynamics can be mapped to specific substrate-enzyme interactions. Our study reveals how subtle changes in substrate binding affect global changes in motion and flexibility extending throughout the protein.

Molecular Characterization of a Mesophilic Cellobiose 2-Epimerase That Maintains a High Catalytic Efficiency at Low Temperatures

Cellobiose 2-epimerase (CE) can catalyze bioconversion of lactose to its prebiotic derivative epilactose. The catalytic property of a novel CE from Treponema brennaborense (Trbr-CE) was investigated. Trbr-CE showed the highest catalytic efficiency of epimerization toward lactose among all of the previously reported CEs. This enzyme's specific activity could reach as high as 208.5 ± 5.3 U/mg at its optimum temperature, which is 45 °C. More importantly, this enzyme demonstrated a considerably high activity at low temperatures, suggesting Trbr-CE as a promising enzyme for industrial low-temperature production of epilactose. This structurally flexible enzyme exhibited a comparatively high binding affinity toward substrates, which was confirmed by both experimental verification and computational analysis. Molecular dynamics (MD) simulations and binding free energy calculations were applied to provide insights into molecular recognition upon temperature changes. Compared with thermophilic CEs, Trbr-CE presents a more negative enthalpy change and a higher entropy change when the temperature drops.

Diversity, Ecological Role and Biotechnological Potential of Antarctic Marine Fungi

The Antarctic Ocean is one of the most remote and inaccessible environments on our planet and hosts potentially high biodiversity, being largely unexplored and undescribed. Fungi have key functions and unique physiological and morphological adaptations even in extreme conditions, from shallow habitats to deep-sea sediments. Here, we summarized information on diversity, the ecological role, and biotechnological potential of marine fungi in the coldest biome on Earth. This review also discloses the importance of boosting research on Antarctic fungi as hidden treasures of biodiversity and bioactive molecules to better understand their role in marine ecosystem functioning and their applications in different biotechnological fields.

A meta-analysis of the activity, stability, and mutational characteristics of temperature-adapted enzymes

Understanding the characteristics that define temperature-adapted enzymes has been a major goal of extremophile enzymology in recent decades. In the present study, we explore these characteristics by comparing psychrophilic, mesophilic, and thermophilic enzymes. Through a meta-analysis of existing data, we show that psychrophilic enzymes exhibit a significantly larger gap (T_g) between their optimum and melting temperatures compared with mesophilic and thermophilic enzymes. These results suggest that T_g may be a useful indicator as to whether an enzyme is psychrophilic or not and that models of psychrophilic enzyme catalysis need to account for this gap. Additionally, by using predictive protein stability software, HoTMuSiC and PoPMuSiC, we show that the deleterious nature of amino acid substitutions to protein stability increases from psychrophiles to thermophiles. How this ultimately affects the mutational tolerance and evolutionary rate of temperature adapted organisms is currently unknown.

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.

Cold-Active β-Galactosidases: Insight into Cold Adaption Mechanisms and Biotechnological Exploitation

β-galactosidases (EC 3.2.1.23) catalyze the hydrolysis of β-galactosidic bonds in oligosaccharides and, under certain conditions, transfer a sugar moiety from a glycosyl donor to an acceptor. Cold-active β-galactosidases are identified in microorganisms endemic to permanently low-temperature environments. While mesophilic β-galactosidases are broadly studied and employed for biotechnological purposes, the cold-active enzymes are still scarcely explored, although they may prove very useful in biotechnological processes at low temperature. This review covers several issues related to cold-active β-galactosidases, including their classification, structure and molecular mechanisms of cold adaptation. Moreover, their applications are discussed, focusing on the production of lactose-free dairy products as well as on the valorization of cheese whey and the synthesis of glycosyl building blocks for the food, cosmetic and pharmaceutical industries.

Hexamerization and thermostability emerged very early during geranylgeranylglyceryl phosphate synthase evolution

A large number of archaea live in hyperthermophilic environments. In consequence, their proteins need to adopt to these harsh conditions, including the enzymes that catalyze the synthesis of their membrane ether lipids. The enzyme that catalyzes the formation of the first ether bond in these lipids, geranylgeranylglyceryl phosphate synthase (GGGPS), exists as a hexamer in many hyperthermophilic archaea, and a recent study suggested that hexamerization serves for a fine‐tuning of the flexibility – stability trade‐off under hyperthermophilic conditions. We have recently reconstructed the sequences of ancestral group II GGGPS enzymes and now present a detailed biochemical characterization of nine of these predecessors, which allowed us to trace back the evolution of hexameric GGGPS and to draw conclusions about the properties of extant GGGPS branches that were not accessible to experiments up to now. Almost all ancestral GGGPS proteins formed hexamers, which demonstrates that hexamerization is even more widespread among the GGGPS family than previously assumed. Furthermore, all experimentally studied ancestral proteins showed high thermostability. Our results indicate that the hexameric oligomerization state and thermostability were present very early during the evolution of group II GGGPS, while the fine tuning of the flexibility – stability trade‐off developed very late, independent of the emergence of hexamerization.

Taking Advantage of Promiscuity of Cold-Active Enzymes

Cold-active enzymes increase their catalytic efficiency at low-temperature, introducing structural flexibility at or near the active sites. Inevitably, this feat seems to be accompanied by lower thermal stability. These characteristics have made cold-active enzymes into attractive targets for the industrial applications, since they could reduce the energy cost in the reaction, attenuate side-reactions, and simply be inactivated. In addition, the increased structural flexibility could result in broad substrate specificity for various non-native substrates, which is called substrate promiscuity. In this perspective, we deal with a less addressed aspect of cold-active enzymes, substrate promiscuity, which has enormous potential for semi-synthesis or enzymatic modification of fine chemicals and drugs. Further structural and directed-evolutional studies on substrate promiscuity of cold-active enzymes will provide a new workhorse in white biotechnology.

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.

The Energetic Origin of Different Catalytic Activities in Temperature-Adapted Trypsins

Psychrophilic enzymes were always observed to have higher catalytic activity (k cat) than their mesophilic homologs at room temperature, while the origin of this phenomenon remains obscure. Here, we used two different temperature-adapted trypsins, the psychrophilic Atlantic cod trypsin (ACT) and the mesophilic bovine trypsin (BT), as a model system to explore the energetic origin of their different catalytic activities using computational methods. The results reproduce the characteristic changing trends in the activation free energy, activation enthalpy, and activation entropy between the psychrophilic and mesophilic enzymes, where, in particular, the slightly decreased activation free energy of ACT is determined by its considerably reduced activation enthalpy rather than by its more negative activation entropy compared to BT. The calculated electrostatic contributions to the solvation free energies in the reactant state/ground sate (RS/GS) and transition state (TS) show that, going from BT to ACT, the TS stabilization has a predominant effect over the RS stabilization on lowering the activation enthalpy of ACT. Comparison between the solvation energy components reveals a more optimized electrostatic preorganization to the TS in ACT, which provides a larger stabilization to the TS through reducing the reorganization energy, thus resulting in the lower activation enthalpy and hence lower activation free energy of ACT. Thus, it can be concluded that it is the difference in the protein electrostatic environment, and hence its different stabilizing effects on the TS, that brings about the different catalytic activities of different temperature-adapted trypsins.