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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

Hyperexpression of a Bacillus circulans xylanase gene in Escherichia coli and characterization of the gene product

A 4.0-kilobase (kb) fragment of Bacillus circulans genomic DNA inserted into pUC19 and encoding endoxylanase activity was subjected to a series of subclonings. A 1.0-kb HindIII-HincII subfragment was found to code for xylanase activity. Maximum expression levels were observed with a subclone that contained an additional 0.3-kb sequence upstream from the coding region. Enhancer sequences in the upstream region are thought to be responsible for these high expression levels. Southern hybridization analyses revealed that the cloned gene hybridized with genomic DNA from Bacillus subtilis and Bacillus polymyxa. Xylanase activity expressed by Escherichia coli harboring the cloned gene was located primarily in the intracellular fraction. Levels of up to 7 U/ml or 35 mg/liter were obtained. The protein product was purified by ion exchange and gel permeation chromatography. The xylanase had a molecular weight of 20,500 and an isoelectric point of 9.0.     

Complete measurement of the pKa values of the carboxyl and imidazole groups in Bacillus circulans xylanase.

Electrostatic interactions in proteins can be dissected experimentally by determining the pKa values of their constituent ionizable amino acids. To complement previous studies of the glutamic acid and histidine residues in Bacillus circulans xylanase (BCX), we have used NMR methods to measure the pKa s of the seven aspartic acids and the C-terminus of this protein. The pKa s of these carboxyls are all less than the corresponding values observed with random coil polypeptides, indicating that their ionization contributes favorably to the stability of the folded enzyme. In general, the aspartic acids with the most reduced pKa s are those with limited exposure to the solvent and a high degree of conservation among homologous xylanases. Most dramatically, Asp 83 and Asp 101 have pKa s < 2 and thus remain deprotonated in native BCX under all conditions examined. Asp 83 is completely buried, forming a strong salt bridge with Arg 136. In contrast, Asp 101 is located on the surface of the protein, stabilized in the deprotonated form by an extensive network of hydrogen bonds involving an internal water molecule and the neutral side-chain and main-chain atoms of Ser 100 and Thr 145. These data provide a complete experimental database for theoretical studies of the ionization behavior of BCX under acidic conditions.     

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.     

Hyperexpression of a Bacillus circulans xylanase gene in Escherichia coli and characterization of the gene product.

A 4.0-kilobase (kb) fragment of Bacillus circulans genomic DNA inserted into pUC19 and encoding endoxylanase activity was subjected to a series of subclonings. A 1.0-kb HindIII-HincII subfragment was found to code for xylanase activity. Maximum expression levels were observed with a subclone that contained an additional 0.3-kb sequence upstream from the coding region. Enhancer sequences in the upstream region are thought to be responsible for these high expression levels. Southern hybridization analyses revealed that the cloned gene hybridized with genomic DNA from Bacillus subtilis and Bacillus polymyxa. Xylanase activity expressed by Escherichia coli harboring the cloned gene was located primarily in the intracellular fraction. Levels of up to 7 U/ml or 35 mg/liter were obtained. The protein product was purified by ion exchange and gel permeation chromatography. The xylanase had a molecular weight of 20,500 and an isoelectric point of 9.0.     

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.     

Thermal stabilization of Bacillus subtilis family-11 xylanase by directed evolution

We used directed evolution to enhance the thermostability of glycosyl hydrolase family-11 xylanase from Bacillus subtilis. By combining random point mutagenesis, saturation mutagenesis, and DNA shuffling, a thermostable variant, Xyl(st), was identified which contained three amino acid substitutions: Q7H, N8F, and S179C. The half-inactivation temperature (the midpoint of the melting curves) for the Xyl(st) variant compared with the wild-type enzyme after incubation for 10 min was elevated from 58 to 68 degrees C. At 60 degrees C the wild-type enzyme was inactivated within 5 min, but Xyl(st) retained full activity for at least 2 h. The stabilization was accompanied by evidence of thermophilicity; that is, an increase in the optimal reaction temperature from 55 to 65 degrees C and lower activity at low temperatures and higher activity at higher temperatures relative to wild type. To elucidate the mechanism of thermal stabilization, three-dimensional structures were determined for the wild-type and Xyl(st) enzymes. A cavity was identified around Gln-7/Asn-8 in wild type that was filled with bulky, hydrophobic residues in Xyl(st). This site was not identified by previous approaches, but directed evolution identified the region as a weak point. Formation of an intermolecular disulfide bridge via Cys-179 was observed between monomers in Xyl(st). However, the stability was essentially the same in the presence and absence of a reducing agent, indicating that the increased hydrophobicity around the Cys-179 accounted for the stability.     

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.     

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.     

Scan-rate dependence in protein calorimetry: the reversible transitions of Bacillus circulans xylanase and a disulfide-bridge mutant.

The stabilities of Bacillus circulans xylanase and a disulfide-bridge-containing mutant (S100C/N148C) were investigated by differential scanning calorimetry (DSC) and thermal inactivation kinetics. The thermal denaturation of both proteins was found to be irreversible, and the apparent transition temperatures showed a considerable dependence upon scanning rate. In the presence of low (nondenaturing) concentrations of urea, calorimetric transitions were observed for both proteins in the second heating cycle, indicating reversible denaturation occurs under those conditions. However, even for these reversible processes, the DSC curves for the wild-type protein showed a scan-rate dependence that was similar to that in the absence of urea. Calorimetric thermograms for the disulfide mutant were significantly less scan-rate dependent in the presence of urea than in the urea-free buffer. The present data show that, just as for irreversible transitions, the apparent transition temperature for the reversible denaturation of proteins can be scan-rate dependent, confirming the prediction of Lepock et al. (Lepock JR, Rithcie KP, Kolios MC, Rodahl AM, Heinz KA, Kruuf J, 1992, Biochemistry 31:12706-12712). The kinetic factors responsible for scan-rate dependence may lead to significant distortions and asymmetry of endotherms, especially at higher scanning rates. This points to the need to check for scan-rate dependence, even in the case of reversible denaturation, before any attempt is made to analyze asymmetric DSC curves by standard thermodynamic procedures. Experiments with the disulfide-bridge-containing mutant indicate that the introduction of the disulfide bond provides additional stabilization of xylanase by changing the rate-limiting step on the thermal denaturation pathway.