Leap‐dynamics: efficient sampling of conformational space of proteins and peptides in solution

A molecular simulation scheme, called Leap‐dynamics, that provides efficient sampling of protein conformational space in solution is presented. The scheme is a combined approach using a fast sampling method, imposing conformational ‘leaps’ to force the system over energy barriers, and molecular dynamics (MD) for refinement. The presence of solvent is approximated by a potential of mean force depending on the solvent accessible surface area. The method has been successfully applied to N‐acetyl‐L‐alanine‐N‐methylamide (alanine dipeptide), sampling experimentally observed conformations inaccessible to MD alone under the chosen conditions. The method predicts correctly the increased partial flexibility of the mutant Y35G compared to native bovine pancreatic trypsin inhibitor. In particular, the improvement over MD consists of the detection of conformational flexibility that corresponds closely to slow motions identified by nuclear magnetic resonance techniques.

[1]  M. Karplus,et al.  Effective energy function for proteins in solution , 1999, Proteins.

[2]  H. Hope.  Cryocrystallography of biological macromolecules: a generally applicable method. , 1988, Acta crystallographica. Section B, Structural science.

[3]  K. Wüthrich,et al.  Conformational sampling by NMR solution structures calculated with the program DIANA evaluated by comparison with long‐time molecular dynamics calculations in explicit water , 1996, Proteins.

[4]  W F van Gunsteren,et al.  A structure refinement method based on molecular dynamics in four spatial dimensions. , 1993, Journal of molecular biology.

[5]  A Kitao,et al.  Energy landscape of a native protein: Jumping‐among‐minima model , 1998, Proteins.

[6]  M. E. Karpen,et al.  Comparing short protein substructures by a method based on backbone torsion angles , 1989, Proteins.

[7]  W. V. van Gunsteren,et al.  An efficient mean solvation force model for use in molecular dynamics simulations of proteins in aqueous solution. , 1996, Journal of molecular biology.

[8]  J. Thornton,et al.  PROCHECK: a program to check the stereochemical quality of protein structures , 1993 .

[9]  P. Kollman,et al.  Pathways to a protein folding intermediate observed in a 1-microsecond simulation in aqueous solution. , 1998, Science.

[10]  T. Oas,et al.  Determinants of backbone dynamics in native BPTI: cooperative influence of the 14-38 disulfide and the Tyr35 side-chain. , 1998, Journal of molecular biology.

[11]  T. Oas,et al.  Enhanced protein flexibility caused by a destabilizing amino acid replacement in BPTI. , 1997, Journal of molecular biology.

[12]  S Doniach,et al.  Protein dynamics simulations from nanoseconds to microseconds. , 1999, Current opinion in structural biology.

[13]  C. Dobson,et al.  Comparison of MD simulations and NMR experiments for hen lysozyme. Analysis of local fluctuations, cooperative motions, and global changes. , 1995, Biochemistry.

[14]  G Vriend,et al.  WHAT IF: a molecular modeling and drug design program. , 1990, Journal of molecular graphics.

[15]  D L Caspar,et al.  Problems in simulating macromolecular movements. , 1995, Structure.

[16]  Charles L. Brooks,et al.  Theoretical study of blocked glycine and alanine peptide analogs , 1991 .

[17]  G. Vriend,et al.  Prediction of protein conformational freedom from distance constraints , 1997, Proteins.

[18]  A T Brünger,et al.  Influence of internal dynamics on accuracy of protein NMR structures: derivation of realistic model distance data from a long molecular dynamics trajectory. , 1999, Journal of molecular biology.

[19]  A. Mark,et al.  Computational approaches to study protein unfolding: Hen egg white lysozyme as a case study , 1995, Proteins.

[20]  Robert Huber,et al.  Structure of bovine pancreatic trypsin inhibitor , 1984 .

[21]  K. Kopple,et al.  Solvent-dependent conformational distributions of some dipeptides , 1980 .

[22]  Franca Fraternali,et al.  Conformational transitions of a dipeptide in water: Effects of imposed pathways using umbrella sampling techniques , 1994 .

[23]  K. Wüthrich,et al.  Structure and internal dynamics of the bovine pancreatic trypsin inhibitor in aqueous solution from long‐time molecular dynamics simulations , 1995, Proteins.

[24]  H. Berendsen,et al.  A Glimpse of the Holy Grail? , 1998, Science.

[25]  M. Billeter,et al.  MOLMOL: a program for display and analysis of macromolecular structures. , 1996, Journal of molecular graphics.

[26]  B. Montgomery Pettitt,et al.  Conformational free energy of hydration for the alanine dipeptide: Thermodynamic analysis , 1988 .

[27]  Wilfred F. van Gunsteren,et al.  SWARM-MD: Searching Conformational Space by Cooperative Molecular Dynamics , 1998 .

[28]  W. Kabsch,et al.  Dictionary of protein secondary structure: Pattern recognition of hydrogen‐bonded and geometrical features , 1983, Biopolymers.

[29]  H. Berendsen,et al.  Essential dynamics of proteins , 1993, Proteins.

[30]  A Wlodawer,et al.  Crystal structure of a Y35G mutant of bovine pancreatic trypsin inhibitor. , 1991, Journal of molecular biology.

[31]  Alan E. Mark,et al.  The GROMOS96 Manual and User Guide , 1996 .

[32]  R. Huber,et al.  Simultaneous refinement of the structure of BPTI against NMR data measured in solution and X-ray diffraction data measured in single crystals. , 1994, Journal of molecular biology.

[33]  Tamar Schlick,et al.  Long timestep dynamics of peptides by the dynamics driver approach , 1995, Proteins.