Scaffolding in the Spliceosome via Single α Helices

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.

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.

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.

Blocking Osa-miR1871 enhances rice resistance against Magnaporthe oryzae and yield

MicroRNAs (miRNAs) play vital roles in plant development and defence responses against various stresses. Here, we show that blocking miR1871 improves rice resistance against Magnaporthe oryzae and enhances grain yield simultaneously. The transgenic lines overexpressing miR1871 (OX1871) exhibit compromised resistance, suppressed defence responses and reduced panicle number resulting in slightly decreased yield. In contrast, the transgenic lines blocking miR1871 (MIM1871) show improved resistance, enhanced defence responses and significantly increased panicle number leading to enhanced yield per plant. The RNA-seq assay and defence response assays reveal that blocking miR1871 resulted in the enhancement of PAMP-triggered immunity (PTI). Intriguingly, miR1871 suppresses the expression of LOC_Os06g22850, which encodes a microfibrillar-associated protein (MFAP1) locating nearby the cell wall and positively regulating PTI responses. The mutants of MFAP1 resemble the phenotype of OX1871. Conversely, the transgenic lines overexpressing MFAP1 (OXMFAP1) or overexpressing both MFAP1 and miR1871 (OXMFAP1/OX1871) resemble the resistance of MIM1871. The time-course experiment data reveal that the expression of miR1871 and MFAP1 in rice leaves, panicles and basal internode is dynamic during the whole growth period to manipulate the resistance and yield traits. Our results suggest that miR1871 regulates rice yield and immunity via MFAP1, and the miR8171-MFAP1 module could be used in rice breeding to improve both immunity and yield.

A multi-factor trafficking site on the spliceosome remodeling enzyme BRR2 recruits C9ORF78 to regulate alternative splicing

The intrinsically unstructured C9ORF78 protein was detected in spliceosomes but its role in splicing is presently unclear. We find that C9ORF78 tightly interacts with the spliceosome remodeling factor, BRR2, in vitro. Affinity purification/mass spectrometry and RNA UV-crosslinking analyses identify additional C9ORF78 interactors in spliceosomes. Cryogenic electron microscopy structures reveal how C9ORF78 and the spliceosomal B complex protein, FBP21, wrap around the C-terminal helicase cassette of BRR2 in a mutually exclusive manner. Knock-down of C9ORF78 leads to alternative NAGNAG 3'-splice site usage and exon skipping, the latter dependent on BRR2. Inspection of spliceosome structures shows that C9ORF78 could contact several detected spliceosome interactors when bound to BRR2, including the suggested 3'-splice site regulating helicase, PRPF22. Together, our data establish C9ORF78 as a late-stage splicing regulatory protein that takes advantage of a multi-factor trafficking site on BRR2, providing one explanation for suggested roles of BRR2 during splicing catalysis and alternative splicing.

Tyramide signal amplification mass spectrometry (TSA-MS) ratio identifies nuclear speckle proteins

We present a simple ratio method to infer protein composition within cellular structures using proximity labeling approaches but compensating for the diffusion of free radicals. We used tyramide signal amplification (TSA) and label-free mass spectrometry (MS) to compare proteins in nuclear speckles versus centromeres. Our “TSA-MS ratio” approach successfully identified known nuclear speckle proteins. For example, 96% and 67% of proteins in the top 30 and 100 sorted proteins, respectively, are known nuclear speckle proteins, including proteins that we validated here as enriched in nuclear speckles. We show that MFAP1, among the top 20 in our list, forms droplets under certain circumstances and that MFAP1 expression levels modulate the size, stability, and dynamics of nuclear speckles. Localization of MFAP1 and its binding partner, PRPF38A, in droplet-like nuclear bodies precedes formation of nuclear speckles during telophase. Our results update older proteomic studies of nuclear speckles and should provide a useful reference dataset to guide future experimental dissection of nuclear speckle structure and function.

Large-scale discovery of protein interactions at residue resolution using co-evolution calculated from genomic sequences

Increasing numbers of protein interactions have been identified in high-throughput experiments, but only a small proportion have solved structures. Recently, sequence coevolution-based approaches have led to a breakthrough in predicting monomer protein structures and protein interaction interfaces. Here, we address the challenges of large-scale interaction prediction at residue resolution with a fast alignment concatenation method and a probabilistic score for the interaction of residues. Importantly, this method (EVcomplex2) is able to assess the likelihood of a protein interaction, as we show here applied to large-scale experimental datasets where the pairwise interactions are unknown. We predict 504 interactions de novo in the E. coli membrane proteome, including 243 that are newly discovered. While EVcomplex2 does not require available structures, coevolving residue pairs can be used to produce structural models of protein interactions, as done here for membrane complexes including the Flagellar Hook-Filament Junction and the Tol/Pal complex.

Our understanding of the residue-level details of protein interactions remains incomplete. Here, the authors show sequence coevolution can be used to infer interacting proteins with residue-level details, including predicting 467 interactions de novo in the Escherichia coli cell envelope proteome.

ER/K-link-Leveraging a native protein linker to probe dynamic cellular interactions

ER/K α-helices are a subset of single alpha helical domains, which exhibit unusual stability as isolated protein secondary structures. They adopt an elongated structural conformation, while regulating the frequency of interactions between proteins or polypeptides fused to their ends. Here we review recent advances on the structure, stability and function of ER/K α-helices as linkers (ER/K linkers) in native proteins. We describe methodological considerations in the molecular cloning, protein expression and measurement of interaction strengths, using sensors incorporating ER/K linkers. We highlight biological insights obtained over the last decade by leveraging distinct biophysical features of ER/K-linked sensors. We conclude with the outlook for the use of ER/K linkers in the selective modulation of dynamic cellular interactions.

Depletion of the MFAP1/SPP381 Splicing Factor Causes R-Loop-Independent Genome Instability

THO/TREX is a conserved complex with a role in messenger ribonucleoprotein biogenesis that links gene expression and genome instability. Here, we show that human THO interacts with MFAP1 (microfibrillar-associated protein 1), a spliceosome-associated factor. Interestingly, MFAP1 depletion impairs cell proliferation and genome integrity, increasing γH2AX foci and DNA breaks. This phenotype is not dependent on either transcription or RNA-DNA hybrids. Mutations in the yeast orthologous gene SPP381 cause similar transcription-independent genome instability, supporting a conserved role. MFAP1 depletion has a wide effect on splicing and gene expression in human cells, determined by transcriptome analyses. MFAP1 depletion affects a number of DNA damage response (DDR) genes, which supports an indirect role of MFAP1 on genome integrity. Our work defines a functional interaction between THO and RNA processing and argues that splicing factors may contribute to genome integrity indirectly by regulating the expression of DDR genes rather than by a direct role.

Molecular structure and function of microfibrillar-associated proteins in skeletal and metabolic disorders and cancers

Microfibrillar-associated proteins (MFAPs) are extracellular matrix glycoproteins, which play a role in microfibril assembly, elastinogenesis, and tissue homeostasis. MFAPs consist of five subfamily members, including MFAP1, MFAP2, MFAP3, MFAP4, and MFAP5. Among these, MFAP2 and MFAP5 are most closely related, and exhibit very limited amino acid sequence homology with MFAP1, MFAP3, and MFAP4. Gene expression profiling analysis reveals that MFAP2, MFAP5, and MFAP4 are specifically expressed in osteoblastic like cells, whereas MFAP1 and MFAP3 are more ubiquitously expressed, indicative of their diverse role in the tropism of tissues. Molecular structural analysis shows that each MFAP family member has distinct features, and functional evidence reveals discrete purposes of individual MFAPs. Animal studies indicate that MFAP2-deficient mice exhibit progressive osteopenia with elevated receptor activator of NF-κB ligand (RANKL) expression, whereas MFAP5-deficient mice are neutropenic, and MFAP4-deficient mice displayed emphysema-like pathology and the impaired formation of neointimal hyperplasia. Emerging data also suggest that MFAPs are involved in cancer progression and fat metabolism. Further understanding of tissue-specific pathophysiology of MFAPs might offer potential novel therapeutic targets for related diseases, such as skeletal and metabolic disorders, and cancers.

Enrichment of Aurora B kinase at the inner kinetochore controls outer kinetochore assembly

Outer kinetochore assembly enables chromosome attachment to microtubules and spindle assembly checkpoint (SAC) signaling in mitosis. Aurora B kinase controls kinetochore assembly by phosphorylating the Mis12 complex (Mis12C) subunit Dsn1. Current models propose Dsn1 phosphorylation relieves autoinhibition, allowing Mis12C binding to inner kinetochore component CENP-C. Using Xenopus laevis egg extracts and biochemical reconstitution, we found that autoinhibition of the Mis12C by Dsn1 impedes its phosphorylation by Aurora B. Our data indicate that the INCENP central region increases Dsn1 phosphorylation by enriching Aurora B at inner kinetochores, close to CENP-C. Furthermore, centromere-bound CENP-C does not exchange in mitosis, and CENP-C binding to the Mis12C dramatically increases Dsn1 phosphorylation by Aurora B. We propose that the coincidence of Aurora B and CENP-C at inner kinetochores ensures the fidelity of kinetochore assembly. We also found that the central region is required for the SAC beyond its role in kinetochore assembly, suggesting that kinetochore enrichment of Aurora B promotes the phosphorylation of other kinetochore substrates.

Remarkable Rigidity of the Single α-Helical Domain of Myosin-VI As Revealed by NMR Spectroscopy

Although the α-helix has long been recognized as an all-important element of secondary structure, it generally requires stabilization by tertiary interactions with other parts of a protein's structure. Highly charged single α-helical (SAH) domains, consisting of a high percentage (>75%) of Arg, Lys, and Glu residues, are exceptions to this rule but have been difficult to characterize structurally. Our study focuses on the 68-residue medial tail domain of myosin-VI, which is found to contain a highly ordered α-helical structure extending from Glu-6 to Lys-63. High hydrogen exchange protection factors (15-150), small (ca. 4 Hz) 3 JHNHα couplings, and a near-perfect fit to an ideal model α-helix for its residual dipolar couplings (RDCs), measured in a filamentous phage medium, support the high regularity of this helix. Remarkably, the hydrogen exchange rates are far more homogeneous than the protection factors derived from them, suggesting that for these transiently broken helices the intrinsic exchange rates derived from the amino acid sequence are not appropriate reference values. 15N relaxation data indicate a very high degree of rotational diffusion anisotropy ( D∥/ D⊥ ≈ 7.6), consistent with the hydrodynamic behavior predicted for such a long, nearly straight α-helix. Alignment of the helix by a paramagnetic lanthanide ion attached to its N-terminal region shows a decrease in alignment as the distance from the tagging site increases. This decrease yields a precise measure for the persistence length of 224 ± 10 Å at 20 °C, supporting the idea that the role of the SAH helix is to act as an extension of the myosin-VI lever arm.

Computational Prediction of MoRFs, Short Disorder-to-order Transitioning Protein Binding Regions

Molecular recognition features (MoRFs) are short protein-binding regions that undergo disorder-to-order transitions (induced folding) upon binding protein partners. These regions are abundant in nature and can be predicted from protein sequences based on their distinctive sequence signatures. This first-of-its-kind survey covers 14 MoRF predictors and six related methods for the prediction of short protein-binding linear motifs, disordered protein-binding regions and semi-disordered regions. We show that the development of MoRF predictors has accelerated in the recent years. These predictors depend on machine learning-derived models that were generated using training datasets where MoRFs are annotated using putative disorder. Our analysis reveals that they generate accurate predictions. We identified eight methods that offer area under the ROC curve (AUC) ≥ 0.7 on experimentally-validated test datasets. We show that modern MoRF predictors accurately find experimentally annotated MoRFs even though they were trained using the putative disorder annotations. They are relatively highly-cited, particularly the methods available as webservers that on average secure three times more citations than methods without this option. MoRF predictions contribute to the experimental discovery of protein-protein interactions, annotation of protein functions and computational analysis of a variety of proteomes, protein families, and pathways. We outline future development and application directions for these tools, stressing the importance to develop novel tools that would target interactions of disordered regions with other types of partners.

Dynamic ion pair behavior stabilizes single α-helices in proteins

Ion pairs are key stabilizing interactions between oppositely charged amino acid side chains in proteins. They are often depicted as single conformer salt bridges (hydrogen-bonded ion pairs) in crystal structures, but it is unclear how dynamic they are in solution. Ion pairs are thought to be particularly important in stabilizing single α-helix (SAH) domains in solution. These highly stable domains are rich in charged residues (such as Arg, Lys, and Glu) with potential ion pairs across adjacent turns of the helix. They provide a good model system to investigate how ion pairs can contribute to protein stability. Using NMR spectroscopy, small-angle X-ray light scattering (SAXS), and molecular dynamics simulations, we provide here experimental evidence that ion pairs exist in a SAH in murine myosin 7a (residues 858-935), but that they are not fixed or long lasting. In silico modeling revealed that the ion pairs within this α-helix exhibit dynamic behavior, rapidly forming and breaking and alternating between different partner residues. The low-energy helical state was compatible with a great variety of ion pair combinations. Flexible ion pair formation utilizing a subset of those available at any one time avoided the entropic penalty of fixing side chain conformations, which likely contributed to helix stability overall. These results indicate the dynamic nature of ion pairs in SAHs. More broadly, thermodynamic stability in other proteins is likely to benefit from the dynamic behavior of multi-option solvent-exposed ion pairs.

Intrinsically Disordered Protein Ntr2 Modulates the Spliceosomal RNA Helicase Brr2

Precursor messenger RNA splicing is mediated by the spliceosome, a large and dynamic molecular machine composed of five small nuclear RNAs and numerous proteins. Many spliceosomal proteins are predicted to be intrinsically disordered or contain large disordered regions, but experimental validation of these predictions is scarce, and the precise functions of these proteins are often unclear. Here, we show via circular dichroism spectroscopy, dynamic light scattering, and NMR spectroscopy that the yeast spliceosomal disassembly factor Ntr2 is largely intrinsically disordered. Peptide SPOT analyses, analytical size-exclusion chromatography, and surface plasmon resonance measurements revealed that Ntr2 uses an N-terminal region to bind the C-terminal helicase unit of the Brr2 RNA helicase, an enzyme involved in spliceosome activation and implicated in splicing catalysis and spliceosome disassembly. NMR analyses suggested that Ntr2 does not adopt a tertiary structure and likely remains disordered upon complex formation. RNA binding and unwinding studies showed that Ntr2 downregulates Brr2 helicase activity in vitro by modulating the fraction of helicase molecules productively bound to the RNA substrate. Our data clarify the nature of a physical link between Brr2 and Ntr2, and point to the possibility of a functional Ntr2-Brr2 interplay during splicing.

Detection of single alpha-helices in large protein sequence sets using hardware acceleration

Single alpha-helices (SAHs) are increasingly recognized as important structural and functional elements of proteins. Comprehensive identification of SAH segments in large protein datasets was largely hindered by the slow speed of the most restrictive prediction tool for their identification, FT_CHARGE on common hardware. We have previously implemented an FPGA-based version of this tool allowing fast analysis of a large number of sequences. Using this implementation, we have set up of a semi-automated pipeline capable of analyzing full UniProt releases in reasonable time and compiling monthly updates of a comprehensive database of SAH segments. Releases of this database, denoted CSAHDB, is available on the CSAHserver 2 website at csahserver.itk.ppke.hu. An overview of human SAH-containing sequences combined with a literature survey suggests specific roles of SAH segments in proteins involved in RNA-based regulation processes as well as cytoskeletal proteins, a number of which is also linked to the development and function of synapses.

Cryo-EM Studies of Pre-mRNA Splicing: From Sample Preparation to Model Visualization

The removal of noncoding introns from pre-messenger RNA (pre-mRNA) is an essential step in eukaryotic gene expression and is catalyzed by a dynamic multi-megadalton ribonucleoprotein complex called the spliceosome. The spliceosome assembles on pre-mRNA substrates by the stepwise addition of small nuclear ribonucleoprotein particles and numerous protein factors. Extensive remodeling is required to form the RNA-based active site and to mediate the pre-mRNA branching and ligation reactions. In the past two years, cryo-electron microscopy (cryo-EM) structures of spliceosomes captured in different assembly and catalytic states have greatly advanced our understanding of its mechanism. This was made possible by long-standing efforts in the purification of spliceosome intermediates as well as recent developments in cryo-EM imaging and computational methodology. The resulting high-resolution densities allow for de novo model building in core regions of the complexes. In peripheral and less ordered regions, the combination of cross-linking, bioinformatics, biochemical, and genetic data is essential for accurate modeling. Here, we summarize these achievements and highlight the critical steps in obtaining near-atomic resolution structures of the spliceosome.

Waggawagga-CLI: A command-line tool for predicting stable single α-helices (SAH-domains), and the SAH-domain distribution across eukaryotes

Stable single-alpha helices (SAH-domains) function as rigid connectors and constant force springs between structural domains, and can provide contact surfaces for protein-protein and protein-RNA interactions. SAH-domains mainly consist of charged amino acids and are monomeric and stable in polar solutions, characteristics which distinguish them from coiled-coil domains and intrinsically disordered regions. Although the number of reported SAH-domains is steadily increasing, genome-wide analyses of SAH-domains in eukaryotic genomes are still missing. Here, we present Waggawagga-CLI, a command-line tool for predicting and analysing SAH-domains in protein sequence datasets. Using Waggawagga-CLI we predicted SAH-domains in 24 datasets from eukaryotes across the tree of life. SAH-domains were predicted in 0.5 to 3.5% of the protein-coding content per species. SAH-domains are particularly present in longer proteins supporting their function as structural building block in multi-domain proteins. In human, SAH-domains are mainly used as alternative building blocks not being present in all transcripts of a gene. Gene ontology analysis showed that yeast proteins with SAH-domains are particular enriched in macromolecular complex subunit organization, cellular component biogenesis and RNA metabolic processes, and that they have a strong nuclear and ribonucleoprotein complex localization and function in ribosome and nucleic acid binding. Human proteins with SAH-domains have roles in all types of RNA processing and cytoskeleton organization, and are predicted to function in RNA binding, protein binding involved in cell and cell-cell adhesion, and cytoskeletal protein binding. Waggawagga-CLI allows the user to adjust the stabilizing and destabilizing contribution of amino acid interactions in i,i+3 and i,i+4 spacings, and provides extensive flexibility for user-designed analyses.

A new role for FBP21 as regulator of Brr2 helicase activity

Splicing of eukaryotic pre-mRNA is carried out by the spliceosome, which assembles stepwise on each splicing substrate. This requires the concerted action of snRNPs and non-snRNP accessory proteins, the functions of which are often not well understood. Of special interest are B complex factors that enter the spliceosome prior to catalytic activation and may alter splicing kinetics and splice site selection. One of these proteins is FBP21, for which we identified several spliceosomal binding partners in a yeast-two-hybrid screen, among them the RNA helicase Brr2. Biochemical and biophysical analyses revealed that an intrinsically disordered region of FBP21 binds to an extended surface of the C-terminal Sec63 unit of Brr2. Additional contacts in the C-terminal helicase cassette are required for allosteric inhibition of Brr2 helicase activity. Furthermore, the direct interaction between FBP21 and the U4/U6 di-snRNA was found to reduce the pool of unwound U4/U6 di-snRNA. Our results suggest FBP21 as a novel key player in the regulation of Brr2.

The spliceosome-associated protein Mfap1 binds to VCP in Drosophila

Posttranscriptional regulation of gene expression contributes to many developmental transitions. Previously, we found that the AAA chaperone Valosin-Containing Protein (VCP) regulates ecdysone-dependent dendrite pruning of Drosophila class IV dendritic arborization (c4da) neurons via an effect on RNA metabolism. In a search for RNA binding proteins associated with VCP, we identified the spliceosome-associated protein Mfap1, a component of the tri-snRNP complex. Mfap1 is a nucleolar protein in neurons and its levels are regulated by VCP. Mfap1 binds to VCP and TDP-43, a disease-associated RNA-binding protein. via distinct regions in its N- and C-terminal halfs. Similar to vcp mutations, Mfap1 overexpression causes c4da neuron dendrite pruning defects and mislocalization of TDP-43 in these cells, but genetic analyses show that Mfap1 is not a crucial VCP target during dendrite pruning. Finally, rescue experiments with a lethal mfap1 mutant show that the VCP binding region is not essential for Mfap1 function, but may act to increase its stability or activity.

Structure of a pre-catalytic spliceosome

Intron removal requires assembly of the spliceosome on precursor mRNA (pre-mRNA) and extensive remodelling to form the spliceosome's catalytic centre. Here we report the cryo-electron microscopy structure of the yeast Saccharomyces cerevisiae pre-catalytic B complex spliceosome at near-atomic resolution. The mobile U2 small nuclear ribonucleoprotein particle (snRNP) associates with U4/U6.U5 tri-snRNP through the U2/U6 helix II and an interface between U4/U6 di-snRNP and the U2 snRNP SF3b-containing domain, which also transiently contacts the helicase Brr2. The 3' region of the U2 snRNP is flexibly attached to the SF3b-containing domain and protrudes over the concave surface of tri-snRNP, where the U1 snRNP may reside before its release from the pre-mRNA 5' splice site. The U6 ACAGAGA sequence forms a hairpin that weakly tethers the 5' splice site. The B complex proteins Prp38, Snu23 and Spp381 bind the Prp8 N-terminal domain and stabilize U6 ACAGAGA stem-pre-mRNA and Brr2-U4 small nuclear RNA interactions. These results provide important insights into the events leading to active site formation.

Human MFAP1 is a cryptic ortholog of the Saccharomyces cerevisiae Spp381 splicing factor

Background

Pre-mRNA splicing involves the stepwise assembly of a pre-catalytic spliceosome, followed by its catalytic activation, splicing catalysis and disassembly. Formation of the pre-catalytic spliceosomal B complex involves the incorporation of the U4/U6.U5 tri-snRNP and of a group of non-snRNP B-specific proteins. While in Saccharomyces cerevisiae the Prp38 and Snu23 proteins are recruited as components of the tri-snRNP, metazoan orthologs of Prp38 and Snu23 associate independently of the tri-snRNP as members of the B-specific proteins. The human spliceosome contains about 80 proteins that lack obvious orthologs in yeast, including most of the B-specific proteins apart from Prp38 and Snu23. Conversely, the tri-snRNP protein Spp381 is one of only five S. cerevisiae splicing factors without a known human ortholog.

Results

Using InParanoid, a state-of-the-art method for ortholog inference between pairs of species, and systematic BLAST searches we identified the human B-specific protein MFAP1 as a putative ortholog of the S. cerevisiae tri-snRNP protein Spp381. Bioinformatics revealed that MFAP1 and Spp381 share characteristic structural features, including intrinsic disorder, an elongated shape, solvent exposure of most residues and a trend to adopt α-helical structures. In vitro binding studies showed that human MFAP1 and yeast Spp381 bind their respective Prp38 proteins via equivalent interfaces and that they cross-interact with the Prp38 proteins of the respective other species. Furthermore, MFAP1 and Spp381 both form higher-order complexes that additionally include Snu23, suggesting that they are parts of equivalent spliceosomal sub-complexes. Finally, similar to yeast Spp381, human MFAP1 partially rescued a growth defect of the temperature-sensitive mutant yeast strain prp38-1.

Conclusions

Human B-specific protein MFAP1 structurally and functionally resembles the yeast tri-snRNP-specific protein Spp381 and thus qualifies as its so far missing ortholog. Our study indicates that the yeast Snu23-Prp38-Spp381 triple complex was evolutionarily reprogrammed from a tri-snRNP-specific module in yeast to the B-specific Snu23-Prp38-MFAP1 module in metazoa, affording higher flexibility in spliceosome assembly and thus, presumably, in splicing regulation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0923-1) contains supplementary material, which is available to authorized users.

Characterization of long and stable de novo single alpha-helix domains provides novel insight into their stability

Naturally-occurring single α-helices (SAHs), are rich in Arg (R), Glu (E) and Lys (K) residues, and stabilized by multiple salt bridges. Understanding how salt bridges promote their stability is challenging as SAHs are long and their sequences highly variable. Thus, we designed and tested simple de novo 98-residue polypeptides containing 7-residue repeats (AEEEXXX, where X is K or R) expected to promote salt-bridge formation between Glu and Lys/Arg. Lys-rich sequences (EK3 (AEEEKKK) and EK2R1 (AEEEKRK)) both form SAHs, of which EK2R1 is more helical and thermo-stable suggesting Arg increases stability. Substituting Lys with Arg (or vice versa) in the naturally-occurring myosin-6 SAH similarly increased (or decreased) its stability. However, Arg-rich de novo sequences (ER3 (AEEERRR) and EK1R2 (AEEEKRR)) aggregated. Combining a PDB analysis with molecular modelling provides a rational explanation, demonstrating that Glu and Arg form salt bridges more commonly, utilize a wider range of rotamer conformations, and are more dynamic than Glu-Lys. This promiscuous nature of Arg helps explain the increased propensity of de novo Arg-rich SAHs to aggregate. Importantly, the specific K:R ratio is likely to be important in determining helical stability in de novo and naturally-occurring polypeptides, giving new insight into how single α-helices are stabilized.