Co-evolution at protein-protein interfaces guides inference of stoichiometry of oligomeric protein complexes by de novo structure prediction.

The quaternary structure with specific stoichiometry is pivotal to the specific function of protein complexes. However, determining the structure of many protein complexes experimentally remains a major bottleneck. Structural bioinformatics approaches, such as the deep learning algorithm Alphafold2-multimer (AF2-multimer), leverage the co-evolution of amino acids and sequence-structure relationships for accurate de novo structure and contact prediction. Pseudo-likelihood maximization direct coupling analysis (plmDCA) has been used to detect co-evolving residue pairs by statistical modeling. Here, we provide evidence that combining both methods can be used for de novo prediction of the quaternary structure and stoichiometry of a protein complex. We achieve this by augmenting the existing AF2-multimer confidence metrics with an interpretable score to identify the complex with an optimal fraction of native contacts of co-evolving residue pairs at intermolecular interfaces. We use this strategy to predict the quaternary structure and non-trivial stoichiometries of Bacillus subtilis spore germination protein complexes with unknown structures. Co-evolution at intermolecular interfaces may therefore synergize with AI-based de novo quaternary structure prediction of structurally uncharacterized bacterial protein complexes.

Methane formation driven by reactive oxygen species across all living organisms.

Methane (CH4), the most abundant hydrocarbon in the atmosphere, originates largely from biogenic sources1 linked to an increasing number of organisms occurring in oxic and anoxic environments. Traditionally, biogenic CH4 has been regarded as the final product of anoxic decomposition of organic matter by methanogenic archaea. However, plants2,3, fungi4, algae5 and cyanobacteria6 can produce CH4 in the presence of oxygen. Although methanogens are known to produce CH4 enzymatically during anaerobic energy metabolism7, the requirements and pathways for CH4 production by non-methanogenic cells are poorly understood. Here, we demonstrate that CH4 formation by Bacillus subtilis and Escherichia coli is triggered by free iron and reactive oxygen species (ROS), which are generated by metabolic activity and enhanced by oxidative stress. ROS-induced methyl radicals, which are derived from organic compounds containing sulfur- or nitrogen-bonded methyl groups, are key intermediates that ultimately lead to CH4 production. We further show CH4 production by many other model organisms from the Bacteria, Archaea and Eukarya domains, including in several human cell lines. All these organisms respond to inducers of oxidative stress by enhanced CH4 formation. Our results imply that all living cells probably possess a common mechanism of CH4 formation that is based on interactions among ROS, iron and methyl donors, opening new perspectives for understanding biochemical CH4 formation and cycling.

Pervasive prophage recombination occurs during evolution of spore-forming Bacilli.

Phages are the main source of within-species bacterial diversity and drivers of horizontal gene transfer, but we know little about the mechanisms that drive genetic diversity of these mobile genetic elements (MGEs). Recently, we showed that a sporulation selection regime promotes evolutionary changes within SPß prophage of Bacillus subtilis, leading to direct antagonistic interactions within the population. Herein, we reveal that under a sporulation selection regime, SPß recombines with low copy number phi3Ts phage DNA present within the B. subtilis population. Recombination results in a new prophage occupying a different integration site, as well as the spontaneous release of virulent phage hybrids. Analysis of Bacillus sp. strains suggests that SPß and phi3T belong to a distinct cluster of unusually large phages inserted into sporulation-related genes that are equipped with a spore-related genetic arsenal. Comparison of Bacillus sp. genomes indicates that similar diversification of SPß-like phages takes place in nature. Our work is a stepping stone toward empirical studies on phage evolution, and understanding the eco-evolutionary relationships between bacteria and their phages. By capturing the first steps of new phage evolution, we reveal striking relationship between survival strategy of bacteria and evolution of their phages.

A spore quality-quantity tradeoff favors diverse sporulation strategies in Bacillus subtilis.

Quality-quantity tradeoffs govern the production of propagules across taxa and can explain variability in life-history traits in higher organisms. A quality-quantity tradeoff was recently discovered in spore forming bacteria, but whether it impacts fitness is unclear. Here we show both theoretically and experimentally that the nutrient supply during spore revival determines the fitness advantage associated with different sporulation behaviors in Bacillus subtilis. By tuning sporulation rates we generate spore-yield and spore-quality strategists that compete with each other in a microscopic life-cycle assay. The quality (yield) strategist is favored when spore revival is triggered by poor (rich) nutrients. We also show that natural isolates from the gut and soil employ different life-cycle strategies that result from genomic variations in the number of rap-phr signaling systems. Taken together, our results suggest that a spore quality-quantity tradeoff contributes to the evolutionary adaptation of sporulating bacteria.

Ratiometric population sensing by a pump-probe signaling system in Bacillus subtilis.

Communication by means of diffusible signaling molecules facilitates higher-level organization of cellular populations. Gram-positive bacteria frequently use signaling peptides, which are either detected at the cell surface or 'probed' by intracellular receptors after being pumped into the cytoplasm. While the former type is used to monitor cell density, the functions of pump-probe networks are less clear. Here we show that pump-probe networks can, in principle, perform different tasks and mediate quorum-sensing, chronometric and ratiometric control. We characterize the properties of the prototypical PhrA-RapA system in Bacillus subtilis using FRET. We find that changes in extracellular PhrA concentrations are tracked rather poorly; instead, cells accumulate and strongly amplify the signal in a dose-dependent manner. This suggests that the PhrA-RapA system, and others like it, have evolved to sense changes in the composition of heterogeneous populations and infer the fraction of signal-producing cells in a mixed population to coordinate cellular behaviors.

Phenotypic memory in Bacillus subtilis links dormancy entry and exit by a spore quantity-quality tradeoff.

Some bacteria, such as Bacillus subtilis, withstand starvation by forming dormant spores that revive when nutrients become available. Although sporulation and spore revival jointly determine survival in fluctuating environments, the relationship between them has been unclear. Here we show that these two processes are linked by a phenotypic "memory" that arises from a carry-over of molecules from the vegetative cell into the spore. By imaging life histories of individual B. subtilis cells using fluorescent reporters, we demonstrate that sporulation timing controls nutrient-induced spore revival. Alanine dehydrogenase contributes to spore memory and controls alanine-induced outgrowth, thereby coupling a spore's revival capacity to the gene expression and growth history of its progenitors. A theoretical analysis, and experiments with signaling mutants exhibiting altered sporulation timing, support the hypothesis that such an intrinsically generated memory leads to a tradeoff between spore quantity and spore quality, which could drive the emergence of complex microbial traits.

Molecular and cellular factors control signal transduction via switchable allosteric modulator proteins (SAMPs).

BACKGROUND: Rap proteins from Bacilli directly target response regulators of bacterial two-component systems and modulate their activity. Their effects are controlled by binding of signaling peptides to an allosteric site. Hence Raps exemplify a class of monomeric signaling receptors, which we call switchable allosteric modulator proteins (SAMPs). These proteins have potential applications in diverse biomedical and biotechnical settings, but a quantitative understanding of the impact of molecular and cellular factors on signal transduction is lacking. Here we introduce mathematical models that elucidate how signals are propagated though the network upon receptor stimulation and control the level of active response regulator. RESULTS: Based on a systematic parameter analysis of the models, we show that key features of the dose-response behavior at steady state are controlled either by the molecular properties of the modulator or the signaling context. In particular, we find that the biochemical activity (i.e. non-enzymatic vs. enzymatic) and allosteric properties of the modulator control the response amplitude. The Hill coefficient and the EC50 are controlled in addition by the relative ligand affinities. By tuning receptor properties, either graded or more switch-like (memory-less) response functions can be fashioned. Furthermore, we show that other contextual factors (e.g. relative concentrations of network components and kinase activity) have a substantial impact on the response, and we predict that there exists a modulator concentration which is optimal for response amplitude. CONCLUSION: We discuss data on Rap-Phr systems in B. subtilis to show how our models can contribute to an integrated view of SAMP signaling by combining biochemical, structural and physiological insights. Our results also suggest that SAMPs could be evolved or engineered to implement diverse response behaviors. However-without additional regulatory controls-they can generate rather variable cellular outputs.

The quorum-sensing regulator ComA from Bacillus subtilis activates transcription using topologically distinct DNA motifs.

ComA-like transcription factors regulate the quorum response in numerous Gram-positive bacteria. ComA proteins belong to the tetrahelical helix-turn-helix superfamily of transcriptional activators, which bind as homodimers to inverted sequence repeats in the DNA. Here, we report that ComA from Bacillus subtilis recognizes a topologically distinct motif, in which the binding elements form a direct repeat. We provide in vitro and in vivo evidence that the canonical and non-canonical site play an important role in facilitating type I and type II promoter activation, respectively, by interacting with different subunits of RNA polymerase. We furthermore show that there is a variety of contexts in which the non-canonical site can occur and identify new direct target genes that are located within the integrative and conjugative element ICEBs1. We therefore suggest that ComA acts as a multifunctional transcriptional activator and provides a striking example for complexity in protein-DNA interactions that evolved in the context of quorum sensing.

Metabolic cross-feeding via intercellular nanotubes among bacteria.

Bacteria frequently exchange metabolites by diffusion through the extracellular environment, yet it remains generally unclear whether bacteria can also use cell-cell connections to directly exchange nutrients. Here we address this question by engineering cross-feeding interactions within and between Acinetobacter baylyi and Escherichia coli, in which two distant bacterial species reciprocally exchange essential amino acids. We establish that in a well-mixed environment E. coli, but likely not A. baylyi, can connect to other bacterial cells via membrane-derived nanotubes and use these to exchange cytoplasmic constituents. Intercellular connections are induced by auxotrophy-causing mutations and cease to establish when amino acids are externally supplied. Electron and fluorescence microscopy reveal a network of nanotubular structures that connects bacterial cells and enables an intercellular transfer of cytoplasmic materials. Together, our results demonstrate that bacteria can use nanotubes to exchange nutrients among connected cells and thus help to distribute metabolic functions within microbial communities.

A modular view of the diversity of cell-density-encoding schemes in bacterial quorum-sensing systems.

Certain environmental parameters are accessible to cells only indirectly and require an encoding step for cells to retrieve the relevant information. A prominent example is the phenomenon of quorum sensing by microorganisms, where information about cell density is encoded by means of secreted signaling molecules. The mapping of cell density to signal molecule concentration and the corresponding network modules involved have been at least partially characterized in many bacteria, and vary markedly between different systems. In this study, we investigate theoretically how differences in signal transport, signal modification, and site of signal detection shape the encoding function and affect the sensitivity and the noise characteristics of the cell-density-encoding process. We find that different modules are capable of implementing both fairly basic as well as more complex encoding schemes, whose qualitative characteristics vary with cell density and are linked to network architecture, providing the basis for a hierarchical classification scheme. We exploit the tight relationship between encoding behavior and network architecture to constrain the network topology of partially characterized natural systems, and verify one such prediction by showing experimentally that Vibrio harveyi is capable of importing Autoinducer 2. The framework developed in this research can serve not only to guide reverse engineering of natural systems but also to stimulate the design of synthetic systems and generally facilitate a better understanding of the complexities arising in the quorum-sensing process because of variations in the physical organization of the encoder network module.

Ectopic integration vectors for generating fluorescent promoter fusions in Bacillus subtilis with minimal dark noise.

Fluorescent protein promoter reporters are important tools that are widely used for diverse purposes in microbiology, systems biology and synthetic biology and considerable engineering efforts are still geared at improving the sensitivity of the reporter systems. Here we focus on dark noise, i.e. the signal that is generated by the empty vector control. We quantitatively characterize the dark noise of a few common bacterial reporter systems by single cell microscopy. All benchmarked reporter systems generated significant amounts of dark noise that exceed the cellular autofluorescence to different extents. We then reengineered a multicolor set of fluorescent ectopic integration vectors for Bacillus subtilis by introducing a terminator immediately upstream of the promoter insertion site, resulting in an up to 2.7-fold reduction of noise levels. The sensitivity and dynamic range of the new high-performance pXFP_Star reporter system is only limited by cellular autofluorescence. Moreover, based on studies of the rapE promoter of B. subtilis we show that the new pXFP_Star reporter system reliably reports on the weak activity of the rapE promoter whereas the original reporter system fails because of transcriptional interference. Since the pXFP_Star reporter system properly isolates the promoter from spurious transcripts, it is a particularly suitable tool for quantitative characterization of weak promoters in B. subtilis.

Effect of adhesion geometry and rigidity on cellular force distributions.

The behavior and fate of tissue cells are controlled by the rigidity and geometry of their adhesive environment, possibly through forces localized to sites of adhesion. We introduce a mechanical model that predicts cellular force distributions for cells adhering to adhesive patterns with different geometries and rigidities. For continuous adhesion along a closed contour, forces are predicted to be localized to the corners. For discrete sites of adhesion, the model predicts the forces to be mainly determined by the lateral pull of the cell contour. With increasing distance between two neighboring sites of adhesion, the adhesion force increases because the cell shape results in steeper pulling directions. Softer substrates result in smaller forces. Our predictions agree well with experimental force patterns measured on pillar assays.

Complexity in bacterial cell-cell communication: quorum signal integration and subpopulation signaling in the Bacillus subtilis phosphorelay.

A common form of quorum sensing in gram-positive bacteria is mediated by peptides that act as phosphatase regulators (Phr) of receptor aspartyl phosphatases (Raps). In Bacillus subtilis, several Phr signals are integrated in sporulation phosphorelay signal transduction. We theoretically demonstrate that the phosphorelay can act as a computational machine performing a sensitive division operation of kinase-encoded signals by quorum-modulated Rap signals, indicative of cells computing a "food per cell" estimate to decide whether to enter sporulation. We predict expression from the rapA-phrA operon to bifurcate as relative environmental signals change in a developing population. We experimentally observe that the rapA-phrA operon is heterogeneously induced in sporulating microcolonies. Uninduced cells sporulate rather synchronously early on, whereas the RapA/PhrA subpopulation sporulates less synchronously throughout later stationary phase. Moreover, we show that cells sustain PhrA expression during periods of active growth. Together with the model, these findings suggest that the phosphorelay may normalize environmental signals by the size of the (sub)population actively competing for nutrients (as signaled by PhrA). Generalizing this concept, the various Phrs could facilitate subpopulation communication in dense isogenic communities to control the physiological strategies followed by differentiated subpopulations by interpreting (environmental) signals based on the spatiotemporal community structure.

Filamentous network mechanics and active contractility determine cell and tissue shape.

For both cells and tissues, shape is closely correlated with function presumably via geometry-dependent distribution of tension. In this study, we identify common shape determinants spanning cell and tissue scales. For cells whose sites of adhesion are restricted to small adhesive islands on a micropatterned substrate, shape resembles a sequence of inward-curved circular arcs. The same shape is observed for fibroblast-populated collagen gels that are pinned to a flat substrate. Quantitative image analysis reveals that, in both cases, arc radii increase with the spanning distance between the pinning points. Although the Laplace law for interfaces under tension predicts circular arcs, it cannot explain the observed dependence on the spanning distance. Computer simulations and theoretical modeling demonstrate that filamentous network mechanics and contractility give rise to a modified Laplace law that quantitatively explains our experimental findings on both cell and tissue scales. Our model in conjunction with actomyosin inhibition experiments further suggests that cell shape is regulated by two different control modes related to motor contractility and structural changes in the actin cytoskeleton.

Collective effects in cellular structure formation mediated by compliant environments: a Monte Carlo study.

Compliant environments can mediate interactions between mechanically active cells like fibroblasts. Starting with a phenomenological model for the behavior of single cells, we use extensive Monte Carlo simulations to predict non-trivial structure formation for cell communities on soft elastic substrates as a function of elastic moduli, cell density, noise and cell position geometry. In general, we find a disordered structure as well as ordered string-like and ring-like structures. The transition between ordered and disordered structures is controlled both by cell density and noise level, while the transition between string- and ring-like ordered structures is controlled by the Poisson ratio. Similar effects are observed in three dimensions. Our results suggest that in regard to elastic effects, healthy connective tissue usually is in a macroscopically disordered state, but can be switched to a macroscopically ordered state by appropriate parameter variations, in a way that is reminiscent of wound contraction or diseased states like contracture.

Focal adhesions as mechanosensors: the two-spring model.

Adhesion-dependent cells actively sense the mechanical properties of their environment through mechanotransductory processes at focal adhesions, which are integrin-based contacts connecting the extracellular matrix to the cytoskeleton. Here we present first steps towards a quantitative understanding of focal adhesions as mechanosensors. It has been shown experimentally that high levels of force are related to growth of and signaling at focal adhesions. In particular, activation of the small GTPase Rho through focal adhesions leads to the formation of stress fibers. Here we discuss one way in which force might regulate the internal state of focal adhesions, namely by modulating the internal rupture dynamics of focal adhesions. A simple two-spring model shows that the stiffer the environment, the more efficient cellular force is built up at focal adhesions by molecular motors interacting with the actin filaments.

Physical determinants of cell organization in soft media.

Cell adhesion is an integral part of many physiological processes in tissues, including development, tissue maintenance, angiogenesis, and wound healing. Recent advances in materials science (including microcontact printing, soft lithography, microfluidics, and nanotechnology) have led to strongly improved control of extracellular ligand distribution and of the properties of the micromechanical environment. As a result, the investigation of cellular response to the physical properties of adhesive surfaces has become a very active area of research. Sophisticated use of elastic substrates has revealed that cell organization in soft media is determined by active mechanosensing at cell-matrix adhesions. In order to determine the underlying mechanisms, quantification and biophysical modelling are essential. In tissue engineering, theory might help to design new environments for cells.