Modeling of protein-anion exchange resin interaction for the human growth hormone charge variants

Modeling ion exchange chromatography (IEC) behavior has generated significant interest because of the wide use of IEC as an analytical technique as well as a preparative protein purification process; indeed there is a need for better understanding of what drives the unique behavior of protein charge variants. We hypothesize that a complex protein molecule, which contains both hydrophobic and charged moieties, would interact strongly with an in silico designed resin through charged electrostatic patches on the surface of the protein. In the present work, variants of recombinant human growth hormone that mimic naturally-occurring deamidation products were produced and characterized in silico. The study included these four variants: rhGH, N149D, N152D, and N149D/N152D. Poisson-Boltzmann calculations were used to determine surface electrostatic potential. Metropolis Monte Carlo simulations were carried out with the resulting variants to simulate IEC systems, examining the free energy of the interaction of the protein with an in silico anion exchange column represented by polylysine polypeptide. The results show that the charge variants have different average binding energies and the free energy of interaction can be used to predict the retention time for the different variants.