Context-specific functions of chromatin remodellers in development and disease.

Chromatin remodellers were once thought to be highly redundant and nonspecific in their actions. However, recent human genetic studies demonstrate remarkable biological specificity and dosage sensitivity of the thirty-two adenosine triphosphate (ATP)-dependent chromatin remodellers encoded in the human genome. Mutations in remodellers produce many human developmental disorders and cancers, motivating efforts to investigate their distinct functions in biologically relevant settings. Exquisitely specific biological functions seem to be an emergent property in mammals, and in many cases are based on the combinatorial assembly of subunits and the generation of stable, composite surfaces. Critical interactions between remodelling complex subunits, the nucleosome and other transcriptional regulators are now being defined from structural and biochemical studies. In addition, in vivo analyses of remodellers at relevant genetic loci have provided minute-by-minute insights into their dynamics. These studies are proposing new models for the determinants of remodeller localization and function on chromatin.

Borrowing Transcriptional Kinases to Activate Apoptosis.

Protein kinases are disease drivers whose therapeutic targeting traditionally centers on inhibition of enzymatic activity. Here chemically induced proximity is leveraged to convert kinase inhibitors into context-specific activators of therapeutic genes. Bivalent molecules that link ligands of the transcription factor B-cell lymphoma 6 (BCL6) to ATP-competitive inhibitors of cyclin-dependent kinases (CDKs) were developed to re-localize CDK to BCL6-bound loci on chromatin and direct phosphorylation of RNA Pol II. The resulting BCL6-target proapoptotic gene expression translated into killing of diffuse large B-cell lymphoma (DLBCL) cells at 72 h with EC50s of 0.9 - 10 nM and highly specific ablation of the BCL6-regulated germinal center response in mice. The molecules exhibited 10,000-fold lower cytotoxicity in normal lymphocytes and are well tolerated in mice. Genomic and proteomic evidence corroborated a gain-of-function mechanism where, instead of global enzyme inhibition, a fraction of total kinase activity is borrowed and re-localized to BCL6-bound loci. The strategy demonstrates how kinase inhibitors can be used to context-specifically activate transcription, accessing new therapeutic space.

Rewiring cancer drivers to activate apoptosis.

Genes that drive the proliferation, survival, invasion and metastasis of malignant cells have been identified for many human cancers1-4. Independent studies have identified cell death pathways that eliminate cells for the good of the organism5,6. The coexistence of cell death pathways with driver mutations suggests that the cancer driver could be rewired to activate cell death using chemical inducers of proximity (CIPs). Here we describe a new class of molecules called transcriptional/epigenetic CIPs (TCIPs) that recruit the endogenous cancer driver, or a downstream transcription factor, to the promoters of cell death genes, thereby activating their expression. We focused on diffuse large B cell lymphoma, in which the transcription factor B cell lymphoma 6 (BCL6) is deregulated7. BCL6 binds to the promoters of cell death genes and epigenetically suppresses their expression8. We produced TCIPs by covalently linking small molecules that bind BCL6 to those that bind to transcriptional activators that contribute to the oncogenic program, such as BRD4. The most potent molecule, TCIP1, increases binding of BRD4 by 50% over genomic BCL6-binding sites to produce transcriptional elongation at pro-apoptotic target genes within 15 min, while reducing binding of BRD4 over enhancers by only 10%, reflecting a gain-of-function mechanism. TCIP1 kills diffuse large B cell lymphoma cell lines, including chemotherapy-resistant, TP53-mutant lines, at EC50 of 1-10 nM in 72 h and exhibits cell-specific and tissue-specific effects, capturing the combinatorial specificity inherent to transcription. The TCIP concept also has therapeutic applications in regulating the expression of genes for regenerative medicine and developmental disorders.

BAF subunit switching regulates chromatin accessibility to control cell cycle exit in the developing mammalian cortex.

mSWI/SNF or BAF chromatin regulatory complexes are dosage-sensitive regulators of human neural development frequently mutated in autism spectrum disorders and intellectual disability. Cell cycle exit and differentiation of neural stem/progenitor cells is accompanied by BAF subunit switching to generate neuron-specific nBAF complexes. We manipulated the timing of BAF subunit exchange in vivo and found that early loss of the npBAF subunit BAF53a stalls the cell cycle to disrupt neurogenesis. Loss of BAF53a results in decreased chromatin accessibility at specific neural transcription factor binding sites, including the pioneer factors SOX2 and ASCL1, due to Polycomb accumulation. This results in repression of cell cycle genes, thereby blocking cell cycle progression and differentiation. Cell cycle block upon Baf53a deletion could be rescued by premature expression of the nBAF subunit BAF53b but not by other major drivers of proliferation or differentiation. WNT, EGF, bFGF, SOX2, c-MYC, or PAX6 all fail to maintain proliferation in the absence of BAF53a, highlighting a novel mechanism underlying neural progenitor cell cycle exit in the continued presence of extrinsic proliferative cues.

A CRISPR/Cas9-Engineered ARID1A-Deficient Human Gastric Cancer Organoid Model Reveals Essential and Nonessential Modes of Oncogenic Transformation.

Mutations in ARID1A rank among the most common molecular aberrations in human cancer. However, oncogenic consequences of ARID1A mutation in human cells remain poorly defined due to lack of forward genetic models. Here, CRISPR/Cas9-mediated ARID1A knockout (KO) in primary TP53-/- human gastric organoids induced morphologic dysplasia, tumorigenicity, and mucinous differentiation. Genetic WNT/ß-catenin activation rescued mucinous differentiation, but not hyperproliferation, suggesting alternative pathways of ARID1A KO-mediated transformation. ARID1A mutation induced transcriptional regulatory modules characteristic of microsatellite instability and Epstein-Barr virus-associated subtype human gastric cancer, including FOXM1-associated mitotic genes and BIRC5/survivin. Convergently, high-throughput compound screening indicated selective vulnerability of ARID1A-deficient organoids to inhibition of BIRC5/survivin, functionally implicating this pathway as an essential mediator of ARID1A KO-dependent early-stage gastric tumorigenesis. Overall, we define distinct pathways downstream of oncogenic ARID1A mutation, with nonessential WNT-inhibited mucinous differentiation in parallel with essential transcriptional FOXM1/BIRC5-stimulated proliferation, illustrating the general utility of organoid-based forward genetic cancer analysis in human cells. SIGNIFICANCE: We establish the first human forward genetic modeling of a commonly mutated tumor suppressor gene, ARID1A. Our study integrates diverse modalities including CRISPR/Cas9 genome editing, organoid culture, systems biology, and small-molecule screening to derive novel insights into early transformation mechanisms of ARID1A-deficient gastric cancers.See related commentary by Zafra and Dow, p. 1327.This article is highlighted in the In This Issue feature, p. 1307.

Alternative Splicing of Opioid Receptor Genes Shows a Conserved Pattern for 6TM Receptor Variants.

The opioid receptor (OPR) family comprises the mu-, delta-, and kappa-opioid, and nociceptin receptors that belong to the superfamily of 7-transmembrane spanning G protein-coupled receptors (GPCRs). The mu-opioid receptor is the main target for clinically used opioid analgesics, and its biology has been extensively studied. The N-terminally truncated 6TM receptors isoform produced through alternative splicing of the OPRM1 gene displays unique signaling and analgesic properties, but it is unclear if other OPRs have the same ability. In this study, we have built a comprehensive map of alternative splicing events that produce 6TM receptor variants in all the OPRs and demonstrated their evolutionary conservation. We then obtained evidence for their translation through ribosomal footprint analysis. We discovered that N-terminally truncated 6TM GPCRs are rare in the human genome and OPRs are overrepresented in this group. Finally, we also observed a significant enrichment of 6TM GPCR genes among genes associated with pain, psychiatric disorders, and addiction. Understanding the biology of 6TM receptors and leveraging this knowledge for drug development should pave the way for novel therapies.

CHD8 dosage regulates transcription in pluripotency and early murine neural differentiation.

The chromatin remodeler CHD8 is among the most frequently mutated genes in autism spectrum disorder (ASD). CHD8 has a dosage-sensitive role in ASD, but when and how it becomes critical to human social function is unclear. Here, we conducted genomic analyses of heterozygous and homozygous Chd8 mouse embryonic stem cells and differentiated neural progenitors. We identify dosage-sensitive CHD8 transcriptional targets, sites of regulated accessibility, and an unexpected cooperation with SOX transcription factors. Collectively, our findings reveal that CHD8 negatively regulates expression of neuronal genes to maintain pluripotency and also during differentiation. Thus, CHD8 is essential for both the maintenance of pluripotency and neural differentiation, providing mechanistic insight into its function with potential implications for ASD.

RNA-Puzzles Round IV: 3D structure predictions of four ribozymes and two aptamers.

RNA-Puzzles is a collective endeavor dedicated to the advancement and improvement of RNA 3D structure prediction. With agreement from crystallographers, the RNA structures are predicted by various groups before the publication of the crystal structures. We now report the prediction of 3D structures for six RNA sequences: four nucleolytic ribozymes and two riboswitches. Systematic protocols for comparing models and crystal structures are described and analyzed. In these six puzzles, we discuss (i) the comparison between the automated web servers and human experts; (ii) the prediction of coaxial stacking; (iii) the prediction of structural details and ligand binding; (iv) the development of novel prediction methods; and (v) the potential improvements to be made. We show that correct prediction of coaxial stacking and tertiary contacts is essential for the prediction of RNA architecture, while ligand binding modes can only be predicted with low resolution and simultaneous prediction of RNA structure with accurate ligand binding still remains out of reach. All the predicted models are available for the future development of force field parameters and the improvement of comparison and assessment tools.

Loss of the neural-specific BAF subunit ACTL6B relieves repression of early response genes and causes recessive autism.

Synaptic activity in neurons leads to the rapid activation of genes involved in mammalian behavior. ATP-dependent chromatin remodelers such as the BAF complex contribute to these responses and are generally thought to activate transcription. However, the mechanisms keeping such "early activation" genes silent have been a mystery. In the course of investigating Mendelian recessive autism, we identified six families with segregating loss-of-function mutations in the neuronal BAF (nBAF) subunit ACTL6B (originally named BAF53b). Accordingly, ACTL6B was the most significantly mutated gene in the Simons Recessive Autism Cohort. At least 14 subunits of the nBAF complex are mutated in autism, collectively making it a major contributor to autism spectrum disorder (ASD). Patient mutations destabilized ACTL6B protein in neurons and rerouted dendrites to the wrong glomerulus in the fly olfactory system. Humans and mice lacking ACTL6B showed corpus callosum hypoplasia, indicating a conserved role for ACTL6B in facilitating neural connectivity. Actl6b knockout mice on two genetic backgrounds exhibited ASD-related behaviors, including social and memory impairments, repetitive behaviors, and hyperactivity. Surprisingly, mutation of Actl6b relieved repression of early response genes including AP1 transcription factors (Fos, Fosl2, Fosb, and Junb), increased chromatin accessibility at AP1 binding sites, and transcriptional changes in late response genes associated with early response transcription factor activity. ACTL6B loss is thus an important cause of recessive ASD, with impaired neuron-specific chromatin repression indicated as a potential mechanism.

Computationally Guided Design of Single-Chain Variable Fragment Improves Specificity of Chimeric Antigen Receptors.

Chimeric antigen receptor (CAR)-T cell-based immunotherapy of malignant disease relies on the specificity and association constant of single-chain variable fragments (scFvs). The latter are synthesized from parent antibodies by fusing their light (VL) and heavy (VH)-chain variable domains into a single chain using a flexible linker peptide. The fusion of VL and VH domains can distort their relative orientation, thereby compromising specificity and association constant of scFv, and reducing the lytic efficacy of CAR-T cells. Here, we circumvent the complications of domains' fusion by designing scFv mutants that stabilize interaction between scFv and its target, thereby rescuing scFv efficacy. We employ an iterative approach, based on structural modeling and mutagenesis driven by computational protein design. To demonstrate the power of this approach, we use the scFv derived from an antibody specific to a human leukocyte antigen A2 (HLA-A2)-HER2-derived peptide complex. Whereas the parental antibody is highly specific to its target, the scFv showed reduced specificity. Using our approach, we design mutations into scFvs that restore specificity of the original antibody.

Distinct Binding Modes of Vinculin Isoforms Underlie Their Functional Differences.

Vinculin and its splice isoform metavinculin play key roles in regulating cellular morphology, motility, and force transduction. Vinculin is distinct from metavinculin in its ability to bundle filamentous actin (F-actin). To elucidate the molecular basis for these differences, we employed computational and experimental approaches. Results from these analyses suggest that the C terminus of both vinculin and metavinculin form stable interactions with the F-actin surface. However, the metavinculin tail (MVt) domain contains a 68 amino acid insert, with helix 1 (H1) sequestered into a globular subdomain, which protrudes from the F-actin surface and prevents actin bundling by sterically occluding actin filaments. Consistent with our model, deletion and selective point mutations within the MVt H1 disrupt this protruding structure, and facilitate actin bundling similar to vinculin tail (Vt) domain.

Limits in accuracy and a strategy of RNA structure prediction using experimental information.

RNA structural complexity and flexibility present a challenge for computational modeling efforts. Experimental information and bioinformatics data can be used as restraints to improve the accuracy of RNA tertiary structure prediction. Regarding utilization of restraints, the fundamental questions are: (i) What is the limit in prediction accuracy that one can achieve with arbitrary number of restraints? (ii) Is there a strategy for selection of the minimal number of restraints that would result in the best structural model? We address the first question by testing the limits in prediction accuracy using native contacts as restraints. To address the second question, we develop an algorithm based on the distance variation allowed by secondary structure (DVASS), which ranks restraints according to their importance to RNA tertiary structure prediction. We find that due to kinetic traps, the greatest improvement in the structure prediction accuracy is achieved when we utilize only 40-60% of the total number of native contacts as restraints. When the restraints are sorted by DVASS algorithm, using only the first 20% ranked restraints can greatly improve the prediction accuracy. Our findings suggest that only a limited number of strategically selected distance restraints can significantly assist in RNA structure modeling.

Cardiomyopathy Mutations in Metavinculin Disrupt Regulation of Vinculin-Induced F-Actin Assemblies.

Debilitating heart conditions, notably dilated and hypertrophic cardiomyopathies (CMs), are associated with point mutations in metavinculin, a larger isoform of the essential cytoskeletal protein vinculin. Metavinculin is co-expressed with vinculin at sub-stoichiometric ratios in cardiac tissues. CM mutations in the metavinculin tail domain (MVt) occur within the extra 68-residue insert that differentiates it from the vinculin tail domain (Vt). Vt binds actin filaments (F-actin) and promotes vinculin dimerization to bundle F-actin into thick fibers. While MVt binds to F-actin in a similar manner to Vt, MVt is incapable of F-actin bundling and inhibits Vt-mediated F-actin bundling. We performed F-actin co-sedimentation and negative-stain EM experiments to dissect the coordinated roles of metavinculin and vinculin in actin fiber assembly and the effects of three known metavinculin CM mutations. These CM mutants were found to weakly induce the formation of disordered F-actin assemblies. Notably, they fail to inhibit Vt-mediated F-actin bundling and instead promote formation of large assemblies embedded with linear bundles. Computational models of MVt bound to F-actin suggest that MVt undergoes a conformational change licensing the formation of a protruding sub-domain incorporating the insert, which sterically prevents dimerization and bundling of F-actin by Vt. Sub-domain formation is destabilized by CM mutations, disrupting this inhibitory mechanism. These findings provide new mechanistic insights into the ability of metavinculin to tune actin organization by vinculin and suggest that dysregulation of this process by CM mutants could underlie their malfunction in disease.

Computational design of chemogenetic and optogenetic split proteins.

Controlling protein activity with chemogenetics and optogenetics has proven to be powerful for testing hypotheses regarding protein function in rapid biological processes. Controlling proteins by splitting them and then rescuing their activity through inducible reassembly offers great potential to control diverse protein activities. Building split proteins has been difficult due to spontaneous assembly, difficulty in identifying appropriate split sites, and inefficient induction of effective reassembly. Here we present an automated approach to design effective split proteins regulated by a ligand or by light (SPELL). We develop a scoring function together with an engineered domain to enable reassembly of protein halves with high efficiency and with reduced spontaneous assembly. We demonstrate SPELL by applying it to proteins of various shapes and sizes in living cells. The SPELL server (spell.dokhlab.org) offers an automated prediction of split sites.

RNA Fibers as Optimized Nanoscaffolds for siRNA Coordination and Reduced Immunological Recognition.

RNA is a versatile biomaterial that can be used to engineer nanoassemblies for personalized treatment of various diseases. Despite promising advancements, the design of RNA nanoassemblies with minimal recognition by the immune system remains a major challenge. Here, an approach is reported to engineer RNA fibrous structures to operate as a customizable platform for efficient coordination of siRNAs and for maintaining low immunostimulation. Functional RNA fibers are studied in silico and their formation is confirmed by various experimental techniques and visualized by atomic force microscopy (AFM). It is demonstrated that the RNA fibers offer multiple advantages among which are: i) programmability and modular design that allow for simultaneous controlled delivery of multiple siRNAs and fluorophores, ii) reduced immunostimulation when compared to other programmable RNA nanoassemblies, and iii) simple production protocol for endotoxin-free fibers with the option of their cotranscriptional assembly. Furthermore, it is shown that functional RNA fibers can be efficiently delivered with various organic and inorganic carriers while retaining their structural integrity in cells. Specific gene silencing triggered by RNA fibers is assessed in human breast cancer and melanoma cell lines, with the confirmed ability of functional fibers to selectively target single nucleotide mutations.

Predicting the functional consequences of non-synonymous single nucleotide polymorphisms in IL8 gene.

Here we report an in-silico approach for identification, characterization and validation of deleterious non-synonymous SNPs (nsSNPs) in the interleukin-8 gene using three steps. In first step, sequence homology-based genetic analysis of a set of 50 coding SNPs associated with 41 rsIDs using SIFT (Sorting Intolerant from Tolerant) and PROVEAN (Protein Variation Effect Analyzer) identified 23 nsSNPs to be putatively damaging/deleterious in at least one of the two tools used. Subsequently, structure-homology based PolyPhen-2 (Polymorphism Phenotyping) analysis predicted 9 of 23 nsSNPs (K4T, E31A, E31K, S41Y, I55N, P59L, P59S, L70P and V88D) to be damaging. According to the conditional hypothesis for the study, only nsSNPs that score damaging/deleterious prediction in both sequence and structural homology-based approach will be considered as 'high-confidence' nsSNPs. In step 2, based on conservation of amino acid residues, stability analysis, structural superimposition, RSMD and docking analysis, the possible structural-functional relationship was ascertained for high-confidence nsSNPs. Finally, in a separate analysis (step 3), the IL-8 deregulation has also appeared to be an important prognostic marker for detection of patients with gastric and lung cancer. This study, for the first time, provided in-depth insights on the effects of amino acid substitutions on IL-8 protein structure, function and disease association.

RNA-Puzzles Round III: 3D RNA structure prediction of five riboswitches and one ribozyme.

RNA-Puzzles is a collective experiment in blind 3D RNA structure prediction. We report here a third round of RNA-Puzzles. Five puzzles, 4, 8, 12, 13, 14, all structures of riboswitch aptamers and puzzle 7, a ribozyme structure, are included in this round of the experiment. The riboswitch structures include biological binding sites for small molecules (S-adenosyl methionine, cyclic diadenosine monophosphate, 5-amino 4-imidazole carboxamide riboside 5'-triphosphate, glutamine) and proteins (YbxF), and one set describes large conformational changes between ligand-free and ligand-bound states. The Varkud satellite ribozyme is the most recently solved structure of a known large ribozyme. All puzzles have established biological functions and require structural understanding to appreciate their molecular mechanisms. Through the use of fast-track experimental data, including multidimensional chemical mapping, and accurate prediction of RNA secondary structure, a large portion of the contacts in 3D have been predicted correctly leading to similar topologies for the top ranking predictions. Template-based and homology-derived predictions could predict structures to particularly high accuracies. However, achieving biological insights from de novo prediction of RNA 3D structures still depends on the size and complexity of the RNA. Blind computational predictions of RNA structures already appear to provide useful structural information in many cases. Similar to the previous RNA-Puzzles Round II experiment, the prediction of non-Watson-Crick interactions and the observed high atomic clash scores reveal a notable need for an algorithm of improvement. All prediction models and assessment results are available at http://ahsoka.u-strasbg.fr/rnapuzzles/.