Design of a pH-dependent cellulose-binding domain

Protein-carbohydrate interactions typically rely on aromatic stacking interactions of tyrosine, phenylalanine and tryptophan side chains with the sugar rings whereas histidine residues are rarely involved. The small cellulose-binding domain of the Cel7A cellobiohydrolase (formerly CBHI) from Trichoderma reesei binds to crystalline cellulose primarily using a planar strip of three tyrosine side chains. Binding of the wild-type Cel7A CBD is practically insensitive to pH. Here we have investigated how histidine residues mediate the binding interaction and whether the protonation of a histidine side chain makes the binding sensitive to pH. Protein engineering of the Cel7A CBD was thus used to replace the tyrosine residues in two different positions with histidine residues. All of the mutants exhibited a clear pH-dependency of the binding, in clear contrast to the wild-type. Although the binding of the mutants at optimal pH was less than for the wild-type, in one case, Y31H, this binding almost reached the wild-type level.     

Insight into the binding of the wild type and mutated alginate lyase (AlyVI) with its substrate: a computational and experimental study

The homology model of the wild type alginate lyase (AlyVI) marine bacterium Vibrio sp. protein, was built using the crystal structure of the Family 7 alginate lyase from Sphingomonas sp. A1. To rationalize the observed structure-affinity relationships of aliginate lyase alyVI with its (GGG) substrate, molecular docking, MD imulations and binding free energy calculations followed by site-directed mutagenesis and alyVI activity assays were carried out. Per-residue decomposition of the (GGG) binding energy revealed that the most important contributions were from polar and charged residues, such as Asn138, Arg143, Asn217, and Lys308, while van der Waals interactions were responsible for binding with the catalytic His200 and Tyr312 residues. The mutants H200A, K308A, Y312A, Y312F, and W165A were found to be inactive or almost inactive. However, the catalytic efficiency (k(cat)/K(m)) of the double mutant L224V/D226G increased by two-fold compared to the wild type enzyme. This first structural model with its substrate binding mode and the agreement with experimental results provide a suitable base for the future rational design of new mutated alyVI structures with improved catalytic activity.     

Structural Insights into the Affinity of Cel7A Carbohydrate-binding Module for Lignin

The high cost of hydrolytic enzymes impedes the commercial production of lignocellulosic biofuels. High enzyme loadings are required in part due to their non-productive adsorption to lignin, a major component of biomass. Despite numerous studies documenting cellulase adsorption to lignin, few attempts have been made to engineer enzymes to reduce lignin binding. In this work, we used alanine-scanning mutagenesis to elucidate the structural basis for the lignin affinity of Trichoderma reesei Cel7A carbohydrate binding module (CBM). T. reesei Cel7A CBM mutants were produced with a Talaromyces emersonii Cel7A catalytic domain and screened for their binding to cellulose and lignin. Mutation of aromatic and polar residues on the planar face of the CBM greatly decreased binding to both cellulose and lignin, supporting the hypothesis that the cellulose-binding face is also responsible for lignin affinity. Cellulose and lignin affinity of the 31 mutants were highly correlated, although several mutants displayed selective reductions in lignin or cellulose affinity. Four mutants with increased cellulose selectivity (Q2A, H4A, V18A, and P30A) did not exhibit improved hydrolysis of cellulose in the presence of lignin. Further reduction in lignin affinity while maintaining a high level of cellulose affinity is thus necessary to generate an enzyme with improved hydrolysis capability. This work provides insights into the structural underpinnings of lignin affinity, identifies residues amenable to mutation without compromising cellulose affinity, and informs engineering strategies for family one CBMs.     

Energetic analysis of binding of progesterone and 5 beta-androstane-3,17-dione to anti-progesterone antibody DB3 using molecular dynamics and free energy calculations.

Molecular dynamics simulations and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) free energy calculations were used to study the energetics of the binding of progesterone (PRG) and 5 beta-androstane-3,17-dione (5AD) to anti-PRG antibody DB3. Although the two steroids bind to DB3 in different orientations, their binding affinities are of the same magnitude, 1 nM for PRG and 8 nM for 5AD. The calculated relative binding free energy of the steroids, 8.8 kJ/mol, is in fair agreement with the experimental energy, 5.4 kJ/mol. In addition, computational alanine scanning was applied to study the role of selected amino acid residues of the ligand-binding site on the steroid cross-reactivity. The electrostatic and van der Waals components of the total binding free energies were found to favour more the binding of PRG, whereas solvation energies were more favourable for the binding of 5AD. The differences in the free energy components are due to the binding of the A rings of the steroids to different binding pockets: PRG is bound to a pocket in which electrostatic antibody-steroid interactions are dominating, whereas 5AD is bound to a pocket in which van der Waals and hydrophobic interactions dominate.     

Binding site analysis of cellulose binding domain CBD(N1) from endoglucanse C of Cellulomonas fimi by site-directed mutagenesis

Endoglucanase C (CenC), a beta1,4 glucanase from the soil bacterium Cellulomonas fimi, binds to amorphous cellulose via two homologous cellulose binding domains, termed CBD(N1) and CBD(N2). In this work, the contributions of 10 amino acids within the binding cleft of CBD(N1) were evaluated by single site-directed mutations to alanine residues. Each isolated domain containing a single mutation was analyzed for binding to an insoluble amorphous preparation of cellulose, phosphoric acid swollen Avicel (PASA), and to a soluble glucopyranoside polymer, barley beta-glucan. The effect of any given mutation on CBD binding was similar for both substrates, suggesting that the mechanism of binding to soluble and insoluble substrates is the same. Tyrosines 19 and 85 were essential for tight binding by CBD(N1) as their replacement by alanine results in affinity decrements of approximately 100-fold on PASA, barley beta-glucan, and soluble cellooligosaccharides. The tertiary structures of unbound Y19A and Y85A were assessed by heteronuclear single quantum coherence (HSQC) spectroscopy. These studies indicated that the structures of both mutants were perturbed but that all perturbations are very near to the site of mutation.     

Probing carbohydrate product expulsion from a processive cellulase with multiple absolute binding free energy methods

Understanding the enzymatic mechanism that cellulases employ to degrade cellulose is critical to efforts to efficiently utilize plant biomass as a sustainable energy resource. A key component of cellulase action on cellulose is product inhibition from monosaccharide and disaccharides in the product site of cellulase tunnel. The absolute binding free energy of cellobiose and glucose to the product site of the catalytic tunnel of the Family 7 cellobiohydrolase (Cel7A) of Trichoderma reesei (Hypocrea jecorina) was calculated using two different approaches: steered molecular dynamics (SMD) simulations and alchemical free energy perturbation molecular dynamics (FEP/MD) simulations. For the SMD approach, three methods based on Jarzynski's equality were used to construct the potential of mean force from multiple pulling trajectories. The calculated binding free energies, -14.4 kcal/mol using SMD and -11.2 kcal/mol using FEP/MD, are in good qualitative agreement. Analysis of the SMD pulling trajectories suggests that several protein residues (Arg-251, Asp-259, Asp-262, Trp-376, and Tyr-381) play key roles in cellobiose and glucose binding to the catalytic tunnel. Five mutations (R251A, D259A, D262A, W376A, and Y381A) were made computationally to measure the changes in free energy during the product expulsion process. The absolute binding free energies of cellobiose to the catalytic tunnel of these five mutants are -13.1, -6.0, -11.5, -7.5, and -8.8 kcal/mol, respectively. The results demonstrated that all of the mutants tested can lower the binding free energy of cellobiose, which provides potential applications in engineering the enzyme to accelerate the product expulsion process and improve the efficiency of biomass conversion.     

Probing the structural and energetic basis of kinesin-microtubule binding using computational alanine-scanning mutagenesis

Kinesin-microtubule (MT) binding plays a critical role in facilitating and regulating the motor function of kinesins. To obtain a detailed structural and energetic picture of kinesin-MT binding, we performed large-scale computational alanine-scanning mutagenesis based on long-time molecular dynamics (MD) simulations of the kinesin-MT complex in both ADP and ATP states. First, we built three all-atom kinesin-MT models: human conventional kinesin bound to ADP and mouse KIF1A bound to ADP and ATP. Then, we performed 30 ns MD simulations followed by kinesin-MT binding free energy calculations for both the wild type and mutants obtained after substitution of each charged residue of kinesin with alanine. We found that the kinesin-MT binding free energy is dominated by van der Waals interactions and further enhanced by electrostatic interactions. The calculated mutational changes in kinesin-MT binding free energy are in excellent agreement with results of an experimental alanine-scanning study with a root-mean-square error of ~0.32 kcal/mol [Woehlke, G., et al. (1997) Cell 90, 207-216]. We identified a set of important charged residues involved in the tuning of kinesin-MT binding, which are clustered on several secondary structural elements of kinesin (including well-studied loops L7, L8, L11, and L12, helices α4, α5, and α6, and less-explored loop L2). In particular, we found several key residues that make different contributions to kinesin-MT binding in ADP and ATP states. The mutations of these residues are predicted to fine-tune the motility of kinesin by modulating the conformational transition between the ADP state and the ATP state of kinesin.     

Alpha-amylase inhibitors selected from a combinatorial library of a cellulose binding domain scaffold

A disulfide bridge-constrained cellulose binding domain (CBD(WT)) derived from the cellobiohydrolase Cel7A from Trichoderma reesei has been investigated for use in scaffold engineering to obtain novel binding proteins. The gene encoding the wild-type 36 aa CBD(WT) domain was first inserted into a phagemid vector and shown to be functionally displayed on M13 filamentous phage as a protein III fusion protein with retained cellulose binding activity. A combinatorial library comprising 46 million variants of the CBD domain was constructed through randomization of 11 positions located at the domain surface and distributed over three separate beta-sheets of the domain. Using the enzyme porcine alpha-amylase (PPA) as target in biopannings, two CBD variants showing selective binding to the enzyme were characterized. Reduction and iodoacetamide blocking of cysteine residues in selected CBD variants resulted in a loss of binding activity, indicating a conformation dependent binding. Interestingly, further studies showed that the selected CBD variants were capable of competing with the binding of the amylase inhibitor acarbose to the enzyme. In addition, the enzyme activity could be partially inhibited by addition of soluble protein, suggesting that the selected CBD variants bind to the active site of the enzyme.     

Chitin binding by Thermobifida fusca cellulase catalytic domains

Cellulose is a linear homopolymer of beta 1-4 linked glucose residues. Chitin is similar to cellulose in structure, and can be described as cellulose with the hydroxyl group on the C2 carbon replaced by an acetylamine group. Both cellulose and chitin form tightly packed, extensively hydrogen-bonded micro-fibrils. Up to now, binding of cellulase catalytic domains (CDs) to chitin has not been reported. In this article, binding of the CDs of Thermobifida fusca Cel6A, Cel6B, Cel48A, Cel5A, and Cel9A to alpha-chitin was investigated. The CDs of endocellulases, Cel6A and Cel5A did not bind to alpha-chitin; one exocellulase, Cel48A CD bound alpha-chitin moderately well; and the exocellulase Cel6B CD and the processive endocellulase Cel9A CD bound extremely tightly to alpha-chitin. Only mutations of Cel6B W329C, W332A and G234S and Cel9A Y206F, Y206S and D261A/R378K caused weaker binding to alpha-chitin than wild-type, and all these mutations were of residues near the catalytic center. One mutant enzyme, Cel9A D261A/R378K had weak chitinase activity, but no soluble products were detected. Chitotriose and chitotetraose were docked successfully to the catalytic cleft of Cel9A. In general, the positioning of the sugar residues in the model structures matched the cellooligosaccharides in the X-ray structure. Our results show that the binding of chitin by a cellulase can provide additional information about its binding to cellulose.     

Three-dimensional structures of three engineered cellulose-binding domains of cellobiohydrolase I from Trichoderma reesei.

Three-dimensional solution structures for three engineered, synthetic CBDs (Y5A, Y31A, and Y32A) of cellobiohydrolase I (CBHI) from Trichoderma reesei were studied with nuclear magnetic resonance (NMR) and circular dichroism (CD) spectroscopy. According to CD measurements the antiparallel beta-sheet structure of the CBD fold was preserved in all engineered peptides. The three-dimensional NMR-based structures of Y31A and Y32A revealed only small local changes due to mutations in the flat face of CBD, which is expected to bind to crystalline cellulose. Therefore, the structural roles of Y31 and Y32 are minor, but their functional importance is obvious because these mutants do not bind strongly to cellulose. In the case of Y5A, the disruption of the structural framework at the N-terminus and the complete loss of binding affinity implies that Y5 has both structural and functional significance. The number of aromatic residues and their precise spatial arrangement in the flat face of the type I CBD fold appears to be critical for specific binding. A model for the CBD binding in which the three aligned aromatic rings stack onto every other glucose ring of the cellulose polymer is discussed.     

A model for glutathione binding and activation in the fosfomycin resistance protein, FosA

The genomically encoded fosfomycin resistance protein from Pseudomonas aeruginosa (FosA(PA)) utilizes Mn(II) and K(+) to catalyze the addition of glutathione (GSH) to C1 of the antibiotic rendering it inactive. Although this protein has been structurally and kinetically characterized with respect to the substrate, fosfomycin, questions remain regarding how the enzyme binds the thiol substrate, GSH. Computational studies have revealed a potential GSH binding site in FosA(PA) that involves six electrostatic or hydrogen-bonding interactions with protein side-chains as well as six additional residues that contribute van der Waals interactions. A strategically placed tyrosine residue, Y39, appears to be involved in the ionization of GSH during catalysis. The Y39F mutant exhibits a 13-fold reduction of catalytic activity (k(cat)=14+/-2s(-1)), suggesting a role in the ionization of GSH. Mutation of five other residues (W34, Q36, S50, K90, and R93) implicated in ionic of hydrogen-bonding interactions resulted in enzymes with reduced catalytic efficiency, affinity for GSH, or both. The mutant enzymes were also found to be less effective resistant proteins in the biological context of Escherichia coli. The more conservative W34H mutant has native-like catalytic efficiency suggesting that the imidazole NH group can replace the indole group of W34 that is important for GSH binding. In the absence of co-crystal structural data with the thiol substrate, these results provide important insights into the role of GSH in catalysis.     

Hypocrea jecorina (Trichoderma reesei) Cel7A as a molecular machine: A docking study

Hypocrea jecorina (formerly Trichoderma reesei) Cel7A has a catalytic domain (CD) and a cellulose-binding domain (CBD) separated by a highly glycosylated linker. Very little is known of how the 2 domains interact to degrade crystalline cellulose. Based on the interaction energies and forces on cello-oligosaccharides computationally docked to the CD and CBD, we propose a molecular machine model, where the CBD wedges itself under a free chain end on the crystalline cellulose surface and feeds it to the CD active site tunnel. Enzyme-substrate interactions produce the forces required to pull cellulose chains from the surface and also to help the enzyme move on the cellulose chain for processive hydrolysis. The energy to generate these forces is ultimately derived from the chemical energy of glycosidic bond breakage.     

Structural basis for enantiomer binding and separation of a common beta-blocker: crystal structure of cellobiohydrolase Cel7A with bound (S)-propranolol at 1.9 A resolution

Cellobiohydrolase Cel7A (previously called CBH 1), the major cellulase produced by the mould fungus Trichoderma reesei, has been successfully exploited as a chiral selector for separation of stereo-isomers of some important pharmaceutical compounds, e.g. adrenergic beta-blockers. Previous investigations, including experiments with catalytically deficient mutants of Cel7A, point unanimously to the active site as being responsible for discrimination of enantiomers.In this work the structural basis for enantioselectivity of basic drugs by Cel7A has been studied by X-ray crystallography. The catalytic domain of Cel7A was co-crystallised with the (S)-enantiomer of a common beta-blocker, propranolol, at pH 7, and the structure of the complex was determined and refined at 1. 9 A resolution. Indeed, (S)-propranolol binds at the active site, in glucosyl-binding subsites -1/+1. The catalytic residues Glu212 and Glu217 make tight salt links with the secondary amino group of (S)-propranolol. The oxygen atom attached to the chiral centre of (S)-propranolol forms hydrogen bonds to the nucleophile Glu212 O(epsilon1) and to Gln175 N(epsilon2), whereas the aromatic naphthyl moiety stacks with the indole ring of Trp376 in site +1. The bidentate charge interaction with the catalytic glutamate residues is apparently crucial, since no enantioselectivity has been obtained with the catalytically deficient mutants E212Q and E217Q.Activity inhibition experiments with wild-type Cel7A were performed in conditions close to those used for crystallisation. Competitive inhibition constants for (R)- and (S)-propranolol were determined at 220 microM and 44 microM, respectively, corresponding to binding free energies of 20 kJ/mol and 24 kJ/mol, respectively. The K(i) value for (R)-propranolol was 57-fold lower than the highest concentration, 12.5 mM, used in co-crystallisation experiments. Still several attempts to obtain a complex with the (R)-enantiomer have failed.By using cellobiose as a selective competing ligand, the retention of the enantiomers of propranolol on the chiral stationary phase (CSP) based on Cel7A mutant D214N were resolved into enantioselective and non- selective binding. The enantioselective binding was weaker for both enantiomers on D214N-CSP than on wild-type-CSP.     

Functional Analysis of the Carbohydrate-Binding Domains of Erwinia chrysanthemi Cel5 (Endoglucanase Z) and an Escherichia coli Putative Chitinase

The Cel5 cellulase (formerly known as endoglucanase Z) from Erwinia chrysanthemi is a multidomain enzyme consisting of a catalytic domain, a linker region, and a cellulose binding domain (CBD). A three-dimensional structure of the CBD(Cel5) has previously been obtained by nuclear magnetic resonance. In order to define the role of individual residues in cellulose binding, site-directed mutagenesis was performed. The role of three aromatic residues (Trp18, Trp43, and Tyr44) in cellulose binding was demonstrated. The exposed potential hydrogen bond donors, residues Gln22 and Glu27, appeared not to play a role in cellulose binding, whereas residue Asp17 was found to be important for the stability of Cel5. A deletion mutant lacking the residues Asp17 to Pro23 bound only weakly to cellulose. The sequence of CBD(Cel5) exhibits homology to a series of five repeating domains of a putative large protein, referred to as Yheb, from Escherichia coli. One of the repeating domains (Yheb1), consisting of 67 amino acids, was cloned from the E. coli chromosome and purified by metal chelating chromatography. While CBD(Cel5) bound to both cellulose and chitin, Yheb1 bound well to chitin, but only very poorly to cellulose. The Yheb protein contains a region that exhibits sequence homology with the catalytic domain of a chitinase, which is consistent with the hypothesis that the Yheb protein is a chitinase.     

Quantification of helix-helix binding affinities in micelles and lipid bilayers

A theoretical approach for estimating association free energies of alpha-helices in nonpolar media has been developed. The parameters of energy functions have been derived from DeltaDeltaG values of mutants in water-soluble proteins and partitioning of organic solutes between water and nonpolar solvents. The proposed approach was verified successfully against three sets of published data: (1) dissociation constants of alpha-helical oligomers formed by 27 hydrophobic peptides; (2) stabilities of 22 bacteriorhodopsin mutants, and (3) protein-ligand binding affinities in aqueous solution. It has been found that coalescence of helices is driven exclusively by van der Waals interactions and H-bonds, whereas the principal destabilizing contributions are represented by side-chain conformational entropy and transfer energy of atoms from a detergent or lipid to the protein interior. Electrostatic interactions of alpha-helices were relatively weak but important for reproducing the experimental data. Immobilization free energy, which originates from restricting rotational and translational rigid-body movements of molecules during their association, was found to be less than 1 kcal/mole. The energetics of amino acid substitutions in bacteriorhodopsin was complicated by specific binding of lipid and water molecules to cavities created in certain mutants.     

Role of hydrogen bonding in the interaction between a xylan binding module and xylan

NMR studies of the internal family 2b carbohydrate binding module (CBM2b-1) of Cellulomonas fimi xylanase 11A have identified six polar residues and two aromatic residues that interact with its target ligand, xylan. To investigate the importance of the various interactions, free energy and enthalpy changes have been measured for the binding of xylan to native and mutant forms of CBM2b-1. The data show that the two aromatic residues, Trp 259 and Trp 291, play a critical role in the binding, and similarly that mutants N264A and T316A have no affinity for the xylose polymer. Interestingly, mutations E257A, Q288A, N292A, E257A/Q288A, E257A/N292A, and E257A/N292A/Q288A do not significantly diminish the affinity of CBM2b-1 for the xylose polymers, but do influence the thermodynamics driving the protein-carbohydrate interactions. These thermodynamic parameters have been interpreted in light of a fresh understanding of enthalpy-entropy compensation and show the following. (1) For proteins whose ligands are bound on an exposed surface, hydrogen bonding confers little specificity or affinity. It also displays little cooperativity. Most specificity and affinity derive from binding between the face of sugar rings and aromatic rings. (2) Loss of hydrogen bonding interactions leads to a redistribution of the remaining bonding interactions such that the entropic mobility of the ligand is maximized, at the expense (if necessary) of enthalpically favorable bonds. (3) Changes in entropy and enthalpy in the binding between polysaccharide and a range of mutants can be interpreted by considering changes in binding and flexibility, without any need to consider solvent reorganization.     

Free energy calculations for the relative binding affinity between DNA and lambda-repressor

The change in the binding free energy between DNA and lambda-repressor resulting from a base substitution, thymine (T)-->deoxyuracil (abbreviated as U), was evaluated by the free energy perturbation method on the basis of molecular dynamics simulations for the DNA-lambda-repressor complex in water with all degrees of freedom and including long-range Coulomb interactions. The binding free energy change that we calculated (1.47 +/- 0.40 kcal/mol) was in good agreement with an experimental value (1.8 kcal/mol). We clarified why the small difference between T and U (CH(3) in T is replaced with H in U) caused such a significant change in the binding free energy: The substitution of CH(3) in T with H in U lowered the dissociated-state free energy level due to the gain of the hydration free energy. Furthermore, the T-->U substitution raised the free energy level in the associated state due to the loss of the favored van der Waals (vdW) interactions with the lambda-repressor amino acid residues. In other words, the amino acid residues of lambda-repressor can recognize the CH(3) in T through the vdW interactions with the CH(3). This recognition is enhanced in a water environment, because the hydrophobic CH(3) prefers the amino acid residues of lambda-repressor to water molecules.     

A kinetic study of Trichoderma reesei Cel7B catalyzed cellulose hydrolysis

One prominent feature of Trichoderma reesei (Tr) endoglucanases catalyzed cellulose hydrolysis is that the reaction slows down quickly after it starts (within minutes). But the mechanism of the slowdown is not well understood. A structural model of Tr- Cel7B catalytic domain bound to cellulose was built computationally and the potentially important binding residues were identified and tested experimentally. The 13 tested mutants show different binding properties in the adsorption to phosphoric acid swollen cellulose and filter paper. Though the partitioning parameter to filter paper is about 10 times smaller than that to phosphoric acid swollen cellulose, a positive correlation is shown for two substrates. The kinetic studies show that the reactions slow down quickly for both substrates. This slowdown is not correlated to the binding constant but anticorrelated to the enzyme initial activity. The amount of reducing sugars released after 24 h by Cel7B in phosphoric acid swollen cellulose, Avicel and filter paper cellulose hydrolysis is correlated with the enzyme activity against a soluble substrate p-nitrophenyl lactoside. Six of the 13 tested mutants, including N47A, N52D, S99A, N323D, S324A, and S346A, yield ∼15–35% more reducing sugars than the wild type (WT) Cel7B in phosphoric acid swollen cellulose and filter paper hydrolysis. This study reveals that the slowdown of the reaction is not due to the binding of the enzyme to cellulose. The activity of Tr- Cel7B against the insoluble substrate cellulose is determined by the enzyme’s capability in hydrolyzing the soluble substrate.     

Properties and mutation analysis of the CelK cellulose-binding domain from the Clostridium thermocellum cellulosome

The family IV cellulose-binding domain of Clostridium thermocellum CelK (CBD(CelK)) was expressed in Escherichia coli and purified. It binds to acid-swollen cellulose (ASC) and bacterial microcrystalline cellulose (BMCC) with capacities of 16.03 and 3.95 micromol/g of cellulose and relative affinities (K(r)) of 2.33 and 9.87 liters/g, respectively. The CBD(CelK) is the first representative of family IV CBDs to exhibit an affinity for BMCC. The CBD(CelK) also binds to the soluble polysaccharides lichenin, glucomannan, and barley beta-glucan, which are substrates for CelK. It does not bind to xylan, galactomannan, and carboxymethyl cellulose. The CBD(CelK) contains 1 mol of calcium per mol. The CBD(CelK) has three thiol groups and one disulfide, reduction of which results in total loss of cellulose-binding ability. To reveal amino acid residues important for biological function of the domain and to investigate the role of calcium in the CBD(CelK) four highly conserved aromatic residues (Trp(56), Trp(94), Tyr(111), and Tyr(136)) and Asp(192) were mutated into alanines, giving the mutants W56A, W94A, Y111A, Y136A, and D192A. In addition 14 N-terminal amino acids were deleted, giving the CBD-N(CelK). The CBD-N(CelK) and D192A retained binding parameters close to that of the intact CBD(CelK), W56A and W94A totally lost the ability to bind to cellulose, Y136A bound to both ASC and BMCC but with significantly reduced binding capacity and K(r) and Y111A bound weakly to ASC and did not bind to BMCC. Mutations of the aromatic residues in the CBD(CelK) led to structural changes revealed by studying solubility, circular-dichroism spectra, dimer formation, and aggregation. Calcium content was drastically decreased in D192A. The results suggest that Asp192 is in the calcium-binding site of the CBD(CelK) and that calcium does not affect binding to cellulose. The 14 amino acids from the N terminus of the CBD(CelK) are not important for binding. Tyr136, corresponding to Cellulomonas fimi CenC CBD(N1) Y85, located near the binding cleft, might be involved in the formation of the binding surface, while Y111, W56A, and W94A are essential for the binding process by keeping the CBD(CelK) correctly folded.