Improving Trichoderma reesei Cel7B Thermostability by Targeting the Weak Spots

For proteins that denature irreversibly, the denaturation is typically triggered by a partial unfolding, followed by a permanent change (e.g., aggregation). The regions that initiate the partial unfolding are named "weak spots". In this work, a molecular dynamics (MD) simulation and data analysis protocol is developed to identify the weak spots of Trichoderma reesei Cel7B, an important endoglucanase in cellulose hydrolysis, through assigning the local melting temperature (Tmp) to individual residue pairs. To test the predicted weak spots, a total of eight disulfide bonds were designed in these regions and all enhanced the enzyme thermostability. The increased stability, quantified by ΔT50 (which is the T50 difference between the mutant and the wild type enzyme), is negatively correlated with the MD-predicted Tmp, demonstrating the effectiveness of the protocol and highlighting the importance of the weak spots. Strengthening interactions in these regions proves to be a useful strategy in improving the thermostability of Tr. Cel7B.

[1]  Valerie Daggett,et al.  Unifying features in protein-folding mechanisms , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[2]  V. Eijsink,et al.  Rational engineering of enzyme stability. , 2004, Journal of biotechnology.

[3]  Frits Goedegebuur,et al.  Hypocrea jecorina CEL6A protein engineering , 2010, Biotechnology for biofuels.

[4]  V. Daggett,et al.  Staphylococcal protein A: unfolding pathways, unfolded states, and differences between the B and E domains. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[5]  Gerrit Groenhof,et al.  GROMACS: Fast, flexible, and free , 2005, J. Comput. Chem..

[6]  Y. Duan,et al.  Probing the stability-limiting regions of an antibody single-chain variable fragment: a molecular dynamics simulation study. , 2011, Protein engineering, design & selection : PEDS.

[7]  Indira Wu,et al.  Engineered thermostable fungal Cel6A and Cel7A cellobiohydrolases hydrolyze cellulose efficiently at elevated temperatures , 2013, Biotechnology and bioengineering.

[8]  S. Grzesiek,et al.  Temperature-dependence of protein hydrogen bond properties as studied by high-resolution NMR. , 2002, Journal of molecular biology.

[9]  A. Fersht,et al.  Simulation and experiment at high temperatures: ultrafast folding of a thermophilic protein by nucleation-condensation. , 2005, Journal of molecular biology.

[10]  Dick B Janssen,et al.  Molecular dynamics simulations as a tool for improving protein stability. , 2002, Protein engineering.

[11]  L. Viikari,et al.  Cellulases without carbohydrate-binding modules in high consistency ethanol production process , 2014, Biotechnology for Biofuels.

[12]  Alan A. Dombkowski,et al.  Disulfide by DesignTM: a computational method for the rational design of disulfide bonds in proteins , 2003, Bioinform..

[13]  N. Pokala,et al.  Energy functions for protein design: adjustment with protein-protein complex affinities, models for the unfolded state, and negative design of solubility and specificity. , 2005, Journal of molecular biology.

[14]  L. Viikari,et al.  The role of carbohydrate binding module (CBM) at high substrate consistency: comparison of Trichoderma reesei and Thermoascus aurantiacus Cel7A (CBHI) and Cel5A (EGII). , 2013, Bioresource technology.

[15]  V. Daggett,et al.  Temperature dependence of the flexibility of thermophilic and mesophilic flavoenzymes of the nitroreductase fold. , 2010, Protein engineering, design & selection : PEDS.

[16]  Martin J. Baumann,et al.  Origin of Initial Burst in Activity for Trichoderma reesei endo-Glucanases Hydrolyzing Insoluble Cellulose* , 2011, The Journal of Biological Chemistry.

[17]  T. Poulos,et al.  Using molecular dynamics to probe the structural basis for enhanced stability in thermal stable cytochromes P450. , 2010, Biochemistry.

[18]  L. Yao,et al.  The protein amide 1HN chemical shift temperature coefficient reflects thermal expansion of the N–H···O=C hydrogen bond , 2012, Journal of Biomolecular NMR.

[19]  D. Baker,et al.  Role of conformational sampling in computing mutation‐induced changes in protein structure and stability , 2011, Proteins.

[20]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[21]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[22]  D. Baker,et al.  High-resolution Structural and Thermodynamic Analysis of Extreme Stabilization of Human Procarboxypeptidase by Computational Protein Design , 2007, Journal of molecular biology.

[23]  V. Eijsink,et al.  Directed evolution of enzyme stability. , 2005, Biomolecular engineering.

[24]  H. Berendsen,et al.  Transport Properties Computed by Linear Response through Weak Coupling to a Bath , 1991 .

[25]  Yu Zhang,et al.  A mechanistic study of Trichoderma reesei Cel7B catalyzed glycosidic bond cleavage. , 2013, The journal of physical chemistry. B.

[26]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[27]  Huan-Xiang Zhou,et al.  Similarity and difference in the unfolding of thermophilic and mesophilic cold shock proteins studied by molecular dynamics simulations. , 2006, Biophysical journal.

[28]  F. Arnold,et al.  Directed evolution of a thermostable esterase. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[29]  P. Kollman,et al.  Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models , 1992 .

[30]  G J Kleywegt,et al.  The crystal structure of the catalytic core domain of endoglucanase I from Trichoderma reesei at 3.6 A resolution, and a comparison with related enzymes. , 1997, Journal of molecular biology.

[31]  Chie Motono,et al.  Thermodynamic and kinetic determinants of Thermotoga maritima cold shock protein stability: A structural and dynamic analysis , 2008, Proteins.

[32]  Amanda L. Jonsson,et al.  The role of the turn in beta-hairpin formation during WW domain folding. , 2007, Protein Science.

[33]  David K. Johnson,et al.  Biomass Recalcitrance: Engineering Plants and Enzymes for Biofuels Production , 2007, Science.

[34]  Jens Meiler,et al.  ROSETTA3: an object-oriented software suite for the simulation and design of macromolecules. , 2011, Methods in enzymology.

[35]  V. Eijsink,et al.  Structural determinants of the stability of thermolysin-like proteinases , 1995, Nature Structural Biology.

[36]  A. Gronenborn,et al.  Insights into conformation and dynamics of protein GB1 during folding and unfolding by NMR. , 2004, Journal of molecular biology.

[37]  A. Bommarius,et al.  The cellulose-binding domain of cellobiohydrolase Cel7A from Trichoderma reesei is also a thermostabilizing domain. , 2011, Journal of biotechnology.

[38]  M. Lever.  Carbohydrate determination with 4-hydroxybenzoic acid hydrazide (PAHBAH): effect of bismuth on the reaction. , 1977, Analytical biochemistry.

[39]  Piero Fariselli,et al.  I-Mutant2.0: predicting stability changes upon mutation from the protein sequence or structure , 2005, Nucleic Acids Res..

[40]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[41]  Holger Gohlke,et al.  Thermostabilizing mutations preferentially occur at structural weak spots with a high mutation ratio. , 2012, Journal of biotechnology.

[42]  Frits Goedegebuur,et al.  Hypocrea jecorina CEL 6 A protein engineering , 2010 .

[43]  D. Bolam,et al.  Carbohydrate-binding modules: fine-tuning polysaccharide recognition. , 2004, The Biochemical journal.

[44]  G. Schreiber,et al.  Assessing computational methods for predicting protein stability upon mutation: good on average but not in the details. , 2009, Protein engineering, design & selection : PEDS.

[45]  M. Parrinello,et al.  Polymorphic transitions in single crystals: A new molecular dynamics method , 1981 .

[46]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[47]  Amanda L. Jonsson,et al.  The role of the turn in β‐hairpin formation during WW domain folding , 2007 .

[48]  Xiangfei Song,et al.  Overexpression of β-glucosidase from Thermotoga maritima for the production of highly purified aglycone isoflavones from soy flour , 2009 .

[49]  Martin A. Nowaka,et al.  ‘ Protein folding and unfolding in microseconds to nanoseconds by experiment and simulation , 2001 .

[50]  S. Grzesiek,et al.  Key stabilizing elements of protein structure identified through pressure and temperature perturbation of its hydrogen bond network. , 2012, Nature chemistry.

[51]  Holger Gohlke,et al.  Constraint Network Analysis (CNA): A Python Software Package for Efficiently Linking Biomacromolecular Structure, Flexibility, (Thermo-)Stability, and Function , 2013, J. Chem. Inf. Model..

[52]  T. Reinikainen,et al.  The three-dimensional crystal structure of the catalytic core of cellobiohydrolase I from Trichoderma reesei. , 1994, Science.

[53]  Gert Vriend,et al.  Extreme Stabilization of a Thermolysin-like Protease by an Engineered Disulfide Bond* , 1997, The Journal of Biological Chemistry.

[54]  G. L. Miller.  Use of Dinitrosalicylic Acid Reagent for Determination of Reducing Sugar , 1959 .

[55]  Berk Hess,et al.  LINCS: A linear constraint solver for molecular simulations , 1997, J. Comput. Chem..

[56]  B. L. de Groot,et al.  Predicting free energy changes using structural ensembles. , 2009, Nature methods.

[57]  Andreas Vogel,et al.  Iterative saturation mutagenesis on the basis of B factors as a strategy for increasing protein thermostability. , 2006, Angewandte Chemie.

[58]  L. Viikari,et al.  Carbohydrate-binding modules (CBMs) revisited: reduced amount of water counterbalances the need for CBMs , 2013, Biotechnology for Biofuels.

[59]  Hein J. Wijma,et al.  Computationally designed libraries for rapid enzyme stabilization , 2014, Protein engineering, design & selection : PEDS.

[60]  Wolfgang Gärtner,et al.  First evidence for phototropin-related blue-light receptors in prokaryotes. , 2002, Biophysical journal.

[61]  L. Serrano,et al.  Predicting changes in the stability of proteins and protein complexes: a study of more than 1000 mutations. , 2002, Journal of molecular biology.

[62]  Jan H. Jensen,et al.  PROPKA3: Consistent Treatment of Internal and Surface Residues in Empirical pKa Predictions. , 2011, Journal of chemical theory and computation.

[63]  G Vriend,et al.  Early Steps in the Unfolding of Thermolysin-like Proteases* , 1998, The Journal of Biological Chemistry.

[64]  Tong Li,et al.  Engineering a More Thermostable Blue Light Photo Receptor Bacillus subtilis YtvA LOV Domain by a Computer Aided Rational Design Method , 2013, PLoS Comput. Biol..

[65]  Tuula T. Teeri,et al.  The roles and function of cellulose-binding domains , 1997 .

[66]  J. Buchert,et al.  Trichoderma reesei cellulases and their core domains in the hydrolysis and modification of chemical pulp , 2000 .

[67]  S. Nosé,et al.  Constant pressure molecular dynamics for molecular systems , 1983 .

[68]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[69]  Holger Gohlke,et al.  Protein rigidity and thermophilic adaptation , 2011, Proteins.

[70]  H. van Tilbeurgh,et al.  Studies of the cellulolytic system of the filamentous fungus Trichoderma reesei QM 9414. Substrate specificity and transfer activity of endoglucanase I. , 1990, The Biochemical journal.