Hydrodynamic properties of macromolecular complexes. V. Improved calculation of rotational diffusion coefficient and intrinsic viscosity

In our previous calculations of rotational diffusion coefficients and intrinsic viscosities of macromolecular complexes modeled by arrays of spherical subunits [J. G. de la Torre & V. A. Bloomfield, Biopolymers 16 , 1765, 1779 (1977); 17 , 1605 (1978)], results were poor when the dominant subunit was located near the center of frictional resistance. A simple means of correcting this flaw, which gives satisfactorily accurate results with little increase in computation time, is to replace the single large subunit with an array of smaller ones. We have examined trigonal bipyramidal, octahedral, and cubic arrays of spheres whose radii were chosen to give the same total volume or the same rotational diffusion coefficient as the parent sphere. These all give similar results, so the details of the modeling are not important. Results obtained using this stratagem are in much better agreement with the theories of Perrin and Simha for short prolate ellipsoids of revolution, and with experimental measurements of rotational diffusion coefficients of T-even bacteriophage without fibers or with fibers retracted. We have also extended our previous calculations to consider phage with various numbers of fibers attached.     

Hydrodynamic properties of macromolecular complexes. IV. Intrinsic viscosity theory,with applications to once-broken rods and multisubunit proteins

We have extended our previous theories of the translational and rotational frictional properties of multisubunit complexes to calculate the intrinsic viscosity of such structures. Our theory is similar to those recently construced by McCammon and Deutch, and by Nakajima and Wada, in that it uses a modified hydrodynamic interaction tensor and solves the system of simultaneous interaction equations by digital computation rather than by successive approximations. However, there are some differences in the formulation and averaging of these equations. Extensive numerical comparison is made between this theory and others that are available—associated with the names of Hearst and Tagami, Abdel-Khalik and Bird, and Tsuda—using as a basis exact results for prolate ellipsoids of revolution. For large axial ratios, only our theory asymptotically approaches the correct limit; but for small axial ratios, only the Tsuda “shell-model” theory is adequate, because the other theories neglect the preponderant influence of the sphere located at the center of rotation. Intrinsic viscosities, translational frictional coefficients, and Scheraga-Mandelkern β parameters, are tabulated for a large number of polygonal and polyhedral subunit structures, with up to eight elements, using both our theory and Tsuda's. Particular application is made to hemerythrin and aspartate transcarbamylase. Finally, the viscosities and friction coefficients o once-broken rods are calculated and compared with an approximate theory by Wilenski.     

The macromolecular properties of peanut agglutinin

The dependence upon solution conditions of the quaternary structure and gross conformation of peanut agglutinin was examined by sedimentation equilibrium, sedimentation velocity, gel chromatography, and circular dichroism. At pH 8, the protein exists as a compactly folded tetramer of molecular weight 98,000. Between pH 4.75 and pH 3.0, the molecular reversibly dissociates to a (still globular) dimer. In the presence of denaturants such as SDS or guanidinium chloride, the protein dissociates to its four equal-sized constituents polypeptide chains. The circular dichroic spectrum of peanut lectin exhibits changes in the near ultraviolet upon binding of lactose, whereas the far ultraviolet spectrum remains unchanged. Dissociation to the dimeric state produces subtle changes in both the near and far ultraviolet circular dichroic spectrum.     

Anomalous properties of water in macromolecular gels.

Low molecular weight solutes often exhibit elution characteristics on gel filtration columns which deviate from ideal behaviour. In many previous studies this anomalous behaviour was attributed to the existence of extremely narrow pores in the gel, inaccessible even to very small solute molecules, to explain Kd values lower than unity. Kd values of small solutes higher than unity were usually ascribed to adsorption of the solute to the gel matrix. In the present paper several observations are presented that contradict these suggestions. Experimental evidence indicates that with small solute molecules Kd values differing from unity can be fully explained by the anomalous properties of vicinal water layers at the gel matrix-water interface.     

Hydrodynamic properties of macromolecular complexes. I. Translation

We have developed an improved theory for calculating the translational frictional coefficients of rigid macromolecular complexes composed of unequal spherical subunits. The Yamakawa hydrodynamic interaction tensor, which improves on the Oseen tensor by taking account of the finite sizes of the frictional subunits, has been generalized to accomodate nonidentical subunits. Iterative numerical methods are described for solving the set of simultaneous hydrodynamic interaction equations, thus avoiding preaveraging. The theory is applied to prolate ellipsoids of revolution, to lollipops, and to dumbbells, and comparison is made with earlier, more approximate theories.     

Microscopic viscosity and rotational diffusion of proteins in a macromolecular environment

The Stokes-Einstein-Debye equation is currently used to obtain information on protein size or on local viscosity from the measurement of the rotational correlation time. However, the implicit assumptions of a continuous and homogeneous solvent do not hold either in vivo, because of the high density of macromolecules, or in vitro, where viscosity is adjusted by adding viscous cosolvents of various size. To quantify the consequence of nonhomogeneity, we have measured the rotational Brownian motion of three globular proteins with molecular mass from 66 to 4000 kD in presence of 1.5 to 2000 kD dextrans as viscous cosolvents. Our results indicate that the linear viscosity dependence of the Stokes-Einstein relation must be replaced by a power law to describe the rotational Brownian motion of proteins in a macromolecular environment. The exponent of the power law expresses the fact that the protein experiences only a fraction of the hydrodynamic interactions of macromolecular cosolvents. An explicit expression of the exponent in terms of protein size and cosolvent's mass is obtained, permitting definition of a microscopic viscosity. Experimental data suggest that a similar effective microviscosity should be introduced in Kramers' equation describing protein reaction rates.     

Statistical interpretation of fluorescence energy transfer measurements in macromolecular systems.

A statistical method is presented for the interpretation of intramolecular distance measurements by the fluorescence energy transfer technique in systems for which the detailed geometries of the donor-acceptor pairs are unknown. This method enables calculation of the probability that a specified distance range corresponds to the actual distance to be measured. It makes use of the numerically calculated probability density function for the distance of interest. The two general systems considered are the single donor-acceptor pair and the multi-donor-single-acceptor transfer. In both systems, the statistical method incorporates the uncertainty in the orientation of the donor and acceptor dipoles. In addition, it can take into account the rotational mobility of the donor dipoles determined by time-dependent emission anisotropy measurements. When more than one donor is involved in the transfer process, the uncertainties associated with the number and location of individual donors and the size and shape of the donor distribution are also incorporated in calculating the distance ranges. Application of the method was demonstrated for a wide range of transfer efficiency and Ro values for the single donor-acceptor system. Specific examples are also presented for interpretation of both single donor-acceptor and multi-donor-single-acceptor energy transfer measurements performed in order to reveal the spatial relationship of the sigma subunit and the rifampicin binding site in the Escherichia coli RNA polymerase (see Wu, C.-W., Yarbrough, L. R., Wu, F. Y.-H., and Hillel, Z. (1976), Biochemistry, preceding paper in this issue). Analysis of these energy transfer data by methods which use average values of the unknown geometrical parameters of the system yielded results similar to those obtained by the statistical method. However, the statistical method represents a more realistic approach to the interpretation of energy transfer measurements since it provides information concerning the entire range of possible distances and their relative likelihood.