J. Mꢁller, F. M. Bickelhaupt et al.
biguously a 0.06 ratio from a 0.11 ratio, a different approach has been
used to obtain interpretable Job plots. Instead of defining one duplex as
one of the components, we have defined one nucleobase as this compo-
nent. Hence, a 1:1 stoichiometry in the Job plot (one Ag+ ion per nucleo-
base) corresponds to two Ag+ ions per base pair. Accordingly, the incor-
poration of only one Ag+ ion per base pair would be represented by a
2:1 stoichiometry in the Job plot (two nucleobases per Ag+ correspond-
ing to one Ag+ per base pair). For example, for the Job plot shown in
Figure 3d, the total concentration was 150 mm. Hence, at cnucleobase =1, the
nucleobase concentration was 150 mm (corresponding to 8.3 mm oligonu-
corrected DFT calculations have corroborated these experi-
mental findings, revealing the favored formation of a two-
fold-metalated base pair bearing one single positive charge.
Bond analyses have indicated a substantial “argentophilic”
Ag-d10–Ag-d10 interaction of about 16 kcalmolꢀ1 caused by
charge transfer from silver 4d to 4s orbitals, which contrib-
utes to the stability of doubly metal-mediated base pairs
such as bp3.
The observations presented herein allow important con-
clusions to be drawn regarding the generation and use of
metal-mediated base pairs in the context of nucleic acid
functionalization. 1) Subtle changes of a nucleobase in-
volved in metal-mediated base pairing may have far-reach-
ing consequences with respect to the metal-binding proper-
ties. 2) The oligonucleotide sequence can also play a signifi-
cant role in determining the stability of the metal-mediated
base pair. 3) The possibility of incorporating two metal ions
into one base pair represents a significant step forward to-
wards obtaining highly functionalized nucleic acids.
4) Mechanistic insights gained from calculations suggest that
the metal ions are inserted in one of the late steps of base-
pair formation. Interestingly, even though Hoogsteen du-
plexes comprising nonmetalated 1,3-dideazaadenine–thy-
mine base pairs are not stable, this base pair might still be
involved as an intermediate during the formation of its
Ag+-mediated analogue.
cleotide single strand or 4.2 mm double helix). Accordingly, at cnucleobase
=
0.1, the nucleobase concentration was 15 mm (corresponding to 0.83 mm
oligonucleotide single strand or 0.42 mm double helix), and the Ag+ con-
centration was 135 mm.
DLS experiments were performed on a DynaPro Titan instrument from
Wyatt Technology. The oligonucleotide samples (c=0.3 mm, 150 mm
NaClO4, 5 mm MOPS, pH 6.8) were incubated in the presence of the re-
spective amounts of AgNO3 by heating to 908C and then slowly cooled
to 108C (cooling rate: 0.5 Kminꢀ1). The samples were then centrifuged at
13500 rpm for at least 2 h. The supernatant was transferred to a 12 mL
cuvette and measured.
LILBID experiments were performed on an in-house-constructed time-
of-flight mass spectrometer (TOF-MS) of the Wiley–McLaren type. De-
tails of this method have been published previously.[33] Briefly, the ion
source contains a commercial droplet dispenser (Microdrop), which in-
jects on demand tiny micro droplets (Øꢂ50 mm; Vꢂ65 pL) from
300 Torr through
a pressure-reduction aperture into high vacuum
(10ꢀ6 Torr). There, they are individually irradiated by high-intensity mid-
IR laser pulses (lꢂ3 mm) generated by an in-house-constructed
Nd:YAG-pumped LiNbO3 optical parametric oscillator. Due to the ab-
sorption of the laser light by the solvent, beyond a certain intensity
threshold the droplets are disrupted, ejecting preformed analyte ions into
the vacuum. For the LILBID MS experiments, oligonucleotide samples
(20 mm oligonucleotide, 0 or 1 mm AgNO3 in 100 mm NH4HCO3/1 mm
Mg
ACHTUNGTRENNUNG
Experimental Section
buffer exchange to 20 mm NH4HCO3/1 mm MgACHTUNGTRENNNUG
All calculations were performed using dispersion-corrected density func-
tional theory (DFT-D) as implemented in the Amsterdam Density Func-
tional program (ADF 2009.1).[34] Geometry optimizations were per-
formed using analytical gradient techniques and without any symmetry
or geometrical constraints. All geometries were verified as true energy
minima, having no imaginary frequencies, by means of vibrational analy-
ses. The dispersion-corrected GGA functional BLYP-D,[35] which has
proved to be adequate for calculations of DNA base dimers,[36] was used
with the doubly polarized, triple-z-quality basis set TZ2P.[37]
In BLYP-D, the BLYP functional[35c,d] is augmented with an empirical cor-
rection for long-range dispersion effects, described by a sum of damped
interatomic potentials of the form C6Rꢀ6 added to the usual DFT ener-
gy.[35a,d] The TZ2P basis set used herein is a large uncontracted set of
Slater-type orbitals (STOs) containing diffuse functions (no Gaussian
functions are involved).[34] The basis set is of triple-z quality for all atoms
and is augmented with two sets of polarization functions, that is, 3d and
4f on C, N, O and 2p, 3d on H. The 1s core shells of C, N, and O were
treated by the frozen-core approximation. An auxiliary set of s, p, d, f,
and g STOs was used to fit the molecular density and to accurately repre-
sent the Coulomb and exchange potentials in each self-consistent field
cycle. The basis set superposition error (BSSE) in the bond energy has
previously been shown to be only 1 kcalmolꢀ1 or less in GGA calcula-
tions involving the TZ2P basis set. Moreover, the dispersion correction in
the present calculations has been developed such that these small BSSE
effects are absorbed into the empirical potential.[35a] Therefore, no explic-
it counterpoise corrections had to be carried out. In addition, scalar rela-
tivistic effects were taken into account by using the zeroth-order regula-
tor approximation (ZORA).[38]
1,3-Dideaza-2’-deoxyadenosine 5 and its FMOC-, DMT-, and CEDIP-
protected analogue 8 were synthesized according to (partially modified)
previously published procedures.[5b,17–19] DNA syntheses were performed
in the DMT-off mode on a K&A Laborgerꢂte H8 DNA/RNA synthesizer
following standard protocols (except for a threefold increase in coupling
time for the artificial phosphoramidite). The oligonucleotides were
cleaved from the solid support and deprotected by exposing them to tert-
butylamine/methanol/water (1:2:1) for 3 h at 658C. They were purified by
polyacrylamide gel electrophoresis under denaturating conditions (gel so-
lution: 7m urea, 1m TBE buffer; 8–18% polyacrylamide/bisacrylamide
(29:1); loading buffer: 11.8m urea; 42 mm Tris·HCl, pH 7.5; 0.83 mm
EDTA, pH 8.0; 8% sucrose). After purification, the oligonucleotides
were desalted by passage through NAP 10 columns. The desalted oligo-
nucleotides were identified by MALDI-TOF mass spectrometry (I: calcd
for [M+H]+: 5475 Da; found: 5476 Da; II: calcd for [M+H]+: 5475 Da;
found: 5474 Da). MALDI-TOF mass spectra were recorded on a Bruker
Reflex IV instrument using a 3-hydroxypicolinic acid/ammonium citrate
matrix. During quantification of the oligonucleotides, a molar extinction
coefficient e260 =7.5 cm2 mmolꢀ1 was used for 5.
UV/Vis spectra were recorded on a Varian CARY BIO 100 spectropho-
tometer. CD spectra were recorded at 108C on a Jasco J-815 spectrome-
ter. UV melting curves were recorded at 260 nm with a heating rate of
1 Kminꢀ1 and a data interval of 0.2 K. Melting temperatures were deter-
mined as the maxima of the first derivatives of the melting curves. Prior
to each measurement, the samples were equilibrated by heating to 908C
and then slowly cooling to 108C (cooling rate: 0.5 Kminꢀ1).
For the generation of a Job plot, the sum of the molar fractions of both
components needs to be constant during all measurements. When using a
double helix (with 18 base pairs) as one component and an Ag+ ion as
the other component, a 1:18 stoichiometry (one Ag+ ion per base pair)
would have to be distinguished from a 1:36 stoichiometry (two Ag+ ions
per base pair). As it would be almost impossible to differentiate unam-
For calculations in water, solvent effects were estimated using the con-
ductor-like screening model (COSMO).[39] In the COSMO calculations, a
corrected ion radius of Ag+ and a corrected solvation energy of H+ were
used. Parameterization was accomplished by performing
a series of
6542
ꢃ 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim
Chem. Eur. J. 2011, 17, 6533 – 6544