Design and synthesis of new enzymes based on the lactate dehydrogenase framework

LDHA: The Obstacle to T cell responses against tumor

Immunotherapy has become a successful therapeutic strategy in certain solid tumors and hematological malignancies. However, this efficacy of immunotherapy is impeded by limited success rates. Cellular metabolic reprogramming determines the functionality and viability in both cancer cells and immune cells. Extensive research has unraveled that the limited success of immunotherapy is related to immune evasive metabolic reprogramming in tumor cells and immune cells. As an enzyme that catalyzes the final step of glycolysis, lactate dehydrogenase A (LDHA) has become a major focus of research. Here, we have addressed the structure, localization, and biological features of LDHA. Furthermore, we have discussed the various aspects of epigenetic regulation of LDHA expression, such as histone modification, DNA methylation, N6-methyladenosine (m6A) RNA methylation, and transcriptional control by noncoding RNA. With a focus on the extrinsic (tumor cells) and intrinsic (T cells) functions of LDHA in T-cell responses against tumors, in this article, we have reviewed the current status of LDHA inhibitors and their combination with T cell-mediated immunotherapies and postulated different strategies for future therapeutic regimens.

Conformational Change and Activity Enhancement of Rabbit Muscle Lactate Dehydrogenase Induced by Polyethyleneimine

For a better understanding on the interaction between polyethyleneimine (PEI) and proteins, spectroscopic studies including UV-vis absorption, resonance Rayleigh scattering, fluorescence, and circular dichroism were conducted to reveal the conformational change of rabbit muscle lactate dehydrogenase (rmLDH) and related to the bioactivity of the enzyme. Regardless of the electrostatic repulsion, PEI could bind on the surface of rmLDH, a basic protein, via hydrogen binding of the dense amine groups and hydrophobic interaction of methyl groups. The competitive binding by PEI led to a reduction of the binding efficiency of rmLDH toward β-nicotinamide adenine dinucleotide, the coenzyme, and sodium pyruvate, the substrate. However, the complex formation with PEI induced a less ordered conformation and an enhanced surface hydrophobicity of rmLDH, facilitating the turnover of the enzyme and generally resulting in an increased activity. PEI of higher molecular weight was more efficient to induce alteration in the conformation and catalytic activity of the enzyme.

Living on the Edge: Physiological and Kinetic Trade-Offs Shape Thermal Tolerance in Intertidal Crabs From Tropical to Sub-Antarctic South America

Temperature is an important abiotic factor that drives the evolution of ectotherms owing to its pervasive effects at all levels of organization. Although a species' thermal tolerance is environmentally driven within a spatial cline, it may be constrained over time due to differential phylogenetic inheritance. At the limits of thermal tolerance, hemolymph oxygen is reduced and lactate formation is increased due to mismatch between oxygen supply and demand; imbalance between enzyme flexibility/stability also impairs the ability to generate energy. Here, we characterized the effects of lower (LL₅₀) and upper (UL₅₀) critical thermal limits on selected descriptors of aerobic and anaerobic metabolism in 12 intertidal crab species distributed from northern Brazil (≈7.8°S) to southern Patagonia (≈53.2°S), considering their phylogeny. We tested for (i) functional trade-offs regarding aerobic and anaerobic metabolism and LDH kinetics in shaping thermal tolerance; (ii) influence of shared ancestry and thermal province on metabolic evolution; and (iii) presence of evolutionary convergences and adaptive peaks in the crab phylogeny. The tropical and subtropical species showed similar systemic and kinetic responses, both differing from the sub-Antarctic crabs. The lower UL₅₀'s of the sub-Antarctic crabs may reflect mismatch between the evolution of aerobic and anaerobic metabolism since these crabs exhibit lower oxygen consumption but higher lactate formation than tropical and subtropical species also at their respective UL₅₀'s. LDH activity increased with temperature increase, while K_m ^Pyr remained fairly constant; catalytic coefficient correlated negatively with thermal niche. Thermal tolerance may rely on a putative evolutionary trade-off between aerobic and anaerobic metabolism regarding energy supply, while temperature compensation of kinetic performance is driven by thermal habitat as revealed by the LDH affinity/efficiency equilibrium. The overall physiological evolution revealed two homoplastic adaptive peaks in the sub-Antarctic crabs with a further shift in the tropical/subtropical clade. The physiological traits at UL₅₀ have evolved in a phylogenetic manner while all others were more plastic. Thus, shared inheritance and thermal environment have driven the crabs' thermal tolerance and metabolic evolution, revealing physiological transformations that have arisen in both colder and warmer climes, especially at higher levels of biological organization and phylogenetic diversity.

Structural and functional regulation of lactate dehydrogenase-A in cancer

Dysregulated metabolism is one of the hallmarks of cancer. Under normal physiological conditions, ATP is primarily generated by oxidative phosphorylation. Cancers commonly undergo a dramatic shift toward glycolysis, despite the presence of oxygen. This phenomenon is known as the Warburg effect, and requires the activity of LDHA. LDHA converts pyruvate to lactate in the final step of glycolysis and is often upregulated in cancer. LDHA inhibitors present a promising therapeutic option, as LDHA blockade leads to apoptosis in cancer cells. Despite this, existing LDHA inhibitors have shown limited clinical efficacy. Here, we review recent progress in LDHA structure, function and regulation as well as strategies to target this critical enzyme.

Crystal structures of Babesia microti lactate dehydrogenase BmLDH reveal a critical role for Arg99 in catalysis

The tick- and transfusion-transmitted human pathogen Babesia microti infects host erythrocytes to cause the pathologic symptoms associated with human babesiosis, an emerging disease with worldwide distribution and potentially fatal clinical outcome. Drugs currently recommended for the treatment of babesiosis are associated with a high failure rate and significant adverse events, highlighting the urgent need for more-effective and safer babesiosis therapies. Unlike other apicomplexan parasites, B. microti lacks a canonical lactate dehydrogenase (LDH) but instead expresses a unique enzyme, B. microti LDH (BmLDH), acquired through evolution by horizontal transfer from a mammalian host. Here, we report the crystal structures of BmLDH in apo state and ternary complex (enzyme-NADH-oxamate) solved at 2.79 and 1.89 Å. Analysis of these structures reveals that upon binding to the coenzyme and substrate, the active pocket of BmLDH undergoes a major conformational change from an opened and disordered to a closed and stabilized state. Biochemical assays using wild-type and mutant B. microti and human LDHs identified Arg99 as a critical residue for the catalytic activity of BmLDH but not its human counterpart. Interestingly, mutation of Arg99 to Ala had no impact on the overall structure and affinity of BmLDH to NADH but dramatically altered the closure of the enzyme's active pocket. Together, these structural and biochemical data highlight significant differences between B. microti and human LDH enzymes and suggest that BmLDH could be a suitable target for the development of selective antibabesial inhibitors.-Yu, L., Shen, Z., Liu, Q., Zhan, X., Luo, X., An, X., Sun, Y., Li, M., Wang, S., Nie, Z., Ao, Y., Zhao, Y., Peng, G., Ben Mamoun, C., He, L., Zhao, J. Crystal structures of Babesia microti lactate dehydrogenase BmLDH reveal a critical role for Arg99 in catalysis.

Roles of Closed- and Open-Loop Conformations in Large-Scale Structural Transitions of l-Lactate Dehydrogenase

The mechanism of l-lactate generation from pyruvate by l-lactate dehydrogenase (LDH) from the rabbit muscle was studied theoretically by the multistructural microiteration (MSM) method combined with the quantum mechanics/molecular mechanics (QM/MM)-ONIOM method, where the MSM method describes the MM environment as a weighted average of multiple different structures that are fully relaxed during geometry optimization or a reaction path calculation for the QM part. The results showed that the substrate binding and product states were stabilized only in the open-loop conformation of LDH and the reaction occurred in the closed-loop conformation. In other words, before and after the chemical reaction, a large-scale structural transition from the open-loop conformation to the closed-loop conformation and vice versa occurred. The closed-loop conformation stabilized the transition state of the reaction. In contrast, the open-loop conformation stabilized the substrate binding and final states. In other words, the closed- to open-loop transition at the substrate binding state urges capture of the substrate molecule, the subsequent open- to closed-loop transition promotes the product generation, and the final closed- to open-loop transition at the final state prevents the reverse reaction going back to the substrate binding state. It is thus suggested that the exchange of stability between the closed- and open-loop conformations at different states promotes the catalytic cycle.

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.

L-2-Hydroxyglutarate production arises from noncanonical enzyme function at acidic pH

The metabolite 2-hydroxyglutarate (2HG) can be produced as either a D(R)- or L(S)- enantiomer, each of which inhibits alpha-ketoglutarate (αKG)-dependent enzymes involved in diverse biologic processes. Oncogenic mutations in isocitrate dehydrogenase produce D-2HG, which causes a pathologic blockade in cell differentiation. On the other hand, oxygen limitation leads to accumulation of L-2HG, which can facilitate physiologic adaptation to hypoxic stress in both normal and malignant cells. Here we demonstrate that purified lactate dehydrogenase (LDH) and malate dehydrogenase (MDH) catalyze stereospecific production of L-2HG via ‘promiscuous’ reduction of the alternative substrate αKG. Acidic pH enhances production of L-2HG by promoting a protonated form of αKG that binds to a key residue in the substrate-binding pocket of LDHA. Acid-enhanced production of L-2HG leads to stabilization of hypoxia-inducible factor 1 alpha (HIF-1α) in normoxia. These findings offer insights into mechanisms whereby microenvironmental factors influence production of metabolites that alter cell fate and function.

Lactate Dehydrogenase C Produces S-2-Hydroxyglutarate in Mouse Testis

Metabolomics is a valuable tool for studying tissue- and organism-specific metabolism. In normal mouse testis, we found 70 μM S-2-hydroxyglutarate (2HG), more than 10-fold greater than in other tissues. S-2HG is a competitive inhibitor of α-ketoglutarate-dependent demethylation enzymes and can alter histone or DNA methylation. To identify the source of testis S-2HG, we fractionated testis extracts and identified the fractions that actively produced S-2HG. Through a combination of ion exchange and size exclusion chromatography, we enriched a single active protein, the lactate dehydrogenase isozyme LDHC, which is primarily expressed in testis. At neutral pH, recombinant mouse LDHC rapidly converted both pyruvate into lactate and α-ketoglutarate into S-2HG, whereas recombinant human LDHC only produced lactate. Rapid S-2HG production by LDHC depends on amino acids 100-102 being Met-Val-Ser, a sequence that occurs only in the rodent protein. Other mammalian LDH can also produce some S-2HG, but at acidic pH. Thus, polymorphisms in the Ldhc gene control testis levels of S-2HG, and thereby epigenetics, across mammals.

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.

Phase Space Bottlenecks in Enzymatic Reactions

The definition of a transition state on an individual reactive trajectory is made via a committor analysis. In the past, the bottleneck definition has often been applied in configuration space. This is an approximation, and in order to expand this definition, we are revisiting an enzyme in which we had identified a fast subpicosecond motion that makes the reaction possible. First we used a time-series analysis method to identify the exact time when this motion initiates donor-acceptor compression. Then we modified the standard committor analysis of transition path sampling to identify events in phase space and found that there is a dividing surface in phase space significantly earlier than the configurationally defined transition-state crossing.

Adaptations of protein structure and function to temperature: there is more than one way to ‘skin a cat’

Sensitivity to temperature helps determine the success of organisms in all habitats, and is caused by the susceptibility of biochemical processes, including enzyme function, to temperature change. A series of studies using two structurally and catalytically related enzymes, A4-lactate dehydrogenase (A4-LDH) and cytosolic malate dehydrogenase (cMDH) have been especially valuable in determining the functional attributes of enzymes most sensitive to temperature, and identifying amino acid substitutions that lead to changes in those attributes. The results of these efforts indicate that ligand binding affinity and catalytic rate are key targets during temperature adaptation: ligand affinity decreases during cold adaptation to allow more rapid catalysis. Structural changes causing these functional shifts often comprise only a single amino acid substitution in an enzyme subunit containing approximately 330 residues; they occur on the surface of the protein in or near regions of the enzyme that move during catalysis, but not in the active site; and they decrease stability in cold-adapted orthologs by altering intra-molecular hydrogen bonding patterns or interactions with the solvent. Despite these structure–function insights, we currently are unable to predict a priori how a particular substitution alters enzyme function in relation to temperature. A predictive ability of this nature might allow a proteome-wide survey of adaptation to temperature and reveal what fraction of the proteome may need to adapt to temperature changes of the order predicted by global warming models. Approaches employing algorithms that calculate changes in protein stability in response to a mutation have the potential to help predict temperature adaptation in enzymes; however, using examples of temperature-adaptive mutations in A4-LDH and cMDH, we find that the algorithms we tested currently lack the sensitivity to detect the small changes in flexibility that are central to enzyme adaptation to temperature.

Screening of novel inhibitors targeting lactate dehydrogenase A via four molecular docking strategies and dynamics simulations

Lactate dehydrogenase A (LDHA) is a metabolic enzyme which catalyzes the interconversion of lactate and pyruvate in the glycolysis pathway, thus playing key roles in aerobic glycolysis. The inhibition of LDHA by small molecules has become an attractive strategy for anticancer therapy in recent years. However, very few LDHA inhibitors have been reported, even though a great deal of effort has directed into identifying LDHA inhibitors using structure-based approaches. Therefore, high-throughput and high-accuracy screening approaches are still urgently needed in order to target LDHA effectively. In the present work, after establishing that our docking strategies performed well using test datasets, we screened 32791 Specs products for their docking scores with the substrate-binding pocket and, separately, the cofactor-binding pocket of LDHA. We subsequently identified 76 hits (i.e., ligands that show low docking scores) for the cofactor-binding pocket and 27 hits for the substrate-binding pocket. Two representative compounds, ZINC20036549 and ZINC19369718, were then chosen for further MD simulation analysis, and we found that these compounds maintained their inhibitory activity during the MD simulations. Meanwhile, we found that ZINC19369718 interacts with a novel binding site close to the active site, and that this interaction may inhibit the catalytic activity of LDHA. Together, these results offer not only a new paradigm for identifying Specs drug-like products for novel therapeutic use but they also provide further opportunity to adopt LDHA inhibition as a strategy for cancer therapy.

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.

Structures of lactate dehydrogenase A (LDHA) in apo, ternary and inhibitor-bound forms

Lactate dehydrogenase (LDH) is an essential metabolic enzyme that catalyzes the interconversion of pyruvate and lactate using NADH/NAD(+) as a co-substrate. Many cancer cells exhibit a glycolytic phenotype known as the Warburg effect, in which elevated LDH levels enhance the conversion of glucose to lactate, making LDH an attractive therapeutic target for oncology. Two known inhibitors of the human muscle LDH isoform, LDHA, designated 1 and 2, were selected, and their IC50 values were determined to be 14.4 ± 3.77 and 2.20 ± 0.15 µM, respectively. The X-ray crystal structures of LDHA in complex with each inhibitor were determined; both inhibitors bind to a site overlapping with the NADH-binding site. Further, an apo LDHA crystal structure solved in a new space group is reported, as well as a complex with both NADH and the substrate analogue oxalate bound in seven of the eight molecules and an oxalate only bound in the eighth molecule in the asymmetric unit. In this latter structure, a kanamycin molecule is located in the inhibitor-binding site, thereby blocking NADH binding. These structures provide insights into LDHA enzyme mechanism and inhibition and a framework for structure-assisted drug design that may contribute to new cancer therapies.

Protein Conformational Landscapes and Catalysis. Influence of Active Site Conformations in the Reaction Catalyzed by L-Lactate Dehydrogenase

In the last decade L-Lactate Dehydrogenase (LDH) has become an extremely useful marker in both clinical diagnosis and in monitoring the course of many human diseases. It has been assumed from the 80s that the full catalytic process of LDH starts with the binding of the cofactor and the substrate followed by the enclosure of the active site by a mobile loop of the protein before the reaction to take place. In this paper we show that the chemical step of the LDH catalyzed reaction can proceed within the open loop conformation, and the different reactivity of the different protein conformations would be in agreement with the broad range of rate constants measured in single molecule spectrometry studies. Starting from a recently solved X-ray diffraction structure that presented an open loop conformation in two of the four chains of the tetramer, QM/MM free energy surfaces have been obtained at different levels of theory. Depending on the level of theory used to describe the electronic structure, the free energy barrier for the transformation of pyruvate into lactate with the open conformation of the protein varies between 12.9 and 16.3 kcal/mol, after quantizing the vibrations and adding the contributions of recrossing and tunneling effects. These values are very close to the experimentally deduced one (14.2 kcal·mol^-1) and ~2 kcal·mol^-1 smaller than the ones obtained with the closed loop conformer. Calculation of primary KIEs and IR spectra in both protein conformations are also consistent with our hypothesis and in agreement with experimental data. Our calculations suggest that the closure of the active site is mainly required for the inverse process; the oxidation of lactate to pyruvate. According to this hypothesis H4 type LDH enzyme molecules, where it has been propose that lactate is transformed into pyruvate, should have a better ability to close the mobile loop than the M4 type LDH molecules.

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.

Human lactate dehydrogenase a inhibitors: a molecular dynamics investigation

Lactate dehydrogenase A (LDHA) is an important enzyme in fermentative glycolysis, generating most energy for cancer cells that rely on anaerobic respiration even under normal oxygen concentrations. This renders LDHA a promising molecular target for the treatment of various cancers. Several efforts have been made recently to develop LDHA inhibitors with nanomolar inhibition and cellular activity, some of which have been studied in complex with the enzyme by X-ray crystallography. In this work, we present a molecular dynamics (MD) study of the binding interactions of selected ligands with human LDHA. Conventional MD simulations demonstrate different binding dynamics of inhibitors with similar binding affinities, whereas steered MD simulations yield discrimination of selected LDHA inhibitors with qualitative correlation between the in silico unbinding difficulty and the experimental binding strength. Further, our results have been used to clarify ambiguities in the binding modes of two well-known LDHA inhibitors.

A hybrid elastic band string algorithm for studies of enzymatic reactions

A common challenge in theoretical biophysics is the identification of a minimum energy path (MEP) for the rearrangement of a group of atoms from one stable configuration to another. The structure with maximum energy along the MEP approximates the transition state for the process and the energy profile itself permits estimation of the transition rates. In this work we describe a computationally efficient algorithm for the identification of minimum energy paths in complicated biosystems. The algorithm is a hybrid of the nudged elastic band (NEB) and string methods. It has been implemented in the pDynamo simulation program and tested by examining elementary steps in the reaction mechanisms of three enzymes: citrate synthase, RasGAP, and lactate dehydrogenase. Good agreement is found for the energies and geometries of the species along the reaction profiles calculated using the new algorithm and previous versions of the NEB and string techniques, and also those obtained by the common method of adiabatic exploration of the potential energy surface as a function of predefined reaction coordinates. Precisely refined structures of the saddle points along the paths may be subsequently obtained with the climbing image variant of the NEB algorithm. Directions in which the utility of the methods that we have implemented can be further improved are discussed.

Design and synthesis of novel lactate dehydrogenase A inhibitors by fragment-based lead generation

Lactate dehydrogenase A (LDHA) catalyzes the conversion of pyruvate to lactate, utilizing NADH as a cofactor. It has been identified as a potential therapeutic target in the area of cancer metabolism. In this manuscript we report our progress using fragment-based lead generation (FBLG), assisted by X-ray crystallography to develop small molecule LDHA inhibitors. Fragment hits were identified through NMR and SPR screening and optimized into lead compounds with nanomolar binding affinities via fragment linking. Also reported is their modification into cellular active compounds suitable for target validation work.

Crowding induces differences in the diffusion of thermophilic and mesophilic proteins: a new look at neutron scattering results

The dynamical basis underlying the increased thermal stability of thermophilic proteins remains uncertain. Here, we challenge the new paradigm established by neutron scattering experiments in solution, in which the adaptation of thermophilic proteins to high temperatures lies in the lower sensitivity of their flexibility to temperature changes. By means of a combination of molecular dynamics and Brownian dynamics simulations, we report a reinterpretation of those experiments and show evidence that under crowding conditions, such as in vivo, thermophilic and homolog mesophilic proteins have diffusional properties with different thermal behavior.

Efficient conversion of phenylpyruvic acid to phenyllactic acid by using whole cells of Bacillus coagulans SDM

Background

Phenyllactic acid (PLA), a novel antimicrobial compound with broad and effective antimicrobial activity against both bacteria and fungi, can be produced by many microorganisms, especially lactic acid bacteria. However, the concentration and productivity of PLA have been low in previous studies. The enzymes responsible for conversion of phenylpyruvic acid (PPA) into PLA are equivocal.

Methodology/Principal Findings

A novel thermophilic strain, Bacillus coagulans SDM, was isolated for production of PLA. When the solubility and dissolution rate of PPA were enhanced at a high temperature, whole cells of B. coagulans SDM could effectively convert PPA into PLA at a high concentration (37.3 g l⁻¹) and high productivity (2.3 g l⁻¹ h⁻¹) under optimal conditions. Enzyme activity staining and kinetic studies identified NAD-dependent lactate dehydrogenases as the key enzymes that reduced PPA to PLA.

Conclusions/Significance

Taking advantage of the thermophilic character of B. coagulans SDM, a high yield and productivity of PLA were obtained. The enzymes involved in PLA production were identified and characterized, which makes possible the rational design and construction of microorganisms suitable for PLA production with metabolic engineering.

Toward Identification of the reaction coordinate directly from the transition state ensemble using the kernel PCA method

We propose a new method for analyzing an ensemble of transition states to extract components of the reaction coordinate. We use the kernel principal component analysis (kPCA), which is a generalization of the ordinary PCA that does not make a linearization approximation We applied this method to a TPS study of human LDH we had previously published [Quaytman, S.; Schwartz, S. D. Proc. Natl. Acad. Sci. U.S.A.2007, 104, 12253] and extracted a reasonable representation for the reaction coordinate.

Comparison studies of the human heart and Bacillus stearothermophilus lactate dehydrogreanse by transition path sampling

Transition path sampling is a well-known technique that generates reactive paths ensembles. Due to the atomic detail of these reactive paths, information about chemical mechanisms can be obtained. We present here a comparative study of Bacillus stearothermophilus and human heart homologues of lactate dehydrogenase (LDH). A comparison of the transition path ensemble of both enzymes revealed that small differences in the active site reverses the order of the particle transfer of the chemical step. Whereas the hydride transfer preceded the proton transfer in the human heart LDH, the order is reversed in the Bacillus stearothermophilis homologue (in the direction of pyruvate to lactate). In addition, transition state analysis revealed that the dividing region that separates reactants and products, the separatrix, is likely wider for B. stearothermophilis LDH as compared to human heart LDH. This would indicate a more variable transition process in the Bacillus enzyme.

Comprehensive structural classification of ligand-binding motifs in proteins

Comprehensive knowledge of protein-ligand interactions should provide a useful basis for annotating protein functions, studying protein evolution, engineering enzymatic activity, and designing drugs. To investigate the diversity and universality of ligand-binding sites in protein structures, we conducted the all-against-all atomic-level structural comparison of over 180,000 ligand-binding sites found in all the known structures in the Protein Data Bank by using a recently developed database search and alignment algorithm. By applying a hybrid top-down-bottom-up clustering analysis to the comparison results, we determined approximately 3000 well-defined structural motifs of ligand-binding sites. Apart from a handful of exceptions, most structural motifs were found to be confined within single families or superfamilies, and to be associated with particular ligands. Furthermore, we analyzed the components of the similarity network and enumerated more than 4000 pairs of structural motifs that were shared across different protein folds.

Theoretical site-directed mutagenesis: Asp168Ala mutant of lactate dehydrogenase

Molecular simulations based on the use of hybrid quantum mechanics/molecular mechanics methods are able to provide detailed information about the complex enzymatic reactions and the consequences of specific mutations on the activity of the enzyme. In this work, the reduction of pyruvate to lactate catalysed by wild-type and Asp168Ala mutant lactate dehydrogenase (LDH) has been studied by means of simulations using a very flexible molecular model consisting of the full tetramer of the enzyme, together with the cofactor NADH, the substrate and solvent water molecules. Our results indicate that the Asp168Ala mutation provokes a shift in the pKa value of Glu199 that becomes unprotonated at neutral pH in the mutant enzyme. This change compensates the loss of the negative charge of Asp168, rendering a still active enzyme. Thus, our methodology gives a calculated barrier height for the Asp168Ala mutant 3 kcal mol-1 higher than that for wild-type LDH, which is in very good agreement with the experiment. The computed potential energy surfaces reveal the reaction pathways and transition structures for the wild-type and mutant enzymes. Hydride transfer is less advanced and the proton transfer is more advanced in the Asp168Ala mutant than in the wild type. This approach provides a very powerful tool for the analysis of the roles of key active-site residues.

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.

Thermal limits and adaptation in marine Antarctic ectotherms: an integrative view

A cause and effect understanding of thermal limitation and adaptation at various levels of biological organization is crucial in the elaboration of how the Antarctic climate has shaped the functional properties of extant Antarctic fauna. At the same time, this understanding requires an integrative view of how the various levels of biological organization may be intertwined. At all levels analysed, the functional specialization to permanently low temperatures implies reduced tolerance of high temperatures, as a trade-off. Maintenance of membrane fluidity, enzyme kinetic properties (Km and k(cat)) and protein structural flexibility in the cold supports metabolic flux and regulation as well as cellular functioning overall. Gene expression patterns and, even more so, loss of genetic information, especially for myoglobin (Mb) and haemoglobin (Hb) in notothenioid fishes, reflect the specialization of Antarctic organisms to a narrow range of low temperatures. The loss of Mb and Hb in icefish, together with enhanced lipid membrane densities (e.g. higher concentrations of mitochondria), becomes explicable by the exploitation of high oxygen solubility at low metabolic rates in the cold, where an enhanced fraction of oxygen supply occurs through diffusive oxygen flux. Conversely, limited oxygen supply to tissues upon warming is an early cause of functional limitation. Low standard metabolic rates may be linked to extreme stenothermy. The evolutionary forces causing low metabolic rates as a uniform character of life in Antarctic ectothermal animals may be linked to the requirement for high energetic efficiency as required to support higher organismic functioning in the cold. This requirement may result from partial compensation for the thermal limitation of growth, while other functions like hatching, development, reproduction and ageing are largely delayed. As a perspective, the integrative approach suggests that the patterns of oxygen- and capacity-limited thermal tolerance are linked, on one hand, with the capacity and design of molecules and membranes, and, on the other hand, with life-history consequences and lifestyles typically seen in the permanent cold. Future research needs to address the detailed aspects of these interrelationships.

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.

Interaction between mammalian glyceraldehyde-3-phosphate dehydrogenase and L-lactate dehydrogenase from heart and muscle

The exceptionally high protein concentration in living cells can favor functional protein-protein interactions that can be difficult to detect with purified proteins. In this study we describe specific interactions between mammalian D-glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and L-lactate dehydrogenase (LDH) isozymes from heart and muscle. We use poly(ethylene-glycol) (PEG)-induced coprecipitation and native agarose electrophoresis as two independent methods uniquely suited to mimic some of the conditions that can favor protein-protein interaction in living cells. We found that GAPDH interacts with heart or muscle isozymes of LDH with approximately one-to-one stoichiometry. The interaction is specific; GAPDH shows interaction with two LDH isozymes that have very different net charge and solubility in PEG solution, while no interaction is observed with GAPDH from other species, other NAD(H) dehydrogenases, or other proteins that have very similar net charge and molecular mass. Analytical ultracentrifugation showed that the LDH and GAPDH complex is insoluble in PEG solution. The interaction is abolished by saturation with NADH, but not by saturation with NAD(+) in correlation with GAPDH solubility in PEG solution. The crystal structures show that GAPDH and LDH isozymes share complementary size, shape, and electric potential surrounding the active sites. The presented results suggest that GAPDH and LDH have a functional interaction that can affect NAD(+)/NADH metabolism and glycolysis in living cells.

Dependence of enzyme reaction mechanism on protonation state of titratable residues and QM level description: lactate dehydrogenase

We have studied the dependence of the chemical reaction mechanism of L-lactate dehydrogenase (LDH) on the protonation state of titratable residues and on the level of the quantum mechanical (QM) description by means of hybrid quantum-mechanical/molecular-mechanical (QM/MM) methods; this methodology has allowed clarification of the timing of the hydride transfer and proton transfer components that hitherto had not been possible to state definitively.

The βαβαβ elementary supersecondary structure of the Rossmann fold from porcine lactate dehydrogenase exhibits characteristics of a molten globule

Protein classifications show that the Rossmann fold, which consists of two betaalphabetaalphabeta motifs (BABAB) related by a rough twofold axis, is the most populated alphabeta fold, and that the betaalphabeta submotif (BAB) is a widespread elementary structural arrangement. Herein, we report MD simulations, circular dichroism and NMR analyses on BAB and BABAB from porcine lactate dehydrogenase to evaluate their intrinsic stability. Our results demonstrate that BAB is not stable in solution and is not a folding nucleus. We also find that BABAB, despite its appearance of a functional and structural unit, is not an independent and thermodynamically stable folding unit. Rather, we show that BABAB retains most native secondary structure but very little tertiary structure, thus displaying characteristics of a molten globule.

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.

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

The dynamical nature of the binding of a substrate surrogate to lactate dehydrogenase is examined on the nanoseconds to milliseconds timescale by laser-induced temperature-jump relaxation spectroscopy. Fluorescence emission of the nicotinamide group of bound NADH is used to define the pathway and kinetics of substrate binding. Assignment of specific kinetic states and elucidation of their structures are accomplished using isotope edited infrared absorption spectroscopy. Such studies are poised to yield a detailed picture of the coupling of protein dynamics to function.

Adaptation of enzymes to temperature: searching for basic “strategies”

The pervasive influence of temperature on biological systems necessitates a suite of temperature--compensatory adaptations that span all levels of biological organization--from behavior to fine-scale molecular structure. Beginning about 50 years ago, physiological studies conducted with whole organisms or isolated tissues, by such pioneers of comparative thermal physiology as V.Ya. Alexandrov, T.H. Bullock, F.E.J. Fry, H. Precht, C.L. Prosser, and P.F. Scholander, began to document in detail the abilities of ectothermic animals to sustain relatively similar rates of metabolic activity at widely different temperatures of adaptation or acclimation. These studies naturally led to investigation of the roles played by enzymatic proteins in metabolic temperature compensation. Peter Hochachka's laboratory became an epicenter of this new focus in comparative physiology. The studies of the enzyme lactate dehydrogenase (LDH) that he initiated as a PhD student at Duke University in the mid-1960s and continued for several years at the University of British Columbia laid much of the foundation for subsequent studies of protein adaptation to temperature. Studies of orthologs of LDH have revealed the importance of conserving kinetic properties (catalytic rate constants (kcat) and Michaelis-Menten constants (Km) and structural stability during adaptation to temperature, and recently have identified the types of amino acid substitutions causing this adaptive variation. The roles of pH and low-molecular-mass organic solutes (osmolytes) in conserving the functional and structural properties of enzymes also have been elucidated using LDH. These studies, begun in Peter Hochachka's laboratory almost 40 years ago, have been instrumental in the development of a conceptual framework for the study of biochemical adaptation, a field whose origin can be traced largely to his creative influences. This framework emphasizes the complementary roles of three "strategies" of adaptation: (1) changes in amino acid sequence that cause adaptive variation in the kinetic properties and stabilities of proteins, (2) shifts in concentrations of proteins, which are mediated through changes in gene expression and protein turnover; and (3) changes in the milieu in which proteins function, which conserve the intrinsic properties of proteins established by their primary structure and modulate protein activity in response to physiological needs. This theoretical framework has helped guide research in adaptational biochemistry for many years and now stands poised to play a critical role in the post-genomic era, as physiologists grapple with the challenge of integrating the wealth of new data on gene sequences (genome), gene expression (transcriptome and proteome), and metabolic profiles (metabolome) into a realistic physiological context that takes into account the evolutionary histories and environmental relationships of species.