Chiral evasion and stereospecific antifolate resistance in Staphylococcus aureus

Antimicrobial resistance presents a significant health care crisis. The mutation F98Y in Staphylococcus aureus dihydrofolate reductase (SaDHFR) confers resistance to the clinically important antifolate trimethoprim (TMP). Propargyl-linked antifolates (PLAs), next generation DHFR inhibitors, are much more resilient than TMP against this F98Y variant, yet this F98Y substitution still reduces efficacy of these agents. Surprisingly, differences in the enantiomeric configuration at the stereogenic center of PLAs influence the isomeric state of the NADPH cofactor. To understand the molecular basis of F98Y-mediated resistance and how PLAs’ inhibition drives NADPH isomeric states, we used protein design algorithms in the osprey protein design software suite to analyze a comprehensive suite of structural, biophysical, biochemical, and computational data. Here, we present a model showing how F98Y SaDHFR exploits a different anomeric configuration of NADPH to evade certain PLAs’ inhibition, while other PLAs remain unaffected by this resistance mechanism.     

Structure-based design of new DHFR-based antibacterial agents: 7-aryl-2,4-diaminoquinazolines

Dihydrofolate reductase (DHFR) inhibitors such as trimethoprim (TMP) have long played a significant role in the treatment of bacterial infections. Not surprisingly, after decades of use there is now bacterial resistance to TMP and therefore a need to develop novel antibacterial agents with expanded spectrum including these resistant strains. In this study, we investigated the optimization of 2,4-diamnoquinazolines for antibacterial potency and selectivity. Using structure-based drug design, several 7-aryl-2,4-diaminoquinazolines were discovered that have excellent sub-100 picomolar potency against bacterial DHFR. These compounds have good antibacterial activity especially on gram-positive pathogens including TMP-resistant strains.     

Crystallographic analysis reveals a novel second binding site for trimethoprim in active site double mutants of human dihydrofolate reductase

In order to produce a more potent replacement for trimethoprim (TMP) used as a therapy for Pneumocystis pneumonia and targets dihydrofolate reductase from Pneumocystis jirovecii (pjDHFR), it is necessary to understand the determinants of potency and selectivity against DHFR from the mammalian host and fungal pathogen cells. To this end, active site residues in human (h) DHFR were replaced with those from pjDHFR. Structural data are reported for two complexes of TMP with the double mutants Gln35Ser/Asn64Phe (Q35S/N64F) and Gln35Lys/Asn64Phe (Q35K/N64F) of hDHFR that unexpectedly show evidence for the binding of two molecules of TMP: one molecule that binds in the normal folate binding site and the second molecule that binds in a novel subpocket site such that the mutated residue Phe64 is involved in van der Waals contacts to the trimethoxyphenyl ring of the second TMP molecule. Kinetic data for the binding of TMP to hDHFR and pjDHFR reveal an 84-fold selectivity of TMP against pjDHFR (K_i 49 nM) compared to hDHFR (K_i 4093 nM). Two mutants that contain one substitution from pj- and one from the closely related Pneumocystis carinii DHFR (pcDHFR) (Q35K/N64F and Q35S/N64F) show K_i values of 593 and 617 nM, respectively; these K_i values are well above both the K_i for pjDHFR and are similar to pcDHFR (Q35K/N64F and Q35S/N64F) (305 nM). These results suggest that active site residues 35 and 64 play key roles in determining selectivity for pneumocystis DHFR, but that other residues contribute to the unique binding of inhibitors to these enzymes.     

Pressure dependence of activity and stability of dihydrofolate reductases of the deep-sea bacterium Moritella profunda and Escherichia coli

To understand the pressure-adaptation mechanism of deep-sea enzymes, we studied the effects of pressure on the enzyme activity and structural stability of dihydrofolate reductase (DHFR) of the deep-sea bacterium Moritella profunda (mpDHFR) in comparison with those of Escherichia coli (ecDHFR). mpDHFR exhibited optimal enzyme activity at 50 MPa whereas ecDHFR was monotonically inactivated by pressure, suggesting inherent pressure-adaptation mechanisms in mpDHFR. The secondary structure of apo-mpDHFR was stable up to 80 °C, as revealed by circular dichroism spectra. The free energy changes due to pressure and urea unfolding of apo-mpDHFR, determined by fluorescence spectroscopy, were smaller than those of ecDHFR, indicating the unstable structure of mpDHFR against pressure and urea despite the three-dimensional crystal structures of both DHFRs being almost the same. The respective volume changes due to pressure and urea unfolding were − 45 and − 53 ml/mol at 25 °C for mpDHFR, which were smaller (less negative) than the corresponding values of − 77 and − 85 ml/mol for ecDHFR. These volume changes can be ascribed to the difference in internal cavity and surface hydration of each DHFR. From these results, we assume that the native structure of mpDHFR is loosely packed and highly hydrated compared with that of ecDHFR in solution.     

Diagnostic chemical shift markers for loop conformation and substrate and cofactor binding in dihydrofolate reductase complexes

Heteronuclear NMR methods have been used to probe the conformation of four complexes of Escherichia coli dihydrofolate reductase (DHFR) in solution. (1)H(N), (15)N, and (13)C(alpha) resonance assignments have been made for the ternary complex with folate and oxidized NADP(+) cofactor and the ternary complex with folate and a reduced cofactor analog, 5,6-dihydroNADPH. The backbone chemical shifts have been compared with those of the binary complex of DHFR with the substrate analog folate and the binary complex with NADPH (the holoenzyme). Analysis of (1)H(N) and (15)N chemical shifts has led to the identification of marker resonances that report on the active site conformation of the enzyme. Other backbone amide resonances report on the presence of ligands in the pterin binding pocket and in the adenosine and nicotinamide-ribose binding sites of the NADPH cofactor. The chemical shift data indicate that the enzyme populates two dominant structural states in solution, with the active site loops in either the closed or occluded conformations defined by X-ray crystallography; there is no evidence that the open conformation observed in some X-ray structures of E. coli DHFR are populated in solution.     

In vitro trials of some antimicrobial combinations against Staphylococcus aureus and Pseudomonas aeruginosa

Staphylococcus aureus and Pseudomonas aeruginosa are rapidly increasing as multidrug resistant strains worldwide. In nosocomial settings because of heavy exposure of different antimicrobials, resistance in these pathogens turned into a grave issue in both developed and developing countries. The aim of this study was to investigate in vitro antibiotic synergism of combinations of β-lactam–β-lactam and β-lactam–aminoglycoside against clinical isolates of S. aureus and P. aeruginosa. Synergy was determined by checkerboard double dilution method. The combination of amoxicillin and cefadroxil was found to be synergistic against 47 S. aureus isolates, in the FICI range of 0.14–0.50 (81.03%) followed by the combination of streptomycin and cefadroxil synergistic against 44 S. aureus isolates in the FICI range of 0.03–0.50 (75.86%). The combination of streptomycin and cefadroxil was observed to be synergistic against 39 P. aeruginosa isolates in the FICI range of 0.16–0.50 (81.28%). Further actions are needed to characterize the possible interaction mechanism between these antibiotics. Moreover, the combination of streptomycin and cefadroxil may lead to the development of a new and vital antimicrobial against simultaneous infections of S. aureus and P. aeruginosa.     

The role of tryptophan-48 in catalysis and binding of inhibitors of Plasmodium falciparum dihydrofolate reductase

Dihydrofolate reductases (DHFRs) from Plasmodium falciparum (Pf) and various species of both prokaryotic and eukaryotic organisms have a conserved tryptophan (Trp) at position 48 in the active site. The role in catalysis and binding of inhibitors of the conserved Trp48 of PfDHFR has been analysed by site-specific mutagenesis, enzyme kinetics and use of a bacterial surrogate system. All 19 mutant enzymes showed undetectable or very low specific activities, with the highest value of k(cat)/K(m) from the Tyr48 (W48Y) mutant (0.12 versus 11.94M(-1)s(-1)), of about 1% of the wild-type enzyme. The inhibition constants for pyrimethamine, cycloguanil and WR99210 of the W48Y mutants are 2.5-5.3 times those of the wild-type enzyme. All mutants, except W48Y, failed to support the growth of Escherichia coli transformed with the parasite gene in the presence of trimethoprim, indicating the loss of functional activity of the parasite enzyme. Hence, Trp48 plays a crucial role in catalysis and inhibitor binding of PfDHFR. Interestingly, W48Y with an additional mutation at Asn188Tyr (N188Y) was found to promote bacterial growth and yielded a higher amount of purified enzyme. However, the kinetic parameters of the purified W48Y+N188Y enzyme were comparable with W48Y and the binding affinities for DHFR inhibitors were also similar to the wild-type enzyme. Due to its conserved nature, Trp48 of PfDHFR is a potential site for interaction with antimalarial inhibitors which would not be compromised by its mutations.     

Cloning,recombinant expression and inhibitor profiles of dihydrofolate reductase from the Australian sheep blow fly,Lucilia cuprina

While dihydrofolate reductase (DHFR) is an important drug target in mammals, bacteria and protozoa, no inhibitors of this enzyme have been developed as commercial insecticides. We therefore examined the potential of this enzyme as a drug target in an important ectoparasite of livestock, the Australian sheep blow fly, Lucilia cuprina (Diptera: Calliphoridae) (Wiedemann). The non‐specific DHFR inhibitors aminopterin and methotrexate significantly inhibited the growth of L. cuprina larvae, with IC₅₀ values at µg levels. Trimethoprim and pyrimethamine were 5–30‐fold less active. Relative IC₅₀ values for the inhibition of recombinant L. cuprina DHFR by various inhibitors were in accordance with their relative effects on larval growth. The active‐site amino acid residues of L. cuprina DHFR differed by between 34% and 50% when compared with two mammalian species, as well as two bacteria and two protozoa. There were significant charge and size differences in specific residues between the blow fly and human DHFR enzymes, notably the L. cuprina Asn21, Lys31 and Lys63 residues. This study provides bioassay evidence to highlight the potential of blow fly DHFR as an insecticide target, and describes differences in active site residues between blow flies and other organisms which could be exploited in the design of blow fly control chemicals.     

Structures of S. aureus thymidylate kinase reveal an atypical active site configuration and an intermediate conformational state upon substrate binding

Methicillin-resistant Staphylococcus aureus (MRSA) poses a major threat to human health, particularly through hospital acquired infection. The spread of MRSA means that novel targets are required to develop potential inhibitors to combat infections caused by such drug-resistant bacteria. Thymidylate kinase (TMK) is attractive as an antibacterial target as it is essential for providing components for DNA synthesis. Here, we report crystal structures of unliganded and thymidylate-bound forms of S. aureus thymidylate kinase (SaTMK). His-tagged and untagged SaTMK crystallize with differing lattice packing and show variations in conformational states for unliganded and thymidylate (TMP) bound forms. In addition to open and closed forms of SaTMK, an intermediate conformation in TMP binding is observed, in which the site is partially closed. Analysis of these structures indicates a sequence of events upon TMP binding, with helix alpha3 shifting position initially, followed by movement of alpha2 to close the substrate site. In addition, we observe significant conformational differences in the TMP-binding site in SaTMK as compared to available TMK structures from other bacterial species, Escherichia coli and Mycobacterium tuberculosis as well as human TMK. In SaTMK, Arg 48 is situated at the base of the TMP-binding site, close to the thymine ring, whereas a cis-proline occupies the equivalent position in other TMKs. The observed TMK structural differences mean that design of compounds highly specific for the S. aureus enzyme looks possible; such inhibitors could minimize the transfer of drug resistance between different bacterial species.     

Structural and functional studies of fatty acyl adenylate ligases from E. coli and L. pneumophila

Fatty acyl-AMP ligase (FAAL) is a new member of a family of adenylate-forming enzymes that were recently discovered in Mycobacterium tuberculosis. They are similar in sequence to fatty acyl-coenzyme A (CoA) ligases (FACLs). However, while FACLs perform a two-step catalytic reaction, AMP ligation followed by CoA ligation using ATP and CoA as cofactors, FAALs produce only the acyl adenylate and are unable to perform the second step. We report X-ray crystal structures of full-length FAAL from Escherichia coli (EcFAAL) and FAAL from Legionella pneumophila (LpFAAL) bound to acyl adenylate, determined at resolution limits of 3.0 and 1.85 Å, respectively. The structures share a larger N-terminal domain and a smaller C-terminal domain, which together resemble the previously determined structures of FAAL and FACL proteins. Our two structures occur in quite different conformations. EcFAAL adopts the adenylate-forming conformation typical of FACLs, whereas LpFAAL exhibits a unique intermediate conformation. Both EcFAAL and LpFAAL have insertion motifs that distinguish them from the FACLs. Structures of EcFAAL and LpFAAL reveal detailed interactions between this insertion motif and the interdomain hinge region and with the C-terminal domain. We suggest that the insertion motifs support sufficient interdomain motions to allow substrate binding and product release during acyl adenylate formation, but they preclude CoA binding, thereby preventing CoA ligation.     

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.     

Crystal Structure of Glyceraldehyde-3-Phosphate Dehydrogenase 1 from Methicillin-Resistant Staphylococcus aureus MRSA252 Provides Novel Insights into Substrate Binding and Catalytic Mechanism

The dreaded pathogen Staphylococcus aureus is one of the causes of morbidity and mortality worldwide. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), one of the key glycolytic enzymes, is irreversibly oxidized under oxidative stress and is responsible for sustenance of the pathogen inside the host. With an aim to elucidate the catalytic mechanism and identification of intermediates involved, we describe in this study different crystal structures of GAPDH1 from methicillin-resistant S. aureus MRSA252 (SaGAPDH1) in apo and holo forms of wild type, thioacyl intermediate, and ternary complexes of active-site mutants with physiological substrate d-glyceraldehyde-3-phosphate (G3P) and coenzyme NAD⁺. A new phosphate recognition site, “new P_i” site, similar to that observed in GAPDH from Thermotoga maritima, is reported here, which is 3.40 Å away from the “classical P_i” site. Ternary complexes discussed are representatives of noncovalent Michaelis complexes in the ground state. d-G3P is bound to all the four subunits of C151S.NAD and C151G.NAD in more reactive hydrate (gem-di-ol) form. However, in C151S + H178N.NAD, the substrate is bound to two chains in aldehyde form and in gem-di-ol form to the other two. This work reports binding of d-G3P to the C151G mutant in an inverted manner for the very first time. The structure of the thiaocyl complex presented here is formed after the hydride transfer. The C3 phosphate of d-G3P is positioned at the “P_s” site in the ternary complexes but at the “new P_i” site in the thioacyl complex and C1-O1 bond points opposite to His178 disrupting the alignment between itself and NE2 of His178. A new conformation (Conformation I) of the 209-215 loop has also been identified, where the interaction between phosphate ion at the “new P_i” site and conserved Gly212 is lost. Altogether, inferences drawn from the kinetic analyses and crystal structures suggest the “flip-flop” model proposed for the enzyme mechanism.     

Influence of environment on the antifolate drug trimethoprim: energy minimization studies

Environmental effects on trimethoprim (TMP), an inhibitor of bacterial dihydrofolate reductase (DHFR), were investigated with energy minimizations in vacuo, in the crystal, and in aqueous solution. The conformations, harmonic dynamics, and energetics of the antibacterial drug calculated in these environments were compared with each other and with those of two enzyme-bound drugs. Valence and torsion angles and their energies and overall intra- and intermolecular energies compensated one another in the minimized TMP structures. The conformations of the isolated and aqueous molecules were similar to that of TMP bound to chicken liver DHFR, while the structures from the TMP crystal and from the Escherichia coli DHFR complex were unique. Since neither the small-molecule crystal nor a local minimum of the isolated molecule gave the conformation of TMP bound to the bacterial enzyme, a combination of several experimental and theoretical techniques may be necessary to probe accessible conformations of a molecule.     

Crystal structure of the Clostridium limosum C3 exoenzyme

C3-like toxins ADP-ribosylate and inactivate Rho GTPases. Seven C3-like ADP-ribosyltransferases produced by Clostridium botulinum, Clostridium limosum, Bacillus cereus and Staphylococcus aureus were identified and two representatives - C3bot from C. botulinum and C3stau2 from S. aureus - were crystallized. Here we present the 1.8 Å structure of C. limosum C3 transferase C3lim and compare it to the structures of other family members. In contrast to the structure of apo-C3bot, the canonical ADP-ribosylating turn turn motif is observed in a primed conformation, ready for NAD binding. This suggests an impact on the binding mode of NAD and on the transferase reaction. The crystal structure explains why auto-ADP-ribosylation of C3lim at Arg41 interferes with the ADP-ribosyltransferase activity of the toxin.     

Kinetic folding of Haloferax volcanii and Escherichia coli dihydrofolate reductases: haloadaptation by unfolded state destabilization at high ionic strength

Salts affect protein stability by multiple mechanisms (e.g., the Hofmeister effect, preferential hydration, electrostatic effects and weak ion binding). These mechanisms can affect the stability of both the native state and the unfolded state. Previous equilibrium stability studies demonstrated that KCl stabilizes dihydrofolate reductases (DHFRs) from Escherichia coli (ecDHFR, E. coli DHFR) and Haloferax volcanii (hvDHFR1, H. volcanii DHFR encoded by the hdrA gene) with similar efficacies, despite adaptation to disparate physiological ionic strengths (0.2 M versus 2 M). Kinetic studies can provide insights on whether equilibrium effects reflect native state stabilization or unfolded state destabilization. Similar kinetic mechanisms describe the folding of urea-denatured ecDHFR and hvDHFR1: a 5-ms stopped-flow burst-phase species that folds to the native state through two sequential intermediates with relaxation times of 0.1-3 s and 25-100 s. The latter kinetic step is very similar to that observed for the refolding of hvDHFR1 from low ionic strength. The unfolding of hvDHFR1 at low ionic strength is relatively slow, suggesting kinetic stabilization as observed for some thermophilic enzymes. Increased KCl concentrations slow the urea-induced unfolding of ecDHFR and hvDHFR1, but much less than expected from equilibrium studies. Unfolding rates extrapolated to 0 M denaturant, k_unf(H₂O), are relatively independent of ionic strength, demonstrating that the KCl-induced stabilization of ecDHFR and hvDHFR1 results predominantly from destabilization of the unfolded state. This supports the hypothesis from previous equilibrium studies that haloadaptation harnesses the effects of elevated salt concentrations on the properties of the aqueous solvent to enhance protein stability.     

Targeting species specific amino acid residues: Design,synthesis and biological evaluation of 6-substituted pyrrolo[2,3-d]pyrimidines as dihydrofolate reductase inhibitors and potential anti-opportunistic infection agents

To combine the potency of trimetrexate (TMQ) or piritrexim (PTX) with the species selectivity of trimethoprim (TMP), target based design was carried out with the X-ray crystal structure of human dihydrofolate reductase (hDHFR) and the homology model of Pneumocystis jirovecii DHFR (pjDHFR). Using variation of amino acids such as Met33/Phe31 (in pjDHFR/hDHFR) that affect the binding of inhibitors due to their distinct positive or negative steric effect at the active binding site of the inhibitor, we designed a series of substituted-pyrrolo[2,3-d]pyrimidines. The best analogs displayed better potency (IC₅₀) than PTX and high selectivity for pjDHFR versus hDHFR, with 4 exhibiting a selectivity for pjDHFR of 24-fold.     

Identification of potential targets in Staphylococcus aureus N315 using computer aided protein data analysis

Staphylococcus aureus is a gram positive bacterium, responsible for both community-acquired and hospital-acquired infection, resulting in a mortality rate of 39%. 43.2% resistance to methicilin and emerging resistance to Fluroquinolone and Oxazolidinone, have evoked the necessity of the establishment of alternative and effective therapeutic approach to treat this bacteria. In this computational study, various database and online software are used to determine some specific targets of Staphylococcus aureus N315 other than those used by Penicillin, Quinolone and Oxazolidinone. For this purpose, among 302 essential proteins, 101 nonhomologous proteins were accrued and 64 proteins which are unique in several metabolic pathways of S. aureus were isolated by using metabolic pathway analysis tools. Furthermore, 7 essentially unique enzymes involved in exclusive metabolic pathways were revealed by this research, which can be potential drug target. Along with these important enzymes, 15 non-homologous proteins located on membrane were identified, which can play a vital role as potential therapeutic targets for the future researchers.     

Mechanistic analysis of allosteric and non-allosteric effects arising from nanobody binding to two epitopes of the dihydrofolate reductase of Escherichia coli

Although allosteric effector antibodies are used widely as modulators of receptors and enzymes, experimental analysis of their mechanism remains highly challenging. Here, we investigate the molecular mechanisms of allosteric and non-allosteric effector antibodies in an experimentally tractable system, consisting of single-domain antibodies (nanobodies) that target the model enzyme dihydrofolate reductase (DHFR) from Escherichia coli. A panel of thirty-five nanobodies was isolated using several strategies to increase nanobody diversity. The nanobodies exhibit a variety of effector properties, including partial inhibition, strong inhibition and stimulation of DHFR activity. Despite these diverse effector properties, chemical shift perturbation NMR epitope mapping identified only two epitope regions: epitope α is a new allosteric site that is over 10 Å from the active site, while epitope β is located in the region of the Met20 loop. The structural basis for DHFR allosteric inhibition or activation upon nanobody binding to the α epitope was examined by solving the crystal structures of DHFR in complex with Nb113 (an allosteric inhibitor) and Nb179 (an allosteric activator). The structures suggest roles for conformational constraint and altered protein dynamics, but not epitope distortion, in the observed allosteric effects. The crystal structure of a β epitope region binder (ca1698) in complex with DHFR is also reported. Although CDR3 of ca1698 occupies the substrate binding site, ca1698 displays linear mixed inhibition kinetics instead of simple competitive inhibition kinetics. Two mechanisms are proposed to account for this apparent anomaly. Evidence for structural convergence of ca1698 and Nb216 during affinity maturation is also presented.     

SsrA (tmRNA) acts as an antisense RNA to regulate Staphylococcus aureus pigment synthesis by base pairing with crtMN mRNA

SsrA RNA (small stable RNA A), also known as tmRNA and 10Sa RNA, functions both as tRNA and mRNA through its unique structure. The carotenoid pigment is the eponymous feature of human pathogen Staphylococcus aureus. Here, we found that the pigment of the mutant strain with ssrA deletion was increased. Furthermore, it was demonstrated that ssrA could act as an antisense RNA aside from its well-known biological function, and crtMN, encoding two essential enzymes for the pigment synthesis, was identified as target of ssrA. Further investigation showed ssrA could specifically base pair with the RBS (ribosomal binding site) region of the crtMN mRNA. Our results revealed a new mechanism by which ssrA regulated the biosynthesis of pigment in S. aureus.     

Functional role for Tyr 31 in the catalytic cycle of chicken dihydrofolate reductase

Despite much work, many key aspects of the mechanism of the dihydrofolate reductase (DHFR) catalyzed reduction of dihydrofolate remain unresolved. In bacterial forms of DHFR both substrate and water access to the active site are controlled by the conformation of the mobile M20 loop. In vertebrate DHFRs only one conformation of the residues corresponding to the M20 loop has been observed. Access to the active site was proposed to be controlled by residue 31. MD simulations of chicken DHFR complexed with substrates and cofactor revealed a closing of the side chain of Tyr 31 over the active site on binding of dihydrofolate. This conformational change was dependent on the presence of glutamate on the para-aminobenzoylamide moiety of dihydrofolate. In its absence, the conformation remained open. Although water could enter the active site and hydrogen bond to N5 of dihydrofolate, indicating the feasibility of water as the proton donor, this was not controlled by the conformation of Tyr 31. The water accessibility of the active site was low for both conformations of Tyr 31. However, when hydride was transferred from NADPH to C6 of dihydrofolate before protonation, the average time during which water was found in hydrogen bonding distance to N5 of dihydrofolate in the active site increased almost fivefold. These results indicated that water can serve as the Broensted acid for the protonation of N5 of dihydrofolate during the DHFR catalyzed reduction.