Dominant Pathways in Protein Folding

P.Faccioli,M.Sega,F.Pederiva,H.Orland
DOI: https://doi.org/10.1103/PhysRevLett.97.108101
2006-07-27
Abstract:We present a method to investigate the kinetics of protein folding on a long time-scale and the dynamics underlying the formation of secondary and tertiary structures during the entire reaction. The approach is based on the formal analogy between thermal and quantum diffusion: by writing the solution of the Fokker-Planck equation for the time-evolution of a protein in a viscous heat-bath in terms of a path integral, we derive a Hamilton-Jacobi variational principle from which we are able to compute the most probable pathway of folding. The method is applied to the folding of the Villin Headpiece Subdomain, in the framework of a Go-model. We have found that, in this model, the transition occurs through an initial collapsing phase driven by the starting coil configuration and a later rearrangement phase, in which secondary structures are formed and all computed paths display strong similarities. This method is completely general, does not require the prior knowledge of any reaction coordinate and represents an efficient tool to perfom ab-initio simulations of the entire folding process with available computers.
Biomolecules,Soft Condensed Matter,High Energy Physics - Phenomenology
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to study the dynamic mechanism of protein folding on a long - time scale, especially the dynamic process when forming secondary and tertiary structures during protein folding. Due to high computational complexity, traditional methods are difficult to simulate the time evolution of typical proteins (about 100 residues) for more than dozens of nanoseconds. Therefore, this paper proposes a new method. By making a formal analogy between Langevin diffusion and quantum propagation, a Hamilton - Jacobi variational principle is derived using the path integral method, so that the most probable folding path can be calculated. This method does not need to know any reaction coordinates in advance and can be used as an effective tool to simulate the entire folding process from scratch on existing computers. Specifically, the paper focuses on: 1. **Protein folding dynamics on a long - time scale**: Overcome the computational problems caused by the huge difference in time scales in traditional methods. 2. **Folding path calculation without prior knowledge**: It does not depend on the selection of specific reaction coordinates and can calculate the most probable folding path. 3. **Efficient simulation of protein folding**: Through the Hamilton - Jacobi variational principle, the calculation is simplified, making it possible to simulate the entire folding process on existing computers. To verify the effectiveness of this method, the author applied it to the study of the folding transition of the 36 - mer Villin Headpiece sub - domain. The results show that this method can effectively describe the two main stages of protein folding: the initial collapse stage and the subsequent rearrangement stage, in which the latter forms the secondary structure. These results are qualitatively consistent with the experimental data, indicating that this method has certain practical value.