Synchronized motion of interface residues for evaluating protein–RNA complex binding affinity: application to aptamer‐mediated inhibition of TDP‐43 aggregates

Investigating the binding between proteins and aptamers, such as peptides or RNA molecules, is of crucial importance both for understanding the molecular mechanisms that regulate cellular activities and for therapeutic applications in several pathologies. Here, a new computational procedure, employing mainly docking, clustering analysis, and molecular dynamics simulations, was designed to estimate the binding affinities between a protein and some RNA aptamers, through the investigation of the dynamical behavior of the predicted molecular complex. Using the state-of-the-art software catRAPID, we computationally designed a set of RNA aptamers interacting with the TAR DNA-binding protein 43 (TDP-43), a protein associated with several neurodegenerative diseases, including amyotrophic lateral sclerosis (ALS). We thus devised a computational protocol to predict the RNA–protein molecular complex, so that the structural and dynamical behavior of such a complex can be investigated through extensive molecular dynamics simulation. We hypothesized that the coordinated and synchronized motion of the protein-binding residues, when in contact with RNA molecule, is a critical requisite in order to have a stable binding. Indeed, we calculated the motion covariance exhibited by the interface residues during molecular dynamics simulation: we tested the results against experimental measurements of binding affinity (in this case, the dissociation constant) for six RNA molecules, resulting in a linear correlation of about 0.9. Our findings suggest that the synchronized movement of interface residues plays a pivotal role in ensuring the stability within RNA–protein complexes, moreover providing insights into the contribution of each interface residue. This promising pipeline could thus contribute to the design of RNA aptamers interacting with proteins.