Effect of cofactor binding and loop conformation on side chain methyl dynamics in dihydrofolate reductase

Dihydrofolate reductase (DHFR) has several flexible active site loops that facilitate ligand binding and catalysis. Previous studies of backbone dynamics in several complexes of DHFR indicate that the time scale and amplitude of motion depend on the conformation of the active site loops. In this study, information on dynamics is extended to methyl-containing side chains. To understand the role of side chain dynamics in ligand binding and loop conformation, methyl deuterium relaxation rates of Escherichia coli DHFR in binary folate and ternary folate:NADP+ complexes have been measured, together with chi(1) rotamer populations for threonine, isoleucine, and valine residues, determined from measurements of 3J(CgammaCO) and 3J(CgammaN) coupling constants. The results indicate that, in addition to backbone motional restriction in the adenosine-binding site, side chain flexibility in the active site and the surrounding active site loops is diminished upon binding NADP+. Resonances for several methyls in the active site and the surrounding active site loops were severely broadened in the folate:NADP+ ternary complex, suggesting the presence of motion on the chemical shift time scale. The side chains of Ile14 and Ile94, which pack against the nicotinamide and pterin rings of the cofactor and substrate, respectively, exhibit rotamer disorder in the ternary folate:NADP+ complex. Conformational fluctuations of these side chains may play a role in transition state stabilization; the observed line broadening for Ile14 suggests motions on a microsecond/millisecond time scale.     

Conformational relaxation following hydride transfer plays a limiting role in dihydrofolate reductase catalysis

The catalytic cycle of an enzyme is frequently associated with conformational changes that may limit maximum catalytic throughput. In Escherichia coli dihydrofolate reductase, release of the tetrahydrofolate (THF) product is the rate-determining step under physiological conditions and is associated with an "occluded" to "closed" conformational change. In this study, we demonstrate that in dihydrofolate reductase the closed to occluded conformational change in the product ternary complex (E.THF.NADP (+)) also gates progression through the catalytic cycle. Using NMR relaxation dispersion, we have measured the temperature and pH dependence of microsecond to millisecond time scale backbone dynamics of the occluded E.THF.NADP (+) complex. Our studies indicate the presence of three independent dynamic regions, associated with the active-site loops, the cofactor binding cleft, and the C-terminus and an adjacent loop, which fluctuate into discrete conformational substates with different kinetic and thermodynamic parameters. The dynamics of the C-terminally associated region is pH-dependent (p K a < 6), but the dynamics of the active-site loops and cofactor binding cleft are pH-independent. The active-site loop dynamics access a closed conformation, and the accompanying closed to occluded rate constant is comparable to the maximum pH-independent hydride transfer rate constant. Together, these results strongly suggest that the closed to occluded conformational transition in the product ternary complex is a prerequisite for progression through the catalytic cycle and that the rate of this process places an effective limit on the maximum rate of the hydride transfer step.     

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.     

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.     

Structure and dynamics of the G121V dihydrofolate reductase mutant: lessons from a transition-state inhibitor complex

It is well known that enzyme flexibility is critical for function. This is due to the observation that the rates of intramolecular enzyme motions are often matched to the rates of intermolecular events such as substrate binding and product release. Beyond this role in progression through the reaction cycle, it has been suggested that enzyme dynamics may also promote the chemical step itself. Dihydrofolate reductase (DHFR) is a model enzyme for which dynamics have been proposed to aid in both substrate flux and catalysis. The G121V mutant of DHFR is a well studied form that exhibits a severe reduction in the rate of hydride transfer yet there remains dispute as to whether this defect is caused by altered structure, dynamics, or both. Here we address this by presenting an NMR study of the G121V mutant bound to reduced cofactor and the transition state inhibitor, methotrexate. NMR chemical shift markers demonstrate that this form predominantly adopts the closed conformation thereby allowing us to provide the first glimpse into the dynamics of a catalytically relevant complex. Based on (15)N and (2)H NMR spin relaxation, we find that the mutant complex has modest changes in ps-ns flexibility with most affected residues residing in the distal adenosine binding domain rather than the active site. Thus, aberrant ps-ns dynamics are likely not the main contributor to the decreased catalytic rate. The most dramatic effect of the mutation involves changes in μs-ms dynamics of the F-G and Met20 loops. Whereas loop motion is quenched in the wild type transition state inhibitor complex, the F-G and Met20 loops undergo excursions from the closed conformation in the mutant complex. These excursions serve to decrease the population of conformers having the correct active site configuration, thus providing an explanation for the G121V catalytic defect.     

Backbone dynamics of the ribonuclease binase active site area using multinuclear ((15)N and (13)CO) NMR relaxation and computational molecular dynamics

The nano-pico second backbone dynamics of the ribonuclease binase, homologous to barnase, is investigated with (15)N, (13)C NMR relaxation at 11.74 and 18.78 T and with a 1.1 ns molecular dynamics simulation. The data are compared with the temperature factors reported for the X-ray structure of this enzyme. The molecular dynamics and X-ray data correspond well and predict motions in the loops 56-61 and 99-104 that contain residues that specifically recognize substrate and are catalytic (His101), respectively. In contrast, the (15)N relaxation data indicate that these loops are mostly ordered at the nano-pico second time scale. Nano-pico second motions in the recognition loop 56-61 are evident from (13)CO-(13)C cross relaxation data, but the mobility of the catalytic loop 99-104 is not detected by (13)CO cross relaxation either. From the results of this and previous work [Wang, L., Pang, Y., Holder, T., Brender, J. R., Kurochkin, A., and Zuiderweg, E. R. P. (2001) Proc. Natl. Acad. Sci. U.S.A., 98, 7684-7689], the following dynamical characterization of the active site area of binase emerges: a beta sheet, rigid at all probed time scales, supports the catalytic residue Glu 72. Both substrate-encapsulating loops are mobile on both fast and slow time scales, but the fast motions of the loop which contains the other catalytic residue, His 101, as predicted by B-factors and computational molecular dynamics is not detected by NMR relaxation. This work strongly argues for the use of several measures in the study of protein dynamics.     

Conformational change of the methionine 20 loop of Escherichia coli dihydrofolate reductase modulates pKa of the bound dihydrofolate

We evaluate the pK(a) of dihydrofolate (H(2)F) at the N(5) position in three ternary complexes with Escherichia coli dihydrofolate reductase (ecDHFR), namely ecDHFR(NADP(+):H(2)F) in the closed form (1), and the Michaelis complexes ecDHFR(NADPH:H(2)F) in the closed (2) and occluded (3) forms, by performing free energy perturbation with molecular dynamics simulations (FEP/MD). Our simulations suggest that in the Michaelis complex the pK(a) is modulated by the Met20 loop fluctuations, providing the largest pK(a) shift in substates with a "tightly closed" loop conformation; in the "partially closed/open" substates, the pK(a) is similar to that in the occluded complex. Conducive to the protonation, tightly closing the Met20 loop enhances the interactions of the cofactor and the substrate with the Met20 side chain and aligns the nicotinamide ring of the cofactor coplanar with the pterin ring of the substrate. Overall, the present study favors the hypothesis that N(5) is protonated directly from solution and provides further insights into the mechanism of the substrate protonation.     

Functional role of a mobile loop of Escherichia coli dihydrofolate reductase in transition-state stabilization.

The function of a highly mobile loop in Escherichia coli dihydrofolate reductase was studied by constructing a mutant (DL1) using cassette mutagenesis that had four residues deleted in the middle section of the loop (Met16-Ala19) and a glycine inserted to seal the gap. This part of the loop involves residues 16-20 and is disordered in the X-ray crystal structures of the apoprotein and the NADP+ binary complex but forms a hairpin turn that folds over the nicotinamide moiety of NADP+ and the pteridine moiety of folate in the ternary complex [Bystroff, C., & Kraut, J. (1991) Biochemistry 30, 2227-2239]. The steady-state and pre-steady-state kinetics and two-dimensional 1H NMR spectra were analyzed and compared to the wild-type protein. The kinetics on the DL1 mutant enzyme show that the KM value for NADPH (5.3 microM), the KM for dihydrofolate (2 microM), the rate constant for the release of the product tetrahydrofolate (10.3 s-1), and the intrinsic pKa value (6.2) are similar to those exhibited by the wild-type enzyme. However, the hydride-transfer rate declines markedly from the wild-type value of 950 s-1 to 1.7 s-1 for the DL1 mutant and when taken with data for substrate binding indicates that the loop contributes to substrate flux by a factor of 3.5 x 10(4). Thus, the mobility of loop I may provide a mechanism of recruiting hydrophobic residues which can properly align the nicotinamide and pteridine rings for the hydride-transfer process (a form of transition-state stabilization).(ABSTRACT TRUNCATED AT 250 WORDS)     

Dynamics of the active site loops in catalyzing aminoacylation reaction in seryl and histidyl tRNA synthetases

Aminoacylation reaction is the first step of protein biosynthesis. The catalytic reorganization at the active site of aminoacyl tRNA synthetases (aaRSs) is driven by the loop motions. There remain lacunae of understanding concerning the catalytic loop dynamics in aaRSs. We analyzed the functional loop dynamics in seryl tRNA synthetase from Methanopyrus kandleri (^mkSerRS) and histidyl tRNA synthetases from Thermus thermophilus (^ttHisRS), respectively, using molecular dynamics. Results confirm that the motif 2 loop and other active site loops are flexible spots within the catalytic domain. Catalytic residues of the loops form a network of interaction with the substrates to form a reactive state. The loops undergo transitions between closed state and open state and the relaxation of the constituent residues occurs in femtosecond to nanosecond time scale. Order parameters are higher for constituent catalytic residues which form a specific network of interaction with the substrates to form a reactive state compared to the Gly residues within the loop. The development of interaction is supported from mutation studies where the catalytic domain with mutated loop exhibits unfavorable binding energy with the substrates. During the open-close motion of the loops, the catalytic residues make relaxation by ultrafast librational motion as well as fast diffusive motion and subsequently relax rather slowly via slower diffusive motion. The Gly residues act as a hinge to facilitate the loop closing and opening by their faster relaxation behavior. The role of bound water is analyzed by comparing implicit solvent-based and explicit solvent-based simulations. Loops fail to form catalytically competent geometry in absence of water. The present result, for the first time reveals the nature of the active site loop dynamics in aaRS and their influence on catalysis.     

NMR studies of the interaction of a type II dihydrofolate reductase with pyridine nucleotides reveal unexpected phosphatase and reductase activity

The interaction of type II R67 dihydrofolate reductase (DHFR) with its cofactor nicotinamide adenine dinucleotide phosphate (NADP(+)) has been studied using nuclear magnetic resonance (NMR). Doubly labeled [U-(13)C,(15)N]DHFR was obtained from Escherichia coli grown on a medium containing [U-(13)C]-D-glucose and (15)NH(4)Cl, and the 16 disordered N-terminal amino acids were removed by treatment with chymotrypsin. Backbone and side chain NMR assignments were made using triple-resonance experiments. The degeneracy of the amide (1)H and (15)N shifts of the tetrameric DHFR was preserved upon addition of NADP(+), consistent with kinetic averaging among equivalent binding sites. Analysis of the more titration-sensitive DHFR amide resonances as a function of added NADP(+) gave a K(D) of 131 +/- 50 microM, consistent with previous determinations using other methodology. We have found that the (1)H spectrum of NADP(+) in the presence of the R67 DHFR changes as a function of time. Comparison with standard samples and mass spectrometric analysis indicates a slow conversion of NADP(+) to NAD(+), i.e., an apparent NADP(+) phosphatase activity. Studies of this activity in the presence of folate and a folate analogue support the conclusion that this activity results from an interaction with the DHFR rather than a contaminating phosphatase. (1)H NMR studies of a mixture of NADP(+) and NADPH in the presence of the enzyme reveal that a ternary complex forms in which the N-4A and N-4B nuclei of the NADPH are in the proximity of the N-4 and N-5 nuclei of NADP(+). Studies using the NADP(+) analogue acetylpyridine adenosine dinucleotide phosphate (APADP(+)) demonstrated a low level of enzyme-catalyzed hydride transfer from NADPH. Analysis of DHFR backbone dynamics revealed little change upon binding of NADP(+). These additional catalytic activities and dynamic behavior are in marked contrast to those of type I DHFR.     

Changes in structure and dynamics of the Fv fragment of a catalytic antibody upon binding of inhibitor

Binding of the product inhibitor p-nitrophenol to the monoclonal esterolytic antibody NPN43C9 has been investigated by performing NMR spectroscopy of the heterodimeric variable-domain fragment (Fv) of the antibody in the presence and absence of inhibitor. Structural information from changes in chemical shift upon binding has been related to the changes in local dynamics in the active site of the catalytic antibody using NMR relaxation measurements. Significant changes in the chemical shifts of the backbone resonances upon binding extend beyond the immediate vicinity of the antigen binding site into the interface between the two associated polypeptides that form the Fv heterodimer, a possible indication that the binding of ligand causes a change in the relative orientations of the component light (V(L)) and heavy (V(H)) chain polypeptides. Significant differences in backbone dynamics were observed between the free Fv and the complex with p-nitrophenol. A number of resonances, including almost all of the third hypervariable loop of the light chain (L3), were greatly broadened in the free form of the protein. Other residues in the antigen-binding site showed less broadening of resonances, but still required exchange terms (R(ex)) in the model-free dynamics analysis, consistent with motion on a slow timescale in the active site region of the free Fv. Binding of p-nitrophenol caused these resonances to sharpen, but some R(ex) terms are still required in the analysis of the backbone dynamics. We conclude that the slow timescale motions in the antigen-binding site are very different in the bound and free forms of the Fv, presumably due to the damping of large-amplitude motions by the bound inhibitor.     

PhosphoThr peptide binding globally rigidifies much of the FHA domain from Arabidopsis receptor kinase-associated protein phosphatase

A net increase in the backbone rigidity of the kinase-interacting FHA domain (KI-FHA) from the Arabidopsis receptor kinase-associated protein phosphatase (KAPP) accompanies the binding of a phosphoThr peptide from its CLV1 receptor-like kinase partner, according to (15)N NMR relaxation at 11.7 and 14.1 T. All of the loops of free KI-FHA display evidence of nanosecond-scale motions. Many of these same residues have residual dipolar couplings that deviate from structural predictions. Binding of the CLV1 pT868 peptide seems to reduce nanosecond-scale fluctuations of all loops, including half of the residues of recognition loops. Residues important for affinity are found to be rigid, i.e., conserved residues and residues of the subsite for the key pT+3 peptide position. This behavior parallels SH2 and PTB domain recognition of pTyr peptides. PhosphoThr peptide binding increases KI-FHA backbone rigidity (S(2)) of three recognition loops, a loop nearby, seven strands from the beta-sandwich, and a distal loop. Compensating the trend of increased rigidity, binding enhances fast mobility at a few sites in four loops on the periphery of the recognition surface and in two loops on the far side of the beta-sandwich. Line broadening evidence of microsecond- to millisecond-scale fluctuations occurs across the six-stranded beta-sheet and nearby edges of the beta-sandwich; this forms a network connected by packing of interior side chains and H-bonding. A patch of the slowly fluctuating residues coincides with the site of segment-swapped dimerization in crystals of the FHA domain of human Chfr. Phosphopeptide binding introduces microsecond- to millisecond-scale fluctuations to more residues of the long 8/9 recognition loop of KI-FHA. The rigidity of this FHA domain appears to couple as a whole to pThr peptide binding.     

Tissue inhibitor of metalloproteinases-1 undergoes microsecond to millisecond motions at sites of matrix metalloproteinase-induced fit

The N-terminal, matrix metalloproteinase (MMP)-inhibitory fragment of recombinant, human tissue inhibitor of metalloproteinases (TIMP-1) exhibits varied backbone dynamics and rigidity. Most striking is the presence of chemical exchange in the MMP-binding ridge reported to undergo conformational change upon MMP binding. Conformational exchange fluctuations in microseconds to milliseconds map to the sites of MMP-induced fit at residues Val29 through Leu34 of the AB loop and to the Ala65 and Cys70 "hinges" of the CD loop of TIMP-1. Slow chemical exchange is also present at the type I turn of the EF loop at the base of the MMP-binding ridge. These functional slow motions and other fast internal motions are evident from backbone (15)N spin relaxation at 500 and 750 MHz, whether interpreted by the model-free formalism with axial diffusion anisotropy or by the reduced spectral density approach. The conformational exchange is confirmed by its deviation from the trend between R(2) and the cross-correlation rate eta. The magnetic field-dependence indicates that the chemical exchange broadening in the AB and CD loops is fast on the time-scale of chemical shift differences. The conformational exchange rates for most of these exchanging residues, which can closely approach MMP, appear to be a few thousand to several thousand per second. The slow dynamics of the TIMP-1 AB loop contrast the picosecond to nanosecond dynamics reported in the longer TIMP-2 AB loop.     

Dynamics of the metallo-beta-lactamase from Bacteroides fragilis in the presence and absence of a tight-binding inhibitor

A significant determinant for the broad substrate specificity of the metallo-beta-lactamases from Bacteroides fragilis and other similar organisms is the presence of a plastic substrate binding site that is nevertheless capable of tight substrate binding in the Michaelis complex. To achieve these two competing ends, the molecule apparently employs a flexible flap that closes over the active site in the presence of substrate. These characteristics imply that dynamic changes are an important component of the mechanism of action of these enzymes. The backbone and tryptophan side chain dynamics of the metallo-beta-lactamase from B. fragilis have been examined using (15)N NMR relaxation measurements. Two states of the protein were examined, in the presence and absence of a tight-binding inhibitor. Relaxation measurements were analyzed by the model-free method. Overall, the metallo-beta-lactamase molecule is rigid and shows little flexibility except in loops. The flexibility of the loop that covers the active site is not unusually great as compared to the other loops of the protein. Local motion on a picosecond time scale was found to be very similar throughout the protein in the presence and absence of the inhibitor, but a significant difference was observed in the motions on a nanosecond time scale (tau(e)). Large-amplitude motions with a time constant of about 1.3 ns were observed for the flexible flap region (residues 45-55) in the absence of the inhibitor. These motions were completely damped out in the presence of the inhibitor. In addition, the motion of a tryptophan side chain at the tip of the beta-hairpin of the flap shows a very significant difference in motion on the ps time scale. These results indicate that the motions of the polypeptide chain in the flap region can be invoked to explain both the wide substrate specificity (the free form has considerable amplitude of motion in this region) and the catalytic efficiency of the metallo-beta-lactamase (the motions are damped out when the inhibitor and by implication a substrate binds in the active site).     

Gating of the active site of triose phosphate isomerase: Brownian dynamics simulations of flexible peptide loops in the enzyme.

The enzyme triose phosphate isomerase has flexible peptide loops at its active sites. The loops close over these sites upon substrate binding, suggesting that the dynamics of the loops could be of mechanistic and kinetic importance. To investigate these issues, the loop motions in the dimeric enzyme were simulated by Brownian dynamics. The two loops, one on each monomer, were represented by linear chains of appropriately parameterized spheres, each sphere corresponding to an amino acid residue. The loops moved in the electrostatic field of the rest of the enzyme, which was held rigid in its crystallographically observed conformation. In the absence of substrate, the loops exhibited gating of the active site with a period of about 1 ns and occupied "closed" conformations for about half of the time. As the period of gating is much shorter than the enzyme-substrate relaxation time, the motion of the loops does not reduce the rate constant for the approach of substrate from its simple diffusion-controlled value. This suggests that the flexible loops may have evolved to create the appropriate environment for catalysis while, at the same time, minimizing the kinetic penalty for gating the active site.     

Insight derived from molecular dynamics simulation into dynamics and molecular motions of cuticle-degrading serine protease Ver112

Cuticle-degrading serine protease Ver112, which derived from a nematophagous fungus Lecanicillium psalliotae, has been exhibited to have high cuticle-degrading and nematicidal activities. We have performed molecular dynamics (MD) simulation based on the crystal structure of Ver112 to investigate its dynamic properties and large-scale concerted motions. The results indicate that the structural core of Ver112 shows a small fluctuation amplitude, whereas the substrate binding sites, and the regions close to and opposite the substrate binding sites experience significant conformational fluctuations. The large concerted motions obtained from essential dynamics (ED) analysis of MD trajectory can lead to open or close of the substrate binding sites, which are proposed to be linked to the functional properties of Ver112, such as substrate binding, orientation, catalytic, and release. The significant motion in the loop regions that is located opposite the binding sites are considered to play an important role in modulating the dynamics of the substrate binding sites. Furthermore, the bottom of free energy landscape (FEL) of Ver112 are rugged, which is mainly caused by the fluctuations of substrate binding regions and loops located opposite the binding site. In addition, the mechanism underlying the high flexibility and catalytic activity of Ver112 was also discussed. Our simulation study complements the biochemical and structural studies, and provides insight into the dynamics-function relationship of cuticle-degrading serine protease Ver112.     

Coupling between Catalytic Loop Motions and Enzyme Global Dynamics

Catalytic loop motions facilitate substrate recognition and binding in many enzymes. While these motions appear to be highly flexible, their functional significance suggests that structure-encoded preferences may play a role in selecting particular mechanisms of motions. We performed an extensive study on a set of enzymes to assess whether the collective/global dynamics, as predicted by elastic network models (ENMs), facilitates or even defines the local motions undergone by functional loops. Our dataset includes a total of 117 crystal structures for ten enzymes of different sizes and oligomerization states. Each enzyme contains a specific functional/catalytic loop (10–21 residues long) that closes over the active site during catalysis. Principal component analysis (PCA) of the available crystal structures (including apo and ligand-bound forms) for each enzyme revealed the dominant conformational changes taking place in these loops upon substrate binding. These experimentally observed loop reconfigurations are shown to be predominantly driven by energetically favored modes of motion intrinsically accessible to the enzyme in the absence of its substrate. The analysis suggests that robust global modes cooperatively defined by the overall enzyme architecture also entail local components that assist in suitable opening/closure of the catalytic loop over the active site.     

Crystal structures of Escherichia coli dihydrofolate reductase: the NADP+ holoenzyme and the folate.NADP+ ternary complex. Substrate binding and a model for the transition state

The crystal structure of dihydrofolate reductase (EC 1.5.1.3) from Escherichia coli has been solved as the binary complex with NADP+ (the holoenzyme) and as the ternary complex with NADP+ and folate. The Bragg law resolutions of the structures are 2.4 and 2.5 A, respectively. The new crystal forms are nonisomorphous with each other and with the methotrexate binary complex reported earlier [Bolin, J. T., Filman, D. J., Matthews, D. A., Hamlin, R. C., & Kraut, J. (1982) J. Biol. Chem. 257, 13650-13662]. In general, NADP+ and folate binding conform to predictions, but the nicotinamide moiety of NADP+ is disordered in the holoenzyme and ordered in the ternary complex. A mobile loop (residues 16-20) involved in binding the nicotinamide is also disordered in the holoenzyme. We report a detailed analysis of the binding interactions for both ligands, paying special attention to several apparently strained interactions that may favor the transition state for hydride transfer. Hypothetical models are presented for the binding of 7,8-dihydrofolate in the Michaelis complex and for the transition-state complex.     

Role of water in the catalytic cycle of E. coli dihydrofolate reductase

Dihydrofolate reductase (DHFR) catalyzes the nicotinamide adenine dinucleotide phosphate (NADPH)-dependent reduction of 7,8-dihydrofolate (H2F) to 5,6,7,8-tetrahydrofolate (H4F). Because of the absence of any ionizable group in the vicinity of N5 of dihydrofolate it has been proposed that N5 could be protonated directly by a water molecule at the active site in the ternary complex of the Escherichia coli enzyme with cofactor and substrate. However, in the X-ray structures representing the Michaelis complex of the E. coli enzyme, a water molecule has never been observed in a position that could allow protonation of N5. In fact, the side chain of Met 20 blocks access to N5. Energy minimization reported here revealed that water could be placed in hydrogen bonding distance of N5 with only minor conformational changes. The r.m.s. deviation between the conformation of the M20 loop observed in the crystal structures of the ternary complexes and the conformation adopted after energy minimization was only 0.79 A. We performed molecular dynamics simulations to determine the accessibility by water of the active site of the Michaelis complex of DHFR. Water could access N5 relatively freely after an equilibration time of approximately 300 psec during which the side chain of Met 20 blocked water access. Protonation of N5 did not increase the accessibility by water. Surprisingly the number of near-attack conformations, in which the distance between the pro-R hydrogen of NADPH and C6 of dihydrofolate was less than 3.5 A and the angle between C4 and the pro-R hydrogen of NADPH and C6 of dihydrofolate was greater than 120 degrees, did not increase after protonation. However, when the hydride was transferred from NADPH to C6 of dihydrofolate before protonation, the side chain of Met 20 moved away from N5 after approximately 100 psec thereby providing water access. The average time during which water was found in hydrogen bonding distance to N5 was significantly increased. These results suggest that hydride transfer might occur early to midway through the reaction followed by protonation. Such a mechanism is supported by the very close contact between C4 of NADP+ and C6 of folate observed in the crystal structures of the ternary enzyme complexes, when the M20 loop is in its closed conformation.     

Solution structure of human sorting nexin 22

The sorting nexins (SNXs) constitute a large group of PX domain-containing proteins that play critical roles in protein trafficking. We report here the solution structure of human sorting nexin 22 (SNX22). Although SNX22 has <30% sequence identity with any PX domain protein of known structure, it was found to contain the alpha/beta fold and compact structural core characteristic of PX domains. Analysis of the backbone dynamics of SNX22 by NMR relaxation measurements revealed that the two walls of the ligand binding cleft undergo internal motions: on the picosecond timescale for the beta1/beta2 loop and on the micro- to millisecond timescale for the loop between the polyproline motif and helix alpha2. Regions of the SNX22 structure that differ from those of other PX domains include the loop connecting strands beta1 and beta2 and the loop connecting helices alpha1 and alpha2, which appear to be more mobile than corresponding loops in other known structures. The interaction of dibutanoyl-phosphatidylinositol-3-phosphate (dibutanoyl-PtdIns(3)P) with SNX22 was investigated by an NMR titration experiment, which identified the binding site in a basic cleft and indicated that ligand binding leads only to a local structural rearrangement as has been found with other PX domains. Because motions in the loops are damped out when dibutanoyl-PtdIns(3)P binds, entropic effects could contribute to the lower affinity of SNX22 for this ligand compared to other PX domains.