A Comparative Study of Disordered and Ordered Protein Folding Dynamics Using Computational Simulation

Rickie Xian
DOI: https://doi.org/10.48550/arXiv.2210.04557
2022-10-10
Abstract:Folding protein dynamics has been an area of high interest for quite some time, especially given the increased focus on the field of Biophysics. Because folding dynamics occur on such short time scales, empirical techniques developed for more "static" protein events, such as X-ray crystallography, nuclear magnetic resonance, and green fluorescent protein (GFP) labelling, aren't as applicable. Instead, computational methods must often be used to simulate these short lived yet highly dynamic events. One such computational method that is proven to provide much valuable insight into protein folding dynamics is Molecular Dynamics Simulation (MD Simulation). This simulation method is both highly computationally demanding, yet highly accurate in its modelling of a proteins physical behaviour. Besides MD Simulation, simulations in general are quite applicable in the context of these protein events. For example, the simple Gillespie algorithm, a computational technique which can be executed on almost any personal computer, provides quite the robust view into protein dynamics given its computational simplicity. This paper will compare the results of two simulations, an MD simulation of a disordered, six-residue, carcinogenic protein fragment, and a Gillespie algorithm based simulation of an ordered folding protein: the mathematically identical nature of the Gillespie algorithm time series of the asymptotically stochastic hyperbolic tangent dynamics for the wild type predicting the exact behaviour of the carcinogenic protein system time series will show the computational power simulations provide for analyzing both disordered and ordered protein systems.
Soft Condensed Matter,Mesoscale and Nanoscale Physics,Biological Physics,Medical Physics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to compare the folding kinetics of disordered proteins and ordered proteins through computational simulation. Specifically, the research aims to explore the following issues: 1. **Kinetic differences between disordered and ordered proteins**: The paper simulates disordered hexapeptide fragments through molecular dynamics (MD) and simulates ordered folding proteins using the Gillespie algorithm to compare the dynamic behaviors of these two types of proteins during the folding process. This helps to understand the structural and functional differences between disordered proteins (such as intrinsically disordered proteins, IDPs) and ordered proteins. 2. **Effectiveness of computational methods**: The paper evaluates the effectiveness and applicability of different computational methods in simulating protein - folding kinetics. In particular, MD simulation is used to simulate the complex dynamics of disordered proteins due to its high computational requirements and high precision; while the Gillespie algorithm is used to simulate ordered proteins because it is computationally simple yet still can provide profound insights into protein dynamics. 3. **Mathematical description of kinetic behaviors**: The paper uses mathematical models (such as the master equation of the Gillespie algorithm and the hyperbolic tangent function) to describe the kinetic behaviors of protein folding and explores the applicability of these mathematical models in explaining the dynamic behaviors of disordered and ordered proteins. ### Main findings - **Similarity of RMSD graphs**: The RMSD graph of the disordered octapeptide system shows similar behavior to the graph of the concentration of the ordered folding protein over time, indicating that although there are significant structural differences between disordered and ordered proteins, they may exhibit similar patterns in some kinetic characteristics. - **Influence of kinetic parameters**: When the folding rate \(k_{\text{fold}}\) is increased while the unfolding rate \(k_{\text{unfold}}\) remains unchanged, the time required for the system to reach equilibrium decreases, and the fluctuation amplitude of the equilibrium part increases. When \(k_{\text{fold}}\) is much larger than \(k_{\text{unfold}}\), the kinetic behavior of the system is similar to that of disordered proteins. - **Characteristics of disordered proteins**: The radius of gyration and the number of hydrogen bonds in the disordered octapeptide system show a high degree of randomness and fluctuation during the simulation process, which reflects the structural characteristics of disordered proteins. ### Conclusions - By adjusting the folding rate \(k_{\text{fold}}\) and the unfolding rate \(k_{\text{unfold}}\), the two - state system of ordered proteins can exhibit kinetic behaviors similar to those of disordered proteins. - Disorder proteins and ordered proteins have similarities in some kinetic characteristics, indicating that the folding kinetics of proteins can be simplified into a two - state system for modeling and analysis. - The choice of computational method has an important impact on the simulation results. MD simulation is suitable for high - precision simulation of complex systems, while the Gillespie algorithm is suitable for rapid simulation of simple systems. Through these studies, the paper provides a new perspective for understanding protein - folding kinetics and provides valuable references for future experimental and theoretical research.