Understanding the impacts of self-shuffling approach on structure and function of shuffled endoglucanase enzyme via MD simulations

Aslı YenenlerORCID iD: https://orcid.org/0000-0002-9169-388X, Umut Gerlevik
  • Department of Biostatistics and Medical Informatics, School of Medicine, Acibadem Mehmet Ali Aydinlar University, Istanbul, Turkey
  • Department of Molecular Biology and Genetics, Istanbul Technical University, Istanbul, Turkey
  • Email
  • Search for other articles:
  • degruyter.comGoogle Scholar
and Ugur Sezerman
  • Department of Biostatistics and Medical Informatics, School of Medicine, Acibadem Mehmet Ali Aydinlar University, Istanbul, Turkey
  • Email
  • Search for other articles:
  • degruyter.comGoogle Scholar



We identify the impacts of structural differences on functionality of EG3_S2 endoglucanase enzyme with MD studies. The results of previous experimental studies have been explained in details with computational approach. The objective of this study is to explain the functional differences between shuffled enzyme (EG3_S2) and its native counterpart (EG3_nat) from Trichoderma reseei, via Molecular Dynamics approach.

Materials and methods

For this purpose, we performed MD simulations along 30 ns at three different reaction temperatures collected as NpT ensemble, and then monitored the backbone motion, flexibilities of residues, and intramolecular interactions of EG3_S2 and EG3_nat enzymes.


According to MD results, we conclude that EG3_S2 and EG3_nat enzymes have unique RMSD patterns, e.g. RMSD pattern of EG3_S2 is more dynamic than that of EG3_nat at all temperatures. In addition to this dynamicity, EG3_S2 establishes more salt bridge interactions than EG3_nat.


By taking these results into an account with the preservation of catalytic Glu residues in a proper manner, we explain the structural basis of differences between shuffled and native enzyme via molecular dynamic studies.

    • Supplementary material
  • 1.

    Kurtzman AL, Govindarajan S, Vahle K, Jones JT, Heinrichs V, Patten PA. Advances in directed protein evolution by recursive genetic recombination: applications to therapeutic proteins. Curr Opin Biotechnol 2001;12:361–70.

  • 2.

    Yenenler A, Sezerman OU. Design and characterizations of two novel cellulases through single-gene shuffling of Cel12A (EG3) gene from Trichoderma reseei. Protein Eng Des Sel 2016;29:219–29.

  • 3.

    Stemmer WP. Molecular breedings of genes, pathways and genomes by DNA shuffling. Biotechnol Bioprocess Eng 2002;7:121–9.

  • 4.

    Jackel C, Kast P, Hilvert D. Protein design by directed evolution. Annu Rev Biophys 2008;37:153–73.

  • 5.

    Song JK, Rhee S. Simultaneous enhancement of thermostability and catalytic activity of phospholipase A(1) by evolutionary molecular engineering. Appl Environ Microbiol 2000;66:890–4.

  • 6.

    Akbulut N, Ozturk MT, Pijning T, Ozturk SI, Gumusel F. Improved activity and thermostability of Bacillus pumilus lipase by directed evolution. J Biotechnol 2013;164:123–9.

  • 7.

    Crameri A, Whitehorn EA, Tate E, Stemmer WP. Improved green fluorescent protein by molecular evolution using DNA shuffling. Nat Biotechnol 1996;14:315–9.

  • 8.

    Yano T, Oue S, Kagamiyama H. Directed evolution of an aspartate aminotransferase with new substrate specificities. Proc Natl Acad Sci USA 1998;95:5511–5.

  • 9.

    Chockalingam K, Chen Z, Katzenellenbogen JA, Zhao H. Directed evolution of specific receptor-ligand pairs for use in the creation of gene switches. Proc Natl Acad Sci USA 2005;102:5691–6.

  • 10.

    Wilson DB. Cellulases and biofuels. Curr Opin Biotechnol 2009;20:295–9.

  • 11.

    Bajpai P. Application of enzymes in the pulp and paper industry. Biotechnol Prog 1999;15:147–57.

  • 12.

    Bayer EA, Morag E, Lamed R. The cellulosome – a treasure-trove for biotechnology. Trends Biotechnol 1994;12:379–86.

  • 13.

    Ulker A, Sprey B. Characterization of an unglycosylated low molecular weight 1,4-beta-glucan-glucanohydrolase of Trichoderma reesei. FEMS Microbiol Lett 1990;57:215–9.

  • 14.

    Sandgren M, Shaw A, Ropp TH, Wu S, Bott R, Cameron AD, et al. The X-ray crystal structure of the Trichoderma reesei family 12 endoglucanase 3, Cel12A, at 1.9 A resolution. J Mol Biol 2001;308:295–310.

  • 15.

    Nidetzky B, Steiner W, Hayn M, Claeyssens M. Cellulose hydrolysis by the cellulases from Trichoderma reesei: a new model for synergistic interaction. Biochem J 1994;298:705–10.

  • 16.

    Zhang Y. I-TASSER server for protein 3D structure prediction. BMC Bioinformatics 2008;9:40.

  • 17.

    Yang J, Zhang Y. I-TASSER server: new development for protein structure and function predictions. Nucleic Acids Res 2015;43:174–81.

  • 18.

    Humphrey W, Dalke A, Schulten K. VMD: visual molecular dynamics. J Mol Graph 1996;14:33–8, 27–8.

  • 19.

    Phillips JC, Braun R, Wang W, Gumbart J, Tajkhorshid E, Villa E, et al. Scalable molecular dynamics with NAMD. J Comput Chem 2005;26:1781–802.

  • 20.

    Brooks BR, Brooks CL, Mackerell AD, Nilsson L, Petrella RJ, Roux B, et al. CHARMM: the biomolecular simulation program. J Comput Chem 2009;30:1545–614.

  • 21.

    MacKerell AD, Bashford D, Bellott M, Dunbrack RL, Evanseck JD, Field MJ, et al. All-atom empirical potential for molecular modeling and dynamics studies of proteins. J Phys Chem B 1998;102:3586–616.

  • 22.

    MacKerell AD, Feig M, Brooks CL. Improved treatment of the protein backbone in empirical force fields. J Am Chem Soc 2004;126:698–9.

  • 23.

    Jorgensen WL, Chandrasekhar J, Madura JD. Comparison of simple potential functions for simulating liquid water. J Chem Phys 1983;79:926–35.

  • 24.

    Darden T, York D, Pedersen L. Particle mesh Ewald: An Nlog(N) method for Ewald sums in large systems. J Chem Phys B 1993;98:10089.

  • 25.

    Essman U, Perera L, Berkowitz M, Darden T, Lee H, Pedersen GG. A smooth particle mesh Ewald method. J Chem Phys 1995;19:8577–93.

  • 26.

    Luthy R, Bowie JU, Eisenberg D. Assessment of protein models with three-dimensional profiles. Nature 1992;356:83–5.

  • 27.

    Bowie JU, Lüthy R, Eisenberg D. A method to identify protein sequences that fold into a known three-dimensional structure. Science (New York, N.Y.) 1991;253:164–70.

  • 28.

    Colovos CT, Yeates O. Verification of protein structures: patterns of nonbonded atomic interactions. Protein Sci 1993;2:1511–9.

  • 29.

    Pontius J, Richelle J, Wodak SJ. Deviations from standard atomic volumes as a quality measure for protein crystal structures. J Mol Biol 1996;264:121–36.

  • 30.

    Laskowski RA, MacArthur MW, Moss DS, Thornton JM. PROCHECK: a program to check the stereochemical quality of protein structures. J Appl Crystallogr 1993;26:283–91.

  • 31.

    Hooft RW, Vriend G, Sander C, Abola EE. Errors in protein structures. Nature 1996;381:272.

  • 32.

    Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, et al. The Protein Data Bank. Nucleic Acids Res 2000;28:235–42.

  • 33.

    Vriend G. WHAT IF: a molecular modeling and drug design program. J Mol Graph 1990;8:52–6.

  • 34.

    Kleywegt GJ. Jones TA. Phi/psi-chology: Ramachandran revisited. Structure 1996;4:1395–400.

  • 35.

    Yenenler A, Venturini A, Burduroglu HC, Sezerman OU. Investigating the structural properties of the active conformation BTL2 of a lipase from Geobacillus thermocatenulatus in toluene using molecular dynamic simulations and engineering BTL2 via in-silico mutation. J Mol Model 2018;24:229.

  • 36.

    Chu X, Wang J. Specificity and affinity quantification of flexible recognition from underlying energy landscape topography. PLoS Comput Biol 2014;10:1003782.

  • 37.

    Spolar RS, Record MT. Coupling of local folding to site-specific binding of proteins to DNA. Science 1994;263:777–84.

  • 38.

    Grunberg R, Nilges M, Leckner J. Flexibility and conformational entropy in protein-protein binding. Structure 2006;14:683–93.

  • 39.

    Marti DN, Bosshard HR. Electrostatic interactions in leucine zippers: thermodynamic analysis of the contributions of Glu and His residues and the effect of mutating salt bridges. J Mol Biol 2003;330:621–37.

  • 40.

    Anderson DE, Becktel WJ, Dahlquist FW. pH-induced denaturation of proteins: a single salt bridge contributes 3–5 kcal/mol to the free energy of folding of T4 lysozyme. Biochemistry 1990;29:2403–8.

  • 41.

    Jaenicke R, Bohm G. The stability of proteins in extreme environments. Curr Opin Struct Biol 1998;8:738–48.

  • 42.

    Kirino H, Aoki M, Aoshima M, Hayashi Y, Ohba M, Yamagishi A, et al. Hydrophobic interaction at the subunit interface contributes to the thermostability of 3-isopropylmalate dehydrogenase from an extreme thermophile, Thermus thermophilus. Eur J Biochem 1994;220:275–81.

  • 43.

    Tanaka S, Igarashi S, Ferri S, Sode K. Increasing stability of water-soluble PQQ glucose dehydrogenase by increasing hydrophobic interaction at dimeric interface. BMC Biochem 2005;6:1.

Purchase article
Get instant unlimited access to the article.
Price including VAT
Log in
Already have access? Please log in.

Journal + Issues

Turkish Journal of Biochemistry (TJB), official journal of Turkish Biochemical Society, is issued electronically every 2 months. The main aim of the journal is to support the research and publishing culture by ensuring that every published manuscript has an added value and thus providing international acceptance of the “readability” of the manuscripts published in the journal.