Vibrational Structure of NAD(P) Cofactors Bound to Three NAD(P) Dependent Enzymes: an Investigation of Ground State Activation

Allostery and Evolution: A Molecular Journey Through the Structural and Dynamical Landscape of an Enzyme Super Family

Allosteric regulation is a powerful mechanism for controlling the efficiency of enzymes. Deciphering the evolutionary mechanisms by which allosteric properties have been acquired in enzymes is of fundamental importance. We used the malate (MalDH) and lactate deydrogenases (LDHs) superfamily as model to elucidate this phenomenon. By introducing a few of mutations associated to the emergence of allosteric LDHs into the non-allosteric MalDH from Methanopyrus kandleri, we have gradually shifted its enzymatic profile toward that typical of allosteric LDHs. We first investigated the process triggering homotropic activation. The structures of the resulting mutants show the typical compact organization of the R-active state of LDHs, but a distorted (T-like) catalytic site demonstrating that they corresponds to hybrid states. Molecular dynamics simulations and free energy calculations confirmed the capability of these mutants to sample the T-inactive state. By adding a final single mutation to fine-tune the flexibility of the catalytic site, we obtained an enzyme with both sigmoid (homotropic) and hyperbolic (heterotropic) substrate activation profiles. Its structure shows a typical extended T-state as in LDHs, whereas its catalytic state has as a restored configuration favorable for catalysis. Free energy calculations indicate that the T and R catalytic site configurations are in an equilibrium that depends on solvent conditions. We observed long-range communication between monomers as required for allosteric activation. Our work links the evolution of allosteric regulation in the LDH/MDH superfamily to the ensemble model of allostery at molecular level, and highlights the important role of the underlying protein dynamics.

NAD(P)H binding configurations revealed by time-resolved fluorescence and two-photon absorption

NADH and NADPH play key roles in the regulation of metabolism. Their endogenous fluorescence is sensitive to enzyme binding, allowing changes in cellular metabolic state to be determined using fluorescence lifetime imaging microscopy (FLIM). However, to fully uncover the underlying biochemistry, the relationships between their fluorescence and binding dynamics require greater understanding. Here we accomplish this through time- and polarization-resolved fluorescence and polarized two-photon absorption measurements. Two lifetimes result from binding of both NADH to lactate dehydrogenase and NADPH to isocitrate dehydrogenase. The composite fluorescence anisotropy indicates the shorter (1.3-1.6 ns) decay component to be accompanied by local motion of the nicotinamide ring, pointing to attachment solely via the adenine moiety. For the longer lifetime (3.2-4.4 ns), the nicotinamide conformational freedom is found to be fully restricted. As full and partial nicotinamide binding are recognized steps in dehydrogenase catalysis, our results unify photophysical, structural, and functional aspects of NADH and NADPH binding and clarify the biochemical processes that underlie their contrasting intracellular lifetimes.

Mechanism of Thermal Adaptation in the Lactate Dehydrogenases

The mechanism of thermal adaptation of enzyme function at the molecular level is poorly understood but is thought to lie within the structure of the protein or its dynamics. Our previous work on pig heart lactate dehydrogenase (phLDH) has determined very high resolution structures of the active site, via isotope edited IR studies, and has characterized its dynamical nature, via laser-induced temperature jump (T-jump) relaxation spectroscopy on the Michaelis complex. These particular probes are quite powerful at getting at the interplay between structure and dynamics in adaptation. Hence, we extend these studies to the psychrophilic protein cgLDH (Champsocephalus gunnari; 0 °C) and the extreme thermophile tmLDH (Thermotoga maritima LDH; 80 °C) for comparison to the mesophile phLDH (38-39 °C). Instead of the native substrate pyruvate, we utilize oxamate as a nonreactive substrate mimic for experimental reasons. Using isotope edited IR spectroscopy, we find small differences in the substate composition that arise from the detailed bonding patterns of oxamate within the active site of the three proteins; however, we find these differences insufficient to explain the mechanism of thermal adaptation. On the other hand, T-jump studies of reduced β-nicotinamide adenine dinucleotide (NADH) emission reveal that the most important parameter affecting thermal adaptation appears to be enzyme control of the specific kinetics and dynamics of protein motions that lie along the catalytic pathway. The relaxation rate of the motions scale as cgLDH > phLDH > tmLDH in a way that faithfully matches kcat of the three isozymes.

The dynamical nature of enzymatic catalysis

As is well-known, enzymes are proteins designed to accelerate specific life essential chemical reactions by many orders of magnitude. A folded protein is a highly dynamical entity, best described as a hierarchy or ensemble of interconverting conformations on all time scales from femtoseconds to minutes. We are just beginning to learn what role these dynamics play in the mechanism of chemical catalysis by enzymes due to extraordinary difficulties in characterizing the conformational space, that is, the energy landscape, of a folded protein. It seems clear now that their role is crucially important. Here we discuss approaches, based on vibrational spectroscopies of various sorts, that can reveal the energy landscape of an enzyme–substrate (Michaelis) complex and decipher which part of the typically very complicated landscape is relevant to catalysis. Vibrational spectroscopy is quite sensitive to small changes in bond order and bond length, with a resolution of 0.01 Å or less. It is this sensitivity that is crucial to its ability to discern bond reactivity.

Using isotope edited IR approaches, we have studied in detail the role of conformational heterogeneity and dynamics in the catalysis of hydride transfer by LDH (lactate dehydrogenase). Upon the binding of substrate, the LDH·substrate system undergoes a search through conformational space to find a range of reactive conformations over the microsecond to millisecond time scale. The ligand is shuttled to the active site via first forming a weakly bound enzyme·ligand complex, probably consisting of several heterogeneous structures. This complex undergoes numerous conformational changes spread throughout the protein that shuttle the enzyme·substrate complex to a range of conformations where the substrate is tightly bound. This ensemble of conformations all have a propensity toward chemistry, but some are much more facile for carrying out chemistry than others. The search for these tightly bound states is clearly directed by the forces that the protein can bring to bear, very much akin to the folding process to form native protein in the first place. In fact, the conformational subspace of reactive conformations of the Michaelis complex can be described as a “collapse” of reactive substates compared with that found in solution, toward a much smaller and much more reactive set.

These studies reveal how dynamic disorder in the protein structure can modulate the on-enzyme reactivity. It is very difficult to account for how the dynamical nature of the ground state of the Michaelis complex modulates function by transition state concepts since dynamical disorder is not a starting feature of the theory. We find that dynamical disorder may well play a larger or similar sized role in the measured Gibbs free energy of a reaction compared with the actual energy barrier involved in the chemical event. Our findings are broadly compatible with qualitative concepts of evolutionary adaptation of function such as adaptation to varying thermal environments. Our work suggests a methodology to determine the important dynamics of the Michaelis complex.

Energy landscape of the Michaelis complex of lactate dehydrogenase: relationship to catalytic mechanism

Lactate dehydrogenase (LDH) catalyzes the interconversion between pyruvate and lactate with nicotinamide adenine dinucleotide (NAD) as a cofactor. Using isotope-edited difference Fourier transform infrared spectroscopy on the “live” reaction mixture (LDH·NADH·pyruvate ⇌ LDH·NAD⁺·lactate) for the wild-type protein and a mutant with an impaired catalytic efficiency, a set of interconverting conformational substates within the pyruvate side of the Michaelis complex tied to chemical activity is revealed. The important structural features of these substates include (1) electronic orbital overlap between pyruvate’s C₂=O bond and the nicotinamide ring of NADH, as shown from the observation of a delocalized vibrational mode involving motions from both moieties, and (2) a characteristic hydrogen bond distance between the pyruvate C₂=O group and active site residues, as shown by the observation of at least four C₂=O stretch bands indicating varying degrees of C₂=O bond polarization. These structural features form a critical part of the expected reaction coordinate along the reaction path, and the ability to quantitatively determine them as well as the substate population ratios in the Michaelis complex provides a unique opportunity to probe the structure–activity relationship in LDH catalysis. The various substates have a strong variance in their propensity toward on enzyme chemistry. Our results suggest a physical mechanism for understanding the LDH-catalyzed chemistry in which the bulk of the rate enhancement can be viewed as arising from a stochastic search through an available phase space that, in the enzyme system, involves a restricted ensemble of more reactive conformational substates as compared to the same chemistry in solution.

Enzyme active site interactions by Raman/FTIR, NMR, and ab initio calculations

Characterization of enzyme active site structure and interactions at high resolution is important for the understanding of the enzyme catalysis. Vibrational frequency and NMR chemical shift measurements of enzyme-bound ligands are often used for such purpose when X-ray structures are not available or when higher resolution active site structures are desired. This review is focused on how ab initio calculations may be integrated with vibrational and NMR chemical shift measurements to quantitatively determine high-resolution ligand structures (up to 0.001 Å for bond length and 0.01 Å for hydrogen bonding distance) and how interaction energies between bound ligand and its surroundings at the active site may be determined. Quantitative characterization of substrate ionic states, bond polarizations, tautomeric forms, conformational changes and its interactions with surroundings in enzyme complexes that mimic ground state or transition state can provide snapshots for visualizing the substrate structural evolution along enzyme-catalyzed reaction pathway. Our results have shown that the integration of spectroscopic studies with theoretical computation greatly enhances our ability to interpret experimental data and significantly increases the reliability of the theoretical analysis.

Conformational heterogeneity within the Michaelis complex of lactate dehydrogenase

A series of isotope edited IR measurements, both static as well as temperature jump relaxation spectroscopy, are performed on lactate dehydrogenase (LDH) to determine the ensemble of structures available to its Michaelis complex. There clearly has been a substantial reduction in the number of states available to the pyruvate substrate (as modeled by the substrate mimic, oxamate) and NADH when bound to protein compared to dissolved in solution, as determined by the bandwidths and positions of the critical C(2)═O band of the bound substrate mimic and the C(4)-H stretch of the NADH reduced nicotinamide group. Moreover, it is found that a strong ionic bond (characterized by a signature IR band discovered in this study) is formed between the carboxyl group of bound pyruvate with (presumably) Arg171, forming a strong "anchor" within the protein matrix. However, conformational heterogeneity within the Michaelis complex is found that has an impact on both catalytic efficiency and thermodynamics of the enzyme.

Probing the role of dynamics in hydride transfer catalyzed by lactate dehydrogenase

The dynamic nature of the interconversion of pyruvate to lactate as catalyzed by lactate dehydrogenase (LDH) is characterized by laser-induced temperature jump relaxation spectroscopy with a resolution of 20 ns. An equilibrium system of LDH.NADH plus pyruvate and LDH.NAD+ plus lactate is perturbed by a sudden T-jump, and the relaxation of the system is monitored by NADH emission and absorption changes. The substrate binding pathway is observed to be similar, although not identical, to previous work on substrate mimics: an encounter complex is formed between LDH.NADH and pyruvate, which collapses to the active Michaelis complex. The previously unresolved hydride transfer event is characterized and separated from other unimolecular isomerizations of the protein important for the catalytic mechanism, such as loop closure, a slower step, and faster events on the nanosecond-microsecond timescales whose structural basis is not understood. The results of this study show that this approach can be applied quite generally to enzyme systems and report on the dynamic nature of proteins over a very wide time range.

On the pathway of forming enzymatically productive ligand-protein complexes in lactate dehydrogenase

We have carried out a series of studies on the binding of a substrate mimic to the enzyme lactate dehydrogenase (LDH) using advanced kinetic approaches, which begin to provide a molecular picture of the dynamics of ligand binding for this protein. Binding proceeds via a binding-competent subpopulation of the nonligated form of the protein (the LDH/NADH binary complex) to form a protein-ligand encounter complex. The work here describes the collapse of the encounter complex to form the catalytically competent Michaelis complex. Isotope-edited static Fourier transform infrared studies on the bound oxamate protein complex reveal two kinds of oxamate environments: 1), a major populated structure wherein all significant hydrogen-bonding patterns are formed at the active site between protein and bound ligand necessary for the catalytically productive Michaelis complex and 2), a minor structure in a configuration of the active site that is unfavorable to carry out catalyzed chemistry. This latter structure likely simulates a dead-end complex in the reaction mixture. Temperature jump isotope-edited transient infrared studies on the binding of oxamate with LDH/NADH suggest that the evolution of the encounter complex between LDH/NADH and oxamate collapses via a branched reaction pathway to form the major and minor bound species. The production of the catalytically competent protein-substrate complex has strong similarities to kinetic pathways found in two-state protein folding processes. Once the encounter complex is formed between LDH/NADH and substrate, the ternary protein-ligand complex appears to "fold" to form a compact productive complex in an all or nothing like fashion with all the important molecular interactions coming together at the same time.

The approach to the Michaelis complex in lactate dehydrogenase: the substrate binding pathway

We examine here the dynamics of forming the Michaelis complex of the enzyme lactate dehydrogenase by characterizing the binding kinetics and thermodynamics of oxamate (a substrate mimic) to the binary lactate dehydrogenase/NADH complex over multiple timescales, from nanoseconds to tens of milliseconds. To access such a wide time range, we employ standard stopped-flow kinetic approaches (slower than 1 ms) and laser-induced temperature-jump relaxation spectroscopy (10 ns-10 ms). The emission from the nicotinamide ring of NADH is used as a marker of structural transformations. The results are well explained by a kinetic model that has binding taking place via a sequence of steps: the formation of an encounter complex in a bimolecular step followed by two unimolecular transformations on the microsecond/millisecond timescales. All steps are well described by single exponential kinetics. It appears that the various key components of the catalytically competent architecture are brought together as separate events, with the formation of strong hydrogen bonding between active site His(195) and substrate early in binding and the closure of the catalytically necessary protein surface loop over the bound substrate as the final event of the binding process. This loop remains closed during the entire period that chemistry takes place for native substrates; however, motions of other key molecular groups bringing the complex in and out of catalytic competence appear to occur on faster timescales. The on-enzyme K(d) values (the ratios of the microscopic rate constants for each unimolecular step) are not far from one. Either substantial, approximately 10-15%, transient melting of the protein or rearrangements of hydrogen bonding and solvent interactions of a number of water molecules or both appear to take place to permit substrate access to the protein binding site. The nature of activating the various steps in the binding process seems to be one overall involving substantial entropic changes.