Добавил:
Sekretar
kiopkiopkiop18@yandex.ru
t.me/Prokururor I Вовсе не секретарь, но почту проверяю
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:Ординатура / Хирургия / Библиотека им академика М.И. Перельмана / Книга_5644_Библиотеки_им_академика_М_И_Перельмана
.pdf
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 modied as well. These low-frequency correlated modication 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 suer three main limitations [60]: (i) the
preference of psoralen for pyrimidines, particularly uracils, can result in the preferential 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 signicantly 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-modied psoralen derivatives to enable the streptavidin-mediated enrichment 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 nonspecic ligation events, the COMRADES approach
generates a control sample by reversing psoralen cross-links and melting RNA
duplexes prior to ligation.
Concerning the identication of RNA–RNA interactions via chemical probing and
MaP [35, 48], the main limitation is that equilibrium uctuations are very rare, so
are correlated modication 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 modeling, achieving signicant 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 signicantly 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 alternative 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 specic 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 experiments. Indeed, while RNA–RNA interaction and spatial proximity mapping
methods capture individual helices or internucleotide distances, chemical probing
experiments typically only provide a reactivity prole 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 modication
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 conformation). 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 proles, enabling the
reconstruction of the reactivity proles 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 dierent 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 identied
clusters corresponding to the number of conformations making up the ensemble.
Although RT drop-o-based probing experiments do not preserve any information about the relationships between individual nucleotides in the RNA strand, as
each cDNA molecule can only capture a single modication site, several computational methods have been devised in an attempt to mitigate this limitation [71–73].
These algorithms take advantage of a partition function-based approach to sample a large number of possible structures from the theoretical ensemble a given
RNA can populate and then use the experimentally determined reactivity prole 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 structures 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 specic tertiary folds in the ribosomal RNAs (rRNA) of bacteria, as well as from natural ligands of bacterial riboswitches, binding of small molecules to RNA often involves
hydrogen bonding with the Watson-Crick faces of exposed nucleobases and aromatic 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 specifically 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. Nonspecic binders can be readily identied as they will result in
SHAPE reactivity changes for both the target and control structures, while specic
binders will only show signicant 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 denitely 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 identication. 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 eciently 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 deconvolution from MaP-based chemical probing experiments likely only provides a
coarse-grained representation of the actual ensemble, as small local structural
dierences might not be eciently 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). Identication 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). Identication 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 dierent 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 condence 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.
Соседние файлы в папке Библиотека им академика М.И. Перельмана
