The contribution of electrostatic and van der Waals interactions to the stereospecificity of the reaction catalyzed by lactate dehydrogenase

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.

Molecular dynamics studies of modular polyketide synthase ketoreductase stereospecificity

Ketoreductases (KRs) from modular polyketide synthases (PKSs) can perform stereospecific catalysis, selecting a polyketide with a D- or L-α-methyl substituent for NADPH-mediated reduction. In this report, molecular dynamics (MD) simulations were performed to investigate the interactions that control stereospecificity. We studied the A1-type KR from the second module of the amphotericin PKS (A1), which is known to be stereospecific for a D-α-methyl-substituted diketide substrate (dkD). MD simulations of two ternary complexes comprised of the enzyme, NADPH, and either the correct substrate, dkD, or its enantiomer (dkL) were performed. The coordinates for the A1/NADPH binary complex were obtained from a crystal structure (PDB entry 3MJS), and substrates were modeled in the binding pocket in conformations appropriate for reduction. Simulations were intended to reproduce the initial weak binding of the polyketide substrate to the enzyme. Long (tens of nanoseconds) MD simulations show that the correct substrate is retained in a conformation closer to the reactive configuration. Many short (up to a nanosecond) MD runs starting from the initial structures display evidence that Q364, three residues N-terminal to the catalytic tyrosine, forms a hydrogen bond to the incorrect dkL substrate to yield an unreactive conformation that is more favorable than the reactive configuration. This interaction is not as strong for dkD, as the D-α-methyl substituent is positioned between the glutamine and the reactive site. This result correlates with experimental findings [Zheng, J., et al. (2010) Structure 18, 913-922] in which a Q364H mutant was observed to lose stereospecificity.

Steric effect for proton, hydrogen-atom, and hydride transfer reactions with geometric isomers of NADH-model ruthenium complexes

Two isomers, [Ru(1)]2+ (Ru = Ru(bpy)2, bpy = 2,2'-bipyridine, 1 = 2-(pyrid-2'-yl)-1-azaacridine) and [Ru(2)]2+ (2 = 3-(pyrid-2'-yl)-4-azaacridine), are bioinspired model compounds containing the nicotinamide functionality and can serve as precursors for the photogeneration of C-H hydrides for studying reactions pertinent to the photochemical reduction of metal-C1 complexes and/or carbon dioxide. While it has been shown that the structural differences between the azaacridine ligands of [Ru(1)]2+ and [Ru(2)]2+ have a significant effect on the mechanism of formation of the hydride donors, [Ru(1HH)]2+ and [Ru(2HH)]2+, in aqueous solution, we describe the steric implications for proton, net-hydrogen-atom and net-hydride transfer reactions in this work. Protonation of [Ru(2*-)] in aprotic and even protic media is slow compared to that of [Ru(1*-)]+. The net hydrogen-atom transfer between *[Ru(1)]2+ and hydroquinone (H2Q) proceeds by one-step EPT, rather than stepwise electron-proton transfer. Such a reaction was not observed for *[Ru(2)]2+ because the non-coordinated N atom is not easily available for an interaction with H2Q. Finally, the rate of the net hydride ion transfer from [Ru(1HH)]2+ to [Ph3C]+ is significantly slower than that of [Ru (2HH)]2+ owing to steric congestion at the donor site.

Modeling electrostatic effects in proteins

Electrostatic energies provide what is perhaps the most effective tool for structure-function correlation of biological molecules. This review considers the current state of simulations of electrostatic energies in macromolecules as well as the early developments of this field. We focus on the relationship between microscopic and macroscopic models, considering the convergence problems of the microscopic models and the fact that the dielectric 'constants' in semimacroscopic models depend on the definition and the specific treatment. The advances and the challenges in the field are illustrated considering a wide range of functional properties including pK(a)'s, redox potentials, ion and proton channels, enzyme catalysis, ligand binding and protein stability. We conclude by pointing out that, despite the current problems and the significant misunderstandings in the field, there is an overall progress that should lead eventually to quantitative descriptions of electrostatic effects in proteins and thus to quantitative descriptions of the function of proteins.

Dynamics of protein ligand binding on multiple time scales: NADH binding to lactate dehydrogenase

Although the importance of atomic motion to how proteins function has been conjectured for several decades, the characterization of protein dynamics on multiple time scales is scant. This is because of severe experimental and theoretical difficulties, particularly characterizing the nanosecond to millisecond time scales. Here, we apply advanced laser-induced temperature-jump relaxation spectroscopic techniques to examine the kinetics of NADH binding to lactate dehydrogenase over this time scale. The bimolecular rate process, at about 290 micros, is easily observed as are multiple faster events (with relaxation times of 200 ns, 3.5 micros, and 24 micros), revealing a rich dynamical nature of the binding step. The results show that there are multiple structures of bound enzyme-ligand complexes, some of which are likely to be far from the catalytically productive structure. The results have important implications for interpretations of the binding thermodynamics of ligands to LDH and, by extension, to other proteins. The observed processes likely play a role in the dynamics of the chemistry that is catalyzed by lactate dehydrogenase.

Hot spots in cold adaptation: localized increases in conformational flexibility in lactate dehydrogenase A4 orthologs of Antarctic notothenioid fishes

To elucidate mechanisms of enzymatic adaptation to extreme cold, we determined kinetic properties, thermal stabilities, and deduced amino acid sequences of lactate dehydrogenase A4 (A4-LDH) from nine Antarctic (-1.86 to 1 degree C) and three South American (4 to 10 degree C) notothenioid teleosts. Higher Michaelis-Menten constants (Km) and catalytic rate constants (kcat) distinguish orthologs of Antarctic from those of South American species, but no relationship exists between adaptation temperature and the rate at which activity is lost because of heat denaturation. In all species, active site residues are conserved fully, and differences in kcat and Km are caused by substitutions elsewhere in the molecule. Within geographic groups, identical kinetic properties are generated by different substitutions. By combining our data with A4-LDH sequences for other vertebrates and information on roles played by localized conformational changes in setting kcat, we conclude that notothenioid A4-LDHs have adapted to cold temperatures by increases in flexibility in small areas of the molecule that affect the mobility of adjacent active-site structures. Using these findings, we propose a model that explains linked temperature-adaptive variation in Km and kcat. Changes in sequence that increase flexibility of regions of the enzyme involved in catalytic conformational changes may reduce energy (enthalpy) barriers to these rate-governing shifts in conformation and, thereby, increase kcat. However, at a common temperature of measurement, the higher configurational entropy of a cold-adapted enzyme may foster conformations that bind ligands poorly, leading to high Km values relative to warm-adapted orthologs.

Electrostatic effects in macromolecules: fundamental concepts and practical modeling

The past few years have seen an exponential growth in the calculations of electrostatic energies of macromolecules and an increased recognition of the crucial role of electrostatic effects. This review considers the current state of the field. Focus is placed on calculations of pKas, redox potentials and binding energies in macromolecules and clarification of the fact that the value of the dielectric 'constant' of a protein depends on its definition and that small dielectric constants should not be used in describing charge-charge interactions by current continuum models.