Previous Article | Next Article ![]()
Journal of Bacteriology, September 2003, p. 5519-5526, Vol. 185, No. 18
0021-9193/03/$08.00+0 DOI: 10.1128/JB.185.18.5519-5526.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
J. M. Wiame Research Institute, Microbiology, Free University of Brussels (VUB), B-1070 Brussels,1 Laboratory of Biochemistry, Institute of Chemistry B6, University of Liège, Liège, Belgium2
Received 6 January 2003/ Accepted 23 June 2003
|
|
|---|
|
|
|---|
At the other end of the temperature scale, psychrophilic (24) enzymes face the necessity of being efficient catalysts at a low energy cost (4, 9, 14). On the basis of a much more limited set of data than for thermophilic enzymes (see recent reviews in references 33, 35, and 45), two basic strategies for psychrophily were distinguished: either an overall increase of flexibility generally associated with a marked decrease of stability, or a more subtle strategy which keeps one domain rigid enough to provide enough affinity for the substrate while holding the catalytic domain in a more flexible state (2, 21, 45). In the former case thermodynamic and genetic data even suggest that some enzymes, such as the
-amylase from a psychrophilic bacterium (10), have achieved the maximum flexibility compatible with the maintenance of a correctly folded structure. This in turn indicates that the lowest temperature compatible with active life (by opposition to mere survival) may not be determined by the freezing point of the aqueous environment only. The challenge of stability at low temperature is however not the only one faced by cold-active enzymes. Kinetic optimization is another one, especially for intracellular biosynthetic enzymes since (i) improving kcat by decreasing the activation enthalpy may have a cost in terms of affinity due to a concomitant increase in conformational entropy (11, 21); (ii) such enzymes usually have to cope with low substrate concentrations. Investigations on the catalytic efficiency of cold-active enzymes acting at key metabolic steps are thus important.
In this work we report the characterization of dihydrofolate reductase (DHFR) (EC 1.5.1.3) from Moritella profunda (DHFRMp), a novel, strictly psychrophilic and moderately piezophilic species isolated from the deep sea (44). To our knowledge this is the first report on a cold-active DHFR.
DHFR catalyzes the NADPH-linked reduction of 7,8-dihydrofolate (DHF) to 5,6,7,8-tetrahydrofolate (THF). It is a metabolically and clinically important enzyme since it is involved in over 20 different one-carbon transfer reactions and is the target of several antibacterial, antiprotozoal and anticancer drugs. It is also an interesting model system because of its small size, thus amenable to nuclear magnetic resonance studies (with a molecular mass of about 19 kDa when monomeric).
|
|
|---|
Cloning of the DHFR gene.
Restriction enzymes and T4 ligase were purchased from Boehringer Mannheim and used according to the manufacturer's instructions. M. profunda genomic DNA partially digested with enzyme Sau3A was used to construct a genomic
ZAP DNA library in pBK-CMV (Stratagene) vector (43). This genomic library was used for transformation of E. coli strain XL1-Blue MRF (Stratagene) and positive clones were selected for resistance to 5 mM trimethoprim (a DHFR inhibitor) at 30°C on plates of medium 853. The corresponding genomic DNA was retrieved by the polymerase-chain reaction and sequenced to check that no mutant gene had been isolated.
Overexpression and purification of recombinant DHFRMp. The DHFR gene of M. profunda was amplified by PCR to produce NdeI and BamHI restriction sites at the ends of the fragment. The NdeI restriction site was designed to overlap the ATG start codon. The fragment was cloned into the pET24a vector (Novagen) which was used to transform E. coli strains and the nucleotide sequence was checked. Exponentially growing cells of E. coli BL21-codonPlus(DE3)-RIL (Stratagene) bearing the plasmid were induced with 1 mM isopropyl-1-thio-ß-D-galactopyranoside (IPTG) for 20 h at 16°C, harvested by centrifugation and washed in 0.9% NaCl. Cells were disrupted with a high-pressure cell disruptor (basic Z model; Constant Systems Ltd.) at a pressure of 28,000 lb/in2 in 50 mM K phosphate buffer, pH 7.0. The extracts were freed of cell debris by centrifugation, first for 30 min at 12,000 x g and then for 1 h at 130,000 x g. The supernatant was applied to a methotrexate-agarose (Sigma) column (XK26/20; Pharmacia) equilibrated with buffer A (50 mM K phosphate buffer, 3 mM 2-mercaptoethanol, pH 7.5). The column was then washed three times until no more protein was detected at 280 nm in the effluent, first with the same buffer, then with buffer A containing 2.3 M KCl, and finally with buffer B (50 mM Tris HCl, 1 M KCl, 3 mM 2-mercaptoethanol, 5% glycerol [vol/vol], pH 8.5). The enzyme was then eluted using 15 ml buffer B containing 15 mM folate at pH 9.0 followed by buffer B until no absorbance at 280 nm was detected in the effluent. Fractions containing the eluted peak were combined and concentrated by ultrafiltration (YM-3 membrane; Amicon) until DHFR activity began to be detected in the flowthrough. Concentrated enzyme was applied to a Hiload 16/60 Superdex 75 prep-grade column equilibrated with 20 mM phosphate buffer, 0.15 M NaCl, 6 mM 2-mercaptoethanol, 5% glycerol (vol/vol), pH 6.5 and calibrated using the Pharmacia molecular weight standards. The activity of the eluted fractions was tested and the fractions pooled.
Fluorescence.
Fluorescence of DHFR samples was recorded using an SML-AMINCO Model 8100 spectrofluorimeter (Spectronic Instruments) at an excitation wavelength of 280 nm and at an emission wavelength of 330 nm (for heat denaturation) or 340 nm (for urea denaturation). In order to reduce photobleaching of the samples, the excitation band-pass was kept at 1 nm (4 nm for emission) and was closed between records. Samples (1 ml) in 20 mM K phosphate buffer at pH 7.0 were adjusted to an A280 of
0.1 and were heated in the fluorometer cell at a rate of 2°C/min using a programmed water bath Lauda Ecoline RE306. Temperature of the sample was recorded in the cell by a Cooper-Constantan thermocouple connected to a calibrated digital thermometer (Thermo Electric). Unfolding by urea was monitored at 15°C after 1 h of incubation of the samples at this temperature. Data were normalized using the pre- and posttransition baseline slopes as described (26). In the linear portion of the transition between the native and unfolded states, the equilibrium constant, K and the free-energy change, G can be calculated using the relations K = fraction unfolded/fraction folded and
G = -RT lnK. Least-squares analysis of G values as a function of urea concentrations allowed to estimate the conformational stability in the absence of denaturant,
GH2O according to the following equation:
G =
GH2O - m [urea].
DSC.
Differential scanning calorimetry (DSC) was performed using a MicroCal MCS-DSC instrument at a scan rate of 2°C/min and under a nitrogen pressure of 2 atm. Samples were desalted by gel filtration through a PD10 column (Pharmacia) equilibrated with 20 mM K phosphate buffer, pH 7.0. This buffer was used in the reference cell and for buffer baseline determination. Thermograms of DHFR were analyzed according to a single non-two-state transition model in which Tm,
Hcal, and
Heff were fitted independently using the MicroCal Origin software (version 2.9).
Enzyme assays. The standard reaction mixture contained 50 mM K phosphate buffer at pH 7.0, 1 mM 2-mercaptoethanol (or 10 mM 2-mercaptoethanol for saturation studies), 0.12 mM NADPH (Roche), 0.05 mM DHF (Sigma) and about 30 ng enzyme. Enzyme was incubated with NADPH, the reaction initiated by addition of DHF and monitored by recording absorption decrease at 340 nM using a PU8740 UV/VIS scanning spectrophotometer (Philips), in a temperature-controlled cuvette. A unit of enzyme is defined as the amount that will reduce 1 µmol of substrate per min, based on a molar extinction coefficient of 12,000 (17). Protein concentrations were determined by the method of Lowry.
In order to record thermal inactivation, DHFR (30 µg/ml) was incubated in 20 mM phosphate buffer (pH 7.0) in the absence of substrates for 15 min at different temperatures and then kept at 0°C. Residual activity was measured in the same buffer under standard conditions at 30°C.
Half-life of activity was determined by incubating enzyme (30 µg/ml) in 20 mM phosphate buffer at pH 7.0 with or without 0.5 mM 2-mercaptoethanol and 0.5% glycerol (vol/vol) at a fixed temperature. Residual activity was measured in the same buffer under standard conditions at 30°C.
The kinetic parameters kcat and Km determined by the initial velocity method using a nonlinear regression fit of the saturation curves on the Michaelis-Menten equation performed by SigmaPlot 4.0. The IC50s (inhibitor concentration required for 50% inhibition of the activity) were obtained by fitting the inhibition data on an nonexponential decay.
|
|
|---|
ZAP library containing M. profunda genomic DNA was used to transform E. coli XL1-Blue recipient cells, selecting for resistance to trimethoprim as described in Materials and Methods. The copy number of the recombinant DHFR gene was expected to be sufficient to produce this phenotype since overexpressing (25- to 30-fold) the wild type E. coli DHFR gene in such a way confers trimethoprim resistance to the host (31). One such clone, found to contain an insert of about 2.7 kb, was analyzed further. Sequencing about 1,000 bp from this fragment revealed two open reading frames homologous to the E. coli DHFR and proline permease genes, respectively. We checked that the cloned DHFR gene was identical to the native one after retrieving the cognate genomic DNA by the polymerase chain-reaction. Enzymological studies (see "Kinetic parameters and temperature dependence") showed that DHFRMp is actually less sensitive to trimethoprim than E. coli DHFR; this may have contributed to the success of the selection procedure. DHFR sequence data have been submitted to the EMBL database under accession number AJ487535.
The DHFRMp open reading frame encodes 162 amino acids and is clearly homologous to all known DHFRs despite a low general level of identity. It presents the highest score (more than 55% identity) with E. coli and Salmonella enterica serovar Typhimurium DHFRs and relatively high values are still obtained with other Proteobacteria belonging to the
-3 group (Fig. 1). Codons with A in the third position, such as AGA, AUA, and CUA, are rare in E. coli but not in Moritella (43). This explains the necessity of using an ad hoc host strain for overexpression (see below).
![]() View larger version (73K): [in a new window] |
FIG. 1. Comparison of DHFR primary structures: alignment of amino acid sequences from M. profunda (MORPR), S. enterica serovar Typhimurium (SALTY), E. coli (ECOLI), Citrobacter freundii (CITFR), Enterobacter aerogenes (ENTAE), Klebsiella aerogenes (KLEAE), and Haemophilus influenzae (HAEIN) (all above 50% identity) and from T. maritima (TEMA). DHFRMp K20 (A19 in E. coli) and E27 (D26 in E. coli) are in bold and underlined. A dotted line above the alignment indicates the extent of the M20 loop (residues 10 to 24) as in E. coli DHFR. Secondary structures, as defined for E. coli DHFR (3) are indicated above the alignment ( -helices and ß-strands). Asterisks indicate identical amino acids and dots indicate similar amino acids conserved in all eight proteins.
|
Expression, purification, and molecular mass of DHFRMp. The DHFRMp gene was inserted into expression vector pET24a, which was used to transform E. coli strain BL21(DE3) (Novagen). However, the expression level of the recombinant protein was low. Since rare codons could have hampered expression (as the sequence analysis confirmed), BL21-codonPlus(DE3)-RIL (Stratagene) competent cells were used. After IPTG induction the expression level was so high that the affinity chromatography step followed by gel filtration yielded as much as 300 mg of pure DHFR (as shown by SDS and native PAGE) from 1 liter of culture. The enzyme was eluted as the major peak from the calibrated gel filtration column at about 19 kDa. Mass spectrometry indicated a molecular mass of 18,290 Da (courtesy of V. Villeret), thus in the expected range for a monomer and in excellent agreement with the value deduced from the amino acid sequence derived from the gene (18,291.8 Da).
Resistance to thermodenaturation. As shown in Fig. 2, DHFRMp withstood incubation up to about 35°C for 15 min but became inactivated at higher temperatures, illustrating its psychrophilic character. The half-life of activity was 25 min at 45°C but in the presence of 0.5 mM 2-mercaptoethanol and 0.5% glycerol no activity loss was observed and the half-live at 50°C was 18 min. Heat-induced unfolding of DHFRMp recorded by fluorescence (Fig. 3) also demonstrated that addition of glycerol and 2-mercaptoethanol shifted the midpoint of the thermal transition upwards by about 25°C. This pronounced stabilization obtained by combining a solvophobic effect and reducing conditions allows proper storage and handling of this relatively unstable enzyme.
![]() View larger version (14K): [in a new window] |
FIG. 2. Effect of temperature on DHFRMp activity and stability. A: residual activity after 15 min incubation at the temperature indicated; assay carried out at 30°C. B: thermodependence of the activity; assay carried out for 5 min at the temperature indicated.
|
![]() View larger version (16K): [in a new window] |
FIG. 3. Stabilization of the heat-labile DHFRMp: heat-induced denaturation was performed in phosphate buffer (open symbols) and in phosphate buffer containing 0.1 mM. 2-mercaptoethanol and 0.1% glycerol (closed symbols).
|
|
View this table: [in a new window] |
TABLE 1. Kinetic parameters of psychrophilic, mesophilic, and thermophilic DHFRs
|
Conformational stability.
In order to compare the structural stability of DHFRMp with that of its mesophilic and thermophilic homologues from E. coli and T. maritima, the reversible urea-induced unfolding curves were recorded by fluorescence at 15°C. As shown in Fig. 4 and Table 2, moderate concentrations of urea induced denaturation of the psychrophilic enzyme, with a concentration for half-denaturation of about 1.6 M. The m value of DHFRMp, which is a measure of the dependence of
G on denaturant concentration, is similar to that of E. coli DHFR. By contrast, linear extrapolation of the
G values for the main transition (Fig. 4) indicates that the conformational stability in the absence of denaturant (
GH2O) is much lower than that of both mesophilic and thermophilic homologs.
![]() View larger version (19K): [in a new window] |
FIG. 4. Urea-induced unfolding curves of DHFR. The unfolded fraction of DHFRMp is plotted for increasing urea concentrations (closed symbols). The sigmoidal transitions of E. coli DHFR (28) and of T. maritima DHFR (6) are included for comparison. The inset shows linear extrapolation of G values in the main transition region of DHFRMp to zero denaturant concentration.
|
|
View this table: [in a new window] |
TABLE 2. Thermodynamic parameters of urea-induced unfolding of DHFRs at 15°C
|
![]() |
Hcal, corresponding to the total amount of heat absorbed during unfolding, and the Tm value were obtained from thermograms of DHFRMp (Fig. 5). The
Cp value, which determine the plot curvature, was taken from data for proteins of similar size (24). Figure 5 illustrates the resulting stability curve of DHFRMp with
Cp values ranging from 1.5 to 2.0 kcal mol-1 K-1, as well as the profiles of the mesophilic E. coli DHFR (5, 19) and of the thermophilic T. maritima DHFR in the presence of 2.9 M guanidinium chloride (6). These stability curves represent the energy required to disrupt the native state (30). The increase in stability in the order (from most to least stable) psychrophile, mesophile, and thermophile is clearly illustrated by the progressively higher
Gmax (top of the curve) and Tm (for
G = 0). Although
Cp has not been determined experimentally for DHFRMp, this parameter is not prone to drastic variations between protein homologues (23) but mainly affects the predicted temperature of cold unfolding at subzero temperatures.
![]() View larger version (19K): [in a new window] |
FIG. 5. Stability profiles of DHFRs. The Gibbs free energy of denaturation for DHFRMp was calculated by equation 1 with Hcal = 50 kcal mol-1, Tm = 45.6°C, and Cp = 1.5 kcal mol-1 K-1(upper curve) or 2.0 kcal mol-1 K-1(lower curve). The stability curves for E. coli DHFR are from reference 19 for the lower curve and from reference 5 for the upper curve. Data for T. maritima DHFR in the presence of 2.9 M guanidinium chloride are from reference 6. The inset shows a differential scanning calorimetry thermogram of DHFRMp.
|
|
|
|---|
|
View this table: [in a new window] |
TABLE 3. Amino acid content in DHFRs from M. profunda, E. coli, and T. maritima
|
G than E. coli DHFR but at much lower temperatures), by increasing the steepness of the stability curve (with similar maximal value of
G than E. coli DHFR at the same temperature but with a lower Tm), or by a global collapse of the stability curve. Each of these alternatives corresponds to a combination of alterations in the
Hcal, Tm, or
Cp values in equation 1. Obviously, the stability curve of DHFRMp follows the third alternative: the weak conformational energy arises from low values of both
Hcal and Tm. This demonstrates that the native state of DHFRMp is stabilized by a lower number of enthalpically driven interactions when compared with its mesophilic and thermophilic homologues. This is in agreement with the above-mentioned lower number of charged residues in DHFRMp. By contrast, the stability profile of T. maritima DHFR is shifted to higher temperatures with a sharply increased
Gmax value (6). This indicates that the weak interactions responsible for the enhanced stability of the thermophilic enzyme are strengthened at high temperatures, suggesting the involvement of nonpolar groups (25). The relatively low number of hydrophobic amino acids in DHFRMp in comparison with T. maritima DHFR is also in agreement with this aspect.
Two further features revealed by previous work on a psychrophilic
-amylase (10) are also confirmed for DHFRMp. First the cold-adapted enzyme is not cold-stable as judged by the stabilization energy value at 0°C and by the extrapolated temperature of cold-unfolding at a subzero temperature. Secondly, the physiological temperature for DHFRMp lies on the left side of the bell-shaped stability curve, in sharp contrast with E. coli and T. maritima DHFRs. Accordingly, the thermal dissipative force is a main determinant of the structural mobility for the latter enzymes. By contrast, group hydration responsible for the decrease of
G values at low temperatures (23) significantly contributes to the conformational flexibility of DHFRMp which is required for enzyme activity around 0°C. It will be interesting to determine whether DHFRMp, like the above-mentioned
-amylase, approaches the maximal flexibility still compatible with stable folding. Indeed, the temperature of cold denaturation of some critical enzymes, and not the freezing point of the aqueous environment, may actually set the lower limit of bacterial development.
E. coli DHFR catalyzes reaction cycles through five detectable kinetic intermediates which have been depicted by crystallographic studies (3, 34). It was shown that the M20 loop (from residues 10 to 24) is predominantly closed in the E. coli holoenzyme. This loop moves over and between the reactants in the Michaelis and transition state complexes but during the remaining steps, the loop protrudes into the nicotinamide-ribose binding pocket. Thus, this loop is expected to assist the formation of the catalytic intermediates by driving the substrates into the active site. When compared to other DHFRs, the psychrophilic enzyme shows extensive conservation of catalytically important residues but differs by the A19K substitution in the M20 loop (Fig. 1). As a consequence, the basic side chain of Lys19 could interact with the oppositely charged folate during the reaction and could be responsible for the high kcat values. It should be mentioned that improved electrostatic fields at the active site of some psychrophilic enzymes seem involved in the adaptation to catalysis at low temperatures (35).
The apparent Km value for DHF of DHFRMp (Table 2) is more than 50 times the value currently reported by several authors for the cognate enzyme in E. coli (27, 29, 36) and is also considerably higher than for L. casei (8) and T. maritima (41). Detailed examination of published data does not suggest that assay conditions influenced actual values to an extent in anyway comparable to the differences we observe. Moreover, the sensitivity toward the inhibitors methotrexate and trimethoprim is much lower than for E. coli DHFR, indicating that the affinity of the enzyme has also been affected. In this respect, the specific replacement of Asp27 in the folate binding site by Glu in DHFRMp (Fig. 1) can possibly decrease the affinity as a result of steric hindrances arising from the larger acidic side chain. As regards NADPH, the Kmapp of DHFRMp is higher than in other homologues but remains in the same order of magnitude. Because these apparent Km values are not substantially different between 5°C and 30°C, the catalytic efficiency (kcat/Km) is definitely lower for DHFRMp than for E. coli DHFR, especially with respect to DHF. Therefore, DHFRMp appears to have adapted to low temperatures by maintaining a relatively high kcat at the expense of its affinity for one of the substrates of the reaction. The fact that this trade-off involves DHF rather than NADPH could be understood by comparison with E. coli DHFR since (i) efficient binding of NADPH appears necessary to promote the release of THF, product of the reaction, and (ii) the enzyme must keep a higher affinity for NADPH than for NADP because the latter is less efficient in THF release than NADPH and thus could slow down the reaction (34). Maintenance of a sufficient affinity for NADPH in DHFRMp was probably selected during the course of evolution.
From the kinetic point of view, DHFRMp thus appears to remain well below optimization, supporting the concept that complete metabolic adaptation to cold may remain out of reach. We have obtained similar data for ornithine carbamoyltransferase from a closely related species of Moritella (42). Since M. profunda is moderately piezophilic it will be interesting to determine whether increased hydrostatic pressure has any effect on the catalytic efficiency of the enzyme. Moreover, since the data presented in this work have been obtained with pure DHFR, one can wonder whether, in vivo, the relatively high Km is not partly alleviated by compartmentalization effects such as obtained in metabolic channeling.
We thank V. Villeret (Institut de Biologie de Lille, France) for helpful assistance in structural interpretation and for the molecular mass determination by mass spectrometry; we thank D. Georlette for her kind help.
|
|
|---|
-amylase from the Antarctic bacterium Alteromonas haloplanctis. Biochemistry 38:4613-4619.[CrossRef][Medline]
-acetyl ornithinase. J. Bacteriol. 182:1609-1615.
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»