Extending the Size of Protein-RNA Complexes Studied by Nuclear Magnetic Resonance Spectroscopy

PRPI-SC: an ensemble deep learning model for predicting plant lncRNA-protein interactions

BACKGROUND

Plant long non-coding RNAs (lncRNAs) play vital roles in many biological processes mainly through interactions with RNA-binding protein (RBP). To understand the function of lncRNAs, a fundamental method is to identify which types of proteins interact with the lncRNAs. However, the models or rules of interactions are a major challenge when calculating and estimating the types of RBP.

RESULTS

In this study, we propose an ensemble deep learning model to predict plant lncRNA-protein interactions using stacked denoising autoencoder and convolutional neural network based on sequence and structural information, named PRPI-SC. PRPI-SC predicts interactions between lncRNAs and proteins based on the k-mer features of RNAs and proteins. Experiments proved good results on Arabidopsis thaliana and Zea mays datasets (ATH948 and ZEA22133). The accuracy rates of ATH948 and ZEA22133 datasets were 88.9% and 82.6%, respectively. PRPI-SC also performed well on some public RNA protein interaction datasets.

CONCLUSIONS

PRPI-SC accurately predicts the interaction between plant lncRNA and protein, which plays a guiding role in studying the function and expression of plant lncRNA. At the same time, PRPI-SC has a strong generalization ability and good prediction effect for non-plant data.

In vitro iCLIP-based modeling uncovers how the splicing factor U2AF2 relies on regulation by cofactors

Alternative splicing generates distinct mRNA isoforms and is crucial for proteome diversity in eukaryotes. The RNA-binding protein (RBP) U2AF2 is central to splicing decisions, as it recognizes 3′ splice sites and recruits the spliceosome. We establish “in vitro iCLIP” experiments, in which recombinant RBPs are incubated with long transcripts, to study how U2AF2 recognizes RNA sequences and how this is modulated by trans-acting RBPs. We measure U2AF2 affinities at hundreds of binding sites and compare in vitro and in vivo binding landscapes by mathematical modeling. We find that trans-acting RBPs extensively regulate U2AF2 binding in vivo, including enhanced recruitment to 3′ splice sites and clearance of introns. Using machine learning, we identify and experimentally validate novel trans-acting RBPs (including FUBP1, CELF6, and PCBP1) that modulate U2AF2 binding and affect splicing outcomes. Our study offers a blueprint for the high-throughput characterization of in vitro mRNP assembly and in vivo splicing regulation.

Recognition of the 3' splice site RNA by the U2AF heterodimer involves a dynamic population shift

An essential early step in the assembly of human spliceosomes onto pre-mRNA involves the recognition of regulatory RNA cis elements in the 3' splice site by the U2 auxiliary factor (U2AF). The large (U2AF65) and small (U2AF35) subunits of the U2AF heterodimer contact the polypyrimidine tract (Py-tract) and the AG-dinucleotide, respectively. The tandem RNA recognition motif domains (RRM1,2) of U2AF65 adopt closed/inactive and open/active conformations in the free form and when bound to bona fide Py-tract RNA ligands. To investigate the molecular mechanism and dynamics of 3' splice site recognition by U2AF65 and the role of U2AF35 in the U2AF heterodimer, we have combined single-pair FRET and NMR experiments. In the absence of RNA, the RRM1,2 domain arrangement is highly dynamic on a submillisecond time scale, switching between closed and open conformations. The addition of Py-tract RNA ligands with increasing binding affinity (strength) gradually shifts the equilibrium toward an open conformation. Notably, the protein-RNA complex is rigid in the presence of a strong Py-tract but exhibits internal motion with weak Py-tracts. Surprisingly, the presence of U2AF35, whose UHM domain interacts with U2AF65 RRM1, increases the population of the open arrangement of U2AF65 RRM1,2 in the absence and presence of a weak Py-tract. These data indicate that the U2AF heterodimer promotes spliceosome assembly by a dynamic population shift toward the open conformation of U2AF65 to facilitate the recognition of weak Py-tracts at the 3' splice site. The structure and RNA binding of the heterodimer was unaffected by cancer-linked myelodysplastic syndrome mutants.

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.

NMR approaches for structural analysis of multidomain proteins and complexes in solution

NMR spectroscopy is a key method for studying the structure and dynamics of (large) multidomain proteins and complexes in solution. It plays a unique role in integrated structural biology approaches as especially information about conformational dynamics can be readily obtained at residue resolution. Here, we review NMR techniques for such studies focusing on state-of-the-art tools and practical aspects. An efficient approach for determining the quaternary structure of multidomain complexes starts from the structures of individual domains or subunits. The arrangement of the domains/subunits within the complex is then defined based on NMR measurements that provide information about the domain interfaces combined with (long-range) distance and orientational restraints. Aspects discussed include sample preparation, specific isotope labeling and spin labeling; determination of binding interfaces and domain/subunit arrangements from chemical shift perturbations (CSP), nuclear Overhauser effects (NOEs), isotope editing/filtering, cross-saturation, and differential line broadening; and based on paramagnetic relaxation enhancements (PRE) using covalent and soluble spin labels. Finally, the utility of complementary methods such as small-angle X-ray or neutron scattering (SAXS, SANS), electron paramagnetic resonance (EPR) or fluorescence spectroscopy techniques is discussed. The applications of NMR techniques are illustrated with studies of challenging (high molecular weight) protein complexes.

A fully enzymatic method for site-directed spin labeling of long RNA

Site-directed spin labeling is emerging as an essential tool to investigate the structural and dynamical features of RNA. We propose here an enzymatic method, which allows the insertion of a paramagnetic center at a specific position in an RNA molecule. The technique is based on a segmental approach using a ligation protocol with T4 RNA ligase 2. One transcribed acceptor RNA is ligated to a donor RNA in which a thio-modified nucleotide is introduced at its 5′-end by in vitro transcription with T7 RNA polymerase. The paramagnetic thiol-specific reagent is subsequently attached to the RNA ligation product. This novel strategy is demonstrated by introducing a paramagnetic probe into the 55 nucleotides long RNA corresponding to K-turn and Specifier Loop domains from the Bacillus subtilis tyrS T-Box leader RNA. The efficiency of the coupling reaction and the quality of the resulting spin-labeled RNA were assessed by Mass Spectrometry, Electron Paramagnetic Resonance (EPR) and Nuclear Magnetic Resonance (NMR). This method enables various combinations of isotopic segmental labeling and spin labeling schemes, a strategy that will be of particular interest to investigate the structural and dynamical properties of large RNA complexes by NMR and EPR spectroscopies.

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.

The dynamic duo: combining NMR and small angle scattering in structural biology

Structural biology provides essential information for elucidating molecular mechanisms that underlie biological function. Advances in hardware, sample preparation, experimental methods, and computational approaches now enable structural analysis of protein complexes with increasing complexity that more closely represent biologically entities in the cellular environment. Integrated multidisciplinary approaches are required to overcome limitations of individual methods and take advantage of complementary aspects provided by different structural biology techniques. Although X-ray crystallography remains the method of choice for structural analysis of large complexes, crystallization of flexible systems is often difficult and does typically not provide insights into conformational dynamics present in solution. Nuclear magnetic resonance spectroscopy (NMR) is well-suited to study dynamics at picosecond to second time scales, and to map binding interfaces even of large systems at residue resolution but suffers from poor sensitivity with increasing molecular weight. Small angle scattering (SAS) methods provide low resolution information in solution and can characterize dynamics and conformational equilibria complementary to crystallography and NMR. The combination of NMR, crystallography, and SAS is, thus, very useful for analysis of the structure and conformational dynamics of (large) protein complexes in solution. In high molecular weight systems, where NMR data are often sparse, SAS provides additional structural information and can differentiate between NMR-derived models. Scattering data can also validate the solution conformation of a crystal structure and indicate the presence of conformational equilibria. Here, we review current state-of-the-art approaches for combining NMR, crystallography, and SAS data to characterize protein complexes in solution.

NMR and small-angle scattering-based structural analysis of protein complexes in solution

Structural analysis of multi-domain protein complexes is a key challenge in current biology and a prerequisite for understanding the molecular basis of essential cellular processes. The use of solution techniques is important for characterizing the quaternary arrangements and dynamics of domains and subunits of these complexes. In this respect solution NMR is the only technique that allows atomic- or residue-resolution structure determination and investigation of dynamic properties of multi-domain proteins and their complexes. As experimental NMR data for large protein complexes are sparse, it is advantageous to combine these data with additional information from other solution techniques. Here, the utility and computational approaches of combining solution state NMR with small-angle X-ray and Neutron scattering (SAXS/SANS) experiments for structural analysis of large protein complexes is reviewed. Recent progress in experimental and computational approaches of combining NMR and SAS are discussed and illustrated with recent examples from the literature. The complementary aspects of combining NMR and SAS data for studying multi-domain proteins, i.e. where weakly interacting domains are connected by flexible linkers, are illustrated with the structural analysis of the tandem RNA recognition motif (RRM) domains (RRM1-RRM2) of the human splicing factor U2AF65 bound to a nine-uridine (U9) RNA oligonucleotide.

The impact of window functions on NMR-based paramagnetic relaxation enhancement measurements in membrane proteins

Though challenging, solution NMR spectroscopy allows fundamental interrogation of the structure and dynamics of membrane proteins. One major technical hurdle in studies of helical membrane proteins by NMR is the difficulty of obtaining sufficient long range NOEs to determine tertiary structure. For this reason, long range distance information is sometimes sought through measurement of paramagnetic relaxation enhancements (PRE) of NMR nuclei as a function of distance from an introduced paramagnetic probe. Current PRE interpretation is based on the assumption of Lorentzian resonance lineshapes. However, in order to optimize spectral resolution, modern multidimensional NMR spectra are almost always subjected to resolution-enhancement, leading to distortions in the Lorentizian peak shape. Here it is shown that when PREs are derived using peak intensities (i.e., peak height) and linewidths from both real and simulated spectra that were produced using a wide range of apodization/window functions, that there is little variation in the distances determined (<1 A at the extremes). This indicates that the high degree of resolution enhancement required to obtain well-resolved spectra from helical membrane proteins is compatible with the use of PRE data as a source of distance restraints. While these conclusions are particularly important for helical membrane proteins, they are generally applicable to all PRE measurements made using resolution-enhanced data.

Preparation of selective and segmentally labeled single-stranded DNA for NMR by self-primed PCR and asymmetrical endonuclease double digestion

We demonstrate a new, efficient and easy-to-use method for enzymatic synthesis of (stereo-)specific and segmental (13)C/(15)N/(2)H isotope-labeled single-stranded DNA in amounts sufficient for NMR, based on the highly efficient self-primed PCR. To achieve this, new approaches are introduced and combined. (i) Asymmetric endonuclease double digestion of tandem-repeated PCR product. (ii) T4 DNA ligase mediated ligation of two ssDNA segments. (iii) In vitro dNTP synthesis, consisting of in vitro rNTP synthesis followed by enzymatic stereo-selective reduction of the C2' of the rNTP, and a one-pot add-up synthesis of dTTP from dUTP. The method is demonstrated on two ssDNAs: (i) a 36-nt three-way junction, selectively (13)C(9)/(15)N(3)/(2)H((1',2'',3',4',5',5''))-dC labeled and (ii) a 39-nt triple-repeat three-way junction, selectively (13)C(9)/(15)N(3)/(2)H((1',2'',3',4',5',5''))-dC and (13)C(9)/(15)N(2)/(2)H((1',2'',3',4',5',5''))-dT labeled in segment C20-C39. Their NMR spectra show the spectral simplification, while the stereo-selective (2)H-labeling in the deoxyribose of the dC-residues, straightforwardly provided assignment of their C1'-H2' and C2'-H2' resonances. The labeling protocols can be extended to larger ssDNA molecules and to more than two segments.

The spliceosome: design principles of a dynamic RNP machine

Ribonucleoproteins (RNPs) mediate key cellular functions such as gene expression and its regulation. Whereas most RNP enzymes are stable in composition and harbor preformed active sites, the spliceosome, which removes noncoding introns from precursor messenger RNAs (pre-mRNAs), follows fundamentally different strategies. In order to provide both accuracy to the recognition of reactive splice sites in the pre-mRNA and flexibility to the choice of splice sites during alternative splicing, the spliceosome exhibits exceptional compositional and structural dynamics that are exploited during substrate-dependent complex assembly, catalytic activation, and active site remodeling.

Use of residual dipolar couplings in structural analysis of protein-ligand complexes by solution NMR spectroscopy

Investigation of structure-function relationships in protein complexes, specifically protein-ligand interactions, carry great significance in elucidating the structural and mechanistic bases of molecular recognition events and their role in regulating cell processes. Nuclear magnetic resonance (NMR) spectroscopy is one of the leading structural and analytical techniques in in-depth studies of protein-ligand interactions. Recent advances in NMR methodology such as transverse relaxation-optimized spectroscopy (TROSY) and residual dipolar couplings (RDCs) measured in liquid crystalline alignment medium, offer a viable alternative to traditional nuclear Overhauser enhancement (NOE)-based approaches for structure determination of large protein complexes. RDCs provide a way to constrain the relative orientation of two molecules in complex with each other by aligning their independently determined order tensors. The potential for utilization of RDCs can be extended to proteins with multiple ligands or even multimeric protein-ligand complexes, where symmetry properties of the protein can be taken advantage of. Availability of effective RDC data collection and analysis protocols can certainly aid this process by their incorporation into structure calculation protocols using intramolecular and intermolecular orientational restraints. This chapter discusses in detail some of these protocols including methods for sample preparation in liquid crystalline media, NMR experiments for RDC data collection, as well as software tools for RDC data analysis and protein-ligand complex structure determination.

Solution NMR structure of putidaredoxin-cytochrome P450cam complex via a combined residual dipolar coupling-spin labeling approach suggests a role for Trp106 of putidaredoxin in complex formation

The 58-kDa complex formed between the [2Fe-2S] ferredoxin, putidaredoxin (Pdx), and cytochrome P450cam (CYP101) from the bacterium Pseudomonas putida has been investigated by high-resolution solution NMR spectroscopy. Pdx serves as both the physiological reductant and effector for CYP101 in the enzymatic reaction involving conversion of substrate camphor to 5-exo-hydroxycamphor. In order to obtain an experimental structure for the oxidized Pdx-CYP101 complex, a combined approach using orientational data on the two proteins derived from residual dipolar couplings and distance restraints from site-specific spin labeling of Pdx has been applied. Spectral changes for residues in and near the paramagnetic metal cluster region of Pdx in complex with CYP101 have also been mapped for the first time using (15)N and (13)C NMR spectroscopy, leading to direct identification of the residues strongly affected by CYP101 binding. The new NMR structure of the Pdx-CYP101 complex agrees well with results from previous mutagenesis and biophysical studies involving residues at the binding interface such as formation of a salt bridge between Asp38 of Pdx and Arg112 of CYP101, while at the same time identifying key features different from those of earlier modeling studies. Analysis of the binding interface of the complex reveals that the side chain of Trp106, the C-terminal residue of Pdx and critical for binding to CYP101, is located across from the heme-binding loop of CYP101 and forms non-polar contacts with several residues in the vicinity of the heme group on CYP101, pointing to a potentially important role in complex formation.

A structure refinement protocol combining NMR residual dipolar couplings and small angle scattering restraints

We present the implementation of a target function based on Small Angle Scattering data (Gabel et al. Eur Biophys J 35(4):313-327, 2006) into the Crystallography and NMR Systems (CNS) and demonstrate its utility in NMR structure calculations by simultaneous application of small angle scattering (SAS) and residual dipolar coupling (RDC) restraints. The efficiency and stability of the approach are demonstrated by reconstructing the structure of a two domain region of the 31 kDa nuclear export factor TAP (TIP-associated protein). Starting with the high resolution X-ray structures of the two individual TAP domains, the translational and orientational domain arrangement is refined simultaneously. We tested the stability of the protocol against variations of the SAS target parameters and the number of RDCs and their uncertainties. The activation of SAS restraints results in an improved translational clustering of the domain positions and lifts part of the fourfold degeneracy of their orientations (associated with a single alignment tensor). The resulting ensemble of structures reflects the conformational space that is consistent with the experimental SAS and RDC data. The SAS target function is computationally very efficient. SAS restraints can be activated at different levels of precision and only a limited SAS angular range is required. When combined with additional data from chemical shift perturbation, paramagnetic relaxation enhancement or mutational analysis the SAS refinement is an efficient approach for defining the topology of multi-domain and/or multimeric biomolecular complexes in solution based on available high resolution structures (NMR or X-ray) of the individual domains.

A target function for quaternary structural refinement from small angle scattering and NMR orientational restraints

We present a novel target function based on atomic coordinates that permits quaternary structural refinement of multi-domain protein-protein or protein-RNA complexes. It requires that the high-resolution structures of the individual domains are known and that small angle scattering (SAS) data as well as NMR orientational restraints from residual dipolar couplings (RDCs) of the complex are available. We show that, when used in combination, the translational and rotational restraints contained in SAS intensities and RDCs, respectively, define a target potential function that permits to determine the overall topology of complexes made up of domains with low internal symmetry. We apply the target function on a modestly anisotropic model system, the Barnase/Barstar complex, and discuss factors that influence the structural refinement such as data errors and the geometrical properties of the individual domains.