Most recent works on Brownian dynamics simulation employ a first-order algorithm developed by Ermak and McCammon [J. Chem. Phys. 69, 1352 (1978)]. In this work we propose the use of a second-order algorithm in which the step is a combination of two first-order steps, like in the second-order Runge–Kutta method for differential equations. Although the computer time per step is roughly doubled, the second-order algorithm is more efficient than the previous one because a given accuracy in the results can be achieved with less than half the number of steps. The new algorithm also allows for longer time steps without divergence. The advantage of the new procedure is illustrated in the simulation of four macromolecular systems: A quasirigid dumbbell, a semiflexible trumbbell, a semiflexible hinged rod, and a Gaussian polymer chain.

The translational diffusion coefficients, rotational relaxation times and intrinsic viscosities of rigid bent rods, composed by two rodlike arms joined rigidly at an angle alpha, have been evaluated for varying conformation using the latest advances in hydrodynamic theory. We have considered semiflexible rods in which the joint is an elastic hinge or swivel, with a potential V(alpha) = 1/2Q alpha 2 with constant Q. Accepting the rigid-body treatment, we calculate properties of broken rods by averaging alpha-dependent values for rigid rods. The results are finally used to interpret literature values of the properties of myosin rod. Q is regarded as an adjustable parameter, and the value fitted is such that the average bending angle of myosin rod is approximately 60 degrees.

SynopsisA general treatment for the solution dynamics of segmentally flexible macromolecules having two subunits is presented. Bead modeling allows for a complete inclusion of hydrodynamic interactions in this treatment. "he finite size of the beads is also considered, so that it is therefore possible to account properly for torsional motions of the subunits. Expressions for the components of the resistance matrix are derived. From them, the translational and rotational diffusion coefEcients can be calculated. Distinction is made between hinged macromolecules, whose only internal motion is bending, and swivel-jointed macromolecules, for which torsions of the subunits are also allowed. Numerical results are presented for broken rods with the two types of flexibility.The effects of hydrodynamic interaction between arms of broken rods are about 25% for translation and under 10% for rotation. These findings give support to the treatments of Harvey, Wegener, and co-workers in which interactions were neglected. The rotational dynamics of hinged and swivel-jointed rods are compared. Although there are differences in the short-time behavior, the longest relaxation time is the same for the two cases. Finally, the validity of Wegener's rotational diffusion constants is discussed.

A simple bead model is proposed for the antibody molecule immunoglobulin IgG1. The partial flexibility of the hinge is represented by a quadratic potential associated to the angles between arms. Conformational and hydrodynamic properties are calculated using Monte Carlo (rigid-body) and Brownian dynamics simulations. Comparison of experimental and calculated values for some overall properties allows the assignment of dimensions and other model parameters. The Brownian dynamics technique is used next to simulate a rotational correlation function that is comparable with the decay of fluorescence emission anisotropy. This is done with varying flexibility at the hinge. The longest relaxation time shows a threefold decrease when going from the rigid Y-shaped conformation to the completely flexible case. The calculations are in good agreement with the decay times observed for IgG1. A flexibility analysis of the latter indicates that a variability of +/- 55 degrees (standard deviation) in the angle between the Fab arms.

No abstract