Purification and properties of the catalytic domain of the thermostable pullulanase type II from Thermococcus hydrothermalis

A pullulanase type II was produced in Escherichia coli using the relevant gene from Thermococcus hydrothermalis. This protein was purified and its pullulanolytic and amylolytic activities were characterised. The optimum temperature and Ca²⁺ concentration for each activity were identical (105 °C and 0.09 mM), whereas the optimum pH (pH_pullulan 5.75, pH_amylose 5) and the influence of Ca²⁺ ions on the kinetic parameters were different. Further analyses revealed that this enzyme exhibits an endo-processive-like action and specifically cleaves -1,6 bonds in pullulan.     

Xylanase XynA from the hyperthermophilic bacterium Thermotoga maritima: structure and stability of the recombinant enzyme and its isolated cellulose-binding domain.

The hyperthermophilic bacterium Thermotoga maritima is capable of gaining metabolic energy utilizing xylan. XynA, one of the corresponding hydrolases required for its degradation, is a 120-kDa endo-1,4-D-xylanase exhibiting high intrinsic stability and a temperature optimum approximately 90 degrees C. Sequence alignments with other xylanases suggest the enzyme to consist of five domains. The C-terminal part of XynA was previously shown to be responsible for cellulose binding (Winterhalter C, Heinrich P, Candussio A, Wich G, Liebl W. 1995. Identification of a novel cellulose-binding domain within the multi-domain 120 kDa Xylanase XynA of the hyperthermophilic bacterium Thermotoga maritima. Mol Microbiol 15:431-444). In order to characterize the domain organization and the stability of XynA and its C-terminal cellulose-binding domain (CBD), the two separate proteins were expressed in Escherichia coli. CBD, because of its instability in its ligand-free form, was expressed as a glutathione S-transferase fusion protein with a specific thrombin cleavage site as linker. XynA and CBD were compared regarding their hydrodynamic and spectral properties. As taken from analytical ultracentrifugation and gel permeation chromatography, both are monomers with 116 and 22 kDa molecular masses, respectively. In the presence of glucose as a ligand, CBD shows high intrinsic stability. Denaturation/renaturation experiments with isolated CBD yield > 80% renaturation, indicating that the domain folds independently. Making use of fluorescence emission and far-UV circular dichroism in order to characterize protein stability, guanidine-induced unfolding of XynA leads to biphasic transitions, with half-concentrations c1/2 (GdmCl) approximately 4 M and > 5 M, in accordance with the extreme thermal stability. At acid pH, XynA exhibits increased stability, indicated by a shift of the second guanidine-transition from 5 to 7 M GdmCl. This can be tentatively attributed to the cellulose-binding domain. Differences in the transition profiles monitored by fluorescence emission and dichroic absorption indicate multi-state behavior of XynA. In the case of CBD, a temperature-induced increase in negative ellipticity at 217 nm is caused by alterations in the environment of aromatic residues that contribute to the far-UV CD in the native state.     

The structure of protein evolution and the evolution of protein structure

The observed distribution of protein structures can give us important clues about the underlying evolutionary process, imposing important constraints on possible models. The availability of results from an increasing number of genome projects has made the development of these models an active area of research. Models explaining the observed distribution of structures have focused on the inherent functional capabilities and structural properties of different folds and on the evolutionary dynamics. Increasingly, these elements are being combined.     

Electrostatic interactions contribute to reduced heat capacity change of unfolding in a thermophilic ribosomal protein l30e

The origin of reduced heat capacity change of unfolding (DeltaC(p)) commonly observed in thermophilic proteins is controversial. The established theory that DeltaC(p) is correlated with change of solvent-accessible surface area cannot account for the large differences in DeltaC(p) observed for thermophilic and mesophilic homologous proteins, which are very similar in structures. We have determined the protein stability curves, which describe the temperature dependency of the free energy change of unfolding, for a thermophilic ribosomal protein L30e from Thermococcus celer, and its mesophilic homologue from yeast. Values of DeltaC(p), obtained by fitting the free energy change of unfolding to the Gibbs-Helmholtz equation, were 5.3 kJ mol(-1) K(-1) and 10.5 kJ mol(-1) K(-1) for T.celer and yeast L30e, respectively. We have created six charge-to-neutral mutants of T.celer L30e. Removal of charges at Glu6, Lys9, and Arg92 decreased the melting temperatures of T.celer L30e by approximately 3-9 degrees C, and the differences in melting temperatures were smaller with increasing concentration of salt. These results suggest that these mutations destabilize T.celer L30e by disrupting favorable electrostatic interactions. To determine whether electrostatic interactions contribute to the reduced DeltaC(p) of the thermophilic protein, we have determined DeltaC(p) for wild-type and mutant T.celer L30e by Gibbs-Helmholtz and by van't Hoff analyses. A concomitant increase in DeltaC(p) was observed for those charge-to-neutral mutants that destabilize T.celer L30e by removing favorable electrostatic interactions. The crystal structures of K9A, E90A, and R92A, were determined, and no structural change was observed. Taken together, our results support the conclusion that electrostatic interactions contribute to the reduced DeltaC(p) of T.celer L30e.     

Crystal structure of an archaeal DNA sliding clamp: proliferating cell nuclear antigen from Pyrococcus furiosus

The proliferating cell nuclear antigen (PCNA) is now recognized as one of the key proteins in DNA metabolic events because of its direct interactions with many proteins involved in important cellular processes. We have determined the crystal structure of PCNA from a hyperthermophilic archaeon, Pyrococcus furiosus (pfuPCNA), at 2.1 A resolution. pfuPCNA forms a toroidal, ring-shaped structure consisting of homotrimeric molecules, which is also observed in the PCNA crystals from human and yeast. The overall structure of pfuPCNA is highly conserved with other PCNA proteins, as well as with the bacterial ss clamp and the bacteriophage gp45. This result shows that the three-dimensional structure of the sliding clamp is conserved in the three domains of life. pfuPCNA has two remarkable features compared with the human and yeast PCNA molecules: it has more ion pairs and fewer intermolecular main chain hydrogen bonds. The former may contribute to the thermal stability of pfuPCNA, and the latter may be the cause of the stimulatory effect of pfuPCNA on the DNA synthesizing activity of P. furiosus DNA polymerases in the absence of the clamp loader replication factor C in vitro.     

Full-sequence computational design and solution structure of a thermostable protein variant

Computational protein design procedures were applied to the redesign of the entire sequence of a 51 amino acid residue protein, Drosophila melanogaster engrailed homeodomain. Various sequence optimization algorithms were compared and two resulting designed sequences were experimentally evaluated. The two sequences differ by 11 mutations and share 22% and 24% sequence identity with the wild-type protein. Both computationally designed proteins were considerably more stable than the naturally occurring protein, with midpoints of thermal denaturation greater than 99 degrees C. The solution structure was determined for one of the two sequences using multidimensional heteronuclear NMR spectroscopy, and the structure was found to closely match the original design template scaffold.     

Combining protein evolution and secondary structure

An evolutionary model that combines protein secondary structure and amino acid replacement is introduced. It allows likelihood analysis of aligned protein sequences and does not require the underlying secondary (or tertiary) structures of these sequences to be known. One component of the model describes the organization of secondary structure along a protein sequence and another specifies the evolutionary process for each category of secondary structure. A database of proteins with known secondary structures is used to estimate model parameters representing these two components. Phylogeny, the third component of the model, can be estimated from the data set of interest. As an example, we employ our model to analyze a set of sucrose synthase sequences. For the evolution of sucrose synthase, a parametric bootstrap approach indicates that our model is statistically preferable to one that ignores secondary structure.      

Alternative splicing and protein structure evolution

Alternative splicing is thought to be one of the major sources for functional diversity in higher eukaryotes. Interestingly, when mapping splicing events onto protein structures, about half of the events affect structured and even highly conserved regions i.e. are non-trivial on the structure level. This has led to the controversial hypothesis that such splice variants result in nonsense-mediated mRNA decay or non-functional, unstructured proteins, which do not contribute to the functional diversity of an organism. Here we show in a comprehensive study on alternative splicing that proteins appear to be much more tolerant to structural deletions, insertions and replacements than previously thought. We find literature evidence that such non-trivial splicing isoforms exhibit different functional properties compared to their native counterparts and allow for interesting regulatory patterns on the protein network level. We provide examples that splicing events may represent transitions between different folds in the protein sequence–structure space and explain these links by a common genetic mechanism. Taken together, those findings hint to a more prominent role of splicing in protein structure evolution and to a different view of phenotypic plasticity of protein structures.     

Relation between protein stability, evolution and structure, as probed by carboxylic acid mutations

Native proteins are marginally stable. Low thermodynamic stability may actually be advantageous, although the accumulation of neutral, destabilizing mutations may have also contributed to it. In any case, once marginal stability has been reached, it appears plausible that mutations at non-constrained positions become fixed in the course of evolution (due to random drift) with frequencies that roughly reflect the mutation effects on stability ("pseudo-equilibrium hypothesis"). We have found that all glutamate-->aspartate mutations in wild-type Escherichia coli thioredoxin are destabilizing, as well as most of the aspartate-->glutamate mutations. Furthermore, the effect of these mutations on thioredoxin thermodynamic stability shows a robust correlation with the frequencies of occurrence of the involved residues in several-hundred sequence alignments derived from a BLAST search. These results provide direct and quantitative experimental evidence for the pseudo-equilibrium hypothesis and should have general consequences for the interpretation of mutation effects on protein stability, as they suggest that residue environments in proteins may be optimized for stabilizing interactions to a remarkable degree of specificity. We also provide evidence that such stabilizing interactions may be detected in sequence alignments, and briefly discuss the implications of this possibility for the derivation of structural information (on native and denatured states) from comparative sequence analyses.     

A thermodynamic model of protein structure evolution explains empirical amino acid substitution matrices

Proteins evolve under a myriad of biophysical selection pressures that collectively control the patterns of amino acid substitutions. These evolutionary pressures are sufficiently consistent over time and across protein families to produce substitution patterns, summarized in global amino acid substitution matrices such as BLOSUM, JTT, WAG, and LG, which can be used to successfully detect homologs, infer phylogenies, and reconstruct ancestral sequences. Although the factors that govern the variation of amino acid substitution rates have received much attention, the influence of thermodynamic stability constraints remains unresolved. Here we develop a simple model to calculate amino acid substitution matrices from evolutionary dynamics controlled by a fitness function that reports on the thermodynamic effects of amino acid mutations in protein structures. This hybrid biophysical and evolutionary model accounts for nucleotide transition/transversion rate bias, multi‐nucleotide codon changes, the number of codons per amino acid, and thermodynamic protein stability. We find that our theoretical model accurately recapitulates the complex yet universal pattern observed in common global amino acid substitution matrices used in phylogenetics. These results suggest that selection for thermodynamically stable proteins, coupled with nucleotide mutation bias filtered by the structure of the genetic code, is the primary driver behind the global amino acid substitution patterns observed in proteins throughout the tree of life.     

一种可促进重组蛋白表达量和稳定性的多功能纯化标签的开发与利用

自组装双亲短肽(Self-assembling amphipathic peptides,SAPs)是一类亲疏水氨基酸按一定规律分布、具有自聚合效应的短肽,融合在酶蛋白N端时,具有促进表达和稳定化的功能。根据前期研究结论,设计一条全新的基于SAPs (S1vw,HNANARARHNANARARHNANARARHNARARAR)的可促进融合蛋白表达和稳定性,并可用于镍柱亲和层析的多功能短肽标签,在大肠杆菌表达系统中,将S1vw以PT-linker(PTPPTTPTPPTTPTP)融合在碱性果胶酶(Alkaline polygalacturonate lyase,PGL)、脂肪氧合酶(Lipoxygenase,LOX)及绿色荧光蛋白(Green fluorescent protein,GFP)的N末端时,与对应的野生型相比,PGL及LOX的粗酶活分别提高了3.1倍和1.89倍,GFP的荧光强度提高了16.22倍。S1vw的3种融合酶均可用镍柱进行亲和纯化,并具有较高的回收率。PGL及LOX在对应的热处理条件下,与野生型相比,半衰期分别提高了2.16倍和3.2倍。将GFP-S1vw在枯草芽孢杆菌及毕赤酵母表达系统中表达,发现在枯草芽孢杆菌中融合蛋白表达量提高明显,但在毕赤酵母中表达量几乎没有改变。说明在原核表达体系中,S1vw可作为一种新型的促表达、稳定化及可纯化的多功能标签。     

Cloning, overexpression, purification, and physicochemical characterization of a cold shock protein homolog from the hyperthermophilic bacterium Thermotoga maritima

Thermotoga maritima (Tm) expresses a 7 kDa monomeric protein whose 18 N-terminal amino acids show 81% identity to N-terminal sequences of cold shock proteins (Csps) from Bacillus caldolyticus and Bacillus stearothermophilus. There were only trace amounts of the protein in Thermotoga cells grown at 80 degrees C. Therefore, to perform physicochemical experiments, the gene was cloned in Escherichia coli. A DNA probe was produced by PCR from genomic Tm DNA with degenerated primers developed from the known N-terminus of TmCsp and the known C-terminus of CspB from Bacillus subtilis. Southern blot analysis of genomic Tm DNA allowed to produce a partial gene library, which was used as a template for PCRs with gene- and vector-specific primers to identify the complete DNA sequence. As reported for other csp genes, the 5' untranslated region of the mRNA was anomalously long; it contained the putative Shine-Dalgarno sequence. The coding part of the gene contained 198 bp, i.e., 66 amino acids. The sequence showed 61% identity to CspB from B. caldolyticus and high similarity to all other known Csps. Computer-based homology modeling allowed the conclusion that TmCsp represents a beta-barrel similar to CspB from B. subtilis and CspA from E. coli. As indicated by spectroscopic analysis, analytical gel permeation chromatography, and mass spectrometry, overexpression of the recombinant protein yielded authentic TmCsp with a molecular weight of 7,474 Da. This was in agreement with the results of analytical ultracentrifugation confirming the monomeric state of the protein. The temperature-induced equilibrium transition at 87 degrees C exceeds the maximum growth temperature of Tm and represents the maximal Tm-value reported for Csps so far.     

Directed evolution of a thermostable l-aminoacylase biocatalyst

Enzymes from extreme environments possess highly desirable traits of activity and stability for application under process conditions. One such example is l-aminoacylase (E.C. 3.5.1.14) from Thermococcus litoralis (TliACY), which catalyzes the enantioselective amide hydrolysis of N-protected l-amino acids, useful for resolving racemic mixtures in the preparation of chiral intermediates. Variants of this enzyme with improved activity and altered substrate preference are highly desirable. We have created a structural homology model of the enzyme and applied various two different directed evolution strategies to identify improved variants. Mutants P237S and F251Y were 2.4-fold more active towards N-benzoyl valine relative to the wild type at 65 °C. F251 mutations to basic residues resulted in 4.5-11-fold shifts in the substrate preference towards N-benzoyl phenylalanine relative to N-benzoyl valine. The substrate preference of wild type decreases with increasingly branched and sterically hindered substrates. However, the mutant S100T/M106K disrupted this simple trend by selectively improving the substrate preference for N-benzoyl valine, with a >30-fold shift in the ratio of k_cat values for N-benzoyl valine and N-benzoyl phenylalanine. Mutations that favoured N-benzoyl-phenylalanine appeared at the active site entrance, whereas those improving activity towards N-benzoyl-valine occurred in the hinge region loops linking the dimerization and zinc-binding domains in each monomer. These observations support a previously proposed substrate induced conformational transition between open and closed forms of aminoacylases.     

Rapid cooperative two-state folding of a miniature alpha-beta protein and design of a thermostable variant

There is a great deal of interest in developing small stably folded miniature proteins. A limited number of these molecules have been described, however they typically have not been characterized in depth. In particular, almost no detailed studies of the thermodynamics and folding kinetics of these proteins have been reported. Here we describe detailed studies of the thermodynamics and kinetics of folding of a 39 residue mixed alpha-beta protein (NTL9(1-39)) derived from the N-terminal domain of the ribosomal protein L9. The protein folds cooperatively and rapidly in a two-state fashion to a native state typical of those found for normal globular proteins. At pH 5.4 in 20mM sodium acetate, 100mM NaCl the temperature of maximum stability is 6 degrees C, the t(m) is 65.3 degrees C, deltaH degrees (t(m)) is between 24.6 kcalmol(-1) and 26.3 kcalmol(-1), and deltaC(p) degrees is 0.38 kcalmol(-1)deg(-1). The thermodynamic parameters are in the range expected on the basis of per residue values determined from databases of globular proteins. H/2H exchange measurements reveal a set of amides that exchange via global unfolding, exactly as expected for a normal cooperatively folded globular protein. Kinetic measurements show that folding is two-state folding. The folding rate is 640 s(-1) and the value of deltaG degrees calculated from the folding and unfolding rates is in excellent agreement with the equilibrium value. A designed thermostable variant, generated by mutating K12 to M, was characterized and found to have a t(m) of 82 degrees C. Equilibrium and kinetic measurements demonstrate that its folding is cooperative and two-state.     

Dissecting protein structure and function using directed evolution

To characterize the contributions of individual amino acids to the structure or function of a protein, researchers have adopted directed evolution approaches, which use iterated cycles of mutagenesis and selection or screening to search vast areas of sequence space for sets of mutations that provide insights into the protein of interest.     

Simulating protein evolution in sequence and structure space

Naturally occurring proteins comprise a special subset of all plausible sequences and structures selected through evolution. Simulating protein evolution with simplified and all-atom models has shed light on the evolutionary dynamics of protein populations, the nature of evolved sequences and structures, and the extent to which today's proteins are shaped by selection pressures on folding, structure and function. Extensive mapping of the native structure, stability and folding rate in sequence space using lattice proteins has revealed organizational principles of the sequence/structure map important for evolutionary dynamics. Evolutionary simulations with lattice proteins have highlighted the importance of fitness landscapes, evolutionary mechanisms, population dynamics and sequence space entropy in shaping the generic properties of proteins. Finally, evolutionary-like simulations with all-atom models, in particular computational protein design, have helped identify the dominant selection pressures on naturally occurring protein sequences and structures.     

The effects of temperature, pH and sulphide on the community structure of hyperthermophilic streamers in hot springs of northern Thailand

Hyperthermophilic community diversity was assessed in hot-spring streamers along gradients of temperature, pH and sulphide in northern Thailand. A hierarchical sampling design was employed to obtain biomass for culture-independent estimates of 16S rRNA gene-defined prokaryotic diversity. All springs supported several archaeal and bacterial phylotypes, including novel phylotypes that expand the known phylogenetic diversity of terrestrial hyperthermophiles. Diversity appeared significantly greater than that observed for several other geographic locations. Phylotypes belonging to the Aquificales were ubiquitous, further supporting the hypothesis that these chemolithoautotrophs are key members of all hyperthermophilic communities. The chemoorganotrophic genus Thermus was also represented by phylotypes in all springs. Other bacterial taxa represented by environmental sequences included Bacillus, Thermotoga and various unidentified Alphaproteobacteria and Betaproteobacteria. Archaeal phylotypes included the Crenarchaea Desulfurococcus, Pyrobaculum, plus several unidentified hyperthermophilic lineages. A Methanothermococcus-like Euryarchaeon was also identified, with this genus not previously known from streamer communities. A multivariate approach to the analysis of biotic and abiotic data revealed that diversity patterns were best explained by a combination of temperature and sulphide rather than by any other abiotic variable either individually or in combination.     

Unusual features of cereal seed protein structure and evolution

The alcohol-soluble (prolamin) storage proteins of barley, wheat and rye vary in their structures, but all have two features in common: the presence of distinct structural domains differing in amino acid compositions, and of repeats within one of these domains. Detailed comparisons of amino acid sequences show that all appear to have evolved from a single ancestral gene consisting of three short related regions (called A, B and C). Regions related to A, B and C are also present in the minor prolamins of maize and in three other groups of seed proteins: inhibitors of alpha-amylase and/or trypsin from cereals. 25 storage globulins from several dicotyledonous species and a 2S albumin from sunflower. It is suggested that these proteins together constitute a protein superfamily with limited sequence homology.     

The interface of protein structure, protein biophysics, and molecular evolution

Abstract The interface of protein structural biology, protein biophysics, molecular evolution, and molecular population genetics forms the foundations for a mechanistic understanding of many aspects of protein biochemistry. Current efforts in interdisciplinary protein modeling are in their infancy and the state-of-the art of such models is described. Beyond the relationship between amino acid substitution and static protein structure, protein function, and corresponding organismal fitness, other considerations are also discussed. More complex mutational processes such as insertion and deletion and domain rearrangements and even circular permutations should be evaluated. The role of intrinsically disordered proteins is still controversial, but may be increasingly important to consider. Protein geometry and protein dynamics as a deviation from static considerations of protein structure are also important. Protein expression level is known to be a major determinant of evolutionary rate and several considerations including selection at the mRNA level and the role of interaction specificity are discussed. Lastly, the relationship between modeling and needed high-throughput experimental data as well as experimental examination of protein evolution using ancestral sequence resurrection and in vitro biochemistry are presented, towards an aim of ultimately generating better models for biological inference and prediction.