Direct imaging of DNA motif sequences with encoded nanoparticles

We present a method for encoded tagging and imaging of short nucleic acid motif chains (oligomotifs) using selective hybridization of heterogeneous Au nanoparticles (Au-NP). The resulting encoded NP string is thus representative of the underlying motif sequence. As the NPs are much more massive than the motifs, the motif chain order can be directly observed using scanning electron microscopy. Using this technique we demonstrate direct sequencing of oligomotifs in single DNA molecules consisting of four 100-nt motif chains tagged with four different types of NPs. The method outlined is a precursor for a high density direct sequencing technology.     

Autophosphorylation of phosphorylase kinase. Divalent metal cation and nucleotide dependency

This study reports on the divalent metal ion specificity for phosphorylase kinase autophosphorylation and, in particular, provides a comparison between the efficacy of Mg2+ and Mn2+ in this role. As well as requiring Ca2+ plus divalent metal ion-ATP2- as substrate, both phosphorylase kinase autoactivation and phosphorylase conversion are additionally modulated by divalent cations. However, these reactions are affected differently by different ions. Phosphorylase kinase-catalyzed phosphorylase conversion is maximally enhanced by a 4- to 10-fold lower concentration of Mg2+ than is autocatalysis and, whereas both reactions are stimulated by Mg2+, autophosphorylation is activated by Mn2+, Co2+, and Ni2+ while phosphorylase a formation is inhibited. This difference may be due to an effect of free Mn2+ on phosphorylase rather than the inability of phosphorylase kinase to use MnATP as a substrate when catalyzing phosphorylase conversion since Mn2+, when added at a level which minimally decreases [MgATP], greatly inhibits phosphorylase phosphorylation. The interactions of Mn2+ with phosphorylase kinase are different from those of Mg2+. Not only are the effects of these ions on phosphorylase activation opposite, but they also provoke different patterns of subunit phosphorylation during phosphorylase kinase autocatalysis. With Mn2+, the time lag of phosphorylation of both the alpha and beta subunits of phosphorylase kinase in autocatalysis is diminished in comparison to what is observed with Mg2+, and the beta subunit is only phosphorylated to a maximum of 1 mol/mol of subunit. With both Mg2+ and Mn2+ the alpha subunit is phosphorylated to a level in excess of 3 mol/mol, a level similar to that obtained for beta subunit phosphorylation in the presence of Mg2+. The support of autophosphorylation by both Co2+ and Ni2+ has characteristics similar to those observed with Mn2+. Although Mn2+ stimulation of autophosphorylation occurs at levels much higher than normal physiological levels, the possible potential of phosphorylase kinase autophosphorylation as a control mechanism is illustrated by the 80- to 100-fold activation that occurs in the presence of Mn2+, a level far in excess of the enzyme activity change normally seen with covalent modification. Autophosphorylation of phosphorylase kinase demonstrates a Km for Mg X ATP2- of 27.7 microM and a Ka for Mg2+ of 3.1 mM. The reaction mechanism of autophosphorylation is intramolecular. This latter observation may indicate that phosphorylase kinase autocatalysis could be of potential physiological relevance and could occur with equal facility in cells containing either constitutively high or low levels of this enzyme.     

A new method of global search for functional DNA segments using a fractal representation of nucleotide texts]

New global method for computer prediction of functional sites in nucleotide sequences, based on the fractal representation, is presented. Fractal representation of set of sequences (FRS) provides simple way for generating recognitions matrix of functionally similar sequences and simple estimations of its efficiency for searching homologous regions in new sequences. Other advantages of the method are absence of the necessity of sequences alignment during generating based set and searching new homologous regions and small CPU time. Usage of the method illustrated for searching globin and histone genes, for ALU repeats in human genome and long terminal repeats in virus genome.     

DNA sequence representation without degeneracy

Graphical representation of DNA sequence provides a simple way of viewing, sorting and comparing various gene structures. A new two-dimensional graphical representation method using a two- quadrant Cartesian coordinates system has been derived for mathematical denotation of DNA sequence. The two-dimensional graphic representation resolves sequences’ degeneracy and is mathematically proven to eliminate circuit formation. Given x-projection and y-projection of any point on the graphical representation, the number of A, G, C and T from the beginning of the sequence to that point could be found. Compared with previous methods, this graphical representation is more in-line with the conventional recognition of linear sequences by molecular biologists, and also provides a metaphor in two dimensions for local and global DNA sequence comparison.     

DNA motif elucidation using belief propagation

Protein-binding microarray (PBM) is a high-throughout platform that can measure the DNA-binding preference of a protein in a comprehensive and unbiased manner. A typical PBM experiment can measure binding signal intensities of a protein to all the possible DNA k-mers (k = 8 ∼10); such comprehensive binding affinity data usually need to be reduced and represented as motif models before they can be further analyzed and applied. Since proteins can often bind to DNA in multiple modes, one of the major challenges is to decompose the comprehensive affinity data into multimodal motif representations. Here, we describe a new algorithm that uses Hidden Markov Models (HMMs) and can derive precise and multimodal motifs using belief propagations. We describe an HMM-based approach using belief propagations (kmerHMM), which accepts and preprocesses PBM probe raw data into median-binding intensities of individual k-mers. The k-mers are ranked and aligned for training an HMM as the underlying motif representation. Multiple motifs are then extracted from the HMM using belief propagations. Comparisons of kmerHMM with other leading methods on several data sets demonstrated its effectiveness and uniqueness. Especially, it achieved the best performance on more than half of the data sets. In addition, the multiple binding modes derived by kmerHMM are biologically meaningful and will be useful in interpreting other genome-wide data such as those generated from ChIP-seq. The executables and source codes are available at the authors’ websites: e.g. http://www.cs.toronto.edu/∼wkc/kmerHMM.     

DNA motif associated with meiotic double-strand break regions in Saccharomyces cerevisiae

Meiotic recombination in yeast is initiated by DNA double-strand breaks (DSBs) that occur at preferred sites, distributed along the chromosomes. These DSB sites undergo changes in chromatin structure early in meiosis, but their common features at the level of DNA sequence have not been defined until now. Alignment of 1 kb sequences flanking six well-mapped DSBs has allowed us to define a flexible sequence motif, the CoHR profile, which predicts the great majority of meiotic DSB locations. The 50 bp profile contains a poly(A) tract in its centre and may have several gaps of unrelated sequences over a total length of up to 250 bp. The major exceptions to the correlation between CoHRs and preferred DSB sites are at telomeric regions, where DSBs do not occur. The CoHR sequence may provide the basis for understanding meiosis-induced chromatin changes that enable DSBs to occur at defined chromosomal sites.     

phi 29 DNA polymerase residue Leu384, highly conserved in motif B of eukaryotic type DNA replicases,is involved in nucleotide insertion fidelity

Replicative DNA polymerases achieve insertion fidelity by geometric selection of a complementary nucleotide followed by induced fit: movement of the fingers subdomain toward the active site to enclose the incoming and templating nucleotides generating a binding pocket for the nascent base pair. Several residues of motif B of DNA polymerases from families A and B, localized in the fingers subdomain, have been described to be involved in template/primer binding and dNTP selection. Here we complete the analysis of this motif, which has the consensus "KLX2NSXYG" in DNA polymerases from family B, characterized by mutational analysis of conserved leucine, Leu384 of phi 29 DNA polymerase. Mutation of Leu384 into Arg resulted in a phi 29 DNA polymerase with reduced nucleotide insertion fidelity during DNA-primed polymerization and protein-primed initiation reactions. However, the mutation did not alter the intrinsic affinity for the different dNTPs, as shown in the template-independent terminal protein-deoxynucleotidylation reaction. We conclude that Leu384 of phi 29 DNA polymerase plays an important role in positioning the templating nucleotide at the polymerization active site and in controlling nucleotide insertion fidelity. This agrees with the localization of the corresponding residue in the closed ternary complexes of family A and family B DNA polymerases, contributing to form the binding pocket for the nascent base pair. As an additional effect, mutant polymerase L384R was strongly reduced in DNA binding, resulting in reduced processivity during polymerization.     

Using structural motif templates to identify proteins with DNA binding function

This work describes a method for predicting DNA binding function from structure using 3-dimensional templates. Proteins that bind DNA using small contiguous helix–turn–helix (HTH) motifs comprise a significant number of all DNA-binding proteins. A structural template library of seven HTH motifs has been created from non-homologous DNA-binding proteins in the Protein Data Bank. The templates were used to scan complete protein structures using an algorithm that calculated the root mean squared deviation (rmsd) for the optimal superposition of each template on each structure, based on C_α backbone coordinates. Distributions of rmsd values for known HTH-containing proteins (true hits) and non-HTH proteins (false hits) were calculated. A threshold value of 1.6 Å rmsd was selected that gave a true hit rate of 88.4% and a false positive rate of 0.7%. The false positive rate was further reduced to 0.5% by introducing an accessible surface area threshold value of 990 Ų per HTH motif. The template library and the validated thresholds were used to make predictions for target proteins from a structural genomics project.     

Compensation for nucleotide bias in a genome by representation as a discrete channel with noise

MOTIVATION: Calculation of the information content of motifs in genomes highly biased in nucleotide composition is likely to lead to overestimates of the amount of useful information in the motif. Calculating relative information can compensate for biases, however the resulting information content is the amount seen by an observer and not by a macromolecule binding to the motif. The latter is needed to calculate the discriminatory power of the motif and to compare motifs between species. RESULTS: By treating a biased genome as a discrete channel with noise, in accordance with Shannon Information Theory, we were able to remove both 'Distortion' and 'Noise' from the motif and recover a more instructive biological 'signal.' A Java application, LogoPaint, was developed to remove nucleotide bias distortion and triplet frequency noise from motifs, calculate information content and present the motif as a logo. We demonstrate how this technique can 'unmask' motifs in the translation initiation regions of bacteria that are obscured by strong sequence biases. AVAILABILITY: LogoPaint is available to all users from the authors as an executable JAR file. Source code is available by arrangement.     

DNA 序列中模式发现的一种快速算法

模式发现是生物信息学的一个重要研究方向,但目前的大部分算法还不能保证获得最优的模式．文章推导了针对三个序列片段相似性关系的判据,将其作为剪枝规则,提出并实现了一种深度优先的穷举搜索算法——判据搜索算法(criterion search algorithm,CRISA),理论分析表明,对绝大多数模式发现问题,CRISA具有多项式的计算时间复杂度和线性的空间复杂度。对仿真的和实际的生物序列数据的测试也表明,CRISA能够快速而完全地识别出序列中所有的模式,具有优于其它算法的总体评价,能够应用于实际的模式发现问题。     

Pyrimidine motif triplexes containing polypurine RNA or DNA with oligo 2'-O-methyl or DNA triplex forming oligonucleotides

Triplex forming oligonucleotides (TFOs) are potentially useful in targeting RNA for antisense therapeutic applications. To determine the feasibility of targeting polypurine RNA with nuclease-resistant oligonucleotides, TFOs containing 2'-deoxy or 2'-O-methyl (2'-OMe) backbones, designed to form pyrimidine motif triplexes with RNA, were synthesized. TFOs were made which can form trimolecular triplexes, or bimolecular, 'clamp' triplexes with polypurine RNA and DNA. It was found that the relative stabilities of the triplexes formed followed the order: M.DM(clamp)>D.DD approximately M.DD>M. RM>D.DM>M.RD approximately M.DM, where M is a 2'-OMe, D is a DNA and R is an RNA backbone. The third strand is listed first, separated by a dot from the purine strand of the Watson-Crick duplex, followed by the pyrimidine strand of the duplex. The results described here provide insight into the feasibility of using TFOs containing a 2'-OMe backbone as antisense agents.     

DNA homoduplexes containing no pyrimidine nucleotide

We show using polyacrylamide gel electrophoresis that guanine+adenine repeat strands of DNA associate into homoduplexes at neutral pH and moderate ionic strength. The homoduplexes melt in a cooperative way like the Watson-Crick duplex, although they contain no Watson-Crick base pair. Guanine is absolutely needed for the homoduplex formation and the homoduplex stability increases with the guanine content of the repeat. The present results have implications for the nature of the first replicators, as well as regarding forces stabilizing the duplexes of DNA.     

Polymerase zeta dependency of increased adaptive mutation frequencies in nucleotide excision repair-deficient yeast strains

Reversions of an auxotrophy-causing frameshift allele during prolonged starvation of yeast cells were used as a means to elucidate the mechanisms concerned with the generation of spontaneous adaptive mutations in cell cycle-arrested cells. Whereas about 50% of these reversions were previously shown to depend on the non-homologous end joining pathway of DNA double-strand break repair, the origin of the residual 50% remains unknown. In search for a mechanism for generation of the latter fraction of reversions we examined the role of the translesion synthesis (TLS) polymerases zeta, eta and Rev1p in cells with wild-type or impaired nucleotide excision repair (NER) capacity. The basal level of adaptive mutations in the repair-proficient wild type was not influenced by disruptions of the genes coding for these three TLS polymerases. Intriguingly, a deficiency in NER by disruption of RAD14, RAD16 or RAD26 resulted in a significantly higher frequency of adaptive mutation, yet this increase was strictly dependent on an intact REV3 gene, coding for the catalytic subunit of polymerase zeta. Furthermore, we observed that intact REV3 was also required for the occurrence of increased frequencies of adaptive mutants in the NER-proficient wild type following UV irradiation. While in proliferating cells the translesion synthesis function of polymerase zeta is connected to DNA replication, our data suggest that in cell cycle-arrested cells this enzyme is able to carry out either TLS or error-prone polymerization along an undamaged template in the course of repair processes. Such a hitherto unappreciated activity of polymerase zeta in non-replicating cells may contribute to the incidence of mutations in evolution, aging and cancer.     

DNA binding sites: representation and discovery

The purpose of this article is to provide a brief history of the development and application of computer algorithms for the analysis and prediction of DNA binding sites. This problem can be conveniently divided into two subproblems. The first is, given a collection of known binding sites, develop a representation of those sites that can be used to search new sequences and reliably predict where additional binding sites occur. The second is, given a set of sequences known to contain binding sites for a common factor, but not knowing where the sites are, discover the location of the sites in each sequence and a representation for the specificity of the protein.