NAC4ED: A high-throughput computational platform for the rational design of enzyme activity and substrate selectivity

Abstract

In silico computational methods have been widely utilized to study enzyme catalytic mechanisms and design enzyme performance, including molecular docking, molecular dynamics, quantum mechanics, and multiscale QM/MM approaches. However, the manual operation associated with these methods poses challenges for simulating enzymes and enzyme variants in a high-throughput manner. We developed the NAC4ED, a high-throughput enzyme mutagenesis computational platform based on the “near-attack conformation” design strategy for enzyme catalysis substrates. This platform circumvents the complex calculations involved in transition-state searching by representing enzyme catalytic mechanisms with parameters derived from near-attack conformations. NAC4ED enables the automated, high-throughput, and systematic computation of enzyme mutants, including protein model construction, complex structure acquisition, molecular dynamics simulation, and analysis of active conformation populations. Validation of the accuracy of NAC4ED demonstrated a prediction accuracy of 92.5% for 40 mutations, showing strong consistency between the computational predictions and experimental results. The time required for automated determination of a single enzyme mutant using NAC4ED is 1/764th of that needed for experimental methods. This has significantly enhanced the efficiency of predicting enzyme mutations, leading to revolutionary breakthroughs in improving the performance of high-throughput screening of enzyme variants. NAC4ED facilitates the efficient generation of a large amount of annotated data, providing high-quality data for statistical modeling and machine learning. NAC4ED is currently available at http://lujialab.org.cn/software/.

Keywords

high-throughput screening / near-attack conformation / protein engineering / rational design

References

[1] Hauer B. Embracing nature’s catalysts: a viewpoint on the future of biocatalysis. ACS Catal. 2020;10:8418–8427.

[2] Winkler CK, Schrittwieser JH, Kroutil W. Power of biocatalysis for organic synthesis. ACS Cent Sci. 2021;7:55–71.

[3] Chen K, Arnold FH. Engineering new catalytic activities in enzymes. Nat Catal. 2020;3:203–213.

[4] Nödling AR, Santi N, Williams TL, Tsai YH, Luk LYP. Enabling protein-hosted organocatalytic transformations. RSC Adv. 2020;10:16147–16161.

[5] CHEN K, Arnold FH. Tuning the activity of an enzyme for unusual environments: sequential random mutagenesis of subtilisin E for catalysis in dimethylformamide. Proc Natl Acad Sci USA. 1993;90:5618–5622.

[6] Arango Gutierrez E, Mundhada H, Meier T, Duefel H, Bocola M, Schwaneberg U. Reengineered glucose oxidase for amperometric glucose determination in diabetes analytics. Biosens Bioelectron. 2013;50:84–90.

[7] Stemmer WPC. Rapid evolution of a protein in vitro by DNA shuffling. Nature. 1994;370:389–391.

[8] Wrenbeck EE, Azouz LR, Whitehead TA. Single-mutation fitness landscapes for an enzyme on multiple substrates reveal specificity is globally encoded. Nat Commun. 2017;8:15695.

[9] Bornscheuer UT, Huisman GW, Kazlauskas RJ, Lutz S, Moore JC, Robins K. Engineering the third wave of biocatalysis. Nature. 2012;485:185–194.

[10] Doerr S, Harvey MJ, Noé F, De Fabritiis G. HTMD: high-throughput molecular dynamics for molecular discovery. J Chem Theory Comput. 2016;12:1845–1852.

[11] Kiss G, Çelebi-Ölçüm N, Moretti R, Baker D, Houk KN. Computational enzyme design. Angew Chem Int Ed. 2013;52:5700–5725.

[12] Bunzel HA, Garrabou X, Pott M, Hilvert D. Speeding up enzyme discovery and engineering with ultrahigh-throughput methods. Curr Opin Struct Biol. 2018;48:149–156.

[13] Cui Q, Karplus M. Molecular properties from combined QM/MM methods. I. Analytical second derivative and vibrational calculations. J Chem Phys. 2000;112:1133–1149.

[14] Lightstone FC, Bruice TC. Ground state conformations and entropic and enthalpic factors in the efficiency of intramolecular and enzymatic reactions. 1. Cyclic anhydride formation by substituted glutarates, succinate, and 3,6-endoxo-Δ₄-tetrahydrophthalate monophenyl esters. J Am Chem Soc. 1996;118:2595–2605.

[15] Bruice TC. A view at the millennium: the efficiency of enzymatic catalysis. Acc Chem Res. 2002;35:139–148.

[16] Hur S, Bruice TC. The near attack conformation approach to the study of the chorismate to prephenate reaction. Proc Natl Acad Sci USA. 2003;100:12015–12020.

[17] Boral S, Sen S, Kushwaha T, Inampudi KK, De S. Extein residues regulate the catalytic function of Spl DnaX intein enzyme by restricting the near-attack conformations of the active-site residues. Prot Sci. 2023;32:e4699.

[18] Yang W, Zhuang J, Li C, Cheng GJ. Unveiling the methyl transfer mechanisms in the epigenetic machinery DNMT3A-3 L: a comprehensive study integrating assembly dynamics with catalytic reactions. Comput Struct Biotechnol J. 2023;21:2086–2099.

[19] Ji S, Luo L, Li C, Liu M, Liu Y, Jiang L. Rational modulation of the enzymatic intermediates for tuning the phosphatase activity of histidine kinase HK853. Biochem Biophys Res Commun. 2020;523:733–738.

[20] Sinko G. Modeling of a near-attack conformation of oxime in phosphorylated acetylcholinesterase via a reactivation product, a phosphorylated oxime. Chem Biol Interact. 2023;383:110656.

[21] Ashworth MA, Bombino E, de Jong RM, Wijma HJ, Janssen DB, McLean KJ, et al. Computation-aided engineering of cytochrome P450 for the production of pravastatin. ACS Catal. 2022;12:15028–15044.

[22] Vieira LA, Almeida JSFD, De Koning MC, LaPlante SR, Borges Jr. I, et al. Molecular modeling of Mannich phenols as reactivators of human acetylcholinesterase inhibited by A-series nerve agents. Chem Biol Interact. 2023;382:110622.

[23] Cicek E, Monard G, Sungur FA. Molecular mechanism of protein arginine deiminase 2: a study involving multiple microsecond long molecular dynamics simulations. Biochemistry. 2022;61:1286–1297.

[24] Ramírez-Palacios C, Wijma HJ, Thallmair S, Marrink SJ, Janssen DB. Computational prediction of ω-transaminase specificity by a combination of docking and molecular dynamics simulations. J Chem Inf Model. 2021;61:5569–5580.

[25] Doerr M, Romero A, Daza MC. Effect of the acyl-group length on the chemoselectivity of the lipase-catalyzed acylation of propranolol—a computational study. J Mol Model. 2021;27:198.

[26] Liu Y, Xu G, Zhou J, Ni J, Zhang L, Hou X, et al. Structure-guided engineering of d-carbamoylase reveals a key loop at substrate entrance tunnel. ACS Catal. 2020;10:12393–12402.

[27] Wang X, Jing X, Deng Y, Nie Y, Xu F, Xu Y, et al. Evolutionary coupling saturation mutagenesis: coevolution-guided identification of distant sites influencing Bacillus naganoensis pullulanase activity. FEBS Lett. 2020;594:799–812.

[28] Cooper AM, Kästner J. Averaging techniques for reaction barriers in QM/MM simulations. Chemphyschem. 2014;15:3264–3269.

[29] Ryde U. How many conformations need to be sampled to obtain converged QM/MM energies? The curse of exponential averaging. J Chem Theory Comput. 2017;13:5745–5752.

[30] Xu T, Gao B, Zhang L, Lin J, Wang X, Wei D. Template-based modeling of a psychrophilic lipase: conformational changes, novel structural features and its application in predicting the enantioselectivity of lipase catalyzed transesterification of secondary alcohols. Biochim Biophys Acta. 2010;1804:2183–2190.

[31] Shi T, Liu L, Tao W, Luo S, Fan S, Wang X-L, et al. Theoretical studies on the catalytic mechanism and substrate diversity for macrocyclization of pikromycin thioesterase. ACS Catal. 2018;8:4323–4332.

[32] Chen D, Li Y, Li X, Hong X, Fan X, Savidge T. Key difference between transition state stabilization and ground state destabilization: increasing atomic charge densities before or during enzyme–substrate binding. Chem Sci. 2022;13:8193–8202.

[33] Song Z, Zhang Q, Wu W, Pu Z, Yu H. Rational design of enzyme activity and enantioselectivity. Front Bioeng Biotechnol. 2023;11:1129149.

[34] Reetz MT, Wang LW, Bocola M. Directed evolution of enantioselective enzymes: iterative cycles of CASTing for probing protein-sequence space. Angew Chem Int Ed. 2006;45:1236–1241.

[35] Reetz MT, Sanchis J. Constructing and analyzing the fitness landscape of an experimental evolutionary process. ChemBioChem. 2008;9:2260–2267.

[36] Reetz MT, Bocola M, Wang LW, Sanchis J, Cronin A, Arand M, et al. Directed evolution of an enantioselective epoxide hydrolase: uncovering the source of enantioselectivity at each evolutionary stage. J Am Chem Soc. 2009;131:7334–7343.

[37] Berman HM, Battistuz T, Bhat TN, Bluhm WF, Bourne PE, Burkhardt K, et al. The protein data bank. Acta Crystallogr D. 2002;58:899–907.

[38] Jumper J, Evans R, Pritzel A, Green T, Figurnov M, Ronneberger O, et al. Highly accurate protein structure prediction with AlphaFold. Nature. 2021;596:583–589.

[39] Baek M, DiMaio F, Anishchenko I, Dauparas J, Ovchinnikov S, Lee GR, et al. Accurate prediction of protein structures and interactions using a three-track neural network. Science. 2021;373:871–876.

[40] Bianco G, Holcomb M, Santos-Martins D, Tillack A, Hansel-Harris A, Forli S. Reactive docking: a computational method for high-throughput virtual screenings of reactive species. J Chem Inf Model. 2023;63:5631–5640.

[41] Shen HC, Hammock BD. Discovery of inhibitors of soluble epoxide hydrolase: a target with multiple potential therapeutic indications. J Med Chem. 2012;55:1789–1808.

[42] Johnston RC, Yao K, Kaplan Z, Chelliah M, Leswing K, Seekins S, et al. Epik: pK(a) and protonation state prediction through machine learning. J Chem Theory Comput. 2023;19:2380–2388.

[43] Friesner RA, Banks JL, Murphy RB, Halgren TA, Klicic JJ, Mainz DT, et al. Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. J Med Chem. 2004;47:1739–1749.

[44] Reetz MT. Recent advances in directed evolution of stereoselective enzymes. In: Alcalde M, editor. Directed enzyme evolution: Advances and applications. Cham: Springer; 2017. p. 69–99.

[45] Li D, Partin AC, Zhao L, Chen I, Michaels ML, Wang Z, et al. Protocol for high-throughput cloning, expression, purification, and evaluation of bispecific antibodies. STAR Protocols. 2022;3:101428.

[46] Bezsudnova EY, Popov VO, Boyko KM. Structural insight into the substrate specificity of PLP fold type IV transaminases. Appl Microbiol Biotechnol. 2020;104:2343–2357.

[47] Xiang C, Ao YF, Höhne M, Bornscheuer UT. Shifting the pH optima of (R)-selective transaminases by protein engineering. Int J Mol Sci. 2022;23:15347.

RIGHTS & PERMISSIONS

2024 The Author(s). mLife published by John Wiley & Sons Australia, Ltd on behalf of Institute of Microbiology, Chinese Academy of Sciences.