Basins of attraction for cascading maps

We study a finite uni-directional array of "cascading" or "threshold coupled" chaotic maps. Such systems have been proposed for use in nonlinear computing and have been applied to classification problems in bioinformatics. We describe some of the attractors for such systems and prove general results about their basins of attraction. In particular, we show that the basins of attraction have infinitely many path components. We show that these components always accumulate at the corners of the domain of the system. For all threshold parameters above a certain value, we show that they accumulate at a Cantor set in the interior of the domain. For certain ranges of the threshold, we prove that the system has many attractors.Comment: 15 pages, 9 figures. To appear in International Journal of Bifurcations and Chao

Characteristic Temperatures of Folding of a Small Peptide

We perform a generalized-ensemble simulation of a small peptide taking the interactions among all atoms into account. From this simulation we obtain thermodynamic quantities over a wide range of temperatures. In particular, we show that the folding of a small peptide is a multi-stage process associated with two characteristic temperatures, the collapse temperature T_{\theta} and the folding temperature T_f. Our results give supporting evidence for the energy landscape picture and funnel concept. These ideas were previously developed in the context of studies of simplified protein models, and here for the first time checked in an all-atom Monte Carlo simulation.Comment: Latex, 6 Figure

Three-body Interactions Improve the Prediction of Rate and Mechanism in Protein Folding Models

Here we study the effects of many-body interactions on rate and mechanism in protein folding, using the results of molecular dynamics simulations on numerous coarse-grained C-alpha-model single-domain proteins. After adding three-body interactions explicitly as a perturbation to a Go-like Hamiltonian with native pair-wise interactions only, we have found 1) a significantly increased correlation with experimental phi-values and folding rates, 2) a stronger correlation of folding rate with contact order, matching the experimental range in rates when the fraction of three-body energy in the native state is ~ 20%, and 3) a considerably larger amount of 3-body energy present in Chymotripsin inhibitor than other proteins studied.Comment: 9 pages, 2 tables and 5 figure

Non-Markovian Configurational Diffusion and Reaction Coordinates for Protein Folding

The non-Markovian nature of polymer motions is accounted for in folding kinetics, using frequency-dependent friction. Folding, like many other problems in the physics of disordered systems, involves barrier crossing on a correlated energy landscape. A variational transition state theory (VTST) that reduces to the usual Bryngelson-Wolynes Kramers approach when the non-Markovian aspects are neglected is used to obtain the rate, without making any assumptions regarding the size of the barrier, or the memory time of the friction. The transformation to collective variables dependent on the dynamics of the system allows the theory to address the controversial issue of what are ``good'' reaction coordinates for folding.Comment: 9 pages RevTeX, 3 eps-figures included, submitted to PR

Molecular dynamics simulation of polymer helix formation using rigid-link methods

Molecular dynamics simulations are used to study structure formation in simple model polymer chains that are subject to excluded volume and torsional interactions. The changing conformations exhibited by chains of different lengths under gradual cooling are followed until each reaches a state from which no further change is possible. The interactions are chosen so that the true ground state is a helix, and a high proportion of simulation runs succeed in reaching this state; the fraction that manage to form defect-free helices is a function of both chain length and cooling rate. In order to demonstrate behavior analogous to the formation of protein tertiary structure, additional attractive interactions are introduced into the model, leading to the appearance of aligned, antiparallel helix pairs. The simulations employ a computational approach that deals directly with the internal coordinates in a recursive manner; this representation is able to maintain constant bond lengths and angles without the necessity of treating them as an algebraic constraint problem supplementary to the equations of motion.Comment: 15 pages, 14 figure