Assessing protein loop flexibility by hierarchical Monte Carlo sampling.

Loop flexibility is often crucial to protein biological function in solution. We report a new Monte Carlo method for generating conformational ensembles for protein loops and cyclic peptides. The approach incorporates the triaxial loop closure method which addresses the inverse kinematic problem for generating backbone move sets that do not break the loop. Sidechains are sampled together with the backbone in a hierarchical way, making it possible to make large moves that cross energy barriers. As an initial application, we apply the method to the flexible loop in triosephosphate isomerase that caps the active site, and demonstrate that the resulting loop ensembles agree well with key observations from previous structural studies. We also demonstrate, with 3 other test cases, the ability to distinguish relatively flexible and rigid loops within the same protein.

[1]  R. Friesner,et al.  Generalized Born Model Based on a Surface Integral Formulation , 1998 .

[2]  Arvin C Dar,et al.  Small molecule recognition of c-Src via the Imatinib-binding conformation. , 2008, Chemistry & biology.

[3]  Jack Snoeyink,et al.  Probik: Protein Backbone Motion by Inverse Kinematics , 2005 .

[4]  B. Berne,et al.  Multiple "time step" Monte Carlo , 2002 .

[5]  Chaok Seok,et al.  A kinematic view of loop closure , 2004, J. Comput. Chem..

[6]  Hagai Meirovitch,et al.  Efficiency of the Local Torsional Deformations Method for Identifying the Stable Structures of Cyclic Molecules , 1997 .

[7]  Jerome Nilmeier,et al.  Multiscale Monte Carlo Sampling of Protein Sidechains: Application to Binding Pocket Flexibility. , 2008, Journal of chemical theory and computation.

[8]  Thierry Siméon,et al.  Geometric algorithms for the conformational analysis of long protein loops , 2004, J. Comput. Chem..

[9]  H. Meirovitch.  Calculation of entropy with computer simulation methods , 1977 .

[10]  D. Theodorou,et al.  A concerted rotation algorithm for atomistic Monte Carlo simulation of polymer melts and glasses , 1993 .

[11]  M Karplus,et al.  Structure of the triosephosphate isomerase-phosphoglycolohydroxamate complex: an analogue of the intermediate on the reaction pathway. , 1991, Biochemistry.

[12]  William L. Jorgensen,et al.  Monte Carlo backbone sampling for polypeptides with variable bond angles and dihedral angles using concerted rotations and a Gaussian bias , 2003 .

[13]  D. Wishart,et al.  Rapid and accurate calculation of protein 1H, 13C and 15N chemical shifts , 2003, Journal of Biomolecular NMR.

[14]  J. Douglas Faires,et al.  Numerical Analysis , 1981 .

[15]  Hagai Meirovitch,et al.  Absolute free energy and entropy of a mobile loop of the enzyme acetylcholinesterase. , 2009, The journal of physical chemistry. B.

[16]  T Schlick,et al.  The loop opening/closing motion of the enzyme triosephosphate isomerase. , 1998, Biophysical journal.

[17]  C. Tanford.  Macromolecules , 1994, Nature.

[18]  Ronald M. Levy,et al.  The SGB/NP hydration free energy model based on the surface generalized born solvent reaction field and novel nonpolar hydration free energy estimators , 2002, J. Comput. Chem..

[19]  Ernst-Walter Knapp,et al.  Polypeptide folding with off-lattice Monte Carlo dynamics: the method , 1996, European Biophysics Journal.

[20]  Quincy Teng,et al.  Protein Structure Determination from NMR Data , 2013 .

[21]  E. Coutsias,et al.  Sub-angstrom accuracy in protein loop reconstruction by robotics-inspired conformational sampling , 2009, Nature Methods.

[22]  M. Deem,et al.  Analytical rebridging Monte Carlo: Application to cis/trans isomerization in proline-containing, cyclic peptides , 1999, physics/9904057.

[23]  A. Sali,et al.  Modeling of loops in protein structures , 2000, Protein science : a publication of the Protein Society.

[24]  Jerome Nilmeier,et al.  Monte Carlo Sampling with Hierarchical Move Sets: POSH Monte Carlo. , 2009, Journal of chemical theory and computation.

[25]  David C. Jones,et al.  Progress in protein structure prediction. , 1997, Current opinion in structural biology.

[26]  Aaron R. Dinner,et al.  Local deformations of polymers with nonplanar rigid main‐chain internal coordinates , 2000 .

[27]  M Karplus,et al.  Anatomy of a conformational change: hinged "lid" motion of the triosephosphate isomerase loop. , 1990, Science.

[28]  B. Honig,et al.  A hierarchical approach to all‐atom protein loop prediction , 2004, Proteins.

[29]  Adrian A Canutescu,et al.  Cyclic coordinate descent: A robotics algorithm for protein loop closure , 2003, Protein science : a publication of the Protein Society.

[30]  Michael J. Cahill,et al.  On the kinematics of protein folding , 2002, J. Comput. Chem..

[31]  宁北芳,et al.  疟原虫var基因转换速率变化导致抗原变异[英]／Paul H, Robert P, Christodoulou Z, et al//Proc Natl Acad Sci U S A , 2005 .

[32]  J. Xu,et al.  Diversity in the CDR3 region of V(H) is sufficient for most antibody specificities. , 2000, Immunity.

[33]  Steven Hayward,et al.  The effect of end constraints on protein loop kinematics. , 2010, Biophysical journal.

[34]  Harianto Tjong,et al.  Spontaneous conformational change and toxin binding in α7 acetylcholine receptor: Insight into channel activation and inhibition , 2008, Proceedings of the National Academy of Sciences.

[35]  Sergio E. Wong,et al.  Competition between intramolecular hydrogen bonds and solvation in phosphorylated peptides: simulations with explicit and implicit solvent. , 2005, The journal of physical chemistry. B.

[36]  Michael W. Deem,et al.  A configurational bias Monte Carlo method for linear and cyclic peptides , 1996, cond-mat/9709330.

[37]  H. Scheraga,et al.  Exact analytical loop closure in proteins using polynomial equations , 1999 .

[38]  F. Young.  Biochemistry , 1955, The Indian Medical Gazette.

[39]  C. Rapp,et al.  Crystal packing effects on protein loops , 2005, Proteins.

[40]  F. C. Hartman,et al.  Structure of yeast triosephosphate isomerase at 1.9-A resolution. , 1990, Biochemistry.

[41]  G. Favrin,et al.  Monte Carlo update for chain molecules: Biased Gaussian steps in torsional space , 2001, cond-mat/0103580.

[42]  Michael W. Deem,et al.  Efficient Monte Carlo methods for cyclic peptides , 1999 .

[43]  Robert E. Bruccoleri,et al.  Chain closure with bond angle variations , 1985 .

[44]  Victor H Hernandez,et al.  Nature Methods , 2007 .

[45]  Chih-Cheng Chen,et al.  A combined optimization method for solving the inverse kinematics problems of mechanical manipulators , 1991, IEEE Trans. Robotics Autom..

[46]  Francesca Massi,et al.  Solution NMR and computer simulation studies of active site loop motion in triosephosphate isomerase. , 2006, Biochemistry.

[47]  O. Brock,et al.  A methodology for efficiently sampling the conformation space of molecular structures , 2005, Physical Biology.

[48]  N. Go,et al.  Ring Closure and Local Conformational Deformations of Chain Molecules , 1970 .

[49]  A. McDermott,et al.  Triosephosphate isomerase: 15N and 13C chemical shift assignments and conformational change upon ligand binding by magic-angle spinning solid-state NMR spectroscopy. , 2010, Journal of molecular biology.

[50]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[51]  K. Dill,et al.  Resultants and Loop Closure , 2006 .

[52]  R. Friesner,et al.  Evaluation and Reparametrization of the OPLS-AA Force Field for Proteins via Comparison with Accurate Quantum Chemical Calculations on Peptides† , 2001 .

[53]  Michael Czisch,et al.  Efficient methods for obtaining phase-sensitive gradient-enhanced HMQC spectra , 1993 .

[54]  P. Tait.  Vector Analysis , 1893, Nature.

[55]  Matthew P Jacobson,et al.  Conformational selection in silico: Loop latching motions and ligand binding in enzymes , 2008, Proteins.

[56]  E. Padlan,et al.  X-ray crystallography of antibodies. , 1996, Advances in protein chemistry.

[57]  Jean-Claude Latombe,et al.  On the structure of the inverse kinematics map of a fragment of protein backbone , 2008, J. Comput. Chem..

[58]  A. McDermott,et al.  Dynamics of the flexible loop of triosephosphate isomerase: the loop motion is not ligand gated. , 1995, Biochemistry.

[59]  G A Petsko,et al.  Crystallographic analysis of the complex between triosephosphate isomerase and 2-phosphoglycolate at 2.5-A resolution: implications for catalysis. , 1990, Biochemistry.