Fast free-energy-based neutral set size estimates for the RNA genotype-phenotype map

Area of Science:

• Computational Biology
• Evolutionary Biology
• Bioinformatics

Background:

• The genotype-phenotype (GP) map of RNA secondary structure connects RNA sequences to their structures.
• Large-scale properties of GP maps, like neutral set size, impact evolutionary trajectories.
• Efficient computation of neutral set sizes is necessary for their application in evolutionary studies.

Purpose of the Study:

• To develop a computationally efficient and accurate method for calculating RNA neutral set sizes.
• To investigate the factors influencing variations in neutral set sizes.
• To extend neutral set size calculations to a many-to-many framework.

Main Methods:

• A novel method based on free energy estimates for computing neutral set sizes.
• Comparison with existing sample-based methods for speed and accuracy.
• Generalization of calculations from a many-to-one to a many-to-many GP map framework.

Main Results:

• The proposed free energy-based method is significantly faster than sample-based approaches.
• The method provides insights into variations in neutral set sizes, linking them to structural features like stack counts.
• Neutral set sizes are consistent across many-to-one and many-to-many frameworks, indicating robustness.

Conclusions:

• A new, efficient computational method for RNA neutral set size calculation is presented.
• Structural properties, such as fewer stacks, are associated with larger neutral set sizes.
• The choice of genotype-phenotype map framework does not fundamentally alter which RNA structures exhibit the largest neutral set sizes.