Self‐guided Langevin dynamics study of regulatory interactions in NtrC

Multiple self‐guided Langevin dynamics (SGLD) simulations were performed to examine structural and dynamical properties of the receiver domain of nitrogen regulatory protein C (NtrCr). SGLD and MD simulations of the phosphorylated active form structure suggest a mostly stable but broad structural ensemble of this protein. The finite difference Poisson–Boltzmann calculations of the pKa values of the active site residues suggest an increase in the pKa of His‐84 on phosphorylation of Asp‐54. In SGLD simulations of the phosphorylated active form with charged His‐84, the average position of the regulatory helix α4 is found closer to the starting structure than in simulations with the neutral His‐84. To model the transition pathway, the phosphate group was removed from the simulations. After 7 ns of simulations, the regulatory helix α4 was found approximately halfway between positions in the NMR structures of the active and inactive forms. Removal of the phosphate group stimulated loss of helix α4, suggesting that the pathway of conformational transition may involve partial unfolding mechanism. The study illustrates the potential utility of the SGLD method in studies of the coupling between ligand binding and conformational transitions. Proteins 2009. © 2009 Wiley‐Liss, Inc.

[1]  Paula Rahal,et al.  Molecular models of NS3 protease variants of the Hepatitis C virus , 2005, BMC Structural Biology.

[2]  M Karplus,et al.  Improving the accuracy of protein pKa calculations: Conformational averaging versus the average structure , 1998, Proteins.

[3]  Bernard R. Brooks,et al.  Interfacing Q‐Chem and CHARMM to perform QM/MM reaction path calculations , 2007, J. Comput. Chem..

[4]  David D. Thomas,et al.  Molecular dynamics simulations reveal a disorder-to-order transition on phosphorylation of smooth muscle myosin. , 2007, Biophysical journal.

[5]  Ron Elber,et al.  Long-timescale simulation methods. , 2005, Current opinion in structural biology.

[6]  Ray Luo,et al.  Enhanced ab initio protein folding simulations in Poisson-Boltzmann molecular dynamics with self-guiding forces. , 2004, Journal of molecular graphics & modelling.

[7]  E. Nogales,et al.  The structural basis for regulated assembly and function of the transcriptional activator NtrC. , 2006, Genes & development.

[8]  G. Vriend,et al.  Optimizing the hydrogen‐bond network in Poisson–Boltzmann equation‐based pKa calculations , 2001, Proteins.

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

[10]  A. Ninfa,et al.  Protein phosphorylation and regulation of adaptive responses in bacteria. , 1989, Microbiological reviews.

[11]  W. Kabsch,et al.  How good are predictions of protein secondary structure? , 1983, FEBS letters.

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

[13]  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.

[14]  J. Gibrat,et al.  Protein secondary structure assignment revisited: a detailed analysis of different assignment methods , 2005, BMC Structural Biology.

[15]  Ming S. Liu,et al.  Coarse‐grained dynamics of the receiver domain of NtrC: Fluctuations, correlations and implications for allosteric cooperativity , 2008, Proteins.

[16]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[17]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[18]  S. Pantano,et al.  The role of phosphorylation on the structure and dynamics of phospholamban: A model from molecular simulations , 2006, Proteins.

[19]  P. Cohen,et al.  The origins of protein phosphorylation , 2002, Nature Cell Biology.

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

[21]  A. Warshel,et al.  What are the dielectric “constants” of proteins and how to validate electrostatic models? , 2001, Proteins.

[22]  L. Johnson,et al.  Structural basis for control by phosphorylation. , 1997, Chemical reviews.

[23]  B. Brooks,et al.  Self-guided Langevin dynamics simulation method , 2003 .

[24]  Q. Cui,et al.  Reconciling the “old” and “new” views of protein allostery: A molecular simulation study of chemotaxis Y protein (CheY) , 2006, Proteins.

[25]  Shawn T. Brown,et al.  Advances in methods and algorithms in a modern quantum chemistry program package. , 2006, Physical chemistry chemical physics : PCCP.

[26]  B. Brooks,et al.  Constant pressure molecular dynamics simulation: The Langevin piston method , 1995 .

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

[28]  M. Gilson,et al.  The determinants of pKas in proteins. , 1996, Biochemistry.

[29]  Bernard R Brooks,et al.  Backbone relaxation coupled to the ionization of internal groups in proteins: a self-guided Langevin dynamics study. , 2008, Biophysical journal.

[30]  Eaton E Lattman,et al.  Experimental pK(a) values of buried residues: analysis with continuum methods and role of water penetration. , 2002, Biophysical journal.

[31]  L. R. Scott,et al.  Electrostatics and diffusion of molecules in solution: simulations with the University of Houston Brownian dynamics program , 1995 .

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

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

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

[35]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[36]  D. Wemmer,et al.  Beryllofluoride mimics phosphorylation of NtrC and other bacterial response regulators. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

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

[38]  W. L. Jorgensen,et al.  The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. , 1988, Journal of the American Chemical Society.

[39]  D. Wemmer,et al.  High-resolution solution structure of the beryllofluoride-activated NtrC receiver domain. , 2003, Biochemistry.