Attractant and repellent induce opposing changes in the four-helix bundle ligand-binding domain of a bacterial chemoreceptor.

Motile bacteria navigate toward favorable conditions and away from unfavorable environments using chemotaxis. Mechanisms of sensing attractants are well understood; however, molecular aspects of how bacteria sense repellents have not been established. Here, we identified malate as a repellent recognized by the MCP2201 chemoreceptor in a bacterium Comamonas testosteroni and showed that it binds to the same site as an attractant citrate. Binding determinants for a repellent and an attractant had only minor differences, and a single amino acid substitution in the binding site inverted the response to malate from a repellent to an attractant. We found that malate and citrate affect the oligomerization state of the ligand-binding domain in opposing way. We also observed opposing effects of repellent and attractant binding on the orientation of an alpha helix connecting the sensory domain to the transmembrane helix. We propose a model to illustrate how positive and negative signals might be generated.

MiST 4.0: a new release of the microbial signal transduction database, now with a metagenomic component.

Signal transduction systems in bacteria and archaea link environmental stimuli to specific adaptive cellular responses. They control gene expression, motility, biofilm formation, development and other processes that are vital to survival. The microbial signal transduction (MiST) database is an online resource that stores tens of thousands of genomes and allows users to explore their signal transduction profiles, analyze genomes in bulk using the database application programming interface (API) and make testable hypotheses about the functions of newly identified signaling systems. However, signal transduction in metagenomes remained completely unexplored. To lay the foundation for research in metagenomic signal transduction, we have prepared a new release of the MiST database, MiST 4.0, which features over 10 000 metagenome-assembled genomes (MAGs), a scaled representation of proteins and detailed BioSample information. In addition, several thousands of new genomes have been processed and stored in the database. A new interface has been developed that allows users to seamlessly switch between genomes and MAGs. MiST 4.0 is freely available at https://mistdb.com; metagenomes and MAGs can also be explored using the API available on the same page.

Amine-recognizing domain in diverse receptors from bacteria and archaea evolved from the universal amino acid sensor.

Bacteria possess various receptors that sense different signals and transmit information to enable an optimal adaptation to the environment. A major limitation in microbiology is the lack of information on the signal molecules that activate receptors. Signals recognized by sensor domains are poorly reflected in overall sequence identity, and therefore, the identification of signals from the amino acid sequence of the sensor alone presents a challenge. Biogenic amines are of great physiological importance for microorganisms and humans. They serve as substrates for aerobic and anaerobic growth and play a role of neurotransmitters and osmoprotectants. Here, we report the identification of a sequence motif that is specific for amine-sensing sensor domains that belong to the Cache superfamily of the most abundant extracellular sensors in prokaryotes. We identified approximately 13,000 sensor histidine kinases, chemoreceptors, receptors involved in second messenger homeostasis and Ser/Thr phosphatases from 8,000 bacterial and archaeal species that contain the amine-recognizing motif. The screening of compound libraries and microcalorimetric titrations of selected sensor domains confirmed their ability to specifically bind biogenic amines. Mutants in the amine-binding motif or domains that contain a single mismatch in the binding motif had either no or a largely reduced affinity for amines. We demonstrate that the amine-recognizing domain originated from the universal amino acid-sensing Cache domain, thus providing insight into receptor evolution. Our approach enables precise "wet"-lab experiments to define the function of regulatory systems and therefore holds a strong promise to enable the identification of signals stimulating numerous receptors.

Evolutionary history of MEK1 illuminates the nature of deleterious mutations.

Mutations in signal transduction pathways lead to various diseases including cancers. MEK1 kinase, encoded by the human MAP2K1 gene, is one of the central components of the MAPK pathway and more than a hundred somatic mutations in the MAP2K1 gene were identified in various tumors. Germline mutations deregulating MEK1 also lead to congenital abnormalities, such as the cardiofaciocutaneous syndrome and arteriovenous malformation. Evaluating variants associated with a disease is a challenge, and computational genomic approaches aid in this process. Establishing evolutionary history of a gene improves computational prediction of disease-causing mutations; however, the evolutionary history of MEK1 is not well understood. Here, by revealing a precise evolutionary history of MEK1, we construct a well-defined dataset of MEK1 metazoan orthologs, which provides sufficient depth to distinguish between conserved and variable amino acid positions. We matched known and predicted disease-causing and benign mutations to evolutionary changes observed in corresponding amino acid positions and found that all known and many suspected disease-causing mutations are evolutionarily intolerable. We selected several variants that cannot be unambiguously assessed by automated prediction tools but that are confidently identified as "damaging" by our approach, for experimental validation in Drosophila. In all cases, evolutionary intolerant variants caused increased mortality and severe defects in fruit fly embryos confirming their damaging nature. We anticipate that our analysis will serve as a blueprint to help evaluate known and novel missense variants in MEK1 and that our approach will contribute to improving automated tools for disease-associated variant interpretation.

Amino acid sensor conserved from bacteria to humans.

SignificanceAmino acids are the building blocks of life and important signaling molecules. Despite their common structure, no universal mechanism for amino acid recognition by cellular receptors is currently known. We discovered a simple motif, which binds amino acids in various receptor proteins from all major life-forms. In humans, this motif is found in subunits of calcium channels that are implicated in pain and neurodevelopmental disorders. Our findings suggest that γ-aminobutyric acid-derived drugs bind to the same motif in human proteins that binds natural ligands in bacterial receptors, thus enabling future improvement of important drugs.

The structural analysis of the periplasmic domain of Sinorhizobium meliloti chemoreceptor McpZ reveals a novel fold and suggests a complex mechanism of transmembrane signaling.

Chemotaxis is a fundamental process whereby bacteria seek out nutrient sources and avoid harmful chemicals. For the symbiotic soil bacterium Sinorhizobium meliloti, the chemotaxis system also plays an essential role in the interaction with its legume host. The chemotactic signaling cascade is initiated through interactions of an attractant or repellent compound with chemoreceptors or methyl-accepting chemotaxis proteins (MCPs). S. meliloti possesses eight chemoreceptors to mediate chemotaxis. Six of these receptors are transmembrane proteins with periplasmic ligand-binding domains (LBDs). The specific functions of McpW and McpZ are still unknown. Here, we report the crystal structure of the periplasmic domain of McpZ (McpZPD) at 2.7 Å resolution. McpZPD assumes a novel fold consisting of three concatenated four-helix bundle modules. Through phylogenetic analyses, we discovered that this helical tri-modular domain fold arose within the Rhizobiaceae family and is still evolving rapidly. The structure, offering a rare view of a ligand-free dimeric MCP-LBD, reveals a novel dimerization interface. Molecular dynamics calculations suggest ligand binding will induce conformational changes that result in large horizontal helix movements within the membrane-proximal domains of the McpZPD dimer that are accompanied by a 5 Å vertical shift of the terminal helix toward the inner cell membrane. These results suggest a mechanism of transmembrane signaling for this family of MCPs that entails both piston-type and scissoring movements. The predicted movements terminate in a conformation that closely mirrors those observed in related ligand-bound MCP-LBDs.

Biallelic CACNA2D1 loss-of-function variants cause early-onset developmental epileptic encephalopathy.

Voltage-gated calcium (CaV) channels form three subfamilies (CaV1-3). The CaV1 and CaV2 channels are heteromeric, consisting of an α1 pore-forming subunit, associated with auxiliary CaVß and α2δ subunits. The α2δ subunits are encoded in mammals by four genes, CACNA2D1-4. They play important roles in trafficking and function of the CaV channel complexes. Here we report biallelic variants in CACNA2D1, encoding the α2δ-1 protein, in two unrelated individuals showing a developmental and epileptic encephalopathy. Patient 1 has a homozygous frameshift variant c.818_821dup/p.(Ser275Asnfs*13) resulting in nonsense-mediated mRNA decay of the CACNA2D1 transcripts, and absence of α2δ-1 protein detected in patient-derived fibroblasts. Patient 2 is compound heterozygous for an early frameshift variant c.13_23dup/p.(Leu9Alafs*5), highly probably representing a null allele and a missense variant c.626G>A/p.(Gly209Asp). Our functional studies show that this amino-acid change severely impairs the function of α2δ-1 as a calcium channel subunit, with strongly reduced trafficking of α2δ-1G209D to the cell surface and a complete inability of α2δ-1G209D to increase the trafficking and function of CaV2 channels. Thus, biallelic loss-of-function variants in CACNA2D1 underlie the severe neurodevelopmental disorder in these two patients. Our results demonstrate the critical importance and non-interchangeability of α2δ-1 and other α2δ proteins for normal human neuronal development.

Photoactive Yellow Protein Represents a Distinct, Evolutionarily Novel Family of PAS Domains.

Photoactive yellow protein (PYP) is a model photoreceptor. It binds a p-coumaric acid as a chromophore, thus enabling blue light sensing. The first discovered single-domain PYP from Halorhodospira halophila has been studied thoroughly in terms of its structural dynamics and photochemical properties. However, the evolutionary origins and biological role of PYP homologs are not well understood. Here, we show that PYP is an evolutionarily novel domain family of the ubiquitous PAS (Per-Arnt-Sim) superfamily. It likely originated from the phylum Myxococcota and was then horizontally transferred to representatives of a few other bacterial phyla. We show that PYP is associated with signal transduction either by domain fusion or by genome context. Key cellular functions modulated by PYP-initiated signal transduction pathways likely involve gene expression, motility, and biofilm formation. We identified three clades of the PYP family, one of which is poorly understood and potentially has novel functional properties. The Tyr42, Glu46, and Cys69 residues that are involved in p-coumaric acid binding in the model PYP from H. halophila are well conserved in the PYP family. However, we also identified cases where substitutions in these residues might have led to neofunctionalization, such as the proposed transition from light to redox sensing. Overall, this study provides definition, a newly built hidden Markov model, and the current genomic landscape of the PYP family and sets the stage for the future exploration of its signaling mechanisms and functional diversity. IMPORTANCE Photoactive yellow protein is a model bacterial photoreceptor. For many years, it was considered a prototypical model of the ubiquitous PAS domain superfamily. Here, we show that, in fact, the PYP family is evolutionarily novel, restricted to a few bacterial phyla and distinct from other PAS domains. We also reveal the diversity of PYP-containing signal transduction proteins and their potential mechanisms.

TREND: a platform for exploring protein function in prokaryotes based on phylogenetic, domain architecture and gene neighborhood analyses.

Key steps in a computational study of protein function involve analysis of (i) relationships between homologous proteins, (ii) protein domain architecture and (iii) gene neighborhoods the corresponding proteins are encoded in. Each of these steps requires a separate computational task and sets of tools. Currently in order to relate protein features and gene neighborhoods information to phylogeny, researchers need to prepare all the necessary data and combine them by hand, which is time-consuming and error-prone. Here, we present a new platform, TREND (tree-based exploration of neighborhoods and domains), which can perform all the necessary steps in automated fashion and put the derived information into phylogenomic context, thus making evolutionary based protein function analysis more efficient. A rich set of adjustable components allows a user to run the computational steps specific to his task. TREND is freely available at http://trend.zhulinlab.org.

MiST 3.0: an updated microbial signal transduction database with an emphasis on chemosensory systems.

Bacteria and archaea employ dedicated signal transduction systems that modulate gene expression, second-messenger turnover, quorum sensing, biofilm formation, motility, host-pathogen and beneficial interactions. The updated MiST database provides a comprehensive classification of microbial signal transduction systems. This update is a result of a substantial scaling to accommodate constantly growing microbial genomic data. More than 125 000 genomes, 516 million genes and almost 100 million unique protein sequences are currently stored in the database. For each bacterial and archaeal genome, MiST 3.0 provides a complete signal transduction profile, thus facilitating theoretical and experimental studies on signal transduction and gene regulation. New software infrastructure and distributed pipeline implemented in MiST 3.0 enable regular genome updates based on the NCBI RefSeq database. A novel MiST feature is the integration of unique profile HMMs to link complex chemosensory systems with corresponding chemoreceptors in bacterial and archaeal genomes. The data can be explored online or via RESTful API (freely available at https://mistdb.com).

A di-iron protein recruited as an Fe[II] and oxygen sensor for bacterial chemotaxis functions by stabilizing an iron-peroxy species.

Many bacteria contain cytoplasmic chemoreceptors that lack sensor domains. Here, we demonstrate that such cytoplasmic receptors found in 8 different bacterial and archaeal phyla genetically couple to metalloproteins related to ß-lactamases and nitric oxide reductases. We show that this oxygen-binding di-iron protein (ODP) acts as a sensor for chemotactic responses to both iron and oxygen in the human pathogen Treponema denticola (Td). The ODP di-iron site binds oxygen at high affinity to reversibly form an unusually stable µ-peroxo adduct. Crystal structures of ODP from Td and the thermophile Thermotoga maritima (Tm) in the Fe[III]2-O22-, Zn[II], and apo states display differences in subunit association, conformation, and metal coordination that indicate potential mechanisms for sensing. In reconstituted systems, iron-peroxo ODP destabilizes the phosphorylated form of the receptor-coupled histidine kinase CheA, thereby providing a biochemical link between oxygen sensing and chemotaxis in diverse prokaryotes, including anaerobes of ancient origin.

The ligand-binding domain of a chemoreceptor from Comamonas testosteroni has a previously unknown homotrimeric structure.

Transmembrane chemoreceptors are widely present in Bacteria and Archaea. They play a critical role in sensing various signals outside and transmitting to the cell interior. Here, we report the structure of the periplasmic ligand-binding domain (LBD) of the transmembrane chemoreceptor MCP2201, which governs chemotaxis to citrate and other organic compounds in Comamonas testosteroni. The apo-form LBD crystal revealed a typical four-helix bundle homodimer, similar to previously well-studied chemoreceptors such as Tar and Tsr of Escherichia coli. However, the citrate-bound LBD revealed a four-helix bundle homotrimer that had not been observed in bacterial chemoreceptor LBDs. This homotrimer was further confirmed with size-exclusion chromatography, analytical ultracentrifugation and cross-linking experiments. The physiological importance of the homotrimer for chemotaxis was demonstrated with site-directed mutations of key amino acid residues in C. testosteroni mutants.

Assigning chemoreceptors to chemosensory pathways in Pseudomonas aeruginosa.

In contrast to Escherichia coli, a model organism for chemotaxis that has 5 chemoreceptors and a single chemosensory pathway, Pseudomonas aeruginosa PAO1 has a much more complex chemosensory network, which consists of 26 chemoreceptors feeding into four chemosensory pathways. While several chemoreceptors were rigorously linked to specific pathways in a series of experimental studies, for most of them this information is not available. Thus, we addressed the problem computationally. Protein-protein interaction network prediction, coexpression data mining, and phylogenetic profiling all produced incomplete and uncertain assignments of chemoreceptors to pathways. However, comparative sequence analysis specifically targeting chemoreceptor regions involved in pathway interactions revealed conserved sequence patterns that enabled us to unambiguously link all 26 chemoreceptors to four pathways. Placing computational evidence in the context of experimental data allowed us to conclude that three chemosensory pathways in P. aeruginosa utilize one chemoreceptor per pathway, whereas the fourth pathway, which is the main system controlling chemotaxis, utilizes the other 23 chemoreceptors. Our results show that while only a very few amino acid positions in receptors, kinases, and adaptors determine their pathway specificity, assigning receptors to pathways computationally is possible. This requires substantial knowledge about interacting partners on a molecular level and focusing comparative sequence analysis on the pathway-specific regions. This general principle should be applicable to resolving many other receptor-pathway interactions.

Comparative study of the effect of disease causing and benign mutations in position Q92 on cholesterol binding by the NPC1 n-terminal domain.

The Niemann-Pick type C1 (NPC1) protein is a large transmembrane protein located in lysosomes/endosomes. NPC1 binds cholesterol (CLR) and transports it to cellular membrane and endoplasmic reticulum. Mutations in NPC1 cause Niemann-Pick type C (NPC) disease, a rare autosomal disorder characterized by intracellular accumulations of CLR and subsequent neurodegeneration leading to premature death. Among known disease-causing mutations in NPC1, Q92R is the one that is located in the N-terminal cholesterol-binding domain [NTD]. Here we study the effect of the mutation on the ability of NPC1 (NTD) to bind and retain CLR in the binding pocket using structural analysis. We compare characteristics of the Q92R and Q92S mutant type (MT) protein, which is predicted to be benign. We provide detailed investigation of the CLR-NPC1 (NTD) binding process; and propose the mechanism, by which Q92R mutation causes NPC disease. We show that although Q92 residue neither directly participates in catalytic activity of the NPC1 (NTD), nor defines its CLR-binding specificity - it is important for the overall protein structure as well as for providing favorable electrostatic environment for CLR transfer. Our results suggest that a negative electrostatic potential of the CLR binding site (the S-opening) might promote NPC2 interaction with NPC1 (NTD) and/or proper CLR orientation and its enforced transfer. We show that in contrast to the benign Q92S mutation, Q92R significantly reduces electrostatic potential around S-opening, and thus likely affects NPC1 (NTD)-NPC2 interaction and/or CLR transfer from NPC2 to NPC1.

Class III Histidine Kinases: a Recently Accessorized Kinase Domain in Putative Modulators of Type IV Pilus-Based Motility.

Histidine kinases are key components of regulatory systems that enable bacteria to respond to environmental changes. Two major classes of histidine kinases are recognized on the basis of their modular design: classical (HKI) and chemotaxis specific (HKII). Recently, a new type of histidine kinase that appeared to have features of both HKIs and HKIIs was identified and termed HKIII; however, the details of HKIII's relationship to other two classes of histidine kinases, their function, and evolutionary history remain unknown. Here, we carried out genomic, phylogenetic, and protein sequence analyses that allowed us to reveal the unusual evolutionary history of this protein family, formalize its distinctive features, and propose its putative function. HKIIIs are characterized by the presence of sensory domains and the lack of a dimerization domain, which is typically present in all histidine kinases. In addition to a single-domain response regulator, HKIII signal transduction systems utilize CheX phosphatase and, in many instances, an unorthodox soluble chemoreceptor that are usual components of chemotaxis signal transduction systems. However, many HKIII genes are found in genomes completely lacking chemotaxis genes, thus decoupling their function from chemotaxis. By contrast, all HKIII-containing genomes also contain pilT, a marker gene for bacterial type IV pilus-based motility, whose regulation is proposed as a putative function for HKIII. These signal transduction systems have a narrow phyletic distribution but are present in many emerging and opportunistic pathogens, thus offering an attractive potential target for future antimicrobial drug design.IMPORTANCE Bacteria adapt to their environment and their hosts by detecting signals and regulating their cellular functions accordingly. Here, we describe a largely unexplored family of signal transduction histidine kinases, called HKIII, that have a unique modular design. While they are currently identified in a relatively short list of bacterial species, this list contains many emerging pathogens. We show that HKIIIs likely control bacterial motility across solid surfaces, which is a key virulence factor in many bacteria, including those causing severe infections. Full understanding of this putative function may help in designing effective drugs against pathogens that will not affect the majority of the beneficial human microbiome.

Aquerium: A web application for comparative exploration of domain-based protein occurrences on the taxonomically clustered genome tree.

Gene duplication and loss are major driving forces in evolution. While many important genomic resources provide information on gene presence, there is a lack of tools giving equal importance to presence and absence information as well as web platforms enabling easy visual comparison of multiple domain-based protein occurrences at once. Here, we present Aquerium, a platform for visualizing genomic presence and absence of biomolecules with a focus on protein domain architectures. The web server offers advanced domain organization querying against the database of pre-computed domains for ∼26,000 organisms and it can be utilized for identification of evolutionary events, such as fusion, disassociation, duplication, and shuffling of protein domains. The tool also allows alternative inputs of custom entries or BLASTP results for visualization. Aquerium will be a useful tool for biologists who perform comparative genomic and evolutionary analyses. The web server is freely accessible at http://aquerium.utk.edu. Proteins 2016; 85:72-77. © 2016 Wiley Periodicals, Inc.

Evolutionary Genomics Suggests That CheV Is an Additional Adaptor for Accommodating Specific Chemoreceptors within the Chemotaxis Signaling Complex.

Escherichia coli and Salmonella enterica are models for many experiments in molecular biology including chemotaxis, and most of the results obtained with one organism have been generalized to another. While most components of the chemotaxis pathway are strongly conserved between the two species, Salmonella genomes contain some chemoreceptors and an additional protein, CheV, that are not found in E. coli. The role of CheV was examined in distantly related species Bacillus subtilis and Helicobacter pylori, but its role in bacterial chemotaxis is still not well understood. We tested a hypothesis that in enterobacteria CheV functions as an additional adaptor linking the CheA kinase to certain types of chemoreceptors that cannot be effectively accommodated by the universal adaptor CheW. Phylogenetic profiling, genomic context and comparative protein sequence analyses suggested that CheV interacts with specific domains of CheA and chemoreceptors from an orthologous group exemplified by the Salmonella McpC protein. Structural consideration of the conservation patterns suggests that CheV and CheW share the same binding spot on the chemoreceptor structure, but have some affinity bias towards chemoreceptors from different orthologous groups. Finally, published experimental results and data newly obtained via comparative genomics support the idea that CheV functions as a "phosphate sink" possibly to off-set the over-stimulation of the kinase by certain types of chemoreceptors. Overall, our results strongly suggest that CheV is an additional adaptor for accommodating specific chemoreceptors within the chemotaxis signaling complex.

Cache Domains That are Homologous to, but Different from PAS Domains Comprise the Largest Superfamily of Extracellular Sensors in Prokaryotes.

Cellular receptors usually contain a designated sensory domain that recognizes the signal. Per/Arnt/Sim (PAS) domains are ubiquitous sensors in thousands of species ranging from bacteria to humans. Although PAS domains were described as intracellular sensors, recent structural studies revealed PAS-like domains in extracytoplasmic regions in several transmembrane receptors. However, these structurally defined extracellular PAS-like domains do not match sequence-derived PAS domain models, and thus their distribution across the genomic landscape remains largely unknown. Here we show that structurally defined extracellular PAS-like domains belong to the Cache superfamily, which is homologous to, but distinct from the PAS superfamily. Our newly built computational models enabled identification of Cache domains in tens of thousands of signal transduction proteins including those from important pathogens and model organisms. Furthermore, we show that Cache domains comprise the dominant mode of extracellular sensing in prokaryotes.