Abstract:Alchemical free energy calculations are widely used to predict the binding affinity of small molecule ligands to protein targets; however, the application of these methods to RNA targets has not been deeply explored. We systematically investigated how modeling decisions affect the performance of absolute binding free energy calculations for a relatively simple RNA model system: theophylline-binding RNA aptamer with theophylline and five analogs. The goal of this investigation was twofold: (1) understanding the performance levels we can expect from absolute free energy calculations for a simple RNA complex and (2) learning about practical modeling considerations that impact the success of RNA binding predictions, which may be different than the best practices established for protein targets. We learned that magnesium ion (Mg ) placement is a critical decision that impacts affinity predictions. When information regarding Mg positions is lacking, implementing RNA backbone restraints is an alternative way of stabilizing RNA structure that recapitulates prediction accuracy. Since mistakes in Mg placement can be detrimental, omitting magnesium ions entirely and using RNA backbone restraints is attractive as a risk-mitigating approach. We found that predictions are sensitive to modeling experimental buffer conditions correctly, including salt type and ionic strength. We explored the effects of sampling in the alchemical protocol, choice of the ligand force field (GAFF2/OpenFF Sage), and water model (TIP3P/OPC) on predictions, which allowed us to give practical advice for the application of free energy methods to RNA targets. By capturing experimental buffer conditions and implementing RNA backbone restraints, we were able to compute binding affinities accurately (MAE = 2.2 kcal/mol, Pearson's correlation coefficient = 0.9, Kendall's tau = 0.7). We believe there is much to learn about how to apply free energy calculations for RNA targets and how to enhance their performance in prospective predictions. This study is an important first step for learning best practices and special considerations for RNA-ligand free energy calculations. Future studies will consider increasingly complicated ligands and diverse RNA systems and help the development of general protocols for therapeutically relevant RNA targets.

Molecular simulations to investigate the impact of N6-methylation in RNA recognition: Improving accuracy and precision of binding free energy prediction

Preparation of Human Nuclear RNA m⁶A Methyltransferases and Demethylases and Biochemical Characterization of Their Catalytic Activity.

Molecular simulations matching denaturation experiments for N6-Methyladenosine

An RNA Methylation-Sensitive AIEgen-Aptamer Reporting System for Quantitatively Evaluating M6a Methylase and Demethylase Activities

Structural and Dynamic Properties of the YTH Domain in Complex with N6-methyladenosine RNA Studied by Accelerated Molecular Dynamics Simulations

Targeted Manipulation of Cellular RNA m6A Methylation at theSingle-Base Level br

Atomistic and Thermodynamic Analysis of N6-Methyladenosine (m6a) Recognition by the Reader Domain of YTHDC1

Mapping Messenger RNA Methylations at Single Base Resolution

Recognition of N6-methyladenosine by the YTHDC1 YTH domain studied by molecular dynamics and NMR spectroscopy: The role of hydration

Structure and thermodynamics of N6-methyladenosine in RNA: a spring-loaded base modification.

Structural and Functional Characterization of the Proteins Responsible for N(6)-Methyladenosine Modification and Recognition.

The catalytic mechanism of the RNA methyltransferase METTL3

Adenine Methylation Enhances the Conformational Flexibility of an RNA Hairpin Tetraloop

Essential Considerations for Free Energy Calculations of RNA-Small Molecule Complexes: Lessons from Theophylline-Binding RNA Aptamer

Theoretical study on the binding mechanism between N6-methyladenine and natural DNA bases

Structural Basis for the Discriminative Recognition of N6-Methyladenosine RNA by the Human YT521-B Homology Domain Family of Proteins.

Identifying RNA N6-Methyladenine Sites in Three Species Based on a Markov Model

Exploring the Binding Process of Cognate Ligand to Add Adenine Riboswitch Aptamer by Using Explicit Solvent Molecular Dynamics (MD) Simulation.

Readers, writers and erasers of N6-methylated adenosine modification.

Empirical corrections to the Amber RNA force field with Target Metadynamics

Identification of a Selective DNA Ligase for Accurate Recognition and Ultrasensitive Quantification of N-6-methyladenosine in RNA at One-Nucleotide Resolution