Proteins with selected sequences: a heteropolymeric study

Equilibrium properties of realistic random heteropolymers and their relevance for globular and naturally unfolded proteins

Random heteropolymers do not display the typical equilibrium properties of globular proteins, but are the starting point to understand the physics of proteins and, in particular, to describe their non-native states. So far, they have been studied with mean-field models in the thermodynamic limit, or with computer simulations of very small chains on lattice. After describing a self-adjusting parallel-tempering technique to sample efficiently the low-energy states of frustrated systems without the need of tuning the system-dependent parameters of the algorithm, we apply it to random heteropolymers moving in continuous space. We show that if the mean interaction between monomers is negative, the usual description through the random-energy model is nearly correct, provided that it is extended to account for noncompact conformations. If the mean interaction is positive, such a simple description breaks out and the system behaves in a way more similar to Ising spin glasses. The former case is a model for the denatured state of globular proteins, the latter of naturally unfolded proteins, whose equilibrium properties thus result as qualitatively different.

Untangling influences of hydrophobicity on protein sequences and structures

We perform a statistical analysis of solvent accessibility and hydrophobicity profiles of a representative set of proteins. The joint probability distribution is well fitted to a multivariable Gaussian, which takes a relatively simple form when expressed in terms of the Fourier transforms of the profiles. This allows us to quantify the asymmetric manner by which these profiles influence each other. For example, the alpha-helix periodicity in sequence hydrophobicity is dictated by the solvent accessibility of structures, and not vice versa, possibly indicating the faster evolution of sequences compared to structures. The decorrelated hydrophobicity and solvent accessibility profiles show distinct behaviors at long periods, where sequence hydrophobicity fluctuates less, while solvent accessibility fluctuates more than average. The correlations between the two profiles can be interpreted as the Boltzmann weight of the solvation energy at room temperature, consistent with earlier observations.

Embedding a native state into a random heteropolymer model: the dynamic approach

We study a random heteropolymer model with Langevin dynamics, in the supersymmetric formulation. Employing a procedure similar to one that has been used in static calculations, we construct an ensemble in which the affinity of the system for a native state is controlled by a "selection temperature" T0. In the limit of high T0, the model reduces to a random heteropolymer, while for T0-->0 the system is forced into the native state. Within the Gaussian variational approach that we employed previously for the random heteropolymer, we explore the phases of the system for high and low T0. For high T0, the system exhibits a (dynamical) spin-glass phase, like that found for the random heteropolymer, below a temperature T(g). For low T0, we find an ordered phase, characterized by a nonzero overlap with the native state, below a temperature T(n) proportional to 1/T(0)>T(g). However, the random-globule phase remains locally stable below T(n), down to the dynamical glass transition at T(g). Thus, in this model, folding is rapid for temperatures between T(g) and T(n), but below T(g) the system can get trapped in conformations uncorrelated with the native state. At a lower temperature, the ordered phase can also undergo a dynamical glass transition, splitting into substates separated by large barriers.

Random heteropolymer dynamics

We study the Langevin dynamics of the standard random heteropolymer model by mapping the problem to a supersymmetric field theory using the Martin-Siggia-Rose formalism. The resulting model is solved nonperturbatively employing a Gaussian variational approach. In constructing the solution, we assume that the chain is very long and impose the translational invariance which is expected to be present in the bulk of the globule by averaging over the center of mass coordinate. In this way we derive equations of motion for the correlation and response functions C(t,t') and R(t,t'). The order parameters are extracted from the asymptotic behavior of these functions. We find a dynamical phase diagram with frozen (glassy) and melted (ergodic) phases. In the glassy phase the system fails to reach equilibrium and exhibits aging of the type found in p-spin glasses. Within the approximations used in this study, the random heteropolymer model can be mapped to the problem of a manifold in a random potential with power law correlations.