Malate dehydrogenase of the cytosol. A kinetic investigation of the reaction mechanism and a comparison with lactate dehydrogenase

Catalytic mechanism and kinetics of malate dehydrogenase

Malate dehydrogenase (MDH) is a ubiquitous and central enzyme in cellular metabolism, found in all kingdoms of life, where it plays vital roles in the cytoplasm and various organelles. It catalyzes the reversible NAD+-dependent reduction of L-malate to oxaloacetate. This review describes the reaction mechanism for MDH and the effects of mutations in and around the active site on catalytic activity and substrate specificity, with a particular focus on the loop that encloses the active site after the substrates have bound. While MDH exhibits selectivity for its preferred substrates, mutations can alter the specificity of MDH for each cosubstrate. The kinetic characteristics and similarities of a variety of MDH isozymes are summarized, and they illustrate that the KM values are consistent with the relative concentrations of the substrates in cells. As a result of its existence in different cellular environments, MDH properties vary, making it an attractive model enzyme for studying enzyme activity and structure under different conditions.

Molecular Characteristics of the Malate Dehydrogenase (MDH) Gene Family in Spirometra mansoni (Cestoda: Diphyllobothriidea)

The plerocercoid larva of Spirometra mansoni can cause a parasitic zoonosis-sparganosis. Malate dehydrogenase (MDH) plays a very important role in the life activities of parasites. However, little is known about the MDH family in S. mansoni. We identified eight new MDH members in S. mansoni in this study. Clustering analysis divided SmMDHs into two groups and revealed patterns similar to the conserved motif organization. RT-qPCR suggested that five MDHs were highly expressed in the mature proglottid and that three MDHs were highly expressed in the gravid proglottid. Phylogenetic analysis revealed that SmMDHs contain both conserved family members and members in the process of further diversification. rSmMDH has an NAD binding domain, a dimer interface and a substrate binding domain. Natural SmMDH was immunolocalized in the tissues and follicles around the uterus in the mature or gravid proglottid and eggshells. The maximum forward and reverse reaction activities of rSmMDH were observed at pH 8.5 and 9.0, respectively. The optimum temperature for enzyme activity was 37 °C in the forward reaction and 40 °C in the reverse reaction. These results lay the foundation for studying the molecular functions and mechanisms of MDHs in S. mansoni and related taxa.

Characterization of the kinetics of cardiac cytosolic malate dehydrogenase and comparative analysis of cytosolic and mitochondrial isoforms

Because the mitochondrial inner membrane is impermeable to pyridine nucleotides, transport of reducing equivalents between the mitochondrial matrix and the cytoplasm relies on shuttle mechanisms, including the malate-aspartate shuttle and the glycerol-3-phosphate shuttle. These shuttles are needed for reducing equivalents generated by metabolic reactions in the cytosol to be oxidized via aerobic metabolism. Two isoenzymes of malate dehydrogenase (MDH) operate as components of the malate-aspartate shuttle, in which a reducing equivalent is transported via malate, which when oxidized to oxaloacetate, transfers an electron pair to reduce NAD to NADH. Several competing mechanisms have been proposed for the MDH-catalyzed reaction. This study aims to identify the pH-dependent kinetic mechanism for cytoplasmic MDH (cMDH) catalyzed oxidation/reduction of MAL/OAA. Experiments were conducted assaying the forward and reverse directions with products initially present, varying pH between 6.5 and 9.0. By fitting time-course data to various mechanisms, it is determined that an ordered bi-bi mechanism with coenzyme binding first followed by the binding of substrate is able to explain the kinetic data. The proposed mechanism is similar to, but not identical to, the mechanism recently determined for the mitochondrial isoform, mMDH. cMDH and mMDH mechanisms are also shown to both be reduced versions of a common, more complex mechanism that can explain the kinetic data for both isoforms. Comparing the simulated activity (ratio of initial velocity to the enzyme concentration) under physiological conditions, the mitochondrial MDH (mMDH) activity is predicted to be higher than cMDH activity under mitochondrial matrix conditions while the cMDH activity is higher than mMDH activity under cytoplasmic conditions, suggesting that the functions of the isoforms are kinetically tuned to their individual physiological roles.

Intrinsic noise analyzer: a software package for the exploration of stochastic biochemical kinetics using the system size expansion

The accepted stochastic descriptions of biochemical dynamics under well-mixed conditions are given by the Chemical Master Equation and the Stochastic Simulation Algorithm, which are equivalent. The latter is a Monte-Carlo method, which, despite enjoying broad availability in a large number of existing software packages, is computationally expensive due to the huge amounts of ensemble averaging required for obtaining accurate statistical information. The former is a set of coupled differential-difference equations for the probability of the system being in any one of the possible mesoscopic states; these equations are typically computationally intractable because of the inherently large state space. Here we introduce the software package intrinsic Noise Analyzer (iNA), which allows for systematic analysis of stochastic biochemical kinetics by means of van Kampen's system size expansion of the Chemical Master Equation. iNA is platform independent and supports the popular SBML format natively. The present implementation is the first to adopt a complementary approach that combines state-of-the-art analysis tools using the computer algebra system Ginac with traditional methods of stochastic simulation. iNA integrates two approximation methods based on the system size expansion, the Linear Noise Approximation and effective mesoscopic rate equations, which to-date have not been available to non-expert users, into an easy-to-use graphical user interface. In particular, the present methods allow for quick approximate analysis of time-dependent mean concentrations, variances, covariances and correlations coefficients, which typically outperforms stochastic simulations. These analytical tools are complemented by automated multi-core stochastic simulations with direct statistical evaluation and visualization. We showcase iNA's performance by using it to explore the stochastic properties of cooperative and non-cooperative enzyme kinetics and a gene network associated with circadian rhythms. The software iNA is freely available as executable binaries for Linux, MacOSX and Microsoft Windows, as well as the full source code under an open source license.

Communication: limitations of the stochastic quasi-steady-state approximation in open biochemical reaction networks

It is commonly believed that, whenever timescale separation holds, the predictions of reduced chemical master equations obtained using the stochastic quasi-steady-state approximation are in very good agreement with the predictions of the full master equations. We use the linear noise approximation to obtain a simple formula for the relative error between the predictions of the two master equations for the Michaelis-Menten reaction with substrate input. The reduced approach is predicted to overestimate the variance of the substrate concentration fluctuations by as much as 30%. The theoretical results are validated by stochastic simulations using experimental parameter values for enzymes involved in proteolysis, gluconeogenesis, and fermentation.

Spatially resolved fluorescence lifetime mapping of enzyme kinetics in living cells

Traditional cuvette-based enzyme studies lack spatial information and do not allow real-time monitoring of the effects of modulating enzyme functions in vivo. In order to probe the realistic timescales of steric modifications in enzyme-substrate complexes and functional binding-unbinding kinetics in living cells without losing spatial information, it is imperative to develop sensitive imaging strategies that can report enzyme kinetics in real time over a wide dynamic range of timescales. Here we present a multi-photon excitation-based, ultra-fast photon detection using a streak camera and Laguerre expansion-based fast deconvolution approach for achieving high spatio-temporal resolution in monitoring real-time enzyme kinetics in single cells. In particular, we report spatially resolved, nanosecond-scale fluorescence dynamics associated with binding-unbinding kinetics of endogenous metabolic co-factor nicotinamide adenine dinucleotide with enzymes in intact living cells. By monitoring real-time kinetics of NAD(P)H-enzyme kinetics in primary hepatocytes isolated from young and aged mouse models, we observed that the mechanism of inhibition of mitochondrial respiration at complex I site is mediated by redistribution of free and protein-bound nicotinamide adenine dinucleotide pools and that this equilibrium redistribution is affected by age-related modifications in mitochondrial function. We describe unique advantages of Laguerre deconvolution algorithm in comparison with conventional lifetime analysis approaches. Non-invasive monitoring of metabolic dysfunctions in intact animal models is an attractive strategy for gaining insight into the dynamics of tissue metabolism in health and in various metabolic syndromes such as cancer, diabetes and aging-induced metabolic dysfunctions. Besides the example demonstrated above, we envisage that the proposed method can find applications in a variety of other situations where intensity-based approaches fall short owing to spectroscopic artefacts.

Simulation of the enzyme reaction mechanism of malate dehydrogenase

A hybrid numerical method, which employs molecular mechanics to describe the bulk of the solvent-protein matrix and a semiempirical quantum-mechanical treatment for atoms near the reactive site, was utilized to simulate the minimum energy surface and reaction pathway for the interconversion of malate and oxaloacetate catalyzed by the enzyme malate dehydrogenase (MDH). A reaction mechanism for proton and hydride transfers associated with MDH and cofactor nicotinamide adenine dinucleotide (NAD) is deduced from the topology of the calculated energy surface. The proposed mechanism consists of (1) a sequential reaction with proton transfer preceding hydride transfer (malate to oxaloacetate direction), (2) the existence of two transition states with energy barriers of approximately 7 and 15 kcal/mol for the proton and hydride transfers, respectively, and (3) reactant (malate) and product (oxaloacetate) states that are nearly isoenergetic. Simulation analysis of the calculated energy profile shows that solvent effects due to the protein matrix dramatically alter the intrinsic reactivity of the functional groups involved in the MDH reaction, resulting in energetics similar to that found in aqueous solution. An energy decomposition analysis indicates that specific MDH residues (Arg-81, Arg-87, Asn-119, Asp-150, and Arg-153) in the vicinity of the substrate make significant energetic contributions to the stabilization of proton transfer and destabilization of hydride transfer. This suggests that these amino acids play an important role in the catalytic properties of MDH.

D-2-hydroxy-4-methylvalerate dehydrogenase from Lactobacillus delbrueckii subsp. bulgaricus. I. Kinetic mechanism and pH dependence of kinetic parameters, coenzyme binding and substrate inhibition

The steady-state kinetics of D-2-hydroxy-4-methylvalerate dehydrogenase have been studied at pH 8.0 by initial velocity, product inhibition, and dead-end inhibition techniques. The mechanism is rapid-equilibrium ordered in the NAD+ plus D-2-hydroxy-4-methylvalerate direction, and steady-state ordered in the other direction. In both cases coenzyme is the first substrate added and both the E-NADH-D-2-hydroxy-4-methylvalerate and E-NAD+-2-oxo-4-methylvalerate give rise to abortive complexes which cause excess substrate inhibition. Steady-state measurements show that the rate-limiting step in both directions at pH 8.0 is between formation of the enzyme-coenzyme-substrate ternary complex and the release of the first product of the reaction. Transient kinetics combined with primary kinetic deuterium isotope effects show that in the NADH-->NAD+ direction there is a slow, rate-limiting rearrangement of the E-NADH-oxoacid complex while hydride transfer is very fast. The release of NAD+ at pH 8.0 is 200-times faster than Kcat (NADH-->NAD+) whereas the release of NADH is only 5-times faster than Kcat (NAD+-->NADH). The pH dependence of NADH binding depends upon the presence of two ionizable residues with a pKa of about 5.9. The pH dependence of kinetic parameters is explained by a third ionizable residue with pKa values 7.2 (in the E-NADH complex) and < or = 6.4 (in the E-NAD+ complex) which may be the proton donor and acceptor for the chemical reaction. At pH 6.5 the mechanism changes in the NADH-->NAD+ direction to be partly limited by the chemical step with a measured primary kinetic isotope effect of 5.7 and partly by an only slightly faster dissociation of NAD+. In addition the inhibition by excess oxo-4-methylvalerate is more pronounced. The mechanism implies that removing the positive charges created by the two groups which control coenzyme affinity could both enhance the catalytic rate at pH 6.5 and diminish excess substrate inhibition to provide an enzyme better suited to the bulk synthesis of D-2-hydroxyacids.

Catalytic-rate improvement of a thermostable malate dehydrogenase by a subtle alteration in cofactor binding

The nucleotide-binding fold of many NAD(+)-dependent dehydrogenases contains a conserved acidic amino acid residue which hydrogen-bonds with the 2'- and 3'-hydroxy groups of the adenine-ribose of the cofactor. This residue is highly conserved as aspartate in malate dehydrogenases, except in the thermophilic enzyme from Thermus aquaticus B (TaqMDH), which has glutamic acid-41 in the equivalent position. The catalytic mechanism was dissected to investigate the functional significance of this difference in TaqMDH with respect to a mutant enzyme where glutamic acid-41 was replaced by aspartic acid. The mutant enzyme was found to retain a high degree of protein structural stability to both thermal and chemical denaturation. When compared with the wild-type enzyme the mutant had a higher Km and Kd for both reduced and oxidized cofactors (NADH and NAD+) and a 2-3-fold increase in steady-state kcat in both assay directions. The rate-determining step for the reduction of oxaloacetate by wild-type TaqMDH was shown to be the rate of NAD+ release, which was about 2.5-fold higher for the mutant enzyme. This correlates well with the 1.8-fold higher steady-state kcat of the mutant enzyme and represents an improvement in the steady-state kcat of a thermophilic enzyme at moderate temperature by a conservative amino acid substitution which increases the rate of product release.

Determinants of coenzyme specificity in glyceraldehyde-3-phosphate dehydrogenase: role of the acidic residue in the fingerprint region of the nucleotide binding fold

On the basis of the three-dimensional structure of the glycolytic NAD-dependent glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and of sequence comparison with the photosynthetic NAD(P)-dependent GAPDH of the chloroplast, a series of mutants of GAPDH from Bacillus stearothermophilus have been constructed. The results deduced from kinetic and binding studies suggest that the absence of activity of the wild-type GAPDH with NADP as a cofactor is the consequence of at least three factors: (1) steric hindrance, (2) electrostatic repulsion between the charged carboxyl group of Asp32 and the 2'PO4, and (3) structural determinants at the subunit interface of the tetramer. The best value for kcat/KM and KD for NADP was observed for the D32A-L187A-P188S mutant. This triple mutation leads to a switch in favor of NADP specificity but with a kcat/KM ratio 50- and 80-fold less than that observed for the wild type with NAD and for the chloroplast GAPDH with NADP, respectively. Substituting the invariant chloroplastic Thr33-Gly34-Gly35 for the B. stearothermophilus Leu33-Thr34-Asp35 residues on the double mutant Ala187-Ser188 does not improve significantly the affinity for NADP while substituting Ala32 for Asp32 on the double mutant does. Clearly, other subtle adjustments in the adenosine subsite are needed to reconcile the presence of the carboxylate group of Asp32 and the 2'-phosphate of NADP. Kinetic studies indicate a change of the rate-limiting step for the mutants. This could be the consequence of an incomplete apo-holo transition.(ABSTRACT TRUNCATED AT 250 WORDS)

A specific, highly active malate dehydrogenase by redesign of a lactate dehydrogenase framework

Three variations to the structure of the nicotinamide adenine dinucleotide (NAD)-dependent L-lactate dehydrogenase from Bacillus stearothermophilus were made to try to change the substrate specificity from lactate to malate: Asp197----Asn, Thr246----Gly, and Gln102----Arg). Each modification shifts the specificity from lactate to malate, although only the last (Gln102----Arg) provides an effective and highly specific catalyst for the new substrate. This synthetic enzyme has a ratio of catalytic rate (kcat) to Michaelis constant (Km) for oxaloacetate of 4.2 x 10(6)M-1 s-1, equal to that of native lactate dehydrogenase for its natural substrate, pyruvate, and a maximum velocity (250 s-1), which is double that reported for a natural malate dehydrogenase from B. stearothermophilus.

Lactate and malate dehydrogenase binding to the microsomal fraction from chicken liver

Chicken liver microsomal fractions show lactate and malate dehydrogenase activities which behave differently with respect to successive extractions by sonication in 0.15 M NaCl, 0.2% Triton X-100 and 0.15 M NaCl, respectively. The Triton X-100-treated pellet did not show malate dehydrogenase activity but exhibited a 10-fold increase in lactate dehydrogenase activity with respect to the sonicated pellet. Total extracted lactate and malate dehydrogenase activities were, respectively, 7.5 and 1.7 times higher than that in the initial pellet. Different isoenzyme compositions were observed for cytosoluble and microsomal extracted lactate and malate dehydrogenases. When the ionic strength (0-500 mM) or the pH values (6.1-8.7) of the media were increased, an efficient release of lactate dehydrogenase was found at NaCl 30-70 mM and pH 6.6-7.3. Malate dehydrogenase solubilization under the same conditions was very small, even at NaCl 500 mM, but it attained a maximum in the 7.3-8.7 pH range. Cytosoluble lactate dehydrogenase bound in vitro to 0.15 M NaCl-treated (M2) and sonicated (M3) microsomal fractions but not to the crude microsomal fraction (M1). Particle saturation by lactate dehydrogenase occurred with M2 and M3, which contained binding sites with different affinities. Cytosoluble malate dehydrogenase did not bind to M1, M2 and M3 fractions, however, a little binding was found when purified basic malate dehydrogenase was incubated with M2 or M3 fractions.

Site-directed mutagenesis reveals role of mobile arginine residue in lactate dehydrogenase catalysis

The binding of substrates to lactate dehydrogenases induces a marked rearrangement of the protein structure in which a 'loop' of polypeptide (residues 98-110) closes over the active site of the enzyme. In this rearrangement, arginine 109 (a basic residue conserved in all known lactate dehydrogenase sequences and in the homologous malate dehydrogenases) moves 0.8 nm from a position in the solvent to one in the active site where its guanidinium group resides within hydrogen bonding distance of both the reactive carbonyl of pyruvate and imidazole ring of the catalytic histidine 195 (see Fig. 1). Whilst this feature of the enzyme has been commented upon previously, the function of this mobile arginine residue during catalysis has not been tested experimentally. The advent of protein engineering has now enabled us to define the role of this basic residue by substituting it with the neutral glutamine. Transient kinetic and equilibrium studies of the mutant enzyme indicate that arginine 109 enhances the polarization of the pyruvate carbonyl group in the ground state and stabilizes the transition state. The gross active-site structure of the enzyme is not altered by the mutation since an alternative catalytic function of the enzyme (rate of addition of sulphite to NAD+), which does not require hydride transfer, is insensitive to the arginine----glutamine substitution.

Computer simulation of metabolism in palmitate-perfused rat heart. II. Behavior of complete model

Intermediary metabolism in rat hearts perfused with 11 mM glucose plus 1 mM palmitate was simulated by a computer model. Several enzyme submodels in a previous version of the isolated rat heart computer model were improved, and a new fatty acid oxidation pathway model was added. Compartmentation of metabolites in a pseudo-stationary state was calculated, and its implications are discussed, e.g., citrate level may not regulate glycolysis because it is mostly mitochondrial. Citrate synthetase, controlled largely by its inhibitors, is of key importance in regulating fatty acid metabolism. The response of aconitase to the mitochondrial Mg2+ level is of major importance in setting both the mitochondrial citrate and isocitrate levels. Pyruvate dehydrogenase is about 96% in the inactive phosphorylated form, and the active form is also 15% inhibited by products, severely limiting pyruvate oxidation and causing preferential utilization of palmitate as the metabolic fuel. The simulation is consistent with a creatine phosphate shuttle which delivers high energy phosphate to the site of its utilization for mechanical work.

Kinetic mechanism of the molecular forms of chicken liver mitochondrial malate dehydrogenase

1. The reaction kinetic mechanism (pH 7.4) of the molecular forms of chicken liver m-MDH is of the order bi-bi ternary complex type with the existence of the E-oxaloacetate, E-L-malate, E-NAD+-oxaloacetate, E-NADH-L-malate, E-NAD+-NADH, E-NAD+-NAD+, E-NADH-NAD+ and E-NAD-NADH abortive complexes. 2. The saturating concentration values of the substrates are notably modified, in certain cases, in the presence of the reaction products.