U2AF65 adapts to diverse pre-mRNA splice sites through conformational selection of specific and promiscuous RNA recognition motifs

Exonic splicing code and coordination of divalent metals in proteins

Exonic sequences contain both protein-coding and RNA splicing information but the interplay of the protein and splicing code is complex and poorly understood. Here, we have studied traditional and auxiliary splicing codes of human exons that encode residues coordinating two essential divalent metals at the opposite ends of the Irving-Williams series, a universal order of relative stabilities of metal-organic complexes. We show that exons encoding Zn2+-coordinating amino acids are supported much less by the auxiliary splicing motifs than exons coordinating Ca2+. The handicap of the former is compensated by stronger splice sites and uridine-richer polypyrimidine tracts, except for position -3 relative to 3' splice junctions. However, both Ca2+ and Zn2+ exons exhibit close-to-constitutive splicing in multiple tissues, consistent with their critical importance for metalloprotein function and a relatively small fraction of expendable, alternatively spliced exons. These results indicate that constraints imposed by metal coordination spheres on RNA splicing have been efficiently overcome by the plasticity of exon-intron architecture to ensure adequate metalloprotein expression.

Assessing the Binding Mode of a Splicing Modulator Stimulating Pre-mRNA Binding to the Plastic U2AF2 Splicing Factor

RNA recognition motifs (RRMs) play a pivotal role in RNA metabolism and the regulation of gene expression. Owing to their plasticity and fuzziness, targeting RRM/RNA interfaces with small molecules is a daunting challenge for drug discovery campaigns. The U2AF2 splicing factor, which recognizes the polypyrimidine (polyPy) sequence of premature messenger (pre-m)RNA, exhibits a dynamic architecture consisting of two RRMs joined by a disordered linker. An inhibitor, NSC-194308, was shown to enhance the binding of pre-mRNA to U2AF2, selectively triggering cell death in leukemia cell lines containing spliceosome mutations. The NSC-194308 binding mode remains elusive; yet, unraveling its knowledge may offer intriguing insights for effectively targeting U2AF2 and other flexible protein/protein/RNA interfaces with small molecules. To infer plausible NSC-194308 binding poses to U2AF2, here, we applied and benchmarked the performance of static and dynamic docking approaches, elucidating the molecular basis of the NSC-194308-induced pre-mRNA stabilization on U2AF2. We demonstrate that introducing dynamic effects is mandatory to assess the binding mode of the inhibitors when they target plastic and modular architectures, such as those formed by interacting RRMs. The latter are widespread across RNA binding proteins; therefore, this mechanism may be broadly applicable to discover new therapeutics aimed at selectively modulating the RNA function by targeting protein/protein/RNA interfaces.

Functional analysis of 3'-UTR hairpins supports a two-tiered model for posttranscriptional regulation of MAT2A by METTL16

S-adenosylmethionine (SAM) is the methyl donor for nearly all cellular methylation events, so cells need to carefully control SAM levels. MAT2A encodes the only SAM synthetase expressed in the majority of human cells, and its 3'-UTR has six conserved regulatory hairpins (hp1-6) that can be methylated by the N6-methyladenosine methyltransferase METTL16. Hp1 begins 8 nt from the stop codon, whereas hp2-6 are clustered further downstream (∼800 nt). These hairpins have been proposed to regulate MAT2A mRNA levels in response to intracellular SAM levels by regulating intron detention of the last intron of MAT2A and by modulating the stability of the fully spliced mRNA. However, a dissection of these two posttranscriptional mechanisms has not been previously reported. Using a modular reporter system, we show that hp1 functions primarily when the detained intron is included in the reporter and when that intron has a suboptimal polypyrimidine tract. In contrast, the hp2-6 cluster modulates mRNA stability independent of the detained intron, although hp1 may make a minor contribution to the regulation of decay as well. Taken with previously published reports, these data support a two-tiered model for MAT2A posttranscriptional regulation by METTL16 through its interactions with hp1 and hp2-6. In the upstream tier, hp1 and METTL16 control MAT2A intron detention, whereas the second tier involves METTL16-dependent methylation of hp2-6 to control MAT2A mRNA stability. Thus, cells use a similar set of molecular factors to achieve considerable complexity in the posttranscriptional regulation of SAM homeostasis.

Why to target G-quadruplexes using peptides: Next-generation G4-interacting ligands

Guanine-rich oligonucleotides existing in both DNA and RNA are able to fold into four-stranded DNA secondary structures via Hoogsteen type hydrogen-bonding, where four guanines self-assemble into a square planar arrangement, which, when stacked upon each other, results in the formation of higher-order structures called G-quadruplexes. Their distribution is not random; they are more frequently present at telomeres, proto-oncogenic promoters, introns, 5'- and 3'-untranslated regions, stem cell markers, ribosome binding sites and so forth and are associated with various biological functions, all of which play a pivotal role in various incurable diseases like cancer and cellular ageing. Several studies have suggested that G-quadruplexes could not regulate biological processes by themselves; instead, various proteins take part in this regulation and can be important therapeutic targets. There are certain limitations in using whole G4-protein for therapeutics purpose because of its high manufacturing cost, laborious structure prediction, dynamic nature, unavailability for oral administration due to its degradation in the gut and inefficient penetration to reach the target site because of the large size. Hence, biologically active peptides can be the potential candidates for therapeutic intervention instead of the whole G4-protein complex. In this review, we aimed to clarify the biological roles of G4s, how we can identify them throughout the genome via bioinformatics, the proteins interacting with G4s and how G4-interacting peptide molecules may be the potential next-generation ligands for targeting the G4 motifs located in biologically important regions.

Early Splicing Complexes and Human Disease

Over the last decade, our understanding of spliceosome structure and function has significantly improved, refining the study of the impact of dysregulated splicing on human disease. As a result, targeted splicing therapeutics have been developed, treating various diseases including spinal muscular atrophy and Duchenne muscular dystrophy. These advancements are very promising and emphasize the critical role of proper splicing in maintaining human health. Herein, we provide an overview of the current information on the composition and assembly of early splicing complexes-commitment complex and pre-spliceosome-and their association with human disease.

Monovalent Ionic Atmosphere Modulates the Selection of Suboptimal RNA Sequences by Splicing Factors' RNA Recognition Motifs

The U2AF2 splicing factor is involved in the RNA recognition of the pre-mRNA poly-pyrimidine signaling sequence. This protein contains two RRM domains connected by a flexible linker, which ensure the preferential selection of a poly-uridine sequence over a poly-cytosine one. In this work, all-atom simulations provide insights into the U2AF2 recognition mechanism and on the features underlying its selectivity. Our outcomes show that U2AF2's RNA recognition is driven by cooperative events modulated by RNA-protein and RNA-ion interactions. Stunningly, monovalent ions contribute to mediating the binding of the weakly binding polyC strand, thus contributing to the selection of suboptimal poly-pyrimidine tracts. This finding broadens our understanding of the diverse traits tuning splicing factors' selectivity and adaptability to precisely handle and process diverse pre-mRNA sequences.

All-Atom Simulations Elucidate the Impact of U2AF2 Cancer-Associated Mutations on Pre-mRNA Recognition

The U2AF2 splicing factor, made of two tandem RNA recognition motifs (RRMs) joined by a flexible linker, selects the intronic polypyrimidine sequence of premature mRNA, thus ensuring splicing fidelity. Increasing evidence links mutations of key splicing factors, including U2AF2, to a variety of cancers. Nevertheless, the impact of U2AF2 cancer-associated mutations on polypyrimidine recognition remains unclear. Here, we combined extensive (18 μs-long) all-atom molecular dynamics simulations and dynamical network theory analysis (NWA) of U2AF2, in its wild-type form and in the presence of the six most frequent cancer-associated mutations, bound to a poly-U strand. Our results reveal that the selected mutations affect the pre-mRNA binding at two hot spot regions, irrespectively of where these mutants are placed on the distinct U2AF2 domains. Complementarily, NWA traced the existence of cross-communication pathways, connecting each mutation site to these recognition hot spots, whose strength is altered by the mutations. Our outcomes suggest the existence of a structural/dynamical interplay of the two U2AF2's RRMs underlying the recognition of the polypyrimidine tract and reveal that the cancer-associated mutations affect the polypyrimidine selection by altering the RRMs' cooperativity. This mechanism may be shared by other RNA binding proteins hallmarked, like U2AF2, by multidomain architecture and high plasticity.

Insight into the Tumor Suppression Mechanism from the Structure of Human Polypyrimidine Splicing Factor (PSF/SFPQ) Complexed with a 30mer RNA from Murine Virus-like 30S Transcript-1

Human polypyrimidine-binding splicing factor (PSF/SFPQ) is a tumor suppressor protein that regulates the gene expression of several proto-oncogenes and binds to the 5'-polyuridine negative-sense template (5'-PUN) of some RNA viruses. The activity of PSF is negatively regulated by long-noncoding RNAs, human metastasis associated in lung adenocarcinoma transcript-1 and murine virus-like 30S transcript-1 (VL30-1). PSF is a 707-amino acid protein that has a DNA-binding domain and two RNA recognition motifs (RRMs). Although the structure of the apo-truncated PSF is known, how PSF recognizes RNA remains elusive. Here, we report the 2.8 Å and 3.5 Å resolution crystal structures of a biologically active truncated construct of PSF (sPSF, consisting of residues 214-598) alone and in a complex with a 30mer fragment of VL30-1 RNA, respectively. The structure of the complex reveals how the 30mer RNA is recognized at two U-specific induced-fit binding pockets, located at the previously unrecognized domain-swapped, inter-subunit RRM1 (of the first subunit)-RRM2 (of the second subunit) interfaces that do not exist in the apo structure. Thus, the sPSF dimer appears to have two conformations in solution: one in a low-affinity state for RNA binding, as seen in the apo-structure, and the other in a high-affinity state for RNA binding, as seen in the sPSF-RNA complex. PSF undergoes an all or nothing transition between having two or no RNA-binding pockets. We predict that the RNA binds with a high degree of positive cooperativity. These structures provide an insight into a new regulatory mechanism that is likely involved in promoting malignancies and other human diseases.

Pre-mRNA splicing factor U2AF2 recognizes distinct conformations of nucleotide variants at the center of the pre-mRNA splice site signal

The essential pre-mRNA splicing factor U2AF2 (also called U2AF65) identifies polypyrimidine (Py) tract signals of nascent transcripts, despite length and sequence variations. Previous studies have shown that the U2AF2 RNA recognition motifs (RRM1 and RRM2) preferentially bind uridine-rich RNAs. Nonetheless, the specificity of the RRM1/RRM2 interface for the central Py tract nucleotide has yet to be investigated. We addressed this question by determining crystal structures of U2AF2 bound to a cytidine, guanosine, or adenosine at the central position of the Py tract, and compared U2AF2-bound uridine structures. Local movements of the RNA site accommodated the different nucleotides, whereas the polypeptide backbone remained similar among the structures. Accordingly, molecular dynamics simulations revealed flexible conformations of the central, U2AF2-bound nucleotide. The RNA binding affinities and splicing efficiencies of structure-guided mutants demonstrated that U2AF2 tolerates nucleotide substitutions at the central position of the Py tract. Moreover, enhanced UV-crosslinking and immunoprecipitation of endogenous U2AF2 in human erythroleukemia cells showed uridine-sensitive binding sites, with lower sequence conservation at the central nucleotide positions of otherwise uridine-rich, U2AF2-bound splice sites. Altogether, these results highlight the importance of RNA flexibility for protein recognition and take a step towards relating splice site motifs to pre-mRNA splicing efficiencies.

Somatic Mutations in Core Spliceosome Components Promote Tumorigenesis and Generate an Exploitable Vulnerability in Human Cancer

Simple Summary

High throughput exome sequencing approaches have disclosed recurrent cancer-associated mutations in spliceosomal components, which drive aberrant pre-mRNA processing events and support the tumor phenotype. At the same time, mutations in spliceosome genes and aberrant splicing regulation establish a selective vulnerability of cancer cells to splicing-targeting approaches, which could be exploited therapeutically. It is conceivable that a better understanding of the mechanisms and roles of abnormal splicing in tumor metabolism will facilitate the development of a novel generation of tumor-targeting drugs. In this review, we describe recent advances in the elucidation of the biological impact and biochemical effects of somatic mutations in core spliceosome components on splicing choices and their associated targetable vulnerabilities.

Abstract

Alternative pre-mRNA processing enables the production of distinct mRNA and protein isoforms from a single gene, thus greatly expanding the coding potential of eukaryotic genomes and fine-tuning gene expression programs. Splicing is carried out by the spliceosome, a complex molecular machinery which assembles step-wise on mRNA precursors in the nucleus of eukaryotic cells. In the last decade, exome sequencing technologies have allowed the identification of point mutations in genes encoding splicing factors as a recurrent hallmark of human cancers, with higher incidence in hematological malignancies. These mutations lead to production of splicing factors that reduce the fidelity of the splicing process and yield splicing variants that are often advantageous for cancer cells. However, at the same time, these mutations increase the sensitivity of transformed cells to splicing inhibitors, thus offering a therapeutic opportunity for novel targeted strategies. Herein, we review the recent literature documenting cancer-associated mutations in components of the early spliceosome complex and discuss novel therapeutic strategies based on small-molecule spliceosome inhibitors that exhibit strong anti-tumor effects, particularly against cancer cells harboring mutations in spliceosomal components.

Splicing factor mutations in hematologic malignancies

Mutations in genes encoding RNA splicing factors were discovered nearly ten years ago and are now understood to be amongst the most recurrent genetic abnormalities in patients with all forms of myeloid neoplasms and several types of lymphoproliferative disorders as well as subjects with clonal hematopoiesis. These discoveries implicate aberrant RNA splicing, the process by which precursor RNA is converted into mature messenger RNA, in the development of clonal hematopoietic conditions. Both the protein as well as the RNA components of the splicing machinery are affected by mutations at highly specific residues and a number of these mutations alter splicing in a manner distinct from loss of function. Importantly, cells bearing these mutations have now been shown to generate mRNA species with novel aberrant sequences, some of which may be critical to disease pathogenesis and/or novel targets for therapy. These findings have opened new avenues of research to understand biological pathways disrupted by altered splicing. In parallel, multiple studies have revealed that cells bearing change-of-function mutation in splicing factors are preferentially sensitized to any further genetic or chemical perturbations of the splicing machinery. These discoveries are now being pursued in several early phase clinical trials using molecules with diverse mechanisms of action. Here we review the molecular effects of splicing factor mutations on splicing, mechanisms by which these mutations drive clonal transformation of hematopoietic cells, and the development of new therapeutics targeting these genetic subsets of hematopoietic malignancies.

A Combined Experimental and Computational Study of Halogen and Hydrogen Bonding in Molecular Salts of 5-Bromocytosine

Although natural or artificial modified pyrimidine nucleobases represent important molecules with valuable properties as constituents of DNA and RNA, no systematic analyses of the structural aspects of bromo derivatives of cytosine have appeared so far in the literature. In view of the biochemical and pharmaceutical relevance of these compounds, six different crystals containing proton-transfer derivatives of 5-bromocytosine are prepared and analyzed in the solid-state by single crystal X-ray diffraction. All six compounds are organic salts, with proton transfer occurring to the N_imino atom of the pyridine ring. Experimental results are then complemented with Hirshfeld surface analysis to quantitively evaluate the contribution of different intermolecular interactions in the crystal packing. Furthermore, theoretical calculations, based on different arrangements of molecules extracted from the crystal structure determinations, are carried out to analyze the formation mechanism of halogen bonds (XBs) in these compounds and provide insights into the nature and strength of the observed interactions. The results show that the supramolecular architectures of the six molecular salts involve extensive classical intermolecular hydrogen bonds. However, in all but one proton-transfer adducts, weak to moderate XBs are revealed by C-Br^…O short contacts between the bromine atom in the fifth position, which acts as XB donor (electron acceptor). Moreover, the lone pair electrons of the oxygen atom of adjacent pyrimidine nucleobases and/or counterions or water molecules, which acts as XB acceptor (electron donor).

SAM homeostasis is regulated by CFIm-mediated splicing of MAT2A

S-adenosylmethionine (SAM) is the methyl donor for nearly all cellular methylation events. Cells regulate intracellular SAM levels through intron detention of MAT2A, the only SAM synthetase expressed in most cells. The N6-adenosine methyltransferase METTL16 promotes splicing of the MAT2A detained intron by an unknown mechanism. Using an unbiased CRISPR knock-out screen, we identified CFIm25 (NUDT21) as a regulator of MAT2A intron detention and intracellular SAM levels. CFIm25 is a component of the cleavage factor Im (CFIm) complex that regulates poly(A) site selection, but we show it promotes MAT2A splicing independent of poly(A) site selection. CFIm25-mediated MAT2A splicing induction requires the RS domains of its binding partners, CFIm68 and CFIm59 as well as binding sites in the detained intron and 3´ UTR. These studies uncover mechanisms that regulate MAT2A intron detention and reveal a previously undescribed role for CFIm in splicing and SAM metabolism.

Inducible nuclear import by TetR aptamer-controlled 3' splice site selection

Correct cellular localization is essential for the function of many eukaryotic proteins and hence cell physiology. Here, we present a synthetic genetic device that allows the control of nuclear and cytosolic localization based on controlled alternative splicing in human cells. The device is based on the fact that an alternative 3' splice site is located within a TetR aptamer that in turn is positioned between the branch point and the canonical splice site. The novel splice site is only recognized when the TetR repressor is bound. Addition of doxycycline prevents TetR aptamer binding and leads to recognition of the canonical 3' splice site. It is thus possible to produce two independent splice isoforms. Since the terminal loop of the aptamer may be replaced with any sequence of choice, one of the two isoforms may be extended by the respective sequence of choice depending on the presence of doxycycline. In a proof-of-concept study, we fused a nuclear localization sequence to a cytosolic target protein, thus directing the protein into the nucleus. However, the system is not limited to the control of nuclear localization. In principle, any target sequence can be integrated into the aptamer, allowing not only the production of a variety of different isoforms on demand, but also to study the function of mislocalized proteins. Moreover, it also provides a valuable tool for investigating the mechanism of alternative splicing in human cells.

A splice site-sensing conformational switch in U2AF2 is modulated by U2AF1 and its recurrent myelodysplasia-associated mutation

An essential heterodimer of the U2AF1 and U2AF2 pre-mRNA splicing factors nucleates spliceosome assembly at polypyrimidine (Py) signals preceding the major class of 3′ splice sites. U2AF1 frequently acquires an S34F-encoding mutation among patients with myelodysplastic syndromes (MDS). The influence of the U2AF1 subunit and its S34F mutation on the U2AF2 conformations remains unknown. Here, we employ single molecule Förster resonance energy transfer (FRET) to determine the influence of wild-type or S34F-substituted U2AF1 on the conformational dynamics of U2AF2 and its splice site RNA complexes. In the absence of RNA, the U2AF1 subunit stabilizes a high FRET value, which by structure-guided mutagenesis corresponds to a closed conformation of the tandem U2AF2 RNA recognition motifs (RRMs). When the U2AF heterodimer is bound to a strong, uridine-rich splice site, U2AF2 switches to a lower FRET value characteristic of an open, side-by-side arrangement of the RRMs. Remarkably, the U2AF heterodimer binds weak, uridine-poor Py tracts as a mixture of closed and open U2AF2 conformations, which are modulated by the S34F mutation. Shifts between open and closed U2AF2 may underlie U2AF1-dependent splicing of degenerate Py tracts and contribute to a subset of S34F-dysregulated splicing events in MDS patients.

Quantifying Intramolecular Halogen Bonds in Nucleic Acids: A Combined Protein Data Bank and Theoretical Study

In this study, we report experimental (Protein Data Bank (PDB) search) and theoretical (RI-MP2/def2-TZVP level of theory) evidence of the nature, stability, and directionality properties of intramolecular halogen bonding interactions (HaBs) between 5-bromo/5-iodoracil bases and backbone phosphate groups in nucleic acids (NAs). A PDB survey revealed relevant examples where intramolecular HaBs are undertaken and serve as a structural source of stability in RNA and DNA molecules. In order to develop suitable energy predictors, we started this investigation by calculating the interaction energy values and both the potential V(r) and kinetic G(r) energy densities (using Bader's "atoms in molecules" theory) of several halogen bond complexes involving 5-bromo/5-iodoracil molecules and biologically relevant electron donors. Once the energy predictors based on V(r)/G(r) energy densities were developed, we analyzed the HaBs observed in the biological examples retrieved from the PDB search in order to estimate the strength of the interaction. As far as our knowledge extends, intramolecular halogen bonds in NAs have not been previously quantified in the literature using this methodology and may be of great importance in understanding their structural properties as well as in the construction of molecular materials with DNA and other biological macromolecules.

How does temperature affect splicing events? Isoform switching of splicing factors regulates splicing of LATE ELONGATED HYPOCOTYL (LHY)

One of the ways in which plants can respond to temperature is via alternative splicing (AS). Previous work showed that temperature changes affected the splicing of several circadian clock gene transcripts. Here, we investigated the role of RNA-binding splicing factors (SFs) in temperature-sensitive AS of the clock gene LATE ELONGATED HYPOCOTYL (LHY). We characterized, in wild type plants, temperature-associated isoform switching and expression patterns for SF transcripts from a high-resolution temperature and time series RNA-seq experiment. In addition, we employed quantitative RT-PCR of SF mutant plants to explore the role of the SFs in cooling-associated AS of LHY. We show that the splicing and expression of several SFs responds sufficiently, rapidly, and sensitively to temperature changes to contribute to the splicing of the 5'UTR of LHY. Moreover, the choice of splice site in LHY was altered in some SF mutants. The splicing of the 5'UTR region of LHY has characteristics of a molecular thermostat, where the ratio of transcript isoforms is sensitive to temperature changes as modest as 2 °C and is scalable over a wide dynamic range of temperature. Our work provides novel insight into SF-mediated coupling of the perception of temperature to post-transcriptional regulation of the clock.

H3B-8800, an orally available small-molecule splicing modulator, induces lethality in spliceosome-mutant cancers

Genomic analyses of cancer have identified recurrent point mutations in the RNA splicing factor-encoding genes SF3B1, U2AF1, and SRSF2 that confer an alteration of function. Cancer cells bearing these mutations are preferentially dependent on wild-type (WT) spliceosome function, but clinically relevant means to therapeutically target the spliceosome do not currently exist. Here we describe an orally available modulator of the SF3b complex, H3B-8800, which potently and preferentially kills spliceosome-mutant epithelial and hematologic tumor cells. These killing effects of H3B-8800 are due to its direct interaction with the SF3b complex, as evidenced by loss of H3B-8800 activity in drug-resistant cells bearing mutations in genes encoding SF3b components. Although H3B-8800 modulates WT and mutant spliceosome activity, the preferential killing of spliceosome-mutant cells is due to retention of short, GC-rich introns, which are enriched for genes encoding spliceosome components. These data demonstrate the therapeutic potential of splicing modulation in spliceosome-mutant cancers.

Stoichiometries of U2AF35, U2AF65 and U2 snRNP reveal new early spliceosome assembly pathways

The selection of 3΄ splice sites (3΄ss) is an essential early step in mammalian RNA splicing reactions, but the processes involved are unknown. We have used single molecule methods to test whether the major components implicated in selection, the proteins U2AF35 and U2AF65 and the U2 snRNP, are able to recognize alternative candidate sites or are restricted to one pre-specified site. In the presence of adenosine triphosphate (ATP), all three components bind in a 1:1 stoichiometry with a 3΄ss. Pre-mRNA molecules with two alternative 3΄ss can be bound concurrently by two molecules of U2AF or two U2 snRNPs, so none of the components are restricted. However, concurrent occupancy inhibits splicing. Stoichiometric binding requires conditions consistent with coalescence of the 5΄ and 3΄ sites in a complex (I, initial), but if this cannot form the components show unrestricted and stochastic association. In the absence of ATP, when complex E forms, U2 snRNP association is unrestricted. However, if protein dephosphorylation is prevented, an I-like complex forms with stoichiometric association of U2 snRNPs and the U2 snRNA is base-paired to the pre-mRNA. Complex I differs from complex A in that the formation of complex A is associated with the loss of U2AF65 and 35.

Cancer-Associated Mutations Mapped on High-Resolution Structures of the U2AF2 RNA Recognition Motifs

Acquired point mutations of pre-mRNA splicing factors recur among cancers, leukemias, and related neoplasms. Several studies have established that somatic mutations of a U2AF1 subunit, which normally recognizes 3' splice site junctions, recur among myelodysplastic syndromes. The U2AF2 splicing factor recognizes polypyrimidine signals that precede most 3' splice sites as a heterodimer with U2AF1. In contrast with those of the well-studied U2AF1 subunit, descriptions of cancer-relevant U2AF2 mutations and their structural relationships are lacking. Here, we survey databases of cancer-associated mutations and identify recurring missense mutations in the U2AF2 gene. We determine ultra-high-resolution structures of the U2AF2 RNA recognition motifs (RRM1 and RRM2) at 1.1 Å resolution and map the structural locations of the mutated U2AF2 residues. Comparison with prior, lower-resolution structures of the tandem U2AF2 RRMs in the RNA-bound and apo states reveals clusters of cancer-associated mutations at the U2AF2 RRM-RNA or apo-RRM1-RRM2 interfaces. Although the role of U2AF2 mutations in malignant transformation remains uncertain, our results show that cancer-associated mutations correlate with functionally important surfaces of the U2AF2 splicing factor.

Structural study of the Fox-1 RRM protein hydration reveals a role for key water molecules in RRM-RNA recognition

The Fox-1 RNA recognition motif (RRM) domain is an important member of the RRM protein family. We report a 1.8 Å X-ray structure of the free Fox-1 containing six distinct monomers. We use this and the nuclear magnetic resonance (NMR) structure of the Fox-1 protein/RNA complex for molecular dynamics (MD) analyses of the structured hydration. The individual monomers of the X-ray structure show diverse hydration patterns, however, MD excellently reproduces the most occupied hydration sites. Simulations of the protein/RNA complex show hydration consistent with the isolated protein complemented by hydration sites specific to the protein/RNA interface. MD predicts intricate hydration sites with water-binding times extending up to hundreds of nanoseconds. We characterize two of them using NMR spectroscopy, RNA binding with switchSENSE and free-energy calculations of mutant proteins. Both hydration sites are experimentally confirmed and their abolishment reduces the binding free-energy. A quantitative agreement between theory and experiment is achieved for the S155A substitution but not for the S122A mutant. The S155 hydration site is evolutionarily conserved within the RRM domains. In conclusion, MD is an effective tool for predicting and interpreting the hydration patterns of protein/RNA complexes. Hydration is not easily detectable in NMR experiments but can affect stability of protein/RNA complexes.

Integrated structural biology to unravel molecular mechanisms of protein-RNA recognition

Recent advances in RNA sequencing technologies have greatly expanded our knowledge of the RNA landscape in cells, often with spatiotemporal resolution. These techniques identified many new (often non-coding) RNA molecules. Large-scale studies have also discovered novel RNA binding proteins (RBPs), which exhibit single or multiple RNA binding domains (RBDs) for recognition of specific sequence or structured motifs in RNA. Starting from these large-scale approaches it is crucial to unravel the molecular principles of protein-RNA recognition in ribonucleoprotein complexes (RNPs) to understand the underlying mechanisms of gene regulation. Structural biology and biophysical studies at highest possible resolution are key to elucidate molecular mechanisms of RNA recognition by RBPs and how conformational dynamics, weak interactions and cooperative binding contribute to the formation of specific, context-dependent RNPs. While large compact RNPs can be well studied by X-ray crystallography and cryo-EM, analysis of dynamics and weak interaction necessitates the use of solution methods to capture these properties. Here, we illustrate methods to study the structure and conformational dynamics of protein-RNA complexes in solution starting from the identification of interaction partners in a given RNP. Biophysical and biochemical techniques support the characterization of a protein-RNA complex and identify regions relevant in structural analysis. Nuclear magnetic resonance (NMR) is a powerful tool to gain information on folding, stability and dynamics of RNAs and characterize RNPs in solution. It provides crucial information that is complementary to the static pictures derived from other techniques. NMR can be readily combined with other solution techniques, such as small angle X-ray and/or neutron scattering (SAXS/SANS), electron paramagnetic resonance (EPR), and Förster resonance energy transfer (FRET), which provide information about overall shapes, internal domain arrangements and dynamics. Principles of protein-RNA recognition and current approaches are reviewed and illustrated with recent studies.

Protein Ensembles: How Does Nature Harness Thermodynamic Fluctuations for Life? The Diverse Functional Roles of Conformational Ensembles in the Cell

All soluble proteins populate conformational ensembles that together constitute the native state. Their fluctuations in water are intrinsic thermodynamic phenomena, and the distributions of the states on the energy landscape are determined by statistical thermodynamics; however, they are optimized to perform their biological functions. In this review we briefly describe advances in free energy landscape studies of protein conformational ensembles. Experimental (nuclear magnetic resonance, small-angle X-ray scattering, single-molecule spectroscopy, and cryo-electron microscopy) and computational (replica-exchange molecular dynamics, metadynamics, and Markov state models) approaches have made great progress in recent years. These address the challenging characterization of the highly flexible and heterogeneous protein ensembles. We focus on structural aspects of protein conformational distributions, from collective motions of single- and multi-domain proteins, intrinsically disordered proteins, to multiprotein complexes. Importantly, we highlight recent studies that illustrate functional adjustment of protein conformational ensembles in the crowded cellular environment. We center on the role of the ensemble in recognition of small- and macro-molecules (protein and RNA/DNA) and emphasize emerging concepts of protein dynamics in enzyme catalysis. Overall, protein ensembles link fundamental physicochemical principles and protein behavior and the cellular network and its regulation.

An extended U2AF(65)-RNA-binding domain recognizes the 3' splice site signal

How the essential pre-mRNA splicing factor U2AF⁶⁵ recognizes the polypyrimidine (Py) signals of the major class of 3′ splice sites in human gene transcripts remains incompletely understood. We determined four structures of an extended U2AF⁶⁵–RNA-binding domain bound to Py-tract oligonucleotides at resolutions between 2.0 and 1.5 Å. These structures together with RNA binding and splicing assays reveal unforeseen roles for U2AF⁶⁵ inter-domain residues in recognizing a contiguous, nine-nucleotide Py tract. The U2AF⁶⁵ linker residues between the dual RNA recognition motifs (RRMs) recognize the central nucleotide, whereas the N- and C-terminal RRM extensions recognize the 3′ terminus and third nucleotide. Single-molecule FRET experiments suggest that conformational selection and induced fit of the U2AF⁶⁵ RRMs are complementary mechanisms for Py-tract association. Altogether, these results advance the mechanistic understanding of molecular recognition for a major class of splice site signals.

The pre-mRNA splicing factor U2AF65 recognizes 3′ splice sites in human gene transcripts, but the details are not fully understood. Here, the authors report U2AF65 structures and single molecule FRET that reveal mechanistic insights into splice site recognition.

The 3.8 Å structure of the U4/U6.U5 tri-snRNP: Insights into spliceosome assembly and catalysis

Splicing of precursor messenger RNA is accomplished by a dynamic megacomplex known as the spliceosome. Assembly of a functional spliceosome requires a preassembled U4/U6.U5 tri-snRNP complex, which comprises the U5 small nuclear ribonucleoprotein (snRNP), the U4 and U6 small nuclear RNA (snRNA) duplex, and a number of protein factors. Here we report the three-dimensional structure of a Saccharomyces cerevisiae U4/U6.U5 tri-snRNP at an overall resolution of 3.8 angstroms by single-particle electron cryomicroscopy. The local resolution for the core regions of the tri-snRNP reaches 3.0 to 3.5 angstroms, allowing construction of a refined atomic model. Our structure contains U5 snRNA, the extensively base-paired U4/U6 snRNA, and 30 proteins including Prp8 and Snu114, which amount to 8495 amino acids and 263 nucleotides with a combined molecular mass of ~1 megadalton. The catalytic nucleotide U80 from U6 snRNA exists in an inactive conformation, stabilized by its base-pairing interactions with U4 snRNA and protected by Prp3. Pre-messenger RNA is bound in the tri-snRNP through base-pairing interactions with U6 snRNA and loop I of U5 snRNA. This structure, together with that of the spliceosome, reveals the molecular choreography of the snRNAs in the activation process of the spliceosomal ribozyme.

The Spliceosome: The Ultimate RNA Chaperone and Sculptor

The spliceosome, one of the most complex machineries of eukaryotic cells, removes intronic sequences from primary transcripts to generate functional messenger and long noncoding RNAs (lncRNA). Genetic, biochemical, and structural data reveal that the spliceosome is an RNA-based enzyme. Striking mechanistic and structural similarities strongly argue that pre-mRNA introns originated from self-catalytic group II ribozymes. However, in the spliceosome, protein components organize and activate the catalytic-site RNAs, and recognize and pair together splice sites at intron boundaries. The spliceosome is a dynamic, reversible, and flexible machine that chaperones small nuclear (sn) RNAs and a variety of pre-mRNA sequences into conformations that enable intron removal. This malleability likely contributes to the regulation of alternative splicing, a prevalent process contributing to cell differentiation, homeostasis, and disease.

Correlated Evolution of Nucleotide Positions within Splice Sites in Mammals

Splice sites (SSs)--short nucleotide sequences flanking introns--are under selection for spliceosome binding, and adhere to consensus sequences. However, non-consensus nucleotides, many of which probably reduce SS performance, are frequent. Little is known about the mechanisms maintaining such apparently suboptimal SSs. Here, we study the correlations between strengths of nucleotides occupying different positions of the same SS. Such correlations may arise due to epistatic interactions between positions (i.e., a situation when the fitness effect of a nucleotide in one position depends on the nucleotide in another position), their evolutionary history, or to other reasons. Within both the intronic and the exonic parts of donor SSs, nucleotides that increase (decrease) SS strength tend to co-occur with other nucleotides increasing (respectively, decreasing) it, consistent with positive epistasis. Between the intronic and exonic parts of donor SSs, the correlations of nucleotide strengths tend to be negative, consistent with negative epistasis. In the course of evolution, substitutions at a donor SS tend to decrease the strength of its exonic part, and either increase or do not change the strength of its intronic part. In acceptor SSs, the situation is more complicated; the correlations between adjacent positions appear to be driven mainly by avoidance of the AG dinucleotide which may cause aberrant splicing. In summary, both the content and the evolution of SSs is shaped by a complex network of interdependences between adjacent nucleotides that respond to a range of sometimes conflicting selective constraints.

Evidence for cooperative tandem binding of hnRNP C RRMs in mRNA processing

The human hnRNP C is a ubiquitous cellular protein involved in mRNA maturation. Recently, we have shown that this protein specifically recognizes uridine (U) pentamers through its single RNA recognition motif (RRM). However, a large fraction of natural RNA targets of hnRNP C consists of much longer contiguous uridine stretches. To understand how these extended sites are recognized, we studied the binding of the RRM to U-tracts of 8-11 bases. In vivo investigation of internal translation activation of unr (upstream of N-ras) mRNA indicates that the conservation of the entire hnRNP C binding site, UC(U)8, is required for hnRNP C-dependent IRES activation. The assays further suggest a synergistic interplay between hnRNP C monomers, dependent on the protein's ability to oligomerize. In vitro spectroscopic and thermodynamic analyses show that isolated RRMs bind to (U)11 oligomers as dimers. Structural modeling of a ternary double-RRM/RNA complex indicates additionally that two RRM copies can be accommodated on the canonical sequence UC(U)8. The proposed tandem RRM binding is in very good agreement with the transcriptome-wide recognition of extended U-tracts by full-length hnRNP C, which displays a cross-linking pattern consistent with a positively cooperative RRM dimer binding model.

Retention and splicing complex (RES) - the importance of cooperativity

One of the great challenges to structural biologists lies in explaining the complexities of the spliceosome – a ribosome-sized molecular machine that is assembled in a step-wise manner and is responsible for pre-mRNA splicing. The spliceosome is both fascinating and difficult to work with, because of its dynamic nature. At each discrete step of splicing tens of proteins come and go orchestrating the functional transition through massive structural rearrangements. The retention and splicing complex (RES) is an important splicing factor interacting with pre-mRNA at the onset of the first transesterification reaction. RES is a specific splicing factor for a number of genes and is important for controlling pre-mRNA retention in the nucleus. RES is a 71 kDa heterotrimer composed of the 3 proteins Pml1p, Bud13p and Snu17p. We solved the 3-dimensional structure of the core of the RES complex as well as the 2 dimers, Snu17p-Bud13p and Snu17p-Pml1p. Further biophysical analysis revealed an astounding cooperativity that governs the assembly of this trimeric complex as well as its interaction with pre-mRNA. The more than 100-fold cooperativity originates from the progressive rigidification of Snu17p upon coupled binding-and-folding of protein regions, which are disordered in the unbound state. Our work highlights the role of cooperativity in the spliceosome and poses new questions about the structure and assembly of the spliceosome.

Integration of the proteome and transcriptome reveals multiple levels of gene regulation in the rice dl2 mutant

Leaf vascular system differentiation and venation patterns play a key role in transporting nutrients and maintaining the plant shape, which is an important agronomic trait for improving photosynthetic efficiency. However, there is little knowledge about the regulation of leaf vascular specification and development. Here we utilized the rice midribless mutant (dl2) to investigate the molecular changes in transcriptome and proteome profiles during leaf vascular specification and differentiation. Using isobaric tags for relative and absolute quantification (iTRAQ) with digital gene expression (DGE) techniques, a nearly complete catalog of expressed protein and mRNA was acquired. From the catalog, we reliably identified 3172 proteins and 9,865,230 tags mapped to genes, and subsets of 141 proteins and 98 mRNAs, which were differentially expressed between the dl2 mutant and wild type. The correlation analysis between the abundance of differentially expressed mRNA and DEPs (differentially expressed proteins) revealed numerous discordant changes in mRNA/protein pairs and only a modest correlation was observed, indicative of divergent regulation of transcription and translational processes. The DEPs were analyzed for their involvement in biological processes and metabolic pathways. Up- or down- regulation of some key proteins confirmed that the physiological process of vascular differentiation is an active process. These key proteins included those not previously reported to be associated with vascular differentiation processes, and included proteins that are involved in the spliceosome pathway. Together, our results show that the developmental and physiological process of the leaf vascular system is a thoroughly regulated and complicated process and this work has identified potential targets for genetic modification that could be used to regulate the development of the leaf vasculature.

One, Two, Three, Four! How Multiple RRMs Read the Genome Sequence

RRM-containing proteins are involved in most of the RNA metabolism steps. Their functions are closely related to their mode of RNA recognition, which has been studied by structural biologists for more than 20 years. In this chapter, we report on high-resolution structures of single and multi RRM-RNA complexes to explain the numerous strategies used by these domains to interact specifically with a large repertoire of RNA sequences. We show that multiple variations of their canonical fold can be used to adapt to different single-stranded sequences with a large range of affinities. Furthermore, we describe the consequences on RNA binding of the different structural arrangements found in tandem RRMs and higher order RNPs. Importantly, these structures also reveal with very high accuracy the RNA motifs bound specifically by RRM-containing proteins, which correspond very often to consensus sequences identified with genome-wide approaches. Finally, we show how structural and cellular biology can benefit from each other and pave a way for understanding, defining, and predicting a code of RNA recognition by the RRMs.

Structure-guided U2AF65 variant improves recognition and splicing of a defective pre-mRNA

Purine interruptions of polypyrimidine (Py) tract splice site signals contribute to human genetic diseases. The essential splicing factor U2AF(65) normally recognizes a Py tract consensus sequence preceding the major class of 3' splice sites. We found that neurofibromatosis- or retinitis pigmentosa-causing mutations in the 5' regions of Py tracts severely reduce U2AF(65) affinity. Conversely, we identified a preferred binding site of U2AF(65) for purine substitutions in the 3' regions of Py tracts. Based on a comparison of new U2AF(65) structures bound to either A- or G-containing Py tracts with previously identified pyrimidine-containing structures, we expected to find that a D231V amino acid change in U2AF(65) would specify U over other nucleotides. We found that the crystal structure of the U2AF(65)-D231V variant confirms favorable packing between the engineered valine and a target uracil base. The D231V amino acid change restores U2AF(65) affinity for two mutated splice sites that cause human genetic diseases and successfully promotes splicing of a defective retinitis pigmentosa-causing transcript. We conclude that reduced U2AF(65) binding is a molecular consequence of disease-relevant mutations, and that a structure-guided U2AF(65) variant is capable of manipulating gene expression in eukaryotic cells.

Exon first nucleotide mutations in splicing: evaluation of in silico prediction tools

Mutations in the first nucleotide of exons (E⁺¹) mostly affect pre-mRNA splicing when found in AG-dependent 3′ splice sites, whereas AG-independent splice sites are more resistant. The AG-dependency, however, may be difficult to assess just from primary sequence data as it depends on the quality of the polypyrimidine tract. For this reason, in silico prediction tools are commonly used to score 3′ splice sites. In this study, we have assessed the ability of sequence features and in silico prediction tools to discriminate between the splicing-affecting and non-affecting E⁺¹ variants. For this purpose, we newly tested 16 substitutions in vitro and derived other variants from literature. Surprisingly, we found that in the presence of the substituting nucleotide, the quality of the polypyrimidine tract alone was not conclusive about its splicing fate. Rather, it was the identity of the substituting nucleotide that markedly influenced it. Among the computational tools tested, the best performance was achieved using the Maximum Entropy Model and Position-Specific Scoring Matrix. As a result of this study, we have now established preliminary discriminative cut-off values showing sensitivity up to 95% and specificity up to 90%. This is expected to improve our ability to detect splicing-affecting variants in a clinical genetic setting.

Multiple conformational selection and induced fit events take place in allosteric propagation

The fact that we observe a single conformational selection event during binding does not necessarily mean that only a single conformational selection event takes place, even though this is the common assumption. Here we suggest that conformational selection takes place not once in a given binding/allosteric event, but at every step along the allosteric pathway. This view generalizes conformational selection and makes it applicable also to other allosteric events, such as post-translational modifications (PTMs) and photon absorption. Similar to binding, at each step along a propagation pathway, conformational selection is coupled with induced fit which optimizes the interactions. Thus, as in binding, the allosteric effects induced by PTMs and light relate not only to population shift; but to conformational selection as well. Conformational selection and population shift take place conjointly.

Coregulation of HIV-1 dependency factors in individuals heterozygous to the CCR5-delta32 deletion

Background

CCR5-delta32 heterozygous individuals are susceptible to HIV-1. However, it is not clear if there is a relevant protective effect against transmission and a beneficial effect in terms of HIV progression which cannot be attributed to CCR5 surface density alone. Therefore we investigated HIV-1 dependency factors (HDF) which might be differently regulated in CCR5 wild type (WT) and CCR5-delta32 heterozygous individuals.

Methods

We examined CD34+ hematopoietic progenitor cells derived from bone marrow samples from 19 healthy volunteers, 12 individuals with CCR5 WT and 7 with heterozygous CCR5-delta32 deletion. Samples were analyzed using a global gene expression oligonucleotide microarray (HG-U133plus 2.0, Affymetrix Inc.).

Results

A total of 205 genes were found with altered expression (3fold difference, present call rate of 75%, p < 0.05) and 7 of these had a connection to HIV-1 pathogenesis. In 4 genes: TOP1, CXCR2, SREBF2, and TAP we found a different regulation which was consistent with a supposed beneficial effect for CCR5-delta32 heterozygotes.

Conclusion

The CCR5-delta32 deletion is associated with other HDFs in HIV-1 pathogenesis as a possible explanation for beneficial effects regarding the deletion leading to a variant expression profile in heterozygous carriers of this mutation.