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.     

Structural and functional adaptations to extreme temperatures in psychrophilic,mesophilic, and thermophilic DNA ligases

Psychrophiles, host of permanently cold habitats, display metabolic fluxes comparable to those exhibited by mesophilic organisms at moderate temperatures. These organisms have evolved by producing, among other peculiarities, cold-active enzymes that have the properties to cope with the reduction of chemical reaction rates induced by low temperatures. The emerging picture suggests that these enzymes display a high catalytic efficiency at low temperatures through an improved flexibility of the structural components involved in the catalytic cycle, whereas other protein regions, if not implicated in catalysis, may be even more rigid than their mesophilic counterparts. In return, the increased flexibility leads to a decreased stability of psychrophilic enzymes. In order to gain further advances in the analysis of the activity/flexibility/stability concept, psychrophilic, mesophilic, and thermophilic DNA ligases have been compared by three-dimensional-modeling studies, as well as regards their activity, surface hydrophobicity, structural permeability, conformational stabilities, and irreversible thermal unfolding. These data show that the cold-adapted DNA ligase is characterized by an increased activity at low and moderate temperatures, an overall destabilization of the molecular edifice, especially at the active site, and a high conformational flexibility. The opposite trend is observed in the mesophilic and thermophilic counterparts, the latter being characterized by a reduced low temperature activity, high stability and reduced flexibility. These results strongly suggest a complex relationship between activity, flexibility and stability. In addition, they also indicate that in cold-adapted enzymes, the driving force for denaturation is a large entropy change.     

Catalytic efficiency and some structural properties of cold-active protein-tyrosine-phosphatase

A procedure was established for expression and purification of abundant recombinant cold-active protein-tyrosine-phosphatase (RCPTPase), which showed identical enzymatic characteristics to the native enzyme (NCPTPase). The purified RCPTPase showed high catalytic activity at low temperature and maximal activity at 30 degrees C. RCPTPase has a thermodynamic characteristic in that its activation enthalpy was determined to be low, 4.3 kcal/mol, at temperatures below 19.3 degrees C, where the Arrhenius relationship exhibited an inflection point, in comparison with 20.3 kcal/mol above 19.3 degrees C. Also, the thermostability, DeltaG(water), of the catalytic site in the RCPTPase molecule was increased with a decrease in temperature. It was considered that cold-active protein-tyrosine-phosphatase could maintain its catalytic site in a stable conformation for eliciting high catalytic activity with low activation enthalpy at low temperature.     

Kinetic and thermodynamic studies of peptidyltransferase in ribosomes from the extreme thermophile Thermus thermophilus

Throughout evolution, emerging organisms survived by adapting existing biochemical processes to new reaction conditions. Simple protein enzymes balanced changes in structural stability with changes that permitted optimal catalysis by adjustments in both entropic and enthalpic contributions to the free energy of activation for the reaction. Study of adaptive mechanisms by large multicomponent enzymes such as the ribosome has been largely unexplored. Here we have determined the kinetic and thermodynamic parameters of peptidyltransferase in ribosomes from the extreme thermophile Thermus thermophilus. Activity of thermophilic enzymes can be assayed over a wide range of temperatures, enabling one to measure accurate catalytic rates and determine enthalpic and entropic contributions to the free energy of activation of the reaction. Differences in the reaction conditions used here and in published studies on mesophilic ribosomes prevent direct comparison, but our data on Thermus ribosomes suggest that these ribosomes have adapted to changing environments using the same strategies as simple protein enzymes, balancing stability and flexibility without loss of catalytic rate. This strategy must be a very ancient process, perhaps first used by primitive ribosomes in the RNA World.     

Molecular basis of cold adaptation

Cold-adapted, or psychrophilic, organisms are able to thrive at low temperatures in permanently cold environments, which in fact characterize the greatest proportion of our planet. Psychrophiles include both prokaryotic and eukaryotic organisms and thus represent a significant proportion of the living world. These organisms produce cold-evolved enzymes that are partially able to cope with the reduction in chemical reaction rates induced by low temperatures. As a rule, cold-active enzymes display a high catalytic efficiency, associated however, with a low thermal stability. In most cases, the adaptation to cold is achieved through a reduction in the activation energy that possibly originates from an increased flexibility of either a selected area or of the overall protein structure. This enhanced plasticity seems in turn to be induced by the weak thermal stability of psychrophilic enzymes. The adaptation strategies are beginning to be understood thanks to recent advances in the elucidation of the molecular characteristics of cold-adapted enzymes derived from X-ray crystallography, protein engineering and biophysical methods. Psychrophilic organisms and their enzymes have, in recent years, increasingly attracted the attention of the scientific community due to their peculiar properties that render them particularly useful in investigating the possible relationship existing between stability, flexibility and specific activity and as valuable tools for biotechnological purposes.     

What are the roles of substrate-assisted catalysis and proximity effects in peptide bond formation by the ribosome?

The action of the peptidyl transferase center of the large ribosomal unit presents a fundamental step in the evolution from the RNA world to the protein world. Thus, it is important to understand the origin of the catalytic power of this ancient enzyme. Earlier studies suggested that the ribosome catalyzes peptide bond formation by using one of its groups as a general base, while more recent works have proposed that the catalysis is due to proximity effects or to substrate-assisted catalysis. However, the actual nature of the catalytic mechanism remains controversial. This work addresses the origin of the catalytic power of the ribosome by using computer simulation approaches and comparing the energetics of the peptide bond formation in the ribosome and in water. It is found that a significant part of the observed activation entropy of the reference solution reaction is due to solvation entropy, and that the proximity effect is smaller than previously thought. It is also found that the 2'-OH of the A76 ribose, which is associated with a large rate acceleration in the ribosome reaction, does not catalyze peptide bond formation in water. Thus, the catalytic effect cannot be attributed to substrate-assisted catalysis but rather to the effect of the ribosome on the reacting system. Overall, our calculations indicate that the reduction of the activation free energy is mainly due to electrostatic effects. The nature of these effects and their relationship to catalytic factors in modern enzymes is analyzed and discussed.     

Enzymatic activity assay of D-hydantoinase by isothermal titration calorimetry. Determination of the thermodynamic activation parameters for the hydrolysis of several substrates

Isothermal titration calorimetry (ITC) has been applied to the determination of the activity of D-hydantoinase (EC 3.5.2.2) with several substrates by monitoring the heat released during the reaction. The method is based on the proportionality between the reaction rate and the thermal power (heat/time) generated. Microcalorimetric assays carried out at different temperatures provided the dependence of the catalytic rate constant on temperature. We show that ITC assay is a nondestructive method that allows the determination of the catalytic rate constant (kcat), Michaelis constant (KM), activation energy and activation Gibbs energy, enthalpy and entropy of this reaction.     

The role of group bulkiness in the catalytic activity of psychrophile cold-active protein tyrosine phosphatase

The cold-active protein tyrosine phosphatase found in psychrophilic Shewanella species exhibits high catalytic efficiency at low temperatures as well as low thermostability, both of which are characteristics shared by many cold-active enzymes. The structure of cold-active protein tyrosine phosphatase is notable for the presence of three hydrophobic sites (termed the CA, Zn-1 and Zn-2 sites) behind the loop structures comprising the catalytic region. To identify the structural components responsible for specific enzyme characteristics, we determined the structure of wild-type cold-active protein tyrosine phosphatase at high resolution (1.1 A) and measured the catalytic efficiencies of enzymes containing mutations in the three hydrophobic sites. The bulkiness of the amino acid side chains in the core region of the Zn-1 site strongly affects the thermostability and the catalytic efficiency at low temperatures. The mutant enzyme I115M possessed a higher kcat at low temperatures. Elucidation of the crystal structure of I115M at a resolution of 1.5 A revealed that the loop structures involved in retaining the nucleophilic group and the acid catalyst are more flexible than in the wild-type enzyme.     

Catalysis and linear free energy relationships in aspartic proteases

Aspartic proteases are receiving considerable attention as potential drug targets in several serious diseases, such as AIDS, malaria, and Alzheimer's disease. These enzymes cleave polypeptide chains, often between specific amino acid residues, but despite the common reaction mechanism, they exhibit large structural differences. Here, the catalytic mechanism of aspartic proteases plasmepsin II, cathepsin D, and HIV-1 protease is examined by computer simulations utilizing the empirical valence bond approach in combination with molecular dynamics and free energy perturbation calculations. Free energy profiles are established for four different substrates, each six amino acids long and containing hydrophobic side chains in the P1 and P1' positions. Our simulations reproduce the catalytic effect of these enzymes, which accelerate the reaction rate by a factor of approximately 10(10) compared to that of the corresponding uncatalyzed reaction in water. The calculations elucidate the origin of the catalytic effect and allow a rationalization of the fact that, despite large structural differences between plasmepsin II/cathepsin D and HIV-1 protease, the magnitude of their rate enhancement is very similar. Amino acid residues surrounding the active site together with structurally conserved water molecules are found to play an important role in catalysis, mainly through dipolar (electrostatic) stabilization. A linear free energy relationship for the reactions in the different enzymes is established that also demonstrates the reduced reorganization energy in the enzymes compared to that in the uncatalyzed water reaction.     

The catalytic effect of dihydrofolate reductase and its mutants is determined by reorganization energies

The effect of distant mutations on the catalytic reaction of dihydrofolate reductase (DHFR) is reexamined by empirical valence bond simulations. The simulations reproduce for the first time the changes in the observed rate constants (without the use of adjustable parameters for this purpose) and show that the changes in activation barriers are strongly correlated with the corresponding changes in the reorganization energy. The preorganization of the polar groups of enzymes is the key catalytic factor, and anticatalytic mutations destroy this preorganization. Some anticatalytic mutations in DHFR also increase the distance between the donor and acceptor, but this effect is not directly related to catalysis since the native enzyme and the uncatalyzed reaction in water have similar average donor-acceptor distances. Insight into the effect of a mutation is provided by constructing the relevant free energy surfaces in terms of the generalized solute-solvent coordinates. It is shown how the mutations change the reaction coordinate and the activation barrier, and it is clarified that the corresponding changes do not reflect dynamical effects. It is also pointed out that all reactions in a condensed phase involve correlated motions (both in enzymes and in solution) and that the change of such motions upon mutations is a result of the change in the shape of the multidimensional reaction path on the solute-solvent surface, rather than the reason for the change in rate constant. Thus, as far as catalysis is concerned, the change in the activation barrier is due to the change in the electrostatic preorganization energy.     

Thermodynamic analysis of fast stages of specific lesion recognition by DNA repair enzymes

The methodology of determination of the thermodynamic parameters of fast stages of recognition and cleavage of DNA substrates is described for the enzymatic processes catalyzed by DNA glycosylases Fpg and hOGG1 and AP endonuclease APE1 during base excision repair (BER) pathway. For this purpose, stopped-flow pre-steady-state kinetic analysis of tryptophan fluorescence intensity changes in proteins and fluorophores in DNA substrates was performed at various temperatures. This approach made it possible to determine the changes of standard Gibbs free energy, enthalpy, and entropy of sequential steps of DNA-substrate binding, as well as activation enthalpy and entropy for the transition complex formation of the catalytic stage. The unified features of mechanism for search and recognition of damaged DNA sites by various enzymes of the BER pathway were discovered.     

A rigidifying salt-bridge favors the activity of thermophilic enzyme at high temperatures at the expense of low-temperature activity

Background

Thermophilic enzymes are often less active than their mesophilic homologues at low temperatures. One hypothesis to explain this observation is that the extra stabilizing interactions increase the rigidity of thermophilic enzymes and hence reduce their activity. Here we employed a thermophilic acylphosphatase from Pyrococcus horikoshii and its homologous mesophilic acylphosphatase from human as a model to study how local rigidity of an active-site residue affects the enzymatic activity.

Methods and Findings

Acylphosphatases have a unique structural feature that its conserved active-site arginine residue forms a salt-bridge with the C-terminal carboxyl group only in thermophilic acylphosphatases, but not in mesophilic acylphosphatases. We perturbed the local rigidity of this active-site residue by removing the salt-bridge in the thermophilic acylphosphatase and by introducing the salt-bridge in the mesophilic homologue. The mutagenesis design was confirmed by x-ray crystallography. Removing the salt-bridge in the thermophilic enzyme lowered the activation energy that decreased the activation enthalpy and entropy. Conversely, the introduction of the salt-bridge to the mesophilic homologue increased the activation energy and resulted in increases in both activation enthalpy and entropy. Revealed by molecular dynamics simulations, the unrestrained arginine residue can populate more rotamer conformations, and the loss of this conformational freedom upon the formation of transition state justified the observed reduction in activation entropy.

Conclusions

Our results support the conclusion that restricting the active-site flexibility entropically favors the enzymatic activity at high temperatures. However, the accompanying enthalpy-entropy compensation leads to a stronger temperature-dependency of the enzymatic activity, which explains the less active nature of the thermophilic enzymes at low temperatures.     

Enzyme-catalyzed decomposition of dibenzoyl peroxide in organic solvents

Catalytic activity of catalase (CAT, EC 1.11.1.6), immobilized on carbon black NORIT and soot PM-100, with respect to decomposition of dibenzoyl peroxide (BPO) in non-aqueous media (acetonitrile and tetrachloromethane), was investigated with a quantitative UV-spectrophotometrical approach. Progress of the above reaction was controlled by selected kinetic parameters: the apparent Michaelis constant (Km(app)), the specific rate constant (k(sp)), the activation energy (Ea), the maximum reaction rate (Vmax), and the Arrhenius' pre-exponential factor (Z0). Conclusions on the tentative mechanism of the catalytic process observed were drawn from the calculated values of the Gibbs energy of activation (deltaG*), the enthalpy of activation (deltaH*), and entropy of activation (deltaS*).     

Cold-active enzymes studied by comparative molecular dynamics simulation

Enzymes from cold-adapted species are significantly more active at low temperatures, even those close to zero Celsius, but the rationale of this adaptation is complex and relatively poorly understood. It is commonly stated that there is a relationship between the flexibility of an enzyme and its catalytic activity at low temperature. This paper gives the results of a study using molecular dynamics simulations performed for five pairs of enzymes, each pair comprising a cold-active enzyme plus its mesophilic or thermophilic counterpart. The enzyme pairs included α-amylase, citrate synthase, malate dehydrogenase, alkaline protease and xylanase. Numerous sites with elevated flexibility were observed in all enzymes; however, differences in flexibilities were not striking. Nevertheless, amino acid residues common in both enzymes of a pair (not present in insertions of a structure alignment) are generally more flexible in the cold-active enzymes. The further application of principle component analysis to the protein dynamics revealed that there are differences in the rate and/or extent of opening and closing of the active sites. The results indicate that protein dynamics play an important role in catalytic processes where structural rearrangements, such as those required for active site access by substrate, are involved. They also support the notion that cold adaptation may have evolved by selective changes in regions of enzyme structure rather than in global change to the whole protein. Figure Collective motions in C^α atoms of the active site of cold-active xylanase Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Kinetic and thermodynamic properties of wild-type and engineered mutants of tyrosyl-tRNA synthetase analyzed by pyrophosphate-exchange kinetics.

The first step of the reaction catalyzed by the aminoacyl-tRNA synthetases is the formation of enzyme-bound aminoacyl adenylate. The steady-state kinetics of this step has conventionally been studied by measuring the rate of isotopic exchange between pyrophosphate and ATP. A simple kinetic analysis of the pyrophosphate-exchange reaction catalyzed by the tyrosyl-tRNA synthetase from Bacillus stearothermophilus is given in which all the observed rate and binding constants can be assigned to identifiable physical processes under a variety of limiting conditions. The free energies of binding to the enzyme of tyrosine, ATP, and the transition state for tyrosyl adenylate formation can be measured in relatively straightforward experiments. The excellent agreement between parameters measured in these experiments and those from earlier pre-steady-state kinetics confirms that the intermediates isolated in the presteady state are kinetically competent. The dissociation constant of ATP from the unligated enzyme, a constant that has previously been experimentally inaccessible, has been measured for wild-type and several mutant enzymes. The changes in enthalpy and entropy of activation on mutation have been measured by a rapid procedure for mutants that have altered contacts with tyrosine and ATP. Those mutants that have large changes of enthalpy and entropy of binding are likely to have structural changes and so warrant further examination by protein crystallography.     

Importin-beta: structural and dynamic determinants of a molecular spring

The beta-karyopherin/RanGTP system constitutes the largest known family of cellular cargo transporters. The flexibility of the karyopherin transport receptors is the key to their versatility in binding cargoes of different shape and size. Despite strong binding of the Ran complex, the comparably low energy associated with GTP hydrolysis suffices to drive dissociation and fuel the transport cycle. Here, we elucidate the drastic structural dynamics of the prototypic karyopherin, importin-beta, and show that its flexibility also solves this energetic puzzle. Our nonequilibrium atomistic simulations reveal fast conformational changes, validated by small-angle X-ray scattering data, and unusually large structural fluctuations. The characteristic dynamic patterns of importin-beta and the observed unfolding pathway of the IBB domain suggest a cooperative mechanism of importin-beta function in the nucleus. We propose a molecular model in which the stored energy and structural dynamics account for an exchange pathway that explains the high observed rates of nucleocytoplasmic transport. Karyopherins utilize a mechanism of entropy/enthalpy control that might be a general feature of highly flexible proteins involved in protein-protein interactions.     

Effects of anions on the activation thermodynamics and fluorescence emission spectrum of alkaline phosphatase: evidence for enzyme hydration changes during catalysis

The effect of anions on the thermodynamic activation functions for a model enzyme, calf intestinal alkaline phosphatase (EC 3.1.3.1), have been studied in order to examine the role of protein hydration changes in establishing the energetics of enzyme catalysis. The influences of these anions on the activation volume (delta V) and activation free energy (delta G) reflected clear Hofmeister (lyotropic) series effects, in the order F- greater than Cl- greater than Br- greater than I- (order of increasing salting-out potential). A pronounced covariation was observed between the influences of these anions on Vmax, which is proportional to delta G, and on the negative activation volume of the reaction. Fluoride was able to counteract the influences of Br- and I- on both Vmax and delta V when combinations of these anions were employed. The effects of Br- and I- on Vmax and delta V were more pronounced at lower temperatures. The control delta V was increasingly negative at reduced temperatures. The effects of the neutral salts and propanol on delta V and delta G, as well as the effects of salting-in anions on the activation enthalpy and the negative activation entropy of the reaction, are consistent with a model which proposes that peptide groups or polar side chains on the native enzyme exergonically increase their exposure to solvent during the catalytic activation event. These conclusions are in accord with the known free energy, enthalpy, entropy, and volume changes which occur when model peptide groups are transferred between water and concentrated salt solutions. Consistent with the kinetic results, the fluorescence emission wavelength maximum of alkaline phosphatase increased in the presence of anions in the order F- greater than Cl- greater than Br- greater than I-. The salting-out ion (F-) and the salting-in ions (Br- and I-) shifted lambda max in different directions, and these lambda max shifts could be counterbalanced by using equimolar combinations of salting-in and salting-out anions. Control experiments with a model compound, N-acetyltryptophanamide, showed that the spectra shifts caused by the salts did not result solely from differential quenching by the anions of the solvent-exposed tryptophan(s) on the enzyme. Hofmeister additivity phenomena indicated that the solvent is at the basis of these salt-induced enzyme structural changes. It is concluded that changes in protein solvation during enzymic reactions contribute significantly to the thermodynamic activation parameters in both the native and the salt-perturbed enzyme.     

Entropy-driven mechanism of an E3 ligase

Ubiquitin-like modifications are macromolecular chemistry for which our understanding of the enzymatic mechanisms is lacking. Most E3 ligases in ubiquitin-like modifications do not directly participate in chemistry but are thought to confer allosteric effects; however, the nature of the allosteric effects has been elusive. Recent molecular dynamics simulations suggested that an E3 binding enhances the population of the conformational states of the E2·SUMO thioester that favor reactions. In this study, we conducted the first temperature-dependent enzyme kinetic analysis to investigate the role of an E3 on activation entropy and enthalpy. The small ubiquitin-like modifier (SUMO) E3, RanBP2, confers unusually large, favorable activation entropy to lower the activation energy of the reaction. Mutants of RanBP2, designed to alter the flexibilities of the E2·SUMO thioester, showed a direct correlation of their favorable entropic effects with their ability to restrict the conformational flexibility of the E2·SUMO thioester. While the more favorable activation entropy is consistent with the previously suggested role of E3 in conformational selection, the large positive entropy suggests a significant role of solvent in catalysis. Indeed, molecular dynamics simulations in explicit water revealed that the more stable E2·SUMO thioester upon E3 binding results in stabilization of a large number of bound water molecules. Liberating such structured water at the transition state can result in large favorable activation entropy but unfavorable activation enthalpy. The entropy-driven mechanism of the E3 is consistent with the lack of structural conservation among E3s despite their similar functions. This study also illustrates how proteins that bind both SUMO and E2 can function as E3s and how intrinsically unstructured proteins can enhance macromolecular chemistry in addition to their known advantages in protein--protein interactions.     

Kinetics of the interaction of dihydroalprenolol with beta-adrenergic receptors in rat cerebral cortex

Time and temperature dependence of the binding of 3H-dihydroalprenolol (3H-DHA) to beta-adrenergic receptors in rat cerebral cortex is described. The kinetic data obtained suggest that 3H-DHA binding proceeds through a two-step reaction scheme consisting of a bimolecular association step followed by an unimolecular internal conversion of the radioligand receptor complex (isomerisation). Equilibrium thermodynamic analysis provided evidence that the over-all binding process is associated with a small decrease in enthalpy and a substantial increase in entropy. Within the framework of the two-step binding kinetics, the evaluation of the temperature dependence by the van't Hoff analysis resulted in values for thermodynamic parameters for the single equilibrium steps. The data suggest that the association step can be considered as a bimolecular hydrophobic interaction which is mainly entropy-driven due to the release of structural water, while the isomerisation step is accompanied by a large negative change in both enthalpy and entropy. The large negative change in the activation entropy for the forward reaction of the isomerisation step, obtained from evaluation of Arrhenius plots, indicates an internal conversion to a highly ordered receptor-ligand complex, while the low activation energy points to a small threshold energy for reaching this structure. Thus, these result support a previous assumption that the hydrophobic center of an adrenergic antagonist interacts with the receptor by entering a pocket (Cherksey et al. 1981).