Добавил:
kiopkiopkiop18@yandex.ru t.me/Prokururor I Вовсе не секретарь, но почту проверяю Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Ординатура / Хирургия / Библиотека им академика М.И. Перельмана / Книга_5644_Библиотеки_им_академика_М_И_Перельмана

.pdf
Скачиваний:
0
Добавлен:
15.09.2026
Размер:
14 Мб
Скачать
☆
2.2 Experimentally Guided RNA Structure Modeling 15
https://t.me/med1917
Figure 2.4 General outline of the protocol for the mapping of RNA–RNA interactions using psoralen (or its derivatives). After having stabilized RNA duplexes by psoralen cross-linking, the RNA is partially fragmented (or digested) to both remove single-stranded loops and to shorten long RNA duplexes. The resulting short RNA duplexes are joined by intramolecular ligation, generating chimeric RNA molecules. Chimeras are reverse transcribed and sequenced. Following mapping of the reads to the reference transcriptome, the two “halves” of the chimeras will map to distal regions of one (or two, in the case of intermolecular duplexes) RNA, hence enabling the direct detection of RNA duplexes.
Psoralen
cross-link
5′
RNA partial digestion, intramolecular ligation,
crosslink reversal and reverse transcription
Sequencing
RNA
duplex
Mapping to reference and
inference of RNA–RNA interactions
16 2 RNA Structure Probing, Dynamics, and Folding
https://t.me/med1917
Probe-induced
chemical modification
Equilibrium fluctuation (structure “breathing”)
Figure 2.5 General outline of the principle exploited by RING-MaP and PAIR-MaP to indirectly infer RNA–RNA interactions from chemical probing experiments. Unpaired bases (i.e. loops) are readily modified by the chemical probe with high efficiency. In rare occasions, base-paired bases can become transiently accessible to the probe due to equilibrium fluctuations. Modification of one partner in the pair permanently destabilizes base pairing, making also the other partner available for modification. Following reverse transcription under mutational profiling (MaP) conditions, sites of modifications are recorded as mutations in the cDNA, and sites of correlated modification due to structure “breathing” can be identified and used to detect putative base-pairing interactions.
Mutational profiling (MaP)
RT
5′ 5′ 5′ 5′ 5′ 5′
Correlated modification events
due to structure “breathing” (rare)
of one of the bases in a pair results in the permanent exposure of its partner, which can then be modied as well. These low-frequency correlated modication events can then be simultaneously detected as mutations within a single cDNA molecule by MaP (Figure 2.5).
2.2.2.1 Limits of RNA–RNA Interaction Mapping
Methods based on psoralen cross-linking suer three main limitations [60]: (i) the preference of psoralen for pyrimidines, particularly uracils, can result in the prefer­ential stabilization and capture of pyrimidine-rich duplexes, hence making these methods not quantitative; (ii) although intramolecular ligation is favored between RNA strands in a duplex, any two RNA strands in a sample can be intermolecularly ligated, resulting in a number of artifactual RNA–RNA interactions being mapped; and (iii) the vast majority of the adapter-ligated RNA fragments are non-chimeric, hence signicantly increasing the required sequencing depth of these experiments. Potential workarounds have been proposed for some of these issues. For instance, SPLASH [54] and COMRADES [57], respectively, take advantage of biotin- or azido-modied psoralen derivatives to enable the streptavidin-mediated enrich­ment of RNA duplexes, while PARIS [53] exploits a native-denatured 2-dimensional gel electrophoresis approach to enrich the cross-linked RNA duplexes. Additionally, in order to account for nonspecic ligation events, the COMRADES approach generates a control sample by reversing psoralen cross-links and melting RNA duplexes prior to ligation.
Concerning the identication of RNA–RNA interactions via chemical probing and MaP [35, 48], the main limitation is that equilibrium uctuations are very rare, so are correlated modication events, hence requiring extreme sequencing depths to be detected, making this type of approach only suited for the targeted analysis of individual RNAs.
2.2 Experimentally Guided RNA Structure Modeling 17
https://t.me/med1917
Spatially proximal and flexible nucleotides
O
Figure 2.6 Reaction of TBIA and SHARC reagents with the ribose 2′-OH of spatially proximal and structurally flexible nucleotides.
Base
O
OHO
Base
O
OHO
O
TBIA, SHARC
O
Base
O
OO
Base
R
O
O
O
OO
O
2.2.3 Mapping Spatially Proximal Nucleotides in RNA molecules
An exciting recent evolution of SHAPE probes has led to the development of bifunctional reagents. These compounds bear two moieties capable of forming adducts with the ribose 2′-OH of two spatially proximal and structurally exible nucleotides. These include trans-bis-isatoic anhydride (TBIA) [61] and spatial 2′-hydroxyl acylation reversible cross-linking (SHARC) reagents [62] (Figure 2.6). The length of the exible linker between the two adduct-forming moieties of SHARC reagents can in principle be adjusted to create “molecular rulers” that can be used to measure approximate internucleotide distances.
Two readout strategies have been proposed for these experiments: selective 2′-hydroxyl acylation analyzed by primer extension and juxtaposed merged pairs (SHAPE-JuMP) and SHARC-exo (Figure 2.7). SHAPE-JuMP [61] exploits an engineered RT, RT-C8 [63], which has the ability to traverse (or jump across) chemical cross-links, hence recording the cross-linked sites as deletions in the cDNA. SHARC-exo [62], instead, involves partial RNA digestion with RNase III, followed by denatured-denatured 2-dimensional gel electrophoresis to enrich cross-linked RNA fragments. Fragments are then digested by RNase R to trim the 3′ends up to ∼5 nucleotides from the cross-linking site, followed by intramolecular ligation to yield chimeric RNA fragments and cross-linking reversal under mild alkaline conditions.
Data from these experiments can be used to constrain RNA 3D structure mod­eling, achieving signicant improvements as compared to unconstrained structure models [61, 62, 64].
2.2.3.1 Limits of Methods for Spatial Proximity Mapping
The biggest limitation of these methods is, analogously to methods for RNA–RNA interaction capture, the low abundance of chimeric RNA fragments in the nal library, which has been reported to be ∼3–15% for SHARC-exo [62]. This is possibly the consequence of the higher frequency at which TBIA/SHARC mono-adducts, which occur when one moiety of the probe reacts with a nucleotide while the other is hydrolyzed in water [61], are formed as compared to di-adducts. This limitation signicantly increases the required sequencing depths, hence posing an important challenge to the feasibility of transcriptome-wide analyses.
18 2 RNA Structure Probing, Dynamics, and Folding
https://t.me/med1917
TBIA/SHARC
cross-link
SHAPE-JuMP
Reverse
transcription
5′
5′
3′
5′
Jump
(cross-linking site)
Sequencing
Mapping to reference and
inference of through-space interactions
Through-space
interaction
RNase R
3′
5′
SHARC-exo
Cross-linking sites
5′
}
~5nt
RNA partial digestion, intramolecular ligation, cross-link reversal and
reverse transcription
Through-space
interaction
}
~5nt
Figure 2.7 General outline of the SHAPE-JuMP and SHARC-exo protocols for the mapping of through-space interactions, following cross-linking of spatially proximal nucleotides by using either trans-bis-isatoic anhydride (TBIA) or spatial 2′-hydroxyl acylation reversible cross-linking (SHARC) reagents. The SHAPE-JuMP approach (left) uses a special RT enzyme capable of “jumping” across the reagent-induced cross-links, hence permanently recording the cross-linking site as a deletion in the cDNA. The SHARC-exo approach, instead, involves a partial digestion of the RNA, followed by enrichment of the cross-linked RNA fragments by denatured-denatured 2-dimensional gel electrophoresis. The recovered fragments are shortened down to ∼5 nucleotides from the site of cross-linking at their 3′ends by RNase R digestion and joined by intramolecular ligation. Cross-linking is then reversed, allowing reverse transcription of the chimeric RNA fragments. Following mapping of the reads to the reference transcriptome, the two “halves” of the chimeras will map to distal regions of the RNA. For SHAPE-JuMP, as the RT “jumps” across the cross-link, the coordinates of the cross-linked nucleotides can be determined from the start and end positions of the resulting deletion. For SHARC-exo, as the cross-linked RNA fragments are trimmed by RNase R digestion, the coordinates of the cross-linked nucleotides can be approximately inferred as the positions ∼5 nucleotides upstream of the 3′end of the two “halves” of the chimeras.
2.3 Dealing with RNA Structure Heterogeneity 19
https://t.me/med1917
2.3 Dealing with RNA Structure Heterogeneity
As previously mentioned, an important feature of RNA molecules is their ability to sample a potentially very large conformational space. Many cellular RNAs are likely to populate multiple alternative structures that can coexist as part of a heterogeneous and dynamic ensemble [3, 60]. The ability to switch between alter­native conformations has been reported to be crucial for several RNAs. Bacterial riboswitches [65] and RNA thermometers [66], which are capable of altering their structure in response to the presence of specic metabolites or due to temperature changes, respectively, are possibly the better characterized examples.
It is therefore paramount to point out that data generated using any of the above-described methods must not be interpreted just as a function of a single RNA structure, but rather as an aggregate of all the structures populating the ensemble for a given RNA. This point is particularly relevant for chemical probing exper­iments. Indeed, while RNA–RNA interaction and spatial proximity mapping methods capture individual helices or internucleotide distances, chemical probing experiments typically only provide a reactivity prole representing a weighted average of the reactivities of all the structures in the ensemble.
The recent introduction of the MaP-based readout provided an important workaround to this problem. As discussed earlier, the key advantage of MaP over the traditional RT drop-o-based readout is that multiple sites of modication within a single RNA molecule can be all recorded as mutations in a single cDNA product. Therefore, this type of approach allows preserving information about the relationship between the individual nucleotides in the RNA strand (i.e. nucleotides that were simultaneously unpaired, or structurally exible, in a certain conforma­tion). By analyzing the co-mutation patterns in the cDNA molecules, it is possible to learn these relationships and to determine how many alternative structures populate the ensemble for a given RNA (Figure 2.8). For instance, given a set of three nucleotides (e.g. 1, 2, and 3) on an RNA, one would expect to observe any of the three possible pairs of co-mutations (i.e. 1-2, 2-3, and 1-3) to occur with comparable probability for an RNA folding into a single structure.
Deviations from this expected behavior would suggest the presence of multiple structures. For example, if only 1-2 and 2-3 were observed to co-mutate, this would indicate that nucleotides 1 and 3 tend not to be simultaneously reactive to the chemical probe, indicating the presence of two mutually exclusive structures. After having determined the number of conformations making up the ensemble, the cDNAs can be clustered based on their co-mutation proles, enabling the reconstruction of the reactivity proles for the individual conformations. To date, three main computational methods exploiting this approach to reconstruct RNA structure ensembles from MaP chemical probing experiments have been developed, namely detection of RNA folding ensembles using expectation–maximization (DREEM) [67], deconvolution of RNA alternative conformations (DRACO) [68], and deconvolution and annotation of ribonucleic conformational ensembles (DANCE-MaP) [69]. A fourth method, dubbed determination of the variation of the RNA structure conformation through stochastic context-free grammar
20 2 RNA Structure Probing, Dynamics, and Folding
https://t.me/med1917
Co-mutation patterns observed in cDNAs
18 36 50
9
Co-mutation patterns absent in cDNAs
18
9
50
36
Conf. #A
Conf. #B
Both #A & #B
Neither #A or #B
18
50
36
9
60
1
62
Conformation #A
Unpaired only in one conformation
1
62
Conformation #B
Figure 2.8 Outline of the principle exploiting chemical probing and mutational profiling (MaP) analysis for the deconvolution of coexisting alternative RNA structures. In the depicted example, the same RNA populates an ensemble of two conformations (A and B). Both conformations share a number of unpaired bases, for example, bases 9 and 36 (in green). Some bases, however , are unpaired only in one conformation and base paired in the other (in red). Among these, for example, base 18 is unpaired in conformation A and base paired in conformation B, while base 50 is base paired in conformation A and unpaired in conformation B. Chemical probing and MaP analysis of this ensemble will result in a very specific set of cDNA co-mutation patterns. For instance, co-mutations of bases 9 and 36 can be observed both for conformations A and B. Co-mutations of either (or both) bases with base 18 can only be observed for conformation A, while co-mutations with base 50 can only be observed for conformation B. However, as base pairing of bases 18 and 50 is mutually exclusive between the two conformations, reads carrying co-mutations of bases 18 and 50 will (almost) never be observed.
(DaVinci) [70], has been recently introduced. DaVinci uses a dierent approach, combining thermodynamics-based structure prediction and MaP data (Figure 2.9). As each sequencing read provides information on the pairing status of the nucleotides in a single RNA molecule, DaVinci predicts a structure from each
Unpaired in both conformations
Mutational profiling (MaP)
https://t.me/med1917
RT
Dimension 2
2.3 Dealing with RNA Structure Heterogeneity 21
MaP-constrained thermodynamics-driven structure prediction
5′ 5′ 5′ 5′ 5′ 5′
Clustering
Predicted ensemble
Dimension 1
Figure 2.9 Outline of the DaVinci stochastic context-free grammar approach. Following chemical probing, mutational profiling (MaP) analysis, and sequencing, the mutational profile of each read, corresponding to the bases being simultaneously unpaired in a single original RNA molecule, is used to constrain the thermodynamics-driven structure prediction. This results in N predicted secondary structures, each derived from a single read (so, theoretically, from a single original RNA molecule). This predicted ensemble is then subjected to clustering. The number of clusters identified by the analysis corresponds to the number of alternative conformations making up the ensemble for the RNA in analysis.
read by constraining the mutated bases in the read to be unpaired in the predicted structure. The resulting structures are then clustered, with the number of identied clusters corresponding to the number of conformations making up the ensemble.
Although RT drop-o-based probing experiments do not preserve any informa­tion about the relationships between individual nucleotides in the RNA strand, as each cDNA molecule can only capture a single modication site, several computa­tional methods have been devised in an attempt to mitigate this limitation [71–73]. These algorithms take advantage of a partition function-based approach to sam­ple a large number of possible structures from the theoretical ensemble a given RNA can populate and then use the experimentally determined reactivity prole to select a parsimonious set of structures that better explains the observed experimental
22 2 RNA Structure Probing, Dynamics, and Folding
https://t.me/med1917
reactivities. All these methods have been thoroughly reviewed elsewhere [74], so they will not be discussed in detail in this chapter.
2.4 Querying RNA–Small Molecule Interactions with Chemical Probing
Determining the structure of RNA molecules in the context of living cells has recently attracted huge interest in industry [75–77]. Just like proteins, RNA struc­tures possess ligandable pockets,amenable for targeting small molecules [78]. In this perspective, mapping which structures are functional and druggable is paramount for the development of innovative RNA-targeted therapeutic strategies, possibly opening the way to targeting proteins currently considered to be undruggable.
As learned from antibiotics such as aminoglycosides, which target specic ter­tiary folds in the ribosomal RNAs (rRNA) of bacteria, as well as from natural lig­ands of bacterial riboswitches, binding of small molecules to RNA often involves hydrogen bonding with the Watson-Crick faces of exposed nucleobases and aro­matic stacking interactions [79]. These interactions can be captured and mapped by chemical probing, providing a powerful tool to investigate RNA-small molecule target engagement, as well as o-target interactions, and the impact of RNA-small molecule interactions on the structure of RNA molecules.
For example, the antibiotic spectinomycin has been previously shown to specif­ically interact with residues G1064 and C1192 in the E. coli 16S rRNA [80]. Accordingly, DMS chemical probing readout by MaP of the 16S rRNA from living bacteria treated with spectinomycin shows strong protection of C1192 from DMS due to hydrogen bonding between C1192 and spectinomycin [81]. More recently, SHAPE probing has been proposed as a rapid way to conduct high-throughput fragment-based screening of small-molecule RNA ligands [82]. In this method, a synthetic RNA molecule is designed to contain both the target RNA structure and a control structure. Nonspecic binders can be readily identied as they will result in SHAPE reactivity changes for both the target and control structures, while specic binders will only show signicant reactivity changes on the target structure. An alternative SHAPE-based approach has also been proposed to map the interaction of small molecules with RNA on a transcriptome-wide scale [83]. In this approach, the small molecule of interest is functionalized with an acylimidazole-substituted linker to enable adduct formation with the ribose 2′-OH at structurally exible nucleotides, further allowing the detection of small molecule–RNA interactions via RT drop-o or MaP-based readouts.
2.5 Conclusions and Future Prospects
Despite substantial advances, we are still far from being able to accurately map the structure of cellular transcriptomes. Of all the potential limitations and challenges of the methods discussed in this chapter, dealing with RNA structure heterogeneity
References 23
https://t.me/med1917
in the cellular context is denitely the biggest. Although several proof-of-concept studies demonstrated the targeted deconvolution of the RNA structure ensembles of individual transcripts, to date no one has attempted exploring RNA structure heterogeneity on a transcriptome-wide scale. Achieving such resolution of the RNA structurome is essential in the perspective of using RNA as a drug target. Pathologically relevant RNAs might populate multiple conformations, and only one of these conformations might actually be responsible for the pathological phenotype, hence making the deconvolution of the RNA structure ensemble a crucial step in target identication. Meeting this challenge requires the cooperative development of improved experimental and computational strategies. For instance, ensemble deconvolution from MaP-based chemical probing experiments has so far been attempted mostly in conjunction with DMS probing, owing to the typically higher signal-to-noise ratio of DMS as compared to other probes, such as SHAPE reagents. This, however, poses substantial limits to our ability to eciently identify structurally heterogeneous regions as DMS can only probe roughly half of the bases in RNA molecules. The increased signal-to-noise ratio of recently developed SHAPE probes [42], possibly combined with the development of engineered RTs with improved read-through apabilities [84], might help advance with respect to the current state of the art. Furthermore, although powerful, ensemble decon­volution from MaP-based chemical probing experiments likely only provides a coarse-grained representation of the actual ensemble, as small local structural dierences might not be eciently captured. In this regard, the ability to combine the information from orthogonal probing reagents (e.g. DMS and SHAPE), or from orthogonal experimental approaches (e.g. chemical probing and RNA duplex mapping), might help improve the resolution over each individual method, possibly providing more holistic solutions to map cellular RNA structure ensembles.
References
1 ENCODE Project Consortium, Moore, J.E., Purcaro, M.J. et al. (2020). Expanded
encyclopaedias of DNA elements in the human and mouse genomes. Nature 583 (7818): 699–710.
2 Frankish, A., Diekhans, M., Jungreis, I. et al. (2021). GENCODE 2021. Nucleic
Acids Res. 49 (D1): D916–D923.
3 Ganser, L.R., Kelly, M.L., Herschlag, D., and Al-Hashimi, H.M. (2019). The roles
of structural dynamics in the cellular functions of RNAs. Nat. Rev. Mol. Cell Biol. 20 (8): 474–489.
4 Schroeder, R., Barta, A., and Semrad, K. (2004). Strategies for RNA folding and
assembly. Nat. Rev. Mol. Cell Biol. 5 (11): 908–919.
5 Mitterer, V. and Pertschy, B. RNA folding and functions of RNA helicases in
ribosome biogenesis. RNA Biol. 19 (1): 781–810.
6 Leppek, K., Das, R., and Barna, M. (2018). Functional 5’ UTR mRNA structures
in eukaryotic translation regulation and how to nd them. Nat. Rev. Mol. Cell Biol. 19 (3): 158–174.
24 2 RNA Structure Probing, Dynamics, and Folding
https://t.me/med1917
7 Halvorsen, M., Martin, J.S., Broadaway, S., and Laederach, A. (2010).
Disease-associated mutations that alter the RNA structural ensemble. PLos Genet. 6 (8): e1001074.
8 Napierała, M. and Krzyzosiak, W.J. (1997). CUG repeats present in myotonin
kinase RNA form metastable “slippery” hairpins. J Biol Chem 272 (49): 31079–31085.
9 Leppert, J., Urbinati, C.R., Häfner, S. et al. (2004). Identication of NH...N
hydrogen bonds by magic angle spinning solid state NMR in a double-stranded RNA associated with myotonic dystrophy. Nucleic Acids Res. 32 (3): 1177–1183.
10 Timchenko, L.T., Miller, J.W., Timchenko, N.A. et al. (1996). Identication of
a (CUG)n triplet repeat RNA-binding protein and its expression in myotonic dystrophy. Nucleic Acids Res. 24 (22): 4407–4414.
11 Roberts, R., Timchenko, N.A., Miller, J.W. et al. (1997). Altered phosphorylation
and intracellular distribution of a (CUG)n triplet repeat RNA-binding protein in patients with myotonic dystrophy and in myotonin protein kinase knockout mice. Proc. Natl. Acad. Sci. U.S.A. 94 (24): 13221–13226.
12 Philips, A.V., Timchenko, L.T., and Cooper, T.A. (1998). Disruption of splicing
regulated by a CUG-binding protein in myotonic dystrophy. Science 280 (5364): 737–741.
13 Boerneke, M.A., Ehrhardt, J.E., and Weeks, K.M. (2019). Physical and functional
analysis of viral RNA genomes by SHAPE. Annu. Rev. Virol. 6 (1): 93–117.
14 Manfredonia, I. and Incarnato, D. (2020). Structure and regulation of coronavirus
genomes: state-of-the-art and novel insights from SARS-CoV-2 studies. Biochem. Soc. Trans. 49 (1): 341–352.
15 Lee, E., Bujalowski, P.J., Teramoto, T. et al. (2021). Structures of avivirus RNA
promoters suggest two binding modes with NS5 polymerase. Nat. Commun. 12 (1): 2530.
16 Sherpa, C., Rausch, J.W., Le Grice, S.F.J. et al. (2015). The HIV-1 Rev response
element (RRE) adopts alternative conformations that promote dierent rates of virus replication. Nucleic Acids Res. 43 (9): 4676–4686.
17 Zuker, M. and Sanko, D. (1984). RNA secondary structures and their prediction.
Bull. Math. Biol. 46 (4): 591–621.
18 Chen, S.-J. (2008). RNA folding: conformational statistics, folding kinetics, and
ion electrostatics. Annu. Rev. Biophys. 37: 197–214.
19 Herschlag, D. (1995). RNA chaperones and the RNA folding problem. J. Biol.
Chem. 270 (36): 20871–20874.
20 Xia, T., SantaLucia, J., Burkard, M.E. et al. (1998). Thermodynamic parameters
for an expanded nearest-neighbor model for formation of RNA duplexes with Watson-Crick base pairs. Biochemistry 37 (42): 14719–14735.
21 Mathews, D.H. (2004). Using an RNA secondary structure partition function to
determine condence in base pairs predicted by free energy minimization. RNA 10 (8): 1178–1190.
22 Zhang, J., Fei, Y., Sun, L., and Zhang, Q.C. (2022). Advances and opportunities
in RNA structure experimental determination and computational modeling. Nat. Methods 19 (10): 1193–1207.