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

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.

Murburn precepts for lactic-acidosis, Cori cycle, and Warburg effect: Interactive dynamics of dehydrogenases, protons, and oxygen

It is unresolved why lactate is transported to the liver for further utilization within the physiological purview of Cori cycle, when muscles have more lactate dehydrogenase (LDH) than liver. We point out that the answer lies in thermodynamics/equilibriums. While the utilization of NADH for the reduction of pyruvate to lactate can be mediated via the classical mechanism, the oxidation of lactate (with/without the uphill reduction of NAD⁺ ) necessitates alternative physiological approaches. The latter pathway occurs via interactive equilibriums involving the enzyme, protons and oxygen or diffusible reactive oxygen species (DROS). Since liver has high DROS, the murburn activity at LDH would enable the cellular system to tide over the unfavorable energy barriers of the forward reaction (~476 kJ/mol; earlier miscalculated as ~26 kJ/mole). Further, the new mechanism does not necessitate any "smart decision-making" or sophisticated control by/of proteins. The DROS-based murburn theory explains the invariant active-site structure of LDH isozymes and their multimeric nature. The theoretical insights, in silico evidence and analyses of literature herein also enrich our understanding of the underpinnings of "lactic acidosis" (lowering of physiological pH accompanied by lactate production), Warburg effect (increased lactate production at high pO₂ by cancer cells) and approach for cancer therapy.

Structure, Function, and Thermodynamics of Lactate Dehydrogenases from Humans and the Malaria Parasite P. falciparum

Temperature adaptation is ubiquitous among all living organisms, yet the molecular basis for this process remains poorly understood. It can be assumed that for parasite-host systems, the same enzymes found in both organisms respond to the same selection factor (human body temperature) with similar structural changes. Herein, we report the existence of a reversible temperature-dependent structural transition for the glycolytic enzyme lactate dehydrogenase (LDH) from the malaria parasite Plasmodium falciparum (pfLDH) and human heart (hhLDH) occurring in the temperature range of human fever. This transition is observed for LDHs from psychrophiles, mesophiles, and moderate thermophiles in their operating temperature range. Thermodynamic analysis reveals unique thermodynamic signatures of the LDH-substrate complexes defining a specific temperature range to which human LDH is adapted and parasite LDH is not, despite their common mesophilic nature. The results of spectroscopic analysis combined with the available crystallographic data reveal the existence of an active center within pfLDH that imparts psychrophilic structural properties to the enzyme. This center consists of two pockets, one formed by the five amino acids (5AA insert) within the substrate specificity loop and the other by the active site, that mutually regulate one another in response to temperature and induce structural and functional changes in the Michaelis complex. Our findings pave the way toward a new strategy for malaria treatments and drug design using therapeutic agents that inactivate malarial LDH selectively at a specific temperature range of the cyclic malaria paroxysm.

Lactate dehydrogenases amplify reactive oxygen species in cancer cells in response to oxidative stimuli

Previous studies demonstrated that superoxide could initiate and amplify LDH-catalyzed hydrogen peroxide production in aqueous phase, but its physiological relevance is unknown. Here we showed that LDHA and LDHB both exhibited hydrogen peroxide-producing activity, which was significantly enhanced by the superoxide generated from the isolated mitochondria from HeLa cells and patients’ cholangiocarcinoma specimen. After LDHA or LDHB were knocked out, hydrogen peroxide produced by Hela or 4T1 cancer cells were significantly reduced. Re-expression of LDHA in LDHA-knockout HeLa cells partially restored hydrogen peroxide production. In HeLa and 4T1 cells, LDHA or LDHB knockout or LDH inhibitor FX11 significantly decreased ROS induction by modulators of the mitochondrial electron transfer chain (antimycin, oligomycin, rotenone), hypoxia, and pharmacological ROS inducers piperlogumine (PL) and phenethyl isothiocyanate (PEITC). Moreover, the tumors formed by LDHA or LDHB knockout HeLa or 4T1 cells exhibited a significantly less oxidative state than those formed by control cells. Collectively, we provide a mechanistic understanding of a link between LDH and cellular hydrogen peroxide production or oxidative stress in cancer cells in vitro and in vivo.

Effect of Protein Isotope Labeling on the Catalytic Mechanism of Lactate Dehydrogenase

We investigate how isotopic labeling of the enzyme lactate dehydrogenase (LDH) affects its function. LDH is of special interest because there exists a line of residues spanning the protein that are involved in the transition state (TS) of the chemical reaction coordinate (so-called promoting vibration). Hence, studies have been carried out on this protein (as well as others) using labeled protein (so-called heavy protein) along with measurements of single turnover k_cat yielding a KIE (=k_cat^light/k_cat^heavy) aimed at understanding the effect of labeling generally and more specifically this line of residues. Here, it is shown that ¹³C, ¹⁵N, and ²H atom labeling of hhLDH (human heart) affects its internal structure which in turn affects its dynamics and catalytic mechanism. Spectral studies employing advanced FTIR difference spectroscopy show that the height of the electronic potential surface of the TS is lowered (probably by ground state destabilization) by labeling. Moreover, laser-induced T-jump relaxation kinetic spectroscopy shows that the microsecond to millisecond nuclear motions internal to the protein are affected by labeling. While the effects are small, they are sufficient to contribute to the observed KIE values as well or even more than promoting vibration effects.

Engineering D-Lactate Dehydrogenase from Pediococcus acidilactici for Improved Activity on 2-Hydroxy Acids with Bulky C3 Functional Group

Engineering D-lactic acid dehydrogenases for higher activity on various 2-oxo acids is important for the synthesis of 2-hydroxy acids that can be utilized in a wide range of industrial fields including the production of biopolymers, pharmaceuticals, and cosmetic compounds. Although there are many D-lactate dehydrogenases (D-LDH) available from a diverse range of sources, there is a lack of biocatalysts with high activities for 2-oxo acids with large functional group at C₃. In this study, the D-LDH from Pediococcus acidilactici was rationally designed and further engineered by controlling the intermolecular interactions between substrates and the surrounding residues via analysis of the active site structure of D-LDH. As a result, Y51L mutant with the catalytic efficiency on phenylpyruvate of 2200 s^-1 mM^-1 and Y51F mutant on 2-oxobutryate and 3-methyl-2-oxobutyrate of 37.2 and 23.2 s^-1 mM^-1 were found, which were 138-, 8.5-, and 26-fold increases than the wild type on the substrates, respectively. Structural analysis revealed that the distance and the nature of the interactions between the side chain of residue 51 and the substrate C₃ substituent group significantly affected the kinetic parameters. Bioconversion of phenyllactate as a practical example of production of the 2-hydroxy acids was investigated, and the Y51F mutant presented the highest productivity in in vitro conversion of D-PLA.

Small molecule cores demonstrate non-competitive inhibition of lactate dehydrogenase

Lactate dehydrogenase (LDH) has recently garnered attention as an attractive target for cancer therapies, owing to the enzyme's critical role in cellular metabolism. Current inhibition strategies, employing substrate or cofactor analogues, are insufficiently specific for use as pharmaceutical agents. The possibility of allosteric inhibition of LDH was postulated on the basis of theoretical docking studies of a small molecule inhibitor to LDH. The present study examined structural analogues of this proposed inhibitor to gauge its potency and attempt to elucidate the molecular mechanism of action. These analogues display encouraging in vitro inhibition of porcine heart LDH, including micromolar K_i values and a maximum inhibition of up to 50% in the steady state. Furthermore, Michaelis-Menten kinetics and fluorescence data both suggest the simple, acetaminophen derivatives are non-competitive in binding to the enzyme. Kinetic comparisons of a panel of increasingly decorated structural analogues imply that the binding is specific, and the small molecule core provides a privileged scaffold for further pharmaceutical development of a novel, allosteric drug.

Trehalose Mediated Inhibition of Lactate Dehydrogenase from Rabbit Muscle. The Application of Kramers' Theory in Enzyme Catalysis

Lactate dehydrogenase (LDH) catalyzes the reduction of pyruvate to lactate by using NADH. LDH kinetics has been proposed to be dependent on the dynamics of a loop over the active site. Kramers' theory has been useful in the study of enzyme catalysis dependent on large structural dynamics. In this work, LDH kinetics was studied in the presence of trehalose and at different temperatures. In the absence of trehalose, temperature increase raised exponentially the LDH V_max and revealed a sigmoid transition of K_m toward a low-affinity state similar to protein unfolding. Notably, LDH V_max diminished when in the presence of trehalose, while pyruvate affinity increased and the temperature-mediated binding site transition was hindered. The effect of trehalose on k_cat was viscosity dependent as described by Kramers' theory since V_max correlated inversely with the viscosity of the medium. As a result, activation energy ( E_a) for pyruvate reduction was dramatically increased by trehalose presence. This work provides experimental evidence that the dynamics of a structural component in LDH is essential for catalysis, i.e., the closing of the loop on the active site. While the trehalose mediated-increased of pyruvate affinity is proposed to be due to the compaction and/or increase of structural order at the binding site.

Mechanistic Analysis of Fluorescence Quenching of Reduced Nicotinamide Adenine Dinucleotide by Oxamate in Lactate Dehydrogenase Ternary Complexes

Fluorescence of Reduced Nicotinamide Adenine Dinucleotide (NADH) is extensively employed in studies of oxidoreductases. A substantial amount of static and kinetic work has focused on the binding of pyruvate or substrate mimic oxamate to the binary complex of lactate dehydrogenase (LDH)-NADH where substantial fluorescence quenching is typically observed. However, the quenching mechanism is not well understood limiting structural interpretation. Based on time-dependent density functional theory (TDDFT) computations with cam-B3LYP functional in conjunction with the analysis of previous experimental results, we propose that bound oxamate acts as an electron acceptor in the quenching of fluorescence of NADH in the ternary complex, where a charge transfer (CT) state characterized by excitation from the highest occupied molecular orbital (HOMO) of the nicotinamide moiety of NADH to the lowest unoccupied molecular orbital (LUMO) of oxamate exists close to the locally excited (LE) state involving only the nicotinamide moiety. Efficient quenching in the encounter complex like in pig heart LDH requires that oxamate forms a salt bridge with Arg-171 and hydrogen bonds with His-195, Thr-246 and Asn-140. Further structural rearrangement and loop closure, which also brings about another hydrogen bond between oxamate and Arg-109, will increase the rate of fluorescence quenching as well.

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.

Triple Isotope Effects Support Concerted Hydride and Proton Transfer and Promoting Vibrations in Human Heart Lactate Dehydrogenase

Transition path sampling simulations have proposed that human heart lactate dehydrogenase (LDH) employs protein promoting vibrations (PPVs) on the femtosecond (fs) to picosecond (ps) time scale to promote crossing of the chemical barrier. This chemical barrier involves both hydride and proton transfers to pyruvate to form l-lactate, using reduced nicotinamide adenine dinucleotide (NADH) as the cofactor. Here we report experimental evidence from three types of isotope effect experiments that support coupling of the promoting vibrations to barrier crossing and the coincidence of hydride and proton transfer. We prepared the native (light) LDH and a heavy LDH labeled with ¹³C, ¹⁵N, and nonexchangeable ²H (D) to perturb the predicted PPVs. Heavy LDH has slowed chemistry in single turnover experiments, supporting a contribution of PPVs to transition state formation. Both the [4-²H]NADH (NADD) kinetic isotope effect and the D₂O solvent isotope effect were increased in dual-label experiments combining both NADD and D₂O, a pattern maintained with both light and heavy LDHs. These isotope effects support concerted hydride and proton transfer for both light and heavy LDHs. Although the transition state barrier-crossing probability is reduced in heavy LDH, the concerted mechanism of the hydride-proton transfer reaction is not altered. This study takes advantage of triple isotope effects to resolve the chemical mechanism of LDH and establish the coupling of fs-ps protein dynamics to barrier crossing.

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.

Direct evidence of catalytic heterogeneity in lactate dehydrogenase by temperature jump infrared spectroscopy

Protein conformational heterogeneity and dynamics are known to play an important role in enzyme catalysis, but their influence has been difficult to observe directly. We have studied the effects of heterogeneity in the catalytic reaction of pig heart lactate dehydrogenase using isotope edited infrared spectroscopy, laser-induced temperature jump relaxation, and kinetic modeling. The isotope edited infrared spectrum reveals the presence of multiple reactive conformations of pyruvate bound to the enzyme, with three major reactive populations having substrate C2 carbonyl stretches at 1686, 1679, and 1674 cm(-1), respectively. The temperature jump relaxation measurements and kinetic modeling indicate that these substates form a heterogeneous branched reaction pathway, and each substate catalyzes the conversion of pyruvate to lactate with a different rate. Furthermore, the rate of hydride transfer is inversely correlated with the frequency of the C2 carbonyl stretch (the rate increases as the frequency decreases), consistent with the relationship between the frequency of this mode and the polarization of the bond, which determines its reactivity toward hydride transfer. The enzyme does not appear to be optimized to use the fastest pathway preferentially but rather accesses multiple pathways in a search process that often selects slower ones. These results provide further support for a dynamic view of enzyme catalysis where the role of the enzyme is not just to bring reactants together but also to guide the conformational search for chemically competent interactions.

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.

Common enzymological experiments allow free energy profile determination

The determination of a complete set of rate constants [free energy profiles (FEPs)] for a complex kinetic mechanism is challenging. Enzymologists have devised a variety of informative steady-state kinetic experiments (e.g., Michaelis-Menten kinetics, viscosity dependence of kinetic parameters, kinetic isotope effects, etc.) that each provide distinct information regarding a particular kinetic system. A simple method for combining steady-state experiments in a single analysis is presented here, which allows microscopic rate constants and intrinsic kinetic isotope effects to be determined. It is first shown that Michaelis-Menten kinetic parameters (kcat and Km values), kinetic isotope efffets, solvent viscosity effects, and intermediate partitioning measurements are sufficient to define the rate constants for a reversible uni-uni mechanism with an intermediate, EZ, between the ES and EP complexes. Global optimization provides the framework for combining the independent experimental measurements, and the search for rate constants is performed using algorithms implemented in the biochemical software COPASI. This method is applied to the determination of FEPs for both alanine racemase and triosephosphate isomerase. The FEPs obtained from global optimization agree with those in the literature, with important exceptions. The method opens the door to routine and large-scale determination of FEPs for enzymes.

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.

Investigation of catalytic loop structure, dynamics, and function relationship of Yersinia protein tyrosine phosphatase by temperature-jump relaxation spectroscopy and X-ray structural determination

Yersinia protein tyrosine phosphatase (YopH) is the most efficient enzyme among all PTPases. YopH is hyperactive compared to human PTPases, interfering with mammalian cellular pathways to achieve the pathogenicity of Yersinia. Two properties related to the catalytic loop structure differences have been proposed to affect its dynamics and enzyme efficiency. One is the ability of the loop to form stabilizing interactions to bound ligand after loop closure, which has long been recognized. In addition, the loop flexibility/mobility was suggested in a previous study to be a factor as well, based on the observation that incremental changes in PTPase loop structure by single point mutations to alanine often induce incremental changes in enzyme catalytic efficiency. In this study, the temperature jump relaxation spectroscopy (T-jump) has been used to discern the subtle changes of the loop dynamics due to point loop mutations. As expected, our results suggest a correlation between loop dynamics and the size of the residue on the catalytic loop. The stabilization of the enzyme-ligand complex is often enthalpy driven, achieved by formation of additional favorable hydrogen bonding/ionic interactions after loop closure. Interestingly, our T-jump and X-ray crystallography studies on YopH suggest that the elimination of some ligand-protein interactions by mutation does not necessarily destabilize the ligand-enzyme complex after loop closure, since the increased entropy in the forms of more mobile protein residues may be sufficient to compensate the free energy loss due to lost interactions and may even lead to enhanced efficiency of the enzyme catalysis. How these competing loop properties may affect loop dynamics and enzyme function are discussed.

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.

Oxamic acid analogues as LDH-C4-specific competitive inhibitors

We performed kinetic studies to determine whether oxamate analogues are selective inhibitors of LDH-C4, owing to their potential usefulness in fertility control and treatment of some cancers. These substances were shown to be competitive inhibitors of LDH isozymes and are able to discriminate among subtle differences that differentiate the active sites of LDH-A4, LDH-B4 and LDH-C4. N-Ethyl oxamate was the most potent inhibitor showing the highest affinity for LDH-C4. However, N-propyl oxamate was the most selective inhibitor showing a high degree of selectivity towards LDH-C4. Non-polar four carbon atoms chains, linear or branched, dramatically diminished the affinity and selectivity towards LDH-C4. N-Propyl oxamate significantly reduced ATP levels, capacitation and mouse sperm motility, in line with results shown by others, suggesting that LDH-C4 plays an essential role in mouse fertility.

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.

Can hydridic-to-protonic hydrogen bonds catalyze hydride transfers in biological systems?

Catalysis of hydride transfer by hydridic-to-protonic hydrogen (HHH) bonding in α-hydroxy carbonyl isomerization reactions was examined computationally in the lithium salts of 7-substituted endo-3-hydroxybicyclo[2.2.1]hept-5-en-2-ones. The barrier for intramolecular hydride transfer in the parent system was calculated to be 17.2 kcal/mol. Traditional proton donors, such as OH and NH(3)(+), stabilized the metal cation-bridged transition state by 1.4 and 3.3 kcal/mol, respectively. Moreover, among the conformers of the OH systems, the one in which the proton donor is able to interact with the migrating hydride (H(m)) has an activation barrier lower by 1.3 and 1.7 kcal/mol than the other possible OH conformers. By contrast, the presence of an electronegative group such as F, which disfavors the migration electronically by opposing development of hydridic charge, destabilizes the hydride migration by 1.5 kcal/mol relative to the epimeric exo system. In both ground and transition states the H(m)···H distance decreased with increasing acidity of the proton donor, reaching a minimum of 1.58 Å at the transition state for NH(3)(+). Both Mulliken and NPA charges show enhancement of negative character of the migrating hydride in the cases in which HHH bonding is possible.

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.

Mechanistic role of movement and strain sensitivity in muscle contraction

Tension generation can be studied by applying step perturbations to contracting muscle fibers and subdividing the mechanical response into exponential phases. The de novo tension-generating isomerization is associated with one of these phases. Earlier work has shown that a temperature jump perturbs the equilibrium constant directly to increase tension. Here, we show that a length jump functions quite differently. A step release (relative movement of thick and thin filaments) appears to release a steric constraint on an ensemble of noncompetent postphosphate release actomyosin cross-bridges, enabling them to generate tension, a concentration jump in effect. Structural studies [Taylor KA, et al. (1999) Tomographic 3D reconstruction of quick-frozen, Ca(2+)-activated contracting insect flight muscle. Cell 99:421-431] that map to these kinetics indicate that both catalytic and lever arm domains of noncompetent myosin heads change angle on actin, whereas lever arm movement alone mediates the power stroke. Together, these kinetic and structural observations show a 13-nm overall interaction distance of myosin with actin, including a final 4- to 6-nm power stroke when the catalytic domain is fixed on actin. Raising fiber temperature with both perturbation techniques accelerates the forward, but slows the reverse rate constant of tension generation, kinetics akin to the unfolding/folding of small proteins. Decreasing strain, however, causes both forward and reverse rate constants to increase. Despite these changes in rate, the equilibrium constant is strain-insensitive. Activation enthalpy and entropy data show this invariance to be the result of enthalpy-entropy compensation. Reaction amplitudes confirm a strain-invariant equilibrium constant and thus a strain-insensitive ratio of pretension- to tension-generating states as work is done.