Structural perturbations in the Ala --> Val polymorphism of methylenetetrahydrofolate reductase: how binding of folates may protect against inactivation

In human methylenetetrahydrofolate reductase (MTHFR) the Ala222Val (677C-->T) polymorphism encodes a heat-labile gene product that is associated with elevated levels of homocysteine and possibly with risk for cardiovascular disease. Generation of the equivalent Ala to Val mutation in Escherichia coli MTHFR, which is 30% identical to the catalytic domain of the human enzyme, creates a protein with enhanced thermolability. In both human and E. coli MTHFR, the A --> V mutation increases the rate of dissociation of FAD, and in both enzymes, loss of FAD is linked to changes in quaternary structure [Yamada, K., Chen, Z., Rozen, R., and Matthews, R. G. (2001) Proc. Natl. Acad. Sci. U.S.A. 98, 14853-14858; Guenther, B. D., Sheppard, C. A., Tran, P., Rozen, R., Matthews, R. G., and Ludwig, M. L. (1999) Nat. Struct. Biol. 6, 359-365]. Folates have been shown to protect both human and bacterial enzymes from loss of FAD. Despite its effect on affinity for FAD, the A --> V mutation is located at the bottom of the (betaalpha)(8) barrel of the catalytic domain in a position that does not contact the bound FAD prosthetic group. Here we report the structures of the Ala177Val mutant of E. coli MTHFR and of its complex with the 5,10-dideazafolate analogue, LY309887, and suggest mechanisms by which the mutation may perturb FAD binding. Helix alpha5, which immediately precedes the loop bearing the mutation, carries several residues that interact with FAD, including Asn168, Arg171, and Lys172. In the structures of the mutant enzyme this helix is displaced, perturbing protein-FAD interactions. In the complex with LY309887, the pterin-like ring of the analogue stacks against the si face of the flavin and is secured by hydrogen bonds to residues Gln183 and Asp120 that adjoin this face. The direct interactions of bound folate with the cofactor provide one mechanism for linkage between binding of FAD and folate binding that could account in part for the protective action of folates. Conformation changes induced by folate binding may also suppress dissociation of FAD.     

The structure and properties of methylenetetrahydrofolate reductase from Escherichia coli suggest how folate ameliorates human hyperhomocysteinemia

Elevated plasma homocysteine levels are associated with increased risk for cardiovascular disease and neural tube defects in humans. Folate treatment decreases homocysteine levels and dramatically reduces the incidence of neural tube defects. The flavoprotein methylenetetrahydrofolate reductase (MTHFR) is a likely target for these actions of folate. The most common genetic cause of mildly elevated plasma homocysteine in humans is the MTHFR polymorphism A222V (base change C677-->T). The X-ray analysis of E. coli MTHFR, reported here, provides a model for the catalytic domain that is shared by all MTHFRs. This domain is a beta8alpha8 barrel that binds FAD in a novel fashion. Ala 177, corresponding to Ala 222 in human MTHFR, is near the bottom of the barrel and distant from the FAD. The mutation A177V does not affect Km or k(cat) but instead increases the propensity for bacterial MTHFR to lose its essential flavin cofactor. Folate derivatives protect wild-type and mutant E. coli enzymes against flavin loss, and protect human MTHFR and the A222V mutant against thermal inactivation, suggesting a mechanism by which folate treatment reduces homocysteine levels.     

A Structured-based Model for the Decreased Activity of Ala222Val and Glu429Ala Methylenetetrahydrofolate Reductase (MTHFR) Mutants

The structure of human Methylenetetrahydrofolate Reductase (MTHFR) is not known either by NMR or by X-ray methods. Phosphorylation seems to play an important role in the functioning of this flavoprotein. MTHFR catalyzes an irreversible reaction in homocysteine metabolism. Phosphorylation decreases the activity of MTHFR by enhancing the sensitivity of the enzyme to SAdenosylmethione. Two common polymorphisms in MTHFR, Ala222Val and Glu429Ala, can result in a number of vascular diseases. Effects of the Glu429Ala polymorphism on the structure of human MTHFR remain undetermined due to limited structural information. Hence, structural models of the MTHFR mutants were constructed using I-TASSER and assessed by PROCHECK, DFIRE and Verify3D tools. A mechanism is further suggested for the decreased activity of the Ala222Val and Glu429Ala mutants due to a decrease in number of serine phosphorylation sites using information gleaned from the molecular models. This provides insights for the understanding of structure-function relationship for MTHFR.     

Properties and crystal structure of methylenetetrahydrofolate reductase from Thermus thermophilus HB8

Background

Methylenetetrahydrofolate reductase (MTHFR) is one of the enzymes involved in homocysteine metabolism. Despite considerable genetic and clinical attention, the reaction mechanism and regulation of this enzyme are not fully understood because of difficult production and poor stability. While recombinant enzymes from thermophilic organisms are often stable and easy to prepare, properties of thermostable MTHFRs have not yet been reported.

Methodology/Principal Findings

MTHFR from Thermus thermophilus HB8, a homologue of Escherichia coli MetF, has been expressed in E. coli and purified. The purified MTHFR was chiefly obtained as a heterodimer of apo- and holo-subunits, that is, one flavin adenine dinucleotide (FAD) prosthetic group bound per dimer. The crystal structure of the holo-subunit was quite similar to the β₈α₈ barrel of E. coli MTHFR, while that of the apo-subunit was a previously unobserved closed form. In addition, the intersubunit interface of the dimer in the crystals was different from any of the subunit interfaces of the tetramer of E. coli MTHFR. Free FAD could be incorporated into the apo-subunit of the purified Thermus enzyme after purification, forming a homodimer of holo-subunits. Comparison of the crystal structures of the heterodimer and the homodimer revealed different intersubunit interfaces, indicating a large conformational change upon FAD binding. Most of the biochemical properties of the heterodimer and the homodimer were the same, except that the homodimer showed ≈50% activity per FAD-bound subunit in folate-dependent reactions.

Conclusions/Significance

The different intersubunit interfaces and rearrangement of subunits of Thermus MTHFR may be related to human enzyme properties, such as the allosteric regulation by S-adenosylmethionine and the enhanced instability of the Ala222Val mutant upon loss of FAD. Whereas E. coli MTHFR was the only structural model for human MTHFR to date, our findings suggest that Thermus MTHFR will be another useful model for this important enzyme.     

A proteomic approach for the characterization of C677T mutation of the human gene methylenetetrahydrofolate reductase

Methylenetetrahydrofolate reductase (MTHFR) catalyzes the conversion of methylenetetrahydrofolate (CH2H4folate) to methyltetrahydrofolate (CH3H4folate). The C677T mutation is a common polymorphism of the human enzyme that leads to the replacement of Ala222Val, thermolability of MTHFR, and mild elevation of plasma homocysteine levels. A mild hyperhomocysteinemia is known to be risk factor for cardiovascular and thrombotic diseases, ischemic stroke, neural tube defects, late on-set dementia, and pregnancy complications. Human plasma of subjects carrying the C677T mutation in the MTHFR gene has been investigated for their protein pattern in order to identify novel molecular hallmarks. 2-D analysis of the plasma protein allowed the identification of a specific pattern associated with the TT mutant genotype. Noteworthy, we found one spot shifted to a more basic pI in mutant individuals, and MS identification corresponded to vitamin D-binding protein (DBP or group component (Gc) globulin). MS/MS peptide sequencing allowed to discriminate different allelic variants in the investigated clinical groups. These data confirmed by molecular genetic analysis highlight the novel association between the C677T MTHFR genotype with the Gc2 polymorphism of the DBP. Moreover, we found a quantitative reduction of Apolipoprotein A-I in mutant individuals, which was associated, in previous studies by others to an increased cardiovascular risk.     

Molecular dynamics-based analyses of the structural instability and secondary structure of the fibrinogen gamma chain protein with the D356V mutation

Mutations in the fibrinogen gamma chain (FGG) gene have been associated with various disorders, such as dysfibrinogenemia, thrombophilia, and hypofibrinogenemia. A literature survey showed that a residue exchange in fibrinogen Milano I from γ Asp to Val at position 330 impairs fibrin polymerization. The D356V (D330V) mutation located in the C-terminus was predicted to be highly deleterious and to affect the function of the protein. The pathogenicity of the altered gene and changes in protein functions were predicted using in silico methods, such as SIFT, PolyPhen 2, I-Mutant 3.0, Align GV–GD, PhD–SNP, and SNPs&GO. The secondary structure of the mutant protein was unwound by the end of the 50-ns simulation period, and a structural change in the helix-turn transition of the alpha-helical (352–356) region residues was observed. Moreover, a change in the length of the helical region was visualized in the mutant trajectory file, indicating the local transient unfolding of the protein. The obtained computational results suggest that the substitution of the neutral amino acid valine for the acidic amino acid aspartic acid at position 356 results in an unwound conformation within 50 ns, which might contribute to defective polymerization. Our analysis also provides insights into the effect of the conformational change in the D356V (D330V) mutant on protein structure and function.     

Activity of human dihydrolipoamide dehydrogenase is reduced by mutation at threonine-44 of FAD-binding region to valine

Dihydrolipoamide dehydrogenase (E3) is a member of the pyridine nucleotide-disulfide oxidoreductase family. Thr residues are highly conserved. They are at the active site disulfide-bond regions of most E3s and other oxidoreductases. The crystal structure of Azotobacter vinelandii E3 suggests that the hydroxyl group of Thr that are involved in the FAD binding interact with the adenosine phosphate of FAD. However, several prokaryotic E3s have Val instead of Thr. To investigate the meaning and importance of the Thr conservation in many E3s, the corresponding residue, Thr-44, in human E3 was substituted to Val by site-directed mutagenesis. The mutant's E3 activity showed about a 2.2-fold decrease. Its UV-visible and fluorescence spectra indicated that the mutant might have a slightly different microenvironment at the FAD-binding region.     

Role of surface hydrophobic residues in the conformational stability of human lysozyme at three different positions

To evaluate the contribution of the amino acid residues on the surface of a protein to its stability, a series of hydrophobic mutant human lysozymes (Val to Gly, Ala, Leu, Ile, Met, and Phe) modified at three different positions on the surface, which are located in the alpha-helix (Val 110), the beta-sheet (Val 2), and the loop (Val 74), were constructed. Their thermodynamic parameters of denaturation and crystal structures were examined by calorimetry and by X-ray crystallography at 100 K, respectively. Differences in the denaturation Gibbs energy change between the wild-type and the hydrophobic mutant proteins ranged from 4.6 to -9.6 kJ/mol, 2.7 to -1.5 kJ/mol, and 3.6 to -0.2 kJ/mol at positions 2, 74, and 110, respectively. The identical substitution at different positions and different substitutions at the same position resulted in different degrees of stabilization. Changes in the stability of the mutant proteins could be evaluated by a unique equation considering the conformational changes due to the substitutions [Funahashi et al. (1999) Protein Eng. 12, 841-850]. For this calculation, secondary structural propensities were newly considered. However, some mutant proteins were not adapted to the equation. The hydration structures around the mutation sites of the exceptional mutant proteins were affected due to the substitutions. The stability changes in the exceptional mutant proteins could be explained by the formation or destruction of the hydration structures. These results suggest that the hydration structure mediated via hydrogen bonds covering the protein surface plays an important role in the conformational stability of the protein.     

The role of the Val57 amino‐acid residue in the hinge loop of the human cystatin C. Conformational studies of the beta2‐L1‐beta3 segments of wild‐type human cystatin C and its mutants

Human cystatin C (HCC) is one of the amyloidogenic proteins to be shown to oligomerize via a three‐dimensional domain swapping mechanism. This process precedes the formation of a stable dimer and proceeds particularly easily in the case of the L68Q mutant. According to the proposed mechanism, dimerization of the HCC precedes conformational changes within the β2 and β3 strands. In this article, we present conformational studies, using circular dichroism and MD methods, of the β2‐L1‐β3 (His43‐Thr72) fragment of the HCC involved in HCC dimer formation. We also carried out studies of the β2‐L1‐β3 peptide, in which the Val57 residue was replaced by residues promoting β‐turn structure formation (Asp, Asn, or Pro). The present study established that point mutation could modify the structure of the L1 loop in the β‐hairpin peptide. Our results showed that the L1 loop in the peptide excised from human cystatin C is broader than that in cystatin C. In the HCC protein, broadening of the L1 loop together with the unfavorable L68Q mutation in the hydrophobic pocket could be a force sufficient to cause the partial unfolding and then the opening of HCC or its L68Q mutant structure for further dimerization. We presume further that the Asp57 and Asn57 mutations in the L1 loop of HCC could stabilize the closed form of HCC, whereas the Pro57 mutation could lead to the opening of the HCC structure and then to dimer/oligomer formation. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 373–383, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Role of tyrosine hot‐spot residues at the interface of colicin E9 and immunity protein 9: A comparative free energy simulation study

The endonuclease activity of the bacterial colicin 9 enzyme is controlled by the specific and high‐affinity binding of immunity protein 9 (Im9). Molecular dynamics simulation studies in explicit solvent were used to investigate the free energy change associated with the mutation of two hot‐spot interface residues [tyrosine (Tyr): Tyr54 and Tyr55] of Im9 to Ala. In addition, the effect of several other mutations (Leu33Ala, Leu52Ala, Val34Ala, Val37Ala, Ser48Ala, and Ile53Ala) with smaller influence on binding affinity was also studied. Good qualitative agreement of calculated free energy changes and experimental data on binding affinity of the mutations was observed. The simulation studies can help to elucidate the molecular details on how the mutations influence protein–protein binding affinity. The role of solvent and conformational flexibility of the partner proteins was studied by comparing the results in the presence or absence of solvent and with or without positional restraints. Restriction of the conformational mobility of protein partners resulted in significant changes of the calculated free energies but of similar magnitude for isolated Im9 and for the complex and therefore in only modest changes of binding free energy differences. Although the overall binding free energy change was similar for the two Tyr–Ala mutations, the physical origin appeared to be different with solvation changes contributing significantly to the Tyr55Ala mutation and to a loss of direct protein–protein interactions dominating the free energy change due to the Tyr54Ala mutation. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

MTHFR 677 CT/MTHFR 1298 CC genotypes are associated with increased risk of hypertension in Indians

The goals of our present study were to measure plasma homocysteine levels and determine their association with methylenetetrahydrofolate reductase (MTHFR) gene polymorphisms (C677T and A1298C) in essential hypertensive subjects. Plasma total homocysteine and folic acid levels were measured in essential hypertensive patients (n = 153) before and after oral supplementation with either 5 mg folic acid tablet/day or 5 mg placebo/day for 4 weeks and compared with age and sex matched normotensive controls (n = 133). MTHFR gene polymorphisms (C677T and A1298C) were studied by restriction fragment length polymorphism and correlated with plasma homocysteine levels. Homocysteine levels were significantly higher in hypertensive patients as compared to controls and showed a negative correlation with plasma folate levels. Folic acid supplementation (5 mg/day) for 4 weeks resulted in a significant decrease in plasma homocysteine concentrations in these patients. Patients carrying MTHFR 677T allele (OR = 1.90; 95%CI: 1.14–3.19) or MTHFR 1298C (OR = 2.6, 95%CI: 1.55–4.40) allele were at increased risk of hypertension. The frequency of co-occurrence of MTHFR 677 CT/1298 CC genotypes was significantly higher in the patients compared to controls (P < 0.05) and was associated with increased risk of hypertension (OR = 3.54, 95%CI: 0.37–4.30). Subjects with MTHFR 1298 CC genotype had significantly higher homocysteine levels compared to those with MTHFR 1298 AA genotype (P < 0.05). Our results indicate that MTHFR 677T and 1298C alleles and co-occurrence of MTHFR 677 CT/MTHFR 1298 CC genotypes are associated with increased risk of hypertension and MTHFR 1298 CC genotype is associated with higher homocysteine levels in our subjects.     

Molecular genetics of MTHFR: polymorphisms are not all benign

Methylenetetrahydrofolate reductase (MTHFR) is a key regulatory enzyme in folate and homocysteine metabolism. Research performed during the past decade has clarified our understanding of MTHFR deficiencies that cause homocystinuria or mild hyperhomocysteinemia. Our cloning of the MTHFR coding sequence was initially followed by the identification of the first deleterious mutations in MTHFR, in patients with homocystinuria and marked hyperhomocysteinemia. Shortly thereafter, we identified the 677C-->T variant and showed that it encoded a thermolabile enzyme with reduced activity. Currently, a total of 41 rare but deleterious mutations in MTHFR, as well as about 60 polymorphisms have been reported. The 677C-->T (Ala222Val) variant has been particularly noteworthy since it has become recognized as the most common genetic cause of hyperhomocysteinemia. The disruption of homocysteine metabolism by this polymorphism influences risk for several complex disorders, including cardiovascular disease, neural tube defects and some cancers. We describe here the complex structure of the MTHFR gene, summarize the current state of knowledge on rare and common mutations in MTHFR and discuss some relevant findings in a mouse model for MTHFR deficiency.     

Engineering stability in NADPH oxidases: A common strategy for enzyme production

NADPH oxidases (NOXs) are membrane enzymes whose sole function is the generation of reactive oxygen species. Humans have seven NOX isoenzymes that feature distinct functions in immune response and cell signaling but share the same catalytic core comprising a FAD-binding dehydrogenase domain and a heme-binding transmembrane domain. We previously described a mutation that stabilizes the dehydrogenase domain of a prokaryotic homolog of human NOX5. The thermostable mutant exhibited a large 19?°C increase in the apparent melting temperature (app T_m) and a much tighter binding of the FAD cofactor, which allowed the crystallization and structure determination of the domain holo-form. Here, we analyze the transferability of this mutation onto prokaryotic and eukaryotic full-length NOX enzymes. We found that the mutation exerts a significative stabilizing effect on the full-length NOX5 from both Cylindrospermum stagnale (app T_m increase of 8?°C) and Homo sapiens (app ΔT_m of 2?°C). Enhanced thermal stability resulted in more homogeneous preparations of the bacterial NOX5 with less aggregation problems. Moreover, we also found that the mutation increases the overall expression of recombinant human NOX4 and NOX5 in mammalian cells. Such a 2–5-fold increase is mainly due to the lowered cell toxicity, which leads to higher biomasses. Because of the high sequence identity of the catalytic core within this family of enzymes, this strategy can be a general tool to boost the production of all NOXs.     

Organization of the multiple coenzymes and subunits and role of the covalent flavin link in the complex heterotetrameric sarcosine oxidase

Heterotetrameric (alphabetagammadelta) sarcosine oxidase from Corynebacterium sp. P-1 (cTSOX) contains noncovalently bound FAD and NAD(+) and covalently bound FMN, attached to beta(His173). The beta(His173Asn) mutant is expressed as a catalytically inactive, labile heterotetramer. The beta and delta subunits are lost during mutant enzyme purification, which yields a stable alphagamma complex. Addition of stabilizing agents prevents loss of the delta but not the beta subunit. The covalent flavin link is clearly a critical structural element and essential for TSOX activity or preventing FMN loss. The alpha subunit was expressed by itself and purified by affinity chromatography. The alpha and beta subunits each contain an NH(2)-terminal ADP-binding motif that could serve as part of the binding site for NAD(+) or FAD. The alpha subunit and the alphagamma complex were each found to contain 1 mol of NAD(+) but no FAD. Since NAD(+) binds to alpha, FAD probably binds to beta. The latter could not be directly demonstrated since it was not possible to express beta by itself. However, FAD in TSOX from Pseudomonas maltophilia (pTSOX) exhibits properties similar to those observed for the covalently bound FAD in monomeric sarcosine oxidase and N-methyltryptophan oxidase, enzymes that exhibit sequence homology with beta. A highly conserved glycine in the ADP-binding motif of the alpha(Gly139) or beta(Gly30) subunit was mutated in an attempt to generate NAD(+)- or FAD-free cTSOX, respectively. The alpha(Gly139Ala) mutant is expressed only at low temperature (t(optimum) = 15 degrees C), but the purified enzyme exhibited properties indistinguishable from the wild-type enzyme. The much larger barrier to NAD(+) binding in the case of the alpha(Gly139Val) mutant could not be overcome even by growth at 3 degrees C, suggesting that NAD(+) binding is required for TSOX expression. The beta(Gly30Ala) mutant exhibited subunit expression levels similar to those of the wild-type enzyme, but the mutation blocked subunit assembly and covalent attachment of FMN, suggesting that both processes require a conformational change in beta that is induced upon FAD binding. About half of the covalent FMN in recombinant preparations of cTSOX or pTSOX is present as a reversible covalent 4a-adduct with a cysteine residue. Adduct formation is not prevented by mutating any of the three cysteine residues in the beta subunit of cTSOX to Ser or Ala. Since FMN is attached via its 8-methyl group to the beta subunit, the FMN ring must be located at the interface between beta and another subunit that contains the reactive cysteine residue.     

Biosynthesis of Biotin-vitamers from Oleic Acid by Cryptococcus neoformans

Chitosanase (ChoA) from Mitsuaria chitosanitabida 3001 was successfully evolved with secretion efficiency and thermal stability. The inactive ChoA mutant (G151D) gene was used to mutate by an error-prone PCR technique and mutant genes that restored chitosanase activity were isolated. Two desirable mutants, designated M5S and M7T, were isolated. Two amino acids, Leu74 and Val75, in the signal peptide of ChoA were changed to Gln and Ile respectively in the M7T mutant, in addition to the G151D mutation. The L74Q/V75I double ChoA mutant was 1.5-fold higher in specific activity than wild-type ChoA due to efficient secretion of ChoA. One amino acid Asn222 was changed to Ser in the M5S mutant in addition to the G151D mutation. The N222S single ChoA mutant was 1.2-fold higher in specific activity and showed a 17% increase in thermal stability at 50 °C as compared with wild-type ChoA. This is the first study to achieve an evolutional increase in enzyme capability among chitosanses.     

Protein stability and mutations in the axial methionine loop of a minimal cytochrome<Emphasis Type="Italic"> c</Emphasis>

The minimal mono-heme ferricytochrome c from Bacillus pasteurii, containing 71 amino acids, has been further investigated through mutagenesis of different positions in the loop containing the iron ligand Met71. These mutations have been designed to sample different aspects of the loop structure, in order to obtain insights into the determinants of the stability of the iron(III) environment. In particular, positions 68, 72 and 75 have been essayed. Gln68 has been mutated to Lys to provide a suitable alternate ligand that can displace Met71 under denaturing conditions. Pro72 has been mutated to Gly and Ala to modify the range of allowed backbone conformations. Ile75, which is in van der Waals contact with Met71 and partly shields a long-lived water molecule in a protein cavity, has been substituted by Val and Ala to affect the network of inter-residue interactions around the metal site. The different contributions of the above amino acids to protein parameters such as structure, redox potential and the overall stability against unfolding with guanidinium hydrochloride are analyzed. While the structure remains essentially the same, the stability decreases with mutations. The comparison with mitochondrial c-type cytochromes is instructive.Abbreviations Bpcytc soluble fragment of cytochrome c₅₅₃ from Bacillus pasteurii - GdmCl guanidinium chloride - I75A Ile75 to Ala mutant - I75V Ile75 to Val mutant - P72A Pro72 to Ala mutant - P72G Pro72 to Gly mutant - Q68K Gln75 to Lys mutant - WT wild type     

Altered Activation of Protein Kinase PKR and Enhanced Apoptosis in Dystonia Cells Carrying a Mutation in PKR Activator Protein PACT

PACT is a stress-modulated activator of the interferon-induced double-stranded RNA-activated protein kinase (PKR). Stress-induced phosphorylation of PACT is essential for PACT''s association with PKR leading to PKR activation. PKR activation leads to phosphorylation of translation initiation factor eIF2α inhibition of protein synthesis and apoptosis. A recessively inherited form of early-onset dystonia DYT16 has been recently identified to arise due to a homozygous missense mutation P222L in PACT. To examine if the mutant P222L protein alters the stress-response pathway, we examined the ability of mutant P222L to interact with and activate PKR. Our results indicate that the substitution mutant P222L activates PKR more robustly and for longer duration albeit with slower kinetics in response to the endoplasmic reticulum stress. In addition, the affinity of PACT-PACT and PACT-PKR interactions is enhanced in dystonia patient lymphoblasts, thereby leading to intensified PKR activation and enhanced cellular death. P222L mutation also changes the affinity of PACT-TRBP interaction after cellular stress, thereby offering a mechanism for the delayed PKR activation in response to stress. Our results demonstrate the impact of a dystonia-causing substitution mutation on stress-induced cellular apoptosis.     

Aspartate 120 of Escherichia coli methylenetetrahydrofolate reductase: evidence for major roles in folate binding and catalysis and a minor role in flavin reactivity

Escherichia coli methylenetetrahydrofolate reductase (MTHFR) catalyzes the NADH-linked reduction of 5,10-methylenetetrahydrofolate (CH(2)-H(4)folate) to 5-methyltetrahydrofolate (CH(3)-H(4)folate) using flavin adenine dinucleotide (FAD) as cofactor. MTHFR is unusual among flavin oxidoreductases because it contains a conserved, negatively rather than positively charged amino acid (aspartate 120) near the N1-C2=O position of the flavin. At this location, Asp 120 is expected to influence the redox properties of the enzyme-bound FAD. Modeling of the CH(3)-H(4)folate product into the enzyme active site suggests that Asp 120 may also play crucial roles in folate binding and catalysis. We have replaced Asp 120 with Asn, Ser, Ala, Val, and Lys and have characterized the mutant enzymes. Consistent with a loss of negative charge near the flavin, the midpoint potentials of the mutants increased from 17 to 30 mV. A small kinetic effect on the NADH reductive half-reaction was also observed as the mutants exhibited a 1.2-1.5-fold faster reduction rate than the wild-type enzyme. Catalytic efficiency (k(cat)/K(m)) in the CH(2)-H(4)folate oxidative half-reaction was decreased significantly (up to 70000-fold) and in a manner generally consistent with the negative charge density of position 120, supporting a major role for Asp 120 in electrostatic stabilization of the putative 5-iminium cation intermediate during catalysis. Asp 120 is also intimately involved in folate binding as increases in the apparent K(d) of up to 15-fold were obtained for the mutants. Examining the E(red) + CH(2)-H(4)folate reaction at 4 degrees C, we obtained, for the first time, evidence for the rapid formation of a reduced enzyme-folate complex with wild-type MTHFR. The more active Asp120Ala mutant, but not the severely impaired Asp120Lys mutant, demonstrated the species, suggesting a connection between the extent of complex formation and catalytic efficiency.     

XLOS-Observed Mutations of MID1 Bbox1 Domain Cause Domain Unfolding

MID1 catalyzes the ubiquitination of the protein alpha4 and the catalytic subunit of protein phosphatase 2A. Mutations within the MID1 Bbox1 domain are associated with X-linked Opitz G syndrome (XLOS). Our functional assays have shown that mutations of Ala130 to Val or Thr, Cys142 to Ser and Cys145 to Thr completely disrupt the polyubiquitination of alpha4. Using NMR spectroscopy, we characterize the effect of these mutations on the tertiary structure of the Bbox1 domain by itself and in tandem with the Bbox2 domain. The mutation of either Cys142 or Cys145, each of which is involved in coordinating one of the two zinc ions, results in the collapse of signal dispersion in the HSQC spectrum of the Bbox1 domain indicating that the mutant protein structure is unfolded. Each mutation caused the coordination of both zinc ions, which are ∼13 Å apart, to be lost. Although Ala130 is not involved in the coordination of a zinc ion, the Ala130Thr mutant Bbox1 domain yields a poorly dispersed HSQC spectrum similar to those of the Cys142Ser and Cys145Thr mutants. Interestingly, neither cysteine mutation affects the structure of the adjacent Bbox2 domain when the two Bbox domains are engineered in their native tandem Bbox1-Bbox2 protein construct. Dynamic light scattering measurements suggest that the mutant Bbox1 domain has an increased propensity to form aggregates compared to the wild type Bbox1 domain. These studies provide insight into the mechanism by which mutations observed in XLOS affect the structure and function of the MID1 Bbox1 domain.