Accurate in silico predictions of modified RNA interactions to a prototypical RNA-binding protein with λ-dynamics

RNA-binding proteins shape biology through their widespread functions in RNA biochemistry. Their function requires the recognition of specific RNA motifs for targeted binding. These RNA binding elements can be composed of both unmodified and chemically modified RNAs, of which over 170 chemical modifications have been identified in biology. Unmodified RNA sequence preferences for RNA-binding proteins have been widely studied, with numerous methods available to identify their preferred sequence motifs. However, only a few techniques can detect preferred RNA modifications, and no current method can comprehensively screen the vast array of hundreds of natural RNA modifications. Prior work demonstrated that λ-dynamics is an accurate in silico method to predict RNA base binding preferences of an RNA-binding antibody. This work extends that effort by using λ-dynamics to predict unmodified and modified RNA binding preferences of human Pumilio, a prototypical RNA-binding protein. A library of RNA modifications was screened at eight nucleotide positions along the RNA to identify modifications predicted to affect Pumilio binding. Computed binding affinities were compared with experimental data to reveal high predictive accuracy. In silico force field accuracies were also evaluated between CHARMM36 and Amber RNA force fields to determine the best parameter set to use in binding calculations. This work demonstrates that λ-dynamics can predict RNA interactions to a bona fide RNA-binding protein without the requirements of chemical reagents or new methods to experimentally test binding at the bench. Advancing in silico methods like λ-dynamics will unlock new frontiers in understanding how RNA modifications shape RNA biochemistry.