Invariant Amino Acid Research Articles

Ancestral role of Pax6 in chordate brain regionalization.

The Pax6 gene is essential for eye and brain development across various animal species. Here, we investigate the function of Pax6 in the development of the anterior central nervous system (CNS) of the invertebrate chordate amphioxus using CRISPR/Cas9-induced genome editing. Specifically, we examined Pax6 mutants featuring a 6bp deletion encompassing two invariant amino acids in the conserved paired domain, hypothesized to impair Pax6 DNA-binding capacity and gene regulatory functions. Although this mutation did not result in gross morphological changes in amphioxus larvae, it demonstrated a reduced ability to activate Pax6-responsive reporter gene, suggesting a hypomorphic effect. Expression analysis in mutant larvae revealed changes in gene expression within the anterior CNS, supporting the conserved role of Pax6 gene in brain regionalization across chordates. Additionally, our findings lend support to the hypothesis of a zona limitans intrathalamica (ZLI)-like region in amphioxus, suggesting evolutionary continuity in brain patterning mechanisms. ZLI region, found in both hemichordates and vertebrates, functions as a key signaling center and serves as a restrictive boundary between major thalamic regions.

Selective Peptide Binders to the Perfluorinated Sulfonic Acid Ionomer Nafion

AbstractFuel cells used for transport applications hold polymer membranes as a key element. Their efficiency can be significantly increased if structured ion channels are implemented at the molecular level into the proton‐conducting membrane. New functional molecules with selective affinity for ionomers are needed to obtain such a membrane design. This study presents a method to screen for selective peptide binders to perfluorinated sulfonic acid ionomers, e.g., Nafion using ultra‐high density peptide arrays with a spot size of up to 30 µm. First, the ionomer dispersion is incubated with the peptide chip containing 56014 randomly chosen 6‐mer peptides. Afterward, the peptide chip is washed. The peptide WIWHCW with the helix structure is identified as a selective binder to Nafion. The invariant amino acids responsible for binding are also determined using a peptide chip approach. The specific binding pocket of WIWHCW is formed by histidine and tryptophan. Its dissociation constant to the ionomer is measured at ≈140 µM.

An insight into the active sites of the catalytic basic protein subunit PB1 of the RNA polymerase of human influenza viruses

RNA polymerase from human influenza viruses is a heterotrimeric enzyme and performs the crucial functions of both transcription and replication for multiplication of the viruses in human cells. The heterotrimeric enzyme is made up of two basic protein subunits (PB1 and PB2) and an acidic protein subunit (PA). All the three subunits perform well-defined function(s) in the transcription and replication processes in the human cells. The basic protein subunit PB1 is shown to possess the polymerization activity, whereas the PB2 and the PA subunits are found to be involved in the cap-snatching and proofreading (PR) activities, respectively. The polymerase activity in the catalytic subunit, PB1, is found to be an RNA-dependent RNA polymerase (RdRp). Multiple sequence alignment (MSA) analysis of the PB1 subunits from all the three human influenza viruses, A, B and C shows large number of highly conserved peptides, amino acid motifs and invariant amino acids. Site-directed mutagenesis (SDM) analysis and X-ray crystallographic data have shown that two completely conserved motifs, viz. –GDN- and –SDD-, are involved in binding to the catalytic Mg2+ ion. These data are in close agreement with the MSA analysis data of the polymerases from all the three human influenza viruses. Furthermore, two highly conserved polymerase catalytic regions are identified in the PB1 subunits by sequence similarity to other DNA/RNA polymerases and hence, are proposed to function in the nucleotidyl transfer activities. Presence of the two catalytic regions suggest that the polymerase may function in a dual mode, i.e., in phase I, in association with the cap-snatching subunit PB2, it could be involved in the synthesis of mRNAs (transcription mode) and once enough proteins are made from the mRNAs, in the second phase, in association with PR exonuclease subunit PA, it could switch to the replication mode to synthesize error-free, exact copies of the viral genome. For both the activities, it could use the same invariant catalytic Mg2+-binding –GDN- and –SDD- motifs.

Resemblance-Ranking Peptide Library to Screen for Binders to Antibodies on a Peptidomic Scale.

A novel resemblance-ranking peptide library with 160,000 10-meric peptides was designed to search for selective binders to antibodies. The resemblance-ranking principle enabled the selection of sequences that are most similar to the human peptidome. The library was synthesized with ultra-high-density peptide arrays. As proof of principle, screens for selective binders were performed for the therapeutic anti-CD20 antibody rituximab. Several features in the amino acid composition of antibody-binding peptides were identified. The selective affinity of rituximab increased with an increase in the number of hydrophobic amino acids in a peptide, mainly tryptophan and phenylalanine, while a total charge of the peptide remained relatively small. Peptides with a higher affinity exhibited a lower sum helix propensity. For the 30 strongest peptide binders, a substitutional analysis was performed to determine dissociation constants and the invariant amino acids for binding to rituximab. The strongest selective peptides had a dissociation constant in the hundreds of the nano-molar range. The substitutional analysis revealed a specific hydrophobic epitope for rituximab. To show that conformational binders can, in principle, be detected in array format, cyclic peptide substitutions that are similar to the target of rituximab were investigated. Since the specific binders selected via the resemblance-ranking peptide library were based on the hydrophobic interactions that are widespread in the world of biomolecules, the library can be used to screen for potential linear epitopes that may provide information about the cross-reactivity of antibodies.

Spatially and functionally distinct pools of calmodulin mRNA in cardiac myocytes

Calmodulin (CaM) is a multifunctional calcium-binding protein that modulates activity of many different ion channels, enzymes and other proteins. In vertebrates, CaM represents a unique case, where the exact same, invariant amino acid sequence is encoded by multiple genes on different chromosomes (three in case of the mouse and human). These multiple CaM genes have been hypothesized to afford high spatiotemporal resolution in the control of numerous CaM-dependent processes within highly specialized cells, such as neurons and cardiac myocytes.

Genome-Wide Analysis of Purple Acid Phosphatase Genes in Brassica rapa and Their Association with Pollen Development and Phosphorus Deprivation Stress

PAPs (purple acid phosphatases) belong to the metallo-phosphoesterase superfamily and play important roles in developmental processes, phosphorus foraging, and recycling. However, the specific functions of BrPAPs in Brassica rapa are poorly understood. In this study, 39 BrPAPs were identified and divided into three major clades and nine subgroups. In 8 of the 39 BrPAPs, some invariant amino acid residues were lost or shifted. Based on an expression profiling analysis, BrPAP11, 14, 20, 24, 29, and 34 were specifically expressed in fertile floral buds, indicating their critical roles during pollen development. A total of 21 BrPAPs responded to Pi deprivation in either shoots or roots. Of these, BrPAP4, 5, 19, and 21 were upregulated in roots under Pi depravation conditions, while BrPAP12 was upregulated in the roots in normal conditions. BrPAP28 was upregulated in shoots under Pi depravation conditions, indicating its function shifted compared with its Arabidopsis homolog, AtPAP26. The present work contributes to further investigation of BrPAPs as candidate genes for genetic improvement studies of low phosphorus tolerance as well as for creating male sterile lines based on gene editing methods in Brassica rapa.

Genome-wide identification and evolutionary analysis of MLO gene family in Rosaceae plants

Mutants lacking wild-type MLO (Mildew resistance Locus O) proteins show broad-spectrum resistance to the powdery mildew fungus, and dysregulated cell death control, with spontaneous cell death in response to developmental or abiotic stimuli. In order to understand the evolution and divergence patterns of the MLO gene family in Rosaceae plants, we analysed systematically genome-wide data from Fragaria vesca, Prunus persica, Prunus mume, Malus domestica, Pyrus bretschneideri and Rubus occidentalis based on bioinformatics methods. Using three phylogenetic methods (the neighbour-joining, maximum likelihood, and Bayesian methods), we identified 117 MLO genes from 6 Rosaceae species. The results of all three phylogenetic analysis methods supported that these genes were divided into six clades. Conserved motif analysis found that only motif 2 was present in all MLO proteins and had 3 nearly invariant amino acid residues. The findings indicated that motif 2 might be shared by the MLO gene family. The structural features of these genes showed large variations in sequence length among different species, although the lengths and the numbers of exons exhibited high degrees of similarity. Selective pressure analysis showed extremely significant differences in all 6 clades, with 2, 1, and 1 site(s) under significant positive selection detected in clades III, IV, and VI, respectively. These positive selection sites were important driving forces for the promotion of the functional differentiation of the MLO genes. Functional divergence analysis showed that the significantly divergent sites were located within the domains of the MLO genes. Functional distance analysis showed that the clade V had more conservative functions and might have retained more original functions during the evolutionary process. However, clade I may have undergone extensive altered functional constraints as a specialised functional role. Moreover, the most original function of the MLO genes in Rosaceae could be related to the evolution of their resistance to powdery mildew, which then gradually evolved into functions such as the regulation of flower development, the control of root morphology, and seed evolution due to the different evolutionary rates after gene duplication. These results provide a theoretical basis for further studies of the molecular evolutionary patterns of the plant MLO gene family.

In silico Evolutionary Divergence Analysis Suggests the Potentiality of Capsid Protein VP2 in Serotype-Independent Foot-and-Mouth Disease Virus Detection.

Foot-and-mouth disease (FMD) is an economically devastating disease of the livestock worldwide and caused by the FMD virus (FMDV), which has seven immunologically distinct serotypes (O, A, Asia1, C, and SAT1–SAT3). Studies suggest that VP2 is relatively conserved among three surface-exposed capsid proteins (VP1–VP3) of FMDV, but the level of conservation has not yet been reported. Here we analyzed the comparative evolutionary divergence of VP2 and VP1 to determine the level of conservation in VP2 at different hierarchical levels of three FMDV serotypes (O, A, and Asia1) currently circulating in Asia through an in-depth computational analysis of 14 compiled datasets and designed a consensus VP2 protein that can be used for the development of a serotype-independent FMDV detection tool. The phylogenetic analysis clearly represented a significant level of conservation in VP2 over VP1 at each subgroup level. The protein variability analysis and mutational study showed the presence of 67.4% invariant amino acids in VP2, with the N-terminal end being highly conserved. Nine inter-serotypically conserved fragments located on VP2 have been identified, among which four sites showed promising antigenicity value and surface exposure. The designed 130 amino acid long consensus VP2 protein possessed six surface-exposed B cell epitopes, which suggests the possible potentiality of the protein for the development of a serotype-independent FMDV detection tool in Asia. Conclusively, this is the first study to report the comparative evolutionary divergence between VP2 and VP1, along with proposing the possible potentiality of a designed protein candidate in serotype-independent FMDV detection.

The Phylogenetic Functional Conservation of  <i>Drosophila</i> Seven-in-Absentia (SINA) E3 Ligase and Its Two Human Paralogs, SIAH1 and SIAH2, in  <i>Drosophila</i> Eye Development

Seven-IN-Absentia (SINA) is the most downstream signaling gatekeeper identified in the Drosophila RAS/EGFR signaling pathway. Underscoring the central importance of Drosophila SINA and human SINA homologs (SIAHs) is its phylogenetic conservation across metazoans. SIAH is a major tumor vulnerability in human cancer. To further delineate the SINA function in Drosophila, a genetic modifier screen was conducted, and 28 new sina mutant alleles were isolated. By sequencing analyses, we identified three invariable amino acid residues in SINA’s RING-domain whose single point mutation ablates SINA function . To demonstrate the functional conservation of this medically important family of RING domain E3 ligases, we compared expression of either wild type (WT) or proteolysis-deficient (PD) SINA/SIAH inhibitors of Drosophila SINAWT/PD and human SIAH1WT/PD /2WT/PD under tissue-specific GAL4-drivers in Drosophila. Our results showed that Drosophila SINA and human SIAH1/2 are functionally conserved. The SINAPD/SIAHPD inhibitors are effective in blocking the RAS-dependent neuronal cell fate determination in Drosophila.

Full-length human cytomegalovirus terminase pUL89 adopts a two-domain structure specific for DNA packaging.

A key step in replication of human cytomegalovirus (HCMV) in the host cell is the generation and packaging of unit-length genomes into preformed capsids. The enzymes involved in this process are the terminases. The HCMV terminase complex consists of two terminase subunits, the ATPase pUL56 and the nuclease pUL89. A potential third component pUL51 has been proposed. Even though the terminase subunit pUL89 has been shown to be essential for DNA packaging and interaction with pUL56, it is not known how pUL89 mechanistically achieves sequence-specific DNA binding and nicking. To identify essential domains and invariant amino acids vis-a-vis nuclease activity and DNA binding, alanine substitutions of predicted motifs were analyzed. The analyses indicated that aspartate 463 is an invariant amino acid for the nuclease activity, while argine 544 is an invariant aa for DNA binding. Structural analysis of recombinant protein using electron microscopy in conjunction with single particle analysis revealed a curvilinear monomer with two distinct domains connected by a thinner hinge-like region that agrees well with the predicted structure. These results allow us to model how the terminase subunit pUL89’s structure may mediate its function.

The Universally Conserved Residues Are Not Universally Required for Stable Protein Expression or Functions of Cryptochromes.

Universally conserved residues (UCRs) are invariable amino acids evolutionarily conserved among members of a protein family across diverse kingdoms of life. UCRs are considered important for stability and/or function of protein families, but it has not been experimentally examined systematically. Cryptochromes are photoreceptors in plants or light-independent components of the circadian clocks in mammals. We experimentally analyzed 51 UCRs of Arabidopsis cryptochrome 2 (CRY2) that are universally conserved in eukaryotic cryptochromes from Arabidopsis to human. Surprisingly, we found that UCRs required for stable protein expression of CRY2 in plants are not similarly required for stable protein expression of human hCRY1 in human cells. Moreover, 74% of the stably expressed CRY2 proteins mutated in UCRs retained wild-type-like activities for at least one photoresponses analyzed. Our finding suggests that the evolutionary mechanisms underlying conservation of UCRs or that distinguish UCRs from non-UCRs determining the same functions of individual cryptochromes remain to be investigated.

Probing the role of an invariant active site His in family GH1 β-glycosidases

The reaction mechanism of glycoside hydrolases belonging to family 1 (GH1) of carbohydrate-active enzymes classification, hydrolysing β-O-glycosidic bonds, is well characterised. This family includes several thousands of enzymes with more than 20 different EC numbers depending on the sugar glycone recognised as substrate. Most GH1 β-glycosidases bind their substrates with similar specificity through invariant amino acid residues. Despite extensive studies, the clear identification of the roles played by each of these residues in the recognition of different glycones is not always possible. We demonstrated here that a histidine residue, completely conserved in the active site of the enzymes of this family, interacts with the C2-OH of the substrate in addition to the C3-OH as previously shown by 3 D-structure determination.

Genome-Wide Assessment of Putative Superoxide Dismutases in Unicellular and Filamentous Cyanobacteria

Abstract Cyanobacteria are photoautotrophic prokaryotes capable to grow in diverse ecological habitats, originated 2.5-3.5 billion years ago and were first to produce oxygen. Since then superoxide dismutases (SOD) acquired great significance due to their ability to catalyze detoxification of byproducts of oxygenic photosynthesis i.e. superoxide radicals. In the present study, we extracted information regarding SODs from species of sequenced cyanobacteria and investigated their diversity, conservation, domain structure, and evolution. 144 putative SOD homologs were identified. Unlike other protein families (ex. serine-threonine kinases) SODs are present in all cyanobacterial species reflecting their significant role in survival. However, their distribution varies fewer (0.01%-0.09%) found in unicellular marine strains whereas abundant (0.02%-0.07%) in filamentous nitrogen-fixing cyanobacteria. They were classified into three major subfamilies according to their domain structures: Fe/MnSOD, Cu/ZnSOD and NiSOD. Interestingly, they lack additional domains as found in proteins of other families however motifs and invariant amino acids typical in eukaryotic SODs were conserved well in these proteins indicating similar catalytic mechanism as eukaryotic SODs. Phylogenetic relationships correspond well with phylogenies based on 16S rRNA and clustering occurs on the basis of structural characteristics such as domain organization. Gene gain-and-loss is insignificant during SOD evolution as evidenced by the absence of additional domain. This study has not only examined an overall background of sequence-structure-function interactions for the SOD gene family but also revealed variation among SOD distribution based on ecophysiological and morphological characters.

Crystal structure of glycosyltrehalose synthase from Sulfolobus shibatae DSM5389.

Glycosyltrehalose synthase (GTSase) converts the glucosidic bond between the last two glucose residues of amylose from an α-1,4 bond to an α-1,1 bond, generating a nonreducing glycosyl trehaloside, in the first step of the biosynthesis of trehalose. To better understand the structural basis of the catalytic mechanism, the crystal structure of GTSase from the hyperthermophilic archaeon Sulfolobus shibatae DSM5389 (5389-GTSase) has been determined to 2.4 Å resolution by X-ray crystallography. The structure of 5389-GTSase can be divided into five domains. The central domain contains the (β/α)8-barrel fold that is conserved as the catalytic domain in the α-amylase family. Three invariant catalytic carboxylic amino acids in the α-amylase family are also found in GTSase at positions Asp241, Glu269 and Asp460 in the catalytic domain. The shape of the catalytic cavity and the pocket size at the bottom of the cavity correspond to the intramolecular transglycosylation mechanism proposed from previous enzymatic studies.

Structural determinants for protein unfolding and translocation by the Hsp104 protein disaggregase.

The ring-forming Hsp104 ATPase cooperates with Hsp70 and Hsp40 molecular chaperones to rescue stress-damaged proteins from both amorphous and amyloid-forming aggregates. The ability to do so relies upon pore loops present in the first ATP-binding domain (AAA-1; loop-1 and loop-2 ) and in the second ATP-binding domain (AAA-2; loop-3) of Hsp104, which face the protein translocating channel and couple ATP-driven changes in pore loop conformation to substrate translocation. A hallmark of loop-1 and loop-3 is an invariable and mutational sensitive aromatic amino acid (Tyr257 and Tyr662) involved in substrate binding. However, the role of conserved aliphatic residues (Lys256, Lys258, and Val663) flanking the pore loop tyrosines, and the function of loop-2 in protein disaggregation has not been investigated. Here we present the crystal structure of an N-terminal fragment of Saccharomyces cerevisiae Hsp104 exhibiting molecular interactions involving both AAA-1 pore loops, which resemble contacts with bound substrate. Corroborated by biochemical experiments and functional studies in yeast, we show that aliphatic residues flanking Tyr257 and Tyr662 are equally important for substrate interaction, and abolish Hsp104 function when mutated to glycine. Unexpectedly, we find that loop-2 is sensitive to aspartate substitutions that impair Hsp104 function and abolish protein disaggregation when loop-2 is replaced by four aspartate residues. Our observations suggest that Hsp104 pore loops have non-overlapping functions in protein disaggregation and together coordinate substrate binding, unfolding, and translocation through the Hsp104 hexamer.

Phylogenetic analysis of the SINA/SIAH ubiquitin E3 ligase family in Metazoa

BackgroundThe RAS signaling pathway is a pivotal developmental pathway that controls many fundamental biological processes including cell proliferation, differentiation, movement and apoptosis. Drosophila Seven-IN-Absentia (SINA) is a ubiquitin E3 ligase that is the most downstream signaling “gatekeeper” whose biological activity is essential for proper RAS signal transduction. Vertebrate SINA homologs (SIAHs) share a high degree of amino acid identity with that of Drosophila SINA. SINA/SIAH is the most conserved signaling component in the canonical EGFR/RAS/RAF/MAPK signal transduction pathway.ResultsVertebrate SIAH1, 2, and 3 are the three orthologs to invertebrate SINA protein. SINA and SIAH1 orthologs are found in all major taxa of metazoans. These proteins have four conserved functional domains, known as RING (Really Interesting New Gene), SZF (SIAH-type zinc finger), SBS (substrate binding site) and DIMER (Dimerization). In addition to the siah1 gene, most vertebrates encode two additional siah genes (siah2 and siah3) in their genomes. Vertebrate SIAH2 has a highly divergent and extended N-terminal sequence, while its RING, SZF, SBS and DIMER domains maintain high amino acid identity/similarity to that of SIAH1. But unlike vertebrate SIAH1 and SIAH2, SIAH3 lacks a functional RING domain, suggesting that SIAH3 may be an inactive E3 ligase. The SIAH3 subtree exhibits a high degree of amino acid divergence when compared to the SIAH1 and SIAH2 subtrees. We find that SIAH1 and SIAH2 are expressed in all human epithelial cell lines examined thus far, while SIAH3 is only expressed in a limited subset of cancer cell lines.ConclusionThrough phylogenetic analyses of metazoan SINA and SIAH E3 ligases, we identified many invariant and divergent amino acid residues, as well as the evolutionarily conserved functional motifs in this medically relevant gene family. Our phylomedicinal study of this unique metazoan SINA/SIAH protein family has provided invaluable evolution-based support towards future effort to design logical, potent, and durable anti-SIAH-based anticancer strategies against oncogenic K-RAS-driven metastatic human cancers. Thus, this method of evolutionary study should be of interest in cancer biology.

Correlation between antigenicity and variability in the vls antigenic variation system of Borrelia burgdorferi

Many parasites have evolved antigenic variation systems that alter surface proteins in order to evade recognition by presently expressed antibodies and subsequent death. Although the amino acid positions in antigens to which antibodies most commonly target are expected to be the most variable, this assumption has not been investigated. Using the vls antigenic variation system of Borrelia burgdorferi as a model, we first investigated this assumption computationally and then developed a sensitive immunoassay to experimentally validate the computational results. There was a strong correlation between variability at an amino acid position and each of the computational metrics associated with antibody reactivity. However, empirical measures of antibody reactivity were not consistently greater at the variable amino acid positions than at the invariant amino acid positions. The inconsistent experimental support for this hypothesis suggests that the biological effect of variability at an amino acid position is obfuscated by other factors.

Insight into the phylogenetic relationship and structural features of vertebrate myoglobin family

Myoglobin (Mb) is studied to clarify the structure-function relationships in protein science. In this work, we report the results of a comparative analysis of amino acid sequences from 298 vertebrate Mbs. Forty-one high conserved residues were identified and seven of them were invariants [E18, G25, F43, V68, L72, H93 (proximal histidine) and H97]. E18 is the only invariant amino acid residue located out of the heme-pocket and Xe-cavities playing a role in interaction between the A and E-helices. A comparative analysis of several parameters related to amino acid composition shows an increase of average mass, accessible surface area and volume per residue from Actinopterygii to Mammalia and Aves. This may be due to an increased number of bulky residues reducing the non-specific cavities volume and thus improving the oxygen flow between the heme site and the outside of the protein. Finally, the phylogenetic analyses of Mb in vertebrates are consistent with an evolution that runs with the diversification of the species, but in which several episodes of gene duplication and lost have occurred, less frequently in the ancestors of great taxons, cartilaginous fishes and non-avian reptiles, most frequently in ray-finned fishes and mammals, and very frequently in birds.

The role of WRKY transcription factors in plants

The WRKY proteins are one of the largest families of transcriptional regulators in plants. Their biosynthesis is induced during certain stages of plant development and upon pathogen infection. A single WRKY transcription factor may be involved in the regulation of several apparently disparate processes. Their hallmark is strong conservation of the DNA binding domain which contains an invariant WRKYGQK amino acid sequence and a zinc binding motif. However, the overall sequences of individual representatives are highly divergent. Little is known about the 3D-structure of the WRKY proteins. Up to date there have only been results of structural strudies of DNA binding domain available. In this review, the biological function as well the structural studies of the WRKY proteins have been recapitulated.

A bioinformatic survey of distribution, conservation, and probable functions of LuxR solo regulators in bacteria.

LuxR solo transcriptional regulators contain both an autoinducer binding domain (ABD; N-terminal) and a DNA binding Helix-Turn-Helix domain (HTH; C-terminal), but are not associated with a cognate N-acyl homoserine lactone (AHL) synthase coding gene in the same genome. Although a few LuxR solos have been characterized, their distributions as well as their role in bacterial signal perception and other processes are poorly understood. In this study we have carried out a systematic survey of distribution of all ABD containing LuxR transcriptional regulators (QS domain LuxRs) available in the InterPro database (IPR005143), and identified those lacking a cognate AHL synthase. These LuxR solos were then analyzed regarding their taxonomical distribution, predicted functions of neighboring genes and the presence of complete AHL-QS systems in the genomes that carry them. Our analyses reveal the presence of one or multiple predicted LuxR solos in many proteobacterial genomes carrying QS domain LuxRs, some of them harboring genes for one or more AHL-QS circuits. The presence of LuxR solos in bacteria occupying diverse environments suggests potential ecological functions for these proteins beyond AHL and interkingdom signaling. Based on gene context and the conservation levels of invariant amino acids of ABD, we have classified LuxR solos into functionally meaningful groups or putative orthologs. Surprisingly, putative LuxR solos were also found in a few non-proteobacterial genomes which are not known to carry AHL-QS systems. Multiple predicted LuxR solos in the same genome appeared to have different levels of conservation of invariant amino acid residues of ABD questioning their binding to AHLs. In summary, this study provides a detailed overview of distribution of LuxR solos and their probable roles in bacteria with genome sequence information.