Abstract
xrRNAs from flaviviruses survive in host cells because of their exceptional dichotomic response to the unfolding action of different enzymes. They can be unwound, and hence copied, by replicases, and yet can resist degradation by exonucleases. How the same stretch of xrRNA can encode such diverse responses is an open question. Here, by using atomistic models and translocation simulations, we uncover an elaborate and directional mechanism for how stress propagates when the two xrRNA ends, \({5}^{\prime}\) and \({3}^{\prime}\), are driven through a pore. Pulling the \({3}^{\prime}\) end, as done by replicases, elicits a progressive unfolding; pulling the \({5}^{\prime}\) end, as done by exonucleases, triggers a counterintuitive molecular tightening. Thus, in what appears to be a remarkable instance of intra-molecular tensegrity, the very pulling of the \({5}^{\prime}\) end is what boosts resistance to translocation and consequently to degradation. The uncovered mechanistic principle might be co-opted to design molecular meta-materials.
Similar content being viewed by others
Introduction
Subgenomic flavivirus RNAs (sfRNAs) are noncoding RNAs that result from the partial degradation of the viral genomic RNAs by cellular Xrn1 exoribonucleases. They accumulate in cells infected by flaviviruses, such as Zika, disrupting several molecular processes of the host. While sfRNAs cannot be degraded by exonucleases, they can still be processed by polymerases and reverse transcriptases1,2,3,4,5.
As part of the ongoing endeavor to explain sfRNAs dichotomic response to processive enzymes, which is at the basis of flaviviruses pathogenicity, a series of targeted in vitro assays, including mutagenesis experiments have been carried out. They clarified that resistance to exonucleases is provided by a ca 71-nucleotide-long subsequence, xrRNA in short, that can withstand degradation at its \({5}^{\prime}\) end by various other exonucleases besides Xrn15. Inspection of xrRNA structureâthe one of Zika is shown in Fig. 1âreveals multiple pseudoknots and a ring-like motif that appears well-suited to dock onto the surface of the exonuclease2,3.
The elaborate architecture of the relatively-short xrRNAs poses general questions about the compliance of xrRNAs to being translocated through the lumen of enzymes working from the \({5}^{\prime}\) to \({3}^{\prime}\) ends, as done by exonucleases, or from the \({3}^{\prime}\) to \({5}^{\prime}\) ones, as done by replicases and reverse transcriptases. What are the mechanistic and thermodynamic underpinnings of xrRNA resistance? Would the same directional resistance be observed in single-molecule setups forcing translocation through narrow pores? And, can one learn general tensegrity principles transferable to molecular design?
Prompted by these questions, we set out to study the microscopic mechanical underpinnings of sfRNAs resistance by using an atomistic model of the Zika xrRNA and stochastic simulations of its driven translocation through a narrow cylindrical pore. The pore-translocation setup provides a convenient abstraction of the resistance of sfRNAs to being driven through various types of exonucleases. In a similar spirit, our model system is exclusively informed with the native three-dimensional organization of the Zika xrRNA, thus discounting sequence-dependent effects and specific interactions with nucleases. More in general, it allows for a comparison with other biologically-motivated contexts where the intra-molecular organization has been investigated in connection with the compliance to translocate through enzymes or synthetic pores6,7,8,9,10,11,12,13,14,15,16.
In our systematic study, we perform hundreds of force-ramped translocation simulations from the \({3}^{\prime}\) and \({5}^{\prime}\) ends. Structural effectors of translocation compliance are established by analysing the network of strained interactions as well as by selectively removing groups of them. BellâEvans analysis and especially metadynamics simulations are then used to obtain a characterization of the free-energy landscape, of the barriers associated with translocation resistance and of the pathway to the effective transition state.
Throughout the explored combination of ramping protocols and temperatures, we observe that translocation from the \({5}^{\prime}\) end requires forces that significantly exceed those needed for the \({3}^{\prime}\) end. This shows that Zika xrRNA architecture, which is the sole input of our model, suffices to encode a significant directional response to translocation. The enhanced hindrance at the \({5}^{\prime}\) end is shown to originate from an elaborate stress redistribution mechanism that has no analog at the \({3}^{\prime}\) end. The much higher free-energy barrier encountered at the \({5}^{\prime}\) end reflects an orders of magnitude difference of the time required to activate translocation from the two ends at constant force.
The results illuminate the structure-based contribution of xrRNA directional resistance and provide a single quantitative framework that recapitulates the resistance to \({5}^{\prime}\to {3}^{\prime}\) exonuclease degradation and \({3}^{\prime}\to {5}^{\prime}\) processability by reverse transcriptases and polymerases. It also indicates that pore-translocation setups, which could verify our simulation results, could be used as valuable probes of xrRNA resistance.
Results
Structure of Zika xrRNA
Figure 1a shows the structure of the Zika xrRNA, a 71-nucleotide long noncoding subgenomic RNA that can resists degradation by Xrn1 and other \({5}^{\prime}\to {3}^{\prime}\) ribonucleases3. Its canonical and noncanonical base pairings are summarized in Fig. 1b. Heavy atoms contacts, including those of sugar and phosphate moieties, are shown in the native contact map of Fig. 1c. The secondary and tertiary organization is articulated over four helices (P1âP4), and multiple pseudoknot-like elements (Pk1âPk5).
Response to pore translocation
We integrated atomistic modeling, molecular dynamics (MD) simulations, force-spectroscopy analysis, and thermodynamic sampling techniques to understand how exactly the molecular architecture of the Zika xrRNA underpins its resistance to being driven through the lumen of exonucleases from the \({5}^{\prime}\) end, while still allowing to be processed and unwound by enzymes, such as polymerases and reverse transcriptases, which work from the \({3}^{\prime}\) end.
Inspired by force-spectroscopy approaches17,18, we adopted a force-ramping protocol where the Zika xrRNA is driven through a cylindrical pore. The pore is 11.7âà wide and runs perpendicularly through a parallelepiped slab that is 19.5âà thick, Fig. 2a. The pore geometry approximates the Xrn1 nuclease lumen19 and allows the passage of only a single RNA strand at a time. Mimicking an electrophoretic setup11,20, the driving force was applied exclusively to the xrRNA P atoms inside the pore, whose longitudinal axis we take as the vertical or z axis of a Cartesian coordinate system, with the xây plane coinciding with the cis surface of the slab.
The xrRNA structure was represented with SMOG. The model retains the full atomistic details and has been validated for various molecular systems, including nucleic acids21,22,23,24. SMOG employs a native-centric force field that allows us to highlight structure-dependent properties while filtering out details more related to sequence or specific RNAâprotein interactions.
We carried out hundreds of ramped translocation simulations of the molecule after priming it at the pore entrance from both ends, Fig. 2a. A representative set of translocation trajectories from \({5}^{\prime}\) and \({3}^{\prime}\) ends is shown in Fig. 2b. The curves show the temporal trace of the translocated fraction of the molecule at the reference temperature Tâ=â300âK and force-ramping rate r0â=â270âpN/μs. Translocation initiates only after a long loading stage, after which it proceeds precipitously. Major differences in translocation resistance emerge between the \({5}^{\prime}\) and \({3}^{\prime}\) ends.
The unlocking, or triggering of translocation, occurs when the vertical position of the leading P atom, measured relative to the cis plane, falls below zâ=ââ19.5âà or zâ=ââ15âà for \({5}^{\prime}\) or \({3}^{\prime}\) entries, respectively, and is followed by a cascade of ruptured contacts without significant pauses between the breaking of various secondary elements (Supplementary Figs. 2, 3). The typical order of rupture is given Fig. 2c. For the inherent stochasticity of the process, the triggering force is not unique but varies in a range, following a unimodal probability distribution as shown in Fig. 2b. Compared with the \({3}^{\prime}\) case, the most probable triggering force at the \({5}^{\prime}\) end is more than double; it is also much larger than the standard deviation of the triggering forces at the same end, 8âpN, or at the opposite one, 17âpN.
Comparison with xrRNA stretching response
The unusually high resistance at the \({5}^{\prime}\) end is further highlighted by contrasting the translocation forces at this end with those required to unfold the xrRNA by pulling the two termini apart, as in force spectroscopy. As shown in Fig. 2b, unfolding by stretching, which occurs when the longitudinal end-to-end distance, \({R}_{ee}^{| | }\) exceeds 16ânm, requires forces three times smaller than those needed to initiate \({5}^{\prime}\) translocations. Interestingly, the typical sequence of secondary elements rupture by stretching is similar to the case of \({3}^{\prime}\) translocations. The main difference is that the P4 helix, which yields when \({R}_{ee}^{| | } \,\sim\)â9ânm, now unfolds well before the rest of the structure (Fig. 2b and Supplementary Fig. 2).
In summary, at similar conditions, the characteristic force needed to initiate xrRNA translocation from the \({5}^{\prime}\) end, which resists exonuclease degradation, is distinctly higher than the forces that causes unfolding by stretching or by translocation from the \({3}^{\prime}\), which is processable by replicases. The same results are seen for the several hundreds of trajectories gathered over 19 combinations of T and r, from lower to higher temperatures and faster or slower ramping protocols than the reference values, see Supplementary Figs. 4, 7.
BellâEvans analysis
The extensive set of translocations at all explored (T, r) combinations can be recapitulated by the BellâEvans (BE) analysis. In force-spectroscopy contexts it is customary to resort to the following BE expression for the most probable rupture force17,18,25:
where ν is a kinetic coefficient and ET and Î, the key quantities of interest, are the effective height and width of the thermodynamic barrier associated to the rupture event at zero force.
Figure 3a shows that expression (1) provides a good fit of the most probable translocation forces, yielding a Ï2â~â3, hence of order unity, for both the \({5}^{\prime}\) and \({3}^{\prime}\) ends when all 19 (T, r) combinations are considered simultaneously. For clarity, only a subset of the considered datapoints are shown in Fig. 3âsee Supplementary Fig. 7 for the complete setâand the best fit to the BE expression is shown for Tâ=â300âK only.
The best-fit parameters are \({E}_{T}^{{\rm{BE}}}\,=\,36\,\pm\, 2\)âkcalâmolâ1 and ÎBEâ=â4.4â±â0.4âà for \({5}^{\prime}\) entries and \({E}_{T}^{{\rm{BE}}}\,=\,24\,\pm\, 1\)âkcalâmolâ1 and ÎBEâ=â2.0â±â0.1âà for \({3}^{\prime}\) ones. At room temperature, Tâ=â300âK, these barriers are approximately equal to 60âKBT and 41âKBT for the \({5}^{\prime}\) and \({3}^{\prime}\) ends, respectively.
These barriers describe the system in the absence of a driving force and reflect that translocation is very unlikely to be activated by thermal fluctuations only, especially at the \({5}^{\prime}\) end. In the presence of a constant driving force, an estimate of the activation time can be obtained via the so-called BellâEvans lifetimes18 and the best-fit parameterization of Eq. (1), see âMethodsâ. For constant driving forces of 50â100âpN, a range accessible to molecular motors26, the estimated activation time from the \({5}^{\prime}\) end exceeds by more than two orders of magnitude the one at the \({3}^{\prime}\) end.
Directional resistance and stress redistribution
The precipitous unlocking of translocation occurs when the P atom at the \({5}^{\prime}\) end is just about at the trans edge of the pore, zâ~ââ19.5âà (Supplementary Fig. 3). This triggering event is preceded by a long loading phase during which the driving force, which is exclusively imparted to the P atoms inside the pore, propagates an increasing tensile disturbance to the xrRNA that is still on the cis side.
The ensuing effects are best discussed by considering how the network of native contacts is distorted during the loading stage. The latter can be quantified in terms of the strain, s, of the native interactions for each nucleotide. Negative or positive values of s indicate that contacting atoms are respectively drawn closer or further apart than in the unstrained xrRNA structure, see âMethodsâ. Kymoplots of the s profile, during the loading stage of translocation from both ends, are shown in Fig. 4a and Supplementary Fig. 6.
The data and accompanying structures of Fig. 4 clarify that, as the initial transmission of tension along the RNA backbone draws the first nucleotides closer to or into the pore, their network of interacting nucleotides is increasingly distorted. In fact, regions G21-U28 and G38-A45 wrap more tightly around the \({5}^{\prime}\) end, particularly nucleotide G3 (Fig. 4b). Thanks to this tightening, which has no analog at the \({3}^{\prime}\) end, the increasing load is distributed over a broad portion of the molecule, which can hence better withstand the dragging force.
Increasing the load at the \({5}^{\prime}\) end eventually weakens the native interactions in several regions. These include xrRNA portions directly exposed to the traction of the \({5}^{\prime}\) end (U4-G14) and their partner strands (A17-C22 and C47-A52), see Fig. 4c. Note that the distal region A45-C48 is further disrupted by being pressed against the outer surface of the pore (Fig. 4c and Supplementary Fig. 5). At the \({3}^{\prime}\) end, instead, no noticeable tightening occurs around the pulled portion. The terminal tract of the P4 helix and their interacting partners in Pk4 and P4 (G31 and C58-C60) are, instead, directly exposed to the unmitigated pulling into the pore. (Fig. 4a, d). The lack of a mechanism for distributing and dissipating stress is what causes the \({3}^{\prime}\) end, the one processed by replicases, to yield at much lower forces than the \({5}^{\prime}\) end, the one resisting to exonucleases.
The tightened portion at the \({5}^{\prime}\) end, highlighted in blue in Fig. 4, is disconnected along the sequence, but is structurally compact. The emergence of a tightened structural element in response to a pulling force, is a distinctive feature of systems with mechanical tensegrity. In these systems, mechanical resistance results from an equilibrium between elements put under tension (positive strain) and other elements that are put under compression (negative strain) by the same applied force. Thus, the \({5}^{\prime}\) architecture provides a heretofore rare example of intra-molecular tensegrity, where the very pulling of the \({5}^{\prime}\) end is what boosts resistance to translocation and consequently to degradation.
Native contacts network and \({5}^{\prime}\) resistance
We next investigated which set of native contacts and structural elements at the \({5}^{\prime}\) end are most conducive to the redistribution of tensile stress and translocation resistance. Our model allows for directly pinpointing these critical contributors, namely by switching off selected combinations of native secondary contacts from the interaction potential. In experiments this is not possible, and similar information can be obtained by performing mutation scans where the identity of nucleotides is affected.
Our analysis indicates that the observed resistance is virtually ascribable to the sole contacts in pseudoknots Pk1 and Pk2. These, we recall, include contacts formed by sugar and phosphate moieties, Fig. 1c. The results are illustrated in Fig. 3b as boxâwhisker plots for the critical translocation forces of the native xrRNA chain, and three variants lacking the native contacts in either of the two pseudoknots, ÎPk1 and ÎPk2, or both, ÎPk1âPk2.
In variant ÎPk1 the critical force is appreciably lowered, but it still is almost twice as large as for the \({3}^{\prime}\) end. A stronger variation is seen for ÎPk2, whose critical force is reduced so much that it overlaps with the critical forces for the \({3}^{\prime}\) end. A complete loss of directional resistance is found for ÎPk1âPk2, which yields practically the same forces encountered at the \({3}^{\prime}\) end. This implies a fallback to a baseline hindrance once Pk1 and Pk2 are disrupted. Of the two sets of pseudoknotted contacts, those in Pk2 are noticeably the most important ones.
\({5}^{\prime}\) resistance, free-energy profile and transition state
We finally used the metadynamics framework27 to characterize \({5}^{\prime}\) resistance in more detail and beyond the BE analysis of out-of-equilibrium trajectories. In metadynamics, reversible transitions across significant barriers are achieved with an affordable computational expenditure by using history-dependent biases on multiple order parameters at the same time. Building on the centrality of Pk1 and Pk2 contacts on translocation resistance (Fig. 3b) we used the corresponding fraction of native contacts, QPk1 and QPk2, as natural order parameters, together with the amount of pore insertion of the \({5}^{\prime}\) P atom, z. To avoid effects related to the refolding of the RNA portion that has passed to the trans side, we extended the length of the pore and profiled the free energy up to zâ=ââ22Ã .
A first result of the metadynamics analysis is the dominant translocation pathway shown in Fig. 5a. The pathway projections in the (QPk1, z) and (QPk2, z) planes clarify that during the initial stages of translocation, z from 0 to â10âà , contacts in Pk1 and Pk2 are only modestly disrupted. Further insertion into the pore, z from â10âà to â15âà , causes an appreciable loss of Pk2 contacts. The final insertion tract, z from â15âà to â19.5âà , produces a concurrent loss of contacts in both Pk1 and Pk2. Representative structures along the pathway to the force-induced transition state are shown next to the pathway, along with two-dimensional projections of the free-energy profile at the same z values, Fig. 5bâd.
The results clarify that Pk2 interactions yield at zâ~ââ14âà , when the contacts involving the sugar and phosphate moieties of nucleotides G3-C22 are sheared and disrupted by the drawing of G3 into the pore while C22 is sustained and held in place by the elaborate secondary and tertiary architecture surrounding the \({5}^{\prime}\) end. The disruption of the nonlocal contact G3-C22 opens the way for the disruption of other contacts in Pk2 and also Pk1, such as G3 with C44 and A45, for zâ~ââ18âà .
We note that the dominant pathway established from metadynamics analysis has direct bearings on out-of-equilibrium, force-ramped translocations too. This is clarified by the scattered datapoints in Fig. 5a which represent the same order parameters for conformations sampled from the force-ramped translocation trajectories at Tâ=â300âK and râ=âr0. The agreement between the dominant pathway and the trajectories has three implications: it corroborates both methods, it indicates that translocation-induced unfolding occurs via a single pathway, and shows that the latter is the same in and out-of-equilibrium. The results thus cover slower ramping rates than considered in Fig. 3, including those accessible experimentally.
The salient thermodynamics of the process is aptly conveyed by the one-dimensional free-energy profile of Fig. 5e, obtained by thermodynamic integration of the metadynamics free-energy at fixed z. The free-energy increase up to zââââ10âà reflects the loss of conformational entropy and native contacts during the loading stage. The steeper rise beyond zââââ10âà reflects the progressive disruption of Pk2 and Pk1 contacts. The free-energy cost for reaching the depth zâ=ââ19.5âà , the triggering condition in force-ramped simulations, is about 46âkcalâmolâ1.
Applying constant translocation forces tilts the free-energy profile as shown in Fig. 5e. One notes that at sufficiently high force the slope of the free-energy changes sign at about zâ~ââ20âà . This value of z pleasingly matches the aforementioned triggering insertion depth obtained in force-ramp translocations (Supplementary Fig. 3). The pore insertion depth zâ~ââ20âà can be interpreted as an effective transition state, embodying the rate-limiting step for force-induced translocations.
In summary, the metadynamics results aptly complement those from force-ramped translocations, giving detailed indications of the free-energy barrier associated with the triggering of translocation from the \({5}^{\prime}\) end, and of the associated pathway. The combined analysis shows that, in the loading stage, regions G21-U28 and G38-A45, which surround the \({5}^{\prime}\) end, become more tightly wrapped around it. This mechanism, that has no analog at the \({3}^{\prime}\) end, stabilizes the network of long-range native contacts around nucleotide G3, hindering their entrance into pore until contacts in Pk1 and Pk2, including G3-C22, are eventually disrupted.
Discussion
Nanopore translocation is a powerful method to probe biopolymers as well as understanding how they interact with processive enzymes8,15,24,28,29,30,31,32,33. Here we used it to clarify the structural determinants of xrRNAs resistance to exonuclease degradation. To this end we carried out hundreds of simulations where an atomistic, native-centric model of Zika xrRNA was translocated through narrow a pore from either end. The setup provides a transparent abstraction of the biological process where the \({5}^{\prime}\) and \({3}^{\prime}\) xrRNA termini are engaged and driven through the lumen of enzymes, such as nucleases or polymerases, encountering different directional resistance.
First, we observed that translocation resistance is significantly higher at the \({5}^{\prime}\) end, the one resisting exonuclease degradation, than at the \({3}^{\prime}\) one, the one yielding to the unwinding operated by replicases and reverse transcriptases. This holds systematically throughout the many considered combinations of temperatures and force-ramping protocols. It should not go unnoticed that, in our ramping simulations, \({3}^{\prime}\)-entry translocations not only initiate, but also complete at forces well below those needed to start translocation from the \({5}^{\prime}\) end. The xrRNA resistance to unfold when pulled through a pore by the \({5}^{\prime}\) end was also shown to be significantly higher than to unfold by pulling away the two ends (stretching). BellâEvans analysis indicates that, in the absence of a driving force, the translocation free-energy barriers at the two ends differ by about 12âkcalâmolâ1. At constant forces of 50â100âpN, the estimated difference of activation times exceeds two orders of magnitude, and ought to be detectable experimentally. Incidentally, we note that the directional effect, as well as the different response to translocation and stretching, suggest that experiments where the two ends are pulled apart might lead to results that do not report directly on the resistance offered to electrophoretic translocation and processive molecular motors.
Secondly, the directional effect originates from the different ways in which the pulling effects are propagated from the termini to the rest of the molecule. During the early loading stages, pulling the \({5}^{\prime}\) end causes regions G21-U28 and G38-A45 to wrap more tightly around it. At later loading stages, interactions between the pulled terminal tract U4-G14 and partner strands A17-C22 and A45-A52 start to yield. At the same time, region A45-C48 also distorts due to the tight contact with the pore surface. Both mechanisms reflect in a compression of the region surrounding the pulled \({5}^{\prime}\) terminal; this region is compact in space but not in sequence. The tension-induced compression, and the ensuing resistance, make the effect entirely analogous to that of systems with mechanical tensegrity, which oppose tensile forces with compressive ones elicited by the applied tension itself.
As in tensegrity systems, the resistance of the \({5}^{\prime}\) end is encoded in the architectural organization of the molecule and how it interacts with the pore. By contrast, no such mechanism is available at the \({3}^{\prime}\) end for its simpler architecture. By selectively removing native secondary interactions, we established that \({5}^{\prime}\) translocation resistance is primarily encoded in pseudoknot Pk2, and secondarily in Pk1. Once Pk1 and Pk2 contacts are disrupted, no other outstanding barriers are encountered at the \({5}^{\prime}\) end, and the resistance becomes comparable to that of the \({3}^{\prime}\) end. The latter fact also rules out the helical asymmetry of nucleic acids, see e.g., refs. 23,34, as a cause for xrRNAs directional resistance.
Lastly, analysis of free-energy calculations and out-of-equilibrium simulations are consistent at indicating that translocation and the concomitant unfolding from the \({5}^{\prime}\) end occur via a single dominant pathway. Inspection of the pathway shows that a critical step is the disruption of the interactions of nucleotides G3 and C22, which unlocks the breaking of the other Pk2 and Pk1 contacts eventually triggering the irreversible translocation. Also, comparison with metadynamics results indicates that free-energy barriers estimated with customary BellâEvans analysis are susceptible to simplifying assumptions on the barrier widths. This finding might be relevant in the analysis of force-ramp experiments too.
In conclusion, these results demonstrate that the structural architecture of xrRNA, which is the sole feature informing our model, is predisposed to generate a profoundly asymmetric resistance to translocation from the two ends. This occurs because mechanical stress is distributed and accumulated in the network of contacts in a way that depends not only on the end that is directly pulled, but also on the regions interacting with it, that become more tightly wound around it. This pulling-induced tightening protects the \({5}^{\prime}\) end and allows it to better withstand further pulling of the end itself.
Note that this resistance mechanism is altogether different from the catch-bonding behavior35 of knotted molecules translocating through narrow pores, which is qualitatively opposite to what observed here. Knotted filaments, particularly those with twist knots, can translocate at sufficiently low forces but jam at high ones36. For such systems, jamming happens when the knot, which is pulled against the pore entrance, becomes tight enough that any stochastic sliding motion of the chain along its knotted contour is suppressed9,36. In xrRNA, instead, translocation is stalled at low forces but is unjammed at sufficiently high ones, and moreover in a directional-dependent manner. For xrRNA the initial resistance is offered by the strained network of native interactions that compactify the molecule around the region at the pore entrance, but that eventually yields ad sufficiently high forces.
The results help rationalize and recapitulate in a single structural and quantitative framework not only xrRNA resistance to \({5}^{\prime}\) degradation from various exonucleases5, but also the concomitant much lower hindrance to being fully processable from the \({3}^{\prime}\) in order to be copied by virally encoded RNA-dependent RNA polymerases. Both properties thus find a natural explanation as different manifestations of the same architecture-dependent response. The results also indicate that present-day experimental setups for pore translocation ought to be adequate to capture and reveal the discussed signature effects of the xrRNA directional resistance. In perspective, the structural organization principles that Zika xrRNAs have acquired by evolutionary selection could be co-opted for designing synthetic molecules or meta-materials37 with directional resistance to translocation and context-dependent tensegrity response.
Methods
Zika xrRNA and system setup
The reference, native conformation of Zika xrRNA was taken from the protein databank38 entry 5TPY3. Simulations of the 71-nucleotide long RNA (1527 heavy atoms) were based on SMOG, a native-centric model retaining the full atomistic structural details39,40. The SMOG potential includes terms for bonded, nonbonded, angular, and dihedral interactions. The model is designed for exposing structurally-encoded properties by neglecting effects, such as transient formation of non-native contacts and details of interaction with water and ions. These effects are not expected to be crucial in this specific mechanistic application where, besides, the faster ramping protocols called by explicit solvent treatments would antagonize directional resistance, Fig. 3a. Default values for the relative strength of the potentials and the cutoff contact map were used40,41. The absolute temperature was calibrated similarly to ref. 23 by matching the experimental condition that mechanical unfolding of RNA helices42,43 occur at about 15âpN at 300âK, see Supplementary Methods and Supplementary Table 1.
The P atom at the \({5}^{\prime}\) or \({3}^{\prime}\) end was primed at the entrance of a cylindrical pore and prevented from retracting from it by a ratchet-like restraint. The pore is 11.7âà wide and transverses perpendicularly a 19.5âà thick slab, which separates the cis and trans sides. The pore geometry approximates the lumen of the Xrn1 exoribonucleases, PDB entry 2Y3519. The same SMOG inter-atom excluded volume potential was used for steric interactions of the RNA atoms with the surface of the slab and pore.
Molecular dynamics simulations
Langevin MD simulations were carried out with the LAMMPS package44, with proper atomic masses, default values for the friction coefficient and with integration timestep equal to 8.7âÃâ10â4ÏMD, where ÏMD is the characteristic MD timescale, see Supplementary Methods. Chain translocation from the cis to the trans sides are driven by a force, F, parallel to the pore axis, and that is exclusively applied to the P atoms inside the pore, among which it is equally subdivided at any given time. The driving force is ramped linearly with time t, Fâ=âr â t. In stretching simulations, the \({5}^{\prime}\) and \({3}^{\prime}\) P atoms are pulled in opposite directions with the same loading protocol. We considered seven equispaced temperature in the range 281â319âK, and as many pulling protocols, r/r0â=â0.3, 1, 5, 10, 20, 50, 100, with r0â=â6.2âÃâ10â4âpN/ÏMD. Conventional mapping to physical units yields ÏMDâ~â2.3âps, see Supplementary Methods, which yields r0â=â270âpN/μs, which is conservative for molecular dynamics simulations though still faster than experiments. Note, however, due to the lack of explicit solvent in the model, the effective simulated timescale (pulling rate) reported here represents a lower (upper) bound estimate45, implying a better agreement with experimental pulling rates. A total of 19 combinations of parameters were explored, see Supplementary Table 2, and for any of them 20 translocation simulations were typically collected for each orientation (\({5}^{\prime}\to {3}^{\prime}\) and vice versa). At the reference temperature and ramping rate (Tâ=â300K, râ=âr0) the number of timesteps needed to trigger translocation from the \({5}^{\prime}\) end is about 7âÃâ108. Including stretching simulations, more than 1000 trajectories were run on the Ulysses supercomputer at SISSA.
BellâEvans analysis
For each orientation, the most probable translocation forces at the 19 combinations of temperature and loading rate, F*(T, r), were simultaneously fitted with the BellâEvans expression of Eq. (1). The fit consisted of a Ï2 minimization over the free parameters ET, Î and ν. F*(T, r) was established from a Gaussian kernel-density analysis of translocation triggering events, corresponding to a pore insertion depth of the terminal P atom equal to zâ=ââ19.5âà and â15âà for \({5}^{\prime}\) and \({3}^{\prime}\) entries, respectively. The estimated error on F*(T, r) was set equal to the error of the mean translocation force. The BE waiting time for triggering translocation at a fixed force, F, is Ïâ=âFβ/rF, where Fβâ=âKB T/Î and rF is the loading rate for which F*â=âF18.
Strain analysis
We considered the relative strain of a nucleotide, sâ=â(q0âââqF)/q0. q0 is the average contact fraction (overlap) of all native interactions of the nucleotide computed at a reference range of low forces, 0â10âpN. qF is analogously computed at force F. The contact fraction of two atoms at distance d is \(q(d)\,=\,1/(1\,+\,{[d/(1.5\,\cdot\, {d}^{0})]}^{6})\), d0 being their native distance (default cutoff distance of 4âà ).
Metadynamics
Metadynamics simulations27 were carried out with the COLVAR module46 using three collective variables: z, QPk1 and QPk2, which are respectively the depth of insertion of the \({5}^{\prime}\) P atom into the pore, and the fraction of native contacts in pseudoknots Pk1 and Pk2. QPk was computed as \(\frac{1}{{n}_{c}}\mathop{\sum }\nolimits_{j \,=\, 1}^{{n}_{c}}q({d}_{j})\), where j runs over the nc native contacts of the pseudoknot. Unbiased simulations, see Supplementary Fig. 8, were used to set the widths of the metadynamics multivariate Gaussians (Ïzâ=â0.25âà , \({\sigma }_{{Q}_{{\rm{Pk1}}}}\,=\,{\sigma }_{{Q}_{{\rm{Pk1}}}}\,=\,0.0125\)), their height (0.15âkcalâmolâ1) and deposition interval (103 timesteps). Retraction of chain from the pore was prevented by a restraining potential, \(V\,=\,\frac{k}{2}\,\cdot\, {z}^{2}\,\cdot\, \theta (-z)\), where kâ=â15.6âkcalâmolâ1âà â2 and θ is the Heaviside function. For the metadynamics analysis, the pore was extended so to reconstructs the free-energy landscape up to zâ=ââ22âà , which required runs of 1.5âÃâ109 timesteps, Supplementary Fig. 9. A constant force bias of 110âpN was used to reduce the time required for filling the basin. The bias was then discounted with a suitable reweighting a posteriori. The dominant pathway was obtained by discretizing z with intervals of width Îzâ=â0.8âà and computing the corresponding canonical averages of QPK1 and QPK2.
Reporting summary
Further information on research design is available in the Nature Research Reporting Summary linked to this article.
Data availability
Data supporting the findings of this paper are available from the corresponding author upon reasonable request. A reporting summary for this Article is available as a Supplementary Information file. The source data underlying Figs. 1c, 2b, 3a, b, 4a, 5 are provided as Source data file. Source data are provided with this paper.
Code availability
Typical setup files for the LAMMPS simulations are provided as Supplementary Software 1. Setup files for SMOG simulations, natively created for GROMACS, were converted for LAMMPS using custom extensions and conversion scripts. The latter are made freely available for academic use at the "SMOG-converter" GitHub repository, github.com/CristianMicheletti/SMOG-converter. Source data are provided with this paper.
References
Pijlman, G. P. et al. A highly structured, nuclease-resistant, noncoding RNA produced by flaviviruses is required for pathogenicity. Cell Host Microbe 4, 579â591 (2008).
Chapman, E. G. et al. The structural basis of pathogenic subgenomic flavivirus RNA (sfRNA) production. Science 344, 307â310 (2014).
Akiyama, B. M. et al. Zika virus produces noncoding RNAs using a multi-pseudoknot structure that confounds a cellular exonuclease. Science 354, 1148â1152 (2016).
Slonchak, A. & Khromykh, A. A. Subgenomic flaviviral RNAs: what do we know after the first decade of research. Antivir. Res. 159, 13â25 (2018).
MacFadden, A. et al. Mechanism and structural diversity of exoribonuclease-resistant RNA structures in flaviviral RNAs. Nat. Commun. 9, 119 (2018).
Luo, K., Ala-Nissila, T., Ying, S.-C. & Bhattacharya, A. Sequence dependence of DNA translocation through a nanopore. Phys. Rev. Lett. 100, 058101 (2008).
Chen, G., Chang, K.-Y., Chou, M.-Y., Bustamante, C. & Tinoco, I. Triplex structures in an RNA pseudoknot enhance mechanical stability and increase efficiency ofâ1 ribosomal frameshifting. Proc. Natl Acad. Sci. U.S.A. 106, 12706â12711 (2009).
Comer, J., Dimitrov, V., Zhao, Q., Timp, G. & Aksimentiev, A. Microscopic mechanics of hairpin DNA translocation through synthetic nanopores. Biophys. J. 96, 593â608 (2009).
Suma, A., Rosa, A. & Micheletti, C. Pore translocation of knotted polymer chains: how friction depends on knot complexity. ACS Macro Lett. 4, 1420â1424 (2015).
Ziegler, F. et al. Knotting and unknotting of a protein in single molecule experiments. Proc. Natl Acad. Sci. U.S.A. 113, 7533â7538 (2016).
Plesa, C. et al. Direct observation of DNA knots using a solid-state nanopore. Nat. Nanotechnol. 11, 1093â1097 (2016).
Suma, A. & Micheletti, C. Pore translocation of knotted DNA rings. Proc. Natl Acad. Sci. U.S.A. 114, E2991âE2997 (2017).
San MartÃn, Ã. et al. Knots can impair protein degradation by atp-dependent proteases. Proc. Natl Acad. Sci. U.S.A. 114, 9864â9869 (2017).
Sriramoju, M. K., Chen, Y., Lee, Y.-T. C. & Hsu, S.-T. D. Topologically knotted deubiquitinases exhibit unprecedented mechanostability to withstand the proteolysis by an aaa+ protease. Sci. Rep. 8, 7076 (2018).
Cao, C. et al. Mapping the sensing spots of aerolysin for single oligonucleotides analysis. Nat. Commun. 9, 1â9 (2018).
Sivertsson, E. M., Jackson, S. E. & Itzhaki, L. S. The aaa+ protease clpxp can easily degrade a 3 1 and a 5 2-knotted protein. Sci. Rep. 9, 1â14 (2019).
Dudko, O. K., Hummer, G. & Szabo, A. Theory, analysis, and interpretation of single-molecule force spectroscopy experiments. Proc. Natl Acad. Sci. U.S.A. 105, 15755â15760 (2008).
Evans, E. Probing the relation between force-lifetime-and chemistry in single molecular bonds. Annu. Rev. Biophys. 30, 105â128 (2001).
Jinek, M., Coyle, S. M. & Doudna, J. A. Coupled 5â nucleotide recognition and processivity in xrn1-mediated mRNA decay. Mol. Cell 41, 600â608 (2011).
Kowalczyk, S. W., Wells, D. B., Aksimentiev, A. & Dekker, C. Slowing down DNA translocation through a nanopore in lithium chloride. Nano Lett. 12, 1038â1044 (2012).
Whitford, P. C. & Sanbonmatsu, K. Y. Simulating movement of tRNA through the ribosome during hybrid-state formation. J. Chem. Phys. 139, 121919 (2013).
Nguyen, K. & Whitford, P. C. Steric interactions lead to collective tilting motion in the ribosome during mRNAâtRNA translocation. Nat. Commun. 7, 1â9 (2016).
Colizzi, F. et al. Asymmetric base-pair opening drives helicase unwinding dynamics. Proc. Natl Acad. Sci. U.S.A. 116, 22471â22477 (2019).
Bandarkar, P., Yang, H., Henley, R. Y., Wanunu, M. & Whitford, P. C. How nanopore translocation experiments can measure RNA unfolding. Biophys. J. 118, 1612â1220 (2020).
Björnham, O. & Andersson, M. Theory for nonlinear dynamic force spectroscopy. Eur. Biophys. J. 46, 225â233 (2017).
Smith, D. E. et al. The bacteriophage Ï29 portal motor can package DNA against a large internal force. Nature 413, 748â752 (2001).
Laio, A. & Parrinello, M. Escaping free-energy minima. Proc. Natl Acad. Sci. U.S.A. 99, 12562â12566 (2002).
Kantor, Y. & Kardar, M. Anomalous dynamics of forced translocation. Phys. Rev. E 69, 021806 (2004).
Branton, D. et al. The potential and challenges of nanopore sequencing. Nat. Biotechnol. 26, 1146â1153 (2008).
Henley, R. Y., Carson, S. & Wanunu, M. Studies of RNA sequence and structure using nanopores. Prog. Mol. Biol. Transl. Sci. 139, 73â99 (2016).
Meller, A., Nivon, L. & Branton, D. Voltage-driven DNA translocations through a nanopore. Phys. Rev. Lett. 86, 3435 (2001).
Trepagnier, E. H., Radenovic, A., Sivak, D., Geissler, P. & Liphardt, J. Controlling DNA capture and propagation through artificial nanopores. Nano Lett. 7, 2824â2830 (2007).
Cao, C., Yu, J., Wang, Y.-Q., Ying, Y.-L. & Long, Y.-T. Driven translocation of polynucleotides through an aerolysin nanopore. Anal. Chem. 88, 5046â5049 (2016).
Mathé, J., Aksimentiev, A., Nelson, D. R., Schulten, K. & Meller, A. Orientation discrimination of single-stranded DNA inside the α-hemolysin membrane channel. Proc. Natl. Acad. Sci. U.S.A. 102, 12377â12382 (2005).
Pereverzev, Y. V., Prezhdo, O. V., Forero, M., Sokurenko, E. V. & Thomas, W. E. The two-pathway model for the catch-slip transition in biological adhesion. Biophys. J. 89, 1446â1454 (2005).
Rosa, A., Di Ventra, M. & Micheletti, C. Topological jamming of spontaneously knotted polyelectrolyte chains driven through a nanopore. Phys. Rev. Lett. 109, 118301 (2012).
Paulose, J., Chen, B. G.-g & Vitelli, V. Topological modes bound to dislocations in mechanical metamaterials. Nat. Phys. 11, 153â156 (2015).
Berman, H. M., Bourne, P. E., Westbrook, J. & Zardecki, C. The protein data bank. In Protein Structure, 394â410 (CRC Press, 2003).
Whitford, P. C. et al. An all-atom structure-based potential for proteins: bridging minimal models with all-atom empirical forcefields. Proteins 75, 430â441 (2009).
Noel, J. K. et al. Smog 2: a versatile software package for generating structure-based models. PLOS Comput. Biol. 12, e1004794 (2016).
Hayes, R. L. et al. Reduced model captures Mg2+-RNA interaction free energy of riboswitches. Biophys. J. 106, 1508â1519 (2014).
Manosas, M., Collin, D. & Ritort, F. Force-dependent fragility in RNA hairpins. Phys. Rev. Lett. 96, 218301 (2006).
Wen, J.-D. et al. Force unfolding kinetics of RNA using optical tweezers. Biophys. J. 92, 2996â3009 (2007).
Plimpton, S. Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 117, 1â19 (1995).
Yang, H. et al. Diffusion of tRNA inside the ribosome is position-dependent. J. Chem. Phys. 151, 085102 (2019).
Fiorin, G., Klein, M. L. & Hénin, J. Using collective variables to drive molecular dynamics simulations. Mol. Phys. 111, 3345â3362 (2013).
Bottaro, S. et al. Barnaba: software for analysis of nucleic acid structures and trajectories. RNA 25, 219â231 (2019).
Humphrey, W., Dalke, A. & Schulten, K. Visual molecular dynamics. J. Mol. Graph. 14, 33â38 (1996).
Acknowledgements
This work was partially supported by MIUR, the Italian Ministry of Education. We thank Antonio De Simone for useful discussions.
Author information
Authors and Affiliations
Contributions
A.S., L.C., G.B., and C.M. designed research; A.S., L.C., and C.M. carried out the simulations; A.S. and C.M did the analysis and wrote the paper with inputs from G.B.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information Nature Communications thanks Aleksei Aksimentiev, and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Source data
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Suma, A., Coronel, L., Bussi, G. et al. Directional translocation resistance of Zika xrRNA. Nat Commun 11, 3749 (2020). https://doi.org/10.1038/s41467-020-17508-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-020-17508-7
This article is cited by
-
Driven translocation of a semiflexible polymer through a conical channel in the presence of attractive surface interactions
Scientific Reports (2022)
-
Enhancement of prime editing via xrRNA motif-joined pegRNA
Nature Communications (2022)
-
Mechanical strength of RNA knot in Zika virus protects against cellular defenses
Nature Chemical Biology (2021)
-
Pseudoknot length modulates the folding, conformational dynamics, and robustness of Xrn1 resistance of flaviviral xrRNAs
Nature Communications (2021)
-
Molecular mechanisms underlying the extreme mechanical anisotropy of the flaviviral exoribonuclease-resistant RNAs (xrRNAs)
Nature Communications (2020)