Crystallographic insights into DNA minor groove recognition by drugs

This review surveys the crystal structures between minor groove drugs and oligonucleotides, of which over thirty have now been determined. The various factors that are involved in the observed A/T sequence selectivity of these drugs are examined in structural terms. The roles of, in particular, hydrogen-bond recognition and sequence-dependent groove width, are assessed, and as a consequence the minor groove drugs have been classified into two categories, dependent on the relative roles played by these two factors in sequence recognition. Implications for the recognition of non-A/T sequences are discussed. © 1997 John Wiley & Sons, Inc. Biopoly 44: 105–121, 1997     

DNA recognition by lexitropsins,minor groove binding agents

Consideration is given to alternative approaches to the development of DNA sequences selective binding agents because of their potential applications in diagnosis and treatment of cancer as well as in molecular biology. The concept of lexitropsins, or information-reading molecules, is introduced within the antigene strategy as an alternative to, and complementary with, the antigene approach for cellular intervention and gene control The chemical, physical and paharmacological factors involved in the design of effective lexitropsins are discussed and illustrated with experimental results. Among the factors contributing to the molecular recognition processes are: the presence and disposition of hydrogen bond accepting and donating groups, ligand shape, chirality, stereochemistry, flexibility and charge. For longer ligands, such as are required to target unique sequences in biological systems (14–16 base pairs), the critical feature is the phasing or spatial corresponding between repeat units in the ligand and receptor. The recently discovered 2:1 lexitropsin-DNA binding motif provides a further refinement in molecular recognition in permitting discrimination between GC and CG base pairs. The application of these factors in the design and synthesis of novel agents which exhibits anticancer, antiviral and antitretroviral properties, and inhibition of critical cellular enzymes including topoisomerases is discussed. The emerging evidence of a relationship between sequence selectivity of the new agents and the biological responses they invoke is also described.     

DNA recognition by the EcoRV restriction endonuclease probed using base analogues

The EcoRV restriction endonuclease recognises palindromic GATATC sequences and cuts between the central T and dA bases in a reaction that has an absolute requirement for a divalent metal ion, physiologically Mg(2+). Use has been made of base analogues, which delete hydrogen bonds between the protein and DNA (or hydrophobic interactions in the case of the 5-CH(3) group of thymine), to evaluate the roles of the outer two base-pairs (GATATC) in DNA recognition. Selectivity arises at both the binding steps leading to the formation of the enzyme-DNA-metal ion ternary complex (assayed by measuring the dissociation constant in the presence of the non-reactive metal Ca(2+)) and the catalytic step (evaluated using single-turnover hydrolysis in the presence of Mg(2+)), with each protein-DNA contact contributing to recognition. With the A:T base-pair, binding was reduced by the amount expected for the simple loss of a single contact; much more severe effects were observed with the G:C base-pair, suggesting additional conformational perturbation. Most of the modified bases lowered the rate of hydrolysis; furthermore, the presence of an analogue in one strand of the duplex diminished cutting at the second, unmodified strand, indicative of communication between DNA binding and the active site. The essential metal ion Mg(2+) plays a key role in mediating interactions between the DNA binding site and active centre and in many instances rescue of hydrolysis was seen with Mn(2+). It is suggested that contacts between the GATATC site are required for tight binding and for the correct assembly of metal ions and bound water at the catalytic site, functions important in providing acid/base catalysis and transition state stabilisation.     

Characterisation of the hydrophobic collapse of polystyrene in water using free energy techniques

We characterise the hydrophobic collapse of single polystyrene chains in water using molecular dynamics simulations. Specifically, we calculate the potential of mean force for the collapse of a single polystyrene chain in water using metadynamics, comparing the results between all atomistic with coarse-grained (CG) molecular simulation. We next explore the scaling behaviour of the collapsed globular shape at the minimum energy configuration, characterised by the radius of gyration, as a function of chain length. The exponent is close to one third, consistent with that predicted for a polymer chain in bad solvent. We also explore the scaling behaviour of the solvent accessible surface area (SASA) as a function of chain length, finding a similar exponent for both all atomistic and CG simulations. Furthermore, calculation of the local water density as a function of chain length near the minimum energy configuration suggests that intermediate chain lengths are more likely to form dewetted states, as compared to shorter or longer chain lengths.     

DNA sequence recognition in the minor groove by hairpin microgonotropens

Two novel microgonotropens (MGTs) comprised of hairpin N-propylaminepyrrole polyamides linked to a Hoechst 33258 (Ht) analogue (3 and 4) were synthesized on solid phase by adopting an Fmoc technique using a series of HOBt mediated coupling reactions. The dsDNA-binding properties of MGTs 3 and 4 were determined by thermal denaturation experiments. Both MGTs were found to be selective for their nine-bp match dsDNA sequence 9 and were less tolerant of G/C bp substitutions in the binding region than linear progenitor MGT 1. MGT 3 was intolerant of a G/C substitution located in the middle of the binding region and did not bind to sequences 13 and 14. MGT 4 also did not bind to sequence 13, and its linker-bound Ht moiety was found to be more sensitive to a G/C substitution in the Ht-binding target, as demonstrated by the lack of binding to sequence 16.     

Cooperative stabilization of Zn2+:DNA complexes through netropsin binding in the minor groove of FdU-substituted DNA

The simultaneous binding of netropsin in the minor groove and Zn²⁺ in the major groove of a DNA hairpin that includes 10 consecutive FdU nucleotides at the 3′-terminus (3′FdU) was demonstrated based upon NMR spectroscopy, circular dichroism (CD), and computational modeling studies. The resulting Zn²⁺/netropsin: 3′FdU complex had very high thermal stability with aspects of the complex intact at 85?°C, conditions that result in complete dissociation of Mg²⁺ complexes. CD and ¹⁹F NMR spectroscopy were consistent with Zn²⁺ binding in the major groove of the DNA duplex and utilizing F5 and O4 of consecutive FdU nucleotides as ligands with FdU nucleotides hemi-deprotonated in the complex. Netropsin is bound in the minor groove of the DNA duplex based upon 2D NOESY data demonstrating contacts between AH2 ¹H and netropsin ¹H resonances. The Zn²⁺/netropsin: 3′FdU complex displayed increased cytotoxicity towards PC3 prostate cancer (PCa) cells relative to the constituent components or separate complexes (e.g. Zn²⁺:3′FdU) indicating that this new structural motif may be therapeutically useful for PCa treatment.

An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:32     

Thermodynamic and kinetic basis for the relaxed DNA sequence specificity of "promiscuous" mutant EcoRI endonucleases

Promiscuous mutant EcoRI endonucleases produce lethal to sublethal effects because they cleave Escherichia coli DNA despite the presence of the EcoRI methylase. Three promiscuous mutant forms, Ala138Thr, Glu192Lys and His114Tyr, have been characterized with respect to their binding affinities and first-order cleavage rate constants towards the three classes of DNA sites: specific, miscognate (EcoRI*) and non-specific. We have made the unanticipated and counterintuitive observations that the mutant restriction endonucleases that exhibit relaxed specificity in vivo nevertheless bind more tightly than the wild-type enzyme to the specific recognition sequence in vitro, and show even greater preference for binding to the cognate GAATTC site over miscognate sites. Binding preference for EcoRI* over non-specific DNA is also improved. The first-order cleavage rate constants of the mutant enzymes are normal for the cognate site GAATTC, but are greater than those of the wild-type enzyme at EcoRI* sites. Thus, the mutant enzymes use two mechanisms to partially bypass the multiple fail-safe mechanisms that protect against cleavage of genomic DNA in cells carrying the wild-type EcoRI restriction-modification system: (a) binding to EcoRI* sites is more probable than for wild-type enzyme because non-specific DNA is less effective as a competitive inhibitor; (b) the combination of increased affinity and elevated cleavage rate constants at EcoRI* sites makes double-strand cleavage of these sites a more probable outcome than it is for the wild-type enzyme. Semi-quantitative estimates of rates of EcoRI* site cleavage in vivo, predicted using the binding and cleavage constants measured in vitro, are in accord with the observed lethal phenotypes associated with the three mutations.     

A contact energy function considering residue hydrophobic environment and its application in protein fold recognition

The three-dimensional （3D） structure prediction of proteins ：is an important task in bioinformatics. Finding energy functions that can better represent residue-residue and residue-solvent interactions is a crucial way to improve the prediction accu- racy. The widely used contact energy functions mostly only consider the contact frequency between different types of residues; however, we find that the contact frequency also relates to the residue hydrophobic environment. Accordingly, we present an improved contact energy function to integrate the two factors, which can reflect the influence of hydrophobic interaction on the stabilization of protein 3D structure more effectively. Furthermore, a fold recognition （threading） approach based on this energy function is developed. The testing results obtained with 20 randomly selected proteins demonstrate that, compared with common contact energy functions, the proposed energy function can improve the accuracy of the fold template prediction from 20% to 50%, and can also improve the accuracy of the sequence-template alignment from 35% to 65%.     

DNA minor groove electrostatic potential: influence of sequence-specific transitions of the torsion angle gamma and deoxyribose conformations

The structural adjustments of the sugar-phosphate DNA backbone (switching of the γ angle (O5′–C5′–C4′–C3′) from canonical to alternative conformations and/or C2′-endo → C3′-endo transition of deoxyribose) lead to the sequence-specific changes in accessible surface area of both polar and non-polar atoms of the grooves and the polar/hydrophobic profile of the latter ones. The distribution of the minor groove electrostatic potential is likely to be changing as a result of such conformational rearrangements in sugar-phosphate DNA backbone. Our analysis of the crystal structures of the short free DNA fragments and calculation of their electrostatic potentials allowed us to determine: (1) the number of classical and alternative γ angle conformations in the free B-DNA; (2) changes in the minor groove electrostatic potential, depending on the conformation of the sugar-phosphate DNA backbone; (3) the effect of the DNA sequence on the minor groove electrostatic potential. We have demonstrated that the structural adjustments of the DNA double helix (the conformations of the sugar-phosphate backbone and the minor groove dimensions) induce changes in the distribution of the minor groove electrostatic potential and are sequence-specific. Therefore, these features of the minor groove sizes and distribution of minor groove electrostatic potential can be used as a signal for recognition of the target DNA sequence by protein in the implementation of the indirect readout mechanism.     

Exploring the role of chirality in nucleic acid recognition

The study of the base-pairing properties of nucleic acids with sugar moieties in the backbone belonging to the L-series (β-L-DNA, β-L-RNA, and their analogs) are reviewed. The major structural factors underlying the formation of stable heterochiral complexes obtained by incorporation of modified nucleotides into natural duplexes, or by hybridization between homochiral strands of opposite sense of chirality are highlighted. In addition, the perspective use of L-nucleic acids as candidates for various therapeutic applications, or as tools for both synthetic biology and etiology-oriented investigations on the structure and stereochemistry of natural nucleic acids is discussed.     

Application of electrospray ionization mass spectrometry to study the hydrophobic interaction between the epsilon and theta subunits of DNA polymerase III

The interactions between the N-terminal domain of the epsilon (epsilon186) and theta subunits of DNA polymerase III of Escherichia coli were investigated using electrospray ionization mass spectrometry. The epsilon186-theta complex was stable in 9 M ammonium actetate (pH 8), suggesting that hydrophobic interactions have a predominant contribution to the stability of the complex. Addition of primary alkanols to epsilon186-theta in 0.1 M ammonium acetate (pH 8), led to dissociation of the complex, as observed in the mass spectrometer. The concentrations of methanol, ethanol, and 1-propanol required to dissociate 50% of the complex were 8.9 M, 4.8 M, and 1.7 M, respectively. Closer scrutiny of the effect of alkanols on epsilon186, theta, and epsilon186-theta showed that epsilon186 formed soluble aggregates prior to precipitation, and that the association of epsilon186 with theta stabilized epsilon186. In-source collision-induced dissociation experiments and other results suggested that the epsilon186-theta complex dissociated in the mass spectrometer, and that the stability (with respect to dissociation) of the complex in vacuo was dependent on the solution from which it was sampled.     

DNA sequence recognition in the minor groove by hairpin pyrrole polyamide-Hoechst 33258 analogue conjugate

A hairpin pyrrole polyamide conjugated to a Hoechst 33258 (Ht) analogue, PyPyPy-gamma-PyPyPy-gamma-Ht, was synthesized on solid-phase by adaptation of an Fmoc technique using a series of PyBOP/HOBt mediated coupling reactions. Sequence selectivity and complex stabilities were characterized by spectrofluorometric titrations and thermal melting studies. The polyamide of the conjugate was observed to bind in a hairpin motif forming 1:1 conjugate:dsDNA complexes. The conjugate is able to recognize nine contiguous A/T bps, discriminating from the sequences containing fewer than nine contiguous A/T bps.     

Studies on the retention of plasmid DNA and Escherichia coli nucleic acids by hydrophobic interaction chromatography

This work presents studies on the interactions of supercoiled plasmid DNA and Escherichia coli genomic DNA (gDNA) and RNA, with an hydrophobic interaction chromatography (HIC) gel, obtained by derivatisation of Sepharose CL-6B with 1,4-butanediol diglycidyl ether. Nucleic acids purified from E. coli were injected separately in the above HIC column and eluted with 1.5 M (NH₄)₂SO₄ in the buffer. The column was able to separate single-stranded from double-stranded nucleic acids. RNA and denatured gDNA were retarded in a different way due to the interactions of the exposed hydrophobic bases with the ligands. Supercoiled plasmid DNA, on the contrary, eluted in the flowthrough. This revised version was published online in July 2006 with corrections to the Cover Date.     

ATP-dependent minor groove recognition of TA base pairs is required for template melting by the E1 initiator protein

Template melting is an essential step in the initiation of DNA replication, but the mechanism of template melting is unknown for any replicon. Here we demonstrate that melting of the bovine papillomavirus type 1 ori is a sequence-dependent process which relies on specific recognition of TA base pairs in the minor groove by the E1 initiator. We show that correct template melting is a prerequisite for the formation of a stable double hexamer with helicase activity and that ori mutants that fail to melt correctly are defective for ori unwinding and DNA replication in vivo. Our results also indicate that melting of the DNA is achieved by destabilization of the double helix along its length through multiple interactions with E1, each of which is responsible for melting of a few base pairs, resulting in the extensive melting that is required for initiation of DNA replication.     

Expanding the recognition of the minor groove of DNA by incorporation of beta-alanine in hairpin polyamides

In order to expand the recognition code by hairpin polyamides to include DNA sequences of the type 5'-CWWC-3' two polyamides, PyPyPyPy-(R)(H2N)gamma-ImPyPyIm-beta-Dp (1) and PyPyPyPy-(R)(H2N)gamma-ImPy-beta-Im-beta-Dp (2) were synthesized which have in common an Py/Im pair in the terminal position for targeting C x G but differ with respect to internal placement of a beta-alanine residue. The equilibrium association constants (Ka) were determined at four DNA sites which differ at a single common position, 5'-TNTACA-3' (N = T, A, G, C). Quantitative DNase I footprint titration experiments reveal that the eight-ring hairpin PyPyPyPy-(R)(H2N)gamma-ImPyPyIm-beta-Dp (1) binds the four binding sites with similar affinities, Ka = 1.3-1.9 x 10(10) M(-1) indicating that there is no preference for the position N. In contrast, a redesigned polyamide PyPyPyPy-(R)(H2N)gamma-ImPy-beta-Im-beta-Dp (2) that places an internal flexible aliphatic beta-alanine to the 5'-side of a key imidazole group bound the match site 5'-TCTACA-3' with high affinity and good sequence discrimination (Ka(match) = 4.9 x 10(10) M(-1) and the single base pair mismatch sites with 5- to 25-fold lower affinity). These results expand the repertoire of sequences targetable by hairpins and emphasize the importance of beta-alanine as a key element for minor groove recognition.     

A chromatographic approach to the determination of relative free energies of interaction between hydrophobic and amphiphilic amino acid side chains.

A liquid chromatographic stationary phase was prepared by covalently binding to the surface of microparticulate silica gel functionality (benzylsilane), which mimics the side chain of the amino acid phenylalanine. The chromatographic retentions of the N-acetyl C-(N'-methyl) amides of various hydrophobic and amphiphilic amino acids on this stationary phase were measured using an aqueous mobile phase. A retention order of Gly < Ala < Cys < Val < Met < Pro < Ile < Leu < Tyr < Phe < Trp is seen at room temperature. Chromatographic retentions were used to derive free energies of adsorption of the amino acid derivatives on the chromatographic support relative to that of the glycine derivative. The temperature dependencies of the retention of aromatic and aliphatic amino acid derivatives differ in curvature, indicating a qualitative difference in the absorption mechanism. An adsorption model for retention is proposed, and arguments are made as to the suitability of an adsorption model for describing the contacts between amino acid side chains during the initial steps of protein folding.     

A novel method reveals that solvent water favors polyproline II over beta-strand conformation in peptides and unfolded proteins: conditional hydrophobic accessible surface area (CHASA)

In aqueous solution, the ensemble of conformations sampled by peptides and unfolded proteins is largely determined by their interaction with water. It has been a long-standing goal to capture these solute-water energetics accurately and efficiently in calculations. Historically, accessible surface area (ASA) has been used to estimate these energies, but this method breaks down when applied to amphipathic peptides and proteins. Here we introduce a novel method in which hydrophobic ASA is determined after first positioning water oxygens in hydrogen-bonded orientations proximate to all accessible peptide/protein backbone N and O atoms. This conditional hydrophobic accessible surface area is termed CHASA. The CHASA method was validated by predicting the polyproline-II (P(II)) and beta-strand conformational preferences of non-proline residues in the coil library (i.e., non-alpha-helix, non-beta-strand, non-beta-turn library derived from X-ray elucidated structures). Further, the method successfully rationalizes the previously unexplained solvation energies in polyalanyl peptides and compares favorably with published experimentally determined P(II) residue propensities. We dedicate this paper to Frederic M. Richards.