Modeling complex biological macromolecules: reduction of multibead models.

Abstract

The shape of simple and complex biological macromolecules can be approximated by bead modeling procedures. Such approaches are required, for example, for the analysis of the scattering and hydrodynamic behavior of the models under analysis and the prediction of their molecular properties. Using the atomic coordinates of proteins for modeling inevitably leads to models composed of a multitude of beads. In particular, for hydrodynamic modeling, a drastic reduction of the bead number may become unavoidable to enable computation. A systematic investigation of different approaches and computation modes shows that the 'running mean', 'cubic grid,' and 'hexagonal grid' approaches are successful, provided that the extent of reduction does not exceed a factor of 100 and the grid approaches use beads of unequal size and the beads are located at the centers of gravity. Further precautions to be taken include usage of appropriate interaction tensors for overlapping beads of unequal size and appropriate volume corrections when calculating intrinsic viscosities. The applied procedures were tested with the small protein lysozyme in a case study and were then applied to the huge capsid of the phage fr and its trimeric building block. The appearance of the models and the agreement of molecular properties and distance distribution functions of unreduced and reduced models can be used as evaluation criteria.