Domains are the three-dimensional building blocks of proteins. An individual domain can be found in a variety of protein architectures that perform unique functions and are subject to different evolutionary selective pressures. We describe an approach to evaluate the variability in amino acid sequences of a single domain across architectural contexts. The ability to distinguish the different evolutionary paths of one protein domain can help determine whether existing knowledge about a specific domain will apply to an uncharacterized protein. Such discrimination can lead to insights and hypotheses about function, as well as guide experimental priorities.We developed and tested our approach on CheW-like domains (PF01584), which mediate protein/protein interactions and are difficult to compare experimentally. CheW-like domains occur in CheW scaffolding proteins, CheA kinases, and CheV proteins that regulate bacterial chemotaxis. We chose 16 protein Architectures that included 94% of all CheW-like domains found in nature. Because some Architectures had more than one CheW-like domain, CheW-like domains were divided into 21 distinct Contexts. The CheW-like domain sequences were closely related within most Contexts; however, one Context was subdivided into three Types. The resulting 23 sequence Types coalesced into five or six Classes of CheW-like domains, which we described in detail.In addition, we created SimpLogo, an innovative method for visualizing amino acid composition across large sets of multiple sequence alignments of arbitrary length. SimpLogo offers substantial advantages over standard sequence logos for comparison and analysis of related protein sequences. The R package for SimpLogo is freely available.

ABSTRACT Microorganisms and plants utilize two-component systems to regulate adaptive responses to changing environmental conditions. Sensor kinases detect stimuli and alter their autophosphorylation activity accordingly. Signal propagation occurs via the transfer of phosphoryl groups from upstream kinases to downstream response regulator proteins. Removal of phosphoryl groups from the response regulator typically resets the system. Members of the same protein family may catalyze phosphorylation and dephosphorylation reactions with different efficiencies, exhibiting rate constants spanning many orders of magnitude to accommodate response time scales from milliseconds to days. We previously found that variable positions one or two residues to the C-terminal side of the conserved Asp phosphorylation site (D+2) or Thr/Ser (T+1/T+2) in response regulators alter reaction kinetics by direct interaction with phosphodonor or phosphoacceptor molecules. Here, we explore the kinetic effects of amino acid substitutions at the two positions immediately C-terminal to the conserved Lys (K+1/K+2) in the model Escherichia coli response regulator CheY. We measured CheY autophosphorylation and autodephosphorylation rate constants for 27 pairs of K+1/K+2 residues that represent 84% of naturally occurring response regulators. Effects on autodephosphorylation were modest, but autophosphorylation rate constants varied by 2 orders of magnitude, suggesting that the K+1/K+2 positions influence reaction kinetics by altering the conformational spectrum sampled by CheY at equilibrium. Additional evidence supporting this indirect mechanism includes the following: the effect on autophosphorylation rate constants is independent of the phosphodonor, the autophosphorylation rate constants and dissociation constants for the phosphoryl group analog BeF3− are inversely correlated, and the K+1/K+2 positions are distant from the phosphorylation site. IMPORTANCE We have identified five variable positions in response regulators that allow the rate constants of autophosphorylation and autodephosporylation reactions each to be altered over 3 orders of magnitude in CheY. The distributions of variable residue combinations across response regulator subfamilies suggest that distinct mechanisms associated with different variable positions allow reaction rates to be tuned independently during evolution for diverse biological purposes. This knowledge could be used in synthetic-biology applications to adjust the properties (e.g., background noise and response duration) of biosensors and may allow prediction of response regulator reaction kinetics from the primary amino acid sequence.

Azorhizobium caulinodans is a nitrogen-fixing bacterium that forms root nodules on its host legume, Sesbania rostrata . This agriculturally significant symbiotic relationship is important in lowland rice cultivation, and allows for nitrogen fixation under flood conditions. Chemotaxis plays an important role in bacterial colonization of the rhizosphere. Plant roots release chemical compounds that are sensed by bacteria, triggering chemotaxis along a concentration gradient toward the roots. This gives motile bacteria a significant competitive advantage during root surface colonization. Although plant-associated bacterial genomes often encode multiple chemotaxis systems, A. caulinodans appears to encode only one. The che cluster on the A. caulinodans genome contains cheA , cheW , cheY2 , cheB , and cheR . Two other chemotaxis genes, cheY1 and cheZ , are located independently from the che operon. Both CheY1 and CheY2 are involved in chemotaxis, with CheY1 being the predominant signaling protein. A. caulinodans CheA contains an unusual set of C-terminal domains: a CheW-like/Receiver pair (termed W2-Rec), follows the more common single CheW-like domain. W2-Rec impacts both chemotaxis and CheA function. We found a preference for transfer of phosphoryl groups from CheA to CheY2, rather than to W2-Rec or CheY1, which appears to be involved in flagellar motor binding. Furthermore, we observed increased phosphoryl group stabilities on CheY1 compared to CheY2 or W2-Rec. Finally, CheZ enhanced dephosphorylation of CheY2 substantially more than CheY1, but had no effect on the dephosphorylation rate of W2-Rec. This network of phosphotransfer reactions highlights a previously uncharacterized scheme for regulation of chemotactic responses. IMPORTANCE Chemotaxis allows bacteria to move towards nutrients and away from toxins in their environment. Chemotactic movement provides a competitive advantage over non-specific motion. CheY is an essential mediator of the chemotactic response with phosphorylated and unphosphorylated forms of CheY differentially interacting with the flagellar motor to change swimming behavior. Previously established schemes of CheY dephosphorylation include action of a phosphatase and/or transfer of the phosphoryl group to another receiver domain that acts as a sink. Here, we propose A. caulinodans uses a concerted mechanism in which the Hpt domain of CheA, CheY2, and CheZ function together as a dual sink system to rapidly reset chemotactic signaling. To the best of our knowledge, this mechanism is unlike any that have previously been evaluated. Chemotaxis systems that utilize both receiver and Hpt domains as phosphate sinks likely occur in other bacterial species.

Domains are the three-dimensional building blocks of proteins. An individual domain can occur in a variety of domain architectures that perform unique functions and are subject to different evolutionary selective pressures. We describe an approach to evaluate the variability in amino acid sequences of a single domain across architectural contexts. The ability to distinguish different evolutionary outcomes of one protein domain can help determine whether existing knowledge about a specific domain will apply to an uncharacterized protein, lead to insights and hypotheses about function, and guide experimental priorities. We developed and tested our approach on CheW-like domains (PF01584), which mediate protein/protein interactions and are difficult to compare experimentally. CheW-like domains occur in CheW scaffolding proteins, CheA kinases, and CheV proteins that regulate bacterial chemotaxis. We analyzed 16 domain architectures that included 94% of all CheW-like domains found in nature. We identified six Classes of CheW-like domains with presumed functional differences. CheV and most CheW proteins contained Class 1 domains, whereas some CheW proteins contained Class 6 ($20%) or Class 2 ($1%) domains instead.Most CheA proteins contained Class 3 domains. CheA proteins with multiple Hpt domains contained Class 4 domains. CheA proteins with two CheW-like domains contained one Class 3 and one Class 5. We also created SimpLogo, an innovative method for visualizing amino acid composition across large sets of multiple sequence alignments of arbitrary length. SimpLogo offers substantial advantages over standard sequence logos for comparison and analysis of related protein sequences. The R package for SimpLogo is freely available.

No abstract

Two-component signaling is a primary method by which microorganisms interact with their environments. A kinase detects stimuli and modulates autophosphorylation activity. The signal propagates by phosphotransfer from the kinase to a response regulator, eliciting a response. Response regulators operate over a range of time scales, corresponding to their related biological processes. Response regulator active site chemistry is highly conserved, but certain variable residues can influence phosphorylation kinetics. An Ala-to-Pro substitution (K+4, residue 113) in the Escherichia coli response regulator CheY triggers a constitutively active phenotype; however, the A113P substitution is too far from the active site to directly affect phosphochemistry. To better understand the activating mechanism(s) of the substitution, we analyzed receiver domain sequences to characterize the evolutionary role of the K+4 position. Although most featured Pro, Leu, Ile, and Val residues, chemotaxis-related proteins exhibited atypical Ala, Gly, Asp, and Glu residues at K+4. Structural and in silico analyses revealed that CheY A113P adopted a partially active configuration. Biochemical data showed that A113P shifted CheY toward a more activated state, enhancing autophosphorylation. By characterizing CheY variants, we determined that this functionality was transmitted through a hydrophobic network bounded by the β5α5 loop and the α1 helix of CheY. This region also interacts with the phosphodonor CheA P1 , suggesting that binding generates an activating perturbation similar to the A113P substitution. Atypical residues like Ala at the K+4 position likely serve two purposes. First, restricting autophosphorylation may minimize background noise generated by intracellular phosphodonors such as acetyl phosphate. Second, optimizing interactions with upstream partners may help prime the receiver domain for phosphorylation.

The data that support the findings of this study were derived from the following resources available in the public domain: Pfam (release v.33), http://pfam.xfam.org; MiST, https://mistdb.com. The data we compiled are available as Datasets S1 through S4 in the Supporting Information of this article.

The diacylglycerol kinases (DGKs) are a family of enzymes responsible for the conversion of diacylglycerol (DAG) to phosphatidic acid (PA). In addition to their primary function in lipid metabolism, DGKs have recently been identified as potential therapeutic targets in multiple cancers, including glioblastoma (GBM) and melanoma. Aside from its tumorigenic properties, DGKα is also a known promoter of T-cell anergy, supporting a role as a recently-recognized T cell checkpoint. In fact, the only significant phenotype previously observed in Dgka knockout (KO) mice is the enhancement of T-cell activity. Herein we reveal a novel, macrophage-specific, immune-regulatory function of DGKα. In bone marrow-derived macrophages (BMDMs) cultured from wild-type (WT) and KO mice, we observed increased responsiveness of KO macrophages to diverse stimuli that yield different phenotypes, including LPS, IL-4, and the chemoattractant MCP-1. Knockdown (KD) of Dgka in a murine macrophage cell line resulted in similar increased responsiveness. Demonstrating in vivo relevance, we observed significantly smaller wounds in Dgka-/- mice with full-thickness cutaneous burns, a complex wound healing process in which macrophages play a key role. The burned area also demonstrated increased numbers of macrophages. In a cortical stab wound model, Dgka-/- brains show increased Iba1+ cell numbers at the needle track versus that in WT brains. Taken together, these findings identify a novel immune-regulatory checkpoint function of DGKα in macrophages with potential implications for wound healing, cancer therapy, and other settings.