The diversity of RNA structural elements and their documented role in human diseases make RNA an attractive therapeutic target. However, progress in drug discovery and development has been hindered by challenges in the determination of high-resolution RNA structures and a limited understanding of the parameters that drive RNA recognition by small molecules, including a lack of validated quantitative structure-activity relationships (QSAR). Herein, we developed QSAR models that quantitatively predict both thermodynamic and kinetic-based binding parameters of small molecules and the HIV-1 TAR model RNA system. A set of small molecules bearing diverse scaffolds was screened against the HIV-1-TAR construct using surface plasmon resonance, which provided the binding kinetics and affinities. The data was then analyzed using multiple linear regression (MLR) combined with feature selection to afford robust models for binding of diverse RNA-targeted scaffolds. The predictivity of the model was validated on untested small molecules. The QSAR models presented herein represent the first application of validated and predictive 2D-QSAR using multiple scaffolds against an RNA target. We expect the workflow to be generally applicable to other RNA structures, ultimately providing essential insight into the small molecule descriptors that drive selective binding interactions and, consequently, providing a platform that can exponentially increase the efficiency of ligand design and optimization without the need for high-resolution RNA structures.
A comparative analysis of the triloops in all high-resolution RNA structures reveals sequence structure relationships