Embedding a Native State into a Random Heteropolymer Model: The Dynamic Approach

Zoran Konkoli,John Hertz

DOI: https://doi.org/10.1103/PhysRevE.67.051915

2003-01-24

Abstract:We study a random heteropolymer model with Langevin dynamics, in the supersymmetric formulation. Employing a procedure similar to one that has been used in static calculations, we construct an ensemble in which the affinity of the system for a native state is controlled by a "selection temperature" T0. In the limit of high T0, the model reduces to a random heteropolymer, while for T0-->0 the system is forced into the native state. Within the Gaussian variational approach that we employed previously for the random heteropolymer, we explore the phases of the system for large and small T0. For large T0, the system exhibits a (dynamical) spin glass phase, like that found for the random heteropolymer, below a temperature Tg. For small T0, we find an ordered phase, characterized by a nonzero overlap with the native state, below a temperature Tn \propto 1/T0 > Tg. However, the random-globule phase remains locally stable below Tn, down to the dynamical glass transition at Tg. Thus, in this model, folding is rapid for temperatures between Tg and Tn, but below Tg the system can get trapped in conformations uncorrelated with the native state. At a lower temperature, the ordered phase can also undergo a dynamical glass transition, splitting into substates separated by large barriers.

Statistical Mechanics

What problem does this paper attempt to address?

The core problem that this paper attempts to solve is to understand the kinetic behavior in the protein - folding process, especially to explore how to embed a native state in the random heteropolymer model and study the influence of this embedding on the phase behavior of the system. Specifically, the author constructs a random heteropolymer model containing the native state by introducing a "selection temperature" \(T_0\) to control the system's affinity for the native state. The following is the specific description of this problem:

1. **Background and Motivation**:
- Protein folding is a crucial process in biology that determines the function of proteins.
- The random heteropolymer model (RHP) is widely used to study the general characteristics of protein folding, but these models usually do not contain the native state.
- In order to understand why proteins do not exhibit glassy - state behavior, the author hopes to study a model that can both exhibit the glassy state and fold into the native state.

2. **Research Objectives**:
- Construct a random heteropolymer model containing the native state and control the system's affinity for the native state by adjusting the selection temperature \(T_0\).
- Study the phase behavior of the system at different temperatures, especially the phase transitions at high \(T_0\) and low \(T_0\).

3. **Methods**:
- Use Langevin dynamics and supersymmetric formulation to describe the kinetic behavior of the model.
- Introduce a selection temperature \(T_0\). When \(T_0\rightarrow\infty\), the model degenerates into a random heteropolymer; when \(T_0\rightarrow0\), the system is forced into the native state.
- Study the phase diagram of the system through the Gaussian variational approach (GVA), especially the kinetic behavior at high \(T_0\) and low \(T_0\).

4. **Main Findings**:
- For high \(T_0\), the system exhibits a dynamic spin - glass phase, similar to the random heteropolymer model.
- For low \(T_0\), the system enters an ordered phase below a certain temperature \(T_n\propto1/T_0\) and shows a non - zero overlap with the native state.
- At low temperatures, the ordered phase can further undergo a dynamic glass transition and split into multiple metastable states.

In summary, this paper aims to study the embedding of the native state in the random heteropolymer model and its influence on the phase behavior of the system by introducing the selection temperature \(T_0\), so as to better understand the kinetic behavior in the protein - folding process.