Conformational and Reaction Dynamic Coupling in Histidine Kinases: Insights from Hybrid QM/MM Simulations

Structural and biophysical characterization of the cytoplasmic domains of HprS kinase and its interactions with the cognate regulator HprR

The HprSR constitutes the bacterial two-component regulatory system engaged by Escherichia coli to reduce the damaging effects of reactive chlorine and oxygen species present in its cytosol. Hypochlorous acid (HOCl) has been shown to be the molecule capable of activating of the HprSR system. HOCl is produced upon pathogen invasion by phagocytic cells of the human innate immune system, particularly neutrophils, to take advantage of its powerful antimicrobial attributes. Therefore, comprehensive studies concerning bacterial sensing and regulatory HprSR system are indispensable in understanding and effectively eliminating pathogens. Here we present the first crystal structure, solved at 1.7 Å resolution, of the HprS cytoplasmic domains arranged as a homodimer. In both protomers, the catalytic ATP-binding domain contains a non-hydrolysable ATP analog coordinated by a magnesium ion. This structure allowed us to provide a detailed characterization of kinase-substrate interaction. Furthermore, the structural data are supported by biophysical studies of kinase interaction with cognate response regulator HprR and substrate ATP. The kinase activity is also assessed in the presence or absence of HprR.

Mechanism of activation and autophosphorylation of a histidine kinase

Histidine kinases (HK) are one of the main prokaryotic signaling systems. Two structurally conserved catalytic domains inside the HK enable autokinase, phosphotransfer, and phosphatase activities. Here, we focus on a detailed mechanistic understanding of the functional cycle of the WalK HK by a multi-scale simulation approach, consisting of classical as well as hybrid QM/MM molecular dynamics simulation. Strikingly, a conformational transition induced solely in DHp leads to the correct activated conformation in CA crucial for autophosphorylation. This finding explains how variable sensor domains induce the transition from inactive to active state. The subsequent autophosphorylation inside DHp proceeds via a penta-coordinated transition state to a protonated phosphohistidine intermediate. This intermediate is consequently deprotonated by a suitable nearby base. The reaction energetics are controlled by the final proton acceptor and presence of a magnesium cation. The slow rates of the process result from the high energy barrier of the conformational transition between inactive and active states. The phosphorylation step exhibits a lower barrier and down-the-hill energetics. Thus, our work suggests a detailed mechanistic model for HK autophosphorylation.

An allosteric switch ensures efficient unidirectional information transmission by the histidine kinase DesK from Bacillus subtilis

Phosphorylation carries chemical information in biological systems. In two-component systems (TCSs), the sensor histidine kinase and the response regulator are connected through phosphoryl transfer reactions that may be uni- or bidirectional. Directionality enables the construction of complex regulatory networks that optimize signal propagation and ensure the forward flow of information. We combined x-ray crystallography, hybrid quantum mechanics/molecular mechanics (QM/MM) simulations, and systems-integrative kinetic modeling approaches to study phosphoryl flow through the Bacillus subtilis thermosensing TCS DesK-DesR. The allosteric regulation of the histidine kinase DesK was critical to avoid back transfer of phosphoryl groups and futile phosphorylation-dephosphorylation cycles by isolating phosphatase, autokinase, and phosphotransferase activities. Interactions between the kinase's ATP-binding domain and the regulator's receiver domain placed the regulator in two distinct positions in the phosphotransferase and phosphatase complexes, thereby determining whether a key glutamine residue in DesK was properly situated to assist in the dephosphorylation reaction. Moreover, an energetically unfavorable phosphotransferase conformation when DesK was not phosphorylated minimized reverse phosphoryl transfer. DesR dimerization and a dissociative phosphoryl transfer reaction also enforced the direction of phosphoryl flow. Shorter or longer distances between the phosphoryl acceptor and donor residues shifted the phosphoryl transfer equilibrium by modulating the stabilizing effect of the Mg²⁺ cofactor. These mechanisms control the directionality of signal transmission and show how structure-encoded allostery stores and transmits information in signaling systems.

TAB1 binding induced p38α conformation change: an accelerated molecular dynamics simulation study

p38α mitogen-activated protein kinase (MAPK) undergoes autophosphorylation induced by the binding of TGFβ-activated kinase 1 binding protein 1 (TAB1) in myocardial ischemia. Investigation of the conformational transformations in p38α triggered by TAB1 binding is motivated by the need to find selective p38α activation inhibitors to treat myocardial ischemia. Herein, the conformational transformations of p38α were studied via all-atom accelerated molecular dynamics simulations and principal component analysis. With the binding of TAB1, the conformational changes of p38α auto-activation were characterized by the movement of the activation loop (A-loop) away from the αG helix toward the αF, αE helixes and L16-loop. In addition, a diverse intermediate state with an extensional and phosphorylated A-loop different from the transition intermediate state was explored. The conformational changes, including the A-loop alpha-structure breaking and the stronger hydrogen bond network formation, are accompanied by the extension of the A-loop and more intramolecular interactions in p38α. TAB1 correlates with other regions of p38α that are distal from the TAB1-binding site, including the A-loop, αC helix, and L16-loop, which regulates the intramolecular correlation of p38α. And, the phosphorylation further enhances the correlations between the A-loop and the other regions of p38α. The correlation results imply the regulation process of p38α conformational transformations. These findings will improve our understanding of the autophosphorylation of kinase and facilitate the development of selective inhibitors for the treatment of ischemic injury.

Specificity and Reactivity of Mycobacterium tuberculosis Serine/Threonine Kinases PknG and PknB

Mycobacterium tuberculosis (Mtb), the causative agent of Tuberculosis, has 11 eukaryotic-like serine/threonine protein kinases, which play essential roles in cell growth, signal transduction, and pathogenesis. Protein kinase G (PknG) regulates the carbon and nitrogen metabolism by phosphorylation of the glycogen accumulation regulator (GarA) protein at Thr21. Protein kinase B (PknB) is involved in cell wall synthesis and cell shape, as well as phosphorylates GarA but at Thr22. While PknG seems to be constitutively activated and recognition of GarA requires phosphorylation in its unstructured tail, PknB activation is triggered by phosphorylation of its activation loop, which allows binding of the forkhead-associated domain of GarA. In the present work, we used molecular dynamics and quantum-mechanics/molecular mechanics simulations of the catalytically competent complex and kinase activity assays to understand PknG/PknB specificity and reactivity toward GarA. Two hydrophobic residues in GarA, Val24 and Phe25, seem essential for PknG binding and allow specificity for Thr21 phosphorylation. On the other hand, phosphorylated residues in PknB bind Arg26 in GarA and regulate its specificity for Thr22. We also provide a detailed analysis of the free energy profile for the phospho-transfer reaction and show why PknG has a constitutively active conformation not requiring priming phosphorylation in contrast to PknB. Our results provide new insights into these two key enzymes relevant for Mtb and the mechanisms of serine/threonine phosphorylation in bacteria.

Microsolvation of Histidine—A Theoretical Study of Intermolecular Interactions Based on AIM and SAPT Approaches

Histidine is unique among amino acids because of its rich tautomeric properties. It participates in essential enzymatic centers, such as catalytic triads. The main aim of the study is the modeling of the change of molecular properties between the gas phase and solution using microsolvation models. We investigate histidine in its three protonation states, microsolvated with 1:6 water molecules. These clusters are studied computationally, in the gas phase and with water as a solvent (Polarizable Continuum Model, PCM) within the Density Functional Theory (DFT) framework. The structural analysis reveals the presence of intra- and intermolecular hydrogen bonds. The Atoms-in-Molecules (AIM) theory is employed to determine the impact of solvation on the charge flow within the histidine, with emphasis on the similarity of the two imidazole nitrogen atoms—topologically not equivalent, they are revealed as electronically similar due to the heterocyclic ring aromaticity. Finally, the Symmetry-Adapted Perturbation Theory (SAPT) is used to examine the stability of the microsolvation clusters. While electrostatic and exchange terms dominate in magnitude over polarization and dispersion, the sum of electrostatic and exchange term is close to zero. This makes polarization the factor governing the actual interaction energy. The most important finding of this study is that even with microsolvation, the polarization induced by the presence of implicit solvent is still significant. Therefore, we recommend combined approaches, mixing explicit water molecules with implicit models.