Kerr constant of multi-subunit particles and semiflexible, wormlike chains.

A formalism is developed for the Kerr constant of particles composed of subunits whose electro-optical properties have axial symmetry. A protocol is devised for the calculation of the Kerr constant from the dipole moments and the electrical and optical polarizabilities of the subunits. The conformational average is required if the particle is flexible. Based on that formalism, the Kerr constant of semiflexible, wormlike chains is calculated with the help of a discrete version to which the previous formalism is applied. The required conformational averages are evaluated by means of Monte Carlo simulation. Thus we obtain expressions for the Kerr constant of wormlike particles over the whole range of conformations.

MULTIHYDRO and MONTEHYDRO: conformational search and Monte Carlo calculation of solution properties of rigid or flexible bead models.

A computer program, MULTIHYDRO, has been constructed for the calculation of hydrodynamic coefficients and other solution properties of multiple possible conformations of a bead model. With minimal additional programming to describe the model under study, this program interfaces efficiently with HYDRO for the calculation of solution properties, including hydrodynamic coefficients, radius of gyration, covolume, etc. A useful application is the conformation search of rigid macromolecules, because many possible conformations can be evaluated in a single run of the program. In this paper we also pay attention to the properties of flexible macromolecules, in the so-called Monte Carlo rigid-body approximation, which is virtually exact for the simpler solution properties. The theoretical aspects of the procedure are described, and we show how MULTIHYDRO can be employed for this calculation. However, for flexible molecules, a more general simulation scheme is importance-sampling Monte Carlo generation. We describe how this procedure is implemented in another computer program, MONTEHYDRO. Examples of the usage of these tools are provided.

Transient electric birefringence of wormlike macromolecules in electric fields of arbitrary strength: a computer simulation study.

We have studied the birefringence decay of linear models of macromolecules for two different types of flexibility, the broken-rod chain and the wormlike chain, using a computer simulation of a transient electric birefringence experiment. We have paid particular attention to the influence of the intensity of the orienting field, including two orienting mechanisms, the induced dipole, and the permanent dipole. We have compared wormlike and broken-rod models of the same radius of gyration, finding that they present a different decay curve under the influence of the same intensity of the field. We have seen that these differences are due to the faster relaxation times (smaller in the wormlike chain model) and amplitudes, because, regardless of the type of flexibility, the overall size of a molecule (measured by the radius of gyration) essentially determines the longest relaxation time. We have also analyzed how the relaxation process is affected by the degree of flexibility, the orientation mechanisms, and the intensity of the field. Studying a different aspect, we have paid attention to the deformation of a molecule in a transient electric birefringence experiment as a source of information. In this work we have developed equations to characterize this deformation in terms of one of the components of the gyration tensor, if a dynamic light scattering experiment under the influence of an electric field could be performed. To develop this work we have simulated the Brownian dynamics of the different models, relaxing after the removal of an orienting external electric field of arbitrary strength. A comparison with other methods such a the rigid body treatment or the correlation analysis of Brownian trajectories has also been included. We have seen that differences between the two Brownian dynamics methods are small and that the rigid-body treatment is only an acceptable approximation to obtain the longest relaxation time.

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.