Initial DNA Interactions of the Binuclear Threading Intercalator Λ,Λ‐[μ‐bidppz(bipy)4Ru2]4+: An NMR Study with [d(CGCGAATTCGCG)]2

Binuclear polypyridine ruthenium compounds have been shown to slowly intercalate into DNA, following a fast initial binding on the DNA surface. For these compounds, intercalation requires threading of a bulky substituent, containing one Ru^II, through the DNA base‐pair stack, and the accompanying DNA duplex distortions are much more severe than with intercalation of mononuclear compounds. Structural understanding of the process of intercalation may greatly gain from a characterisation of the initial interactions between binuclear Ru^II compounds and DNA. We report a structural NMR study on the binuclear Ru^II intercalator Λ,Λ‐B (Λ,Λ‐[μ‐bidppz(bipy)₄Ru₂]⁴⁺; bidppz=11,11′‐bis(dipyrido[3,2‐a:2′,3′‐c]phenazinyl, bipy = 2,2′‐bipyridine) mixed with the palindromic DNA [d(CGCGAATTCGCG)]₂. Threading of Λ,Λ‐B depends on the presence and length of AT stretches in the DNA. Therefore, the latter was selected to promote initial binding, but due to the short stretch of AT base pairs, final intercalation is prevented. Structural calculations provide a model for the interaction: Λ,Λ‐B is trapped in a well‐defined surface‐bound state consisting of an eccentric minor‐groove binding. Most of the interaction enthalpy originates from electrostatic and van der Waals contacts, whereas intermolecular hydrogen bonds may help to define a unique position of Λ,Λ‐B. Molecular dynamics simulations show that this minor‐groove binding mode is stable on a nanosecond scale. To the best of our knowledge, this is the first structural study by NMR spectroscopy on a binuclear Ru compound bound to DNA. In the calculated structure, one of the positively charged Ru²⁺ moieties is near the central AATT region; this is favourable in view of potential intercalation as observed by optical methods for DNA with longer AT stretches. Circular dichroism (CD) spectroscopy suggests that a similar binding geometry is formed in mixtures of Λ,Λ‐B with natural calf thymus DNA. The present minor‐groove binding mode is proposed to represent the initial surface interactions of binuclear Ru^II compounds prior to intercalation into AT‐rich DNA.     

Binding of piano‐stool Ru(II) complexes to DNA; QM/MM study

Ru(II) “piano‐stool” complexes belong to group of biologically active metallocomplexes with promising anticancer activity. In this study, we investigate the reaction mechanism of [(η⁶‐benzene)Ru(II)(en)(H₂O)]²⁺ (en = ethylenediamine) complex binding to DNA by hybrid QM/MM computational techniques. The reaction when the Ru(II) complex is coordinated on N7‐guanine from major groove is explored. Two reaction pathways, direct binding to N7 position and two‐step mechanism passing through O6 position, are considered. It was found that the reaction is exothermic and the direct binding process is preferred kinetically. In analogy to cisplatin, we also explored the possibility of intrastrand cross‐link formation where the Ru(II) complex makes a bridge between two adjacent guanines. Two different pathways were found, leading to a final structure with released benzene ligand. This process is exothermic; however, one pathway is blocked by relatively high initial activation barrier. Geometries, energies, and electronic properties analyzed by atoms in molecules and natural population analysis methods are discussed. © 2012 Wiley Periodicals, Inc.     

诺氟沙星-DNA复合物的分子动力学模拟

采用分子模建的方法构建了诺氟沙星-DNA复合物的初始结构, 通过2 ns的分子动力学(MD)模拟研究表明: 诺氟沙星能够和双螺旋d[ATATCGATAT]₂形成稳定的复合物, 药物分子可紧密结合在DNA的小沟区域, 并且能够与DNA的鸟嘌呤碱基形成两个稳定的氢键. 在分子水平上提供了诺氟沙星直接与双螺旋DNA相互作用的结构及复合物的动态变化情况.     

DNA小沟结合二脒:水分子的识别作用

采用分子动力学模拟了DB921-DNA复合物, 通过7 ns的模拟研究表明: DB921一端的氨基氮原子与一个水分子形成氢键, 同时, 水分子又与DNA的5位A碱基的氮原子形成一个氢键. 水分子在DB921与DNA小沟结合中起了桥连的作用, 使得直线型的芳香二脒化合物DB921通过水桥与DNA小沟结合, 水分子诱导DB921分子与DNA的小沟域构型相适应, 与DNA小沟域的AATTC碱基有较强的结合作用. 在分子水平上提供了DB921与双螺旋DNA相互作用的结构及复合物的动态变化情况, 指出水分子在DNA小沟结合二脒化合物中的识别作用, 为设计出更高生物活性的DNA小沟结合剂提供一定的理论依据.     

DNA小沟结合药物DB818的分子动力学模拟和结合自由能分析

采用分子动力学模拟了DNA小沟结合芳香二脒药物DB818形成的复合物. 通过5 ns的模拟研究表明: DB818药物分子可紧密结合在DNA的AATTC小沟区域, 和双螺旋d[CGCGAATTCGCG]2形成稳定的复合物. 由于噻吩硫原子的弱电负性, 使DB818能够以更大的伸展程度与DNA的小沟结合, 形成更强的结合力. DB818苯并咪唑的氮原子能够与DNA 7位和19位T碱基上的氧原子形成两个稳定的氢键, 同时, DB818末端氨基氮原子分别与DNA 的20位T碱基的氧原子和9位C碱基的氧原子形成两个氢键. 另外, 运用MM_PBSA方法计算了DB293-DNA和DB818-DNA复合物的结合自由能, 计算结合能与实验值能较好的吻合, 通过比较其结合自由能, 从热力学能量角度说明了DB818有较大的熵值与较小的焓值贡献, 从而与DNA小沟结合的结合力比DB293强. 本文在分子水平上提供了DB818直接与双螺旋DNA相互作用的结构及复合物的动态变化情况, 为设计出更高生物活性的DNA小沟结合剂提供一定的理论依据.     

Monitoring the DNA binding kinetics of a binuclear ruthenium complex by energy transfer: evidence for slow shuffling

The semirigid binuclear ruthenium complex Delta,Delta-[mu-(11,11'-bidppz)(phen)(4)Ru(2)](4+) has been shown to rearrange slowly from an initial groove-bound nonluminescent state to a final intercalated emissive state by threading one of its bulky Ru(phen)(2) moieties through the DNA base stack. When this complex binds to poly[d(A-T)(2)], a further increase in emission from the complex is observed after completion of the intercalation, assigned to reorganization of the intercalated complex. We here report a study of the threading process in poly[d(A-T)(2)], in which the minor groove binding dye DAPI is used as an energy transfer probe molecule to assess the distribution of ruthenium complex during and also after the actual threading phase. The emission from DAPI is found to change with the same rate as the emission from the ruthenium complex, and furthermore, DAPI does not disturb the binding kinetics of the latter, justifying it as a good probe of both the threading and the reorganization processes. We conclude from the change in the emission from both DAPI and the ruthenium complex with time that DAPI-ruthenium interactions are most pronounced during the process of threading of the complex, suggesting that the complexes are initially threaded slightly anticooperatively and thereafter redistribute along the DNA to reach their thermodynamically most favorable distribution. The final distribution is characterized by a small but significant binding cooperativity, probably as a result of hydrophobic interactions between the complex ions despite their tetravalent positive charges. The mechanism of "shuffling" the complex along the DNA chain is discussed, i.e., whether the ruthenium complex remains threaded (requiring sequential base-pair openings) or if unthreading followed by lateral diffusion within the ionic atmosphere of the DNA and rethreading occurs.     

DNA与二脒DB293复合物的分子动力学模拟

采用分子动力学模拟了DNA小沟与芳香二脒化合物DB293结合形成的复合物,通过5ns的模拟研究表明,DB293分子可紧密结合在DNA的AATT小沟区域,和双螺旋d[CGCGAATTCGCG]2形成稳定的复合物。DB293苯并咪唑的氮原子N2能够与DNA胸腺嘧啶碱基T7的O2原子和T19的O2原子形成两个较强的氢键,同时,其末端氨基的N3原子和T20的O2原子形成一个较弱的氢键。本文在分子水平上提供了DB293直接与双螺旋DNA相互作用的结构及复合物的动态变化情况,为设计出更高活性的芳香二脒类DNA小沟结合剂提供一定的理论依据。     

DNA sequence dependent monomer-dimer binding modulation of asymmetric benzimidazole derivatives

A number of studies indicate that DNA sequences such as AATT and TTAA have significantly different physical and interaction properties. To probe these interaction differences in detail and determine the influence of charge, we have synthesized three bisbenzimidazole derivatives, a diamidine, DB185, and monoamidines, DB183 and DB210, that are related to the well-known minor groove agent, Hoechst 33258. Footprinting studies with several natural and designed DNA fragments indicate that the synthetic compounds bind at AT sequences in the minor groove and interact more weakly at sites with TpA steps relative to sites without such steps. Circular dichroism spectroscopy also indicates that the compounds bind in the DNA minor groove. Surprisingly, Tm studies as a function of ratio indicate that the monoamidines bind to TTAA sequences as dimers, whereas the diamidine binds as a monomer. Biosensor-surface plasmon resonance (SPR) studies allowed us to quantitate the interaction differences in more detail. SPR results clearly show that the monoamidine compounds bind to the TTAA sequence in a cooperative 2:1 complex but bind as monomers to AATT. The dication binds to both sequences in monomer complexes but the binding to AATT is significantly stronger than binding to TTAA. Molecular dynamics simulations indicate that the AATT sequence has a narrow time-average minor groove width that is a very good receptor site for the bisbenzimidazole compounds. The groove in TTAA sequences is wider and the width must be reduced to form a favorable monomer complex. The monocations thus form cooperative dimers that stack in an antiparallel orientation and closely fit the structure of the TTAA minor groove. The amidine groups in the dimer are oriented in the 5' direction of the strand to which they are closest. Charge repulsion in the dication apparently keeps it from forming the dimer. It instead reduces the TTAA groove width, in an induced fit process, sufficiently to form a minor groove complex. The dimer-binding mode of DB183 and DB210 is a new DNA recognition motif and offers novel design concepts for selective targeting of DNA sequences with a wider minor groove, including those with TpA steps.     

Azole-bridged diplatinum anticancer compounds. Modulating DNA flexibility to escape repair mechanism and avoid cross resistance

Dinuclear azole-bridged Pt compounds bind to DNA helices, forming intrastrand crosslinks between adjacent guanines in a similar way to cisplatin. Their cytotoxic profile is, however, different from that of first and second generation Pt drugs in that they lack cross resistance in cisplatin-resistant cell lines. In contrast to cisplatin, which induces a large kink in DNA duplex, structural NMR studies and molecular dynamics simulations have shown that azole-bridged diplatinum compounds induce only small structural changes in double-stranded DNA. These structural differences have been invoked to explain the different cytotoxic profile of these compounds. Here, we show that in addition to the small structural changes in DNA, dinuclear Pt compounds also affect DNA minor groove flexibility in a different way than cisplatin. Free-energy calculations on azole-bridged diplatinum DNA adducts reveal that opening of the minor groove requires a higher free-energy cost (DeltaG ~ 7-15 kcal/mol) than in the corresponding cisplatin-DNA adduct (DeltaG ~ 0 kcal/mol). This could prevent minor groove binding proteins from binding to diplatinum-DNA adducts thus leading to a different cellular response than cisplatin and possibly decreasing the activity of excision repair enzymes. Although the development of drug resistance is a highly complex mechanism, our findings provide an additional rationale for the improved cytotoxic activity of these compounds in cell lines resistant to cisplatin.     

Molecular modeling of cationic porphyrin-anthraquinone hybrids as DNA topoisomerase IIβ inhibitors

Human DNA Topoisomerase II has been regarded as a promising target in anticancer drug discovery. In the present study, we designed six porphyrin-anthraquinone hybrids bearing pyrazole or pyridine group as meso substituents and evaluated their potentials as DNA Topoisomerase IIβ inhibitor. First, we investigated the binding orientation of porphyrin hybrids into DNA topoisomerase IIβ employing AutoDock 4.2 and then performed 20-ns molecular dynamics simulations to see the dynamic stability of each porphyrin-Topo IIβ complex using Amber 14. We found that the binding of porphyrin hybrids occured through intercalation and groove binding mode in addition interaction with the amino acid residues constituting the active cavity of Topo IIβ. Each porphyrin-Topo IIβ complex was stabilized during 20-ns dynamics simulations. The MM-PBSA free energy calculation shows that the binding affinities of porphyrin hybrids were modified with the number of meso substituent. Interestingly, the affinity of all porphyrin hybrids to Topo IIβ was stronger than that of native ligand (EVP), indicating the potential of the designed porphyrin to be considered in experimental research.     

Electrostatic recognition and induced fit in the kappa-PVIIA toxin binding to Shaker potassium channel

Brownian dynamics (BD) and molecular dynamics (MD) simulations and electrostatic calculations were performed to study the binding process of kappa-PVIIA to the Shaker potassium channel and the structure of the resulting complex. BD simulations, guided by electrostatic interactions, led to an initial alignment between the toxin and the channel protein. MD simulations were then carried out to allow for rearrangements from this initial structure. After approximately 4 ns, a critical "induced fit" process was observed to last for approximately 2 ns. In this process, the interface was reorganized, and side chains were moved so that favorable atomic contacts were formed or strengthened, while unfavorable contacts were eliminated. The final complex structure was stabilized through electrostatic interactions with the positively charged side chain of Lys7 of kappa-PVIIA deeply inserted into the channel pore and other hydrogen bonds and by hydrophobic interactions involving Phe9 and Phe23 of the toxin. The validity of the predicted structure for the complex was assessed by calculating the effects of mutating charged and polar residues of both the toxin and the channel protein, with the calculated effects correlating reasonably well with experimental data. The present study suggests a general binding mechanism, whereby proteins are pre-aligned in their diffusional encounter by long-range electrostatic attraction, and nanosecond-scale rearrangements within the initial complex then lead to a specifically bound complex.     

DNA-binding of semirigid binuclear ruthenium complex delta,delta-[mu-(11,11'-bidppz)(phen)(4)ru(2)](4+): extremely slow intercalation kinetics

We here report a remarkably slow rearrangement of binding modes for a binuclear ruthenium(II) complex upon interaction with DNA. It has been previously shown that Delta,Delta-[mu-(11,11'-bidppz)(phen)4Ru2]4+ binds to DNA in one of the grooves. However, we find that this is only an initial, metastable, binding mode, which is extremely slowly reorganized into an intercalative binding geometry. The slow rearrangement and dissociation, revealed by flow linear dichroism and fluorescence spectroscopy, are concluded to be a result from the complex being threaded through the DNA, with one of the bridging aromatic dppz ligands intercalated between the base pairs of the DNA, placing one metal center in the minor groove and one in the major groove. A negative LD, a high luminescence quantum yield, and long luminescence lifetimes, similar to the intercalating complex Delta-[Ru(phen)2dppz]2+, indicate intercalation of the bidppz moiety. The unique slow dissociation of the complex in its final DNA-binding mode suggests that this class of threading, partially intercalated binuclear complexes may be interesting in the context of cancer therapy. Also, their unique optical and photophysical properties could make such complexes, either alone or scaffolded by DNA structures, of interest for the development of nanometer-sized molecular optoelectronic devices.     

DNA recognition by the anthracycline antibiotic respinomycin D: NMR structure of the intercalation complex with d(AGACGTCT)2

Respinomycin D is a member of the anthracycline family of antitumour antibiotics that interact with double stranded DNA through intercalation. The clinical agents daunomycin and doxorubicin are the most well-studied of this class but have a relatively simple molecular architecture in which the pendant daunosamine sugar resides in the DNA minor groove. Respinomycin D, which belongs to the nogalamycin group of anthracyclines, possesses additional sugar residues at either end of the aglycone chromophore that modulate the biological activity but whose role in molecular recognition is unknown. We report the NMR structure of the respinomycin D-d(AGACGTCT)2 complex in solution derived from NOE restraints and molecular dynamics simulations. We show that the drug threads through the DNA double helix forming stabilising interactions in both the major and minor groove, the latter through a different binding geometry to that previously reported. The bicycloaminoglucose sugar resides in the major groove and makes specific contacts with guanine at the 5'-CpG intercalation site, however, the disaccharide attached at the C4 position plays little part in drug binding and DNA recognition and is largely solvent exposed.     

AFM‐Correlated CSM‐Coupled Raman/Fluorescence Studies on Selective Interactions of Water Soluble Porphyrins with DNA

The Raman and fluorescence spectroscopic properties of water‐soluble oxo‐titanium(IV) mesotetrakis (1‐methyl pyridium‐4‐yl) porphyrin (O=Ti(TMPyP)⁴⁺) bound with calf thymus DNA and artificial DNAs such as double stranded poly[d(A‐T)₂] and poly[d(G‐C)₂] have been investigated on the single DNA molecule basis by AFM‐correlated confocal scanning microscope (CSM)‐coupled Raman and fluorescence spectroscopic techniques as well as the ensemble‐averaged spectroscopy. The ensemble‐averaged spectroscopic studies imply that the porphyrin interacts with DNA in different groove binding patterns depending on the base pairs. AFM‐images of the different DNAs bound with O=Ti(TMPyP)⁴⁺ were measured, and their morphologies are found to depend on kind of base pairs interacting with O=Ti(TMPyP)⁴⁺. Being correlated with the AFM images, the CSM‐coupled Raman and fluorescence spectral properties of the three different single O=Ti(TMPyP)⁴⁺‐DNA complexes were observed to be highly resolved and sensitive to base pair‐dependent axial ligation of Ti‐O bond as compared to the corresponding ensemble‐averaged spectral properties, which affect the groove binding and its strength of the O=Ti(TMPyP)⁴⁺ with DNA. The axial ligation was found to be accompanied by vibration structural change of the porphyrin ring, leading to keep the shape of double stranded poly[d(A‐T)₂] rigid while poly‐[d(G‐C)₂] and calf thymus DNA flexible after binding with the oxo‐titanyl porphyrin. The base pair dependence of the fluorescence decay times of the DNA‐bound porphyrins was also observed, implying that an excited‐state charge transfer takes place in the G‐C rich major groove in calf thymus DNA. These results suggest that binding of O=Ti(TMPyP)⁴⁺ is more preferential with the G‐C rich major groove than with the A‐T rich minor groove in calf thymus DNA so that the morphology of DNA is changed.     

Prediction of protein–protein complexes using replica exchange with repulsive scaling

The realistic prediction of protein–protein complex structures is import to ultimately model the interaction of all proteins in a cell and for the design of new protein–protein interactions. In principle, molecular dynamics (MD) simulations allow one to follow the association process under realistic conditions including full partner flexibility and surrounding solvent. However, due to the many local binding energy minima at the surface of protein partners, MD simulations are frequently trapped for long times in transient association states. We have designed a replica-exchange based scheme employing different levels of a repulsive biasing between partners in each replica simulation. The bias acts only on intermolecular interactions based on an increase in effective pairwise van der Waals radii (repulsive scaling (RS)-REMD) without affecting interactions within each protein or with the solvent. For a set of five protein test cases (out of six) the RS-REMD technique allowed the sampling of near-native complex structures even when starting from the opposide site with respect to the native binding site for one partner. Using the same start structures and same computational demand regular MD simulations sampled near native complex structures only for one case. The method showed also improved results for the refinement of docked structures in the vicinity of the native binding geometry compared to regular MD refinement.     

New transition metal complexes of 9,10‐phenanthrenequinone p‐toluyl hydrazone Schiff base: Synthesis,spectroscopy, DNA and HSA interactions,antimicrobial, DFT and docking studies

The new complexes of Cu (II) and Ni (II) of a tridentate Schiff base ligand derived from 9,10‐phenanthrenequinone and p‐toluic hydrazide have been synthesized and characterized by elemental analysis, electrical conductometry, FT‐IR, Mass, NMR and UV–Vis. The DFT calculations were carried out at B3LYP/6‐31G*(d) level for the determination of the optimized structure of the ligand and its complexes. The as‐synthesized compounds were screened for their antimicrobial activity. Also, their binding behavior with fish salmon‐DNA (FS‐DNA) and human serum albumin (HSA) were studied by different kinds of spectroscopic and molecular modeling techniques. The fluorescence data at different temperatures were applied in order to estimate the thermodynamics parameters of interactions of ligand and its complexes with DNA and HSA. The results showed that the as‐made compounds could bind to FS‐DNA and HSA via the groove binding as the major binding mode. According to molecular docking calculation and competitive binding experiments, these compounds bind to the minor groove of DNA and hydrophobic residues located in the subdomain IB of HSA. In addition, the molecular docking results kept in good consistence with experimental data.     

Sequence‐Specific DNA Recognition by Cyclic Pyrrole–Imidazole Cysteine‐Derived Polyamide Dimers

Pyrrole–imidazole (PI) polyamides bind to the minor groove of the DNA duplex in a sequence‐specific manner and thus have the potential to regulate gene expression. To date, various types of PI polyamides have been designed as sequence‐specific DNA binding ligands. One of these, cysteine cyclic PI polyamides containing two β‐alanine molecules, were designed to recognize a 7 bp DNA sequence with high binding affinity. In this study, an efficient cyclization reaction between a cysteine and a chloroacetyl residue was used for dimerization in the synthesis of a unit that recognizes symmetrical DNA sequences. To evaluate specific DNA binding properties, dimeric PI polyamide binding was measured by using a surface plasmon resonance (SPR) method. Extending this molecular design, we synthesized a large dimeric PI polyamide that can recognize a 14 bp region in duplex DNA.