-
Protein function is related to its structure and dynamic change. Molecular dynamics simulation is an important tool for studying protein dynamics by exploring its conformational space, however, conformational sampling is a nontrivial issue, because of the risk of missing key details during sampling. In recent years, deep learning methods, such as auto-encoder, can couple with MD to explore conformational space of protein. After being trained with the MD trajectories, auto-encoder can generate new conformations quickly by inputting random numbers in low dimension space. However, some problems still exist, such as requirements for the quality of the training set, the limitation of explorable area and the undefined sampling direction. In this work, we build a supervised auto-encoder, in which some reaction coordinates are used to guide conformational exploration along certain directions. We also try to expand the explorable area by training through the data generated by the model. Two multi-domain proteins, bacteriophage T4 lysozyme and adenylate kinase, are used to illustrate the method. In the case of the training set consisting of only under-sampled simulated trajectories, the supervised auto-encoder can still explore along the given reaction coordinates. The explored conformational space can cover all the experimental structures of the proteins and be extended to regions far from the training sets. Having been verified by molecular dynamics and secondary structure calculations, most of the conformations explored are found to be plausible. The supervised auto-encoder provides a way to efficiently expand the conformational space of a protein with limited computational resources, although some suitable reaction coordinates are required. By integrating appropriate reaction coordinates or experimental data, the supervised auto-encoder may serve as an efficient tool for exploring conformational space of proteins.
-
Keywords:
- protein conformational space /
- molecular dynamics simulation /
- machine learning /
- auto-encoder
[1] Chu X, Gan L, Wang E, Wang J 2013 Proc. Natl. Acad. Sci. U.S.A. 110 E2342Google Scholar
[2] Smyth M S, Martin J H 2000 Mol. Pathol. 53 8Google Scholar
[3] Danev R, Yanagisawa H, Kikkawa M 2019 Trends Biochem. Sci. 44 837Google Scholar
[4] Vincenzi M, Mercurio F A, Leone M 2021 Curr. Med. Chem. 28 2729Google Scholar
[5] Kachala M, Valentini E, Svergun D I 2015 Adv. Exp. Med. Biol. 870 261Google Scholar
[6] Chu F, Thornton D T, Nguyen H T 2018 Methods 144 53Google Scholar
[7] Bhaumik S R 2021 Emerg. Top Life Sci. 5 49Google Scholar
[8] Jumper J, Evans R, Pritzel A, Green T, Figurnov M, Ronneberger O, Tunyasuvunakool K, Bates R, Žídek A, Potapenko A, Bridgland A, Meyer C, Kohl S A A, Ballard A J, Cowie A, Romera-Paredes B, Nikolov S, Jain R, Adler J, Back T, Petersen S, Reiman D, Clancy E, Zielinski M, Steinegger M, Pacholska M, Berghammer T, Bodenstein S, Silver D, Vinyals O, Senior A W, Kavukcuoglu K, Kohli P, Hassabis D 2021 Nature 596 583Google Scholar
[9] Baek M, DiMaio F, Anishchenko I, Dauparas J, Ovchinnikov S, Lee G R, Wang J, Cong Q, Kinch L N, Schaeffer R D, Millán C, Park H, Adams C, Glassman C R, DeGiovanni A, Pereira J H, Rodrigues A V, van Dijk A A, Ebrecht A C, Opperman D J, Sagmeister T, Buhlheller C, Pavkov-Keller T, Rathinaswamy M K, Dalwadi U, Yip C K, Burke J E, Garcia K C, Grishin N V, Adams P D, Read R J, Baker D 2021 Science 373 871Google Scholar
[10] Karplus M, Kuriyan J 2005 Proc. Natl. Acad. Sci. 102 6679Google Scholar
[11] Bernardi R C, Melo M C R, Schulten K 2015 Biochim. Biophys. Acta 1850 872Google Scholar
[12] Mu J, Liu H, Zhang J, Luo R, Chen H F 2021 J. Chem. Inf. Model. 61 1037Google Scholar
[13] Lemke T, Peter C 2019 J. Chem. Theory Comput. 15 1209Google Scholar
[14] Zhu J, Wang J, Han W, Xu D 2022 Nat. Commun. 13 1661Google Scholar
[15] Hinton G E, Salakhutdinov R R 2006 Science 313 504Google Scholar
[16] Degiacomi M T 2019 Structure 27 1034Google Scholar
[17] Wen B, Peng J, Zuo X, Gong Q, Zhang Z 2014 Biophysical J. 107 956Google Scholar
[18] Giri Rao V V H, Gosavi S 2014 PLOS Computational Biology 10 e1003938Google Scholar
[19] Abraham M J, Murtola T, Schulz R, Páll S, Smith J C, Hess B, Lindahl E 2015 SoftwareX 1–2 19Google Scholar
[20] Weaver L H, Matthews B W 1987 J. Mol. Biol. 193 189Google Scholar
[21] Zhang X J, Wozniak J A, Matthews B W 1995 J. Mol. Biol. 250 527Google Scholar
[22] Müller C W, Schulz G E 1992 J. Mol. Biol. 224 159Google Scholar
[23] Müller C W, Schlauderer G J, Reinstein J, Schulz G E 1996 Structure 4 147Google Scholar
[24] Hornak V, Abel R, Okur A, Strockbine B, Roitberg A, Simmerling C 2006 Proteins Struct. Funct. Bioinf. 65 712Google Scholar
[25] Izadi S, Anandakrishnan R, Onufriev A V 2014 J. Phys. Chem. Lett. 5 3863Google Scholar
[26] Huang J, Rauscher S, Nawrocki G, Ran T, Feig M, de Groot B L, Grubmüller H, MacKerell A D 2017 Nat. Methods 14 71Google Scholar
[27] Bussi G, Donadio D, Parrinello M 2007 J. Chem. Phys. 126 014101Google Scholar
[28] Essmann U, Perera L E, Berkowitz M L, Darden T A, Lee H C, Pedersen L G 1995 J. Chem. Phys. 103 8577Google Scholar
[29] Kingma D P, Ba J 2014 arXiv:1412.6980 [cs.LG
[30] Lovell S C, Davis I W, Arendall III W B, de Bakker P I W, Word J M, Prisant M G, Richardson J S, Richardson D C 2003 Proteins Struct. Funct. Bioinf. 50 437Google Scholar
[31] Eastman P, Swails J, Chodera J D, McGibbon R T, Zhao Y, Beauchamp K A, Wang L P, Simmonett A C, Harrigan M P, Stern C D, Wiewiora R P, Brooks B R, Pande V S 2017 PLoS Comput. Biol. 13 e1005659Google Scholar
[32] Shirts M R, Klein C, Swails J M, Yin J, Gilson M K, Mobley D L, Case D A, Zhong E D 2017 J. Comput. -Aided Mol. Des. 31 147Google Scholar
[33] Touw W G, Baakman C, Black J, te Beek T A, Krieger E, Joosten R P, Vriend G 2015 Nucleic Acids Res. 43 D364Google Scholar
-
图 2 本研究中使用的两种蛋白质分子的不同结构 (a) T4L的闭合(不透明)和打开(透明)结构, 紫色为α螺旋, 黄色为β折叠; (b) AdK的闭合(不透明)和打开(透明)结构, 不同颜色表示不同的结构域
Figure 2. Different structures of the two proteins in the work. (a) The close (opaque) and open (transparent) state of T4L. α-helix is colored in purple and β-sheet is colored in yellow. (b) The close (opaque) and open (transparent) state of AdK. Different domains are colored in different colors.
图 4 探索到的不同T4L构象 (a) PDB编号173L的晶体结构(不透明)与探索到的相似结构(透明); (b) 开合程度不同的两个构象; (c) 扭动情况不同的两个构象; 紫色为α螺旋, 黄色为β折叠
Figure 4. Different T4L conformations explored: (a) PDB:173L (opaque) and a similar structure explored; (b) two conformations with different degrees of opening and closing; (c) two conformations with different degrees of twisting. α-helix is colored in purple and β-sheet is colored in yellow.
图 5 T4L构象探索结果的合理性检验 (a) 使用AMBER99SB力场/OPC水模型; (b) 使用CHARMM36m力场/TIP3P水模型; (c) 修复后各代表构象的二级结构含量, 参考值为模拟轨迹的平均值
Figure 5. Plausibility check of T4L conformational exploration results: (a) With AMBER99SB/OPC; (b) with CHARMM36m/TIP3P; (c) secondary structure counts of each representative conformation after fixing, the reference is the average value of the simulated trajectory.
图 9 AdK构象探索结果的合理性检验 (a) 使用AMBER99SB力场/OPC水模型; (b)使用CHARMM36m力场/TIP3P水模型; (c) 修复后各代表构象的二级结构含量, 参考值为模拟轨迹的平均值
Figure 9. Plausibility check of AdK conformational exploration results: (a) With AMBER99SB/OPC; (b) with CHARMM36m/TIP3P; (c) secondary structure counts of each representative conformation after fixing, the reference is the average value of the simulated trajectory.
-
[1] Chu X, Gan L, Wang E, Wang J 2013 Proc. Natl. Acad. Sci. U.S.A. 110 E2342Google Scholar
[2] Smyth M S, Martin J H 2000 Mol. Pathol. 53 8Google Scholar
[3] Danev R, Yanagisawa H, Kikkawa M 2019 Trends Biochem. Sci. 44 837Google Scholar
[4] Vincenzi M, Mercurio F A, Leone M 2021 Curr. Med. Chem. 28 2729Google Scholar
[5] Kachala M, Valentini E, Svergun D I 2015 Adv. Exp. Med. Biol. 870 261Google Scholar
[6] Chu F, Thornton D T, Nguyen H T 2018 Methods 144 53Google Scholar
[7] Bhaumik S R 2021 Emerg. Top Life Sci. 5 49Google Scholar
[8] Jumper J, Evans R, Pritzel A, Green T, Figurnov M, Ronneberger O, Tunyasuvunakool K, Bates R, Žídek A, Potapenko A, Bridgland A, Meyer C, Kohl S A A, Ballard A J, Cowie A, Romera-Paredes B, Nikolov S, Jain R, Adler J, Back T, Petersen S, Reiman D, Clancy E, Zielinski M, Steinegger M, Pacholska M, Berghammer T, Bodenstein S, Silver D, Vinyals O, Senior A W, Kavukcuoglu K, Kohli P, Hassabis D 2021 Nature 596 583Google Scholar
[9] Baek M, DiMaio F, Anishchenko I, Dauparas J, Ovchinnikov S, Lee G R, Wang J, Cong Q, Kinch L N, Schaeffer R D, Millán C, Park H, Adams C, Glassman C R, DeGiovanni A, Pereira J H, Rodrigues A V, van Dijk A A, Ebrecht A C, Opperman D J, Sagmeister T, Buhlheller C, Pavkov-Keller T, Rathinaswamy M K, Dalwadi U, Yip C K, Burke J E, Garcia K C, Grishin N V, Adams P D, Read R J, Baker D 2021 Science 373 871Google Scholar
[10] Karplus M, Kuriyan J 2005 Proc. Natl. Acad. Sci. 102 6679Google Scholar
[11] Bernardi R C, Melo M C R, Schulten K 2015 Biochim. Biophys. Acta 1850 872Google Scholar
[12] Mu J, Liu H, Zhang J, Luo R, Chen H F 2021 J. Chem. Inf. Model. 61 1037Google Scholar
[13] Lemke T, Peter C 2019 J. Chem. Theory Comput. 15 1209Google Scholar
[14] Zhu J, Wang J, Han W, Xu D 2022 Nat. Commun. 13 1661Google Scholar
[15] Hinton G E, Salakhutdinov R R 2006 Science 313 504Google Scholar
[16] Degiacomi M T 2019 Structure 27 1034Google Scholar
[17] Wen B, Peng J, Zuo X, Gong Q, Zhang Z 2014 Biophysical J. 107 956Google Scholar
[18] Giri Rao V V H, Gosavi S 2014 PLOS Computational Biology 10 e1003938Google Scholar
[19] Abraham M J, Murtola T, Schulz R, Páll S, Smith J C, Hess B, Lindahl E 2015 SoftwareX 1–2 19Google Scholar
[20] Weaver L H, Matthews B W 1987 J. Mol. Biol. 193 189Google Scholar
[21] Zhang X J, Wozniak J A, Matthews B W 1995 J. Mol. Biol. 250 527Google Scholar
[22] Müller C W, Schulz G E 1992 J. Mol. Biol. 224 159Google Scholar
[23] Müller C W, Schlauderer G J, Reinstein J, Schulz G E 1996 Structure 4 147Google Scholar
[24] Hornak V, Abel R, Okur A, Strockbine B, Roitberg A, Simmerling C 2006 Proteins Struct. Funct. Bioinf. 65 712Google Scholar
[25] Izadi S, Anandakrishnan R, Onufriev A V 2014 J. Phys. Chem. Lett. 5 3863Google Scholar
[26] Huang J, Rauscher S, Nawrocki G, Ran T, Feig M, de Groot B L, Grubmüller H, MacKerell A D 2017 Nat. Methods 14 71Google Scholar
[27] Bussi G, Donadio D, Parrinello M 2007 J. Chem. Phys. 126 014101Google Scholar
[28] Essmann U, Perera L E, Berkowitz M L, Darden T A, Lee H C, Pedersen L G 1995 J. Chem. Phys. 103 8577Google Scholar
[29] Kingma D P, Ba J 2014 arXiv:1412.6980 [cs.LG
[30] Lovell S C, Davis I W, Arendall III W B, de Bakker P I W, Word J M, Prisant M G, Richardson J S, Richardson D C 2003 Proteins Struct. Funct. Bioinf. 50 437Google Scholar
[31] Eastman P, Swails J, Chodera J D, McGibbon R T, Zhao Y, Beauchamp K A, Wang L P, Simmonett A C, Harrigan M P, Stern C D, Wiewiora R P, Brooks B R, Pande V S 2017 PLoS Comput. Biol. 13 e1005659Google Scholar
[32] Shirts M R, Klein C, Swails J M, Yin J, Gilson M K, Mobley D L, Case D A, Zhong E D 2017 J. Comput. -Aided Mol. Des. 31 147Google Scholar
[33] Touw W G, Baakman C, Black J, te Beek T A, Krieger E, Joosten R P, Vriend G 2015 Nucleic Acids Res. 43 D364Google Scholar
Catalog
Metrics
- Abstract views: 2960
- PDF Downloads: 184
- Cited By: 0