Feature-incorporated alignment based ligand-binding residue prediction for carbohydrate-binding modules.

MOTIVATION: Carbohydrate-binding modules (CBMs) share similar secondary and tertiary topology, but their primary sequence identity is low. Computational identification of ligand-binding residues allows biologists to better understand the protein-carbohydrate binding mechanism. In general, functional characterization can be alternatively solved by alignment-based manners. As alignment accuracy based on conventional methods is often sensitive to sequence identity, low sequence identity among query sequences makes it difficult to precisely locate small portions of relevant features. Therefore, we propose a feature-incorporated alignment (FIA) to flexibly align conserved signatures in CBMs. Then, an FIA-based target-template prediction model was further implemented to identify functional ligand-binding residues. RESULTS: Arabidopsis thaliana CBM45 and CBM53 were used to validate the FIA-based prediction model. The predicted ligand-binding residues residing on the surface in the hypothetical structures were verified to be ligand-binding residues. In the absence of 3D structural information, FIA demonstrated significant improvement in the estimation of sequence similarity and identity for a total of 808 sequences from 11 different CBM families as compared with six leading tools by Friedman rank test.

CBM21 starch-binding domain: a new purification tag for recombinant protein engineering.

The use of protein fusion tag technology simplifies and facilitates purification of recombinant proteins. In this article, we have found that the starch-binding domain derived from Rhizopus oryzae glucoamylase (RoSBD), a member of carbohydrate-binding module family 21 (CBM21) with raw starch-binding activity, is favorable to be applied as an affinity tag for fusion protein engineering and purification in Escherichia coli and Pichia pastoris systems. To determine suitable spatial arrangement of RoSBD as a fusion handle, enhanced green fluorescent protein (eGFP) was fused to either the N- or C-terminus of the SBD, expressed by E. coli, and purified for yield assessment and functional analysis. Binding assays showed that the ligand-binding capacity was fully retained when the RoSBD was engineered at either the N-terminal or the C-terminal end. Similar results have been obtained with the RoSBD-conjugated phytase secreted by P. pastoris. The effective adsorption onto raw starch and low cost of starch make RoSBD practically applicable in terms of development of a new affinity fusion tag for recombinant protein engineering in an economic manner.

Crystal structures of the starch-binding domain from Rhizopus oryzae glucoamylase reveal a polysaccharide-binding path.

GA (glucoamylase) hydrolyses starch and polysaccharides to beta-D-glucose. RoGA (Rhizopus oryzae GA) consists of two functional domains, an N-terminal SBD (starch-binding domain) and a C-terminal catalytic domain, which are connected by an O-glycosylated linker. In the present study, the crystal structures of the SBD from RoGA (RoGACBM21) and the complexes with beta-cyclodextrin (SBD-betaCD) and maltoheptaose (SBD-G7) were determined. Two carbohydrate binding sites, I (Trp(47)) and II (Tyr(32)), were resolved and their binding was co-operative. Besides the hydrophobic interaction, two unique polyN loops comprising consecutive asparagine residues also participate in the sugar binding. A conformational change in Tyr(32) was observed between unliganded and liganded SBDs. To elucidate the mechanism of polysaccharide binding, a number of mutants were constructed and characterized by a quantitative binding isotherm and Scatchard analysis. A possible binding path for long-chain polysaccharides in RoGACBM21 was proposed.

Identification and characterization of a novel fibril forming peptide in fungal starch binding domain.

Scanty information is available regarding the chemical basis for structural alterations of the carbohydrate-binding modules (CBMs). The N-terminal starch binding domain (SBD) of Rhizopus oryzae glucoamylase (GA) forms fibrils under thermal stress, presenting an unusual conformational change from immunoglobulin-like to beta-sheet-rich structure. Site-directed mutagenesis revealed that the C-terminal Lys of SBD played a crucial role in the fibril formation. The synthetic peptide (DNNNSANYQVSTSK) representing the C-terminal 14 amino acid residues of SBD was further demonstrated to act as a fibril-forming segment, in which terminal charges and an internal NNNxxNYQ motif were key fibril-forming determinants. The formation of fibril structure in a fungal SBD, caused by its chemical and biophysical requirements, was demonstrated for the first time.

Solution structure of family 21 carbohydrate-binding module from Rhizopus oryzae glucoamylase.

CBMs (carbohydrate-binding modules) function independently to assist carbohydrate-active enzymes. Family 21 CBMs contain approx. 100 amino acid residues, and some members have starchbinding functions or glycogen-binding activities. We report here the first structure of a family 21 CBM from the SBD (starch-binding domain) of Rhizopus oryzae glucoamylase (RoCBM21) determined by NMR spectroscopy. This CBM has a beta-sandwich fold with an immunoglobulin-like structure. Ligand-binding properties of RoCBM21 were analysed by chemical-shift perturbations and automated docking. Structural comparisons with previously reported SBDs revealed two types of topologies, namely type I and type II, with CBM20, CBM25, CBM26 and CBM41 showing type I topology, with CBM21 and CBM34 showing type II topology. According to the chemical-shift perturbations, RoCBM21 contains two ligand-binding sites. Residues in site II are similar to those found in the family 20 CBM from Aspergillus niger glucoamylase (AnCBM20). Site I, however, is embedded in a region with unique sequence motifs only found in some members of CBM21s. Additionally, docking of beta-cyclodextrin and malto-oligosaccharides highlights that side chains of Y83 and W47 (one-letter amino acid code) form the central part of the conserved binding platform in the SBD. The structure of RoCBM21 provides the first direct evidence of the structural features and the basis for protein-carbohydrate recognition from an SBD of CBM21.

REMUS: a tool for identification of unique peptide segments as epitopes.

We provide a 'R(E)MUS' (reinforced merging techniques for unique peptide segments) web server for identification of the locations and compositions of unique peptide segments from a set of protein family sequences. Different levels of uniqueness are determined according to substitutional relationship in the amino acids, frequency of appearance and biological properties such as priority for serving as candidates for epitopes where antibodies recognize. R(E)MUS also provides interactive visualization of 3D structures for allocation and comparison of the identified unique peptide segments. Accuracy of the algorithm was found to be 70% in terms of mapping a unique peptide segment as an epitope. The R(E)MUS web server is available at http://biotools.cs.ntou.edu.tw/REMUS and the PC version software can be freely downloaded either at http://bioinfo.life.nthu.edu.tw/REMUS or http://spider.cs.ntou.edu.tw/BioTools/REMUS. User guide and working examples for PC version are available at http://spider.cs.ntou.edu.tw/BioTools/REMUS-DOCS.html, and details of the proposed algorithm can be referred to the documents as described previously [H. T. Chang, T. W. Pai, T. C. Fan, B. H. Su, P. C. Wu, C. Y. Tang, C. T. Chang, S. H. Liu and M. D. T. Chang (2006) BMC Bioinformatics, 7, 38 and T. W. Pai, B. H. Su, P. C. Wu, M. D. T. Chang, H. T. Chang, T. C. Fan and S. H. Liu (2006) J. Bioinform. Comput. Biol., 4, 75-92].

Role of the linker region in the expression of Rhizopus oryzae glucoamylase.

BACKGROUND: Rhizopus oryzae glucoamylase (RoGA) consists of three domains: an amino (N)-terminal raw starch-binding domain (SBD), a glycosylated linker domain, and a carboxy (C)-terminal catalytic domain. The 36-amino-acid linker region (residues 132-167) connects the two functional domains, but its structural and functional roles are unclear. RESULTS: To characterize the linker sequences of RoGA and its involvement in protein expression, a number of RoGA variants containing deletions and mutations were constructed and expressed in Saccharomyces cerevisiae. Deletion analyses demonstrate that the linker region, especially within residues 161 to 167, is required for protein expression. In addition, site-directed mutagenesis and deglycosylation studies reveal that the linker region of RoGA contains both N- and O-linked carbohydrate moieties, and the N-linked oligosaccharides play a major role in the formation of active enzyme. Although the linker segment itself appears to have no ordered secondary structural conformation, the flexible region indeed contributes to the stabilization of functional N- and C-terminal domains. CONCLUSION: Our data provide direct evidence that the length, composition, and glycosylation of the interdomain linker play a central role in the structure and function of RoGA.

The family 21 carbohydrate-binding module of glucoamylase from Rhizopus oryzae consists of two sites playing distinct roles in ligand binding.

The starch-hydrolysing enzyme GA (glucoamylase) from Rhizopus oryzae is a commonly used glycoside hydrolase in industry. It consists of a C-terminal catalytic domain and an N-terminal starch-binding domain, which belong to the CBM21 (carbohydrate-binding module, family 21). In the present study, a molecular model of CBM21 from R. oryzae GA (RoGACBM21) was constructed according to PSSC (progressive secondary structure correlation), modified structure-based sequence alignment, and site-directed mutagenesis was used to identify and characterize potential ligand-binding sites. Our model suggests that RoGACBM21 contains two ligand-binding sites, with Tyr32 and Tyr67 grouped into site I, and Trp47, Tyr83 and Tyr93 grouped into site II. The involvement of these aromatic residues has been validated using chemical modification, UV difference spectroscopy studies, and both qualitative and quantitative binding assays on a series of RoGACBM21 mutants. Our results further reveal that binding sites I and II play distinct roles in ligand binding, the former not only is involved in binding insoluble starch, but also facilitates the binding of RoGACBM21 to long-chain soluble polysaccharides, whereas the latter serves as the major binding site mediating the binding of both soluble polysaccharide and insoluble ligands. In the present study we have for the first time demonstrated that the key ligand-binding residues of RoGACBM21 can be identified and characterized by a combination of novel bioinformatics methodologies in the absence of resolved three-dimensional structural information.

Constrained multiple sequence alignment tool development and its application to RNase family alignment.

In this paper, we design a heuristic algorithm of computing a constrained multiple sequence alignment (CMSA for short) for guaranteeing that the generated alignment satisfies the user-specified constraints that some particular residues should be aligned together. If the number of residues needed to be aligned together is a constant alpha, then the time-complexity of our CMSA algorithm for aligning K sequences is O(alphaKn(4)), where n is the maximum of the lengths of sequences. In addition, we have built up such a CMSA software system and made several experiments on the RNase sequences, which mainly function in catalyzing the degradation of RNA molecules. The resulting alignments illustrate the practicability of our method.

Two unique ligand-binding clamps of Rhizopus oryzae starch binding domain for helical structure disruption of amylose.

The N-terminal starch binding domain of Rhizopus oryzae glucoamylase (RoSBD) has a high binding affinity for raw starch. RoSBD has two ligand-binding sites, each containing a ligand-binding clamp: a polyN clamp residing near binding site I is unique in that it is expressed in only three members of carbohydrate binding module family 21 (CBM21) members, and a Y32/F58 clamp located at binding site II is conserved in several CBMs. Here we characterized different roles of these sites in the binding of insoluble and soluble starches using an amylose-iodine complex assay, atomic force microscopy, isothermal titration calorimetry, site-directed mutagenesis, and structural bioinformatics. RoSBD induced the release of iodine from the amylose helical cavity and disrupted the helical structure of amylose type III, thereby significantly diminishing the thickness and length of the amylose type III fibrils. A point mutation in the critical ligand-binding residues of sites I and II, however, reduced both the binding affinity and amylose helix disruption. This is the first molecular model for structure disruption of the amylose helix by a non-hydrolytic CBM21 member. RoSBD apparently twists the helical amylose strands apart to expose more ligand surface for further SBD binding. Repeating the process triggers the relaxation and unwinding of amylose helices to generate thinner and shorter amylose fibrils, which are more susceptible to hydrolysis by glucoamylase. This model aids in understanding the natural roles of CBMs in protein-glycan interactions and contributes to potential molecular engineering of CBMs.

Hydrophilic aromatic residue and in silico structure for carbohydrate binding module.

Carbohydrate binding modules (CBMs) are found in polysaccharide-targeting enzymes and increase catalytic efficiency. Because only a relatively small number of CBM structures have been solved, computational modeling represents an alternative approach in conjunction with experimental assessment of CBM functionality and ligand-binding properties. An accurate target-template sequence alignment is the crucial step during homology modeling. However, low sequence identities between target/template sequences can be a major bottleneck. We therefore incorporated the predicted hydrophilic aromatic residues (HARs) and secondary structure elements into our feature-incorporated alignment (FIA) algorithm to increase CBM alignment accuracy. An alignment performance comparison for FIA and six others was made, and the greatest average sequence identities and similarities were achieved by FIA. In addition, structure models were built for 817 representative CBMs. Our models possessed the smallest average surface-potential z scores. Besides, a large true positive value for liagnd-binding aromatic residue prediction was obtained by HAR identification. Finally, the pre-simulated CBM structures have been deposited in the Database of Simulated CBM structures (DS-CBMs). The web service is publicly available at http://dscbm.life.nthu.edu.tw/ and http://dscbm.cs.ntou.edu.tw/.

Improved secretory production of glucoamylase in Pichia pastoris by combination of genetic manipulations.

Rhizopus oryzae glucoamylase (GA) has been genetically engineered with modified signal peptide (MSP), increased copy number of the gene, and coexpression of SEC4, a gene encoding a Rab protein associated with secretory vesicles, and its secretion level has been successfully raised up to 100-fold in Pichia pastoris. The MSP was designed to contain the signal peptide of mouse salivary alpha-amylase (S8L) fused to the pro-region of the signal peptide of Saccharomyces cerevisiae alpha-mating factor to replace the wild type signal peptide (WTSP) of GA. The P. pastoris transformant MSPGA-1 containing a single copy of MSPGA gene showed a 3.6-fold increase in GA secretion as compared to that of WTSPGA-1. Moreover, the P. pastoris transformant MSPGA-7 harboring seven copies of the MSPGA inserts was identified and showed 56-fold higher secreted GA than WTSPGA-1. In addition, we found that overexpression of SEC4 further doubled the secretion level of GA in each MSPGA/P. pastoris transformant. Taken together, the MSPGA-7-SEC4 clone showed as much as 100-fold secretion level of GA when compared to WTSPGA-1. In summary, we have demonstrated that combination of the aforementioned genetic manipulations resulted in high level secretion of R. oryzae GA in P. pastoris.

Molecular genetic manipulation of Pichia pastoris SEC4 governs cell growth and glucoamylase secretion.

We have previously engineered a recombinant Pichia pastoris GS115 transformant, MSPGA-7, harboring seven copies of glucoamylase (GA) fused with modified signal peptide. High yield secretion of GA was achieved as an extra copy of SEC4 was integrated to the transformant. To elucidate the physiological role of SEC4, a dominant-negative mutant of SEC4, SEC4(S28N), was overexpressed under the control of alchohol oxidase 1 (AOX1) promoter in P. pastoris strain MSPGA-7 as well as a set of host cells harboring multi-copy of wild type SEC4. We found that SEC4(S28N) mutation in the key guanine nucleotide binding domain reduced guanine nucleotide binding affinity, hence it blocked the transport of vesicles required for targeting and fusion to the plasma membrane. The inhibitory levels of cell growth and GA secretion were correlated with the dosage of SEC4(S28N) gene. In addition, overexpression of SEC4 driven by AOX1 promoter in MSPGA-7 improved the secretory production of GA, but demonstrated the delay of cell growth by increased gene dosage of SEC4. Interestingly, a limited level of Sec4p did not disturb the cell growth. It was because expression of only one copy of SEC4 resulted in delay of cell growth at an early stage while still maintaining high level Sec4p at long-term incubation. Accordingly, as glyceraldehyde-3-phosphate dehydrogenase promoter was used to substitute AOX1 promoter to drive the SEC4 expression, enhanced GA secretion but not inhibition of cell growth was achieved. Taken together, our results demonstrate that SEC4 is essential for P. pastoris in regulating cell growth and heterologous protein secretion in a dosage-dependent manner.

Constrained multiple sequence alignment tool development and its application to RNase family alignment.

In this paper, we design an algorithm of computing a constrained multiple sequence alignment (CMSA for short) for guaranteeing that the generated alignment satisfies the user-specified constraints that some particular residues should be aligned together. If the number of residues needed to be aligned together is a constant alpha, then the time-complexity of our CMSA algorithm for aligning K sequences is O (alphaKn4), where n is the maximum of the lengths of sequences. In addition, we have build up such a CMSA software system and made several experiments on the RNase sequences, which mainly function in catalyzing the degradation of RNA molecules. The resulting alignments illustrate the practicability of our method.