Stereochemical aspects of interaction of DNA binding domain of human progesterone receptor with d(AGGTCATGCT)2.

Structures of (i) 66 amino-acid fragment (residues 567-633) from DNA binding domain of human progesterone receptor (hPR), (ii) a ten base pair DNA sequence d(AGGTCATGCT)2 from hormone responsive element (HRE) and (iii) a complex of these two are optimised by computer modelling and molecular mechanics technique using extensive steric constraints from secondary structure predictions, comparison with the structures of known metalloproteins, geometric constraints imposed by tetrahedral coordination with the zinc ion and comparison with structures of DNA binding domains of human glucocorticoid and estrogen receptors (hGR and hER). Structure of the complex was obtained using genetic modification data on steroid receptors and general consensus about protein-DNA interaction. DNA is in distorted B conformation. Sequence dependent as well as protein-induced conformation changes are noticed. There is change in propeller twist, buckle and angle between glycosyl bonds. However, H-bonding network is preserved. The complex is stabilized with eighteen hydrogen-bonds, mainly between peptide side-chains and backbone phosphate. There are five specific H-bonds between basic amino acid side chains, Lys 22, Lys 26 and Arg 27, and DNA bases, A1, G3, G16 and A17. Gly 19, Ser 20 and Val 23 are in close proximity of DNA.     

Molecular mechanics simulation of the interaction of estrogen receptor with estrogen regulatory element.

A model for the interaction of 31 amino acid fragment (protein) from DNA binding domain of human estrogen receptor (hER) with a five base pair DNA sequence 5'GGTCA 3' from estrogen regulatory element (ERE) has been obtained using a step-wise procedure based on structural data on model peptides, DNA binding domain of hER, steric constrains imposed by tetrahedral coordination of the Cys sulphurs with zinc ion and classical secondary structural elements. Structure of the protein as well as its complex with DNA is obtained by energy minimization followed by refinement by molecular mechanics. The complex is stabilized by H-bonds between Lys22, Lys26 and Arg27 with DNA bases G2, T3 and T6. Lys22 also made H-bond with the backbone of G2. The backbone of Cys18 H-bonded with N7 of G1. DNA was in distorted B form and showed evidence of protein-induced conformational changes.     

A study of the interaction of Cr(III) complexes and their selective binding with B-DNA: a molecular modeling approach

Molecular modeling and energy minimisation calculations have been used to investigate the interaction of chromium(III) complexes in different ligand environments with various sequences of B-DNA. The complexes are [Cr(salen)(H(2)O)(2)](+); salen denotes 1, 2 bis-salicylideneaminoethane, [Cr(salprn)(H(2)O)(2)](+); salprn denotes 1, 3 bis- salicylideneaminopropane, [Cr(phen)(3)](3+); phen denotes 1, 10 phenanthroline and [Cr(en)(3)](3+); en denotes ethylenediamine. All the chromium(III) complexes are interacted with the minor groove and major groove of d(AT)(12), d(CGCGAATTCGCG)(2) and d(GC)(12) sequences of DNA. The binding energy and hydrogen bond parameters of DNA-Cr complex adduct in both the groove have been determined using molecular mechanics approach. The binding energy and formation of hydrogen bonds between chromium(III) complex and DNA has shown that all complexes of chromium(III) prefer minor groove interaction as the favourable binding mode.     

Molecular modelling of the interaction of carbocyclic analogues of netropsin and distamycin with d(CGCGAATTCGCG)2

A molecular mechanics and molecular dynamics approach was used to examine the structure of complexes formed between the d(CGCGAATTCGCG)2 duplex and netropsin, distamycin, and four carbocyclic analogues of netropsin and distamycin (1-4). The resulting structures of the ligand-DNA model complexes and their energetics were examined. It is predicted that the compounds 1-4 should have a decreased affinity for the minor groove of AT-rich regions in comparison to netropsin and distamycin. From the energetic analysis it appears that van der Waals and electrostatic interactions are more important than specific hydrogen bonds in stabilizing the ligand-duplex complexes. We predict that compounds 1 and 2 are effectively isohelical with the DNA minor groove. The superior DNA-binding afforded by 1 and 2 in comparison to 3 and 4 results from their more effective penetration into the minor groove and smaller perturbation of molecular structure upon complex formation.     

A Study of the Interaction of Cr(III) Complexes and Their Selective Binding with B-DNA: A Molecular Modeling Approach

Abstract

Molecular modeling and energy minimisation calculations have been used to investigate the interaction of chromium(III) complexes in different ligand environments with various sequences of B-DNA. The complexes are [Cr(salen)(H₂O)₂]⁺; salen denotes 1, 2 bis-salicylideneaminoethane, [Cr(salprn)(H₂O)₂]⁺; salprn denotes 1, 3 bis- salicylideneamino-propane, [Cr(phen)₃]³⁺; phen denotes 1, 10 phenanthroline and [Cr(en)₃]³⁺; en denotes eth- ylenediamine. All the chromium(III) complexes are interacted with the minor groove and major groove of d(AT)₁₂, d(CGCGAATTCGCG)₂ and d(GC)₁₂ sequences of DNA. The binding energy and hydrogen bond parameters of DNA-Cr complex adduct in both the groove have been determined using molecular mechanics approach. The binding energy and formation of hydrogen bonds between chromium(III) complex and DNA has shown that all complexes of chromium(III) prefer minor groove interaction as the favourable binding mode.     

Minor groove deformability of DNA: a molecular dynamics free energy simulation study

The conformational deformability of nucleic acids can influence their function and recognition by proteins. A class of DNA binding proteins including the TATA box binding protein binds to the DNA minor groove, resulting in an opening of the minor groove and DNA bending toward the major groove. Explicit solvent molecular dynamics simulations in combination with the umbrella sampling approach have been performed to investigate the molecular mechanism of DNA minor groove deformations and the indirect energetic contribution to protein binding. As a reaction coordinate, the distance between backbone segments on opposite strands was used. The resulting deformed structures showed close agreement with experimental DNA structures in complex with minor groove-binding proteins. The calculated free energy of minor groove deformation was approximately 4-6 kcal mol(-1) in the case of a central TATATA sequence. A smaller equilibrium minor groove width and more restricted minor groove mobility was found for the central AAATTT and also a significantly ( approximately 2 times) larger free energy change for opening the minor groove. The helical parameter analysis of trajectories indicates that an easier partial unstacking of a central TA versus AT basepair step is a likely reason for the larger groove flexibility of the central TATATA case.     

Structure of the Wilms tumor suppressor protein zinc finger domain bound to DNA

The zinc finger domain of the Wilms tumor suppressor protein (WT1) contains four canonical Cys(2)His(2) zinc fingers. WT1 binds preferentially to DNA sequences that are closely related to the EGR-1 consensus site. We report the structure determination by both X-ray crystallography and NMR spectroscopy of the WT1 zinc finger domain in complex with DNA. The X-ray structure was determined for the complex with a cognate 14 base-pair oligonucleotide, and composite X-ray/NMR structures were determined for complexes with both the 14 base-pair and an extended 17 base-pair DNA. This combined approach allowed unambiguous determination of the position of the first zinc finger, which is influenced by lattice contacts in the crystal structure. The crystal structure shows the second, third and fourth zinc finger domains inserted deep into the major groove of the DNA where they make base-specific interactions. The DNA duplex is distorted in the vicinity of the first zinc finger, with a cytidine twisted and tilted out of the base stack to pack against finger 1 and the tip of finger 2. By contrast, the composite X-ray/NMR structures show that finger 1 continues to follow the major groove in the solution complexes. However, the orientation of the helix is non-canonical, and the fingertip and the N terminus of the helix project out of the major groove; as a consequence, the zinc finger side-chains that are commonly involved in base recognition make no contact with the DNA. We conclude that finger 1 helps to anchor WT1 to the DNA by amplifying the binding affinity although it does not contribute significantly to binding specificity. The structures provide molecular level insights into the potential consequences of mutations in zinc fingers 2 and 3 that are associated with Denys-Drash syndrome and nephritic syndrome. The mutations are of two types, and either destabilize the zinc finger structure or replace key base contact residues.     

Zinc finger proteins: getting a grip on RNA

C2H2 (Cys-Cys-His-His motif) zinc finger proteins are members of a large superfamily of nucleic-acid-binding proteins in eukaryotes. On the basis of NMR and X-ray structures, we know that DNA sequence recognition involves a short alpha helix bound to the major groove. Exactly how some zinc finger proteins bind to double-stranded RNA has been a complete mystery for over two decades. This has been resolved by the long-awaited crystal structure of part of the TFIIIA-5S RNA complex. A comparison can be made with identical fingers in a TFIIIA-DNA structure. Additionally, the NMR structure of TIS11d bound to an AU-rich element reveals the molecular details of the interaction between CCCH fingers and single-stranded RNA. Together, these results contrast the different ways that zinc finger proteins bind with high specificity to their RNA targets.     

Induction of DNA bending by bifunctional intercalating agents of the 7H-pyridocarbazole family

The structures and binding energetics of selected complexes formed between the deoxynucleotides d(CpGpGpCpG).d(CpGpCpCpG), d(CpGpApTpCpG)2, d(GpCpGpCpCpG).d(CpGpGpCpGpC), and d(CpGpCpCpCpG)2 with the DNA bifunctional intercalating agent ditercalinium and three of its rigid linking chain derivatives have been investigated theoretically by means of a molecular mechanics approach that takes into account nucleic acid flexibility, ligand flexibility and solvent dielectric effects (R. Lavery, in: Unusual DNA structures, eds S. Harvey and R. Wells (Pergamon, New York, 1988) p. 189; R. Lavery, in: DNA bending and curvature, eds W.K. Olson et al. (Adenine Press, New York, 1988) p. 191). The piperidinium chains of the bis-intercalating ligands are always located in the major groove of DNA. For the energy-minimized complexes the ligand proceeds to bind following preferentially the 5'-pyrimidine-purine-3' alternating sequence, thus dictating the number of internal exclusion sites. The complexes with three exclusion sites will present (i) a bending of the structure towards the major groove, and (ii) a non-ideal distribution of unwinding angles; complexes with less than three exclusion sites will remain essentially linear. The absence of a bend does not preclude other types of local deformations of the base-pairs such as inclination, buckle and tip. The proposed structures of the d(CpGpApTpCpG)2 complexes are in agreement with NMR structural results. The possible relevance of these findings to a previously proposed mode of interaction for ditercalinium-like DNA ligands is discussed.     

Structural basis for the binding affinity of a homologous series of synthetic phenoxazone drugs with DNA: NMR and molecular mechanics analysis

The molecular basis of the marked structure-activity relationship for a homologous series of DNA-binding phenoxazone drugs (ActII-ActIV) has been investigated by NMR spectroscopy and molecular mechanics. The spatial structures of the complexes between the drugs and a model deoxytetranucleotide, 5'-d(TpGpCpA), have been determined by molecular mechanics methods using homonuclear (1)H-(1)H 2D-NOESY and heteronuclear (1)H-(31)P (HMBC) NMR spectroscopic data. Observed intermolecular NOE contacts and equilibrium binding studies confirm that the binding affinity of the synthetic phenoxazone derivatives with d(TGCA) decreases with an increase in the number of CH(2) groups in the dimethylaminoalkyl side chains, i.e., ActII > ActIII > ActIV, in agreement with the observed biological activity of these compounds. Molecular mechanics calculations of the spatial structures of the intercalated complexes of ActII-ActIV with d(TGCA) indicate that the different binding constants of the phenoxazone derivatives with the DNA oligomer are due to the different degrees of intercalation of the chromophore and the different steric arrangements of aminoalkyl side chains in the minor groove of the tetramer duplex; this results in different distances between the negatively-charged phosphates of the DNA duplex and the terminal positively-charged N(CH(3))(2) groups of the side chains.     

Partial B-to-A DNA transition upon minor groove binding of protein Sac7d monitored by Raman spectroscopy

Members of the Sso7d/Sac7d protein family and other related proteins are believed to play an important role in DNA packaging and maintenance in archeons. Sso7d/Sac7d are small, abundant, basic, and nonspecific DNA-binding proteins of the hyperthermophilic archeon Sulfolobus. Structures of several complexes of Sso7d/Sac7d with DNA octamers are known. These structures are characterized by sequence unspecific minor groove binding of the proteins and sharp kinking of the double helix. Corresponding Raman vibrational signatures have been identified in this study. A Raman spectroscopic analysis of Sac7d binding to the oligonucleotide decamer d(GAGGCGCCTC)(2) reveals large conformational perturbations in the DNA structure upon complex formation. Perturbed Raman bands are associated with the vibrational modes of the sugar phosphate backbone and frequency shifts of bands assigned to nucleoside vibrations. Large changes in the DNA backbone and partial B- to A-form DNA transitions are indicated that are closely associated with C2'-endo/anti to C3'-endo/anti conversion of the deoxyadenosyl moiety upon Sac7d binding. The major spectral feature of Sac7d binding is kinking of the DNA. Raman markers of minor groove binding do not largely contribute to spectral differences; however, clear indications for minor groove binding come from G-N2 and G-N3 signals that are supported by Trp24 features. Trp24 is the only tryptophan present in Sac7d and binds to guanine N3, as has been demonstrated clearly in X-ray structures of Sac7d-DNA complexes. No changes of the Sac7d secondary structure have been detected upon DNA binding.     

Estrogen receptor impairs interleukin-6 expression by preventing protein binding on the NF-kappaB site.

Interleukin-6 (IL-6) is a multifunctional cytokine thought to be a key factor in post-menopausal osteoporosis, given its ability to induce osteoclast maturation and its down regulation by estrogens. We have previously shown that the effects of TNFalphaand estradiol on the human IL-6 promoter were dependent on a region of the promoter containing a C/EBP site and a NF-kappaB site. To define the molecular mode of action of estrogens, we performed gel shift assays with this DNA fragment as a probe, and nuclear extracts from TNFalpha-induced HeLa, MCF7 and Saos2 cells. Several induced complexes specifically bound the probe. The use of various competitor DNA suggested that most of the complexes detected contained NF-kappaB factors, and that C/EBP site binding factors were important for the overall binding to the probe. Addition of in vitro translated human estrogen receptor (hER) impaired the binding of three complexes in HeLa cells and two complexes in MCF7 and Saos2 cells. Competition experiments suggested that the NF-kappaB site was necessary for the effect of hER. The use of antisera against NF-kappaB and C/EBP proteins showed that the target complexes of hER contained the c-rel proto-oncogene product and to a lesser extent, the RelA protein. Taken together, these data show that hER impairs TNFalphainduction of IL-6 by preventing c-rel and, to a lesser extent, RelA proteins binding to the NF-kappaB site of the IL-6 promoter.     

DNA-induced alpha-helix capping in conserved linker sequences is a determinant of binding affinity in Cys(2)-His(2) zinc fingers

High-affinity, sequence-specific DNA binding by Cys(2)-His(2) zinc finger proteins is mediated by both specific protein-base interactions and non-specific contacts between charged side-chains and the phosphate backbone. In addition, in DNA complexes of multiple zinc fingers, protein-protein interactions between the finger units contribute to the binding affinity. We present NMR evidence for another contribution to high- affinity binding, a highly specific DNA-induced helix capping involving residues in the linker sequence between fingers. Capping at the C terminus of the alpha-helix in each zinc finger, incorporating a consensus TGEKP linker sequence that follows each finger, provides substantial binding energy to the DNA complexes of zinc fingers 1-3 of TFIIIA (zf1-3) and the four zinc fingers of the Wilms' tumor suppressor protein (wt1-4). The same alpha-helix C-capping motif is observed in the X-ray structures of four other protein-DNA complexes. The structures of each of the TGEKP linkers in these complexes can be superimposed on the linker sequences in the zf1-3 complex, revealing a remarkable similarity in both backbone and side-chain conformations. The canonical linker structures from the zinc-finger-DNA complexes have been compared to the NMR structure of the TGEKP linker connecting fingers 1 and 2 in zf1-3 in the absence of DNA. This comparison reveals that additional stabilization likely arises in the DNA complexes from hydrogen bonding between the backbone amide of E3 and the side-chain O(gamma) of T1 in the linker. We suggest that these DNA-induced C-capping interactions provide a means whereby the multiple-finger complex, which must necessarily be domain-flexible in the unbound state as it searches for the correct DNA sequence, can be "snap-locked" in place once the correct DNA sequence is encountered. These observations provide a rationale for the high conservation of the TGEKP linker sequences in Cys(2)-His(2) zinc finger proteins.     

A thermodynamic and structural analysis of DNA minor-groove complex formation

As part of an effort to develop a better understanding of the structural and thermodynamic principles of DNA minor groove recognition, we have investigated complexes of three diphenylfuran dications with the d(CGCGAATTCGCG)(2) duplex. The parent compound, furamidine (DB75), has two amidine substituents while DB244 has cyclopentyl amidine substituents and DB226 has 3-pentyl amidines. The structure for the DB244-DNA complex is reported here and is compared to the structure of the DB75 complex. Crystals were not obtained with DB226 but information from the DB75 and DB244 structures as well as previous NMR results on DB226 indicate that all three compounds bind in the minor groove at the AATT site of the duplex. DB244 and DB75 penetrate to the floor of the groove and form hydrogen bonds with T8 on one strand and T20 on the opposite strand while DB226 forms a complex with fewer interactions. Binding studies by surface plasmon resonance (SPR) yield -delta G degrees values in the order DB244>DB75>DB226 that are relatively constant with temperature. The equilibrium binding constants for DB244 are 10-20 times greater than that for DB226. Isothermal titration calorimetric (ITC) experiments indicate that, in contrast to delta G degrees, delta H degrees varies considerably with temperature to yield large negative delta Cp degrees values. The thermodynamic results, analyzed in terms of structures of the DNA complexes, provide an explanation of why DB244 binds more strongly to DNA than DB75, while DB266 binds more weakly. All three compounds have a major contribution to binding from hydrophobic interactions but the hydrophobic term is most favorable for DB244. DB244 also has strong contributions from molecular interactions in its DNA complex and all of these factors combine to give it the largest-delta G degrees for binding. Although the factors that influence the energetics of minor groove interactions are varied and complex, results from the literature coupled with those on the furan derivatives indicate that there are some common characteristics for minor groove recognition by unfused heterocyclic cations that can be used in molecular design.     

Molecular details of quinolone-DNA interactions: solution structure of an unusually stable DNA duplex with covalently linked nalidixic acid residues and non-covalent complexes derived from it

Quinolones are antibacterial drugs that are thought to bind preferentially to disturbed regions of DNA. They do not fall into the classical categories of intercalators, groove binders or electrostatic binders to the backbone. We solved the 3D structure of the DNA duplex (ACGCGU-NA)2, where NA denotes a nalidixic acid residue covalently linked to the 2'-position of 2'-amino-2'-deoxyuridine, by NMR and restrained torsion angle molecular dynamics (MD). In the complex, the quinolones stack on G:C base pairs of the core tetramer and disrupt the terminal A:U base pair. The displaced dA residues can stack on the quinolones, while the uracil rings bind in the minor groove. The duplex-bridging interactions of the drugs and the contacts of the displaced nucleotides explain the high UV-melting temperature for d(ACGCGU-NA)2 of up to 53 degrees C. Further, non-covalently linked complexes between quinolones and DNA of the sequence ACGCGT can be generated via MD using constraints obtained for d(ACGCGU-NA)2. This is demonstrated for unconjugated nalidixic acid and its 6-fluoro derivative. The well-ordered and tightly packed structures thus obtained are compatible with a published model for the quinolone-DNA complex in the active site of gyrases.     

Structural Basis for the Binding Affinity of a Homologous Series of Synthetic Phenoxazone Drugs with DNA: NMR and Molecular Mechanics Analysis

Abstract

The molecular basis of the marked structure-activity relationship for a homologous series of DNA-binding phenoxazone drugs (ActII-ActIV) has been investigated by NMR spectroscopy and molecular mechanics. The spatial structures of the complexes between the drugs and a model deoxytetranucleotide, 5′-d(TpGpCpA), have been determined by molecular mechanics methods using homonuclear ¹H-¹H 2D-NOESY and heteronuclear ¹H-³¹P (HMBC) NMR spectroscopic data. Observed intermolecular NOE contacts and equilibrium binding studies confirm that the binding affinity of the synthetic phenoxazone derivatives with d(TGCA) decreases with an increase in the number of CH₂ groups in the dimethylami- noalkyl side chains, i.e., ActII > ActIII > ActIV, in agreement with the observed biological activity of these compounds. Molecular mechanics calculations of the spatial structures of the intercalated complexes of ActII-ActIV with d(TGCA) indicate that the different binding constants of the phenoxazone derivatives with the DNA oligomer are due to the different degrees of intercalation of the chromophore and the different steric arrangements of aminoalkyl side chains in the minor groove of the tetramer duplex; this results in different distances between the negatively-charged phosphates of the DNA duplex and the terminal positively-charged N(CH₃)₂ groups of the side chains.     

Flexible docking of DNA fragments and actinocin derivatives

We have applied molecular docking methods to systems containing nucleic acids as targets and biologically active substances as ligands. The complexes of DNA fragments and actinocin derivatives with different lengths of aminoalkyl side chains were obtained by molecular docking. It was observed that actinocin derivatives could form energetically favourable complexes with DNA both as intercalators and minor groove binders. It was shown that small changes in the binding energy (～1?kcal/mol) could result in complexes with substantially different structure. The complexes of actinocin derivatives and DNA fragments were stabilized by hydrogen bonding upon intercalation and minor groove binding. It was found that the change of solvent-accessible surface area upon binding of the actinocin derivative to DNA linear increased with the growth of methylene groups' number in ligand side chains. The solvation energy change upon binding of actinocin derivatives to DNA calculated by the WSAS method was favourable in the case of small uncharged ligands and unfavourable for positively charged ligands.     

Solution structure of [d(A-T)5]2 via complete relaxation matrix analysis of two-dimensional nuclear Overhauser effect spectra and molecular mechanics calculations: evidence for a hydration tunnel

Pure absorption phase proton two-dimensional nuclear Overhauser effect (2D NOE) spectra at 500 MHz have been obtained for [d(5'ATATATATAT3')]2 in deuterium oxide solution at several mixing times. The 100 nonexchangeable proton resonances have been assigned. The experimental 2D NOE spectra were compared with theoretical spectra calculated by using the complete relaxation matrix analysis method [Keepers, J. W., & James, T. L. (1984) J. Magn. Reson. 57, 404-426] and x-ray diffraction determined molecular coordinates of A, B, alternating B, left-handed B, C, D, and wrinkled D forms of DNA and of energy-minimized structures calculated from the most promising X-ray crystal structures by using the molecular mechanics program AMBER, in which all hydrogens, counterions, and hydration water molecules were included. The analysis of all features of the 2D NOE spectra played an important role in extracting the promising structures, and it was concluded that the wrinkled D form yields the best fit for the 2D NOE data of the A-T decamer. The molecular mechanics calculation indicated that this model structure, whose minor groove is comparatively deep and narrow, may be energetically more stable than the B form for alternating d(A-T) DNA. Interesting features of the structure include possible intra- and interchain sugar-phosphate attractions and a hydration tunnel inside the minor groove capable of accommodating three types of water molecules that aid in helix stabilization via hydrogen bonding. Counterions (sodium) serve to reduce interchain phosphate-phosphate repulsive effects.     

Molecular dynamics of spermine-DNA interactions: sequence specificity and DNA bending for a simple ligand.

We used molecular dynamics to model interactions between the physiologically important polyamine spermine and two B-DNA oligomers, the homopolymer (dG)10-(dC)10 and the heteropolymer (dGdC)5-(dGdC)5. Water and counterions were included in the simulation. Starting coordinates for spermine-DNA complexes were structures obtained by molecular mechanics modeling of spermine with the two oligomers; in these models, spermine binding induced a bend in the heteropolymer but not in the homopolymer. During approximately 40 psec of molecular dynamics simulation, spermine moves away from the floor of the major groove and interacts nospecifically with d(G)10-d(C)10. In contrast, a spermine-induced bend in the helix of (dGdC)5-(dGdC)5 is maintained throughout the simulation and spermine remains closely associated with the major groove. These results provide further evidence that the binding of spermine to nucleic acids can be sequence specific and that bending of alternating purine-pyrimidine sequences may be a physiologically important result of spermine binding.