Prediction of atomic structure from sequence for double helical DNA oligomers

DNA can adopt different conformations depending on the base sequence, solvent, electrolyte composition and concentration, pH, temperature, and interaction with proteins. Here we present a model for calculating the three-dimensional atomic structure of double-stranded DNA oligomers. A theoretical energy function is used for calculating the interactions within the base steps and an empirical backbone function is used to restrict the conformational space accessible to the bases and to account for the conformational coupling of neighboring steps in a sequence. Conformational searching on large structures or a large number of structures is possible, because each base step can be described by just two primary degrees of freedom (slide and shift). A genetic algorithm is used to search for low-energy structures in slide-shift space, and this allows very rapid optimization of DNA oligomers. The other base step parameters have been previously optimized for all possible slide-shift sequence combinations, and a heuristic algorithm is used to add the atomic details of the backbone conformation in the final step of the calculation. The structures obtained by this method are very similar to the corresponding X-ray crystal structures observed experimentally. The average RMSD is 2.24 Angstroms for a set of 20 oligomer structures. For 15 of these sequences, the X-ray crystal structure is the global energy minimum. The other 5 are bistable sequences that have B-form global energy minima but crystallize as A-DNA.     

A fractal method to distinguish coding and non-coding sequences in a complete genome based on a number sequence representation

A fractal method to distinguish coding and non-coding sequences in a complete genome is proposed, based on different statistical behaviors between these two kinds of sequences. We first propose a number sequence representation of DNA sequences. Multifractal analysis is then performed on the measure representation of the obtained number sequence. The three exponents C(-1), C1 and C2 are selected from the result of multifractal analysis. Each DNA may be represented by a point in the three-dimensional space generated by these three-component vectors. It is shown that points corresponding to coding and non-coding sequences in the complete genome of many prokaryotes are roughly distributed in different regions. Fisher's discriminant algorithm can be used to separate these two regions in the spanned space. If the point (C(-1),C1,C2) for a DNA sequence is situated in the region corresponding to coding sequences, the sequence is discriminated as a coding sequence; otherwise, the sequence is classified as a non-coding one. For all 51 prokaryotes we considered , the average discriminant accuracies pc,pnc,qc and qnc reach 72.28%, 84.65%, 72.53% and 84.18%, respectively.     

A structural similarity analysis of double-helical DNA

A database of the structural properties of all 32,896 unique DNA octamer sequences has been calculated, including information on stability, the minimum energy conformation and flexibility. The contents of the database have been analysed using a variety of Euclidean distance similarity measures. A global comparison of sequence similarity with structural similarity shows that the structural properties of DNA are much less diverse than the sequences, and that DNA sequence space is larger and more diverse than DNA structure space. Thus, there are many very different sequences that have very similar structural properties, and this may be useful for identifying DNA motifs that have similar functional properties that are not apparent from the sequences. On the other hand, there are also small numbers of almost identical sequences that have very different structural properties, and these could give rise to false-positives in methods used to identify function based on sequence alignment. A simple validation test demonstrates that structural similarity can differentiate between promoter and non-promoter DNA. Combining structural and sequence similarity improves promoter recall beyond that possible using either similarity measure alone, demonstrating that there is indeed information available in the structure of double-helical DNA that is not readily apparent from the sequence.     

A prediction of the amino acids and structures involved in DNA recognition by type I DNA restriction and modification enzymes.

The S subunits of type I DNA restriction/modification enzymes are responsible for recognising the DNA target sequence for the enzyme. They contain two domains of approximately 150 amino acids, each of which is responsible for recognising one half of the bipartite asymmetric target. In the absence of any known tertiary structure for type I enzymes or recognisable DNA recognition motifs in the highly variable amino acid sequences of the S subunits, it has previously not been possible to predict which amino acids are responsible for sequence recognition. Using a combination of sequence alignment and secondary structure prediction methods to analyse the sequences of S subunits, we predict that all of the 51 known target recognition domains (TRDs) have the same tertiary structure. Furthermore, this structure is similar to the structure of the TRD of the C5-cytosine methyltransferase, Hha I, which recognises its DNA target via interactions with two short polypeptide loops and a beta strand. Our results predict the location of these sequence recognition structures within the TRDs of all type I S subunits.     

DNA序列信息的一种新的测度

根据信息理论给出了测度ＤＮＡ序列信息的一种新的方法,获得ＤＮＡ序列４个层次的信息量测度：Ｉｂ,Ｉｆ（１）,Ｉｆ（２）ａｎｄＩｆ（３）,这４种信息测度可分别用来测度ＤＮＡ的碱基序列、密码子序列、编码蛋白质序列和功能蛋白质序列的信息量。从Ｍ．ｅｄｕｌｉｓ的线粒体基因组中两个较短的编码蛋白质的ＤＮＡ序列和使用具有不同倍性的间并密码子组组成的模拟ＤＮＡ序列中所获得计算结果表明,这些信息测度确实能用来揭示所     

Rec A-independent homologous recombination induced by a putative fold-back tetraplex DNA

We have recently reported that a GC-rich palindromic repeat sequence presumably adopts a stable fold-back tetraplex DNA structure under supercoiling. To establish the biological significance of this structure, we inserted this sequence between two direct repeat sequences, separated by 200 bp, in a plasmid. We then investigated the effect of this sequence on homologous recombination events. Here we report that the putative fold-back DNA tetraplex structure induces homologous recombination between direct repeat sequences. Interestingly, this recombination event is independent of recA, a major driving force for homologous recombination. We think that the fold-back structure forces the repeat sequences to come into close proximity and therefore leads to strand exchange. Although triplex-induced recombination has been well documented, our results for the first time directly establish the potential of a tetraplex structure to induce recA-independent homologous recombination in vivo. This finding might have a significant implication for site-directed gene deletion in the context of the correction of genetic defects.     

Computer modeling studies of the structure of a repressor

Due to advances in molecular biology the DNA sequences of structural genes coding for proteins are often known before a protein is characterized or even isolated. The function of a protein whose amino acid sequence has been deduced from a DNA sequence may not even be known. This has created greater interest in the development of methods to predict the tertiary structures of proteins. The a priori prediction of a protein's structure from its amino acid sequence is not yet possible. However, since proteins with similar amino acid sequences are observed to have similar three-dimensional structures, it is possible to use an analogy with a protein of known structure to draw some conclusions about the structure and properties of an uncharacterized protein. The process of predicting the tertiary structure of a protein relies very much upon computer modeling and analysis of the structure. The prediction of the structure of the bacteriophage 434 cro repressor is used as an example illustrating current procedures.     

Plant defensins: Common fold,multiple functions

Plant defensins represent a large class of structurally similar peptides found throughout the plant kingdom. Despite a conserved cysteine spacing pattern and three-dimensional structure, their sequences are highly divergent and they display a range of activities including antifungal and antibacterial activities, enzyme inhibitory activities as well as roles in heavy metal tolerance and development. The vast number of sequences along with their diverse range of activities makes it impossible to test the activity and assign function to all plant defensins. However, as the number of characterized defensins increases, in depth sequence analysis may allow us to predict the function of newly identified peptides. In this review, we analyze the sequences of defensins whose activities have been described and group these based on similarity using a maximum-likelihood phylogenetic tree. We also compare the amino acids that have been described as essential for the activity of various plant defensins between these groups. While many more plant defensins will need to be characterized before we can develop rules to predict the activity of novel sequences, this approach may prove useful in identifying structure–function relationships.     

Nucleosomal locations of dominant DNA sequence motifs for histone-DNA interactions and nucleosome positioning

DNA sequence is an important determinant of the positioning, stability, and activity of nucleosomes, yet the molecular basis of these effects remains elusive. A "consensus DNA sequence" for nucleosome positioning has not been reported and, while certain DNA sequence preferences or motifs for nucleosome positioning have been discovered, how they function is not known. Here, we report that an unexpected observation concerning the reassembly of nucleosomes during salt gradient dialysis has allowed a breakthrough in our efforts to identify the nucleosomal locations of the DNA sequence motifs that dominate histone-DNA interactions and nucleosome positioning. We conclude that a previous selection experiment for high-affinity, nucleosome-forming DNA sequences exerted selective pressure chiefly on the central stretch of the nucleosomal DNA. This observation implies that algorithms for aligning the selected DNA sequences should seek to optimize the alignment over much less than the full 147 bp of nucleosomal DNA. A new alignment calculation implemented these ideas and successfully aligned 19 of the 41 sequences in a non-redundant database of selected high-affinity, nucleosome-positioning sequences. The resulting alignment reveals strong conservation of several stretches within a central 71 bp of the nucleosomal DNA. The alignment further reveals an inherent palindromic symmetry in the selected DNAs; it makes testable predictions of nucleosome positioning on the aligned sequences and for the creation of new positioning sequences, both of which are upheld experimentally; and it suggests new signals that may be important in translational nucleosome positioning.     

Plant defensins: Common fold,multiple functions

Plant defensins represent a large class of structurally similar peptides found throughout the plant kingdom. Despite a conserved cysteine spacing pattern and three-dimensional structure, their sequences are highly divergent and they display a range of activities including antifungal and antibacterial activities, enzyme inhibitory activities as well as roles in heavy metal tolerance and development. The vast number of sequences along with their diverse range of activities makes it impossible to test the activity and assign function to all plant defensins. However, as the number of characterized defensins increases, in depth sequence analysis may allow us to predict the function of newly identified peptides. In this review, we analyze the sequences of defensins whose activities have been described and group these based on similarity using a maximum-likelihood phylogenetic tree. We also compare the amino acids that have been described as essential for the activity of various plant defensins between these groups. While many more plant defensins will need to be characterized before we can develop rules to predict the activity of novel sequences, this approach may prove useful in identifying structure–function relationships.     

Site-directed mutagenesis studies of tn5 transposase residues involved in synaptic complex formation

Transposition (the movement of discrete segments of DNA, resulting in rearrangement of genomic DNA) initiates when transposase forms a dimeric DNA-protein synaptic complex with transposon DNA end sequences. The synaptic complex is a prerequisite for catalytic reactions that occur during the transposition process. The transposase-DNA interactions involved in the synaptic complex have been of great interest. Here we undertook a study to verify the protein-DNA interactions that lead to synapsis in the Tn5 system. Specifically, we studied (i) Arg342, Glu344, and Asn348 and (ii) Ser438, Lys439, and Ser445, which, based on the previously published cocrystal structure of Tn5 transposase bound to a precleaved transposon end sequence, make cis and trans contacts with transposon end sequence DNA, respectively. By using genetic and biochemical assays, we showed that in all cases except one, each of these residues plays an important role in synaptic complex formation, as predicted by the cocrystal structure.     

在AppleⅡ型微机上实现核酸数据处理的工作程序

本文报道了在AppleⅡ型微机上实现核酸数据处理的一系列工作程序。应用这些程序,可进行核酸数据的贮存、对指定的核酸数据结构的改造、限制性内切酶识别位点的检索、核酸序列至蛋白序列的翻译、相关核酸序列及蛋白序列的同源性比较、氨基酸密码使用频率的统计和基因的启动子结构的初步探索等方面的工作。     

Reproducible but dynamic positioning of DNA in chromosomes during mitosis

How DNA is folded into chromosomes is unknown. Mitotic chromosome banding shows reproducibility in longitudinal compaction at a resolution of several megabase pairs, but it is less clear whether DNA sequences are targeted laterally to specific locations. The in vitro chromosome assembly of prokaryotic DNA suggests that there is a lack of sequence requirements for chromosome condensation, implying an absence of DNA targeting. Protein extraction experiments indicate, however, that specific DNA sequences may bind to a chromosome scaffold. Chromosome banding patterns, using dyes with differential sequence specificity, have been interpreted to result from the alignment of AT-rich sequences in a partially helically folded chromosome scaffold. But fluorescence in situ hybridization experiments, perhaps owing to technical limitations, have shown at best only slight deviation from a random, lateral sequence distribution. Here we show that there is highly reproducible targeting of specific chromosome segments to the metaphase chromatid axis, but that these segments localize to the periphery of prophase and telophase chromosomes. Unfolding intermediates during anaphase and telophase suggest that sequence repositioning occurs through the global uncoiling of an underlying chromatid structure.     

Stability, structure and complexity of yeast chromosome III.

The complete sequence of yeast chromosome III provides a model for studies relating DNA sequence and structure at different levels of organisation in eukaryotic chromosomes. DNA helical stability, intrinsic curvature and sequence complexity have been calculated for the complete chromosome. These features are compartmentalised at different levels of organisation. Compartmentalisation of thermal stability is observed from the level delineating coding/non-coding sequences, to higher levels of organisation which correspond to regions varying in G + C content. The three-dimensional path reveals a symmetrical structure for the chromosome, with a densely packed central region and more diffuse and linear subtelomeric regions. This interspersion of regions of high and low curvature is reflected at lower levels of organisation. Complexity of n-tuplets (n = 1 to 6) also reveals compartmentalisation of the chromosome at different levels of organisation, in many cases corresponding to the structural features. DNA stability, conformation and complexity delineate telomeres, centromere, autonomous replication sequences (ARS), transposition hotspots, recombination hotspots and the mating-type loci.     

The role of internal packing interactions in determining the structure and stability of a protein

Cassette mutagenesis has been used to investigate how internal packing interactions help to specify a protein's three-dimensional structure and stability. Three interacting residues in the hydrophobic core of the N-terminal domain of lambda repressor were randomized combinatorially. The randomization was restricted to the five amino acids Val, Leu, Ile, Met and Phe, thereby generating a sterically diverse set of core sequences composed solely of hydrophobic residues. We have isolated 78 of the 125 possible sequences generated by this randomization. Approximately 70% of the isolated sequences show some level of biological activity, and thus still carry sufficient information to encode the basic structure of lambda repressor. An assay based on the temperature dependence of activity in vivo has been used to estimate the relative activities and thermal stabilities of the set of mutants. In addition, nine mutants have been purified and their stabilities and DNA binding activities characterized in vitro. Of the 56 active sequences, only two, in addition to the wild-type, maintain the wild-type level of stability and activity. All three of these proteins satisfy stringent requirements for specifically shaped residues at each position. All of the remaining active sequences have reduced stabilities and/or reduced DNA binding affinities. These and previous results suggest that there are two levels of structural information encoded in core residues. At the first level, the basic structural information appears to reside largely in the hydrophobic character of these residues. The majority of sequences that simply maintain hydrophobicity at core positions are able to adopt the overall lambda repressor fold and maintain moderate stability. At the second, more detailed level, specific steric features of these residues and their packing interactions clearly act as important determinants of the protein's precise structure and stability. These results imply that many of the basic structural features of a protein could be predicted from relatively simple, degenerate sequence patterns.     

Sequence comparison of single-stranded DNA binding proteins and its structural implications

The primary sequences were compared among several proteins: gene product 5 protein (GP5) from phage M13; PIKE from phage Ike; gene product 32 protein (GP32) from phage T4; RecA, SSB and SSF from Escherichia coli. These proteins bind strongly and cooperatively to single-stranded DNA with no sequence specificity. GP5 is the smallest in this group and its three-dimensional structure is well-characterized. Using the entire sequence of GP5 as a template we searched for the regions in other single-stranded DNA binding proteins yielding the best alignment of aromatic and basic residues. The identified domains show alignment of five aromatic and four charged residues in these proteins. The domains in PIKE, GP32 and RecA exhibit statistically significant sequence homology with GP5. These observations strongly favor the hypothesis that the protein-single-stranded DNA complex in this class of proteins is stabilized by the stacking interaction of the aromatic residues with the bases of the DNA, and by the electrostatic interaction of the basic residues with the phosphate groups of the DNA. We also find that the DNA binding domains of these proteins have similar secondary structural preferences, mainly beta structures. The triple-stranded beta-sheet may be a common motif in the DNA binding domains of these proteins.     

A comparison of two cloned mouse beta-globin genes and their surrounding and intervening sequences

The BALC/c mouse has two nonallelic beta-globin genes that appear to reside on two different Eco R1 fragments of genomic DNA. We have already cloned one of these fragments and shown that the gene encoded within it is interrupted by at least one large intervening sequence of DNA. We have now cloned and characterized the second beta-globin gene-containing fragment. The coding sequence of its gene is also interrupted by an intervening sequence of DNA that occurs in about the same position, relative to the coding sequence, as does the first. Because some shared features of the structure of these two genes might be responsible for their coordinate expression and the elimination of their intervening sequences, we have compared their surrounding, coding and intervening sequences by restriction endonuclease analysis and by visualization of the heteroduplex structures formed between them. Of the 7000 bp of sequence compared in this way, we find only a few hundred base pairs of homology in addition to the coding sequence. These shared sequences flank the coding sequence and appear to include only those portions of the intervening sequence immediately adjacent to the interrupted structural gene.