Role of Lys-32 residues in R67 dihydrofolate reductase probed by asymmetric mutations

R67 dihydrofolate reductase (R67 DHFR) is a novel protein encoded by an R-plasmid that confers resistance to the antibiotic, trimethoprim. This homotetrameric enzyme possesses 222 symmetry, which imposes numerous constraints on the single active site pore, including a "one-site-fits-both" strategy for binding its ligands, dihydrofolate (DHF) and NADPH. Previous studies uncovered salt effects on binding and catalysis (Hicks, S. N., Smiley, R. D., Hamilton, J. B., and Howell, E. E. (2003) Biochemistry 42, 10569-10578), however the one or more residues that participate in ionic contacts with the negatively charged tail of DHF as well as the phosphate groups in NADPH were not identified. Several studies predict that Lys-32 residues were involved, however mutations at this residue destabilize the R67 DHFR homotetramer. To study the role of Lys-32 in binding and catalysis, asymmetric K32M mutations have been utilized. To create asymmetry, individual mutations were added to a tandem array of four in-frame gene copies. These studies show one K32M mutation is tolerated quite well, whereas addition of two mutations has variable effects. Two double mutants, K32M:1+2 and K32M: 1+4, which place the mutations on opposite sides of the pore, reduce kcat. However a third double mutant, K32M: 1+3, that places two mutations on the same half pore, enhances kcat 4- to 5-fold compared with the parent enzyme, albeit at the expense of weaker binding of ligands. Because the kcat/Km values for this double mutant series are similar, these mutations appear to have uncovered some degree of non-productive binding. This non-productive binding mode likely arises from formation of an ionic interaction that must be broken to allow access to the transition state. The K32M:1+3 mutant data suggest this interaction is an ionic interaction between Lys-32 and the charged tail of dihydrofolate. This unusual catalytic scenario arises from the 222 symmetry imposed on the single active site pore.     

Hydrophobic interactions via mutants of Escherichia coli dihydrofolate reductase: separation of binding and catalysis

The strictly conserved residue leucine-54 of Escherichia coli dihydrofolate reductase forms part of the hydrophobic wall which binds the p-aminobenzoyl side chain of dihydrofolate. In addition to the previously reported glycine-54 mutant, isoleucine-54 and asparagine-54 substitutions have been constructed and characterized with regard to their effects on binding and catalysis. NADP+ and NADPH binding is virtually unaffected with the exception of a 15-fold decrease in NADPH dissociation from the Gly-54 mutant. The synergistic effect of NADPH on tetrahydrofolate dissociation seen in the wild-type enzyme is lost in the isoleucine-54 mutant: little acceleration is seen in tetrahydrofolate dissociation when cofactor is bound, and there is no discrimination between reduced and oxidized cofactor. The dissociation constants for dihydrofolate and methotrexate increase in the order Leu less than Ile less than Asn less than Gly, varying by a maximum factor of 1700 for dihydrofolate and 6300 for methotrexate. Despite these large changes in binding affinity, the hydride transfer rate of 950 s-1 in the wild-type enzyme is decreased by a constant factor of ca. 30 (2 kcal/mol) regardless of the mutant. Thus, the contributions of residue 54 to binding and catalysis appear to have been separated.     

A balancing act between net uptake of water during dihydrofolate binding and net release of water upon NADPH binding in R67 dihydrofolate reductase

R67 dihydrofolate reductase (DHFR) catalyzes the reduction of dihydrofolate (DHF) to tetrahydrofolate using NADPH as a cofactor. This enzyme is a homotetramer possessing 222 symmetry, and a single active site pore traverses the length of the protein. A promiscuous binding surface can accommodate either DHF or NADPH, thus two nonproductive complexes can form (2NADPH or 2DHF) as well as a productive complex (NADPH.DHF). The role of water in binding was monitored using a number of different osmolytes. From isothermal titration calorimetry (ITC) studies, binding of NADPH is accompanied by the net release of 38 water molecules. In contrast, from both steady state kinetics and ITC studies, binding of DHF is accompanied by the net uptake of water. Although different osmolytes have similar effects on NADPH binding, variable results are observed when DHF binding is probed. Sensitivity to water activity can also be probed by an in vivo selection using the antibacterial drug, trimethoprim, where the water content of the media is decreased by increasing concentrations of sorbitol. The ability of wild type and mutant clones of R67 DHFR to allow host Escherichia coli to grow in the presence of trimethoprim plus added sorbitol parallels the catalytic efficiency of the DHFR clones, indicating water content strongly correlates with the in vivo function of R67 DHFR.     

Structure and hydride transfer mechanism of a moderate thermophilic dihydrofolate reductase from Bacillus stearothermophilus and comparison to its mesophilic and hyperthermophilic homologues

Dihydrofolate reductase (DHFR) from a moderate thermophilic organism, Bacillus stearothermophilus, has been cloned and expressed. Physical characterization of the protein (BsDHFR) indicates that it is a monomeric protein with a molecular mass of 18,694.6 Da (0.8), coincident with the mass of 18 694.67 Da calculated from the primary sequence. Determination of the X-ray structure of BsDHFR provides the first structure for a monomeric DHFR from a thermophilic organism, indicating a high degree of conservation of structure in relation to all chromosomal DHFRs. Structurally based sequence alignment of DHFRs indicates the following levels of sequence identity and similarity for BsDHFR: 38 and 58% with Escherichia coli, 35 and 56% with Lactobacillus casei, and 23 and 40% with Thermotoga maritima, respectively. Steady state kinetic isotope effect studies indicate an ordered kinetic mechanism at elevated temperatures, with NADPH binding first to the enzyme. This converts to a more random mechanism at reduced temperatures, reflected in a greatly reduced K(m) for dihydrofolate at 20 degrees C in relation to that at 60 degrees C. A reduction in either temperature or pH reduces the degree to which the hydride transfer step is rate-determining for the second-order reaction of DHF with the enzyme-NADPH binary complex. Transient state kinetics have been used to study the temperature dependence of the isotope effect on hydride transfer at pH 9 between 10 and 50 degrees C. The data support rate-limiting hydride transfer with a moderate enthalpy of activation (E(a) = 5.5 kcal/mol) and a somewhat greater temperature dependence for the kinetic isotope effect than predicted from classical behavior [A(H)/A(D) = 0.57 (0.15)]. Comparison of kinetic parameters for BsDHFR to published data for DHFR from E. coli and T. maritima shows a decreasing trend in efficiency of hydride transfer with increasing thermophilicity of the protein. These results are discussed in the context of the capacity of each enzyme to optimize H-tunneling from donor (NADPH) to acceptor (DHF) substrates.     

Catalysis by entropic effects: the action of cytidine deaminase on 5,6-dihydrocytidine

In neutral solution, 5,6-dihydrocytidine undergoes spontaneous deamination (k25 approximately 3.2 x 10(-5) s(-1)) much more rapidly than does cytidine (k25 approximately 3.0 x 10(-10) s(-1)), with a more favorable enthalpy of activation (DeltaDeltaH# = -8.7 kcal/mol) compensated by a less favorable entropy of activation (TDeltaDeltaS# = -1.8 kcal/mol at 25 degrees C). E. coli cytidine deaminase enhances the rate of deamination of 5,6-dihydrocytidine (kcat/k(non) = 4.4 x 10(5)) by enhancing the entropy of activation (DeltaDeltaH# = 0 kcal/mol; TDeltaDeltaS# = +7.6 kcal/mol, at 25 degrees C). Binding of the competitive inhibitor 3,4,5,6-tetrahydrouridine (THU), a stable analogue of 5,6-dihydrocytidine in the transition state for its deamination, is accompanied by a release of enthalpy (DeltaH = -7.1 kcal/mol, TDeltaDeltaS = +2.2 kcal/mol) that approaches the estimated enthalpy of binding of the actual substrate in the transition state for deamination of 5,6-dihydrocytidine (DeltaH = -8.1 kcal/mol, TDeltaDeltaS = +6.0 kcal/mol). Thus, the shortcomings of THU in capturing all of the binding affinity expected of an ideal transition-state analogue reflect a less favorable entropy of association. That difference may arise from the analogue's inability to displace a water molecule from the "leaving group site" at which ammonia is generated in the normal reaction. The effect on binding of removing the 4-OH group from the transition-state analogue THU, to form 3,4,5,6-tetrahydrozebularine (THZ) (DeltaDeltaH = -2.1 kcal/mol, TDeltaDeltaS = -4.4 kcal/mol), is mainly entropic, consistent with the inability of THZ to displace water from the "attacking group site". These results are consistent with earlier indications [Snider, M. J., and Wolfenden, R. (2001) Biochemistry 40, 11364] that site-bound water plays a prominent role in substrate activation and inhibitor binding by cytidine deaminase.     

R67, the other dihydrofolate reductase: rational design of an alternate active site configuration

R67 dihydrofolate reductase (DHFR) bears no sequence or structural homologies with chromosomal DHFRs. The gene for this enzyme produces subunits that are 78 amino acids long, which assemble into a homotetramer possessing 222 symmetry. More recently, a tandem array of four gene copies linked in-frame was constructed, which produces a monomer containing 312 amino acids named Quad3. Asymmetric mutations in Quad3 have also been constructed to probe the role of Q67 and K32 residues in catalysis. This present study mixes and matches mutations to determine if the Q67H mutation, which tightens binding approximately 100-fold to both dihydrofolate (DHF) and NADPH, can help rescue the K32M mutation. While the latter mutation weakens DHF binding over 60-fold, it concurrently increases kcat by a factor of 5. Two Q67H mutations were added to gene copies 1 and 4 in conjunction with the K32M mutation in gene copies 1 and 3. Addition of these Q67H mutations tightens binding 40-fold, and the catalytic efficiency (kcat/Km(DHF)) of the resulting protein is similar to that of Quad3. Since these Q67H mutations can mostly compensate for the K32M lesion, K32 must not be necessary for DHF binding. Another multimutant combines the K32M mutation in gene copies 1 and 3 with the Q67H mutation in all gene copies. This mutant is inhibited by DHF but not NADPH, indicating that NADPH binds only to the wild type half of the pore, while DHF can bind to either the wild type or mutant half of the pore. This inhibition pattern contrasts with the mutant containing only the Q67H substitution in all four gene copies, which is severely inhibited by both NADPH and substrate. Since gene duplication and divergence are evolutionary tools for gaining function, these constructs are a first step toward building preferences for NADPH and DHF in each half of the active site pore of this primitive enzyme.     

Interaction of NADP(H) with oxidized and reduced P450 reductase during catalysis. Studies with nucleotide analogues

Previous studies have shown that the interaction of P450 reductase with bound NADP(H) is essential to ensure fast electron transfer through the two flavin cofactors. In this study we investigated in detail the interaction of the house fly flavoprotein with NADP(H) and a number of nucleotide analogues. 1,4,5,6-Tetrahydro-NADP, an analogue of NADPH, was used to characterize the interaction of P450 reductase with the reduced nucleotide. This analogue is inactive as electron donor, but its binding affinity and rate constant of release are very close to those for NADPH. The 2'-phosphate contributes about 5 kcal/mol of the binding energy of NADP(H). Oxidized nicotinamide does not interact with the oxidized flavoprotein, while reduced nicotinamide contributes 1.3 kcal/mol of the binding energy. Oxidized P450 reductase binds NADPH with a K(d) of 0.3 microM, while the affinity of the reduced enzyme is considerably lower, K(d) = 1.9 microM. P450 reductase catalyzes a transhydrogenase reaction between NADPH and oxidized nucleotides, such as thionicotinamide-NADP(+), acetylpyridine-NADP(+), or [(3)H]NADP(+). The reverse reaction, reduction of [(3)H]NADP(+) by the reduced analogues, is also catalyzed by P450 reductase. We define the mechanism of the transhydrogenase reaction as follows: NADPH binding, hydride ion transfer, and release of the NADP(+) formed. An NADP(+) or its analogue binds to the two-electron-reduced flavoprotein, and the electron-transfer steps reverse to transfer hydride ion to the oxidized nucleotide, which is released. Measurements of the flavin semiquinone content, rate constant for NADPH release, and transhydrogenase turnover rates allowed us to estimate the steady-state distribution of P450 reductase species during catalysis, and to calculate equilibrium constants for the interconversion of catalytic intermediates. Our results demonstrate that equilibrium redox potentials of the flavin cofactors are not the sole factor governing rapid electron transfer during catalysis, but conformational changes must be considered to understand P450 reductase catalysis.     

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)     

Effects of distal point-site mutations on the binding and catalysis of dihydrofolate reductase from Escherichia coli

The importance of salt bridge interactions at the NADPH binding site of dihydrofolate reductase has been studied by using site-directed mutagenesis. The mutations R44L and H45Q respectively disrupt the ionic contacts made between the 2'-phosphate and pyrophosphoryl moiety of the coenzyme and the N-terminal region of helix C. Equilibrium fluorescence experiments indicate that while the overall binding of NADPH to both free mutants is weakened by 1.1 and 1.5 kcal/mol (H45Q and R44L, respectively), the binding of dihydrofolate and tetrahydrofolate is unaffected. Despite the similar binding energies for both mutants, the transition state for the chemical hydride step is differentially destabilized relative to wild type (0.6 and 1.8 kcal/mol for H45Q and R44L, respectively). Both stopped-flow and pre-steady-state experiments suggest that the root of this effect may lie in multiple conformations for the E-NADPH complex of R44L. The ability of both mutants to transmit their effects beyond the local environment of the NADPH pocket is manifested in several details: (1) the pKa of Asp-27 (25 A away from the sites of mutation) is elevated from 6.5 in the wild type to 7.5 and 8.4 in H45Q and R44L, respectively; (2) NADPH elevates the off rates for tetrahydrofolate from 12 s-1 in the wild type to greater than 45 s-1 in R44L; and (3) bound tetrahydrofolate decreases the affinity of the enzymes for NADPH as reflected in the Km from 2 to 40 microM for H45Q (similar to wild type) but from 8 to 5000 microM for R44L.     

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.     

Vibrational structure of dihydrofolate bound to R67 dihydrofolate reductase.

R67 is a Type II dihydrofolate reductase (DHFR) that catalyzes the reduction of dihydrofolate (DHF) to tetrahydrofolate by facilitating the addition of a proton to N5 of DHF and the transfer of a hydride ion from NADPH to C6. Because this enzyme is a plasmid-encoded DHFR from trimethoprim-resistant bacteria, extensive studies on R67 with various methods have been performed to elucidate its reaction mechanism. Here, Raman difference measurements, conducted on the ternary complex of R67.NADP(+).DHF believed to be an accurate mimic of the productive DHFR.NADPH.DHF complex, show that the pK(a) of N5 in the complex is less than 4. This is in clear contrast to the behavior observed in Escherichia coli DHFR, a substantially more efficient enzyme, where the pK(a) of bound DHF at N5 is increased to 6.5 compared with its solution value of 2.6. A comparison of the ternary complexes in R67 and E. coli DHFRs suggests that enzymic raising of the pK(a) at N5 can significantly increase the catalytic efficiency of the hydride transfer step. However, R67 shows that even without such a strategy an effective DHFR can still be designed.     

Breaking symmetry: mutations engineered into R67 dihydrofolate reductase,a D2 symmetric homotetramer possessing a single active site pore

R67 dihydrofolate reductase (DHFR) is an enzyme, encoded by an R-plasmid, that confers resistance to the antibacterial agent trimethoprim. This homotetramer possesses a single active site pore and exact 222 symmetry. The symmetry imposes constraints on the ability of the enzyme to optimize binding of the substrate, dihydrofolate (DHF), and the cofactor, NADPH, resulting in a "one site fits both ligands" approach. This approach allows formation of either a NADPH.NADPH, dihydrofolate.dihydrofolate, or NADPH.dihydrofolate complex. The first two complexes are nonproductive, while the third is the productive catalytic species. To break the symmetry of the active site, a tandem array of four R67 DHFR genes has been linked in frame, allowing individual manipulation of each gene copy. Various numbers and combinations of asymmetric Q67H mutations have been engineered into the tandem gene array. The Q67H mutation was chosen for investigation as it was previously found to tighten binding to both dihydrofolate and NADPH by approximately 100-fold in homotetrameric R67 DHFR [Park, H., Bradrick, T. D., and Howell, E. E. (1997) Protein Eng. 10, 1415-1424]. Nonadditive effects on ligand binding are observed when one to four mutations are inserted, indicating either conformational changes in the protein or different cooperativity patterns in the ligand-ligand interactions. From steady state kinetics, addition of Q67H mutations does not drastically affect formation of the NADPH.dihydrofolate complex; however, a large energy difference between the productive and nonproductive complexes is no longer maintained. A role for Q67 in discriminating between these various states is proposed. Since theories of protein evolution suggest gene duplication followed by accumulation of mutations can lead to divergence of activity, this study is a first step toward asking if introduction of asymmetric mutations in the quadrupled R67 DHFR gene can lead to optimization of ligand binding sites.     

Novel mechanism-based substrates of dihydrofolate reductase and the thermodynamics of ligand binding: A comparison of theory and experiment for 8-methylpterin and 6,8-dimethylpterin

Molecular dynamics simulation and free energy perturbation techniques have been used to study the relative binding free energies of the designed mechanism-based pterins, 8-methylpterin and 6,8-dimethylpterin, to dihydrofolate reductase (DHFR), with co-factor nicotinamide adenine dinucleotide phosphate (NADPH). The calculated free energy differences suggest that DHFR.NADPH.6,8-dimethylpterin is thermodynamically more stable than DHFR.NADPH.8-methylpterin by 2.4 kcal/mol when the substrates are protonated and by 1.3 kcal/mol when neutral. The greater binding strength of 6,8-dimethylpterin may be attributed largely to hydration effects. In terms of an appropriate model for the pH-dependent kinetic mechanism, these differences can be interpreted consistently with experimental data obtained from previous kinetic studies, i.e., 6,8-dimethylpterin is a more efficient substrate of vertebrate DHFRs than 8-methylpterin. The kinetic data suggest a value of 6.6 ± 0.2 for the pK_a of the active site Glu-30 in DHFR.NADPH. We have also used experimental data to estimate absolute values for thermodynamic dissociation constants of the active (i.e., protonated) forms of the substrates: these are of the same order as for the binding of folate (0.1–10 μM). The relative binding free energy calculated from the empirically derived dissociation constants for the protonated forms of 8-methylpterin and 6,8-dimethylpterin is 1.4 kcal/mol, a value which compares reasonably well with the theoretical value of 2.4 kcal/mol. © 1993 Wiley-Liss, Inc.     

Analysis of hydride transfer and cofactor fluorescence decay in mutants of dihydrofolate reductase: possible evidence for participation of enzyme molecular motions in catalysis.

A remarkable correlation has been discovered between fluorescence lifetimes of bound NADPH and rates of hydride transfer among mutants of dihydrofolate reductase (DHFR) from Escherichia coli. Rates of hydride transfer from NADPH to dihydrofolate change by a factor of 1,000 for the series of mutant enzymes. Since binding constants for the initial complex between coenzyme and DHFR change by only a factor of 10, the major portion of the change in hydride transfer must be attributed to losses in transition-state stabilization. The time course of fluorescence decay for NADPH bound to DHFR is biphasic. Lifetimes ranging from 0.3 to 0.5 ns are attributed to a solvent-exposed dihydronicotinamide conformation of bound coenzyme which is presumably not active in catalysis, while decay times (tau 2) in the range of 1.3 to 2.3 ns are assigned to a more tightly bound species of NADPH in which dihydronicotinamide is sequestered from solvent. It is this slower component that is of interest. Ternary complexes with three different inhibitors, methotrexate, 5-deazafolate, and trimethoprim, were investigated, along with the holoenzyme complex; 3-acetylNADPH was also investigated. Fluorescence polarization decay, excitation polarization spectra, the temperature variation of fluorescence lifetimes, fluorescence amplitudes, and wavelength of absorbance maxima were measured. We suggest that dynamic quenching or internal conversion promotes decay of the excited state in NADPH-DHFR. When rates of hydride transfer are plotted against the fluorescence lifetime (tau 2) of tightly bound NADPH, an unusual correlation is observed. The fluorescence lifetime becomes longer as the rate of catalysis decreases for most mutants studied. However, the fluorescence lifetime is unchanged for those mutations that principally alter the binding of dihydrofolate while leaving most dihydronicotinamide interactions relatively undisturbed. The data are interpreted in terms of possible dynamic motions of a flexible loop region in DHFR which closes over both substrate and coenzyme binding sites. These motions could lead to faster rates of fluorescence decay in holoenzyme complexes and, when correlated over time, may be involved in other motions which give rise to enhanced rates of catalysis in DHFR.     

Selective probing of a NADPH site controlled light‐induced enzymatic catalysis

Achieving molecular recognition of NADPH binding sites is a compelling strategy to control many redox biological processes. The NADPH sites recognize the ubiquitous NADPH cofactor via highly conserved binding interactions, despite differences in the regulation of the hydride transfer in redox active proteins. We recently developed a photoactive NADPH substitute, called nanotrigger NT synchronizing the initiation of enzymatic catalysis of the endothelial NO‐synthase (eNOS) with a laser pulse. Spatial and temporal control of enzymatic activity by such a designed light‐driven activator would benefit from achieving molecular selectivity, i.e. activation of a single NADPH‐mediated enzyme. In this work, we probe the ability of NT to discriminate between two NADPH sites with light. The selected NADPH sites belong to dihydrofolate reductase dihydrofolate reductase enzyme (DHFR) and endothelial NO‐synthase (eNOS). Ultrafast kinetics showed that NT could not activate DHFR catalysis with a laser pulse in contrast with the observed trigger of eNOS catalysis leading to NO formation. Homology modelling, molecular dynamics simulations showed that NT discriminated between the two NADPH sites by different donor to acceptor distances and by local steric effects hindering light activation of DHFR catalysis. The data suggested that the narrow NADPH site required a tight fit of the nanotrigger at a suitable distance/angle to the electron acceptor for a specific activation of the catalysis. The ability of the nanotrigger to activate eNOS combined with a low reactivity in unfavourable NADPH sites makes NT a highly promising tool for targeting eNOS in endothelial cells with a laser pulse. Copyright © 2009 John Wiley & Sons, Ltd.     

"Catch 222," the effects of symmetry on ligand binding and catalysis in R67 dihydrofolate reductase as determined by mutations at Tyr-69

R67 dihydrofolate reductase (R67 DHFR) catalyzes the transfer of a hydride ion from NADPH to dihydrofolate, generating tetrahydrofolate. The homotetrameric enzyme provides a unique environment for catalysis as both ligands bind within a single active site pore possessing 222 symmetry. Mutation of one active site residue results in concurrent mutation of three additional symmetry-related residues, causing large effects on binding of both ligands as well as catalysis. For example, mutation of symmetry-related tyrosine 69 residues to phenylalanine (Y69F), results in large increases in Km values for both ligands and a 2-fold rise in the kcat value for the reaction (Strader, M. B., Smiley, R. D., Stinnett, L. G., VerBerkmoes, N. C., and Howell, E. E. (2001) Biochemistry 40, 11344-11352). To understand the interactions between specific Tyr-69 residues and each ligand, asymmetric Y69F mutants were generated that contain one to four Y69F mutations. A general trend observed from isothermal titration calorimetry and steady-state kinetic studies of these asymmetric mutants is that increasing the number of Y69F mutations results in an increase in the Kd and Km values. In addition, a comparison of steady-state kinetic values suggests that two Tyr-69 residues in one half of the active site pore are necessary for NADPH to exhibit a wild-type Km value. A tyrosine 69 to leucine mutant was also generated to approach the type(s) of interaction(s) occurring between Tyr-69 residues and the ligands. These studies suggest that the hydroxyl group of Tyr-69 is important for interactions with NADPH, whereas both the hydroxyl group and hydrophobic ring atoms of the Tyr-69 residues are necessary for proper interactions with dihydrofolate.     

Role of ionic interactions in ligand binding and catalysis of R67 dihydrofolate reductase

R67 dihydrofolate reductase (DHFR), which catalyzes the NADPH dependent reduction of dihydrofolate to tetrahydrofolate, belongs to a type II family of R-plasmid encoded DHFRs that confer resistance to the antibacterial drug trimethoprim. Crystal structure data reveals this enzyme is a homotetramer that possesses a single active site pore. Only two charged residues in each monomer are located near the pore, K32 and K33. Site-directed mutants were constructed to probe the role of these residues in ligand binding and/or catalysis. As a result of the 222 symmetry of this enzyme, mutagenesis of one residue results in modification at four related sites. All mutants at K32 affected the quaternary structure, producing an inactive dimer. The K33M mutant shows only a 2-4-fold effect on K(m) values. Salt effects on ligand binding and catalysis for K33M and wildtype R67 DHFRs were investigated to determine if these lysines are involved in forming ionic interactions with the negatively charged substrates, dihydrofolate (overall charge of -2) and NADPH (overall charge of -3). Binding studies indicate that two ionic interactions occur between NADPH and R67 DHFR. In contrast, the binding of folate, a poor substrate, to R67 DHFR.NADPH appears weak as a titration in enthalpy is lost at low ionic strength. Steady-state kinetic studies for both wild type (wt) and K33M R67 DHFRs also support a strong electrostatic interaction between NADPH and the enzyme. Interestingly, quantitation of the observed salt effects by measuring the slopes of the log of ionic strength versus the log of k(cat)/K(m) plots indicates that only one ionic interaction is involved in forming the transition state. These data support a model where two ionic interactions are formed between NADPH and symmetry related K32 residues in the ground state. To reach the transition state, an ionic interaction between K32 and the pyrophosphate bridge is broken. This unusual scenario likely arises from the constraints imposed by the 222 symmetry of the enzyme.     

Protein-Ligand interaction studies on 2, 4, 6- trisubstituted triazine derivatives as anti-malarial DHFR agents using AutoDock

The dihydrofolate reductase (DHFR) domain of P. falciparum is one of the few well defined targets in malarial chemotherapy. The enzyme catalyzes the nicotinamide adenine dinucleotide phosphate (NADPH) dependent reduction of dihydrofolate to tetrahydrofolate. Protein-ligand interactions were studied using DHFR protein 2BL9, extracted from PDB to evaluate the strength of affinity of various molecules towards ligand binding site and to study the extent of correlation between experimental values and computational dock scores. AutoDock runs resulted in binding energy scores from -7.14 to -10.72 kcal/mol. Among the five inhibitors (Bioorganic and Medicinal Chemistry Letters 15 2005 531-533) selected for docking studies, an excellent correlation was observed in all cases, for instance, experimentally reported most active molecule 2a (MIC: 1μg/ml) showed a high dock score (-10.72 kcal/mol) than the remaining inhibitors. Therefore, molecular docking using AutoDock suggests the importance of evaluating the prediction accuracy of various molecules as evidenced by a correlation coefficient of 0.961 between experimental activities and AutoDock binding energies.     

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

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

Insight into the molecular mechanism about lowered dihydrofolate binding affinity to dihydrofolate reductase-like 1 (DHFRL1)

Human dihydrofolate reductase-like 1 (DHFRL1) has been identified as a second human dihydrofolate reductase (DHFR) enzyme. Although DHFRL1 have high sequence homology with human DHFR, dihydrofolate (DHF) exhibits a lowered binding affinity to DHFRL1 and the corresponding molecular mechanism is still unknown. To address this question, we studied the binding of DHF to DHFRL1 and DHFR by using molecular dynamics simulation. Moreover, to investigate the role the 24th residue of DHFR/DHFRL1 plays in DHF binding, R24W DHFRL1 mutant was also studied. The van der Waals interaction are more crucial for the total DHF binding energies, while the difference between the DHF binding energies of human DHFR and DHFRL1 can be attributed to the electrostatic interaction and the polar desolvation free energy. More specifically, lower DHF affinity to DHFRL1 can be mainly attributed to the reduction of net electrostatic interactions of residues Arg32 and Gln35 of DHFRL1 with DHF as being affected by Arg24. The side chain of Arg24 in DHFRL1 can extend deeply into the binding sites of DHF and NADPH, and disturb the DHF binding by steric effect, which rarely happens in human DHFR and R24W DHFRL1 mutant. Additionally, the conformation of loop I in DHFRL1 was also studied in this work. Interestingly, the loop conformation resemble to normal closed state of Escherichia coli DHFR other than the closed state of human DHFR. We hope this work will be useful to understand the general characteristics of DHFRL1.