Rotational Brownian dynamics of semiflexible broken rods

Dynamic electro-optic properties of macromolecules and nanoparticles in solution: A review of computational and simulation methodologies

This paper reviews some theories, and computational and simulation procedures available for the calculation of the time-course of electro-optic properties of particles in solution. For rigid particles, the time evolution of the properties is directly related to their rotational diffusion; therefore, the computational procedures for the calculation of hydrodynamic properties find a direct application in electro-optics. Several of such computational procedures, based on bead models, are reviewed. For flexible particles, the simultaneous effects the external field and the flexibility can be treated with Brownian dynamics simulation. We illustrate the various procedures, with applications to rigid bent rods and flexible, wormlike or hinged rods, trying to show how the absence or presence of flexibility, and its kind, influences the dynamic electro-optic properties, which are therefore valuable sources of information about the conformation of macromolecules and nanoparticles.

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.

Hydrodynamics of segmentally flexible macromolecules

Segmentally flexible macromolecules are composed of a few rigid subunits linked by joints which are more or less flexible. The dynamics in solution of this type of macromolecule present special aspects that are reviewed here. Three alternative approaches are described. One is the rigid-body treatment, which is shown to be valid for overall dynamic properties such as translational diffusion and intrinsic viscosity. Another approach is the Harvey-Wegener treatment, which is particularly suited for rotational diffusion. The simplest version of this treatment, which ignores hydrodynamic interaction (HI) effects, is found to be quite accurate when compared to a more rigorous version including HI. A third approach is the Brownian dynamics simulation that, albeit at some computational cost, might describe rigorously cases of arbitrary complexity. This technique has been used to test the approximations in the rigid-body and Harvey-Wegener treatments, thus allowing a better understanding of their validity. Brownian trajectories of simplified models such as the trumbbell and the broken rod have been simulated. The comparison of the decay rates of some correlation functions with the predictions of the two treatments leads to a general conclusion: the Harvey-Wegener treatment determines the initial rate, while the long-time behavior is dominated by the rigid-body relaxation time. As an example of application to a specific biological macromolecule, we present a simulation of an immunoglobulin molecule, showing how Brownian Dynamics can be used to predict rotational and internal dynamics. Another typical example is myosin. Literature data of hydrodynamic properties of whole myosin and the myosin rod are compared with predictions from the Harvey-Wegener and rigid-body treatments. The present situation of the problem on myosin flexibility is analyzed, and some indications are given for future experimental and simulation work.