Evidence for two interconverting protein isomers in the methotrexate complex of dihydrofolate reductase from Escherichia coli

Distal Regions Regulate Dihydrofolate Reductase-Ligand Interactions

Protein motions play a fundamental role in enzyme catalysis and ligand binding. The relationship between protein motion and function has been extensively investigated in the model enzyme dihydrofolate reductase (DHFR). DHFR is an essential enzyme that catalyzes the reduction of dihydrofolate to tetrahydrofolate. Numerous experimental and computational methods have been used to probe the motions of DHFR through the catalytic cycle and to investigate the effect of distal mutations on DHFR motions and ligand binding. These experimental investigations have pushed forward the study of protein motions and their role in protein-ligand interactions. The introduction of mutations distal to the active site has been shown to have profound effects on ligand binding, hydride transfer rates and catalytic efficacy and these changes are captured by enzyme kinetics measurements. Distal mutations have been shown to exert their effects through a network of correlated amino acids and these effects have been investigated by NMR, protein dynamics, and analysis of coupled amino acids. The experimental methods and the findings that are reviewed here have broad implications for our understanding of enzyme mechanisms, ligand binding and for the future design and discovery of enzyme inhibitors.

Target identification by chromatographic co-elution: monitoring of drug-protein interactions without immobilization or chemical derivatization

Bioactive molecules typically mediate their biological effects through direct physical association with one or more cellular proteins. The detection of drug-target interactions is therefore essential for the characterization of compound mechanism of action and off-target effects, but generic label-free approaches for detecting binding events in biological mixtures have remained elusive. Here, we report a method termed target identification by chromatographic co-elution (TICC) for routinely monitoring the interaction of drugs with cellular proteins under nearly physiological conditions in vitro based on simple liquid chromatographic separations of cell-free lysates. Correlative proteomic analysis of drug-bound protein fractions by shotgun sequencing is then performed to identify candidate target(s). The method is highly reproducible, does not require immobilization or derivatization of drug or protein, and is applicable to diverse natural products and synthetic compounds. The capability of TICC to detect known drug-protein target physical interactions (K(d) range: micromolar to nanomolar) is demonstrated both qualitatively and quantitatively. We subsequently used TICC to uncover the sterol biosynthetic enzyme Erg6p as a novel putative anti-fungal target. Furthermore, TICC identified Asc1 and Dak1, a core 40 S ribosomal protein that represses gene expression, and dihydroxyacetone kinase involved in stress adaptation, respectively, as novel yeast targets of a dopamine receptor agonist.

The solution structure of Bacillus anthracis dihydrofolate reductase yields insight into the analysis of structure-activity relationships for novel inhibitors

There is a significant need for new therapeutics to treat infections caused by the biodefense agent Bacillus anthracis. In pursuit of drug discovery against this organism, we have developed novel propargyl-linked inhibitors that target the essential enzyme dihydrofolate reductase (DHFR) from B. anthracis. Previously, we reported an initial series of these inhibitors and a high-resolution crystal structure of the ternary complex of the enzyme bound to its cofactor and one of the most potent inhibitors, UCP120B [Beierlein, J., Frey, K., Bolstad, D., Pelphrey, P., Joska, T., Smith, A., Priestley, N., Wright, D., and Anderson, A. (2008) J. Med. Chem. 51, 7532-7540]. Herein, we describe a three-dimensional solution structure of the ternary complex as determined by NMR. A comparison of this solution structure to the crystal structure reveals a general conservation of the DHFR fold and cofactor interactions as well as differences in the location of an active site helix and specific ligand interactions. In addition to data for the fully assigned ternary complex, data for the binary (enzyme-cofactor) complex were collected, providing chemical shift comparisons and revealing perturbations in residues that accommodate ligand binding. Dynamics of the protein, measured using (15)N T(1) and T(2) relaxation times and {(1)H}-(15)N heteronuclear NOEs, reveal residue flexibility at the active site that explains enzyme inhibition and structure-activity relationships for two different series of these propargyl-linked inhibitors. The information obtained from the solution structure regarding active site flexibility will be especially valuable in the design of inhibitors with increased potency.

Carbon-deuterium bonds as probes of dihydrofolate reductase

Much effort has been directed toward understanding the contributions of electrostatics and dynamics to protein function and especially to enzyme catalysis. Unfortunately, these studies have been limited by the absence of direct experimental probes. We have been developing the use of carbon-deuterium bonds as probes of proteins and now report the application of the technique to the enzyme dihydrofolate reductase, which catalyzes a hydride transfer and has served as a paradigm for biological catalysis. We observe that the stretching absorption frequency of (methyl- d 3) methionine carbon-deuterium bonds shows an approximately linear dependence on solvent dielectric. Solvent and computational studies support the empirical interpretation of the stretching frequency in terms of local polarity. To begin to explore the use of this technique to study enzyme function and mechanism, we report a preliminary analysis of (methyl- d 3) methionine residues within dihydrofolate reductase. Specifically, we characterize the IR absorptions at Met16 and Met20, within the catalytically important Met20 loop, and Met42, which is located within the hydrophobic core of the enzyme. The results confirm the sensitivity of the carbon-deuterium bonds to their local protein environment, demonstrate that dihydrofolate reductase is electrostatically and dynamically heterogeneous, and lay the foundation for the direct characterization protein electrostatics and dynamics and, potentially, their contribution to catalysis.

Substrate and inhibitor specificity of Mycobacterium avium dihydrofolate reductase

Dihydrofolate reductase (EC 1.5.1.3) is a key enzyme in the folate biosynthetic pathway. Information regarding key residues in the dihydrofolate-binding site of Mycobacterium avium dihydrofolate reductase is lacking. On the basis of previous information, Asp31 and Leu32 were selected as residues that are potentially important in interactions with dihydrofolate and antifolates (e.g. trimethoprim), respectively. Asp31 and Leu32 were modified by site-directed mutagenesis, giving the mutants D31A, D31E, D31Q, D31N and D31L, and L32A, L32F and L32D. Mutated proteins were expressed in Escherichia coli BL21(DE3)pLysS and purified using His-Bind resin; functionality was assessed in comparison with the recombinant wild type by a standard enzyme assay, and growth complementation and kinetic parameters were evaluated. All Asp31 substitutions affected enzyme function; D31E, D31Q and D31N reduced activity by 80-90%, and D31A and D31L by > 90%. All D31 mutants had modified kinetics, ranging from three-fold (D31N) to 283-fold (D31L) increases in K(m) for dihydrofolate, and 12-fold (D31N) to 223 077-fold (D31L) decreases in k(cat)/K(m). Of the Leu32 substitutions, only L32D caused reduced enzyme activity (67%) and kinetic differences from the wild type (seven-fold increase in K(m); 21-fold decrease in k(cat)/K(m)). Only minor variations in the K(m) for NADPH were observed for all substitutions. Whereas the L32F mutant retained similar trimethoprim affinity as the wild type, the L32A mutation resulted in a 12-fold decrease in affinity and the L32D mutation resulted in a seven-fold increase in affinity for trimethoprim. These findings support the hypotheses that Asp31 plays a functional role in binding of the substrate and Leu32 plays a functional role in binding of trimethoprim.

Loop anchor modification causes the population of an alternative native state in an SH3-like domain

Many stably folded proteins are proposed to contain long, unstructured loops. A series of hybrid proteins (EbE1-4) containing the folded scaffold of photosystem I accessory protein E (PsaE), an SH3-like protein, and the 40-residue heme-binding loop of cytochrome b(5) was created to inspect the dependence of thermodynamic and kinetic parameters on the residues at the interface of folded and flexible regions. Compared to the simplest hybrid (EbE1), the chimeras differed by Gly insertions (EbE2, EbE3) or an asymmetric four-residue restructuring of loop termini (EbE4). NMR spectroscopy indicated that the chimeras retained the PsaE topology; native and unfolded state solubilities, however, were affected to varying degrees. Thermal and chemical denaturation experiments revealed that the EbE2 and EbE1 constructs resulted in a modest destabilization of the PsaE core, whereas apparent stability was increased by >5 kJ/mol in EbE4. EbE3 aggregated at microM concentrations and was not studied in detail. EbE4 populated two native states (N1 and N2), which differed by hydrophobic core packing and C-terminal interactions. At room temperature, the population ratio ( approximately 3-4:1) favored the state whose spectroscopic properties most resembled those of PsaE (N1). EbE4 also demonstrated altered folding kinetics, displaying multiple slow phases related to the population of intermediates and possibly N2. It was concluded that loop anchors can affect protein properties, including stability, via short-range effects on local structure and long-range communication with the packed hydrophobic core. Modification of the attachment points appears to be a possible stepping stone in the transition from one three-dimensional structure to another.

The relationship between chain connectivity and domain stability in the equilibrium and kinetic folding mechanisms of dihydrofolate reductase from E.coli

The role of domains in defining the equilibrium and kinetic folding properties of dihydrofolate reductase (DHFR) from Escherichia coli was probed by examining the thermodynamic and kinetic properties of a set of variants in which the chain connectivity in the discontinuous loop domain (DLD) and the adenosine-binding domain (ABD) was altered by permutation. To test the concept that chain cleavage can selectively destabilize the domain in which the N- and C-termini are resident, permutations were introduced at one position within the ABD, one within the DLD and one at a boundary between the domains. The results demonstrated that a continuous ABD is required for a stable thermal intermediate and a continuous DLD is required for a stable urea intermediate. The permutation at the domain interface had both a thermal and urea intermediate. Strikingly, the observable kinetic folding responses of all three permuted proteins were very similar to the wild-type protein. These results demonstrate a crucial role for stable domains in defining the energy surface for the equilibrium folding reaction of DHFR. If domain connectivity affects the kinetic mechanism, the effects must occur in the sub-millisecond time range.

Direct NMR observation of a substrate protein bound to the chaperonin GroEL

The reaction cycle and the major structural states of the molecular chaperone GroEL and its cochaperone, GroES, are well characterized. In contrast, very little is known about the nonnative states of the substrate polypeptide acted on by the chaperonin machinery. In this study, we investigated the substrate protein human dihydrofolate reductase (hDHFR) while bound to GroEL or to a single-ring analog, SR1, by NMR spectroscopy in solution under conditions where hDHFR was efficiently recovered as a folded, enzymatically active protein from the stable complexes upon addition of ATP and GroES. By using the NMR techniques of transverse relaxation-optimized spectroscopy (TROSY), cross-correlated relaxation-induced polarization transfer (CRIPT), and cross-correlated relaxation-enhanced polarization transfer (CRINEPT), bound hDHFR could be observed directly. Measurements of the buildup of hDHFR NMR signals by different magnetization transfer mechanisms were used to characterize the dynamic properties of the NMR-observable parts of the bound substrate. The NMR data suggest that the bound state includes random coil conformations devoid of stable native-like tertiary contacts and that the bound hDHFR might best be described as a dynamic ensemble of randomly structured conformers.

Structure, dynamics, and catalytic function of dihydrofolate reductase

Molecular motions are widely regarded as contributing factors in many aspects of protein function. The enzyme dihydrofolate reductase (DHFR), and particularly that from Escherichia coli, has become an important system for investigating the linkage between protein dynamics and catalytic function, both because of the location and timescales of the motions observed and because of the availability of a large amount of structural and mechanistic data that provides a detailed context within which the motions can be interpreted. Changes in protein dynamics in response to ligand binding, conformational change, and mutagenesis have been probed using numerous experimental and theoretical approaches, including X-ray crystallography, fluorescence, nuclear magnetic resonance (NMR), molecular dynamics simulations, and hybrid quantum/classical dynamics methods. These studies provide a detailed map of changes in conformation and dynamics throughout the catalytic cycle of DHFR and give new insights into the role of protein motions in the catalytic activity of this enzyme.

The coordination of the isomerization of a conserved non-prolyl cis peptide bond with the rate-limiting steps in the folding of dihydrofolate reductase

The propensity for peptide bonds to adopt the trans configuration in native and unfolded proteins, and the relatively slow rates of cis-trans isomerization reactions, imply that the formation of cis peptide bonds in native conformations are likely to limit folding reactions. The role of the conserved cis Gly95-Gly96 peptide bond in dihydrofolate reductase (DHFR) from Escherichia coli was examined by replacing Gly95 with alanine. The introduction of a beta carbon at position 95 is expected to increase the propensity for the trans isomer and perturb the isomerization reaction required to reach the native conformation. Although G95A DHFR is 1.30 kcal mol(-1) less stable than the wild-type protein, it adopts a well-folded structure that can be chemically denatured in a cooperative fashion. The mutant protein also retains the complex refolding kinetic pattern attributed to a parallel-channel mechanism in wild-type DHFR. The spectroscopic response upon refolding monitored by Trp fluorescence and the absence of a Trp/Trp exciton coupling apparent in the far-UV CD spectrum of the wild-type protein, however, indicated that the tertiary structure of the folded state for G95A DHFR is altered. The addition of methotrexate (MTX), a tight-binding inhibitor, to folded G95A DHFR restored the exciton coupling and the fluorescence properties through five slow kinetic events whose relaxation times are independent of the ligand and the denaturant concentrations. The results were interpreted to mean that MTX-binding drives the formation of the cis isomer of the peptide bond between Ala95 and Gly96 in five compact and stable but not wild-type-like conformations that contain the trans isomer. Folding studies in the presence of MTX for both wild-type and G95A DHFR support the notion that the cis peptide bond between Gly95 and Gly96 in the wild-type protein forms during four parallel rate-limiting steps, which are primarily controlled by folding reactions, and lead directly to a set of native, or native-like, conformers. The isomerization of the cis peptide bond is not a source of the parallel channels that characterize the complex folding mechanism for DHFR.

BetaCore, a designed water soluble four-stranded antiparallel beta-sheet protein

BetaCore is a designed approximately 50-residue protein in which two BPTI-derived core modules, CM I and CM II, are connected by a 22-atom cross-link. At low temperature and pH 3, homo- and heteronuclear NMR data report a dominant folded ('f') conformation with well-dispersed chemical shifts, i, i+1 periodicity, numerous long-range NOEs, and slowed amide hydrogen isotope exchange patterns that is a four-stranded antiparallel beta-sheet with nonsymmetrical and specific association of CM I and CM II. BetaCore 'f' conformations undergo reversible, global, moderately cooperative, non-two-state thermal transitions to an equilibrium ensemble of unfolded 'u' conformations. There is a significant energy barrier between 'f' and 'u' conformations. This is the first designed four-stranded antiparallel beta-sheet that folds in water.

Testing the role of chain connectivity on the stability and structure of dihydrofolate reductase from E. coli: fragment complementation and circular permutation reveal stable, alternatively folded forms

The effects of chain cleavage and circular permutation on the structure, stability, and activity of dihydrofolate reductase (DHFR) from Escherichia coli were investigated by various spectroscopic and biochemical methods. Cleavage of the backbone after position 86 resulted in two fragments, (1--86) and (87--159) each of which are poorly structured and enzymatically inactive. When combined in a 1 : 1 molar ratio, however, the fragments formed a high-affinity (K(a) = 2.6 x 10(7) M(-1)) complex that displays a weakly cooperative urea-induced unfolding transition at micromolar concentrations. The retention of about 15% of the enzymatic activity of full-length DHFR is surprising, considering that the secondary structure in the complex is substantially reduced from its wild-type counterpart. In contrast, a circularly permuted form with its N-terminus at position 86 has similar overall stability to full-length DHFR, about 50% of its activity, substantial secondary structure, altered side-chain packing in the adenosine binding domain, and unfolds via an equilibrium intermediate not observed in the wild-type protein. After addition of ligand or the tight-binding inhibitor methotrexate, both the fragment complex and the circular permutant adopt more native-like secondary and tertiary structures. These results show that changes in the backbone connectivity can produce alternatively folded forms and highlight the importance of protein-ligand interactions in stabilizing the active site architecture of DHFR.

Effect of ligand binding on the flexibility of dihydrofolate reductase as revealed by compressibility

The partial specific volume, v, and adiabatic compressibility, beta(s), of Escherichia coli dihydrofolate reductase were measured at 30 degrees C in the presence of various ligands (folate, dihydrofolate, tetrahydrofolate, NADPH, NADP, methotrexate, and KCl). Binding of these ligands (binary and ternary complexes) brought about large changes of v (0.734-0.754 cm(3) g(-1)) and beta(s) (6. 6x10(-6)-9.8x10(-6) bar(-1)), keeping a linear relationship between the two parameters. The values of v and beta(s) increased with an increase in internal cavity, V(cav), and a decrease in accessible surface area, ASA, which were calculated from the X-ray crystal structures of the complexes. A large variation of V(cav) relative to ASA by ligand binding suggested that the cavity is a dominant factor and the effect of hydration might be small for the ligand-induced changes of v and beta(s). The beta(s) values of the binary and ternary complexes suggested a characteristic conformational flexibility of the kinetic intermediates in the enzyme reaction coordinate. Comparison of beta(s) with the cavity distribution in the crystal structures revealed that the flexibility of the intermediates was mainly determined by the total cavity volume with minor contributions of the number, position, and size of cavities. These results demonstrate that the compressibility is a useful measure of the conformational flexibility of the intermediates in the enzyme reaction and that the combined study of compressibility and X-ray crystallography gives new insight into the protein dynamics through the behavior of the cavities.

Cyclophilin-promoted folding of mouse dihydrofolate reductase does not include the slow conversion of the late-folding intermediate to the active enzyme

Cyclophilins accelerate slow protein folding reactions in vitro by catalyzing the cis/trans isomerization of peptidyl-prolyl bonds. Cyclophilins were reported to be involved in a variety of cellular functions, including the promotion of protein folding by use of the substrate mouse dihydrofolate reductase (DHFR). The interaction of cyclophilin with DHFR has only been studied under limited conditions so far, not taking into account that native DHFR exists in equilibrium with a non-native late-folding intermediate. Here we report a systematic analysis of catalysis of DHFR folding by cyclophilins. The specific ligand methotrexate traps DHFR in its native state, permitting a specific analysis of the action of cyclophilin on both denatured DHFR with non-native prolyl bonds and denatured DHFR with all-native prolyl bonds. Cyclophilins from yeast and Neurospora crassa as well as the related prolyl isomerase b from Escherichia coli promote the folding of different forms of DHFR to the enzymatically active form, demonstrating the generality of cyclophilin-catalyzed folding of DHFR. The slow equilibrium between the late-folding intermediate and native DHFR suggests that prolyl isomerization may be required for this final phase of conversion to native DHFR. However, by reversible trapping of the intermediate, we analyze the slow interconversion between native and late-folding conformations in the backward and forward reactions and show a complete independence of cyclophilin. We conclude that cyclophilin catalyzes folding of DHFR, but surprisingly not in the last slow folding step.

High expression and steady-state kinetic characterization of methionine site-directed mutants of Escherichia coli methionyl- and selenomethionyl-dihydrofolate reductase

A high expression system that produces Escherichia coli dihydrofolate reductase (DHFR) at 30% total cellular protein was constructed. This expression vector, named pCOCK, allowed for the purification of nearly 100 mg of homogeneous DHFR from a 11 bacterial culture. A simple, single Q-Sepharose anion exchange column purification was developed on an FPLC instrument. Methionine site-directed mutants were constructed in DHFR to assess the role of Met within the enzymes. These mutants consisted of a Met16leucine (Leu), Met20Leu, Met42Leu, Met92Leu, Met16,20Leu and Met16,20,42Leu. Steady-state kinetic studies showed that the Met16Leu, Met42Leu and Met92Leu mutants possessed essentially the same kcat, Km(DHF) and Km(NADPH) as that of wild-type (wt) DHFR (13.7 s-1, 0.97 microM and 2.52 microM, respectively). Mutants which contained a Leu at position 20 possessed substantially elevated specific activity and kcat values. The specific activity and kcat of wt, Met20Leu, Met16,20Leu and Met16,20,42Leu were 45.9, 92.7, 90.2 and 172 mumol/min/mg and 13.7, 24.6, 25.2 and 52.7 s-1, respectively. Upon substitution of Met by selenomethionine (SeMet) in the aforementioned mutants, further information as to the effect of SeMet incorporation into proteins was ascertained. Steady-state kinetic parameters of the SeMet substituted Met16Leu, Met20Leu, Met42Leu and Met92Leu mutants were nearly identical to those of their Met containing counterparts. These data indicate that Met apparently has a limited role in the protein structure and function of DHFR and that SeMet incorporation has no effect on the steady-state kinetic constants of DHFR.

Dynamics of the dihydrofolate reductase-folate complex: catalytic sites and regions known to undergo conformational change exhibit diverse dynamical features

Backbone and tryptophan side-chain dynamics of uniformly 15N-labeled Escherichia coli dihydrofolate reductase were determined for the binary folate complex. The 15N T1 and T2 relaxation times and [1H]-15N heteronuclear NOEs were measured for 118 protonated backbone nitrogen atoms. The generalized order parameter (S2), the effective correlation time for internal motions (tau e), and the contribution to spin-spin relaxation through 15N exchange broadening (Rex) were determined for each residue by model-free analysis. Back-calculation of the relaxation rates for each resonance showed that the calculated dynamical parameters accurately predict the experimental data. Diverse dynamical features were evident in the DHFR backbone. Six sites exhibited order parameters significantly below the weighted mean S2 value (for the complex) of 0.81 +/- 0.002: residues G67 and D69 of the adenosine binding domain, and "hinge" residues K38 and V88, exhibited low S2 (0.29 < or = S2 < or = 0.6) and high tau e values (700 ps < or = tau e < or = 2 ns), as did sites within the beta A-alpha B loop and the beta F-beta G loop. Thus, large amplitude backbone motions, on the picosecond and nanosecond time scales, occurred at regions implicated in transition-state stabilization and in ligand-dependent conformational change. Significant Rex values (> or = 1 s-1) were determined for 45% of assigned resonances, many of which arise from residues surrounding the folate binding site. The mean S2 values of the occupied folate binding site and the unoccupied NADPH binding site were similar, indicating the backbone of the latter is at least as conformationally restricted as that of the occupied folate site. We conclude that the observed time-dependent structural fluctuations of the binary complex are in fact associated with catalytic properties of the molecule.

Detection of a stable intermediate in the thermal unfolding of a cysteine-free form of dihydrofolate reductase from Escherichia coli

The reversible temperature-induced unfolding of a cysteine-free mutant (C85S/C152E, des-Cys) of dihydrofolate reductase from Escherichia coli has been studied by absorbance and by both far- and near-ultraviolet circular dichroism spectroscopies. The non-coincidence of all three transition curves demonstrated the existence of a highly populated partially-folded form near 39 degrees C at pH 7.8. This intermediate retains substantial secondary structure and partially excludes one or more of the five tryptophans from solvent; however, the intermediate has lost specific tertiary packing around its aromatic residues. Increases in enthalpy, entropy, and heat capacity are observed for both the native/intermediate and intermediate/unfolded transitions; the majority of the changes in these parameters occurs in the first transition. These results suggest that the thermal unfolding reaction of des-Cys dihydrofolate reductase involves a stable intermediate whose properties resemble those of a molten globule.

Early intermediates in the folding of dihydrofolate reductase from Escherichia coli detected by hydrogen exchange and NMR

The kinetic folding mechanism for Escherichia coli dihydrofolate reductase postulates two distinct types of transient intermediates. The first forms within 5 ms and has substantial secondary structure but little stability. The second is a set of four species that appear over the course of several hundred milliseconds and have secondary structure, specific tertiary structure, and significant stability (Jennings PA, Finn BE, Jones BE, Matthews CR, 1993, Biochemistry 32:3783-3789). Pulse labeling hydrogen exchange experiments were performed to determine the specific amide hydrogens in alpha-helices and beta-strands that become protected from exchange through the formation of stable hydrogen bonds during this time period. A significant degree of protection was observed for two subsets of the amide hydrogens within the dead time of this experiment (6 ms). The side chains of one subset form a continuous nonpolar strip linking six of the eight strands in the beta-sheet. The other subset corresponds to a nonpolar cluster on the opposite face of the sheet and links three of the strands and two alpha-helices. Taken together, these data demonstrate that the complex strand topology of this eight-stranded sheet can be formed correctly within 6 ms. Measurement of the protection factors at three different folding times (13 ms, 141 ms, and 500 ms) indicates that, of the 13 amide hydrogens displaying significant protection within 6 ms, 8 exhibit an increase in their protection factors from approximately 5 to approximately 50 over this time range; the remaining five exhibit protection factors > 100 at 13 ms. Only approximately half of the population of molecules form this set of stable hydrogen bonds. Thirteen additional hydrogens in the beta-sheet become protected from exchange as the set of native conformers appear, suggesting that the stabilization of this network reflects the global cooperativity of the folding reaction.

1H, 15N and 13C resonance assignments, secondary structure, and the conformation of substrate in the binary folate complex of Escherichia coli dihydrofolate reductase

By using fully 15N- and 15N/13C-labeled Escherichia coli dihydrofolate reductase, the sequence-specific 1H and 15N NMR assignments were achieved for 95% of the backbone resonances and for 90% of the 13C alpha resonances in the binary folate complex. These assignments were made through a variety of three-dimensional proton-detected 15N and 13C experiments. A smaller but significant subset of side-chain 1H and 13C assignments were also determined. In this complex, only one 15N or 13C resonance was detected per 15N or 13C protein nucleus, which indicated a single conformation. Proton-detected 13C experiments were also performed with unlabeled DHFR, complexed with 13C-7/13C-9 folate to probe for multiple conformations of the substrate in its binary complex. As was found for the protein resonances, only a single bound resonance corresponding to a productive conformation could be detected for C-7. These results are consistent with an earlier report based on 1H NMR data [Falzone, C.J. et al. (1990) Biochemistry, 29, 9667-9677] and suggest that the E. coli enzyme is not involved in any catalytically unproductive binding modes in the binary complex. This feature of the E. coli enzyme seems to be unique among the bacterial forms of DHFR that have been studied to date.

Dynamics of a flexible loop in dihydrofolate reductase from Escherichia coli and its implication for catalysis

Apo-dihydrofolate reductase from Escherichia coli samples two distinct environments slowly on the NMR time scale at room temperature. Several assigned resonances belong to residues in, or proximal to, a loop (loop I) which is comprised of residues 9-24. This exchange process was altered (either removed or made fast on the NMR time scale) by deleting three hairpin turn forming residues from the loop and filling the gap with a single glycine [Li, L., Falzone, C. J., Wright, P. E., & Benkovic, S. J. (1992) Biochemistry 31, 7826-7833]. An approximate value of 35 s-1 for the exchange rate associated with loop I in apo-DHFR was obtained in two-dimensional nuclear Overhauser spectra by analyzing the time dependence of the cross-peak volume for N epsilon H of Trp-22, a residue which is located in this loop and which has resolved cross-peaks. Owing to the critical role that this loop plays in catalysis, the correspondence between this rate of conformational exchange and off-rates for tetrahydrofolate and the reduced nicotinamide cofactor from product and substrate complexes suggests that loop movement may be a limiting factor in substrate turnover.

Collapse of parallel folding channels in dihydrofolate reductase from Escherichia coli by site-directed mutagenesis

The rate-limiting steps in the folding of dihydrofolate reductase from Escherichia coli have been shown to involve the conversion of a set of four intermediates to a corresponding set of native conformers via four parallel channels [Jennings et al. (1993) Biochemistry 32, 3783-3789]. Fluorescence and absorbance studies of the unfolding and refolding of the C85S/C152E double mutant at various final urea concentrations reveal two slow folding reactions, two fewer than observed in the wild-type protein. Refolding in the presence of substoichiometric levels of the inhibitor methotrexate shows that the two remaining slow reactions correspond to two parallel channels which lead to a pair of native conformers capable of binding the inhibitor. A combination of stopped-flow circular dichroism and cofactor binding studies confirms that the four parallel channels observed in the wild-type protein have collapsed into two channels in the mutant. Kinetic and equilibrium studies of the single cysteine mutants suggest that replacements of Cysteine-85 which perturb the hydrophobic core containing this side chain are responsible for the simplification of the kinetic mechanism. These results demonstrate that at least two of the parallel folding channels in dihydrofolate reductase arise when tertiary structure develops and are not dependent upon cis/trans isomerization at prolyl peptide bonds.

How do mutations at phenylalanine-153 and isoleucine-155 partially suppress the effects of the aspartate-27-->serine mutation in Escherichia coli dihydrofolate reductase?

Several second-site suppressors of the D27S lesion in Escherichia coli dihydrofolate reductase (DHFR) have been identified. The activity of the primary mutant, D27S DHRF, was found to be greatly decreased at pH 7.0, consistent with aspartic acid-27 being critically involved in proton donation during catalysis. Partial suppressors of the D27S mutation have been selected by their ability to confer an increased resistance to trimethoprim upon host E. coli; the suppressors have been identified as F153S or I155N substitutions. D27S+F153S and D27S+I155N DHFRs display 2-3-fold increases in kcat over D27S DHFR values, but only the F153S mutation decreases the Km for dihydrofolate by a factor of 2. Neither double mutant approaches wild-type DHFR activity. Unexpectedly, Phe153 and Ile155 occur on the surface of the protein and are approximately 8 and 14 A distant from the active site. Ile155 is a member of a beta-bulge. A previously identified suppressing mutation, F137S, occurs nearby and is also a member of the same beta-bulge [Howell et al. (1990) Biochemistry 29, 8561-8569]. Clustering of these three second-site mutations indicates this area of the structure may be important in protein function. Conformational changes due to the presence of these suppressing mutations are likely as the F153S and I155N mutations do not affect hydride-transfer rates upon introduction in wild-type DHFR and alterations in circular dichroism spectra are associated with the double-mutant DHFRs.

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.

Functional role of a mobile loop of Escherichia coli dihydrofolate reductase in transition-state stabilization

The function of a highly mobile loop in Escherichia coli dihydrofolate reductase was studied by constructing a mutant (DL1) using cassette mutagenesis that had four residues deleted in the middle section of the loop (Met16-Ala19) and a glycine inserted to seal the gap. This part of the loop involves residues 16-20 and is disordered in the X-ray crystal structures of the apoprotein and the NADP+ binary complex but forms a hairpin turn that folds over the nicotinamide moiety of NADP+ and the pteridine moiety of folate in the ternary complex [Bystroff, C., & Kraut, J. (1991) Biochemistry 30, 2227-2239]. The steady-state and pre-steady-state kinetics and two-dimensional 1H NMR spectra were analyzed and compared to the wild-type protein. The kinetics on the DL1 mutant enzyme show that the KM value for NADPH (5.3 microM), the KM for dihydrofolate (2 microM), the rate constant for the release of the product tetrahydrofolate (10.3 s-1), and the intrinsic pKa value (6.2) are similar to those exhibited by the wild-type enzyme. However, the hydride-transfer rate declines markedly from the wild-type value of 950 s-1 to 1.7 s-1 for the DL1 mutant and when taken with data for substrate binding indicates that the loop contributes to substrate flux by a factor of 3.5 x 10(4). Thus, the mobility of loop I may provide a mechanism of recruiting hydrophobic residues which can properly align the nicotinamide and pteridine rings for the hydride-transfer process (a form of transition-state stabilization).(ABSTRACT TRUNCATED AT 250 WORDS)

15N NMR studies of the conformation of E. coli dihydrofolate reductase in complex with folate or methotrexate

We have employed 15N NMR to characterize the conformations of Escherichia coli dihydrofolate reductase (ECDHFR) in complex with [5-15N]folate or [5-15N]methotrexate (MTX). Two 15N resonances were observed for DHFR/MTX binary complex. The relative population of these two conformations is pH dependent. Addition of NADP+ or NADPH results in the disappearance of the low field resonance. In contrast, only one conformation was observed for both the DHFR/folate and DHFR/folate/NADP+ complexes. However, the 15N chemical shift of [5-15N]folate in the binary DHFR/folate complex is 7.28 ppm upfield from that of the ternary complex, suggesting the possible loss of a hydrogen bonding to N5 of folate in the ternary complex.

Methotrexate binds in a non-productive orientation to human dihydrofolate reductase in solution, based on NMR spectroscopy

Dihydrofolate reductase (DHFR) is an intracellular target enzyme for folate antagonist drugs, including methotrexate. In order to compare the binding of methotrexate to human DHFR in solution with that observed in the crystalline state, NMR spectroscopy has been used to determine the conformation of the drug bound to human DHFR in solution. In agreement with what has been observed in the crystalline state, NOE's identified protein and methotrexate protons indicate that methotrexate binds in a non-productive orientation. In contrast to what has been reported for E. coli DHFR in solution, only one bound conformation of methotrexate is observed.