G. Vernizzi,P. Ribeca,H. Orland,A. Zee
Abstract:We consider the folding of a self-avoiding homopolymer on a lattice, with saturating hydrogen bond interactions. Our goal is to numerically evaluate the statistical distribution of the topological genus of pseudoknotted configurations. The genus has been recently proposed for classifying pseudoknots (and their topological complexity) in the context of RNA folding. We compare our results on the distribution of the genus of pseudoknots, with the theoretical predictions of an existing combinatorial model for an infinitely flexible and stretchable homopolymer. We thus obtain that steric and geometric constraints considerably limit the topological complexity of pseudoknotted configurations, as it occurs for instance in real RNA molecules. We also analyze the scaling properties at large homopolymer length, and the genus distributions above and below the critical temperature between the swollen phase and the compact-globule phase, both in two and three dimensions.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to evaluate the topological complexity of pseudoknots in homopolymers and study the influence of geometric and steric constraints on this complexity. Specifically, the author hopes to determine the topological type distribution of pseudoknots, especially the statistical distribution of its topological genus, through numerical simulation. These pseudoknots play an important role in RNA folding, so understanding their topological properties is of great significance for predicting the tertiary structure of RNA.
### Research Background and Problem Description
1. **Importance of RNA Pseudoknots**
- The tertiary structure of RNA molecules contains many important structural motifs, one of which is called pseudoknots. Pseudoknots are non - planar structures, that is, their disk diagrams contain intersecting arcs.
- Pseudoknots are found in almost all organisms and play important roles in regulation and function. However, due to their complex topological structures, existing RNA structure prediction algorithms usually ignore the existence of pseudoknots.
2. **Definition of Topological Genus**
- The topological genus \( g \) is an integer used to describe the complexity of pseudoknots. It is defined by the Euler characteristic \( \chi \) :
\[
\chi = 1 - 2g
\]
- where \( \chi \) is the Euler characteristic of the disk diagram. The genus \( g \) represents the number of handles that need to be added to the disk surface in order to make all arcs non - intersecting.
3. **Limitations of Existing Models**
- Previous research proposed a combinatorial model, assuming that homopolymers are infinitely flexible and stretchable, ignoring the actual three - dimensional conformation and geometric constraints. This model predicts that the average genus of homopolymers is close to the maximum value \( L/4 \), where \( L \) is the length of the homopolymer.
- However, real RNA molecules are not infinitely flexible, but are limited by factors such as bending rigidity and helical persistence length. Therefore, the actual complexity of pseudoknots should be lower than theoretically predicted.
### Research Objectives
The main objective of this paper is to study the influence of geometric and steric constraints on the topological complexity of pseudoknots, especially the change in genus distribution, through numerical simulation. Specifically:
- **Evaluate the Influence of Geometric and Steric Constraints**: By modeling self - avoiding random walks (SAW) on a cubic lattice and considering saturated hydrogen - bond interactions, study how these constraints affect the genus distribution of pseudoknots.
- **Compare Behaviors at Different Temperatures**: Analyze the change in genus distribution above and below the critical temperature \( T_\theta \).
- **Verify Theoretical Predictions**: Compare the numerical results with previous combinatorial models to verify whether geometric and steric constraints significantly reduce the complexity of pseudoknots.
Through these studies, the author hopes to better understand the topological properties of pseudoknots in real RNA molecules and provide a theoretical basis for improving RNA structure prediction algorithms.