The importance of arginine 102 for the substrate specificity of Escherichia coli malate dehydrogenase

The structural biology and dynamics of malate dehydrogenases

Malate dehydrogenase (MDH) enzymes catalyze the reversible oxidoreduction of malate to oxaloacetate using NAD(P) as a cofactor. This reaction is vital for metabolism and the exchange of reducing equivalents between cellular compartments. There are more than 100 structures of MDH in the Protein Data Bank, representing species from archaea, bacteria, and eukaryotes. This conserved family of enzymes shares a common nucleotide-binding domain, substrate-binding domain, and subunits associate to form a dimeric or a tetrameric enzyme. Despite the variety of crystallization conditions and ligands in the experimental structures, the conformation and configuration of MDH are similar. The quaternary structure and active site dynamics account for most conformational differences in the experimental MDH structures. Oligomerization appears essential for activity despite each subunit having a structurally independent active site. There are two dynamic regions within the active site that influence substrate binding and possibly catalysis, with one of these regions adjoining the subunit interface. In this review, we introduce the reader to the general structural framework of MDH highlighting the conservation of certain features and pointing out unique differences that regulate MDH enzyme activity.

Malate dehydrogenase: a story of diverse evolutionary radiation

Malate dehydrogenase (MDH) is a ubiquitous enzyme involved in cellular respiration across all domains of life. MDH's ubiquity allows it to act as an excellent model for considering the history of life and how the rise of aerobic respiration and eukaryogenesis influenced this evolutionary process. Here, we present the diversity of the MDH family of enzymes across bacteria, archaea, and eukarya, the relationship between MDH and lactate dehydrogenase (LDH) in the formation of a protein superfamily, and the connections between MDH and endosymbiosis in the formation of mitochondria and chloroplasts. The development of novel and powerful DNA sequencing techniques has challenged some of the conventional wisdom underlying MDH evolution and suggests a history dominated by gene duplication, horizontal gene transfer, and cryptic endosymbiosis events and adaptation to a diverse range of environments across all domains of life over evolutionary time. The data also suggest a superfamily of proteins that do not share high levels of sequential similarity but yet retain strong conservation of core function via key amino acid residues and secondary structural components. As DNA sequencing and 'big data' analysis techniques continue to improve in the life sciences, it is likely that the story of MDH will continue to refine as more examples of superfamily diversity are recovered from nature and analyzed.

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.

Cysteine oxidation as a regulatory mechanism of Arabidopsis plastidial NAD-dependent malate dehydrogenase

Malate dehydrogenases (MDHs) catalyze a reversible NAD(P)-dependent-oxidoreductase reaction that plays an important role in central metabolism and redox homeostasis of plant cells. Recent studies suggest a moonlighting function of plastidial NAD-dependent MDH (plNAD-MDH; EC 1.1.1.37) in plastid biogenesis, independent of its enzyme activity. In this study, redox effects on activity and conformation of recombinant plNAD-MDH from Arabidopsis thaliana were investigated. We show that reduced plNAD-MDH is active while it is inhibited upon oxidation. Interestingly, the presence of its cofactors NAD+ and NADH could prevent oxidative inhibition of plNAD-MDH. In addition, a conformational change upon oxidation could be observed via non-reducing SDS-PAGE. Both effects, its inhibition and conformational change, were reversible by re-reduction. Further investigation of single cysteine substitutions and mass spectrometry revealed that oxidation of plNAD-MDH leads to oxidation of all four cysteine residues. However, cysteine oxidation of C129 leads to inhibition of plNAD-MDH activity and oxidation of C147 induces its conformational change. In contrast, oxidation of C190 and C333 does not affect plNAD-MDH activity or structure. Our results demonstrate that plNAD-MDH activity can be reversibly inhibited, but not inactivated, by cysteine oxidation and might be co-regulated by the availability of its cofactors in vivo.

Lactate Dehydrogenase Superfamily in Rice and Arabidopsis: Understanding the Molecular Evolution and Structural Diversity

Lactate/malate dehydrogenases (Ldh/Maldh) are ubiquitous enzymes involved in the central metabolic pathway of plants and animals. The role of malate dehydrogenases in the plant system is very well documented. However, the role of its homolog L-lactate dehydrogenases still remains elusive. Though its occurrence is experimentally proven in a few plant species, not much is known about its role in rice. Therefore, a comprehensive genome-wide in silico investigation was carried out to identify all Ldh genes in model plants, rice and Arabidopsis, which revealed Ldh to be a multigene family encoding multiple proteins. Publicly available data suggest its role in a wide range of abiotic stresses such as anoxia, salinity, heat, submergence, cold and heavy metal stress, as also confirmed by our qRT-PCR analysis, especially in salinity and heavy metal mediated stresses. A detailed protein modelling and docking analysis using Schrodinger Suite reveals the presence of three putatively functional L-lactate dehydrogenases in rice, namely OsLdh3, OsLdh7 and OsLdh9. The analysis also highlights the important role of Ser-219, Gly-220 and His-251 in the active site geometry of OsLdh3, OsLdh7 and OsLdh9, respectively. In fact, these three genes have also been found to be highly upregulated under salinity, hypoxia and heavy metal mediated stresses in rice.

Synthetic anaplerotic modules for the direct synthesis of complex molecules from CO2

Anaplerosis is an essential feature of metabolism that allows the continuous operation of natural metabolic networks, such as the citric acid cycle, by constantly replenishing drained intermediates. However, this concept has not been applied to synthetic in vitro metabolic networks, thus far. Here we used anaplerotic strategies to directly access the core sequence of the CETCH cycle, a new-to-nature in vitro CO₂-fixation pathway that features several C₃-C₅ biosynthetic precursors. We drafted four different anaplerotic modules that use CO₂ to replenish the CETCH cycle's intermediates and validated our designs by producing 6-deoxyerythronolide B (6-DEB), the C₂₁-macrolide backbone of erythromycin. Our best design allowed the carbon-positive synthesis of 6-DEB via 54 enzymatic reactions in vitro at yields comparable to those with isolated 6-DEB polyketide synthase (DEBS). Our work showcases how new-to-nature anaplerotic modules can be designed and tailored to enhance and expand the synthetic capabilities of complex catalytic in vitro reaction networks.

Structural Comparison of hMDH2 Complexed with Natural Substrates and Cofactors: The Importance of Phosphate Binding for Active Conformation and Catalysis

Malate dehydrogenase (MDH), which catalyzes a reversible conversion of _L-malate to oxaloacetate, plays essential roles in common metabolic processes, such as the tricarboxylic acid cycle, the oxaloacetate–malate shuttle, and the glyoxylate cycle. MDH2 has lately been recognized as a promising anticancer target; however, the structural information for the human homologue with natural ligands is very limited. In this study, various complex structures of hMDH2, with its substrates and/or cofactors, were solved by X-ray crystallography, which could offer knowledge about the molecular and enzymatic mechanism of this enzyme and be utilized to design novel inhibitors. The structural comparison suggests that phosphate binds to the substrate binding site and brings the conformational change of the active loop to a closed state, which can secure the substate and cofactor to facilitate enzymatic activity.

Identification and characterization of cytosolic malate dehydrogenase from the liver fluke Fasciola gigantica

The liver fluke zoonoses, Fasciola spp. are parasitic helminths infecting humans and animals globally. Recent sequencing of the genome of Fasciola gigantica has provided a basis to understand the biochemistry of this parasite. Here, we identified the cytosolic malate dehydrogenase in F. gigantica (FgMDH) and characterized the enzyme biochemically and structurally. F. gigantica encodes a single cytosolic MDH, a key enzyme of the citric acid cycle. It catalyzes the reversible oxidation of malate to oxaloacetate using NAD⁺. The Fgmdh gene was amplified and cloned for expression of the recombinant protein. The purified protein showed a molecular weight of ~ 36 kDa that existed in a dimeric form in solution. The recombinant enzyme was catalytically active as it catalyzed both forward and reverse reactions efficiently. The kinetic parameters were determined for both directions. The structure of FgMDH and human MDH were modeled and validated. The superimposition of both the model structures showed overall structural similarity in the active site loop region, however, the conformation of the residues was different. Molecular docking elucidated the binding sites and affinities of the substrates and cofactors to the enzyme. Simulation of molecular dynamics and principal component analysis indicated the stability of the systems and collective motions, respectively. Understanding the structural and functional properties of MDH is important to better understand the roles of this enzyme in the biochemistry of the parasite.

Modeling genome-wide enzyme evolution predicts strong epistasis underlying catalytic turnover rates

Systems biology describes cellular phenotypes as properties that emerge from the complex interactions of individual system components. Little is known about how these interactions have affected the evolution of metabolic enzymes. Here, we combine genome-scale metabolic modeling with population genetics models to simulate the evolution of enzyme turnover numbers (k_cats) from a theoretical ancestor with inefficient enzymes. This systems view of biochemical evolution reveals strong epistatic interactions between metabolic genes that shape evolutionary trajectories and influence the magnitude of evolved k_cats. Diminishing returns epistasis prevents enzymes from developing higher k_cats in all reactions and keeps the organism far from the potential fitness optimum. Multifunctional enzymes cause synergistic epistasis that slows down adaptation. The resulting fitness landscape allows k_cat evolution to be convergent. Predicted k_cat parameters show a significant correlation with experimental data, validating our modeling approach. Our analysis reveals how evolutionary forces shape modern k_cats and the whole of metabolism.

Functional and Structural Resilience of the Active Site Loop in the Evolution of Plasmodium Lactate Dehydrogenase

The malarial pathogen Plasmodium falciparum ( Pf) is a member of the Apicomplexa, which independently evolved a highly specific lactate dehydrogenase (LDH) from an ancestral malate dehydrogenase (MDH) via a five-residue insertion in a key active site loop. PfLDH is widely considered an attractive drug target because of its unique active site. The conservation of the apicomplexan loop suggests that a precise insertion sequence was required for the evolution of LDH specificity. Aside from a single critical tryptophan, W107f, the functional and structural roles of residues in the loop are currently unknown. Here we show that the loop is remarkably robust to mutation, as activity is resilient to radical perturbations of both loop identity and length. Thus, alternative insertions could have evolved LDH specificity as long as they contained a tryptophan in the proper location. PfLDH likely has great potential to develop resistance to drugs designed to target its distinctive active site loop.

Structurally Linked Dynamics in Lactate Dehydrogenases of Evolutionarily Distinct Species

We present new findings about how primary and secondary structure affects the role of fast protein motions in the reaction coordinates of enzymatic reactions. Using transition path sampling and committor distribution analysis, we examined the difference in the role of these fast protein motions in the reaction coordinate of lactate dehydrogenases (LDHs) of Apicomplexa organisms Plasmodium falciparum and Cryptosporidium parvum. Having evolved separately from a common malate dehydrogenase ancestor, the two enzymes exhibit several important structural differences, notably a five-amino acid insertion in the active site loop of P. falciparum LDH. We find that these active site differences between the two organisms' LDHs likely cause a decrease in the contribution of the previously determined LDH rate-promoting vibration to the reaction coordinate of P. falciparum LDH compared to that of C. parvum LDH, specifically in the coupling of the rate-promoting vibration and the hydride transfer. This effect, while subtle, directly shows how changes in structure near the active site of LDH alter catalytically important motions. Insights provided by studying these alterations would prove to be useful in identifying LDH inhibitors that specifically target the isozymes of these parasitic organisms.

Comparative Genome-Wide Analysis of the Malate Dehydrogenase Gene Families in Cotton

Malate dehydrogenases (MDHs) play crucial roles in the physiological processes of plant growth and development. In this study, 13 and 25 MDH genes were identified from Gossypium raimondii and Gossypium hirsutum, respectively. Using these and 13 previously reported Gossypium arboretum MDH genes, a comparative molecular analysis between identified MDH genes from G. raimondii, G. hirsutum, and G. arboretum was performed. Based on multiple sequence alignments, cotton MDHs were divided into five subgroups: mitochondrial MDH, peroxisomal MDH, plastidial MDH, chloroplastic MDH and cytoplasmic MDH. Almost all of the MDHs within the same subgroup shared similar gene structure, amino acid sequence, and conserved motifs in their functional domains. An analysis of chromosomal localization suggested that segmental duplication played a major role in the expansion of cotton MDH gene families. Additionally, a selective pressure analysis indicated that purifying selection acted as a vital force in the evolution of MDH gene families in cotton. Meanwhile, an expression analysis showed the distinct expression profiles of GhMDHs in different vegetative tissues and at different fiber developmental stages, suggesting the functional diversification of these genes in cotton growth and fiber development. Finally, a promoter analysis indicated redundant but typical cis-regulatory elements for the potential functions and stress activity of many MDH genes. This study provides fundamental information for a better understanding of cotton MDH gene families and aids in functional analyses of the MDH genes in cotton fiber development.

Function, kinetic properties, crystallization, and regulation of microbial malate dehydrogenase*

Malate dehydrogenase (MDH) is an enzyme widely distributed among living organisms and is a key protein in the central oxidative pathway. It catalyzes the interconversion between malate and oxaloacetate using NAD+ or NADP+ as a cofactor. Surprisingly, this enzyme has been extensively studied in eukaryotes but there are few reports about this enzyme in prokaryotes. It is necessary to review the relevant information to gain a better understanding of the function of this enzyme. Our review of the data generated from studies in bacteria shows much diversity in their molecular properties, including weight, oligomeric states, cofactor and substrate binding affinities, as well as differences in the direction of the enzymatic reaction. Furthermore, due to the importance of its function, the transcription and activity of this enzyme are rigorously regulated. Crystal structures of MDH from different bacterial sources led to the identification of the regions involved in substrate and cofactor binding and the residues important for the dimer-dimer interface. This structural information allows one to make direct modifications to improve the enzyme catalysis by increasing its activity, cofactor binding capacity, substrate specificity, and thermostability. A comparative analysis of the phylogenetic reconstruction of MDH reveals interesting facts about its evolutionary history, dividing this superfamily of proteins into two principle clades and establishing relationships between MDHs from different cellular compartments from archaea, bacteria, and eukaryotes.

Gradual neofunctionalization in the convergent evolution of trichomonad lactate and malate dehydrogenases

Lactate and malate dehydrogenases (LDH and MDH) are homologous, core metabolic enzymes common to nearly all living organisms. LDHs have evolved convergently from MDHs at least four times, achieving altered substrate specificity by a different mechanism each time. For instance, the LDH of anaerobic trichomonad parasites recently evolved independently from an ancestral trichomonad MDH by gene duplication. LDH plays a central role in trichomonad metabolism by catalyzing the reduction of pyruvate to lactate, thereby regenerating the NAD+ required for glycolysis. Using ancestral reconstruction methods, we identified the biochemical and evolutionary mechanisms responsible for this convergent event. The last common ancestor of these enzymes was a highly specific MDH, similar to modern trichomonad MDHs. In contrast, the LDH lineage evolved promiscuous activity by relaxing specificity in a gradual process of neofunctionalization involving one highly detrimental substitution at the “specificity residue” (R91L) and many additional mutations of small effect. L91 has different functional consequences in LDHs and in MDHs, indicating a prominent role for epistasis. Crystal structures of modern‐day and ancestral enzymes show that the evolution of substrate specificity paralleled structural changes in dimerization and α‐helix orientation. The relatively small “specificity residue” of the trichomonad LDHs can accommodate a range of substrate sizes and may permit solvent to access the active site, both of which promote substrate promiscuity. The trichomonad LDHs present a multi‐faceted counterpoint to the independent evolution of LDHs in other organisms and illustrate the diverse mechanisms by which protein function, structure, and stability coevolve.

PDB Code(s): 4UUL; 4UUM; 4UUN; 4UUO; 4UUP; 5A1T

Reassessment of the transhydrogenase/malate shunt pathway in Clostridium thermocellum ATCC 27405 through kinetic characterization of malic enzyme and malate dehydrogenase

Clostridium thermocellum produces ethanol as one of its major end products from direct fermentation of cellulosic biomass. Therefore, it is viewed as an attractive model for the production of biofuels via consolidated bioprocessing. However, a better understanding of the metabolic pathways, along with their putative regulation, could lead to improved strategies for increasing the production of ethanol. In the absence of an annotated pyruvate kinase in the genome, alternate means of generating pyruvate have been sought. Previous proteomic and transcriptomic work detected high levels of a malate dehydrogenase and malic enzyme, which may be used as part of a malate shunt for the generation of pyruvate from phosphoenolpyruvate. The purification and characterization of the malate dehydrogenase and malic enzyme are described in order to elucidate their putative roles in malate shunt and their potential role in C. thermocellum metabolism. The malate dehydrogenase catalyzed the reduction of oxaloacetate to malate utilizing NADH or NADPH with a kcat of 45.8 s(-1) or 14.9 s(-1), respectively, resulting in a 12-fold increase in catalytic efficiency when using NADH over NADPH. The malic enzyme displayed reversible malate decarboxylation activity with a kcat of 520.8 s(-1). The malic enzyme used NADP(+) as a cofactor along with NH4 (+) and Mn(2+) as activators. Pyrophosphate was found to be a potent inhibitor of malic enzyme activity, with a Ki of 0.036 mM. We propose a putative regulatory mechanism of the malate shunt by pyrophosphate and NH4 (+) based on the characterization of the malate dehydrogenase and malic enzyme.

An atomic-resolution view of neofunctionalization in the evolution of apicomplexan lactate dehydrogenases

Malate and lactate dehydrogenases (MDH and LDH) are homologous, core metabolic enzymes that share a fold and catalytic mechanism yet possess strict specificity for their substrates. In the Apicomplexa, convergent evolution of an unusual LDH from MDH produced a difference in specificity exceeding 12 orders of magnitude. The mechanisms responsible for this extraordinary functional shift are currently unknown. Using ancestral protein resurrection, we find that specificity evolved in apicomplexan LDHs by classic neofunctionalization characterized by long-range epistasis, a promiscuous intermediate, and few gain-of-function mutations of large effect. In canonical MDHs and LDHs, a single residue in the active-site loop governs substrate specificity: Arg102 in MDHs and Gln102 in LDHs. During the evolution of the apicomplexan LDH, however, specificity switched via an insertion that shifted the position and identity of this 'specificity residue' to Trp107f. Residues far from the active site also determine specificity, as shown by the crystal structures of three ancestral proteins bracketing the key duplication event. This work provides an unprecedented atomic-resolution view of evolutionary trajectories creating a nascent enzymatic function.

Proteomic analysis of the seed development in Jatropha curcas: from carbon flux to the lipid accumulation

To characterize the metabolic signatures of lipid accumulation in Jatropha curcas seeds, comparative proteomic technique was employed to profile protein changes during the seed development. Temporal changes in comparative proteome were examined using gels-based proteomic technique at six developmental stages for lipid accumulation. And 104 differentially expressed proteins were identified by MALDI-TOF/TOF tandem mass spectrometry. These protein species were classified into 10 functional categories, and the results demonstrated that protein species related to energy and metabolism were notably accumulated and involved in the carbon flux to lipid accumulation that occurs primarily from early to late stage in seed development. Glycolysis and oxidative pentose phosphate pathways were the major pathways of producing carbon flux, and the glucose-6-phosphate and triose-phosphate are the major carbon source for fatty acid synthesis. Lipid analysis revealed that fatty acid accumulation initiated 25days after flowering at the late stage of seed development of J. curcas. Furthermore, C16:0 was initially synthesized as the precursor for the elongation to C18:1 and C18:2 in the developing seeds of J. curcas. Together, the metabolic signatures on protein changes in seed development provide profound knowledge and perspective insights into understanding lipid network in J. curcas.

BIOLOGICAL SIGNIFICANCE

Due to the abundant oil content in seeds, Jatropha curcas seeds are being considered as the ideal materials for biodiesel. Although several studies had carried out the transcriptomic project to study the genes expression profiles in seed development of J. curcas, these ESTs hadn't been confirmed by qRT-PCR. Yet, the seed development of J. curcas had been described for a pool of developing seeds instead of being characterized systematically. Moreover, cellular metabolic events are also controlled by protein-protein interactions, posttranslational protein modifications, and enzymatic activities which cannot be described by transcriptional profiling approaches alone. In this study, within the overall objective of profiling differential protein abundance in developing J. curcas seeds, we provide a setting of physiological data with dynamic proteomic and qRT-PCR analysis to characterize the metabolic pathways and the relationship between mRNA and protein patterns from early stage to seed filling during the seed development of J. curcas. The construction of J. curcas seed development proteome profiles will significantly increase our understanding of the process of seed development and provide a foundation to examine the dynamic changes of the metabolic network during seed development process and certainly suggest some clues to improve the lipid content of J. curcas seeds.

Analysis of quaternary structure of a [LDH-like] malate dehydrogenase of Plasmodium falciparum with oligomeric mutants

L-Malate dehydrogenase (PfMDH) from Plasmodium falciparum, the causative agent for the most severe form of malaria, has shown remarkable similarities to L: -lactate dehydrogenase (PfLDH). PfMDH is more closely related to [LDH-like] MDHs characterized in archae and other prokaryotes. Initial sequence analysis and identification of critical amino acid residues involved in inter-subunit salt-bridge interactions predict tetrameric structure for PfMDH. The catalytically active recombinant PfMDH was characterized as a tetramer. The enzyme is localized primarily in the parasites cytosol. To gain molecular insights into PfMDH/PfLDH relationships and to understand the quaternary structure of PfMDH, dimers were generated by mutation to the potential salt-bridge interacting sites. The R183A and R214G mutations, which snapped the salt bridges between the dimers and resulted in lower dimeric state, did not affect catalytic properties of the enzyme. The mutant dimers of PfMDH were active equally as the wild-type PfMDH. The studies reveal structure of PfMDH as a dimer of dimers. The tetrameric state of PfMDH was not essential for catalytic functions of the enzyme but may be an evolutionary adaptation for cytosolic localization to support its role in NAD/NADH coupling, an important metabolic function for survival of the malaria parasite.

Identification of functional paralog shift mutations: conversion of Escherichia coli malate dehydrogenase to a lactate dehydrogenase

Five positions in the Escherichia coli malate dehydrogenase (eMDH) sequence, which distinguish MDH from lactate dehydrogenase (LDH) activity, were identified through a combination of Venn diagrams constructed from whole genomic data and from unbiased representative sequences from terminal clades. Incorporation of the five changes in eMDH sufficed to convert the enzyme from one with (k(cat)/K(m)(pyruvate))/(k(cat)/K(m)(oxaloacetate)) = 6.1 x 10(-9) to one with that ratio = 28. The substrate specificity was thus changed by a factor of 4.6 x 10(9). The k(cat)/K(m)(pyruvate) value for the pentamutant (eMDH I12V/R81Q/M85E/G210A/V214I) is 3,500 M(-1).s(-1), which is approximately equal 1/1,000 of the values found for typical wild-type LDHs. The procedure isolates an intersection of "strong forcing sets" that should prove to be of general use in switching paralog function.

Purification, characterization, and overexpression of psychrophilic and thermolabile malate dehydrogenase of a novel antarctic psychrotolerant, Flavobacterium frigidimaris KUC-1

We purified the psychrophilic and thermolabile malate dehydrogenase to homogeneity from a novel psychrotolerant, Flavobacterium frigidimaris KUC-1, isolated from Antarctic seawater. The enzyme was a homotetramer with a molecular weight of about 123 k and that of the subunit was about 32 k. The enzyme required NAD(P)(+) as a coenzyme and catalyzed the oxidation of L-malate and the reduction of oxalacetate specifically. The reaction proceeded through an ordered bi-bi mechanism. The enzyme was highly susceptible to heat treatment, and the half-life time at 40 degrees C was estimated to be 3.0 min. The k(cat)/K(m) (microM(-1).s(-1)) values for L-malate and NAD(+) at 30 degrees C were 289 and 2,790, respectively. The enzyme showed pro-R stereospecificity for hydrogen transfer at the C4 position of the nicotinamide moiety of the coenzyme. The enzyme contained 311 amino acid residues and much lower numbers of proline and arginine residues than other malate dehydrogenases.

Probing the specificity of a trypanosomal aromatic alpha-hydroxy acid dehydrogenase by site-directed mutagenesis

The aromatic l-alpha-hydroxy acid dehydrogenase (AHDAH) from Trypanosoma cruzi has over 50% sequence identity with cytosolic malate dehydrogenases (cMDHs), yet it is unable to reduce oxaloacetate. Molecular modeling of the three-dimensional structure of AHADH using the pig cMDH as template directed the construction of several mutants. AHADH shares with MDHs the essential catalytic residues H195 and R171 (using Eventoff's numbering). The AHADH A102R mutant became able to reduce oxaloacetate, while remaining fully active towards aromatic alpha-oxoacids. The Y237G mutant diminished its affinity for all of the natural substrates, whereas the double mutant A102R/Y237G was more active than Y237G and had similar activity with oxaloacetate and with aromatic substrates. The present results reinforce our proposal that AHADH arose by a moderate number of point mutations from a cMDH no longer present in the parasite.

Alteration of the specificity of malate dehydrogenase by chemical modulation of an active site arginine

Malate dehydrogenase from Escherichia coli is highly specific for the oxidation of malate to oxaloacetate. The technique of site-specific modulation has been used to alter the substrate binding site of this enzyme. Introduction of a cysteine in place of the active site binding residue arginine 153 results in a mutant enzyme with diminished catalytic activity, but with K(m) values for malate and oxaloacetate that are surprisingly unaffected. Reaction of this introduced cysteine with a series of amino acid analog reagents leads to the incorporation of a range of functional groups at the active site of malate dehydrogenase. The introduction of a positively charged group such as an amine or an amidine at this position results in improved affinity for several inhibitors over that observed with the native enzyme. However, the recovery of catalytic activity is less dramatic, with less than one third of the native activity achieved with the optimal reagents. These modified enzymes do have altered substrate specificity, with alpha-ketoglutarate and hydroxypyruvate no longer functioning as alternative substrates.

Structural analyses of a malate dehydrogenase with a variable active site

Malate dehydrogenase specifically oxidizes malate to oxaloacetate. The specificity arises from three arginines in the active site pocket that coordinate the carboxyl groups of the substrate and stabilize the newly forming hydroxyl/keto group during catalysis. Here, the role of Arg-153 in distinguishing substrate specificity is examined by the mutant R153C. The x-ray structure of the NAD binary complex at 2.1 A reveals two sulfate ions bound in the closed form of the active site. The sulfate that occupies the substrate binding site has been translated approximately 2 A toward the opening of the active site cavity. Its new location suggests that the low catalytic turnover observed in the R153C mutant may be due to misalignment of the hydroxyl or ketone group of the substrate with the appropriate catalytic residues. In the NAD.pyruvate ternary complex, the monocarboxylic inhibitor is bound in the open conformation of the active site. The pyruvate is coordinated not by the active site arginines, but through weak hydrogen bonds to the amide backbone. Energy minimized molecular models of unnatural analogues of R153C (Wright, S. K., and Viola, R. E. (2001) J. Biol. Chem. 276, 31151-31155) reveal that the regenerated amino and amido side chains can form favorable hydrogen-bonding interactions with the substrate, although a return to native enzymatic activity is not observed. The low activity of the modified R153C enzymes suggests that precise positioning of the guanidino side chain is essential for optimal orientation of the substrate.

The role of active site arginines of sorghum NADP-malate dehydrogenase in thioredoxin-dependent activation and activity

The activation of sorghum NADP-malate dehydrogenase is initiated by thiol/disulfide interchanges with reduced thioredoxin followed by the release of the C-terminal autoinhibitory extension and a structural modification shaping the active site into a high efficiency and high affinity for oxaloacetate conformation. In the present study, the role of the active site arginines in the activation and catalysis was investigated by site-directed mutagenesis and arginyl-specific chemical derivatization using butanedione. Sequence and mass spectrometry analysis were used to identify the chemically modified groups. Taken together, our data reveal the involvement of Arg-134 and Arg-204 in oxaloacetate coordination, suggest an indirect role for Arg-140 in substrate binding and catalysis, and clearly confirm that Arg-87 is implicated in cofactor binding. In contrast with NAD-malate dehydrogenase, no lactate dehydrogenase activity could be promoted by the R134Q mutation. The decreased susceptibility of the activation of the R204K mutant to NADP and its increased sensitivity to the histidine-specific reagent diethylpyrocarbonate indicated that Arg-204 is involved in the locking of the active site. These results are discussed in relation with the recently published NADP-MDH three-dimensional structures and the previously established three-dimensional structures of NAD-malate dehydrogenase and lactate dehydrogenase.

From malate dehydrogenase to phenyllactate dehydrogenase. Incorporation of unnatural amino acids to generate an improved enzyme-catalyzed activity

Malate dehydrogenase (MDH) from Escherichia coli is highly specific for its keto acid substrate. The placement of the active site-binding groups in MDH effectively discriminates against both the shorter and the longer keto dicarboxylic acids that could potentially serve as alternative substrates. A notable exception to this specificity is the alternative substrate phenylpyruvate. This aromatic keto acid can be reduced by MDH, albeit at a somewhat slower rate and with greatly diminished affinity, despite the presence of several substrate-binding arginyl residues and the absence of a hydrophobic pocket in the active site. The specificity of MDH for phenylpyruvate has now been enhanced, and that for the physiological substrate oxaloacetate has been diminished, through the replacement of one of the binding arginyl residues with several unnatural alkyl and aryl amino acid analogs. This approach, called site-specific modulation, incorporates systematic structural variations at a site of interest. Molecular modeling studies have suggested a structural basis for the affinity of native MDH for phenylpyruvate and a rationale for the improved catalytic activity that is observed with these new, modified phenyllactate dehydrogenases.

A family of highly conserved glycosomal 2-hydroxyacid dehydrogenases from Phytomonas sp

Phytomonas sp. contains two malate dehydrogenase isoforms, a mitochondrial isoenzyme with a high specificity for oxaloacetate and a glycosomal isozyme that acts on a broad range of substrates (Uttaro, A. D., and Opperdoes, F.R. (1997) Mol. Biochem. Parasitol. 89, 51-59). Here, we show that the low specificity of the latter isoenzyme is the result of a number of recent gene duplications that gave rise to a family of glycosomal 2-hydroxyacid dehydrogenase genes. Two of these genes were cloned, sequenced, and overexpressed in Escherichia coli. Although both gene products have 322 amino acids, share 90.4% identical residues, and have a similar hydrophobicity profile and net charge, their kinetic properties were strikingly different. One isoform behaved as a real malate dehydrogenase with a high specificity for oxaloacetate, whereas the other showed no activity with oxaloacetate but was able to reduce other oxoacids, such as phenyl pyruvate, 2-oxoisocaproate, 2-oxovalerate, 2-oxobutyrate, 2-oxo-4-methiolbutyrate, and pyruvate.

The NAD-linked aromatic alpha-hydroxy acid dehydrogenase from Trypanosoma cruzi. A new member of the cytosolic malate dehydrogenases group without malate dehydrogenase activity

Trypanosoma cruzi, the protozoan parasite causing Chagas disease, contains a novel aromatic alpha-hydroxy acid dehydrogenase. This enzyme is responsible, together with tyrosine aminotransferase, for the catabolism of aromatic amino acids, which leads to the excretion of aromatic lactate derivatives into the culture medium. The gene encoding the aromatic alpha-hydroxy acid dehydrogenase has been cloned through a combined approach using screening of an expression genomic library with antibodies, peptide sequencing and PCR amplification. Its sequence shows high similarity to the cytosolic malate dehydrogenases. However, the enzyme has no malate dehydrogenase activity. The gene seems to be present in a single copy per haploid genome and is differentially expressed throughout the parasite's life cycle, the highest levels being found in the insect forms of T. cruzi. The purified recombinant enzyme, expressed in Escherichia coli, was unable to reduce oxaloacetate and had kinetic constants similar to those of the natural aromatic alpha-hydroxy acid dehydrogenase. Sequence comparisons suggest that the aromatic alpha-hydroxy acid dehydrogenase derives from a cytosolic malate dehydrogenase no longer present in the parasite, made redundant by the presence of a glycosomal malate dehydrogenase as a member of a shuttle device involving the mitochondrial isoenzyme.

Convergent evolution of Trichomonas vaginalis lactate dehydrogenase from malate dehydrogenase

Lactate dehydrogenase (LDH) is present in the amitochondriate parasitic protist Trichomonas vaginalis and some but not all other trichomonad species. The derived amino acid sequence of T. vaginalis LDH (TvLDH) was found to be more closely related to the cytosolic malate dehydrogenase (MDH) of the same species than to any other LDH. A key difference between the two T. vaginalis sequences was that Arg91 of MDH, known to be important in coordinating the C-4 carboxyl of oxalacetate/malate, was replaced by Leu91 in LDH. The change Leu91Arg by site-directed mutagenesis converted TvLDH into an MDH. The reverse single amino acid change Arg91Leu in TvMDH, however, gave a product with no measurable LDH activity. Phylogenetic reconstructions indicate that TvLDH arose from an MDH relatively recently.

Malate/lactate dehydrogenase in mollicutes: evidence for a multienzyme protein

The malate (MDH) and lactate (LDH) dehydrogenases belong to the homologous class of 2-ketoacid dehydrogenases. The specificity for their respective substrates depends on residues differing at two or three regions within each molecule. Theoretical peptide-mass fingerprinting and PROSITE analysis of nine MDH and six LDH molecules were used to describe conserved sites related to function. A unique LDH is described which probably also confers MDH activity within the 580 kbp genome of Mycoplasma genitalium (class: Mollicutes). A single hydrophilic arginine residue was found in the active site of the M. genitalium LDH enzyme, differing from an hydrophobic residue normally present in these molecules. The effect of this residue may be to alter active site substrate specificity, allowing the enzyme to perform two closely related tasks. Evidence for a single gene affording dual enzymatic function is discussed in terms of genome size reduction in the simplest of free-living organisms. Since Mollicutes are thought to lack enzymes of the tricarboxylic acid cycle that would otherwise bind and interact with MDH in bacterial species possessing this pathway, active site modification of M. genitalium LDH is the sole requirement for MDH activity of this molecule. The closely related helical Mollicute, Spiroplasma melliferum, was shown to possess two distinct gene products for MDH/LDH activity.

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.

Tetrameric malate dehydrogenase from a thermophilic Bacillus: cloning, sequence and overexpression of the gene encoding the enzyme and isolation and characterization of the recombinant enzyme

The gene encoding the tetrameric malate dehydrogenase (MDH) in a thermophilic Bacillus species (BI) has been cloned in an Escherichia coli plasmid. The nucleotide sequence of the gene, the first to be elucidated for a tetrameric MDH, shows the MDH subunit to contain 312 amino acids and have a molecular mass of 33648 Da, which confirms the experimentally determined value of about 35 kDa. Like the genomic DNA of BI, the MDH gene is relatively AT-rich; this contrasts with the generally GC-rich nature of the DNA of thermophilic Bacillus species. Comparison of amino acid sequences reveals that BI MDH bears greater structural similarity to lactate dehydrogenases (LDHs) than to other (dimeric) MDHs. MDHs and LDHs resemble each other in catalytic mechanism and several other respects. However, whereas MDHs in the majority of organisms are dimers, the tetrameric structure is favoured among LDHs. The stronger structural resemblance that BI MDH has to LDHs than to the dimeric MDHs provides some explanation as to why Bacillus MDH, unlike most other MDHs, is tetrameric. A 1 kb fragment containing the BI MDH gene, produced in a PCR, has been cloned into a high-expression E. coli plasmid vector. BI MDH synthesized from this clone constitutes about 47% of the total protein in cell extracts of the E. coli strain carrying the clone. MDH purified from BI and that purified from the E. coli strain carrying the MDH gene clone appear to be identical proteins by several criteria. A number of characteristics of the MDH have been elucidated, including the molecular masses of the native enzyme and the subunit, N-terminal amino acid sequence, isoelectric point, pH optimum for activity, thermostability, stability to pH, urea and guanidinium chloride and several kinetic parameters. Whereas the MDH is a stable tetramer in the pH range 5-7, it appears to be converted into a stable dimer at pH 3.5. This suggests that the dimer is a stable intermediate in the dissociation of the tetramer to monomers at low pH.

A highly active decarboxylating dehydrogenase with rationally inverted coenzyme specificity

The isocitrate dehydrogenase of Escherichia coli, which lacks the Rossmann fold common to other dehydrogenases, displays a 7000-fold preference for NADP over NAD (calculated as the ratio of kcat/Km). Guided by x-ray crystal structures and molecular modeling, site-directed mutagenesis has been used to introduce six substitutions in the adenosine binding pocket that systematically shift coenzyme preference toward NAD. The engineered enzyme displays an 850-fold preference for NAD over NADP, which exceeds the 140-fold preference displayed by a homologous NAD-dependent enzyme. Of the six mutations introduced, only one is identical in all related NAD-dependent enzyme sequences--strict adherence to homology as a criterion for replacing these amino acids impairs function. Two additional mutations at remote sites improve performance further, resulting in a final mutant enzyme with kinetic characteristics and coenzyme preference comparable to naturally occurring homologous NAD-dependent enzymes.

Malate dehydrogenase: a model for structure, evolution, and catalysis

Malate dehydrogenases are widely distributed and alignment of the amino acid sequences show that the enzyme has diverged into 2 main phylogenetic groups. Multiple amino acid sequence alignments of malate dehydrogenases also show that there is a low degree of primary structural similarity, apart from in several positions crucial for nucleotide binding, catalysis, and the subunit interface. The 3-dimensional structures of several malate dehydrogenases are similar, despite their low amino acid sequence identity. The coenzyme specificity of malate dehydrogenase may be modulated by substitution of a single residue, as can the substrate specificity. The mechanism of catalysis of malate dehydrogenase is similar to that of lactate dehydrogenase, an enzyme with which it shares a similar 3-dimensional structure. Substitution of a single amino acid residue of a lactate dehydrogenase changes the enzyme specificity to that of a malate dehydrogenase, but a similar substitution in a malate dehydrogenase resulted in relaxation of the high degree of specificity for oxaloacetate. Knowledge of the 3-dimensional structures of malate and lactate dehydrogenases allows the redesign of enzymes by rational rather than random mutation and may have important commercial implications.

A single amino acid mutation enhances the thermal stability of Escherichia coli malate dehydrogenase

The stability of wild-type Escherichia coli malate dehydrogenase was compared with a mutant form of the enzyme with the amino acid residue at position 102 changed from arginine to glutamine. The mutation occurs on the underside of a mobile loop which closes over the active-site cleft on formation of the enzyme/cofactor/substrate ternary complex. The mutant enzyme is kinetically compromised while the wild-type enzyme is highly specific for oxaloacetate. The mutant enzyme was shown to be more resistant to irreversible thermal denaturation by thermal inactivation experiments and high-sensitivity differential scanning calorimetry than the wild-type enzyme. In contrast, resistance of both enzymes to reversible unfolding in guanidinium chloride was similar. Circular dichroic spectropolarimetry shows the secondary structures of the enzymes are similar but there is a demonstrable difference in tertiary structure. From the position of the mutation, it is conjectured that the substitution on a mobile surface loop results in partial closure of the loop and greater resistance to thermal inactivation of the mutant enzyme. However, molecular modelling combined with circular dichroic spectropolarimetry indicate that the mutation may have a more widespread effect on the structure than simply partial closure of the mobile surface loop as the environment of distant tyrosine residues is altered. Resistance of the wild-type enzyme to thermal inactivation can be increased by cofactor addition, which may have the effect of partial closure of the mobile surface loop, but has little effect on the mutant enzyme.

Substitution of the amino acid at position 102 with polar and aromatic residues influences substrate specificity of lactate dehydrogenase

The Gln residue at amino acid position 102 of Bacillus stearothermophilus lactate dehydrogenase was replaced with Ser, Thr, Tyr, or Phe to investigate the effect on substrate recognition. The Q102S and Q102T mutant enzymes were found to have a broader range of substrate specificity (measured by kcat/Km) than the wild-type enzyme. However, it is evident that either Ser or Thr at position 102 are of a size able to accommodate a wide variety of substrates in the active site and substrate specificity appears to rely largely on size discrimination in these mutants. The Q102F and Q102Y mutant enzymes have low catalytic efficiency and do not show this relaxed substrate specificity. However, their activities are restored by the presence of an aromatic substrate. All of the enzymes have a very low catalytic efficiency with branched chain aliphatic substrates.