1
|
Zayats V, Sikora M, Perlinska AP, Stasiulewicz A, Gren BA, Sulkowska JI. Conservation of knotted and slipknotted topology in transmembrane transporters. Biophys J 2023; 122:4528-4541. [PMID: 37919904 PMCID: PMC10719070 DOI: 10.1016/j.bpj.2023.10.031] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2023] [Revised: 08/25/2023] [Accepted: 10/30/2023] [Indexed: 11/04/2023] Open
Abstract
The existence of nontrivial topology is well accepted in globular proteins but not in membrane proteins. Our comprehensive topological analysis of the Protein Data Bank structures reveals 18 families of transmembrane proteins with nontrivial topology, showing that they constitute a significant number of membrane proteins. Moreover, we found that they comprise one of the largest groups of secondary active transporters. We classified them based on their knotted fingerprint into four groups: three slipknotted and one knotted. Unexpectedly, we found that the same protein can possess two distinct slipknot motifs that correspond to its outward- and inward-open conformational state. Based on the analysis of structures and knotted fingerprints, we show that slipknot topology is directly involved in the conformational transition and substrate transfer. Therefore, entanglement can be used to classify proteins and to find their structure-function relationship. Furthermore, based on the topological analysis of the transmembrane protein structures predicted by AlphaFold, we identified new potentially slipknotted protein families.
Collapse
Affiliation(s)
- Vasilina Zayats
- Centre of New Technologies, University of Warsaw, Warsaw, Poland
| | - Maciej Sikora
- Centre of New Technologies, University of Warsaw, Warsaw, Poland; Faculty of Mathematics, Informatics and Mechanics, University of Warsaw, Warsaw, Poland
| | | | - Adam Stasiulewicz
- Centre of New Technologies, University of Warsaw, Warsaw, Poland; Department of Drug Chemistry, Faculty of Pharmacy, Medical University of Warsaw, Warsaw, Poland
| | - Bartosz A Gren
- Centre of New Technologies, University of Warsaw, Warsaw, Poland
| | | |
Collapse
|
2
|
Salicari L, Baiesi M, Orlandini E, Trovato A. Folding kinetics of an entangled protein. PLoS Comput Biol 2023; 19:e1011107. [PMID: 37956216 PMCID: PMC10681328 DOI: 10.1371/journal.pcbi.1011107] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2023] [Revised: 11/27/2023] [Accepted: 11/02/2023] [Indexed: 11/15/2023] Open
Abstract
The possibility of the protein backbone adopting lasso-like entangled motifs has attracted increasing attention. After discovering the surprising abundance of natively entangled protein domain structures, it was shown that misfolded entangled subpopulations might become thermosensitive or escape the homeostasis network just after translation. To investigate the role of entanglement in shaping folding kinetics, we introduce a novel indicator and analyze simulations of a coarse-grained, structure-based model for two small single-domain proteins. The model recapitulates the well-known two-state folding mechanism of a non-entangled SH3 domain. However, despite its small size, a natively entangled antifreeze RD1 protein displays a rich refolding behavior, populating two distinct kinetic intermediates: a short-lived, entangled, near-unfolded state and a longer-lived, non-entangled, near-native state. The former directs refolding along a fast pathway, whereas the latter is a kinetic trap, consistently with known experimental evidence of two different characteristic times. Upon trapping, the natively entangled loop folds without being threaded by the N-terminal residues. After trapping, the native entangled structure emerges by either backtracking to the unfolded state or threading through the already formed but not yet entangled loop. Along the fast pathway, trapping does not occur because the native contacts at the closure of the lasso-like loop fold after those involved in the N-terminal thread, confirming previous predictions. Despite this, entanglement may appear already in unfolded configurations. Remarkably, a longer-lived, near-native intermediate, with non-native entanglement properties, recalls what was observed in cotranslational folding.
Collapse
Affiliation(s)
- Leonardo Salicari
- Department of Physics and Astronomy “G. Galilei”, University of Padova, Padova, Italy
- National Institute of Nuclear Physics (INFN), Padova Section, Padova, Italy
| | - Marco Baiesi
- Department of Physics and Astronomy “G. Galilei”, University of Padova, Padova, Italy
- National Institute of Nuclear Physics (INFN), Padova Section, Padova, Italy
| | - Enzo Orlandini
- Department of Physics and Astronomy “G. Galilei”, University of Padova, Padova, Italy
- National Institute of Nuclear Physics (INFN), Padova Section, Padova, Italy
| | - Antonio Trovato
- Department of Physics and Astronomy “G. Galilei”, University of Padova, Padova, Italy
- National Institute of Nuclear Physics (INFN), Padova Section, Padova, Italy
| |
Collapse
|
3
|
Especial JNC, Faísca PFN. Effects of sequence-dependent non-native interactions in equilibrium and kinetic folding properties of knotted proteins. J Chem Phys 2023; 159:065101. [PMID: 37551809 DOI: 10.1063/5.0160886] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2023] [Accepted: 07/24/2023] [Indexed: 08/09/2023] Open
Abstract
Determining the role of non-native interactions in folding dynamics, kinetics, and mechanisms is a classic problem in protein folding. More recently, this question has witnessed a renewed interest in light of the hypothesis that knotted proteins require the assistance of non-native interactions to fold efficiently. Here, we conduct extensive equilibrium and kinetic Monte Carlo simulations of a simple off-lattice C-alpha model to explore the role of non-native interactions in the thermodynamics and kinetics of three proteins embedding a trefoil knot in their native structure. We find that equilibrium knotted conformations are stabilized by non-native interactions that are non-local, and proximal to native ones, thus enhancing them. Additionally, non-native interactions increase the knotting frequency at high temperatures, and in partially folded conformations below the transition temperatures. Although non-native interactions clearly enhance the efficiency of transition from an unfolded conformation to a partially folded knotted one, they are not required to efficiently fold a knotted protein. Indeed, a native-centric interaction potential drives the most efficient folding transition, provided that the simulation temperature is well below the transition temperature of the considered model system.
Collapse
Affiliation(s)
- João N C Especial
- Departamento de Física, Faculdade de Ciências, Ed. C8, Universidade de Lisboa, Campo Grande, Lisboa, Portugal
- BioISI-Biosystems and Integrative Sciences Institute, Faculdade de Ciências, Universidade de Lisboa, Campo Grande, Lisboa, Portugal
| | - Patrícia F N Faísca
- Departamento de Física, Faculdade de Ciências, Ed. C8, Universidade de Lisboa, Campo Grande, Lisboa, Portugal
- BioISI-Biosystems and Integrative Sciences Institute, Faculdade de Ciências, Universidade de Lisboa, Campo Grande, Lisboa, Portugal
| |
Collapse
|
4
|
A Note on the Effects of Linear Topology Preservation in Monte Carlo Simulations of Knotted Proteins. Int J Mol Sci 2022; 23:ijms232213871. [PMID: 36430350 PMCID: PMC9695063 DOI: 10.3390/ijms232213871] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2022] [Revised: 11/07/2022] [Accepted: 11/08/2022] [Indexed: 11/12/2022] Open
Abstract
Monte Carlo simulations are a powerful technique and are widely used in different fields. When applied to complex molecular systems with long chains, such as those in synthetic polymers and proteins, they have the advantage of providing a fast and computationally efficient way to sample equilibrium ensembles and calculate thermodynamic and structural properties under desired conditions. Conformational Monte Carlo techniques employ a move set to perform the transitions in the simulation Markov chain. While accepted conformations must preserve the sequential bonding of the protein chain model and excluded volume among its units, the moves themselves may take the chain across itself. We call this a break in linear topology preservation. In this manuscript, we show, using simple protein models, that there is no difference in equilibrium properties calculated with a move set that preserves linear topology and one that does not. However, for complex structures, such as those of deeply knotted proteins, the preservation of linear topology provides correct equilibrium results but only after long relaxation. In any case, to analyze folding pathways, knotting mechanisms and folding kinetics, the preservation of linear topology may be an unavoidable requirement.
Collapse
|
5
|
Kumari K, Ravi Prakash J, Padinhateeri R. Heterogeneous interactions and polymer entropy decide organization and dynamics of chromatin domains. Biophys J 2022; 121:2794-2812. [PMID: 35672951 PMCID: PMC9382282 DOI: 10.1016/j.bpj.2022.06.008] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2022] [Revised: 04/28/2022] [Accepted: 06/01/2022] [Indexed: 11/02/2022] Open
Abstract
Chromatin is known to be organized into multiple domains of varying sizes and compaction. While these domains are often imagined as static structures, they are highly dynamic and show cell-to-cell variability. Since processes such as gene regulation and DNA replication occur in the context of these domains, it is important to understand their organization, fluctuation, and dynamics. To simulate chromatin domains, one requires knowledge of interaction strengths among chromatin segments. Here, we derive interaction-strength parameters from experimentally known contact maps and use them to predict chromatin organization and dynamics. Taking two domains on the human chromosome as examples, we investigate its three-dimensional organization, size/shape fluctuations, and dynamics of different segments within a domain, accounting for hydrodynamic effects. Considering different cell types, we quantify changes in interaction strengths and chromatin shape fluctuations in different epigenetic states. Perturbing the interaction strengths systematically, we further investigate how epigenetic-like changes can alter the spatio-temporal nature of the domains. Our results show that heterogeneous weak interactions are crucial in determining the organization of the domains. Computing effective stiffness and relaxation times, we investigate how perturbations in interactions affect the solid- and liquid-like nature of chromatin domains. Quantifying dynamics of chromatin segments within a domain, we show how the competition between polymer entropy and interaction energy influence the timescales of loop formation and maintenance of stable loops.
Collapse
Affiliation(s)
- Kiran Kumari
- IITB-Monash Research Academy, Indian Institute of Technology Bombay, Mumbai, Maharashtra, 400076, India; Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Mumbai, Maharashtra 400076, India; Department of Chemical Engineering, Monash University, Melbourne, VIC 3800, Australia.
| | - J Ravi Prakash
- Department of Chemical Engineering, Monash University, Melbourne, VIC 3800, Australia.
| | - Ranjith Padinhateeri
- Department of Biosciences and Bioengineering, Indian Institute of Technology Bombay, Mumbai, Maharashtra 400076, India.
| |
Collapse
|
6
|
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
|
7
|
Slipknotted and unknotted monovalent cation-proton antiporters evolved from a common ancestor. PLoS Comput Biol 2021; 17:e1009502. [PMID: 34648493 PMCID: PMC8562792 DOI: 10.1371/journal.pcbi.1009502] [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: 06/30/2020] [Revised: 11/02/2021] [Accepted: 09/28/2021] [Indexed: 11/20/2022] Open
Abstract
While the slipknot topology in proteins has been known for over a decade, its evolutionary origin is still a mystery. We have identified a previously overlooked slipknot motif in a family of two-domain membrane transporters. Moreover, we found that these proteins are homologous to several families of unknotted membrane proteins. This allows us to directly investigate the evolution of the slipknot motif. Based on our comprehensive analysis of 17 distantly related protein families, we have found that slipknotted and unknotted proteins share a common structural motif. Furthermore, this motif is conserved on the sequential level as well. Our results suggest that, regardless of topology, the proteins we studied evolved from a common unknotted ancestor single domain protein. Our phylogenetic analysis suggests the presence of at least seven parallel evolutionary scenarios that led to the current diversity of proteins in question. The tools we have developed in the process can now be used to investigate the evolution of other repeated-domain proteins. In proteins with the slipknot topology, the polypeptide chain forms a slipknot—a structure that is not necessarily manifest to a naked eye, but it can be detected using mathematical methods. Slipknots are conserved motifs often found at catalytic sites and are directly involved in molecular transport. Although the first proteins with slipknots were found in 2007, many questions remain unanswered, e.g. how these proteins appeared, or whether the slipknotted proteins evolved from unknotted ones or vice versa. Here we provide the first analysis of homologous slipknotted and unknotted transmembrane proteins in order to elucidate their evolutionary relationship. We show that two-domain slipknotted and unknotted membrane transporters share the same one-domain unknotted protein as an ancestor. The ancestor gene duplicated and underwent various diversification and fusion events during the evolution, which have led to the appearance of a large superfamily of secondary active transporters. The slipknot motif seems to have been created by chance after a fusion of two single domain genes. Therefore, we show here that the slipknotted transporter evolved from an unknotted one-domain protein and that there are at least seven different evolutionary scenarios that gave rise to this large superfamily of transporters.
Collapse
|
8
|
Fonseka HYY, Javidi A, Oliveira LFL, Micheletti C, Stan G. Unfolding and Translocation of Knotted Proteins by Clp Biological Nanomachines: Synergistic Contribution of Primary Sequence and Topology Revealed by Molecular Dynamics Simulations. J Phys Chem B 2021; 125:7335-7350. [PMID: 34110163 DOI: 10.1021/acs.jpcb.1c00898] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We use Langevin dynamics simulations to model, at an atomistic resolution, how various natively knotted proteins are unfolded in repeated allosteric translocating cycles of the ClpY ATPase. We consider proteins representative of different topologies, from the simplest knot (trefoil 31), to the three-twist 52 knot, to the most complex stevedore, 61, knot. We harness the atomistic detail of the simulations to address aspects that have so far remained largely unexplored, such as sequence-dependent effects on the ruggedness of the landscape traversed during knot sliding. Our simulations reveal the combined effect on translocation of the knotted protein structure, i.e., backbone topology and geometry, and primary sequence, i.e., side chain size and interactions, and show that the latter can dominate translocation hindrance. In addition, we observe that due to the interplay between the knotted topology and intramolecular contacts the transmission of tension along the polypeptide chain occurs very differently from that of homopolymers. Finally, by considering native and non-native interactions, we examine how the disruption or formation of such contacts can affect the translocation processivity and concomitantly create multiple unfolding pathways with very different activation barriers.
Collapse
Affiliation(s)
| | - Alex Javidi
- Data Sciences, Janssen Research and Development, Spring House, Pennsylvania 19477, United States
| | - Luiz F L Oliveira
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221, United States
| | - Cristian Micheletti
- Molecular and Statistical Biophysics, Scuola Internazionale Superiore di Studi Avanzati (SISSA), 34136 Trieste, Italy
| | - George Stan
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221, United States
| |
Collapse
|
9
|
Especial JNC, Faísca PFN. A Specific Set of Heterogeneous Native Interactions Yields Efficient Knotting in Protein Folding. J Phys Chem B 2021; 125:7359-7367. [PMID: 34197706 DOI: 10.1021/acs.jpcb.1c03127] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Native interactions are crucial for folding, and non-native interactions appear to be critical for efficiently knotting proteins. Therefore, it is important to understand both their roles in the folding of knotted proteins. It has been proposed that non-native interactions drive the correct order of contact formation, which is essential to avoid backtracking and efficiently self-tie. In this study, we ask if non-native interactions are strictly necessary to tangle a protein or if the correct order of contact formation can be assured by a specific set of native, but otherwise heterogeneous (i.e., having distinct energies), interactions. In order to address this problem, we conducted extensive Monte Carlo simulations of lattice models of protein-like sequences designed to fold into a preselected knotted conformation embedding a trefoil knot. We were able to identify a specific set of heterogeneous native interactions that drives efficient knotting and is able to fold the protein when combined with the remaining native interactions modeled as homogeneous. This specific set of heterogeneous native interactions is strictly enough to efficiently self-tie. A distinctive feature of these native interactions is that they do not backtrack because their energies ensure the correct order of contact formation. Furthermore, they stabilize a knotted intermediate state, which is en route to the native structure. Our results thus show that-at least in the context of the adopted model-non-native interactions are not necessary to knot a protein. However, when they are taken into account in protein energetics, it is possible to find specific, nonlocal non-native interactions that operate as a scaffold that assists the knotting step.
Collapse
Affiliation(s)
- João N C Especial
- Departamento de Física, Universidade de Lisboa, Campo Grande, Ed. C8, 1749-016 Lisboa, Portugal.,BioISI-Biosystems and Integrative Sciences Institute, Faculdade de Ciências, Universidade de Lisboa, Campo Grande, Ed. C8, 1749-016 Lisboa, Portugal
| | - Patrícia F N Faísca
- Departamento de Física, Universidade de Lisboa, Campo Grande, Ed. C8, 1749-016 Lisboa, Portugal.,BioISI-Biosystems and Integrative Sciences Institute, Faculdade de Ciências, Universidade de Lisboa, Campo Grande, Ed. C8, 1749-016 Lisboa, Portugal
| |
Collapse
|
10
|
Michieletto D, Orlandini E, Turner MS, Micheletti C. Separation of Geometrical and Topological Entanglement in Confined Polymers Driven out of Equilibrium. ACS Macro Lett 2020; 9:1081-1085. [PMID: 35653213 DOI: 10.1021/acsmacrolett.0c00366] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
We use Brownian dynamics simulations and advanced topological profiling methods to characterize the out-of-equilibrium evolution of self-entanglement in linear polymers confined into nanochannels and under periodic compression. By introducing suitable observables, we can distinguish two main forms of entanglement that we term geometrical and topological. The latter is measured by the number of (essential) crossings of the physical knot detected after a suitable bridging of the chain termini. The former is instead measured as the average number of times a linear chain appears to cross itself when viewed under all projections and is irrespective of the physical knotted state. The key discovery of our work is that these two forms of entanglement are uncoupled and evolve with distinct dynamics. While geometrical entanglement is typically in phase with the compression-elongation cycles and it is primarily sensitive to its force f, the topological measure is mildly sensitive to cyclic modulation but strongly depends on both compression force f and duration k. The findings could assist the interpretation of experiments using fluorescence molecular tracers to track physical knots in polymers. Furthermore, we identify optimal regions in the experimentally controllable parameter space in which to obtain more/less topological and geometrical entanglement; this may help designing polymers with targeted topology.
Collapse
Affiliation(s)
- Davide Michieletto
- School of Physics and Astronomy, University of Edinburgh, Peter Guthrie Tait Road, Edinburgh, EH9 3FD, United Kingdom.,MRC Human Genetics Unit, Institute of Genetics and Molecular Medicine, University of Edinburgh, North Crewe Road, Edinburgh, EH4 2XU, United Kingdom
| | - Enzo Orlandini
- Dipartimento di Fisica e Astronomia and Sezione INFN, Universitá degli Studi di Padova, I-35131 Padova, Italy
| | - Matthew S Turner
- Department of Physics and Centre for Complexity Science, University of Warwick, Coventry, CV4 7AL, U.K.,Department of Chemical Engineering, Kyoto University, Kyoto, Japan
| | - Cristian Micheletti
- SISSA (Scuola Internazionale Superiore di Studi Avanzati), Via Bonomea 265, 34136 Trieste, Italy
| |
Collapse
|
11
|
Piejko M, Niewieczerzal S, Sulkowska JI. The Folding of Knotted Proteins: Distinguishing the Distinct Behavior of Shallow and Deep Knots. Isr J Chem 2020. [DOI: 10.1002/ijch.202000036] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Maciej Piejko
- Faculty of ChemistryUniversity of Warsaw Pasteura 1 Warsaw 02-093 Poland
- Centre of New TechnologiesUniversity of Warsaw Banacha 2c Warsaw 02-097 Poland
| | | | - Joanna I. Sulkowska
- Faculty of ChemistryUniversity of Warsaw Pasteura 1 Warsaw 02-093 Poland
- Centre of New TechnologiesUniversity of Warsaw Banacha 2c Warsaw 02-097 Poland
| |
Collapse
|
12
|
Rivera M, Hao Y, Maillard RA, Baez M. Mechanical unfolding of a knotted protein unveils the kinetic and thermodynamic consequences of threading a polypeptide chain. Sci Rep 2020; 10:9562. [PMID: 32533020 PMCID: PMC7292828 DOI: 10.1038/s41598-020-66258-5] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2019] [Accepted: 05/12/2020] [Indexed: 12/21/2022] Open
Abstract
Knots are remarkable topological features in nature. The presence of knots in crystallographic structures of proteins have stimulated considerable research to determine the kinetic and thermodynamic consequences of threading a polypeptide chain. By mechanically manipulating MJ0366, a small single domain protein harboring a shallow trefoil knot, we allow the protein to refold from either the knotted or the unknotted denatured state to characterize the free energy profile associated to both folding pathways. By comparing the stability of the native state with reference to the knotted and unknotted denatured state we find that knotting the polypeptide chain of MJ0366 increase the folding energy barrier in a magnitude close to the energy cost of forming a knot randomly in the denatured state. These results support that a protein knot can be formed during a single cooperative step of folding but occurs at the expenses of a large increment on the free energy barrier.
Collapse
Affiliation(s)
- Maira Rivera
- Departamento de Bioquímica y Biología Molecular, Facultad de Ciencias Químicas y Farmacéuticas, Universidad de Chile, Santiago, Chile
| | - Yuxin Hao
- Department of Chemistry, Georgetown University, Washington, DC, 20057, USA
| | - Rodrigo A Maillard
- Department of Chemistry, Georgetown University, Washington, DC, 20057, USA.
| | - Mauricio Baez
- Departamento de Bioquímica y Biología Molecular, Facultad de Ciencias Químicas y Farmacéuticas, Universidad de Chile, Santiago, Chile.
| |
Collapse
|
13
|
Sulkowska JI. On folding of entangled proteins: knots, lassos, links and θ-curves. Curr Opin Struct Biol 2020; 60:131-141. [PMID: 32062143 DOI: 10.1016/j.sbi.2020.01.007] [Citation(s) in RCA: 36] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2019] [Revised: 01/02/2020] [Accepted: 01/12/2020] [Indexed: 12/15/2022]
Abstract
Around 6% of protein structures deposited in the PDB are entangled, forming knots, slipknots, lassos, links, and θ-curves. In each of these cases, the protein backbone weaves through itself in a complex way, and at some point passes through a closed loop, formed by other regions of the protein structure. Such a passing can be interpreted as crossing a topological barrier. How proteins overcome such barriers, and therefore different degrees of frustration, challenged scientists and has shed new light on the field of protein folding. In this review, we summarize the current knowledge about the free energy landscape of proteins with non-trivial topology. We describe identified mechanisms which lead proteins to self-tying. We discuss the influence of excluded volume, such as crowding and chaperones, on tying, based on available data. We briefly discuss the diversity of topological complexity of proteins and their evolution. We also list available tools to investigate non-trivial topology. Finally, we formulate intriguing and challenging questions at the boundary of biophysics, bioinformatics, biology, and mathematics, which arise from the discovery of entangled proteins.
Collapse
Affiliation(s)
- Joanna Ida Sulkowska
- Centre of New Technologies, University of Warsaw, Warsaw, Poland; Faculty of Chemistry, University of Warsaw, Warsaw, Poland.
| |
Collapse
|
14
|
Perego C, Potestio R. Computational methods in the study of self-entangled proteins: a critical appraisal. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2019; 31:443001. [PMID: 31269476 DOI: 10.1088/1361-648x/ab2f19] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
The existence of self-entangled proteins, the native structure of which features a complex topology, unveils puzzling, and thus fascinating, aspects of protein biology and evolution. The discovery that a polypeptide chain can encode the capability to self-entangle in an efficient and reproducible way during folding, has raised many questions, regarding the possible function of these knots, their conservation along evolution, and their role in the folding paradigm. Understanding the function and origin of these entanglements would lead to deep implications in protein science, and this has stimulated the scientific community to investigate self-entangled proteins for decades by now. In this endeavour, advanced experimental techniques are more and more supported by computational approaches, that can provide theoretical guidelines for the interpretation of experimental results, and for the effective design of new experiments. In this review we provide an introduction to the computational study of self-entangled proteins, focusing in particular on the methodological developments related to this research field. A comprehensive collection of techniques is gathered, ranging from knot theory algorithms, that allow detection and classification of protein topology, to Monte Carlo or molecular dynamics strategies, that constitute crucial instruments for investigating thermodynamics and kinetics of this class of proteins.
Collapse
Affiliation(s)
- Claudio Perego
- Max Panck Institute for Polymer Research, Ackermannweg 10, Mainz 55128, Germany
| | | |
Collapse
|
15
|
Xu Y, Li S, Yan Z, Ge B, Huang F, Yue T. Revealing Cooperation between Knotted Conformation and Dimerization in Protein Stabilization by Molecular Dynamics Simulations. J Phys Chem Lett 2019; 10:5815-5822. [PMID: 31525988 DOI: 10.1021/acs.jpclett.9b02209] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
The topological knot is thought to play a stabilizing role in maintaining the global fold and nature of proteins with the underlying mechanism yet to be elucidated. Given that most proteins containing trefoil knots exist and function as homodimers with a large part of the dimer interface occupied by the knotted region, we reason that the knotted conformation cooperates with dimerization in protein stabilization. Here, we take YbeA from Escherichia coli as the knotted protein model, using molecular dynamics (MD) simulations to compare the stability of two pairs of dimeric proteins having the same sequence and secondary structures but differing in the presence or absence of a trefoil knot in each subunit. The dimer interface of YbeA is identified to involve favorable contacts among three α-helices (α1, α3, and α5), one of which (α5) is threaded through a loop connected with α3 to form the knot. Upon removal of the knot by appropriate change of the knot-making crossing of the polypeptide chain, relevant domains are less constrained and exhibit enhanced fluctuations to decrease contacts at the interface. Unknotted subunits are less compact and undergo structural changes to ease the dimer separation. Such a stabilizing effect is evidenced by steered MD simulations, showing that the mechanical force required for dimer separation is significantly reduced by removing the knot. In addition to the knotted conformation, dimerization further improves the protein stability by restricting the α1-α5 separation, which is defined as a leading step for protein unfolding. These results provide important insights into the structure-function relationship of dimerization in knotted proteins.
Collapse
Affiliation(s)
- Yan Xu
- State Key Laboratory of Heavy Oil Processing, Center for Bioengineering and Biotechnology, College of Chemical Engineering , China University of Petroleum (East China) , Qingdao 266580 , China
- College of Electronic Engineering and Automation , Shandong University of Science and Technology , Qingdao 266590 , China
| | - Shixin Li
- State Key Laboratory of Heavy Oil Processing, Center for Bioengineering and Biotechnology, College of Chemical Engineering , China University of Petroleum (East China) , Qingdao 266580 , China
| | - Zengshuai Yan
- State Key Laboratory of Heavy Oil Processing, Center for Bioengineering and Biotechnology, College of Chemical Engineering , China University of Petroleum (East China) , Qingdao 266580 , China
| | - Baosheng Ge
- State Key Laboratory of Heavy Oil Processing, Center for Bioengineering and Biotechnology, College of Chemical Engineering , China University of Petroleum (East China) , Qingdao 266580 , China
| | - Fang Huang
- State Key Laboratory of Heavy Oil Processing, Center for Bioengineering and Biotechnology, College of Chemical Engineering , China University of Petroleum (East China) , Qingdao 266580 , China
| | - Tongtao Yue
- State Key Laboratory of Heavy Oil Processing, Center for Bioengineering and Biotechnology, College of Chemical Engineering , China University of Petroleum (East China) , Qingdao 266580 , China
| |
Collapse
|
16
|
Especial J, Nunes A, Rey A, Faísca PF. Hydrophobic confinement modulates thermal stability and assists knotting in the folding of tangled proteins. Phys Chem Chem Phys 2019; 21:11764-11775. [PMID: 31114834 DOI: 10.1039/c9cp01701a] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
There is growing support for the idea that the in vivo folding process of knotted proteins is assisted by chaperonins, but the mechanism of chaperonin assisted folding remains elusive. Here, we conduct extensive Monte Carlo simulations of lattice and off-lattice models to explore the effects of confinement and hydrophobic intermolecular interactions with the chaperonin cage in the folding and knotting processes. We find that moderate to high protein-cavity interactions (which are likely to be established in the beginning of the chaperonin working cycle) cause an energetic destabilization of the protein that overcomes the entropic stabilization driven by excluded volume, and leads to a decrease of the melting temperature relative to bulk conditions. Moreover, mild-to-moderate hydrophobic interactions with the cavity (which would be established later in the cycle) lead to a significant enhancement of knotting probability in relation to bulk conditions while simultaneously moderating the effect of steric confinement in the enhancement of thermal stability. Our results thus indicate that the chaperonin may be able to assist knotting without simultaneously thermally stabilizing potential misfolded states to a point that would hamper productive folding thus compromising its functional role.
Collapse
Affiliation(s)
- João Especial
- Departamento de Física, Universidade de Lisboa, Campo Grande, Ed. C8, 1749-016 Lisboa, Portugal.
| | | | | | | |
Collapse
|
17
|
Perego C, Potestio R. Searching the Optimal Folding Routes of a Complex Lasso Protein. Biophys J 2019; 117:214-228. [PMID: 31235180 PMCID: PMC6700606 DOI: 10.1016/j.bpj.2019.05.025] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2019] [Revised: 04/29/2019] [Accepted: 05/30/2019] [Indexed: 10/27/2022] Open
Abstract
Understanding how polypeptides can efficiently and reproducibly attain a self-entangled conformation is a compelling biophysical challenge that might shed new light on our general knowledge of protein folding. Complex lassos, namely self-entangled protein structures characterized by a covalent loop sealed by a cysteine bridge, represent an ideal test system in the framework of entangled folding. Indeed, because cysteine bridges form in oxidizing conditions, they can be used as on/off switches of the structure topology to investigate the role played by the backbone entanglement in the process. In this work, we have used molecular dynamics to simulate the folding of a complex lasso glycoprotein, granulocyte-macrophage colony-stimulating factor, modeling both reducing and oxidizing conditions. Together with a well-established Gō-like description, we have employed the elastic folder model, a coarse-grained, minimalistic representation of the polypeptide chain driven by a structure-based angular potential. The purpose of this study is to assess the kinetically optimal pathways in relation to the formation of the native topology. To this end, we have implemented an evolutionary strategy that tunes the elastic folder model potentials to maximize the folding probability within the early stages of the dynamics. The resulting protein model is capable of folding with high success rate, avoiding the kinetic traps that hamper the efficient folding in the other tested models. Employing specifically designed topological descriptors, we could observe that the selected folding routes avoid the topological bottleneck by locking the cysteine bridge after the topology is formed. These results provide valuable insights on the selection of mechanisms in self-entangled protein folding while, at the same time, the proposed methodology can complement the usage of established minimalistic models and draw useful guidelines for more detailed simulations.
Collapse
Affiliation(s)
- Claudio Perego
- Polymer Theory Department, Max Planck Institute for Polymer Research, Mainz, Germany.
| | - Raffaello Potestio
- Department of Physics, University of Trento, Trento, Italy; INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, Trento, Italy
| |
Collapse
|
18
|
Sivertsson EM, Jackson SE, Itzhaki LS. The AAA+ protease ClpXP can easily degrade a 3 1 and a 5 2-knotted protein. Sci Rep 2019; 9:2421. [PMID: 30787316 PMCID: PMC6382783 DOI: 10.1038/s41598-018-38173-3] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2018] [Accepted: 12/04/2018] [Indexed: 12/16/2022] Open
Abstract
Knots in proteins are hypothesized to make them resistant to enzymatic degradation by ATP-dependent proteases and recent studies have shown that whereas ClpXP can easily degrade a protein with a shallow 31 knot, it cannot degrade 52-knotted proteins if degradation is initiated at the C-terminus. Here, we present detailed studies of the degradation of both 31- and 52-knotted proteins by ClpXP using numerous constructs where proteins are tagged for degradation at both N- and C-termini. Our results confirm and extend earlier work and show that ClpXP can easily degrade a deeply 31-knotted protein. In contrast to recently published work on the degradation of 52-knotted proteins, our results show that the ClpXP machinery can also easily degrade these proteins. However, the degradation depends critically on the location of the degradation tag and the local stability near the tag. Our results are consistent with mechanisms in which either the knot simply slips along the polypeptide chain and falls off the free terminus, or one in which the tightened knot enters the translocation pore of ClpXP. Results of experiments on knotted protein fusions with a highly stable domain show partial degradation and the formation of degradation intermediates.
Collapse
Affiliation(s)
- Elin M Sivertsson
- Department of Pharmacology, Tennis Court Road, Cambridge, CB2 1PD, UK
| | - Sophie E Jackson
- Department of Chemistry, Lensfield Road, Cambridge, CB2 1EW, UK.
| | - Laura S Itzhaki
- Department of Pharmacology, Tennis Court Road, Cambridge, CB2 1PD, UK.
| |
Collapse
|
19
|
|
20
|
Sulkowska JI, Sułkowski P. Entangled Proteins: Knots, Slipknots, Links, and Lassos. SPRINGER SERIES IN SOLID-STATE SCIENCES 2018. [DOI: 10.1007/978-3-319-76596-9_8] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
21
|
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
|
22
|
Jackson SE, Suma A, Micheletti C. How to fold intricately: using theory and experiments to unravel the properties of knotted proteins. Curr Opin Struct Biol 2016; 42:6-14. [PMID: 27794211 DOI: 10.1016/j.sbi.2016.10.002] [Citation(s) in RCA: 56] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2016] [Revised: 09/06/2016] [Accepted: 10/03/2016] [Indexed: 11/15/2022]
Abstract
Over the years, advances in experimental and computational methods have helped us to understand the role of thermodynamic, kinetic and active (chaperone-aided) effects in coordinating the folding steps required to achieving a knotted native state. Here, we review such developments by paying particular attention to the complementarity of experimental and computational studies. Key open issues that could be tackled with either or both approaches are finally pointed out.
Collapse
Affiliation(s)
- Sophie E Jackson
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom.
| | - Antonio Suma
- SISSA, International School for Advanced Studies, via Bonomea 265, I-34136 Trieste, Italy
| | - Cristian Micheletti
- SISSA, International School for Advanced Studies, via Bonomea 265, I-34136 Trieste, Italy.
| |
Collapse
|
23
|
Najafi S, Potestio R. Folding of small knotted proteins: Insights from a mean field coarse-grained model. J Chem Phys 2016; 143:243121. [PMID: 26723606 DOI: 10.1063/1.4934541] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
A small but relevant number of proteins whose native structure is known features nontrivial topology, i.e., they are knotted. Understanding the process of folding from a swollen unknotted state to the biologically relevant native conformation is, for these proteins, particularly difficult, due to their rate-limiting topological entanglement. To shed some light into this conundrum, we introduced a structure-based coarse-grained model of the protein, where the information about the folded conformation is encoded in bonded angular interactions only, which do not favor the formation of native contacts. A stochastic search scheme in parameter space is employed to identify a set of interactions that maximizes the probability to attain the knotted state. The optimal knotting pathways of the two smallest knotted proteins, obtained through this approach, are consistent with the results derived by means of coarse-grained as well as full atomistic simulations.
Collapse
Affiliation(s)
- Saeed Najafi
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Raffaello Potestio
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| |
Collapse
|
24
|
Wang F, Cazzolli G, Wintrode P, Faccioli P. Folding Mechanism of Proteins Im7 and Im9: Insight from All-Atom Simulations in Implicit and Explicit Solvent. J Phys Chem B 2016; 120:9297-307. [PMID: 27532482 DOI: 10.1021/acs.jpcb.6b05819] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
Im7 and Im9 are evolutionary related proteins with almost identical native structures. In spite of their structural similarity, experiments show that Im7 folds through a long-lived on-pathway intermediate, while Im9 folds according to two-state kinetics. In this work, we use a recently developed enhanced path sampling method to generate many folding trajectories for these proteins, using realistic atomistic force fields, in both implicit and explicit solvent. Overall, our results are in good agreement with the experimental ϕ values and with the result of ϕ-value-restrained molecular dynamics (MD) simulations. However, our implicit solvent simulations fail to predict a qualitative difference in the folding pathways of Im7 and Im9. In contrast, our simulations in explicit solvent correctly reproduce the fact that only protein Im7 folds through a on-pathway intermediate. By analyzing our atomistic trajectories, we provide a physical picture which explains the observed difference in the folding kinetics of these chains.
Collapse
Affiliation(s)
- F Wang
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy , Baltimore, Maryland 21201, United States
| | - G Cazzolli
- Physics Department, University of Trento , via Sommarive 14 Povo, Trento 38128, Italy
| | - P Wintrode
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy , Baltimore, Maryland 21201, United States
| | - P Faccioli
- Physics Department, University of Trento , via Sommarive 14 Povo, Trento 38128, Italy
| |
Collapse
|
25
|
Folding analysis of the most complex Stevedore's protein knot. Sci Rep 2016; 6:31514. [PMID: 27527519 PMCID: PMC4985754 DOI: 10.1038/srep31514] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2016] [Accepted: 07/21/2016] [Indexed: 12/21/2022] Open
Abstract
DehI is a homodimeric haloacid dehalogenase from Pseudomonas putida that contains the most complex 61 Stevedore's protein knot within its folding topology. To examine how DehI attains such an intricate knotted topology we combined far-UV circular dichroism (CD), intrinsic fluorescence spectroscopy and small angle X-ray scattering (SAXS) to investigate its folding mechanism. Equilibrium unfolding of DehI by chemical denaturation indicated the presence of two highly populated folding intermediates, I and I'. While the two intermediates vary in secondary structure contents and tertiary packing according to CD and intrinsic fluorescence, respectively, their overall dimension and compactness are similar according to SAXS. Three single-tryptophan variants (W34, W53, and W196) were generated to probe non-cooperative unfolding events localized around the three fluorophores. Kinetic fluorescence measurements indicated that the transition from the intermediate I' to the unfolded state is rate limiting. Our multiparametric folding analyses suggest that DehI unfolds through a linear folding pathway with two distinct folding intermediates by initial hydrophobic collapse followed by nucleation condensation, and that knotting precedes the formation of secondary structures.
Collapse
|
26
|
Abstract
Spontaneous folding of a polypeptide chain into a knotted structure remains one of the most puzzling and fascinating features of protein folding. The folding of knotted proteins is on the timescale of minutes and thus hard to reproduce with atomistic simulations that have been able to reproduce features of ultrafast folding in great detail. Furthermore, it is generally not possible to control the topology of the unfolded state. Single-molecule force spectroscopy is an ideal tool for overcoming this problem: by variation of pulling directions, we controlled the knotting topology of the unfolded state of the 52-knotted protein ubiquitin C-terminal hydrolase isoenzyme L1 (UCH-L1) and have therefore been able to quantify the influence of knotting on its folding rate. Here, we provide direct evidence that a threading event associated with formation of either a 31 or 52 knot, or a step closely associated with it, significantly slows down the folding of UCH-L1. The results of the optical tweezers experiments highlight the complex nature of the folding pathway, many additional intermediate structures being detected that cannot be resolved by intrinsic fluorescence. Mechanical stretching of knotted proteins is also of importance for understanding the possible implications of knots in proteins for cellular degradation. Compared with a simple 31 knot, we measure a significantly larger size for the 52 knot in the unfolded state that can be further tightened with higher forces. Our results highlight the potential difficulties in degrading a 52 knot compared with a 31 knot.
Collapse
|
27
|
Potestio R, Tubiana L. Discretized knot motion on a tensioned fiber induced by transverse waves. SOFT MATTER 2016; 12:669-673. [PMID: 26510521 DOI: 10.1039/c5sm01766a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
Topological entanglement is a ubiquitous feature of many biological as well as artificial polymers and fibers. While the equilibrium properties of entangled chains have been the subject of several studies, little is known about their out-of-equilibrium behavior. Here, we address the problem of a stretched knotted fiber driven by a periodic force applied to one of its termini. We show that the onset of standing waves kinetically traps the knot in spatially localized states where the amplitude of the oscillations is maximal, while the knot normal diffusive dynamics is replaced by a discrete jump dynamics.
Collapse
Affiliation(s)
- Raffaello Potestio
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany.
| | | |
Collapse
|
28
|
Dabrowski-Tumanski P, Jarmolinska AI, Sulkowska JI. Prediction of the optimal set of contacts to fold the smallest knotted protein. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2015; 27:354109. [PMID: 26291339 DOI: 10.1088/0953-8984/27/35/354109] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Knotted protein chains represent a new motif in protein folds. They have been linked to various diseases, and recent extensive analysis of the Protein Data Bank shows that they constitute 1.5% of all deposited protein structures. Despite thorough theoretical and experimental investigations, the role of knots in proteins still remains elusive. Nonetheless, it is believed that knots play an important role in mechanical and thermal stability of proteins. Here, we perform a comprehensive analysis of native, shadow-specific and non-native interactions which describe free energy landscape of the smallest knotted protein (PDB id 2efv). We show that the addition of shadow-specific contacts in the loop region greatly enhances folding kinetics, while the addition of shadow-specific contacts along the C-terminal region (H3 or H4) results in a new folding route with slower kinetics. By means of direct coupling analysis (DCA) we predict non-native contacts which also can accelerate kinetics. Next, we show that the length of the C-terminal knot tail is responsible for the shape of the free energy barrier, while the influence of the elongation of the N-terminus is not significant. Finally, we develop a concept of a minimal contact map sufficient for 2efv protein to fold and analyze properties of this protein using this map.
Collapse
Affiliation(s)
- P Dabrowski-Tumanski
- Centre of New Technologies, University of Warsaw, Banacha 2c, 02-097 Warsaw, Poland. Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | | | | |
Collapse
|
29
|
Lim NCH, Jackson SE. Molecular knots in biology and chemistry. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2015; 27:354101. [PMID: 26291690 DOI: 10.1088/0953-8984/27/35/354101] [Citation(s) in RCA: 99] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Knots and entanglements are ubiquitous. Beyond their aesthetic appeal, these fascinating topological entities can be either useful or cumbersome. In recent decades, the importance and prevalence of molecular knots have been increasingly recognised by scientists from different disciplines. In this review, we provide an overview on the various molecular knots found in naturally occurring biological systems (DNA, RNA and proteins), and those created by synthetic chemists. We discuss the current knowledge in these fields, including recent developments in experimental and, in some cases, computational studies which are beginning to shed light into the complex interplay between the structure, formation and properties of these topologically intricate molecules.
Collapse
Affiliation(s)
- Nicole C H Lim
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK. Faculty of Sciences, Universiti Brunei Darussalam, Gadong BE 1410, Brunei Darussalam
| | | |
Collapse
|
30
|
Faísca PF. Knotted proteins: A tangled tale of Structural Biology. Comput Struct Biotechnol J 2015; 13:459-68. [PMID: 26380658 PMCID: PMC4556803 DOI: 10.1016/j.csbj.2015.08.003] [Citation(s) in RCA: 58] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2015] [Revised: 07/31/2015] [Accepted: 08/07/2015] [Indexed: 01/19/2023] Open
Abstract
Knotted proteins have their native structures arranged in the form of an open knot. In the last ten years researchers have been making significant efforts to reveal their folding mechanism and understand which functional advantage(s) knots convey to their carriers. Molecular simulations have been playing a fundamental role in this endeavor, and early computational predictions about the knotting mechanism have just been confirmed in wet lab experiments. Here we review a collection of simulation results that allow outlining the current status of the field of knotted proteins, and discuss directions for future research.
Collapse
|
31
|
Computing Reaction Pathways of Rare Biomolecular Transitions using Atomistic Force-Fields. Biophys Chem 2015; 208:62-7. [PMID: 26320390 DOI: 10.1016/j.bpc.2015.06.014] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2015] [Revised: 06/23/2015] [Accepted: 06/23/2015] [Indexed: 02/02/2023]
Abstract
The Dominant Reaction Pathway (DRP) method is an approximate variational scheme which can be used to compute reaction pathways in conformational transitions undergone by large biomolecules (up to ~10(3) amino-acids) using realistic all-atom force fields. We first review the status of development of this method. Next, we discuss its validation against the results of plain MD protein folding simulations performed by the DE-Shaw group using the Anton supercomputer. Finally, we review a few representative applications of the DRP approach to study reactions which are far too complex and rare to be investigated by plain MD, even on the Anton machine.
Collapse
|
32
|
Chen T, Song J, Chan HS. Theoretical perspectives on nonnative interactions and intrinsic disorder in protein folding and binding. Curr Opin Struct Biol 2014; 30:32-42. [PMID: 25544254 DOI: 10.1016/j.sbi.2014.12.002] [Citation(s) in RCA: 64] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2014] [Revised: 12/02/2014] [Accepted: 12/02/2014] [Indexed: 11/29/2022]
Abstract
The diverse biological functions of intrinsically disordered proteins (IDPs) have markedly raised our appreciation of protein conformational versatility, whereas the existence of energetically favorable yet functional detrimental nonnative interactions underscores the physical limitations of evolutionary optimization. Here we survey recent advances in using biophysical modeling to gain insight into experimentally observed nonnative behaviors and IDP properties. Simulations of IDP interactions to date focus mostly on coupled folding-binding, which follows essentially the same organizing principle as the local-nonlocal coupling mechanism in cooperative folding of monomeric globular proteins. By contrast, more innovative theories of electrostatic and aromatic interactions are needed for the conceptually novel but less-explored 'fuzzy' complexes in which the functionally bound IDPs remain largely disordered.
Collapse
Affiliation(s)
- Tao Chen
- Department of Biochemistry, University of Toronto, Toronto, Ontario M5S 1A8, Canada; Department of Molecular Genetics, University of Toronto, Toronto, Ontario M5S 1A8, Canada; Department of Physics, University of Toronto, Toronto, Ontario M5S 1A7, Canada
| | - Jianhui Song
- Department of Biochemistry, University of Toronto, Toronto, Ontario M5S 1A8, Canada; Department of Molecular Genetics, University of Toronto, Toronto, Ontario M5S 1A8, Canada; Department of Physics, University of Toronto, Toronto, Ontario M5S 1A7, Canada
| | - Hue Sun Chan
- Department of Biochemistry, University of Toronto, Toronto, Ontario M5S 1A8, Canada; Department of Molecular Genetics, University of Toronto, Toronto, Ontario M5S 1A8, Canada; Department of Physics, University of Toronto, Toronto, Ontario M5S 1A7, Canada.
| |
Collapse
|
33
|
Di Stefano M, Tubiana L, Di Ventra M, Micheletti C. Driving knots on DNA with AC/DC electric fields: topological friction and memory effects. SOFT MATTER 2014; 10:6491-6498. [PMID: 25048107 DOI: 10.1039/c4sm00160e] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Abstract
The dynamical properties of entangled polyelectrolytes are investigated theoretically and computationally for a proposed novel micromanipulation setup. Specifically, we investigate the effects of DC and AC electric fields acting longitudinally on knotted DNA chains, modelled as semiflexible chains of charged beads, under mechanical tension. We consider various experimentally accessible values of the field amplitude and frequency as well as several of the simplest knot types. In particular, we consider both torus and twist knots because they are respectively known to be able or unable to slide along macroscopic threads and ropes. Strikingly, this qualitative distinction disappears in this microscopic context because all the considered knot types acquire a systematic drift in the direction of the electric force. Notably, the knot drift velocity and diffusion coefficient in zero field (both measurable also experimentally) can be used to define a characteristic "frictional" lengthscale for the various knot types. This previously unexplored length provides valuable information on the extent of self-interactions in the nominal knotted region. It is finally observed that the motion of a knot can effectively follow the AC field only if the driving period is larger than the knot relaxation time (for which the self-diffusion time provides an upper bound). These results suggest that salient aspects of the intrinsic dynamics of knots in DNA chains could be probed experimentally by means of external, time-dependent electric fields.
Collapse
Affiliation(s)
- Marco Di Stefano
- SISSA - Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, 34136 Trieste, Italy.
| | | | | | | |
Collapse
|
34
|
Mashaghi A, van Wijk R, Tans S. Circuit Topology of Proteins and Nucleic Acids. Structure 2014; 22:1227-1237. [DOI: 10.1016/j.str.2014.06.015] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2013] [Revised: 05/10/2014] [Accepted: 06/17/2014] [Indexed: 01/19/2023]
|
35
|
Yadahalli S, Hemanth Giri Rao VV, Gosavi S. Modeling Non-Native Interactions in Designed Proteins. Isr J Chem 2014. [DOI: 10.1002/ijch.201400035] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|
36
|
Haglund E, Sulkowska JI, Noel JK, Lammert H, Onuchic JN, Jennings PA. Pierced Lasso Bundles are a new class of knot-like motifs. PLoS Comput Biol 2014; 10:e1003613. [PMID: 24945798 PMCID: PMC4063663 DOI: 10.1371/journal.pcbi.1003613] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2014] [Accepted: 03/26/2014] [Indexed: 01/11/2023] Open
Abstract
A four-helix bundle is a well-characterized motif often used as a target for designed pharmaceutical therapeutics and nutritional supplements. Recently, we discovered a new structural complexity within this motif created by a disulphide bridge in the long-chain helical bundle cytokine leptin. When oxidized, leptin contains a disulphide bridge creating a covalent-loop through which part of the polypeptide chain is threaded (as seen in knotted proteins). We explored whether other proteins contain a similar intriguing knot-like structure as in leptin and discovered 11 structurally homologous proteins in the PDB. We call this new helical family class the Pierced Lasso Bundle (PLB) and the knot-like threaded structural motif a Pierced Lasso (PL). In the current study, we use structure-based simulation to investigate the threading/folding mechanisms for all the PLBs along with three unthreaded homologs as the covalent loop (or lasso) in leptin is important in folding dynamics and activity. We find that the presence of a small covalent loop leads to a mechanism where structural elements slipknot to thread through the covalent loop. Larger loops use a piercing mechanism where the free terminal plugs through the covalent loop. Remarkably, the position of the loop as well as its size influences the native state dynamics, which can impact receptor binding and biological activity. This previously unrecognized complexity of knot-like proteins within the helical bundle family comprises a completely new class within the knot family, and the hidden complexity we unraveled in the PLBs is expected to be found in other protein structures outside the four-helix bundles. The insights gained here provide critical new elements for future investigation of this emerging class of proteins, where function and the energetic landscape can be controlled by hidden topology, and should be take into account in ab initio predictions of newly identified protein targets.
Collapse
Affiliation(s)
- Ellinor Haglund
- Center for Theoretical Biological Physics (CTBP) and Department of Physics, University of California at San Diego (UCSD), La Jolla, California, United States of America
- Center for Theoretical Biological Physics (CTBP) and Departments of Physics and Astronomy, Chemistry and Biochemistry and Cell Biology, Rice University, Houston, Texas, United States of America
| | | | - Jeffrey K. Noel
- Center for Theoretical Biological Physics (CTBP) and Departments of Physics and Astronomy, Chemistry and Biochemistry and Cell Biology, Rice University, Houston, Texas, United States of America
| | - Heiko Lammert
- Center for Theoretical Biological Physics (CTBP) and Departments of Physics and Astronomy, Chemistry and Biochemistry and Cell Biology, Rice University, Houston, Texas, United States of America
| | - José N. Onuchic
- Center for Theoretical Biological Physics (CTBP) and Departments of Physics and Astronomy, Chemistry and Biochemistry and Cell Biology, Rice University, Houston, Texas, United States of America
| | - Patricia A. Jennings
- Departments of Chemistry and Biochemistry, University of California at San Diego (UCSD), La Jolla, California, United States of America
- * E-mail:
| |
Collapse
|
37
|
Tubiana L. Computational study on the progressive factorization of composite polymer knots into separated prime components. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 89:052602. [PMID: 25353821 DOI: 10.1103/physreve.89.052602] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/28/2013] [Indexed: 06/04/2023]
Abstract
Using Monte Carlo simulations and advanced knot localization methods, we analyze the length and distribution of prime components in composite knots tied on freely jointed rings. For increasing contour length, we observe the progressive factorization of composite knots into separated prime components. However, we observe that a complete factorization, equivalent to the "decorated ring" picture, is not obtained even for rings of contour lengths N ≃ 3 N(0), about tens of times the most probable length of the prime knots tied on the rings. The decorated ring hypothesis has been used in the literature to justify the factorization of composite knot probabilities into the knotting probabilities of their prime components. Following our results, we suggest that such a hypothesis may not be necessary to explain the factorization of the knotting probabilities, at least when polymers excluding volume is not relevant. We rationalize the behavior of the system through a simple one-dimensional model in which prime knots are replaced by slip links randomly placed on a circle, with the only constraint being that the length of the loops has the same distribution as that of the length of the corresponding prime knots.
Collapse
Affiliation(s)
- Luca Tubiana
- Department of Theoretical Physics, Jožef Stefan Institute, Jamova cesta 39, SI-1000 Ljubljana (Slovenia) and SISSA-Scuola Internazionale Superiore di Studi Avanzati Via Bonomea 265, 34136 Trieste, Italy
| |
Collapse
|
38
|
Hills RD. Balancing bond, nonbond, and gō-like terms in coarse grain simulations of conformational dynamics. Methods Mol Biol 2014; 1084:123-140. [PMID: 24061919 DOI: 10.1007/978-1-62703-658-0_7] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/02/2023]
Abstract
Characterization of the protein conformational landscape remains a challenging problem, whether it concerns elucidating folding mechanisms, predicting native structures or modeling functional transitions. Coarse-grained molecular dynamics simulation methods enable exhaustive sampling of the energetic landscape at resolutions of biological interest. The general utility of structure-based models is reviewed along with their differing levels of approximation. Simple Gō models incorporate attractive native interactions and repulsive nonnative contacts, resulting in an ideal smooth landscape. Non-Gō coarse-grained models reduce the parameter set as needed but do not include bias to any desired native structure. While non-Gō models have achieved limited success in protein coarse-graining, they can be combined with native structured-based potentials to create a balanced and powerful force field. Recent applications of such Gō-like models have yielded insight into complex folding mechanisms and conformational transitions in large macromolecules. The accuracy and usefulness of reduced representations are also revealed to be a function of the mathematical treatment of the intrinsic bonded topology.
Collapse
Affiliation(s)
- Ronald D Hills
- Department of Pharmaceutical Sciences, University of New England, Portland, ME, USA
| |
Collapse
|
39
|
Covino R, Skrbić T, Beccara SA, Faccioli P, Micheletti C. The role of non-native interactions in the folding of knotted proteins: insights from molecular dynamics simulations. Biomolecules 2013; 4:1-19. [PMID: 24970203 PMCID: PMC4030985 DOI: 10.3390/biom4010001] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2013] [Revised: 12/10/2013] [Accepted: 12/20/2013] [Indexed: 12/14/2022] Open
Abstract
For several decades, the presence of knots in naturally-occurring proteins was largely ruled out a priori for its supposed incompatibility with the efficiency and robustness of folding processes. For this very same reason, the later discovery of several unrelated families of knotted proteins motivated researchers to look into the physico-chemical mechanisms governing the concerted sequence of folding steps leading to the consistent formation of the same knot type in the same protein location. Besides experiments, computational studies are providing considerable insight into these mechanisms. Here, we revisit a number of such recent investigations within a common conceptual and methodological framework. By considering studies employing protein models with different structural resolution (coarse-grained or atomistic) and various force fields (from pure native-centric to realistic atomistic ones), we focus on the role of native and non-native interactions. For various unrelated instances of knotted proteins, non-native interactions are shown to be very important for favoring the emergence of conformations primed for successful self-knotting events.
Collapse
Affiliation(s)
- Roberto Covino
- Department of Physics, University of Trento, Via Sommarive 14, Trento 38123, Italy.
| | - Tatjana Skrbić
- SISSA-Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, Trieste 34136, Italy.
| | - Silvio A Beccara
- Interdisciplinary Laboratory for Computational Science, FBK-CMM and University of Trento,Trento 38123, Italy.
| | - Pietro Faccioli
- Department of Physics, University of Trento, Via Sommarive 14, Trento 38123, Italy.
| | - Cristian Micheletti
- SISSA-Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, Trieste 34136, Italy.
| |
Collapse
|
40
|
Unfolding thermodynamics of cysteine-rich proteins and molecular thermal-adaptation of marine ciliates. Biomolecules 2013; 3:967-85. [PMID: 24970199 PMCID: PMC4030967 DOI: 10.3390/biom3040967] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2013] [Revised: 10/28/2013] [Accepted: 10/29/2013] [Indexed: 11/17/2022] Open
Abstract
Euplotes nobilii and Euplotes raikovi are phylogenetically closely allied species of marine ciliates, living in polar and temperate waters, respectively. Their evolutional relation and the sharply different temperatures of their natural environments make them ideal organisms to investigate thermal-adaptation. We perform a comparative study of the thermal unfolding of disulfide-rich protein pheromones produced by these ciliates. Recent circular dichroism (CD) measurements have shown that the two psychrophilic (E. nobilii) and mesophilic (E. raikovi) protein families are characterized by very different melting temperatures, despite their close structural homology. The enhanced thermal stability of the E. raikovi pheromones is realized notwithstanding the fact that these proteins form, as a rule, a smaller number of disulfide bonds. We perform Monte Carlo (MC) simulations in a structure-based coarse-grained (CG) model to show that the higher stability of the E. raikovi pheromones is due to the lower locality of the disulfide bonds, which yields a lower entropy increase in the unfolding process. Our study suggests that the higher stability of the mesophilic E. raikovi phermones is not mainly due to the presence of a strongly hydrophobic core, as it was proposed in the literature. In addition, we argue that the molecular adaptation of these ciliates may have occurred from cold to warm, and not from warm to cold. To provide a testable prediction, we identify a point-mutation of an E. nobilii pheromone that should lead to an unfolding temperature typical of that of E. raikovi pheromones.
Collapse
|
41
|
Abstract
Most proteins, in order to perform their biological function, have to fold to a compact native state. The increasing number of knotted and slipknotted proteins identified suggests that proteins are able to manoeuvre around topological barriers during folding. In the present article, we review the current progress in elucidating the knotting process in proteins. Although we concentrate on theoretical approaches, where a knotted topology can be unambiguously detected, comparison with experiments is also reviewed. Numerical simulations suggest that the folding process for small knotted proteins is composed of twisted loop formation and then threading by either slipknot geometries or flipping. As the size of the knotted proteins increases, particularly for more deeply threaded termini, the prevalence of traps in the free energy landscape also increases. Thus, in the case of longer knotted and slipknotted proteins, the folding mechanism is probably supported by chaperones. Overall, results imply that knotted proteins can be folded efficiently and survive evolutionary pressure in order to perform their biological functions.
Collapse
|
42
|
Abstract
This work explores the impact of knots, knot depth and motif of the threading terminus in protein folding properties (kinetics, thermodynamics and mechanism) via extensive Monte Carlo simulations of lattice models. A knotted backbone has no effect on protein thermodynamic stability but it may affect key aspects of folding kinetics. In this regard, we found clear evidence for a functional advantage of knots: knots enhance kinetic stability because a knotted protein unfolds at a distinctively slower rate than its unknotted counterpart. However, an increase in knot deepness does not necessarily lead to more effective changes in folding properties. In this regard, a terminus with a non-trivial conformation (e.g. hairpin) can have a more dramatic effect in enhancing kinetic stability than knot depth. Nevertheless, our results suggest that the probability of the denatured ensemble to keep knotted is higher for proteins with deeper knots, indicating that knot depth plays a role in determining the topology of the denatured state. Refolding simulations starting from denatured knotted conformations show that not every knot is able to nucleate folding and further indicate that the formation of the knotting loop is a key event in the folding of knotted trefoils. They also show that there are specific native contacts within the knotted core that are crucial to keep a native knotting loop in denatured conformations which otherwise have no detectable structure. The study of the knotting mechanism reveals that the threading of the knotting loop generally occurs towards late folding in conformations that exhibit a significant degree of structural consolidation.
Collapse
Affiliation(s)
- Miguel A. Soler
- Centro de Física da Matéria Condensada, Universidade de Lisboa, Lisboa, Portugal
- Departamento de Física, Universidade de Lisboa, Lisboa, Portugal
| | - Patrícia F. N. Faísca
- Centro de Física da Matéria Condensada, Universidade de Lisboa, Lisboa, Portugal
- Departamento de Física, Universidade de Lisboa, Lisboa, Portugal
| |
Collapse
|
43
|
Contessoto VG, Lima DT, Oliveira RJ, Bruni AT, Chahine J, Leite VBP. Analyzing the effect of homogeneous frustration in protein folding. Proteins 2013; 81:1727-37. [PMID: 23609962 DOI: 10.1002/prot.24309] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2012] [Revised: 03/13/2013] [Accepted: 03/18/2013] [Indexed: 11/09/2022]
Abstract
The energy landscape theory has been an invaluable theoretical framework in the understanding of biological processes such as protein folding, oligomerization, and functional transitions. According to the theory, the energy landscape of protein folding is funneled toward the native state, a conformational state that is consistent with the principle of minimal frustration. It has been accepted that real proteins are selected through natural evolution, satisfying the minimum frustration criterion. However, there is evidence that a low degree of frustration accelerates folding. We examined the interplay between topological and energetic protein frustration. We employed a Cα structure-based model for simulations with a controlled nonspecific energetic frustration added to the potential energy function. Thermodynamics and kinetics of a group of 19 proteins are completely characterized as a function of increasing level of energetic frustration. We observed two well-separated groups of proteins: one group where a little frustration enhances folding rates to an optimal value and another where any energetic frustration slows down folding. Protein energetic frustration regimes and their mechanisms are explained by the role of non-native contact interactions in different folding scenarios. These findings strongly correlate with the protein free-energy folding barrier and the absolute contact order parameters. These computational results are corroborated by principal component analysis and partial least square techniques. One simple theoretical model is proposed as a useful tool for experimentalists to predict the limits of improvements in real proteins.
Collapse
Affiliation(s)
- Vinícius G Contessoto
- Departamento de Física, Instituto de Biociências, Letras e Ciências Exatas, Universidade Estadual Paulista, Sao José do Rio Preto, São Paulo, 15054-000, Brazil
| | | | | | | | | | | |
Collapse
|
44
|
a Beccara S, Škrbić T, Covino R, Micheletti C, Faccioli P. Folding pathways of a knotted protein with a realistic atomistic force field. PLoS Comput Biol 2013; 9:e1003002. [PMID: 23555232 PMCID: PMC3605060 DOI: 10.1371/journal.pcbi.1003002] [Citation(s) in RCA: 66] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2012] [Accepted: 02/04/2013] [Indexed: 11/20/2022] Open
Abstract
We report on atomistic simulation of the folding of a natively-knotted protein, MJ0366, based on a realistic force field. To the best of our knowledge this is the first reported effort where a realistic force field is used to investigate the folding pathways of a protein with complex native topology. By using the dominant-reaction pathway scheme we collected about 30 successful folding trajectories for the 82-amino acid long trefoil-knotted protein. Despite the dissimilarity of their initial unfolded configuration, these trajectories reach the natively-knotted state through a remarkably similar succession of steps. In particular it is found that knotting occurs essentially through a threading mechanism, involving the passage of the C-terminal through an open region created by the formation of the native -sheet at an earlier stage. The dominance of the knotting by threading mechanism is not observed in MJ0366 folding simulations using simplified, native-centric models. This points to a previously underappreciated role of concerted amino acid interactions, including non-native ones, in aiding the appropriate order of contact formation to achieve knotting. It has been recently observed that the native structure of proteins can contain knots. These are formed during the folding process and are tightened in a specific (i.e. native) location, along the poly-peptide chain. The existence of knots hence implies a high degree coordination of local and global conformational changes, during the folding reaction. In this work we investigate how the knot is formed and what are the dynamical mechanisms which drive the self-entanglement process. To this end, we report on the first atomistically detailed numerical simulation of the folding of a knotted protein, based on a realistic description of the inter-atomic forces. These simulations show that the knot is formed by following a specific sequence of contacts. The comparison of the findings with those based on simplified folding models suggest that the productive succession of contacts is aided by a concerted interplay of amino acid interactions, arguably including non-native ones.
Collapse
Affiliation(s)
| | | | - Roberto Covino
- Physics Department, University of Trento, Trento, Italy
- INFN, Gruppo Collegato di Trento, Trento, Italy
| | | | - Pietro Faccioli
- Physics Department, University of Trento, Trento, Italy
- INFN, Gruppo Collegato di Trento, Trento, Italy
- * E-mail:
| |
Collapse
|
45
|
Andrews BT, Capraro DT, Sulkowska JI, Onuchic JN, Jennings PA. Hysteresis as a Marker for Complex, Overlapping Landscapes in Proteins. J Phys Chem Lett 2013; 4:180-188. [PMID: 23525263 PMCID: PMC3601837 DOI: 10.1021/jz301893w] [Citation(s) in RCA: 45] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/21/2023]
Abstract
Topologically complex proteins fold by multiple routes as a result of hard-to-fold regions of the proteins. Oftentimes these regions are introduced into the protein scaffold for function and increase frustration in the otherwise smooth-funneled landscape. Interestingly, while functional regions add complexity to folding landscapes, they may also contribute to a unique behavior referred to as hysteresis. While hysteresis is predicted to be rare, it is observed in various proteins, including proteins containing a unique peptide cyclization to form a fluorescent chromophore as well as proteins containing a knotted topology in their native fold. Here, hysteresis is demonstrated to be a consequence of the decoupling of unfolding events from the isomerization or hula-twist of a chromophore in one protein and the untying of the knot in a second protein system. The question now is- can hysteresis be a marker for the interplay of landscapes where complex folding and functional regions overlap?
Collapse
Affiliation(s)
| | - Dominique T. Capraro
- Department of Chemistry and Biochemistry, University of California, San Diego, San Diego, CA
| | | | - José N. Onuchic
- Center for Theoretical Biological Physics, Rice University, 6100 Main Street, Houston TX 77005
| | - Patricia A. Jennings
- Department of Chemistry and Biochemistry, University of California, San Diego, San Diego, CA
- Department of Chemistry and Biochemistry, University of California at San Diego (UCSD), 9500 Gilman Drive, La Jolla, CA 92093-0375, USA. Tel: 858-534-6417, fax: 858-534-7042,
| |
Collapse
|
46
|
Nickson AA, Wensley BG, Clarke J. Take home lessons from studies of related proteins. Curr Opin Struct Biol 2012; 23:66-74. [PMID: 23265640 PMCID: PMC3578095 DOI: 10.1016/j.sbi.2012.11.009] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2012] [Revised: 11/26/2012] [Accepted: 11/27/2012] [Indexed: 11/30/2022]
Abstract
The 'Fold Approach' involves a detailed analysis of the folding of several topologically, structurally and/or evolutionarily related proteins. Such studies can reveal determinants of the folding mechanism beyond the gross topology, and can dissect the residues required for folding from those required for stability or function. While this approach has not yet matured to the point where we can predict the native conformation of any polypeptide chain in silico, it has been able to highlight, amongst others, the specific residues that are responsible for nucleation, pathway malleability, kinetic intermediates, chain knotting, internal friction and Paracelsus switches. Some of the most interesting discoveries have resulted from the attempt to explain differences between homologues.
Collapse
Affiliation(s)
- Adrian A Nickson
- Department of Chemistry, University of Cambridge, Lensfield Rd, Cambridge CB2 1EW, UK.
| | | | | |
Collapse
|
47
|
Soler MA, Faísca PFN. How difficult is it to fold a knotted protein? In silico insights from surface-tethered folding experiments. PLoS One 2012; 7:e52343. [PMID: 23284997 PMCID: PMC3527535 DOI: 10.1371/journal.pone.0052343] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2012] [Accepted: 11/14/2012] [Indexed: 11/25/2022] Open
Abstract
We explore the effect of surface tethering on the folding process of a lattice protein that contains a trefoil knot in its native structure via Monte Carlo simulations. We show that the outcome of the tethering experiment depends critically on which terminus is used to link the protein to a chemically inert plane. In particular, if surface tethering occurs at the bead that is closer to the knotted core the folding rate becomes exceedingly slow and the protein is not able to find the native structure in all the attempted folding trajectories. Such low folding efficiency is also apparent from the analysis of the probability of knot formation, p(knot), as a function of nativeness. Indeed, p(knot) increases abruptly from ∼0 to ∼1 only when the protein has more than 80% of its native contacts formed, showing that a highly compact conformation must undergo substantial structural re-arrangement in order to get effectively knotted. When the protein is surface tethered by the bead that is placed more far away from the knotted core p(knot) is higher than in the other folding setups (including folding in the bulk), especially if conformations are highly native-like. These results show that the mobility of the terminus closest to the knotted core is critical for successful folding of trefoil proteins, which, in turn, highlights the importance of a knotting mechanism that is based on a threading movement of this terminus through a knotting loop. The results reported here predict that if this movement is blocked, knotting occurs via an alternative mechanism, the so-called spindle mechanism, which is prone to misfolding. Our simulations show that in the three considered folding setups the formation of the knot is typically a late event in the folding process. We discuss the implications of our findings for co-translational folding of knotted trefoils.
Collapse
Affiliation(s)
- Miguel A. Soler
- Centro de Física da Matéria Condensada, Universidade de Lisboa, Lisboa, Portugal
- Departamento de Física, Universidade de Lisboa, Lisboa, Portugal
| | - Patrícia F. N. Faísca
- Centro de Física da Matéria Condensada, Universidade de Lisboa, Lisboa, Portugal
- Departamento de Física, Universidade de Lisboa, Lisboa, Portugal
| |
Collapse
|