Diagnostic chemical shift markers for loop conformation and substrate and cofactor binding in dihydrofolate reductase complexes

Heteronuclear NMR methods have been used to probe the conformation of four complexes of Escherichia coli dihydrofolate reductase (DHFR) in solution. (1)H(N), (15)N, and (13)C(alpha) resonance assignments have been made for the ternary complex with folate and oxidized NADP(+) cofactor and the ternary complex with folate and a reduced cofactor analog, 5,6-dihydroNADPH. The backbone chemical shifts have been compared with those of the binary complex of DHFR with the substrate analog folate and the binary complex with NADPH (the holoenzyme). Analysis of (1)H(N) and (15)N chemical shifts has led to the identification of marker resonances that report on the active site conformation of the enzyme. Other backbone amide resonances report on the presence of ligands in the pterin binding pocket and in the adenosine and nicotinamide-ribose binding sites of the NADPH cofactor. The chemical shift data indicate that the enzyme populates two dominant structural states in solution, with the active site loops in either the closed or occluded conformations defined by X-ray crystallography; there is no evidence that the open conformation observed in some X-ray structures of E. coli DHFR are populated in solution.     

Backbone dynamics in dihydrofolate reductase complexes: role of loop flexibility in the catalytic mechanism

To elucidate the influence of local motion of the polypeptide chain on the catalytic mechanism of an enzyme, we have measured (15)N relaxation data for Escherichia coli dihydrofolate reductase in three different complexes, representing different stages in the catalytic cycle of the enzyme. NMR relaxation data were analyzed by the model-free approach, corrected for rotational anisotropy, to provide insights into the backbone dynamics. There are significant differences in the backbone dynamics in the different complexes. Complexes in which the cofactor binding site is occluded by the Met20 loop display large amplitude motions on the picosecond/nanosecond time scale for residues in the Met20 loop, the adjacent betaF-betaG loop and for residues 67-69 in the adenosine binding loop. Formation of the closed Met20 loop conformation in the ternary complex with folate and NADP(+), results in attenuation of the motions in the Met20 loop and the betaF-betaG loop but leads to increased flexibility in the adenosine binding loop. New fluctuations on a microsecond/millisecond time scale are observed in the closed E:folate:NADP(+) complex in regions that form hydrogen bonds between the Met20 and the betaF-betaG loops. The data provide insights into the changes in backbone dynamics during the catalytic cycle and point to an important role of the Met20 and betaF-betaG loops in controlling access to the active site. The high flexibility of these loops in the occluded conformation is expected to promote tetrahydrofolate-assisted product release and facilitate binding of the nicotinamide ring to form the Michaelis complex. The backbone fluctuations in the Met20 loop become attenuated once it closes over the active site, thereby stabilizing the nicotinamide ring in a geometry conducive to hydride transfer. Finally, the relaxation data provide evidence for long-range motional coupling between the adenosine binding loop and distant regions of the protein.     

Conformational changes in the active site loops of dihydrofolate reductase during the catalytic cycle

Escherichia coli dihydrofolate reductase (DHFR) has several flexible loops surrounding the active site that play a functional role in substrate and cofactor binding and in catalysis. We have used heteronuclear NMR methods to probe the loop conformations in solution in complexes of DHFR formed during the catalytic cycle. To facilitate the NMR analysis, the enzyme was labeled selectively with [(15)N]alanine. The 13 alanine resonances provide a fingerprint of the protein structure and report on the active site loop conformations and binding of substrate, product, and cofactor. Spectra were recorded for binary and ternary complexes of wild-type DHFR bound to the substrate dihydrofolate (DHF), the product tetrahydrofolate (THF), the pseudosubstrate folate, reduced and oxidized NADPH cofactor, and the inactive cofactor analogue 5,6-dihydroNADPH. The data show that DHFR exists in solution in two dominant conformational states, with the active site loops adopting conformations that closely approximate the occluded or closed conformations identified in earlier X-ray crystallographic analyses. A minor population of a third conformer of unknown structure was observed for the apoenzyme and for the disordered binary complex with 5,6-dihydroNADPH. The reactive Michaelis complex, with both DHF and NADPH bound to the enzyme, could not be studied directly but was modeled by the ternary folate:NADP(+) and dihydrofolate:NADP(+) complexes. From the NMR data, we are able to characterize the active site loop conformation and the occupancy of the substrate and cofactor binding sites in all intermediates formed in the extended catalytic cycle. In the dominant kinetic pathway under steady-state conditions, only the holoenzyme (the binary NADPH complex) and the Michaelis complex adopt the closed loop conformation, and all product complexes are occluded. The catalytic cycle thus involves obligatory conformational transitions between the closed and occluded states. Parallel studies on the catalytically impaired G121V mutant DHFR show that formation of the closed state, in which the nicotinamide ring of the cofactor is inserted into the active site, is energetically disfavored. The G121V mutation, at a position distant from the active site, interferes with coupled loop movements and appears to impair catalysis by destabilizing the closed Michaelis complex and introducing an extra step into the kinetic pathway.     

Interligand Overhauser effects in type II dihydrofolate reductase

R67 dihydrofolate reductase (DHFR) is a type II DHFR produced by bacteria as a resistance mechanism to the increased clinical use of the antibacterial drug trimethoprim. Type II DHFRs are not homologous in either sequence or structure with chromosomal DHFRs. The type II enzymes contain four identical subunits which form a homotetramer containing a single active site pore accessible from either end. Although the crystal structure of the complex of R67 DHFR with folate has been reported [Narayana et al. (1995) Nat. Struct. Biol. 2, 1018], the nature of the ternary complex which must form with substrate and cofactor is unclear. We have performed transferred NOE and interligand NOE (ILOE) studies to analyze the ternary complexes formed from NADP(+) and folate in order to probe the structure of the ternary complex. Consistent with previous studies of the binary complex formed from another type II DHFR, the ribonicotinamide bond of NADP(+) was found to adopt a syn conformation, while the adenosine moiety adopts an anti conformation. Large ILOE peaks connecting NADP(+) H4 and H5 with folate H9 protons are observed, while the absence of a large ILOE connecting NADP(+) H4 and H5 with folate H7 indicates that the relative orientation of the two ligands differs significantly from the orientation in the chromosomal enzyme. To obtain more detailed insight, we prepared and studied the folate analogue 2-deamino-2-methyl-5,8-dideazafolate (DMDDF) which contains additional protons in order to provide additional NOEs. For this analogue, the exchange characteristics of the corresponding ternary complex were considerably poorer, and it was necessary to utilize higher enzyme concentrations and higher temperature in order to obtain ILOE information. The results support a structure in which the NADP(+) and folate/DMDDF molecules extend in opposite directions parallel to the long axis of the pore, with the nicotinamide and pterin ring systems approximately stacked at the center. Such a structure leads to a ternary complex which is in many respects similar to the gas-phase theoretical calculations of the dihydrofolate-NADPH transition state by Andres et al. [(1996) Bioorg. Chem. 24, 10-18]. Analogous NMR studies performed on folate, DMDDF, and R67 DHFR indicate formation of a ternary complex in which two symmetry-related binding sites are occupied by folate and DMDDF.     

Crystal structure of chicken liver dihydrofolate reductase complexed with NADP+ and biopterin.

The 2.2-A crystal structure of chicken liver dihydrofolate reductase (EC 1.5.1.3, DHFR) has been solved as a ternary complex with NADP+ and biopterin (a poor substrate). The space group and unit cell are isomorphous with the previously reported structure of chicken liver DHFR complexed with NADPH and phenyltriazine [Volz, K. W., Matthews, D. A., Alden, R. A., Freer, S. T., Hansch, C., Kaufman, B. T., & Kraut, J. (1982) J. Biol. Chem. 257, 2528-2536]. The structure contains an ordered water molecule hydrogen-bonded to both hydroxyls of the biopterin dihydroxypropyl group as well as to O4 and N5 of the biopterin pteridine ring. This water molecule, not observed in previously determined DHFR structures, is positioned to complete a proposed route for proton transfer from the side-chain carboxylate of E30 to N5 of the pteridine ring. Protonation of N5 is believed to occur during the reduction of dihydropteridine substrates. The positions of the NADP+ nicotinamide and biopterin pteridine rings are quite similar to the nicotinamide and pteridine ring positions in the Escherichia coli DHFR.NADP+.folate complex [Bystroff, C., Oatley, S. J., & Kraut, J. (1990) Biochemistry 29, 3263-3277], suggesting that the reduction of biopterin and the reduction of folate occur via similar mechanisms, that the binding geometry of the nicotinamide and pteridine rings is conserved between DHFR species, and that the p-aminobenzoylglutamate moiety of folate is not required for correct positioning of the pteridine ring in ground-state ternary complexes. Instead, binding of the p-aminobenzoylglutamate moiety of folate may induce the side chain of residue 31 (tyrosine or phenylalanine) in vertebrate DHFRs to adopt a conformation in which the opening to the pteridine binding site is too narrow to allow the substrate to diffuse away rapidly. A reverse conformational change of residue 31 is proposed to be required for tetrahydrofolate release.     

NMR studies of the interaction of a type II dihydrofolate reductase with pyridine nucleotides reveal unexpected phosphatase and reductase activity

The interaction of type II R67 dihydrofolate reductase (DHFR) with its cofactor nicotinamide adenine dinucleotide phosphate (NADP(+)) has been studied using nuclear magnetic resonance (NMR). Doubly labeled [U-(13)C,(15)N]DHFR was obtained from Escherichia coli grown on a medium containing [U-(13)C]-D-glucose and (15)NH(4)Cl, and the 16 disordered N-terminal amino acids were removed by treatment with chymotrypsin. Backbone and side chain NMR assignments were made using triple-resonance experiments. The degeneracy of the amide (1)H and (15)N shifts of the tetrameric DHFR was preserved upon addition of NADP(+), consistent with kinetic averaging among equivalent binding sites. Analysis of the more titration-sensitive DHFR amide resonances as a function of added NADP(+) gave a K(D) of 131 +/- 50 microM, consistent with previous determinations using other methodology. We have found that the (1)H spectrum of NADP(+) in the presence of the R67 DHFR changes as a function of time. Comparison with standard samples and mass spectrometric analysis indicates a slow conversion of NADP(+) to NAD(+), i.e., an apparent NADP(+) phosphatase activity. Studies of this activity in the presence of folate and a folate analogue support the conclusion that this activity results from an interaction with the DHFR rather than a contaminating phosphatase. (1)H NMR studies of a mixture of NADP(+) and NADPH in the presence of the enzyme reveal that a ternary complex forms in which the N-4A and N-4B nuclei of the NADPH are in the proximity of the N-4 and N-5 nuclei of NADP(+). Studies using the NADP(+) analogue acetylpyridine adenosine dinucleotide phosphate (APADP(+)) demonstrated a low level of enzyme-catalyzed hydride transfer from NADPH. Analysis of DHFR backbone dynamics revealed little change upon binding of NADP(+). These additional catalytic activities and dynamic behavior are in marked contrast to those of type I DHFR.     

Structure of the Q67H mutant of R67 dihydrofolate reductase-NADP+ complex reveals a novel cofactor binding mode

Plasmid-encoded bacterial R67 dihydrofolate reductase (DHFR) is a NADPH-dependent enzyme unrelated to chromosomal DHFR in amino acid sequence and structure. R67 DHFR is insensitive to the bacterial drug trimethoprim in contrast to chromosomal DHFR. The crystal structure of Q67H mutant of R67 DHFR bound to NADP(+) has been determined at 1.15 angstroms resolution. The cofactor assumes an extended conformation with the nicotinamide ring bound near the center of the active site pore, the ribose and pyrophosphate group (PP(i)) extending toward the outer pore. The ribonicotinamide exhibits anti conformation as in chromosomal DHFR complexes. The relative orientation between the PP(i) and the nicotinamide ribose differs from that observed in chromosomal DHFR-NADP(+) complexes. The coenzyme displays symmetrical binding mode with several water-mediated hydrogen bonds with the protein besides ionic, stacking, and van der Waals interactions. The structure provides a molecular basis for the observed stoichiometry and cooperativity in ligand binding. The ternary model based on the present structure and the previous R67 DHFR-folate complex provides insight into the catalytic mechanism and indicates that the relative orientation of the reactants in plasmid DHFR is different from that seen in chromosomal DHFRs.     

Ligand-induced conformational changes in the crystal structures of Pneumocystis carinii dihydrofolate reductase complexes with folate and NADP+

Structural data from two independent crystal forms (P212121 and P21) of the folate (FA) binary complex and from the ternary complex with the oxidized coenzyme, NADP+, and recombinant Pneumocystis carinii dihydrofolate reductase (pcDHFR) refined to an average of 2.15 A resolution, show the first evidence of ligand-induced conformational changes in the structure of pcDHFR. These data are also compared with the crystal structure of the ternary complex of methotrexate (MTX) with NADPH and pcDHFR in the monoclinic lattice with data to 2.5 A resolution. Comparison of the data for the FA binary complex of pcDHFR with those for the ternary structures reveals significant differences, with a >7 A movement of the loop region near residue 23 that results in a new "flap-open" position for the binary complex, and a "closed" position in the ternary complexes, similar to that reported for Escherichia coli (ec) DHFR complexes. In the orthorhombic lattice for the binary FA pcDHFR complex, there is also an unwinding of a short helical region near residue 47 that places hydrophobic residues Phe-46 and Phe-49 toward the outer surface, a conformation that is stabilized by intermolecular packing contacts. The pyrophosphate moiety of NADP+ in the ternary folate pcDHFR complexes shows significant differences in conformation compared with that observed in the MTX-NADPH-pcDHFR ternary complex. Additionally, comparison of the conformations among these four pcDHFR structures reveals evidence for subdomain movement that correlates with cofactor binding states. The larger binding site access in the new "flap-open" loop 23 conformation of the binary FA complex is consistent with the rapid release of cofactor from the product complex during catalysis as well as the more rapid release of substrate product from the binary complex as a result of the weaker contacts of the closed loop 23 conformation, compared to ecDHFR.     

Role of water in the catalytic cycle of E. coli dihydrofolate reductase

Dihydrofolate reductase (DHFR) catalyzes the nicotinamide adenine dinucleotide phosphate (NADPH)-dependent reduction of 7,8-dihydrofolate (H2F) to 5,6,7,8-tetrahydrofolate (H4F). Because of the absence of any ionizable group in the vicinity of N5 of dihydrofolate it has been proposed that N5 could be protonated directly by a water molecule at the active site in the ternary complex of the Escherichia coli enzyme with cofactor and substrate. However, in the X-ray structures representing the Michaelis complex of the E. coli enzyme, a water molecule has never been observed in a position that could allow protonation of N5. In fact, the side chain of Met 20 blocks access to N5. Energy minimization reported here revealed that water could be placed in hydrogen bonding distance of N5 with only minor conformational changes. The r.m.s. deviation between the conformation of the M20 loop observed in the crystal structures of the ternary complexes and the conformation adopted after energy minimization was only 0.79 A. We performed molecular dynamics simulations to determine the accessibility by water of the active site of the Michaelis complex of DHFR. Water could access N5 relatively freely after an equilibration time of approximately 300 psec during which the side chain of Met 20 blocked water access. Protonation of N5 did not increase the accessibility by water. Surprisingly the number of near-attack conformations, in which the distance between the pro-R hydrogen of NADPH and C6 of dihydrofolate was less than 3.5 A and the angle between C4 and the pro-R hydrogen of NADPH and C6 of dihydrofolate was greater than 120 degrees, did not increase after protonation. However, when the hydride was transferred from NADPH to C6 of dihydrofolate before protonation, the side chain of Met 20 moved away from N5 after approximately 100 psec thereby providing water access. The average time during which water was found in hydrogen bonding distance to N5 was significantly increased. These results suggest that hydride transfer might occur early to midway through the reaction followed by protonation. Such a mechanism is supported by the very close contact between C4 of NADP+ and C6 of folate observed in the crystal structures of the ternary enzyme complexes, when the M20 loop is in its closed conformation.     

13C NMR studies of complexes of Escherichia coli dihydrofolate reductase formed with methotrexate and with folic acid.

13C NMR studies of 13C-labelled ligands bound to dihydrofolate reductase provide (DHFR) a powerful means of detecting and characterizing multiple bound conformations. Such studies of complexes of Escherichia coli DHFR with [4,7,8a,9-13C]- and [2,4a,6-13C]methotrexate (MTX) and [4,6,8a-13C]- and [2,4a,7,9-13C]folic acid confirm that in the binary complexes, MTX binds in two conformational forms and folate binds as a single conformation. Earlier studies on the corresponding complexes with Lactobacillus casei DHFR indicated that, in this case, MTX binds as a single conformation whereas folate binds in multiple conformational forms (both in its binary complex and ternary complex with NADP+); two of the bound conformational states for the folate complexes are very different from each other in that there is a 180 degrees difference in their pteridine ring orientation. In contrast, the two different conformational states observed for MTX bound to E. coli DHFR do not show such a major difference in ring orientation and bind with N1 protonated in both forms. The major difference appears to involve the manner in which the 4-NH2 group of MTX binds to the enzyme (although the same protein residues are probably involved in both interactions). Addition of either NADP+ or NADPH to the E. coli DHFR-MTX complex results in a single set of 13C signals for bound methotrexate consistent with only one conformational form in the ternary complexes.     

Mass spectrometry on hydrogen/deuterium exchange of dihydrofolate reductase: effects of ligand binding

To address the effects of ligand binding on the structural fluctuations of Escherichia coli dihydrofolate reductase (DHFR), the hydrogen/deuterium (H/D) exchange kinetics of its binary and ternary complexes formed with various ligands (folate, dihydrofolate, tetrahydrofolate, NADPH, NADP(+), and methotrexate) were examined using electrospray ionization mass spectrometry. The kinetic parameters of H/D exchange reactions, which consisted of two phases with fast and slow rates, were sensitively influenced by ligand binding, indicating that changes in the structural fluctuation of the DHFR molecule are associated with the alternating binding and release of the cofactor and substrate. No additivity was observed in the kinetic parameters between a ternary complex and its constitutive binary complexes, indicating that ligand binding cooperatively affects the structural fluctuation of the DHFR molecule via long-range interactions. The local H/D exchange profile of pepsin digestion fragments was determined by matrix-assisted laser desorption/ionization mass spectrometry, and the helix and loop regions that appear to participate in substrate binding, largely fluctuating in the apo-form, are dominantly influenced by ligand binding. These results demonstrate that the structural fluctuation of kinetic intermediates plays an important role in enzyme function, and that mass spectrometry on H/D exchange coupled with ligand binding and protease digestion provide new insight into the structure-fluctuation-function relationship of enzymes.     

Conformational relaxation following hydride transfer plays a limiting role in dihydrofolate reductase catalysis

The catalytic cycle of an enzyme is frequently associated with conformational changes that may limit maximum catalytic throughput. In Escherichia coli dihydrofolate reductase, release of the tetrahydrofolate (THF) product is the rate-determining step under physiological conditions and is associated with an "occluded" to "closed" conformational change. In this study, we demonstrate that in dihydrofolate reductase the closed to occluded conformational change in the product ternary complex (E.THF.NADP (+)) also gates progression through the catalytic cycle. Using NMR relaxation dispersion, we have measured the temperature and pH dependence of microsecond to millisecond time scale backbone dynamics of the occluded E.THF.NADP (+) complex. Our studies indicate the presence of three independent dynamic regions, associated with the active-site loops, the cofactor binding cleft, and the C-terminus and an adjacent loop, which fluctuate into discrete conformational substates with different kinetic and thermodynamic parameters. The dynamics of the C-terminally associated region is pH-dependent (p K a < 6), but the dynamics of the active-site loops and cofactor binding cleft are pH-independent. The active-site loop dynamics access a closed conformation, and the accompanying closed to occluded rate constant is comparable to the maximum pH-independent hydride transfer rate constant. Together, these results strongly suggest that the closed to occluded conformational transition in the product ternary complex is a prerequisite for progression through the catalytic cycle and that the rate of this process places an effective limit on the maximum rate of the hydride transfer step.     

Important role of the cys-191 cys-220 disulfide bond in thrombin function and allostery

Little is known on the role of disulfide bonds in the catalytic domain of serine proteases. The Cys-191-Cys-220 disulfide bond is located between the 190 strand leading to the oxyanion hole and the 220-loop that contributes to the architecture of the primary specificity pocket and the Na+ binding site in allosteric proteases. Removal of this bond in thrombin produces an approximately 100-fold loss of activity toward several chromogenic and natural substrates carrying Arg or Lys at P1. Na+ activation is compromised, and no fluorescence change can be detected in response to Na+ binding. A 1.54-A resolution structure of the C191A/C220A mutant in the free form reveals a conformation similar to the Na+-free slow form of wild type. The lack of disulfide bond exposes the side chain of Asp-189 to solvent, flips the backbone O atom of Gly-219, and generates disorder in portions of the 186 and 220 loops defining the Na+ site. This conformation, featuring perturbation of the Na+ site but with the active site accessible to substrate, offers a possible representation of the recently identified E* form of thrombin. Disorder in the 186 and 220 loops and the flip of Gly-219 are corrected by the active site inhibitor H-D-Phe-Pro-Arg-CH(2)Cl, as revealed by the 1.8-A resolution structure of the complex. We conclude that the Cys-191-Cys-220 disulfide bond confers stability to the primary specificity pocket by shielding Asp-189 from the solvent and orients the backbone O atom of Gly-219 for optimal substrate binding. In addition, the disulfide bond stabilizes the 186 and 220 loops that are critical for Na+ binding and activation.     

Advantages of fine-grained side chain conformer libraries

We compare the modelling accuracy of two common rotamer libraries, the Dunbrack-Cohen and the 'Penultimate' rotamer libraries, with that of a novel library of discrete side chain conformations extracted from the Protein Data Bank. These side chain conformer libraries are extracted automatically from high-quality protein structures using stringent filters and maintain crystallographic bond lengths and angles. This contrasts with traditional rotamer libraries defined in terms of chi angles under the assumption of idealized covalent geometry. We demonstrate that side chain modelling onto native and near-native main chain conformations is significantly more successful with the conformer libraries than with the rotamer libraries when solely considering excluded-volume interactions. The rotamer libraries are inadequate to model side chains without atomic clashes on over 20% of targets if the backbone is held fixed in the native conformation. An algorithm is described for simultaneously modelling both main chain and side chain atoms during discrete ab initio sampling. The resulting models have equivalent root mean square deviations from the experimentally determined protein loops as models from backbone-only ensembles, indicating that all-atom modelling does not detract from the accuracy of conformational sampling.     

A structural basis for the mechanism of aspartate-beta-semialdehyde dehydrogenase from Vibrio cholerae

L-Aspartate-beta-semialdehyde dehydrogenase (ASADH) catalyzes the reductive dephosphorylation of beta-aspartyl phosphate to L-aspartate-beta-semialdehyde in the aspartate biosynthetic pathway of plants and micro-organisms. The aspartate pathway produces fully one-quarter of the naturally occurring amino acids, but is not found in humans or other eukaryotic organisms, making ASADH an attractive target for the development of new antibacterial, fungicidal, or herbicidal compounds. We have determined the structure of ASADH from Vibrio cholerae in two states; the apoenzyme and a complex with NADP, and a covalently bound active site inhibitor, S-methyl-L-cysteine sulfoxide. Upon binding the inhibitor undergoes an enzyme-catalyzed reductive demethylation leading to a covalently bound cysteine that is observed in the complex structure. The enzyme is a functional homodimer, with extensive intersubunit contacts and a symmetrical 4-amino acid bridge linking the active site residues in adjacent subunits that could serve as a communication channel. The active site is essentially preformed, with minimal differences in active site conformation in the apoenzyme relative to the ternary inhibitor complex. The conformational changes that do occur result primarily from NADP binding, and are localized to the repositioning of two surface loops located on the rim at opposite sides of the NADP cleft.     

Crystal structures of Escherichia coli dihydrofolate reductase: the NADP+ holoenzyme and the folate.NADP+ ternary complex. Substrate binding and a model for the transition state

The crystal structure of dihydrofolate reductase (EC 1.5.1.3) from Escherichia coli has been solved as the binary complex with NADP+ (the holoenzyme) and as the ternary complex with NADP+ and folate. The Bragg law resolutions of the structures are 2.4 and 2.5 A, respectively. The new crystal forms are nonisomorphous with each other and with the methotrexate binary complex reported earlier [Bolin, J. T., Filman, D. J., Matthews, D. A., Hamlin, R. C., & Kraut, J. (1982) J. Biol. Chem. 257, 13650-13662]. In general, NADP+ and folate binding conform to predictions, but the nicotinamide moiety of NADP+ is disordered in the holoenzyme and ordered in the ternary complex. A mobile loop (residues 16-20) involved in binding the nicotinamide is also disordered in the holoenzyme. We report a detailed analysis of the binding interactions for both ligands, paying special attention to several apparently strained interactions that may favor the transition state for hydride transfer. Hypothetical models are presented for the binding of 7,8-dihydrofolate in the Michaelis complex and for the transition-state complex.     

Dismutation of dihydrofolate by dihydrofolate reductase

Degradation of 7,8-dihydrofolate (H2folate) in the presence of dihydrofolate reductase (DHFR) has been shown due not to an oxygenase activity of the reductase as previously reported but to dismutation of H2folate to folate and 5,6,7,8-tetrahydrofolate (H4folate). The reaction can be followed spectrophotometrically or by analysis of the reaction mixture by high-performance liquid chromatography (HPLC). The products have also been isolated and characterized. Oxygen uptake during the reaction is much less than stoichiometric with H2folate disappearance and is attributed to autoxidation of the H4folate formed. The dismutation activity is a property of highly purified Streptococcus faecium DHFR isoenzyme 2 (but not isoenzyme 1) and of Lactobacillus casei DHFR, but not of bovine liver DHFR. The activity is dependent on tightly bound NADP+ and/or NADPH. Removal of the nucleotide results in loss of dismutation activity, which is restored by adding NADP+ or NADPH. Maximum activity is obtained when approximately 1 mol equiv of nucleotide is added per mol of DHFR. It is proposed that in the dismutation reaction bound NADP(H) is alternately reduced and oxidized by incoming molecules of H2folate with release of folate and H4folate, respectively. The relatively slow rate of folate formation presumably limits the rate of the overall reaction. The equilibrium constant for the dismutation reaction is 19.4 +/- 7.4 at 22 degrees C and pH 7.0. Calculation of standard oxidation-reduction potentials at pH 7 gave values of -0.230 V for the H2folate/H4 folate pair and -0.268 V for the folate/H2folate pair. The mechanism by which NADP+ is retained by the enzyme from some sources during purification procedures is unclear.(ABSTRACT TRUNCATED AT 250 WORDS)     

Site-directed deletion mutants of a carboxyl-terminal region of human dihydrofolate reductase.

Substrate and inhibitor binding to dihydrofolate reductase (DHFR) primarily involves residues in the amino-terminal half of the enzyme; however, antibody binding studies performed in this laboratory suggested that the loop region located in the carboxyl terminus of human DHFR (hDHFR; residues 140-186) is involved in conformational changes that occur upon ligand binding and affect enzyme function (Ratnam, M., Tan, X., Prendergast, N.J., Smith, P.L. & Freisheim, J.H. (1988) Biochemistry 27, 4800-4804). To investigate this observation further, site-directed mutagenesis was used to construct deletion mutants of hDHFR missing 1 (del-1), 2 (del-2), 4 (del-4), and 6 (del-6) residues from loops in the carboxyl terminus of the enzyme. The del-1 mutant enzyme has a two-amino acid substitution in addition to the one-amino acid deletion. Deletion of only one amino acid resulted in a 35% decrease in the specific activity of the enzyme. The del-6 mutant enzyme was inactive. Surprisingly, the del-4 mutant enzyme retained a specific activity almost 33% that of the wild type. The specific activity of the del-2 mutant enzyme was slightly higher (38% wild-type activity) than that of the del-4 mutant. All three active deletion mutants were much less stable than the wild-type enzyme, and all three showed at least a 10-fold increase in Km values for both substrates. The del-1 and del-2 mutants exhibited a similar increase in KD values for both substrate and cofactor. The three active deletion mutants lost activity at concentrations of activating agents such as KCl, urea, and p-hydroxymercuribenzoate that continued to stimulate the wild-type enzyme. Antibody binding studies revealed conformational differences between the wild-type and mutant enzymes both in the absence and presence of bound folate. Thus, although the loops near the carboxyl terminus are far removed from the active site, small deletions of this region significantly affect DHFR function, indicating that the loop structure in mammalian DHFR plays an important functional role in its conformation and catalysis.     

Toward better refinement of comparative models: predicting loops in inexact environments

Achieving atomic-level accuracy in comparative protein models is limited by our ability to refine the initial, homolog-derived model closer to the native state. Despite considerable effort, progress in developing a generalized refinement method has been limited. In contrast, methods have been described that can accurately reconstruct loop conformations in native protein structures. We hypothesize that loop refinement in homology models is much more difficult than loop reconstruction in crystal structures, in part, because side-chain, backbone, and other structural inaccuracies surrounding the loop create a challenging sampling problem; the loop cannot be refined without simultaneously refining adjacent portions. In this work, we single out one sampling issue in an artificial but useful test set and examine how loop refinement accuracy is affected by errors in surrounding side-chains. In 80 high-resolution crystal structures, we first perturbed 6-12 residue loops away from the crystal conformation, and placed all protein side chains in non-native but low energy conformations. Even these relatively small perturbations in the surroundings made the loop prediction problem much more challenging. Using a previously published loop prediction method, median backbone (N-Calpha-C-O) RMSD's for groups of 6, 8, 10, and 12 residue loops are 0.3/0.6/0.4/0.6 A, respectively, on native structures and increase to 1.1/2.2/1.5/2.3 A on the perturbed cases. We then augmented our previous loop prediction method to simultaneously optimize the rotamer states of side chains surrounding the loop. Our results show that this augmented loop prediction method can recover the native state in many perturbed structures where the previous method failed; the median RMSD's for the 6, 8, 10, and 12 residue perturbed loops improve to 0.4/0.8/1.1/1.2 A. Finally, we highlight three comparative models from blind tests, in which our new method predicted loops closer to the native conformation than first modeled using the homolog template, a task generally understood to be difficult. Although many challenges remain in refining full comparative models to high accuracy, this work offers a methodical step toward that goal.     

The electrostatic potential of Escherichia coli dihydrofolate reductase

Escherichia coli dihydrofolate reductase (DHFR) carries a net charge of -10 electrons yet it binds ligands with net charges of -4 (NADPH) and -2 (folate or dihydrofolate). Evaluation and analysis of the electrostatic potential of the enzyme give insight as to how this is accomplished. The results show that the enzyme is covered by an overall negative potential (as expected) except for the ligand binding sites, which are located inside "pockets" of positive potential that enable the enzyme to bind the negatively charged ligands. The electrostatic potential can be related to the asymmetric distribution of charged residues in the enzyme. The asymmetric charge distribution, along with the dielectric boundary that occurs at the solvent-protein interface, is analogous to the situation occurring in superoxide dismutase. Thus DHFR is another case where the shape of the active site focuses electric fields out into solution. The positive electrostatic potential at the entrance of the ligand binding site in E. coli DHFR is shown to be a direct consequence of the presence of three positively charged residues at positions 32, 52, and 57--residues which have also been shown recently to contribute significantly to electronic polarization of the ligand folate. The latter has been postulated to be involved in the catalytic process. A similar structural motif of three positively charged amino acids that gives rise to a positive potential at the entrance to the active site is also found in DHFR from chicken liver, and is suggested to be a common feature in DHFRs from many species. It is noted that, although the net charges of DHFRs from different species vary from +3 to -10, the enzymes are able to bind the same negatively charged ligands, and perform the same catalytic function.