Structure of a pre-catalytic spliceosome

Cryo-EM heterogeneity analysis using regularized covariance estimation and kernel regression

Proteins and the complexes they form are central to nearly all cellular processes. Their flexibility, expressed through a continuum of states, provides a window into their biological functions. Cryogenic electron microscopy (cryo-EM) is an ideal tool to study these dynamic states as it captures specimens in noncrystalline conditions and enables high-resolution reconstructions. However, analyzing the heterogeneous distributions of conformations from cryo-EM data is challenging. We present RECOVAR, a method for analyzing these distributions based on principal component analysis (PCA) computed using a REgularized COVARiance estimator. RECOVAR is fast, robust, interpretable, expressive, and competitive with state-of-the-art neural network methods on heterogeneous cryo-EM datasets. The regularized covariance method efficiently computes a large number of high-resolution principal components that can encode rich heterogeneous distributions of conformations and does so robustly thanks to an automatic regularization scheme. The reconstruction method based on adaptive kernel regression resolves conformational states to a higher resolution than all other tested methods on extensive independent benchmarks while remaining highly interpretable. Additionally, we exploit favorable properties of the PCA embedding to estimate the conformational density accurately. This density allows for better interpretability of the latent space by identifying stable states and low free-energy motions. Finally, we present a scheme to navigate the high-dimensional latent space by automatically identifying these low free-energy trajectories. We make the code freely available at https://github.com/ma-gilles/recovar.

Dynamics and evolutionary conservation of B complex protein recruitment during spliceosome activation

Spliceosome assembly and catalytic activation involve dozens of protein and snRNA binding and unbinding events. The B-complex specific proteins (Prp38, Snu23, and Spp381) have critical roles in stabilizing the spliceosome during rearrangements essential for activation. While these proteins are conserved, different mechanisms have been proposed for their recruitment to spliceosomes. To visualize recruitment directly, we used Colocalization Single Molecule Spectroscopy (CoSMoS) to study the dynamics of Prp38, Snu23, and Spp381 during splicing in real time. These proteins can bind to and release from spliceosomes simultaneously and are likely associated with one another. We designate the assembly minimally containing Prp38, Snu23, and Spp381 as a potential B complex protein (BCP) subcomplex of the spliceosome. Under splicing conditions, BCP proteins associate with pre-mRNA after tri-snRNP binding. BCP protein release predominantly occurs after U4 snRNP dissociation and after NineTeen Complex (NTC) association. Under low concentrations of ATP, BCP proteins preassociate with the tri-snRNP resulting in their simultaneous binding to pre-mRNA. Together, our results reveal that BCP protein recruitment to the spliceosome is conserved between Saccharomyces cerevisiae and humans. The binding of BCP proteins to the tri-snRNP when ATP is limiting may result in the formation of unproductive complexes that could be used to regulate splicing.

CryoSTAR: leveraging structural priors and constraints for cryo-EM heterogeneous reconstruction

Resolving conformational heterogeneity in cryogenic electron microscopy datasets remains an important challenge in structural biology. Previous methods have often been restricted to working exclusively on volumetric densities, neglecting the potential of incorporating any preexisting structural knowledge as prior or constraints. Here we present cryoSTAR, which harnesses atomic model information as structural regularization to elucidate such heterogeneity. Our method uniquely outputs both coarse-grained models and density maps, showcasing the molecular conformational changes at different levels. Validated against four diverse experimental datasets, spanning large complexes, a membrane protein and a small single-chain protein, our results consistently demonstrate an efficient and effective solution to conformational heterogeneity with minimal human bias. By integrating atomic model insights with cryogenic electron microscopy data, cryoSTAR represents a meaningful step forward, paving the way for a deeper understanding of dynamic biological processes.

DynaMight: estimating molecular motions with improved reconstruction from cryo-EM images

How to deal with continuously flexing molecules is one of the biggest outstanding challenges in single-particle analysis of proteins from cryogenic-electron microscopy (cryo-EM) images. Here, we present DynaMight, a software tool that estimates a continuous space of conformations in a cryo-EM dataset by learning three-dimensional deformations of a Gaussian pseudo-atomic model of a consensus structure for every particle image. Inversion of the learned deformations is then used to obtain an improved reconstruction of the consensus structure. We illustrate the performance of DynaMight for several experimental cryo-EM datasets. We also show how error estimates on the deformations may be obtained by independently training two variational autoencoders on half sets of the cryo-EM data, and how regularization of the three-dimensional deformations through the use of atomic models may lead to important artifacts due to model bias. DynaMight is distributed as free, open-source software, as part of RELION-5.

Building molecular model series from heterogeneous CryoEM structures using Gaussian mixture models and deep neural networks

Cryogenic electron microscopy (CryoEM) produces structures of macromolecules at near-atomic resolution. However, building molecular models with good stereochemical geometry from those structures can be challenging and time-consuming, especially when many structures are obtained from datasets with conformational heterogeneity. Here we present a model refinement protocol that automatically generates series of molecular models from CryoEM datasets, which describe the dynamics of the macromolecular system and have near-perfect geometry scores.

Structure and RNA-binding of the helically extended Roquin CCCH-type zinc finger

Zinc finger (ZnF) domains appear in a pool of structural contexts and despite their small size achieve varying target specificities, covering single-stranded and double-stranded DNA and RNA as well as proteins. Combined with other RNA-binding domains, ZnFs enhance affinity and specificity of RNA-binding proteins (RBPs). The ZnF-containing immunoregulatory RBP Roquin initiates mRNA decay, thereby controlling the adaptive immune system. Its unique ROQ domain shape-specifically recognizes stem-looped cis-elements in mRNA 3'-untranslated regions (UTR). The N-terminus of Roquin contains a RING domain for protein-protein interactions and a ZnF, which was suggested to play an essential role in RNA decay by Roquin. The ZnF domain boundaries, its RNA motif preference and its interplay with the ROQ domain have remained elusive, also driven by the lack of high-resolution data of the challenging protein. We provide the solution structure of the Roquin-1 ZnF and use an RBNS-NMR pipeline to show that the ZnF recognizes AU-rich RNAs. We systematically refine the contributions of adenines in a poly(U)-background to specific complex formation. With the simultaneous binding of ROQ and ZnF to a natural target transcript of Roquin, our study for the first time suggests how Roquin integrates RNA shape and sequence features through the ROQ-ZnF tandem.

Cryo-EM Heterogeneity Analysis using Regularized Covariance Estimation and Kernel Regression

Proteins and the complexes they form are central to nearly all cellular processes. Their flexibility, expressed through a continuum of states, provides a window into their biological functions. Cryogenic electron microscopy (cryo-EM) is an ideal tool to study these dynamic states as it captures specimens in non-crystalline conditions and enables high-resolution reconstructions. However, analyzing the heterogeneous distributions of conformations from cryo-EM data is challenging. We present RECOVAR, a method for analyzing these distributions based on principal component analysis (PCA) computed using a REgularized COVARiance estimator. RECOVAR is fast, robust, interpretable, expressive, and competitive with the state-of-art neural network methods on heterogeneous cryo-EM datasets. The regularized covariance method efficiently computes a large number of high-resolution principal components that can encode rich heterogeneous distributions of conformations and does so robustly thanks to an automatic regularization scheme. The novel reconstruction method based on adaptive kernel regression resolves conformational states to a higher resolution than all other tested methods on extensive independent benchmarks while remaining highly interpretable. Additionally, we exploit favorable properties of the PCA embedding to estimate the conformational density accurately. This density allows for better interpretability of the latent space by identifying stable states and low free-energy motions. Finally, we present a scheme to navigate the high-dimensional latent space by automatically identifying these low free-energy trajectories. We make the code freely available at https://github.com/ma-gilles/recovar.

Characterization of Cwc2, U6 snRNA, and Prp8 interactions destabilized by Prp16 ATPase at the transition between the first and second steps of splicing

The spliceosome performs two consecutive transesterification reactions using one catalytic center, thus requiring its rearrangement between the two catalytic steps of splicing. The Prp16 ATPase facilitates exit from the first-step conformation of the catalytic center by destabilizing some interactions important for catalysis. To better understand rearrangements within the Saccharomyces cerevisiae catalytic center, we characterize factors that modulate the function of Prp16: Cwc2, N-terminal domain of Prp8, and U6-41AACAAU46 region. Alleles of these factors were identified through genetic screens for mutants that correct cs defects of prp16-302 alleles. Several of the identified U6, cwc2, and prp8 alleles are located in close proximity of each other in cryo-EM structures of the spliceosomal catalytic conformations. Cwc2 and U6 interact with the intron sequences in the first step, but they do not seem to contribute to the stability of the second-step catalytic center. On the other hand, the N-terminal segment of Prp8 not only affects intron positioning for the first step, but it also makes important contacts in the proximity of the active site for both the first and second steps of splicing. By identifying interactions important for the stability of catalytic conformations, our genetic analyses indirectly inform us about features of the transition-state conformation of the spliceosome.

Dynamics and Evolutionary Conservation of B Complex Protein Recruitment During Spliceosome Activation

Spliceosome assembly and catalytic site formation (called activation) involve dozens of protein and snRNA binding and unbinding events. The B-complex specific proteins Prp38, Snu23, and Spp381 have critical roles in stabilizing the spliceosome during conformational changes essential for activation. While these proteins are conserved, different mechanisms have been proposed for their recruitment to spliceosomes. To visualize recruitment directly, we used Colocalization Single Molecule Spectroscopy (CoSMoS) to study the dynamics of Prp38, Snu23, and Spp381 during splicing in real time. These proteins bind to and release from spliceosomes simultaneously and are likely associated with one another. We designate the complex of Prp38, Snu23, and Spp381 as the B Complex Protein (BCP) subcomplex. Under splicing conditions, the BCP associates with pre-mRNA after tri-snRNP binding. BCP release predominantly occurs after U4 snRNP dissociation and after NineTeen Complex (NTC) association. Under low concentrations of ATP, the BCP pre-associates with the tri-snRNP resulting in their simultaneous binding to pre-mRNA. Together, our results reveal that the BCP recruitment pathway to the spliceosome is conserved between S. cerevisiae and humans. Binding of the BCP to the tri-snRNP when ATP is limiting may result in formation of unproductive complexes that could be used to regulate splicing.

Mechanism for the initiation of spliceosome disassembly

Precursor-mRNA (pre-mRNA) splicing requires the assembly, remodelling and disassembly of the multi-megadalton ribonucleoprotein complex called the spliceosome1. Recent studies have shed light on spliceosome assembly and remodelling for catalysis2-6, but the mechanism of disassembly remains unclear. Here we report cryo-electron microscopy structures of nematode and human terminal intron lariat spliceosomes along with biochemical and genetic data. Our results uncover how four disassembly factors and the conserved RNA helicase DHX15 initiate spliceosome disassembly. The disassembly factors probe large inner and outer spliceosome surfaces to detect the release of ligated mRNA. Two of these factors, TFIP11 and C19L1, and three general spliceosome subunits, SYF1, SYF2 and SDE2, then dock and activate DHX15 on the catalytic U6 snRNA to initiate disassembly. U6 therefore controls both the start5 and end of pre-mRNA splicing. Taken together, our results explain the molecular basis of the initiation of canonical spliceosome disassembly and provide a framework to understand general spliceosomal RNA helicase control and the discard of aberrant spliceosomes.

Molecular basis for the activation of human spliceosome

The spliceosome executes pre-mRNA splicing through four sequential stages: assembly, activation, catalysis, and disassembly. Activation of the spliceosome, namely remodeling of the pre-catalytic spliceosome (B complex) into the activated spliceosome (Bact complex) and the catalytically activated spliceosome (B* complex), involves major flux of protein components and structural rearrangements. Relying on a splicing inhibitor, we have captured six intermediate states between the B and B* complexes: pre-Bact, Bact-I, Bact-II, Bact-III, Bact-IV, and post-Bact. Their cryo-EM structures, together with an improved structure of the catalytic step I spliceosome (C complex), reveal how the catalytic center matures around the internal stem loop of U6 snRNA, how the branch site approaches 5'-splice site, how the RNA helicase PRP2 rearranges to bind pre-mRNA, and how U2 snRNP undergoes remarkable movement to facilitate activation. We identify a previously unrecognized key role of PRP2 in spliceosome activation. Our study recapitulates a molecular choreography of the human spliceosome during its catalytic activation.

Mechanisms and regulation of spliceosome-mediated pre-mRNA splicing in Saccharomyces cerevisiae

Pre-mRNA splicing, the removal of introns and ligation of flanking exons, is a crucial step in eukaryotic gene expression. The spliceosome, a macromolecular complex made up of five small nuclear RNAs (snRNAs) and dozens of proteins, assembles on introns via a complex pathway before catalyzing the two transesterification reactions necessary for splicing. All of these steps have the potential to be highly regulated to ensure correct mRNA isoform production for proper cellular function. While Saccharomyces cerevisiae (yeast) has a limited set of intron-containing genes, many of these genes are highly expressed, resulting in a large number of transcripts in a cell being spliced. As a result, splicing regulation is of critical importance for yeast. Just as in humans, yeast splicing can be influenced by protein components of the splicing machinery, structures and properties of the pre-mRNA itself, or by the action of trans-acting factors. It is likely that further analysis of the mechanisms and pathways of splicing regulation in yeast can reveal general principles applicable to other eukaryotes. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Processing > Splicing Regulation/Alternative Splicing.

Data-driven regularization lowers the size barrier of cryo-EM structure determination

Macromolecular structure determination by electron cryo-microscopy (cryo-EM) is limited by the alignment of noisy images of individual particles. Because smaller particles have weaker signals, alignment errors impose size limitations on its applicability. Here, we explore how image alignment is improved by the application of deep learning to exploit prior knowledge about biological macromolecular structures that would otherwise be difficult to express mathematically. We train a denoising convolutional neural network on pairs of half-set reconstructions from the electron microscopy data bank (EMDB) and use this denoiser as an alternative to a commonly used smoothness prior. We demonstrate that this approach, which we call Blush regularization, yields better reconstructions than do existing algorithms, in particular for data with low signal-to-noise ratios. The reconstruction of a protein-nucleic acid complex with a molecular weight of 40 kDa, which was previously intractable, illustrates that denoising neural networks will expand the applicability of cryo-EM structure determination for a wide range of biological macromolecules.

Growth-regulated co-occupancy of Mediator and Lsm3 at intronic ribosomal protein genes

Mediator is a well-known transcriptional co-regulator and serves as an adaptor between gene-specific regulatory proteins and RNA polymerase II. Studies on the chromatin-bound form of Mediator revealed interactions with additional protein complexes involved in various transcription-related processes, such as the Lsm2-8 complex that is part of the spliceosomal U6 small nuclear ribonucleoprotein complex. Here, we employ Chromatin Immunoprecipitation sequencing (ChIP-seq) of chromatin associated with the Lsm3 protein and the Med1 or Med15 Mediator subunits. We identify 86 genes co-occupied by both Lsm3 and Mediator, of which 73 were intron-containing ribosomal protein genes. In logarithmically growing cells, Mediator primarily binds to their promoter regions but also shows a second, less pronounced occupancy at their 3'-exons. During the late exponential phase, we observe a near-complete transition of Mediator from these promoters to a position in their 3'-ends, overlapping the Lsm3 binding sites ∼250 bp downstream of their last intron-exon boundaries. Using an unbiased RNA sequencing approach, we show that transition of Mediator from promoters to the last exon of these genes correlates to reduction of both their messenger RNA levels and splicing ratios, indicating that the Mediator and Lsm complexes cooperate to control growth-regulated expression of intron-containing ribosomal protein genes at the levels of transcription and splicing.

Conserved intronic secondary structures with concealed branch sites regulate alternative splicing of poison exons

Alternative splicing (AS) generates multiple RNA isoforms and increases the complexities of transcriptomes and proteomes. However, it remains unclear how RNA structures contribute to AS regulation. Here, we systematically search transcriptomes for secondary structures with concealed branch sites (BSs) in the alternatively spliced introns and predict thousands of them from six organisms, of which many are evolutionarily conserved. Intriguingly, a highly conserved stem-loop structure with concealed BSs is found in animal SF3B3 genes and colocalizes with a downstream poison exon (PE). Destabilization of this structure allows increased usage of the BSs and results in enhanced PE inclusion in human and Drosophila cells, leading to decreased expression of SF3B3. This structure is experimentally validated using an in-cell SHAPE-MaP assay. Through RNA interference screens of 28 RNA-binding proteins, we find that this stem-loop structure is sensitive to U2 factors. Furthermore, we find that SF3B3 also facilitates DNA repair and protects genome stability by enhancing interaction between ERCC6/CSB and arrested RNA polymerase II. Importantly, both Drosophila and human cells with the secondary structure mutated by genome editing exhibit altered DNA repair in vivo. This study provides a novel and common mechanism for AS regulation of PEs and reveals a physiological function of SF3B3 in DNA repair.

Deep generative priors for biomolecular 3D heterogeneous reconstruction from cryo-EM projections

Cryo-electron microscopy has become a powerful tool to determine three-dimensional (3D) structures of rigid biological macromolecules from noisy micrographs with single-particle reconstruction. Recently, deep neural networks, e.g., CryoDRGN, have demonstrated conformational and compositional heterogeneity of complexes. However, the lack of ground-truth conformations poses a challenge to assess the performance of heterogeneity analysis methods. In this work, variational autoencoders (VAE) with three types of deep generative priors were learned for latent variable inference and heterogeneous 3D reconstruction via Bayesian inference. More specifically, VAEs with "Variational Mixture of Posteriors" priors (VampPrior-SPR), non-parametric exemplar-based priors (ExemplarPrior-SPR) and priors from latent score-based generative models (LSGM-SPR) were quantitatively compared with CryoDRGN. We built four simulated datasets composed of hypothetical continuous conformation or discrete states of the hERG K + channel. Empirical and quantitative comparisons of inferred latent representations were performed with affine-transformation-based metrics. These models with more informative priors gave better regularized, interpretable factorized latent representations with better conserved pairwise distances, less deformed latent distributions and lower within-cluster variances. They were also tested on experimental datasets to resolve compositional and conformational heterogeneity (50S ribosome assembly, cowpea chlorotic mottle virus, and pre-catalytic spliceosome) with comparable high resolution. Codes and data are available: https://github.com/benjamin3344/DGP-SPR.

Molecular mechanisms of non-genetic aberrant peptide production in cancer

The cancer peptidome has long been known to be altered by genetic mutations. However, more recently, non-genetic polypeptide mutations have also been related to cancer cells. These non-genetic mutations occur post-t30ranscriptionally, leading to the modification of the peptide primary structure, while the corresponding genes remain unchanged. Three main processes participate in the production of these aberrant proteins: mRNA alternative splicing, mRNA editing, and mRNA aberrant translation. In this review, we summarize the molecular mechanisms underlying these processes and the recent findings on the functions of the aberrant proteins, as well as their exploitability as new therapeutic targets due to their specific enrichment in cancer cells. These non-genetic aberrant polypeptides represent a source of novel cancer cell targets independent from their level of mutational burden, still to be exhaustively explored.

Environmentally relevant concentrations of 2,3,7,8-TCDD induced inhibition of multicellular alternative splicing and transcriptional dysregulation

Alternative splicing (AS), found in approximately 95 % of human genes, significantly amplifies protein diversity and is implicated in disease pathogenesis when dysregulated. However, the precise involvement of AS in the toxic mechanisms induced by TCDD (2,3,7,8-tetrachlorodibenzo-p-dioxin) remains incompletely elucidated. This study conducted a thorough global AS analysis in six human cell lines following TCDD exposure. Our findings revealed that environmentally relevant concentration (0.1 nM) of TCDD significantly suppressed AS events in all cell types, notably inhibiting diverse splicing events and reducing transcript diversity, potentially attributed to modifications in the splicing patterns of the inhibitory factor family, particularly hnRNP. And we identified 151 genes with substantial AS alterations shared among these cell types, particularly enriched in immune and metabolic pathways. Moreover, TCDD induced cell-specific changes in splicing patterns and transcript levels, with increased sensitivity notably in THP-1 monocyte, potentially linked to aberrant expression of pivotal genes within the spliceosome pathway (DDX5, EFTUD2, PUF60, RBM25, SRSF1, and CRNKL1). This study extends our understanding of disrupted alternative splicing and its relation to the multisystem toxicity of TCDD. It sheds light on how environmental toxins affect post-transcriptional regulatory processes, offering a fresh perspective for toxicology and disease etiology investigations.

Cryo-EM analyses of dimerized spliceosomes provide new insights into the functions of B complex proteins

The B complex is a key intermediate stage of spliceosome assembly. To improve the structural resolution of monomeric, human spliceosomal B (hB) complexes and thereby generate a more comprehensive hB molecular model, we determined the cryo-EM structure of B complex dimers formed in the presence of ATP γ S. The enhanced resolution of these complexes allows a finer molecular dissection of how the 5' splice site (5'ss) is recognized in hB, and new insights into molecular interactions of FBP21, SNU23 and PRP38 with the U6/5'ss helix and with each other. It also reveals that SMU1 and RED are present as a heterotetrameric complex and are located at the interface of the B dimer protomers. We further show that MFAP1 and UBL5 form a 5' exon binding channel in hB, and elucidate the molecular contacts stabilizing the 5' exon at this stage. Our studies thus yield more accurate models of protein and RNA components of hB complexes. They further allow the localization of additional proteins and protein domains (such as SF3B6, BUD31 and TCERG1) whose position was not previously known, thereby uncovering new functions for B-specific and other hB proteins during pre-mRNA splicing.

OPUS-DSD: deep structural disentanglement for cryo-EM single-particle analysis

Cryo-electron microscopy (cryo-EM) captures snapshots of dynamic macromolecules, collectively illustrating the involved structural landscapes. This provides an exciting opportunity to explore the structural variations of macromolecules under study. However, traditional cryo-EM single-particle analysis often yields static structures. Here we describe OPUS-DSD, an algorithm capable of efficiently reconstructing the structural landscape embedded in cryo-EM data. OPUS-DSD uses a three-dimensional convolutional encoder-decoder architecture trained with cryo-EM images, thereby encoding structural variations into a smooth and easily analyzable low-dimension space. This space can be traversed to reconstruct continuous dynamics or clustered to identify distinct conformations. OPUS-DSD can offer meaningful insights into the structural variations of macromolecules, filling in the gaps left by traditional cryo-EM structural determination, and potentially improves the reconstruction resolution by reliably clustering similar particles within the dataset. These functionalities are especially relevant to the study of highly dynamic biological systems. OPUS-DSD is available at https://github.com/alncat/opusDSD .

The SMN complex drives structural changes in human snRNAs to enable snRNP assembly

Spliceosomal snRNPs are multicomponent particles that undergo a complex maturation pathway. Human Sm-class snRNAs are generated as 3'-end extended precursors, which are exported to the cytoplasm and assembled together with Sm proteins into core RNPs by the SMN complex. Here, we provide evidence that these pre-snRNA substrates contain compact, evolutionarily conserved secondary structures that overlap with the Sm binding site. These structural motifs in pre-snRNAs are predicted to interfere with Sm core assembly. We model structural rearrangements that lead to an open pre-snRNA conformation compatible with Sm protein interaction. The predicted rearrangement pathway is conserved in Metazoa and requires an external factor that initiates snRNA remodeling. We show that the essential helicase Gemin3, which is a component of the SMN complex, is crucial for snRNA structural rearrangements during snRNP maturation. The SMN complex thus facilitates ATP-driven structural changes in snRNAs that expose the Sm site and enable Sm protein binding.

Role of the RNA-binding protein family in gynecologic cancers

Gynecological cancers pose a threat to women's health. Although early-stage gynecological cancers show good outcomes after standardized treatment, the prognosis of patients with advanced, met-astatic, and recurrent cancers is poor. RNA-binding proteins (RBPs) are important cellular proteins that interact with RNA through RNA-binding domains and participate extensively in post-transcriptional regulatory processes, such as mRNA alternative splicing, polyadenylation, intracellular localization and stability, and translation. Abnormal RBP expression affects the normal function of oncogenes and tumor suppressor genes in many malignancies, thus leading to the occurrence or progression of cancers. Similarly, RBPs play crucial roles in gynecological carcinogenesis. We summarize the role of RBPs in gynecological malignancies and explore their potential in the diagnosis and treatment of cancers. The findings summarized in this review may provide a guide for future research on the functions of RBPs.

A network of DZF proteins controls alternative splicing regulation and fidelity

Proteins containing DZF (domain associated with zinc fingers) modules play important roles throughout gene expression, from transcription to translation. Derived from nucleotidyltransferases but lacking catalytic residues, DZF domains serve as heterodimerization surfaces between DZF protein pairs. Three DZF proteins are widely expressed in mammalian tissues, ILF2, ILF3 and ZFR, which form mutually exclusive ILF2-ILF3 and ILF2-ZFR heterodimers. Using eCLIP-Seq, we find that ZFR binds across broad intronic regions to regulate the alternative splicing of cassette and mutually exclusive exons. ZFR preferentially binds dsRNA in vitro and is enriched on introns containing conserved dsRNA elements in cells. Many splicing events are similarly altered upon depletion of any of the three DZF proteins; however, we also identify independent and opposing roles for ZFR and ILF3 in alternative splicing regulation. Along with widespread involvement in cassette exon splicing, the DZF proteins control the fidelity and regulation of over a dozen highly validated mutually exclusive splicing events. Our findings indicate that the DZF proteins form a complex regulatory network that leverages dsRNA binding by ILF3 and ZFR to modulate splicing regulation and fidelity.

Therapeutic Targeting of RNA Splicing in Cancer

RNA splicing is a key regulatory step in the proper control of gene expression. It is a highly dynamic process orchestrated by the spliceosome, a macro-molecular machinery that consists of protein and RNA components. The dysregulation of RNA splicing has been observed in many human pathologies ranging from neurodegenerative diseases to cancer. The recent identification of recurrent mutations in the core components of the spliceosome in hematologic malignancies has advanced our knowledge of how splicing alterations contribute to disease pathogenesis. This review article will discuss our current understanding of how aberrant RNA splicing regulation drives tumor initiation and progression. We will also review current therapeutic modalities and highlight emerging technologies designed to target RNA splicing for cancer treatment.

Integrating Molecular Models Into CryoEM Heterogeneity Analysis Using Scalable High-resolution Deep Gaussian Mixture Models

Resolving the structural variability of proteins is often key to understanding the structure-function relationship of those macromolecular machines. Single particle analysis using Cryogenic electron microscopy (CryoEM), combined with machine learning algorithms, provides a way to reveal the dynamics within the protein system from noisy micrographs. Here, we introduce an improved computational method that uses Gaussian mixture models for protein structure representation and deep neural networks for conformation space embedding. By integrating information from molecular models into the heterogeneity analysis, we can analyze continuous protein conformational changes using structural information at the frequency of 1/3 Å-1, and present the results in a more interpretable form.

Spliceosome assembly and regulation: insights from analysis of highly reduced spliceosomes

Premessenger RNA splicing is catalyzed by the spliceosome, a multimegadalton RNA-protein complex that assembles in a highly regulated process on each intronic substrate. Most studies of splicing and spliceosomes have been carried out in human or S. cerevisiae model systems. There exists, however, a large diversity of spliceosomes, particularly in organisms with reduced genomes, that suggests a means of analyzing the essential elements of spliceosome assembly and regulation. In this review, we characterize changes in spliceosome composition across phyla, describing those that are most frequently observed and highlighting an analysis of the reduced spliceosome of the red alga Cyanidioschyzon merolae We used homology modeling to predict what effect splicing protein loss would have on the spliceosome, based on currently available cryo-EM structures. We observe strongly correlated loss of proteins that function in the same process, for example, in interacting with the U1 snRNP (which is absent in C. merolae), regulation of Brr2, or coupling transcription and splicing. Based on our observations, we predict splicing in C. merolae to be inefficient, inaccurate, and post-transcriptional, consistent with the apparent trend toward its elimination in this lineage. This work highlights the striking flexibility of the splicing pathway and the spliceosome when viewed in the context of eukaryotic diversity.

Conformational space exploration of cryo-EM structures by variability refinement

Cryo-EM observation of biological samples enables visualization of sample heterogeneity, in the form of discrete states that are separable, or continuous heterogeneity as a result of local protein motion before flash freezing. Variability analysis of this continuous heterogeneity describes the variance between a particle stack and a volume, and results in a map series describing the various steps undertaken by the sample in the particle stack. While this observation is absolutely stunning, it is very hard to pinpoint structural details to elements of the maps. In order to bridge the gap between observation and explanation, we designed a tool that refines an ensemble of structures into all the maps from variability analysis. Using this bundle of structures, it is easy to spot variable parts of the structure, as well as the parts that are not moving. Comparison with molecular dynamics simulations highlights the fact that the movements follow the same directions, albeit with different amplitudes. Ligand can also be investigated using this method. Variability refinement is available in the Phenix software suite, accessible under the program name phenix.varref.

Post-Transcriptional Modification by Alternative Splicing and Pathogenic Splicing Variants in Cardiovascular Development and Congenital Heart Defects

Advancements in genomics, bioinformatics, and genome editing have uncovered new dimensions in gene regulation. Post-transcriptional modifications by the alternative splicing of mRNA transcripts are critical regulatory mechanisms of mammalian gene expression. In the heart, there is an expanding interest in elucidating the role of alternative splicing in transcriptome regulation. Substantial efforts were directed toward investigating this process in heart development and failure. However, few studies shed light on alternative splicing products and their dysregulation in congenital heart defects (CHDs). While elegant reports showed the crucial roles of RNA binding proteins (RBPs) in orchestrating splicing transitions during heart development and failure, the impact of RBPs dysregulation or genetic variation on CHDs has not been fully addressed. Herein, we review the current understanding of alternative splicing and RBPs' roles in heart development and CHDs. Wediscuss the impact of perinatal splicing transition and its dysregulation in CHDs. We further summarize the discoveries made of causal splicing variants in key transcription factors that are implicated in CHDs. An improved understanding of the roles of alternative splicing in heart development and CHDs may potentially inform novel preventive and therapeutic advancements for newborn infants with CHDs.

Estimating conformational landscapes from Cryo-EM particles by 3D Zernike polynomials

The new developments in Cryo-EM Single Particle Analysis are helping us to understand how the macromolecular structure and function meet to drive biological processes. By capturing many states at the particle level, it is possible to address how macromolecules explore different conformations, information that is classically extracted through 3D classification. However, the limitations of classical approaches prevent us from fully understanding the complete conformational landscape due to the reduced number of discrete states accurately reconstructed. To characterize the whole structural spectrum of a macromolecule, we propose an extension of our Zernike3D approach, able to extract per-image continuous flexibility information directly from a particle dataset. Also, our method can be seamlessly applied to images, maps or atomic models, opening integrative possibilities. Furthermore, we introduce the ZART reconstruction algorithm, which considers the Zernike3D deformation fields to revert particle conformational changes during the reconstruction process, thus minimizing the blurring induced by molecular motions.

Frozen in time: analyzing molecular dynamics with time-resolved cryo-EM

Molecular machines, such as polymerases, ribosomes, or proteasomes, fulfill complex tasks requiring the thermal energy of their environment. They achieve this by restricting random motion along a path of possible conformational changes. These changes are often directed through engagement with different cofactors, which can best be compared to a Brownian ratchet. Many molecular machines undergo three major steps throughout their functional cycles, including initialization, repetitive processing, and termination. Several of these major states have been elucidated by cryogenic electron microscopy (cryo-EM). However, the individual steps for these machines are unique and multistep processes themselves, and their coordination in time is still elusive. To measure these short-lived intermediate events by cryo-EM, the total reaction time needs to be shortened to enrich for the respective pre-equilibrium states. This approach is termed time-resolved cryo-EM (trEM). In this review, we sum up the methodological development of trEM and its application to a range of biological questions.

Recruitment of a splicing factor to the nuclear lamina for its inactivation

Precursor messenger RNA splicing is a highly regulated process, mediated by a complex RNA-protein machinery, the spliceosome, that encompasses several hundred proteins and five small nuclear RNAs in humans. Emerging evidence suggests that the spatial organization of splicing factors and their spatio-temporal dynamics participate in the regulation of splicing. So far, methods to manipulate the spatial distribution of splicing factors in a temporally defined manner in living cells are missing. Here, we describe such an approach that takes advantage of a reversible chemical dimerizer, and outline the requirements for efficient, reversible re-localization of splicing factors to selected sub-nuclear compartments. In a proof-of-principle study, the partial re-localization of the PRPF38A protein to the nuclear lamina in HEK293T cells induced a moderate increase in intron retention. Our approach allows fast and reversible re-localization of splicing factors, has few side effects and can be applied to many splicing factors by fusion of a protein tag through genome engineering. Apart from the systematic analysis of the spatio-temporal aspects of splicing regulation, the approach has a large potential for the fast induction and reversal of splicing switches and can reveal mechanisms of splicing regulation in native nuclear environments.

Through the use of a reversible chemical dimerizer, the splicing factor PRPF38A is re-localized to the nuclear lamina, paving the way for a systematic analysis of spatio-temporal splicing regulation.

Amortized Inference for Heterogeneous Reconstruction in Cryo-EM

Cryo-electron microscopy (cryo-EM) is an imaging modality that provides unique insights into the dynamics of proteins and other building blocks of life. The algorithmic challenge of jointly estimating the poses, 3D structure, and conformational heterogeneity of a biomolecule from millions of noisy and randomly oriented 2D projections in a computationally efficient manner, however, remains unsolved. Our method, cryoFIRE, performs ab initio heterogeneous reconstruction with unknown poses in an amortized framework, thereby avoiding the computationally expensive step of pose search while enabling the analysis of conformational heterogeneity. Poses and conformation are jointly estimated by an encoder while a physics-based decoder aggregates the images into an implicit neural representation of the conformational space. We show that our method can provide one order of magnitude speedup on datasets containing millions of images without any loss of accuracy. We validate that the joint estimation of poses and conformations can be amortized over the size of the dataset. For the first time, we prove that an amortized method can extract interpretable dynamic information from experimental datasets.

Identification of transient intermediates during spliceosome activation by single molecule fluorescence microscopy

Spliceosome activation is the process of creating the catalytic site for RNA splicing and occurs de novo on each intron following spliceosome assembly. Dozens of factors bind to or are released from the activating spliceosome including the Lsm2-8 heteroheptameric ring that binds the U6 small nuclear RNA 3'-end. Lsm2-8 must be released to permit active site stabilization by the Prp19-containing complex (NineTeen Complex, NTC); however, little is known about the temporal order of events and dynamic interactions that lead up to and follow Lsm2-8 release. We have used colocalization single molecule spectroscopy (CoSMoS) to visualize Lsm2-8 dynamics during activation of Saccharomyces cerevisiae spliceosomes in vitro. Lsm2-8 is recruited as a component of the tri-snRNP and is released after integration of the Prp19-containing complex (NTC). Despite Lsm2-8 and the NTC being mutually exclusive in existing cryo-EM structures of yeast B complex spliceosomes, we identify a transient intermediate containing both ([Formula: see text]) and provide a kinetic framework for its formation and transformation during activation. Prior to [Formula: see text] assembly, the NTC rapidly and reversibly samples the spliceosome suggesting a mechanism for preventing NTC sequestration by defective spliceosomes that fail to properly activate. In complementary ensemble assays, we show that a base-pairing-dependent ternary complex can form between Lsm2-8 and U2 and U6 helix II RNAs. We propose that this interaction may play a role in formation of transient spliceosome intermediates formed during activation.

Structural and functional investigation of the human snRNP assembly factor AAR2 in complex with the RNase H-like domain of PRPF8

The crystal structure of human AAR2 bound to the central spliceosomal factor PRPF8 and in vitro functional data yield insights into the structural basis of snRNP assembly in humans.

Small nuclear ribonucleoprotein complexes (snRNPs) represent the main subunits of the spliceosome. While the assembly of the snRNP core particles has been well characterized, comparably little is known of the incorporation of snRNP-specific proteins and the mechanisms of snRNP recycling. U5 snRNP assembly in yeast requires binding of the the Aar2 protein to Prp8p as a placeholder to preclude premature assembly of the SNRNP200 helicase, but the role of the human AAR2 homolog has not yet been investigated in detail. Here, a crystal structure of human AAR2 in complex with the RNase H-like domain of the U5-specific PRPF8 (PRP8F RH) is reported, revealing a significantly different interaction between the two proteins compared with that in yeast. Based on the structure of the AAR2–PRPF8 RH complex, the importance of the interacting regions and residues was probed and AAR2 variants were designed that failed to stably bind PRPF8 in vitro. Protein-interaction studies of AAR2 with U5 proteins using size-exclusion chromatography reveal similarities and marked differences in the interaction patterns compared with yeast Aar2p and imply phosphorylation-dependent regulation of AAR2 reminiscent of that in yeast. It is found that in vitro AAR2 seems to lock PRPF8 RH in a conformation that is only compatible with the first transesterification step of the splicing reaction and blocks a conformational switch to the step 2-like, Mg²⁺-coordinated conformation that is likely during U5 snRNP biogenesis. These findings extend the picture of AAR2 PRP8 interaction from yeast to humans and indicate a function for AAR2 in the spliceosomal assembly process beyond its role as an SNRNP200 placeholder in yeast.

CryoAI: Amortized Inference of Poses for Ab Initio Reconstruction of 3D Molecular Volumes from Real Cryo-EM Images

Cryo-electron microscopy (cryo-EM) has become a tool of fundamental importance in structural biology, helping us understand the basic building blocks of life. The algorithmic challenge of cryo-EM is to jointly estimate the unknown 3D poses and the 3D electron scattering potential of a biomolecule from millions of extremely noisy 2D images. Existing reconstruction algorithms, however, cannot easily keep pace with the rapidly growing size of cryo-EM datasets due to their high computational and memory cost. We introduce cryoAI, an ab initio reconstruction algorithm for homogeneous conformations that uses direct gradient-based optimization of particle poses and the electron scattering potential from single-particle cryo-EM data. CryoAI combines a learned encoder that predicts the poses of each particle image with a physics-based decoder to aggregate each particle image into an implicit representation of the scattering potential volume. This volume is stored in the Fourier domain for computational efficiency and leverages a modern coordinate network architecture for memory efficiency. Combined with a symmetrized loss function, this framework achieves results of a quality on par with state-of-the-art cryo-EM solvers for both simulated and experimental data, one order of magnitude faster for large datasets and with significantly lower memory requirements than existing methods.

Emerging Methods and Applications to Decrypt Allostery in Proteins and Nucleic Acids

Many large protein-nucleic acid complexes exhibit allosteric regulation. In these systems, the propagation of the allosteric signaling is strongly coupled to conformational dynamics and catalytic function, challenging state-of-the-art analytical methods. Here, we review established and innovative approaches used to elucidate allosteric mechanisms in these complexes. Specifically, we report network models derived from graph theory and centrality analyses in combination with molecular dynamics (MD) simulations, introducing novel schemes that implement the synergistic use of graph theory with enhanced simulations methods and ab-initio MD. Accelerated MD simulations are used to construct "enhanced network models", describing the allosteric response over long timescales and capturing the relation between allostery and conformational changes. "Ab-initio network models" combine graph theory with ab-initio MD and quantum mechanics/molecular mechanics (QM/MM) simulations to describe the allosteric regulation of catalysis by following the step-by-step dynamics of biochemical reactions. This approach characterizes how the allosteric regulation changes from reactants to products and how it affects the transition state, revealing a tense-to-relaxed allosteric regulation along the chemical step. Allosteric models and applications are showcased for three paradigmatic examples of allostery in protein-nucleic acid complexes: (i) the nucleosome core particle, (ii) the CRISPR-Cas9 genome editing system and (iii) the spliceosome. These methods and applications create innovative protocols to determine allosteric mechanisms in protein-nucleic acid complexes that show tremendous promise for medicine and bioengineering.

Cell-cycle and apoptosis related and proteomics-based signaling pathways of human hepatoma Huh-7 cells treated by three currently used multi-RTK inhibitors

Sorafenib, lenvatinib and regorafenib, the multi-RTK inhibitors with potent anti-angiogenesis effects, are currently therapeutic drugs generally recommended for the patients with advanced hepatocellular carcinoma (HCC). To date, however, there have been no published studies on the mechanism underling differential effects of the three drugs on HCC cell proliferation, and the proteomic analysis in HCC cell lines treated by regorafenib or lenvatinib. The present study for the first time performed a direct comparison of the cell cycle arrest and apoptosis induction in the Huh-7 cells caused by sorafenib, regorafenib and lenvatinib at respective IC₅₀ using flow cytometry technique, as well as their pharmacological interventions for influencing whole cell proteomics using tandem mass tag-based peptide-labeling coupled with the nLC-HRMS technique. Sorafenib, regorafenib and lenvatinib at respective IC₅₀ drove the remaining surviving Huh-7 cells into a G₀/G₁ arrest, but lenvatinib and regorafenib were much more effective than sorafenib. Lenvatinib produced a much stronger induction of Huh-7 cells into early apoptosis than sorafenib and regorafenib, while necrotic cell proportion induced by regorafenib was 2.4 times as large as that by lenvatinib. The proteomic study revealed 419 proteins downregulated commonly by the three drugs at respective IC₅₀. KEGG pathway analysis of the downregulated proteins indicated the ranking of top six signaling pathways including the spliceosome, DNA replication, cell cycle, mRNA surveillance, P53 and nucleotide excision repair involved in 33 proteins, all of which were directly related to their pharmacological effects on cell cycle and cell apoptosis. Notably, lenvatinib and regorafenib downregulated the proteins of PCNA, Cyclin B1, BCL-xL, TSP1, BUD31, SF3A1 and Mad2 much more strongly than sorafenib. Moreover, most of the proteins in the P53 signaling pathway were downregulated with lenvatinib and regorafenib by more than 36% at least. In conclusion, lenvatinib and regorafenib have much stronger potency against Huh-7 cell proliferation than sorafenib because of their more potent effects on cell cycle arrest and apoptosis induction. The underling mechanism may be at least due to the 33 downregulated proteins centralizing the signal pathways of cell cycle, p53 and DNA synthesis based on the present proteomics study.

Visualizing Conformational Space of Functional Biomolecular Complexes by Deep Manifold Learning

The cellular functions are executed by biological macromolecular complexes in nonequilibrium dynamic processes, which exhibit a vast diversity of conformational states. Solving the conformational continuum of important biomolecular complexes at the atomic level is essential to understanding their functional mechanisms and guiding structure-based drug discovery. Here, we introduce a deep manifold learning framework, named AlphaCryo4D, which enables atomic-level cryogenic electron microscopy (cryo-EM) reconstructions that approximately visualize the conformational space of biomolecular complexes of interest. AlphaCryo4D integrates 3D deep residual learning with manifold embedding of pseudo-energy landscapes, which simultaneously improves 3D classification accuracy and reconstruction resolution via an energy-based particle-voting algorithm. In blind assessments using simulated heterogeneous datasets, AlphaCryo4D achieved 3D classification accuracy three times those of alternative methods and reconstructed continuous conformational changes of a 130-kDa protein at sub-3 Å resolution. By applying this approach to analyze several experimental datasets of the proteasome, ribosome and spliceosome, we demonstrate its potential generality in exploring hidden conformational space or transient states of macromolecular complexes that remain hitherto invisible. Integration of this approach with time-resolved cryo-EM further allows visualization of conformational continuum in a nonequilibrium regime at the atomic level, thus potentially enabling therapeutic discovery against highly dynamic biomolecular targets.

Extended DNA binding interfaces beyond the canonical SAP domain contribute to the function of replication stress regulator SDE2 at DNA replication forks

Elevated DNA replication stress causes instability of the DNA replication fork and increased DNA mutations, which underlies tumorigenesis. The DNA replication stress regulator silencing-defective 2 (SDE2) is known to bind to TIMELESS (TIM), a protein of the fork protection complex, and enhances its stability, thereby supporting replisome activity at DNA replication forks. However, the DNA-binding activity of SDE2 is not well defined. Here, we structurally and functionally characterize a new conserved DNA-binding motif related to the SAP (SAF-A/B, Acinus, PIAS) domain in human SDE2 and establish its preference for ssDNA. Our NMR solution structure of the SDE2^SAP domain reveals a helix-extended loop-helix core with the helices aligned parallel to each other, consistent with known canonical SAP folds. Notably, we have shown that the DNA interaction of this SAP domain extends beyond the core SAP domain and is augmented by two lysine residues in the C-terminal tail, which is uniquely positioned adjacent to the SAP motif and conserved in the pre-mRNA splicing factor SF3A3. Furthermore, we found that mutation in the SAP domain and extended C terminus not only disrupts ssDNA binding but also impairs TIM localization at replication forks, thus inhibiting efficient fork progression. Taken together, our results establish SDE2^SAP as an essential element for SDE2 to exert its role in preserving replication fork integrity via fork protection complex regulation and highlight the structural diversity of the DNA–protein interactions achieved by a specialized DNA-binding motif.

Developing Graphene Grids for Cryoelectron Microscopy

Cryogenic electron microscopy (cryo-EM) single particle analysis has become one of the major techniques used to study high-resolution 3D structures of biological macromolecules. Specimens are generally prepared in a thin layer of vitrified ice using a holey carbon grid. However, the sample quality using this type of grid is not always ideal for high-resolution imaging even when the specimens in the test tube behave ideally. Various problems occur during a vitrification procedure, including poor/nonuniform distribution of particles, preferred orientation of particles, specimen denaturation/degradation, high background from thick ice, and beam-induced motion, which have become important bottlenecks in high-resolution structural studies using cryo-EM in many projects. In recent years, grids with support films made of graphene and its derivatives have been developed to efficiently solve these problems. Here, the various advantages of graphene grids over conventional holey carbon film grids, functionalization of graphene support films, production methods of graphene grids, and origins of pristine graphene contamination are reviewed and discussed.

A UHM - ULM interface with unusual structural features contributes to U2AF2 and SF3B1 association for pre-mRNA splicing

During spliceosome assembly, the 3′ splice site is recognized by sequential U2AF2 complexes, first with Splicing Factor 1 (SF1) and second by the SF3B1 subunit of the U2 small nuclear ribonuclear protein particle. The U2AF2–SF1 interface is well characterized, comprising a U2AF homology motif (UHM) of U2AF2 bound to a U2AF ligand motif (ULM) of SF1. However, the structure of the U2AF2–SF3B1 interface and its importance for pre-mRNA splicing are unknown. To address this knowledge gap, we determined the crystal structure of the U2AF2 UHM bound to a SF3B1 ULM site at 1.8-Å resolution. We discovered a distinctive trajectory of the SF3B1 ULM across the U2AF2 UHM surface, which differs from prior UHM/ULM structures and is expected to modulate the orientations of the full-length proteins. We established that the binding affinity of the U2AF2 UHM for the cocrystallized SF3B1 ULM rivals that of a nearly full-length U2AF2 protein for an N-terminal SF3B1 region. An additional SF3B6 subunit had no detectable effect on the U2AF2–SF3B1 binding affinities. We further showed that key residues at the U2AF2 UHM–SF3B1 ULM interface contribute to coimmunoprecipitation of the splicing factors. Moreover, disrupting the U2AF2–SF3B1 interface changed splicing of representative human transcripts. From analysis of genome-wide data, we found that many of the splice sites coregulated by U2AF2 and SF3B1 differ from those coregulated by U2AF2 and SF1. Taken together, these findings support distinct structural and functional roles for the U2AF2—SF1 and U2AF2—SF3B1 complexes during the pre-mRNA splicing process.

Zinc finger structure determination by NMR: Why zinc fingers can be a handful

Zinc fingers can be loosely defined as protein domains containing one or more tetrahedrally-co-ordinated zinc ions whose role is to stabilise the structure rather than to be involved in enzymatic chemistry; such zinc ions are often referred to as "structural zincs". Although structural zincs can occur in proteins of any size, they assume particular significance for very small protein domains, where they are often essential for maintaining a folded state. Such small structures, that sometimes have only marginal stability, can present particular difficulties in terms of sample preparation, handling and structure determination, and early on they gained a reputation for being resistant to crystallisation. As a result, NMR has played a more prominent role in structural studies of zinc finger proteins than it has for many other types of proteins. This review will present an overview of the particular issues that arise for structure determination of zinc fingers by NMR, and ways in which these may be addressed.

Hierarchical autoclassification of cryo-EM samples and macromolecular energy landscape determination

BACKGROUND AND OBJECTIVE

Cryo-electron microscopy using single particle analysis is a powerful technique for obtaining 3D reconstructions of macromolecules in near native conditions. One of its major advances is its capacity to reveal conformations of dynamic molecular complexes. Most popular and successful current approaches to analyzing heterogeneous complexes are founded on Bayesian inference. However, these 3D classification methods require the tuning of specific parameters by the user and the use of complicated 3D re-classification procedures for samples affected by extensive heterogeneity. Thus, the success of these approaches highly depends on the user experience. We introduce a robust approach to identify many different conformations presented in a cryo-EM dataset based on Bayesian inference through Relion classification methods that does not require tuning of parameters and reclassification strategies.

METHODS

The algorithm allows both 2D and 3D classification and is based on a hierarchical clustering approach that runs automatically without requiring typical inputs, such as the number of conformations present in the dataset or the required classification iterations. This approach is applied to robustly determine the energy landscapes of macromolecules.

RESULTS

We tested the performance of the methods proposed here using four different datasets, comprising structurally homogeneous and highly heterogeneous cases. In all cases, the approach provided excellent results. The routines are publicly available as part of the CryoMethods plugin included in the Scipion package.

CONCLUSIONS

Our results show that the proposed method can be used to align and classify homogeneous and heterogeneous datasets without requiring previous alignment information or any prior knowledge about the number of co-existing conformations. The approach can be used for both 2D and 3D autoclassification and only requires an initial volume. In addition, the approach is robust to the "attractor" problem providing many different conformations/views for samples affected by extensive heterogeneity. The obtained 3D classes can render high resolution 3D structures, while the obtained energy landscapes can be used to determine structural trajectories.

Sequence-specific RNA recognition by an RGG motif connects U1 and U2 snRNP for spliceosome assembly

In mammals, the structural basis for the interaction between U1 and U2 small nuclear ribonucleoproteins (snRNPs) during the early steps of splicing is still elusive. The binding of the ubiquitin-like (UBL) domain of SF3A1 to the stem-loop 4 of U1 snRNP (U1-SL4) contributes to this interaction. Here, we determined the 3D structure of the complex between the UBL of SF3A1 and U1-SL4 RNA. Our crystallography, NMR spectroscopy, and cross-linking mass spectrometry data show that SF3A1-UBL recognizes, sequence specifically, the GCG/CGC RNA stem and the apical UUCG tetraloop of U1-SL4. In vitro and in vivo mutational analyses support the observed intermolecular contacts and demonstrate that the carboxyl-terminal arginine-glycine-glycine-arginine (RGGR) motif of SF3A1-UBL binds sequence specifically by inserting into the RNA major groove. Thus, the characterization of the SF3A1-UBL/U1-SL4 complex expands the repertoire of RNA binding domains and reveals the capacity of RGG/RG motifs to bind RNA in a sequence-specific manner.

The spliceosome component Usp39 controls B cell development by regulating immunoglobulin gene rearrangement

The spliceosome is a large ribonucleoprotein complex responsible for pre-mRNA splicing and genome stability maintenance. Disruption of the spliceosome activity may lead to developmental disorders and tumorigenesis. However, the physiological role that the spliceosome plays in B cell development and function is still poorly defined. Here, we demonstrate that ubiquitin-specific peptidase 39 (Usp39), a spliceosome component of the U4/U6.U5 tri-snRNP complex, is essential for B cell development. Ablation of Usp39 in B cell lineage blocks pre-pro-B to pro-B cell transition in the bone marrow, leading to a profound reduction of mature B cells in the periphery. We show that Usp39 specifically regulates immunoglobulin gene rearrangement in a spliceosome-dependent manner, which involves modulating chromatin interactions at the Igh locus. Moreover, our results indicate that Usp39 deletion reduces the pre-malignant B cells in Eμ-Myc transgenic mice and significantly improves their survival.

Structural basis of branch site recognition by the human spliceosome

Recognition of the intron branch site (BS) by the U2 small nuclear ribonucleoprotein (snRNP) is a critical event during spliceosome assembly. In mammals, BS sequences are poorly conserved, and unambiguous intron recognition cannot be achieved solely through a base-pairing mechanism. We isolated human 17S U2 snRNP and reconstituted in vitro its adenosine 5´-triphosphate (ATP)–dependent remodeling and binding to the pre–messenger RNA substrate. We determined a series of high-resolution (2.0 to 2.2 angstrom) structures providing snapshots of the BS selection process. The substrate-bound U2 snRNP shows that SF3B6 stabilizes the BS:U2 snRNA duplex, which could aid binding of introns with poor sequence complementarity. ATP-dependent remodeling uncoupled from substrate binding captures U2 snRNA in a conformation that competes with BS recognition, providing a selection mechanism based on branch helix stability.

Network theory reveals principles of spliceosome structure and dynamics

Cryoelectron microscopy has revolutionized spliceosome structural biology, and structures representing much of the splicing process have been determined. Comparison of these structures is challenging due to extreme dynamics of the splicing machinery and the thousands of changing interactions during splicing. We have used network theory to analyze splicing factor interactions by constructing structure-based networks from protein-protein, protein-RNA, and RNA-RNA interactions found in eight different spliceosome structures. Our analyses reveal that connectivity dynamics result in step-specific impacts of factors on network topology. The spliceosome's connectivity is focused on the active site, in part due to contributions from nonglobular proteins. Many essential factors exhibit large shifts in centralities during splicing. Others show transiently high betweenness centralities at certain stages, thereby suggesting mechanisms for regulating splicing by briefly bridging otherwise poorly connected network nodes. These observations provide insights into organizing principles of the spliceosome and provide frameworks for comparative analysis of other macromolecular machines.

Single-particle cryo-EM reveals conformational variability of the oligomeric VCC β-barrel pore in a lipid bilayer

Vibrio cholerae cytolysin (VCC) is a water-soluble, membrane-damaging, pore-forming toxin (PFT) secreted by pathogenic V. cholerae, which causes eukaryotic cell death by altering the plasma membrane permeability. VCC self-assembles on the cell surface and undergoes a dramatic conformational change from prepore to heptameric pore structure. Over the past few years, several high-resolution structures of detergent-solubilized PFTs have been characterized. However, high-resolution structural characterization of small β-PFTs in a lipid environment is still rare. Therefore, we used single-particle cryo-EM to characterize the structure of the VCC oligomer in large unilamellar vesicles, which is the first atomic-resolution cryo-EM structure of VCC. From our study, we were able to provide the first documented visualization of the rim domain amino acid residues of VCC interacting with lipid membrane. Furthermore, cryo-EM characterization of lipid bilayer–embedded VCC suggests interesting conformational variabilities, especially in the transmembrane channel, which could have a potential impact on the pore architecture and assist us in understanding the pore formation mechanism.

U2 snRNA structure is influenced by SF3A and SF3B proteins but not by SF3B inhibitors

U2 snRNP is an essential component of the spliceosome. It is responsible for branch point recognition in the spliceosome A-complex via base-pairing of U2 snRNA with an intron to form the branch helix. Small molecule inhibitors target the SF3B component of the U2 snRNP and interfere with A-complex formation during spliceosome assembly. We previously found that the first SF3B inhibited-complex is less stable than A-complex and hypothesized that SF3B inhibitors interfere with U2 snRNA secondary structure changes required to form the branch helix. Using RNA chemical modifiers, we probed U2 snRNA structure in A-complex and SF3B inhibited splicing complexes. The reactivity pattern for U2 snRNA in the SF3B inhibited-complex is indistinguishable from that of A-complex, suggesting that they have the same secondary structure conformation, including the branch helix. This observation suggests SF3B inhibited-complex instability does not stem from an alternate RNA conformation and instead points to the inhibitors interfering with protein component interactions that normally stabilize U2 snRNP’s association with an intron. In addition, we probed U2 snRNA in the free U2 snRNP in the presence of SF3B inhibitor and again saw no differences. However, increased protection of nucleotides upstream of Stem I in the absence of SF3A and SF3B proteins suggests a change of secondary structure at the very 5′ end of U2 snRNA. Chemical probing of synthetic U2 snRNA in the absence of proteins results in similar protections and predicts a previously uncharacterized extension of Stem I. Because this stem must be disrupted for SF3A and SF3B proteins to stably join the snRNP, the structure has the potential to influence snRNP assembly and recycling after spliceosome disassembly.