- Research article
- Open Access
Disturbance of DNA conformation by the binding of testosterone-based platinum drugs via groove-face and intercalative interactions: a molecular dynamics simulation study
BMC Structural Biology volume 13, Article number: 4 (2013)
To explore novel platinum-based anticancer agents that are distinct from the structure and interaction mode of the traditional cisplatin by forming the bifunctional intrastrand 1,2 GpG adduct, the monofunctional platinum + DNA adducts with extensive non-covalent interactions had been studied. It was reported that the monofunctional testosterone-based platinum(II) agents present the high anticancer activity. Moreover, it was also found that the testosterone-based platinum agents could cause the DNA helix to undergo significant unwinding and bending over the non-testosterone-based platinum agents. However, the interaction mechanisms of these platinum agents with DNA at the atomic level are not yet clear so far.
In the present work, we used molecular dynamics (MD) simulations and DNA conformational dynamics calculations to study the DNA distortion properties of the testosterone-based platinum + DNA, the improved testosterone-based platinum + DNA and the non-testosterone-based platinum + DNA adducts. The results show that the intercalative interaction of the improved flexible testosterone-based platinum agent with DNA molecule could cause larger DNA conformational distortion than the groove-face interaction of the rigid testosterone-based platinum agent with DNA molecule. Further investigations for the non-testosterone-based platinum agent reveal the occurrence of insignificant change of DNA conformation due to the absence of testosterone ligand in such agent. Based on the DNA dynamics analysis, the DNA base motions relating to DNA groove parameter changes and hydrogen bond destruction of DNA base pairs were also discussed in this work.
The flexible linker in the improved testosterone-based platinum agent causes an intercalative interaction with DNA in the improved testosterone-based platinum + DNA adduct, which is different from the groove-face interaction caused by a rigid linker in the testosterone-based platinum agent. The present investigations provide useful information of DNA conformation affected by a testosterone-based platinum complex at the atomic level.
Since the discovery of cisplatin [cis-(NH3)2PtCl2], as an anticancer agent, a series of platinum anticancer drugs were used in treatment of various cancers in clinical chemotherapy [1–6]. The platinum(II) center of traditional cisplatin reacts with DNA forming two covalent bonds to N7 atoms of two adjacent guanine (G) bases with the bifunctional intrastrand 1,2 GpG adduct, to prevent the replication and transcription of cancer genes, and ultimately to induce tumor cell apoptosis [6–9]. However, the efficacy of traditional cisplatin is often compromised because of the intrinsic and acquired resistance, as well as the toxic side effects [3, 10–14]. Much effort has been devoted to the developments of novel platinum-based anticancer agents which might form monofunctional platinum + DNA adducts with extensive non-covalent interactions to circumvent such drawbacks [11, 15, 16]. The structures of these monofunctional platinum + DNA adducts with extensive non-covalent interactions are different from that produced by cisplatin [8, 9]. To the best of our knowledge, the non-covalent interactions in platinum + DNA adducts, mostly resulting from the interactions between the ligands of platinum agents and DNA molecules, include the electrostatic interaction, groove-face interaction and intercalative interaction [17–22]. Especially, the groove-face and intercalative interactions could greatly stabilize the distorted DNA molecule via hydrogen bonds and hydrophobic interactions, etc [17, 18, 20, 21]. Therefore, some pioneering strategies toward improving the ligand properties of monofunctional platinum agents have emerged [15, 23–27]. However, the systematic studies on the relationship between the ligand properties of platinum agents and interaction modes in the platinum + DNA adducts have not yet been clearly detailed so far.
Recently, it has been shown that the monofunctional platinum complexes of nitrogen-containing heterocyclic amines, such as pyridine, pyrimidine, purine, piperidine, picoline, and their derivatives could enhance cytotoxicity, perhaps because of formations of monofunctional platinum + DNA adducts rather than those of bifunctional ones [15, 23]. Especially, Hannon and co-workers reported that the linkage of a testosterone to aromatic N-heteroatomic monofunctional platinum(II) agents confers relatively high activity to otherwise non-active platinum(II) agents. Moreover, they indicated that the conjugation of testosterone enhances the delivery ability into the tumor cell, and inherent antitumor activity of testosterone-based platinum agents . Importantly, it was also proved that the testosterone-based platinum agents cause the DNA helix to undergo significant unwinding and bending over the non-testosterone-based platinum agents. This might be caused by the steric bulk of testosterone which requires greater unwinding/bending of DNA helix to accommodate the agent in the DNA double helix . Theoretically, the structures of some new designed drugs were usually optimized by the density functional theory (DFT) method . Some molecular dynamics (MD) simulations were used for investigating the interaction properties between DNA molecule and different platinum agents, such as cisplatin, oxaliplatin, [PtCl(en)(ACRAMTU-S)](NO3)2 and so on [29–36]. Moreover, the MD simulations were also used for investigating the molecular interactions between the androgen receptor protein and testosterones in the platinum agents [24, 37]. However, few theoretical studies devoted to the interactions between the monofunctional N-heteroatomic platinum agents with the testosterone ligands and DNA molecules.
Delineating the structural details of DNA bound by platinum agents will help us to understand the features that are responsible for the remarkable potency of these platinum agents. Nevertheless, how the testosterone ligands affect the interaction modes of platinum agents with DNA molecules is still unknown. In the present work, we used molecular dynamics simulations and DNA dynamics analysis to study the conformational properties of DNA disturbed by different monofunctional N-heteroatomic platinum agents. Due to the potent antitumor activity of cis-[Pt(Testo)Cl][NO3] ( Testo = (NH3)2(17α-pyridyl-3-ethynyltestosterone)) as a rigid testosterone-based platinum agent studied by Hannon and co-workers , the first MD simulation was performed on the Pt(Testo)(II) + DNA adduct to investigate interaction properties of this platinum agent with DNA, and effects of the testosterone ligand on DNA conformation. To better understand the effects of the testosterone ligand on interaction modes in the platinum + DNA adducts, the second MD simulation was performed on the improved Pt(Testo)(II) agent with flexible linker interacting with DNA molecule. To compare with these two testosterone-based platinum agents, a non-testosterone-based platinum agent adding to DNA molecule has also been simulated.
On the basis of previous experimental studies [26, 27], we choose the initial structure of the cis-[ Pt(Testo)Cl][NO3] ( Testo = (NH3)2(17α-pyridyl-3-ethynyltestosterone)) from the geometry optimized by the DFT method at the B3LYP/6-31G** + LanL2DZ level of theory due to the inexistence of its crystal structure (see Figure 1(a)). The corresponding geometry parameters with the available experimental data from similar molecules are shown in Additional file 1: Table S1 and Additional file 2: Figure S1. The 20bp DNA fragment, 5′-d(CGGTGAAAACCTCTGACACA)-3′, was chosen from the primer PUC19 F studied by the experiment of polymerase-chain reaction (PCR)  and shown in Figure 1(b). To build a Pt(Testo)(II) + DNA adduct, its initial coordinates used in our simulations were generated by manually docking [38–40]. Namely, the platinum center is bound to the N7 atom of G15 base in the major groove of DNA due to the experimental prediction for a similar platinum agent . The plane of platinum center coordinated by four ligands is perpendicular to the plane of G15 base; ultimately, the linked pyridyl ring is approximately perpendicular to the base pair planes of DNA, and is parallel to the wall of major groove of DNA; the testosterone ligand with four 5/6-membered rings linked to the pyridyl ring is approximately localized around the T14:A27 base pair across the major groove of DNA, minimising the steric overlap between the agent and the DNA molecule. Since several different orientations of platinum center plane to the G15 base of DNA for the Pt(Testo)(II) + DNA adduct were chosen as the starting structures for our MD simulations, this provides a good test of whether our MD simulations are capable of driving significantly distinct starting structures to a non-distinguishable one when simulations reach equilibrium. To further investigate the influence of testosterone ligand of platinum agent on interaction modes, we improved the flexibility of linker between the pyridyl and testosterone ligands in the Pt(Testo)(II) agent with changing the linker of ethynyl –C ≡ C– to that of ethyl –CH2–CH2– (assigned as Im-Pt(Testo)(II) agent); simultaneously a non-testosterone-based platinum agent, cis-[Pt(NH3)2(pyridine)Cl][NO3] (assigned as Pt(Py)(II) agent), was also chosen to perform a comparison. The Im-Pt(Testo)(II) and Pt(Py)(II) agents were fully optimized by the same method used for the Pt(Testo)(II) agent, and shown in Figure 1(c) and (d) [26, 27]. The Im-Pt(Testo)(II) + DNA and Pt(Py) (II) + DNA adducts were similarly built by using the same method for the Pt(Testo)(II) + DNA adduct described above. To compare the differences between the DNA adduct conformations and an undamaged DNA molecule, a bare-DNA (B-DNA) molecule simulation was also performed, in which an idealized B-DNA with the same base-pair sequence in these adducts was used as a starting structure for the simulation. Given that each strand of DNA has some phosphate groups, Na+ counterions were added to each system to achieve electroneutrality. The systems were explicitly solvated using the TIP3P water potential inside a box large enough to ensure the solvent shell extended to 10 Å in all directions.
Force field parameter preparation
The atom types for the studied platinum agents, except for the platinum atoms, were generated using the ANTECHAMBER module in the AMBER9 program . The electrostatic potentials of the platinum agents used for RESP charge calculations were calculated at the B3LYP/6-31G** + LanL2DZ [28, 31, 42, 43] level of theory using the Gaussian09 program . The RESP charges of platinum agents were derived by the RESP program based on the calculated electrostatic potentials. The force field parameters around the Pt center were generated by quantum chemical calculations, which was reported particularly in our previous work [45, 46]. Other force field parameters of the platinum agents were generated from the gaff force field in the AMBER9 program.
Molecular dynamics simulations protocols
All MD simulations were carried out using the AMBER9 package  with the parm99 force field [47, 48], the parmbsc0 refinement  and gaff  force field parameters. The details of MD protocols are given in Additional file 3: Methodology.
Principal component analysis
Principal component analysis can be used to segregate large-scale correlated motions from random thermal fluctuations, thereby probing the essential dynamics of the system. The details of this analysis method are available in Additional file 3: Methodology.
DNA groove parameter analyses
The frequency distributions (fraction of the time spent in each conformation) from the trajectories of simulations for the models and a canonical B-DNA were calculated using the CURVES program  to investigate the distortion of DNA. To account for the distortion of whole DNA backbone, the overall bend, tilt and roll angles of the DNA time-averaged structures for the studied models were calculated by using the MadBend program from the CURVES outputs . The details of the calculation method are available in Additional file 3: Methodology.
The root-mean-square deviation (RMSD) values of all backbone atoms referenced to the corresponding starting structures over all three trajectories for the Pt(Testo)(II) + DNA, Im-Pt(Testo)(II) + DNA and Pt(Py)(II) + DNA models were examined to determine if each system had attained equilibrium. It is often considered that small RMSD values of a simulation indicate a stable state of the system. Plots of RMSDs of three system simulations over times are shown in Figure 2. It can be seen from Figure 2 that each platinum + DNA adduct reached equilibrium after 30 ns, and their energies were found to be stable during the remainder of each simulation. Therefore, the trajectory analysis for each of three systems has extracted the equilibrated conformation between 30 ns and 50 ns of simulation time, recording 10000 snapshots at every 2 ps time-interval of each trajectory.
Disturbance of DNA conformation by the binding of Pt(Testo)(II) agent via groove-face interaction
Conformation analysis of Pt(Testo)(II) + DNA adduct
Based on the results of previous experiment , the testosterone-based platinum agent, cis-[Pt(NH3)2(17α-pyridyl-3-ethynyltestosterone)Cl]+ (called Pt(Testo)(II)), which consists of testosterone, ethynyl linker and platinum(II) center coordinated by one chlorine atom and three nitrogen atoms from two ammonias and one pyridyl, as an efficient anticancer agent, binding to a DNA molecule (called Pt(Testo)(II) + DNA adduct), was tested by MD simulation. The average structure of Pt(Testo)(II) + DNA adduct was obtained by the trajectory analysis extracting the equilibrated conformations between 30 ns and 50 ns of simulation times and shown in Figure 3 along with that of undamaged DNA molecule. It is demonstrated from this average structure that the testosterone with a groove-face interaction mode locates at the major groove of DNA molecule; the platinum center binds to the N7 atom of G15 base of DNA molecule; the pyridyl ring is perpendicular to the plane of G15 base of DNA; the plane of testosterone with the rigid ethynyl (–C ≡ C–) linker is perpendicular to the pyridyl ring and plane of T14:A27 base pair of DNA molecule.
Because the testosterone ligand locates at the DNA major groove via the groove-face interaction mode, both major groove and minor groove of DNA molecule tend to widen and shoal. Figure 4 shows the DNA groove parameters around the connection sites, C11:G30 ~ G15:C26 base pairs of DNA, which are near to the testosterone position, for the Pt(Testo)(II) + DNA and B-DNA models. Namely, the average major groove width and minor groove width increase from 12.13 to 13.68 Å and from 5.75 to 8.04 Å, respectively, compared to a normal B-DNA; then their depths are shoaled from 6.00 to 2.04 Å and from 4.73 to 3.63 Å, respectively (Figure 4(a), (b), (c), (d) and (e); see the red lines with circles and the black lines with squares for the Pt(Testo)(II) + DNA and B-DNA models, respectively). In addition, the measured degrees of overall DNA bend, tilt and roll angles for the bound-DNA backbone conformation of Pt(Testo)(II) + DNA adduct are shown in Table 1. It can be seen that the average deviation percentage of helical angles with respect to the undamaged B-DNA molecule is 261.64% (260.02% for bend; 270.54% for tilt; 254.37% for roll) for the Pt(Testo)(II) + DNA adduct. Interestingly, the bend angle toward the major groove was observed, which suggests there is great stability of this adduct for this orientation with the groove-face interaction mode. As expected, the testosterone-based platinum agent could cause a great DNA helical conformation distortion, which is consistent with the previous experimental study .
Hydrogen bond and hydrophobic interaction analyses of Pt(Testo)(II) + DNA adduct
To address the extent of DNA distortion through the groove-face interaction mode, the destruction of hydrogen bonds in the DNA molecule and the formation of new hydrogen bonds between the Pt agent and DNA molecule have been analyzed and shown in Figure 5, Figure 6(a) and Table 2, respectively. A distance of less than 3.5 Å and an angle of greater than 120° between the potential hydrogen bond donor and acceptor were used as the criteria for a hydrogen bond formation [33, 53]. Due to the interaction of testosterone ligand with DNA bases at the major groove, the hydrogen bonds at T14:A27 base pair in DNA molecule are completely destroyed compared with the undamaged B-DNA (see Figure 5). The destructions of hydrogen bonds result from the formations of some new hydrogen bonds between the C-H groups of testosterone ligand and O6 atom of G28 base, and between the C-H groups of pyridyl ligand and O6 atom of G15 base with 20.82% and 13.78% occupancy times, respectively (see Table 2 and Figure 6(a)). Moreover, hydrophobic interactions, defined as a distance between carbon atoms shorter than 4.5 Å, have been analyzed and shown in Table 3. Namely, it was found that the hydrophobic interactions between the carbon atoms of testosterone ligand and carbon atoms of C26, A27 and T14 bases of DNA occur with the occupancy times of 211.69%, 36.08% and 23.44%, respectively. Both hydrogen bond destruction in DNA molecule and new interaction formation of Pt(Testo)(II)-DNA interface facilitate the DNA distortion, which might make DNA repair more difficult.
Disturbance of DNA conformation by the binding of Im-Pt(Testo)(II) agent via intercalative interaction
It has been shown from the structure analysis of Pt(Testo)(II) + DNA adduct discussed above that the Pt(Testo)(II) agent interacts with DNA molecule adopting the groove-face interaction mode due to the rigid ethynyl –C ≡ C– linker between the testosterone and pyridyl ligands. To produce great effects of Pt(Testo)(II) agent on DNA conformation distortion via intercalative interaction mode, the rigid ethynyl –C ≡ C– linker in the Pt(Testo)(II) agent was changed to a flexible linker, ethyl –CH2–CH2–, in order to build a new modified platinum agent, called Im-Pt(Testo)(II). The Im-Pt(Testo)(II) agent, as an improved anticancer agent, binding to a DNA molecule (called Im-Pt(Testo)(II) + DNA adduct) was also tested by MD simulation. The average structure of Im-Pt(Testo)(II) + DNA adduct are shown in Figure 3 along with that of undamaged DNA molecule, and demonstrates that the plane of testosterone ligand embeds into the middle of T12 and C13 bases of DNA molecule at the DNA major groove via an intercalative interaction mode, and parallels to both T12 and C13 bases; the platinum center also binds to the N7 atom of G15 base of DNA molecule; the pyridyl ring is still perpendicular to the plane of G15 base of DNA. The DNA conformation distortion greatly increases due to the modified platinum agent interacting with DNA molecule via the intercalative interaction mode.
Because the testosterone ligand embeds into the bases of DNA molecule at the major groove via the intercalative interaction mode, the major and minor grooves of DNA molecule in the Im-Pt(Testo)(II) + DNA adduct become wider and shallower than those in the Pt(Testo)(II) + DNA adduct. Namely, the average width of DNA major groove in the Im-Pt(Testo)(II) + DNA adduct increases by 74.27% (from 13.68 Å to 23.84 Å) compared to the Pt(Testo)(II) + DNA adduct; then the average major groove depth and minor groove depth are shoaled by 14.71% (from 2.04 Å to 1.74 Å) and by 18.73% (from 3.63 Å to 2.95 Å), respectively (Figure 4(b), (c) and (e); see the red lines with circles and the Magenta lines with down-triangles for Pt(Testo)(II) + DNA and Im-Pt(Testo)(II) + DNA adducts, respectively). In addition, the average deviation percentage of helical angles for the Im-Pt(Testo)(II) + DNA adduct is larger by 43.63% than that for the Pt(Testo)(II) + DNA adduct (see Table 1). Similarly, the occupancy percentages of formations of new hydrogen bonds and hydrophobic interactions in the Im-Pt(Testo)(II) + DNA adduct are larger than those in the Pt(Testo)(II) + DNA adduct due to the intercalative interaction mode, except for the same extent destruction of hydrogen bonds at the T12:A29 base pair of DNA molecule compared with the Pt(Testo)(II) + DNA adduct. Namely, the Im-Pt(Testo)(II) + DNA adduct spends more occupancy times by 38.14% (from 34.60% to 72.74%) forming new hydrogen bonds between the C-H groups of testosterone ligand and the O2/N3/N3 atom of T12/A29/G30 base than the Pt(Testo)(II) + DNA adduct (shown in Figure 6(b)); simultaneously, it still spends more occupancy times by 178.5% (from 271.21% to 449.71%) forming more hydrophobic interactions between C atoms of testosterone and C atoms of T12/C13/A29/G30 base, and between C atoms of ethyl linker and C atoms of T14 base (shown in Table 3). These results predict that the distortion of DNA conformation in the Im-Pt(Testo)(II) + DNA adduct via the intercalative interaction mode is greater than that in the Pt(Testo)(II) + DNA adduct via the groove-face interaction mode.
Binding of a non-testosterone-based platinum agent to DNA
Based on the previous experiments [26, 27], a non-testosterone-based platinum complex displays little anticancer activity. We also performed MD simulation on the non-testosterone-based platinum agent, cis-[Pt(NH3)2(pyridine)Cl]+ (assigned as Pt(Py)(II)), in which a platinum(II) center is coordinated by one chlorine atom and three nitrogen atoms from two ammonias and only one pyridine, binding to DNA molecule (assigned as Pt(Py)(II) + DNA adduct). It can be found that there is hardly any change in DNA conformation and DNA groove parameters in the Pt(Py)(II) + DNA adduct compared with the undamaged B-DNA (Figure 4(b), (c), (d) and (e); see the blue lines with up-triangles and the black lines with squares for the Pt(Py)(II) + DNA and B-DNA models, respectively), except that the average deviation percentage of helical angle is 143.93% with respect to the undamaged B-DNA molecule (see Table 1). Simultaneously, there are no any destruction of hydrogen bond in the DNA molecule and formation of new hydrogen bond/hydrophobic interaction between the Pt(Py)(II) agent and DNA molecule in the Pt(Py)(II) + DNA adduct due to no testosterone ligand in the Pt(Py)(II) agent. These results suggest that the DNA conformation changes caused by the binding of the non-testosterone-based platinum agent are unconspicuous, which predicts that the non-testosterone-based Pt(Py)(II) agent is an inefficient anticancer agent reported by the previous experiment .
Principal component analysis of major conformational dynamics
Principal component analysis was used to analyze the trajectories from the corresponding simulations to examine the dominant DNA dynamic motions in the studied Pt + DNA adducts. The first three principal components (PC1, PC2 and PC3) described about 80% of all the motion modes of dynamics for Pt(Testo)(II) + DNA, Im-Pt(Testo)(II) + DNA and Pt(Py)(II) + DNA adducts (see Table 4). It was found that the first three components of conformational motions roughly correspond to a superposition of bending, unwinding and twisting motions (see Additional file 4: Movie S1, Additional file 5: Movie S2 and Additional file 6: Movie S3, respectively, for the Im-Pt(Testo)(II) + DNA adduct). The overall average centroid structures of these adducts were analyzed and shown in Figure 3. Visual analyses of average DNA structures from the trajectories support the PCA dominant motions. It can be seen that the DNA conformations of Pt + DNA adducts with the groove-face and intercalative interaction modes including the bending, unwinding and twisting motions show obvious differences from the undamaged B-DNA conformation. For example, the percentages of occupancy times of PC1 and PC2 with the bending and unwinding motions, respectively, for both the Pt(Testo)(II) + DNA and Im-Pt(Testo)(II) + DNA adducts are larger than those for the Pt(Py)(II) + DNA adduct, which predicts that the platinum agents with the big-size testosterone ligand binding to DNA molecule via either groove-face interaction mode or intercalative interaction mode lead to greater bending and unwinding for DNA double helix than that without the testosterone ligand. Otherwise, the twisting motions for all three adducts make little contribution to the distortion of DNA conformation.
DNA conformational dynamics in the Pt + DNA adducts with different interaction modes
The Pt(Testo)(II) and Im-Pt(Testo)(II) agents respectively via the groove-face and intercalative interaction modes locate at the major groove of DNA and disturb the DNA conformation greatly. The frequency distributions of DNA helical dynamics parameters at the base pairs near the binding site and testosterone position have been analyzed and shown in Figure 7 and Additional file 7: Figure S2 for the Pt(Testo)(II) + DNA and Im-Pt(Testo)(II) + DNA adducts, respectively, along with the distribution patterns of the undamaged B-DNA. Generally, the groove width and depth were measured mainly by the opening, twist and rise motions, and the shift motion for DNA bases, respectively, which has been defined by the Curves analysis [54–56].
In the Pt(Testo)(II) + DNA model, the opening parameter of DNA base pair helix near the testosterone position T14:A27 base pair is changed by ~ -50° away from the B-DNA molecule (see Figure 7); the deviations of twist angle values are ~ -35º, ~25º, ~-25º and ~20º for the C13 · T14, T14 · G15, C26 · A27 and A27 · G28 base pair steps with respect to ~ 36º for a base pair step of B-DNA molecule, in which the value of 25º is agreement with the unwinding angle of 21º determined by using gel mobility shift assays for the Pt(Testo)(II) agent in plasmid DNA solution ; moreover, the deviations of rise values are ~ -1.5Å, ~5Å and ~4.5Å for the T12 · C13, C13 · T14 and C26 · A27 base pair steps, respectively, with respect to ~ 3.5 Å for the B-DNA molecule (see Figure 7), which results in a widening of major groove by 1.55 Å and minor groove by 2.29 Å around the C11:G30 ~ G15:C26 base pairs (see Figure 4(a), (b) (red line) and (d) (red line)). The deviations of shift values are ~2Å, ~1.5Å, ~-2.5Å and ~ -2Å for the T12 · C13, T14 · G15, C26 · A27and A27 · G28 base pair steps, respectively, away from the B-DNA molecule (see Figure 7), which results in a shoaling of major groove by 3.96 Å and minor groove by 1.10 Å around the C11:G30 ~ G15:C26 base pairs (see Figure 4(a), (c) (red line) and (e) (red line)). In addition, the deviations of roll angle values are ~15°, ~-25º, ~40º and ~ -20° for the T12 · C13, C13 · T14, C26 · A27 and A27 · G28 base pair steps, respectively, away from the B-DNA molecule (see Figure 7), which contributes to the significant bend of DNA helix toward the major groove by 30.50°. Moreover, the conformational dynamics of DNA helical parameters in the vicinity of the testosterone position are associated with the pattern of hydrogen bond destruction and formation. For example, except of the opening parameter mentioned above, the propel, stagger, and stretch parameters for the DNA molecule in the Pt(Testo)(II) + DNA model near the testosterone position T14:A27 base pair are also changed by ~ -20°, ~4Å and ~ -2.5Å, respectively, away from the B-DNA molecule, which results in the destruction of hydrogen bonds at the base pair T14:A27. In addition, the DNA base motions of shift, twist at T14 · G15 and A27 · G28 base pair steps, and roll at A27 · G28 base pair step with the large deviations mentioned above (see Figure 7) are associated with the formations of new hydrogen bonds between the C-H groups of testosterone ligand and O6 atom of G28 base, and between the C-H groups of pyridyl ligand and O6 atom of G15 base with 20.82% and 13.78% occupancy times, respectively, for the Pt(Testo)(II) + DNA adduct (see Figure 6(a) and Table 2 ). These calculated results suggest that the DNA base dynamics motions caused by DNA-bound testosterone-based platinum agents induce DNA conformational distortion via the rearrangement of hydrogen bonds at the vicinity of testosterone position.
The deviation extents of DNA dynamics parameters in the Im-Pt(Testo)(II) + DNA model away from the B-DNA molecule are generally larger than those in the Pt(Testo)(II) + DNA model. For example, the deviations of twist angle values for the C26 · A27 base pair step in the Im-Pt(Testo)(II) + DNA model increases by 80% (from ~ -25º to ~ -45º) with respect to that in the Pt(Testo)(II) + DNA model. Especially, the deviation directions of some DNA dynamics parameters in the Im-Pt(Testo)(II) + DNA model are opposite to those in the Pt(Testo)(II) + DNA model. For example, the deviation direction of opening motion with ~65° at the T12:A29 base pair in the Im-Pt(Testo)(II) + DNA model is opposite to that with ~ -50° at the T14:A27 base pair in the Pt(Testo)(II) + DNA model, which suggests the difference of groove-face and intercalative interaction modes in the platinum + DNA adducts. These observations result from the intercalative interaction mode making the base pair step opening to the DNA major groove by a positive value, and oppositely the groove-face interaction mode making that opening to the minor groove by a negative value due to the testosterone ligand locating at the DNA major groove surface. These observations suggest that the different interaction modes of platinum agents with DNA molecules caused by the flexibility of platinum agents might greatly affect the DNA base dynamics motions and DNA conformational damage. Moreover, it can be seen in the present study that the intercalative interaction mode of platinum agent with DNA molecule might induce extensive DNA conformation distortion over the groove-face interaction mode in the platinum agent + DNA adducts.
Molecular dynamics simulations and DNA dynamics analyses for a series of the Pt + DNA adducts were carried out to examine the distortions of DNA double-helical structure perturbed by the binding of testosterone-based platinum agent Pt(Testo)(II) (cis-[Pt(NH3)2(17α-pyridyl-3-ethyltestosterone)](II)), improved testosterone-based platinum agent Im-Pt(Testo)(II) (cis-[Pt(NH3)2(17α-pyridyl-3-ethynyltestosterone)](II)) and non-testosterone platinum agent Pt(Py)(II) (cis-[Pt(NH3)2(pyridine)](II)). It has been found that the rigid testosterone-based platinum agent Pt(Testo)(II) interacts with the DNA molecule via the groove-face interaction mode at the major groove of DNA; however, the improved flexible testosterone-based platinum agent Im-Pt(Testo)(II) interacts with DNA molecule via the intercalative interaction mode with the testosterone ligand inserting the middle of T12 and C13 bases. The distortion of DNA helical conformation caused by Im-Pt(Testo)(II) agent with intercalative interaction mode is larger than that caused by the Pt(Testo)(II) agent with groove-face interaction mode, which is supported by the formations of more hydrogen bonds and hydrophobic interactions at the platinum agent-DNA interface in the Im-Pt(Testo)(II) + DNA adduct. It can be found through the DNA dynamics analysis that the DNA base motions of opening, shift, rise, roll and twist away from the DNA groove caused by the binding of Pt(Testo)(II) and Im-Pt(Testo)(II) agents result in the widening and shoaling of DNA major/minor groove, and hydrogen bond destruction of DNA base pairs. Moreover, the non-testosterone-based platinum agent Pt(Py)(II) might cause insignificant change of DNA conformation due to the absence of testosterone ligand. Our simulation results might provide useful insights into understanding how a DNA conformation is affected by a testosterone-based platinum complex at the atomic level, and might be helpful for anticancer agent design.
Klein AV, Hambley TW: Platinum drug distribution in cancer cells and tumors. Chem Rev 2009, 109(10):4911–4920. 10.1021/cr9001066
Lovejoy KS, Lippard SJ: Non-traditional platinum compounds for improved accumulation, oral bioavailability, and tumor targeting. Dalton Trans 2009, 48: 10651–10659.
Wheate NJ, Walker S, Craig GE, Oun R: The status of platinum anticancer drugs in the clinic and in clinical trials. Dalton Trans 2010, 39(35):8113–8127. 10.1039/c0dt00292e
Wong E, Giandomenico CM: ChemInform abstract: current status of platinum-based antitumor drugs. ChemInform 1999, 30(48):no-no.
Rodger A, Patel KK, Sanders KJ, Datt M, Sacht C, Hannon MJ: Anti-tumour platinum acylthiourea complexes and their interactions with DNA. J Chem Soc Dalton Trans 2002, 19: 3656–3663.
Hannon MJ: Metal-based anticancer drugs: from a past anchored in platinum chemistry to a post-genomic future of diverse chemistry and biology. Pure Appl Chem 2007, 79(12):2243–2261. 10.1351/pac200779122243
Kostova I: Platinum complexes as anticancer agents. Recent Pat Anticancer Drug Discov 2006, 1(1):1–22.
Reedijk J: Platinum anticancer coordination compounds: study of dna binding inspires new drug design. Eur J Inorg Chem 2009, 2009(10):1303–1312. 10.1002/ejic.200900054
Sherman SE, Gibson D, Wang A, Lippard SJ: X-ray structure of the major adduct of the anticancer drug cisplatin with DNA: cis-[Pt (NH3) 2 (d (pGpG))]. Science (New York, NY) 1985, 230(4724):412. 10.1126/science.4048939
Abu-Surrah AS, Kettunen M: Platinum group antitumor chemistry: design and development of new anticancer drugs complementary to cisplatin. Curr Med Chem 2006, 13(11):1337–1357. 10.2174/092986706776872970
Bruijnincx PC, Sadler PJ: New trends for metal complexes with anticancer activity. Curr Opin Chem Biol 2008, 12(2):197–206. 10.1016/j.cbpa.2007.11.013
Harper BW, Krause-Heuer AM, Grant MP, Manohar M, Garbutcheon-Singh KB, Aldrich-Wright JR: Advances in platinum chemotherapeutics. Chemistry 2010, 16(24):7064–7077.
McWhinney SR, Goldberg RM, McLeod HL: Platinum neurotoxicity pharmacogenetics. Mol Cancer Ther 2009, 8(1):10–16.
Kandala PK, Srivastava SK: Diindolylmethane suppresses ovarian cancer growth and potentiates the effect of cisplatin in tumor mouse model by targeting signal transducer and activator of transcription 3 (STAT3). BMC Med 2012, 10(1):9. 10.1186/1741-7015-10-9
Zhang J, Wang X, Tu C, Lin J, Ding J, Lin L, Wang Z, He C, Yan C, You X: Monofunctional platinum complexes showing potent cytotoxicity against human liver carcinoma cell line BEL-7402. J Med Chem 2003, 46(16):3502–3507. 10.1021/jm020593j
Sanchez-Cano C, Hannon MJ: Cytotoxicity, cellular localisation and biomolecular interaction of non-covalent metallo-intercalators with appended sex hormone steroid vectors. Dalton Trans 2009, 48: 10765–10773.
Coll M, Frederick CA, Wang A, Rich A: A bifurcated hydrogen-bonded conformation in the d (AT) base pairs of the DNA dodecamer d (CGCAAATTTGCG) and its complex with distamycin. Proc Natl Acad Sci 1987, 84(23):8385–8389. 10.1073/pnas.84.23.8385
Gessner RV, Quigley GJ, Wang AHJ, Van der Marel GA, Van Boom JH, Rich A: Structural basis for stabilization of Z-DNA by cobalt hexaammine and magnesium cations. Biochemistry 1985, 24(2):237–240. 10.1021/bi00323a001
Łęczkowska A, Vilar R: Interaction of metal complexes with nucleic acids. Annu Rep Prog Chem Sect A: Inorg Chem 2012, 108(1):330–349.
Lipscomb LA, Zhou FX, Presnell SR, Woo RJ, Peek ME, Plaskon RR, Williams LD: Structure of a DNA-porphyrin complex. Biochemistry 1996, 35(9):2818–2823. 10.1021/bi952443z
Liu HK, Sadler PJ: Metal complexes as DNA intercalators. Accounts Chem Res 2011, 44(5):349–359. 10.1021/ar100140e
Nair RB, Teng ES, Kirkland SL, Murphy CJ: Synthesis and DNA-binding properties of [Ru (NH3) 4dppz] 2+. Inorg Chem 1998, 37(1):139–141. 10.1021/ic970432j
Lovejoy KS, Todd RC, Zhang S, McCormick MS, D’Aquino JA, Reardon JT, Sancar A, Giacomini KM, Lippard SJ: cis-Diammine(pyridine)chloroplatinum(II), a monofunctional platinum(II) antitumor agent: uptake, structure, function, and prospects. Proc Natl Acad Sci USA 2008, 105(26):8902–8907. 10.1073/pnas.0803441105
Gagnon V, St-Germain ME, Descoteaux C, Provencher-Mandeville J, Parent S, Mandal SK, Asselin E, Berube G: Biological evaluation of novel estrogen-platinum(II) hybrid molecules on uterine and ovarian cancers-molecular modeling studies. Bioorg Med Chem Lett 2004, 14(23):5919–5924. 10.1016/j.bmcl.2004.09.015
Hannon MJ, Green PS, Fisher DM, Derrick PJ, Beck JL, Watt SJ, Ralph SF, Sheil MM, Barker PR, Alcock NW: An estrogen-platinum terpyridine conjugate: DNA and protein binding and cellular delivery. Chemistry 2006, 12(31):8000–8013. 10.1002/chem.200501012
Huxley M, Sanchez-Cano C, Browning MJ, Navarro-Ranninger C, Quiroga AG, Rodger A, Hannon MJ: An androgenic steroid delivery vector that imparts activity to a non-conventional platinum(II) metallo-drug. Dalton Trans 2010, 39(47):11353–11364. 10.1039/c0dt00838a
Sanchez-Cano C, Huxley M, Ducani C, Hamad AE, Browning MJ, Navarro-Ranninger C, Quiroga AG, Rodger A, Hannon MJ: Conjugation of testosterone modifies the interaction of mono-functional cationic platinum(II) complexes with DNA, causing significant alterations to the DNA helix. Dalton Trans 2010, 39(47):11365–11374. 10.1039/c0dt00839g
Ruiz J, Rodriguez V, Cutillas N, Espinosa A, Hannon MJ, Novel C: N-chelate platinum(II) antitumor complexes bearing a lipophilic ethisterone pendant. J Inorg Biochem 2011, 105(4):525–531. 10.1016/j.jinorgbio.2010.12.005
Baruah H, Wright MW, Bierbach U: Solution structural study of a DNA duplex containing the guanine-N7 adduct formed by a cytotoxic platinum-acridine hybrid agent. Biochemistry 2005, 44(16):6059–6070. 10.1021/bi050021b
Chaney SG, Ramachandran S, Sharma S, Dokholyan NV, Temple B, Bhattacharyya D, Wu Y, Campbell S: Differences in conformation and conformational dynamics between cisplatin and oxaliplatin DNA adducts. In Platinum and Other Heavy Metal Compounds in Cancer Chemotherapy. Edited by: Bonetti A, Leone R, Muggia F, Howell S. Humana Press; 2009:157–169.
Mantri Y, Lippard SJ, Baik MH: Bifunctional binding of cisplatin to DNA: why does cisplatin form 1, 2-intrastrand cross-links with AG but not with GA? J Am Chem Soc 2007, 129(16):5023–5030. 10.1021/ja067631z
Ramachandran S, Temple BR, Chaney SG, Dokholyan NV: Structural basis for the sequence-dependent effects of platinum-DNA adducts. Nucleic Acids Res 2009, 37(8):2434–2448. 10.1093/nar/gkp029
Sharma S, Gong P, Temple B, Bhattacharyya D, Dokholyan NV, Chaney SG: Molecular dynamic simulations of cisplatin- and oxaliplatin-d(GG) intrastand cross-links reveal differences in their conformational dynamics. J Mol Biol 2007, 373(5):1123–1140. 10.1016/j.jmb.2007.07.079
Spingler B, Whittington DA, Lippard SJ: 2.4 Å crystal structure of an oxaliplatin 1, 2-d (GpG) intrastrand cross-link in a DNA dodecamer duplex. Inorg Chem 2001, 40(22):5596–5602. 10.1021/ic010790t
Zhu J, Zhao Y, Zhu Y, Wu Z, Lin M, He W, Wang Y, Chen G, Dong L, Zhang J: DNA cross-linking patterns induced by an antitumor-active trinuclear platinum complex and comparison with its dinuclear analogue. Chemistry 2009, 15(21):5245–5253. 10.1002/chem.200900217
Zhu Y, Wang Y, Chen G: Differences in conformational dynamics of [Pt3(HPTAB)]6 + -DNA adducts with various cross-linking modes. Nucleic Acids Res 2009, 37(17):5930–5942. 10.1093/nar/gkp618
Marhefka CA, Moore BM II, Bishop TC, Kirkovsky L, Mukherjee A, Dalton JT, Miller DD: Homology modeling using multiple molecular dynamics simulations and docking studies of the human androgen receptor ligand binding domain bound to testosterone and nonsteroidal ligands. J Med Chem 2001, 44(11):1729–1740. 10.1021/jm0005353
Zhou H, Singh N, Kim KS: Homology modeling and molecular dynamics study of chorismate synthase from Shigella flexneri . J Mol Graph Model 2006, 25(4):434–441. 10.1016/j.jmgm.2006.02.013
Prusis P, Schiöth HB, Muceniece R, Herzyk P, Afshar M, Hubbard RE, Wikberg JES: Modeling of the three-dimensional structure of the human melanocortin 1 receptor, using an automated method and docking of a rigid cyclic melanocyte-stimulating hormone core peptide. J Mol Graph Model 1997, 15(5):307–317. 10.1016/S1093-3263(98)00004-7
Yang X, Wang Z, Dong W, Ling L, Yang H, Chen R: Modeling and docking of the three-dimensional structure of the human melanocortin 4 receptor. J Protein Chem 2003, 22(4):335–344. 10.1023/A:1025386022852
Case D, Darden T, Cheatham T, Simmerling C, Wang J, Duke R, Luo R, Merz K, Pearlman D, Crowley M: AMBER 9. San Francisco: University of California; 2006.
Kobayashi A, Ohbayashi K, Aoki R, Chang HC, Kato M: Synthesis, structure and photophysical properties of a flavin-based platinum (II) complex. Dalton Trans 2011, 40(14):3484–3489. 10.1039/c0dt01139h
Zhou L: Theoretical analysis on the transition state of the anticancer drug trans-[PtCl2 (isopropylamine) 2] and its cis isomer binding to DNA purine bases. J Phys Chem B 2009, 113(7):2110–2127. 10.1021/jp806661g
Frisch MJ, Trucks G, Schlegel H, Scuseria G, Robb M, Cheeseman J, Scalmani G, Barone V, Mennucci B, Petersson G: Gaussian 09; Gaussian. CT: Inc, Wallingford; 2009.
Zhu Y, Su Y, Li X, Wang Y, Chen G: Evaluation of Amber force field parameters for copper (II) with pyridylmethyl-amine and benzimidazolylmethyl-amine ligands: a quantum chemical study. Chem Phys Lett 2008, 455(4):354–360. 10.1016/j.cplett.2008.03.004
Zhu Y, Wang Y, Chen G, Zhan CG: A three-point method for evaluations of AMBER force field parameters: an application to copper-based artificial nucleases. Theor Chem Acc: Theory, Computation, and Modeling (Theoretica Chimica Acta) 2009, 122(3):167–178. 10.1007/s00214-008-0496-6
Duan Y, Wu C, Chowdhury S, Lee MC, Xiong G, Zhang W, Yang R, Cieplak P, Luo R, Lee T: A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculations. J Comput Chem 2003, 24(16):1999–2012. 10.1002/jcc.10349
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 2004, 55(3):620–634. 10.1002/prot.10470
Yao S, Plastaras JP, Marzilli LG: A molecular mechanics AMBER-type force field for modeling platinum complexes of guanine derivatives. Inorg Chem 1994, 33(26):6061–6077. 10.1021/ic00104a015
Wang J, Wolf RM, Caldwell JW, Kollman PA, Case DA: Development and testing of a general amber force field. J Comput Chem 2004, 25(9):1157–1174. 10.1002/jcc.20035
Lavery R, Sklenar H: The definition of generalized helicoidal parameters and of axis curvature for irregular nucleic acids. J Biomol Struct Dyn 1988, 6(1):63–91. 10.1080/07391102.1988.10506483
Strahs D, Schlick T: A-tract bending: insights into experimental structures by computational models. J Mol Biol 2000, 301(3):643–664. 10.1006/jmbi.2000.3863
Yang B, Zhu Y, Wang Y, Chen G: Interaction identification of Zif268 and TATAZF proteins with GC‐/AT‐rich DNA sequence: a theoretical study. J Comput Chem 2010, 32(3):416–428.
Goodsell DS, Dickerson RE: Bending and curvature calculations in B-DNA. Nucleic Acids Res 1994, 22(24):5497. 10.1093/nar/22.24.5497
Poncin M, Piazzola D, Lavery R: DNA flexibility as a function of allomorphic conformation and of base sequence. Biopolymers 1992, 32(8):1077–1103. 10.1002/bip.360320817
Stofer E, Lavery R: Measuring the geometry of DNA grooves. Biopolymers 1994, 34(3):337–346. 10.1002/bip.360340305
The authors acknowledge research support from the National Science Foundation of China (No. 21271029, 21131003, 21073015 and 20973024), the Major State Basic Research Development Programs (grant No. 2011CB808500).
The authors declare that they have no competing interests.
SC carried out the computation and analysis of data, and drafted the manuscript. YW conceived this work and critically revised the manuscript. GC have made substantial contributions to the interpretation and evaluation of the results, and helped in construction of the manuscript. We also wish to thank the advice of our reviewers. All authors read and approved the manuscript.
Electronic supplementary material
Additional file 7: Figure S2. Frequency distributions of DNA duplex base-pair helical/step parameters for B-DNA and Im-Pt(Testo)(II) + DNA adduct. Selected frequency distributions of the representative DNA duplex base-pair helical/step parameters for at the base pairs near the binding site and testosterone position for B-DNA (black line) and Im-Pt(Testo)(II) + DNA adduct (Magenta line). (DOC 744 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
- Molecular dynamics simulations
- Groove-face and intercalative interactions
- Testosterone-based platinum agent
- Pt + DNA adducts
- DNA conformation distortion