Systematic Dissection and Trajectory-Scanning Mutagenesis of the Molecular Interface That Ensures Specificity of Two-Component Signaling Pathways

Two-component signal transduction systems enable bacteria to sense and respond to a wide range of environmental stimuli. Sensor histidine kinases transmit signals to their cognate response regulators via phosphorylation. The faithful transmission of information through two-component pathways and the avoidance of unwanted cross-talk require exquisite specificity of histidine kinase-response regulator interactions to ensure that cells mount the appropriate response to external signals. To identify putative specificity-determining residues, we have analyzed amino acid coevolution in two-component proteins and identified a set of residues that can be used to rationally rewire a model signaling pathway, EnvZ-OmpR. To explore how a relatively small set of residues can dictate partner selectivity, we combined alanine-scanning mutagenesis with an approach we call trajectory-scanning mutagenesis, in which all mutational intermediates between the specificity residues of EnvZ and another kinase, RstB, were systematically examined for phosphotransfer specificity. The same approach was used for the response regulators OmpR and RstA. Collectively, the results begin to reveal the molecular mechanism by which a small set of amino acids enables an individual kinase to discriminate amongst a large set of highly-related response regulators and vice versa. Our results also suggest that the mutational trajectories taken by two-component signaling proteins following gene or pathway duplication may be constrained and subject to differential selective pressures. Only some trajectories allow both the maintenance of phosphotransfer and the avoidance of unwanted cross-talk

Microfluidics with fluid walls

Microfluidics has great potential, but the complexity of fabricating and operating devices has limited its use. Here we describe a method - Freestyle Fluidics - that overcomes many key limitations. In this method, liquids are confined by fluid (not solid) walls. Aqueous circuits with any 2D shape are printed in seconds on plastic or glass Petri dishes; then, interfacial forces pin liquids to substrates, and overlaying an immiscible liquid prevents evaporation. Confining fluid walls are pliant and resilient; they self-heal when liquids are pipetted through them. We drive flow through a wide range of circuits passively by manipulating surface tension and hydrostatic pressure, and actively using external pumps. Finally, we validate the technology with two challenging applications - triggering an inflammatory response in human cells and chemotaxis in bacterial biofilms. This approach provides a powerful and versatile alternative to traditional microfluidics.The complexity of fabricating and operating microfluidic devices limits their use. Walsh et al. describe a method in which circuits are printed as quickly and simply as writing with a pen, and liquids in them are confined by fluid instead of solid walls

Robust Signal Processing in Living Cells

Cellular signaling networks have evolved an astonishing ability to function reliably and with high fidelity in uncertain environments. A crucial prerequisite for the high precision exhibited by many signaling circuits is their ability to keep the concentrations of active signaling compounds within tightly defined bounds, despite strong stochastic fluctuations in copy numbers and other detrimental influences. Based on a simple mathematical formalism, we identify topological organizing principles that facilitate such robust control of intracellular concentrations in the face of multifarious perturbations. Our framework allows us to judge whether a multiple-input-multiple-output reaction network is robust against large perturbations of network parameters and enables the predictive design of perfectly robust synthetic network architectures. Utilizing the Escherichia coli chemotaxis pathway as a hallmark example, we provide experimental evidence that our framework indeed allows us to unravel the topological organization of robust signaling. We demonstrate that the specific organization of the pathway allows the system to maintain global concentration robustness of the diffusible response regulator CheY with respect to several dominant perturbations. Our framework provides a counterpoint to the hypothesis that cellular function relies on an extensive machinery to fine-tune or control intracellular parameters. Rather, we suggest that for a large class of perturbations, there exists an appropriate topology that renders the network output invariant to the respective perturbations

Coupling of protein localization and cell movements by a dynamically localized response regulator in Myxococcus xanthus

Myxococcus xanthus cells harbor two motility machineries, type IV pili (Tfp) and the A-engine. During reversals, the two machineries switch polarity synchronously. We present a mechanism that synchronizes this polarity switching. We identify the required for motility response regulator (RomR) as essential for A-motility. RomR localizes in a bipolar, asymmetric pattern with a large cluster at the lagging cell pole. The large RomR cluster relocates to the new lagging pole in parallel with cell reversals. Dynamic RomR localization is essential for cell reversals, suggesting that RomR relocalization induces the polarity switching of the A-engine. The analysis of RomR mutants shows that the output domain targets RomR to the poles and the receiver domain is essential for dynamic localization. The small GTPase MglA establishes correct RomR polarity, and the Frz two-component system regulates dynamic RomR localization. FrzS localizes with Tfp at the leading pole and relocates in an Frz-dependent manner to the opposite pole during reversals; FrzS and RomR localize and oscillate independently. The Frz system synchronizes these oscillations and thus the synchronous polarity switching of the motility machineries

The genetic organisation of prokaryotic two-component system signalling pathways

<p>Abstract</p> <p>Background</p> <p>Two-component systems (TCSs) are modular and diverse signalling pathways, involving a stimulus-responsive transfer of phosphoryl groups from transmitter to partner receiver domains. TCS gene and domain organisation are both potentially informative regarding biological function, interaction partnerships and molecular mechanisms. However, there is currently little understanding of the relationships between domain architecture, gene organisation and TCS pathway structure.</p> <p>Results</p> <p>Here we classify the gene and domain organisation of TCS gene loci from 1405 prokaryotic replicons (>40,000 TCS proteins). We find that 200 bp is the most appropriate distance cut-off for defining whether two TCS genes are functionally linked. More than 90% of all TCS gene loci encode just one or two transmitter and/or receiver domains, however numerous other geometries exist, often with large numbers of encoded TCS domains. Such information provides insights into the distribution of TCS domains between genes, and within genes. As expected, the organisation of TCS genes and domains is affected by phylogeny, and plasmid-encoded TCS exhibit differences in organisation from their chromosomally-encoded counterparts.</p> <p>Conclusions</p> <p>We provide here an overview of the genomic and genetic organisation of TCS domains, as a resource for further research. We also propose novel metrics that build upon TCS gene/domain organisation data and allow comparisons between genomic complements of TCSs. In particular, '<it>percentage orphaned TCS genes</it>' (or 'Dissemination') and '<it>percentage of complex loci</it>' (or 'Sophistication') appear to be useful discriminators, and to reflect mechanistic aspects of TCS organisation not captured by existing metrics.</p

Exopolysaccharide-Independent Social Motility of Myxococcus xanthus

Social motility (S motility), the coordinated movement of large cell groups on agar surfaces, of Myxococcus xanthus requires type IV pili (TFP) and exopolysaccharides (EPS). Previous models proposed that this behavior, which only occurred within cell groups, requires cycles of TFP extension and retraction triggered by the close interaction of TFP with EPS. However, the curious observation that M. xanthus can perform TFP-dependent motility at a single-cell level when placed onto polystyrene surfaces in a highly viscous medium containing 1% methylcellulose indicated that “S motility” is not limited to group movements. In an apparent further challenge of the previous findings for S motility, mutants defective in EPS production were found to perform TFP-dependent motility on polystyrene surface in methylcellulose-containing medium. By exploring the interactions between pilin and surface materials, we found that the binding of TFP onto polystyrene surfaces eliminated the requirement for EPS in EPS- cells and thus enabled TFP-dependent motility on a single cell level. However, the presence of a general anchoring surface in a viscous environment could not substitute for the role of cell surface EPS in group movement. Furthermore, EPS was found to serve as a self-produced anchoring substrate that can be shed onto surfaces to enable cells to conduct TFP-dependent motility regardless of surface properties. These results suggested that in certain environments, such as in methylcellulose solution, the cells could bypass the need for EPS to anchor their TPF and conduct single-cell S motility to promote exploratory movement of colonies over new specific surfaces

Inference of Co-Evolving Site Pairs: an Excellent Predictor of Contact Residue Pairs in Protein 3D structures

Residue-residue interactions that fold a protein into a unique three-dimensional structure and make it play a specific function impose structural and functional constraints on each residue site. Selective constraints on residue sites are recorded in amino acid orders in homologous sequences and also in the evolutionary trace of amino acid substitutions. A challenge is to extract direct dependences between residue sites by removing indirect dependences through other residues within a protein or even through other molecules. Recent attempts of disentangling direct from indirect dependences of amino acid types between residue positions in multiple sequence alignments have revealed that the strength of inferred residue pair couplings is an excellent predictor of residue-residue proximity in folded structures. Here, we report an alternative attempt of inferring co-evolving site pairs from concurrent and compensatory substitutions between sites in each branch of a phylogenetic tree. First, branch lengths of a phylogenetic tree inferred by the neighbor-joining method are optimized as well as other parameters by maximizing a likelihood of the tree in a mechanistic codon substitution model. Mean changes of quantities, which are characteristic of concurrent and compensatory substitutions, accompanied by substitutions at each site in each branch of the tree are estimated with the likelihood of each substitution. Partial correlation coefficients of the characteristic changes along branches between sites are calculated and used to rank co-evolving site pairs. Accuracy of contact prediction based on the present co-evolution score is comparable to that achieved by a maximum entropy model of protein sequences for 15 protein families taken from the Pfam release 26.0. Besides, this excellent accuracy indicates that compensatory substitutions are significant in protein evolution.Comment: 17 pages, 4 figures, and 4 tables with supplementary information of 5 figure

Temporal Controls of the Asymmetric Cell Division Cycle in Caulobacter crescentus

The asymmetric cell division cycle of Caulobacter crescentus is orchestrated by an elaborate gene-protein regulatory network, centered on three major control proteins, DnaA, GcrA and CtrA. The regulatory network is cast into a quantitative computational model to investigate in a systematic fashion how these three proteins control the relevant genetic, biochemical and physiological properties of proliferating bacteria. Different controls for both swarmer and stalked cell cycles are represented in the mathematical scheme. The model is validated against observed phenotypes of wild-type cells and relevant mutants, and it predicts the phenotypes of novel mutants and of known mutants under novel experimental conditions. Because the cell cycle control proteins of Caulobacter are conserved across many species of alpha-proteobacteria, the model we are proposing here may be applicable to other genera of importance to agriculture and medicine (e.g., Rhizobium, Brucella)

Validation of Coevolving Residue Algorithms via Pipeline Sensitivity Analysis: ELSC and OMES and ZNMI, Oh My!

Correlated amino acid substitution algorithms attempt to discover groups of residues that co-fluctuate due to either structural or functional constraints. Although these algorithms could inform both ab initio protein folding calculations and evolutionary studies, their utility for these purposes has been hindered by a lack of confidence in their predictions due to hard to control sources of error. To complicate matters further, naive users are confronted with a multitude of methods to choose from, in addition to the mechanics of assembling and pruning a dataset. We first introduce a new pair scoring method, called ZNMI (Z-scored-product Normalized Mutual Information), which drastically improves the performance of mutual information for co-fluctuating residue prediction. Second and more important, we recast the process of finding coevolving residues in proteins as a data-processing pipeline inspired by the medical imaging literature. We construct an ensemble of alignment partitions that can be used in a cross-validation scheme to assess the effects of choices made during the procedure on the resulting predictions. This pipeline sensitivity study gives a measure of reproducibility (how similar are the predictions given perturbations to the pipeline?) and accuracy (are residue pairs with large couplings on average close in tertiary structure?). We choose a handful of published methods, along with ZNMI, and compare their reproducibility and accuracy on three diverse protein families. We find that (i) of the algorithms tested, while none appear to be both highly reproducible and accurate, ZNMI is one of the most accurate by far and (ii) while users should be wary of predictions drawn from a single alignment, considering an ensemble of sub-alignments can help to determine both highly accurate and reproducible couplings. Our cross-validation approach should be of interest both to developers and end users of algorithms that try to detect correlated amino acid substitutions

Cellular Levels and Binding of c-di-GMP Control Subcellular Localization and Activity of the Vibrio cholerae Transcriptional Regulator VpsT

The second messenger, cyclic diguanylate (c-di-GMP), regulates diverse cellular processes in bacteria. C-di-GMP is produced by diguanylate cyclases (DGCs), degraded by phosphodiesterases (PDEs), and receptors couple c-di-GMP production to cellular responses. In many bacteria, including Vibrio cholerae, multiple DGCs and PDEs contribute to c-di-GMP signaling, and it is currently unclear whether the compartmentalization of c-di-GMP signaling components is required to mediate c-di-GMP signal transduction. In this study we show that the transcriptional regulator, VpsT, requires c-di-GMP binding for subcellular localization and activity. Only the additive deletion of five DGCs markedly decreases the localization of VpsT, while single deletions of each DGC do not impact VpsT localization. Moreover, mutations in residues required for c-di-GMP binding, c-di-GMP-stabilized dimerization and DNA binding of VpsT abrogate wild type localization and activity. VpsT does not co-localize or interact with DGCs suggesting that c-di-GMP from these DGCs diffuses to VpsT, supporting a model in which c-di-GMP acts at a distance. Furthermore, VpsT localization in a heterologous host, Escherichia coli, requires a catalytically active DGC and is enhanced by the presence of VpsT-target sequences. Our data show that c-di-GMP signaling can be executed through an additive cellular c-di-GMP level from multiple DGCs affecting the localization and activity of a c-di-GMP receptor and furthers our understanding of the mechanisms of second messenger signaling