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
|Authors (as registered SciPost users):||Regina Ricotta|
|Preprint Link:||scipost_202212_00043v2 (pdf)|
|Date submitted:||2023-05-12 05:18|
|Submitted by:||Ricotta, Regina|
|Submitted to:||SciPost Physics Proceedings|
|Proceedings issue:||34th International Colloquium on Group Theoretical Methods in Physics (GROUP2022)|
Recently, a mathematical method to solve the Fokker Plank equation (FPE) enabled the analysis of the protein folding kinetics, through the construction of the temporal evolution of the probability density. A symmetric tri-stable potential function was used to describe the unfolded and folded states of the protein as well as an intermediate state of the protein. In this paper, the main points of the methodology are reviewed, based on the algebraic Supersymmetric Quantum Mechanics (SQM) formalism, and new results on the kinetics of the evolution of the system characterized in terms of the diffusion parameter are presented.
For Journal SciPost Physics Proceedings: Publish
(status: Editorial decision fixed and (if required) accepted by authors)
Author comments upon resubmission
List of changes
1. The article presents a summary of a theoretical methodology already published (Reference ) and new results (Section 3). It is a mathematical model built to understand the protein folding process by solving the Fokker Plank equation. Through a known mapping between this equation and a Schrödinger-type equation, both equations have the same spectrum. Thus, the approximate spectrum is solved through the algebraic method of supersymmetric quantum mechanics (SQM) which allows obtaining the time dependent probability distribution, solution of the the Fokker Planck equation for a given free energy V(x). This allows the subsequent analysis of the (theoretical) kinetics of the protein folding. It is therefore an initial model for the protein folding process. Thus, in future work (in progress) the free energy with specific parameters and units obtained from the computational simulation will be treated and the results compared.
2. The Abstract and Introduction were rewritten in order to make clear the points raised by the reviewers. New references were added, References  and . We were careful to ilustrate the model in Section 2 with results different from those published in reference  (different free energy V(x)).
3. Questions concerning the choice of parameters in Equation (10) are explained below it, in Section 2.1 with the addition of the following: "The choice of parameters was made in order to have several symmetrical potentials with V(0)=0 and different lateral depths of the wells, as in Reference . Thus, only the variation of the parameter "a" in each tri-stable potential is enough to deal with the depth of the lateral minima."
4. Question concerning the choice of parameter Q is explained in Section 2.2 with the addition of the following: "To illustrate the model, we choose the potential V (x) with a = 3.90456 and fixed diffusion (arbitrary) constant, Q = 0.5, as in Figure 2. The value of Q is arbitrary but it has to be fixed in order to apply the SQM methodology to obtain the spectrum, as it can be seen from the superpotential W1 in equation (7). In Section 3 we vary the diffusion constant and evaluate its impact on the folding kinetics. It should be mentioned that in other works, , the quantities are given in units 1/Q."
5. The new results are presented in Section 3. The first paragraph was rewritten, to avoid confusion between the time t and the passage time (τ). "Using the methodology developed in  for the protein folding process, the characteristic passage time τ was evaluated for different values of the diffusion parameter Q for the free energy of Figure 2, V (x) = 3.90456x6 − 8.93851x4 + 5.42373x2. For each fixed value of Q in the interval 0.4 < Q < 20, the passage time τ for the evolution of the population to the right well was calculated, starting from the initial position x0 = xmin = −1.05282, as shown in Figure 6 (dotted line). Figure 6 also shows the best fitting for the results (solid line), given as function of 1/Q."
Submission & Refereeing History
You are currently on this page