Structural Flexibility of the d(CCAGTACTGG)2B-DNA Decamer and Its Complex with Two Polyamides

Dynamic hydrogen bonding and DNA flexibility in minor groove binders: molecular dynamics simulation of the polyamide f-ImPyIm bound to the Mlu1 (MCB) sequence 5'-ACGCGT-3' in 2:1 motif

Molecular dynamics simulations of the DNA 10-mer 5'-CCACGCGTGG-3' alone and complexed with the formamido-imidazole-pyrrole-imidazole (f-ImPyIm) polyamide minor groove binder in a 2:1 fashion were conducted for 50 ns using the pbsc0 parameters within the AMBER 12 software package. The change in DNA structure upon binding of f-ImPyIm was evaluated via minor groove width and depth, base pair parameters of Slide, Twist, Roll, Stretch, Stagger, Opening, Propeller, and x-displacement, dihedral angle distributions of ζ, ε, α, and γ determined using the Curves+ software program, and hydrogen bond formation. The dynamic hydrogen bonding between the f-ImPyIm and its cognate DNA sequence was compared to the static image used to predict sequence recognition by polyamide minor groove binders. Many of the predicted hydrogen bonds were present in less than 50% of the simulation; however, persistent hydrogen bonds between G5/15 and the formamido group of f-ImPyIm were observed. It was determined that the DNA is wider in the Complex than without the polyamide binder; however, there is flexibility in this particular sequence, even in the presence of the f-ImPyIm as evidenced by the range of minor groove widths the DNA exhibits and the dynamics of the hydrogen bonding that binds the two f-ImPyIm ions to the minor groove. The Complex consisting of the DNA and the 2 f-ImPyIm binders shows slight fraying of the 5' end of the 10-mer at the end of the simulation, but the portion of the oligomer responsible for recognition and binding is stable throughout the simulation. Several structural changes in the Complex indicate that minor groove binders may have a more active role in inhibiting transcription than just preventing binding of important transcription factors.

Modifying the N-terminus of polyamides: PyImPyIm has improved sequence specificity over f-ImPyIm

Seven N-terminus modified derivatives of a previously published minor-groove binding polyamide (f-ImPyIm, 1) were synthesized and the biochemical and biophysical chemistry evaluated. These compounds were synthesized with the aim of attaining a higher level of sequence selectivity over f-ImPyIm (1), a previously published strong minor-groove binder. Two compounds possessing a furan or a benzofuran moiety at the N-terminus showed a footprint of 0.5microM at the cognate ACGCGT site (determined by DNase I footprinting); however, the specificity of these compounds was not improved. In contrast, PyImPyIm (4) produced a footprint of 0.5microM but showed a superior specificity using the same technique. When evaluated by thermal melting experiments and circular dichroism using ACGCGT and the non-cognate AAATTT sequence, all compounds were shown to bind in the minor-groove of DNA and stabilize the cognate sequence much better than the non-cognate (except for the non-amido-compound that did not bind either sequence, as expected). PyImPyIm (4) was interesting as the DeltaT(m) for this compound was only 4 degrees C but the footprint was very selective. No binding was observed for this compound with a third DNA (non-cognate, ACCGGT). ITC studies on compound 4 showed exothermic binding with ACGCGT and no heat change was observed for titrating the compound to the other two DNA sequences. The heat capacity (DeltaC(p)) of the PIPI/ACGCGT complex calculated from the hydrophobic interactions and SASA calculations was comparable to the experimental value obtained from ITC (-146calmol(-1)K(-1)). SPR results provided confirmation of the sequence specificity of PyImPyIm (4), with a K(eq) value determined to be 7.1x10(6) M(-1) for the cognate sequence and no observable binding to AAATTT and ACCGGT. Molecular dynamic simulations affirmed that PyImPyIm (4) binds as a dimer in an overlapped conformation, and it fits snugly in the minor-groove of the ACGCGT oligonucleotide. PyImPyIm (4) is an especially interesting molecule, because although the binding affinity is slightly reduced, the specificity with respect to f-ImPyIm (1) is significantly improved.

Physical and Structural Basis for the Strong Interactions of the -ImPy- Central Pairing Motif in the Polyamide f-ImPyIm

The polyamide f-ImPyIm has a higher affinity for its cognate DNA than either the parent analogue, distamycin A (10-fold), or the structural isomer, f-PyImIm (250-fold), has for its respective cognate DNA sequence. These findings have led to the formulation of a two-letter polyamide "language" in which the -ImPy- central pairings associate more strongly with Watson-Crick DNA than -PyPy-, -PyIm-, and -ImIm-. Herein, we further characterize f-ImPyIm and f-PyImIm, and we report thermodynamic and structural differences between -ImPy- (f-ImPyIm) and -PyIm- (f-PyImIm) central pairings. DNase I footprinting studies confirmed that f-ImPyIm is a stronger binder than distamycin A and f-PyImIm and that f-ImPyIm preferentially binds CGCG over multiple competing sequences. The difference in the binding of f-ImPyIm and f-PyImIm to their cognate sequences was supported by the Na(+)-dependent nature of DNA melting studies, in which significantly higher Na(+) concentrations were needed to match the ability of f-ImPyIm to stabilize CGCG with that of f-PyImIm stabilizing CCGG. The selectivity of f-ImPyIm beyond the four-base CGCG recognition site was tested by circular dichroism and isothermal titration microcalorimetry, which shows that f-ImPyIm has marginal selectivity for (A.T)CGCG(A.T) over (G.C)CGCG(G.C). In addition, changes adjacent to this 6 bp binding site do not affect f-ImPyIm affinity. Calorimetric studies revealed that binding of f-ImPyIm, f-PyImIm, and distamycin A to their respective hairpin cognate sequences is exothermic; however, changes in enthalpy, entropy, and heat capacity (DeltaC(p)) contribute differently to formation of the 2:1 complexes for each triamide. Experimental and theoretical determinations of DeltaC(p) for binding of f-ImPyIm to CGCG were in good agreement (-142 and -177 cal mol(-)(1) K(-)(1), respectively). (1)H NMR of f-ImPyIm and f-PyImIm complexed with their respective cognate DNAs confirmed positively cooperative formation of distinct 2:1 complexes. The NMR results also showed that these triamides bind in the DNA minor groove and that the oligonucleotide retains the B-form conformation. Using minimal distance restraints from the NMR experiments, molecular modeling and dynamics were used to illustrate the structural complementarity between f-ImPyIm and CGCG. Collectively, the NMR and ITC experiments show that formation of the 2:1 f-ImPyIm-CGCG complex achieves a structure more ordered and more thermodynamically favored than the structure of the 2:1 f-PyImIm-CCGG complex.

Towards an understanding of DNA recognition by the methyl-CpG binding domain 1

CpG methylation determines a variety of biological functions of DNA. The methylation signal is interpreted by proteins containing a methyl-CpG binding domain (MBDs). Based on the NMR structure of MBD1 complexed with methylated DNA we analysed the recognition mode by means of molecular dynamics simulations. As the protein is monomeric and recognizes a symmetrically methylated CpG step, the recognition mode is an asymmetric one. We find that the two methyl groups do not contribute equally to the binding energy. One methyl group is associated with the major part of the binding energy and the other one nearly does not contribute at all. The contribution of the two cytosine methyl groups to binding energy is calculated to be -3.6 kcal/mol. This implies a contribution of greater than two orders of magnitude to the binding constant. The conserved amino acid Asp32 is known to be essential for DNA binding by MBD1, but so far no direct contact with DNA has been observed. We detected a direct DNA base contact to Asp32. This could be the main reason for the importance of this amino acid. MBD contacts DNA exclusively in the major groove, the minor groove is reserved for histone contacts. We found a deformation of the minor groove shape due to complexation by MBD1, which indicates an information transfer between the major and the minor groove.

Molecular Docking Studies of Natural Cholinesterase-Inhibiting Steroidal Alkaloids from Sarcococca saligna

Alkaloids isolated from Sarcococca saligna significantly inhibit acetyl- and butyrylcholinesterase enzyme, suggesting discovery of inhibitors for nervous-system disorders. Studying interactions with the active site of the AChE enzyme from Torpedo californica, we have identified hydrophobic interactions inside the aromatic gorge area as the major stabilizing factor in enzyme-inhibitor complexes of these alkaloids. Molecular Dynamics simulation of a predicted complex indicates that ligand binding does not extensively alter enzyme structure, but reduces flexibility at the gorge.

Molecular dynamics simulations and thermodynamics analysis of DNA-drug complexes. Minor groove binding between 4',6-diamidino-2-phenylindole and DNA duplexes in solution

An extended set of nanosecond-scale molecular dynamics simulations of DNA duplex sequences in explicit solvent interacting with the minor groove binding drug 4',6-diamidino-2-phenylindole (DAPI) are investigated for four different and sequence specific binding modes. Force fields for DAPI have been parametrized to properly reflect its internal nonplanarity. Sequences investigated include the binding modes observed experimentally, that is, AATT in d(CGCGAATTCGCG)(2) and ATTG in d(GGCCAATTGG)(2) and alternative shifted binding modes ATTC and AATT, respectively. In each case, stable MD simulations are obtained, well reproducing specific hydration patterns seen in the experiments. In contrast to the 2.4 A d(CGCGAATTCGCG)(2) crystal structure, the DAPI is nonplanar, consistent with its gas-phase geometry and the higher resolution crystal structure. The simulations also suggest that the DAPI molecule is able to adopt different conformational substates accompanied by specific hydration patterns that include long-residing waters. The MM_PBSA technology for estimating relative free energies was utilized. The most consistent free energy results were obtained with an approach that uses a single trajectory of the DNA-DAPI complex to estimate all free energy terms. It is demonstrated that explicit inclusion of a subset of bound water molecules shifts the calculated relative binding free energies in favor of both crystallographically observed binding modes, underlining the importance of structured hydration.

Hydration of hydroxypyrrole influences binding of ImHpPyPy-beta-Dp polyamide to DNA

Ligands which are able to recognize DNA sequence specifically are of fundamental interest as transcription controlling drugs. Recently a polyamide ligand was developed (ImHpPyPy-beta-Dp) which differentiates in a dimeric arrangement between all four possible base pair steps in the minor groove. This is a landmark for the design of DNA binding drugs because it was believed that such a recognition could only be possible in the major groove of DNA. Although the OH groups of the hydroxypyrrole (Hp) moieties of the ligands are responsible for this sequence discrimination, experiments showed that this OH group also reduces the absolute binding constant. We performed a free energy calculation by means of thermodynamic integration in order to find out the influence of this single hydroxyl on DNA binding. In our simulation, we found that the hydroxyl group reduces binding by about 1.3 kcal/mol, which is in excellent agreement with the experimentally determined value of 1.2 kcal/mol. In further MD simulations, the structural reasons for this reduction was estimated. The results of these simulations qualitatively agree with the X-ray structures, but in contrast, in the simulations both (ImHpPyPy-beta-Dp and ImPyPyPy-beta-Dp) ligand-DNA (d(CCAGTACTGG)(2)) complexes exhibit only slight structural differences. This is consistent with a recently published second pair of similar polyamide DNA crystal structures. Thus, we believe that the explanations resulting from the X-ray structures must be modified. We attribute the large structural differences between the two polyamide DNA complexes to a buffer molecule which binds only in the case of the ImHpPyPy-beta-Dp-DNA complex at the region of interest. We propose that the differential hydration of both ligands in the unbound state is responsible for the reduction of the binding constant. Additionally, we suggest an indirect readout of DNA, because of a lengthening of the Watson-Crick base pairs, which possibly contributes to the differentiation between T.A, A.T from G.C, C.G base pairs.

PvuII-endonuclease induces structural alterations at the scissile phosphate group of its cognate DNA

We investigated the PvuII endonuclease with its cognate DNA by means of molecular dynamics simulations. Comparing the complexed DNA with a reference simulation of free DNA, we saw structural changes at the scissile phosphodiester bond. At this GpC step, the enzyme induces the highest twist and axial rise, inclination is increased and the minor groove widened. The distance between the scissile phosphate group and the phosphate group of the following thymine base is shortened significantly, indicating a substrate-assisted catalysis. A feasible reason for this vicinity is the catalytically important amino acid residue lysine 70, which bridges the free oxygen atoms of the successive phosphate groups. Due to this geometry, a compact reaction pocket is formed where a water molecule can be held, thus bringing the reaction partners for hydrolysis into contact. The O1-P-O2 angle of the scissile nucleotide is decreased, probably due to a complexation of the negative oxygen atoms through protein and solvent contacts.

Influence of netropsin's charges on the minor groove width of d(CGCGAATTCGCG)2

The exact understanding of the interaction of minor groove binding drugs with DNA is of interest due to their importance as transcription controlling drugs. In this study we performed four molecular dynamics simulations, one of the uncomplexed d(CGCGAATTCGCG)(2) dodecamer and three simulations of the DNA complexed with the minor groove binder netropsin. The charged guanidinium and amidinium ends of the small ligand were in one simulation formally uncharged, in the second one normally charged, and in the third simulation we doubled the charges of the two ends. So we are able to filter out the influence the charges exert on the DNA structure. The positive charges reduce the width of the minor groove showing that charges are able to modify the groove width by charge neutralization of the negative phosphate groups. The quality of the used force field was successfully tested by comparing the results of the uncomplexed dodecamer with already reported NMR and x-ray studies. Thus our simulations should be able to describe the minor groove width of DNA in a correct manner underlying the validity of the results.

Simulations of nucleic acids and their complexes

Recent years have seen considerable progress in simulations of nucleic acids. Improvements in force fields, simulation techniques and protocols, and increasing computer power have all contributed to making nanosecond-scale simulations of both DNA and RNA commonplace. The results are already helping to explain how nucleic acids respond to their environment and to their base sequence and to reveal the factors underlying recognition processes by probing biologically important nucleic acid-protein interactions and medically important nucleic acid-drug complexation. This Account summarizes methodological progress and applications of molecular dynamics to nucleic acids over the past few years and tries to identify remaining challenges.

Significance of ligand tails for interaction with the minor groove of B-DNA

Minor groove binding ligands are of great interest due to their extraordinary importance as transcription controlling drugs. We performed three molecular dynamics simulations of the unbound d(CGCGAATTCGCG)(2) dodecamer and its complexes with Hoechst33258 and Netropsin. The structural behavior of the piperazine tail of Hoechst33258, which has already been shown to be a contributor in sequence-specific recognition, was analyzed. The simulations also reveal that the tails of the ligands are able to influence the width of the minor groove. The groove width is even sensitive for conformational transitions of these tails, indicating a high adaptability of the minor groove. Furthermore, the ligands also exert an influence on the B(I)/B(II) backbone conformational substate behavior. All together these results are important for the understanding of the binding process of sequence-specific ligands.

Complex of B-DNA with polyamides freezes DNA backbone flexibility

The development of sequence-specific minor groove binding ligands is a modern and rapidly growing field of research because of their extraordinary importance as transcription-controlling drugs. We performed three molecular dynamics simulations in order to clarify the influence of minor groove binding of two ImHpPyPy-beta-Dp polyamides to the d(CCAGTACTGG)(2) decamer in the B-form. This decamer contains the recognition sequence for the trp repressor (5'-GTACT-3'), and it was investigated recently by X-ray crystallography. On one hand we are able to reproduce X-ray-determined DNA--drug contacts, and on the other hand we provide new contact information which is important for the development of potential ligands. The new insights show how the beta-tail of the polyamide ligands contributes to binding. Our simulations also indicate that complexation freezes the DNA backbone in a specific B(I) or B(II) substate conformation and thus optimizes nonbonded contacts. The existence of this distinct B(I)/B(II) substate pattern also allows the formation of water-mediated contacts. Thus, we suggest the B(I) <==> B(II) substate behavior to be an important part of the indirect readout of DNA.