The transcription of genomic information in eukaryotes is regulated in large part by chromatin. How a diverse array of chromatin regulator (CR) proteins with different functions and genomic localization patterns coordinates chromatin activity to control transcription remains unclear. Here, we take a synthetic biology approach to decipher the complexity of chromatin regulation by studying emergent transcriptional behaviors from engineered combinatorial, spatial, and temporal patterns of individual CRs. We fuse 223 yeast CRs to programmable zinc finger proteins. Site-specific and combinatorial recruitment of CRs to distinct intralocus locations reveals a range of transcriptional logic and behaviors, including synergistic activation, long-range and spatial regulation, and gene expression memory. Comparing these transcriptional behaviors with annotated CR complex and function terms provides design principles for the engineering of transcriptional regulation. This work presents a bottom-up approach to investigating chromatin-mediated transcriptional regulation and introduces chromatin-based components and systems for synthetic biology and cellular engineering.

It is generally held that random-coil polypeptide chains undergo a barrier-less continuous collapse when the solvent conditions are changed to favor the fully folded native conformation. We test this hypothesis by probing intramolecular distance distributions during folding in one of the paradigms of folding reactions, that of cytochrome c. The Trp59-to-heme distance was probed by time-resolved Förster resonance energy transfer in the microsecond time range of refolding. Contrary to expectation, a state with a Trp59-heme distance close to that of the guanidinium hydrochloride (GdnHCl) denatured state is present after ~27 μs of folding. A concomitant decrease in the population of this state and an increase in the population of a compact high-FRET (Förster resonance energy transfer) state (efficiency>90%) show that the collapse is barrier limited. Small-angle X-ray scattering (SAXS) measurements over a similar time range show that the radius of gyration under native favoring conditions is comparable to that of the GdnHCl denatured unfolded state. An independent comprehensive global thermodynamic analysis reveals that marginally stable partially folded structures are also present in the nominally unfolded GdnHCl denatured state. These observations suggest that specifically collapsed intermediate structures with low stability in rapid equilibrium with the unfolded state may contribute to the apparent chain contraction observed in previous fluorescence studies using steady-state detection. In the absence of significant dynamic averaging of marginally stable partially folded states and with the use of probes sensitive to distance distributions, barrier-limited chain contraction is observed upon transfer of the GdnHCl denatured state ensemble to native-like conditions.

Chemical modifications to DNA and histone proteins are involved in epigenetic programs underlying cellular differentiation and development. Regulatory networks involving molecular writers and readers of chromatin marks are thought to control these programs. Guided by this common principle, we established an orthogonal epigenetic regulatory system in mammalian cells using N6-methyladenine (m6A), a DNA modification not commonly found in metazoan epigenomes. Our system utilizes synthetic factors that write and read m6A and consequently recruit transcriptional regulators to control reporter loci. Inspired by models of chromatin spreading and epigenetic inheritance, we used our system and mathematical models to construct regulatory circuits that induce m6A-dependent transcriptional states, promote their spatial propagation, and maintain epigenetic memory of the states. These minimal circuits were able to program epigenetic functions de novo, conceptually validating "read-write" architectures. This work provides a toolkit for investigating models of epigenetic regulation and encoding additional layers of epigenetic information in cells.

Genetic, epidemiological and experimental evidence implicate lysosomal dysfunction in Parkinson's disease (PD) and related synucleinopathies. Investigate several mouse models of lysosomal storage diseases (LSDs) and evaluate pathologies reminiscent of synucleinopathies. We obtained brain tissue from symptomatic mouse models of Gaucher, Fabry, Sandhoff, Niemann-Pick A (NPA), Hurler, Pompe and Niemann-Pick C (NPC) diseases and assessed for the presence of Lewy body-like pathology (proteinase K-resistant α-synuclein and tau aggregates) and neuroinflammation (microglial Iba1 and astrocytic GFAP) by immunofluorescence. All seven LSD models exhibited evidence of proteinopathy and/or inflammation in the central nervous system (CNS). However, these phenotypes were divergent. Gaucher and Fabry mouse models displayed proteinase K-resistant α-synuclein and tau aggregates but no neuroinflammation; whereas Sandhoff, NPA and NPC showed marked neuroinflammation and no overt proteinopathy. Pompe disease animals uniquely displayed widespread distribution of tau aggregates accompanied by moderate microglial activation. Hurler mice also demonstrated proteinopathy and microglial activation. The present study demonstrated additional links between LSDs and pathogenic phenotypes that are hallmarks of synucleinopathies. The data suggest that lysosomal dysregulation can contribute to brain region-specific protein aggregation and induce widespread neuroinflammation in the brain. However, only a few LSD models examined exhibited phenotypes consistent with synucleinopathies. While no model can recapitulate the complexity of PD, they can enable the study of specific pathways and mechanisms contributing to disease pathophysiology. The present study provides evidence that there are existing, previously unutilized mouse models that can be employed to study pathogenic mechanisms and gain insights into potential PD subtypes, helping to determine if they are amenable to pathway-specific therapeutic interventions.

The gut microbiome is essential for processing complex food compounds and synthesizing nutrients that the host cannot digest or produce, respectively. New model systems are needed to study how the metabolic capacity provided by the gut microbiome impacts the nutritional status of the host, and to explore possibilities for altering host metabolic capacity via the microbiome. Here, we colonized the nematode Caenorhabditis elegans gut with cellulolytic bacteria that enabled C. elegans to utilize cellulose, an otherwise indigestible substrate, as a carbon source. Cellulolytic bacteria as a community component in the worm gut can also support additional bacterial species with specialized roles, which we demonstrate by using Lactobacillus plantarum to protect C. elegans against Salmonella enterica infection. This work shows that engineered microbiome communities can be used to endow host organisms with novel functions, such as the ability to utilize alternate nutrient sources or to better fight pathogenic bacteria.

Despite the availability of Cas9 variants with varied protospacer-adjacent motif (PAM) compatibilities, some genomic loci-especially those with pyrimidine-rich PAM sequences-remain inaccessible by high-activity Cas9 proteins. Moreover, broadening PAM sequence compatibility through engineering can increase off-target activity. With directed evolution, we generated four Cas9 variants that together enable targeting of most pyrimidine-rich PAM sequences in the human genome. Using phage-assisted noncontinuous evolution and eVOLVER-supported phage-assisted continuous evolution, we evolved Nme2Cas9, a compact Cas9 variant, into variants that recognize single-nucleotide pyrimidine-PAM sequences. We developed a general selection strategy that requires functional editing with fully specified target protospacers and PAMs. We applied this selection to evolve high-activity variants eNme2-T.1, eNme2-T.2, eNme2-C and eNme2-C.NR. Variants eNme2-T.1 and eNme2-T.2 offer access to N4TN PAM sequences with comparable editing efficiencies as existing variants, while eNme2-C and eNme2-C.NR offer less restrictive PAM requirements, comparable or higher activity in a variety of human cell types and lower off-target activity at N4CN PAM sequences.

Synthetic biology toolkits are one of the core foundations on which the field has been built, facilitating and accelerating efforts to reprogram cells and organisms for diverse biotechnological applications. The yeast Saccharomyces cerevisiae, an important model and industrial organism, has benefited from a wide range of toolkits. In particular, the MoClo Yeast Toolkit (YTK) enables the fast and straightforward construction of multigene plasmids from a library of highly characterized parts for programming new cellular behavior in a more predictable manner. While YTK has cultivated a strong parts ecosystem and excels in plasmid construction, it is limited in the extent and flexibility with which it can create new strains of yeast. Here, we describe a new and improved toolkit, the Multiplex Yeast Toolkit (MYT), that extends the capabilities of YTK and addresses strain engineering limitations. MYT provides a set of new integration vectors and selectable markers usable across common laboratory strains, as well as additional assembly cassettes to increase the number of transcriptional units in multigene constructs, CRISPR-Cas9 tools for highly efficient multiplexed vector integration, and three orthogonal and inducible promoter systems for conditional programming of gene expression. With these tools, we provide yeast synthetic biologists with a powerful platform to take their engineering ambitions to exciting new levels.

A ubiquitous feature of eukaryotic transcriptional regulation is cooperative self-assembly between transcription factors (TFs) and DNA cis-regulatory motifs. It is thought that this strategy enables specific regulatory connections to be formed in gene networks between otherwise weakly interacting, low-specificity molecular components. Here, using synthetic gene circuits constructed in yeast, we find that high regulatory specificity can emerge from cooperative, multivalent interactions among artificial zinc-finger-based TFs. We show that circuits "wired" using the strategy of cooperative TF assembly are effectively insulated from aberrant misregulation of the host cell genome. As we demonstrate in experiments and mathematical models, this mechanism is sufficient to rescue circuit-driven fitness defects, resulting in genetic and functional stability of circuits in long-term continuous culture. Our naturally inspired approach offers a simple, generalizable means for building high-fidelity, evolutionarily robust gene circuits that can be scaled to a wide range of host organisms and applications.

CRISPR-based gene drives have the potential to spread within populations and are considered as promising vector control tools. A doublesex-targeting gene drive was able to suppress laboratory Anopheles mosquito populations in small and large cages, and it is considered for field application. Challenges related to the field-use of gene drives and the evolving regulatory framework suggest that systems able to modulate or revert the action of gene drives, could be part of post-release risk-mitigation plans. In this study, we challenge an AcrIIA4-based anti-drive to inhibit gene drive spread in age-structured Anopheles gambiae population under complex feeding and behavioural conditions. A stochastic model predicts the experimentally-observed genotype dynamics in age-structured populations in medium-sized cages and highlights the necessity of large-sized cage trials. These experiments and experimental-modelling framework demonstrate the effectiveness of the anti-drive in different scenarios, providing further corroboration for its use in controlling the spread of gene drive in Anopheles.

Alpha-synuclein (αSyn) protein levels correlate with the risk and severity of Parkinson's disease and related neurodegenerative diseases. Lowering αSyn is being actively investigated as a therapeutic modality. Here, we systematically map the regulatory network that controls endogenous αSyn using sequential CRISPR-knockout and -interference screens in an αSyn gene (SNCA)-tagged cell line and induced pluripotent stem cell-derived neurons (iNeurons). We uncover αSyn modifiers at multiple regulatory layers, with amino-terminal acetyltransferase B (NatB) enzymes being the most potent endogenous αSyn modifiers in both cell lines. Amino-terminal acetylation protects the cytosolic αSyn from rapid degradation by the proteasome in a Ube2w-dependent manner. Moreover, we show that pharmacological inhibition of methionyl-aminopeptidase 2, a regulator of NatB complex formation, attenuates endogenous αSyn in iNeurons carrying SNCA triplication. Together, our study reveals several gene networks that control endogenous αSyn, identifies mechanisms mediating the degradation of nonacetylated αSyn, and illustrates potential therapeutic pathways for decreasing αSyn levels in synucleinopathies.

Precise control over microbial cell growth conditions could enable detection of minute phenotypic changes, which would improve our understanding of how genotypes are shaped by adaptive selection. Although automated cell-culture systems such as bioreactors offer strict control over liquid culture conditions, they often do not scale to high-throughput or require cumbersome redesign to alter growth conditions. We report the design and validation of eVOLVER, a scalable do-it-yourself (DIY) framework, which can be configured to carry out high-throughput growth experiments in molecular evolution, systems biology, and microbiology. High-throughput evolution of yeast populations grown at different densities reveals that eVOLVER can be applied to characterize adaptive niches. Growth selection on a genome-wide yeast knockout library, using temperatures varied over different timescales, finds strains sensitive to temperature changes or frequency of temperature change. Inspired by large-scale integration of electronics and microfluidics, we also demonstrate millifluidic multiplexing modules that enable multiplexed media routing, cleaning, vial-to-vial transfers and automated yeast mating.

Clonal populations of cells exhibit cell-to-cell variation in the transcription of individual genes. In addition to this noise in gene expression, heterogeneity in the proteome and the proteostasis network expands the phenotypic diversity of a population. Heat shock factor 1 (Hsf1) regulates chaperone gene expression, thereby coupling transcriptional noise to proteostasis. Here we show that cell-to-cell variation in Hsf1 activity is an important determinant of phenotypic plasticity. Budding yeast cells with high Hsf1 activity were enriched for the ability to acquire resistance to an antifungal drug, and this enrichment depended on Hsp90, a known phenotypic capacitor and canonical Hsf1 target. We show that Hsf1 phosphorylation promotes cell-to-cell variation, and this variation, rather than absolute Hsf1 activity, promotes antifungal resistance. We propose that Hsf1 phosphorylation enables differential tuning of the proteostasis network in individual cells, allowing populations to access a range of phenotypic states.

Models for regulation of the eukaryotic heat shock response typically invoke a negative feedback loop consisting of the transcriptional activator Hsf1 and a molecular chaperone. Previously we identified Hsp70 as the chaperone responsible for Hsf1 repression and constructed a mathematical model that recapitulated the yeast heat shock response (Zheng et al., 2016). The model was based on two assumptions: dissociation of Hsp70 activates Hsf1, and transcriptional induction of Hsp70 deactivates Hsf1. Here we validate these assumptions. First, we severed the feedback loop by uncoupling Hsp70 expression from Hsf1 regulation. As predicted by the model, Hsf1 was unable to efficiently deactivate in the absence of Hsp70 transcriptional induction. Next, we mapped a discrete Hsp70 binding site on Hsf1 to a C-terminal segment known as conserved element 2 (CE2). In vitro, CE2 binds to Hsp70 with low affinity (9 µM), in agreement with model requirements. In cells, removal of CE2 resulted in increased basal Hsf1 activity and delayed deactivation during heat shock, while tandem repeats of CE2 sped up Hsf1 deactivation. Finally, we uncovered a role for the N-terminal domain of Hsf1 in negatively regulating DNA binding. These results reveal the quantitative control mechanisms underlying the heat shock response.

Alpha-synuclein (αS) is a conformationally plastic protein that reversibly binds to cellular membranes. It aggregates and is genetically linked to Parkinson's disease (PD). Here, we show that αS directly modulates processing bodies (P-bodies), membraneless organelles that function in mRNA turnover and storage. The N terminus of αS, but not other synucleins, dictates mutually exclusive binding either to cellular membranes or to P-bodies in the cytosol. αS associates with multiple decapping proteins in close proximity on the Edc4 scaffold. As αS pathologically accumulates, aberrant interaction with Edc4 occurs at the expense of physiologic decapping-module interactions. mRNA decay kinetics within PD-relevant pathways are correspondingly disrupted in PD patient neurons and brain. Genetic modulation of P-body components alters αS toxicity, and human genetic analysis lends support to the disease-relevance of these interactions. Beyond revealing an unexpected aspect of αS function and pathology, our data highlight the versatility of conformationally plastic proteins with high intrinsic disorder.

Amino acid replacements at dozens of positions in the dimeric protein human, Cu,Zn superoxide dismutase (SOD1) can cause amyotrophic lateral sclerosis (ALS). Although it has long been hypothesized that these mutations might enhance the populations of marginally-stable aggregation-prone species responsible for cellular toxicity, there has been little quantitative evidence to support this notion. Perturbations of the folding free energy landscapes of metal-free versions of five ALS-inducing variants, A4V, L38V, G93A, L106V and S134N SOD1, were determined with a global analysis of kinetic and thermodynamic folding data for dimeric and stable monomeric versions of these variants. Utilizing this global analysis approach, the perturbations on the global stability in response to mutation can be partitioned between the monomer folding and association steps, and the effects of mutation on the populations of the folded and unfolded monomeric states can be determined. The 2- to 10-fold increase in the population of the folded monomeric state for A4V, L38V and L106V and the 80- to 480-fold increase in the population of the unfolded monomeric states for all but S134N would dramatically increase their propensity for aggregation through high-order nucleation reactions. The wild-type-like populations of these states for the metal-binding region S134N variant suggest that even wild-type SOD1 may also be prone to aggregation in the absence of metals.

Mutations in GBA, the gene encoding the lysosomal enzyme glucocerebrosidase (GCase), represent the greatest genetic risk factor for developing synucleinopathies including Parkinson's disease (PD). Additionally, PD patients harboring a mutant GBA allele present with an earlier disease onset and an accelerated disease progression of both motor and non-motor symptoms. Preclinical studies in mouse models of synucleinopathy suggest that modulation of the sphingolipid metabolism pathway via inhibition of glucosylceramide synthase (GCS) using a CNS-penetrant small molecule may be a potential treatment for synucleinopathies. Here, we aim to alleviate the lipid storage burden by inhibiting the de novo synthesis of the primary glycosphingolipid substrate of GCase, glucosylceramide (GlcCer). We have previously shown that systemic GCS inhibition reduced GlcCer and glucosylsphingosine (GlcSph) accumulation, slowed α-synuclein buildup in the hippocampus, and improved cognitive deficits. Here, we studied the efficacy of a brain-penetrant clinical candidate GCS inhibitor, venglustat, in mouse models of GBA-related synucleinopathy, including a heterozygous Gba mouse model which more closely replicates the typical GBA-PD patient genotype. Collectively, these data support the rationale for modulation of GCase-related sphingolipid metabolism as a therapeutic strategy for treating GBA-related synucleinopathies.

Protein aggregation is a hallmark of many diseases but also underlies a wide range of positive cellular functions. This phenomenon has been difficult to study because of a lack of quantitative and high-throughput cellular tools. Here, we develop a synthetic genetic tool to sense and control protein aggregation. We apply the technology to yeast prions, developing sensors to track their aggregation states and employing prion fusions to encode synthetic memories in yeast cells. Utilizing high-throughput screens, we identify prion-curing mutants and engineer "anti-prion drives" that reverse the non-Mendelian inheritance pattern of prions and eliminate them from yeast populations. We extend our technology to yeast RNA-binding proteins (RBPs) by tracking their propensity to aggregate, searching for co-occurring aggregates, and uncovering a group of coalescing RBPs through screens enabled by our platform. Our work establishes a quantitative, high-throughput, and generalizable technology to study and control diverse protein aggregation processes in cells.

Gene tagging with fluorescent proteins is commonly applied to investigate the localization and dynamics of proteins in their cellular environment. Ideally, a fluorescent tag is genetically inserted at the endogenous locus at the N- or C- terminus of the gene of interest without disrupting regulatory sequences including the 5' and 3' untranslated region (UTR) and without introducing any extraneous unwanted "scar" sequences, which may create unpredictable transcriptional or translational effects. We present a reliable, low-cost, and highly efficient method for the construction of such scarless C-terminal and N-terminal fusions with fluorescent proteins in yeast. The method relies on sequential positive and negative selection and uses an integration cassette with long flanking regions, which is assembled by two-step PCR, to increase the homologous recombination frequency. The method also enables scarless tagging of essential genes with no need for a complementing plasmid. To further ease high-throughput strain construction, we have computationally automated design of the primers, applied the primer design code to all open reading frames (ORFs) of the budding yeast Saccharomyces cerevisiae (S. cerevisiae) and the fission yeast Schizosaccharomyces pombe (S. pombe), and provide here the computed sequences. To illustrate the scarless N- and C-terminal gene tagging methods in S. cerevisiae, we tagged various genes including the E3 ubiquitin ligase RSP5, the proteasome subunit PRE1, and the eleven Rab GTPases with yeast codon-optimized mNeonGreen or mCherry; several of these represent essential genes. We also implemented the scarless C-terminal gene tagging method in the distantly related organism S. pombe using kanMX6 and HSV1tk as positive and negative selection markers, respectively, as well as ura4. The scarless gene tagging methods presented here are widely applicable to visualize and investigate the functional roles of proteins in living cells.

Heat shock factor (Hsf1) regulates the expression of molecular chaperones to maintain protein homeostasis. Despite its central role in stress resistance, disease and aging, the mechanisms that control Hsf1 activity remain unresolved. Here we show that in budding yeast, Hsf1 basally associates with the chaperone Hsp70 and this association is transiently disrupted by heat shock, providing the first evidence that a chaperone repressor directly regulates Hsf1 activity. We develop and experimentally validate a mathematical model of Hsf1 activation by heat shock in which unfolded proteins compete with Hsf1 for binding to Hsp70. Surprisingly, we find that Hsf1 phosphorylation, previously thought to be required for activation, in fact only positively tunes Hsf1 and does so without affecting Hsp70 binding. Our work reveals two uncoupled forms of regulation - an ON/OFF chaperone switch and a tunable phosphorylation gain - that allow Hsf1 to flexibly integrate signals from the proteostasis network and cell signaling pathways.

As synthetic biology approaches are extended to diverse applications throughout medicine, biotechnology and basic biological research, there is an increasing need to engineer yeast, plant and mammalian cells. Eukaryotic genomes are regulated by the diverse biochemical and biophysical states of chromatin, which brings distinct challenges, as well as opportunities, over applications in bacteria. Recent synthetic approaches, including 'epigenome editing', have allowed the direct and functional dissection of many aspects of physiological chromatin regulation. These studies lay the foundation for biomedical and biotechnological engineering applications that could take advantage of the unique combinatorial and spatiotemporal layers of chromatin regulation to create synthetic systems of unprecedented sophistication.