A programmable reaction-diffusion system for spatiotemporal cell signaling circuit design.

Cells self-organize molecules in space and time to generate complex behaviors, but we lack synthetic strategies for engineering spatiotemporal signaling. We present a programmable reaction-diffusion platform for designing protein oscillations, patterns, and circuits in mammalian cells using two bacterial proteins, MinD and MinE (MinDE). MinDE circuits act like "single-cell radios," emitting frequency-barcoded fluorescence signals that can be spectrally isolated and analyzed using digital signal processing tools. We define how to genetically program these signals and connect their spatiotemporal dynamics to cell biology using engineerable protein-protein interactions. This enabled us to construct sensitive reporter circuits that broadcast endogenous cell signaling dynamics on a frequency-barcoded imaging channel and to build control signal circuits that synthetically pattern activities in the cell, such as protein condensate assembly and actin filamentation. Our work establishes a paradigm for visualizing, probing, and engineering cellular activities at length and timescales critical for biological function.

Engineering Customized Cell Sensing and Response Behaviors Using Synthetic Notch Receptors.

The Notch protein is one of the most mechanistically direct transmembrane receptors-the intracellular domain contains a transcriptional regulator that is released from the membrane when engagement of the cognate extracellular ligand induces intramembrane proteolysis. We find that chimeric forms of Notch, in which both the extracellular sensor module and the intracellular transcriptional module are replaced with heterologous protein domains, can serve as a general platform for generating novel cell-cell contact signaling pathways. Synthetic Notch (synNotch) pathways can drive user-defined functional responses in diverse mammalian cell types. Because individual synNotch pathways do not share common signaling intermediates, the pathways are functionally orthogonal. Thus, multiple synNotch receptors can be used in the same cell to achieve combinatorial integration of environmental cues, including Boolean response programs, multi-cellular signaling cascades, and self-organized cellular patterns. SynNotch receptors provide extraordinary flexibility in engineering cells with customized sensing/response behaviors to user-specified extracellular cues.

Product binding enforces the genomic specificity of a yeast polycomb repressive complex.

We characterize the Polycomb system that assembles repressive subtelomeric domains of H3K27 methylation (H3K27me) in the yeast Cryptococcus neoformans. Purification of this PRC2-like protein complex reveals orthologs of animal PRC2 components as well as a chromodomain-containing subunit, Ccc1, which recognizes H3K27me. Whereas removal of either the EZH or EED ortholog eliminates H3K27me, disruption of mark recognition by Ccc1 causes H3K27me to redistribute. Strikingly, the resulting pattern of H3K27me coincides with domains of heterochromatin marked by H3K9me. Indeed, additional removal of the C. neoformans H3K9 methyltransferase Clr4 results in loss of both H3K9me and the redistributed H3K27me marks. These findings indicate that the anchoring of a chromatin-modifying complex to its product suppresses its attraction to a different chromatin type, explaining how enzymes that act on histones, which often harbor product recognition modules, may deposit distinct chromatin domains despite sharing a highly abundant and largely identical substrate-the nucleosome.

Exploitation of latent allostery enables the evolution of new modes of MAP kinase regulation.

Allosteric interactions provide precise spatiotemporal control over signaling proteins, but how allosteric activators and their targets coevolve is poorly understood. Here, we trace the evolution of two allosteric activator motifs within the yeast scaffold protein Ste5 that specifically target the mating MAP kinase Fus3. One activator (Ste5-VWA) provides pathway insulation and dates to the divergence of Fus3 from its paralog, Kss1; a second activator (Ste5-FBD) that tunes mating behavior is, in contrast, not conserved in most lineages. Surprisingly, both Ste5 activator motifs could regulate MAP kinases that diverged from Fus3 prior to the emergence of Ste5, suggesting that Ste5 activators arose by exploiting latent regulatory features already present in the MAPK ancestor. The magnitude of this latent allosteric potential drifts widely among pre-Ste5 MAP kinases, providing a pool of hidden phenotypic diversity that, when revealed by new activators, could lead to functional divergence and to the evolution of distinct signaling behaviors.

Coupled 5' nucleotide recognition and processivity in Xrn1-mediated mRNA decay.

Messenger RNA decay plays a central role in the regulation and surveillance of eukaryotic gene expression. The conserved multidomain exoribonuclease Xrn1 targets cytoplasmic RNA substrates marked by a 5' monophosphate for processive 5'-to-3' degradation by an unknown mechanism. Here, we report the crystal structure of an Xrn1-substrate complex. The single-stranded substrate is held in place by stacking of the 5'-terminal trinucleotide between aromatic side chains while a highly basic pocket specifically recognizes the 5' phosphate. Mutations of residues involved in binding the 5'-terminal nucleotide impair Xrn1 processivity. The substrate recognition mechanism allows Xrn1 to couple processive hydrolysis to duplex melting in RNA substrates with sufficiently long single-stranded 5' overhangs. The Xrn1-substrate complex structure thus rationalizes the exclusive specificity of Xrn1 for 5'-monophosphorylated substrates, ensuring fidelity of mRNA turnover, and posits a model for translocation-coupled unwinding of structured RNA substrates.

Ers1 links HP1 to RNAi.

Pericentromeric heterochromatin formation is mediated by repressive histone H3 lysine 9 methylation (H3K9Me) and its recognition by HP1 proteins. Intriguingly, in many organisms, RNAi is coupled to this process through poorly understood mechanisms. In Schizosaccharomyces pombe, the H3-K9 methyltransferase Clr4 and the heterochromatin protein 1 (HP1) ortholog Swi6 are critical for RNAi, whereas RNAi stimulates H3K9Me. In addition to the endoribonuclease Dcr1, RNAi in S. pombe requires two interacting protein complexes, the RITS complex, which contains an Argonaute subunit, and the RDRC complex, which contains an RNA-dependent RNA polymerase subunit. We previously identified Ers1 (essential for RNAi-dependent silencing) as an orphan protein that genetically acts in the RNAi pathway. Using recombinant proteins, we show here that Ers1 directly and specifically interacts with HP1/Swi6. Two-hybrid assays indicate that Ers1 also directly interacts with several RNAi factors. Consistent with these interactions, Ers1 associates in vivo with the RITS complex, the RDRC complex, and Dcr1, and it promotes interactions between these factors. Ers1, like Swi6, is also required for RNAi complexes to associate with pericentromeric noncoding RNAs. Overexpression of Ers1 results in a dominant-negative phenotype that can be specifically suppressed by increasing levels of the RDRC subunit Hrr1 or of Dcr1, further supporting a functional role for Ers1 in promoting the assembly of the RNAi machinery. Through the interactions described here, Ers1 may promote RNAi by tethering the corresponding enzyme complexes to HP1-coated chromatin, thereby placing them in proximity to the nascent noncoding RNA substrate.

Comparative profiling of cellular gait on adhesive micropatterns defines statistical patterns of activity that underlie native and cancerous cell dynamics.

Cell dynamics are powered by patterns of activity, but it is not straightforward to quantify these patterns or compare them across different environmental conditions or cell-types. Here we digitize the long-term shape fluctuations of metazoan cells grown on micropatterned fibronectin islands to define and extract statistical features of cell dynamics without the need for genetic modification or fluorescence imaging. These shape fluctuations generate single-cell morphological signals that can be decomposed into two major components: a continuous, slow-timescale meandering of morphology about an average steady-state shape; and short-lived "events" of rapid morphology change that sporadically occur throughout the timecourse. By developing statistical metrics for each of these components, we used thousands of hours of single-cell data to quantitatively define how each axis of cell dynamics was impacted by environmental conditions or cell-type. We found the size and spatial complexity of the micropattern island modulated the statistics of morphological events-lifetime, frequency, and orientation-but not its baseline shape fluctuations. Extending this approach to profile a panel of triple negative breast cancer cell-lines, we found that different cell-types could be distinguished from one another along specific and unique statistical axes of their behavior. Our results suggest that micropatterned substrates provide a generalizable method to build statistical profiles of cell dynamics to classify and compare emergent cell behaviors.

Regulatable assembly of synthetic microtubule architectures using engineered MAP-IDR condensates.

Microtubules filaments are assembled into higher-order structures and machines critical for cellular processes using microtubule-associated proteins (MAPs). However, the design of synthetic MAPs that direct the formation of new structures in cells is challenging, as nanoscale biochemical activities must be organized across micron length-scales. Here we develop synthetic MAP-IDR condensates (synMAPs) that provide tunable and regulatable assembly of higher-order microtubule structures in vitro and in mammalian cells. synMAPs harness a small microtubule-binding domain from oligodendrocytes (TPPP) whose activity can be synthetically rewired by interaction with condensate-forming IDR sequences. This combination allows synMAPs to self-organize multivalent structures that bind and bridge microtubules into synthetic architectures. Regulating the connection between the microtubule-binding and condensate-forming components allows synMAPs to act as nodes in more complex cytoskeletal circuits in which the formation and dynamics of the microtubule structure can be controlled by small molecules or cell-signaling inputs. By systematically testing a panel of synMAP circuit designs, we define a two-level control scheme for dynamic assembly of microtubule architectures at the nanoscale (via microtubule-binding) and microscale (via condensate formation). synMAPs provide a compact and rationally engineerable starting point for the design of more complex microtubule architectures and cellular machines.

Structural and biochemical studies of a fluoroacetyl-CoA-specific thioesterase reveal a molecular basis for fluorine selectivity.

We have initiated a broad-based program aimed at understanding the molecular basis of fluorine specificity in enzymatic systems, and in this context, we report crystallographic and biochemical studies on a fluoroacetyl-coenzyme A (CoA) specific thioesterase (FlK) from Streptomyces cattleya. Our data establish that FlK is competent to protect its host from fluoroacetate toxicity in vivo and demonstrate a 10(6)-fold discrimination between fluoroacetyl-CoA (k(cat)/K(M) = 5 × 107 M⁻¹ s⁻¹) and acetyl-CoA (k(cat)/K(M) = 30 M⁻¹ s⁻¹) based on a single fluorine substitution that originates from differences in both substrate reactivity and binding. We show that Thr 42, Glu 50, and His 76 are key catalytic residues and identify several factors that influence substrate selectivity. We propose that FlK minimizes interaction with the thioester carbonyl, leading to selection against acetyl-CoA binding that can be recovered in part by new C═O interactions in the T42S and T42C mutants. We hypothesize that the loss of these interactions is compensated by the entropic driving force for fluorinated substrate binding in a hydrophobic binding pocket created by a lid structure, containing Val 23, Leu 26, Phe 33, and Phe 36, that is not found in other structurally characterized members of this superfamily. We further suggest that water plays a critical role in fluorine specificity based on biochemical and structural studies focused on the unique Phe 36 "gate" residue, which functions to exclude water from the active site. Taken together, the findings from these studies offer molecular insights into organofluorine recognition and design of fluorine-specific enzymes.

Ciliate behavior: blueprints for dynamic cell biology and microscale robotics.

Place a drop of pond water under the microscope, and you will likely find an ocean of extraordinary and diverse single-celled organisms called ciliates. This remarkable group of single-celled organisms wield microtubules, active systems, electrical signaling, and chemical sensors to build intricate geometrical structures and perform complex behaviors that can appear indistinguishable from those of macroscopic animals. Advances in computer vision and machine learning are making it possible to completely digitize and track the dynamics of complex ciliates and mine these data for the hidden structure, patterns, and motifs that are responsible for their behaviors. By deconstructing the diversity of ciliate behaviors in the natural world, themes for organizing and controlling matter at the microscale are beginning to take hold, suggesting new modular approaches for the design of autonomous molecular machines that emulate nature's finest examples.

Starvation induces genomic rearrangements and starvation-resilient phenotypes in yeast.

Evolution has shaped a wide variety of genomes across eukaryotic taxa. However, the forces that shape the genomes are generally unknown. Because organisms in nature commonly experience prolonged periods of nutrient depletion, we posit that diverse demographic, physiological, and genomic responses to starvation can occur. To test for these possibilities, we subjected replicate yeast populations to prolonged starvation. We observed that clones repeatedly gave rise to descendants that were karyotypically diverse. After a 1-month starvation period, approximately 70% of randomly isolated members of starved populations harbored one or more genomic rearrangements. Further, we found that 5 of 16 karyotypically differentiated groups of isolates from starved populations were more resilient to starvation than nonstarved clones and their common ancestor. Phylogenetic analysis of these isolates suggests that genomic rearrangements that arose during starvation can be adaptive in the context of a nutrient-depleted environment. Altogether our data illustrate the profound influence of environmental conditions on adaptive genome evolution in eukaryotes.

Coupled Active Systems Encode an Emergent Hunting Behavior in the Unicellular Predator Lacrymaria olor.

Many single-celled protists use rapid morphology changes to perform fast animal-like behaviors. To understand how such behaviors are encoded, we analyzed the hunting dynamics of the predatory ciliate Lacrymaria olor, which locates and captures prey using the tip of a slender "neck" that can rapidly extend more than seven times its body length (500 µm from its body) and retract in seconds. By tracking single cells in real-time over hours and analyzing millions of sub-cellular postures, we find that these fast extension-contraction cycles underlie an emergent hunting behavior that comprehensively samples a broad area within the cell's reach. Although this behavior appears complex, we show that it arises naturally as alternating sub-cellular ciliary and contractile activities rearrange the cell's underlying helical cytoskeleton to extend or retract the neck. At short timescales, a retracting neck behaves like an elastic filament under load, such that compression activates a series of buckling modes that reorient the head and scramble its extensile trajectory. At longer timescales, the fundamental length of this filament can change, altering the location in space where these transitions occur. Coupling these fast and slow dynamics together, we present a simple model for how Lacrymaria samples the range of geometries and orientations needed to ensure dense stochastic sampling of the immediate environment when hunting to locate and strike at prey. More generally, coupling active mechanical and chemical signaling systems across different timescales may provide a general strategy by which mechanically encoded emergent cell behaviors can be understood or engineered.

Reverse engineering GTPase programming languages with reconstituted signaling networks.

The Ras superfamily GTPases represent one of the most prolific signaling currencies used in Eukaryotes. With these remarkable molecules, evolution has built GTPase networks that control diverse cellular processes such as growth, morphology, motility and trafficking. (1-4) Our knowledge of the individual players that underlie the function of these networks is deep; decades of biochemical and structural data has provided a mechanistic understanding of the molecules that turn GTPases ON and OFF, as well as how those GTPase states signal by controlling the assembly of downstream effectors. However, we know less about how these different activities work together as a system to specify complex dynamic signaling outcomes. Decoding this molecular "programming language" would help us understand how different species and cell types have used the same GTPase machinery in different ways to accomplish different tasks, and would also provide new insights as to how mutations to these networks can cause disease. We recently developed a bead-based microscopy assay to watch reconstituted H-Ras signaling systems at work under arbitrary configurations of regulators and effectors. (5) Here we highlight key observations and insights from this study and propose extensions to our method to further study this and other GTPase signaling systems.

Mapping the functional versatility and fragility of Ras GTPase signaling circuits through in vitro network reconstitution.

The Ras-superfamily GTPases are central controllers of cell proliferation and morphology. Ras signaling is mediated by a system of interacting molecules: upstream enzymes (GEF/GAP) regulate Ras's ability to recruit multiple competing downstream effectors. We developed a multiplexed, multi-turnover assay for measuring the dynamic signaling behavior of in vitro reconstituted H-Ras signaling systems. By including both upstream regulators and downstream effectors, we can systematically map how different network configurations shape the dynamic system response. The concentration and identity of both upstream and downstream signaling components strongly impacted the timing, duration, shape, and amplitude of effector outputs. The distorted output of oncogenic alleles of Ras was highly dependent on the balance of positive (GAP) and negative (GEF) regulators in the system. We found that different effectors interpreted the same inputs with distinct output dynamics, enabling a Ras system to encode multiple unique temporal outputs in response to a single input. We also found that different Ras-to-GEF positive feedback mechanisms could reshape output dynamics in distinct ways, such as signal amplification or overshoot minimization. Mapping of the space of output behaviors accessible to Ras provides a design manual for programming Ras circuits, and reveals how these systems are readily adapted to produce an array of dynamic signaling behaviors. Nonetheless, this versatility comes with a trade-off of fragility, as there exist numerous paths to altered signaling behaviors that could cause disease.

Specificity in endoplasmic reticulum-stress signaling in yeast entails a step-wise engagement of HAC1 mRNA to clusters of the stress sensor Ire1.

Insufficient protein-folding capacity in the endoplasmic reticulum (ER) induces the unfolded protein response (UPR). In the ER lumen, accumulation of unfolded proteins activates the transmembrane ER-stress sensor Ire1 and drives its oligomerization. In the cytosol, Ire1 recruits HAC1 mRNA, mediating its non-conventional splicing. The spliced mRNA is translated into Hac1, the key transcription activator of UPR target genes that mitigate ER-stress. In this study, we report that oligomeric assembly of the ER-lumenal domain is sufficient to drive Ire1 clustering. Clustering facilitates Ire1's cytosolic oligomeric assembly and HAC1 mRNA docking onto a positively charged motif in Ire1's cytosolic linker domain that tethers the kinase/RNase to the transmembrane domain. By the use of a synthetic bypass, we demonstrate that mRNA docking per se is a pre-requisite for initiating Ire1's RNase activity and, hence, splicing. We posit that such step-wise engagement between Ire1 and its mRNA substrate contributes to selectivity and efficiency in UPR signaling.

Starvation-associated genome restructuring can lead to reproductive isolation in yeast.

Knowledge of the mechanisms that lead to reproductive isolation is essential for understanding population structure and speciation. While several models have been advanced to explain post-mating reproductive isolation, experimental data supporting most are indirect. Laboratory investigations of this phenomenon are typically carried out under benign conditions, which result in low rates of genetic change unlikely to initiate reproductive isolation. Previously, we described an experimental system using the yeast Saccharomyces cerevisiae where starvation served as a proxy to any stress that decreases reproduction and/or survivorship. We showed that novel lineages with restructured genomes quickly emerged in starved populations, and that these survivors were more fit than their ancestors when re-starved. Here we show that certain yeast lineages that survive starvation have become reproductively isolated from their ancestor. We further demonstrate that reproductive isolation arises from genomic rearrangements, whose frequency in starving yeast is several orders of magnitude greater than an unstarved control. By contrast, the frequency of point mutations is less than 2-fold greater. In a particular case, we observe that a starved lineage becomes reproductively isolated as a direct result of the stress-related accumulation of a single chromosome. We recapitulate this result by demonstrating that introducing an extra copy of one or several chromosomes into naïve, i.e. unstarved, yeast significantly diminishes their fertility. This type of reproductive barrier, whether arising spontaneously or via genetic manipulation, can be removed by making a lineage euploid for the altered chromosomes. Our model provides direct genetic evidence that reproductive isolation can arise frequently in stressed populations via genome restructuring without the precondition of geographic isolation.

Conformational control of the Ste5 scaffold protein insulates against MAP kinase misactivation.

Cells reuse signaling proteins in multiple pathways, raising the potential for improper cross talk. Scaffold proteins are thought to insulate against such miscommunication by sequestering proteins into distinct physical complexes. We show that the scaffold protein Ste5, which organizes the yeast mating mitogen-activated protein kinase (MAPK) pathway, does not use sequestration to prevent misactivation of the mating response. Instead, Ste5 appears to use a conformation mechanism: Under basal conditions, an intramolecular interaction of the pleckstrin homology (PH) domain with the von Willebrand type A (VWA) domain blocks the ability to coactivate the mating-specific MAPK Fus3. Pheromone-induced membrane binding of Ste5 triggers release of this autoinhibition. Thus, in addition to serving as a conduit guiding kinase communication, Ste5 directly receives input information to decide if and when signal can be transmitted to mating output.

Structural insights into the human GW182-PABC interaction in microRNA-mediated deadenylation.

GW182-family proteins are essential for microRNA-mediated translational repression and deadenylation in animal cells. Here we show that a conserved motif in the human GW182 paralog TNRC6C interacts with the C-terminal domain of polyadenylate binding protein 1 (PABC) and present the crystal structure of the complex. Mutations at the complex interface impair mRNA deadenylation in mammalian cell extracts, suggesting that the GW182-PABC interaction contributes to microRNA-mediated gene silencing.

Direct link between RACK1 function and localization at the ribosome in vivo.

The receptor for activated C-kinase (RACK1), a conserved protein implicated in numerous signaling pathways, is a stoichiometric component of eukaryotic ribosomes located on the head of the 40S ribosomal subunit. To test the hypothesis that ribosome association is central to the function of RACK1 in vivo, we determined the 2.1-A crystal structure of RACK1 from Saccharomyces cerevisiae (Asc1p) and used it to design eight mutant versions of RACK1 to assess roles in ribosome binding and in vivo function. Conserved charged amino acids on one side of the beta-propeller structure were found to confer most of the 40S subunit binding affinity, whereas an adjacent conserved and structured loop had little effect on RACK1-ribosome association. Yeast mutations that confer moderate to strong defects in ribosome binding mimic some phenotypes of a RACK1 deletion strain, including increased sensitivity to drugs affecting cell wall biosynthesis and translation elongation. Furthermore, disruption of RACK1's position at the 40S ribosomal subunit results in the failure of the mRNA binding protein Scp160 to associate with actively translating ribosomes. These results provide the first direct evidence that RACK1 functions from the ribosome, implying a physical link between the eukaryotic ribosome and cell signaling pathways in vivo.

Structural basis for DNase activity of a conserved protein implicated in CRISPR-mediated genome defense.

Acquired immunity in prokaryotes is achieved by integrating short fragments of foreign nucleic acids into clustered regularly interspaced short palindromic repeats (CRISPRs). This nucleic acid-based immune system is mediated by a variable cassette of up to 45 protein families that represent distinct immune system subtypes. CRISPR-associated gene 1 (cas1) encodes the only universally conserved protein component of CRISPR immune systems, yet its function is unknown. Here we show that the Cas1 protein is a metal-dependent DNA-specific endonuclease that produces double-stranded DNA fragments of approximately 80 base pairs in length. The 2.2 A crystal structure of the Cas1 protein reveals a distinct fold and a conserved divalent metal ion-binding site. Mutation of metal ion-binding residues, chelation of metal ions, or metal-ion substitution inhibits Cas1-catalyzed DNA degradation. These results provide a foundation for understanding how Cas1 contributes to CRISPR function, perhaps as part of the machinery for processing foreign nucleic acids.