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.     

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.     

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.     

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.     

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.     

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.     

Structural disulfide bonds in the Bacillus thuringiensis subsp. israelensis protein crystal.

We examined disulfide bonds in mosquito larvicidal crystals produced by Bacillus thuringiensis subsp. israelensis. Intact crystals contained 2.01 X 10(-8) mol of free sulfhydryls and 3.24 X 10(-8) mol of disulfides per mg of protein. Reduced samples of alkali-solubilized crystals resolved into several proteins, the most prominent having apparent molecular sizes of 28, 70, 135, and 140 kilodaltons (kDa). Nonreduced samples contained two new proteins of 52 and 26 kDa. When reduced, both the 52- and 26-kDa proteins were converted to 28-kDa proteins. Furthermore, both bands reacted with antiserum prepared against reduced 28-kDa protein. Approximately 50% of the crystal proteins could be solubilized without disulfide cleavage. These proteins were 70 kDa or smaller. Solubilization of the 135- and 140-kDa proteins required disulfide cleavage. Incubation of crystals at pH 12.0 for 2 h cleaved 40% of the disulfide bonds and solubilized 83% of the crystal protein. Alkali-stable disulfides were present in both the soluble and insoluble portions. The insoluble pellet contained 12 to 14 disulfides per 100 kDa of protein and was devoid of sulfhydryl groups. Alkali-solubilized proteins contained both intrachain and interchain disulfide bonds. Despite their structural significance, it is unlikely that disulfide bonds are involved in the formation or release of the larvicidal toxin.     

Growth dynamics of Bacillus circulans colony

We have investigated the growth dynamics of Bacillus circulans colony exhibiting the knotted-branching pattern by swarming on a hard agar medium. The knotted-branching pattern consists of many circular clusters, so-called subcolonies, and their trajectories. We analysed the processes of a subcolony because they are presumably the key elements for the formation of knotted-branching pattern. It was found that a subcolony has three processes, i.e. "generation", "growth", and "migration" by microscopic and time-resolved observations. An embryonic small subcolony (child subcolony) formed around an existing subcolony (parent subcolony) grows larger and migrates away from the parent subcolony. We proposed a simple model to explain the migration and the growth processes. It is assumed that the internal part of the subcolony is unfavorable for the bacteria and that the motion of the child subcolony on the agar medium can be modeled using a frictional force. The experimental data were quantitatively analysed in order to compare with models. Our models are consistent with the experimental results on following three points: (1) the radius of a subcolony increases linearly with the incubation time, (2) a subcolony stops just after formation and then starts to migrate suddenly, and (3) the trajectory of a subcolony predicted by the model agrees with the experimental one.     

Acid stabilization of Bacillus licheniformis alpha amylase through introduction of mutations

This paper provided further understanding of the relationships between acid resistance and structural features of different mutants in Bacillus licheniformis alpha amylase (BLA) due to the changes of two crucial positions Leu134 and Ser320. In order to investigate effect of the two positions on the acid stability, we described the detailed characterization of wild-type and the single mutants L134R and S320A as well as the double mutant L134R/S320A. The highest k _cat /Km with pH 4.5, approximately 14 times that of wild type, was observed in L134R/S320A. The k _cat /Km corresponding to L134R and S320A were at an intermediate values between those for wild type and L134R/S320A. In addition, compared with wild type, which had a rapid decline of the activity, L134R/S320A could maintain its activity strongly in low pH. Meanwhile, lower tolerance of L134R and S320A in acidic conditions than that of L134R/S320A was determined. Surprisingly, the acid-resistant capability of L134R/S320A was significantly enhanced by directed evolution. These results, combined with three-dimensional structure analysis, show that the electrostatic effects play a significant role in determining the stability of BLA at two crucial positions, 134 and 320.     

Profile of the disulfide bonds in acetylcholinesterase

The inter- and intrasubunit disulfide bridges for the 11 S form of acetylcholinesterase isolated from Torpedo californica have been identified. Localized within the basal lamina of the synapse, the dimensionally asymmetric forms of acetylcholinesterase contain either two (13 S) or three (17 S) sets of catalytic subunits linked to collagenous and noncollagenous structural subunits. Limited proteolysis of these molecules yields a tetramer of catalytic subunits that sediments at 11 S. Each catalytic subunit contains 8 cysteine residues which were identified following tryptic digestion of the reduced, 14C-carboxymethylated protein. The tryptic peptides were purified by gel filtration followed by reverse-phase high performance liquid chromatography (HPLC) and then sequenced. The disulfide bonding profile was determined by treating the native, nonreduced 11 S form of acetylcholinesterase with a fluorescent, sulfhydryl-specific reagent, monobromobimane, prior to tryptic digestion. Peptides again were resolved by gel filtration and reverse-phase HPLC. One fluorescent cysteine-containing peptide was identified, indicating that a single sulfhydryl residue, Cys231, was present in its reduced form. Three pairs of disulfide-bonded peptides were identified. These were localized in the polypeptide chain based on the cDNA-deduced sequence of the protein and were identified as Cys67-Cys94, Cys254-Cys265, and Cys402-Cys521. Hence, three loops are found in the secondary structure. Cys572, located in the carboxyl-terminal tryptic peptide, was disulfide-bonded to an identical peptide and most likely forms an intersubunit cross-link. Since the 6 cysteine residues in acetylcholinesterase that are involved in intrachain disulfide bonds are conserved in the sequence of the homologous protein thyroglobulin, it is likely that both proteins share a common folding pattern in their respective tertiary structures. Cys231 and the carboxyl-terminal cysteine residue Cys572 are not conserved in thyroglobulin.     

Kinetic stabilization of Bacillus licheniformis alpha-amylase through introduction of hydrophobic residues at the surface

It is generally assumed that in proteins hydrophobic residues are not favorable at solvent-exposed sites, and that amino acid substitutions on the surface have little effect on protein thermostability. Contrary to these assumptions, we have identified hyperthermostable variants of Bacillus licheniformis alpha-amylase (BLA) that result from the incorporation of hydrophobic residues at the surface. Under highly destabilizing conditions, a variant combining five stabilizing mutations unfolds 32 times more slowly and at a temperature 13 degrees C higher than the wild-type. Crystal structure analysis at 1.7 A resolution suggests that stabilization is achieved through (a) extension of the concept of increased hydrophobic packing, usually applied to cavities, to surface indentations, (b) introduction of favorable aromatic-aromatic interactions on the surface, (c) specific stabilization of intrinsic metal binding sites, and (d) stabilization of a beta-sheet by introducing a residue with high beta-sheet forming propensity. All mutated residues are involved in forming complex, cooperative interaction networks that extend from the interior of the protein to its surface and which may therefore constitute "weak points" where BLA unfolding is initiated. This might explain the unexpectedly large effect induced by some of the substitutions on the kinetic stability of BLA. Our study shows that substantial protein stabilization can be achieved by stabilizing surface positions that participate in underlying cooperatively formed substructures. At such positions, even the apparently thermodynamically unfavorable introduction of hydrophobic residues should be explored.     

Maltodextrin-dependent crystallization of cyclomaltodextrin glucanotransferase from Bacillus circulans

Crystals of cyclomaltodextrin glucanotransferase from Bacillus circulans (EC 2.4.1.19) suitable for high-resolution X-ray analysis were obtained by vapor diffusion against 60% (v/v) 2-methyl 2,4-pentanediol buffered with 100 mM-sodium Hepes, pH 7.55. The crystals have P2(1)2(1)2(1) space group symmetry, with a = 120.4 A, b = 110.9 A and c = 66.4 A, and contain one molecule of 68,000 in the asymmetric unit. Growth of single enzyme crystals was found to require the presence of either alpha-cyclodextrin, beta-cyclodextrin, gamma-cyclodextrin, or maltose in high molar excess, a requirement that could not be fulfilled by glucose, the basic building block of these compounds. Although the exact role of cyclic and linear maltodextrins in enzyme crystallization is not yet known, we have preliminary evidence that these compounds are degraded by the enzyme in the crystallization droplet.