Determinants of homodimerization specificity in histidine kinases

Symmetry facilitated the evolution of heterospecificity and high-order stoichiometry in vertebrate hemoglobin

Many proteins form paralogous multimers-molecular complexes in which evolutionarily related proteins are arranged into specific quaternary structures. Little is known about the mechanisms by which they acquired their stoichiometry (the number of total subunits in the complex) and heterospecificity (the preference of subunits for their paralogs rather than other copies of the same protein). Here, we use ancestral protein reconstruction and biochemical experiments to study historical increases in stoichiometry and specificity during the evolution of vertebrate hemoglobin (Hb), an α2β2 heterotetramer that evolved from a homodimeric ancestor after a gene duplication. We show that the mechanisms for this evolutionary transition were simple. One hydrophobic substitution in subunit β after the gene duplication was sufficient to cause the ancestral dimer to homotetramerize with high affinity across a new interface. During this same interval, a single-residue deletion in subunit α at the older interface conferred specificity for the heterotetrameric form and the trans-orientation of subunits within it. These sudden transitions in stoichiometry and specificity were possible because the interfaces in Hb are isologous, binding via the same surface patch on interacting subunits, but rotated 180° relative to each other. This architecture amplifies the impacts of individual mutations on stoichiometry and specificity, especially in higher-order complexes, and allows single substitutions to differentially affect heteromeric and homomeric interactions. Our findings suggest that elaborate and specific symmetrical molecular complexes may often evolve via simple genetic and physical mechanisms.

Symmetry facilitated the evolution of heterospecificity and high-order stoichiometry in vertebrate hemoglobin

Many proteins form paralogous multimers - molecular complexes in which evolutionarily related proteins are arranged into specific quaternary structures. Little is known about the mechanisms by which they acquired their stoichiometry (the number of total subunits in the complex) and heterospecificity (the preference of subunits for their paralogs rather than other copies of the same protein). Here we use ancestral protein reconstruction and biochemical experiments to study historical increases in stoichiometry and specificity during the evolution of vertebrate hemoglobin (Hb), aα 2 β 2 heterotetramer that evolved from a homodimeric ancestor after a gene duplication. We show that the mechanisms for this evolutionary transition were simple. One hydrophobic substitution in subunit β after the gene duplication was sufficient to cause the ancestral dimer to homotetramerize with high affinity across a new interface. During this same interval, a single-residue deletion in subunit α at the older interface conferred specificity for the heterotetrameric form and the trans-orientation of subunits within it. These sudden transitions in stoichiometry and specificity were possible because the interfaces in Hb are isologous - involving the same surface patch on interacting subunits, rotated 180° relative to each other. This architecture amplifies the impacts of individual mutations on stoichiometry and specificity, especially in higher-order complexes, and allows single substitutions to differentially affect heteromeric vs homomeric interactions. Our findings suggest that elaborate and specific symmetrical molecular complexes may often evolve via simple genetic and physical mechanisms.

Mutational biases favor complexity increases in protein interaction networks after gene duplication

Biological systems can gain complexity over time. While some of these transitions are likely driven by natural selection, the extent to which they occur without providing an adaptive benefit is unknown. At the molecular level, one example is heteromeric complexes replacing homomeric ones following gene duplication. Here, we build a biophysical model and simulate the evolution of homodimers and heterodimers following gene duplication using distributions of mutational effects inferred from available protein structures. We keep the specific activity of each dimer identical, so their concentrations drift neutrally without new functions. We show that for more than 60% of tested dimer structures, the relative concentration of the heteromer increases over time due to mutational biases that favor the heterodimer. However, allowing mutational effects on synthesis rates and differences in the specific activity of homo- and heterodimers can limit or reverse the observed bias toward heterodimers. Our results show that the accumulation of more complex protein quaternary structures is likely under neutral evolution, and that natural selection would be needed to reverse this tendency.

Crosstalk involving two-component systems in Staphylococcus aureus signaling networks

Staphylococcus aureus poses a serious global threat to human health due to its pathogenic nature, adaptation to environmental stress, high virulence, and the prevalence of antimicrobial resistance. The signaling network in S. aureus coordinates and integrates various internal and external inputs and stimuli to adapt and formulate a response to the environment. Two-component systems (TCSs) of S. aureus play a central role in this network where surface-expressed histidine kinases (HKs) receive and relay external signals to their cognate response regulators (RRs). Despite the purported high fidelity of signaling, crosstalk within TCSs, between HK and non-cognate RR, and between TCSs and other systems has been detected widely in bacteria. The examples of crosstalk in S. aureus are very limited, and there needs to be more understanding of its molecular recognition mechanisms, although some crosstalk can be inferred from similar bacterial systems that share structural similarities. Understanding the cellular processes mediated by this crosstalk and how it alters signaling, especially under stress conditions, may help decipher the emergence of antibiotic resistance. This review highlights examples of signaling crosstalk in bacteria in general and S. aureus in particular, as well as the effect of TCS mutations on signaling and crosstalk.

N-terminal domain of ARF-GEF GNOM prevents heterodimerization with functionally divergent GNL1 in Arabidopsis

Evolutionary change following gene duplication can lead to functionally divergent paralogous proteins. If comprising identical subunits their random assortment would also form potentially detrimental heteromeric proteins. In Arabidopsis, the ARF GTPase guanine-nucleotide exchange factor GNOM is essential for polar recycling of auxin-efflux transporter PIN1 from endosomes to the basal plasma membrane whereas its paralog GNL1 mediates retrograde Golgi-endoplasmic reticulum traffic. Here we show that both GNOM and GNL1 form homodimers but no heterodimers. To assess the biological significance of this, we generated transgenic plants expressing engineered heterodimer-compatible GNOM variants. Those plants showed developmental defects such as the failure to produce lateral roots. To identify mechanisms underlying heterodimer prevention, we analyzed interactions of the N-terminal dimerization and cyclophilin-binding (DCB) domain. Each DCB domain interacted with the complementary fragment (ΔDCB) both of their own and of the paralogous protein. However, only DCB^GNOM interacted with itself whereas DCB^GNL1 failed to interact with itself and with DCB^GNOM . GNOM variants in which the DCB domain was removed or replaced by DCB^GNL1 revealed a role for DCB-DCB interaction in the prevention of GNOM-GNL1 heterodimers whereas DCB-ΔDCB interaction was essential for dimer formation and GNOM function. Our data suggest a model of early DCB-DCB interaction that facilitates GNOM homodimer formation, indirectly precluding formation of detrimental heterodimers.

Robust and tunable signal processing in mammalian cells via engineered covalent modification cycles

Engineered signaling networks can impart cells with new functionalities useful for directing differentiation and actuating cellular therapies. For such applications, the engineered networks must be tunable, precisely regulate target gene expression, and be robust to perturbations within the complex context of mammalian cells. Here, we use bacterial two-component signaling proteins to develop synthetic phosphoregulation devices that exhibit these properties in mammalian cells. First, we engineer a synthetic covalent modification cycle based on kinase and phosphatase proteins derived from the bifunctional histidine kinase EnvZ, enabling analog tuning of gene expression via its response regulator OmpR. By regulating phosphatase expression with endogenous miRNAs, we demonstrate cell-type specific signaling responses and a new strategy for accurate cell type classification. Finally, we implement a tunable negative feedback controller via a small molecule-stabilized phosphatase, reducing output expression variance and mitigating the context-dependent effects of off-target regulation and resource competition. Our work lays the foundation for establishing tunable, precise, and robust control over cell behavior with synthetic signaling networks.

Allosteric mechanism of signal transduction in the two-component system histidine kinase PhoQ

Transmembrane signaling proteins couple extracytosolic sensors to cytosolic effectors. Here, we examine how binding of Mg²⁺ to the sensor domain of an E. coli two component histidine kinase (HK), PhoQ, modulates its cytoplasmic kinase domain. We use cysteine-crosslinking and reporter-gene assays to simultaneously and independently probe the signaling state of PhoQ’s sensor and autokinase domains in a set of over 30 mutants. Strikingly, conservative single-site mutations distant from the sensor or catalytic site strongly influence PhoQ’s ligand-sensitivity as well as the magnitude and direction of the signal. Data from 35 mutants are explained by a semi-empirical three-domain model in which the sensor, intervening HAMP, and catalytic domains can adopt kinase-promoting or inhibiting conformations that are in allosteric communication. The catalytic and sensor domains intrinsically favor a constitutively ‘kinase-on’ conformation, while the HAMP domain favors the ‘off’ state; when coupled, they create a bistable system responsive to physiological concentrations of Mg²⁺. Mutations alter signaling by locally modulating domain intrinsic equilibrium constants and interdomain couplings. Our model suggests signals transmit via interdomain allostery rather than propagation of a single concerted conformational change, explaining the diversity of signaling structural transitions observed in individual HK domains.

Manipulation of Bacterial Signaling Using Engineered Histidine Kinases

Two-component systems allow bacteria to respond to changes in environmental or cytosolic conditions through autophosphorylation of a histidine kinase (HK) and subsequent transfer of the phosphate group to its downstream cognate response regulator (RR). The RR then elicits a cellular response, commonly through regulation of transcription. Engineering two-component system signaling networks provides a strategy to study bacterial signaling mechanisms related to bacterial cell survival, symbiosis, and virulence, and to develop sensory devices in synthetic biology. Here we focus on the principles for engineering the HK to identify unknown signal inputs, test signal transmission mechanisms, design small molecule sensors, and rewire two-component signaling networks.

Artificial signaling in mammalian cells enabled by prokaryotic two-component system

Augmenting live cells with novel signal transduction capabilities is a key objective in genetic engineering and synthetic biology. We showed earlier that two-component signaling pathways could function in mammalian cells, albeit while losing their ligand sensitivity. Here we show how to transduce small molecule ligands in a dose-dependent fashion into gene expression in mammalian cells using two-component signaling machinery. First, we engineer mutually complementing truncated mutants of a histidine kinase unable to dimerize and phosphorylate the response regulator. Next, we fuse these mutants to protein domains capable of ligand-induced dimerization, which restores the phosphoryl transfer in a ligand-dependent manner. Cytoplasmic ligands are transduced by facilitating mutant dimerization in the cytoplasm, while extracellular ligands trigger dimerization at the inner side of a plasma membrane. These findings point to the potential of two-component regulatory systems as enabling tools for orthogonal signaling pathways in mammalian cells.

The role of structural pleiotropy and regulatory evolution in the retention of heteromers of paralogs

Gene duplication is a driver of the evolution of new functions. The duplication of genes encoding homomeric proteins leads to the formation of homomers and heteromers of paralogs, creating new complexes after a single duplication event. The loss of these heteromers may be required for the two paralogs to evolve independent functions. Using yeast as a model, we find that heteromerization is frequent among duplicated homomers and correlates with functional similarity between paralogs. Using in silico evolution, we show that for homomers and heteromers sharing binding interfaces, mutations in one paralog can have structural pleiotropic effects on both interactions, resulting in highly correlated responses of the complexes to selection. Therefore, heteromerization could be preserved indirectly due to selection for the maintenance of homomers, thus slowing down functional divergence between paralogs. We suggest that paralogs can overcome the obstacle of structural pleiotropy by regulatory evolution at the transcriptional and post-translational levels.

The role of polyproline motifs in the histidine kinase EnvZ

Although distinct amino acid motifs containing consecutive prolines (polyP) cause ribosome stalling, which necessitates recruitment of the translation elongation factor P (EF-P), they occur strikingly often in bacterial proteomes. For example, polyP motifs are found in more than half of all histidine kinases in Escherichia coli K-12, which raises the question of their role(s) in receptor function. Here we have investigated the roles of two polyP motifs in the osmosensor and histidine kinase EnvZ. We show that the IPPPL motif in the HAMP domain is required for dimerization of EnvZ. Moreover, replacement of the prolines in this motif by alanines disables the receptor’s sensor function. The second motif, VVPPA, which is located in the periplasmic domain, was found to be required for interaction with the modulator protein MzrA. Our study also reveals that polyP-dependent stalling has little effect on EnvZ levels. Hence, both polyP motifs in EnvZ are primarily involved in protein-protein interaction. Furthermore, while the first motif occurs in almost all EnvZ homologues, the second motif is only found in species that have MzrA, indicating co-evolution of the two proteins.

Sodium Lactate Negatively Regulates Shewanella putrefaciens CN32 Biofilm Formation via a Three-Component Regulatory System (LrbS-LrbA-LrbR)

The capability of biofilm formation has a major impact on the industrial and biotechnological applications of Shewanella putrefaciens CN32. However, the detailed regulatory mechanisms underlying biofilm formation in this strain remain largely unknown. In the present report, we describe a three-component regulatory system which negatively regulates the biofilm formation of S. putrefaciens CN32. This system consists of a histidine kinase LrbS (Sputcn32_0303) and two cognate response regulators, including a transcription factor, LrbA (Sputcn32_0304), and a phosphodiesterase, LrbR (Sputcn32_0305). LrbS responds to the signal of the carbon source sodium lactate and subsequently activates LrbA. The activated LrbA then promotes the expression of lrbR, the gene for the other response regulator. The bis-(3'-5')-cyclic dimeric GMP (c-di-GMP) phosphodiesterase LrbR, containing an EAL domain, decreases the concentration of intracellular c-di-GMP, thereby negatively regulating biofilm formation. In summary, the carbon source sodium lactate acts as a signal molecule that regulates biofilm formation via a three-component regulatory system (LrbS-LrbA-LrbR) in S. putrefaciens CN32.IMPORTANCE Biofilm formation is a significant capability used by some bacteria to survive in adverse environments. Numerous environmental factors can affect biofilm formation through different signal transduction pathways. Carbon sources are critical nutrients for bacterial growth, and their concentrations and types significantly influence the biomass and structure of biofilms. However, knowledge about the underlying mechanism of biofilm formation regulation by carbon source is still limited. This work elucidates a modulation pattern of biofilm formation negatively regulated by sodium lactate as a carbon source via a three-component regulatory system in S. putrefaciens CN32, which may serve as a good example for studying how the carbon sources impact biofilm development in other bacteria.

Structure and function of HWE/HisKA2-family sensor histidine kinases

Sensor histidine kinases regulate adaptive cellular responses to changes in the chemical or physical state of the environment. HWE/HisKA2-family kinases comprise a subset of histidine kinases that is defined by unique sequence motifs in both the catalytic and non-catalytic regions. Recent crystal structures have defined conserved intramolecular interactions that inform models of kinase regulation that are unique to the HWE/HisKA2 superfamily. Emerging genetic, biochemical and genomic data indicate that, unlike typical histidine kinases, HWE/HisKA2 kinases do not generally signal via classical DNA-binding response regulators. Rather, these unusual kinases are often part of atypical regulatory pathways that control changes in gene expression via modulation of protein-protein interactions or transcription anti-termination.

Atypical modes of bacterial histidine kinase signaling

The environment of a cell has a profound influence on its physiology, development and evolution. Accordingly, the capacity to sense and respond to physical and chemical signals in the environment is an important feature of cellular biology. In bacteria, environmental sensory perception is often regulated by two-component signal transduction systems (TCSTs). Canonical TCST entails signal-induced autophosphorylation of a sensor histidine kinase (HK) followed by phosphoryl transfer to a cognate response regulator (RR) protein, which may affect gene expression at multiple levels. Recent studies provide evidence for systems that do not adhere to this archetypal TCST signaling model. We present selected examples of atypical modes of signal transduction including inactivation of HK activity via homo- and hetero oligomerization, and cross-phosphorylation between HKs. These examples highlight mechanisms bacteria use to integrate environmental signals to control complex adaptive processes.

Evidence that Autophosphorylation of the Major Sporulation Kinase in Bacillus subtilis Is Able To Occur in trans

Entry into sporulation in Bacillus subtilis is governed by a multicomponent phosphorelay, a complex version of a two-component system which includes at least three histidine kinases (KinA to KinC), two phosphotransferases (Spo0F and Spo0B), and a response regulator (Spo0A). Among the three histidine kinases, KinA is known as the major sporulation kinase; it is autophosphorylated with ATP upon starvation and then transfers a phosphoryl group to the downstream components in a His-Asp-His-Asp signaling pathway. Our recent study demonstrated that KinA forms a homotetramer, not a dimer, mediated by the N-terminal domain, as a functional unit. Furthermore, when the N-terminal domain was overexpressed in the starving wild-type strain, sporulation was impaired. We hypothesized that this impairment of sporulation could be explained by the formation of a nonfunctional heterotetramer of KinA, resulting in the reduced level of phosphorylated Spo0A (Spo0A∼P), and thus, autophosphorylation of KinA could occur in trans. To test this hypothesis, we generated a series of B. subtilis strains expressing homo- or heterogeneous KinA protein complexes consisting of various combinations of the phosphoryl-accepting histidine point mutant protein and the catalytic ATP-binding domain point mutant protein. We found that the ATP-binding-deficient protein was phosphorylated when the phosphorylation-deficient protein was present in a 1:1 stoichiometry in the tetramer complex, while each of the mutant homocomplexes was not phosphorylated. These results suggest that ATP initially binds to one protomer within the tetramer complex and then the γ-phosphoryl group is transmitted to another in a trans fashion. We further found that the sporulation defect of each of the mutant proteins is complemented when the proteins are coexpressed in vivo. Taken together, these in vitro and in vivo results reinforce the evidence that KinA autophosphorylation is able to occur in a trans fashion.

IMPORTANCE

Autophosphorylation of histidine kinases is known to occur by either the cis (one subunit of kinase phosphorylating itself within the multimer) or the trans (one subunit of the multimer phosphorylates the other subunit) mechanism. The present study provided direct in vivo and in vitro evidence that autophosphorylation of the major sporulation histidine kinase (KinA) is able to occur in trans within the homotetramer complex. While the physiological and mechanistic significance of the trans autophosphorylation reaction remains obscure, understanding the detailed reaction mechanism of the sporulation kinase is the first step toward gaining insight into the molecular mechanisms of the initiation of sporulation, which is believed to be triggered by unknown factors produced under conditions of nutrient depletion.

Protein evolution. Pervasive degeneracy and epistasis in a protein-protein interface

Mapping protein sequence space is a difficult problem that necessitates the analysis of 20(N) combinations for sequences of length N. We systematically mapped the sequence space of four key residues in the Escherichia coli protein kinase PhoQ that drive recognition of its substrate PhoP. We generated a library containing all 160,000 variants of PhoQ at these positions and used a two-step selection coupled to next-generation sequencing to identify 1659 functional variants. Our results reveal extensive degeneracy in the PhoQ-PhoP interface and epistasis, with the effect of individual substitutions often highly dependent on context. Together, epistasis and the genetic code create a pattern of connectivity of functional variants in sequence space that likely constrains PhoQ evolution. Consequently, the diversity of PhoQ orthologs is substantially lower than that of functional PhoQ variants.

Data-driven prediction and design of bZIP coiled-coil interactions

Selective dimerization of the basic-region leucine-zipper (bZIP) transcription factors presents a vivid example of how a high degree of interaction specificity can be achieved within a family of structurally similar proteins. The coiled-coil motif that mediates homo- or hetero-dimerization of the bZIP proteins has been intensively studied, and a variety of methods have been proposed to predict these interactions from sequence data. In this work, we used a large quantitative set of 4,549 bZIP coiled-coil interactions to develop a predictive model that exploits knowledge of structurally conserved residue-residue interactions in the coiled-coil motif. Our model, which expresses interaction energies as a sum of interpretable residue-pair and triplet terms, achieves a correlation with experimental binding free energies of R = 0.68 and significantly out-performs other scoring functions. To use our model in protein design applications, we devised a strategy in which synthetic peptides are built by assembling 7-residue native-protein heptad modules into new combinations. An integer linear program was used to find the optimal combination of heptads to bind selectively to a target human bZIP coiled coil, but not to target paralogs. Using this approach, we designed peptides to interact with the bZIP domains from human JUN, XBP1, ATF4 and ATF5. Testing more than 132 candidate protein complexes using a fluorescence resonance energy transfer assay confirmed the formation of tight and selective heterodimers between the designed peptides and their targets. This approach can be used to make inhibitors of native proteins, or to develop novel peptides for applications in synthetic biology or nanotechnology.

Conservation weighting functions enable covariance analyses to detect functionally important amino acids

The explosive growth in the number of protein sequences gives rise to the possibility of using the natural variation in sequences of homologous proteins to find residues that control different protein phenotypes. Because in many cases different phenotypes are each controlled by a group of residues, the mutations that separate one version of a phenotype from another will be correlated. Here we incorporate biological knowledge about protein phenotypes and their variability in the sequence alignment of interest into algorithms that detect correlated mutations, improving their ability to detect the residues that control those phenotypes. We demonstrate the power of this approach using simulations and recent experimental data. Applying these principles to the protein families encoded by Dscam and Protocadherin allows us to make testable predictions about the residues that dictate the specificity of molecular interactions.

Increasing the affinity of selective bZIP-binding peptides through surface residue redesign

The coiled-coil dimer is a prevalent protein interaction motif that is important for many cellular processes. The basic leucine-zipper (bZIP) transcription factors are one family of proteins for which coiled-coil mediated dimerization is essential for function, and misregulation of bZIPs can lead to disease states including cancer. This makes coiled coils attractive protein-protein interaction targets to disrupt using engineered molecules. Previous work designing peptides to compete with native coiled-coil interactions focused primarily on designing the core residues of the interface to achieve affinity and specificity. However, folding studies on the model bZIP GCN4 show that coiled-coil surface residues also contribute to binding affinity. Here we extend a prior study in which peptides were designed to bind tightly and specifically to representative members of each of 20 human bZIP families. These "anti-bZIP" peptides were designed with an emphasis on target-binding specificity, with contributions to design-target specificity and affinity engineered considering only the coiled-coil core residues. High-throughput testing using peptide arrays indicated many successes. We have now measured the binding affinities and specificities of anti-bZIPs that bind to FOS, XBP1, ATF6, and CREBZF in solution and tested whether redesigning the surface residues can increase design-target affinity. Incorporating residues that favor helix formation into the designs increased binding affinities in all cases, providing low-nanomolar binders of each target. However, changes in surface electrostatic interactions sometimes changed the binding specificity of the designed peptides.

Vestigialization of an allosteric switch: genetic and structural mechanisms for the evolution of constitutive activity in a steroid hormone receptor

An important goal in molecular evolution is to understand the genetic and physical mechanisms by which protein functions evolve and, in turn, to characterize how a protein's physical architecture influences its evolution. Here we dissect the mechanisms for an evolutionary shift in function in the mollusk ortholog of the steroid hormone receptors (SRs), a family of biologically essential transcription factors. In vertebrates, the activity of SRs allosterically depends on binding a hormonal ligand; in mollusks, however, the SR ortholog (called ER, because of high sequence similarity to vertebrate estrogen receptors) activates transcription in the absence of ligand and does not respond to steroid hormones. To understand how this shift in regulation evolved, we combined evolutionary, structural, and functional analyses. We first determined the X-ray crystal structure of the ER of the Pacific oyster Crassostrea gigas (CgER), and found that its ligand pocket is filled with bulky residues that prevent ligand occupancy. To understand the genetic basis for the evolution of mollusk ERs' unique functions, we resurrected an ancient SR progenitor and characterized the effect of historical amino acid replacements on its functions. We found that reintroducing just two ancient replacements from the lineage leading to mollusk ERs recapitulates the evolution of full constitutive activity and the loss of ligand activation. These substitutions stabilize interactions among key helices, causing the allosteric switch to become “stuck” in the active conformation and making activation independent of ligand binding. Subsequent changes filled the ligand pocket without further affecting activity; by degrading the allosteric switch, these substitutions vestigialized elements of the protein's architecture required for ligand regulation and made reversal to the ancestral function more complex. These findings show how the physical architecture of allostery enabled a few large-effect mutations to trigger a profound evolutionary change in the protein's function and shaped the genetics of evolutionary reversibility.

An important goal in evolutionary genetics is to understand how genetic mutations cause the evolution of new protein functions and how a protein's structure shapes its evolution. Here we address these questions by studying a dramatic lineage-specific shift in function in steroid hormone receptors (SRs), a physiologically important family of transcription factors. In vertebrates, SRs bind hormones and then undergo a structural change that allows them to activate gene expression. In mollusks, SRs do not bind hormone and are always active. We identified the genetic and structural mechanisms for the evolution of constitutive activity in the mollusk SRs by using X-ray crystallography, ancestral sequence reconstruction, and experimental studies of the effects of ancient mutations on protein structure and function. We found that constitutive activity evolved due to just two historical substitutions that subtly stabilized elements of the active conformation, and subsequent mutations filled the hormone-binding cavity. The structural characteristics required for a hormone-sensitive activator were thus vestigialized, much the same way that a whale's hindlimbs became vestiges of their ancestral form after they became dispensable. Our findings show how the architecture of a protein can shape its evolution, allowing radically different functions to evolve by a few large-effect mutations.

The MprB extracytoplasmic domain negatively regulates activation of the Mycobacterium tuberculosis MprAB two-component system

Mycobacterium tuberculosis is an acid-fast pathogen of humans and the etiological agent of tuberculosis (TB). It is estimated that one-third of the world's population is latently (persistently) infected with M. tuberculosis. M. tuberculosis persistence is regulated, in part, by the MprAB two-component signal transduction system, which is activated by and mediates resistance to cell envelope stress. Here we identify MprAB as part of an evolutionarily conserved cell envelope stress response network and demonstrate that MprAB-mediated signal transduction is negatively regulated by the MprB extracytoplasmic domain (ECD). In particular, we report that deregulated production of the MprB sensor kinase, or of derivatives of this protein, negatively impacts M. tuberculosis growth. The observed growth attenuation is dependent on MprAB-mediated signal transduction and is exacerbated in strains of M. tuberculosis producing an MprB variant lacking its ECD. Interestingly, full-length MprB, and the ECD of MprB specifically, immunoprecipitates the Hsp70 chaperone DnaK in vivo, while overexpression of dnaK inhibits MprAB-mediated signal transduction in M. tuberculosis grown in the absence or presence of cell envelope stress. We propose that under nonstress conditions, or under conditions in which proteins present in the extracytoplasmic space are properly folded, signaling through the MprAB system is inhibited by the MprB ECD. Following exposure to cell envelope stress, proteins present in the extracytoplasmic space become unfolded or misfolded, leading to removal of the ECD-mediated negative regulation of MprB and subsequent activation of MprAB.

Networks of bZIP protein-protein interactions diversified over a billion years of evolution

Differences in biomolecular sequence and function underlie dramatic ranges of appearance and behavior among species. We studied the basic region-leucine zipper (bZIP) transcription factors and quantified bZIP dimerization networks for five metazoan and two single-cell species, measuring interactions in vitro for 2891 protein pairs. Metazoans have a higher proportion of heteromeric bZIP interactions and more network complexity than the single-cell species. The metazoan bZIP interactomes have broadly similar structures, but there has been extensive rewiring of connections compared to the last common ancestor, and each species network is highly distinct. Many metazoan bZIP orthologs and paralogs have strikingly different interaction specificities, and some differences arise from minor sequence changes. Our data show that a shifting landscape of biochemical functions related to signaling and gene expression contributes to species diversity.

Missense substitutions reflecting regulatory control of transmitter phosphatase activity in two-component signalling

Negative control in two-component signal transduction results from sensor transmitter phosphatase activity for phospho-receiver dephosphorylation. A hypothetical mechanism for this reaction involves a catalytic residue in the H-box active-site region. However, a complete understanding of transmitter phosphatase regulation is hampered by the abundance of kinase-competent, phosphatase-defective missense substitutions (K(+) P(-) phenotype) outside of the active-site region. For the Escherichia coli NarX sensor, a model for the HisKA_3 sequence family, DHp domain K(+) P(-) mutants defined two classes. Interaction mutants mapped to the active site-distal base of the DHp helix 1, whereas conformation mutants were affected in the X-box region of helix 2. Thus, different types of perturbations can influence transmitter phosphatase activity indirectly. By comparison, K(+) P(-) substitutions in the HisKA sensors EnvZ and NtrB additionally map to a third region, at the active site-proximal top of the DHp helix 1, independently identified as important for DHp-CA domain interaction in this sensor class. Moreover, the NarX transmitter phosphatase activity was independent of nucleotides, in contrast to the activity for many HisKA family sensors. Therefore, distinctions involving both the DHp and the CA domains suggest functional diversity in the regulation of HisKA and HisKA_3 transmitter phosphatase activities.

Determinants of specificity in two-component signal transduction

Maintaining the faithful flow of information through signal transduction pathways is critical to the survival and proliferation of organisms. This problem is particularly challenging as many signaling proteins are part of large, paralogous families that are highly similar at the sequence and structural levels, increasing the risk of unwanted cross-talk. To detect environmental signals and process information, bacteria rely heavily on two-component signaling systems comprised of sensor histidine kinases and their cognate response regulators. Although most species encode dozens of these signaling pathways, there is relatively little cross-talk, indicating that individual pathways are well insulated and highly specific. Here, we review the molecular mechanisms that enforce this specificity. Further, we highlight recent studies that have revealed how these mechanisms evolve to accommodate the introduction of new pathways by gene duplication.

Helix bundle loops determine whether histidine kinases autophosphorylate in cis or in trans

Bacteria frequently use two-component signal transduction pathways to sense and respond to environmental and intracellular stimuli. Upon receipt of a stimulus, a homodimeric sensor histidine kinase autophosphorylates and then transfers its phosphoryl group to a cognate response regulator. The autophosphorylation of histidine kinases has been reported to occur both in cis and in trans, but the molecular determinants dictating which mechanism is employed are unknown. Based on structural considerations, one model posits that the handedness of a loop at the base of the helical dimerization domain plays a critical role. Here, we tested this model by replacing the loop from Escherichia coli EnvZ, which autophosphorylates in trans, with the loop from three PhoR orthologs that autophosphorylate in cis. These chimeric kinases autophosphorylated in cis, indicating that this small loop is sufficient to determine autophosphorylation mechanism. Further, we report that the mechanism of autophosphorylation is conserved in orthologous sets of histidine kinases despite highly dissimilar loop sequences. These findings suggest that histidine kinases are under selective pressure to maintain their mode of autophosphorylation, but they can do so with a wide range of sequences.

Using analyses of amino Acid coevolution to understand protein structure and function

Determining which residues of a protein contribute to a specific function is a difficult problem. Analyses of amino acid covariation within a protein family can serve as a useful guide by identifying residues that are functionally coupled. Covariation analyses have been successfully used on several different protein families to identify residues that work together to promote folding, enable protein-protein interactions, or contribute to an enzymatic activity. Covariation is a statistical signal that can be measured in a multiple sequence alignment of homologous proteins. As sequence databases have expanded dramatically, covariation analyses have become easier and more powerful. In this chapter, we describe how functional covariation arises during the evolution of proteins and how this signal can be distinguished from various background signals. We discuss the basic methodology for performing amino acid covariation analysis, using bacterial two-component signal transduction proteins as an example. We provide practical suggestions for each step of the process including assembly of protein sequences, construction of a multiple sequence alignment, measurement of covariation, and analysis of results.

Statistical analyses of protein sequence alignments identify structures and mechanisms in signal activation of sensor histidine kinases

Statistical analyses of genome sequence-derived protein sequence data can identify amino acid residues that interact between proteins or between domains of a protein. These statistical methods are based on evolution-directed amino acid variation responding to structural and functional constraints in proteins. The identified residues form a basis for determining structure and folding of proteins as well as inferring mechanisms of protein function. When applied to two-component systems, several research groups have shown they can be used to identify the amino acid interactions between response regulators and histidine kinases and the specificity therein. Recently, statistical studies between the HisKA and HATPase-ATP-binding domains of histidine kinases identified amino acid interactions for both the inactive and the active catalytic states of such kinases. The identified interactions generated a model structure for the domain conformation of the active state. This conformation requires an unwinding of a portion of the C-terminal helix of the HisKA domain that destroys the inactive state residue contacts and suggests how signal-binding determines the equilibrium between the inactive and active states of histidine kinases. The rapidly accumulating protein sequence databases from genome, metagenome and microbiome studies are an important resource for functional and structural understanding of proteins and protein complexes in microbes.

Intramolecular arrangement of sensor and regulator overcomes relaxed specificity in hybrid two-component systems

Cellular processes require specific interactions between cognate protein partners and concomitant discrimination against noncognate partners. Signal transduction by classical two-component regulatory systems typically entails an intermolecular phosphoryl transfer between a sensor kinase (SK) and a cognate response regulator (RR). Interactions between noncognate partners are rare because SK/RR pairs coevolve unique interfaces that dictate phosphotransfer specificity. Here we report that the in vitro phosphotransfer specificity is relaxed in hybrid two-component systems (HTCSs) from the human gut symbiont Bacteroides thetaiotaomicron, which harbor both the SK and RR in a single polypeptide. In contrast, phosphotransfer specificity is retained in classical two-component regulatory systems from this organism. This relaxed specificity enabled us to rewire a HTCS successfully to transduce signals between noncognate SK/RR pairs. Despite the relaxed specificity between SK and RRs, HTCSs remained insulated from cross-talk with noncognate proteins in vivo. Our data suggest that the high local concentration of the SK and RR present in the same polypeptide maintains specificity while relaxing the constraints on coevolving unique contact interfaces.

bZIP transcription factors affecting secondary metabolism, sexual development and stress responses in Aspergillus nidulans

The eukaryotic basic leucine zipper (bZIP) transcription factors play critical roles in the organismal response to the environment. Recently, a novel YAP-like bZIP, restorer of secondary metabolism A (RsmA), was found in a suppressor screen of an Aspergillus nidulans secondary metabolism (SM) mutant in which overexpression of rsmA was found to partially remediate loss of SM in Velvet Complex mutants. The Velvet Complex is a conserved fungal transcriptional heteromer that couples SM with sexual development in fungi. Here we characterized and contrasted SM in mutants of RsmA and four other A. nidulans bZIP proteins (NapA, ZipA, ZipB and ZipC) with predicted DNA binding motifs similar to RsmA. Only two overexpression mutants exhibited both SM and sexual abnormalities that were noteworthy: OE : : rsmA resulted in a 100-fold increase in sterigmatocystin and a near loss of meiotic spore production. OE : : napA displayed decreased production of sterigmatocystin, emericellin, asperthecin, shamixanthone and epishamixanthone, coupled with a shift from sexual to asexual development. Quantification of bZIP homodimer and heterodimer formation using fluorescence resonance energy transfer (FRET) suggested that these proteins preferentially self-associate.

Structural insights into ChpT, an essential dimeric histidine phosphotransferase regulating the cell cycle in Caulobacter crescentus

Two-component and phosphorelay signal-transduction proteins are crucial for bacterial cell-cycle regulation in Caulobacter crescentus. ChpT is an essential histidine phosphotransferase that controls the activity of the master cell-cycle regulator CtrA by phosphorylation. Here, the 2.2 Å resolution crystal structure of ChpT is reported. ChpT is a homodimer and adopts the domain architecture of the intracellular part of class I histidine kinases. Each subunit consists of two distinct domains: an N-terminal helical hairpin domain and a C-terminal α/β domain. The two N-terminal domains are adjacent within the dimer, forming a four-helix bundle. The ChpT C-terminal domain adopts an atypical Bergerat ATP-binding fold.

Evolution of two-component signal transduction systems

To exist in a wide range of environmental niches, bacteria must sense and respond to a variety of external signals. A primary means by which this occurs is through two-component signal transduction pathways, typically composed of a sensor histidine kinase that receives the input stimuli and then phosphorylates a response regulator that effects an appropriate change in cellular physiology. Histidine kinases and response regulators have an intrinsic modularity that separates signal input, phosphotransfer, and output response; this modularity has allowed bacteria to dramatically expand and diversify their signaling capabilities. Recent work has begun to reveal the molecular basis by which two-component proteins evolve. How and why do orthologous signaling proteins diverge? How do cells gain new pathways and recognize new signals? What changes are needed to insulate a new pathway from existing pathways? What constraints are there on gene duplication and lateral gene transfer? Here, we review progress made in answering these questions, highlighting how the integration of genome sequence data with experimental studies is providing major new insights.

Designing specific protein-protein interactions using computation, experimental library screening, or integrated methods

Given the importance of protein-protein interactions for nearly all biological processes, the design of protein affinity reagents for use in research, diagnosis or therapy is an important endeavor. Engineered proteins would ideally have high specificities for their intended targets, but achieving interaction specificity by design can be challenging. There are two major approaches to protein design or redesign. Most commonly, proteins and peptides are engineered using experimental library screening and/or in vitro evolution. An alternative approach involves using protein structure and computational modeling to rationally choose sequences predicted to have desirable properties. Computational design has successfully produced novel proteins with enhanced stability, desired interactions and enzymatic function. Here we review the strengths and limitations of experimental library screening and computational structure-based design, giving examples where these methods have been applied to designing protein interaction specificity. We highlight recent studies that demonstrate strategies for combining computational modeling with library screening. The computational methods provide focused libraries predicted to be enriched in sequences with the properties of interest. Such integrated approaches represent a promising way to increase the efficiency of protein design and to engineer complex functionality such as interaction specificity.

CitA/CitB two-component system regulating citrate fermentation in Escherichia coli and its relation to the DcuS/DcuR system in vivo

Citrate fermentation by Escherichia coli requires the function of the citrate/succinate antiporter CitT (citT gene) and of citrate lyase (citCDEFXG genes). Earlier experiments suggested that the two-component system CitA/CitB, consisting of the membrane-bound sensor kinase CitA and the response regulator CitB, stimulates the expression of the genes in the presence of citrate, similarly to CitA/CitB of Klebsiella pneumoniae. In this study, the expression of a chromosomal citC-lacZ gene fusion was shown to depend on CitA/CitB and citrate. CitA/CitB is related to the DcuS/DcuR two-component system which induces the expression of genes for fumarate respiration in response to C(4)-dicarboxylates and citrate. Unlike DcuS, CitA required none of the cognate transporters (CitT, DcuB, or DcuC) for function, and the deletion of the corresponding genes showed no effect on the expression of citC-lacZ. The citAB operon is preceded by a DcuR binding site. Phosphorylated DcuR bound specifically to the promoter region, and the deletion of dcuS or dcuR reduced the expression of citC. The data indicate the presence of a regulatory cascade consisting of DcuS/DcuR modulating citAB expression (and CitA/CitB levels) and CitA/CitB controlling the expression of the citCDEFXGT gene cluster in response to citrate. In vivo fluorescence resonance energy transfer (FRET) and the bacterial two-hybrid system (BACTH) showed interaction between the DcuS and CitA proteins. However, BACTH and expression studies demonstrated the lack of interaction and cross-regulation between CitA and DcuR or DcuS and CitB. Therefore, there is only linear phosphoryl transfer (DcuS→DcuR and CitA→CitB) without cross-regulation between DcuS/DcuR and CitA/CitB.