A Broad range of conformations contribute to the solution ensemble of the essential splicing factor U2AF(65)

Mapping the Constrained Coding Regions in the Human Genome to Their Corresponding Proteins

Constrained Coding Regions (CCRs) in the human genome have been derived from DNA sequencing data of large cohorts of healthy control populations, available in the Genome Aggregation Database (gnomAD) [1]. They identify regions depleted of protein-changing variants and thus identify segments of the genome that have been constrained during human evolution. By mapping these DNA-defined regions from genomic coordinates onto the corresponding protein positions and combining this information with protein annotations, we have explored the distribution of CCRs and compared their co-occurrence with different protein functional features, previously annotated at the amino acid level in public databases. As expected, our results reveal that functional amino acids involved in interactions with DNA/RNA, protein-protein contacts and catalytic sites are the protein features most likely to be highly constrained for variation in the control population. More surprisingly, we also found that linear motifs, linear interacting peptides (LIPs), disorder-order transitions upon binding with other protein partners and liquid-liquid phase separating (LLPS) regions are also strongly associated with high constraint for variability. We also compared intra-species constraints in the human CCRs with inter-species conservation and functional residues to explore how such CCRs may contribute to the analysis of protein variants. As has been previously observed, CCRs are only weakly correlated with conservation, suggesting that intraspecies constraints complement interspecies conservation and can provide more information to interpret variant effects.

Representative cancer-associated U2AF2 mutations alter RNA interactions and splicing

High-throughput sequencing of hematologic malignancies and other cancers has revealed recurrent mis-sense mutations of genes encoding pre-mRNA splicing factors. The essential splicing factor U2AF2 recognizes a polypyrimidine-tract splice-site signal and initiates spliceosome assembly. Here, we investigate representative, acquired U2AF2 mutations, namely N196K or G301D amino acid substitutions associated with leukemia or solid tumors, respectively. We determined crystal structures of the wild-type (WT) compared with N196K- or G301D-substituted U2AF2 proteins, each bound to a prototypical AdML polypyrimidine tract, at 1.5, 1.4, or 1.7 Å resolutions. The N196K residue appears to stabilize the open conformation of U2AF2 with an inter-RNA recognition motif hydrogen bond, in agreement with an increased apparent RNA-binding affinity of the N196K-substituted protein. The G301D residue remains in a similar position as the WT residue, where unfavorable proximity to the RNA phosphodiester could explain the decreased RNA-binding affinity of the G301D-substituted protein. We found that expression of the G301D-substituted U2AF2 protein reduces splicing of a minigene transcript carrying prototypical splice sites. We further show that expression of either N196K- or G301D-substituted U2AF2 can subtly alter splicing of representative endogenous transcripts, despite the presence of endogenous, WT U2AF2 such as would be present in cancer cells. Altogether, our results demonstrate that acquired U2AF2 mutations such as N196K and G301D are capable of dysregulating gene expression for neoplastic transformation.

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.

Treating Biomacromolecular Conformational Variability

The function of a biomacromolecule is related not only to its structure but also to the different conformations that its structural elements can sample. It is therefore important to determine the extent of the structural fluctuations and to identify the states that are actually populated as a result of the rearrangement. However, this accomplishment is undermined by an intrinsic limitation: the amount of experimental data is by and large inferior to the number of the states that a biomacromolecule can actually sample. This means that additional, a priori information must be applied in order to derive the most from the available experimental data but not to run into overinterpretation. In this chapter we will give a summary of the experimental observables that can be used towards the reconstruction of structural ensembles, how the data can be profitably combined and to what extent the data are affected by error; finally we will give an overview of the computational methods that have been developed to model structural ensembles, highlighting their difference and similarities, advantages and disadvantages.

Capturing dynamic conformational shifts in protein–ligand recognition using integrative structural biology in solution

In recent years, a dynamic view of the structure and function of biological macromolecules is emerging, highlighting an essential role of dynamic conformational equilibria to understand molecular mechanisms of biological functions. The structure of a biomolecule, i.e. protein or nucleic acid in solution, is often best described as a dynamic ensemble of conformations, rather than a single structural state. Strikingly, the molecular interactions and functions of the biological macromolecule can then involve a shift between conformations that pre-exist in such an ensemble. Upon external cues, such population shifts of pre-existing conformations allow gradually relaying the signal to the downstream biological events. An inherent feature of this principle is conformational dynamics, where intrinsically disordered regions often play important roles to modulate the conformational ensemble. Unequivocally, solution-state NMR spectroscopy is a powerful technique to study the structure and dynamics of such biomolecules in solution. NMR is increasingly combined with complementary techniques, including fluorescence spectroscopy and small angle scattering. The combination of these techniques provides complementary information about the conformation and dynamics in solution and thus affords a comprehensive description of biomolecular functions and regulations. Here, we illustrate how an integrated approach combining complementary techniques can assess the structure and dynamics of proteins and protein complexes in solution.

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.

RNA and Proteins: Mutual Respect

Proteins and RNA are often found in ribonucleoprotein particles (RNPs), where they function in cellular processes to synthesize proteins (the ribosome), chemically modify RNAs (small nucleolar RNPs), splice pre-mRNAs (the spliceosome), and, on a larger scale, sequester RNAs, degrade them, or process them (P bodies, Cajal bodies, and nucleoli). Each RNA–protein interaction is a story in itself, as both molecules can change conformation, compete for binding sites, and regulate cellular functions. Recent studies of Xist long non-coding RNP, the U4/5/6 tri-small nuclear RNP complex, and an activated state of a spliceosome reveal new features of RNA interactions with proteins, and, although their stories are incomplete, they are already fascinating.

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.

Structural Analysis of Protein-RNA Complexes in Solution Using NMR Paramagnetic Relaxation Enhancements

Biological activity in the cell is predominantly mediated by large multiprotein and protein-nucleic acid complexes that act together to ensure functional fidelity. Nuclear magnetic resonance (NMR) spectroscopy is the only method that can provide information for high-resolution three-dimensional structures and the conformational dynamics of these complexes in solution. Mapping of binding interfaces and molecular interactions along with the characterization of conformational dynamics is possible for very large protein complexes. In contrast, de novo structure determination by NMR becomes very time consuming and difficult for protein complexes larger than 30 kDa as data are noisy and sparse. Fortunately, high-resolution structures are often available for individual domains or subunits of a protein complex and thus sparse data can be used to define their arrangement and dynamics within the assembled complex. In these cases, NMR can therefore be efficiently combined with complementary solution techniques, such as small-angle X-ray or neutron scattering, to provide a comprehensive description of the structure and dynamics of protein complexes in solution. Particularly useful are NMR-derived paramagnetic relaxation enhancements (PREs), which provide long-range distance restraints (ca. 20Å) for structural analysis of large complexes and also report on conformational dynamics in solution. Here, we describe the use of PREs from sample production to structure calculation, focusing on protein-RNA complexes. On the basis of recent examples from our own research, we demonstrate the utility, present protocols, and discuss potential pitfalls when using PREs for studying the structure and dynamic features of protein-RNA complexes.

Transient electrostatic interactions dominate the conformational equilibrium sampled by multidomain splicing factor U2AF65: a combined NMR and SAXS study

Multidomain proteins containing intrinsically disordered linkers exhibit large-scale dynamic modes that play key roles in a multitude of molecular recognition and signaling processes. Here, we determine the conformational space sampled by the multidomain splicing factor U2AF65 using complementary nuclear magnetic resonance spectroscopy and small-angle scattering data. Available degrees of conformational freedom are initially stochastically sampled and experimental data then used to delineate the potential energy landscape in terms of statistical probability. The spatial distribution of U2AF65 conformations is found to be highly anisotropic, comprising significantly populated interdomain contacts that appear to be electrostatic in origin. This hypothesis is supported by the reduction of signature PREs reporting on expected interfaces with increasing salt concentration. The described spatial distribution reveals the complete spectrum of the unbound forms of U2AF65 that coexist with the small percentage of a preformed RNA-bound domain arrangement required for polypyrimidine-tract recognition by conformational selection. More generally, the proposed approach to describing conformational equilibria of multidomain proteins can be further combined with other experimental data that are sensitive to domain dynamics.

Iron availability modulates aberrant splicing of ferrochelatase through the iron- and 2-oxoglutarate dependent dioxygenase Jmjd6 and U2AF(65.)

Erythropoietic protoporphyria (EPP) results from partial deficiency of ferrochelatase (FECH). Genetically, EPP patients differ from asymptomatic mutation carriers at the unmutated FECH allele, the expression of which is modulated by single nucleotide polymorphism IVS3-48C/T. The IVS3-48C genotype, which is present among patients, leads to correct splicing of 60% of the pre-mRNA and to alternative splicing of 40%, the latter mRNA-product being destroyed by nonsense-mediated decay. An IVS3-48T genotype generates 80% correct and 20% aberrant products. Our study demonstrated that under iron deficient conditions, the aberrant splice product was increased to 56% and 50% of total FECH mRNA in erythroleukemic K562 and lymphoblastoid cell lines, respectively, both being homozygous for IVS3-48T. Concomitantly, FECH protein was decreased. Iron deficiency had less effect on the FECH splice ratio in an IVS3-48C/C lymphoblastoid cell line. Effects similar to iron deficiency were generated by siRNA knockdown of either splicing factor U2AF(65) or Fe(II)- and 2-oxoglutarate-dependent dioxygenase Jumonji domain-containing protein 6 (Jmjd6), which interacts with U2AF(65) by lysyl-hydroxylation. Based on these results, we propose that the availability of iron, a co-factor of Jmjd6, modulates U2AF(65)-lysyl-hydroxylation. This in turn, influences the relative amounts of correct and aberrant FECH mRNA splice products and thus, regulates the FECH enzyme activity.

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

Degenerate splice site sequences mark the intron boundaries of pre-mRNA transcripts in multicellular eukaryotes. The essential pre-mRNA splicing factor U2AF⁶⁵ is faced with the paradoxical tasks of accurately targeting polypyrimidine (Py) tracts preceding 3′ splice sites while adapting to both cytidine and uridine nucleotides with nearly equivalent frequencies. To understand how U2AF⁶⁵ recognizes degenerate Py tracts, we determined six crystal structures of human U2AF⁶⁵ bound to cytidine-containing Py tracts. As deoxy-ribose backbones were required for co-crystallization with these Py tracts, we also determined two baseline structures of U2AF⁶⁵ bound to the deoxy-uridine counterparts and compared the original, RNA-bound structure. Local structural changes suggest that the N-terminal RNA recognition motif 1 (RRM1) is more promiscuous for cytosine-containing Py tracts than the C-terminal RRM2. These structural differences between the RRMs were reinforced by the specificities of wild-type and site-directed mutant U2AF⁶⁵ for region-dependent cytosine- and uracil-containing RNA sites. Small-angle X-ray scattering analyses further demonstrated that Py tract variations select distinct inter-RRM spacings from a pre-existing ensemble of U2AF⁶⁵ conformations. Our results highlight both local and global conformational selection as a means for universal 3′ splice site recognition by U2AF⁶⁵.

RRM-RNA recognition: NMR or crystallography…and new findings

To characterize protein-RNA recognition at the molecular level, structural biology has turned out to be an indispensable approach. Detailed and direct insights into the mechanism of RNA binding and specificity have emerged from protein-RNA structures, especially from the most abundant RNA recognition motif (RRM). Although this protein domain has a very conserved α-β fold, it can recognize a large number of different RNA sequences and shapes and can be involved in a multitude of biological processes. Here, we report on recent single and multiple RRM-RNA structures and point out those features that provide novel insights into the mechanism of RNA recognition by RRMs. We further outline inherent problems to both NMR spectroscopy and X-ray crystallography methods and review recent strategies that emphasize the need to use both methodologies for more rapid and accurate structure determinations.