New concepts in protein-DNA recognition: sequence-directed DNA bending and flexibility

DNA structural features of eukaryotic TATA-containing and TATA-less promoters

Eukaryotic genes can be broadly classified as TATA‐containing and TATA‐less based on the presence of TATA box in their promoters. Experiments on both classes of genes have revealed a disparity in the regulation of gene expression and cellular functions between the two classes. In this study, we report characteristic differences in promoter sequences and associated structural properties of the two categories of genes in six different eukaryotes. We have analyzed three structural features, DNA duplex stability, bendability, and curvature along with the distribution of A‐tracts, G‐quadruplex motifs, and CpG islands. The structural feature analyses reveal that while the two classes of gene promoters are distinctly different from each other, the properties are also distinguishable across the six organisms.

A theoretical investigation of DNA dynamics and desolvation kinetics for zinc finger proteinZif268

Background

Transcription factors, regulating the expression inventory of a cell, interact with its respective DNA subjugated by a specific recognition pattern, which if well exploited may ensure targeted genome engineering. The mostly widely studied transcription factors are zinc finger proteins that bind to its target DNA via direct and indirect recognition levels at the interaction interface. Exploiting the binding specificity and affinity of the interaction between the zinc fingers and the respective DNA can help in generating engineered zinc fingers for therapeutic applications. Experimental evidences lucidly substantiate the effect of indirect interaction like DNA deformation and desolvation kinetics, in empowering ZFPs to accomplish partial sequence specificity functioning around structural properties of DNA. Exploring the structure-function relationships of the existing zinc finger-DNA complexes at the indirect recognition level can aid in predicting the probable zinc fingers that could bind to any target DNA. Deformation energy, which defines the energy required to bend DNA from its native shape to its shape when bound to the ZFP, is an effect of indirect recognition mechanism. Water is treated as a co-reactant for unfurling the affinity studies in ZFP-DNA binding equilibria that takes into account the unavoidable change in hydration that occurs when these two solvated surfaces come into contact.

Results

Aspects like desolvation and DNA deformation have been theoretically investigated based on simulations and free energy perturbation data revealing a consensus in correlating affinity and specificity as well as stability for ZFP-DNA interactions. Greater loss of water at the interaction interface of the DNA calls for binding with higher affinity, eventually distorting the DNA to a greater extent accounted by the change in major groove width and DNA tilt, stretch and rise.

Conclusion

Most prediction algorithms for ZFPs do not account for water loss at the interface. The above findings may significantly affect these algorithms. Further the sequence dependent deformation in the DNA upon complexation with our prototype as well as preference of bases at the 2^nd and 3^rd position of the repeating triplet provide an absolutely new insight about the indirect interactions undergoing a change that have not been probed yet.

Structural waters in the minor and major grooves of DNA--a major factor governing structural adjustments of the A-T mini-helix

The role of microhydration in structural adjustments of the AT-tract in B-DNA was studied at the B97-D/def2-SV(P) level. The (dA:dT)5 complexes with 10 water molecules in minor and 15 water molecules in major grooves were studied. The obtained network of hydrogen bonds revealed the dependence between the groove width and the types of water patterns. In the minor groove, the following patterns were observed: interstrand one-water bridges similar to that of the Dickerson "water spine" and interstrand two-water bridges. The network of structural waters in the major groove is more diverse than that in the minor groove, which agrees with crystallographic data. As the major groove is wider, it is enriched by water molecules forming two- and three-water bridges. Results suggest the nucleobase-water interactions in both grooves prevent AT-tract twisting and its "collapse" along the minor groove. Whereby, a helix structure with narrow minor and wide major grooves is formed. The structural waters affect the polynucleotide conformation so that it becomes similar to poly(dA)·poly(dT) in fibers and acquires features of the A-tracts in DNA in solution. We suggest that formation of specific water patterns in both grooves is the factor responsible for stabilization of A-tracts with a narrowed minor groove, leading in turn to their strong intrinsic bending in DNA.

Computational approaches for predicting the binding sites and understanding the recognition mechanism of protein-DNA complexes

Protein-DNA recognition plays an important role in the regulation of gene expression. Understanding the influence of specific residues for protein-DNA interactions and the recognition mechanism of protein-DNA complexes is a challenging task in molecular and computational biology. Several computational approaches have been put forward to tackle these problems from different perspectives: (i) development of databases for the interactions between protein and DNA and binding specificity of protein-DNA complexes, (ii) structural analysis of protein-DNA complexes, (iii) discriminating DNA-binding proteins from amino acid sequence, (iv) prediction of DNA-binding sites and protein-DNA binding specificity using sequence and/or structural information, and (v) understanding the recognition mechanism of protein-DNA complexes. In this review, we focus on all these issues and extensively discuss the advancements on the development of comprehensive bioinformatics databases for protein-DNA interactions, efficient tools for identifying the binding sites, and plausible mechanisms for understanding the recognition of protein-DNA complexes. Further, the available online resources for understanding protein-DNA interactions are collectively listed, which will serve as ready-to-use information for the research community.

Scaling behavior of nucleotide cluster in DNA sequences

In this paper we study the scaling behavior of nucleotide cluster in 11 chromosomes of Encephalitozoon cuniculi Genome. The statistical distribution of nucleotide clusters for 11 chromosomes is characterized by the scaling behavior of P(S) proportional, variant e(-alphaS), where S represents nucleotide cluster size. The cluster-size distribution P(S(1)+S(2)) with the total size of sequential C-G cluster and A-T cluster S(1)+S(2) were also studied. P(S(1)+S(2)) follows exponential decay. There does not exist the case of large C-G cluster following large A-T cluster or large A-T cluster following large C-G cluster. We also discuss the relatively random walk length function L(n) and the local compositional complexity of nucleotide sequences based on a new model. These investigations may provide some insight into nucleotide cluster of DNA sequence.

DNA sequence-specific recognition by a transcriptional regulator requires indirect readout of A-tracts

The bacteriophage Ø29 transcriptional regulator p4 binds to promoters of different intrinsic activities. The p4–DNA complex contains two identical protomers that make similar interactions with the target sequence 5′-AACTTTTT-15 bp-AAAATGTT-3′. To define how the various elements in the target sequence contribute to p4's affinity, we studied p4 binding to a series of mutated binding sites. The binding specificity depends critically on base pairs of the target sequence through both direct as well as indirect readout. There is only one specific contact between a base and an amino acid residue; other contacts take place with the phosphate backbone. Alteration of direct amino acid–base contacts, or mutation of non-contacted A·T base pairs at A-tracts abolished binding. We generated three 5 ns molecular dynamics (MD) simulations to investigate the basis for the p4–DNA complex specificity. Recognition is controlled by the protein and depends on DNA dynamic properties. MD results on protein–DNA contacts and the divergence of p4 affinity to modified binding sites reveal an inherent asymmetry, which is required for p4-specific binding and may be crucial for transcription regulation.

DNA deformation energy as an indirect recognition mechanism in protein-DNA interactions

Proteins that bind to specific locations in genomic DNA control many basic cellular functions. Proteins detect their binding sites using both direct and indirect recognition mechanisms. Deformation energy, which models the energy required to bend DNA from its native shape to its shape when bound to a protein, has been shown to be an indirect recognition mechanism for one particular protein, Integration Host Factor (IHF). This work extends the analysis of deformation to two other DNA-binding proteins, CRP and SRF, and two endonucleases, I-CreI and I-PpoI. Known binding sites for all five proteins showed statistically significant differences in mean deformation energy as compared to random sequences. Binding sites for the three DNA-binding proteins and one of the endonucleases had mean deformation energies lower than random sequences. Binding sites for I-PpoI had mean deformation energy higher than random sequences. Classifiers that were trained using the deformation energy at each base pair step showed good cross-validated accuracy when classifying unseen sequences as binders or nonbinders. These results support DNA deformation energy as an indirect recognition mechanism across a wider range of DNA-binding proteins. Deformation energy may also have a predictive capacity for the underlying catalytic mechanism of DNA-binding enzymes.

Structure-function analysis of SWI2/SNF2 enzymes

Biochemical and structural progress over the last years has revealed that SWI2/SNF2 family chromatin remodeling or DNA repair enzymes are molecular motors that transport duplex DNA along a helicase-like domain using ATP-hydrolysis. The screw motion of DNA along the active site probably generates the force to disrupt chromatin or other protein:DNA complexes. In this chapter, we describe biochemical and structural approaches to study the molecular mechanism of SWI2/SNF2 enzymes. In particular, we describe assays to monitor DNA dependent ATPase activity, translocation on duplex DNA, and DNA distortion activity. We also describe recent progress in the crystallization and structure determination of SWI2/SNF2 enzymes in complex with duplex DNA.

Effect of DNA structural flexibility on promoter strength--molecular dynamics studies of E. coli promoter sequences

To study possible correlations between promoter activity and the structural flexibility of the DNA helix, we have carried out unrestrained molecular dynamics simulations of the -10 consensus region sequence and five variants forming the -10 region of various Escherichia coli promoter sequences. Analyses of the trajectories obtained from the simulations show that the consensus sequence has a pattern of two structurally flexible nucleotide steps sandwiched between two stiff steps. In the other sequences, this pattern varies in consonance with the change in the sequence. The variations in the patterns show correlation with the promoter strength.

Role of inter and intramolecular interactions in protein-DNA recognition

Protein-DNA recognition plays an essential role in the regulation of gene expression. Regulatory proteins are known to recognize specific DNA sequences directly through atomic contacts between protein and DNA, and/or indirectly through the conformational properties of the DNA. In this work, we have analyzed the specificity of intermolecular interactions by statistical analysis of base-amino acid interactions within protein-DNA complexes as well as the computer simulations of base-amino acid interactions. The specificity of the intramolecular interactions was studied by statistical analysis of the sequence-dependent DNA conformational parameters and the elastic properties of DNA. Systematic comparison of these specificities in a large number of protein-DNA complexes revealed that both intermolecular and intramolecular interactions contribute to the specificity of protein-DNA recognition, and their relative contributions vary depending upon the protein-DNA complex. We demonstrated that combination of the intermolecular and intramolecular energies leads to enhanced specificity and the combined energy could explain experimental data on binding affinity changes caused by base mutations. These results provided new insight into the relationship between specificity and structure in the process of protein-DNA recognition, which would lead to prediction of specific protein-DNA binding sites.

Free energy contributions to direct readout of a DNA sequence

The energetic contributions of individual DNA-contacting side chains to specific DNA recognition in the human papillomavirus 16 E2C-DNA complex is small (less than 1.0 kcal mol(-1)), independent of the physical and chemical nature of the interaction, and is strictly additive. The sum of the individual contributions differs 1.0 kcal mol(-1) from the binding energy of the wild-type protein. This difference corresponds to the contribution from the deformability of the DNA, known as "indirect readout." Thus, we can dissect the energetic contribution to DNA binding into 90% direct and 10% indirect readout components. The lack of high energy interactions indicates the absence of "hot spots," such as those found in protein-protein interfaces. These results are compatible with a highly dynamic and "wet" protein-DNA interface, yet highly specific and tight, where individual interactions are constantly being formed and broken.

Comparative binding of p53 to its promoter and DNA recognition elements

Tumor suppressor p53 is a transcription factor that transactivates a wide range of genes, including those in DNA repair, cell cycle arrest, apoptosis and its own degradation. To estimate the role of selectivity in binding to its promoters, we measured the binding affinities of a tetrameric p53 construct (p53CT) in vitro with 20 of its recognition elements from a variety of representative genes. The binding of full length p53 to four representative sequences exactly paralleled the affinities to p53CT. The binding of p53 to different recognition elements was co-operative and the affinities varied by up to 50-fold. p53 bound with high affinity to the recognition elements of all the genes involved in cell cycle arrest and some of the genes in apoptosis. All of the lower affinity-binding sites were in genes involved in apoptosis. Our quantitative-binding data were in agreement with published cell-based assays. The regulation of p53 activity is in part determined through the specificity of its DNA-binding interactions.

Investigating local conformations of double-stranded DNA by low-energy circular dichroism of pyrrolo-cytosine

Local base stacking and conformational mobility play a major role in the structure and function of nucleic acids. We have recently shown that the low-energy CD spectrum of 2-aminopurine (2-AP), i.e., the CD spectral region above 300 nm, can be used to monitor conformational changes in polynucleotides at or near mono- and dinucleotide 2-AP residues that replace adenine residues in DNA and RNA. Here, we extend this technique to pyrrolo-cytosine (PC), a fluorescent analogue of cytosine. The low-energy CD spectrum of a PC dinucleotide in dsDNA exhibits an exciton couplet with two bands of opposite sign centered at 350 nm. This signal is characteristic of base stacking between adjacent PC residues in a right helical conformation. Isolated PC nucleotide residues inserted into polynucleotide chains also display chirality that reflects the asymmetric environment of their sequence context. Thus, we show that the low-energy CD spectra of C(PC)A and A(PC)C sequences in dsDNA have opposite signs. It appears that the measurement of the low-energy CD spectra of PC residues will usefully complement 2-AP measurements by serving to characterize the local conformations and dynamics of nucleic acids near C residues and G-C base pairs.

Influence of DNA stiffness in protein–DNA recognition

Protein-DNA recognition plays an essential role in the regulation of gene expression. The protein-DNA binding specificity is based on direct atomic contacts between protein and DNA and/or the conformational properties of DNA. In this work, we have analyzed the influence of DNA stiffness (E) to the specificity of protein-DNA complexes. The average DNA stiffness parameters for several protein-DNA complexes have been computed using the structure based sequence dependent stiffness scale. The relationship between DNA stiffness and experimental protein-DNA binding specificity has been brought out. We have investigated the importance of DNA stiffness with the aid of experimental free energy changes (DeltaDeltaG) due to binding in several protein-DNA complexes, such as, ETS proteins, 434, lambda, Mnt and trp repressors, 434 cro protein, EcoRV endonuclease V and zinc fingers. We found a correlation in the range 0.65-0.97 between DeltaDeltaG and E in these examples. Further, we have qualitatively analyzed the effect of mutations in the target sequence of lambda repressor and we observed that the DNA stiffness could correctly identify 70% of the correct bases among the considered nine positions.

New elements of the termination of transcription in prokaryotes

We report here that phased runs of adenines and thymines are very frequent in the neighborhood of 3' of the coding regions of Escherichia coli and Bacillus subtilis. These findings suggest that the DNA curvature could affect transcription termination either directly, through contacts with RNA polymerase, or indirectly, via contacts with some regulatory proteins.

The role of water in protein-DNA recognition

Is it by design or by default that water molecules are observed at the interfaces of some protein-DNA complexes? Both experimental and theoretical studies on the thermodynamics of protein-DNA binding overwhelmingly support the extended hydrophobic view that water release from interfaces favors binding. Structural and energy analyses indicate that the waters that remain at the interfaces of protein-DNA complexes ensure liquid-state packing densities, screen the electrostatic repulsions between like charges (which seems to be by design), and in a few cases act as linkers between complementary charges on the biomolecules (which may well be by default). This review presents a survey of the current literature on water in protein-DNA complexes and a critique of various interpretations of the data in the context of the role of water in protein-DNA binding and principles of protein-DNA recognition in general.

Binding of phage Phi29 architectural protein p6 to the viral genome: evidence for topological restriction of the phage linear DNA

Bacillus subtilis phage Phi29 protein p6 is required for DNA replication and promotes the switch from early to late transcription. In vivo it binds all along the viral linear DNA, which suggests a global role as an architectural protein; in contrast, binding to bacterial DNA is negligible. This specificity could be due to the p6 binding preference for less negatively supercoiled DNA, as is presumably the case with viral (with respect to bacterial) DNA. Here we demonstrate that p6 binding to Phi29 DNA is greatly increased when negative supercoiling is decreased by novobiocin; in addition, gyrase is required for DNA replication. This indicates that, although non-covalently closed, the viral genome is topologically constrained in vivo. We also show that the p6 binding to different Phi29 DNA regions is modulated by the structural properties of their nucleotide sequences. The higher affinity for DNA ends is possibly related to the presence of sequences in which their bendability properties favor the formation of the p6-DNA complex, whereas the lower affinity for the transcription control region is most probably due to the presence of a rigid intrinsic DNA curvature.

Intermolecular and Intramolecular Readout Mechanisms in Protein–DNA Recognition

Protein-DNA recognition plays an essential role in the regulation of gene expression. Regulatory proteins are known to recognize specific DNA sequences directly through atomic contacts (intermolecular readout) and/or indirectly through the conformational properties of the DNA (intramolecular readout). However, little is known about the respective contributions made by these so-called direct and indirect readout mechanisms. We addressed this question by making use of information extracted from a structural database containing many protein-DNA complexes. We quantified the specificity of intermolecular (direct) readout by statistical analysis of base-amino acid interactions within protein-DNA complexes. The specificity of the intramolecular (indirect) readout due to DNA was quantified by statistical analysis of the sequence-dependent DNA conformation. Systematic comparison of these specificities in a large number of protein-DNA complexes revealed that both intermolecular and intramolecular readouts contribute to the specificity of protein-DNA recognition, and that their relative contributions vary depending upon the protein-DNA complexes. We demonstrated that combination of the intermolecular and intramolecular energies derived from the statistical analyses lead to enhanced specificity, and that the combined energy could explain experimental data on binding affinity changes caused by base mutations. These results provided new insight into the relationship between specificity and structure in the process of protein-DNA recognition, which would lead to prediction of specific protein-DNA binding sites.

Sequence-dependent dynamics of duplex DNA: the applicability of a dinucleotide model

The short-time (submicrosecond) bending dynamics of duplex DNA were measured to determine the effect of sequence on dynamics. All measurements were obtained from a single site on duplex DNA, using a single, site-specific modified base containing a rigidly tethered, electron paramagnetic resonance active spin probe. The observed dynamics are interpreted in terms of single-step sequence-dependent bending force constants, determined from the mean squared amplitude of bending relative to the end-to-end vector using the modified weakly bending rod model. The bending dynamics at a single site are a function of the sequence of the nucleotides constituting the duplex DNA. We developed and examined several dinucleotide-based models for flexibility. The models indicate that the dominant feature of the dynamics is best explained in terms of purine- and pyrimidine-type steps, although distinction is made among all 10 unique steps: It was found that purine-purine steps (which are the same as pyrimidine-pyrimidine steps) were near average in flexibility, but the pyrimidine-purine steps (5' to 3') were nearly twice as flexible, whereas purine-pyrimidine steps were more than half as flexible as average DNA. Therefore, the range of stepwise flexibility is approximately fourfold and is characterized by both the type of base pair step (pyrimidine/purine combination) and the identity of the bases within the pair (G, A, T, or C). All of the four models considered here underscore the complexity of the dependence of dynamics on DNA sequence with certain sequences not satisfactorily explainable in terms of any dinucleotide model. These findings provide a quantitative basis for interpreting the dynamics and kinetics of DNA-sequence-dependent biological processes, including protein recognition and chromatin packaging.

The papillomavirus E2 proteins: structure, function, and biology

Nearly twenty years after the first high-resolution crystal structures of specific protein-DNA complexes were determined, the stereo-chemical basis for protein-DNA recognition remains an active area of investigation. One outstanding question is, how are proteins able to detect noncontacted sequences in their binding sites? The papillomavirus E2 proteins represent a particularly suitable group of proteins in which to examine the mechanisms of "indirect readout." Coordinated structural and thermodynamic studies of the E2-DNA interaction conducted over the past five years are summarized in this review. The data support a model in which the electrostatic properties of the individual E2 proteins correlate with their affinities for intrinsically flexible or rigidly prebent DNA targets.

Complex local dynamics in DNA on the picosecond and nanosecond time scales

Time-resolved Stokes shift measurements of the local structural relaxation of three DNA oligonucleotides are presented. Logarithmic relaxation is seen for over three decades in time (40 ps-40 ns), indicating a complex relaxation among a large number of conformational substates. The observed relaxation is the same in all the sequences. Sequence dependence of the localized dynamics of DNA does not appear within this time range. We infer that 30%-50% of the relaxation is faster than 40 ps, has a nonlogarithmic decay and has a sequence dependent amplitude.

Free-energy component analysis of 40 protein-DNA complexes: a consensus view on the thermodynamics of binding at the molecular level

Noncovalent association of proteins to specific target sites on DNA--a process central to gene expression and regulation--has thus far proven to be idiosyncratic and elusive to generalizations on the nature of the driving forces. The spate of structural information on protein--DNA complexes sets the stage for theoretical investigations on the molecular thermodynamics of binding aimed at identifying forces responsible for specific macromolecular recognition. Computation of absolute binding free energies for systems of this complexity transiting from structural information is a stupendous task. Adopting some recent progresses in treating atomic level interactions in proteins and nucleic acids including solvent and salt effects, we have put together an energy component methodology cast in a phenomenological mode and amenable to systematic improvements and developed a computational first atlas of the free energy contributors to binding in approximately 40 protein-DNA complexes representing a variety of structural motifs and functions. Illustrating vividly the compensatory nature of the free energy components contributing to the energetics of recognition for attaining optimal binding, our results highlight unambiguously the roles played by packing, electrostatics including hydrogen bonds, ion and water release (cavitation) in protein-DNA binding. Cavitation and van der Waals contributions without exception favor complexation. The electrostatics is marginally unfavorable in a consensus view. Basic residues on the protein contribute favorably to binding despite the desolvation expense. The electrostatics arising from the acidic and neutral residues proves unfavorable to binding. An enveloping mode of binding to short stretches of DNA makes for a strong unfavorable net electrostatics but a highly favorable van der Waals and cavitation contribution. Thus, noncovalent protein-DNA association is a system-specific fine balancing act of these diverse competing forces. With the advances in computational methods as applied to macromolecular recognition, the challenge now seems to be to correlate the differential (initial vs. final) energetics to substituent effects in drug design and to move from affinity to specificity.

Transcriptional activation by LR1 at the Emu enhancer and switch region sites

LR1 is a B cell-specific, sequence-specific duplex DNA binding activity which is induced in B cells carrying out class switch recombination. Here we identify several properties of LR1 which enable it to function in transcriptional regulation. We show that LR1 contributes to transcriptional activation by the Emu immunoglobulin heavy chain intron enhancer by binding to a site within the enhancer core. We further show that LR1 bends DNA upon binding. In addition, we show that LR1 is itself a bona fide transcriptional activator, as multimerized LR1 sites produce an element which can enhance transcription from a minimal promoter. In order for class switch recombination to occur, an activating signal must be transmitted via the Emu core, and both S regions targeted for recombination must be actively transcribed. The properties of LR1 that we have identified suggest distinct potential functions of LR1 duplex DNA binding activity in class switch recombination. First, LR1 may contribute to recombinational activation by the Emu core. Second, there are multiple potential LR1 duplex binding sites in each of the G-rich switch regions, and LR1 bound at contiguous sites may enhance recombination by stimulating transcription of the S regions.

Protein-DNA recognition complexes: conservation of structure and binding energy in the transition state

This paper considers how enzymes that catalyze reactions at specific DNA sites have been engineered to overcome the problem of competitive inhibition by excess nonspecific binding sites on DNA. The formation of a specific protein-DNA recognition complex is discussed from both structural and thermodynamic perspectives, and contrasted with formation of nonspecific complexes. Evidence (from EcoRI and BamHI endonucleases) is presented that a wide variety of perturbations of the DNA substrate alter binding free energy but do not affect the free energy of activation for the chemical step; that is, many energetic factors contribute equally to the recognition complex and the transition-state complex. This implies that the specific recognition complex bears a close resemblance to the transition-state complex, such that very tight binding to the recognition site on the DNA substrate does not inhibit catalysis, but instead provides energy that is efficiently utilized along the path to the transition state. It is suggested that this view can be usefully extended to "noncatalytic" site-specific DNA-binding proteins like transcriptional activators and general transcription factors.

Sequence-dependent dynamics in duplex DNA

The submicrosecond bending dynamics of duplex DNA were measured at a single site, using a site-specific electron paramagnetic resonance active spin probe. The observed dynamics are interpreted in terms of the mean squared amplitude of bending relative to the end-to-end vector defined by the weakly bending rod model. The bending dynamics monitored at the single site varied when the length and position of a repeated AT sequence, distant from the spin probe, were changed. As the distance between the probe and the AT sequence was increased, the mean squared amplitude of bending seen by the probe due to that sequence decreased. A model for the sequence-dependent internal flexural motion of duplex DNA, which casts the mean squared bending amplitudes in terms of sequence-dependent bending parameters, has been developed. The best fit of the data to the model occurs when the (AT)(n) basepairs are assumed to be 20% more flexible than the average of the basepairs within the control sequence. These findings provide a quantitative basis for interpreting the kinetics of biological processes that depend on duplex DNA flexibility, such as protein recognition and chromatin packaging.

Cytoplasmic ribosomal protein genes of the fission yeast Schizosaccharomyces pombe display a unique promoter type: a suggestion for nomenclature of cytoplasmic ribosomal proteins in databases

We identified 34 new ribosomal protein genes in the Schizosaccharomyces pombe database at the Sanger Centre coding for 30 different ribosomal proteins. All contain the Homol D-box in their promoter. We have shown that Homol D is, in this promoter type, the TATA-analogue. Many promoters contain the Homol E-box, which serves as a proximal activation sequence. Furthermore, comparative sequence analysis revealed a ribosomal protein gene encoding a protein which is the equivalent of the mammalian ribosomal protein L28. The budding yeast Saccharomyces cerevisiae has no L28 equivalent. Over the past 10 years we have isolated and characterized nine ribosomal protein (rp) genes from the fission yeast S.pombe . This endeavor yielded promoters which we have used to investigate the regulation of rp genes. Since eukaryotic ribosomal proteins are remarkably conserved and several rp genes of the budding yeast S.cerevisiae were sequenced in 1985, we probed DNA fragments encoding S.cerevisiae ribosomal proteins with genomic libraries of S.pombe . The deduced amino acid sequence of the different isolated rp genes of fission yeast share between 65 and 85% identical amino acids with their counterparts of budding yeast.

Submillimolar levels of calcium regulates DNA structure at the dinucleotide repeat (TG/AC)n

Submillimolar levels of calcium, similar to the physiological total (bound + free) intranuclear concentration (0.01-1 mM), induced a conformational change within d(TG/AC)n, one of the frequent dinucleotide repeats of the mammalian genome. This change is calcium-specific, because no other tested cation induced it and it was detected as a concentration-dependent transition from B- to a non-B-DNA conformation expanding from 3' end toward the 5' of the repeat. Genomic footprinting of various rat brain regions revealed the existence of similar non-B-DNA conformation within a d(TG/AC)28 repeat of the endogenous enkephalin gene only in enkephalin-expressing caudate nucleus and not in the nonexpressing thalamus. Binding assays demonstrated that DNA could bind calcium and can compete with calmodulin for calcium.

Interaction of a bZip oligopeptide model with oligodeoxyribonucleotides modelling DNA binding sites. The effect of flanking sequences

A leucine zipper (bZip) binding peptide BP1 was constructed based on the DNA binding sequence of the GCN4 protein, slightly modified to make it more similar to the sequence of other bZip proteins (Jun) with related DNA binding specificity. Self-complementary DNA hexadecanucleotides containing ATF/CRE, AP-1 and C/EPB target sites were used to study peptide-DNA complex formation. Conformation changes in both components that occur on complex formation were studied by circular dichroism (CD) spectroscopy. The results show that the amount of alpha-helix formed in the peptide strongly depends not only on the target site present, but also on the type of the sequence flanking the ATF/CRE target site. Highest amount of the alpha-helix induced in the peptide was observed when homopurine homopyrimidine flanking sequences were present, whereas the presence of alternating sequences, especially of the CA/TG type, showed considerably lower effects. The change in DNA conformation on complex formation was generally small, but also depended on the type of the flanking sequence. It appears that the sequences flanking the target site can considerably modify the ability of the target sequence to bind specifically the bZip peptide, probably by slightly varying the overall DNA conformation.

Proposed mechanism for the stabilization of nuclear receptor DNA binding via protein dimerization

Hepatocyte nuclear factor 4 (HNF-4) defines a new subgroup of nuclear receptors that exist in solution and bind DNA exclusively as homodimers. We recently showed that the putative ligand binding domain (LBD) of HNF-4 is responsible for dimerization in solution and prevents heterodimerization with other receptors. In this report, the role of the LBD in DNA binding by HNF-4 is further investigated by using electrophoretic mobility shift analysis. A comparison of constructs containing either the DNA binding domain (DBD) alone or the DBD plus the LBD of HNF-4 showed that dimerization via the DBD was sufficient to provide nearly the full DNA binding affinity of the full-length HNF-4. In contrast, dimerization via the DBD was not sufficient to produce a stable protein-DNA complex, whereas dimerization via the LBD increased the half-life of the complex by at least 100-fold. Circular permutation analysis showed that full-length HNF-4 bent DNA by approximately 80 degrees while the DBD bent DNA by only 24 degrees. Nonetheless, analysis of other constructs indicated that the increase in stability afforded by the LBD could be explained only partially by an increased ability to bend DNA. Coimmunoprecipitation studies, on the other hand, showed that dimerization via the LBD produced a protein-protein complex that was much more stable than the corresponding protein-DNA complex. These results led us to propose a model in which dimerization via the LBD stabilizes the receptor on DNA by converting an energetically favorable two-step dissociation event into an energetically unfavorable single-step event. Implications of this one-step model for other nuclear receptors are discussed.

Strained DNA is kinked by low concentrations of Zn2+

A novel atomic force microscope with a magnetically oscillated tip has provided unprecedented resolution of small DNA fragments spontaneously adsorbed to mica and imaged in situ in the presence of divalent ions. Kinks (localized bends of average angle 78 degrees) were observed in axially strained minicircles consisting of tandemly repeated d(A)5 and d(GGGCC[C]) sequences. The frequency of kinks in identical minicircles increased 4-fold in the presence of 1 mM Zn2+ compared with 1 mM Mg2+. Kinking persisted in mixed Mg2+/Zn2+ electrolytes until the Zn2+ concentration dropped below 100 microM, indicating that this type of kinking may occur under physiological conditions. Kinking appears to replace intrinsic bending, and statistical analysis shows that kinks are not localized within any single sequence element. A surprisingly small free energy is associated with kink formation.

Does TATA matter? A structural exploration of the selectivity determinants in its complexes with TATA box-binding protein

The binding of the TATA box-binding protein (TBP) to a TATA sequence in DNA is essential for eukaryotic basal transcription. TBP binds in the minor groove of DNA, causing a large distortion of the DNA helix. Given the apparent stereochemical equivalence of AT and TA basepairs in the minor groove, DNA deformability must play a significant role in binding site selection, because not all AT-rich sequences are bound effectively by TBP. To gain insight into the precise role that the properties of the TATA sequence have in determining the specificity of the DNA substrates of TBP, the solution structure and dynamics of seven DNA dodecamers have been studied by using molecular dynamics simulations. The analysis of the structural properties of basepair steps in these TATA sequences suggests a reason for the preference for alternating pyrimidine-purine (YR) sequences, but indicates that these properties cannot be the sole determinant of the sequence specificity of TBP. Rather, recognition depends on the interplay between the inherent deformability of the DNA and steric complementarity at the molecular interface.

Flexible DNA: genetically unstable CTG.CAG and CGG.CCG from human hereditary neuromuscular disease genes

The properties of duplex CTG.CAG and CGG.CCG, which are involved in the etiology of several hereditary neurodegenerative diseases, were investigated by a variety of methods, including circularization kinetics, apparent helical repeat determination, and polyacrylamide gel electrophoresis. The bending moduli were 1.13 x 10(-19) erg.cm for CTG and 1.27 x 10(-19) erg.cm for CGG, approximately 40% less than for random B-DNA. Also, the persistence lengths of the triplet repeat sequences were approximately 60% the value for random B-DNA. However, the torsional moduli and the helical repeats were 2.3 x 10(-19) erg.cm and 10.4 base pairs (bp)/turn for CTG and 2.4 x 10(-19) erg.cm and 10.3 bp/turn for CGG, respectively, all within the range for random B-DNA. Determination of the apparent helical repeat by the band shift assay indicated that the writhe of the repeats was different from that of random B-DNA. In addition, molecules of 224-245 bp in length (64-71 triplet repeats) were able to form topological isomers upon cyclization. The low bending moduli are consistent with predictions from crystallographic variations in slide, roll, and tilt. No unpaired bases or non-B-DNA structures could be detected by chemical and enzymatic probe analyses, two-dimensional agarose gel electrophoresis, and immunological studies. Hence, CTG and CGG are more flexible and highly writhed than random B-DNA and thus would be expected to act as sinks for the accumulation of superhelical density.

DNA bending is essential for the site-specific recognition of DNA response elements by the DNA binding domain of the tumor suppressor protein p53

We have used circular permutation assays to determine the extent and location of the DNA bend induced by the DNA binding domain of human wild type p53 (p53DBD) upon binding to several naturally occurring DNA response elements. We have found that p53DBD binding induces axial bending in all of the response elements investigated. In particular, response elements having a d(CATG) sequence at the junction of two consensus pentamers in each half-site favor highly bent complexes (bending angle is approximately 50 degrees ), whereas response elements having d(CTTG) bases at this position are less bent (bending angles from approximately 37 to approximately 25 degrees ). Quantitative electrophoretic mobility shift assays of different complexes show a direct correlation between the DNA bending angle and the binding affinity of the p53DBD with the response elements, i.e. the greater the stability of the complex, the more the DNA is bent by p53DBD binding. The study provides evidence that the energetics of DNA bending, as determined by the presence or absence of flexible sites in the response elements, may contribute significantly to the overall binding affinity of the p53DBD for different sequences. The results therefore suggest that both the structure and the stability of the p53-DNA complex may vary with different response elements. This variability may be correlated with variability in p53 function.

Architectural accommodation in the complex of four p53 DNA binding domain peptides with the p21/waf1/cip1 DNA response element

High resolution chemical footprinting and cross-linking experiments have provided a basis for elucidating the overall architecture of the complex between the core DNA binding domain of p53 (p53DBD, amino acids 98-309) and the p21/waf1/cip1 DNA response element implicated in the G1/S phase cell cycle checkpoint. These studies complement both a crystal structure and earlier biophysical studies and provide the first direct experimental evidence that four subunits of p53DBD bind to the response element in a regular staggered array having pseudodyad symmetry. The invariant guanosines in the highly conserved C(A/T)|(T/A)G parts of the consensus half-sites are critical to the p53DBD-DNA binding. Molecular modeling of the complex using the observed peptide-DNA contacts shows that when four subunits of p53DBD bind the response element, the DNA has to bend approximately 50 degrees to relieve steric clashes among different subunits, consistent with recent DNA cyclization studies. The overall lateral arrangement of the four p53 subunits with respect to the DNA loop comprises a novel nucleoprotein assembly that has not been reported previously in other complexes. We suggest that this kind of nucleoprotein superstructure may be important for p53 binding to response elements packed in chromatin and for subsequent transactivation of p53-mediated genes.

Polymorphic curvature of satellite DNA in three subspecies of the beetle Pimelia sparsa

The curvature of the monomeric repeats of satellite DNAs from three subspecies of the beetle Pimelia sparsa (Coleoptera, Tenebrionidae) has been analysed. Evidence of curvature was inferred from their retarded migration in native polyacrylamide gels, which was confirmed by direct electron microscopy visualisation. Sequence-comparison analysis, which included sequence alignments and modelling studies, was used to reveal the patterns of local bending and curvature. The effects of the minor-groove-binding drugs distamycin and berenil on the curvature were assayed, and analysed with respect to their (A+T)-rich sequence preferences. Since our study deals with satellite DNAs from closely related organisms (three subspecies), we correlated the differences in sequence, and also the high similarity conserved in the (A+T)-rich regions, with the changes in the patterns of curvature.

DNA curvature controls termination of plus strand DNA synthesis at the centre of HIV-1 genome

In vivo and in vitro, reverse transcriptase (RT) from human immunodeficiency virus type 1 (HIV-1) terminates plus strand synthesis at the centre of the viral genome. The central termination sequence (CTS) contains curved DNA fragments located upstream of each terminator site. Two different models, relying either on the A-tract or general sequence roll assumptions, were used to predict the extent and the direction of this curvature as well as to design mutants, which abolished it. Straightening of each curved element abolished termination at the site located immediately downstream from the curvature. When synthesis was performed on the other strand and in the opposite direction, the two curved elements C1 and C2 associated with the two termination sites Ter1 and Ter2, led again to termination of DNA synthesis. Therefore, termination occurred as a nascent bent duplex was synthesized within the template primer binding cleft of RT, even when putative strand-specific motifs have been removed by the inversion. Computation of DNA paths upstream of other known arrest sites suggested that this feature was of general relevance for termination. At the CTS, termination occurred more precisely at the 3' end of an AnTm motif (n + m = 7). The possible structures, adopted by this motif, are discussed and confronted with the present crystallographic and biochemical data obtained on HIV-1 RT-DNA interactions and on HIV-1 RT processivity.

A hot spot for hydrogen peroxide-induced damage in the human hypoxia-inducible factor 1 binding site of the PGK 1 gene

Using ligation-mediated polymerase chain reaction to separately map the distribution of induced oxidized bases and strand breaks along the human PGK1 promoter at nucleotide resolution, we previously described the pattern of oxidative DNA damage induced in vitro by Cu(II)/ascorbate/H2O2 [J. Biol. Chem. 270, 17633-17640 (1995)]. Here we report that the pattern of in vivo base damage caused by H2O2 is almost identical to that of the previously used in vitro system with the exception of transcription factor-associated footprints. An unusually strong positive footprint for both strand breaks and oxidized bases is associated with binding of the hypoxia-inducible transcription factor-1. Base damage at this footprint was 52-91% repaired in 24 h, which was similar to the global base damage repair rate. However, strand breaks at this footprint were only 39-55% repaired in 24 h or approximately 100-fold slower than the global strand break repair rate.

Formation and processing of UV photoproducts: effects of DNA sequence and chromatin environment

Cyclobutane pyrimidine dimers and (6-4) photoproducts are the two major classes of lesions produced in DNA by UVB and UVC irradiation. Their distribution along genes is nucleotide sequence-dependent. In vivo, the frequency of these lesions at specific sites is modulated by nucleosomes and other DNA binding proteins. Repair of UV photoproducts is dependent on the transcriptional status of the sequences to be repaired and on the chromatin environment. The formation of DNA photolesions by UV light is responsible for the induction of mutations and the development of skin cancer. To understand the mechanisms of UV mutagenesis, it is important to know how these lesions are formed, by which cellular pathways they are repaired and how they are dealt with by DNA polymerases.

Evolution of promoter sequences: elements of a canonical promoter for prespore genes of Dictyostelium

An attempt is made to define a minimal prespore promoter which contains all elements essential for correct regulation of expression of a prespore gene. The prespore genes of Dictyostelium are coregulated during development. Most begin transcription at the same early stage, and activity of all is restricted to prespore tissue during the later slug stage. Sequences 5' to the coding sequences of eight prespore genes were searched for all elements proposed to control transcription and for new elements. The meaningfulness of occurrences of elements and pairs of elements in prespore promoters was evaluated by comparison with frequencies of occurrences in promoters of other, nonprespore genes. These comparisons resulted in definition of a canonical prespore promoter, a stretch of about 200 nucleotides containing at least one of each of three elements. Certain limitations were found on the spacing of elements. Orientation of elements with respect to each other appeared unrestricted. All elements often occurred in multiple copies. This structure suggests that individual copies of each element are not conserved during evolution, but instead continually appear and disappear.

The effects of sequence context on DNA curvature

Recent experiments have exposed significant discrepancies between experimental data and predictive models for DNA structure. These results strongly suggest that DNA structural parameters incorporated in the models are not always sufficient to account for the influence of sequence context and of specific ion effects. In an attempt to evaluate these two effects, we have investigated repetitive DNA sequences with the sequence motif GAGAG.CTCTC located in different helical phasing arrangements with respect to poly(A) tracts and GGGCCC.GGGCCC sequence motifs. Methods used are ligase-mediated cyclization and gel mobility experiments along with DNase I cutting and chemical probe studies. The results provide new evidence for curvature in poly(A) tracts. They also show that the sequence context in which bending and flexible sequence elements are found is an important aspect of sequence-dependent DNA conformation. Although dinucleotide models generally have good predictive power, this work demonstrates that in some instances sequence elements larger than the dinucleotide must be taken into account, and hence it provides a starting point for the appropriate modification and refinement of existing structural models for DNA.

Large-amplitude picosecond anisotropy decay of the intrinsic fluorescence of double-stranded DNA

The conformational flexibility of the DNA double helix is of great interest because of its potential role in protein recognition, packaging into chromosomes, formation of photodefects, and interaction with drugs. Theory finds that DNA is very flexible; however, there is a scarcity of experimental results that examine intrinsic properties of the DNA bases for the inherent flexibility in solution. We have studied the dynamics of poly(dA).poly(dT) and (dA)20.(dT)20 in a 50 mM cacodylate, 0.1 M NaCl, pH 7 buffer by using the time-correlated picosecond fluorescence anisotropy of thymine selectively excited at 293 nm. For both nucleic acids, a large-amplitude biphasic decrease in the anisotropy is observed that has a very fast, large-amplitude component on the picosecond time scale and a slower, smaller-amplitude component on the nanosecond time scale. These modes are sensitive to sucrose concentration, and are greatly attenuated at 77% sucrose by volume. This observation suggests that motions of the bases make a significant contribution to the observed fluorescence depolarization (in the absence of sucrose). Measurements on the single-stranded systems poly(dT) and (dT)20 reveal a much smaller amplitude of the very fast depolarization mode. These observations are consistent with a mechanism that involves concerted motions in the interior of the double-stranded systems.

UV damage and repair mechanisms in mammalian cells

The formation of DNA photoproducts by ultraviolet (UV) light is responsible for induction of mutations and development of skin cancer. To understand UV mutagenesis, it is important to know the mechanisms of formation and repair of these lesions. Cyclobutane pyrimidine dimers and (6-4)photoproducts are the two major classes of UV-induced DNA lesions. Their distribution along DNA sequences in vivo is strongly influenced by nucleosomes and other DNA binding proteins. Repair of UV photoproducts is dependent on the transcriptional status of the sequences to be repaired and on the chromatin environment. Sensitive techniques are now available to study repair of UV damage at the level of nucleotide resolution in mammalian cells. With the aid of in vitro systems, the entire nucleotide excision repair process has been reconstituted from purified protein components with naked DNA as a substrate. Future work will focus on the development of in vitro assays for transcription-coupled repair and repair in chromatin.

Bending and torsional flexibility of G/C-rich sequences as determined by cyclization assays

The structural polymorphism of DNA is a vital aspect of its biological function. However, it has become increasingly apparent in recent years that DNA polymorphism is a complicated, multidimensional phenomenon that includes not only static sequence-directed structures but dynamic effects as well, including influences of counterions and sequence context. In order to address some of these additional factors that govern DNA conformation, we have used T4 ligase-mediated cyclization to investigate bending in a series of DNA sequences containing the GGGCCC.GGGCCC motif in different sequence contexts including various helical phasings with (A)5-tracts. We present evidence for curvature in GGGCCC.GGGCCC and (A)5-tract motifs in the presence of physiological levels of Mg2+ and show that these motifs curve through similar but oppositely directed bending angles under these ionic strength conditions. Although these two sequence motifs appear to bend similarly, our results suggest significant differences in stiffness and stability of curvature between them. We also show that under the same experimental conditions, the CTAG-CTAG sequence element possesses unusual torsional flexibility and that this appears to be associated with the central TA.TA dinucleotide. The results underscore the need to include sequence context and specific ion effects as well as a dynamic basis in more complete predictive models for functionally related DNA polymorphism.

The tandem repeat AGGGTAGGGT is, in the fission yeast, a proximal activation sequence and activates basal transcription mediated by the sequence TGTGACTG

Ribosomal protein (rp) genes in the fission yeast Schizosaccharomyces pombe display two highly conserved sequence elements in the promoter region. The molecular dissection of these promoters revealed that basal transcription is not based on a TATA element. The sequence which promotes basal transcription is the conserved sequence CAGTCACA or the inverted form TGTGACTG, called the homol D box. Upstream of the homol D box a tandem repeat AGGGTAGGGT or the inverted form ACCCTACCCT appears in some promoters, called homol E. This element functions in the proximal arrangement with homol D as an activation sequence. A compilation of homol D and homol E sequences identified in other S.pombe promoters revealed that several putative polymerase II and polymerase III promoters display a homol D box or the homol E/homol D arrangement.

Members of a new family of DNA-binding proteins bind to a conserved cis-element in the promoters of alpha-Amy2 genes

The promoters of wheat, barley and wild oat alpha-Amy2 genes contain a number of conserved cis-acting elements that bind nuclear protein, we report here the isolation of two cDNAs encoding proteins (ABF1 and ABF2) that bind specifically to one of these elements, Box 2 (ATTGACTTGACCGTCATCGG). The two proteins are unrelated to each other except for a conserved region of 56-58 amino acids that consists of 25 highly conserved amino acids followed by a putative zinc finger motif, C-X4-5-C-X22-23-H-X1-H. ABF1 contains two such conserved regions, whereas ABF2 possesses only one but also contains a potential leucine zipper motif, suggesting that it could form homo- or heterodimers. ABF1 and ABF2 expressed in Escherichia coli bound specifically to Box 2 probes in gel retardation experiments; this binding was abolished by the transition-metal-chelating agent, 1,10-o-phenanthroline and by EDTA. We propose that ABF1 and ABF2 are representatives of two classes of a new family of plant sequence-specific DNA-binding proteins.

DNA loops induced by cooperative binding of transcriptional activator proteins and preinitiation complexes

DNA conformational changes are essential for the assembly of multiprotein complexes that contact several DNA sequence elements. An approach based on atomic force microscopy was chosen to visualize specific protein-DNA interactions occurring on eukaryotic class II nuclear gene promoters. Here we report that binding of the transcription regulatory protein Jun to linearized plasmid DNA containing the consensus AP-1 binding site upstream of a class II gene promoter leads to bending of the DNA template. This binding of Jun was found to be essential for the formation of preinitiation complexes (PICs). The cooperative binding of Jun and PIC led to looping of DNA at the protein binding sites. These loops were not seen in the absence of either PICs, Jun, or the AP-1 binding site, suggesting a direct interaction between DNA-bound Jun homodimers and proteins bound to the core promoter. This direct visualization of functional transcriptional complexes confirms the theoretical predictions for the mode of gene regulation by trans-activating proteins.