|
Advertisement | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 283, Issue 8, 4690-4698, February 22, 2008
A Balancing Act between Net Uptake of Water during Dihydrofolate Binding and Net Release of Water upon NADPH Binding in R67 Dihydrofolate Reductase*
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
R67 DHFR is a homotetramer with a single active site pore; the overall structure possesses 222 symmetry as seen in Fig. 1 (5). The symmetry of the active site results in overlapping binding sites for substrate, DHF, and cofactor, NADPH. This can be observed as R67 DHFR binds a total of two ligands as follows: either two NADPH molecules or two folate/DHF molecules or one NADPH plus one folate/DHF molecule (6). The first two complexes are dead-end (binary) complexes, whereas the third is the productive ternary complex. Because of the 222 symmetry, binding to either ligand is unlikely to be optimal.
The active site pore of R67 DHFR is unusual in its hourglass shape as well as its large size (2938 Å3 for apoR67 DHFR with hydrogens added, calculated by CASTp (4, 7)). Because of this large volume, DHF and NADPH cannot occupy all the space in the pore and must use water to mediate some contacts with the protein.
How do substrate and cofactor bind to R67 DHFR? From NMR, crystallography, and docking studies, the pteridine ring of DHF/folate binds at the center of the pore, but the para-aminobenzoylglutamate tail of dihydrofolate/folate is disordered (5, 8, 9). Interactions have been predicted between symmetry-related Lys-32 residues (at either edge of the pore) and the glutamate tail of substrate. Also, increasing concentrations of sucrose or trehalose (10) increase Km (DHF) as well as Kd (DHF). As polyols can have a variety of effects, including perturbation of solvent structure, a possibility is that these compounds affect the water structure and, in turn, DHF binding.
If water is involved in an interaction, perturbation of water content should affect binding. Thus we probed the role of water in DHF/folate binding to R67 DHFR, using osmolytes with different properties to determine whether the primary effect was because of water or some other variable.
To provide a complete picture, the effects of osmolytes on the interaction of NADPH with R67 DHFR were also monitored. Previous crystallography and NMR studies (8, 9, 11, 12) show NADPH binds in an extended conformation, with numerous specific interactions, including H-bonds between the nicotinamide carboxamide of NADP+ with backbone NH and O atoms from Ile-68. Ionic interactions also form between the nicotinamide phosphate and the adenosyl-2'-phosphate with symmetry-related Lys-32 residues (13–16).
Because of the 222 symmetry of the active site pore, DHF and NADPH bind to the same promiscuous surface. Using various osmolytes, this study finds water release accompanies NADPH binding, whereas water uptake accompanies DHF binding.
|
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
Steady State Kinetics—Steady state kinetic data were obtained using a PerkinElmer Life Sciences
35 spectrophotometer interfaced with an IBM personal computer as described previously (19). MTA (MTA buffer: 100 mM Tris, 50 mM MES, 50 mM acetic acid polybuffer plus 1 mM β-mercaptoethanol) polybuffer was used. This buffering system maintains a constant ionic strength from pH 4.5 to 9.5 (20). The steady state kinetic rates were measured at 30 °C by the addition of substrate and cofactor, followed by the addition of enzyme to initiate the reaction. To obtain kcat and Km values, the concentration of NADPH was held constant at a subsaturating level, whereas the concentration of DHF was varied. This process was repeated using four additional subsaturating concentrations of NADPH. The data were fit globally to the nonlinear bisubstrate Michaelis-Menten equation utilizing SAS (statistical analysis software (21, 22)). The NLINEK macro for use in SAS is available on line. For those conditions where Km (NADPH) did not change, saturating concentrations of NADPH were used, and DHF concentrations varied to obtain kcat and Km (DHF) values.
Water Activity Measurements—A Westcor 5500 vapor pressure osmometer was used to obtain the osmolality of the solutions. This value was converted into water activity using Equation 1,
![]() |
The binding of DHF to R67 DHFR·NADPH in a buffer containing neutral osmolytes (S) can be described by Equation 2,
![]() |
H2O and
S are the stoichiometric coefficients of water and osmolyte, respectively. From Wyman (24), it can be calculated that ln Ka (or ln kcat/Km (DHF)) is related to ln as according to the following:
H2O and
S, the stoichiometric coefficients of water and osmolyte, as shown in Equation 3,
![]() |
H2O would contribute more significantly to Equation 3 than
S. This scenario can be tested by use of osmolytes from different classes, which differentially associate with proteins (and presumably ligands as well (25, 26)).
Isothermal Titration Calorimetry—Affinities, stoichiometries, as well as
H values associated with binding were determined using isothermal titration calorimetry (ITC) as described previously (6). Measurements were performed on a VP-ITC microcalorimeter from MicroCal interfaced to a Gateway personal computer for data acquisition and analysis. Origin version 5 scientific software was used to analyze the data. This instrument has been described previously (28). R67 DHFR concentrations typically ranged from 60 to 150 µM in MTA buffer (pH 7). For titrations with osmolyte present, MTA buffer plus osmolyte was used in the reference cell. The "c value" (=[Ptotal]/Kd) for binding was from 3 to 127, within the suggested range of 1–1000 (28).
The heat capacity (
Cp) or the change in
H as a function of temperature is described by Equation 4,
![]() |
Escherichia coli Growth Conditions—The ability of the E. coli strain DH5
to grow on M9 minimal media (29) containing 0.02% casamino acids with 20 µg of TMP/ml was assessed. Another layer of screening added increasing concentrations of sorbitol to the media. The water activity of the plates was measured using an Aqualab meter (Decagon Devices) and converted into osmolality using Equation 1.
| RESULTS |
|---|
|
|
|---|
With this background in mind, the steady state kinetic behavior of R67 DHFR was monitored in the presence of the neutral osmolytes glycerol, ethylene glycol, trimethylamine N-oxide (TMAO), dimethyl sulfoxide (Me2SO), glycine betaine, and sucrose. Little to no effect (<1.5-fold) is observed on the catalytic rate constant, kcat, whereas the Km for DHF increases in the presence of all these osmolytes. A possible decrease in Km for NADPH is also observed. If a linear relationship is observed in plots of either osmolality or ln water activity versus ln kcat/Km (DHF), then effects on water are involved (30, 38, 40, 41). Fig. 2 shows the linear relationships associated with plots of osmolality versus ln kcat/Km (DHF). As a precaution, we also plotted effects on solution viscosity or solution dielectric and overlapping data were not observed. These figures are available as supplemental Figs. 1 and 2. Another control experiment monitored the effect of osmolytes on the pH-dependent dissociation of tetramer to two protonated dimers. Addition of 10% ethylene glycol, 12.5% Me2SO, or 20% PEG400 did not greatly alter the pH titration data (experimental design in supplemental material as well as supplemental Fig. 3), indicating that osmolyte addition minimally perturbs the oligomeric state of R67 DHFR.
|
nW, involved in DHF binding to the enzyme·NADPH complex and going to the transition state. The slopes associated with the individual osmolytes for this type of plot are given in Table 1. The slopes associated with TMAO, glycerol and ethylene glycol are low, around 16–25, and the data overlap reasonably well. A fit to these three combined data sets yields a slope of 17 ± 2. However, the slopes for the glycine betaine, Me2SO, and sucrose data are higher, ranging from 40 to 60. Variable slopes are common in osmolality studies, but their origin is not clear (37, 38, 41–54). This issue will be broached under the "Discussion."
|
nW/55.6, the slope of the combined betaine and sucrose data were converted to a
nW of 29 ± 3. This value is lower than the effect of sucrose and betaine on kcat/Km (DHF) and suggests that Km may contain kinetic terms. We also present the data as a plot of total heat (Qtotal) versus [DHF] at increasing concentrations of betaine (see Fig. 4A).
Although the data for betaine and sucrose overlap in Fig. 3A, steeper slopes were observed using PEG400 (
nW = 77 ± 11), indicating DHF binding remains sensitive to osmolyte identity. We also performed titrations with Me2SO, glycerol, and ethylene glycol, and at low cosolvent concentrations, decreases in Ka were observed. However, at higher osmolyte concentrations, good fits were not obtained; thus these data are not presented.
From Fig. 4A, as the osmolyte concentration increases, the observed enthalpy value becomes smaller. Less negative
H values associated with water release have been observed previously as well as more negative
H values associated with water uptake (55, 56). Water-mediated enthalpy-entropy compensation may be the mechanism by which changes in enthalpy occur (57–60).
We additionally monitored osmolyte effects on folate binding to enzyme (i.e. formation of the enzyme·2folate complex). Large effects are observed, particularly on the enthalpy change, making fitting difficult. Thus the data are presented in Fig. 4B as a plot of Qtotal versus [folate].
To determine whether osmolytes affect NADPH binding, several ITC titrations were performed in varying concentrations of betaine (0–20%). Fig. 4C shows a plot of Qtotal versus [NADPH]. An opposite trend compared with DHF/folate binding is noted in these plots, with increasing osmolyte concentrations tightening binding. Similar trends are noted when ethylene glycol, Me2SO, or PEG400 are used. The data were fit to a sequential sites model, and a plot of ln Ka1 versus osmolality is shown in Fig. 3B. All the data overlay reasonably well. Converting the slope of this plot to
nW yields a value of –38 ± 6. This result indicates that release of water accompanies binding of NADPH, in contrast to water uptake upon binding of DHF. These opposite results are surprising as R67 DHFR possesses 222 symmetry and DHF shares a related binding site with NADPH. Also of note, we find that different osmolytes provide overlapping data when NADPH binding is probed whereas addition of different solutes yields variable slopes for DHF binding.
|
1.9 osM (67, 68) in minimal medium. Our above studies indicate that increasing osmolality negatively affects the ability of R67 DHFR to bind DHF and proceed to the transition state. Does this have any impact with respect to function of R67 DHFR in bacteria? We addressed this question by adding increasing concentrations of sorbitol (69, 70) to M9 minimal media containing TMP and asked whether this osmolyte impacted growth of E. coli strain DH5
, which has been transformed with wild type or mutant R67 DHFR clones. This approach appeared feasible as well as unique because few reports describe uptake of water upon binding, and those that do may involve conformational changes (36, 56, 71–76). (Note: This enzyme is a homotetramer; when a single residue is mentioned, all four related residues are implied.)
|
did not grow on minimal M9 control plates containing 20 µg of TMP/ml as this antibacterial agent selectively inhibits chromosomal DHFR. Transformation of pUC8 carrying the wild type (WT) R67 DHFR gene into DH5
allows growth on this selective media. Transformation of mutant R67 DHFR clones into DH5
allowed a range of growth patterns on TMP plates. The WT gene and mutants with high kcat/Km (DHF) values allowed confluent growth overnight. In contrast, the Y69L mutant (with a lower catalytic efficiency) shows growth after 3 days. The K32M mutant has a very low catalytic efficiency and low protein expression levels; it is unable to rescue DH5
from TMP-selective pressure. In the second set of control plates, all cultures were streaked on M9 minimal agar containing increasing concentrations of sorbitol to determine the effect of osmolality on cell growth. Cells showed confluent growth overnight until 1.95 osmol was reached, at which point 2 days were required for growth. Higher concentrations of sorbitol were unable to support growth.
Next, we asked whether sorbitol could impact the ability of DH5
carrying various R67 DHFR clones to grow on minimal M9 plates containing 20 µg of TMP/ml. The number of days required to obtain visible growth on these plates at 37 °C is listed in supplemental Table S1. The WT R67 DHFR clone allows growth on all sorbitol conditions until the osmolyte concentration becomes too high (or cell water content becomes too low (42)) and cell growth stops. The Q67H mutant, which has a reasonable catalytic efficiency, but displays substrate and cofactor inhibition (77), allows growth up to 1.81 osmol conditions. The I68M and Y69L mutants, with lower kcat/Km (DHF) values, allow growth up to 1.44 and 0.81 osmol conditions, respectively. Because DHFR activity is required to restore folate end product prototrophy and as cell growth parallels to some degree the efficiency of the mutants, it appears that lowering the intracellular water activity by sorbitol addition confers another level to the selection process. To our knowledge, this is the first time osmolyte addition has been successfully used as a selection tool to probe for the correlation between in vitro water uptake and in vivo function.
Heat Capacity—Because water uptake has been shown to be involved in DHF binding and water release associated with NADPH binding, we asked whether water reorganization might contribute to the binding enthalpy. As Chervenak and Toone (57) have shown that water reorganization correlates with heat capacity values, we monitored the change in enthalpy as a function of temperature. For NADPH binding, the Kd value decreases, accompanied by a
H change. The data for
H for the first NADPH site are shown in Fig. 5A. The slope (
Cp) is –178 ± 15 cal/K mol. (There are changes in
H for the second NADPH site as well, but the errors are larger.) We note that NADPH binding to R67 DHFR occurs without release or uptake of protons, as does binding of DHF to enzyme·NADP+, so protonation effects do not contribute to the
H change (14).
|
H are observed. This plot is depicted in Fig. 5B, and heat capacity of –199 ± 16 cal/K mol is calculated. Both
Cp values for DHF and NADPH binding are small, negative, and within error of each other. Because a number of the interactions between NADPH and protein or DHF and protein are similar (8, 15, 78, 79), the convergence of the
Cp values may describe at some level the ability of symmetry-related sites to accommodate the two different ligands.
From Fig. 5A, NADPH binding is enthalpy- and entropy-driven and the entropy term is positive (favorable) until 281 K (8 °C). Above this temperature, the entropy term becomes negative, and binding becomes enthalpy-driven. A similar scenario occurs for DHF binding, although the point at which
S changes sign occurs at a higher temperature, 291 K (18 °C). The changes in
G are much smaller than those in
H and T
S, consistent with enthalpy-entropy compensation. Dunitz (58) has proposed that enthalpy-entropy compensation commonly occurs in weak interactions. Also this type of compensation is a hallmark of water involvement (80, 81) as calculations indicate that the
H value associated with cavity formation for accommodating a solvent molecule is exactly balanced by the entropy of the cavity (82, 83). Thus, contributions to
H and
S can occur but will not necessarily show up in
G as
Gsolvation equals zero (=
Hsolvation – T
Ssolvation). Similar results have been observed previously in R67 DHFR upon titrating folate into R67 DHFR·NADPH in the presence of increasing salt concentrations. A titration in the
H value was observed but not in
G (15). Solvent-accessible area calculations have been performed and are provided in the supplemental material.
| DISCUSSION |
|---|
|
|
|---|
Variable
nW Values—Because the slope values describing kcat/Km (DHF) vary in Fig. 2, several potential scenarios could be occurring. The first possibility is that the different slopes arise by preferential binding of the osmolyte upon the protein surface (37, 42). From Equation 3, if an osmolyte is preferentially bound by the protein, then the net change in water molecules associated with ligand binding will be overestimated if water is taken up (or underestimated if water is released). In other words, the slope can be generated by any combination of numbers for
H2O and
solute with the caveat that
solute is a positive value (37). For example, a "
nW"of +17 could be generated by 17H2O + 0solute or 11H2O + 6solute or 5H2O + 12solute, etc.
At first glance, binding of solutes does seem possible as R67 DHFR possesses a promiscuous binding surface and could be imagined to weakly bind other molecules as well. However, NADPH binding serves as an internal control. If osmolytes do bind and provide a range of slopes for the DHF-binding site, why isn't the same phenomenon observed for NADPH? In other words, because DHF and NADPH bind at symmetry-related sites, we might expect both slopes in Fig. 2 and Fig. 3B to be sensitive to osmolyte identity. At 1–2 osmol, Fig. 2 shows the kcat/Km (DHF) values clearly separating with different osmolytes, whereas at
2 osmol in Fig. 3B, the Ka values describing NADPH binding overlap.
A second option that could describe the variable slopes is crowding or volume exclusion (43–45). This term describes the volume of the protein probed plus the volume swept out by a sphere with radius r rolled over the surface. Because of different radii, the excluded volume varies for each cosolvent with larger osmolytes probing larger surface features. PEGs of varying sizes were used to analyze this possibility. The average radius of gyration for PEG400 is 8.1 Å; however, its radius increases because of hydration; the radii probed vary from 9 to 15 Å, depending on the macromolecule studied (84). This range of values contrasts with a maximal pore radius in R67 DHFR of 12 Å (5). When PEG400, PEG3350, and PEG8000 are used to probe DHF binding, the slopes of the ln aH2O versus ln kcat/Km(DHF) plot increase to 78, 145, and 353, respectively (Table 1). The increase in slope indicates the PEGs are probing larger surface features as the osmolyte becomes too big to fit in the active site pore. Supplemental Fig. 3 shows a plot of ln molar volume (Vmol) versus the slopes listed in Table 1, indicating crowding plays some role, particularly with larger molecular weight PEGs. However, if crowding causes the variable slopes associated with DHF binding in Fig. 2, then one might expect crowding to affect NADPH binding in a similar fashion. Because this is not observed, crowding is probably not the origin of the behavior.
Another possibility is that conformational changes occur in the protein and/or ligands such that release of water upon ligand binding is masked by water uptake associated with the conformational change, which results in a larger surface area. Examples of this phenomenon include glucose binding to hexokinase (85), aspartate to aspartate transcarbamoylase (71), DNA to repressors (30), etc. However NMR studies of R67 DHFR indicate little change in the generalized order parameter upon binding NADP+ (
Sfree2
= 0.89,
Sbound2
0.86), consistent with a fairly rigid structure (12). Small changes may occur when the apo and ternary complex structures are compared, as introduction of asymmetry, particularly at Gln-67 and Tyr-69 residues, is observed (8).
Another origin of the different slope effects could lie in the fact that DHF binds to an enzyme·NADPH complex, whereas NADPH binds to the apoenzyme. However, binding of a second NADPH molecule to enzyme·NADPH still shows water release and overlapping data with varying osmolytes (
nW =–20 ± 8; data not shown), and binding of two folate molecules to the free enzyme shows very large effects (Fig. 4B). Both of these observations suggest it is the identity of the ligands rather than the position in the binding sequence (i.e. first versus second molecule) that causes the different behavior.
Because NADPH and DHF compete for binding to related sites in R67 DHFR, it seems likely that osmolyte sensitivity is either reporting on differences in how the ligands bind, differences in solvation of the two ligands (86, 87), differences in solvent that may be induced by osmolyte addition, and/or differences in protein solvation levels (88, 89). For example, for the latter possibility, Dzingeleski and Wolfenden (72) saw an uptake of 9 water molecules associated with substrate and inhibitor binding in adenosine deaminase. They proposed that a hydrated form of this enzyme is in equilibrium with a dehydrated form and that only the hydrated form can bind substrate. A similar mechanism could exist in R67 with NADPH preferentially binding to the dehydrated form and DHF preferentially interacting with the hydrated form. At the moment, we are unable to differentiate between the models proposed in this paragraph.
Correlation of in Vitro Results with in Vivo Function?—Given the uncertainty in interpreting the variable slopes associated with DHF binding, we note that what is important is that the results of our in vitro experiments correlate with the in vivo efficacy of R67 DHFR. For example, if we consider the Y69L mutant clone, its catalytic efficiency and intracellular protein concentration appear barely sufficient to allow DH5
to grow in the presence of TMP. This is not surprising as its Km (DHF) equals 180 µM (90); however, the normal intracellular DHF concentration in E. coli is
300 µM (91, 92); thus it is likely working under kcat/Km conditions. From Fig. 2, these are exactly the conditions where we might expect to observe the greatest osmolyte effect. If betaine is one of the intracellular osmolytes produced in response to osmotic stress (42), there should be a large inhibitory effect on R67 DHFR efficiency. Thus addition of sorbitol increases the osmotic pressure, cellular water content decreases, and the Y69L mutant is no longer able to support growth on TMP plates.
Extrapolating from our in vitro studies to in vivo conditions, we also might expect increasing osmolyte concentrations to tighten NADPH binding to both sites, potentially leading to inhibition by the two NADPH complex (see Fig. 5C). Concurrently, for the two substrate complex, we would predict that it is less likely to form in the cell (Fig. 5B). As we have shown, formation of the transition state is compromised upon addition of osmolytes. This sensitivity to water content suggests R67 DHFR is not very efficient in vivo, supporting our previous hypothesis that R67 DHFR is a primitive enzyme (4).
Conclusion—R67 DHFR is an unusual enzyme that uses symmetry-related sites to bind two different ligands. NADPH interacts with R67 DHFR through specific contacts, and a unique bound conformation exists (8). This structural description is consistent with results from this study, where removing surface water by osmotic stress results in tighter binding of NADPH. In contrast, DHF utilizes specific contacts for binding of its pteridine ring (8); however, its para-aminobenzoylglutamate tail of dihydrofolate/folate is disordered. This different mode of binding may correlate with its different response to osmolyte addition. Because NADPH serves as an internal control, whatever effect causes the different behavior associated with DHF binding must lie in either the free ligand, free protein, or in the different mode of binding in the complexes. Clearly interfacial water contributes to binding plasticity in R67 DHFR and can provide hydrogen bonds either alone or in networks between ligand and protein (93–97). In other words, these osmolyte studies indicate that NADPH and R67 DHFR interact more directly, utilizing more protein contacts (i.e. a "dry interface" according to Janin (98)), whereas the DHF interaction with protein uses more water contacts ("wet interface"). This view proposes water is intimately involved in binding specificity as well as affinity in R67 DHFR.
| FOOTNOTES |
|---|
The on-line version of this article (available at http://www.jbc.org) contains supplemental text, equations, references, Figs. 1–4, and Table S1. ![]()
1 To whom correspondence should be addressed: Dept. of Biochemistry, Cellular and Molecular Biology, University of Tennessee, Knoxville, TN 37996-0840. Tel.: 865-974-4507; Fax: 865-974-6306; E-mail: lzh{at}utk.edu.
2 The abbreviations used are: DHFR, dihydrofolate reductase; DHF, dihydrofolate; ITC, isothermal titration calorimetry; MES, 4-morpholineethanesulfonic acid; PEG, polyethylene glycol; TMAO, trimethylamine N-oxide; TMP, trimethoprim. ![]()
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| All ASBMB Journals | Molecular and Cellular Proteomics |
| Journal of Lipid Research | ASBMB Today |