Large-scale purification and characterization of dihydrofolate reductase from a methotrexate-resistant strain of Lactobacillus casei

The impact of the human gut microbiome on the treatment of autoimmune disease

Autoimmune (or rheumatic) diseases are increasing in prevalence but selecting the best therapy for each patient proceeds in trial-and-error fashion. This strategy can lead to ineffective therapy resulting in irreversible damage and suffering; thus, there is a need to bring the promise of precision medicine to patients with autoimmune disease. While host factors partially determine the therapeutic response to immunosuppressive drugs, these are not routinely used to tailor therapy. Thus, non-host factors likely contribute. Here, we consider the impact of the human gut microbiome in the treatment of autoimmunity. We propose that the gut microbiome can be manipulated to improve therapy and to derive greater benefit from existing therapies. We focus on the mechanisms by which the human gut microbiome impacts treatment response, provide a framework to interrogate these mechanisms, review a case study of a widely-used anti-rheumatic drug, and discuss challenges with studying multiple complex systems: the microbiome, the human immune system, and autoimmune disease. We consider open questions that remain in the field and speculate on the future of drug-microbiome-autoimmune disease interactions. Finally, we present a blue-sky vision for how the microbiome can be used to bring the promise of precision medicine to patients with rheumatic disease.

Reduction of Folate by Dihydrofolate Reductase from Thermotoga maritima

Mammalian dihydrofolate reductases (DHFRs) catalyze the reduction of folate more efficiently than the equivalent bacterial enzymes do, despite typically having similar efficiencies for the reduction of their natural substrate, dihydrofolate. In contrast, we show here that DHFR from the hyperthermophilic bacterium Thermotoga maritima can catalyze reduction of folate to tetrahydrofolate with an efficiency similar to that of reduction of dihydrofolate under saturating conditions. Nuclear magnetic resonance and mass spectrometry experiments showed no evidence of the production of free dihydrofolate during either the EcDHFR- or TmDHFR-catalyzed reductions of folate, suggesting that both enzymes perform the two reduction steps without release of the partially reduced substrate. Our results imply that the reaction proceeds more efficiently in TmDHFR than in EcDHFR because the more open active site of TmDHFR facilitates protonation of folate. Because T. maritima lives under extreme conditions where tetrahydrofolate is particularly prone to oxidation, this ability to salvage folate may impart an advantage to the bacterium by minimizing the squandering of a valuable cofactor.

NMR structures of apo L. casei dihydrofolate reductase and its complexes with trimethoprim and NADPH: contributions to positive cooperative binding from ligand-induced refolding, conformational changes, and interligand hydrophobic interactions

In order to examine the origins of the large positive cooperativity (ΔG(0)(coop) = -2.9 kcal mol(-1)) of trimethoprim (TMP) binding to a bacterial dihydrofolate reductase (DHFR) in the presence of NADPH, we have determined and compared NMR solution structures of L. casei apo DHFR and its binary and ternary complexes with TMP and NADPH and made complementary thermodynamic measurements. The DHFR structures are generally very similar except for the A-B loop region and part of helix B (residues 15-31) which could not be directly detected for L. casei apo DHFR because of line broadening from exchange between folded and unfolded forms. Thermodynamic and NMR measurements suggested that a significant contribution to the cooperativity comes from refolding of apo DHFR on binding the first ligand (up to -0.95 kcals mol(-1) if 80% of A-B loop requires refolding). Comparisons of Cα-Cα distance differences and domain rotation angles between apo DHFR and its complexes indicated that generally similar conformational changes involving domain movements accompany formation of the binary complexes with either TMP or NADPH and that the binary structures are approaching that of the ternary complex as would be expected for positive cooperativity. These favorable ligand-induced structural changes upon binding the first ligand will also contribute significantly to the cooperative binding. A further substantial contribution to cooperative binding results from the proximity of the bound ligands in the ternary complex: this reduces the solvent accessible area of the ligand and provides a favorable entropic hydrophobic contribution (up to -1.4 kcal mol(-1)).

Kinetic and structural characterization of dihydrofolate reductase from Streptococcus pneumoniae

Drug resistance associated with dihydrofolate reductase (DHFR) has emerged as a critical issue in the treatment of bacterial infections. In our efforts to understand the mechanism of a drug-resistant dihydrofolate reductase (DHFR) from a pathogenic bacterial source, we report the first kinetic characterization of Streptococcus pneumoniae DHFR (spDHFR) along with its X-ray structure. This study revealed that the kinetic properties of spDHFR were significantly different from those of Escherichia coli DHFR. The product (tetrahydrofolate) dissociation step that is the rate-limiting step in E. coli DHFR is significantly accelerated in spDHFR so that hydride transfer or a preceding step is rate-limiting. Comparison of the binding parameters of this enzyme to those of a mutant spDHFR (Sp9) confirmed that the Leu100 residue in spDHFR is the critical element for the trimethoprim (TMP) resistance. Steady-state kinetics exhibited a pH dependence in k(cat), which prompted us to elucidate the role of the new catalytic residue (His33) in the active site of spDHFR. Structural data of the Sp9 mutant in complex with NADPH and methotrexate confirmed the participation of His33 in a hydrogen bonding network involving a water molecule, the hydroxyl group of Thr119, and the carboxylate ion of Glu30. Sequence analysis of the DHFR superfamily revealed that the His residue is the major amino acid component at this position and is found mostly in pathogenic bacterial DHFRs. A mutation of Val100 to Leu demonstrated a steric clash of the leucine side chain with the side chains of Ile8 and Phe34, rationalizing weaker binding of trimethoprim to Leu100 DHFR. Understanding the role of specific amino acids in the active site coupled with detailed structural analysis will inform us on how to better design inhibitors targeting drug-resistant pathogenic bacterial DHFRs.

Crystallization and preliminary crystallographic studies of dihydrofolate reductase-thymidylate synthase from Trypanosoma cruzi, the Chagas disease pathogen

Trypanosoma cruzi dihydrofolate reductase-thymidylate synthase (TcDHFR-TS) was crystallized in complexes with the dihydrotriazine-based or quinazoline-based antifolates C-448, cycloguanil (CYC) and Q-8 in order to gain insight into the interactions of this DHFR enzyme with classical and novel inhibitors. The TcDHFR-TS-C-448-NDP-dUMP crystal belonged to space group C222(1) with two molecules per asymmetric unit and diffracted to 2.37 angstrom resolution. The TcDHFR-TS-CYC, TcDHFR-TS-CYC-NDP and TcDHFR-TS-Q-8-NDP crystals belonged to space group P2(1) with four molecules per asymmetric unit and diffracted to 2.1, 2.6 and 2.8 angstrom resolution, respectively. Crystals belonging to the two different space groups were suitable for structure determination.

Crystal structure of Bacillus anthracis dihydrofolate reductase with the dihydrophthalazine-based trimethoprim derivative RAB1 provides a structural explanation of potency and selectivity

Bacillus anthracis possesses an innate resistance to the antibiotic trimethoprim due to poor binding to dihydrofolate reductase (DHFR); currently, there are no commercial antibacterials that target this enzyme in B. anthracis. We have previously reported a series of dihydrophthalazine-based trimethoprim derivatives that are inhibitors for this target. In the present work, we have synthesized one compound (RAB1) displaying favorable 50% inhibitory concentration (54 nM) and MIC (< or =12.8 microg/ml) values. RAB1 was cocrystallized with the B. anthracis DHFR in the space group P2(1)2(1)2(1), and X-ray diffraction data were collected to a 2.3-A resolution. Binding of RAB1 causes a conformational change of the side chain of Arg58 and Met37 to accommodate the dihydrophthalazine moiety. Unlike the natural substrate or trimethoprim, the dihydrophthalazine group provides a large hydrophobic anchor that embeds within the DHFR active site and accounts for its selective inhibitory activity against B. anthracis.

Origins of specificity in the binding of small molecules to dihydrofolate reductase

Dihydrofolate reductase is the target for the therapeutically important 'anti-folate' drugs such as methotrexate and trimethoprim. Methotrexate is a powerful inhibitor of the enzyme, binding up to 10,000 times more tightly than the structurally similar substrate, folate. Two contributions to this striking difference in affinity have been identified: the two ligands bind in different charge states, and there are conformational differences between the two complexes. The origins of the tight binding of methotrexate have been explored further by studying the binding of 2,4-diaminopyrimidine and p-aminobenzoyl-L-glutamate, which may be considered as 'fragments' of methotrexate. These two compounds bind simultaneously but also cooperatively, the binding of one 'fragment' leading to a 50-fold increase in the affinity for the other. Studies of structural analogues of these fragments show that the specificity as well as the strength of binding can be altered by the presence of the other 'fragment'; both positive and negative cooperativity are observed. The relation of these observations to methotrexate binding, and the notion of intramolecular cooperativity in ligand binding are discussed.

Discovery of potent pteridine reductase inhibitors to guide antiparasite drug development

Pteridine reductase (PTR1) is essential for salvage of pterins by parasitic trypanosomatids and is a target for the development of improved therapies. To identify inhibitors of Leishmania major and Trypanosoma cruzi PTR1, we combined a rapid-screening strategy using a folate-based library with structure-based design. Assays were carried out against folate-dependent enzymes including PTR1, dihydrofolate reductase (DHFR), and thymidylate synthase. Affinity profiling determined selectivity and specificity of a series of quinoxaline and 2,4-diaminopteridine derivatives, and nine compounds showed greater activity against parasite enzymes compared with human enzymes. Compound 6a displayed a K(i) of 100 nM toward LmPTR1, and the crystal structure of the LmPTR1:NADPH:6a ternary complex revealed a substrate-like binding mode distinct from that previously observed for similar compounds. A second round of design, synthesis, and assay produced a compound (6b) with a significantly improved K(i) (37 nM) against LmPTR1, and the structure of this complex was also determined. Biological evaluation of selected inhibitors was performed against the extracellular forms of T. cruzi and L. major, both wild-type and overexpressing PTR1 lines, as a model for PTR1-driven antifolate drug resistance and the intracellular form of T. cruzi. An additive profile was observed when PTR1 inhibitors were used in combination with known DHFR inhibitors, and a reduction in toxicity of treatment was observed with respect to administration of a DHFR inhibitor alone. The successful combination of antifolates targeting two enzymes indicates high potential for such an approach in the development of previously undescribed antiparasitic drugs.

Synthesis of N-(5,7-diamino-3-phenyl-quinoxalin-2-yl)-3,4,5-substituted anilines and N-[4[(5,7-diamino-3-phenylquinoxalin-2-yl)amino]benzoyl]-l-glutamic acid diethyl ester: Evaluation of in vitro anti-cancer and anti-folate activities

Several diamino quinoxalines were designed, synthesized and evaluated as anti-tumor agents. Two compounds showed the most potent cytotoxic activities against the leukemia CCRF-CEM cell line (GI(50)<0.01microM) and the ovarian cancer cell line OVCAR-4 (GI(50)=0.03microM), respectively, with comparable/better activities than Methotrexate (MTX). Docking calculations of the complexes of hDHFR with the most active compounds identified the binding mode of the described molecules with respect to MTX.

Cloning and characterization of dihydrofolate reductase from a facultative alkaliphilic and halotolerant bacillus strain

Elucidation of the molecular basis of the stability of enzymes from extremophilic organisms is of fundamental importance for various industrial applications. Due to the wealth of structural data from various species, dihydrofolate reductase (DHFR, EC 1.5.1.3) provides an excellent model for systematic investigations. In this report, DHFR from alkaliphilic Bacillus halodurans C-125 was cloned and expressed in E. coli. Functional analyses revealed that BhDHFR exhibits the most alkali-stable phenotype of DHFRs characterized so far. Optimal enzyme activity was observed in a slightly basic pH region ranging from 7.25 to 8.75. Alkali-stability is associated with a remarkable resistance to elevated temperatures (half-life of 60 min at 52.5 degrees C) and to high concentrations of urea (up to 3 M). Although the secondary structure shows distinct similarities to those of mesophilic DHFR molecules, BhDHFR exhibits molecular features contributing to its alkaliphilic properties. Interestingly, the unique phenotype is diminished by C-terminal addition of a His-tag sequence. Therefore, His-tag-derivatized BhDHFR offers the opportunity to obtain deeper insights into the specific mechanisms of alkaliphilic adaption by comparison of the three dimensional structure of both BhDHFR molecules.

Effects of co-operative ligand binding on protein amide NH hydrogen exchange

Amide protection factors have been determined from NMR measurements of hydrogen/deuterium amide NH exchange rates measured on assigned signals from Lactobacillus casei apo-DHFR and its binary and ternary complexes with trimethoprim (TMP), folinic acid and coenzymes (NADPH/NADP(+)). The substantial sizes of the residue-specific DeltaH and TDeltaS values for the opening/closing events in NH exchange for most of the measurable residues in apo-DHFR indicate that sub-global or global rather than local exchange mechanisms are usually involved. The amide groups of residues in helices and sheets are those most protected in apo-DHFR and its complexes, and the protection factors are generally related to the tightness of ligand binding. The effects of ligand binding that lead to changes in amide protection are not localised to specific binding sites but are spread throughout the structure via a network of intramolecular interactions. Although the increase in protein stability in the DHFR.TMP.NADPH complex involves increased ordering in the protein structure (requiring TDeltaS energy) this is recovered, to a large extent, by the stronger binding (enthalpic DeltaH) interactions made possible by the reduced motion in the protein. The ligand-induced protection effects in the ternary complexes DHFR.TMP.NADPH (large positive binding co-operativity) and DHFR.folinic acid.NADPH (large negative binding co-operativity) mirror the co-operative effects seen in the ligand binding. For the DHFR.TMP.NADPH complex, the ligand-induced protection factors result in DeltaDeltaG(o) values for many residues being larger than the DeltaDeltaG(o) values in the corresponding binary complexes. In contrast, for DHFR.folinic acid.NADPH, the DeltaDeltaG(o) values are generally smaller than many of those in the corresponding binary complexes. The results indicate that changes in protein conformational flexibility on formation of the ligand complex play an important role in determining the co-operativity in the ligand binding.

Biochemical characterization of a dihydromethanopterin reductase involved in tetrahydromethanopterin biosynthesis in Methylobacterium extorquens AM1

During growth on one-carbon (C1) compounds, the aerobic alpha-proteobacterium Methylobacterium extorquens AM1 synthesizes the tetrahydromethanopterin (H4MPT) derivative dephospho-H4MPT as a C1 carrier in addition to tetrahydrofolate. The enzymes involved in dephospho-H4MPT biosynthesis have not been identified in bacteria. In archaea, the final step in the proposed pathway of H4MPT biosynthesis is the reduction of dihydromethanopterin (H2MPT) to H4MPT, a reaction analogous to the reaction of the bacterial dihydrofolate reductase. A gene encoding a dihydrofolate reductase homolog has previously been reported for M. extorquens and assigned as the putative H2MPT reductase gene (dmrA). In the present work, we describe the biochemical characterization of H2MPT reductase (DmrA), which is encoded by dmrA. The gene was expressed with a six-histidine tag in Escherichia coli, and the recombinant protein was purified by nickel affinity chromatography and gel filtration. Purified DmrA catalyzed the NAD(P)H-dependent reduction of H2MPT with a specific activity of 2.8 micromol of NADPH oxidized per min per mg of protein at 30 degrees C and pH 5.3. Dihydrofolate was not a substrate for DmrA at the physiological pH of 6.8. While the existence of an H2MPT reductase has been proposed previously, this is the first biochemical evidence for such an enzyme in any organism, including archaea. Curiously, no DmrA homologs have been identified in the genomes of known methanogenic archaea, suggesting that bacteria and archaea produce two evolutionarily distinct forms of dihydromethanopterin reductase. This may be a consequence of different electron donors, NAD(P)H versus reduced F420, used, respectively, in bacteria and methanogenic archaea.

Purification and characterization of NAD-specific 6-phosphogluconate dehydrogenase from Leuconostoc lactis SHO-54

The 6-phosphogluconate dehydrogenase (EC 1.1.1.44) from Leuconostoc lactis SHO-54 was purified with an overall yield of 38% and a specific activity of 140.0 units/mg protein. The enzyme had a tetrameric structure and a molecular mass of 32.8 kDa. The amino acid composition of the purified enzyme was determined, and the enzyme contained no sulfhydryl amino acids. The K(m) values for 6-phosphogluconate and NAD were 0.95 mM and 0.32 mM, respectively.

Moritella cold-active dihydrofolate reductase: are there natural limits to optimization of catalytic efficiency at low temperature?

Adapting metabolic enzymes of microorganisms to low temperature environments may require a difficult compromise between velocity and affinity. We have investigated catalytic efficiency in a key metabolic enzyme (dihydrofolate reductase) of Moritella profunda sp. nov., a strictly psychrophilic bacterium with a maximal growth rate at 2 degrees C or less. The enzyme is monomeric (Mr=18,291), 55% identical to its Escherichia coli counterpart, and displays Tm and denaturation enthalpy changes much lower than E. coli and Thermotoga maritima homologues. Its stability curve indicates a maximum stability above the temperature range of the organism, and predicts cold denaturation below 0 degrees C. At mesophilic temperatures the apparent Km value for dihydrofolate is 50- to 80-fold higher than for E. coli, Lactobacillus casei, and T. maritima dihydrofolate reductases, whereas the apparent Km value for NADPH, though higher, remains in the same order of magnitude. At 5 degrees C these values are not significantly modified. The enzyme is also much less sensitive than its E. coli counterpart to the inhibitors methotrexate and trimethoprim. The catalytic efficiency (kcat/Km) with respect to dihydrofolate is thus much lower than in the other three bacteria. The higher affinity for NADPH could have been maintained by selection since NADPH assists the release of the product tetrahydrofolate. Dihydrofolate reductase adaptation to low temperature thus appears to have entailed a pronounced trade-off between affinity and catalytic velocity. The kinetic features of this psychrophilic protein suggest that enzyme adaptation to low temperature may be constrained by natural limits to optimization of catalytic efficiency.

Quinoxaline chemistry. Part 15. 4-[2-Quinoxalylmethylenimino]-benzoylglutamates and -benzoates, 4-[2-quinoxalylmethyl-N-methylamino]-benzoylglutamates as analogues of classical antifolate agents. Synthesis, elucidation of structures and in vitro evaluation of antifolate and anticancer activities

We report on an extension of our previous discovery of in vitro anticancer activity of trifluoromethylquinoxalines as analogues of classical and non-classical antifolic methotrexate and trimetrexate. In this case a small number of Schiff bases were obtained from the reaction of 2-bromethyl-3-R-6(7)trifluoromethylquinoxaline and ethyl p-aminobenzoylglutamate, ethyl p-aminobenzoate, p-toluidine instead of the expected 4-[2-quinoxalyl]methyl-N-methylanilino derivatives, which in turn formed with N-methylanilino derivatives. The reaction mechanism has been put forward. Structure elucidation of both Schiff bases and N-methylanilino analogues was achieved by a combination of 1H and 13C NMR spectra and hetcor experiments. Compounds 3a, 3b, 3c, 8, 11, 12, 13, Ie were tested in antifolic enzyme assay [Lactobacillus casei (LcTS), Leishmania major (LmTs), human Thymidylate synthase (hTs), human TS, human dihydrofolate reductase (hDHFR)] while compounds 3a, 3b, 3c were tested for anticancer activity. These results seem to indicate that the Schiff bases are somewhat active either as anticancer or as folate inhibitors, while compound Ie was selectively active against hDHFR with an inhibition constant (Ki) of 200 nM with a specificity of about 1000-folds with respect to hTS.

10-Formyl-dihydrofolic acid is bioactive in human leukemia cells

The bioactivity of 10-formyl-7,8-dihydrofolic acid and 10-formyl-folic acid was determined in human leukemia (CCRF-CEM) cells grown in a folate-depleted medium containing methotrexate. Excess 10-formyl-7,8-dihydrofolic acid, (but not 10-formyl folic acid) supported the growth of these cells, but it was less potent than5-formyl-5,6,7,8-tetrahydrofolic acid (a control). 10-formyl-7, 8-dihydrofolic acid (not 10-formyl folic acid) was active as substrate for aminoimidazole carboxamide ribotide transformylase and dihydrofolate reductase. This is the first experimental evidence that 10-formyl-7,8-dihydrofolic acid is a bioactive folate in mammalian cells. These experiments and several other lines of evidence in the literature suggest that 10-formyl-folic acid must be metabolized to bioactive folate by enteric bacteria before it can be utilized by the vertebrate host.

Structure and dynamics in solution of the complex of Lactobacillus casei dihydrofolate reductase with the new lipophilic antifolate drug trimetrexate

We have determined the three-dimensional solution structure of the complex of Lactobacillus casei dihydrofolate reductase and the anticancer drug trimetrexate. Two thousand seventy distance, 345 dihedral angle, and 144 hydrogen bond restraints were obtained from analysis of multidimensional NMR spectra recorded for complexes containing 15N-labeled protein. Simulated annealing calculations produced a family of 22 structures fully consistent with the constraints. Several intermolecular protein-ligand NOEs were obtained by using a novel approach monitoring temperature effects of NOE signals resulting from dynamic processes in the bound ligand. At low temperature (5 degrees C) the trimethoxy ring of bound trimetrexate is flipping sufficiently slowly to give narrow signals in slow exchange, which give good NOE cross peaks. At higher temperature these broaden and their NOE cross peaks disappear thus allowing the signals in the lower-temperature spectrum to be identified as NOEs involving ligand protons. The binding site for trimetrexate is well defined and this was compared with the binding sites in related complexes formed with methotrexate and trimethoprim. No major conformational differences were detected between the different complexes. The 2,4-diaminopyrimidine-containing moieties in the three drugs bind essentially in the same binding pocket and the remaining parts of their molecules adapt their conformations such that they can make effective van der Waals interactions with essentially the same set of hydrophobic amino acids, the side-chain orientations and local conformations of which are not greatly changed in the different complexes (similar chi1 and chi2 values).

The solution structure of the complex of Lactobacillus casei dihydrofolate reductase with methotrexate

We have determined the three-dimensional solution structure of the complex of Lactobacillus casei dihydrofolate reductase (18.3 kDa, 162 amino acid residues) formed with the anticancer drug methotrexate using 2531 distance, 361 dihedral angle and 48 hydrogen bond restraints obtained from analysis of multidimensional NMR spectra. Simulated annealing calculations produced a family of 21 structures fully consistent with the constraints. The structure has four alpha-helices and eight beta-strands with two other regions, comprising residues 11 to 14 and 126 to 127, also interacting with each other in a beta-sheet manner. The methotrexate binding site is very well defined and the structure around its glutamate moiety was improved by including restraints reflecting the previously determined specific interactions between the glutamate alpha-carboxylate group with Arg57 and the gamma-carboxylate group with His28. The overall fold of the binary complex in solution is very similar to that observed in the X-ray studies of the ternary complex of L. casei dihydrofolate reductase formed with methotrexate and NADPH (the structures of the binary and ternary complexes have a root-mean-square difference over the backbone atoms of 0.97 A). Thus no major conformational change takes place when NADPH binds to the binary complex. In the binary complex, the loop comprising residues 9 to 23 which forms part of the active site has been shown to be in the "closed" conformation as defined by M. R. Sawaya & J. Kraut, who considered the corresponding loops in crystal structures of complexes of dihydrofolate reductases from several organisms. Thus the absence of the NADPH does not result in the "occluded" form of the loop as seen in crystal studies of some other dihydrofolate reductases in the absence of coenzyme. Some regions of the structure in the binary complex which form interaction sites for NADPH are less well defined than other regions. However, in general terms, the NADPH binding site appears to be essentially pre-formed in the binary complex. This may contribute to the tighter binding of coenzyme in the presence of methotrexate.

NMR detection of arginine-ligand interactions in complexes of Lactobacillus casei dihydrofolate reductase

1H-NMR and 15N-NMR signal assignments have been made for the eight arginine residues in Lactobacillus casei dihydrofolate reductase in its binary complex with methotrexate and in its ternary complex with methotrexate and NADPH. 1H-NMR chemical shifts for the guanidino groups of two of the arginines (Arg57 and Arg43) were sensitive to different modes of binding of the guanidino groups with charged oxygen atoms of the ligands. In the complexes formed with methotrexate, Arg57 showed four non-equivalent NH eta proton signals indicating hindered rotation about the N epsilon-C zeta and C zeta-N eta bonds. The NH eta 12 and NH eta 22 protons showed large downfield shifts, which would be expected for a symmetric end-on interaction of these protons with the charged oxygen atoms of a carboxylate group in methotrexate. These effects were not observed for the complex formed with trimethoprim, which does not contain any carboxylate groups. In the complex formed with NADPH present, Arg43 showed a large downfield chemical shift for its NH epsilon proton and a retardation of its rate of exchange with water. This pattern of deshielding contrasts with that detected for Arg57 and is that expected for a side-on interaction of the guanidino group protons with charged oxygen atoms of the ribose 2'-phosphate group of NADPH.

31P solid-state NMR measurements used to detect interactions between NADPH and water and to determine the ionisation state of NADPH in a protein-ligand complex subjected to low-level hydration

31P-NMR spectra of NADPH and NADPH bound to Lactobacillus casei dihydrofolate reductase have been recorded using the techniques of cross-polarization, magic-angle spinning and high-power proton-decoupling on both lyophilized and hydrated samples. Previous studies on the lyophilized complex of L. casei dihydrofolate reductase with NADPH and methotrexate, measuring the isotropic shifts and principal components of the chemical shift tensors, have shown that the 2'-phosphate group of bound NADPH exists as a mixture of the dianionic and monoanionic states [Gerothanassis, I. P, Barrie, P. J., Birdsall, B. & Feeney, J. (1994) Eur J. Biochem. 226, 211-218]. In the present study on hydrated samples, the characterization of the isotropic shift and chemical shift tensors of the 2'-phosphate signal indicates that the 2'-phosphate is almost exclusively in the dianionic state. This is in agreement with earlier 31P-NMR studies in solution [Feeney, J., Birdsall, B., Roberts, G. C. K. & Burgen, A. S. V. (1975) Nature 257, 564-566]. In experiments examining progressively hydrated (6%, 12%, 15%, by mass) samples, the observed signals become increasingly narrower probably because the microenvironments of the 31P nuclei become more homogeneous upon sample hydration. Chemical exchange between mobile water molecules and bound protons close to individual sites on NADPH has been indirectly monitored on a hydrated sample (15% water, by mass) using a pulse sequence proposed by Harbison and coworkers [Harbison, G. S., Roberts, J. E., Herzfeld, J. & Griffin, R. G. (1988) J. Am. Chem. Soc. 110, 7221-7223]. In this experiment, the two diphosphate signals are totally suppressed while the 2'-phosphate phosphorus signal remains: this indicates a significant polarization of the 2'-phosphate nuclei from protons in exchange with those of mobile water molecules.

31P-NMR studies of NADPH, NADP+ and the complex of NADPH and methotrexate with Lactobacillus casei dihydrofolate reductase in the solid state

31P-NMR spectra on solid samples of NADP+, NADPH and NADPH bound to Lactobacillus casei dihydrofolate reductase have been recorded using the techniques of cross polarisation, magic angle spinning and high power proton decoupling. The isotropic chemical shifts, the principal components of the shielding tensors and the asymmetry parameters for the 31P nuclei in the 2'-phosphate and pyrophosphate groups have been measured. The isotropic shifts show similar trends to the chemical shifts measured in solution. The isotropic shifts and the shielding tensors for the dianionic and monoanionic states of the 2'-phosphate group have been determined and the presence of both ionisation states has been detected in a solid sample of the lyophilised complex of L. casei dihydrofolate reductase with NADPH and methotrexate. This contrasts with the behaviour in solution, where only the dianionic form is bound to the enzyme. The signals from the two pyrophosphates 31P nuclei in bound NADPH were resolved and identified. The asymmetry parameters in the different ionisation states and the orientations of the shielding tensors within the molecular framework are considered in the context of previous 31P studies on phosphate-containing compounds.

3H-n.m.r. studies of multiple conformations and dynamic processes in complexes of folate and methotrexate with Lactobacillus casei dihydrofolate reductase

[7,3',5'-3H3]- and [7,9-3H3]-folic acid and [7,3',5'-3H3]methotrexate (MTX) have been prepared and 3H-n.m.r. spectra obtained for their complexes with Lactobacillus casei dihydrofolate reductase (DHFR). The 3H results confirm the presence of three pH-dependent different conformational forms in the complex DHFR.NADP+.folate. The folate benzoyl ring could be shown to be in essentially the same environment in the different forms, with the major differences being associated with the pterin ring. The appearance of a single resonance for the 3',5'-tritons showed that the benzoyl ring is flipping rapidly in all three forms. In contrast, the MTX complex was shown to exist as a single conformational state with the benzoyl ring flipping rate being too low to give a single averaged signal for the 3',5'-nuclei over the temperature range 283-313 K.

High expression in Escherichia coli of the gene coding for dihydrofolate reductase of the extremely halophilic archaebacterium Haloferax volcanii. Reconstitution of the active enzyme and mutation studies

The gene coding for the enzyme dihydrofolate reductase of the extremely halophilic archaebacterium Haloferax volcanii was recombined into the Escherichia coli expression vector pET11d. Following induction, the enzyme was produced in large quantities and accumulated in the cells in an insoluble form. The enzymic activity could be efficiently reconstituted by dissolving the aggregate in 6 M guanidine hydrochloride followed by dilution into salt solutions. Mutants were produced in which Lys30 was converted to Leu (K30L), Lys31 was converted to Ala (K31A) and a double mutant in which both lysines were converted (K30L, K31A). The mutated enzymes were produced in E. coli, activated and purified to homogeneity. The effect of the salt concentration on the steady-state kinetic parameters was determined. It was found that the salt concentration affects the Km but not kcat of the various mutants.

Stereospecific assignments of the leucine methyl resonances in the 1H NMR spectrum of Lactobacillus casei dihydrofolate reductase

A general method is described for the stereospecific assignment of methyl resonances in protein NMR spectra based on selective deuteration procedures. A selectively deuterated dihydrofolate reductase from L. casei was prepared by incorporating stereoselectively deuterated L-leucine, (2S,4R)[5,5,5-2H3]leucine. By comparing the COSY spectra of the dihydrofolate reductase-methotrexate complexes formed using deuterated and non-deuterated enzyme the stereospecific assignments for resonances of all 13 leucine residues were obtained by noting the absence of cross-peaks in spectra from the deuterated proteins.

Nuclear magnetic resonance detection of bound water molecules in the active site of Lactobacillus casei dihydrofolate reductase in aqueous solution

Proton nuclear magnetic resonance spectroscopy has been used to detect two water molecules bound to residues in the active site of the Lactobacillus casei dihydrofolate reductase (DHFR). Their presence was detected by measuring nuclear Overhauser effects between NH protons in protein residues and protons in the individual bound water molecules in two-dimensional nuclear Overhauser effect spectroscopy (NOESY), in nuclear Overhauser effect spectroscopy in the rotating frame (ROESY) and three-dimensional 1H-15N ROESY-heteronuclear multiple quantum coherence spectra recorded on samples containing appropriately 15N-labelled DHFR. For the DHFR-methotrexate-NADPH complex, two bound molecules were found, one close to the Trp5 amide NH proton and the other near to the Trp21 indole HE1 proton: these correspond to two of the water molecules (Wat201 and Wat253) detected in the crystal structure studies described by Bolin and co-workers. However, the nuclear magnetic resonance experiments did not detect any of the other bound water molecules observed in the X-ray studies. The nuclear magnetic resonance results indicate that the two bound water molecules that were detected have lifetimes in the solution state that are longer than approximately two nanoseconds. This is of considerable interest, since one of these water molecules (Wat253) has been implicated as the likely proton donor in the catalytic reduction of dihydrofolate to tetrahydrofolate.

31P-NMR assignment and conformational study of NADPH bound to Lactobacillus casei dihydrofolate reductase based on two-dimensional 1H-31P-heteronuclear and 1H-detected 1H-31P-shift-correlation experiments

For any detailed NMR conformational study of a protein-ligand complex it is essential to have specific resonance assignments. We have now assigned the pyrophosphate 31P resonances in spectra of NADPH bound to Lactobacillus casei dihydrofolate reductase (DHFR) by using a combination of 1H-31P-heteronuclear shift-correlation (HETCOR), 1H-31P-heteronuclear multiple-quantum-coherence correlation spectroscopy (HMQC-COSY), 1H-1H COSY, homonuclear Hartmann-Hahn (HOHAHA) and NOE spectroscopy (NOESY) experiments. The nicotinamide pyrophosphate phosphorus, P(n), has been unequivocally assigned to a signal (-14.07 ppm) which shows a large 3JP-O-C-H coupling constant. Such a coupling constant when combined with the appropriate Karplus relationship provides conformational information about the P-O-C-H torsion angle. The torsion angle changes by 65 degrees +/- 10 degrees for the binary complex compared with the value in free NADPH. The observed coupling constants for the binary (DHFR--NADPH) and ternary (DHFR--NADPH--methotrexate) complexes (12.3 and 10.5 +/- 0.6 Hz, respectively) indicate that the pyrophosphate group has a similar conformation in the two complexes.

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.

Production of mutant dihydrofolate reductases of Lactobacillus casei for nuclear magnetic resonance spectroscopy

Seven mutations (L4P, W21L, D26E, D26N, R57H, R57K and T63Q) affecting residues of dihydrofolate reductase of Lactobacillus casei, suspected of being important in substrate, inhibitor, or cofactor binding, were made by gapped-duplex site-directed mutagenesis. Expression of the L. casei dhfr gene required the removal of nucleotide sequences flanking the coding region. Temperature-inducible expression from the lambda pL promoter of plasmid pPLc28 allowed synthesis and subsequent affinity purification of five mutant proteins in amounts and purity sufficient for nuclear magnetic resonance (NMR) spectroscopic analysis (100 mg or more) from 10-liter cultures. W21L required the growth of 40-liter batches, and L4P was not found. Using a two-plasmid system with pcI857 providing lambda repressor and pMAC5-14 expressing the mutant gene, any auxotrophic strain of Escherichia coli can be used as a host, allowing isotopic labelling of each amino acid of any protein for rapid NMR peak assignment.

Heterologous expression of active thymidylate synthase-dihydrofolate reductase from Plasmodium falciparum

The coding sequence of the bifunctional thymidylate synthase-dihydrofolate reductase (TS-DHFR) from a moderately pyrimethamine-resistant strain (HB3) of Plasmodium falciparum was assembled in a pUC expression vector. The coding sequence possesses unique Nco1 and Xba1 sites which flank 243 bp of the DHFR gene that include all point mutations thus far linked to pyrimethamine resistance. Wild-type (3D7) and highly pyrimethamine-resistant (7G8) TS-DHFRs were made from this vector by cassette mutagenesis using Nco1-Xba1 fragments from the corresponding cloned TS-DHFR genes. Catalytically active recombinant TS-DHFRs were expressed in Escherichia coli, albeit at low levels. Both TS and DHFR coeluted upon gel filtration and copurified upon affinity and anion exchange chromatography. Gel filtration and SDS-PAGE indicated that the enzyme was a dimer with identical 67-kDa subunits, characteristic of protozoan TS-DHFRs. Amino-terminal sequencing gave 10 amino acids which perfectly matched the sequence predicted from the nucleotide sequence. The recombinant TS-DHFR was purified to homogeneity by 10-formylfolate affinity chromatography followed by Mono Q FPLC. The inhibition properties of pyrimethamine toward the purified recombinant enzymes show that the point mutations are the molecular basis of pyrimethamine resistance in P. falciparum.

Thermodynamic study of the interaction of methotrexate, its metabolites, and new antifolates with thymidylate synthase: influence of FdUMP

A microcalorimetric method was used for the direct study of the interaction of methotrexate, its metabolites, and new antifolates N10-propargyl-5,8-dideazafolate (CB 3717) and 2-methyl,2-desamino N10-propargyl-5,8-dideazafolate (CB 3819), with thymidylate synthase. We show that 7-hydroxymethotrexate and dideazafolates require the prior binding of dUMP or its fluorinated derivative FdUMP to bind to thymidylate synthase, as does methotrexate. Conversely, we show that methotrexate-G2 can interact directly with the enzyme alone. On the other hand, both dUMP and FdUMP exhibited a large cooperative effect on the affinity for thymidylate synthase of the inhibitors, and surprisingly, no significant difference was shown at this level between the natural substrate dUMP and its fluorinated derivative. It was demonstrated that this cooperative effect had an enthalpic origin. In the presence of FdUMP or dUMP, all the studied compounds except 7-hydroxymethotrexate exhibited a large negative enthalpy variation when binding to thymidylate synthase (from -44 to -91 kJ/mol). CB 3717 and methotrexate-G2 are competitors for the same protein binding site. Polyglutamation of methotrexate lead to compounds with higher affinity (association constants were 6.6 x 10(3) M-1 and 2.3 x 10(6) M-1 for methotrexate and methotrexate-G2 respectively) while hydroxylation has an unfavourable effect (association constant of 7-hydroxymethotrexate inferior to 500 M-1). Evidence for the influence of polyglutamation was also provided by the relatively low affinity of dideazofolates for thymidylate synthase (association constant equal to 1.4 and 1.7 x 10(7) M-1 for CB 3717 and CB 3819, respectively), whereas these compounds are known to be strong inhibitors of the enzyme in cells in their polyglutamated forms.

A novel method of preparing totally alpha-deuterated amino acids for selective incorporation into proteins. Application to assignment of 1H resonances of valine residues in dihydrofolate reductase

The pyridoxal/2H2O exchange reaction of the alpha-CH of amino acids is known to be accompanied by racemisation: Thus by using a D-amino acid as the starting material any L-amino acid formed in the reaction will be essentially fully deuterated at its alpha-position. We have used this method to prepare alpha-deuterated L-valine and incorporated this biosynthetically into L. casei dihydrofolate reductase. A comparison of the alpha CH-NH fingerprint regions of COSY spectra of deuterated and normal DHFR complexes allows one to identify cross-peaks from 15 of the 16 valine residues.

The 1H-NMR assignments of the aromatic resonances in complexes of Lactobacillus casei dihydrofolate reductase and the origins of their chemical shifts

All the aromatic proton resonances in the 500-MHz NMR spectra of Lactobacillus casei dihydrofolate reductase have been assigned for several of its complexes with inhibitors. For the complexes with methotrexate and trimethoprim this was achieved by using a combination of NMR techniques in conjunction with a selectively deuterated protein designed to simplify the spectra such that nuclear Overhauser effect (NOE) connections could be detected with greater ease and certainty. By correlating these NOE data with crystal structure data on related complexes it was possible to assign all the aromatic resonances and to extend these assignments to spectra of other complexes of dihydrofolate reductase. The conformation-dependent chemical shifts observed for many of the resonances could be explained qualitatively, but not quantitatively, in terms of ring-current shifts. The discrepancies between calculated ring-current shifts and the observed conformation-dependent shifts could not in general be accounted for satisfactorily in terms of carbonyl-group anisotropic shielding contributions calculated using presently available models. In the case of the H delta 1, delta 2 protons of Phe30 some of the discrepancy probably results from a difference in the conformation of the Phe ring between the solution and crystal states.

The purification of dihydrofolate reductase from Drosophila melanogaster

Dihydrofolate reductase (DHFR) has been purified over 30,000-fold from Drosophila adults with a yield of 35%, using a combination of low pH extraction, (NH4)2SO4 precipitation, Sephadex gel filtration, Affi-Gel blue affinity chromatography, ion exchange and gel filtration FPLC. The Drosophila enzyme is a soluble, 17-22 kDa monomeric protein displaying the two pH optima characteristic of eukaryotic DHFRs. The sequence of the first 23 amino acids from the amino-terminal end of the protein shows that Drosophila DHFR is more homologous to the mosquito and vertebrate DHFRs than to the prokaryotic enzymes. However, the percent similarity between the two insect enzymes is not as close as expected when compared to the virtually identical initial sequence conservation of mammalian DHFRs.

Biospecific interactions: their quantitative characterization and use for solute purification

Biospecificity is due largely to the formation and dissociation of non-covalent complexes between small molecules and macromolecules, or between two macromolecules. The first part of this review is concerned with the use of such biospecificity in the fractionation and identification of solutes. Major emphasis is given to affinity chromatography, especially in regard to the practical considerations inherent in an experimental situation and to the wide range of specific interactions that can be utilized. The second part of the review considers the quantitative characterization of biospecific complex formation. The merits of frontal gel chromatography, electrophoretic methods and affinity chromatography are discussed, and special consideration is given to the effects of ligand and/or acceptor multivalency because of its relevance to many biospecific interactions. Finally attention is drawn to the feasibility of employing quantitative affinity chromatographic theory for the determination of association constants for antigen-antibody systems by radioimmunoassay and enzyme-linked immunosorbent assay techniques.

Interactions between inhibitors of dihydrofolate reductase

The binding of substrates and inhibitors to dihydrofolate reductase was studied by steady-state kinetics and high-field 1H-n.m.r. spectroscopy. A series of 5-substituted 2,4-diaminopyrimidines were examined and were found to be 'tightly binding' inhibitors of the enzyme (Ki less than 10(-9) M). Studies on the binding of 4-substituted benzenesulphonamides and benzenesulphonic acids also established the existence of a 'sulphonamide-binding site' on the enzyme. Subsequent n.m.r. experiments showed that there are two binding sites for the sulphonamides on the enzyme, one of which overlaps the coenzyme (NADPH) adenine-ring-binding site. An examination of the pH-dependence of the binding of sulphonamides to the enzyme indicated the influence of an ionizable group on the enzyme that was not directly involved in the sulphonamide binding. The change in pKa value from 6.7 to 7.2 observed on sulphonamide binding suggests the involvement of a histidine residue, which could be histidine-28.

Logarithmic timebase for stopped-flow data acquisition and analysis

A method for capturing stopped-flow and other rapid reaction records in logarithmic time, using a microcomputer, is described. Apart from the ability to record processes over several decades in time in a single experiment, the method shows distinct advantages in subsequent nonlinear regression analysis of multiexponential processes. The method is illustrated by a study of the binding of NADPH to dihydrofolate reductase and the reaction between formycin triphosphate and heavy meromyosin.

Direct 19F NMR observation of the conformational selection of optically active rotamers of the antifolate compound fluoronitropyrimethamine bound to the enzyme dihydrofolate reductase

The molecular basis of the binding of the lipophilic antifolate compound fluoronitropyrimethamine [2,4-diamino-5-(4-fluoro-3-nitrophenyl)-6-ethylpyrimidine] to its target enzyme dihydrofolate reductase has been investigated using a combination of 19F NMR spectroscopy and molecular mechanical calculations. 19F NMR reveals the presence of two different conformational states for the fluoronitropyrimethamine-Lactobacillus casei enzyme complex. MM2 molecular mechanical calculations predict restricted rotation about the C5-C1' bond of the ligand and this gives rise to two slowly interconverting rotamers which are an enantiomeric pair. The results of 19F NMR spectroscopy reveal that both these isomers bind to the enzyme, with different affinities. There is no detectable interconversion of the bound rotamers themselves on the NMR timescale. The effect of the addition of co-enzyme to the sample is to reverse the preference the enzyme has for each rotamer.

Dynamics of trimethoprim bound to dihydrofolate reductase

The conformation of a small molecule in its binding site on a protein is a major factor in the specificity of the interaction between them. In this paper, we report the use of 1H and 13C NMR spectroscopy to study the fluctuations in conformation of the anti-bacterial drug trimethoprim when it is bound to its "target," dihydrofolate reductase. 13C relaxation measurements reveal dihedral angle changes of +/- 25 degrees to +/- 35 degrees on the subnanosecond time scale, while 13C line-shape analysis demonstrates dihedral angle changes of at least +/- 65 degrees on the millisecond time scale. 1H NMR shows that a specific hydrogen bond between the inhibitor and enzyme, which is believed to make an important contribution to binding, makes and breaks rapidly at room temperature.

Direct experimental evidence for competitive inhibition of dihydrofolate reductase by methotrexate

Steady-state enzyme kinetic techniques at very low enzyme concentration (0.4 nM) were used successfully to measure the inhibition constant (53 pM) for the dissociation of methotrexate from the ternary complex of methotrexate, NADPH and dihydrofolate reductase from Lactobacillus casei; and to demonstrate unequivocally that the inhibition was, indeed, competitive with respect to dihydrofolate.

Thermodynamic characterization of the interactions of methotrexate with dihydrofolate reductase by quantitative affinity chromatography

Affinity chromatography on methotrexate-Sepharose has been used to evaluate dissociation constants for interactions of methotrexate with dihydrofolate reductase from Lactobacillus casei. Equilibrium constants of 0.25 microM and 0.6 nM were obtained for dissociation of the inhibitor from the enzyme-methotrexate and enzyme-NADPH-methotrexate complexes, respectively, these estimates being in good agreement with the corresponding published values for dihydrofolate reductase from Streptococcus faecium. By employing a different method for evaluating the thermodynamic dissociation constant for the enzyme-NADPH-methotrexate interaction, this investigation provides independent support for the inference drawn from published fluorescence quenching studies that the interaction of methotrexate with dihydrofolate reductase-NADPH complex is governed by a dissociation constant in the vicinity of 600 pM, a value slightly higher than, but approaching, the inhibition constant of 50-60 pM obtained by enzyme kinetic techniques.