Molecular dynamics simulation and experimental validation of the effect of pH on protein desorption in hydrophobic charge induction chromatography

The ligands in hydrophobic charge induction chromatography (HCIC) are hydrophobic and ionisable. Thus, the pH is crucial for the separation performance in HCIC, especially for elution. However, it is difficult to obtain the microscopic information in HCIC through experimental means. In this work, molecular dynamics simulations are performed to examine the effect of pH on elution and protein conformational transition in HCIC, using a 46-bead -barrel coarse-grained model protein and an HCIC adsorbent pore model constructed in an earlier work. Corresponding experiments are carried out for the validation of simulation results, using lysozyme and MEP Hypercel. Both the activities and fluorescence of lysozyme are examined to evaluate the conformational transition. The simulations indicate that the elution efficiency of protein increases with decreasing pH value in a non-linear manner. This is qualitatively consistent with the experimental results. MD simulations indicate that protein unfolding occurs in elution at all pH values. However, the experimental data show that the activity and conformation of lysozyme is independent of pH of the elution buffer. The microscopic information from simulation shows that protein unfolding is mainly observed on the adsorbent surface, but it cannot be detected in the experiments that only probe the proteins in the bulk liquid.