Calculation of the solution properties of flexible macromolecules: methods and applications

While the prediction of hydrodynamic properties of rigid particles is nowadays feasible using simple and efficient computer programs, the calculation of such properties and, in general, the dynamic behavior of flexible macromolecules has not reached a similar situation. Although the theories are available, usually the computational work is done using solutions specific for each problem. We intend to develop computer programs that would greatly facilitate the task of predicting solution behavior of flexible macromolecules. In this paper, we first present an overview of the two approaches that are most practical: the Monte Carlo rigid-body treatment, and the Brownian dynamics simulation technique. The Monte Carlo procedure is based on the calculation of properties for instantaneous conformations of the macromolecule that are regarded as if they were instantaneously rigid. We describe how a Monte Carlo program can be interfaced to the programs in the HYDRO suite for rigid particles, and provide an example of such calculation, for a hypothetical particle: a protein with two domains connected by a flexible linker. We also describe briefly the essentials of Brownian dynamics, and propose a general mechanical model that includes several kinds of intramolecular interactions, such as bending, internal rotation, excluded volume effects, etc. We provide an example of the application of this methodology to the dynamics of a semiflexible, wormlike DNA.