Crystal structures and properties of de novo circularly permuted 1,3-1,4-beta-glucanases

Fusion then fission: splitting and reassembly of an artificial fusion-protein nanocage

A split-protein system is a simple approach to introduce new termini which are useful as modification sites in protein engineering, but has been adapted mainly for monomeric proteins. Here we demonstrate the design of split subunits of the 60-mer artificial fusion-protein nanocage TIP60. The subunit fragments successfully reformed the cage structure in the same manner as prior to splitting. One of the newly introduced terminals at the interior surface can be modified using a tag peptide and green fluorescent protein. Therefore, the termini could serve as a versatile modification site for incorporating a wide variety of functional peptides and proteins.

SeqCP: A sequence-based algorithm for searching circularly permuted proteins

Circular permutation (CP) is a protein sequence rearrangement in which the amino- and carboxyl-termini of a protein can be created in different positions along the imaginary circularized sequence. Circularly permutated proteins usually exhibit conserved three-dimensional structures and functions. By comparing the structures of circular permutants (CPMs), protein research and bioengineering applications can be approached in ways that are difficult to achieve by traditional mutagenesis. Most current CP detection algorithms depend on structural information. Because there is a vast number of proteins with unknown structures, many CP pairs may remain unidentified. An efficient sequence-based CP detector will help identify more CP pairs and advance many protein studies. For instance, some hypothetical proteins may have CPMs with known functions and structures that are informative for functional annotation, but existing structure-based CP search methods cannot be applied when those hypothetical proteins lack structural information. Despite the considerable potential for applications, sequence-based CP search methods have not been well developed. We present a sequence-based method, SeqCP, which analyzes normal and duplicated sequence alignments to identify CPMs and determine candidate CP sites for proteins. SeqCP was trained by data obtained from the Circular Permutation Database and tested with nonredundant datasets from the Protein Data Bank. It shows high reliability in CP identification and achieves an AUC of 0.9. SeqCP has been implemented into a web server available at: http://pcnas.life.nthu.edu.tw/SeqCP/.

Preliminary X-ray diffraction analysis of a thermophilic β-1,3-1,4-glucanase from Clostridium thermocellum

β-1,3-1,4-Glucanases catalyze the specific hydrolysis of internal β-1,4-glycosidic bonds adjacent to the 3-O-substituted glucose residues in mixed-linked β-glucans. The thermophilic glycoside hydrolase CtGlu16A from Clostridium thermocellum exhibits superior thermal profiles, high specific activity and broad pH adaptability. Here, the catalytic domain of CtGlu16A was expressed in Escherichia coli, purified and crystallized in the trigonal space group P3121, with unit-cell parameters a=b=74.5, c=182.9 Å, by the sitting-drop vapour-diffusion method and diffracted to 1.95 Å resolution. The crystal contains two protein molecules in an asymmetric unit. Further structural determination and refinement are in progress.

Improved biocatalysts from a synthetic circular permutation library of the flavin-dependent oxidoreductase old yellow enzyme

Members of the old yellow enzyme (OYE) family are widely used, effective biocatalysts for the stereoselective trans-hydrogenation of activated alkenes. To further expand their substrate scope and improve catalytic performance, we have applied a protein engineering strategy called circular permutation (CP) to enhance the function of OYE1 from Saccharomyces pastorianus. CP can influence a biocatalyst's function by altering protein backbone flexibility and active site accessibility, both critical performance features because the catalytic cycle for OYE1 is thought to involve rate-limiting conformational changes. To explore the impact of CP throughout the OYE1 protein sequence, we implemented a highly efficient approach for cell-free cpOYE library preparation by combining whole-gene synthesis with in vitro transcription/translation. The versatility of such an ex vivo system was further demonstrated by the rapid and reliable functional evaluation of library members under variable environmental conditions with three reference substrates ketoisophorone, cinnamaldehyde, and (S)-carvone. Library analysis identified over 70 functional OYE1 variants with several biocatalysts exhibiting over an order of magnitude improved catalytic activity. Although catalytic gains of individual cpOYE library members vary by substrate, the locations of new protein termini in functional variants for all tested substates fall within the same four distinct loop/lid regions near the active site. Our findings demonstrate the importance of these structural elements in enzyme function and support the hypothesis of conformational flexibility as a limiting factor for catalysis in wild type OYE.

Thermal stabilization of an endoglucanase by cyclization

An intein-driven protein splicing approach allowed for the covalent linkage between the N- and C-termini of a polypeptide chain to create circular variants of the endo-β-1,3-1,4-glucanase, LicA, from Bacillus licheniformis. Two circular variants, LicA-C1 and LicA-C2, which have connecting loops of 20 and 14 amino acids, respectively, showed catalytic activities that are approximately two and three times higher, respectively, compared to that of the linear LicA (LicA-L1). The thermal stability of the circular variants was significantly increased compared to the linear form. Whereas the linear glucanase lost half of its activity after 3 min at 65 °C, the two circular variants have 6-fold (LicA-C1) and 16-fold (LicA-C2) increased half-life time of inactivation. In agreement with this, fluorescence spectroscopy and differential scanning calorimetry studies revealed that circular enzymes undergo structural changes at higher temperatures compared to that of the linear form. The effect of calcium on the conformational stability and function of the circular LicAs was also investigated, and we observed that the presence of calcium ions results in increased thermal stability. The impact of the length of the designed loops on thermal stability of the circular proteins is discussed, and it is suggested that cyclization may be an efficient strategy for the increased stability of proteins.

Structural redesign of lipase B from Candida antarctica by circular permutation and incremental truncation

Circular permutation of Candida antarctica lipase B yields several enzyme variants with substantially increased catalytic activity. To better understand the structural and functional consequences of protein termini reorganization, we have applied protein engineering and x-ray crystallography to cp283, one of the most active hydrolase variants. Our initial investigation has focused on the role of an extended surface loop, created by linking the native N- and C-termini, on protein integrity. Incremental truncation of the loop partially compensates for observed losses in secondary structure and the permutants' temperature of unfolding. Unexpectedly, the improvements are accompanied by quaternary-structure changes from monomer to dimer. The crystal structures of one truncated variant (cp283 Delta 7) in the apo-form determined at 1.49 A resolution and with a bound phosphonate inhibitor at 1.69 A resolution confirmed the formation of a homodimer by swapping of the enzyme's 35-residue N-terminal region. Separately, the new protein termini at amino acid positions 282/283 convert the narrow access tunnel to the catalytic triad into a broad crevice for accelerated substrate entry and product exit while preserving the native active-site topology for optimal catalytic turnover.

Random amino acid mutations and protein misfolding lead to Shannon limit in sequence-structure communication

The transmission of genomic information from coding sequence to protein structure during protein synthesis is subject to stochastic errors. To analyze transmission limits in the presence of spurious errors, Shannon's noisy channel theorem is applied to a communication channel between amino acid sequences and their structures established from a large-scale statistical analysis of protein atomic coordinates. While Shannon's theorem confirms that in close to native conformations information is transmitted with limited error probability, additional random errors in sequence (amino acid substitutions) and in structure (structural defects) trigger a decrease in communication capacity toward a Shannon limit at 0.010 bits per amino acid symbol at which communication breaks down. In several controls, simulated error rates above a critical threshold and models of unfolded structures always produce capacities below this limiting value. Thus an essential biological system can be realistically modeled as a digital communication channel that is (a) sensitive to random errors and (b) restricted by a Shannon error limit. This forms a novel basis for predictions consistent with observed rates of defective ribosomal products during protein synthesis, and with the estimated excess of mutual information in protein contact potentials.

A comprehensive analysis of non-sequential alignments between all protein structures

Background

The majority of relations between proteins can be represented as a conventional sequential alignment. Nevertheless, unusual non-sequential alignments with different connectivity of the aligned fragments in compared proteins have been reported by many researchers. It is interesting to understand those non-sequential alignments; are they unique, sporadic cases or they occur frequently; do they belong to a few specific folds or spread among many different folds, as a common feature of protein structure. We present here a comprehensive large-scale study of non-sequential alignments between available protein structures in Protein Data Bank.

Results

The study has been conducted on a non-redundant set of 8,865 protein structures aligned with the aid of the TOPOFIT method. It has been estimated that between 17.4% and 35.2% of all alignments are non-sequential depending on variations in the parameters. Analysis of the data revealed that non-sequential relations between proteins do occur systematically and in large quantities. Various sizes and numbers of non-sequential fragments have been observed with all possible complexities of fragment rearrangements found for alignments consisting of up to 12 fragments. It has been found that non-sequential alignments are not limited to proteins of any particular fold and are present in more than two hundred of them. Moreover, many of them are found between proteins with different fold assignments. It has been shown that protein structure symmetry does not explain non-sequential alignments. Therefore, compelling evidences have been provided that non-sequential alignments between proteins are systematic and widespread across the protein universe.

Conclusion

The phenomenon of the widespread occurrence of non-sequential alignments between proteins might represent a missing rule of protein structure organization. More detailed study of this phenomenon will enhance our understanding of protein stability, folding, and evolution.

The relationship between chain connectivity and domain stability in the equilibrium and kinetic folding mechanisms of dihydrofolate reductase from E.coli

The role of domains in defining the equilibrium and kinetic folding properties of dihydrofolate reductase (DHFR) from Escherichia coli was probed by examining the thermodynamic and kinetic properties of a set of variants in which the chain connectivity in the discontinuous loop domain (DLD) and the adenosine-binding domain (ABD) was altered by permutation. To test the concept that chain cleavage can selectively destabilize the domain in which the N- and C-termini are resident, permutations were introduced at one position within the ABD, one within the DLD and one at a boundary between the domains. The results demonstrated that a continuous ABD is required for a stable thermal intermediate and a continuous DLD is required for a stable urea intermediate. The permutation at the domain interface had both a thermal and urea intermediate. Strikingly, the observable kinetic folding responses of all three permuted proteins were very similar to the wild-type protein. These results demonstrate a crucial role for stable domains in defining the energy surface for the equilibrium folding reaction of DHFR. If domain connectivity affects the kinetic mechanism, the effects must occur in the sub-millisecond time range.

Circular permutation as a tool to reduce surface entropy triggers crystallization of the signal recognition particle receptor beta subunit

The production of diffraction-quality crystals remains a difficult obstacle on the road to high-resolution structural characterization of proteins. This is primarily a result of the empirical nature of the process. Although crystallization is not predictable, factors inhibiting it are well established. First, crystal formation is always entropically unfavorable. Reducing the entropic cost of crystallizing a given protein is thus desirable. It is common practice to map boundaries and remove unstructured regions surrounding the folded protein domain. However, a problem arises when flexible regions are not at the boundaries but within a domain. Such regions cannot be deleted without adding new restraints to the domain. We encountered this problem during an attempt to crystallize the beta subunit of the eukaryotic signal recognition particle (SRbeta), bearing a long and flexible internal loop. Native SRbeta did not crystallize. However, after circularly permuting the protein by connecting the spatially close N and C termini with a short heptapeptide linker GGGSGGG and removing 26 highly flexible loop residues within the domain, we obtained diffraction-quality crystals. This protein-engineering method is simple and should be applicable to other proteins, especially because N and C termini of protein domains are often close in space. The success of this method profits from prior knowledge of the domain fold, which is becoming increasingly common in today's postgenomic era.

Effect of linker sequence on the stability of circularly permuted variants of ribonuclease T1

Circularly permuted variants of ribonuclease T1 were constructed with a library of residues covalently linking the original amino and carboxyl terminal ends of the wild-type protein. The library of linking peptides consisted of three amino acids containing any combination of proline, aspartate, asparagine, serine, threonine, tyrosine, alanine, and histidine. Forty two unique linker sequences were isolated and 10 of these mutants were further characterized with regard to catalytic activity and overall thermodynamic stability. The 10 mutants with the different linking sequences (HPD, TPH, DTD, TPD, PYH, PAT, PHP, DSS, SPP, and TPS), in addition to GGG and GPG, were 4.0-6.2 kcal/mol less stable than the wild-type ribonuclease T1. However, these circular permuted variants were only 0.4-2.6 kcal/mol less stable than the direct parent protein that is missing the disulfide bond connecting residues 2 and 10. The most stable linking peptide was HPD.

Crystal structure of a natural circularly permuted jellyroll protein: 1,3-1,4-beta-D-glucanase from Fibrobacter succinogenes

The 1,3-1,4-beta-D-glucanase from Fibrobacter succinogenes (Fsbeta-glucanase) is classified as one of the family 16 glycosyl hydrolases. It hydrolyzes the glycosidic bond in the mixed-linked glucans containing beta-1,3- and beta-1,4-glycosidic linkages. We constructed a truncated form of recombinant Fsbeta-glucanase containing the catalytic domain from amino acid residues 1-258, which exhibited a higher thermal stability and enzymatic activity than the full-length enzyme. The crystal structure of the truncated Fsbeta-glucanase was solved at a resolution of 1.7A by the multiple wavelength anomalous dispersion (MAD) method using the anomalous signals from the seleno-methionine-labeled protein. The overall topology of the truncated Fsbeta-glucanase consists mainly of two eight-stranded anti-parallel beta-sheets arranged in a jellyroll beta-sandwich, similar to the fold of many glycosyl hydrolases and carbohydrate-binding modules. Sequence comparison with other bacterial glucanases showed that Fsbeta-glucanase is the only naturally occurring circularly permuted beta-glucanase with reversed sequences. Structural comparison shows that the engineered circular-permuted Bacillus enzymes are more similar to their parent enzymes with which they share approximately 70% sequence identity, than to the naturally occurring Fsbeta-glucanase of similar topology with 30% identity. This result suggests that protein structure relies more on sequence identity than topology. The high-resolution structure of Fsbeta-glucanase provides a structural rationale for the different activities obtained from a series of mutant glucanases and a basis for the development of engineered enzymes with increased activity and structural stability.

Testing the relationship between foldability and the early folding events of dihydrofolate reductase from Escherichia coli

A "folding element" is a contiguous peptide segment crucial for a protein to be foldable and is a new concept that could assist in our understanding of the protein-folding problem. It is known that the presence of the complete set of folding elements of dihydrofolate reductase (DHFR) from Escherichia coli is essential for the protein to be foldable. Since almost all of the amino acid residues known to be involved in the early folding events of DHFR are located within the folding elements, a close relationship between the folding elements and early folding events is hypothesized. In order to test this hypothesis, we have investigated whether or not the early folding events are preserved in circular permutants and topological mutants of DHFR, in which the order of the folding elements is changed but the complete set of folding elements is present. The stopped-flow circular dichroism (CD) measurements show that the CD spectra at the early stages of folding are similar among the mutants and the wild-type DHFR, indicating that the presence of the complete set of folding elements is sufficient to preserve the early folding events. We have further examined whether or not sequence perturbation on the folding elements by a single amino acid substitution affects the early folding events of DHFR. The results show that the amino acid substitutions inside of the folding elements can affect the burst-phase CD spectra, whereas the substitutions outside do not. Taken together, these results indicate that the above hypothesis is true, suggesting a close relationship between the foldability of a protein and the early folding events. We propose that the folding elements interact with each other and coalesce to form a productive intermediate(s) early in the folding, and these early folding events are important for a protein to be foldable.

Intersubunit circular permutation of human hemoglobin

For many years, human hemoglobin (Hb) isolated from erythrocytes has been investigated as a potential oxygen delivery therapeutic. Advantages with respect to the need for blood typing were balanced with various undesirable properties of cell-free Hb, including cost, overall oxygen affinity, alterations in cooperativity, and ready dissociation into toxic dimeric species. The use of total gene synthesis has resulted in very high levels of functional human Hb expression in Escherichia coli, but there remains a desire for effecting the crosslinking of the hemoglobin tetramer and providing for ready means for increasing the globular molecular weight. In this communication, we report a novel method for linking alpha chains. By circularly permuting one alpha sequence, the second alpha chain in the Hb tetramer can be linked with glycine residues to form 2 bridges across the central cavity. The second alpha chain thus presents its amino and carboxyl termini on a solvent exposed surface, providing for additional polymerization of oxygen-carrying subunits or attachment of any other peptide-based therapeutic.

Four-helix bundle topology re-engineered: monomeric Rop protein variants with different loop arrangements

We converted the small homodimeric four-helix bundle repressor of primer protein (Rop) into a monomeric four-helix bundle by introduction of connecting loops. Both left- and right-handed four-helix bundles were produced. The left-handed bundles were more stable and were used to introduce biologically interesting peptides in one of the loops.

Fold change in evolution of protein structures

Typically, protein spatial structures are more conserved in evolution than amino acid sequences. However, the recent explosion of sequence and structure information accompanied by the development of powerful computational methods led to the accumulation of examples of homologous proteins with globally distinct structures. Significant sequence conservation, local structural resemblance, and functional similarity strongly indicate evolutionary relationships between these proteins despite pronounced structural differences at the fold level. Several mechanisms such as insertions/deletions/substitutions, circular permutations, and rearrangements in beta-sheet topologies account for the majority of detected structural irregularities. The existence of evolutionarily related proteins that possess different folds brings new challenges to the homology modeling techniques and the structure classification strategies and offers new opportunities for protein design in experimental studies.

Testing the role of chain connectivity on the stability and structure of dihydrofolate reductase from E. coli: fragment complementation and circular permutation reveal stable, alternatively folded forms

The effects of chain cleavage and circular permutation on the structure, stability, and activity of dihydrofolate reductase (DHFR) from Escherichia coli were investigated by various spectroscopic and biochemical methods. Cleavage of the backbone after position 86 resulted in two fragments, (1--86) and (87--159) each of which are poorly structured and enzymatically inactive. When combined in a 1 : 1 molar ratio, however, the fragments formed a high-affinity (K(a) = 2.6 x 10(7) M(-1)) complex that displays a weakly cooperative urea-induced unfolding transition at micromolar concentrations. The retention of about 15% of the enzymatic activity of full-length DHFR is surprising, considering that the secondary structure in the complex is substantially reduced from its wild-type counterpart. In contrast, a circularly permuted form with its N-terminus at position 86 has similar overall stability to full-length DHFR, about 50% of its activity, substantial secondary structure, altered side-chain packing in the adenosine binding domain, and unfolds via an equilibrium intermediate not observed in the wild-type protein. After addition of ligand or the tight-binding inhibitor methotrexate, both the fragment complex and the circular permutant adopt more native-like secondary and tertiary structures. These results show that changes in the backbone connectivity can produce alternatively folded forms and highlight the importance of protein-ligand interactions in stabilizing the active site architecture of DHFR.

Bacterial 1,3-1,4-beta-glucanases: structure, function and protein engineering

1,3-1,4-beta-Glucanases (or lichenases, EC 3.2.1.73) hydrolyse linear beta-glucans containing beta-1,3 and beta-1,4 linkages such as cereal beta-glucans and lichenan, with a strict cleavage specificity for beta-1,4 glycosidic bonds on 3-O-substituted glucosyl residues. The bacterial enzymes are retaining glycosyl hydrolases of family 16 with a jellyroll beta-sandwich fold and a substrate binding cleft composed of six subsites. The present paper reviews the structure-function aspects of the enzymatic action including mechanistic enzymology, protein engineering and X-ray crystallographic studies.

Characterization of a beta-1,3-glucanase encoded by chlorella virus PBCV-1

Sequence analysis of the 330-kb chlorella virus PBCV-1 genome revealed an open-reading frame, A94L, that encodes a protein with significant amino acid identity to Glycoside Hydrolase Family 16 beta-1,3-glucanases. The a94l gene was cloned and the protein was expressed as a GST-A94L fusion protein in Escherichia coli. The recombinant A94L protein hydrolyzed the beta-1,3-glucose polymer laminarin and had slightly less hydrolytic activity on beta-1,3-1, 4-glucose polymers, lichenan and barley beta-glucan. The recombinant enzyme had the highest activity at 65 degrees C and pH 8. We predicted that the a94l-encoded beta-1,3-glucanase is involved in degrading the host cell wall either during virus release and/or is packaged in the virion particle and involved in virus entry. Therefore, we expected a94l to be expressed late in virus infection. However, contrary to expectations, both the a94l mRNA and the A94L protein appeared 15 min after PBCV-1 infection and disappeared 60- and 120-min p.i. postinfection, respectively, indicating that a94l is an early gene. Twenty-seven of 42 chlorella viruses contained the a94l gene. To our knowledge, this is the first report of a virus-encoded beta-1,3-glucanase.

Changing the surface of a virus shell fusion of an enzyme to polyoma VP1

Recent developments on virus-like particles have demonstrated their potential in transfecting eucaryotic cells. In the case of particles based on the major coat protein VP1 of polyoma virus, transfection occurs via binding of VP1 to sialic acids. Since sialic acid is present on almost every eucaryotic cell line, this results in an unspecific cell targeting. Generation of a cell-type specificity of this system would imply the presentation of a new function on the surface of VP1. To analyze whether a new functional protein can be placed on VP1, we inserted dihydrofolate reductase from Escherichia coli as a model protein. The effect of such an insertion on both VP1 and the inserted protein was investigated, respectively. The function of VP1, like the formation of pentameric capsomers and its ability to assemble into capsids, was not influenced by the insertion. The inserted dihydrofolate reductase showed major changes when compared to the wild-type form. The thermal stability of the enzyme was dramatically reduced in the fusion protein; nevertheless, the dihydrofolate reductase proved to be a fully active enzyme with only slightly increased K(M) values for its substrates. This model system provides the basis for further modifications of the VP1 protein to achieve an altered surface of VP1 with new properties.

Protein redesign

It has been demonstrated that, for a number of proteins, it is possible to dramatically alter the connectivities between elements of secondary structure. Remarkably large loop insertions are tolerated and many redesigns have generated proteins that successfully fold to stable, active structures. Some redesigns have been entirely the choice of the investigators, whereas others have incorporated a randomization and selection step to identify optimal sequences. These studies have provided basic guidelines for the rational manipulation of protein structure and stability, they have allowed the dissection of folding pathways and they have generated proteins with the potential for practical therapeutic applications.

A simple model for evolution of proteins towards the global minimum of free energy

BACKGROUND

Proteins seem to have their native structure in a global minimum of free energy. No mechanism is known, however, for ensuring this property. Furthermore, computational complexity studies suggest that such a mechanism is not feasible. These seemingly contradictory observations can be reconciled by the suggestion that evolutionary selection can yield proteins whose native conformation is in the global minimum of free energy. The aim of this study is to investigate such evolutionary processes in a simple model of protein folding.

RESULTS

Three possible evolutionary processes are explored.First, if the free energy of the chain is kept below a fixed threshold there is no improvement towards the global minimum. Second, if free energy is minimized directly, sequences emerge whose native conformation is in the global minimum of free energy. Third, when evolutionary pressure is applied within a small set of close homologs, sequences emerge whose functional conformation is in the global minimum of free energy.

CONCLUSIONS

Although minimizing free energy does select for sequences whose functional conformation is in the global free energy minimum, we argue that for most proteins, which typically have free energy values of only 5-15 kcal/mol, such evolutionary pressure cannot be considered biologically plausible. In contrast, by repeatedly forcing sequences to avoid drifting towards competing "non-native" conformations, sequences emerge whose native conformation becomes very close to the global minimum of free energy. We argue that such a mechanism is both efficient and biologically plausible.

Structure and function of the Bacillus hybrid enzyme GluXyn-1: native-like jellyroll fold preserved after insertion of autonomous globular domain

The 1,3-1,4-beta-glucanase from Bacillus macerans (wtGLU) and the 1, 4-beta-xylanase from Bacillus subtilis (wtXYN) are both single-domain jellyroll proteins catalyzing similar enzymatic reactions. In the fusion protein GluXyn-1, the two proteins are joined by insertion of the entire XYN domain into a surface loop of cpMAC-57, a circularly permuted variant of wtGLU. GluXyn-1 was generated by protein engineering methods, produced in Escherichia coli and shown to fold spontaneously and have both enzymatic activities at wild-type level. The crystal structure of GluXyn-1 was determined at 2.1 A resolution and refined to R = 17.7% and R(free) = 22.4%. It shows nearly ideal, native-like folding of both protein domains and a small, but significant hinge bending between the domains. The active sites are independent and accessible explaining the observed enzymatic activity. Because in GluXyn-1 the complete XYN domain is inserted into the compact folding unit of GLU, the wild-type-like activity and tertiary structure of the latter proves that the folding process of GLU does not depend on intramolecular interactions that are short-ranged in the sequence. Insertion fusions of the GluXyn-1 type may prove to be an easy route toward more stable bifunctional proteins in which the two parts are more closely associated than in linear end-to-end protein fusions.