1
|
da Silva FB, Martins de Oliveira V, de Oliveira Junior AB, Contessoto VDG, Leite VBP. Probing the Energy Landscape of Spectrin R15 and R16 and the Effects of Non-native Interactions. J Phys Chem B 2023; 127:1291-1300. [PMID: 36723393 DOI: 10.1021/acs.jpcb.2c06178] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
Understanding the details of a protein folding mechanism can be a challenging and complex task. One system with an interesting folding behavior is the α-spectrin domain, where the R15 folds three-orders of magnitude faster than its homologues R16 and R17, despite having similar structures. The molecular origins that explain these folding rate differences remain unclear, but our previous work revealed that a combined effect produced by non-native interactions could be a reasonable cause for these differences. In this study, we explore further the folding process by identifying the molecular paths, metastable states, and the collective motions that lead these unfolded proteins to their native state conformation. Our results uncovered the differences between the folding pathways for the wild-type R15 and R16 and an R16 mutant. The metastable ensembles that speed down the folding were identified using an energy landscape visualization method (ELViM). These ensembles correspond to similar experimentally reported configurations. Our observations indicate that the non-native interactions are also associated with secondary structure misdocking. This computational methodology can be used as a fast, straightforward protocol for shedding light on systems with unclear folding or conformational traps.
Collapse
Affiliation(s)
- Fernando Bruno da Silva
- Department of Physics, São Paulo State University (UNESP), Institute of Biosciences, Humanities and Exact Sciences, São José do Rio Preto, São Paulo15054-000, Brazil
| | - Vinícius Martins de Oliveira
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, Maryland21201, United States
| | | | | | - Vitor B P Leite
- Department of Physics, São Paulo State University (UNESP), Institute of Biosciences, Humanities and Exact Sciences, São José do Rio Preto, São Paulo15054-000, Brazil
| |
Collapse
|
2
|
Contessoto VG, Ferreira PHB, Chahine J, Leite VBP, Oliveira RJ. Small Neutral Crowding Solute Effects on Protein Folding Thermodynamic Stability and Kinetics. J Phys Chem B 2021; 125:11673-11686. [PMID: 34644091 DOI: 10.1021/acs.jpcb.1c07663] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Molecular crowding is a ubiquitous phenomenon in biological systems, with significant consequences on protein folding and stability. Small compounds, such as the osmolyte trimethylamine N-oxide (TMAO), can also present similar effects. To analyze the effects arising from these solute-like molecules, we performed a series of crowded coarse-grained folding simulations. Two well-known proteins were chosen, CI2 and SH3, modeled by the alpha-carbon-structure-based model. In the simulations, the crowding molecules were represented by low-sized neutral atom beads in different concentrations. The results show that a low level of the volume fraction occupied by neutral agents can change protein stability and folding kinetics for the two systems. However, the kinetics were shown to be unaffected in their respective folding temperatures. The faster kinetics correlates with changes in the folding route for systems at the same temperature, where non-native contacts were shown to be relevant. The transition states of the two systems with and without crowders are similar. In the case of SH3, there are differences in the structuring of two strands, which may be associated with the increase in its folding rate, in addition to the destabilization of the denatured ensemble. The present study also detected a crossover in the thermodynamic stability behavior, previously observed experimentally and theoretically. As the temperature increases, crowders change from destabilizing to stabilizing agents.
Collapse
Affiliation(s)
- Vinícius G Contessoto
- Department of Physics, Institute of Biosciences, Letters and Exact Sciences, São Paulo State University, São José do Rio Preto 15054-000, Brazil
| | - Paulo H B Ferreira
- Laboratório de Biofísica Teórica, Departamento de Física, Instituto de Ciências Exatas, Naturais e Educação, Universidade Federal do Triângulo Mineiro, Uberaba 38064-200, Brazil
| | - Jorge Chahine
- Department of Physics, Institute of Biosciences, Letters and Exact Sciences, São Paulo State University, São José do Rio Preto 15054-000, Brazil
| | - Vitor B P Leite
- Department of Physics, Institute of Biosciences, Letters and Exact Sciences, São Paulo State University, São José do Rio Preto 15054-000, Brazil
| | - Ronaldo J Oliveira
- Laboratório de Biofísica Teórica, Departamento de Física, Instituto de Ciências Exatas, Naturais e Educação, Universidade Federal do Triângulo Mineiro, Uberaba 38064-200, Brazil
| |
Collapse
|
3
|
Wang S, Gong W, Deng X, Liu Y, Li C. Exploring the dynamics of RNA molecules with multiscale Gaussian network model. Chem Phys 2020. [DOI: 10.1016/j.chemphys.2020.110820] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
|
4
|
Hattori LT, Pinheiro BA, Frigori RB, Benítez CMV, Lopes HS. PathMolD-AB: Spatiotemporal pathways of protein folding using parallel molecular dynamics with a coarse-grained model. Comput Biol Chem 2020; 87:107301. [PMID: 32554177 DOI: 10.1016/j.compbiolchem.2020.107301] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2020] [Revised: 05/25/2020] [Accepted: 05/28/2020] [Indexed: 10/24/2022]
Abstract
Solving the protein folding problem (PFP) is one of the grand challenges still open in computational biophysics. Globular proteins are believed to evolve from initial configurations through folding pathways connecting several thermodynamically accessible states in a free energy landscape until reaching its minimum, inhabited by the stable native structures. Despite its huge computational burden, molecular dynamics (MD) is the leading approach in the PFP studies by preserving the Newtonian temporal evolution in the canonical ensemble. Non-trivial improvements are provided by highly parallel implementations of MD in cost-effective GPUs, concomitant to multiscale descriptions of proteins by coarse-grained minimalist models. In this vein, we present the PathMolD-AB framework, a comprehensive software package for massively parallel MD simulations using the canonical ensemble, structural analysis, and visualization of the folding pathways using the minimalist AB-model. It has, also, a tool to compare the results with proteins re-scaled from the PDB. We simulate and analyze, as case studies, the folding of four proteins: 13FIBO, 2GB1, 1PLC and 5ANZ, with 13, 55, 99 and 223 amino acids, respectively. The datasets generated from simulations correspond to the MD evolution of 3500 folding pathways, encompassing 35×106 states, which contains the spatial amino acid positions, the protein free energies and radii of gyration at each time step. Results indicate that the speedup of our approach grows logarithmically with the protein length and, therefore, it is suited for most of the proteins in the PDB. The predicted structures simulated by PathMolD-AB were similar to the re-scaled biological structures, indicating that it is promising for the study of the PFP study.
Collapse
Affiliation(s)
- Leandro Takeshi Hattori
- Bioinformatics and Computational Intelligence Laboratory (LABIC), Federal University of Technology Paraná (UTFPR), Av. 7 de Setembro, 3165, 80230-901 Curitiba, PR, Brazil.
| | - Bruna Araujo Pinheiro
- Bioinformatics and Computational Intelligence Laboratory (LABIC), Federal University of Technology Paraná (UTFPR), Av. 7 de Setembro, 3165, 80230-901 Curitiba, PR, Brazil.
| | - Rafael Bertolini Frigori
- Bioinformatics and Computational Intelligence Laboratory (LABIC), Federal University of Technology Paraná (UTFPR), Av. 7 de Setembro, 3165, 80230-901 Curitiba, PR, Brazil.
| | - César Manuel Vargas Benítez
- Bioinformatics and Computational Intelligence Laboratory (LABIC), Federal University of Technology Paraná (UTFPR), Av. 7 de Setembro, 3165, 80230-901 Curitiba, PR, Brazil
| | - Heitor Silvério Lopes
- Bioinformatics and Computational Intelligence Laboratory (LABIC), Federal University of Technology Paraná (UTFPR), Av. 7 de Setembro, 3165, 80230-901 Curitiba, PR, Brazil.
| |
Collapse
|
5
|
Freitas FC, Lima AN, Contessoto VDG, Whitford PC, Oliveira RJD. Drift-diffusion (DrDiff) framework determines kinetics and thermodynamics of two-state folding trajectory and tunes diffusion models. J Chem Phys 2019; 151:114106. [DOI: 10.1063/1.5113499] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Affiliation(s)
- Frederico Campos Freitas
- Laboratório de Biofísica Teórica, Departamento de Física, Instituto de Ciências Exatas, Naturais e Educação, Universidade Federal do Triângulo Mineiro, Uberaba, MG, Brazil
| | - Angelica Nakagawa Lima
- Laboratório de Biofísica Teórica, Departamento de Física, Instituto de Ciências Exatas, Naturais e Educação, Universidade Federal do Triângulo Mineiro, Uberaba, MG, Brazil
- Laboratório de Biologia Computacional e Bioinformática, Universidade Federal do ABC, Santo André, SP, Brazil
| | - Vinícius de Godoi Contessoto
- Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, USA
- Departamento de Física, Universidade Estadual Paulista, São José do Rio Preto, SP, Brazil
- Brazilian Biorenewables National Laboratory - LNBR, Brazilian Center for Research in Energy and Materials - CNPEM, Campinas, SP, Brazil
| | - Paul C. Whitford
- Department of Physics, Northeastern University, Boston, Massachusetts 02115, USA
| | - Ronaldo Junio de Oliveira
- Laboratório de Biofísica Teórica, Departamento de Física, Instituto de Ciências Exatas, Naturais e Educação, Universidade Federal do Triângulo Mineiro, Uberaba, MG, Brazil
| |
Collapse
|
6
|
Bruno da Silva F, Contessoto VG, de Oliveira VM, Clarke J, Leite VBP. Non-Native Cooperative Interactions Modulate Protein Folding Rates. J Phys Chem B 2018; 122:10817-10824. [DOI: 10.1021/acs.jpcb.8b08990] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Fernando Bruno da Silva
- Department of Physics, Institute of Biosciences, Humanities and Exact Sciences, São Paulo State University (UNESP), São José do Rio Preto - São Paulo 15054-000, Brazil
| | - Vinícius G. Contessoto
- Brazilian Bioethanol Science and Technology Laboratory - CTBE, Campinas - São Paulo 13083-100, Brazil
| | - Vinícius M. de Oliveira
- Department of Physics, Institute of Biosciences, Humanities and Exact Sciences, São Paulo State University (UNESP), São José do Rio Preto - São Paulo 15054-000, Brazil
| | - Jane Clarke
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Vitor B. P. Leite
- Department of Physics, Institute of Biosciences, Humanities and Exact Sciences, São Paulo State University (UNESP), São José do Rio Preto - São Paulo 15054-000, Brazil
| |
Collapse
|
7
|
Martins de Oliveira V, Godoi Contessoto VD, Bruno da Silva F, Zago Caetano DL, Jurado de Carvalho S, Pereira Leite VB. Effects of pH and Salt Concentration on Stability of a Protein G Variant Using Coarse-Grained Models. Biophys J 2018; 114:65-75. [PMID: 29320697 PMCID: PMC5984902 DOI: 10.1016/j.bpj.2017.11.012] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2017] [Revised: 10/20/2017] [Accepted: 11/13/2017] [Indexed: 11/18/2022] Open
Abstract
The importance of charge-charge interactions in the thermal stability of proteins is widely known. pH and ionic strength play a crucial role in these electrostatic interactions, as well as in the arrangement of ionizable residues in each protein-folding stage. In this study, two coarse-grained models were used to evaluate the effect of pH and salt concentration on the thermal stability of a protein G variant (1PGB-QDD), which was chosen due to the quantity of experimental data exploring these effects on its stability. One of these coarse-grained models, the TKSA, calculates the electrostatic free energy of the protein in the native state via the Tanford-Kirkwood approach for each residue. The other one, CpHMD-SBM, uses a Coulomb screening potential in addition to the structure-based model Cα. Both models simulate the system in constant pH. The comparison between the experimental stability analysis and the computational results obtained by these simple models showed a good agreement. Through the TKSA method, the role of each charged residue in the protein's thermal stability was inferred. Using CpHMD-SBM, it was possible to evaluate salt and pH effects throughout the folding process. Finally, the computational pKa values were calculated by both methods and presented a good level of agreement with the experiments. This study provides, to our knowledge, new information and a comprehensive description of the electrostatic contribution to protein G stability.
Collapse
Affiliation(s)
- Vinícius Martins de Oliveira
- São Paulo State University (UNESP), Institute of Biosciences, Humanities and Exact Sciences (Ibilce), Campus São José do Rio Preto, Brazil
| | - Vinícius de Godoi Contessoto
- São Paulo State University (UNESP), Institute of Biosciences, Humanities and Exact Sciences (Ibilce), Campus São José do Rio Preto, Brazil; Brazilian Bioethanol Science and Technology Laboratory- (CTBE), Campinas, Brazil
| | - Fernando Bruno da Silva
- São Paulo State University (UNESP), Institute of Biosciences, Humanities and Exact Sciences (Ibilce), Campus São José do Rio Preto, Brazil
| | - Daniel Lucas Zago Caetano
- São Paulo State University (UNESP), Institute of Biosciences, Humanities and Exact Sciences (Ibilce), Campus São José do Rio Preto, Brazil
| | - Sidney Jurado de Carvalho
- São Paulo State University (UNESP), Institute of Biosciences, Humanities and Exact Sciences (Ibilce), Campus São José do Rio Preto, Brazil
| | - Vitor Barbanti Pereira Leite
- São Paulo State University (UNESP), Institute of Biosciences, Humanities and Exact Sciences (Ibilce), Campus São José do Rio Preto, Brazil.
| |
Collapse
|
8
|
Mouro PR, de Godoi Contessoto V, Chahine J, Junio de Oliveira R, Pereira Leite VB. Quantifying Nonnative Interactions in the Protein-Folding Free-Energy Landscape. Biophys J 2017; 111:287-293. [PMID: 27463131 DOI: 10.1016/j.bpj.2016.05.041] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2015] [Revised: 05/10/2016] [Accepted: 05/17/2016] [Indexed: 11/27/2022] Open
Abstract
Protein folding is a central problem in biological physics. Energetic roughness is an important aspect that controls protein-folding stability and kinetics. The roughness is associated with conflicting interactions in the protein and is also known as frustration. Recent studies indicate that an addition of a small amount of energetic frustration may enhance folding speed for certain proteins. In this study, we have investigated the conditions under which frustration increases the folding rate. We used a Cα structure-based model to simulate a group of proteins. We found that the free-energy barrier at the transition state (ΔF) correlates with nonnative-contact variation (ΔA), and the simulated proteins are clustered according to their fold motifs. These findings are corroborated by the Clementi-Plotkin analytical model. As a consequence, the optimum frustration regime for protein folding can be predicted analytically.
Collapse
Affiliation(s)
- Paulo Ricardo Mouro
- Departamento de Física, Instituto de Biociências, Letras e Ciências Exatas, Universidade Estadual Paulista, São José do Rio Preto, São Paulo, Brazil
| | - Vinícius de Godoi Contessoto
- Departamento de Física, Instituto de Biociências, Letras e Ciências Exatas, Universidade Estadual Paulista, São José do Rio Preto, São Paulo, Brazil
| | - Jorge Chahine
- Departamento de Física, Instituto de Biociências, Letras e Ciências Exatas, Universidade Estadual Paulista, São José do Rio Preto, São Paulo, Brazil
| | - Ronaldo Junio de Oliveira
- Laboratório de Biofísica Teórica, Departamento de Física, Instituto de Ciências Exatas, Naturais e Educação, Universidade Federal do Triângulo Mineiro, Uberaba, Minas Gerais, Brazil
| | - Vitor Barbanti Pereira Leite
- Departamento de Física, Instituto de Biociências, Letras e Ciências Exatas, Universidade Estadual Paulista, São José do Rio Preto, São Paulo, Brazil.
| |
Collapse
|
9
|
Contessoto VG, de Oliveira VM, de Carvalho SJ, Oliveira LC, Leite VBP. NTL9 Folding at Constant pH: The Importance of Electrostatic Interaction and pH Dependence. J Chem Theory Comput 2016; 12:3270-7. [PMID: 27327651 DOI: 10.1021/acs.jctc.6b00399] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The folding process of the N-terminal domain of ribosomal protein L9 (NTL9) was investigated at constant-pH computer simulations. Evaluation of the role of electrostatic interaction during folding was carried out by including a Debye-Hückel potential into a Cα structure-based model (SBM). In this study, the charges of the ionizable residues and the electrostatic potential are susceptible to the solution conditions, such as pH and ionic strength, as well as to the presence of charged groups. Simulations were performed under different pHs, and the results were validated by comparing them with experimental values of pKa and with denaturation experiment data. Also, the free energy profiles, Φ-values, and folding routes were calculated for each condition. It was shown how charges vary along the folding under different pH, which is subject to different scenarios. This study reveals how simplified models can capture essential physical features, reproducing experimental results, and presenting the role of electrostatic interactions before, during, and after the transition state.
Collapse
Affiliation(s)
- Vinícius G Contessoto
- Departamento de Física, Instituto de Biociências, Letras e Ciências Exatas, Universidade Estadual Paulista (UNESP) , São José do Rio Preto, São Paulo 15054-000, Brazil
| | - Vinícius M de Oliveira
- Departamento de Física, Instituto de Biociências, Letras e Ciências Exatas, Universidade Estadual Paulista (UNESP) , São José do Rio Preto, São Paulo 15054-000, Brazil
| | - Sidney J de Carvalho
- Departamento de Física, Instituto de Biociências, Letras e Ciências Exatas, Universidade Estadual Paulista (UNESP) , São José do Rio Preto, São Paulo 15054-000, Brazil
| | - Leandro C Oliveira
- Departamento de Física, Instituto de Biociências, Letras e Ciências Exatas, Universidade Estadual Paulista (UNESP) , São José do Rio Preto, São Paulo 15054-000, Brazil
| | - Vitor B P Leite
- Departamento de Física, Instituto de Biociências, Letras e Ciências Exatas, Universidade Estadual Paulista (UNESP) , São José do Rio Preto, São Paulo 15054-000, Brazil
| |
Collapse
|