A novel approach for large-scale polypeptide folding based on elastic networks using continuous optimization.

We present a new computationally efficient method for large-scale polypeptide folding using coarse-grained elastic networks and gradient-based continuous optimization techniques. The folding is governed by minimization of energy based on Miyazawa-Jernigan contact potentials. Using this method we are able to substantially reduce the computation time on ordinary desktop computers for simulation of polypeptide folding starting from a fully unfolded state. We compare our results with available native state structures from Protein Data Bank (PDB) for a few de-novo proteins and two natural proteins, Ubiquitin and Lysozyme. Based on our simulations we are able to draw the energy landscape for a small de-novo protein, Chignolin. We also use two well known protein structure prediction software, MODELLER and GROMACS to compare our results. In the end, we show how a modification of normal elastic network model can lead to higher accuracy and lower time required for simulation.

[1]  G. Chirikjian,et al.  Elastic models of conformational transitions in macromolecules. , 2002, Journal of molecular graphics & modelling.

[2]  Andrej ⩽ali,et al.  Comparative protein modeling by satisfaction of spatial restraints , 1995 .

[3]  Burak Erman,et al.  Gaussian model of protein folding , 2000 .

[4]  Burak Erman,et al.  The elastic net algorithm and protein structure prediction , 2002, J. Comput. Chem..

[5]  K. Dill,et al.  From Levinthal to pathways to funnels , 1997, Nature Structural Biology.

[6]  J. Onuchic,et al.  Protein folding funnels: a kinetic approach to the sequence-structure relationship. , 1992, Proceedings of the National Academy of Sciences of the United States of America.

[7]  M. Levitt,et al.  Funnel sculpting for in silico assembly of secondary structure elements of proteins , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[8]  Gideon Schreiber,et al.  Folding and binding: an extended family business , 2005 .

[9]  Robert L. Jernigan,et al.  Dynamics of large proteins through hierarchical levels of coarse‐grained structures , 2002, J. Comput. Chem..

[10]  T. Blundell,et al.  Comparative protein modelling by satisfaction of spatial restraints. , 1993, Journal of molecular biology.

[11]  G J Williams,et al.  The Protein Data Bank: a computer-based archival file for macromolecular structures. , 1978, Archives of biochemistry and biophysics.

[12]  H. Scheraga,et al.  Packing helices in proteins by global optimization of a potential energy function , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[13]  R. Jernigan,et al.  Anisotropy of fluctuation dynamics of proteins with an elastic network model. , 2001, Biophysical journal.

[14]  G. Schulz,et al.  The Covalent Structure of Proteins , 1979 .

[15]  R. Jernigan,et al.  Residue-residue potentials with a favorable contact pair term and an unfavorable high packing density term, for simulation and threading. , 1996, Journal of molecular biology.

[16]  M Karplus,et al.  Protein folding dynamics: The diffusion‐collision model and experimental data , 1994, Protein science : a publication of the Protein Society.

[17]  T. Creighton,et al.  Protein Folding , 1992 .

[18]  D. van der Spoel,et al.  GROMACS: A message-passing parallel molecular dynamics implementation , 1995 .

[19]  G J Williams,et al.  The Protein Data Bank: a computer-based archival file for macromolecular structures. , 1977, Journal of molecular biology.

[20]  I. Bahar,et al.  Coarse-grained normal mode analysis in structural biology. , 2005, Current opinion in structural biology.

[21]  Gregory S. Chirikjian,et al.  Normal mode analysis of proteins: a comparison of rigid cluster modes with Cα coarse graining , 2004 .

[22]  Shinya Honda,et al.  10 residue folded peptide designed by segment statistics. , 2004, Structure.

[23]  K. Dill,et al.  Iterative assembly of helical proteins by optimal hydrophobic packing. , 2008, Structure.

[24]  David de Sancho,et al.  Evolutionary method for the assembly of rigid protein fragments , 2005, J. Comput. Chem..

[25]  K. Dill Theory for the folding and stability of globular proteins. , 1985, Biochemistry.

[26]  I. Bahar,et al.  Gaussian Dynamics of Folded Proteins , 1997 .

[27]  K Yue,et al.  Constraint-based assembly of tertiary protein structures from secondary structure elements. , 2000, Protein science : a publication of the Protein Society.

[28]  R. Cook,et al.  Concepts and Applications of Finite Element Analysis , 1974 .

[29]  Robert L Jernigan,et al.  Rigid-cluster models of conformational transitions in macromolecular machines and assemblies. , 2005, Biophysical journal.

[30]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[31]  R. Jernigan,et al.  Inter-residue potentials in globular proteins and the dominance of highly specific hydrophilic interactions at close separation. , 1997, Journal of molecular biology.

[32]  Georg E. Schulz,et al.  Principles of Protein Structure , 1979 .

[33]  R. Jernigan,et al.  Equilibrium states of rigid bodies with multiple interaction sites: Application to protein helices , 1997 .