Modeling of Isotope Effects on Binding Oxamate to Lactic Dehydrogenase

Potential Mechanisms of Lactate Dehydrogenase and Bovine Serum Albumin Proteins as Antioxidants: A Mixed Experimental-Computational Study

Proteins have shown varying degrees of antioxidant activity. This study examined the potential mechanisms of interactions between proteins and radicals using chemical kinetics and computational methods. The study quantified the antioxidant activity of lactate dehydrogenase (LDH) and bovine serum albumin (BSA) through Trolox equivalent antioxidant capacity (TEAC) and oxygen radical absorbance capacity (ORAC) assays. BSA was about seven times and LDH 12 times more potent as antioxidants for 2,2'-azino-bis(3-ethylbenzothiazoline-6-sulfonate) (ABTS•-) than they were for peroxyl radicals. According to the evaluation of Trolox equivalents (TE) of 20 proteinogenic amino acids, tryptophan (with a TE value of 101 μmol TE/μmol) exhibited the highest antioxidant activity for ABTS•-, followed by tyrosine (38.7 μmol TE/μmol) and cysteine (30.5 μmol TE/μmol), lysine (0.193 μmol TE/μmol), arginine (0.0325 μmol TE/μmol), valine (0.0280 μmol TE/μmol), histidine (0.00689 μmol TE/μmol), and leucine (0.00560 μmol TE/μmol). The EC50 showed a similar order with a swap between valine and histidine. The antioxidant activity of the amino acids and proteins was temperature dependent. The rate laws, activation energy, and pre-exponential factor A of these reactions provided information on the reaction mechanisms, i.e., a biomolecular elementary step for the reaction of ABTS•- with amino acids tryptophan, tyrosine, cysteine, or protein LDH, and a more complicated mechanism for BSA. The presence of -NH- or hydroxyl groups on aromatic rings enhanced the antioxidant ability of tryptophan and tyrosine. LDH's antioxidant activity did not affect its enzymatic activity, indicating that the radical reaction likely happened on the protein's surface without significantly altering its conformation. The molecular modeling and visualization showed potential reaction sites on the proteins' accessible tryptophan and tyrosine residues. However, the mere surface exposure of tryptophan and tyrosine does not guarantee their antioxidant activities.

The catalytic action of human d-lactate dehydrogenase is severely inhibited by oxalate and is impaired by mutations triggering d-lactate acidosis

d-lactate dehydrogenases are known to be expressed by prokaryotes and by eukaryotic invertebrates, and over the years the functional and structural features of some bacterial representatives of this enzyme ensemble have been investigated quite in detail. Remarkably, a human gene coding for a putative d-lactate dehydrogenase (DLDH) was identified and characterized, disclosing the occurrence of alternative splicing of its primary transcript. This translates into the expression of two human DLDH (hDLDH) isoforms, the molecular mass of which is expected to differ by 2.7 kDa. However, no information on these two hDLDH isoforms is available at the protein level. Here we report on the catalytic action of these enzymes, along with a first analysis of their structural features. In particular, we show that hDLDH is strictly stereospecific, with the larger isoform (hDLDH-1) featuring higher activity at the expense of d-lactate when compared to its smaller counterpart (hDLDH-2). Furthermore, we found that hDLDH is strongly inhibited by oxalate, as indicated by a Ki equal to 1.2 μM for this dicarboxylic acid. Structurally speaking, hDLDH-1 and hDLDH-2 were determined, by means of gel filtration and dynamic light scattering experiments, to be a hexamer and a tetramer, respectively. Moreover, in agreement with previous studies performed with human mitochondria, we identified FAD as the cofactor of hDLDH, and we report here a model of FAD binding by the human d-lactate dehydrogenase. Interestingly, the mutations W323C and T412 M negatively affect the activity of hDLDH, most likely by impairing the enzyme electron-acceptor site.

Activation and friction in enzymatic loop opening and closing dynamics

Protein loop dynamics have recently been recognized as central to enzymatic activity, specificity and stability. However, the factors controlling loop opening and closing kinetics have remained elusive. Here, we combine molecular dynamics simulations with string-method determination of complex reaction coordinates to elucidate the molecular mechanism and rate-limiting step for WPD-loop dynamics in the PTP1B enzyme. While protein conformational dynamics is often represented as diffusive motion hindered by solvent viscosity and internal friction, we demonstrate that loop opening and closing is activated. It is governed by torsional rearrangement around a single loop peptide group and by significant friction caused by backbone adjustments, which can dynamically trap the loop. Considering both torsional barrier and time-dependent friction, our calculated rate constants exhibit very good agreement with experimental measurements, reproducing the change in loop opening kinetics between proteins. Furthermore, we demonstrate the applicability of our results to other enzymatic loops, including the M20 DHFR loop, thereby offering prospects for loop engineering potentially leading to enhanced designs.

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.

Non-Covalent Isotope Effects

In this Perspective, we present examples of isotope effects that originate from noncovalent interactions, mainly hydrogen bonding, electrostatics, and confinement. They are traditionally widely used in isotopic enrichment processes, as well as in studies of mechanisms of different (bio)chemical and physical phenomena. We then show the emerging areas of their applications, mainly medical and material sciences. We stress that these emerging applications require either high enrichment or isotopic substitution, which requires the development of new effective techniques of isotopic purification.

Effects of freezing rate on structural changes in L-lactate dehydrogenase during the freezing process

Freezing is a common method for improving enzyme storage stability. During the freezing process, the freezing rate is an important parameter that can affect protein stability. However, there is limited information on the denaturation mechanisms and protein conformational changes associated with the freezing rate. In this study, the effects of freezing rate on activity loss and conformational changes in a model enzyme, l-lactate dehydrogenase, were evaluated. Enzyme solutions were frozen at various rates, from 0.2 to 70.6 °C/min, and ice seeding was conducted to reduce supercooling. The results demonstrated that fast freezing results in activity loss, structural changes, and aggregation. The residual activities at freezing rates of 0.2, 12.8, and 70.6 °C/min were 77.6 ± 0.9%, 64.1 ± 0.4%, and 44.8 ± 2.0%, respectively. As the freezing rate increased, the degree of dissociation and unfolding increased significantly, as determined using blue native-polyacrylamide gel electrophoresis and fluorescence spectroscopy. Moreover, a large number of amyloid aggregates were detected in samples frozen at a fast freezing rate (70.6 °C/min). The enzyme inactivation mechanism induced by fast freezing was proposed in terms of increased dehydration at the enzyme surface and an ice/unfroze solution interface, which could be helpful to establish a common understanding of enzyme inactivation during the freezing process.

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.

Predicting protein-ligand binding affinity and correcting crystal structures with quantum mechanical calculations: lactate dehydrogenase A

Accurately computing the geometry and energy of host-guest and protein-ligand interactions requires a physically accurate description of the forces in action. Quantum mechanics can provide this accuracy but the calculations can require a prohibitive quantity of computational resources. The size of the calculations can be reduced by including only the atoms of the receptor that are in close proximity to the ligand. We show that when combined with log P values for the ligand (which can be computed easily) this approach can significantly improve the agreement between computed and measured binding energies. When the approach is applied to lactate dehydrogenase A, it can make quantitative predictions about conformational, tautomeric and protonation state preferences as well as stereoselectivity and even identifies potential errors in structures deposited in the Protein Data Bank for this enzyme. By broadening the evidence base for these structures from only the diffraction data, more chemically realistic structures can be proposed.

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.

The self-inhibitory nature of metabolic networks and its alleviation through compartmentalization

Metabolites can inhibit the enzymes that generate them. To explore the general nature of metabolic self-inhibition, we surveyed enzymological data accrued from a century of experimentation and generated a genome-scale enzyme-inhibition network. Enzyme inhibition is often driven by essential metabolites, affects the majority of biochemical processes, and is executed by a structured network whose topological organization is reflecting chemical similarities that exist between metabolites. Most inhibitory interactions are competitive, emerge in the close neighbourhood of the inhibited enzymes, and result from structural similarities between substrate and inhibitors. Structural constraints also explain one-third of allosteric inhibitors, a finding rationalized by crystallographic analysis of allosterically inhibited L-lactate dehydrogenase. Our findings suggest that the primary cause of metabolic enzyme inhibition is not the evolution of regulatory metabolite-enzyme interactions, but a finite structural diversity prevalent within the metabolome. In eukaryotes, compartmentalization minimizes inevitable enzyme inhibition and alleviates constraints that self-inhibition places on metabolism.

Ligand binding phenomena that pertain to the metabolic function of renalase

Renalase catalyzes the oxidation of isomers of β-NAD(P)H that carry the hydride in the 2 or 6 positions of the nicotinamide base to form β-NAD(P)⁺. This activity is thought to alleviate inhibition of multiple β-NAD(P)-dependent enzymes of primary and secondary metabolism by these isomers. Here we present evidence for a variety of ligand binding phenomena relevant to the function of renalase. We offer evidence of the potential for primary metabolism inhibition with structures of malate dehydrogenase and lactate dehydrogenase bound to the 6-dihydroNAD isomer. The previously observed preference of renalase from Pseudomonas for NAD-derived substrates over those derived from NADP is accounted for by the structure of the enzyme in complex with NADPH. We also show that nicotinamide nucleosides and mononucleotides reduced in the 2- and 6-positions are renalase substrates, but bind weakly. A seven-fold enhancement of acquisition (k_red/K_d) for 6-dihydronicotinamide riboside was observed for human renalase in the presence of ADP. However, generally the addition of complement ligands, AMP for mononucleotide or ADP for nucleoside substrates, did not enhance the reductive half-reaction. Non-substrate nicotinamide nucleosides or nucleotides bind weakly suggesting that only β-NADH and β-NADPH compete with dinucleotide substrates for access to the active site.

Impact of cyclometalated ruthenium(II) complexes on lactate dehydrogenase activity and cytotoxicity in gastric and colon cancer cells

Lactate dehydrogenase (LDH) is a redox enzyme often overexpressed in cancer cells allowing their survival in stressful metabolic tumor environment. Ruthenium(II) complexes have been shown to impact on the activity of purified horseradish peroxidase and glucose oxidase but the physiological relevance remains unclear. In this study we investigated how ruthenium complexes impact on the activity of LDH in vitro and in cancer cells and performed a comparative study using polypyridine ruthenium(II) complex [Ru(bpy)₃]²⁺ (1) and its structurally related cyclometalated 2-phenylpyridinato counterpart [Ru(phpy)(bpy)₂]⁺ (2) (bpy=2,2'-bipyridine, phpyH=2-phenylpyridine). We show that the cytotoxicity in gastric and colon cancer cells induced by 2 is significantly higher compared to 1. The kinetic inhibition mechanisms on purified LDH and the corresponding inhibition constants K_i or i_0.5 values were calculated. Though complexes 1 and 2 are structurally very similar (one Ru-C bond in 2 replaces one Ru-N bond in 1), their inhibition modes are different. Cyclometalated complex 2 behaves exclusively as a non-competitive inhibitor of LDH from rabbit muscle (LDH_rm), strongly suggesting that 2 does not interact with LDH in the vicinities of either lactate/pyruvate or NAD⁺/NADH binding sites. Sites of interaction of 1 and 2 with LDH_rm were revealed theoretically through computational molecular docking. Inhibition of LDH activity by 2 was confirmed in cancer cells. Altogether, these results revealed an inhibition of LDH activity by ruthenium complex through a direct interaction structurally tuned by a Ru-C bond.

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.

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.

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.

Binding ligands and cofactor to L-lactate dehydrogenase from human skeletal and heart muscles

Binding affinities of cofactor and ligands to the active site of two different isoforms of lactate dehydrogenase (LDH), heart and skeletal muscles (H4 and M4, respectively), can be used for medical and biological applications. Herein, a hybrid QM/MM computational approach based on free energy perturbation methods has been carried out to estimate binding affinities and binding isotope effects (BIEs) for NADH/NAD(+) and oxamate, pyruvate, L-lactate, and D-lactate ligands to the M4 and H4 isoforms of L-LDH. Here, we show that determining how cofactor and ligands interact with the active site of LDH isoforms advanced the still open discussion on the intracellular lactate shuttle hypothesis. In our discussion we deny the key concept of this hypothesis showing, based on interaction energy values, that there is no evidence that the M4 type of LDH in the skeletal muscles cells served as a catalyst of the conversion of lactate to pyruvate. Additionally, theoretical determination of BIEs for H4 and M4 types of LDH shows that there is a way of using the BIEs as a tool capable to distinguish these isoforms, and for this purpose D-lactate labeled with deuterium in positions 11 or 7, 8, 9 ([11-2H]-BIE and [7,8,9-2H3]-BIE) or L-lactate labeled only in position 11 ([11-2H]-BIE) could be used. We propose the BIEs as a useful tool which can be applied in order to experimentally determine the types of LDH.

Discovery of N-hydroxyindole-based inhibitors of human lactate dehydrogenase isoform A (LDH-A) as starvation agents against cancer cells

Highly invasive tumor cells are characterized by a metabolic switch, known as the Warburg effect, from "normal" oxidative phosphorylation to increased glycolysis even under sufficiently oxygenated conditions. This dependence on glycolysis also confers a growth advantage to cells present in hypoxic regions of the tumor. One of the key enzymes involved in glycolysis, the muscle isoform of lactate dehydrogenase (LDH-A), is overexpressed by metastatic cancer cells and is linked to the vitality of tumors in hypoxia. This enzyme may be considered as a potential target for new anticancer agents, since its inhibition cuts cancer energetic and anabolic supply, thus reducing the metastatic and invasive potential of cancer cells. We have discovered new and efficient N-hydroxyindole-based inhibitors of LDH-A, which are isoform-selective (over LDH-B) and competitive with both the substrate (pyruvate) and the cofactor (NADH). The antiproliferative activity of these compounds was confirmed on a series of cancer cell lines, and they proved to be particularly effective under hypoxic conditions. Moreover, NMR experiments showed that these compounds are able to reduce the glucose-to-lactate conversion inside the cell.

A new scheme to calculate isotope effects

We present a new scheme to calculate isotope effects. Only selected frequencies at the target level of theory are calculated. The frequencies are selected by an analysis of the Hessian from a lower level of theory. We obtain accurate isotope effects without calculating the full Hessian at the target level of theory. The calculated frequencies are very accurate. The scheme converges to the correct isotope effect.

Importance of the lactate dehydrogenase quaternary structure in theoretical calculations

Using the example of lactate dehydrogenase, we show that enzyme quaternary structure has an important influence on the structure of the active site and that models that comprise all amino acids in the vicinity of an active site, but are missing this structural information, can lead to incorrect results. We also show that binding isotope effects are very sensitive to the geometric parameters, and thus one should be very cautious when interpreting results obtained with models that are too coarse. In terms of the type of hydrogen bonds, our results indicate that binding isotope effects are pronounced only when a hydrogen bond exhibits some covalent character.