C. M.
Grytz
a,
A.
Marko
a,
P.
Cekan
b,
S. Th.
Sigurdsson
*b and
T. F.
Prisner
*a
aInstitute of Physical and Theoretical Chemistry and Center of Biomolecular Magnetic Resonance, Goethe University, 60438 Frankfurt am Main, Max-von-Laue-Str. 9, Hessen, Germany. E-mail: prisner@chemie.uni-frankfurt.de
bUniversity of Iceland, Department of Chemistry, Science Institute, Dunhaga 3, 107 Reykjavik, Iceland. E-mail: snorrisi@hi.is
First published on 18th December 2015
The cocaine aptamer is a DNA three-way junction that binds cocaine at its helical junction. We studied the global conformation and overall flexibility of the aptamer in the absence and presence of cocaine by pulsed electron–electron double resonance (PELDOR) spectroscopy, also called double electron–electron resonance (DEER). The rigid nitroxide spin label Ç was incorporated pairwise into two helices of the aptamer. Multi-frequency 2D PELDOR experiments allow the determination of the mutual orientation and the distances between two Çs. Since Ç is rigidly attached to double-stranded DNA, it directly reports on the aptamer dynamics. The cocaine-bound and the non-bound states could be differentiated by their conformational flexibility, which decreases upon binding to cocaine. We observed a small change in the width and mean value of the distance distribution between the two spin labels upon cocaine binding. Further structural insights were obtained by investigating the relative orientation between the two spin-labeled stems of the aptamer. We determined the bend angle between this two stems. By combining the orientation information with a priori knowledge about the secondary structure of the aptamer, we obtained a molecular model describing the global folding and flexibility of the cocaine aptamer.
The DNA-based cocaine aptamer was discovered and subsequently engineered as an aptamer-based fluorescent sensor for cocaine by Stojanovic and coworkers (Fig. 1B).10 This particular construct was designed to be a single strand in the absence of the ligand. The aptamer folds up into a three-way junction, as detected by fluorescence spectroscopy, upon binding to cocaine11 with the dissociation constant in the micromolar range.11–13 We have previously incorporated the rigid cytidine–analogue, spin-labeled nitroxide Ç (Fig. 1A), into single sites of this cocaine aptamer construct and studied its folding upon cocaine binding.12 Since Ç becomes fluorescent upon reduction,14–16 it offered a unique opportunity to study the folding of this aptamer by two spectroscopic techniques, continuous-wave (CW) electron paramagnetic resonance (EPR) and fluorescence spectroscopies. The interpretation of our data was consistent with a model in which the aptamer was partially folded up into two co-axially stacked helices in the absence of cocaine (helices II and III, Fig. 1B). Our data indicated that the co-axially stacked helices would tilt, relative to each other, upon concomitant binding to cocaine and formation of helix I. More recently, Johnson and coworkers described a stem-length dependent binding mechanism of cocaine to the aptamer using a combination of nuclear magnetic resonance (NMR) spectroscopy and isothermal titration calorimetry.13 In contrast to our previous results,12 they did not observe a pre-formed helix III with the construct having a short helix I (Fig. 1B), while all three helices were formed when helix I was extended to seven base-pairs. For all the constructs that they investigated, binding to cocaine led to a well-defined bound-state structure with all three helices formed. However, these NMR studies were performed in the absence of salts, which affects single-strand to duplex equilibria.
PELDOR (pulsed electron–electron double resonance),17–19 also called DEER (double electron–electron resonance), is an EPR technique that has successfully been used to study the structure and conformational changes of biopolymers.20,21 By recording the PELDOR signal, which oscillates with the frequency νdip of the dipole–dipole interaction between two electron spin centers, one can determine their distance R
R = 3.736·νdip−1/3|3·cos2Θ − 1|1/3, | (1) |
The spin label Ç is an ideal probe for PELDOR studies since it contains a nitroxide fused to the nucleobase14,16 and, therefore, fixes the location and orientation of Ç relative to the framework of the nucleic acid helix14 unlike other PELDOR studies of DNA using flexible spin lables.22–25 Consequentially, a well-defined distance between such a pair of spin labels and their mutual orientations can be extracted from orientation-selective 2D PELDOR experiments (evolution time t vs. pump–probe frequency offset Δν = νA − νB), as we have recently shown for double-stranded DNA molecules.26 The frequency offset averaged PELDOR data reflect only the conformational flexibility of the DNA molecule itself without the need to consider any conformational flexibility of the spin label which is negligibly small in comparison to the DNA mobility.27 Therefore, the orientation between two spin-labeled double-stranded DNA helices can also be extracted from such 2D PELDOR data sets, allowing determination of the global folding and tertiary structure of nucleic acid molecules containing secondary structure elements, such as stems, bulges and hairpin loops.
In this work, we utilize the spin label Ç and 2D PELDOR to investigate conformational changes of the cocaine aptamer. Both constructs (with a short or elongated helix I) undergo a cocaine-induced conformational change, which presumably originates at the junction where cocaine is known to bind.10,11 This is directly visible in the PELDOR time traces by more pronounced oscillations and a small distance change, from the offset-averaged PELDOR time trace applied Tikhonov regularization. A detailed simulation of the 2D PELDOR time traces revealed that there is considerable conformational flexibility in the helical junction, even after binding to cocaine. In most cases so far, PELDOR spectroscopy has been used as a supporting technique to obtain long-range distance restraints for structural studies of large biomolecules. This study underlines the potential of using rigid spin labels to determine directly the flexibility, conformational changes of oligonucleotide molecules and the relative orientations of secondary structural elements which is useful for structure determination.
Two peaks can be observed in the distance distributions obtained by Tikhonov regularization of aptamer 1 in the absence of cocaine (Fig. 3B, black curve), whereas all other samples show a single peak. These two peaks might correspond to two different conformations: the major peak is at 2.9 nm, while another broad and weak peak appears at around 4.8 nm. The observation of two peaks in the distance distributions would be consistent with a previously proposed two-state model for folding of the cocaine aptamer consisting of T- and Y-shaped junctions,12 that are expected to give different distances. However, the reliability of distance distributions depends strongly on the maximum dipolar evolution time tmax.30 Due to this reason, we recorded time traces with an extended dipolar evolution time up to 5 μs with samples of aptamer 1 in the presence and the absence of cocaine in a deuterated solvent. These measurements were performed at the Q-band to avoid artifacts induced by a poor signal-to-noise ratio. Measurements at Q-band frequencies have higher sensitivity compared to X-band and a somewhat reduced orientation selectivity. Fig. S2 in the ESI,† shows that we again observe a second peak in the distance distribution of aptamer 1 without cocaine. The PELDOR data for the cocaine-bound state display also a second peak, but with considerably reduced intensity. A model distribution consisting of one and two Gaussian peaks were applied to fit the averaged PELDOR time traces. Two Gaussians have to be used to fit the averaged PELDOR time trace of aptamer 1 in the absence of cocaine which also supports a second conformer with a larger distance (data not shown). Nevertheless because of the low intensity of the second broad peak the assignment to the proposed two-state model remains somewhat speculative.
Another small but significant effect can be observed upon binding of cocaine for both constructs 1 and 2, namely a shift ΔR of the mean distance of about 0.3 nm, indicating a conformational change upon binding to cocaine. This shift of the maximum peak of the distance distribution is clearly visible for both aptamer samples and for the PELDOR experiments performed at X- and Q-band frequencies in non-deuterated as well as deuterated solvent (Fig. S2 and S11, ESI†). For flexible spin labels such a shift could also be induced by changes of the rotameric state of the spin label itself but that can be excluded for the rigid spin label Ç, again demonstrating the ideal properties of this spin label with respect to quantitative distance determination. However, the change in distance can also be due to orientation selection.
The dependence of the oscillation frequency on the angle Θ can easily be visualized by the dipolar oscillation formula, whereas the dependence of the modulation depth on the specific geometry and offset frequency requires a quantitative modelling of pumped and probed (A and B) spins for a given spin pair geometry. The modulation depth parameter is also sensitive to experimental parameters such as the MW power, resonator bandwidth, pulse length and B1 homogeneity over the sample length. We have shown that the modulation depth is reproducible within 5–10% error from experiment to experiment, if the experiments are performed under identical conditions.32 Within a 2D PELDOR data set, where each offset measurement is performed immediately after the other, this error will be even lower. Thus, the different variations of the modulation depth as a function of the offset frequency in the experiments on aptamers 1 and 2, with and without cocaine, (Fig. 2) reflect different conformational distributions of relative orientations between the two spin labels of the different samples, rather than experimental errors.
Fitting of the experimental 2D PELDOR data set of the cocaine-bound aptamer 1, the spin label conformers of the PELDOR database library reproduced all PELDOR time traces very well including the modulation depth (Fig. 4A and Fig. S5, ESI†). In total, 100 iterative fitting steps where performed yielding 53 non-redundant spin-label conformers for the cocaine-bound state of this aptamer (Fig. 4B). A clustering of all spin label conformers in the β angle range 60–90° is obvious (as mentioned above β1 and β2 can be interchanged). The same procedure has been applied to the 2D PELDOR data set of aptamer 1 in the non-cocaine bound state. This data show no clear offset dependence and indeed, as can be expected, a much broader orientation distribution results from the fit procedure (see the ESI† for more details about the fit procedure especially with respect to the modulation depth). This confirms that aptamer 1 shows high conformational flexibility in the absence of the ligand. In contrast, the cocaine-bound state of the aptamer is more ordered.
As explained above, a priori knowledge about the secondary structure of the three-way junction was subsequently used to construct physically possible solutions for the relative orientation between helix II and III from this angular information. The procedure yielded an ensemble of 1125 possible structural models from 549504 possible helix orientations that are in agreement with the experimental PELDOR data (Fig. S7, ESI†).
For visualization of the extracted relative conformational flexibility between aptamer helices II and III, their geometry was described by the model sketched in Fig. 5A, defined by the bend angle φ (angle between both helix axes) and the two twist angles ψ and θ (rotation of the two nitroxide N–O vectors xx,II and xx,III of the respective helix axis). The bend angle can be directly extracted from the angles β1 and β2, whereas the twist angles ψ and θ are only restricted by the geometrical restraints imposed by the a priori knowledge from the secondary structure of the aptamer.
Fig. 5 (A) Definition of bend φ and twist angles (θ and ψ). Both helices are represented as cylinders. More details can be found in the ESI.† (B) A density plot of twist angle θ against bend angle φ for the conformers found from the simulation of the 2D PELDOR data set for aptamer 1 in the cocaine-bound state. (C) A density plot of twist angle ψ against twist angle θ. Color bar specifies the density of the conformers to the solution. The data points are interpolated. |
As can be seen for the cocaine-bound aptamer 1 the bend angle is very well conserved within the conformational ensemble. The analysis gives a weighted mean bend angle of 86 ± 19°. The twist angles θ and ψ scatter much more but with some correlation between them. The molecular models using the helices as rigid bodies can be found in the ESI† (Fig. S9).
The geometrical parameters resulting from the algorithm fit to the PELDOR data might not be unique, especially for rather flexible systems, where different ensembles might lead to almost identical good fits.33 To test the robustness of our fit we removed the area of the spin label conformational space where the solutions were concentrated (Fig. 6B, blue squared region) from the library and tried to fit the experimental PELDOR data set by such a restricted library. As can be seen in Fig. 6A, the fit gets significantly worse. Interestingly, the new conformers found from this fit border on the edges of the cut-region (Fig. 6B) approaching the missing structures as much as possible. These results confirm that our found spin label conformations from the PELDOR database fit are essential and not a solution replaceable by another ensemble of spin label conformations.
The rigid spin label Ç, used in this EPR study, contains a nitroxide rigidly fused to the nucleobase (Fig. 1A). Therefore, Ç is sensitive to the motion of the base and thereby to helical motions when the spin label is paired with guanine in a duplex region.14 Thus, in addition to accurate distance measurements between two Ç-labels, the damping of the dipolar oscillations in the PELDOR time traces give detailed information on the conformational flexibility of the DNA molecule.27 We have previously shown that PELDOR signals from duplex DNA that contains two Çs, show well resolved dipolar oscillations (depending on the labelling position and the mutual orientation of the two spin labels).26 Furthermore, the oscillation frequency, modulation depth and damping of the PELDOR signals vary strongly as a function of the probe frequency offset Δν.29 In comparison, the oscillations in the PELDOR signals for both aptamers 1 and 2 are more strongly damped (Fig. 2 and 3) and show much less variation in the oscillation frequency as a function of the frequency offset (Fig. 2), indicating that both have considerable conformational flexibility, albeit less so in the cocaine-bound state. In summary, these data show that the cocaine aptamer can be considered as a fairly flexible junction, especially in the absence of cocaine.
Our previous studies of the cocaine aptamer using fluorescence and CW-EPR spectroscopy indicated that helices II and III are already formed in the absence of cocaine.12 Furthermore, upon addition of cocaine an increase in fluorescence of Çf, located in helix II immediately at the junction (Fig. 1B), was taken as evidence for helical stacking of aptamer 1 before addition of cocaine.12 A NMR study has also been performed by Johnson and coworkers on a cocaine aptamer that was constructed similar to aptamer 1 of this study, in which helix I contained three base pairs and helix III was a four base-pair stem loop. According to their data, “only a small amount of secondary structure” was formed in the absence of cocaine.13 However, the NMR measurements were performed in the absence of salt, which reduces duplex stability, and may explain the lack of secondary structure. The PELDOR studies of aptamers 1 and 2, performed in the presence of salt, show a well-defined distance at 2.9 nm without cocaine (Fig. 3B), confirming that helix III already forms in the absence of cocaine. The additional peak at higher distances could reflect the presence of two conformations; speculatively T- and Y-shaped junctions12 with a difference in the spin–spin distance between the two spin labels of around 1.9 nm.
The averaged time traces and distance distributions are almost identical for aptamers 1 and 2 in the presence of cocaine, which shows that elongation of helix I in aptamer 2 does not influence the structure of the cocaine-bound state (Fig. 3 and Fig. S1, ESI†). This is consistent with the previous results that have shown that cocaine stabilizes the formation of helix I and the aptamer.10,12,13
The mean distance between the spin labels in both aptamer constructs increased by about 0.3 nm upon addition of cocaine. This small change in the distance is at the edge of what is detectable by PELDOR. However, since this shift is observed for both aptamer samples and in all our control PELDOR experiments at X- and Q-band frequencies (Fig. S2, ESI†), we judge it to be significant. It is doubtful that this relatively short change in the distance between the two spin labels could have been determined using flexible spin labels, although small differences in the distance have been reported.40 This underscores the advantage of using rigid spin labels such as Ç. This change in the inter-spin distance indicates a cocaine-induced conformational change, which presumably originates at the junction where cocaine is known to bind.10,11 The approximation to average out the orientation information is qualitatively very good.27 Nevertheless, the small shift observed upon cocaine binding is at the limit of resolution of the PELOR method.
Aptamer 1 consists of a 13-mer with the sequence 5′-d(GACAAGAAGÇGAC) which is paired with a 24-mer with the sequence 5′-d(GTCGCTTCTTCAAÇGAAGTGGGTC). Aptamer 2 contains a longer helix I (seven base pairs instead of three, Fig. 1B). A GT-mismatch in helix III of the original construct12 was substituted by a Watson–Crick base pair GÇ which results in a tighter binding of cocaine.13 Helix II, which is known to form in the absence of cocaine,12 was elongated to enable spin labels to be incorporated into helix II further away from the junction. A two-strand construct was used to give more flexibility in preparation of spin-labeled oligonucleotides and higher yields.
The spin-labeled oligonucleotides were synthesized as previously described.14,16 10 nmol of spin-labeled DNA strands were annealed in 100 μl of 10 mM phosphate buffer (pH 7), 100 mM NaCl and 0.1 mM EDTA. Aptamer samples with the cocaine ligand (5 mM) were prepared in a similar manner. After evaporation of the solvent, each sample was dissolved in 100 μl of 20% non-deuterated or deuterated ethylene glycol/water (H2O or D2O) mixture. Ethylene glycol was used as a cryoprotectant. Before the PELDOR experiments at X-band frequencies, 25 μl of the spin-labeled samples were filled in EPR tubes (2.8 mm OD) and frozen in liquid nitrogen.
π/2(νA) − τ1 − π(νA) − (τ1 + t) − π(νB) − (τ2 − t) − π(νA) − τ2 − echo signal |
For the samples with non-deuterated solvent, the detection pulses (νA) were set to 32 ns for both the π and the π/2 pulses. The pump pulse (νB) had a length of 12 ns and was located at the maximum of the anisotropic powder spectrum selecting the nitrogen spin mI = 0 sublevel with random orientations of the nitroxide radical with respect to the external magnetic field (Fig. 7c). This frequency was set as the central resonator frequency. The amplitude of the pump pulse was optimized by the inversion of a Hahn-echo at the frequency νB. The frequency offset Δν = νA − νB was varied between 40 and 90 MHz in steps of 10 MHz, selecting different orientations of the nitroxide radicals with respect to the external magnetic field through the anisotropy of the nitrogen hyperfine tensor A. Therefore, the probe pulses (νA) were slightly off-resonant, but still within the bandwidth of the microwave resonator. Proton modulation was suppressed by averaging of 8 spectra of variable τ1 ranging from τ1 = 132 to 196 ns in steps Δτ1 of 8 ns. The time delay between the second and the third detection pulse (τ1 + τ2) was set to 2.5 μs. In order to eliminate receiver offsets, the π/2 pulse was phase-cycled (0°, 180°). A video amplifier bandwidth of 25 MHz was chosen. The spectra were recorded with a repetition time of the pulse sequence of 4000 μs or 5000 μs. Depending on the frequency offset Δν, the approximate accumulation time for each time trace was usually 2–3 h (for Δν = 40 − 70 MHz) and about 8–15 h (for Δν = 80 and 90 MHz).
Quantitative analysis of the orientation dependency of the PELDOR time traces recorded at X-band frequencies was performed using a numerical algorithm to globally fit the 2D PELDOR data set; the algorithm was developed by our group for the application on rigid spin labels.27 In the first step, a library of all 2D PELDOR data sets that can occur at X-band frequencies for nitroxide radical pairs with arbitrary orientations and inter-spin distances in the experimentally accessible range was created. In the second step, an iterative fitting was applied to find an optimal combination of these pre-simulated PELDOR time traces to reproduce the experimental data set. The fitting algorithm searches for an ensemble of spin label conformers that minimize the mean squared deviation between the experimental and set of library data sets:
∥n−1(Sr,1 + … + Sr,n−1 + Sr,n) − Sexp∥2 → min | (2) |
At X-band frequencies, only the z-component of the nitrogen hyperfine tensor is resolved as shown in Fig. 7C. Due this reason, only the angles β1, α2, β2 and the distance R (Fig. 7B) have to be varied to create the library of data sets needed for the fitting procedure.
The γ1 and γ2 angles cannot be determined from X-band PELDOR measurements, because of the similar values of Axx and Ayy such that the asymmetry in x and y directions cannot be determined (Fig. 7C). It should be noted however that γ1 and γ2 are required for unique description of the mutual orientation between both Ç labels. This means that while it is possible to obtain the bend angle, φ, between both helices for the X-band 2D PELDOR datasets, the twist angle, θ, of one helix relative to the other cannot be determined. Therefore, it is not possible to define a unique solution from our experimental data. This problem could in principle be solved by performing orientation dependent PELDOR measurements at higher magnetic fields, where the x and y components of the anisotropic g-tensor are spectrally resolved. To achieve this, a frequency of 95 GHz or above would be needed. However, unfortunately the inherent flexibility of the aptamer would lead to a very small modulation depth at this frequency due to variations in the Lamor frequency of the spin pairs. This made these experiments unfeasible. Instead, we used geometrical knowledge about the helices to restrict the conformational space and provide further refinement of the simulation.
Two B-DNA helices were constructed using the web interface software w3dna.46 For helix II the sequence 5′d-(AAGAAGCGAC) paired with 5′d-(GTCGCTTCTT) and for helix III the sequence 5′-d(GTTG) was paired with 5′-d(CAAC). The GA mismatch in helix III was not considered. The rigid nitroxide spin label Ç was modeled into both DNA helices. The coordinates of all atoms were loaded into a custom written Matlab program, which orients the DNA according to the nitroxide spin label orientation found by the fitting algorithm (for nitroxide 1 the Euler angle β1 and for nitroxide 2 the Euler angles α2, β2) and the varied γ1 and γ2 angles. The second Ç labeled helix was additionally translated by the distance R between the two electron spins each located in the middle of the NO-bond. In addition, because the two nitroxide spin labels cannot be distinguished, we considered all possible solutions for spin label 1 in helix II and spin label 2 in helix II and vice versa. Similar, we considered that the nitroxide radical normal (zz) of each spin can point up or down, which turns the helix upside down. Within the analysis, the Ç spin labeled DNA helices were treated as rigid bodies. With the iterative fit algorithm, the non-redundant parameter sets (α2, β1, β2, R), corresponding to spin label conformers were found to best fit the experimental 2D PELDOR data set. From this data set all possible helix orientations were calculated, which yielded all possible helix orientations.
Fortunately, not all of these structure models for the cocaine-bound state are compatible with the basic geometrical knowledge of the aptamer structure and can therefore be discarded. All structures were filtered by three conditions, which are based on the secondary structure of the cocaine aptamer. The criterion for choosing the threshold conditions was that the experimental PELDOR data could still be satisfactorily simulated with the remaining bundle of conformers. From the secondary structure of the aptamer it is known that helices II and III have to be covalently bound. One can imagine that the junction might stretch the pitch height. Therefore, we chose a threshold value that allows stretching of the helix by maximum 0.8 nm. Structural models with a separation distance between the P-atom at the 5′ end of Cytosin24 and 3′OH at Thymin23 longer than 0.8 nm were filtered out. Consideration of the maximum and the minimum space required by the non-labeled helix I (diameter of a B-DNA helix 2 nm) completed the geometry of the three-way junction. All structures with a separation between the 3′OH at the 3′ end of helix III and the P-atom at the 5′ end of helix II (deoxyadenosine 14), beyond 1.5 nm and 2.7 nm were filtered out. Moreover, steric clashes between helices II and III were prevented. The plugin, “Show bumps”, was implemented in PyMOL to calculate the van-der Waals repulsion between the helices (http://www.pymolwiki.org/index.php/Show_bumps). Again all structural models with steric conflicts caused by the repulsion between the two helices were discarded. Using these three filter threshold conditions reduces the number of possible structures.
This study is the first example of applying PELDOR with rigid spin labels to study the changes in the structure and conformational flexibility of a nucleic acid junction upon ligand binding. The results clearly show the advantage of using rigid spin labels to investigate such systems, since the flexibility of the spin label itself does not need to be taken into account, allowing more accurate distance determination and unravelling additional information about the relative orientation of the two spin labels with respect to each other. This yielded conclusive information on the structure and flexibility of an aptamer from a single pair of spin label positions. In spite of no available atomic resolution structure of the investigated cocaine aptamer, a molecular model describing the global folding and flexibility of the cocaine aptamer using the helices as rigid bodies was obtained. This work lays the foundation to develop an efficient EPR structure calculation program for large nucleic acids or nucleic acid–protein complexes.
Footnote |
† Electronic supplementary information (ESI) available: X-and Q-band data for aptamer 1 in deuterated solvent, details about the analysis, structural models, and twist and bend angle calculations. See DOI: 10.1039/c5cp06158j |
This journal is © the Owner Societies 2016 |