Conformational properties of B-Z junctions in DNA.

The structural consequences of specific base sequences in DNA can exert a strong influence on the function of DNA. It has previously been reported that the presence of multiple B-Z conformational junctions in constructed DNA oligomers results in unusually enhanced electrophoretic gel mobilities of these oligomers [Winkle, S. A., & Sheardy, R. D. (1990) Biochemistry 29, 6514-6521]. In order to investigate this phenomenon further, we designed and synthesized several DNA oligomers capable of pure Z or B-Z junction formation for polyacrylamide gel electrophoresis studies. The results indicate that both pure Z-DNA and polymorphic B-Z-DNA oligomers exhibit unusual gel migratory properties. The results of gel mobility studies in the absence and presence of cobalt hexamine indicate that a B-Z junction corresponds to a stiff bend of the helix axis, with two or more conformers accessible at the junction site. This is a different bend and mechanism than that in oligo(A) tracts.     

The role of material properties for the mechanical adaptation at branch junctions

Branch junctions are mechanically particularly interesting areas of trees, because they have to withstand a combination of static and dynamic loads, from the stem as well as from the branch. In the present work, the local adaptation of material properties at branch junctions was assessed by mapping microfibril angle and tissue density. Images of the density distribution were obtained by computer tomography (CT). Wide angle X-ray scattering (WAXS) was used to determine the microfibril angle distribution with high-resolution around the junctions. The stem tissue around the junctions showed increased density and microfibril angle, which points towards an optimisation for fracture toughness. The tissue at the branch bases showed low density combined with high MFA, which provides deformability and flexibility and might act as protection of the stem against load transmission from the branch.     

Conformational changes in phosphoglucose isomerase induced by ligand binding

Phosphoglucose isomerase (PGI; EC 5.3.1.9) is the second enzyme in glycolysis, where it catalyzes the isomerization of D-glucose-6-phosphate to D-fructose-6-phosphate. It is the same protein as autocrine motility factor, differentiation and maturation mediator, and neuroleukin. Here, we report a new X-ray crystal structure of rabbit PGI (rPGI) without ligands bound in its active site. The structure was solved at 1.8A resolution by isomorphous phasing with a previously solved X-ray crystal structure of the rPGI dimer containing 6-phosphogluconate in its active site. Comparison of the new structure to previously reported structures enables identification of conformational changes that occur during binding of substrate or inhibitor molecules. Ligand binding causes an induced fit of regions containing amino acid residues 209-215, 245-259 and 385-389. This conformational change differs from the change previously reported to occur between the ring-opening and isomerization steps, in which the helix containing residues 513-521 moves toward the bound substrate. Differences between the liganded and unliganded structures are limited to the region within and close to the active-site pocket.     

The kinetic properties of cruciform extrusion are determined by DNA base-sequence.

The extrusion kinetics of two cruciforms derived from unrelated DNA sequences differ markedly. Kinetic barriers exist for both reactions, necessitating elevated temperatures before extrusion proceeds at measureable speeds, but the dependence upon temperature and ionic strength is quite different for the two sequences. One, the ColE1 inverted repeat, exhibits a remarkably great temperature dependence of reaction rate and is suppressed by moderate amounts of NaCl or MgCl2. In contrast, the other, a synthetic inverted repeat present in pIRbke8, shows more modest temperature dependence and has a requirement for the presence of salt, with optimal concentrations being 50 mM NaCl or 100 microM MgCl2. Under optimal conditions, cruciform extrusion rates are fast (t1/2 less than 60m) at 37 degrees C for both sequences at native superhelix densities. In 50 mM NaCl the pIRbke8 inverted repeat is characterised by an Arrhenius activation energy of 42.4 +/- 3.2 kcal mole -1. The differences in kinetic properties between the two sequences indicate that DNA base sequence is itself an important factor in determining cruciform kinetics, and possibly even in the selection of the mechanistic pathway.     

DNA binding affinity and sequence permutation preference of the telomere protein from Euplotes crassus

Telomere end binding proteins from diverse organisms use various forms of an ancient protein structure to recognize and bind with single-strand DNA found at the ends of telomeres. To further understand the biochemistry and evolution of these proteins, we have characterized the DNA binding properties of the telomere end binding protein from Euplotes crassus (EcTEBP). EcTEBP and its predicted amino-terminal DNA-binding domain, EcTEBP-N, were expressed in Escherichia coli and purified. Each protein formed stoichiometric (1:1) complexes with single-strand DNA oligos derived from the precisely defined d(TTTTGGGGTTTTGG) sequence found at DNA termini in Euplotes. Dissociation constants for DNA x EcTEBP and DNA x EcTEBP-N complexes were comparable: K(D-DNA) = 38 +/- 2 nM for the full-length protein and K(D-DNA) = 60 +/- 4 nM for the N-terminal domain, indicating that the N-terminal domain retains a high affinity for DNA even in the absence of potentially stabilizing moieties located in the C-terminal domain. Rate constants for DNA association and DNA dissociation corroborated a slightly improved DNA binding performance for the full-length protein (ka = 45 +/- 4 microM(-1) s(-1), kd = 0.10 +/- 0.02 s(-1)) relative to that of the N-terminal domain (ka = 18 +/- 1 microM(-1) s(-1), kd = 0.15 +/- 0.01 s(-1)). Equilibrium dissociation constants measured for sequence permutations of the telomere repeat spanned the range of 55-1400 nM, with EcTEBP and EcTEBP-N binding most tightly to d(TTGGGGTTTTGG), the sequence corresponding to that of mature DNA termini. Additionally, competition experiments showed that EcTEBP recognizes and binds the telomere-derived 14-nucleotide DNA in preference to shorter 5'-truncation variants. Compared with the results for multisubunit complexes assembled with telomere single-strand DNA from Oxytricha nova, our results highlight the relative simplicity of the E. crassus system where a telomere end binding protein has biochemical properties indicating one protein subunit caps the single-strand DNA.     

Single molecule fluorescence analysis of branch migration of holliday junctions: effect of DNA sequence

The Holliday junction is a central intermediate in various genetic processes including homologous, site-specific recombination and DNA replication. Recent single molecule FRET experiments led to the model for branch migration as a stepwise stochastic process in which the branch migration hop is terminated by the folding of the junction. In this article, we studied the effect of the sequence on Holliday junction dynamics and branch migration process. We show that a GC pair placed at the border of the homologous region almost prevents the migration into this position. At the same time, insertion of a GC pair into the middle of the AT tract does not show this effect, however when the junction folds at this position, it resides at this position much longer time in comparison to the folding at AT pairs. Two contiguous GC pairs do not block migration as well and generally manifest the same effect as one GC pair—the junction when it folds resides at these positions for a relatively long time. The same elevated residence time was obtained for the design with the homology region that consists of only GC pairs. These data suggest a model for branch migration in which the sequence modulates the overall stochastic process of the junction dynamics and branch migration by the variability of the time that the junction dwells before making a migration hop.     

RuvAB-directed branch migration of individual Holliday junctions is impeded by sequence heterology

The Holliday junction, the key intermediate of recombination, is generated by strand exchange resulting in a covalent connection between two recombining DNA molecules. Translocation of a Holliday junction along DNA, or branch migration, progressively exchanges one DNA strand for another and determines the amount of information that is transferred between two recombining partners. In Escherichia coli, the RuvAB protein complex promotes rapid and unidirectional branch migration of Holliday junctions. We have studied translocation of Holliday junctions using a quantitative biochemical system together with a 'single-molecule' branch migration assay. We demonstrate that RuvAB translocates the junctions through identical DNA sequences in a processive manner with a broad distribution of individual branch migration rates. However, when the complex encounters short heterologous sequences, translocation of the Holliday junctions is impeded. We conclude that translocation of the junctions through a sequence heterology occurs with a probability of bypass being determined both by the length of the heterologous region and the lifetime of the stalled RuvAB complex.     

Conformational properties of topoisomerase II inhibitors and sequence specificity of DNA cleavage

The sequence specificity of topoisomerase-II-mediated DNA cleavage, stimulated by 2-methyl-9-hydroxy ellipticinium and 4′, 5′,7-trihydroyflavone (genistein) was investigated by sequencing analysis of DNA cleavage sites and molecular modeling techniques. The former drug exhibits a marked preference for a T base at the position immediately preceding the cleavage site (?1). The latter shares the preference for the same base, with an additional preference for a thymine at position +1. The cleavage intensity patterns for the two drugs differ considerably. From a conformational point of view, ellipticinium and genistein exhibit similar overall shape and dimensions. However, the fused ring system in the former generates a planar structure whereas the single bond, connecting the two aromatic portions in the latter, allows internal rotation. The most stable conformation of genistein corresponds to a deviation of about 40° from planarity. A computer-assisted analysis was carried out to compare the steric and electrostatic properties of the two compounds. Two types of preferred (energetically almost degenerate) alignment for the two molecules were found. One corresponds to overlapping of the 9-hydroxyl containing ring of ellipticinium with the 4′-hydroxyphenyl moiety of genistein, the other envisages the same moiety of ellipticine superimposed to the hydroxyl-benzopyrone portion of genistein. The structural similarities of the test drugs might account for the common preference for stimulation of DNA cleavage at position +1, whereas the different possible arrangements of genistein in the cleavable complex could explain both the additional +1 specificity exhibited by this compound and the differences in cleavage intensity patterns observed in comparison to ellipticinium.     

DNA binding of dinuclear iron(II) metallosupramolecular cylinders. DNA unwinding and sequence preference

[Fe₂L₃]⁴⁺ (L = C₂₅H₂₀N₄) is a synthetic tetracationic supramolecular cylinder (with a triple helical architecture) that targets the major groove of DNA and can bind to DNA Y-shaped junctions. To explore the DNA-binding mode of [Fe₂L₃]⁴⁺, we examine herein the interactions of pure enantiomers of this cylinder with DNA by biochemical and molecular biology methods. The results have revealed that, in addition to the previously reported bending of DNA, the enantiomers extensively unwind DNA, with the M enantiomer being the more efficient at unwinding, and exhibit preferential binding to regular alternating purine–pyrimidine sequences, with the M enantiomer showing a greater preference. Also, interestingly, the DNA binding of bulky cylinders [Fe₂(L-CF₃)₃]⁴⁺ and [Fe₂(L-Ph)₃]⁴⁺ results in no DNA unwinding and also no sequence preference of their DNA binding was observed. The observation of sequence-preference in the binding of these supramolecular cylinders suggests that a concept based on the use of metallosupramolecular cylinders might result in molecular designs that recognize the genetic code in a sequence-dependent manner with a potential ability to affect the processing of the genetic code.     

DNA recombination: holliday junctions dynamics and branch migration

Holliday junctions are critical intermediates for homologous, site-specific recombination, DNA repair, and replication. A wealth of structural information is available for immobile four-way junctions, but the controversy on the mechanism of branch migration of Holliday junctions remains unsolved. Two models for the mechanism of branch migration were suggested. According to the early model of Alberts-Meselson-Sigal (Sigal, N., and Alberts, B. (1972) J. Mol. Biol. 71, 789-793 and Meselson, M. (1972) J. Mol. Biol. 71, 795-798), exchanging DNA strands around the junction remain parallel during branch migration. Kinetic studies of branch migration (Panyutin, I. G., and Hsieh, P. (1994) Proc. Natl. Acad. Sci. U. S. A. 91, 2021-2025) suggest an alternative model in which the junction adopts an extended conformation. We tested these models using a Holliday junction undergoing branch migration and time-lapse atomic force microscopy, an imaging technique capable of imaging DNA dynamics. The single molecule atomic force microscopy experiments performed in the presence and in the absence of divalent cations show that mobile Holliday junctions adopt an unfolded conformation during branch migration that is retained despite a broad range of motion in the arms of the junction. This conformation of the junction remains unchanged until strand separation. The data obtained support the model for branch migration having the extended conformation of the Holliday junction.     

The relative ligand binding preference of the murine ileal lipid binding protein

The ileal lipid binding protein (ILBP), a member of the intracellular lipid binding protein family, is a 14-kDa protein that has bile and fatty acids as possible physiological ligands. The ligand binding specificity of this protein is not well characterized. Therefore, we studied the lipid binding activity of purified recombinant murine ILBP (mILBP) in vitro. These studies demonstrated by direct analysis the interaction of mILBP with naturally occurring bile and fatty acids. The rank order of binding preference for fatty acids, or unconjugated and conjugated bile acids, was assessed. Among fatty acids, mILBP preferred species that had longer chain length and increased saturation, similar to other members of the intracellular lipid binding protein family. Among the bile acids, mILBP showed the greatest preference for conjugated species that contained a doubly hydroxylated steroid moiety. The results demonstrate that mILBP exhibits a preference for certain species of bile and fatty acids.     

G-quadruplex conformation and dynamics are determined by loop length and sequence

The quadruplex forming G-rich sequences are unevenly distributed throughout the human genome. Their enrichment in oncogenic promoters and telomeres has generated interest in targeting G-quadruplex (GQ) for an anticancer therapy. Here, we present a quantitative analysis on the conformations and dynamics of GQ forming sequences measured by single molecule fluorescence. Additionally, we relate these properties to GQ targeting ligands and G4 resolvase 1 (G4R1) protein binding. Our result shows that both the loop (non-G components) length and sequence contribute to the conformation of the GQ. Real time single molecule traces reveal that the folding dynamics also depend on the loop composition. We demonstrate that GQ-stabilizing small molecules, N-methyl mesoporphyrin IX (NMM), its analog, NMP and the G4R1 protein bind selectively to the parallel GQ conformation. Our findings point to the complexity of GQ folding governed by the loop length and sequence and how the GQ conformation determines the small molecule and protein binding propensity.     

Molecular recognition between oligopeptides and nucleic acids. DNA sequence specificity and binding properties of an acridine-linked netropsin hybrid ligand

The binding to DNA of a mixed function ligand (NETGA) is described, in which a potential intercalating group, an acridine moiety, is incorporated at the carboxyl terminus of the minor groove binding oligopeptide netropsin skeleton. Scatchard analysis of absorption data provided evidence of two modes of binding to DNA with K1 = 9.1 x 10(5) M-1 at low r values (0.003-0.1), and a binding site size n = 10, indicative of binding of both moeities. At high binding ratios (greater than 0.1), K2 = 0.9 x 10(5) M-1 and n = 5 corresponding to external binding. Complementary strand MPE footprinting on a pBR322 restriction fragment showed NETGA binds to 5'-AAAT like netropsin. It causes enhanced cleavage by MPE, particularly at G-C rich sequences and remote from the preferred binding sites. Viscometry measurements provided evidence for biphasic modes of the two binding portions of NETGA. Fluorescence polarization and linear dichroism measurements were in accord with distinct modes of interaction of the acridine (intercalation) and oligopeptide (minor groove binding) portions of NETGA. LD measurements on NETGA indicate that the oligopeptide moiety (netropsin-like) has an orientation typical of minor groove binders, whereas the degree of intercalation of the acridine group is decreased by association of the oligopeptide moiety.     

Conformational constraint in protein ligand design and the inconsistency of binding entropy

It is an accepted practice in ligand design to introduce conformational constraint with the expectation of improving affinity, justified by the theoretical possibility that an unfavorable change in binding entropy will be reduced. This rationale of minimizing the entropic penalty through imposing structural constraints upon a ligand, however, has been voiced more often than verified. Here we examine three modified cyclic peptides, along with multiple versions of their linear control analogs, and determine their thermodynamic parameters when binding the same host, the third PDZ domain (PDZ3) of the mammalian postsynaptic density-95 (PSD-95) protein. To begin a two-stage investigation, the initial evaluation involved solution binding studies with isothermal titration calorimetry (ITC), which provided the changes in Gibbs free energy (DeltaG), enthalpy (DeltaH), and entropy (TDeltaS) upon formation of the protein-ligand complex. In the second stage, a selected macrocycle along with two matched linear controls were subjected to more rigorous analysis by ITC, which included (1) change in heat of buffer ionization (DeltaH(ion)) titrations, to examine the role of proton transfer events; (2) change in heat capacity (DeltaC(p)) determinations, to indirectly probe the nature of the binding surface; and (3) osmotic stress experiments, to evaluate desolvation effects and quantitate water release. Together, these demonstrate that the entropic relationship between a macrocyclic ligand and a linear counterpart can be a complex one that is difficult to rationalize. Further, the addition of constraint can, counterintuitively, lead to a less favorable change in binding entropy. This underscores the need to use matched linear control ligands to assure that comparisons are made in a meaningful manner.     

Modeling of DNA condensation and decondensation caused by ligand binding

A theoretical method for computer modeling of DNA condensation caused by ligand binding is developed. In the method, starting (s) and condensed (c) states are characterized by different free energies for ligand free DNA (F(s) and F(c) respectively), ligand binding constants (K(s) and K(c)) and stoichiometry dependent parameters (c(sm) and c(cm) - maximum relative concentration of bound ligands (per base pair) for starting and condensed state respectively). The method allows computation of the dependence of the degree of condensation (the fraction of condensed DNA molecules) on ligand concentration. Calculations demonstrate that condensation transition occurs under an increase in ligand concentration if F(s) < F(c) (i.e. S(sc) = exp [- (F(c) - F(s)) / (RT)], the equilibrium constant of the s-c transition, is low (S(sc) < 1)) and K(s) < K(c). It was also found that condensation is followed by decondensation at high ligand concentration if the condensed DNA state provides the number of sites for ligand binding less than the starting state (c(sm) > c(cm)). A similar condensation-decondensation effect was found in recent experimental studies. We propose its simple explanation.