Abstract
The replacement of ionizable functional groups that are predominantly charged at physiological pH with neutral bioisosteres is a common strategy in medicinal chemistry; however, its impact on binding affinity is often context-dependent. Here, we investigated a series of amide derivatives of a glycomimetic E-selectin ligand, in which the carboxylate group of the lead compound is substituted with a range of amide and isosteric analogs. Despite the expected loss of the salt-bridge interaction with Arg97, several amides retained or even improved the binding affinity. Co-crystal structures revealed conserved binding poses across the series, with consistent interactions involving the carbonyl oxygen of the amide and the key residues Tyr48 and Arg97. High-level quantum chemical calculations ruled out a direct correlation between carbonyl partial charges and affinity. Instead, a moderate correlation was observed between ligand binding and the out-of-plane pyramidality of the amide nitrogen, suggesting a favorable steric adaptation within the binding site. Molecular dynamics (MD) simulations revealed that high-affinity ligands exhibit enhanced solution-phase pre-organization toward the bioactive conformation, likely reducing the entropic penalty upon binding. Further analysis of protein-ligand complexes using Molecular mechanics/Generalized born surface area (MM-GB/SA) decomposition suggested minor lipophilic contributions from amide substituents. Taken together, this work underscores the importance of geometric and conformational descriptors, beyond classical electrostatics, in driving affinity in glycomimetic ligand design and provides new insights into the nuanced role of amides as carboxylate isosteres in protein-ligand recognition.