Two-component signaling systems are ubiquitous in bacteria, Archaea and plants and play important roles in sensing and responding to environmental stimuli. To propagate a signaling response the typical system employs a sensory histidine kinase that phosphorylates a Receiver (REC) domain on a conserved aspartate (Asp) residue. Although it is known that some REC domains are missing this Asp residue, it remains unclear as to how many of these divergent REC domains exist, what their functional roles are and how they are regulated in the absence of the conserved Asp. Here we have compiled all deposited REC domains missing their phosphorylatable Asp residue, renamed here as the Aspartate-Less Receiver (ALR) domains. Our data show that ALRs are surprisingly common and are enriched for when attached to more rare effector outputs. Analysis of our informatics and the available ALR atomic structures, combined with structural, biochemical and genetic data of the ALR archetype RitR from Streptococcus pneumoniae presented here suggest that ALRs have reorganized their active pockets to instead take on a constitutive regulatory role or accommodate input signals other than Asp phosphorylation, while largely retaining the canonical post-phosphorylation mechanisms and dimeric interface. This work defines ALRs as an atypical REC subclass and provides insights into shared mechanisms of activation between ALR and REC domains.
In a changing environment organisms must have means to effectively respond to external stimuli or perish. In bacteria, although so-called One-Component Signaling systems numerically dominate microbial genomes [1], Two-Component Signaling (TCS) systems also play crucial roles in adaptation to the changing external environment [2,3]. Classical TCS systems consist of a membrane-bound histidine kinase that upon sensing an external stimulus, autophosphorylates using ATP as a phosphate donor and then subsequently transfers the phosphate to a conserved Asp residue within a cognate Receiver (REC) domain-containing partner in the cytoplasm. REC domains adopt the flavodoxin fold, a very common (α/β) 5 architecture with a central 5-parallell β-sheet. The typical REC active pocket structure consists of an ‘acidic triad’ that includes the phospho-accepting Asp residue at the end of the β3 strand and two other acidic amino acids within the β1-α1 loop, an ‘invariant’ lysine (Lys) residue at the end of the β5 strand and a metal ion (usually Mg2+) that are coordinated through hydrogen bonding [4]. Crystalline and solution structure comparisons of inactive (unphosphorylated) and active (phosphorylated) REC states have revealed a conserved signaling mechanism that centers around an equilibrium shift [5,6] and reorientation of a threonine/serine (Thr/Ser) and tyrosine/phenylalanine (Tyr/Phe) pair toward the phosphorylated Asp residue [7–11]. First described by Zhu et al. [12], this so-called “Y/T-coupling” results in a repositioning of the quaternary structure, usually through the α4-α5-α5 interface of the REC protein architecture, to allow the precise alignment of ionic and hydrophobic residues to form the normally observed active homodimer [13,14], although now a possible alternative dimeric interface centered around the α1-α5 face has also been described [15,16]. After these events the C-terminal Effector Domain (ED) [17], which is often a DNA-binding domain but can also take the form of many other output modules [4,18], is then freed of the physical restraint implemented by close contact with the REC domain to enable downstream signaling [4,15,19]. With a plethora of available genome sequences to sample, a second class of more divergent REC domains has emerged. These sequences typically lack one or more of the aforementioned standard ‘invariant’ features of the canonical variety, including loss of secondary structural features and key residues involved in coordination of the catalytic pocket. Some of these ‘atypical’ REC domains have diversified to accommodate signaling inputs that include binding small molecules such as antibiotics [20], and interestingly Ser-Thr phosphorylation through bacterial eukaryotic-like Ser-Thr kinases and their cognate PP2C phosphatases, which regulate crucial cell envelope-related functions and virulence in many microbes [21,22]. A key feature of many atypical REC domains is the lack of conservation of the phosphorylatable Asp residue, indicating they are regulated by input signals other than histidine phosphorylation. Described examples identified thus far come from bacterial pathogens and environmental microorganisms including RitR from Streptococcus pneumoniae [23,24], AmiR from Pseudomonas aeruginosa [17], RedZ from Streptomyces coelicolor [20], HP1043 from Helicobacter pylori [25], ChxR from Chlamydia trachomatis [26,27], FrzS from Myxococcus xanthus [28], KiaA from Synechococcus elongatus [29], and several have also been described in plants [30], all of which are involved in important cellular processes. As the importance of these non-canonical REC domains are becoming increasingly apparent, we wanted to know how prevalent they are in nature, if they bear functional and structural similarities to canonical REC domains, and importantly how they transduce a signal in the absence of phosphorylation. To gain insight into these questions we first designed a custom program to extract all deposited REC domain sequences that are missing the predicted phosphorylatable Asp residue. Surprisingly, we found that these sequences comprised ~4% (or in bacteria ~2 per completed genome) of all REC domains. Given their representation, we then defined them as the Aspartate-Less Receiver (ALR) domain subclass of atypical RECs. These data reveal that although the largest category of ALR EDs consists of DNA-interacting modules, based on their ED appendages the complete ALR dataset suggests a functional consolidation into more rarely observed specialized roles such as secondary messenger signaling, RNA-binding, Ser phosphorylation and other enzymatic activities. Structural and biochemical analyses of the ALR RitR [23,24], which regulates iron and oxidative stress in the human pathogen Streptococcus pneumoniae, showed that in the absence of typical Asp phospho-regulation RitR likely retains the conserved Y/T-coupling mechanism of activation. Activation is shown to be prevented by an extensive ‘Hydrophobic Gate’ barrier comprised of residues within the canonical α4-β5 dimer interface. Changing the R-group in Gate residues to a methyl group (Ala) enabled RitR to actively dimerize through its α4-β5 face, bind DNA and modulate its target promoter activity (the Pneumococcal Iron Uptake (Piu) transporter) in vivo in the absence of an inducer. Collectively, this work presents the ALR domains and gives insight as to how they might work in the absence of typical phospho-regulatory mechanisms.
To reveal the extent to which the REC domain family was missing the conserved phospho-accepting Asp, we downloaded the available (103,233) REC sequences from the Pfam database, and from this used a series of custom programs to extract a subset of 3,484 sequences lacking the phospho-regulated Asp residue. When redundant sequences were removed 74,816 unique REC-containing proteins were identified, 2,976 of which were missing or had substitutions at the conserved Asp position and indicated that as much as 4% of REC domains do not possess this conserved phosphorylatable residue. Given the frequency of these substitutions we renamed this subset of REC-type sequences Aspartate-Less Receivers (or ALRs). The full dataset of ALR sequences and their accession numbers are given in S5 Table. ALR ED architectures, their accompanying Pfam ED accession numbers and phyla associated with these specific EDs are listed in S6 Table. Examination of our ALR dataset revealed several substitutions within predicted acidic triad positions, which includes the phosphorylatable Asp residue (phospho-Asp). In typical REC sequences these three acidic residues facilitate the phospho-transfer reaction and then help to coordinate the newly formed phosphorylated active pocket. We found that the former phosphorylated Asp position is most frequently substituted with a Glu (26%; Fig 1a). In typical REC sequences when the phospho-Asp is replaced with Glu this can result in a constitutive phosphate-independent activation [31–33]. However, often Asp to Glu substitutions alone are not sufficient to produce such an effect, and therefore must be combined with other mutations to result in constitutive activity (e. g. CheY [34], PhoB [35] and Spo0A [36]). Similar to canonical REC domains, ALR constitutive activity might also be enabled by a change in the phosphorylatable Asp to Glu along with other key functional residues. For example in the case of ChxR, an ALR from Chlamydia trachomatis that contains a Glu in this position, when tested in vitro substitutions with Asp or Ala alone were not sufficient to abrogate its observed constitutive DNA-binding activity [37]. After Glu, the next most frequent substitutions at the conserved phospho-Asp are Asn (16%), Ser (15%), Gly (11%) and Ala (9%) (Fig 1a), all of which have a representative characterized ALR and accompanying atomic structure (except for Gly; Figs 1a and S2a; S3 Table). The remaining 23% of ALR phospho-Asp substitutions are far less frequent, some of which would be predicted to significantly change the typical REC active pocket hydrogen-bonding network. For example, we found 21 ALR sequences with a Pro substitution and another 23 sequences with either a Phe or Tyr aromatic. It is also worth noting that a substantial number of deletions in ALRs are observed in and around the former phospho-Asp position and other pocket residues, which would be predicted to result in changes to the typical REC domain pocket and/or overall structure (Fig 1a). Although only two examples thus far of an ALR display an absent portion of its secondary structure within the pocket region (ChxR [37] and HP1043 [38]), our data presented here suggest that such structural features might be more common in the ALR domain family than in typical RECs. Other more rarely used residues at the former phospho-Asp position include the positively charged (basic) amino acids of Arg, Lys and His, and when combined were found to constitute 4% of ALR phospho-Asp substitutions. One possibility is that such an alteration in charge or a similar change in any of the acidic triad residues might suffice to replace the missing active pocket metal ion (yet to be identified in any ALR). The only amino acid that we did not find in place of the phospho-Asp was Trp, presumably due to the large steric clashes and instability that would likely ensue on the conserved REC three-dimensional architecture. In fact, only one Trp substitution was discovered among all three acid triad residues in ALRs. In general other hydrophobic residues (Val, Ile, Leu, and Met) were also rarely observed in place of any of the acidic triad amino acids (Fig 1a). One exception to this rule appears to be Ala substitutions, where its single methyl side chain might disrupt the overall pocket structure to a lesser extent. In support, the described ALR FrzS structure from M. xanthus takes on a typical REC (α/β) 5 fold where the phospho-Asp has been substituted with an Ala residue (ref. [28] and S2a Fig). Another noteworthy observation was that the N-terminal acidic triad residue (referred to here on as “acidic triad-1”, in RitR it is coordinate Glu9) is generally retained in ALRs (65%), whereas the second acid triad residue (referred to here on as “acidic triad-2”, in RitR it is coordinate Lys10) is only an acidic residue in 29% of ALR sequences (Fig 1a). Furthermore, the conserved Lys that helps coordinate the typical active pocket in REC domains was also largely retained in our ALR dataset (70% of sequences were basic at this position, in RitR it is coordinate Lys103) (Fig 1a), suggesting that this residue and at least one acidic triad residue (more often acidic triad-1) are maintained in ALR sequences. Whether these conservations are present to maintain structural integrity, protein function or both is at present unclear. We also examined the retention of Y/T-coupling residues in ALRs. Based on sequence alignment, available structural information and computational methods, Figs 1a, S1, S2, S3 and S3 Table collectively show that in the majority of cases Tyr/Phe (Y/F) and Thr/Ser (S/T) residues are present in ALRs at the expected positions (66% and 61% of the time, respectively). We also calculated that both a Y/F and S/T residue were present within the same ALR domain 44. 5% of the time, which indicates that in ALRs the pairing of these residues is important, but to a lesser extent than canonical REC domains. Interestingly, if an S/T residue is present in an ALR sequence (61%) then a Y/F is enriched for and appears 77% of the time, whereas if a Y/F is present (66%) then an S/T residue appears at the same rate of 66% (S3 Fig). We also observe that if one residue of the Y/T pair is present, then the other appears with a greater frequency than observed in the total ALR population (for Y/F residues 11% more, and for S/T residues 5% more; S3 Fig). Combined, these observations suggest that in ALRs there is selective pressure to retain Y/T coupling residues, and to a certain extent (44. 5%) to keep them together in the same ALR. Importantly, our data also show that in ALRs it is more important to retain the Y/F residue than the S/T. One explanation for these data could be that to accommodate alternative functions assigned to the sometimes drastically changed ALR “active pockets”, the S/T residue might have to change more frequently as it (but not the Y/T partner) is actually part of the altered catalytic core. Thus the S/T residue could be subjected to a greater evolutionary pressure to accommodate changes in the active pocket specific to a particular ALR function (e. g. atomistic coordination of a bound small molecule). We observe that if one of the Y/T residues is present but not the other, then the substitution tends to be hydrophobic (e. g. Val, Ile, Ala, Leu), and is infrequently charged (RitR is one exception to this rule as it has an Asp residue in place of the Ser/Thr). The reason for these preferences remains to be determined. Taken together our results shown here allude to selective pressure in ALRs to retain Y/T-coupling residues, whose exact roles in signaling will likely have to be determined on an individual experimental basis. Intriguingly, we noticed that many of the substitutions in ALRs are capable of receiving PTMs other than Asp phosphorylation. The most common ALR substitution Glu (Fig 1a) has been shown as phosphorylatable [39], however to the best of our knowledge there are no known examples of such modifications contributing to bacterial signaling. A more likely Glu-specific PTM would be methylation, as this modification has already been convincingly shown to alter bacterial protein signaling function [40]. The next most common substitution, Asn (as seen in RitR; Fig 1a), has the potential to become modified by deamidation to a phosphorylatable Asp. Indeed it has been demonstrated that a REC domain harboring an Asn at the phosphorylatable position can undergo rapid reversion back to Asp [41]. How prevalent this reaction is in Asn-harboring ALRs, and if there are functional consequences of such chemical reversions has yet to be explored with RitR and other ALRs. We also noticed that the third most common ALR substitution at the conserved phospho-Asp site is Ser, and when combined with Thr and Tyr replacements these phosphorylatable residues together constitute 12% (556 sequences) of all identified ALRs (Fig 1a). These amino acids are all well known to hold potential for phospho-regulation by bacterial Ser-Thr [42] and Tyr [43] kinases. Additionally 63 ALR sequences harbor an Arg at the typical REC phospho-Asp position, an amino acid that also has the potential for phospho-regulation in prokaryotes (by arginine-specific kinases and their cognate phosphatases [44,45]). We also noted that the seventh most common substitution at the phospho-Asp position was Cys, a residue capable of sensing environmental changes in the cellular redox state [46]. The 26 identified ALR sequences with a Cys substitution (Cys-ALRs) were further examined by generating an alignment (Fig 2a) and accompanying phylogenetic tree (Fig 2b). Results revealed a distribution among both pathogenic and environmental microbes, many of which require strict oxygen conditions for growth (e. g. Neptuniibacter caesariensis, Bifidobacterium gallicum and Geobacter metallireducens). Cys-ALR EDs were found to cluster into defined outputs and include a DsbA module from the marine aerobe N. caesariensis (Fig 2b), a domain which can act as a disulfide oxidoreductase to ensure proper folding of proteins, especially excreted toxins and virulence factors [47]. Examination of the N. caesariensis genome revealed its encoded Cys-ALR to be adjacent to an operon containing glutathione S-transferase and a predicted NADPH ferrodoxin, enzymes that Dsb family members are often associated with [47] (S7 Fig). Two other predicted genes encoding for Cys-ALRs from the marine bacteria Saccharophagus degradans and human pathogen Stenotrophomonas maltophilia were found adjacent to operons encoding iron/heme uptake systems (S7 Fig). At present PTMs other than histidine kinase phosphorylation have not been described at the conserved REC domain phospho-Asp position. However, there is a plethora of literature describing alternative PTM regulation at other REC and associated ED residues by bacterial kinases [42,43,45,48], as well as by Cys-mediated oxidation [46]. To the best of our knowledge the only example thus far of an atypical REC or ALR domain-containing protein being regulated through PTM modification at any residue is RitR (Ser phosphorylation of the DNA-binding domain [24]). Although speculative, it is exciting to think that some ALRs might be regulated at their former phospho-Asp site by other PTMs (e. g. Ser, Thr, Tyr, Arg phosphorylation or oxidation), whereby the newly-formed modified amino acid would be accommodated by a novel restructured ALR catalytic pocket. Future work will be required to make a definitive verdict on this subject. The phylogenetic numeric distribution of ALRs spanning all three domains of life is depicted in Fig 1b, and also displayed in Fig 1c as the average number of ALRs per Completed Genome (CG) categorized by phyla. Expanded bioinformatic statistics of typical REC domains versus that of ALRs are displayed in S1 Table, and ALR EDs and their association with specific phyla are given in S6 Table. For bacteria we found on average that any given CG possesses 2 ALR sequences, making ALRs a significant portion of the REC-like sequences in prokaryotes. In fact, only 5 out of the 45 known bacterial classes lacked ALRs (the Mollicutes, Thermotogae, Dictyoglomia, Elusimicrobia and Synergistia; S1 Table). Although ALRs have thus far been predominantly characterized in pathogens, our data show they are numerically dominated by environmental bacteria, where greater selective pressures might be responsible for driving REC divergence to accommodate new input signals. We found that the highest average number of ALRs per CG were in the arsenic utilizing Chrystiogenetes (7 per CG), followed by the largely unstudied and mainly soil-dwelling Acidobacteria, the budding Plantomycetacia, the nitrite-oxidizing Nitrosprira and the phylum of Verrucomicrobia (Clamydiae). Photosynthetic bacteria also were enriched for ALRs including the oxygenic photosynthetic Cyanobacteria and the Chloroflexi green non-sulfur bacteria (Fig 1c and S1 Table). In general, the expansive phylum of Proteobacteria was found to contain an average of 2. 2 ALRs per CG. However, the class of δ-proteobacteria, which is largely composed of the related pathogenic species of Helicobacter and Campylobacter, was the most highly enriched in this taxon possessing approximately 5 ALRs per CG (S1 Table). ALRs were found to comprise as much as 10% of the total REC sequences in Archaea, with the classes of Halobacteria and Methanomicrobia containing the highest numbers (Fig 1b and S1 Table). Eukaryotes were also enriched in ALR sequences, which were found in yeast, fungi, ciliated protozoa, land plants and both green and red algae (Fig 1b and S1 Table), but were most abundant in the Eurotiomycetes class of fungi (58 ALRs), the plant phylum of Streptophyta that includes vascular plants (88 ALRs), and finally the highly-evolved Liliopsida (or lilly) class of flowering plants (75 ALRs). A surprising finding was the presence of REC domains and some ALRs in ticks (Arachnida, 2 RECs, 1 ALR), and the primitive invertebrate marine taxa of Anthozoa (mainly sea urchins, 17 RECs, 1 ALR) and Placozoa (4 RECs, 1 ALR), Placozoa being the most basal forms of invertebrates known [49]. As far as we are aware REC sequences and TCS systems have yet to be described in insects and animals. These sequences could be artefacts of bacterial contamination and/or endosymbionts [49], but given the potential for expansion of TCS into the animal kingdom further investigation is warranted. To shed light on the enrichment of ALR ED sequences and how they might correlate with specific taxa, we extracted all ED statistics associated with canonical RECs and ALRs and cross-referenced these results with their taxonomic distribution. Results shown in Fig 1d and S2 Table indicate that the majority (53%) of canonical REC sequences are either ‘stand alone’ REC domains, or possess C-terminal DNA-binding domain extensions, with the largest portion being OmpR/PhoB-type transcription factors such as RitR [50], followed by the LuxR/GerE family. For most DNA-binding effectors, their percentage association with ALRs roughly paralleled the percentage associated with total REC sequences (green versus black bars in Fig 1d). On the other hand, the less commonly seen DNA-binding EDs of YcbB (or GlnL) that possesses a novel helix-turn-helix motif, Myb domains that are dominant in plants, and the β-proteobacterial transcriptional regulator (BetR) family whose function (s) is currently unknown, were all considerably enriched in ALRs (Fig 1d and S2 Table). The remaining 47% of ED sequences are diverse in their functional roles and phylogenetic distribution (a complete listing can be seen in S6 Table). In particular, our data show that ALRs were highly enriched for when attached to RNA binding (ANTAR) motifs, cyclic-di GMP (GGDEF) signaling domains, the Constans-Constans–like TOC1 (CCT) domains involved in controlling plant circadian rhythms, and modules containing both Ser-Thr kinase and phosphatase domains [17,18,30] (Fig 1d and S2 Table). Also highly associated with ALRs were several enzymatic EDs that include pyruvate phosphate dikinases involved in plant C4 carbon metabolism, and HDOD phosphohydrolase domains whose exact functional role (s) is currently unknown [18]. These findings reveal that although ALRs exist in greater overall numbers within the more common ED families (e. g. OmpR-type and LysR), they represent a larger percentage of more rare and less studied REC-associated effector outputs. To further illustrate this point the OmpR-type EDs are associated with the highest raw number of ALRs (609 sequences), yet this figure only constitutes 2% of all REC-[OmpR/PhoB winged helix-turn-helix] architectures. In contrast, out of only 24 REC-CbiA (colbamin binding domain) architectures known to exist, 11 are ALRs (or 46%; Fig 1d; S2 Table). To the best of our knowledge, there are currently five available functionally-characterized prokaryotic ALR structures: (i) AmiR from Pseudomonas aeruginosa [17], (ii) HP1043 from Helicobacter pylori [25], (iii) ChxR from Chlamydia trachomatis [26,27], (iv) FrzS from Myxococcus xanthus [28] and (v) KiaA from Synechococcus elongatus [29]. However, our computational searches revealed an additional 5 unpublished and uncharacterized ALR structures that have been deposited in the RCSB protein databank (www. rcsb. org; PDB IDs 2B4A, 2QZJ, 2ZAY, 3HV2, and 3KTO). For comparison available ALR structures and their known attributes are listed in S3 Table, and an alignment of these sequences and structural comparisons are shown in S2 Fig. To add to the understanding of how ALRs are able to function in the absence of a phospho-Asp signal, and to decipher if RitR parallels these previously solved ALR as well as canonical REC structures, we obtained a crystal structure to 1. 6 Å resolution of the ALR (REC) domain of RitR (RitRALR; PDB ID 4LZL; Table 1), responsible for regulating oxidative stress and iron uptake in the important human pathogen S. pneumoniae [23,24]. As implied by its crystallization in the P1 space group, which lacks crystallographic symmetry elements, RitRALR is monomeric and suggests that this structure could represent an inactive state of the protein [4]. The overall fold of the RitR structure shares the (α/β) 5 topology common to the OmpR/PhoB family of response regulators [4,50,51] (Fig 3a). One notable exception is a defined kink in the α4 helix, which we have subdivided into α4a that interacts with the would-be catalytic pocket in RitR, and α4b that likely interacts with helix 8 of the DBD to preclude its association with DNA [15] (Fig 3b and 3c). A search of available ALR and typical REC atomic structures revealed only two others that exhibit a similar α4 helical bend: PhoP from Mycobacterium tuberculosis [52] (PDB ID 3r0j) and DrrD from Thermotoga maritima [53] (PDB ID 1KGS). Similar to other ALRs and atypical REC sequences, the would-be “active site” of RitR harbors several divergent substitutions relative to classical REC domains (Figs 1a, S1 and S2), all of which gave well-defined electron density for analysis (Fig 4a). First, the phosphorylatable Asp position in RitR is substituted with asparagine (Asn53), thus defining it as a bona fide ALR. Although Asn53 is in approximately the same position as its more canonical Asp-containing counterparts, the conformation differs by a roughly 90° rotation of the χ1 angle where it now hydrogen bonds to acidic triad-1 (Glu9; Figs 3b and S2), which in typical REC domains binds the Mg2+ ion required for the phospho-transfer. When compared with other ALR structures we see that the acid triad-1 position and conserved would-be active pocket Lys (RitR coordinate Lys103) almost without exception make electrostatic interactions to presumably help stabilize the ALR ‘active’ pocket, with the residue that replaces the phospho-Asp site also participating when side chain hydrogen bonding is possible (S2 Fig). Conversely, the acidic triad-2 position side chain only forms a potential hydrogen bond with other ALR acid triad residues in the case of AmiR (S2 Fig). These observations suggest that ALRs have largely retained an acidic residue at the acidic triad-1 position (66%) and a basic residue at the conserved Lys position (70%) to maintain the overall structure of a typical REC domain, while the acid triad-2 remains much more variable (only 29% conserved in ALRs; Fig 1a). Of note, in many ALRs such as RitR, HP1043 and 3KTO acidic triad-2 is replaced by a positively charged (or basic) residue, for example in RitR by Lys10 (Figs 4a and 4b and S2). This normally invariant position of the classical REC domain, located between the β1 strand and α1 helix (Lβ1α1), is crucial in coordinating the active pocket metal and phospho-Asp residue in typical REC sequences. One hypothesis is that such a change to a positive charge might act to take the place of the normally present positively-charged Mg2+ ion, and therefore serve to maintain the structural integrity of the overall REC/ALR fold. In RitR we see that the Lys substitution at acidic triad-2 dramatically alters the position of the entire Lβ1α1 loop relative to typical REC domains such that the cleft, which would normally receive the metal ion, collapses and interrupts what would usually be in typical REC structures a continuous electronegative surface (Fig 4c and 4d). Whatever the reason for these observed changes, our computational and structural data presented here point to ALRs retaining at least one of the two remaining acidic triad residues (usually acidic triad-1), and the conserved Lys-Pro motif, which together help to maintain the classical three-dimensional REC structure in the absence of the usually bound metal ion and typical Asp phosphorylation. Conversely, variation seen within the other two acidic triad residues (the changed phospho-Asp position and acidic triad-2) might be tailored to individual ALR mechanisms. A distinguishing feature of the RitR structure is the α4 helix, which is broken at residues 88–90 into two smaller helices (α4a and α4b; Fig 3a–3c). A careful examination of the available REC structures reveals variations in this region from a straight helical extension. However, few REC α4 structures are as profusely interrupted as the RitR α4, where the measured angle between the helical axes of α4a and α4b is approximately 100°. Situated on the α4 are three hydrophobic residues (Leu86, Leu90, and Val93) that pack against the central β-sheet and make van der Waals contacts with Val54 (Lβ3α3), Ile79 (β4), Ala97 (Lα4β5), Tyr100 (β5) and Tyr102 (β5) to form what we name here the ‘Hydrophobic Gate’ (or ‘Gate’ residues; Fig 3b and 3c). The Gate creates an obvious mechanical barrier to impede the reorientation of Tyr100 (i. e. Y/T-coupling), which in REC domains usually results in α4- β5- α5 dimer association and DNA binding [10,11]. In RitR we noticed that the position of residue Leu90 appeared especially important in that it directly blocks Tyr100 from entering the Gate and also holds the helix 4 kink together through a hydrogen bond (2. 7 Å) to the carbonyl of Glu87. As a result, the side chain of the Tyr100 dimerization (Y/T-coupling) trigger is oriented outward to extend into the solvent—away from the α4-β5- α5 face of the protein as observed in the typical monomeric, inactive REC conformation [7,9]. This positioning might explain why in the case of RitR we observe a predominantly monomeric state in both crystalline and solution environments (see below). The outward-facing Tyr100 rotomer is also enforced by Gate residue Leu86, which in RitR occupies the space normally filled by the Y/T-coupling Thr or Ser residue in the unphosphorylated state. Instead, the equivalent of this residue in RitR (Asp81) orients itself towards Asn53, which is reminiscent of the active or inward-facing rotomer when typical REC domains are phosphorylated. Thus in our RitR structure Y/T-coupling is effectively trapped in an ‘uncoupled’ intermediate state, whereby Asp81 is in the “in” orientation and Tyr100 in the “out” orientation. A survey of the available ALR structures reveals that the RitR Tyr100 outward-facing rotomer is not the norm. In fact, all but one currently available ALR structure exhibits an inward-facing conformation, which could explain why several ALR structures with this rotomer are dimeric (i. e. AmiR (PDB ID 1QO0), HP1043 (PDB ID 2PLN), PDB ID 2QZJ, PDB ID 2ZAY and PDB ID 3HV2; S2 Fig), with some demonstrating DNA binding and activity in vivo without phospho-Asp driven changes influencing their conformational equilibrium states (e. g. HP1043 (PDB: 2PLN) [38] and FrzS (PDB ID 2I6F) [28]). However, in the case of ChxR this trend is contradicted. The available structure reveals that the conserved ChxR Tyr is not only found facing outward, but also participates in dimer formation (S2 Fig; [37]), reinforcing the notion that an observed outward Tyr rotomer in REC/ALR domains does not always result in a monomeric quarternary state [19]. Interestingly, recent data from Kern and colleagues show that the conserved Tyr/Phe and Ser/Thr pair in the response regulator NtrC operate on different timescales in solution, and specifically that the Tyr/Phe does not participate in the active conformation of this protein [54]. Although we did observe that in RitR Tyr100 was important for DNA-binding and in vivo repression of piu (Fig 5), its precise role in RitR activation still remains to be determined. To date it remains enigmatic as to how both atypical REC domains and ALRs are activated in the absence of phosphorylation. Do these proteins exhibit equilibrium between active and inactive states and are driven to activation by an input signal through their Y/T-coupling residues? [5] Or are they instead mainly constitutive and impervious to signaling inputs and resulting downstream structural changes? Recent data has revealed that the activity of several ALRs as well as canonical REC domains can be modulated through ligand or PTM mediated signals other than Asp phosphorylation, which includes Ser/Thr phosphorylation [22,42] and direct binding by small molecule activators [20,55]. The observed conservation of the canonical REC Y/T-coupling residues in ALRs shown here (Figs 1a, S1, and S3) suggests that although the ALR inducer can be variable, after this signal induction takes place in many cases some version of the Y/T-coupling relay and dimerization mechanism may work to modulate ALR activity. At present, it still remains unclear as to why most ALRs have retained these residues and have thus far shown to be functionally important in vivo in the few examples available (ref. [28] and these studies). The unique uncoupled state of the RitR Y/T (Y/D in the case of RitR) pair in our structure presented an opportunity to explore the activation of an ALR in the absence of the inducer. To attempt to force RitR to dimerize in the absence an upstream signal we introduced substitutions within divergent acidic triad residues Lys10 and Asn53 (Fig 5a in cyan), Y/T-coupling residues Tyr100 and Asp81 (Fig 5a in orange), and Gate residues Leu86, Leu90, Val93, Tyr102 and Phe96 (Fig 5a in green). Purified RitR mutant and wild-type proteins (S6 Fig) were then assessed for their oligomeric state via Size Exclusion Chromatography (SEC) and the ability to bind to pre-defined RitR target DNA [23] by Electrophoretic Mobility Shift Assays (EMSAs). Furthermore, to examine the contributions of these individual residues to the RitR mediated repression of piu iron transporter expression in vivo [23,24], we created β-galactosidase reporter strains in S. pneumoniae where we could measure Piu promoter (Ppiu) activity in response to the same RitR isogenic variants (see Methods section for details). Results shown in Fig 5b–5e indicate that mutagenesis of the RitR acidic triad residue coordinates Lys10 and Asn53 had no effect on the ability of RitR to bind DNA, form a dimer in vitro or modulate Ppiu activity in vivo. These data suggest that similar to other functionally-characterized ALRs [25,27,28] and also atypical REC domains that possess the phosphorylatable Asp [9], that at least under these tested conditions RitR does not rely on its acidic triad residue positions to function. Conversely, when Y/T-coupling residues Tyr100 or Asp81 were replaced with Ala we observed DNA binding with as little as 0. 22 μM of RitR in EMSA experiments (Fig 5c and 5d). This was significant in that previous to these experiments DNA binding with full-length RitR (RitRFL) was unobtainable without first removing the ALR (REC) domain (wild-type sample Fig 5c and 5d and ref. [23]), a phenomenon likely due the restraint of the helix-turn-helix motif by intramolecular ALR-DBD domain contacts in the absence of an inducer [19]. SEC experimentation of the Asp81 mutant revealed a split peak profile, neither of which perfectly aligned with the expected monomer or dimer molecular weights, whereas the Tyr100 mutant was mostly monomeric (Fig 5b). However, in vivo results show that Asp81 and Tyr100 are required for full repression of piu transcription (Fig 5e). As REC/ALR dimers are normally the ‘activated’ form of the protein that binds DNA, the SEC results are apparently in conflict with the EMSA and in vivo data. A likely explanation is that in the absence of DNA mutations within Asp81 and Tyr100 destabilize RitR such that it is unable to form a proper dimer in vitro. Experimental support for this hypothesis includes our SEC analysis of Y/T-coupling residues Asp81 and Tyr100 in both RitRFL and RitRALR constructs, as well as our accompanying RitRALR HSQC spectra, where mutations within these residues produces protein with enhanced aggregation and heterogeneic properties (S4 Fig). In contrast to mutagenesis of the RitR Y/T-coupling residues, substitutions within Gate residues Leu86 and Leu90 yielded well-behaved protein that resulted in DNA binding, clear dimer formation in vitro and an enhancement effect on Ppiu activity in vivo (Fig 5). An Ala substitution of Gate residue Val93 alone did not influence RitR DNA binding, but was sufficient to produce more dimer in solution (Fig 5b). However, when the Val93 mutant was combined with a Leu86 mutation, we then observed a large shift in equilibrium towards the dimeric form of RitRFL when compared to the SEC profile of either of these Gate mutations alone (Fig 5b). The Leu86/Val93 double mutant also exhibited DNA binding and modified Ppiu activity in vivo (Fig 5e). Mutation of Gate residue Phe96 produced a marked loss of promoter repression in vivo, however in vitro this mutant displayed a heterogeneous SEC profile and no observed DNA binding activity, suggesting that this residue might be crucial to protein stability. On the other hand a substitution at the Gate residue Tyr102 did not have an effect on RitR DNA binding, modulation of Ppiu activity in vivo or result in a clear dimeric state in solution (Fig 5). From these experiments we can conclude that substitutions within key Gate residues Leu86, Leu90, and Val93 results in the enhancement of dimerization and DNA-binding in vitro, and in vivo these same mutations cause a derepression and/or activation of Ppiu (Fig 5e). However, these in vitro and in vivo data are seemingly in conflict, as RitR is classically known as a repressor of piu transcription [23,24,56], and thus enhanced DNA binding would be expected to result in maintained Ppiu repression in vivo, and not the observed increase seen in Ppiu activity (Fig 5e). In response, recent data from our lab has shown that RitR regulation is much more complex than this simple paradigm, where other PTMs play a pivotal role in its regulatory effects, including Ser-Thr phosphorylation and oxidation ([24] and Wright et al. under review). Moreover, we have found that RitR is also able to act as an activator of piu transcription when occupying different Ppiu transcriptional regulatory sites, which here we have shown can be artificially induced through Y/T-coupling and Gate mutations (Fig 5e). Future studies will have to determine how PTMs, possibly other transcription factors and changing environmental conditions precisely contribute to the complexity of RitR promoter occupancy and transcriptional control in vivo. Although we observed clear RitR dimerization with substitutions in Gate residues, we wanted to know if the RitR dimer interface was facilitated through the typical α4-β5-α5 face used by regular REC domains and other known ALRs ([27,38]; S2 Fig; S3 Table). To answer this question we collected three-dimensional NMR data with [U-15N, 13C] RitR. The NMR spectra displayed excellent peak quality and distribution that would be predicted of a monomeric RitRALR, and were assigned to completion. Next we produced purified 15N-labeled protein of the Ala substitution mutants (ALR domain only, or RitRALR) of Leu86, Leu90, Asp81 and Tyr100. When SEC was used to analyze their oligomeric states as expected the wild-type was predominantly monomeric (Figs 6a and S4a). Conversely, when Gate residue mutants were similarly examined (Leu86 and Leu90) they were almost entirely dimeric, with the Leu90 mutant showing a small monomeric peak, and the Leu86 mutant displaying more heterogeneity in this region with two peaks (Figs 6a and S4a). Paralleling the data of full-length protein SEC experiments, the Y/T-coupling residue mutants were more poorly behaved than Gate substituted versions. SEC data showed the Tyr100 RitRALR mutant to be highly unstable, where most of the protein was found within the void volume (S4a Fig). Moreover the Asp81 mutant formed mainly aggregates after initial Ni2+ affinity purification and was therefore not examined further. These results suggest that the RitR DNA-binding domain aids in the stability of the protein in the active form, which could be further stabilized by the presence of its target DNA sequence as is often seen with helix-turn-helix containing transcription factors. NMR data generally reflected what we observed by SEC, where a new set of peaks likely representing the chemically distinct dimeric state of RitR could be detected in the Leu86 and Leu90 mutants (Figs 6b and S4b). Although not identical, the Leu86 and Leu90 mutant spectra exhibited considerable spectral overlap (S4b Fig), which suggests that these two mutants both form similar dimeric structures. Differences between these spectra could be due the chemical shift contributions within the immediate chemical environment of the mutated residues. When compared to the wild-type HSQC plot, some of the amides from the Leu86/Leu90 mutant samples were found to completely disappear, which we deemed to be in very different chemical environments and contributing to the dimer interface, while other shifts remained in similar chemical environments to varying degrees (Figs 6b and S4b). To better understand what amino acids were changing their immediate chemical environments within the RitR dimer, we superimposed the wild-type HSQC plot over the most homogeneous mutant sample (Leu90; Fig 6b). Residues that experienced the largest differences in the HSQC spectra were then mapped onto the RitRALR crystal structure using a color-coded chemical shift intensity scheme (Fig 6c). The largest changes in chemical shifts (i. e. the peaks in question were either severely diminished in intensity or had completely disappeared in the Leu90 sample) were mostly centered around the α4-β5 face of the protein, which participates in the canonical dimerization region in other ALR/REC structures. In addition to the backbone amide chemical shifts from residues Gln87 and Gln94 being completely absent from the Leu86 and Leu90 mutant HSQC spectra, their respective NH2 side chain chemical shifts were also missing (Fig 6b). Indeed, the Gln87 and Gln94 side chains, both of which project into the solvent in the ‘inactive’ monomeric crystal structure, are predicted to form a component of the α4-β5 dimer interface according to other dimeric ALR and typical REC available atomic models, and in our dimeric RitR atomic model (see below). Collectively these results show that Gate residues Leu90 and Leu86 stabilize the inactive form of the α4, and suggest that Gln87/Gln94 participate in the RitR dimer. Other regions of the protein experienced little change save the N-termini of the α1 and α2 helices (Fig 6c). Our results show that similar to other REC and ALR domains, the α4-β5/Gate region of RitR is likely the major component involved in dimerization of the protein. The observed changes in the α1 and α2 could be explained by another minor and unidentified dimer species in slow exchange with the major α4- β5 dimeric species, or possible intra-protein rearrangements resulting from the larger chemical shift changes in the α4-β5/Gate region. Interestingly, in the Leu86 or Leu90 mutants the α5 helix that would normally complete a canonical dimer face exhibited very little changes in chemical environment, suggesting that this helix does not play a major role in RitR dimer formation. Indeed, the degree of the REC/ALR dimeric interface can vary considerably [15,16,37]. In support of this RitR “α4-β5 only” dimer hypothesis, we found one other ALR structure (PDB ID 3HV2; S5a Fig) and two canonical FixJ regulators ([57,58]; S5a Fig), which all use a α4-β5 dimer interface without α5 participation. Furthermore, when we generated an atomic model of the active RitR dimer based on the structure of the FixJ homolog DctD [57,58] we observed that the amino acids which exhibited the greatest changes in chemical shift/intensities aligned at the α4-β5 dimer interface and were predicted to use the expected crucial Gate residues Leu86, Leu90 and Val93 help form the outward facing inter-protomor contacts (S5b Fig). These contacts also included the two Gln residues (Gln87 and Gln94) whose chemical shifts completely disappeared from their wild-type positions in our Leu86 and Leu90 mutant RitR constructs (Figs 6b, S4b, and S5b). Future three-dimensional structures and further biochemical and genetic analyses will be required to determine the precise contribution of individual residues to dimer formation. Taken together our structural, biochemical and genetic experiments shown here suggest that similar to other REC and ALR domains, RitR uses its conserved Y/T-coupling residues in concert with its Hydrophobic Gate to initiate and form an active dimer (Figs 6 and S5). However, although RitR likely dimerizes at the same face as typical REC domains, we observe that similar to the ALR 3HV2 and FixJ-type transcriptions factors it uses only the α4- β5 to facilitate this interaction. We found that Gate residues Leu86 and Leu90 were especially important in activation of the dimeric form of RitR, and likely involves a repositioning of the conserved Tyr residue from the outward to inward rotomer. Previous studies of typical REC domains have noted the conservation of hydrophobic residues within the α4-β5 Gate, most notably at the RitR Leu90 and Val93 positions, with Leu86 being more variable (see S1 and S2b Figs). In fact in one study it was found that 261 of 269 surveyed OmpR/PhoB-type response regulators possessed a hydrophobic residue at the RitR Leu90 coordinate [53]. In ALRs this residue is also highly conserved as a hydrophobic amino acid (S1 Fig). A survey of the available ALR atomic structures reveals that Leu90 and Val93 are hydrophobic residues in 7/10 and 9/10 of the sequences, respectively (S2b Fig). Moreover, we found that in ALR structures where dimers were available (ChxR (PDB ID 3Q7R), HP1043 (PDB ID 2PLN), AmiA (PDB ID 1QO0), PDB ID 2QZJ and PBD ID 2ZAY), the equivalent coordinates of Leu86, Leu90 and Val93 in all cases participated in their respective dimer interfaces, including in our RitR activated dimer model structure based on DctD (S5b Fig). Combined with our computational data showing their conservation (Figs 1 and S3), these observations further suggest that many ALRs use some form of Y/T-coupling to modulate their activity, and that Gate residues might have a more general function in both canonical REC and ALR domains to ‘fine tune’ the rotomeric states of the conserved Tyr trigger. Additional support comes through studies from Alber and colleagues [28]. They show that similar to most other ALR structures, the ALR FrzS has an inward or “activated” Tyr rotomer. When mutated little difference was observed in protein structure, however in vivo data demonstrated a similar phenotype to a full frzS deletion and suggests that FrzS might use a ‘reverse’ version of Y/T coupling [28]. Considering that most ALR structures solved thus far show their Tyr trigger in the inward active state (S2 Fig), this mechanism could be more widespread within the ALR domain family. At present we can only speculate that ALRs use a version of the Y/T-coupling mechanism of activation, which is for now based on the conservation of Gate and Y/T coupling residues, and data from only two existing examples (i. e. RitR and FtzS). Now that the entire available ALR dataset has been identified (S5 Table), we will be able to more comprehensively study the ALR family of signaling domains to aid in our understanding of their diverse mechanistic roles they play in nature, and importantly how they function in the absence of histidine phosphorylation.
Cultures of pneumococcus were grown overnight from frozen (-80°C) stocks in CAT medium (0. 5% (w/v) tryptone, 1% (w/v) casein digest, 0. 1% (w/v) yeast extract, and 5 mg/L choline) containing 0. 02% (w/v) glucose and the appropriate antibiotics at 37°C in a humidified 5% CO2 incubator. The following day cultures were then diluted 1: 10 in Todd Hewitt broth medium (Becton Dickinson) supplemented with 1% (w/v) yeast extract (THY) containing appropriate antibiotics, and the cell density measured periodically at 600 nm using a Biomate 3 spectrophotometer (Thermo Scientific, Waltham, MA). DNA was transformed using the Competence Stimulating Peptide CSP-1, a generous gift from Donald Morrison, University of Illinois-Chicago. For transformations, S. pneumoniae cells were inoculated from frozen stocks into THY broth containing the appropriate selection antibiotics and cultured to early exponential phase (i. e. an A590 of 0. 01–0. 03). At this time 100 ng of CSP-1 and 100–200 ng of DNA were then added to a 0. 5 ml volume of the culture. The transformation reactions were then placed back into the CO2 chamber and incubated for at least another 2 hours before being plated onto Tryptic Soy Blood Agar (TSBA) plates (Becton Dickinson, Oxford, UK) containing 4% (v/v) defibrinated sheep blood (Rockland Immunochemicals, Gilbertsville, PA) with appropriate antibiotics and grown at 37°C under 5% CO2 until resistant colonies appeared. The kanamycin (Kan) cassette to delete ritR was selected with 200–300 μg/ml Kan, and the erythromycin (Erm) cassette to complement with ritR variants was selected with 1 μg/ml Erm. Strains of E. coli were cultured for plasmid purification overnight with aeration in a 37°C incubator in standard lysogeny broth (LB) medium supplemented with appropriate antibiotics: 50 μg/ml spectinomycin, 34 μg/ml chloramphenicol (Cam), or 50 μg/ml ampicillin (Amp) or the more stable Amp substitute carbenicillin (Carb). Crystals of RitRALR were grown by the hanging-drop, vapor diffusion method. Drops were comprised of equal parts protein solution (5–10 mg/ml RitRALR in 10 mM BIS-TRIS propane, pH 7. 5) and crystallization solution (20–25% (w/v) polyethylene glycol (PEG) 3350,2. 5 mM magnesium formate, and 20 mM TRIS, pH 8. 5). Long, narrow rods appeared after several days. Crystals were prepared for flash-cooling by sequential soaks in solution containing 30% PEG (w/v) 3350,5 mM magnesium formate, 20 mM TRIS (pH 8. 5), and 5,10, or 20% glycerol (v/v). Flash-cooling was accomplished by plunging glycerol-soaked crystals into liquid nitrogen. Diffraction data were collected from a 20 x 20 x 200 μm crystal at beam line 21-ID-D of the Life Science Collaborative Access Team at the Advanced Photon Source (APS), Chicago, Illinois. Data were processed with HKL2000 [67]. Initial phases were obtained by molecular replacement using the program PHASER [68] in the CCP4 suite [69] and the structure of the PhoP REC domain (PDB ID 1MVO, ref. [70]) as the starting model. The initial electron density maps were poor. Automatic rebuilding using the AutoBuild program [71] in the PHENIX suite [72] improved the phases considerably, allowing approximately 95% of the polypeptide to be modeled, albeit with several discontinuities. This improved model was then subjected to iterative cycles of maximum likelihood refinement in the program PHENIX. refine [73] and manual model building in COOT [74]. After completing the polypeptide, riding hydrogen atoms were added to the model with REDUCE [75], and ordered solvent molecules were added with PHENIX. refine. At this point, the crystallographic and free R factors converged at 0. 160 and 0. 195, respectively. Subsequent refinement of anisotropic motion in the form of TLS parameters for groups identified by the TLSMD server [76] resulted in significant reductions in both R factors (Rcryst = 0. 146, Rfree = 0. 167). After final adjustments to the model, including the addition of glycerol molecules from the cryo-protectant solution, the R factors converged at their final values of Rcryst = 0. 136 and Rfree = 0. 157. Model quality was assessed using the comprehensive validation tools implemented in the PHENIX suite, including MolProbity [77]. Data collection and refinement statistics are presented in Table 1. [U-13C, 15N] RitRALR or [U-15N] RitR variants were resuspended in a buffer containing 40 mM 2H-Tris-HCl pH 7. 6,300 mM NaCl, 10% (v/v) D2O, and 0. 02% (v/v) NaN3 and then used to collect 15N-edited, 13C-aliphatic-edited, and 13C-aromatic-edited NOESY spectra on a 600MHz Bruker Advance spectrometer at 25°C. Additional spectra were collected to enable the three-dimensional unambiguous assignments of backbone and side chain resonances of the wild-type sample that included 1H-15N-HSQC, HNCO, HNCA, HN (CO) CA, HN (CO) CACB, HNCACB, C (CO) NH, H (CCO) NH, HBHA (CO) NH, and HCCH-TOCSY experiments. 1H-15N-heteronuclear NOE experiments were also collected to determine the relative order or flexibility of protein on the nanosecond-picosecond timescales (for further discussion of many of these experiments see ref. [78]. Backbone resonances were initially assigned by automation using the program GARANT [79] and were confirmed by manual inspection in XEASY [80]. Side chains were assigned using GARANT’s genPeaks command for the HCCH-TOCSY spectrum, and side chain assignments were edited and corrected manually in XEASY. Aromatic protons were assigned manually in XEASY from the 13C-aromatic-edited NOESY spectrum. SEC experiments were carried out using an Agilent 1220 Compact HPLC equipped with a 250 x 4. 6 mm BioBasic SEC-300A column equilibrated with 50 mM TRIS pH 8. 0,150 mM NaCl and 10 mM DTT. The column was calibrated with the Gel Filtration Molecular Weight Marker kit from Sigma-Aldrich (Cytochrome C (12. 4 kDa), Carbonic Anhydrase (29 kDa), Bovine Serum Albumin (66 kDa), and Sweet Potato Amylase (200 kDa) ). The wild-type and mutant forms of RitRFL protein were pre-treated with 10 mM DTT, injected (5 μl) onto the column and separated at a flow rate of 0. 5 ml/min at ambient temperature. EMSA analysis was conducted using HEX (Hexachlorofluorocein) -labeled double stranded DNA representing the RitR 33-mer Binding Site 2 (BS2) within the Piu promoter [23]. The probes were made by combining at room temperature a forward primer containing the HEX label (BS2-HEX-F) with a complementary and unlabeled reverse primer (BS2–R) (S4 Table). EMSAs were carried out in a final volume of 20 μl containing: 20 mM HEPES pH 7. 2,5 mM MgCl2,1 mM CaCl2,0. 1 mM EDTA, 10 mM DTT, 10% (v/v) glycerol, 0. 5 μM of the double stranded BS2 or control DNA oligos, 800 ng of Poly (deoxyinosinic-deoxycytidylic) acid (poly dI-dC) and RitRFL protein at concentrations between 0 to 6. 6 μM. EMSA reaction mixtures were incubated at room temperature for 10 minutes prior to resolution by 4% non-denaturing 1XTAE (Tris-Acetate EDTA) PAGE. Gels were visualized at 560 nm using a FLA3000 FujiFilm imager. For quantification of RitRFL DNA affinity to BS2 (Fig 4d), a fixed concentration of RitRFL protein (2. 2 μM) was used in a reaction with HEX-labeled BS2 as above. EMSAs were done in triplicate and quantified using ImageJ [81].
Two-Component Signaling (TCS) systems are a primary means by which bacteria sense their constantly changing external environment. For pathogens, these relays are key to their host survival. A typical TCS pair consists of a sensor protein, which upon an environmental stimulation will initiate a cellular response by transferring a phosphate group onto a crucial aspartate amino acid within a secondary receiver (REC) protein. However, more recently it has come to light that some REC proteins are missing key amino acids involved in this signaling relay process, including the aspartate residue required for the signaling relay, hinting that a subset of REC proteins might function in different ways. Here we use custom programing to extract all deposited Aspartate-Less Receivers (ALRs) for further examination. Surprisingly, we found that on average there are two ALRs present in every microbial species, making them a significant portion of the bacterial signaling family. Combining computational data with structural, biochemical and genetic examination of a founding member of the ALR family, Repressor of Iron Transport Regulator (RitR) from the human pathogen Streptococcus pneumoniae, we demonstrate that ALRs might have evolved to accommodate more diverse environmental signals, while largely retaining their time-tested ancestral post-input signaling mechanisms.
The Aspartate-Less Receiver (ALR) Domains: Distribution, Structure and Function