ABSTRACT
To study the adsorption of individual peptides in implicit solvent, we propose a version of the Wang-Landau Monte Carlo algorithm that uses a single surface, with no need for a confining wall or grafting. Our "wall-free" method is both more efficient than the traditional ones and free of additional assumptions or approximations. We illustrate it by simulating an HP-model lattice peptide on planar surfaces with a variety of patterns of adsorption sites, discovering a temperature-induced switch of surface selection which is due to a balance of energetic and entropic effects.
Subject(s)
Models, Chemical , Models, Statistical , Peptides/chemistry , Adsorption , Binding Sites , Computer Simulation , Protein Binding , Surface PropertiesABSTRACT
The phase diagram for a lattice ring polymer under applied force, with variable solvent quality, for different topological knot states, is determined for the first time. In addition to eliminating pseudophases where the polymer is flattened into a single layer, it is found that nontrivial knots result in additional pseudophases under tensile force conditions.
Subject(s)
Algorithms , Models, Chemical , Models, Molecular , Phase Transition , Polymers/chemistry , Computer SimulationABSTRACT
The 1/t Wang-Landau algorithm is tested on simple models of polymers and proteins. It is found that this method resolves the problem of the saturation of the error present in the original algorithm for lattice polymers. However, for lattice proteins, which have a rough energy landscape with an unknown energy minimum, it is found that the density of states does not converge in all runs. A new variant of the Wang-Landau algorithm that appears to solve this problem is described and tested. In the new variant, the optimum modification factor is calculated in the same straightforward way throughout the simulation. There is only one free parameter for which a value of unity appears to give near optimal convergence for all run lengths for lattice homopolymers when pull moves are used. For lattice proteins, a much smaller value of the parameter is needed to ensure rapid convergence of the density of states for energies discovered late in the simulation, which unfortunately results in poor convergence early on in the run.
Subject(s)
Algorithms , Polymers/chemistry , Proteins/chemistry , Thermodynamics , Hydrophobic and Hydrophilic Interactions , Monte Carlo MethodABSTRACT
We propose several improvements to the Monte Carlo simulation techniques for lattice peptide adsorption on surfaces. Firstly, we examine the implementation of "pull" moves and discuss the most efficient way of selecting them. Secondly, we explicitly show how Wang-Landau sampling may be used to calculate the appropriate density of states for a peptide chain in contact with a single surface, and how the information from such a simulation may be used to calculate results for slit geometry with a range of wall separations. Lastly, we consider further possible modifications of the simulation method and its application to adsorption on structured and patterned surfaces.