Copper nucleases as a famous class of artificial metallonucleases have attracted considerable interest in relation to their diverse potentials not only as therapeutic agents but also in genomic researches. Copper nucleases present high efficient oxidative cleavage of DNA, in which DNA strand scission occurs generally after hydrogen atom abstracted from a sugar moiety. In order to achieve the selective cleavage of DNA sequences by copper nucleases, the DNA specific recognition agents of the Dervan-type hairpin and cyclic polyamides can be considered as proper carriers of copper nucleases. Investigation of the DNA cleavage selectivity of copper nucleases assisted by the hairpin and cyclic polyamides at the molecular level has not yet been elucidated.
We carried out a series of molecular dynamics simulations for the nuclease [Cu(BPA)]2+ or [Cu(IDB)]2+ bound to the hairpin/cyclic polyamide and associated with DNA to investigate the selective DNA cleavage properties of Cu(II)-based artificial nucleases. The simulated results demonstrate that the DNA cleavage selectivity of the two nucleases assisted by the hairpin polyamide is improved efficiently. The [Cu(BPA)]2+ or [Cu(IDB)]2+ nuclease with a substrate OOH- bound to the hairpin polyamide can be stably located at the minor groove of DNA, and possibly abstracts H atom from the sugar of DNA. However, the DNA cleavage properties of the two nucleases assisted by the cyclic polyamide are significantly poor due to the rigidity of linking region between the cyclic polyamide and nuclease. With introduction of the flexible linker -CH2CH2CH2NH2, the modified cyclic polyamide can assist the two copper nucleases to improve the selective DNA cleavage properties efficiently.
A flexible linker and a proper binding site of the polyamide-type recognition agents play an important role in improving the DNA cleavage selectivity of copper nucleases. Current investigations provide an insight into the DNA cleavage specificities of chemical nucleases assisted by an appropriate nucleic acid recognition agent.
Artificial nucleases have attracted considerable attention for their diverse applications not only as therapeutic agents but also in biochemical research [1–4]. Among them, transition metal complexes have been studied extensively due to their diversity in structure and reactivity [5–8]. Copper complexes with the biological and accessible oxidative/reductive potentials have become a class of the most frequently studied artificial metallonucleases [9–16]. The first copper nuclease, [Cu(OP)2]2+ (OP = 1,10-phenanthroline) [17, 18], with high efficiency DNA cleavage property, motivated decade-long studies on mononuclear copper nucleases [3, 5, 19–22]. Especially, a great number of mononuclear copper nucleases[23–26] with BPA (BPA = bis(2-pyridylmethyl) amine) and IDB (IDB = N,N-bis(2-benzimidazolylmethyl) amine) ligands have attracted much attention due to their synthetic accessibility, low molecular weight and efficient DNA cleavage ability. Moreover, some experiments and theoretical calculations have predicted that the mononuclear copper nucleases present high efficient oxidative cleavage of DNA [7, 23, 29–33], in which DNA strand scission occurs generally after hydrogen atom abstractions of deoxyribose sugar moiety [2, 27, 34]. However, the oxidative DNA cleavage process only involves the hydrogen atom abstraction from any sugar in DNA, which makes most copper nucleases do not selectively recognize and cleave DNA sequences [35–37].
In order to improve the selective cleavage of copper nucleases to DNA molecules, specific DNA sequence recognition agents were introduced as carriers of copper nucleases [23, 38]. In the case of DNA recognition agents, polyamides have presented promising characteristics due to its simple structure, synthetic accessibility and sequence specific affinity to base pairs in the DNA minor groove [39–42]. Recently, the experimental observations reported by Aldrich-Wright and co-workers demonstrated that the complexes of the Dervan-type polyamides combined with platinum compounds can successfully recognize the relevant sequences of DNA [43, 44]. Moreover, among the Dervan-type polyamides, the hairpin and cyclic polyamides have presented higher affinities and specificities to DNA sequences than those of single-chain and antiparallel-double polyamides [1, 40, 45, 46]. Specifically, eight-ring hairpin polyamide with the flexible β-alanine (β) residue that binds to DNA as a "turn-to-tail" model possesses excellent affinity and sequence specificity of DNA molecule [6, 47]. The cyclic polyamides with γ-butyric acid (γ) can also efficiently recognize DNA sequences by a "turn-to-turn" model [48–50]. In the present work, we have performed a molecular dynamics (MD) study on the interactions of DNA with the ligands formed by each of the high efficient DNA-cleavage[22, 51] copper nucleases, [Cu(BPA)]2+ and [Cu(IDB)]2+, and by either the hairpin Dervan-type polyamide, ImImPyPy-γ-PyPyPyPy-β-Dp, or the cyclic Dervan-type polyamide, . Namely, six independent simulations were carried out, of which first four simulations were performed on hairpin-polyamide + [Cu(BPA)]2+- DNA, cyclic-polyamide + [Cu(BPA)]2+- DNA, hairpin-polyamide + [Cu(IDB)]2+- DNA and cyclic-polyamide + [Cu(IDB)]2+- DNA, to study the influences of different recognition agents and nuclease cleavage agents on the hydrogen abstractions from the sugar moiety of DNA; and the last two MD simulations were performed on modified-cyclic-polyamide and modified-cyclic-polyamide to investigate an improvement of the DNA cleavage selectivity of metal nucleases assisted by the modified cyclic polyamide with a flexible linker -CH2CH2CH2NH2.
The two polyamide-DNA systems are taken from the X-ray crystal structures of polyamide-DNA complexes. One is hairpin polyamide-DNA complex, ImImPyPy-γ-PyPyPyPy-β-Dp-d(AATATCCACCTGCA)2 (PDB code: 1M1A) ; and the other is cyclic polyamide-DNA complex, - d(CGCTAACAGGC)2 (PDB code: 1PQQ) , assigned as HPD and CPD, respectively. The branch C atom of [Cu(BPA)]2+ or [Cu(IDB)]2+ nuclease was bound manually to the tail N atom (sp3 hybridization) of each polyamide chain in HPD and CPD, which employed the protocols studied by Aldrich-Wright and co-workers for the polyamide - platinum(II) complexes [43, 44]. Then the substrate OOH- was introduced to the [Cu(BPA)]2+ or [Cu(IDB)]2+ nuclease to form an oxygen bridge for addressing the cleavage possibility. Finally, each whole system was explored to AutoDock 3.0 program for selecting the initial structure for the MD simulation. However, the position of N-tail in CPD lies in the distal Py residue of cyclic polyamide, which makes the orientation of the linking site point away from the DNA. For selecting an appropriate linking site, the group of -NH2 (as a linker) on the C atom (sp3 hybridization) of γ region was introduced to form the cyclic polyamide of associated with DNA, assigned as CPDγ, which was suggested by the previous experimental observations [55, 56]. Moreover, to improve the flexibility of liking region between the cyclic polyamide and nuclease, and to increase better the possibility of sugar hydrogen abstraction, the flexible linker of -CH2CH2CH2NH2 on the C atom (sp3 hybridization) of γ region was introduced to form the modified cyclic polyamide of associated with DNA, assigned as CPDγTail. The component sketches of six polyamide+copper nuclease+OOH- ligands are shown in the additional files 1A, B, C, D, E and 1F. Therefore, HPD + [Cu(BPA)OOH]+, CPDγ+ [Cu(BPA)OOH]+ and CPDγTail + [Cu(BPA)OOH]+ were assigned as HPD-BPA, CPDγ-BPA and CPDγTail-BPA, respectively. Similarly, HPD + [Cu(IDB)OOH]+, CPDγ+ [Cu(IDB)OOH]+ and CPDγTail + [Cu(IDB)OOH]+ were assigned as HPD-IDB, CPDγ-IDB and CPDγTail-IDB, respectively. Given that each strand of DNA has some phosphate groups, sodium ion counterions were separately added to each system to achieve electroneutrality. The systems were solvated explicitly by using the TIP3P water potential inside a central simulation box. The box dimensions ensure the solvent shell extended to10 Å in all directions of each system studied.
Molecular dynamics simulations
To take advantage of previous extensive simulations [44, 54, 57–60], the protocols employed therein were directly adapted in this study. The six MD simulations were carried out using the AMBER9 package with the AMBER force fields of parm99[61, 62] and gaff . The atomic types for the studied polyamides, OOH-, [Cu(BPA)]2+ and [Cu(IDB)]2+, except for copper atoms surrounding of [Cu(BPA)]2+ and [Cu(IDB)]2+ nucleases, were generated by using the ANTECHAMBER module included in AMBER9 program package. The force field parameters of copper center of [Cu(BPA)OOH]+ and [Cu(IDB)OOH]+ were evaluated based on the quantum chemical calculations in our previous work reported elsewhere [30, 33, 57]. The electrostatic potentials for RESP calculations have been calculated at the HF/6-31G** level of theory[64–66] using Gaussian03 program . The RESP charges of the polyamides, OOH-, [Cu(BPA)]2+ and [Cu(IDB)]2+ were derived by the RESP program based on the calculated electrostatic potentials.
The protocol for all MD simulations was described herein as follows: (1) The systems were energetically minimized to remove unfavorable contacts. Four cycles of minimizations were performed as follows, i.e. 5000 steps of each minimization with harmonic restraints from 100, 75, 50 to 25 kcal·mol-1·Å-2, which means that the restraints were relaxed stepwise by 25 kcal·mol-1·Å-2 per cycle, on DNA, polyamide and nuclease positions. The fifth cycle consists of 10000 steps of unrestrained minimization before heating process. The cutoff distance used for the nonbonded interactions was 10 Å. The SHAKE algorithm  was used to restrain the bonds containing hydrogen atoms. (2) Each energy-minimized structure was heated over 120 ps from 0 to 300 K (with a temperature coupling of 0.2 ps), while the positions of DNA, polyamide and nuclease were restrained with a small value of 25 kcal·mol-1·Å -2. The constant volume was maintained during the processes. (3) The unrestrained equilibration of 200 ps with constant pressure and temperature conditions was carried out for each system. The temperature and pressure were allowed to fluctuate around 300 K and 1 bar with the corresponding coupling of 0.2 ps, respectively. For each simulation, the integration step of 2 fs was used. (4) Finally, production runs of 30 ns were carried out by following the same protocol. A time of 200 ps after thermal warm-up in each simulation was selected as a starting point for data collections. During the production run, 15000 structures for each simulation were saved for post-processing by uniformly sampling the trajectory.
Results and Discussion
With the help of thirty nanosecond-long MD simulations, it was noted that these systems in water environment can be successfully simulated by using the protocol described above. The energies and root-mean-square deviations (RMSDs) for each simulation have been examined to clarify if each system had attained equilibrium. The RMSDs of each studied system with respect to its starting structure have attained equilibrium after 5 ns. The energies were found to be stable during the course of each remainder simulation. Hence, a time of 5 ns after thermal warm-up in each simulation was selected as a starting point for data collection and further analysis. The trajectory analysis for each system involves extracting the equilibrated conformations between 5 ns and 30 ns of simulation time, recording 12500 snapshots at every 2 ps time-interval of each trajectory. For the analyses of trajectories, the PTRAJ module of the AMBER9 program was used to extract production conformations.
Simulations of polyamide + [Cu(BPA)OOH]+-DNA complexes (HPD-BPA and CPDγ-BPA)
To investigate the DNA cleavage specificities of the metal nuclease assisted by different recognition agents, the [Cu(BPA)OOH]+ ligand was bound to hairpin and cyclic polyamides, respectively. The two polyamide + [Cu(BPA)OOH]+ - DNA complexes are assigned as HPD-BPA and CPDγ-BPA. The small RMSD values of these two systems predict that the interactions between the polyamide ligand and the DNA molecule in each system are strong enough to hold firmly the [Cu(BPA)OOH]+ ligand in the minor groove of DNA. Namely, the average RMSD values of two systems derived from the simulations over 12500 structures in the computed trajectories for HPD-BPA and CPDγ-BPA are respective 2.38 and 2.31 Å, as shown in Figures 1A and 1B. The mean-square fluctuations correlated with temperature factors for the systems support that each system reaches a relatively stable situation with small RMSD changes during the course of the simulation. The radial distribution functions (RDF) between the Od atom of OOH- ligand and the hydrogen atoms of sugar of DNA may supply primary information regarding DNA strand scission possibility. The RDFs, probability distributions and integration plots of probabilities of the Od-C1'H and Od-C4'H distances of the nearest cytosine and thymine bases, are shown in Figures 1C, E and 1G, respectively, for HPD-BPA. The RDFs and the probability distributions of the Od-C1'H and Od-C4'H distances of the nearest cytosine and thymine bases are shown in Figures 1D, and 1F, respectively, for CPDγ-BPA. Figures 1C and 1D show the first sharp RDF peak centered at 2.85 Å corresponding to a direct interaction of OOH- with C4'H atom for HPD-BPA, and that centered at 5.98 Å with the C4'H atom for CPDγ-BPA, respectively. In addition, Figures 1E and 1F present the probability distributions of the Od-C4'H distance less than 3.5 Å with 76.5% of the simulation time versus those of the Od-C1'H distance less than 3.5 Å with 16.5% of the simulation time for HPD-BPA, and the probability distributions of the Od-sugar hydrogen atom distances less than 3.5 Å all with 0% of simulation times for CPDγ-BPA, respectively. The integration plots of probabilities in Figure 1G show that the C4'H abstraction is more favorable than the C1'H one for HPD-BPA. However, Figure 1F presents the properties of the poor hydrogen abstractions from the sugars of DNA for CPDγ-BPA.
Analyses of the average structures of the trajectories obtained from the simulations of HPD-BPA and CPDγ-BPA reveal spatial details of the interactions of the entire ligands with DNA molecules, which supports the general observations described above. Specifically, (1) for HPD-BPA, the cleavage agent of [Cu(BPA) OOH]+ is located at the minor grove of DNA during the course of the simulation; the equatorial plane of the [Cu(BPA)]2+ nuclease is parallel to the wall of the minor groove of DNA with the angle of 12° between the plane of nuclease and the wall of minor groove, as shown in Figure 2A. However, for CPDγ-BPA, the cleavage agent of [Cu(BPA) OOH]+ cannot remain in the minor groove, whose orientation is removed from the minor groove of DNA, resulting from the rigidity of the N-tail of cyclic polyamide that constrains the flexibility of the [Cu(BPA)]2+ nuclease, as shown in Figure 2B. (2) The chains of polyamide for each system also remain almost parallel to the wall of the minor groove of DNA; the N-H groups nearby Im/Py of polyamides play an important role in acting as H-bonds donors to the N or O atoms of nucleobases located on the floor of the minor groove of DNA. Moreover, the H atom of the branch C atom in [Cu(BPA)]2+ and the N or O atom of bases of neighboring DNA can also form H-bonds. The affinity of polyamides with DNA are enhanced mainly by an increase in the number and strength of observed H-bonds between the polyamide-nuclease ligand and the DNA molecule, i.e., the H-bonds formed between the C-H groups of [Cu(BPA)]2+ and O4 atom of T8 base in DNA with C-O distance of about 3.1 Å. (3) The average structure of each system presents a relatively stable spatial conformation during the course of simulation. Figures 2A and 2B show the average structures from the simulations for HPD-BPA and CPDγ-BPA, respectively. The average structure with the average Od-C4'H distance of 2.78 Å for HPD-BPA in Figure 2C predicts the DNA cleavage possibility at the tail region of the polyamide, because that the Cu atom of [Cu(BPA)]2+ approaches the nearest sugar C4'H of C10 base. However, from Figure 2B for CPDγ-BPA, it is obvious that the structural characteristics reveal theoretically the poor cleavage ability of [Cu(BPA)]2+ nuclease due to the restraint of the rigid linker of the cyclic polyamide bound to the [Cu(BPA)]2+ nuclease. Namely, the restraint of DNA cleavage of the [Cu(BPA)]2+ nuclease results from the short linker between the branch C atom of nuclease and backbone of cyclic polyamide.
Simulations of polyamide + [Cu(IDB)OOH] +- DNA complexes (HPD-IDB and CPDγ-IDB)
For the [Cu(IDB)OOH]+ ligand, the average RMSD values derived from the simulations over 12500 structures in the computed trajectories for HPD-IDB and CPDγ-IDB are 2.70 Å and 2.40 Å, respectively, as shown in Figures 3A and 3B. It has been tested that the mean-square fluctuations correlated with temperature factors for the systems are small, which predicts that the systems have attained simulation equilibriums. Figures 3C and 3D show the first sharp RDF peak, centered at 2.85 Å and 6.05Å, corresponding to direct and nonexistent interaction of OOH- with the C4'H atom of sugar for HPD-IDB and CPDγ-IDB, respectively. In addition, Figures 3E and 3F present the probability distributions of Od-C4'H distance less than 3.5 Å with 34.7% of the simulation time versus those of the Od-C1'H distance less than 3.5 Å with 0% of the simulation time for HPD-IDB, and the probability distributions of the Od-sugar hydrogen atom distances of less than 3.5 Å all with 0% of the simulation time for CPDγ-IDB, respectively. The integration plots of probabilities in Figure 3G show that the C4'H abstraction is more favorable than the C1'H one for HPD-IDB. However, Figure 3F shows the properties of poor hydrogen abstractions for CPDγ-IDB.
The average structure of each trajectory obtained from the corresponding simulation supports the general observations described above. The equatorial plane of [Cu(IDB)]2+ nuclease in HPD-IDB is located perpendicularly at the wall of the minor groove of DNA with the angle of 80° between the plane of nuclease and the wall of minor groove, as shown in Figure 4A. However, for CPDγ-IDB as shown in Figure 4B, the [Cu(IDB)OOH]+ ligand cannot be located at the minor groove of DNA, whose orientation is away from the minor groove of DNA, resulting from the rigidity of the N-tail of cyclic polyamide constraining the flexibility of [Cu(IDB) OOH]+ ligand. This observation is similar to the properties of the cleavage agent of [Cu(BPA)OOH]+ discussed above. There are also H-bonds formed between the N-H groups nearby Im/Py from the polyamides and the N or O atoms of nucleobases located on the floor of the minor groove of DNA (i.e., C...N distance of 2.52 Å), and between the H atom of the branch C atom in [Cu(IDB)]2+ and the N/O atom of neighboring bases of DNA (i.e., C...O distance of 3.05 Å ). The average structures of HPD-IDB and CPDγ-IDB are shown in Figures 4A and 4B. The average structure from the simulation for HPD-IDB, as shown in Figure 4C, predicts the DNA cleavage possibility at the tail region of the polyamide, resulting from the Cu atom of [Cu(BPA)]2+ approaching the nearest sugar C4'H of C10 base with the average Od-C4'H distance of 2.75 Å. However, it is obvious that the structural characteristics of CPDγ-IDB reveal theoretically the poor cleavage ability of [Cu(IDB)OOH]+ ligand, which is similar to that of [Cu(BPA)OOH]+ ligand discussed above.
Based on the discussion above, the two systems of HPD-BPA and HPD-IDB with different nuclease cleavage agents present similar cleavage properties to DNA. Namely, the [Cu(BPA)OOH]+ and [Cu(IDB)OOH]+ nucleases can be stably located at the minor groove of DNA, though the orientation of each cleavage agent in the minor groove of DNA is different. The [Cu(BPA)OOH]+ nuclease is parallel to the wall of the minor groove of DNA, and the [Cu(IDB)OOH]+ nuclease is located perpendicularly at the minor groove of DNA. On the other hand, for hairpin-polyamide recognition, the probability of hydrogen abstraction of the copper nuclease ligand [Cu(BPA) OOH]+ from DNA is similar to that of the [Cu(IDB)OOH]+ ligand, i.e., the C4'H abstraction from the sugar of DNA is favorable. However, for the two systems of CPDγ-BPA and CPDγ-IDB, both nuclease ligands cannot cause DNA strand scissions due to the inaccessibility of any hydrogen of the sugar moiety of DNA, which suggests that the DNA cleavage properties of [Cu(BPA)]2+ and [Cu(IDB)]2+ nucleases assisted by the cyclic polyamide are poor. In addition, though the recognition property of the cyclic polyamide with DNA is better than that of the hairpin polyamide,[48, 53, 56] the [Cu(BPA)]2+ or [Cu(IDB)]2+ nuclease assisted by the cyclic polyamide cannot achieve the specific DNA strand scission.
Simulations of modified-cyclic-polyamide + [Cu(BPA)OOH]+/[Cu(IDB)OOH]+-DNA complexes (CPDγTail-BPA and CPDγTail -IDB)
The discussions above indicate that the nuclease ligands for both systems of CPDγ-BPA and CPDγ-IDB cannot cause any hydrogen abstractions from the sugar moiety of DNA due to the rigidity of the N-tail of cyclic polyamide. It is obvious that the linker between the branch C atom of nuclease and the backbone of polyamide is very important for specific cleavage of metal nucleases to DNA. To improve the selective cleavage properties of metal nucleases assisted by cyclic polyamide, i.e., to increase the possibility of hydrogen abstraction from the sugar moiety, the flexible linker of -CH2CH2CH2NH2 was introduced to the cyclic polyamide in order to form the modified cyclic polyamide bound by each of the two nucleases, as shown in the additional files 1E and 1F. The corresponding complexes, modified-cyclic-polyamide + [Cu(BPA)OOH] +/[Cu(IDB)OOH] + - DNA, assigned as CPDγTail-BPA and CPDγTail-IDB, are built for the simulation calculations. The analyses of RMSDs of the two systems with respect to their initial structures, RDFs, probability distributions and integration plots of probabilities of Od-C1'H and Od -C4'H distances of the nearest bases (i.e., the cytosine in CPDγTail-BPA and the thymine in CPDγTail-IDB) are shown in Figures 5A, C, E and 5G for CPDγTail-BPA and Figures 5B, D, F and 5H for CPDγTail-IDB, respectively. Figures 5A and 5B show that the simulations initiated from the docked structures correspond closely to the average structures for both CPDγTail-BPA and CPDγTail-IDB systems. The small RMSD values predict that the interaction between the polyamide ligand and DNA molecule in each system is strong enough to hold the small [Cu(BPA)OOH]+ and [Cu(IDB)OOH]+ ligands in the minor groove of DNA. The average RMSDs of two systems derived from the simulations over 12500 structures in the computed trajectories for CPDγTail-BPA and CPDγTail-IDB are respective 2.27 and 2.14 Å (see Figures 5A and 5B), which supports that each system reaches a relatively stable situation during the course of the simulation. Figures 5C and 5D show the first sharp RDF peaks centered at 2.65 Å and 2.58 Å corresponding to direct interactions of Od atoms of OOH- ligands with the C1'H atoms of sugars for CPDγTail-BPA and CPDγTail-IDB, respectively. In addition, Figure 5E presents the probability distributions of the Od -C1'H and Od-C4'H distances less than 3.5 Å with 72.3% and 15.8% of the simulation time for CPDγTail-BPA, respectively. Figure 5F presents the probability distributions of the Od -C1'H and Od-C4'H distances less than 3.5 Å with 62.5% and 25.2% of the simulation time for CPDγTail-IDB, respectively. The integration plots of probabilities in Figures 5G and 5H show that the numbers of C1'H atoms presented around Od-H distance of 3.5 Å are more than the numbers of C4'H atoms for these two systems. The results of probability distributions indicate that C1'H abstractions are favorable for both CPDγTail-BPA and CPDγTail-IDB systems. However, the possibility of C4'H or other sugar hydrogen atom abstractions is quite small.
Visual analyses of the average structures obtained from trajectories of CPDγTail-BPA and CPDγTail-IDB reveal spatial details of the interactions of the entire ligands with DNA molecules, which supports the general observations described above. Specifically, (1) the modified cyclic polyamide + [Cu(BPA)OOH]+ ligand or the modified cyclic polyamide + [Cu(IDB)OOH]+ ligand remains stable in the minor groove of the corresponding DNA molecule; (2) the equatorial plane of [Cu(BPA)]2+ nuclease is located perpendicularly at the wall of the minor groove of DNA with the angle of 72° between the plane of nuclease and the wall of minor groove of DNA for CPDγTail-BPA (see Figure 6A). The corresponding plane in CPDγTail-IDB is parallel to the wall of the minor groove of DNA with the angle of 28° between them (see Figure 6B); (3) for CPDγTail-BPA, the [Cu(BPA)OOH]+ ligand is located at the small region of C8/G8 dyad based on the connection of the branch C atom of [Cu(BPA)]2+ with the tail N atom of . The position of the Od atom of [Cu(BPA) OOH]+ is close to the nearest sugar C1'H of G8 base with the average Od-C1'H distance of 2.52 Å, as shown in Figure 6C. At the same time, H-bonds between the C-H groups of [Cu(BPA)OOH)]+ and O4/N2 atoms of C8/G8 bases in DNA with C...O or C...N distance of about 3.12 Å are observed in the data analysis. For the simulation of CPDγTail-IDB, the Od atom of the [Cu(IDB)OOH]+ ligand approaches the nearest sugar C1'H atom of T4 base at the tail of with the average Od -C1'H distance of 2.75 Å, as presented in Figure 6D. At the same time, H-bond events between the C-H groups of [Cu(IDB)OOH]+ and O2/N2 atoms of T4/A4 bases in DNA occur with the C...O or C...N distance of about 3.25 Å.
In summary, the simulated results of CPDγTail-BPA and CPDγTail-IDB present the stable conformations of the [Cu(BPA)OOH]+ and [Cu(IDB)OOH]+ ligands respectively bound to the N-tail of the modified cyclic polyamide associated with DNA. With the -CH2CH2CH2NH2 linker introduced, the [Cu(BPA)OOH]+ or[Cu(IDB)OOH]+ ligand can stably lie in the minor groove of DNA. There are some contacts between the linker and the minor groove of DNA, which leads to the increase of binding affinity, resulting in enhancing the probability of DNA strand scission. These observations suggest that a flexibility linker between the recognition and cleavage agent can provide crucial contributions to specific DNA cleavage of metal nucleases. In addition, it can be seen from the comparison of the HPD-BPA(IDB) and CPDγTail-BPA(IDB) systems, the two nucleases assisted by both hairpin polyamide and modified cyclic polyamide with a linker -CH2CH2CH2NH2 all can achieve selective DNA cleavage via the hydrogen abstraction from the sugar moiety of DNA, although the abstracted hydrogen type is different, i.e., the nucleases bound to the hairpin polyamide can efficiently abstract the C4'H atom of sugar; however, the nucleases bound to the modified cyclic polyamide with a linker -CH2CH2CH2NH2 can efficiently abstract the C1'H atom of sugar. That is to say the chemical nucleases bound to an appropriate polyamide-type recognition agent can improve efficiently the DNA cleavage selectivity.
A series of molecular dynamics simulations was performed to investigate the DNA cleavage specificities of the nucleases [Cu(BPA)]2+ and [Cu(IDB)]2+ assisted by the DNA recognition agents, hairpin and cyclic polyamides, via the analysis of hydrogen atom abstractions from the sugar rings in DNA. The simulated results indicate that the studied polyamide bound by the nuclease [Cu(BPA)]2+ or [Cu(IDB)]2+, presents a stable conformation located in the minor groove of DNA, which is consistent with the conformation of the corresponding polyamide bound to DNA presented in X-ray structure. For the different nucleases assisted by hairpin polyamide in the two systems of hairpin-polyamide + [Cu(BPA)OOH] + - DNA and hairpin-polyamide + [Cu(IDB)OOH] + - DNA, the nucleases present the similar DNA cleavage properties. That is to say, the [Cu(BPA) ]2+ or [Cu(IDB) ]2+ nuclease with a substrate OOH- bound to the hairpin polyamide can be stably located at the minor groove of DNA, and possibly abstracts the C4'H atom from the sugar moiety of DNA. However, for the [Cu(BPA)]2+ or [Cu(BPA)]2+ nuclease assisted with the cyclic polyamide in the systems of cyclic-polyamide + [Cu(BPA)OOH] + - DNA and cyclic-polyamide + [Cu(IDB)OOH] + - DNA, the probabilities of any hydrogen abstraction of the two nucleases from the sugar moiety of DNA are very small, which suggests that the nucleases cannot achieve DNA strand scission, due to the short and rigid linker between the branch C atom of nucleases and backbone of the cyclic polyamide. Therefore, the flexible linker of -CH2CH2CH2NH2 of γ region was introduced to the cyclic polyamide to improve DNA cleavage selectivity of the nucleases. It is of note that with the flexible linker introduction, the simulated results for the systems of modified-cyclic-polyamide + [Cu(BPA)OOH] + - DNA and modified-cyclic-polyamide + [Cu(IDB)OOH] + - DNA exhibit the high possibility of DNA cleavage by the [Cu(BPA)]2+ or [Cu(IDB)]2+ nuclease bound at the tail of the linker of the modified cyclic polyamide via C1'H atom abstraction from the sugar ring in DNA. Current observations suggest that the polyamide-type recognition agents with a flexible linker (γ-amino butyric acid) can efficiently improve the DNA cleavage properties of copper nucleases.
The authors acknowledge research support from the National Science Foundation of China (No. 20631020, 20771017, 20973024 and 21073015), the Major State Basic Research Development Programs (grant No. 2011CB808500). We also thank the HPSC of Beijing Normal University for providing computer resources.
College of Chemistry, Beijing Normal University
Geierstanger BH, Mrksich M, Dervan PB, Wemmer DE: Design of a G·C-specific DNA minor groove-binding peptide.Science 1994, 266: 646–650. 10.1126/science.7939719View ArticlePubMed
Pratviel G, Bernadou J, Meuminer B: Carbon-hydrogen bonds of DNA sugar units as targets for chemical nucleases and drugs.Angew Chem Int Ed Engl 1995, 34: 746–769. 10.1002/anie.199507461View Article
Jiang Q, Xiao N, Shi PF, Zhu YG, Guo ZJ: Design of artificial metallonucleases with oxidative mechanism.Coord Chem Rev 2007, 251: 1951–1972. 10.1016/j.ccr.2007.02.013View Article
Sigman DS, Bruice TW, Mazumder A, Sutton CL: Targeted chemical nucleases.Acc Chem Res 1993, 26: 98–104. 10.1021/ar00027a004View Article
Sigman DS, Mazumder A, Perrin DM: Chemical nucleases.Chem Rev 1993, 93: 2295–2316. 10.1021/cr00022a011View Article
Erkkila KE, Odom DT, Barton JK: Recognition and reaction of metallointercalators with DNA.Chem Rev 1999, 99: 2777–2795. 10.1021/cr9804341View ArticlePubMed
Liu CL, Wang M, Zhang TL, Sun HZ: DNA hydrolysis promoted by di- and multi-nuclear metal complexes.Coord Chem Rev 2004, 248: 147–168. 10.1016/j.cct.2003.11.002View Article
McMilllin DR, McNett KM: Photoprocesses of copper complexes that bind to DNA.Chem Rev 1998, 98: 1201–1219. 10.1021/cr9601167View Article
Humphreys KJ, Johnson AE, Karlin KD, Rokita SE: Oxidative strand scission of nucleic acids by a multinuclear copper(II) complex.J Biol Inorg Chem 2002, 7: 835–842. 10.1007/s00775-002-0369-8View ArticlePubMed
Li L, Karlin KD, Rokita SE: Changing selectivity of DNA oxidation from deoxyribose to guanine by ligand design and a new binuclear copper complex.J Am Chem Soc 2005, 127: 520–521. 10.1021/ja044209eView ArticlePubMed
Humphreys KJ, Karlin KD, Rokita SE: Targeted strand scission of DNA substrates by a tricopper(II) coordination complex.J Am Chem Soc 2002, 124: 8055–8066. 10.1021/ja012539iView ArticlePubMed
Humphreys KJ, Karlin KD, Rokita SE: Efficient and specific strand scission of DNA by a dinuclear copper complex: comparative reactivity of complexes with linked tris(2-pyridylmethyl)amine moieties.J Am Chem Soc 2002, 124: 6009–6019. 10.1021/ja020039zView ArticlePubMed
Humphreys KJ, Karlin KD, Rokita SE: Recognition and strand scission at junctions between single- and double-stranded DNA by a trinuclear copper complex.J Am Chem Soc 2001, 123: 5588–5589. 10.1021/ja010403eView ArticlePubMed
Li L, Murthy NN, Telser J, Zakharov LN, Yap GPA, Rheingold AL, Karlin KD, Rokita SE: Targeted guanine oxidation by a dinuclear copper(II) complex at single stranded/double stranded DNA junctions.Inorg Chem 2006, 45: 7144–7159. 10.1021/ic0605930View ArticlePubMed
Zhu Q, Lian Y, Thyagarajan S, Rokita SE, Karlin KD, Blough NV: Hydrogen peroxide and dioxygen activation by dinuclear copper complexes in aqueous solution: hydroxyl radical production initiated by internal electron transfer.J Am Chem Soc 2008, 130: 6304–6305. 10.1021/ja800080zView ArticlePubMed
Sigman DS, Graham DR, D'Aurora V, Stern AM: Oxygen-dependent cleavage of DNA by the 1,10-phenanthroline-cuprous complex. Inhibition of E. coli DNA polymerase I.J Biol Chem 1979, 254: 12269–12272.PubMed
Graham DR, Marshall LE, Reich KA, Sigman DS: Cleavage of DNA by coordination complexes. superoxide formation in the oxidation of 1,10-phenanthroline-cuprous complexes by oxygen - relevance to DNA-cleavage reaction.J Am Chem Soc 1980, 102: 5419–5421. 10.1021/ja00536a063View Article
Pratviel G, Duarte V, Bernadou J, Meunier B: Nonenzymic cleavage and ligation of DNA at a three A·T base pair site. A two-step pseudohydrolysis of DNA.J Am Chem Soc 1993, 115: 7939–7943. 10.1021/ja00071a002View Article
Keck MV, Lippard SJ: Unwinding of supercoiled DNA by platinum-ethidium and related complexes.J Am Chem Soc 1992, 114: 3386–3390. 10.1021/ja00035a033View Article
An Y, Liu SD, Deng SY, Ji LN, Mao ZW: Cleavage of double-strand DNA by linear and triangular trinuclear copper complexes.J Inorg Biochem 2006, 100: 1586–1593. 10.1016/j.jinorgbio.2006.05.002View ArticlePubMed
Thyagarajan S, Murthy NN, Narducci Sarjeant AA, Karlin KD, Rokita SE: Selective DNA strand scission with binuclear copper complexes: implications for an active Cu2-O2 species.J Am Chem Soc 2006, 128: 7003–7008. 10.1021/ja061014tPubMed CentralView ArticlePubMed
Dervan PB: Reagents for the site-specific cleavage of megabase DNA.Nature 1992, 359: 87–88. 10.1038/359087a0View ArticlePubMed
Sreedhara A, Freed JD, Cowan JA: Efficient inorganic deoxyribonucleases greater than 50-million-fold rate enhancement in enzyme-like DNA cleavage.J Am Chem Soc 2000, 122: 8814–8824. 10.1021/ja994411vView Article
Pogozelski WK, Tullius TD: Oxidative strand scission of nucleic acids: routes initiated by hydrogen abstraction from the sugar moiety.Chem Rev 1998, 98: 1089–1108. 10.1021/cr960437iView ArticlePubMed
Frank BL, Leroy W Jr, Christner DF, Kozarich JW, Stubbe J, Kappen LS, Goldberg IH: J Am Chem Soc. 1991, 113: 2271–2275. 10.1021/ja00006a054View Article
Liang Q, Ananias DC, Long EC: Ni(II)·Xaa-Xaa-His induced DNA cleavage: deoxyribose modification by a common "activated" intermediate derived from KHSO5, MMPP, or H2O2.J Am Chem Soc 1998, 120: 248–257. 10.1021/ja9720218View Article
Mirica LM, Ottenwaelder X, Stack TDP: Structure and spectroscopy of copper-dioxygen complexes.Chem Rev 2004, 104: 1013–1045. 10.1021/cr020632zView ArticlePubMed
Zhu YY, Wang Y, Chen GJ, Zhan CG: A three-point method for evaluations of AMBER force field parameters: an application to copper-based artificial nucleases.Theor Chem Account 2009, 122: 167–178. 10.1007/s00214-008-0496-6View Article
Zhu YY, Su YW, Qi Y, Tan HW, Wang Y, Chen GJ: Molecular docking studies on the metallonucleases and polyamides with double-strand DNA molecules.Chem J Chinese Univ 2009, 30: 781–785.
Zhu YY, Chen ZF, Guo ZJ, Wang Y, Chen GJ: Systematic characterization on electronic structures and spectra for a series of complexes, M(IDB)Cl2 (M = Mn, Fe, Co, Ni, Cu and Zn): a theoretical study.J Mol Model 2009, 15: 469–479. 10.1007/s00894-008-0432-7View ArticlePubMed
Zhu YY, Su YW, Li XC, Wang Y, Chen GJ: Evaluation of AMBER force field parameters for copper(II) with pyridylmethyl-amine and benzimidazolylmethyl-amine ligands: a quantum chemical study.Chem Phys Letter 2008, 455: 354–360. 10.1016/j.cplett.2008.03.004View Article
Jin Y, Cowan JA: DNA cleavage by copper-ATCUN complexes. factors influencing cleavage mechanism and linearization of dsDNA.J Am Chem Soc 2005, 127: 8408–8415. 10.1021/ja0503985View ArticlePubMed
Meng X, Liu L, Zhou C, Wang L, Liu C: Dinuclear copper(II) complexes of a polybenzimidazole ligand: their structures and inductive roles in DNA condensation.Inorg Chem 2008, 47: 6572–6574. 10.1021/ic800532qView ArticlePubMed
Chen TQ, Greenberg MM: Model studies indicate that copper phenanthroline induces direct strand breaks via β-elimination of the 2'-deoxyribonolactone intermediate observed in enediyne mediated DNA damage.J Am Chem Soc 1998, 120: 3815–3816. 10.1021/ja9800075View Article
Pitie M, Burrows CJ, Meunier B: Mechanisms of DNA cleavage by copper complexes of 3-Clip-Phen and of its conjugate with a distamycin analogue.Nucleic Acids Res 2000, 28: 4856–4864. 10.1093/nar/28.24.4856PubMed CentralView ArticlePubMed
Uma V, Kanthimathi M, Subramanian J, Unni Nair B: A new dinuclear biphenylene bridged copper(II) complex: DNA cleavage under hydrolytic conditions.Biochim Biophy Acta 2006, 1760: 814–819.View Article
White S, Szewczyk JW, Turner JM, Baird EE, Dervan PB: Recognition of the four Watson-Crick base pairs in the DNA minor groove by synthetic ligands.Nature 1998, 391: 468–471. 10.1038/36116View ArticlePubMed
Dervan PB, Edelson BS: Recognition of the DNA minor groove by pyrrole-imidazole polyamides.Curr Opin Struct Biol 2003, 13: 284–299. 10.1016/S0959-440X(03)00081-2View ArticlePubMed
Ellervik U, Wang CCC, Dervan PB: Hydroxybenzamide/pyrrole pair distinguishes T·A from A·T base pairs in the minor groove of DNA.J Am Chem Soc 2000, 122: 9354–9360. 10.1021/ja001692uView Article
Jaramillo D, Wheate NJ, Ralph SF, Howard WA, Tor Y, Aldrich-Wright JR: Polyamide platinum anticancer complexes designed to target specific DNA sequences.Inorg Chem 2006, 45: 6004–6013. 10.1021/ic060383nView ArticlePubMed
Kielkopf CL, Bremer RE, White S, Szewczyk JW, Turner JM, Baird EE, Dervan PB, Rees DC: Structural effects of DNA sequence on T·A recognition by hydroxypyrrole/pyrrole pairs in the minor groove.J Mol Biol 2000, 295: 557–567. 10.1006/jmbi.1999.3364View ArticlePubMed
Parks ME, Baird EE, Dervan PB: Optimization of the hairpin polyamide design for recognition of the minor groove of DNA.J Am Chem Soc 1996, 118: 6147–6152. 10.1021/ja960721rView Article
Swalley SE, Baird EE, Dervan PB: Discrimination of 5'-GGGG-3', 5'-GCGC-3', and 5'-GGCC-3' sequences in the minor groove of DNA by eight-ring hairpin polyamides.J Am Chem Soc 1997, 119: 6953–6961. 10.1021/ja971184tView Article
Chenoweth DM, Harki DA, Phillips JW, Dose C, Dervan PB: Cyclic pyrrole-imidazole polyamides targeted to the androgen response element.J Am Chem Soc 2009, 131: 7182–7188. 10.1021/ja901309zPubMed CentralView ArticlePubMed
Wemmer DE: Designed sequence-specific minor groove ligands.Annu Rev Biophys Biomol Struct 2000, 29: 439–461. 10.1146/annurev.biophys.29.1.439View ArticlePubMed
Mrksich M, Parks ME, Dervan PB: Hairpin peptide motif. a new class of oligopeptides for sequence-specific recognition in the minor groove of double-helical DNA.J Am Chem Soc 1994, 116: 7983–7988. 10.1021/ja00097a004View Article
Zhao YM, Zhu JH, He WJ, Yang Z, Zhu YG, Li YZ, Zhang JF, Guo ZJ: Oxidative DNA cleavage promoted by multinuclear copper complexes: activity dependence on the complex structure.Chem Eur J 2006, 12: 6621–6629. 10.1002/chem.200600044View ArticlePubMed
Suto RK, Edayathumangalam RS, White CL, Melander C, Gottesfeld JM, Dervan PB, Luger K: Crystal structures of nucleosome core particles in complex with minor groove DNA-binding ligands.J Mol Biol 2003, 326: 371–380. 10.1016/S0022-2836(02)01407-9View ArticlePubMed
Zhang Q, Dwyer TJ, Tsui V, Case DA, Cho J, Dervan PB, Wemmer DE: NMR structure of a cyclic polyamide-DNA complex.J Am Chem Soc 2004, 126: 7958–7966. 10.1021/ja0373622View ArticlePubMed
Case DA, Darden TA, Cheatham TE, Simmerling CL, Wang J, Duke RE, Luo RM, K M, Pearlman DA, Crowley M, Walker RC, et al.: AMBER9. University of California, San Francisco; 2006.
Herman DM, Turner JM, Baird EE, Dervan PB: Cycle polyamide motif for recognition of the minor groove of DNA.J Am Chem Soc 1999, 121: 1121–1129. 10.1021/ja983206xView Article
Zhu YY, Wang Y, Chen GJ: Molecular dynamics simulations on binding models of Dervan-type polyamide + Cu(II) nuclease ligands to DNA.J Phys Chem B 2009, 113: 839–848. 10.1021/jp8091545View ArticlePubMed
Young MA, Jayaram B, Beveridge DL: Intrusion of counterions into the spine of hydration in the minor groove of B-DNA: fractional occupancy of electronegative pockets.J Am Chem Soc 1997, 119: 59–69. 10.1021/ja960459mView Article
Young MA, Ravishanker G, Beveridge DL: A 5-nanosecond molecular dynamics trajectory for B-DNA: analysis of structure, motions, and solvation.Biophys J 1997, 73: 2313–2336. 10.1016/S0006-3495(97)78263-8PubMed CentralView ArticlePubMed
De Souza ON, Ornstein RL: Effect of warmup protocol and sampling time on convergence of molecular dynamics simulations of a DNA dodecamer using AMBER 4.1 and particle-mesh Ewald method.J Biomol Struct Dyn 1997, 14: 607–611.View Article
Duan Y, Wu C, Chowdhury S, Lee MC, Xiong G, Zhang W, Yang R, Cieplak P, Luo R, Lee T, et al.: A point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum mechanical calculations.J Comput Chem 2003, 24: 1999–2012. 10.1002/jcc.10349View ArticlePubMed
Lee MC, Duan Y: Distinguish protein decoys by using a scoring function based on a new AMBER force field, short molecular dynamics simulations, and the generalized Born solvent model.Proteins: Structure, Function, and Bioinformatics 2004, 55: 620–634. 10.1002/prot.10470View Article
Wang JM, Wolf RM, Caldwell JW, Kollman PA, Case DA: Development and testing of a general Amber force field.J Comput Chem 2004, 25: 1157–1174. 10.1002/jcc.20035View ArticlePubMed
Roothaan CCJ: New developments in molecular orbital theory.Rev Modern Phys 1951, 23: 69–89. 10.1103/RevModPhys.23.69View Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.