Computational design of 5-carboxyl-2-pentenoyl-CoA reductase from Thermobifida fusca to enhance adipic acid production
YANG Ju1, MAO Yin1, HUANG Xiaoqiang2, ZHOU Shenghu1, DENG Yu1*
1(School of Bioengineering,Jiangnan University,Wuxi 214122,China) 2(Department of Computational Medicine & Bioinformatics,University of Michigan,Ann Arbor 48109,USA)
Abstract: In order to enhance the biosynthesis of adipic acid, the active pocket of 5-carboxy-2-pentenoyl-CoA reductase was modified using computational enzyme design.Based on the substrate binding model,fourteen residues including Ser 88,Leu 89,Ile 90,Pro 91,Ala 92,Val 93,Lys 95,Leu 96,Thr 161,Thr 246,Thr 249,Ile 250,Gln 253 and Tyr 367 were designed to improve the combination between the enzyme and 5-carboxy-2-pentenoyl-CoA and the catalytic activity of 5-carboxy-2-pentenoyl-CoA reductase by introducing hydrogen bonds network.In the 10 designs (Des0-Des9), the carboxyl of 5-carboxy-2-pentenoyl-CoA could form hydrogen bonds with Gln253Arg and Ile250Gln in Des0, Des3, Des4 and Des9.And Gln253Arg could form hydrogen bonds with Leu89Ser (Des0 and Des9) or Leu89Thr (Des3 and Des4) with Thr364.In order to test these designs, the stability of hydrogen bonds were analyzed by molecular dynamics simulation.The results showed that the four hydrogen bonds designed in Des0 were stable during the process of 16 ns molecular dynamics simulation.It indicated that Des0 may have a strong binding effect with 5-carboxy-2-pentenoyl-CoA.Accordingly, we speculated that Des0 could improve the catalytic activity of 5-carboxy-2-pentenoyl-CoA reductase, which should be verified by subsequent experiments.
BABU T,YUN E J,KIM S,et al.Engineering Escherichia coli for the production of adipic acid through the reversed β-oxidation pathway[J].Process Biochemistry,2015,50:2 066-2 071.
[2]
NOACK H,GEORGIEV V,BLOMBERG M R A,et al.Theoretical insights into heme-catalyzed oxidation of cyclohexane to adipic acid[J].Inorganic Chemistry,2011,50:1 194-1 202.
[3]
LI X K,WU D,LU T,et al.Highly efficient chemical process to convert mucic acid into adipic acid and DFT studies of the mechanism of the rhenium-catalyzed deoxydehydration[J].Angewandte Chemie,2014,53:4 200-4 204.
[4]
BLACH P,BÖSTROM Z,FRANCESCHI-MESSANT S,et al.Recyclable process for sustainable adipic acid production in microemulsions[J].Tetrahedron,2010,66:7 124-7 128.
[5]
SATO K,AOKI M,NOYORI R.A “green” route to adipic acid:Direct oxidation of cyclohexenes with 30 percent hydrogen peroxide[J].Science,1998,281:1 646-1 647.
[6]
NIU W,DRATHS K M,FROST J W.Benzene-free synthesis of adipic acid[J].Biotechnology Progress,2002,18:201-211.
[7]
DENG Y,MAO Y.Production of adipic acid by the native-occurring pathway in Thermobifida fusca B6[J].Journal of Applied Microbiology,2015,119:1 057-1 063.
[8]
ZHAO M,HUANG D X,ZHANG X J,et al.Metabolic engineering of Escherichia coli for producing adipic acid through the reverse adipate-degradation pathway[J].Metabolic Engineering,2018,47:254-262.
[9]
TIAN Y,HUANG X Q,LI Q,et al.Computational design of variants for cephalosporin C acylase from Pseudomonas strain N176 with improved stability and activity[J].Applied Microbiology Biotechnology,2017,101:621-632.
[10]
HUANG X Q,HAN K H,ZHU Y S.Systematic optimization model and algorithm for binding sequence selection in computational enzyme design[J].Protein Science,2013,22:929-941.
[11]
HE J W,HUANG X Q,XUE J,et al.Computational redesign of penicillin acylase for cephradine synthesis with high kinetic selectivity[J].Green Chemistry,2018,20:5 484-5 490.
[12]
HUANG X Q,XUE J,ZHU Y S.Computational design of cephradine synthase in a new scaffold identified from structural databases[J].Chemical Communications,2017,53:7 604-7 607.
[13]
HUANG X Q,YANG J,ZHU Y S.A solvated ligand rotamer approach and its application in computational protein design[J].Journal of Molecular Modeling,2013,19:1 355-1 367.
[14]
TIAN Y,XU Z B,HUANG X Q,et al.Computational design to improve catalytic activity of cephalosporin C acylase from Pseudomonas strain N176[J].RSC Advances,2017,7:30 370-30 375.
[15]
JIANG L,ALTHOFF E A,Clemente F R,et al.De novo computational design of retro-aldol enzymes[J].Science,2008,319:1 387-1 391.
[16]
RÖTHLISBERGER D,KHERSONSKY O,WOLLACOTT A M,et al.Kemp elimination catalysts by computational enzyme design[J].Nature,2008,453:190-195.
[17]
SIEGEL J B,ZANGHELLINI A,LOVICK H M,et al.Computational design of an enzyme catalyst for a stereoselective bimolecular Diels-Alder reaction[J].Science,2010,329:309-313.
[18]
BERMAN H M,WESTBROOK,FENG Z K,et al.The protein data bank[J].Nucleic Acids Research,2000,28:235-242.
[19]
ALTSCHUL S F,MADDEN T L,SCHÄFFER A A,et al.Gapped BLAST and PSI-BLAST:A new generation of protein database search programs[J].Nucleic Acids Research,1997,25:3 389-3 402.
[20]
BIASINI M,BIENERT S,WATERHOUSE A,et al.SWISS-MODEL:Modelling protein tertiary and quaternary structure using evolutionary information[J].Nucleic Acids Research,2014,42:W252-W258.
[21]
CHOWDHURY N P,MOWAFY A M,DEMMER J K,et al.Studies on the mechanism of electron bifurcation catalyzed by electron transferring flavoprotein (Etf) and butyryl-CoA dehydrogenase (Bcd) of Acidaminococcus fermentans[J].The Journal of Biological Chemistry,2014,289(8):5 145-5 157.
[22]
ABRAHAM M J,MURTOLA T,SCHULZ R,et al.GROMACS:High performance molecular simulations through multi-level parallelism from laptops to supercomputers[J].SoftwareX,2015,1:19-25.
[23]
YANG J,WEI Y F,LI G H,et al.Computer-aided engineering of adipyl-CoA synthetase for enhancing adipic acid synthesis[J].Biotechnology Letters,2020,42:2 693-2 701.
[24]
XUE J,HUANG X Q,ZHU Y S.Using molecular dynamics simulations to evaluate active designs of cephradine hydrolase by molecular mechanics/Poisson-Boltzmann surface area and molecular mechanics/generalized Born surface area methods[J].RSC Advances,2019,9:13 868-13 877.
[25]
LI Q,HUANG X Q,ZHU Y S.Evaluation of active designs of cephalosporin C acylase by molecular dynamics simulation and molecular docking[J].Journal of molecular modeling,2014,20:2 314.
[26]
HUANG X Q,PEARCE R,ZHANG Y.EvoEF2:Accurate and fast energy function for computational protein design[J].Bioinformatics,2020,36:1 135-1 142.
[27]
ABENDROTH J,GARDBERG A S,ROBINSON J I,et al.SAD phasing using iodide ions in a high-throughput structural genomics environment[J].Journal of Structural And Functional Genomics,2011,12:83-95.
[28]
PEARCE R,HUANG X Q,SETIAWAN D,et al.EvoDesign:Designing protein-protein binding interactions using evolutionary interface profiles in conjunction with an optimized physical energy function[J].Journal of Molecular Biology,2019,431:2 467-2 476.
[29]
SHAPOVALOV M V,DUNBRACK R L.A smoothed backbone-dependent rotamer library for proteins derived from adaptive kernel density estimates and regressions[J].Structure,2011,19(6):844-858.