Post-transcriptional modifications are crucial for RNA
function
and can affect its structure and dynamics. Force-field-based classical
molecular dynamics simulations are a fundamental tool to characterize
biomolecular dynamics, and their application to RNA is flourishing.
Here, we show that the set of force-field parameters for N6-methyladenosine (m6A) developed for the commonly used
AMBER force field does not reproduce duplex denaturation experiments
and, specifically, cannot be used to describe both paired and unpaired
states. Then, we use reweighting techniques to derive new parameters
matching available experimental data. The resulting force field can
be used to properly describe paired and unpaired m6A in
both syn and anti conformation,
which thus opens the way to the use of molecular simulations to investigate
the effects of N6 methylations on RNA structural dynamics.