Different roles of electrostatics in heat and in cold: adaptation by citrate synthase

Cold-Active Lipases and Esterases: A Review on Recombinant Overexpression and Other Essential Issues

Cold environments characterised by diverse temperatures close to or below the water freezing point dominate about 80% of the Earth's biosphere. One of the survival strategies adopted by microorganisms living in cold environments is their expression of cold-active enzymes that enable them to perform an efficient metabolic flux at low temperatures necessary to thrive and reproduce under those constraints. Cold-active enzymes are ideal biocatalysts that can reduce the need for heating procedures and improve industrial processes' quality, sustainability, and cost-effectiveness. Despite their wide applications, their industrial usage is still limited, and the major contributing factor is the lack of complete understanding of their structure and cold adaptation mechanisms. The current review looked at the recombinant overexpression, purification, and recent mechanism of cold adaptation, various approaches for purification, and three-dimensional (3D) crystal structure elucidation of cold-active lipases and esterase.

Coevolution of the Ess1-CTD axis in polar fungi suggests a role for phase separation in cold tolerance

Most of the world's biodiversity lives in cold (-2° to 4°C) and hypersaline environments. To understand how cells adapt to such conditions, we isolated two key components of the transcription machinery from fungal species that live in extreme polar environments: the Ess1 prolyl isomerase and its target, the carboxy-terminal domain (CTD) of RNA polymerase II. Polar Ess1 enzymes are conserved and functional in the model yeast, Saccharomyces cerevisiae. By contrast, polar CTDs diverge from the consensus (YSPTSPS)26 and are not fully functional in S. cerevisiae. These CTDs retain the critical Ess1 Ser-Pro target motifs, but substitutions at Y1, T4, and S7 profoundly affected their ability to undergo phase separation in vitro and localize in vivo. We propose that environmentally tuned phase separation by the CTD and other intrinsically disordered regions plays an adaptive role in cold tolerance by concentrating enzymes and substrates to overcome energetic barriers to metabolic activity.

Structure and in silico simulations of a cold-active esterase reveals its prime cold-adaptation mechanism

Here we determined the structure of a cold active family IV esterase (EstN7) cloned from Bacillus cohnii strain N1. EstN7 is a dimer with a classical α/β hydrolase fold. It has an acidic surface that is thought to play a role in cold-adaption by retaining solvation under changed water solvent entropy at lower temperatures. The conformation of the functionally important cap region is significantly different to EstN7's closest relatives, forming a bridge-like structure with reduced helical content providing greater access to the active site through more than one substrate access tunnel. However, dynamics do not appear to play a major role in cold adaption. Molecular dynamics at different temperatures, rigidity analysis, normal mode analysis and geometric simulations of motion confirm the flexibility of the cap region but suggest that the rest of the protein is largely rigid. Rigidity analysis indicates the distribution of hydrophobic tethers is appropriate to colder conditions, where the hydrophobic effect is weaker than in mesophilic conditions due to reduced water entropy. Thus, it is likely that increased substrate accessibility and tolerance to changes in water entropy are important for of EstN7's cold adaptation rather than changes in dynamics.

Computational Study on Temperature Driven Structure-Function Relationship of Polysaccharide Producing Bacterial Glycosyl Transferase Enzyme

Glycosyltransferase (GTs) is a wide class of enzymes that transfer sugar moiety, playing a key role in the synthesis of bacterial exopolysaccharide (EPS) biopolymer. In recent years, increased demand for bacterial EPSs has been observed in pharmaceutical, food, and other industries. The application of the EPSs largely depends upon their thermal stability, as any industrial application is mainly reliant on slow thermal degradation. Keeping this in context, EPS producing GT enzymes from three different bacterial sources based on growth temperature (mesophile, thermophile, and hyperthermophile) are considered for in silico analysis of the structural–functional relationship. From the present study, it was observed that the structural integrity of GT increases significantly from mesophile to thermophile to hyperthermophile. In contrast, the structural plasticity runs in an opposite direction towards mesophile. This interesting temperature-dependent structural property has directed the GT–UDP-glucose interactions in a way that thermophile has finally demonstrated better binding affinity (−5.57 to −10.70) with an increased number of hydrogen bonds (355) and stabilizing amino acids (Phe, Ala, Glu, Tyr, and Ser). The results from this study may direct utilization of thermophile-origin GT as best for industrial-level bacterial polysaccharide production.

X-ray crystallographic structural studies of α-amylase I from Eisenia fetida

The earthworm Eisenia fetida possesses several cold-active enzymes, including α-amylase, β-glucanase and β-mannanase. E. fetida possesses two isoforms of α-amylase (Ef-Amy I and II) to digest raw starch. Ef-Amy I retains its catalytic activity at temperatures below 10°C. To identify the molecular properties of Ef-Amy I, X-ray crystal structures were determined of the wild type and of the inactive E249Q mutant. Ef-Amy I has structural similarities to mammalian α-amylases, including the porcine pancreatic and human pancreatic α-amylases. Structural comparisons of the overall structures as well as of the Ca2+-binding sites of Ef-Amy I and the mammalian α-amylases indicate that Ef-Amy I has increased structural flexibility and more solvent-exposed acidic residues. These structural features of Ef-Amy I may contribute to its observed catalytic activity at low temperatures, as many cold-adapted enzymes have similar structural properties. The structure of the substrate complex of the inactive mutant of Ef-Amy I shows that a maltohexaose molecule is bound in the active site and a maltotetraose molecule is bound in the cleft between the N- and C-terminal domains. The recognition of substrate molecules by Ef-Amy I exhibits some differences from that observed in structures of human pancreatic α-amylase. This result provides insights into the structural modulation of the recognition of substrates and inhibitors.

Insights into the role of electrostatics in temperature adaptation: a comparative study of psychrophilic, mesophilic, and thermophilic subtilisin-like serine proteases

To investigate the role of electrostatics in different temperature adaptations, we performed a comparative study on subtilisin-like serine proteases from psychrophilic Vibrio sp. PA-44 (VPR), mesophilic Engyodontium album (Tritirachium album) (PRK), and thermophilic Thermus aquaticus (AQN) using multiple-replica molecular dynamics (MD) simulations combined with continuum electrostatics calculations. The results reveal that although salt bridges are not a crucial factor in determining the overall thermostability of these three proteases, they on average provide the greatest, moderate, and least electrostatic stabilization to AQN, PRK, and VPR, respectively, at the respective organism growth temperatures. Most salt bridges in AQN are effectively stabilizing and thus contribute to maintaining the overall structural stability at 343 K, while nearly half of the salt bridges in VPR interconvert between being stabilizing and being destabilizing, likely aiding in enhancing the local conformational flexibility at 283 K. The individual salt bridges, salt-bridge networks, and calcium ions contribute differentially to local stability and flexibility of these three enzyme structures, depending on their spatial distributions and electrostatic strengths. The shared negatively charged surface potential at the active center of the three enzymes may provide the active-center flexibility necessary for nucleophilic attack and proton transfer. The differences in distributions of the electro-negative, electro-positive, and electro-neutral potentials, particularly over the back surfaces of the three proteases, may modulate/affect not only protein solubility and thermostability but also structural stability and flexibility/rigidity. These results demonstrate that electrostatics contributes to both heat and cold adaptation of subtilisin-like serine proteases through fine-tuning, either globally or locally, the structural stability and conformational flexibility/rigidity, thus providing a foundation for further engineering and mutagenesis studies.

Differently charged surface patches contribute to temperature adaptation of subtilisin-like serine proteases through affecting/modulating the protein solubility and thermostability and the structural flexibility/rigidity/stability.

Characterization of a cold-adapted DNA photolyase from C. psychrerythraea 34H

The phrB gene encoding a putative cold-adapted DNA photolyase was cloned from the bacterial genomic DNA of Colwellia psychrerythraea 34H, a psychrophilic bacterium. Recombinant DNA photolyase, rCpPL, was overexpressed and purified from three different vectors. rCpPL binds its DNA substrate by flipping a cyclobutane pyrimidine dimer (CPD) into its active site and repairs CPD-containing DNA in vitro. rCpPL contains one catalytic flavin adenine dinucleotide (FAD) cofactor, but displays promiscuity in cofactor binding, in which either a flavin mononucleotide (FMN) or a methenyltetrahydrofolate (MTHF) molecule is bound as an antenna molecule and found in sub-stoichiometric amounts. The UV/Vis spectrum of oxidized rCpPL shows that the FAD_OX absorption maximum is the most red-shifted reported for a PL, suggesting a unique cavity electrostatic environment. Modest FAD vibronic structure suggests that the binding pocket is more flexible than warmer PLs, corroborating the hypothesis that psychrophilic proteins must be highly flexible to function at low temperatures. Fluorescence excitation data show that the freshly purified flavin cofactor is in its fully reduced state (FADH¯). A homology analysis of PL protein structures spanning 70 °C in growth temperature supports the data that the structure of CpPL is quite different from its warmer cousins.

Scrutiny of electrostatic-driven conformational ordering of polypeptide chains in DMSO: a study with a model oligopeptide

The molecular mechanism of DMSO-induced stabilisation of β-sheets is attributed to the combination of polar electrostatic interactions among side chains, and backbone desolvation through bulky side chains which promotes backbone hydrogen bonding.

Discovery, Molecular Mechanisms, and Industrial Applications of Cold-Active Enzymes

Cold-active enzymes constitute an attractive resource for biotechnological applications. Their high catalytic activity at temperatures below 25°C makes them excellent biocatalysts that eliminate the need of heating processes hampering the quality, sustainability, and cost-effectiveness of industrial production. Here we provide a review of the isolation and characterization of novel cold-active enzymes from microorganisms inhabiting different environments, including a revision of the latest techniques that have been used for accomplishing these paramount tasks. We address the progress made in the overexpression and purification of cold-adapted enzymes, the evolutionary and molecular basis of their high activity at low temperatures and the experimental and computational techniques used for their identification, along with protein engineering endeavors based on these observations to improve some of the properties of cold-adapted enzymes to better suit specific applications. We finally focus on examples of the evaluation of their potential use as biocatalysts under conditions that reproduce the challenges imposed by the use of solvents and additives in industrial processes and of the successful use of cold-adapted enzymes in biotechnological and industrial applications.

Crystal Structures of Two Isozymes of Citrate Synthase from Sulfolobus tokodaii Strain 7

Thermoacidophilic archaeon Sulfolobus tokodaii strain 7 has two citrate synthase genes (ST1805-CS and ST0587-CS) in the genome with 45% sequence identity. Because they exhibit similar optimal temperatures of catalytic activity and thermal inactivation profiles, we performed structural comparisons between these isozymes to elucidate adaptation mechanisms to high temperatures in thermophilic CSs. The crystal structures of ST1805-CS and ST0587-CS were determined at 2.0 Å and 2.7 Å resolutions, respectively. Structural comparison reveals that both of them are dimeric enzymes composed of two identical subunits, and these dimeric structures are quite similar to those of citrate synthases from archaea and eubacteria. ST0587-CS has, however, 55 ion pairs within whole dimer structure, while having only 36 in ST1805-CS. Although the number and distributions of ion pairs are distinct from each other, intersubunit ion pairs between two domains of each isozyme are identical especially in interterminal region. Because the location and number of ion pairs are in a trend with other CSs from thermophilic microorganisms, the factors responsible for thermal adaptation of ST-CS isozymes are characterized by ion pairs in interterminal region.

Structure Prediction of a Novel Exo-β-1,3-Glucanase: Insights into the Cold Adaptation of Psychrophilic Yeast Glaciozyma antarctica PI12

We report a detailed structural analysis of the psychrophilic exo-β-1,3-glucanase (GaExg55) from Glaciozyma antarctica PI12. This study elucidates the structural basis of exo-1,3-β-1,3-glucanase from this psychrophilic yeast. The structural prediction of GaExg55 remains a challenge because of its low sequence identity (37 %). A 3D model was constructed for GaExg55. Threading approach was employed to determine a suitable template and generate optimal target-template alignment for establishing the model using MODELLER9v15. The primary sequence analysis of GaExg55 with other mesophilic exo-1,3-β-glucanases indicated that an increased flexibility conferred to the enzyme by a set of amino acids substitutions in the surface and loop regions of GaExg55, thereby facilitating its structure to cold adaptation. A comparison of GaExg55 with other mesophilic exo-β-1,3-glucanases proposed that the catalytic activity and structural flexibility at cold environment were attained through a reduced amount of hydrogen bonds and salt bridges, as well as an increased exposure of the hydrophobic side chains to the solvent. A molecular dynamics simulation was also performed using GROMACS software to evaluate the stability of the GaExg55 structure at varying low temperatures. The simulation result confirmed the above findings for cold adaptation of the psychrophilic GaExg55. Furthermore, the structural analysis of GaExg55 with large catalytic cleft and wide active site pocket confirmed the high activity of GaExg55 to hydrolyze polysaccharide substrates.

Physicochemical nature of interfaces controlling ferredoxin NADP(+) reductase activity through its interprotein interactions with ferredoxin

Although acidic residues of ferredoxin (Fd) are known to be essential for activities of various Fd-dependent enzymes, including ferredoxin NADP(+) reductase (FNR) and sulfite reductase (SiR), through electrostatic interactions with basic residues of partner enzymes, non-electrostatic contributions such as hydrophobic forces remain largely unknown. We herein demonstrated that intermolecular hydrophobic and charge-charge interactions between Fd and enzymes were both critical for enzymatic activity. Systematic site-directed mutagenesis, which altered physicochemical properties of residues on the interfaces of Fd for FNR /SiR, revealed various changes in activities of both enzymes. The replacement of serine 43 of Fd to a hydrophobic residue (S43W) and charged residue (S43D) increased and decreased FNR activity, respectively, while S43W showed significantly lower SiR activity without affecting SiR activity by S43D, suggesting that hydrophobic and electrostatic interprotein forces affected FNR activity. Enzyme kinetics revealed that changes in FNR activity by mutating Fd correlated with Km, but not with kcat or activation energy, indicating that interprotein interactions determined FNR activity. Calorimetry-based binding thermodynamics between Fd and FNR showed different binding modes of FNR to wild-type, S43W, or S43D, which were controlled by enthalpy and entropy, as shown by the driving force plot. Residue-based NMR spectroscopy of (15)N FNR with Fds also revealed distinct binding modes of each complex based on different directions of NMR peak shifts with similar overall chemical shift differences. We proposed that subtle adjustments in both hydrophobic and electrostatic forces were critical for enzymatic activity, and these results may be applicable to protein-based electron transfer systems.

Structural prediction of a novel laminarinase from the psychrophilic Glaciozyma antarctica PI12 and its temperature adaptation analysis

Here, we present a novel psychrophilic β-glucanase from Glaciozyma antarctica PI12 yeast that has been structurally modeled and analyzed in detail. To our knowledge, this is the first attempt to model a psychrophilic laminarinase from yeast. Because of the low sequence identity (<40%), a threading method was applied to predict a 3D structure of the enzyme using the MODELLER9v12 program. The results of a comparative study using other mesophilic, thermophilic, and hyperthermophilic laminarinases indicated several amino acid substitutions on the surface of psychrophilic laminarinase that totally increased the flexibility of its structure for efficient catalytic reactions at low temperatures. Whereas several structural factors in the overall structure can explain the weak thermal stability, this research suggests that the psychrophilic adaptation and catalytic activity at low temperatures were achieved through existence of longer loops and shorter or broken helices and strands, an increase in the number of aromatic and hydrophobic residues, a reduction in the number of hydrogen bonds and salt bridges, a higher total solvent accessible surface area, and an increase in the exposure of the hydrophobic side chains to the solvent. The results of comparative molecular dynamics simulation and principal component analysis confirmed the above strategies adopted by psychrophilic laminarinase to increase its catalytic efficiency and structural flexibility to be active at cold temperature.

Cold denaturation of α-synuclein amyloid fibrils

Although amyloid fibrils are associated with numerous pathologies, their conformational stability remains largely unclear. Herein, we probe the thermal stability of various amyloid fibrils. α-Synuclein fibrils cold-denatured to monomers at 0-20 °C and heat-denatured at 60-110 °C. Meanwhile, the fibrils of β2-microglobulin, Alzheimer's Aβ1-40/Aβ1-42 peptides, and insulin exhibited only heat denaturation, although they showed a decrease in stability at low temperature. A comparison of structural parameters with positive enthalpy and heat capacity changes which showed opposite signs to protein folding suggested that the burial of charged residues in fibril cores contributed to the cold denaturation of α-synuclein fibrils. We propose that although cold-denaturation is common to both native proteins and misfolded fibrillar states, the main-chain dominated amyloid structures may explain amyloid-specific cold denaturation arising from the unfavorable burial of charged side-chains in fibril cores.

Structural and functional analysis of a novel psychrophilic β-mannanase from Glaciozyma antarctica PI12

The structure of a novel psychrophilic β-mannanase enzyme from Glaciozyma antarctica PI12 yeast has been modelled and analysed in detail. To our knowledge, this is the first attempt to model a psychrophilic β-mannanase from yeast. To this end, a 3D structure of the enzyme was first predicted using a threading method because of the low sequence identity (<30%) using MODELLER9v12 and simulated using GROMACS at varying low temperatures for structure refinement. Comparisons with mesophilic and thermophilic mannanases revealed that the psychrophilic mannanase contains longer loops and shorter helices, increases in the number of aromatic and hydrophobic residues, reductions in the number of hydrogen bonds and salt bridges and numerous amino acid substitutions on the surface that increased the flexibility and its efficiency for catalytic reactions at low temperatures.

Circular dichroism and fluorescence spectroscopy of cysteinyl-tRNA synthetase from Halobacterium salinarum ssp. NRC-1 demonstrates that group I cations are particularly effective in providing structure and stability to this halophilic protein

Proteins from extremophiles have the ability to fold and remain stable in their extreme environment. Here, we investigate the presence of this effect in the cysteinyl-tRNA synthetase from Halobacterium salinarum ssp. NRC-1 (NRC-1), which was used as a model halophilic protein. The effects of salt on the structure and stability of NRC-1 and of E. coli CysRS were investigated through far-UV circular dichroism (CD) spectroscopy, fluorescence spectroscopy, and thermal denaturation melts. The CD of NRC-1 CysRS was examined in different group I and group II chloride salts to examine the effects of the metal ions. Potassium was observed to have the strongest effect on NRC-1 CysRS structure, with the other group I salts having reduced strength. The group II salts had little effect on the protein. This suggests that the halophilic adaptations in this protein are mediated by potassium. CD and fluorescence spectra showed structural changes taking place in NRC-1 CysRS over the concentration range of 0-3 M KCl, while the structure of E. coli CysRS was relatively unaffected. Salt was also shown to increase the thermal stability of NRC-1 CysRS since the melt temperature of the CysRS from NRC-1 was increased in the presence of high salt, whereas the E. coli enzyme showed a decrease. By characterizing these interactions, this study not only explains the stability of halophilic proteins in extremes of salt, but also helps us to understand why and how group I salts stabilize proteins in general.

Protein adaptations in archaeal extremophiles

Extremophiles, especially those in Archaea, have a myriad of adaptations that keep their cellular proteins stable and active under the extreme conditions in which they live. Rather than having one basic set of adaptations that works for all environments, Archaea have evolved separate protein features that are customized for each environment. We categorized the Archaea into three general groups to describe what is known about their protein adaptations: thermophilic, psychrophilic, and halophilic. Thermophilic proteins tend to have a prominent hydrophobic core and increased electrostatic interactions to maintain activity at high temperatures. Psychrophilic proteins have a reduced hydrophobic core and a less charged protein surface to maintain flexibility and activity under cold temperatures. Halophilic proteins are characterized by increased negative surface charge due to increased acidic amino acid content and peptide insertions, which compensates for the extreme ionic conditions. While acidophiles, alkaliphiles, and piezophiles are their own class of Archaea, their protein adaptations toward pH and pressure are less discernible. By understanding the protein adaptations used by archaeal extremophiles, we hope to be able to engineer and utilize proteins for industrial, environmental, and biotechnological applications where function in extreme conditions is required for activity.

Function and biotechnology of extremophilic enzymes in low water activity

Enzymes from extremophilic microorganisms usually catalyze chemical reactions in non-standard conditions. Such conditions promote aggregation, precipitation, and denaturation, reducing the activity of most non-extremophilic enzymes, frequently due to the absence of sufficient hydration. Some extremophilic enzymes maintain a tight hydration shell and remain active in solution even when liquid water is limiting, e.g. in the presence of high ionic concentrations, or at cold temperature when water is close to the freezing point. Extremophilic enzymes are able to compete for hydration via alterations especially to their surface through greater surface charges and increased molecular motion. These properties have enabled some extremophilic enzymes to function in the presence of non-aqueous organic solvents, with potential for design of useful catalysts. In this review, we summarize the current state of knowledge of extremophilic enzymes functioning in high salinity and cold temperatures, focusing on their strategy for function at low water activity. We discuss how the understanding of extremophilic enzyme function is leading to the design of a new generation of enzyme catalysts and their applications to biotechnology.

Molecular dynamics of mesophilic-like mutants of a cold-adapted enzyme: insights into distal effects induced by the mutations

Networks and clusters of intramolecular interactions, as well as their "communication" across the three-dimensional architecture have a prominent role in determining protein stability and function. Special attention has been dedicated to their role in thermal adaptation. In the present contribution, seven previously experimentally characterized mutants of a cold-adapted α-amylase, featuring mesophilic-like behavior, have been investigated by multiple molecular dynamics simulations, essential dynamics and analyses of correlated motions and electrostatic interactions. Our data elucidate the molecular mechanisms underlying the ability of single and multiple mutations to globally modulate dynamic properties of the cold-adapted α-amylase, including both local and complex unpredictable distal effects. Our investigation also shows, in agreement with the experimental data, that the conversion of the cold-adapted enzyme in a warm-adapted variant cannot be completely achieved by the introduction of few mutations, also providing the rationale behind these effects. Moreover, pivotal residues, which are likely to mediate the effects induced by the mutations, have been identified from our analyses, as well as a group of suitable candidates for protein engineering. In fact, a subset of residues here identified (as an isoleucine, or networks of mesophilic-like salt bridges in the proximity of the catalytic site) should be considered, in experimental studies, to get a more efficient modification of the features of the cold-adapted enzyme.

Contribution of inter-subunit interactions to the thermostability of Pyrococcus furiosus citrate synthase

Using citrate synthase from the hyperthermophile Pyrococcus furiosus (PfCS) as our test molecule, we show through guanidine hydrochloride-induced unfolding that the dimer separates into folded, but inactive, monomers before individual subunit unfolding takes place. Given that forces across the dimer interface are vital for thermostability, a robust computational method was derived that uses the University of Houston Brownian Dynamics (UHBD) program to calculate both the hydrophobic and electrostatic contribution to the dimerisation energy at 100°C. The results from computational and experimental determination of the lowered stability of interface mutants were correlated, being both of the same order of magnitude and placing the mutant proteins in the same order of stability. This computational method, optimised for hyperthermophilic molecules and tested in the laboratory, after further testing on other examples, could be of widespread use in the prediction of thermostabilising mutations in other oligomeric proteins for which dissociation is the first step in unfolding.

A novel method for high-level production of psychrophilic TAB5 alkaline phosphatase

Heat labile alkaline phosphatases (APs) are widely used in biomedical research for they can easily be heat inactivated once they have done their job. Here we reported a novel method for high-level production of recombinant psychrophilic Antarctic bacterium strain TAB5 alkaline phosphatase (TAP) in Escherichia coli. We synthesized the whole TAP gene according to the synonymous codon choice that is optimal for the E. coli translational system. Then the gene was cloned into pT7 expression vector, expressed in BL21 (DE3) pLysS strain by auto-induction system. The recombinant protein was purified by Ni-NTA affinity chromatography and anion exchange chromatography. The typical yield was 90.9 mg protein from 16.2 g wet cells in 1L culture medium with the purity over 99%, 340 times as many mg/L (and 21 times the mg/g cells) compared to previous methods. The dephosphorylation activity assay showed that the purified recombinant TAP has similar activity to calf intestinal alkaline phosphatase in room temperature, and it can be totally inactivated by treatment at 60°C for 15 min. Our research provides a novel method for high-level expression, purification and characterization of TAP which is sufficient for high throughput genome analysis and may replace the widely used shrimp AP because of its low cost.

Environmental adaptation of proteins: regression models with simple physicochemical properties

Bio-sequences from ortholog proteins are well suited for statistical inference when the sequences can be divided into ordinal groups based on known environmental features or traits of the host organisms. In this paper two new regression models are described for extracting proteomic trends of extreme environments. The approach is based on physicochemical properties of the amino acids, and may also utilise stratification of the data. We are especially looking for connections of temperature adaptation between the organism and its molecular level. To show the applicability of the methods, we present analyses of genomic data from proteobacteria orders, where we examine the cold adaptation of membrane proteins, intracellular proteins, and the enzyme endonuclease I. Our results confirm earlier findings that redistribution of charge and increase of surface hydrophobicity might be some of the most important signatures for cold adaptation.

Characteristics of mutants designed to incorporate a new ion pair into the structure of a cold adapted subtilisin-like serine proteinase

Structural comparisons of VPR, a subtilisin-like serine proteinase from a psychrotrophic Vibrio species and a thermophilic homologue, aqualysin I, have led us to hypothesize about the roles of different residues in the temperature adaptation of the enzymes. Some of these hypotheses are now being examined by analysis of mutants of the enzymes. The selected substitutions are believed to increase the stability of the cold adapted enzyme based on structural analysis of the thermostable structure. We report here on mutants, which were designed to incorporate an ion pair into the structure of VPR. The residues Asp17 and Arg259 are assumed to form an ion pair in aqualysin I. The cold adapted VPR contains Asn (Asn15) and Lys (Lys257) at corresponding sites in its structure. In VPR, Asn 15 is located on a surface loop with its side group pointing towards the side chain of Lys257. By substituting Asn15 by Asp (N15D) it was considered feasible that a salt bridge would form between the oppositely charged groups. To mimic further the putative salt bridge from the thermophile enzyme the corresponding double mutant (N15D/K257R) was also produced. The N15D mutation increased the thermal stability of VPR by approximately 3 degrees C, both in T(50%) and T(m). Addition of the K257R mutation did not however, increase the stability of the double mutant any further. Despite this stabilization of the VPR mutants the catalytic activity (k(cat)) against the substrate Suc-AAPF-NH-Np was increased in the mutants. Molecular dynamics simulations on wild type and the two mutant proteins suggested that indeed a salt bridge was formed in both cases. Furthermore, a truncated form of the N15D mutant (N15DDeltaC) was produced, lacking a 15 residue long C-terminal extended sequence not present in the thermophilic enzyme. In wild type VPR this supposedly moveable, negatively charged arm on the protein molecule might interfere with the new salt bridge introduced as a result of the N15D mutation. Removal of the C-terminal arm improved the thermal stability (T(m) approximately +1.5 degrees C) of the truncated enzyme (VPRDeltaC) as compared to the wild type VPR. Introduction of the N15D substitution into VPRDeltaC improved the thermal stability further by about 3 degrees C, or to about the same extent as in the wild type. However, contrary to what was observed for the wild type, the introduction of the putative salt bridge did not affect the catalytic properties (k(cat)) of the C-terminal truncated enzyme.

Cold adaptation of enzyme reaction rates

A major issue for organisms living at extreme temperatures is to preserve both stability and activity of their enzymes. Cold-adapted enzymes generally have a reduced thermal stability, to counteract freezing, and show a lower enthalpy and a more negative entropy of activation compared to mesophilic and thermophilic homologues. Such a balance of thermodynamic activation parameters can make the reaction rate decrease more linearly, rather than exponentially, as the temperature is lowered, but the structural basis for rate optimization toward low working temperatures remains unclear. In order to computationally address this problem, it is clear that reaction simulations rather than standard molecular dynamics calculations are needed. We have thus carried out extensive computer simulations of the keto-enol(ate) isomerization steps in differently adapted citrate synthases to explore the structure-function relationships behind catalytic rate adaptation to different temperatures. The calculations reproduce the absolute rates of the psychrophilic and mesophilic enzymes at 300 K, as well as the lower enthalpy and more negative entropy of activation of the cold-adapted enzyme, where the latter simulation result is obtained from high-precision Arrhenius plots. The overall catalytic effect originates from electrostatic stabilization of the transition state and enolate and the reduction of reorganization free energy. The simulations, however, show psychrophilic, mesophilic, and hyperthermophilic citrate synthases to have increasingly stronger electrostatic stabilization of the transition state, while the energetic penalty in terms of internal protein interactions follows the reverse order with the cold-adapted enzyme having the most favorable energy term. The lower activation enthalpy and more negative activation entropy observed for cold-adapted enzymes are found to be associated with a decreased protein stiffness. The origin of this effect is, however, not localized to the active site but to other regions of the protein structure.

Large-scale modulation of thermodynamic protein folding barriers linked to electrostatics

Protein folding barriers, which range from zero to the tens of RT that result in classical two-state kinetics, are primarily determined by protein size and structural topology [Plaxco KW, Simons KT, Baker D (1998) J Mol Biol 277:985-994]. Here, we investigate the thermodynamic folding barriers of two relatively large proteins of the same size and topology: bovine alpha-lactalbumin (BLA) and hen-egg-white lysozyme (HEWL). From the analysis of differential scanning calorimetry experiments with the variable-barrier model [Muñoz V, Sanchez-Ruiz JM (2004) Proc Natl Acad Sci USA 101:17646-17651] we obtain a high barrier for HEWL and a marginal folding barrier for BLA. These results demonstrate a remarkable tuning range of at least 30 kJ/mol (i.e., five to six orders of magnitude in population) within a unique protein scaffold. Experimental and theoretical analyses on these proteins indicate that the surprisingly small thermodynamic folding barrier of BLA arises from the stabilization of partially unfolded conformations by electrostatic interactions. Interestingly, there is clear reciprocity between the barrier height and the biological function of the two proteins, suggesting that the marginal barrier of BLA is a product of natural selection. Electrostatic surface interactions thus emerge as a mechanism for the modulation of folding barriers in response to special functional requirements within a given structural fold.

Cooperativity of complex salt bridges

The energetic contribution of complex salt bridges, in which one charged residue (anchor residue) forms salt bridges with two or more residues simultaneously, has been suggested to have importance for protein stability. Detailed analysis of the net energetics of complex salt bridge formation using double- and triple-mutant cycle analysis revealed conflicting results. In two cases, it was shown that complex salt bridge formation is cooperative, i.e., the net strength of the complex salt bridge is more than the sum of the energies of individual pairs. In one case, it was reported that complex salt bridge formation is anti-cooperative. To resolve these different findings, we performed analysis of the geometries of salt bridges in a representative set of structures from the PDB and found that over 87% of all complex salt bridges anchored by Arg/Lys have a geometry such that the angle formed by their Calpha atoms, Theta, is <90 degrees . This preferred geometry is observed in the two reported instances when the energetics of complex salt bridge formation is cooperative, while in the reported anti-cooperative complex salt bridge, Theta is close to 160 degrees . Based on these observations, we hypothesized that complex salt bridges are cooperative for Theta < 90 degrees and anti-cooperative for 90 degrees < Theta < 180 degrees . To provide a further experimental test for this hypothesis, we engineered a complex salt bridge with Theta = 150 degrees into a model protein, the activation domain of human procarboxypeptidase A2 (ADA2h). Experimentally derived stabilities of the ADA2h variants allowed us to show that the complex salt bridge in ADA2h is anti-cooperative.

Electrostatic interactions play an essential role in DNA repair and cold-adaptation of uracil DNA glycosylase

Life has adapted to most environments on earth, including low and high temperature niches. The increased catalytic efficiency and thermoliability observed for enzymes from organisms living in constantly cold regions when compared to their mesophilic and thermophilic cousins are poorly understood at the molecular level. Uracil DNA glycosylase (UNG) from cod (cUNG) catalyzes removal of uracil from DNA with an increased k(cat) and reduced K(m) relative to its warm-active human (hUNG) counterpart. Specific issues related to DNA repair and substrate binding/recognition (K(m)) are here investigated by continuum electrostatics calculations, MD simulations and free energy calculations. Continuum electrostatic calculations reveal that cUNG has surface potentials that are more complementary to the DNA potential at and around the catalytic site when compared to hUNG, indicating improved substrate binding. Comparative MD simulations combined with free energy calculations using the molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) method show that large opposing energies are involved when forming the enzyme-substrate complexes. Furthermore, the binding free energies obtained reveal that the Michaelis-Menten complex is more stable for cUNG, primarily due to enhanced electrostatic properties, suggesting that energetic fine-tuning of electrostatics can be utilized for enzymatic temperature adaptation. Energy decomposition pinpoints the residual determinants responsible for this adaptation.

Dispersion interactions govern the strong thermal stability of a protein

Rubredoxin from the hyperthermophile Pyrococcus furiosus (Pf Rd) is an extremely thermostable protein, which makes it an attractive subject of protein folding and stability studies. A fundamental question arises as to what the reason for such extreme stability is and how it can be elucidated from a complex set of interatomic interactions. We addressed this issue first theoretically through a computational analysis of the hydrophobic core of the protein and its mutants, including the interactions taking place inside the core. Here we show that a single mutation of one of phenylalanine's residues inside the protein's hydrophobic core results in a dramatic decrease in its thermal stability. The calculated unfolding Gibbs energy as well as the stabilization energy differences between a few core residues follows the same trend as the melting temperature of protein variants determined experimentally by microcalorimetry measurements. NMR spectroscopy experiments have shown that the only part of the protein affected by mutation is the reasonably rearranged hydrophobic core. It is hence concluded that stabilization energies, which are dominated by London dispersion, represent the main source of stability of this protein.

Structural and Functional Properties of Isocitrate Dehydrogenase from the Psychrophilic Bacterium Desulfotalea psychrophila Reveal a Cold-active Enzyme with an Unusual High Thermal Stability

Isocitrate dehydrogenase (IDH) has been studied extensively due to its central role in the Krebs cycle, catalyzing the oxidative NAD(P)(+)-dependent decarboxylation of isocitrate to alpha-ketoglutarate and CO(2). Here, we present the first crystal structure of IDH from a psychrophilic bacterium, Desulfotalea psychrophila (DpIDH). The structural information is combined with a detailed biochemical characterization and a comparative study with IDHs from the mesophilic bacterium Desulfitobacterium hafniense (DhIDH), porcine (PcIDH), human cytosolic (HcIDH) and the hyperthermophilic Thermotoga maritima (TmIDH). DpIDH was found to have a higher melting temperature (T(m)=66.9 degrees C) than its mesophilic homologues and a suboptimal catalytic efficiency at low temperatures. The thermodynamic activation parameters indicated a disordered active site, as seen also for the drastic increase in K(m) for isocitrate at elevated temperatures. A methionine cluster situated at the dimeric interface between the two active sites and a cluster of destabilizing charged amino acids in a region close to the active site might explain the poor isocitrate affinity. On the other hand, DpIDH was optimized for interacting with NADP(+) and the crystal structure revealed unique interactions with the cofactor. The highly acidic surface, destabilizing charged residues, fewer ion pairs and reduced size of ionic networks in DpIDH suggest a flexible global structure. However, strategic placement of ionic interactions stabilizing the N and C termini, and additional ionic interactions in the clasp domain as well as two enlarged aromatic clusters might counteract the destabilizing interactions and promote the increased thermal stability. The structure analysis of DpIDH illustrates how psychrophilic enzymes can adjust their flexibility in dynamic regions during their catalytic cycle without compromising the global stability of the protein.

Structure of phenylalanine hydroxylase from Colwellia psychrerythraea 34H, a monomeric cold active enzyme with local flexibility around the active site and high overall stability

The characteristic of cold-adapted enzymes, high catalytic efficiency at low temperatures, is often associated with low thermostability and high flexibility. In this context, we analyzed the catalytic properties and solved the crystal structure of phenylalanine hydroxylase from the psychrophilic bacterium Colwellia psychrerythraea 34H (CpPAH). CpPAH displays highest activity with tetrahydrobiopterin (BH(4)) as cofactor and at 25 degrees C (15 degrees C above the optimal growth temperature). Although the enzyme is monomeric with a single L-Phe-binding site, the substrate binds cooperatively. In comparison with PAH from mesophilic bacteria and mammalian organisms, CpPAH shows elevated [S(0.5)](L-Phe) (= 1.1 +/- 0.1 mm) and K(m)(BH(4))(= 0.3 +/- 0.1 mm), as well as high catalytic efficiency at 10 degrees C. However, the half-inactivation and denaturation temperature is only slightly lowered (T(m) approximately 52 degrees C; where T(m) is half-denaturation temperature), in contrast to other cold-adapted enzymes. The crystal structure shows regions of local flexibility close to the highly solvent accessible binding sites for BH(4) (Gly(87)/Phe(88)/Gly(89)) and l-Phe (Tyr(114)-Pro(118)). Normal mode and COREX analysis also detect these and other areas with high flexibility. Greater mobility around the active site and disrupted hydrogen bonding abilities for the cofactor appear to represent cold-adaptive properties that do not markedly affect the thermostability of CpPAH.

Characterizing structural features of cuticle-degrading proteases from fungi by molecular modeling

Background

Serine proteases secreted by nematode and insect pathogenic fungi are bio-control agents which have commercial potential for developing into effective bio-pesticides. A thorough understanding of the structural and functional features of these proteases would significantly assist with targeting the design of efficient bio-control agents.

Results

Structural models of serine proteases PR1 from entomophagous fungus, Ver112 and VCP1 from nematophagous fungi, have been modeled using the homology modeling technique based on the crystal coordinate of the proteinase K. In combination with multiple sequence alignment, these models suggest one similar calcium-binding site and two common disulfide bridges in the three cuticle-degrading enzymes. In addition, the predicted models of the three cuticle-degrading enzymes present an essentially identical backbone topology and similar geometric properties with the exception of a limited number of sites exhibiting relatively large local conformational differences only in some surface loops and the N-, C termini. However, they differ from each other in the electrostatic surface potential, in hydrophobicity and size of the S4 substrate-binding pocket, and in the number and distribution of hydrogen bonds and salt bridges within regions that are part of or in close proximity to the S2-loop.

Conclusion

These differences likely lead to variations in substrate specificity and catalytic efficiency among the three enzymes. Amino acid polymorphisms in cuticle-degrading enzymes were discussed with respect to functional effects and host preference. It is hoped that these structural models would provide a further basis for exploitation of these serine proteases from pathogenic fungi as effective bio-control agents.

Discovery of a thermophilic protein complex stabilized by topologically interlinked chains

A growing number of organisms have been discovered inhabiting extreme environments, including temperatures in excess of 100 degrees C. How cellular proteins from such organisms retain their native folds under extreme conditions is still not fully understood. Recent computational and structural studies have identified disulfide bonding as an important mechanism for stabilizing intracellular proteins in certain thermophilic microbes. Here, we present the first proteomic analysis of intracellular disulfide bonding in the hyperthermophilic archaeon Pyrobaculum aerophilum. Our study reveals that the utilization of disulfide bonds extends beyond individual proteins to include many protein-protein complexes. We report the 1.6 A crystal structure of one such complex, a citrate synthase homodimer. The structure contains two intramolecular disulfide bonds, one per subunit, which result in the cyclization of each protein chain in such a way that the two chains are topologically interlinked, rendering them inseparable. This unusual feature emphasizes the variety and sophistication of the molecular mechanisms that can be achieved by evolution.