A new simulation of branched polymers

Flory theory of randomly branched polymers

Randomly branched polymer chains (or trees) are a classical subject of polymer physics with connections to the theory of magnetic systems, percolation and critical phenomena. More recently, the model has been reconsidered for RNA, supercoiled DNA and the crumpling of topologically-constrained polymers. While solvable in the ideal case, little is known exactly about randomly branched polymers with volume interactions. Flory theory provides a simple, unifying description for a wide range of branched systems, including isolated trees in good and θ-solvent, and tree melts. In particular, the approach provides a common framework for the description of randomly branched polymers with quenched connectivity and for randomly branching polymers with annealed connectivity. Here we review the Flory theory for interacting trees in the asymptotic limit of very large polymerization degree for good solvent, θ-solutions and melts, and report its predictions for annealed connectivity in θ-solvents. We compare the predictions of Flory theory for randomly branched polymers to a wide range of available analytical and numerical results and conclude that they are qualitatively excellent and quantitatively good in most cases.

Mean Wiener numbers and other mean extensions for alkane trees

The "Wiener number" (the sum over intersite graph distances of a structure) as averaged over all alkane structural isomers of a fixed number N of carbon atoms is considered. This and several other measures of average graphical "extension" of N-site alkanes are computed for N up to 90 (where there are over 10(35) such isomers). Fits are then made for several surmised or derived asymptotic forms, and a heuristic argument is made relating these results to geometric extensions of a random mix of (N-site) alkanes.

Computer simulation of the free energy of peptides with the local states method: analogues of gonadotropin releasing hormone in the random coil and stable states

The Helmholtz free energy F (rather than the energy) is the correct criterion for stability; therefore, calculation of F is important for peptides and proteins that can populate a large number of metastable states. The local states (LS) method proposed by H. Meirovitch [(1977) Chemical Physics Letters, Vol. 45, p. 389] enables one to obtain upper and lower bounds of the conformational free energy, FB (b,l) and FA (b,l), respectively, from molecular dynamics (MD) or Monte Carlo samples. The correlation parameter b is the number of consecutive dihedral or valence angles along the chain that are taken into account explicitly. The continuum angles are approximated by a discretization parameter l; the larger are b and l, the better the approximations; while FA can be estimated efficiently, it is more difficult to estimate FB. The method is further developed here by applying it to MD trajectories of a relatively large molecule (188 atoms), the potent "Asp4-Dpr10" antagonist [cyclo(4/10)-(Ac-delta 3Pro1-D-pFPhe2-D-Trp3-Asp4-Tyr5-D-Nal6-Leu7-Arg8 -Pro9- Dpr10-NH2)] of gonadotropin releasing hormone (GnRH). The molecule was simulated in vacuo at T = 300 K in two conformational states, previously investigated [J. Rizo et al. Journal of the American Chemical Society, (1992) Vol. 114, p. 2860], which differ by the orientation of the N-terminal tail, above (tail up, TU) and below (tail down, TD) the cyclic heptapeptide ring. As in previous applications of the LS method, we have found the following: (1) While FA is a crude approximation for the correct F, results for the difference, delta FA = FA (TD)-FA (TU) converge rapidly to 5.6 (1) kcal/mole as the approximation is improved (i.e., as b and l are increased), which suggests that this is the correct value for delta F; therefore TD is more stable than TU. (The corresponding difference in entropy, T delta SA = 1.3(2) kcal/mole, is equal to the value obtained by the harmonic approximation.) (2) The lowest approximation, which has the minimal number of local states, i.e., based on b = 0 (no correlations) and l = 1 (the angle values are distributed homogeneously), also leads to the correct value of delta F, within the error bars. This is important since the lowest approximation can be applied even to large proteins. (3) The method enables one to define the entropy of a part of the molecule and thus to measure the flexibility of this part.(ABSTRACT TRUNCATED AT 400 WORDS)