Molecular dynamic investigation of the interaction of supported affinity ligands with monoclonal antibodies

Diagnostics and therapeutic treatments based on monoclonal antibodies have been attaining an increasing importance in the past decades, but their large scale employment requires the optimization of purification processes. To obtain this goal, research is focusing on affinity chromatography techniques and the development of new synthetic ligands. In this work we present a computational investigation aimed at obtaining some guidelines for the rational design of affinity ligands, through the study of their interactions with both monoclonal antibodies (modeled as the FC domain of human IgG) and a model support material (agarose). The study was carried out performing molecular dynamics simulations of the support-spacer-ligand-IgG complex in explicit water. Binding energies between IgG and two supported ligands, a disubstituted derivative of trichlorotriazine and a tetrameric peptide, were determined with the linear interaction energy and MM-GBSA approaches. A detailed study of the possible binding sites of the considered ligands was performed exploiting docking protocols and MD simulations. It was found that both ligands bind IgG in the same site as protein A, which is the hinge region between the CH2 and CH3 domains of IgG. However this site is not easily accessible and requires a high mobility of the ligands. The energetic analysis revealed that van der Waals and electrostatic energies of interaction of the triazine ligand with the support are significant and comparable to those with the protein, so that they limit its capability to reach the protein binding site. A similar result was found also for the tetrameric peptide, which is however able to circumvent the problem; for steric reasons only two of its arms can interact at the same time with the agarose support, thus leaving the remaining two available to bind the protein. These results indicate that the interaction between ligand and support material is an important parameter, which should be considered in the computational and experimental design of ligands for affinity chromatography.