Ying
He‡
a,
Yanbin
Zhang‡
b,
Lukasz
Wojtas
a,
Novruz G.
Akhmedov
c,
David
Thai
a,
Heng
Wang
a,
Xiaopeng
Li
a,
Hao
Guo
*b and
Xiaodong
Shi
*a
aDepartment of Chemistry, University of South Florida, 4202 E. Fowler Avenue, Tampa, Florida 33620, USA. E-mail: xmshi@usf.edu
bDepartment of Chemistry, Fudan University, 2005 Songhu Road, Shanghai, 200438, People's Republic of China
cDepartment of Chemistry, West Virginia University, Morgantown, WV 26505, USA
First published on 6th March 2019
Methanol soluble and stable guanosine octamers were successfully achieved via H-bond self-assembly. Through structural conformational design, we developed a new class of guanosine derivatives with modification on guanine (8-aryl) and ribose (2′,3′-isopropylidene). This unique design led to the formation of the first discrete G8-octamer with its structure characterized by single crystal X-ray diffraction, MS and NMR spectroscopy. The G8-octamer showed unique cation recognition properties, including the formation of a stable Rb+ templated G-quadruplex. Based on this observation, further modification on the 8-aryl moiety was performed to incorporate a cross-layer H-bond or covalent linkage. Similar G-octamers were obtained in both cases with structures confirmed by single crystal X-ray diffraction. Furthermore, the covalently linked G-quadruplex exhibited excellent stability even in MeOH and DMSO, suggesting a promising future for this new H-bond self-assembly system in biological and material applications.
A G-quartet is an interesting supramolecular scaffold formed by H-bonds.13,14 As shown in Scheme 1, with an approximately 90-degree angle between the H-bond donor and acceptor, four guanine units are held together to form a G-quartet. Through ion–dipole interactions, alkali and alkaline earth metal cations can enhance the process by serving as the template to coordinate with central oxygen atoms.15–17 Stacking of G-quartets gives G-quadruplexes as bioactive building blocks found in DNA and RNA folding.18,19 In this case, the extent of G-quartet stacking in a G-quadruplex will be determined by the phosphate backbone, which is often associated with the formation of a counter folding subunit, such as the i-motif.20
Inspired by this unique H-bond assembly, researchers have been devoted to developing guanine derivatives to achieve controllable G-quadruplex formation from small molecules.21–23 Some interesting applications have been identified with various G-quartet assemblies, including lipophilic ion channels,24–26 supramolecular hydrogels,27–29 nanomaterials,30,31 potential targets for cancer therapy,32,33 and more.34–37 Although many examples of G-quartet formation through various modified guanine derivatives have been reported, studies on controlling G-quartet stacking molecularity are relatively rare.38,39 Factors to be taken into consideration include the cation concentration,39 the solvent,40 the anion,41 and so on.42 In many cases, mixtures of various “stacking isomers” (G8, G12 or G16) were observed, which highlights the significant challenges associated with controlling the vertical stacking.43–46 Moreover, the assembly of a specific and stable G-quadruplex in H-bond competitive solvents remains a challenging task.48 Herein, we report the construction of the first G-octamer with structures characterized by single crystal X-ray diffraction through monomer conformational design. Moreover, with this new system, stable G-quadruplexes were formed with significantly improved stability. Through the design of cross-layer H-bonds and covalent linkage, G-octamers were prepared with excellent stability in MeOH (no dissociation) and even in 50% DMSO, which offers a potential opportunity to extend the H-bond assembly system into biosystems for future applications.
The simplest plausible G-quadruplex would be a G-octamer which is likely to adopt either top-to-top (T–T) or bottom-to-bottom (B–B) stacking patterns (Scheme 2A).50 In previous studies, Meijer et al. employed different concentrations of a guanosine derivative resulting in the formation of tunable G-octamers.41 Spada has reported G-octamer formation upon exposure to UV light through alkene isomerization.51 Wu presented several excellent examples of G-octamer formation through π–π interactions using NMR and MS studies.52 Hirao and coworkers applied Au(I)–Au(I) interaction between two G-quartet layers to achieve a G-octamer confirmed by NMR and CD spectra.53 However, to the best of our knowledge, no single crystal structure of a G-octamer has ever been reported, implying the challenging nature of preparing a stable and discrete G-octamer in a dynamic equilibrium.
To tackle this problem, we set out to design a G-derivative where the structure is rigid and predisposed for the conformation of a potential octamer structure so as to minimize the entropy cost involved in the self-assembly process. As shown in Scheme 2B, modification of guanosine often occurs on two positions: C-8 of purine and the hydroxyl of ribose. The Sessler group first reported on 8-aryl substituted guanosine in the formation of a G-quartet without templating cation in both solution and solid state.54
This seminal work initiated the concept of conformational control for G-quartet formation: the steric effect between the aryl substituent and protected ribose helped guanosine to adopt a syn conformation, preventing the ribbon formation and giving a tetramer as the dominant conformation. On the other hand, Davis's group developed lipophilic guanosine with bulky ribose to form a G-hexadecamer (Scheme 2B) both in solution and solid state (Mn+ = K+, Ba2+, Sr2+, and Pb2+).55 Notably, a picrate anion bridge played a crucial role: as revealed by single crystal X-ray diffraction (Scheme 2C), four picrate anions linked two G-octamers through H-bonds between anion and two inner G-quartets. The two G-octamers (from adjacent inner and outer G4) gave top-to-bottom stacking with ribose interdigitated between the adjacent layers. Interestingly, the two inner layers adopted bottom-to-bottom (B–B) stacking, which is more favorable than the T–T mode with the cation binding on the more “naked” convex face between the two layers. This result aroused our interest in developing a G-octamer through similar B–B stacking.
Considering the steric interaction between the C-8 substituent and ribose, we postulated that incorporation of C-8 aryl and the rigid ribose ring might provide a new system with steric hindrance between the G-quartet to force the G4 bowls to stack in a bottom-to-bottom manner, while obstructing ribose interdigitation at the top-face (Scheme 2C). To confirm this idea, compound 1b was designed, prepared and applied to assemble with various alkali and alkaline earth metal cations. 1H NMR spectra were obtained and selected regions of the 1H NMR spectra of these G-quadruplexes were compared with the G-hexadecamer from 1a as shown in Fig. 1.
Fig. 1 1H NMR spectra of G-quadruplexes (A) [(1a)16Ba2]4+·(Pic−)4; (B) [(1b)8Ba]2+·(Pic−)2; (C) [(1b)8K]+·(Pic−); (D) [(1b)8K]+·(PF6−); (E) [(1b)8Sr]2+·(Pic−)2 in CDCl3. |
As previously reported, treating 1a with alkali and alkaline earth metal salts (K+, Ba2+, Sr2+ and so on) gave two sets of signals in the 1H NMR spectra, corresponding to the inner and outer G-quartet.56 Conducting similar cation binding experiments with 1b in CDCl3 gave a single set of 1H NMR signals in all cases (M = K+, Ba2+, and Sr2+; A− = Picrate− or PF6−). Furthermore, ESI-MS demonstrated a clear doubly charged peak at m/z = 2123.01, corresponding to a mol. wt of 4246.68 for [(1b)8Ba]2+. The experimental and calculated isotope patterns further suggested an octameric composition. In addition, traveling wave ion mobility-mass spectrometry (TWIM-MS)57 confirmed no formation of stacking isomers, which excluded the formation of random aggregates in gas phase (see ESI† for details).
Finally, single crystal structures were obtained and unambiguously verified the G8-octamer formation with the proposed bottom-to-bottom stacking (Fig. 2A). The top view of the crystal structure (Fig. 2B) shows the G-quartet self-assembly in the tail to tail orientation. The five crystal structures obtained with ligand 1b include monovalent cations (K+ and Rb+) and divalent cations (Ba2+ and Sr2+). Picrate anion showed no clear binding with the G-quartet, consistent with what was observed in the 1H NMR spectra (see Fig. S2†). A complex with non-coordinated PF6− anion was also successfully obtained, [(1b)8K]+·(PF6−), confirming the “anion-free” binding mode of this new type of G8-quadruplex. The distances of the H-bond within the G4-quartet and between the G4 layers are compared in Table 1.
Fig. 2 X-ray single crystal structure of the G8-octamer (A) [(1b)8Rb]+·(Pic−); (B) [(1b)8K]+·(PF6−). Similar structures for [(1b)8K]+·(Pic−), [(1b)8Ba]2+·(Pic−)2 and [(1b)8Sr]2+·(Pic−)2 were also obtained.§ |
G-quadruplex | d(N1⋯O6) | d(N2⋯N7) | d(O–M) | d(G4–G4) | |
---|---|---|---|---|---|
a See ref. 46. | |||||
[(1a)16Ba2]4+·(Pic−)4a | Inner | 2.92 ± 0.01 | 2.91 ± 0.07 | 2.75 ± 0.02 | 3.06 (i–o) |
Outer | 2.86 ± 0.01 | 2.89 ± 0.04 | 2.79 ± 0.03 | 3.58 (i–i) | |
[(1b)8K]+·(Pic−) | 2.82 ± 0.08 | 2.87 ± 0.04 | 2.72 ± 0.17 | 2.85 | |
[(1b)8K]+·(PF6−) | 2.82 ± 0.04 | 2.89 ± 0.05 | 2.77 ± 0.03 | 2.96 | |
[(1b)8Ba]2+·(Pic−)2 | 2.89 ± 0.04 | 2.89 ± 0.04 | 2.72 ± 0.05 | 2.89 | |
[(1b)8Sr]2+·(Pic−)2 | 2.83 ± 0.03 | 2.86 ± 0.04 | 2.62 ± 0.05 | 2.75 | |
[(1b)8Rb]+·(Pic−) | 2.85 ± 0.03 | 2.89 ± 0.03 | 2.81 ± 0.05 | 3.04 |
According to these crystal structures, all G-octamers from 1b gave G-quartets with an (N1⋯O6) and (N2⋯N7) H-bond distance of around 2.9 Å, similar to that of the inner and outer layer in the G16-hexadecamer formed from 1a.46 These results indicated that both cations and the C-8 phenyl substituent had little influence on the H-bond in the G4-quartet. However, the size of the G8 was influenced by the average O–M distances of all these complexes from 1b, which might follow the trend where higher ionic potential (Z/r) resulted in shorter O–M distance (2.75–2.81 Å). Interestingly, in comparison with the G16-hexadecamer [(1a)16Ba2]4+·(Pic−)4, the G8-octamers [(1b)8Ba2]2+·(Pic−)2 gave a slightly shorter distance between the two G4 layers (2.89 Å vs. 3.06 Å and 3.58 Å). This result implied the stronger cation interaction of the (1b)4-quartet than the (1a)4-quartet. This improved cation interaction has been supported by G4-binding studies with Rb+ using the 1b ligand and led to the first crystal structures of Rb+ coordinated G-quadruplexes. Among all the G8 crystals, [(1b)8Rb]+·(Pic−) gave the longest O–M and G4–G4 distance due to its large radii58 and low valency. Very few examples of Rb coordinated G-quartets have been reported so far, indicating how challenging it is for guanosine to bind with Rb to form a discrete G-quadruplex.59 To the best of our knowledge, this is the first single crystal structure of a G-quadruplex containing Rb+, clearly suggesting the promising cation binding ability of guanosine derivative 1b.
Having successfully confirmed a new concise G8-quadruplex structure in solution (NMR), solid state (XRD) and gas phase (ESI-MS and TWIM-MS), we evaluated its stability in MeOH. Dissolving octamer [(1b)8K]+·(Pic−) in CD3OD gave a mixture of two sets of signals in NMR spectra, suggesting partial decomposition of this G8-octamer (vide infra). To obtain MeOH stable G-quadruplexes, further modification is still needed.
Notably, the Rivera group have reported the formation of an intralayer H-bond between carbonyl oxygen with N(2)H within the same G-quartet by using 2′-deoxy guanosine derivatives without rigid ribose functionalization.34,36–38 This work suggested the possibility of forming extra H-bonds by using both hydrogens of the N(2)–NH2 group. Inspired by this work, a carbonyl group was introduced at the meta-position of the 8-aryl position of 1b as illustrated in Fig. 3A. Based on this design, we hypothesized that N(2)–HA would form a H-bond with neighboring guanosine within the G-quartet, while the N(2)–HB could interact with the carbonyl group by forming a cross-layer H-bond.
Fig. 3 (A) General design of establishing a cross-layer H-bond; (B) single crystal structure and (C) H-bond information of [(1c)8Ba]2+·(Pic−)2.§ |
To confirm this idea, compound 1c was synthesized and applied to G-quadruplex construction upon interacting with metal cations. According to the 1H NMR spectra, treating 1c with Ba(Pic)2 led to the formation of a new G-quadruplex with one set of signals, similar to the G8-octamer obtained from 1b. Analysis of the NMR sample (in CDCl3) by ESI-MS gave a dominant, doubly charged peak with m/z at 2291.34, corresponding to [(1c)8Ba]2+ (mw = 4582.68).
Having confirmed the G8-octamer [(1b)8Ba]2+·(Pic−)2 formation, our next goal was to determine if there was a cross-layer H-bond as designed above. The 1H NMR spectra of [(1b)8Ba]2+·(Pic−)2 and [(1c)8Ba]2+·(Pic−)2 did not show peaks corresponding to N(2)H at room temperature. This is likely due to the rapid exchange between the two NH2 protons, even with the formation of a H-bond. Thus, the exchange rate of the two protons provides a direct indication of the H-bond strength in the G4-quartet. To explore the dynamic structure, variable temperature (VT) NMR experiments with [(1b)8Ba]2+·(Pic−)2 and [(1c)8Ba]2+·(Pic−)2 were performed and are summarized in Fig. 4.
Fig. 4 (A) VT NMR spectra of [(1b)8Ba]2+·(Pic−)2 and (B) VT NMR spectra of [(1c)8Ba]2+·(Pic−)2 confirmed the cross-layer H-bond design in establishing a cross-layer H-bond. |
For complex [(1c)8Ba]2+·(Pic−)2, the NH2 protons started appearing as broad peaks at 0 °C with the chemical shift at 10.28 ppm (HA) and 7.25 ppm (HB). In contrast, the VT NMR spectra of [(1b)8Ba]2+·(Pic−)2 did not show apparent peaks of the N(2)–H signals until further cooling the sample to −40 °C. The results indicated that there might be an extra H-bond in the 1c complex to lock the N(2)–NH2 from rapid exchange. Furthermore, a significantly downfield shifted chemical shift (7.25 ppm) was ascribed to the N(2)–HB proton in the (1c)8-octamer compared with the (1b)8-octamer (5.98 ppm). These observations provide clear evidence of the formation of a cross-layer H-bond in the (1c)8-octamer as designed.
Finally, the G8-octamer was verified by X-ray crystallography as shown in Fig. 3B. The crystal structure also confirmed the presence of the cross-layer interactions with a mean H-bond distance of 3.09 Å and a bond angle of 161.5°, suggesting a weak cross-layer H-bond present in solid state (Fig. 3C). This makes the structure a “self-assembled molecular-cuboid” purely constructed by H-bond linkage with all eight guanosine units. On the other hand, the distance of the two G-quartet layers (2.90 Å) in [(1c)8Ba]2+·(Pic−)2 remained similar to the 1b complex. Considering the specific “cage” size in the G8-octamer shown in Table 1, the cross-layer H-bond was not strong enough to generate extra enthalpy gain to balance the entropy cost caused by holding the two layers tighter.
With the confirmed cross-layer H-bonds, we evaluated the complex stability of [(1c)8K]+·(Pic−) in methanol. The results showed a similar stability to the complex formed with 1b (vide infra). Although the cross-layer H-bond approach could not boost G8-quadruplex stability in MeOH as anticipated, it provided an effective and novel approach to enhance G-quadruplex stability from monomer conformational design. Structural amendment was required to further improve the stability of the G8-octamer.
Fig. 5 (A) Synthesis of bis-guanosine derivatives. (B) Single crystal structure of [(1d′)4K]+·(Pic−).§ |
Self-assembly of 1d with K+ and Ba2+ cations in CDCl3 gave a similar one-set of signals in the 1H NMR spectra, consistent with the formation of a G8-octamer. ESI-MS of complexes from 1d and Ba2+ gave a doubly charged peak at m/z = 2175.55 as the dominant signal, indicating a mol. wt of 4351.1 for the supermolecule as [(1d)4Ba]2+. An attempt to obtain a single crystal of the 1d complex failed initially, resulting in a rather thin, film-like solid formation. Fortunately, the single crystal structure was successfully obtained by switching the monomer to 1d′ using DMSO as a co-solvent, confirming the cross-layer covalent linked structure as proposed (Fig. 5B). Notably, for complex [(1d′)4Ba]2+·(Pic−)2, the dihedral angle between 8-aryl and guanine is 40.7 degrees, similar to the dihedral angles in complex [(1b)8Ba]2+·(Pic−)2 (42.5 degrees). Overall, through G-monomer conformational analysis, a series of G8-octamers was successfully prepared with functionalization at 8-phenyl (1b), a cross-layer H-bond linker (1c) and a covalent linker (1d).
Fig. 6 (A) 1H NMR spectra showing G-quadruplex stability in CD3OD. (B) MS of [(1d)4Ba]2+·(Pic−)2 in MeOH. |
Fig. 7 DMSO-d6 titration: 1H NMR spectra of (A) [(1a)16K4]4+·(Pic−)4, (B) [(1b)8K]+·(Pic−), (C) [(1c)8K]+·(Pic−), (D) [(1d)4K]+·(Pic−) in CDCl3:DMSO-d6 with DMSO-d6 fractions of 10%, 20% and 50%. |
As shown in the 1H NMR spectra, G-quadruplexes from 1a, 1b and 1c started to dissociate in mixed solvents containing 20% DMSO-d6. Compared with the reported G16-hexadecamer from 1a,46 the G8-octamer formed by 1b and 1c showed comparable stability in 20% DMSO-d6. Eventually, all three G-quadruplexes gave complete dissociation in 50% DMSO-d6 solution with only one set of signals corresponding to the monomer. In contrast, [(1d)4K]+·(Pic−) showed significantly improved stability, with only 2% complex dissociation in 50% DMSO-d6. This result demonstrates the significantly enhanced stability of G-quadruplexes constructed by 1d.
To quantify the thermodynamic stability of G-quadruplexes,47 VT NMR experiments were carried out for pure complexes of [(1a)16K4]4+·(Pic−)4, [(1b)8K]+·(Pic−), [(1c)8K]+·(Pic−), and [(1d)4K]+·(Pic−) in CDCl3/DMSO-d6 with a fraction of 20% DMSO-d6. The values of complex dissociation enthalpy and entropy for each G-quadruplex were calculated from van't Hoff plots and are compared in Table S6 (see detailed discussion in the ESI†). For complex [(1d)4K]+·(Pic−), no significant increase in monomer concentration was observed with the increase in temperature. This might be attributed to the high kinetic barrier for G-quadruplex [(1d)4K]+·(Pic−) dissociation. To confirm this hypothesis, a NOESY experiment at 50 °C was performed for [(1b)8K]+·(Pic−), [(1c)8K]+·(Pic−) and [(1d)4K]+·(Pic−) (Fig. S14–S16†). The results suggested that kinetic exchange between the complex and monomer for [(1d)4K]+·(Pic−) was too slow to be recorded by NMR spectroscopy.
In addition to DMSO titration, the stability of G8 and G16 complexes could also be evaluated using tandem-MS by increasing the operating voltage. The cation fragments of G-quadruplexes were separated and treated with increasing voltage. The operating voltage for G-quadruplex cation fragment decomposition are summarized in Table 2.
It is noteworthy that [(1a)16Ba2]4+·(Pic−)4 only showed a doubly charged peak at m/z = 2123.31 corresponding to [(1a)8Ba]2+, indicating that the picrate bridge dissociated under the MS conditions. Further comparison of all the cation fragments of the G-quadruplexes revealed a clear stability trend as [(1d)4Ba]2+ > [(1c)8Ba]2+ > [(1b)8Ba]2+ > [(1a)8Ba]2+. Overall, the covalent linking strategy significantly helped to stabilize the G8-octamer, both in H-bond competitive solvents and gas phase.
Footnotes |
† Electronic supplementary information (ESI) available: Experimental section, NMR spectra, ESI-MS spectra and crystallographic data. CCDC 1871565–1871570 and 1871754. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c9sc00190e |
‡ Equal contribution. |
§ All the structures reported in this article have been deposited with the Cambridge Crystallographic Data Centre. The accession numbers for [(1b)8Ba]2+·(Pic−)2; [(1b)8K]+·(PF6−); [(1b)8K]+·(Pic−); [(1b)8Sr]2+·(Pic−)2; [(1b)8Rb]+·(Pic−); [(1c)8Ba]2+·(Pic−)2; [(1d′)4K]+·(Pic−) reported in this paper are: 1871565, 1871566, 1871567, 1871568, 1871569, 1871570, 1871754, correspondingly. |
This journal is © The Royal Society of Chemistry 2019 |