Coarse-grained dynamics of the receiver domain of NtrC: fluctuations, correlations and implications for allosteric cooperativity

Small phospho-donors phosphorylate MorR without inducing protein conformational changes

Phosphorylation is an essential mechanism of protein control and plays an important role in biology. The two-component system (TCS) is a bacterial regulation mechanism mediated by a response regulator (RR) protein and a kinase protein, which synchronize the regulatory circuit according to the environment. Phosphorylation is a key element in TCS function as it controls RR activity. In the present study, we characterize the behavior of MorR, an RR associated with Mo homeostasis, upon acetylphosphate and phosphoramidate treatment in vitro. Our results show that MorR was phosphorylated by both phospho-donors. Fluorescence experiments showed that MorR tryptophan emission is quenched by phosphoramidate. Furthermore, theoretical and computational results demonstrate that phosphorylation by phosphoramidate is more favorable than that by acetylphosphate. In conclusion, phosphorylated MorR is a monomeric protein and phosphorylation does not appear to induce observable conformational changes in the protein structure.

Molecular dynamics reveal a novel kinase-substrate interface that regulates protein translation

A key control point in gene expression is the initiation of protein translation, with a universal stress response being constituted by inhibitory phosphorylation of the eukaryotic initiation factor 2α (eIF2α). In humans, four kinases sense diverse physiological stresses to regulate eIF2α to control cell differentiation, adaptation, and survival. Here we develop a computational molecular model of eIF2α and one of its kinases, the protein kinase R, to simulate the dynamics of their interaction. Predictions generated by coarse-grained dynamics simulations suggest a novel mode of action. Experimentation substantiates these predictions, identifying a previously unrecognized interface in the protein complex, which is constituted by dynamic residues in both eIF2α and its kinases that are crucial to regulate protein translation. These findings call for a reinterpretation of the current mechanism of action of the eIF2α kinases and demonstrate the value of conducting computational analysis to evaluate protein function.

Evidence against the "Y-T coupling" mechanism of activation in the response regulator NtrC

The dominant theory on the mechanism of response regulators activation in two-component bacterial signaling systems is the "Y-T coupling" mechanism, wherein the χ1 rotameric state of a highly conserved aromatic residue correlates with the activation of the protein via structural rearrangements coupled to a conserved tyrosine. In this paper, we present evidence that, in the receiver domain of the response regulator nitrogen regulatory protein C (NtrC(R)), the interconversion of this tyrosine (Y101) between its rotameric states is actually faster than the rate of inactive/active conversion and is not correlated to the activation process. Data gathered from NMR relaxation dispersion experiments show that a subset of residues surrounding the conserved tyrosine sense a process that is occurring at a faster rate than the inactive/active conformational transition. We show that this process is related to χ1 rotamer exchange of Y101 and that mutation of this aromatic residue to a leucine eliminated this second faster process without affecting activation. Computational simulations of NtrC(R) in its active conformation further demonstrate that the rotameric state of Y101 is uncorrelated with the global conformational transition during activation. Moreover, the tyrosine does not appear to be involved in the stabilization of the active form upon phosphorylation and is not essential in propagating the signal downstream for ATPase activity of the central domain. Our data provide experimental evidence against the generally accepted "Y-T coupling" mechanism of activation in NtrC(R).

Entropic mechanism of large fluctuation in allosteric transition

A statistical mechanical model of allosteric transitions in proteins is developed by extending the structure-based model of protein folding to cases of multiple native conformations. The partition function is calculated exactly within the model and the free-energy surface reflecting allostery is derived. This approach is applied to an example protein, the receiver domain of the bacterial enhancer-binding protein NtrC. The model predicts the large entropy associated with a combinatorial number of preexisting transition routes. This large entropy lowers the free-energy barrier of the allosteric transition, which explains the large structural fluctuation observed in the NMR data of NtrC. The global allosteric transformation of NtrC is explained by the shift of preexisting distribution of conformations upon phosphorylation, but the local structural adjustment around the phosphorylation site is explained by the complementary induced-fit mechanism. Structural disordering accompanied by fluctuating interactions specific to two allosteric conformations underlies a large number of routes of allosteric transition.

Receiver domain structure and function in response regulator proteins

During signal transduction by two-component regulatory systems, sensor kinases detect and encode input information while response regulators (RRs) control output. Most receiver domains function as phosphorylation-mediated switches within RRs, but some transfer phosphoryl groups in multistep phosphorelays. Conserved features of receiver domain amino acid sequence correlate with structure and hence function. Receiver domains catalyze their own phosphorylation and dephosphorylation in reactions requiring a divalent cation. Molecular dynamics simulations are supplementing structural investigation of the conformational changes that underlie receiver domain switch function. As understanding of features shared by all receiver domains matures, factors conferring differences (e.g. in reaction rate or specificity) are receiving increased attention. Numerous examples of atypical receiver or pseudo-receiver domains that function without phosphorylation have recently been characterized.

Allosteric Conformational Transition in Adenylate Kinase: Dynamic Correlations and Implication for Allostery

An essential aspect of protein science is to determine the deductive relationship between structure, dynamics, and various sets of functions. The role of dynamics is currently challenging our understanding of protein functions, both experimentally and theoretically. To verify the internal fluctuations and dynamics correlations in an enzyme protein undergoing conformational transitions, we have applied a coarse-grained dynamics algorithm using the elastic network model for adenylate kinase. Normal mode analysis reveals possible dynamical and allosteric pathways for the transition between the open and the closed states of adenylate kinase. As the ligands binding induces significant flexibility changes of the nucleotides monophosphate (NMP) domain and adenosine triphosphate (ATP) domain, the diagonalized correlation between different structural transition states shows that most correlated motions occur between the NMP domain and the helices surrounding the ATP domain. The simultaneous existence of positive and negative correlations indicates that the conformational changes of adenylate kinase take place in an allosteric manner. Analyses of the cumulated normal mode overlap coefficients and long-range correlated motion provide new insights of operating mechanisms and dynamics of adenylate kinase. They also suggest a quantitative dynamics criterion for determining the allosteric cooperativity, which may be applicable to other proteins.

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 (NtrC(r)). 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 pK(a) values of the active site residues suggest an increase in the pK(a) 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 alpha4 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 alpha4 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 alpha4, 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.