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

Computational Analysis of Heat Capacity Effects in Protein-Ligand Binding

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

Cooperative Conformational Transitions Underpin the Activation Heat Capacity in the Temperature Dependence of Enzyme Catalysis

Many enzymes display non-Arrhenius behavior with curved Arrhenius plots in the absence of denaturation. There has been significant debate about the origin of this behavior and recently the role of the activation heat capacity (ΔCP⧧) has been widely discussed. If enzyme-catalyzed reactions occur with appreciable negative values of ΔCP⧧ (arising from narrowing of the conformational space along the reaction coordinate), then curved Arrhenius plots are a consequence. To investigate these phenomena in detail, we have collected high precision temperature-rate data over a wide temperature interval for a model glycosidase enzyme MalL, and a series of mutants that change the temperature-dependence of the enzyme-catalyzed rate. We use these data to test a range of models including macromolecular rate theory (MMRT) and an equilibrium model. In addition, we have performed extensive molecular dynamics (MD) simulations to characterize the conformational landscape traversed by MalL in the enzyme-substrate complex and an enzyme-transition state complex. We have crystallized the enzyme in a transition state-like conformation in the absence of a ligand and determined an X-ray crystal structure at very high resolution (1.10 Å). We show (using simulation) that this enzyme-transition state conformation has a more restricted conformational landscape than the wildtype enzyme. We coin the term "transition state-like conformation (TLC)" to apply to this state of the enzyme. Together, these results imply a cooperative conformational transition between an enzyme-substrate conformation (ES) and a transition-state-like conformation (TLC) that precedes the chemical step. We present a two-state model as an extension of MMRT (MMRT-2S) that describes the data along with a convenient approximation with linear temperature dependence of the activation heat capacity (MMRT-1L) that can be used where fewer data points are available. Our model rationalizes disparate behavior seen for MalL and previous results for a thermophilic alcohol dehydrogenase and is consistent with a raft of data for other enzymes. Our model can be used to characterize the conformational changes required for enzyme catalysis and provides insights into the role of cooperative conformational changes in transition state stabilization that are accompanied by changes in heat capacity for the system along the reaction coordinate. TLCs are likely to be of wide importance in understanding the temperature dependence of enzyme activity and other aspects of enzyme catalysis.

Glycosylation increases active site rigidity leading to improved enzyme stability and turnover

Glycosylation is the most prevalent protein post-translational modification, with a quarter of glycosylated proteins having enzymatic properties. Yet, the full impact of glycosylation on the protein structure-function relationship, especially in enzymes, is still limited. Here, we show that glycosylation rigidifies the important commercial enzyme horseradish peroxidase (HRP), which in turn increases its turnover and stability. Circular dichroism spectroscopy revealed that glycosylation increased holo-HRP's thermal stability and promoted significant helical structure in the absence of haem (apo-HRP). Glycosylation also resulted in a 10-fold increase in enzymatic turnover towards o-phenylenediamine dihydrochloride when compared to its nonglycosylated form. Utilising a naturally occurring site-specific probe of active site flexibility (Trp117) in combination with red-edge excitation shift fluorescence spectroscopy, we found that glycosylation significantly rigidified the enzyme. In silico simulations confirmed that glycosylation largely decreased protein backbone flexibility, especially in regions close to the active site and the substrate access channel. Thus, our data show that glycosylation does not just have a passive effect on HRP stability but can exert long-range effects that mediate the 'native' enzyme's activity and stability through changes in inherent dynamics.

Immobilization of Kluyveromyces lactis β-Galactosidase on Meso-macroporous Silica: Use of Infrared Spectroscopy to Rationalize the Catalytic Efficiency

Enzyme immobilization on adequate carriers is a challenging strategy. Understanding the enzyme-carrier interactions and their effects on enzyme conformation and bioactivity is critical. In this study, a meso-macropores silica (MMS) was used to immobilize β-galactosidase from the yeast Kluyveromyces lactis (β-gal-KL) by physical adsorption. The bioactivity of the immobilized β-gal-KL was altered, evidenced by the increased K_m , decreased V_max and k_cat , and increased activity at alkaline values. By performing infrared spectroscopy analysis and subsequent secondary structure assessment from the amide I band, the immobilized β-gal-KL suffered a β-sheet (∼31-35 %) to α-helix (∼15-19 %) transition with increased turns (∼21-22 %) with respect to the free β-gal-KL having ∼12 % α-helix, ∼42 % β-sheet, and ∼17 % turns. These findings led us to correlate the observed bioactivity performance to structural alterations to a non-native conformation. The presented line of thought can lead to a better understanding of the reasons causing bioactivity alterations upon enzyme immobilization.

Functional roles of enzyme dynamics in accelerating active site chemistry: Emerging techniques and changing concepts

With the growing acceptance of the contribution of protein conformational ensembles to enzyme catalysis and regulation, research in the field of protein dynamics has shifted toward an understanding of the atomistic properties of protein dynamical networks and the mechanisms and time scales that control such behavior. A full description of an enzymatic reaction coordinate is expected to extend beyond the active site and include site-specific networks that communicate with the protein/water interface. Advances in experimental tools for the spatial resolution of thermal activation pathways are being complemented by biophysical methods for visualizing dynamics in real time. An emerging multidimensional model integrates the impacts of bound substrate/effector on the distribution of protein substates that are in rapid equilibration near room temperature with reaction-specific protein embedded heat transfer conduits.