1
|
Languin-Cattoën O, Laborie E, Yurkova DO, Melchionna S, Derreumaux P, Belyaev AV, Sterpone F. Exposure of Von Willebrand Factor Cleavage Site in A1A2A3-Fragment under Extreme Hydrodynamic Shear. Polymers (Basel) 2021; 13:polym13223912. [PMID: 34833213 PMCID: PMC8625202 DOI: 10.3390/polym13223912] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2021] [Revised: 11/05/2021] [Accepted: 11/06/2021] [Indexed: 12/22/2022] Open
Abstract
Von Willebrand Factor (vWf) is a giant multimeric extracellular blood plasma involved in hemostasis. In this work we present multi-scale simulations of its three-domains fragment A1A2A3. These three domains are essential for the functional regulation of vWf. Namely the A2 domain hosts the site where the protease ADAMTS13 cleavages the multimeric vWf allowing for its length control that prevents thrombotic conditions. The exposure of the cleavage site follows the elongation/unfolding of the domain that is caused by an increased shear stress in blood. By deploying Lattice Boltzmann molecular dynamics simulations based on the OPEP coarse-grained model for proteins, we investigated at molecular level the unfolding of the A2 domain under the action of a perturbing shear flow. We described the structural steps of this unfolding that mainly concerns the β-strand structures of the domain, and we compared the process occurring under shear with that produced by the action of a directional pulling force, a typical condition of single molecule experiments. We observe, that under the action of shear flow, the competition among the elongational and rotational components of the fluid field leads to a complex behaviour of the domain, where elongated structures can be followed by partially collapsed melted globule structures with a very different degree of exposure of the cleavage site. Our simulations pose the base for the development of a multi-scale in-silico description of vWf dynamics and functionality in physiological conditions, including high resolution details for molecular relevant events, e.g., the binding to platelets and collagen during coagulation or thrombosis.
Collapse
Affiliation(s)
- Olivier Languin-Cattoën
- Laboratoire de Biochimie Théorique, CNRS, Université de Paris, UPR 9080, 13 rue Pierre et Marie Curie, F-75005 Paris, France; (O.L.-C.); (E.L.); (P.D.)
| | - Emeline Laborie
- Laboratoire de Biochimie Théorique, CNRS, Université de Paris, UPR 9080, 13 rue Pierre et Marie Curie, F-75005 Paris, France; (O.L.-C.); (E.L.); (P.D.)
| | - Daria O. Yurkova
- Faculty of Physics, Lomonosov Moscow State University, 119991 Moscow, Russia;
| | - Simone Melchionna
- Dipartimento di Fisica, Università Sapienza, P.le A. Moro 5, 00185 Rome, Italy;
| | - Philippe Derreumaux
- Laboratoire de Biochimie Théorique, CNRS, Université de Paris, UPR 9080, 13 rue Pierre et Marie Curie, F-75005 Paris, France; (O.L.-C.); (E.L.); (P.D.)
| | - Aleksey V. Belyaev
- Faculty of Physics, Lomonosov Moscow State University, 119991 Moscow, Russia;
- Correspondence: (A.V.B.); (F.S.)
| | - Fabio Sterpone
- Laboratoire de Biochimie Théorique, CNRS, Université de Paris, UPR 9080, 13 rue Pierre et Marie Curie, F-75005 Paris, France; (O.L.-C.); (E.L.); (P.D.)
- Correspondence: (A.V.B.); (F.S.)
| |
Collapse
|
2
|
Tapia-Rojo R, Mazo JJ, Falo F. Thermal versus mechanical unfolding in a model protein. J Chem Phys 2019; 151:185105. [PMID: 31731855 DOI: 10.1063/1.5126071] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Abstract
Force spectroscopy techniques are often used to learn about the free energy landscape of single biomolecules, typically by recovering free energy quantities that, extrapolated to zero force, are compared to those measured in bulk experiments. However, it is not always clear how the information obtained from a mechanically perturbed system can be related to the information obtained using other denaturants since tensioned molecules unfold and refold along a reaction coordinate imposed by the force, which is not likely to be meaningful in its absence. Here, we explore this dichotomy by investigating the unfolding landscape of a model protein, which is unfolded first mechanically through typical force spectroscopy-like protocols and next thermally. When unfolded by nonequilibrium force extension and constant force protocols, we recover a simple two-barrier landscape as the protein reaches the extended conformation through a metastable intermediate. Interestingly, folding-unfolding equilibrium simulations at low forces suggested a totally different scenario, where this metastable state plays little role in the unfolding mechanism, and the protein unfolds through two competing pathways [R. Tapia-Rojo et al., J. Chem. Phys. 141, 135102 (2014)]. Finally, we use Markov state models to describe the configurational space of the unperturbed protein close to the critical temperature. The thermal dynamics is well understood by a one-dimensional landscape along an appropriate reaction coordinate, however it is very different from the mechanical picture. In this sense, the results of our protein model for the mechanical and thermal descriptions provide incompatible views of the folding/unfolding landscape of the system, and the estimated quantities to zero force result are hard to interpret.
Collapse
Affiliation(s)
- Rafael Tapia-Rojo
- Departamento de Física de la Materia Condensada, Instituto de Biocomputación y Física de Sistemas Complejos, Universidad de Zaragoza, 50009 Zaragoza, Spain
| | - Juan J Mazo
- Departamento de Física de la Materia Condensada, Instituto de Ciencia de Materiales de Aragón, CSIC-Universidad de Zaragoza, 50009 Zaragoza, Spain
| | - Fernando Falo
- Departamento de Física de la Materia Condensada, Instituto de Biocomputación y Física de Sistemas Complejos, Universidad de Zaragoza, 50009 Zaragoza, Spain
| |
Collapse
|
3
|
Szymczak P. Periodic forces trigger knot untying during translocation of knotted proteins. Sci Rep 2016; 6:21702. [PMID: 26996878 PMCID: PMC4800218 DOI: 10.1038/srep21702] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2015] [Accepted: 01/29/2016] [Indexed: 11/19/2022] Open
Abstract
Proteins need to be unfolded when translocated through the pores in mitochondrial and other cellular membranes. Knotted proteins, however, might get stuck during this process, jamming the pore, since the diameter of the pore is smaller than the size of maximally tightened knot. The jamming probability dramatically increases as the magnitude of the driving force exceeds a critical value, Fc. In this numerical study, we show that for deep knots Fc lies below the force range over which molecular import motors operate, which suggest that in these cases the knots will tighten and block the pores. Next, we show how such topological traps might be prevented by using a pulling protocol of a repetitive, on-off character. Such a repetitive pulling is biologically relevant, since the mitochondrial import motor, like other molecular motors transforms chemical energy into directed motions via nucleotide-hydrolysis-mediated conformational changes, which are cyclic in character.
Collapse
Affiliation(s)
- Piotr Szymczak
- Institute of Theoretical Physics, Faculty of Physics, University of Warsaw, Warsaw, Poland
| |
Collapse
|
4
|
Wojciechowski M, Szymczak P, Carrión-Vázquez M, Cieplak M. Protein unfolding by biological unfoldases: insights from modeling. Biophys J 2015; 107:1661-8. [PMID: 25296319 DOI: 10.1016/j.bpj.2014.07.035] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2014] [Revised: 06/09/2014] [Accepted: 07/15/2014] [Indexed: 10/24/2022] Open
Abstract
The molecular determinants of the high efficiency of biological machines like unfoldases (e.g., the proteasome) are not well understood. We propose a model to study protein translocation into the chamber of biological unfoldases represented as a funnel. It is argued that translocation is a much faster way of unfolding a protein than end-to-end stretching, especially in a low-force regime, because it allows for a conformational freedom while concentrating local tension on consecutive regions of a protein chain and preventing refolding. This results in a serial unfolding of the protein structures dominated by unzipping. Thus, pulling against the unfoldase pore is an efficient catalyst of the unfolding reaction. We also show that the presence of the funnel makes the tension along the backbone of the substrate protein nonuniform even when the protein gets unfolded. Hence, the stalling force measured by single-molecule force spectroscopy techniques may be smaller than the traction force of the unfoldase motor.
Collapse
Affiliation(s)
| | - Piotr Szymczak
- Institute of Theoretical Physics, Faculty of Physics, University of Warsaw, Warsaw, Poland
| | - Mariano Carrión-Vázquez
- Instituto Cajal, Consejo Superior de Investigaciones Cientificas and Instituto Madrileño de Estudios Avanzados en Nanociencia, Madrid, Spain
| | - Marek Cieplak
- Institute of Physics, Polish Academy of Sciences, Warsaw, Poland.
| |
Collapse
|
5
|
Tapia-Rojo R, Arregui S, Mazo JJ, Falo F. Mechanical unfolding of a simple model protein goes beyond the reach of one-dimensional descriptions. J Chem Phys 2014; 141:135102. [PMID: 25296840 DOI: 10.1063/1.4896620] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We study the mechanical unfolding of a simple model protein. The Langevin dynamics results are analyzed using Markov-model methods which allow to describe completely the configurational space of the system. Using transition-path theory we also provide a quantitative description of the unfolding pathways followed by the system. Our study shows a complex dynamical scenario. In particular, we see that the usual one-dimensional picture: free-energy vs end-to-end distance representation, gives a misleading description of the process. Unfolding can occur following different pathways and configurations which seem to play a central role in one-dimensional pictures are not the intermediate states of the unfolding dynamics.
Collapse
Affiliation(s)
- R Tapia-Rojo
- Instituto de Biocomputación y Física de Sistemas Complejos and Departamento de Física de la Materia Condensada, Universidad de Zaragoza, 50009 Zaragoza, Spain
| | - S Arregui
- Instituto de Biocomputación y Física de Sistemas Complejos and Departamento de Física de la Materia Condensada, Universidad de Zaragoza, 50009 Zaragoza, Spain
| | - J J Mazo
- Instituto de Ciencia de Materiales de Aragón and Departamento de Física de la Materia Condensada, CSIC-Universidad de Zaragoza, 50009 Zaragoza, Spain
| | - F Falo
- Instituto de Biocomputación y Física de Sistemas Complejos and Departamento de Física de la Materia Condensada, Universidad de Zaragoza, 50009 Zaragoza, Spain
| |
Collapse
|
6
|
|
7
|
Mazo JJ, Fajardo OY, Zueco D. Thermal activation at moderate-to-high and high damping: Finite barrier effects and force spectroscopy. J Chem Phys 2013; 138:104105. [DOI: 10.1063/1.4793983] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
|
8
|
Wales DJ, Head-Gordon T. Evolution of the potential energy landscape with static pulling force for two model proteins. J Phys Chem B 2012; 116:8394-411. [PMID: 22432920 DOI: 10.1021/jp211806z] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
The energy landscape is analyzed for off-lattice bead models of protein L and protein G as a function of a static pulling force. Two different pairs of attachment points (pulling directions) are compared in each case, namely, residues 1/56 and 10/32. For the terminal residue pulling direction 1/56, the distinct global minimum structures are all extended, aside from the compact geometry that correlates with zero force. The helical turns finally disappear at the highest pulling forces considered. For the 10/32 pulling direction, the changes are more complicated, with a variety of competing arrangements for beads outside the region where the force is directly applied. These alternatives produce frustrated energy landscapes, with low-lying minima separated by high barriers. The calculated folding pathways in the absence of force are in good agreement with previous work. The N-terminal hairpin folds first for protein L and the C-terminal hairpin for protein G, which exhibits an intermediate. However, for a relatively low static force, where the global minimum retains its structure, the folding mechanisms change, sometimes dramatically, depending on the protein and the attachment points. The scaling relations predicted by catastrophe theory are found to hold in the limit of short path lengths.
Collapse
Affiliation(s)
- David J Wales
- University Chemical Laboratories, Lensfield Road, Cambridge CB2 1EW, UK.
| | | |
Collapse
|
9
|
Lichter S, Rafferty B, Flohr Z, Martini A. Protein high-force pulling simulations yield low-force results. PLoS One 2012; 7:e34781. [PMID: 22529933 PMCID: PMC3329509 DOI: 10.1371/journal.pone.0034781] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2011] [Accepted: 03/09/2012] [Indexed: 11/20/2022] Open
Abstract
All-atom explicit-solvent molecular dynamics simulations are used to pull with extremely large constant force (750–3000 pN) on three small proteins. The introduction of a nondimensional timescale permits direct comparison of unfolding across all forces. A crossover force of approximately 1100 pN divides unfolding dynamics into two regimes. At higher forces, residues sequentially unfold from the pulling end while maintaining the remainder of the protein force-free. Measurements of hydrodynamic viscous stresses are made easy by the high speeds of unfolding. Using an exact low-Reynolds-number scaling, these measurements can be extrapolated to provide, for the first time, an estimate of the hydrodynamic force on low-force unfolding. Below 1100 pN, but surprisingly still at extremely large applied force, intermediate states and cooperative unfoldings as seen at much lower forces are observed. The force-insensitive persistence of these structures indicates that decomposition into unfolded fragments requires a large fluctuation. This finding suggests how proteins are constructed to resist transient high force. The progression of helix and sheet unfolding is also found to be insensitive to force. The force-insensitivity of key aspects of unfolding opens the possibility that numerical simulations can be accelerated by high applied force while still maintaining critical features of unfolding.
Collapse
Affiliation(s)
- Seth Lichter
- Department of Mechanical Engineering, Northwestern University, Evanston, Illinois, United States of America.
| | | | | | | |
Collapse
|
10
|
Caraglio M, Imparato A, Pelizzola A. Direction-dependent mechanical unfolding and green fluorescent protein as a force sensor. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2011; 84:021918. [PMID: 21929030 DOI: 10.1103/physreve.84.021918] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/08/2011] [Revised: 07/05/2011] [Indexed: 05/31/2023]
Abstract
An Ising-like model of proteins is used to investigate the mechanical unfolding of the green fluorescent protein along different directions. When the protein is pulled from its ends, we recover the major and minor unfolding pathways observed in experiments. Upon varying the pulling direction, we find the correct order of magnitude and ranking of the unfolding forces. Exploiting the direction dependence of the unfolding force at equilibrium, we propose a force sensor whose luminescence depends on the applied force.
Collapse
Affiliation(s)
- M Caraglio
- Dipartimento di Fisica and CNISM, Politecnico di Torino, c. Duca degli Abruzzi 24, Torino, Italy.
| | | | | |
Collapse
|
11
|
Luccioli S, Imparato A, Lepri S, Piazza F, Torcini A. Discrete breathers in a realistic coarse-grained model of proteins. Phys Biol 2011; 8:046008. [PMID: 21670494 DOI: 10.1088/1478-3975/8/4/046008] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
We report the results of molecular dynamics simulations of an off-lattice protein model featuring a physical force-field and amino-acid sequence. We show that localized modes of nonlinear origin, discrete breathers (DBs), emerge naturally as continuations of a subset of high-frequency normal modes residing at specific sites dictated by the native fold. DBs are time-periodic, space-localized vibrational modes that exist generically in nonlinear discrete systems and are known for their resilience and ability to concentrate energy for long times. In the case of the small β-barrel structure that we consider, DB-mediated localization occurs on the turns connecting the strands. At high energies, DBs stabilize the structure by concentrating energy on a few sites, while their collapse marks the onset of large-amplitude fluctuations of the protein. Furthermore, we show how breathers develop as energy-accumulating centres following perturbations even at distant locations, thus mediating efficient and irreversible energy transfers. Remarkably, due to the presence of angular potentials, the breather induces a local static distortion of the native fold. Altogether, the combination of these two nonlinear effects may provide a ready means for remotely controlling local conformational changes in proteins.
Collapse
Affiliation(s)
- Stefano Luccioli
- CNR-Consiglio Nazionale delle Ricerche, Istituto dei Sistemi Complessi, via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy.
| | | | | | | | | |
Collapse
|
12
|
Cérou F, Guyader A, Lelièvre T, Pommier D. A multiple replica approach to simulate reactive trajectories. J Chem Phys 2011; 134:054108. [DOI: 10.1063/1.3518708] [Citation(s) in RCA: 46] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
13
|
Caraglio M, Imparato A, Pelizzola A. Pathways of mechanical unfolding of FnIII10: Low force intermediates. J Chem Phys 2010; 133:065101. [DOI: 10.1063/1.3464476] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
14
|
Matsunaga Y, Li CB, Komatsuzaki T. Cooperativity at different space and time scales in multiscale protein dynamics. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2010; 82:016213. [PMID: 20866713 DOI: 10.1103/physreve.82.016213] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/02/2010] [Indexed: 05/29/2023]
Abstract
A method proposed by Matsunaga [Phys. Rev. Lett. 99, 238103 (2007)] is applied to simple stochastic models and two model proteins composed of 46 amino beads with three different kinds of residues. The method, which is based on the combination of the principal component analysis and the finite size Lyapunov exponent, characterize the coarse-grained dynamics in different spatiotemporal hierarchies in protein dynamics. The application of the method to model proteins reveals that the low-indexed (large-variance) principal components carry less-divergent, regularized dynamics at the coarse-grained scales on a less-frustrated energy landscape, whereas this less-divergent nature is less pronounced for a protein model with a more frustrated energy landscape. It is also revealed that our technique can differentiate the collective motions on the projected principal component space inherent to the system and the apparent collective behavior which can appear even in high-dimensional stochastic systems.
Collapse
|