Structures of the apo and holo forms of formate dehydrogenase from the bacteriumMoraxellasp. C-1: towards understanding the mechanism of the closure of the interdomain cleft

Effect of Additional Amino Acid Replacements on the Properties of Multi-point Mutant Bacterial Formate Dehyderogenase PseFDH SM4S

Formate dehydrogenase from Pseudomonas sp. 101 bacterium (PseFDH, EC 1.2.1.2) is a research model for the elucidation of the catalytic mechanism of 2-oxyacid D-specific dehydrogenases enzyme superfamily. The enzyme is actively used for regeneration of the reduced form of NAD(P)H in chiral synthesis with oxidoreductases. A multi-point mutant PseFDH SM4S with an improved thermal and chemical stability has been prepared earlier in this laboratory. To further improve the properties of the mutant, additional single-point replacements have been introduced to generate five new PseFDH mutants. All new enzymes have been highly purified, and their kinetic properties and thermal stability studied using analysis of thermal inactivation kinetics and differential scanning calorimetry. The E170D amino acid change in PseFDH SM4S shows an increase in thermal stability 1.76- and 10-fold compared to the starting mutant and the wild-type enzyme, respectively.

Role of Protein Motions in Catalysis by Formate Dehydrogenase

We have analyzed the reaction catalyzed by formate dehydrogenase using transition path sampling. This system has recently received experimental attention using infrared spectroscopy and heavy-enzyme studies. Some of the experimental results point to the possible importance of protein motions that are coupled to the chemical step. We found that the residue Val123 that lies behind the nicotinamide ring occasionally comes into van der Waals contact with the acceptor and that in all reactive trajectories, the barrier-crossing event is preceded by this contact, meaning that the motion of Val123 is part of the reaction coordinate. Experimental results have been interpreted with a two-dimensional formula for the chemical rate, which cannot capture effects such as the one we describe.

Effect of Metal Ions on the Activity of Ten NAD-Dependent Formate Dehydrogenases

NAD-dependent formate dehydrogenase (FDH) enzymes are frequently used in industrial and scientific applications. FDH is a reversible enzyme that reduces the NAD molecule to NADH and produces CO₂ by oxidation of the formate ion, whereas it causes CO₂ reduction in the reverse reaction. Some transition metal elements - Fe³⁺, Mo⁶⁺ and W^{6 +} - can be found in the FDH structure of anaerobic and archaeal microorganisms, and these enzymes require cations and other redox-active cofactors for their FDH activity. While NAD-dependent FDHs do not necessarily require any metal cations, the presence of various metal cations can still affect FDH activities. To study the effect of 11 different metal ions, NAD-dependent FDH enzymes from ten different microorganisms were tested: Ancylobacter aquaticus (AaFDH), Candida boidinii (CboFDH), Candida methylica (CmFDH), Ceriporiopsis subvermispora (CsFDH), Chaetomium thermophilum (CtFDH), Moraxella sp. (MsFDH), Myceliophthora thermophila (MtFDH), Paracoccus sp. (PsFDH), Saccharomyces cerevisiae (ScFDH) and Thiobacillus sp. (TsFDH). It was found that metal ions (mainly Cu²⁺ and Zn²⁺) could have quite strong inhibition effects on several enzymes in the forward reaction, whereas several cations (Li⁺, Mg²⁺, Mn²⁺, Fe³⁺ and W⁶⁺) could increase the forward reaction of two FDHs. The highest activity increase (1.97 fold) was caused by Fe³⁺ in AaFDH. The effect on the reverse reaction was minimal. The modelled structures of ten FDHs showed that the active site is formed by 15 highly conserved amino acid residues spatially settling around the formate binding site in a conserved way. However, the residue differences at some of the sites close to the substrate do not explain the activity differences. The active site space is very tight, excluding water molecules, as observed in earlier studies. Structural examination indicated that smaller metal ions might be spaced close to the active site to affect the reaction. Metal ion size showed partial correlation to the effect on inhibition or activation. Affinity of the substrate may also affect the sensitivity to the metal's effect. In addition, amino acid differences on the protein surface may also be important for the metal ion effect.

Classification, substrate specificity and structural features of D-2-hydroxyacid dehydrogenases: 2HADH knowledgebase

BACKGROUND

The family of D-isomer specific 2-hydroxyacid dehydrogenases (2HADHs) contains a wide range of oxidoreductases with various metabolic roles as well as biotechnological applications. Despite a vast amount of biochemical and structural data for various representatives of the family, the long and complex evolution and broad sequence diversity hinder functional annotations for uncharacterized members.

RESULTS

We report an in-depth phylogenetic analysis, followed by mapping of available biochemical and structural data on the reconstructed phylogenetic tree. The analysis suggests that some subfamilies comprising enzymes with similar yet broad substrate specificity profiles diverged early in the evolution of 2HADHs. Based on the phylogenetic tree, we present a revised classification of the family that comprises 22 subfamilies, including 13 new subfamilies not studied biochemically. We summarize characteristics of the nine biochemically studied subfamilies by aggregating all available sequence, biochemical, and structural data, providing comprehensive descriptions of the active site, cofactor-binding residues, and potential roles of specific structural regions in substrate recognition. In addition, we concisely present our analysis as an online 2HADH enzymes knowledgebase.

CONCLUSIONS

The knowledgebase enables navigation over the 2HADHs classification, search through collected data, and functional predictions of uncharacterized 2HADHs. Future characterization of the new subfamilies may result in discoveries of enzymes with novel metabolic roles and with properties beneficial for biotechnological applications.

Temperature dependence of dynamic, tunnelling and kinetic isotope effects in formate dehydrogenase

The origin of the catalytic power of enzymes has been a question of debate for a long time. In this regard, the possible contribution of protein dynamics in enzymatic catalysis has become one of the most controversial topics. In the present work, the hydride transfer step in the formate dehydrogenase (FDH EC 1.2.1.2) enzyme is studied by means of molecular dynamic (MD) simulations with quantum mechanics/molecular mechanics (QM/MM) potentials in order to explore any correlation between dynamics, tunnelling effects and the rate constant. The temperature dependence of the kinetic isotope effects (KIEs), which is one of the few tests that can be studied by experiments and simulations to shed light on this debate, has been computed and the results have been compared with previous experimental data. The classical mechanical free energy barrier and the number of recrossing trajectories seem to be temperature-independent while the quantum vibrational corrections and the tunnelling effects are slightly temperature-dependent over the interval of 5-45 °C. The computed primary KIEs are in very good agreement with previous experimental data, being almost temperature-independent within the standard deviations. The modest dependence on the temperature is due to just the quantum vibrational correction contribution. These results, together with the analysis of the evolution of the collective variables such as the electrostatic potential or the electric field created by the protein on the key atoms involved in the reaction, confirm that while the protein is well preorganised, some changes take place along the reaction that favour the hydride transfer and the product release. Coordinates defining these movements are, in fact, part of the real reaction coordinate.

Structural, Biochemical, and Evolutionary Characterizations of Glyoxylate/Hydroxypyruvate Reductases Show Their Division into Two Distinct Subfamilies

The d-2-hydroxyacid dehydrogenase (2HADH) family illustrates a complex evolutionary history with multiple lateral gene transfers and gene duplications and losses. As a result, the exact functional annotation of individual members can be extrapolated to a very limited extent. Here, we revise the previous simplified view on the classification of the 2HADH family; specifically, we show that the previously delineated glyoxylate/hydroxypyruvate reductase (GHPR) subfamily consists of two evolutionary separated GHRA and GHRB subfamilies. We compare two representatives of these subfamilies from Sinorhizobium meliloti (SmGhrA and SmGhrB), employing a combination of biochemical, structural, and bioinformatics approaches. Our kinetic results show that both enzymes reduce several 2-ketocarboxylic acids with overlapping, but not equivalent, substrate preferences. SmGhrA and SmGhrB show highest activity with glyoxylate and hydroxypyruvate, respectively; in addition, only SmGhrB reduces 2-keto-d-gluconate, and only SmGhrA reduces pyruvate (with low efficiency). We present nine crystal structures of both enzymes in apo forms and in complexes with cofactors and substrates/substrate analogues. In particular, we determined a crystal structure of SmGhrB with 2-keto-d-gluconate, which is the biggest substrate cocrystallized with a 2HADH member. The structures reveal significant differences between SmGhrA and SmGhrB, both in the overall structure and within the substrate-binding pocket, offering insight into the molecular basis for the observed substrate preferences and subfamily differences. In addition, we provide an overview of all GHRA and GHRB structures complexed with a ligand in the active site.

Structural basis for double cofactor specificity in a new formate dehydrogenase from the acidobacterium Granulicella mallensis MP5ACTX8

Formate dehydrogenases (FDHs) are considered particularly useful enzymes in biocatalysis when the regeneration of the cofactor NAD(P)H is required, that is, in chiral synthesis with dehydrogenases. Their utilization is however limited to the recycling of NAD(+), since all (apart one) of the FDHs characterized so far are strictly specific for this cofactor, and this is a major drawback for their otherwise wide applicability. Despite the many attempts performed to modify cofactor specificity by protein engineering different NAD(+)-dependent FDHs, in the general practice, glucose or phosphite dehydrogenases are chosen for the recycling of NADP(+). We report on the functional and structural characterization of a new FDH, GraFDH, identified by mining the genome of the extremophile prokaryote Granulicella mallensis MP5ACTX8. The new enzyme displays a valuable stability in the presence of many organic cosolvents as well as double cofactor specificity, with NADP(+) preferred over NAD(+) at acidic pH values, at which it also shows the highest stability. The quite low affinities for both cofactors as well as for the substrate formate indicate, however, that the native enzyme requires optimization to be applied as biocatalytic tool. We also determined the crystal structure of GraFDH both as apoprotein and as holoprotein, either in complex with NAD(+) or NADP(+). Noticeably, the latter represents the first structure of an FDH enzyme in complex with NADP(+). This fine picture of the structural determinants involved in cofactor selectivity will possibly boost protein engineering of the new enzyme or other homolog FDHs in view of their biocatalytic exploitation for NADP(+) recycling.

Improvement of the soy formate dehydrogenase properties by rational design

Previous experiments on substitution of the residue Phe290 to Asp, Asn and Ser in NAD(+)-dependent formate dehydrogenase from soya Glycine max (SoyFDH) showed important role of the residue in enzyme thermal stability and catalytic properties (Alekseeva et al. Prot. Eng. Des. Sel., 2012a; 25: :781-88). In this work, we continued site-directed mutagenesis experiments of the Phe290 and the residue was changed to Ala, Thr, Tyr, Glu and Gln. All amino acid changes resulted in increase of catalytic constant from 2.9 to 3.5-4.7 s(-1). The substitution Phe290Ala led to KM (NAD+) decrease from 13.3 to 8.6 μM, and substitutions Phe290Tyr and Phe290Glu resulted in decrease and increase of KM (HCOO-) from 1.5 to 0.9 and -2.9 mM, respectively. The highest improvement of catalytic properties was observed for SoyFDH Phe290Ala which showed 2-fold higher catalytic efficiency with both substrates. Stability of mutants was examined by study of thermal inactivation kinetics and differential scanning calorimetry (DSC). All five amino acids provided increase of thermal stability of mutant SoyFDH in comparison with wild-type enzyme. Mutant SoyFDH Phe290Glu showed the highest improvement-the stabilization effect was 43 at 56°C. The DSC data agree with results of thermal inactivation kinetics. Substitutions Phe290Tyr, Phe290Thr, Phe290Gln and Phe290Glu provided Tm value increase 0.6°-6.6°. SoyFDH Phe290Glu and previously prepared SoyFDH Phe290Asp show similar thermal stability as enzymes from Candida boidinii and Mycobacterium vaccae N10 and have higher catalytic efficiency with NAD(+) compared with all described FDHs. Therefore, these mutants are very perspective enzymes for coenzyme regeneration in processes of chiral synthesis with dehydrogenases.

Structural insights into the efficient CO2-reducing activity of an NAD-dependent formate dehydrogenase from Thiobacillus sp. KNK65MA

CO2 fixation is thought to be one of the key factors in mitigating global warming. Of the various methods for removing CO2, the NAD-dependent formate dehydrogenase from Candida boidinii (CbFDH) has been widely used in various biological CO2-reduction systems; however, practical applications of CbFDH have often been impeded owing to its low CO2-reducing activity. It has recently been demonstrated that the NAD-dependent formate dehydrogenase from Thiobacillus sp. KNK65MA (TsFDH) has a higher CO2-reducing activity compared with CbFDH. The crystal structure of TsFDH revealed that the biological unit in the asymmetric unit has two conformations, i.e. open (NAD(+)-unbound) and closed (NAD(+)-bound) forms. Three major differences are observed in the crystal structures of TsFDH and CbFDH. Firstly, hole 2 in TsFDH is blocked by helix α20, whereas it is not blocked in CbFDH. Secondly, the sizes of holes 1 and 2 are larger in TsFDH than in CbFDH. Thirdly, Lys287 in TsFDH, which is crucial for the capture of formate and its subsequent delivery to the active site, is an alanine in CbFDH. A computational simulation suggested that the higher CO2-reducing activity of TsFDH is owing to its lower free-energy barrier to CO2 reduction than in CbFDH.

Engineering catalytic properties and thermal stability of plant formate dehydrogenase by single-point mutations

The analysis of the 3D model structure of the ternary complex of recombinant formate dehydrogenase from soya Glycine max (EC 1.2.1.2., SoyFDH) with bound NAD+ and an inhibitor azide ion revealed the presence of hydrophobic Phe290 in the coenzyme-binding domain. This residue should shield the enzyme active site from solvent. On the basis of the alignment of plant FDHs sequences, Asp, Asn and Ser were selected as candidates to substitute Phe290. Computer modeling indicated the formation of two (Ser and Asn) or three (Asp) new hydrogen bonds in such mutants. The mutant SoyFDHs were expressed in Escherichia coli, purified and characterized. All amino acid substitutions increased K(м)(HCOO-) from 1.5 to 4.1-5.0 mM, whereas the K(м)(NAD+) values remained almost unchanged in the range from 9.1 to 14.0 μM, which is close to wt-SoyFDH (13.3 μM). The catalytic constants for F290N, F290D and F290S mutants of SoyFDH equaled 2.8, 5.1 and 4.1 s⁻¹, respectively; while that of the wild-type enzyme was 2.9 s⁻¹. The thermal stability of all mutant SoyFDHs was much higher compared with the wild-type enzyme. The differential scanning calorimetry data were in agreement with the results of thermal inactivation kinetics. The mutations F290S, F290N and F290D introduced into SoyFDH increased the T(m) values by 2.9°C, 4.3°C and 7.8°C, respectively. The best mutant F290D exhibited thermal stability similar to that of FDH from the plant Arabidopsis thaliana and exceeded that of the enzymes from the yeast Candida boidinii and the bacterium Moraxella sp. C1.

Investigation of formate transport through the substrate channel of formate dehydrogenase by steered molecular dynamics simulations

Steered molecular dynamics simulation has revealed the mechanism of formate transport via the substrate channel of formate dehydrogenase. It is shown that the structural organization of the channel promotes the transport of formate anion in spite of the fact that the channel is too narrow even for such a small molecule. The conformational mobility of Arg284 residue, one of the residues forming the wall of the substrate channel, provides for the binding and delivery of formate to the active site.

Formate formation and formate conversion in biological fuels production

Biomethanation is a mature technology for fuel production. Fourth generation biofuels research will focus on sequestering CO(2) and providing carbon-neutral or carbon-negative strategies to cope with dwindling fossil fuel supplies and environmental impact. Formate is an important intermediate in the methanogenic breakdown of complex organic material and serves as an important precursor for biological fuels production in the form of methane, hydrogen, and potentially methanol. Formate is produced by either CoA-dependent cleavage of pyruvate or enzymatic reduction of CO(2) in an NADH- or ferredoxin-dependent manner. Formate is consumed through oxidation to CO(2) and H(2) or can be further reduced via the Wood-Ljungdahl pathway for carbon fixation or industrially for the production of methanol. Here, we review the enzymes involved in the interconversion of formate and discuss potential applications for biofuels production.

Crystallization and preliminary X-ray analysis of formate oxidase, an enzyme of the glucose-methanol-choline oxidoreductase family

Formate oxidase (FOD), which catalyzes the oxidation of formate to yield carbon dioxide and hydrogen peroxide, belongs to the glucose-methanol-choline oxidoreductase (GMCO) family. FOD from Aspergillus oryzae RIB40, which has a modified FAD as a cofactor, was crystallized at 293 K by the hanging-drop vapour-diffusion method. The crystal was orthorhombic and belonged to space group C222(1). Diffraction data were collected from a single crystal to 2.4 A resolution.