Computer simulation of antibody binding specificity

A Monte Carlo algorithm that searches for the optimal docking configuration of hen egg white lysozyme to an antibody is developed. Both the lysozyme and the antibody are kept rigid. Unlike the work of other authors, our algorithm does not attempt to explicitly maximize surface contact, but minimizes the energy computed using coarse-grained pair potentials. The final refinement of our best solutions using all-atom OPLS potentials (Jorgensen and Tirado-Rives8) consistently yields the native conformation as the preferred solution for three different antibodies. We find that the use of an exponential distance-dependent dielectric function is an improvement over the more commonly used linear form.