A rigid network of long-range contacts increases thermostability in a mutant endoglucanase

Thermodynamic stability of a protein at elevated temperatures is a key factor for thermostable enzymes to catalyze their specific reactions. Yet our understanding of biological determinants of thermostability is far from complete. Many different atomistic factors have been suggested as possible means for such proteins to preserve their activity at high temperatures. Among these factors are specific local interatomic interactions or enrichment of specific amino acid types. The case of glycosyl hydrolase family endoglucanase of Trichoderma reesei defies current hypotheses for thermostability because a single mutation far from the active site (A35?V) converts this mesostable protein into a thermostable protein without significant change in the protein structure. This substantial change in enzymatic activity cannot be explained on the basis of local intramolecular interactions alone. Here we present a more global view of the induced thermostability and show that the A35?V mutation affects the underlying structural rigidity of the whole protein via a number of long-range, non-local interactions. Our analysis of this structure reveals a precisely tuned, rigid network of atomic interactions. This cooperative, allosteric effect promotes the transformation of this mesostable protein into a thermostable one.     

Environment-dependent long-range structural distortion in a temperature-sensitive point mutant

Extensive environment-dependent rearrangement of the helix-turn-helix DNA recognition region and adjacent L-tryptophan binding pocket is reported in the crystal structure of dimeric E. coli trp aporepressor with point mutation Leu75Phe. In one of two subunits, the eight residues immediately C-terminal to the mutation are shifted forward in helical register by three positions, and the five following residues form an extrahelical loop accommodating the register shift. In contrast, the second subunit has wildtype-like conformation, as do both subunits in an isomorphous wildtype control structure. Treated together as an ensemble pair, the distorted and wildtype-like conformations of the mutant apoprotein agree more fully than either conformation alone with previously reported NOE measurements, and account more completely for its diverse biochemical and biophysical properties. The register-shifted segment Ile79-Ala80-Thr81-Ile82-Thr83 is helical in both conformations despite low helical propensity, suggesting an important structural role for the steric constraints imposed by β-branched residues in helical conformation.     

Dimerisation and an increase in active site aromatic groups as adaptations to high temperatures: X-ray solution scattering and substrate-bound crystal structures of Rhodothermus marinus endoglucanase Cel12A

Cellulose, a polysaccharide consisting of beta-1,4-linked glucose, is the major component of plant cell walls and consequently one of the most abundant biopolymers on earth. Carbohydrate polymers such as cellulose are molecules with vast diversity in structure and function, and a multiplicity of hydrolases operating in concert are required for depolymerisation. The bacterium Rhodothermus marinus, isolated from shallow water marine hot springs, produces a number of carbohydrate-degrading enzymes including a family 12 cellulase Cel12A. The structure of R.marinus Cel12A in the ligand-free form (at 1.54 angstroms) and structures of RmCel12A after crystals were soaked in cellopentaose for two different lengths of time, have been determined. The shorter soaked complex revealed the conformation of unhydrolysed cellotetraose, while cellopentaose had been degraded more completely during the longer soak. Comparison of these structures with those of mesophilic family 12 cellulases in complex with inhibitors and substrate revealed that RmCel12A has a more extensive aromatic network in the active site cleft which ejects products after hydrolysis. The substrate structure confirms that during hydrolysis by family 12 cellulases glucose does not pass through a (2,5)B conformation. Small-angle X-ray scattering analysis of RmCel12A showed that the enzyme forms a loosely associated antiparallel dimer in solution, which may target the enzyme to the antiparallel polymer strands in cellulose.     

The structure of Rhodothermus marinus Cel12A,a highly thermostable family 12 endoglucanase,at 1.8 A resolution

Cellulose is one of the most abundant polysaccharides in nature and microorganisms have developed a comprehensive system for enzymatic breakdown of this ubiquitous carbon source, a subject of much interest in the biotechnology industry. Rhodothermus marinus produces a hyperthermostable cellulase, with a temperature optimum of more than 90 degrees C, the structure of which is presented here to 1.8 A resolution. The enzyme has been classified into glycoside hydrolase family 12; this is the first structure of a thermophilic member of this family to have been solved. The beta-jelly roll fold observed has identical topology to those of the two mesophilic members of the family whose structures have been elucidated previously. A Hepes buffer molecule bound in the active site may have triggered a conformational change to an active configuration as the two catalytic residues Glu124 and Glu207, together with dependent residues, are observed in a conformation similar to that seen in the structure of Streptomyces lividans CelB2 complexed with an inhibitor. The structural similarity between this cellulase and the mesophilic enzymes serves to highlight features that may be responsible for its thermostability, chiefly an increase in ion pair number and the considerable stabilisation of a mobile region seen in S. lividans CelB2. Additional aromatic residues in the active site region may also contribute to the difference in thermophilicity.     

Removal of distal protein-water hydrogen bonds in a plant epoxide hydrolase increases catalytic turnover but decreases thermostability

A putative proton wire in potato soluble epoxide hydrolase 1, StEH1, was identified and investigated by means of site-directed mutagenesis, steady-state kinetic measurements, temperature inactivation studies, and X-ray crystallography. The chain of hydrogen bonds includes five water molecules coordinated through backbone carbonyl oxygens of Pro(186), Leu(266), His(269), and the His(153) imidazole. The hydroxyl of Tyr(149) is also an integrated component of the chain, which leads to the hydroxyl of Tyr(154). Available data suggest that Tyr(154) functions as a final proton donor to the anionic alkylenzyme intermediate formed during catalysis. To investigate the role of the putative proton wire, mutants Y149F, H153F, and Y149F/H153F were constructed and purified. The structure of the Y149F mutant was solved by molecular replacement and refined to 2.0 A resolution. Comparison with the structure of wild-type StEH1 revealed only subtle structural differences. The hydroxyl group lost as a result of the mutation was replaced by a water molecule, thus maintaining a functioning hydrogen bond network in the proton wire. All mutants showed decreased catalytic efficiencies with the R,R-enantiomer of trans-stilbene oxide, whereas with the S,S-enantiomer, k (cat)/K (M) was similar or slightly increased compared with the wild-type reactions. k (cat) for the Y149F mutant with either TSO enantiomer was increased; thus the lowered enzyme efficiencies were due to increases in K (M). Thermal inactivation studies revealed that the mutated enzymes were more sensitive to elevated temperatures than the wild-type enzyme. Hence, structural alterations affecting the hydrogen bond chain caused increases in k (cat) but lowered thermostability.     

Importance of hydrophobic cluster formation through long-range contacts in the folding transition state of two-state proteins

Understanding the folding pathways of proteins is a challenging task. The Phi value approach provides a detailed understanding of transition-state structures of folded proteins. In this work, we have computed the hydrophobicity associated with each residue in the folded state of 16 two-state proteins and compared the Phi values of each mutant residue. We found that most of the residues with high Phi value coincide with local maximum in surrounding hydrophobicity, or have nearby residues that show such maximum in hydrophobicity, indicating the importance of hydrophobic interactions in the transition state. We have tested our approach to different structural classes of proteins, such as alpha-helical, SH3 domains of all-beta proteins, beta-sandwich, and alpha/beta proteins, and we observed a good agreement with experimental results. Further, we have proposed a hydrophobic contact network pattern to relate the Phi values with long-range contacts, which will be helpful to understand the transition-state structures of folded proteins. The present approach could be used to identify potential hydrophobic clusters that may form through long-range contacts during the transition state.     

Role of glycine 256 residue in improving the catalytic efficiency of mutant endoglucanase of family 5 glycoside hydrolase from Bacillus amyloliquefaciens SS35

Wild-type, BaGH5-WT and mutant, BaGH5-UV2 (aspartate residue mutated to glycine), endoglucanases belonging to glycoside hydrolase family 5 (GH5), from wild-type, and UV2 mutant strain of Bacillus amyloliquefaciens SS35, respectively, were earlier cloned in pHTP0 cloning vector. In this study, genes encoding BaGH5-WT or BaGH5-UV2 were cloned into pET28a(+) expression-vector and expressed in Escherichia coli BL-21(DE3)pLysS cells. BaGH5-UV2 showed 10-fold (43.6 U/mg) higher specific activity against carboxymethylcellulose sodium salt (CMC-Na), higher optimal temperature by 10°C at 65°C, and 22-fold higher catalytic efficiency against CMC-Na, than BaGH5-WT. BaGH5-UV2 showed stability in wider acidic pH range (5.0–7.0) unlike BaGH5-WT in narrow basic pH range (7.0–7.5). BaGH5-UV2 displayed a mutation, Asp256Gly in L11 loop, connecting β₆-sheet with α₆-helix, near active site toward the domain surface of (α/β)₈-TIM barrel fold. Molecular dynamics simulation studies showed more stable structure, accessibility of substrate for a catalytic site, and increased flexibility of loop L11 of BaGH5-UV2 than the wild type, suggesting enhanced catalysis by BaGH5-UV2. Molecular docking analysis displayed enhanced hydrogen bond interactions of cello-oligosaccharides with BaGH5-UV2, unlike BaGH5-WT. Thus, Gly256 residue of loop L11 plays an important role in enhancing catalytic efficiency, and pH stability of GH5 endoglucanase. Therefore, these results help in protein engineering of GH5 endoglucanase for improved biochemical properties.     

Ethylene glycol and the thermostability of trypsin in a reverse micelle system

The influence of ethylene glycol (EG) on the kinetics of hydrolysis of N-alpha-benzoyl-L-arginine ethyl ether catalyzed by trypsin encapsulated in sodium bis-(2-ethylhexyl)sulfosuccinate (AOT)-based reverse micelles was studied at different temperatures. Ethylene glycol was shown to shift the range of the trypsin activity in the reverse micelles towards higher temperatures. Infrared spectroscopy showed a stabilizing effect of EG on the secondary structure of the protein in the system of reverse micelles. Electron spin resonance spectroscopy showed that the solubilized protein affected the interactions of EG with the polar head groups of AOT and altered the rigidity of the micellar matrix. The results indicate that EG increases the thermostability of the solubilized enzyme in microemulsion media by two mechanisms.     

木糖异构酶序列结构特点、耐热机理及分子改造研究进展

近年来,随着发展低碳经济的迫切需要,可再生资源利用研究方兴未艾,其中,构建充分利用木质纤维素水解产物木糖生产乙醇的重组菌成为研究热点.木糖异构酶由于不需要辅酶,成为构建利用木糖重组酵母的首选途径.文中对近年来木糖异构酶的研究进展进行了综述.首先介绍了木糖异构酶的基本性质、序列、结构和功能特性,然后对其耐热机理进行了总结归纳;重点阐述了基于序列及结构进行的酶分子改造研究,包括底物特异性改造、热稳定性改造等;同时,结合作者的研究经历,对如何提高嗜热木糖异构酶在常温下的活性进行了探讨.最后,对木糖异构酶的研究进展进行了总结和展望,对基于结构改善木糖异构酶催化活性及构建新型高效利用木糖生产乙醇的重组菌具有重要的指导意义.     

Potential hydrophobic interaction between two cysteines in interior hydrophobic region improves thermostability of a family 11 xylanase from Neocallimastix Patriciarum

In this study, we employed directed evolution and site‐directed mutagenesis to screen thermostable mutants of a family 11 xylanase from Neocallimastix patriciarum, and found that the thermostability and specific activity are both enhanced when mutations (G201C and C60A) take place in the interior hydrophobic region of the enzyme. Far‐ultraviolet circular dichroism analysis showed that the melting temperatures (T_m) of the G201C and C60A–G201C mutants are higher than that of the wild type by about 10 and 12°C, respectively. At 72°C, their specific activities are about 4 and 6 times as that of the wild type, respectively. Homology modeling and site‐directed mutagenesis demonstrated that the enhanced thermostability of the G201C and C60A–G201C mutants may be mainly attributed to a potential stronger hydrophobic interaction between the two well‐packed cysteines at sites 50 and 201, rather than the disulfide bond formation which was ruled out by thiol titration with dithionitrobenzoic acid (DTNB). And the strength of such interaction depends on the packing of the side‐chain and hydrophobicity of residues at these two sites. This suggests that cysteine could stabilize a protein not only by forming a disulfide bond, but also by the strong hydrophobicity itself. Biotechnol. Bioeng. 2010;105: 861–870. © 2009 Wiley Periodicals, Inc.     

Understanding the alpha-helix to coil transition in polypeptides using network rigidity: predicting heat and cold denaturation in mixed solvent conditions

Thermodynamic stability in polypeptides is described using a novel Distance Constraint Model (DCM). Here, microscopic interactions are represented as constraints. A topological arrangement of constraints define a mechanical framework. Each constraint in the framework is associated with an enthalpic and entropic contribution. All accessible topological arrangements of distance constraints form an ensemble of mechanical frameworks, each representing a microstate of the polypeptide. A partition function is calculated exactly using a transfer matrix approach, where in many respects the DCM is similar to the Lifson-Roig model. The crucial difference is that the effect of network rigidity is explicitly calculated for each mechanical framework in the ensemble. Network rigidity is a mechanical interaction that provides a mechanism for long-range molecular cooperativity and enables a proper treatment of the nonadditivity of a microscopic free energy decomposition. Accounting for (1) helix <--> coil conformation changes along the backbone similar to the Lifson-Roig model, (2) i to i + 4 hydrogen-bond formation <--> breaking similar to the Zimm-Bragg model, and (3) structured <--> unstructured solvent interaction (hydration effects), a six-parameter DCM describes normal and inverted helix-coil transitions in polypeptides. Under suitable mixed solvent conditions heat and cold denaturation is predicted. Model parameters are fitted to experimental data showing different degrees of cold denaturation in monomeric polypeptides in aqueous hexafluoroisopropanol (HFIP) solution at various HFIP concentrations. By assuming a linear HFIP concentration dependence (up to 6% by mole fraction) on model parameters, all essential experimentally observed features are captured.     

An enhanced elastic network model to represent the motions of domain-swapped proteins

Domain swapping is a process where two (or more) protein molecules form a dimer (or higher oligomer) by exchanging an identical domain. In this article, based on the observation that domains are rigid and hinge loops are highly flexible, we propose a new Elastic Network Model, domain-ENM, for domain-swapped proteins. In this model, the rigidity of domains is taken into account by using a larger spring constant for intradomain contacts. The large-scale transition of domain swapping is then novelly decomposed into the relative motion between the rigid domains (only 6 degrees of freedom) plus the internal fluctuations of each domain. Consequently, this approach has the potential to produce much more meaningful transition pathways than other simulation approaches that try to find pathways in a search space of large numbers of dimensions. In this article, we also propose a new way to define the overlap measure. Past approaches used an inappropriate comparison of the large-scale conformation displacement against the computed infinitesimal motions of modes. Here, we propose an infinitesimal version of the large-scale conformation change and then compare it with the modes of motions. As a result, we obtain much better overlap values. Using this new overlap definition, we are also able for the first time to give a clear, intuitive explanation why "open" forms tend to produce better overlap values than "closed" forms with traditional ENMs. Finally, as an application, we present a simple approach to show how domain-ENM can be used to generated transition pathways for domain-swapped proteins.     

The Limited Role of Nonnative Contacts in the Folding Pathways of a Lattice Protein

Models of protein energetics that neglect interactions between amino acids that are not adjacent in the native state, such as the Gō model, encode or underlie many influential ideas on protein folding. Implicit in this simplification is a crucial assumption that has never been critically evaluated in a broad context: Detailed mechanisms of protein folding are not biased by nonnative contacts, typically argued to be a consequence of sequence design and/or topology. Here we present, using computer simulations of a well-studied lattice heteropolymer model, the first systematic test of this oft-assumed correspondence over the statistically significant range of hundreds of thousands of amino acid sequences that fold to the same native structure. Contrary to previous conjectures, we find a multiplicity of folding mechanisms, suggesting that Gō-like models cannot be justified by considerations of topology alone. Instead, we find that the crucial factor in discriminating among topological pathways is the heterogeneity of native contact energies: The order in which native contacts accumulate is profoundly insensitive to omission of nonnative interactions, provided that native contact heterogeneity is retained. This robustness holds over a surprisingly wide range of folding rates for our designed sequences. Mirroring predictions based on the principle of minimum frustration, fast-folding sequences match their Gō-like counterparts in both topological mechanism and transit times. Less optimized sequences dwell much longer in the unfolded state and/or off-pathway intermediates than do Gō-like models. For dynamics that bridge unfolded and unfolded states, however, even slow folders exhibit topological mechanisms and transit times nearly identical with those of their Gō-like counterparts. Our results do not imply a direct correspondence between folding trajectories of Gō-like models and those of real proteins, but they do help to clarify key topological and energetic assumptions that are commonly used to justify such caricatures.     

Increase in the thermostability of Bacillus sp. strain TAR-1 xylanase using a site saturation mutagenesis library

Site saturation mutagenesis library is a recently developed technique, in which any one out of all amino acid residues in a target region is substituted into other 19 amino acid residues. In this study, we used this technique to increase the thermostability of a GH10 xylanase, XynR, from Bacillus sp. strain TAR-1. We hypothesized that the substrate binding region of XynR is flexible, and that the thermostability of XynR will increase if the flexibility of the substrate binding region is decreased without impairing the substrate binding ability. Site saturation mutagenesis libraries of amino acid residues Tyr43–Lys115 and Ala300–Asn325 of XynR were constructed. By screening 480 clones, S92E was selected as the most thermostable one, exhibiting the residual activity of 80% after heat treatment at 80°C for 15 min in the hydrolysis of Remazol Brilliant Blue-xylan. Our results suggest that this strategy is effective for stabilization of GH10 xylanase.

Abbreviations: DNS: 3,5-dinitrosalicylic acid; RBB-xylan: Remazol Brilliant Blue-xylan     

Primary recovery of a genetically engineered Trichoderma reesei endoglucanase I (Cel 7B) fusion protein in cloud point extraction systems

Here we present data to demonstrate how partitioning of a hydrophilic enzyme can be directed to the hydrophobic detergent-enriched phase of an aqueous two-phase system by addition of short stretches of amino acid residues to the protein molecule. The target enzyme was the industrially important endoglucanase I, EGI (endo-1,4-beta-D-glucan-4-glucanohydrolase, EC 3.2.1.4, Cel7B) of Trichoderma reesei. We investigated the partitioning of three EGI variants containing various C-terminal peptide extensions including Trp-Pro motifs of different lengths and localizations. Additionally, a recently developed system composed of the thermoseparating copolymer HM-EOPO was utilized to study the effects of fusion tags. The addition of peptides containing tryptohan residues enhanced the partitioning of EGI to the HM-EOPO-rich phase. The system composed of a nonionic detergent (Agrimul NRE1205) resulted in the highest partition coefficient (K = 31) and yield (90%) with the construct EGI(core-P5)(WP)(4) containing (Trp-Pro)(4) after a short linker stretch. A recombinant strain of T. reesei Rut-C30 for large-scale production was constructed in which the fusion protein EGI(core-P5)(WP)(4) was expressed from the strong promoter of the cellulase gene cbh1. The fusion protein was successfully expressed and secreted from the fungus during shake-flask cultivations. Cultivation in a 28-L bioreactor however, revealed that the fusion protein is sensitive to proteases. Consequently, only low production levels were obtained in large-scale production trials.     

Crystal structure of heat-labile enterotoxin from Escherichia coli with increased thermostability introduced by an engineered disulfide bond in the A subunit.

Cholera toxin (CT) produced by Vibrio cholerae and heat-labile enterotoxin (LT-I), produced by enterotoxigenic Escherichia coli, are AB5 heterohexamers with an ADP-ribosylating A subunit and a GM1 receptor binding B pentamer. These toxins are among the most potent mucosal adjuvants known and, hence, are of interest both for the development of anti-diarrheal vaccines against cholera or enterotoxigenic Escherichia coli diarrhea and also for vaccines in general. However, the A subunits of CT and LT-I are known to be relatively temperature sensitive. To improve the thermostability of LT-I an additional disulfide bond was introduced in the A1 subunit by means of the double mutation N40C and G166C. The crystal structure of this double mutant of LT-I has been determined to 2.0 A resolution. The protein structure of the N40C/G166C double mutant is very similar to the native structure except for a few local shifts near the new disulfide bond. The introduction of this additional disulfide bond increases the thermal stability of the A subunit of LT-I by 6 degrees C. The enhancement in thermostability could make this disulfide bond variant of LT-I of considerable interest for the design of enterotoxin-based vaccines.     

Staying in touch with the endocytic network: The importance of contacts for cholesterol transport

Cholesterol homeostasis is critical for cell function and human health. Cholesterol is heterogeneously distributed among cellular membranes, with the redistribution of endocytosed dietary cholesterol playing a pivotal role in the regulation of cholesterol homeostasis. While gaps remain in our understanding of intracellular dietary cholesterol transport, a highly complex network of pathways is starting to emerge, often involving inter‐dependent vesicular and non‐vesicular transport mechanisms. The last decade has seen a surge in interest in non‐vesicular transport and inter‐organellar communication at membrane contact sites. By providing platforms for protein interactions, signalling events, lipid exchange and calcium flux, membrane contact sites (MCS) are now appreciated as controlling the fate of large amounts of lipid and play central roles in the regulation and co‐ordination of endocytic trafficking. Here, we review the role of MCS in multiple pathways for cholesterol export from the endocytic pathway and highlight the intriguing interplay between vesicular and non‐vesicular transport mechanisms and relationship with neurodegenerative disease.     

Random main effects of treatment: A case study with a network meta‐analysis

If the number of treatments in a network meta‐analysis is large, it may be possible and useful to model the main effect of treatment as random, that is to say as random realizations from a normal distribution of possible treatment effects. This then constitutes a third sort of random effect that may be considered in connection with such analyses. The first and most common models treatment‐by‐trial interaction as being random and the second, rather rarer, models the main effects of trial as being random and thus permits the recovery of intertrial information. Taking the example of a network meta‐analysis of 44 similar treatments in 10 trials, we illustrate how a hierarchical approach to modeling a random main effect of treatment can be used to produce shrunk (toward the overall mean) estimates of effects for individual treatments. As a related problem, we also consider the issue of using a random‐effect model for the within‐trial variances from trial to trial. We provide a number of possible graphical representations of the results and discuss the advantages and disadvantages of such an approach.     

The evolutionary dynamics of evolvability in a gene network model

Evolvability, the ability of populations to adapt, has recently emerged as a major unifying concept in biology. Although the study of evolvability offers new insights into many important biological questions, the conceptual bases of evolvability, and the mechanisms of its evolution, remain controversial. We used simulated evolution of a model of gene network dynamics to test the contentious hypothesis that natural selection can favour high evolvability, in particular in sexual populations. Our results conclusively demonstrate that fluctuating natural selection can increase the capacity of model gene networks to adapt to new environments. Detailed studies of the evolutionary dynamics of these networks establish a broad range of validity for this result and quantify the evolutionary forces responsible for changes in evolvability. Analysis of the genotype–phenotype map of these networks also reveals mechanisms connecting evolvability, genetic architecture and robustness. Our results suggest that the evolution of evolvability can have a pervasive influence on many aspects of organisms.