NMR studies of multiple conformations in complexes of Lactobacillus casei dihydrofolate reductase with analogues of pyrimethamine

1H and 19F NMR signals from bound ligands have been assigned in one- and two-dimensional NMR spectra of complexes of Lactobacillus casei dihydrofolate reductase with various pyrimethamine analogues (including pyrimethamine [1, 2,4-diamino-5-(4'-chlorophenyl)-6-ethylpyrimidine], fluoropyrimethamine [2, 2,4-diamino-5-(4'-fluorophenyl)-6-ethylpyrimidine], fluoronitropyrimethamine [3, 2,4-diamino-5-(4'-fluoro-3'-nitrophenyl) -6-ethylpyrimidine], and methylbenzoprim [4, 2,4-diamino-5-[4'- (methylbenzylamino)-3'-nitrophenyl]-6-ethylpyrimidine]). The signals were identified mainly by correlating signals from bound and free ligands by using 2D exchange experiments. Analogues (such as 1 and 2) with symmetrically substituted phenyl rings give rise to 1H signals from four nonequivalent aromatic protons, clearly indicating the presence of hindered rotation about the pyrimidine-phenyl bond. Analogues containing asymmetrically substituted aromatic rings (such as 3 and 4) exist as mixtures of two rotational isomers (an enantiomeric pair) because of this hindered rotation and the NMR spectra revealed that both isomers (forms A and B) bind to the enzyme with comparable, though unequal, binding energies. In this case two complete sets of bound proton signals were observed. The phenyl ring protons in each of the two forms experience essentially the same protein environment (same shielding) as that experienced by the corresponding protons in bound pyrimethamine: this confirms that forms A and B correspond to two rotational isomers resulting from approximately 180 degrees rotation about the pyrimidine-phenyl bond, with the 2,4-diaminopyrimidine ring being bound similarly in both forms. The relative orientations of the two forms have been determined from NOE through-space connections between protons on the ligand and protein.(ABSTRACT TRUNCATED AT 250 WORDS)     

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.     

The kinetic mechanism of wild-type and mutant mouse dihydrofolate reductases

A kinetic mechanism is presented for mouse dihydrofolate reductase that predicts all the steady-state parameters and full time-course kinetics. This mechanism was derived from association and dissociation rate constants and pre-steady-state transients by using stopped-flow fluorescence and absorbance measurements. The major features of this kinetic mechanism are as follows: (1) the two native enzyme conformers, E1 and E2, bind ligands with varying affinities although only one conformer, E1, can support catalysis in the forward direction, (2) tetrahydrofolate dissociation is the rate-limiting step under steady-state turnover at low pH, and (3) the pH-independent rate of hydride transfer from NADPH to dihydrofolate is fast (khyd = 9000 s-1) and favorable (Keq = 100). The overall mechanism is similar in form to the Escherichia coli kinetic scheme (Fierke et al., 1987), although several differences are observed: (1) substrates and products predominantly bind the same form of the E. coli enzyme, and (2) the hydride transfer rate from NADPH to either folate or dihydrofolate is considerably faster for the mouse enzyme. The role of Glu-30 (Asp-27 in E. coli) in mouse DHFR has also been examined by using site-directed mutagenesis as a potential source of these differences. While aspartic acid is strictly conserved in all bacterial DHFRs, glutamic acid is conserved in all known eucaryotes. The two major effects of substituting Asp for Glu-30 in the mouse enzyme are (1) a decreased rate of folate reduction and (2) an increased rate of hydride transfer from NADPH to dihydrofolate.(ABSTRACT TRUNCATED AT 250 WORDS)     

NMR studies of ligand carboxylate group interactions with arginine residues in complexes of Lactobacillus casei dihydrofolate reductase with substrates and substrate analogues

In a series of complexes of Lactobacillus casei dihydrofolate reductase (DHFR) formed with substrates and substrate analogues, the (1)H/(15)N NMR chemical shifts for the guanidino group of the conserved Arg 57 residue were found to be sensitive to the mode of binding of their H(eta) protons to the charged oxygen atoms in ligand carboxylate groups. In all cases, Arg 57 showed four nonequivalent H(eta) signals indicating hindered rotation about the N(epsilon)-C(zeta) and C(zeta)-N(eta) bonds. The H(eta)(12) and H(eta)(22) protons have large downfield shifts as expected for a symmetrical end-on interaction with the ligand carboxylate group. The chemical shifts are essentially the same in the complexes with folate and p-aminobenzoyl-L-glutamate (PABG) and similar to those found previously for the methotrexate complex reflecting the strong and similar hydrogen bonds formed with the carboxylate oxygens. Interestingly, the rates of rotation about the N(epsilon)-C(zeta) bond for the complexes containing the weakly binding PABG fragment are almost identical to those measured in the complex with methotrexate, which binds 10(7) times more tightly. In the methotrexate complex, this rotation depends on correlated rotations about the N(epsilon)-C(zeta) bond of Arg 57 and the C(alpha)-C' bond of the ligand glutamate alpha-carboxylate group. Thus, even in a fragment such as PABG, which has a much faster off-rate, the carboxylate group binds to the enzyme in a similar way to that in a parent molecule such as folate and methotrexate with the rotation about the N(epsilon)-C(zeta) bond of Arg 57 being essentially the same in all the different complexes.     

A kinetic study of wild-type and mutant dihydrofolate reductases from Lactobacillus casei

A kinetic scheme is presented for Lactobacillus casei dihydrofolate reductase that predicts steady-state kinetic parameters. This scheme was derived from measuring association and dissociation rate constants and pre-steady-state transients by using stopped-flow fluorescence and absorbance spectroscopy. Two major features of this kinetic scheme are the following: (i) product dissociation is the rate-limiting step for steady-state turnover at low pH and follows a specific, preferred pathway in which tetrahydrofolate (H4F) dissociation occurs after NADPH replaces NADP+ in the ternary complex; (ii) the rate constant for hydride transfer from NADPH to dihydrofolate (H2F) is rapid (khyd = 430 s-1), favorable (Keq = 290), and pH dependent (pKa = 6.0), reflecting ionization of a single group. Not only is this scheme identical in form with the Escherichia coli kinetic scheme [Fierke et al. (1987) Biochemistry 26, 4085] but moreover none of the rate constants vary by more than 40-fold despite there being less than 30% amino acid homology between the two enzymes. This similarity is consistent with their overall structural congruence. The role of Trp-21 of L. casei dihydrofolate reductase in binding and catalysis was probed by amino acid substitution. Trp-21, a strictly conserved residue near both the folate and coenzyme binding sites, was replaced by leucine. Two major effects of this substitution are on (i) the rate constant for hydride transfer which decreases 100-fold, becoming the rate-limiting step in steady-state turnover, and (ii) the affinities for NADPH and NADP+ which decrease by approximately 3.5 and approximately 0.5 kcal mol-1, respectively.(ABSTRACT TRUNCATED AT 250 WORDS)     

1H/15N HSQC NMR studies of ligand carboxylate group interactions with arginine residues in complexes of brodimoprim analogues and Lactobacillus casei dihydrofolate reductase

1H and 15N NMR studies have been undertaken on complexes of Lactobacillus casei dihydrofolate reductase (DHFR) formed with analogues of the antibacterial drug brodimoprim (2,4-diamino-5-(3', 5'-dimethoxy-4'-bromobenzyl)pyrimidine) in order to monitor interactions between carboxylate groups on the ligands and basic residues in the protein. These analogues had been designed by computer modeling with carboxylated alkyl chains introduced at the 3'-O position in order to improve their binding properties by making additional interactions with basic groups in the protein. Specific interactions between ligand carboxylate groups and the conserved Arg57 residue have been detected in studies of 1H/15N HSQC spectra of complexes of DHFR with both the 4-carboxylate and the 4, 6-dicarboxylate brodimoprim analogues. The spectra from both complexes showed four resolved signals for the four NHeta protons of the guanidino group of Arg57, and this is consistent with hindered rotation in the guanidino group resulting from interactions with the 4-carboxylate group in each analogue. In the spectra of each complex, one of the protons from each of the two NH2 groups and both nitrogens are considerably deshielded compared to the shielding values normally observed for such nuclei. This pattern of deshielding is that expected for a symmetrical end-on interaction of the carboxylate oxygens with the NHeta12 and NHeta22 guanidino protons. The differences in the degree of deshielding between the complexes of the two structurally similar brodimoprim analogues and the methotrexate indicates that the shielding is very sensitive to geometry, most probably to hydrogen bond lengths. The 1H/15N HSQC spectrum of the DHFR complex with the brodimoprim-6-carboxylate analogue does not feature any deshielded Arg NHeta protons and this argues against a similar interaction with the Arg57 in this case. It has not proved possible to determine whether the 6-carboxylate in this analogue is interacting directly with any residue in the protein. 1H/15N HSQC spectra have been fully assigned for the complexes with the three brodimoprim analogues and chemical shift mapping used to explore interactions in the binding site. The 1H signals of the bound ligands for all three brodimoprim analogues have been assigned. Their 1H chemical shifts were found to be fairly similar in the different complexes indicating that the 2, 4-diaminopyrimidine and the benzyl ring are binding in essentially the same binding sites and with the same overall conformation in the different complexes. The rotation rate about the NepsilonCzeta bond in the brodimoprim-4,6-dicarboxylate complex with DHFR has been determined from a zz-HSQC exchange experiment, and its value is quite similar to that observed in the DHFR.methotrexate complex (24 +/- 10 s-1 at 8 degrees C and 50 +/- 10 s-1 at 15 degrees C, respectively). The 1H and 15N chemical shift differences of selected amide and guanidino NH groups, measured between the DHFR complexes, provided further evidence about the interactions involving Arg57 with the 4-carboxylate and 4,6-dicarboxylate brodimoprim analogues.     

Comparison of solution structures of dihydrofolate reductases and enzyme-ligand complexes using cross-reacting antibodies

Polyclonal antibodies against dihydrofolate reductase (DHFR) from the human lymphoblastoid cell line WIL-2/M4 were used as probes to compare the antigenic structures in solution of native DHFRs obtained from a broad range of species and their complexes with substrate, cofactor, and folate antagonist inhibitors. All these antibodies could bind to the denatured human DHFR, indicating that they were specific for the primary structure of this enzyme. Denatured chicken liver and L1210 murine leukemic DHFRs competed for all of the antibodies that bound to the human enzyme, although less effectively than the denatured human enzyme, showing the presence of similar epitopes among the vertebrate enzymes. However, both direct binding and competition experiments showed low antibody cross-reactivities with native chicken liver (8%) and murine (10%) DHFRs, suggesting differences in the disposition of similar epitopes in these enzymes. The lactobacillus casei DHFR showed a low amount (less than 2%) of cross-reactivity with the antibodies while the same antibodies did not cross-react with the Escherichia coli enzyme. DHFR from soybean seedlings competed for a large proportion (70%) of the anti-human DHFR antibodies, indicating a close similarity in the antigenic structures of plant and animal DHFRs. Binary complexes of the L. casei, avian, murine, and human DHFRs with dihydrofolate, methotrexate (MTX), trimethoprim (TMP), NADPH, and NADP+ all showed significantly lower antibody binding capacity as compared with the corresponding free enzymes. Further, these ligands inhibited antibody binding to the enzyme to varying degrees.(ABSTRACT TRUNCATED AT 250 WORDS)     

Calculations of the conformations of Iturin A in relation with NMR studies

Iturin A is an antifungal antibiotic which was isolated from a strain of Bacillus subtilis, and contains a lipophilic beta amino acid closing an heptapeptide cycle with polar L and D residues. Iturin A belongs to a lipopeptide family of which the LDDLLDL sequence is kept constant. NMR spectroscopy and semi-empirical energy calculations are combined to design the conformations of Iturin A in pyridine solution. J coupling constants and nOes (nuclear Overhauser enhancements) are used as guiding line for energy calculations. This preliminary study shows that Iturin A in pyridine appears as rather rigid, especially in the L Pro 5-D Asn 6 region, probably involved in a beta turn. The polar side chains can form different networks of intramolecular hydrogen bonds. The Tyr side chain, relatively mobile, could be involved in interactions with an hydrophobic environment as the beta amino acid side chain found away from the peptide cycle.     

Comparative studies on dihydrofolate reductases from Plasmodium falciparum and Aotus trivirgatus.

Dihydrofolate reductase (E.C. 1.5.1.3) from Plasmodium falciparum and from its host, the owl monkey (Aotus trivirgatus), were partially purified and characterized. The molecular weight of the parasite enzyme was estimated to be over 10 times as high as that of the host enzyme. The host enzyme had 2 pH optima whereas the parasite enzyme only one. The activity of the host enzyme was greatly stimulated by KCl and urea, while that of the parasite enzyme was inhibited at high concentrations of such chaotropic agents. Km of the parasite enzyme was significantly higher than that of the host enzyme. The parasite enzyme had much lower Ki for pyrimethamine than the host enzyme. Dihydrofolate reductases isolated from pyrimethamine-resistant and pyrimethamine sensitive strains of P. falciparum were found to be similar.     

Calorimetric studies of ligand binding in R67 dihydrofolate reductase

R67 dihydrofolate reductase (DHFR) is a novel bacterial protein that possesses 222 symmetry and a single active site pore. Although the 222 symmetry implies that four symmetry-related binding sites must exist for each substrate as well as for each cofactor, various studies indicate only two molecules bind. Three possible combinations include two dihydrofolate molecules, two NADPH molecules, or one substrate plus one cofactor. The latter is the productive ternary complex. To explore the role of various ligand substituents during binding, numerous analogues, inhibitors, and fragments of NADPH and/or folate were used in both isothermal titration calorimetry (ITC) and K(i) studies. Not surprisingly, as the length of the molecule is shortened, affinity is lost, indicating that ligand connectivity is important in binding. The observed enthalpy change in ITC measurements arises from all components involved in the binding process, including proton uptake. As a buffer dependence for binding of folate was observed, this likely correlates with perturbation of the bound N3 pK(a), such that a neutral pteridine ring is preferred for pairwise interaction with the protein. Of interest, there is no enthalpic signal for binding of folate fragments such as dihydrobiopterin where the p-aminobenzoylglutamate tail has been removed, pointing to the tail as providing most of the enthalpic signal. For binding of NADPH and its analogues, the nicotinamide carboxamide is quite important. Differences between binary (binding of two identical ligands) and ternary complex formation are observed, indicating interligand pairing preferences. For example, while aminopterin and methotrexate both form binary complexes, albeit weakly, neither readily forms ternary complexes with the cofactor. These observations suggest a role for the O4 atom of folate in a pairing preference with NADPH, which ultimately facilitates catalysis.     

NMR studies of the heme pocket conformations of monomeric hemoglobins from Glycera dibranchiata. Implications for ligand binding

Two-dimensional 1H-NMR methods have been used to assign side-chain resonances for the tryptophan residues and for several amino acids located in the heme pockets of the carbon monoxide complexes of the major monomeric hemoglobins from Glycera dibranchiata. The NMR spectra reveal a high degree of conservation of the heme pocket structure in the different hemoglobins. However some conformational differences are evident and residues at positions B10 and G8 on the distal side of the heme pocket are not conserved. From the present NMR studies it appears that the monomeric G. dibranchiata hemoglobin examined by X-ray crystallography [Padlan, E. A. & Love, W. (1974) J. Biol. Chem. 249, 4067-4078] corresponds to HbC. Except that the orientation of the heme in solution is the reverse of that reported in the crystal structure, there is a close correspondence between the heme pocket structure in the crystal and in solution. The proximal histidine coordination geometry is almost identical in the CO complexes of the three monomeric hemoglobins studied. Distal residues are strongly implicated in determining the observed kinetic differences in ligand binding reactions. In particular, steric crowding of the ligand binding site in hemoglobin A is probably a major factor in the slower kinetics of this component.     

1H NMR studies of homo and mixed ligand complexes of Tl+ ion with several polyazamacrocycles

Proton NMR was used as a probe to study the interaction of the Tl(+) ion with 9-18-membered macromonocyclic tri-, tetra-, and hexaamines in dimethylformamide (DMF) solution. A study of proton chemical shift of ligands as a function of Tl(+) ion to ligand mole ratio revealed that the complexation reactions occur in a stepwise manner. Formation of a 1:1 complex is followed by the addition of a second complexant molecule to form a homo-sandwich complex for triazamacrocycle ligands and a mixed ligand complex in the case of hexamethylhexacyclen (HMHCY) and 1,4,7-triazacyclononane ([9]aneN(3)). The formation constants of resulting 1:1 and 1:2 (homo and mixed ligand sandwich) complexes in DMF solution were evaluated from computer fitting of the chemical shift-mole ratio data. The mixed ligand complexes may be more stable than the parent complex in which both ligands are the same. The influence of cavity size and substitution of methyl groups on nitrogen atoms of the macrocyclic ring the stability of the resulting complexes is discussed. The geometries of the tri- and tetraazamacrocycle ligands and their Tl(+) ion complexes were optimized by an ab initio method, and the calculated binding energies of resulting complexes were compared. Both the experimental and theoretical studies revealed that, in the presence of methyl groups, the stability of triazamacrocycle complexes with Tl(+) ion was decreased.     

NMR studies of the conformations of leghemoglobins from soybean and lupin

Phase-sensitive two-dimensional NMR methods have been used to obtain extensive proton resonance assignments for the carbon monoxide complexes of lupin leghemoglobins I and II and soybean leghemoglobin a. The assigned resonances provide information on the solution conformations of the proteins, particularly in the vicinity of the heme. The structure of the CO complex of lupin leghemoglobin II in solution is compared with the X-ray crystal structure of the cyanide complex by comparison of observed and calculated ring current shifts. The structures are generally very similar but significant differences are observed for the ligand contact residues, Phe30, His63 and Val67, and for the proximal His97 ligand. Certain residues are disordered and adopt two interconverting conformations in lupin leghemoglobin II in solution. The proximal heme pocket structure is closely conserved in the lupin leghemoglobins I and II but small differences in conformation in the distal heme pocket are apparent. Larger conformational differences are observed when comparisons are made with the CO complex of soybean leghemoglobin. Altered protein-heme packing is indicated on the proximal side of the heme and some conformational differences are evident in the distal heme pocket. The small conformational differences between the three leghemoglobins probably contribute to the known differences in their O2 and CO association and dissociation kinetics. The heme pocket conformations of the three leghemoglobins are more closely related to each other than to sperm whale myoglobin. The most notable differences between the leghemoglobins and myoglobin are: (a) reduced steric crowding of the ligand binding site in the leghemoglobins, (b) different orientations of the distal histidine, and (c) small but significant differences in proximal histidine coordination geometry. These changes probably contribute to the large differences in ligand binding kinetics between the leghemoglobins and myoglobin.     

13C NMR evidence of the slow exchange of tryptophans in dihydrofolate reductase between stable conformations.

¹³C NMR spectra are reported for dihydrofolate reductase of $\text{Streptococcus}\text{faecium}$ labeled with [γ-¹³C]tryptophan. Two of the four tryptophans generate unusual resonances indicating slow exchange of the residues between alternative stable conformations. Since 3′, 5′-dichloromethotrexate sharpens two of the resonances, it apparently locks the corresponding residues into one conformation.     

Volume and compressibility differences between protein conformations revealed by high-pressure NMR

Proteins often interconvert between different conformations in ways critical to their function. Although manipulating such equilibria for biophysical study is often challenging, the application of pressure is a potential route to achieve such control by favoring the population of lower volume states. Here, we use this feature to study the interconversion of ARNT PAS-B Y456T, which undergoes a dramatic +3 slip in the β-strand register as it switches between two stably folded conformations. Using high-pressure biomolecular NMR approaches, we obtained the first, to our knowledge, quantitative data testing two key hypotheses of this process: the slipped conformation is both smaller and less compressible than the wild-type equivalent, and the interconversion proceeds through a chiefly unfolded intermediate state. Data collected in steady-state pressure and time-resolved pressure-jump modes, including observed pressure-dependent changes in the populations of the two conformers and increased rate of interconversion between conformers, support both hypotheses. Our work exemplifies how these approaches, which can be generally applied to protein conformational switches, can provide unique information that is not easily accessible through other techniques.     

Interactions of 5-deazapteridine derivatives with Mycobacterium tuberculosis and with human dihydrofolate reductases

There are major differences between the structures of human dihydrofolate reductase (hDHFR) and Mycobacterium tuberculosis dihydrofolate reductase (mtDHFR). These differences may allow us to design more selective mtDHFR inhibitors. In this paper we study the reactions of six different compounds derived from 5-deazapteridine with human and bacterial enzymes. Results suggest that the addition of hydrophobic groups to the aminophenyl ring would increase mtDHFR-inhibitor affinity and selectivity.     

15N and 1H NMR evidence for multiple conformations of the complex of dihydrofolate reductase with its substrate, folate

The binding of folate to Lactobacillus casei dihydrofolate reductase in the presence and absence of NADP+ has been studied by 15N NMR, using [5-15N]folate. In the presence of NADP+, three separate signals were observed for the single 15N atom, in agreement with our earlier evidence from 1H and 13C NMR for multiple conformations of this complex [(1982) Biochemistry 21, 5831-5838]. The 15N spectra of the binary enzyme-folate complex provide evidence for the first time that this complex also exists in at least two conformational states. This is confirmed by the observation of two separate resonances for the 7-proton of bound folate, located by two-dimensional exchange spectroscopy.     

Solution conformations and dynamics of ABL kinase-inhibitor complexes determined by NMR substantiate the different binding modes of imatinib/nilotinib and dasatinib

Current structural understanding of kinases is largely based on x-ray crystallographic studies, whereas very little data exist on the conformations and dynamics that kinases adopt in the solution state. ABL kinase is an important drug target in the treatment of chronic myelogenous leukemia. Here, we present the first characterization of ABL kinase in complex with three clinical inhibitors (imatinib, nilotinib, and dasatinib) by modern solution NMR techniques. Structural and dynamical results were derived from complete backbone resonance assignments, experimental residual dipolar couplings, and (15)N relaxation data. Residual dipolar coupling data on the imatinib and nilotinib complexes show that the activation loop adopts the inactive conformation, whereas the dasatinib complex preserves the active conformation, which does not support contrary predictions based upon molecular modeling. Nanosecond as well as microsecond dynamics can be detected for certain residues in the activation loop in the inactive and active conformation complexes.