A new intrinsic thermal parameter for enzymes reveals true temperature optima

Two established thermal properties of enzymes are the Arrhenius activation energy and thermal stability. Arising from anomalies found in the variation of enzyme activity with temperature, a comparison has been made of experimental data for the activity and stability properties of five different enzymes with theoretical models. The results provide evidence for a new and fundamental third thermal parameter of enzymes, T(eq), arising from a subsecond timescale-reversible temperature-dependent equilibrium between the active enzyme and an inactive (or less active) form. Thus, at temperatures above its optimum, the decrease in enzyme activity arising from the temperature-dependent shift in this equilibrium is up to two orders of magnitude greater than what occurs through thermal denaturation. This parameter has important implications for our understanding of the connection between catalytic activity and thermostability and of the effect of temperature on enzyme reactions within the cell. Unlike the Arrhenius activation energy, which is unaffected by the source ("evolved") temperature of the enzyme, and enzyme stability, which is not necessarily related to activity, T(eq) is central to the physiological adaptation of an enzyme to its environmental temperature and links the molecular, physiological, and environmental aspects of the adaptation of life to temperature in a way that has not been described previously. We may therefore expect the effect of evolution on T(eq) with respect to enzyme temperature/activity effects to be more important than on thermal stability. T(eq) is also an important parameter to consider when engineering enzymes to modify their thermal properties by both rational design and by directed enzyme evolution.     

Temperature and the catalytic activity of enzymes: A fresh understanding

The discovery of an additional step in the progression of an enzyme from the active to inactive state under the influence of temperature has led to a better match with experimental data for all enzymes that follow Michaelis–Menten kinetics, and to an increased understanding of the process. The new model of the process, the Equilibrium Model, describes an additional mechanism by which temperature affects the activity of enzymes, with implications for ecological, metabolic, structural, and applied studies of enzymes.     

High thermal stability of 3-isopropylmalate dehydrogenase from Thermus thermophilus resulting from low DeltaC(p) of unfolding.

To characterize the thermal stability of 3-isopropylmalate dehydrogenase (IPMDH) from an extreme thermophile, Thermus thermophilus, urea-induced unfolding of the enzyme and of its mesophilic counterpart from Escherichia coli was investigated at various temperatures. The unfolding curves were analyzed with a three-state model for E.coli IPMDH and with a two-state model for T.thermophilus IPMDH, to obtain the free energy change DeltaG degrees of each unfolding process. Other thermodynamic parameters, enthalpy change DeltaH, entropy change DeltaS and heat capacity change DeltaC(p), were derived from the temperature dependence of DeltaG degrees. The main feature of the thermophilic enzyme was its lower dependence of DeltaG degrees on temperature resulting from a low DeltaC(p). The thermophilic IPMDH had a significantly lower DeltaC(p), 1.73 kcal/mol.K, than that of E.coli IPMDH (20.7 kcal/mol.K). The low DeltaC(p) of T.thermophilus IPMDH could not be predicted from its change in solvent-accessible surface area DeltaASA. The results suggested that there is a large structural difference between the unfolded state of T.thermophilus and that of E.coli IPMDH. Another responsible factor for the higher thermal stability of T.thermophilus IPMDH was the increase in the most stable temperature T(s). The DeltaG degrees maximum of T.thermophilus IPMDH was much smaller than that of E.coli IPMDH. The present results clearly demonstrated that a higher melting temperature T(m) is not necessarily accompanied by a higher DeltaG degrees maximum.     

Allostery in trypsin-like proteases suggests new therapeutic strategies

Trypsin-like proteases (TLPs) are a large family of enzymes responsible for digestion, blood coagulation, fibrinolysis, development, fertilization, apoptosis and immunity. A current paradigm posits that the irreversible transition from an inactive zymogen to the active protease form enables productive interaction with substrate and catalysis. Analysis of the entire structural database reveals two distinct conformations of the active site: one fully accessible to substrate (E) and the other occluded by the collapse of a specific segment (E*). The allosteric E*-E equilibrium provides a reversible mechanism for activity and regulation in addition to the irreversible zymogen to protease conversion and points to new therapeutic strategies aimed at inhibiting or activating the enzyme. In this review, we discuss relevant examples, with emphasis on the rational engineering of anticoagulant thrombin mutants.     

Radioiodinated antibodies,a tool in studies on the presence and role of inactive enzyme forms: Regulation of chalcone synthase in parsley cell suspension cultures

A new technique, the quantitative determination of total enzyme concentrations by specific immunoprecipitation with purified, radioiodinated antibodies, was used to investigate the presence and possible roles of inactive enzyme in the regulation of chalcone synthase. Dark-grown cell suspension cultures from parsley (Petroselinum hortense) contained neither catalytically active nor detectable amounts of immunoprecipitable chalcone synthase. Irradiation induced large increases and subsequent decreases of both. Significant differences in the peak positions and in the half-lives of active and total chalcone synthase indicated that induced cells contained inactive as well as active enzyme forms. The presence of inactive enzyme could be explained by two different modes of regulation, (i) simultaneous de novo synthesis of active and inactive enzyme (“Simultaneous Model”), or (ii) de novo synthesis of active enzyme only, with sequential steps of inactivation and degradation (“Sequential Model”). Both models were compatible with experimental results, as analyzed mathematically by investigating the relations between curves for rate of enzyme synthesis, enzyme activity, total enzyme, and half-lives of active and total enzyme. However, the “Simultaneous Model” postulated that de novo synthesis of inactive enzyme represented always the vast majority of total enzyme synthesis, while the Sequential Model integrated inactive enzyme with facility in a sequence of irreversible inactivation and degradation of active enzyme. Experiments with repeated induction indicated that cells containing large amounts of inactive enzyme increased enzyme activity by de novo synthesis rather than by activation of preexisting inactive enzyme.     

Effect of phospholipids on the thermal stability of microsomal UDP-glucuronosyltransferase

The GT2P isoform of microsomal UDP-glucuronosyltransferase from pig liver is a lipid-dependent enzyme. The data in the present work indicate that, in addition to regulation of activity, the thermal stability of the enzyme also is modulated by the acyl chain composition of phosphatidylcholines (PC) used to reconstitute the activity of pure enzyme. There was a reversible, temperature-dependent change in the state of the pure enzyme to an inactive form with onset at T greater than 38 degrees C, depending on the environment of the enzyme. The midpoint for the transition shifted from 39.8 degrees C for enzyme in a bilayer of distearoylphosphatidylcholine (DSPC) to 47.5 degrees C for enzyme in a bilayer of 1-stearoyl-2-oleoylphosphatidylcholine (SOPC). For all lipids, the transition from a catalytically active to an inactive form of the enzyme was associated with large compensating changes in H and S. Lipid-induced stabilization of the active form of UDP-glucuronosyltransferase at T greater than 37 degrees C was associated with decreases in delta H and delta S, but the decreases in delta S were larger, indicating that lipid-induced stabilization of the active form of the enzyme was entropic. The transition between the active and inactive forms of the enzyme was too rapid in either direction to measure in a standard spectrophotometer. In addition to reversible inactivation of the enzyme, there was a slower irreversible, temperature-dependent inactivation. The rate of this process depended on the acyl chains of the phosphocholines interacting with the enzyme. However, there was no obvious correlation between the structures of lipids that stabilized the different inactivation reactions.     

Modeling the inhibitor activity and relative binding affinities in enzyme-inhibitor-protein systems: application to developmental regulation in a PG-PGIP system

Within a number of classes of hydrolytic enzymes are certain enzymes whose activity is modulated by a specific inhibitor-protein that binds to the enzyme and forms an inactive complex. One unit of a specific inhibitor-protein activity is often defined as the amount necessary to inhibit one unit of its target enzyme by 50 %. No objective quantitative means is available to determine this point of 50 % inhibition in crude systems such as those encountered during purification. Two models were derived: the first model is based on an irreversible binding approximation, and the second, or equilibrium, model is based on reversible binding. The two models were validated using the inhibition data for the polygalacturonase-polygalacturonase-inhibiting protein (PG-PGIP) system. Theory and experimental results indicate that the first model can be used for inhibitor protein activity determination and the second model can be used for inhibitor protein activity determination as well as for comparison of association constants among enzymes and their inhibitor-proteins from multiple sources. The models were used to identify and further clarify the nature of a differential regulation of expression of polygalacturonase-inhibiting protein in developing cantaloupe fruit. These are the first relations that provide for an objective and quantitative determination of inhibitor-protein activity in both pure and crude systems. Application of these models should prove valuable in gaining insights into regulatory mechanisms and enzyme-inhibitor-protein interactions.     

Mutational analysis of the subunit interface of Vibrio harveyi bacterial luciferase

Bacterial luciferase is a heterodimeric (alphabeta) enzyme which catalyzes a light-producing reaction in Vibrio harveyi. In addition to the alphabeta enzyme, the beta subunit can self-associate to form a stable but inactive homodimer [Sinclair, J. F., Ziegler, M. M., and Baldwin, T. O. (1994) Nat. Struct. Biol. 1, 320-326]. The studies reported here were undertaken to explore the role of the subunit interface in the conformational stability of the enzyme. To this end, we constructed four mutant heterodimers in which residues at the subunit interface were changed in an effort to alter the volume of an apparent solvent accessible channel at the interface or to alter H-bonding groups. Equilibrium unfolding data for the heterodimer have been interpreted in terms of a three-state mechanism [Clark, C. A., Sinclair, J. F., and Baldwin, T. O. (1993) J. Biol. Chem. 268, 10773-10779]. However, we found that unfolding for the wild-type and mutant luciferases is better described by a four-state model. This change in the proposed mechanism of unfolding is based on observation of residual structure in the subunits following dissociation of the heterodimeric intermediate. All of the mutants display modest reductions in activity but, surprisingly, no change in the DeltaG2H2O value for subunit dissociation and no measurable change in the equilibrium dissociation constant relative to that of the wild-type heterodimer. However, the DeltaG1H2O value for the formation of the dimeric intermediate that precedes subunit dissociation is reduced for three of the mutants, indicating that mutations at the interface can alter the stability of a region of the alpha subunit that is distant from the interface. We conclude that the interface region communicates with the distal domains of this subunit, probably through the active center region of the enzyme.     

Mechanism of the phospholipid transfer protein-mediated transfer of phospholipids from model lipid vesicles to high density lipoproteins

To study the effects of the phospholipid transfer protein (PLTP) on the thermodynamic parameters governing the transfer of phospholipids (PL) from single bilayer vesicles (SBV) to high density lipoprotein (HDL), we performed transfer measurements at various temperatures between 4 and 65 degrees C, using a pyrenylphosphatidylcholine (Pyr-PC) as probe. The proportion of excimer (E) to monomer (M) fluorescence of a pyrenyl moiety constitutes a direct measure of its local concentration. The transfers of Pyr-PC were monitored by following the decrease of E/M. The data were used to calculate the rate constants K(+1) for the transfer from SBV to HDL and to generate the corresponding Arrhenius plots. The equilibrium constants, K(eq), for the same reactions were also determined and used to generate Van't Hoff plots. From these data, we calculated the thermodynamic parameters for both the whole transfer reaction and the transition state. Both K(+1) and K(eq) values clearly varied with temperature. PLTP induced very similar decreases in the free energy for the whole reaction (DeltaG) and in that for the transition state (DeltaG(#)). At 37 degrees C, the decreases were of 0.37 and 0.29 kcal/mol, respectively. We studied the thermal denaturation of PLTP between 37 and 65 degrees C, and the effects of denatured PLTP samples on the PL transfer reaction were then determined. In all cases, the changes of DeltaG remained comparable to those of DeltaG(#). Thus the essential action of PLTP is to facilitate the first step of the reaction, which can be considered as the desorption of PL molecules from the surface of donor particles.     

Thermal stability of acetylcholinesterase from Bungarus fasciatus venom as investigated by capillary electrophoresis

Previous studies on the conformation of the monomeric acetylcholinesterase (AChE) from the krait (Bungarus fasciatus) venom showed that the protein possesses a large permanent dipole moment. These studies predicted that thermal irreversible denaturation must occur via partially unfolded states. The thermal stability of Bungarus AChE was determined using capillary electrophoresis (CE) with optimized conditions. Runs performed at convenient temperature scanning rates provided evidence for an irreversible denaturation process according to the Lumry and Eyring model. The mid-transition temperature, T(m), and the effective enthalpy change, DeltaH(m) were determined at different pH. The temperature dependence of the free energy, DeltaG, of Bungarus AChE unfolding was drawn using values of T(m), DeltaH(m) and DeltaC(p) determined by CE. The thermodynamic parameters for the thermal denaturation of the monomeric snake enzyme were compared with those of different dimeric and tetrameric ChEs. It was shown that the changes in the ratio of DeltaH(cal/)DeltaH(vH) and DeltaC(p) reflect the oligomerization state of these proteins. All these results indicate that wild-type monomeric Bungarus AChE is a stable enzyme under standard conditions. However, designed mutants of this enzyme capable of degrading organophosphates have to be engineered to enhance their thermostability.     

A thermodynamic molecular switch in biological systems: ribonuclease S' fragment complementation reactions

It is well known that essentially all biological systems function over a very narrow temperature range. Most typical macromolecular interactions show DeltaH degrees (T) positive (unfavorable) and a positive DeltaS degrees (T) (favorable) at low temperature, because of a positive (DeltaCp degrees /T). Because DeltaG degrees (T) for biological systems shows a complicated behavior, wherein DeltaG degrees (T) changes from positive to negative, then reaches a negative value of maximum magnitude (favorable), and finally becomes positive as temperature increases, it is clear that a deeper-lying thermodynamic explanation is required. This communication demonstrates that the critical factor is a temperature-dependent DeltaCp degrees (T) (heat capacity change) of reaction that is positive at low temperature but switches to a negative value at a temperature well below the ambient range. Thus the thermodynamic molecular switch determines the behavior patterns of the Gibbs free energy change and hence a change in the equilibrium constant, K(eq), and/or spontaneity. The subsequent, mathematically predictable changes in DeltaH degrees (T), DeltaS degrees (T), DeltaW degrees (T), and DeltaG degrees (T) give rise to the classically observed behavior patterns in biological reactivity, as may be seen in ribonuclease S' fragment complementation reactions.     

Structural equilibrium fluctuations in mesophilic and thermophilic alpha-amylase

By comparing a mesophilic alpha-amylase with its thermophilic homolog, we investigated the relationship between thermal stability and internal equilibrium fluctuations. Fourier transform infrared spectroscopy monitoring hydrogen/deuterium (H/D) exchange kinetics and incoherent neutron scattering measuring picosecond dynamics were used to study dynamic features of the folded state at room temperature. Fairly similar rates of slowly exchanging amide protons indicate about the same free energy of stabilization DeltaG(stab) for both enzymes at room temperature. With respect to motions on shorter time scales, the thermophilic enzyme is characterized by an unexpected higher structural flexibility as compared to the mesophilic counterpart. In particular, the picosecond dynamics revealed a higher degree of conformational freedom for the thermophilic alpha-amylase. The mechanism proposed for increasing thermal stability in the present case is characterized by entropic stabilization and by flattening of the curvature of DeltaG(stab) as a function of temperature.     

Stability of butyrylcholinesterase: thermal inactivation in water and deuterium oxide

Irreversible thermal inactivation of the tetrameric form of human plasma butyrylcholinesterase (cholinesterase; EC 3.1.1.8) was studied in water and in deuterium oxide at pH 7 in the temperature range 53-65 degrees C. The enzyme inactivation follows a complex kinetics that may be described by the sum of two apparent first-order processes. The Eyring plot for enzyme inactivation exhibits a wavelike discontinuity over a span of 2 C degrees around 58 degrees C. This transition was interpreted in terms of equilibrium between two temperature-dependent conformational states. Though 2H2O does not alter the overall multistep inactivation process, a slight solvent isotope effect was observed: a stabilizing effect and a shift in the transition temperature. A comparison between several enzyme preparations revealed differences in thermodynamic activation parameters of inactivation suggesting microheterogeneity in enzyme structures. Kinetics of inactivation of usual (E1uE1u) and atypical (E1aEa1a++) enzymes were compared. The atypical enzyme was found to be more stable than the usual phenotype.     

Conformational equilibrium of an enzyme catalytic site in the allosteric transition.

The dynamic equilibrium of a catalytic site between active and inactive conformations, the missing link between the structure and function of allosteric enzymes, was identified using protein engineering and NMR techniques. Kinetic analyses of the wild-type and three mutants of Thermus L-lactate dehydrogenase established that the allosteric property of the enzyme is associated with a concerted transition between the high-affinity (R) and low-affinity (T) states. By introducing mutations, we prepared an enzyme in which the R and T states were balanced. The conformation of the enzyme-bound coenzyme, NAD+, which interacts directly with the substrate, was analyzed using NMR spectroscopy. NAD+ bound to the mutant enzyme was in a conformational mixture of the active and inactive forms, while NAD+ took on predominantly one of the two forms when it was bound to the other enzymes we had analyzed. We interpret this to mean that the catalytic site is in equilibrium between the two conformations. The ratio of the conformers of each enzyme agreed with the [T]/[R] ratio as determined by kinetic analyses. Therefore, it is the identified conformational equilibrium of the catalytic site that governs the allosteric regulation of the enzyme activity.     

Purification, kinetic and thermodynamic studies of a new ribonuclease from a mutant of Aspergillus niger

Ribonuclease was purified from Aspergillus niger SA-13-20 to homogeneity level by using (NH(4))(2)SO(4) precipitation, DEAE-cellulose anion-exchange chromatography, ultrafiltration and Sephacryl HR-200 chromatography. The molecular weight and isoelectric point of the enzyme was 40.1kDa and 5.3, respectively. The pH- and temperature-dependent kinetic parameters were determined. The RNase showed the strongest affinity with RNA as the substrate, and the highest catalytic efficiency for hydrolysis of the substrate at pH 3.5 and 65 degrees C. It exhibited Michaelis-Menten Kinetics with k(cat) of 118.1s(-1) and K(m) of 57.0 microg ml(-1), respectively. Thermodynamic parameters for catalysis and thermal denaturation were also determined. Activation energy (E(a)) for catalysis of A. niger SA-13-20 RNase was 50.31 kJ mol(-1) and free energy (DeltaG(#)), enthalpy (DeltaH(#)) and entropy (DeltaS(#)) of activation for catalysis of the enzyme at 65 degrees C were 69.76, 47.50 and -65.83 Jmol(-1)K(-1), respectively. Activation energy (E(a,d)) for denaturation of the enzyme was 200.53 kJ mol(-1) and free energy (DeltaG(d)(#)), enthalpy (DeltaH(d)(#)) and entropy (DeltaS(d)(#)) of activation for denaturation of the enzyme at 45 degrees C were 79.18 kJ mol(-1), 197.88 and 373.09 Jmol(-1)K(-1), respectively.     

Protein ubiquitination is regulated by phosphorylation. An in vitro study.

Protein ubiquitination has been implicated in ATP-dependent protein turnover and in a number of biological processes in eukaryotic cells. The ubiquitination activating enzyme, E1, and ubiquitin carrier protein, E2, are two essential enzymes in the protein ubiquitination machinery. Using purified E1 and E2 from rabbit reticulocytes and various protein kinases, which include cAMP-dependent protein kinase, protein kinase C, and protein tyrosine kinase, we demonstrated that E1 is phosphorylated by protein kinase C, with a stoichiometry of 0.65 mol of phosphate/mol of E1, and one of the E2 isoforms, E2(32kDa), is phosphorylated by protein tyrosine kinase to 2 eq of phosphate/mol of protein. Phosphorylation of E1 causes a 2-fold enhancement of its activity as monitored by ubiquitin-dependent ATP in equilibrium PPi exchange. When 1 eq of phosphate was incorporated into E2(32kDa), a 2.4-fold activation was also observed for its activity to catalyze the ubiquitination of histone H2A. The regulatory significance of this finding is discussed.     

Thermostability of irreversible unfolding alpha-amylases analyzed by unfolding kinetics

For most multidomain proteins the thermal unfolding transitions are accompanied by an irreversible step, often related to aggregation at elevated temperatures. As a consequence the analysis of thermostabilities in terms of equilibrium thermodynamics is not applicable, at least not if the irreversible process is fast with respect the structural unfolding transition. In a comparative study we investigated aggregation effects and unfolding kinetics for five homologous alpha-amylases, all from mesophilic sources but with rather different thermostabilities. The results indicate that for all enzymes the irreversible process is fast and the precedent unfolding transition is the rate-limiting step. In this case the kinetic barrier toward unfolding, as measured by unfolding rates as function of temperature, is the key feature in thermostability. The investigated enzymes exhibit activation energies (E(a)) between 208 and 364 kJmol(-1) and pronounced differences in the corresponding unfolding rates. The most thermostable alpha-amylase from Bacillus licheniformis (apparent transition temperature, T(1/2) approximately 100 degrees C) shows an unfolding rate which is four orders of magnitude smaller as compared with the alpha-amylase from pig pancreas (T(1/2) approximately 65 degrees C). Even with respect to two other alpha-amylases from Bacillus species (T(1/2) approximately 86 degrees C) the difference in unfolding rates is still two orders of magnitude.     

Inactivation kinetics of a new target of beta-lactam antibiotics

Peptidoglycan is predominantly cross-linked by serine DD-transpeptidases in most bacterial species. The enzymes are the essential targets of β-lactam antibiotics. However, unrelated cysteine LD-transpeptidases have been recently recognized as a predominant mode of peptidoglycan cross-linking in Mycobacterium tuberculosis and as a bypass mechanism conferring resistance to all β-lactams, except carbapenems such as imipenem, in Enterococcus faecium. Investigation of the mechanism of inhibition of this new β-lactam target showed that acylation of the E. faecium enzyme (Ldt(fm)) by imipenem is irreversible. Using fluorescence kinetics, an original approach was developed to independently determine the catalytic constants for imipenem binding (k(1) = 0.061 μM(-1) min(-1)) and acylation (k(inact) = 4.5 min(-1)). The binding step was limiting at the minimal drug concentration required for bacterial growth inhibition. The Michaelis complex was committed to acylation because its dissociation was negligible. The emergence of imipenem resistance involved substitutions in Ldt(fm) that reduced the rate of formation of the non-covalent complex but only marginally affected the efficiency of the acylation step. The methods described in this study will facilitate development of new carbapenems active on extensively resistant M. tuberculosis.     

Determination of the free-energy change for repair of a DNA phosphodiester bond

The repair of phosphodiester bonds in nicked DNA is catalyzed by DNA ligases. Ligation is coupled to cleavage of a phosphoanhydride bond in a nucleotide cofactor resulting in a thermodynamically favorable process. A free energy value for phosphodiester bond formation was calculated using the reversibility of the T4 DNA ligase reaction. The relative number of DNA nicks to phosphodiester bonds in a circular plasmid DNA, formed during this reaction at fixed concentrations of ATP to AMP and PP(i), was quantified. At 25 degrees C, pH 7, the equilibrium constant (K(eq)) for the ligation reaction is 3.89 x 10(4) m. This value corresponds to a standard free energy (DeltaG degrees ') of -6.3 kcal mol(-1). By subtracting the known energy contribution due to hydrolysis of ATP to AMP and PP(i), DeltaG degrees ' for the hydrolysis of a DNA phosphodiester bond is -5.3 kcal mol(-1).