Two-state protein-like folding of a homopolymer chain

All-or-none folding of a flexible polymer chain in cylindrical nanoconfinement

Geometric confinement of a polymer chain results in a loss of conformational entropy. For a chain that can fold into a compact native state via a first-order-like transition, as is the case for many small proteins, confinement typically provides an entropic stabilization of the folded state, thereby shifting the location of the transition. This allows for the possibility of confinement (entropy) driven folding. Here, we investigate such confinement effects for a flexible square-well-sphere N-mer chain (monomer diameter σ) confined within a long cylindrical pore (diameter D) or a closed cylindrical box (height H = D). We carry out Wang-Landau simulations to construct the density of states, which provides access to the complete thermodynamics of the system. For a wide pore, an entropic stabilization of the folded state is observed. However, as the pore diameter approaches the size of the folded chain (D ∼ N^1/3σ), we find a destabilization effect. For pore diameters smaller than the native ground-state, the chain folds into a different, higher energy, ground state ensemble and the T vs D phase diagram displays non-monotonic behavior as the system is forced into different ground states for different ranges of D. In this regime, isothermal reduction of the confinement dimension can induce folding, unfolding, or crystallite restructuring. For the cylindrical box, we find a monotonic stabilization effect with decreasing D. Scaling laws for the confinement free energy in the athermal limit are also investigated.

Configurational contribution to the Soret effect of a protein ligand system : An investigation with density-of-states simulations

Many of the biological functions of proteins are closely associated with their ability to bind ligands and change conformations in response to changing conditions. Since binding state and conformation of a protein affect its response to a temperature gradient, they may be probed with thermophoresis. In recent years, thermophoretic techniques to investigate biomolecular interactions, quantify ligand binding, and probe conformational changes have become established. To develop a better understanding of the mechanisms underlying the thermophoretic behavior of proteins and ligands, we employ a simple, off-lattice model for a protein and ligand in explicit solvent. To investigate the partitioning of the particles in a temperature gradient, we perform Wang-Landau-type simulations in a divided simulation box and construct the density of states over a two-dimensional state space. This method gives us access to the entropy and energy of the divided system and allows us to estimate the configurational contribution to the Soret coefficient. In this work, we focus on dilute solutions of hydrophobic proteins and investigate the effect of ligand binding on their thermophoretic behavior. We find that our simple model captures important aspects of protein-ligand interactions and allows us to relate the binding energy to the change in Soret coefficient upon ligand binding.

On the polymer physics origins of protein folding thermodynamics

A remarkable feature of the spontaneous folding of many small proteins is the striking similarity in the thermodynamics of the folding process. This process is characterized by simple two-state thermodynamics with large and compensating changes in entropy and enthalpy and a funnel-like free energy landscape with a free-energy barrier that varies linearly with temperature. One might attribute the commonality of this two-state folding behavior to features particular to these proteins (e.g., chain length, hydrophobic/hydrophilic balance, attributes of the native state) or one might suspect that this similarity in behavior has a more general polymer-physics origin. Here we show that this behavior is also typical for flexible homopolymer chains with sufficiently short range interactions. Two-state behavior arises from the presence of a low entropy ground (folded) state separated from a set of high entropy disordered (unfolded) states by a free energy barrier. This homopolymer model exhibits a funneled free energy landscape that reveals a complex underlying dynamics involving competition between folding and non-folding pathways. Despite the presence of multiple pathways, this simple physics model gives the robust result of two-state thermodynamics for both the cases of folding from a basin of expanded coil states and from a basin of compact globule states.

Scaling Laws in Spatial Network Formation

Geometric constraints impact the formation of a broad range of spatial networks, from amino acid chains folding to proteins structures to rearranging particle aggregates. How the network of interactions dynamically self-organizes in such systems is far from fully understood. Here, we analyze a class of spatial network formation processes by introducing a mapping from geometric to graph-theoretic constraints. Combining stochastic and mean field analyses yields an algebraic scaling law for the extent (graph diameter) of the resulting networks with system size, in contrast to logarithmic scaling known for networks without constraints. Intriguingly, the exponent falls between that of self-avoiding random walks and that of space filling arrangements, consistent with experimentally observed scaling of the radius of gyration of protein tertiary structures with their chain length.

Wang-Landau simulation of Gō model molecules

Gō-like models are one of the oldest protein modeling concepts in computational physics and have proven their value over and over for forty years. The essence of a Gō model is to define a native contact matrix for a well-defined low-energy polymer configuration, e.g., the native state in the case of proteins or peptides. Many different potential shapes and many different cut-off distances in the definition of this native contact matrix have been proposed and applied. We investigate here the physical consequences of the choice for this cut-off distance in the Gō models derived for a square-well tangent sphere homopolymer chain. For this purpose we are performing flat-histogram Monte Carlo simulations of Wang-Landau type, obtaining the thermodynamic and structural properties of such models over the complete temperature range. Differences and similarities with Gō models for proteins and peptides are discussed.

Crumpled globule formation during collapse of a long flexible and semiflexible polymer in poor solvent

By introducing explicit solvent particles and hydrodynamic interactions we demonstrate that crumpled globules are formed after the collapse of long polymer chains (N = 10(4)) in a poor solvent. During the collapse crumples of all sizes form sequentially, but small crumples are not stable and convert to blobs with Gaussian statistics. The observed effective mean squared distance R(2)(n) ∼ n(0.38) at n > Ne and contact probability index p(n) ∼ n(-0.5) at n ≫ Ne, which is not following either the model of a fractal globule, or the predictions for an equilibrium globule. Polymer chain stiffness pushes the system to form globular crystallite, and this freezes crumpled structure with R(2)(n) ∼ n(0.33) at n > Ne as a stable state. We note that there is some similarity to crumple globule formation and crystallization of polymer melt.

Theoretical volume profiles as a tool for probing transition states: folding kinetics

The mechanism by which conformational changes, particularly folding and unfolding, occur in proteins and other biopolymers has been widely discussed in the literature. Molecular dynamics (MD) simulations of protein folding present a formidable challenge since these conformational changes occur on a time scale much longer than what can be afforded at the current level of computational technology. Transition state (TS) theory offers a more economic description of kinetic properties of a reaction system by relating them to the properties of the TS, or for flexible systems, the TS ensemble (TSE). The application of TS theory to protein folding is limited by ambiguity in the definition of the TSE for this process. We propose to identify the TSE for conformational changes in flexible systems by comparison of its experimentally determined volumetric property, known as the volume of activation, to the structure-specific volume profile of the process calculated using MD. We illustrate this approach by its successful application to unfolding of a model chain system.

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.

Collapse of a two-state sticky hard sphere chain

The collapse of a homopolymer gaussian chain into a globule is represented as a transition between two states, viz., extended and collapsed. Appropriately, this model has been labeled as the all-or-none view of chain collapse. In the collapsed state, the single polymer partition function is expressed by a single Mayer diagram with the maximum number of f-bonds arising from nonbonded square well interactions. Our target is the dependence of the transition temperature on chain length and the interaction range of the square well, as indicated through the behavior of the radius of gyration and the constant volume heat capacity. Properties of the collapse transition are calculated exactly for chains with three to six backbone atoms and heuristically for long chains using arguments derived from the small chains and from conditions of integrability. Comparison with simulation studies is made.