Skip to main content

Dinuclear complex-induced DNA melting


Dinuclear copper complexes have been designed for molecular recognition in order to selectively bind to two neighboring phosphate moieties in the backbone of double strand DNA. Associated biophysical, biochemical and cytotoxic effects on DNA were investigated in previous works, where atomic force microscopy (AFM) in ambient conditions turned out to be a particular valuable asset, since the complexes influence the macromechanical properties and configurations of the strands. To investigate and scrutinize these effects in more depth from a structural point of view, cutting-edge preparation methods and scanning force microscopy under ultra-high vacuum (UHV) conditions were employed to yield submolecular resolution images. DNA strand mechanics and interactions could be resolved on the single base pair level, including the amplified formation of melting bubbles. Even the interaction of singular complex molecules could be observed. To better assess the results, the appearance of treated DNA is also compared to the behavior of untreated DNA in UHV on different substrates. Finally, we present data from a statistical simulation reasoning about the nanomechanics of strand dissociation. This sort of quantitative experimental insights paralleled by statistical simulations impressively shade light on the rationale for strand dissociations of this novel DNA interaction process, that is an important nanomechanistic key and novel approach for the development of new chemotherapeutic agents.


The design of new cytostatic compounds with novel mechanisms of action for treatment of cancer has been an ongoing task, in which maximum apoptosis effectiveness for various tumor types is desired, as well as to limit the damage to the surrounding tissue [1,2,3]. Therefore, a new class of molecules had been synthesized, imitating the hydrolytic cleaving mechanism of metalloenzymes to the backbone of DNA, but regulated by the steric structure of the complex and even tunable by the choice of metal cores and terminal donors [4,5,6,7]. This way, cytotoxic effects by blocking DNA polymerase could be anticipated [8].

The complexes have been synthesized by coordination chemistry with metal ions such as NiII [5, 8], albeit in the following investigations, complexes with CuII were chosen, since their impact on DNA was shown to be less radical and more feasible to observe [4]. Metal ions display differing hydrolytic cleaving reactivities and can therefore be used to tune the effectiveness to the targeted cancer type [9,10,11].

The interaction of these complexes with dsDNA has been investigated in depth employing various methods, although exact details of the binding modalities could not finally be uncovered. Therefore, this report aims at employing AFM under UHV conditions to ultimately shed light on the binding scheme, with its unmatched high-resolution capabilities.

In the previous work of Jany et al. [4], investigations using techniques including polymerase chain reaction (PCR), gel electrophoresis, cell arrays and even optical and magnetic tweezers showed that the complexes indeed bind successfully to DNA and cause cell toxicity. The binding was shown to be irreversible and to cause the strands to form intra- and interstrand aggregation.

Since the detailed structural peculiarities of DNA strands are far below the resolution limit of conventional microscopy, AFM is the method of choice for the structural elucidation without further treatments. To get access to details below the grooves, one needs to get rid of adverse influences of ambient conditions, i.e., moisture adsorbate layers, which is only provided by environments of liquids or vacuum. Indeed, measurements under liquid conditions showed remarkable detailed images [12, 13]. There have also been a few investigations of DNA under vacuum [14], though these were mostly limited to the use of substrates that provide strong adhesion of biomolecules, like mica [15,16,17], since the focus mostly lays on the properties of the strands under physiological conditions before adsorption. Also, the sample preparations relied on the coating and drying of solutions still under ambient conditions, since DNA is not suited for thermal evaporation.

In surface science, a relatively new method of UHV in-situ sample preparation employs the principles of electrospray ionization, originating from mass spectrometry. Hereby, a solution of the desired molecules gets inserted into vacuum where droplet Coulomb fission is provoked by high electric fields [18]. Ideally, the entire solvent evaporates before the molecules reach the sample, where they adhere gently to the clean surface, without contaminants of residual solvent [19,20,21,22]. With this new and innovative sample deposition method available, DNA can finally be investigated with extraordinary high resolution [23].


Figure 1a, b illustrate the design of the [(HtomMe){Cu(OAc)}2]+ complexes, with two metal cores that are positioned by the naphthalene backbone at a distance of 0.6–0.7 nm, providing the precise molecular recognition ability of two neighboring phosphate moieties of the DNA backbone. Each of the two pyridyl groups coordinate around the metal atoms in such a way that they sterically limit the hydrolytic cleavage reaction and hamper the nicking of the strand [24].

Fig. 1
figure 1

a Chemical structure of the [(HtomMe){Cu(OAc)}2]+ complexes. The descriptor for the complexes changes according to the chemical surroundings, for clarification see Additional file 1: Figure S1. b Model of binding modes of the complex molecules to adsorb to DNA and to each other. c UHV-AFM topographic image of the agglomeration of multiple strands (frequency shift setpoint df =− 2.5 Hz, oscillation amplitude A = 14.1 nm, resonance frequency f0 = 259.0 Hz) and the dedicated interpretation d: two double strands (yellow, height 0.6 nm), one singular ds (orange, height 0.4 nm), melting bubble dissociation in two single strands (red, height 0.1 nm)

A solution of λ-phage DNA solution (Lambda DNA Hind III Digest, Sigma-Aldrich, 12.5 ng/µl) (λ-DNA) was treated with a 0.4 nM solution of the dinuclear copper complexes [(HtomMe){Cu(OH2)}2]3+, with an incubation time of 45 h, as this was determined to display an indicated onset of expected effects to the DNA. Regarding the concentration dependent influence of the complexes to DNA, a dilution series can be found in the Additional File 1: (Figure S2).

Backside bonding

A result that was already expected from investigations under ambient conditions was the interstrand agglomeration [4]. The high-resolution images taken under UHV conditions in Fig. 1c, d show how multiple strands align in order to agglomerate, which can be geometrically traced to individual strands.

As Jany et al. already speculated, this effect results from π-stacking effects of the outwards-oriented naphthalene und pyridyl residue moieties of [(HtomMe){Cu(OAc)}2]+. X-ray diffraction showed that the freely rotatable pyridyl groups orient in such a way that, together with the center naphthalene, they form three parallel ridges, which would allow two molecules to attach by coercing multiple π-stacking interactions.

Melting bubbles

Another remarkable feature is introduced in Fig. 2. Along the strand appear several partial fissions of various length into thinner strands, with these sectors ranging in length from 1 to 8 nm. They are not only characterized by a reduced height (0.1–0.2 nm in contrast to 0.4 nm of the main strand), but also by a string of bead-like appearance with equidistant protrusions along a strand.

Fig. 2
figure 2

a Dissociated melting bubbles along a double strand of λ-DNA (taken under UHV, df =− 2.0 Hz, A = 14.1 nm, f0 = 259.0 Hz). b Closer look at one of the larger bubbles, including an area of supercoiled buildup. c Distribution of interval distances between the protrusion features along the bubbles, with a mean of 1.42 nm. d Distances and antiphasic patterning of features. The larger than average distance of 0.78 nm can be attributed to observations where complex coated ssDNA regions appear elongated [25]. e Overlay of a possible molecular model to match the observed topography [DNA model from [26] (1BNA)]

The phenomenon of sectoral dehybridization of DNA double strands, known as melting-, denaturation- or breathing bubbles, is well-known in DNA science [27]. It is crucial to DNA replication, by introducing insertion point for single strand interacting polymerase enzymes [28, 29] and is often connected to AT-rich (adenine, thymine) DNA regions, which are, due to two instead of three hydrogen bonds, slightly less strongly bound than their GC-counterparts [30] (guanine, cytosine).

They were observed indirectly (by fluorescence quenching [31, 32]) and directly (TEM [25] and AFM [33]), mainly in circular DNA and plasmids, where the openings are induced by twist stress. They appear also in linear DNA [34], caused statistically by molecular motion [27, 35,36,37].

The height of the measured protrusions lies with 0.2 nm well below the one measured for a double strand in air (1 nm) and vacuum (0.4 nm). In addition, the average interval of the protrusions is with 1.42 nm too large for a single nucleobase length and too low for a double helix pitch, but nicely matches the twofold distance between two neighboring phosphate groups. Therefore, we conclude that the “melted bubble” sections indeed consist of dehybridized double DNA strands into two single strands, where the beads representing the coordinatively bonded [(HtomMe)Cu2]3+ molecules. It shall also be noted that upon closer inspection, a zigzag type arrangement of antiphasic binding sites along the single strands can readily be discerned (see Fig. 2d).

Reasoning of the melting bubbles

With the dissociated sections making up approximately 12% of the contour length of a strand, which is even above the usual occurrence rate of single strand regions in linear DNA of 4% in supercoiled form [25], the observation of the formation here does not appear to be a purely statistical process.

Melting bubbles do not appear frequently under ambient conditions and are obviously not caused by the UHV conditions, rough handling or freezing as additional measurements have shown. To the best of our knowledge, it has not been shown in literature where DNA adsorbed to low adhesion substrates was subjected to vacuum [38, 39].

Since they do not structurally change over time, it can also be excluded that they are the result of tip-induced damage by the scanning motion. Therefore, a correlation to the presence of the metal-coordinated backbone binding complex molecules is comprehensible.

DNA behavior in vacuum

To correctly assess these findings, it is essential to take into account the behavior of dsDNA in vacuum environment, as well as the influence of the substrate.

It is widely known that dsDNA in environments of low humidity transitions into a more densely packed formation, the so-called A-form. This is caused by the lower availability of stabilizing water molecules and ions surrounding the negatively charged backbone. A lower helical pitch distance and a more pronounced major groove are characteristic for the A-form [40, 41].

Since it is the most common choice of substrate for investigations under ambient conditions, most AFM-works of DNA in vacuum conditions feature the DNA adsorbed onto functionalized mica. This heavily influences the structural properties of the adsorbed strands, whether it conserves the physiological B-form [15] or changes the elastical properties derived from persistence lengths and equilibrated geometries [42]. The reason is that in order to attach DNA to mica, it is necessary to use a functionalization, to avoid the repelling negative DNA charges of the phosphate backbone as well as of the mica (\({\text{SiO}}_{x}^{-}\)). Therefore, different methods are employed, that offer different affinities to the DNA. This ranges from less invasive approaches by using divalent cations [12] like Mg2+ or Ni2+, to the use of strong bonding surface modifications like APTES [42], ODA [34], PLL [43] or PLO [44]. The strength of the bonding influences the observed geometry, since the stronger the bonding, the less dynamics is allowed after adsorption. This leads to 2D equilibrated strands on ionic functionalization and 3D projected strands on covalent bonding surfaces. Electrically neutral substrates like Au or HOPG offer little to no attraction to DNA, disregarding mirror-dipoles and associated van-der-Waals forces, thus, making them mostly unsuitable for use under ambient and liquid conditions. In UHV on the other hand, gold is of great advantage, since it can be ultra-cleanly prepared and atomically controlled in-situ, while mica introduces a distinct roughness that is highly disadvantageous to high-resolution measurements.

Since literature is sparse, regarding the interaction of DNA with gold, its behavior during electrospray ionization and even its behavior in vacuum conditions, an appropriate series of tests of deposition on different substrates has been conducted.


A string of bead-like appearance characterizes the pristine DNA-strands on gold (Fig. 3a). The grooves between the features correlate to the major groove [38] while smaller features could not be resolved, due to instability of the sample. Another explanation features for the beads to be correlated to electrostatic forces from residual ions, which are known to reside primarily in the grooves [45, 46]. Caused by the low interaction between strand and substrate, the strands are easily moved during scanning, which prevents higher resolution. Figure 3d shows a distribution of measured groove distances. The mean value of 2.38 nm is in agreement to the expected pitch value range for A-form transformed DNA (2.82 nm [47]). The associated strand height measures to about 0.5 nm.

Fig. 3
figure 3

λ-DNA subjected to UHV conditions on different substrates. a Untreated DNA on Au(111) (df =− 6.2 Hz, A = 12.6 nm, f0 = 286.3 Hz). b [(HtomMe)Cu2]3+-treated DNA on Au(111) (df =− 2.0 Hz, A = 14.1 nm, f0 = 259.0 Hz). The residues on the otherwise flat atomic planes are likely due to residual complex molecules. c DNA on potassium bromide (df =− 8.6 Hz, A = 4.7 nm, f0 = 302.0 Hz, inset: df =− 20.2 Hz, A = 1.9 nm, f0 = 302.0 Hz). d Distribution of groove distances, along with corresponding mean values for different substrates (Au, KBr)

Potassium bromide

As a reference, pristine λ-DNA strands were additionally deposited on an in-situ cleaved surface of KBr. Strands can clearly be identified, with features that resemble the major groove, even despite signs of minor rearrangement of ion plateaus, possibly caused by residual solvent hitting the surface (Fig. 3c).

A statistical assessment of groove distances shows high agreement to the expected full turn distance of 3.4 nm of the B-form DNA-configuration.

For most parts, the right-handedness can be readily observed. However, upon closer look, segments of differing groove distance and even left-handedness, resembling A- and even the Z-form of dsDNA [41], can be additionally discerned. This is known to occur under high salt concentrations at GC-rich sections [48] and even on Ni-functionalized mica [49].

The expected vacuum-induced transition to the A-form is apparently hampered for the most part, due to an abundance of ions that act stabilizing to the physiological B-form and being able to retain a commensurate number of water molecules from evaporating [45].

The measured dsDNA height is about 0.6 nm, due to flattening from adhesion forces and embedding into residual solvent.

[(Htom Me )Cu 2 ] 3+ -bound DNA

Aside from the already discussed strand aggregations and bubble formations, the DNA bonded with [(HtomMe)Cu2]3+ exhibits a similar string of bead-like appearance, away from the larger aggregations (Fig. 3b). A statistical analysis of their groove distances leads to a mean value of 2.94 nm. While the pitch can slightly be influenced by the base sequence [48], this is far in between the untreated DNA on gold that has taken on the A-form and the presumably conserved B-form on KBr.

There were no differences in persistence lengths noticeable for each of the three cases, residing each time around 10 nm.


Freezing of B-conformation

As stated above, the complexes are designed in such way that they recognize and selectively bind to two consecutive phosphate moieties of the backbone of a DNA strand. For such a specific level of molecular recognition to a specific binding site geometry, not only the distance of the metal atoms is important, but also their metal-coordination site angles. Since the acetate groups, which are present in the unsolvated state, point outward at 110° (see Additional file 1: Figure S4), the binding could favor sites larger than the copper distance.

The measurements of untreated DNA showed that DNA tends to transition into a different configuration when not subjected to stabilizing factors. One of these stabilizing factors was shown to be salt ions, while they may not work directly but rather by binding residual water molecules necessary to stabilize the B-form. Here, the complexes appear to act in a similar fashion to the backbone of the DNA. Upon conformational change, a smaller phosphate distance would be highly unfavorable for the binding kinetics. Since the binding already occurred during the incubation in liquid environment and was shown to be irreversible [4], the adsorbed complexes act like a brace, forcing the phosphates to the former larger distance, and therefore partially conserving the B-conformation, at least to some degree, an effect well described for other compounds [50]. Incidentally, this effect might also be related to the high occurrence of melting bubbles.

Melting bubbles

The reason for the B-to-A conformational change can be found in the puckering of the furanose rings of the backbone. In the hydrated form they take on a C2’-endo form, whether without the influence of water molecules, the energetically favorable conformation is the C3’-endo form, which is characterized by a switched angle between the 2’- and 3’-carbon atoms (illustrated in Additional file 1: Figure S5). This altered angle leads to a reduced geometrical distance between the phosphates from 6.72 [26] to 5.85 Å [51], thus the denser packing in A-configuration [47, 52, 53].

Without the stabilizing influences of water molecules, the former state is now energetically unfavorable, but nonetheless enforced by the complexes. This way, an angular mechanical strain is introduced onto the base pairs, which are prone to dissociating due to rotational motions [54]. We argue that it is actually favored for specific segments to break up and melt into single strands in order to relief this mechanical stress.

Since there are hints that residual solvent hits the surface before evaporating completely, the molecules are already in contact with the substrate when complete evaporation occurs. For a transformation to happen now, it would require extensive on-surface dynamics, which the employed gold substrate indeed allows for. When bubbles open up, the overall twist sum of helical turns needs to be conserved. Still, after adsorption, the strands are somewhat restrained in their freedom of motion, despite the low adhesive force of the substrate. This would cause a local supercoiling buildup, which can indeed be observed locally (Fig. 2b).

To determine possible causes, a statistical analysis correlated the dissociation rate with the occurrence of specific conditions and geometric adsorption patterns. It showed that the percentage of passages of more than two subsequent AT-pairs in the λ-phage genome is 37.5% (Fig. 4a). Two base pairs corresponding here to the smallest observed gap distances, so the phenomenon cannot purely be related to the AT-pair occurrence.

Fig. 4
figure 4

a Contour length percentage along a strand when summing up dissociated bubble segments from a certain threshold value (purple). Segments with consequential AT base paring along the λ-phage genome sequence (blue). The percentage of segments with odd complex adsorption pattern, derived from a statistical simulation (yellow). Percentage of areas where both conditions are met at the same time (green), fitting well to the actual observed bubble lengths. The dip at 1 bp accounts for the visual counting threshold. Inset: Distribution of melting bubble lengths. The relation between base pair and nm of 1 bp \(\widehat{=}\) 0.71 nm was determined from the interval distance in Fig. 2c. b Explanation model for the bubble formation. Top: Zigzag adsorption pattern of complex molecules. Middle: Reduced phosphate distance by partial A-to-B transformation introduces tensile stress on the base pairing. Bottom: Dehybridization of weaker bound AT-base pairs

In contrast, with a purely statistical adsorption model, the occurrence of passages of more than two adversely packed complex molecules (equaling three base pairs) is roughly at 38.4%, which is far higher than the observed dissociated areas. When combining these two factors, simulating the adsorption of complex molecules onto the λ-phage genome, the percentage where both of these conditions occur at the same time is about 10%. This fits well with the 12% observed dissociation contour length.

This lets us conclude that an interplay of adsorption pattern and the base pair sequence is responsible for the observed formation of melting bubbles. Having both conditions occur simultaneously makes a strand segment energetically prone for dissociation, with tensile stress on the base pairs that would otherwise not be present when the adsorption pattern was parallel or random, as illustrated in Fig. 4b. Although this statistical method only gives a rough estimation, since it leaves out more complex effects like anisotropic stacking interaction [55] and cooperative area opening [37], it nicely reflects the overall findings.


Double stranded λ-phage DNA was successfully treated with dinuclear [(HtomMe)Cu2]3+ complexes upon selective metal coordination and deposited in UHV by electrospray deposition onto in-situ prepared Au(111) and investigated by AFM. Binding to the phosphate backbone could be confirmed, as well as the effect of backside π-stacking mediated strand agglomeration.

Interestingly, high amounts of melting bubbles were evidenced, where even individual complex molecules could be identified. The unfolding confirms the expected bonding scheme in contrast to other feasible bonding modes such a quasi-minor groove strand bridging. Aside from coverage of every two phosphates, the binding can take on even the closest expected packing, without any non-cooperative dynamics.

Compared to our previous series of experiments, in which we characterized binding phenomena of dinuclear complexes under ambient conditions, the UHV data presented here exhibit a resolution down to the single base pair level. The observed melting bubbles were explained by the bound complexes blocking the vacuum-induced B-A-transition by shortening of backbones, which leads to the mechanically preferable melting of AT-rich segments coupled with areas of odd packing, where single strand stability is energetically favored over conservation of the double strand conformation.

These quantitative experimental results, paralleled by statistical simulations impressively shade light on the rationale for strand dissociations of a novel DNA interaction process.

Besides these findings, the presented study demonstrates the value of such an approach of observing the behavior of individual molecules, employing the high-resolution capabilities of AFM, complementing common ensemble-based techniques. Going forward, this methodology can be applied for various chemotherapeutical agents, in order to gain detailed insights into their nanomechanistic properties, assess their effectiveness in blocking DNA-synthesis and potential to cause overall damage. This way aiding not only the development of new medicinal compounds, but also the associated clearance process through the federal drug agencies.


Complex synthesis

The synthesis of the [(HtomMe){Cu(OAc)}2](OAc) compound is described elsewhere [5]. The complex (M = 1158.22 g/mol) was dissolved in distilled water (MilliQ) at a concentration of 200 µM before this solution was added to the DNA solution. Incubation times were \(\sim\) 45 h.

Atomic force microscopy

The scanning force images were recorded on a RHK UHV 7500 system (RHK Technology, Troy, USA) with R9 controller electronics. Cantilevers used were Tap300Al-G (Budget Sensors, Sofia, Bulgaria) and NCH-PPP (Nanosensors, Neuchâtel, Switzerland), which were sputtered before use by Ar+ ions (8 × 10–7 mbar) at 680 eV for 90 s. Images were recorded in AFM non-contact imaging mode, at room temperature and at 90 K (cooled by a liquid nitrogen flow cryostat), at a pressure of 5 × 10–11 mbar. They were corrected for thermal drift [56], using Gwyddion SPM image software [57]. Persistence lengths were measured with Easyworm [58].

Sample preparation

The samples were prepared in-situ on a commercial substrate consisting of a 300 nm layer of Au(111) epitaxially grown onto a sheet of mica (Georg Albert PVD, Silz, Germany). They were precleaned in ethanol and water before transferring to UHV. There they were sputtered for 5 min with Ar+ ions at 8 × 10–7 mbar und subsequently annealed at 523 K for 2 h, for two cycles and checked for cleanliness before deposition of molecules.

Electrospray deposition

A solution of λ-phage DNA in MilliQ water (12.5 ng/µl, 396 µl) was treated with the complex solution (40 nmol/l, 4 µl) with a resulting complex concentration of 0.4 nmol/l and was diluted with methanol in a 2:1 ratio immediately before ESI deposition [59]. The system used for the procedure was a commercial MolecularSpray (MolecularSpray Ltd.) setup, with three individually pumped chambers of gradually decreasing pressure (1 × 10− 1 mbar, 1 × 10− 2 mbar, 1 × 10− 5 mbar). The solution is carried by a syringe through PEEK tubing to the capillary where it gets nebulized by a voltage of 2 kV and injected into the vacuum chamber. The syringe provided a flux of 10 µl/min, approximately 500 µl of sample solution was injected.


The statistical simulation was written in Matlab 2018b. The λ-phage genome sequence originates from [60]. Pseudocode is attached in Additional file 1.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.


  1. Pages BJ, Ang DL, Wright EP, Aldrich-Wright JR. Metal complex interactions with DNA. Dalton Trans. 2015;44:3505–26.

    Article  CAS  Google Scholar 

  2. Santini C, Pellei M, Gandin V, Porchia M, Tisato F, Marzano C. Advances in copper complexes as anticancer agents. Chem Rev. 2014;114:815–62.

    Article  CAS  Google Scholar 

  3. Johnstone TC, Suntharalingam K, Lippard SJ. The next generation of platinum drugs: targeted pt(II) agents, nanoparticle delivery, and pt(IV) prodrugs. Chem Rev. 2016;116:3436–86.

    Article  CAS  Google Scholar 

  4. Jany T, Moreth A, Gruschka C, Sischka A, Spiering A, Dieding M, Wang Y, Samo SH, Stammler A, Bögge H, von Fischer Mollard G, Anselmetti D, Glaser T. Rational design of a cytotoxic dinuclear Cu2 complex that binds by molecular recognition at two neighboring phosphates of the DNA backbone. Inorg Chem. 2015;54:2679–90.

    Article  CAS  Google Scholar 

  5. Jany T, Horstmann née Gruschka C, Bögge H, Stammler A, Glaser T. A series of dinuclear complexes with a flexible naphthalene-spacer and MOM-cleavage by pre-coordinated lewis acids. Z Anorg Allg Chem. 2015;641:2157–68.

    Article  CAS  Google Scholar 

  6. Simon J, Stammler A, Oldengott J, Bögge H, Glaser T. Proof of phosphate diester binding ability of cytotoxic DNA-binding complexes. Inorg Chem. 2020;59:14615–9.

    Article  CAS  Google Scholar 

  7. Simon J, Horstmann née Gruschka C, Mix A, Stammler A, Oldengott J, Bögge H, Glaser T. Evaluation of the binding mode of a cytotoxic dinuclear nickel complex to two neighboring phosphates of the DNA backbone. Dalton Trans. 2022;51:2863–75.

    Article  CAS  Google Scholar 

  8. Schwarzbich S, Horstmann née Gruschka C, Simon J, Siebe L, Moreth A, Wiegand C, Lavrentieva A, Scheper T, Stammler A, Bögge H, von Fischer Mollard G, Glaser T. Stronger cytotoxicity for cancer cells than for fast proliferating human stem cells by rationally designed dinuclear complexes. Inorg Chem. 2020;59:14464–77.

    Article  CAS  Google Scholar 

  9. Dallas A, Vlassov AV, Kazakov SA. Principles of nucleic acid cleavage by metal ions. In: Zenkova MA, editor. Artificial nucleases. Berlin: Springer berlin heidelberg; 2004.

    Google Scholar 

  10. García-Ramos JC, Galindo-Murillo R, Cortés-Guzmán F, Ruiz-Azuara L. Metal-based drug-DNA interactions. J Mex Chem Soc. 2013;57:245–59.

    Google Scholar 

  11. Liu C, Wang L. DNA hydrolytic cleavage catalyzed by synthetic multinuclear metallonucleases. Dalton Trans. 2009.

    Article  Google Scholar 

  12. Heenan PR, Perkins TT. Imaging DNA equilibrated onto mica in liquid using biochemically relevant deposition conditions. ACS Nano. 2019;13:4220–9.

    Article  CAS  Google Scholar 

  13. Pyne A, Thompson R, Leung C, Roy D, Hoogenboom BW. Single-molecule reconstruction of oligonucleotide secondary structure by atomic force microscopy. Small. 2014;10:3257–61.

    Article  CAS  Google Scholar 

  14. Maeda Y, Matsumoto T, Tanaka H, Kawai T. Imaging of the DNA (deoxyribonucleic acid) double helix structure by noncontact atomic force microscopy. Jpn J Appl Phys. 1999;38:L1211.

    Article  CAS  Google Scholar 

  15. Cerreta A, Vobornik D, Di Santo G, Tobenas S, Alonso-Sarduy L, Adamcik J, Dietler G. FM-AFM constant height imaging and force curves: high resolution study of DNA-tip interactions. J Mol Recog. 2012;25:486–93.

    Article  CAS  Google Scholar 

  16. Maeda Y, Matsumoto T, Kawai T. Observation of single- and double-stranded DNA using non-contact atomic force microscopy. Appl Surf Sci. 1999;140:400–5.

    Article  CAS  Google Scholar 

  17. Di Santo G, Tobenas S, Adamcik J, Dietler G. Imaging and force spectroscopy on DNA. Imag Micr. 2008;10:44–7.

    Google Scholar 

  18. Fenn JB, Mann M, Meng CK, Wong SF, Whitehouse CM. Electrospray ionization for mass spectrometry of large biomolecules. Science. 1989;246:64–71.

    Article  CAS  Google Scholar 

  19. Hinaut A, Pawlak R, Meyer E, Glatzel T. Electrospray deposition of organic molecules on bulk insulator surfaces. Beilstein J Nanotechnol. 2015;6:1927–34.

    Article  CAS  Google Scholar 

  20. Satterley CJ, Perdigão LMA, Saywell A, Magnano G, Rienzo A, Mayor LC, Dhanak VR, Beton PH, O’Shea JN. Electrospray deposition of fullerenes in ultra-high vacuum situ scanning tunneling microscopy and photoemission spectroscopy. Nanotechnology. 2007;18:455304.

    Article  Google Scholar 

  21. Wu X, Delbianco M, Anggara K, Michnowicz T, Pardo-Vargas A, Bharate P, Sen S, Pristl M, Rauschenbach S, Schlickum U, Abb S, Seeberger PH, Kern K. Imaging single glycans. Nature. 2020;582:375–8.

    Article  CAS  Google Scholar 

  22. Rauschenbach S, Stadler FL, Lunedei E, Malinowski N, Koltsov S, Costantini G, Kern K. Electrospray ion beam deposition of clusters and biomolecules. Small. 2006;2:540–7.

    Article  CAS  Google Scholar 

  23. Pawlak R, Vilhena JG, Hinaut A, Meier T, Glatzel T, Baratoff A, Gnecco E, Pérez R, Meyer E. Conformations and cryo-force spectroscopy of spray-deposited single-strand DNA on gold. Nat Com. 2019;10:1–7.

    Google Scholar 

  24. Wilcox DE. Binuclear metallohydrolases. Chem Rev. 1996;96:2435–58.

    Article  CAS  Google Scholar 

  25. Brack C, Bickle TA, Yuan R. The relation of single-stranded regions in bacteriophage PM2 supercoiled DNA to the early melting sequences. J Mol Bio. 1975;96:693–702.

    Article  CAS  Google Scholar 

  26. Drew HR, Wing RM, Takano T, Broka C, Tanaka S, Itakura K, Dickerson RE. Structure of a B-DNA dodecamer: conformation and dynamics. PNAS. 1981;78:2179–83.

    Article  CAS  Google Scholar 

  27. Altan-Bonnet G, Libchaber A, Krichevsky O. Bubble dynamics in double-stranded DNA. Phys Rev Lett. 2003;90:138101.

    Article  Google Scholar 

  28. Rieloff E, Nunes SCC, Pais AACC, Skepö M. Structural characterization of bubbles formed in DNA melting: a monte carlo simulation study. ACS Omega. 2017;2:1915–21.

    Article  CAS  Google Scholar 

  29. Tchernaenko V, Halvorson HR, Kashlev M, Lutter LC. DNA bubble formation in transcription initiation. Biochemistry. 2008;47:1871–84.

    Article  CAS  Google Scholar 

  30. Borer PN, Dengler B, Tinoco I, Uhlenbeck OC. Stability of ribonucleic acid double-stranded helices. J Mol Bio. 1974;86:843–53.

    Article  CAS  Google Scholar 

  31. Ambjörnsson T, Banik SK, Krichevsky O, Metzler R. Breathing dynamics in heteropolymer DNA. Biophys J. 2007;92:2674–84.

    Article  Google Scholar 

  32. Paul T, Bera SC, Mishra PP. Direct observation of breathing dynamics at the mismatch induced DNA bubble with nanometre accuracy: a smFRET study. Nanoscale. 2017;9:5835–42.

    Article  CAS  Google Scholar 

  33. Jeon J-H, Adamcik J, Dietler G, Metzler R. Supercoiling induces denaturation bubbles in circular DNA. Phys Rev Lett. 2010;105:208101.

    Article  Google Scholar 

  34. Dubrovin EV, Gerritsen JW, Zivkovic J, Yaminsky IV, Speller S. The effect of underlying octadecylamine monolayer on the DNA conformation on the graphite surface. Coll Surf B. 2010;76:63–9.

    Article  CAS  Google Scholar 

  35. Zeng Y, Montrichok A, Zocchi G. Length and statistical weight of bubbles in DNA melting. Phys Rev Lett. 2003;91:148101.

    Article  Google Scholar 

  36. Jeon J-H, Park PJ, Sung W. The effect of sequence correlation on bubble statistics in double-stranded DNA. J Chem Phys. 2006;125:164901.

    Article  Google Scholar 

  37. Wartell RM, Benight AS. Thermal denaturation of DNA molecules: a comparison of theory with experiment. Phys Rep. 1985;126:67–107.

    Article  CAS  Google Scholar 

  38. Tanaka H, Hamai C, Kanno T, Kawai T. High-resolution scanning tunneling microscopy imaging of DNA molecules on Cu(111) surfaces. Surf Sci. 1999;432:L611–6.

    Article  CAS  Google Scholar 

  39. Davies E, Teng KS, Conlan RS, Wilks SP. Ultra-high resolution imaging of DNA and nucleosomes using non-contact atomic force microscopy. FEBS Lett. 2005;579:1702–6.

    Article  CAS  Google Scholar 

  40. Wahl MC, Sundaralingam M. Crystal structures of A-DNA duplexes. Biopolymers. 1997;44:45–63.

    Article  CAS  Google Scholar 

  41. Ghosh A, Bansal M. A glossary of DNA structures from A to Z. Acta Cryst D. 2003;59:620–6.

    Article  Google Scholar 

  42. Adamcik J, Klinov DV, Witz G, Sekatskii SK, Dietler G. Observation of single-stranded DNA on mica and highly oriented pyrolytic graphite by atomic force microscopy. FEBS Lett. 2006;580:5671–5.

    Article  CAS  Google Scholar 

  43. Brouns T; Keersmaecker de H; Konrad S F; Kodera N; Ando T; Lipfert J; Feyter S de; Vanderlinden W. 2018 Free energy landscape and dynamics of supercoiled dna by high-speed atomic force microscopy. ACS Nano. 12 11907–11916

  44. Podesta A, Imperadori L, Colnaghi W, Finzi L, Milani P, Dunlap D. Atomic force microscopy study of DNA deposited on poly-l-ornithine-coated mica. J Micr. 2004;215:236–40.

    Article  CAS  Google Scholar 

  45. Yoo J, Aksimentiev A. Competitive binding of cations to duplex DNA revealed through molecular dynamics simulations. J Phys Chem B. 2012;116:12946–54.

    Article  CAS  Google Scholar 

  46. Lipfert J, Doniach S, Das R, Herschlag D. Understanding nucleic acid-ion interactions. Annu Rev Biochem. 2014;83:813–41.

    Article  CAS  Google Scholar 

  47. Wood BR. The importance of hydration and DNA conformation in interpreting infrared spectra of cells and tissues. Chem Soc Rev. 2016;45:1980–98.

    Article  CAS  Google Scholar 

  48. Kanno T, Tanaka H, Miyoshi N, Fukuda M, Kawai T. Base sequence dependence of deoxyribonucleic acid studied by scanning tunneling microscopy. Jpn J Appl Phys. 1892;2000:39.

    Google Scholar 

  49. Leung C, Bestembayeva A, Thorogate R, Stinson J, Pyne A, Marcovich C, Yang J, Drechsler U, Despont M, Jankowski T, Tschöpe M, Hoogenboom BW. Atomic force microscopy with nanoscale cantilevers resolves different structural conformations of the DNA double helix. Nano Lett. 2012;12:3846–50.

    Article  CAS  Google Scholar 

  50. Fritzsche H, Brandes R, Rupprecht A, Song Z, Weidlich T, Kearns DR. The formation of A-DNA in NaDNA films is suppressed by netropsin. Nucl Acid Res. 1992;20:1223–8.

    Article  CAS  Google Scholar 

  51. Bingman CA, Zon G, Sundaralingam M. Crystal and molecular structure of the A-DNA dodecamer d(CCGTACGTACGG). J Mol Bio. 1992;227:738–56.

    Article  CAS  Google Scholar 

  52. Saenger W, Hunter WN, Kennard O. DNA conformation is determined by economics in the hydration of phosphate groups. Nature. 1986;324:385–8.

    Article  CAS  Google Scholar 

  53. Sundaralingam M, Westhof E. The, “rigid” nucleotide concept in perspective. Int J Quant Chem. 1979;16:115–30.

    Article  Google Scholar 

  54. Yakushevich LV. Nonlinear physics of DNA. Weinheim: Wiley-VCH; 2004.

    Book  Google Scholar 

  55. Protozanova E, Yakovchuk P, Frank-Kamenetskii MD. Stacked-unstacked equilibrium at the nick site of DNA. J Mol Bio. 2004;342:775–85.

    Article  CAS  Google Scholar 

  56. Rahe P, Bechstein R, Kühnle A. Vertical and lateral drift corrections of scanning probe microscopy images. J Vac Sci Tec B. 2010;28:C4E31-8.

    Article  CAS  Google Scholar 

  57. Nečas D, Klapetek P. Gwyddion: an open-source software for SPM data analysis. Open Physics. 2012;10:99.

    Article  Google Scholar 

  58. Lamour G, Kirkegaard JB, Li H, Knowles TP, Gsponer J. Easyworm: an open-source software tool to determine the mechanical properties of worm-like chains. Source Code Bio Med. 2014;9:16.

    Article  Google Scholar 

  59. Olumee Z, Callahan JH, Vertes A. Droplet dynamics changes in electrostatic sprays of methanol−water mixtures. J Phys Chem A. 1998;102:9154–60.

    Article  CAS  Google Scholar 

  60. Sanger F, Coulson AR, Hong GF, Hill DF, Petersen GB. Nucleotide sequence of bacteriophage λ DNA. J Mol Bio. 1982;162:729–73.

    Article  CAS  Google Scholar 

Download references


Open Access funding enabled and organized by Projekt DEAL. We gratefully acknowledge financial support from the German Science Foundation (DFG), under the grants DA 370/8-1 and GL 295/8-1.

Author information

Authors and Affiliations



SS and TG designed and performed the synthesis of the complex molecules. NB performed all measurements and analyzed the data. DK and VW contributed substantially to sample preparation and interpretation of results. NB and DA wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Niklas Biere.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

: Figure S1. Various stages of the (HtomMe)- complexes. (a) The basic form of the H2tomMe ligand. (b) A methoxymethyl(MOM)-based precursor for synthesis of the metal complexes. (c) The cationic form crystallized with an acetate OAc counter anion. (d) In aqueous solution the acetate groups are replaces by H2O. (e) The copper ions coordinate to the phosphates of the DNA backbone. Figure S2. λ-DNA subjected to dinuclear copper complexes at various concentrations. Images taken in air in tapping mode. Same length and height scales apply to each image. (a) Untreated λ-DNA (12.5 ng/µl), strand height 0.6 nm. (b) Final complex concentration 0.7 µM, increase of strand height to 1 nm. (c) 1.4 µM, height 0.8 nm. (d) 2 µM, increase of coagulation (indicated by arrows), height 1 nm. (e) 7 µM, height 1 nm. (f) 14 µM, complete coagulation, height 0.4 - 1 nm. (g) 0.4 nM with an incubation time of 45 h. (h) 4 nM. Figure S3. Empirical determination of a suitable complex concentrations: measurements of electrospray depositions of complex-treated DNA under UHV. (a) At 2 µM, a massive residual background is present, the strands appear visibly embedded into a layer (df = - 2.2 Hz, A = 36.1 nm, f0 = 257.4 Hz). (b) At 0.4 nM, still a slight residual molecular background is visible, the gold surface appears to be covered (df = - 2.5 Hz, A = 14.1 nm, f0 = 259.0 Hz). Although technically there should be less complexes than binding sites to the DNA, due to the obvious excess of molecules, it can be assumed that the strands are maximally coated. Figure S4. The [(HtomMe){Cu(OAc)}2]+ complex chemical structure (hydrogen left out) from different perspectives, based on data from [1,2]. (a) Angular view, (b) side view, (c) top view. Figure S5. The puckering of the furanose ring of the DNA backbone changes in humid environments. The (*)-marked C2’ atom switches from the front into the back of the ring plane. The designations C2’- and C3’-endo refer to the atom which is positioned on the side of the C5’ atom. Structural X-ray diffraction data from [3] (1BNA, B-DNA) and [4] (116D, A-DNA). Figure S6. Pseudocode for the analysis of the λ-genome sequence. Firstly, the sequence is reduced for easier counting, then the counting of segments in dependence of the threshold value length occurs. Figure S7. Pseudocode for the adsorption simulation and subsequent analysis. The adsorption routine is depicted in black. Two arrays of singular binding spots, representing the two single strands of one double helix, get subsequently and randomly filled up until no free spots remain. Since the molecules require two adjacent spots, full adsorption will leave randomly distributed singular free spots. In yellow and green depicted are the analysis routines, according to their colors in the graph of Figure 4 in the main paper. The adsorbed strand is reduced for easier identification of molecules. Then two counters are iterated over the arrays, each counting one zigzag phase. The exit conditions for the counting of one segment are that both spots of the same index are either empty or both occupied, meaning the zigzag arrangement would be interrupted either way. In addition, when a GC-pair is present (green).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Biere, N., Kreft, D., Walhorn, V. et al. Dinuclear complex-induced DNA melting. J Nanobiotechnol 21, 26 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • DNA
  • AFM
  • Molecular recognition
  • Electrospray ionization
  • Biomolecules