The free energy landscape of a protein-like chain is constructed from exhaustive simulation studies using a combination of discontinuous molecular dynamics and parallel tempering methods. The protein model is a repeating sequence of four kinds of monomers, in which hydrogen bond attraction, electrostatic repulsion, and covalent bond vibrations are modeled by step, shoulder and square-well potentials, respectively. These protein-like chains exhibit a helical structure in their folded states. The model allows a natural definition of a configuration by considering which beads are bonded. In the absence of a solvent, the relative free energy of dominant structures is determined from the relative populations, and the probabilities predicted from the calculated free energies are found to be in excellent agreement with the observed probabilities at different temperatures. The free energy landscape of the protein-like chain is analyzed and confirmed to have funnel-like characteristics, confirmed by the fact that the probability of observing the most common configuration approaches unity at low enough temperatures for chains with fewer than 30 beads. The effect on the free energy landscape of an explicit square-well solvent, where the beads that can form intra-chain bonds can also form (weaker) bonds with solvent molecules while other beads are insoluble, is also examined. Simulations for chains of 15, 20 and 25 beads show that at low temperatures, the most likely structures are collapsed helical structures. The temperature at which collapsed helical structures become dominant is higher than in the absence of a solvent. Finally, the dynamics of the protein-like chain immersed in an implicit hard sphere solvent is studied using a simple model in which the implicit solvent interacts on a fast time scale with the chain beads and provides sufficient friction so that the motion of monomers is governed by the Smoluchowski equation. Using a Markovian model of the kinetics of transitions between conformations, the equilibration process from an ensemble of initially extended configurations to mainly folded configurations is investigated at low effective temperatures for a number of different chain lengths. It was observed that folding profiles appear to be single exponentials and independent of temperature at low temperatures.
Identifer | oai:union.ndltd.org:TORONTO/oai:tspace.library.utoronto.ca:1807/31689 |
Date | 05 January 2012 |
Creators | Bayat Movahed, Hanif |
Contributors | Schofield, Jeremy |
Source Sets | University of Toronto |
Language | en_ca |
Detected Language | English |
Type | Thesis |
Page generated in 0.0013 seconds