Conformational Equilibria of Multimodal Chromatography Ligands in Water and Bound to Protein Surfaces

Multimodal chromatography uses small ligands with multiple modes of interaction, e.g., charged, hydrophobic or hydrogen bonding, to separate proteins from complex mixtures. The mechanism by which multimodal ligands interact with proteins is expected to be affected by ligand conformations, among other factors. Here, we study conformational equilibria of two commercially used multimodal cation exchange ligands, Capto MMC and Nuvia cPrime, in a range of solvents, a Lennard-Jones (LJ) liquid, ethanol, and water, using molecular dynamics (MD) simulations. By mapping ligand conformations onto two key torsion angles, ω and φ, in these solvents and in low and high dielectric media, we quantify the relative importance of intramolecular and solvent-mediated interactions. In a high dielectric medium, Capto MMC preferentially samples three conformations, which are stabilized by a combination of an intramolecular torsion potential (on ω) and LJ interactions. In an LJ liquid, solvent molecules compete with intramolecular interactions while simultaneously providing an osmotic force, stabilizing both closer and farther distances between ligand sites. This has the overall effect of “flattening out” the conformational landscape. Interestingly, in ethanol and water, hydrogen bonding between the amide hydrogen and solvent molecules stabilizes two additional conformations of Capto MMC in which ω takes on less favorable cis-like configurations. MD simulations of ligands in free solution with three therapeutic antibody fragments show that ligand conformational equilibria remain effectively unchanged upon binding to proteins. Although, there is 20–30% dehydration of the overall ligand upon binding, the hydrogen-bonding sites are dehydrated to a much smaller extent, particularly in cis-like configurations. Conformational preferences of Nuvia cPrime are similar to that of Capto MMC, except for the effect of symmetry arising from the absence of an alkyl thiol tail. Characterizing the conformational equilibria of these two ligands in free solution and bound to a protein provides a foundation for developing a mechanistic understanding of protein–multimodal ligand interactions.