N,N-dimethylglycine, a dietary supplement, has been reported to be beneficial in children with autism and pervasive developmental disorder. We examined the effectiveness of dimethylglycine in children with autism and pervasive developmental disorder in a double-blind, placebo-controlled study. Thirty-seven children between 3 and 11 years of age with a diagnosis of autism and/or pervasive developmental disorder were gender and age matched and randomly assigned to receive either placebo or dimethylglycine for 4 weeks. All children were assessed before and after treatment on two behavioral measures, the Vineland Maladaptive Behavior Domain and the Aberrant Behavior Checklist. Standardized neurologic examinations before and after treatment on 33 children showed no change. An overall improvement on all behavioral measures was observed for both the placebo and the dimethylglycine groups. However, the improvement among the children who received dimethylglycine was not statistically different from the improvement observed among the children who received the placebo. The children who participated in this study were a heterogeneous group, and their apparent responses to the dimethylglycine varied. Some children appeared to respond positively to the dimethylglycine, and there was a smaller proportion of negative changes in the dimethylglycine group, but the quantitative changes in the dimethylglycine behavioral assessments were not significantly different from what was observed among children who received placebo. (J Child Neurol 2001;16:169-173).
The human dimethylglycine dehydrogenase (hDMGDH) is a flavin adenine dinucleotide (FAD)‐ and tetrahydrofolate (THF)‐dependent, mitochondrial matrix enzyme taking part in choline degradation, one‐carbon metabolism and electron transfer to the respiratory chain. The rare natural variant H109R causes dimethylglycine dehydrogenase deficiency leading to increased blood and urinary dimethylglycine concentrations. A detailed biochemical and structural characterization of hDMGDH was thus far hampered by insufficient heterologous expression of the protein. In the present study, we report the development of an intracellular, heterologous expression system in Komagataella phaffii (formerly known as Pichia pastoris) providing the opportunity to determine kinetic parameters, spectroscopic properties, thermostability, and the redox potential of hDMGDH. Moreover, we have successfully crystallized the wild‐type enzyme and determined the structure to 3.1‐Å resolution. The structure‐based analysis of our biochemical data provided new insights into the kinetic properties of the enzyme in particular with respect to oxygen reactivity. A comparative study with the H109R variant demonstrated that the variant suffers from decreased protein stability, cofactor saturation, and substrate affinity.
Structural data are available in the PDB database under the accession number 5L46.
Keywords: electron transfer, flavin adenine dinucleotide, genetic disease, recombinant protein expression, X‐ray crystallography
The human flavoproteome comprises 90 enzymes with versatile functions, structures, and protein characteristics. Sixty percent of the known flavoproteins are involved in human diseases and disorders emphasizing their importance in human metabolism 1. The human mitochondrial matrix flavoprotein dimethylglycine dehydrogenase (hDMGDH, EC: 22.214.171.124) catalyzes the oxidative demethylation of dimethylglycine (DMG) to sarcosine, and also to a lesser extent the conversion of sarcosine to glycine as part of choline degradation (Fig. 1) 2, 3. Choline is an essential nutrient and building block in a variety of vital biomolecules such as the membrane phospholipid phosphatidylcholine and the neurotransmitter acetylcholine 4. Degradation of choline proceeds by consecutive oxidations via betaine aldehyde, betaine, dimethylglycine, and sarcosine to the amino acid glycine (Fig. 1) 5. hDMGDH, a key enzyme of this pathway, requires two cofactors: flavin adenine dinucleotide (FAD) and tetrahydrofolate (THF). The FAD is covalently attached via its 8α‐position to the N3 of a histidyl residue and serves as the electron acceptor in the oxidation of DMG. On the other hand, THF is used as the acceptor of the incipient methyl group and thus prevents the release of cell‐toxic formaldehyde during catalysis 3, 6. In the course of this reaction, N‐5,10‐methylene tetrahydrofolate is formed, which plays an important role in one‐carbon metabolism. Regeneration of oxidized FAD is achieved by electron transfer to the human electron‐transferring flavoprotein (hETF), which in turn transfers the electrons to the membrane‐anchored ETF‐ubiquinone oxidoreductase (ETF‐QO) for further utilization in the mitochondrial respiratory chain 7, 8.
Mammalian choline degradation pathway with DMGDH as metabolic branch point. Choline degradation to dimethylglycine is catalyzed by three enzymes, choline dehydrogenase (CHD), betaine aldehyde dehydrogenase (BAD), and betaine homocysteine S‐methyltransferase...
A rare naturally occurring point mutation found in 58 out of 118 656 (0.049%) analyzed human gene sequences (Exome Aggregation Consortium online browser; http://exac.broadinstitute.org/variant/5-78351682-T-C) results in the exchange of histidine 109 to arginine (H109R). This variant reportedly suffers from lower stability and enzymatic activity 9, 10. The observed phenotype was described as dimethylglycine dehydrogenase deficiency with significantly higher levels of dimethylglycine in human body fluids causing muscle fatigue and a fish‐like odor (OMIM: 605850). Interestingly, the mutation is predominantly found in individuals of African descent.
In recent years, also the product of the reaction catalyzed by hDMGDH, sarcosine, attracted attention as a biomarker for aggressive prostate cancer 11, 12. Therefore, investigations on the regulation of sarcosine levels and the enzymes involved are in the focus of research. Furthermore, a recent epidemiological study revealed a possible connection of dimethylglycine dehydrogenase deficiency to the development of diabetes, further emphasizing the importance of the enzyme 13.
Dimethylglycine dehydrogenase (DMGDH) was first identified and purified from rat liver in the 1950s and early 1960s 14, 15. Since then several publications have dealt with the characterization of mammalian dimethylglycine dehydrogenases from rat 2, 3, 16 and pig liver 8 as well as the recombinant rat 17, 18 and human enzyme 19 using Escherichia coli as expression host. Furthermore, dimethylglycine oxidase from Arthrobacter globiformis (AgDMGO) was characterized in biochemical and structural detail 6, 20, 21, 22, 23. However, the available information on the hDMGDH appears insufficient and partly contradictory. In order to overcome these deficiencies, we present a profound analysis of the recombinant hDMGDH and the H109R variant concerning kinetics, redox behavior, and spectral properties. Furthermore, we have elucidated the crystal structure of the human wild‐type enzyme.
In order to characterize hDMGDH, we have established a recombinant expression system for the wild‐type and the H109R variant in Komagataella phaffii (formerly known as Pichia pastoris or Komagataella pastoris24) and a subsequent purification protocol. This approach enabled us to determine key kinetic parameters as well as physical and spectroscopic properties of wild‐type DMGDH. Consequently, we have also employed our expression system to generate the H109R variant to address the putative loss of function that leads to DMGDH‐deficiency in humans. Moreover, successful crystallization of the WT and subsequent crystallographic analysis provided us with the opportunity to analyze the structure with regard to substrate binding and oxygen reactivity. The latter issue is especially interesting as the control of oxygen reactivity in flavoenzymes is still a controversial topic. The direct comparison of a dehydrogenase (hDMGDH) and an oxidase (AgDMGO) allowed new insights into the structural elements involved in oxygen reactivity and supports concepts that are based on gatekeeper residues in the vicinity of the isoalloxazine ring system. In this vein, we show that a previously proposed model to rationalize oxygen reactivity in the vanillyl oxidase family is applicable to the family of sarcosine and dimethylglycine dehydrogenases and oxidases, respectively 25, 26.
Enzyme expression and purification
Initial attempts to express the gene encoding hDMGDH in Escherichia coli (BL21 DE3) yielded largely insoluble protein. On the other hand, heterologous expression in the methanotrophic yeast Komagataella phaffi (formerly known as Pichia pastoris24) was successful. As shown in Fig. 2A, western blot analysis of K. phaffii cell lysates at different time points indicates a stable expression of the protein after induction with methanol (MeOH). Typically, fermentations were stopped after 96 h after MeOH induction resulting in 1.6–1.9 kg of wet cell pellet. A comparison of signal intensities showed that the WT was expressed in higher amounts than the variant under identical fermentation conditions thus leading to a higher yield of WT (Fig. 2A). After cell disruption and Ni‐NTA affinity chromatography, the yield of proteins (Fig. 2B, lanes 3) was approximately 70 and 25 mg of the WT and H109R variant, respectively (40 or 15 μg enzyme per g wet cell weight). In order to achieve higher purity for crystallization trials, the proteins were further purified using anion exchange chromatography (Fig. 2B, lanes 4) resulting in lower protein yields of 30 and 2 mg (17 or 1 μg enzyme per g wet cell weight) of WT and H109R variant, respectively. The protein loss mainly occurred during the necessary buffer change to lower salt concentrations after Ni‐NTA affinity chromatography required for the subsequent anion exchange chromatography. Although we experienced that hDMGDH tends to precipitate at low salt concentrations, other chromatographic methods, like size‐exclusion chromatography or hydrophobic interaction chromatography, were explored but did not give a similar purity.
Heterologous expression and purification of hDMGDH‐WT and hDMGDH‐H109R variant. (A) Western blot of K. phaffii cell lysates taken at different time points after methanol induction. Antibodies are directed against a C‐terminal nona‐histidine...
The A280/A450 ratios of the highly pure protein fractions were usually between 14–16 and 20–25 for the WT and variant protein, respectively. The UV‐Vis absorption spectra show a clear difference between enzyme‐bound FAD and free FAD (Fig. 2C), whereas hDMGDH‐WT and the variant featured similar spectral properties with identical absorption maxima. In comparison with free FAD in solution, the UV‐Vis absorption maxima of FAD bound to either WT or the H109R variant exhibit a hypsochromic shift from 370 to 350 nm and a bathochromic shift from 445 to 460 nm (Fig. 2C).
Determination of kinetic parameters for hDMGDH‐WT and H109R variant
The evaluation of the steady‐state kinetic parameters for hDMGDH‐WT and H109R variant was conducted using either ferrocene or 2,6‐dichlorophenolindophenol (DCPIP) as electron acceptor to oxidize the reduced FAD cofactor. The obtained parameters show significant differences between the two assay systems both in terms of the kcat as well as the KM for DMG (Fig. 3, Table 1). When ferrocene was used as the artificial electron acceptor the kcat was fivefold greater than with DCPIP (213 ± 4 versus 44 ± 1 min−1, Table 1) and the KM was significantly higher (1.4 ± 0.1 versus 0.3 ± 0.002 mm, Table 1). Interestingly, the H109R variant exhibited a twofold higher activity than the WT in both assays. At the same time, the KM of DMG was strongly increased by a factor of 23 and 16 in the ferrocene and DCPIP assay, respectively, resulting in an approximately 10‐fold lower catalytic efficiency of the H109R variant (Table 1). Notably, when DCPIP was used as electron acceptor, data linearization using the method of Eadie and Hofstee revealed a second, much higher KM value of ≈ 30 mm for DMG (inset Fig. 3A). As the rate at which FAD is oxidized by ferrocene is closer to the rate of flavin reduction, it appears that ferrocene is superior to DCPIP as electron acceptor.
Steady‐state kinetics of hDMGDH‐WT and H109R variant. Steady‐state kinetic experiments were performed with DCPIP (A) and ferrocene activity assays (B) by plotting initial reaction velocities (v0) against the substrate concentration....
Summary of kinetic parameters obtained with steady‐state kinetics for hDMGDH‐WT and hDMGDH‐H109R variant with two different activity assays (DCPIP and ferrocene)
Due to much lower enzyme amounts available of the H109R variant, presteady‐state kinetics were solely measured with the wild‐type enzyme. As shown in Fig. 4A, the rate of reduction as a function of DMG concentration fitted to a hyperbolic equation yielding a limiting reductive rate of 17 ± 0.3 s−1 and a dissociation constant of 4.9 ± 0.3 mm for DMG. Thus, the reductive rate is 5–20 times faster than kcat determined in the ferrocene and DCPIP (Table 1) indicating that the electron transfer from the reduced FAD to the artificial electron acceptor is rate‐limiting in the assays used to evaluate steady‐state kinetics. Reoxidation of hDMGDH with air‐saturated buffer proceeds very slowly with an observed rate of reoxidation of 0.006 ± 0.001 s−1 at an oxygen concentration of 135 μm (inset Fig. 4A). Rapid reaction measurements with sarcosine yielded a limiting rate of reduction of 0.7 ± 0.1 s−1 and a dissociation constant of 280 ± 10 mm clearly indicating that sarcosine is a much poorer substrate than DMG (Fig. 4A).
Presteady‐state kinetics of hDMGDH‐WT. Reductive and oxidative rates were measured under anoxic conditions with the stopped flow device. (A) Reductive rates were determined for DMG (black, left axis) and sarcosine (red, right axis). The...
Finally, we also studied the effect of THF on the activity of hDMGDH using the ferrocene assay (Fig. 5). The activity measurements were done similarly to the ferrocene activity assay as described above at 25 °C and pH 7.8 in 50 mm HEPES/NaOH with 150 mm NaCl, 100 nm hDMGDH, 1 mm DMG, 0.1 mm EDTA, and 0.2 mm freshly prepared ferrocene in a concentration range of 0–75 μm THF. As shown in Fig. 5, the activity of WT decreases as a function of the THF concentration.
hDMGDH‐WT activity as a function of the THF concentration. The activity was measured using the ferrocene activity assay. The error bars show the standard deviation of at least three measurements for each THF concentration.
Photoreduction and determination of the redox potential of hDMGDH
Reduction of hDMGDH by DMG shows a monophasic conversion of the oxidized FAD to the two‐electron reduced dihydroquinone state (Fig. 4B) without any transient appearance of a semiquinone radical. In contrast to substrate reduction, photoreduction of hDMGDH leads to the anionic (red) flavin semiquinone with a typical absorption maximum at 372 nm (Fig. 6A). Further photoreduction eventually leads to the fully reduced dihydroquinone (Fig. 6B), which is directly converted to the oxidized state without the occurrence of a semiquinone upon reaction with molecular dioxygen (Fig. 6C). This was also seen when reduced hDMGDH was reacted with air‐saturated buffer in the stopped flow device (data not shown).
Photoreduction and redoxpotential of hDMGDH‐WT. (A) and (B) show selected absorption spectra of the anoxic photoreduction of the wild‐type enzyme. About 20 μm of hDMGDH was photoreduced in the presence of 1 mm EDTA, 2 μ...
In order to determine the redox potential of hDMGDH, we employed the xanthine oxidase/xanthine method described by Massey 27 using indigotrisulfonic acid potassium salt as reference dye. The slope of the double logarithmic plot was close to unity indicating that the dye and the flavin receive an equal number of electrons (i.e. two), although the anionic semiquinone radical is clearly observable during the redox titration (Fig. 6D). The redox potential for the reduction of FAD bound to hDMGDH was calculated from six independent measurements to −93 ± 1 mV.
Thermal stability of hDMGDH‐WT and hDMGDH‐H109R
The thermal stability of both enzymes was determined using Thermofluor® and circular dichroism spectroscopy. Thermofluor® temperature scans from 25 to 95 °C showed that the WT and the H109R variant have melting points (Tm) of 52 °C and 47 °C, respectively (Fig. 7A). A similar result was obtained by means of CD spectroscopy yielding a Tm of 53 °C and 47 °C, for WT and the variant, respectively (Fig. 7B). Thus, both independent methods confirm that the variant exhibits reduced thermal stability.
Thermal stability of hDMGDH‐WT and the H109R variant. The melting temperatures of wild‐type and the variant enzyme differ significantly when measured with Thermofluor® (A) or CD spectroscopy (B). Thermofluor® was done with...
Protein structure analysis
The X‐ray crystal structure of hDMGDH was determined at 3.1 Å resolution (Table 2). The crystal belongs to the monoclinic space group (P21) and contains two dehydrogenase molecules in the asymmetric unit. The two molecules in the asymmetric unit are very similar to each other with an rmsd of 0.22 Å calculated after the superposition of 810 out of 828 Cα‐atoms using the program pymol (DeLano Scientific). As previously observed in rDMGDH and also after superposition with the rat enzyme (rmsd of 0.37 Å calculated for 739 out of 809 Cα‐atoms), the human enzyme consists of two domains: the FAD‐binding domain (N‐terminal domain with residues 46–466) and the folate‐binding domain (C‐terminal domain with residues 467–855, Fig. 8A). The FAD cofactor is covalently attached to the protein through a linkage between N3 of His91 and the 8α‐methyl group of the isoalloxazine moiety as confirmed by the crystal structure. The folate‐binding domain can be divided into three subdomains positioned in a ‘cloverleaf’‐like arrangement as previously described for the rat enzyme 18. Subdomain 1 (residues 467–537 and 626–728) contains a Greek‐key motif surrounded by α‐helices. Subdomain 2 (residues 538–625 and 729–770) is defined by a five‐stranded antiparallel β‐sheet with flanking α‐helices. Subdomain 3 (residues 771–850) adopts a jellyroll fold. The structure of the rDMGDH was determined in complex with THF 18. We also observed some difference electron density in the same area in the folate‐binding domain. Due to the lower resolution of the diffraction data and/or an incomplete occupancy of the ligand, however, this density was not clear enough to unequivocally position a folate molecule. Based on the structural similarity of the two enzymes, it is still very likely that folate binds in the same position (Fig. 8B).
Crystal structure of hDMGDH. (A) The overall structure of hDMGDH with the FAD‐binding domain is shown in deep blue and the folate‐binding domain in raspberry. The THF‐binding site is depicted with an arrow and the FAD cofactor...
Data collection and refinement statistics
Heterologous enzyme expression and purification
In this study, we demonstrated that intracellular expression of hDMGDH in K. phaffii is feasible and provides enough hDMGDH to perform a detailed biochemical and structural characterization. The expression of the full‐length protein and a Δ50 truncated version as well as secretory expression were unsuccessful. However, an N‐terminal Δ28 truncation of the enzyme was expressed successfully representing the mature form of the enzyme lacking its mitochondrial targeting sequence 28. The intracellular coexpression of the Saccharomyces cerevisiae protein disulfide isomerase enzyme resulted in improved protein yields. The expression of hDMGDH in E. coli as reported before 19 was unsuccessful and failed to yield detectable amounts of protein as judged by western blotting. A C‐terminal nona‐histidine‐tag was added to the gene in order to facilitate purification by means of affinity chromatography. According to the previously published structure of the rat enzyme 18, this C‐terminal tag should not interfere with the native fold of the protein or catalytic activity. Western blot analysis showed that the H109R variant is expressed in lower amounts compared to the WT (Fig. 2A).
A theoretical A280/A450 ratio for hDMGDH fully saturated with FAD cofactor can be roughly calculated to 14. The observed A280/A450 ratios of the highly pure proteins, 14–16 and 20–24 for WT and H109R variant, respectively, revealed that the WT is almost fully loaded with the covalently attached FAD. In the case of the variant, the ratio was much higher indicating a much lower saturation compared to the WT (ca. 40–50% apoprotein), similar to previously reported observations 19. The higher fraction of apoprotein may be responsible for the recurring protein precipitation during purification of the variant. Overall, this led to substantially reduced yields of the variant protein and limited the scope of experiments that could be conducted.
Spectral characteristics and steady‐state kinetics of hDMGDH
Recombinant hDMGDH features characteristic UV‐Vis absorption properties that clearly distinguish the bound FAD from FAD free in solution (Fig. 2C). These spectral characteristics are very similar to those reported earlier for DMGDH isolated from rat and pig liver as well as for recombinant AgDMGO 8, 16, 20. In contrast, the UV‐Vis absorption spectra for heterologously expressed rat DMGDH 17 and human DMGDH 19 are either of very poor quality or apparently lack the spectral features that are typically found in these enzymes. In fact, the only published spectrum of hDMGDH expressed in E. coli is very similar to free FAD in solution (compare spectra shown in Fig. 9).
Protein spectrum of this study (black line) compared with the published hDMGDH spectrum (gray line, 19) and free FAD in solution (dotted blue line). The published spectrum resembles free FAD in solution and not our spectral findings.
A summary of our data and the previously published kinetic steady‐state data is given in Table 3. These values were obtained using two different enzymatic activity assays, on one hand with ferrocenium hexafluorophosphate as electron acceptor and on the other hand, the combination of phenazine methosulfate (PMS) as electron mediator and DCPIP as terminal electron acceptor. Both assays are physiologically irrelevant as all mammalian DMGDHs interact with the electron‐transferring flavoprotein (hETF) for reoxidation of the reduced cofactor. With the ferrocene assay, we have obtained comparable results as reported for the rat enzyme 18 and also for the bacterial DMGO 22. However, analysis of the data obtained with the PMS/DCPIP yielded two Michaelis–Menten constants similar to previously published data for the rat enzyme 16, 17. Interestingly, aside from the results reported by McAndrew et al. 19, two KM values were only observed for the two component DCPIP activity assay, with PMS as electron mediator, but not in the one component ferrocene or H2O2 assay suggesting that it is an experimental artifact. In view of this apparent inconsistency, speculations that the presence of two KM values reflect the presence of activating or regulatory sites in the enzyme are unwarranted 16, 17. Moreover, we have shown that the rate of enzyme reduction is faster than kcat indicating that reoxidation is rate‐limiting for turnover in both assays. Thus, the affinity of DMG to hDMGDH is more reliably expressed as the dissociation constant determined in presteady‐state experiments.
Comparison of steady‐state kinetic parameters of published results for the DMGDH from Homo sapiens (human), Rattus norvegicus (rat), Sus scrofa (pig), as well as the dimethylglycine oxidase from A. globiformis (Ag)
Effects of a naturally occurring H109R enzyme variation
The rare naturally occurring enzyme variant H109R causes a nonfatal disease called dimethylglycine dehydrogenase deficiency (OMIM: 605850). The defect results in an accumulation of DMG in blood serum and urine, and patients suffer from a mild phenotype including muscle fatigue and fish‐like odor 9, 10. During heterologous expression in K. phaffii, we noticed decreased expression levels resulting in much lower yields of the variant protein. Further characterization of the variant protein demonstrated lower cofactor saturation and thermal stability. Interestingly, the variant showed elevated turnover rates in both of our assay systems although KM values were 15–25 times higher and the catalytic efficiency 10‐fold lower than for the WT (Table 1). Hence, our results are clearly in contradiction with previously published data 19 that claimed significantly decreased catalytic activity of the H109R variant. As the spectral properties of DMGDH preparations obtained from the E. coli expression system lacks the characteristic features typically observed for DMGDH isolated from natural sources (e.g., rat and pig liver), it is conceivable that recombinant DMGDH generated in bacterial host cells suffers from an altered FAD‐binding mode, such as noncovalent binding. In fact, the absorption spectrum reported for wild‐type DMGDH and its H109R variant is clearly more similar to that of free FAD (Fig. 9). In our hands, the bacterial expression system produced solely insoluble protein and therefore analysis of FAD binding to the protein could not be inspected further.
Previous reports on the structure of bacterial DMGO and rat DMGDH identified a large internal cavity in the protein 6, 18. This cavity serves as ‘reaction chamber’ responsible for substrate delivery to the active site, substrate conversion, as well as intermediate channeling to THF and product formation. The binding site of THF is the only entry‐exit point of the protein, which was revealed by MD simulations for the structurally related DMGO 6 and also confirmed by cavity analysis of our structure (Fig. 10A). In accordance with this interpretation, increasing THF concentrations exhibited a negative effect on the turnover rate (Fig. 5) as the substrate DMG competes with THF for entrance into the cavity at the THF‐binding site.
hDMGDH internal cavity analysis and effect of H109R variation. (A) Cavity analysis of hDMGDH done with CASoX in pymol without a bound THF cofactor. The protein cavity is shown in green surface view and the THF‐binding site indicated by an arrow....
The consistently (16‐ to 23‐fold) higher KM observed in both steady‐state assays with the variant suggests that replacement of histidine by the bulkier arginine impedes access of DMG to the binding site (Fig. 10B). On the other hand, the H109R variant possessed higher kcat values indicating that the variant is more permissible in transferring electrons from the reduced FAD to the artificial electron acceptors.
In conclusion, our analysis of the disease‐related H109R variant has revealed a lower thermal stability and an approximately 10‐fold lower catalytic efficiency (Table 1) as major effects of the single amino acid exchange. The diminished protein expression and cofactor saturation found in our expression system may also be relevant for homologous expression of the variant in human cells and thus contribute to the overall reduced DMGDH activity manifested in affected human subjects.
Substrate specificity, redox behavior, and active site composition
The overall protein structure and the active site of hDMGDH with the covalently bound FAD as well as the THF‐binding site are very similar to the rat enzyme (PDB ID: 4PAA; 18), and the structure of AgDMGO 21. In 2003, AgDMGO could be crystallized with a bound folic acid (PDB ID: 1PJ6) but additionally also with a bound acetate ion that mimics the negatively charged carboxylate group of DMG (PDB ID: 1PJ5). Most residues in the active site of AgDMGO described to play a role in catalysis and substrate binding are also present in hDMGDH (Fig. 11). Docking of DMG and sarcosine to the active site of the hDMGDH structure results in a similar positioning of both substrates (Fig. 11A,B) and thus do not rationalize the differences observed in affinity and catalytic activity. This binding mode is similar to the acetate moiety in the crystal structure of DMGO (Fig. 11C). Presteady‐state measurements comparing DMG and sarcosine as substrate for hDMGDH‐WT show a 25‐fold decrease of the reductive rate for sarcosine as the substrate and a more than 50‐fold higher KD. These results are in very good agreement with data obtained for the bacterial DMGO (Table 4) 22. Therefore, sarcosine is a poor but possibly physiologically relevant substrate for the enzyme as already shown previously for the pig enzyme 3 and AgDMGO 22. On the other hand, sarcosine is apparently oxidized by a homologous sarcosine dehydrogenase (SARDH) 16 and also by peroxisomal sarcosine oxidase (PIPOX) 29. Therefore, the relative turnover of sarcosine by these three FAD‐dependent enzymes must await further studies.
Docking and active site composition of hDMGDH in comparison with bound acetate in DMGO structure. DMG (A) and sarcosine (B) were docked into the active site of hDMGDH with autodock/vina47. The active site of DMGO with bound acetate (C) was retrieved...
Summary and comparison of transient kinetic data of hDMGDH (this study) with published results from A. globiformis DMGO 22, 23. Measurements were done in a stopped flow at pH 7.8 and 25 °C for hDMGDH and 25 °C and pH 8.5 for the bacterial...
In the bacterial enzyme, a catalytic dyad, consisting of a histidine and a tyrosine, is responsible for deprotonation of the substrate amine 21. This catalytic dyad is also present in the human enzyme (Fig. 11). Furthermore, residues interacting with DMG (and sarcosine) were identified in the structure of hDMGDH (Fig. 11A). After binding and conversion of the substrate at the active site, an intermediate, most likely a cyclic lactone is channeled through the protein cavity to the THF‐binding site 6. The THF‐binding sites of mammalian and bacterial enzymes are again very similar and were already extensively discussed by Luka et al. 18.
Basran et al. reported a multistep kinetic model for the reduction of DMGO by DMG. In contrast, the reduction of hDMGDH by DMG was adequately fit by a single exponential equation. The actual measurement of the intermediate cyclic lactone could not be observed, most likely due to its short lifetime as described elsewhere 6.
The two dehydrogenases in choline degradation belong to a family of enzymes that delivers electrons to the electron transferring flavoprotein (ETF). Among these clients, only the redox potential for the two‐electron reduction of human medium‐chain acyl‐CoA dehydrogenases (hMCAD) was previously determined 30, 31. This prompted us to determine the redox potential for the covalently bound FAD in hDMGDH. The obtained value of −93 ± 1 mV is 40 mV more positive than that of hMCAD (−135 mV 30, 31) in agreement with the fact that covalent linkages to the 8‐methyl group increase the redox potential 32. In any case, the much higher redox potential of hETF (+37 mV for the oxidized/semiquinone couple, 33) ensures that electrons from hDMGDH are readily transferred to the ETF‐bound FAD.
The rates of reduction determined for hDMGDH are in good agreement with those observed for bacterial DMGO (Table 4). In contrast to that, the oxygen reactivity of reduced hDMGDH is more than 300 times lower and apparently physiologically irrelevant as reoxidation by its cognate partner protein hETF proceeds over 30 times faster (P. Augustin, E. Gerstmann and P. Macheroux, unpublished results).
Recently, a study concerning steric control of dioxygen reduction in enzymes of the vanillyl‐alcohol oxidase (VAO) family discovered a new ‘gatekeeper’ residue in berberine bridge enzyme and the pollen allergen Phl p 4 25. This gatekeeper appears to control access to an oxyanion hole essential to stabilize reaction intermediates during oxidation of the reduced isoalloxazine ring by dioxygen. In this family, valine was found to grant access to the oxyanion hole, whereas isoleucine denies access 25. Similarly, Leferink et al. discovered that a large number of oxidases in the VAO family contain either a glycine or a proline at a structurally conserved position near the isoalloxazine ring 26. Although hDMGDH does not belong to the VAO protein superfamily, comparison of the hDMGDH to the AgDMGO structure reveals that hDMGDH features residues typical for dehydrogenases suppressing oxygen reactivity, i.e., an alanine and isoleucine, whereas AgDMGO possesses residues compatible with high oxygen reactivity, i.e., proline and valine (Fig. 12). Nevertheless, in order to further support, if the emerging concept of how oxygen reactivity is controlled in flavoproteins is applicable to DMGDH, a member of the d‐amino acid oxidase (DAO) family of FAD‐dependent oxidoreductases (pfam: 01266), a detailed mutagenesis study of the concerned residues has to be conducted.
Active site composition of the oxygen‐binding site of hDMGDH and AgDMGO. The human structure is shown in panel A with the most significant alternated residues on the si‐face of the flavin isoalloxazine ring shown in green sticks. The...
Enzymes and reagents
Restriction enzymes and Phusion DNA polymerase were from Thermo Fisher Scientific (Waltham, MA, USA), purification columns and materials from GE Healthcare (Chalfont St Giles, UK) and salt‐free purified oligos for site‐directed mutagenesis from VBC‐Biotech (Vienna, Austria). Antibodies used for western blot were from Cell Signaling Technology (Cambridge, UK). The redox dye for the determination of the redox potential, indigotrisulfonic acid potassium salt, was from TCI Europe N.V. (Zwijndrecht, Belgium). All other chemicals and media were from Carl Roth GmbH (Karlsruhe, Germany) and Sigma‐Aldrich (St. Louis, MO, USA) and were of the highest grade available.
hDMGDH expression strain generation
hDMGDH and hDMGDH‐H109R expression plasmids were designed for intracellular enzyme expression in Komagataella phaffii KM71H strain (formerly known as Pichia pastoris or Komagataella pastoris24). The hDMGDH sequence was codon optimized for expression in K. phaffii using GeneOptimizer® (Thermo Fisher Scientific). In accordance to Binzak et al., the sequence of the mature form of hDMGDH (lacking the first 28 amino acids, 28) was used, flanked with XhoI and NotI restriction sites. For enzyme expression and purification, a Kozak sequence in front of the start codon (AAAA) and a C‐terminal nona‐histidine tag were added. The designed gene was cloned into a pPICZ B expression vector (Thermo Fisher Scientific) and verified by automated sequencing. The recombinant plasmid was transformed to K. phaffii KM71H strain (Thermo Fisher Scientific), harboring a pPICK‐PDI vector containing the gene for the protein disulfide isomerase from Saccharomyces cerevisiae following the guidelines of the Invitrogen EasySelect™Pichia Expression Kit (Thermo Fisher Scientific). Positive clones were selected by Zeocin resistance and verified by western blot. The H109R variant was constructed by a two‐step site‐directed mutagenesis method employing Phusion DNA polymerase and 5′‐AACTTGAAGAAGATCAGATACGACTCCATCAAG‐3′ as forward and 5′‐CTTGATGGAGTCGTATCTGATCTTCTTCAAGTT‐3′ as reverse primer (alternated codon underlined). For the first reaction step, two separate PCR reactions using either the forward or the reverse primer were run with the PCR program: 98 °C (2 min) – [98 °C (50 s) – 60 °C (20 s) – 68 °C (16.5 min)] × 5 – 4 °C ∞. Afterwards, the two PCR reactions were combined and the same program was employed for another 20 cycles. Then, after a 2‐h DpnI digestion step, the plasmid was transformed to Escherichia coli TOP10 cells for amplification before transformation to K. phaffii KM71H [pPICK‐PDI]. Strain selection was done as described above.
hDMGDH and hDMGDH‐H109R expression was carried out in a 7‐L glass fermenter BBI CT5‐2 system (Sartorius, Göttingen, Germany). Preparation of inoculum, and preparation of the fermenter, batch, and fed‐batch were carried out as described by Schrittwieser et al. 34. The fermentation medium was a minimal basal salt medium (MgSO4·7 H2O – 1.2 g·L−1; K2SO4 – 14.9 g·L−1; KOH – 18.2 g·L−1; NaCl – 4.13 g·L−1; glycerol – 40 g·L−1; 85% H3PO4 – 27 mL·L−1). After the fed‐batch, the MeOH induction was started by injection of 5 g MeOH directly into the fermenter and afterwards by a stepwise increase in MeOH feed from 0 to 6 g·h−1 over 2 h. Glycerol feed was reduced from 15 to 6 g·h−1 over 2 h and a mixed MeOH/glycerol feed was maintained until the end of the fermentation. In total, 400–500 g MeOH was fed to the cultures before harvest after 96 h of induced growth. Samples were taken for western blot analysis at different time points. The cell pellet was collected by centrifugation (2704 g, 30 min, 4 °C) and stored at −20 °C.
Cell lysates were prepared by Merckenschlager (Braun Biotech International, Melsungen, Germany) glass beads cell homogenization under CO2 cooling (3‐min disruption steps). For every 1 g harvested wet cell pellet, 2‐mL cell lysis buffer (50 mm HEPES/NaOH, 150 mm NaCl, 35 mm imidazole, 1 mm PMSF, pH 7.8) and a spatula tip of FAD were added before cell disruption. The lysates were cleared by centrifugation (38 720 g, 30 min, 4 °C) and filtered through a paper filter. Nickel ion affinity chromatography was done by applying the lysates onto self‐packed nickel‐Sepharose 6 Fast Flow columns (GE Healthcare) at 4 °C. Afterwards, the material was washed with washing buffer (50 mm HEPES/NaOH, 150 mm NaCl, 75 mm imidazole, pH 7.8) and the enzyme stripped with elution buffer (50 mm HEPES/NaOH, 150 mm NaCl, 300 mm imidazole, pH 7.8). The purification was monitored by SDS/PAGE and fractions containing hDMGDH were concentrated and rebuffered to storage buffer (50 mm HEPES/NaOH, 150 mm NaCl, pH 7.8) using Amicon ultracentrifugal filter units (50 kDa cut‐off; Merck‐Milipore, Darmstadt, Germany). The obtained enzyme purity was sufficient for kinetic and photometric studies of the enzymes. If a higher enzyme purity was needed, especially for protein crystallization, the eluted fractions of the affinity chromatography were rebuffered to buffer A (25 mm HEPES/NaOH, 25 mm NaCl, pH 7.8), and anion exchange chromatography with a MonoQ 5/50 GL column (GE Healthcare) was done on an ÄKTA FPLC system (GE Healthcare) at 4 °C. The protein was eluted with a gradient from 0% to 50% buffer B (25 mm HEPES/NaOH, 1 m NaCl, pH 7.8). The chromatography was controlled by SDS/PAGE and fractions containing hDMGDH were concentrated and rebuffered to storage buffer as described above. Purification was done identically for the wild‐type and the variant enzyme.
Enzyme samples were separated by SDS/PAGE, with 12.5% separation gels (for western blot 10%), under reducing conditions as described by Laemmli 35. Gels were either used for western blot analysis or stained with Coomassie brilliant blue R‐250 for purification control. As protein standard, a PageRuler Prestained protein ladder (Thermo Fisher Scientific) was employed.
K. phaffii cell lysates for western blot screening and control of fermentation were prepared by glass bead disruption following the Invitrogen EasySelect™Pichia Expression Kit manual (p. 42, Thermo Fisher Scientific). Western blot analysis was done on nitrocellulose membranes essentially following the General Protocol for western blotting from Bio‐Rad (Hercules, CA, USA; Bulletin 6376). As antibodies, a rabbit anti‐histidine IgG antibody at 1 : 2000 overnight at 4 °C as primary and an HRP‐linked goat anti‐rabbit IgG antibody at 1 : 5000 for 1 h at room temperature as secondary antibody were used. As loading control, a rabbit anti‐GAPDH antibody was used in a 1 : 1000 dilution overnight. Immunoreactive signals were obtained using SuperSignal West Pico Chemiluminescent Substrate (Thermo Fisher Scientific) and detection in a G:BOX (Syngene, Cambridge, UK).
UV‐Vis absorption spectroscopy
UV‐Vis absorption spectra to assess protein concentration, activity, purity, quality, and cofactor saturation, as well as for steady‐state kinetic measurements and photoreduction were recorded with a Specord 210 spectrophotometer (Analytik Jena, Jena, Germany).
Enzyme quantification and calculation of the extinction coefficient
Protein concentrations of purified hDMGDH enzymes were determined according to the characteristic absorption of bound FAD at 450 nm. The molar extinction coefficient for hDMGDH was calculated using the method described by Macheroux 36 to 11 600 m−1·cm−1.
Steady‐state kinetic parameters were determined using DCPIP according to Okamura‐Ikeda et al. 37 or ferrocenium hexafluorophosphate (Fc+) as described by Lehman et al. 38 as terminal electron acceptor and measurable dimension in the spectrophotometer. In short, DCPIP assays were performed at 25 °C and pH 7.8 in 50 mm HEPES/NaOH, 150 mm NaCl, 135 μm DCPIP, and 100 nm hDMGDH with freshly prepared 3 mm phenazine methosulfate (PMS) as intermediate electron mediator and 0–100 mm dimethylglycine (DMG) following the change of absorption spectrophotometrically at 600 nm over 3 min. Ferrocenium assays were performed at 25 °C at pH 7.8 in 50 mm HEPES/NaOH, 150 mm NaCl, 0.1 mm EDTA, 100 nm hDMGDH, 0–100 mm DMG, and freshly prepared 0.2 mm Fc+ following the change of absorption at 300 nm over 3 min. For each concentration, at least a triplicate measurement was performed, the initial velocities were determined, and KM and kcat were assessed by nonlinear hyperbolic fit in origin 8.6 (OriginLab Corp., Northampton, MA, USA).
Presteady‐state reaction kinetics were measured anaerobically with a Hi‐Tech stopped flow instrument (SF‐61DX2; TgK Scientific Limited, Bradford‐on‐Avon, UK) in a glove box (Belle Technology, Weymouth, UK) at 25 °C. For determination of the reductive rates of enzyme‐bound FAD, 20 μm purified hDMGDH in 50 mm HEPES/NaOH, pH 7.8 containing 150 mm NaCl was shot against different concentrations of DMG (0.25–250 mm) or sarcosine (0–500 mm) dissolved in 50 mm HEPES/NaOH, 150 mm NaCl pH 7.8 and monitored at 460 nm with a KinetaScanT diode array detector (MG‐6560; TgK Scientific Limited). Each concentration was measured in triplicates and the observed rate constants (kobs) for different substrate concentrations were calculated using an exponential fitting function in the kinetic studio software (TgK Scientific Limited). The dissociation constant (KD) was determined by a hyperbolic fitting curve employing origin 8.6 (OriginLab Corp.). The oxidative rate (kox) of the enzyme was measured three times at an oxygen concentration of 10.5% O2 (135 μm O2), by mixing substrate‐reduced hDMGDH with air equilibrated buffer (50 mm HEPES/NaOH, 150 mm NaCl, pH 7.8).
Determination of the redox potential
The redox potential (E0