Emergence of stable and fast folding protein structures

The number of protein structures is far less than the number of sequences. By imposing simple generic features of proteins (low energy and compaction) on all possible sequences we show that the structure space is sparse compared to the sequence space. Even though the sequence space grows exponentially with N (the number of amino acids) we conjecture that the number of low energy compact structures only scales as ln N. This implies that many sequences must map onto countable number of basins in the structure space. The number of sequences for which a given fold emerges as a native structure is further reduced by the dual requirements of stability and kinetic accessibility. The factor that determines the dual requirement is related to the sequence dependent temperatures, T_\theta (collapse transition temperature) and T_F (folding transition temperature). Sequences, for which \sigma =(T_\theta-T_F)/T_\theta is small, typically fold fast by generically collapsing to the native-like structures and then rapidly assembling to the native state. Such sequences satisfy the dual requirements over a wide temperature range. We also suggest that the functional requirement may further reduce the number of sequences that are biologically competent. The scheme developed here for thinning of the sequence space that leads to foldable structures arises naturally using simple physical characteristics of proteins. The reduction in sequence space leading to the emergence of foldable structures is demonstrated using lattice models of proteins.Comment: latex, 18 pages, 8 figures, to be published in the conference proceedings "Stochastic Dynamics and Pattern Formation in Biological Systems

Dependence of folding rates on protein length

Using three-dimensional Go lattice models with side chains for proteins, we investigate the dependence of folding times on protein length. In agreement with previous theoretical predictions, we find that the folding time grows as a power law with the chain length N with exponent $\lambda \approx 3.6$ for the Go model, in which all native interactions (i.e., between all side chains and backbone atoms) are uniform. If the interactions between side chains are given by pairwise statistical potentials, which introduce heterogeneity in the contact energies, then the power law fits yield large $\lambda$ values that typically signifies a crossover to an underlying activated process. Accordingly, the dependence of folding time is best described by the stretched exponential \exp(\sqrt{N}). The study also shows that the incorporation of side chains considerably slows down folding by introducing energetic and topological frustration.Comment: 6 pages, 5 eps figure

Probing the Mechanisms of Fibril Formation Using Lattice Models

Using exhaustive Monte Carlo simulations we study the kinetics and mechanism of fibril formation using lattice models as a function of temperature and the number of chains. While these models are, at best, caricatures of peptides, we show that a number of generic features thought to govern fibril assembly are present in the toy model. The monomer, which contains eight beads made from three letters (hydrophobic, polar, and charged), adopts a compact conformation in the native state. The kinetics of fibril assembly occurs in three distinct stages. In each stage there is a cascade of events that transforms the monomers and oligomers to ordered structures. In the first "burst" stage highly mobile oligomers of varying sizes form. The conversion to the aggregation-prone conformation occurs within the oligomers during the second stage. As time progresses, a dominant cluster emerges that contains a majority of the chains. In the final stage, the aggregation-prone conformation particles serve as a template onto which smaller oligomers or monomers can dock and undergo conversion to fibril structures. The overall time for growth in the latter stages is well described by the Lifshitz-Slyazov growth kinetics for crystallization from super-saturated solutions.Comment: 27 pages, 6 figure

Multiexcitons confined within a sub-excitonic volume: Spectroscopic and dynamical signatures of neutral and charged biexcitons in ultrasmall semiconductor nanocrystals

The use of ultrafast gating techniques allows us to resolve both spectrally and temporally the emission from short-lived neutral and negatively charged biexcitons in ultrasmall (sub-10 nm) CdSe nanocrystals (nanocrystal quantum dots). Because of forced overlap of electronic wave functions and reduced dielectric screening, these states are characterized by giant interaction energies of tens (neutral biexcitons) to hundreds (charged biexcitons) of meV. Both types of biexcitons show extremely short lifetimes (from sub-100 picoseconds to sub-picosecond time scales) that rapidly shorten with decreasing nanocrystal size. These ultrafast relaxation dynamics are explained in terms of highly efficient nonradiative Auger recombination.Comment: 5 pages, 4 figures, to be published in Phys. Rev.

Finite size effects on thermal denaturation of globular proteins

Finite size effects on the cooperative thermal denaturation of proteins are considered. A dimensionless measure of cooperativity, Omega, scales as N^zeta, where N is the number of amino acids. Surprisingly, we find that zeta is universal with zeta = 1 + gamma, where the exponent gamma characterizes the divergence of the susceptibility for a self-avoiding walk. Our lattice model simulations and experimental data are consistent with the theory. Our finding rationalizes the marginal stability of proteins and substantiates the earlier predictions that the efficient folding of two-state proteins requires the folding transition temperature to be close to the collapse temperature.Comment: 3 figures. Physical Review Letters (in press