A multiscale study to characterize SHAPE probe/ RNA molecules pre-reactive complex and its role in initiating the SHAPE reaction
Abstract
Ribonucleic acid (RNA) molecules play a crucial role in nearly every cellular process, with their function closely tied to their three-dimensional (3D) structure. As a result, determining the precise 3D structure of RNAs is essential to understand their biological functions. However, obtaining high-resolution 3D structures remains a significant challenge using traditional biophysical techniques. To address this, chemical probing methods such as SHAPE (Selective 2'-Hydroxyl Acylation analyzed by Primer Extension) have gained widespread popularity. SHAPE reactivities have been introduced in 2D RNA predictors as soft constraints on nucleotide base pairing although the acylation reaction involves the ribose moiety. Little is known about the physical chemistry behind this reaction, leaving several reactivities unexplained. In this context, our aim is to unveil the complex relationship between the local structure of RNA, its dynamics, and the SHAPE chemical reactivity. In this study, using a multiscale approach based on biased molecular dynamics simulations and quantum mechanics/molecular mechanics calculations on the well characterized GAAA RNA tetraloop, we provide new molecular insights on the pre-reactive complex and its binding mode. Our findings emphasize the critical role of the local environment in facilitating the recruitment and proper accommodation of the SHAPE probe. All-atom umbrella sampling (US) molecular dynamics (MD) simulations underscored the significance of the binding angle in forming the pre-reactive complex, which was thoroughly analyzed using geometrical descriptors. In US QM/MM simulations, we used different initial structures to demonstrate the impact of favorable binding on the acylation reaction by characterizing the first step. Moreover, we investigated the deprotonation of the hydroxyl and its role in the acylation reaction, highlighting the role of the local environment and the probe proximity in this process. Our results suggest that the formation of oxyanion prior to the binding is not required to initiate the chemical reaction.
Related
- Substrate-bound structures of yeast ALA synthase reveal protein dynamics initiating heme biosynthesis
- Coupling Raman Spectroscopy and Acoustic Microscopy to Characterize Glioblastoma Resistance
- Integrated non-targeted metabolomics, volatile profiling and targeted quantification characterize chemical changes in Jinhua jasmine tea during the scenting process
- A sensitive Interactome-based approach to characterize drug activity in living human cells
- Using Bayesian inference to characterize cell-size regulation across organisms
BibTeX
@article{hognon2026multiscale,
title = {{A multiscale study to characterize SHAPE probe/ RNA molecules pre-reactive complex and its role in initiating the SHAPE reaction}},
author = {Cécilia Hognon and Ameni Ben Abdeljaoued and Pierre Hardouin and Mélanie Etheve-Quelquejeu and Bruno Sargueil and Damien Laage and Elisa Frezza},
journal = {Journal of Chemical Information and Modeling},
year = {2026},
month = sep,
doi = {10.1021/acs.jcim.5c00812},
url = {https://hal.science/hal-05226372v2},
note = {CC BY 4.0},
}
RIS
TY - JOUR
TI - A multiscale study to characterize SHAPE probe/ RNA molecules pre-reactive complex and its role in initiating the SHAPE reaction
AU - Cécilia Hognon
AU - Ameni Ben Abdeljaoued
AU - Pierre Hardouin
AU - Mélanie Etheve-Quelquejeu
AU - Bruno Sargueil
AU - Damien Laage
AU - Elisa Frezza
PY - 2026
DA - 2026/09/16
JO - Journal of Chemical Information and Modeling
DO - 10.1021/acs.jcim.5c00812
UR - https://hal.science/hal-05226372v2
AB - Ribonucleic acid (RNA) molecules play a crucial role in nearly every cellular process, with their function closely tied to their three-dimensional (3D) structure. As a result, determining the precise 3D structure of RNAs is essential to understand their biological functions. However, obtaining high-resolution 3D structures remains a significant challenge using traditional biophysical techniques. To address this, chemical probing methods such as SHAPE (Selective 2'-Hydroxyl Acylation analyzed by Primer Extension) have gained widespread popularity. SHAPE reactivities have been introduced in 2D RNA predictors as soft constraints on nucleotide base pairing although the acylation reaction involves the ribose moiety. Little is known about the physical chemistry behind this reaction, leaving several reactivities unexplained. In this context, our aim is to unveil the complex relationship between the local structure of RNA, its dynamics, and the SHAPE chemical reactivity. In this study, using a multiscale approach based on biased molecular dynamics simulations and quantum mechanics/molecular mechanics calculations on the well characterized GAAA RNA tetraloop, we provide new molecular insights on the pre-reactive complex and its binding mode. Our findings emphasize the critical role of the local environment in facilitating the recruitment and proper accommodation of the SHAPE probe. All-atom umbrella sampling (US) molecular dynamics (MD) simulations underscored the significance of the binding angle in forming the pre-reactive complex, which was thoroughly analyzed using geometrical descriptors. In US QM/MM simulations, we used different initial structures to demonstrate the impact of favorable binding on the acylation reaction by characterizing the first step. Moreover, we investigated the deprotonation of the hydroxyl and its role in the acylation reaction, highlighting the role of the local environment and the probe proximity in this process. Our results suggest that the formation of oxyanion prior to the binding is not required to initiate the chemical reaction.
SN - 1549-9596
SN - 1549-960X
ER -