Validating Small-Molecule Force Fields for Macrocyclic Compounds Using NMR Data in Different Solvents

Abstract

Macrocycles are a promising class of compounds as therapeutics for difficult drug targets due to a favorable combination of properties: They often exhibit improved binding affinity compared to their linear counterparts due to their reduced conformational flexibility, while still being able to adapt to environments of different polarity. To assist in the rational design of macrocyclic drugs, there is need for computational methods that can accurately predict conformational ensembles of macrocycles in different environments. Molecular dynamics (MD) simulations remain one of the most accurate methods to predict ensembles quantitatively, although the accuracy is governed by the underlying force field. In this work, we benchmark four different force fields for their application to macrocycles by performing replica exchange with solute tempering (REST2) simulations of 11 macrocyclic compounds and comparing the obtained conformational ensembles to nuclear Overhauser effect (NOE) upper distance bounds from NMR experiments. Especially, the modern force fields OpenFF 2.0 and XFF yield good results, outperforming force fields like GAFF2 and OPLS/AA. We conclude that REST2 in combination with modern force fields can often produce accurate ensembles of macrocyclic compounds. However, we also highlight examples for which all examined force fields fail to produce ensembles that fulfill the experimental constraints.