Foldamer dynamics expressed via Markov state models. I. Explicit solvent molecular-dynamics simulations in acetonitrile, chloroform, methanol, and water.

In this article, we analyze the folding dynamics of an all-atom model of a polyphenylacetylene (pPA) 12-mer in explicit solvent for four common organic and aqueous solvents: acetonitrile, chloroform, methanol, and water. The solvent quality has a dramatic effect on the time scales in which pPA 12-mers fold. Acetonitrile was found to manifest ideal folding conditions as suggested by optimal folding times on the order of approximately 100-200 ns, depending on temperature. In contrast, chloroform and water were observed to hinder the folding of the pPA 12-mer due to extreme solvation conditions relative to acetonitrile; chloroform denatures the oligomer, whereas water promotes aggregation and traps. The pPA 12-mer in a pure methanol solution folded in approximately 400 ns at 300 K, compared relative to the experimental 12-mer folding time of approximately 160 ns measured in a 1:1 v/v THF/methanol solution. Requisite in drawing the aforementioned conclusions, analysis techniques based on Markov state models are applied to multiple short independent trajectories to extrapolate the long-time scale dynamics of the 12-mer in each respective solvent. We review the theory of Markov chains and derive a method to impose detailed balance on a transition-probability matrix computed from simulation data.

[1]  Michael R. Shirts,et al.  Extremely precise free energy calculations of amino acid side chain analogs: Comparison of common molecular mechanics force fields for proteins , 2003 .

[2]  E. W. Meijer,et al.  Self-assembly of folded m-phenylene ethynylene oligomers into helical columns. , 2001, Journal of the American Chemical Society.

[3]  A. Fersht On the simulation of protein folding by short time scale molecular dynamics and distributed computing , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[4]  V. Pande,et al.  Foldamer dynamics expressed via Markov state models. II. State space decomposition. , 2005, The Journal of chemical physics.

[5]  N. Kampen,et al.  Stochastic processes in physics and chemistry , 1981 .

[6]  Jeffrey S. Moore,et al.  Foldamer-Based Molecular Recognition , 2000 .

[7]  Vijay S Pande,et al.  Length dependent folding kinetics of phenylacetylene oligomers: structural characterization of a kinetic trap. , 2005, The Journal of chemical physics.

[8]  Michael R. Shirts,et al.  Solvation free energies of amino acid side chain analogs for common molecular mechanics water models. , 2005, The Journal of chemical physics.

[9]  J S Moore,et al.  Solvophobically driven folding of nonbiological oligomers. , 1997, Science.

[10]  Michael Heylin SALARY SURVEY: Chemists are swept along by the longest decline in jobs since World War II, bringing a record high level of unemployment for ACS members , 2003 .

[11]  H. Stanley,et al.  Statistical physics of macromolecules , 1995 .

[12]  William Swope,et al.  Describing Protein Folding Kinetics by Molecular Dynamics Simulations. 2. Example Applications to Alanine Dipeptide and a β-Hairpin Peptide† , 2004 .

[13]  Jeffery G. Saven,et al.  Cooperative Conformational Transitions in Phenylene Ethynylene Oligomers: Chain-Length Dependence , 1999 .

[14]  Gin,et al.  Helical twist sense bias in oligo(phenylene ethynylene)s induced by an optically active flexible tether , 2000, Organic letters.

[15]  David J. Hill,et al.  Helicogenicity of solvents in the conformational equilibrium of oligo(m-phenylene ethynylene)s: Implications for foldamer research , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[16]  A New Twist on the Helix-Coil Transition: A Non-Biological Helix with Protein-Like Intermediates and Traps , 2001 .

[17]  George H. Weiss,et al.  Stochastic Processes in Chemical Physics: The Master Equation , 1977 .

[18]  J. Saven,et al.  Simulation Studies of a Helical m-Phenylene Ethynylene Foldamer , 2004 .

[19]  Moore,et al.  Twist Sense Bias Induced by Chiral Side Chains in Helically Folded Oligomers. , 2000, Angewandte Chemie.

[20]  Eric J. Sorin,et al.  Does native state topology determine the RNA folding mechanism? , 2004, Journal of molecular biology.

[21]  Martin Gruebele,et al.  Transition from Exponential to Nonexponential Kinetics during Formation of a Nonbiological Helix , 2000 .

[22]  Jeffrey S. Moore,et al.  Folding-driven synthesis of oligomers , 2001, Nature.

[23]  P. Deuflhard,et al.  Identification of almost invariant aggregates in reversible nearly uncoupled Markov chains , 2000 .

[24]  Jeffrey S. Moore,et al.  Reversible polymerization driven by folding. , 2002, Journal of the American Chemical Society.

[25]  Matthew J. Mio,et al.  A field guide to foldamers. , 2001, Chemical reviews.

[26]  Jeffrey S. Moore Shape-Persistent Molecular Architectures of Nanoscale Dimension , 1997 .

[27]  N. Mikami,et al.  Electronic spectra of tolane in a supersonic free jet: Large-amplitude torsional motion , 1984 .

[28]  T. S. Hughes,et al.  Foldamers as dynamic receptors: probing the mechanism of molecular association between helical oligomers and rodlike ligands. , 2002, Angewandte Chemie.

[29]  Jeffrey S. Moore,et al.  A water-soluble m-phenylene ethynylene foldamer. , 2004, Organic letters.

[30]  Jeffrey S. Moore,et al.  Hexagonal Packing of Oligo(m-phenylene ethynylene)s in the Solid State: Helical Nanotubules , 2000 .

[31]  Vijay S Pande,et al.  Foldamer simulations: novel computational methods and applications to poly-phenylacetylene oligomers. , 2004, The Journal of chemical physics.

[32]  K. Shuler,et al.  Stochastic processes in chemical physics , 1969 .

[33]  Vijay S Pande,et al.  Using path sampling to build better Markovian state models: predicting the folding rate and mechanism of a tryptophan zipper beta hairpin. , 2004, The Journal of chemical physics.

[34]  William Swope,et al.  Describing Protein Folding Kinetics by Molecular Dynamics Simulations. 1. Theory , 2004 .

[35]  Matthew J. Mio,et al.  Chain length-dependent affinity of helical foldamers for a rodlike guest. , 2001, Journal of the American Chemical Society.