Foldamer simulations: Novel computational methods and applications to poly-phenylacetylene oligomers

Abstract

We apply several methods to probe the ensemble kinetic and structural properties of a model system of poly-phenylacetylene (pPA) oligomer folding trajectories. The kinetic methods employed included a brute force accounting of conformations, a Markovian state matrix method, and a nonlinear least squares fit to a minimalist kinetic model used to extract the folding time. Each method gave similar measures for the folding time of the 12-mer chain, calculated to be on the order of 7 ns for the complete folding of the chain from an extended conformation. Utilizing both a linear and a nonlinear scaling relationship between the viscosity and the folding time to correct for a low simulation viscosity, we obtain an upper and a lower bound for the approximate folding time within the range 70 ns<tau<350 ns. This is in agreement with the experimentally measured folding time on the order of 160 ns. The kinetic model used to fit the kinetic behavior of the ensemble of trajectories provides a framework to describe the bulk folding mechanism. We were able to identify two unique clusters of conformations that provide a structural basis to account for the appearance of a kinetic intermediate in the mechanism. We discuss the implications of these findings in the context of helix-coil theory.

How Well Can Simulation Predict Protein Folding Kinetics and Thermodynamics?

“…These include, but are not limited to, the characterization of protein folding rates in pores (58), in chaperonins (2), and under force (71), as well as rate predictions for small RNAs (106) and nonbiological polymers (27). Additionally, several coarse-grained calculations have been employed to study folding and unfolding rates (11,50,78).…”

Network models for molecular kinetics and their initial applications to human health

“…Many of the initial applications of MSMs have been to small peptides due to the challenges of simulating larger proteins [12, 13, 17, 20–22, 34, 35]. For example, some of the fundamental methods were laid out by Singhal et al [36] in an early work on a tryptophan zipper β hairpin.…”

Hydrogen Bond Migration between Molecular Sites Observed with Ultrafast 2D IR Chemical Exchange Spectroscopy

“…Phenylacetylene can be polymerized to form poly(phenylacetylene), a highly conjugated polymer which forms folded structures in solution 17. This polymer has been used as model system for investigating polymer folding 17.…”