1
|
Oliveira RJD. Coordinate-Dependent Drift-Diffusion Reveals the Kinetic Intermediate Traps of Top7-Based Proteins. J Phys Chem B 2022; 126:10854-10869. [PMID: 36519977 DOI: 10.1021/acs.jpcb.2c07031] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
The computer-designed Top7 served as a scaffold to produce immunoreactive proteins by grafting of the 2F5 HIV-1 antibody epitope (Top7-2F5) followed by biotinylation (Top7-2F5-biotin). The resulting nonimmunoglobulin affinity proteins were effective in inducing and detecting the HIV-1 antibody. However, the grafted Top7-2F5 design led to protein aggregation, as opposed to the soluble biotinylated Top7-2F5-biotin. The structure-based model predicted that the thermodynamic cooperativity of Top7 increases after grafting and biotin-labeling, reducing their intermediate state populations. In this work, the folding kinetic traps that might contribute to the aggregation propensity are investigated by the diffusion theory. Since the engineered proteins have similar sequence and structural homology, they served as protein models to study the kinetic intermediate traps that were uncovered by characterizing the position-dependent drift-velocity (v(Q)) and the diffusion (D(Q)) coefficients. These coordinate-dependent coefficients were taken into account to obtain the folding and transition path times over the free energy transition states containing the intermediate kinetic traps. This analysis may be useful to predict the aggregated kinetic traps of scaffold-epitope proteins that might compose novel diagnostic and therapeutic platforms.
Collapse
Affiliation(s)
- 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, MG38064-200, Brazil
| |
Collapse
|
2
|
Sarkar A, Gasic AG, Cheung MS, Morrison G. Effects of Protein Crowders and Charge on the Folding of Superoxide Dismutase 1 Variants: A Computational Study. J Phys Chem B 2022; 126:4458-4471. [PMID: 35686856 DOI: 10.1021/acs.jpcb.2c00819] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The neurodegenerative disease amyotrophic lateral sclerosis (ALS) is associated with the misfolding and aggregation of the metalloenzyme protein superoxide dismutase 1 (SOD1) via mutations that destabilize the monomer-dimer interface. In a cellular environment, crowding and electrostatic screening play essential roles in the folding and aggregation of the SOD1 monomers. Despite numerous studies on the effects of mutations on SOD1 folding, a clear understanding of the interplay between crowding, folding, and aggregation in vivo remains lacking. Using a structure-based minimal model for molecular dynamics simulations, we investigate the role of self-crowding and charge on the folding stability of SOD1 and the G41D mutant where experimentalists were intrigued by an alteration of the folding mechanism by a single point mutation from glycine to charged aspartic acid. We show that unfolded SOD1 configurations are significantly affected by charge and crowding, a finding that would be extremely costly to achieve with all-atom simulations, while the native state is not significantly altered. The mutation at residue 41 alters the interactions between proteins in the unfolded states instead of those within a protein. This paper suggests electrostatics may play an important role in the folding pathway of SOD1 and modifying the charge via mutation and ion concentration may change the dominant interactions between proteins, with potential impacts for aggregation of the mutants. This work provides a plausible reason for the alteration of the unfolded states to address why the mutant G41D causes the changes to the folding mechanism of SOD1 that have intrigued experimentalists.
Collapse
Affiliation(s)
- Atrayee Sarkar
- Department of Physics, University of Houston, Houston, Texas 77204, United States.,Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, United States
| | - Andrei G Gasic
- Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, United States
| | - Margaret S Cheung
- Department of Physics, University of Houston, Houston, Texas 77204, United States.,Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, United States.,Pacific Northwest National Laboratory, Seattle Research Center, Seattle, Washington 98109, United States
| | - Greg Morrison
- Department of Physics, University of Houston, Houston, Texas 77204, United States.,Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, United States
| |
Collapse
|
3
|
Freitas FC, Maldonado M, Oliveira Junior AB, Onuchic JN, Oliveira RJD. Biotin-painted proteins have thermodynamic stability switched by kinetic folding routes. J Chem Phys 2022; 156:195101. [PMID: 35597640 DOI: 10.1063/5.0083875] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
Abstract
Biotin-labeled proteins are widely used as tools to study protein-protein interactions and proximity in living cells. Proteomic methods broadly employ proximity-labeling technologies based on protein biotinylation in order to investigate the transient encounters of biomolecules in subcellular compartments. Biotinylation is a post-translation modification in which the biotin molecule is attached to lysine or tyrosine residues. So far, biotin-based technologies proved to be effective instruments as affinity and proximity tags. However, the influence of biotinylation on aspects such as folding, binding, mobility, thermodynamic stability, and kinetics needs to be investigated. Here, we selected two proteins [biotin carboxyl carrier protein (BCCP) and FKBP3] to test the influence of biotinylation on thermodynamic and kinetic properties. Apo (without biotin) and holo (biotinylated) protein structures were used separately to generate all-atom structure-based model simulations in a wide range of temperatures. Holo BCCP contains one biotinylation site, and FKBP3 was modeled with up to 23 biotinylated lysines. The two proteins had their estimated thermodynamic stability changed by altering their energy landscape. In all cases, after comparison between the apo and holo simulations, differences were observed on the free-energy profiles and folding routes. Energetic barriers were altered with the density of states clearly showing changes in the transition state. This study suggests that analysis of large-scale datasets of biotinylation-based proximity experiments might consider possible alterations in thermostability and folding mechanisms imposed by the attached biotins.
Collapse
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 38064-200, Brazil
| | - Michelli Maldonado
- Departamento de Matemática, Instituto de Ciências Exatas, Naturais e Educação, Universidade Federal do Triângulo Mineiro, Uberaba, MG 38064-200, Brazil
| | - Antonio Bento Oliveira Junior
- Center for Theoretical Biological Physics, Rice University, BioScience Research Collaborative, 6566 Main St., Houston, Texas 77030, USA
| | - José Nelson Onuchic
- Center for Theoretical Biological Physics, Rice University, BioScience Research Collaborative, 6566 Main St., Houston, Texas 77030, 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 38064-200, Brazil
| |
Collapse
|
4
|
Beltrán HI, Alas-Guardado SJ, González-Pérez PP. Improving coarse-grained models of protein folding through weighting of polar-polar/hydrophobic-hydrophobic interactions into crowded spaces. J Mol Model 2022; 28:87. [PMID: 35262807 DOI: 10.1007/s00894-022-05071-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2021] [Accepted: 02/26/2022] [Indexed: 10/18/2022]
Abstract
Herein were tested 7 hydrophobic-polar sequences in two types of 2D-square space lattices, homogeneous and correlated, the latter simulating molecular crowding included as a geometric boundary restriction. Optimization of 2D structures was carried out using a variant of Dill's model, inspired by convex function, taking into account both hydrophobic (Dill's model) and polar interactions, including more structural information to reach better folding solutions. While using correlated networks, degrees of freedom in the folding of sequences were limited; as a result in all cases, more successful structural trials were found in comparison to a homogeneous lattice. The majority of employed sequences were designed by our workgroup, two of them were folded with other approaches, and another is a modified version of a previous sequence, initial forms of the other two have been employed but without taking into account polar-polar contributions. Three of them are newly proposed, intended to test the conjoint hydrophobic-hydrophobic and polar-polar contributions in crowded spaces. One sequence turned out to be the most difficult of the seven folded, this perhaps due to intrinsic (i) degrees of freedom and (ii) motifs of the expected 2D HP structure. Meanwhile two-sequence, although optimal folding was not achieved for neither of the two approaches, folding with correlated network approach not only produced better results than homogeneous space, but for them the best values found with crowding were very close to the expected optimal fitness. In general, five sequences were better folded with medium lattice units for correlated media; instead, another two sequences were better folded with a bit larger degree of lattice unit, revealing that depending on the degrees of freedom and particular folding, motifs in each sequence would require tuned crowding to achieve better folding. Therefore, the main goal herein was to obtain a modified 2D HP lattice model to mimic folding of proteins or secondary structures, like β-sheets, taking into account both hydrophobic-hydrophobic and polar-polar interactions, and fold them in a crowded environment. This simple but enough construction would be conducted to determine the needed information to fold sequences in a sort of a minimal but complete heuristic model. Finally, we claim that all folded sequences into crowded spaces achieve better results than homogeneous ones.
Collapse
Affiliation(s)
- Hiram Isaac Beltrán
- Departamento de Ciencias Básicas, Universidad Autónoma Metropolitana, Unidad Azcapotzalco, CDMX 02200, Mexico, Mexico
| | - Salomón J Alas-Guardado
- Departamento de Ciencias Naturales, Universidad Autónoma Metropolitana Unidad Cuajimalpa, CDMX 05300, Mexico, Mexico.
| | - Pedro Pablo González-Pérez
- Departamento de Matemáticas Aplicadas y Sistemas, Universidad Autónoma Metropolitana, Unidad Cuajimalpa, CDMX 05300, Mexico, Mexico.
| |
Collapse
|
5
|
The folding and misfolding mechanisms of multidomain proteins. MEDICINE IN DRUG DISCOVERY 2022. [DOI: 10.1016/j.medidd.2022.100126] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
|
6
|
Structural and thermodynamical insights into the binding and inhibition of FIH-1 by the N-terminal disordered region of Mint3. J Biol Chem 2021; 297:101304. [PMID: 34655613 PMCID: PMC8571082 DOI: 10.1016/j.jbc.2021.101304] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2021] [Revised: 10/09/2021] [Accepted: 10/12/2021] [Indexed: 11/29/2022] Open
Abstract
Mint3 is known to enhance aerobic ATP production, known as the Warburg effect, by binding to FIH-1. Since this effect is considered to be beneficial for cancer cells, the interaction is a promising target for cancer therapy. However, previous research has suggested that the interacting region of Mint3 with FIH-1 is intrinsically disordered, which makes investigation of this interaction challenging. Therefore, we adopted thermodynamic and structural studies in solution to clarify the structural and thermodynamical changes of Mint3 binding to FIH-1. First, using a combination of circular dichroism, nuclear magnetic resonance, and hydrogen/deuterium exchange–mass spectrometry (HDX-MS), we confirmed that the N-terminal half, which is the interacting part of Mint3, is mostly disordered. Next, we revealed a large enthalpy and entropy change in the interaction of Mint3 using isothermal titration calorimetry (ITC). The profile is consistent with the model that the flexibility of disordered Mint3 is drastically reduced upon binding to FIH-1. Moreover, we performed a series of ITC experiments with several types of truncated Mint3s, an effective approach since the interacting part of Mint3 is disordered, and identified amino acids 78 to 88 as a novel core site for binding to FIH-1. The truncation study of Mint3 also revealed the thermodynamic contribution of each part of Mint3 to the interaction with FIH-1, where the core sites contribute to the affinity (ΔG), while other sites only affect enthalpy (ΔH), by forming noncovalent bonds. This insight can serve as a foothold for further investigation of intrinsically disordered regions (IDRs) and drug development for cancer therapy.
Collapse
|
7
|
Nie QM, Sun LZ, Li HB, Chu X, Wang J. Effects of electrostatic interactions on global folding and local conformational dynamics of a multidomain Y-family DNA polymerase. Phys Chem Chem Phys 2021; 23:20841-20847. [PMID: 34533560 DOI: 10.1039/d1cp02832d] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
The Y-family DNA polymerases specialize in translesion DNA synthesis, which is essential for replicating damaged DNA. The Y-family polymerases, which are made up of four stable domains, exhibit extensive distributions of charged residues, and are responsible for the tight formation of the protein-DNA complex. However, it is still unclear how the electrostatic interactions influence the conformational dynamics of the polymerases. Here, we focus on the case of a prototype Y-family DNA polymerase, Dpo4. Using coarse-grained models including a salt-dependent electrostatic potential, we investigate the effects of the electrostatic interactions on the folding process of Dpo4. Our simulations show that strong electrostatic interactions result in a three-state folding of Dpo4, consistent with the experimental observations. This folding process exhibits low cooperativity led by low salt concentration, where the individual domains fold one by one through one single pathway. Since the refined folding order of domains in multidomain proteins can shrink the configurational space, we suggest that the electrostatic interactions facilitate the Dpo4 folding. In addition, we study the local conformational dynamics of Dpo4 in terms of fluctuation and frustration analyses. We show that the electrostatic interactions can exaggerate the local conformational properties, which are in favor of the large-scale conformational transition of Dpo4 during the functional DNA binding. Our results underline the importance of electrostatic interactions in the conformational dynamics of Dpo4 at both the global and local scale, providing useful guidance in protein engineering at the multidomain level.
Collapse
Affiliation(s)
- Qing-Miao Nie
- Department of Applied Physics, Zhejiang University of Technology, 288, Liuhe Road, Hangzhou 310023, P. R. China
| | - Li-Zhen Sun
- Department of Applied Physics, Zhejiang University of Technology, 288, Liuhe Road, Hangzhou 310023, P. R. China
| | - Hai-Bin Li
- Department of Applied Physics, Zhejiang University of Technology, 288, Liuhe Road, Hangzhou 310023, P. R. China
| | - Xiakun Chu
- Department of Chemistry, State University of New York at Stony Brook, Stony Brook, New York 11794, USA.
| | - Jin Wang
- Department of Chemistry, State University of New York at Stony Brook, Stony Brook, New York 11794, USA.
| |
Collapse
|
8
|
Giulini M, Rigoli M, Mattiotti G, Menichetti R, Tarenzi T, Fiorentini R, Potestio R. From System Modeling to System Analysis: The Impact of Resolution Level and Resolution Distribution in the Computer-Aided Investigation of Biomolecules. Front Mol Biosci 2021; 8:676976. [PMID: 34164432 PMCID: PMC8215203 DOI: 10.3389/fmolb.2021.676976] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2021] [Accepted: 05/06/2021] [Indexed: 12/18/2022] Open
Abstract
The ever increasing computer power, together with the improved accuracy of atomistic force fields, enables researchers to investigate biological systems at the molecular level with remarkable detail. However, the relevant length and time scales of many processes of interest are still hardly within reach even for state-of-the-art hardware, thus leaving important questions often unanswered. The computer-aided investigation of many biological physics problems thus largely benefits from the usage of coarse-grained models, that is, simplified representations of a molecule at a level of resolution that is lower than atomistic. A plethora of coarse-grained models have been developed, which differ most notably in their granularity; this latter aspect determines one of the crucial open issues in the field, i.e. the identification of an optimal degree of coarsening, which enables the greatest simplification at the expenses of the smallest information loss. In this review, we present the problem of coarse-grained modeling in biophysics from the viewpoint of system representation and information content. In particular, we discuss two distinct yet complementary aspects of protein modeling: on the one hand, the relationship between the resolution of a model and its capacity of accurately reproducing the properties of interest; on the other hand, the possibility of employing a lower resolution description of a detailed model to extract simple, useful, and intelligible information from the latter.
Collapse
Affiliation(s)
- Marco Giulini
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| | - Marta Rigoli
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| | - Giovanni Mattiotti
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| | - Roberto Menichetti
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| | - Thomas Tarenzi
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| | - Raffaele Fiorentini
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| | - Raffaello Potestio
- Physics Department, University of Trento, Trento, Italy.,INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| |
Collapse
|
9
|
Zhong Q, Li G. Adaptively Iterative Multiscale Switching Simulation Strategy and Applications to Protein Folding and Structure Prediction. J Phys Chem Lett 2021; 12:3151-3162. [PMID: 33755493 DOI: 10.1021/acs.jpclett.1c00618] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Structure prediction is an important means to quickly understand new protein functions. However, the prediction of effects of proteins that have no detectable templates is still to be improved. Molecular dynamics simulation is supposed to be the primary research tool for structure predictions, but it still has limitations of huge computational cost in all-atom (AA) models and rough accuracy in coarse-grained (CG) models. We propose a universal multiscale simulation strategy named AIMS in which simulations can iteratively switch among multiple resolutions in order to adaptively trade off AA accuracy and CG high-efficiency. AIMS follows the idea of CG-guided enhanced sampling so that final results always keep AA accuracy. We successfully achieve four ab initio and four data-assisted protein structure predictions using AIMS. The prediction result is an ensemble rather than a structure and provides special insights on folding metastable states. AIMS is estimated to achieve a computational speed about 40 times faster than that of conventional AA simulations.
Collapse
Affiliation(s)
- Qinglu Zhong
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Guohui Li
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| |
Collapse
|
10
|
Alas-Guardado SJ, González-Pérez PP, Beltrán HI. Contributions of topological polar-polar contacts to achieve better folding stability of 2D/3D HP lattice proteins: An in silico approach. AIMS BIOPHYSICS 2021. [DOI: 10.3934/biophy.2021023] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
<abstract>
<p>Many of the simplistic hydrophobic-polar lattice models, such as Dill's model (called <bold>Model 1</bold> herein), are aimed to fold structures through hydrophobic-hydrophobic interactions mimicking the well-known hydrophobic collapse present in protein structures. In this work, we studied 11 designed hydrophobic-polar sequences, S<sub>1</sub>-S<sub>8</sub> folded in 2D-square lattice, and S<sub>9</sub>-S<sub>11</sub> folded in 3D-cubic lattice. And to better fold these structures we have developed <bold>Model 2</bold> as an approximation to convex function aimed to weight hydrophobic-hydrophobic but also polar-polar contacts as an augmented version of <bold>Model 1</bold>. In this partitioned approach hydrophobic-hydrophobic ponderation was tuned as <italic>α</italic>-1 and polar-polar ponderation as <italic>α</italic>. This model is centered in preserving required hydrophobic substructure, and at the same time including polar-polar interactions, otherwise absent, to reach a better folding score now also acquiring the polar-polar substructure. In all tested cases the folding trials were better achieved with <bold>Model 2</bold>, using <italic>α</italic> values of 0.05, 0.1, 0.2 and 0.3 depending of sequence size, even finding optimal scores not reached with <bold>Model 1</bold>. An important result is that the better folding score, required the lower <italic>α</italic> weighting. And when <italic>α</italic> values above 0.3 are employed, no matter the nature of the hydrophobic-polar sequence, banning of hydrophobic-hydrophobic contacts started, thus yielding misfolding of sequences. Therefore, the value of <italic>α</italic> to correctly fold structures is the result of a careful weighting among hydrophobic-hydrophobic and polar-polar contacts.</p>
</abstract>
Collapse
|
11
|
Chu X, Suo Z, Wang J. Investigating the trade-off between folding and function in a multidomain Y-family DNA polymerase. eLife 2020; 9:60434. [PMID: 33079059 PMCID: PMC7641590 DOI: 10.7554/elife.60434] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2020] [Accepted: 10/16/2020] [Indexed: 01/01/2023] Open
Abstract
The way in which multidomain proteins fold has been a puzzling question for decades. Until now, the mechanisms and functions of domain interactions involved in multidomain protein folding have been obscure. Here, we develop structure-based models to investigate the folding and DNA-binding processes of the multidomain Y-family DNA polymerase IV (DPO4). We uncover shifts in the folding mechanism among ordered domain-wise folding, backtracking folding, and cooperative folding, modulated by interdomain interactions. These lead to ‘U-shaped’ DPO4 folding kinetics. We characterize the effects of interdomain flexibility on the promotion of DPO4–DNA (un)binding, which probably contributes to the ability of DPO4 to bypass DNA lesions, which is a known biological role of Y-family polymerases. We suggest that the native topology of DPO4 leads to a trade-off between fast, stable folding and tight functional DNA binding. Our approach provides an effective way to quantitatively correlate the roles of protein interactions in conformational dynamics at the multidomain level.
Collapse
Affiliation(s)
- Xiakun Chu
- Department of Chemistry, State University of New York at Stony Brook, New York, United States
| | - Zucai Suo
- Department of Biomedical Sciences, College of Medicine, Florida State University, Tallahassee, United States
| | - Jin Wang
- Department of Chemistry, State University of New York at Stony Brook, New York, United States
| |
Collapse
|
12
|
Aumpuchin P, Hamaue S, Kikuchi T. Prediction of the initial folding sites and the entire folding processes for Ig-like beta-sandwich proteins. Proteins 2019; 88:740-758. [PMID: 31833097 DOI: 10.1002/prot.25862] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2019] [Revised: 11/13/2019] [Accepted: 12/06/2019] [Indexed: 11/12/2022]
Abstract
Describing the whole story of protein folding is currently the main enigmatic problem in molecular bioinformatics study. Protein folding mechanisms have been intensively investigated with experimental as well as simulation techniques. Since a protein folds into its specific 3D structure from a unique amino acid sequence, it is interesting to extract as much information as possible from the amino acid sequence of a protein. Analyses based on inter-residue average distance statistics and a coarse-grained Gō-model simulation were conducted on Ig and FN3 domains of a titin protein to decode the folding mechanisms from their sequence data and native structure data, respectively. The central region of all domains was predicted to be an initial folding unit, that is, stable in an early state of folding. This common feature coincides well with the experimental results and underscores the significance of the β-sandwich proteins' common structure, namely, the key strands for folding and the Greek-key motif, which is located in the central region. We confirmed that our sequence-based techniques were able to predict the initial folding event just next to the denatured state and that a 3D-based Gō-model simulation can be used to investigate the whole process of protein folding.
Collapse
Affiliation(s)
- Panyavut Aumpuchin
- Department of Bioinformatics, College of Life Sciences, Ritsumeikan University, Kusatsu, Shiga, Japan
| | - Shoya Hamaue
- Department of Bioinformatics, College of Life Sciences, Ritsumeikan University, Kusatsu, Shiga, Japan
| | - Takeshi Kikuchi
- Department of Bioinformatics, College of Life Sciences, Ritsumeikan University, Kusatsu, Shiga, Japan
| |
Collapse
|
13
|
Gasic AG, Boob MM, Prigozhin MB, Homouz D, Wirth AJ, Daugherty CM, Gruebele M, Cheung MS. Critical phenomena in the temperature-pressure-crowding phase diagram of a protein. PHYSICAL REVIEW. X 2019; 9:041035. [PMID: 32642303 PMCID: PMC7343146 DOI: 10.1103/physrevx.9.041035] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
In the cell, proteins fold and perform complex functions through global structural rearrangements. Function requires a protein to be at the brink of stability to be susceptible to small environmental fluctuations, yet stable enough to maintain structural integrity. These apparently conflicting behaviors are exhibited by systems near a critical point, where distinct phases merge-a concept beyond previous studies indicating proteins have a well-defined folded/unfolded phase boundary in the pressure-temperature plane. Here, by modeling the protein phosphoglycerate kinase (PGK) on the temperature (T), pressure (P), and crowding volume-fraction (ϕ) phase diagram, we demonstrate a critical transition where phases merge, and PGK exhibits large structural fluctuations. Above the critical point, the difference between the intermediate and unfolded phases disappears. When ϕ increases, the critical point moves to lower T c. We verify the calculations with experiments mapping the T-P-ϕ space, which likewise reveal a critical point at 305 K and 170 MPa that moves to lower T c as ϕ increases. Crowding places PGK near a critical line in its natural parameter space, where large conformational changes can occur without costly free energy barriers. Specific structures are proposed for each phase based on simulation.
Collapse
Affiliation(s)
- Andrei G. Gasic
- University of Houston, Department of Physics, Houston, Texas, 77204, United States
- Center for Theoretical Biological Physics, Rice University, 77005, United States
| | - Mayank M. Boob
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Champaign, IL 61801, United States
| | - Maxim B. Prigozhin
- Department of Chemistry, University of Illinois at Urbana-Champaign, Champaign, IL, 61801, United States
| | - Dirar Homouz
- University of Houston, Department of Physics, Houston, Texas, 77204, United States
- Center for Theoretical Biological Physics, Rice University, 77005, United States
- Khalifa University of Science and Technology, Department of Physics, P.O. Box 127788, Abu Dhabi, United Arab Emirates
| | - Anna Jean Wirth
- Department of Chemistry, University of Illinois at Urbana-Champaign, Champaign, IL, 61801, United States
| | - Caleb M. Daugherty
- University of Houston, Department of Physics, Houston, Texas, 77204, United States
- Center for Theoretical Biological Physics, Rice University, 77005, United States
| | - Martin Gruebele
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Champaign, IL 61801, United States
- Department of Chemistry, University of Illinois at Urbana-Champaign, Champaign, IL, 61801, United States
- Department of Physics and Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Champaign, IL 61801, United States
| | - Margaret S. Cheung
- University of Houston, Department of Physics, Houston, Texas, 77204, United States
- Center for Theoretical Biological Physics, Rice University, 77005, United States
| |
Collapse
|
14
|
Observation of Continuous Contraction and a Metastable Misfolded State during the Collapse and Folding of a Small Protein. J Mol Biol 2019; 431:3814-3826. [DOI: 10.1016/j.jmb.2019.07.024] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2019] [Revised: 07/12/2019] [Accepted: 07/12/2019] [Indexed: 01/22/2023]
|
15
|
Alas SDJ, González-Pérez PP, Beltrán HI. In silico minimalist approach to study 2D HP protein folding into an inhomogeneous space mimicking osmolyte effect: First trial in the search of foldameric backbones. Biosystems 2019; 181:31-43. [PMID: 31029589 DOI: 10.1016/j.biosystems.2019.04.005] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2018] [Revised: 04/01/2019] [Accepted: 04/08/2019] [Indexed: 12/22/2022]
Abstract
We have employed our bioinformatics workbench, named Evolution, a Multi-Agent System based architecture with lattice-bead-models, evolutionary-algorithms, and correlated-networks as inhomogeneous spaces, with different correlation lengths, mimicking osmolyte effect (molecular crowding), to in silico survey protein folding. Resolution is with hydrophobic-polar (H-P) sequences in inhomogeneous 2D square lattices, since general biophysicochemical trends consider i) that the backbone is one of the major components responsible for protein folding and ii) osmolyte effect plays an important role to better folding kinetics and reach deeper optima. We have designed foldamers, as square n × n (n = 3, 4, 5, 6) arrays of hydrophobic cores stabilized by H⋯H contacts, attached through short PP (P2) or long PPPP (P4) loops, giving rise to 8 sequences (S1 to S8) with known optimal scores. Designed sequences were folded into different inhomogeneous spaces and indeed crowded media induced deeper optima, being crowding necessary to best fold, but the space should be enough constrained to induce folding without banning chain movement. The constrained space plays an important role to reach the optimal structure, depending on designed foldamer sequence size, for an optimal correlation length, implying that media affects the folding pathways as happens in real systems. Designed structures were found, moreover, they undergo to degenerated states, both folding states could survey considering i) backbone information and ii) osmolyte effect. In nature, the proteins fold in different structures aiming to reach a global minimum, but a local minimum could be enough to the protein to be functional or dysfunctional.
Collapse
Affiliation(s)
- Salomón de Jesús Alas
- Departamento de Ciencias Naturales, Universidad Autónoma Metropolitana Unidad Cuajimalpa, Ciudad de México, 05300, Mexico.
| | - Pedro Pablo González-Pérez
- Departamento de Matemáticas Aplicadas y Sistemas, Universidad Autónoma Metropolitana Unidad Cuajimalpa, 05300, Ciudad de Mexico, Mexico
| | - Hiram Isaac Beltrán
- Departamento de Ciencias Básicas, Universidad Autónoma Metropolitana Unidad Azcapotzalco, Ciudad de México, 02200, Mexico.
| |
Collapse
|
16
|
Aumpuchin P, Kikuchi T. Prediction of folding mechanisms for Ig-like beta sandwich proteins based on inter-residue average distance statistics methods. Proteins 2018; 87:120-135. [PMID: 30520530 DOI: 10.1002/prot.25637] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2018] [Revised: 10/05/2018] [Accepted: 11/29/2018] [Indexed: 11/11/2022]
Abstract
To understand the folding mechanism of a protein is one of the goals in bioinformatics study. Nowadays, it is enigmatic and difficult to extract folding information from amino acid sequence using standard bioinformatics techniques or even experimental protocols which can be time consuming. To overcome these problems, we aim to extract the initial folding unit for titin protein (Ig and fnIII domains) by means of inter-residue average distance statistics, Average Distance Map (ADM) and contact frequency analysis (F-value). TI I27 and TNfn3 domains are used to represent the Ig-domain and fnIII-domain, respectively. Beta-strands 2, 3, 5, and 6 are significant for the initial folding processes of TI I27. The central strands of TNfn3 were predicted as a primary folding segment. Known 3D structure and unknown 3D structure domains were investigated by structure or non-structure based multiple sequence alignment, respectively, to learn the conserved hydrophobic residues and predicted compact region relevant to evolution. Our results show good correspondence to experimental data, phi-value and protection factor from H-D exchange experiments. The significance of conserved hydrophobic residues near F-value peaks for structural stability using hydrophobic packing is confirmed. Our prediction methods once again could extract a folding mechanism only knowing the amino acid sequence.
Collapse
Affiliation(s)
- Panyavut Aumpuchin
- Department of Bioinformatics, College of Life Sciences, Ritsumeikan University, Kusatsu, Shiga, Japan
| | - Takeshi Kikuchi
- Department of Bioinformatics, College of Life Sciences, Ritsumeikan University, Kusatsu, Shiga, Japan
| |
Collapse
|
17
|
Wang Y, Tian P, Boomsma W, Lindorff-Larsen K. Monte Carlo Sampling of Protein Folding by Combining an All-Atom Physics-Based Model with a Native State Bias. J Phys Chem B 2018; 122:11174-11185. [DOI: 10.1021/acs.jpcb.8b06335] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Yong Wang
- Structural Biology and NMR Laboratory, Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Ole Maaløes Vej 5, DK-2200 Copenhagen N, Denmark
| | - Pengfei Tian
- Structural Biology and NMR Laboratory, Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Ole Maaløes Vej 5, DK-2200 Copenhagen N, Denmark
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Wouter Boomsma
- Department of Computer Science, University of Copenhagen, 2100 Copenhagen Ø, Denmark
| | - Kresten Lindorff-Larsen
- Structural Biology and NMR Laboratory, Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Ole Maaløes Vej 5, DK-2200 Copenhagen N, Denmark
| |
Collapse
|
18
|
Cheng C, Wu J, Liu G, Shi S, Chen T. Effects of Non-native Interactions on Frustrated Proteins Folding under Confinement. J Phys Chem B 2018; 122:7654-7667. [DOI: 10.1021/acs.jpcb.8b04147] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Chenqian Cheng
- Key Laboratory of Synthetic and Natural Functional Molecular Chemistry of the Ministry of Education, College of Chemistry and Materials Science, Northwest University, Xi’an 710127, P. R. China
| | - Jing Wu
- Key Laboratory of Synthetic and Natural Functional Molecular Chemistry of the Ministry of Education, College of Chemistry and Materials Science, Northwest University, Xi’an 710127, P. R. China
| | - Gaoyuan Liu
- Key Laboratory of Synthetic and Natural Functional Molecular Chemistry of the Ministry of Education, College of Chemistry and Materials Science, Northwest University, Xi’an 710127, P. R. China
| | - Suqing Shi
- Key Laboratory of Synthetic and Natural Functional Molecular Chemistry of the Ministry of Education, College of Chemistry and Materials Science, Northwest University, Xi’an 710127, P. R. China
| | - Tao Chen
- Key Laboratory of Synthetic and Natural Functional Molecular Chemistry of the Ministry of Education, College of Chemistry and Materials Science, Northwest University, Xi’an 710127, P. R. China
| |
Collapse
|
19
|
Wu J, Cheng C, Liu G, Zhang P, Chen T. The folding pathways and thermodynamics of semiflexible polymers. J Chem Phys 2018; 148:184901. [PMID: 29764123 DOI: 10.1063/1.5018114] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
Inspired by the protein folding and DNA packing, we have systematically studied the thermodynamic and kinetic behaviors of single semiflexible homopolymers by Langevin dynamics simulations. In line with experiments, a rich variety of folding products, such as rod-like bundles, hairpins, toroids, and a mixture of them, are observed in the complete diagram of states. Moreover, knotted structures with a significant population are found in a certain range of bending stiffness in thermal equilibrium. As the solvent quality becomes poorer, the population of the intermediate occurring in the folding process increases, which leads to a severe chevron rollover for the folding arm. However, the population of the intermediates in the unfolding process is very low, insufficient to induce unfolding arm rollover. The total types of folding pathways from the coil state to the toroidal state for a semiflexible polymer chain remain unchanged by varying the solvent quality or temperature, whereas the kinetic partitioning into different folding events can be tuned significantly. In the process of knotting, three types of mechanisms, namely, plugging, slipknotting, and sliding, are discovered. Along the folding evolution, a semiflexible homopolymer chain can knot at any stage of folding upon leaving the extended coil state, and the probability to find a knot increases with chain compactness. In addition, we find rich types of knotted topologies during the folding of a semiflexible homopolymer chain. This study should be helpful in gaining insight into the general principles of biopolymer folding.
Collapse
Affiliation(s)
- Jing Wu
- Key Laboratory of Synthetic and Natural Functional Molecular Chemistry of the Ministry of Education, College of Chemistry and Materials Science, Northwest University, Xi'an 710127, People's Republic of China
| | - Chenqian Cheng
- Key Laboratory of Synthetic and Natural Functional Molecular Chemistry of the Ministry of Education, College of Chemistry and Materials Science, Northwest University, Xi'an 710127, People's Republic of China
| | - Gaoyuan Liu
- Key Laboratory of Synthetic and Natural Functional Molecular Chemistry of the Ministry of Education, College of Chemistry and Materials Science, Northwest University, Xi'an 710127, People's Republic of China
| | - Ping Zhang
- Key Laboratory of Synthetic and Natural Functional Molecular Chemistry of the Ministry of Education, College of Chemistry and Materials Science, Northwest University, Xi'an 710127, People's Republic of China
| | - Tao Chen
- Key Laboratory of Synthetic and Natural Functional Molecular Chemistry of the Ministry of Education, College of Chemistry and Materials Science, Northwest University, Xi'an 710127, People's Republic of China
| |
Collapse
|
20
|
Gopi S, Paul S, Ranu S, Naganathan AN. Extracting the Hidden Distributions Underlying the Mean Transition State Structures in Protein Folding. J Phys Chem Lett 2018; 9:1771-1777. [PMID: 29565127 DOI: 10.1021/acs.jpclett.8b00538] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
The inherent conflict between noncovalent interactions and the large conformational entropy of the polypeptide chain forces folding reactions and their mechanisms to deviate significantly from chemical reactions. Accordingly, measures of structure in the transition state ensemble (TSE) are strongly influenced by the underlying distributions of microscopic folding pathways that are challenging to discern experimentally. Here, we present a detailed analysis of 150,000 folding transition paths of five proteins at three different thermodynamic conditions from an experimentally consistent statistical mechanical model. We find that the underlying TSE structural distributions are rarely unimodal, and the average experimental measures arise from complex underlying distributions. Unfolding pathways also exhibit subtle differences from folding counterparts due to a combination of Hammond behavior and native-state movements. Local interactions and topological complexity, to a lesser extent, are found to determine pathway heterogeneity, underscoring the importance of the balance between local and nonlocal energetics in protein folding.
Collapse
|
21
|
Ouyang Y, Zhao L, Zhang Z. Characterization of the structural ensembles of p53 TAD2 by molecular dynamics simulations with different force fields. Phys Chem Chem Phys 2018. [DOI: 10.1039/c8cp00067k] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
The conformations of p53 TAD2 in complexes and sampled in simulations with five force fields.
Collapse
Affiliation(s)
- Yanhua Ouyang
- College of Life Science, University of Chinese Academy of Sciences
- Beijing
- China
| | - Likun Zhao
- College of Life Science, University of Chinese Academy of Sciences
- Beijing
- China
| | - Zhuqing Zhang
- College of Life Science, University of Chinese Academy of Sciences
- Beijing
- China
| |
Collapse
|
22
|
Song J, Gomes GN, Shi T, Gradinaru CC, Chan HS. Conformational Heterogeneity and FRET Data Interpretation for Dimensions of Unfolded Proteins. Biophys J 2017; 113:1012-1024. [PMID: 28877485 DOI: 10.1016/j.bpj.2017.07.023] [Citation(s) in RCA: 48] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2017] [Revised: 07/26/2017] [Accepted: 07/31/2017] [Indexed: 12/28/2022] Open
Abstract
A mathematico-physically valid formulation is required to infer properties of disordered protein conformations from single-molecule Förster resonance energy transfer (smFRET). Conformational dimensions inferred by conventional approaches that presume a homogeneous conformational ensemble can be unphysical. When all possible-heterogeneous as well as homogeneous-conformational distributions are taken into account without prejudgment, a single value of average transfer efficiency 〈E〉 between dyes at two chain ends is generally consistent with highly diverse, multiple values of the average radius of gyration 〈Rg〉. Here we utilize unbiased conformational statistics from a coarse-grained explicit-chain model to establish a general logical framework to quantify this fundamental ambiguity in smFRET inference. As an application, we address the long-standing controversy regarding the denaturant dependence of 〈Rg〉 of unfolded proteins, focusing on Protein L as an example. Conventional smFRET inference concluded that 〈Rg〉 of unfolded Protein L is highly sensitive to [GuHCl], but data from SAXS suggested a near-constant 〈Rg〉 irrespective of [GuHCl]. Strikingly, our analysis indicates that although the reported 〈E〉 values for Protein L at [GuHCl] = 1 and 7 M are very different at 0.75 and 0.45, respectively, the Bayesian Rg2 distributions consistent with these two 〈E〉 values overlap by as much as 75%. Our findings suggest, in general, that the smFRET-SAXS discrepancy regarding unfolded protein dimensions likely arise from highly heterogeneous conformational ensembles at low or zero denaturant, and that additional experimental probes are needed to ascertain the nature of this heterogeneity.
Collapse
Affiliation(s)
- Jianhui Song
- School of Polymer Science and Engineering, Qingdao University of Science and Technology, Qingdao, Shandong, China; Departments of Biochemistry and Molecular Genetics, University of Toronto, Toronto, Ontario, Canada
| | - Gregory-Neal Gomes
- Department of Chemical and Physical Sciences, University of Toronto Mississauga, Mississauga, Ontario, Canada; Department of Physics, University of Toronto, Toronto, Ontario, Canada
| | - Tongfei Shi
- State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, Jilin, China
| | - Claudiu C Gradinaru
- Department of Chemical and Physical Sciences, University of Toronto Mississauga, Mississauga, Ontario, Canada; Department of Physics, University of Toronto, Toronto, Ontario, Canada
| | - Hue Sun Chan
- Departments of Biochemistry and Molecular Genetics, University of Toronto, Toronto, Ontario, Canada.
| |
Collapse
|
23
|
Wu J, Chen G, Zhang Z, Zhang P, Chen T. The low populated folding intermediate of a mutant of the Fyn SH3 domain identified by a simple model. Phys Chem Chem Phys 2017; 19:22321-22328. [DOI: 10.1039/c7cp04139j] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
The low populated on-pathway folding intermediate of the A39V/N53P/V55L Fyn SH3 domain is captured by a native-centric model augmented by sequence-dependent nonnative hydrophobic interactions.
Collapse
Affiliation(s)
- Jing Wu
- Key Laboratory of Synthetic and Natural Functional Molecular Chemistry of the Ministry of Education
- College of Chemistry and Materials Science
- Northwest University
- Xi'an
- P. R. China
| | - Guojun Chen
- Key Laboratory of Synthetic and Natural Functional Molecular Chemistry of the Ministry of Education
- College of Chemistry and Materials Science
- Northwest University
- Xi'an
- P. R. China
| | - Zhuqing Zhang
- College of Life Sciences
- University of Chinese Academy of Sciences
- Beijing
- P. R. China
| | - Ping Zhang
- Key Laboratory of Synthetic and Natural Functional Molecular Chemistry of the Ministry of Education
- College of Chemistry and Materials Science
- Northwest University
- Xi'an
- P. R. China
| | - Tao Chen
- Key Laboratory of Synthetic and Natural Functional Molecular Chemistry of the Ministry of Education
- College of Chemistry and Materials Science
- Northwest University
- Xi'an
- P. R. China
| |
Collapse
|