Academia.eduAcademia.edu

Outline

Title
Abstract
Introduction
Results and Discussion
Analysis of a Mutation
Conclusions
Methods
References
All Topics
Physics
Thermodynamics

Visualization of protein folding funnels in lattice models

Profile image of Vbp LeiteVbp Leite

2014, PloS one

https://doi.org/10.1371/JOURNAL.PONE.0100861
visibility

…

description

9 pages

link

1 file

Abstract

Protein folding occurs in a very high dimensional phase space with an exponentially large number of states, and according to the energy landscape theory it exhibits a topology resembling a funnel. In this statistical approach, the folding mechanism is unveiled by describing the local minima in an effective one-dimensional representation. Other approaches based on potential energy landscapes address the hierarchical structure of local energy minima through disconnectivity graphs. In this paper, we introduce a metric to describe the distance between any two conformations, which also allows us to go beyond the one-dimensional representation and visualize the folding funnel in 2D and 3D. In this way it is possible to assess the folding process in detail, e.g., by identifying the connectivity between conformations and establishing the paths to reach the native state, in addition to regions where trapping may occur. Unlike the disconnectivity maps method, which is based on the kinetic con...

Related papers

Topography of funneled landscapes determines the thermodynamics and kinetics of protein folding
Vitor Barbanti Pereira Leite

Proceedings of the National Academy of Sciences, 2012

The energy landscape approach has played a fundamental role in advancing our understanding of protein folding. Here, we quantify protein folding energy landscapes by exploring the underlying density of states. We identify three quantities essential for characterizing landscape topography: the stabilizing energy gap between the native and nonnative ensembles δ E , the energetic roughness Δ E , and the scale of landscape measured by the entropy S . We show that the dimensionless ratio between the gap, roughness, and entropy of the system accurately predicts the thermodynamics, as well as the kinetics of folding. Large Λ implies that the energy gap (or landscape slope towards the native state) is dominant, leading to more funneled landscapes. We investigate the role of topological and energetic roughness for proteins of different sizes and for proteins of the same size, but with different structural topologies. The landscape topography ratio Λ is shown to be monotonically correlated wi...

View PDFchevron_right
Locating transition path region in the free energy landscape of protein folding
amar nath gupta

arXiv: Biological Physics, 2017

Protein folding processes are generally described statistically with the help of multidimensional free energy landscape, typically reduced to a 1-D free energy profile along good reaction co-ordinate. There are many physical parameters which are responsible for protein molecule to hop between the native and unfolded states. The transition path region across the barrier is the region corresponding to the minimum fluctuation. The extent to which this transition region can extend beyond the obvious 1/2 kBT region has been a question of interest for a long time. We propose a new method to locate this transition path region and to study its dependence on the asymmetry of the transition state for a given free energy landscape. We have performed Brownian dynamics simulations with Gaussian white noise and Monte-Carlo simulation by sampling ten thousand successful transitions across the barrier for three different energy landscape having fixed barrier height with asymmetry in their curvature...

View PDFchevron_right
Folding Funnels and Frustration in Off-Lattice Minimalist Protein Landscapes
Angel Garcia

Proceedings of The National Academy of Sciences, 1998

A full quantitative understanding of the protein folding problem is now becoming possible with the help of the energy landscape theory and the protein folding funnel concept. Good folding sequences have a landscape that resembles a rough funnel where the energy bias towards the native state is larger than its ruggedness. Such a landscape leads not only to fast folding and stable native conformations but, more importantly, to sequences that are robust to variations in the protein environment and to sequence mutations. In this paper, an off-lattice model of sequences that fold into a ␤-barrel native structure is used to describe a framework that can quantitatively distinguish good and bad folders. The two sequences analyzed have the same native structure, but one of them is minimally frustrated whereas the other one exhibits a high degree of frustration.

View PDFchevron_right
Network structure of protein folding pathways
Zoltán Toroczkai

2007

The classical approach to protein folding inspired by statistical mechanics avoids the high dimensional structure of the conformation space by using effective coordinates. Here we introduce a network approach to capture the statistical properties of the structure of conformation spaces. Conformations are represented as nodes of the network, while links are transitions via elementary rotations around a chemical bond. Self-avoidance of a polypeptide chain introduces degree correlations in the conformation network, which in turn lead to energy landscape correlations. Folding can be interpreted as a biased random walk on the conformation network. We show that the folding pathways along energy gradients organize themselves into scale free networks, thus explaining previous observations made via molecular dynamics simulations. We also show that these energy landscape correlations are essential for recovering the observed connectivity exponent, which belongs to a different universality class than that of random energy models. In addition, we predict that the exponent and therefore the structure of the folding network fundamentally changes at high temperatures, as verified by our simulations on the AK peptide.

View PDFchevron_right
Folding domain B of protein A on a dynamically partitioned free energy landscape
Erik D Nelson

Proceedings of the National Academy of Sciences, 2008

View PDFchevron_right
A topological order parameter for describing folding free energy landscapes of proteins
pham lan

The Journal of Chemical Physics, 2018

We studied the refolding free energy landscape of 26 proteins using the Go-like model. The distance between the denaturated state and the transition state, X F , was calculated using the Bell theory and the nonlinear Dudko-Hummer-Szabo theory, and its relation to the geometrical properties of the native state was considered in detail. We showed that none of the structural parameters, such as the contact order, protein length, and radius of cross section, correlate with X F for all classes of proteins. To overcome this problem, we have introduced the nematic order parameter P 02 , which describes the ordering of the structured elements of the native state. Due to its topologically global nature, P 02 is better than other structural parameters in describing the folding free energy landscape. In particular, P 02 displays a good correlation with X F extracted from the nonlinear theory for all three classes of proteins. Therefore, this parameter can be used to predict X F for any protein, if its native structure is known.

View PDFchevron_right
Analyzing energy landscapes for folding model proteins
Roy L Johnston

The Journal of Chemical Physics, 2006

A new benchmark 20-bead HP model protein sequence ͑on a square lattice͒, which has 17 distinct but degenerate global minimum ͑GM͒ energy structures, has been studied using a genetic algorithm ͑GA͒. The relative probabilities of finding particular GM conformations are determined and related to the theoretical probability of generating these structures using a recoil growth constructor operator. It is found that for longer successful GA runs, the GM probability distribution is generally very different from the constructor probability, as other GA operators have had time to overcome any initial bias in the originally generated population of structures. Structural and metric relationships ͑e.g., Hamming distances͒ between the 17 distinct GM are investigated and used, in conjunction with data on the connectivities of the GM and the pathways that link them, to explain the GM probability distributions obtained by the GA. A comparison is made of searches where the sequence is defined in the normal ͑forward͒ and reverse directions. The ease of finding mirror image solutions are also compared. Finally, this approach is applied to rationalize the ease or difficulty of finding the GM for a number of standard benchmark HP sequences on the square lattice. It is shown that the relative probabilities of finding particular members of a set of degenerate global minima depend critically on the topography of the energy landscape in the vicinity of the GM, the connections and distances between the GM, and the nature of the operators used in the chosen search method.

View PDFchevron_right

References (59)

  1. Anfinsen CB (1973) Principles that govern the folding of protein chains. Science (New York, NY) 181: 223-230.
  2. Levinthal C (1968) Are there pathways for protein folding? Extrait du Journal de Chimie Physique 65.
  3. Bryngelson JD, Onuchic JN, Socci ND, Wolynes PG (1995) Funnels, pathways, and the energy landscape of protein folding: A synthesis. Proteins: Structure, Function, and Bioinformatics 21: 167-195.
  4. Leopold PE, Montal M, Onuchic JN (1992) Protein folding funnels: a kinetic approach to the sequence-structure relationship. Proceedings of the National Academy of Sciences of the United States of America 89: 8721-8725.
  5. Thirumalai D, O'Brien EP, Morrison G, Hyeon C (2010) Theoretical perspectives on protein folding. Annual Review of Biophysics 39: 159-183.
  6. Dill KA, Ozkan SB, Shell MS, Weikl TR (2008) The protein folding problem. Annual Review of Biophysics 37: 289-316.
  7. Onuchic JN, Luthey-Schulten Z, Wolynes PG (1997) Theory of protein folding: the energy landscape perspective. Annual review of physical chemistry 48: 545- 600.
  8. Klimov DK, Thirumalai D (1998) Linking rates of folding in lattice models of proteins with underlying thermodynamic characteristics. The Journal of Chemical Physics 109: 4119-4125.
  9. Sabelko J, Ervin J, Gruebele M (1999) Observation of strange kinetics in protein folding. Proceedings of the National Academy of Sciences of the United States of America 96: 6031-6036.
  10. Nymeyer H, Garcı ´a AE, Onuchic JN (1998) Folding funnels and frustration in off-lattice minimalist protein landscapes. Proceedings of the National Academy of Sciences 95: 5921-5928.
  11. Onuchic JN, Nymeyer H, Garcı ´a AE, Chahine J, Socci ND (2000) The energy landscape theory of protein folding: insights into folding mechanisms and scenarios. Advances in protein chemistry 53: 87-152.
  12. Schuler B, Lipman EA, Eaton WA (2002) Probing the free-energy surface for protein folding with single-molecule uorescence spectroscopy. Nature 419: 743- 747.
  13. Lee CL, Stell G, Wang J (2003) First-passage time distribution and non- markovian diffusion dynamics of protein folding. The Journal of Chemical Physics 118: 959-968.
  14. Chavez LL, Onuchic JN, Clementi C (2004) Quantifying the roughness on the free energy landscape: entropic bottlenecks and protein folding rates. Journal of the American Chemical Society 126: 8426-8432.
  15. Wang J, Oliveira RJ, Chu X, Whitford PC, Chahine J, et al. (2012) Topography of funneled landscapes determines the thermodynamics and kinetics of protein folding. Proceedings of the National Academy of Sciences 109: 15763-15768.
  16. Zhuravlev PI, Papoian GA (2010) Protein functional landscapes, dynamics, allostery: a tortuous path towards a universal theoretical framework. Quarterly Reviews of Biophysics 43: 295-332.
  17. Potoyan DA, Papoian GA (2011) Energy landscape analyses of disordered histone tails reveal special organization of their conformational dynamics. Journal of the American Chemical Society 133: 7405-7415.
  18. Itzhaki LS, Otzen DE, Fersht AR (1995) The structure of the transition state for folding of chymotrypsin inhibitor 2 analysed by protein engineering methods: evidence for a nucleationcondensation mechanism for protein folding. Journal of molecular biology 254: 260-288.
  19. Clementi C, Nymeyer H, Onuchic JN (2000) Topological and energetic factors: what determines the structural details of the transition state ensemble and ''en- route'' intermediates for protein folding? an investigation for small globular proteins. Journal of molecular biology 298: 937-953.
  20. Wales DJ (2010) Energy landscapes: some new horizons. Current Opinion in Structural Biology 20: 3-10.
  21. Wales D (2003) Energy Landscapes: Applications to Clusters, Biomolecules and Glasses. Cambridge University Press.
  22. Shan Y, Arkhipov A, Kim ET, Pan AC, Shaw DE (2013) Transitions to catalytically inactive conformations in EGFR kinase. Proceedings of the National Academy of Sciences of the United States of America 110: 7270-7275.
  23. Dobson CM (2003) Protein folding and misfolding. Nature 426: 884-890.
  24. Reddy AS, Wang L, Singh S, Ling YL, Buchanan L, et al. (2010) Stable and metastable states of human amylin in solution. Biophysical Journal 99: 2208- 2216.
  25. Wales DJ (2012) Decoding the energy landscape: extracting structure, dynamics and thermodynamics. Philosophical transactions Series A, Mathematical, physical, and engineering sciences 370: 2877-2899.
  26. Wales DJ, Bogdan TV (2006) Potential energy and free energy landscapes. The Journal of Physical Chemistry B 110: 20765-20776.
  27. Becker OM, Karplus M (1997) The topology of multidimensional potential energy surfaces: Theory and application to peptide structure and kinetics. The Journal of Chemical Physics 106: 1495-1517.
  28. Becker OM (1997) Quantitative visualization of a macromolecular potential energy ''funnel''. Journal of Molecular Structure: THEOCHEM 398-399: 507- 516.
  29. Wales DJ, Miller MA, Walsh TR (1998) Archetypal energy landscapes. Nature 394: 758-760.
  30. Miller MA, Doye JPK, Wales DJ (1999) Structural relaxation in atomic clusters: Master equation dynamics. Physical Review E 60: 3701-3718.
  31. Doye J, Miller M, Wales D (1999) The double-funnel energy landscape of the 38-atom lennard-jones cluster. The Journal of Chemical Physics 110: 6896.
  32. Doye J, Miller M, Wales D (1999) Evolution of the potential energy surface with size for lennardjones clusters. The Journal of Chemical Physics 111: 8417.
  33. Wales DJ, Dewsbury PEJ (2004) Effect of salt bridges on the energy landscape of a model protein. The Journal of Chemical Physics 121: 10284-10290.
  34. Miller MA,Wales DJ (1999) Energy landscape of a model protein. The Journal of Chemical Physics 111: 6610-6616.
  35. Evans DA, Wales DJ (2003) The free energy landscape and dynamics of met- enkephalin. The Journal of Chemical Physics 119: 9947-9955.
  36. Krivov SV, Karplus M (2002) Free energy disconnectivity graphs: Application to peptide models. The Journal of Chemical Physics 117: 10894-10903.
  37. Noe ´F, Fischer S (2008) Transition networks for modeling the kinetics of conformational change in macromolecules. Current Opinion in Structural Biology 18: 154-162.
  38. Prada-Gracia D, Go ´mez-Garden ˜es J, Echenique P, Falo F (2009) Exploring the free energy landscape: From dynamics to networks and back. PLoS Comput Biol 5: e1000415.
  39. Noe ´F, Horenko I, Schutte C, Smith J (2007) Hierarchical analysis of conformational dynamics in biomolecules: Transition networks of metastable states. The Journal of Chemical Physics 126.
  40. Rao F, Caisch A (2004) The protein folding network. Journal of molecular biology 342: 299-306.
  41. Dickson A, Brooks CL (2013) Native states of fast-folding proteins are kinetic traps. Journal of the American Chemical Society 135: 4729-4734.
  42. Socci ND, Onuchic JN (1995) Kinetic and thermodynamic analysis of proteinlike heteropolymers: Monte carlo histogram technique. The Journal of Chemical Physics 103: 4732-4744.
  43. Socci ND, Onuchic JN, Wolynes PG (1998) Protein folding mechanisms and the multidimensional folding funnel. Proteins 32: 136-158.
  44. Garstecki P, Hoang TX, Cieplak M (1999) Energy landscapes, supergraphs, and folding funnel in spin systems. Physical Review E 60: 3219-3226.
  45. Tejada E, Minghim R, Nonato LG (2003) On improved projection techniques to support visual exploration of multi-dimensional data sets. Information Visualization 2: 218-231.
  46. Choi S, Cha S, Tappert C (2010) A survey of binary similarity and distance measures. Journal on Systemics, Cybernetics and Informatics 8: 43-48.
  47. Tan PN, Steinbach M, Kumar V (2005) Introduction to data mining. Boston: Pearson Addison Wesley.
  48. Wales DJ (2006) Energy landscapes: calculating pathways and rates. Interna- tional Reviews in Physical Chemistry 25: 237-282.
  49. Wales DJ (2002) Discrete path sampling. Molecular Physics 100: 3285-3305.
  50. Wales DJ (2004) Some further applications of discrete path sampling to cluster isomerization. Molecular Physics 102: 891-908.
  51. Zheng W, Rohrdanz MA, Clementi C (2013) Rapid exploration of configuration space with diffusion-map-directed molecular dynamics. The journal of physical chemistry B.
  52. Laio A, Parrinello M (2002) Escaping free-energy minima. Proceedings of the National Academy of Sciences 99: 12562-12566.
  53. Otzen DE, Oliveberg M (1999) Salt-induced detour through compact regions of the protein folding landscape. Proceedings of the National Academy of Sciences 96: 11746-11751.
  54. Kurnik M, Hedberg L, Danielsson J, Oliveberg M (2012) Folding without charges. Proceedings of the National Academy of Sciences 109: 5705-5710.
  55. Chahine J, Nymeyer H, Leite VBP, Socci ND, Onuchic JN (2002) Specific and nonspecific collapse in protein folding funnels. Physical review letters 88: 168101.
  56. Cox TF, Cox MAA (2010) Multidimensional Scaling, Second Edition. CRC Press.
  57. Ware C (2001) Designing with a 2 1/2d attitude. Information Design Journal 10: 2001.
  58. Dima RI, Banavar JR, Cieplak M, Maritan A (1999) Statistical mechanics of protein-like heteropolymers. Proceedings of the National Academy of Sciences 96: 4904-4907.
  59. Abkevich VI, Gutin AM, Shakhnovich EI (1994) Free energy landscape for protein folding kinetics: Intermediates, traps, and multiple pathways in theory and lattice model simulations. The Journal of Chemical Physics 101: 6052-6062.

Related papers

Cunning simplicity of protein folding landscapes
A. Finkelstein

Protein Engineering, Design and Selection, 2001

Funnel-like landscapes are widely used to visualize protein folding. It might seem that any funnel-like energy landscape helps to avoid the 'Levinthal paradox', i.e. to avoid sampling the impossibly large number of conformations for a folding protein. This cunning suggestion, reinforced by beautiful drawings of the energy funnels, stimulated some simple models of protein folding; one of them [D.J. Bicout and A. Szabo (2000) Protein Sci., 9, 452-465] is especially straightforward and instructive. A thorough analysis of this strict funnel model (which does not consider a nucleation of phase separation in the course of folding) shows that it cannot provide a simultaneous explanation for both major features observed for protein folding: (i) folding within non-astronomical time, and (ii) coexistence of the native and the unfolded states during the folding process. On the contrary, the nucleation mechanism of protein folding can account for both these major features simultaneously. Keywords: coexistence of the native and the unfolded phases/ folding nucleus/protein folding/rate of folding/thermodynamic mid-transition/transition state/two-state kinetics

View PDFchevron_right
Protein folding funnels: the nature of the transition state ensemble
Nicholas Socci, Zaida Luthey-schulten

Folding and Design, 1996

Background: Energy landscape theory predicts that the folding funnel for a small fast-folding ␣-helical protein will have a transition state halfway to the native state. Estimates of the position of the transition state along an appropriate reaction coordinate can be obtained from linear free energy relationships observed for folding and unfolding rate constants as a function of denaturant concentration. The experimental results of Huang and Oas for repressor, Fersht and collaborators for CI2, and Gray and collaborators for cytochrome c indicate a free energy barrier midway between the folded and unfolded regions. This barrier arises from an entropic bottleneck for the folding process. Results: In keeping with the experimental results, lattice simulations based on the folding funnel description show that the transition state is not just a single conformation, but rather an ensemble of a relatively large number of configurations that can be described by specific values of one or a few order parameters (e.g. the fraction of native contacts). Analysis of this transition state or bottleneck region from our lattice simulations and from atomistic models for small ␣-helical proteins by Boczko and Brooks indicates a broad distribution for native contact participation in the transition state ensemble centered around 50%. Importantly, however, the lattice-simulated transition state ensemble does include some particularly hot contacts, as seen in the experiments, which have been termed by others a folding nucleus. Conclusions: Linear free energy relations provide a crude spectroscopy of the transition state, allowing us to infer the values of a reaction coordinate based on the fraction of native contacts. This bottleneck may be thought of as a collection of delocalized nuclei where different native contacts will have different degrees of participation. The agreement between the experimental results and the theoretical predictions provides strong support for the landscape analysis.

View PDFchevron_right
Topological methods for exploring low-density states in biomolecular folding pathways
Gunnar Carlsson

The Journal of Chemical Physics, 2009

View PDFchevron_right
Visualization of the nature of protein folding by a study of a distance constraint approach in two-dimensional models
Hiroshi Wako

Biopolymers, 1982

A distance constraint approach is applied to two-dimensional models of proteins in order to visualize the nature of protein folding and to examine the relative roles of different ranges of interaction. Three different native structures (I, 11, and 111) are considered; they have two different kinds of residues, viz., hydrophobic and hydrophilic, and different sequences of these residues. We examine how the distance constraint approach functions in the prediction of protein folding when we know the sequence of the residues, the (fixed) bond lengths, the mean distances between residues i and i + 2, and i and i + 3, and the mean distances for hydrophobic-hydrophobic, hydrophobic-hydrophilic, and hydrophilic-hydrophilic contacts between residues i and i + j , where j Z 4. This approach involves optimization of an object function with respect to 98 variables and is not free of the multiple-minimum problem. The optimization is always terminated if the chain is entangled and/or the segments (residues) are packed too compactly to move. In order to escape from such situations and to take the excluded-volume effect into account, a Monte Carlo method is used after the optimization is trapped in local minima. Success in the prediction of folding is found to depend on the starting conformations and on the native conformations. Fair success is obtained in predicting the helix-like structure in protein I and the overall structure of protein 111, but not the P-like structures of proteins I and 11. Insofar as the prediction of the structure of protein I11 is reasonable, it appears that some sequences of residues produce greater constraints on their conformations than others, if one considers only the hydrophobic and hydrophilic nature of the residues. These results imply that, in the folding of real proteins in three dimensions, the competition for hydrophobic (and hydrophilic) residues for inside (outside) positions in the molecule probably constitutes a necessary but not a sufficient condition to form and stabilize the native structure. The failure to predict the structure of protein 11, and part of that of protein I, suggests that there are two types of long-range interactions. One (which we considered here) is nonspecific (i.e., is defined only in terms of contacts between residues of the same or different polarity) and acts at any stage of protein folding; the other (which we did not consider here) is a specific interaction between residues in pairs and contributes only when the residues in the specific pair take on the native conformation. Presumably, incorporation of such specific long-range interactions, together with the nonspecific ones, is necessary for successful protein folding, using the distance constraint approach.

View PDFchevron_right
Topological determinants of protein folding
Nikolay Dokholyan

Proceedings of the National Academy of Sciences, 2002

The folding of many small proteins is kinetically a two-state process that represents overcoming the major free-energy barrier. A kinetic characteristic of a conformation, its probability to descend to the native state domain in the amount of time that represents a small fraction of total folding time, has been introduced to determine to which side of the free-energy barrier a conformation belongs. However, which features make a protein conformation on the folding pathway become committed to rapidly descending to the native state has been a mystery. Using two small, well characterized proteins, CI2 and C-Src SH3, we show how topological properties of protein conformations determine their kinetic ability to fold. We use a macroscopic measure of the protein contact network topology, the average graph connectivity, by constructing graphs that are based on the geometry of protein conformations. We find that the average connectivity is higher for conformations with a high folding probability than for those with a high probability to unfold. Other macroscopic measures of protein structural and energetic properties such as radius of gyration, rms distance, solvent-accessible surface area, contact order, and potential energy fail to serve as predictors of the probability of a given conformation to fold.

View PDFchevron_right

Related topics

  • Engineering
  • Physics
  • Thermodynamics
  • Chemistry
  • Biology
  • Medicine
  • Protein Conformation

    • Cited by

    Intrinsically disordered proteins: Ensembles at the limits of Anfinsen's dogma
    Atish Mohanty

    Biophysics Reviews, 2022

    Intrinsically disordered proteins (IDPs) are proteins that lack rigid 3D structure. Hence, they are often misconceived to present a challenge to Anfinsen's dogma. However, IDPs exist as ensembles that sample a quasi-continuum of rapidly interconverting conformations and, as such, may represent proteins at the extreme limit of the Anfinsen postulate. IDPs play important biological roles and are key components of the cellular protein interaction network (PIN). Many IDPs can interconvert between disordered and ordered states as they bind to appropriate partners. Conformational dynamics of IDPs contribute to conformational noise in the cell. Thus, the dysregulation of IDPs contributes to increased noise and “promiscuous” interactions. This leads to PIN rewiring to output an appropriate response underscoring the critical role of IDPs in cellular decision making. Nonetheless, IDPs are not easily tractable experimentally. Furthermore, in the absence of a reference conformation, discern...

    View PDFchevron_right
    Exploring Energy Landscapes of Intrinsically Disordered Proteins: Insights into Functional Mechanisms
    Vitor Barbanti Pereira Leite

    2021

    Intrinsically disordered proteins (IDPs) lack a rigid 3D structure and populate a polymorphic ensemble of conformations. Because of the lack of a reference conformation, their energy landscape representation in terms of reaction coordinates presents a daunting challenge. Here, our newly developed Energy Landscape Visualization Method (ELViM), a reaction coordinate-free approach, shows its prime application to explore frustrated energy landscapes of an intrinsically disordered protein, Prostate-Associated Gene 4 (PAGE4). PAGE4 is a transcriptional coactivator that potentiates the oncogene c-Jun. Two kinases, namely HIPK1 and CLK2, phosphorylate PAGE4 generating variants phosphorylated at different serine/threonine residues (HIPK1-PAGE4 and CLK2-PAGE4, respectively) with opposing functions. While HIPK1-PAGE4 predominantly phosphorylates Thr51 and potentiates c-Jun, CLK2-PAGE4 hyper-phosphorylates PAGE4 and attenuates transactivation. To understand the underlying mechanisms of conforma...

    View PDFchevron_right
    Drift-diffusion (DrDiff) framework determines kinetics and thermodynamics of two-state folding trajectory and tunes diffusion models
    Vinicius Contessoto

    The Journal of Chemical Physics, 2019

    The stochastic drift-diffusion (DrDiff) theory is an approach used to characterize the dynamical properties of simulation data. With new features in transition times analyses, the framework characterized the thermodynamic free-energy profile [F(Q)], the folding time (τ f), and transition path time (τ TP) by determining the coordinate-dependent drift-velocity [v(Q)] and diffusion [D(Q)] coefficients from trajectory time traces. In order to explore the DrDiff approach and to tune it with two other methods (Bayesian analysis and fep1D algorithm), a numerical integration of the Langevin equation with known D(Q) and F(Q) was performed and the inputted coefficients were recovered with success by the diffusion models. DrDiff was also applied to investigate the prion protein (PrP) kinetics and thermodynamics by analyzing folding/unfolding simulations. The protein structure-based model, the well-known Gō-model, was employed in a coarse-grained Cα level to generate long constant-temperature time series. PrP was chosen due to recent experimental single-molecule studies in D and τ TP that stressed the importance and the difficulty of probing these quantities and the rare transition state events related to prion misfolding and aggregation. The PrP thermodynamic double-well F(Q) profile, the "X" shape of τ f (T), and the linear shape of τ TP (T) were predicted with v(Q) and D(Q) obtained by the DrDiff algorithm. With the advance of single-molecule techniques, the DrDiff framework might be a useful ally for determining kinetic and thermodynamic properties by analyzing time observables of biomolecular systems. The code is freely available at https://github.com/ronaldolab/DrDiff.

    View PDFchevron_right
    580 California St., Suite 400
    San Francisco, CA, 94104
    © 2025 Academia. All rights reserved