David
Fischermeier
a,
Christian
Steinmetzger
b,
Claudia
Höbartner
b and
Roland
Mitrić
*a
aInstitut für Physikalische und Theoretische Chemie, Universität Würzburg, Am Hubland, 97074 Würzburg, Germany. E-mail: roland.mitric@uni-wuerzburg.de
bInstitut für Organische Chemie, Universität Würzburg, Am Hubland, 97074 Würzburg, Germany
First published on 24th November 2023
Förster resonant energy transfer (FRET) can be utilized in the study of tertiary structures of RNA aptamers, which bind specific fluorophoric ligands to form a fluorogenic aptamer complex. By introducing the emissive nucleobase analog 4-cyanoindole into the fluorogenic Chili RNA aptamer a FRET pair was established. The interpretation of studies aiming to investigate those tertiary structures using FRET, however, relies on prior knowledge about conformational properties of the nucleobase, which govern exciton transfer capabilities. Herein we employed classical molecular dynamics combined with Förster exciton theory to elucidate the preferred orientation relative to proximate bases and the influence on exciton transfer efficiency in multiple substitution sites. We did this by comparing the chromophoric distances emergent from MD simulations with experimental FRET data based on structural data of the native aptamer. We present the outlined methodology as a means to reliably evaluate future nucleobase analogue candidates in terms of their structural behavior and emergent exciton transfer properties as exemplified in the study of the preferred orientation of 4-cyanoindole in the Chili RNA aptamer.
The exploration of FRET between nucleobase-modified RNA and bound ligands is a new addition to fluorogenic aptamers that was demonstrated for the Chili RNA aptamer, which activates fluorescence of DMHBI+ and DMHBO+. The universal nucleobase analogue 4-cyanoindole (CI) was identified as a suitable FRET partner of the Chili ligands and used to map their binding site within the RNA aptamer.17,25 However, quantitative interpretation of the results was hampered by the lack of information on the relative orientation of CI and the ligand.
Recently, the structure of the Chili aptamer in complex with two cationic ligands, DMHBI+ and DMHBO+,15,16 was solved by X-ray crystallography26 (Fig. 1). The RNA forms a central two-tiered G-quadruplex as the ligand-binding element, which is anchored by a basal stem and an apical stem-loop. A long-range base pair helps to establish a preorganized proton transfer pathway to N7 of G15, which enables efficient ligand-to-RNA excited state proton transfer (ESPT) in less than 130 fs and causes fluorescence emission with a large Stokes shift. Site-specific modification of Chili with the universal fluorescent ribonucleoside analog r4CI was used to engineer a supramolecular base-to-ligand FRET system. By incorporating r4CI at eight positions (G6, A11, G17, U23, U34, G36, G45, G46) based on previously reported Tb3+ probing and mutagenesis experiments,15,16 the exciton transfer efficiency was mapped for the different domains of the Chili aptamer earlier experimentally.17
The goal of this study is assessing the effect of those substitutions on the structural integrity of the wild-type aptamer as well as the orientational preferences of r4CI. The orientation of FRET paired bases is vital to the exciton transfer process and thus the interpretation of experimental FRET data. Due to a lack of knowledge regarding the orientation of the fluorescent base one often assumes isotropic behavior, which can lead to errors of orders of magnitude when trying to assess the pair distance from FRET experiments as already shown for DNA analogues.27,28 This orientation may deviate from the orientation of the unsubstituted base in the native aptamer's crystal structure due to differences in binding affinities between original base and substitute. Insight into the orientations of the base can be gained from MD simulations. While this study focuses heavily on the Chili aptamer complex, the outlined methodology is not restricted to this field of research and is widely usable in any context where FRET experiments are utilized to gain insights into the structure of large systems.
Using classical MD simulations we were able to gain some insights into the influence of the newly introduced base on pairing, stacking and the general structure of the immediate neighbouring bases. We were also able to analyze the preferred orientations of the substituted bases and their effect on the exciton transfer and quenching. To this end, these classical MD simulations were combined with quantum chemical TD-DFT calculations to theoretically study the exciton transfer based on the classical Förster model.
The outlined methodology can help scan a plethora of different bases and select for stability of the tertiary structure as well as the desired orientation patterns in synthetic nucleoside strands, which allows for a faster and more accurate interpretation of future experimental studies and selection of modification sites and bases. As the employed method only relies on a wild-type crystal structure, the often prohibitively arduous – or in case of in vivo FRET studies impossible – process of isolating and crystallizing a multitude of mutations can be partially or fully replaced by these simulations.
κ = A·D − 3(A·)(D·), | (1) |
κ = cosθD–A − 3cosθA–RcosθD–R, | (2) |
(3) |
(4) |
Fig. 2 Depiction of the angular relation between the transition dipole moment (D and A) and distance vectors (). |
To assess the stability of the system the root mean squared deviation (RMSD) of the atomic positions for parts and the entirety of the system have been calculated according to
(5) |
To evaluate the the orientation factor κ, the orientations of the correct transition dipole moments are paramount. Fig. 3 shows the orientation of the transition dipole moments of both chromophores' and experimental data for comparison. For CI the transition dipole calculated by TD-DFT is in excellent agreement with experimental data.49 While the DMHBO+ transition dipole moment has not been determined yet, several similar species are already studied, which show similar angles (70–100°) between the CO bond and the transition dipole moment50 suggesting that both transition dipole vectors are accurately simulated. The FRET donor r4CI was individually introduced into the wild-type Chili RNA aptamer at multiple sites (G6, A11, G17, U23, U34, G36, G45, G46). However, not all orientations of these bases, e.g. stacked, unstacked and different rotational degrees of freedom, are guaranteed to be visited within a single, limited trajectory due to the energy barriers involved likely to lead to trapping a replica in a local minimum. Thus for every substitution several starting orientations were propagated. Unless not permitted by the position of the original base (e.g. sterical constraints) these starting orientations consisted of two varieties with r4CI occupying the intrahelical space of the original base, one of which features a 180° rotation (I1,2). Furthermore several orientations with extrahelical r4CI pointing outwards of the complex (E1–3, see ESI.† Fig. S1–S35) were also simulated.
To assess the movement and stability of the system the root mean squared deviation (RMSD) of the atomic positions was calculated for the structures within the trajectories for the complete systems as well as the chromophores as shown in Fig. 4. The U23 species was chosen as a good baseline example, since U23 is located at the aptamer's loop and as such can't pair or stack as bases in the helix/tetraloop could and as such can reasonably only take a single orientation, where the cyanoindole can freely rotate around the glycosidic bond (see ESI.† Fig. S36–S43 for RMSD data of the remaining species). To exclude the influence of rotational and translational motion of the complete system, the structures of the complete system were aligned with each other along the backbone of the aptamer. The U23 E1 modified system is equilibrated within roughly 1 ns, after which the RMSD of atomic positions of the system lies between 3 and 5 Å relative to the starting orientation. The atomic positions of DMHBO+ fluctuate roughly around 2 Å when aligning the complete system. An illustration of the positional changes of the complete system and DMHBO+ is shown in Fig. 5. This positional rigidity implies a strong interaction between DMHBO+ and the aptamer. r4CI shows the highest positional flexibility with an RMSD between 2 and 7 Å. In the same way the rotation and translation of the system were excluded by aligning its structures, the change in the local geometry of the chromophores can be characterized by aligning the chromophore's structures. CI shows nearly no change in geometry, while DMHBO+ also shows very little change in structure along the trajectory. This shows, that no rotation of the rings is observed for DMHBO+. The RMSD data for r4CI and the aligned CI show, that CI is very rigid as expected, but can rotate rather freely around the glycosidic bond for the U23 species. This is also true for all sterically unrestricted orientations.
From MD simulations the chromophore distance RDA (MD) is readily available for each orientation of all substituted species. Using the exciton model this chromophoric distance (RDA (Exc.)) can also be calculated with eqn (4) from κ, which is obtainable from eqn (2). If the chromophoric distance RDA (Exc.) based on experimental data and the orientation factor is sufficiently congruent with RDA (MD), it suggests that the orientation and chromophoric distance observed during the MD simulation is indeed dominant in solution. Fig. 6 shows a comparison of the donor–acceptor distance, RDA (MD), over the course of the MD simulation as well as RDA (Exc.), which was obtained from experimental data17 shown in Table 1 in conjunction with the MD simulations. We show the data for the trajectories, which yield the best agreement for the chromophoric distances as well as other noteworthy behavior affecting those. A table with all values for all trajectories can be found in the ESI.† We note here, that for all species multiple starting orientations give good agreements for the RDA values as they cross over and converge geometrically over the course of the trajectory, which is also shown and discussed in the ESI.†
Fig. 6 Comparison of the donor–acceptor distance, RDA, over the course of the MD simulation (MD, blue) as well as calculated based on the exciton model (exc.,red) and experimental values shown in Table 1. For visibility the 10 value running average for the exciton model is shown due to the high fluctuations. For G45 I2 the angle between the donor and acceptor transition dipoles (θD–A) is also displayed in yellow. |
Species | Quantum yield | E τ | J | R 0 (Å) | κ 2 | R DA (Å) | |
---|---|---|---|---|---|---|---|
MD | Exc. | ||||||
a 1014 l nm4 cm−1 mol−1. | |||||||
G6 E3 | 0.48 | 0.36 | 4.40 | 29.7 | 0.13 | 25.0 | 28.2 |
A11 E2 | 0.25 | 0.67 | 4.90 | 22.5 | 0.09 | 18.1 | 20.0 |
G17 E3 | 0.07 | 0.86 | 4.10 | 21.0 | 0.65 | 20.6 | 19.4 |
U23 E1 | 0.49 | 0.80 | 4.80 | 35.8 | 0.50 | 31.5 | 28.1 |
U34 E2 | 0.15 | 0.70 | 4.50 | 22.0 | 0.15 | 21.1 | 19.5 |
G36 E3 | 0.28 | 0.82 | 4.40 | 24.0 | 0.95 | 22.0 | 28.2 |
G45 E3 | 0.10 | 0.37 | 4.60 | 15.4 | 0.15 | 22.5 | 23.2 |
G46 I1 | 0.28 | 0.62 | 4.80 | 18.9 | 0.03 | 17.1 | 18.3 |
Overall we were able to identify at least one trajectory for every mutation, which produced chromophoric distances very similar to the distances calculated with the exciton model based on the MD and experimental spectroscopic data. A selection of those will be discussed herein, while the datasets for all simulations are shown in the ESI† (Fig. S44–S50). Some external orientations for each species did cross over into each other at some points, for example G6 E2/E3. The G46 mutation was the only species, for which the orientation located inside the helix (G46 I1) produced the best agreement. This confirms the conclusions found by thermal melting studies17 suggesting that substituted sites possess more conformational freedom compared to their unsubstituted variant and that the modified aptamer's r4CI exhibits a significant unstacked population. The I2 orientation of U34 also produced very good agreement for the majority of the trajectory until the orientation changes substantially due to movement in the RNA backbone, which leads to a divergence between the RDA values. This suggests that its dominant orientation might be facing inward or an equilibrium of I2 and E2. The G45 I2 orientation also produces two close RDA curves initially, but diverges rapidly due to multiple subtle geometry changes in the aptamer complex, which result in a steady change of θD–A effecting the change in RDA. This, however, could indicate that r4CI in the G45 species might also be in an inwards facing orientation or in an equilibrium of G45 I2 and G45 E3. Notably the U34, G45 and G46 moieties are very close in proximity at the very bottom of the quadruplex. Thus it seems that this region, where the helix opens up to transition into the quadruplex, allows r4CI to pair with neighbouring bases possibly due to the lack of sterical constraint.
The G17 E3 and U23 E1 species exhibit orientation factors close to the isotropic limit of 2/3. In Fig. 7 we show that θD–A, which should show high fluctuations for isotropic behavior, is nearly constant for U23 E1 and in most other trajectories. Although r4CI exhibits considerable movement when located outside of the helix as seen in Fig. 4, the most freely accessible movement for any unpaired r4CI nucleotide is the rotation around the glycosidic bond. Fig. 3 shows the almost parallel alignment of r4CI's transition dipole moment with the glycosidic bond, which results in the invariance of θD–A in regards to this rotation. Thus the chromophores' relative transition dipole orientation is mostly constant and the value for the orientation factor is not a result of isotropic behavior, but rather a rather rigid orientation resulting in the same value as the isotropic limit.
Fig. 7 θ D–A over the course the selected MD trajectories. The constant behavior shows a rigid orientation of the transition dipole moment. |
For the G17 E3 species the angular relation between the transition dipole moments shows different behavior with a large change at roughly 5 ns suggesting a more flexible environment. In fact the same trend can be observed in Fig. 6, where the agreement of both RDA values is very close before 5 ns with discrepancies arising afterwards. This shift in behavior was found to be a result of a sort of “flipping” motion in the RNA backbone at the G17 position (see Fig. 8). Whereas r4CI “leans” on the helix at the beginning of the simulation, it unlatches from that position and the backbone rotates r4CI further away from the helix tilting the loop portion of the complex during the motion. Following this motion the congruence in RDA fluctuates, but still shows very good agreement at times. This could suggest, that the flipping motion might represent real behavior of the aptamer complex. Some other sites like G6 and G45, where the backbone seems less rigid, also undergo frequent, albeit mostly minor changes in orientation, which can, however, result in large spikes for the calculated RDA.
Fig. 8 The initial leaning orientation of r4CI (red) in the G17 position starting from the G17 E3 orientation is illustrated in (a), while (b) shows r4CI facing away from the complex. |
While we have shown that the chromophores exhibit a rather rigid geometry due to their extensive aromatic, planar rings, these geometric changes might still affect the transition dipole vectors, which would affect the orientation factors and thus fail to produce accurate values for RDA. To assess whether our method of constructing the transition dipole vectors from transition charges is sufficiently accurate, the transition dipole vectors of both chromophores were directly calculated using TD-DFT in GAUSSIAN16 along the trajectory of the U23 species. Furthermore the chromophore distance was calculated subsequently based on those vectors as seen in Fig. 6.
The average chromophoric distances of U23 with 25.2 Å for DFT and 28.1 Å for the exciton method show a good agreement and both curves appear similar in shape especially after the first 4 ns, showing, that both react very similarly to changes in relative orientations as seen in Fig. 6. This shows that the dipole moments constructed using atomic transition charges are in good agreement with the transition dipole moments obtained directly from DFT and thus the exciton model derived transition dipole vectors and orientation factors should give an accurate description of the exciton transfer process.
Even localized outside of the aptamer complex the angular relation of the r4CI and DMHBO+ is similarly constant as it is in an orientation facing inwards to the complex. This suggests that the isotropic limit of κ2 = 2/3 should not be assumed when working with r4CI as nucleoside analog.
Studying which characteristics allow for and enforce intrahelical, well defined extrahelical orientations or enable free movement could prove important to finding suitable sites for base substitution in the future.
The method employed was shown to be comparatively accurate in obtaining transition dipoles and orientation factors to full DFT at a fraction of the computational cost and produced results widely in agreement with prior experimental data. Obtaining crystals of modified biomacromolecules is often virtually impossible due to the prohibitive effort required especially while maintaining the precise tertiary structure present in solution. The outlined methodology could be widely employed to gain very detailed insight into the topology of a variety of artificial nucleobases and their influence on the tertiary structural stability as well as any other chromophoric biomolecules used in FRET applications. Provided a crystal structure for the unmodified system alongside FRET data is provided, both FRET chromophores can be studied to obtain properties like pairing, stacking or true isotropic behavior.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3cp04704k |
This journal is © the Owner Societies 2024 |