Correlated Motions and Residual Frustration in Thrombin

Thrombin is the central protease in the cascade of blood coagulation proteases. The structure of thrombin consists of a double β-barrel core surrounded by connecting loops and helices. Compared to chymotrypsin, thrombin has more extended loops that are thought to have arisen from insertions in the serine protease that evolved to impart greater specificity. Previous experiments showed thermodynamic coupling between ligand binding at the active site and distal exosites. We present a combined approach of molecular dynamics (MD), accelerated molecular dynamics (AMD), and analysis of the residual local frustration of apo-thrombin and active-site-bound (PPACK-thrombin). Community analysis of the MD ensembles identified changes upon active site occupation in groups of residues linked through correlated motions and physical contacts. AMD simulations, calibrated on measured residual dipolar couplings, reveal that upon active site ligation, correlated loop motions are quenched, but new ones connecting the active site with distal sites where allosteric regulators bind emerge. Residual local frustration analysis reveals a striking correlation between frustrated contacts and regions undergoing slow time scale dynamics. The results elucidate a motional network that probably evolved through retention of frustrated contacts to provide facile conversion between ensembles of states.

[1]  Aric Hagberg,et al.  Exploring Network Structure, Dynamics, and Function using NetworkX , 2008, Proceedings of the Python in Science Conference.

[2]  P. Kollman,et al.  Molecular Dynamics Simulations on Solvated Biomolecular Systems: The Particle Mesh Ewald Method Leads to Stable Trajectories of DNA, RNA, and Proteins , 1995 .

[3]  James J. Feng,et al.  How malaria parasites reduce the deformability of infected red blood cells. , 2012, Biophysical journal.

[4]  Zaida Luthey-Schulten,et al.  Exploring residue component contributions to dynamical network models of allostery. , 2012, Journal of chemical theory and computation.

[5]  Peter G Wolynes,et al.  The spectrum of biomolecular states and motions , 2008, HFSP journal.

[6]  Robert Huber,et al.  The refined 1.9 A crystal structure of human alpha‐thrombin: interaction with D‐Phe‐Pro‐Arg chloromethylketone and significance of the Tyr‐Pro‐Pro‐Trp insertion segment. , 1989 .

[7]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[8]  Peter G Wolynes,et al.  Localizing frustration in native proteins and protein assemblies , 2007, Proceedings of the National Academy of Sciences.

[9]  R. Huber,et al.  The refined 1.9 A crystal structure of human alpha‐thrombin: interaction with D‐Phe‐Pro‐Arg chloromethylketone and significance of the Tyr‐Pro‐Pro‐Trp insertion segment. , 1989, The EMBO journal.

[10]  J Andrew McCammon,et al.  The dynamic structure of thrombin in solution. , 2012, Biophysical journal.

[11]  E. Komives,et al.  Allosteric changes in solvent accessibility observed in thrombin upon active site occupation. , 2004, Biochemistry.

[12]  R. Nussinov,et al.  Allostery: absence of a change in shape does not imply that allostery is not at play. , 2008, Journal of molecular biology.

[13]  Oliver F. Lange,et al.  Generalized correlation for biomolecular dynamics , 2005, Proteins.

[14]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[15]  Michael Nilges,et al.  Computational approaches to the interpretation of NMR data for studying protein dynamics , 2012 .

[16]  J. Huntington,et al.  Ligand Binding Shuttles Thrombin along a Continuum of Zymogen- and Proteinase-like States* , 2010, The Journal of Biological Chemistry.

[17]  J. Mccammon,et al.  Allosteric networks in thrombin distinguish procoagulant vs. anticoagulant activities , 2012, Proceedings of the National Academy of Sciences.

[18]  Z. Luthey-Schulten,et al.  Dynamical networks in tRNA:protein complexes , 2009, Proceedings of the National Academy of Sciences.

[19]  J. Onuchic,et al.  Theory of protein folding: the energy landscape perspective. , 1997, Annual review of physical chemistry.

[20]  J. Mccammon,et al.  Functional Dynamics of the Folded Ankyrin Repeats of IκBα Revealed by Nuclear Magnetic Resonance† , 2009, Biochemistry.

[21]  J. Mccammon,et al.  Sampling of slow diffusive conformational transitions with accelerated molecular dynamics. , 2007, The Journal of chemical physics.

[22]  Mohammad M. Sultan,et al.  Allosteric pathways in imidazole glycerol phosphate synthase , 2012, Proceedings of the National Academy of Sciences.

[23]  D. Dryden,et al.  Allostery without conformational change , 1984, European Biophysics Journal.

[24]  R. MacGillivray,et al.  Evolution of prothrombin: Isolation and characterization of the cDNAs encoding chicken and hagfish prothrombin , 1994, Journal of Molecular Evolution.

[25]  D. Kern,et al.  The role of dynamics in allosteric regulation. , 2003, Current opinion in structural biology.

[26]  C. Esmon.  Regulation of blood coagulation. , 2000, Biochimica et biophysica acta.

[27]  E. Komives,et al.  Thermodynamic compensation upon binding to exosite 1 and the active site of thrombin. , 2011, Biochemistry.

[28]  J. Onuchic,et al.  Nonlinear elasticity, proteinquakes, and the energy landscapes of functional transitions in proteins , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[29]  E. Komives,et al.  Amide H/2H exchange reveals a mechanism of thrombin activation. , 2006, Biochemistry.

[30]  Adrian Gustavo Turjanski,et al.  Protein frustratometer: a tool to localize energetic frustration in protein molecules , 2012, Nucleic Acids Res..

[31]  J Andrew McCammon,et al.  Studying functional dynamics in bio-molecules using accelerated molecular dynamics. , 2011, Physical chemistry chemical physics : PCCP.

[32]  Berk Hess,et al.  GROMACS 3.0: a package for molecular simulation and trajectory analysis , 2001 .

[33]  Charalampos G. Kalodimos,et al.  Protein dynamics and allostery: an NMR view. , 2011, Current opinion in structural biology.

[34]  Rafael Brüschweiler,et al.  Protein conformational flexibility from structure-free analysis of NMR dipolar couplings: quantitative and absolute determination of backbone motion in ubiquitin. , 2009, Angewandte Chemie.

[35]  P. Wolynes,et al.  On the role of frustration in the energy landscapes of allosteric proteins , 2011, Proceedings of the National Academy of Sciences.

[36]  J. Changeux,et al.  ON THE NATURE OF ALLOSTERIC TRANSITIONS: A PLAUSIBLE MODEL. , 1965, Journal of molecular biology.

[37]  Samuel Genheden,et al.  Starting-Condition Dependence of Order Parameters Derived from Molecular Dynamics Simulations. , 2010, Journal of chemical theory and computation.

[38]  S. Freund,et al.  NMR resonance assignments of thrombin reveal the conformational and dynamic effects of ligation , 2010, Proceedings of the National Academy of Sciences.