Ligand binding and protein dynamics in lactate dehydrogenase

L-Lactate dehydrogenase from Cyanidioschyzon merolae shows high catalytic efficiency for pyruvate reduction and is inhibited by ATP

L-Lactate is a commodity chemical used in various fields. Microorganisms have produced L-lactate via lactic fermentation using saccharides derived from crops as carbon sources. Recently, L-lactate production using microalgae, whose carbon source is carbon dioxide, has been spotlighted because the prices of the crops have increased. A red alga Cyanidioschyzon merolae produce L-lactate via lactic fermentation under dark anaerobic conditions. The L-lactate titer of C. merolae is higher than those of other microalgae but lower than those of heterotrophic bacteria. Therefore, an increase in the L-lactate titer is required in C. merolae. L-Lactate dehydrogenase (L-LDH) catalyzes the reduction of pyruvate to L-lactate during lactic fermentation. C. merolae possesses five isozymes of L-LDH. The results of previous transcriptome analysis suggested that L-LDHs are the key enzymes in the lactic fermentation of C. merolae. However, their biochemical characteristics, such as catalytic efficiency and tolerance for metabolites, have not been revealed. We compared the amino acid sequences of C. merolae L-LDHs (CmLDHs) and characterized one of the isozymes, CmLDH1. BLAST analysis revealed that the sequence similarities of CmLDH1 and the other isozymes were above 99%. The catalytic efficiency of CmLDH1 under its optimum conditions was higher than those of L-LDHs of other organisms. ATP decreased the affinity and turnover number of CmLDH1 for NADH. These findings contribute to understanding the characteristics of L-LDHs of microalgae and the regulatory mechanisms of lactic fermentation in C. merolae.

[Identification of the binding proteins of organic acid metabolites by matrix thermal shift assay]

Organic acid metabolites exhibit acidic properties. These metabolites serve as intermediates in major carbon metabolic pathways and are involved in several biochemical pathways, including the tricarboxylic acid (TCA) cycle and glycolysis. They also regulate cellular activity and play crucial roles in epigenetics, tumorigenesis, and cellular signal transduction. Knowledge of the binding proteins of organic acid metabolites is crucial for understanding their biological functions. However, identifying the binding proteins of these metabolites has long been a challenging task owing to the transient and weak nature of their interactions. Moreover, traditional methods are unsuitable for the structural modification of the ligands of organic acid metabolites because these metabolites have simple and similar structures. Even minor structural modifications can significantly affect protein interactions. Thermal proteome profiling (TPP) provides a promising avenue for identifying binding proteins without the need for structural modifications. This approach has been successfully applied to the identification of the binding proteins of several metabolites. In this study, we investigated the binding proteins of two TCA cycle intermediates, i.e., succinate and fumarate, and lactate, an end-product of glycolysis, using the matrix thermal shift assay (mTSA) technique. This technique involves combining single-temperature (52 ℃) TPP and dose-response curve analysis to identify ligand-binding proteins with high levels of confidence and determine the binding affinity between ligands and proteins. To this end, HeLa cells were lysed, followed by protein desalting to remove endogenous metabolites from the cell lysates. The desalted cell lysates were treated with fumarate or succinate at final concentrations of 0.004, 0.04, 0.4, and 2 mmol/L in the experimental groups or 2 mmol/L sodium chloride in the control group. Considering that the cellular concentration of lactate can be as high as 2-30 mmol/L, we then applied lactate at final concentrations of 0.2, 1, 5, 10, and 25 mmol/L in the experimental groups or 25 mmol/L sodium chloride in the control group. Using high-sensitivity mass spectrometry coupled with data-independent acquisition (DIA) quantification, we quantified 5870, 5744, and 5816 proteins in succinate, fumarate, and lactate mTSA experiments, respectively. By setting stringent cut-off values (i.e., significance of changes in protein thermal stability (p-value)<0.001 and quality of the dose-response curve fitting (square of Pearson's correlation coefficient, R2)>0.95), multiple binding proteins for these organic acid metabolites from background proteins were confidently determined. Several known binding proteins were identified, notably fumarate hydratase (FH) as a binding protein for fumarate, and α-ketoglutarate-dependent dioxygenase (FTO) as a binding protein for both fumarate and succinate. Additionally, the affinity data for the interactions between these metabolites and their binding proteins were obtained, which closely matched those reported in the literature. Interestingly, ornithine aminotransferase (OAT), which is involved in amino acid biosynthesis, and 3-mercaptopyruvate sulfurtransferase (MPST), which acts as an antioxidant in cells, were identified as lactate-binding proteins. Subsequently, an orthogonal assay technique developed in our laboratory, the solvent-induced precipitation (SIP) technique, was used to validate the mTSA results. SIP identified OAT as the top target candidate, validating the mTSA-based finding that OAT is a novel lactate-binding protein. Although MPST was not identified as a lactate-binding protein by SIP, statistical analysis of MPST in the mTSA experiments with 10 or 25 mmol/L lactate revealed that MPST is a lactate-binding protein with a high level of confidence. Peptide-level empirical Bayes t-tests combined with Fisher's exact test also supported the conclusion that MPST is a lactate-binding protein. Lactate is structurally similar to pyruvate, the known binding protein of MPST. Therefore, assuming that lactate could potentially occupy the binding site of pyruvate on MPST. Overall, the novel binding proteins identified for lactate suggest their potential involvement in amino acid synthesis and redox balance regulation.

The mechanistic insights into different aspects of promiscuity in metalloenzymes

Enzymes are nature's ultimate machinery to catalyze complex reactions. Though enzymes are evolved to catalyze specific reactions, they also show significant promiscuity in reactions and substrate selection. Metalloenzymes contain a metal ion or metal cofactor in their active site, which is crucial in their catalytic activity. Depending on the metal and its coordination environment, the metal ion or cofactor may function as a Lewis acid or base and a redox center and thus can catalyze a plethora of natural reactions. In fact, the versatility in the oxidation state of the metal ions provides metalloenzymes with a high level of catalytic adaptability and promiscuity. In this chapter, we discuss different aspects of promiscuity in metalloenzymes by using several recent experimental and theoretical works as case studies. We start our discussion by introducing the concept of promiscuity and then we delve into the mechanistic insight into promiscuity at the molecular level.

Deciphering Evolutionary Trajectories of Lactate Dehydrogenases Provides New Insights into Allostery

Lactate dehydrogenase (LDH, EC.1.1.127) is an important enzyme engaged in the anaerobic metabolism of cells, catalyzing the conversion of pyruvate to lactate and NADH to NAD+. LDH is a relevant enzyme to investigate structure-function relationships. The present work provides the missing link in our understanding of the evolution of LDHs. This allows to explain (i) the various evolutionary origins of LDHs in eukaryotic cells and their further diversification and (ii) subtle phenotypic modifications with respect to their regulation capacity. We identified a group of cyanobacterial LDHs displaying eukaryotic-like LDH sequence features. The biochemical and structural characterization of Cyanobacterium aponinum LDH, taken as representative, unexpectedly revealed that it displays homotropic and heterotropic activation, typical of an allosteric enzyme, whereas it harbors a long N-terminal extension, a structural feature considered responsible for the lack of allosteric capacity in eukaryotic LDHs. Its crystallographic structure was solved in 2 different configurations typical of the R-active and T-inactive states encountered in allosteric LDHs. Structural comparisons coupled with our evolutionary analyses helped to identify 2 amino acid positions that could have had a major role in the attenuation and extinction of the allosteric activation in eukaryotic LDHs rather than the presence of the N-terminal extension. We tested this hypothesis by site-directed mutagenesis. The resulting C. aponinum LDH mutants displayed reduced allosteric capacity mimicking those encountered in plants and human LDHs. This study provides a new evolutionary scenario of LDHs that unifies descriptions of regulatory properties with structural and mutational patterns of these important enzymes.

Protein Dynamics and Enzymatic Catalysis

This Perspective presents a review of our work and that of others in the highly controversial topic of the coupling of protein dynamics to reaction in enzymes. We have been involved in studying this topic for many years. Thus, this perspective will naturally present our own views, but it also is designed to present an overview of the variety of viewpoints of this topic, both experimental and theoretical. This is obviously a large and contentious topic.

Design and Synthesis of New Pyrimidine-Quinolone Hybrids as Novel hLDHA Inhibitors

A battery of novel pyrimidine-quinolone hybrids was designed by docking scaffold replacement as lactate dehydrogenase A (hLDHA) inhibitors. Structures with different linkers between the pyrimidine and quinolone scaffolds (10-21 and 24−31) were studied in silico, and those with the 2-aminophenylsulfide (U-shaped) and 4-aminophenylsulfide linkers (24−31) were finally selected. These new pyrimidine-quinolone hybrids (24−31)(a−c) were easily synthesized in good to excellent yields by a green catalyst-free microwave-assisted aromatic nucleophilic substitution reaction between 3-(((2/4-aminophenyl)thio)methyl)quinolin-2(1H)-ones 22/23(a−c) and 4-aryl-2-chloropyrimidines (1−4). The inhibitory activity against hLDHA of the synthesized hybrids was evaluated, resulting IC50 values of the U-shaped hybrids 24−27(a−c) much better than the ones of the 1,4-linked hybrids 28−31(a−c). From these results, a preliminary structure−activity relationship (SAR) was established, which enabled the design of novel 1,3-linked pyrimidine-quinolone hybrids (33−36)(a−c). Compounds 35(a−c), the most promising ones, were synthesized and evaluated, fitting the experimental results with the predictions from docking analysis. In this way, we obtained novel pyrimidine-quinolone hybrids (25a, 25b, and 35a) with good IC50 values (<20 μM) and developed a preliminary SAR.

Proteomic profiling of cisplatin-resistant and cisplatin-sensitive germ cell tumour cell lines using quantitative mass spectrometry

Purpose

Advanced testicular germ cell tumours (GCT) generally have a good prognosis owing to their unique sensitivity towards cisplatin-based chemotherapies. However, cisplatin-resistant GCT have a poor outcome. Further studies are mandatory to better understand resistance mechanisms and develop therapeutic strategies for refractory GCTs.

Methods

Protein levels in cisplatin-resistant GCT cell lines of NTERA-2, NCCIT and 2102EP were analyzed by quantitative proteomic mass spectrometry (MS) in combination with stable isotope labelling by amino acids in cell culture (SILAC). Differentially abundant protein markers of acquired cisplatin resistance were validated by Western blotting. Comprehensive bioinformatical annotation using gene set enrichment analyses (GSEA) and STRING interaction analysis were performed to identify commonly affected pathways in cisplatin resistance and the data were compared to the GCT cohort of the ‘The Cancer Genome Atlas’.

Results

A total of 4375 proteins were quantified by MS, 144 of which were found to be differentially abundant between isogenic resistant and sensitive cell line pairs (24 proteins for NTERA-2, 60 proteins for NCCIT, 75 proteins for 2102EP). Western blotting confirmed regulation of key resistance-associated proteins (CBS, ANXA1, LDHA, CTH, FDXR). GSEA revealed a statistically significant enrichment of DNA repair-associated proteins in all three resistant cell lines and specific additional processes for individual cell lines.

Conclusion

High resolution MS combined with SILAC is a powerful tool and 144 significantly deregulated proteins were found in cisplatin-resistant GCT cell lines. Our study provides the largest proteomic in vitro library for cisplatin resistance in GCT, yet, enabling further studies to develop new treatment options for patients with refractory GCT.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00345-022-03936-1.

An Allosteric Inhibitory Potential of Triterpenes from Combretum racemosum on the Structural and Functional Dynamics of Plasmodium falciparum Lactate Dehydrogenase Binding Landscape

Multidrug resistance is a significant drawback in malaria treatment, and mutations in the active sites of the many critical antimalarial drug targets have remained challenging. Therefore, this has necessitated the global search for new drugs with new mechanisms of action. Plasmodium falciparum lactate dehydrogenase (pfLHD), a glycolytic enzyme, has emerged as a potential target for developing new drugs due to the parasite reliance on glycolysis for energy. Strong substrate-binding is required in pfLDH enzymatic catalysis; however, there is a lack of information on small molecules' inhibitory mechanism bound to the substrate-binding pocket. Therefore, this study investigated a potential allosteric inhibition of pfLDH by targeting the substrate-binding site. The structural and functional behaviour of madecassic acid (MA), the most promising among the six triterpenes bound to pfLDH, were unravelled using molecular dynamic simulations at 300 ns to gain insights into its mechanism of binding and inhibition and chloroquine as a standard drug. The docking studies identified that the substrate site has the preferred position for the compounds even in the absence of a co-factor. The bound ligands showed comparably higher binding affinity at the substrate site than at the co-factor site. Mechanistically, a characteristic loop implicated in the enzyme catalytic activity was identified at the substrate site. This loop accommodates key interacting residues (LYS174, MET175, LEU177 and LYS179) pivotal in the MA binding and inhibitory action. The MA-bound pfLHD average RMSD (1.60 Å) relative to chloroquine-bound pfLHD RMSD (2.00 Å) showed higher stability for the substrate pocket, explaining the higher binding affinity (-33.40 kcal/mol) observed in the energy calculations, indicating that MA exhibited profound inhibitory activity. The significant pfLDH loop conformational changes and the allostery substrate-binding landscape suggested inhibiting the enzyme function, which provides an avenue for designing antimalarial compounds in the future studies of pfLDH protein as a target.

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.

LNA oligonucleotide mediates an anti-inflammatory effect in autoimmune myocarditis via targeting lactate dehydrogenase B

Treatment of myocarditis is often limited to symptomatic treatment due to unknown pathomechanisms. In order to identify new therapeutic approaches, the contribution of locked nucleic acid antisense oligonucleotides (LNA ASOs) in autoimmune myocarditis was investigated. Hence, A/J mice were immunized with cardiac troponin I (TnI) to induce experimental autoimmune myocarditis (EAM) and treated with LNA ASOs. The results showed an unexpected anti‐inflammatory effect for one administered LNA ASO MB_1114 by reducing cardiac inflammation and fibrosis. The target sequence of MB_1114 was identified as lactate dehydrogenase B (mLDHB). For further analysis, mice received mLdhb‐specific GapmeR during induction of EAM. Here, mice receiving the mLdhb‐specific GapmeR showed increased protein levels of cardiac mLDHB and a reduced cardiac inflammation and fibrosis. The effect of increased cardiac mLDHB protein level was associated with a downregulation of genes of reactive oxygen species (ROS)‐associated proteins, indicating a reduction in ROS. Here, the suppression of murine pro‐apoptotic Bcl‐2‐associated X protein (mBax) was also observed. In our study, an unexpected anti‐inflammatory effect of LNA ASO MB_1114 and mLdhb‐specific GapmeR during induction of EAM could be demonstrated in vivo. This effect was associated with increased protein levels of cardiac mLDHB, mBax suppression and reduced ROS activation. Thus, LDHB and LNA ASOs may be considered as a promising target for directed therapy of myocarditis. Nevertheless, further investigations are necessary to clarify the mechanism of action of anti‐inflammatory LDHB‐triggered effects.

Temperature-Dependent Low-Frequency Modes in the Active Site of Bovine Carbonic Anhydrase II Probed by 2D-IR Spectroscopy

Enzyme catalysis achieves tremendous rate accelerations. Enzyme reaction centers provide a constraint geometry that preferentially binds an activated form of the substrate and thus lowers the energy barrier. However, this transition state picture neglects the flexibility of proteins and its role in enzymatic catalysis. Especially for proton transfer reactions, it has been suggested that motions of the protein modulate the donor-acceptor distance and prepare a tunneling-ready state. We report the detection of frequency fluctuations of an azide anion (N₃^-) bound in the active site of the protein carbonic anhydrase II, where a low-frequency mode of the protein has been proposed to facilitate proton transfer over two water molecules during the catalyzed reaction. 2D-IR spectroscopy resolves an underdamped low-frequency mode at about 1 THz (30 cm^-1). We find its frequency to be viscosity- and temperature-dependent and to decrease by 6 cm^-1 between 230 and 320 K, reporting the softening of the mode's potential.

Acceleration of catalysis in dihydrofolate reductase by transient, site-specific photothermal excitation

We have studied the role of protein dynamics in chemical catalysis in the enzyme dihydrofolate reductase (DHFR), using a pump-probe method that employs pulsed-laser photothermal heating of a gold nanoparticle (AuNP) to directly excite a local region of the protein structure and transient absorbance to probe the effect on enzyme activity. Enzyme activity is accelerated by pulsed-laser excitation when the AuNP is attached close to a network of coupled motions in DHFR (on the FG loop, containing residues 116-132, or on a nearby alpha helix). No rate acceleration is observed when the AuNP is attached away from the network (distal mutant and His-tagged mutant) with pulsed excitation, or for any attachment site with continuous wave excitation. We interpret these results within an energy landscape model in which transient, site-specific addition of energy to the enzyme speeds up the search for reactive conformations by activating motions that facilitate this search.

Pressure tolerance of deep-sea enzymes can be evolved through increasing volume changes in protein transitions: a study with lactate dehydrogenases from abyssal and hadal fishes

We explore the principles of pressure tolerance in enzymes of deep-sea fishes using lactate dehydrogenases (LDH) as a case study. We compared the effects of pressure on the activities of LDH from hadal snailfishes Notoliparis kermadecensis and Pseudoliparis swirei with those from a shallow-adapted Liparis florae and an abyssal grenadier Coryphaenoides armatus. We then quantified the LDH content in muscle homogenates using mass-spectrometric determination of the LDH-specific conserved peptide LNLVQR. Existing theory suggests that adaptation to high pressure requires a decrease in volume changes in enzymatic catalysis. Accordingly, evolved pressure tolerance must be accompanied with an important reduction in the volume change associated with pressure-promoted alteration of enzymatic activity ( Δ V PP ∘ ). Our results suggest an important revision to this paradigm. Here, we describe an opposite effect of pressure adaptation-a substantial increase in the absolute value of Δ V PP ∘ in deep-living species compared to shallow-water counterparts. With this change, the enzyme activities in abyssal and hadal species do not substantially decrease their activity with pressure increasing up to 1-2 kbar, well beyond full-ocean depth pressures. In contrast, the activity of the enzyme from the tidepool snailfish, L. florae, decreases nearly linearly from 1 to 2500 bar. The increased tolerance of LDH activity to pressure comes at the expense of decreased catalytic efficiency, which is compensated with increased enzyme contents in high-pressure-adapted species. The newly discovered strategy is presumably used when the enzyme mechanism involves the formation of potentially unstable excited transient states associated with substantial changes in enzyme-solvent interactions.

Ferrihydrite nanoparticles insights: Structural characterization, lactate dehydrogenase binding and virtual screening assay

The binding between the enzyme lactate dehydrogenase (LDH) and ferrihydrite nanoparticles (Fh-NPs) was investigated by means of small-angle neutron scattering (SANS), Fourier-transform infrared (FTIR) spectroscopy, fluorescence and Förster resonance energy transfer (FRET) and molecular docking. Fh-NPs - LDH compounds of dimensions under 100 nm are formed. The conformational changes and the mechanism of interaction between LDH and Fh-NPs simple and doped with Cu and Co, and the effect of these NPs on the thermal denaturation of LDH were monitored. The quenching mechanism is static, the binding occurring with moderate affinity, being mainly driven by hydrogen bonding and van der Waals forces. FRET occurs at a minimal distance of 2.55 nm. Thermal denaturation of LDH in the presence of simple and doped Fh-NPs shows that the thermodynamic parameters of protein unfolding are significantly changed with temperature. The denaturation temperature of LDH shifts to higher values in the presence of all Fh-NPs, than in the case of simple LDH. The docking approach estimates the energy corresponding to the best fit of the ferrihydrite in the LDH binding site near Trp. These results have direct implications on the uses of the complex of LDH with Fh-NPs in various biochemical, biological, or clinical applications.

Substrate Channeling via a Transient Protein-Protein Complex: The case of D-Glyceraldehyde-3-Phosphate Dehydrogenase and L-Lactate Dehydrogenase

Substrate channeling studies have frequently failed to provide conclusive results due to poor understanding of this subtle phenomenon. We analyzed the mechanism of NADH-channeling from D-glyceraldehyde-3-phosphate dehydrogenase (GAPDH) to L-lactate Dehydrogenase (LDH) using enzymes from different cells. Enzyme kinetics studies showed that LDH activity with free NADH and GAPDH-NADH complex always take place in parallel. The channeling is observed only in assays that mimic cytosolic conditions where free NADH concentration is negligible and the GAPDH-NADH complex is dominant. Molecular dynamics and protein-protein interaction studies showed that LDH and GAPDH can form a leaky channeling complex only at the limiting NADH concentrations. Surface calculations showed that positive electric field between the NAD(H) binding sites on LDH and GAPDH tetramers can merge in the LDH-GAPDH complex. NAD(H)-channeling within the LDH-GAPDH complex can be an extension of NAD(H)-channeling within each tetramer. In the case of a transient LDH-(GAPDH-NADH) complex, the relative contribution from the channeled and the diffusive paths depends on the overlap between the off-rates for the LDH-(GAPDH-NADH) complex and the GAPDH-NADH complex. Molecular evolution or metabolic engineering protocols can exploit substrate channeling for metabolic flux control by fine-tuning substrate-binding affinity for the key enzymes in the competing reaction paths.

Crowding-Induced Uncompetitive Inhibition of Lactate Dehydrogenase: Role of Entropic Pushing

The cell is an extremely complex environment, notably highly crowded, segmented, and confining. Overall, there is overwhelming and ever-growing evidence that to understand how biochemical reactions proceed in vivo, one cannot separate the biochemical actors from their environment. Effects such as excluded volume, obstructed diffusion, weak nonspecific interactions, and fluctuations all team up to steer biochemical reactions often very far from what is observed in ideal conditions. In this paper, we use Ficoll PM70 and PEG 6000 to build an artificial crowded milieu of controlled composition and density in order to assess how such environments influence the biocatalytic activity of lactate dehydrogenase (LDH). Our measurements show that the normalized apparent affinity and maximum velocity decrease in the same fashion, a behavior reminiscent of uncompetitive inhibition, with PEG resulting in the largest reduction. In line with previous studies on other enzymes of the same family, and in agreement with the known role of a surface loop involved in enzyme isomerization and regulation of access to the active site, we suggest that the crowding matrix interferes with the conformational ensemble of the enzyme. This likely results in both impaired enzyme-complex isomerization and thwarted product release. Molecular dynamics simulations confirm that excluded-volume effects lead to an entropic force that effectively tends to push the loop closed, thereby effectively shifting the conformational ensemble of the enzyme in favor of a more stable complex isoform. Overall, our study substantiates the idea that most biochemical kinetics cannot be fully explained without including the subtle action of the environment where they take place naturally, in particular accounting for important factors such as excluded-volume effects and also weak nonspecific interactions when present, confinement, and fluctuations.

Osmolytes modify protein dynamics and function of tetrameric lactate dehydrogenase upon pressurization

We present a study of the combined effects of natural cosolvents (TMAO, glycine, urea) and pressure on the activity of the tetrameric enzyme lactate dehydrogenase (LDH). To this end, high-pressure stopped-flow methodology in concert with fast UV/Vis spectroscopic detection of product formation was applied. To reveal possible pressure effects on the stability and dynamics of the enzyme, FTIR spectroscopic and neutron scattering measurements were carried out. In neat buffer solution, the catalytic turnover number of the enzyme, k_cat, increases up to 1000 bar, the pressure range where dissociation of the tetrameric species to dimers sets in. Accordingly, we obtain a negative activation volume, ΔV^# = -45.3 mL mol^-1. Further, the enzyme substrate complex has a larger volume compared to the enzyme and substrate in the unbound state. The neutron scattering data show that changes in the fast internal dynamics of the enzyme are not responsible for the increase of k_cat upon compression. Whereas the magnitude of k_cat is similar in the presence of the osmolytes, the pressure of deactivation is modulated by the addition of cosolvents. TMAO and glycine increase the pressure of deactivation, and in accordance with the observed stabilizing effect both cosolvents exhibit against denaturation and/or dissociation of proteins. While urea does not markedly affect the magnitude of the Michaelis constant, K_M, both 1 M TMAO and 1 M glycine exhibit smaller K_M values of about 0.07 mM and 0.05 mM below about 1 kbar. Such positive effect on the substrate affinity could be rationalized by the effect the two cosolutes impose on the thermodynamic activities of the reactants, which reflect changes in water-mediated intermolecular interactions. Our data show that the intracellular milieu, i.e., the solution conditions that have evolved, may be sufficient to maintain enzymatic activity under extreme environmental conditions, including the whole pressure range encountered on Earth.

Glycation profile of minor abundant erythrocyte proteome across varying glycemic index in diabetes mellitus

PURPOSE

Long-term glycemic index in patients with diabetes mellitus (DM) is measured by glycated hemoglobin (HbA_1c) besides blood glucose. In DM, the primary amino groups of proteins get glycated via non-enzymatic post-translational modification. This study aims at identifying and characterizing site-specific glycation of erythrocyte proteome across varying glycemic index in patients with DM.

EXPERIMENTS

We isolated the glycated erythrocyte proteome devoid of hemoglobin from control and diabetic samples using boronate affinity chromatography. Proteomic analysis was performed using nanoLC/ESI-MS proteomics platform. The site-specific modification on different proteins was deciphered using a customized database.

RESULTS

We report 37 glycated proteins identified and characterized from samples with HbA_1c of 6%, 8%, 12%, and 16%. Our results show that both extent and site-specific modification of proteins increased with increasing HbA_1c. The observed residue-specific modifications of catalase, peroxiredoxin, carbonic anhydrase, lactate dehydrogenase B and delta-aminolevulinic acid dehydratase were correlated with the literature report on their functional disorder in DM.

CONCLUSIONS

and clinical relevance: 37 glycated erythrocyte proteins apart from hemoglobin were characterized from DM patient samples with varying HbA_1c values. We correlated the site-specific glycation and associated functional disorder of five representative proteins. However, the clinical correlation with the observed modifications needs further investigation.

Profiling of Methylglyoxal Blood Metabolism and Advanced Glycation End-Product Proteome Using a Chemical Probe

Methylglyoxal (MG) is quantitatively the most important precursor to advanced glycation end-products (AGEs), and evidence is accumulating that it is also a causally linked to diabetes and aging related diseases. Living systems primarily reside on the glyoxalase system to detoxify MG into benign d-lactate. The flux to either glycation or detoxification, accordingly, is a key parameter for how well a system handles the ubiquitous glyoxal burden. Furthermore, insight into proteins and in particular their individual modification sites are central to understanding the involvement of MG and AGE in diabetes and aging related diseases. Here, we present a simple method to simultaneously monitor the flux of MG both to d-lactate and to protein AGE formation in a biological sample by employing an alkyne-labeled methylglyoxal probe. We apply the method to blood and plasma to demonstrate the impact of blood cell glyoxalase activity on plasma protein AGE formation. We move on to isolate proteins modified by the MG probe and accordingly can present the first general inventory of more than 100 proteins and 300 binding sites of the methylglyoxal probe on plasma as well as erythrocytic proteins. Some of the data could be validated against a number of in vivo and in vitro targets for advanced glycation previously known from the literature; the majority of proteins and specific sites however were previously unknown and may guide future research into MG and AGE to elucidate how these are functionally linked to diabetic disease and aging.

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.

Promoting Vibrations and the Function of Enzymes. Emerging Theoretical and Experimental Convergence

A complete understanding of enzyme catalysis requires knowledge of both transition state features and the detailed motions of atoms that cause reactant molecules to form and traverse the transition state. The seeming intractability of the problem arises from the femtosecond lifetime of chemical transition states, preventing most experimental access. Computational chemistry is admirably suited to short time scale analysis but can be misled by inappropriate starting points or by biased assumptions. Kinetic isotope effects provide an experimental approach to transition state structure and a method for obtaining transition state analogues but, alone, do not inform how that transition state is reached. Enzyme structures with transition state analogues provide computational starting points near the transition state geometry. These well-conditioned starting points, combined with the unbiased computational method of transition path sampling, provide realistic atomistic motions involved in transition state formation and passage. In many, but not all, enzymatic systems, femtosecond local protein motions near the catalytic site are linked to transition state formation. These motions are not inherently revealed by most approaches of transition state theory, because transition state theory replaces dynamics with the statistics of the transition state. Experimental and theoretical convergence of the link between local catalytic site vibrational modes and catalysis comes from heavy atom ("Born-Oppenheimer") enzymes. Fully labeled and catalytic site local heavy atom labels perturb the probability of finding enzymatic transition states in ways that can be analyzed and predicted by transition path sampling. Recent applications of these experimental and computational approaches reveal how subpicosecond local catalytic site protein modes play important roles in creating the transition state.

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.

Retention of functional characteristics of glutathione-S-transferase and lactate dehydrogenase-A in fusion protein

A paradigm shift toward fusion proteins to render multiple functionalities and applications on a single platform has been incurred in enzyme based diagnosis. Herein, we report development and systematic characterizations of glutathione-S-transferase (GST) and human lactate dehydrogenase A (hLDHA) in a fusion protein (GST-hLDHA) to achieve functional activities of GST and hLDHA simultaneously. The GST-pGEX-4T-2 vector system was used for cloning and purification of hLDHA, utilizing the affinity based interaction between GST and GSH in column chromatography. Bacterially purified protein was subjected to the Western blot analysis and structural analysis by circular dichroism spectroscopy, which revealed intact structural framework of the fusion construct. Kinetic characterization of the fusion GST-hLDHA protein toward GSH and NADH, suggested retention of functional activities of GST and hLDHA in fused protein as indicated by the kinetic parameters k_m and k_cat/k_m. Further analysis of effect of temperature and pH on GST-hLDHA activity revealed maximum activity around human physiological conditions (37°C and pH 8). Preservation of the structural and functional characteristics of the fusion enzyme paves the way for potential application for the detection of NADH and GSH in conjunction as biomarkers for cancer diagnosis.

Neural Stem Cell Death Mechanisms Induced by Amyloid Beta

Background and Purpose

Amyloid beta (Aβ) is the main component of amyloid plaques, which are deposited in the brains of patients with Alzheimer's disease (AD). Biochemical and animal studies support the central role of Aβ in AD pathogenesis. Despite several investigations focused on the pathogenic mechanisms of Aβ, it is still unclear how Aβ accumulates in the central nervous system and subsequently initiates the disease at the cellular level. In this study, we investigated the pathogenic mechanisms of Aβ using proteomics and antibody microarrays.

Methods

To evaluate the effect of Aβ on neural stem cells (NSCs), we treated primary cultured cortical NSCs with several doses of Aβ for 48 h. A 3-[4,5-dimethylthiazol-2-yl]-2,5-diphenyltetrazolium bromide assay, trypan blue staining, and bromodeoxyuridine cell proliferation assay were performed. We detected several intracellular proteins that may be associated with Aβ by proteomics and Western blotting analysis.

Results

Various viability tests showed that Aβ decreased NSCs viability and cell proliferation in a concentration-dependent manner. Aβ treatment significantly decreased lactate dehydrogenase B, high-mobility group box 1, aldolase C, Ezrin, and survival signals including phosphorylated phosphoinositide 3-kinase, Akt, and glycogen synthase kinase-3β.

Conclusions

These results suggest that several factors determined by proteomics and Western blot hold the clue to Aβ pathogenesis. Further studies are required to investigate the role of these factors.

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.

Thermodynamic and Structural Adaptation Differences between the Mesophilic and Psychrophilic Lactate Dehydrogenases

The thermodynamics of substrate binding and enzymatic activity of a glycolytic enzyme, lactate dehydrogenase (LDH), from both porcine heart, phLDH (Sus scrofa; a mesophile), and mackerel icefish, cgLDH (Chamapsocephalus gunnari; a psychrophile), were investigated. Using a novel and quite sensitive fluorescence assay that can distinguish protein conformational changes close to and distal from the substrate binding pocket, a reversible global protein structural transition preceding the high-temperature transition (denaturation) was surprisingly found to coincide with a marked change in enzymatic activity for both LDHs. A similar reversible structural transition of the active site structure was observed for phLDH but not for cgLDH. An observed lower substrate binding affinity for cgLDH compared to that for phLDH was accompanied by a larger contribution of entropy to ΔG, which reflects a higher functional plasticity of the psychrophilic cgLDH compared to that of the mesophilic phLDH. The natural osmolyte, trimethylamine N-oxide (TMAO), increases stability and shifts all structural transitions to higher temperatures for both orthologs while simultaneously reducing catalytic activity. The presence of TMAO causes cgLDH to adopt catalytic parameters like those of phLDH in the absence of the osmolyte. Our results are most naturally understood within a model of enzyme dynamics whereby different conformations of the enzyme that have varied catalytic parameters (i.e., binding and catalytic proclivity) and whose population profiles are temperature-dependent and influenced by osmolytes interconvert among themselves. Our results also show that adaptation can be achieved by means other than gene mutations and complements the synchronic evolution of the cellular milieu.

Multiple helical conformations of the helix-turn-helix region revealed by NOE-restrained MD simulations of tryptophan aporepressor, TrpR

The nature of flexibility in the helix-turn-helix region of E. coli trp aporepressor has been unexplained for many years. The original ensemble of nuclear magnetic resonance (NMR structures showed apparent disorder, but chemical shift and relaxation measurements indicated a helical region. Nuclear Overhauser effect (NOE) data for a temperature-sensitive mutant showed more helical character in its helix-turn-helix region, but nevertheless also led to an apparently disordered ensemble. However, conventional NMR structure determination methods require all structures in the ensemble to be consistent with every NOE simultaneously. This work uses an alternative approach in which some structures of the ensemble are allowed to violate some NOEs to permit modeling of multiple conformational states that are in dynamic equilibrium. Newly measured NOE data for wild-type aporepressor are used as time-averaged distance restraints in molecular dynamics simulations to generate an ensemble of helical conformations that is more consistent with the observed NMR data than the apparent disorder in the previously reported NMR structures. The results indicate the presence of alternating helical conformations that provide a better explanation for the flexibility of the helix-turn-helix region of trp aporepressor. Structures representing these conformations have been deposited with PDB ID: 5TM0. Proteins 2017; 85:731-740. © 2016 Wiley Periodicals, Inc.

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.

Structural and dynamical correlations in PfHGXPRT oligomers: A molecular dynamics simulation study

PfHGXPRT is a key enzyme involved in purine nucleotide salvage pathway of the malarial parasite, Plasmodium falciparum. Atomistic molecular dynamics simulations have been performed on two types of PfHGXPRT dimers (D1 and D3) and its tetramer in their apo and ligand-bound states. A significant event in the catalytic cycle is the dynamics of a gate that provides access for the ligand molecules to the reaction center. The gate is formed by loops II and IV, the former being the most flexible. Large amplitude conformational changes have been observed in active site loop II. Upon complete occupancy of the active site, loop II gets stabilized due to specific interactions between its residues and the ligand molecules. Remote loop, X, is seen to be less fluxional in the D3 dimer than in D1 which is rationalized as due to the greater number of inter-subunit contacts in the former. The presence of ligand molecules in subunits of the tetramer further reduces the flexibility of loop X epitomizing a communication between this region and the active sites in the tetramer. These observations are in accordance with the outcomes of several experimental investigations. Participation of loop X in the oligomerization process has also been discerned. Between the two types of dimers in solution, D1 tetramerizes readily and thus would not be present as free dimers. We conjecture an equilibrium to exist between D3 and the tetramer in solution; upon binding of the ligand molecules to the D3 dimer, this equilibrium shifts toward the tetramer.

Lactate modulates the intracellular pH sensitivity of human TREK1 channels

Tissue acidosis and high lactate concentrations are associated with cerebral ischaemia. The degree of acidosis is dependent on circulating glucose concentration, hyperglycaemia being associated with increased acidosis. Among other agents, lactate and protons have been shown to activate the leak potassium channel; TREK1 (TWIK related potassium channel 1) from the intracellular side and its increased activity is implicated in tolerance towards ischaemic cell damage. In the present study, we show that ischaemic concentrations of lactate (30 mM) at pH 7.0 and 6.5, commonly observed during ischemia, cause robust potentiation of human TREK1 (hTREK1) activity at single-channel level in cell-free inside-out membrane patches, while 30 mM lactate at pH 6.0 to 5.5, commonly observed during hyperglycaemic ischemia, reduces hTREK1 channel activity significantly. The biphasic effect of 30 mM lactate (ischaemic concentrations) on modulation of hTREK1 by varying pH conditions is specific since basal concentrations of lactate (3 mM) and 30 mM pyruvate at pH 7.0 and 5.5 failed to show similar effect as lactate. Experiments with deletion and point mutants of hTREK1 channel suggest that lactate changes the pH modulation of hTREK1 by interacting differently with the histidine residue at 328th position (H328) above and below its pKa (∼6.0) in the intracellular carboxyl-terminal domain of TREK1. This lactate-induced pH modulation of hTREK1 is absent in C-terminal deletion mutant, CTDΔ100, and is similar in E321A-hTREK1 mutant as in wild-type hTREK1 suggesting that it is independent of pH-sensitive glutamate residue at 321st position. Such a differential pH-dependent effect of lactate on an ion channel function has not been reported earlier and has important implications in different stages of ischaemia.

Inhibitory effects of ionic liquids on the lactic dehydrogenase activity

Ionic liquids (ILs) were widely used in scientific and industrial application and have been reported to possess potential toxicity to the environment and human health. The effects of six typical N-methylimidazolium-based ILs ([Cnmim]X, n=4, 6, 8; X=Br(-), Cl(-), BF4(-), CF3SO3(-)) on the lactic dehydrogenase (LDH) activity and the molecular interaction mechanism of ILs and the LDH were investigated with the aid of spectroscopic techniques. Experimental results showed that the LDH activity was inhibited in the presence of ILs. For the ILs with the same anion but different cations, their inhibitory ability on the LDH activity increased with increasing the alkyl chain length on the IL cation. Thermodynamic parameters, enthalpy change (ΔH) and entropy change (ΔS) were obtained by analyzing the fluorescence behavior of LDH with the addition of ILs. Both positive ΔH and ΔS suggested that hydrophobicity was the major driven force in the interaction process as expected.

Ischaemic concentrations of lactate increase TREK1 channel activity by interacting with a single histidine residue in the carboxy terminal domain

KEY POINTS

The physiological metabolite, lactate and the two-pore domain leak potassium channel, TREK1 are known neuroprotectants against cerebral ischaemia. However, it is not known whether lactate interacts with TREK1 channel to provide neuroprotection. In this study we show that lactate increases TREK1 channel activity and hyperpolarizes CA1 stratum radiatum astrocytes in hippocampal slices. Lactate increases open probability and decreases longer close time of the human (h)TREK1 channel in a concentration dependent manner. Lactate interacts with histidine 328 (H328) in the carboxy terminal domain of hTREK1 channel to decrease its dwell time in the longer closed state. This interaction was dependent on the charge on H328. Lactate-insensitive mutant H328A hTREK1 showed pH sensitivity similar to wild-type hTREK1, indicating that the effect of lactate on hTREK1 is independent of pH change. A rise in lactate concentration and the leak potassium channel TREK1 have been independently associated with cerebral ischaemia. Recent literature suggests lactate to be neuroprotective and TREK1 knockout mice show an increased sensitivity to brain and spinal cord ischaemia; however, the connecting link between the two is missing. Therefore we hypothesized that lactate might interact with TREK1 channels. In the present study, we show that lactate at ischaemic concentrations (15-30 mm) at pH 7.4 increases TREK1 current in CA1 stratum radiatum astrocytes and causes membrane hyperpolarization. We confirm the intracellular action of lactate on TREK1 in hippocampal slices using monocarboxylate transporter blockers and at single channel level in cell-free inside-out membrane patches. The intracellular effect of lactate on TREK1 is specific since other monocarboxylates such as pyruvate and acetate at pH 7.4 failed to increase TREK1 current. Deletion and point mutation experiments suggest that lactate decreases the longer close dwell time incrementally with increase in lactate concentration by interacting with the histidine residue at position 328 (H328) in the carboxy terminal domain of the TREK1 channel. The interaction of lactate with H328 is dependent on the charge on the histidine residue since isosteric mutation of H328 to glutamine did not show an increase in TREK1 channel activity with lactate. This is the first demonstration of a direct effect of lactate on ion channel activity. The action of lactate on the TREK1 channel signifies a separate neuroprotective mechanism in ischaemia since it was found to be independent of the effect of acidic pH on channel activity.

The Regulation and Function of Lactate Dehydrogenase A: Therapeutic Potential in Brain Tumor

There are over 120 types of brain tumor and approximately 45% of primary brain tumors are gliomas, of which glioblastoma multiforme (GBM) is the most common and aggressive with a median survival rate of 14 months. Despite progress in our knowledge, current therapies are unable to effectively combat primary brain tumors and patient survival remains poor. Tumor metabolism is important to consider in therapeutic approaches and is the focus of numerous research investigations. Lactate dehydrogenase A (LDHA) is a cytosolic enzyme, predominantly involved in anaerobic and aerobic glycolysis (the Warburg effect); however, it has multiple additional functions in non‐neoplastic and neoplastic tissues, which are not commonly known or discussed. This review summarizes what is currently known about the function of LDHA and identifies areas that would benefit from further exploration. The current knowledge of the role of LDHA in the brain and its potential as a therapeutic target for brain tumors will also be highlighted. The Warburg effect appears to be universal in tumors, including primary brain tumors, and LDHA (because of its involvement with this process) has been identified as a potential therapeutic target. Currently, there are, however, no suitable LDHA inhibitors available for tumor therapies in the clinic.

Free energy surface of the Michaelis complex of lactate dehydrogenase: a network analysis of microsecond simulations

It has long been recognized that the structure of a protein creates a hierarchy of conformations interconverting on multiple time scales. The conformational heterogeneity of the Michaelis complex is of particular interest in the context of enzymatic catalysis in which the reactant is usually represented by a single conformation of the enzyme/substrate complex. Lactate dehydrogenase (LDH) catalyzes the interconversion of pyruvate and lactate with concomitant interconversion of two forms of the cofactor nicotinamide adenine dinucleotide (NADH and NAD(+)). Recent experimental results suggest that multiple substates exist within the Michaelis complex of LDH, and they show a strong variance in their propensity toward the on-enzyme chemical step. In this study, microsecond-scale all-atom molecular dynamics simulations were performed on LDH to explore the free energy landscape of the Michaelis complex, and network analysis was used to characterize the distribution of the conformations. Our results provide a detailed view of the kinetic network of the Michaelis complex and the structures of the substates at atomistic scales. They also shed light on the complete picture of the catalytic mechanism of LDH.

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.

Understanding interactions of functionalized nanoparticles with proteins: a case study on lactate dehydrogenase

Nanomaterials in biological solutions are known to interact with proteins and have been documented to affect protein function, such as enzyme activity. Understanding the interactions of nanoparticles with biological components at the molecular level will allow for rational designs of nanomaterials for use in medical technologies. Here we present the first detailed molecular mechanics model of functionalized gold nanoparticle (NP) interacting with an enzyme (L-lactate dehydrogenase (LDH) enzyme). Molecular dynamics (MD) simulations of the response of LDH to the NP binding demonstrate that although atomic motions (dynamics) of the main chain exhibit only a minor response to the binding, the dynamics of side chains are significantly constrained in all four active sites that predict alteration in kinetic properties of the enzyme. It is also demonstrated that the 5 nm gold NPs cause a decrease in the maximal velocity of the enzyme reaction (V(max)) and a trend towards a reduced affinity (increased K(m)) for the β-NAD binding site, while pyruvate enzyme kinetics (K(m) and V(max)) are not significantly altered in the presence of the gold NPs. These results demonstrate that modeling of NP:protein interactions can be used to understand alterations in protein function.

Structural characterization of the apo form and NADH binary complex of human lactate dehydrogenase

Lactate dehydrogenase A (LDH-A) is a key enzyme in anaerobic respiration that is predominantly found in skeletal muscle and catalyses the reversible conversion of pyruvate to lactate in the presence of NADH. LDH-A is overexpressed in many tumours and has therefore emerged as an attractive target for anticancer drug discovery. Crystal structures of human LDH-A in the presence of inhibitors have been described, but currently no structures of the apo or binary NADH-bound forms are available for any mammalian LDH-A. Here, the apo structure of human LDH-A was solved at a resolution of 2.1 Å in space group P4122. The active-site loop adopts an open conformation and the packing and crystallization conditions suggest that the crystal form is suitable for soaking experiments. The soaking potential was assessed with the cofactor NADH, which yielded a ligand-bound crystal structure in the absence of any inhibitors. The structures show that NADH binding induces small conformational changes in the active-site loop and an adjacent helix. A comparison with other eukaryotic apo LDH structures reveals the conservation of intra-loop interactions. The structures provide novel insight into cofactor binding and provide the foundation for soaking experiments with fragments and inhibitors.

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

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

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

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

Conformational stability of leukocyte lactate dehydrogenase in healthy men of different age

Comparative analysis of cytosolic fraction distribution and conformational stability in lactate dehydrogenase (LDH) isozymes in healthy men of different age revealed significant differences indicating age-related characteristics of the enzyme. Activity of cytosolic fraction of LDH-1 and conformational stability of LDH-2 and LDH-3 in 40-year-old men were higher than in 20-year-olds. These data suggest age-related determination of conformational stability of individual LDH isozymes and functional heterogeneity of the enzyme within each isoform.

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.