Structural transformations in the dynamics of Michaelis complex formation in lactate dehydrogenase

Resolution of Submillisecond Kinetics of Multiple Reaction Pathways for Lactate Dehydrogenase

Enzymes are known to exhibit conformational flexibility. An important consequence of this flexibility is that the same enzyme reaction can occur via multiple reaction pathways on a reaction landscape. A model enzyme for the study of reaction landscapes is lactate dehydrogenase. We have previously used temperature-jump (T-jump) methods to demonstrate that the reaction landscape of lactate dehydrogenase branches at multiple points creating pathways with varied reactivity. A limitation of this previous work is that the T-jump method makes only small perturbations to equilibrium and may not report conclusively on all steps in a reaction. Therefore, interpreting T-jump results of lactate dehydrogenase kinetics has required extensive computational modeling work. Rapid mixing methods offer a complementary approach that can access large perturbations from equilibrium; however, traditional enzyme mixing methods like stopped-flow do not allow for the observation of fast protein dynamics. In this report, we apply a microfluidic rapid mixing device with a mixing time of <100 μs that allows us to study these fast dynamics and the catalytic redox step of the enzyme reaction. Additionally, we report UV absorbance and emission T-jump results with improved signal-to-noise ratio at fast times. The combination of mixing and T-jump results yields an unprecedented view of lactate dehydrogenase enzymology, confirming the timescale of substrate-induced conformational change and presence of multiple reaction pathways.

Conformational Heterogeneity in the Michaelis Complex of Lactate Dehydrogenase: An Analysis of Vibrational Spectroscopy Using Markov and Hidden Markov Models

Lactate dehydrogenase (LDH) catalyzes the interconversion of pyruvate and lactate. Recent isotope-edited IR spectroscopy suggests that conformational heterogeneity exists within the Michaelis complex of LDH, and this heterogeneity affects the propensity toward the on-enzyme chemical step for each Michaelis substate. By combining molecular dynamics simulations with Markov and hidden Markov models, we obtained a detailed kinetic network of the substates of the Michaelis complex of LDH. The ensemble-average electric fields exerted onto the vibrational probe were calculated to provide a direct comparison with the vibrational spectroscopy. Structural features of the Michaelis substates were also analyzed on atomistic scales. Our work not only clearly demonstrates the conformational heterogeneity in the Michaelis complex of LDH and its coupling to the reactivities of the substates, but it also suggests a methodology to simultaneously resolve kinetics and structures on atomistic scales, which can be directly compared with the vibrational spectroscopy.

Active-Loop Dynamics within the Michaelis Complex of Lactate Dehydrogenase from Bacillus stearothermophilus

Laser-induced temperature-jump relaxation spectroscopy was used to study the active site mobile-loop dynamics found in the binding of the NADH nucleotide cofactor and oxamate substrate mimic to lactate dehydrogenase in Bacillus stearothermophilus thermophilic bacteria (bsLDH). The kinetic data can be best described by a model in which NADH can bind only to the open-loop apoenzyme, oxamate can bind only to the bsLDH·NADH binary complex in the open-loop conformation, and oxamate binding is followed by closing of the active site loop preventing oxamate unbinding. The open and closed states of the loop are in dynamic equilibrium and interconvert on the submillisecond time scale. This interconversion strongly accelerates with an increase in temperature because of significant enthalpy barriers. Binding of NADH to bsLDH results in minor changes of the loop dynamics and does not shift the open-closed equilibrium, but binding of the oxamate substrate mimic shifts this equilibrium to the closed state. At high excess oxamate concentrations where all active sites are nearly saturated with the substrate mimic, all active site mobile loops are mainly closed. The observed active-loop dynamics for bsLDH is very similar to that previously observed for pig heart LDH.

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 cis conformation of proline leads to weaker binding of a p53 peptide to MDM2 compared to trans

The cis and trans conformations of the Xaa-Pro (Xaa: any amino acid) peptide bond are thermodynamically stable while other peptide bonds strongly prefer trans. The effect of proline cis-trans isomerization on protein binding has not been thoroughly investigated. In this study, computer simulations were used to calculate the absolute binding affinity for a p53 peptide (residues 17-29) to MDM2 for both cis and trans isomers of the p53 proline in position 27. Results show that the cis isomer of p53(17-29) binds more weakly to MDM2 than the trans isomer, and that this is primarily due to the difference in the free energy cost associated with the loss of conformational entropy of p53(17-29) when it binds to MDM2. The population of cis p53(17-29) was estimated to be 0.8% of the total population in the bound state. The stronger binding of trans p53(17-29) to MDM2 compared to cis may leave a minimal level of p53 available to respond to cellular stress. This study demonstrates that it is feasible to estimate the absolute binding affinity for an intrinsically disordered protein fragment binding to an ordered protein that are in good agreement with experimental results.

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.

Nuclear lactate dehydrogenase modulates histone modification in human hepatocytes

It is becoming increasingly apparent that the nucleus harbors metabolic enzymes that affect genetic transforming events. Here, we describe a nuclear isoform of lactate dehydrogenase (nLDH) and its ability to orchestrate histone deacetylation by controlling the availability of nicotinamide adenine dinucleotide (NAD(+)), a key ingredient of the sirtuin-1 (SIRT1) deacetylase system. There was an increase in the expression of nLDH concomitant with the presence of hydrogen peroxide (H2O2) in the culture medium. Under oxidative stress, the NAD(+) generated by nLDH resulted in the enhanced deacetylation of histones compared to the control hepatocytes despite no discernable change in the levels of SIRT1. There appeared to be an intimate association between nLDH and SIRT1 as these two enzymes co-immunoprecipitated. The ability of nLDH to regulate epigenetic modifications by manipulating NAD(+) reveals an intricate link between metabolism and the processing of genetic information.

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.

Changes in protein architecture and subpicosecond protein dynamics impact the reaction catalyzed by lactate dehydrogenase

We have previously established the importance of a promoting vibration, a subpicosecond protein motion that propagates through a specific axis of residues, in the reaction coordinate of lactate dehydrogenase (LDH). To test the effect that perturbation of this motion would have on the enzymatic reaction, we employ transition path sampling to obtain transition path ensembles for four independent LDH enzymatic systems: the wild type enzyme, a version of the enzyme expressing heavy isotopic substitution, and two enzymes with mutations in the promoting vibration axis. We show that even slight changes in the promoting vibration of LDH result in dramatic changes in enzymatic chemistry. In the "heavy" version of the enzyme, we find that the dampening of the subpicosecond dynamics from heavy isotopic substitution leads to a drastic increase in the time of barrier crossing. Furthermore, we see that mutation of the promoting vibration axis causes a decrease in the variability of transition paths available to the enzymatic reaction. The combined results reveal the importance of the protein architecture of LDH in enzymatic catalysis by establishing how the promoting vibration is finely tuned to facilitate chemistry.

Large scale dynamics of the Michaelis complex in Bacillus stearothermophilus lactate dehydrogenase revealed by a single-tryptophan mutant study

Large scale dynamics within the Michaelis complex mimic of Bacillus stearothermophilus thermophilic lactate dehydrogenase, bsLDH·NADH·oxamate, were studied with site specific resolution by laser-induced temperature jump relaxation spectroscopy with a time resolution of 20 ns. NADH emission and Trp emission from the wild type and a series of single-tryptophan bsLDH mutants, with the tryptophan positions different distances from the active site, were used as reporters of evolving structure in response to the rapid change in temperature. Several distinct dynamical events were observed on the millisecond to microsecond time scale involving motion of atoms spread over the protein, some occurring concomitantly or nearly concomitantly with structural changes at the active site. This suggests that a large portion of the protein-substrate complex moves in a rather concerted fashion to bring about catalysis. The catalytically important surface loop undergoes two distinct movements, both needed for a competent enzyme. Our results also suggest that what is called "loop motion" is not just localized to the loop and active site residues. Rather, it involves the motion of atoms spread over the protein, even some quite distal from the active site. How these results bear on the catalytic mechanism of bsLDH is discussed.

Homology model and targeted mutagenesis identify critical residues for arachidonic acid inhibition of Kv4 channels

Polyunsaturated fatty acids such as arachidonic acid (AA) exhibit inhibitory modulation of Kv4 potassium channels. Molecular docking approaches using a Kv4.2 homology model predicted a membrane-embedded binding pocket for AA comprised of the S4-S5 linker on one subunit and several hydrophobic residues within S3, S5 and S6 from an adjacent subunit. The pocket is conserved among Kv4 channels. We tested the hypothesis that modulatory effects of AA on Kv4.2/KChIP channels require access to this site. Targeted mutation of a polar residue (K318) and a nonpolar residue (G314) within the S4-S5 linker as well as a nonpolar residue in S3 (V261) significantly impaired the effects of AA on K (+) currents in Xenopus oocytes. These residues may be important in stabilizing (K318) or regulating access to (V261, G314) the negatively charged carboxylate moiety on the fatty acid. Structural specificity was supported by the lack of disruption of AA effects observed with mutations at residues located near, but not within the predicted binding pocket. Furthermore, we found that the crystal structure of the related Kv1.2/2.1 chimera lacks the structural features present in the proposed AA docking site of Kv4.2 and the Kv1.2/2.1 K (+) currents were unaffected by AA. We simulated the mutagenic substitutions in our Kv4.2 model to demonstrate how specific mutations may disrupt the putative AA binding pocket. We conclude that AA inhibits Kv4 channel currents and facilitates current decay by binding within a hydrophobic pocket in the channel in which K318 within the S4-S5 linker is a critical residue for AA interaction.

Fragment growing and linking lead to novel nanomolar lactate dehydrogenase inhibitors

Lactate dehydrogenase A (LDH-A) catalyzes the interconversion of lactate and pyruvate in the glycolysis pathway. Cancer cells rely heavily on glycolysis instead of oxidative phosphorylation to generate ATP, a phenomenon known as the Warburg effect. The inhibition of LDH-A by small molecules is therefore of interest for potential cancer treatments. We describe the identification and optimization of LDH-A inhibitors by fragment-based drug discovery. We applied ligand based NMR screening to identify low affinity fragments binding to LDH-A. The dissociation constants (K(d)) and enzyme inhibition (IC(50)) of fragment hits were measured by surface plasmon resonance (SPR) and enzyme assays, respectively. The binding modes of selected fragments were investigated by X-ray crystallography. Fragment growing and linking, followed by chemical optimization, resulted in nanomolar LDH-A inhibitors that demonstrated stoichiometric binding to LDH-A. Selected molecules inhibited lactate production in cells, suggesting target-specific inhibition in cancer cell lines.

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.

Effect of osmolytes on protein dynamics in the lactate dehydrogenase-catalyzed reaction

Laser-induced temperature jump relaxation spectroscopy was used to probe the effect of osmolytes on the microscopic rate constants of the lactate dehydrogenase-catalyzed reaction. NADH fluorescence and absorption relaxation kinetics were measured for the lactate dehydrogenase (LDH) reaction system in the presence of varying amounts of trimethylamine N-oxide (TMAO), a protein-stabilizing osmolyte, or urea, a protein-destabilizing osmolyte. Trimethylamine N-oxide (TMAO) at a concentration of 1 M strongly increases the rate of hydride transfer, nearly nullifies its activation energy, and also slightly increases the enthalpy of hydride transfer. In 1 M urea, the hydride transfer enthalpy is almost nullified, but the activation energy of the step is not affected significantly. TMAO increases the preference of the closed conformation of the active site loop in the LDH·NAD(+)·lactate complex; urea decreases it. The loop opening rate in the LDH·NADH·pyruvate complex changes its temperature dependence to inverse Arrhenius with TMAO. In this complex, urea accelerates the loop motion, without changing the loop opening enthalpy. A strong, non-Arrhenius decrease in the pyruvate binding rate in the presence of TMAO offers a decrease in the fraction of the open loop, pyruvate binding competent form at higher temperatures. The pyruvate off rate is not affected by urea but decreases with TMAO. Thus, the osmolytes strongly affect the rates and thermodynamics of specific events along the LDH-catalyzed reaction: binding of substrates, loop closure, and the chemical event. Qualitatively, these results can be understood as an osmolyte-induced change in the energy landscape of the protein complexes, shifting the conformational nature of functional substates within the protein ensemble.

Loop-tryptophan human purine nucleoside phosphorylase reveals submillisecond protein dynamics

Human PNP is a homotrimer containing three tryptophan residues at positions 16, 94, and 178, all remote from the catalytic site. The catalytic sites of PNP are located near the subunit-subunit interfaces where F159 is a catalytic site residue donated from an adjacent subunit. F159 covers the top (beta) surface of the ribosyl group at the catalytic site. QM/MM calculations of human PNP have shown that F159 is the center of the most mobile region of the protein providing access to the substrate in the active site. F159 is also the key residue in a cluster of hydrophobic residues that shield catalytic site ligands from bulk solvent. Trp-free human PNP (Leuko-PNP) was previously engineered by replacing the three Trp residues of native PNP with Tyr. From this active construct, a single Trp residue was placed in the catalytic site loop (F159W-Leuko-PNP) as a reporter group for the ribosyl region of the catalytic site. The F159W-Leuko-PNP fluorescence is red shifted compared to native PNP, suggesting a solvent-exposed Trp residue. Upon ligand binding (hypoxanthine), the 3-fold fluorescence quench confirms conformational packing of the catalytic site pocket hydrophobic cluster. F159W-Leuko-PNP has an on-enzyme thermodynamic equilibrium constant (Keq) near unity in the temperature range between 20 and 30 degrees C and nonzero enthalpic components, making it suitable for laser-induced T-jump analyses. T-jump relaxation kinetics of F159W-Leuko-PNP in equilibrium with substrates and/or products indicate the conformational equilibria of at least two ternary complex intermediates in the nano- to millisecond time scale (1000-10000 s-1) that equilibrate prior to the slower chemical step (approximately 200 s-1). F159W-Leuko-PNP provides a novel protein platform to investigate the protein conformational dynamics occurring prior to transition state formation.

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.

Tryptophan-free human PNP reveals catalytic site interactions

Human purine nucleoside phosphorylase (PNP) is a homotrimer, containing three nonconserved tryptophan residues at positions 16, 94, and 178, all remote from the catalytic site. The Trp residues were replaced with Tyr to produce Trp-free PNP (Leuko-PNP). Leuko-PNP showed near-normal kinetic properties. It was used (1) to determine the tautomeric form of guanine that produces strong fluorescence when bound to PNP, (2) for thermodynamic binding analysis of binary and ternary complexes with substrates, (3) in temperature-jump perturbation of complexes for evidence of multiple conformational complexes, and (4) to establish the ionization state of a catalytic site tyrosine involved in phosphate nucleophile activation. The (13)C NMR spectrum of guanine bound to Leuko-PNP, its fluorescent properties, and molecular orbital electronic transition analysis establish that its fluorescence originates from the lowest singlet excited state of the N1H, 6-keto, N7H guanine tautomer. Binding of guanine and phosphate to PNP and Leuko-PNP are random, with decreased affinity for formation of ternary complexes. Pre-steady-state kinetics and temperature-jump studies indicate that the ternary complex (enzyme-substrate-phosphate) forms in single binding steps without kinetically significant protein conformational changes as monitored by guanine fluorescence. Spectral changes of Leuko-PNP upon phosphate binding establish that the hydroxyl of Tyr88 is not ionized to the phenolate anion when phosphate is bound. A loop region (residues 243-266) near the purine base becomes highly ordered upon substrate/inhibitor binding. A single Trp residue was introduced into the catalytic loop of Leuko-PNP (Y249W-Leuko-PNP) to determine effects on catalysis and to introduce a fluorescence catalytic site probe. Although Y249W-Leuko-PNP is highly fluorescent and catalytically active, substrate binding did not perturb the fluorescence. Thermodynamic boxes, constructed to characterize the binding of phosphate, guanine, and hypoxanthine to native, Leuko-, and Y249W-Leuko-PNPs, establish that Leuko-PNP provides a versatile protein scaffold for introduction of specific Trp catalytic site probes.

Lactate dehydrogenase undergoes a substantial structural change to bind its substrate

Employing temperature-jump relaxation spectroscopy, we investigate the kinetics and thermodynamics of the formation of a very early ternary binding intermediate formed when lactate dehydrogenase (LDH) binds a substrate mimic on its way to forming the productive LDH/NADH.substrate Michaelis complex. Temperature-jump scans show two distinct submillisecond processes are involved in the formation of this ternary binding intermediate, called the encounter complex here. The on-rate of the formation of the encounter complex from LDH/NADH with oxamate (a substrate mimic) is determined as a function of temperature and in the presence of small concentrations of a protein destabilizer (urea) and protein stabilizer (TMAO). It shows a strong temperature dependence with inverse Arrhenius behavior and a temperature-dependent enthalpy (heat capacity of 610 +/- 84 cal/Mol K), is slowed in the presence of TMAO and speeded up in the presence of urea. These results suggest that LDH/NADH occupies a range of conformations, some competent to bind substrate (open structure; a minority population) and others noncompetent (closed), in fast equilibrium with each other in accord with a select fit model of binding. From the thermodynamic results, the two species differ in the rearrangement of low energy hydrogen bonds as would arise from changes in internal hydrogen bonding and/or increases in the solvation of the protein structure. The binding-competent species can bind ligand at or very near diffusion-limited speeds, suggesting that the binding pocket is substantially exposed to solvent in these species. This would be in contrast to the putative closed structure where the binding pocket resides deep within the protein interior.

Ligand binding and protein dynamics in lactate dehydrogenase

Recent experimental studies suggest that lactate dehydrogenase (LDH) binds its substrate via the formation of a LDH/NADH.substrate encounter complex through a select-fit mechanism, whereby only a minority population of LDH/NADH is binding-competent. In this study, we perform molecular dynamics calculations to explore the variations in structure accessible to the binary complex with a focus on identifying structures that seem likely to be binding-competent and which are in accord with the known experimental characterization of forming binding-competent species. We find that LDH/NADH samples quite a range of protein conformations within our 2.148 ns calculations, some of which yield quite facile access of solvent to the active site. The results suggest that the mobile loop of LDH is perhaps just partially open in these conformations and that multiple open conformations, yielding multiple binding pathways, are likely. These open conformations do not require large-scale unfolding/melting of the binary complex. Rather, open versus closed conformations are due to subtle protein and water rearrangements. Nevertheless, the large heat capacity change observed between binding-competent and binding-incompetent can be explained by changes in solvation and an internal rearrangement of hydrogen bonds. We speculate that such a strategy for binding may be necessary to get a ligand efficiently to a binding pocket that is located fairly deep within the protein's interior.

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.