Generic folding and transition hierarchies for surface adsorption of hydrophobic-polar lattice model proteins

Thermodynamics of hydrophobic-polar model proteins on the face-centered cubic lattice

The HP model, a coarse-grained protein representation with only hydrophobic (H) and polar (P) amino acids, has already been extensively studied on the simple cubic (SC) lattice. However, this geometry severely restricts possible bond angles, and a simple improvement is to instead use the face-centered cubic (fcc) lattice. In this paper, the density of states and ground state energies are calculated for several benchmark HP sequences on the fcc lattice using the replica-exchange Wang-Landau algorithm and a powerful set of Monte Carlo trial moves. Results from the fcc lattice proteins are directly compared with those obtained from a previous lattice protein folding study with a similar methodology on the SC lattice. A thermodynamic analysis shows comparable folding behavior between the two lattice geometries, but with a greater rate of hydrophobic-core formation persisting into lower temperatures on the fcc lattice.

Simulation study on the critical adsorption and diffusion of polymer chains on heterogeneous surfaces with random adsorption sites

The critical adsorption and diffusion of a linear polymer chain on a heterogeneous surface with randomly distributed adsorption sites are studied using dynamic Monte Carlo simulations. Results show that the critical fraction of the adsorption sites at which critical adsorption takes place decreases exponentially with the increasing polymer-surface attraction strength and, at the same time, decreases with the increasing intra-polymer attraction strength. For adsorbed polymers with large intra-polymer attraction strength, we also find an adsorption-induced structural transition from a three-dimensional compact globule to a two-dimensional compacted pancake with an increasing fraction of adsorption sites. Anomalous sub-diffusion is observed for the adsorbed polymer diffusion on heterogeneous surfaces, in contrast to the normal diffusion on a homogeneous surface. The polymer on heterogeneous surfaces shows larger fluctuation in the total surface attraction energy and a longer waiting time.

Elucidating thermal behavior, native contacts, and folding funnels of simple lattice proteins using replica exchange Wang-Landau sampling

We studied the folding behavior of two coarse-grained, lattice models, the HP (hydrophobic-polar) model and the semi-flexible H0P model, whose 124 monomer long sequences were derived from the protein Ribonuclease A. Taking advantage of advanced parallel computing techniques, we applied replica exchange Wang-Landau sampling and calculated the density of states over the models entire energy ranges to high accuracy. We then determined both energetic and structural quantities in order to elucidate the folding behavior of each model completely. As a result of sufficiently long sequences and model complexity, yet computational accessibility, we were able to depict distinct free energy folding funnels for both models. In particular, we found that the HP model folds in a single-step process with a very highly degenerate native state and relatively flat low temperature folding funnel minimum. By contrast, the semi-flexible H0P model folds via a multi-step process and the native state is almost four orders of magnitude less degenerate than that for the HP model. In addition, for the H0P model, the bottom of the free energy folding funnel remains rough, even at low temperatures.

The role of chain-stiffness in lattice protein models: A replica-exchange Wang-Landau study

Using Monte Carlo simulations, we investigate simple, physically motivated extensions to the hydrophobic-polar lattice protein model for the small (46 amino acid) protein Crambin. We use two-dimensional replica-exchange Wang-Landau sampling to study the effects of a bond angle stiffness parameter on the folding and uncover a new step in the collapse process for particular values of this stiffness parameter. A physical interpretation of the folding is developed by analysis of changes in structural quantities, and the free energy landscape is explored. For these special values of stiffness, we find non-degenerate ground states, a property that is consistent with behavior of real proteins, and we use these unique ground states to elucidate the formation of native contacts during the folding process. Through this analysis, we conclude that chain-stiffness is particularly influential in the low energy, low temperature regime of the folding process once the lattice protein has partially collapsed.

Polymer adsorption on curved surfaces

The conformational behavior of a coarse-grained finite polymer chain near an attractive spherical surface was investigated by means of multicanonical Monte Carlo computer simulations. In a detailed analysis of canonical equilibrium data over a wide range of sphere radius and temperature, we have constructed entire phase diagrams both for nongrafted and end-grafted polymers. For the identification of the conformational phases, we have calculated several energetic and structural observables such as gyration tensor based shape parameters and their fluctuations by canonical statistical analysis. Despite the simplicity of our model, it qualitatively represents in the considered parameter range real systems that are studied in experiments. The work discussed here could have experimental implications from protein-ligand interactions to designing nanosmart materials.

Macroscopically constrained Wang-Landau method for systems with multiple order parameters and its application to drawing complex phase diagrams

A generalized approach to Wang-Landau simulations, macroscopically constrained Wang-Landau, is proposed to simulate the density of states of a system with multiple macroscopic order parameters. The method breaks a multidimensional random-walk process in phase space into many separate, one-dimensional random-walk processes in well-defined subspaces. Each of these random walks is constrained to a different set of values of the macroscopic order parameters. When the multivariable density of states is obtained for one set of values of fieldlike model parameters, the density of states for any other values of these parameters can be obtained by a simple transformation of the total system energy. All thermodynamic quantities of the system can then be rapidly calculated at any point in the phase diagram. We demonstrate how to use the multivariable density of states to draw the phase diagram, as well as order-parameter probability distributions at specific phase points, for a model spin-crossover material: an antiferromagnetic Ising model with ferromagnetic long-range interactions. The fieldlike parameters in this model are an effective magnetic field and the strength of the long-range interaction.

Subphase transitions in first-order aggregation processes

In this paper, we investigate the properties of aggregation transitions in the context of generic coarse-grained homopolymer systems. By means of parallel replica-exchange Monte Carlo methods, we perform extensive simulations of systems consisting of up to 20 individual oligomer chains with five monomers each. Using the tools of the versatile microcanonical inflection-point analysis, we show that the aggregation transition is a first-order process consisting of a sequence of subtransitions between intermediate structural phases. We unravel the properties of these intermediate phases by collecting and analyzing their individual contributions towards the density of states of the system. The central theme of this systematic study revolves around translational entropy and its role in the striking phenomena of missing intermediate phases. We conclude with a brief discussion of the scaling properties of the transition temperature and the latent heat.

Characterizing folding funnels with replica exchange Wang-Landau simulation of lattice proteins

We have studied the folding of ribonuclease A by mapping it onto coarse-grained lattice protein models. With replica exchange Wang-Landau sampling, we calculated the free energy vs end-to-end distance as a function of temperature. A mapping to the famous hydrophobic-polar (HP) model shows a relatively shallow folding funnel and flat free energy minimum, reflecting the high degeneracy of the ground state. In contrast, extending the HP model with an additional "neutral" monomer type (i.e., a mapping to the three-letter H0P model) has a well developed, rough free energy funnel with a low degeneracy ground state. In both cases, folding funnels are asymmetric with temperature dependent shape.

Dynamics of adsorbed polymers on attractive homogeneous surfaces

Dynamic behaviors of polymer chains adsorbed on an attractive, homogeneous surface are studied by using dynamic Monte Carlo simulations. The translational diffusion coefficient D_xy parallel to the surface decreases as the intra-polymer attraction strength E_PP or the polymer-surface attraction strength E_PS increases. The rotational relaxation time τ_R increases with E_PS, but the dependence of τ_R on E_PP is dependent on the adsorption state of the polymer. We find that τ_R decreases with increasing E_PP for a partially adsorbed polymer but it increases with E_PP for a fully adsorbed polymer. Scaling relations D_xy ~ N^−α and τ_R ~ N^β are found for long polymers. The scaling exponent α is independent of E_PS for long polymers but increases with E_PP from α = 1.06 at E_PP = 0. While β ≈ 2.7 is also roughly independent of E_PS for the adsorbed polymer at E_PP = 0, but β increases with E_PS at E_PP > 0. Moreover, we find that β always decreases with increasing E_PP. Our results reveal different effects of the attractive surface on the diffusion and rotation of adsorbed polymers.

Interlocking order parameter fluctuations in structural transitions between adsorbed polymer phases

By means of contact-density chain-growth simulations of a simple coarse-grained lattice model for a polymer grafted at a solid homogeneous substrate, we investigate the complementary behavior of the numbers of surface-monomer and monomer-monomer contacts under various solvent and thermal conditions. This pair of contact numbers represents an appropriate set of order parameters that enables the distinct discrimination of significantly different compact phases of polymer adsorption. Depending on the transition scenario, these order parameters can interlock in perfect cooperation. The analysis helps understand the transitions from compact filmlike adsorbed polymer conformations into layered morphologies and dissolved adsorbed structures, respectively, in more detail.

Confinement effects upon the separation of structural transitions in linear systems with restricted bond fluctuation ranges

By means of advanced parallel replica-exchange Monte Carlo methods we examine the influence of elasticity and confinement on the structural transitions of linear systems with restricted bonded interaction. For this purpose, we adopt a model for coarse-grained flexible polymers of finite length in the dilute regime. Hyperphase diagrams are constructed using energy-dependent canonical quantities to demonstrate the effects of the changes in the range of the confined interaction on the liquid and solid structural phases. With increasing bonded interaction range we observe the disappearance of the liquid phase and the fusion of the gas-liquid (or Θ) and the liquid-solid transitions. One of the most remarkable features, the liquid-gas transition, changes from second to first order if the confined interaction range exceeds a threshold that separates polymeric from nonpolymeric systems. The notoriously difficult sampling of the entropically suppressed conformations in the region of very strong first-order transitions is improved by using multiple Gaussian modified ensembles.

Simple model of pH-induced protein denaturation

The pH-induced conformational changes of proteins are systematically studied in the framework of a hydrophobic-polar (HP) model, in which proteins are dramatically simplified as chains of hydrophobic (H) and polar (P) beads on a lattice. We express the electrostatic interaction, the principal driving force of pH-induced unfolding that is not included in the conventional HP model, as the repulsive energy term between P monomers. As a result of the exact enumeration of all of the 14- to 18-mers, it is found that lowest-energy states in many sequences change from single "native" conformations to multiple sets of "denatured" conformations with an increase in the electrostatic repulsion. The switching of the lowest-energy states occurs in quite a similar way to real proteins: it is almost always between two states, while in a small fraction of ≥16-mers it is between three states. We also calculate the structural fluctuations for all of the denatured states and find that the denatured states contain a broad range of incompletely unfolded conformations, similar to "molten globule" states referred to in acid or alkaline denatured real proteins. These results show that the proposed model provides a simple physical picture of pH-induced protein denaturation.

Interplay of Coil-Globule Transition and Surface Adsorption of a Lattice HP Protein Model

An end-grafted hydrophobic-polar (HP) model protein chain with alternating H and P monomers is studied to examine interactions between the critical adsorption transition due to surface attraction and the collapse transition due to pairwise attractive H–H interactions. We find that the critical adsorption phenomenon can always be observed; however, the critical adsorption temperature T_CAP is influenced by the attractive H–H interactions in some cases. When the collapse temperature T_c is lower than T_CAP, the critical adsorption of the HP chain is similar to that of a homopolymer without intrachain attractions and T_CAP remains unchanged, whereas the collapse transition is suppressed by the adsorption. In contrast, for cases where T_c is close to or higher than T_CAP, T_CAP of the HP chain is increased, indicating that a collapsed chain is more easily adsorbed on the surface. The strength of the H–H attraction also influences the statistical size and shape of the polymer, with strong H–H attractions resulting in adsorbed and collapsed chains adopting two-dimensional, circular conformations.

Effect of single-site mutations on hydrophobic-polar lattice proteins

We developed a heuristic method for determining the ground-state degeneracy of hydrophobic-polar (HP) lattice proteins, based on Wang-Landau and multicanonical sampling. It is applied during comprehensive studies of single-site mutations in specific HP proteins with different sequences. The effects in which we are interested include structural changes in ground states, changes of ground-state energy, degeneracy, and thermodynamic properties of the system. With respect to mutations, both extremely sensitive and insensitive positions in the HP sequence have been found. That is, ground-state energies and degeneracies, as well as other thermodynamic and structural quantities, may be either largely unaffected or may change significantly due to mutation.

Scalable replica-exchange framework for Wang-Landau sampling

We investigate a generic, parallel replica-exchange framework for Monte Carlo simulations based on the Wang-Landau method. To demonstrate its advantages and general applicability for massively parallel simulations of complex systems, we apply it to lattice spin models, the self-assembly process in amphiphilic solutions, and the adsorption of molecules on surfaces. While of general current interest, the latter phenomena are challenging to study computationally because of multiple structural transitions occurring over a broad temperature range. We show how the parallel framework facilitates simulations of such processes and, without any loss of accuracy or precision, gives a significant speedup and allows for the study of much larger systems and much wider temperature ranges than possible with single-walker methods.

Adsorption and pattern recognition of polymers at complex surfaces with attractive stripelike motifs

We construct the complete structural phase diagram of polymer adsorption at substrates with attractive stripelike patterns in the parameter space spanned by the adsorption affinity of the stripes and temperature. Results were obtained by extensive generalized-ensemble Monte Carlo simulations of a generic model for the hybrid organic-inorganic system. By comparing with adhesion properties at homogeneous substrates, we find substantial differences in the formation of adsorbed polymer structures if translational invariance at the surface is broken by a regular pattern. Beside a more specific understanding of polymer adsorption processes, our results are potentially relevant for the design of macromolecular pattern recognition devices such as sensors.

Dynamics of a polymer adsorbed to an attractive homogeneous flat surface

Polymer conformation and statistical sizes change with the surface contact number during the adsorption.

Effects of knots on protein folding properties

This work explores the impact of knots, knot depth and motif of the threading terminus in protein folding properties (kinetics, thermodynamics and mechanism) via extensive Monte Carlo simulations of lattice models. A knotted backbone has no effect on protein thermodynamic stability but it may affect key aspects of folding kinetics. In this regard, we found clear evidence for a functional advantage of knots: knots enhance kinetic stability because a knotted protein unfolds at a distinctively slower rate than its unknotted counterpart. However, an increase in knot deepness does not necessarily lead to more effective changes in folding properties. In this regard, a terminus with a non-trivial conformation (e.g. hairpin) can have a more dramatic effect in enhancing kinetic stability than knot depth. Nevertheless, our results suggest that the probability of the denatured ensemble to keep knotted is higher for proteins with deeper knots, indicating that knot depth plays a role in determining the topology of the denatured state. Refolding simulations starting from denatured knotted conformations show that not every knot is able to nucleate folding and further indicate that the formation of the knotting loop is a key event in the folding of knotted trefoils. They also show that there are specific native contacts within the knotted core that are crucial to keep a native knotting loop in denatured conformations which otherwise have no detectable structure. The study of the knotting mechanism reveals that the threading of the knotting loop generally occurs towards late folding in conformations that exhibit a significant degree of structural consolidation.

Generic, hierarchical framework for massively parallel Wang-Landau sampling

We introduce a parallel Wang-Landau method based on the replica-exchange framework for Monte Carlo simulations. To demonstrate its advantages and general applicability for simulations of complex systems, we apply it to different spin models including spin glasses, the Ising model, and the Potts model, lattice protein adsorption, and the self-assembly process in amphiphilic solutions. Without loss of accuracy, the method gives significant speed-up and potentially scales up to petaflop machines.