Mutation of non-conserved amino acids surrounding catalytic site to shift pH optimum of Bacillus circulans xylanase

Improvement of enzyme function by engineering pH dependence of enzymatic activity is of importance for industrial application of Bacillus circulans xylanases. Target mutation sites were selected by structural alignment between B. circulans xylanase and other xylanases having different pH optima. We selected non-conserved mutant sites within 8 Å from the catalytic residues, to see whether these residues have some role in modulating pK_as of the catalytic residues. We hypothesized that the non-conserved residues which may not have any role in enzyme catalysis might perturb pK_as of the catalytic residues. Change in pK_a of a titratable group due to change in electrostatic potential of a mutation was calculated and the change in pH optimum was predicted from the change in pK_a of the catalytic residues. Our strategy is proved to be useful in selection of promising mutants to shift the pH optimum of the xylanases towards desired side.     

Shifting pH optimum of Bacillus circulans xylanase based on molecular modeling

Although hydrolases are used in several industrial processes, its industrial applications have some limitations in specific cases since some industrial processes are carried out at pH value which is different from optimum pH of the enzyme. Alkaline side optimum pH of hydrolases is always desirable, and it is proved difficult to achieve that by mutation only. Hence, molecular modeling was applied to select the promising mutants. The changes in electrostatic potential, which was calculated using Delphi, were compared to the changes in pH optimum of four hydolases and their mutants. The results showed that the change in electrostatic potential can be used as an indicator for selecting relevant candidates of mutation. Bacillus circulans xylanase was selected as a model system, and the promising mutants were picked up by the molecular modeling. Q167M and R73V, had a higher pH optimum than the wild type, while K175Q had a similar pH-activity profile of the wild type.     

Shifting the optimum pH of Bacillus circulans xylanase towards acidic side by introducing arginine

Electrostatic interactions are important in protein folding, binding, flexibility, stability and function. The pH at which the enzyme is maximally active is determined by the pK_as of the active site residues, which are modulated by several factors including the change in electrostatics in its vicinity. As the acidic xylanases are important in food and animal feed industries, electrostatic interactions are introduced in Bacillus circulans xylanase to shift their pH optima towards the acidic side. Arg substitutions are made to modulate the pKas of the active site residues. Neutral residues are substituted by Arg in such a way that the substituted residue can make direct interaction with the catalytic residues. However, the mutations with other titratable residues (Asp, Arg, Lys, His, Tyr, and Ser) present in between the catalytic sites and the substituted sites are avoided. Site directed mutagenesis was conducted to confirm the strategy. The results show the shift in pH optima of the mutants towards the acidic side by 0.5–1.5 unit. Molecular dynamics simulation of the mutant V37R reveals that the decrease in activity is due to the increase in distance between the substrate oxygen atoms and catalytic glutamates.     

Identification of catalytic amino acids of cyclodextran glucanotransferase from Bacillus circulans T-3040

In glycoside hydrolase family 66 (see http://afmb.cnrs-mrs.fr/CAZY/), cyclodextran glucanotransferase (CITase) is the only transglycosylation enzyme, all the other family 66 enzymes being dextranases. To analyze the catalytic amino acids of CITase, we modified CITase chemically from the T-3040 strain of Bacillus circulans with 1-ethyl-3-(3-dimethylaminopropyl)-carbodiimide (EDC). EDC inactivated the enzyme by following pseudo-first order kinetics. In addition, the substrates of an isomaltooligosaccharide and a cyclodextran inhibited EDC-induced enzyme inactivation, implicating the carboxyl groups of CITase as the catalytic amino acids of the enzyme. When two conserved aspartic acid residues, Asp145 and Asp270, were replaced with Asn in T-3040 mature CITase, CIT-D270N was completely inactive, and CIT-D145N had reduced activity. The V(max) of CIT-D145N was 1% of that of wild-type CITase, whereas the K(m) of CIT-D145N was about the same as that of the wild-type enzyme. These findings indicate that Asp145 and Asp270 play an important role in the enzymatic reaction of T-3040 CITase.     

Dissecting electrostatic interactions in Bacillus circulans xylanase through NMR-monitored pH titrations

NMR-monitored pH titration curves of proteins provide a rich source of structural and electrostatic information. Although relatively straightforward to measure, interpreting pH-dependent chemical shift changes to obtain site-specific acid dissociation constants (pK (A) values) is challenging. In order to analyze the biphasic titrations exhibited by the side chain (13)C(γ) nuclei of the nucleophilic Glu78 and general acid/base Glu172 in Bacillus circulans xylanase, we have revisited the formalism for the ionization equilibria of two coupled acidic residues. In general, fitting NMR-monitored pH titration curves for such a system will only yield the two macroscopic pK (A) values that reflect the combined effects of both deprotonation reactions. However, through the use of mutations complemented with ionic strength-dependent measurements, we are able to extract the four microscopic pK (Ai) values governing the branched acid/base equilibria of Glu78 and Glu172 in BcX. These data, confirmed through theoretical calculations, help explain the pH-dependent mechanism of this model GH11 xylanase by demonstrating that the kinetically determined pK (A) values and hence catalytic roles of these two residues result from their electrostatic coupling.     

Expression of Bacillus circulans Teri-42 xylanase gene in Bacillus subtilis

The xylanase gene of Bacillus circulans Teri-42 was cloned in both B. subtilis and Escherichia coli. The enzyme activity was almost 87% higher in B. subtilis (pBA7) than in E. coli (pAQ4). No cellulase activity was detected in the clones, B. subtilis (pBA7) and E. coli (pAQ4). Approximately 1120 U (80%) of the xylanase was secreted extracellularly by the clone B. subtilis (pBA7) as compared to 79 U (88%) excreted in E. coli (pAQ4). In B. subtilis (pBA7) the optimal xylanase activity was at pH 7.0 and 50 degrees C, which was the same as that of the parent B. circulans Teri-42. The recombinant xylanase in B. subtilis was more stable at higher temperatures than the parent B. circulans Teri-42. Purification of xylanase from the clone B. subtilis (pBA7) showed a 71 kDa polypeptide similar to that observed in B. circulans Teri-42.     

Mutational and crystallographic analyses of the active site residues of the Bacillus circulans xylanase.

Using site-directed mutagenesis we have investigated the catalytic residues in a xylanase from Bacillus circulans. Analysis of the mutants E78D and E172D indicated that mutations in these conserved residues do not grossly alter the structure of the enzyme and that these residues participate in the catalytic mechanism. We have now determined the crystal structure of an enzyme-substrate complex to 108 A resolution using a catalytically incompetent mutant (E172C). In addition to the catalytic residues, Glu 78 and Glu 172, we have identified 2 tyrosine residues, Tyr 69 and Tyr 80, which likely function in substrate binding, and an arginine residue, Arg 112, which plays an important role in the active site of this enzyme. On the basis of our work we would propose that Glu 78 is the nucleophile and that Glu 172 is the acid-base catalyst in the reaction.     

Secondary structure and NMR assignments of Bacillus circulans xylanase.

Bacillus circulans xylanase (BCX) is a member of the family of low molecular weight endo-beta-(1,4)-xylanases. The main-chain 1H, 13C, and 15N resonances of this 20.4-kDa enzyme were assigned using heteronuclear NMR experiments recorded on a combination of selectively and uniformly labeled protein samples. Using chemical shift, NOE, J coupling, and amide hydrogen exchange information, 14 beta-strands, arranged in a network of three beta-sheets, and a single alpha-helix were identified in BCX. The NMR-derived secondary structure and beta-sheet topology agree closely with that observed in the crystal structure of this protein. The HN of Ile 118 has a strongly upfield-shifted resonance at 4.03 ppm, indicative of a potential amide-aromatic hydrogen bond to the indole ring of Trp 71. This interaction, which is conserved in all low molecular weight xylanases of known structure, may play an important role in establishing the active site conformation of these enzymes. Following hen egg white and bacteriophage T4 lysozymes, B. circulans xylanase represents the third family of beta-glycanases for which extensive NMR assignments have been reported. These assignments provide the background for detailed studies of the mechanism of carbohydrate recognition and hydrolysis by this bacterial xylanase.     

Mechanism of stabilization of Bacillus circulans xylanase upon the introduction of disulfide bonds

The introduction of disulfide bonds has been used as a strategy to enhance the stability of Bacillus circulans xylanase. The transition temperature of the S100C/N148C (DS1), V98C/A152C (DS2), and A1GC/G187,C188 (cXl) in comparison to the wild type was increased by 5.0, 4.1 and 3.8 degrees C, respectively. Interestingly, a combination of two disulfide bonds of DS1 and cXl (cDS1, circular disulfide 1) led to a 12 degrees C increase in the transition temperature. Importantly, an increase in the melting point and DeltaDeltaG values of the cDS1 mutant was cooperative. These results suggest that the mechanism of stabilization by disulfide bonds under irreversible denaturation condition is achieved through: (1) a change in the rate-limiting step on the denaturation pathway; (2) destabilizing the unfolded state without affecting the relative rate constants on the denaturation pathway (like cXl mutant); and (3) or combination of the two (cDS1 mutant).     

Refolding the unfoldable: A systematic approach for renaturation of Bacillus circulans xylanase

Xylanases are important polysaccharide‐cleaving catalysts for the pulp and paper, animal feeds and biofuels industries. They have also proved to be valuable model systems for understanding enzymatic catalysis, with one of the best studied being the GH11 xylanase from Bacillus circulans (Bcx). However, proteins from this class are very recalcitrant to refolding in vitro. This both limits their high level expression in heterologous hosts, and prevents experimental approaches, such as peptide ligation or chemical modifications, to probe and engineer their stability and function. To solve this problem, a systematic screening approach was employed to identify suitable buffer conditions for renaturing Bcx in vitro. The fractional factorial screen employed identified starting conditions for refolding, which were then refined and developed into a generic protocol for renaturing preparative amounts of active Bcx in a 50–60% yield from inclusion bodies. The method is robust and proved equally proficient at refolding circularly permuted versions that carry cysteine mutations. This general approach should be applicable to related GH11 xylanases, as well as proteins adopting a similar β‐jellyroll fold, that are otherwise recalcitrant to refolding in vitro.     

Statistical optimization of thermo-tolerant xylanase activity from Amazon isolated Bacillus circulans on solid-state cultivation

A 2(2) factorial design was performed to find the best conditions of pH and temperature for xylanolytic activity of Bacillus circulans BL53 isolated from the Amazon environment. Solid-state cultivation was carried out on an inexpensive, abundant agro-industrial soybean residue. The central composite design (CCD) used for the analysis of treatment combinations showed that a second-order polynomial regression model was in good agreement with experimental results, with R(2) = 0.9369 (P < 0.05). The maximum activity was obtained at a high temperature (80 degrees C) and over a large pH range (4.0-7.0). Enzymatic activity was maintained in heated extracts up to 50 degrees C, suggesting that the xylanases of B. circulans BL53 are thermo-tolerant biocatalysts, being of interest for industrial processes. The crude enzyme extract hydrolyzed rice straw, sugar cane bagasse and soybean fiber and its activity was stimulated by Co(2+), Fe(3+), and beta-mercaptoethanol but inhibited by Mn(2+), Cu(2+), Ca(2+), Zn(2+), Ba(2+), Mg(2+) and by EDTA.     

A single amino acid substitution enhances the catalytic activity of family 11 xylanase at alkaline pH

Random mutagenesis of the gene encoding family 11 xylanase was used to obtain alkalophilic mutants. The catalytic domain of the chimeric enzyme Stx15, which was constructed from Streptomyces lividans xylanase B and Thermobifida fusca xylanase A, was mutated using error-prone PCR and screened for halo formation on dye-linked xylan plates and activity toward soluble xylan. A positive mutant, M1011, was isolated, and it was found that mutation A49V was responsible for the alkalophilicity of the mutant. Mutation A49V increased the specific activity at pH 9.1 and the stability of mutant A49V was not significantly different from that of Stx15 at 60 degrees C. Both enzymes retained more than 90% of their relative activity from pH 4.7 to 9.1 after 1 h of incubation at 60 degrees C. Analysis of the kinetic parameters at various pH values showed that the A49V mutation reduced the Km in the alkaline pH range, resulting in the higher specific activity of the A49V mutant enzyme.     

Enhancing the activity of Bacillus circulans xylanase by modulating the flexibility of the hinge region

Enzymes undergo multiple conformational changes in solution, and these dynamics are considered to play a critical role in enzyme activity. Hinge-bending motions, resulting from reciprocal movements of dynamical quasi-rigid bodies, are thought to be related to turnover rate and are affected by the physical properties of the hinge regions. In this study, hinge identification and flexibility modification of the regions by mutagenesis were conducted to explore the relationship between hinge flexibility and catalytic activity. Bacillus circulans xylanase was selected for the identification and mutation of the hinge regions. As a result, turnover rate (V _max) was improved approximately twofold in mutants that have more rigid hinge structure, despite the decrease in K _m and V _max/K _m. This result indicates that the rigidly mutated hinge has positive effects on B. circulans xylanase activity.     

Characterization of a buried neutral histidine residue in Bacillus circulans xylanase: NMR assignments, pH titration, and hydrogen exchange.

Bacillus circulans xylanase contains two histidines, one of which (His 156) is solvent exposed, whereas the other (His 149) is buried within its hydrophobic core. His 149 is involved in a network of hydrogen bonds with an internal water and Ser 130, as well as a potential weak aromatic-aromatic interaction with Tyr 105. These three residues, and their network of interactions with the bound water, are conserved in four homologous xylanases. To probe the structural role played by His 149, NMR spectroscopy was used to characterize the histidines in BCX. Complete assignments of the 1H, 13C, and 15N resonances and tautomeric forms of the imidazole rings were obtained from two-dimensional heteronuclear correlation experiments. An unusual spectroscopic feature of BCX is a peak near 12 ppm arising from the nitrogen bonded 1H epsilon 2 of His 149. Due to its solvent inaccessibility and hydrogen bonding to an internal water molecule, the exchange rate of this proton (4.0 x 10(-5) s-1 at pH*7.04 and 30 degrees C) is retarded by > 10(6)-fold relative to an exposed histidine. The pKa of His 156 is unperturbed at approximately 6.5, as measured from the pH dependence of the 15N- and 1H-NMR spectra of BCX. In contrast, His 149 has a pKa < 2.3, existing in the neutral N epsilon 2H tautomeric state under all conditions examined. BCX unfolds at low pH and 30 degrees C, and thus His 149 is never protonated significantly in the context of the native enzyme. The structural importance of this buried histidine is confirmed by the destablizing effect of substituting a phenylalanine or glutamine at position 149 in BCX.     

Production of a thermostable alkali-tolerant xylanase from Bacillus circulans AB 16 grown on wheat straw

Bacillus circulans AB 16 was able to produce 50 IU/ml of xylanase, with negligible cellulase activity when grown on untreated wheat straw. The pH optimum of the crude enzyme was 6–7 with a temperature optimum of 80 C. The enzyme showed high pH and thermal stability retaining 100% activity at 60 C, pH 8 and 9 after 2.5 h of incubation. The residual activity at 70 C after 2.5 h was 62% and 45% at pH 8 and 9, respectively. At 75 C only 22.2% activity remained at pH 8 after 1 h incubation. Since Kraft pulp is alkaline this enzyme could be used for prebleaching of pulp at temperatures up to 70 C without pH adjustment.     

Analysis of the dynamic properties of Bacillus circulans xylanase upon formation of a covalent glycosyl-enzyme intermediate

NMR spectroscopy was used to search for mechanistically significant differences in the local mobility of the main-chain amides of Bacillus circulans xylanase (BCX) in its native and catalytically competent covalent glycosyl-enzyme intermediate states. 15N T1, T2, and 15N[1H] NOE values were measured for approximately 120 out of 178 peptide groups in both the apo form of the protein and in BCX covalently modified at position Glu78 with a mechanism-based 2-deoxy-2-fluoro-beta-xylobioside inactivator. Employing the model-free formalism of Lipari and Szabo, the measured relaxation parameters were used to calculate a global correlation time (tau(m)) for the protein in each form (9.2 +/- 0.2 ns for apo-BCX; 9.8 +/- 0.3 ns for the modified protein), as well as individual order parameters for the main-chain NH bond vectors. Average values of the order parameters for the protein in the apo and complexed forms were S2 = 0.86 +/- 0.04 and S2 = 0.91 +/- 0.04, respectively. No correlation is observed between these order parameters and the secondary structure, solvent accessibility, or hydrogen bonding patterns of amides in either form of the protein. These results demonstrate that the backbone of BCX is well ordered in both states and that formation of the glycosyl-enzyme intermediate leads to little change, in any, in the dynamic properties of BCX on the time scales sampled by 15N-NMR relaxation measurements.     

Identification of catalytic amino acids in the human GTP fucose pyrophosphorylase active site

GTP-l-fucose pyrophosphorylase(GFPP) catalyzes the reversible formation of the nucleotide-sugar GDP-beta-l-fucose from guanosine triphosphate and beta-l-fucose-1-phosphate. The enzyme functions primarily in the mammalian liver and kidney to salvage free fucose during the breakdown of glycoproteins and glycolipids. GFPP shares little primary sequence identity with other nucleotide-sugar metabolizing enzymes, and the three-dimensional structure of the protein is unknown. The enzyme does contain several sequences that could be nucleotide binding sites, but none of them are an exact match to consensus sequences. Using a combination of site-directed mutagenesis and UV photoaffinity cross-linking, we have identified five amino acid residues that are critical for catalysis. Some of these amino acids are found within the poorly conserved nucleotide binding consensus structures, while others represent new motifs. Two active site lysines can be cross-linked to photoaffinity probes. The site of cross-linking depends on the probe used. The identification of these critical residues highlights how distinct GFPP is from other nucleotide-sugar pyrophosphorylases.     

Detection and characterization of serine and threonine hydroxyl protons in Bacillus circulans xylanase by NMR spectroscopy

Hydroxyl protons on serine and threonine residues are not well characterized in protein structures determined by both NMR spectroscopy and X-ray crystallography. In the case of NMR spectroscopy, this is in large part because hydroxyl proton signals are usually hidden under crowded regions of ¹H-NMR spectra and remain undetected by conventional heteronuclear correlation approaches that rely on strong one-bond ¹H–¹⁵N or ¹H–¹³C couplings. However, by filtering against protons directly bonded to ¹³C or ¹⁵N nuclei, signals from slowly-exchanging hydroxyls can be observed in the ¹H-NMR spectrum of a uniformly ¹³C/¹⁵N-labeled protein. Here we demonstrate the use of a simple selective labeling scheme in combination with long-range heteronuclear scalar correlation experiments as an easy and relatively inexpensive way to detect and assign these hydroxyl proton signals. Using auxtrophic Escherichia coli strains, we produced Bacillus circulans xylanase (BcX) labeled with ¹³C/¹⁵N-serine or ¹³C/¹⁵N-threonine. Signals from two serine and three threonine hydroxyls in these protein samples were readily observed via ³J_C–OH couplings in long-range ¹³C-HSQC spectra. These scalar couplings (~5–7 Hz) were measured in a sample of uniformly ¹³C/¹⁵N-labeled BcX using a quantitative ¹³C/¹⁵N-filtered spin-echo difference experiment. In a similar approach, the threonine and serine hydroxyl hydrogen exchange kinetics were measured using a ¹³C/¹⁵N-filtered CLEANEX-PM pulse sequence. Collectively, these experiments provide insights into the structural and dynamic properties of several serine and threonine hydroxyls within this model protein.     

Thermostabilization of Bacillus circulans xylanase: computational optimization of unstable residues based on thermal fluctuation analysis

Low thermostability often hampers the applications of xylanases in industrial processes operated at high temperature, such as degradation of biomass or pulp bleaching. Thermostability of enzymes can be improved by the optimization of unstable residues via protein engineering. In this study, computational modeling instead of random mutagenesis was used to optimize unstable residues of Bacillus circulans xylanase (Bcx). The thermal fluctuations of unstable residues known as important to the thermal unfolding of Bcx were investigated by the molecular dynamics (MD) simulations at 300 K and 330 K to identify promising residues. The N52 site in unstable regions showed the highest thermal fluctuations. Subsequently, computational design was conducted to predict the optimal sequences of unstable residues. Five optimal single mutants were predicted by the computational design, and the N52Y mutation showed the thermostabilization effect. The N52 residue is conserved in Bacillus species xylanases and the structure analysis revealed that the N52Y mutation introduced more hydrophobic clusters for thermostability, as well as a more favorable aromatic stacking environment for substrate binding. We confirm that flexible residues at high temperature in unstable regions can be promising targets to improve thermostability of enzymes.