Torsional directed walks, entropic elasticity, and DNA twist stiffness

DNA and other biopolymers differ from classical polymers because of their torsional stiffness. This property changes the statistical character of their conformations under tension from a classical random walk to a problem we call the “torsional directed walk.” Motivated by a recent experiment on single lambda-DNA molecules [Strick, T. R., Allemand, J.-F., Bensimon, D., Bensimon, A. & Croquette, V. (1996) Science 271, 1835–1837], we formulate the torsional directed walk problem and solve it analytically in the appropriate force regime. Our technique affords a direct physical determination of the microscopic twist stiffness C and twist-stretch coupling D relevant for DNA functionality. The theory quantitatively fits existing experimental data for relative extension as a function of overtwist over a wide range of applied force; fitting to the experimental data yields the numerical values C = 120 nm and D = 50 nm. Future experiments will refine these values. We also predict that the phenomenon of reduction of effective twist stiffness by bend fluctuations should be testable in future single-molecule experiments, and we give its analytic form.

Shape anisotropy of a single random-walk polymer

Random walks have been used to describe a wide variety of systems ranging from cell colonies to polymers. Sixty-five years ago, Kuhn [Kuhn, W. (1934) Kolloid-Z. 68, 2–11] made the prediction, backed later by computer simulations, that the overall shape of a random-walk polymer is aspherical, yet no experimental work has directly tested Kuhn's general idea and subsequent computer simulations. By using fluorescence microscopy, we monitored the conformation of individual, long, random-walk polymers (fluorescently labeled DNA molecules) at equilibrium. We found that a polymer most frequently adopts highly extended, nonfractal structures with a strongly anisotropic shape. The ensemble-average ratio of the lengths of the long and short axes of the best-fit ellipse of the polymer was much larger than unity.

A modified gambler's ruin model of polyethylene chains in the amorphous region.

Polyethylene chains in the amorphous region between two crystalline lamellae M unit apart are modeled as random walks with one-step memory on a cubic lattice between two absorbing boundaries. These walks avoid the two preceding steps, though they are not true self-avoiding walks. Systems of difference equations are introduced to calculate the statistics of the restricted random walks. They yield that the fraction of loops is (2M - 2)/(2M + 1), the fraction of ties 3/(2M + 1), the average length of loops 2M - 0.5, the average length of ties 2/3M2 + 2/3M - 4/3, the average length of walks equals 3M - 3, the variance of the loop length 16/15M3 + O(M2), the variance of the tie length 28/45M4 + O(M3), and the variance of the walk length 2M3 + O(M2).

Evidence for nonrandom hydrophobicity structures in protein chains.

The question of whether proteins originate from random sequences of amino acids is addressed. A statistical analysis is performed in terms of blocked and random walk values formed by binary hydrophobic assignments of the amino acids along the protein chains. Theoretical expectations of these variables from random distributions of hydrophobicities are compared with those obtained from functional proteins. The results, which are based upon proteins in the SWISS-PROT data base, convincingly show that the amino acid sequences in proteins differ from what is expected from random sequences in a statistically significant way. By performing Fourier transforms on the random walks, one obtains additional evidence for nonrandomness of the distributions. We have also analyzed results from a synthetic model containing only two amino acid types, hydrophobic and hydrophilic. With reasonable criteria on good folding properties in terms of thermodynamical and kinetic behavior, sequences that fold well are isolated. Performing the same statistical analysis on the sequences that fold well indicates similar deviations from randomness as for the functional proteins. The deviations from randomness can be interpreted as originating from anticorrelations in terms of an Ising spin model for the hydrophobicities. Our results, which differ from some previous investigations using other methods, might have impact on how permissive with respect to sequence specificity protein folding process is-only sequences with nonrandom hydrophobicity distributions fold well. Other distributions give rise to energy landscapes with poor folding properties and hence did not survive the evolution.