1
|
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
|
2
|
Salicari L, Trovato A. Entangled Motifs in Membrane Protein Structures. Int J Mol Sci 2023; 24:ijms24119193. [PMID: 37298146 DOI: 10.3390/ijms24119193] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2023] [Revised: 05/18/2023] [Accepted: 05/20/2023] [Indexed: 06/12/2023] Open
Abstract
Entangled motifs are found in one-third of protein domain structures, a reference set that contains mostly globular proteins. Their properties suggest a connection with co-translational folding. Here, we wish to investigate the presence and properties of entangled motifs in membrane protein structures. From existing databases, we build a non-redundant data set of membrane protein domains, annotated with the monotopic/transmembrane and peripheral/integral labels. We evaluate the presence of entangled motifs using the Gaussian entanglement indicator. We find that entangled motifs appear in one-fifth of transmembrane and one-fourth of monotopic proteins. Surprisingly, the main features of the distribution of the values of the entanglement indicator are similar to the reference case of general proteins. The distribution is conserved across different organisms. Differences with respect to the reference set emerge when considering the chirality of entangled motifs. Although the same chirality bias is found for single-winding motifs in both membrane and reference proteins, the bias is reversed, strikingly, for double-winding motifs only in the reference set. We speculate that these observations can be rationalized in terms of the constraints exerted on the nascent chain by the co-translational bio-genesis machinery, which is different for membrane and globular proteins.
Collapse
Affiliation(s)
- Leonardo Salicari
- Department of Physics and Astronomy 'Galileo Galilei', University of Padova, Via Marzolo 8, 35031 Padova, PD, Italy
- National Institute of Nuclear Physics (INFN), Padova Section, Via Marzolo 8, 35131 Padova, PD, Italy
| | - Antonio Trovato
- Department of Physics and Astronomy 'Galileo Galilei', University of Padova, Via Marzolo 8, 35031 Padova, PD, Italy
- National Institute of Nuclear Physics (INFN), Padova Section, Via Marzolo 8, 35131 Padova, PD, Italy
| |
Collapse
|
3
|
Barkataki K, Panagiotou E. The Jones polynomial of collections of open curves in 3-space. Proc Math Phys Eng Sci 2022. [DOI: 10.1098/rspa.2022.0302] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022] Open
Abstract
Measuring the entanglement complexity of collections of open curves in 3-space has been an intractable, yet pressing mathematical problem, relevant to a plethora of physical systems, such as in polymers and biopolymers. In this manuscript, we give a novel definition of the Jones polynomial that generalizes the classic Jones polynomial to collections of open curves in 3-space. More precisely, first we provide a novel definition of the Jones polynomial of linkoids (open link diagrams) and show that this is a well-defined single variable polynomial that is a topological invariant, which, for link-type linkoids, coincides with that of the corresponding link. Using the framework introduced in (Panagiotou E, Kauffman L. 2020
Proc. R. Soc. A
476
, 20200124. ((
doi:10.1098/rspa.2020.0124
)), this enables us to define the Jones polynomial of collections of open and closed curves in 3-space. For collections of open curves in 3-space, the Jones polynomial has real coefficients and it is a continuous function of the curves’ coordinates. As the endpoints of the curves tend to coincide, the Jones polynomial tends to that of the resultant link. We demonstrate with numerical examples that the novel Jones polynomial enables us to characterize the topological/geometrical complexity of collections of open curves in 3-space for the first time.
Collapse
Affiliation(s)
- Kasturi Barkataki
- School of Mathematical and Statistical Sciences, Arizona State University, Tempe, AZ 85287-1804, USA
| | - Eleni Panagiotou
- Department of Mathematics and SimCenter, University of Tennessee at Chattanooga, 613 McCallie Avenue, Chattanooga,TN 37403, USA
| |
Collapse
|
4
|
Baldwin Q, Sumpter B, Panagiotou E. The Local Topological Free Energy of the SARS-CoV-2 Spike Protein. Polymers (Basel) 2022; 14:polym14153014. [PMID: 35893978 PMCID: PMC9332627 DOI: 10.3390/polym14153014] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2022] [Revised: 07/15/2022] [Accepted: 07/19/2022] [Indexed: 02/04/2023] Open
Abstract
The novel coronavirus SARS-CoV-2 infects human cells using a mechanism that involves binding and structural rearrangement of its Spike protein. Understanding protein rearrangement and identifying specific amino acids where mutations affect protein rearrangement has attracted much attention for drug development. In this manuscript, we use a mathematical method to characterize the local topology/geometry of the SARS-CoV-2 Spike protein backbone. Our results show that local conformational changes in the FP, HR1, and CH domains are associated with global conformational changes in the RBD domain. The SARS-CoV-2 variants analyzed in this manuscript (alpha, beta, gamma, delta Mink, G614, N501) show differences in the local conformations of the FP, HR1, and CH domains as well. Finally, most mutations of concern are either in or in the vicinity of high local topological free energy conformations, suggesting that high local topological free energy conformations could be targets for mutations with significant impact of protein function. Namely, the residues 484, 570, 614, 796, and 969, which are present in variants of concern and are targeted as important in protein function, are predicted as such from our model.
Collapse
Affiliation(s)
- Quenisha Baldwin
- Department of Biology, Tuskegee University, Tuskegee, AL 36088, USA;
| | - Bobby Sumpter
- Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA;
| | - Eleni Panagiotou
- Department of Mathematics and SimCenter, University of Tennessee at Chattanooga, Chattanooga, TN 37403, USA
- Correspondence: or
| |
Collapse
|
5
|
Sleiman JL, Burton RH, Caraglio M, Gutierrez Fosado YA, Michieletto D. Geometric Predictors of Knotted and Linked Arcs. ACS POLYMERS AU 2022; 2:341-350. [PMID: 36254317 PMCID: PMC9562465 DOI: 10.1021/acspolymersau.2c00021] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
Inspired by how certain proteins “sense”
knots and
entanglements in DNA molecules, here, we ask if local geometric features
that may be used as a readout of the underlying topology of generic
polymers exist. We perform molecular simulations of knotted and linked
semiflexible polymers and study four geometric measures to predict
topological entanglements: local curvature, local density, local 1D
writhe, and nonlocal 3D writhe. We discover that local curvature is
a poor predictor of entanglements. In contrast, segments with maximum
local density or writhe correlate as much as 90% of the time with
the shortest knotted and linked arcs. We find that this accuracy is
preserved across different knot types and also under significant spherical
confinement, which is known to delocalize essential crossings in knotted
polymers. We further discover that nonlocal 3D writhe is the best
geometric readout of the knot location. Finally, we discuss how these
geometric features may be used to computationally analyze entanglements
in generic polymer melts and gels.
Collapse
Affiliation(s)
- Joseph L. Sleiman
- School of Physics and Astronomy, University of Edinburgh, Peter Guthrie Tait Road, Edinburgh EH9 3FD, United Kingdom
| | - Robin H. Burton
- School of Physics and Astronomy, University of Edinburgh, Peter Guthrie Tait Road, Edinburgh EH9 3FD, United Kingdom
| | - Michele Caraglio
- Institut für Theoretische Physik, Universität Innsbruck, Technikerstraße 21A, A-6020 Innsbruck, Austria
| | - Yair Augusto Gutierrez Fosado
- School of Physics and Astronomy, University of Edinburgh, Peter Guthrie Tait Road, Edinburgh EH9 3FD, United Kingdom
| | - 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 Cancer, University of Edinburgh, Edinburgh EH4 2XU, United Kingdom
| |
Collapse
|
6
|
Panagiotou E, Kauffman LH. Vassiliev measures of complexity of open and closed curves in 3-space. Proc Math Phys Eng Sci 2021. [DOI: 10.1098/rspa.2021.0440] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
In this article, we define Vassiliev measures of complexity for open curves in 3-space. These are related to the coefficients of the enhanced Jones polynomial of open curves in 3-space. These Vassiliev measures are continuous functions of the curve coordinates; as the ends of the curve tend to coincide, they converge to the corresponding Vassiliev invariants of the resulting knot. We focus on the second Vassiliev measure from the enhanced Jones polynomial for closed and open curves in 3-space. For closed curves, this second Vassiliev measure can be computed by a Gauss code diagram and it has an integral formulation, the double alternating self-linking integral. The double alternating self-linking integral is a topological invariant of closed curves and a continuous function of the curve coordinates for open curves in 3-space. For polygonal curves, the double alternating self-linking integral obtains a simpler expression in terms of geometric probabilities.
Collapse
Affiliation(s)
- Eleni Panagiotou
- Department of Mathematics and SimCenter, University of Tennessee at Chattanooga, Chattanooga, TN 37403, USA
| | - Louis H. Kauffman
- Department of Mathematics, Statistics and Computer Science, University of Illinois at Chicago, Chicago, IL 60607, USA
- Department of Mechanics and Mathematics, Novosibirsk State University, Novosibirsk, Russia
| |
Collapse
|
7
|
Herschberg T, Carrillo JMY, Sumpter BG, Panagiotou E, Kumar R. Topological Effects Near Order–Disorder Transitions in Symmetric Diblock Copolymer Melts. Macromolecules 2021. [DOI: 10.1021/acs.macromol.1c00780] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Tom Herschberg
- Department of Computer Science and Engineering, University of Tennessee at Chattanooga, Chattanooga, Tennessee 37403, United States
| | - Jan-Michael Y. Carrillo
- Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Bobby G. Sumpter
- Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Eleni Panagiotou
- Department of Mathematics, University of Tennessee at Chattanooga, Chattanooga, Tennessee 37403, United States
| | - Rajeev Kumar
- Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| |
Collapse
|
8
|
Panagiotou E, Kauffman LH. Knot polynomials of open and closed curves. Proc Math Phys Eng Sci 2020; 476:20200124. [PMID: 32922152 DOI: 10.1098/rspa.2020.0124] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2020] [Accepted: 07/01/2020] [Indexed: 01/05/2023] Open
Abstract
In this manuscript, we introduce a method to measure entanglement of curves in 3-space that extends the notion of knot and link polynomials to open curves. We define the bracket polynomial of curves in 3-space and show that it has real coefficients and is a continuous function of the curve coordinates. This is used to define the Jones polynomial in a way that it is applicable to both open and closed curves in 3-space. For open curves, the Jones polynomial has real coefficients and it is a continuous function of the curve coordinates and as the endpoints of the curve tend to coincide, the Jones polynomial of the open curve tends to that of the resulting knot. For closed curves, it is a topological invariant, as the classical Jones polynomial. We show how these measures attain a simpler expression for polygonal curves and provide a finite form for their computation in the case of polygonal curves of 3 and 4 edges.
Collapse
Affiliation(s)
- Eleni Panagiotou
- Department of Mathematics and SimCenter, University of Tennessee at Chattanooga, Chattanooga, TN 37403, USA
| | - Louis H Kauffman
- Department of Mathematics, Statistics and Computer Science University of Illinois at Chicago, Chicago, IL 60607-7045, USA.,Department of Mechanics and Mathematics, Novosibirsk State University, Novosibirsk, Russia
| |
Collapse
|
9
|
Grønbæk C, Hamelryck T, Røgen P. GISA: using Gauss Integrals to identify rare conformations in protein structures. PeerJ 2020; 8:e9159. [PMID: 32566389 PMCID: PMC7293858 DOI: 10.7717/peerj.9159] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2019] [Accepted: 04/18/2020] [Indexed: 12/13/2022] Open
Abstract
The native structure of a protein is important for its function, and therefore methods for exploring protein structures have attracted much research. However, rather few methods are sensitive to topologic-geometric features, the examples being knots, slipknots, lassos, links, and pokes, and with each method aimed only for a specific set of such configurations. We here propose a general method which transforms a structure into a ”fingerprint of topological-geometric values” consisting in a series of real-valued descriptors from mathematical Knot Theory. The extent to which a structure contains unusual configurations can then be judged from this fingerprint. The method is not confined to a particular pre-defined topology or geometry (like a knot or a poke), and so, unlike existing methods, it is general. To achieve this our new algorithm, GISA, as a key novelty produces the descriptors, so called Gauss integrals, not only for the full chains of a protein but for all its sub-chains. This allows fingerprinting on any scale from local to global. The Gauss integrals are known to be effective descriptors of global protein folds. Applying GISA to sets of several thousand high resolution structures, we first show how the most basic Gauss integral, the writhe, enables swift identification of pre-defined geometries such as pokes and links. We then apply GISA with no restrictions on geometry, to show how it allows identifying rare conformations by finding rare invariant values only. In this unrestricted search, pokes and links are still found, but also knotted conformations, as well as more highly entangled configurations not previously described. Thus, an application of the basic scan method in GISA’s tool-box revealed 10 known cases of knots as the top positive writhe cases, while placing at the top of the negative writhe 14 cases in cis-trans isomerases sharing a spatial motif of little secondary structure content, which possibly has gone unnoticed. Possible general applications of GISA are fold classification and structural alignment based on local Gauss integrals. Others include finding errors in protein models and identifying unusual conformations that might be important for protein folding and function. By its broad potential, we believe that GISA will be of general benefit to the structural bioinformatics community. GISA is coded in C and comes as a command line tool. Source and compiled code for GISA plus read-me and examples are publicly available at GitHub (https://github.com).
Collapse
Affiliation(s)
- Christian Grønbæk
- Department of Biology, University of Copenhagen, Copenhagen, Denmark.,Current affiliation: Department of Biology, Novo Nordisk Foundation Center for Basic Metabolic Research, University of Copenhagen, Copenhagen, Denmark
| | - Thomas Hamelryck
- Department of Biology, Department of Computer Science, University of Copenhagen, Copenhagen, Denmark
| | - Peter Røgen
- DTU COMPUTE, Technical University of Denmark, Kgs. Lyngby, Denmark
| |
Collapse
|
10
|
Panyukov S. Theory of Flexible Polymer Networks: Elasticity and Heterogeneities. Polymers (Basel) 2020; 12:E767. [PMID: 32244601 PMCID: PMC7240557 DOI: 10.3390/polym12040767] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2020] [Revised: 03/19/2020] [Accepted: 03/21/2020] [Indexed: 11/16/2022] Open
Abstract
A review of the main elasticity models of flexible polymer networks is presented. Classical models of phantom networks suggest that the networks have a tree-like structure. The conformations of their strands are described by the model of a combined chain, which consists of the network strand and two virtual chains attached to its ends. The distribution of lengths of virtual chains in real polydisperse networks is calculated using the results of the presented replica model of polymer networks. This model describes actual networks having strongly overlapping and interconnected loops of finite sizes. The conformations of their strands are characterized by the generalized combined chain model. The model of a sliding tube is represented, which describes the general anisotropic deformations of an entangled network in the melt. I propose a generalization of this model to describe the crossover between the entangled and phantom regimes of a swollen network. The obtained dependence of the Mooney-Rivlin parameters C 1 and C 2 on the polymer volume fraction is in agreement with experiments. The main results of the theory of heterogeneities in polymer networks are also discussed.
Collapse
Affiliation(s)
- Sergey Panyukov
- P. N. Lebedev Physics Institute, Russian Academy of Sciences, Moscow 117924, Russia
| |
Collapse
|
11
|
Ahmad R, Paul S, Basu S. Characterization of entanglements in glassy polymeric ensembles using the Gaussian linking number. Phys Rev E 2020; 101:022503. [PMID: 32168721 DOI: 10.1103/physreve.101.022503] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2019] [Accepted: 01/21/2020] [Indexed: 11/06/2022]
Abstract
We propose a method for enumerating entanglements between long chained, linear polymers that is based on the Gaussian linking number. The linking number is calculated between closely approaching segments of the macromolecular chains. Topological features of an entanglement, i.e., the extent to which one open segment winds around another, are reflected by the linking number. We show that using this measure, we can track disentanglement events through a deformation history and gain insights into how large scale disentanglements lead to failure. Incorporating an additional step where the topological entanglements identified along each chain are optimally clustered using standard clustering algorithms, we can also obtain a measure of the average number of rheological constraints that exist along each chain in an ensemble. Comparisons with other methods of enumerating entanglements, especially the primitive path analysis, are also made. Our results indicate that the linking number between two entangled segments in the undeformed state is a good indicator of the strength of the entanglement. Also, disentanglements occurring overwhelmingly around chain ends are an important cause of failure when a triaxial stress state exists in the polymer.
Collapse
Affiliation(s)
- Rasool Ahmad
- Department of Mechanical Engineering, Indian Institute of Technology Kanpur, Kanpur 208016, Uttar Pradesh, India
| | - Saptarshi Paul
- Department of Mechanical Engineering, Indian Institute of Technology Kanpur, Kanpur 208016, Uttar Pradesh, India
| | - Sumit Basu
- Department of Mechanical Engineering, Indian Institute of Technology Kanpur, Kanpur 208016, Uttar Pradesh, India
| |
Collapse
|
12
|
Norbiato F, Seno F, Trovato A, Baiesi M. Folding Rate Optimization Promotes Frustrated Interactions in Entangled Protein Structures. Int J Mol Sci 2019; 21:ijms21010213. [PMID: 31892272 PMCID: PMC6981561 DOI: 10.3390/ijms21010213] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2019] [Revised: 12/16/2019] [Accepted: 12/25/2019] [Indexed: 12/15/2022] Open
Abstract
Many native structures of proteins accomodate complex topological motifs such as knots, lassos, and other geometrical entanglements. How proteins can fold quickly even in the presence of such topological obstacles is a debated question in structural biology. Recently, the hypothesis that energetic frustration might be a mechanism to avoid topological frustration has been put forward based on the empirical observation that loops involved in entanglements are stabilized by weak interactions between amino-acids at their extrema. To verify this idea, we use a toy lattice model for the folding of proteins into two almost identical structures, one entangled and one not. As expected, the folding time is longer when random sequences folds into the entangled structure. This holds also under an evolutionary pressure simulated by optimizing the folding time. It turns out that optmized protein sequences in the entangled structure are in fact characterized by frustrated interactions at the closures of entangled loops. This phenomenon is much less enhanced in the control case where the entanglement is not present. Our findings, which are in agreement with experimental observations, corroborate the idea that an evolutionary pressure shapes the folding funnel to avoid topological and kinetic traps.
Collapse
Affiliation(s)
- Federico Norbiato
- Department of Physics and Astronomy, University of Padova, Via Marzolo 8, I-35131 Padova, Italy; (F.N.); (F.S.); (A.T.)
| | - Flavio Seno
- Department of Physics and Astronomy, University of Padova, Via Marzolo 8, I-35131 Padova, Italy; (F.N.); (F.S.); (A.T.)
- INFN, Sezione di Padova, Via Marzolo 8, I-35131 Padova, Italy
| | - Antonio Trovato
- Department of Physics and Astronomy, University of Padova, Via Marzolo 8, I-35131 Padova, Italy; (F.N.); (F.S.); (A.T.)
- INFN, Sezione di Padova, Via Marzolo 8, I-35131 Padova, Italy
| | - Marco Baiesi
- Department of Physics and Astronomy, University of Padova, Via Marzolo 8, I-35131 Padova, Italy; (F.N.); (F.S.); (A.T.)
- INFN, Sezione di Padova, Via Marzolo 8, I-35131 Padova, Italy
- Correspondence:
| |
Collapse
|
13
|
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
|
14
|
Baiesi M, Orlandini E, Seno F, Trovato A. Sequence and structural patterns detected in entangled proteins reveal the importance of co-translational folding. Sci Rep 2019; 9:8426. [PMID: 31182755 PMCID: PMC6557820 DOI: 10.1038/s41598-019-44928-3] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2019] [Accepted: 05/23/2019] [Indexed: 11/09/2022] Open
Abstract
Proteins must fold quickly to acquire their biologically functional three-dimensional native structures. Hence, these are mainly stabilized by local contacts, while intricate topologies such as knots are rare. Here, we reveal the existence of specific patterns adopted by protein sequences and structures to deal with backbone self-entanglement. A large scale analysis of the Protein Data Bank shows that loops significantly intertwined with another chain portion are typically closed by weakly bound amino acids. Why is this energetic frustration maintained? A possible picture is that entangled loops are formed only toward the end of the folding process to avoid kinetic traps. Consistently, these loops are more frequently found to be wrapped around a portion of the chain on their N-terminal side, the one translated earlier at the ribosome. Finally, these motifs are less abundant in natural native states than in simulated protein-like structures, yet they appear in 32% of proteins, which in some cases display an amazingly complex intertwining.
Collapse
Affiliation(s)
- Marco Baiesi
- Department of Physics and Astronomy, University of Padova, Via Marzolo 8, I-35131, Padova, Italy
- INFN, Sezione di Padova, Via Marzolo 8, I-35131, Padova, Italy
| | - Enzo Orlandini
- Department of Physics and Astronomy, University of Padova, Via Marzolo 8, I-35131, Padova, Italy
- INFN, Sezione di Padova, Via Marzolo 8, I-35131, Padova, Italy
| | - Flavio Seno
- Department of Physics and Astronomy, University of Padova, Via Marzolo 8, I-35131, Padova, Italy.
- INFN, Sezione di Padova, Via Marzolo 8, I-35131, Padova, Italy.
| | - Antonio Trovato
- Department of Physics and Astronomy, University of Padova, Via Marzolo 8, I-35131, Padova, Italy
- INFN, Sezione di Padova, Via Marzolo 8, I-35131, Padova, Italy
| |
Collapse
|
15
|
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
|
16
|
Panagiotou E, Millett KC, Atzberger PJ. Topological Methods for Polymeric Materials: Characterizing the Relationship Between Polymer Entanglement and Viscoelasticity. Polymers (Basel) 2019; 11:E437. [PMID: 30960421 PMCID: PMC6473770 DOI: 10.3390/polym11030437] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2019] [Revised: 02/12/2019] [Accepted: 02/26/2019] [Indexed: 12/31/2022] Open
Abstract
We develop topological methods for characterizing the relationship between polymer chain entanglement and bulk viscoelastic responses. We introduce generalized Linking Number and Writhe characteristics that are applicable to open linear chains. We investigate the rheology of polymeric chains entangled into weaves with varying topologies and levels of chain density. To investigate viscoelastic responses, we perform non-equilibrium molecular simulations over a range of frequencies using sheared Lees⁻Edwards boundary conditions. We show how our topological characteristics can be used to capture key features of the polymer entanglements related to the viscoelastic responses. We find there is a linear relation over a significant range of frequencies between the mean absolute Writhe W r and the Loss Tangent tan ( δ ) . We also find an approximate inverse linear relationship between the mean absolute Periodic Linking Number L K P and the Loss Tangent tan ( δ ) . Our results show some of the ways topological methods can be used to characterize chain entanglements to better understand the origins of mechanical responses in polymeric materials.
Collapse
Affiliation(s)
- Eleni Panagiotou
- Department of Mathematics and SimCenter, University of Tennessee at Chattanooga, Chattanooga, TN 37403, USA.
| | - Kenneth C Millett
- Department of Mathematics, University of California Santa Barbara, Santa Barbara, CA 93106-3080, USA.
| | - Paul J Atzberger
- Department of Mathematics and Department of Mechanical Engineering, University of California Santa Barbara, Santa Barbara, CA 93106-3080, USA.
| |
Collapse
|
17
|
Dai L, Renner CB, Doyle PS. Origin of metastable knots in single flexible chains. PHYSICAL REVIEW LETTERS 2015; 114:037801. [PMID: 25659023 DOI: 10.1103/physrevlett.114.037801] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/28/2014] [Indexed: 06/04/2023]
Abstract
Recent theoretical progress has explained the physics of knotting of semiflexible polymers, yet knotting of flexible polymers is relatively unexplored. We herein develop a new theory for the size distribution of knots on a flexible polymer and the existence of metastable knots. We show the free energy of a flexible molecule in a tube can be mapped to quantitatively reproduce the free energy distribution of a knot on a flexible chain. The size distribution of knots on flexible chains is expected to be universal and might be observed at a macroscopic scale, such as a string of hard balls.
Collapse
Affiliation(s)
- Liang Dai
- BioSystems and Micromechanics IRG, Singapore-MIT Alliance for Research and Technology Centre, Singapore 117543, Singapore
| | - C Benjamin Renner
- Department of Chemical Engineering, Massachusetts Institute of Technology (MIT), Cambridge, Massachusetts 02139, USA
| | - Patrick S Doyle
- BioSystems and Micromechanics IRG, Singapore-MIT Alliance for Research and Technology Centre, Singapore 117543, Singapore and Department of Chemical Engineering, Massachusetts Institute of Technology (MIT), Cambridge, Massachusetts 02139, USA
| |
Collapse
|
18
|
Foteinopoulou K, Karayiannis NC, Laso M. Monte Carlo simulations of densely-packed athermal polymers in the bulk and under confinement. Chem Eng Sci 2015. [DOI: 10.1016/j.ces.2014.08.021] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
19
|
Panagiotou E, Kröger M. Pulling-force-induced elongation and alignment effects on entanglement and knotting characteristics of linear polymers in a melt. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 90:042602. [PMID: 25375516 DOI: 10.1103/physreve.90.042602] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/12/2014] [Indexed: 06/04/2023]
Abstract
We employ a primitive path (PP) algorithm and the Gauss linking integral to study the degree of entanglement and knotting characteristics of linear polymer model chains in a melt under the action of a constant pulling force applied to selected chain ends. Our results for the amount of entanglement, the linking number, the average crossing number, the writhe of the chains and their PPs and the writhe of the entanglement strands all suggest a different response at the length scale of entanglement strands than that of the chains themselves and of the corresponding PPs. Our findings indicate that the chains first stretch at the level of entanglement strands and next the PP (tube) gets oriented with the "flow." These two phases of the extension and alignment of the chains coincide with two phases related to the disentanglement of the chains. Soon after the onset of external force the PPs attain a more entangled conformation, and the number of nontrivially linked end-to-end closed chains increases. Next, the chains disentangle continuously to attain an almost unentangled conformation. Using the linking matrix of the chains in the melt, we furthermore show that these phases are accompanied by a different scaling of the homogeneity of the global entanglement in the system. The homogeneity of the end-to-end closed chains first increases to a maximum and then decreases slowly to a value characterizing a completely unlinked system.
Collapse
Affiliation(s)
- E Panagiotou
- Department of Mathematics, University of California, Santa Barbara, California 93106, USA
| | - M Kröger
- Polymer Physics, Department of Materials, ETH Zurich, CH-8093 Zurich, Switzerland
| |
Collapse
|