Aukkrapon Dangsopona,
Nattawee Poomsukb,
Khatcharin Siriwongb,
Tirayut Vilaivanc and
Chaturong Suparpprom*a
aDepartment of Chemistry and Center of Excellence for Innovation in Chemistry, Faculty of Science, Naresuan University, Muang, Phitsanulok 65000, Thailand. E-mail: chaturongs@nu.ac.th
bMaterials Chemistry Research Center, Department of Chemistry and Center of Excellence for Innovation in Chemistry, Faculty of Science, Khon Kaen University, Khon Kaen 40002, Thailand
cOrganic Synthesis Research Unit, Department of Chemistry, Faculty of Science, Chulalongkorn University, Phayathai Road, Patumwan, Bangkok 10330, Thailand
First published on 27th July 2016
This work aims to explore a new clickable carbazole-based fluorescent label, its incorporation into pyrrolidinyl peptide nucleic acid (acpcPNA), and the interactions between the labeled PNA probe and DNA by fluorescence spectrophotometry and molecular dynamics simulation. A carbazole derivative, namely 3,6-diaminocarbazole (DAC), was synthesized and incorporated into the internal and terminal positions of azide-modified acpcPNA via a sequential reductive alkylation and click reaction previously developed by our group. The DAC-modified acpcPNA can form a stable hybrid with complementary DNA with somewhat lower stability compared to unmodified acpcPNA. Most importantly, the DAC-modified acpcPNA exhibits a remarkable fluorescence increase in the presence of DNA (up to 35.5 fold with complementary DNA). Non-complementary as well as single mismatched DNA targets gave a smaller fluorescence increase (1.1 to 18.6 fold), and the discrimination could be further improved by increasing the temperature to dissociate the mismatched hybrids. Molecular dynamic simulations revealed that the DAC interacts with adjacent nucleobases in single stranded PNA, resulting in quenching of the fluorescence signal. When the PNA formed a hybrid with DNA, the DAC was pushed away from the duplex, resulting in a fluorescence increase. Thus, the DAC-labeled acpcPNA is a potential candidate as a self-reporting probe for determination of DNA sequences.
Carbazole is an electron rich aromatic molecule that showed many interesting properties such as high thermal stability, excellent photophysical and electroconductive properties and the ease of synthesis and chemical modification at 3-, 6- and 9-positions. Carbazole derivatives, both in monomeric and oligomeric forms, have found widespread applications as photoconductivity polymers,7–10 polymeric light emitting devices,11,12 high temperature thermostable polymers,13 electrochromic materials,14–16 dye-sensitized solar cell and dendrimeric electro phosphorescent devices.17–19 In addition, carbazole derivatives have been used as fluorophores in sensor applications, including DNA detection.20 Carbazole derivatives with cationic groups at 3,6-positions have been reported to bind in the DNA groove of AT-rich sequences.21–24 However, the use of carbazole derivatives in combination with an oligonucleotide or analogues for sequence-specific detection of DNA sequences is unknown.
Peptide nucleic acid or PNA is a DNA mimic consisting of nucleobases attached onto a peptide-like backbone.25,26 PNA can hybridize with DNA according to Watson–Crick base pairing rules with high affinity and sequence specificity. The absence of phosphate group in the PNA backbone results in decreasing of the electrostatic repulsion with DNA, and therefore PNA can form a more stable duplex with DNA than DNA itself. In addition, the binding of PNA and DNA is highly sequence specific and is not sensitive to ionic strength. These properties, together with the excellent chemical and biological stabilities, make PNA a promising tool for several applications in medical and biotechnology fields. Recently, we reported a series of new PNA system with restricted conformation flexibility as a result of incorporation of cyclic structures into the PNA backbone.27–30 Pyrrolidinyl PNA probes labeled with environment sensitive dyes such as pyrene,31,32 Nile red,32,33 or other DNA binding dyes such as thiazole orange34 have been successfully applied as self-reporting fluorescence probes for DNA sequence detection.
From the remarkable photophysical properties and DNA binding ability mentioned above, carbazole has a potential to be developed into a probe for DNA detection by fluorescence spectroscopy. In the present work, an electron rich carbazole derivative, namely 3,6-diaminocarbazole (DAC), was incorporated at an internal position of acpcPNA in order to develop a new self-reporting fluorescence hybridization probe. Here we report the synthesis, thermal stabilities and fluorescence properties of DAC-modified acpcPNA probes, as well as their ability to discriminate between complementary and single-mismatched DNA targets. Moreover, molecular dynamics simulation was performed to explain the behavior of DAC-labeled acpcPNA probes.
AcpcPNA sequences were synthesized manually by Fmoc solid phase peptide synthesis as described previously.27,35 In order to provide a handle for attachment of the DAC via the previously reported reductive alkylation-Click chemistry strategy,34 one of the 2-amino cyclopentanecarboxylic acid (ACPC) spacer in the middle of the acpcPNA strand was replaced with 3-aminopyrrolidine-4-carboxylic acid (APC) for internal labeling.35 For terminal labeling, no modification with APC was required, but the last ACPC residue was omitted from the PNA sequence. After end-capping by acetylation (this step was skipped for the terminal modification) and removal of all protecting groups from the nucleobase and the APC spacer, the azide group was introduced by reductive alkylation with 4-azidobutanal. Finally, the clickable DAC label (5) was attached to the azide-modified acpcPNA on the solid support by Cu-catalyzed click chemistry. Four DAC-modified acpcPNA sequences were designed to investigate the DNA binding and fluorescent properties under different sequence context. These include an internally-labeled homothymine (T9DAC), two internally-labeled mix sequences with different flanking bases (A/T and C/G pairs) (M10ATDAC and M10GCDAC) and a terminally-labeled mix sequence (DAC-M10). All modified PNAs were purified by reverse phase HPLC (to >90% purity). The identities of the obtained PNA sequences were confirmed by MALDI-TOF mass spectrometry, which showed the mass signals that are in a good agreement with the expected values. The sequences and characterization data of all PNA involved in this studies are shown in Table 1.
| PNA | Sequence (N → C) | tRa (min) | Yieldb (%) | m/z (calcd) | m/z (found) |
|---|---|---|---|---|---|
a HPLC condition: C18 column 4.6 × 50 mm, particle size 3 μm, gradient 0.1% TFA in H2O : MeOH 90 : 10 for 5 min then linear gradient to 10 : 90 over 30 min, flow rate 0.5 mL min−1.b Isolated yield after HPLC. |
|||||
| T9DAC | Ac-TTTT(DAC)TTTTT-LysNH2 | 32.0 | 5.8 | 3513.6 | 3512.4 |
| M10ATDAC | Ac-GTAGA(DAC)TCACT-Lys NH2 | 29.5 | 5.2 | 3892.8 | 3894.1 |
| M10GCDAC | Ac-GTAGC(DAC)GCACT-Lys NH2 | 27.7 | 4.4 | 3893.8 | 3893.2 |
| DAC-M10 | DAC-GTAGATCACT-Lys NH2 | 29.0 | 5.4 | 3738.8 | 3738.3 |
| Entry | PNA | DNA (3′ → 5′) | Tm (°C) | F/F0 | ΦF | ΦF (ds)/ΦF (ss)b | Notesc |
|---|---|---|---|---|---|---|---|
| a All measurements were performed in 10 mM sodium phosphate buffer (pH 7.0), [PNA] = 1.0 μM, [DNA] = 1.2 μM, λex 315 nm, λem 450 nm.b Quantum yields were measured in phosphate buffer using quinine sulfate (0.540) as a reference.c ss = single stranded PNA, comp = complementary, sm = single mismatched. | |||||||
| 1 | T9DAC | None | — | — | 0.009 | — | ss |
| 2 | AAAAAAAAA | 65.2 | 35.5 | 0.242 | 26.7 | comp | |
| 3 | AAAACAAAA | 33.7 | 18.6 | 0.139 | 15.3 | smC | |
| 4 | AAAAGAAAA | 23.6 | 6.0 | 0.041 | 4.5 | smG | |
| 5 | AAAATAAAA | 38.2 | 15.7 | 0.160 | 17.6 | smT | |
| 6 | M10ATDAC | None | — | — | 0.007 | — | ss |
| 7 | AGTGATCTAC | 53.4 | 13.2 | 0.083 | 11.9 | comp | |
| 8 | AGTGCTCTAC | 29.0 | 2.7 | 0.012 | 1.7 | smC | |
| 9 | AGTGGTCTAC | 24.5 | 1.7 | 0.014 | 2.0 | smG | |
| 10 | AGTGTTCTAC | 35.0 | 6.5 | 0.034 | 4.8 | smT | |
| 11 | M10GCDAC | None | — | — | 0.011 | — | ss |
| 12 | AGTGCGCTAC | 59.0 | 10.2 | 0.117 | 11.1 | comp | |
| 13 | AGTGATCTAC | 30.8 | 4.9 | 0.099 | 9.4 | smA | |
| 14 | AGTGGTCTAC | 33.4 | 7.0 | 0.085 | 8.0 | smG | |
| 15 | AGTGTTCTAC | 39.2 | 9.1 | 0.104 | 9.9 | smT | |
Next, fluorescence experiments were performed with all PNA sequences and their DNA hybrids. The responsiveness of the PNA probe was reported in terms of fluorescence intensity ratio (F/F0; where F = fluorescence intensity of PNA·DNA hybrid, F0 = fluorescence intensity of single stranded PNA; λex 315 nm; λem 450 nm). While free 3,6-diaminocarbazole showed a high ΦF of 0.489, all internally-labeled single stranded PNA probes exhibited very low ΦF values of 0.01 or less, regardless of the identity of the neighboring nucleobase (Table 2, entries 1, 6, 11). This suggests that the DAC label is quenched upon incorporation into the PNA, probably by interactions with the nucleobases similar to pyrene.31 The sequences M10ATDAC and T9DAC showed lower quantum yields than M10GCDAC, which means that A/T are better quenchers than C/G bases for DAC label. Nevertheless, judging from the uniformly low quantum yields of all PNA probes, the difference in quenching efficiency by different nucleobases is much less than pyrene, which emphasizes the important feature of the DAC label.31
A remarkable increase in fluorescence signal was observed upon hybridization with complementary DNA, as shown by the fluorescence intensity ratio (F/F0) between 10.2 and 35.5 (Table 2, entries 2, 7, 12). The quantum yield of the DAC-labeled PNA·DNA duplexes was also consistently increased by >10-fold to around 0.08–0.24, which is still considerably lower than the free DAC (3). This indicates that the interaction between the DAC and the neighboring nucleobases is reduced in PNA·DNA duplexes due to the formation of base pairs. This behavior is similar to pyrene label and a similar mechanism of fluorescence increase is proposed.31 It should be noted that free DAC also interact with single stranded DNA, resulting in fluorescence quenching, which is consistent with the proposed mechanism (see ESI Fig. S22†). CD spectra of M10GCDAC and its complementary DNA hybrid (see ESI Fig. S23†) also clearly confirms the formation of PNA·DNA hybrids as shown by the change in CD signals at 210 and 260 nm. No CD signal was observed in the UV absorption region of DAC, suggesting that the DAC label was not well-oriented and thus did not appreciably interact with the PNA·DNA duplex.
Single mismatched PNA·DNA hybrids generally showed smaller fluorescence increase, but some variations in the sequence context were observed (Fig. 1). The M10ATDAC showed a better mismatch discrimination than the T9DAC and M10GCDAC probes. Possible explanation is that certain mismatched PNA·DNA hybrids (Table 2, entries 3, 5, 10, 13–15) are still sufficiently stable at 20 °C – the temperature at which the fluorescent experiments were carried out. Completely unrelated DNA target provided a signal that was indistinguishable from the single stranded probes in all cases. For terminal labeled PNA probe DAC-M10, no fluorescence change could be observed for both complementary and mismatched DNA targets (Fig. 1b). Since the single stranded PNA probe exhibited low fluorescence similar to internally labeled PNA probes, it can be concluded that the quenching by neighboring nucleobases still operates. Tm measurement confirms the strong binding between the DAC-M10 probe and DNA in a sequence specific fashion (58 °C for complementary duplexes and <20 °C for mismatched duplexes). The absence of fluorescence increase upon duplex formation was therefore explained by the interaction between the DAC and the terminal GC base pair of the PNA·DNA duplex similar to the pyrene label.31
These results clearly demonstrate that the fluorescence of internally DAC-labeled acpcPNA is low in single stranded PNA probes due to quenching by neighboring nucleobases, especially T and A bases. This interaction diminished – resulting in a large fluorescence increase – when the PNA·DNA duplexes are formed, provided that the DAC label does not locate at the end of the duplex which allows interactions with the terminal base pairs.
Since the maximal difference will be obtained at the temperature whereby the mismatched duplexes were dissociated as much as possible and the complementary duplexes are still stable, it is important to vary the temperature to find the optimal value case by case. The specific temperature that provided maximum differentiation in the case of T9DAC, M10ATDAC and M10GCDAC was 55, 45 and 50 °C, respectively. The mismatched hybrids of T9DAC and M10GCDAC were more stable than those of M10ATDAC and therefore required a higher temperature to complete the separation. At the aforementioned optimal temperature, the fluorescence ratio between complementary and single mismatch hybrids improved from 1.0–2.5 folds to 6.0–21.1 folds (Fig. 2b and c and Table 3). The results confirm that the specificity of DAC-labeled acpcPNA in distinguishing between complementary and single-mismatched DNA by DAC-labeled PNA probes can be improved by simply increasing the temperature without requiring more complex operations such as enzymatic digestion.
| Entry | PNA | DNA (3′ → 5′) | F/F0 (25 °C) | Fcomp/Fmm (25 °C) | F/F0 (high temp)b | Fcomp/Fmm (high temp)b |
|---|---|---|---|---|---|---|
| a All measurements were performed in 10 mM sodium phosphate buffer (pH 7.0), [PNA] = 1.0 μM, [DNA] = 1.2 μM, λex 315 nm, λem 450 nm.b High temperature: T9DAC = 55 °C; M10ATDAC = 45 °C; M10GCDAC = 50 °C. | ||||||
| 1 | T9DAC | AAAAAAAAA | 35.4 | — | 35.7 | — |
| 2 | AAAACAAAA | 18.6 | 1.9 | 2.5 | 14.2 | |
| 3 | AAAAGAAAA | 6.0 | 5.9 | — | — | |
| 4 | AAAATAAAA | 15.7 | 2.3 | 1.7 | 21.1 | |
| 5 | TTTTTTTTT | 1.1 | 32.4 | — | — | |
| 6 | M10ATDAC | AGTGATCTAC | 13.2 | — | 12.9 | — |
| 7 | AGTGCTCTAC | 2.6 | 5.2 | 1.3 | 10.0 | |
| 8 | AGTGGTCTAC | 1.6 | 8.2 | — | — | |
| 9 | AGTGTTCTAC | 5.7 | 2.3 | 1.2 | 10.6 | |
| 10 | GGGGGGGGGG | 1.2 | 10.9 | — | — | |
| 11 | M10GCDAC | AGTGCGCTAC | 10.1 | — | 10.2 | — |
| 12 | AGTGATCTAC | 5.3 | 1.9 | — | — | |
| 13 | AGTGGTCTAC | 7.1 | 1.4 | 1.0 | 10.3 | |
| 14 | AGTGTTCTAC | 9.0 | 1.1 | 1.7 | 6.0 | |
| 15 | GGGGGGGGGG | 0.9 | 11.1 | — | — | |
Next, to support the proposed mechanism of fluorescence change of the DAC-labeled PNA probe, molecular dynamics simulation and energy calculation were performed on both single stranded and duplex PNA probe.
| System | ΔEbinding (kcal mol−1) |
|---|---|
| a The calculation is expressed as ΔEbinding = EDAC–PNA − EDAC − EPNA, where EDAC–PNA, EDAC and EPNA are the total energy of DAC and three PNA nucleotides, DAC moiety and three PNA nucleotides, respectively. | |
| ssPNA | −11.40 |
| PNA·DNA | −3.86 |
C), 1720, 1702 (C
O), 1524, 1490, 1048, 1023, 864, 792, 578 cm−1; HRMS (ESI+): m/z calcd for C25H29N3O4Na+: 458.2056 [M + Na]+ found: 458.2092, m/z calcd for C25H29N3O4K+: [M + K]+ 474.1795 found: 474.1826.
:
1 aqueous ammonia/dioxane at 60 °C overnight. The apc-modified acpcPNA (0.5 μmol) was alkylated at the pyrrolidine ring nitrogen of the APC spacer with 4-azidobutanal (15 μmol, 30 equiv.) in the presence of sodium cyanoborohydride (30 μmol, 60 equiv.) and acetic acid (30 μmol, 60 equiv.) in methanol (100 μL) at room temperature overnight.34 The azide-modified acpcPNA was further reacted with N9-propargyl-3,6-bis(Boc-amino)carbazole (5) (7.5 μmol, 15 equiv.) while still attached on the solid support in the presence of tris[(1-benzyl-1H-1,2,3-triazol-4yl)methyl]amine (TBTA, 30 μmol, 60 equiv.), tetrakis(acetonitrile)copper(I) hexafluorophosphate (15 μmol, 30 equiv.) and (+)-sodium-L-ascorbate (60 μmol, 120 equiv.) in 3
:
1 (v/v) DMSO
:
tBuOH at room temperature overnight.34 After the reaction was completed as monitored by MALDI-TOF analyses (12 h), the diaminocarbazole-modified acpcPNA was cleaved from the resin with 50% trifluoroacetic acid (TFA) in dichloromethane and 10% triisopropylsilane (TIS) (500 μL × 30 min × 3). The combined cleavage solution was dried under a stream of nitrogen gas and the crude PNA was precipitated by addition of diethyl ether (Scheme 2). The crude PNA was purified by reverse-phase HPLC with using mobile phase A
:
B (A: 0.1% TFA in H2O; B: 0.1% TFA in MeOH) 9
:
1 for 5 min followed by a linear gradient to A
:
B 1
:
9 over 55 min. The purified PNA was characterized by MALDI-TOF mass spectrometry (Microflex, Bruker Daltonics) by using α-cyano-4-hydroxy-cinnamic acid (CCA) as a matrix.
| ΦF(sample) = ΦF(standard)[m(sample)/m(standard)][η(sample)2/η(standard)2] | (1) |
![[A with combining low line]](https://www.rsc.org/images/entities/char_0041_0332.gif)
CACT whereby the DAC was linked to the apc spacer at the underlined position, and the corresponding acpcPNA·DNA duplex. Because the three-dimensional structures of ssPNA and PNA·DNA duplex were not available, the starting structure of PNA·DNA duplex was built up based on the corresponding sequence of DNA·DNA duplex (dGTAGATCACT·dAGTGATCTAC) which was generated using the NAB module of the AMBER 12 package.37 The first strand of the duplex was replaced with PNA strand in an antiparallel direction yielding the PNA·DNA duplex.38 The structure of diaminocarbazole (DAC) linked to backbone of PNA was first optimized with B3LYP/6-31G* calculation using Gaussian 09 (ref. 39) before adding to the PNA strand. The initial structure of single stranded acpcPNA was obtained in a similar manner. The force field parameters of the PNA and DAC were generated using the standard procedure.38 To setup the simulations, the starting structure of each system was immersed in a rectangular water box with TIP3P water molecules40 extended by 10 Å in each direction from the solute and then neutralized by adding Na+ counterions (for PNA·DNA system). The MD simulations were carried out using AMBER 12 suite of program with parm10 force field complemented with the prepared parameters of PNA and DAC, SHAKE algorithm41 employed for all hydrogen atoms, a 9 Å cutoff applied to the non-bonding Lennard-Jones interactions, the particle mesh Ewald method used to account for long-range interactions,42 time step of 1 fs and NPT ensemble (P = 1 atm and T = 300 K). The simulation protocol consists of heating, equilibration and production steps, and was described in details elsewhere.38 The MD trajectories were generated for 10 ns and the coordinates were stored every 1 ps, totally 10
000 structures.Footnote |
| † Electronic supplementary information (ESI) available. See DOI: 10.1039/c6ra15861g |
| This journal is © The Royal Society of Chemistry 2016 |