Previous Article | Next Article ![]()
Applied and Environmental Microbiology, May 2007, p. 3165-3172, Vol. 73, No. 10
0099-2240/07/$08.00+0 doi:10.1128/AEM.02960-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Field of Microbiology, Cornell University, Ithaca, New York 14850,1 Department of Molecular Biology and Genetics, Cornell University, Ithaca, New York 148532
Received 20 December 2006/ Accepted 10 March 2007
|
|
|---|
|
|
|---|
Thermobifida fusca is a filamentous soil bacterium that grows at 50°C in defined medium and can utilize cellulose as its sole carbon source. It is a major degrader of plant cell walls in heated organic materials such as compost piles and rotting hay. Six different cellulases have been purified and characterized in our lab using standard purification methods. Among these six cellulases, three of them are endocellulases, Cel9B, Cel6A, and Cel5A (7, 8); two of them are exocellulases, Cel6B and Cel48A (6, 20); and one is a processive endocellulase, Cel9A (5, 7).
Thermobifida fusca Cel9A-90 consists of a family 9 catalytic domain (CD), a family 3c cellulose binding module (CBM3c), a fibronectin III-like domain, and a family 2 CBM. It has the highest activity of any individual T. fusca enzyme on crystalline substrates. Cel9A has relatively high activity on carboxymethyl cellulose (CMC), which is a feature of an endocellulase. However, it produces 87% soluble products from filter paper (FP), a form of crystalline cellulose, and only 13% insoluble products (19). True endocellulases produce 30 to 40% insoluble products from FP, and exocellulases produce only 5 to 8% insoluble products; thus, Cel9A is quite processive (19). Another interesting feature of Cel9A is that it can synergize with both reducing and nonreducing end-directed exocellulases and with endocellulases, while true endocellulases do not synergize with each other. This evidence confirms that Cel9A possesses the properties of both endo- and exocellulases (19). Cel9A-68 consists of the family 9 CD rigidly attached to the CBM3c by a short linker. The crystal structure of this enzyme was determined by X-ray crystallography at 1.9-Å resolution (Protein Data Bank code 4TF4) (13). The CD has a shallow open active site cleft with six known glucose binding sites, numbered from 4 to +2, aligned with a flat binding surface of the CBM3c (Fig. 1). The linker region between the CD and this CBM (G445 to P461) lies along the opposite side of the CD from the active cleft with many contacts with the CD and the CBM. The CBM consists of a 10ß-stranded sandwich motif (13). It binds weakly to crystalline cellulose and is thought to play a role in catalysis by feeding a single cellulose molecule into the catalytic cleft. The CBM is critical for hydrolysis and for processivity, as a cloned enzyme construction of the family 9 CD with the linker Cel9A-51 but without this CBM had very low activity on CMC, bacterial cellulose (BC), and swollen cellulose (SWC) and showed no processivity (4) (Table 1).
![]() View larger version (71K): [in a new window] |
FIG. 1. Space-filled structure of Cel9A-68 (4TF4) with six glucose molecules in the catalytic cleft. This enzyme product structure was obtained after soaking the crystals in cellopentaose (G5). G5 was cleaved to G4 and G1 or G3 and G2, and density for six glucose residues can be seen in the structure because G5 can bind from the position of Glc(3) to Glc(+2) or from Glc(4) to Glc(+1). The CBM residues chosen for mutation are shown in dark gray.
|
|
View this table: [in a new window] |
TABLE 1. Activity and binding characteristics of mutant enzymes
|
|
|
|---|
Mutagenesis and DNA manipulations.
Standard nucleic acid techniques were used as described by Sambrook et al. (14). All mutant genes were generated using the QuikChange site-directed mutagenesis method (Stratagene). Two oligonucleotide primers for each mutation were designed, complementary to each other on opposite strands at the site of the desired mutation. The length of both primers was 40 to 45 bases; G or C was on each terminus, and one restriction enzyme site was created or deleted for screening. QuikChange PCR was performed under the following conditions: (i) denaturation at 95°C for 2 min; (ii) addition of Pfu polymerase (Stratagene); (iii) 18 cycles of denaturation at 95°C for 60 s, annealing at 60°C for 50 s, and extension at 68°C for 8 min; and (iv) a final step at 68°C for 7 min. The sample (25 µl) was then incubated with 20 units of DpnI for 3 h at 37°C to remove the methylated template, and the sample was concentrated by ethanol precipitation at 20°C overnight. The PCR product was transformed into E. coli XL10, and plasmid minipreparations of individual colonies were screened by restriction enzyme digestion for the altered site. The DNA sequences of the candidate plasmids were determined by the Cornell Biotechnology Resource Center using a Perkin-Elmer/Applied Biosystems automated sequencer to confirm the presence of the desired mutation and to confirm that no other mutations were present. Plasmid DNA was transformed into E. coli BL21(DE3) for protein expression.
Enzyme production and purification.
E. coli BL21(DE3) strains were grown at 37°C overnight in 30 ml of Luria broth with 60-µg/ml kanamycin. The overnight culture was transferred into 1 liter of M9 medium containing 0.8% glucose and 60-µg/ml kanamycin. After growth at 30°C to an optical density at 600 nm (OD600) of about 0.8, isopropylthio-ß-D-galactoside was added to the culture to 0.8 mM and growth was continued at 30°C overnight. The proteins were purified from the culture supernatants on a phenyl-Sepharose column, followed by a Q-Sepharose column, as described in previous work (4). Protein purity was assessed on sodium dodecyl sulfate gels, and the protein concentrations were measured by absorbance at OD280 using molar extinction coefficients of 165,480 for Cel9A-68 and 129,110 for Cel9A-51, which were determined from the predicted amino acid compositions using the DNA Star Protean program.
Activity assays.
The activities of Cel9A-68 wild-type and mutant enzymes were determined concurrently on CMC (10 mg/ml), BC (2.5 mg/ml), phosphoric acid-SWC (2.5 mg/ml), and Whatman FP no. 1 (8 mg/ml). All assays were run in a reaction volume of 0.4 ml in triplicate at 50°C for 16 h in 0.05 M sodium acetate, 0.015 M CaCl2 (pH 5.5), using from 0.5 to 600 pmol of protein as needed. Reducing sugars were measured using the dinitrosalicylic acid (DNS) reagent (3) and a glucose standard curve using Sigma 635-100 glucose standard solution. Thin-layer chromatography of the BC assay products under these conditions showed that the major product was cellobiose along with approximately 10% glucose and a small amount of cellotriose for some mutants (data not shown). To calculate the activity, the OD600 values were converted to micromoles of cellobiose by comparing a cellobiose standard curve and a glucose standard curve. A major advantage of the DNS reagent is that it does not have a blank with the amount of enzyme used and it fits the assay range well. The disadvantage of the DNS reagent, used with a cellobiose standard, is that glucose will be about 30% underestimated and cellotriose will be slightly overestimated. Mutant enzyme activities were determined concurrently with those of the wild type; the nanomoles of protein used (X) were plotted versus the OD600 (Y) and fitted to the equation Y = m1X/(m2 + X) using the program KaleidaGraph from Synergy Software. This curve was used to determine the amount of enzyme required to digest the substrate to the target value of 5% for CMC and FP (0.585 µmol and 0.467 µmol cellobiose, respectively) and 10% for SWC and BC (0.292 µmol cellobiose). For those enzymes which could not reach the target digestion, the activity was calculated from the micromoles of cellobiose produced by 0.6 nmol of enzyme. Although the estimated overall error of CMC and FP activities is around 5% and that for SWC and BC is around 10% from previous work, the calculated coefficients of variation for a given assay are below 5% as noted in Table 1.
Processivity was determined from the distribution of reducing ends between the FP (insoluble) and the supernatant (soluble) following a published procedure (5). Assays were performed on FP using 0.6 nmol of enzyme per reaction; the supernatant was separated from the FP circle, and the reducing sugar content of each fraction was determined with the DNS reagent (3).
Binding assays.
Binding to BC was measured in 2.2-ml "low-adhesion" tubes (USA Scientific) by adding 2 nmol of enzyme to 2.5 mg/ml BC in a total volume of 0.5 ml of 50 mM sodium acetate buffer, pH 5.5, plus 10% glycerol. A rack of tubes was placed on a Nutator rocking table (Clay-Adams) at 4°C for 1 h. The tubes were then centrifuged for 5 min at 13,000 rpm, and the top 300 µl of liquid was transferred to 1.5-ml Eppendorf tubes with a gel loading pipette tip. The tubes were centrifuged again for 5 min at 13,000 rpm, and the top 150 µl of liquid was transferred to clean 0.5-ml Eppendorf tubes with a gel loading pipette tip. The amount of unbound protein was determined by the protein concentration calculated from the OD280; control reactions without BC and also with BC alone were performed for every assay.
Ligand binding studies.
Dissociation constants for the binding of 4-methylumbelliferyl ß-cellotrioside (MUG3) to Cel9A-68 and to mutant enzymes were determined by direct fluorescence quenching titrations using an Aminco SLM 8000C spectrofluorimeter and were calculated as described previously (1). The initial concentration of MUG3 was 1.8 µM in 903 µl of 50 mM sodium acetate buffer, pH 5.5, and 100 µM protein was added at 3.5 µl/min. The temperature was maintained at 6.5°C, and excitation was at 316 nm, while emission was measured at 360 nm.
Azide rescue studies.
The azide activity rescue studies were run in a reaction volume of 0.4 ml in triplicate at 50°C for 16 h in 0.05 M sodium acetate (pH 5.5), with 0.4 nmol of protein and addition of 0.075 M to 0.6 M sodium azide; 200 µl of 5.0-mg/ml SWC was used as substrate in the reactions. Reducing sugars were measured using the DNS method, and the activity was determined from the micromoles of cellobiose produced by 0.4 nmol of enzyme.
Accession numbers.
The sequence files for Cel9A (formerly called E4) are L20093 and P26221. The structure files for Cel9A-68 are 1JS4, 1TF4, 3TF4, and 4TF4.
|
|
|---|
Activity of CD mutant enzymes.
The selection of CD target residues focused on the catalytic site, especially those residues with hydrogen bonding or hydrophobic stacking to one of the six sugar residues in the active cleft. The residues chosen are listed in Table 1 and illustrated in Fig. 2. Data for some mutated residues previously published but discussed in this paper are included in Table 1 and referenced.
![]() View larger version (26K): [in a new window] |
FIG. 2. Cel9A-68 product complex structure (4TF4) showing amino acids chosen for mutation around the catalytic center in light gray and the six sugar residues in dark gray. Hydrogen bonds are shown as dashed lines. Cleavage occurs between Glc(1) and Glc(+1).
|
![]() View larger version (22K): [in a new window] |
FIG. 3. Sodium azide rescue test. The samples were treated with various concentrations of azide (0 M to 0.6 M) for 16 h, and then SWC activity was assayed and calculated [µmol of cellobiose/(min·µmol of enzyme)]. The average coefficients of variation were 4%, 5%, 3%, and 4% for wild-type and E424, D55, and D58 mutant enzymes, respectively.
|
![]() View larger version (26K): [in a new window] |
FIG. 4. Enzyme product complex (4TF4) showing the hydrogen bond network between Y206, D55, Glc(1), D58, and H125. The hydrogen bonds from the catalytic acid E424 to Glc(1) and Glc(+1) are also shown.
|
F205 and Y429 (Fig. 2) do not have direct interactions with any sugar residues, even though they are adjacent to the catalytic site. The F205A mutant enzyme has lower activity on all substrates, lower processivity, and significantly lower oligosaccharide binding. F205 has one hydrogen bond from its backbone carbonyl group shared between Glc(2) O6 and Glc(1) O3, which is consistent with lower binding. The Y429A mutant has normal binding but low crystalline cellulose activity and impaired processivity.
W209 appears to stack against Glc(3), and W256 stacks against Glc(4) and is hydrogen bonded to the Glc(3) O6 (Fig. 2). D261 is at the bottom of the substrate cleft and is hydrogen bonded to Glc(4) O2 and to S210 OG. D261 and S210 work together to bind Ca2+. R317 is hydrogen bonded to Glc(2) O2, and a previously studied residue, Y318, appears to be hydrogen bonded to Glc(2) O3, Glc(3) O5, or Glc(3) O6. All of the enzymes with mutations of the above residues, Y318A (21), Y318F (21), W209S, W256A, R317K, D261A, and D261N, have almost-normal activity on SWC and severalfold-higher activity on CMC, yet only 15% of wild-type activity on BC and reduced processivity. All of these residues are part of the Glc(2) to Glc(4) subsites, and it appears that each of them participates in binding crystalline substrates such as BC in the active cleft. The higher CMC activities could be due to the smaller side chains of the mutants allowing more room for the carboxymethyl groups to enter and leave the catalytic cleft.
Another interesting mutation is R378K, which makes two hydrogen bonds, NH1 and NH2, to Glc(+1) O2 (Fig. 2). The mutation does not change enzyme activity very much, but processivity is improved significantly. So far this is the only mutation in the CD that results in obviously higher processivity. The D261A mutant enzyme has lower processivity than the wild type, and it is located near Glc(4). An D261A/R378K double mutant enzyme was made, and its processivity decreased dramatically, even lower than that of the D261A mutant enzyme, indicating that processivity requires a precise balance between catalytic cleft binding on the two sides of the cleavage site.
The length of one cellobiose unit is about 10.3 Å (17), and the distance from the top of the catalytic cleft to the end of the CBM3c is about 70 to
80 Å, so mathematically as many as 16 glucose units could span the CD and the CBM. Amino acids which interact with sugar residues other than those in the six known subsites could play important roles in activity and processivity. It is apparent from the structure that W128 (Fig. 2) could stack or interact with a possible Glc(+3) subsite. The W128A mutant enzyme lost activity on all substrates, and processivity dropped significantly. It is also interesting that this tryptophan side chain is predicted to have several energetically favorable positions based on the rotamer scores (three out of seven positions have scores of 1) reported by Swiss-PdbViewer.
Activity of CBM3c mutant enzymes.
The initial step in cellulose degradation is the binding of the enzyme to the substrate, and this process clearly involves the CBM3c although Cel9A-68 binds only weakly to BC. There is no direct evidence for which CBM residues bind cellulose, but Tormo et al. (18) hypothesized that the flat face of the family 3a CBM, which is made up of four ß-strands, could interact with three chains of cellulose: a strip of mostly aromatic residues, D, H, Y, R, and W, could interact with one chain; two well-conserved polar residues, N and Q, could interact with the middle chain; and nonconserved polar residues could interact with the third chain. The Cel9A CBM3c does not contain the strip of aromatic residues which exists in family 3a and 3b CBMs, even though it has the flat face which is made up of strands 1, 2, 4, and 7 (13). Six residues in the Cel9A-68 CBM that could potentially interact with a cellulose molecule were chosen for mutation based on sequence alignment (Fig. 5), a model (2), and consideration of the structure (Table 1). BC activity was lowered significantly in the R557A/E559A, Y520A, and R563A mutant enzymes. These amino acids are "in line" with the catalytic cleft (Fig. 1). SWC activity was only slightly lowered by these mutations, and CMC activity was not significantly affected. In spite of the lower BC activity, processivity was not decreased, which is an unexpected result. Mutation of the highly conserved F476 residue to A resulted in 21% decreased BC activity but increased processivity. The D513A mutation increased BC activity by 21%. D513 is located at the end of the CBM3c and appears to be not quite aligned with the catalytic cleft. Another interesting residue is I514, which is conserved neither in family 3c nor on the ß strands of the flat face. However, sequence alignment shows that this position is conserved, as an H in family 3a and a W in family 3b. In a family 3a CBM from Clostridium cellulovorans, the H is in a planar strip and plays an important role in cellulose binding (9). The I514H mutant enzyme has higher activity on BC and almost-wild-type activity on the other substrates. The BC binding of all of the tested CBM3c mutant enzymes was as good as or better than that of the wild type.
![]() View larger version (93K): [in a new window] |
FIG. 5. Structure-based sequence alignment of family 3 CBMs. Secondary structural ß-strands are indicated by arrows and enumerated. Proposed cellulose binding residues are shown in dashed boxes, and for T. fusca Cel9A-68, they are F476, D513, Y520, Q561, and R563. Residues shown in solid boxes are conserved, and those T. fusca residues are V465, Y501, Y538, E559, N574, D575, S577, Y594, V600, and P605. Sequences were obtained from the respective GenBank accession codes (included parenthetically) as follows: Cel9A_Thefu, from Thermobifida fusca (M73322); CelQ_Cloth, from Clostridium thermocellum (AB047845); Cel9G_Cloce, from Clostridium cellulolyticum (M87018); EngH_Cloce, from Clostridium cellulovorans (U34793); CelZ2_Clost, from Clostridium stercorarium (X55299); Cel9I2_Cloth, from Clostridium thermocellum (L04735); Cel9B_Cloth, from Clostridium thermocellum (X60545); Cel9A_Celfi, from Cellulomonas fimi (M64644); CipC_Cloce, from Clostridium cellulolyticum (U40345); CipA1_Cloth, from Clostridium thermocellum (X67406); CbpA_Cloce, from Clostridium cellulovorans (M73817); CipA_Clojo, from Clostridium josui (AB004845); CelZ1_Clost, from Clostridium stercorarium (X55299); CelA_Calsa, from Caldicellulosiruptor saccharolyticus (M36063); CelA_Paela, from Paenibacillus lautus (M76588); and Cel9I1_Cloth, from Clostridium thermocellum (L04735).
|
|
|
|---|
An interesting finding is that all of the mutations that increase CMC activity are in the catalytic cleft from Glc(2) to Glc(4). The W209S, W256A, R317K, and D261A mutants all have smaller side chains than those in the wild type, which may allow the carboxymethyl-modified groups on CMC to enter and leave the catalytic cleft more easily to achieve more efficient cutting. This hypothesis does not explain why the D261N mutant enzyme has the highest CMC activity of any mutant enzyme studied so far, as the side chain is not changed in size by mutation. There is one hydrogen bond formed between D261 OD2 and Glc(4) O2, and the distance between the two atoms is 2.76 Å. A hydrogen bond is also predicted between N261 ND2 and Glc(4) O2, but the distance between the two atoms is 3.14 Å; therefore, this new hydrogen bond should be weaker than the original one so that the fragment of CMC after hydrolysis might leave the catalytic cleft more easily.
The ability of a cellulase to move along the substrate is critical for crystalline cellulose degradation. However, the mechanism of processivity is not yet completely understood. The R378K mutant is the first mutant with obviously higher processivity. R378 might have the function of moving the cellulose chain toward the catalytic center and into the right subsites. The R378K mutant enzyme loses a hydrogen bond to Glc(+1); thus, it might be easier for the enzyme to slide along the substrate. In contrast, the mutations in residues near subsites Glc(1) to Glc(4) (F205, Y429, W313, W256, R317, and D261) decrease processivity markedly. It is reasonable to propose that these mutations disturb the coordination between the release of the product and the subsequent binding of the cellulose molecule into these subsites as it moves from CBM3c.
Cel9A-51, which retains the linker but none of CBM3c, lost most of the activity on CMC and SWC and had barely measurable crystalline cellulose activity (Table 1). The CD alone has a ratio of soluble to insoluble reducing sugars of only 0.6 compared to 3.1 for Cel9A-68 (4). Clearly the CBM3c domain is essential for cellulose hydrolysis. Cel9A-68 D55C has 32% binding to BC compared to only 10% for Cel9A-51 D55C, indicating that the family 3c domain is also important for cellulose binding. Surprisingly, none of the CBM3c mutations reduced processivity very much; in fact most of them had improved processivity. Even the R557A/E559A double mutant, which has only 15% BC activity, has wild-type processivity. In contrast, all but one of the catalytic domain mutants had reduced processivity; some of them were drastically reduced. This strongly implies that it is the equilibrium binding into the catalytic cleft that controls processivity. It was proposed by Tormo et al. (18) that polar residues can interact with cellulose to contribute to affinity and may replace some of the interchain hydrogen bonds of cellulose so as to disrupt the tight cellulose structure. Four of the CBM3c mutants showed that BC activity was more severely affected than SWC activity, supporting this hypothesis. Based on these studies it appears that the function of CBM3c is to loosely anchor the enzyme on the substrate and to disrupt the hydrogen bond network in crystalline cellulose.
The detailed mechanism of how CBMs disrupt the hydrogen bond network in the crystalline substrate is still unclear. Preliminary data show that single mutations of N470, F476, K480, Y520, Q522, Q553, and R563, which are all suspected binding residues in the CBM3c, do not reduce the enzyme activity very much (the current data and data not shown). One reason could be that the pattern of interaction between the CBM and crystalline cellulose is so complicated, involving many residues simultaneously, that double or triple mutations are required to see changes in activity.
This research was supported by grant DE-FG02-ER15356 from the U.S Department of Energy.
Published ahead of print on 16 March 2007. ![]()
|
|
|---|
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»