Dihydrofolate reductase: sequential resonance assignments using 2D and 3D NMR and secondary structure determination in solution

Three-dimensional (3D) heteronuclear NMR techniques have been used to make sequential 1H and 15N resonance assignments for most of the residues of Lactobacillus casei dihydrofolate reductase (DHFR), a monomeric protein of molecular mass 18,300 Da. A uniformly 15N-labeled sample of the protein was prepared and its complex with methotrexate (MTX) studied by 3D 15N/1H nuclear Overhauser-heteronuclear multiple quantum coherence (NOESY-HMQC), Hartmann-Hahn-heteronuclear multiple quantum coherence (HOHAHA-HMQC), and HMQC-NOESY-HMQC experiments. These experiments overcame most of the spectral overlap problems caused by chemical shift degeneracies in 2D spectra and allowed the 1H-1H through-space and through-bond connectivities to be identified unambiguously, leading to the resonance assignments. The novel HMQC-NOESY-HMQC experiment allows NOE cross peaks to be detected between NH protons even when their 1H chemical shifts are degenerate as long as the amide 15N chemical shifts are nondegenerate. The 3D experiments, in combination with conventional 2D NOESY, COSY, and HOHAHA experiments on unlabelled and selectively deuterated DHFR, provide backbone assignments for 146 of the 162 residues and side-chain assignments for 104 residues of the protein. Data from the NOE-based experiments and identification of the slowly exchanging amide protons provide detailed information about the secondary structure of the binary complex of the protein with methotrexate. Sequential NHi-NHi+1 NOEs define four regions with helical structure. Two of these regions, residues 44-49 and 79-89, correspond to within one amino acid to helices C and E in the crystal structure of the DHFR.methotrexate.NADPH complex [Bolin et al. (1982) J. Biol. Chem. 257, 13650-13662], while the NMR-determined helix formed by residues 26-35 is about one turn shorter at the N-terminus than helix B in the crystal structure, which spans residues 23-34. Similarly, the NMR-determined helical region comprising residues 102-110 is somewhat offset from the crystal structure's helix F, which encompasses residues 97-107. Regions of beta-sheet structure were characterized in the binary complex by strong alpha CHi-NHi+1 NOEs and by slowly exchanging amide protons. In addition, several long-range NOEs were identified linking together these stretches to form a beta-sheet. These elements align perfectly with corresponding elements in the crystal structure of the DHFR.methotrexate.NADPH complex, which contains an eight-stranded beta-sheet, indicating that the main body of the beta-sheet is preserved in the binary complex in solution.     

Structural studies of alpha-bungarotoxin. 1. Sequence-specific 1H NMR resonance assignments

We report the complete sequence-specific assignment of the backbone resonances and most of the side-chain resonances in the 1H NMR spectrum of alpha-bungarotoxin by two-dimensional NMR. Problems with resonance overlap were resolved with the assistance of the HRNOESY experiment described in an accompanying paper [Basus, V.J., & Scheek, R.M. (1988) Biochemistry (second paper of three in this issue)]. Significant differences exist between the solution structure described here and the crystal structure of alpha-bungarotoxin, on the basis of the proton to proton distances obtained by nuclear Overhauser enhancement spectroscopy (NOESY) and the corresponding distances from the X-ray crystal structure [Love, R.A., & Stroud, R.M. (1986) Protein Eng. 1, 37]. These differences include a larger beta-sheet in solution and a different orientation of the invariant tryptophan, Trp-28, making the solution structure more consistent with the crystal structure of the homologous neurotoxin alpha-cobratoxin. Four errors in the order of the amino acids in the primary sequence were indicated by the NMR data. These errors were confirmed by chemical means, as described in an accompanying paper [Kosen, P.A., Finer-Moore, J., McCarthy, M.P., & Basus, V.J. (1988) Biochemistry (third paper of three in this issue)].     

19F-n.m.r. studies of 3'',5''-difluoromethotrexate binding to Lactobacillus casei dihydrofolate reductase. Molecular motion and coenzyme-induced conformational changes.

19F-n.m.r. spectroscopy was used to study the binding of 3',5'-difluoromethotrexate to dihydrofolate reductase (tetrahydrofolate dehydrogenase) from Lactobacillus casei. The benzoyl ring of the bound difluoromethotrexate was found to 'flip' about its symmetry axis, and the rate (7.3 X 10(3) s-1 at 298 K) and activation parameters for this process were determined by lineshape analysis of the 19F-n.m.r. spectrum at a series of temperatures in the range 273-308 K. The contributions to the barrier for this process are discussed. Addition of NADP+ or NADPH to form the enzyme-difluoromethotrexate-coenzyme ternary complex led to an increase in the rate of benzoyl ring flipping by a factor of 2.6-2.8-fold, and to substantial changes in the 19F-n.m.r. chemical shifts. The possible nature of the coenzyme-induced conformational changes responsible for these effects is discussed.     

1H, 13C and 15N resonance assignments of the arsenate reductase from Synechocystis sp. strain PCC 6803

Arsenate reductases (ArsC) are a group of enzymes that play essential roles in biological arsenic detoxification pathways by catalyzing the intracellular reduction of arsenate to arsenite, which is subsequently extruded from the cells by specific transport systems. The ArsC protein from cyanobacterium Synechocystis sp. strain PCC 6803 (SynArsC) is related to the thioredoxin-dependent ArsC family, but uses the glutathione/glutaredoxin system for arsenate reduction. Therefore, it is classified to a novel thioredoxin/glutaredoxin hybrid arsenate reductase family. Herein we report the chemical shift assignments of ¹H, ¹³C and ¹⁵N atoms for the reduced form of SynArsC, which provides a starting point for further structural analysis and elucidation of its enzymatic mechanism.     

Nuclear magnetic resonance studies of the snake toxin echistatin. 1H resonance assignments and secondary structure.

The 1H-NMR spectrum of the snake toxin echistatin has been assigned using homonuclear two-dimensional methods. Consideration of the NOE patterns, coupling constants and putative hydrogen bonds enabled two regular features of secondary structure to be deduced: a beta-sheet/turn between residues 8 and 13 and a small anti-parallel beta-sheet and bulge linking residues 16-20 with residues 30-33. The recognition region of the protein containing the residues RGD lies in a loop joining the two strands of the beta-sheet. The beta-bulge and the loop containing the RGD sequence undergo pH-dependent conformational interconversion, modulated by the side chain of Asp29.     

Sequential resonance assignments in protein 1H nuclear magnetic resonance spectra. Basic pancreatic trypsin inhibitor

The assignment of the ¹H nuclear magnetic resonance spectrum of the basic pancreatic trypsin inhibitor with the use of two-dimensional ¹H nuclear magnetic resonance techniques at 500 MHz is described. The assignments are based entirely on the known amino acid sequence and the nuclear magnetic resonance data. Individual resonance assignments were obtained for all backbone and C^β protons, with the exception of those of Arg1, Pro2, Pro13 and the amide proton of Gly37. The side-chain resonance assignments are complete, with the exception of Pro2 and Pro13, the N^δ protons of Asn44 and the peripheral protons of the lysine residues and all but two of the arginine residues.     

1H NMR studies of echistatin in solution. Sequential resonance assignments and secondary structure.

Two-dimensional 1H-NMR methods have been used to obtain complete proton resonance assignments for the 49-residue protein echistatin from the viper Echis carinatus. The protein in solution contains only a small amount of regular secondary structure with four very short beta-strands. These beta-strands form two short segments of antiparallel beta-sheet, as evidenced by the observed cross-strand NOE. The first two strands are connected with a tight reverse turn, whereas the remaining two strands are linked together by an 11-residue loop forming a so-called hairpin. The tripeptide unit Arg-Gly-Asp, responsible for the binding of echistatin to the fibrinogen receptor glycoprotein GPIIb/IIIa, is located at the tip of this very hydrophilic loop.     

1H NMR studies of lambda cro repressor. 2. Sequential resonance assignments of the 1H NMR spectrum

The cro repressor protein from bacteriophage lambda has been studied in solution by two-dimensional nuclear magnetic resonance spectroscopy (2D NMR). Following the approach of Wüthrich and co-workers [Wüthrich, K., Wider, G., Wagner, G., & Braun, W. (1982) J. Mol. Biol. 155, 311-319], individual spin systems were identified by J-correlated spectroscopy (COSY) supplemented, where necessary, by relayed coherence transfer spectroscopy (RELAY). Nuclear Overhauser effect spectroscopy (NOESY) was used to obtain sequence-specific assignments. From the two-dimensional spectra, the peptide backbone resonances (NH and C alpha H) for 65 of the 66 amino acids were assigned, as well as most of the side chain resonances. The chemical shifts for the assigned protons are reported at 35 degrees C in 10 mM potassium phosphate, pH 6.8, and in 10 mM potassium phosphate, pH 4.6, 0.2 M KCl, and 0.1 mM EDTA. Small shifts were observed for some resonances upon addition of salt, but no major changes in the spectrum were seen, indicating that no global structural change occurs between these ionic strengths. NOE patterns characteristic of alpha-helices, beta-strands, and turns are seen in various regions of the primary sequence. From the location of these regions the secondary structure of cro in solution appears to be virtually identical with the crystal structure [Anderson, W. F., Ohlendorf, D. H., Takeda, Y., & Matthews, B. W. (1981) Nature (London) 290, 754-758]. Missing assignments include the Pro-59 resonances and the peripheral protons of the eight lysine, the three arginine, and three of the five isoleucine residues.     

Two-dimensional 1H NMR studies of histidine-containing protein from Escherichia coli. 1. Sequential resonance assignments

Two-dimensional NMR studies at 500 MHz have been performed on the histidine-containing protein (HPr) from Escherichia coli. HPr is one of the phosphocarrier proteins involved in the bacterial phosphoenolpyruvate:sugar phosphotransferase system that is responsible for the concomitant phosphorylation and translocation of a number of sugars. Sequential resonance assignments of HPr are complete. The conventional method of sequential assignments involving J-correlated spectroscopy (COSY) and nuclear Overhauser spectroscopy (NOESY) has been supplemented by optimized relayed coherence transfer spectroscopy (RELAY) to help overcome the spectral overlap that is inevitable in the spectra of proteins the size of HPr. RELAY experiments were performed in H2O to obtain NH-C beta H connectivities and in D2O to obtain C alpha H-C gamma H connectivities. The abundance of relayed coherence transfer peaks in the two experiments greatly aided in the assignment process of the complicated protein spectrum. The assignments lay the groundwork for the determination of the solution structure of HPr, as described in the accompanying paper [Klevit, R. E., & Waygood, E. B. (1986) Biochemistry (third paper of three in this issue)].     

Dihydrofolate reductase. The stereochemistry of inhibitor selectivity

X-ray structural results are reported for 10 triazine and pyrimidine inhibitors of dihydrofolate reductase, each one studied as a ternary complex with NADPH and chicken dihydrofolate reductase. Analysis of these data and comparison with structural results from the preceding paper (Matthews, D.A., Bolin, J.T., Burridge, J.M., Filman, D.J., Volz, K.W., Kaufman, B. T., Beddell, C.R., Champness, J.N., Stammers, D.K., and Kraut, J. (1985) J. Biol. Chem. 260, 381-391) in which we contrasted binding of the antibiotic trimethoprim (TMP) to chicken dihydrofolate reductase on the one hand with its binding to Escherichia coli dihydrofolate reductase on the other, permit identification of differences that are important in accounting for TMP's selectivity. The crystallographic evidence strongly suggests that loss of a potential hydrogen bond between the 4-amino group of TMP and the backbone carbonyl of Val-115 when TMP binds to chicken dihydrofolate reductase but not when it binds to the E. coli reductase is the major factor responsible for this drug's more potent inhibition of bacterial dihydrofolate reductase. A key finding of the current study which is important in understanding why TMP binds differently to chicken and E. coli dihydrofolate reductases is that residues on opposite sides of the active-site cleft in chicken dihydrofolate reductase are about 1.5-2.0 A further apart than are structurally equivalent residues in the E. coli enzyme.     

Determining stereo-specific 1H nuclear magnetic resonance assignments from distance geometry calculations

Stereo-specific 1H nuclear magnetic resonance assignments can be obtained following distance geometry structure calculations. The key to this method is to allow stereo-related atoms or methyls to float between pro-R and pro-S configurations, the final configuration being determined by the experimental constraints. Resonances from stereo-related pairs are given initial random assignments (either pro-R or pro-S) for identifying nuclear Overhauser effects (NOEs). A list of distance constraints using these assignments is compiled and a series of structures calculated where the chirality of non-C alpha chiral centers is not constrained; no pseudoatom corrections are required. Calculated structures are both locally and globally well-determined since the assignments rely upon the structure determination rather than the structure quality relying upon stereo-specific assignments. The method represents a global approach to determining stereo-specific assignments versus previously reported methods where only intraresidue NOEs and J-coupling information are used.     

Dihydrofolate reductase from bovine liver. Enzymatic and structural properties.

Dihydrofolate reductase from bovine liver has been purified 5000-fold employing conventional techniques and methotrexate/aminohexyl/Sepharose affinity chromatography. Electrophoresis of the isolated enzyme on polyacrylamide gels resulted in the separation of two enzymatically active protein components which were not interconvertible by treatment with dihydrofolate and/or the coenzyme. The two forms, present in a ratio of 20:1, were found by isoelectric focusing to have isoelectric points of 7.15 and 5.94. They had identical specific activities toward dihydrofolate (26.1-27.0 U/mg) and folate (1.3-2.2 U/mg), and had identical molecular weights (23500) and amino acid compositions. Due to the small quantity of the acidic form and the similarity of the two forms, the amino-terminal sequence (19 residues) was determined on a mixture of carboxymethylated reductase. The single sulfhydryl group of the enzyme can be modified by several sulfhydryl reagents in the native enzyme without loss of activity. Modification of the same residue occurs in the denaturated state and partially inhibits renaturation to the fully acitve enzyme. One disulfide bridge was detected by reduction and alkylation. The cleavage of this bond did not effect the enzymatic activity.     

1H-NMR stereospecific assignments by conformational data-base searches

A search procedure is described for making stereospecific assignments at prochiral centers in proteins on the basis of nuclear Overhauser enhancement and coupling constant data derived from nmr experiments. A data base comprising torsion angles, associated 1H-1H coupling constants and interproton distances is searched by a computer algorithm for sets of values that match the experimental data within specified error limits. Two different data bases are used. The first is a crystallographic data base derived from 34 well-refined crystal structures; the second is a systematic data base derived from conformations of a short peptide fragment with idealized geometry by systematically varying the phi, psi, and chi 1 torsion angles. Both approaches are tested for beta-methylene groups with model data obtained from 20 crystal structures. The results for the two methods are similar though not identical, so that a combination of the two methods appears to be useful. With an appropriate choice of error estimates, around 80% of the beta-methylene groups could be assigned in the test calculations. In addition, results with experimental nmr data indicate that a similar percentage of stereospecific assignments can be made in practical situations.     

Dihydrofolate reductase and aminopterin resistance in Pneumococcus

Summary Wild-type pneumococci derived from Avery's strain R36A are sensitive to extracellular concentrations of the folate antimetabolite aminopterin exceeding 1.0x10^-6 M. Three classes of resistant strains are phenotypically distinguishable: amiB-r, amiA-r and amiD-r strains are resistant to low (1.5x10^-6 M), intermediate (0.5–4.0×10^-5 M) and high (4.5x10^-4 M) aminopterin levels respectively. The amiA and amiB regions are weakly linked, but linkage has not been established between either of these loci and the amiD region.Consistent with the maximum resistance conferred by mutations in the amiA locus, dihydrofolate (FH₂) reductase in cell-free extracts (CFE) of amiA-r strains has a two- to six-fold greater affinity for the substrate than dose the enzyme in wild-type CFE (Table 1); FH₂ reductase from amiA-r strains may also have reduced affinity for aminopterin. Specific activity of the enzyme is not affected by mutation in the amiA locus (Table 1) and its affinity for the cofactor (NADPH) is probably unaffected by mutation in this locus (Table 4). Dihydrofolate reductase activity in amiA5 CFE is considerably more thermolabile than that in wild-type CFE (Table 2).The enzyme in CFE of the high resistance strain amiD1 has the same affinity for the substrate, cofactor and antimetabolite as FH₂ reductase in wild-type CFE (Figs. 1–4, 8 and 9; Table 4). However, specific activity of the enzyme in amiD1 CFE is 11-fold higher than that in wild-type CFE (Table 1) and it is much more heat stable (Table 2).Some properties of FH₂ reductase in CFE of the high resistance recombinant strain amiA5amiD1 are intermediate between those in CFE of wild-type and amiD1.Preliminary results suggest that CFE of wild-type and amiA5 contain a factor, which is neither dialyzable nor heat sensitive, that has an inhibitory effect upon activity and stability of FH₂ reductase in amiD1 CFE (Tables 2 and 3).