1
|
Wakabayashi T, Oide M, Nakasako M. CryoEM-sampling of metastable conformations appearing in cofactor-ligand association and catalysis of glutamate dehydrogenase. Sci Rep 2024; 14:11165. [PMID: 38750092 PMCID: PMC11096400 DOI: 10.1038/s41598-024-61793-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2023] [Accepted: 05/09/2024] [Indexed: 05/18/2024] Open
Abstract
Kinetic aspects of enzymatic reactions are described by equations based on the Michaelis-Menten theory for the initial stage. However, the kinetic parameters provide little information on the atomic mechanism of the reaction. In this study, we analyzed structures of glutamate dehydrogenase in the initial and steady stages of the reaction using cryoEM at near-atomic resolution. In the initial stage, four metastable conformations displayed different domain motions and cofactor/ligand association modes. The most striking finding was that the enzyme-cofactor-substrate complex, treated as a single state in the enzyme kinetic theory, comprised at least three different metastable conformations. In the steady stage, seven conformations, including derivatives from the four conformations in the initial stage, made the reaction pathway complicated. Based on the visualized conformations, we discussed stage-dependent pathways to illustrate the dynamics of the enzyme in action.
Collapse
Grants
- JPMJPR22E2 Japan Science and Technology Agency
- jp13480214 Japan Society for the Promotion of Science
- jp19204042 Japan Society for the Promotion of Science
- jp22244054 Japan Society for the Promotion of Science
- jp21H01050 Japan Society for the Promotion of Science
- jp26800227 Japan Society for the Promotion of Science
- 18J11653 Japan Society for the Promotion of Science
- jp15076210 Ministry of Education, Culture, Sports, Science and Technology of Japan
- jp20050030 Ministry of Education, Culture, Sports, Science and Technology of Japan
- jp22018027 Ministry of Education, Culture, Sports, Science and Technology of Japan
- jp23120525, jp25120725 Ministry of Education, Culture, Sports, Science and Technology of Japan
- jp15H01647 Ministry of Education, Culture, Sports, Science and Technology of Japan
- jp17H05891 Ministry of Education, Culture, Sports, Science and Technology of Japan
Collapse
Affiliation(s)
- Taiki Wakabayashi
- Department of Physics, Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoko-Ku, Yokohama, Kanagawa, 223-8522, Japan
- RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-Cho, Sayo-Gun, Hyogo, 679-5148, Japan
| | - Mao Oide
- Department of Physics, Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoko-Ku, Yokohama, Kanagawa, 223-8522, Japan
- RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-Cho, Sayo-Gun, Hyogo, 679-5148, Japan
- PRESTO, Japan Science and Technology Agency, Chiyoda-Ku, Tokyo, 102-0076, Japan
- Protein Research Institute, Osaka University, Yamadaoka, Suita, Osaka, 565-0871, Japan
| | - Masayoshi Nakasako
- Department of Physics, Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoko-Ku, Yokohama, Kanagawa, 223-8522, Japan.
- RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-Cho, Sayo-Gun, Hyogo, 679-5148, Japan.
| |
Collapse
|
2
|
Ichinomiya T. Topological data analysis gives two folding paths in HP35(nle-nle), double mutant of villin headpiece subdomain. Sci Rep 2022; 12:2719. [PMID: 35177744 PMCID: PMC8854739 DOI: 10.1038/s41598-022-06682-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2021] [Accepted: 02/04/2022] [Indexed: 11/16/2022] Open
Abstract
The folding dynamics of proteins is a primary area of interest in protein science. We carried out topological data analysis (TDA) of the folding process of HP35(nle-nle), a double-mutant of the villin headpiece subdomain. Using persistent homology and non-negative matrix factorization, we reduced the dimension of protein structure and investigated the flow in the reduced space. We found this protein has two folding paths, distinguished by the pairings of inter-helix residues. Our analysis showed the excellent performance of TDA in capturing the formation of tertiary structure.
Collapse
Affiliation(s)
- Takashi Ichinomiya
- Department of Systems Biology, Gifu University School of Medicine, Yanagido 1-1, Gifu, 501-1194, Japan. .,The United Graduate School of Drug Discovery and Medical Information Sciences of Gifu University, Yanagido 1-1, Gifu, 501-1194, Japan.
| |
Collapse
|
3
|
Tao P, Wang E, Xiao Y. Pathway regulation mechanism revealed by cotranslational folding of villin headpiece subdomain HP35. Phys Rev E 2021; 101:052403. [PMID: 32575289 DOI: 10.1103/physreve.101.052403] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2019] [Accepted: 04/13/2020] [Indexed: 01/07/2023]
Abstract
Cotranslational folding is one of the most important features of protein folding in vivo. Although many studies have shown that the folding pathways of cotranslational folding are different from free folding in vitro, the detailed mechanism of folding dynamics is lacking. Here we combine all-atom molecular simulations with an ideal ribosome tunnel model to investigate the cotranslational folding of villin headpiece subdomain HP35. By comparing the folding dynamics between cotranslational folding and free folding, we found that cotranslational folding tends to fold along the pathway that is easier to fold into native state in the latter. In addition, the roles of the ribosome tunnel and sequential folding are analyzed separately. Our results show that the ribosome can prevent the untimely folding of the C segment of HP35 to reduce the non-native interactions, while the translation speed can regulate the amounts of native and non-native interactions and the balance between them. Overall, these results give insights into the general mechanisms of cotranslational protein folding.
Collapse
Affiliation(s)
- Peng Tao
- School of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| | - Ercheng Wang
- School of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China.,College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Yi Xiao
- School of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| |
Collapse
|
4
|
Yasuda T, Shigeta Y, Harada R. The Folding of Trp-cage is Regulated by Stochastic Flip of the Side Chain of Tryptophan. CHEM LETT 2021. [DOI: 10.1246/cl.200699] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022]
Affiliation(s)
- Takunori Yasuda
- College of Biological Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-0821, Japan
| | - Yasuteru Shigeta
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-0821, Japan
| | - Ryuhei Harada
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-0821, Japan
| |
Collapse
|
5
|
Heilmann N, Wolf M, Kozlowska M, Sedghamiz E, Setzler J, Brieg M, Wenzel W. Sampling of the conformational landscape of small proteins with Monte Carlo methods. Sci Rep 2020; 10:18211. [PMID: 33097750 PMCID: PMC7585447 DOI: 10.1038/s41598-020-75239-7] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2020] [Accepted: 10/12/2020] [Indexed: 12/24/2022] Open
Abstract
Computer simulation provides an increasingly realistic picture of large-scale conformational change of proteins, but investigations remain fundamentally constrained by the femtosecond timestep of molecular dynamics simulations. For this reason, many biologically interesting questions cannot be addressed using accessible state-of-the-art computational resources. Here, we report the development of an all-atom Monte Carlo approach that permits the modelling of the large-scale conformational change of proteins using standard off-the-shelf computational hardware and standard all-atom force fields. We demonstrate extensive thermodynamic characterization of the folding process of the α-helical Trp-cage, the Villin headpiece and the β-sheet WW-domain. We fully characterize the free energy landscape, transition states, energy barriers between different states, and the per-residue stability of individual amino acids over a wide temperature range. We demonstrate that a state-of-the-art intramolecular force field can be combined with an implicit solvent model to obtain a high quality of the folded structures and also discuss limitations that still remain.
Collapse
Affiliation(s)
- Nana Heilmann
- Institute of Nanotechnology (INT), Karlsruhe Institute of Technology (KIT), Hermann-von-Helmholtz-Platz 1, 76344, Eggenstein-Leopoldshafen, Germany
| | - Moritz Wolf
- Institute of Nanotechnology (INT), Karlsruhe Institute of Technology (KIT), Hermann-von-Helmholtz-Platz 1, 76344, Eggenstein-Leopoldshafen, Germany
| | - Mariana Kozlowska
- Institute of Nanotechnology (INT), Karlsruhe Institute of Technology (KIT), Hermann-von-Helmholtz-Platz 1, 76344, Eggenstein-Leopoldshafen, Germany
| | - Elaheh Sedghamiz
- Institute of Nanotechnology (INT), Karlsruhe Institute of Technology (KIT), Hermann-von-Helmholtz-Platz 1, 76344, Eggenstein-Leopoldshafen, Germany
| | - Julia Setzler
- Institute of Nanotechnology (INT), Karlsruhe Institute of Technology (KIT), Hermann-von-Helmholtz-Platz 1, 76344, Eggenstein-Leopoldshafen, Germany
| | - Martin Brieg
- Institute of Nanotechnology (INT), Karlsruhe Institute of Technology (KIT), Hermann-von-Helmholtz-Platz 1, 76344, Eggenstein-Leopoldshafen, Germany
| | - Wolfgang Wenzel
- Institute of Nanotechnology (INT), Karlsruhe Institute of Technology (KIT), Hermann-von-Helmholtz-Platz 1, 76344, Eggenstein-Leopoldshafen, Germany.
| |
Collapse
|
6
|
Tao P, Xiao Y. Using the generalized Born surface area model to fold proteins yields more effective sampling while qualitatively preserving the folding landscape. Phys Rev E 2020; 101:062417. [PMID: 32688556 DOI: 10.1103/physreve.101.062417] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2020] [Accepted: 06/01/2020] [Indexed: 11/07/2022]
Abstract
Protein folding is a long-standing problem and has been widely investigated using molecular dynamics simulations with both explicit and implicit solvents. However, to what extent the folding mechanisms observed in two water models agree remains an open question. In this study, ab initio folding simulations of ten proteins with different topologies are performed in two combinations of force fields and water models (ff14SB+TIP3P and ff14SBonlysc+GB-Neck2). Interestingly, the latter combination not only folds more proteins but also provides a better balance of different secondary structures than the former in the same number of integration time steps. More importantly, the folding pathways found in the two types of simulations are conserved and they may only differ in their weights. Our results suggest that simulations with an implicit solvent may also be suitable for the investigation of the mechanism of protein folding.
Collapse
Affiliation(s)
- Peng Tao
- School of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| | - Yi Xiao
- School of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| |
Collapse
|
7
|
Oide M, Kato T, Oroguchi T, Nakasako M. Energy landscape of domain motion in glutamate dehydrogenase deduced from cryo-electron microscopy. FEBS J 2020; 287:3472-3493. [PMID: 31976609 DOI: 10.1111/febs.15224] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2019] [Revised: 11/09/2019] [Accepted: 01/21/2020] [Indexed: 11/28/2022]
Abstract
Analysis of the conformational changes of protein is important to elucidate the mechanisms of protein motions correlating with their function. Here, we studied the spontaneous domain motion of unliganded glutamate dehydrogenase from Thermococcus profundus using cryo-electron microscopy and proposed a novel method to construct free-energy landscape of protein conformations. Each subunit of the homo-hexameric enzyme comprises nucleotide-binding domain (NAD domain) and hexamer-forming core domain. A large active-site cleft is situated between the two domains and varies from open to close according to the motion of a NAD domain. A three-dimensional map reconstructed from all cryo-electron microscopy images displayed disordered volumes of NAD domains, suggesting that NAD domains in the collected images adopted various conformations in domain motion. Focused classifications on NAD domain of subunits provided several maps of possible conformations in domain motion. To deduce what kinds of conformations appeared in EM images, we developed a novel analysis method that describe the EM maps as a linear combination of representative conformations appearing in a 200-ns molecular dynamics simulation as reference. The analysis enabled us to estimate the appearance frequencies of the representative conformations, which illustrated a free-energy landscape in domain motion. In the open/close domain motion, two free-energy basins hindered the direct transformation from open to closed state. Structure models constructed for representative EM maps in classifications demonstrated the correlation between the energy landscape and conformations in domain motion. Based on the results, the domain motion in glutamate dehydrogenase and the analysis method to visualize conformational changes and free-energy landscape were discussed. DATABASE: The EM maps of the four conformations were deposited to Electron Microscopy Data Bank (EMDB) as accession codes EMD-9845 (open), EMD-9846 (half-open1), EMD-9847 (half-open2), and EMD-9848 (closed), respectively. In addition, the structural models built for the four conformations were deposited to the Protein Data Bank (PDB) as accession codes 6JN9 (open), 6JNA (half-open1), 6JNC (half-open2), and 6JND (closed), respectively.
Collapse
Affiliation(s)
- Mao Oide
- Department of Physics, Faculty of Science and Technology, Keio University, Yokohama, Japan.,RIKEN SPring-8 Center, Sayo-gun, Japan
| | - Takayuki Kato
- Graduate School of Frontier Biosciences, Osaka University, Suita, Japan
| | - Tomotaka Oroguchi
- Department of Physics, Faculty of Science and Technology, Keio University, Yokohama, Japan.,RIKEN SPring-8 Center, Sayo-gun, Japan
| | - Masayoshi Nakasako
- Department of Physics, Faculty of Science and Technology, Keio University, Yokohama, Japan.,RIKEN SPring-8 Center, Sayo-gun, Japan
| |
Collapse
|
8
|
Abstract
Comprehensive data about the composition and structure of cellular components have enabled the construction of quantitative whole-cell models. While kinetic network-type models have been established, it is also becoming possible to build physical, molecular-level models of cellular environments. This review outlines challenges in constructing and simulating such models and discusses near- and long-term opportunities for developing physical whole-cell models that can connect molecular structure with biological function.
Collapse
Affiliation(s)
- Michael Feig
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824, USA;
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Kobe, Hyogo 650-0047, Japan
| | - Yuji Sugita
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, Kobe, Hyogo 650-0047, Japan
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama 351-0198, Japan
| |
Collapse
|
9
|
Wang E, Tao P, Wang J, Xiao Y. A novel folding pathway of the villin headpiece subdomain HP35. Phys Chem Chem Phys 2019; 21:18219-18226. [PMID: 31389931 DOI: 10.1039/c9cp01703h] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The villin headpiece subdomain (HP35) is a fast-folding protein with 35 residues and its folding pathways have been extensively studied experimentally and theoretically but remain controversial. While experiments showed that HP35 might have multiple folding pathways, most theoretical studies only found one major pathway, although a few theoretical studies revealed two. Here we report our results of molecular dynamics simulations of HP35 folding by using the newest AMBER ff14SB force field and show that HP35 has a novel folding pathway in addition to the two pathways shown previously. We also study the mechanism of determining the folding pathways and found that the dynamics of Helix2 may play a special role in the folding of HP35. Our results may be helpful to understand the folding mechanism of HP35 further.
Collapse
Affiliation(s)
- Ercheng Wang
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | | | | | | |
Collapse
|
10
|
Harada R, Shigeta Y. Hybrid Cascade-Type Molecular Dynamics with a Markov State Model for Efficient Free Energy Calculations. J Chem Theory Comput 2018; 15:680-687. [PMID: 30468705 DOI: 10.1021/acs.jctc.8b00802] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A protocol for calculating free energy landscapes (FELs) is proposed based on a combination of two cascade-type molecular dynamics (MD) methods, parallel cascade selection MD (PaCS-MD) and outlier flooding method (OFLOOD), with the help of a Markov state model (MSM). The former rapidly generates approximated transition paths directly connecting reactants with products, and the latter complementary resamples marginal conformational subspaces. Trajectories obtained by them give reliable microstates in MSM providing accurate FEL with low computational costs. As a demonstration, the present method was applied to a miniprotein (Chignolin and Trp-cage) in explicit water and successfully elucidated multiple folding paths on their free energy landscapes. Our method could be applicable to a wide variety of biological systems to estimate their free energy profiles.
Collapse
Affiliation(s)
- Ryuhei Harada
- Center for Computational Sciences , University of Tsukuba , 1-1-1 Tennodai , Tsukuba , Ibaraki 305-8577 , Japan
| | - Yasuteru Shigeta
- Center for Computational Sciences , University of Tsukuba , 1-1-1 Tennodai , Tsukuba , Ibaraki 305-8577 , Japan
| |
Collapse
|
11
|
Harada R, Aida H, Shigeta Y. The Formation of Hydrophobic Core Regulates the Protein Folding of Villin Elucidated with Parallel Cascade Selection Molecular Dynamics. CHEM LETT 2018. [DOI: 10.1246/cl.180596] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Ryuhei Harada
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-0821, Japan
| | - Hayato Aida
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-0821, Japan
| | - Yasuteru Shigeta
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-0821, Japan
| |
Collapse
|
12
|
Harada R. Simple, yet Efficient Conformational Sampling Methods for Reproducing/Predicting Biologically Rare Events of Proteins. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2018. [DOI: 10.1246/bcsj.20180170] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Affiliation(s)
- Ryuhei Harada
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8571, Japan
| |
Collapse
|
13
|
Hayashi T, Yasuda S, Škrbić T, Giacometti A, Kinoshita M. Unraveling protein folding mechanism by analyzing the hierarchy of models with increasing level of detail. J Chem Phys 2017; 147:125102. [DOI: 10.1063/1.4999376] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Affiliation(s)
- Tomohiko Hayashi
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611-0011, Japan
| | - Satoshi Yasuda
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611-0011, Japan
- Graduate School of Science, Chiba University, 1-33 Yayoi-cho, Inage, Chiba 263-8522, Japan
- Molecular Chirality Research Center, Chiba University, 1-33 Yayoi-cho, Inage, Chiba 263-8522, Japan
| | - Tatjana Škrbić
- Dipartimento di Scienze Molecolari e Nanosistemi, Università Ca’ Foscari Venezia, Edificio Alfa Campus Scientifico, Via Torino 155, Venezia-Mestre I-3010, Italy
| | - Achille Giacometti
- Dipartimento di Scienze Molecolari e Nanosistemi, Università Ca’ Foscari Venezia, Edificio Alfa Campus Scientifico, Via Torino 155, Venezia-Mestre I-3010, Italy
| | - Masahiro Kinoshita
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611-0011, Japan
| |
Collapse
|
14
|
Harada R, Takano Y, Shigeta Y. Common folding processes of mini-proteins: Partial formations of secondary structures initiate the immediate protein folding. J Comput Chem 2017; 38:790-797. [DOI: 10.1002/jcc.24748] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2016] [Revised: 12/20/2016] [Accepted: 01/12/2017] [Indexed: 11/10/2022]
Affiliation(s)
- Ryuhei Harada
- Center for Computational Sciences, University of Tsukuba; 1-1-1 Tennodai Tsukuba Ibaraki 305-8577 Japan
| | - Yu Takano
- Department of Biomedical Information Sciences; Graduate School of Information Sciences, Hiroshima City University; 3-4-1 Ozuka-Higashi, Asa-Minami-Ku Hiroshima 731-3194 Japan
| | - Yasuteru Shigeta
- Center for Computational Sciences, University of Tsukuba; 1-1-1 Tennodai Tsukuba Ibaraki 305-8577 Japan
| |
Collapse
|
15
|
Harada R, Nakamura T, Shigeta Y. A Fast Convergent Simulated Annealing Algorithm for Protein-Folding: Simulated Annealing Outlier FLOODing (SA-OFLOOD) Method. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2016. [DOI: 10.1246/bcsj.20160244] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
|
16
|
Harada R, Takano Y, Baba T, Shigeta Y. Simple, yet powerful methodologies for conformational sampling of proteins. Phys Chem Chem Phys 2016; 17:6155-73. [PMID: 25659594 DOI: 10.1039/c4cp05262e] [Citation(s) in RCA: 46] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Several biological functions, such as molecular recognition, enzyme catalysis, signal transduction, allosteric regulation, and protein folding, are strongly related to conformational transitions of proteins. These conformational transitions are generally induced as slow dynamics upon collective motions, including biologically relevant large-amplitude fluctuations of proteins. Although molecular dynamics (MD) simulation has become a powerful tool for extracting conformational transitions of proteins, it might still be difficult to reach time scales of the biological functions because the accessible time scales of MD simulations are far from biological time scales, even if straightforward conventional MD (CMD) simulations using massively parallel computers are employed. Thus, it is desirable to develop efficient methods to achieve canonical ensembles with low computational costs. From this perspective, we review several enhanced conformational sampling techniques of biomolecules developed by us. In our methods, multiple independent short-time MD simulations are employed instead of single straightforward long-time CMD simulations. Our basic strategy is as follows: (i) selection of initial seeds (initial structures) for the conformational sampling in restarting MD simulations. Here, the seeds should be selected as candidates with high potential to transit. (ii) Resampling from the selected seeds by initializing velocities in restarting short-time MD simulations. A cycle of these simple protocols might drastically promote the conformational transitions of biomolecules. (iii) Once reactive trajectories extracted from the cycles of short-time MD simulations are obtained, a free energy profile is evaluated by means of umbrella sampling (US) techniques with the weighted histogram analysis method (WHAM) as a post-processing technique. For the selection of the initial seeds, we proposed four different choices: (1) Parallel CaScade molecular dynamics (PaCS-MD), (2) Fluctuation Flooding Method (FFM), (3) Outlier FLOODing (OFLOOD) method, and (4) TaBoo SeArch (TBSA) method. We demonstrate applications of our methods to several biological systems, such as domain motions of proteins with large-amplitude fluctuations, conformational transitions upon ligand binding, and protein folding/refolding to native structures of proteins. Finally, we show the conformational sampling efficiencies of our methods compared with those by CMD simulations and other previously developed enhanced conformational sampling methods.
Collapse
Affiliation(s)
- Ryuhei Harada
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8571, Japan.
| | | | | | | |
Collapse
|
17
|
Harada R, Kitao A. Nontargeted Parallel Cascade Selection Molecular Dynamics for Enhancing the Conformational Sampling of Proteins. J Chem Theory Comput 2015; 11:5493-502. [PMID: 26574337 DOI: 10.1021/acs.jctc.5b00723] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Nontargeted parallel cascade selection molecular dynamics (nt-PaCS-MD) is proposed as an efficient conformational sampling method to enhance the conformational transitions of proteins, which is an extension of the original targeted PaCS-MD (t-PaCS-MD). The original PaCS-MD comprises cycles of (i) selection of initial structures for multiple independent MD simulations toward a predetermined target and (ii) conformational sampling by the independent MDs. In nt-PaCS-MD, structures that significantly deviate from an average are regarded as candidates that have high potential to address other metastable states and are chosen as the initial structures in the selection. To select significantly deviated structures, we examine the root-mean-square deviation (RMSD) of snapshots generated from the average structure based on Gram-Schmidt orthogonalization. nt-PaCS-MD was applied to the folding of the mini-protein chignolin in implicit solvent and to the open-closed conformational transitions of T4 lysozyme (T4L) and glutamine binding protein (QBP) in explicit solvent. We show that nt-PaCS-MD can reach chignolin's native state and can also cause the open-closed transition of T4L and QBP on a nanosecond time scale, which are very efficient in terms of conformational sampling and comparable to that with t-PaCS-MD.
Collapse
Affiliation(s)
| | - Akio Kitao
- Institute of Molecular and Cellular Bioscience, The University of Tokyo , 1-1-1 Yayoi, Bunkyo-ku, Tokyo 113-0032, Japan
| |
Collapse
|
18
|
Miao Y, Feixas F, Eun C, McCammon JA. Accelerated molecular dynamics simulations of protein folding. J Comput Chem 2015; 36:1536-49. [PMID: 26096263 DOI: 10.1002/jcc.23964] [Citation(s) in RCA: 110] [Impact Index Per Article: 12.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2015] [Revised: 05/11/2015] [Accepted: 05/19/2015] [Indexed: 02/02/2023]
Abstract
Folding of four fast-folding proteins, including chignolin, Trp-cage, villin headpiece and WW domain, was simulated via accelerated molecular dynamics (aMD). In comparison with hundred-of-microsecond timescale conventional molecular dynamics (cMD) simulations performed on the Anton supercomputer, aMD captured complete folding of the four proteins in significantly shorter simulation time. The folded protein conformations were found within 0.2-2.1 Å of the native NMR or X-ray crystal structures. Free energy profiles calculated through improved reweighting of the aMD simulations using cumulant expansion to the second-order are in good agreement with those obtained from cMD simulations. This allows us to identify distinct conformational states (e.g., unfolded and intermediate) other than the native structure and the protein folding energy barriers. Detailed analysis of protein secondary structures and local key residue interactions provided important insights into the protein folding pathways. Furthermore, the selections of force fields and aMD simulation parameters are discussed in detail. Our work shows usefulness and accuracy of aMD in studying protein folding, providing basic references in using aMD in future protein-folding studies.
Collapse
Affiliation(s)
- Yinglong Miao
- Howard Hughes Medical Institute, University of California at San Diego, La Jolla, California
| | - Ferran Feixas
- Department of Chemistry and Biochemistry, University of California at San Diego, La Jolla, California.,Department of Pharmacology, University of California at San Diego, La Jolla, California
| | - Changsun Eun
- Howard Hughes Medical Institute, University of California at San Diego, La Jolla, California
| | - J Andrew McCammon
- Howard Hughes Medical Institute, University of California at San Diego, La Jolla, California.,Department of Chemistry and Biochemistry, University of California at San Diego, La Jolla, California.,Department of Pharmacology, University of California at San Diego, La Jolla, California
| |
Collapse
|
19
|
Efficient conformational sampling of proteins based on a multi-dimensional TaBoo SeArch algorithm: An application to folding of chignolin in explicit solvent. Chem Phys Lett 2015. [DOI: 10.1016/j.cplett.2015.04.039] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
20
|
Harada R, Kitao A. Parallel Cascade Selection Molecular Dynamics (PaCS-MD) to generate conformational transition pathway. J Chem Phys 2014; 139:035103. [PMID: 23883057 DOI: 10.1063/1.4813023] [Citation(s) in RCA: 91] [Impact Index Per Article: 9.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Parallel Cascade Selection Molecular Dynamics (PaCS-MD) is proposed as a molecular simulation method to generate conformational transition pathway under the condition that a set of "reactant" and "product" structures is known a priori. In PaCS-MD, the cycle of short multiple independent molecular dynamics simulations and selection of the structures close to the product structure for the next cycle are repeated until the simulated structures move sufficiently close to the product. Folding of 10-residue mini-protein chignolin from the extended to native structures and open-close conformational transition of T4 lysozyme were investigated by PaCS-MD. In both cases, tens of cycles of 100-ps MD were sufficient to reach the product structures, indicating the efficient generation of conformational transition pathway in PaCS-MD with a series of conventional MD without additional external biases. Using the snapshots along the pathway as the initial coordinates, free energy landscapes were calculated by the combination with multiple independent umbrella samplings to statistically elucidate the conformational transition pathways.
Collapse
Affiliation(s)
- Ryuhei Harada
- Department of Physics, Graduate School of Science, The University of Tokyo, Tokyo, 7-3-1, Hongo, Japan
| | | |
Collapse
|
21
|
Yamamori Y, Kitao A. MuSTAR MD: Multi-scale sampling using temperature accelerated and replica exchange molecular dynamics. J Chem Phys 2013; 139:145105. [DOI: 10.1063/1.4823743] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
|
22
|
Harada R, Tochio N, Kigawa T, Sugita Y, Feig M. Reduced native state stability in crowded cellular environment due to protein-protein interactions. J Am Chem Soc 2013; 135:3696-701. [PMID: 23402619 DOI: 10.1021/ja3126992] [Citation(s) in RCA: 130] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The effect of cellular crowding environments on protein structure and stability is a key issue in molecular and cellular biology. The classical view of crowding emphasizes the volume exclusion effect that generally favors compact, native states. Here, results from molecular dynamics simulations and NMR experiments show that protein crowders may destabilize native states via protein-protein interactions. In the model system considered here, mixtures of villin head piece and protein G at high concentrations, villin structures become increasingly destabilized upon increasing crowder concentrations. The denatured states observed in the simulation involve partial unfolding as well as more subtle conformational shifts. The unfolded states remain overall compact and only partially overlap with unfolded ensembles at high temperature and in the presence of urea. NMR measurements on the same systems confirm structural changes upon crowding based on changes of chemical shifts relative to dilute conditions. An analysis of protein-protein interactions and energetic aspects suggests the importance of enthalpic and solvation contributions to the crowding free energies that challenge an entropic-centered view of crowding effects.
Collapse
Affiliation(s)
- Ryuhei Harada
- RIKEN Advanced Institute for Computational Science, 7-1-26 minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047 Japan
| | | | | | | | | |
Collapse
|