SciPost Submission Page
The impact of the diffusion parameter on the passage time of the folding process
by Marcelo Tozo Araujo, Jorge Chahine, Elso Drigo Filho and Regina Maria Ricotta
This is not the latest submitted version.
This Submission thread is now published as
Submission summary
Authors (as registered SciPost users): | Regina Ricotta |
Submission information | |
---|---|
Preprint Link: | scipost_202212_00043v1 (pdf) |
Date submitted: | 2022-12-18 04:27 |
Submitted by: | Ricotta, Regina |
Submitted to: | SciPost Physics Proceedings |
Proceedings issue: | 34th International Colloquium on Group Theoretical Methods in Physics (GROUP2022) |
Ontological classification | |
---|---|
Academic field: | Physics |
Specialties: |
|
Approach: | Theoretical |
Abstract
This work shows a mathematical method to solve the diffusion equation, which enabled the analysis of the protein folding kinetics, through the construction of the temporal evolution of the probability density inspired by the funnel concept of energy landscape. A symmetric tri-stable potential function is used to describe the unfolded and folded states of the protein as well as a set of intermediate states of the protein. The kinetics of the evolution of the system is characterized in terms of the diffusion parameter.
Current status:
Reports on this Submission
Report #2 by Anonymous (Referee 1) on 2023-4-14 (Invited Report)
- Cite as: Anonymous, Report on arXiv:scipost_202212_00043v1, delivered 2023-04-14, doi: 10.21468/SciPost.Report.7047
Strengths
The mathematical formalism and biophysical approach
Weaknesses
lack of discussion, and references
Report
The Authors have applied the Fokker-Plank equation to the diffusion model during the protein folding process using a random potential V that represents 3 protein energy state profiles: unfolded, intermediate/transition, and folded, according to my understanding. The mathematical formalism and biophysical approach in this work is very sophisticated and highlighted the need of research in this area.
However, a few points should be addressed:
Introduction: There is not many references for the topic, especially around the protein folding issue and its deregulation and the consequences in the medical field. The author could explore this area in more detail and show data for some diseases and how this work can help to address those biological questions.
Methodology and Formalism:
In general there is a lack of physics parameters units, such as for V, Q, t and f.
On the equation 10, where do the constants -8.93851 and +5.42373 come from? Are they arbitrary or coming from a simulation? Same for the Figure 2 where “a” was fixed as 3.90456. It should be shown how those numbers were fixed, or mentioned if it was arbitrary.
Is “a” and “” (alpha) the same variable? Could that be a typo?
What is the definition of N(t) shown in the equation 11, is it the diffusion process and its distribution? A better variable definition is required.
t and t’ were mentioned on the section 2.2 and defined only in the section 2.3 as passing time and then the same variable is defined as passage time. Please define the variable and its units before it appears in the equations, and keep it consistent throughout the paper.
All the figures (graphics) need to have a bigger x and y variables allocated to its axis for better visualisation. Such as “t x Q” or “t x Xo” or “Population x t”.
In this work there is a lack of discussion, for example the authors could explore how this model could address the Levinthal’s protein folding paradox for (B Bagchi, 1992). Also, there is room to discuss how this Fokker-Plank equation to diffusion model can impact biomedical or biophysics research or techniques, as a novelty proposed in the introduction.
Overall, this is good paper and addresses a part of biophysics research where research is lacking and could bring significant breakthrough in science and medical research. I would recommend it for publication after the authors checked the points raised above.
Requested changes
All the figures (graphics) need to have a bigger x and y variables allocated to its axis, for a better visualisation. Such as “t x Q” or “t x Xo” or “Population x t”.
Report #1 by Anonymous (Referee 2) on 2023-3-22 (Invited Report)
- Cite as: Anonymous, Report on arXiv:scipost_202212_00043v1, delivered 2023-03-22, doi: 10.21468/SciPost.Report.6944
Report
In this work, the authors present an application of the Fokker-Planck equation to model the diffusion process during protein folding. The authors use a potential V that mimics a protein free energy profile with 3 states (unfolded, intermediate, and folded). Nonetheless, there is a lack of discussion in associating this work with protein folding. Comparison with experimental results should be performed. The language and writing should be revised. Breaking sentences would help the reader. There are also major scientific points that should be discussed:
There is a lack of references, especially discussing seminal work from Dmitrii E. Makarov, Robert B. Best, Gerhard Hummer, and others.
There is also a lack of references about the protein folding problem, particularly in the introduction.
Section 2.1 - How do the authors justify the choice of parameters in equation 10?
Section 2.2:
The same, Is there any argumentation for using those values of a and Q?
Why choose a potential V(x) where the intermediate is more stable than the native state? I recommend doing the same calculations for a scenario where the protein's native state is more stable than the intermediate.
What happens when d(V(x)/dx = 0? In this case, the force is zero. How does the protein decide to go right or left in the potential? Do authors consider the case of half o the time, when reaching this inflection point, the protein returns to the previous state?
Why the potential in the native state goes up? What exists after the native state?
The authors mentioned that the reaction coordinate x is the number of native contacts (page 2). How is possible to have negative values of x? The minimum value must be 0. A negative number of contacts is inconsistent.
How these calculations are different from a partible diffusing under the influence of the potential V(x)?
Section 3:
There is a lack of physical units in all calculations and figures. What are the units of V, f, Q, t?
What are the criteria to consider the protein folded when calculating the passage time?
Section 4:
There is a lack of discussion, and the conclusion seems like a textbook exercise. I recommend exploring different parameters and connecting the results with experiments or even simulations.
Anonymous on 2023-03-22 [id 3502]
In this work, the authors present an application of the Fokker-Planck equation to model the diffusion process during protein folding. The authors use a potential V that mimics a protein free energy profile with 3 states (unfolded, intermediate, and folded). Nonetheless, there is a lack of discussion in associating this work with protein folding. Comparison with experimental results should be performed. The language and writing should be revised. Breaking sentences would help the reader. There are also major scientific points that should be discussed:
There is a lack of references, especially discussing seminal work from Dmitrii E. Makarov, Robert B. Best, Gerhard Hummer, and others.
There is also a lack of references about the protein folding problem, particularly in the introduction.
Section 2.1 - How do the authors justify the choice of parameters in equation 10?
Section 2.2:
The same, Is there any argumentation for using those values of a and Q?
Why choose a potential V(x) where the intermediate is more stable than the native state? I recommend doing the same calculations for a scenario where the protein's native state is more stable than the intermediate.
What happens when d(V(x)/dx = 0? In this case, the force is zero. How does the protein decide to go right or left in the potential? Do authors consider the case of half o the time, when reaching this inflection point, the protein returns to the previous state?
Why the potential in the native state goes up? What exists after the native state?
The authors mentioned that the reaction coordinate x is the number of native contacts (page 2). How is possible to have negative values of x? The minimum value must be 0. A negative number of contacts is inconsistent.
How these calculations are different from a partible diffusing under the influence of the potential V(x)?
Section 3:
There is a lack of physical units in all calculations and figures. What are the units of V, f, Q, t?
What are the criteria to consider the protein folded when calculating the passage time?
Section 4:
There is a lack of discussion, and the conclusion seems like a textbook exercise. I recommend exploring different parameters and connecting the results with experiments or even simulations.
Anonymous on 2023-03-10 [id 3469]
This work was presented in the context of the 34th International Colloquium on Group Theoretical Methods in Physics Strasbourg, 18-22 July 2022, according to the authors.
My comments are as follows: * There is neither mention or relevance to group theory * The connection to protein folding is vague, at least much less than the work quoted in reference: Reference 2 by Wolynes. * The paper is very similar to Réference [1] by the same authors. * The scientific level is weak, a simple exercise based on the diffusion equation (1), no new finding.
The paper should not be accepted in SciPost.