Free energy landscape of activation in a signaling protein at atomic resolution

The interconversion between inactive and active protein states, traditionally described by two static structures, is at the heart of signaling. However, how folded states interconvert is largely unknown due to the inability to experimentally observe transition pathways. Here we explore the free energy landscape of the bacterial response regulator NtrC by combining computation and NMR, and discover unexpected features underlying efficient signaling. We find that functional states are defined purely in kinetic and not structural terms. The need of a well-defined conformer, crucial to the active state, is absent in the inactive state, which comprises a heterogeneous collection of conformers. The transition between active and inactive states occurs through multiple pathways, facilitated by a number of nonnative transient hydrogen bonds, thus lowering the transition barrier through both entropic and enthalpic contributions. These findings may represent general features for functional conformational transitions within the folded state.

[1]  Frank Noé,et al.  Crystal structure of nucleotide-free dynamin , 2011, Nature.

[2]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[3]  G. Ciccotti,et al.  String method in collective variables: minimum free energy paths and isocommittor surfaces. , 2006, The Journal of chemical physics.

[4]  Donald Hilvert,et al.  Precision is essential for efficient catalysis in an evolved Kemp eliminase , 2013, Nature.

[5]  R. Dror,et al.  Systematic Validation of Protein Force Fields against Experimental Data , 2012, PloS one.

[6]  D. Kern,et al.  Evidence against the "Y-T coupling" mechanism of activation in the response regulator NtrC. , 2014, Journal of molecular biology.

[7]  S. L. Mayo,et al.  Enzyme-like proteins by computational design , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[8]  S. Grzesiek,et al.  A simple and sensitive experiment for measurement of JCC couplings between backbone carbonyl and methyl carbons in isotopically enriched proteins , 1993, Journal of Biomolecular NMR.

[9]  Vijay S. Pande,et al.  Screen Savers of the World Unite! , 2000, Science.

[10]  Frank Noé,et al.  An Introduction to Markov State Models and Their Application to Long Timescale Molecular Simulation , 2014, Advances in Experimental Medicine and Biology.

[11]  Peter G. Bolhuis,et al.  Predicting the reaction coordinates of millisecond light-induced conformational changes in photoactive yellow protein , 2010, Proceedings of the National Academy of Sciences.

[12]  Bernard R Brooks,et al.  Self‐guided Langevin dynamics study of regulatory interactions in NtrC , 2009, Proteins.

[13]  Dong Wang,et al.  A Two-State Model for the Dynamics of the Pyrophosphate Ion Release in Bacterial RNA Polymerase , 2013, PLoS Comput. Biol..

[14]  P. Wolynes,et al.  Conformational switching upon phosphorylation: a predictive framework based on energy landscape principles. , 2008, Biochemistry.

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

[16]  Diwakar Shukla,et al.  Activation pathway of Src kinase reveals intermediate states as targets for drug design , 2014 .

[17]  D E Wemmer,et al.  Two-state allosteric behavior in a single-domain signaling protein. , 2001, Science.

[18]  Albert C. Pan,et al.  Finding transition pathways using the string method with swarms of trajectories. , 2008, The journal of physical chemistry. B.

[19]  D E Wemmer,et al.  Three-dimensional solution structure of the N-terminal receiver domain of NTRC. , 1995, Biochemistry.

[20]  C. Schütte,et al.  Supplementary Information for “ Constructing the Equilibrium Ensemble of Folding Pathways from Short Off-Equilibrium Simulations ” , 2009 .

[21]  Vijay S Pande,et al.  Protein folded states are kinetic hubs , 2010, Proceedings of the National Academy of Sciences.

[22]  Bruce A. Johnson,et al.  NMR View: A computer program for the visualization and analysis of NMR data , 1994, Journal of biomolecular NMR.

[23]  L. Kay,et al.  Gradient-Enhanced Triple-Resonance Three-Dimensional NMR Experiments with Improved Sensitivity , 1994 .

[24]  Diwakar Shukla,et al.  A network of molecular switches controls the activation of the two-component response regulator NtrC , 2015, Nature Communications.

[25]  Michael R. Shirts,et al.  COMPUTING: Screen Savers of the World Unite! , 2000, Science.

[26]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[27]  Yang Shen,et al.  Identification of helix capping and β-turn motifs from NMR chemical shifts , 2012, Journal of biomolecular NMR.

[28]  S. Kustu,et al.  The phosphorylated form of the enhancer-binding protein NTRC has an ATPase activity that is essential for activation of transcription , 1991, Cell.

[29]  Pathways for conformational change in nitrogen regulatory protein C from discrete path sampling. , 2008, The journal of physical chemistry. B.

[30]  N. C. Price,et al.  3D HCCH3 -TOCSY for resonance assignment of methyl-containing side chains in 13C-labeled proteins , 2000 .

[31]  Eric A. Althoff,et al.  Kemp elimination catalysts by computational enzyme design , 2008, Nature.

[32]  S. Grzesiek,et al.  Two-Dimensional NMR Methods for Determining χ1 Angles of Aromatic Residues in Proteins from Three-Bond JC‘Cγ and JNCγ Couplings , 1997 .

[33]  Molecular Dynamic Simulations of the N-Terminal Receiver Domain of NtrC Reveal Intrinsic Conformational Flexibility in the Inactive State , 2006, Journal of biomolecular structure & dynamics.

[34]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[35]  L. Johnson,et al.  Active and Inactive Protein Kinases: Structural Basis for Regulation , 1996, Cell.

[36]  Vijay S Pande,et al.  Simple few-state models reveal hidden complexity in protein folding , 2012, Proceedings of the National Academy of Sciences.

[37]  K. Wüthrich,et al.  Torsion angle dynamics for NMR structure calculation with the new program DYANA. , 1997, Journal of molecular biology.

[38]  M. Klein,et al.  Constant pressure molecular dynamics algorithms , 1994 .

[39]  Xuhui Huang,et al.  Using generalized ensemble simulations and Markov state models to identify conformational states. , 2009, Methods.

[40]  M. Parrinello,et al.  Polymorphic transitions in single crystals: A new molecular dynamics method , 1981 .

[41]  Eric Vanden-Eijnden,et al.  Revisiting the finite temperature string method for the calculation of reaction tubes and free energies. , 2009, The Journal of chemical physics.

[42]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[43]  Hiroshi Fujisaki,et al.  Minimum Free Energy Path of Ligand-Induced Transition in Adenylate Kinase , 2012, PLoS Comput. Biol..

[44]  Joseph A. Bank,et al.  Supporting Online Material Materials and Methods Figs. S1 to S10 Table S1 References Movies S1 to S3 Atomic-level Characterization of the Structural Dynamics of Proteins , 2022 .

[45]  B. Roux.  The calculation of the potential of mean force using computer simulations , 1995 .

[46]  A. Bax,et al.  TALOS+: a hybrid method for predicting protein backbone torsion angles from NMR chemical shifts , 2009, Journal of biomolecular NMR.

[47]  Martin Karplus,et al.  Segmented transition pathway of the signaling protein nitrogen regulatory protein C. , 2009, Journal of molecular biology.

[48]  P. Deuflhard,et al.  A Direct Approach to Conformational Dynamics Based on Hybrid Monte Carlo , 1999 .

[49]  R. Altman,et al.  Cloud-based simulations on Google Exacycle reveal ligand-modulation of GPCR activation pathways , 2013, Nature chemistry.

[50]  L. Kay,et al.  Two-dimensional NMR experiments for correlating carbon-13.beta. and proton.delta./.epsilon. chemical shifts of aromatic residues in 13C-labeled proteins via scalar couplings , 1993 .

[51]  S. Grzesiek,et al.  Correlation of Backbone Amide and Aliphatic Side-Chain Resonances in 13C/15N-Enriched Proteins by Isotropic Mixing of 13C Magnetization , 1993 .

[52]  D. Kahn,et al.  Conformational changes induced by phosphorylation of the FixJ receiver domain. , 1999, Structure.

[53]  Brian F. Volkman,et al.  Structure of a transiently phosphorylated switch in bacterial signal transduction , 2000, Nature.

[54]  Masaki Sasai,et al.  Entropic mechanism of large fluctuation in allosteric transition , 2010, Proceedings of the National Academy of Sciences.

[55]  Oliver F. Lange,et al.  Evaluation and optimization of discrete state models of protein folding. , 2012, The journal of physical chemistry. B.

[56]  Vijay S Pande,et al.  Improvements in Markov State Model Construction Reveal Many Non-Native Interactions in the Folding of NTL9. , 2013, Journal of chemical theory and computation.

[57]  S. Grzesiek,et al.  NMRPipe: A multidimensional spectral processing system based on UNIX pipes , 1995, Journal of biomolecular NMR.

[58]  C. Brooks,et al.  Native states of fast-folding proteins are kinetic traps. , 2013, Journal of the American Chemical Society.

[59]  Jeremy C. Smith,et al.  Hierarchical analysis of conformational dynamics in biomolecules: transition networks of metastable states. , 2007, The Journal of chemical physics.

[60]  Wilfred Pinfold,et al.  Proceedings of the Conference on High Performance Computing Networking, Storage and Analysis , 2009, HiPC 2009.

[61]  Jeremy C. Smith,et al.  Transition Networks for the Comprehensive Characterization of Complex Conformational Change in Proteins. , 2006, Journal of chemical theory and computation.

[62]  Wladimir Labeikovsky,et al.  Transient Non-native Hydrogen Bonds Promote Activation of a Signaling Protein , 2009, Cell.

[63]  Frank Noé,et al.  Markov state models of biomolecular conformational dynamics. , 2014, Current opinion in structural biology.

[64]  C. Griesinger,et al.  Novel strategies for sensitivity enhancement in heteronuclear multi—dimensional NMR experiments employing pulsed field gradients , 1995, Journal of biomolecular NMR.

[65]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[66]  Frank Noé,et al.  Kinetic characterization of the critical step in HIV-1 protease maturation , 2012, Proceedings of the National Academy of Sciences.

[67]  John L. Klepeis,et al.  Millisecond-scale molecular dynamics simulations on Anton , 2009, Proceedings of the Conference on High Performance Computing Networking, Storage and Analysis.

[68]  G. Bowman,et al.  Equilibrium fluctuations of a single folded protein reveal a multitude of potential cryptic allosteric sites , 2012, Proceedings of the National Academy of Sciences.

[69]  Massimiliano Bonomi,et al.  PLUMED: A portable plugin for free-energy calculations with molecular dynamics , 2009, Comput. Phys. Commun..

[70]  Kyle A. Beauchamp,et al.  Markov state model reveals folding and functional dynamics in ultra-long MD trajectories. , 2011, Journal of the American Chemical Society.